Showing 1–50 of 295 results for author: Nakanishi, K

Robust off-policy Reinforcement Learning via Soft Constrained Adversary

Authors: Kosuke Nakanishi, Akihiro Kubo, Yuji Yasui, Shin Ishii

Abstract: Recently, robust reinforcement learning (RL) methods against input observation have garnered significant attention and undergone rapid evolution due to RL's potential vulnerability. Although these advanced methods have achieved reasonable success, there have been two limitations when considering adversary in terms of long-term horizons. First, the mutual dependency between the policy and its corre… ▽ More Recently, robust reinforcement learning (RL) methods against input observation have garnered significant attention and undergone rapid evolution due to RL's potential vulnerability. Although these advanced methods have achieved reasonable success, there have been two limitations when considering adversary in terms of long-term horizons. First, the mutual dependency between the policy and its corresponding optimal adversary limits the development of off-policy RL algorithms; although obtaining optimal adversary should depend on the current policy, this has restricted applications to off-policy RL. Second, these methods generally assume perturbations based only on the $L_p$-norm, even when prior knowledge of the perturbation distribution in the environment is available. We here introduce another perspective on adversarial RL: an f-divergence constrained problem with the prior knowledge distribution. From this, we derive two typical attacks and their corresponding robust learning frameworks. The evaluation of robustness is conducted and the results demonstrate that our proposed methods achieve excellent performance in sample-efficient off-policy RL. △ Less

Submitted 31 August, 2024; originally announced September 2024.

Comments: 33 pages, 12 figures, 2 tables

Stable beam operation of approximately 1 mA beam under highly efficient energy recovery conditions at compact energy-recovery linac

Authors: Hiroshi Sakai, Dai Arakawa, Takaaki Furuya, Kaiichi Haga, Masayuki Hagiwara, Kentaro Harada, Yosuke Honda, Teruya Honma, Eiji Kako, Ryukou Kato, Yuuji Kojima, Taro Konomi, Hiroshi Matsumura, Taichi Miura, Takako Miura, Shinya Nagahashi, Hirotaka Nakai, Norio Nakamura, Kota Nakanishi, Kazuyuki Nigorikawa, Takashi Nogami, Takashi Obina, Feng Qiu, Hidenori Sagehashi, Shogo Sakanaka , et al. (15 additional authors not shown)

Abstract: A compact energy-recovery linac (cERL) has been un-der construction at KEK since 2009 to develop key technologies for the energy-recovery linac. The cERL began operating in 2013 to create a high-current beam with a low-emittance beam with stable continuous wave (CW) superconducting cavities. Owing to the development of critical components, such as the DC gun, superconducting cavities, and the desi… ▽ More A compact energy-recovery linac (cERL) has been un-der construction at KEK since 2009 to develop key technologies for the energy-recovery linac. The cERL began operating in 2013 to create a high-current beam with a low-emittance beam with stable continuous wave (CW) superconducting cavities. Owing to the development of critical components, such as the DC gun, superconducting cavities, and the design of ideal beam transport optics, we have successfully established approximately 1 mA stable CW operation with a small beam emittance and extremely small beam loss. This study presents the details of our key technologies and experimental results for achieving 100% energy recovery operation with extremely small beam loss during a stable, approximately 1 mA CW beam operation. △ Less

Submitted 24 August, 2024; originally announced August 2024.

Comments: 11 pages, 16 figures

Measuring 60-pc-scale Star Formation Rate of the Nearby Seyfert Galaxy NGC 1068 with ALMA, HST, VLT/MUSE, and VLA

Authors: Yuzuki Nagashima, Toshiki Saito, Soh Ikarashi, Shuro Takano, Kouichiro Nakanishi, Nanase Harada, Taku Nakajima, Akio Taniguchi, Tomoka Tosaki, Kazuharu Bamba

Abstract: Star formation rate (SFR) is a fundamental parameter for describing galaxies and inferring their evolutionary course. HII regions yield the best measure of instantaneous SFR in galaxies, although the derived SFR can have large uncertainties depending on tracers and assumptions. We present an SFR calibration for the entire molecular gas disk of the nearby Seyfert galaxy NGC 1068, based on our new h… ▽ More Star formation rate (SFR) is a fundamental parameter for describing galaxies and inferring their evolutionary course. HII regions yield the best measure of instantaneous SFR in galaxies, although the derived SFR can have large uncertainties depending on tracers and assumptions. We present an SFR calibration for the entire molecular gas disk of the nearby Seyfert galaxy NGC 1068, based on our new high-sensitivity ALMA 100GHz continuum data at 55pc (=0."8) resolution in combination with the HST Paαあるふぁ line data. In this calibration, we account for the spatial variations of dust extinction, electron temperature of HII regions, AGN contamination, and diffuse ionized gas (DIG) based on publicly available multi-wavelength data. Especially, given the extended nature and the possible non-negligible contribution to the total SFR, a careful consideration of DIG is essential. With a cross-calibration between two corrected ionized gas tracers (free-free continuum&Paαあるふぁ), the total SFR of the NGC 1068 disk is estimated to be 3.2\pm0.5 Msol/yr, one-third of the SFR without accounting for DIG (9.1\pm1.4 Msol/yr). We confirmed high SFR around the southern bar-end and the corotation radius, which is consistent with the previous SFR measurements. In addition, our total SFR exceeds the total SFR based on 8μみゅーm dust emission by a factor of 1.5. We attribute this discrepancy to the differences in the young stars at different stages of evolution traced by each tracer and their respective timescales. This study provides an example to address the various uncertainties in conventional SFR measurements and their potential to lead to significant SFR miscalculations. △ Less

Submitted 24 July, 2024; v1 submitted 22 July, 2024; originally announced July 2024.

Comments: 25 pages, 19 figures, 5 tables, accepted for publication in ApJ

An ALCHEMI inspection of sulphur-bearing species towards the central molecular zone of NGC 253

Authors: M. Bouvier, S. Viti, E. Behrens, J. Butterworth, K. -Y. Huang, J. G. Mangum, N. Harada, S. Martín, V. M. Rivilla, S. Muller, K. Sakamoto, Y. Yoshimura, K. Tanaka, K. Nakanishi, R. Herrero-Illana, L. Colzi, M. D. Gorski, C. Henkel, P. K. Humire, D. S. Meier, P. P. van der Werf, Y. T. Yan

Abstract: Sulphur-bearing species are detected in various environments within Galactic star-forming regions and are particularly abundant in the gas phase of outflows and shocks, and photo-dissociation regions. In this work, we aim to investigate the nature of the emission from the most common sulphur-bearing species observable at millimetre wavelengths towards the nuclear starburst of the galaxy NGC 253. W… ▽ More Sulphur-bearing species are detected in various environments within Galactic star-forming regions and are particularly abundant in the gas phase of outflows and shocks, and photo-dissociation regions. In this work, we aim to investigate the nature of the emission from the most common sulphur-bearing species observable at millimetre wavelengths towards the nuclear starburst of the galaxy NGC 253. We intend to understand which type of regions are probed by sulphur-bearing species and which process(es) dominate(s) the release of sulphur into the gas phase. We used the high-angular resolution (1.6" or 27 pc) observations from the ALCHEMI ALMA Large Program to image several sulphur-bearing species towards the central molecular zone (CMZ) of NGC 253. We performed local thermodynamic equilibrium (LTE) and non-LTE large velocity gradient (LVG) analyses to derive the physical conditions of the gas in which S-bearing species are emitted, and their abundance ratios across the CMZ. Finally, we compared our results with previous ALCHEMI studies and a few selected Galactic environments. We found that not all sulphur-bearing species trace the same type of gas: strong evidence indicates that H2S and part of the emission of OCS, H2CS, and SO, are tracing shocks whilst part of SO and CS emission rather trace the dense molecular gas. For some species, such as CCS and SO2, we could not firmly conclude on their origin of emission. The present analysis indicates that the emission from most sulphur-bearing species throughout the CMZ is likely dominated by shocks associated with ongoing star formation. In the inner part of the CMZ where the presence of super star clusters was previously indicated, we could not distinguish between shocks or thermal evaporation as the main process releasing the S-bearing species. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: 44 pages, 20 figures, Accepted for publication in A&A

Journal ref: A&A 689, A64 (2024)

Zone-selection effect of photoelectron intensity distributions in a nonsymmorphic system RAlSi (R : Ce and Nd)

Authors: Yusei Morita, K. Nakanishi, T. Iwata, K. Ohwada, Y. Nishioka, T. Kousa, M. Nurmamat, K. Yamagami, A. Kimura, T. Yamada, H. Tanida, Kenta Kuroda

Abstract: We investigate the electronic structures of noncentrosymmetric Weyl semimetals RAlSi (R: Ce and Nd) using soft x-ray angle-resolved photoemission spectroscopy. We find that the photoelectron intensity distribution observed in the momentum-resolved electronic bands is highly sensitive to the covered Brillouin zone (BZ) due to the zone-selection effect arising from the nonsymmorphic crystal structur… ▽ More We investigate the electronic structures of noncentrosymmetric Weyl semimetals RAlSi (R: Ce and Nd) using soft x-ray angle-resolved photoemission spectroscopy. We find that the photoelectron intensity distribution observed in the momentum-resolved electronic bands is highly sensitive to the covered Brillouin zone (BZ) due to the zone-selection effect arising from the nonsymmorphic crystal structure of RAlSi. Our data reconstruct the photoelectron distributions varied according to the zone-selection effect, and reveal comprehensive information about the electronic band structures reproduced by band calculations. This detailed information enables us to experimentally trace the Weyl-cone dispersion throughout three-dimensional momentum space, providing valuable insights into the unique properties of RAlSi. △ Less

Submitted 7 May, 2024; originally announced May 2024.

Merging gas-rich galaxies that harbor low-luminosity twin quasars at z = 6.05: a promising progenitor of the most luminous quasars

Authors: Takuma Izumi, Yoshiki Matsuoka, Masafusa Onoue, Michael A. Strauss, Hideki Umehata, John D. Silverman, Tohru Nagao, Masatoshi Imanishi, Kotaro Kohno, Yoshiki Toba, Kazushi Iwasawa, Kouichiro Nakanishi, Mahoshi Sawamura, Seiji Fujimoto, Satoshi Kikuta, Toshihiro Kawaguchi, Kentaro Aoki, Tomotsugu Goto

Abstract: We present ALMA [CII] 158 $μみゅー$m line and underlying far-infrared (FIR) continuum emission observations ($0''.57 \times 0''.46$ resolution) toward a quasar-quasar pair system recently discovered at $z = 6.05$ (Matsuoka et al. 2024). The quasar nuclei (C1 and C2) are faint ($M_{\rm 1450} \gtrsim -23$ mag), but we detect very bright [CII] emission bridging the 12 kpc between the two objects and extend… ▽ More We present ALMA [CII] 158 $μみゅー$m line and underlying far-infrared (FIR) continuum emission observations ($0''.57 \times 0''.46$ resolution) toward a quasar-quasar pair system recently discovered at $z = 6.05$ (Matsuoka et al. 2024). The quasar nuclei (C1 and C2) are faint ($M_{\rm 1450} \gtrsim -23$ mag), but we detect very bright [CII] emission bridging the 12 kpc between the two objects and extending beyond them (total luminosity $L_{\rm [CII]} \simeq 6 \times 10^9~L_\odot$). The [CII]-based total star formation rate of the system is $\sim 550~M_\odot$ yr$^{-1}$ (IR-based dust-obscured SFR is $\sim 100~M_\odot$ yr$^{-1}$), with a [CII]-based total gas mass of $\sim 10^{11}~M_\odot$. The dynamical masses of the two galaxies are large ($\sim 9 \times 10^{10}~M_\odot$ for C1 and $\sim 5 \times 10^{10}~M_\odot$ for C2). There is a smooth velocity gradient in [CII], indicating that these quasars are a tidally interacting system. We identified a dynamically distinct, fast [CII] component around C1: detailed inspection of the line spectrum there reveals the presence of a broad wing component, which we interpret as the indication of fast outflows with a velocity of $\sim 600$ km s$^{-1}$. The expected mass loading factor of the outflows, after accounting for multiphase gas, is $\gtrsim 2-3$, which is intermediate between AGN-driven and starburst-driven outflows. Hydrodynamic simulations in the literature predicted that this pair will evolve to a luminous ($M_{\rm 1450} \lesssim -26$ mag), starbursting ($\gtrsim 1000~M_\odot$ yr$^{-1}$) quasar after coalescence, one of the most extreme populations in the early universe. △ Less

Submitted 22 August, 2024; v1 submitted 3 May, 2024; originally announced May 2024.

Comments: 16 pages incl. Appendix, 4 main figures, 1 table, 2 appendix figures. Accepted for publication in the Astrophysical Journal on 2024 June 11

Physical Properties of the Southwest Outflow Streamer in the Starburst Galaxy NGC 253 with ALCHEMI

Authors: Min Bao, Nanase Harada, Kotaro Kohno, Yuki Yoshimura, Fumi Egusa, Yuri Nishimura, Kunihiko Tanaka, Kouichiro Nakanishi, Sergio Martín, Jeffrey G. Mangum, Kazushi Sakamoto, Sébastien Muller, Mathilde Bouvier, Laura Colzi, Kimberly L. Emig, David S. Meier, Christian Henkel, Pedro Humire, Ko-Yun Huang, Víctor M. Rivilla, Paul van der Werf, Serena Viti

Abstract: The physical properties of galactic molecular outflows are important as they could constrain outflow formation mechanisms. We study the properties of the southwest (SW) outflow streamer including gas kinematics, optical depth, dense gas fraction, and shock strength in the central molecular zone of the starburst galaxy NGC 253. We image the molecular emission at a spatial resolution of $\sim$27 pc… ▽ More The physical properties of galactic molecular outflows are important as they could constrain outflow formation mechanisms. We study the properties of the southwest (SW) outflow streamer including gas kinematics, optical depth, dense gas fraction, and shock strength in the central molecular zone of the starburst galaxy NGC 253. We image the molecular emission at a spatial resolution of $\sim$27 pc based on data from the ALCHEMI program. We trace the kinematics of molecular gas with CO(1-0) line. We constrain the optical depth of CO emission with CO/$^{13}$CO(1-0) ratio, the dense gas fraction with HCN/CO(1-0) ratio, as well as the shock strength with SiO(2-1)/$^{13}$CO(1-0) ratio. The CO/$^{13}$CO(1-0) integrated intensity ratio is $\sim$21 in the SW streamer region, which approximates the C/$^{13}$C isotopic abundance ratio. The higher integrated intensity ratio compared to the disk can be attributed to the optically thinner environment for CO(1-0) emission inside the SW streamer. The HCN/CO(1-0) and SiO(2-1)/$^{13}$CO(1-0) integrated intensity ratios both approach $\sim$0.2 in three giant molecular clouds (GMCs) at the base of the outflow streamers, which implies the higher dense gas fraction and enhanced strength of fast shocks in those GMCs than in the disk. The contours of those two integrated intensity ratios are extended towards the directions of outflow streamers, which connects the enhanced dense gas fraction and shock strength with molecular outflow. Moreover, the molecular gas with enhanced dense gas fraction and shock strength located at the base of the SW streamer shares the same velocity with the outflow. These phenomena suggest that the star formation inside the GMCs can trigger the shocks and further drive the molecular outflow. △ Less

Submitted 6 April, 2024; originally announced April 2024.

Comments: Accepted for publication in A&A, 14 pages, 11 figures

Sub-Parsec-Scale Jet-Driven Water Maser with Possible Gravitational Acceleration in the Radio Galaxy NGC 1052

Authors: Seiji Kameno, Yuichi Harikane, Satoko Sawada-Satoh, Tsuyoshi Sawada, Toshiki Saito, Kouichiro Nakanishi, Elizabeth Humphreys, C. M. Violette Impellizzeri

Abstract: We report sub-pc-scale observations of the 321-GHz H$_2$O emission line in the radio galaxy NGC 1052. The H$_2$O line emitter size is constrained in $< 0.6$ milliarcsec distributed on the continuum core component. The brightness temperature exceeding $10^6$ K and the intensity variation indicate certain evidence for maser emission. The maser spectrum consists of redshifted and blueshifted velocity… ▽ More We report sub-pc-scale observations of the 321-GHz H$_2$O emission line in the radio galaxy NGC 1052. The H$_2$O line emitter size is constrained in $< 0.6$ milliarcsec distributed on the continuum core component. The brightness temperature exceeding $10^6$ K and the intensity variation indicate certain evidence for maser emission. The maser spectrum consists of redshifted and blueshifted velocity components spanning $\sim 400$ km s$^{-1}$, separated by a local minimum around the systemic velocity of the galaxy. Spatial distribution of maser components show velocity gradient along the jet direction, implying that the population-inverted gas is driven by the jets interacting with the molecular torus. We identified significant change of the maser spectra between two sessions separated by 14 days. The maser profile showed a radial velocity drift of $127 \pm 13$ km s$^{-1}$ yr$^{-1}$ implying inward gravitational acceleration at 5000 Schwarzschild radii. The results demonstrate feasibility of future VLBI observations to resolve the jet-torus interacting region. △ Less

Submitted 8 February, 2024; originally announced February 2024.

Comments: 14 pages, 9 figures, accepted for publication in Publications of the Astronomical Society of Japan

Global wellposedness of general nonlinear evolution equations for distributions on the Fourier half space

Authors: Kenji Nakanishi, Baoxiang Wang

Abstract: The Cauchy problem is studied for very general systems of evolution equations, where the time derivative of solution is written by Fourier multipliers in space and analytic nonlinearity, with no other structural requirement. We construct a function space for the Fourier transform embedded in the space of distributions, and establish the global wellposedness with no size restriction. The major rest… ▽ More The Cauchy problem is studied for very general systems of evolution equations, where the time derivative of solution is written by Fourier multipliers in space and analytic nonlinearity, with no other structural requirement. We construct a function space for the Fourier transform embedded in the space of distributions, and establish the global wellposedness with no size restriction. The major restriction on the initial data is that the Fourier transform is supported on the half space, decaying at the boundary in the sense of measure. We also require uniform integrability for the orthogonal directions in the distribution sense, but no other condition. In particular, the initial data may be much more rough than the tempered distributions, and may grow polynomially at the spatial infinity. A simpler argument is also presented for the solutions locally integrable in the frequency. When the Fourier support is slightly more restricted to a conical region, the generality of equations is extremely wide, including those that are even locally illposed in the standard function spaces, such as the backward heat equations, as well as those with infinite derivatives and beyond the natural boundary of the analytic nonlinearity. As more classical examples, our results may be applied to the incompressible and compressible Navier-Stokes and Euler equations, the nonlinear diffusion and wave equations, and so on. The major drawback of the Fourier support restriction is that the solutions cannot be real valued. △ Less

Submitted 18 January, 2024; originally announced January 2024.

Comments: 52 Pages

MSC Class: 35A01; 35A02; 35B60; 35F55; 35G55; 35K55; 35L60; 35L70; 35Q30; 35Q31; 35Q55; 46F05; 46F10

The ALCHEMI atlas: principal component analysis reveals starburst evolution in NGC 253

Authors: Nanase Harada, David S. Meier, Sergio Martín, Sebastien Muller, Kazushi Sakamoto, Toshiki Saito, Mark D. Gorski, Christian Henkel, Kunihiko Tanaka, Jeffrey G. Mangum, Susanne Aalto, Rebeca Aladro, Mathilde Bouvier, Laura Colzi, Kimberly L. Emig, Rubén Herrero-Illana, Ko-Yun Huang, Kotaro Kohno, Sabine König, Kouichiro Nakanishi, Yuri Nishimura, Shuro Takano, Víctor M. Rivilla, Serena Viti, Yoshimasa Watanabe , et al. (2 additional authors not shown)

Abstract: Molecular lines are powerful diagnostics of the physical and chemical properties of the interstellar medium (ISM). These ISM properties, which affect future star formation, are expected to differ in starburst galaxies from those of more quiescent galaxies. We investigate the ISM properties in the central molecular zone of the nearby starburst galaxy NGC 253 using the ultra-wide millimeter spectral… ▽ More Molecular lines are powerful diagnostics of the physical and chemical properties of the interstellar medium (ISM). These ISM properties, which affect future star formation, are expected to differ in starburst galaxies from those of more quiescent galaxies. We investigate the ISM properties in the central molecular zone of the nearby starburst galaxy NGC 253 using the ultra-wide millimeter spectral scan survey from the ALMA Large Program ALCHEMI. We present an atlas of velocity-integrated images at a 1".6 resolution of 148 unblended transitions from 44 species, including the first extragalactic detection of HCNH$^+$ and the first interferometric images of C$_3$H$^+$, NO, HCS$^+$. We conduct a principal component analysis (PCA) on these images to extract correlated chemical species and to identify key groups of diagnostic transitions. To the best of our knowledge, our dataset is currently the largest astronomical set of molecular lines to which PCA has been applied. The PCA can categorize transitions coming from different physical components in NGC 253 such as i) young starburst tracers characterized by high-excitation transitions of HC$_3$N and complex organic molecules (COMs) versus tracers of on-going star formation (radio recombination lines) and high-excitation transitions of CCH and CN tracing PDRs, ii) tracers of cloud-collision-induced shocks (low-excitation transitions of CH$_3$OH, HNCO, HOCO$^+$, and OCS) versus shocks from star-formation-induced outflows (high-excitation transitions of SiO), as well as iii) outflows showing emission from HOC$^+$, CCH, H$_3$O$^+$, CO isotopologues, HCN, HCO$^+$, CS, and CN. Our findings show these intensities vary with galactic dynamics, star formation activities, and stellar feedback. △ Less

Submitted 4 January, 2024; originally announced January 2024.

Comments: 65 pages, 39 figures. Accepted for publication in ApJS

Dust Polarization of Prestellar and Protostellar Sources in OMC-3

Authors: Yuhua Liu, Satoko Takahashi, Masahiro Machida, Kohji Tomisaka, Josep Miquel Girart, Paul T. P. Ho, Kouichiro Nakanishi, Asako Sato

Abstract: We present the Atacama Large Millimeter/submillimeter Array (ALMA) observations of linearly polarized 1.1 mm continuum emission at $\sim$0.14" (55 au) resolution and CO ($J$=2$-$1) emission at $\sim$1.5" (590 au) resolution towards one prestellar (MMS 4), four Class 0 (MMS$\,$1, MMS$\,$3, MMS$\,$5, and MMS$\,$6), one Class I (MMS$\,$7), and one flat-spectrum (MMS$\,$2) sources in the Orion Molecul… ▽ More We present the Atacama Large Millimeter/submillimeter Array (ALMA) observations of linearly polarized 1.1 mm continuum emission at $\sim$0.14" (55 au) resolution and CO ($J$=2$-$1) emission at $\sim$1.5" (590 au) resolution towards one prestellar (MMS 4), four Class 0 (MMS$\,$1, MMS$\,$3, MMS$\,$5, and MMS$\,$6), one Class I (MMS$\,$7), and one flat-spectrum (MMS$\,$2) sources in the Orion Molecular Cloud$\,$3 region. The dust disk-like structures and clear CO outflows are detected towards all sources except for MMS$\,$4. The diameters of these disk-like structures, ranging from 16 au to 97 au, are estimated based on the deconvolved full width half maximum (FWHM) values obtained from the multi-Gaussian fitting. Polarized emissions are detected towards MMS$\,$2, MMS$\,$5, MMS$\,$6, and MMS$\,$7, while no polarized emission is detected towards MMS$\,$1, MMS$\,$3, and MMS$\,$4. MMS$\,$2, MMS$\,$5, and MMS$\,$7 show organized polarization vectors aligned with the minor axes of the disk-like structures, with mean polarization fractions ranging from 0.6$\%$ to 1.2$\%$. The strongest millimeter source, MMS$\,$6, exhibits complex polarization orientations and a remarkably high polarization fraction of $\sim$10$\%$ around the Stokes $I$ peak, and 15$-$20$\%$ on the arm-like structure, as reported by Takahashi et al. (2019). The origins of the polarized emission, such as self-scattering and dust alignment due to the magnetic field or radiative torque, are discussed for individual sources. Some disk-like sources exhibit a polarized intensity peak shift towards the nearside of the disk, which supports that the polarized emission originates from self-scattering. △ Less

Submitted 20 December, 2023; originally announced December 2023.

Comments: 46 pages, 19 figures, accepted for publication in Astrophysical Journal

Volume density structure of the NGC253 CMZ through ALCHEMI excitation analysis

Authors: Kunihiko Tanaka, Jeffrey G. Mangum, Serena Viti, Sergio Martin, Nanase Harada, Kazushi Sakamoto, Sebastien Muller, Yuki Yoshimura, Kouichiro Nakanishi, Ruben Herrero Illana, Kimberly L. Emig, S. Muhle, Hiroyuki Kaneko, Tomoka Tosaki, Erica Behrens, Victor M. Rivilla, Laura Colzi, Yuri Nishimura, P. K. Humire, Mathilde Bouvier, Ko-Yun Huang, Joshua Butterworth, David S. Meier, Paul P. van der Werf

Abstract: We present a spatially-resolved excitation analysis for the central molecular zone (CMZ) of the starburst galaxy NGC 253 using the data from the ALMA Large program ALCHEMI, whereby we explore parameters distinguishing NGC 253 from the quiescent Milky Way's Galactic Center (GC). Non-LTE analyses employing a hierarchical Bayesian framework are applied to Band 3-7 transitions from nine molecular spec… ▽ More We present a spatially-resolved excitation analysis for the central molecular zone (CMZ) of the starburst galaxy NGC 253 using the data from the ALMA Large program ALCHEMI, whereby we explore parameters distinguishing NGC 253 from the quiescent Milky Way's Galactic Center (GC). Non-LTE analyses employing a hierarchical Bayesian framework are applied to Band 3-7 transitions from nine molecular species to delineate the position-position-velocity distributions of column density ($N_\mathrm{H_2}$), volume density ($n_\mathrm{H_2}$), and temperature ($T_\mathrm{kin}$) at 27 pc resolution. Two distinct components are detected: a low-density component with $(n_\mathrm{H_2},\ T_\mathrm{kin})\sim(10^{3.3}\ \mathrm{cm}^{-3}, 85 K)$ and a high-density component with $(n_\mathrm{H_2},\ T_\mathrm{kin})\sim (10^{4.4}\ \mathrm{cm}^{-3}, 110\ \mathrm{K})$, separated at $n_\mathrm{H_2}\sim10^{3.8}\ \mathrm{cm}^{-3}$. NGC 253 has $\sim10$ times the high-density gas mass and $\sim3$ times the dense-gas mass fraction of the GC. These properties are consistent with their HCN/CO ratio but cannot alone explain the factor of $\sim30$ difference in their star formation efficiencies (SFEs), contradicting the dense-gas mass to star formation rate scaling law. The $n_\mathrm{H_2}$ histogram toward NGC 253 exhibits a shallow declining slope up to $n_\mathrm{H_2}\sim10^6\ \mathrm{cm}^{-3}$, while that of the GC steeply drops in $n_\mathrm{H_2}\gtrsim10^{4.5}\ \mathrm{cm}^{-3}$ and vanishes at $10^5\ \mathrm{cm}^{-3}$. Their dense-gas mass fraction ratio becomes consistent with their SFEs when the threshold $n_\mathrm{H_2}$ for the dense gas is taken at $\sim 10^{4.2\mbox{-}4.6}\ \mathrm{cm}^{-3}$. The rich abundance of gas above this density range in the NGC 253 CMZ, or its scarcity in the GC, is likely to be the critical difference characterizing the contrasting star formation in the centers of the two galaxies. △ Less

Submitted 20 November, 2023; originally announced November 2023.

Comments: 49 pages, 27 figures, 7 tables, accepted for publication in the Astrophysical Journal

ALMA 0.5 kpc Resolution Spatially Resolved Investigations of Nuclear Dense Molecular Gas Properties in Nearby Ultraluminous Infrared Galaxies Based on HCN and HCO$^{+}$ Three Transition Line Data

Authors: Masatoshi Imanishi, Shunsuke Baba, Kouichiro Nakanishi, Takuma Izumi

Abstract: We present the results of our ALMA $\lesssim$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 ([U]LIRGs). 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, 2 kpc circular regions, and… ▽ More We present the results of our ALMA $\lesssim$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 ([U]LIRGs). 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, 2 kpc circular regions, and 0.5-1 and 1-2 kpc annular regions. Based on non-LTE model calculations, we quantitatively confirm that the innermost ($\lesssim$0.5 kpc) molecular gas is very dense ($\gtrsim$10$^{5}$ cm$^{-3}$) and warm ($\gtrsim$300 K) in ULIRGs, and that in one LIRG is also modestly dense (10$^{4-5}$ cm$^{-3}$) and warm ($\sim$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 decreasing trend of these ratios from the innermost ($\lesssim$0.5 kpc) to outer nuclear (0.5-2 kpc) region is discernible in a significant fraction of the observed ULIRGs. For two AGN-hosting ULIRGs 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. △ Less

Submitted 27 July, 2023; originally announced July 2023.

Comments: 62 pages, 29 figures. ApJ accepted. Higher resolution version is available at http://www2.nao.ac.jp/~masaimanishi/Aph/ImanishiULIRG500pcH.pdf

Cosmic-ray-driven enhancement of the C$^0$/CO abundance ratio in W51C

Authors: Mitsuyoshi Yamagishi, Kenji Furuya, Hidetoshi Sano, Natsuko Izumi, Tatsuya Takekoshi, Hidehiro Kaneda, Kouichiro Nakanishi, Takashi Shimonishi

Abstract: We examine spatial variations of the C$^0$/CO abundance ratio ($X_{\mathrm{C/CO}}$) in the vicinity of the $γがんま$-ray supernova remnant W51C, based on [CI] ($^3P_1$-$^3P_0$), $^{12}$CO(1-0), and $^{13}$CO(1-0) observations with the ASTE and Nobeyama 45-m telescopes. We find that $X_{\mathrm{C/CO}}$ varies in a range of 0.02-0.16 (0.05 in median) inside the molecular clouds of $A_V>$100 mag, where pho… ▽ More We examine spatial variations of the C$^0$/CO abundance ratio ($X_{\mathrm{C/CO}}$) in the vicinity of the $γがんま$-ray supernova remnant W51C, based on [CI] ($^3P_1$-$^3P_0$), $^{12}$CO(1-0), and $^{13}$CO(1-0) observations with the ASTE and Nobeyama 45-m telescopes. We find that $X_{\mathrm{C/CO}}$ varies in a range of 0.02-0.16 (0.05 in median) inside the molecular clouds of $A_V>$100 mag, where photodissociation of CO by the interstellar UV is negligible. Furthermore, $X_{\mathrm{C/CO}}$ is locally enhanced up to by a factor of four near the W51C center, depending on the projected distance from the W51C center. In high-$A_V$ molecular clouds, $X_{\mathrm{C/CO}}$ is determined by the ratio of the cosmic-ray (CR) ionization rate to the H$_2$ density, and we find no clear spatial variation of the H$_2$ density against the projected distance. Hence, the high CR ionization rate may locally enhance $X_{\mathrm{C/CO}}$ near the W51C center. We also find that the observed spatial extent of the enhanced $X_{\mathrm{C/CO}}$ ($\sim$17 pc) is consistent with the diffusion distance of CRs with the energy of 100 MeV. The fact suggests that the low-energy CRs accelerated in W51C enhance $X_{\mathrm{C/CO}}$. The CR ionization rate at the $X_{\mathrm{C/CO}}$-enhanced cloud is estimated to be 3$\times$10$^{-16}$ s$^{-1}$ on the basis of time-dependent PDR simulations of $X_{\mathrm{C/CO}}$, the value of which is 30 times higher than that in the standard Galactic environment. These results demonstrate that [CI] is a powerful probe to investigate the interaction between CRs and the interstellar medium for a wide area in the vicinity of supernova remnants. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 17 pages, 8 figures, accepted for publication in PASJ

Digital quantum simulator for the time-dependent Dirac equation using discrete-time quantum walks

Authors: Shigetora Miyashita, Takahiko Satoh, Michihiko Sugawara, Naphan Benchasattabuse, Ken M. Nakanishi, Michal Hajdušek, Hyensoo Choi, Rodney Van Meter

Abstract: We introduce a quantum algorithm for simulating the time-dependent Dirac equation in 3+1 dimensions using discrete-time quantum walks. Thus far, promising quantum algorithms have been proposed to simulate quantum dynamics in non-relativistic regimes efficiently. However, only some studies have attempted to simulate relativistic dynamics due to its theoretical and computational difficulty. By lever… ▽ More We introduce a quantum algorithm for simulating the time-dependent Dirac equation in 3+1 dimensions using discrete-time quantum walks. Thus far, promising quantum algorithms have been proposed to simulate quantum dynamics in non-relativistic regimes efficiently. However, only some studies have attempted to simulate relativistic dynamics due to its theoretical and computational difficulty. By leveraging the convergence of discrete-time quantum walks to the Dirac equation, we develop a quantum spectral method that approximates smooth solutions with exponential convergence. This mitigates errors in implementing potential functions and reduces the overall gate complexity that depends on errors. We demonstrate that our approach does not require additional operations compared to the asymptotic gate complexity of non-relativistic real-space algorithms. Our findings indicate that simulating relativistic dynamics is achievable with quantum computers and can provide insights into relativistic quantum physics and chemistry. △ Less

Submitted 31 May, 2023; originally announced May 2023.

Comments: 11 pages, 5 figures

Supermassive black hole feeding and feedback observed on sub-parsec scales

Authors: Takuma Izumi, Keiichi Wada, Masatoshi Imanishi, Kouichiro Nakanishi, Kotaro Kohno, Yuki Kudoh, Taiki Kawamuro, Shunsuke Baba, Naoki Matsumoto, Yutaka Fujita, Konrad R. W. Tristram

Abstract: Active galaxies contain a supermassive black hole at their center, which grows by accreting matter from the surrounding galaxy. The accretion process in the central ~10 parsecs has not been directly resolved in previous observations, due to the small apparent angular sizes involved. We observed the active nucleus of the Circinus Galaxy using sub-millimeter interferometry. A dense inflow of molecul… ▽ More Active galaxies contain a supermassive black hole at their center, which grows by accreting matter from the surrounding galaxy. The accretion process in the central ~10 parsecs has not been directly resolved in previous observations, due to the small apparent angular sizes involved. We observed the active nucleus of the Circinus Galaxy using sub-millimeter interferometry. A dense inflow of molecular gas is evident on sub-parsec scales. We calculate 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. The observations also reveal a dense gas disk surrounding the inflow; the disk is gravitationally unstable which drives the accretion into the central ~1 parsec. △ Less

Submitted 13 November, 2023; v1 submitted 6 May, 2023; originally announced May 2023.

Comments: First release on Nov 3, 2023 in Science. 32 pages in one column = Main (13 pages, 4 figures) + Supplement (19 pages, 9 figures + 2 tables). This is the accepted version after peer review

Journal ref: Science, Volume 382, Year 2023, Issue 6670, Pages 554-559

Diverse Molecular Structures Across The Whole Star-Forming Disk of M83: High fidelity Imaging at 40pc Resolution

Authors: Jin Koda, Akihiko Hirota, Fumi Egusa, Kazushi Sakamoto, Tsuyoshi Sawada, Mark Heyer, Junichi Baba, Samuel Boissier, Daniela Calzetti, Jennifer Donovan Meyer, Bruce G. Elmegreen, Armando Gil de Paz, Nanase Harada, Luis C. Ho, Masato I. N. Kobayashi, Nario Kuno, Amanda M Lee, Barry F. Madore, Fumiya Maeda, Sergio Martin, Kazuyuki Muraoka, Kouichiro Nakanishi, Sachiko Onodera, Jorge L. Pineda, Nick Scoville , et al. (1 additional authors not shown)

Abstract: We present high-fidelity CO(1-0) imaging of molecular gas across the full star-forming disk of M83, using ALMA's 12m, 7m, and TP arrays and the MIRIAD package. The data have a mass sensitivity and resolution of 10^4Msun and 40 pc. The full disk coverage shows that the characteristics of molecular gas change radially from the center to outer disk. The molecular gas distribution shows coherent large… ▽ More We present high-fidelity CO(1-0) imaging of molecular gas across the full star-forming disk of M83, using ALMA's 12m, 7m, and TP arrays and the MIRIAD package. The data have a mass sensitivity and resolution of 10^4Msun and 40 pc. The full disk coverage shows that the characteristics of molecular gas change radially from the center to outer disk. The molecular gas distribution shows coherent large-scale structures in the inner part, including the central concentration, bar offset ridges, and prominent molecular spiral arms. In the outer disk, the spiral arms appear less spatially coherent, and even flocculent. Massive filamentary gas concentrations are abundant even in the interarm regions. Building up these structures in the interarm regions would require a very long time (~>100Myr). Instead, they must have formed within stellar spiral arms and been released into the interarm regions. For such structures to survive through the dynamical processes, the lifetimes of these structures and their constituent molecules and molecular clouds must be long (~>100Myr). These interarm structures host little or no star formation traced by Halpha. The new map also shows extended CO emission, which likely represents an ensemble of unresolved molecular clouds. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: Accepted for publication in ApJ

Dense Molecular Gas Properties of the Central Kpc of Nearby Ultraluminous Infrared Galaxies Constrained by ALMA Three Transition-line Observations

Authors: Masatoshi Imanishi, Shunsuke Baba, Kouichiro Nakanishi, Takuma Izumi

Abstract: We report the results of 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 ten 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 lin… ▽ More We report the results of 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 ten 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: (1) Molecular gas at the (U)LIRGs' nuclei is dense ($\gtrsim$10$^{3-4}$ cm$^{-3}$) and warm ($\gtrsim$100 K). (2) Molecular gas density and temperature in nine ULIRGs' nuclei are significantly higher than that of one LIRG's nucleus. (3) 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. △ Less

Submitted 14 March, 2023; originally announced March 2023.

Comments: 42 pages, 13 figures. ApJ accepted

ALMA Detection of 321 GHz water maser emission in the radio galaxy NGC 1052

Authors: Seiji Kameno, Yuichi Harikane, Satoko Sawada-Satoh, Tsuyoshi Sawada, Toshiki Saito, Kouichiro Nakanishi, Elizabeth Humphreys

Abstract: The Atacama Large Millimeter/submillimeter Array (ALMA) serendipitously detected H$_2$O $J_{Ka, Kc} = 10_{2,9} - 9_{3,6}$ emission at 321 GHz in NGC 1052. This is the first submillimeter maser detection in a radio galaxy and the most luminous 321-GHz H$_2$O maser known to date with the isotropic luminosity of 1090 $L_{\odot}$. The line profile consists of a broad velocity component with FWHM… ▽ More The Atacama Large Millimeter/submillimeter Array (ALMA) serendipitously detected H$_2$O $J_{Ka, Kc} = 10_{2,9} - 9_{3,6}$ emission at 321 GHz in NGC 1052. This is the first submillimeter maser detection in a radio galaxy and the most luminous 321-GHz H$_2$O maser known to date with the isotropic luminosity of 1090 $L_{\odot}$. The line profile consists of a broad velocity component with FWHM $= 208 \pm 12$ km s$^{-1}$ straddling the systemic velocity and a narrow component with FWHM $= 44 \pm 3$ km s$^{-1}$ blueshifted by 160 km s$^{-1}$. The profile is significantly different from the known 22-GHz $6_{1,6} - 5_{2,3}$ maser which shows a broad profile redshifted by 193 km s$^{-1}$. The submillimeter maser is spatially unresolved with a synthesized beam of $0^{\prime \prime}.68 \times 0^{\prime \prime}.56$ and coincides with the continuum core position within 12 pc. These results indicate amplification of the continuum emission through high-temperature ($>1000$ K) and dense ($n({\rm H}_2{\rm O}) > 10^4$ cm$^{-3}$) molecular gas in front of the core △ Less

Submitted 8 February, 2023; originally announced February 2023.

Comments: 3 figures, Accepted for publication in Publications of Astronomical Society of Japan

CO($J$=1-0) mapping survey of 64 galaxies in the Fornax cluster with the ALMA Morita array

Authors: Kana Morokuma-Matsui, Kenji Bekki, Jing Wang, Paolo Serra, Yusei Koyama, Tomoki Morokuma, Fumi Egusa, Bi-Qing For, Kouichiro Nakanishi, Bäbel S. Koribalski, Takashi Okamoto, Tadayuki Kodama, Bumhyun Lee, Filippo M. Maccagni, Rie E. Miura, Daniel Espada, Tsutomu T. Takeuchi, Dong Yang, Minju M. Lee, Masaki Ueda, Kyoko Matsushita

Abstract: We conduct a $^{12}$C$^{16}$O($J$=1-0) (hereafter CO) mapping survey of 64 galaxies in the Fornax cluster using the ALMA Morita array in cycle 5. CO emission is detected from 23 out of the 64 galaxies. Our sample includes dwarf, spiral and elliptical galaxies with stellar masses of $M_{\rm star}\sim10^{6.3-11.6}$~M$_\odot$. The achieved beam size and sensitivity are $15''\times8''$ and $\sim12$~mJ… ▽ More We conduct a $^{12}$C$^{16}$O($J$=1-0) (hereafter CO) mapping survey of 64 galaxies in the Fornax cluster using the ALMA Morita array in cycle 5. CO emission is detected from 23 out of the 64 galaxies. Our sample includes dwarf, spiral and elliptical galaxies with stellar masses of $M_{\rm star}\sim10^{6.3-11.6}$~M$_\odot$. The achieved beam size and sensitivity are $15''\times8''$ and $\sim12$~mJy~beam$^{-1}$ at the velocity resolution of $\sim10$~km~s$^{-1}$, respectively. We study the cold-gas (molecular- and atomic-gas) properties of 38 subsamples with $M_{\rm star}>10^9$~M$_\odot$ combined with literature HI data. We find that: (1) the low star-formation (SF) activity in the Fornax galaxies is caused by the decrease in the cold-gas mass fraction with respect to stellar mass (hereafter, gas fraction) rather than the decrease of the SF efficiency from the cold gas; (2) the atomic-gas fraction is more heavily reduced than the molecular-gas fraction of such galaxies with low SF activity. A comparison between the cold-gas properties of the Fornax galaxies and their environmental properties suggests that the atomic gas is stripped tidally and by the ram pressure, which leads to the molecular gas depletion with an aid of the strangulation and consequently SF quenching. Pre-processes in the group environment would also play a role in reducing cold-gas reservoirs in some Fornax galaxies. △ Less

Submitted 16 October, 2022; originally announced October 2022.

Comments: 53 pages, 41 figures, accepted for publication in ApJS

Optimization of RF phase and beam loading distribution among RF stations in SuperKEKB

Authors: S. Ogasawara, K. Akai, T. Kobayashi, K. Nakanishi, M. Nishiwaki

Abstract: SuperKEKB is the e-/e+ collider which targets the world highest luminosity. In recent operation, SuperKEKB achieved a new world record of $4.71\times 10^{34} cm^{-2}s^{-1}$ for luminosity with the beam current of 1.4 A. In the future, the beam current will be increased further to aim at the design value of 3.6 A and much higher luminosity. The RF system consists of 38 cavities (30 klystron station… ▽ More SuperKEKB is the e-/e+ collider which targets the world highest luminosity. In recent operation, SuperKEKB achieved a new world record of $4.71\times 10^{34} cm^{-2}s^{-1}$ for luminosity with the beam current of 1.4 A. In the future, the beam current will be increased further to aim at the design value of 3.6 A and much higher luminosity. The RF system consists of 38 cavities (30 klystron stations), which share the huge beam loading brought by high current beam with each other cavities. For beam stability and power efficiency, it is important to distribute the beam loading properly among RF cavities. It is equivalent to adjust the acceleration phase of each cavity. However, it is difficult to evaluate the acceleration phase using only the pickup signal. Therefore, we established a method to evaluate the beam loading balance among RF stations from the RF power measurement for each cavity, and to adjust the acceleration phase. This report introduces the method for evaluating and optimizing the beam loading (acceleration phase) among stations in SuperKEKB, which has a large number of RF stations, and its operation. △ Less

Submitted 10 October, 2022; originally announced October 2022.

Comments: Talk presented at LLRF Workshop 2022 (LLRF2022, arXiv:2208.13680)

Report number: LLRF2022/8

Tracing Interstellar Heating: An ALCHEMI Measurement of the HCN Isomers in NGC 253

Authors: Erica Behrens, Jeffrey G. Mangum, Jonathan Holdship, Serena Viti, Nanase Harada, Sergio Martin, Kazushi Sakamoto, Sebastien Muller, Kunihiko Tanaka, Kouichiro Nakanishi, Ruben Herrero-Illana, Yuki Yoshimura, Rebeca Aladro, Laura Colzi, Kimberly L. Emig, Christian Henkel, Ko-Yun Huang, P. K. Humire, David S. Meier, Victor M. Rivilla

Abstract: We analyze HCN and HNC emission in the nearby starburst galaxy NGC 253 to investigate its effectiveness in tracing heating processes associated with star formation. This study uses multiple HCN and HNC rotational transitions observed using ALMA via the ALCHEMI Large Program. To understand the conditions and associated heating mechanisms within NGC 253's dense gas, we employ Bayesian nested samplin… ▽ More We analyze HCN and HNC emission in the nearby starburst galaxy NGC 253 to investigate its effectiveness in tracing heating processes associated with star formation. This study uses multiple HCN and HNC rotational transitions observed using ALMA via the ALCHEMI Large Program. To understand the conditions and associated heating mechanisms within NGC 253's dense gas, we employ Bayesian nested sampling techniques applied to chemical and radiative transfer models which are constrained using our HCN and HNC measurements. We find that the volume density $n_{\text{H}_{2}}$ and cosmic ray ionization rate (CRIR) $ζぜーた$ are enhanced by about an order of magnitude in the galaxy's central regions as compared to those further from the nucleus. In NGC 253's central GMCs, where observed HCN/HNC abundance ratios are lowest, $n \sim 10^{5.5}$ cm$^{-3}$ and $ζぜーた\sim 10^{-12}$ s$^{-1}$ (greater than $10^4$ times the average Galactic rate). We find a positive correlation in the association of both density and CRIR with the number of star formation-related heating sources (supernova remnants, HII regions, and super hot cores) located in each GMC, as well as a correlation between CRIRs and supernova rates. Additionally, we see an anticorrelation between the HCN/HNC ratio and CRIR, indicating that this ratio will be lower in regions where $ζぜーた$ is higher. Though previous studies suggested HCN and HNC may reveal strong mechanical heating processes in NGC 253's CMZ, we find cosmic ray heating dominates the heating budget, and mechanical heating does not play a significant role in the HCN and HNC chemistry. △ Less

Submitted 8 November, 2022; v1 submitted 13 September, 2022; originally announced September 2022.

Comments: 33 pages, 23 figures, accepted for publication by the Astrophysical Journal

ALCHEMI finds a "shocking" carbon footprint in the starburst galaxy NGC~253

Authors: Nanase Harada, Sergio Martin, Jeff Mangum, Kazushi Sakamoto, Sebastian Muller, Victor Rivilla, Christian Henkel, David Meier, Laura Colzi, Mitsuyoshi Yamagishi, Kunihiko Tanaka, Kouichiro Nakanishi, Ruben Herrero-Illana, Yuki Yoshimura, Pedro Humire, Rebeca Aladro, Paul van der Werf, Kim Emig

Abstract: Centers of starburst galaxies may be characterized by a specific gas and ice chemistry due to their gas dynamics and the presence of various ice desorption mechanisms. This may result in a peculiar observable composition. We analyze abundances of $CO_2$, a reliable tracer of ice chemistry, from data collected as part of the ALMA large program ALCHEMI, a wide-frequency spectral scan toward the star… ▽ More Centers of starburst galaxies may be characterized by a specific gas and ice chemistry due to their gas dynamics and the presence of various ice desorption mechanisms. This may result in a peculiar observable composition. We analyze abundances of $CO_2$, a reliable tracer of ice chemistry, from data collected as part of the ALMA large program ALCHEMI, a wide-frequency spectral scan toward the starburst galaxy NGC~253 with an angular resolution of 1.6$''$. We constrain the $CO_2$ abundances in the gas phase using its protonated form $HOCO^+$. The distribution of $HOCO^+$ is similar to that of methanol, which suggests that $HOCO^+$ is indeed produced from the protonation of $CO_2$ sublimated from ice. The $HOCO^+$ fractional abundances are found to be $(1-2)\times10^{-9}$ at the outer part of the central molecular zone (CMZ), while they are lower ($\sim10^{-10}$) near the kinematic center. This peak fractional abundance at the outer CMZ is comparable to that in the Milky Way CMZ, and orders of magnitude higher than that in Galactic disk star-forming regions. From the range of $HOCO^+/CO_2$ ratios suggested from chemical models, the gas-phase $CO_2$ fractional abundance is estimated to be $(1-20)\times10^{-7}$ at the outer CMZ, and orders of magnitude lower near the center. We estimate the $CO_2$ ice fractional abundances at the outer CMZ to be $(2-5)\times10^{-6}$ from the literature. A comparison between the ice and gas $CO_2$ abundances suggests an efficient sublimation mechanism. This sublimation is attributed to large-scale shocks at the orbital intersections of the bar and CMZ. △ Less

Submitted 30 August, 2022; originally announced August 2022.

Comments: 22 pages, 9 figures. Accepted for publication in the Astrophysical Journal

An AGN with an ionized gas outflow in a massive quiescent galaxy in a protocluster at $\bf z=3.09$

Authors: Mariko Kubo, Hideki Umehata, Yuichi Matsuda, Masaru Kajisawa, Charles C. Steidel, Toru Yamada, Ichi Tanaka, Bunyo Hatsukade, Yoichi Tamura, Kouichiro Nakanishi, Kotaro Kohno, Kianhong Lee, Keiichi Matsuda, Yiping Ao, Tohru Nagao, Min S. Yun

Abstract: We report the detection of an ionized gas outflow from an $X$-ray active galactic nucleus (AGN) hosted in a massive quiescent galaxy in a protocluster at $z=3.09$ (J221737.29+001823.4). It is a type-2 QSO with broad ($W_{80}>1000$ km s$^{-1}$) and strong ($\log (L_{\rm [OIII]}$ / erg s$^{-1})\approx43.4$) [O {\footnotesize III}]$λらむだλらむだ$4959,5007 emission lines detected by slit spectroscopy in three-po… ▽ More We report the detection of an ionized gas outflow from an $X$-ray active galactic nucleus (AGN) hosted in a massive quiescent galaxy in a protocluster at $z=3.09$ (J221737.29+001823.4). It is a type-2 QSO with broad ($W_{80}>1000$ km s$^{-1}$) and strong ($\log (L_{\rm [OIII]}$ / erg s$^{-1})\approx43.4$) [O {\footnotesize III}]$λらむだλらむだ$4959,5007 emission lines detected by slit spectroscopy in three-position angles using Multi-Object Infra-Red Camera and Spectrograph (MOIRCS) on the Subaru telescope and the Multi-Object Spectrometer For Infra-Red Exploration (MOSFIRE) on the Keck-I telescope. In the all slit directions, [O {\footnotesize III}] emission is extended to $\sim15$ physical kpc and indicates a powerful outflow spreading over the host galaxy. The inferred ionized gas mass outflow rate is $\rm 22\pm3~M_{\odot}~yr^{-1}$. Although it is a radio source, according to the line diagnostics using H$βべーた$, [O {\footnotesize II}], and [O {\footnotesize III}], photoionization by the central QSO is likely the dominant ionization mechanism rather than shocks caused by radio jets. On the other hand, the spectral energy distribution of the host galaxy is well characterized as a quiescent galaxy that has shut down star formation by several hundred Myr ago. Our results suggest a scenario that QSOs are powered after the shut-down of the star formation and help to complete the quenching of massive quiescent galaxies at high redshift. △ Less

Submitted 7 July, 2022; originally announced July 2022.

Comments: 18 pages 6 figures, accepted for publication in ApJ

Methanol masers in NGC 253 with ALCHEMI

Authors: P. K. Humire, C. Henkel, A. Hernández-Gómez, S. Martín, J. Mangum, N. Harada, S. Muller, K. Sakamoto, K. Tanaka, Y. Yoshimura, K. Nakanishi, S. Mühle, R. Herrero-Illana, D. S. Meier, E. Caux, R. Aladro, R. Mauersberger, S. Viti, L. Colzi, V. M. Rivilla, M. Gorski, K. M. Menten, K. -Y. Huang, S. Aalto, P. P. van der Werf , et al. (1 additional authors not shown)

Abstract: Context: Methanol masers of Class I (collisionally-pumped) and Class II (radiatively-pumped) have been studied in great detail in our Galaxy in a variety of astrophysical environments such as shocks and star-forming regions and are helpful to analyze the properties of the dense interstellar medium. However, the study of methanol masers in external galaxies is still in its infancy. Aims: Our main g… ▽ More Context: Methanol masers of Class I (collisionally-pumped) and Class II (radiatively-pumped) have been studied in great detail in our Galaxy in a variety of astrophysical environments such as shocks and star-forming regions and are helpful to analyze the properties of the dense interstellar medium. However, the study of methanol masers in external galaxies is still in its infancy. Aims: Our main goal is to search for methanol masers in the central molecular zone (CMZ; inner 500 pc) of the nearby starburst galaxy NGC 253. Methods: Covering a frequency range between 84 and 373 GHz ($λらむだ$ = 3.6 to 0.8 mm) at high angular (1.6"$\sim$27 pc) and spectral ($\sim$8--9 km s$^{-1}$) resolution with the ALMA large program ALCHEMI, we have probed different regions across the CMZ of NGC 253. In order to look for methanol maser candidates, we employed the rotation diagram method and a set of radiative transfer models. Results: We detect for the first time masers above 84 GHz in NGC 253, covering an ample portion of the $J_{-1}\rightarrow(J-$ 1)$_{0}-E$ line series (at 84, 132, 229, and 278 GHz) and the $J_{0}\rightarrow(J-$ 1)$_{1}-A$ series (at 95, 146, and 198 GHz). This confirms the presence of the Class I maser line at 84 GHz, already reported but now being detected in more than one location. For the $J_{-1}\rightarrow(J-$ 1)$_{0}-E$ line series, we observe a lack of Class I maser candidates in the central star-forming disk. Conclusions: The physical conditions for maser excitation in the $J_{-1}\rightarrow(J-$ 1)$_{0}-E$ line series can be weak shocks and cloud-cloud collisions as suggested by shock tracers (SiO and HNCO) in bi-symmetric shock/active regions located in the outskirts of the CMZ. On the other hand, the presence of photodissociation regions due to a high star-formation rate would be needed to explain the lack of Class I masers in the very central regions. △ Less

Submitted 6 May, 2022; originally announced May 2022.

Comments: Accepted for publication in A&A. 29 pages, 17 figures (4 in Appendix)

Energizing Star Formation: The Cosmic Ray Ionization Rate in NGC 253 Derived From ALCHEMI Measurements of H$_3$O$^+$ and SO

Authors: Jonathan Holdship, Jeffrey G. Mangum, Serena Viti, Erica Behrens, Nanase Harada, Sergio Martín, Kazushi Sakamoto, Sebastien Muller, Kunihiko Tanaka, Kouichiro Nakanishi, Rubén Herrero-Illana, Yuki Yoshimura, Rebeca Aladro, Laura Colzi, Kimberly L. Emig, Christian Henkel, Yuri Nishimura, Víctor M. Rivilla, Paul P. van der Werf

Abstract: The cosmic ray ionization rate (CRIR) is a key parameter in understanding the physical and chemical processes in the interstellar medium. Cosmic rays are a significant source of energy in star formation regions, which impacts the physical and chemical processes which drive the formation of stars. Previous studies of the circum-molecular zone (CMZ) of the starburst galaxy NGC 253 have found evidenc… ▽ More The cosmic ray ionization rate (CRIR) is a key parameter in understanding the physical and chemical processes in the interstellar medium. Cosmic rays are a significant source of energy in star formation regions, which impacts the physical and chemical processes which drive the formation of stars. Previous studies of the circum-molecular zone (CMZ) of the starburst galaxy NGC 253 have found evidence for a high CRIR value; $10^3-10^6$ times the average cosmic ray ionization rate within the Milky Way. This is a broad constraint and one goal of this study is to determine this value with much higher precision. We exploit ALMA observations towards the central molecular zone of NGC 253 to measure the CRIR. We first demonstrate that the abundance ratio of H$_3$O$^+$ and SO is strongly sensitive to the CRIR. We then combine chemical and radiative transfer models with nested sampling to infer the gas properties and CRIR of several star-forming regions in NGC 253 due to emission from their transitions. We find that each of the four regions modelled has a CRIR in the range $(1-80)\times10^{-14}$ s$^{-1}$ and that this result adequately fits the abundances of other species that are believed to be sensitive to cosmic rays including C$_2$H, HCO$^+$, HOC$^+$, and CO. From shock and PDR/XDR models, we further find that neither UV/X-ray driven nor shock dominated chemistry are a viable single alternative as none of these processes can adequately fit the abundances of all of these species. △ Less

Submitted 11 May, 2022; v1 submitted 7 April, 2022; originally announced April 2022.

Comments: 24 pages, 15 figures, accepted for publication in ApJ, updated to include changes made for proofs

Multi-wavelength and Multi-CO View of The Minor Merger Driven Star Formation in the Nearby LIRG NGC 3110

Authors: Yuka Kawana, Toshiki Saito, Sachiko K. Okumura, Ryohei Kawabe, Daniel Espada, Daisuke Iono, Hiroyuki Kaneko, Minju M. Lee, Tomonari Michiyama, Kentaro Motohara, Kouichiro Nakanishi, Alex R. Pettitt, Zara Randriamanakoto, Junko Ueda, Takuji Yamashita

Abstract: We present Atacama Large Millimeter/submillimeter Array observations of multiple CO(1-0), $^{13}$CO(1-0), and C$^{18}$O(1-0) lines and 2.9 mm and 1.3 mm continuum emission toward the nearby interacting luminous infrared galaxy NGC 3110, supplemented with similar spatial resolution H$αあるふぁ$, 1.4GHz continuum, and $K$-band data. We estimate the typical CO-to-H$_2$ conversion factor of 1.7 $M_{\odot}$ (K… ▽ More We present Atacama Large Millimeter/submillimeter Array observations of multiple CO(1-0), $^{13}$CO(1-0), and C$^{18}$O(1-0) lines and 2.9 mm and 1.3 mm continuum emission toward the nearby interacting luminous infrared galaxy NGC 3110, supplemented with similar spatial resolution H$αあるふぁ$, 1.4GHz continuum, and $K$-band data. We estimate the typical CO-to-H$_2$ conversion factor of 1.7 $M_{\odot}$ (K km s$^{-1}$ pc$^2$)$^{-1}$ within the disk using LTE-based and dust-based H$_2$ column densities, and measure the 1-kpc scale surface densities of star formation rate ($Σしぐま_{\rm SFR}$), super star clusters ($Σしぐま_{\rm SSC}$), molecular gas mass, and star formation efficiency (SFE) toward the entire gas disk. These parameters show a peak at the southern part of the southern spiral arm (SFE $\sim$ 10$^{-8.2}$ yr$^{-1}$, $Σしぐま_{\rm SFR}$ $\sim$ 10$^{-0.6}$ $M_{\odot}$ kpc$^{-2}$ yr$^{-1}$, $Σしぐま_{\rm SSC}$ $\sim$ 6.0 kpc$^{-2}$), which is likely attributed to the on-going tidal interaction with the companion galaxy MCG-01-26-013, as well as toward the circumnuclear region. We also find that thermal free-free emission contributes to a significant fraction of the millimeter continuum emission at the southern peak position. Those measurements imply that the peak of the southern arm is an active and young star-forming region, whereas the central part of NGC 3110 is a site of long-continued star formation. We suggest that, during the early stage of the galaxy-galaxy interaction with large mass ratio that in NGC 3110, fragmentation along the main galaxy's arms is an important driver of merger-induced star formation and massive gas inflow results in dusty nuclear starbursts. △ Less

Submitted 19 March, 2022; originally announced March 2022.

Comments: 17 pages, 9 figures, 2 tables, accepted for publication in The Astrophysical Journal

High Dimensional Statistical Analysis and its Application to ALMA Map of NGC 253

Authors: Tsutomu T. Takeuchi, Kazuyoshi Yata, Kento Egashira, Makoto Aoshima, Aki Ishii, Suchetha Cooray, Kouichiro Nakanishi, Kotaro Kohno, Kai T. Kono

Abstract: In astronomy, if we denote the dimension of data as $d$ and the number of samples as $n$, we often meet a case with $n \ll d$. Traditionally, such a situation is regarded as ill-posed, and there was no choice but to throw away most of the information in data dimension to let $d < n$. The data with $n \ll d$ is referred to as high-dimensional low sample size (HDLSS). {}To deal with HDLSS problems,… ▽ More In astronomy, if we denote the dimension of data as $d$ and the number of samples as $n$, we often meet a case with $n \ll d$. Traditionally, such a situation is regarded as ill-posed, and there was no choice but to throw away most of the information in data dimension to let $d < n$. The data with $n \ll d$ is referred to as high-dimensional low sample size (HDLSS). {}To deal with HDLSS problems, a method called high-dimensional statistics has been developed rapidly in the last decade. In this work, we first introduce the high-dimensional statistical analysis to the astronomical community. We apply two representative methods in the high-dimensional statistical analysis methods, the noise-reduction principal component analysis (NRPCA) and regularized principal component analysis (RPCA), to a spectroscopic map of a nearby archetype starburst galaxy NGC 253 taken by the Atacama Large Millimeter/Submillimeter Array (ALMA). The ALMA map is a typical HDLSS dataset. First we analyzed the original data including the Doppler shift due to the systemic rotation. The high-dimensional PCA could describe the spatial structure of the rotation precisely. We then applied to the Doppler-shift corrected data to analyze more subtle spectral features. The NRPCA and RPCA could quantify the very complicated characteristics of the ALMA spectra. Particularly, we could extract the information of the global outflow from the center of NGC 253. This method can also be applied not only to spectroscopic survey data, but also any type of data with small sample size and large dimension. △ Less

Submitted 20 February, 2024; v1 submitted 9 March, 2022; originally announced March 2022.

Comments: 33 pages, 21 figures, accepted for publication in ApJS (Jan. 31, 2024)

Detection of nitrogen and oxygen in a galaxy at the end of reionization

Authors: Ken-ichi Tadaki, Akiyoshi Tsujita, Yoichi Tamura, Kotaro Kohno, Bunyo Hatsukade, Daisuke Iono, Minju M. Lee, Yuichi Matsuda, Tomonari Michiyama, Tohru Nagao, Kouichiro Nakanishi, Yuri Nishimura, Toshiki Saito, Hideki Umehata, Jorge Zavala

Abstract: We present observations of [NII] 205 $μみゅー$m, [OIII] 88 $μみゅー$m and dust emission in a strongly-lensed, submillimeter galaxy (SMG) at $z=6.0$, G09.83808, with the Atacama Large Millimeter/submillimeter Array (ALMA). Both [NII] and [OIII] line emissions are detected at $>12σしぐま$ in the 0.8$"$-resolution maps. Lens modeling indicates that the spatial distribution of the dust continuum emission is well charac… ▽ More We present observations of [NII] 205 $μみゅー$m, [OIII] 88 $μみゅー$m and dust emission in a strongly-lensed, submillimeter galaxy (SMG) at $z=6.0$, G09.83808, with the Atacama Large Millimeter/submillimeter Array (ALMA). Both [NII] and [OIII] line emissions are detected at $>12σしぐま$ in the 0.8$"$-resolution maps. Lens modeling indicates that the spatial distribution of the dust continuum emission is well characterized by a compact disk with an effective radius of 0.64$\pm$0.02 kpc and a high infrared surface brightness of $Σしぐま_\mathrm{IR}=(1.8\pm0.3)\times10^{12}~L_\odot$ kpc$^{-2}$. This result supports that G09.83808 is the progenitors of compact quiescent galaxies at $z\sim4$, where the majority of its stars are expected to be formed through a strong and short burst of star formation. G09.83808 and other lensed SMGs show a decreasing trend of the [NII] line to infrared luminosity ratio with increasing continuum flux density ratio between 63 $μみゅー$m and 158 $μみゅー$m, as seen in local luminous infrared galaxies (LIRGs). The decreasing trend can be reproduced by photoionization models with increasing ionization parameters. Furthermore, by combining the [NII]/[OIII] luminosity ratio with far-infrared continuum flux density ratio in G09.83808, we infer that the gas phase metallicity is already $Z\approx 0.5-0.7~Z_\odot$. G09.83808 is likely one of the earliest galaxies that has been chemically enriched at the end of reionization. △ Less

Submitted 20 February, 2022; originally announced February 2022.

Comments: 8 pages, 5 figures, accepted for publication in PASJ

Properties of Molecular Gas in Galaxies in Early and Mid Stages of Interaction. III. Resolved Kennicutt-Schmidt Law

Authors: Hiroyuki Kaneko, Nario Kuno, Daisuke Iono, Yoichi Tamura, Tomoka Tosaki, Koichiro Nakanishi, Tsuyoshi Sawada

Abstract: We study properties of the interstellar medium, an ingredient of stars, and star formation activity, in four nearby galaxy pairs in the early and mid stages of interaction for both a galaxy scale and a kpc scale. The galaxy-scale Kennicutt-Schmidt law shows that seven of eight interacting galaxies have a star formation rate within a factor of three compared with the best-fit of the isolated galaxi… ▽ More We study properties of the interstellar medium, an ingredient of stars, and star formation activity, in four nearby galaxy pairs in the early and mid stages of interaction for both a galaxy scale and a kpc scale. The galaxy-scale Kennicutt-Schmidt law shows that seven of eight interacting galaxies have a star formation rate within a factor of three compared with the best-fit of the isolated galaxies, although we have shown that molecular hydrogen gas is efficiently produced from atomic hydrogen during the interaction in the previous paper. The galaxy-scale specific star formation rate (sSFR) and star formation efficiency (SFE) in interacting galaxies are comparable to those in isolated galaxies. We also investigate SFE and the Kennicutt-Schmidt law on a kpc scale. The spatial distributions of SFE reveal that SFE is locally enhanced, and the enhanced regions take place asymmetrically or at off-centre regions. The local enhancement of SFE could be induced by shock. We find that the index of the Kennicutt-Schmidt law for the interacting galaxies in the early stage is 1.30$\pm$0.04, which is consistent with that of the isolated galaxies. Since CO emission, which is used in the Kennicutt-Schmidt law, is a tracer of the amount of molecular gas, this fact suggests that dense gas, which is more directly connected to star formation, is not changed at the early stage of interaction. △ Less

Submitted 6 January, 2022; originally announced January 2022.

Comments: 40 pages, 15 figures, 10 tables; accepted for publication in PASJ

First extragalactic detection of a phosphorus-bearing molecule with ALCHEMI: phosphorus nitride (PN)

Authors: D. Haasler, V. M. Rivilla, S. Martín, J. Holdship, S. Viti, N. Harada, J. Mangum, K. Sakamoto, S. Muller, K. Tanaka, Y. Yoshimura, K. Nakanishi, L. Colzi, L. Hunt, K. L. Emig, R. Aladro, P. Humire, C. Henkel, P. van der Werf

Abstract: Phosphorus (P) is a crucial element for life given its central role in several biomolecules. P-bearing molecules have been discovered in different regions of the Milky Way, but not yet towards an extragalactic environment. We have searched for P-bearing molecules towards the nearby starburst Galaxy NGC 253. Using observations from the ALMA Comprehensive High-resolution Extragalactic Molecular Inve… ▽ More Phosphorus (P) is a crucial element for life given its central role in several biomolecules. P-bearing molecules have been discovered in different regions of the Milky Way, but not yet towards an extragalactic environment. We have searched for P-bearing molecules towards the nearby starburst Galaxy NGC 253. Using observations from the ALMA Comprehensive High-resolution Extragalactic Molecular Inventory (ALCHEMI) project, we used the MADCUBA package to model the emission of P-bearing molecules assuming Local Thermodynamic Equilibrium (LTE) conditions. We have also performed a non-LTE analysis using SpectralRadex. We report the detection of a P-bearing molecule, phosphorus nitride (PN), for the first time in an extragalactic environment, towards two giant molecular clouds (GMCs) of NGC 253. The LTE analysis yields total PN beam-averaged column densities $N$=(1.20$\pm$0.09)$\times$10$^{13}$ cm$^{-2}$ and $N$=(6.5$\pm$1.6)$\times$10$^{12}$ cm$^{-2}$, which translate into abundances with respect to H$_2$ of $χかい$=(8.0$\pm$1.0)$\times$10$^{-12}$ and $χかい$=(4.4$\pm$1.2)$\times$10$^{-12}$. We derived a low excitation temperature of $T_{\rm ex}$=(4.4$\pm$1.3) K towards the GMC with the brightest PN emission, which indicates that PN is sub-thermally excited. The non-LTE analysis results in column densities consistent with the LTE values. We have also searched for other P-bearing molecules (PO, PH$_{3}$, CP and CCP), and upper limits were derived. The derived PO/PN ratios are $<$1.3 and $<$1.7. The abundance ratio between PN and the shock-tracer SiO derived towards NGC 253 follows the same trend previously found towards Galactic sources. Comparison of the observations with chemical models indicates that the derived molecular abundances of PN in NGC 253 can be explained by shock-driven chemistry followed by cosmic-ray-driven photochemistry. △ Less

Submitted 15 January, 2022; v1 submitted 9 December, 2021; originally announced December 2021.

Comments: 14 pages, 14 figures, 4 tables. Accepted for publication in A&A

Journal ref: A&A 659, A158 (2022)

ALMA Sub-arcsec-resolution 183 GHz H2O and Dense Molecular Line Observations of Nearby Ultraluminous Infrared Galaxies

Authors: Masatoshi Imanishi, Kouichiro Nakanishi, Takuma Izumi, Shunsuke Baba

Abstract: We present the results of ALMA $\sim$2 mm, $\lesssim$1$''$-resolution observations of ten (ultra)luminous infrared galaxies ([U]LIRGs; infrared luminosity $\gtrsim$ 10$^{11.7}$L$_{\odot}$) at $z <$ 0.15, targeting dense ($>$10$^{4}$ cm$^{-3}$) molecular (HCN, HCO$^{+}$, and HNC J=2-1) and 183 GHz H$_{2}$O 3$_{1,3}$-2$_{2,0}$ emission lines. Higher HCN to HCO$^{+}$ J=2-1 flux ratios are observed in… ▽ More We present the results of ALMA $\sim$2 mm, $\lesssim$1$''$-resolution observations of ten (ultra)luminous infrared galaxies ([U]LIRGs; infrared luminosity $\gtrsim$ 10$^{11.7}$L$_{\odot}$) at $z <$ 0.15, targeting dense ($>$10$^{4}$ cm$^{-3}$) molecular (HCN, HCO$^{+}$, and HNC J=2-1) and 183 GHz H$_{2}$O 3$_{1,3}$-2$_{2,0}$ emission lines. Higher HCN to HCO$^{+}$ J=2-1 flux ratios are observed in some, but not all, AGN-important ULIRGs than in starburst-classified sources. We detect 183 GHz H$_{2}$O emission in almost all AGN-important ULIRGs, and elevated H$_{2}$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$_{2}$O emission largely comes from the entire nuclear regions ($\sim$1 kpc), rather than AGN-origin megamaser at the very center ($<<$1 kpc). Nuclear ($\sim$1 kpc) dense molecular gas mass derived from HCO$^{+}$ J=2-1 luminosity is $\gtrsim$a few $\times$ 10$^{8}$M$_{\odot}$, and its depletion time is estimated to be $\gtrsim$10$^{6}$ 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 $\sim$2 mm data of ten (U)LIRGs, two continuum sources are serendipitously detected within $\sim$10$''$, which are likely to be an infrared luminous dusty galaxy at $z >$ 1 and a blazar. △ Less

Submitted 17 November, 2021; originally announced November 2021.

Comments: 50 pages, 16 Figures. ApJ accepted

Quantum-gate decomposer

Authors: Ken M. Nakanishi, Takahiko Satoh, Synge Todo

Abstract: Efficient decompositions of multi-qubit gates are essential in NISQ applications, where the number of gates or the circuit depth is limited. This paper presents efficient decompositions of CCZ and CCCZ gates, typical multi-qubit gates, under several qubit connectivities. We can construct the CCZ gate with only four CZ-depth when the qubit is square-shaped, including one auxiliary qubit. In T-shape… ▽ More Efficient decompositions of multi-qubit gates are essential in NISQ applications, where the number of gates or the circuit depth is limited. This paper presents efficient decompositions of CCZ and CCCZ gates, typical multi-qubit gates, under several qubit connectivities. We can construct the CCZ gate with only four CZ-depth when the qubit is square-shaped, including one auxiliary qubit. In T-shaped qubit connectivity, which has no closed loop, we can decompose the CCCZ gate with 17 CZ gates. While previous studies have shown a CCCZ gate decomposition with 14 CZ gates for the fully connected case, we found only four connections are sufficient for 14 CZ gates' implementation. The search for constraint-sufficient decompositions is aided by an optimization method we devised to bring the parameterized quantum circuit closer to the target quantum gate. We can apply this scheme to decompose any quantum gates, not only CCZ and CCCZ. Such decompositions of multi-qubit gates, together with the newly found CCZ and CCCZ decompositions, shorten the execution time of quantum circuits and improve the accuracy of complex quantum algorithms on near future QPUs. △ Less

Submitted 27 September, 2021; originally announced September 2021.

Comments: 7 pages, 4 figures

ALCHEMI: an ALMA Comprehensive High-resolution Extragalactic Molecular Inventory. Survey presentation and first results from the ACA array

Authors: S. Martín, J. G. Mangum, N. Harada, F. Costagliola, K. Sakamoto, S. Muller, R. Aladro, K. Tanaka, Y. Yoshimura, K. Nakanishi, R. Herrero-Illana, S. Mühle, S. Aalto, E. Behrens, L. Colzi, K. L. Emig, G. A. Fuller, S. García-Burillo, T. R. Greve, C. Henkel, J. Holdship, P. Humire, L. Hunt, T. Izumi, K. Kohno , et al. (10 additional authors not shown)

Abstract: We used the Atacama Large Millimeter/submillimeter Array (ALMA), covering a nearly contiguous 289 GHz frequency range between 84.2 and 373.2 GHz, to image the continuum and spectral line emission at 1.6\arcsec ($\sim 28$ pc) resolution down to a sensitivity of $30-50$ mK. This article describes the ALMA Comprehensive High-resolution Extragalactic Molecular Inventory (ALCHEMI) Large Program. We foc… ▽ More We used the Atacama Large Millimeter/submillimeter Array (ALMA), covering a nearly contiguous 289 GHz frequency range between 84.2 and 373.2 GHz, to image the continuum and spectral line emission at 1.6\arcsec ($\sim 28$ pc) resolution down to a sensitivity of $30-50$ mK. This article describes the ALMA Comprehensive High-resolution Extragalactic Molecular Inventory (ALCHEMI) Large Program. We focus on the analysis of the spectra extracted from the $15''$ ($\sim255$ pc) resolution ALMA Compact Array data. We model the molecular emission assuming local thermodynamic equilibrium with 78 species detected. Additionally, multiple hydrogen and helium recombination lines are identified. Spectral lines contribute 5 to 36\% of the total emission in frequency bins of 50 GHz. We report the first extragalactic detections of C$_2$H$_5$OH, HOCN, HC$_3$HO, and several rare isotopologues. Isotopic ratios of carbon, oxygen, sulfur, nitrogen and silicon were measure with multiple species. Infrared pumped vibrationaly excited HCN, HNC, and HC$_3$N emission, originating in massive star formation locations, is clearly detected at low resolution, while we do not detect it for HCO$^+$. We suggest high temperature conditions in these regions driving a seemingly "carbon-rich" chemistry which may also explain the observed high abundance of organic species close to those in Galactic hot cores. The $L_{vib}/L_{IR}$ ratio is used as a proxy to estimate a $3\%$ contribution from proto super star cluster to the global infrared emission. Measured isotopic ratios with high dipole moment species agree with those within the central kiloparsec of the Galaxy, while those derived from $\rm^{13}C^{18}O$ are a factor of 5 larger, confirming the existence of multiple ISM components within NGC 253 with different degrees of nucleosynthesis enrichment. ALCHEMI provides a template for early Universe galaxies. △ Less

Submitted 17 September, 2021; originally announced September 2021.

Comments: 60 pages, 32 figures, 15 tables, Accepted for publication at A&A

Journal ref: A&A 656, A46 (2021)

Starburst Energy Feedback Seen Through HCO$^+$/HOC$^+$ Emission in NGC 253 from ALCHEMI

Authors: Nanase Harada, Sergio Martín, Jeffrey G. Mangum, Kazushi Sakamoto, Sebastien Muller, Kunihiko Tanaka, Kouichiro Nakanishi, Rubén Herrero-Illana, Yuki Yoshimura, Stefanie Mühle, Rebeca Aladro, Laura Colzi, Víctor M. Rivilla, Susanne Aalto, Erica Behrens, Christian Henkel, Jonathan Holdship, P. K. Humire, David S. Meier, Yuri Nishimura, Paul P. van der Werf, Serena Viti

Abstract: Molecular abundances are sensitive to UV-photon flux and cosmic-ray ionization rate. In starburst environments, the effects of high-energy photons and particles are expected to be stronger. We examine these astrochemical signatures through multiple transitions of HCO$^+$ and its metastable isomer HOC$^+$ in the center of the starburst galaxy NGC 253 using data from the ALMA large program ALCHEMI.… ▽ More Molecular abundances are sensitive to UV-photon flux and cosmic-ray ionization rate. In starburst environments, the effects of high-energy photons and particles are expected to be stronger. We examine these astrochemical signatures through multiple transitions of HCO$^+$ and its metastable isomer HOC$^+$ in the center of the starburst galaxy NGC 253 using data from the ALMA large program ALCHEMI. The distribution of the HOC$^+$(1-0) integrated intensity shows its association with "superbubbles", cavities created either by supernovae or expanding HII regions. The observed HCO$^+$/HOC$^+$ abundance ratios are $\sim 10-150$, and the fractional abundance of HOC$^+$ relative to H$_2$ is $\sim 1.5\times 10^{-11} - 6\times 10^{-10}$, which implies that the HOC$^+$ abundance in the center of NGC 253 is significantly higher than in quiescent spiral-arm dark clouds in the Galaxy and the Galactic center clouds. Comparison with chemical models implies either an interstellar radiation field of $G_0\gtrsim 10^3$ if the maximum visual extinction is $\gtrsim 5$, or a cosmic-ray ionization rate of $ζぜーた\gtrsim 10^{-14}$ s$^{-1}$ (3-4 orders of magnitude higher than that within clouds in the Galactic spiral-arms) to reproduce the observed results. From the difference in formation routes of HOC$^+$, we propose that a low-excitation line of HOC$^+$ traces cosmic-ray dominated regions, while high-excitation lines trace photodissociation regions. Our results suggest that the interstellar medium in the center of NGC 253 is significantly affected by energy input from UV-photons and cosmic rays, sources of energy feedback. △ Less

Submitted 14 September, 2021; originally announced September 2021.

Comments: 57 pages, 32 figures, accepted for publication in ApJ

Global dynamics around 2-solitons for the nonlinear damped Klein-Gordon equations

Authors: Kenjiro Ishizuka, Kenji Nakanishi

Abstract: Global behavior of solutions is studied for the nonlinear Klein-Gordon equation with a focusing power nonlinearity and a damping term in the energy space on the Euclidean space. We give a complete classification of solutions into 5 types of global behavior for all initial data in a small neighborhood of each superposition of two ground states (2-solitons) with the opposite signs and sufficient spa… ▽ More Global behavior of solutions is studied for the nonlinear Klein-Gordon equation with a focusing power nonlinearity and a damping term in the energy space on the Euclidean space. We give a complete classification of solutions into 5 types of global behavior for all initial data in a small neighborhood of each superposition of two ground states (2-solitons) with the opposite signs and sufficient spatial distance. The neighborhood contains, for each sign of the ground state, the manifold with codimension one in the energy space, consisting of solutions that converge to the ground state at time infinity. The two manifolds are joined at their boundary by the manifold with codimension two of solutions that are asymptotic to 2-solitons moving away from each other. The connected union of these three manifolds separates the rest of the neighborhood into the open set of global decaying solutions and that of blow-up. The main ingredient in the proof is a difference estimate on two solutions starting near 2-solitons and asymptotic to 1-solitons. The main difficulty is in controlling the direction of the two unstable modes attached to 2-solitons, while the soliton interactions are not uniformly integrable in time. It is resolved by showing that the non-scalar part of the interaction between the unstable modes is uniformly integrable due to the symmetry of the equation and the 2-solitons. △ Less

Submitted 8 September, 2021; originally announced September 2021.

Comments: 62 pages

MSC Class: 35L71(Primary); 37L25; 35B40; 35B30(Secondary)

ALMA Measurement of 10 kpc-scale Lensing Power Spectra towards the Lensed Quasar MG J0414+0534

Authors: Kaiki Taro Inoue, Takeo Minezaki, Satoki Matsushita, Koichiro Nakanishi

Abstract: The lensing power spectra for gravitational potential, astrometric shift, and convergence perturbations are powerful probes to investigate dark matter structures on small scales. We report the first lower and upper bounds of these lensing power spectra on angular scale $\sim 1$arcsec towards the anomalous quadruply lensed quasar MG$\,$J0414+0534 at a redshift $z=2.639$. To obtain the spectra, we c… ▽ More The lensing power spectra for gravitational potential, astrometric shift, and convergence perturbations are powerful probes to investigate dark matter structures on small scales. We report the first lower and upper bounds of these lensing power spectra on angular scale $\sim 1$arcsec towards the anomalous quadruply lensed quasar MG$\,$J0414+0534 at a redshift $z=2.639$. To obtain the spectra, we conducted observations of MG$\,$J0414+0534 using the Atacama Large Millimeter/submillimeter Array (ALMA) with high angular resolution (0.02 arcsec to 0.05 arcsec). We developed a new partially non-parametric method in which Fourier coefficients of potential perturbation are adjusted to minimize the difference between linear combinations of weighted mean de-lensed images. Using positions of radio jet components, extended dust emission on scales $>1\,$kpc, and mid-infrared flux ratios, the range of measured convergence, astrometric shift, and potential powers at an angular scale of $\sim 1.1$arcsec (corresponding to an angular wave number of $l=1.2\times 10^6$ or $\sim 9\,$kpc in the primary lens plane) within $1\,σしぐま$ are $\varDelta_κかっぱ=0.021-0.028$, $\varDelta_αあるふぁ=7-9\,$mas, and $\varDelta_ψぷさい=1.2-1.6\,$$\textrm{mas}^2$, respectively. Our result is consistent with the predicted abundance of halos in the line of sight and subhalos in cold dark matter models. Our partially non-parametric lens models suggest a presence of a clump in the vicinity of object Y, a possible dusty dwarf galaxy and some small clumps in the vicinity of other lensed quadruple images. Although much fainter than the previous report, we detected weak continuum emission possibly from object Y with a peak flux of $\sim 100\,μみゅー\textrm{Jy}\, \textrm{beam} ^{-1}$ at the $\sim 4\,σしぐま$ level. △ Less

Submitted 26 July, 2023; v1 submitted 2 September, 2021; originally announced September 2021.

Comments: 54 pages, 35 figures, version accepted in ApJ. Analyses of models with a possible dwarf galaxy, object Y and that of chi^2 fit in the visibility plane were added. A new section 7 about CDM consistency was added. The angular scale of the smallest angular bin in Figure 24 was corrected

Journal ref: The Astrophysical Journal (2023) 954, 197 (33pp) (Open Access)

An ACA Survey of [CI] $^3P_1-^3P_0$, CO $J=4-3$, and Dust Continuum in Nearby U/LIRG

Authors: Tomonari Michiyama, Toshiki Saito, Ken-ichi Tadaki, Junko Ueda, Ming-Yang Zhuang, Juan Molina, Bumhyun Lee, Ran Wang, Alberto Bolatto, Daisuke Iono, Kouichiro Nakanishi, Takuma Izumi, Takuji Yamashita, Luis C. Ho

Abstract: We present the results of surveying [CI] $^3P_1-^3P_0$, $^{12}$CO $J=4-3$, and 630 $μみゅー$m dust continuum emission for 36 nearby ultra/luminous infrared galaxies (U/LIRGs) using the Band 8 receiver mounted on the Atacama Compact Array (ACA) of the Atacama Large Millimeter/submillimeter Array. We describe the survey, observations, data reduction, and results; the main results are as follows. (i) We co… ▽ More We present the results of surveying [CI] $^3P_1-^3P_0$, $^{12}$CO $J=4-3$, and 630 $μみゅー$m dust continuum emission for 36 nearby ultra/luminous infrared galaxies (U/LIRGs) using the Band 8 receiver mounted on the Atacama Compact Array (ACA) of the Atacama Large Millimeter/submillimeter Array. We describe the survey, observations, data reduction, and results; the main results are as follows. (i) We confirmed that [CI] $^3P_1-^3P_0$ has a linear relationship with both the $^{12}$CO $J=4-3$and 630 $μみゅー$m continuum. (ii) In NGC 6052 and NGC 7679, $^{12}$CO $J=4-3$ was detected but [CI] $^3P_1-^3P_0$ was not detected with a [CI] $^3P_1-^3P_0$/ $^{12}$CO $J=4-3$ ratio of $\lesssim0.08$. Two possible scenarios of weak [CI] $^3P_1-^3P_0$ emission are C$^0$-poor/CO-rich environments or an environment with an extremely large [CI] $^3P_1-^3P_0$ missing flux. (iii) There is no clear evidence showing that galaxy mergers, AGNs, and dust temperatures control the ratios of [CI] $^3P_1-^3P_0$/ $^{12}$CO $J=4-3$ and $L'_{\rm [CI](1-0)}/L_{\rm 630μみゅーm}$. (iv) We compare our nearby U/LIRGs with high-z galaxies, such as galaxies on the star formation main sequence (MS) at z$\sim1$ and submillimeter galaxies (SMGs) at $z=2-4$. We found that the mean value for the [CII] $^3P_1$--$^3P_0$/ $^{12}$CO $J=4-3$ ratio of U/LIRGs is similar to that of SMGs but smaller than that of galaxies on the MS. △ Less

Submitted 26 July, 2021; originally announced July 2021.

Comments: 11 pages, 7 figures, and 42 pages for the appendix, accepted for publication in ApJS

The Distribution and Origin of C$_2$H in NGC 253 from ALCHEMI

Authors: J. Holdship, S. Viti, S. Martín, N. Harada, J. Mangum, K. Sakamoto, S. Muller, K. Tanaka, Y. Yoshimura, K. Nakanishi, R. Herrero-Illana, S. Mühle, R. Aladro, L. Colzi, K. L. Emig, S. García-Burillo, C. Henkel, P. Humire, D. S. Meier, V. M. Rivilla, P. van der Werf

Abstract: Observations of chemical species can provide an insight into the physical conditions of the emitting gas but it is important to understand how their abundances and excitation vary within different heating environments. C$_2$H is a molecule typically found in PDR regions of our own Galaxy but there is evidence to suggest it also traces other regions undergoing energetic processing in extragalactic… ▽ More Observations of chemical species can provide an insight into the physical conditions of the emitting gas but it is important to understand how their abundances and excitation vary within different heating environments. C$_2$H is a molecule typically found in PDR regions of our own Galaxy but there is evidence to suggest it also traces other regions undergoing energetic processing in extragalactic environments. As part of the ALCHEMI ALMA large program, the emission of C$_2$H in the central molecular zone of the nearby starburst galaxy NGC 253 was mapped at 1.6 " (28 pc) resolution and characterized to understand its chemical origins. Spectral modelling of the N=1-0 through N=4-3 rotational transitions of C$_2$H was used to derive the C$_2$H column densities towards the dense clouds in NGC 253. Chemical modelling, including PDR, dense cloud, and shock models were then used to investigate the chemical processes and physical conditions that are producing the molecular emission. We find high C$_2$H column densities of $\sim 10^{15} cm^{-3}$ detected towards the dense regions of NGC 253. We further find that these column densities cannot be reproduced by assuming that the emission arises from the PDR regions at the edge of the clouds. Instead, we find that the C$_2$H abundance remains high even in the high visual extinction interior of these clouds and that this is most likely caused by a high cosmic-ray ionization rate. △ Less

Submitted 16 July, 2021; v1 submitted 9 July, 2021; originally announced July 2021.

Comments: 18 pages, 13 figures, accepted for publication in A&A

Journal ref: A&A 654, A55 (2021)

ALMA Observations of Lyman-alpha Blob 1: Multiple major-mergers and widely distributed interstellar media

Authors: Hideki Umehata, Ian Smail, Charles C. Steidel, Matthew Hayes, Douglas Scott, A. M. Swinbank, R. J. Ivison, Toru Nagao, Mariko Kubo, Kouichiro Nakanishi, Yuichi Matsuda, Soh Ikarashi, Yoichi Tamura, J. E. Geach

Abstract: We present observations of a giant Lyman-alpha blob in the SSA22 proto-cluster at z=3.1, SSA22-LAB1, taken with the Atacama Large Millimeter/submillimeter Array (ALMA). Dust continuum, along with [C II]158um, and CO(4-3) line emission have been detected in LAB1, showing complex morphology and kinematics across a ~100 kpc central region. Seven galaxies at z=3.0987-3.1016 in the surroundings are ide… ▽ More We present observations of a giant Lyman-alpha blob in the SSA22 proto-cluster at z=3.1, SSA22-LAB1, taken with the Atacama Large Millimeter/submillimeter Array (ALMA). Dust continuum, along with [C II]158um, and CO(4-3) line emission have been detected in LAB1, showing complex morphology and kinematics across a ~100 kpc central region. Seven galaxies at z=3.0987-3.1016 in the surroundings are identified in [C II] and dust continuum emission, with two of them potential companions or tidal structures associated with the most massive galaxies. Spatially resolved [C II] and infrared luminosity ratios for the widely distributed media (L[C II]/LIR~0.01-0.001) suggest that the observed extended interstellar media are likely to have originated from star-formation activity and the contribution from shocked gas is probably not dominant. LAB1 is found to harbour a total molecular gas mass Mmol=(8.7+/-2.0)e+10 Msun, concentrated in the core region of the Ly-alpha-emitting area. While (primarily obscured) star-formation activity in the LAB1 core is one of the most plausible power sources for the Ly-alpha emission, multiple major-mergers found in the core may also play a role in making LAB1 exceptionally bright and extended in Ly-alpha as a result of cooling radiation induced by gravitational interactions. △ Less

Submitted 2 July, 2021; originally announced July 2021.

Comments: 23 pages, 11 figures, accepted for publication in ApJ

A massive quiescent galaxy confirmed in a protocluster at z=3.09

Authors: Mariko Kubo, Hideki Umehata, Yuichi Matsuda, Masaru Kajisawa, Charles C. Steidel, Toru Yamada, Ichi Tanaka, Bunyo Hatsukade, Yoichi Tamura, Kouichiro Nakanishi, Kotaro Kohno, Chien-Feng Lee, Keiichi Matsuda

Abstract: We report a massive quiescent galaxy at $z_{\rm spec}=3.0922^{+0.008}_{-0.004}$ spectroscopically confirmed at a protocluster in the SSA22 field by detecting the Balmer and Ca {\footnotesize II} absorption features with multi-object spectrometer for infrared exploration (MOSFIRE) on the Keck I telescope. This is the most distant quiescent galaxy confirmed in a protocluster to date. We fit the opti… ▽ More We report a massive quiescent galaxy at $z_{\rm spec}=3.0922^{+0.008}_{-0.004}$ spectroscopically confirmed at a protocluster in the SSA22 field by detecting the Balmer and Ca {\footnotesize II} absorption features with multi-object spectrometer for infrared exploration (MOSFIRE) on the Keck I telescope. This is the most distant quiescent galaxy confirmed in a protocluster to date. We fit the optical to mid-infrared photometry and spectrum simultaneously with spectral energy distribution (SED) models of parametric and nonparametric star formation histories (SFH). Both models fit the observed SED well and confirm that this object is a massive quiescent galaxy with the stellar mass of $\log(\rm M_{\star}/M_{\odot}) = 11.26^{+0.03}_{-0.04}$ and $11.54^{+0.03}_{-0.00}$, and star formation rate of $\rm SFR/M_{\odot}~yr^{-1} <0.3$ and $=0.01^{+0.03}_{-0.01}$ for parametric and nonparametric models, respectively. The SFH from the former modeling is described as an instantaneous starburst while that of the latter modeling is longer-lived but both models agree with a sudden quenching of the star formation at $\sim0.6$ Gyr ago. This massive quiescent galaxy is confirmed in an extremely dense group of galaxies predicted as a progenitor of a brightest cluster galaxy formed via multiple mergers in cosmological numerical simulations. We newly find three plausible [O III]$λらむだ$5007 emitters at $3.0791\leq z_{\rm spec}\leq3.0833$ happened to be detected around the target. Two of them just between the target and its nearest massive galaxy are possible evidence of their interactions. They suggest the future strong size and stellar mass evolution of this massive quiescent galaxy via mergers. △ Less

Submitted 20 June, 2021; originally announced June 2021.

Comments: 13 pages, 7 figures (in main text. 21 pages including appendices), accepted for publication in ApJ

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). XIII. Large-scale Feedback and Star Formation in a Low-Luminosity Quasar at z = 7.07 on the Local Black Hole to Host Mass Relation

Authors: Takuma Izumi, Yoshiki Matsuoka, Seiji Fujimoto, Masafusa Onoue, Michael A. Strauss, Hideki Umehata, Masatoshi Imanishi, Kotaro Kohno, Toshihiro Kawaguchi, Taiki Kawamuro, Shunsuke Baba, Tohru Nagao, Yoshiki Toba, Kohei Inayoshi, John D. Silverman, Akio K. Inoue, Soh Ikarashi, Kazushi Iwasawa, Nobunari Kashikawa, Takuya Hashimoto, Kouichiro Nakanishi, Yoshihiro Ueda, Malte Schramm, Chien-Hsiu Lee, Hyewon Suh

Abstract: We present ALMA [CII] 158 $μみゅー$m line and underlying far-infrared (FIR) continuum emission observations ($0''.70 \times 0''.56$ resolution) toward HSC J124353.93$+$010038.5 (J1243$+$0100) at $z = 7.07$, the only low-luminosity ($M_{\rm 1450} > -25$ mag) quasar currently known at $z > 7$. The FIR continuum is bright (1.52 mJy) and resolved with a total luminosity of… ▽ More We present ALMA [CII] 158 $μみゅー$m line and underlying far-infrared (FIR) continuum emission observations ($0''.70 \times 0''.56$ resolution) toward HSC J124353.93$+$010038.5 (J1243$+$0100) at $z = 7.07$, the only low-luminosity ($M_{\rm 1450} > -25$ mag) quasar currently known at $z > 7$. The FIR continuum is bright (1.52 mJy) and resolved with a total luminosity of $L_{\rm FIR} = 3.5 \times 10^{12}~L_\odot$. The spatially extended component is responsible for $\sim 40\%$ of the emission. The area-integrated [CII] spectrum shows a broad wing (${\rm FWHM} = 997$ km s$^{-1}$, $L_{\rm [CII]} = 1.2 \times 10^9~L_\odot$) as well as a bright core (${\rm FWHM} = 235$ km s$^{-1}$, $L_{\rm [CII]} = 1.9 \times 10^9~L_\odot$). This wing is the first detection of a galactic-scale quasar-driven outflow (atomic outflow rate $> 447~M_\odot$ yr$^{-1}$) at $z > 7$. The estimated large mass loading factor of the total outflow (e.g., $\gtrsim 9$ relative to the [CII]-based SFR) suggests that this outflow will soon quench the star-formation of the host. The core gas dynamics are governed by rotation, with a rotation curve suggestive of a compact bulge ($\sim 3.3 \times 10^{10}~M_\odot$), although it is not yet spatially resolved. Finally, we found that J1243$+$0100 has a black hole mass-to-dynamical mass ratio (and -to-bulge mass ratio) of $\sim 0.4\%$ ($\sim 1\%$), consistent with the local value within uncertainties. Our results therefore suggest that the black hole-host co-evolution relation is already in place at $z \sim 7$ for this object. △ Less

Submitted 20 May, 2021; v1 submitted 12 April, 2021; originally announced April 2021.

Comments: Slightly revised version after the proof correction. 22 pages, 13 figures, 6 tables. Accepted for publication in ApJ

A phase-space view of cold-gas properties of Virgo-cluster galaxies: multiple quenching processes at work?

Authors: Kana Morokuma-Matsui, Tadayuki Kodama, Tomoki Morokuma, Kouichiro Nakanishi, Yusei Koyama, Takuji Yamashita, Shuhei Koyama, Takashi Okamoto

Abstract: We investigate the cold-gas properties of massive Virgo galaxies ($>10^9$ M$_\odot$) at $<3R_{200}$ ($R_{200}$ is the radius where the mean interior density is 200 times the critical density) on the projected phase-space diagram (PSD) with the largest archival dataset to date to understand the environmental effect on galaxy evolution in the Virgo cluster. We find: lower HI and H$_2$ mass fractions… ▽ More We investigate the cold-gas properties of massive Virgo galaxies ($>10^9$ M$_\odot$) at $<3R_{200}$ ($R_{200}$ is the radius where the mean interior density is 200 times the critical density) on the projected phase-space diagram (PSD) with the largest archival dataset to date to understand the environmental effect on galaxy evolution in the Virgo cluster. We find: lower HI and H$_2$ mass fractions and higher star-formation efficiencies (SFEs) from HI and H$_2$ in the Virgo galaxies than the field galaxies for matched stellar masses; the Virgo galaxies generally follow the field relationships between the offset from the main sequence of the star-forming galaxies [$Δでるた$(MS)] with gas fractions and SFEs but slightly offset to lower gas fractions or higher SFEs than field galaxies at $Δでるた({\rm MS})< 0$; lower gas fractions in galaxies with smaller clustocentric distance and velocity; lower gas fractions in the galaxies in the W cloud, a substructure of the Virgo cluster. Our results suggest the cold-gas properties of some Virgo galaxies are affected by their environment at least at $3 R_{200}$ maybe via strangulation and/or pre-processes and HI and H$_2$ in some galaxies are removed by ram pressure at $<1.5 R_{200}$. Our data cannot rule the possibility of the other processes such as strangulation and galaxy harassment accounting for the gas reduction in some galaxies at $<1.5 R_{200}$. Future dedicated observations of a mass-limited complete sample are required for definitive conclusions. △ Less

Submitted 24 March, 2021; v1 submitted 9 March, 2021; originally announced March 2021.

Comments: 32 pages, 16 figures, accepted for publication in ApJ; Updated stellar masses of galaxies after the first acceptance, relevant texts, and figures

Dense and warm neutral gas in BR1202-0725 at z = 4.7 as traced by the [O I] 145 um line

Authors: Minju M. Lee, Tohru Nagao, Carlos De Breuck, Stefano Carniani, Giovanni Cresci, Bunyo Hatsukade, Ryohei Kawabe, Kotaro Kohno, Roberto Maiolino, Filippo Mannucci, Alessandro Marconi, Kouichiro Nakanishi, Paulina Troncoso, Hideki Umehata

Abstract: We report the detection of [O I]145.5um in the BR 1202-0725 system, a compact group at z=4.7 consisting of a quasar (QSO), a submillimeter-bright galaxy (SMG), and three faint Lya emitters. By taking into account the previous detections and upper limits, the [O I]/[C II] line ratios of the now five known high-z galaxies are higher than or on the high-end of the observed values in local galaxies ([… ▽ More We report the detection of [O I]145.5um in the BR 1202-0725 system, a compact group at z=4.7 consisting of a quasar (QSO), a submillimeter-bright galaxy (SMG), and three faint Lya emitters. By taking into account the previous detections and upper limits, the [O I]/[C II] line ratios of the now five known high-z galaxies are higher than or on the high-end of the observed values in local galaxies ([O I]/[C II]$\gtrsim$0.13). The high [O I]/[C II] ratios and the joint analysis with the previous detection of [N II] lines for both the QSO and the SMG suggest the presence of warm and dense neutral gas in these highly star-forming galaxies. This is further supported by new CO (12-11) line detections and a comparison with cosmological simulations. There is a possible positive correlation between the [NII]122/205 line ratio and the [O I]/[C II] ratio when all local and high-z sources are taken into account, indicating that the denser the ionized gas, the denser and warmer the neutral gas (or vice versa). The detection of the [O I] line in the BR1202-0725 system with a relatively short amount of integration with ALMA demonstrates the great potential of this line as a dense gas tracer for high-z galaxies. △ Less

Submitted 18 February, 2021; originally announced February 2021.

Comments: 11 pages, 7 figures, 1 table, Accepted for publication in ApJ, figures 2 and 4 show the main results

ALMA detection of millimetre 183 GHz H2O maser emission in the Superantennae galaxy at z ~ 0.06

Authors: Masatoshi Imanishi, Yoshiaki Hagiwara, Shinji Horiuchi, Takuma Izumi, Kouichiro Nakanishi

Abstract: We present the results of ALMA band-5 (~170 GHz) observations of the merging ultraluminous infrared galaxy, the "Superantennae" (IRAS 19254-7245) at z=0.0617, which has been diagnosed as containing a luminous obscured active galactic nucleus (AGN). In addition to dense molecular line emission (HCN, HCO+, and HNC J = 2-1), we detect a highly luminous (~6e4Lsun) 183 GHz H2O 3(1,3)-2(2,0) emission li… ▽ More We present the results of ALMA band-5 (~170 GHz) observations of the merging ultraluminous infrared galaxy, the "Superantennae" (IRAS 19254-7245) at z=0.0617, which has been diagnosed as containing a luminous obscured active galactic nucleus (AGN). In addition to dense molecular line emission (HCN, HCO+, and HNC J = 2-1), we detect a highly luminous (~6e4Lsun) 183 GHz H2O 3(1,3)-2(2,0) emission line. We interpret the strong H2O emission as largely originating in maser amplification in AGN-illuminated dense and warm molecular gas, based on (1) the spatially compact (<220 pc) nature of the H2O emission, unlike spatially resolved (>500 pc) dense molecular emission, and (2) a strikingly different velocity profile from, and (3) significantly elevated flux ratio relative to, dense molecular emission lines. H2O maser emission, other than the widely studied 22 GHz 6(1,6)-5(2,3) line, has been expected to provide important information on the physical properties of gas in the vicinity of a central mass-accreting supermassive black hole (SMBH), because of different excitation energy. We here demonstrate that with highly sensitive ALMA, millimetre 183 GHz H2O maser detection is feasible out to >270 Mpc, opening a new window to scrutinize molecular gas properties around a mass-accreting SMBH far beyond the immediately local universe. △ Less

Submitted 13 January, 2021; originally announced January 2021.

Comments: 6 pages, 2 figures. MNRAS Letters accepted

Improving ALMA's data processing effciency using a holistic approach

Authors: Theodoros Nakos, Harold Francke, Kouichiro Nakanishi, Dirk Petry, Thomas Stanke, Catarina Ubach, Luciano Cerrigone, Erica Keller, Alfonso Trejo, Junko Ueda

Abstract: ALMA (Atacama Large Millimeter/submillimeter Array) is the world's largest ground-based facility for observations in the millimeter/submillimeter regime. One of ALMA's outstanding characteristics is the large effort dedicated to the quality assurance (QA) of the calibrated and imaged data products offered to the astronomical community. The Data Management Group (DMG), in charge of the data process… ▽ More ALMA (Atacama Large Millimeter/submillimeter Array) is the world's largest ground-based facility for observations in the millimeter/submillimeter regime. One of ALMA's outstanding characteristics is the large effort dedicated to the quality assurance (QA) of the calibrated and imaged data products offered to the astronomical community. The Data Management Group (DMG), in charge of the data processing, review, and delivery of the ALMA data, consists of approximately 60 experts in data reduction, from the ALMA Regional Centers (ARCs) and the Joint ALMA Observatory (JAO), distributed in fourteen countries. With a throughput of more than 3,000 datasets per year, meeting the goal of delivering the pipeline-able data products within 30 days after data acquisition is a huge challenge. This paper presents (a) the history of data processing at ALMA, (b) the challenges our team had and is still facing, (c) the methodology followed to mitigate the operational risks, (d) the ongoing optimization initiatives, (e) the current data processing status, (f) the strategy which is being followed so that, in a few Cycles from now, a team of approximately 10 data reducers (DRs) at JAO can process and review some 80% of the datasets collected during an observing cycle, and, finally, (g) the important role of the ARCs for processing the remaining datasets. △ Less

Submitted 9 January, 2021; originally announced January 2021.

Comments: 18 pages, 5 figures, paper presented at the December 2020 SPIE virtual conference on Astronomical Telescopes and Instrumentation

Journal ref: Proc. SPIE 11449, Observatory Operations: Strategies, Processes, and Systems VIII, 114491T (15 December 2020)

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) XII. Extended [C II] Structure (Merger or Outflow) in a z = 6.72 Red Quasar

Authors: Takuma Izumi, Masafusa Onoue, Yoshiki Matsuoka, Michael A. Strauss, Seiji Fujimoto, Hideki Umehata, Masatoshi Imanishi, Taiki Kawamuro, Tohru Nagao, Yoshiki Toba, Kotaro Kohno, Nobunari Kashikawa, Kohei Inayoshi, Toshihiro Kawaguchi, Kazushi Iwasawa, Akio K. Inoue, Tomotsugu Goto, Shunsuke Baba, Malte Schram, Hyewon Suh, Yuichi Harikane, Yoshihiro Ueda, John D. Silverman, Takuya Hashimoto, Yasuhiro Hashimoto , et al. (9 additional authors not shown)

Abstract: We present ALMA [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 the beam size of $\sim 0''.8 \times 0''.5$ (or 4.1 kpc $\times$ 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 unobsc… ▽ More We present ALMA [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 the beam size of $\sim 0''.8 \times 0''.5$ (or 4.1 kpc $\times$ 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 which blows out a cocoon of interstellar medium (ISM). The FIR continuum of J1205$-$0000 is bright, with an estimated luminosity of $L_{\rm FIR} \sim 3 \times 10^{12}~L_\odot$. The [C II] line emission is extended on scales of $r \sim 5$ kpc, greater than the FIR continuum. The line profiles at the extended regions are complex and broad (FWHM $\sim 630-780$ km s$^{-1}$). 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 ($\sim 10~M_\odot$ yr$^{-1}$), but are not detected by our Subaru optical observations ($y_{\rm AB,5σしぐま} = 24.4$ mag). For the case of (ii), our lower-limit to the cold neutral outflow rate is $\sim 100~M_\odot$ yr$^{-1}$. The outflow kinetic energy and momentum are both much smaller than what predicted in energy-conserving wind models, suggesting that the AGN feedback in this quasar is not capable of completely suppressing its star formation. △ Less

Submitted 4 January, 2021; originally announced January 2021.

Comments: 12 pages, 4 figures, 2 tables, accepted for publication in ApJ

Qulacs: a fast and versatile quantum circuit simulator for research purpose

Authors: Yasunari Suzuki, Yoshiaki Kawase, Yuya Masumura, Yuria Hiraga, Masahiro Nakadai, Jiabao Chen, Ken M. Nakanishi, Kosuke Mitarai, Ryosuke Imai, Shiro Tamiya, Takahiro Yamamoto, Tennin Yan, Toru Kawakubo, Yuya O. Nakagawa, Yohei Ibe, Youyuan Zhang, Hirotsugu Yamashita, Hikaru Yoshimura, Akihiro Hayashi, Keisuke Fujii

Abstract: To explore the possibilities of a near-term intermediate-scale quantum algorithm and long-term fault-tolerant quantum computing, a fast and versatile quantum circuit simulator is needed. Here, we introduce Qulacs, a fast simulator for quantum circuits intended for research purpose. We show the main concepts of Qulacs, explain how to use its features via examples, describe numerical techniques to s… ▽ More To explore the possibilities of a near-term intermediate-scale quantum algorithm and long-term fault-tolerant quantum computing, a fast and versatile quantum circuit simulator is needed. Here, we introduce Qulacs, a fast simulator for quantum circuits intended for research purpose. We show the main concepts of Qulacs, explain how to use its features via examples, describe numerical techniques to speed-up simulation, and demonstrate its performance with numerical benchmarks. △ Less

Submitted 5 October, 2021; v1 submitted 26 November, 2020; originally announced November 2020.

Comments: 34 pages, 12 figures

Journal ref: Quantum 5, 559 (2021)

The Explosion in Orion-KL as Seen by Mosaicking the Magnetic Field with ALMA

Authors: Paulo C. Cortes, Valentin J. M. Le Gouellec, Charles L. H. Hull, Josep M. Girart, Fabien Louvet, Edward B. Fomalont, Seiji Kameno, George A. Moellenbrock, Hiroshi Nagai, Kouichiro Nakanishi, Eric Villard

Abstract: We present the first linear-polarization mosaicked observations performed by the Atacama Large Millimeter/submillimeter Array (ALMA). We mapped the Orion-KLeinmann-Low (Orion-KL) nebula using super-sampled mosaics at 3.1 and 1.3 mm as part of the ALMA Extension and Optimization of Capabilities (EOC) program. We derive the magnetic field morphology in the plane of the sky by assuming that dust grai… ▽ More We present the first linear-polarization mosaicked observations performed by the Atacama Large Millimeter/submillimeter Array (ALMA). We mapped the Orion-KLeinmann-Low (Orion-KL) nebula using super-sampled mosaics at 3.1 and 1.3 mm as part of the ALMA Extension and Optimization of Capabilities (EOC) program. We derive the magnetic field morphology in the plane of the sky by assuming that dust grains are aligned with respect to the ambient magnetic field. At the center of the nebula, we find a quasi-radial magnetic field pattern that is aligned with the explosive CO outflow up to a radius of approximately 12 arc-seconds (~ 5000 au), beyond which the pattern smoothly transitions into a quasi-hourglass shape resembling the morphology seen in larger-scale observations by the James-Clerk-Maxwell Telescope (JCMT). We estimate an average magnetic field strength $\langle B\rangle = 9.4$ mG and a total magnetic energy of 2 x 10^45 ergs, which is three orders of magnitude less than the energy in the explosive CO outflow. We conclude that the field has been overwhelmed by the outflow and that a shock is propagating from the center of the nebula, where the shock front is seen in the magnetic field lines at a distance of ~ 5000 au from the explosion center. △ Less

Submitted 15 November, 2020; originally announced November 2020.

Comments: Accepted for publication in ApJ

Post-Hartree-Fock method in Quantum Chemistry for Quantum Computer

Authors: Yutaka Shikano, Hiroshi C. Watanabe, Ken M. Nakanishi, Yu-ya Ohnishi

Abstract: Quantum computational chemistry is a potential application of quantum computers that is expected to effectively solve several quantum-chemistry problems, particularly the electronic structure problem. Quantum computational chemistry can be compared to the conventional computational devices. This review comprehensively investigates the applications and overview of quantum computational chemistry, i… ▽ More Quantum computational chemistry is a potential application of quantum computers that is expected to effectively solve several quantum-chemistry problems, particularly the electronic structure problem. Quantum computational chemistry can be compared to the conventional computational devices. This review comprehensively investigates the applications and overview of quantum computational chemistry, including a review of the Hartree-Fock method for quantum information scientists. Quantum algorithms, quantum phase estimation, and variational quantum eigensolver, have been applied to the post-Hartree-Fock method. △ Less

Submitted 3 November, 2020; originally announced November 2020.

Comments: 31 pages, 6 figures

Journal ref: Eur. Phys. J. Spec. Top. 230, 1037-1051 (2021)