Showing 1–50 of 194 results for author: Akiyama, M

Ongoing and fossil large-scale outflows detected in a high-redshift radio galaxy: [C II] observations of TN J0924$-$2201 at $z=5.174$

Authors: Kianhong Lee, Masayuki Akiyama, Kotaro Kohno, Daisuke Iono, Masatoshi Imanishi, Bunyo Hatsukade, Hideki Umehata, Tohru Nagao, Yoshiki Toba, Xiaoyang Chen, Fumi Egusa, Kohei Ichikawa, Takuma Izumi, Naoki Matsumoto, Malte Schramm, Kenta Matsuoka

Abstract: We present Atacama Large Millimeter/submillimeter Array observations of the [C II] 158 $μ$m line and the underlying continuum emission of TN J0924$-$2201, which is one of the most distant known radio galaxies at $z>5$. The [C II] line and 1-mm continuum emission are detected at the host galaxy. The systemic redshift derived from the [C II] line is $z_{\rm [C II]}=5.1736\pm0.0002$, indicating that… ▽ More We present Atacama Large Millimeter/submillimeter Array observations of the [C II] 158 $μ$m line and the underlying continuum emission of TN J0924$-$2201, which is one of the most distant known radio galaxies at $z>5$. The [C II] line and 1-mm continuum emission are detected at the host galaxy. The systemic redshift derived from the [C II] line is $z_{\rm [C II]}=5.1736\pm0.0002$, indicating that the Ly$α$ line is redshifted by a velocity of $1035\pm10$ km s$^{-1}$, marking the largest velocity offset between the [C II] and Ly$α$ lines recorded at $z>5$ to date. In the central region of the host galaxy, we identified a redshifted substructure of [C II] with a velocity of $702\pm17$ km s$^{-1}$, which is close to the CIV line with a velocity of $500\pm10$ km s$^{-1}$. The position and the velocity offsets align with a model of an outflowing shell structure, consistent with the large velocity offset of Ly$α$. The non-detection of [C II] and dust emission from the three CO(1--0)-detected companions indicates their different nature compared to dwarf galaxies based on the photodissociation region model. Given their large velocity of $\sim1500$ km s$^{-1}$, outflowing molecular clouds induced by the AGN is the most plausible interpretation, and they may exceed the escape velocity of a $10^{13}\,M_{\odot}$ halo. These results suggest that TN J0924$-$2201, with the ongoing and fossil large-scale outflows, is in a distinctive phase of removing molecular gas from a central massive galaxy in an overdense region in the early universe. A dusty HI absorber at the host galaxy is an alternative interpretation. △ Less

Submitted 1 July, 2024; originally announced July 2024.

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

Discovery of merging twin quasars at z = 6.05

Authors: Yoshiki Matsuoka, Takuma Izumi, Masafusa Onoue, Michael A. Strauss, Kazushi Iwasawa, Nobunari Kashikawa, Masayuki Akiyama, Kentaro Aoki, Junya Arita, Masatoshi Imanishi, Rikako Ishimoto, Toshihiro Kawaguchi, Kotaro Kohno, Chien-Hsiu Lee, Tohru Nagao, John D. Silverman, Yoshiki Toba

Abstract: We report the discovery of two quasars at a redshift of $z$ = 6.05, in the process of merging. They were serendipitously discovered from the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The quasars, HSC $J$121503.42$-$014858.7 (C1) and HSC $J$121503.55$-$014859.3 (C2), both have luminous ($>$10$^{43}$ erg s$^{-1}$) Ly$α$ emission with a cle… ▽ More We report the discovery of two quasars at a redshift of $z$ = 6.05, in the process of merging. They were serendipitously discovered from the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The quasars, HSC $J$121503.42$-$014858.7 (C1) and HSC $J$121503.55$-$014859.3 (C2), both have luminous ($>$10$^{43}$ erg s$^{-1}$) Ly$α$ emission with a clear broad component (full width at half maximum $>$1000 km s$^{-1}$). The rest-frame ultraviolet (UV) absolute magnitudes are $M_{1450} = -23.106 \pm 0.017$ (C1) and $-22.662 \pm 0.024$ (C2). Our crude estimates of the black hole masses provide $\log (M_{\rm BH}/M_\odot) = 8.1 \pm 0.3$ in both sources. The two quasars are separated by 12 kpc in projected proper distance, bridged by a structure in the rest-UV light suggesting that they are undergoing a merger. This pair is one of the most distant merging quasars reported to date, providing crucial insight into galaxy and black hole build-up in the hierarchical structure formation scenario. A companion paper will present the gas and dust properties captured by Atacama Large Millimeter/submillimeter Array observations, which provide additional evidence for and detailed measurements of the merger and also demonstrate that the two sources are not gravitationally-lensed images of a single quasar. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Comments: Published from ApJL. A companion paper (T. Izumi et al.) presents the associated ALMA observations

Journal ref: The Astrophysical Journal Letters, Volume 965, Issue 1, id.L4, 8 pp (2024)

Katachi: Decoding the Imprints of Past Star Formation on Present Day Morphology in Galaxies with Interpretable CNNs

Authors: Juan Pablo Alfonzo, Kartheik G. Iyer, Masayuki Akiyama, Greg L. Bryan, Suchetha Cooray, Eric Ludwig, Lamiya Mowla, Kiyoaki C. Omori, Camilla Pacifici, Joshua S. Speagle, John F. Wu

Abstract: The physical processes responsible for shaping how galaxies form and quench over time leave imprints on both the spatial (galaxy morphology) and temporal (star formation history; SFH) tracers that we use to study galaxies. While the morphology-SFR connection is well studied, the correlation with past star formation activity is not as well understood. To quantify this we present Katachi, an interpr… ▽ More The physical processes responsible for shaping how galaxies form and quench over time leave imprints on both the spatial (galaxy morphology) and temporal (star formation history; SFH) tracers that we use to study galaxies. While the morphology-SFR connection is well studied, the correlation with past star formation activity is not as well understood. To quantify this we present Katachi, an interpretable convolutional neural network (CNN) framework that learns the connection between the factors regulating star formation in galaxies on different spatial and temporal scales. Katachi is trained on 9904 galaxies at 0.02$<$z$<$0.1 in the SDSS-IV MaNGA DR17 sample to predict stellar mass (M$_*$; RMSE 0.22 dex), current star formation rate (SFR; RMSE 0.31 dex) and half-mass time (t$_{50}$; RMSE 0.23 dex). This information allows us to reconstruct non-parametric SFHs for each galaxy from \textit{gri} imaging alone. To quantify the morphological features informing the SFH predictions we use SHAP (SHapley Additive exPlanations). We recover the expected trends of M$_*$ governed by the growth of galaxy bulges, and SFR correlating with spiral arms and other star-forming regions. We also find the SHAP maps of D4000 are more complex than those of M$_*$ and SFR, and that morphology is correlated with t$_{50}$ even at fixed mass and SFR. Katachi serves as a scalable public framework to predict galaxy properties from large imaging surveys including Rubin, Roman, and Euclid, with large datasets of high SNR imaging across limited photometric bands. △ Less

Submitted 7 April, 2024; originally announced April 2024.

Comments: 35 pages, 25 figures. Accepted for publication in ApJ, comments are welcomed

Observational Properties of AGN Obscuration During the Peak of Accretion Growth

Authors: Bovornpratch Vijarnwannaluk, Masayuki Akiyama, Malte Schramm, Yoshihiro Ueda, Yoshiki Matsuoka, Yoshiki Toba, Naoki Matsumoto, Angel Ruiz, Ioannis Georgantopoulos, Ektoras Pouliasis, Elias Koulouridis, Kohei Ichikawa, Marcin Sawicki, Stephen Gwyn

Abstract: We investigated the gas obscuration and host galaxy properties of active galactic nuclei (AGN) during the peak of cosmic accretion growth of supermassive black holes (SMBHs) at redshift 0.8-1.8 using X-ray detected AGN with mid-infrared and far-infrared detection. The sample was classified as type-1 and type-2 AGN using optical spectral and morphological classification while the host galaxy proper… ▽ More We investigated the gas obscuration and host galaxy properties of active galactic nuclei (AGN) during the peak of cosmic accretion growth of supermassive black holes (SMBHs) at redshift 0.8-1.8 using X-ray detected AGN with mid-infrared and far-infrared detection. The sample was classified as type-1 and type-2 AGN using optical spectral and morphological classification while the host galaxy properties were estimated with multiwavelength SED fitting. For type-1 AGN, the black hole mass was determined from MgII emission lines while the black hole mass of type-2 AGN was inferred from the host galaxy's stellar mass. Based on the derived parameters, the distribution of the sample in the absorption hydrogen column density ($N_{\rm H}$) vs. Eddington ratio diagram is examined. Among the type-2 AGN, $28\pm5$\% are in the forbidden zone, where the obscuration by dust torus cannot be maintained due to radiation pressure on dusty material. The fraction is higher than that observed in the local universe from the BAT AGN Spectroscopic Survey (BASS) data release 2 ($11\pm3$\%). The higher fraction implies that the obscuration of the majority of AGN is consistent with the radiation pressure regulated unified model but with an increased incidence of interstellar matter (ISM) obscured AGN. We discuss the possibility of dust-free absorption in type-1 AGN and heavy ISM absorption in type-2 AGN. We also find no statistical difference in the star-formation activity between type-1 and type-2 AGN which may suggest that obscuration triggered by a gas-rich merging is not common among X-ray detected AGN in this epoch. △ Less

Submitted 10 March, 2024; originally announced March 2024.

Comments: Accepted for publication in MNRAS

AGN X-ray luminosity function and absorption function in the Early Universe ($3\leq z \leq 6$)

Authors: E. Pouliasis, A. Ruiz, I. Georgantopoulos, F. Vito, R. Gilli, C. Vignali, Y. Ueda, E. Koulouridis, M. Akiyama, S. Marchesi, B. Laloux, T. Nagao, S. Paltani, M. Pierre, Y. Toba, M. Habouzit, B. Vijarnwannaluk, C. Garrel

Abstract: The XLF of AGN offers a robust tool to study the evolution and the growth of SMBHs over cosmic time. Owing to the limited area probed by X-ray surveys, optical surveys are routinely used to probe the accretion in the high redshift Universe $z\geq 3$. However, optical surveys may be incomplete because they are strongly affected by dust redenning. In this work, we derive the XLF and its evolution at… ▽ More The XLF of AGN offers a robust tool to study the evolution and the growth of SMBHs over cosmic time. Owing to the limited area probed by X-ray surveys, optical surveys are routinely used to probe the accretion in the high redshift Universe $z\geq 3$. However, optical surveys may be incomplete because they are strongly affected by dust redenning. In this work, we derive the XLF and its evolution at high redshifts using a large sample of AGNs selected in different fields with various areas and depths covering a wide range of luminosities. Additionally, we put the tightest yet constraints on the absorption function in this redshift regime. In particular, we use more than 600 soft X-ray selected high-z sources in the Chandra Deep fields, the Chandra COSMOS Legacy survey and the XMM-XXL northern field. We derive the X-ray spectral properties for all sources via spectral fitting, using a consistent technique and model. For modeling the parametric form of the XLF and the absorption function, we use a Bayesian methodology allowing us to correctly propagate the uncertainties for the observed X-ray properties of our sources and also the absorption effects. The evolution of XLF is in agreement with a pure density evolution model similar to what is witnessed at optical wavelengths, although a luminosity dependent density evolution model cannot be securely ruled out. A large fraction ($60\%)$ of our sources are absorbed by column densities of $\rm N_H \geq 10^{23} cm^{-2} $, while $17$\% of the sources are CTK. Our results favor a scenario where both the ISM of the host and the AGN torus contribute to the obscuration. The derived BHAD is in agreement with the simulations, if one takes into account that the X-ray AGN are hosted by massive galaxies, while it differs from the one derived using JWST data. The latter could be due to the differences in the AGN and host-galaxy properties. △ Less

Submitted 24 January, 2024; originally announced January 2024.

Comments: Accepted for publication in A&A. 25 pages, 15 figures (+4 in Appendix), 2 tables (+1 in Appendix)

Journal ref: A&A 685, A97 (2024)

Black hole mass and Eddington ratio distributions of less-luminous quasars at $z\sim4$ in the Subaru Hyper Suprime-Cam Wide field

Authors: Wanqiu He, Masayuki Akiyama, Motohiro Enoki, Kohei Ichikawa, Kohei Inayoshi, Nobunari Kashikawa, Toshihiro Kawaguchi, Yoshiki Matsuoka, Tohru Nagao, Masafusa Onoue, Taira Oogi, Andreas Schulze, Yoshiki Toba, Yoshihiro Ueda

Abstract: We investigate the black hole mass function (BHMF) and Eddington ratio distribution function (ERDF) of broad-line AGNs at z=4, based on a sample of 52 quasars with i<23.2 at 3.50 < z < 4.25 from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) S16A-Wide2 dataset, and 1,462 quasars with i<20.2 in the same redshift range from the Sloan Digital Sky Survey (SDSS) DR7 quasar catalog. Virial BH… ▽ More We investigate the black hole mass function (BHMF) and Eddington ratio distribution function (ERDF) of broad-line AGNs at z=4, based on a sample of 52 quasars with i<23.2 at 3.50 < z < 4.25 from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) S16A-Wide2 dataset, and 1,462 quasars with i<20.2 in the same redshift range from the Sloan Digital Sky Survey (SDSS) DR7 quasar catalog. Virial BH masses of quasars are estimated using the width of the CIV 1549Å line and the continuum luminosity at 1350Å. To obtain the intrinsic broad-line AGN BHMF and ERDF, we correct for the incompleteness in the low-mass and/or low-Eddington-ratio ranges caused by the flux-limited selection. The resulting BHMF is constrained down to $\log M_{\rm BH}/M_{\odot}\sim7.5$. In comparison with broad-line AGN BHMFs at z=2 in the literature, we find that the number density of massive SMBHs peaks at higher redshifts, consistent with the "down-sizing" evolutionary scenario. Additionally, the resulting ERDF shows a negative dependence on BH mass, suggesting more massive SMBHs tend to accrete at lower Eddington ratios at z=4. With the derived intrinsic broad-line AGN BHMF, we also evaluate the active fraction of broad-line AGNs among the entire SMBH population at z=4. The resulting active fraction may suggest a positive dependence on BH mass. Finally, we examine the time evolution of broad-line AGN BHMF between z=4 and 6 through solving the continuity equation. The results suggest that the broad-line AGN BHMFs at z=4-6 only show evolution in their normalization, but with no significant changes in their shape. △ Less

Submitted 15 November, 2023; originally announced November 2023.

Comments: 38 pages, 26 figures. Submitted to ApJ. Data is available in https://github.com/wanqqq31/BHMF_hsc_z4

Subaru High-$z$ Exploration of Low-luminosity Quasars (SHELLQs). XVII. Black Hole Mass Distribution at $z \sim 6$ Estimated via Spectral Comparison with Low-$z$ Quasars

Authors: Ayumi Takahashi, Yoshiki Matsuoka, Masafusa Onoue, Michael A. Strauss, Nobunari Kashikawa, Yoshiki Toba, Kazushi Iwasawa, Masatoshi Imanishi, Masayuki Akiyama, Toshihiro Kawaguchi, Akatoki Noboriguchi, Chien-Hsiu Lee

Abstract: We report the distribution of black hole (BH) masses and Eddingont ratios estimated for a sample of 131 low luminosity quasars in the early cosmic epoch ($5.6 < z < 7.0$). Our work is based on Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, which has constructed a low luminosity quasar sample down to $M_{1450} \sim - 21$ mag, exploiting the survey data of Hyper Suprime-Cam… ▽ More We report the distribution of black hole (BH) masses and Eddingont ratios estimated for a sample of 131 low luminosity quasars in the early cosmic epoch ($5.6 < z < 7.0$). Our work is based on Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, which has constructed a low luminosity quasar sample down to $M_{1450} \sim - 21$ mag, exploiting the survey data of Hyper Suprime-Cam installed on Subaru Telescope. The discovery spectra of these quasars are limited to the rest-frame wavelengths of $\sim$ 1200 -- 1400 Å, which contains no emission lines that can be used as BH mass estimators. In order to overcome this problem, we made use of low-$z$ counterpart spectra from the Sloan Digital Sky Survey, which are spectrally matched to the high-$z$ spectra in overlapping wavelengths. We then combined the C~{\sc iv} emission line widths of the counterparts with the continuum luminosity from the SHELLQs data to estimate BH masses. The resulting BH mass distribution has a range of $\sim 10^{7-10} M_{\odot}$, with most of the quasars having BH masses $\sim 10^{8} M_{\odot}$ with sub-Eddington accretion. The present study provides not only a new insight into normal quasars in the reionization epoch, but also a new promising way to estimate BH masses of high-$z$ quasars without near-infrared spectroscopy. △ Less

Submitted 18 October, 2023; originally announced October 2023.

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

SHARPEST: The atmospheric turbulence profiling experiment using Shack-Hartmann sensor at the Subaru telescope

Authors: Hajime Ogane, Yoshito Ono, Yosuke Minowa, Shin Oya, Koki Terao, Takumi Akasawa, Riki Homan, Masayuki Akiyama

Abstract: Atmospheric turbulence profile plays an important role in designing and operating adaptive optics (AO) systems with multiple laser guide stars. To obtain representative free atmospheric profiles and resolved ground layer profiles for future AO systems at the Subaru telescope, we are conducting the SHARPEST (Shack-Hartmann Atmospheric tuRbulence Profiling Experiment at the Subaru Telescope) project… ▽ More Atmospheric turbulence profile plays an important role in designing and operating adaptive optics (AO) systems with multiple laser guide stars. To obtain representative free atmospheric profiles and resolved ground layer profiles for future AO systems at the Subaru telescope, we are conducting the SHARPEST (Shack-Hartmann Atmospheric tuRbulence Profiling Experiment at the Subaru Telescope) project. In this project, we develop a turbulence profiler comprising two Shack-Hartmann (SH) sensors to observe a pair of bright stars through the Subaru telescope with high spatial sampling by 2 cm subapertures. We perform two analyses on the SH spot data: variance analysis on the spot scintillation for free atmospheric profiles, and on the spot slope for ground layer profiles. This paper introduces the initial results of free atmospheric profiles as well as total seeing values and wind profiles obtained by the first two engineering runs. The free atmospheric profiles reconstructed by the two independent SH sensors show good agreement. The results are also consistent with simultaneous measurements by another profiler except for turbulence strength at ~1 km, which is explained by an overestimation problem of scintillation-based profilers. Measured total seeing values are also smaller than the simultaneous measurements, possibly due to the difference in ground layer turbulence between the two sites. The wind profiles show good consistency with the direct measurements by a rawinsonde. Through this study, we establish a method to constrain the free atmospheric profile, the total seeing, and the wind profile by analysing data from a single SH sensor with fine subapertures. △ Less

Submitted 14 October, 2023; originally announced October 2023.

Comments: 16 pages, 17 figures, accepted for publication in MNRAS

Reconstruction of Cosmic Black Hole Growth and Mass Distribution from Quasar Luminosity Functions at $z>4$: Implications for Faint and Low-mass Populations in JWST

Authors: Wenxiu Li, Kohei Inayoshi, Masafusa Onoue, Wanqiu He, Yoshiki Matsuoka, Zhiwei Pan, Masayuki Akiyama, Takuma Izumi, Tohru Nagao

Abstract: The evolution of the quasar luminosity function (QLF) is fundamental to understanding the cosmic evolution of black holes (BHs) through their accretion phases. In the era of the James Webb Space Telescope (JWST), Euclid, and Nancy Grace Roman Space Telescope, their unprecedented detection sensitivity and wide survey area can unveil the low-luminosity quasar and low-mass BH population, and provide… ▽ More The evolution of the quasar luminosity function (QLF) is fundamental to understanding the cosmic evolution of black holes (BHs) through their accretion phases. In the era of the James Webb Space Telescope (JWST), Euclid, and Nancy Grace Roman Space Telescope, their unprecedented detection sensitivity and wide survey area can unveil the low-luminosity quasar and low-mass BH population, and provide new insights into quasar host galaxies. We present a theoretical model describing BH growth from initial seeding at $z>20$ to $z\sim 4$,incorporating the duration of accretion episodes, the distribution of Eddington ratios, and the mass dependency of BH accretion rates. By constraining the model parameters with the observed QLFs at $4\leq z\leq6$ across a wide UV luminosity range ($-29<M_{\rm 1450}<-24$), we find that the high-redshift BH population grows rapidly at $z\gtrsim6$, and decelerates the pace in subsequent epochs. Toward lower redshifts ($z<6$), mass-dependent accretion inhibits the growth of high-mass BHs with $M_{\bullet}>10^8~M_\odot$, leading to mass saturation at $M_\bullet\gtrsim 10^{10}~M_\odot$. We predict the BH mass function down to $M_{\bullet}\sim 10^6~M_\odot$ for both unobscured and obscured quasar populations at $4\leq z \leq 11$, offering a benchmark for future observational tests. Our model accounts for the presence of both bright and faint quasars at $z>4$, including those discovered by JWST. Furthermore, our findings suggest two distinct pathways for the early assembly of the BH-galaxy mass correlation: the population with a BH-to-stellar mass ratio near the local value of $M_\bullet/M_{\star}\simeq5\times10^{-3}$ maintains a proximity to the relation through its evolution via moderate growth, while the population that begins to grow above the local relation accretes mass rapidly and becomes as overmassive as $M_\bullet/M_\star \sim 0.01-0.1$ by $z\sim 6$. △ Less

Submitted 30 April, 2024; v1 submitted 9 June, 2023; originally announced June 2023.

Comments: 16 pages, 8 figures; Accepted by ApJ

Quasar Luminosity Function at z = 7

Authors: Yoshiki Matsuoka, Masafusa Onoue, Kazushi Iwasawa, Michael A. Strauss, Nobunari Kashikawa, Takuma Izumi, Tohru Nagao, Masatoshi Imanishi, Masayuki Akiyama, John D. Silverman, Naoko Asami, James Bosch, Hisanori Furusawa, Tomotsugu Goto, James E. Gunn, Yuichi Harikane, Hiroyuki Ikeda, Kohei Inayoshi, Rikako Ishimoto, Toshihiro Kawaguchi, Satoshi Kikuta, Kotaro Kohno, Yutaka Komiyama, Chien-Hsiu Lee, Robert H. Lupton , et al. (19 additional authors not shown)

Abstract: We present the quasar luminosity function (LF) at $z = 7$, measured with 35 spectroscopically confirmed quasars at $6.55 < z < 7.15$. The sample of 22 quasars from the Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, combined with 13 brighter quasars in the literature, covers an unprecedentedly wide range of rest-frame ultraviolet magnitudes over $-28 < M_{1450} < -23$. We… ▽ More We present the quasar luminosity function (LF) at $z = 7$, measured with 35 spectroscopically confirmed quasars at $6.55 < z < 7.15$. The sample of 22 quasars from the Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, combined with 13 brighter quasars in the literature, covers an unprecedentedly wide range of rest-frame ultraviolet magnitudes over $-28 < M_{1450} < -23$. We found that the binned LF flattens significantly toward the faint end populated by the SHELLQs quasars. A maximum likelihood fit to a double power-law model has a break magnitude $M^*_{1450} = -25.60^{+0.40}_{-0.30}$, a characteristic density $Φ^* = 1.35^{+0.47}_{-0.30}$ Gpc$^{-3}$ mag$^{-1}$, and a bright-end slope $β= -3.34^{+0.49}_{-0.57}$, when the faint-end slope is fixed to $α= -1.2$ as observed at $z \le 6$. The overall LF shape remains remarkably similar from $z = 4$ to $7$, while the amplitude decreases substantially toward higher redshifts, with a clear indication of an accelerating decline at $z \ge 6$. The estimated ionizing photon density, $10^{48.2 \pm 0.1}$ s$^{-1}$ Mpc$^{-3}$, is less than 1 % of the critical rate to keep the intergalactic medium ionized at $z = 7$, and thus indicates that quasars are not a major contributor to cosmic reionization. △ Less

Submitted 18 May, 2023; originally announced May 2023.

Comments: The Astrophysical Journal Letters, in press

How WEIRD is Usable Privacy and Security Research? (Extended Version)

Authors: Ayako A. Hasegawa, Daisuke Inoue, Mitsuaki Akiyama

Abstract: In human factor fields such as human-computer interaction (HCI) and psychology, researchers have been concerned that participants mostly come from WEIRD (Western, Educated, Industrialized, Rich, and Democratic) countries. This WEIRD skew may hinder understanding of diverse populations and their cultural differences. The usable privacy and security (UPS) field has inherited many research methodolog… ▽ More In human factor fields such as human-computer interaction (HCI) and psychology, researchers have been concerned that participants mostly come from WEIRD (Western, Educated, Industrialized, Rich, and Democratic) countries. This WEIRD skew may hinder understanding of diverse populations and their cultural differences. The usable privacy and security (UPS) field has inherited many research methodologies from research on human factor fields. We conducted a literature review to understand the extent to which participant samples in UPS papers were from WEIRD countries and the characteristics of the methodologies and research topics in each user study recruiting Western or non-Western participants. We found that the skew toward WEIRD countries in UPS is greater than that in HCI. Geographic and linguistic barriers in the study methods and recruitment methods may cause researchers to conduct user studies locally. In addition, many papers did not report participant demographics, which could hinder the replication of the reported studies, leading to low reproducibility. To improve geographic diversity, we provide the suggestions including facilitate replication studies, address geographic and linguistic issues of study/recruitment methods, and facilitate research on the topics for non-WEIRD populations. △ Less

Submitted 8 October, 2023; v1 submitted 8 May, 2023; originally announced May 2023.

Comments: This paper is the extended version of the paper presented at USENIX SECURITY 2024

J1406+0102: Dust Obscured Galaxy Hiding Super Eddington Accretion System with Bright Radio Emission

Authors: Hikaru Fukuchi, Kohei Ichikawa, Masayuki Akiyama, Shigeo Kimura, Yoshiki Toba, Kohei Inayoshi, Akatoki Noboriguchi, Toshihiro Kawaguchi, Xiaoyang Chen, Itsna K. Fitriana

Abstract: Recent high-$z$ quasar observations strongly indicate that super-Eddington accretion is a crucial phase to describe the existence of supermassive black holes (SMBHs) with $M_\mathrm{BH} \gtrsim 10^9 M_\odot$ at $z \gtrsim 7$. Motivated by the theoretical suggestion that the super-Eddington phase efficiently produces outflows and jets bright in radio bands, we search and find a super-Eddington radi… ▽ More Recent high-$z$ quasar observations strongly indicate that super-Eddington accretion is a crucial phase to describe the existence of supermassive black holes (SMBHs) with $M_\mathrm{BH} \gtrsim 10^9 M_\odot$ at $z \gtrsim 7$. Motivated by the theoretical suggestion that the super-Eddington phase efficiently produces outflows and jets bright in radio bands, we search and find a super-Eddington radio-loud dust-obscured galaxy (DOG) J1406+0102 at $z=0.236$, through cross-matching of the infrared-bright DOGs of Noboriguchi et al. (2019) with the VLA/FIRST 1.4 GHz radio and the SDSS optical spectral catalog. DOG J1406+0102 shows broad components in the Balmer lines. Assuming those lines are from the broad line region, it gives BH mass estimation of $\log\ (M_\mathrm{BH}/M_\odot)=7.30 \pm 0.25$, and AGN luminosity of $\log (L_\mathrm{bol,[OIII]}/\mathrm{erg}~\mathrm{s}^{-1}) = 45.91\pm0.38$ estimated from the intrinsic [OIII] luminosity, resulting in super-Eddington accretion of $λ_\mathrm{Edd}\simeq 3$. We show that 1) DOG J1406+0102 is operating strong AGN feedback: the [OIII] outflow velocity exceeds the escape velocity of the host galaxy halo and the kinetic efficiency is obtained as $\approx$ 8% that can be sufficient to quench the host galaxy, 2) the expected future growth pathway of DOG J1406+0102 would join an over-massive BH trajectory and 3) radio-loud DOGs can provide a significant contribution to the high-energy ($\gtrsim$ 100 TeV) cosmic neutrino background if we assume DOG J1406+0102 as a representative of radio-loud DOGs. △ Less

Submitted 9 March, 2023; originally announced March 2023.

Comments: 10 pages, 5 figures, submitted to ApJL

eROSITA Final Equatorial-Depth Survey (eFEDS): eFEDS X-ray view of WERGS Radio Galaxies selected by the Subaru/HSC and VLA/FIRST survey

Authors: Kohei Ichikawa, Takuji Yamashita, Andrea Merloni, Junyao Li, Teng Liu, Mara Salvato, Masayuki Akiyama, Riccardo Arcodia, Tom Dwelly, Xiaoyang Chen, Masatoshi Imanishi, Kohei Inayoshi, Toshihiro Kawaguchi, Taiki Kawamuro, Mitsuru Kokubo, Yoshiki Matsuoka, Tohru Nagao, Malte Schramm, Hyewon Suh, Masayuki Tanaka, Yoshiki Toba, Yoshihiro Ueda

Abstract: We construct the eROSITA X-ray catalog of radio galaxies discovered by the WERGS survey that is made by the cross-matching of the wide-area Subaru/HSC optical survey and VLA/FIRST 1.4 GHz radio survey. We find 393 eROSITA detected radio galaxies in the 0.5--2 keV band in the eFEDS field covering 140~deg$^2$. Thanks to the wide and medium depth eFEDS X-ray survey, the sample contains the rare and m… ▽ More We construct the eROSITA X-ray catalog of radio galaxies discovered by the WERGS survey that is made by the cross-matching of the wide-area Subaru/HSC optical survey and VLA/FIRST 1.4 GHz radio survey. We find 393 eROSITA detected radio galaxies in the 0.5--2 keV band in the eFEDS field covering 140~deg$^2$. Thanks to the wide and medium depth eFEDS X-ray survey, the sample contains the rare and most X-ray luminous radio galaxies above the knee of the X-ray luminosity function, spanning 44<log L(0.5-2keV,abs)<46.5 at $1<z<4$. Based on the X-ray properties obtained by the spectral fitting, 37 sources show obscured AGN signature with $\log (N_\mathrm{H}/\mathrm{cm}^{-2})>22$. Those obscured and radio AGN reside in $0.4<z<3.2$, indicating that they are obscured counterparts of the radio-loud quasar, which are missed in the previous optical quasar surveys. By combining radio and X-ray luminosities, the jet production efficiency $η_\mathrm{jet}$ is investigated, and we find 14 sources with extremely high jet production efficiency at $η_\mathrm{jet}\approx1$. This high $η_\mathrm{jet}$ value might be a result of 1) the decreased radiation efficiency of $η_\mathrm{rad}<0.1$ due to the low accretion rate for those sources and/or 2) the boosting due to the decline of $L_\mathrm{bol}$ by a factor of 10--100 by keeping $P_\mathrm{jet}$ constant in the previous Myr, indicating the experience of the AGN feedback. Finally, inferring the BH masses from the stellar-mass, we find that X-ray luminous sources show the excess of the radio emission with respect to the value estimated from the fundamental plane. Such radio emission excess cannot be explained by the Doppler booming alone, and therefore disk-jet connection of X-ray luminous eFEDS-WERGS is fundamentally different from the conventional fundamental plane which mainly covers low accretion regime. △ Less

Submitted 15 February, 2023; originally announced February 2023.

Comments: 20 pages, 11 figures, accepted for publication in A&A

Journal ref: A&A 672, A171 (2023)

The cosmic web of X-ray active galactic nuclei seen through the eROSITA Final Equatorial Depth Survey (eFEDS)

Authors: Johan Comparat, Wentao Luo, Andrea Merloni, Surhud More, Mara Salvato, Mirko Krumpe, Takamitsu Miyaji, William Brandt, Antonis Georgakakis, Masayuki Akiyama, Johannes Buchner, Tom Dwelly, Toshihiro Kawaguchi, Teng Liu, Tohru Nagao, Kirpal Nandra, John Silverman, Yoshiki Toba, Scott F. Anderson, Juna Kollmeier

Abstract: Which galaxies in the general population turn into active galactic nuclei (AGNs) is a keystone of galaxy formation and evolution. Thanks to SRG/eROSITA's contiguous 140 square degree pilot survey field, we constructed a large, complete, and unbiased soft X-ray flux-limited ($F_X>6.5\times 10^{-15}$ erg s$^{-1}$ cm$^{-2}$) AGN sample at low redshift, $0.05<z<0.55$. Two summary statistics, the clust… ▽ More Which galaxies in the general population turn into active galactic nuclei (AGNs) is a keystone of galaxy formation and evolution. Thanks to SRG/eROSITA's contiguous 140 square degree pilot survey field, we constructed a large, complete, and unbiased soft X-ray flux-limited ($F_X>6.5\times 10^{-15}$ erg s$^{-1}$ cm$^{-2}$) AGN sample at low redshift, $0.05<z<0.55$. Two summary statistics, the clustering using spectra from SDSS-V and galaxy-galaxy lensing with imaging from HSC, are measured and interpreted with halo occupation distribution and abundance matching models. Both models successfully account for the observations. We obtain an exceptionally complete view of the AGN halo occupation distribution. The population of AGNs is broadly distributed among halos with a mean mass of $3.9 _{- 2.4 }^{+ 2.0 }\times10^{12}M_\odot$. This corresponds to a large-scale halo bias of $b(z=0.34)= 0.99 ^{+0.08}_{-0.10}$. The central occupation has a large transition parameter, $σ_{\log_{10}(M)}=1.28\pm0.2$. The satellite occupation distribution is characterized by a shallow slope, $α_{\rm sat}=0.73\pm0.38$. We find that AGNs in satellites are rare, with $f_{\rm sat}<20\%$. Most soft X-ray-selected AGNs are hosted by central galaxies in their dark matter halo. A weak correlation between soft X-ray luminosity and large-scale halo bias is confirmed (3.3$σ$). We discuss the implications of environmental-dependent AGN triggering. This study paves the way toward fully charting, in the coming decade, the coevolution of X-ray AGNs, their host galaxies, and dark matter halos by combining eROSITA with SDSS-V, 4MOST, DESI, LSST, and \textit{Euclid} data. △ Less

Submitted 14 April, 2023; v1 submitted 3 January, 2023; originally announced January 2023.

Comments: 17 pages, 9 figures, 2 tables. A&A: Received: 19 December 2022 / Accepted: 21 March 2023

Journal ref: A&A 673, A122 (2023)

H1821+643: The most X-ray and infrared luminous AGN in the Swift/BAT survey in the process of rapid stellar and supermassive black hole mass assembly

Authors: Hikaru Fukuchi, Kohei Ichikawa, Masayuki Akiyama, Claudio Ricci, Sunmyon Chon, Mitsuru Kokubo, Ang Liu, Takuya Hashimoto, Takuma Izumi

Abstract: H1821+643 is the most X-ray luminous non-beamed AGN of $L_\mathrm{14-150 keV}= 5.2\times 10^{45}$ erg s$^{-1}$ in the Swift/BAT ultra-hard X-ray survey and it is also a hyper-luminous infrared (IR) galaxy $L_\mathrm{IR} = 10^{13.2} L_\odot$ residing in the center of a massive galaxy cluster, which is a unique environment achieving the rapid mass assembly of black holes (BH) and host galaxies in th… ▽ More H1821+643 is the most X-ray luminous non-beamed AGN of $L_\mathrm{14-150 keV}= 5.2\times 10^{45}$ erg s$^{-1}$ in the Swift/BAT ultra-hard X-ray survey and it is also a hyper-luminous infrared (IR) galaxy $L_\mathrm{IR} = 10^{13.2} L_\odot$ residing in the center of a massive galaxy cluster, which is a unique environment achieving the rapid mass assembly of black holes (BH) and host galaxies in the local universe. We decompose the X-ray to IR spectral energy distribution (SED) into the AGN and starburst component using the SED fitting tool CIGALE-2022.0 and show that H1821+643 consumes a large amount of cold gas ($\dot{M}_\mathrm{con}$) with star-formation rate of $\log ( \mathrm{SFR}/M_{\odot}~\mathrm{yr}^{-1}) = 3.01 \pm 0.04$ and BH accretion rate of $\log (\dot{M}_\mathrm{BH}/M_{\odot}~\mathrm{yr}^{-1}) = 1.20 \pm 0.05$. This high $\dot{M}_\mathrm{con}$ is larger than the cooling rate ($\dot{M}_\mathrm{cool}$) of the intra-cluster medium (ICM), $\dot{M}_\mathrm{con}/\dot{M}_\mathrm{cool} \gtrsim 1$, which is one to two order magnitude higher than the typical value of other systems, indicating that H1821 provides the unique and extreme environment of rapid gas consumption. We also show that H1821+643 has an efficient cooling path achieving from $10^7$ K to $10^2$ K thanks to [OIII] 63 $μ\mathrm{m}$, which is a main coolant in low temperature range ($10^4$ K to $10^2$ K) with a cooling rate of $\dot{M}_{\mathrm{cool}}=3.2\times 10^5\ M_{\odot}\mathrm{~yr^{-1}}$, and the star-forming region extends over 40 kpc scale. △ Less

Submitted 19 September, 2022; originally announced September 2022.

Comments: 20 pages, 8 figures, accepted for publication in ApJ

The Obscured Fraction of Quasars at Cosmic Noon

Authors: Bovornpratch Vijarnwannaluk, Masayuki Akiyama, Malte Schramm, Yoshihiro Ueda, Yoshiki Matsuoka, Yoshiki Toba, Marcin Sawicki, Stephen Gwyn, Janek Pflugradt

Abstract: Statistical studies of X-ray selected Active Galactic Nuclei (AGN) indicate that the fraction of obscured AGN increases with increasing redshift, and the results suggest that a significant part of the accretion growth occurs behind obscuring material in the early universe. We investigate the obscured fraction of highly accreting X-ray AGN at around the peak epoch of supermassive black hole growth… ▽ More Statistical studies of X-ray selected Active Galactic Nuclei (AGN) indicate that the fraction of obscured AGN increases with increasing redshift, and the results suggest that a significant part of the accretion growth occurs behind obscuring material in the early universe. We investigate the obscured fraction of highly accreting X-ray AGN at around the peak epoch of supermassive black hole growth utilizing the wide and deep X-ray and optical/IR imaging datasets. A unique sample of luminous X-ray selected AGNs above $z>2$ was constructed by matching the XMM-SERVS X-ray point-source catalog with a PSF-convolved photometric catalog covering from $u^*$ to 4.5$μ\mathrm{m}$ bands. Photometric redshift, hydrogen column density, and 2-10 keV AGN luminosity of the X-ray selected AGN candidates were estimated. Using the sample of 306 2-10 keV detected AGN at above redshift 2, we estimate the fraction of AGN with $\log N_{\rm H}\ (\rm cm^{-2})>22$, assuming parametric X-ray luminosity and absorption functions. The results suggest that $76_{-3}^{+4}\%$ of luminous quasars ($\log L_X\ (\rm erg\ s^{-1}) >44.5$) above redshift 2 are obscured. The fraction indicates an increased contribution of obscured accretion at high redshift than that in the local universe. We discuss the implications of the increasing obscured fraction with increasing redshift based on the AGN obscuration scenarios, which describe obscuration properties in the local universe. Both the obscured and unobscured $z>2$ AGN show a broad range of SEDs and morphology, which may reflect the broad variety of host galaxy properties and physical processes associated with the obscuration. △ Less

Submitted 30 September, 2022; v1 submitted 16 September, 2022; originally announced September 2022.

Comments: 25 Pages, 19 Figure, ApJ accepted. Updated version contains corrections based on comments from the community

Finding of a population of active galactic nuclei showing a significant luminosity decline in the past $\sim 10^{3-4}$ yrs

Authors: J. Pflugradt, K. Ichikawa, M. Akiyama, M. Kokubo, B. Vijarnwannaluk, H. Noda, X. Chen

Abstract: Recent observations have revealed an interesting active galactic nuclei (AGN) subclass that shows strong activity at large scales ($\sim1$ kpc) but weaker at small scales ($<10$ pc), suggesting a strong change in the mass accretion rate of the central engine in the past $10^{3-4}$ yr. We systematically search for such declining or fading AGN by cross-matching the SDSS type-1 AGN catalog at… ▽ More Recent observations have revealed an interesting active galactic nuclei (AGN) subclass that shows strong activity at large scales ($\sim1$ kpc) but weaker at small scales ($<10$ pc), suggesting a strong change in the mass accretion rate of the central engine in the past $10^{3-4}$ yr. We systematically search for such declining or fading AGN by cross-matching the SDSS type-1 AGN catalog at $z<0.4$, covering the [OIII]$\lambda5007$ emission line which is a tracer for the narrow-line region (NLR) emission, with the WISE mid-infrared (MIR) catalog covering the emissions from the dusty tori. Out of the 7,653 sources, we found 57 AGN whose bolometric luminosities estimated from the MIR band are at least one order of magnitude fainter than those estimated from the [OIII]$\lambda5007$ emission line. This luminosity declining AGN candidate population shows four important properties: 1) the past AGN activity estimated from the [OIII]$\lambda5007$ line reaches around the Eddington-limit, 2) more than 30% of the luminosity declining AGN candidates show a large absolute variability of $ΔW1 > 0.45$ mag in the previous $\sim10$ yr at the WISE 3.4 $μ$m band, 3) the median ratio of $\log$([NII]$\lambda6584/ \mathrm{H}α\lambda6563)=-0.52$, suggesting a lower gas metallicity and/or higher ionization parameter compared to other AGN populations. 4) the second epoch spectra of the population indicate a spectral type change for 15% of the sources. This population provides insights on the possible connection between the luminosity decline which started $\sim10^{3-4}$ yr ago and the decline in the recent $10$ yr. △ Less

Submitted 25 August, 2022; originally announced August 2022.

Comments: 21 pages, 10 figures, 2 tables (1 MRT), accepted for publication in ApJ

Optical IFU Observations of GOALS Sample with KOOLS-IFU on Seimei Telescope: Initial results of 9 U/LIRGs at $z <$ 0.04

Authors: Yoshiki Toba, Satoshi Yamada, Kazuya Matsubayashi, Koki Terao, Aoi Moriya, Yoshihiro Ueda, Kouji Ohta, Aoi Hashiguchi, Kazuharu G. Himoto, Hideyuki Izumiura, Kazuma Joh, Nanako Kato, Shuhei Koyama, Hiroyuki Maehara, Rana Misato, Akatoki Noboriguchi, Shoji Ogawa, Naomi Ota, Mio Shibata, Nozomu Tamada, Anri Yanagawa, Naoki Yonekura, Tohru Nagao, Masayuki Akiyama, Masaru Kajisawa , et al. (1 additional authors not shown)

Abstract: We present ionized gas properties of 9 local ultra/luminous infrared galaxies (U/LIRGs) at $z <$ 0.04 through IFU observations with KOOLS-IFU on Seimei Telescope. The observed targets are drawn from the Great Observatories All-sky LIRG Survey (GOALS), covering a wide range of merger stages. We successfully detect emission lines such as H$β$, [OIII]$λ$5007, H$α$, [NII]$λλ$6549,6583, and [SII]$λλ$67… ▽ More We present ionized gas properties of 9 local ultra/luminous infrared galaxies (U/LIRGs) at $z <$ 0.04 through IFU observations with KOOLS-IFU on Seimei Telescope. The observed targets are drawn from the Great Observatories All-sky LIRG Survey (GOALS), covering a wide range of merger stages. We successfully detect emission lines such as H$β$, [OIII]$λ$5007, H$α$, [NII]$λλ$6549,6583, and [SII]$λλ$6717,6731 with a spectral resolution of $R$ = 1500-2000, which provides (i) spatially-resolved ($\sim$200-700 pc) moment map of ionized gas and (ii) diagnostics for active galactic nucleus (AGN) within the central $\sim$3--11 kpc in diameter for our sample. We find that [OIII] outflow that is expected to be driven by AGN tends to be stronger (i) towards the galactic center and (ii) as a sequence of merger stage. In particular, the outflow strength in the late-stage (stage D) mergers is about 1.5 times stronger than that in the early-state (stage B) mergers, which indicates that galaxy mergers could induce AGN-driven outflow and play an important role in the co-evolution of galaxies and supermassive black holes. △ Less

Submitted 24 August, 2022; originally announced August 2022.

Comments: 12 pages, 8 figures, and 2 tables, accepted for publication in PASJ

High Contrast and High Angular Imaging at Subaru Telescope

Authors: Olivier Guyon, Kyohoon Ahn, Masayuki Akiyama, Thayne Currie, Vincent Deo, Takashi Hattori, Tomoyuki Kudo, Julien Lozi, Yosuke Minowa, Yoshito Ono, Nour Skaf, Motohide Tamura, Vincent Vievard

Abstract: Adaptive Optics projects at Subaru Telescope span a wide field of capabilities ranging from ground-layer adaptive optics (GLAO) providing partial correction over a 20 arcmin FOV to extreme adaptive optics (ExAO) for exoplanet imaging. We describe in this paper current and upcoming narrow field-of-view capabilities provided by the Subaru Extreme Adaptive Optics Adaptive Optics (SCExAO) system and i… ▽ More Adaptive Optics projects at Subaru Telescope span a wide field of capabilities ranging from ground-layer adaptive optics (GLAO) providing partial correction over a 20 arcmin FOV to extreme adaptive optics (ExAO) for exoplanet imaging. We describe in this paper current and upcoming narrow field-of-view capabilities provided by the Subaru Extreme Adaptive Optics Adaptive Optics (SCExAO) system and its instrument modules, as well as the upcoming 3000-actuator upgrade of the Nasmyth AO system. △ Less

Submitted 2 August, 2022; originally announced August 2022.

Comments: 11 pages, to appear in SPIE Proceedings of Astronomical Telescopes + Instrumentation, 2022

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

XXL-HSC: The link between AGN activity and star formation in the Early Universe ($z\geqslant3.5$)

Authors: E. Pouliasis, G. Mountrichas, I. Georgantopoulos, A. Ruiz, R. Gilli, E. Koulouridis, M. Akiyama, Y. Ueda, C. Garrel, T. Nagao, S. Paltani, M. Pierre, Y. Toba, C. Vignali

Abstract: In this work, we aimed at investigating the star formation rate of active galactic nuclei host galaxies in the early Universe. To this end, we constructed a sample of 149 luminous ($\rm L_{2-10keV} > 10^{44}\,erg\,s^{-1}$) X-ray AGNs at $\rm z \geq3.5$ selected in three fields with different depths and observed areas (Chandra COSMOS Legacy survey, XMM-XXL North and eFEDS). We built their spectral… ▽ More In this work, we aimed at investigating the star formation rate of active galactic nuclei host galaxies in the early Universe. To this end, we constructed a sample of 149 luminous ($\rm L_{2-10keV} > 10^{44}\,erg\,s^{-1}$) X-ray AGNs at $\rm z \geq3.5$ selected in three fields with different depths and observed areas (Chandra COSMOS Legacy survey, XMM-XXL North and eFEDS). We built their spectral energy distributions (SED) using available multi-wavelength photometry from X-rays up to far-IR. Then, we estimated the stellar mass, M$_{*}$, and the SFR of the AGNs using the X-CIGALE SED fitting algorithm. After applying several quality criteria, we ended up with 89 high-z sources. More than half (55\%) of the X-ray sample have spectroscopic redshifts. Based on our analysis, our high-z X-ray AGNs live in galaxies with median $\rm M_{*}=5.6 \times10^{10}~M_\odot$ and $\rm SFR_{*}\approx240\,M_\odot yr^{-1}$. The majority of the high-z sources ($\sim89$\%) were found inside or above the main sequence (MS) of star-forming galaxies. Estimation of the normalised SFR, $\rm SFR_{NORM}$, defined as the ratio of the SFR of AGNs to the SFR of MS galaxies, showed that the SFR of AGNs is enhanched by a factor of $\sim 1.8$ compared to non-AGN star-forming systems. Combining our results with previous studies at lower redshifts, we confirmed that $\rm SFR_{NORM}$ does not evolve with redshift. Using the specific BHAR (i.e., $\rm L_X$ divided by $\rm M_{*}$), $\rm λ_{BHAR}$, that can be used as a tracer of the Eddington ratio, we found that the bulk of AGNs that lie inside or above the MS have higher specific accretion rates compared to sources below the MS. Finally, we found indications that the SFR of the most massive AGN host galaxies ($\rm log\,(M_{*}/ M_\odot) >10^{11.5-12}$) remains roughly constant as a function of M$_*$, in agreement with the SFR of MS star-forming galaxies. △ Less

Submitted 28 June, 2022; originally announced June 2022.

Comments: Accepted for publication in A&A. 19 pages, 12 figures (+6 in Appendix), 4 tables

Journal ref: A&A 667, A56 (2022)

Interpreting Graph-based Sybil Detection Methods as Low-Pass Filtering

Authors: Satoshi Furutani, Toshiki Shibahara, Mitsuaki Akiyama, Masaki Aida

Abstract: Online social networks (OSNs) are threatened by Sybil attacks, which create fake accounts (also called Sybils) on OSNs and use them for various malicious activities. Therefore, Sybil detection is a fundamental task for OSN security. Most existing Sybil detection methods are based on the graph structure of OSNs, and various methods have been proposed recently. However, although almost all methods h… ▽ More Online social networks (OSNs) are threatened by Sybil attacks, which create fake accounts (also called Sybils) on OSNs and use them for various malicious activities. Therefore, Sybil detection is a fundamental task for OSN security. Most existing Sybil detection methods are based on the graph structure of OSNs, and various methods have been proposed recently. However, although almost all methods have been compared experimentally in terms of detection performance and noise robustness, theoretical understanding of them is still lacking. In this study, we show that existing graph-based Sybil detection methods can be interpreted in a unified framework of low-pass filtering. This framework enables us to theoretically compare and analyze each method from two perspectives: filter kernel properties and the spectrum of shift matrices. Our analysis reveals that the detection performance of each method depends on how well low-pass filtering can extract low frequency components and remove noisy high frequency components. Furthermore, on the basis of the analysis, we propose a novel Sybil detection method called SybilHeat. Numerical experiments on synthetic graphs and real social networks demonstrate that SybilHeat performs consistently well on graphs with various structural properties. This study lays a theoretical foundation for graph-based Sybil detection and leads to a better understanding of Sybil detection methods. △ Less

Submitted 22 June, 2022; originally announced June 2022.

On the Connection between Supermassive Black Hole and Galaxy Growth in the Reionization Epoch

Authors: Junyao Li, John D. Silverman, Takuma Izumi, Wanqiu He, Masayuki Akiyama, Kohei Inayoshi, Yoshiki Matsuoka, Masafusa Onoue, Yoshiki Toba

Abstract: The correlation between the mass of supermassive black holes (SMBHs; $\mathcal{M}_{\rm BH}$) and their host galaxies ($\mathcal{M}_\star$) in the reionization epoch provides valuable constraints on their early growth. High-redshift quasars typically have a $\mathcal{M}_{\rm BH}$/$\mathcal{M}_\star$ ratio significantly elevated in comparison to the local value. However, the degree to which this app… ▽ More The correlation between the mass of supermassive black holes (SMBHs; $\mathcal{M}_{\rm BH}$) and their host galaxies ($\mathcal{M}_\star$) in the reionization epoch provides valuable constraints on their early growth. High-redshift quasars typically have a $\mathcal{M}_{\rm BH}$/$\mathcal{M}_\star$ ratio significantly elevated in comparison to the local value. However, the degree to which this apparent offset is driven by observational biases is unclear for the most distant quasars. To address this issue, we model the sample selection and measurement biases for a compilation of 20 quasars at $z\sim6$ with host properties based on ALMA observations. We find that the observed distribution of quasars in the $\mathcal{M}_{\rm BH} - \mathcal{M}_\star$ plane can be reproduced by assuming that the underlying SMBH population at $z\sim6$ follows the relationship in the local universe. However, a positive or even a negative evolution in $\mathcal{M}_{\rm BH}$/$\mathcal{M}_\star$ can also explain the data, depending on whether the intrinsic scatter evolves and the strength of various systematic uncertainties. To break these degeneracies, an improvement in the accuracy of mass measurements and an expansion of the current sample to lower $\mathcal{M}_{\rm BH}$ limits are needed. Furthermore, assuming a radiative efficiency of 0.1 and quasar duty cycles estimated from the active SMBH fraction, significant outliers in $\mathcal{M}_{\rm BH}$/$\mathcal{M}_\star$ tend to move toward the local relation given their instantaneous BH growth rate and star formation rate. This may provide evidence for a self-regulated SMBH-galaxy coevolution scenario that is in place at $z\sim6$, with AGN feedback being a possible driver. △ Less

Submitted 29 May, 2022; v1 submitted 20 March, 2022; originally announced March 2022.

Comments: 10 pages, 5 figures, published in ApJL

Multi-line assessment of narrow-line regions in $z \sim$ 3 radio galaxies

Authors: Koki Terao, Tohru Nagao, Kyoko Onishi, Kenta Matsuoka, Masayuki Akiyama, Yoshiki Matsuoka, Takuji Yamashita

Abstract: In this paper, we utilize high-quality rest-UV spectra of three radio galaxies at $z \sim$ 3 observed with VLT/FORS2 to measure the flux of several emission lines including relatively faint ones, such as NIV]$λ$1486, OIII]$λ$1663, and [NeIV]$λ$2424. Additionally, we collect fluxes of faint rest-UV emission lines in 12 $z \sim$ 3 radio galaxies from the literature. Previously, physical and chemical… ▽ More In this paper, we utilize high-quality rest-UV spectra of three radio galaxies at $z \sim$ 3 observed with VLT/FORS2 to measure the flux of several emission lines including relatively faint ones, such as NIV]$λ$1486, OIII]$λ$1663, and [NeIV]$λ$2424. Additionally, we collect fluxes of faint rest-UV emission lines in 12 $z \sim$ 3 radio galaxies from the literature. Previously, physical and chemical properties of narrow-line regions (NLRs) in high-$z$ active galactic nuclei (AGNs) have been investigated mostly by using only strong rest-UV emission-lines (e.g., NV$λ$1240, CIV$λ$1549, HeII$λ$1640, and CIII]$λ$1909). Such strong-line diagnostics are based on various assumptions due to the limitation in the number of available emission-line constraints. In this work, both physical and chemical properties of NLR clouds in each object are estimated by fitting detailed photoionization models to the measured emission-line fluxes. We confirm that the metallicity of NLRs in AGNs at $z \sim$ 3 is solar or super-solar, without assuming the gas density and ionization parameter thanks to the constrains from the faint emission lines. This result suggests that high-$z$ radio galaxies are already chemically matured at $z \sim$ 3. △ Less

Submitted 7 March, 2022; originally announced March 2022.

Comments: 17 pages, 10 figures, accepted for publication in the ApJ

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). XVI. 69 New Quasars at 5.8 < z < 7.0

Authors: Yoshiki Matsuoka, Kazushi Iwasawa, Masafusa Onoue, Takuma Izumi, Nobunari Kashikawa, Michael A. Strauss, Masatoshi Imanishi, Tohru Nagao, Masayuki Akiyama, John D. Silverman, Naoko Asami, James Bosch, Hisanori Furusawa, Tomotsugu Goto, James E. Gunn, Yuichi Harikane, Hiroyuki Ikeda, Rikako Ishimoto, Toshihiro Kawaguchi, Nanako Kato, Satoshi Kikuta, Kotaro Kohno, Yutaka Komiyama, Chien-Hsiu Lee, Robert H. Lupton , et al. (19 additional authors not shown)

Abstract: We present the spectroscopic discovery of 69 quasars at 5.8 < z < 7.0, drawn from the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) imaging survey data. This is the 16th publication from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, and completes identification of all but the faintest candidates (i.e., i-band dropouts with zAB < 24 and y-band detections, and z… ▽ More We present the spectroscopic discovery of 69 quasars at 5.8 < z < 7.0, drawn from the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) imaging survey data. This is the 16th publication from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, and completes identification of all but the faintest candidates (i.e., i-band dropouts with zAB < 24 and y-band detections, and z-band dropouts with yAB < 24) with Bayesian quasar probability Pq > 0.1 in the HSC-SSP third public data release (PDR3). The sample reported here also includes three quasars with Pq < 0.1 at z ~ 6.6, which we selected in an effort to completely cover the reddest point sources with simple color cuts. The number of high-z quasars discovered in SHELLQs has now grown to 162, including 23 type-II quasar candidates. This paper also presents identification of seven galaxies at 5.6 < z < 6.7, an [O III] emitter at z = 0.954, and 31 Galactic cool stars and brown dwarfs. High-z quasars and galaxies comprise 75 % and 16 % respectively of all the spectroscopic SHELLQs objects that pass our latest selection algorithm with the PDR3 photometry. That is, a total of 91 % of the objects lie at z > 5.6. This demonstrates that the algorithm has very high efficiency, even though we are probing an unprecedentedly low-luminosity population down to M1450 ~ -21 mag. △ Less

Submitted 24 November, 2021; originally announced November 2021.

Comments: Accepted for publication in The Astrophysical Journal Supplement Series

XXL-HSC: An updated catalogue of high-redshift ($z\geqslant3.5$) X-ray AGN in the XMM-XXL northern field. Constraints on the bright end of the soft logN-logS

Authors: E. Pouliasis, I. Georgantopoulos, A. Ruiz, R. Gilli, E. Koulouridis, M. Akiyama, Y. Ueda, L. Chiappetti, C. Garrel, C. Horellou, T. Nagao, S. Paltani, M. Pierre, Y. Toba, C. Vignali

Abstract: X-rays offer a reliable method to identify Active Galactic Nuclei (AGN). However, in the high-redshift Universe X-ray AGN are poorly sampled due to their relatively low space density and the small areas covered by X-ray surveys. In addition to wide-area X-ray surveys, it is important to have deep optical data in order to locate the optical counterparts and determine their redshifts. In this work,… ▽ More X-rays offer a reliable method to identify Active Galactic Nuclei (AGN). However, in the high-redshift Universe X-ray AGN are poorly sampled due to their relatively low space density and the small areas covered by X-ray surveys. In addition to wide-area X-ray surveys, it is important to have deep optical data in order to locate the optical counterparts and determine their redshifts. In this work, we build a high-redshift ($z\geqslant3.5$) X-ray selected AGN sample in the XMM-XXL northern field using the most updated [0.5-2 keV] catalogue along with a plethora of new spectroscopic and multi-wavelength catalogues, including the deep optical Subaru Hyper Suprime-Cam (HSC) data reaching magnitude limits i$\sim$26 mag. We select all the spectroscopically confirmed AGN and complement this sample with high-redshift candidates that are HSC g- and r-band dropouts. To confirm the dropouts, we derive their photometric redshifts using spectral energy distribution techniques. We end up with a sample of 54 high-z sources (28 with spec-z), the largest in this field so far (almost three times larger than in previous studies), and we estimate the possible contamination and completeness. We calculate the number counts (logN-logS) in different redshift bins and compare our results with previous studies and models. We provide the strongest high-redshift AGN constraints yet at bright fluxes ($\rm f_{0.5-2~keV} > 10^{-15} erg~s^{-1}~cm^{-2}$). The samples of $z\geqslant3.5$, $z\geqslant4$ and $z\geqslant5$ are in agreement with an exponential decline model similar to that witnessed at optical wavelengths. Our work emphasizes the importance of using wide area X-ray surveys with deep optical data to uncover high redshift AGN. △ Less

Submitted 8 November, 2021; originally announced November 2021.

Comments: Accepted for publication in A&A. 16 pages, 12 figures, 3 tables

Journal ref: A&A 658, A175 (2022)

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

Ionizing radiation from AGNs at z>3.3 with the Subaru Hyper Suprime-Cam Survey and the CFHT Large Area U-band Deep Survey (CLAUDS)

Authors: Ikuru Iwata, Marcin Sawicki, Akio K. Inoue, Masayuki Akiyama, Genoveva Micheva, Toshihiro Kawaguchi, Nobunari Kashikawa, Stephen Gwyn, Stephane Arnouts, Jean Coupon, Guillaume Desprez

Abstract: We use deep and wide imaging data from the CFHT Large Area U-band Deep Survey (CLAUDS) and the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) to constrain the ionizing radiation (Lyman Continuum; LyC) escape fraction from AGNs at $z \sim 3 - 4$. For 94 AGNs with spectroscopic redshifts at $3.3 < z < 4.0$, we use their U-band / i-band flux ratios to estimate LyC transmission of individual AGN… ▽ More We use deep and wide imaging data from the CFHT Large Area U-band Deep Survey (CLAUDS) and the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) to constrain the ionizing radiation (Lyman Continuum; LyC) escape fraction from AGNs at $z \sim 3 - 4$. For 94 AGNs with spectroscopic redshifts at $3.3 < z < 4.0$, we use their U-band / i-band flux ratios to estimate LyC transmission of individual AGNs. The distribution of their LyC transmission shows values lower than the range of LyC transmission values for IGM of the same redshift range, which suggests that LyC escape fraction of AGNs at $z>3.3$ is considerably lower than unity in most cases. We do not find any trend in LyC transmission values depending on their UV luminosities. Based on the photometry of stacked images we find the average flux ratio of LyC and non-ionizing UV photons escaping from the objects $(f_{LyC}/f_{UV})^{out} = 0.182 \pm 0.043$ for AGNs at $3.3<z<3.6$, which corresponds to LyC escape fraction $f_{esc} = 0.303 \pm 0.072$ if we assume a fiducial intrinsic SED of AGN. Based on the estimated LyC escape fraction and the UV luminosity function of AGNs, we argue that UV-selected AGNs' contribution to the LyC emissivity at the epoch is minor, although the size of their contribution largely depends on the shape of the UV luminosity function. △ Less

Submitted 21 September, 2021; originally announced September 2021.

Comments: 19 pages, 12 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

A Wide and Deep Exploration of Radio Galaxies with Subaru HSC (WERGS). IV. Rapidly Growing (Super-)Massive Black Holes in Extremely Radio-Loud Galaxies

Authors: Kohei Ichikawa, Takuji Yamashita, Yoshiki Toba, Tohru Nagao, Kohei Inayoshi, Maria Charisi, Wanqiu He, Alexander Y. Wagner, Masayuki Akiyama, Bovornpratch Vijarnwannaluk, Xaioyang Chen, Masaru Kajisawa, Taiki Kawamuro, Chien-Hsiu Lee, Yoshiki Matsuoka, Malte Schramm, Hyewon Suh, Masayuki Tanaka, Hisakazu Uchiyama, Yoshihiro Ueda, Janek Pflugradt, Hikaru Fukuchi

Abstract: We present the optical and infrared properties of 39 extremely radio-loud galaxies discovered by cross-matching the Subaru/Hyper Suprime-Cam (HSC) deep optical imaging survey and VLA/FIRST 1.4 GHz radio survey. The recent Subaru/HSC strategic survey revealed optically-faint radio galaxies (RG) down to $g_\mathrm{AB} \sim 26$, opening a new parameter space of extremely radio-loud galaxies (ERGs) wi… ▽ More We present the optical and infrared properties of 39 extremely radio-loud galaxies discovered by cross-matching the Subaru/Hyper Suprime-Cam (HSC) deep optical imaging survey and VLA/FIRST 1.4 GHz radio survey. The recent Subaru/HSC strategic survey revealed optically-faint radio galaxies (RG) down to $g_\mathrm{AB} \sim 26$, opening a new parameter space of extremely radio-loud galaxies (ERGs) with radio-loudness parameter of $\log \mathcal{R}_\mathrm{rest} = \log (f_{1.4 \mathrm{GHz,rest}}/f_{g,\mathrm{rest}}) >4$. Because of their optical faintness and small number density of $\sim1~$deg$^{-2}$, such ERGs were difficult to find in the previous wide but shallow, or deep but small area optical surveys. ERGs show intriguing properties that are different from the conventional RGs: (1) most ERGs reside above or on the star-forming main-sequence, and some of them might be low-mass galaxies with $\log (M_\star/M_\odot) < 10$. (2) ERGs exhibit a high specific black hole accretion rate, reaching the order of the Eddington limit. The intrinsic radio-loudness ($\mathcal{R}_\mathrm{int}$), defined by the ratio of jet power over bolometric radiation luminosity, is one order of magnitude higher than that of radio quasars. This suggests that ERGs harbor a unique type of active galactic nuclei (AGN) that show both powerful radiations and jets. Therefore, ERGs are prominent candidates of very rapidly growing black holes reaching Eddington-limited accretion just before the onset of intensive AGN feedback. △ Less

Submitted 5 August, 2021; originally announced August 2021.

Comments: 20 pages, 10 figures, accepted for publication in ApJ

The eROSITA Final Equatorial-Depth Survey (eFEDS): A multiwavelength view of WISE mid-infrared galaxies/active galactic nuclei

Authors: Yoshiki Toba, Teng Liu, Tanya Urrutia, Mara Salvato, Junyao Li, Yoshihiro Ueda, Marcella Brusa, Naomichi Yutani, Keiichi Wada, Atsushi J. Nishizawa, Johannes Buchner, Tohru Nagao, Andrea Merloni, Masayuki Akiyama, Riccardo Arcodia, Bau-Ching Hsieh, Kohei Ichikawa, Masatoshi Imanishi, Kaiki T. Inoue, Toshihiro Kawaguchi, Georg Lamer, Kirpal Nandra, John D. Silverman, Yuichi Terashima

Abstract: We investigate the physical properties--such as the stellar mass, SFR, IR luminosity, X-ray luminosity, and hydrogen column density--of MIR galaxies and AGN at $z < 4$ in the 140 deg$^2$ field observed by SRG/eROSITA through the eFEDS survey. By cross-matching the WISE 22 $μ$m (W4)-detected sample and the eFEDS X-ray point-source catalog, we find that 692 extragalactic objects are detected by eROS… ▽ More We investigate the physical properties--such as the stellar mass, SFR, IR luminosity, X-ray luminosity, and hydrogen column density--of MIR galaxies and AGN at $z < 4$ in the 140 deg$^2$ field observed by SRG/eROSITA through the eFEDS survey. By cross-matching the WISE 22 $μ$m (W4)-detected sample and the eFEDS X-ray point-source catalog, we find that 692 extragalactic objects are detected by eROSITA. We have compiled a multiwavelength dataset. We have also performed (i) an X-ray spectral analysis, (ii) SED fitting using X-CIGALE, (iii) 2D image-decomposition analysis using Subaru HSC images, and (iv) optical spectral fitting with QSFit to investigate the AGN and host-galaxy properties. For 7,088 WISE W4 objects that are undetected by eROSITA, we have performed an X-ray stacking analysis to examine the typical physical properties of these X-ray faint and/or probably obscured objects. We find that (i) 82% of the eFEDS-W4 sources are classified as X-ray AGN with $\log\,L_{\rm X} >$ 42 erg s$^{-1}$; (ii) 67% and 24% of the objects have $\log\,(L_{\rm IR}/L_{\odot}) > 12$ and 13, respectively; (iii) the relationship between $L_{\rm X}$ and the 6 $μ$m luminosity is consistent with that reported in previous works; and (iv) the relationship between the Eddington ratio and $N_{\rm H}$ for the eFEDS-W4 sample and a comparison with a model prediction from a galaxy-merger simulation indicates that approximately 5% of the eFEDS-W4 sources in our sample are likely to be in an AGN-feedback phase, in which strong radiation pressure from the AGN blows out the surrounding material from the nuclear region. Thanks to the wide area coverage of eFEDS, we have been able to constrain the ranges of the physical properties of the WISE W4 sample of AGNs at $z < 4$, providing a benchmark for forthcoming studies on a complete census of MIR galaxies selected from the full-depth eROSITA all-sky survey. △ Less

Submitted 24 August, 2021; v1 submitted 28 June, 2021; originally announced June 2021.

Comments: 18 pages, 19 figures, and 3 tables, accepted to appear on A&A, Special Issue: The Early Data Release of eROSITA and Mikhail Pavlinsky ART-XC on the SRG Mission

Journal ref: A&A 661, A15 (2022)

The eROSITA Final Equatorial-Depth Survey (eFEDS): An X-ray bright, extremely luminous infrared galaxy at z = 1.87

Authors: Yoshiki Toba, Marcella Brusa, Teng Liu, Johannes Buchner, Yuichi Terashima, Tanya Urrutia, Mara Salvato, Masayuki Akiyama, Riccardo Arcodia, Andy D. Goulding, Yuichi Higuchi, Kaiki T. Inoue, Toshihiro Kawaguchi, Georg Lamer, Andrea Merloni, Tohru Nagao, Yoshihiro Ueda, Kirpal Nandra

Abstract: In this study, we investigate the X-ray properties of WISE J090924.01+000211.1 (WISEJ0909+0002), an extremely luminous infrared (IR) galaxy (ELIRG) at $z_{\rm spec}$= 1.871 in the eROSITA final equatorial depth survey (eFEDS). WISEJ0909+0002 is a WISE 22 $μ$m source, located in the GAMA-09 field, detected by eROSITA during the performance and verification phase. The corresponding optical spectrum… ▽ More In this study, we investigate the X-ray properties of WISE J090924.01+000211.1 (WISEJ0909+0002), an extremely luminous infrared (IR) galaxy (ELIRG) at $z_{\rm spec}$= 1.871 in the eROSITA final equatorial depth survey (eFEDS). WISEJ0909+0002 is a WISE 22 $μ$m source, located in the GAMA-09 field, detected by eROSITA during the performance and verification phase. The corresponding optical spectrum indicates that this object is a type-1 active galactic nucleus (AGN). Observations from eROSITA combined with Chandra and XMM-Newton archival data indicate a very luminous ($L$ (2--10 keV) = ($2.1 \pm 0.2) \times 10^{45}$ erg s$^{-1}$) unobscured AGN with a power-law photon index of $Γ$ = 1.73$_{-0.15}^{+0.16}$, and an absorption hydrogen column density of $\log\,(N_{\rm H}/{\rm cm}^{-2}) < 21.0$. The IR luminosity was estimated to be $L_{\rm IR}$ = (1.79 $\pm$ 0.09) $\times 10^{14}\, L_{\odot}$ from spectral energy distribution modeling based on 22 photometric data (X-ray to far-IR) with X-CIGALE, which confirmed that WISEJ0909+0002 is an ELIRG. A remarkably high $L_{\rm IR}$ despite very low $N_{\rm H}$ would indicate that we are witnessing a short-lived phase in which hydrogen gas along the line of sight is blown outwards, whereas warm and hot dust heated by AGNs still exist. As a consequence of eROSITA all-sky survey, $6.8_{-5.6}^{+16}\times 10^2$ such X-ray bright ELIRGs are expected to be discovered in the entire extragalactic sky ($|b| > 10^\circ$). This can potentially be the key population to constrain the bright-end of IR luminosity functions. △ Less

Submitted 6 April, 2021; originally announced April 2021.

Comments: 10 pages, 5 figures, and 3 tables, accepted for publication in A&A Letters (special Issue: First science highlights from SRG/eROSITA)

Journal ref: A&A 649, L11 (2021)

A First Look at COVID-19 Domain Names: Origin and Implications

Authors: Ryo Kawaoka, Daiki Chiba, Takuya Watanabe, Mitsuaki Akiyama, Tatsuya Mori

Abstract: This work takes a first look at domain names related to COVID-19 (Cov19doms in short), using a large-scale registered Internet domain name database, which accounts for 260M of distinct domain names registered for 1.6K of distinct top-level domains. We extracted 167K of Cov19doms that have been registered between the end of December 2019 and the end of September 2020. We attempt to answer the follo… ▽ More This work takes a first look at domain names related to COVID-19 (Cov19doms in short), using a large-scale registered Internet domain name database, which accounts for 260M of distinct domain names registered for 1.6K of distinct top-level domains. We extracted 167K of Cov19doms that have been registered between the end of December 2019 and the end of September 2020. We attempt to answer the following research questions through our measurement study: RQ1: Is the number of Cov19doms registrations correlated with the COVID-19 outbreaks?, RQ2: For what purpose do people register Cov19doms? Our chief findings are as follows: (1) Similar to the global COVID-19 pandemic observed around April 2020, the number of Cov19doms registrations also experienced the drastic growth, which, interestingly, pre-ceded the COVID-19 pandemic by about a month, (2) 70 % of active Cov19doms websites with visible content provided useful information such as health, tools, or product sales related to COVID-19, and (3) non-negligible number of registered Cov19doms was used for malicious purposes. These findings imply that it has become more challenging to distinguish domain names registered for legitimate purposes from others and that it is crucial to pay close attention to how Cov19doms will be used/misused in the future. △ Less

Submitted 10 February, 2021; originally announced February 2021.

Comments: 9 pages, 4 figures, 4 tables. Accepted at the Passive and Active Measurement Conference 2021 (PAM 2021)

Black Hole and Galaxy Coevolution in Moderately Luminous Active Galactic Nuclei at z~1.4 in SXDF

Authors: Kenta Setoguchi, Yoshihiro Ueda, Yoshiki Toba, Masayuki Akiyama

Abstract: We investigate the relation of black hole mass versus host stellar mass and that of mass accretion rate versus star formation rate (SFR) in moderately luminous ($\log L_{\rm bol} \sim 44.5-46.5\ {\rm erg\ s^{-1}}$), X-ray selected broad-line active galactic nuclei (AGNs) at $z=1.18-1.68$ in the Subaru/XMM-Newton Deep Field. The far-infrared to far-ultraviolet spectral energy distributions of 85 AG… ▽ More We investigate the relation of black hole mass versus host stellar mass and that of mass accretion rate versus star formation rate (SFR) in moderately luminous ($\log L_{\rm bol} \sim 44.5-46.5\ {\rm erg\ s^{-1}}$), X-ray selected broad-line active galactic nuclei (AGNs) at $z=1.18-1.68$ in the Subaru/XMM-Newton Deep Field. The far-infrared to far-ultraviolet spectral energy distributions of 85 AGNs are reproduced with the latest version of Code Investigating GALaxy Emission ({\tt CIGALE}), where the AGN clumpy torus model {\tt SKIRTOR} is implemented. Most of their hosts are confirmed to be main-sequence star-forming galaxies. We find that the mean ratio of the black hole mass ($M_{\rm BH}$) to the total stellar mass ($M_{\rm stellar}$) is $\log M_{\rm BH}/M_{\rm stellar} = -2.2$, which is similar to the local black hole-to-bulge mass ratio. This suggests that if the host galaxies of these moderately luminous AGNs at $z\sim1.4$ are dominated by bulges, they already established the local black hole mass-bulge mass relation; if they are disk dominant, their black holes are overmassive relative to the bulges. AGN bolometric luminosities and SFR show a good correlation with ratios higher than that expected from the local black hole-to-bulge mass relation, suggesting that these AGNs are in a SMBH-growth dominant phase. △ Less

Submitted 11 March, 2021; v1 submitted 27 January, 2021; originally announced January 2021.

Comments: 9 pages, 5 figures, 1 table, accepted for publication in ApJ; the error in the figure corrected, sentences changed a little

Introducing piXedfit -- a Spectral Energy Distribution Fitting Code Designed for Resolved Sources

Authors: Abdurro'uf, Yen-Ting Lin, Po-Feng Wu, Masayuki Akiyama

Abstract: We present piXedfit, pixelized spectral energy distribution (SED) fitting, a Python package that provides tools for analyzing spatially resolved properties of galaxies using multiband imaging data alone or in combination with integral field spectroscopy (IFS) data. piXedfit has six modules that can handle all tasks in the spatially resolved SED fitting. The SED fitting module uses the Bayesian inf… ▽ More We present piXedfit, pixelized spectral energy distribution (SED) fitting, a Python package that provides tools for analyzing spatially resolved properties of galaxies using multiband imaging data alone or in combination with integral field spectroscopy (IFS) data. piXedfit has six modules that can handle all tasks in the spatially resolved SED fitting. The SED fitting module uses the Bayesian inference technique with two kinds of posteriors sampling methods: Markov Chain Monte Carlo (MCMC) and random densely-sampling of parameter space (RDSPS). We test the performance of the SED fitting module using mock SEDs of simulated galaxies from IllustrisTNG. The SED fitting with both posteriors sampling methods can recover physical properties and star formation histories of the IllustrisTNG galaxies well. We further test the performance of piXedfit modules by analyzing 20 galaxies observed by the CALIFA and MaNGA surveys. The data comprises of 12-band imaging data from GALEX, SDSS, 2MASS, and WISE, and the IFS data from CALIFA or MaNGA. piXedfit can spatially match (in resolution and sampling) of the imaging and IFS data. By fitting only the photometric SEDs, piXedfit can predict the spectral continuum, $\text{D}_{\rm n}4000$, $H_α$, and $H_β$ well. The star formation rate (SFR) derived by piXedfit is consistent with that derived from $H_α$ emission. The RDSPS method gives equally good fitting results as the MCMC and it is much faster than the MCMC. piXedfit is a versatile tool equipped with a parallel computing module for efficient analysis of large datasets, and will be made publicly available (https://github.com/aabdurrouf/piXedfit). △ Less

Submitted 2 March, 2021; v1 submitted 24 January, 2021; originally announced January 2021.

Comments: 40 pages, 21 figures, 1 table, Accepted for publication in The Astrophysical Journal Supplement Series (ApJS). piXedfit will be available at https://github.com/aabdurrouf/piXedfit and documented at https://pixedfit.readthedocs.io/en/latest/index.html

Atmospheric turbulence profiling with multi-aperture scintillation of a Shack-Hartmann sensor

Authors: Hajime Ogane, Masayuki Akiyama, Shin Oya, Yoshito Ono

Abstract: Adaptive optics (AO) systems using tomographic estimation of three-dimensional structure of atmospheric turbulence requires vertical atmospheric turbulence profile, which describes turbulence strength as a function of altitude as a prior information. We propose a novel method to reconstruct the profile by applying Multi Aperture Scintillation Sensor (MASS) method to scintillation data obtained by… ▽ More Adaptive optics (AO) systems using tomographic estimation of three-dimensional structure of atmospheric turbulence requires vertical atmospheric turbulence profile, which describes turbulence strength as a function of altitude as a prior information. We propose a novel method to reconstruct the profile by applying Multi Aperture Scintillation Sensor (MASS) method to scintillation data obtained by a Shack-Hartmann wavefront sensor (SH-WFS). Compared to the traditional MASS, which uses atmospheric scintillation within 4 concentric annular apertures, the new method utilizes scintillation in several hundreds of spatial patterns, which are created by combinations of SH-WFS subapertures. Accuracy of the turbulence profile reconstruction is evaluated with Bayesian inference, and it is confirmed that turbulence profile with more than 10 layers can be reconstructed thanks to the large number of constraints. We demonstrate the new method with a SH-WFS attached to the 50 cm telescope at Tohoku university and confirm that general characteristics of atmospheric turbulence profile is reproduced. △ Less

Submitted 11 January, 2021; originally announced January 2021.

Comments: 11 pages, 14 figures. Accepted for publication in MNRAS

Faint Quasars Live in the Same Number Density Environments as Lyman Break Galaxies at z ~ 4

Authors: Hisakazu Uchiyama, Masayuki Akiyama, Jun Toshikawa, Nobunari Kashikawa, Roderik Overzier, Tohru Nagao, Kohei Ichikawa, Murilo Marinello, Masatoshi Imanishi, Masayuki Tanaka, Yoshiki Matsuoka, Yutaka Komiyama, Shogo Ishikawa, Masafusa Onoue, Mariko Kubo, Yuichi Harikane, Kei Ito, Shigeru Namiki, Yongming Liang

Abstract: Characterizing high-z quasar environments is key to understanding the co-evolution of quasars and the surrounding galaxies. To restrict their global picture, we statistically examine the g-dropout galaxy overdensity distribution around 570 faint quasar candidates at z ~ 4, based on the Hyper Suprime-Cam Subaru Strategic Program survey. We compare the overdensity significances of g-dropout galaxies… ▽ More Characterizing high-z quasar environments is key to understanding the co-evolution of quasars and the surrounding galaxies. To restrict their global picture, we statistically examine the g-dropout galaxy overdensity distribution around 570 faint quasar candidates at z ~ 4, based on the Hyper Suprime-Cam Subaru Strategic Program survey. We compare the overdensity significances of g-dropout galaxies around the quasars with those around g-dropout galaxies, and find no significant difference between their distributions. A total of 4 (22) out of the 570 faint quasars, 0.7_{-0.4}^{+0.4} (3.9_{-0.8}^{+0.8}) %, are found to be associated with the > 4 sigma overdense regions within an angular separation of 1.8 (3.0) arcmin, which is the typical size of protoclusters at this epoch. This is similar to the fraction of g-dropout galaxies associated with the > 4 sigma overdense regions. This result is consistent with our previous work that 1.3_{-0.9}^{+0.9} % and 2.0_{-1.1}^{+1.1} % of luminous quasars detected in the Sloan Digital Sky Survey exist in the > 4 sigma overdense regions within 1.8 and 3.0 arcmin separations, respectively. Therefore, we suggest that the galaxy number densities around quasars are independent of their luminosity, and most quasars do not preferentially appear in the richest protocluster regions at z ~ 4. The lack of an apparent positive correlation between the quasars and the protoclusters implies that: i) the gas-rich major merger rate is relatively low in the protocluster regions, ii) most high-z quasars may appear through secular processes, or iii) some dust-obscured quasars exist in the protocluster regions. △ Less

Submitted 19 November, 2020; originally announced November 2020.

Comments: 12 pages, 4 figures

NuSTAR Non-detection of a Faint Active Galactic Nucleus in an Ultraluminous IR Galaxy with Kpc-scale Fast Wind

Authors: Xiaoyang Chen, Kohei Ichikawa, Hirofumi Noda, Taiki Kawamuro, Toshihiro Kawaguchi, Yoshiki Toba, Masayuki Akiyama

Abstract: Large-scale outflows are generally considered as a possible evidence that active galactic nuclei (AGNs) can severely affect their host galaxies. Recently an ultraluminous IR galaxy (ULIRG) at $z=0.49$, AKARI J0916248+073034, was found to have a galaxy-scale [OIII] $λ$5007 outflow with one of the highest energy-ejection rates at $z<1.6$. However, the central AGN activity estimated from its torus mi… ▽ More Large-scale outflows are generally considered as a possible evidence that active galactic nuclei (AGNs) can severely affect their host galaxies. Recently an ultraluminous IR galaxy (ULIRG) at $z=0.49$, AKARI J0916248+073034, was found to have a galaxy-scale [OIII] $λ$5007 outflow with one of the highest energy-ejection rates at $z<1.6$. However, the central AGN activity estimated from its torus mid-IR (MIR) radiation is weak relative to the luminous [OIII] emission. In this work we report the first NuSTAR hard X-ray follow-up of this ULIRG to constrain its current AGN luminosity. The intrinsic 2-10 keV luminosity shows a 90% upper-limit of $3.0\times10^{43}$ erg s$^{-1}$ assuming Compton-thick obscuration ($N_{\rm H}=1.5\times10^{24}$ cm$^{-2}$), which is only 3.6% of the luminosity expected from the extinction corrected [OIII] luminosity. With the NuSTAR observation, we succeed to identify that this ULIRG has a most extreme case of X-ray deficit among local ULIRGs. A possible scenario to explain the drastic declining in both of the corona (X-ray) and torus (MIR) is that the primary radiation from the AGN accretion disk is currently in a fading status, as a consequence of a powerful nuclear wind suggested by its powerful ionized outflow in the galaxy scale. △ Less

Submitted 13 November, 2020; originally announced November 2020.

Comments: 8 pages, 2 figures, accepted for publication in The Astrophysical Journal Letters

The Faint End of the Quasar Luminosity Function at $z \sim 5$ from the Subaru Hyper Suprime-Cam Survey

Authors: Mana Niida, Tohru Nagao, Hiroyuki Ikeda, Masayuki Akiyama, Yoshiki Matsuoka, Wanqiu He, Kenta Matsuoka, Yoshiki Toba, Masafusa Onoue, Masakazu A. R. Kobayashi, Yoshiaki Taniguchi, Hisanori Furusawa, Yuichi Harikane, Masatoshi Imanishi, Nobunari Kashikawa, Toshihiro Kawaguchi, Yutaka Komiyama, Hikari Shirakata, Yuichi Terashima, Yoshihiro Ueda

Abstract: We present the quasar luminosity function at $z \sim 5$ derived from the optical wide-field survey data obtained as a part of the Subaru strategic program (SSP) with Hyper Suprime-Cam (HSC). From $\sim$81.8 deg$^2$ area in the Wide layer of the HSC-SSP survey, we selected 224 candidates of low-luminosity quasars at $z \sim 5$ by adopting the Lyman-break method down to $i = 24.1$ mag. Based on our… ▽ More We present the quasar luminosity function at $z \sim 5$ derived from the optical wide-field survey data obtained as a part of the Subaru strategic program (SSP) with Hyper Suprime-Cam (HSC). From $\sim$81.8 deg$^2$ area in the Wide layer of the HSC-SSP survey, we selected 224 candidates of low-luminosity quasars at $z \sim 5$ by adopting the Lyman-break method down to $i = 24.1$ mag. Based on our candidates and spectroscopically-confirmed quasars from the Sloan Digital Sky Survey (SDSS), we derived the quasar luminosity function at $z \sim 5$ covering a wide luminosity range of $-28.76 < M_{\rm 1450} < -22.32$ mag. We found that the quasar luminosity function is fitted by a double power-law model with a break magnitude of $M^{*}_{1450} = -25.05^{+0.10}_{-0.24}$ mag. The inferred number density of low-luminosity quasars is lower, and the derived faint-end slope, $-1.22^{+0.03}_{-0.10}$, is flatter than those of previous studies at $z \sim 5$. A compilation of the quasar luminosity function at $4 \leq z \leq 6$ from the HSC-SSP suggests that there is little redshift evolution in the break magnitude and in the faint-end slope within this redshift range, although previous studies suggest that the faint-end slope becomes steeper at higher redshifts. The number density of low-luminosity quasars decreases more rapidly from $z \sim 5$ to $z \sim 6$ than from $z \sim 4$ to $z \sim 5$. △ Less

Submitted 1 October, 2020; originally announced October 2020.

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

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). IX. Identification of Two Red Quasars at z > 5.6

Authors: Nanako Kato, Yoshiki Matsuoka, Masafusa Onoue, Shuhei Koyama, Yoshiki Toba, Masayuki Akiyama, Seiji Fujimoto, Masatoshi Imanishi, Kazushi Iwasawa, Takuma Izumi, Nobunari Kashikawa, Toshihiro Kawaguchi, Chien-Hsiu Lee, Takeo Minezaki, Tohru Nagao, Akatoki Noboriguchi, Michael A. Strauss

Abstract: We present the first discovery of dust-reddened quasars (red quasars) in the high-z universe (z >5.6). This is a result from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which is based on the sensitive multi-band optical imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. We identified four red quasar candidates from the spectroscopic… ▽ More We present the first discovery of dust-reddened quasars (red quasars) in the high-z universe (z >5.6). This is a result from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which is based on the sensitive multi-band optical imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. We identified four red quasar candidates from the spectroscopically confirmed 93 high-z quasars in the SHELLQs sample, based on detections in the Wide-field Infrared Survey Explorer (WISE) data at 3.4 and 4.6 um (rest-frame ~5000-6500 A). The amount of dust reddening was estimated with spectral energy distribution (SED) fits over optical and mid-infrared wavelengths. Two of the four candidates were found to be red quasars with dust reddening of E(B-V) > 0.1. The remaining SHELLQs quasars without individual WISE detections are significantly fainter in the WISE bands and bluer than the red quasars, although we did detect them in the W1 band in a stacked image. We also conducted the same SED fits for high-z optically-luminous quasars, but no red quasar was found. This demonstrates the power of Subaru HSC to discover high-z red quasars, which are fainter than the limiting magnitudes of past surveys in the rest-frame ultraviolet, due to dust extinction. △ Less

Submitted 28 July, 2020; v1 submitted 16 July, 2020; originally announced July 2020.

Comments: Accepted for publication in PASJ

Wide and Deep Exploration of Radio Galaxies with Subaru HSC (WERGS). III. Discovery of a z = 4.72 Radio Galaxy with Lyman Break Technique

Authors: Takuji Yamashita, Tohru Nagao, Hiroyuki Ikeda, Yoshiki Toba, Masaru Kajisawa, Yoshiaki Ono, Masayuki Tanaka, Masayuki Akiyama, Yuichi Harikane, Kohei Ichikawa, Toshihiro Kawaguchi, Taiki Kawamuro, Kotaro Kohno, Chien-Hsiu Lee, Kianhong Lee, Yoshiki Matsuoka, Mana Niida, Kazuyuki Ogura, Masafusa Onoue, Hisakazu Uchiyama

Abstract: We report a discovery of $z = 4.72$ radio galaxy, HSC J083913.17+011308.1, by using the Lyman break technique with the Hyper Suprime-Cam Subaru Strategic Survey (HSC-SSP) catalog for VLA FIRST radio sources. The number of known high-$z$ radio galaxies (HzRGs) at $z > 3$ is quite small to constrain the evolution of HzRGs so far. The deep and wide-area optical survey by HSC-SSP enables us to apply t… ▽ More We report a discovery of $z = 4.72$ radio galaxy, HSC J083913.17+011308.1, by using the Lyman break technique with the Hyper Suprime-Cam Subaru Strategic Survey (HSC-SSP) catalog for VLA FIRST radio sources. The number of known high-$z$ radio galaxies (HzRGs) at $z > 3$ is quite small to constrain the evolution of HzRGs so far. The deep and wide-area optical survey by HSC-SSP enables us to apply the Lyman break technique to a large search for HzRGs. For an HzRG candidate among pre-selected $r$-band dropouts with a radio detection, a follow-up optical spectroscopy with GMOS/Gemini has been performed. The obtained spectrum presents a clear Ly$α$ emission line redshifted to $z=4.72$. The SED fitting analysis with the rest-frame UV and optical photometries suggests the massive nature of this HzRG with $\log{M_*/M_{\odot}} = 11.4$. The small equivalent width of Ly$α$ and the moderately red UV colors indicate its dusty host galaxy, implying a chemically evolved and dusty system. The radio spectral index does not meet a criterion for an ultra-steep spectrum: $α^{325}_{1400}$ of $-1.1$ and $α^{150}_{1400}$ of $-0.9$, demonstrating that the HSC-SSP survey compensates for a sub-population of HzRGs which are missed in surveys focusing on an ultra-steep spectral index. △ Less

Submitted 8 June, 2020; originally announced June 2020.

Comments: 10 pages, 5 figures, accepted for publication in AJ

Properties of environment around AGN and luminous galaxy pairs through HSC wide survey

Authors: Yuji Shirasaki, Masayuki Akiyama, Yoshiki Toba, Wanqiu He, Tomotsugu Goto

Abstract: We investigated the properties of AGN environments, particularly environments where the association of luminous galaxies (LGs) is found within 4~Mpc from AGNs with redshifts of 0.8 -- 1.1. For comparison, three additional AGN environments, (namely, AGNs of all types, type~1 AGNs with X-ray and/or radio detection, and type~2 AGNs) and an environment of blue $M_{*}$, characteristic luminosity of the… ▽ More We investigated the properties of AGN environments, particularly environments where the association of luminous galaxies (LGs) is found within 4~Mpc from AGNs with redshifts of 0.8 -- 1.1. For comparison, three additional AGN environments, (namely, AGNs of all types, type~1 AGNs with X-ray and/or radio detection, and type~2 AGNs) and an environment of blue $M_{*}$, characteristic luminosity of the Schechter function, galaxies were investigated. The cross-correlation function with the surrounding galaxies was measured and compared between the AGN and blue galaxy samples. We also compared the distributions of color, absolute magnitude, and stellar mass of the galaxies around such target objects. The properties of clusters detected using surrounding galaxies selected based on a photometric redshift were examined and compared for different samples. The target AGNs were drawn from the Million Quasars (MILLIQUAS) catalog, and the blue galaxies were drawn from six redshift survey catalogs (SDSS, WiggleZ, DEEP2, VVDS, VIPERS, and PRIMUS). The galaxies used as a measure of the environment around the targets are drawn from S18a internal data released by the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). We found that, among the five AGN and blue galaxy samples considered, the environment of AGN-LG pairs is the most enriched with luminous galaxies. We also found an enhancement in the number of mass-selected clusters in the AGN-LG pair sample against those in the other samples. The results obtained in this study indicate that existence of multiple clusters is the major driver in the association of AGNs and LGs, rather than a single large-mass dark matter halo hosting the AGN. △ Less

Submitted 13 April, 2020; originally announced April 2020.

Comments: 46 pages, 21 figures, accepted for publication in PASJ

Tracing the co-evolution path of super massive black holes and spheroids with AKARI-selected ultra-luminous IR galaxies at intermediate redshifts

Authors: Xiao-Yang Chen, Masayuki Akiyama, Kohei Ichikawa, Hirofumi Noda, Yoshiki Toba, Issei Yamamura, Toshihiro Kawaguchi, Abdurro'uf, Mitsuru Kokubo

Abstract: We present the stellar population and ionized-gas outflow properties of ultra-luminous IR galaxies (ULIRGs) at $z=$ 0.1-1.0, which are selected from AKARI FIR all-sky survey. We construct a catalog of 1077 ULIRGs to examine feedback effect after major mergers. 202 out of the 1077 ULIRGs are spectroscopically identified by SDSS and Subaru/FOCAS observations. Thanks to deeper depth and higher resolu… ▽ More We present the stellar population and ionized-gas outflow properties of ultra-luminous IR galaxies (ULIRGs) at $z=$ 0.1-1.0, which are selected from AKARI FIR all-sky survey. We construct a catalog of 1077 ULIRGs to examine feedback effect after major mergers. 202 out of the 1077 ULIRGs are spectroscopically identified by SDSS and Subaru/FOCAS observations. Thanks to deeper depth and higher resolution of AKARI compared to the previous IRAS survey, and reliable identification from WISE MIR pointing, the sample is unique in identifying optically-faint (i$\sim$20) IR-bright galaxies, which could be missed in previous surveys. A self-consistent spectrum-SED decomposition method, which constrains stellar population properties in SED modeling based on spectral fitting results, has been employed for 149 ULIRGs whose optical continua are dominated by host galaxies. They are massive galaxies ($M_{\rm star}\sim10^{11}$-$10^{12}$ M$_{\odot}$), associated with intense star formation activities (SFR $\sim$ 200-2000 M$_{\odot}$ yr$^{-1}$). The sample covers a range of AGN bolometric luminosity of $10^{10}$-$10^{13}$ L$_{\odot}$, and the outflow velocity measured from [OIII] 5007A line shows a correlation with AGN luminosity. Eight galaxies show extremely fast outflows with velocity up to 1500-2000 km s$^{-1}$. However, the co-existence of vigorous starbursts and strong outflows suggests the star formation has not been quenched during the ULIRG phase. By deriving stellar mass and mass fraction of young stellar population, we find no significant discrepancies between stellar properties of ULIRGs with weak and powerful AGNs. The results are not consistent with the merger-induced evolutionary scenario, which predicts that SF-dominated ULIRGs show smaller stellar mass and younger stellar populations compared to AGN-dominated ULIRGs. △ Less

Submitted 13 July, 2020; v1 submitted 11 November, 2019; originally announced November 2019.

Comments: 28 pages, 25 figures, 7 tables; accepted by ApJ

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). X. Discovery of 35 Quasars and Luminous Galaxies at 5.7 $\le$ z $\le$ 7.0

Authors: Yoshiki Matsuoka, Kazushi Iwasawa, Masafusa Onoue, Nobunari Kashikawa, Michael A. Strauss, Chien-Hsiu Lee, Masatoshi Imanishi, Tohru Nagao, Masayuki Akiyama, Naoko Asami, James Bosch, Hisanori Furusawa, Tomotsugu Goto, James E. Gunn, Yuichi Harikane, Hiroyuki Ikeda, Takuma Izumi, Toshihiro Kawaguchi, Nanako Kato, Satoshi Kikuta, Kotaro Kohno, Yutaka Komiyama, Shuhei Koyama, Robert H. Lupton, Takeo Minezaki , et al. (22 additional authors not shown)

Abstract: We report the discovery of 28 quasars and 7 luminous galaxies at 5.7 $\le$ z $\le$ 7.0. This is the tenth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The total number of spectroscopically identified objects in SHELLQs has… ▽ More We report the discovery of 28 quasars and 7 luminous galaxies at 5.7 $\le$ z $\le$ 7.0. This is the tenth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The total number of spectroscopically identified objects in SHELLQs has now grown to 93 high-z quasars, 31 high-z luminous galaxies, 16 [O III] emitters at z ~ 0.8, and 65 Galactic cool dwarfs (low-mass stars and brown dwarfs). These objects were found over 900 deg2, surveyed by HSC between 2014 March and 2018 January. The full quasar sample includes 18 objects with very strong and narrow Ly alpha emission, whose stacked spectrum is clearly different from that of other quasars or galaxies. While the stacked spectrum shows N V 1240 emission and resembles that of lower-z narrow-line quasars, the small Ly alpha width may suggest a significant contribution from the host galaxies. Thus these objects may be composites of quasars and star-forming galaxies. △ Less

Submitted 19 August, 2019; originally announced August 2019.

Comments: ApJ, in press. arXiv admin note: text overlap with arXiv:1803.01861

Journal ref: ApJ 883:183, 2019

Analysis and synthesis of feature map for kernel-based quantum classifier

Authors: Yudai Suzuki, Hiroshi Yano, Qi Gao, Shumpei Uno, Tomoki Tanaka, Manato Akiyama, Naoki Yamamoto

Abstract: A method for analyzing the feature map for the kernel-based quantum classifier is developed; that is, we give a general formula for computing a lower bound of the exact training accuracy, which helps us to see whether the selected feature map is suitable for linearly separating the dataset. We show a proof of concept demonstration of this method for a class of 2-qubit classifier, with several 2-di… ▽ More A method for analyzing the feature map for the kernel-based quantum classifier is developed; that is, we give a general formula for computing a lower bound of the exact training accuracy, which helps us to see whether the selected feature map is suitable for linearly separating the dataset. We show a proof of concept demonstration of this method for a class of 2-qubit classifier, with several 2-dimensional dataset. Also, a synthesis method, that combines different kernels to construct a better-performing feature map in a lager feature space, is presented. △ Less

Submitted 17 August, 2020; v1 submitted 25 June, 2019; originally announced June 2019.

Comments: 9 pages, 10 figures, 4 tables

Journal ref: Quantum Machine Intelligence 2, 1-9, 2020

An observational link between AGN Eddington ratio and [NII]λ6583/Ha at 0.6<z<1.7

Authors: Kyuseok Oh, Yoshihiro Ueda, Masayuki Akiyama, Hyewon Suh, Michael J. Koss, Daichi Kashino, Günther Hasinger

Abstract: We present an observed relationship between Eddington ratio (λEdd) and optical narrow-emission-line ratio ([N II]λ6583/Hα) of X-ray-selected broad-line active galactic nuclei (AGN) at 0.6 < z < 1.7. We use 27 near-infrared spectra from the Fiber Multi-Object Spectrograph along with 26 sources from the literature. We show that the λEdd and [N II]λ6583/Hα ratio at 0.6 < z < 1.7 exhibits a similar an… ▽ More We present an observed relationship between Eddington ratio (λEdd) and optical narrow-emission-line ratio ([N II]λ6583/Hα) of X-ray-selected broad-line active galactic nuclei (AGN) at 0.6 < z < 1.7. We use 27 near-infrared spectra from the Fiber Multi-Object Spectrograph along with 26 sources from the literature. We show that the λEdd and [N II]λ6583/Hα ratio at 0.6 < z < 1.7 exhibits a similar anti-correlation distribution of λEdd-[N II]λ6583/Hα as has been found for local (<z> = 0.036), ultra-hard X-ray selected AGN. The observed distribution implies that there is a consistent relationship from local to z \sim 1.7 which corresponds from the present time to 4 Gyr old. Further study of high redshift low Eddington ratio AGN (logλEdd < -2) is necessary to determine fully whether the λEdd-[N II]λ6583/Hα anti-correlation still holds in high-redshift AGN at low Eddington ratios. △ Less

Submitted 10 June, 2019; originally announced June 2019.

Comments: 15 pages, 7 figures, and 2 tables; Accepted for publication in ApJ

A Wide and Deep Exploration of Radio Galaxies with Subaru HSC (WERGS). II. Physical Properties derived from the SED Fitting with Optical, Infrared, and Radio Data

Authors: Yoshiki Toba, Takuji Yamashita, Tohru Nagao, Wei-Hao Wang, Yoshihiro Ueda, Kohei Ichikawa, Toshihiro Kawaguchi, Masayuki Akiyama, Bau-Ching Hsieh, Masaru Kajisawa, Chien-Hsiu Lee, Yoshiki Matsuoka, Akatoki Noboriguchi, Masafusa Onoue, Malte Schramm, Masayuki Tanaka, Yutaka Komiyama

Abstract: We present physical properties of radio galaxies (RGs) with $f_{\rm 1.4 GHz} >$ 1 mJy discovered by Subaru Hyper Supreme-Cam (HSC) and VLA Faint Images of the Radio Sky at Twenty-Centimeters (FIRST) survey. For 1056 FIRST RGs at $0 < z \leq 1.7$ with HSC counterparts in about 100 deg$^2$, we compiled multi-wavelength data of optical, near-infrared (IR), mid-IR, far-IR, and radio (150 MHz). We deri… ▽ More We present physical properties of radio galaxies (RGs) with $f_{\rm 1.4 GHz} >$ 1 mJy discovered by Subaru Hyper Supreme-Cam (HSC) and VLA Faint Images of the Radio Sky at Twenty-Centimeters (FIRST) survey. For 1056 FIRST RGs at $0 < z \leq 1.7$ with HSC counterparts in about 100 deg$^2$, we compiled multi-wavelength data of optical, near-infrared (IR), mid-IR, far-IR, and radio (150 MHz). We derived their color excess ($E (B-V)_{*}$), stellar mass, star formation rate (SFR), IR luminosity, the ratio of IR and radio luminosity ($q_{\rm IR}$), and radio spectral index ($α_{\rm radio}$) that are derived from the SED fitting with CIGALE. We also estimated Eddington ratio based on stellar mass and integration of the best-fit SEDs of AGN component. We found that $E (B-V)_{*}$, SFR, and IR luminosity clearly depend on redshift while stellar mass, $q_{\rm IR}$, and $α_{\rm radio}$ do not significantly depend on redshift. Since optically-faint ($i_{\rm AB} \geq 21.3$) RGs that are newly discovered by our RG survey tend to be high redshift, they tend to not only have a large dust extinction and low stellar mass but also have high SFR and AGN luminosity, high IR luminosity, and high Eddington ratio compared to optically-bright ones. The physical properties of a fraction of RGs in our sample seem to differ from a classical view of RGs with massive stellar mass, low SFR, and low Eddington ratio, demonstrating that our RG survey with HSC and FIRST provides us curious RGs among entire RG population. △ Less

Submitted 3 May, 2019; originally announced May 2019.

Comments: 30 pages, 20 figures, and 4 tables, accepted for publication in ApJS. The catalog and SED template of radio galaxies will be accessible through an online service

Spatially resolved stellar mass buildup and quenching in massive disk galaxies over the last 10 Gyr revealed with spatially resolved SED fitting

Authors: Abdurro'uf, Masayuki Akiyama

Abstract: Despite decreasing cosmic star formation rate density over the last 10 Gyr, the stellar mass ($M_{*}$) buildups in galaxies were still progressing during this epoch. About 50\% of the current $M_{*}$ density in the universe was built over the last $\sim 8.7$ Gyr. In this research, we investigated the stellar mass buildup and quenching of spatially resolved regions within massive disk galaxies over… ▽ More Despite decreasing cosmic star formation rate density over the last 10 Gyr, the stellar mass ($M_{*}$) buildups in galaxies were still progressing during this epoch. About 50\% of the current $M_{*}$ density in the universe was built over the last $\sim 8.7$ Gyr. In this research, we investigated the stellar mass buildup and quenching of spatially resolved regions within massive disk galaxies over the last 10 Gyr. We apply the spectral energy distribution (SED) fitting method to SEDs of sub-galactic regions in galaxies to derive the spatially resolved distributions of SFR and $M_{*}$ in the galaxies. This namely \textit{pixel-to-pixel SED fitting} method is applied to massive disk galaxies at $0.01<z<0.02$ and $0.8<z<1.8$. We found that massive disk galaxies tend to build their $M_{*}$ and quench their star formation progressively from the central region to the outskirts, i.e. \textit{inside-out} stellar mass buildup and quenching. △ Less

Submitted 20 February, 2019; originally announced February 2019.

Comments: 5 pages; to appear in the IAU proceeding PanModel2018: "Challenges in panchromatic galaxy modelling with next generation facilities"

Enhancement of HI absorption associated with the $z=3.1$ large-scale proto-cluster and characteristic structures with AGNs sculptured over Gpc scale in the SSA22 field

Authors: T. Hayashino, A. K. Inoue, K. Kousai, N. Kashikawa, K. Mawatari, Y. Matsuda, N. Tejos, J. X. Prochaska, I. Iwata, S. Noll, D. Burgarella, T. Yamada, M. Akiyama

Abstract: In the SSA22 field which exhibits a large-scale proto-cluster at $z=3.1$, we carried out a spectroscopic survey for Lyman Break Galaxies (LBGs) with the VLT/VIMOS and identified 78 confident LBGs at $z=2.5$--4. We stacked their spectra in the observer's frame by using a sophisticated method. Analyzing the composite spectrum, we have revealed that the large-scale proto-cluster at $z=3.1$ has a stro… ▽ More In the SSA22 field which exhibits a large-scale proto-cluster at $z=3.1$, we carried out a spectroscopic survey for Lyman Break Galaxies (LBGs) with the VLT/VIMOS and identified 78 confident LBGs at $z=2.5$--4. We stacked their spectra in the observer's frame by using a sophisticated method. Analyzing the composite spectrum, we have revealed that the large-scale proto-cluster at $z=3.1$ has a strong HI absorption dip of rest-frame equivalent width of $-1.7$ A. Another strong absorption dip found at $z=3.28$ is associated with a modestly high-density LBG peak, similar to that at $z=3.1$. We have also detected an HI transparency peak at $z=2.98$ in the composite spectrum, coincident with a void in the LBG distribution. In this paper, we also investigated the relation between LBGs, HI gas and AGNs at $z=3$--4 in the SSA22 field. Two AGNs at $z=3.353$ and 3.801 are, respectively, associated with the LBG concentration of an overdensity factor $δ_{\rm LBG}\simeq2$ in the present statistics. Another structure at $z=3.453$ is remarkable: 20 comoving Mpc scale dense HI gas which is not associated with any apparent LBG overdensity but involving a pair of AGNs. Such structure may be a new type of the AGN-matter correlation. If the inhomogeneous structures over a comoving Gpc scale found in this paper are confirmed with sufficient statistics in the future, the SSA22 field will become a key region to test the standard cold dark matter structure formation scenario. △ Less

Submitted 31 January, 2019; originally announced January 2019.

Comments: MNRAS accepted, 24 pages, 20 figures, 8 tables

Discovery of the First Low-Luminosity Quasar at z > 7

Authors: Yoshiki Matsuoka, Masafusa Onoue, Nobunari Kashikawa, Michael A. Strauss, Kazushi Iwasawa, Chien-Hsiu Lee, Masatoshi Imanishi, Tohru Nagao, Masayuki Akiyama, Naoko Asami, James Bosch, Hisanori Furusawa, Tomotsugu Goto, James E. Gunn, Yuichi Harikane, Hiroyuki Ikeda, Takuma Izumi, Toshihiro Kawaguchi, Nanako Kato, Satoshi Kikuta, Kotaro Kohno, Yutaka Komiyama, Shuhei Koyama, Robert H. Lupton, Takeo Minezaki , et al. (23 additional authors not shown)

Abstract: We report the discovery of a quasar at z = 7.07, which was selected from the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. This quasar, HSC J124353.93+010038.5, has an order of magnitude lower luminosity than do the other known quasars at z > 7. The rest-frame ultraviolet absolute magnitude is M1450 = -24.13 +/- 0.08 mag and the bolometric l… ▽ More We report the discovery of a quasar at z = 7.07, which was selected from the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. This quasar, HSC J124353.93+010038.5, has an order of magnitude lower luminosity than do the other known quasars at z > 7. The rest-frame ultraviolet absolute magnitude is M1450 = -24.13 +/- 0.08 mag and the bolometric luminosity is Lbol = (1.4 +/- 0.1) x 10^{46} erg/s. Its spectrum in the optical to near-infrared shows strong emission lines, and shows evidence for a fast gas outflow, as the C IV line is blueshifted and there is indication of broad absorption lines. The Mg II-based black hole mass is Mbh = (3.3 +/- 2.0) x 10^8 Msun, thus indicating a moderate mass accretion rate with an Eddington ratio 0.34 +/- 0.20. It is the first z > 7 quasar with sub-Eddington accretion, besides being the third most distant quasar, known to date. The luminosity and black hole mass are comparable to, or even lower than, those measured for the majority of low-z quasars discovered by the Sloan Digital Sky Survey, and thus this quasar likely represents a z > 7 counterpart to quasars commonly observed in the low-z universe. △ Less

Submitted 29 January, 2019; originally announced January 2019.

Comments: Accepted for publication in ApJ Letters

Journal ref: ApJL 872:L2, 2019

Discovery of a strong ionized-gas outflow in an AKARI-selected Ultra-luminous Infrared Galaxy at z = 0.5

Authors: Xiaoyang Chen, Masayuki Akiyama, Hirofumi Noda, Abdurrouf, Yoshiki Toba, Issei Yamamura, Toshihiro Kawaguchi, Mitsuru Kokubo, Kohei Ichikawa

Abstract: In order to construct a sample of ultra-luminous infrared galaxies (ULIRGs, with infrared luminosity, $L_{\rm IR} > 10^{12}$ L$_{\odot}$) at 0.5 < z < 1, we are conducting an optical follow-up program for bright 90-$μ$m FIR sources with a faint optical (i < 20) counterpart selected in the AKARI Far-Infrared Surveyor (FIS) Bright Source catalog (Ver.2). AKARI-FIS-V2 J0916248+073034, identified as a… ▽ More In order to construct a sample of ultra-luminous infrared galaxies (ULIRGs, with infrared luminosity, $L_{\rm IR} > 10^{12}$ L$_{\odot}$) at 0.5 < z < 1, we are conducting an optical follow-up program for bright 90-$μ$m FIR sources with a faint optical (i < 20) counterpart selected in the AKARI Far-Infrared Surveyor (FIS) Bright Source catalog (Ver.2). AKARI-FIS-V2 J0916248+073034, identified as a ULIRG at z = 0.49 in the spectroscopic follow-up observation, indicates signatures of an extremely strong outflow in its emission line profiles. Its [OIII] 5007Å emission line shows FWHM of 1830 km s$^{-1}$ and velocity shift of -770 km s$^{-1}$ in relative to the stellar absorption lines. Furthermore, low-ionization [OII] 3726Å 3729Å doublet also shows large FWHM of 910 km s$^{-1}$ and velocity shift of -380 km s$^{-1}$. After the removal of an unresolved nuclear component, the long-slit spectroscopy 2D image possibly shows that the outflow extends to radius of 4 kpc. The mass outflow and energy ejection rates are estimated to be 500 M$_{\odot}$ yr$^{-1}$ and $4\times10^{44}$ erg s$^{-1}$, respectively, which imply that the outflow is among the most powerful ones observed in ULIRGs and QSOs at 0.3 < z < 1.6. The co-existence of the strong outflow and intense star formation (star formation rate of 990 M$_{\odot}$ yr$^{-1}$) indicates that the feedback of the strong outflow has not severely affect the star-forming region of the galaxy. △ Less

Submitted 16 January, 2019; originally announced January 2019.

Comments: 20 pages, 18 figures, accepted to be published in PASJ