Search

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 present the first results on the spatial distribution of star formation in 454 star-forming galaxies at $4.8

We present follow-up spectroscopy on a protocluster candidate selected from the wide-field imaging of the Hyper SuprimeCam Subaru Strategic Programme. The target protocluster candidate was identified as a $4.5\sigma$ overdense region of $g$-dropout galaxies, and the redshifts of $g$-dropout galaxies are determined by detecting their Ly$\alpha$ emission. Thirteen galaxies, at least, are found to be clustering in the narrow redshift range of $\Delta z<0.05$ at $z=3.699$. This is clear evidence of the presence of a protocluster in the target region. Following the discovery of the protocluster at $z=3.699$, the physical properties and three-dimensional distribution of its member galaxies are investigated. Based on spectroscopically-confirmed $g$-dropout galaxies, we find an overabundance of rest-frame ultraviolet (UV) bright galaxies in the protocluster. The UV brightest protocluster member turns out to be an active galactic nucleus, and the other UV brighter members tend to show smaller Ly$\alpha$ equivalent widths than field counterparts. The member galaxies tend to densely populate near the centre of the protocluster, but the separation from the nearest neighbour rather than the distance from the centre of the protocluster is more tightly correlated to galaxy properties, implying that the protocluster is still in an early phase of cluster formation and only close neighbours have a significant impact on the physical properties of protocluster members. The number density of massive galaxies, selected from an archival photometric-redshift catalogue, is higher near the centre of the protocluster, while dusty starburst galaxies are distributed on the outskirts. These observational results suggest that the protocluster consists of multiple galaxy populations, whose spatial distributions may hint at the developmental phase of the galaxy cluster., Comment: 11 pages, 8 figures, submitted to MNRAS

We present results from the JWST First Reionization Epoch Spectroscopically Complete Observations survey (FRESCO) on the star forming sequence of galaxies at $1.09.5 M_\odot$ that are lower than found in many earlier studies by up to 0.6 dex, but in good agreement with recent results obtained with the Prospector fitting framework. The difference log(SFR(Pa$\alpha$)-SFR(Prospector)) is -0.09 $\pm$ 0.04 dex at $10^{10-11} M_\odot$. We also measure the empirical relation between Paschen-$\alpha$ luminosity and rest-frame H band magnitude and find that the scatter is only 0.04 dex lower than that of the SFR-M* relation and is much lower than the systematic differences among relations in the literature due to various methods of converting observed measurements to physical properties. We additionally identify examples of sources -- that, with standard cutoffs via the UVJ diagram, would be deemed quiescent -- with significant, typically extended, Paschen-$\alpha$ emission. Our results may be indicative of the potential unification of methods used to derive the star forming sequence with careful selection of star forming galaxies and independent star formation rate and stellar mass indicators., Comment: 17 pages, 9 figures. Accepted for publication in ApJ

This paper utilizes the JWST MIRI multi-band imaging data from the SMILES survey (5-25micron), complemented with HST and NIRCam photometric and spectroscopic data from the JADES and FRESCO surveys for 443 star-forming (non-AGN) galaxies at z=0.7-2.0 to extend the study of dust and PAH emission to a new mass and SFR parameter space beyond our local universe. We find a strong correlation between the fraction of dust in PAHs (PAH fraction, q_PAH) with stellar mass. Moreover, the PAH fraction behavior as a function of gas-phase metallicity is similar to that at z~0 from previous studies, suggesting a universal relation: q_PAH is constant (~3.4%) above a metallicity of ~ 0.5$Z_{\odot}$ and decreases to <1% at metallicities $<0.3Z_{\odot}$. This indicates that metallicity is a good indicator of the ISM properties that affect the balance between the formation and destruction of PAHs. The lack of a redshift evolution from z~0-2 also implies that above $0.5\,Z_{\odot}$, the PAH emission effectively traces obscured luminosity and the previous locally-calibrated PAH-SFR calibrations remain applicable in this metallicity regime. We observe a strong correlation between obscured UV luminosity fraction (ratio of obscured to total luminosity) and stellar mass. Above the stellar mass of $>5\times 10^9M_{\odot}$, on average, more than half of the emitted luminosity is obscured, while there exists a non-negligible population of lower mass galaxies with >50% obscured fractions. At a fixed mass, the obscured fraction correlates with SFR surface density. This is a result of higher dust covering fractions in galaxies with more compact star forming regions. Similarly, galaxies with high IRX (IR to UV luminosity) at a given mass or UV continuum slope tend to have higher SFR surface density and shallower attenuation curves, owing to their higher effective dust optical depths and more compact star forming regions., Comment: A&A accepted. Supplementary material on https://zenodo.org/records/12671075

We present a 24 um-selected spectroscopic sample z > 0.13 (median z = 0.41) in the Lockman Hole field, consisting of 4035 spectra. Our aim is to identify AGNs and determine their fraction in this mid-infrared selected sample. In this work, we use the [Ne V]3426 emission line to spectroscopically identify AGNs. Combined with broad-line Type I AGNs selected in our previous study, our sample consists of 887 (22%) spectroscopically confirmed AGNs. We perform a stacking analysis on the remaining spectra, and find that in various MIR-wedge-selected AGN candidates, the stacked spectra still show significant [Ne V]3426 emission, In contrast, no clear [Ne V]3426 signal is detected in non-AGN candidates falling outside the wedges. Assuming a range of AGN mid-IR SED slope of 0.3< alpha <0.7, and an average star-forming relation derived from 65 star-forming templates, we develop a robust method to separate the AGN and star-forming contributions to the mid-IR SEDs using the rest-frame L12 /L1.6 vs L4.5 /L1.6 diagram. We separate the objects into bins of L12 , and find that AGN fraction increases with increasing L12. We also find that the stacked [Ne V]3426 strength scales with L12 . The pure AGN luminosity at 12 um exhibits a positive correlation with the star formation rates, indicating possible co-evolution and common gas supply between the AGN and their host galaxies. Varying population properties across the redshift range explored contribute to the observed correlation., Comment: 16 pages, 14 figures. Accepted for publication in ApJ

The majority of massive disk galaxies in the local Universe show a stellar barred structure in their central regions, including our Milky Way. Bars are supposed to develop in dynamically cold stellar disks at low redshift, as the strong gas turbulence typical of disk galaxies at high redshift suppresses or delays bar formation. Moreover, simulations predict bars to be almost absent beyond $z = 1.5$ in the progenitors of Milky Way-like galaxies. Here we report observations of ceers-2112, a barred spiral galaxy at redshift $z_{\rm phot} \sim 3$, which was already mature when the Universe was only 2 Gyr old. The stellar mass ($M_{\star} = 3.9 \times 10^9 M_{\odot}$) and barred morphology mean that ceers-2112 can be considered a progenitor of the Milky Way, in terms of both structure and mass-assembly history in the first 2 Gyr of the Universe, and was the closest in mass in the first 4 Gyr. We infer that baryons in galaxies could have already dominated over dark matter at $z \sim 3$, that high-redshift bars could form in approximately 400 Myr and that dynamically cold stellar disks could have been in place by redshift $z = 4-5$ (more than 12 Gyrs ago)., Comment: Published in Nature. Accepted on September 12

The relatively red wavelength range (4800-9300{\AA}) of the VLT Multi Unit Spectroscopic Explorer (MUSE) limits which metallicity diagnostics can be used; in particular excluding those requiring the [O ii]{\lambda}{\lambda}3726,29 doublet. We assess various strong line diagnostics by comparing to sulphur Te-based metallicity measurements for a sample of 671 HII regions from 36 nearby galaxies from the MUSE Atlas of Disks (MAD) survey. We find that the O3N2 and N2 diagnostics return a narrower range of metallicities which lie up to ~0.3 dex below Te-based measurements, with a clear dependence on both metallicity and ionisation parameter. The N2S2H{\alpha} diagnostic shows a near-linear relation with the Te-based metallicities, although with a systematic downward offset of ~0.2 dex, but no clear dependence on ionisation parameter. These results imply that the N2S2H{\alpha} diagnostic produces the most reliable results when studying the distribution of metals within galaxies with MUSE. On sub-HII region scales, the O3N2 and N2 diagnostics measure metallicity decreasing towards the centres of HII regions, contrary to expectations. The S-calibration and N2S2H{\alpha} diagnostics show no evidence of this, and show a positive relationship between ionisation parameter and metallicity at 12 + log(O/H)> 8.4, implying the relationship between ionisation parameter and metallicity differs on local and global scales. We also present HIIdentify, a python tool developed to identify HII regions within galaxies from H{\alpha} emission maps. All segmentation maps and measured emission line strengths for the 4408 HII regions identified within the MAD sample are available to download., Comment: Accepted for publication in MNRAS. 21 pages, 13 figures

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

We present a protocluster search covering $z\sim3$ to $z\sim5$ based on the combination of the Hyper SuprimeCam Subaru Strategic Programme and the CFHT Large Area $U$-band Deep Survey. We identify about 30 protocluster candidates per unit redshift over the $\sim25\,\mathrm{deg^2}$ area of the Deep/Ultra-Deep layer. Protocluster candidates are selected as regions with a significantly enhanced surface density of dropout galaxies. With this large sample, we characterise the properties of their individual member galaxies. We compare the number counts of dropout galaxies in protocluster candidates with that of coeval field galaxies. Rest-frame UV bright galaxies are over-abundant in protocluster candidates, a trend seen across the full redshift range studied. We do not find evidence for their spatial distribution within protocluster candidates to be distinct from their fainter counterparts, nor for their UV colour to be different from that of field galaxies with the same brightness. Cosmological simulations predict this bright-end excess, with the main cause being a richer population of massive galaxies, with only a minor contribution from an enhancement in star formation activity (and therefore UV emission) at fixed mass. $U$-to-$K$ SED modelling of our observed samples supports this interpretation. This environmental differentiation in number counts is already in place at $z\sim5$, with no significant redshift dependence over the range in lookback times probed. These observational results and model predictions suggest that the cosmic clock is ahead in high-density environments., Comment: 17 pages, 13 figures, accepted for publication in MNRAS

With the remarkable sensitivity and resolution of JWST in the infrared, measuring rest-optical kinematics of galaxies at $z>5$ has become possible for the first time. This study pilots a new method for measuring galaxy dynamics for highly multiplexed, unbiased samples by combining FRESCO NIRCam grism spectroscopy and JADES medium-band imaging. Here we present one of the first JWST kinematic measurements for a galaxy at $z>5$. We find a significant velocity gradient, which, if interpreted as rotation yields $V_{rot} = 240\pm50$km/s and we hence refer to this galaxy as Twister-z5. With a rest-frame optical effective radius of $r_e=2.25$kpc, the high rotation velocity in this galaxy is not due to a compact size as may be expected in the early universe but rather a high total mass, ${\rm log(M}_{dyn}/{\rm M}_\odot)=11.0\pm0.2$. This is a factor of roughly 4x higher than the stellar mass within the effective radius. We also observe that the radial H$\alpha$ equivalent width profile and the specific star formation rate map from resolved stellar population modeling is centrally depressed by a factor of $\sim1.5$ from the center to $r_e$. Combined with the morphology of the line-emitting gas in comparison to the continuum, this centrally suppressed star formation is consistent with a star-forming disk surrounding a bulge growing inside-out. While large, rapidly rotating disks are common to z~2, the existence of one after only 1Gyr of cosmic time, shown for the first time in ionized gas, adds to the growing evidence that some galaxies matured earlier than expected in the history of the universe., Comment: Fig. 3 shows the main result

We study 10 spectroscopically confirmed Ly$\alpha$ emitters (LAEs) at $z\approx3.1$ in the UDS field, covered by JWST/NIRCam in the PRIMER program. All LAEs are detected in all NIRCam bands from F090W to F444W, corresponding to restframe 2200\AA--1.2$\mathrm{\mu m}$. Based on morphological analysis of the F200W images, three out of the 10 targets are resolved into pair-like systems with separations of $<0.9''$, and another three show asymmetric structures. We then construct the spectral energy distributions (SEDs) of these LAEs. All sources, including the pairs, show similar SED shapes, with a prominent flux excess in the F200W band, corresponding to extremely strong [O III]+H$\beta$ emission lines (${\rm EW_{rest}}=740$--$6500\,$\AA). The median effective radii, stellar mass, and UV slope of our sample are 0.36$\,$kpc, $3.8\times10^7\,M_\odot$, and --2.48, respectively. The average burst age, estimated by stellar mass over star formation rate, is $<40\,$Myr. These measurements reveal an intriguing starbursting dwarf galaxy population lying off the extrapolations of the $z \sim 3$ scaling relations to the low-mass end: $\sim 0.7$ dex above the star-forming main sequence, $\sim 0.35$ dex below the mass--size relation, and bluer in the UV slope than typical high-z galaxies at similar UV luminosities. We speculate that these numbers may require a larger main sequence scatter or tail in the dwarf galaxy regime towards the starburst outliers., Comment: 17 pages, 6 figures, ApJ accepted

We report the detection of a 13$\sigma$ H$\alpha$ emission line from HDF850.1 at $z=5.188\pm0.001$ using the FRESCO NIRCam F444W grism observations. Detection of H$\alpha$ in HDF850.1 is noteworthy, given its high far-IR luminosity, substantial dust obscuration, and the historical challenges in deriving its redshift. HDF850.1 shows a clear detection in the F444W imaging data, distributed between a northern and southern component, mirroring that seen in [CII] from the Plateau de Bure Interferometer. Modeling the SED of each component separately, we find that the northern component has a higher mass, star formation rate (SFR), and dust extinction than the southern component. The observed H$\alpha$ emission appears to arise entirely from the less-obscured southern component and shows a similar $\Delta$v$\sim$+130 km/s velocity offset to that seen for [CII] relative to the source systemic redshift. Leveraging H$\alpha$-derived redshifts from FRESCO observations, we find that HDF850.1 is forming in one of the richest environments identified to date at $z>5$, with 100 $z=5.17-5.20$ galaxies distributed across 10 structures and a $\sim$(15 cMpc)$^3$ volume. Based on the evolution of analogous structures in cosmological simulations, the $z=5.17-5.20$ structures seem likely to collapse into a single $>$10$^{14}$ $M_{\odot}$ cluster by $z\sim0$. Comparing galaxy properties forming within this overdensity with those outside, we find the masses, SFRs, and $UV$ luminosities inside the overdensity to be clearly higher. The prominence of H$\alpha$ line emission from HDF850.1 and other known highly-obscured $z>5$ galaxies illustrates the potential of NIRCam-grism programs to map both the early build-up of IR-luminous galaxies and overdense structures., Comment: Submitted to MNRAS. 20 pages, 10 figures and 8 tables (including appendices)

Recent JWST observations have revealed an unexpected abundance of massive galaxy candidates in the early Universe, extending further in redshift and to lower luminosity than what had previously been found by sub-millimeter surveys. These JWST candidates have been interpreted as challenging the $\Lambda$CDM cosmology, but, so far, they have mostly relied only on rest-frame ultraviolet data and lacked spectroscopic confirmation of their redshifts. Here we report a systematic study of 36 massive dust-obscured galaxies with spectroscopic redshifts between $z_{\rm spec}=5-9$ from the JWST FRESCO survey. We find no tension with the $\Lambda$CDM model in our sample. However, three ultra-massive galaxies (log$M_{\star}/M_{\odot}$ $\gtrsim11.0$) require an exceptional fraction of 50% of baryons converted into stars -- two to three times higher than even the most efficient galaxies at later epochs. The contribution from an active nucleus is unlikely because of their extended emission. Ultra-massive galaxies account for as much as 17% of the total cosmic star formation rate density at $z\sim5-6$., Comment: Nature in press. Updated to the accepted version. 24 pages, 4 main figures, 7 supplementary figures, 3 supplementary tables

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

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

[Abridged] We combined HST images from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey with JWST images from the Cosmic Evolution Early Release Science (CEERS) survey to measure the stellar and dust-obscured star formation distributions of a mass-complete ($>10^{10}M_\odot$) sample of 69 star-forming galaxies (SFGs) at $0.175$; 35 galaxies). At lower $S/N$, simulations show that simultaneous measurements of the size and S\'ersic index become unreliable. We extended our study to fainter sources ($S/N>10$; 69 galaxies) by fixing their S\'ersic index to unity. The S\'ersic index of bright galaxies ($S/N>75$) has a median value of 0.7, which, together with their axis ratio distribution, suggests a disk-like morphology in the rest-MIR. Galaxies above the main sequence (MS; i.e., starbursts) have rest-MIR sizes that are a factor 2 smaller than their rest-optical sizes. The median rest-optical to rest-MIR size ratio of MS galaxies increases with stellar mass, from 1.1 at $10^{9.8}M_\odot$ to 1.6 at $10^{11}M_\odot$. This mass-dependent trend resembles the one found in the literature between the rest-optical and rest-near-infrared sizes of SFGs, suggesting that it is due to radial color gradients affecting rest-optical sizes and that the sizes of the stellar and star-forming components of SFGs are, on average, consistent at all masses. There is, however, a small population of SFGs (15%) with a compact star-forming component embedded in a larger stellar structure. This could be the missing link between galaxies with an extended stellar component and those with a compact stellar component, the so-called blue nuggets., Comment: Published in A&A. 18 pages, 14 figures

In recent years, observations have uncovered a population of massive galaxies that are invisible or very faint in deep optical/near-infrared (near-IR) surveys but brighter at longer wavelengths. However, the nature of these optically dark or faint galaxies (OFGs; one of several names given to these objects) is highly uncertain. In this work, we investigate the drivers of dust attenuation in the JWST era. In particular, we study the role of stellar mass, size, and orientation in obscuring star-forming galaxies (SFGs) at $3 < z < 7.5$, focusing on the question of why OFGs and similar galaxies are so faint at optical/near-IR wavelengths. We find that stellar mass is the primary proxy for dust attenuation, among the properties studied. Effective radius and axis ratio do not show a clear link with dust attenuation, with the effect of orientation being close to random. However, there is a subset of highly dust attenuated ($A_V > 1$, typically) SFGs, of which OFGs are a specific case. For this subset, we find that the key distinctive feature is their compact size (for massive systems with $\log (M_{*}/M_{\odot}) > 10$); OFGs exhibit a 30% smaller effective radius than the average SFG at the same stellar mass and redshift. On the contrary, OFGs do not exhibit a preference for low axis ratios (i.e., edge-on disks). The results in this work show that stellar mass is the primary proxy for dust attenuation and compact stellar light profiles behind the thick dust columns obscuring typical massive SFGs., Comment: Accepted by A&A. 13 pages, 9 figures

We present an investigation into the first 500 Myr of galaxy evolution from the Cosmic Evolution Early Release Science (CEERS) survey. CEERS, one of 13 JWST ERS programs, targets galaxy formation from z~0.5 to z>10 using several imaging and spectroscopic modes. We make use of the first epoch of CEERS NIRCam imaging, spanning 35.5 sq. arcmin, to search for candidate galaxies at z>9. Following a detailed data reduction process implementing several custom steps to produce high-quality reduced images, we perform multi-band photometry across seven NIRCam broad and medium-band (and six Hubble broadband) filters focusing on robust colors and accurate total fluxes. We measure photometric redshifts and devise a robust set of selection criteria to identify a sample of 26 galaxy candidates at z~9-16. These objects are compact with a median half-light radius of ~0.5 kpc. We present an early estimate of the z~11 rest-frame ultraviolet (UV) luminosity function, finding that the number density of galaxies at M_UV ~ -20 appears to evolve very little from z~9 to z~11. We also find that the abundance (surface density [arcmin^-2]) of our candidates exceeds nearly all theoretical predictions. We explore potential implications, including that at z>10 star formation may be dominated by top-heavy initial mass functions, which would result in an increased ratio of UV light per unit halo mass, though a complete lack of dust attenuation and/or changing star-formation physics may also play a role. While spectroscopic confirmation of these sources is urgently required, our results suggest that the deeper views to come with JWST should yield prolific samples of ultra-high-redshift galaxies with which to further explore these conclusions., Comment: Replaced with published version

We present the data release and data reduction process for the Epoch 1 NIRCam observations for the Cosmic Evolution Early Release Science Survey (CEERS). These data consist of NIRCam imaging in six broadband filters (F115W, F150W, F200W, F277W, F356W and F444W) and one medium band filter (F410M) over four pointings, obtained in parallel with primary CEERS MIRI observations (Yang et al. in prep). We reduced the NIRCam imaging with the JWST Calibration Pipeline, with custom modifications and reduction steps designed to address additional features and challenges with the data. Here we provide a detailed description of each step in our reduction and a discussion of future expected improvements. Our reduction process includes corrections for known pre-launch issues such as 1/f noise, as well as in-flight issues including snowballs, wisps, and astrometric alignment. Many of our custom reduction processes were first developed with pre-launch simulated NIRCam imaging over the full 10 CEERS NIRCam pointings. We present a description of the creation and reduction of this simulated dataset in the Appendix. We provide mosaics of the real images in a public release, as well as our reduction scripts with detailed explanations to allow users to reproduce our final data products. These represent one of the first official public datasets released from the Directors Discretionary Early Release Science (DD-ERS) program., Comment: 27 pages, 14 figures, submitted to ApJ. Accompanying CEERS public Data Release 0.5 available at ceers.github.io/releases.html

The new capabilities that JWST offers in the near- and mid-infrared (IR) are used to investigate in unprecedented detail the nature of optical/near-IR faint, mid-IR bright sources, HST-dark galaxies among them. We gather JWST data from the CEERS survey in the EGS, jointly with HST data, and analyze spatially resolved optical-to-mid-IR spectral energy distributions (SEDs) to estimate both photometric redshifts in 2 dimensions and stellar populations properties in a pixel-by-pixel basis. We select 138 galaxies with F150W-F356W>1.5 mag, F356W<27.5 mag. The nature of these sources is threefold: (1) 71% are dusty star-forming galaxies at 2100 Gyr^-1); (2) 18% are quiescent/dormant (i.e., subject to reignition and rejuvenation) galaxies at 3

We present a comprehensive analysis of the evolution of the morphological and structural properties of a large sample of galaxies at z=3-9 using early JWST CEERS NIRCam observations. Our sample consists of 850 galaxies at z>3 detected in both CANDELS HST imaging and JWST CEERS NIRCam images to enable a comparison of HST and JWST morphologies. Our team conducted a set of visual classifications, with each galaxy in the sample classified by three different individuals. We also measure quantitative morphologies using the publicly available codes across all seven NIRCam filters. Using these measurements, we present the fraction of galaxies of each morphological type as a function of redshift. Overall, we find that galaxies at z>3 have a wide diversity of morphologies. Galaxies with disks make up a total of 60\% of galaxies at z=3 and this fraction drops to ~30% at z=6-9, while galaxies with spheroids make up ~30-40% across the whole redshift range and pure spheroids with no evidence for disks or irregular features make up ~20%. The fraction of galaxies with irregular features is roughly constant at all redshifts (~40-50%), while those that are purely irregular increases from ~12% to ~20% at z>4.5. We note that these are apparent fractions as many selection effects impact the visibility of morphological features at high redshift. The distributions of S\'ersic index, size, and axis ratios show significant differences between the morphological groups. Spheroid Only galaxies have a higher S\'ersic index, smaller size, and higher axis ratio than Disk/Irregular galaxies. Across all redshifts, smaller spheroid and disk galaxies tend to be rounder. Overall, these trends suggest that galaxies with established disks and spheroids exist across the full redshift range of this study and further work with large samples at higher redshift is needed to quantify when these features first formed., Comment: Accepted for publication to ApJL, 24 pages, 14 figures

Stellar bars are key drivers of secular evolution in galaxies and can be effectively studied using rest-frame near-infrared (NIR) images, which trace the underlying stellar mass and are less impacted by dust and star formation than rest-frame UV or optical images. We leverage the power of {\it{JWST}} CEERS NIRCam images to present the first quantitative identification and characterization of stellar bars at $z>1$ based on rest-frame NIR F444W images of high resolution (~1.3 kpc at z ~ 1-3). We identify stellar bars in these images using quantitative criteria based on ellipse fits. For this pilot study, we present six examples of robustly identified bars at $z>1$ with spectroscopic redshifts, including the two highest redshift bars at ~2.136 and 2.312 quantitatively identified and characterized to date. The stellar bars at $z$ ~ 1.1-2.3 presented in our study have projected semi-major axes of ~2.9-4.3 kpc and projected ellipticities of ~0.41-0.53 in the rest-frame NIR. The barred host galaxies have stellar masses ~ $ 1 \times 10^{10}$ to $2 \times 10^{11}$ $M_{\odot}$, star formation rates of ~ 21-295 $M_{\odot}$ yr$^{-1}$, and several have potential nearby companions. Our finding of bars at $z$ ~1.1-2.3 demonstrates the early onset of such instabilities and supports simulations where bars form early in massive dynamically cold disks. It also suggests that if these bars at lookback times of 8-10 Gyr survive out to present epochs, bar-driven secular processes may operate over a long time and have a significant impact on some galaxies by z ~ 0., Comment: 16 pages, 5 figures. Accepted for Publication in Astrophysical Journal Letters

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 comprehensive release of photometric redshifts (photo-z's) from the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS) team. We use statistics based upon the Quantile-Quantile (Q--Q) plot to identify biases and signatures of underestimated or overestimated errors in photo-z probability density functions (PDFs) produced by six groups in the collaboration; correcting for these effects makes the resulting PDFs better match the statistical definition of a PDF. After correcting each group's PDF, we explore three methods of combining the different groups' PDFs for a given object into a consensus curve. Two of these methods are based on identifying the minimum f-divergence curve, i.e., the PDF that is closest in aggregate to the other PDFs in a set (analogous to the median of an array of numbers). We demonstrate that these techniques yield improved results using sets of spectroscopic redshifts independent of those used to optimize PDF modifications. The best photo-z PDFs and point estimates are achieved with the minimum f-divergence using the best 4 PDFs for each object (mFDa4) and the Hierarchical Bayesian (HB4) methods, respectively. The HB4 photo-z point estimates produced $\sigma_{\rm NMAD} = 0.0227/0.0189$ and $|\Delta z/(1+z)| > 0.15$ outlier fraction = 0.067/0.019 for spectroscopic and 3D-HST redshifts, respectively. Finally, we describe the structure and provide guidance for the use of the CANDELS photo-z catalogs, which are available at https://archive.stsci.edu/hlsp/candels., Comment: 35 pages, 19 figures, published in ApJ, data available at https://archive.stsci.edu/hlsp/candels

We report on the host properties of five X-ray luminous Active Galactic Nuclei (AGN) identified at $3 < z < 5$ in the first epoch of imaging from the Cosmic Evolution Early Release Science Survey (CEERS). Each galaxy has been imaged with the \textit{James Webb Space Telescope} (\jwst) Near-Infrared Camera (NIRCam), which provides spatially resolved, rest-frame optical morphologies at these redshifts. We also derive stellar masses and star formation rates for each host galaxy by fitting its spectral energy distribution using a combination of galaxy and AGN templates. The AGN hosts have an average stellar mass of ${\rm log}(M_{*}/{\rm M_{\odot}} )= 11.0$, making them among the most massive galaxies detected at this redshift range in the current CEERS pointings, even after accounting for nuclear light from the AGN. We find that three of the AGN hosts have spheroidal morphologies, one is a bulge-dominated disk and one host is dominated by point-like emission. None are found to show strong morphological disturbances that might indicate a recent interaction or merger event. Notably, all four of the resolved hosts have rest-frame optical colors consistent with a quenched or post-starburst stellar population. The presence of AGN in passively evolving galaxies at $z>3$ is significant because a rapid feedback mechanism is required in most semi-analytic models and cosmological simulations to explain the growing population of massive quiescent galaxies observed at these redshifts. Our findings are in general agreement with this picture and show that AGN can continue to inject energy into these systems after their star formation is curtailed, possibly helping to maintain their quiescent state., Comment: Submitted to ApJ; 12 Pages, 6 Figures

Lyman Break Galaxy (LBG) candidates at z>10 are rapidly being identified in JWST/NIRCam observations. Due to the (redshifted) break produced by neutral hydrogen absorption of rest-frame UV photons, these sources are expected to drop out in the bluer filters while being well detected in redder filters. However, here we show that dust-enshrouded star-forming galaxies at lower redshifts (z<7) may also mimic the near-infrared (near-IR) colors of z>10 LBGs, representing potential contaminants in LBG candidate samples. First, we analyze CEERS-DSFG-1, a NIRCam dropout undetected in the F115W and F150W filters but detected at longer wavelengths. Combining the JWST data with (sub)millimeter constraints, including deep NOEMA interferometric observations, we show that this source is a dusty star-forming galaxy (DSFG) at z~5.1. We also present a tentative 2.6sigma SCUBA-2 detection at 850um around a recently identified z~16 LBG candidate in the same field and show that, if the emission is real and associated with this candidate, the available photometry is consistent with a z~5 dusty galaxy with strong nebular emission lines despite its blue near-IR colors. Further observations on this candidate are imperative to mitigate the low confidence of this tentative submillimeter emission and its positional uncertainty. Our analysis shows that robust (sub)millimeter detections of NIRCam dropout galaxies likely imply z=4-6 redshift solutions, where the observed near-IR break would be the result of a strong rest-frame optical Balmer break combined with high dust attenuation and strong nebular line emission, rather than the rest-frame UV Lyman break. This provides evidence that DSFGs may contaminate searches for ultra high-redshift LBG candidates from JWST observations., Comment: Published in The Astrophysical Journal Letters (updated to match the published version)

We report the discovery of a candidate galaxy with a photo-z of z~12 in the first epoch of the JWST Cosmic Evolution Early Release Science (CEERS) Survey. Following conservative selection criteria we identify a source with a robust z_phot = 11.8^+0.3_-0.2 (1-sigma uncertainty) with m_F200W=27.3, and >7-sigma detections in five filters. The source is not detected at lambda < 1.4um in deep imaging from both HST and JWST, and has faint ~3-sigma detections in JWST F150W and HST F160W, which signal a Ly-alpha break near the red edge of both filters, implying z~12. This object (Maisie's Galaxy) exhibits F115W-F200W > 1.9 mag (2-sigma lower limit) with a blue continuum slope, resulting in 99.6% of the photo-z PDF favoring z > 11. All data quality images show no artifacts at the candidate's position, and independent analyses consistently find a strong preference for z > 11. Its colors are inconsistent with Galactic stars, and it is resolved (r_h = 340 +/- 14 pc). Maisie's Galaxy has log M*/Msol ~ 8.5 and is highly star-forming (log sSFR ~ -8.2 yr^-1), with a blue rest-UV color (beta ~ -2.5) indicating little dust though not extremely low metallicity. While the presence of this source is in tension with most predictions, it agrees with empirical extrapolations assuming UV luminosity functions which smoothly decline with increasing redshift. Should followup spectroscopy validate this redshift, our Universe was already aglow with galaxies less than 400 Myr after the Big Bang., Comment: 17 pages, 6 figures, 3 tables, ApJL in press. Summary of changes from original submission: Improvements in astrometry generated a weak detection in F150W that reduces the photo-z to 11.8 but does not increase the likelihood of lower-z solutions. A full discussion of changes from the original version is available at: https://web.corral.tacc.utexas.edu/ceersdata/papers/Maisie_update.pdf

We present rest-frame optical emission-line flux ratio measurements for five $z>5$ galaxies observed by the JWST Near-Infared Spectrograph (NIRSpec) in the SMACS 0723 Early Release Observations. We add several quality-control and post-processing steps to the NIRSpec pipeline reduction products in order to ensure reliable relative flux calibration of emission lines that are closely separated in wavelength, despite the uncertain \textit{absolute} spectrophotometry of the current version of the reductions. Compared to $z\sim3$ galaxies in the literature, the $z>5$ galaxies have similar [OIII]$\lambda$5008/H$\beta$ ratios, similar [OIII]$\lambda$4364/H$\gamma$ ratios, and higher ($\sim$0.5 dex) [NeIII]$\lambda$3870/[OII]$\lambda$3728 ratios. We compare the observations to MAPPINGS V photoionization models and find that the measured [NeIII]$\lambda$3870/[OII]$\lambda$3728, [OIII]$\lambda$4364/H$\gamma$, and [OIII]$\lambda$5008/H$\beta$ emission-line ratios are consistent with an interstellar medium that has very high ionization ($\log(Q) \simeq 8-9$, units of cm~s$^{-1}$), low metallicity ($Z/Z_\odot \lesssim 0.2$), and very high pressure ($\log(P/k) \simeq 8-9$, units of cm$^{-3}$). The combination of [OIII]$\lambda$4364/H$\gamma$ and [OIII]$\lambda$(4960+5008)/H$\beta$ line ratios indicate very high electron temperatures of $4.1<\log(T_e/{\rm K})<4.4$, further implying metallicities of $Z/Z_\odot \lesssim 0.2$ with the application of low-redshift calibrations for ``$T_e$-based'' metallicities. These observations represent a tantalizing new view of the physical conditions of the interstellar medium in galaxies at cosmic dawn., Comment: Accepted for publication in AAS Journals. 14 pages, 6 figures, 3 tables

The 3D-Drift And SHift (3D-DASH) program is a \textit{Hubble Space Telescope} WFC3 F160W imaging and G141 grism survey of the equatorial COSMOS field. 3D-DASH extends the legacy of HST near-infrared imaging and spectroscopy to degree-scale swaths of the sky, enabling the identification and study of distant galaxies ($z>2$) that are rare or in short-lived phases of galaxy evolution at rest-frame optical wavelengths. Furthermore, when combined with existing ACS/F814W imaging, the program facilitates spatially-resolved studies of the stellar populations and dust content of intermediate-redshift ($0.5

We study the intrinsic 3D shapes of quiescent galaxies over the last half of cosmic history based on their axial ratio distribution. To this end, we construct a sample of unprecedented size, exploiting multi-wavelength $u$-to-$K_s$ photometry from the deep wide area surveys KiDS+VIKING paired with high-quality $i$-band imaging from HSC-SSP. Dependencies of the shapes on mass, redshift, photometric bulge prominence and environment are considered. For comparison, the intrinsic shapes of quenched galaxies in the IllustrisTNG simulations are analyzed and contrasted to their formation history. We find that over the full $0 10^{11}\ M_{\odot}$, with the effect being most pronounced at lower redshifts. In TNG, the most massive galaxies feature the highest ex-situ stellar mass fractions, pointing to violent relaxation via mergers as the mechanism responsible for their 3D shape transformation. Larger differences between observed and simulated shapes are found at low to intermediate masses. At any mass, the most spheroidal quiescent galaxies in TNG feature the highest bulge mass fractions, and conversely observed quiescent galaxies with the highest bulge-to-total ratios are found to be intrinsically the roundest. Finally, we detect an environmental influence on galaxy shape, at least at the highest masses, such that at fixed mass and redshift quiescent galaxies tend to be rounder in denser environments., Comment: Accepted by MNRAS. 18 pages, 11 figures, 2 tables

This paper investigates the neutral gas phase of galactic winds via the Na I D$\lambda\lambda 5890,5895${\AA} feature within $z \sim 0.04$ MaNGA galaxies, and directly compares their incidence and strength to the ionized winds detected within the same parent sample. We find evidence for neutral outflows in 127 galaxies ($\sim 5$ per cent of the analysed line-emitting sample). Na I D winds are preferentially seen in galaxies with dustier central regions and both wind phases are more often found in systems with elevated SFR surface densities, especially when there has been a recent upturn in the star formation activity according to the SFR$_{5Myr}$/SFR$_{800Myr}$ parameter. We find the ionized outflow kinematics to be in line with what we measure in the neutral phase. This demonstrates that, despite their small contributions to the total outflow mass budget, there is value to collecting empirical measurements of the ionized wind phase to provide information on the bulk motion in the outflow. Depending on dust corrections applied to the ionized gas diagnostics, the neutral phase has $\sim 1.2 - 1.8$ dex higher mass outflow rates ($\dot{M}_{out}$), on average, compared to the ionized phase. We quantify scaling relations between $\dot{M}_{out}$ and the strengths of the physical wind drivers (SFR, $L_{AGN}$). Using a radial-azimuthal stacking method, and by considering inclination dependencies, we find results consistent with biconical outflows orthogonal to the disk plane. Our work complements other multi-phase outflow studies in the literature which consider smaller samples, more extreme objects, or proceed via stacking of larger samples., Comment: This is a pre-copyedited, author-produced PDF of an article accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS) following peer review

While the slope of the dust attenuation curve ($\delta$) is found to correlate with effective dust attenuation ($A_V$) as obtained through spectral energy distribution (SED) fitting, it remains unknown how the fitting degeneracies shape this relation. We examine the degeneracy effects by fitting SEDs of a sample of local star-forming galaxies (SFGs) selected from the Galaxy And Mass Assembly survey, in conjunction with mock galaxy SEDs of known attenuation parameters. A well-designed declining starburst star formation history is adopted to generate model SED templates with intrinsic UV slope ($\beta_0$) spanning over a reasonably wide range. The best-fitting $\beta_0$ for our sample SFGs shows a wide coverage, dramatically differing from the limited range of $\beta_0<-2.2$ for a starburst of constant star formation. Our results show that strong degeneracies between $\beta_0$, $\delta$, and $A_V$ in the SED fitting induce systematic biases leading to a false $A_V$--$\delta$ correlation. Our simulation tests reveal that this relationship can be well reproduced even when a flat $A_V$--$\delta$ relation is taken to build the input model galaxy SEDs. The variations in best-fitting $\delta$ are dominated by the fitting errors. We show that assuming a starburst with constant star formation in SED fitting will result in a steeper attenuation curve, smaller degeneracy errors, and a stronger $A_V$--$\delta$ relation. Our findings confirm that the $A_V$--$\delta$ relation obtained through SED fitting is likely driven by the systematic biases induced by the fitting degeneracies between $\beta_0$, $\delta$, and $A_V$., Comment: 21 pages, 13 figures, accepted for publication in the MNRAS, Comments welcome!

Radial metallicity trends provide a key indicator of physical processes such as star formation and radial gas migration within a galaxy. Large IFU surveys allow for detailed studies of these radial variations, with recent observations detecting central dips in the metallicity, which may trace the impact of various evolutionary processes. However, the origin of these dips has not been conclusively determined, with suggestions that they may be diagnostic dependent. In this paper, we use the SDSS-IV MaNGA survey to investigate whether the observed dips represent genuine decreases in the central metallicity, or if they could be an artefact of the diagnostic used. Using a sub-sample of 758 local star-forming galaxies at low inclinations, we investigate in detail the impact of using different strong line diagnostics on the shapes of the returned profiles, and the prevalence of dips. We find no clear evidence of the dips being caused by changing values of the ionisation parameter within galaxies. To investigate physical causes, we explore both global and spatially-resolved parameters, finding that galaxies exhibiting central dips in the O3N2 metallicity profile have on average lower H$\alpha$EW values out to $R/R_\rm{e} \sim 1.5$, and higher values of D$_N$(4000) in the central regions. We additionally find a higher prevalence of dips in galaxies with high stellar mass, and lower values of global specific star formation rate, suggesting a possible link to central quenching. Nevertheless, these results are dependent on the diagnostic used, suggesting caution should be taken when interpreting observed features in galaxy metallicity gradients., Comment: Accepted for publication in MNRAS. 23 pages; 21 figures

In this work, we investigate the strength and impact of ionised gas outflows within $z \sim 0.04$ MaNGA galaxies. We find evidence for outflows in 322 galaxies ($12\%$ of the analysed line-emitting sample), 185 of which show evidence for AGN activity. Most outflows are centrally concentrated with a spatial extent that scales sublinearly with $R_{\rm e}$. The incidence of outflows is enhanced at higher masses, central surface densities and deeper gravitational potentials, as well as at higher SFR and AGN luminosity. We quantify strong correlations between mass outflow rates and the mechanical drivers of the outflow of the form $\dot{M}_{\rm out} \propto \rm SFR^{0.97}$ and $\dot{M}_{\rm out} \propto L_{\rm AGN}^{0.55}$. We derive a master scaling relation describing the mass outflow rate of ionised gas as a function of $M_{\star}$, SFR, $R_{\rm e}$ and $L_{\rm AGN}$. Most of the observed winds are anticipated to act as galactic fountains, with the fraction of galaxies with escaping winds increasing with decreasing potential well depth. We further investigate the physical properties of the outflowing gas finding evidence for enhanced attenuation in the outflow, possibly due to metal-enriched winds, and higher excitation compared to the gas in the galactic disk. Given that the majority of previous studies have focused on more extreme systems with higher SFRs and/or more luminous AGN, our study provides a unique view of the non-gravitational gaseous motions within `typical' galaxies in the low-redshift Universe, where low-luminosity AGN and star formation contribute jointly to the observed outflow phenomenology., Comment: Accepted for publication in MNRAS, 27 pages, Fig 7 & 8 for scaling wind strength with drivers, Fig 10 for master scaling

We compare the star forming main sequence (SFMS) -- both integrated and resolved on 1kpc scales -- between the high-resolution TNG50 simulation of IllustrisTNG and observations from the 3D-HST slitless spectroscopic survey at z~1. Contrasting integrated star formation rates (SFRs), we find that the slope and normalization of the star-forming main sequence in TNG50 are quantitatively consistent with values derived by fitting observations from 3D-HST with the Prospector Bayesian inference framework. The previous offsets of 0.2-1dex between observed and simulated main sequence normalizations are resolved when using the updated masses and SFRs from Prospector. The scatter is generically smaller in TNG50 than in 3D-HST for more massive galaxies with M_*>10^10Msun, even after accounting for observational uncertainties. When comparing resolved star formation, we also find good agreement between TNG50 and 3D-HST: average specific star formation rate (sSFR) radial profiles of galaxies at all masses and radii below, on, and above the SFMS are similar in both normalization and shape. Most noteworthy, massive galaxies with M_*>10^10.5Msun, which have fallen below the SFMS due to ongoing quenching, exhibit a clear central SFR suppression, in both TNG50 and 3D-HST. In TNG this inside-out quenching is due to the supermassive black hole (SMBH) feedback model operating at low accretion rates. In contrast, the original Illustris simulation, without this same physical SMBH mechanism, does not reproduce the central SFR profile suppression seen in data. The observed sSFR profiles provide support for the TNG quenching mechanism and how it affects gas on kiloparsec scales in the centers of galaxies., Comment: Main results in Figs 1 and 3. The TNG50 data is now fully publicly available at https://www.tng-project.org/

The Cepheid Leavitt Law (LL), also known as the Period-Luminosity relation, is a crucial tool for assembling the cosmic distance ladder. By combining data from the OGLE-IV catalogue with mid-infrared photometry from the Spitzer Space Telescope, we have determined the $3.6$ $\mu$m and $4.5$ $\mu$m LLs for the Magellanic Clouds using around 5000 fundamental-mode Classical Cepheids. Mean magnitudes were determined using a Monte Carlo Markov Chain (MCMC) template fitting procedure, with template light curves constructed from a subsample of these Cepheids with fully-phased, well-sampled light curves. The dependence of the Large Magellanic Cloud LL coefficients on various period cuts was tested, in addition to the linearity of the relationship. The zero point of the LL was calibrated using the parallaxes of Milky Way Cepheids from the Hubble Space Telescope and Gaia Data Release 2. Our final calibrated relations are $M_{[3.6]} = -3.246(\pm 0.008)(\log(P)-1.0)-5.784(\pm 0.030)$ and $M_{[4.5]} = -3.162(\pm 0.008)(\log(P)-1.0)-5.751(\pm 0.030)$., Comment: Accepted by MNRAS, the full Cepheid catalogue accompanying this article will be available soon

We present 0.2arcsec-resolution Atacama Large Millimeter/submillimeter Array observations at 870 $\mu$m in a stellar mass-selected sample of 85 massive ($M_\mathrm{star}>10^{11}~M_\odot$) star-forming galaxies (SFGs) at z=1.9-2.6 in the 3D-HST/CANDELS fields of UDS and GOODS-S. We measure the effective radius of the rest-frame far-infrared (FIR) emission for 62 massive SFGs. They are distributed over wide ranges of FIR size from $R_\mathrm{e,FIR}=$0.4 kpc to $R_\mathrm{e,FIR}=$6 kpc. The effective radius of the FIR emission is smaller by a factor of 2.3$^{+1.9}_{-1.0}$ than the effective radius of the optical emission and by a factor of 1.9$^{+1.9}_{-1.0}$ smaller than the half-mass radius. Even with taking into account potential extended components, the FIR size would change by ~10%. By combining the spatial distributions of the FIR and optical emission, we investigate how galaxies change the effective radius of the optical emission and the stellar mass within a radius of 1 kpc, $M_\mathrm{1kpc}$. The compact starburst puts most of massive SFGs on the mass--size relation for quiescent galaxies (QGs) at z~2 within 300 Myr if the current star formation activity and its spatial distribution are maintained. We also find that within 300 Myr, ~38% of massive SFGs can reach the central mass of $M_\mathrm{1kpc}=10^{10.5}~M_\odot$, which is around the boundary between massive SFGs and QGs. These results suggest an outside-in transformation scenario in which a dense core is formed at the center of a more extended disk, likely via dissipative in-disk inflows. Synchronized observations at ALMA 870 $\mu$m and JWST 3-4 $\mu$m will explicitly verify this scenario., Comment: 25 pages, 15 figures, 3 tables, accepted for publication in ApJ

We study the dynamical properties of massive quiescent galaxies at $1.4 < z < 2.1$ using deep Hubble Space Telescope WFC3/F160W imaging and a combination of literature stellar velocity dispersion measurements and new near-infrared spectra obtained using KMOS on the ESO VLT. We use these data to show that the typical dynamical-to-stellar mass ratio has increased by $\sim$0.2 dex from $z = 2$ to the present day, and investigate this evolution in the context of possible changes in the stellar initial mass function (IMF) and/or fraction of dark matter contained within the galaxy effective radius, $f_\mathrm{DM}$. Comparing our high-redshift sample to their likely descendants at low-redshift, we find that $f_\mathrm{DM}$ has increased by a factor of more than 4 since $z \approx 1.8$, from $f_\mathrm{DM}$ = $6.6\pm1.0$% to $\sim$24%. The observed increase appears robust to changes in the methods used to estimate dynamical masses or match progenitors and descendants. We quantify possible variation of the stellar IMF through the offset parameter $\alpha$, defined as the ratio of dynamical mass in stars to the stellar mass estimated using a Chabrier IMF. We demonstrate that the correlation between stellar velocity dispersion and $\alpha$ reported among quiescent galaxies at low-redshift is already in place at $z = 2$, and argue that subsequent evolution through (mostly minor) merging should act to preserve this relation while contributing significantly to galaxies overall growth in size and stellar mass., Comment: 40 pages, 21 figures; accepted for publication in ApJ

We compare the infrared excess (IRX) and Balmer decrement (${\rm H\alpha/H\beta }$) as dust attenuation indicators in relation to other galaxy parameters using a sample of $\sim$32 000 local star-forming galaxies (SFGs) carefully selected from SDSS, GALEX and WISE. While at fixed ${\rm H\alpha/H\beta }$, IRX turns out to be independent on galaxy stellar mass, the Balmer decrement does show a strong mass dependence at fixed IRX. We find the discrepancy, parameterized by the color excess ratio $R_{\rm EBV} \equiv E(B-V)_{\rm IRX}/E(B-V)_{\rm H\alpha/H\beta }$, is not dependent on the gas-phase metallicity and axial ratio but on the specific star formation rate (SSFR) and galaxy size ($R_{\rm e}$) following $R_{\rm EBV}=0.79+0.15\log({\rm SSFR}/R_{\rm e}^{2})$. This finding reveals that the nebular attenuation as probed by the Balmer decrement becomes increasingly larger than the global (stellar) attenuation of SFGs with decreasing SSFR surface density. This can be understood in the context of an enhanced fraction of intermediate-age stellar populations that are less attenuated by dust than the HII region-traced young population, in conjunction with a decreasing dust opacity of the diffuse ISM when spreading over a larger spatial extent. Once the SSFR surface density of an SFG is known, the conversion between attenuation of nebular and stellar emission can be well estimated using our scaling relation., Comment: 9 pages, 5 figures, Accepted for publication in ApJ

We present a WFC3 F160W ($H$-band) selected catalog in the CANDELS/GOODS-N field containing photometry from the ultraviolet (UV) to the far-infrared (IR), photometric redshifts and stellar parameters derived from the analysis of the multi-wavelength data. The catalog contains 35,445 sources over the 171 arcmin$^{2}$ of the CANDELS F160W mosaic. The 5$\sigma$ detection limits (within an aperture of radius 0\farcs17) of the mosaic range between $H=27.8$, 28.2 and 28.7 in the wide, intermediate and deep regions, that span approximately 50\%, 15\% and 35\% of the total area. The multi-wavelength photometry includes broad-band data from UV (U band from KPNO and LBC), optical (HST/ACS F435W, F606W, F775W, F814W, and F850LP), near-to-mid IR (HST/WFC3 F105W, F125W, F140W and F160W, Subaru/MOIRCS Ks, CFHT/Megacam K, and \spitzer/IRAC 3.6, 4.5, 5.8, 8.0 $\mu$m) and far IR (\spitzer/MIPS 24$\mu$m, HERSCHEL/PACS 100 and 160$\mu$m, SPIRE 250, 350 and 500$\mu$m) observations. In addition, the catalog also includes optical medium-band data (R$\sim50$) in 25 consecutive bands, $\lambda=500$ to 950~nm, from the SHARDS survey and WFC3 IR spectroscopic observations with the G102 and G141 grisms (R$\sim210$ and 130). The use of higher spectral resolution data to estimate photometric redshifts provides very high, and nearly uniform, precision from $z=0-2.5$. The comparison to 1,485 good quality spectroscopic redshifts up to $z\sim3$ yields $\Delta z$/(1+$z_{\rm spec}$)$=$0.0032 and an outlier fraction of $\eta=$4.3\%. In addition to the multi-band photometry, we release added-value catalogs with emission line fluxes, stellar masses, dust attenuations, UV- and IR- based star formation rates and rest-frame colors., Comment: The catalogs, images and spectra are available in Rainbow-slicer (https://bit.ly/2OxKKx1 ), navigator (https://bit.ly/2GDS180 ) and MAST (https://bit.ly/2YtoBQ4 ). In addition to the photometry and other added-value catalogs we release UV+IR star formation rates based on Spitzer (MIPS) and Herschel (PACS and SPIRE) in the 5 CANDELS fields

We present the $0.6

We investigate dust obscuration as parameterised by the infrared excess IRX$\equiv$$L_{\rm IR}/L_{\rm UV}$ in relation to global galaxy properties, using a sample of $\sim$32$\,$000 local star-forming galaxies (SFGs) selected from SDSS, GALEX and WISE. We show that IRX generally correlates with stellar mass ($M_\ast$), star formation rate (SFR), gas-phase metallicity ($Z$), infrared luminosity ($L_{\rm IR}$) and the half-light radius ($R_{\rm e}$). A weak correlation of IRX with axial ratio (b/a) is driven by the inclination and thus seen as a projection effect. By examining the tightness and the scatter of these correlations, we find that SFGs obey an empirical relation of the form $IRX$=$10^\alpha\,(L_{\rm IR})^{\beta}\,R_{\rm e}^{-\gamma}\,(b/a)^{-\delta}$ where the power-law indices all increase with metallicity. The best-fitting relation yields a scatter of $\sim$0.17$\,$dex and no dependence on stellar mass. Moreover, this empirical relation also holds for distant SFGs out to $z=3$ in a population-averaged sense, suggesting it to be universal over cosmic time. Our findings reveal that IRX approximately increases with $L_{\rm IR}/R_{\rm e}^{[1.3 - 1.5]}$ instead of $L_{\rm IR}/R_{\rm e}^{2}$ (i.e., surface density). We speculate this may be due to differences in the spatial extent of stars versus star formation and/or complex star-dust geometries. We conclude that not stellar mass but IR luminosity, metallicity and galaxy size are the key parameters jointly determining dust obscuration in SFGs., Comment: 22 pages, 14 figures, accepted for publication in the MNRAS

We investigate the resolved relation between local extinction and star formation surface density within nearby star-forming galaxies selected from the MaNGA survey. Balmer decrement measurements imply an extinction of the H{\alpha} line emission which scales approximately linearly with the logarithm of the star formation surface density: $ A_{H{\alpha}} = 0.46 \log(\Sigma_{SFR}) + 1.53$. Secondary dependencies are observed such that, at a given $\Sigma_{SFR}$, regions of lower metallicity and/or enhanced H{\alpha} equivalent width (EW) suffer less obscuration than regions of higher metallicity and/or lower H{\alpha} EW. Spaxels lying above the mean relation also tend to belong to galaxies that are more massive, larger and viewed under higher inclination than average. We present a simple model in which the observed trends can be accounted for by a metallicity-dependent scaling between $\Sigma_{SFR}$ and $\Sigma_{dust}$ via a super-linear Kennicutt-Schmidt relation ($n_{KS} \sim 1.47$) and a dust-to-gas ratio which scales linearly with metallicity (DGR($Z_{\odot}$) = 0.01). The relation between the resulting total dust column and observed effective extinction towards nebular regions requires a geometry for the relative distribution of H{\alpha} emitting regions and dust that deviates from a uniform foreground screen and also from an entirely homogeneous mixture of dust and emitting sources. The best-fit model features an H{\alpha} EW and galactocentric distance dependent fraction of the dust mass in a clumpy foreground screen in front of a homogeneous mixture., Comment: Accepted for publication in ApJ. Figures 2 and 3 show the observed star formation - extinction relation. Figures 9 and 10 show our favored model

We investigate the relationship between star formation activity and outflow properties on kiloparsec scales in a sample of 28 star forming galaxies at $z\sim$ 2-2.6, using adaptive optics assisted integral field observations from SINFONI on the VLT. The narrow and broad components of the H$\alpha$ emission are used to simultaneously determine the local star formation rate surface density ($\Sigma_{\rm SFR}$), and the outflow velocity $v_{\rm out}$ and mass outflow rate $\dot{M}_{\rm out}$, respectively. We find clear evidence for faster outflows with larger mass loading factors at higher $\Sigma_{\rm SFR}$. The outflow velocities scale as $v_{\rm out}$ $\propto$ $\Sigma_{\rm SFR}^{0.34 \pm 0.10}$, which suggests that the outflows may be driven by a combination of mechanical energy released by supernova explosions and stellar winds, as well as radiation pressure acting on dust grains. The majority of the outflowing material does not have sufficient velocity to escape from the galaxy halos, but will likely be re-accreted and contribute to the chemical enrichment of the galaxies. In the highest $\Sigma_{\rm SFR}$ regions the outflow component contains an average of $\sim$45% of the H$\alpha$ flux, while in the lower $\Sigma_{\rm SFR}$ regions only $\sim$10% of the H$\alpha$ flux is associated with outflows. The mass loading factor, $\eta$ = $\dot{M}_{\rm out}$/SFR, is positively correlated with $\Sigma_{\rm SFR}$ but is relatively low even at the highest $\Sigma_{\rm SFR}$: $\eta \lesssim$ 0.5 $\times$ (380 cm$^{-3}$/n$_e$). This may be in tension with the $\eta$ $\gtrsim$ 1 required by cosmological simulations, unless a significant fraction of the outflowing mass is in other gas phases and has sufficient velocity to escape the galaxy halos., Comment: 28 pages, 12 figures, 2 tables. Accepted for publication in ApJ

We present a cross-calibration of CO- and dust-based molecular gas masses at $z \leqslant 0.2$. Our results are based on a survey with the IRAM 30-m telescope collecting CO(1-0) measurements of 78 massive ($\log M_{\star} / M_{\odot} > 10$) galaxies with known gas-phase metallicities, and with IR photometric coverage from WISE (22 $\mu$m ) and Herschel SPIRE (250, 350, 500 $\mu$m). We find a tight relation ($\sim 0.17$ dex scatter) between the gas masses inferred from CO and dust continuum emission, with a minor systematic offset of 0.05 dex. The two methods can be brought into agreement by applying a metallicity-dependent adjustment factor ($\sim 0.13$ dex scatter). We illustrate that the observed offset is consistent with a scenario in which dust traces not only molecular gas, but also part of the ${\rm H \small I}$ reservoir, residing in the ${\rm H_2}$-dominated region of the galaxy. Observations of the CO(2-1) to CO(1-0) line ratio for two thirds of the sample indicate a narrow range in excitation properties, with a median ratio of luminosities $ \left\langle R_{21} \right\rangle \sim 0.64 $. Finally, we find dynamical mass constraints from spectral line profile fitting to agree well with the anticipated mass budget enclosed within an effective radius, once all mass components (stars, gas and dark matter) are accounted for., Comment: Accepted for publication in MNRAS, 31 pages, 17 figures, 2 tables

We present deep observations of a $z=1.4$ massive, star-forming galaxy in molecular and ionized gas at comparable spatial resolution (CO 3-2, NOEMA; H$\alpha$, LBT). The kinematic tracers agree well, indicating that both gas phases are subject to the same gravitational potential and physical processes affecting the gas dynamics. We combine the one-dimensional velocity and velocity dispersion profiles in CO and H$\alpha$ to forward-model the galaxy in a Bayesian framework, combining a thick exponential disk, a bulge, and a dark matter halo. We determine the dynamical support due to baryons and dark matter, and find a dark matter fraction within one effective radius of $f_{\rm DM}(\leq$$R_{e})=0.18^{+0.06}_{-0.04}$. Our result strengthens the evidence for strong baryon-dominance on galactic scales of massive $z\sim1-3$ star-forming galaxies recently found based on ionized gas kinematics alone., Comment: accepted for publication in the Astrophysical Journal Letters

A randomly chosen star in today's Universe is most likely to live in a galaxy with a stellar mass between that of the Milky Way and Andromeda. Yet it remains uncertain how the structural evolution of these bulge-disk systems proceeded. Most of the unobscured star formation we observe building Andromdeda progenitors at 0.790% of their star formation is reprocessed by dust and remains unaccounted for. Here we map 500micron dust continuum emission in an Andromeda progenitor at z=1.25 to probe where it is growing through dust-obscured star formation. Combining resolved dust measurements from the NOEMA interferometer with Hubble Space Telescope Halpha maps and multicolor imaging (including new UV data from the HDUV survey), we find a bulge growing by dust-obscured star formation: while the unobscured star formation is centrally suppressed, the dust continuum is centrally concentrated, filling in the ring-like structures evident in the Halpha and UV emission. Reflecting this, the dust emission is more compact than the optical/UV tracers of star formation with r_e(dust)=3.4kpc, r_e(Halpha)/r_e(dust)=1.4, and r_e(UV)/r_e(dust)=1.8. Crucially, however, the bulge and disk of this galaxy are building simultaneously; although the dust emission is more compact than the rest-optical emission (r_e(optical)/r_e(dust)=1.4), it is somewhat less compact than the stellar mass (r_e(M_*)/r_e(dust)=0.9). Taking the 500micron emission as a tracer of star formation, the expected structural evolution of this galaxy can be accounted for by star formation: it will grow in size by Delta(r_e)/Delta(M_*)~0.3 and central surface density by Delta(Sigma_cen)/Delta(M_*)~0.9. Finally, our observations are consistent with a picture in which merging and disk instabilities drive gas to the center of galaxies, boosting global star formation rates above the main sequence and building bulges., Comment: Submitted to ApJ, key new result of paper shown in Fig. 2

We performed a detailed study of the evolution of the star formation rate (SFR) and stellar mass of the brightest group galaxies (BGGs) and their relative contribution to the total baryon budget within $R_{200}$ ($f^{BGG}_{b,200}$). The sample comprises 407 BGGs selected from X-ray galaxy groups ($M_{200}=10^{12.8}-10^{14} \;M_{\odot}$) out to $z\sim1.3$ identified in the COSMOS, XMM-LSS, and AEGIS fields. We find that BGGs constitute two distinct populations of quiescent and star-forming galaxies and their mean SFR is $\sim2$ dex higher than the median SFR at $ z<1.3 $. Both the mean and the median SFRs decline with time by $>2$ dex. The mean (median) of stellar mass of BGGs has grown by $0.3$ dex since $z=1.3$ to the present day. We show that up to $\sim45\% $ of the stellar mass growth in a star-forming BGG can be due to its star-formation activity. With respect to $f^{BGG}_{b,200}$, we find it to increase with decreasing redshift by $\sim0.35$ dex while decreasing with halo mass in a redshift dependent manner. We show that the slope of the relation between $f^{BGG}_{b,200}$ and halo mass increases negatively with decreasing redshift. This trend is driven by an insufficient star-formation in BGGs, compared to the halo growth rate. We separately show the BGGs with the 20\% highest $f^{BGG}_{b,200}$ are generally non-star-forming galaxies and grow in mass by processes not related to star formation (e.g., dry mergers and tidal striping). We present the $ M_\star-M_h $ and $ M_\star/M_h-M_h $ relations and compare them with semi-analytic model predictions and a number of results from the literature. We quantify the intrinsic scatter in stellar mass of BGGs at fixed halo mass ($\sigma_{log M_{\star}}$) and find that $\sigma_{log M_{\star}}$ increases from 0.3 dex at $ z\sim0.2 $ to 0.5 dex at $ z\sim1.0 $ due to the bimodal distribution of stellar mass., Comment: Accepted for publication in MNRAS, 25 pages, 11 figures

The rate of major galaxy-galaxy merging is theoretically predicted to steadily increase with redshift during the peak epoch of massive galaxy development ($1{\leq}z{\leq}3$). We use close-pair statistics to objectively study the incidence of massive galaxies (stellar $M_{1}{\geq}2{\times}10^{10}M_{\odot}$) hosting major companions ($1{\leq}M_{1}/M_{2}{\leq}4$; i.e., $<$4:1) at six epochs spanning $0{<}z{<}3$. We select companions from a nearly complete, mass-limited ($\geq5{\times}10^{9}M_{\odot}$) sample of 23,696 galaxies in the five CANDELS fields and the SDSS. Using $5-50$ kpc projected separation and close redshift proximity criteria, we find that the major companion fraction $f_{\mathrm{mc}}(z)$ based on stellar mass-ratio (MR) selection increases from 6% ($z{\sim}0$) to 16% ($z{\sim}0.8$), then turns over at $z{\sim}1$ and decreases to 7% ($z{\sim}3$). Instead, if we use a major F160W flux ratio (FR) selection, we find that $f_{\mathrm{mc}}(z)$ increases steadily until $z=3$ owing to increasing contamination from minor (MR$>$4:1) companions at $z>1$. We show that these evolutionary trends are statistically robust to changes in companion proximity. We find disagreements between published results are resolved when selection criteria are closely matched. If we compute merger rates using constant fraction-to-rate conversion factors ($C_{\mathrm{merg,pair}}{=}0.6$ and $T_{\mathrm{obs,pair}}{=}0.65\mathrm{Gyr}$), we find that MR rates disagree with theoretical predictions at $z{>}1.5$. Instead, if we use an evolving $T_{\mathrm{obs,pair}}(z){\propto}(1+z)^{-2}$ from Snyder et al., our MR-based rates agree with theory at $0{<}z{<}3$. Our analysis underscores the need for detailed calibration of $C_{\mathrm{merg,pair}}$ and $T_{\mathrm{obs,pair}}$ as a function of redshift, mass and companion selection criteria to better constrain the empirical major merger history., Comment: 27 pages, 15 figures, 5 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal