We present the results of our ALMA observations of 11 (ultra)luminous infrared galaxies ((U)LIRGs) at J = 4-3 of HCN, HCO+, and HNC and J = 3-2 of HNC. This is an extension of our previously ...published HCN and HCO+ J = 3-2 observations to multiple rotational J-transitions of multiple molecules, to investigate how molecular emission line flux ratios vary at different J-transitions. We confirm that ultraluminous infrared galaxies (ULIRGs) that contain or may contain luminous obscured active galactic nuclei (AGNs) tend to show higher HCN-to-HCO+ flux ratios than starburst galaxies, both at J = 4-3 and J = 3-2. For selected HCN-flux-enhanced AGN-important ULIRGs, our isotopologue H13CN, H13CO+, and HN13C J = 3-2 line observations suggest a higher abundance of HCN than HCO+ and HNC, which is interpreted to be primarily responsible for the elevated HCN flux in AGN-important galaxies. For such sources, the intrinsic HCN-to-HCO+ flux ratios after line opacity correction will be higher than the observed ratios, making the separation between AGNs and starbursts even larger. The signature of the vibrationally excited (v2 = 1f) HCN J = 4-3 emission line is seen in one ULIRG, IRAS 12112−0305 NE. P Cygni profiles are detected in the HCO+ J = 4-3 and J = 3-2 lines toward IRAS 15250+3609, with an estimated molecular outflow rate of ∼250-750 M yr−1. The SiO J = 6-5 line also exhibits a P Cygni profile in IRAS 12112+0305 NE, suggesting the presence of shocked outflow activity. Shock tracers are detected in many sources, suggesting ubiquitous shock activity in the nearby ULIRG population.
We present the results of our ALMA Cycle 4 high-spatial-resolution (0 04-0 07) observations, at HCN J = 3-2 and HCO+ J = 3-2 lines, of the nucleus of NGC 1068, the nearby prototypical type 2 active ...galactic nucleus (AGN). Our previous ALMA observations identified the compact emission of these lines at the putative location of the torus around a mass-accreting supermassive black hole. We now report that we have detected the rotation of this compact emission, with the eastern and western sides being redshifted and blueshifted, respectively. Unlike the previously reported CO J = 6-5 emission, both the morphological and dynamical alignments of the HCN J = 3-2 and HCO+ J = 3-2 emission are roughly along the east-west direction (i.e., the expected torus direction), suggesting that these molecular lines are better probes of a rotating dense molecular gas component in the torus. The western part of the torus exhibits larger velocity dispersion and stronger emission in the HCN J = 3-2 and HCO+ J = 3-2 lines than the eastern part, revealing a highly inhomogeneous molecular torus. The dense molecular gas in the torus and that of the host galaxy at 0 5-2 0 from the AGN along the torus direction are found to be counter-rotating, suggesting an external process happened in the past at the NGC 1068 nucleus.
Abstract
We report on a study of X-ray-irradiated gas in the central ∼100 pc of the Circinus galaxy, which hosts a Compton-thick active galactic nucleus (AGN), at 10 pc resolution using Chandra and ...ALMA. Based on ∼200 ks Chandra/ACIS-S data, we created an image of the Fe Kα line at 6.4 keV, tracing X-ray-irradiated dense gas. The ALMA data in Bands 6 (∼270 GHz) and 7 (∼350 GHz) cover five molecular lines: CO(J = 3–2), HCN(J = 3–2), HCN(J = 4–3), HCO+(J = 3–2), and HCO+(J = 4–3). The detailed spatial distribution of dense molecular gas was revealed, and compared to the iron line image. The molecular gas emission appeared faint in regions with bright iron emission. Motivated by this, we quantitatively discuss the possibility that the molecular gas is efficiently dissociated by AGN X-ray irradiation (i.e., creating an X-ray-dominated region). Based on a non-local thermodynamic equilibrium model, we constrained the molecular gas densities and found that they are as low as interpreted by X-ray dissociation. Furthermore, judging from inactive star formation reported in the literature, we suggest that the X-ray emission has the potential to suppress star formation, particularly in the proximity of the AGN.
We present deep near-infrared spectroscopy of six quasars at 6.1 ≤ z ≤ 6.7 with Very Large Telescope/X-Shooter and Gemini-N/GNIRS. Our objects, originally discovered through a wide-field optical ...survey with the Hyper Suprime-Cam (HSC) Subaru Strategic Program (HSC-SSP), have the lowest luminosities (−25.5 mag ≤ M1450 ≤ −23.1 mag) of the z > 5.8 quasars with measured black hole (BH) masses. From single-epoch mass measurements based on Mg ii λ2798, we find a wide range in BH masses, from MBH = 107.6 to 109.3 M . The Eddington ratios Lbol/LEdd range from 0.16 to 1.1, but the majority of the HSC quasars are powered by MBH ∼ 109 M supermassive black holes (SMBHs) accreting at sub-Eddington rates. The Eddington ratio distribution of the HSC quasars is inclined to lower accretion rates than those of Willott et al., who measured the BH masses for similarly faint z ∼ 6 quasars. This suggests that the global Eddington ratio distribution is wider than has previously been thought. The presence of MBH ∼ 109 M SMBHs at z ∼ 6 cannot be explained with constant sub-Eddington accretion from stellar remnant seed BHs. Therefore, we may be witnessing the first buildup of the most massive BHs in the first billion years of the universe, the accretion activity of which is transforming from active growth to a quiescent phase. Measurements of a larger complete sample of z 6 low-luminosity quasars, as well as deeper observations with future facilities, will enable us to better understand the early SMBH growth in the reionization epoch.
Abstract The nearby ultraluminous infrared galaxy (ULIRG) IRAS F01004−2237 exhibits 100 kpc scale continuum emission at radio wavelengths. The absence of extended X-ray emission in IRAS F01004−2237 ...has suggested an active galactic nucleus (AGN) origin for the extended radio emission, whose properties and role in merging systems still need to be better understood. We present the results of multifrequency observations of IRAS F01004−2237 conducted by the Very Long Baseline Array at 2.3 and 8.4 GHz. Compact 8.4 GHz continuum emission was detected on a 1 pc scale in the nuclear region with an intrinsic brightness temperature of 10 8.1 K suggesting that the radio source originated from an AGN, potentially driving the extended emission. In contrast, no significant emission was observed at 2.3 GHz, indicating the presence of low-frequency absorption. This absorption cannot be attributed solely to synchrotron self-absorption; alternatively, free–free absorption due to thermal plasma is mainly at work in the spectrum. From combined perspectives, including mid-infrared and X-ray data, the AGN is obscured in a dense environment. The kinetic power of the nonthermal jet, as inferred from the extended emission, can play a more important role in dispersing the surrounding medium than the thermal outflow in IRAS F01004−2237. These findings hint that jet activities in ULIRGs may contribute to AGN feedback during galaxy evolution induced by merger events.
Abstract
We report the results of Atacama Large Millimeter/submillimeter Array (ALMA) 1–2 kpc resolution, three rotational transition-line (
J
= 2–1,
J
= 3–2, and
J
= 4–3) observations of multiple ...dense molecular gas tracers (HCN, HCO
+
, and HNC) for 10 nearby (ultra)luminous infrared galaxies ((U)LIRGs). Following the matching of beam sizes to 1–2 kpc for each (U)LIRG, the high-
J
-to-low-
J
transition-line flux ratios of each molecule and the emission-line flux ratios of different molecules at each
J
transition are derived. We conduct RADEX non-LTE model calculations and find that, under a wide range of gas density and kinetic temperature, the observed HCN-to-HCO
+
flux ratios in the overall (U)LIRGs are naturally reproduced with enhanced HCN abundance compared to HCO
+
. Thereafter, molecular gas properties are constrained primarily through the use of HCN and HCO
+
data and the adoption of fiducial values for the HCO
+
column density and HCN-to-HCO
+
abundance ratio. We quantitatively confirm the following: (i) molecular gas at the (U)LIRGs’ nuclei is dense (≳10
3–4
cm
−3
) and warm (≳100 K), (ii) the molecular gas density and temperature in nine ULIRGs’ nuclei are significantly higher than those of one LIRG’s nucleus, (iii) molecular gas in starburst-dominated sources tends to be less dense and cooler than ULIRGs with luminous AGN signatures. For six selected sources, we also apply a Bayesian approach by freeing all parameters and support the above main results. Our ALMA 1–2 kpc resolution, multiple transition-line data of multiple molecules are a very powerful tool for scrutinizing the properties of molecular gas concentrated around luminous energy sources in nearby (U)LIRGs’ nuclei.
ABSTRACT We present the results of our ALMA observations of three active galactic nucleus (AGN)-dominated nuclei in optical Seyfert 1 galaxies (NGC 7469, I Zw 1, and IC 4329 A) and eleven luminous ...infrared galaxies (LIRGs) with various levels of infrared estimated energetic contributions by AGNs at the HCN and HCO+ J = 3 − 2 emission lines. The HCN and HCO+ J = 3 − 2 emission lines are clearly detected at the main nuclei of all sources, except for IC 4329 A. The vibrationally excited (v2 = 1f) HCN J = 3 − 2 and HCO+ J = 3 − 2 emission lines are simultaneously covered, and HCN v2 = 1f J = 3 − 2 emission line signatures are seen in the main nuclei of two LIRGs, IRAS 12112+0305 and IRAS 22491-1808, neither of which shows clear buried AGN signatures in the infrared. If the vibrational excitation is dominated by infrared radiative pumping, through the absorption of infrared 14 m photons, primarily originating from AGN-heated hot dust emission, then these two LIRGs may contain infrared-elusive, but (sub)millimeter-detectable, extremely deeply buried AGNs. These vibrationally excited emission lines are not detected in the three AGN-dominated optical Seyfert 1 nuclei. However, the observed HCN v2 = 1f to v = 0 flux ratios in these optical Seyferts are still consistent with the intrinsic flux ratios in LIRGs with detectable HCN v2 = 1f emission lines. The observed HCN-to-HCO+ J = 3 − 2 flux ratios tend to be higher in galactic nuclei with luminous AGN signatures compared with starburst-dominated regions, as previously seen at J = 1 − 0 and J = 4 − 3.
Abstract
We present the results of our ALMA ≲0.5 kpc resolution dense molecular line (HCN and HCO
+
J
= 2–1,
J
= 3–2, and
J
= 4–3) observations of 12 nearby (ultra)luminous infrared galaxies ...(ULIRGs). After matching beam sizes of all molecular line data to the same values in all (U)LIRGs, we derive molecular line flux ratios by extracting spectra in the central 0.5, 1, and 2 kpc circular regions and in 0.5–1 and 1–2 kpc annular regions. Based on non–local thermal equilibrium model calculations, we quantitatively confirm that the innermost (≲0.5 kpc) molecular gas is very dense (≳10
5
cm
−3
) and warm (≳300 K) in ULIRGs, and that in one LIRG, it is also modestly dense (10
4–5
cm
−3
) and warm (∼100 K). We then investigate the spatial variation of the HCN-to-HCO
+
flux ratios and high-
J
to low-
J
flux ratios of HCN and HCO
+
. A subtle sign of a decreasing trend in these ratios from the innermost (≲0.5 kpc) to the outer nuclear (0.5–2 kpc) region is discernible in a significant fraction of the observed ULIRGs. For two ULIRGS hosting an active galactic nucleus (AGN), which display the trend most clearly, we find based on a Bayesian approach that the HCN-to-HCO
+
abundance ratio and gas kinetic temperature systematically increase from the outer nuclear to the innermost region. We suggest that this trend comes from potential AGN effects because no such spatial variation is found in a starburst-dominated LIRG.
The obscuring structure surrounding active galactic nuclei (AGN) can be explained as a dust and gas flow cycle that fundamentally connects the AGN with their host galaxies. This structure is believed ...to be associated with dusty winds driven by radiation pressure. However, the role of magnetic fields, which are invoked in almost all models for accretion onto a supermassive black hole and outflows, has not been thoroughly studied. Here we report the first detection of polarized thermal emission by means of magnetically aligned dust grains in the dusty torus of NGC 1068 using ALMA Cycle 4 polarimetric dust continuum observations (0 07, 4.2 pc; 348.5 GHz, 860 m). The polarized torus has an asymmetric variation across the equatorial axis with a peak polarization of 3.7% 0.5% and position angle of 109° 2° (B-vector) at ∼8 pc east from the core. We compute synthetic polarimetric observations of magnetically aligned dust grains assuming a toroidal magnetic field and homogeneous grain alignment. We conclude that the measured 860 m continuum polarization arises from magnetically aligned dust grains in an optically thin region of the torus. The asymmetric polarization across the equatorial axis of the torus arises from (1) an inhomogeneous optical depth and (2) a variation of the velocity dispersion, i.e., a variation of the magnetic field turbulence at subparsec scales, from the eastern to the western region of the torus. These observations and modeling constrain the torus properties beyond spectral energy distribution results. This study strongly supports that magnetic fields up to a few parsecs contribute to the accretion flow onto the active nuclei.
Abstract
We present a large sample of 2.5–38
μ
m galaxy spectra drawn from a cross-archival comparison in the AKARI–Spitzer Extragalactic Spectral Survey, and investigate a subset of 113 star-forming ...galaxies with prominent polycyclic aromatic hydrocarbon (PAH) emission spanning a wide range of star formation properties. With AKARI’s extended 2.5–5
μ
m wavelength coverage, we self-consistently model for the first time
all
PAH emission bands using a modified version of
Pahfit
. We find
L
PAH
3.3
/
L
IR
∼ 0.1%, and the 3.3
μ
m PAH feature contributes ∼1.5%–3% to the total PAH power—somewhat less than earlier dust models have assumed. We establish a calibration between 3.3
μ
m PAH emission and star formation rate, but also find regimes where it loses reliability, including at high luminosity and low metallicity. The 3.4
μ
m aliphatic emission and a broad plateau feature centered at 3.47
μ
m are also modeled. As the PAH feature with the shortest wavelength, the one at 3.3
μ
m is susceptible to attenuation, leading to differences of a factor of ∼3 in the inferred star formation rate at high obscuration with different assumed attenuation geometries. Surprisingly,
L
PAH
3.3
/
L
Σ
PAH
shows no sign of decline at high luminosities, and the low-metallicity dwarf galaxy II Zw 40 exhibits an unusually strong 3.3
μ
m band; both results suggest either that the smallest PAHs are better able to survive under intense radiation fields than presumed, or that PAH emission is shifted to shorter wavelengths in intense and high-energy radiation environments. A photometric surrogate for 3.3
μ
m PAH luminosity using JWST/NIRCam is provided and found to be highly reliable at low redshift.