馬場 俊介

Recent submillimeter observations have revealed signs of pc-scale molecular inflow and atomic outflow in the nearest Seyfert 2 galaxy, the Circinus galaxy. To verify the gas kinematics suggested by these observations, we performed molecular and atomic line transfer calculations based on a physics-based 3D radiation-hydrodynamic model, which has been compared with multi-wavelength observations in this paper series. The major axis position-velocity diagram (PVD) of CO(3-2) reproduces the observed faint emission at the systemic velocity, and our calculations confirm that this component originates from failed winds falling back to the disk plane. The minor-axis PVD of [CI]($^3P_1$-$^3P_0$), when created using only the gas with positive radial velocities, presents a sign of blue- and redshifted offset peaks similar to those in the observation, suggesting that the observed peaks indeed originate from the outflow, but that the model may lack outflows as strong as those in the Circinus galaxy. Similar to the observed HCN(3-2), the similar dense gas tracer HCO$^+$(3-2) can exhibit nuclear spectra with inverse P-Cygni profiles with $\sim$0.5 pc beams, but the line shape is azimuthally dependent. The corresponding continuum absorbers are inflowing clumps at 5-10 pc from the center. To detect significant absorption with a high probability, the inclination must be fairly edge-on ($\gtrsim$85$^\circ$), and the beam size must be small ($\lesssim$1 pc). These results suggest that HCN or HCO$^+$ and [CI] lines are effective for observing pc-scale inflows and outflows, respectively.

We introduce a novel model to spectroscopically constrain the mid-infrared (MIR) extinction/attenuation curve from 3--17 um, using Polycyclic Aromatic Hydrocarbon (PAH) emission drawn from an AKARI-Spitzer extragalactic cross-archival dataset. Currently proposed MIR extinction curves vary significantly in their slopes toward the near-infrared, and the variation of the strengths and shapes of the 9.7 um and 18 um silicate absorption features make MIR spectral modeling and interpretation challenging, particularly for heavily obscured galaxies. By adopting the basic premise that PAH bands have relatively consistent intrinsic ratios within dusty starbursting galaxies, we can, for the first time, empirically determine the overall shape of the MIR attenuation curve by measuring the differential attenuation at specific PAH wavelengths. Our attenuation model shows PAH emission in most (U)LIRGs is unambiguously subjected to attenuation, and we find strong evidence that PAH bands undergo differential attenuation as obscuration increases. Compared to pre-existing results, the MIR attenuation curve derived from the model favors relatively gray continuum absorption from 3-8 $\mu$m and silicate features with intermediate strength at 9.7 um but with stronger than typical 18 um opacity.

Active galaxies contain a supermassive black hole at their center that grows by accreting matter from the surrounding galaxy. The accretion process in about the central 10 parsecs has not been directly resolved in previous observations because of the small apparent angular sizes involved. We observed the active nucleus of the Circinus Galaxy using submillimeter interferometry. A dense inflow of molecular gas was evident on subparsec scales. We calculated that less than 3% of this inflow is accreted by the black hole, with the rest being ejected by multiphase outflows, providing feedback to the host galaxy. Our observations also reveal a dense gas disk surrounding the inflow that is gravitationally unstable, which drives the accretion into about the central 1 parsec.

<jats:title>Abstract</jats:title> <jats:p>We present a catalog of the millimeter-wave (mm-wave) continuum properties of 98 nearby (<jats:italic>z</jats:italic> &lt; 0.05) active galactic nuclei (AGNs) selected from the 70 month Swift/BAT hard-X-ray catalog that have precisely determined X-ray spectral properties and subarcsecond-resolution Atacama Large Millimeter/submillimeter Array Band 6 (211–275 GHz) observations as of 2021 April. Due to the hard-X-ray (&gt;10 keV) selection, the sample is nearly unbiased for obscured systems at least up to Compton-thick-level obscuration, and provides the largest number of AGNs with high-physical-resolution mm-wave data (≲100–200 pc). Our catalog reports emission peak coordinates, spectral indices, and peak fluxes and luminosities at 1.3 mm (230 GHz). Additionally, high-resolution mm-wave images are provided. Using the images and creating radial surface brightness profiles of mm-wave emission, we identify emission extending from the central sources and isolated blob-like emission. Flags indicating the presence of these emission features are tabulated. Among 90 AGNs with significant detections of nuclear emission, 37 AGNs (≈41%) appear to have both or one of extended or blob-like components. We, in particular, investigate AGNs that show well-resolved mm-wave components and find that these seem to have a variety of origins (i.e., a jet, radio lobes, a secondary AGN, stellar clusters, a narrow-line region, galaxy disk, active star formation regions, or AGN-driven outflows), and some components have currently unclear origins.</jats:p>

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⁻³) 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 To understand the origin of nuclear (≲100 pc) millimeter-wave (mm-wave) continuum emission in active galactic nuclei (AGNs), we systematically analyzed subarcsecond resolution Band-6 (211–275 GHz) Atacama Large Millimeter/submillimeter Array data of 98 nearby AGNs (z &lt; 0.05) from the 70 month Swift/BAT catalog. The sample, almost unbiased for obscured systems, provides the largest number of AGNs to date with high mm-wave spatial resolution sampling (∼1–200 pc), and spans broad ranges of 14–150 keV luminosity {$40\lt \mathrm{log}[{L}_{14-150}/(\mathrm{erg}\,{ { \rm{s } } }^{-1})]\lt 45$}, black hole mass $[5\lt \mathrm{log}({M}_{\mathrm{BH } }/{M}_{\odot })\lt 10$], and Eddington ratio ($-4\lt \mathrm{log}{\lambda }_{\mathrm{Edd } }\lt 2$). We find a significant correlation between 1.3 mm (230 GHz) and 14–150 keV luminosities. Its scatter is ≈0.36 dex, and the mm-wave emission may serve as a good proxy of the AGN luminosity, free of dust extinction up to N_H ∼ 10²⁶ cm⁻². While the mm-wave emission could be self-absorbed synchrotron radiation around the X-ray corona according to past works, we also discuss different possible origins of the mm-wave emission: AGN-related dust emission, outflow-driven shocks, and a small-scale (&lt;200 pc) jet. The dust emission is unlikely to be dominant, as the mm-wave slope is generally flatter than expected. Also, due to no increase in the mm-wave luminosity with the Eddington ratio, a radiation-driven outflow model is possibly not the common mechanism. Furthermore, we find independence of the mm-wave luminosity on indicators of the inclination angle from the polar axis of the nuclear structure, which is inconsistent with a jet model whose luminosity depends only on the angle.

Abstract Dust-obscured galaxies (DOGs), which are observationally characterized as faint in the optical and bright in the infrared, are the final stage of galaxy mergers and are essential objects in the evolution of galaxies and active galactic nuclei (AGNs). However, the relationship between the torus-scale gas dynamics around AGNs and the DOGs’ lifetime remains unclear. We obtained the evolution of the spectral energy distributions (SEDs) of a galaxy merger system with AGN feedback from postprocessed pseudo-observations based on an N-body/smoothed particle hydrodynamics (SPH) simulation. We focused on a late-stage merger of two identical galaxies with a supermassive black hole (SMBH) of 10⁸M_⊙. We found that the infrared luminosity of the system reaches ultra- and hyperluminous infrared galaxy classes (10¹² and 10¹³L_⊙, respectively). The DOG phase corresponds to a state in which the AGNs are buried in dense gas and dust, with the infrared luminosity exceeding 3.3 × 10¹²L_⊙. We also identified subcategories of DOGs, namely bump and power-law DOGs, from the SEDs and their evolution. The bump DOGs tend to evolve to power-law DOGs over several Myrs. We found that contribution from the hot dust around the nucleus in the infrared radiation is essential for identifying the system as a power-law DOG; the gas and dust are distributed nonspherically around the nucleus, therefore, the observed properties of DOGs depend on the viewing angle. In our model, the lifetime of merger-driven DOGs is less than 4 Myr, suggesting that the observed DOG phase is a brief aspect of galaxy mergers.

Abstract A recent hydrodynamic model, the radiation-driven fountain model (Wada et al. 2016), presented a dynamical picture that active galactic nuclei (AGNs) tori sustain their geometrical thickness by gas circulation around AGNs, and previous papers have confirmed that this picture is consistent with multiwavelength observations of nearby Seyfert galaxies. Recent near-infrared observations implied that CO rovibrational absorption lines (ΔJ = ± 1, v = 0 − 1, λ ∼ 4.7 μm) could probe the physical properties of the inside tori. However, the origin of the CO absorption lines has been under debate. In this paper, we investigate the origin of the absorption lines and conditions for detecting them by performing line radiative transfer calculations based on the radiation-driven fountain model. We find that CO rovibrational absorption lines are detected at inclination angles θ_obs = 50°–80°. At the inclination angle θ_obs = 77°, we observe multi-velocity components: inflow (v_LOS = 30 km s⁻¹), systemic (v_LOS = 0 km s⁻¹), and outflows (v_LOS = −75, − 95, and −105 km s⁻¹). The inflow and outflow components (v_LOS = 30 and −95 km s⁻¹) are collisionally excited at the excitation temperatures of 186 and 380 K up to J = 12 and 4, respectively. The inflow and outflow components originate from the accreting gas on the equatorial plane at 1.5 pc from the AGN center and the outflowing gas driven by AGN radiation pressure at 1.0 pc, respectively. These results suggest that CO rovibrational absorption lines can provide us with the velocities and kinetic temperatures of the inflow and outflow in the inner few parsec region of AGN tori, and the observations can probe the gas circulation inside the tori.

Abstract The ultraluminous infrared galaxy IRAS 17208−0014 is a late-stage merger that hosts a buried active galactic nucleus (AGN). To investigate its nuclear structure, we performed high-spatial-resolution ( ∼ 0.″04 ∼ 32 pc) Atacama Large Millimeter/submillimeter Array (ALMA) observations in Band 9 (∼450 μm or ∼660 GHz), along with near-infrared AKARI spectroscopy in 2.5–5.0 μm. The Band 9 dust continuum peaks at the AGN location, and toward this position CO(J = 6 − 5) and CS(J = 14 − 13) are detected in absorption. Comparison with nonlocal thermal equilibrium calculations indicates that, within the central beam (r ∼ 20 pc), there exists a concentrated component that is dense (10⁷ cm⁻³) and warm (&gt;200 K) and has a large column density (${N}_{ { {\rm{H } } }_{2 } }\gt {10}^{23}\,{\mathrm{cm } }^{-2}$). The AKARI spectrum shows deep and broad CO rovibrational absorption at 4.67 μm. Its band profile is well reproduced with a similarly dense and large column but hotter (∼1000 K) gas. The region observed through absorption in the near-infrared is highly likely in the nuclear direction, as in the submillimeter, but with a narrower beam including a region closer to the nucleus. The central component is considered to possess a hot structure where vibrationally excited HCN emission originates. The most plausible heating source for the gas is X-rays from the AGN. The AKARI spectrum does not show other AGN signs in 2.5–4 μm, but this absence may be usual for AGNs buried in a hot mid-infrared core. Further, based on our ALMA observations, we relate the various nuclear structures of IRAS 17208−0014 that have been proposed in the literature.

Abstract We present the results of ALMA ∼2 mm, ≲1″-resolution observations of 10 (ultra)luminous infrared galaxies ([U]LIRGs; infrared luminosity ≳10^11.7L_⊙) at z &lt; 0.15, targeting dense (&gt;10⁴ cm⁻³) molecular (HCN, HCO⁺, and HNC J = 2–1) and 183 GHz H₂O 3_1,3–2_2,0 emission lines. Active galactic nucleus (AGN)-important ULIRGs tend to show higher HCN/HCO⁺J = 2–1 flux ratios than starburst-classified sources. We detect 183 GHz H₂O emission in almost all AGN-important ULIRGs, and elevated H₂O emission is found in two sources with elevated HCN J = 2–1 emission, relative to HCO⁺J = 2–1. Except one ULIRG (the Superantennae), the H₂O emission largely comes from the entire nuclear regions (∼1 kpc), rather than an AGN-origin megamaser at the very center (≪1 kpc). Nuclear (∼1 kpc) dense molecular gas mass derived from HCO⁺J = 2–1 luminosity is ≳ a few × 10⁸M_⊙, and its depletion time is estimated to be ≳10⁶ yr in all sources. Vibrationally excited J = 2–1 emission lines of HCN and HNC are detected in a few (U)LIRGs, but those of HCO⁺ are not. It is suggested that in mid-infrared-radiation-exposed innermost regions around energy sources, HCO⁺ and HNC are substantially less abundant than HCN. In our ALMA ∼2 mm data of 10 (U)LIRGs, two continuum sources are serendipitously detected within ∼10″, which are likely to be an infrared-luminous dusty galaxy at z &gt; 1 and a blazar.

Abstract To investigate the role of active galactic nucleus (AGN) X-ray irradiation on the interstellar medium (ISM), we systematically analyzed Chandra and Atacama Large Millimeter/submillimeter Array CO (J = 2–1) data for 26 hard X-ray (&gt;10 keV) selected AGNs at redshifts below 0.05. While Chandra unveils the distribution of X-ray-irradiated gas via Fe-Kα emission, the CO (J = 2–1) observations reveal that of cold molecular gas. At high resolutions ≲1″, we derive Fe-Kα and CO (J = 2–1) maps for the nuclear 2″ region and for the external annular region of 2″–4″, where 2″ is ∼100–600 pc for most of our AGNs. First, focusing on the external regions, we find the Fe-Kα emission for six AGNs above 2σ. Their large equivalent widths (≳1 keV) suggest a fluorescent process as their origin. Moreover, by comparing the 6–7 keV/3–6 keV ratio, as a proxy of Fe-Kα, and CO (J = 2–1) images for three AGNs with the highest significant Fe-Kα detections, we find a possible spatial separation. These suggest the presence of X-ray-irradiated ISM and the change in the ISM properties. Next, examining the nuclear regions, we find that (1) the 20–50 keV luminosity increases with the CO (J = 2–1) luminosity; (2) the ratio of CO (J = 2–1)/HCN (J = 1–0) luminosities increases with 20–50 keV luminosity, suggesting a decrease in the dense gas fraction with X-ray luminosity; and (3) the Fe-Kα-to-X-ray continuum luminosity ratio decreases with the molecular gas mass. This may be explained by a negative AGN feedback scenario: the mass accretion rate increases with gas mass, and simultaneously, the AGN evaporates a portion of the gas, which possibly affects star formation.

Abstract We conducted systematic observations of the H i Brα (4.05 μm) and Brβ (2.63 μm) lines in 52 nearby (z &lt; 0.3) ultraluminous infrared galaxies (ULIRGs) with AKARI. Among 33 ULIRGs wherein the lines are detected, 3 galaxies show anomalous Brβ/Brα line ratios (∼1.0), which are significantly higher than those for case B (0.565). Our observations also show that ULIRGs have a tendency to exhibit higher Brβ/Brα line ratios than those observed in Galactic H ii regions. The high Brβ/Brα line ratios cannot be explained by a combination of dust extinction and case B since dust extinction reduces the ratio. We explore possible causes for the high Brβ/Brα line ratios and show that the observed ratios can be explained by a combination of an optically thick Brα line and an optically thin Brβ line. We simulated the H ii regions in ULIRGs with the Cloudy code, and our results show that the high Brβ/Brα line ratios can be explained by high-density conditions, wherein the Brα line becomes optically thick. To achieve a column density large enough to make the Brα line optically thick within a single H ii region, the gas density must be as high as n ∼ 10⁸ cm⁻³. We therefore propose an ensemble of H ii regions, in each of which the Brα line is optically thick, to explain the high Brβ/Brα line ratio.

ABSTRACT As part of the Measuring Black Holes in below Milky Way-mass (M⋆) galaxies (MBHBM⋆) Project, we present a dynamical measurement of the supermassive black hole (SMBH) mass in the nearby lenticular galaxy NGC 3593, using cold molecular gas 12CO(2-1) emission observed at an angular resolution of ≈0${_{.}^{\prime\prime } }$3 (≈10 pc) with the Atacama Large Millimeter/submillimeter Array (ALMA). Our ALMA observations reveal a circumnuclear molecular gas disc (CND) elongated along the galaxy major axis and rotating around the SMBH. This CND has a relatively low-velocity dispersion (≲10 km s−1) and is morphologically complex, with clumps having higher integrated intensities and velocity dispersions (≲25 km s−1). These clumps are distributed along the ridges of a two-arm/bi-symmetric spiral pattern surrounded by a larger ring-like structure (radius r ≈ 10 arcsec or ≈350 pc). This pattern likely plays an important role to bridge the molecular gas reservoirs in the CND and beyond (10 ≲ r ≲ 35 arcsec or 350 pc ≲ r ≲ 1.2 kpc). Using dynamical modelling, the molecular gas kinematics allow us to infer an SMBH mass $M_{\rm BH}=2.40_{-1.05}^{+1.87}\times 10^6$ M⊙ (only statistical uncertainties at the 3σ level). We also detect a massive core of cold molecular gas (CMC) of mass MCMC = (5.4 ± 1.2) × 106 M⊙ and effective (half-mass) radius rCMC,e = 11.2 ± 2.8 pc, co-spatial with a nuclear star cluster (NSC) of mass MNSC = (1.67 ± 0.48) × 107 M⊙ and effective radius rNSC,e = 5.0 ± 1.0 pc (or 0${_{.}^{\prime\prime } }$15 ± 0${_{.}^{\prime\prime } }$03). The mass profiles of the CMC and NSC are well described by Sérsic functions with indices 1−1.4. Our MBH and MNSC estimates for NGC 3593 agree well with the recently compiled MBH–MNSC scaling relation. Although the MNSC uncertainty is twice the inferred MBH, the rapid central rise of the rotation velocities of the CND (as the radius decreases) clearly suggests an SMBH. Indeed, our dynamical models show that even if MNSC is at the upper end of its allowed range, the evidence for a BH does not vanish, but remains with a lower limit of MBH &amp;gt; 3 × 105 M⊙.

ABSTRACT We present a supermassive black hole (SMBH) mass measurement in the Seyfert 1 galaxy NGC 7469 using Atacama Large Millimeter/submillimeter Array (ALMA) observations of the atomic-[CI](1–0) and molecular-12CO(1–0) emission lines at the spatial resolution of ≈0${_{.}^{\prime\prime } }$3 (or ≈100 pc). These emissions reveal that NGC 7469 hosts a circumnuclear gas disc (CND) with a ring-like structure and a two-arm/bi-symmetric spiral pattern within it, surrounded by a starbursting ring. The CND has a relatively low σgas/V ≈ 0.35 (r ≲ 0${_{.}^{\prime\prime } }$5) and ≈0.19 (r &amp;gt; 0${_{.}^{\prime\prime } }$5), suggesting that the gas is dynamically settled and suitable for dynamically deriving the mass of its central source. As is expected from X-ray dominated region (XDR) effects that dramatically increase an atomic carbon abundance by dissociating CO molecules, we suggest that the atomic [CI](1–0) emission is a better probe of SMBH masses than CO emission in active galactic nuclei (AGNs). Our dynamical model using the [CI](1–0) kinematics yields a $M_{\rm BH}=1.78^{+2.69}_{-1.10}\times 10^7$ M⊙ and $M/L_{\rm F547M}=2.25^{+0.40}_{-0.43}$ (M⊙/L⊙). The model using the 12CO(1–0) kinematics also gives a consistent MBH with a larger uncertainty, up to an order of magnitude, i.e. $M_{\rm BH}=1.60^{+11.52}_{-1.45}\times 10^7$ M⊙. This newly dynamical MBH is ≈2 times higher than the mass determined from the reverberation mapped (RM) method using emissions arising in the unresolved broad-line region (BLR). Given this new MBH, we are able to constrain the specific RM dimensionless scaling factor of $f=7.2^{+4.2}_{-3.4}$ for the AGN BLR in NGC 7469. The gas within the unresolved BLR thus has a Keplerian virial velocity component and the inclination of $i\approx {11.0^\circ }_{-2.5}^{+2.2}$, confirming its face-on orientation in a Seyfert 1 AGN by assuming a geometrically thin BLR model.

Abstract We present Atacama Large Millimeter/submillimeter Array [C ii] 158 μm line and far-infrared (FIR) continuum emission observations toward HSC J120505.09−000027.9 (J1205−0000) at z = 6.72 with a beam size of ∼0.″8 × 0.″5 (or 4.1 kpc × 2.6 kpc), the most distant red quasar known to date. Red quasars are modestly reddened by dust and are thought to be in rapid transition from an obscured starburst to an unobscured normal quasar, driven by powerful active galactic nucleus (AGN) feedback that blows out a cocoon of interstellar medium. The FIR continuum of J1205−0000 is bright, with an estimated luminosity of L_FIR ∼ 3 × 10¹²L_⊙. The [C ii] line emission is extended on scales of r ∼ 5 kpc, greater than that of the FIR continuum. The line profiles at the extended regions are complex and broad (FWHM ∼ 630–780 km s⁻¹). Although it is not practical to identify the nature of this extended structure, possible explanations include (i) companion/merging galaxies and (ii) massive AGN-driven outflows. For the case of (i), the companions are modestly star-forming (∼10 M_⊙ yr⁻¹) but are not detected by our Subaru optical observations (y_AB,5σ = 24.4 mag). For the case of (ii), our lower limit to the cold neutral outflow rate is ∼100 M_⊙ yr⁻¹. The outflow kinetic energy and momentum are both much lower than predicted in energy-conserving wind models, suggesting that the AGN feedback in this quasar is not capable of completely suppressing its star formation.

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.