Showing 1–50 of 222 results for author: Silverman, J D

Rich and diverse molecular gas environments of closely-separated dual quasars viewed by ALMA

Authors: Shenli Tang, John D. Silverman, Zhaoxuan Liu, Manda Banerji, Tomoko Suzuki, Seiji Fujimoto, Andy Goulding, Masatoshi Imanishi, Toshihiro Kawaguchi, Connor Bottrell, Tilman Hartwig, Knud Jahnke, Masafusa Onoue, Malte Schramm, Yoshihiro Ueda

Abstract: We present a study of the molecular gas in five closely-spaced ($R_{\perp}<20$ kpc) dual quasars ($L_{\rm bol}\gtrsim10^{44}~\mathrm{erg~s}^{-1}$) at redshifts $0.4<z<0.8$ with the Atacama Large Millimeter/submillimeter Array. The dual quasar phase represents a distinctive stage during the interaction between two galaxies for investigating quasar fueling and feedback effects on the gas reservoir.… ▽ More We present a study of the molecular gas in five closely-spaced ($R_{\perp}<20$ kpc) dual quasars ($L_{\rm bol}\gtrsim10^{44}~\mathrm{erg~s}^{-1}$) at redshifts $0.4<z<0.8$ with the Atacama Large Millimeter/submillimeter Array. The dual quasar phase represents a distinctive stage during the interaction between two galaxies for investigating quasar fueling and feedback effects on the gas reservoir. The dual quasars were selected from the Sloan Digital Sky Survey and Subaru/Hyper Suprime-Cam Subaru Strategic Program, with confirmatory spectroscopic validation. Based on the detection of the CO J=2--1 emission line with Band 4, we derived key properties including CO luminosities, line widths, and molecular gas masses for these systems. Among the ten quasars of the five pairs, eight have line detections exceeding $5σ$. The detected sources prominently harbor substantial molecular gas reservoirs, with molecular gas masses ($M_{\text{molgas}}$) between $10^{9.6-10.5}~\mathrm{M_{\odot}}$, and molecular gas-to-stellar mass ratios ($μ_{\text{molgas}}$) spanning $18-97\%$. The overall $μ_{\text{molgas}}$ of these dual quasars agrees with that of inactive star-forming main-sequence galaxies at comparable redshifts, indicating no clear evidence of quenching. However, intriguing features in each individual system show possible evidence of AGN feedback, matter transfer, and compaction processes. △ Less

Submitted 12 July, 2024; originally announced July 2024.

Prior-Informed AGN-Host Spectral Decomposition Using PyQSOFit

Authors: Wenke Ren, Hengxiao Guo, Yue Shen, John D. Silverman, Colin J. Burke, Shu Wang, Junxian Wang

Abstract: We introduce an improved method for decomposing the emission of active galactic nuclei (AGN) and their host galaxies using templates from principal component analysis (PCA). This approach integrates prior information from PCA with a penalized pixel fitting mechanism which improves the precision and effectiveness of the decomposition process. Specifically, we have reduced the degeneracy and over-fi… ▽ More We introduce an improved method for decomposing the emission of active galactic nuclei (AGN) and their host galaxies using templates from principal component analysis (PCA). This approach integrates prior information from PCA with a penalized pixel fitting mechanism which improves the precision and effectiveness of the decomposition process. Specifically, we have reduced the degeneracy and over-fitting in AGN-host decomposition, particularly for those with low signal-to-noise ratios (SNR), where traditional methods tend to fail. By applying our method to 76,565 SDSS Data Release 16 quasars with $z<0.8$, we achieve a success rate of $\approx$ 94%, thus establishing the largest host-decomposed spectral catalog of quasars to date. Our fitting results consider the impact of the host galaxy on the overestimation of the AGN luminosity and black hole mass ($M_{\rm BH}$). Furthermore, we obtained stellar velocity dispersion ($σ_*$) measurements for 4,137 quasars. The slope of the $M_{\rm BH}-σ_*$ relation in this subsample is generally consistent with previous quasar studies beyond the local universe. Our method provides a robust and efficient approach to disentangle the AGN and host galaxy components across a wide range of SNRs and redshifts. △ Less

Submitted 25 June, 2024; originally announced June 2024.

Comments: 21 pages, 9 figures, 2 tables. Submitted to ApJ, comments welcome

COSMOS-Web: The over-abundance and physical nature of "little red dots"--Implications for early galaxy and SMBH assembly

Authors: Hollis B. Akins, Caitlin M. Casey, Erini Lambrides, Natalie Allen, Irham T. Andika, Malte Brinch, Jaclyn B. Champagne, Olivia Cooper, Xuheng Ding, Nicole E. Drakos, Andreas Faisst, Steven L. Finkelstein, Maximilien Franco, Seiji Fujimoto, Fabrizio Gentile, Steven Gillman, Ghassem Gozaliasl, Santosh Harish, Christopher C. Hayward, Michaela Hirschmann, Olivier Ilbert, Jeyhan S. Kartaltepe, Dale D. Kocevski, Anton M. Koekemoer, Vasily Kokorev , et al. (16 additional authors not shown)

Abstract: JWST has revealed a population of compact and extremely red galaxies at $z>4$, which likely host active galactic nuclei (AGN). We present a sample of 434 ``little red dots'' (LRDs), selected from the 0.54 deg$^2$ COSMOS-Web survey. We fit galaxy and AGN SED models to derive redshifts and physical properties; the sample spans $z\sim5$-$9$ after removing brown dwarf contaminants. We consider two ext… ▽ More JWST has revealed a population of compact and extremely red galaxies at $z>4$, which likely host active galactic nuclei (AGN). We present a sample of 434 ``little red dots'' (LRDs), selected from the 0.54 deg$^2$ COSMOS-Web survey. We fit galaxy and AGN SED models to derive redshifts and physical properties; the sample spans $z\sim5$-$9$ after removing brown dwarf contaminants. We consider two extreme physical scenarios: either LRDs are all AGN, and their continuum emission is dominated by the accretion disk, or they are all compact star-forming galaxies, and their continuum is dominated by stars. If LRDs are AGN-dominated, our sample exhibits bolometric luminosities $\sim10^{45-47}$ erg\,s$^{-1}$, spanning the gap between JWST AGN in the literature and bright, rare quasars. We derive a bolometric luminosity function (LF) $\sim100$ times the (UV-selected) quasar LF, implying a non-evolving black hole accretion density of $\sim10^{-4}$ M$_\odot$ yr$^{-1}$ Mpc$^{-3}$ from $z\sim2$-$9$. By contrast, if LRDs are dominated by star formation, we derive stellar masses $\sim10^{8.5-10}\,M_\odot$. MIRI/F770W is key to deriving accurate stellar masses; without it, we derive a mass function inconsistent with $Λ$CDM. The median stellar mass profile is broadly consistent with the maximal stellar mass surface densities seen in the nearby universe, though the most massive $\sim50$\% of objects exceed this limit, requiring substantial AGN contribution to the continuum. Nevertheless, stacking all available X-ray, mid-IR, far-IR/sub-mm, and radio data yields non-detections. Whether dominated by dusty AGN, compact star-formation, or both, the high masses/luminosities and remarkable abundance of LRDs implies a dominant mode of early galaxy/SMBH growth. △ Less

Submitted 14 June, 2024; originally announced June 2024.

Comments: 27 pages, 13 figures. Submitted to ApJ, comments welcome! Data access at https://github.com/hollisakins/akins24_cw

Systematic Collapse of the Accretion Disc Across the Supermassive Black Hole Population

Authors: Scott Hagen, Chris Done, John D. Silverman, Junyao Li, Teng Liu, Wenke Ren, Johannes Buchner, Andrea Merloni, Tohru Nagao, Mara Salvato

Abstract: The structure of the accretion flow onto supermassive black holes (SMBH) is not well understood. Standard disc models match to zeroth order in predicting substantial energy dissipation within optically-thick material producing a characteristic strong blue/UV continuum. However they fail at reproducing more detailed comparisons to the observed spectral shapes along with their observed variability.… ▽ More The structure of the accretion flow onto supermassive black holes (SMBH) is not well understood. Standard disc models match to zeroth order in predicting substantial energy dissipation within optically-thick material producing a characteristic strong blue/UV continuum. However they fail at reproducing more detailed comparisons to the observed spectral shapes along with their observed variability. Based on stellar mass black holes within our galaxy, accretion discs should undergo a transition into an X-ray hot, radiatively inefficient flow, below a (mass scaled) luminosity of $\sim 0.02\,L_{\rm{Edd}}$. While this has been seen in limited samples of nearby low-luminosity active galactic nuclei (AGN) and a few rare changing-look AGN, it is not at all clear whether this transition is present in the wider AGN population across cosmic time. A key issue is the difficulty in disentangling a change in spectral state from increased dust obscuration and/or host galaxy contamination, effectively drowning out the AGN emission. Here we use the new eROSITA eFEDS Survey to identify unobscured AGN from their X-ray emission, matched to excellent optical imaging from Subaru's Hyper Suprime-Cam; allowing the subtraction of the host galaxy contamination. The resulting, uncontaminated, AGN spectra reveal a smooth transition from a strongly disc dominated state in bright AGN, to the collapse of the disc into an inefficient X-ray plasma in the low luminosity AGN, with the transition occurring at $\sim 0.02\,L_{\rm{Edd}}$; revealing fundamental aspects of accretion physics in AGN. △ Less

Submitted 10 June, 2024; originally announced June 2024.

Comments: 10 pages, 5 figure, 1 appendix - submitted to MNRAS, comments welcome

COSMOS-Web: The Role of Galaxy Interactions and Disk Instabilities in Producing Starbursts at z<4

Authors: A. L. Faisst, M. Brinch, C. M. Casey, N. Chartab, M. Dessauges-Zavadsky, N. E. Drakos, S. Gillman, G. Gonzaliasl, C. C. Hayward, O. Ilbert, P. Jablonka, J. S. Kartaltepe, A. M. Koekemoer, V. Kokorev, E. Lambrides, D. Liu, C. Maraston, C. L. Martin, A. Renzini, B. E. Robertson, D. B. Sanders, Z. Sattari, N. Scoville, C. M. Urry, A. P. Vijayan , et al. (27 additional authors not shown)

Abstract: We study of the role of galaxy-galaxy interactions and disk instabilities in producing starburst activity in galaxies out to z=4. For this, we use a sample of 387 galaxies with robust total star formation rate measurements from Herschel, gas masses from ALMA, stellar masses and redshifts from multi-band photometry, and JWST/NIRCam rest-frame optical imaging. Using mass-controlled samples, we find… ▽ More We study of the role of galaxy-galaxy interactions and disk instabilities in producing starburst activity in galaxies out to z=4. For this, we use a sample of 387 galaxies with robust total star formation rate measurements from Herschel, gas masses from ALMA, stellar masses and redshifts from multi-band photometry, and JWST/NIRCam rest-frame optical imaging. Using mass-controlled samples, we find an increased fraction of interacting galaxies in the starburst regime at all redshifts out to z=4. This increase correlates with star formation efficiency (SFE), but not with gas fraction. However, the correlation is weak (and only significant out to z=2), which could be explained by the short duration of SFE increase during interaction. In addition, we find that isolated disk galaxies make up a significant fraction of the starburst population. The fraction of such galaxies with star-forming clumps ("clumpy disks") is significantly increased compared to the main-sequence disk population. Furthermore, this fraction directly correlates with SFE. This is direct observational evidence for a long-term increase of SFE maintained due to disk instabilities, contributing to the majority of starburst galaxies in our sample and hence to substantial mass growth in these systems. This result could also be of importance for explaining the growth of the most massive galaxies at z>6. △ Less

Submitted 15 May, 2024; originally announced May 2024.

Comments: 14 pages, 5 figures. Submitted to ApJ. Comments welcome

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

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

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

Submitted 3 May, 2024; originally announced May 2024.

Comments: 14 pages incl. Appendix, 4 main figures, 1 table, 1 appendix figure. Submitted to AAS journal

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)

Cosmic Himalayas: The Highest Quasar Density Peak Identified in a 10,000 deg$^2$ Sky with Spatial Discrepancies between Galaxies, Quasars, and IGM HI

Authors: Yongming Liang, Masami Ouchi, Dongsheng Sun, Nobunari Kashikawa, Zheng Cai, Sebastiano Cantalupo, Kentaro Nagamine, Hidenobu Yajima, Takanobu Kirihara, Haibin Zhang, Mingyu Li, Rhythm Shimakawa, Xiaohui Fan, Kei Ito, Masayuki Tanaka, Yuichi Harikane, J. Xavier Prochaska, Andrea Travascio, Weichen Wang, Martin Elvis, Giuseppina Fabbiano, Junya Arita, Masafusa Onoue, John D. Silverman, Dongdong Shi , et al. (5 additional authors not shown)

Abstract: We report the identification of a quasar overdensity in the BOSSJ0210 field, dubbed Cosmic Himalayas, consisting of 11 quasars at $z=2.16-2.20$, the densest overdensity of quasars ($17σ$) in the $\sim$10,000 deg$^2$ of the Sloan Digital Sky Survey. We present the spatial distributions of galaxies and quasars and an HI absorption map of the intergalactic medium (IGM). On the map of 465 galaxies sel… ▽ More We report the identification of a quasar overdensity in the BOSSJ0210 field, dubbed Cosmic Himalayas, consisting of 11 quasars at $z=2.16-2.20$, the densest overdensity of quasars ($17σ$) in the $\sim$10,000 deg$^2$ of the Sloan Digital Sky Survey. We present the spatial distributions of galaxies and quasars and an HI absorption map of the intergalactic medium (IGM). On the map of 465 galaxies selected from the MAMMOTH-Subaru survey, we find two galaxy density peaks that do not fall on the quasar overdensity but instead exist at the northwest and southeast sides, approximately 25 $h^{-1}$ comoving-Mpc apart from the quasar overdensity. With a spatial resolution of 15 $h^{-1}$ comoving Mpc in projection, we produce a three-dimensional HI tomography map by the IGM Ly$α$ forest in the spectra of 23 SDSS/eBOSS quasars behind the quasar overdensity. Surprisingly, the quasar overdensity coincides with neither an absorption peak nor a transmission peak of IGM HI but lies near the border separating opaque and transparent volumes, with the more luminous quasars located in an environment with lesser IGM HI. Hence remarkably, the overdensity region traced by the 11 quasars, albeit all in coherently active states, has no clear coincidence with peaks of galaxies or HI absorption densities. Current physical scenarios with mixtures of HI overdensities and quasar photoionization cannot fully interpret the emergence of Cosmic Himalayas, suggesting this peculiar structure is an excellent laboratory to unveil the interplay between galaxies, quasars, and the IGM. △ Less

Submitted 24 April, 2024; originally announced April 2024.

Comments: 19 pages, 11 figures, submitted to ApJ, comments are welcome

Tip of the iceberg: overmassive black holes at 4<z<7 found by JWST are not inconsistent with the local $\mathcal{M}_{\rm BH}$-$\mathcal{M}_\star$ relation

Authors: Junyao Li, John D. Silverman, Yue Shen, Marta Volonteri, Knud Jahnke, Ming-Yang Zhuang, Matthew T. Scoggins, Xuheng Ding, Yuichi Harikane, Masafusa Onoue, Takumi S. Tanaka

Abstract: JWST is revealing a new remarkable population of high-redshift ($z\gtrsim4$), low-luminosity Active Galactic Nuclei (AGNs) in deep surveys and detecting the host galaxy stellar light in the most luminous and massive quasars at $z\sim 6$ for the first time. Latest results claim supermassive black holes (SMBHs) in these systems to be significantly more massive than expected from the local BH mass -… ▽ More JWST is revealing a new remarkable population of high-redshift ($z\gtrsim4$), low-luminosity Active Galactic Nuclei (AGNs) in deep surveys and detecting the host galaxy stellar light in the most luminous and massive quasars at $z\sim 6$ for the first time. Latest results claim supermassive black holes (SMBHs) in these systems to be significantly more massive than expected from the local BH mass - stellar mass ($\mathcal{M}_{\rm BH} - \mathcal{M}_\star$) relation and that this is not due to sample selection effects. Through detailed statistical modeling, we demonstrate that the coupled effects of selection biases (i.e., finite detection limit and requirements on detecting broad lines) and measurement uncertainties in $\mathcal{M}_{\rm BH}$ and $\mathcal{M}_\star$ can in fact largely account for the reported offset and flattening in the observed $\mathcal{M}_{\rm BH} - \mathcal{M}_\star$ relation toward the upper envelope of the local relation, even for those at $\mathcal{M}_{\rm BH} < 10^8\,M_{\odot}$. We further investigate the possible evolution of the $\mathcal{M}_{\rm BH} - \mathcal{M}_\star$ relation at $z\gtrsim 4$ with careful treatment of observational biases and consideration of the degeneracy between intrinsic evolution and dispersion in this relation. The bias-corrected intrinsic $\mathcal{M}_{\rm BH} - \mathcal{M}_\star$ relation in the low-mass regime suggests that there might be a large population of low-mass BHs (${\rm log}\,\mathcal{M}_{\rm BH} \lesssim 5$), possibly originating from lighter seeds, remaining undetected or unidentified even in the deepest JWST surveys. These results have important consequences for JWST studies of BH seeding and the coevolution between SMBHs and their host galaxies at the earliest cosmic times. △ Less

Submitted 29 February, 2024; originally announced March 2024.

Comments: 21 pages, 9 figures, submitted to ApJ

Steep-spectrum AGN in eROSITA Final Equatorial-Depth Survey (eFEDS): Their host galaxies and multi-wavelength properties

Authors: K. Iwasawa, T. Liu, Th. Boller, J. Buchner, J. Li, T. Kawaguchi, T. Nagao, Y. Terashima, Y. Toba, J. D. Silverman, R. Arcodia, Th. Dauser, M. Krumpe, K. Nandra, J. Wilms

Abstract: We selected sources with a steep soft-X-ray-band spectrum with a photon index larger than 2.5 -- measured by eROSITA on board the Spectrum-Roentgen-Gamma (SRG) -- from the eFEDS AGN catalogue as candidates of highly accreting supermassive black holes, and investigated their multi-wavelength properties. Among 601 bright AGN with 0.2-5 keV counts of greater than 100, 83 sources (~14%) are classified… ▽ More We selected sources with a steep soft-X-ray-band spectrum with a photon index larger than 2.5 -- measured by eROSITA on board the Spectrum-Roentgen-Gamma (SRG) -- from the eFEDS AGN catalogue as candidates of highly accreting supermassive black holes, and investigated their multi-wavelength properties. Among 601 bright AGN with 0.2-5 keV counts of greater than 100, 83 sources (~14%) are classified as steep-spectrum sources. These sources have typical 0.5-2 keV luminosities of L(SX) ~ 1e44 erg/s and the majority of them are found at redshifts below z=1. In comparison with sources with flatter spectra, these sources have, on average, a UV (or optical) to 2 keV luminosity ratio that is larger by ~0.3 dex and bluer optical-to-UV continuum emission. They also appear to be radio quiet based on the detection rate in the FIRST and VLASS surveys. Their host galaxies -- at least in the redshift range of z=0.2-0.8, where the AGN-galaxy decomposition results from the Subaru Hyper Suprime-Cam imaging are available -- tend to be late-type and have smaller stellar masses than those of sources with flatter spectra. These properties are similar to those found in nearby narrow-line Seyfert 1 galaxies, in agreement with the picture that they are AGN with elevated accretion rates and are in the early growth phase of black hole and galaxy co-evolution. However, the steep-spectrum sources are not exclusively narrow-line Seyfert 1 galaxies; indeed many are broad-line Seyfert 1 galaxies, as found by a catalogue search. This suggests that these steep-spectrum sources may be black holes generally with high accretion rates but of a wide mass range, including a few objects emitting at L(SX)>1e45 erg/s, of which black hole masses can be close to 10^9 M_sun. △ Less

Submitted 21 February, 2024; originally announced February 2024.

Comments: 8 pages, accepted for publication in A&A

Near-IR clumps and their properties in high-z galaxies with JWST/NIRCam

Authors: Boris S. Kalita, John D. Silverman, Emanuele Daddi, Wilfried Mercier, Luis C. Ho, Xuheng Ding

Abstract: Resolved stellar morphology of z>1 galaxies was inaccessible before JWST. This limitation, due to the impact of dust on rest-frame UV light, had withheld major observational conclusions required to understand the importance of clumps in galaxy evolution. Essentially independent of this issue, we use the rest-frame near-IR for a stellar-mass dependent clump detection method and determine reliable e… ▽ More Resolved stellar morphology of z>1 galaxies was inaccessible before JWST. This limitation, due to the impact of dust on rest-frame UV light, had withheld major observational conclusions required to understand the importance of clumps in galaxy evolution. Essentially independent of this issue, we use the rest-frame near-IR for a stellar-mass dependent clump detection method and determine reliable estimations of selection effects. We exploit publicly available JWST/NIRCam and HST/ACS imaging data from CEERS, to create a stellar-mass based picture of clumps in a mass-complete sample of 418 galaxies within a wide wavelength coverage of 0.5-4.6$μ$m and a redshift window of 1<z<2. We find that a near-IR detection gives access to a larger set of clumps within galaxies, with those also detected in UV making up only 28%. Whereas, 85% of the UV clumps are found to have a near-IR counterpart. These near-IR clumps closely follow the UVJ classification of their respective host galaxies, with these hosts mainly populating the star-forming regime besides a fraction of them (16%) that can be considered quiescent. The mass of the detected clumps are found to be within the range of $10^{7.5-9.5}\,\rm M_{\odot}$, therefore expected to drive gas into galaxy cores through tidal torques. However, there is likely contribution from blending of smaller unresolved structures. Furthermore, we observe a radial gradient of increasing clump mass towards the centre of galaxies. This trend could be an indication of clump migration, but accurate star-formation measurements would be required to confirm such a scenario. △ Less

Submitted 4 February, 2024; originally announced February 2024.

Comments: Submitted to MNRAS. Comments welcome

The $M_{\rm BH}-M_*$ relation up to $z\sim2$ through decomposition of COSMOS-Web NIRCam images

Authors: Takumi S. Tanaka, John D. Silverman, Xuheng Ding, Knud Jahnke, Benny Trakhtenbrot, Erini Lambrides, Masafusa Onoue, Irham Taufik Andika, Angela Bongiorno, Andreas L. Faisst, Steven Gillman, Christopher C. Hayward, Michaela Hirschmann, Anton Koekemoer, Vasily Kokorev, Zhaoxuan Liu, Georgios E. Magdis, Alvio Renzini, Caitlin Casey, Nicole E. Drakos, Maximilien Franco, Ghassem Gozaliasl, Jeyhan Kartaltepe, Daizhong Liu, Henry Joy McCracken , et al. (3 additional authors not shown)

Abstract: Our knowledge of relations between supermassive black holes and their host galaxies at $z\gtrsim1$ is still limited, even though being actively sought out to $z\sim6$. Here, we use the high resolution and sensitivity of JWST to measure the host galaxy properties for 61 X-ray-selected type-I AGNs at $0.7<z<2.5$ with rest-frame optical/near-infrared imaging from COSMOS-Web and PRIMER. Black hole mas… ▽ More Our knowledge of relations between supermassive black holes and their host galaxies at $z\gtrsim1$ is still limited, even though being actively sought out to $z\sim6$. Here, we use the high resolution and sensitivity of JWST to measure the host galaxy properties for 61 X-ray-selected type-I AGNs at $0.7<z<2.5$ with rest-frame optical/near-infrared imaging from COSMOS-Web and PRIMER. Black hole masses ($\log\left(M_{\rm BH}/M_\odot\right)\sim7.5-9.5$) are available from previous spectroscopic campaigns. We extract the host galaxy components from four NIRCam broadband images and the HST/ACS F814W image by applying a 2D image decomposition technique. We detect the host galaxy for $\sim90\%$ of the sample after subtracting the unresolved AGN emission. With host photometry free of AGN emission, we determine the stellar mass of the host galaxies to be $\log\left(M_*/M_\odot\right)\sim10-11.5$ through SED fitting and measure the evolution of the mass relation between SMBHs and their host galaxies. Considering selection biases and measurement uncertainties, we find that the $M_\mathrm{ BH}/M_*$ ratio evolves as $\left(1+z\right)^{0.37_{-0.60}^{+0.35}}$ thus remains essentially constant or exhibits mild evolution up to $z\sim2.5$. We also see an amount of scatter ($σ_μ=0.28\pm0.13$) is similar to the local relation and consistent with low-$z$ studies; this appears to not rule out non-causal cosmic assembly where mergers contribute to the statistical averaging towards the local relation. We highlight improvements to come with larger samples from JWST and, particularly, Euclid, which will exceed the statistical power of wide and deep surveys such as Subaru Hyper Suprime-Cam. △ Less

Submitted 24 January, 2024; originally announced January 2024.

Comments: 31 pages, 19 figures, submitted to ApJ

Tracing the rise of supermassive black holes: A panchromatic search for faint, unobscured quasars at z > 6 with COSMOS-Web and other surveys

Authors: Irham T. Andika, Knud Jahnke, Masafusa Onoue, John D. Silverman, Itsna K. Fitriana, Angela Bongiorno, Malte Brinch, Caitlin M. Casey, Andreas Faisst, Steven Gillman, Ghassem Gozaliasl, Christopher C. Hayward, Michaela Hirschmann, Dale Kocevski, Anton M. Koekemoer, Vasily Kokorev, Erini Lambrides, Minju M. Lee, R. Michael Rich, Benny Trakhtenbrot, C. Megan Urry, Stephen M. Wilkins, Aswin P. Vijayan

Abstract: We report the identification of 64 new candidates of compact galaxies, potentially hosting faint quasars with bolometric luminosities of $L_\mathrm{bol} = 10^{43}$--10$^{46}$ erg s$^{-1}$, residing in the reionization epoch within the redshift range of $6 \lesssim z \lesssim 8$. These candidates were selected by harnessing the rich multiband datasets provided by the emerging JWST-driven extragalac… ▽ More We report the identification of 64 new candidates of compact galaxies, potentially hosting faint quasars with bolometric luminosities of $L_\mathrm{bol} = 10^{43}$--10$^{46}$ erg s$^{-1}$, residing in the reionization epoch within the redshift range of $6 \lesssim z \lesssim 8$. These candidates were selected by harnessing the rich multiband datasets provided by the emerging JWST-driven extragalactic surveys, focusing on COSMOS-Web, as well as JADES, UNCOVER, CEERS, and PRIMER. Our search strategy includes two stages: applying stringent photometric cuts to catalog-level data and detailed spectral energy distribution fitting. These techniques effectively isolate the quasar candidates while mitigating contamination from low-redshift interlopers, such as brown dwarfs and nearby galaxies. The selected candidates indicate physical traits compatible with low-luminosity active galactic nuclei, likely hosting $\approx10^5$--$10^7~M_\odot$ supermassive black holes (SMBHs) living in galaxies with stellar masses of $\approx10^8$--$10^{10}~M_\odot$. The SMBHs selected in this study, on average, exhibit an elevated mass compared to their hosts, with the mass ratio distribution slightly higher than those of galaxies in the local Universe. As with other high-$z$ studies, this is at least in part due to the selection method for these quasars. An extensive Monte Carlo analysis provides compelling evidence that heavy black hole seeds from the direct collapse scenario appear to be the preferred pathway to mature this specific subset of SMBHs by $z\approx7$. This work underscores the significance of further spectroscopic observations, as the quasar candidates presented here offer exceptional opportunities to delve into the nature of the earliest galaxies and SMBHs that formed during cosmic infancy. △ Less

Submitted 2 February, 2024; v1 submitted 22 January, 2024; originally announced January 2024.

Comments: Accepted for publication in the Astronomy & Astrophysics journal. 19 pages, 10 figures, and 4 tables. We welcome comments from the reader

Fitting pseudo-S${\rm \acute{e}}$rsic(Spergel) light profiles to galaxies in interferometric data: the excellence of the $uv$-plane

Authors: Qing-Hua Tan, Emanuele Daddi, Victor de Souza Magalhães, Carlos Gómez-Guijarro, Jérôme Pety, Boris S. Kalita, David Elbaz, Zhaoxuan Liu, Benjamin Magnelli, Annagrazia Puglisi, Wiphu Rujopakarn, John D. Silverman, Francesco Valentino, Shao-Bo Zhang

Abstract: Modern (sub)millimeter interferometers, such as ALMA and NOEMA, offer high angular resolution and unprecedented sensitivity. This provides the possibility to characterize the morphology of the gas and dust in distant galaxies. To assess the capabilities of current softwares in recovering morphologies and surface brightness profiles in interferometric observations, we test the performance of the Sp… ▽ More Modern (sub)millimeter interferometers, such as ALMA and NOEMA, offer high angular resolution and unprecedented sensitivity. This provides the possibility to characterize the morphology of the gas and dust in distant galaxies. To assess the capabilities of current softwares in recovering morphologies and surface brightness profiles in interferometric observations, we test the performance of the Spergel model for fitting in the $uv$-plane, which has been recently implemented in the IRAM software GILDAS (uv$\_$fit). Spergel profiles provide an alternative to the Sersic profile, with the advantage of having an analytical Fourier transform, making them ideal to model visibilities in the $uv$-plane. We provide an approximate conversion between Spergel index and Sersic index, which depends on the ratio of the galaxy size to the angular resolution of the data. We show through extensive simulations that Spergel modeling in the $uv$-plane is a more reliable method for parameter estimation than modeling in the image-plane, as it returns parameters that are less affected by systematic biases and results in a higher effective signal-to-noise ratio (S/N). The better performance in the $uv$-plane is likely driven by the difficulty of accounting for correlated signal in interferometric images. Even in the $uv$-plane, the integrated source flux needs to be at least 50 times larger than the noise per beam to enable a reasonably good measurement of a Spergel index. We characterise the performance of Spergel model fitting in detail by showing that parameters biases are generally low (< 10%) and that uncertainties returned by uv$\_$fit are reliable within a factor of two. Finally, we showcase the power of Spergel fitting by re-examining two claims of extended halos around galaxies from the literature, showing that galaxies and halos can be successfully fitted simultaneously with a single Spergel model. △ Less

Submitted 8 December, 2023; originally announced December 2023.

Comments: 23 pages, 15 figures, accepted for publication in A&A

JWST and ALMA discern the assembly of structural and obscured components in a high-redshift starburst galaxy

Authors: Zhaoxuan Liu, John D. Silverman, Emanuele Daddi, Annagrazia Puglisi, Alvio Renzini, Boris S. Kalita, Jeyhan S. Kartaltepe, Daichi Kashino, Giulia Rodighiero, Wiphu Rujopakarn, Tomoko L. Suzuki, Takumi S. Tanaka, Francesco Valentino, Irham Taufik Andika, Caitlin M. Casey, Andreas Faisst, Maximilien Franco, Ghassem Gozaliasl, Steven Gillman, Christopher C. Hayward, Anton M. Koekemoer, Vasily Kokorev, Erini Lambrides, Minju M. Lee, Georgios E. Magdis , et al. (5 additional authors not shown)

Abstract: We present observations and analysis of the starburst, PACS-819, at z=1.45 ($M_*=10^{10.7}$ M$_{ \odot}$), using high-resolution ($0^{\prime \prime}.1$; 0.8 kpc) ALMA and multi-wavelength JWST images from the COSMOS-Web program. Dissimilar to HST/ACS images in the rest-frame UV, the redder NIRCam and MIRI images reveal a smooth central mass concentration and spiral-like features, atypical for such… ▽ More We present observations and analysis of the starburst, PACS-819, at z=1.45 ($M_*=10^{10.7}$ M$_{ \odot}$), using high-resolution ($0^{\prime \prime}.1$; 0.8 kpc) ALMA and multi-wavelength JWST images from the COSMOS-Web program. Dissimilar to HST/ACS images in the rest-frame UV, the redder NIRCam and MIRI images reveal a smooth central mass concentration and spiral-like features, atypical for such an intense starburst. Through dynamical modeling of the CO J=5--4 emission with ALMA, PACS-819 is rotation-dominated thus has a disk-like nature. However, kinematic anomalies in CO and asymmetric features in the bluer JWST bands (e.g., F150W) support a more disturbed nature likely due to interactions. The JWST imaging further enables us to map the distribution of stellar mass and dust attenuation, thus clarifying the relationships between different structural components, not discernable in the previous HST images. The CO J = 5 -- 4 and FIR dust continuum emission are co-spatial with a heavily-obscured starbursting core (<1 kpc) which is partially surrounded by much less obscured star-forming structures including a prominent arc, possibly a tidally-distorted dwarf galaxy, and a clump, either a sign of an ongoing violent disk instability or a recently accreted low-mass satellite. With spatially-resolved maps, we find a high molecular gas fraction in the central area reaching $\sim3$ ($M_{\text{gas}}$/$M_*$) and short depletion times ($M_{\text{gas}}/SFR\sim$ 120 Myrs) across the entire system. These observations provide insights into the complex nature of starbursts in the distant universe and underscore the wealth of complementary information from high-resolution observations with both ALMA and JWST. △ Less

Submitted 10 May, 2024; v1 submitted 24 November, 2023; originally announced November 2023.

Comments: 21 pages(including appendix), 13 figures, accepted by ApJ

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

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

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

Submitted 27 September, 2023; originally announced September 2023.

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

A rest-frame near-IR study of clumps in galaxies at 1 < z < 2 using JWST/NIRCam: connection to galaxy bulges

Authors: Boris S. Kalita, John D. Silverman, Emanuele Daddi, Connor Bottrell, Luis C. Ho, Xuheng Ding, Lilan Yang

Abstract: A key question in galaxy evolution has been the importance of the apparent `clumpiness' of high redshift galaxies. Until now, this property has been primarily investigated in rest-frame UV, limiting our understanding of their relevance. Are they short-lived or are associated with more long-lived massive structures that are part of the underlying stellar disks? We use JWST/NIRCam imaging from CEERS… ▽ More A key question in galaxy evolution has been the importance of the apparent `clumpiness' of high redshift galaxies. Until now, this property has been primarily investigated in rest-frame UV, limiting our understanding of their relevance. Are they short-lived or are associated with more long-lived massive structures that are part of the underlying stellar disks? We use JWST/NIRCam imaging from CEERS to explore the connection between the presence of these `clumps' in a galaxy and its overall stellar morphology, in a mass-complete ($log\,M_{*}/M_{\odot} > 10.0$) sample of galaxies at $1.0 < z < 2.0$. Exploiting the uninterrupted access to rest-frame optical and near-IR light, we simultaneously map the clumps in galactic disks across our wavelength coverage, along with measuring the distribution of stars among their bulges and disks. Firstly, we find that the clumps are not limited to rest-frame UV and optical, but are also apparent in near-IR with $\sim 60\,\%$ spatial overlap. This rest-frame near-IR detection indicates that clumps would also feature in the stellar-mass distribution of the galaxy. A secondary consequence is that these will hence be expected to increase the dynamical friction within galactic disks leading to gas inflow. We find a strong negative correlation between how clumpy a galaxy is and strength of the bulge. This firmly suggests an evolutionary connection, either through clumps driving bulge growth, or the bulge stabilizing the galaxy against clump formation, or a combination of the two. Finally, we find evidence of this correlation differing from rest-frame optical to near-IR, which could suggest a combination of varying formation modes for the clumps. △ Less

Submitted 29 November, 2023; v1 submitted 11 September, 2023; originally announced September 2023.

Comments: Accepted for publication in ApJ

IllustrisTNG in the HSC-SSP: image data release and the major role of mini mergers as drivers of asymmetry and star formation

Authors: Connor Bottrell, Hassen M. Yesuf, Gergö Popping, Kiyoaki Christopher Omori, Shenli Tang, Xuheng Ding, Annalisa Pillepich, Dylan Nelson, Lukas Eisert, Hua Gao, Andy D. Goulding, Boris S. Kalita, Wentao Luo, Jenny E. Greene, Jingjing Shi, John D. Silverman

Abstract: At fixed galaxy stellar mass, there is a clear observational connection between structural asymmetry and offset from the star forming main sequence, $Δ$SFMS. Herein, we use the TNG50 simulation to investigate the relative roles of major mergers (stellar mass ratios $μ\geq0.25$), minor ($0.1 \leq μ< 0.25$), and mini mergers ($0.01 \leq μ< 0.1$) in driving this connection amongst star forming galaxi… ▽ More At fixed galaxy stellar mass, there is a clear observational connection between structural asymmetry and offset from the star forming main sequence, $Δ$SFMS. Herein, we use the TNG50 simulation to investigate the relative roles of major mergers (stellar mass ratios $μ\geq0.25$), minor ($0.1 \leq μ< 0.25$), and mini mergers ($0.01 \leq μ< 0.1$) in driving this connection amongst star forming galaxies (SFGs). We use dust radiative transfer post-processing with SKIRT to make a large, public collection of synthetic Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) images of simulated TNG galaxies over $0.1\leq z \leq 0.7$ with $\log M_{\star} / \mathrm{M}_{\odot}\geq9$ ($\sim750$k images). Using their instantaneous SFRs, known merger histories/forecasts, and HSC-SSP asymmetries, we show (1) that TNG50 SFGs qualitatively reproduce the observed trend between $Δ$SFMS and asymmetry and (2) a strikingly similar trend emerges between $Δ$SFMS and the time-to-coalescence for mini mergers. Controlling for redshift, stellar mass, environment, and gas fraction, we show that individual mini merger events yield small enhancements in SFRs and asymmetries that are sustained on long timescales (at least $\sim3$ Gyr after coalescence, on average) -- in contrast to major/minor merger remnants which peak at much greater amplitudes but are consistent with controls only $\sim1$ Gyr after coalescence. Integrating the boosts in SFRs and asymmetries driven by $μ\geq0.01$ mergers since $z=0.7$ in TNG50 SFGs, we show that mini mergers are responsible for (i) $55$ per cent of all merger-driven star formation and (ii) $70$ per cent of merger-driven asymmetric structure. Due to their relative frequency and prolonged boost timescales, mini mergers dominate over their minor and major counterparts in driving star formation and asymmetry in SFGs. △ Less

Submitted 7 October, 2023; v1 submitted 28 August, 2023; originally announced August 2023.

Comments: 31 pages; 19 figures; MNRAS accepted; Image data available via the TNG website: www.tng-project.org/bottrell23

Uncovering a Massive z~7.7 Galaxy Hosting a Heavily Obscured Radio-Loud QSO Candidate in COSMOS-Web

Authors: Erini Lambrides, Marco Chiaberge, Arianna Long, Daizhong Liu, Hollis B. Akins, Andrew F. Ptak, Irham Taufik Andika, Alessandro Capetti, Caitlin M. Casey, Jaclyn B. Champagne, Katherine Chworowsky, Tracy E. Clarke, Olivia R. Cooper, Xuheng Ding, Dillon Z. Dong, Andreas L. Faisst, Jordan Y. Forman, Maximilien Franco, Steven Gillman, Ghassem Gozaliasl, Kirsten R. Hall, Santosh Harish, Christopher C. Hayward, Michaela Hirschmann, Taylor A. Hutchison , et al. (25 additional authors not shown)

Abstract: In this letter, we report the discovery of the highest redshift, heavily obscured, radio-loud AGN candidate selected using JWST NIRCam/MIRI, mid-IR, sub-mm, and radio imaging in the COSMOS-Web field. Using multi-frequency radio observations and mid-IR photometry, we identify a powerful, radio-loud (RL), growing supermassive black hole (SMBH) with significant spectral steepening of the radio SED (… ▽ More In this letter, we report the discovery of the highest redshift, heavily obscured, radio-loud AGN candidate selected using JWST NIRCam/MIRI, mid-IR, sub-mm, and radio imaging in the COSMOS-Web field. Using multi-frequency radio observations and mid-IR photometry, we identify a powerful, radio-loud (RL), growing supermassive black hole (SMBH) with significant spectral steepening of the radio SED ($f_{1.28 \mathrm{GHz}} \sim 2$ mJy, $q_{24μm} = -1.1$, $α_{1.28-3\mathrm{GHz}}=-1.2$, $Δα= -0.4$). In conjunction with ALMA, deep ground-based observations, ancillary space-based data, and the unprecedented resolution and sensitivity of JWST, we find no evidence of AGN contribution to the UV/optical/NIR data and thus infer heavy amounts of obscuration (N$_{\mathrm{H}} > 10^{23}$ cm$^{-2}$). Using the wealth of deep UV to sub-mm photometric data, we report a singular solution photo-z of $z_\mathrm{phot}$ = 7.7$^{+0.4}_{-0.3}$ and estimate an extremely massive host-galaxy ($\log M_{\star} = 11.4 -12\,\mathrm{M}_{\odot}$) hosting a powerful, growing SMBH (L$_{\mathrm{Bol}} = 4-12 \times 10^{46}$ erg s$^{-1}$). This source represents the furthest known obscured RL AGN candidate, and its level of obscuration aligns with the most representative but observationally scarce population of AGN at these epochs. △ Less

Submitted 15 December, 2023; v1 submitted 24 August, 2023; originally announced August 2023.

Comments: Accepted to ApJL

COSMOS-Web: Intrinsically Luminous z$\gtrsim$10 Galaxy Candidates Test Early Stellar Mass Assembly

Authors: Caitlin M. Casey, Hollis B. Akins, Marko Shuntov, Olivier Ilbert, Louise Paquereau, Maximilien Franco, Christopher C. Hayward, Steven L. Finkelstein, Michael Boylan-Kolchin, Brant E. Robertson, Natalie Allen, Malte Brinch, Olivia R. Cooper, Xuheng Ding, Nicole E. Drakos, Andreas L. Faisst, Seiji Fujimoto, Steven Gillman, Santosh Harish, Michaela Hirschmann, Shuowen Jin, Jeyhan S. Kartaltepe, Anton M. Koekemoer, Vasily Kokorev, Daizhong Liu , et al. (17 additional authors not shown)

Abstract: We report the discovery of 15 exceptionally luminous $10\lesssim z\lesssim14$ candidate galaxies discovered in the first 0.28 deg$^2$ of JWST/NIRCam imaging from the COSMOS-Web Survey. These sources span rest-frame UV magnitudes of $-20.5>M_{\rm UV}>-22$, and thus constitute the most intrinsically luminous $z\gtrsim10$ candidates identified by JWST to-date. Selected via NIRCam imaging with Hubble… ▽ More We report the discovery of 15 exceptionally luminous $10\lesssim z\lesssim14$ candidate galaxies discovered in the first 0.28 deg$^2$ of JWST/NIRCam imaging from the COSMOS-Web Survey. These sources span rest-frame UV magnitudes of $-20.5>M_{\rm UV}>-22$, and thus constitute the most intrinsically luminous $z\gtrsim10$ candidates identified by JWST to-date. Selected via NIRCam imaging with Hubble ACS/F814W, deep ground-based observations corroborate their detection and help significantly constrain their photometric redshifts. We analyze their spectral energy distributions using multiple open-source codes and evaluate the probability of low-redshift solutions; we conclude that 12/15 (80%) are likely genuine $z\gtrsim10$ sources and 3/15 (20%) likely low-redshift contaminants. Three of our $z\sim12$ candidates push the limits of early stellar mass assembly: they have estimated stellar masses $\sim5\times10^{9}\,M_\odot$, implying an effective stellar baryon fraction of $ε_{\star}\sim0.2-0.5$, where $ε_{\star}\equiv M_{\star}/(f_{b}M_{halo})$. The assembly of such stellar reservoirs is made possible due to rapid, burst-driven star formation on timescales $<$100\,Myr where the star-formation rate may far outpace the growth of the underlying dark matter halos. This is supported by the similar volume densities inferred for $M_\star\sim10^{10}\,M_\odot$ galaxies relative to $M_\star\sim10^{9}\,M_\odot$ -- both about $10^{-6}$ Mpc$^{-3}$ -- implying they live in halos of comparable mass. At such high redshifts, the duty cycle for starbursts would be of order unity, which could cause the observed change in the shape of the UVLF from a double powerlaw to Schechter at $z\approx8$. Spectroscopic redshift confirmation and ensuing constraints of their masses will be critical to understanding how, and if, such early massive galaxies push the limits of galaxy formation in $Λ$CDM. △ Less

Submitted 21 August, 2023; originally announced August 2023.

Comments: 30 pages, 9 figures; ApJ submitted

Unveiling the distant Universe: Characterizing $z\ge9$ Galaxies in the first epoch of COSMOS-Web

Authors: Maximilien Franco, Hollis B. Akins, Caitlin M. Casey, Steven L. Finkelstein, Marko Shuntov, Katherine Chworowsky, Andreas L. Faisst, Seiji Fujimoto, Olivier Ilbert, Anton M. Koekemoer, Daizhong Liu, Christopher C. Lovell, Claudia Maraston, Henry Joy McCracken, Jed McKinney, Brant E. Robertson, Micaela B. Bagley, Jaclyn B. Champagne, Olivia R. Cooper, Xuheng Ding, Nicole E. Drakos, Andrea Enia, Steven Gillman, Christopher C. Hayward, Michaela Hirschmann , et al. (25 additional authors not shown)

Abstract: We report the identification of 15 galaxy candidates at $z\ge9$ using the initial COSMOS-Web JWST observations over 77 arcmin$^2$ through four NIRCam filters (F115W, F150W, F277W, F444W) with an overlap with MIRI (F770W) of 8.7 arcmin$^2$. We fit the sample using several publicly-available SED fitting and photometric redshift codes and determine their redshifts between $z=9.3$ and $z=10.9$ (… ▽ More We report the identification of 15 galaxy candidates at $z\ge9$ using the initial COSMOS-Web JWST observations over 77 arcmin$^2$ through four NIRCam filters (F115W, F150W, F277W, F444W) with an overlap with MIRI (F770W) of 8.7 arcmin$^2$. We fit the sample using several publicly-available SED fitting and photometric redshift codes and determine their redshifts between $z=9.3$ and $z=10.9$ ($\langle z\rangle=10.0$), UV-magnitudes between M$_{\rm UV}$ = $-$21.2 and $-$19.5 (with $\langle $M$_{\rm UV}\rangle=-20.2$) and rest-frame UV slopes ($\langle β\rangle=-2.4$). These galaxies are, on average, more luminous than most $z\ge9$ candidates discovered by JWST so far in the literature, while exhibiting similar blue colors in their rest-frame UV. The rest-frame UV slopes derived from SED-fitting are blue ($β\sim$[$-$2.0, $-$2.7]) without reaching extremely blue values as reported in other recent studies at these redshifts. The blue color is consistent with models that suggest the underlying stellar population is not yet fully enriched in metals like similarly luminous galaxies in the lower redshift Universe. The derived stellar masses with $\langle \log_{\rm 10} ($M$_\star/$M$_\odot)\rangle\approx8-9$ are not in tension with the standard $Λ$CDM model and our measurement of the volume density of such UV luminous galaxies aligns well with previously measured values presented in the literature at $z\sim9-10$. Our sample of galaxies, although compact, are significantly resolved. △ Less

Submitted 1 August, 2023; originally announced August 2023.

Comments: Submitted to ApJ

Resolving galactic-scale obscuration of X-ray AGN at $z\gtrsim1$ with COSMOS-Web

Authors: John D. Silverman, Vincenzo Mainieri, Xuheng Ding, Daizhong Liu, Knud Jahnke, Michaela Hirschmann, Jeyhan Kartaltepe, Erini Lambrides, Masafusa Onoue, Benny Trakhtenbrot, Eleni Vardoulaki, Angela Bongiorno, Caitlin Casey, Francesca Civano, Andreas Faisst, Maximilien Franco, Steven Gillman, Ghassem Gozaliasl, Christopher C. Hayward, Anton M. Koekemoer, Vasily Kokorev, Georgios Magdis, Stefano Marchesi, Robert Michael Rich, Martin Sparre , et al. (3 additional authors not shown)

Abstract: A large fraction of the accreting supermassive black hole population is shrouded by copious amounts of gas and dust, particularly in the distant ($z\gtrsim1$) Universe. While much of the obscuration is attributed to a parsec-scale torus, there is a known contribution from the larger-scale host galaxy. Using JWST/NIRCam imaging from the COSMOS-Web survey, we probe the galaxy-wide dust distribution… ▽ More A large fraction of the accreting supermassive black hole population is shrouded by copious amounts of gas and dust, particularly in the distant ($z\gtrsim1$) Universe. While much of the obscuration is attributed to a parsec-scale torus, there is a known contribution from the larger-scale host galaxy. Using JWST/NIRCam imaging from the COSMOS-Web survey, we probe the galaxy-wide dust distribution in X-ray selected AGN up to $z\sim2$. Here, we focus on a sample of three AGNs with their host galaxies exhibiting prominent dust lanes, potentially due to their edge-on alignment. These represent 27% (3 out of 11 with early NIRCam data) of the heavily obscured ($N_H>10^{23}$ cm$^{-2}$) AGN population. With limited signs of a central AGN in the optical and near-infrared, the NIRCam images are used to produce reddening maps $E(B-V)$ of the host galaxies. We compare the mean central value of $E(B-V)$ to the X-ray obscuring column density along the line-of-sight to the AGN ($N_H\sim10^{23-23.5}$ cm$^{-2}$). We find that the extinction due to the host galaxy is present ($0.6\lesssim E(B-V) \lesssim 0.9$; $1.9 \lesssim A_V \lesssim 2.8$) and significantly contributes to the X-ray obscuration at a level of $N_H\sim10^{22.5}$ cm$^{-2}$ assuming an SMC gas-to-dust ratio which amounts to $\lesssim$30% of the total obscuring column density. These early results, including three additional cases from CEERS, demonstrate the ability to resolve such dust structures with JWST and separate the different circumnuclear and galaxy-scale obscuring structures. △ Less

Submitted 5 June, 2023; originally announced June 2023.

Comments: 12 pages, 8 figures, Accepted for publication in ApJL

HYPERION. Interacting companion and outflow in the most luminous $z>6$ quasar

Authors: R. Tripodi, J. Scholtz, R. Maiolino, S. Fujimoto, S. Carniani, J. D. Silverman, C. Feruglio, M. Ginolfi, L. Zappacosta, T. Costa, G. C. Jones, E. Piconcelli, M. Bischetti, F. Fiore

Abstract: We present ALMA deep observations of the [CII] 158 $μ$m emission line and the continuum at 253 GHz and 99 GHz towards SDSS J0100+2802 at $z\simeq 6.3$, the most luminous quasi-stellar object (QSO) at z$>$6. J0100+2802 belongs to the HYPERION sample of luminous QSOs at $z\sim 6-7.5$. The observations (at 2.2" resolution in Band 3 and 0.9" resolution in Band 6) are optimized to detect extended emiss… ▽ More We present ALMA deep observations of the [CII] 158 $μ$m emission line and the continuum at 253 GHz and 99 GHz towards SDSS J0100+2802 at $z\simeq 6.3$, the most luminous quasi-stellar object (QSO) at z$>$6. J0100+2802 belongs to the HYPERION sample of luminous QSOs at $z\sim 6-7.5$. The observations (at 2.2" resolution in Band 3 and 0.9" resolution in Band 6) are optimized to detect extended emission around the QSO. We detect an interacting, tidally disrupted companion both in [CII], peaking at $z\sim 6.332$, and in continuum, stretching on scales up to 20 kpc from the quasar, with a knotty morphology. The higher velocity dispersion in the direction of the companion emission and the complex morphology of tidally stretched galaxy suggest a possible ongoing or future merger. For the newly-detected companion we derive a range of dust mass, $M_{\rm dust}=(0.3-2.6)\times 10^7\ \rm M_\odot$, and of star formation rate, SFR$=[35-344]\ \rm M_\odot$. This shows that both the QSO and its companion are gas rich and that a major merger may be at the origin of the boosted star formation. We also detect a broad blueshifted component in the [CII] spectrum, that we interpret as a gaseous outflow for which we estimate a mass outflow rate in the range $\dot{M}_{\rm out}=(118-269)\ \rm M_\odot\ yr^{-1}$. J0100+2802 was recently found to reside in a strong overdensity, however this close companion remained undetected by both previous higher resolution ALMA observations and by JWST-NIRCAM imaging. Our results highlight the importance of deep medium-resolution ALMA observations for the study of QSOs and their environment at the Epoch of Reionization. △ Less

Submitted 9 January, 2024; v1 submitted 2 June, 2023; originally announced June 2023.

Comments: 11 pages, 7 figures, 3 Tables. Accepted by A&A. In press

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

Morphological asymmetries of quasar host galaxies with Subaru Hyper Suprime-Cam

Authors: Shenli Tang, John D. Silverman, Hassen M. Yesuf, Xuheng Ding, Junyao Li, Connor Bottrell, Andy Goulding, Kiyoaki Christopher Omori, Yoshiki Toba, Toshihiro Kawaguchi

Abstract: How does the host galaxy morphology influence a central quasar or vice versa? We address this question by measuring the asymmetries of 2424 SDSS quasar hosts at $0.2<z<0.8$ using broad-band ($grizy$) images from the Hyper Suprime-Cam Subaru Strategic Program. Control galaxies (without quasars) are selected by matching the redshifts and stellar masses of the quasar hosts. A two-step pipeline is run… ▽ More How does the host galaxy morphology influence a central quasar or vice versa? We address this question by measuring the asymmetries of 2424 SDSS quasar hosts at $0.2<z<0.8$ using broad-band ($grizy$) images from the Hyper Suprime-Cam Subaru Strategic Program. Control galaxies (without quasars) are selected by matching the redshifts and stellar masses of the quasar hosts. A two-step pipeline is run to decompose the PSF and \sersic\ components, and then measure asymmetry indices ($A_{\rm CAS}$, $A_{\rm outer}$, and $A_{\rm shape}$) of each quasar host and control galaxy. We find a mild correlation between host asymmetry and AGN bolometric luminosity ($L_{\rm bol}$) for the full sample (spearman correlation of 0.37) while a stronger trend is evident at the highest luminosities ($L_{\rm bol}>45$). This then manifests itself into quasar hosts being more asymmetric, on average, when they harbor a more massive and highly accreting black hole. The merger fraction also positively correlates with $L_{\rm bol}$ and reaches up to 35\% for the most luminous. Compared to control galaxies, quasar hosts are marginally more asymmetric (excess of 0.017 in median at 9.4$σ$ level) and the merger fractions are similar ($\sim 16.5\%$). We quantify the dependence of asymmetry on optical band which demonstrates that mergers are more likely to be identified with the bluer bands and the correlation between $L_{\rm bol}$ and asymmetry is also stronger in such bands. We stress that the band dependence, indicative of a changing stellar population, is an important factor in considering the influence of mergers on AGN activity. △ Less

Submitted 12 April, 2023; originally announced April 2023.

Comments: 27 pages, 28 figures

Journal ref: MNRAS, volume 521, number 4, pages 5272--5297 (2023)

JWST and ALMA imaging of dust-obscured, massive substructures in a typical $z \sim 3$ star-forming disk galaxy

Authors: Wiphu Rujopakarn, Christina C. Williams, Emanuele Daddi, Malte Schramm, Fengwu Sun, Stacey Alberts, George H. Rieke, Qing-Hua Tan, Sandro Tacchella, Mauro Giavalisco, John D. Silverman

Abstract: We present an identification of dust-attenuated star-forming galactic-disk substructures in a typical star-forming galaxy (SFG), UDF2, at $z = 2.696$. To date, substructures containing significant buildup of stellar mass and actively forming stars have yet to be found in typical (i.e., main-sequence) SFGs at $z > 2$. This is due to the strong dust attenuation common in massive galaxies at the epoc… ▽ More We present an identification of dust-attenuated star-forming galactic-disk substructures in a typical star-forming galaxy (SFG), UDF2, at $z = 2.696$. To date, substructures containing significant buildup of stellar mass and actively forming stars have yet to be found in typical (i.e., main-sequence) SFGs at $z > 2$. This is due to the strong dust attenuation common in massive galaxies at the epoch and the scarcity of high-resolution, high-sensitivity extinction-independent imaging. To search for disk substructures, we subtracted the central stellar-mass disk from the JWST/NIRCam rest-frame 1.2 $μ$m image ($0.13''$ resolution) and subtracted, in the visibility plane, the central starburst disk from ALMA rest-frame 240 $μ$m observations ($0.03''$ resolution). The residual images revealed substructures at rest-frame 1.2 $μ$m co-located with those found at rest-frame 240 $μ$m, $\simeq 2$ kpc away from the galactic center. The largest substructure contains $\simeq20$% of the total stellar mass and $\simeq5$% of the total SFR of the galaxy. While UDF2 exhibits a kinematically-ordered velocity field of molecular gas consistent with a secularly evolving disk, more sensitive observations are required to characterize the nature and the origin of this substructure (spiral arms, minor merger, or other types of disk instabilities). UDF2 resides in an overdense region ($N \geqslant 4$ massive galaxies within 70 kpc projected distance at $z=2.690-2.697$) and the substructures may be associated with interaction-induced instabilities. Importantly, a statistical sample of such substructures identified with JWST and ALMA could play a key role in bridging the gap between the bulge-forming starburst and the rest of the galaxy. △ Less

Submitted 10 April, 2023; originally announced April 2023.

Comments: 7 pages, 5 figures; ApJL, accepted

The eROSITA Final Equatorial-Depth Survey (eFEDS): Host-galaxy Demographics of X-ray AGNs with Subaru Hyper Suprime-Cam

Authors: Junyao Li, John D. Silverman, Andrea Merloni, Mara Salvato, Johannes Buchner, Andy Goulding, Teng Liu, Riccardo Arcodia, Johan Comparat, Xuheng Ding, Kohei Ichikawa, Masatoshi Imanishi, Toshihiro Kawaguchi, Lalitwadee Kawinwanichakij, Yoshiki Toba

Abstract: We investigate the physical properties, such as star-forming activity, disk vs. bulge nature, galaxy size, and obscuration of 3796 X-ray selected AGNs at $0.2<z<0.8$ in the eFEDS field. Using Subaru Hyper Suprime-Cam imaging data in the $grizy$ bands for SRG/eROSITA-detected AGNs, we measure the structural parameters for AGN host galaxies by performing a 2D AGN-host image decomposition. We then co… ▽ More We investigate the physical properties, such as star-forming activity, disk vs. bulge nature, galaxy size, and obscuration of 3796 X-ray selected AGNs at $0.2<z<0.8$ in the eFEDS field. Using Subaru Hyper Suprime-Cam imaging data in the $grizy$ bands for SRG/eROSITA-detected AGNs, we measure the structural parameters for AGN host galaxies by performing a 2D AGN-host image decomposition. We then conduct spectral energy distribution fitting to derive stellar mass and rest-frame colors for AGN hosts. We find that (1) AGNs can contribute significantly to the total optical light down to ${\rm log}\,L_{\rm X}\sim 42.5\ \rm erg\,s^{-1}$, thus ignoring the AGN component can significantly bias the structural measurements; (2) AGN hosts are predominately star-forming galaxies at ${\rm log}\,\mathcal{M}_\star \lesssim 11.3\ M_\odot$; (3) the bulk of AGNs (64%) reside in galaxies with significant stellar disks, while their host galaxies become increasingly bulge dominated and quiescent at ${\rm log}\,\mathcal{M}_\star \gtrsim 11.0\ M_\odot$; (4) the size-stellar mass relation of AGN hosts tends to lie between that of inactive star-forming and quiescent galaxies, suggesting that the physical mechanism responsible for building the central stellar density also efficiently fuel the black hole growth; (5) the hosts of X-ray unobscured AGNs are biased towards face-on systems and the average $E(B-V)/N_{\rm H}$ is similar to the galactic dust-to-gas ratio, suggesting that some of the obscuration of the nuclei could come from galaxy-scale gas and dust, which may partly account for (up to 30%) the deficiency of star-forming disks as host galaxies for the most massive AGNs. These results are consistent with a scenario in which the black hole and galaxy grow in mass while transform in structure and star-forming activity, as desired to establish the local scaling relations. △ Less

Submitted 23 February, 2023; originally announced February 2023.

Comments: 15 pages, 13 figures, submitted to MNRAS

Detection of stellar light from quasar host galaxies at redshifts above 6

Authors: Xuheng Ding, Masafusa Onoue, John D. Silverman, Yoshiki Matsuoka, Takuma Izumi, Michael A. Strauss, Knud Jahnke, Camryn L. Phillips, Junyao Li, Marta Volonteri, Zoltan Haiman, Irham Taufik Andika, Kentaro Aoki, Shunsuke Baba, Rebekka Bieri, Sarah E. I. Bosman, Connor Bottrell, Anna-Christina Eilers, Seiji Fujimoto, Melanie Habouzit, Masatoshi Imanishi, Kohei Inayoshi, Kazushi Iwasawa, Nobunari Kashikawa, Toshihiro Kawaguchi , et al. (19 additional authors not shown)

Abstract: The detection of starlight from the host galaxies of quasars during the reionization epoch ($z>6$) has been elusive, even with deep HST observations. The current highest redshift quasar host detected, at $z=4.5$, required the magnifying effect of a foreground lensing galaxy. Low-luminosity quasars from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) mitigate the challenge of detecting the… ▽ More The detection of starlight from the host galaxies of quasars during the reionization epoch ($z>6$) has been elusive, even with deep HST observations. The current highest redshift quasar host detected, at $z=4.5$, required the magnifying effect of a foreground lensing galaxy. Low-luminosity quasars from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) mitigate the challenge of detecting their underlying, previously-undetected host galaxies. Here we report rest-frame optical images and spectroscopy of two HSC-SSP quasars at $z>6$ with JWST. Using NIRCam imaging at 3.6$μ$m and 1.5$μ$m and subtracting the light from the unresolved quasars, we find that the host galaxies are massive (stellar masses of $13\times$ and $3.4\times$ $10^{10}$ M$_{\odot}$, respectively), compact, and disk-like. NIRSpec medium-resolution spectroscopy shows stellar absorption lines in the more massive quasar, confirming the detection of the host. Velocity-broadened gas in the vicinity of these quasars enables measurements of their black hole masses ($1.4\times 10^9$ and $2.0\times$ $10^{8}$ M$_{\odot}$, respectively). Their location in the black hole mass - stellar mass plane is consistent with the distribution at low redshift, suggesting that the relation between black holes and their host galaxies was already in place less than a billion years after the Big Bang. △ Less

Submitted 23 June, 2023; v1 submitted 25 November, 2022; originally announced November 2022.

Comments: Matched to the published Nature version of the article. 27 pages, 4 main figures, 1 table, 6 supplementary figures, 2 supplementary table

COSMOS-Web: An Overview of the JWST Cosmic Origins Survey

Authors: Caitlin M. Casey, Jeyhan S. Kartaltepe, Nicole E. Drakos, Maximilien Franco, Santosh Harish, Louise Paquereau, Olivier Ilbert, Caitlin Rose, Isabella G. Cox, James W. Nightingale, Brant E. Robertson, John D. Silverman, Anton M. Koekemoer, Richard Massey, Henry Joy McCracken, Jason Rhodes, Hollis B. Akins, Aristeidis Amvrosiadis, Rafael C. Arango-Toro, Micaela B. Bagley, Angela Bongiorno, Peter L. Capak, Jaclyn B. Champagne, Nima Chartab, Oscar A. Chavez Ortiz , et al. (60 additional authors not shown)

Abstract: We present the survey design, implementation, and outlook for COSMOS-Web, a 255 hour treasury program conducted by the James Webb Space Telescope in its first cycle of observations. COSMOS-Web is a contiguous 0.54 deg$^2$ NIRCam imaging survey in four filters (F115W, F150W, F277W, and F444W) that will reach 5$σ$ point source depths ranging $\sim$27.5-28.2 magnitudes. In parallel, we will obtain 0.… ▽ More We present the survey design, implementation, and outlook for COSMOS-Web, a 255 hour treasury program conducted by the James Webb Space Telescope in its first cycle of observations. COSMOS-Web is a contiguous 0.54 deg$^2$ NIRCam imaging survey in four filters (F115W, F150W, F277W, and F444W) that will reach 5$σ$ point source depths ranging $\sim$27.5-28.2 magnitudes. In parallel, we will obtain 0.19 deg$^2$ of MIRI imaging in one filter (F770W) reaching 5$σ$ point source depths of $\sim$25.3-26.0 magnitudes. COSMOS-Web will build on the rich heritage of multiwavelength observations and data products available in the COSMOS field. The design of COSMOS-Web is motivated by three primary science goals: (1) to discover thousands of galaxies in the Epoch of Reionization ($6<z<11$) and map reionization's spatial distribution, environments, and drivers on scales sufficiently large to mitigate cosmic variance, (2) to identify hundreds of rare quiescent galaxies at $z>4$ and place constraints on the formation of the Universe's most massive galaxies ($M_\star>10^{10}$\,M$_\odot$), and (3) directly measure the evolution of the stellar mass to halo mass relation using weak gravitational lensing out to $z\sim2.5$ and measure its variance with galaxies' star formation histories and morphologies. In addition, we anticipate COSMOS-Web's legacy value to reach far beyond these scientific goals, touching many other areas of astrophysics, such as the identification of the first direct collapse black hole candidates, ultracool sub-dwarf stars in the Galactic halo, and possibly the identification of $z>10$ pair-instability supernovae. In this paper we provide an overview of the survey's key measurements, specifications, goals, and prospects for new discovery. △ Less

Submitted 8 March, 2023; v1 submitted 14 November, 2022; originally announced November 2022.

Comments: 46 pages, 16 figures, ApJ accepted

A machine learning approach to assessing the presence of substructure in quasar host galaxies using the Hyper Suprime-Cam Subaru Strategic Program

Authors: Chris Nagele, John D. Silverman, Tilman Hartwig, Junyao Li, Connor Bottrell, Xuheng Ding, Yoshiki Toba

Abstract: The conditions under which galactic nuclear regions become active are largely unknown, although it has been hypothesized that secular processes related to galaxy morphology could play a significant role. We investigate this question using optical i-band images of 3096 SDSS quasars and galaxies at 0.3<z<0.6 from the Hyper Suprime-Cam Subaru Strategic Program, which possess a unique combination of a… ▽ More The conditions under which galactic nuclear regions become active are largely unknown, although it has been hypothesized that secular processes related to galaxy morphology could play a significant role. We investigate this question using optical i-band images of 3096 SDSS quasars and galaxies at 0.3<z<0.6 from the Hyper Suprime-Cam Subaru Strategic Program, which possess a unique combination of area, depth and resolution, allowing the use of residual images, after removal of the quasar and smooth galaxy model, to investigate internal structural features. We employ a variational auto-encoder which is a generative model that acts as a form of dimensionality reduction. We analyze the lower dimensional latent space in search of features which correlate with nuclear activity. We find that the latent space does separate images based on the presence of nuclear activity which appears to be associated with more pronounced components (i.e., arcs, rings and bars) as compared to a matched control sample of inactive galaxies. These results suggest the importance of secular processes, and possibly mergers (by their remnant features) in activating or sustaining black hole growth. Our study highlights the breadth of information available in ground-based imaging taken under optimal seeing conditions and having accurate characterization of the point spread function (PSF) thus demonstrating future science to come from the Rubin Observatory. △ Less

Submitted 17 April, 2023; v1 submitted 21 October, 2022; originally announced October 2022.

Journal ref: 2023 ApJ 947 30

Opening the era of quasar host studies at high redshift with JWST

Authors: Xuheng Ding, John D. Silverman, Masafusa Onoue

Abstract: We measure the host galaxy properties of five quasars with $z\sim 1.6 - 3.5$ selected from SDSS and AEGIS, which fall within the JWST/HST CEERS survey area. A PSF library is constructed based on stars in the full field-of-view of the data and used with the 2-dimensional image modeling tool galight to decompose the quasar and its host with multi-band filters available for HST ACS+WFC3 and JWST NIRC… ▽ More We measure the host galaxy properties of five quasars with $z\sim 1.6 - 3.5$ selected from SDSS and AEGIS, which fall within the JWST/HST CEERS survey area. A PSF library is constructed based on stars in the full field-of-view of the data and used with the 2-dimensional image modeling tool galight to decompose the quasar and its host with multi-band filters available for HST ACS+WFC3 and JWST NIRCAM (12 filters covering HST F606W to JWST F444W). As demonstrated, JWST provides the first capability to detect quasar hosts at $z>3$ and enables spatially-resolved studies of the underlying stellar populations at $z\sim2$ within morphological structures (spiral arms, bar) not possible with HST. Overall, we find quasar hosts to be disk-like, lack merger signatures, and have sizes generally more compact than typical star-forming galaxies at their respective stellar mass, thus in agreement with results at lower redshifts. The fortuitous face-on orientation of SDSSJ1420+5300A at $z = 1.646$ enables us to find higher star formation and younger ages in the central $2-4$ kpc region relative to the outskirts, which may help explain the relatively compact nature of quasar hosts and pose a challenge to AGN feedback models. △ Less

Submitted 18 October, 2022; v1 submitted 7 September, 2022; originally announced September 2022.

Comments: 12 pages, 1 table, 6 figures; ApJL accepted; Comments welcome

Journal ref: 2022, ApJL, 939, L28

Concordance between observations and simulations in the evolution of the mass relation between supermassive black holes and their host galaxies

Authors: Xuheng Ding, John D. Silverman, Tommaso Treu, Junyao Li, Aklant K. Bhowmick, Nicola Menci, Marta Volonteri, Laura Blecha, Tiziana Di Matteo, Yohan Dubois

Abstract: We carry out a comparative analysis of the relation between the mass of supermassive black holes (BHs) and the stellar mass of their host galaxies at $0.2<z<1.7$ using well-matched observations and multiple state-of-the-art simulations (e.g., Massive Black II, Horizon-AGN, Illustris, TNG and a semi-analytic model). The observed sample consists of 646 uniformly-selected SDSS quasars (… ▽ More We carry out a comparative analysis of the relation between the mass of supermassive black holes (BHs) and the stellar mass of their host galaxies at $0.2<z<1.7$ using well-matched observations and multiple state-of-the-art simulations (e.g., Massive Black II, Horizon-AGN, Illustris, TNG and a semi-analytic model). The observed sample consists of 646 uniformly-selected SDSS quasars ($0.2 < z < 0.8$) and 32 broad-line active galactic nuclei (AGNs; $1.2<z<1.7$) with imaging from Hyper Suprime-Cam (HSC) for the former and Hubble Space Telescope (HST) for the latter. We first add realistic observational uncertainties to the simulation data and then construct a simulated sample in the same manner as the observations. Over the full redshift range, our analysis demonstrates that all simulations predict a level of intrinsic scatter of the scaling relations comparable to the observations which appear to agree with the dispersion of the local relation. Regarding the mean relation, Horizon-AGN and TNG are in closest agreement with the observations at low and high redshift ($z\sim$ 0.2 and 1.5, respectively) while the other simulations show subtle differences within the uncertainties. For insight into the physics involved, the scatter of the scaling relation, seen in the SAM, is reduced by a factor of two and closer to the observations after adopting a new feedback model that considers the geometry of the AGN outflow. The consistency in the dispersion with redshift in our analysis supports the importance of both quasar- and radio-mode feedback prescriptions in the simulations. Finally, we highlight the importance of increasing the sensitivity (e.g., using the James Webb Space Telescope), thereby pushing to lower masses and minimizing biases due to selection effects. △ Less

Submitted 21 July, 2022; v1 submitted 9 May, 2022; originally announced May 2022.

Comments: 17 pages, 9 figures, Published by ApJ

Journal ref: 2022, ApJ, 933, 132D

Dark matter halos of luminous AGNs from galaxy-galaxy lensing with the HSC Subaru Strategic Program

Authors: Wentao Luo, John D. Silverman, Surhud More, Andy Goulding, Hironao Miyatake, Takahiro Nishimichi, Chiaki Hikage, Lalitwadee Kawinwanichakij, Junyao Li, Xiangchong Li, Elinor Medezinski, Masamune Oguri, Taira Oogi, Cristobal Sifon

Abstract: We assess the dark matter halo masses of luminous AGNs over the redshift range 0.2 to 1.2 using galaxy-galaxy lensing based on imaging data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). We measure the weak lensing signal of a sample of 8882 AGNs constructed using HSC and WISE photometry. The lensing detection around AGNs has a signal-to-noise ratio of 15. As expected, we find that… ▽ More We assess the dark matter halo masses of luminous AGNs over the redshift range 0.2 to 1.2 using galaxy-galaxy lensing based on imaging data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). We measure the weak lensing signal of a sample of 8882 AGNs constructed using HSC and WISE photometry. The lensing detection around AGNs has a signal-to-noise ratio of 15. As expected, we find that the lensing mass profile is consistent with that of massive galaxies ($M_{*}\sim 10.8~M_\odot$). Surprisingly, the lensing signal remains unchanged when the AGN sample is split into low and high stellar mass hosts. Specifically, we find that the excess surface density (ESD) of AGNs, residing in galaxies with high stellar masses, significantly differs from that of the control sample. We further fit a halo occupation distribution model to the data to infer the posterior distribution of parameters including the average halo mass. We find that the characteristic halo mass of the full AGN population lies near the knee ($\rm log(M_h/h^{-1}M_{\odot})=12.0$) of the stellar-to-halo mass relation (SHMR). Illustrative of the results given above, the halo masses of AGNs residing in host galaxies with high stellar masses (i.e., above the knee of the SHMR) falls below the calibrated SHMR while the halo mass of the low stellar mass sample is more consistent with the established SHMR. These results indicate that massive halos with higher clustering bias tend to suppress AGN activity, probably due to the lack of available gas. △ Less

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

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

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

The ALPINE-ALMA [CII] Survey: the population of [CII]-undetected galaxies and their role in the $\mathrm{L_{[CII]}}$-SFR relation

Authors: Michael Romano, L. Morselli, P. Cassata, M. Ginolfi, D. Schaerer, M. Béthermin, P. Capak, A. Faisst, O. Le Fèvre, J. D. Silverman, L. Yan, S. Bardelli, M. Boquien, M. Dessauges-Zavadsky, S. Fujimoto, N. P. Hathi, G. C. Jones, A. M. Koekemoer, B. C. Lemaux, H. Méndez-Hernández, D. Narayanan, M. Talia, D. Vergani, G. Zamorani, E. Zucca

Abstract: The [CII] 158$~μ$m emission line represents so far one of the most profitable tools for the investigation of the high-redshift galaxies in the early Universe. Being one of the brightest cooling lines in the rest-frame far-infrared regime of star-forming galaxies, it has been successfully exploited as a tracer of star-formation rate (SFR) in local sources. The picture is more complex at higher reds… ▽ More The [CII] 158$~μ$m emission line represents so far one of the most profitable tools for the investigation of the high-redshift galaxies in the early Universe. Being one of the brightest cooling lines in the rest-frame far-infrared regime of star-forming galaxies, it has been successfully exploited as a tracer of star-formation rate (SFR) in local sources. The picture is more complex at higher redshifts, where its usability in this context is still under investigation. Recent results from the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) survey suggest that there is no (or weak) evolution of the L$\mathrm{_{[CII]}}$-SFR relation up to $z\sim6$ but their reliability is hampered by the presence of a large population of [CII] non-detected galaxies. In this work, we characterize the population of [CII] non-detections in ALPINE. By stacking their ALMA spectra, we obtain a signal detected at $\sim5.1σ$, resulting in a [CII] luminosity of $\mathrm{log(L_\mathrm{[CII]}}/\mathrm{L_{\odot}})$ $\sim7.8$. When combining this value with those from the [CII] detections, we find a $\mathrm{L_{[CII]}}$-SFR relation with a slope $b=1.14\pm0.11$, in agreement within the uncertainties both with the linear relation found in the local Universe, and with the previous findings from ALPINE at $z\sim5$. This suggests that the [CII] line can be considered a good tracer of star formation up to the distant Universe. Finally, we show that the galaxies of our sample that most deviate from the observed L$_\mathrm{[CII]}$-SFR relation could suffer from a less precise redshift estimation, perhaps artificially reducing their [CII] luminosity. In this respect, we claim that there is no evidence in favour of a deficit of [CII] content in high-z galaxies, in contrast with earlier studies. △ Less

Submitted 4 February, 2022; originally announced February 2022.

Comments: Accepted for publication in A&A. 10 pages, 5 figures

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

Scale Reliant Inference

Authors: Michelle Pistner Nixon, Kyle C. McGovern, Jeffrey Letourneau, Lawrence A. David, Nicole A. Lazar, Sayan Mukherjee, Justin D. Silverman

Abstract: Scientific fields such as genomics, ecology, and political science often collect multivariate count data. In these fields, the data are often sufficiently noisy such that inferences regarding the total size of the measured systems have substantial uncertainty. This uncertainty can hinder downstream analyses, such as differential analysis in case-control studies. There have historically been two ap… ▽ More Scientific fields such as genomics, ecology, and political science often collect multivariate count data. In these fields, the data are often sufficiently noisy such that inferences regarding the total size of the measured systems have substantial uncertainty. This uncertainty can hinder downstream analyses, such as differential analysis in case-control studies. There have historically been two approaches to this problem: one considers the data as compositional and the other as counts that can be normalized. In this article, we use the framework of partially identified models to rigorously study the types of scientific questions (estimands) that can be answered (estimated) using these data. We prove that satisfying Frequentist inferential criteria is impossible for many estimation problems. In contrast, we find that the criteria for Bayesian inference can be satisfied, yet it requires a particular type of model called a Bayesian partially identified model. We introduce Scale Simulation Random Variables as a flexible and computationally efficient form of Bayesian partially identified models for analyzing these data. We use simulations and data analysis to validate our theory. △ Less

Submitted 5 April, 2024; v1 submitted 10 January, 2022; originally announced January 2022.

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

Galaxy shapes of Light (GaLight): a 2D modeling of galaxy images

Authors: Xuheng Ding, Simon Birrer, Tommaso Treu, John D. Silverman

Abstract: Galight is a Python-based open-source package that can be used to perform two-dimensional model fitting of optical and near-infrared images to characterize the light distribution of galaxies with components including a disk, bulge, bar, and quasar. The decomposition of stellar components has been demonstrated in published studies of inactive galaxies and quasar host galaxies observed by the Hubble… ▽ More Galight is a Python-based open-source package that can be used to perform two-dimensional model fitting of optical and near-infrared images to characterize the light distribution of galaxies with components including a disk, bulge, bar, and quasar. The decomposition of stellar components has been demonstrated in published studies of inactive galaxies and quasar host galaxies observed by the Hubble Space Telescope and Subaru's Hyper Suprime-Cam. Galight utilizes the image modeling capabilities of lenstronomy while redesigning the user interface for the analysis of large samples of extragalactic sources. The package is user-friendly with some automatic features such as determining the cutout size of the modeling frame, searching for PSF-stars in field-of-view, estimating the noise map of the data, identifying all the objects to set the initial model, and associated parameters to fit them simultaneously. These features minimize the manpower and allow the automatic fitting tasks. The software is distributed under the MIT license. The source code, installation guidelines, and example notebooks code can be found at https://galight.readthedocs.io/en/latest/ △ Less

Submitted 16 November, 2021; originally announced November 2021.

Comments: 10 pages, 5 figures. This manuscript is posted as a user manual and will be constantly updated in the future with more upcoming features

Hyper Suprime-Cam Subaru Strategic Program: A Mass-Dependent Slope of the Galaxy Size-Mass Relation at $z<1$

Authors: Lalitwadee Kawinwanichakij, John D. Silverman, Xuheng Ding, Angelo George, Ivana Damjanov, Marcin Sawicki, Masayuki Tanaka, Dan S. Taranu, Simon Birrer, Song Huang, Junyao Li, Masato Onodera, Takatoshi Shibuya, Naoki Yasuda

Abstract: We present the galaxy size-mass ($R_{e}-M_{\ast}$) distributions using a stellar-mass complete sample of $\sim1.5$ million galaxies, covering $\sim100$ deg$^2$, with $\log(M_{\ast}/M_{\odot})>10.2~(9.2)$ over the redshift range $0.2<z<1.0$ $(z<0.6)$ from the second public data release of the Hyper Suprime-Cam Subaru Strategic Program. We confirm that, at fixed redshift and stellar mass over the ra… ▽ More We present the galaxy size-mass ($R_{e}-M_{\ast}$) distributions using a stellar-mass complete sample of $\sim1.5$ million galaxies, covering $\sim100$ deg$^2$, with $\log(M_{\ast}/M_{\odot})>10.2~(9.2)$ over the redshift range $0.2<z<1.0$ $(z<0.6)$ from the second public data release of the Hyper Suprime-Cam Subaru Strategic Program. We confirm that, at fixed redshift and stellar mass over the range of $\log(M_{\ast}/M_{\odot})<11$, star-forming galaxies are on average larger than quiescent galaxies. The large sample of galaxies with accurate size measurements, thanks to the excellent imaging quality, also enables us to demonstrate that the $R_{e}-M_{\ast}$ relations of both populations have a form of broken power-law, with a clear change of slopes at a pivot stellar mass $M_{p}$. For quiescent galaxies, below an (evolving) pivot mass of $\log(M_{p}/M_{\odot})=10.2-10.6$ the relation follows $R_{e}\propto M_{\ast}^{0.1}$; above $M_{p}$ the relation is steeper and follows $R_{e}\propto M_{\ast}^{0.6-0.7}$. For star-forming galaxies, below $\log(M_{p}/M_{\odot})\sim10.7$ the relation follows $R_{e}\propto M_{\ast}^{0.2}$; above $M_{p}$ the relation evolves with redshift and follows $R_{e}\propto M_{\ast}^{0.3-0.6}$. The shallow power-law slope for quiescent galaxies below $M_{p}$ indicates that large low-mass quiescent galaxies have sizes similar to those of their counterpart star-forming galaxies. We take this as evidence that large low-mass quiescent galaxies have been recently quenched (presumably through environment-specific process) without significant structural transformation. Interestingly, the pivot stellar mass of the $R_{e}-M_{\ast}$ relations coincides with mass at which half of the galaxy population is quiescent, implied that the pivot mass represents the transition of galaxy growth from being dominated by in-situ star formation to being dominated by (dry) mergers. △ Less

Submitted 20 September, 2021; originally announced September 2021.

Comments: Accepted for publication in ApJ. 45 pages, 17 figures

The stellar mass versus stellar metallicity relation of star-forming galaxies at $1.6\le z\le3.0$ and implications for the evolution of the $α$-enhancement

Authors: Daichi Kashino, Simon J. Lilly, Alvio Renzini, Emanuele Daddi, Giovanni Zamorani, John D. Silverman, Olivier Ilbert, Yingjie Peng, Vincenzo Mainieri, Sandro Bardelli, Elena Zucca, Jeyhan S. Kartaltepe, David B. Sanders

Abstract: We measure the relationship between stellar mass and stellar metallicity, the stellar mass--metallicity relation (MZR), for 1336 star-forming galaxies at $1.6\le z\le3.0$ (<z>=2.2) using rest-frame far-ultraviolet spectra from the zCOSMOS-deep survey. High signal-to-noise composite spectra containing stellar absorption features are fit with population synthesis model spectra of a range of metallic… ▽ More We measure the relationship between stellar mass and stellar metallicity, the stellar mass--metallicity relation (MZR), for 1336 star-forming galaxies at $1.6\le z\le3.0$ (<z>=2.2) using rest-frame far-ultraviolet spectra from the zCOSMOS-deep survey. High signal-to-noise composite spectra containing stellar absorption features are fit with population synthesis model spectra of a range of metallicity. We find stellar metallicities, which mostly reflect iron abundances, scaling as $(Z_{Fe,\ast}/Z_{Fe,\odot})=-(0.81\pm0.01)+(0.32+0.03)\log(M_\ast/10^{10}M_\odot)$ across the mass range of $10^9\lesssim M_\ast/M_\odot\lesssim10^{11}$, being $\approx6\times$ lower than seen locally at the same masses. The instantaneous oxygen-to-iron ratio ($α$-enhancement) inferred using the gas-phase oxygen MZRs, is on average found to be [O/Fe]$\approx0.47$, being higher than the local [O/Fe]$\approx0$. The observed changes in [O/Fe] and [Fe/H] are reproduced in simple flow-through gas-regulator models with steady star-formation histories (SFHs) that follow the evolving main sequence. Our models show that the [O/Fe] is determined almost entirely by the instantaneous specific star formation rate alone while being independent of the SFHs, mass, and the gas-regulation characteristics of the systems. We find that the locations of $\sim10^{10}M_\odot$ galaxies at z~2 in the [O/Fe]--metallicity planes are in remarkable agreement with the sequence of low-metallicity thick-disk stars in our Galaxy. This manifests a beautiful concordance between the results of Galactic archaeology and observations of high-redshift Milky Way progenitors. However, there remains a question of how and when the old metal-rich, low-$α$/Fe stars seen in the bulge had formed by z~2 because such a stellar population is not seen in our data and difficult to explain in the context of our models. △ Less

Submitted 13 September, 2021; originally announced September 2021.

Comments: Submitted to ApJ; 30 pages, 20 figures

The effect of active galactic nuclei on the cold interstellar medium in distant star-forming galaxies

Authors: Francesco Valentino, Emanuele Daddi, Annagrazia Puglisi, Georgios E. Magdis, Vasily Kokorev, Daizhong Liu, Suzanne C. Madden, Carlos Gomez-Guijarro, Min-Young Lee, Isabella Cortzen, Chiara Circosta, Ivan Delvecchio, James R. Mullaney, Yu Gao, Raphael Gobat, Manuel Aravena, Shuowen Jin, Seiji Fujimoto, John D. Silverman, Helmut Dannerbauer

Abstract: In the framework of a systematic ALMA study of IR-selected main-sequence and starburst galaxies at z~1-1.7 at typical ~1" resolution, we report on the effects of mid-IR- and X-ray-detected active galactic nuclei (AGN) on the reservoirs and excitation of molecular gas in a sample of 55 objects. We find detectable nuclear activity in ~30% of the sample. The presence of dusty tori influences the IR S… ▽ More In the framework of a systematic ALMA study of IR-selected main-sequence and starburst galaxies at z~1-1.7 at typical ~1" resolution, we report on the effects of mid-IR- and X-ray-detected active galactic nuclei (AGN) on the reservoirs and excitation of molecular gas in a sample of 55 objects. We find detectable nuclear activity in ~30% of the sample. The presence of dusty tori influences the IR SED of galaxies, as highlighted by the strong correlation among the AGN contribution to the total IR luminosity budget (fAGN = LIR,AGN/LIR), its hard X-ray emission, and the Rayleigh-Jeans to mid-IR (S1.2mm/S24um) observed color, with consequences on the empirical SFR estimates. Nevertheless, we find only marginal effects of AGN on the CO (J=2,4,5,7) or neutral carbon ([CI](1-0), [CI](2-1)) line luminosities and on the derived molecular gas excitation as gauged by line ratios and the full SLEDs. The [CI] and CO emission up to J=5,7 thus primarily traces the properties of the host in typical IR luminous galaxies. However, we highlight the existence of a large variety of line luminosities and ratios despite the homogeneous selection. In particular, we find a sparse group of AGN-dominated sources with the highest LIR,AGN/LIR,SFR ratios, >3, that are more luminous in CO(5-4) than what is predicted by the L'CO(5-4)-LIR,SFR relation, which might be the result of the nuclear activity. For the general population, our findings translate into AGN having minimal effects on quantities such as gas and dust fractions and SFEs. If anything, we find hints of a marginal tendency of AGN hosts to be compact at far-IR wavelengths and to display 1.8x larger dust optical depths. In general, this is consistent with a marginal impact of the nuclear activity on the gas reservoirs and star formation in average star-forming AGN hosts with LIR>5e11 Lsun, typically underrepresented in surveys of quasars and SMGs. △ Less

Submitted 8 September, 2021; originally announced September 2021.

Comments: 16 pages, 10 figures + Appendix, accepted in A&A on August 28th, 2021. The data compilation will be available on Vizier. (Minor) updates for the AGN population with respect to the data release of Valentino et al. 2020, A&A, 641, 155. For an early access, please contact the corresponding author. Abstract slightly modified to adjust to arXiv's requirements

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

Synchronized Co-evolution between Supermassive Black Holes and Galaxies Over the Last Seven Billion Years as Revealed by the Hyper Suprime-Cam

Authors: Junyao Li, John D. Silverman, Xuheng Ding, Michael A. Strauss, Andy Goulding, Malte Schramm, Hassen M. Yesuf, Mouyuan Sun, Yongquan Xue, Simon Birrer, Jingjing Shi, Yoshiki Toba, Tohru Nagao, Masatoshi Imanishi

Abstract: We measure the evolution of the $M_{\rm BH}-M_*$ relation using 584 uniformly-selected SDSS quasars at $0.2<z<0.8$. The black-hole masses ($M_{\rm BH}$) are derived from the single-epoch virial mass estimator using the H$β$ emission line, and span the range $7.0<{\rm log}\,M_{\rm BH}/M_\odot<9.5$. The host-galaxy stellar masses ($M_*$), which cover the interval… ▽ More We measure the evolution of the $M_{\rm BH}-M_*$ relation using 584 uniformly-selected SDSS quasars at $0.2<z<0.8$. The black-hole masses ($M_{\rm BH}$) are derived from the single-epoch virial mass estimator using the H$β$ emission line, and span the range $7.0<{\rm log}\,M_{\rm BH}/M_\odot<9.5$. The host-galaxy stellar masses ($M_*$), which cover the interval $10.0<{\rm log}\,M_* / M_\odot <11.5$, are determined by performing two-dimensional quasar-host decomposition of the Hyper Suprime-Cam images and spectral energy distribution fitting. To quantify sample-selection biases and measurement uncertainties on the mass terms, a mock quasar sample is constructed to jointly constrain the redshift evolution of the $M_{\rm BH}-M_*$ relation and its intrinsic scatter ($σ_μ$) through forward modeling. We find that the level of evolution is degenerate with $σ_μ$, such that both positive mild evolution (i.e, $M_{\rm BH}/M_*$ increases with redshift) with a small $σ_μ$, and negative mild evolution with a larger $σ_μ$ are consistent with our data. The posterior distribution of $σ_μ$ enables us to put a strong constraint on the intrinsic scatter of the $M_{\rm BH}-M_*$ relation, which has a best inference of $0.25_{-0.04}^{+0.03}$ dex, consistent with the local value. The redshift evolution of the $M_{\rm BH}-M_*$ relation relative to the local relation is constrained to be $(1+z)^{0.12_{-0.27}^{+0.28}}$, in agreement with no significant evolution since $z\sim0.8$. The tight and non-evolving $M_{\rm BH}-M_*$ relation is suggestive of a coupling through AGN feedback or/and a common gas supply at work, thus restricting the mass ratio of galaxies and their black holes to a limited range. Given the considerable stellar disk component, the $M_{\rm BH}-M_{\rm bulge}$ relation may evolve as previously seen at higher redshifts. △ Less

Submitted 6 September, 2021; originally announced September 2021.

Comments: 17 pages, 6 figures, accepted by ApJ

The ALPINE-ALMA [CII] survey: The contribution of major mergers to the galaxy mass assembly at z~5

Authors: M. Romano, P. Cassata, L. Morselli, G. C. Jones, M. Ginolfi, A. Zanella, M. Béthermin, P. Capak, A. Faisst, O. Le Fèvre, D. Schaerer, J. D. Silverman, L. Yan, S. Bardelli, M. Boquien, A. Cimatti, M. Dessauges-Zavadsky, A. Enia, S. Fujimoto, C. Gruppioni, N. P. Hathi, E. Ibar, A. M. Koekemoer, B. C. Lemaux, G. Rodighiero , et al. (3 additional authors not shown)

Abstract: Galaxy mergers are thought to be one of the main mechanisms of the mass assembly of galaxies. Recently, many works have suggested a possible increase in the fraction of major mergers in the early Universe, reviving the debate on which processes (e.g., cold accretion, star formation, mergers) most contribute to the mass build-up of galaxies through cosmic time. To estimate the importance of major m… ▽ More Galaxy mergers are thought to be one of the main mechanisms of the mass assembly of galaxies. Recently, many works have suggested a possible increase in the fraction of major mergers in the early Universe, reviving the debate on which processes (e.g., cold accretion, star formation, mergers) most contribute to the mass build-up of galaxies through cosmic time. To estimate the importance of major mergers in this context, we make use of the new data collected by the ALMA Large Program to INvestigate [CII] at Early times (ALPINE), which observed the [CII] 158 $μ$m emission line from a sample of 75 main-sequence star-forming galaxies at 4.4 < z < 5.9. We used, for the first time, the morpho-kinematic information provided by the [CII] emission to obtain the fraction of major mergers ($f_{MM}$) at z~5. By adopting different prescriptions for the merger timescales ($T_{MM}$), we converted this fraction into the merger rate per galaxy ($R_{MM}$) and per volume ($Γ_{MM}$). We then combined our results with those at lower redshifts from the literature, computing the cosmic evolution of the merger fraction. This is described by a rapid increase from z~0 to higher redshifts, a peak at z~3, and a slow decrease towards earlier epochs. Depending on the timescale prescription used, this fraction translates into a merger rate ranging between ~0.1 and ~4.0 Gyr$^{-1}$ at z~5. Finally, we compare the specific star formation and star-formation rate density with the analogous quantities from major mergers. Our new ALPINE data reveal the presence of a significant merging activity in the early Universe. However, whether this population of mergers can provide a relevant contribution to the galaxy mass assembly at these redshifts and through the cosmic epochs is strongly dependent on the assumption of the merger timescale. △ Less

Submitted 22 July, 2021; originally announced July 2021.

Comments: Accepted for publication in section 4. Extragalactic astronomy of Astronomy and Astrophysics

Journal ref: A&A 653, A111 (2021)

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): The first archetypal Quasar in the feedback phase discovered by eROSITA

Authors: M. Brusa, T. Urrutia, Y. Toba, J. Buchner, J. -Y. Li, T. Liu, M. Perna, M. Salvato, A. Merloni, B. Musiimenta, K. Nandra, J. Wolf, R. Arcodia, T. Dwelly, A. Georgakakis, A. Goulding, Y. Matsuoka, T. Nagao, M. Schramm, J. D. Silverman, Y. Terashima

Abstract: Theoretical models of galaxy-AGN co-evolution ascribe an important role for the feedback process to a short, luminous, obscured, and dust-enshrouded phase during which the accretion rate of the SMBH is expected to be at its maximum and the associated AGN-driven winds are also predicted to be maximally developed. To test this scenario, we have isolated a text-book candidate from the eROSITA Final E… ▽ More Theoretical models of galaxy-AGN co-evolution ascribe an important role for the feedback process to a short, luminous, obscured, and dust-enshrouded phase during which the accretion rate of the SMBH is expected to be at its maximum and the associated AGN-driven winds are also predicted to be maximally developed. To test this scenario, we have isolated a text-book candidate from the eROSITA Final Equatorial-Depth Survey (eFEDS) obtained within the Performance and Verification program of the eROSITA telescope on board Spectrum Roentgen Gamma. From an initial catalog of 246 hard X-ray selected sources matched with the photometric and spectroscopic information available within the eROSITA and Hyper Suprime-Cam consortia, three candidates Quasars in the feedback phase have been isolated applying the diagnostic proposed in Brusa et al. (2015). Only one source (eFEDSU J091157.5+014327) has a spectrum already available (from SDSS-DR16, z=0.603) and it unambiguously shows the presence of a broad component (FWHM~1650 km/s) in the [OIII]5007 line. The associated observed L_[OIII] is ~2.6x10^{42} erg/s, one to two orders of magnitude larger than that observed in local Seyferts and comparable to those observed in a sample of z~0.5 Type 1 Quasars. From the multiwavelength data available we derive an Eddington Ratio (L_bol/L_Edd) of ~0.25, and a bolometric correction in the hard X-ray of k_bol~10, lower than those observed for objects at similar bolometric luminosity. The presence of an outflow, the high X-ray luminosity and moderate X-ray obscuration (L_X~10^44.8 erg/s, N_H~2.7x10^22 cm^-2) and the red optical color, all match the prediction of quasars in the feedback phase from merger driven models. Forecasting to the full eROSITA all-sky survey with its spectroscopic follow-up, we predict that by the end of 2024 we will have a sample of few hundreds such objects at z=0.5-2. △ Less

Submitted 28 June, 2021; originally announced June 2021.

Comments: 11 pages, 6 figures, 2 tables. Resubmitted to A&A for the Special Issue: "The Early Data Release of eROSITA and Mikhail Pavlinsky ART-XC on the SRG Mission"

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

GOODS-ALMA 2.0: Source catalog, number counts, and prevailing compact sizes in 1.1 mm galaxies

Authors: C. Gómez-Guijarro, D. Elbaz, M. Xiao, M. Béthermin, M. Franco, B. Magnelli, E. Daddi, M. Dickinson, R. Demarco, H. Inami, W. Rujopakarn, G. E. Magdis, X. Shu, R. Chary, L. Zhou, D. M. Alexander, F. Bournaud, L. Ciesla, H. C. Ferguson, S. L. Finkelstein, M. Giavalisco, D. Iono, S. Juneau, J. S. Kartaltepe, G. Lagache , et al. (14 additional authors not shown)

Abstract: Submillimeter/millimeter observations of dusty star-forming galaxies with the Atacama Large Millimeter/submillimeter Array (ALMA) have shown that dust continuum emission generally occurs in compact regions smaller than the stellar distribution. However, it remains to be understood how systematic these findings are. Studies often lack homogeneity in the sample selection, target discontinuous areas… ▽ More Submillimeter/millimeter observations of dusty star-forming galaxies with the Atacama Large Millimeter/submillimeter Array (ALMA) have shown that dust continuum emission generally occurs in compact regions smaller than the stellar distribution. However, it remains to be understood how systematic these findings are. Studies often lack homogeneity in the sample selection, target discontinuous areas with inhomogeneous sensitivities, and suffer from modest $uv$ coverage coming from single array configurations. GOODS-ALMA is a 1.1mm galaxy survey over a continuous area of 72.42arcmin$^2$ at a homogeneous sensitivity. In this version 2.0, we present a new low resolution dataset and its combination with the previous high resolution dataset from the survey, improving the $uv$ coverage and sensitivity reaching an average of $σ= 68.4μ$Jy beam$^{-1}$. A total of 88 galaxies are detected in a blind search (compared to 35 in the high resolution dataset alone), 50% at $S/N_{peak} \geq 5$ and 50% at $3.5 \leq S/N_{peak} \leq 5$ aided by priors. Among them, 13 out of the 88 are optically dark or faint sources ($H$- or $K$-band dropouts). The sample dust continuum sizes at 1.1mm are generally compact, with a median effective radius of $R_{e} = 0"10 \pm 0"05$ (a physical size of $R_{e} = 0.73 \pm 0.29$kpc at the redshift of each source). Dust continuum sizes evolve with redshift and stellar mass resembling the trends of the stellar sizes measured at optical wavelengths, albeit a lower normalization compared to those of late-type galaxies. We conclude that for sources with flux densities $S_{1.1mm} > 1$mJy, compact dust continuum emission at 1.1mm prevails, and sizes as extended as typical star-forming stellar disks are rare. The $S_{1.1mm} < 1$mJy sources appear slightly more extended at 1.1mm, although they are still generally compact below the sizes of typical star-forming stellar disks. △ Less

Submitted 21 December, 2021; v1 submitted 24 June, 2021; originally announced June 2021.

Comments: Accepted by A&A. 30 pages, 18 figures

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

Optical Spectroscopy of Dual Quasar Candidates from the Subaru HSC-SSP program

Authors: Shenli Tang, John D. Silverman, Xuheng Ding, Junyao Li, Khee-Gan Lee, Michael A. Strauss, Andy Goulding, Malte Schramm, Lalitwadee Kawinwanichakij, J. Xavier Prochaska, Joseph F. Hennawi, Masatoshi Imanishi, Kazushi Iwasawa, Yoshiki Toba, Issha Kayo, Masamune Oguri, Yoshiki Matsuoka, Kohei Ichikawa, Tilman Hartwig, Nobunari Kashikawa, Toshihiro Kawaguchi, Kotaro Kohno, Yuichi Matsuda, Tohru Nagao, Yoshiaki Ono , et al. (8 additional authors not shown)

Abstract: We report on a spectroscopic program to search for dual quasars using Subaru Hyper Suprime-Cam (HSC) images of SDSS quasars which represent an important stage during galaxy mergers. Using Subaru/FOCAS and Gemini-N/GMOS, we identify three new physically associated quasar pairs having projected separations less than 20 kpc, out of 26 observed candidates. These include the discovery of the highest re… ▽ More We report on a spectroscopic program to search for dual quasars using Subaru Hyper Suprime-Cam (HSC) images of SDSS quasars which represent an important stage during galaxy mergers. Using Subaru/FOCAS and Gemini-N/GMOS, we identify three new physically associated quasar pairs having projected separations less than 20 kpc, out of 26 observed candidates. These include the discovery of the highest redshift ($z=3.1$) quasar pair with a separation $<$ 10 kpc. Based on the sample acquired to date, the success rate of identifying physically associated dual quasars is $19\%$ when excluding stars based on their HSC colors. Using the full sample of six spectroscopically confirmed dual quasars, we find that the black holes in these systems have black hole masses ($M_{BH} \sim 10^{8-9}M_{\odot}$) similar to single SDSS quasars as well as their bolometric luminosities and Eddington ratios. We measure the stellar mass of their host galaxies based on 2D image decomposition of the five-band ($grizy$) optical emission and assess the mass relation between supermassive black holes (SMBHs) and their hosts. Dual SMBHs appear to have elevated masses relative to their host galaxies. Thus mergers may not necessarily align such systems onto the local mass relation, as suggested by the Horizon-AGN simulation. This study suggests that dual luminous quasars are triggered prior to the final coalescence of the two SMBHs, resulting in early mass growth of the black holes relative to their host galaxies. △ Less

Submitted 21 May, 2021; originally announced May 2021.

Comments: 36 pages, 21 figures, 4 tables

The Sizes of Quasar Host Galaxies with the Hyper Suprime-Cam Subaru Strategic Program

Authors: Junyao Li, John D. Silverman, Xuheng Ding, Michael A. Strauss, Andy Goulding, Simon Birrer, Hassen M. Yesuf, Yongquan Xue, Lalitwadee Kawinwanichakij, Yoshiki Matsuoka, Yoshiki Toba, Tohru Nagao, Malte Schramm, Kohei Inayoshi

Abstract: The relationship between quasars and their host galaxies provides clues on how supermassive black holes (SMBHs) and massive galaxies are jointly assembled. To elucidate this connection, we measure the structural and photometric properties of the host galaxies of ~5000 SDSS quasars at 0.2<z<1 using five-band (grizy) optical imaging from the Hyper Suprime-Cam Subaru Strategic Program. An automated a… ▽ More The relationship between quasars and their host galaxies provides clues on how supermassive black holes (SMBHs) and massive galaxies are jointly assembled. To elucidate this connection, we measure the structural and photometric properties of the host galaxies of ~5000 SDSS quasars at 0.2<z<1 using five-band (grizy) optical imaging from the Hyper Suprime-Cam Subaru Strategic Program. An automated analysis tool is used to forward-model the blended emission of the quasar as characterized by the point spread function and the underlying host galaxy as a two-dimensional Sersic profile. In agreement with previous studies, quasars are preferentially hosted by massive star-forming galaxies with disk-like light profiles. Furthermore, we find that the size distribution of quasar hosts is broad at a given stellar mass and the average values exhibit a size-stellar mass relation as seen with inactive galaxies. In contrast, the sizes of quasar hosts are more compact than inactive star-forming galaxies on average, but not as compact as quiescent galaxies of similar stellar masses. This is true irrespective of quasar properties including bolometric luminosity, Eddington ratio, and black hole mass. These results are consistent with a scenario in which galaxies are concurrently fueling a SMBH and building their stellar bulge from a centrally-concentrated gas reservoir. Alternatively, quasar hosts may be experiencing a compaction process in which stars from the disk and inflowing gas are responsible for growing the bulge. In addition, we confirm that the host galaxies of type-1 quasars have a bias of being closer towards face-on, suggesting that galactic-scale dust can contribute to obscuring the broad-line region. △ Less

Submitted 13 May, 2021; originally announced May 2021.

Comments: 34 pages, 22 figures, 3 tables, resubmitted to ApJ after addressing referee's comments; catalog will be available after the paper is published

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

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

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

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

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

The ALPINE-ALMA [CII] Survey: Kinematic Diversity & Rotation in Massive Star Forming Galaxies at z~4.4-5.9

Authors: G. C. Jones, D. Vergani, M. Romano, M. Ginolfi, Y. Fudamoto, M. Bethermin, S. Fujimoto, B. C. Lemaux, L. Morselli, P. Capak, P. Cassata, A. Faisst, O. Le Fevre, D. Schaerer, J. D. Silverman, Lin Yan, M. Boquien, A. Cimatti, M. Dessauges-Zavadsky, E. Ibar, R. Maiolino, F. Rizzo, M. Talia, G. Zamorani

Abstract: While the kinematics of galaxies up to z~3 have been characterized in detail, only a handful of galaxies at high redshift (z>4) have been examined in such a way. The Atacama Large Millimeter/submillimeter Array (ALMA) Large Program to INvestigate [CII] at Early times (ALPINE) survey observed a statistically significant sample of 118 star-forming main sequence galaxies at z=4.4-5.9 in [CII]158um em… ▽ More While the kinematics of galaxies up to z~3 have been characterized in detail, only a handful of galaxies at high redshift (z>4) have been examined in such a way. The Atacama Large Millimeter/submillimeter Array (ALMA) Large Program to INvestigate [CII] at Early times (ALPINE) survey observed a statistically significant sample of 118 star-forming main sequence galaxies at z=4.4-5.9 in [CII]158um emission, increasing the number of such observations by nearly 10x. A preliminary qualitative classification of these sources revealed a diversity of kinematic types (i.e., rotators, mergers, and dispersion-dominated systems). In this work, we supplement the initial classification by applying quantitative analyses to the ALPINE data: a tilted ring model (TRM) fitting code (3DBarolo), a morphological classification (Gini-M20), and a set of disk identification criteria. Of the 75 [CII]-detected ALPINE galaxies, 29 are detected at sufficient significance and spatial resolution to allow for TRM fitting and the derivation of morphological and kinematic parameters. These 29 sources constitute a high-mass subset of the ALPINE sample (M_*>10^9.5Msol). We robustly classify 14 of these sources (six rotators, five mergers, and three dispersion-dominated systems); the remaining sources showing complex behaviour. By exploring the G-M20 of z>4 rest-frame FIR and [CII] data for the first time, we find that our 1"~6kpc resolution data alone are insufficient to separate galaxy types. We compare the rotation curves and dynamical mass profiles of the six ALPINE rotators to the two previously detected z~4-6 unlensed main sequence rotators, finding high rotational velocities (~50-250km/s) and a diversity of rotation curve shapes. △ Less

Submitted 2 August, 2021; v1 submitted 7 April, 2021; originally announced April 2021.

Comments: 25 pages, 15 figures. Accepted for publication in MNRAS