No Arabic abstract
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]
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.
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.
The physical state of the gas in the central 500 pc of NGC~5128 (the radio galaxy Centaurus A - Cen A), was investigated using the far-infrared fine-structure lines of carbon, oxygen, and nitrogen, as well as the CO(4-3) molecular line. The circumnuclear disk (CND) is traced by emission from dust and the neutral gas ([CI] and CO). A gas outflow with a line-of-sight velocity of 60 km/s is evident in both species. The center of the CND is bright in [OI], [OIII], and [CII]; [OI]63mu emission dominates that of [CII] even though it is absorbed with optical depths of 1.0-1.5. The outflow is well-traced by the [NII] and [NIII] lines and also seen in the [CII] and [OIII] lines that peak in the center. Ionized gas densities are moderate in the CND and low everywhere else. Neutral gas densities range from 4000 per cm3 (outflow, extended thin disk ETD) to 20 000 per cm3 (CND). The CND radiation field is weak compared to the ETD starburst field. The outflow has a much stronger radiation field. The total mass of all the CND gas is 9 x 10^(7) M(o) and the mass of the outflowing gas is only 15%-30% of that. The outflow most likely originates from the shock-dominated CND cavity surrounding the central black hole. With a factor of three uncertainty, the mass outflow rate is about 2 M(o)/yr, a thousand times higher than the accretion rate of the black hole. Without replenishment, the CND will be depleted in 15-120 million years. However, the outflow velocity is well below the escape velocity.
We present a detailed, photoionization modeling analysis of XMM-Newton/Reflection Grating Spectrometer observations of the Seyfert 2 galaxy NGC 1068. The spectrum, previously analyzed by Kinkhabwala et al. (2002), reveals a myriad of soft-Xray emission lines, including those from H- and He-like carbon, nitrogen, oxygen, and neon, and M- and L-shell iron. As noted in the earlier analysis, based on the narrowness of the radiative recombination continua, the electron temperatures in the emission-line gas are consistent with photoionization, rather than collisional ionization. The strengths of the carbon and nitrogen emission lines, relative to those of oxygen, suggest unusual elemental abundances, which we attribute to star-formation history of the host galaxy. Overall, the emission-lines are blue-shifted with respect to systemic, with radial velocities ~ 160 km/s, similar to that of [O III] 5007, and thus consistent with the kinematics and orientation of the optical emission-line gas and, hence, likely part of an AGN-driven outflow. We were able to achieve an acceptable fit to most of the strong emission-lines with a two-component photoionization model, generated with Cloudy. The two components have ionization parameters and column densities of logU = -0.05 and 1.22, and logN(H) = 20.85 and 21.2, and covering factors of 0.35 and 0.84, respectively. The total mass of the X-ray gas is roughly of an order of magnitude greater than the mass of ionized gas determined from optical and near-IR spectroscopy, which indicates that it may be the dominant component of the narrow line region. Furthermore, we suggest that the medium which produces the scattered/polarized optical emission in NGC~1068 possesses similar physical characteristics to those of the more highly-ionized of the X-ray model components.
We present an analysis of [OI]63, [OIII]88, [NII]122 and [CII]158 far-infrared (FIR) fine-structure line observations obtained with Herschel/PACS, for ~240 local luminous infrared galaxies (LIRGs) in the Great Observatories All-sky LIRG Survey (GOALS). We find pronounced declines -deficits- of line-to-FIR-continuum emission for [NII]122, [OI]63 and [CII]158 as a function of FIR color and infrared luminosity surface density, $Sigma_{rm IR}$. The median electron density of the ionized gas in LIRGs, based on the [NII]122/[NII]205 ratio, is $n_{rm e}$ = 41 cm$^{-3}$. We find that the dispersion in the [CII]158 deficit of LIRGs is attributed to a varying fractional contribution of photo-dissociation-regions (PDRs) to the observed [CII]158 emission, f([CII]PDR) = [CII]PDR/[CII], which increases from ~60% to ~95% in the warmest LIRGs. The [OI]63/[CII]158PDR ratio is tightly correlated with the PDR gas kinetic temperature in sources where [OI]63 is not optically-thick or self-absorbed. For each galaxy, we derive the average PDR hydrogen density, $n_{rm H}$, and intensity of the interstellar radiation field, in units of G$_0$, and find G$_0$/$n_{rm H}$ ratios ~0.1-50 cm$^3$, with ULIRGs populating the upper end of the distribution. There is a relation between G$_0$/$n_{rm H}$ and $Sigma_{rm IR}$, showing a critical break at $Sigma_{rm IR}^{star}$ ~ 5 x 10$^{10}$ Lsun/kpc$^2$. Below $Sigma_{rm IR}^{star}$, G$_0$/$n_{rm H}$ remains constant, ~0.32 cm$^3$, and variations in $Sigma_{rm IR}$ are driven by the number density of star-forming regions within a galaxy, with no change in their PDR properties. Above $Sigma_{rm IR}^{star}$, G$_0$/$n_{rm H}$ increases rapidly with $Sigma_{rm IR}$, signaling a departure from the typical PDR conditions found in normal star-forming galaxies towards more intense/harder radiation fields and compact geometries typical of starbursting sources.