Shunsuke BABA

HSC 120505.09-000027.9 (J1205–0000) is one of the highest redshift (z = 6.72) dust-reddened quasars (red quasars) known to date. We present an improved analysis of Atacama Large Millimeter/submillimeter Array data of the [C ii] 158 μm line and the underlying rest-frame far-infrared (FIR) continuum emission, previously reported in T. Izumi et al. (2021b), toward J1205–0000. Red quasars are thought to be a transitional phase from an obscured starburst to a luminous blue quasar, in some cases associated with massive outflows driven by the active galactic nucleus (AGN). J1205–0000 has a high FIR luminosity, L _FIR = 2.5 × 10¹² L _⊙ and a total IR luminosity of L _TIR = 3.5 × 10¹² L _⊙, corresponding to a star formation rate of  ∼528 M _⊙ yr⁻¹. With the [C ii]-based dynamical mass of  ∼1 × 10¹¹ M _⊙, we conclude that J1205–0000 is hosted by a starburst galaxy. In contradiction to T. Izumi et al., our improved analysis shows no hint of a broad component in the [C ii] line spectrum. Thus there is no evidence for a host galaxy-scale fast [C ii] outflow, despite the fact that J1205–0000 has fast nuclear ionized outflows seen in the rest-frame UV. We explore several scenarios for this discrepancy (e.g., the early phase of AGN feedback, reliability of the [C ii] line as a tracer of outflows), and we claim that it is still too early to conclude that there is no significant negative AGN feedback on star formation in this red quasar.

Determining the inner structure of the molecular torus around an active galactic nucleus is essential for understanding its formation mechanism. However, spatially resolving the torus is difficult because of its small size. To probe the clump conditions in the torus, we therefore perform the systematic velocity-decomposition analyses of the gaseous ¹²CO rovibrational absorption lines (v = 0 → 1, ΔJ = ±1) at λ ∼ 4.67 μm observed toward four (ultra)luminous infrared galaxies using the high-resolution (R ∼ 5000–10,000) spectroscopy from the Subaru Telescope. We find that each transition has two to five distinct velocity components with different line-of-sight (LOS) velocities (V _LOS ∼ −240 to +100 km s⁻¹) and dispersions (σ _V ∼ 15–190 km s⁻¹), i.e., the components (a), (b), ⋯, beginning with the broadest one in each target, indicating that the tori have clumpy structures. By assuming a hydrostatic disk ( ), we find that the tori has dynamic inner structures, with the innermost component (a) outflowing with velocity ∣V _LOS∣ ∼ 160–240 km s⁻¹, and the outer components (b) and (c) outflowing more slowly or infalling with ∣V _LOS∣ ≲ 100 km s⁻¹. In addition, we find that the innermost component (a) can be attributed to collisionally excited hot (≳530 K) and dense ( ) clumps, based on the level populations. Conversely, the outer component (b) can be attributed to cold (∼30–140 K) clumps radiatively excited by a far-infrared-to-submillimeter background with a brightness temperature higher than ∼20–400 K. These observational results demonstrate the clumpy and dynamic structure of tori in the presence of background radiation.

We investigate the physical origins of the Balmer decrement anomalies in GS-NDG-9422 and RXCJ2248-ID galaxies at z ∼ 6 whose Hα/Hβ values are significantly smaller than 2.7, the latter of which also shows anomalous Hγ/Hβ and Hδ/Hβ values beyond the errors. Because the anomalous Balmer decrements are not reproduced under the Case B recombination, we explore the nebulae with optical depths smaller and larger than the Case B recombination by physical modeling. We find two cases quantitatively explaining the anomalies: (1) density-bounded nebulae that are opaque only up to around Lyγ–Ly8 transitions and (2) ionization-bounded nebulae partly/fully surrounded by optically thick excited H i clouds. The case of (1) produces more Hβ photons via Lyγ absorption in the nebulae, requiring fine tuning in optical depth values, while this case helps ionizing photon escape for cosmic reionization. The case of (2) needs the optically thick excited Hi clouds with N ₂ ≃ 10¹²−10¹³ cm⁻², where N ₂ is the column density of the hydrogen atom with the principal quantum number of n = 2. Interestingly, the high N ₂ values qualitatively agree with the recent claims for GS-NDG-9422 with the strong nebular continuum requiring a number of 2s-state electrons and for RXCJ2248-ID with the dense ionized regions likely coexisting with the optically thick clouds. While the physical origin of the optically thick excited H i clouds is unclear, these results may suggest gas clouds with excessive collisional excitation caused by an amount of accretion and supernovae in the high-z galaxies.

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.

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](³P₁–³P₀), 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 ~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 (≳85°), and the beam size must be small (≲1 pc). These results suggest that HCN or HCO⁺ and [CI] lines are effective for observing pc-scale inflows and outflows, respectively.

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.

We present a catalog of the millimeter-wave (mm-wave) continuum properties of 98 nearby (z < 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 (>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.

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⊙.

We present a new calibration for the second-order light contamination in the near-infrared grism spectroscopy with the Infrared Camera aboard AKARI, specifically for the post-cryogenic phase of the satellite (Phase 3). Following our previous work on the cryogenic phase (Phases 1 and 2), the wavelength and spectral response calibrations were revised. Unlike Phases 1 and 2, during Phase 3 the temperature of the instrument was not stable and gradually increased from 40 to 47 K. To assess the effect of the temperature increase, we divided Phase 3 into three sub-phases and performed the calibrations separately. As in Phases 1 and 2, we confirmed that there was contamination due to the wavelength dependence of the refractive index of the grism material in every sub-phase. The wavelength calibration curves for the three sub-phases coincided with each other and did not show any significant temperature dependence. The response decreased with temperature by ∼10% from the beginning to the end of Phase 3. We approximated the temperature dependence of the response at a linear relation and derived a correction factor as a function of temperature. The relative fraction of the second-order light contamination to the first-order light was found to be 25% smaller than that in Phases 1 and 2.