Showing 1–50 of 132 results for author: Feldmann, R

Elevated UV luminosity density at Cosmic Dawn explained by non-evolving, weakly-mass dependent star formation efficiency

Authors: Robert Feldmann, Michael Boylan-Kolchin, James S. Bullock, Onur Çatmabacak, Claude-André Faucher-Giguère, Christopher C. Hayward, Dušan Kereš, Alexandres Lazar, Lichen Liang, Jorge Moreno, Pascal A. Oesch, Eliot Quataert, Xuejian Shen, Guochao Sun

Abstract: Recent observations with the James Webb Space Telescope (JWST) have uncovered unexpectedly high cosmic star formation activity in the early Universe, mere hundreds of millions of years after the Big Bang. These observations are often understood to reflect an evolutionary shift in star formation efficiency (SFE) caused by changing galactic conditions during these early epochs. We present FIREbox-HR… ▽ More Recent observations with the James Webb Space Telescope (JWST) have uncovered unexpectedly high cosmic star formation activity in the early Universe, mere hundreds of millions of years after the Big Bang. These observations are often understood to reflect an evolutionary shift in star formation efficiency (SFE) caused by changing galactic conditions during these early epochs. We present FIREbox-HR, a high-resolution, cosmological hydrodynamical simulation from the Feedback in Realistic Environments project, which offers insights into the SFE of galaxies during the first billion years of cosmic time. FIREbox-HR re-simulates the cosmic volume (L = 22.1 cMpc) of the original FIREbox run with eight times higher mass resolution (m_b ~ 7800 M_sun), but with identical physics, down to z ~ 6. FIREbox-HR predicts ultraviolet (UV) luminosity functions in good agreement with available observational data. The simulation also successfully reproduces the observed cosmic UV luminosity density at z ~ 6 - 14, demonstrating that relatively high star formation activity in the early Universe is a natural outcome of the baryonic processes encoded in the FIRE-2 model. According to FIREbox-HR, the SFE - halo mass relation for intermediate mass halos (M_halo ~ 10^9 - 10^11 M_sun) does not significantly evolve with redshift and is only weakly mass-dependent. These properties of the SFE - halo mass relation lead to a larger contribution from lower mass halos at higher z, driving the gradual evolution of the observed cosmic UV luminosity density. A theoretical model based on the SFE - halo mass relation inferred from FIREbox-HR allows us to explore implications for galaxy evolution. Future observations of UV faint galaxies at z > 12 will provide an opportunity to further test these predictions and deepen our understanding of star formation during Cosmic Dawn. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: 26 pages, 14 figures, 5 tables, submitted to MNRAS, comments welcome

Renormalized Internally-Contracted Multireference Coupled Cluster with Perturbative Triples

Authors: Robin Feldmann, Markus Reiher

Abstract: In this work, we combine the many-body formulation of the internally contracted multireference coupled cluster (ic-MRCC) method with Evangelista's multireference formulation of the driven similarity renormalization group (DSRG). The DSRG method can be viewed as a unitary multireference coupled cluster theory, which renormalizes the amplitudes based on a flow equation approach to eliminate numerica… ▽ More In this work, we combine the many-body formulation of the internally contracted multireference coupled cluster (ic-MRCC) method with Evangelista's multireference formulation of the driven similarity renormalization group (DSRG). The DSRG method can be viewed as a unitary multireference coupled cluster theory, which renormalizes the amplitudes based on a flow equation approach to eliminate numerical instabilities. We extend this approach by demonstrating that the unitary flow equation approach can be adapted for nonunitary transformations, rationalizing the renormalization of ic-MRCC amplitudes. We denote the new approach, the renormalized ic-MRCC (ric-MRCC) method. To achieve high accuracy with a reasonable computational cost, we introduce a new approximation to the Baker-Campbell-Hausdorff expansion. We fully consider the linear commutator while approximating the quadratic commutator, for which we neglect specific contractions involving amplitudes with active indices. Moreover, we introduce approximate perturbative triples to obtain the ric-MRCCSD[T] method. We demonstrate the accuracy of our approaches in comparison to advanced multireference methods for the potential energy curves of H8, F2, H2O, N2, and Cr2. Additionally, we show that ric-MRCCSD and ric-MRCSSD[T] match the accuracy of CCSD(T) for evaluating spectroscopic constants and of full configuration interaction energies for a set of small molecules. △ Less

Submitted 25 May, 2024; originally announced May 2024.

Comments: 59 pages, 4 figures, 6 tables

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/

Complete Active Space Iterative Coupled Cluster Theory

Authors: Robin Feldmann, Max Mörchen, Jakub Lang, Michał\ Lesiuk, Markus Reiher

Abstract: In this work, we investigate the possibility of improving multireference-driven coupled cluster (CC) approaches with an algorithm that iteratively combines complete active space (CAS) calculations with tailored CC and externally corrected CC. This is accomplished by establishing a feedback loop between the CC and CAS parts of a calculation through similarity transformation of the Hamiltonian with… ▽ More In this work, we investigate the possibility of improving multireference-driven coupled cluster (CC) approaches with an algorithm that iteratively combines complete active space (CAS) calculations with tailored CC and externally corrected CC. This is accomplished by establishing a feedback loop between the CC and CAS parts of a calculation through similarity transformation of the Hamiltonian with those CC amplitudes that are not encompassed by the active space. We denote this approach the complete active space iterative coupled cluster (CASiCC) ansatz. We investigate its efficiency and accuracy in the singles and doubles approximation by studying the prototypical molecules H4, H8, H2O, and N2. Our results demonstrate that CASiCC systematically improves on the single-reference CCSD and the ecCCSD methods across entire potential energy curves, while retaining modest computational costs. However, the tailored coupled cluster method shows superior performance in the strong correlation regime suggesting that its accuracy is based on error compensation. We find that the iterative version of externally corrected and tailored coupled cluster methods converge to the same results. △ Less

Submitted 9 April, 2024; originally announced April 2024.

Comments: 38 pages, 7 figures, 1 table

Predictions for CO emission and the CO-to-H$_2$ conversion factor in galaxy simulations with non-equilibrium chemistry

Authors: Oliver A. Thompson, Alexander J. Richings, Brad K. Gibson, Claude-André Faucher-Giguère, Robert Feldmann, Christopher C. Hayward

Abstract: Our ability to trace the star-forming molecular gas is important to our understanding of the Universe. We can trace this gas using CO emission, converting the observed CO intensity into the H$_2$ gas mass of the region using the CO-to-H$_2$ conversion factor (Xco). In this paper, we use simulations to study the conversion factor and the molecular gas within galaxies. We analysed a suite of simulat… ▽ More Our ability to trace the star-forming molecular gas is important to our understanding of the Universe. We can trace this gas using CO emission, converting the observed CO intensity into the H$_2$ gas mass of the region using the CO-to-H$_2$ conversion factor (Xco). In this paper, we use simulations to study the conversion factor and the molecular gas within galaxies. We analysed a suite of simulations of isolated disc galaxies, ranging from dwarfs to Milky Way-mass galaxies, that were run using the FIRE-2 subgrid models coupled to the CHIMES non-equilibrium chemistry solver. We use the non-equilibrium abundances from the simulations, and we also compare to results using abundances assuming equilibrium, which we calculate from the simulation in post-processing. Our non-equilibrium simulations are able to reproduce the relation between CO and H$_2$ column densities, and the relation between Xco and metallicity, seen within observations of the Milky Way. We also compare to the xCOLD GASS survey, and find agreement with their data to our predicted CO luminosities at fixed star formation rate. We also find the multivariate function used by xCOLD GASS overpredicts the H$_2$ mass for our simulations, motivating us to suggest an alternative multivariate function of our fitting, though we caution that this fitting is uncertain due to the limited range of galaxy conditions covered by our simulations. We also find that the non-equilibrium chemistry has little effect on the conversion factor (<5\%) for our high-mass galaxies, though still affects the H$_2$ mass and Lco by $\approx$25\%. △ Less

Submitted 12 June, 2024; v1 submitted 21 March, 2024; originally announced March 2024.

Comments: 19 pages, 16 figures (including appendices). Accepted for publication in MNRAS

The High-Redshift Gas-Phase Mass-Metallicity Relation in FIRE-2

Authors: Andrew Marszewski, Guochao Sun, Claude-André Faucher-Giguère, Christopher C. Hayward, Robert Feldmann

Abstract: The unprecedented infrared spectroscopic capabilities of JWST have provided high-quality interstellar medium (ISM) metallicity measurements and enabled characterization of the gas-phase mass-metallicity relation (MZR) for galaxies at $z \gtrsim 5$ for the first time. We analyze the gas-phase MZR and its evolution in a high-redshift suite of FIRE-2 cosmological zoom-in simulations at $z=5-12$ and f… ▽ More The unprecedented infrared spectroscopic capabilities of JWST have provided high-quality interstellar medium (ISM) metallicity measurements and enabled characterization of the gas-phase mass-metallicity relation (MZR) for galaxies at $z \gtrsim 5$ for the first time. We analyze the gas-phase MZR and its evolution in a high-redshift suite of FIRE-2 cosmological zoom-in simulations at $z=5-12$ and for stellar masses $M_* \sim 10^6-10^{10} \rm{M}_\odot$. These simulations implement a multi-channel stellar feedback model and produce broadly realistic galaxy properties, including when evolved to $z=0$. The simulations predict very weak redshift evolution of the MZR over the redshift range studied, with the normalization of the MZR increasing by less than $0.01$ dex as redshift decreases from $z = 12$ to $z=5$. The median MZR in the simulations is well-approximated as a constant power-law relation across this redshift range given by $\log(Z/Z_\odot) = 0.37\log(M_*/\rm{M}_\odot) - 4.3$. We find good agreement between our best-fit model and recent observations made by JWST at high redshift. The weak evolution of the MZR at $z > 5$ contrasts with the evolution at $z \lesssim 3$, where increasing normalization of the MZR with decreasing redshift is observed and predicted by most models. The FIRE-2 simulations predict increasing scatter in the gas-phase MZR with decreasing stellar mass, in qualitative agreement with some observations. △ Less

Submitted 23 May, 2024; v1 submitted 13 March, 2024; originally announced March 2024.

Comments: 13 pages, 6 figures

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!

Inflow and outflow properties, not total gas fractions, drive the evolution of the mass-metallicity relation

Authors: Luigi Bassini, Robert Feldmann, Jindra Gensior, Claude-André Faucher-Giguère, Elia Cenci, Jorge Moreno, Mauro Bernardini, Lichen Liang

Abstract: Observations show a tight correlation between the stellar mass of galaxies and their gas-phase metallicity (MZR). This relation evolves with redshift, with higher-redshift galaxies being characterized by lower metallicities. Understanding the physical origin of the slope and redshift evolution of the MZR may provide important insight into the physical processes underpinning it: star formation, fee… ▽ More Observations show a tight correlation between the stellar mass of galaxies and their gas-phase metallicity (MZR). This relation evolves with redshift, with higher-redshift galaxies being characterized by lower metallicities. Understanding the physical origin of the slope and redshift evolution of the MZR may provide important insight into the physical processes underpinning it: star formation, feedback, and cosmological inflows. While theoretical models ascribe the shape of the MZR to the lower efficiency of galactic outflows in more massive galaxies, what drives its evolution remains an open question. In this letter, we analyze how the MZR evolves over $z=0-3$, combining results from the FIREbox cosmological volume simulation with analytical models. Contrary to a frequent assertion in the literature, we find that the evolution of the gas fraction does not contribute significantly to the redshift evolution of the MZR. Instead, we show that the latter is driven by the redshift-dependence of the inflow metallicity, outflow metallicity, and mass loading factor, whose relative importance depends on stellar mass. These findings also suggest that the evolution of the MZR is not explained by galaxies moving along a fixed surface in the space spanned by stellar mass, gas phase metallicity, and star formation rate. △ Less

Submitted 24 January, 2024; originally announced January 2024.

Comments: 6 pages, 4 figures. Submitted for publication on MNRAS

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

Two Distinct Classes of Quiescent Galaxies at Cosmic Noon Revealed by JWST PRIMER and UNCOVER

Authors: Sam E. Cutler, Katherine E. Whitaker, John R. Weaver, Bingjie Wang, Richard Pan, Rachel Bezanson, Lukas J. Furtak, Ivo Labbe, Joel Leja, Sedona H. Price, Yingjie Cheng, Maike Clausen, Fergus Cullen, Pratika Dayal, Anna de Graaff, Mark Dickinson, James S. Dunlop, Robert Feldmann, Marijn Franx, Mauro Giavalisco, Karl Glazebrook, Jenny E. Greene, Norman A. Grogin, Garth Illingworth, Anton M. Koekemoer , et al. (9 additional authors not shown)

Abstract: We present a measurement of the low-mass quiescent size-mass relation at Cosmic Noon (1<z<3) from the JWST PRIMER and UNCOVER treasury surveys, which highlights two distinct classes of quiescent galaxies. While the massive population is well studied at these redshifts, the low-mass end has been previously under-explored due to a lack of observing facilities with sufficient sensitivity and spatial… ▽ More We present a measurement of the low-mass quiescent size-mass relation at Cosmic Noon (1<z<3) from the JWST PRIMER and UNCOVER treasury surveys, which highlights two distinct classes of quiescent galaxies. While the massive population is well studied at these redshifts, the low-mass end has been previously under-explored due to a lack of observing facilities with sufficient sensitivity and spatial resolution. We select a conservative sample of low-mass quiescent galaxy candidates using rest-frame UVJ colors and specific star formation rate criteria and measure galaxy morphology in both rest-frame UV/optical wavelengths (F150W) and rest-frame near-infrared (F444W). We confirm an unambiguous flattening of the low-mass quiescent size-mass relation, which results from the separation of the quiescent galaxy sample into two distinct populations at $\log(M_\star/M_\odot)\sim10.3$: low-mass quiescent galaxies that are notably younger and have disky structures, and massive galaxies consistent with spheroidal morphologies and older median stellar ages. These separate populations imply mass quenching dominates at the massive end while other mechanisms, such as environmental or feedback-driven quenching, form the low-mass end. This stellar mass dependent slope of the quiescent size-mass relation could also indicate a shift from size growth due to star formation (low masses) to growth via mergers (massive galaxies). The transition mass between these two populations also corresponds with other dramatic changes and characteristic masses in several galaxy evolution scaling relations (e.g. star-formation efficiency, dust obscuration, and stellar-halo mass ratios), further highlighting the stark dichotomy between low-mass and massive galaxy formation. △ Less

Submitted 23 April, 2024; v1 submitted 22 December, 2023; originally announced December 2023.

Comments: 16 pages, 6 figures, 1 table. Submitted to ApJL. Revised 2024 April 18

Starburst-induced gas-stars kinematic misalignment

Authors: Elia Cenci, Robert Feldmann, Jindra Gensior, James S. Bullock, Jorge Moreno, Luigi Bassini, Mauro Bernardini

Abstract: A kinematic misalignment of the stellar and gas components is a phenomenon observed in a significant fraction of galaxies. However, the underlying physical mechanisms are not well understood. A commonly proposed scenario for the formation of a misaligned component requires any pre-existing gas disc to be removed, via fly-bys or ejective feedback from an active galactic nucleus. In this Letter, we… ▽ More A kinematic misalignment of the stellar and gas components is a phenomenon observed in a significant fraction of galaxies. However, the underlying physical mechanisms are not well understood. A commonly proposed scenario for the formation of a misaligned component requires any pre-existing gas disc to be removed, via fly-bys or ejective feedback from an active galactic nucleus. In this Letter, we study the evolution of a Milky Way mass galaxy in the FIREbox cosmological volume that displays a thin, counter-rotating gas disc with respect to its stellar component at low redshift. In contrast to scenarios involving gas ejection, we find that pre-existing gas is mainly removed via the conversion into stars in a central starburst, triggered by a merging satellite galaxy. The newly-accreted, counter-rotating gas eventually settles into a kinematically misaligned disc. About 4.4 (8 out of 182) of FIREbox galaxies with stellar masses larger than 5e9 Msun at z=0 exhibit gas-star kinematic misalignment. In all cases, we identify central starburst-driven depletion as the main reason for the removal of the pre-existing co-rotating gas component, with no need for feedback from, e.g., a central active black hole. However, during the starburst, the gas is funneled towards the central regions, likely enhancing black hole activity. By comparing the fraction of misaligned discs between FIREbox and other simulations and observations, we conclude that this channel might have a non-negligible role in inducing kinematic misalignment in galaxies. △ Less

Submitted 12 December, 2023; originally announced December 2023.

Comments: 9 pages; 3 figures; submitted to ApJ Letters

The HI covering fraction of Lyman Limit Systems in FIRE haloes

Authors: Lucas Tortora, Robert Feldmann, Mauro Bernardini, Claude-André Faucher-Giguère

Abstract: Atomic hydrogen (HI) serves a crucial role in connecting galactic-scale properties such as star formation with the large-scale structure of the Universe. While recent numerical simulations have successfully matched the observed covering fraction of HI near Lyman Break Galaxies (LBGs) and in the foreground of luminous quasars at redshifts $z \lesssim 3$, the low-mass end remains as-of-yet unexplore… ▽ More Atomic hydrogen (HI) serves a crucial role in connecting galactic-scale properties such as star formation with the large-scale structure of the Universe. While recent numerical simulations have successfully matched the observed covering fraction of HI near Lyman Break Galaxies (LBGs) and in the foreground of luminous quasars at redshifts $z \lesssim 3$, the low-mass end remains as-of-yet unexplored in observational and computational surveys. We employ a cosmological, hydrodynamical simulation (FIREbox) supplemented with zoom-in simulations (MassiveFIRE) from the Feedback In Realistic Environments (FIRE) project to investigate the HI covering fraction of Lyman Limit Systems ($N_{\mathrm{HI}} \gtrsim 10^{17.2}$ cm$^{-2}$) across a wide range of redshifts ($z=0-6$) and halo masses ($10^8-10^{13} M_{\odot}$ at $z=0$, $10^8-10^{11} M_{\odot}$ at $z=6$) in the absence of feedback from active galactic nuclei (AGN). We find that the covering fraction inside haloes exhibits a strong increase with redshift, with only a weak dependence on halo mass for higher-mass haloes. For massive haloes ($M_{\mathrm{vir}} \sim 10^{11}-10^{12} M_{\odot}$), the radial profiles showcase scale-invariance and remain independent of mass. The radial dependence is well-captured by a fitting function. The covering fractions in our simulations are in good agreement with measurements of the covering fraction in LBGs. Our comprehensive analysis unveils a complex dependence with redshift and halo mass for haloes with $M_{\mathrm{vir}} \lesssim 10^{10} M_{\odot}$ that future observations aim to constrain, providing key insights into the physics of structure formation and gas assembly. △ Less

Submitted 29 November, 2023; originally announced November 2023.

Great Balls of FIRE III: Modeling Black Hole Mergers from Massive Star Clusters in Simulations of Galaxies

Authors: Tristan Bruel, Carl L. Rodriguez, Astrid Lamberts, Michael Y. Grudic, Zachary Hafen, Robert Feldmann

Abstract: After the nearly hundred gravitational-wave detections reported by the LIGO-Virgo-KAGRA Collaboration, the question of the cosmological origin of merging binary black holes (BBHs) remains open. The two main formation channels generally considered are from isolated field binaries or via dynamical assembly in dense star clusters. Here, we focus on understanding the dynamical formation of merging BBH… ▽ More After the nearly hundred gravitational-wave detections reported by the LIGO-Virgo-KAGRA Collaboration, the question of the cosmological origin of merging binary black holes (BBHs) remains open. The two main formation channels generally considered are from isolated field binaries or via dynamical assembly in dense star clusters. Here, we focus on understanding the dynamical formation of merging BBHs within massive clusters in galaxies of different masses. To this end, we apply a new framework to consistently model the formation and evolution of massive star clusters in zoom-in cosmological simulations of galaxies. Each simulation, taken from the FIRE project, provides a realistic star formation environment with a unique star formation history and hosts realistic giant molecular clouds that constitute the birthplace of star clusters. Combined with the code for star cluster evolution CMC, we are able to produce populations of dynamically formed merging BBHs across cosmic time in different environments. As the most massive star clusters preferentially form in dense massive clouds of gas, we find that, despite their low metallicities favourable to the creation of black holes, low-mass galaxies contain few massive clusters and therefore have a limited contribution to the global production of dynamically formed merging BBHs. Furthermore, we find that massive clusters can host hierarchical BBH mergers with clear identifiable physical properties. Looking at the evolution of the BBH merger rate in different galaxies, we find strong correlations between BBH mergers and the most extreme episodes of star formation. Finally, we discuss the implications for future LIGO-Virgo-KAGRA gravitational wave observations. △ Less

Submitted 24 November, 2023; originally announced November 2023.

Comments: 14 pages, 9 figures, 3 tables

Journal ref: A&A 686, A106 (2024)

Galaxies Going Bananas: Inferring the 3D Geometry of High-Redshift Galaxies with JWST-CEERS

Authors: Viraj Pandya, Haowen Zhang, Marc Huertas-Company, Kartheik G. Iyer, Elizabeth McGrath, Guillermo Barro, Steven L. Finkelstein, Martin Kuemmel, William G. Hartley, Henry C. Ferguson, Jeyhan S. Kartaltepe, Joel Primack, Avishai Dekel, Sandra M. Faber, David C. Koo, Greg L. Bryan, Rachel S. Somerville, Ricardo O. Amorin, Pablo Arrabal Haro, Micaela B. Bagley, Eric F. Bell, Emmanuel Bertin, Luca Costantin, Romeel Dave, Mark Dickinson , et al. (31 additional authors not shown)

Abstract: The 3D geometry of high-redshift galaxies remains poorly understood. We build a differentiable Bayesian model and use Hamiltonian Monte Carlo to efficiently and robustly infer the 3D shapes of star-forming galaxies in JWST-CEERS observations with $\log M_*/M_{\odot}=9.0-10.5$ at $z=0.5-8.0$. We reproduce previous results from HST-CANDELS in a fraction of the computing time and constrain the mean e… ▽ More The 3D geometry of high-redshift galaxies remains poorly understood. We build a differentiable Bayesian model and use Hamiltonian Monte Carlo to efficiently and robustly infer the 3D shapes of star-forming galaxies in JWST-CEERS observations with $\log M_*/M_{\odot}=9.0-10.5$ at $z=0.5-8.0$. We reproduce previous results from HST-CANDELS in a fraction of the computing time and constrain the mean ellipticity, triaxiality, size and covariances with samples as small as $\sim50$ galaxies. We find high 3D ellipticities for all mass-redshift bins suggesting oblate (disky) or prolate (elongated) geometries. We break that degeneracy by constraining the mean triaxiality to be $\sim1$ for $\log M_*/M_{\odot}=9.0-9.5$ dwarfs at $z>1$ (favoring the prolate scenario), with significantly lower triaxialities for higher masses and lower redshifts indicating the emergence of disks. The prolate population traces out a ``banana'' in the projected $b/a-\log a$ diagram with an excess of low $b/a$, large $\log a$ galaxies. The dwarf prolate fraction rises from $\sim25\%$ at $z=0.5-1.0$ to $\sim50-80\%$ at $z=3-8$. If these are disks, they cannot be axisymmetric but instead must be unusually oval (triaxial) unlike local circular disks. We simultaneously constrain the 3D size-mass relation and its dependence on 3D geometry. High-probability prolate and oblate candidates show remarkably similar Sérsic indices ($n\sim1$), non-parametric morphological properties and specific star formation rates. Both tend to be visually classified as disks or irregular but edge-on oblate candidates show more dust attenuation. We discuss selection effects, follow-up prospects and theoretical implications. △ Less

Submitted 15 January, 2024; v1 submitted 23 October, 2023; originally announced October 2023.

Comments: Accepted version to appear in ApJ, main body is 36 pages of which ~half are full-page figures

HI discs of L$_{\ast}$ galaxies as probes of the baryonic physics of galaxy evolution

Authors: Jindra Gensior, Robert Feldmann, Marta Reina-Campos, Sebastian Trujillo-Gomez, Lucio Mayer, Benjamin W. Keller, Andrew Wetzel, J. M. Diederik Kruijssen, Philip F. Hopkins, Jorge Moreno

Abstract: Understanding what shapes the cold gas component of galaxies, which both provides the fuel for star formation and is strongly affected by the subsequent stellar feedback, is a crucial step towards a better understanding of galaxy evolution. Here, we analyse the HI properties of a sample of 46 Milky Way halo-mass galaxies, drawn from cosmological simulations (EMP-Pathfinder and FIREbox). This set o… ▽ More Understanding what shapes the cold gas component of galaxies, which both provides the fuel for star formation and is strongly affected by the subsequent stellar feedback, is a crucial step towards a better understanding of galaxy evolution. Here, we analyse the HI properties of a sample of 46 Milky Way halo-mass galaxies, drawn from cosmological simulations (EMP-Pathfinder and FIREbox). This set of simulations comprises galaxies evolved self-consistently across cosmic time with different baryonic sub-grid physics: three different star formation models [constant star formation efficiency (SFE) with different star formation eligibility criteria, and an environmentally-dependent, turbulence-based SFE] and two different feedback prescriptions, where only one sub-sample includes early stellar feedback. We use these simulations to assess the impact of different baryonic physics on the HI content of galaxies. We find that the galaxy-wide HI properties agree with each other and with observations. However, differences appear for small-scale properties. The thin HI discs observed in the local Universe are only reproduced with a turbulence-dependent SFE and/or early stellar feedback. Furthermore, we find that the morphology of HI discs is particularly sensitive to the different physics models: galaxies simulated with a turbulence-based SFE have discs that are smoother and more rotationally symmetric, compared to those simulated with a constant SFE; galaxies simulated with early stellar feedback have more regular discs than supernova-feedback-only galaxies. We find that the rotational asymmetry of the HI discs depends most strongly on the underlying physics model, making this a promising observable for understanding the physics responsible for shaping the interstellar medium of galaxies. △ Less

Submitted 7 May, 2024; v1 submitted 2 October, 2023; originally announced October 2023.

Comments: 13 pages, 10 figures + appendices (7 pages, 10 figures); updated to version accepted by MNRAS

The UNCOVER Survey: A First-look HST+JWST Catalog of Galaxy Redshifts and Stellar Population Properties Spanning $0.2 \lesssim z \lesssim 15$

Authors: Bingjie Wang, Joel Leja, Ivo Labbé, Rachel Bezanson, Katherine E. Whitaker, Gabriel Brammer, Lukas J. Furtak, John R. Weaver, Sedona H. Price, Adi Zitrin, Hakim Atek, Dan Coe, Sam E. Cutler, Pratika Dayal, Pieter van Dokkum, Robert Feldmann, Danilo Marchesini, Marijn Franx, Natascha Förster Schreiber, Seiji Fujimoto, Marla Geha, Karl Glazebrook, Anna de Graaff, Jenny E. Greene, Stéphanie Juneau , et al. (19 additional authors not shown)

Abstract: The recent UNCOVER survey with the James Webb Space Telescope (JWST) exploits the nearby cluster Abell 2744 to create the deepest view of our universe to date by leveraging strong gravitational lensing. In this work, we perform photometric fitting of more than 50,000 robustly detected sources out to $z \sim 15$. We show the redshift evolution of stellar ages, star formation rates, and rest-frame c… ▽ More The recent UNCOVER survey with the James Webb Space Telescope (JWST) exploits the nearby cluster Abell 2744 to create the deepest view of our universe to date by leveraging strong gravitational lensing. In this work, we perform photometric fitting of more than 50,000 robustly detected sources out to $z \sim 15$. We show the redshift evolution of stellar ages, star formation rates, and rest-frame colors across the full range of $0.2 \lesssim z \lesssim 15$. The galaxy properties are inferred using the Prospector Bayesian inference framework using informative Prospector-$β$ priors on masses and star formation histories to produce joint redshift and stellar population posteriors, and additionally lensing magnification is performed on-the-fly to ensure consistency with the scale-dependent priors. We show that this approach produces excellent photometric redshifts with $σ_{\rm NMAD} \sim 0.03$, of a similar quality to the established photometric redshift code EAzY. In line with the open-source scientific objective of the Treasury survey, we publicly release the stellar population catalog with this paper, derived from the photometric catalog adapting aperture sizes based on source profiles. This release includes posterior moments, maximum-likelihood spectra, star-formation histories, and full posterior distributions, offering a rich data set to explore the processes governing galaxy formation and evolution over a parameter space now accessible by JWST. △ Less

Submitted 16 April, 2024; v1 submitted 2 October, 2023; originally announced October 2023.

Comments: Corrected typos: Eq.1 should've been (1-kappa)^2, and the lens maps are normalized to D_ds/D_s=1. These errors were only in the writing; no data products or results were affected. The SPS catalogs are accessible via the UNCOVER survey webpage: https://jwst-uncover.github.io/DR2.html#SPSCatalogs, with a copy deposited to Zenodo: https://doi.org/10.5281/zenodo.8401181

Journal ref: The Astrophysical Journal Supplement Series, 270, 12 (2024)

Starbursts driven by central gas compaction

Authors: Elia Cenci, Robert Feldmann, Jindra Gensior, Jorge Moreno, Luigi Bassini, Mauro Bernardini

Abstract: Starburst (SB) galaxies are a rare population of galaxies with star formation rates (SFRs) greatly exceeding those of the majority of star-forming galaxies with similar stellar mass. It is unclear whether these bursts are the result of either especially large gas reservoirs or enhanced efficiencies in converting gas into stars. Tidal torques resulting from gas-rich galaxy mergers are known to enha… ▽ More Starburst (SB) galaxies are a rare population of galaxies with star formation rates (SFRs) greatly exceeding those of the majority of star-forming galaxies with similar stellar mass. It is unclear whether these bursts are the result of either especially large gas reservoirs or enhanced efficiencies in converting gas into stars. Tidal torques resulting from gas-rich galaxy mergers are known to enhance the SFR by funneling gas towards the centre. However, recent theoretical works show that mergers do not always trigger a SB and not all SB galaxies are interacting systems, raising the question of what drives a SB. We analyse a large sample of SB galaxies and a mass- and redshift-matched sample of control galaxies, drawn from the FIREbox cosmological volume at z=0-1. We find that SB galaxies have both larger molecular gas fractions and shorter molecular depletion times than control galaxies, but similar total gas masses. Control galaxies evolve towards the SB regime by gas compaction in their central regions, over timescales of about 70 Myr, accompanied by an increase in the fraction of ultra-dense and molecular gas. The driving mechanism behind the SB varies depending on the mass of the galaxy. Massive (Mstar > 1e10 Msun) galaxies undergoing intense, long-lasting SBs are mostly driven by galaxy interactions. Conversely, SBs in non-interacting galaxies are often triggered by a global gravitational instability, that can result in a breathing mode in low-mass galaxies. △ Less

Submitted 16 September, 2023; originally announced September 2023.

Comments: 21 pages, 18 figures, submitted to MNRAS

DUALZ: Deep UNCOVER-ALMA Legacy High-Z Survey

Authors: Seiji Fujimoto, Rachel Bezanson, Ivo Labbe, Gabriel Brammer, Sedona H. Price, Bingjie Wang, John R. Weaver, Yoshinobu Fudamoto, Pascal A. Oesch, Christina C. Williams, Pratika Dayal, Robert Feldmann, Jenny E. Greene, Joel Leja, Katherine E. Whitaker, Adi Zitrin, Sam E. Cutler, Lukas J. Furtak, Richard Pan, Iryna Chemerynska, Vasily Kokorev, Tim B. Miller, Hakim Atek, Pieter van Dokkum, Stephanie Juneau , et al. (7 additional authors not shown)

Abstract: We present the survey design and initial results of the ALMA Cycle 9 program of DUALZ, which aims to establish a joint ALMA and JWST public legacy field targeting the massive galaxy cluster Abell 2744. DUALZ features a contiguous $4'\times6'$ ALMA 30-GHz-wide mosaic in Band 6, covering areas of $μ>2$ down to a sensitivity of $σ=32.7~μ$Jy. Through a blind search, we identified 69 dust continuum sou… ▽ More We present the survey design and initial results of the ALMA Cycle 9 program of DUALZ, which aims to establish a joint ALMA and JWST public legacy field targeting the massive galaxy cluster Abell 2744. DUALZ features a contiguous $4'\times6'$ ALMA 30-GHz-wide mosaic in Band 6, covering areas of $μ>2$ down to a sensitivity of $σ=32.7~μ$Jy. Through a blind search, we identified 69 dust continuum sources at S/N $\gtrsim5.0$ with median redshift and intrinsic 1.2-mm flux of $z=2.30$ and $S_{\rm 1.2mm}^{\rm int}=0.24$~mJy. Of these, 27 have been spectroscopically confirmed, leveraged by the latest NIRSpec observations, while photometric redshift estimates are constrained by the comprehensive HST, NIRCam, and ALMA data for the remaining sources. With priors, we further identify a [CII]158 $μ$m line emitter at $z=6.3254\pm0.0004$, confirmed by the latest NIRSpec spectroscopy. The NIRCam counterparts of the 1.2-mm continuum exhibit undisturbed morphologies, denoted either by disk or spheroid, implying the triggers for the faint mm emission are less catastrophic than mergers. We have identified 8 HST-dark galaxies (F150W$>$27mag, F150W$-$F444W$>$2.3) and 2 JWST-dark (F444W$>$30mag) galaxy candidates among the ALMA continuum sources. The former includes face-on disk galaxies, hinting that substantial dust obscuration does not always result from inclination. We also detect a marginal dust emission from an X-ray-detected galaxy at $z_{\rm spec}=10.07$, suggesting an active co-evolution of the central black hole and its host. We assess the infrared luminosity function up to $z\sim10$ and find it consistent with predictions from galaxy formation models. To foster diverse scientific outcomes from the community, we publicly release reduced ALMA mosaic maps, cubes, and the source catalog. △ Less

Submitted 16 September, 2023; v1 submitted 14 September, 2023; originally announced September 2023.

Comments: 33 pages, 16 figures, and 5 tables. Submitted to ApJS. The ALMA products are fully available from here: https://jwst-uncover.github.io/DR2.html#DUALZ

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

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)

Most of the photons that reionized the Universe came from dwarf galaxies

Authors: Hakim Atek, Ivo Labbé, Lukas J. Furtak, Iryna Chemerynska, Seiji Fujimoto, David J. Setton, Tim B. Miller, Pascal Oesch, Rachel Bezanson, Sedona H. Price, Pratika Dayal, Adi Zitrin, Vasily Kokorev, John R. Weaver, Gabriel Brammer, Pieter van Dokkum, Christina C. Williams, Sam E. Cutler, Robert Feldmann, Yoshinobu Fudamoto, Jenny E. Greene, Joel Leja, Michael V. Maseda, Adam Muzzin, Richard Pan , et al. (8 additional authors not shown)

Abstract: The identification of sources driving cosmic reionization, a major phase transition from neutral Hydrogen to ionized plasma around 600-800 Myr after the Big Bang (Dayal et al. 2018, Mason et al. 2019, Robertson et al. 2022), has been a matter of intense debate (Robertson et al. 2022). Some models suggest that high ionizing emissivity and escape fractions ($f_{\rm esc}$) from quasars support their… ▽ More The identification of sources driving cosmic reionization, a major phase transition from neutral Hydrogen to ionized plasma around 600-800 Myr after the Big Bang (Dayal et al. 2018, Mason et al. 2019, Robertson et al. 2022), has been a matter of intense debate (Robertson et al. 2022). Some models suggest that high ionizing emissivity and escape fractions ($f_{\rm esc}$) from quasars support their role in driving cosmic reionization (Madau & Haardt 2015, Mitra et al. 2018). Others propose that the high $f_{\rm esc}$ values from bright galaxies generates sufficient ionizing radiation to drive this process (Naidu et al. 2020). Finally, a few studies suggest that the number density of faint galaxies, when combined with a stellar-mass-dependent model of ionizing efficiency and $f_{\rm esc}$, can effectively dominate cosmic reionization (Finkelstein et al. 2019, Dayal et al. 2020). However, so far, low-mass galaxies have eluded comprehensive spectroscopic studies owing to their extreme faintness. Here we report an analysis of eight ultra-faint galaxies (in a very small field) during the epoch of reionization with absolute magnitudes between $M_{\rm UV}$ $\sim -17$ to $-15$ mag (down to 0.005 $L^{\star}$. We find that faint galaxies during the Universe's first billion years produce ionizing photons with log($ξ_{\mathrm{ion}}$/ Hz erg$^{-1}$) =$25.80\pm 0.14$, a factor of 4 higher than commonly assumed values (Robertson et al. 2015). If this field is representative of the large scale distribution of faint galaxies, the rate of ionizing photons exceeds that needed for reionization, even for escape fractions of order five per cent. △ Less

Submitted 30 April, 2024; v1 submitted 16 August, 2023; originally announced August 2023.

Comments: 29 pages, 7 figures, 2 tables. Published in Nature

Journal ref: Nature Volume 626, 2024, 975-978

Symmetry-Projected Nuclear-Electronic Hartree-Fock: Eliminating Rotational Energy Contamination

Authors: Robin Feldmann, Alberto Baiardi, Markus Reiher

Abstract: We present a symmetry projection technique for enforcing rotational and parity symmetries in nuclear-electronic Hartree-Fock wave functions, which treat electrons and nuclei on equal footing. The molecular Hamiltonian obeys rotational and parity-inversion symmetries, which are, however, broken by expanding in Gaussian basis sets that are fixed in space. We generate a trial wave function with the c… ▽ More We present a symmetry projection technique for enforcing rotational and parity symmetries in nuclear-electronic Hartree-Fock wave functions, which treat electrons and nuclei on equal footing. The molecular Hamiltonian obeys rotational and parity-inversion symmetries, which are, however, broken by expanding in Gaussian basis sets that are fixed in space. We generate a trial wave function with the correct symmetry properties by projecting the wave function onto representations of the three-dimensional rotation group, i.e., the special orthogonal group in three dimensions SO(3). As a consequence, the wave function becomes an eigenfunction of the angular momentum operator which (i) eliminates the contamination of the ground state wave function by highly excited rotational states arising from the broken rotational symmetry, and (ii) enables the targeting of specific rotational states of the molecule. We demonstrate the efficiency of the symmetry projection technique by calculating energies of the low-lying rotational states of the H$_2$ and H$_3^+$ molecules. △ Less

Submitted 2 November, 2023; v1 submitted 15 July, 2023; originally announced July 2023.

Comments: 37 pages, 3 tables

Journal ref: J. Phys. Chem. A 2023, 127, 8943-8954

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

JWST constraints on the UV luminosity density at cosmic dawn: implications for 21-cm cosmology

Authors: Sultan Hassan, Christopher C. Lovell, Piero Madau, Marc Huertas-Company, Rachel S. Somerville, Blakesley Burkhart, Keri L. Dixon, Robert Feldmann, Tjitske K. Starkenburg, John F. Wu, Christian Kragh Jespersen, Joseph D. Gelfand, Ankita Bera

Abstract: An unprecedented array of new observational capabilities are starting to yield key constraints on models of the epoch of first light in the Universe. In this Letter we discuss the implications of the UV radiation background at cosmic dawn inferred by recent JWST observations for radio experiments aimed at detecting the redshifted 21-cm hyperfine transition of diffuse neutral hydrogen. Under the ba… ▽ More An unprecedented array of new observational capabilities are starting to yield key constraints on models of the epoch of first light in the Universe. In this Letter we discuss the implications of the UV radiation background at cosmic dawn inferred by recent JWST observations for radio experiments aimed at detecting the redshifted 21-cm hyperfine transition of diffuse neutral hydrogen. Under the basic assumption that the 21-cm signal is activated by the Ly$α$ photon field produced by metal-poor stellar systems, we show that a detection at the low frequencies of the EDGES and SARAS3 experiments may be expected from a simple extrapolation of the declining UV luminosity density inferred at $z\lesssim 14$ from JWST early galaxy data. Accounting for an early radiation excess above the CMB suggests a shallower or flat evolution to simultaneously reproduce low and high-$z$ current UV luminosity density constraints, which cannot be entirely ruled out, given the large uncertainties from cosmic variance and the faint-end slope of the galaxy luminosity function at cosmic dawn. Our findings raise the intriguing possibility that a high star formation efficiency at early times may trigger the onset of intense Ly$α$ emission at redshift $z\lesssim 20$ and produce a cosmic 21-cm absorption signal 200 Myr after the Big Bang. △ Less

Submitted 11 October, 2023; v1 submitted 4 May, 2023; originally announced May 2023.

Comments: 8 pages, 3 figures, accepted to ApJL

$\rm [C_{II}]$ 158 $\rm μm$ emission as an indicator of galaxy star formation rate

Authors: Lichen Liang, Robert Feldmann, Norman Murray, Desika Narayanan, Christopher C. Hayward, Daniel Anglés-Alcázar, Luigi Bassini, Alexander J. Richings, Claude-André Faucher-Giguère, Dongwoo T. Chung, Jennifer Y. H. Chan, Doǧa Tolgay, Onur Çatmabacak, Dušan Kereš, Philip F. Hopkins

Abstract: Observations of local star-forming galaxies (SFGs) show a tight correlation between their singly ionized carbon line luminosity ($L_{\rm [C_{II}]}$) and star formation rate (SFR), suggesting that $L_{\rm [C_{II}]}$ may be a useful SFR tracer for galaxies. Some other galaxy populations, however, are found to have lower $L_{\rm [C_{II}]}{}/{}\rm SFR$ than the local SFGs, including the infrared-lumin… ▽ More Observations of local star-forming galaxies (SFGs) show a tight correlation between their singly ionized carbon line luminosity ($L_{\rm [C_{II}]}$) and star formation rate (SFR), suggesting that $L_{\rm [C_{II}]}$ may be a useful SFR tracer for galaxies. Some other galaxy populations, however, are found to have lower $L_{\rm [C_{II}]}{}/{}\rm SFR$ than the local SFGs, including the infrared-luminous, starburst galaxies at low and high redshifts, as well as some moderately star-forming galaxies at the epoch of re-ionization (EoR). The origin of this `$\rm [C_{II}]$ deficit' is unclear. In this work, we study the $L_{\rm [C_{II}]}$-SFR relation of galaxies using a sample of $z=0-8$ galaxies with $M_*\approx10^7-5\times10^{11}\,M_\odot$ extracted from cosmological volume and zoom-in simulations from the Feedback in Realistic Environments (FIRE) project. We find a simple analytic expression for $L_{\rm [C_{II}]}$/SFR of galaxies in terms of the following parameters: mass fraction of $\rm [C_{II}]$-emitting gas ($f_{\rm [C_{II}]}$), gas metallicity ($Z_{\rm gas}$), gas density ($n_{\rm gas}$) and gas depletion time ($t_{\rm dep}{}={}M_{\rm gas}{}/{}\rm SFR$). We find two distinct physical regimes, where $t_{\rm dep}$ ($Z_{\rm gas}$) is the main driver of the $\rm [C_{II}]$ deficit in $\rm H_2$-rich ($\rm H_2$-poor) galaxies. The observed $\rm [C_{II}]$ deficit of IR-luminous galaxies and early EoR galaxies, corresponding to the two different regimes, is due to short gas depletion time and low gas metallicity, respectively. Our result indicates that $\rm [C_{II}]$ deficit is a common phenomenon of galaxies, and caution needs to be taken when applying a constant $L_{\rm [C_{II}]}$-to-SFR conversion factor derived from local SFGs to estimate cosmic SFR density at high redshifts and interpret data from upcoming $\rm [C_{II}]$ line intensity mapping experiments. △ Less

Submitted 6 December, 2023; v1 submitted 10 January, 2023; originally announced January 2023.

Comments: 44 page, 24 figures, 8 tables and 8 appendices, accepted for publication in MNRAS. Key figures: fig. 16 & 17. Comments are welcome

The UNCOVER Survey: A first-look HST+JWST catalog of 60,000 galaxies near Abell 2744 and beyond

Authors: John R. Weaver, Sam E. Cutler, Richard Pan, Katherine E. Whitaker, Ivo Labbe, Sedona H. Price, Rachel Bezanson, Gabriel Brammer, Danilo Marchesini, Joel Leja, Bingjie Wang, Lukas J. Furtak, Adi Zitrin, Hakim Atek, Dan Coe, Pratika Dayal, Pieter van Dokkum, Robert Feldmann, Natascha Forster Schreiber, Marijn Franx, Seiji Fujimoto, Yoshinobu Fudamoto, Karl Glazebrook, Anna de Graaff, Jenny E. Greene , et al. (19 additional authors not shown)

Abstract: In November 2022, the James Webb Space Telescope (JWST) returned deep near-infrared images of Abell 2744 -- a powerful lensing cluster capable of magnifying distant, incipient galaxies beyond it. Together with the existing Hubble Space Telescope (HST) imaging, this publicly available dataset opens a fundamentally new discovery space to understand the remaining mysteries of the formation and evolut… ▽ More In November 2022, the James Webb Space Telescope (JWST) returned deep near-infrared images of Abell 2744 -- a powerful lensing cluster capable of magnifying distant, incipient galaxies beyond it. Together with the existing Hubble Space Telescope (HST) imaging, this publicly available dataset opens a fundamentally new discovery space to understand the remaining mysteries of the formation and evolution of galaxies across cosmic time. In this work, we detect and measure some 60,000 objects across the 49 arcmin$^2$ JWST footprint down to a $5\,σ$ limiting magnitude of $\sim$30 mag in 0.32" apertures. Photometry is performed using circular apertures on images matched to the point spread function of the reddest NIRCam broad band, F444W, and cleaned of bright cluster galaxies and the related intra-cluster light. To give an impression of the photometric performance, we measure photometric redshifts and achieve a $σ_{\rm NMAD}\approx0.03$ based on known, but relatively small, spectroscopic samples. With this paper, we publicly release our HST and JWST PSF-matched photometric catalog with optimally assigned aperture sizes for easy use, along with single aperture catalogs, photometric redshifts, rest-frame colors, and individual magnification estimates. These catalogs will set the stage for efficient and deep spectroscopic follow-up of some of the first JWST-selected samples in Summer 2023. △ Less

Submitted 2 October, 2023; v1 submitted 6 January, 2023; originally announced January 2023.

Comments: 28 pages, 19 figures, resubmitted to ApJS following significant data product improvements. Comments welcome. Catalogs can be accessed at https://jwst-uncover.github.io/DR2.html#PhotometricCatalogs

UNCOVERing the extended strong lensing structures of Abell 2744 with the deepest JWST imaging

Authors: Lukas J. Furtak, Adi Zitrin, John R. Weaver, Hakim Atek, Rachel Bezanson, Ivo Labbe, Katherine E. Whitaker, Joel Leja, Sedona H. Price, Gabriel B. Brammer, Bingjie Wang, Danilo Marchesini, Richard Pan, Pratika Dayal, Pieter van Dokkum, Robert Feldmann, Seiji Fujimoto, Marijn Franx, Gourav Khullar, Erica J. Nelson, Lamiya A. Mowla

Abstract: We present a new parametric lens model for the massive galaxy cluster Abell~2744 based on the new ultra-deep JWST imaging taken in the framework of the UNCOVER program. These observations constitute the deepest JWST images of a lensing cluster to date, adding to the existing deep Hubble Space Telescope (HST) images and the recent JWST ERS and DDT data taken for this field. The wide field-of-view o… ▽ More We present a new parametric lens model for the massive galaxy cluster Abell~2744 based on the new ultra-deep JWST imaging taken in the framework of the UNCOVER program. These observations constitute the deepest JWST images of a lensing cluster to date, adding to the existing deep Hubble Space Telescope (HST) images and the recent JWST ERS and DDT data taken for this field. The wide field-of-view of UNCOVER ($\sim45$ arcmin$^2$) extends beyond the cluster's well-studied central core and reveals a spectacular wealth of prominent lensed features around two massive cluster sub-structures in the north and north-west, where no multiple images were previously known. The 75 newly uncovered multiple images and candidates of 16 sources allow us, for the first time, to constrain the lensing properties and total mass distribution around these extended cluster structures using strong lensing (SL). Our model yields an effective Einstein radius of $θ_{E,\mathrm{main}}\simeq23''$ for the main cluster core (for $z_{\mathrm{s}}=2$), enclosing a mass of $M(θ<θ_{E,\mathrm{main}})\simeq7.7\times10^{13}$ M$_{\odot}$, and $θ_{E,\mathrm{NW}}\simeq13''$ for the newly discovered north-western SL structure enclosing $M(θ<θ_{E,\mathrm{NW}})\simeq2.2\times10^{13}$ M$_{\odot}$. The northern clump is somewhat less massive with $θ_{E,\mathrm{N}}\simeq7''$ enclosing $M(θ<θ_{E,\mathrm{N}})\simeq8\times10^{12}$ M$_{\odot}$. We find the northern sub-structures of Abell~2744 to broadly agree with the findings from weak lensing (WL) and align with the filamentary structure found by these previous studies. Our model in particular reveals a large area of high magnifications between the various cluster structures, which will be paramount for lensed galaxy studies in the UNCOVER field. The model is made publicly available to accompany the first UNCOVER data release. △ Less

Submitted 28 May, 2023; v1 submitted 8 December, 2022; originally announced December 2022.

Comments: Accepted for publication in MNRAS. Updated to match the published version

The JWST UNCOVER Treasury survey: Ultradeep NIRSpec and NIRCam ObserVations before the Epoch of Reionization

Authors: Rachel Bezanson, Ivo Labbe, Katherine E. Whitaker, Joel Leja, Sedona H. Price, Marijn Franx, Gabe Brammer, Danilo Marchesini, Adi Zitrin, Bingjie Wang, John R. Weaver, Lukas J. Furtak, Hakim Atek, Dan Coe, Sam E. Cutler, Pratika Dayal, Pieter van Dokkum, Robert Feldmann, Natascha Forster Schreiber, Seiji Fujimoto, Marla Geha, Karl Glazebrook, Anna de Graaff, Jenny E. Greene, Stephanie Juneau , et al. (18 additional authors not shown)

Abstract: In this paper we describe the survey design for the Ultradeep NIRSpec and NIRCam ObserVations before the Epoch of Reionization (UNCOVER) Cycle 1 JWST Treasury program, which executed its early imaging component in November 2022. The UNCOVER survey includes ultradeep ($\sim29-30\mathrm{AB}$) imaging of $\sim$45 arcmin$^2$ on and around the well-studied Abell 2744 galaxy cluster at $z=0.308$ and wil… ▽ More In this paper we describe the survey design for the Ultradeep NIRSpec and NIRCam ObserVations before the Epoch of Reionization (UNCOVER) Cycle 1 JWST Treasury program, which executed its early imaging component in November 2022. The UNCOVER survey includes ultradeep ($\sim29-30\mathrm{AB}$) imaging of $\sim$45 arcmin$^2$ on and around the well-studied Abell 2744 galaxy cluster at $z=0.308$ and will follow-up ${\sim}500$ galaxies with extremely deep low-resolution spectroscopy with the NIRSpec/PRISM during the summer of 2023. We describe the science goals, survey design, target selection, and planned data releases. We also present and characterize the depths of the first NIRCam imaging mosaic, highlighting previously unparalleled resolved and ultradeep 2-4 micron imaging of known objects in the field. The UNCOVER primary NIRCam mosaic spans 28.8 arcmin$^2$ in seven filters (F115W, F150W, F200W, F277W, F356W, F410M, F444W) and 16.8 arcmin$^2$ in our NIRISS parallel (F115W, F150W, F200W, F356W, and F444W). To maximize early community use of the Treasury data set, we publicly release full reduced mosaics of public JWST imaging including 45 arcmin$^2$ NIRCam and 17 arcmin$^2$ NIRISS mosaics on and around the Abell 2744 cluster, including the Hubble Frontier Field primary and parallel footprints. △ Less

Submitted 9 December, 2022; v1 submitted 7 December, 2022; originally announced December 2022.

Comments: 17 pages, 8 figures, 4 tables, submitted to ApJ, comments welcome (v2 with full author list in metadata)

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

The inefficiency of stellar feedback in driving galactic outflows in massive galaxies at high redshift

Authors: L. Bassini, R. Feldmann, J. Gensior, C. C. Hayward, C. -A. Faucher-Giguère, E. Cenci, L. Liang, M. Bernardini

Abstract: Recent observations indicate that galactic outflows are ubiquitous in high redshift galaxies, including normal star forming galaxies, quasar hosts, and dusty star forming galaxies (DSFGs). However, the impact of outflows on the evolution of their hosts is still an open question. Here, we analyse the star formation histories (SFH) and galactic outflow properties of galaxies in massive haloes (… ▽ More Recent observations indicate that galactic outflows are ubiquitous in high redshift galaxies, including normal star forming galaxies, quasar hosts, and dusty star forming galaxies (DSFGs). However, the impact of outflows on the evolution of their hosts is still an open question. Here, we analyse the star formation histories (SFH) and galactic outflow properties of galaxies in massive haloes ($10^{12}M_{\odot}<M_{\rm vir} <5\times 10^{12}M_{\odot}$) at $z\gtrsim5.5$ in three zoom-in cosmological simulations from the MassiveFIRE suite, as part of the Feedback In Realistic Environments (FIRE) project. The simulations were run with the FIRE-2 model, which does not include feedback from active galactic nuclei (AGN). The simulated galaxies resemble $z>4$ DSFGs, with SFRs of $\sim 1000\ M_{\odot}\rm yr^{-1}$ and molecular gas masses of $M_{\rm mol}\sim 10^{10}\ M_{\odot}$. However, the simulated galaxies are characterised by higher circular velocities than those observed in high-z DSFGs. The mass loading factors from stellar feedback are of the order of $\sim 0.1$, implying that stellar feedback is inefficient in driving galactic outflows and gas is consumed by star formation on much shorter time-scales than it is expelled from the interstellar medium (ISM). We also find that stellar feedback is highly inefficient in self-regulating star formation in this regime, with an average integrated star formation efficiency (SFE) per dynamical time of $30\%$. Finally, compared to FIRE-2 galaxies hosted in similarly massive haloes at lower redshift, we find lower mass loading factors and higher SFEs in the high redshift sample. We argue that both effects originate from the higher total and gas surface densities that characterise high$-z$ massive systems. △ Less

Submitted 15 November, 2022; originally announced November 2022.

Comments: 17 pages, 12 figures. Comments welcome

Second-order self-consistent field algorithms: from classical to quantum nuclei

Authors: Robin Feldmann, Alberto Baiardi, Markus Reiher

Abstract: This work presents a general framework for deriving exact and approximate Newton self-consistent field (SCF) orbital optimization algorithms by leveraging concepts borrowed from differential geometry. Within this framework, we extend the augmented Roothaan--Hall (ARH) algorithm to unrestricted electronic and nuclear-electronic calculations. We demonstrate that ARH yields an excellent compromise be… ▽ More This work presents a general framework for deriving exact and approximate Newton self-consistent field (SCF) orbital optimization algorithms by leveraging concepts borrowed from differential geometry. Within this framework, we extend the augmented Roothaan--Hall (ARH) algorithm to unrestricted electronic and nuclear-electronic calculations. We demonstrate that ARH yields an excellent compromise between stability and computational cost for SCF problems that are hard to converge with conventional first-order optimization strategies. In the electronic case, we show that ARH overcomes the slow convergence of orbitals in strongly-correlated molecules with the example of several iron-sulfur clusters. For nuclear-electronic calculations, ARH significantly enhances the convergence already for small molecules, as demonstrated for a series of protonated water clusters. △ Less

Submitted 11 January, 2023; v1 submitted 18 October, 2022; originally announced October 2022.

Comments: 40 pages, 4 figures, 3 tables

Journal ref: J. Chem. Theory Comput. 2023, 19, 856-873

Star Formation Suppresion by Tidal Removal of Cold Molecular Gas from an Intermediate-Redshift Massive Post-Starburst Galaxy

Authors: Justin S. Spilker, Katherine A. Suess, David J. Setton, Rachel Bezanson, Robert Feldmann, Jenny E. Greene, Mariska Kriek, Sidney Lower, Desika Narayanan, Margaret Verrico

Abstract: Observations and simulations have demonstrated that star formation in galaxies must be actively suppressed to prevent the formation of over-massive galaxies. Galactic outflows driven by stellar feedback or supermassive black hole accretion are often invoked to regulate the amount of cold molecular gas available for future star formation, but may not be the only relevant quenching processes in all… ▽ More Observations and simulations have demonstrated that star formation in galaxies must be actively suppressed to prevent the formation of over-massive galaxies. Galactic outflows driven by stellar feedback or supermassive black hole accretion are often invoked to regulate the amount of cold molecular gas available for future star formation, but may not be the only relevant quenching processes in all galaxies. We present the discovery of vast molecular tidal features extending up to 64 kpc outside of a massive z=0.646 post-starburst galaxy that recently concluded its primary star-forming episode. The tidal tails contain (1.2 +/- 0.1)x10^10 Msun of molecular gas, 47 +/- 5 % of the total cold gas reservoir of the system. Both the scale and magnitude of the molecular tidal features are unprecedented compared to all known nearby or high-redshift merging systems. We infer that the cold gas was stripped from the host galaxies during the merger, which is most likely responsible for triggering the initial burst phase and the subsequent suppression of star formation. While only a single example, this result shows that galaxy mergers can regulate the cold gas contents in distant galaxies by directly removing a large fraction of the molecular gas fuel, and plausibly suppress star formation directly, a qualitatively different physical mechanism than feedback-driven outflows. △ Less

Submitted 29 August, 2022; originally announced August 2022.

Comments: Published in ApJL, 8 pages, 5 figures

Modeling HI at the field level

Authors: Andrej Obuljen, Marko Simonović, Aurel Schneider, Robert Feldmann

Abstract: We use an analytical forward model based on perturbation theory to predict the neutral hydrogen (HI) overdensity maps at low redshifts. We investigate its performance by comparing it directly at the field level to the simulated HI from the IllustrisTNG simulation TNG300-1 ($L=205\ h^{-1}$ Mpc), in both real and redshift space. We demonstrate that HI is a biased tracer of the underlying matter fiel… ▽ More We use an analytical forward model based on perturbation theory to predict the neutral hydrogen (HI) overdensity maps at low redshifts. We investigate its performance by comparing it directly at the field level to the simulated HI from the IllustrisTNG simulation TNG300-1 ($L=205\ h^{-1}$ Mpc), in both real and redshift space. We demonstrate that HI is a biased tracer of the underlying matter field and find that the cubic bias model describes the simulated HI power spectrum to within 1% up to $k=0.4 \;(0.3) \,h\,{\rm Mpc}^{-1}$ in real (redshift) space at redshifts $z=0,1$. Looking at counts in cells, we find an excellent agreement between the theory and simulations for cells as small as 5 $h^{-1}$ Mpc. These results are in line with expectations from perturbation theory and they imply that a perturbative description of the HI field is sufficiently accurate given the characteristics of upcoming 21cm intensity mapping surveys. Additionally, we study the statistical properties of the model error - the difference between the truth and the model. We show that on large scales this error is nearly Gaussian and that it has a flat power spectrum, with amplitude significantly lower than the standard noise inferred from the HI power spectrum. We explain the origin of this discrepancy, discuss its implications for the HI power spectrum Fisher matrix forecasts and argue that it motivates the HI field-level cosmological inference. On small scales in redshift space we use the difference between the model and the truth as a proxy for the Fingers-of-God effect. This allows us to estimate the nonlinear velocity dispersion of HI and show that it is smaller than for the typical spectroscopic galaxy samples at the same redshift. Finally, we provide a simple prescription based on the perturbative forward model which can be used to efficiently generate accurate HI mock data, in real and redshift space. △ Less

Submitted 30 November, 2023; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: Hi-Fi mocks code available here: https://github.com/andrejobuljen/Hi-Fi_mocks. Updated to match version accepted by PRD. Minor changes in best-fit parameter values. Added appendix C

Realistic HI scale heights of Milky Way-mass galaxies in the FIREbox cosmological volume

Authors: Jindra Gensior, Robert Feldmann, Lucio Mayer, Andrew Wetzel, Philip F. Hopkins, Claude-André Faucher-Giguère

Abstract: Accurately reproducing the thin cold gas discs observed in nearby spiral galaxies has been a long standing issue in cosmological simulations. Here, we present measurements of the radially resolved HI scale height in 22 non-interacting Milky Way-mass galaxies from the FIREbox cosmological volume. We measure the HI scale heights using five different approaches commonly used in the literature: fittin… ▽ More Accurately reproducing the thin cold gas discs observed in nearby spiral galaxies has been a long standing issue in cosmological simulations. Here, we present measurements of the radially resolved HI scale height in 22 non-interacting Milky Way-mass galaxies from the FIREbox cosmological volume. We measure the HI scale heights using five different approaches commonly used in the literature: fitting the vertical volume density distribution with a Gaussian, the distance between maximum and half-maximum of the vertical volume density distribution, a semi-empirical description using the velocity dispersion and the galactic gravitational potential, the analytic assumption of hydrostatic equilibrium, and the distance from the midplane which encloses $\gtrsim$60 per cent of the HI mass. We find median HI scale heights, measured using the vertical volume distribution, that range from ~100 pc in the galactic centres to ~800 pc in the outskirts and are in excellent agreement with recent observational results. We speculate that the presence of a realistic multiphase interstellar medium, including cold gas, and realistic stellar feedback are the drivers behind the realistic HI scale heights. △ Less

Submitted 7 November, 2022; v1 submitted 7 July, 2022; originally announced July 2022.

Comments: 5 pages + appendix, 3 figures, updated to version accepted by MNRAS

The observability of galaxy merger signatures in nearby gas-rich spirals

Authors: Rebecca McElroy, Connor Bottrell, Maan H. Hani, Jorge Moreno, Scott M. Croom, Christopher C. Hayward, Angela Twum, Robert Feldmann, Philip F. Hopkins, Lars Hernquist, Bernd Husemann

Abstract: Galaxy mergers are crucial to understanding galaxy evolution, therefore we must determine their observational signatures to select them from large IFU galaxy samples such as MUSE and SAMI. We employ 24 high-resolution idealised hydrodynamical galaxy merger simulations based on the "Feedback In Realistic Environment" (FIRE-2) model to determine the observability of mergers to various configurations… ▽ More Galaxy mergers are crucial to understanding galaxy evolution, therefore we must determine their observational signatures to select them from large IFU galaxy samples such as MUSE and SAMI. We employ 24 high-resolution idealised hydrodynamical galaxy merger simulations based on the "Feedback In Realistic Environment" (FIRE-2) model to determine the observability of mergers to various configurations and stages using synthetic images and velocity maps. Our mergers cover a range of orbital configurations at fixed 1:2.5 stellar mass ratio for two gas rich spirals at low redshift. Morphological and kinematic asymmetries are computed for synthetic images and velocity maps spanning each interaction. We divide the interaction sequence into three: (1) the pair phase; (2) the merging phase; and (3) the post-coalescence phase. We correctly identify mergers between first pericentre passage and 500 Myr after coalescence using kinematic asymmetry with 66% completeness, depending upon merger phase and the field-of-view of the observation. We detect fewer mergers in the pair phase (40%) and many more in the merging and post-coalescence phases (97%). We find that merger detectability decreases with field-of-view, except in retrograde mergers, where centrally concentrated asymmetric kinematic features enhances their detectability. Using a cut-off derived from a combination of photometric and kinematic asymmetry, we increase these detections to 89% overall, 79% in pairs, and close to 100% in the merging and post-coalescent phases. By using this combined asymmetry cut-off we mitigate some of the effects caused by smaller fields-of-view subtended by massively multiplexed integral field spectroscopy programmes. △ Less

Submitted 15 June, 2022; originally announced June 2022.

Comments: Accepted to MNRAS

FIREbox: Simulating galaxies at high dynamic range in a cosmological volume

Authors: Robert Feldmann, Eliot Quataert, Claude-André Faucher-Giguère, Philip F. Hopkins, Onur Çatmabacak, Dušan Kereš, Luigi Bassini, Mauro Bernardini, James S. Bullock, Elia Cenci, Jindra Gensior, Lichen Liang, Jorge Moreno, Andrew Wetzel

Abstract: We introduce a suite of cosmological volume simulations to study the evolution of galaxies as part of the Feedback in Realistic Environments project. FIREbox, the principal simulation of the present suite, provides a representative sample of galaxies (~1000 galaxies with Mstar > 10^8 Msun at z=0) at a resolution (~20 pc, m_b ~ 6x10^4 Msun) comparable to state-of-the-art galaxy zoom-in simulations.… ▽ More We introduce a suite of cosmological volume simulations to study the evolution of galaxies as part of the Feedback in Realistic Environments project. FIREbox, the principal simulation of the present suite, provides a representative sample of galaxies (~1000 galaxies with Mstar > 10^8 Msun at z=0) at a resolution (~20 pc, m_b ~ 6x10^4 Msun) comparable to state-of-the-art galaxy zoom-in simulations. FIREbox captures the multiphase nature of the interstellar medium in a fully cosmological setting (L=22.1 Mpc) thanks to its exceptionally high dynamic range (~10^6) and the inclusion of multi-channel stellar feedback. Here, we focus on validating the simulation predictions by comparing to observational data. We find that simulated galaxies with Mstar < 10^{10.5-11} Msun have star formation rates, gas masses, and metallicities in broad agreement with observations. These galaxy scaling relations extend to low masses (Mstar ~ 10^7 Msun) and follow a (broken) power-law relationship. Also reproduced are the evolution of the cosmic HI density and the HI column density distribution at z~0-5. At low z, FIREbox predicts a peak in the stellar-mass--halo-mass relation, but also a higher abundance of massive galaxies and a higher cosmic star formation rate density than observed, showing that stellar feedback alone is insufficient to reproduce the properties of massive galaxies at late times. Given its high resolution and sample size, FIREbox offers a baseline prediction of galaxy formation theory in a $Λ$CDM Universe while also highlighting modeling challenges to be addressed in next-generation galaxy simulations. △ Less

Submitted 21 April, 2023; v1 submitted 30 May, 2022; originally announced May 2022.

Comments: 31 pages, 17 figures, 3 tables, accepted for publication in MNRAS

The rotation of planet-hosting stars

Authors: Yves Sibony, Ravit Helled, Robert Feldmann

Abstract: Understanding the distribution of angular momentum during the formation of planetary systems is a key topic in astrophysics. Data from the $\textit{Kepler}$ and $\textit{Gaia}$ missions allow to investigate whether stellar rotation is correlated with the presence of planets around Sun-like stars. Here, we perform a statistical analysis of the rotation period of 493 planet-hosting stars. These are… ▽ More Understanding the distribution of angular momentum during the formation of planetary systems is a key topic in astrophysics. Data from the $\textit{Kepler}$ and $\textit{Gaia}$ missions allow to investigate whether stellar rotation is correlated with the presence of planets around Sun-like stars. Here, we perform a statistical analysis of the rotation period of 493 planet-hosting stars. These are matched to a control sample, without detected planets, with similar effective temperatures, masses, radii, metallicities, and ages. We find that planet-hosting stars rotate on average $1.63 \pm 0.40$ days slower. The difference in rotation is statistically significant both in samples including and not including planets confirmed by radial velocity follow-up observations. We also analyse the dependence of rotation distribution on various stellar and planetary properties. Our results could potentially be explained by planet detection biases depending on the rotation period of their host stars in both RV and transit methods. Alternatively, they could point to a physical link between the existence of planets and stellar rotation, emphasising the need to understand the role of angular momentum in the formation and evolution planetary systems. △ Less

Submitted 4 April, 2022; originally announced April 2022.

Comments: 10 pages, 8 figures, 3 Appendix sections (10 pages, 13 figures). Accepted in MNRAS on April 4 2022

The Compact Structures of Massive $z\sim0.7$ Post-Starburst Galaxies in the SQuIGG$\vec{L}$E Sample

Authors: David J. Setton, Margaret Verrico, Rachel Bezanson, Jenny E. Greene, Katherine A. Suess, Andy D. Goulding, Justin S. Spilker, Mariska Kriek, Robert Feldmann, Desika Narayanan, Khalil Hall-Hooper, Erin Kado-Fong

Abstract: We present structural measurements of 145 spectroscopically selected intermediate-redshift (z$\sim$0.7), massive ($M_\star \sim 10^{11} \ M_\odot$) post-starburst galaxies from the SQuIGG$\vec{L}$E Sample measured using wide-depth Hyper Suprime-Cam i-band imaging. This deep imaging allows us to probe the sizes and structures of these galaxies, which we compare to a control sample of star forming a… ▽ More We present structural measurements of 145 spectroscopically selected intermediate-redshift (z$\sim$0.7), massive ($M_\star \sim 10^{11} \ M_\odot$) post-starburst galaxies from the SQuIGG$\vec{L}$E Sample measured using wide-depth Hyper Suprime-Cam i-band imaging. This deep imaging allows us to probe the sizes and structures of these galaxies, which we compare to a control sample of star forming and quiescent galaxies drawn from the LEGA-C Survey. We find that post-starburst galaxies systematically lie $\sim0.1$ dex below the quiescent mass-size (half-light radius) relation, with a scatter of $\sim0.2$ dex. This finding is bolstered by non-parametric measures, such as the Gini coefficient and the concentration, which also reveal these galaxies to have more compact light profiles than both quiescent and star-forming populations at similar mass and redshift. The sizes of post-starburst galaxies show either negative or no correlation with the time since quenching, such that more recently quenched galaxies are larger or similarly sized. This empirical finding disfavors the formation of post-starburst galaxies via a purely central burst of star formation that simultaneously shrinks the galaxy and shuts off star formation. We show that the central densities of post-starburst and quiescent galaxies at this epoch are very similar, in contrast with their effective radii. The structural properties of z$\sim$0.7 post-starburst galaxies match those of quiescent galaxies that formed in the early universe, suggesting that rapid quenching in the present epoch is driven by a similar mechanism to the one at high redshift. △ Less

Submitted 13 June, 2022; v1 submitted 16 March, 2022; originally announced March 2022.

Comments: 28 pages, 19 figures, 2 appendices, updated to accepted ApJ version

Exploring supermassive black hole physics and galaxy quenching across halo mass in FIRE cosmological zoom simulations

Authors: Sarah Wellons, Claude-André Faucher-Giguère, Philip F. Hopkins, Eliot Quataert, Daniel Anglés-Alcázar, Robert Feldmann, Christopher C. Hayward, Dušan Kereš, Kung-Yi Su, Andrew Wetzel

Abstract: Feedback from accreting supermassive black holes (SMBHs) is thought to be a primary driver of quenching in massive galaxies, but the best way to implement SMBH physics into galaxy formation simulations remains ambiguous. As part of the Feedback in Realistic Environments (FIRE) project, we explore the effects of different modeling choices for SMBH accretion and feedback in a suite of $\sim500$ cosm… ▽ More Feedback from accreting supermassive black holes (SMBHs) is thought to be a primary driver of quenching in massive galaxies, but the best way to implement SMBH physics into galaxy formation simulations remains ambiguous. As part of the Feedback in Realistic Environments (FIRE) project, we explore the effects of different modeling choices for SMBH accretion and feedback in a suite of $\sim500$ cosmological zoom-in simulations across a wide range of halo mass (10^10-10^13 Msun). Within the suite, we vary the numerical schemes for BH accretion and feedback, the accretion efficiency, and the strength of mechanical, radiative, and cosmic ray feedback independently. We then compare the outcomes to observed galaxy scaling relations. We find several models that satisfy the observational constraints, and for which the energetics in different feedback channels are physically plausible. Interestingly, cosmic rays accelerated by SMBHs play an important role in many successful models. However, it is non-trivial to reproduce scaling relations across halo mass, and many model variations produce qualitatively incorrect results regardless of parameter choices. The growth of stellar and BH mass are closely related: for example, over-massive BHs tend to over-quench galaxies. BH mass is most strongly affected by the choice of accretion efficiency in high-mass halos, but by feedback efficiency in low-mass halos. The amount of star formation suppression by SMBH feedback in low-mass halos is determined primarily by the time-integrated feedback energy. For massive galaxies, the "responsiveness" of a model (i.e. how quickly and powerfully the BH responds to gas available for accretion) is an additional important factor for quenching. △ Less

Submitted 11 March, 2022; originally announced March 2022.

Comments: 18 pages, 9 figures. Submitted to MNRAS, comments welcome

FIRE-3: Updated Stellar Evolution Models, Yields, & Microphysics and Fitting Functions for Applications in Galaxy Simulations

Authors: Philip F. Hopkins, Andrew Wetzel, Coral Wheeler, Robyn Sanderson, Michael Y. Grudic, Omid Sameie, Michael Boylan-Kolchin, Matthew Orr, Xiangcheng Ma, Claude-Andre Faucher-Giguere, Dusan Keres, Eliot Quataert, Kung-Yi Su, Jorge Moreno, Robert Feldmann, James S. Bullock, Sarah R. Loebman, Daniel Angles-Alcazar, Jonathan Stern, Lina Necib, Christopher C. Hayward

Abstract: Increasingly, uncertainties in predictions from galaxy formation simulations (at sub-Milky Way masses) are dominated by uncertainties in stellar evolution inputs. In this paper, we present the full set of updates from the FIRE-2 version of the Feedback In Realistic Environments (FIRE) project code, to the next version, FIRE-3. While the transition from FIRE-1 to FIRE-2 focused on improving numeric… ▽ More Increasingly, uncertainties in predictions from galaxy formation simulations (at sub-Milky Way masses) are dominated by uncertainties in stellar evolution inputs. In this paper, we present the full set of updates from the FIRE-2 version of the Feedback In Realistic Environments (FIRE) project code, to the next version, FIRE-3. While the transition from FIRE-1 to FIRE-2 focused on improving numerical methods, here we update the stellar evolution tracks used to determine stellar feedback inputs, e.g. stellar mass-loss (O/B and AGB), spectra (luminosities and ionization rates), and supernova rates (core-collapse and Ia), as well as detailed mass-dependent yields. We also update the low-temperature cooling and chemistry, to enable improved accuracy at $T \lesssim 10^{4}\,$K and densities $n\gg 1\,{\rm cm^{-3}}$, and the meta-galactic ionizing background. All of these synthesize newer empirical constraints on these quantities and updated stellar evolution and yield models from a number of groups, addressing different aspects of stellar evolution. To make the updated models as accessible as possible, we provide fitting functions for all of the relevant updated tracks, yields, etc, in a form specifically designed so they can be directly 'plugged in' to existing galaxy formation simulations. We also summarize the default FIRE-3 implementations of 'optional' physics, including spectrally-resolved cosmic rays and supermassive black hole growth and feedback. △ Less

Submitted 28 February, 2022; originally announced March 2022.

Comments: 26 pages, 12 figures. Submitted to MNRAS. Comments welcome

Public data release of the FIRE-2 cosmological zoom-in simulations of galaxy formation

Authors: Andrew Wetzel, Christopher C. Hayward, Robyn E. Sanderson, Xiangcheng Ma, Daniel Angles-Alcazar, Robert Feldmann, T. K Chan, Kareem El-Badry, Coral Wheeler, Shea Garrison-Kimmel, Farnik Nikakhtar, Nondh Panithanpaisal, Arpit Arora, Alexander B. Gurvich, Jenna Samuel, Omid Sameie, Viraj Pandya, Zachary Hafen, Cameron Hummels, Sarah Loebman, Michael Boylan-Kolchin, James S. Bullock, Claude-Andre Faucher-Giguere, Dusan Keres, Eliot Quataert , et al. (1 additional authors not shown)

Abstract: We describe a public data release of the FIRE-2 cosmological zoom-in simulations of galaxy formation, available at http://flathub.flatironinstitute.org/fire, from the Feedback In Realistic Environments (FIRE) project. FIRE-2 simulations achieve parsec-scale resolution to explicitly model the multi-phase interstellar medium while implementing direct models for stellar evolution and feedback, includ… ▽ More We describe a public data release of the FIRE-2 cosmological zoom-in simulations of galaxy formation, available at http://flathub.flatironinstitute.org/fire, from the Feedback In Realistic Environments (FIRE) project. FIRE-2 simulations achieve parsec-scale resolution to explicitly model the multi-phase interstellar medium while implementing direct models for stellar evolution and feedback, including stellar winds, core-collapse and Ia supernovae, radiation pressure, photoionization, and photoelectric heating. We release complete snapshots from 3 suites of simulations. The first comprises 20 simulations that zoom in on 14 Milky Way-mass galaxies, 5 SMC/LMC-mass galaxies, and 4 lower-mass galaxies including 1 ultra-faint; we release 39 snapshots across z = 0 - 10. The second comprises 4 massive galaxies, with 19 snapshots across z = 1 - 10. Finally, a high-redshift suite comprises 22 simulations, with 11 snapshots across z = 5 - 10. Each simulation also includes dozens of resolved lower-mass (satellite) galaxies in its zoom-in region. Snapshots include all stored properties for all dark matter, gas, and star particles, including 11 elemental abundances for stars and gas, and formation times (ages) of star particles. We also release accompanying (sub)halo catalogs, which include galaxy properties and member star particles. For the simulations to z = 0, including all Milky Way-mass galaxies, we release the formation coordinates and an "ex-situ" flag for all star particles, pointers to track particles across snapshots, catalogs of stellar streams, and multipole basis expansions for the halo mass distributions. We describe publicly available python packages for reading and analyzing these simulations. △ Less

Submitted 29 March, 2023; v1 submitted 14 February, 2022; originally announced February 2022.

Comments: 22 pages. Accepted for publication in ApJS, matches the published version, with updated references. Now includes up to 39 snapshots per simulation and formation coordinates for all star particles. Data available at http://flathub.flatironinstitute.org/fire

Galaxies lacking dark matter produced by close encounters in a cosmological simulation

Authors: Jorge Moreno, Shany Danieli, James S. Bullock, Robert Feldmann, Philip F. Hopkins, Onur Catmabacak, Alexander Gurvich, Alexandres Lazar, Courtney Klein, Cameron B. Hummels, Zachary Hafen, Francisco J. Mercado, Sijie Yu, Fangzhou Jiang, Coral Wheeler, Andrew Wetzel, Daniel Angles-Alcazar, Michael Boylan-Kolchin, Eliot Quataert, Claude-Andre Faucher-Giguere, Dusan Keres

Abstract: The standard cold dark matter plus cosmological constant model predicts that galaxies form within dark-matter haloes, and that low-mass galaxies are more dark-matter dominated than massive ones. The unexpected discovery of two low-mass galaxies lacking dark matter immediately provoked concerns about the standard cosmology and ignited explorations of alternatives, including self-interacting dark ma… ▽ More The standard cold dark matter plus cosmological constant model predicts that galaxies form within dark-matter haloes, and that low-mass galaxies are more dark-matter dominated than massive ones. The unexpected discovery of two low-mass galaxies lacking dark matter immediately provoked concerns about the standard cosmology and ignited explorations of alternatives, including self-interacting dark matter and modified gravity. Apprehension grew after several cosmological simulations using the conventional model failed to form adequate numerical analogues with comparable internal characteristics (stellar masses, sizes, velocity dispersions and morphologies). Here we show that the standard paradigm naturally produces galaxies lacking dark matter with internal characteristics in agreement with observations. Using a state-of-the-art cosmological simulation and a meticulous galaxy-identification technique, we find that extreme close encounters with massive neighbours can be responsible for this. We predict that approximately 30 percent of massive central galaxies (with at least 1e11 solar masses in stars) harbour at least one dark-matter-deficient satellite (with 1e8 - 1e9 solar masses in stars). This distinctive class of galaxies provides an additional layer in our understanding of the role of interactions in shaping galactic properties. Future observations surveying galaxies in the aforementioned regime will provide a crucial test of this scenario. △ Less

Submitted 11 February, 2022; originally announced February 2022.

Comments: 55 pages, 4 figures, 13 supplementary figures. 29 pages, 4 figures. Accepted for publication in Nature Astronomy. To appear on 14-February-2022. Published version: https://www.nature.com/articles/s41550-021-01598-4

Hot-mode accretion and the physics of thin-disk galaxy formation

Authors: Zachary Hafen, Jonathan Stern, James Bullock, Alex B. Gurvich, Sijie Yu, Claude-Andre Faucher-Giguere, Drummond B. Fielding, Daniel Angles-Alcazar, Eliot Quataert, Andrew Wetzel, Tjitske Starkenburg, Michael Boylan-Kolchin, Jorge Moreno, Robert Feldmann, Kareem El-Badry, T. K. Chan, Cameron Trapp, Dusan Keres, Philip F. Hopkins

Abstract: We use FIRE simulations to study disk formation in z~0, Milky Way-mass galaxies, and conclude that a key ingredient for the formation of thin stellar disks is the ability for accreting gas to develop an aligned angular momentum distribution via internal cancellation *prior* to joining the galaxy. Among galaxies with a high fraction (>70%) of their young stars in a thin disk (h/R~0.1) we find that:… ▽ More We use FIRE simulations to study disk formation in z~0, Milky Way-mass galaxies, and conclude that a key ingredient for the formation of thin stellar disks is the ability for accreting gas to develop an aligned angular momentum distribution via internal cancellation *prior* to joining the galaxy. Among galaxies with a high fraction (>70%) of their young stars in a thin disk (h/R~0.1) we find that: (i) hot, virial-temperature gas dominates the inflowing gas mass on halo scales (>~20 kpc), with radiative losses offset by compression heating; (ii) this hot accretion proceeds until angular momentum support slows inward motion, at which point the gas cools to T~10^4 K or less; (iii) prior to cooling, the accreting gas develops an angular momentum distribution that is aligned with the galaxy disk, and while cooling transitions from a quasi-spherical spatial configuration to a more flattened, disk-like configuration. We show that the existence of this "rotating cooling flow" accretion mode is strongly correlated with the fraction of stars forming in a thin disk among a sample of 17 z~0 galaxies spanning a halo mass range of 10^10.5 solar masses to 10^12 solar masses, or a stellar mass range 10^8 solar masses to 10^11 solar masses. Notably, galaxies with a thick disk or irregular morphology do not undergo significant angular momentum alignment of gas prior to accretion and show no correspondence between halo gas cooling and flattening. Our results suggest that rotating cooling flows (or, more generally, rotating subsonic flows) that become coherent and angular momentum-supported prior to accretion onto the galaxy are likely a necessary condition for the formation of thin, star-forming disk galaxies in a LambdaCDM universe. △ Less

Submitted 6 June, 2022; v1 submitted 18 January, 2022; originally announced January 2022.

Comments: 17 pages, 11 figures. Comments welcome. Accompanying interactive visualization at zhafen.github.io/rotating-cooling-flows

The galaxy-halo size relation of low-mass galaxies in FIRE

Authors: Eric Rohr, Robert Feldmann, James Bullock, Onur Çatmabacak, Michael Boylan-Kolchin, Claude-André Faucher-Giguère, Dušan Kereš, Lichen Liang, Jorge Moreno, Andrew Wetzel

Abstract: Galaxy sizes correlate closely with the sizes of their parent dark matter haloes, suggesting a link between halo formation and galaxy growth. However, the precise nature of this relation and its scatter remains to be understood fully, especially for low-mass galaxies. We analyse the galaxy-halo size relation for low-mass ($M_\star \sim 10^{7-9} {\rm M_\odot}$) central galaxies over the past 12.5 b… ▽ More Galaxy sizes correlate closely with the sizes of their parent dark matter haloes, suggesting a link between halo formation and galaxy growth. However, the precise nature of this relation and its scatter remains to be understood fully, especially for low-mass galaxies. We analyse the galaxy-halo size relation for low-mass ($M_\star \sim 10^{7-9} {\rm M_\odot}$) central galaxies over the past 12.5 billion years with the help of cosmological volume simulations (FIREbox) from the Feedback in Realistic Environments (FIRE) project. We find a nearly linear relationship between the half-stellar mass galaxy size $R_{1/2}$ and the parent dark matter halo virial radius $R_{\rm vir}$. This relation evolves only weakly since redshift $z = 5$: $R_{1/2} {\rm kpc} = (0.053\pm0.002)(R_{\rm vir}/35 {\rm kpc})^{0.934\pm0.054}$, with a nearly constant scatter $\langle σ\rangle = 0.084 [{\rm dex}]$. Whilst this ratio is similar to what is expected from models where galaxy disc sizes are set by halo angular momentum, the low-mass galaxies in our sample are not angular momentum supported, with stellar rotational to circular velocity ratios $v_{\rm rot} / v_{\rm circ} \sim 0.15$. Introducing redshift as another parameter to the GHSR does not decrease the scatter. Furthermore, this scatter does not correlate with any of the halo properties we investigate -- including spin and concentration -- suggesting that baryonic processes and feedback physics are instead critical in setting the scatter in the galaxy-halo size relation. Given the relatively small scatter and the weak dependence of the galaxy-halo size relation on redshift and halo properties for these low-mass central galaxies, we propose using galaxy sizes as an independent method from stellar masses to infer halo masses. △ Less

Submitted 16 April, 2024; v1 submitted 9 December, 2021; originally announced December 2021.

Comments: 14 pages, 9 figures, 3 tables + 3 appendices (5 pages, 5 figures, 2 tables). Accepted for publication in MNRAS. Reuploaded after proofing stage

Journal ref: MNRAS, 510, 3967-3985 (2022)

SQuIGGLE: Studying Quenching in Intermediate-z Galaxies -- Gas, AnguLar Momentum, and Evolution

Authors: Katherine A. Suess, Mariska Kriek, Rachel Bezanson, Jenny E. Greene, David Setton, Justin S. Spilker, Robert Feldmann, Andy D. Goulding, Benjamin D. Johnson, Joel Leja, Desika Narayanan, Khalil Hall-Hooper, Qiana Hunt, Sidney Lower, Margaret Verrico

Abstract: We describe the SQuIGGLE survey of intermediate-redshift post-starburst galaxies. We leverage the large sky coverage of the SDSS to select ~1300 recently-quenched galaxies at 0.5<z<~0.9 based on their unique spectral shapes. These bright, intermediate-redshift galaxies are ideal laboratories to study the physics responsible for the rapid quenching of star formation: they are distant enough to be u… ▽ More We describe the SQuIGGLE survey of intermediate-redshift post-starburst galaxies. We leverage the large sky coverage of the SDSS to select ~1300 recently-quenched galaxies at 0.5<z<~0.9 based on their unique spectral shapes. These bright, intermediate-redshift galaxies are ideal laboratories to study the physics responsible for the rapid quenching of star formation: they are distant enough to be useful analogs for high-redshift quenching galaxies, but low enough redshift that multi-wavelength follow-up observations are feasible with modest telescope investments. We use the Prospector code to infer the stellar population properties and non-parametric star formation histories of all galaxies in the sample. We find that SQuIGGLE galaxies are both very massive (M* ~ 10^11.25 Msun) and quenched, with inferred star formation rates <~1Msun/yr, more than an order of magnitude below the star-forming main sequence. The best-fit star formation histories confirm that these galaxies recently quenched a major burst of star formation: >75% of SQuIGGLE galaxies formed at least a quarter of their total stellar mass in the recent burst, which ended just ~200Myr before observation. We find that SQuIGGLE galaxies are on average younger and more burst-dominated than most other z<~1 post-starburst samples. This large sample of bright post-starburst galaxies at intermediate redshift opens a wide range of studies into the quenching process. In particular, the full SQuIGGLE survey will investigate the molecular gas reservoirs, morphologies, kinematics, resolved stellar populations, AGN incidence, and infrared properties of this unique sample of galaxies in order to place definitive constraints on the quenching process. △ Less

Submitted 29 November, 2021; originally announced November 2021.

Comments: 23 pages, 16 figures, accepted to ApJ

Now you see it, now you don't: Star formation truncation precedes the loss of molecular gas by ~100 Myr in massive post-starburst galaxies at z~0.6

Authors: Rachel Bezanson, Justin S. Spilker, Katherine A. Suess, David J. Setton, Robert Feldmann, Jenny E. Greene, Mariska Kriek, Desika Narayanan, Margaret Verrico

Abstract: We use ALMA observations of CO(2-1) in 13 massive ($M_{\star}\gtrsim 10^{11} M_{\odot}$) post-starburst galaxies at $z\sim0.6$ to constrain the molecular gas content in galaxies shortly after they quench their major star-forming episode. The post-starburst galaxies in this study are selected from the Sloan Digital Sky Survey spectroscopic samples (DR14) based on their spectral shapes, as part of t… ▽ More We use ALMA observations of CO(2-1) in 13 massive ($M_{\star}\gtrsim 10^{11} M_{\odot}$) post-starburst galaxies at $z\sim0.6$ to constrain the molecular gas content in galaxies shortly after they quench their major star-forming episode. The post-starburst galaxies in this study are selected from the Sloan Digital Sky Survey spectroscopic samples (DR14) based on their spectral shapes, as part of the SQuIGGLE program. Early results showed that two post-starburst galaxies host large H$_2$ reservoirs despite their low inferred star formation rates. Here we expand this analysis to a larger statistical sample of 13 galaxies. Six of the primary targets (45%) are detected, with $M_{H_2}\gtrsim10^9 M_{\odot}$. Given their high stellar masses, this mass limit corresponds to an average gas fraction of $\langle f_{H_2} \equiv M_{H_2}/M_{\star} \rangle \sim7\%$, or ${\sim}14\%$ using lower stellar masses estimates derived from analytic, exponentially declining star formation histories. The gas fraction correlates with the $D_n4000$ spectral index, suggesting that the cold gas reservoirs decrease with time since burst, as found in local K+A galaxies. Star formation histories derived from flexible stellar population synthesis modeling support this empirical finding: galaxies that quenched $\lesssim 150$ Myr prior to observation host detectable CO(2-1) emission, while older post-starburst galaxies are undetected. The large $\mathrm{H_2}$ reservoirs and low star formation rates in the sample imply that the quenching of star formation precedes the disappearance of the cold gas reservoirs. However, within the following 100-200 Myrs, the SQuIGGLE galaxies require the additional and efficient heating or removal of cold gas to bring their low star formation rates in line with standard $\mathrm{H_2}$ scaling relations. △ Less

Submitted 29 November, 2021; originally announced November 2021.

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

From EMBER to FIRE: predicting high resolution baryon fields from dark matter simulations with Deep Learning

Authors: Mauro Bernardini, Robert Feldmann, Daniel Anglés-Alcázar, Mike Boylan-Kolchin, James Bullock, Lucio Mayer, Joachim Stadel

Abstract: Hydrodynamic simulations provide a powerful, but computationally expensive, approach to study the interplay of dark matter and baryons in cosmological structure formation. Here we introduce the EMulating Baryonic EnRichment (EMBER) Deep Learning framework to predict baryon fields based on dark-matter-only simulations thereby reducing computational cost. EMBER comprises two network architectures, U… ▽ More Hydrodynamic simulations provide a powerful, but computationally expensive, approach to study the interplay of dark matter and baryons in cosmological structure formation. Here we introduce the EMulating Baryonic EnRichment (EMBER) Deep Learning framework to predict baryon fields based on dark-matter-only simulations thereby reducing computational cost. EMBER comprises two network architectures, U-Net and Wasserstein Generative Adversarial Networks (WGANs), to predict two-dimensional gas and HI densities from dark matter fields. We design the conditional WGANs as stochastic emulators, such that multiple target fields can be sampled from the same dark matter input. For training we combine cosmological volume and zoom-in hydrodynamical simulations from the Feedback in Realistic Environments (FIRE) project to represent a large range of scales. Our fiducial WGAN model reproduces the gas and HI power spectra within 10% accuracy down to ~10 kpc scales. Furthermore, we investigate the capability of EMBER to predict high resolution baryon fields from low resolution dark matter inputs through upsampling techniques. As a practical application, we use this methodology to emulate high-resolution HI maps for a dark matter simulation of a L=100 Mpc/h comoving cosmological box. The gas content of dark matter haloes and the HI column density distributions predicted by EMBER agree well with results of large volume cosmological simulations and abundance matching models. Our method provides a computationally efficient, stochastic emulator for augmenting dark matter only simulations with physically consistent maps of baryon fields. △ Less

Submitted 31 March, 2022; v1 submitted 22 October, 2021; originally announced October 2021.

Comments: 21 pages, 12 figures, accepted by MNRAS, comments welcome

Quantum Proton Effects from Density Matrix Renormalization Group Calculations

Authors: Robin Feldmann, Andrea Muolo, Alberto Baiardi, Markus Reiher

Abstract: We recently introduced [J. Chem. Phys. 152 2020, 204103] the nuclear-electronic all-particle density matrix renormalization group method (NEAP-DMRG) to solve the molecular Schrödinger equation, based on a stochastically optimized orbital basis, without invoking the Born-Oppenheimer approximation. In this work, we combine the DMRG with nuclear-electronic Hartree-Fock (NEHF-DMRG), treating nuclei an… ▽ More We recently introduced [J. Chem. Phys. 152 2020, 204103] the nuclear-electronic all-particle density matrix renormalization group method (NEAP-DMRG) to solve the molecular Schrödinger equation, based on a stochastically optimized orbital basis, without invoking the Born-Oppenheimer approximation. In this work, we combine the DMRG with nuclear-electronic Hartree-Fock (NEHF-DMRG), treating nuclei and electrons on the same footing. Inter- and intra-species correlations are described within the DMRG without truncating the excitation degree of the full configuration interaction wave function. We extend the concept of orbital entanglement and mutual information to nuclear-electronic wave functions and demonstrate that they are reliable metrics to detect strong correlation effects. We apply NEHF-DMRG to the HeHHe$^+$ molecular ion, to obtain accurate proton densities, ground-state total energies, and vibrational transition frequencies by comparison with state-of-the-art data obtained with grid-based approaches and modern configuration interaction methods. For HCN, we improve on the accuracy of the latter approaches with respect to both ground-state absolute energy and proton density which is a major challenge for multi-reference nuclear-electronic state-of-the-art methods. △ Less

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

Comments: 55 pages, 8 figures, 5 tables

Journal ref: J. Chem. Theory Comput. 2022, 18, 1, 234-250

Why do more massive stars host larger planets?

Authors: Michael Lozovsky, Ravit Helled, Illaria Pascucci, Caroline Dorn, Julia Venturini, Robert Feldmann

Abstract: It has been suggested that planetary radii increase with the stellar mass, for planets below 6 R$_{\oplus}$ and host below 1 M$_\odot$. In this study, we explore whether this inferred relation between planetary size and the host star's mass can be explained by a larger planetary mass among planets orbiting more massive stars, inflation of the planetary radius due to the difference in stellar irrad… ▽ More It has been suggested that planetary radii increase with the stellar mass, for planets below 6 R$_{\oplus}$ and host below 1 M$_\odot$. In this study, we explore whether this inferred relation between planetary size and the host star's mass can be explained by a larger planetary mass among planets orbiting more massive stars, inflation of the planetary radius due to the difference in stellar irradiation, or different planetary compositions and structures. Using exoplanetary data of planets with measured masses and radii, we investigate the relations between stellar mass and various planetary properties for G- and K- stars, and confirm that more massive stars host larger planets and more massive. We find that the differences in the planetary masses and temperatures are insufficient to explain the measured differences in radii between planets surrounding different stellar types. We show that the larger planetary radii can be explained by a larger fraction of volatile material (H-He atmospheres) among planets surrounding more massive stars. We conclude that planets around more massive stars are larger most probably as a result of larger H-He atmospheres. Our findings imply that planets forming around more massive stars tend to accrete H-He atmospheres more efficiently. △ Less

Submitted 20 July, 2021; originally announced July 2021.

Report number: AA/2021/40563

Journal ref: A&A 652, A110 (2021)