Search

We present new results from the ASTRID simulation from $z=3$ to $z=0.5$, covering the epoch of cosmic noon. The galaxy stellar mass function, as well as the black hole mass and luminosity functions in ASTRID, exhibit good agreement with recent observational constraints. We study the $M_{\rm BH}$-$M_*$ scaling relation and its connections to AGN luminosity, galaxy color, and star formation rate, demonstrating that AGN feedback plays a crucial role in the quenching of massive galaxies ($M_*>10^{10.5} M_{\odot}$). Although AGN feedback suppresses star formation through quenching, AGN-host galaxies still exhibit statistically higher levels of star formation compared to inactive ones, due to the positive correlation between AGN activity and star formation, both fueled by a shared gas reservoir. The fraction of quiescent galaxies in ASTRID increases with both galaxy mass and redshift evolution, aligning well with observational trends. We find that different quenching mechanisms can leave distinct morphological imprints on quenched galaxies. Massive, compact quiescent galaxies typically experience shorter quenching timescales, have younger central regions, and host overmassive black holes. This is usually due to a compaction-like quenching mechanism that funnels gas into the galaxy center, leading to starbursts and triggering AGN kinetic feedback. In contrast, quiescent galaxies with more diffuse morphologies generally experience `inside-out' quenching, which is characterized by older central regions compared to the outskirts. These galaxies typically experience longer quenching timescales due to quenching processes operating on a larger halo scale, which gradually deplete the galactic star-forming gas. Data of the \astrid simulation down to $z=0.5$ is available at \url{https://astrid.psc.edu}., Comment: Summited to APJ, comments are welcome

Understanding the chemical enrichment of different elements is crucial to gaining a complete picture of galaxy chemical evolution. In this study, we present a new sample of 46 low-redshift, low-mass star-forming galaxies at $M_*\sim 10^{8-10}M_{\odot}$ along with two quiescent galaxies at $M_*\sim 10^{8.8}M_{\odot}$ observed with the Keck Cosmic Web Imager (KCWI), aiming to investigate the chemical evolution of galaxies in the transition zone between Local Group satellites and massive field galaxies. We develop a novel method to simultaneously determine stellar abundances of iron and magnesium in star-forming galaxies. With the gas-phase oxygen abundance (O/H)$_{\rm g}$ measured using the strong line method, we are able to make the first-ever apples-to-apples comparison of $\alpha$ elements in the stars and the ISM. We find that the [Mg/H]$_*$-[O/H]$_{\rm g}$ relation is much tighter than the [Fe/H]$_*$-[O/H]$_{\rm g}$ relation, which can be explained by the similar production processes of $\alpha$ elements. Most galaxies in our sample exhibit higher [O/H]$_{\rm g}$ than [Fe/H]$_*$ and [Mg/H]$_*$. In addition, we construct mass-metallicity relations (MZRs) measured as three different elements (Fe$_*$, Mg$_*$, O$_{\rm g}$). Compared to the gas O-MZR, the stellar Fe- and Mg-MZRs show larger scatter driven by variations in specific star formation rates (sSFR), with star-forming galaxies exhibiting higher sSFR and lower stellar abundances at fixed mass. The excess of [O/H]$_{\rm g}$ compared to stellar abundances as well as the anti-correlation between sSFR and stellar abundance suggests that galaxy quenching of intermediate-mass galaxies at $M_*\sim 10^{8-10}M_{\odot}$ is primarily driven by starvation., Comment: 36 pages, 16 figures, 4 tables; accepted for publication in ApJ. Main results in Figure 7, 10 and 11

We analyze ionized gas emission lines in deep rest-frame optical spectra of 16 quiescent galaxies at redshift $1.7

Recent JWST observations have uncovered an unexpectedly large population of massive quiescent galaxies at $z>3$. Using the cosmological simulations IllustrisTNG and ASTRID, we identify analogous galaxies and investigate their abundance, formation, quenching mechanisms, and post-quenching evolution for stellar masses $9.5 < \log_{10}{(M_\star/{\rm M}_\odot)} < 12$. We apply three different quenching definitions and find that both simulations significantly underestimate the comoving number density of quenched massive galaxies at $z \gtrsim 3$ compared to JWST observations by up to $\sim 2$ dex. This fact highlights the necessity for improved physical models of AGN feedback in galaxy formation simulations. In both simulations, the high-$z$ quenched massive galaxies often host overmassive central black holes above the standard $M_{BH}-M_\star$ relation, implying that the AGN feedback plays a crucial role in quenching galaxies in the early Universe. The typical quenching timescales for these galaxies are $\sim 200-600$ Myr. IllustrisTNG primarily employs AGN kinetic feedback, while ASTRID relies on AGN thermal feedback, which is less effective and has a longer quenching timescale. We also study the post-quenching evolution of the high-$z$ massive quiescent galaxies and find that many experience subsequent reactivation of star formation, evolving into primary progenitors of $z=0$ brightest cluster galaxies., Comment: Submitted to The Astrophysical Journal. 13 pages, 12 figures

The complex physics governing nebular emission in galaxies, particularly in the early universe, often defy simple low-dimensional models. This has proven to be a significant barrier in understanding the (often diverse) ionizing sources powering this emission. We present Cue, a highly flexible tool for interpreting nebular emission across a wide range of abundances and ionizing conditions of galaxies at different redshifts. Unlike typical nebular models used to interpret extragalactic nebular emission, our model does not require a specific ionizing spectrum as a source, instead approximating the ionizing spectrum with a 4-part piece-wise power-law. We train a neural net emulator based on the CLOUDY photoionization modeling code and make self-consistent nebular continuum and line emission predictions. Along with the flexible ionizing spectra, we allow freedom in [O/H], [N/O], [C/O], gas density, and total ionizing photon budget. This flexibility allows us to either marginalize over or directly measure the incident ionizing radiation, thereby directly interrogating the source of the ionizing photons in distant galaxies via their nebular emission. Our emulator demonstrates a high accuracy, with $\sim$1% uncertainty in predicting the nebular continuum and $\sim$5% uncertainty in the emission lines. Mock tests suggest Cue is well-calibrated and produces useful constraints on the ionizing spectra when $S/N (\mathrm{H}_\alpha) \gtrsim 10$, and furthermore capable of distinguishing between the ionizing spectra predicted by single and binary stellar models. The compute efficiency of neural networks facilitates future applications of Cue for rapid modeling of the nebular emission in large samples and Monte Carlo sampling techniques., Comment: submitted to ApJ

Massive quiescent galaxies in the young universe are expected to be quenched rapidly, but it is unclear whether they all experience starbursts before quenching and what physical mechanism drives rapid quenching. We study 16 massive quiescent galaxies ($\log(M_\star/M_\odot) > 10$) at $z\sim2$ selected from a representative sample of the Blue Jay survey. We reconstruct their star formation histories by fitting spectral energy distribution models to the JWST/NIRSpec $R\sim1000$ spectra. We find that massive quiescent galaxies can be split into three categories with roughly equal numbers of galaxies according to their SFHs: 1) Relatively old galaxies quenched at early epochs; 2) Galaxies that are rapidly and recently quenched after a flat or bursty formation history (depending on the assumed prior); 3) Galaxies that are rapidly and recently quenched after a major starburst. Most recently quenched galaxies show neutral gas outflows, probed by blueshifted $\rm Na\,I\,D$ absorption, and ionized gas emission, with line ratios consistent with active galactic nucleus (AGN) diagnostics. This suggests that AGN activity drives multi-phase gas outflows, leading to rapid quenching. By tracing back the SFHs of the entire sample, we predict the number density of massive quiescent galaxies at $z=4-6$: $n=3.0\pm1.4\times10^{-5}\,\rm Mpc^{-3}$. The two oldest massive quiescent galaxies in our sample appear to have extremely early formation and quenching ($z\gtrsim6$), possibly descendants of early post-starbursts at $z>3$. These galaxies still show neutral gas reservoirs and low-level star formation, consistent with weak H$\alpha$ emission, perhaps because the ejective AGN feedback that caused rapid quenching has weakened over time., Comment: Submitted to ApJ; 21 pages, 9 figures, 3 tables (excluding appendix)

(Reduced) Using the recent cosmological high-resolution zoom-in simulations, NewHorizon and Galactica, in which the evolution of black hole spin is followed on the fly, we have tracked the cosmic history of a hundred of black holes (BHs) with a mass greater than 2x10^4 Ms. For each of them, we have studied the variations of the three dimensional angle (Psi) subtended between the BH spins and the angular momentum vectors of their host galaxies. The analysis of the individual evolution of the most massive BHs suggests that they are generally passing by three different regimes. First, for a short period after their birth, low mass BHs (<3x10^4 Ms) are rapidly spun up by gas accretion and their spin tends to be aligned with their host galaxy spin. Then follows a second phase in which the accretion of gas onto low mass BHs (<10^5 Ms) is quite chaotic and inefficient, reflecting the complex and disturbed morphologies of forming proto-galaxies at high redshifts. The variations of Psi are rather erratic during this phase and are mainly driven by the rapid changes of the direction of the galaxy angular momentum. Then, in a third and long phase, BHs are generally well settled in the center of galaxies around which the gas accretion becomes much more coherent (>10^5 Ms). In this case, the BH spins tend to be well aligned with the angular momentum of their host galaxy and this configuration is generally stable even though BH merger episodes can temporally induce misalignment. We have also derived the distributions of cos(Psi) at different redshifts and found that BHs and galaxy spins are generally aligned. Finally, based on a Monte Carlo method, we also predict statistics for the 2-d projected spin-orbit angles lambda. In particular, the distribution of lambda traces well the alignment tendency in the 3-d analysis. Such predictions provide an interesting background for future observational analyses., Comment: Accepted for publication in A&A (12/03/2024)

We make an in-depth analysis of different AGN jet models' signatures, inducing quiescence in galaxies with a halo mass of $10^{12} M_\odot$. Three jet models, including cosmic ray-dominant, hot thermal, and precessing kinetic jets, are studied at two energy flux levels each, compared to a jet-free, stellar feedback-only simulation. We examine the distribution of Mg II, O VI, and O VIII ions, alongside gas temperature and density profiles. Low-energy ions, like Mg II, concentrate in the ISM, while higher energy ions, e.g., O VIII, prevail at the AGN jet cocoon's edge. High-energy flux jets display an isotropic ion distribution with lower overall density. High-energy thermal or cosmic ray jets pressurize at smaller radii, significantly suppressing core density. The cosmic ray jet provides extra pressure support, extending cool and warm gas distribution. A break in the ion-to-mass ratio slope in O VI and O VIII is demonstrated in the ISM-to-CGM transition (between 10-30 kpc), growing smoothly towards the CGM at greater distances., Comment: 17 pages, 13 figures

We use deep JWST/NIRSpec R~1000 slit spectra of 113 galaxies at 1.7 < z < 3.5, selected from the mass-complete Blue Jay survey, to investigate the prevalence and typical properties of neutral gas outflows at cosmic noon. We detect excess Na I D absorption (beyond the stellar contribution) in 46% of massive galaxies ($\log$ M$_*$/M$_\odot >$ 10), with similar incidence rates in star-forming and quenching systems. Half of the absorption profiles are blueshifted by at least 100 km/s, providing unambiguous evidence for neutral gas outflows. Galaxies with strong Na I D absorption are distinguished by enhanced emission line ratios consistent with AGN ionization. We conservatively measure mass outflow rates of 3 - 100 $M_\odot$ yr$^{-1}$; comparable to or exceeding ionized gas outflow rates measured for galaxies at similar stellar mass and redshift. The outflows from the quenching systems (log(sSFR)[yr$^{-1}$] $\lesssim$ -10) have mass loading factors of 4 - 360, and the energy and momentum outflow rates exceed the expected injection rates from supernova explosions, suggesting that these galaxies could possibly be caught in a rapid blowout phase powered by the AGN. Our findings suggest that AGN-driven ejection of cold gas may be a dominant mechanism for fast quenching of star formation at z~2., Comment: 17 pages, 8 figures, accepted for publication in MNRAS

We report the serendipitous discovery of an extended stellar halo surrounding the low-mass galaxy Ark 227 ($M_\ast=5\times10^9 M_\odot$; d=35 Mpc) in deep JWST NIRCam imaging from the Blue Jay Survey. The F200W-F444W color provides robust star-galaxy separation, enabling the identification of stars at very low density. By combining resolved stars at large galactocentric distances with diffuse emission from NIRCam and Dragonfly imaging at smaller distances, we trace the surface brightness and color profiles of this galaxy over the entire extent of its predicted dark matter halo, from 0.1-100 kpc. Controlled N-body simulations have predicted that minor mergers create "accretion shelves" in the surface brightness profile at large radius. We observe such a feature in Ark 227 at 10-20 kpc, which, according to models, could be caused by a merger with total mass ratio 1:10. The metallicity declines over this radial range, further supporting the minor merger scenario. There is tentative evidence of a second shelf at $\mu_V\approx 35$ mag arcsec$^{-2}$ extending from 50-100 kpc, along with a corresponding drop in metallicity. The stellar mass in this outermost envelope is $\approx10^7M_\odot$. These results suggest that Ark 227 experienced multiple mergers with a spectrum of lower-mass galaxies -- a scenario that is broadly consistent with the hierarchical growth of structure in a cold dark matter-dominated universe. Finally, we identify an ultra-faint dwarf associated with Ark 227 with $M_\ast\approx10^5 M_\odot$ and $\mu_{V,e}=28.1$ mag arcsec$^{-2}$, demonstrating that JWST is capable of detecting very low-mass dwarfs to distances of at least ~30 Mpc., Comment: 10 pages, 9 figures, submitted to ApJ

Large-scale outflows driven by supermassive black holes are thought to play a fundamental role in suppressing star formation in massive galaxies. However, direct observational evidence for this hypothesis is still lacking, particularly in the young universe where star formation quenching is remarkably rapid, thus requiring effective removal of gas as opposed to slow gas heating. While outflows of ionized gas are commonly detected in massive distant galaxies, the amount of ejected mass is too small to be able to suppress star formation. Gas ejection is expected to be more efficient in the neutral and molecular phases, but at high redshift these have only been observed in starbursts and quasars. Here we report JWST spectroscopy of a massive galaxy experiencing rapid quenching at redshift z=2.445. We detect a weak outflow of ionized gas and a powerful outflow of neutral gas, with a mass outflow rate that is sufficient to quench the star formation. Neither X-ray or radio activity are detected; however, the presence of a supermassive black hole is suggested by the properties of the ionized gas emission lines. We thus conclude that supermassive black holes are able to rapidly suppress star formation in massive galaxies by efficiently ejecting neutral gas., Comment: 21 pages, 9 figures. Published in Nature

Thick disks are a prevalent feature observed in numerous disk galaxies including our own Milky Way. Their significance has been reported to vary widely, ranging from a few to 100% of the disk mass, depending on the galaxy and the measurement method. We use the NewHorizon simulation which has high spatial and stellar mass resolutions to investigate the issue of thick disk mass fraction. We also use the NewHorizon2 simulation that was run on the same initial conditions but additionally traced nine chemical elements. Based on a sample of 27 massive disk galaxies with M* > 10^10 M_{\odot} in NewHorizon, the contribution of the thick disk was found to be 34 \pm 15% in r-band luminosity or 48 \pm 13% in mass to the overall galactic disk, which seems in agreement with observational data. The vertical profiles of 0, 22, and 5 galaxies are best fitted by 1, 2, or 3 sech2 components, respectively. The NewHorizon2 data show that the selection of thick disk stars based on a single [{\alpha}/Fe] cut is severely contaminated by stars of different kinematic properties while missing a bulk of kinematically thick disk stars. Vertical luminosity profile fits recover the key properties of thick disks reasonably well. The majority of stars are born near the galactic mid-plane with high circularity and get heated with time via fluctuation in the force field. Depending on the star formation and merger histories, galaxies may naturally develop thick disks with significantly different properties., Comment: 28 pages, 24 figures. Submitted to ApJS

We identify the first quiescent galaxies in TNG300, the largest volume of the IllustrisTNG cosmological simulation suite, and explore their quenching processes and time evolution to z=0. We find that the first quiescent galaxies with stellar masses M_* > 3 x 10^{10} M_sun and specific star formation rates sSFR < 10^{-11} yr^{-1} emerge at z~4.2 in TNG300. Suppression of star formation in these galaxies begins with a thermal mode of AGN feedback at z~6, and a kinetic feedback mode acts in each galaxy by z~4.7 to complete the quenching process, which occurs on a time-scale of ~0.35 Gyr. Surprisingly, we find that the majority of these galaxies are not the main progenitors of their z=0 descendants; instead, four of the five galaxies fall into more massive galaxies in subsequent mergers at a range of redshifts 2.5 < z < 0.2. By z=0, these descendants are the centres of galaxy clusters with average stellar masses of 8 x 10^{11} M_sun. We make predictions for the first quenched galaxies to be located by the James Webb Space Telescope (JWST)., Comment: 6 pages, 4 figures

The existence of massive quiescent galaxies at high redshift seems to require rapid quenching, but it is unclear whether all quiescent galaxies have gone through this phase and what physical mechanisms are involved. To study rapid quenching, we use rest-frame colors to select 12 young quiescent galaxies at $z \sim 1.5$. From spectral energy distribution fitting, we find that they all experienced intense starbursts prior to rapid quenching. We confirm this with deep Magellan/FIRE spectroscopic observations for a subset of seven galaxies. Broad emission lines are detected for two galaxies and are most likely caused by AGN activity. The other five galaxies do not show any emission features, suggesting that gas has already been removed or depleted. Most of the rapidly quenched galaxies are more compact than normal quiescent galaxies, providing evidence for a central starburst in the recent past. We estimate an average transition time of $300\,\rm Myr$ for the rapid quenching phase. Approximately $4\%$ of quiescent galaxies at $z=1.5$ have gone through rapid quenching; this fraction increases to $23\%$ at $z=2.2$. We identify analogs in the TNG100 simulation and find that rapid quenching for these galaxies is driven by AGN, and for half of the cases, gas-rich major mergers seem to trigger the starburst. We conclude that these massive quiescent galaxies are not just rapidly quenched but also rapidly formed through a major starburst. We speculate that mergers drive gas inflow towards the central regions and grow supermassive black holes, leading to rapid quenching by AGN feedback., Comment: Submitted to ApJ. Comments are welcome

The Milky Way (MW) galaxy is in focus, thanks to new observational data. Here we shed new light on the MW's past by studying the structural evolution of MW progenitors, which we identify from extragalactic surveys. Specifically, we constrain the stellar-mass growth history (SMGH) of the MW with two methods: ($i$) direct measurement of the MW's star formation history, and ($ii$) assuming the MW is a typical star-forming galaxy that remains on the star-forming main sequence. We select MW progenitors based on these two SMGHs at $z=0.2-2.0$ from the CANDELS/3D-HST data. We estimate the structural parameters (including half-mass radius $r_{50}$ and S\'ersic index) from the stellar-mass profiles. Our key finding is that the progenitors of the MW galaxy grow self-similarly on spatially resolved scales with roughly a constant half-mass radius ($\sim2-3$ kpc) over the past 10 Gyr, while their stellar masses increase by about 1 dex, implying little-to-no inside-out growth. We discover that the radius containing $20\%$ of the stellar mass ($r_{20}$) decreases by $60\%$ between redshifts of $z=2.0$ and $z=0.7$, while the central stellar-mass density ($\Sigma_1$) increases by a factor of 1.3 dex over the same time and the S\'ersic index changes as $n\propto(1+z)^{-1.41\pm0.19}$. This is consistent with an early ($z>1$) formation of a thick disk, followed by the formation of a bar that led to an increase in the mass in the core. The formation and evolution of the thin disk had only little impact on the overall half-mass size. We also show that the constant-size evolution of the MW progenitors challenges semiempirical approaches and numerical simulations., Comment: Published in ApJL; 8 pages, 5 figures; doi:10.3847/2041-8213/ac76c8

Observations have shown that the star-formation activity and the morphology of galaxies are closely related, but the underlying physical connection is not well understood. Using the TNG50 simulation, we explore the quenching and the morphological evolution of the 102 massive quiescent galaxies in the mass range of $10.5<\log(M_{\rm stellar}/M_{\odot})<11.5$ selected at $z=0$. The morphology of galaxies is quantified based on their kinematics, and we measure the quenching timescale of individual galaxies directly from star formation history. We show that galaxies tend to be quenched more rapidly if they: (i) are satellites in massive halos, (ii) have lower star-forming gas fractions, or (iii) inject a larger amount of black hole kinetic feedback energy. By following the global evolutionary pathways, we conclude that quiescent discs are mainly disc galaxies that are recently and slowly quenched. Approximately half of the quiescent ellipticals at $z=0$ are rapidly quenched at higher redshifts while still disc-like. While being quiescent, they gradually become more elliptical mostly by disc heating, yet these ellipticals still retain some degree of rotation. The other half of quiescent ellipticals with the most random motion-dominated kinematics build up large spheroidal components before quenching primarily by mergers, or in some cases, misaligned gas accretion. However, the mergers that contribute to morphological transformation do not immediately quench galaxies in many cases. In summary, we find that quenching and morphological transformation are decoupled. We conclude that the TNG black hole feedback -- in combination with the stochastic merger history of galaxies -- leads to a large diversity of quenching timescales and a rich morphological landscape., Comment: Submitted to MNRAS. 9 Figures 17 Pages. Comments are very welcome!

Ever since a thick disk was proposed to explain the vertical distribution of the Milky Way disk stars, its origin has been a recurrent question. We aim to answer this question by inspecting 19 disk galaxies with stellar mass greater than $10^{10}\,\rm M_\odot$ in recent cosmological high-resolution zoom-in simulations: Galactica and NewHorizon. The thin and thick disks are reasonably reproduced by the simulations with scale heights and luminosity ratios as observed. We then spatially classify the thin and thick disks and find that the thick disk stars are older, metal-poorer, kinematically-hotter, and higher in accreted star fraction, while both disks are dominated by the stars formed in situ. Half of the in-situ stars in the thick disks are formed before the galaxies develop their disks, and the rest are formed in spatially and kinematically thinner disks and then thickened with time by heating. However, the 19 galaxies have various properties and evolutionary routes, highlighting the need for statistically-large samples to draw general conclusions. We conclude from our simulations that the thin and thick disk components are not entirely distinct in terms of formation processes, but rather markers of the evolution of galactic disks. Moreover, as the combined result of the thickening of the existing disk stars and the continued formation of young thin-disk stars, the vertical distribution of stars does not change much after the disks settle, pointing to the modulation of both orbital diffusion and star formation by the same confounding factor: the proximity of galaxies to marginal stability., Comment: Accepted for publication in ApJS