Search

The Euclid mission is generating a vast amount of imaging data in four broadband filters at high angular resolution. This will allow the detailed study of mass, metallicity, and stellar populations across galaxies, which will constrain their formation and evolutionary pathways. Transforming the Euclid imaging for large samples of galaxies into maps of physical parameters in an efficient and reliable manner is an outstanding challenge. We investigate the power and reliability of machine learning techniques to extract the distribution of physical parameters within well-resolved galaxies. We focus on estimating stellar mass surface density, mass-averaged stellar metallicity and age. We generate noise-free, synthetic high-resolution imaging data in the Euclid photometric bands for a set of 1154 galaxies from the TNG50 cosmological simulation. The images are generated with the SKIRT radiative transfer code, taking into account the complex 3D distribution of stellar populations and interstellar dust attenuation. We use a machine learning framework to map the idealised mock observational data to the physical parameters on a pixel-by-pixel basis. We find that stellar mass surface density can be accurately recovered with a $\leq 0.130 {\rm \,dex}$ scatter. Conversely, stellar metallicity and age estimates are, as expected, less robust, but still contain significant information which originates from underlying correlations at a sub-kpc scale between stellar mass surface density and stellar population properties.

Context. The UVJ color-color diagram is a widely used diagnostic to separate star-forming and quiescent galaxies. Observational data from photometric surveys reveal a strong stellar mass trend, with higher-mass star-forming galaxies being systematically more dust-reddened. Aims. We analyze the UVJ diagram in the TNG100 cosmological simulation at cosmic noon ($z\approx2$). Specifically, we focus on the trend between UVJ colors and mass which has so far not been reproduced in any cosmological simulation. Methods. We apply the SKIRT dust radiative transfer code to the TNG100 simulation to generate rest-frame UVJ fluxes. These UVJ colors are then compared to observational data from several well-studied extragalactic fields from the CANDELS/3D-HST programs, augmented by recent JWST/NIRCam photometry. Results. Quiescent and low-mass ($M_\star\lesssim10^{10.5}\,\mathrm{M}_\odot$) galaxies at cosmic noon do not require significant levels of dust reddening, as opposed to massive ($M_\star\gtrsim10^{11}\,\mathrm{M}_\odot$) star-forming galaxies. An extensive range of possible dust models fall short of the required dust reddening in V-J color for massive star-forming galaxies, with the simulated galaxies being too blue by $\approx0.9\,\mathrm{mag}$. Conclusions. We find that only variations in the star-to-dust geometries of the simulated galaxies can yield V-J colors that are red enough to match the observations. A toy model with isolated dust screens around younger stellar populations (with ages below $\sim1\,\mathrm{Gyr}$) can reproduce the observational data, while all conventional dust radiative transfer models (where the dust distribution follows the metals in the interstellar medium) fail to achieve the required V-J colors., Comment: Main text 17 pages, 12 figures. Resubmitted to A&A after first referee report. Comments warmly welcome!

We present a study of the size growth of the red sequence between $0.510$, whereas half-mass radii are derived from the F444W profiles together with the F277W-F444W color-$M_*$/L relation. We investigate the dependence of the ratio $r_{e, \mathrm{mass}}/r_{e, \mathrm{light}}$ on redshift, stellar mass, and the wavelength used to measure $r_{e, \mathrm{light}}$, also separating the sample into younger and older quiescent galaxies. Our data demonstrate that rest-frame infrared sizes accurately trace mass-weighted sizes while sizes measured at rest-frame optical wavelengths (0.5-0.7$\mu$m) are 0.1-0.2 dex larger, with only minor variations in redshift. We find that the average size of young quiescent galaxies agrees with that of old quiescent galaxies at intermediate masses, $10<$log($M_*/M_{\odot}$)$<11$, within their respective uncertainties in all observed-frame half-light, rest-frame half-light and half-mass radius measurements. At face value, our results point to a combination of progenitor bias and minor mergers driving the size growth of intermediate-mass quiescent galaxies at $0.5

The James Webb Space Telescope (JWST) is unveiling the rest-frame near-IR structure of galaxies. We measure the evolution with redshift of the rest-frame optical and near-IR S\'ersic index ($n$), and examine the dependence on stellar mass and star-formation activity across the redshift range $0.5\leq z\leq2.5$. We infer rest-frame near-IR S\'ersic profiles for $\approx 15.000$ galaxies in publicly available NIRCam imaging mosaics from the COSMOS-Web and PRIMER surveys. We augment these with rest-frame optical S\'ersic indices, previously measured from HST imaging mosaics. The median S\'ersic index evolves slowly or not at all with redshift, except for very high-mass galaxies ($M_\star > 10^{11}~{\text{M}}_\odot$), which show an increase from $n\approx 2.5$ to $n\approx 4$ at $z<1$. High-mass galaxies have higher $n$ than lower-mass galaxies ($M_\star=10^{9.5}~{\text{M}}_\odot$) at all redshifts, with a stronger dependence in the rest-frame near-IR than in the rest-frame optical at $z>1$. This wavelength dependence is caused by star-forming galaxies that have lower optical than near-IR $n$ at z>1 (but not at z<1). Both at optical and near-IR wavelengths, star-forming galaxies have lower $n$ than quiescent galaxies, fortifying the connection between star-formation activity and radial stellar mass distribution. At $z>1$ the median near-IR $n$ varies strongly with star formation activity, but not with stellar mass. The scatter in near-IR $n$ is higher in the green valley (0.25 dex) than on the star-forming sequence and among quiescent galaxies (0.18 dex) -- this trend is not seen in the optical because dust and young stars contribute to the variety in optical light profiles. Our newly measured rest-frame near-IR radial light profiles motivate future comparisons with radial stellar mass profiles of simulated galaxies as a stringent constraint on processes that govern galaxy formation., Comment: Accepted for publication on A&A

Current galaxy formation simulations often approximate star-formation, necessitating models of star-forming regions to produce observables. In the first paper of the series, we introduced TODDLERS, a time-resolved model of UV-mm emission from star-forming regions implemented in the radiative transfer code SKIRT. This work uses SKIRT-TODDLERS to produce synthetic observations, demonstrating its potential through observables related to star-formation and comparing results with existing models in SKIRT. We calculate broadband and line emission maps for 30 Milky Way-like galaxies from the Auriga simulation at z=0. Analyzing FUV and IR data, we calculate kpc-resolved IR correction factors (k_IR), quantifying the ratio of FUV luminosity absorbed by dust to reprocessed IR luminosity. We use IR maps to calculate kpc-scale MIR (8 micron / 24 micron) and FIR (70 micron / 500 micron) colors. H alpha and H beta line maps are used to study the Balmer decrement and dust correction. We also verify the fidelity of our model FIR fine structure lines as SFR indicators. We find that the Auriga integrated UV-mm SEDs show higher FUV/NUV attenuation and lower 24 micron emission when using TODDLERS instead of the existing models in SKIRT, alleviating tensions with observations. The light-weighted mean k_IR increases with aperture and inclination, while its correlation with kpc-resolved sSFR is weaker than literature values. Kpc-scale MIR-FIR colors agree excellently with local observations, with anti-correlation degree varying by galaxy morphology. The Balmer decrement effectively corrects for dust, with attenuation law varying with dust amount. H alpha emission attenuation levels are comparable to high-density regions in state-of-the-art simulations. FIR fine-structure line emission-based luminosity-SFR relations align with global observations, [CII] showing best agreement., Comment: 22 pages, 20 figures. Accepted for publication in Astronomy and Astrophysics on October 21, 2024

Hubble Space Telescope imaging shows that most star-forming galaxies at cosmic noon -- the peak of cosmic star formation history -- appear disk-dominated, leaving the origin of the dense cores in their quiescent descendants unclear. With the James Webb Space Telescope's (JWST) high-resolution imaging to 5 {\mu}m, we can now map the rest-frame near-infrared emission, a much closer proxy for stellar mass distribution, in these massive galaxies. We selected 70 star-forming galaxies with 10$<$log(M)$<$12 and 1.5$<$z$<$3 in the CEERS survey and compare their morphologies in the rest-frame optical to those in the rest-frame near-IR. While the bulk of these galaxies are disk-dominated in 1.5 {\mu}m (rest-frame optical) imaging, they appear more bulge-dominated at 4.4 {\mu}m (rest-frame near-infrared). Our analysis reveals that in massive star-forming galaxies at z$\sim$2, the radial surface brightness profiles steepen significantly, from a slope of $\sim$0.3/dex at 1.5 {\mu}m to $\sim$1.4/dex at 4.4 {\mu}m within radii $<$ 1 kpc. Additionally, we find their total flux contained within the central 1 kpc is approximately 7 times higher in F444W than in F150W. In rest-optical emission, a galaxy's central surface density appears to be the strongest indicator of whether it is quenched or star-forming. Our most significant finding is that at redder wavelengths, the central surface density ratio between quiescent and star-forming galaxies dramatically decreases from $\sim$10 to $\sim$1. This suggests the high central densities associated with galaxy quenching are already in place during the star-forming phase, imposing new constraints on the transition from star formation to quiescence., Comment: 9 pages, 5 figures, published in ApJL

We present new JWST observations of the nearby, prototypical edge-on, spiral galaxy NGC 891. The northern half of the disk was observed with NIRCam in its F150W and F277W filters. Absorption is clearly visible in the mid-plane of the F150W image, along with vertical dusty plumes that closely resemble the ones seen in the optical. A $\sim 10 \times 3~{\rm kpc}^2$ area of the lower circumgalactic medium (CGM) was mapped with MIRI F770W at 12 pc scales. Thanks to the sensitivity and resolution of JWST, we detect dust emission out to $\sim 4$ kpc from the disk, in the form of filaments, arcs, and super-bubbles. Some of these filaments can be traced back to regions with recent star formation activity, suggesting that feedback-driven galactic winds play an important role in regulating baryonic cycling. The presence of dust at these altitudes raises questions about the transport mechanisms at play and suggests that small dust grains are able to survive for several tens of million years after having been ejected by galactic winds in the disk-halo interface. We lay out several scenarios that could explain this emission: dust grains may be shielded in the outer layers of cool dense clouds expelled from the galaxy disk, and/or the emission comes from the mixing layers around these cool clumps where material from the hot gas is able to cool down and mix with these cool cloudlets. This first set of data and upcoming spectroscopy will be very helpful to understand the survival of dust grains in energetic environments, and their contribution to recycling baryonic material in the mid-plane of galaxies., Comment: Accepted for publication in Astronomy & Astrophysics; 16 pages, 8 figures

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)\lesssim11.8$) quiescent galaxies at $0.6\lesssim z\lesssim1.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., Comment: Accepted for publication in MNRAS; minor typesetting corrections after copyediting

We use deep spectroscopy from the SAMI Galaxy Survey to explore the precision of the fundamental plane of early-type galaxies (FP) as a distance indicator for future single-fibre spectroscopy surveys. We study the optimal trade-off between sample size and signal-to-noise ratio (SNR), and investigate which additional observables can be used to construct hyperplanes with smaller intrinsic scatter than the FP. We add increasing levels of random noise (parametrised as effective exposure time) to the SAMI spectra to study the effect of increasing measurement uncertainties on the FP-and hyperplane-inferred distances. We find that, using direct-fit methods, the values of the FP and hyperplane best-fit coefficients depend on the spectral SNR, and reach asymptotic values for a mean SNR=40 {\AA}$^{-1}$. As additional variables for the FP we consider three stellar-population observables: light-weighted age, stellar mass-to-light ratio and a novel combination of Lick indices (I$_{\rm age}$). For a SNR=45 {\AA}$^{-1}$ (equivalent to 1-hour exposure on a 4-m telescope), all three hyperplanes outperform the FP as distance indicators. Being an empirical spectral index, I$_{\rm age}$ avoids the model-dependent uncertainties and bias underlying age and mass-to-light ratio measurements, yet yields a 10 per cent reduction of the median distance uncertainty compared to the FP. We also find that, as a by-product, the Iage hyperplane removes most of the reported environment bias of the FP. After accounting for the different signal-to-noise ratio, these conclusions also apply to a 50 times larger sample from SDSS-III. However, in this case, only age removes the environment bias., Comment: 24 pages, 18 figures, accepted for publication in MNRAS

We present the galaxy stellar mass - size relation in the rest-frame near-IR ($1.5~\mu{\text{m}}$) and its evolution with redshift up to $z=2.5$. S\'ersic profiles are measured for $\sim$ $26\,000$ galaxies with stellar masses $M_\star > 10^9~{\text{M}}_\odot$ from JWST/NIRCam F277W and F444W imaging provided by the COSMOS-WEB and PRIMER surveys, using coordinates, redshifts, colors and stellar mass estimates from the COSMOS2020 catalog. The new rest-frame near-IR effective radii are generally smaller than previously measured rest-frame optical sizes, on average by 0.14~dex, with no significant dependence on redshift. For quiescent galaxies this size offset does not depend on stellar mass, but for star-forming galaxies the offset increases from -0.1~dex at $M_\star = 10^{9.5}~{\text{M}}_\odot$ to -0.25~dex at $M_\star > 10^{11}~{\text{M}}_\odot$. That is, we find that the near-IR stellar mass - size relation for star-forming galaxies is flatter in the rest-frame near-IR than in the rest-frame optical at all redshifts $0.510^{11}~{\text{M}}_\odot$) star-forming galaxies evolve in size almost as fast as quiescent galaxies. Low-mass ($M_\star<10^{10}~{\text{M}}_\odot$)~quiescent galaxies evolve as slow as star-forming galaxies. Our main conclusion is that the size evolution narrative as it has emerged over the past two decades does not radically change when accessing with JWST the rest-frame near-IR, a better proxy of the underlying stellar mass distribution., Comment: Accepted for publication. V2: includes a new table (table 1) of the data produced and available in machine-readable format on the journal version of the article

Using deep imaging from the CANDELS and HFF surveys, we present bulge+disc decompositions with GalfitM for $\sim$17,000 galaxies over $0.2 \leq z\leq 1.5$. We use various model parameters to select reliable samples of discs and bulges, and derive their stellar masses using an empirically calibrated relation between mass-to-light ratio and colour. Across our entire redshift range, we show that discs follow stellar mass-size relations that are consistent with those of star-forming galaxies, suggesting that discs primarily evolve via star formation. In contrast, the stellar mass-size relations of bulges are mass-independent. Our novel dataset further enables us to separate components into star-forming and quiescent based on their specific star formation rates. We find that both star-forming discs and star-forming bulges lie on stellar mass-size relations that are similar to those of star-forming galaxies, while quiescent discs are typically smaller than star-forming discs and lie on steeper relations, implying distinct evolutionary mechanisms. Similar to quiescent galaxies, quiescent bulges show a flattening in the stellar mass-size relation at $\sim$10$^{10}$M$_\odot$, below which they show little mass dependence. However, their best-fitting relations have lower normalisations, indicating that at a given mass, bulges are smaller than quiescent galaxies. Finally, we obtain rest-frame colours for individual components, showing that bulges typically have redder colours than discs, as expected. We visually derive UVJ criteria to separate star-forming and quiescent components and show that this separation agrees well with component colour. HFF bulge+disc decomposition catalogues used for these analyses are publicly released with this paper., Comment: 30 pages, 18 figures, and 6 tables. Resubmitted to MNRAS after addressing a thorough and constructive referee report

We present a study of the growth of the quiescent galaxy population between 0.5 < z < 3 by tracing the number density and structural evolution of a sample of 4518 old and 583 young quiescent galaxies with log($M_*$/$M_{\odot}$)>10.4, selected from the COSMOS2020 catalog with complementary HST/F160W imaging from the 3D-DASH survey. Among the quiescent population at z$\sim$2, roughly 50% are recently quenched galaxies; these young quiescent galaxies become increasingly rare towards lower redshift, supporting the idea that the peak epoch of massive galaxy quenching occurred at z>2. Our data show that while the effective half-light radii of quiescent galaxies generally increases with time, young quiescent galaxies are significantly smaller than their older counterparts at the same redshift. In this work we investigate the connection between this size difference and other structural properties, including axis ratios, color gradients, stellar mass, and the intrinsic scatter in effective radii. We demonstrate that the size difference is driven by the most massive sub-population (log($M_*$/$M_{\odot}$)>11) and does not persist when restricting the sample to intermediate mass galaxies (10.4

We apply the 3D dust radiative transfer code SKIRT to the low-redshift ($z\leq0.1$) galaxy population in the TNG100 cosmological simulation, the fiducial run of the IllustrisTNG project. We compute global fluxes and spectral energy distributions (SEDs) from the far-ultraviolet to the sub-millimeter for $\approx\,$60 000 galaxies. Our post-processing methodology follows the study of Tr\v{c}ka et al. (2022) of the higher-resolution TNG50 simulation. We verify that TNG100 reproduces observational luminosity functions at low redshifts to excellent precision, unlike TNG50. Additionally, we test the realism of our TNG100 plus SKIRT fluxes by comparing various flux and color relations to data from the GAMA survey. TNG100 broadly reproduces the observed distributions, but we predict ultraviolet colors that are too blue by $\approx\,$0.4 mag, possibly related to the extinction in the star-forming regions subgrid model not being selective enough. Furthermore, we find that the simulated galaxies exhibit mid-infrared fluxes elevated by up to $\approx\,$0.5 mag that we attribute to overly effective stochastic heating of the diffuse dust. All synthetic broadband fluxes and SEDs are made publicly available in three orientations and four apertures, and can readily be used to study TNG100 galaxies in a mock observational fashion., Comment: Main text 14 pages, 5 figures, accepted to MNRAS. Generated data products are publicly available at www.tng-project.org/gebek24

Galaxy sizes correlate with many other important properties of galaxies, and the cosmic evolution of galaxy sizes is an important observational diagnostic for constraining galaxy evolution models. The effective radius is probably the most widely used indicator of galaxy size. We used the TNG50-SKIRT Atlas to investigate the wavelength dependence of the effective radius of galaxies at optical and near-infrared (NIR) wavelengths. We find that, on average, the effective radius in every band exceeds the stellar mass effective radius, and that this excess systematically decreases with increasing wavelength. The optical g-band (NIR Ks-band) effective radius is on average 58% (13%) larger than the stellar mass effective radius. Effective radii measured from dust-obscured images are systematically larger than those measured from dust-free images, although the effect is limited (8.7% in the g-band, 2.1% in the Ks-band). We find that stellar population gradients are the dominant factor (about 80%) in driving the wavelength dependence of the effective radius, and that differential dust attenuation is a secondary factor (20%). Comparing our results to recent observational data, we find offsets in the absolute values of the median effective radii, up to 50% for the population of blue galaxies. We find better agreement in the slope of the wavelength dependence of the effective radius, with red galaxies having a slightly steeper slope than green-blue galaxies. Comparing our effective radii with those of galaxies from the Siena Galaxy Atlas in separate bins in z-band absolute magnitude and g-z colour, we find excellent agreement for the reddest galaxies, but again significant offsets for the blue populations: up to 70% for galaxies around Mz=-21.5. This difference in median effective radius for the bluer galaxies is most probably due to (abridged...)., Comment: 11 pages, 10 figures, accepted for publication in A&A

Galaxy morphology is a powerful diagnostic to assess the realism of cosmological hydrodynamical simulations. Determining the morphology of simulated galaxies requires the generation of synthetic images through 3D radiative transfer post-processing that properly accounts for different stellar populations and interstellar dust attenuation. We use the SKIRT code to generate the TNG50-SKIRT Atlas, a synthetic UV to near-infrared broadband image atlas for a complete stellar-mass selected sample of 1154 galaxies extracted from the TNG50 cosmological simulation at $z=0$. The images have a high spatial resolution (100 pc) and a wide field of view (160 kpc). In addition to the dust-obscured images, we also release dust-free images and physical parameter property maps with matching characteristics. As a sanity check and preview application we discuss the UVJ diagram of the galaxy sample. We investigate the effect of dust attenuation on the UVJ diagram and find that it affects both the star-forming and the quiescent galaxy populations. The quiescent galaxy region is polluted by younger and star-forming highly inclined galaxies, while dust attenuation induces a separation in inclination of the star-forming galaxy population, with low-inclination galaxies remaining at the blue side of the diagram and high-inclination galaxies systematically moving towards the red side. This image atlas can be used for a variety of other applications, including galaxy morphology studies and the investigation of local scaling relations. We publicly release the images and parameter maps, and we invite the community to use them., Comment: 13 pages, 7 figures, accepted for publication in A&A

Theories of spiral structure traditionally separate into tight-winding Lin-Shu spiral density waves and the swing-amplified material patterns of Goldreich & Lynden-Bell and Julian & Toomre. In this paper we consolidate these two types of spirals into a unified description, treating density waves beyond the tight-winding limit, in the regime of shearing and non-steady open spirals. This 'shearing wave' scenario novelly captures swing amplification that enables structure formation above conventional Q thresholds. However, it also highlights the fundamental role of spiral forcing on the amplification process in general, whether the wave is shearing or not. Thus it captures resonant and non-resonant mode growth through the donkey effect described by Lynden-Bell & Kalnajs and, critically, the cessation of growth when donkey behavior is no longer permitted. Our calculations predict growth exclusive to trailing spirals above the Jeans length, the prominence of spirals across a range of orientations that increases with decreasing arm multiplicity, and a critical orientation where growth is fastest that is the same for both modes and material patterns. Predicted structures are consistent with highly regular, high-multiplicity gaseous spur features and long filaments spaced close to the Jeans scale in spirals and bars. Applied to stellar disks, conditions favor low multiplicity (m<5) open trailing spirals with pitch angles in the observed range $10 deg$<$i_p$<$50 deg$. The results of this work serve as a basis for describing spirals as a unified class of transient waves, abundantly stimulated but narrowly selected for growth depending on local conditions., Comment: Accepted for publication in ApJ, 30 pages, 4 figures

We test whether we can predict optical spectra from deep-field photometry of distant galaxies. Our goal is to perform a comparison in data space, highlighting the differences between predicted and observed spectra. The Large Early Galaxy Astrophysics Census (LEGA-C) provides high-quality optical spectra of thousands of galaxies at redshift $0.6

We present and publicly release a new star-forming regions emission library TODDLERS (Time evolution of Observables including Dust Diagnostics and Line Emission from Regions containing young Stars) for the publicly available radiative transfer code SKIRT. The library generation involves the spherical evolution of a homogeneous gas cloud around a young stellar cluster that accounts for stellar feedback processes including stellar winds, supernovae, and radiation pressure, as well as the gravitational forces on the gas. The semi-analytical evolution model is coupled with the photoionization code Cloudy to calculate time-dependent UV-mm spectral energy distributions (SEDs) from star-forming regions of varying metallicity, star-formation efficiency, birth-cloud density, and mass. The calculated SEDs include the stellar, nebular, and dust continuum emission along with a wide range of emission lines originating from H ii, photo-dissociation, and molecular gas regimes tabulated at high resolution. The SEDs incorporated in SKIRT are generated by calculating a stellar-mass normalized luminosity, which assumes that each emission source is composed of a power-law population of star-forming clouds. When compared to the previous treatment of star-forming regions in SKIRT, TODDLERS shows a better agreement with low-redshift observational data in the IR wavelength range while offering a more comprehensive line-emission support. This paves the way for a variety of applications using simulated galaxies at low and high redshift., Comment: Accepted for publication in MNRAS on 27/09/23

We examine the wavelength dependence of radial light profiles based on S\'ersic index $n$ measurements of 1067 galaxies with M$_*\geq$ 10$^{9.5}$M$_\odot$ and in the redshift range $0.5 < z < 3$. The sample and rest-frame optical light profiles are drawn from CANDELS$+$3D-HST; rest-frame near-infrared light profiles are inferred from CEERS JWST/NIRCam imaging. $n$ shows only weak dependence on wavelength, regardless of redshift, galaxy mass and type: on average, star-forming galaxies have $n = 1-1.5$ and quiescent galaxies have $n = 3-4$ in the rest-frame optical and near-infrared. The strong correlation at all wavelengths between $n$ and star-formation activity implies a physical connection between the radial stellar mass profile and star-formation activity. The main caveat is that the current sample is too small to discern trends for the most massive galaxies (M$_* > 10^{11}M_\odot$)., Comment: 15 pages, 7 figures, Submitted to ApJ. Comments welcome

We present a comprehensive multi-wavelength study of 5 poststarburst galaxies with $M_\ast > 10^{11} M_\odot$ at $z\sim 0.7$, examining their stars, gas, and current and past star-formation activities. Using optical images from the Subaru telescope and Hubble Space Telescope, we observe a high incidence of companion galaxies and low surface brightness tidal features, indicating that quenching is closely related to interactions between galaxies. From optical spectra provided by the LEGA-C survey, we model the stellar continuum to derive the star-formation histories and show that the stellar masses of progenitors ranging from $2\times10^9 M_\odot$ to $10^{11} M_\odot$, undergoing a burst of star formation several hundred million years prior to observation, with a decay time scale of $\sim100$ million years. Our ALMA observations detect CO(2-1) emission in four galaxies, with the molecular gas spreading over up to $>1"$, or $\sim10$ kpc, with a mass of up to $\sim2 \times10^{10} M_\odot$. However, star-forming regions are unresolved by either the slit spectra or 3~GHz continuum observed by the Very Large Array. Comparisons between the star-formation rates and gas masses, and the sizes of CO emission and star-forming regions suggest a low star-forming efficiency. We show that the star-formation rates derived from IR and radio luminosities with commonly-used calibrations tend to overestimate the true values because of the prodigious amount of radiation from old stars and the contribution from AGN, as the optical spectra reveal weak AGN-driven outflows., Comment: Accepted by ApJ

We present individual star-formation histories of $\sim3000$ massive galaxies (log($\mathrm{M_*/M_{\odot}}$) > 10.5) from the Large Early Galaxy Astrophysics Census (LEGA-C) spectroscopic survey at a lookback time of $\sim$7 billion years and quantify the population trends leveraging 20hr-deep integrated spectra of these $\sim$ 1800 star-forming and $\sim$ 1200 quiescent galaxies at 0.6 < $z$ < 1.0. Essentially all galaxies at this epoch contain stars of age < 3 Gyr, in contrast with older massive galaxies today, facilitating better recovery of previous generations of star formation at cosmic noon and earlier. We conduct spectro-photometric analysis using parametric and non-parametric Bayesian SPS modeling tools - Bagpipes and Prospector to constrain the median star-formation histories of this mass-complete sample and characterize population trends. A consistent picture arises for the late-time stellar mass growth when quantified as $t_{50}$ and $t_{90}$, corresponding to the age of the universe when galaxies formed 50\% and 90\% of their total stellar mass, although the two sets of models disagree at the earliest formation times (e.g. $t_{10}$). Our results reveal trends in both stellar mass and stellar velocity dispersion as in the local universe - low-mass galaxies with shallower potential wells grow their stellar masses later in cosmic history compared to high-mass galaxies. Unlike local quiescent galaxies, the median duration of late-time star-formation ($\tau_{SF,late}$ = $t_{90}$ - $t_{50}$) does not consistently depend on the stellar mass. This census sets a benchmark for future deep spectro-photometric studies of the more distant universe., Comment: 28 pages, 12 figures, 4 tables, Accepted to ApJ

We use CEERS JWST/NIRCam imaging to measure rest-frame near-IR light profiles of $>$500 $M_\star>10^{10}~M_\odot$ galaxies in the redshift range $0.510^{11}~M_\odot$) quiescent galaxies between $z=0.5$ and $z=2.3$, again comparing $R_{\rm{opt}}$, \rmass~and \rmassd. We conclude that the main tenets of the size evolution narrative established over the past 20 years, based on rest-frame optical light profile analysis, still hold in the era of JWST/NIRCam observations in the rest-frame near-IR., Comment: Submitted to ApJ. Comments welcome

We analyse the joint distribution of dust attenuation and projected axis ratios, together with galaxy size and surface brightness profile information, to infer lessons on the dust content and star/dust geometry within star-forming galaxies at 0 < z <2.5. To do so, we make use of large observational datasets from KiDS+VIKING+HSC-SSP and extend the analysis out to redshift z = 2.5 using the HST surveys CANDELS and 3D-DASH. We construct suites of SKIRT radiative transfer models for idealized galaxies observed under random viewing angles with the aim of reproducing the aforementioned distributions, including the level and inclination dependence of dust attenuation. We find that attenuation-based dust mass estimates are at odds with constraints from far-infrared observations, especially at higher redshifts, when assuming smooth star and dust geometries of equal extent. We demonstrate that UV-to-near-IR and far-infrared constraints can be reconciled by invoking clumpier dust geometries for galaxies at higher redshifts and/or very compact dust cores. We discuss implications for the significant wavelength- and redshift-dependent differences between half-light and half-mass radii that result from spatially varying dust columns within -- especially massive -- star-forming galaxies., Comment: Accepted for publication in MNRAS

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 SAMI Galaxy survey. We perform spectral fitting on our synthetic cubes and compute the resulting rotation velocity ($V_{\rm{rot}}$) and velocity dispersion within 1$R_{\text{e}}$ ($\sigma_{\text{e}}$) for a sub-set 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/\sigma$ by nearly a factor two (after inclination correction). We also demonstrate the utility of our method for high-redshift galaxies by emulating spatially resolved continuum spectra from the LEGA-C survey and, looking forward, E-ELT 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., Comment: 17 pages, 11 figures. Accepted for publication in MNRAS

The massive end of the gas-phase mass--metallicity relation (MZR) is a sensitive probe of active galactic nuclei (AGN) feedback that is a crucial but highly uncertain component of galaxy evolution models. In this paper, we extend the $z\sim0.7$ MZR by $\sim$0.5 dex up to log$(M_\star/\textrm{M}_\odot)\sim11.1$. We use extremely deep VLT VIMOS spectra from the Large Early Galaxy Astrophysics Census (LEGA-C) survey to measure metallicities for 145 galaxies. The LEGA-C MZR matches the normalization of the $z\sim0.8$ DEEP2 MZR where they overlap, so we combine the two to create an MZR spanning from 9.3 to 11.1 log$(M_\star/\textrm{M}_\odot)$. The LEGA-C+DEEP2 MZR at $z\sim0.7$ is offset to slightly lower metallicities (0.05-0.13 dex) than the $z\sim0$ MZR, but it otherwise mirrors the established power law rise at low/intermediate stellar masses and asymptotic flattening at high stellar masses. We compare the LEGA-C+DEEP2 MZR to the MZR from two cosmological simulations (IllustrisTNG and SIMBA), which predict qualitatively different metallicity trends for high-mass galaxies. This comparison highlights that our extended MZR provides a crucial observational constraint for galaxy evolution models in a mass regime where the MZR is very sensitive to choices about the implementation of AGN feedback., Comment: 10 pages, 4 figures, 1 table

The radio jets of an active galactic nucleus (AGN) can heat up the gas around a host galaxy and quench star formation activity. The presence of a radio jet could be related to the evolutionary path of the host galaxy and may be imprinted in the morphology and kinematics of the galaxy. In this work, we use data from the Sloan Digital Sky Survey's Mapping Nearby Galaxies at Apache Point Observatory survey and the Low Frequency Array (LOFAR) Two-Metre Sky Survey as well as the National Radio Astronomy Observatory (NRAO) the Karl G. Jansky Very Large Array (VLA) Sky Survey and the Faint Images of the Radio Sky at Twenty Centimeter survey. We combine these integral field spectroscopic data and radio data to study the link between stellar kinematics and radio AGNs. We find that the luminosity-weighted stellar angular momentum $\lambda_{Re}$ is tightly related to the range of radio luminosity and the fraction of radio AGNs F radio present in galaxies, as high-luminosity radio AGNs are only in galaxies with a small $\lambda_{Re}$, and the $F_{radio}$ at a fixed stellar mass decreases with $\lambda_{Re}$. These results indicate that galaxies with stronger random stellar motions with respect to the ordered motions might be better breeding grounds for powerful radio AGNs. This would also imply that the merger events of galaxies are important in the triggering of powerful radio jets in our sample., Comment: 10 pages,7 figures. Accepted in A&A

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 parametrised as a Gauss-Hermite series. Conservatively using a redshift-independent cut in stellar mass ($M_\star = 10^{11}\,{\rm M}_\odot$), and matching the stellar-mass distributions of our samples, we find 7 $\sigma$ evidence of $h_4$ increasing with cosmic time, from a median value of 0.019$\pm$0.002 at z=0.8 to 0.059$\pm$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 $\sigma$ 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., Comment: 19 pages, 9 figures Accepted for publication in MNRAS

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 > 10^{10.5} {\rm 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 parametrised as a Gauss-Hermite series. We find that at all redshifts where we have large enough samples, $h_4$ anti-correlates with the ratio between rotation and dispersion, highlighting the physical connection between these two kinematic observables. In addition, and independently from the anti-correlation with rotation-to-dispersion ratio, we also find a correlation between $h_4$ and $M_\star$, 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., Comment: 26 pages, 15 figures Accepted for publication in MNRAS

We present elemental abundance patterns (C, N, Mg, Si, Ca, Ti, V, Cr, Fe, Co, and Ni) for a population of 135 massive quiescent galaxies at $z\sim0.7$ with ultra-deep rest-frame optical spectroscopy drawn from the LEGA-C survey. We derive average ages and elemental abundances in four bins of stellar velocity dispersion ($\sigma_v$) ranging from 150$~$km$\,$s$^{-1}$ to 250$~$km$\,$s$^{-1}$ using a full-spectrum hierarchical Bayesian model. The resulting elemental abundance measurements are precise to 0.05$\,$dex. The majority of elements, as well as the total metallicity and stellar age, show a positive correlation with $\sigma_v$. Thus, the highest dispersion galaxies formed the earliest and are the most metal-rich. We find only mild or non-significant trends between [X/Fe] and $\sigma_v$, suggesting that the average star-formation timescale does not strongly depend on velocity dispersion. To first order, the abundance patterns of the $z\sim0.7$ quiescent galaxies are strikingly similar to those at $z\sim0$. However, at the lowest velocity dispersions the $z\sim0.7$ galaxies have slightly enhanced N, Mg, Ti, and Ni abundance ratios and earlier formation redshifts than their $z\sim0$ counterparts. Thus, while the higher-mass quiescent galaxy population shows little evolution, the low-mass quiescent galaxies population has grown significantly over the past six billion years. Finally, the abundance patterns of both $z\sim0$ and $z\sim0.7$ quiescent galaxies differ considerably from theoretical prediction based on a chemical evolution model, indicating that our understanding of the enrichment histories of these galaxies is still very limited., Comment: 24 pages, 12 figures, accepted for publication in ApJ

We explore the properties of the host galaxies of X-ray selected AGN in the COSMOS field using the Chandra Legacy sample and the LEGA-C survey VLT optical spectra. Our main goal is to compare the relative ages of the host galaxies of the obscured and unobscured AGN by means of the calcium break Dn(4000) and the Hdelta Balmer line. The host galaxy ages are examined in conjunction with other properties such as the galaxy stellar mass, and star-formation rate as well as the AGN Eddington ratio. Our sample consists of 50 unobscured or mildly obscured (logN_H (cm-2) < 23) and 23 heavily obscured AGN (log N_H (cm-2) >23) in the redshift range z=0.6-1. We take specific caution to create control samples in order to match the exact luminosity and redshift distributions for the obscured and unobscured AGN. The majority of unobscured AGN appear to live in young galaxies in contrast to the obscured AGN which appear to live in galaxies located between the young and old galaxy populations. This finding may be in contrast to those evolutionary AGN unification models which postulate that the AGN begin their life in a heavy obscuration phase. The host galaxies of the obscured AGN have significantly lower levels of specific star-formation. At the same time the obscured AGN have lower Eddington ratios indicating a link between the star-formation and the black hole accretion. We find that the distribution of the stellar masses of the host galaxies of obscured AGN is skewed towards higher stellar masses in agreement with previous findings. Our results on the relative age of obscured AGN are valid when we match our obscured and unobscured AGN samples according to the stellar mass of their host galaxies. All the above results become less conspicuous when a lower column density (log N_H(cm-2)= 21.5 or 22) is used to separate the obscured and unobscured AGN populations., Comment: Accepted for publication in A&A

JWST/MIRI imaging of the nearby galaxies IC 5332, NGC 628, NGC 1365 and NGC 7496 from PHANGS reveals a richness of gas structures that in each case form a quasi-regular network of interconnected filaments, shells and voids. We examine whether this multi-scale network of structure is consistent with the fragmentation of the gas disk through gravitational instability. We use FilFinder to detect the web of filamentary features in each galaxy and determine their characteristic radial and azimuthal spacings. These spacings are then compared to estimates of the most Toomre-unstable length (a few kpc), the turbulent Jeans length (a few hundred pc) and the disk scale height (tens of pc) reconstructed using PHANGS-ALMA observations of the molecular gas as a dynamical tracer. Our analysis of the four galaxies targeted in this work indicates that Jeans-scale structure is pervasive. Future work will be essential for determining how the structure observed in gas disks impacts not only the rate and location of star formation but also how stellar feedback interacts positively or negatively with the surrounding multi-phase gas reservoir., Comment: 15 pages, 4 figures. Accepted as part of a PHANGS-JWST ApJL Focus Issue

We present a joint analysis of the galaxy S04590 at $z=8.496$ based on NIRSpec, NIRCam, and NIRISS observations obtained through as part of Early Release Observations programme of the James Webb Space Telescope (JWST) and the far-infrared [CII]-$158\mu$m emission line detected by dedicated Atacama Large Millimeter/submillimeter Array (ALMA) observations. We determine the physical properties of S04590 from modelling of the spectral energy distribution (SED) and through the redshifted optical nebular emission lines detected with JWST/NIRSpec. The best-fit SED model reveals a low-mass ($M_\star = 10^{7.2}-10^{8}\,M_{\odot}$) galaxy with a low oxygen abundance of $12+\log{\rm (O/H)} = 7.16^{+0.10}_{-0.12}$ derived from the strong nebular and auroral emission lines. Assuming that [CII] effectively traces the interstellar medium (ISM), we estimate the total gas mass of the galaxy to be $M_{\rm gas} = (8.0\pm 4.0)\times 10^{8}\,M_\odot$ based on the luminosity and spatial extent of [CII]. This yields an exceptionally high gas fraction, $f_{\rm gas} = M_{\rm gas}/(M_{\rm gas} + M_\star) \gtrsim 90\%$, though still consistent within the range expected for its low metallicity. We further derive the metal mass of the galaxy based on the gas mass and gas-phase metallicity, which we find to be consistent with the expected metal production from Type II supernovae. Finally, we make the first constraints on the dust-to-gas (DTG) and dust-to-metals (DTM) ratios of galaxies in the epoch of reionization at $z\gtrsim 6$, showing overall low mass ratios of logDGT $<-3.8$ and logDTM $<-0.5$, though consistent with local scaling relations and in particular the local metal-poor galaxy I Zwicky 18. Our analysis highlights the synergy between ALMA and JWST in characterizing the gas, metal, and stellar content of the first generation of galaxies., Comment: Accepted for publication in ApJL

Over the last two decades, around 300 quasars have been discovered at $z\gtrsim6$, yet only one has identified as being strongly gravitationally lensed. We explore a new approach -- enlarging the permitted spectral parameter space, while introducing a new spatial geometry veto criterion -- which is implemented via image-based deep learning. We first apply this approach to a systematic search for reionization-era lensed quasars, using data from the Dark Energy Survey, the Visible and Infrared Survey Telescope for Astronomy Hemisphere Survey, and the Wide-field Infrared Survey Explorer.Our search method consists of two main parts: (i) the preselection of the candidates based on their spectral energy distributions (SEDs) using catalog-level photometry and (ii) relative probabilities calculation of the candidates being a lens or some contaminant, utilizing a convolutional neural network (CNN) classification. The training data sets are constructed by painting deflected point-source lights over actual galaxy images, to generate realistic galaxy-quasar lens models, optimized to find systems with small image separations, i.e., Einstein radii of $\theta_\mathrm{E} \leq 1$ arcsec. Visual inspection is then performed for sources with CNN scores of $P_\mathrm{lens} > 0.1$, which leads us to obtain 36 newly selected lens candidates, which are awaiting spectroscopic confirmation. These findings show that automated SED modeling and deep learning pipelines, supported by modest human input, are a promising route for detecting strong lenses from large catalogs that can overcome the veto limitations of primarily dropout-based SED selection approaches., Comment: 24 pages, 17 figures, and 2 tables. Accepted for publication in The Astrophysical Journal. We welcome comments from the reader

Rejuvenating galaxies are unusual galaxies that fully quench and then subsequently experience a "rejuvenation" event to become star-forming once more. Rejuvenation rates vary substantially in models of galaxy formation: 10%-70% of massive galaxies are expected to experience rejuvenation by z = 0. Measuring the rate of rejuvenation is therefore important for calibrating the strength of star formation feedback mechanisms. However, these observations are challenging because rejuvenating systems blend in with normal star-forming galaxies in broadband photometry. In this paper, we use the galaxy spectral energy distribution (SED)-fitting code Prospector to search for observational markers that distinguish normal star-forming galaxies from rejuvenating galaxies. We find that rejuvenating galaxies have smaller Balmer absorption line equivalent widths (EWs) than star-forming galaxies. This is analogous to the well-known "K + A" or post-starburst galaxies, which have strong Balmer absorption due to A-stars dominating the light: in this case, rejuvenating systems have a lack of A-stars, instead resembling "O - A" systems. We find star-forming galaxies that have H$\beta$, H$\gamma$, and/or H$\delta$ absorption EWs $\lesssim 3${\AA} corresponds to a highly pure selection of rejuvenating systems. Interestingly, while this technique is highly effective at identifying mild rejuvenation, "strongly" rejuvenating systems remain nearly indistinguishable from star-forming galaxies due to the well-known stellar outshining effect. We conclude that measuring Balmer absorption line EWs in star-forming galaxy populations is an efficient method to identify rejuvenating populations, and discuss several techniques to either remove or resolve the nebular emission which typically lies on top of these absorption lines., Comment: Published in ApJ, 16 pages, 9 figures

We compare ionised gas and stellar kinematics of 16 star-forming galaxies ($\log(M_\star/M_\odot)=9.7-11.2$, SFR=6-86 $M_\odot/yr$) at $z\sim1$ using near-infrared integral field spectroscopy (IFS) of H$\alpha$ emission from the KMOS$^{\rm 3D}$ survey and optical slit spectroscopy of stellar absorption and gas emission from the LEGA-C survey. H$\alpha$ 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$\alpha$ relative to [OII], 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$\alpha$ dynamical mass models agree remarkably well out to ~10 kpc for all but one galaxy (average $\Delta M_{\rm dyn}(R_{e,\rm F814W})<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$\alpha$ kinematic major axis and the photometric position angle, highlighting the added value of IFS observations for dynamics studies. The good agreement between the JAM models and the dynamical models based on H$\alpha$ kinematics at $z\sim1$ corroborates the validity of dynamical mass measurements from H$\alpha$ IFS observations also for higher redshift rotating disc galaxies., Comment: 29 pages, 16 figures, 4 tables; accepted for publication in MNRAS; doi:10.1093/mnras/stad3826

We present the first direct spectroscopic measurement of the stellar velocity dispersion function (VDF) for massive quiescent and star-forming galaxies at $0.6 < z \leq 1.0$. For this analysis we use individual measurements of stellar velocity dispersion from high-S/N spectra from the public Large Early Galaxy Astrophysics Census (LEGA-C) survey. We report a remarkable stability of the VDF for both quiescent and star-forming galaxies within this redshift range, though we note the presence of weak evolution in the number densities of star-forming galaxies. We compare both VDFs with previous direct and inferred measurements at local and intermediate redshifts, with the caveat that previous measurements of the VDF for star-forming galaxies are poorly constrained at all epochs. We emphasize that this work is the first to directly push to low-stellar velocity dispersion ($\sigma_\star > 100$ km s$^{-1}$) and extend to star-forming galaxies. We are largely consistent with the high-sigma tail measured from BOSS, and we find that the VDF remains constant from the median redshift of LEGA-C, $z\sim0.8$, to the present day., Comment: 11 pages, 3 figures, accepted for publication in ApJ

Dynamical models for $673$ galaxies at $z=0.6-1.0$ with spatially resolved (long-slit) stellar kinematic data from LEGA-C are used to calibrate virial mass estimates defined as $M_{\rm{vir}}=K \sigma'^2_{\star,\rm{int}} R$, with $K$ a scaling factor, $\sigma'_{\star,\rm{int}}$ the spatially-integrated stellar velocity second moment from the LEGA-C survey and $R$ the effective radius measured from a S\'ersic profile fit to HST imaging. The sample is representative for $M_{\star}>3\times10^{10}~M_{\odot}$ and includes all types of galaxies, irrespective of morphology and color. We demonstrate that using $R=R_{\rm{sma}}$~(the semi-major axis length of the ellipse that encloses 50\% of the light) in combination with an inclination correction on $\sigma'_{\star,\rm{int}}$~produces an unbiased $M_{\rm{vir}}$. We confirm the importance of projection effects on $\sigma'_{\star,\rm{int}}$ by showing the existence of a similar residual trend between virial mass estimates and inclination for the nearby early-type galaxies in the ATLAS$^{\rm{3D}}$~survey. Also, as previously shown, when using a S\'ersic profile-based $R$ estimate, then a S\'{e}rsic index-dependent correction to account for non-homology in the radial profiles is required. With respect to analogous dynamical models for low-redshift galaxies from the ATLAS$^{\rm{3D}}$~survey we find a systematic offset of 0.1 dex in the calibrated virial constant for LEGA-C, which may be due to physical differences between the galaxy samples or an unknown systematic error. Either way, with our work we establish a common mass scale for galaxies across 8 Gyr of cosmic time with a systematic uncertainty of at most 0.1 dex., Comment: Published in ApJ

$JWST$'s first glimpse of the $z>10$ Universe has yielded a surprising abundance of luminous galaxy candidates. Here we present the most extreme of these systems: CEERS-1749. Based on $0.6-5\mu$m photometry, this strikingly luminous ($\approx$26 mag) galaxy appears to lie at $z\approx17$. This would make it an $M_{\rm{UV}}\approx-22$, $M_{\rm{\star}}\approx5\times10^{9}M_{\rm{\odot}}$ system that formed a mere $\sim220$ Myrs after the Big Bang. The implied number density of this galaxy and its analogues challenges virtually every early galaxy evolution model that assumes $\Lambda$CDM cosmology. However, there is strong environmental evidence supporting a secondary redshift solution of $z\approx5$: all three of the galaxy's nearest neighbors at $<2.5$" have photometric redshifts of $z\approx5$. Further, we show that CEERS-1749 may lie in a $z\approx5$ protocluster that is $\gtrsim5\times$ overdense compared to the field. Intense line emission at $z\approx5$ from a quiescent galaxy harboring ionized gas, or from a dusty starburst, may provide satisfactory explanations for CEERS-1749's photometry. The emission lines at $z\approx5$ conspire to boost the $>2\mu$m photometry, producing an apparent blue slope as well as a strong break in the SED. Such a perfectly disguised contaminant is possible only in a narrow redshift window ($\Delta z\lesssim0.1$), implying that the permitted volume for such interlopers may not be a major concern for $z>10$ searches, particularly when medium-bands are deployed. If CEERS-1749 is confirmed to lie at $z\approx5$, it will be the highest-redshift quiescent galaxy, or one of the lowest mass dusty galaxies of the early Universe detected to-date. Both redshift solutions of this intriguing galaxy hold the potential to challenge existing models of early galaxy evolution, making spectroscopic follow-up of this source critical., Comment: Submitted to ApJL. Fig. 1 shows the "superposed" state, Fig. 2 presents the potential z~5 protocluster, Fig. 3 shows some atypical z~5 SEDs, Fig. 4 studies the morphology, Fig. 5 summarizes the extraordinary situation at z~17, and Fig. 6 argues for why such interlopers may not be a major issue for z>10 searches. Comments warmly welcomed

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 ($<10\%$ 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., Comment: accepted for publication in MNRAS; resolution of some figures reduced for arxiv version

Links between the properties of radio-loud active galactic nuclei (RLAGNs) and the morphology of their hosts may provide important clues for our understanding of how RLAGNs are triggered. In this work, focusing on passive galaxies, we study the shape of the hosts of RLAGNs selected from the Karl G. Jansky Very Large Array Cosmic Evolution Survey (VLA-COSMOS) 3GHz Large Project, and compare them with previous results based on the first data release (DR1) of the LOFAR Two-Metre Sky Survey (LoTSS). We find that, at redshifts of between 0.6 and 1, high-luminosity ($L_{1.4 GHz}\gtrsim10^{24}\rm W Hz^{-1}$) RLAGNs have a wider range of optical projected axis ratios than their low-redshift counterparts, which are essentially all found in round galaxies with axis ratios of higher than 0.7. We construct control samples and show that although the hosts of high-redshift RLAGNs with the highest luminosities still have a rounder shape compared with the non-RLAGNs, they on average have a smaller axis ratio (more elongated) than the local RLAGNs with similar stellar masses and radio luminosities. This evolution can be interpreted as a byproduct of radio luminosity evolution, namely that galaxies at fixed stellar mass are more radio luminous at high redshifts: artificially increasing the radio luminosities of local galaxies ($z\leq$0.3) by a factor of 2 to 4 can remove the observed evolution of the axis ratio distribution. If this interpretation is correct then the implication is that the link between AGN radio luminosity and host galaxy shape is similar at $z\simeq1$ to in the present-day Universe., Comment: 10 pages, 6 figures. Accepted by A&A