Showing 1–50 of 153 results for author: Takeuchi, T T

Density-based clustering algorithm for galaxy group/cluster identification

Authors: Hai-Xia Ma, Tsutomu T. Takeuchi, Suchetha Cooray, Yongda Zhu

Abstract: A direct approach to studying the galaxy-halo connection is the analysis of observed groups and clusters of galaxies that trace the underlying dark matter halos, making identifying galaxy clusters and their associated brightest cluster galaxies (BCGs) crucial. We test and propose a robust density-based clustering algorithm that outperforms the traditional Friends-of-Friends (FoF) algorithm in the… ▽ More A direct approach to studying the galaxy-halo connection is the analysis of observed groups and clusters of galaxies that trace the underlying dark matter halos, making identifying galaxy clusters and their associated brightest cluster galaxies (BCGs) crucial. We test and propose a robust density-based clustering algorithm that outperforms the traditional Friends-of-Friends (FoF) algorithm in the currently available galaxy group/cluster catalogs. Our new approach is a modified version of the Ordering Points To Identify the Clustering Structure (OPTICS) algorithm, which accounts for line-of-sight positional uncertainties due to redshift space distortions by incorporating a scaling factor, and is thereby referred to as sOPTICS. When tested on both a galaxy group catalog based on semi-analytic galaxy formation simulations and observational data, our algorithm demonstrated robustness to outliers and relative insensitivity to hyperparameter choices. In total, we compared the results of eight clustering algorithms. The proposed density-based clustering method, sOPTICS, outperforms FoF in accurately identifying giant galaxy clusters and their associated BCGs in various environments with higher purity and recovery rate, also successfully recovering 115 BCGs out of 118 reliable BCGs from a large galaxy sample. △ Less

Submitted 16 May, 2024; originally announced May 2024.

Near-Ultraviolet Radiation toward Molecular Cloud N4 in W 50/SS 433: Evidence for Direct Interaction of the Jet with Molecular Cloud

Authors: H. Yamamoto, T. Ishikawa, T. T. Takeuchi

Abstract: We compared the molecular clouds in the western part of SS 433 with near-ultraviolet radiation data obtained from GALEX. Near-ultraviolet radiation is prominently confirmed toward only N4, while no near-ultraviolet radiation is detected toward N1, N2, and N3. The radiative region of near-ultraviolet radiation is nearly the same as the CO-emitting region in N4, and does not extend beyond the jet se… ▽ More We compared the molecular clouds in the western part of SS 433 with near-ultraviolet radiation data obtained from GALEX. Near-ultraviolet radiation is prominently confirmed toward only N4, while no near-ultraviolet radiation is detected toward N1, N2, and N3. The radiative region of near-ultraviolet radiation is nearly the same as the CO-emitting region in N4, and does not extend beyond the jet seen in X-ray radiation. Near-ultraviolet radiation cannot be explained solely by broadband continuous radiation and may originate from line emissions. The intensity of near-ultraviolet radiation exhibits an anti-correlation with that of 13CO(J=3-2) emission. This anti-correlation, along with strong far-infrared radiation in the region with weaker near-ultraviolet radiation intensity compared to its surroundings, suggests that near-ultraviolet radiation originates from behind the molecular cloud, heating up the interstellar dust in N4. Subsequently, the dust in N4 re-radiates in the far-infrared band. In the same region, a high peak T_MB ratio of 12CO(J=3-2)/12CO(J=1-0) of ~0.9, and a high kinetic temperature of T_k ~56 K in the molecular cloud indicate that CO molecules are highly excited, and the molecular cloud is heated through photoelectric heating. This heating results from electrons released due to the photoelectric effect caused by the phenomenon where interstellar dust absorbs near-ultraviolet radiation. In terms of the timescale of near-ultraviolet radiation originating from line emissions, near-ultraviolet radiation towards N4 cannot be explained by the shock of the blast wave from a supernova that created W 50. These findings also suggest that N4 directly interacts with the jet from SS 433. As a result of this direct interaction, near-ultraviolet radiation is emitted from an interacting layer between the jet and N4. △ Less

Submitted 12 January, 2024; originally announced January 2024.

Comments: 16 pages, 7 figures, Accepted for publication in PASJ

The fate of the interstellar medium in early-type galaxies. III. The mechanism of ISM removal and quenching of star formation

Authors: Michał J. Michałowski, C. Gall, J. Hjorth, D. T. Frayer, A. -L. Tsai, K. Rowlands, T. T. Takeuchi, A. Leśniewska, D. Behrendt, N. Bourne, D. H. Hughes, M. P. Koprowski, J. Nadolny, O. Ryzhov, M. Solar, E. Spring, J. Zavala, P. Bartczak

Abstract: Understanding how galaxies quench their star formation is crucial for studies of galaxy evolution. Quenching is related to the cold gas decrease. In the first paper we showed that the dust removal timescale in early-type galaxies (ETGs) is about 2.5 Gyr. Here we present carbon monoxide (CO) and 21 cm hydrogen (H I) line observations of these galaxies and measure the timescale of removal of the col… ▽ More Understanding how galaxies quench their star formation is crucial for studies of galaxy evolution. Quenching is related to the cold gas decrease. In the first paper we showed that the dust removal timescale in early-type galaxies (ETGs) is about 2.5 Gyr. Here we present carbon monoxide (CO) and 21 cm hydrogen (H I) line observations of these galaxies and measure the timescale of removal of the cold interstellar medium (ISM). We find that all the cold ISM components (dust, molecular and atomic gas) decline at similar rates. This allows us to rule out a wide range of potential ISM removal mechanisms (including starburst-driven outflows, astration, a decline in the number of asymptotic giant branch stars), and artificial effects like stellar mass-age correlation, environmental influence, mergers, and selection bias, leaving ionization by evolved low-mass stars and ionization/outflows by supernovae Type Ia or active galactic nuclei as viable mechanisms. We also provide evidence for an internal origin of the detected ISM. Moreover, we find that the quenching of star formation in these galaxies cannot be explained by a reduction in gas amount alone, because the star formation rates (SFRs) decrease faster (on a timescale of about 1.8 Gyr) than the amount of cold gas. Furthermore, the star formation efficiency of the ETGs (SFE = SFR/MH2) is lower than that of star-forming galaxies, whereas their gas mass fractions (fH2 = MH2/M*) are normal. This may be explained by the stabilization of gas against fragmentation, for example due to morphological quenching, turbulence, or magnetic fields. △ Less

Submitted 26 March, 2024; v1 submitted 9 January, 2024; originally announced January 2024.

Comments: Astrophysical Journal, published; 29 pages, 13 figures, 6 tables; V2: matched the published version and typos in the SFR and outflow equations fixed - no change in the conclusions

Journal ref: 2024, ApJ, 964, 129

Extracting an Informative Latent Representation of High-Dimensional Galaxy Spectra

Authors: Daiki Iwasaki, Suchetha Cooray, Tsutomu T. Takeuchi

Abstract: To understand the fundamental parameters of galaxy evolution, we investigated the minimum set of parameters that explain the observed galaxy spectra in the local Universe. We identified four latent variables that efficiently represent the diversity of high-dimensional galaxy spectral energy distributions (SEDs) observed by the Sloan Digital Sky Survey. Additionally, we constructed meaningful laten… ▽ More To understand the fundamental parameters of galaxy evolution, we investigated the minimum set of parameters that explain the observed galaxy spectra in the local Universe. We identified four latent variables that efficiently represent the diversity of high-dimensional galaxy spectral energy distributions (SEDs) observed by the Sloan Digital Sky Survey. Additionally, we constructed meaningful latent representation using conditional variational autoencoders trained with different permutations of galaxy physical properties, which helped us quantify the information that these traditionally used properties have on the reconstruction of galaxy spectra. The four parameters suggest a view that complex SED population models with a very large number of parameters will be difficult to constrain even with spectroscopic galaxy data. Through an Explainable AI (XAI) method, we found that the region below 5000\textupÅ and prominent emission lines ([O II], [O III], and H$α$) are particularly informative for predicting the latent variables. Our findings suggest that these latent variables provide a more efficient and fundamental representation of galaxy spectra than conventionally considered galaxy physical properties. △ Less

Submitted 29 November, 2023; originally announced November 2023.

Comments: 5 pages, 6 figures, accepted by NeurIPS 2023

Galaxy mergers in Subaru HSC-SSP: a deep representation learning approach for identification and the role of environment on merger incidence

Authors: Kiyoaki Christopher Omori, Connor Bottrell, Mike Walmsley, Hassen M. Yesuf, Andy D. Goulding, Xuheng Ding, Gergö Popping, John D. Silverman, Tsutomu T. Takeuchi, Yoshiki Toba

Abstract: We take a deep learning-based approach for galaxy merger identification in Subaru HSC-SSP, specifically through the use of deep representation learning and fine-tuning, with the aim of creating a pure and complete merger sample within the HSC-SSP survey. We can use this merger sample to conduct studies on how mergers affect galaxy evolution. We use Zoobot, a deep learning representation learning m… ▽ More We take a deep learning-based approach for galaxy merger identification in Subaru HSC-SSP, specifically through the use of deep representation learning and fine-tuning, with the aim of creating a pure and complete merger sample within the HSC-SSP survey. We can use this merger sample to conduct studies on how mergers affect galaxy evolution. We use Zoobot, a deep learning representation learning model pre-trained on citizen science votes on Galaxy Zoo DeCALS images. We fine-tune Zoobot for the purpose of merger classification of images of SDSS and GAMA galaxies in HSC-SSP PDR 3. Fine-tuning is done using 1200 synthetic HSC-SSP images of galaxies from the TNG simulation. We then find merger probabilities on observed HSC images using the fine-tuned model. Using our merger probabilities, we examine the relationship between merger activity and environment. We find that our fine-tuned model returns an accuracy on the synthetic validation data of 76%. This number is comparable to those of previous studies where convolutional neural networks were trained with simulation images, but with our work requiring a far smaller number of training samples. For our synthetic data, our model is able to achieve completeness and precision values of 80%. In addition, our model is able to correctly classify both mergers and non-mergers of diverse morphologies and structures, including those at various stages and mass ratios, while distinguishing between projections and merger pairs. For the relation between galaxy mergers and environment, we find two distinct trends. Using stellar mass overdensity estimates for TNG simulations and observations using SDSS and GAMA, we find that galaxies with higher merger scores favor lower density environments on scales of 0.5 to 8 h^-1 Mpc. However, below these scales in the simulations, we find that galaxies with higher merger scores favor higher density environments. △ Less

Submitted 27 September, 2023; originally announced September 2023.

Comments: 36 pages, 15 figures, accepted to Astronomy and Astrophysics

Probing Ultra-late Reionization: Direct Measurements of the Mean Free Path over $5<z<6$

Authors: Yongda Zhu, George D. Becker, Holly M. Christenson, Anson D'Aloisio, Sarah E. I. Bosman, Tom Bakx, Valentina D'Odorico, Manuela Bischetti, Christopher Cain, Frederick B. Davies, Rebecca L. Davies, Anna-Christina Eilers, Xiaohui Fan, Prakash Gaikwad, Martin G. Haehnelt, Laura C. Keating, Girish Kulkarni, Samuel Lai, Hai-Xia Ma, Andrei Mesinger, Yuxiang Qin, Sindhu Satyavolu, Tsutomu T. Takeuchi, Hideki Umehata, Jinyi Yang

Abstract: The mean free path of ionizing photons, $λ_{\rm mfp}$, is a critical parameter for modeling the intergalactic medium (IGM) both during and after reionization. We present direct measurements of $λ_{\rm mfp}$ from QSO spectra over the redshift range $5<z<6$, including the first measurements at $z\simeq5.3$ and 5.6. Our sample includes data from the XQR-30 VLT large program, as well as new Keck/ESI o… ▽ More The mean free path of ionizing photons, $λ_{\rm mfp}$, is a critical parameter for modeling the intergalactic medium (IGM) both during and after reionization. We present direct measurements of $λ_{\rm mfp}$ from QSO spectra over the redshift range $5<z<6$, including the first measurements at $z\simeq5.3$ and 5.6. Our sample includes data from the XQR-30 VLT large program, as well as new Keck/ESI observations of QSOs near $z \sim 5.5$, for which we also acquire new [C II] 158$μ$m redshifts with ALMA. By measuring the Lyman continuum transmission profile in stacked QSO spectra, we find $λ_{\rm mfp} = 9.33_{-1.80}^{+2.06}$, $5.40_{-1.40}^{+1.47}$, $3.31_{-1.34}^{+2.74}$, and $0.81_{-0.48}^{+0.73}$ pMpc at $z=5.08$, 5.31, 5.65, and 5.93, respectively. Our results demonstrate that $λ_{\rm mfp}$ increases steadily and rapidly with time over $5<z<6$. Notably, we find that $λ_{\rm mfp}$ deviates significantly from predictions based on a fully ionized and relaxed IGM as late as $z=5.3$. By comparing our results to model predictions and indirect $λ_{\rm mfp}$ constraints based on IGM Ly$α$ opacity, we find that the $λ_{\rm mfp}$ evolution is consistent with scenarios wherein the IGM is still undergoing reionization and/or retains large fluctuations in the ionizing UV background well below redshift six. △ Less

Submitted 8 August, 2023; originally announced August 2023.

Comments: 17 pages, 6 figures, 2 tables; accepted for publication in ApJ

CO Multi-line Imaging of Nearby Galaxies (COMING). XII. CO-to-H$_{2}$ Conversion Factor and Dust-to-Gas Ratio

Authors: Atsushi Yasuda, Nario Kuno, Kazuo Sorai, Kazuyuki Muraoka, Yusuke Miyamoto, Hiroyuki Kaneko, Yoshiyuki Yajima, Takahiro Tanaka, Kana Morokuma-Matsui, Tsutomu T. Takeuchi, Masato I. N. Kobayashi

Abstract: We simultaneously measured the spatially-resolved CO-to-H$_{2}$ conversion factor ($α_\mathrm{CO}$) and dust-to-gas ratio (DGR) in nearby galaxies on a kiloparsec scale. In this study, we used $^{12}$CO($J=1-0$) data obtained by the Nobeyama 45-m radio telescope with HI and dust mass surface densities. We obtained the values of global $α_\mathrm{CO}$ and DGR in 22 nearby spiral galaxies, with aver… ▽ More We simultaneously measured the spatially-resolved CO-to-H$_{2}$ conversion factor ($α_\mathrm{CO}$) and dust-to-gas ratio (DGR) in nearby galaxies on a kiloparsec scale. In this study, we used $^{12}$CO($J=1-0$) data obtained by the Nobeyama 45-m radio telescope with HI and dust mass surface densities. We obtained the values of global $α_\mathrm{CO}$ and DGR in 22 nearby spiral galaxies, with averages of $2.66 \pm 1.36\ M_\odot\ \mathrm{pc}^{-2}\ (\mathrm{K\ km\ s^{-1}})^{-1}$ and $0.0052 \pm 0.0026$, respectively. Furthermore, the radial variations of $α_\mathrm{CO}$ and DGR in four barred spiral galaxies (IC 342, NGC 3627, NGC 5236, and NGC 6946) were obtained by dividing them into the inner and outer regions with a boundary of $0.2R_{25}$, where $R_{25}$ is the isophotal radius at 25 mag arcsec$^{-2}$ in the $B$ band. The averages of $α_\mathrm{CO}$ and DGR in the inner region ($\leq 0.2R_{25}$) are $0.36 \pm 0.08\ M_\odot\ \mathrm{pc}^{-2}\ (\mathrm{K\ km\ s^{-1}})^{-1}$ and $0.0199 \pm 0.0058$, while those in the outer region ($> 0.2R_{25}$) are $1.49 \pm 0.76\ M_\odot\ \mathrm{pc}^{-2}\ (\mathrm{K\ km\ s^{-1}})^{-1}$ and $0.0084 \pm 0.0037$, respectively. The value of $α_\mathrm{CO}$ in the outer region is 2.3 to 5.3 times larger than that of the inner region. When separated into the inner and outer regions, we find that $α_\mathrm{CO}$ and DGR correlate with the metallicity and the star formation rate surface density. The value of $α_\mathrm{CO}$ derived in this study tends to be smaller than those obtained in previous studies for the Milky Way and nearby star-forming galaxies. This fact can be attributed to our measurements being biased toward the inner region; we measured $α_\mathrm{CO}$ at 0.85 and 0.76 times smaller in radius than the previous works for nearby star-forming galaxies and the Milky Way, respectively. △ Less

Submitted 12 May, 2023; originally announced May 2023.

Comments: 64 pages, 32 figures, accepted to PASJ

Photometric Observations and Period Analysis of an SU UMa-type Dwarf Nova, MASTER OT J004527.52+503213.8

Authors: Sena A. Matsui, Tsutomu T. Takeuchi, Kai T. Kono, Suchetha Cooray

Abstract: MASTER OT J004527.52+503213.8 (hereafter MASTER J004527) is a dwarf nova discovered by the MASTER project in 2013. At 18:20 UTC on 24 October 2020, brightening of this object was reported to vsnet-alert (24843 by Denisenko). This was the second report of a superoutburst after its discovery. Photometric observations were made using the 23.5-cm Schmidt-Cassegrain telescope at Okayama University of S… ▽ More MASTER OT J004527.52+503213.8 (hereafter MASTER J004527) is a dwarf nova discovered by the MASTER project in 2013. At 18:20 UTC on 24 October 2020, brightening of this object was reported to vsnet-alert (24843 by Denisenko). This was the second report of a superoutburst after its discovery. Photometric observations were made using the 23.5-cm Schmidt-Cassegrain telescope at Okayama University of Science observatory soon after the alert through 4 November 2020. In this work, we present the photometric data from our observation, and the analysis of the light curves of MASTER J004527 during the 2020 outburst. We propose a method to determine the period of superhumps by polynomial fitting, which can be applied to a light curve with many missing data. In addition to our own data, we incorporate other all sky survey data of the outburst to better understand the properties of the superhumps. Based on our observations, we conclude that MASTER J004527 is an SU UMa-type dwarf nova, since no early superhumps occurred. △ Less

Submitted 25 March, 2023; originally announced March 2023.

Comments: Accepted for publication in JAAVSO, 9 pages, 8 figures

EMPRESS. XII. Statistics on the Dynamics and Gas Mass Fraction of Extremely-Metal Poor Galaxies

Authors: Yi Xu, Masami Ouchi, Yuki Isobe, Kimihiko Nakajima, Shinobu Ozaki, Nicolas F. Bouché, John H. Wise, Eric Emsellem, Haruka Kusakabe, Takashi Hattori, Tohru Nagao, Gen Chiaki, Hajime Fukushima, Yuichi Harikane, Kohei Hayashi, Yutaka Hirai, Ji Hoon Kim, Michael V. Maseda, Kentaro Nagamine, Takatoshi Shibuya, Yuma Sugahara, Hidenobu Yajima, Shohei Aoyama, Seiji Fujimoto, Keita Fukushima , et al. (27 additional authors not shown)

Abstract: We present demography of the dynamics and gas-mass fraction of 33 extremely metal-poor galaxies (EMPGs) with metallicities of $0.015-0.195~Z_\odot$ and low stellar masses of $10^4-10^8~M_\odot$ in the local universe. We conduct deep optical integral-field spectroscopy (IFS) for the low-mass EMPGs with the medium high resolution ($R=7500$) grism of the 8m-Subaru FOCAS IFU instrument by the EMPRESS… ▽ More We present demography of the dynamics and gas-mass fraction of 33 extremely metal-poor galaxies (EMPGs) with metallicities of $0.015-0.195~Z_\odot$ and low stellar masses of $10^4-10^8~M_\odot$ in the local universe. We conduct deep optical integral-field spectroscopy (IFS) for the low-mass EMPGs with the medium high resolution ($R=7500$) grism of the 8m-Subaru FOCAS IFU instrument by the EMPRESS 3D survey, and investigate H$α$ emission of the EMPGs. Exploiting the resolution high enough for the low-mass galaxies, we derive gas dynamics with the H$α$ lines by the fitting of 3-dimensional disk models. We obtain an average maximum rotation velocity ($v_\mathrm{rot}$) of $15\pm3~\mathrm{km~s^{-1}}$ and an average intrinsic velocity dispersion ($σ_0$) of $27\pm10~\mathrm{km~s^{-1}}$ for 15 spatially resolved EMPGs out of the 33 EMPGs, and find that all of the 15 EMPGs have $v_\mathrm{rot}/σ_0<1$ suggesting dispersion dominated systems. There is a clear decreasing trend of $v_\mathrm{rot}/σ_0$ with the decreasing stellar mass and metallicity. We derive the gas mass fraction ($f_\mathrm{gas}$) for all of the 33 EMPGs, and find no clear dependence on stellar mass and metallicity. These $v_\mathrm{rot}/σ_0$ and $f_\mathrm{gas}$ trends should be compared with young high-$z$ galaxies observed by the forthcoming JWST IFS programs to understand the physical origins of the EMPGs in the local universe. △ Less

Submitted 26 January, 2024; v1 submitted 22 March, 2023; originally announced March 2023.

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

Journal ref: The Astrophysical Journal, Volume 961, Number 1, January 2024, Page 49-53

The 300 pc resolution imaging of a z = 8.31 galaxy: Turbulent ionized gas and potential stellar feedback 600 million years after the Big Bang

Authors: Yoichi Tamura, Tom J. L. C. Bakx, Akio K. Inoue, Takuya Hashimoto, Tsuyoshi Tokuoka, Chihiro Imamura, Bunyo Hatsukade, Minju M. Lee, Kana Moriwaki, Takashi Okamoto, Kazuaki Ota, Hideki Umehata, Naoki Yoshida, Erik Zackrisson, Masato Hagimoto, Hiroshi Matsuo, Ikkoh Shimizu, Yuma Sugahara, Tsutomu T. Takeuchi

Abstract: We present the results of 300 pc resolution ALMA imaging of the [OIII] 88 $μ$m line and dust continuum emission from a $z = 8.312$ Lyman break galaxy MACS0416_Y1. The velocity-integrated [OIII] emission has three peaks which are likely associated with three young stellar clumps of MACS0416_Y1, while the channel map shows a complicated velocity structure with little indication of a global velocity… ▽ More We present the results of 300 pc resolution ALMA imaging of the [OIII] 88 $μ$m line and dust continuum emission from a $z = 8.312$ Lyman break galaxy MACS0416_Y1. The velocity-integrated [OIII] emission has three peaks which are likely associated with three young stellar clumps of MACS0416_Y1, while the channel map shows a complicated velocity structure with little indication of a global velocity gradient unlike what was found in [CII] 158 $μ$m at a larger scale, suggesting random bulk motion of ionized gas clouds inside the galaxy. In contrast, dust emission appears as two individual clumps apparently separating or bridging the [OIII]/stellar clumps. The cross correlation coefficient between dust and ultraviolet-related emission (i.e., [OIII] and ultraviolet continuum) is unity on a galactic scale, while it drops at < 1 kpc, suggesting well mixed geometry of multi-phase interstellar media on sub-kpc scales. If the cutoff scale characterizes different stages of star formation, the cutoff scale can be explained by gravitational instability of turbulent gas. We also report on a kpc-scale off-center cavity embedded in the dust continuum image. This could be a superbubble producing galactic-scale outflows, since the energy injection from the 4 Myr starburst suggested by a spectral energy distribution analysis is large enough to push the surrounding media creating a kpc-scale cavity. △ Less

Submitted 16 July, 2023; v1 submitted 20 March, 2023; originally announced March 2023.

Comments: 18 pages, 6 figures, 1 table, published in ApJ

The SCUBA-2 Large eXtragalactic Survey: 850um map, catalogue and the bright-end number counts of the XMM-LSS field

Authors: T. K. Garratt, J. E. Geach, Y. Tamura, K. E. K. Coppin, M. Franco, Y. Ao, C. -C. Chen, C. Cheng, D. L. Clements, Y. S. Dai, H. Dannerbauer, T. R. Greve, B. Hatsukade, H. S. Hwang, L. Jiang, K. Kohno, M. P. Koprowski, M. J. Michalowski, M. Sawicki, D. Scott, H. Shim, T. T. Takeuchi, W. -H. Wang, Y. Q. Xue, C. Yang

Abstract: We present 850um imaging of the XMM-LSS field observed for 170 hours as part of the James Clerk Maxwell Telescope SCUBA-2 Large eXtragalactic Survey (S2LXS). S2LXS XMM-LSS maps an area of 9 square degrees, reaching a moderate depth of 1-sigma ~ 4 mJy/beam. This is the largest contiguous area of extragalactic sky mapped by JCMT at 850um to date. The wide area of the S2LXS XMM-LSS survey allows us t… ▽ More We present 850um imaging of the XMM-LSS field observed for 170 hours as part of the James Clerk Maxwell Telescope SCUBA-2 Large eXtragalactic Survey (S2LXS). S2LXS XMM-LSS maps an area of 9 square degrees, reaching a moderate depth of 1-sigma ~ 4 mJy/beam. This is the largest contiguous area of extragalactic sky mapped by JCMT at 850um to date. The wide area of the S2LXS XMM-LSS survey allows us to probe the ultra-bright (S_850um > 15 mJy), yet rare submillimetre population. We present the S2LXS XMM-LSS catalogue, which comprises 40 sources detected at >5-sigma significance, with deboosted flux densities in the range of 7 mJy to 48 mJy. We robustly measure the bright-end of the 850um number counts at flux densities >7 mJy, reducing the Poisson errors compared to existing measurements. The S2LXS XMM-LSS observed number counts show the characteristic upturn at bright fluxes, expected to be motivated by local sources of submillimetre emission and high-redshift strongly lensed galaxies. We find that the observed 850um number counts are best reproduced by model predictions that include either strong lensing or source blending from a 15 arcsec beam, indicating that both may make an important contribution to the observed over-abundance of bright single dish 850um selected sources. We make the S2LXS XMM-LSS 850um map and >5-sigma catalogue presented here publicly available. △ Less

Submitted 25 January, 2023; originally announced January 2023.

Comments: Accepted for publication in MNRAS

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

Galaxy Manifold: Characterizing and understanding galaxies with two parameters

Authors: Suchetha Cooray, Tsutomu T. Takeuchi, Daichi Kashino, Shuntaro A. Yoshida, Hai-Xia Ma, Kai T. Kono

Abstract: We report the discovery of a two-dimensional Galaxy Manifold within the multi-dimensional luminosity space of local galaxies. The multi-dimensional luminosity space is constructed using 11 bands that span from far ultraviolet to near-infrared for redshift < 0.1 galaxies observed with GALEX, SDSS, and UKIDSS. The two latent parameters are sufficient to express 93.2% of the variance in the galaxy sa… ▽ More We report the discovery of a two-dimensional Galaxy Manifold within the multi-dimensional luminosity space of local galaxies. The multi-dimensional luminosity space is constructed using 11 bands that span from far ultraviolet to near-infrared for redshift < 0.1 galaxies observed with GALEX, SDSS, and UKIDSS. The two latent parameters are sufficient to express 93.2% of the variance in the galaxy sample, suggesting that this Galaxy Manifold is one of the most efficient representations of galaxies. The transformation between the observed luminosities and the manifold parameters as an analytic mapping is provided. The manifold representation provides accurate (85%) morphological classifications with a simple linear boundary, and galaxy properties can be estimated with minimal scatter (0.12 dex and 0.04 dex for star formation rate and stellar mass, respectively) by calibrating with the two-dimensional manifold location. Under the assumption that the manifold expresses the possible parameter space of galaxies, the evolution on the manifold is considered. We find that constant and exponentially decreasing star formation histories form almost orthogonal modes of evolution on the manifold. Through these simple models, we understand that the two modes are closely related to gas content, which suggests the close relationship of the manifold to gas accretion. Without assuming a star formation history, a gas-regulated model reproduces an exponentially declining star formation history with a timescale of $\sim$1.2 Gyrs on the manifold. Lastly, the found manifold suggests a paradigm where galaxies are characterized by their mass/scale and specific SFR, which agrees with previous studies of dimensionality reduction. △ Less

Submitted 11 October, 2022; originally announced October 2022.

Comments: Submitted to MNRAS. 15 pages and 17 figures. All comments are welcome

Deblurring galaxy images with Tikhonov regularization on magnitude domain

Authors: Kazumi Murata, Tsutomu T. Takeuchi

Abstract: We propose a regularization-based deblurring method that works efficiently for galaxy images. The spatial resolution of a ground-based telescope is generally limited by seeing conditions and much worse than space-based telescopes. This circumstance has generated considerable research interest in restoration of spatial resolution. Since image deblurring is a typical inverse problem and often ill-po… ▽ More We propose a regularization-based deblurring method that works efficiently for galaxy images. The spatial resolution of a ground-based telescope is generally limited by seeing conditions and much worse than space-based telescopes. This circumstance has generated considerable research interest in restoration of spatial resolution. Since image deblurring is a typical inverse problem and often ill-posed, solutions tend to be unstable. To obtain a stable solution, much research has adopted regularization-based methods for image deblurring, but the regularization term is not necessarily appropriate for galaxy images. Although galaxies have an exponential or Sersic profile, the conventional regularization assumes the image profiles to behave linear in space. The significant deviation between the assumption and real situation leads to blurring the images and smoothing out the detailed structures. Clearly, regularization on logarithmic, i.e. magnitude domain, should provide a more appropriate assumption, which we explore in this study. We formulate a problem of deblurring galaxy images by an objective function with a Tikhonov regularization term on magnitude domain. We introduce an iterative algorithm minimizing the objective function with a primal-dual splitting method. We investigate the feasibility of the proposed method using simulation and observation images. In the simulation, we blur galaxy images with a realistic point spread function and add both Gaussian and Poisson noises. For the evaluation with the observed images, we use galaxy images taken by the Subaru HSC-SSP. Both of these evaluations show that our method successfully recovers the spatial resolution of the images and significantly outperforms the conventional methods. The code is publicly available from the Github ( https://github.com/kzmurata-astro/PSFdeconv_amag ). △ Less

Submitted 22 August, 2022; originally announced August 2022.

Comments: 14 pages, 9 figures, accepted for publication in PASJ. The code is available at https://github.com/kzmurata-astro/PSFdeconv_amag

Dissecting Nearby Galaxies with piXedfit: II. Spatially Resolved Scaling Relations Among Stars, Dust, and Gas

Authors: Abdurro'uf, Yen-Ting Lin, Hiroyuki Hirashita, Takahiro Morishita, Sandro Tacchella, Po-Feng Wu, Masayuki Akiyama, Tsutomu T. Takeuchi

Abstract: We study spatially resolved scaling relations among stars, dust, and gas in ten nearby spiral galaxies. In a preceding paper Abdurro'uf et al. (2022), we have derived spatially resolved properties of the stellar population and dust by panchromatic spectral energy distribution (SED) fitting using piXedfit. Now, we investigate resolved star formation ($Σ_{\rm H_{2}}$--$Σ_{\rm SFR}$--$Σ_{*}$) and dus… ▽ More We study spatially resolved scaling relations among stars, dust, and gas in ten nearby spiral galaxies. In a preceding paper Abdurro'uf et al. (2022), we have derived spatially resolved properties of the stellar population and dust by panchromatic spectral energy distribution (SED) fitting using piXedfit. Now, we investigate resolved star formation ($Σ_{\rm H_{2}}$--$Σ_{\rm SFR}$--$Σ_{*}$) and dust scaling relations. While the relations with all sub-galactic regions of the galaxies are reasonably tight ($σ\lesssim 0.3$ dex), we find that most of the scaling relations exhibit galaxy-to-galaxy variations in normalization and shape. Only two relations of $Σ_{\rm dust}$--$Σ_{\rm gas}$ and $Σ_{\rm dust}$--$Σ_{\rm H_{2}}$ do not show noticeable galaxy-to-galaxy variations among our sample galaxies. We further investigate correlations among the scaling relations. We find significant correlations among the normalization of the $Σ_{\rm H_{2}}$--$Σ_{\rm SFR}$--$Σ_{*}$ relations, which suggest that galaxies with higher levels of resolved $\text{H}_{2}$ fraction ($f_{\rm H_{2}}$) tend to have higher levels of resolved star formation efficiency (SFE) and specific star formation rate (sSFR). We also observe that galaxies with higher levels of resolved dust-to-stellar mass ratios tend to have higher levels of resolved sSFR, SFE, and $f_{\rm H_{2}}$. Moreover, we find that galaxies with higher global sSFR and less compact morphology tend to have higher levels of the resolved sSFR, SFE, and $f_{\rm H_{2}}$, which can explain the variations in the normalization of the $Σ_{\rm H_{2}}$--$Σ_{\rm SFR}$--$Σ_{*}$ relationships. Overall, we observe indications of the contributions of both global and local factors in governing the star formation process in galaxies. △ Less

Submitted 18 July, 2022; originally announced July 2022.

Comments: 28 pages, 13 figures, accepted for publication in ApJ. piXedfit v1.0 is publicly available at https://github.com/aabdurrouf/piXedfit. Documentation is available at https://pixedfit.readthedocs.io/en/latest/. Some tutorials of practical usages in jupyter notebooks can be found at https://github.com/aabdurrouf/piXedfit/tree/main/examples

EMPRESS. IX. Extremely Metal-Poor Galaxies are Very Gas-Rich Dispersion-Dominated Systems: Will JWST Witness Gaseous Turbulent High-z Primordial Galaxies?

Authors: Yuki Isobe, Masami Ouchi, Kimihiko Nakajima, Shinobu Ozaki, Nicolas F. Bouche, John H. Wise, Yi Xu, Eric Emsellem, Haruka Kusakabe, Takashi Hattori, Tohru Nagao, Gen Chiaki, Hajime Fukushima, Yuichi Harikane, Kohei Hayashi, Yutaka Hirai, Ji Hoon Kim, Michael V. Maseda, Kentaro Nagamine, Takatoshi Shibuya, Yuma Sugahara, Hidenobu Yajima, Shohei Aoyama, Seiji Fujimoto, Keita Fukushima , et al. (27 additional authors not shown)

Abstract: We present kinematics of 6 local extremely metal-poor galaxies (EMPGs) with low metallicities ($0.016-0.098\ Z_{\odot}$) and low stellar masses ($10^{4.7}-10^{7.6} M_{\odot}$). Taking deep medium-high resolution ($R\sim7500$) integral-field spectra with 8.2-m Subaru, we resolve the small inner velocity gradients and dispersions of the EMPGs with H$α$ emission. Carefully masking out sub-structures… ▽ More We present kinematics of 6 local extremely metal-poor galaxies (EMPGs) with low metallicities ($0.016-0.098\ Z_{\odot}$) and low stellar masses ($10^{4.7}-10^{7.6} M_{\odot}$). Taking deep medium-high resolution ($R\sim7500$) integral-field spectra with 8.2-m Subaru, we resolve the small inner velocity gradients and dispersions of the EMPGs with H$α$ emission. Carefully masking out sub-structures originated by inflow and/or outflow, we fit 3-dimensional disk models to the observed H$α$ flux, velocity, and velocity-dispersion maps. All the EMPGs show rotational velocities ($v_{\rm rot}$) of 5--23 km s$^{-1}$ smaller than the velocity dispersions ($σ_{0}$) of 17--31 km s$^{-1}$, indicating dispersion-dominated ($v_{\rm rot}/σ_{0}=0.29-0.80<1$) systems affected by inflow and/or outflow. Except for two EMPGs with large uncertainties, we find that the EMPGs have very large gas-mass fractions of $f_{\rm gas}\simeq 0.9-1.0$. Comparing our results with other H$α$ kinematics studies, we find that $v_{\rm rot}/σ_{0}$ decreases and $f_{\rm gas}$ increases with decreasing metallicity, decreasing stellar mass, and increasing specific star-formation rate. We also find that simulated high-$z$ ($z\sim 7$) forming galaxies have gas fractions and dynamics similar to the observed EMPGs. Our EMPG observations and the simulations suggest that primordial galaxies are gas-rich dispersion-dominated systems, which would be identified by the forthcoming James Webb Space Telescope (JWST) observations at $z\sim 7$. △ Less

Submitted 19 April, 2023; v1 submitted 9 June, 2022; originally announced June 2022.

Comments: ApJ in Press

A new galaxy spectral energy distribution model consistent with the evolution of dust

Authors: Kazuki Y. Nishida, Tsutomu T. Takeuchi, Takuma Nagata, Ryosuke S. Asano

Abstract: The spectral energy distribution (SED) of galaxies provides fundamental information on the related physical processes. However, the SED is significantly affected by dust in its interstellar medium. Dust is mainly produced by asymptotic giant branch stars and Type II supernovae. In addition, the dust mass increases through the metal accretion, and the grain size changes by the collisions between th… ▽ More The spectral energy distribution (SED) of galaxies provides fundamental information on the related physical processes. However, the SED is significantly affected by dust in its interstellar medium. Dust is mainly produced by asymptotic giant branch stars and Type II supernovae. In addition, the dust mass increases through the metal accretion, and the grain size changes by the collisions between the grains. The contribution of each process and the extinction depend on the size distribution. Therefore, the SED model should treat the evolution of the dust mass and size distribution. In spite of the importance of dust evolution, many previous SED models have not considered the evolution of the total mass and size distribution in a physically consistent manner. In this work, we constructed a new radiative transfer SED model, based on our dust evolution model consistent with the chemical evolution. To reduce the computational cost, we adopted the mega-grain and the one-dimensional plane parallel galaxy approximation. As a fiducial case, we calculated Milky Way-like galaxy SEDs at various ages under the closed-box model. We found that a galaxy at the age of 100~Myr does not produce small grains such as polycyclic aromatic hydrocarbons. After 1~Gyr, we observed a drastic increase of infrared emission and attenuation caused by a rapid increase of dust mass. This phenomenon can be treated appropriately for the first time by our new model. This model can be used for the SED fitting to a galaxy at any stage of evolution. △ Less

Submitted 16 May, 2022; originally announced May 2022.

Comments: 19 Pages, 19 Figures, accepted for publication in MNRAS

Estimation of the Star Formation Rate of Galaxies with Radio Continuum Obtained with Murchison Widefield Array

Authors: Tsutomu T. Takeuchi, Shuntaro A. Yoshida, Luca Cortese, O. Ivy Wong, Barbara Catinella, Suchetha Cooray

Abstract: We investigate the correlation between the integrated low-frequency and infrared (IR) emissions of star-forming galaxies extracted from the {\sl Herschel} Reference Survey. By taking advantage of the GaLactic Extragalactic All-sky MWA (GLEAM) survey operated by the Murchison Widefield Array (MWA) we examine how this correlation varies at a function of frequency across the 20 GLEAM narrow bands at… ▽ More We investigate the correlation between the integrated low-frequency and infrared (IR) emissions of star-forming galaxies extracted from the {\sl Herschel} Reference Survey. By taking advantage of the GaLactic Extragalactic All-sky MWA (GLEAM) survey operated by the Murchison Widefield Array (MWA) we examine how this correlation varies at a function of frequency across the 20 GLEAM narrow bands at $72\mbox{--}231\; [\mbox{MHz}]$. These examinations are important for ensuring the reliability of the radio luminosity as a SFR indicator. In this study, we focus on 18 star-forming galaxies whose radio emission is detected by the GLEAM survey. These galaxies show that a single power-law is sufficient to characterise the far-infrared-to-radio correlation across the GLEAM frequency bands and up to $1.5\; [\mbox{GHz}]$. Thus, the radio continuum in this wavelength range can serve as a good dust extinction-free SFR estimator. This is particularly important for future investigation of the cosmic SFR independently from other estimators, since the radio continuum can be detected from $z=0$ to high redshifts ($z \sim 5\mbox{--}10$) in a coherent manner. △ Less

Submitted 2 April, 2022; originally announced April 2022.

Comments: 13 pages, 6 figures, 2 tables, submitted

EMPRESS. VIII. A New Determination of Primordial He Abundance with Extremely Metal-Poor Galaxies: A Suggestion of the Lepton Asymmetry and Implications for the Hubble Tension

Authors: Akinori Matsumoto, Masami Ouchi, Kimihiko Nakajima, Masahiro Kawasaki, Kai Murai, Kentaro Motohara, Yuichi Harikane, Yoshiaki Ono, Kosuke Kushibiki, Shuhei Koyama, Shohei Aoyama, Masahiro Konishi, Hidenori Takahashi, Yuki Isobe, Hiroya Umeda, Yuma Sugahara, Masato Onodera, Kentaro Nagamine, Haruka Kusakabe, Yutaka Hirai, Takashi J. Moriya, Takatoshi Shibuya, Yutaka Komiyama, Keita Fukushima, Seiji Fujimoto , et al. (20 additional authors not shown)

Abstract: The primordial He abundance $Y_\mathrm{P}$ is a powerful probe of cosmology. Currently, $Y_\mathrm{P}$ is best determined by observations of metal-poor galaxies, while there are only a few known local extremely metal-poor ($<0.1 Z_\odot$) galaxies (EMPGs) having reliable He/H measurements with HeI$λ$10830 near-infrared (NIR) emission. Here we present deep Subaru NIR spectroscopy for 10 EMPGs. Comb… ▽ More The primordial He abundance $Y_\mathrm{P}$ is a powerful probe of cosmology. Currently, $Y_\mathrm{P}$ is best determined by observations of metal-poor galaxies, while there are only a few known local extremely metal-poor ($<0.1 Z_\odot$) galaxies (EMPGs) having reliable He/H measurements with HeI$λ$10830 near-infrared (NIR) emission. Here we present deep Subaru NIR spectroscopy for 10 EMPGs. Combining the existing optical data, He/H values of 5 out of the 10 EMPGs are reliably derived by the Markov chain Monte Carlo algorithm. Adding the existing 3 EMPGs and 51 moderately metal-poor ($0.1-0.4 Z_\odot$) galaxies with reliable He/H estimates, we obtain $Y_\mathrm{P}=0.2370^{+0.0034}_{-0.0033}$ by linear regression in the $\mathrm{(He/H)}-\mathrm{(O/H)}$ plane, where we increase the number of EMPGs from 3 to 8 anchoring He/H of the most metal-poor gas in galaxies. Although our $Y_\mathrm{P}$ measurement and previous measurements are consistent, our result is slightly ($\sim 1σ$) smaller due to our EMPGs. With our $Y_\mathrm{P}$ and the existing primordial deuterium $D_\mathrm{P}$ measurement, we constrain the effective number of neutrino species $N_\mathrm{eff}$ and the baryon-to-photon ratio $η$ showing $\gtrsim 1-2σ$ tensions with the Standard Model and Planck Collaboration et al. (2020). Motivated by the tensions, we allow the degeneracy parameter of electron-neutrino $ξ_e$ to vary as well as $N_\mathrm{eff}$ and $η$. We obtain $ξ_e = 0.05^{+0.03}_{-0.02}$, $N_\mathrm{eff}=3.11^{+0.34}_{-0.31}$, and $η\times10^{10}=6.08^{+0.06}_{-0.06}$ from the $Y_\mathrm{P}$ and $D_\mathrm{P}$ measurements with a prior of $η$ taken from Planck Collaboration et al. (2020). Our constraints suggest a lepton asymmetry and allow for a high value of $N_\mathrm{eff}$ within the $1σ$ level, which could mitigate the Hubble tension. △ Less

Submitted 27 November, 2022; v1 submitted 17 March, 2022; originally announced March 2022.

Comments: 22 pages, 9 figures, accepted for publication in ApJ

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)

The bimodality in the mass-metallicity relation in SDSS-MaNGA galaxy pairs

Authors: Kiyoaki Christopher Omori, Tsutomu T. Takeuchi

Abstract: Aims: Interacting galaxies show a metallicity dilution compared to isolated galaxies of similar masses in the mass-metallicity space at the global scale. We investigate the spatially resolved mass-metallicity relation (MZR) of galaxy pairs in the SDSS-MaNGA survey to confirm that the local relation between the stellar mass surface density and the metallicity is consistent with the MZR at the globa… ▽ More Aims: Interacting galaxies show a metallicity dilution compared to isolated galaxies of similar masses in the mass-metallicity space at the global scale. We investigate the spatially resolved mass-metallicity relation (MZR) of galaxy pairs in the SDSS-MaNGA survey to confirm that the local relation between the stellar mass surface density and the metallicity is consistent with the MZR at the global scale. Methods: We investigate the relationship between the stellar mass surface density and the metallicity abundance, 12+log (O/H), for star-forming spaxels belonging to 298 galaxy pairs identified using visual and kinematic indicators in the SDSS-MaNGA survey. We also investigate if a)the location of a spaxel relative to the galaxy centre and b)the galaxy pair separation have any effect on the local MZR. Results: We find that the correlation between mass and metallicity holds for interacting galaxies at the local level. However, we find two peaks in spaxel distribution, one peak with enriched metallicity and the other with diluted metallicity. We find that the spaxels belonging to the galaxy central regions (i.e. at lower R/Reff) are concentrated close to the two peaks. We also find that the metallicity-diluted spaxels belong to galaxy pairs with closer projected separations and that spaxels with enriched metallicity belong to galaxy pairs with greater projected separations. Conclusions: We find two discrete peaks in the spatially resolved MZR for star-forming spaxels that belong to galaxy pairs. The peaks are likely related to the galaxy projected separation or the stage of the interaction process of a galaxy pair. △ Less

Submitted 28 February, 2022; v1 submitted 7 December, 2021; originally announced December 2021.

Comments: 9 pages, 6 figures, accepted by A&A Feb. 17, 2022

Journal ref: A&A 660, A145 (2022)

An Empirical Estimation of the Galaxy Radio Number Counts model at 3~GHz and the Effect of Gravitational Lensing

Authors: K. T. Kono, T. T. Takeuchi

Abstract: We analyzed the evolution of 3GHz luminosity function of radio selected galaxies through cosmic time, up to about 1 Gyr after Big Bang, with making use of the VLA-COSMOS3GHz large project. We estimated the luminosity function for SFGs and AGNs in the VLA-COSMOS field respectively with two independent methods, the non-parametric C-method for accuracy to avoid the effect of cosmic variance and param… ▽ More We analyzed the evolution of 3GHz luminosity function of radio selected galaxies through cosmic time, up to about 1 Gyr after Big Bang, with making use of the VLA-COSMOS3GHz large project. We estimated the luminosity function for SFGs and AGNs in the VLA-COSMOS field respectively with two independent methods, the non-parametric C-method for accuracy to avoid the effect of cosmic variance and parametric MCMC method assuming pure luminosity evolution. We obtained pure luminosity evolution parameter as L*=(1+z)(3.36-0.31)z for SFG and L* = (1+z)(2.91-0.99)z for AGN. The resultant model agreed well with observed number counts down to sub-muJy level. Additionally, we compared the model with cosmic star formation rate density history with multiple observations in various wavelength and confirmed that our result is consistent with literature up to z=3. We further analyzed the statistical effect of gravitational lensing due to the dark matter halos and found that the effect is maximized at S_3GHz = 100muJy with about 1% contribution for SFG and less than 0.5% for AGN for the flux range considered in this study. In order to evaluate the effect of lensing magnification on parameter estimation in the SKA observations, we performed Fisher analysis and find that while the bias is limited, there is a distinguishable effect on SFG LF evolution parameter with SKA I MID wide survey, and all-sky and wide survey for AGN LF evolution parameters. △ Less

Submitted 2 December, 2021; originally announced December 2021.

Comments: 12 pages, 7 figures

Wavelets and sparsity for Faraday tomography

Authors: Suchetha Cooray, Tsutomu T. Takeuchi, Shinsuke Ideguchi, Takuya Akahori, Yoshimitsu Miyashita, Keitaro Takahashi

Abstract: Faraday tomography through broadband polarimetry can provide crucial information on magnetized astronomical objects, such as quasars, galaxies, or galaxy clusters. However, the limited wavelength coverage of the instruments requires that we solve an ill-posed inverse problem when we want to obtain the Faraday dispersion function (FDF), a tomographic distribution of the magnetoionic media along the… ▽ More Faraday tomography through broadband polarimetry can provide crucial information on magnetized astronomical objects, such as quasars, galaxies, or galaxy clusters. However, the limited wavelength coverage of the instruments requires that we solve an ill-posed inverse problem when we want to obtain the Faraday dispersion function (FDF), a tomographic distribution of the magnetoionic media along the line of sight. This paper explores the use of wavelet transforms and the sparsity of the transformed FDFs in the form of wavelet shrinkage (WS) for finding better solutions to the inverse problem. We recently proposed the Constraining and Restoring iterative Algorithm for Faraday Tomography (CRAFT; Cooray et al. 2021), a new flexible algorithm that showed significant improvements over the popular methods such as Rotation Measure Synthesis. In this work, we introduce CRAFT+WS, a new version of CRAFT incorporating the ideas of wavelets and sparsity. CRAFT+WS exhibit significant improvements over the original CRAFT when tested for a complex FDF of realistic Galactic model. Reconstructions of FDFs demonstrate super-resolution in Faraday depth, uncovering previously unseen Faraday complexities in observations. The proposed approach will be necessary for effective cosmic magnetism studies using the Square Kilometre Array and its precursors. △ Less

Submitted 6 May, 2022; v1 submitted 2 December, 2021; originally announced December 2021.

Comments: 10 Pages, 3 Figures, Submitted to PASJ. Code will be made available upon publication

Dissecting Nearby Galaxies with piXedfit: I. Spatially Resolved Properties of Stars, Dust, and Gas as Revealed by Panchromatic SED Fitting

Authors: Abdurro'uf, Yen-Ting Lin, Hiroyuki Hirashita, Takahiro Morishita, Sandro Tacchella, Masayuki Akiyama, Tsutomu T. Takeuchi, Po-Feng Wu

Abstract: We study spatially resolved properties (on spatial scales of $\sim 1-2$ kpc out to at least $3$ effective radii) of the stars, dust, and gas in ten nearby spiral galaxies. The properties of the stellar population and dust are derived by fitting the spatially resolved spectral energy distribution (SED) with more than 20 photometric bands ranging from far-ultraviolet to far-infrared. Our newly devel… ▽ More We study spatially resolved properties (on spatial scales of $\sim 1-2$ kpc out to at least $3$ effective radii) of the stars, dust, and gas in ten nearby spiral galaxies. The properties of the stellar population and dust are derived by fitting the spatially resolved spectral energy distribution (SED) with more than 20 photometric bands ranging from far-ultraviolet to far-infrared. Our newly developed software piXedfit performs point spread function matching of images, pixel binning, and models the stellar light, dust attenuation, dust emission, and emission from a dusty torus heated by an active galactic nucleus simultaneously through the energy balance approach. With this self-consistent analysis, we present the spatially resolved version of the IRX--$β$ relation, finding that it is consistent with the relationship from the integrated photometry. We show that the old stellar populations contribute to the dust heating, which causes an overestimation of star formation rate (SFR) derived from the total ultraviolet and infrared luminosities on kpc scales. With archival high-resolution maps of atomic and molecular gas, we study the radial variation of the properties of the stellar populations (including stellar mass, age, metallicity, and SFR), dust (including dust mass, dust temperature, and abundance of polycyclic aromatic hydrocarbon), gas, as well as dust-to-stellar mass and dust-to-gas mass ratios. We observe a depletion of molecular gas mass fraction in the central region of the majority of the galaxies, suggesting that the lack of available fuel is an important factor in suppressing the specific SFR at the center. △ Less

Submitted 14 December, 2021; v1 submitted 6 October, 2021; originally announced October 2021.

Comments: 35 pages, 19 figures, 4 tables. Accepted by ApJ (results remain the same). Complementary figures to Fig. 4 (for 10 galaxies) and high resolution version of Fig. 8 are available at https://github.com/aabdurrouf/maps-nearby-galaxies

Jets from MRC 0600-399 bent by magnetic fields in the cluster Abell 3376

Authors: James O. Chibueze, Haruka Sakemi, Takumi Ohmura, Mami Machida, Hiroki Akamatsu, Takuya Akahori, Hiroyuki Nakanishi, Viral Parekh, Ruby van Rooyen, Tsutomu T. Takeuchi

Abstract: Galaxy clusters are known to harbour magnetic fields. The nature of the intra-cluster magnetic fields remains an unresolved question. Intra-cluster magnetic field can be observed at the density contact discontinuity formed by cool and dense plasma running into hot ambient plasma, and the discontinuity exists near the 2nd BCG MRC 0600-399 of a merging galaxy cluster Abell 3376 (z=0.0461, hereafter… ▽ More Galaxy clusters are known to harbour magnetic fields. The nature of the intra-cluster magnetic fields remains an unresolved question. Intra-cluster magnetic field can be observed at the density contact discontinuity formed by cool and dense plasma running into hot ambient plasma, and the discontinuity exists near the 2nd BCG MRC 0600-399 of a merging galaxy cluster Abell 3376 (z=0.0461, hereafter as A3376). Elongated X-ray image in the east-west direction with a comet-like structure reaches a Mpc-scale (Fig1.(a)). Previous radio observations detected the bent jets from MRC 0600-399, moving in same direction as the sub-cluster's motion against ram pressure.Here we report a new radio observation of a radio galaxy MRC 0600-399 which is 3.4 times and 11 times higher resolution and sensitivity than the previous results. Contrary to typical jets, the MRC 0600-399 shows a 90deg bend at the contact discontinuity and the collimated jets further extend over 100 kpc from the bend point. Diffuse, elongated emission named "double-scythe" structures were detected for the first time. The spectral index flattens downstream of the bend point, indicating cosmic-ray re-acceleration. High-resolution numerical simulations reveal that the ordered magnetic field along the discontinuity plays a significant role in the change in the jet direction. The morphology of the "double-scythe" bear remarkable similarities with the simulations, which strengthens our understanding of the interaction between relativistic electrons and intra-cluster magnetic field. △ Less

Submitted 24 June, 2021; originally announced June 2021.

Comments: 4 figures

Journal ref: Nature 593, 47-50 (2021). (https://www.nature.com/articles/s41586-021-03434-1)

Constraining the luminosity function of active galactic nuclei through the reionization observations in the SKA era

Authors: Kai T. Kono, Taichi Takeuchi, Hiroyuki Tashiro, Kiyotomo Ichiki, Tsutomu T. Takeuchi

Abstract: Ultraviolet (UV) and X-ray photons from active galactic nuclei (AGNs) can ionize hydrogen in the intergalactic medium (IGM). We solve radiative transfer around AGNs in high redshift to evaluate the 21-cm line emission from the neutral hydrogen in the IGM and obtain the radial profile of the brightness temperature in the epoch of reionization. The ionization profile extends over 10 [Mpc] comoving d… ▽ More Ultraviolet (UV) and X-ray photons from active galactic nuclei (AGNs) can ionize hydrogen in the intergalactic medium (IGM). We solve radiative transfer around AGNs in high redshift to evaluate the 21-cm line emission from the neutral hydrogen in the IGM and obtain the radial profile of the brightness temperature in the epoch of reionization. The ionization profile extends over 10 [Mpc] comoving distance which can be observed in the order of 10 [arcmin]. From estimation of the radio galaxy number counts with high sensitivity observation through the Square Kilometre Array (SKA), we investigate the capability of parameter constrains for AGN luminosity function with Fisher analysis for three evolution model through cosmic time. We find that the errors for each parameter are restricted to a few percent when AGNs are sufficiently bright at high redshifts. We also investigate the possibility of further parameter constraints with future observation beyond the era of SKA. △ Less

Submitted 12 April, 2021; originally announced April 2021.

Comments: 12 pages ,9 figures and 2 tables

Dynamical Evolution of Voids with Surrounding Gravitational Tidal Field

Authors: Mutsumi Minoguchi, Atsushi J. Nishizawa, Tsutomu T. Takeuchi, Naoshi Sugiyama

Abstract: The void ellipticity distribution today can be well explained by the tidal field. Going a step further from the overall distribution, we investigate individuality on the tidal response of void shape in non-linear dynamical evolution. We perform an N-body simulation and trace individual voids using particle ID. The voids are defined based on Voronoi tessellation and watershed algorithm, using publi… ▽ More The void ellipticity distribution today can be well explained by the tidal field. Going a step further from the overall distribution, we investigate individuality on the tidal response of void shape in non-linear dynamical evolution. We perform an N-body simulation and trace individual voids using particle ID. The voids are defined based on Voronoi tessellation and watershed algorithm, using public code VIDE. A positive correlation is found between the time variation of void ellipticity and tidal field around a void if the void maintains its constituent particles. Such voids tend to have smaller mass densities. Conversely, not a few voids significantly deform by particle exchange, rather than the tidal field. Those voids may prevent us from correctly proving a quadrupole field of gravity out of a void shape. △ Less

Submitted 2 March, 2021; originally announced March 2021.

Comments: 10 pages, 12 figures

CO Multi-line Imaging of Nearby Galaxies (COMING). IX. $^{12}$CO($J$=2-1)/$^{12}$CO($J$=1-0) line ratio on kiloparsec scales

Authors: Yoshiyuki Yajima, Kazuo Sorai, Yusuke Miyamoto, Kazuyuki Muraoka, Nario Kuno, Hiroyuki Kaneko, Tsutomu T. Takeuchi, Atsushi Yasuda, Takahiro Tanaka, Kana Morokuma-Matsui, Masato I. N. Kobayashi

Abstract: While molecular gas mass is usually derived from $^{12}$CO($J$=1-0) - the most fundamental line to explore molecular gas - it is often derived from $^{12}$CO($J$=2-1) assuming a constant $^{12}$CO($J$=2-1)/$^{12}$CO($J$=1-0) line ratio ($R_{2/1}$). We present variations of $R_{2/1}$ and effects of the assumption that $R_{2/1}$ is a constant in 24 nearby galaxies using $^{12}$CO data obtained with… ▽ More While molecular gas mass is usually derived from $^{12}$CO($J$=1-0) - the most fundamental line to explore molecular gas - it is often derived from $^{12}$CO($J$=2-1) assuming a constant $^{12}$CO($J$=2-1)/$^{12}$CO($J$=1-0) line ratio ($R_{2/1}$). We present variations of $R_{2/1}$ and effects of the assumption that $R_{2/1}$ is a constant in 24 nearby galaxies using $^{12}$CO data obtained with the Nobeyama 45-m radio telescope and IRAM 30-m telescope. The median of $R_{2/1}$ for all galaxies is 0.61, and the weighted mean of $R_{2/1}$ by $^{12}$CO($J$=1-0) integrated-intensity is 0.66 with a standard deviation of 0.19. The radial variation of $R_{2/1}$ shows that it is high (~0.8) in the inner ~1 kpc while its median in disks is nearly constant at 0.60 when all galaxies are compiled. In the case that the constant $R_{2/1}$ of 0.7 is adopted, we found that the total molecular gas mass derived from $^{12}$CO($J$=2-1) is underestimated/overestimated by ~20%, and at most by 35%. The scatter of a molecular gas surface density within each galaxy becomes larger by ~30%, and at most by 120%. Indices of the spatially resolved Kennicutt-Schmidt relation by $^{12}$CO($J$=2-1) are underestimated by 10-20%, at most 39% in 17 out of 24 galaxies. $R_{2/1}$ has good positive correlations with star-formation rate and infrared color, and a negative correlation with molecular gas depletion time. There is a clear tendency of increasing $R_{2/1}$ with increasing kinetic temperature ($T_{\rm kin}$). Further, we found that not only $T_{\rm kin}$ but also pressure of molecular gas is important to understand variations of $R_{2/1}$. Special considerations should be made when discussing molecular gas mass and molecular gas properties inferred from $^{12}$CO($J$=2-1) instead of $^{12}$CO($J$=1-0). △ Less

Submitted 15 December, 2020; originally announced December 2020.

Comments: 29 pages, 18 figures, 5 tables. Accepted for publication in PASJ. The original resolution version is available here (https://astro3.sci.hokudai.ac.jp/~radio/coming/publications/COMING_IX_org_res.pdf)

Spatially Resolved Properties of Galaxies with a Kinematically Distinct Core

Authors: Kiyoaki Christopher Omori, Tsutomu T. Takeuchi

Abstract: Aims. Interacting galaxies show unique irregularities in their kinematic structure. By investigating the spatially resolved kinematics and stellar population properties of galaxies that show irregularities, we can paint a detailed picture of the formation and evolutionary processes that took place during its lifetimes. Methods. In this work, we focus on galaxies with a specific kinematic irregular… ▽ More Aims. Interacting galaxies show unique irregularities in their kinematic structure. By investigating the spatially resolved kinematics and stellar population properties of galaxies that show irregularities, we can paint a detailed picture of the formation and evolutionary processes that took place during its lifetimes. Methods. In this work, we focus on galaxies with a specific kinematic irregularity, a kinematically distinct stellar core (KDC), in particular, counter-rotating galaxies where the core and main body of the galaxy are rotating in opposite directions. We visually identify eleven MaNGA galaxies with a KDC from their stellar kinematics, and investigate their spatially resolved stellar and gaseous kinematic properties, namely the two-dimensional stellar and gaseous velocity and velocity dispersion (σ) maps. Additionally, we examine the stellar population properties, as well as spatially resolved recent star formation histories using the Dn4000 and Hδ gradients. Results. The galaxies display multiple off-centred, symmetrical peaks in the stellar σ maps. The gaseous velocity and σ maps display regular properties. The stellar population properties and their respective gradients show differing properties depending on the results of the spatially resolved emission line diagnostics of the galaxies, with some galaxies showing inside-out quenching but others not. The star formation histories also largely differ based on the spatially resolved emission line diagnostics, but most galaxies show indications of recent star formation either in their outskirts or core. Conclusions. We find a distinct difference in kinematic and stellar population properties in galaxies with a counter-rotating stellar core, depending on its classification using spatially resolved emission line diagnostics. △ Less

Submitted 14 December, 2020; originally announced December 2020.

Comments: 21 pages, 9 figures, accepted in A&A

Journal ref: A&A 647, A181 (2021)

An Iterative Reconstruction Algorithm for Faraday Tomography

Authors: Suchetha Cooray, Tsutomu T. Takeuchi, Takuya Akahori, Yoshimitsu Miyashita, Shinsuke Ideguchi, Keitaro Takahashi, Kiyotomo Ichiki

Abstract: Faraday tomography offers crucial information on the magnetized astronomical objects, such as quasars, galaxies, or galaxy clusters, by observing its magnetoionic media. The observed linear polarization spectrum is inverse Fourier transformed to obtain the Faraday dispersion function (FDF), providing us a tomographic distribution of the magnetoionic media along the line of sight. However, this tra… ▽ More Faraday tomography offers crucial information on the magnetized astronomical objects, such as quasars, galaxies, or galaxy clusters, by observing its magnetoionic media. The observed linear polarization spectrum is inverse Fourier transformed to obtain the Faraday dispersion function (FDF), providing us a tomographic distribution of the magnetoionic media along the line of sight. However, this transform gives a poor reconstruction of the FDF because of the instrument's limited wavelength coverage. The current Faraday tomography techniques' inability to reliably solve the above inverse problem has noticeably plagued cosmic magnetism studies. We propose a new algorithm inspired by the well-studied area of signal restoration, called the Constraining and Restoring iterative Algorithm for Faraday Tomography (CRAFT). This iterative model-independent algorithm is computationally inexpensive and only requires weak physically-motivated assumptions to produce high fidelity FDF reconstructions. We demonstrate an application for a realistic synthetic model FDF of the Milky Way, where CRAFT shows greater potential over other popular model-independent techniques. The dependence of observational frequency coverage on the various techniques' reconstruction performance is also demonstrated for a simpler FDF. CRAFT exhibits improvements even over model-dependent techniques (i.e., QU-fitting) by capturing complex multi-scale features of the FDF amplitude and polarization angle variations within a source. The proposed approach will be of utmost importance for future cosmic magnetism studies, especially with broadband polarization data from the Square Kilometre Array and its precursors. We make the CRAFT code publicly available. △ Less

Submitted 21 November, 2020; originally announced November 2020.

Comments: Accepted for publication in MNRAS. 13 pages and 12 figures

CO Multi-line Imaging of Nearby Galaxies (COMING). X. Physical conditions of molecular gas and the local SFR-Mass relation

Authors: Kana Morokuma-Matsui, Kazuo Sorai, Yuya Sato, Nario Kuno, Tsutomu T. Takeuchi, Dragan Salak, Yusuke Miyamoto, Yoshiyuki Yajima, Kazuyuki Muraoka, Hiroyuki Kaneko

Abstract: We investigate the molecular gas properties of galaxies across the main sequence of star-forming (SF) galaxies in the local Universe using $^{12}$CO($J=1-0$) (hereafter $^{12}$CO) and $^{13}$CO($J=1-0$) ($^{13}$CO) mapping data of 147 nearby galaxies obtained in the COMING project, a legacy project of the Nobeyama Radio Observatory. In order to improve signal-to-noise ratios of both lines, we stac… ▽ More We investigate the molecular gas properties of galaxies across the main sequence of star-forming (SF) galaxies in the local Universe using $^{12}$CO($J=1-0$) (hereafter $^{12}$CO) and $^{13}$CO($J=1-0$) ($^{13}$CO) mapping data of 147 nearby galaxies obtained in the COMING project, a legacy project of the Nobeyama Radio Observatory. In order to improve signal-to-noise ratios of both lines, we stack all the pixels where $^{12}$CO emission is detected after aligning the line center expected from the first-moment map of $^{12}$CO. As a result, $^{13}$CO emission is successfully detected in 80 galaxies with a signal-to-noise ratio larger than three. The error-weighted mean of integrated-intensity ratio of $^{12}$CO to $^{13}$CO lines ($R_{1213}$) of the 80 galaxies is $10.9$ with a standard deviation of $7.0$. We find that (1) $R_{1213}$ positively correlates to specific star-formation rate (sSFR) with a correlation coefficient of $0.46$, and (2) both flux ratio of IRAS 60~$μ$m to 100~$μ$m ($f_{60}/f_{100}$) and inclination-corrected linewidth of $^{12}$CO stacked spectra ($σ_{{\rm ^{12}CO},i}$) also correlate with sSFR for galaxies with the $R_{1213}$ measurement. Our results support the scenario where $R_{1213}$ variation is mainly caused by the changes in molecular-gas properties such as temperature and turbulence. The consequent variation of CO-to-H$_2$ conversion factor across the SF main sequence is not large enough to completely extinguish the known correlations between sSFR and $M_{\rm mol}/M_{\rm star}$ ($μ_{\rm mol}$) or star-formation efficiency (SFE) reported in previous studies, while this variation would strengthen (weaken) the sSFR-SFE (sSFR-$μ_{\rm mol}$) correlation. △ Less

Submitted 3 August, 2020; originally announced August 2020.

Comments: 16 pages, 7 figures, 1 table, accepted for publication in PASJ

Constructing a multivariate distribution function with a vine copula: toward multivariate luminosity and mass functions

Authors: Tsutomu T. Takeuchi, Kai T. Kono

Abstract: The need for a method to construct multidimensional distribution function is increasing recently, in the era of huge multiwavelength surveys. We have proposed a systematic method to build a bivariate luminosity or mass function of galaxies by using a copula. It allows us to construct a distribution function when only its marginal distributions are known, and we have to estimate the dependence stru… ▽ More The need for a method to construct multidimensional distribution function is increasing recently, in the era of huge multiwavelength surveys. We have proposed a systematic method to build a bivariate luminosity or mass function of galaxies by using a copula. It allows us to construct a distribution function when only its marginal distributions are known, and we have to estimate the dependence structure from data. A typical example is the situation that we have univariate luminosity functions at some wavelengths for a survey, but the joint distribution is unknown. Main limitation of the copula method is that it is not easy to extend a joint function to higher dimensions ($d > 2$), except some special cases like multidimensional Gaussian. Even if we find such a multivariate analytic function in some fortunate case, it would often be inflexible and impractical. In this work, we show a systematic method to extend the copula method to unlimitedly higher dimensions by a vine copula. This is based on the pair-copula decomposition of a general multivariate distribution. We show how the vine copula construction is flexible and extendable. We also present an example of the construction of an stellar mass--atomic gas--molecular gas 3-dimensional mass function. We demonstrate the maximum likelihood estimation of the best functional form for this function, as well as a proper model selection via vine copula. △ Less

Submitted 18 August, 2020; v1 submitted 10 June, 2020; originally announced June 2020.

Comments: 15 pages, 9 figures. Submitted to MNRAS on 7 May 2020, accepted for publication on 17 August 2020

A Study on the Baryon Acoustic Oscillation with Topological Data Analysis

Authors: Kai T. Kono, Tsutomu T. Takeuchi, Suchetha Cooray, Atsushi J. Nishizawa, Koya Murakami

Abstract: (abridged) The scale of the acoustic oscillation of baryons at the baryon-photon decoupling is imprinted on the spatial distribution of galaxies in the Universe, known as the baryon acoustic oscillation (BAO). The correlation functions and power spectrum are used as a central tool for the studies on the BAO analysis. In this work, we analyzed the spatial distribution of galaxies with a method from… ▽ More (abridged) The scale of the acoustic oscillation of baryons at the baryon-photon decoupling is imprinted on the spatial distribution of galaxies in the Universe, known as the baryon acoustic oscillation (BAO). The correlation functions and power spectrum are used as a central tool for the studies on the BAO analysis. In this work, we analyzed the spatial distribution of galaxies with a method from the topological data analysis (TDA), in order to detect and examine the BAO signal in the galaxy distribution. The TDA provides a method to treat various types of "holes" in point set data, by constructing the persistent homology (PH) group from the geometric structure of data points and handling the topological information of the dataset. We can obtain the information on the size, position, and statistical significance of the holes in the data. A particularly strong point of the persistent homology is that it can classify the holes by their spatial dimension, i.e., a 0-dim separation, 1-dim loop, 2-dim shell, etc. We first analyzed the simulation datasets with and without the baryon physics to examine the performance of the PH method. We found that the PH is indeed able to detect the BAO signal: simulation data with baryon physics present a prominent signal from the BAO, while data without baryon physics does not show this signal. Then, we applied the PH to a quasar sample at $z <1.0$ from extended Baryon Oscillation Spectroscopic Survey in Sloan Digital Sky Survey Data Release 14. We discovered a characteristic hole (a hollow shell) at a scaler $\sim150\ [{\rm Mpc}]$. This exactly corresponds to the BAO signature imprinted in the galaxy/quasar distribution. We performed this analysis on a small subsample of 2000 quasars. This clearly demonstrates that the PH analysis is very efficient in finding this type of topological structures even if the sampling is very sparse. △ Less

Submitted 4 June, 2020; originally announced June 2020.

Comments: 29 pages, 13 figures

A Method for Unmasking Incomplete Astronomical Signals: Application to CO Multi-line Imaging of Nearby Galaxies Project

Authors: Suchetha Cooray, Tsutomu T. Takeuchi, Moe Yoda, Kazuo Sorai

Abstract: Photometric surveys have provided incredible amounts of astronomical information in the form of images. However, astronomical images often contain artifacts that can critically hinder scientific analysis by misrepresenting intensities or contaminating catalogs as artificial objects. These affected pixels need to be masked and dealt with in any data reduction pipeline. In this paper, we present a f… ▽ More Photometric surveys have provided incredible amounts of astronomical information in the form of images. However, astronomical images often contain artifacts that can critically hinder scientific analysis by misrepresenting intensities or contaminating catalogs as artificial objects. These affected pixels need to be masked and dealt with in any data reduction pipeline. In this paper, we present a flexible, iterative algorithm to recover (unmask) astronomical images where some pixels are lacking. We demonstrate the application of the method on some intensity calibration source images in CO Multi-line Imaging of Nearby Galaxies (COMING) Project conducted using the 45m telescope at Nobeyama Radio Observatory (NRO). The proposed algorithm restored artifacts due to a detector error in the intensity calibration source images. The restored images were used to calibrate 11 out of 147 observed galaxy maps in the survey. The tests show that the algorithm can restore measured intensities at sub 1% error even for noisy images (SNR = 2.4), despite lacking a significant part of the image. We present the formulation of the reconstruction algorithm, discuss its possibilities and limitations for extensions to other astronomical signals and the results of the COMING application. △ Less

Submitted 15 April, 2020; originally announced April 2020.

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

Observational and theoretical constraints on the formation and early evolution of the first dust grains in galaxies at 5 < z < 10

Authors: D. Burgarella, A. Nanni, H. Hirashita, P. Theule, A. K. Inoue, T. T. Takeuchi

Abstract: The first generation of stars were born a few hundred million years after the big bang. These stars synthesized elements heavier than H and He, that are later expelled into the interstellar medium, initiating the rise of metals. Within this enriched medium, the first dust grains formed. This event is cosmological crucial for molecule formation as dust plays a major role by cooling low-metallicity… ▽ More The first generation of stars were born a few hundred million years after the big bang. These stars synthesized elements heavier than H and He, that are later expelled into the interstellar medium, initiating the rise of metals. Within this enriched medium, the first dust grains formed. This event is cosmological crucial for molecule formation as dust plays a major role by cooling low-metallicity star-forming clouds which can fragment to create lower mass stars. Collecting information on these first dust grains is difficult because of the negative alliance of large distances and low dust masses. We combine the observational information from galaxies at redshifts 5 < z < 10 to constrain their dust emission and theoretically understand the first evolutionary phases of the dust cycle. Spectral energy distributions (SEDs) are fitted with CIGALE and the physical parameters and their evolution are modelled. From this SED fitting, we build a dust emission template for this population of galaxies in the epoch of reionization. Our new models explain why some early galaxies are observed and others are not. We follow in time the formation of the first grains by supernovae later destroyed by other supernova blasts and expelled in the circumgalactic and intergalactic media. We have found evidence for the first dust grains formed in the universe. But, above all, this letter underlines the need to collect more data and to develop new facilities to further constrain the dust cycle in galaxies in the epoch of reionization. △ Less

Submitted 11 February, 2020; v1 submitted 5 February, 2020; originally announced February 2020.

Comments: Paper accepted for publication in section 4. Extragalactic astronomy of Astronomy and Astrophysics (Feb. 2020). Correction typo in caption of Fig. 1

Journal ref: A&A 637, A32 (2020)

The fate of the interstellar medium in early-type galaxies. I. First direct measurement of the timescale of dust removal

Authors: Michał J. Michałowski, J. Hjorth, C. Gall, D. T. Frayer, A. -L. Tsai, H. Hirashita, K. Rowlands, T. T. Takeuchi, A. Leśniewska, D. Behrendt, N. Bourne, D. H. Hughes, E. Spring, J. Zavala, P. Bartczak

Abstract: An important aspect of quenching star formation is the removal of the cold interstellar medium (ISM; non-ionised gas and dust) from a galaxy. In addition, dust grains can be destroyed in a hot or turbulent medium. The adopted timescale of dust removal usually relies on uncertain theoretical estimates. It is tricky to track the dust removal, because usually dust is constantly replenished by consecu… ▽ More An important aspect of quenching star formation is the removal of the cold interstellar medium (ISM; non-ionised gas and dust) from a galaxy. In addition, dust grains can be destroyed in a hot or turbulent medium. The adopted timescale of dust removal usually relies on uncertain theoretical estimates. It is tricky to track the dust removal, because usually dust is constantly replenished by consecutive generations of stars. Our objective is to measure observationally the timescale of dust removal. We here explore an approach to select galaxies which do have detectable amounts of dust and cold ISM but exhibit a low current dust production rate. Any decrease of the dust and gas content as a function of the age of such galaxies therefore must be attributed to processes governing the ISM removal. We used a sample of galaxies detected by Herschel in the far-infrared with visually assigned early-type morphology or spirals with red colours. We also obtained JCMT/SCUBA-2 observations for five of them. We discovered an exponential decline of the dust-to-stellar mass ratio with age, which we interpret as an evolutionary trend of dust removal from these galaxies. For the first time we directly measure the dust removal timescale in such galaxies to be tau=(2.5+-0.4) Gyr (the corresponding half-life time is (1.75+-0.25) Gyr). This quantity may be used in models in which it must be assumed a priori and cannot be derived. Any process which removes dust in these galaxies, such as dust grain destruction, cannot happen on shorter timescales. The timescale is comparable to the quenching timescales found in simulations for galaxies with similar stellar masses. The dust is likely of internal, not external origin. It was either formed in the past directly by supernovae, or from seeds produced by SNe and with grain growth in the ISM contributing substantially to the dust mass accumulation. △ Less

Submitted 14 October, 2019; originally announced October 2019.

Comments: Astronomy & Astrophysics, accepted; 13 pages, 9 figures, 1 table

Journal ref: A&A 632, A43 (2019)

CO Multi-line Imaging of Nearby Galaxies (COMING) IV. Overview of the Project

Authors: Kazuo Sorai, Nario Kuno, Kazuyuki Muraoka, Yusuke Miyamoto, Hiroyuki Kaneko, Hiroyuki Nakanishi, Naomasa Nakai, Kazuki Yanagitani, Takahiro Tanaka, Yuya Sato, Dragan Salak, Michiko Umei, Kana Morokuma-Matsui, Naoko Matsumoto, Saeko Ueno, Hsi-An Pan, Yuto Noma, Tsutomu T. Takeuchi, Moe Yoda, Mayu Kuroda, Atsushi Yasuda, Yoshiyuki Yajima, Nagisa Oi, Shugo Shibata, Masumichi Seta , et al. (14 additional authors not shown)

Abstract: Observations of the molecular gas in galaxies are vital to understanding the evolution and star-forming histories of galaxies. However, galaxies with molecular gas maps of their whole discs having sufficient resolution to distinguish galactic structures are severely lacking. Millimeter wavelength studies at a high angular resolution across multiple lines and transitions are particularly needed, se… ▽ More Observations of the molecular gas in galaxies are vital to understanding the evolution and star-forming histories of galaxies. However, galaxies with molecular gas maps of their whole discs having sufficient resolution to distinguish galactic structures are severely lacking. Millimeter wavelength studies at a high angular resolution across multiple lines and transitions are particularly needed, severely limiting our ability to infer the universal properties of molecular gas in galaxies. Hence, we conducted a legacy project with the 45 m telescope of the Nobeyama Radio Observatory, called the CO Multi-line Imaging of Nearby Galaxies (COMING), which simultaneously observed 147 galaxies with high far-infrared flux in $^{12}$CO, $^{13}$CO, and C$^{18}$O $J=1-0$ lines. The total molecular gas mass was derived using the standard CO-to-H$_2$ conversion factor and found to be positively correlated with the total stellar mass derived from the WISE $3.4 μ$m band data. The fraction of the total molecular gas mass to the total stellar mass in galaxies does not depend on their Hubble types nor the existence of a galactic bar, although when galaxies in individual morphological types are investigated separately, the fraction seems to decrease with the total stellar mass in early-type galaxies and vice versa in late-type galaxies. No differences in the distribution of the total molecular gas mass, stellar mass, and the total molecular gas to stellar mass ratio was observed between barred and non-barred galaxies, which is likely the result of our sample selection criteria, in that we prioritized observing FIR bright (and thus molecular gas-rich) galaxies. △ Less

Submitted 9 October, 2019; originally announced October 2019.

Comments: Accepted for publication in PASJ; 47 pages, 5 tables, 29 figures. On-line supplementary images are available at this URL (https://astro3.sci.hokudai.ac.jp/~radio/coming/publications/). CO data is available at the Japanese Virtual Observatory (JVO) website (https://jvo.nao.ac.jp/portal/nobeyama/coming.do) and the project website (https://astro3.sci.hokudai.ac.jp/~radio/coming/data/)

ADF22: Blind detections of [CII] line emitters shown to be spurious

Authors: Natsuki H. Hayatsu, Rob J. Ivison, Paola Andreani, Hideki Umehata, Yuichi Matsuda, Naoki Yoshida, Kotaro Kohno, Bunyo Hatsukade, Akio K. Inoue, Yoichi Tamura, Tutomu T. Takeuchi, Seiji Fujimoto, Minju M. Lee, Tohru Nagao, Yiping Ao

Abstract: We report Atacama Large Millimetre/submillimeter Array (ALMA) Cycle-5 follow-up observations of two candidate [CII] emitters at z ~ 6 in the ALMA deep field in SSA22 (ADF22). The candidates were detected blindly in a Cycle-2 ALMA survey covering ~ 5 square arcmins, with a single tuning, along with two CO lines associated with galaxies at lower redshifts. Various tests suggested at least one of the… ▽ More We report Atacama Large Millimetre/submillimeter Array (ALMA) Cycle-5 follow-up observations of two candidate [CII] emitters at z ~ 6 in the ALMA deep field in SSA22 (ADF22). The candidates were detected blindly in a Cycle-2 ALMA survey covering ~ 5 square arcmins, with a single tuning, along with two CO lines associated with galaxies at lower redshifts. Various tests suggested at least one of the two > 6-sigma [CII] candidates should be robust (Hayatsu et al. 2017). Nevertheless, our new, deeper observations recover neither candidate, demonstrating a higher contamination rate than expected. The cause of the spurious detections is under investigation but at present it remains unclarified. △ Less

Submitted 12 July, 2019; originally announced July 2019.

Comments: Submitted to AAS research note. One figure, three pages. Comments, suggestions, questions are welcome. The details will be discussed in Hayatsu et al. in prep

Complex distribution and velocity field of molecular gas in NGC 1316 as revealed by Morita Array of ALMA

Authors: Kana Morokuma-Matsui, Paolo Serra, Filippo M. Maccagni, Bi-Qing For, Jing Wang, Kenji Bekki, Tomoki Morokuma, Fumi Egusa, Daniel Espada, Rie Miura, E., Kouichiro Nakanishi, Bärbel S. Koribalski, Tsutomu T. Takeuchi

Abstract: We present the results of $^{12}$CO($J$=1-0) mosaicing observations of the cD galaxy NGC 1316 at kpc-resolution performed with the Morita Array of the Atacama Large Millimeter/submillimeter Array (ALMA). We reveal the detailed distribution of the molecular gas in the central region for the first time: a shell structure in the northwest, a barely resolved blob in the southeast of the center and som… ▽ More We present the results of $^{12}$CO($J$=1-0) mosaicing observations of the cD galaxy NGC 1316 at kpc-resolution performed with the Morita Array of the Atacama Large Millimeter/submillimeter Array (ALMA). We reveal the detailed distribution of the molecular gas in the central region for the first time: a shell structure in the northwest, a barely resolved blob in the southeast of the center and some clumps between them. The total molecular gas mass obtained with a standard Milky-Way CO-to-H$_2$ conversion factor is $(5.62 \pm 0.53)\times10^8$ M$_\odot$, which is consistent with previous studies. The disturbed velocity field of the molecular gas suggests that the molecular gas is injected very recently ($<1$ Gyr) if it has an external origin and is in the process of settling into a rotating disk. Assuming that a low-mass gas-rich galaxy has accreted, the gas-to-dust ratio and H$_2$-to-HI ratio are unusually low ($\sim 28$) and high ($\sim 5.6$), respectively. To explain these ratios, additional processes should be taken into accounts such as an effective dust formation and conversion from atomic to molecular gas during the interaction. We also discuss the interaction between the nuclear jet and the molecular gas. △ Less

Submitted 27 May, 2019; originally announced May 2019.

Comments: 16 pages, 8 figures, 3 tables, accepted for publication in PASJ

AGN selection in the AKARI NEP deep field with the fuzzy SVM algorithm

Authors: Artem Poliszczuk, Aleksandra Solarz, Agnieszka Pollo, Maciej Bilicki, Tsutomu T. Takeuchi, Hideo Matsuhara, Tomotsugu Goto, Toshinobu Takagi, Takehiko Wada, Yoichi Ohyama, Hitoshi Hanami, Takamitsu Miyaji, Nagisa Oi, Matthew Malkan, Kazumi Murata, Helen Kim, Jorge Díaz Tello

Abstract: The aim of this work is to create a new catalog of reliable AGN candidates selected from the AKARI NEP-Deep field. Selection of the AGN candidates was done by applying a fuzzy SVM algorithm, which allows to incorporate measurement uncertainties into the classification process. The training dataset was based on the spectroscopic data available for selected objects in the NEP-Deep and NEP-Wide field… ▽ More The aim of this work is to create a new catalog of reliable AGN candidates selected from the AKARI NEP-Deep field. Selection of the AGN candidates was done by applying a fuzzy SVM algorithm, which allows to incorporate measurement uncertainties into the classification process. The training dataset was based on the spectroscopic data available for selected objects in the NEP-Deep and NEP-Wide fields. The generalization sample was based on the AKARI NEP-Deep field data including objects without optical counterparts and making use of the infrared information only. A high quality catalog of previously unclassified 275 AGN candidates was prepared. △ Less

Submitted 18 March, 2019; v1 submitted 11 February, 2019; originally announced February 2019.

Comments: 17 pages, accepted to PASJ AKARI special issue

Journal ref: Publications of the Astronomical Society of Japan, Volume 71, Issue 3, June 2019, 65

CO Multi-line Imaging of Nearby Galaxies (COMING). III. Dynamical effect on molecular gas density and star formation in the barred spiral galaxy NGC 4303

Authors: Yoshiyuki Yajima, Kazuo Sorai, Nario Kuno, Kazuyuki Muraoka, Yusuke Miyamoto, Hiroyuki Kaneko, Hiroyuki Nakanishi, Naomasa Nakai, Takahiro Tanaka, Yuya Sato, Dragan Salak, Kana Morokuma-Matsui, Naoko Matsumoto, His-An Pan, Yuto Noma, Tsutomu T. Takeuchi, Moe Yoda, Mayu Kuroda, Atsushi Yasuda, Nagisa Oi, Shugo Shibata, Masumichi Seta, Yoshimasa Watanabe, Shoichiro Kita, Ryusei Komatsuzaki , et al. (2 additional authors not shown)

Abstract: We present the results of $^{12}$CO($J$=1-0) and $^{13}$CO($J$=1-0) simultaneous mappings toward the nearby barred spiral galaxy NGC 4303 as a part of the CO Multi-line Imaging of Nearby Galaxies (COMING) project. Barred spiral galaxies often show lower star-formation efficiency (SFE) in their bar region compared to the spiral arms. In this paper, we examine the relation between the SFEs and the v… ▽ More We present the results of $^{12}$CO($J$=1-0) and $^{13}$CO($J$=1-0) simultaneous mappings toward the nearby barred spiral galaxy NGC 4303 as a part of the CO Multi-line Imaging of Nearby Galaxies (COMING) project. Barred spiral galaxies often show lower star-formation efficiency (SFE) in their bar region compared to the spiral arms. In this paper, we examine the relation between the SFEs and the volume densities of molecular gas $n(\rm{H}_2)$ in the eight different regions within the galactic disk with CO data combined with archival far-ultraviolet and 24 $μ$m data. We confirmed that SFE in the bar region is lower by 39% than that in the spiral arms. Moreover, velocity-alignment stacking analysis was performed for the spectra in the individual regions. The integrated intensity ratios of $^{12}$CO to $^{13}$CO ($R_{12/13}$) range from 10 to 17 as the results of stacking. Fixing a kinetic temperature of molecular gas, $n(\rm{H}_2)$ was derived from $R_{12/13}$ via non-local thermodynamic equilibrium (non-LTE) analysis. The density $n(\rm{H}_2)$ in the bar is lower by 31-37% than that in the arms and there is a rather tight positive correlation between SFEs and $n(\rm{H}_2)$, with a correlation coefficient of $\sim 0.8$. Furthermore, we found a dependence of $n(\rm{H}_2)$ on the velocity dispersion of inter-molecular clouds ($ΔV/ \sin i$). Specifically, $n(\rm{H}_2)$ increases as $ΔV/ \sin i$ increases when $ΔV/ \sin i < 100$ km s$^{-1}$. On the other hand, $n(\rm{H}_2)$ decreases as $ΔV/ \sin i$ increases when $ΔV/ \sin i > 100$ km s$^{-1}$. These relations indicate that the variations of SFE could be caused by the volume densities of molecular gas, and the volume densities could be governed by the dynamical influence such as cloud-cloud collisions, shear and enhanced inner-cloud turbulence. △ Less

Submitted 12 February, 2019; originally announced February 2019.

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

Journal ref: 2019, PASJ, 71, S13

CO Multi-line Imaging of Nearby Galaxies (COMING): VI. Radial variations in star formation efficiency

Authors: Kazuyuki Muraoka, Kazuo Sorai, Yusuke Miyamoto, Moe Yoda, Kana Morokuma-matsui, Masato I. N. Kobayashi, Mayu Kuroda, Hiroyuki Kaneko, Nario Kuno, Tsutomu T. Takeuchi, Hiroyuki Nakanishi, Yoshimasa Watanabe, Takahiro Tanaka, Atsushi Yasuda, Yoshiyuki Yajima, Shugo Shibata, Dragan Salak, Daniel Espada, Naoko Matsumoto, Yuto Noma, Shoichiro Kita, Ryusei Komatsuzaki, Ayumi Kajikawa, Yu Yashima, Hsi-An Pan , et al. (3 additional authors not shown)

Abstract: We examined radial variations in molecular-gas based star formation efficiency (SFE), which is defined as star formation rate per unit molecular gas mass, for 80 galaxies selected from the CO Multi-line Imaging of Nearby Galaxies project (Sorai et al. 2019). The radial variations in SFE for individual galaxies are typically a factor of 2 -- 3, which suggests that SFE is nearly constant along galac… ▽ More We examined radial variations in molecular-gas based star formation efficiency (SFE), which is defined as star formation rate per unit molecular gas mass, for 80 galaxies selected from the CO Multi-line Imaging of Nearby Galaxies project (Sorai et al. 2019). The radial variations in SFE for individual galaxies are typically a factor of 2 -- 3, which suggests that SFE is nearly constant along galactocentric radius. We found the averaged SFE in 80 galaxies of $(1.69 \pm 1.1) \times 10^{-9}$ yr$^{-1}$, which is consistent with Leroy et al. 2008 if we consider the contribution of helium to the molecular gas mass evaluation and the difference in the assumed initial mass function between two studies. We compared SFE among different morphological (i.e., SA, SAB, and SB) types, and found that SFE within the inner radii ($r/r_{25} < 0.3$, where $r_{25}$ is $B$-band isophotal radius at 25 mag arcsec$^{-2}$) of SB galaxies is slightly higher than that of SA and SAB galaxies. This trend can be partly explained by the dependence of SFE on global stellar mass, which probably relates to the CO-to-H$_2$ conversion factor through the metallicity. For two representative SB galaxies in our sample, NGC 3367 and NGC 7479, the ellipse of $r/r_{25}$ = 0.3 seems to cover not only the central region but also the inner part of the disk, mainly the bar. These two galaxies show higher SFE in the bar than in spiral arms. However, we found an opposite trend in NGC 4303; SFE is lower in the bar than in spiral arms, which is consistent with earlier studies (e.g., Momose et al. 2010). These results suggest diversity of star formation activities in the bar. △ Less

Submitted 30 January, 2019; originally announced January 2019.

Comments: 28 pages, 12 figures, accepted for publication in PASJ

Radio-Infrared Correlation for Local Dusty Galaxies and Dusty AGNs from the AKARI All Sky Survey

Authors: A. Solarz, A. Pollo, M. Bilicki, A. Pępiak, T. T. Takeuchi, P. Piątek

Abstract: We use the new release of the AKARI Far-Infrared all sky Survey matched with the NVSS radio database to investigate the local ($z<0.25$) far infrared-radio correlation (FIRC) of different types of extragalactic sources. To obtain the redshift information for the AKARI FIS sources we crossmatch the catalogue with the SDSS DR8. This also allows us to use emission line properties to divide sources in… ▽ More We use the new release of the AKARI Far-Infrared all sky Survey matched with the NVSS radio database to investigate the local ($z<0.25$) far infrared-radio correlation (FIRC) of different types of extragalactic sources. To obtain the redshift information for the AKARI FIS sources we crossmatch the catalogue with the SDSS DR8. This also allows us to use emission line properties to divide sources into four categories: i) star-forming galaxies (SFGs), ii) composite galaxies (displaying both star-formation and active nucleus components), iii) Seyfert galaxies, and iv) low-ionization nuclear emission-line region (LINER) galaxies. We find that the Seyfert galaxies have the lowest FIR/radio flux ratios and display excess radio emission when compared to the SFGs. We conclude that FIRC can be used to separate SFGs and AGNs only for the most radio-loud objects. △ Less

Submitted 29 January, 2019; originally announced January 2019.

Comments: 9 pages, accepted to PASJ

CO Multi-line Imaging of Nearby Galaxies (COMING). VII. Fourier decomposition of molecular gas velocity fields and bar pattern speed

Authors: Dragan Salak, Yuto Noma, Kazuo Sorai, Yusuke Miyamoto, Nario Kuno, Alex R. Pettitt, Hiroyuki Kaneko, Takahiro Tanaka, Atsushi Yasuda, Shoichiro Kita, Yoshiyuki Yajima, Shugo Shibata, Naomasa Nakai, Masumichi Seta, Kazuyuki Muraoka, Mayu Kuroda, Hiroyuki Nakanishi, Tsutomu T. Takeuchi, Moe Yoda, Kana Morokuma-Matsui, Yoshimasa Watanabe, Naoko Matsumoto, Nagisa Oi, Hsi-An Pan, Ayumi Kajikawa , et al. (2 additional authors not shown)

Abstract: The $^{12}$CO $(J=1\rightarrow0)$ velocity fields of a sample of 20 nearby spiral galaxies, selected from the CO Multi-line Imaging of Nearby Galaxies (COMING) legacy project of Nobeyama Radio Observatory, have been analyzed by Fourier decomposition to determine their basic kinematic properties, such as circular and noncircular velocities. On average, the investigated barred (SAB and SB) galaxies… ▽ More The $^{12}$CO $(J=1\rightarrow0)$ velocity fields of a sample of 20 nearby spiral galaxies, selected from the CO Multi-line Imaging of Nearby Galaxies (COMING) legacy project of Nobeyama Radio Observatory, have been analyzed by Fourier decomposition to determine their basic kinematic properties, such as circular and noncircular velocities. On average, the investigated barred (SAB and SB) galaxies exhibit a ratio of noncircular to circular velocities of molecular gas larger by a factor of 1.5-2 than non-barred (SA) spiral galaxies at radii within the bar semimajor axis $a_\mathrm{b}$ at 1 kpc resolution, with a maximum at a radius of $R/a_\mathrm{b}\sim0.3$. Residual velocity field images, created by subtracting model velocity fields from the data, reveal that this trend is caused by kpc-scale streaming motions of molecular gas in the bar region. Applying a new method based on radial velocity reversal, we estimated the corotation radius $R_\mathrm{CR}$ and bar pattern speed $Ω_\mathrm{b}$ in seven SAB and SB systems. The ratio of the corotation to bar radius is found to be in a range of $\mathcal{R}\equiv R_\mathrm{CR}/a_\mathrm{b}\sim0.8\mathrm{-}1.6$, suggesting that intermediate (SBb-SBc), luminous barred spiral galaxies host fast and slow rotator bars. Tentative negative correlations are found for $Ω_\mathrm{b}$ vs. $a_\mathrm{b}$ and $Ω_\mathrm{b}$ vs. total stellar mass $M_\ast$, indicating that bars in massive disks are larger and rotate slower, possibly a consequence of angular momentum transfer. The kinematic properties of SAB and SB galaxies, derived from Fourier decomposition, are compared with recent numerical simulations that incorporate various rotation curve models and galaxy interactions. △ Less

Submitted 3 January, 2019; originally announced January 2019.

Comments: Accepted for publication in PASJ; 39 pages, 15 figures

AKARI mid-infrared slitless spectroscopic survey of star-forming galaxies at $z\lesssim 0.5$

Authors: Y. Ohyama, T. Wada, H. Matsuhara, T. Takagi, M. Malkan, T. Goto, E. Egami, H. M. Lee, M. Im, J. -H. Kim, C. Pearson, H. Inami, S. Oyabu, F. Usui, D. Burgarella, F. Mazyed, M. Imanishi, W. -S. Jeong, T. Miyaji, J. Díaz Tello, T. Nakagawa, S. Serjeant, T. T. Takeuchi, Y. Toba, G. J. White , et al. (2 additional authors not shown)

Abstract: Deep MIR surveys have revealed numerous strongly star-forming galaxies at redshift z~<2. Their MIR fluxes are produced by a combination of continuum and PAH emission features. The PAH features can dominate the total MIR flux, but are difficult to measure without spectroscopy. We aim to study star-forming galaxies by using a blind spectroscopic survey at MIR wavelengths to understand evolution of t… ▽ More Deep MIR surveys have revealed numerous strongly star-forming galaxies at redshift z~<2. Their MIR fluxes are produced by a combination of continuum and PAH emission features. The PAH features can dominate the total MIR flux, but are difficult to measure without spectroscopy. We aim to study star-forming galaxies by using a blind spectroscopic survey at MIR wavelengths to understand evolution of their star formation rate (SFR) and specific SFR up to z~=0.5, by paying particular attention to their PAH properties. We conducted a low-resolution (R~=50) slitless spectroscopic survey at 5-13um of 9um flux-selected sources (>0.3 mJy) around the North Ecliptic Pole with the Infrared Camera (IRC) onboard AKARI. After removing 11 AGN candidates by using the IRC photometry, we identified 48 PAH galaxies with PAH 6.2, 7.7, and 8.6um features at z<0.5. The rest-frame optical-MIR spectral energy distributions (SEDs) based on CFHT and AKARI/IRC imaging covering 0.37-18um were produced, and analysed in conjunction with the PAH spectroscopy. We defined the PAH enhancement by using the luminosity ratio of the 7.7um PAH feature over the 3.5um stellar component of the SEDs. The rest-frame SEDs of all PAH galaxies have a universal shape with stellar and 7.7um bumps, except that the PAH enhancement significantly varies as a function of the PAH luminosities. We identified a PAH-enhanced population at z~>0.35, whose SEDs and luminosities are typical of luminous infrared galaxies. They show particularly larger PAH enhancement at high luminosity, implying that they are vigorous star-forming galaxies with elevated specific SFR. Our composite starburst model that combines a very young and optically very thick starburst with a very old population can successfully reproduce most of their SED characteristics, although we could not confirm this optically thick component from our spectral analysis. △ Less

Submitted 14 August, 2018; originally announced August 2018.

Comments: 61 pages, 24 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 618, A101 (2018)

Detection of the Far-infrared [O III] and Dust Emission in a Galaxy at Redshift 8.312: Early Metal Enrichment in the Heart of the Reionization Era

Authors: Yoichi Tamura, Ken Mawatari, Takuya Hashimoto, Akio K. Inoue, Erik Zackrisson, Lise Christensen, Christian Binggeli, Yuichi Matsuda, Hiroshi Matsuo, Tsutomu T. Takeuchi, Ryosuke S. Asano, Kaho Sunaga, Ikkoh Shimizu, Takashi Okamoto, Naoki Yoshida, Minju Lee, Takatoshi Shibuya, Yoshiaki Taniguchi, Hideki Umehata, Bunyo Hatsukade, Kotaro Kohno, Kazuaki Ota

Abstract: We present the Atacama Large Millimeter/submillimeter Array (ALMA) detection of the [O III] 88 $μ$m line and rest-frame 90 $μ$m dust continuum emission in a Y-dropout Lyman break galaxy (LBG), MACS0416_Y1, lying behind the Frontier Field cluster MACS J0416.1-2403. This [O III] detection confirms the LBG with a spectroscopic redshift of $z = 8.3118 \pm 0.0003$, making this object one of the furthes… ▽ More We present the Atacama Large Millimeter/submillimeter Array (ALMA) detection of the [O III] 88 $μ$m line and rest-frame 90 $μ$m dust continuum emission in a Y-dropout Lyman break galaxy (LBG), MACS0416_Y1, lying behind the Frontier Field cluster MACS J0416.1-2403. This [O III] detection confirms the LBG with a spectroscopic redshift of $z = 8.3118 \pm 0.0003$, making this object one of the furthest galaxies ever identified spectroscopically. The observed 850 $μ$m flux density of $137 \pm 26$ $μ$Jy corresponds to a de-lensed total infrared (IR) luminosity of $L_{\rm IR} = (1.7 \pm 0.3) \times 10^{11} L_{\odot}$ if assuming a dust temperature of $T_{\rm dust} = 50$ K and an emissivity index of $β= 1.5$, yielding a large dust mass of $4 \times 10^6 M_{\odot}$. The ultraviolet-to-far IR spectral energy distribution modeling where the [O III] emissivity model is incorporated suggests the presence of a young ($τ_{\rm age} \approx 4$ Myr), star-forming (SFR $\approx 60 M_{\odot}$ yr$^{-1}$), moderately metal-polluted ($Z \approx 0.2 Z_{\odot}$) stellar component with a mass of $M_{\rm star} = 3 \times 10^8 M_{\odot}$. An analytic dust mass evolution model with a single episode of star-formation does not reproduce the metallicity and dust mass in $τ_{\rm age} \approx 4$ Myr, suggesting a pre-existing evolved stellar component with $M_{\rm star} \sim 3 \times 10^9 M_{\odot}$ and $τ_{\rm age} \sim 0.3$ Gyr as the origin of the dust mass. △ Less

Submitted 11 February, 2019; v1 submitted 11 June, 2018; originally announced June 2018.

Comments: Accepted by ApJ. 18 pages, 10 figures, 5 tables

The VIMOS Public Extragalactic Redshift Survey (VIPERS). Unsupervised classification with photometric redshifts: a method to accurately classify large galaxy samples without spectroscopic information

Authors: M. Siudek, K. Małek, A. Pollo, B. R. Granett, M. Scodeggio, T. Moutard, A. Iovino, L. Guzzo, B. Garilli, M. Bolzonella, S. de la Torre, U. Abbas, C. Adami, D. Bottini, A. Cappi, O. Cucciati, I. Davidzon, P. Franzetti, A. Fritz, J. Krywult, V. Le Brun, O. Le Fèvre, D. Maccagni, F. Marulli, M. Polletta , et al. (13 additional authors not shown)

Abstract: Techniques to classify galaxies solely based on photometry will be necessary for future large cosmology missions, such as Euclid or LSST. However, the precision of classification is always lower in photometric surveys and can be systematically biased with respect to classifications based upon spectroscopic data. We verified how precisely the detailed classification scheme introduced by Siudek et a… ▽ More Techniques to classify galaxies solely based on photometry will be necessary for future large cosmology missions, such as Euclid or LSST. However, the precision of classification is always lower in photometric surveys and can be systematically biased with respect to classifications based upon spectroscopic data. We verified how precisely the detailed classification scheme introduced by Siudek et al. (2018, hereafter: S1) for galaxies at z~0.7 could be reproduced if only photometric data are available. We applied the Fisher Expectation-Maximization (FEM) unsupervised clustering algorithm to 54,293 VIPERS galaxies working in a parameter space of reliable photometric redshifts and 12 corresponding rest-frame magnitudes. The FEM algorithm distinguishes four main groups: (1) red, (2) green, (3) blue, and (4) outliers. Each group is further divided into 3, 3, 4, and 2 subclasses, respectively. The accuracy of reproducing galaxy classes using spectroscopic data is high: 92%, 84%, 96% for red, green, and blue classes, respectively, except for dusty star-forming galaxies. The presented verification of the photometric classification demonstrates that large photometric samples can be used to distinguish different galaxy classes at z > 0.5 with an accuracy provided so far only by spectroscopic data except for particular galaxy classes. △ Less

Submitted 18 December, 2018; v1 submitted 24 May, 2018; originally announced May 2018.

Comments: the letter submitted to MNRAS, 6 pages, 4 figures. Please contact Malgorzata Siudek (msiudek@ifae.es) if you are interested in the catalogue

The VIMOS Public Extragalactic Redshift Survey (VIPERS). The complexity of galaxy populations at 0.4< z<1.3 revealed with unsupervised machine-learning algorithms

Authors: M. Siudek, K. Małek, A. Pollo, T. Krakowski, A. Iovino, M. Scodeggio, T. Moutard, G. Zamorani, L. Guzzo, B. Garilli, B. R. Granett, M. Bolzonella, S. de la Torre, U. Abbas, C. Adami, D. Bottini, A. Cappi, O. Cucciati, I. Davidzon, P. Franzetti, A. Fritz, J. Krywult, V. Le Brun, O. Le Fèvre, D. Maccagni , et al. (15 additional authors not shown)

Abstract: Various galaxy classification schemes have been developed so far to constrain the main physical processes regulating evolution of different galaxy types. In the era of a deluge of astrophysical information and recent progress in machine learning, a new approach to galaxy classification becomes imperative. We employ a Fisher Expectation-Maximization unsupervised algorithm working in a parameter s… ▽ More Various galaxy classification schemes have been developed so far to constrain the main physical processes regulating evolution of different galaxy types. In the era of a deluge of astrophysical information and recent progress in machine learning, a new approach to galaxy classification becomes imperative. We employ a Fisher Expectation-Maximization unsupervised algorithm working in a parameter space of 12 rest-frame magnitudes and spectroscopic redshift. The model (DBk) and the number of classes (12) were established based on the joint analysis of standard statistical criteria and confirmed by the analysis of the galaxy distribution with respect to a number of classes and their properties. This new approach allows us to classify galaxies based just on their redshifts and UV-NIR spectral energy distributions. The FEM unsupervised algorithm has automatically distinguished 12 classes: 11 classes of VIPERS galaxies and an additional class of broad-line AGNs. After a first broad division into blue, green and red categories we obtained a further sub-division into three red, three green, and five blue galaxy classes. The FEM classes follow the galaxy sequence from the earliest to the latest types that is reflected in their colours (which are constructed from rest-frame magnitudes used in classification procedure) but also their morphological, physical, and spectroscopic properties (not included in the classification scheme). We demonstrate that the members of each class share similar physical and spectral properties. In particular, we are able to find three different classes of red passive galaxy populations. Thus, we demonstrate the potential of an unsupervised approach to galaxy classification and we retrieve the complexity of galaxy populations at z~0.7, a task that usual simpler colour-based approaches cannot fulfil. △ Less

Submitted 18 December, 2018; v1 submitted 24 May, 2018; originally announced May 2018.

Comments: published in A&A, 27 pages, 10 figures, 4 appendixes. Please contact Malgorzata Siudek (msiudek@ifae.es) if you are interested in the catalogue

Journal ref: A&A 617, A70 (2018)

GAMA/H-ATLAS: The Local Dust Mass Function and Cosmic Density as a Function of Galaxy Type - A Benchmark for Models of Galaxy Evolution

Authors: R. A. Beeston, A. H. Wright, S. Maddox, H. L. Gomez, L. Dunne, S. P. Driver, A. Robotham, C. J. R. Clark, K. Vinsen, T. T. Takeuchi, G. Popping, N. Bourne, M. N. Bremer, S. Phillipps, A. J. Moffett, M. Baes, S. Brough, P. De Vis, S. A. Eales, B. W. Holwerda, J. Loveday, M. W. L. Smith, D. J. B. Smith, C. Vlahakis, L. Wang

Abstract: We present the dust mass function (DMF) of 15,750 galaxies with redshift $z< 0.1$, drawn from the overlapping area of the GAMA and {\it H-}ATLAS surveys. The DMF is derived using the density corrected $V_{\rm max}$ method, where we estimate $V_{\rm max}$ using: (i) the normal photometric selection limit ($pV_{\rm max}$) and (ii) a bivariate brightness distribution (BBD) technique, which accounts f… ▽ More We present the dust mass function (DMF) of 15,750 galaxies with redshift $z< 0.1$, drawn from the overlapping area of the GAMA and {\it H-}ATLAS surveys. The DMF is derived using the density corrected $V_{\rm max}$ method, where we estimate $V_{\rm max}$ using: (i) the normal photometric selection limit ($pV_{\rm max}$) and (ii) a bivariate brightness distribution (BBD) technique, which accounts for two selection effects. We fit the data with a Schechter function, and find $M^{*}=(4.65\pm0.18)\times 10^{7}\,h^2_{70}\, M_{\odot}$, $α=(1.22\pm 0.01)$, $φ^{*}=(6.26\pm 0.28)\times 10^{-3}\,h^3_{70}\,\rm Mpc^{-3}\,dex^{-1}$. The resulting dust mass density parameter integrated down to $10^4\,M_{\odot}$ is $Ω_{\rm d}=(1.11 \pm0.02)\times 10^{-6}$ which implies the mass fraction of baryons in dust is $f_{m_b}=(2.40\pm0.04)\times 10^{-5}$; cosmic variance adds an extra 7-17\,per\,cent uncertainty to the quoted statistical errors. Our measurements have fewer galaxies with high dust mass than predicted by semi-analytic models. This is because the models include too much dust in high stellar mass galaxies. Conversely, our measurements find more galaxies with high dust mass than predicted by hydrodynamical cosmological simulations. This is likely to be from the long timescales for grain growth assumed in the models. We calculate DMFs split by galaxy type and find dust mass densities of $Ω_{\rm d}=(0.88\pm0.03)\times 10^{-6}$ and $Ω_{\rm d}=(0.060\pm0.005)\times 10^{-6}$ for late-types and early-types respectively. Comparing to the equivalent galaxy stellar mass functions (GSMF) we find that the DMF for late-types is well matched by the GMSF scaled by $(8.07\pm0.35) \times 10^{-4}$. △ Less

Submitted 1 June, 2018; v1 submitted 19 December, 2017; originally announced December 2017.

Comments: 25 Pages, 18 Figures. Submitted December 2017

Overview of North Ecliptic Pole Deep multi-wavelength Survey (NEP-Deep)

Authors: H. Matsuhara, T. Wada, N. Oi, T. Takagi, T. Nakagawa, K. Murata, T. Goto, S. Oyabu, T. T. Takeuchi, K. Malek, A. Solarz, Y. Ohyama, T. Miyaji, M. Krumpe, H. M. Lee, M. Im, S. Serjeant, C. P. Pearson, G. J. White, M. A. Malkan, H. Hanami, T. Ishigaki, D. Burgarella, V. Buat, A. Pollo

Abstract: The recent updates of the North Ecliptic Pole deep (0.5~deg$^2$, NEP-Deep) multi-wavelength survey covering from X-ray to radio-wave is presented. The NEP-Deep provides us with several thousands of 15~$μ$m or 18~$μ$m selected sample of galaxies, which is the largest sample ever made at this wavelengths. A continuous filter coverage in the mid-infrared wavelength (7, 9, 11, 15, 18, and 24~$μ$m) is… ▽ More The recent updates of the North Ecliptic Pole deep (0.5~deg$^2$, NEP-Deep) multi-wavelength survey covering from X-ray to radio-wave is presented. The NEP-Deep provides us with several thousands of 15~$μ$m or 18~$μ$m selected sample of galaxies, which is the largest sample ever made at this wavelengths. A continuous filter coverage in the mid-infrared wavelength (7, 9, 11, 15, 18, and 24~$μ$m) is unique and vital to diagnose the contributions from starbursts and AGNs in the galaxies out to $z$=2.The new goal of the project is to resolve the nature of the cosmic star formation history at the violent epoch (e.g. $z$=1--2), and to find a clue to understand its decline from $z$=1 to present universe by utilizing the unique power of the multi-wavelength survey. The progress in this context is briefly mentioned. △ Less

Submitted 28 October, 2017; originally announced October 2017.

Comments: 4 pages, 2 figures, AKARI conference proceedings; accepted to PKAS

Journal ref: 2017, Publications of The Korean Astronomical Society, vol. 32, issue 1, pp. 213-217