Search

We present spatially resolved, simple stellar population equivalent ages, stellar metallicities, and abundance ratios for 456 massive (⁠|$10.3\lesssim \log (\mathrm{M}_*/\mathrm{M}_\odot)\lesssim 11.8$|⁠) quiescent galaxies at |$0.6\lesssim z\lesssim 1.0$| from the Large Early Galaxy Astrophysics Census, derived using full-spectrum models. Typically, we find flat age and [Mg/Fe] gradients, and negative [Fe/H] gradients, implying iron-rich cores. We also estimate intrinsic [Fe/H] gradients via forward modelling. We examine the observed gradients in three age bins. Younger quiescent galaxies typically have negative [Fe/H] gradients and positive age gradients, possibly indicating a recent central starburst. Additionally, this finding suggests that photometrically measured flat colour gradients in young quiescent galaxies are the result of the positive age and negative metallicity gradients cancelling each other. For older quiescent galaxies, the age gradients become flat and [Fe/H] gradients weaken, though remain negative. Thus, negative colour gradients at older ages are likely driven by metallicity gradients. The diminishing age gradient may result from the starburst fading. Furthermore, the persistence of the [Fe/H] gradients may suggest that the outskirts are simultaneously built up by mergers with lower metallicity satellites. On the other hand, the gradients could be inherited from the star-forming phase, in which case mergers may not be needed to explain our findings. This work illustrates the need for resolved spectroscopy, instead of just photometry, to measure stellar population gradients. Extending these measurements to higher redshift is imperative for understanding how stellar populations in quiescent galaxies are assembled over cosmic time. [ABSTRACT FROM AUTHOR]

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

We compare ionized gas and stellar kinematics of 16 star-forming galaxies (log (M ⋆/ M ⊙) = 9.7–11.2, SFR =6 − 86 M⊙ yr−1) at z ∼ 1 using near-infrared integral field spectroscopy (IFS) of Hα emission from the KMOS3D (the K-band Multi-Object Spectrograph 3D) survey and optical slit spectroscopy of stellar absorption and gas emission from the LEGA-C (Large Early Galaxy Astrophysics Census) survey. Hα is dynamically colder than stars, with higher disc rotation velocities (by ∼45 per cent) and lower disc velocity dispersions (by a factor ∼2). This is similar to trends observed in the local Universe. We find higher rotational support for Hα relative to [O  ii ], potentially explaining systematic offsets in kinematic scaling relations found in the literature. Regarding dynamical mass measurements, for six galaxies with cumulative mass profiles from Jeans Anisotropic Multi-Gaussian Expansion (JAM) models the Hα dynamical mass models agree remarkably well out to ∼10 kpc for all but one galaxy (average |$\Delta M_{\rm dyn}(R_{e,\rm F814W})\lt 0.1$| dex). Simpler dynamical mass estimates based on integrated stellar velocity dispersion are less accurate (standard deviation 0.24 dex). Differences in dynamical mass estimates are larger, for example, for galaxies with stronger misalignments of the Hα kinematic major axis and the photometric position angle, highlighting the added value of IFS observations for dynamics studies. The good agreement between the JAM and the dynamical models based on Hα kinematics at z ∼ 1 corroborates the validity of dynamical mass measurements from Hα IFS observations, which can be more easily obtained for higher redshift galaxies. [ABSTRACT FROM AUTHOR]

We present the first study of spatially integrated higher-order stellar kinematics over cosmic time. We use deep rest-frame optical spectroscopy of quiescent galaxies at redshifts z  = 0.05, 0.3, and 0.8 from the SAMI, MAGPI, and LEGA-C surveys to measure the excess kurtosis h 4 of the stellar velocity distribution, the latter parametrized as a Gauss-Hermite series. Conservatively using a redshift-independent cut in stellar mass (⁠|$M_\star = 10^{11}\, \mathrm{M_\odot }$|⁠) and matching the stellar-mass distributions of our samples, we find 7σ evidence of h 4 increasing with cosmic time, from a median value of 0.019 ± 0.002 at z  = 0.8 to 0.059 ± 0.004 at z  = 0.06. Alternatively, we use a physically motivated sample selection based on the mass distribution of the progenitors of local quiescent galaxies as inferred from numerical simulations; in this case, we find 10σ evidence. This evolution suggests that, over the last 7 Gyr, there has been a gradual decrease in the rotation-to-dispersion ratio and an increase in the radial anisotropy of the stellar velocity distribution, qualitatively consistent with accretion of gas-poor satellites. These findings demonstrate that massive galaxies continue to accrete mass and increase their dispersion support after becoming quiescent. [ABSTRACT FROM AUTHOR]

We present the first statistical study of spatially integrated non-Gaussian stellar kinematics spanning 7 Gyr in cosmic time. We use deep, rest-frame optical spectroscopy of massive galaxies (stellar mass |$M_\star \gt 10^{10.5} \, \mathrm{M_\odot }$|⁠) at redshifts z = 0.05, 0.3, and 0.8 from the SAMI, MAGPI, and LEGA-C surveys, to measure the excess kurtosis h 4 of the stellar velocity distribution, the latter parametrized as a Gauss–Hermite series. We find that at all redshifts where we have large enough samples, h 4 anticorrelates with the ratio between rotation and dispersion, highlighting the physical connection between these two kinematic observables. In addition, and independently from the anticorrelation with rotation-to-dispersion ratio, we also find a correlation between h 4 and M ⋆, potentially connected to the assembly history of galaxies. In contrast, after controlling for mass, we find no evidence of independent correlation between h 4 and aperture velocity dispersion or galaxy size. These results hold for both star-forming and quiescent galaxies. For quiescent galaxies, h 4 also correlates with projected shape, even after controlling for the rotation-to-dispersion ratio. At any given redshift, star-forming galaxies have lower h 4 compared to quiescent galaxies, highlighting the link between kinematic structure and star-forming activity. [ABSTRACT FROM AUTHOR]

We present the results of a search for high-redshift (z > 9) galaxy candidates in the JWST UNCOVER survey, using deep NIRCam and NIRISS imaging in seven bands over ∼45 arcmin2 and ancillary Hubble Space Telescope (HST) observations. The NIRCam observations reach a 5σ limiting magnitude of ∼29.2 AB. The identification of high- z candidates relies on a combination of a dropout selection and photometric redshifts. We find 16 candidates at 9 < z < 12 and three candidates at 12 < z < 13, eight candidates are deemed very robust. Their lensing amplification ranges from μ = 1.2 to 11.5. Candidates have a wide range of (lensing corrected) luminosities and young ages, with low stellar masses [6.8 < log(M⋆/M⊙) < 9.5] and low star formation rates (SFR = 0.2–7 M⊙ yr−1), confirming previous findings in early JWST observations of z > 9. A few galaxies at z ∼ 9−10 appear to show a clear Balmer break between the F 356 W and F 444 W / F 410 M bands, which helps constrain their stellar mass. We estimate blue UV continuum slopes between β = −1.8 and −2.3, typical for early galaxies at z > 9 but not as extreme as the bluest recently discovered sources. We also find evidence for a rapid redshift-evolution of the mass-luminosity relation and a redshift evolution of the UV continuum slope for a given range of intrinsic magnitude, in line with theoretical predictions. These findings suggest that deeper JWST observations are needed to reach the fainter galaxy population at those early epochs, and follow-up spectroscopy will help better constrain the physical properties and star formation histories of a larger sample of galaxies. [ABSTRACT FROM AUTHOR]

We present a framework to build realistic mock spectroscopic observations for state-of-the-art hydrodynamical simulations, using high spectral resolution stellar population models and full radiative transfer treatment with skirt. As a first application, we generate stellar continuum mock observations for the Auriga cosmological zoom simulations emulating integral-field observations from the Sydney–AOO Multi-object Integral Field Spectrograph (SAMI) Galaxy Survey. We perform spectral fitting on our synthetic cubes and compute the resulting rotation velocity (V rot) and velocity dispersion within 1 R e (σ e) for a subset of the Auriga sample. We find that the kinematics produced by Auriga are in good agreement with the observations from the SAMI galaxy survey after taking into account the effects of dust and the systematics produced by the observation limitations. We also explore the effects of seeing convolution, inclination, and attenuation on the line-of-sight velocity distribution. For highly inclined galaxies, these effects can lead to an artificial decrease in the measured V / σ by nearly a factor of two (after inclination correction). We also demonstrate the utility of our method for high-redshift galaxies by emulating spatially resolved continuum spectra from the Large Early Galaxy Census (LEGA-C) survey and, looking forward, the European Extremely Large Telescope (E-ELT) High Angular Resolution Monolithic Optical and Near-infrared Integral field spectrograph (HARMONI). Our framework represents a crucial link between the ground truth for stellar populations and kinematics in simulations and the observed stellar continuum observations at low and high redshift. [ABSTRACT FROM AUTHOR]

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 × 1010 M⊙ 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 active galactic nucleus 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 × 1011 M⊙. We make predictions for the first quenched galaxies to be located by the JWST. [ABSTRACT FROM AUTHOR]

We use the EAGLE cosmological simulations to perform a comprehensive and systematic analysis of the z  = 0.1 fundamental plane (FP), the tight relation between galaxy size, mass, and velocity dispersion. We first measure the total mass and velocity dispersion (including both random and rotational motions) within the effective radius to show that simulated galaxies obey a total mass FP that is very close to the virial relation (⁠|$\lt 10{{\ \rm per\ cent}}$| deviation), indicating that the effects of non-homology are weak. When we instead use the stellar mass, we find a strong deviation from the virial plane, which is driven by variations in the dark matter content. The dark matter fraction is a smooth function of the size and stellar mass, and thereby sets the coefficients of the stellar mass FP without substantially increasing the scatter. Hence, both star-forming and quiescent galaxies obey the same FP, with equally low scatter (⁠|$0.02\,$| dex). We employ simulations with a variable stellar initial mass function (IMF) to show that IMF variations have a modest additional effect on this FP. Moreover, when we use luminosity-weighted mock observations of the size and spatially integrated velocity dispersion, the inferred FP changes only slightly. However, the scatter increases significantly, due to the luminosity-weighting and line-of-sight projection of the velocity dispersions, and measurement uncertainties on the half-light radii. Importantly, we find significant differences between the simulated FP and observations, which likely reflects a systematic difference in the stellar mass distributions. Therefore, we suggest the stellar mass FP offers a simple test for cosmological simulations, requiring minimal post-processing of simulation data. [ABSTRACT FROM AUTHOR]

We investigate changes in stellar population age and metallicity ([Z/H]) scaling relations for quiescent galaxies from intermediate redshift (0.60 ≤   |$z$| ≤ 0.76) using the LEGA-C Survey to low redshift (0.014 ≤   |$z$| ≤ 0.10) using the SAMI Galaxy Survey. Specifically, we study how the spatially integrated global age and metallicity of individual quiescent galaxies vary in the mass–size plane, using the stellar mass M * and a dynamical mass proxy derived from the virial theorem M D  ∝ σ2  R e. We find that, similarly to at low redshift, the metallicity of quiescent galaxies at 0.60 ≤   |$z$| ≤ 0.76 closely correlates with M / R e (a proxy for the gravitational potential or escape velocity), in that galaxies with deeper potential wells are more metal-rich. This supports the hypothesis that the relation arises due to the gravitational potential regulating the retention of metals by determining the escape velocity for metal-rich stellar and supernova ejecta to escape the system and avoid being recycled into later stellar generations. Conversely, we find no correlation between age and surface density (⁠|$M/R_\mathrm{e}^2$|⁠) at 0.60 ≤   |$z$| ≤ 0.76, despite this relation being strong at low redshift. We consider this change in the age– |$M/R_\mathrm{e}^2$| relation in the context of the redshift evolution of the star-forming and quiescent mass–size relations, and find our results are consistent with galaxies forming more compactly at higher redshifts and remaining compact throughout their evolution. Furthermore, galaxies appear to quench at a characteristic surface density that decreases with decreasing redshift. The |$z$| ∼  0 age– |$M/R_\mathrm{e}^2$| relation is therefore a result of building up the quiescent and star-forming populations with galaxies that formed at a range of redshifts and therefore a range of surface densities. [ABSTRACT FROM AUTHOR]

We use deep, spatially resolved spectroscopy from the Large Early Galaxy Astrophysics Census Survey to study radial variations in the stellar population of 17 spectroscopically selected post-starburst (PSB) galaxies. We use spectral fitting to measure two Lick indices, H δA and Fe 4383 , and find that, on average, PSB galaxies have radially decreasing H δA and increasing Fe 4383  profiles. In contrast, a control sample of quiescent, non-PSB galaxies in the same mass range shows outwardly increasing H δA and decreasing Fe 4383 . The observed gradients are weak (≈−0.2 Å/ R e), mainly due to seeing convolution. A two-SSP (simple stellar population) model suggests that intrinsic gradients are as strong as observed in local PSB galaxies (≈−0.8 Å/ R e). We interpret these results in terms of inside-out growth (for the bulk of the quiescent population) versus star formation occurring last in the centre (for PSB galaxies). At z ≈ 0.8, central starbursts are often the result of gas-rich mergers, as evidenced by the high fraction of PSB galaxies with disturbed morphologies and tidal features (40 per cent). Our results provide additional evidence for multiple paths to quiescence: a standard path, associated with inside-out disc formation and with gradually decreasing star formation activity, without fundamental structural transformation, and a fast path, associated with centrally concentrated starbursts, leaving an inverse age gradient and smaller half-light radius. [ABSTRACT FROM AUTHOR]