No Arabic abstract
We explore the warm molecular and ionized gas in the centers of two megamaser disk galaxies using K-band spectroscopy. Our ultimate goal is to determine how gas is funneled onto the accretion disk, here traced by megamaser spots on sub-pc scales. We present NIR IFU data with a resolution of ~50 pc for two galaxies: NGC 4388 with VLT/SINFONI and NGC 1194 with Keck/OSIRIS+AO. The high spatial resolution and rich spectral diagnostics allow us to study both the stellar and gas kinematics as well as gas excitation on scales only an order of magnitude larger than the maser disk. We find a drop in the stellar velocity dispersion in the inner ~100 pc of NGC 4388, a common signature of a dynamically cold central component seen in many active nuclei. We also see evidence for non-circular gas motions in the molecular hydrogen on similar scales, with the gas kinematics on 100-pc scales aligned with the megamaser disk. In contrast, the high ionization lines and Br-gamma trace outflow along the 100 pc-scale jet. In NGC 1194, the continuum from the accreting black hole is very strong, making it difficult to measure robust two-dimensional kinematics, but the spatial distribution and line ratios of the molecular hydrogen and Br-gamma have consistent properties between the two galaxies.
We present a 190-307 GHz broadband spectrum obtained with Z-Spec of NGC 1068 with new measurements of molecular rotational transitions. After combining our measurements with those previously published and considering the specific geometry of this Seyfert 2 galaxy, we conduct a multi-species Bayesian likelihood analysis of the density, temperature, and relative molecular abundances of HCN, HNC, CS, and HCO+. We find that these molecules trace warm (T > 100 K) gas of H2 number densities 10^4.2 - 10^4.9 cm^-3. Our models also place strong constraints on the column densities and relative abundances of these molecules, as well as on the total mass in the circumnuclear disk. Using the uniform calibration afforded by the broad Z-Spec bandpass, we compare our line ratios to X-ray dominated region (XDR) and photon-dominated region models. The majority of our line ratios are consistent with the XDR models at the densities indicated by the likelihood analysis, lending substantial support to the emerging interpretation that the energetics in the circumnuclear disk of NGC 1068 are dominated by accretion onto an active galactic nucleus.
Using HST, we identify circumnuclear ($100$-$500$ pc scale) structures in nine new H$_2$O megamaser host galaxies to understand the flow of matter from kpc-scale galactic structures down to the supermassive black holes (SMBHs) at galactic centers. We double the sample analyzed in a similar way by Greene et al. (2013) and consider the properties of the combined sample of 18 sources. We find that disk-like structure is virtually ubiquitous when we can resolve $<200$ pc scales, in support of the notion that non-axisymmetries on these scales are a necessary condition for SMBH fueling. We perform an analysis of the orientation of our identified nuclear regions and compare it with the orientation of megamaser disks and the kpc-scale disks of the hosts. We find marginal evidence that the disk-like nuclear structures show increasing misalignment from the kpc-scale host galaxy disk as the scale of the structure decreases. In turn, we find that the orientation of both the $sim100$ pc scale nuclear structures and their host galaxy large-scale disks is consistent with random with respect to the orientation of their respective megamaser disks.
We report the detection of far-IR CO rotational emission from the prototypical Seyfert 2 galaxy NGC 1068. Using Herschel-PACS, we have detected 11 transitions in the J_upper=14-30 (E_upper/k_B = 580-2565 K) range, all of which are consistent with arising from within the central 10 (700 pc). The detected transitions are modeled as arising from 2 different components: a moderate excitation (ME) component close to the galaxy systemic velocity, and a high excitation (HE) component that is blueshifted by ~80 km s^{-1}. We employ a large velocity gradient (LVG) model and derive n_H2~10^{5.6} cm^{-3}, T_kin~170 K, and M_H2~10^{6.7} M_sun for the ME component, and n_H2~10^{6.4} cm^{-3}, T_kin~570 K, and M_H2~10^{5.6} M_sun for the HE component, although for both components the uncertainties in the density and mass are plus/minus (0.6-0.9) dex. We compare the CO line profiles with those of other molecular tracers observed at higher spatial and spectral resolution, and find that the ME transitions are consistent with these lines arising in the ~200 pc diameter ring of material traced by H_2 1-0 S(1) observations. The blueshift of the HE lines may also be consistent with the bluest regions of this H_2 ring, but a better kinematic match is found with a clump of infalling gas ~40 pc north of the AGN. We discuss the prospects of placing the HE component near the AGN, and conclude that while the moderate thermal pressure precludes an association with the ~1 pc radius H_2O maser disk, the HE component could potentially be located only a few parsecs more distant from the AGN, and might then provide the N_H~10^{25} cm^{-2} column obscuring the nuclear hard X-rays. Finally, we also report sensitive upper limits extending up to J_upper=50, which place constraints on a previous model prediction for the CO emission from the X-ray obscuring torus. [Abridged]
A joint XMM-Newton and NuSTAR observation was conducted for the bright, local Seyfert 1.9 galaxy, NGC 1194. The hard spectral form of this AGN was modeled using the toroidal reprocessor MYTORUS. The decoupled model form provides a good description of the spectrum, with reflection arising from gas with a global average column density > 4 x 10^24 cm^-2 and transmission of the continuum through an order-of-magnitude lower column. In this model, the reflection strength is a factor ~3 higher than expected from a simple torus. Such a result may indicate that much of the intrinsic X-ray continuum is hidden from view. An alternative model is that of a patchy torus, where 85% of sight-lines are obscured by Compton-thick gas and the remaining 15% by Compton-thin gas. The patchy torus model is based on a solar abundance of Fe and is consistent with X-ray partial-covering results found in other AGN. That a patchy torus model would relieve the issue with the strength of the reflection signature is not an intuitive result: such an insight regarding the geometry of the global reprocessing gas could not have been obtained using ad hoc model components to describe the spectral form.
We present the first observations of H$^{13}$CN$(1-0)$, H$^{13}$CO$^+(1-0)$ and SiO$(2-1)$ in NGC,6240, obtained with the IRAM PdBI. Combining a Markov Chain Monte Carlo (MCMC) code with Large Velocity Gradient (LVG) modelling, and with additional data from the literature, we simultaneously fit three gas phases and six molecular species to constrain the physical condition of the molecular gas, including mass$-$luminosity conversion factors. We find $sim10^{10}M_odot$ of dense molecular gas in cold, dense clouds ($T_{rm k}sim10$,K, $n_{{rm H}_2}sim10^6$,cm$^{-3}$) with a volume filling factor $<0.002$, embedded in a shock heated molecular medium ($T_{rm k}sim2000$,K, $n_{{rm H}_2}sim10^{3.6}$,cm$^{-3}$), both surrounded by an extended diffuse phase ($T_{rm k}sim200$,K, $n_{{rm H}_2}sim10^{2.5}$,cm$^{-3}$). We derive a global $alpha_{rm CO}=1.5^{7.1}_{1.1}$ with gas masses $log_{10}left(M / [M_odot]right)=10.1_{10.0}^{10.8}$, dominated by the dense gas. We also find $alpha_{rm HCN} = 32^{89}_{13}$, which traces the cold, dense gas. The [$^{12}$C]/[$^{13}$C] ratio is only slightly elevated ($98^{230}_{65}$), contrary to the very high [CO]/[$^{13}$CO] ratio (300-500) reported in the literature. However, we find very high [HCN]/[H$^{13}$CN] and [HCO$^+$]/[H$^{13}$CO$^+$] abundance ratios $(300^{500}_{200})$ which we attribute to isotope fractionation in the cold, dense clouds.