Showing 1–50 of 136 results for author: Goulding, A

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.

Quasars can Signpost Supermassive Black Hole Binaries

Authors: J. Andrew Casey-Clyde, Chiara M. F. Mingarelli, Jenny E. Greene, Andy D. Goulding, Siyuan Chen, Jonathan R. Trump

Abstract: Supermassive black holes (SMBHs) are found in the centers of massive galaxies, and galaxy mergers should eventually lead to SMBH mergers. Quasar activity has long been associated with galaxy mergers, so here we investigate if supermassive black hole binaries (SMBHBs) are preferentially found in quasars. Our multimessenger investigation folds together a gravitational wave background signal from NAN… ▽ More Supermassive black holes (SMBHs) are found in the centers of massive galaxies, and galaxy mergers should eventually lead to SMBH mergers. Quasar activity has long been associated with galaxy mergers, so here we investigate if supermassive black hole binaries (SMBHBs) are preferentially found in quasars. Our multimessenger investigation folds together a gravitational wave background signal from NANOGrav, a sample of periodic AGN candidates from the Catalina Real-Time Transient Survey, and a quasar mass function, to estimate an upper limit on the fraction of quasars which could host a SMBHB. We find at 95\% confidence that quasars are at most seven times as likely to host a SMBHB as a random galaxy. Quasars should therefore be prioritized as targets for SMBHB searches in pulsar timing arrays. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: 9 pages, 5 figures, 1 appendix

Medium Bands, Mega Science: a JWST/NIRCam Medium-Band Imaging Survey of Abell 2744

Authors: Katherine A. Suess, John R. Weaver, Sedona H. Price, Richard Pan, Bingjie Wang, Rachel Bezanson, Gabriel Brammer, Sam E. Cutler, Ivo Labbe, Joel Leja, Christina C. Williams, Katherine E. Whitaker, Pratika Dayal, Anna de Graaff, Robert Feldmann, Marijn Franx, Yoshinobu Fudamoto, Seiji Fujimoto, Lukas J. Furtak, Andy D. Goulding, Jenny E. Greene, Gourav Khullar, Vasily Kokorev, Mariska Kriek, Brian Lorenz , et al. (17 additional authors not shown)

Abstract: In this paper, we describe the "Medium Bands, Mega Science" JWST Cycle 2 survey (JWST-GO-4111) and demonstrate the power of these data to reveal both the spatially-integrated and spatially-resolved properties of galaxies from the local universe to the era of cosmic dawn. Executed in November 2023, MegaScience obtained ~30 arcmin^2 of deep multiband NIRCam imaging centered on the z~0.3 Abell 2744 c… ▽ More In this paper, we describe the "Medium Bands, Mega Science" JWST Cycle 2 survey (JWST-GO-4111) and demonstrate the power of these data to reveal both the spatially-integrated and spatially-resolved properties of galaxies from the local universe to the era of cosmic dawn. Executed in November 2023, MegaScience obtained ~30 arcmin^2 of deep multiband NIRCam imaging centered on the z~0.3 Abell 2744 cluster, including eleven medium-band filters and the two shortest-wavelength broad-band filters, F070W and F090W. Together, MegaScience and the UNCOVER Cycle 1 treasury program provide a complete set of deep (~28-30 mag) images in all NIRCam medium- and broad-band filters. This unique dataset allows us to precisely constrain photometric redshifts, map stellar populations and dust attenuation for large samples of distant galaxies, and examine the connection between galaxy structures and formation histories. MegaScience also includes ~17 arcmin^2 of NIRISS parallel imaging in two broad-band and four medium-band filters from 0.9-4.8um, expanding the footprint where robust spectral energy distribution (SED) fitting is possible. We provide example SEDs and multi-band cutouts at a variety of redshifts, and use a catalog of JWST spectroscopic redshifts to show that MegaScience improves both the scatter and catastrophic outlier rate of photometric redshifts by factors of 2-3. Additionally, we demonstrate the spatially-resolved science enabled by MegaScience by presenting maps of the [OIII] line emission and continuum emission in three spectroscopically-confirmed z>6 galaxies. We show that line emission in reionization-era galaxies can be clumpy, extended, and spatially offset from continuum emission, implying that galaxy assembly histories are complex even at these early epochs. We publicly release fully reduced mosaics and photometric catalogs for both the NIRCam primary and NIRISS parallel fields. △ Less

Submitted 19 April, 2024; originally announced April 2024.

Comments: 21 pages, 10 figures. Fully reduced imaging, photometric catalogs, and photometric redshift fits publicly available at https://jwst-uncover.github.io/megascience/

Estimating Galaxy Parameters with Self-Organizing Maps and the Effect of Missing Data

Authors: Valentina La Torre, Anna Sajina, Andy D. Goulding, Danilo Marchesini, Rachel Bezanson, Alan N. Pearl, Laerte Sodré Jr

Abstract: The current and upcoming large data volume galaxy surveys require the use of machine learning techniques to maximize their scientific return. This study explores the use of Self-Organizing Maps (SOMs) to estimate galaxy parameters with a focus on handling cases of missing data and providing realistic probability distribution functions for the parameters. We train a SOM with a simulated mass-limite… ▽ More The current and upcoming large data volume galaxy surveys require the use of machine learning techniques to maximize their scientific return. This study explores the use of Self-Organizing Maps (SOMs) to estimate galaxy parameters with a focus on handling cases of missing data and providing realistic probability distribution functions for the parameters. We train a SOM with a simulated mass-limited lightcone assuming a ugrizYJHKs+IRAC dataset, mimicking the Hyper Suprime-Cam (HSC) Deep joint dataset. For parameter estimation, we derive SOM likelihood surfaces considering photometric errors to derive total (statistical and systematic) uncertainties. We explore the effects of missing data including which bands are particular critical to the accuracy of the derived parameters. We demonstrate that the parameter recovery is significantly better when the missing bands are "filled-in" rather than if they are completely omitted. We propose a practical method for such recovery of missing data. △ Less

Submitted 27 March, 2024; originally announced March 2024.

Comments: 22 pages, 12 figures, 1 table, accepted for publication in The Astronomical Journal

Galaxy merger challenge: A comparison study between machine learning-based detection methods

Authors: B. Margalef-Bentabol, L. Wang, A. La Marca, C. Blanco-Prieto, D. Chudy, H. Domínguez-Sánchez, A. D. Goulding, A. Guzmán-Ortega, M. Huertas-Company, G. Martin, W. J. Pearson, V. Rodriguez-Gomez, M. Walmsley, R. W. Bickley, C. Bottrell, C. Conselice, D. O'Ryan

Abstract: Various galaxy merger detection methods have been applied to diverse datasets. However, it is difficult to understand how they compare. We aim to benchmark the relative performance of machine learning (ML) merger detection methods. We explore six leading ML methods using three main datasets. The first one (the training data) consists of mock observations from the IllustrisTNG simulations and allow… ▽ More Various galaxy merger detection methods have been applied to diverse datasets. However, it is difficult to understand how they compare. We aim to benchmark the relative performance of machine learning (ML) merger detection methods. We explore six leading ML methods using three main datasets. The first one (the training data) consists of mock observations from the IllustrisTNG simulations and allows us to quantify the performance metrics of the detection methods. The second one consists of mock observations from the Horizon-AGN simulations, introduced to evaluate the performance of classifiers trained on different, but comparable data. The third one consists of real observations from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) survey. For the binary classification task (mergers vs. non-mergers), all methods perform reasonably well in the domain of the training data. At $0.1<z<0.3$, precision and recall range between $\sim$70\% and 80\%, both of which decrease with increasing $z$ as expected (by $\sim$5\% for precision and $\sim$10\% for recall at $0.76<z<1.0$). When transferred to a different domain, the precision of all classifiers is only slightly reduced, but the recall is significantly worse (by $\sim$20-40\% depending on the method). Zoobot offers the best overall performance in terms of precision and F1 score. When applied to real HSC observations, all methods agree well with visual labels of clear mergers but can differ by more than an order of magnitude in predicting the overall fraction of major mergers. For the multi-class classification task to distinguish pre-, post- and non-mergers, none of the methods offer a good performance, which could be partly due to limitations in resolution and depth of the data. With the advent of better quality data (e.g. JWST and Euclid), it is important to improve our ability to detect mergers and distinguish between merger stages. △ Less

Submitted 15 April, 2024; v1 submitted 22 March, 2024; originally announced March 2024.

Journal ref: A&A 687, A24 (2024)

RUBIES: JWST/NIRSpec Confirmation of an Infrared-luminous, Broad-line Little Red Dot with an Ionized Outflow

Authors: Bingjie Wang, Anna de Graaff, Rebecca L. Davies, Jenny E. Greene, Joel Leja, Andy D. Goulding, Christina C. Williams, Gabriel B. Brammer, Katherine A. Suess, Andrea Weibel, Rachel Bezanson, Leindert A. Boogaard, Nikko J. Cleri, Michaela Hirschmann, Harley Katz, Ivo Labbe, Michael V. Maseda, Jorryt Matthee, Ian McConachie, Rohan P. Naidu, Pascal A. Oesch, Hans-Walter Rix, David J. Setton, Katherine E. Whitaker

Abstract: The JWST discovery of ``little red dots'' (LRDs) is reshaping our picture of the early Universe, yet the physical mechanisms driving their compact size and UV-optical colors remain elusive. Here we report an unusually bright LRD ($z=3.1$) observed as part of the RUBIES program. This LRD exhibits broad emission lines (FWHM $\sim4000$km/s), a blue UV continuum, a clear Balmer break and a red continu… ▽ More The JWST discovery of ``little red dots'' (LRDs) is reshaping our picture of the early Universe, yet the physical mechanisms driving their compact size and UV-optical colors remain elusive. Here we report an unusually bright LRD ($z=3.1$) observed as part of the RUBIES program. This LRD exhibits broad emission lines (FWHM $\sim4000$km/s), a blue UV continuum, a clear Balmer break and a red continuum sampled out to rest 4 $μ$m with MIRI. We develop a new joint galaxy and AGN model within the Prospector Bayesian inference framework and perform spectrophotometric modeling using NIRCam, MIRI, and NIRSpec/Prism observations. Our fiducial model reveals a $M_*\sim 10^9M_\odot$ galaxy alongside a dust-reddened AGN driving the optical emission. Explaining the rest-frame optical color as a reddened AGN requires $A_{\rm v}\gtrsim4$, suggesting that a great majority of the accretion disk energy is re-radiated as dust emission. Yet despite clear AGN signatures, we find a surprising lack of hot torus emission, which implies that either the dust emission in this object must be cold, or the red continuum must instead be driven by a massive, evolved stellar population of the host galaxy -- seemingly inconsistent with the high EW broad lines (H$α$ EW $\sim800$Å). The widths and luminosities of Pa$β$, Pa$δ$, Pa$γ$, and H$α$ imply a modest black hole mass of $M_{\rm BH}\sim10^8M_\odot$. Additionally, we identify a narrow blue-shifted HeI absorption in G395M spectra, signaling an ionized outflow with kinetic energy up to $\sim1$\% the luminosity of the AGN. The low redshift of RUBIES-BLAGN-1 combined with the depth and richness of the JWST imaging and spectroscopic observations provide a unique opportunity to build a physical model for these so-far mysterious LRDs, which may prove to be a crucial phase in the early formation of massive galaxies and their supermassive black holes. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: 20 pages, 8 figures, 4 tables

UNCOVER NIRSpec/PRISM Spectroscopy Unveils Evidence of Early Core Formation in a Massive, Centrally Dusty Quiescent Galaxy at $z_{spec}=3.97$

Authors: David J. Setton, Gourav Khullar, Tim B. Miller, Rachel Bezanson, Jenny E. Greene, Katherine A. Suess, Katherine E. Whitaker, Jacqueline Antwi-Danso, Hakim Atek, Gabriel Brammer, Sam E. Cutler, Pratika Dayal, Robert Feldmann, Lukas J. Furtak, Seiji Fujimoto, Karl Glazebrook, Andy D. Goulding, Vasily Kokorev, Ivo Labbe, Joel Leja, Yilun Ma, Danilo Marchesini, Themiya Nanayakkara, Richard Pan, Sedona H. Price , et al. (6 additional authors not shown)

Abstract: We report the spectroscopic confirmation of a massive ($\log(M_\star/M_\odot)=10.34 \pm_{0.07}^{0.06}$), HST-dark ($m_\mathrm{F150W} - m_\mathrm{F444W} = 3.6$) quiescent galaxy at $z_{spec}=3.97$ in the UNCOVER survey. NIRSpec/PRISM spectroscopy and a non-detection in deep ALMA imaging surprisingly reveals that the galaxy is consistent with a low ($<$10 $M_\odot \ \mathrm{yr^{-1}}$) star formation… ▽ More We report the spectroscopic confirmation of a massive ($\log(M_\star/M_\odot)=10.34 \pm_{0.07}^{0.06}$), HST-dark ($m_\mathrm{F150W} - m_\mathrm{F444W} = 3.6$) quiescent galaxy at $z_{spec}=3.97$ in the UNCOVER survey. NIRSpec/PRISM spectroscopy and a non-detection in deep ALMA imaging surprisingly reveals that the galaxy is consistent with a low ($<$10 $M_\odot \ \mathrm{yr^{-1}}$) star formation rate despite evidence for moderate dust attenuation. The F444W image is well modeled with a two component \sersic fit that favors a compact, $r_e\sim200$ pc, $n\sim2.9$ component and a more extended, $r_e\sim1.6$ kpc, $n\sim1.7$ component. The galaxy exhibits strong color gradients: the inner regions are significantly redder than the outskirts. Spectral energy distribution models that reproduce both the red colors and low star formation rate in the center of UNCOVER 18407 require both significant ($A_v\sim1.4$ mag) dust attenuation and a stellar mass-weighted age of 900 Myr, implying 50\% of the stars in the core already formed by $z=7.5$. Using spatially resolved annular mass-to-light measurements enabled by the galaxy's moderate magnification ($μ=2.12\pm_{0.01}^{0.05}$) to reconstruct a radial mass profile from the best-fitting two-component \sersic model, we infer a total mass-weighted $r_\mathrm{eff} = 0.72 \pm_{0.11}^{0.15}$ kpc and log$(Σ_\mathrm{1 kpc} \ [\mathrm{M_\odot/kpc^2}]) = 9.61 \pm_{0.10}^{0.08}$. The early formation of a dense, low star formation rate, and dusty core embedded in a less attenuated stellar envelope suggests an evolutionary link between the earliest-forming massive galaxies and their elliptical descendants. Furthermore, the disparity between the global, integrated dust properties and the spatially resolved gradients highlights the importance of accounting for radially varying stellar populations when characterizing the early growth of galaxy structure. △ Less

Submitted 12 May, 2024; v1 submitted 8 February, 2024; originally announced February 2024.

Comments: 17 pages, 9 figures, 2 tables. Resubmitted to ApJ after response to referee and update to include new medium band imaging from the JWST MEGASCIENCE program. Comments welcome!

Spectroscopic Confirmation of Obscured AGN Populations from Unsupervised Machine Learning

Authors: Raphael E. Hviding, Kevin N. Hainline, Andy D. Goulding, Jenny E. Greene

Abstract: We present the result of a spectroscopic campaign targeting Active Galactic Nucleus (AGN) candidates selected using a novel unsupervised machine-learning (ML) algorithm trained on optical and mid-infrared (mid-IR) photometry. AGN candidates are chosen without incorporating prior AGN selection criteria and are fainter, redder, and more numerous, $\sim$340 AGN deg$^{-2}$, than comparable photometric… ▽ More We present the result of a spectroscopic campaign targeting Active Galactic Nucleus (AGN) candidates selected using a novel unsupervised machine-learning (ML) algorithm trained on optical and mid-infrared (mid-IR) photometry. AGN candidates are chosen without incorporating prior AGN selection criteria and are fainter, redder, and more numerous, $\sim$340 AGN deg$^{-2}$, than comparable photometric and spectroscopic samples. In this work we obtain 178 rest-optical spectra from two candidate ML-identified AGN classes with the Hectospec spectrograph on the MMT Observatory. We find that our first ML-identified group, is dominated by Type I AGNs (85%) with a $<3$% contamination rate from non-AGNs. Our second ML-identified group is comprised mostly of Type II AGNs (65%) with a moderate contamination rate of 15% primarily from star-forming galaxies. Our spectroscopic analyses suggest that the classes recover more obscured AGNs, confirming that ML techniques are effective at recovering large populations of AGNs at high levels of extinction. We demonstrate the efficacy of pairing existing WISE data with large-area and deep optical/near-infrared photometric surveys to select large populations of AGNs and recover obscured SMBH growth. This approach is well suited to upcoming photometric surveys, such as Euclid, Rubin, and Roman. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: 23 pages, 16 figures. Accepted for publication in The Astronomical Journal

The HMXB Luminosity Functions of Dwarf Galaxies

Authors: Robel Geda, Andy D. Goulding, Bret D. Lehmer, Jenny E. Greene, Anish Kulkarni

Abstract: Drawing from the Chandra archive and using a carefully selected set of nearby dwarf galaxies, we present a calibrated high-mass X-ray binary (HMXB) luminosity function in the low-mass galaxy regime and search for an already hinted at dependence on metallicity. Our study introduces a new sample of local dwarf galaxies (D < 12.5 Mpc and M* < 5 x 10^9 M_sun), expanding the specific star-formation rat… ▽ More Drawing from the Chandra archive and using a carefully selected set of nearby dwarf galaxies, we present a calibrated high-mass X-ray binary (HMXB) luminosity function in the low-mass galaxy regime and search for an already hinted at dependence on metallicity. Our study introduces a new sample of local dwarf galaxies (D < 12.5 Mpc and M* < 5 x 10^9 M_sun), expanding the specific star-formation rates (sSFR) and gas-phase metallicities probed in previous investigations. Our analysis of the observed X-ray luminosity function indicates a shallower power-law slope for the dwarf galaxy HMXB population. In our study, we focus on dwarf galaxies that are more representative in terms of sSFR compared to prior work. In this regime, the HMXB luminosity function exhibits significant stochastic sampling at high luminosities. This likely accounts for the pronounced scatter observed in the galaxy-integrated HMXB population's Lx/SFR versus metallicity for our galaxy sample. Our calibration is necessary to understand the AGN content of low mass galaxies identified in current and future X-ray survey fields and has implications for binary population synthesis models, as well as X-ray driven cosmic heating in the early universe. △ Less

Submitted 25 January, 2024; originally announced January 2024.

Comments: 19 pages, 7 figures

UNCOVERing the contribution of black holes to reionization in the JWST era

Authors: Pratika Dayal, Marta Volonteri, Jenny E. Greene, Vasily Kokorev, Andy D. Goulding, Christina C. Williams, Lukas J. Furtak, Adi Zitrin, Hakim Atek, Iryna Chemerynska, Robert Feldmann, Karl Glazebrook, Ivo Labbe, Themiya Nanayakkara, Pascal A. Oesch, John R. Weaver

Abstract: With its sensitivity in the rest-frame optical, the James Webb Space Telescope (JWST) has uncovered active galactic nuclei (AGN), comprising both intrinsically faint and heavily reddened sources, well into the first billion years of the Universe, at $z \sim 4-11$. In this work, we revisit the AGN contribution to reionization given the high number densities associated with these objects. We use the… ▽ More With its sensitivity in the rest-frame optical, the James Webb Space Telescope (JWST) has uncovered active galactic nuclei (AGN), comprising both intrinsically faint and heavily reddened sources, well into the first billion years of the Universe, at $z \sim 4-11$. In this work, we revisit the AGN contribution to reionization given the high number densities associated with these objects. We use the DELPHI semi-analytic model, base-lined against the latest high-redshift datasets from the JWST and the Atacama Large millimetre Array (ALMA) to model early star forming galaxies and AGN. We calculate the escape fractions of ionizing radiation from both star formation and AGN and include the impact of reionization feeback in suppressing the baryonic content of low-mass galaxies in ionized regions. This model is validated against the key observables for star forming galaxy, AGN and reionization. In our {\it fiducial} model, reionization reaches its mid-point at $z \sim 6.9$ and ends by $z \sim 5.9$. Low stellar mass ($M_*\leq 10^9M_\odot$) star forming galaxies are found to be the key drivers of the reionization process, providing about $77\%$ of the total photon budget. Despite their high numbers, high accretion rates and higher escape fractions compared to star forming galaxies at $z \sim 5$, AGN only provide about $23\%$ of the total reionization budget which is dominated by black holes in high stellar mass systems (with $M_* \geq 10^9M_\odot$). This is because AGN number densities become relevant only at $z \leq 7$ - as a result, AGN contribute as much as galaxies as late as $z \sim 6.2$, when reionization is already in its end stages. Finally, we find that even contrasting models of the AGN ionizing photon escape fraction (increasing or decreasing with stellar mass) do not qualitatively change our results. △ Less

Submitted 20 January, 2024; originally announced January 2024.

Comments: Submitted to A&A; comments welcome

Evidence for Intrinsic X-ray Weakness Among Red Quasars at Cosmic Noon

Authors: Yilun Ma, Andy Goulding, Jenny E. Greene, Nadia L. Zakamska, Dominika Wylezalek, Yan-Fei Jiang

Abstract: Quasar feedback is a key ingredient in shaping galaxy evolution. A rare population of extremely red quasars (ERQs) at $z=2-3$ are often associated with high-velocity [OIII]$\lambda5008$ outflows and may represent sites of strong feedback. In this paper, we present an X-ray study of 50 ERQs to investigate the link between the X-ray and outflow properties of these intriguing objects. Using hardness… ▽ More Quasar feedback is a key ingredient in shaping galaxy evolution. A rare population of extremely red quasars (ERQs) at $z=2-3$ are often associated with high-velocity [OIII]$\lambda5008$ outflows and may represent sites of strong feedback. In this paper, we present an X-ray study of 50 ERQs to investigate the link between the X-ray and outflow properties of these intriguing objects. Using hardness ratio analysis, we confirm that the ERQs are heavily obscured systems with gas column density reaching $N_\mathrm{H}=10^{23-24}\,\mathrm{cm^{-2}}$. We identify 17 X-ray-non-detected ERQs at high mid-infrared luminosities of $νL_\mathrm{ν,6μm}\gtrsim3\times10^{46}\,\mathrm{erg\,s^{-1}}$. By stacking the X-ray observations, we find that the non-detected ERQs are on average underluminous in X-rays by a factor of $\sim5$ for their bolometric luminosities. We consider such X-ray weakness to be due to both gas absorption and intrinsic factors. Moreover, we find that the X-ray-weak sources also display higher-velocity outflows. One option to explain this trend is that weaker X-rays facilitate more vigorous line-driven winds, which then accelerate the [OIII]-emitting gas to kpc-scales. Alternatively, super-Eddington accretion could also lead to intrinsic X-ray weakness and more powerful continuum-driven outflow. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: 22 pages, 10 figures, submitted to ApJ

JWST UNCOVER: The Overabundance of Ultraviolet-luminous Galaxies at $z>9$

Authors: Iryna Chemerynska, Hakim Atek, Lukas J. Furtak, Adi Zitrin, Jenny E. Greene, Pratika Dayal, Andrea Weibel, Vasily Kokorev, Andy D. Goulding, Christina C. Williams, Themiya Nanayakkara, Rachel Bezanson, Gabriel Brammer, Sam E. Cutler, Ivo Labbe, Joel Leja, Richard Pan, Sedona H. Price, Bingjie Wang, John R. Weaver, Katherine E. Whitaker

Abstract: Over the past year, JWST has uncovered galaxies at record-breaking distances up to $z \sim 13$. The JWST UNCOVER (ultra-deep NIRSpec and NIRcam observations before the epoch of reionization) program has obtained ultra-deep multiwavelength NIRCam imaging of the massive galaxy cluster Abell 2744 over $\sim 45$ arcmin$^{2}$ down to $\sim 29.5$ AB mag. Here, we present a robust ultraviolet (UV) lumino… ▽ More Over the past year, JWST has uncovered galaxies at record-breaking distances up to $z \sim 13$. The JWST UNCOVER (ultra-deep NIRSpec and NIRcam observations before the epoch of reionization) program has obtained ultra-deep multiwavelength NIRCam imaging of the massive galaxy cluster Abell 2744 over $\sim 45$ arcmin$^{2}$ down to $\sim 29.5$ AB mag. Here, we present a robust ultraviolet (UV) luminosity function derived through lensing clusters at $9<z<12$. Using comprehensive end-to-end simulations, we account for all lensing effects and systematic uncertainties in deriving both the amplification factors and the effective survey volume. Our results confirm the intriguing excess of UV-bright galaxies ($M_{UV} < -20$ AB mag) previously reported at $z>9$ in recent JWST studies. In particular, a double power-law (DPL) describes better the bright-end of the luminosity function compared to the classical Schechter form. The number density of these bright galaxies is 10-100 times larger than theoretical predictions and previous findings based on Hubble Space Telescope (HST) observations. Additionally, we measure a star formation rate density of $ρ_{\rm SFR} = 10^{-2.64}$ M$_{\odot}$ yr$^{-1}$ Mpc$^{-3}$ at these redshifts, which is 4 to 10 times higher than galaxy formation models that assume a constant star formation efficiency. Future wide-area surveys and accurate modeling of lensing-assisted observations will reliably constrain both the bright and the dim end of the UV luminosity function at $z>9$, which will provide key benchmarks for galaxy formation models. △ Less

Submitted 8 December, 2023; originally announced December 2023.

Comments: Submitted to MNRAS

Constructing Impactful Machine Learning Research for Astronomy: Best Practices for Researchers and Reviewers

Authors: D. Huppenkothen, M. Ntampaka, M. Ho, M. Fouesneau, B. Nord, J. E. G. Peek, M. Walmsley, J. F. Wu, C. Avestruz, T. Buck, M. Brescia, D. P. Finkbeiner, A. D. Goulding, T. Kacprzak, P. Melchior, M. Pasquato, N. Ramachandra, Y. -S. Ting, G. van de Ven, S. Villar, V. A. Villar, E. Zinger

Abstract: Machine learning has rapidly become a tool of choice for the astronomical community. It is being applied across a wide range of wavelengths and problems, from the classification of transients to neural network emulators of cosmological simulations, and is shifting paradigms about how we generate and report scientific results. At the same time, this class of method comes with its own set of best pr… ▽ More Machine learning has rapidly become a tool of choice for the astronomical community. It is being applied across a wide range of wavelengths and problems, from the classification of transients to neural network emulators of cosmological simulations, and is shifting paradigms about how we generate and report scientific results. At the same time, this class of method comes with its own set of best practices, challenges, and drawbacks, which, at present, are often reported on incompletely in the astrophysical literature. With this paper, we aim to provide a primer to the astronomical community, including authors, reviewers, and editors, on how to implement machine learning models and report their results in a way that ensures the accuracy of the results, reproducibility of the findings, and usefulness of the method. △ Less

Submitted 19 October, 2023; originally announced October 2023.

Comments: 14 pages, 3 figures; submitted to the Bulletin of the American Astronomical Society

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

UNCOVER spectroscopy confirms a surprising ubiquity of AGN in red galaxies at $z>5$

Authors: Jenny E. Greene, Ivo Labbe, Andy D. Goulding, Lukas J. Furtak, Iryna Chemerynska, Vasily Kokorev, Pratika Dayal, Christina C. Williams, Bingjie Wang, David J. Setton, Adam J. Burgasser, Rachel Bezanson, Hakim Atek, Gabriel Brammer, Sam E. Cutler, Robert Feldmann, Seiji Fujimoto, Karl Glazebrook, Anna de Graaff, Joel Leja, Danilo Marchesini, Michael V. Maseda, Jorryt Matthee, Tim B. Miller, Rohan P. Naidu , et al. (9 additional authors not shown)

Abstract: JWST is revealing a new population of dust-reddened broad-line active galactic nuclei (AGN) at redshifts $z\gtrsim5$. Here we present deep NIRSpec/Prism spectroscopy from the Cycle 1 Treasury program UNCOVER of 15 AGN candidates selected to be compact, with red continua in the rest-frame optical but with blue slopes in the UV. From NIRCam photometry alone, they could have been dominated by dusty s… ▽ More JWST is revealing a new population of dust-reddened broad-line active galactic nuclei (AGN) at redshifts $z\gtrsim5$. Here we present deep NIRSpec/Prism spectroscopy from the Cycle 1 Treasury program UNCOVER of 15 AGN candidates selected to be compact, with red continua in the rest-frame optical but with blue slopes in the UV. From NIRCam photometry alone, they could have been dominated by dusty star formation or AGN. Here we show that the majority of the compact red sources in UNCOVER are dust-reddened AGN: $60\%$ show definitive evidence for broad-line H$α$ with FWHM$\, >2000$ km/s, for $20\%$ current data are inconclusive, and $20\%$ are brown dwarf stars. We propose an updated photometric criterion to select red $z>5$ AGN that excludes brown dwarfs and is expected to yield $>80\%$ AGN. Remarkably, among all $z_{\rm phot}>5$ galaxies with F277W$-$F444W$>1$ in UNCOVER at least $33\%$ are AGN regardless of compactness, climbing to at least $80\%$ AGN for sources with F277W$-$F444W$>1.6$. The confirmed AGN have black hole masses of $10^7-10^9$ M$_{\odot}$. While their UV-luminosities ($-16>M_{\rm UV}>-20$ AB mag) are low compared to UV-selected AGN at these epochs, consistent with percent-level scattered AGN light or low levels of unobscured star formation, the inferred bolometric luminosities are typical of $10^7-10^9$ M$_{\odot}$ black holes radiating at $\sim 10-40\%$ of Eddington. The number densities are surprisingly high at $\sim10^{-5}$ Mpc$^{-3}$ mag$^{-1}$, 100 times more common than the faintest UV-selected quasars, while accounting for $\sim1\%$ of the UV-selected galaxies. While their UV-faintness suggest they may not contribute strongly to reionization, their ubiquity poses challenges to models of black hole growth. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: 23 pages, 9 figures, 5 tables, submitted to 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

UNCOVER: A NIRSpec Identification of a Broad Line AGN at z = 8.50

Authors: Vasily Kokorev, Seiji Fujimoto, Ivo Labbe, Jenny E. Greene, Rachel Bezanson, Pratika Dayal, Erica J. Nelson, Hakim Atek, Gabriel Brammer, Karina I. Caputi, Iryna Chemerynska, Sam E. Cutler, Robert Feldmann, Yoshinobu Fudamoto, Lukas J. Furtak, Andy D. Goulding, Anna de Graaff, Joel Leja, Danilo Marchesini, Tim B. Miller, Themiya Nanayakkara, Pascal Oesch, Richard Pan, Sedona H. Price, David J. Setton , et al. (7 additional authors not shown)

Abstract: Deep observations with JWST have revealed an emerging population of red point-like sources that could provide a link between the postulated supermassive black hole seeds and observed quasars. In this work we present a JWST/NIRSpec spectrum from the JWST Cycle 1 UNCOVER Treasury survey, of a massive accreting black hole at $z=8.50$, displaying a clear broad-line component as inferred from the H$β$… ▽ More Deep observations with JWST have revealed an emerging population of red point-like sources that could provide a link between the postulated supermassive black hole seeds and observed quasars. In this work we present a JWST/NIRSpec spectrum from the JWST Cycle 1 UNCOVER Treasury survey, of a massive accreting black hole at $z=8.50$, displaying a clear broad-line component as inferred from the H$β$ line with FWHM = $3439\pm413$ km s$^{-1}$, typical of the broad line region of an active galactic nucleus (AGN). The AGN nature of this object is further supported by high ionization, as inferred from emission lines, and a point-source morphology. We compute the black hole mass of log$_{10}(M_{\rm BH}/M_\odot)=8.17\pm0.42$, and a bolometric luminosity of $L_{\rm bol}\sim6.6\times10^{45}$ erg s$^{-1}$. These values imply that our object is accreting at $\sim 40\%$ of the Eddington limit. Detailed modeling of the spectral energy distribution in the optical and near-infrared, together with constraints from ALMA, indicate an upper limit on the stellar mass of log$_{10}(M_{\rm *}/M_\odot)<8.7$, which would lead to an unprecedented ratio of black hole to host mass of at least $\sim 30 \%$. This is orders of magnitude higher compared to the local QSOs, but is consistent with recent AGN studies at high redshift with JWST. This finding suggests that a non-negligible fraction of supermassive black holes either started out from massive seeds and/or grew at a super-Eddington rate at high redshift. Given the predicted number densities of high-$z$ faint AGN, future NIRSpec observations of larger samples will allow us to further investigate the galaxy-black hole co-evolution in the early Universe. △ Less

Submitted 15 October, 2023; v1 submitted 22 August, 2023; originally announced August 2023.

Comments: 14 pages, 6 figures, 2 tables. Accepted for a publication in ApJL

UNCOVER: A NIRSpec Census of Lensed Galaxies at z=8.50-13.08 Probing a High AGN Fraction and Ionized Bubbles in the Shadow

Authors: Seiji Fujimoto, Bingjie Wang, John Weaver, Vasily Kokorev, Hakim Atek, Rachel Bezanson, Ivo Labbe, Gabriel Brammer, Jenny E. Greene, Iryna Chemerynska, Pratika Dayal, Anna de Graaff, Lukas J. Furtak, Pascal A. Oesch, David J. Setton, Sedona H. Price, Tim B. Miller, Christina C. Williams, Katherine E. Whitaker, Adi Zitrin, Sam E. Cutler, Joel Leja, Richard Pan, Dan Coe, Pieter van Dokkum , et al. (11 additional authors not shown)

Abstract: We present JWST NIRSpec prism spectroscopy of gravitationally lensed galaxies at $z\gtrsim9$ found behind the massive galaxy cluster Abell 2744 in the UNCOVER Cycle 1 Treasury Program. We confirm the source redshift via emission lines and/or the Ly$α$ break feature for ten galaxies at z=8.50-13.08 down to $M_{\rm UV}=-17.3$. We achieve a high confirmation rate of 100\% for $z>9$ candidates reporte… ▽ More We present JWST NIRSpec prism spectroscopy of gravitationally lensed galaxies at $z\gtrsim9$ found behind the massive galaxy cluster Abell 2744 in the UNCOVER Cycle 1 Treasury Program. We confirm the source redshift via emission lines and/or the Ly$α$ break feature for ten galaxies at z=8.50-13.08 down to $M_{\rm UV}=-17.3$. We achieve a high confirmation rate of 100\% for $z>9$ candidates reported in Atek et al. (2023). Using six sources with multiple emission line detections, we find that the offset of the redshift estimates between the lines and the Ly$α$ break alone with prism can be as large as $\pm0.2$, raising caution in designing future follow-up spectroscopy for the break-only sources. With spec-$z$ confirmed sources in UNCOVER and the literature, we derive lower limits on the rest-frame ultraviolet (UV) luminosity function (LF) at $z\simeq9$-12 and find these lower limits to be consistent with recent photometric measurements. We identify at least two unambiguous and several possible active galactic nucleus (AGN) systems based on X-ray emission, broad line (BL) H$β$, high ionization line (e.g., NIV]1487, CIV1549) detections, and excess in UVLF. This requires the AGN LFs at $z\simeq$ 9-10 to be comparable or even higher than the X-ray AGN LF estimated at $z\sim6$ and indicates a plausible cause of the high abundance of $z>9$ galaxies claimed in recent photometric studies may be AGNs. One UV-luminous source is confirmed at the same redshift as a dusty BL AGN at $z=8.50$ with a physical separation of 380 kpc in the source plane. These two sources show blueward Ly$α$ line or continuum emission, suggesting that they reside in the same ionized bubble with a radius of $7.69\pm0.18$ pMpc. Our results imply that AGNs have a non-negligible contribution to cosmic reionization. △ Less

Submitted 25 August, 2023; v1 submitted 22 August, 2023; originally announced August 2023.

Comments: 27 pages, 11 figures, 4 tables, submitted to ApJ (See also arXiv:2308.11610)

A supermassive black hole in the early universe growing in the shadows

Authors: Lukas J. Furtak, Ivo Labbé, Adi Zitrin, Jenny E. Greene, Pratika Dayal, Iryna Chemerynska, Vasily Kokorev, Tim B. Miller, Andy D. Goulding, Rachel Bezanson, Gabriel B. Brammer, Sam E. Cutler, Joel Leja, Richard Pan, Sedona H. Price, Bingjie Wang, John R. Weaver, Katherine E. Whitaker, Hakim Atek, Ákos Bogdán, Stéphane Charlot, Emma Curtis-Lake, Pieter van Dokkum, Ryan Endsley, Yoshinobu Fudamoto , et al. (12 additional authors not shown)

Abstract: Early JWST observations have uncovered a new, substantial population of red sources that might represent a previously overlooked phase of actively growing supermassive black holes (Kocevski et al. 2023, Matthee et al. 2023, Labbe et al. 2023). One of the most intriguing examples is an extremely red, point-like object that was found to be triply-imaged by the strong lensing galaxy cluster Abell 274… ▽ More Early JWST observations have uncovered a new, substantial population of red sources that might represent a previously overlooked phase of actively growing supermassive black holes (Kocevski et al. 2023, Matthee et al. 2023, Labbe et al. 2023). One of the most intriguing examples is an extremely red, point-like object that was found to be triply-imaged by the strong lensing galaxy cluster Abell 2744 (Furtak et al. 2023), allowing an unprecedented detailed look into this enigmatic population. Here we present deep spectroscopic JWST/NIRSpec observations of this object, Abell2744-QSO1. The spectroscopy confirms that the three images are of the same object, and that it is a highly reddened ($A_V\sim3$) broad emission-line Active Galactic Nucleus (AGN) at a redshift of $z_{\mathrm{spec}}=7.0451\pm0.0005$. From the width of H$β$ ($\mathrm{FWHM}=2800\pm250\,\frac{\mathrm{km}}{\mathrm{s}}$) we derive a black hole mass of $M_{\mathrm{BH}}=3_{-1}^{+2}\times10^7\,\mathrm{M}_{\odot}$. We infer a very high ratio of black hole to galaxy mass of at least 3% and possibly as high as 100%, an order of magnitude or more than is seen in local galaxies. The lack of strong metal lines in the spectrum together with the high bolometric luminosity ($L_{\mathrm{bol}}=(1.1\pm0.3)\times10^{45}\,\frac{\mathrm{erg}}{\mathrm{s}}$) suggest that we are seeing the black hole in a phase of rapid growth, accreting at 30% of the Eddington limit. Based on early JWST imaging studies we estimate that such heavily reddened, low-mass black holes can be $\sim100$ times more common than UV-selected ones at this epoch. The rapid growth and high black hole to galaxy mass ratio of A2744-QSO1 suggests that it may represent the missing link between black hole seeds (Inayoshi et al. 2020; Greene et al. 2020; Volonteri 2021) and the first luminous quasars (Fan et al. 2023). △ Less

Submitted 10 August, 2023; originally announced August 2023.

Comments: Submitted. Comments welcome!

UNCOVER: Illuminating the Early Universe -- JWST/NIRSpec Confirmation of $z > 12$ Galaxies

Authors: Bingjie Wang, Seiji Fujimoto, Ivo Labbe, Lukas J. Furtak, Tim B. Miller, David J. Setton, Adi Zitrin, Hakim Atek, Rachel Bezanson, Gabriel Brammer, Joel Leja, Pascal A. Oesch, Sedona H. Price, Iryna Chemerynska, Sam E. Cutler, Pratika Dayal, Pieter van Dokkum, Andy D. Goulding, Jenny E. Greene, Y. Fudamoto, Gourav Khullar, Vasily Kokorev, Danilo Marchesini, Richard Pan, John R. Weaver , et al. (2 additional authors not shown)

Abstract: Observations of high-redshift galaxies provide a critical direct test to the theories of early galaxy formation, yet to date, only three have been spectroscopically confirmed at $z>12$. Due to strong gravitational lensing over a wide area, the galaxy cluster field A2744 is ideal for searching for the earliest galaxies. Here we present JWST/NIRSpec observations of two galaxies: a robust detection a… ▽ More Observations of high-redshift galaxies provide a critical direct test to the theories of early galaxy formation, yet to date, only three have been spectroscopically confirmed at $z>12$. Due to strong gravitational lensing over a wide area, the galaxy cluster field A2744 is ideal for searching for the earliest galaxies. Here we present JWST/NIRSpec observations of two galaxies: a robust detection at $z_{\rm spec} = 12.393^{+0.004}_{-0.001}$, and a plausible candidate at $z_{\rm spec} = 13.079^{+0.013}_{-0.001}$. The galaxies are discovered in JWST/NIRCam imaging and their distances are inferred with JWST/NIRSpec spectroscopy, all from the JWST Cycle 1 UNCOVER Treasury survey. Detailed stellar population modeling using JWST NIRCam and NIRSpec data corroborates the primeval characteristics of these galaxies: low mass ($\sim 10^8~{\rm M_\odot}$), young, rapidly-assembling, metal-poor, and star-forming. Interestingly, both galaxies are spatially resolved, having lensing-corrected rest-UV effective radii on the order of 300-400 pc, which are notably larger than other spectroscopically confirmed systems at similar redshifts. The observed dynamic range of $z \gtrsim 10$ sizes spans over 1 order of magnitude, implying a significant scatter in the size-mass relation at early times. Deep into the epoch of reionization, these discoveries elucidate the emergence of the first galaxies. △ Less

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

Comments: ApJL in press; 16 pages, 6 figures, 2 tables

Journal ref: The Astrophysical Journal Letters, 957, L34 (2023)

UNCOVER: The growth of the first massive black holes from JWST/NIRSpec -- spectroscopic redshift confirmation of an X-ray luminous AGN at z=10.1

Authors: Andy D. Goulding, Jenny E. Greene, David J. Setton, Ivo Labbe, Rachel Bezanson, Tim B. Miller, Hakim Atek, Akos Bogdan, Gabriel Brammer, Iryna Chemerynska, Sam E. Cutler, Pratika Dayal, Yoshinobu Fudamoto, Seiji Fujimoto, Lukas J. Furtak, Vasily Kokorev, Gourav Khullar, Joel Leja, Danilo Marchesini, Priyamvada Natarajan, Erica Nelson, Pascal A. Oesch, Richard Pan, Casey Papovich, Sedona H. Price , et al. (5 additional authors not shown)

Abstract: The James Webb Space Telescope is now detecting early black holes (BHs) as they transition from "seeds" to supermassive BHs. Recently Bogdan et al. (2023) reported the detection of an X-ray luminous supermassive BH, UHZ-1, with a photometric redshift at $z > 10$. Such an extreme source at this very high redshift provides new insights on seeding and growth models for BHs given the short time availa… ▽ More The James Webb Space Telescope is now detecting early black holes (BHs) as they transition from "seeds" to supermassive BHs. Recently Bogdan et al. (2023) reported the detection of an X-ray luminous supermassive BH, UHZ-1, with a photometric redshift at $z > 10$. Such an extreme source at this very high redshift provides new insights on seeding and growth models for BHs given the short time available for formation and growth. Harnessing the exquisite sensitivity of JWST/NIRSpec, here we report the spectroscopic confirmation of UHZ-1 at $z = 10.073 \pm 0.002$. We find that the NIRSpec/Prism spectrum is typical of recently discovered z~10 galaxies, characterized primarily by star-formation features. We see no clear evidence of the powerful X-ray source in the rest-frame UV/optical spectrum, which may suggest heavy obscuration of the central BH, in line with the Compton-thick column density measured in the X-rays. We perform a stellar population fit simultaneously to the new NIRSpec spectroscopy and previously available photometry. The fit yields a stellar mass estimate for the host galaxy that is significantly better constrained than prior photometric estimates ($M_*\sim 1.4^{+0.3}_{-0.4} \times 10^8 M_\odot$). Given the predicted BH mass ($M_{\rm BH}\sim10^7-10^8 M_\odot$), the resulting ratio of $M_{\rm BH}/M_*$ remains two to three orders of magnitude higher than local values, thus lending support to the heavy seeding channel for the formation of supermassive BHs within the first billion years of cosmic evolution. △ Less

Submitted 19 September, 2023; v1 submitted 4 August, 2023; originally announced August 2023.

Comments: 10 pages, 4 figures, ApJL in-press. Updated to accepted version incl. X-ray image

First Detection of an Over-Massive Black Hole Galaxy UHZ1: Evidence for Heavy Black Hole Seed Formation from Direct Collapse

Authors: Priyamvada Natarajan, Fabio Pacucci, Angelo Ricarte, Akos Bogdan, Andy D. Goulding, Nico Cappelluti

Abstract: The recent Chandra-JWST discovery of a quasar in the z = 10.1 galaxy UHZ1 reveals that accreting supermassive black holes (SMBHs) were already in place 470 million years after the Big Bang. The Chandra X-ray source detected in UHZ1 is a Compton-thick quasar with a bolometric luminosity of $L_{\rm bol}\sim5\times10^{45}\ \rm{erg\ s^{-1}},$ which corresponds to an estimated BH mass of… ▽ More The recent Chandra-JWST discovery of a quasar in the z = 10.1 galaxy UHZ1 reveals that accreting supermassive black holes (SMBHs) were already in place 470 million years after the Big Bang. The Chandra X-ray source detected in UHZ1 is a Compton-thick quasar with a bolometric luminosity of $L_{\rm bol}\sim5\times10^{45}\ \rm{erg\ s^{-1}},$ which corresponds to an estimated BH mass of $\sim4\times10^7 \ \rm{M_{\odot}}$ assuming accretion at the Eddington rate. JWST NIRCAM and NIRSpec data yield a stellar mass estimate for UHZ1 comparable to its BH mass. These characteristics are in excellent agreement with prior theoretical predictions for a unique class of transient, high-redshift objects, Over-massive Black Hole Galaxies [OBGs] by Natarajan et al. that harbor a heavy initial black hole seed that likely formed from the direct collapse of the gas. Based on the excellent agreement between the observed multi-wavelength properties of UHZ1 with theoretical model template predictions, suggests that UHZ1 is the first detected OBG candidate. Our assertion rests on multiple lines of concordant evidence between model predictions and the following observed properties of UHZ1: its X-ray detection and the estimated ratio of the X-ray flux to the IR flux that is consistent with theoretical expectations for a heavy initial BH seed; its high measured redshift of z = 10.1, as predicted for the transient OBG stage (9 < z< 12); the amplitude and shape of the detected JWST Spectral Energy Distribution (SED) between 1 - 5 microns, which is in very good agreement with simulated template SEDs for OBGs; and the extended JWST morphology of UHZ1 that is suggestive of a recent merger, also expected for the formation of transient OBGs. As the first OBG candidate, UHZ1 provides compelling evidence for the formation of heavy initial seeds from direct collapse in the early Universe. △ Less

Submitted 25 November, 2023; v1 submitted 4 August, 2023; originally announced August 2023.

Comments: 9 pages, 4 figures, accepted ApJ Letters

Modeling the SED of the AGN inside NGC 4395

Authors: Hector Afonso G. Cruz, Andy D. Goulding, Jenny E. Greene

Abstract: We study the broad-band spectral energy distribution (SED) of the prototypical low-mass active galactic nucleus (AGN) in NGC 4395. We jointly model the optical through mid-infrared SED with a combination of galaxy and AGN light, and find that on arcsecond scales, the AGN dominates at most wavelengths. However, there is still some ambiguity about emission from the galaxy, owing partially to the str… ▽ More We study the broad-band spectral energy distribution (SED) of the prototypical low-mass active galactic nucleus (AGN) in NGC 4395. We jointly model the optical through mid-infrared SED with a combination of galaxy and AGN light, and find that on arcsecond scales, the AGN dominates at most wavelengths. However, there is still some ambiguity about emission from the galaxy, owing partially to the strong short-term variability of the black hole. We investigate the use of smooth and clumpy-torus models in order to disentangle the nuclear infrared emission, as well as exploring the use of poloidal wind emission to account for the blue spectral slope observed in the near-IR. Even when simultaneously fitting the full optical-IR spectral range, we find that degeneracies still remain in the best-fit models. We conclude that high spatial resolution and wider wavelength coverage with the James Webb Space Telescope is needed to understand the mid-infrared emission in this complex highly-variable object, which is the best nearby example to provide a blueprint to finding other low-mass AGN via their mid-infrared emission in the future. △ Less

Submitted 27 July, 2023; originally announced July 2023.

Comments: 17 pages, 8 figures

First results from the JWST Early Release Science Program Q3D: Powerful quasar-driven galactic scale outflow at $z=3$

Authors: Andrey Vayner, Nadia L. Zakamska, Yuzo Ishikawa, Swetha Sankar, Dominika Wylezalek, David S. N. Rupke, Sylvain Veilleux, Caroline Bertemes, Jorge K. Barrera-Ballesteros, Hsiao-Wen Chen, Nadiia Diachenko, Andy D. Goulding, Jenny E. Greene, Kevin N. Hainline, Fred Hamann, Timothy Heckman, Sean D. Johnson, Hui Xian Grace Lim, Weizhe Liu, Dieter Lutz, Nora Lutzgendorf, Vincenzo Mainieri, Ryan McCrory, Grey Murphree, Nicole P. H. Nesvadba , et al. (3 additional authors not shown)

Abstract: Quasar-driven galactic outflows are a major driver of the evolution of massive galaxies. We report observations of a powerful galactic-scale outflow in a $z=3$ extremely red, intrinsically luminous ($L_{\rm bol}\simeq 5\times 10^{47}$erg s$^{-1}$) quasar SDSSJ1652+1728 with the Near Infrared Spectrograph (NIRSpec) on board JWST. We analyze the kinematics of rest-frame optical emission lines and id… ▽ More Quasar-driven galactic outflows are a major driver of the evolution of massive galaxies. We report observations of a powerful galactic-scale outflow in a $z=3$ extremely red, intrinsically luminous ($L_{\rm bol}\simeq 5\times 10^{47}$erg s$^{-1}$) quasar SDSSJ1652+1728 with the Near Infrared Spectrograph (NIRSpec) on board JWST. We analyze the kinematics of rest-frame optical emission lines and identify the quasar-driven outflow extending out to $\sim 10$ kpc from the quasar with a velocity offset of ($v_{r}=\pm 500$ km s$^{-1}$) and high velocity dispersion (FWHM$=700-2400$ km s$^{-1}$). Due to JWST's unprecedented surface brightness sensitivity in the near-infrared -- we unambiguously show that the powerful high velocity outflow in an extremely red quasar (ERQ) encompasses a large swath of the host galaxy's interstellar medium (ISM). Using the kinematics and dynamics of optical emission lines, we estimate the mass outflow rate -- in the warm ionized phase alone -- to be at least $2300\pm1400$ $M_{\odot}$ yr$^{-1}$. We measure a momentum flux ratio between the outflow and the quasar accretion disk of $\sim$1 on kpc scale, indicating that the outflow was likely driven in a relatively high ($>10^{23}$cm$^{-2}$) column density environment through radiation pressure on dust grains. We find a coupling efficiency between the bolometric luminosity of the quasar and the outflow of 0.1$\%$, matching the theoretical prediction of the minimum coupling efficiency necessary for negative quasar feedback. The outflow has sufficient energetics to drive the observed turbulence seen in shocked regions of the quasar host galaxy, likely directly responsible for prolonging the time it takes for gas to cool efficiently. △ Less

Submitted 25 July, 2023; originally announced July 2023.

Comments: 14 pages, 7 figures, 1 table, submitted for publication in ApJ

Outlier Detection in the DESI Bright Galaxy Survey

Authors: Yan Liang, Peter Melchior, ChangHoon Hahn, Jeff Shen, Andy Goulding, Charlotte Ward

Abstract: We present an unsupervised search for outliers in the Bright Galaxy Survey (BGS) dataset from the DESI Early Data Release. This analysis utilizes an autoencoder to compress galaxy spectra into a compact, redshift-invariant latent space, and a normalizing flow to identify low-probability objects. The most prominent outliers show distinctive spectral features such as irregular or double-peaked emiss… ▽ More We present an unsupervised search for outliers in the Bright Galaxy Survey (BGS) dataset from the DESI Early Data Release. This analysis utilizes an autoencoder to compress galaxy spectra into a compact, redshift-invariant latent space, and a normalizing flow to identify low-probability objects. The most prominent outliers show distinctive spectral features such as irregular or double-peaked emission lines, or originate from galaxy mergers, blended sources, and rare quasar types, including one previously unknown Broad Absorption Line system. A significant portion of the BGS outliers are stars spectroscopically misclassified as galaxies. By building our own star model trained on spectra from the DESI Milky Way Survey, we have determined that the misclassification likely stems from the Principle Component Analysis of stars in the DESI pipeline. To aid follow-up studies, we make the full probability catalog of all BGS objects and our pre-trained models publicly available. △ Less

Submitted 14 July, 2023; originally announced July 2023.

Comments: 9 pages, 5 figures

ReveaLLAGN 0: First Look at JWST MIRI data of Sombrero and NGC 1052

Authors: K. Goold, A. Seth, M. Molina, D. Ohlson, J. C. Runnoe, T. Boeker, T. A. Davis, A. Dumont, M. Eracleous, J. A. Fernández-Ontiveros, E. Gallo, A. D. Goulding, J. E. Greene, L. C. Ho, S. B. Markoff, N. Neumayer, R. Plotkin, A. Prieto, S. Satyapal, G. Van De Ven, J. L. Walsh, F. Yuan, A. Feldmeier-Krause, K. Gültekin, S. Hoenig , et al. (6 additional authors not shown)

Abstract: We present the first results from the Revealing Low-Luminosity Active Galactic Nuclei (ReveaLLAGN) survey, a JWST survey of seven nearby LLAGN. We focus on two observations with the Mid-Infrared Instrument's (MIRI) Medium Resolution Spectrograph (MRS) of the nuclei of NGC 1052 and Sombrero (NGC 4594 / M104). We also compare these data to public JWST data of a higher-luminosity AGN, NGC 7319 and NG… ▽ More We present the first results from the Revealing Low-Luminosity Active Galactic Nuclei (ReveaLLAGN) survey, a JWST survey of seven nearby LLAGN. We focus on two observations with the Mid-Infrared Instrument's (MIRI) Medium Resolution Spectrograph (MRS) of the nuclei of NGC 1052 and Sombrero (NGC 4594 / M104). We also compare these data to public JWST data of a higher-luminosity AGN, NGC 7319 and NGC 7469. JWST clearly separates the AGN spectrum from the galaxy light even in Sombrero, the faintest target in our survey; the AGN components have very red spectra. We find that the emission-line widths in both NGC 1052 and Sombrero increase with increasing ionization potential, with FWHM > 1000 km/s for lines with ionization potential > 50 eV. These lines are also significantly blue-shifted in both LLAGN. The high ionization potential lines in NGC 7319 show neither broad widths or significant blue shifts. Many of the lower ionization potential emission lines in Sombrero show significant blue wings extending > 1000 km/s. These features and the emission-line maps in both galaxies are consistent with outflows along the jet direction. Sombrero has the lowest luminosity high-ionization potential lines ([Ne V] and [O IV]) ever measured in the mid-IR, but the relative strengths of these lines are consistent with higher luminosity AGN. On the other hand, the [Ne V] emission is much weaker relative to the [Ne III] and [Ne II] lines of higher-luminosity AGN. These initial results show the great promise that JWST holds for identifying and studying the physical nature of LLAGN. △ Less

Submitted 1 March, 2024; v1 submitted 3 July, 2023; originally announced July 2023.

Comments: Accepted to ApJ Feb 28, 2024

First results from the JWST Early Release Science Program Q3D: Benchmark Comparison of Optical and Mid-IR Tracers of a Dusty, Ionized Red Quasar Wind at z=0.435

Authors: D. S. N. Rupke, D. Wylezalek, N. L. Zakamska, S. Veilleux, C. Bertemes, Y. Ishikawa, W. Liu, S. Sankar, A. Vayner, H. X. G. Lim, R. McCrory, G. Murphree, L. Whitesell, L. Shen, G. Liu, J. K. Barrera-Ballesteros, H. -W. Chen, N. Diachenko, A. D. Goulding, J. E. Greene, K. N. Hainline, F. Hamann, T. Heckman, S. D. Johnson, D. Lutz , et al. (5 additional authors not shown)

Abstract: The [OIII] 5007 A emission line is the most common tracer of warm, ionized outflows in active galactic nuclei across cosmic time. JWST newly allows us to use mid-infrared spectral features at both high spatial and spectral resolution to probe these same winds. Here we present a comparison of ground-based, seeing-limited [OIII] and space-based, diffraction-limited [SIV] 10.51 micron maps of the pow… ▽ More The [OIII] 5007 A emission line is the most common tracer of warm, ionized outflows in active galactic nuclei across cosmic time. JWST newly allows us to use mid-infrared spectral features at both high spatial and spectral resolution to probe these same winds. Here we present a comparison of ground-based, seeing-limited [OIII] and space-based, diffraction-limited [SIV] 10.51 micron maps of the powerful, kiloparsec-scale outflow in the Type 1 red quasar SDSS J110648.32+480712.3. The JWST data are from the Mid-InfraRed Instrument (MIRI). There is a close match in resolution between the datasets (0."6), in ionization potential of the O$^{+2}$ and S$^{+3}$ ions (35 eV), and in line sensitivity (1e-17 to 2e-17 erg/s/cm$^2$/arcsec$^2$). The [OIII] and [SIV] line shapes match in velocity and linewidth over much of the 20 kpc outflowing nebula, and [SIV] is the brightest line in the rest-frame 3.5-19.5 micron range, demonstrating its usefulness as a mid-IR probe of quasar outflows. [OIII] is nevertheless intriniscally brighter and provides better contrast with the point-source continuum, which is strong in the mid-IR. There is a strong anticorrelation of [OIII]/[SIV] with average velocity, which is consistent with a scenario of differential obscuration between the approaching (blueshifted) and receding (redshifted) sides of the flow. The dust in the wind may also obscure the central quasar, consistent with models that attribute red quasar extinction to dusty winds. △ Less

Submitted 11 December, 2023; v1 submitted 21 June, 2023; originally announced June 2023.

Comments: resubmitted to match published version; minor changes

Journal ref: 2023ApJ...953L..26R

UNCOVER: Candidate Red Active Galactic Nuclei at 3<z<7 with JWST and ALMA

Authors: Ivo Labbe, Jenny E. Greene, Rachel Bezanson, Seiji Fujimoto, Lukas J. Furtak, Andy D. Goulding, Jorryt Matthee, Rohan P. Naidu, Pascal A. Oesch, Hakim Atek, Gabriel Brammer, Iryna Chemerynska, Dan Coe, Sam E. Cutler, Pratika Dayal, Robert Feldmann, Marijn Franx, Karl Glazebrook, Joel Leja, Danilo Marchesini, Michael Maseda, Themiya Nanayakkara, Erica J. Nelson, Richard Pan, Casey Papovich , et al. (6 additional authors not shown)

Abstract: The James Webb Space Telescope (JWST) is revolutionizing our knowledge of $z>5$ galaxies and their actively accreting black holes. Using the JWST Cycle 1 Treasury program Ultradeep NIRSpec and NIRCam ObserVations before the Epoch of Reionization (UNCOVER) in the lensing field Abell 2744, we report the identification of a sample of little red dots at $3 < z_{\rm{phot}} < 7$ that likely contain high… ▽ More The James Webb Space Telescope (JWST) is revolutionizing our knowledge of $z>5$ galaxies and their actively accreting black holes. Using the JWST Cycle 1 Treasury program Ultradeep NIRSpec and NIRCam ObserVations before the Epoch of Reionization (UNCOVER) in the lensing field Abell 2744, we report the identification of a sample of little red dots at $3 < z_{\rm{phot}} < 7$ that likely contain highly-reddened accreting supermassive black holes. Using a NIRCam-only selection to F444W$<27.7$ mag, we find 26 sources over the $\sim45$ arcmin$^{2}$ field that are blue in F115W$-$F200W$\sim0$ (or $β_{\rm UV}\sim-2.0$ for $f_λ \propto λ^β$), red in F200W$-$F444W = $1-4$ ($β_{\rm opt} \sim +2.0$), and are dominated by a point-source like central component. Of the 20 sources with deep ALMA 1.2-mm coverage, none are detected individually or in a stack. For the majority of the sample, SED fits to the JWST+ALMA observations prefer models with hot dust rather than obscured star-formation to reproduce the red NIRCam colors and ALMA 1.2-mm non-detections. While compact dusty star formation can not be ruled out, the combination of extremely small sizes ($\langle r_e \rangle\approx50$ pc after correction for magnification), red rest-frame optical slopes, and hot dust can by explained by reddened broad-line active galactic nuclei (AGNs). Our targets have faint $M_{\rm 1450} \approx -14\ \, {\rm to} -18$ mag but inferred bolometric luminosities of $L_{\rm bol} = 10^{43}-10^{46}$ erg/s, reflecting their obscured nature. If the candidates are confirmed as AGNs with upcoming UNCOVER spectroscopy, then we have found an abundant population of reddened luminous AGN that are at least ten times more numerous than UV-luminous AGN at the same intrinsic bolometric luminosity. △ Less

Submitted 12 June, 2023; originally announced June 2023.

Comments: submitted to ApJ

Little Red Dots: an abundant population of faint AGN at z~5 revealed by the EIGER and FRESCO JWST surveys

Authors: Jorryt Matthee, Rohan P. Naidu, Gabriel Brammer, John Chisholm, Anna-Christina Eilers, Andy Goulding, Jenny Greene, Daichi Kashino, Ivo Labbe, Simon J. Lilly, Ruari Mackenzie, Pascal A. Oesch, Andrea Weibel, Stijn Wuyts, Mengyuan Xiao, Rongmon Bordoloi, Rychard Bouwens, Pieter van Dokkum, Garth Illingworth, Ivan Kramarenko, Michael V. Maseda, Charlotte Mason, Romain A. Meyer, Erica J. Nelson, Naveen A. Reddy , et al. (3 additional authors not shown)

Abstract: Characterising the prevalence and properties of faint active galactic nuclei (AGN) in the early Universe is key for understanding the formation of supermassive black holes (SMBHs) and determining their role in cosmic reionization. We perform a spectroscopic search for broad H$α$ emitters at $z\approx4-6$ using deep JWST/NIRCam imaging and wide field slitless spectroscopy from the EIGER and FRESCO… ▽ More Characterising the prevalence and properties of faint active galactic nuclei (AGN) in the early Universe is key for understanding the formation of supermassive black holes (SMBHs) and determining their role in cosmic reionization. We perform a spectroscopic search for broad H$α$ emitters at $z\approx4-6$ using deep JWST/NIRCam imaging and wide field slitless spectroscopy from the EIGER and FRESCO surveys. We identify 20 H$α$ lines at $z=4.2-5.5$ that have broad components with line widths from $\sim1200-3700$ km s$^{-1}$, contributing $\sim30-90$ % of the total line flux. We interpret these broad components as being powered by accretion onto SMBHs with implied masses $\sim10^{7-8}$ M$_{\odot}$. In the UV luminosity range M$_{\rm UV}=-21$ to $-18$, we measure number densities of $\approx10^{-5}$ cMpc$^{-3}$. This is an order of magnitude higher than expected from extrapolating quasar UV luminosity functions. Yet, such AGN are found in only $<1$ % of star-forming galaxies at $z\sim5$. The SMBH mass function agrees with large cosmological simulations. In two objects we detect narrow red- and blue-shifted H$α$ absorption indicative, respectively, of dense gas fueling SMBH growth and outflows. We may be witnessing early AGN feedback that will clear dust-free pathways through which more massive blue quasars are seen. We uncover a strong correlation between reddening and the fraction of total galaxy luminosity arising from faint AGN. This implies that early SMBH growth is highly obscured and that faint AGN are only minor contributors to cosmic reionization. △ Less

Submitted 1 February, 2024; v1 submitted 8 June, 2023; originally announced June 2023.

Comments: 26 pages, 20 figures. Accepted for publication in ApJ. Various improvements and updated citations. Main additions: colours (Fig 11), broad Halpha LF (Fig 18, Sect 5.2), comparisons to simulations (Sect 5.3)

Varstrometry for Off-nucleus and Dual sub-Kpc AGN (VODKA): Investigating the Nature of J0823+2418 at $z=1.81$: a Likely Lensed Quasar

Authors: Arran C. Gross, Yu-Ching Chen, Adi Foord, Xin Liu, Yue Shen, Masamune Oguri, Andy Goulding, Hsiang-Chih Hwang, Nadia Zakamska, Yilun Ma, Liam Nolan

Abstract: Dual quasars at small physical separations are an important precursor phase of galaxy mergers, ultimately leading to the coalescence of the two supermassive black holes. Starting from a sample of dual/lensed quasar candidates discovered using astrometric jitter in Gaia data, we present a pilot case study of one of the most promising yet puzzling candidate dual quasars at cosmic noon (z$\sim$1.8).… ▽ More Dual quasars at small physical separations are an important precursor phase of galaxy mergers, ultimately leading to the coalescence of the two supermassive black holes. Starting from a sample of dual/lensed quasar candidates discovered using astrometric jitter in Gaia data, we present a pilot case study of one of the most promising yet puzzling candidate dual quasars at cosmic noon (z$\sim$1.8). Using multi-wavelength imaging and spectroscopy from X-ray to radio, we test whether the J0823+2418 system is two individual quasars in a bound pair at separation$\sim$0.64'', or instead a single quasar being gravitationally lensed by a foreground galaxy. We find consistent flux ratios ($\sim$1.25-1.45) between the two sources in optical, NIR, UV, and radio, and thus similar spectral energy distributions, suggesting a strong lensing scenario. However, differences in the radio spectral index, as well as changing X-ray fluxes, hint at either a dual quasar with otherwise nearly identical properties, or perhaps lensing-based time lag of $\sim$3 days paired with intrinsic variability. We find with lens mass modeling that the relative NIR positions and magnitudes of the two quasars and a marginally detected central galaxy are consistent with strong lensing. Archival SDSS spectra likewise suggest a foreground absorber via Mg II absorption lines. We conclude that J0823+2418 is likely a lensed quasar, and therefore that the VODKA sample contains a population of these lensed systems (perhaps as high as 50%) as well as dual quasars. △ Less

Submitted 6 June, 2023; originally announced June 2023.

Comments: 25 pages, 14 figures, submitted to ApJ

Evidence for heavy seed origin of early supermassive black holes from a z~10 X-ray quasar

Authors: Akos Bogdan, Andy Goulding, Priyamvada Natarajan, Orsolya Kovacs, Grant Tremblay, Urmila Chadayammuri, Marta Volonteri, Ralph Kraft, William Forman, Christine Jones, Eugene Churazov, Irina Zhuravleva

Abstract: Observations of quasars reveal that many supermassive black holes (BHs) were in place less than 700 million years after the Big Bang. However, the origin of the first BHs remains a mystery. Seeds of the first BHs are postulated to be either light (i.e., $10-100~\rm{M_{\odot}})$, remnants of the first stars or heavy (i.e., $10^4-10^5~\rm{M_{\odot}})$, originating from the direct collapse of gas clo… ▽ More Observations of quasars reveal that many supermassive black holes (BHs) were in place less than 700 million years after the Big Bang. However, the origin of the first BHs remains a mystery. Seeds of the first BHs are postulated to be either light (i.e., $10-100~\rm{M_{\odot}})$, remnants of the first stars or heavy (i.e., $10^4-10^5~\rm{M_{\odot}})$, originating from the direct collapse of gas clouds. Harnessing recent data from the Chandra X-ray Observatory, we report the detection of an X-ray-luminous massive BH in a gravitationally-lensed galaxy identified by JWST at $z\approx10.3$ behind the cluster lens Abell 2744. This heavily-obscured quasar with a bolometric luminosity of $L_{\rm bol}\sim5\times10^{45}~\rm{erg\ s^{-1}}$ harbors a $M_{\rm BH}\sim10^7-10^8~\rm{M_{\odot}}$ BH assuming accretion at the Eddington limit. This mass is comparable to the inferred stellar mass of its host galaxy, in contrast to what is found in the local Universe wherein the BH mass is $\sim0.1\%$ of the host galaxy's stellar mass. The combination of such a high BH mass and large BH-to-galaxy stellar mass ratio just $\sim$500 Myrs after the Big Bang was theoretically predicted and is consistent with a picture wherein BHs originated from heavy seeds. △ Less

Submitted 25 September, 2023; v1 submitted 24 May, 2023; originally announced May 2023.

Comments: 27 pages, 6 figures, accepted

The Spitzer Coverage of HSC-Deep with IRAC for Z studies (SHIRAZ) I: IRAC mosaics

Authors: Marianna Annunziatella, Anna Sajina, Mauro Stefanon, Danilo Marchesini, Mark Lacy, Ivo Labbe, Lilianna Houston, Rachel Bezanson, Eiichi Egami, Xiaohui Fan, Duncan Farrah, Jenny Greene, Andy Goulding, Yen-Ting Lin, Xin Liu, Thibaud Moutard, Yoshiaki Ono, Masami Ouchi, Marcin Sawicki, Jason Surace, Katherine Whitaker

Abstract: We present new Spitzer Infrared Array Camera (IRAC) 3.6 and 4.5μm mosaics of three fields, E-COSMOS, DEEP2-F3, and ELAIS-N1. Our mosaics include both new IRAC observations as well as re-processed archival data in these fields. These fields are part of the HSC-Deep grizy survey and have a wealth of additional ancillary data. The addition of these new IRAC mosaics is critical in allowing for improve… ▽ More We present new Spitzer Infrared Array Camera (IRAC) 3.6 and 4.5μm mosaics of three fields, E-COSMOS, DEEP2-F3, and ELAIS-N1. Our mosaics include both new IRAC observations as well as re-processed archival data in these fields. These fields are part of the HSC-Deep grizy survey and have a wealth of additional ancillary data. The addition of these new IRAC mosaics is critical in allowing for improved photometric redshifts and stellar population parameters at cosmic noon and earlier epochs. The total area mapped by this work is {\sim} 17 deg2 with a mean integration time of {\sim}1200s, providing a median 5σ depth of 23.7(23.3) at 3.6(4.5)μm in AB. We perform SExtractor photometry both on the combined mosaics as well as the single-epoch mosaics taken {\sim}6 months apart. The resultant IRAC number counts show good agreement with previous studies. In combination with the wealth of existing and upcoming spectro-photometric data in these fields, our IRAC mosaics will enable a wide range of galactic evolution and AGN studies. With that goal in mind, we make the combined IRAC mosaics and coverage maps of these three fields publicly available. counts show good agreement with previous studies. △ Less

Submitted 22 May, 2023; originally announced May 2023.

Comments: accepted for publication in AJ

JWST UNCOVER: Discovery of $z>9$ Galaxy Candidates Behind the Lensing Cluster Abell 2744

Authors: Hakim Atek, Iryna Chemerynska, Bingjie Wang, Lukas Furtak, Andrea Weibel, Pascal Oesch, John R. Weaver, Ivo Labbé, Rachel Bezanson, Pieter van Dokkum, Adi Zitrin, Pratika Dayal, Christina C. Williams, Themiya Nannayakkara, Sedona H. Price, Gabriel Brammer, Andy D. Goulding, Joel Leja, Danilo Marchesini, Erica J. Nelson, Richard Pan, Katherine E. Whitaker

Abstract: We present the results of a search for high-redshift ($z>9$) galaxy candidates in the JWST UNCOVER survey, using deep NIRCam and NIRISS imaging in 7 bands over $\sim45$ arcmin$^2$ and ancillary HST observations. The NIRCam observations reach a $5-σ$ limiting magnitude of $\sim 29.2$ AB. The identification of high$-z$ candidates relies on a combination of a dropout selection and photometric redshif… ▽ More We present the results of a search for high-redshift ($z>9$) galaxy candidates in the JWST UNCOVER survey, using deep NIRCam and NIRISS imaging in 7 bands over $\sim45$ arcmin$^2$ and ancillary HST observations. The NIRCam observations reach a $5-σ$ limiting magnitude of $\sim 29.2$ AB. The identification of high$-z$ candidates relies on a combination of a dropout selection and photometric redshifts. We find 16 candidates at $9<z<12$ and 3 candidates at $12<z<13$, eight candidates are deemed very robust. Their lensing amplification ranges from $μ=1.2$ to 11.5. Candidates have a wide range of (lensing-corrected) luminosities and young ages, with low stellar masses ($6.8<$ log(M$_{\star}$/M$_{\odot}$) $<9.5$) and low star formation rates (SFR=0.2-7 M$_{\odot}$ yr$^{-1}$), confirming previous findings in early JWST observations of $z>9$. A few galaxies at $z\sim9-10$ appear to show a clear Balmer break between the F356W and F444W/F410M bands, which helps constrain their stellar mass. We estimate blue UV continuum slopes between $β=-1.8$ and $-2.3$, typical for early galaxies at $z>9$ but not as extreme as the bluest recently discovered sources. We also find evidence for a rapid redshift-evolution of the mass-luminosity relation and a redshift-evolution of the UV continuum slope for a given range of intrinsic magnitude, in line with theoretical predictions. These findings suggest that deeper JWST observations are needed to reach the fainter galaxy population at those early epochs, and follow-up spectroscopy will help better constrain the physical properties and star formation histories of a larger sample of galaxies. △ Less

Submitted 4 May, 2023; v1 submitted 2 May, 2023; originally announced May 2023.

Comments: Submitted to MNRAS

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)

First results from the JWST Early Release Science Program Q3D: The Warm Ionized Gas Outflow in z ~ 1.6 Quasar XID 2028 and its Impact on the Host Galaxy

Authors: Sylvain Veilleux, Weizhe Liu, Andrey Vayner, Dominika Wylezalek, David S. N. Rupke, Nadia L. Zakamska, Yuzo Ishikawa, Caroline Bertemes, Jorge K. Barrera-Ballesteros, Hsiao-Wen Chen, Nadiia Diachenko, Andy D. Goulding, Jenny E. Greene, Kevin N. Hainline, Fred Hamann, Timothy Heckman, Sean D. Johnson, Hui Xian Grace Lim, Dieter Lutz, Nora Lutzgendorf, Vincenzo Mainieri, Roberto Maiolino, Ryan McCrory, Grey Murphree, Nicole P. H. Nesvadba , et al. (4 additional authors not shown)

Abstract: Quasar feedback may regulate the growth of supermassive black holes, quench coeval star formation, and impact galaxy morphology and the circumgalactic medium. However, direct evidence for quasar feedback in action at the epoch of peak black hole accretion at z ~ 2 remains elusive. A good case in point is the z = 1.6 quasar WISEA J100211.29+013706.7 (XID 2028) where past analyses of the same ground… ▽ More Quasar feedback may regulate the growth of supermassive black holes, quench coeval star formation, and impact galaxy morphology and the circumgalactic medium. However, direct evidence for quasar feedback in action at the epoch of peak black hole accretion at z ~ 2 remains elusive. A good case in point is the z = 1.6 quasar WISEA J100211.29+013706.7 (XID 2028) where past analyses of the same ground-based data have come to different conclusions. Here we revisit this object with the integral field unit of the Near Infrared Spectrograph (NIRSpec) on board the James Webb Space Telescope (JWST) as part of Early Release Science program Q3D. The excellent angular resolution and sensitivity of the JWST data reveal new morphological and kinematic sub-structures in the outflowing gas plume. An analysis of the emission line ratios indicates that photoionization by the central quasar dominates the ionization state of the gas with no obvious sign for a major contribution from hot young stars anywhere in the host galaxy. Rest-frame near-ultraviolet emission aligned along the wide-angle cone of outflowing gas is interpreted as a scattering cone. The outflow has cleared a channel in the dusty host galaxy through which some of the quasar ionizing radiation is able to escape and heat the surrounding interstellar and circumgalactic media. The warm ionized outflow is not powerful enough to impact the host galaxy via mechanical feedback, but radiative feedback by the AGN, aided by the outflow, may help explain the unusually small molecular gas mass fraction in the galaxy host. △ Less

Submitted 22 June, 2023; v1 submitted 15 March, 2023; originally announced March 2023.

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

First results from the JWST Early Release Science Program Q3D: Ionization cone, clumpy star formation and shocks in a $z=3$ extremely red quasar host

Authors: Andrey Vayner, Nadia L. Zakamska, Yuzo Ishikawa, Swetha Sankar, Dominika Wylezalek, David S. N. Rupke, Sylvain Veilleux, Caroline Bertemes, Jorge K. Barrera-Ballesteros, Hsiao-Wen Chen, Nadiia Diachenko, Andy D. Goulding, Jenny E. Greene, Kevin N. Hainline, Fred Hamann, Timothy Heckman, Sean D. Johnson, Hui Xian Grace Lim, Weizhe Liu, Dieter Lutz, Nora Lutzgendorf, Vincenzo Mainieri, Ryan McCrory, Grey Murphree, Nicole P. H. Nesvadba , et al. (3 additional authors not shown)

Abstract: Massive galaxies formed most actively at redshifts $z=1-3$ during the period known as `cosmic noon.' Here we present an emission-line study of an extremely red quasar SDSSJ165202.64+172852.3 host galaxy at $z=2.94$, based on observations with the Near Infrared Spectrograph (NIRSpec) integral field unit (IFU) on board JWST. We use standard emission-line diagnostic ratios to map the sources of gas i… ▽ More Massive galaxies formed most actively at redshifts $z=1-3$ during the period known as `cosmic noon.' Here we present an emission-line study of an extremely red quasar SDSSJ165202.64+172852.3 host galaxy at $z=2.94$, based on observations with the Near Infrared Spectrograph (NIRSpec) integral field unit (IFU) on board JWST. We use standard emission-line diagnostic ratios to map the sources of gas ionization across the host and a swarm of companion galaxies. The quasar dominates the photoionization, but we also discover shock-excited regions orthogonal to the ionization cone and the quasar-driven outflow. These shocks could be merger-induced or -- more likely, given the presence of a powerful galactic-scale quasar outflow -- these are signatures of wide-angle outflows that can reach parts of the galaxy that are not directly illuminated by the quasar. Finally, the kinematically narrow emission associated with the host galaxy presents as a collection of 1 kpc-scale clumps forming stars at a rate of at least 200 $M_{\odot}$ yr$^{-1}$. The ISM within these clumps shows high electron densities, reaching up to 3,000 cm$^{-3}$ with metallicities ranging from half to a third solar with a positive metallicity gradient and V band extinctions up to 3 magnitudes. The star formation conditions are far more extreme in these regions than in local star-forming galaxies but consistent with that of massive galaxies at cosmic noon. JWST observations reveal an archetypical rapidly forming massive galaxy undergoing a merger, a clumpy starburst, an episode of obscured near-Eddington quasar activity, and an extremely powerful quasar outflow simultaneously. △ Less

Submitted 25 July, 2023; v1 submitted 13 March, 2023; originally announced March 2023.

Comments: 19 pages, 8 figures. Accepted for publication in ApJ

Beyond Ultra-diffuse Galaxies. II. Environmental Quenching of Mass-Size Outliers among the Satellites of Milky Way Analogs

Authors: Jiaxuan Li, Jenny E. Greene, Johnny Greco, Rachael Beaton, Shany Danieli, Andy Goulding, Song Huang, Erin Kado-Fong

Abstract: Recent observations have reignited interest in a population of dwarf galaxies that are large and diffuse for their mass, often called ultra-diffuse galaxies. However, the origin and evolution of these mass-size outliers and the role of the environment are still unclear. Using the exquisitely deep and wide Hyper Suprime-Cam Strategic Survey images, we search for ultra-puffy galaxies, defined as bei… ▽ More Recent observations have reignited interest in a population of dwarf galaxies that are large and diffuse for their mass, often called ultra-diffuse galaxies. However, the origin and evolution of these mass-size outliers and the role of the environment are still unclear. Using the exquisitely deep and wide Hyper Suprime-Cam Strategic Survey images, we search for ultra-puffy galaxies, defined as being $1.5σ$ larger than the average size for their mass, around Milky Way-like galaxies. We present the sizes and radial distributions of mass-size outliers and derive their quenched fraction to explore the impact of the environment. Surprisingly, despite being outliers in size, the ultra-puffy galaxies have a similar quenched fraction as normal-sized satellites of Milky Way analogs in both observations and simulations, suggesting that quenching is not tied to being a mass-size outlier. The quenched fraction is higher for the ultra-puffy galaxies associated with redder hosts, as well as those that are closer to the host in projection. In contrast, the ultra-diffuse galaxies are overall redder and more quiescent compared with normal satellites. We show that the classic definition of ultra-diffuse galaxies is heavily weighted toward quenched galaxies and thus cannot be used for a study of quenching of mass-size outliers. △ Less

Submitted 11 September, 2023; v1 submitted 27 February, 2023; originally announced February 2023.

Comments: 13 pages, 6 figures, ApJ accepted. Key results are shown in Figures 4-6

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

Observational Signatures of Carbon-Oxygen White Dwarf Merger Remnants

Authors: Philippe Z. Yao, Eliot Quataert, Andy Goulding

Abstract: Many double white dwarf (WD) mergers likely do not lead to a prompt thermonuclear explosion. We investigate the prospects for observationally detecting the surviving remnants of such mergers, focusing on the case of mergers of double Carbon-Oxygen WDs. For $\sim 10^4$ yr, the merger remnant is observationally similar to an extreme AGB star evolving to become a massive WD. Identifying merger remnan… ▽ More Many double white dwarf (WD) mergers likely do not lead to a prompt thermonuclear explosion. We investigate the prospects for observationally detecting the surviving remnants of such mergers, focusing on the case of mergers of double Carbon-Oxygen WDs. For $\sim 10^4$ yr, the merger remnant is observationally similar to an extreme AGB star evolving to become a massive WD. Identifying merger remnants is thus easiest in galaxies with high stellar masses (high WD merger rate) and low star formation rates (low birth rate of $\sim 6-10 \,{\rm M_{\odot}}$ stars). Photometrically identifying merger remnants is challenging even in these cases because the merger remnants appear similar to He stars and post-outburst classical novae. We propose that the most promising technique for discovering WD merger remnants is through their unusual surrounding photoionized nebulae. We use CLOUDY photoionization calculations to investigate their unique spectral features. Merger remnants should produce weak hydrogen lines and strong carbon and oxygen recombination and fine-structure lines in the UV, optical and IR. With narrow-band imaging or integral field spectrographs, we predict that multiple candidates are detectable in the bulge of M31, the outskirts of M87 and other nearby massive galaxies, and the Milky Way. Our models roughly reproduce the WISE nebula surrounding the Galactic WD merger candidate IRAS 00500+6713; we predict detectable [Ne\,VI] and [Mg\,VII] lines with JWST but that the mid-IR WISE emission is dominated by dust not fine-structure lines. △ Less

Submitted 29 June, 2023; v1 submitted 15 February, 2023; originally announced February 2023.

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

Autoencoding Galaxy Spectra II: Redshift Invariance and Outlier Detection

Authors: Yan Liang, Peter Melchior, Sicong Lu, Andy Goulding, Charlotte Ward

Abstract: We present an unsupervised outlier detection method for galaxy spectra based on the spectrum autoencoder architecture spender, which reliably captures spectral features and provides highly realistic reconstructions for SDSS galaxy spectra. We interpret the sample density in the autoencoder latent space as a probability distribution, and identify outliers as low-probability objects with a normalizi… ▽ More We present an unsupervised outlier detection method for galaxy spectra based on the spectrum autoencoder architecture spender, which reliably captures spectral features and provides highly realistic reconstructions for SDSS galaxy spectra. We interpret the sample density in the autoencoder latent space as a probability distribution, and identify outliers as low-probability objects with a normalizing flow. However, we found that the latent-space position is not, as expected from the architecture, redshift invariant, which introduces stochasticity into the latent space and the outlier detection method. We solve this problem by adding two novel loss terms during training, which explicitly link latent-space distances to data-space distances, preserving locality in the autoencoding process. Minimizing the additional losses leads to a redshift-invariant, non-degenerate latent space distribution with clear separations between common and anomalous data. We inspect the spectra with the lowest probability and find them to include blends with foreground stars, extremely reddened galaxies, galaxy pairs and triples, and stars that are misclassified as galaxies. We release the newly trained spender model and the latent-space probability for the entire SDSS-I galaxy sample to aid further investigations. △ Less

Submitted 5 February, 2023; originally announced February 2023.

Comments: 12 pages, 7 figures, submitted to AJ. Code available at https://github.com/pmelchior/spender

Merger Signatures are Common, but not Universal, In Massive, Recently-Quenched Galaxies at z~0.7

Authors: Margaret Verrico, David J. Setton, Rachel Bezanson, Jenny E. Greene, Katherine A. Suess, Andy D. Goulding, Justin S. Spilker, Mariska Kriek, Robert Feldmann, Desika Narayanan, Vincenzo Donofrio, Gourav Khullar

Abstract: We present visual classifications of merger-induced tidal disturbances in 143 $\rm{M}_* \sim 10^{11}\rm{M}_\odot$ post-starburst galaxies at z$\sim$0.7 identified in the SQuIGG$\vec{L}$E Sample. This sample spectroscopically selects galaxies from the Sloan Digital Sky Survey that have stopped their primary epoch of star formation within the past $\sim$500 Myrs. Visual classifications are performed… ▽ More We present visual classifications of merger-induced tidal disturbances in 143 $\rm{M}_* \sim 10^{11}\rm{M}_\odot$ post-starburst galaxies at z$\sim$0.7 identified in the SQuIGG$\vec{L}$E Sample. This sample spectroscopically selects galaxies from the Sloan Digital Sky Survey that have stopped their primary epoch of star formation within the past $\sim$500 Myrs. Visual classifications are performed on Hyper Suprime Cam (HSC) i-band imaging. We compare to a control sample of mass- and redshift-matched star-forming and quiescent galaxies from the Large Early Galaxy Census and find that post-starburst galaxies are more likely to be classified as disturbed than either category. This corresponds to a factor of $3.6^{+2.9}_{-1.3}$ times the disturbance rate of older quiescent galaxies and $2.1^{+1.9}_{-.73}$ times the disturbance rate of star-forming galaxies. Assuming tidal features persist for $\lesssim500$ Myr, this suggests merging is coincident with quenching in a significant fraction of these post-starbursts. Galaxies with tidal disturbances are younger on average than undisturbed post-starburst galaxies in our sample, suggesting tidal features from a major merger may have faded over time. This may be exacerbated by the fact that, on average, the undisturbed subset is fainter, rendering low surface brightness tidal features harder to identify. However, the presence of ten young ($\lesssim150$ Myr since quenching) undisturbed galaxies suggests that major mergers are not the only fast physical mechanism that shut down the primary epoch of star formation in massive galaxies at intermediate redshift. △ Less

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

Comments: 13 pages, 8 figures, accepted to ApJ

Autoencoding Galaxy Spectra I: Architecture

Authors: Peter Melchior, Yan Liang, ChangHoon Hahn, Andy Goulding

Abstract: We introduce the neural network architecture SPENDER as a core differentiable building block for analyzing, representing, and creating galaxy spectra. It combines a convolutional encoder, which pays attention to up to 256 spectral features and compresses them into a low-dimensional latent space, with a decoder that generates a restframe representation, whose spectral range and resolution exceeds t… ▽ More We introduce the neural network architecture SPENDER as a core differentiable building block for analyzing, representing, and creating galaxy spectra. It combines a convolutional encoder, which pays attention to up to 256 spectral features and compresses them into a low-dimensional latent space, with a decoder that generates a restframe representation, whose spectral range and resolution exceeds that of the observing instrument. The decoder is followed by explicit redshift, resampling, and convolution transformations to match the observations. The architecture takes galaxy spectra at arbitrary redshifts and is robust to glitches like residuals of the skyline subtraction, so that spectra from a large survey can be ingested directly without additional preprocessing. We demonstrate the performance of SPENDER by training on the entire spectroscopic galaxy sample of SDSS-II; show its ability to create highly accurate reconstructions with substantially reduced noise; perform deconvolution and oversampling for a super-resolution model that resolves the [OII] doublet; introduce a novel method to interpret attention weights as proxies for important spectral features; and infer the main degrees of freedom represented in the latent space. We conclude with a discussion of future improvements and applications. △ Less

Submitted 14 November, 2022; originally announced November 2022.

Comments: 17 pages, 5 figures, submitted to AJ. Code available at https://github.com/pmelchior/spender

Beyond Ultra-diffuse Galaxies. I. Mass--Size Outliers among the Satellites of Milky Way Analogs

Authors: Jiaxuan Li, Jenny E. Greene, Johnny P. Greco, Song Huang, Peter Melchior, Rachael Beaton, Kirsten Casey, Shany Danieli, Andy Goulding, Remy Joseph, Erin Kado-Fong, Ji Hoon Kim, Lauren A. MacArthur

Abstract: Large diffuse galaxies are hard to find, but understanding the environments where they live, their numbers, and ultimately their origins, is of intense interest and importance for galaxy formation and evolution. Using Subaru's Hyper Suprime-Cam Strategic Survey Program, we perform a systematic search for low surface brightness galaxies and present novel and effective methods for detecting and mode… ▽ More Large diffuse galaxies are hard to find, but understanding the environments where they live, their numbers, and ultimately their origins, is of intense interest and importance for galaxy formation and evolution. Using Subaru's Hyper Suprime-Cam Strategic Survey Program, we perform a systematic search for low surface brightness galaxies and present novel and effective methods for detecting and modeling them. As a case study, we surveyed 922 Milky Way analogs in the nearby Universe ($0.01 < z < 0.04$) and build a large sample of satellite galaxies that are outliers in the mass-size relation. These ``ultra-puffy'' galaxies (UPGs), defined to be $1.5σ$ above the average mass-size relation, represent the tail of the satellite size distribution. We find that each MW analog hosts $N_{\rm UPG} = 0.31\pm 0.05$ ultra-puffy galaxies on average, which is consistent with but slightly lower than the observed abundance at this halo mass in the Local Volume. We also construct a sample of ultra-diffuse galaxies (UDGs) in MW analogs and find an abundance of $N_{\rm UDG} = 0.44\pm0.05$ per host. With literature results, we confirm that the UDG abundance scales with the host halo mass following a sublinear power law. We argue that our definition for ultra-puffy galaxies, which is based on the mass-size relation, is more physically-motivated than the common definition of ultra-diffuse galaxies, which depends on surface brightness and size cuts and thus yields different surface mass density cuts for quenched and star-forming galaxies. △ Less

Submitted 11 September, 2023; v1 submitted 26 October, 2022; originally announced October 2022.

Comments: 19 pages, 7 figures, ApJ accepted. Also see the companion paper arXiv:2302.14108

ELVES IV: The Satellite Stellar-to-Halo Mass Relation Beyond the Milky-Way

Authors: Shany Danieli, Jenny E. Greene, Scott Carlsten, Fangzhou Jiang, Rachael Beaton, Andy D. Goulding

Abstract: Quantifying the connection between galaxies and their host dark matter halos has been key for testing cosmological models on various scales. Below $M_\star \sim 10^9\,M_\odot$, such studies have primarily relied on the satellite galaxy population orbiting the Milky Way. Here we present new constraints on the connection between satellite galaxies and their host dark matter subhalos using the larges… ▽ More Quantifying the connection between galaxies and their host dark matter halos has been key for testing cosmological models on various scales. Below $M_\star \sim 10^9\,M_\odot$, such studies have primarily relied on the satellite galaxy population orbiting the Milky Way. Here we present new constraints on the connection between satellite galaxies and their host dark matter subhalos using the largest sample of satellite galaxies in the Local Volume ($D \lesssim 12\,\mathrm{Mpc}$) to date. We use $250$ confirmed and $71$ candidate dwarf satellites around 27 Milky Way (MW)-like hosts from the Exploration of Local VolumE Satellites (ELVES) Survey and use the semi-analytical SatGen model for predicting the population of dark matter subhalos expected in the same volume. Through a Bayesian model comparison of the observed and the forward-modeled satellite stellar mass functions (SSMF), we infer the satellite stellar-to-halo mass relation. We find that the observed SSMF is best reproduced when subhalos at the low mass end are populated by a relation of the form $M_\star \propto M^α_\mathrm{peak}$, with a moderate slope of $α_\mathrm{const}=2.10 \pm 0.01$ and a low scatter, constant as a function of the peak halo mass, of $σ_\mathrm{const}=0.06^{+0.07}_{-0.05}$. A model with a steeper slope ($α_\mathrm{grow}=2.39 \pm 0.06$) and a scatter that grows with decreasing $M_\mathrm{peak}$ is also consistent with the observed SSMF but is not required. Our new model for the satellite-subhalo connection, based on hundreds of Local Volume satellite galaxies, is in line with what was previously derived using only the Milky Way satellites. △ Less

Submitted 13 August, 2023; v1 submitted 25 October, 2022; originally announced October 2022.

Comments: Accepted for publication in ApJ. Figure 8 shows the key result -- the Satellite Stellar to Halo Mass relation obtained in this work, in comparison to previous studies

First results from the JWST Early Release Science Program Q3D: Turbulent times in the life of a $z \sim 3$ extremely red quasar revealed by NIRSpec IFU

Authors: Dominika Wylezalek, Andrey Vayner, David S. N. Rupke, Nadia L. Zakamska, Sylvain Veilleux, Yuzo Ishikawa, Caroline Bertemes, Weizhe Liu, Jorge K. Barrera-Ballesteros, Hsiao-Wen Chen, Andy D. Goulding, Jenny E. Greene, Kevin N. Hainline, Nora Lützgendorf, Fred Hamann, Timothy Heckman, Sean D. Johnson, Dieter Lutz, Vincenzo Mainieri, Roberto Maiolino, Nicole P. H. Nesvadba, Patrick Ogle, Eckhard Sturm

Abstract: Extremely red quasars, with bolometric luminosities exceeding $10^{47}$ erg s$^{-1}$, are a fascinating high-redshift population that is absent in the local universe. They are the best candidates for supermassive black holes accreting at rates at or above the Eddington limit, and they are associated with the most rapid and powerful outflows of ionized gas known to date. They are also hosted by mas… ▽ More Extremely red quasars, with bolometric luminosities exceeding $10^{47}$ erg s$^{-1}$, are a fascinating high-redshift population that is absent in the local universe. They are the best candidates for supermassive black holes accreting at rates at or above the Eddington limit, and they are associated with the most rapid and powerful outflows of ionized gas known to date. They are also hosted by massive galaxies. Here we present the first integral field unit (IFU) observations of a high-redshift quasar obtained by the Near Infrared Spectrograph (NIRSpec) on board the James Webb Space Telescope (JWST), which targeted SDSSJ165202.64+172852.3, an extremely red quasar at $z=2.94$. JWST observations reveal extended ionized gas - as traced by [OIII]$λ$5007Å- in the host galaxy of the quasar, its outflow, and the circumgalactic medium. The complex morphology and kinematics imply that the quasar resides in a very dense environment with several interacting companion galaxies within projected distances of 10-15 kpc. The high density of the environment and the large velocities of the companion galaxies suggest that this system may represent the core of a forming cluster of galaxies. The system is a good candidate for a merger of two or more dark matter halos, each with a mass of a few $10^{13}$ M$_\odot$ and traces potentially one of the densest knots at $z\sim3$. △ Less

Submitted 18 October, 2022; originally announced October 2022.

Comments: ApJL in press, 11 pages, please also see associated ESA / STScI press release from Oct 20th 2022

Ultra-Diffuse Galaxies as Extreme Star-forming Environments I: Mapping Star Formation in HI-Rich UDGs

Authors: Erin Kado-Fong, Jenny E. Greene, Song Huang, Andy Goulding

Abstract: Ultra-Diffuse Galaxies are both extreme products of galaxy evolution and extreme environments in which to test our understanding of star formation. In this work, we contrast the spatially resolved star formation activity of a sample of 22 HI-selected UDGs and 35 low-mass galaxies from the NASA Sloan Atlas (NSA) within 120 Mpc. We employ a new joint SED fitting method to compute star formation rate… ▽ More Ultra-Diffuse Galaxies are both extreme products of galaxy evolution and extreme environments in which to test our understanding of star formation. In this work, we contrast the spatially resolved star formation activity of a sample of 22 HI-selected UDGs and 35 low-mass galaxies from the NASA Sloan Atlas (NSA) within 120 Mpc. We employ a new joint SED fitting method to compute star formation rate and stellar mass surface density maps that leverage the high spatial resolution optical imaging data of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) and the UV coverage of GALEX, along with HI radial profiles estimated from a subset of galaxies that have spatially resolved HI maps. We find that the UDGs have low star formation efficiencies as a function of their atomic gas down to scales of 500 pc. We additionally find that the stellar mass-weighted sizes of our UDG sample are unremarkable when considered as a function of their HI mass -- their stellar sizes are comparable to the NSA dwarfs at fixed HI mass. This is a natural result in the picture where UDGs are forming stars normally, but at low efficiencies. We compare our results to predictions from contemporary models of galaxy formation, and find in particular that our observations are difficult to reproduce in models where UDGs undergo stellar expansion due to vigorous star formation feedback should bursty star formation be required down to $z=0$. △ Less

Submitted 24 October, 2022; v1 submitted 12 September, 2022; originally announced September 2022.

Comments: Accepted to ApJ, 27 pages, 18 figures

An empirical approach to selecting the first growing black hole seeds with JWST/NIRCam

Authors: Andy D. Goulding, Jenny E. Greene

Abstract: The James Webb Space Telescope (JWST) will have the sensitivity to detect early low-mass black holes as they transition from "seeds" to supermassive black holes (BHs). Based on the JAGUAR mock catalog of galaxies, we present a clean color selection that takes advantage of the unique UV slope of accreting supermassive black holes with relatively low mass and high accretion rates. We show for a rang… ▽ More The James Webb Space Telescope (JWST) will have the sensitivity to detect early low-mass black holes as they transition from "seeds" to supermassive black holes (BHs). Based on the JAGUAR mock catalog of galaxies, we present a clean color selection that takes advantage of the unique UV slope of accreting supermassive black holes with relatively low mass and high accretion rates. We show for a range of host stellar masses that those galaxies hosting $\sim 10^{6} M_{\odot}$ BHs radiating at >10% of their Eddington luminosity separate in color space from inactive systems. Here we propose a set of 3-band, 2-color selection boxes (with 90% completeness; 90% purity; balanced purity/completeness) with JWST/NIRCam to identify the most promising growing BH candidates at $z \sim 7-10$. △ Less

Submitted 16 March, 2023; v1 submitted 4 August, 2022; originally announced August 2022.

Comments: Updated to accepted version. 11 pages, 4 figures, 1 table. Published in ApJL

Spectropolarimetric measurements of hidden broad lines in nearby megamaser galaxies: a lack of clear evidence for a correlation between black hole masses and virial products

Authors: Nora B. Linzer, Andy D. Goulding, Jenny E. Greene, Ryan C. Hickox

Abstract: High-accuracy black hole (BH) masses require excellent spatial resolution that is only achievable for galaxies within ~100 Mpc using present-day technology. At larger distances, BH masses are often estimated with single-epoch scaling relations for active galactic nuclei. This method requires only luminosity and the velocity dispersion of the broad line region (BLR) to calculate a virial product, a… ▽ More High-accuracy black hole (BH) masses require excellent spatial resolution that is only achievable for galaxies within ~100 Mpc using present-day technology. At larger distances, BH masses are often estimated with single-epoch scaling relations for active galactic nuclei. This method requires only luminosity and the velocity dispersion of the broad line region (BLR) to calculate a virial product, and an additional virial factor, $f$, to determine BH mass. The accuracy of these single-epoch masses, however, is unknown, and there are few empirical constraints on the variance of $f$ between objects. We attempt to calibrate single-epoch BH masses using spectropolarimetric measurements of nine megamaser galaxies from which we measure the velocity distribution of the BLR. We do not find strong evidence for a correlation between the virial products used for single-epoch masses and dynamical mass, neither for the megamaser sample alone or when combined with dynamical masses from reverberation mapping modeling. Furthermore, we find evidence that the virial parameter $f$ varies between objects, but we do not find strong evidence for a correlation with other observable parameters such as luminosity or broad line width. Although we cannot definitively rule out the existence of any correlation between dynamical mass and virial product, we find tension between allowed $f$ values for masers and those widely used in the literature. We conclude that the single-epoch method requires further investigation if it is to be used successfully to infer BH masses. △ Less

Submitted 18 October, 2022; v1 submitted 2 August, 2022; originally announced August 2022.

Comments: 27 pages, 17 figures, accepted by ApJ

The Nature of Low Surface Brightness Galaxies in the Hyper Suprime-Cam Survey

Authors: Jenny E Greene, Johnny P. Greco, Andy D. Goulding, Song Huang, Erin Kado-Fong, Shany Danieli, Jiaxuan Li, Ji Hoon Kim, Yutaka Komiyama, Alexie Leauthaud, Lauren A. MacArthur, Cristobal Sifon

Abstract: We present the statistical redshift distribution of a large sample of low surface brightness (LSB) galaxies identified in the first 200 deg$^2$ of the Hyper Suprime-Cam Strategic Survey Program. Through cross-correlation with the NASA-SDSS Atlas, we find that the majority of objects lie within z<0.15 or ~ 500 Mpc, yielding a mass range of $M_*$ ~ $10^7-10^9$ M_sun and size range of… ▽ More We present the statistical redshift distribution of a large sample of low surface brightness (LSB) galaxies identified in the first 200 deg$^2$ of the Hyper Suprime-Cam Strategic Survey Program. Through cross-correlation with the NASA-SDSS Atlas, we find that the majority of objects lie within z<0.15 or ~ 500 Mpc, yielding a mass range of $M_*$ ~ $10^7-10^9$ M_sun and size range of $r_{\rm{eff,g}}$ ~ 1-8 kpc. We find a peak in the distance distribution within 100 Mpc, corresponding mostly to ~ $10^7$ M_sun galaxies that fall on the known mass-size relation. There is also a tail in the redshift distribution out to z~0.15, comprising more massive ($M_*=10^8-10^9$ M_sun) galaxies at the larger end of our size range. We see tentative evidence that at the higher-mass end ($M_* > 10^8$ M_sun) the LSB galaxies do not form a smooth extension of the mass-size relation of higher surface brightness galaxies, perhaps suggesting that the LSB galaxy population is distinct in its formation path. △ Less

Submitted 25 April, 2022; originally announced April 2022.

Comments: 15 pages, 6 figures, resubmitted to ApJ after review

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