Search

We present the star formation rate–stellar mass (SFR– M _* ) relation for galaxies in the Cosmic Evolution Early Release Science survey at 4.5 ≤ z ≤ 12. We model the JWST and Hubble Space Telescope rest-UV and rest-optical photometry of galaxies with flexible star formation histories (SFHs) using BAGPIPES. We consider SFRs averaged from the SFHs over 10 Myr (SFR _10 ) and 100 Myr (SFR _100 ), where the photometry probes SFRs on these timescales, effectively tracing nebular emission lines in the rest-optical (on ~10 Myr timescales) and the UV/optical continuum (on ~100 Myr timescales). We measure the slope, normalization and intrinsic scatter of the SFR– M _* relation, taking into account the uncertainty and the covariance of galaxy SFRs and M _* . From z ~ 5 to 9 there is larger scatter in the SFR _10 – M _* relation, with $\sigma (\mathrm{log}{\rm{S}}{\rm{F}}{{\rm{R}}}_{100})=0.4$ dex, compared to the SFR _100 – M _* relation, with $\sigma (\mathrm{log}SF{R}_{10})=0.1$ dex. This scatter increases with redshift and increasing stellar mass, at least out to z ~ 7. These results can be explained if galaxies at higher redshift experience an increase in star formation variability and form primarily in short, active periods, followed by a lull in star formation (i.e., “napping” phases). We see a significant trend in the ratio R _SFR = SFR _10 /SFR _100 in which, on average, R _SFR decreases with increasing stellar mass and increasing redshift. This yields a star formation “duty cycle” of ~40% for galaxies with $\mathrm{log}{M}_{* }/{M}_{\odot }\geqslant 9.3$ at z ~ 5, declining to ~20% at z ~ 9. Galaxies also experience longer lulls in star formation at higher redshift and at higher stellar mass, such that galaxies transition from periods of higher SFR variability at z ≳ 6 to smoother SFR evolution at z ≲ 4.5.

We present measurements of stellar populations properties of a z = 9.1 gravitationally lensed galaxy MACS1149-JD1 using deep James Webb Space Telescope NIRISS slitless spectroscopy as well as NIRISS and NIRCam imaging from the CAnadian NIRISS Unbiased Cluster Survey (CANUCS). The galaxy is split into four components. Three magnified ( μ ∼ 11) star-forming components are unresolved, giving intrinsic sizes of

Photometrically derived stellar masses are known to suffer from systematic uncertainties, particularly due to nebular emission contributions to the spectral energy distribution (SED). Using James Webb Space Telescope NIRCam imaging from the CAnadian NIRISS Unbiased Cluster Survey, we introduce a comparison study of photometrically derived redshifts and stellar masses based on two photometric catalogs of the same field spanning ∼0.4–4.5 μ m: one consisting solely of wide band photometry, and another employing a combination of wide and medium band photometry. We find that ∼70% of likely line emitters have consistent photometric redshifts between both catalogs, with a median stellar mass difference between the two catalogs of 5. These systematic differences caused by the poor spectral resolution of wide bands have implications for both ongoing and future planned observing programs that determine stellar mass and other physical properties of high-redshift galaxies solely via wide band photometry.

We present a sample of 88 candidate z ∼ 8.5–14.5 galaxies selected from the completed NIRCam imaging from the Cosmic Evolution Early Release Science survey. These data cover ∼90 arcmin ^2 (10 NIRCam pointings) in six broadband imaging filters and one medium-band imaging filter. With this sample we confirm at higher confidence early JWST conclusions that bright galaxies in this epoch are more abundant than predicted by most theoretical models. We construct the rest-frame ultraviolet luminosity functions at z ∼ 9, 11, and 14 and show that the space density of bright ( M _UV = −20) galaxies changes only modestly from z ∼ 14 to z ∼ 9, compared to a steeper increase from z ∼ 8 to z ∼ 4. While our candidates are photometrically selected, spectroscopic follow-up has now confirmed 13 of them, with only one significant interloper, implying that the fidelity of this sample is high. Successfully explaining the evidence for a flatter evolution in the number densities of UV-bright z > 10 galaxies may thus require changes to the dominant physical processes regulating star formation. While our results indicate that significant variations of dust attenuation with redshift are unlikely to be the dominant factor at these high redshifts, they are consistent with predictions from models that naturally have enhanced star formation efficiency and/or stochasticity. An evolving stellar initial mass function could also bring model predictions into better agreement with our results. Deep spectroscopic follow-up of a large sample of early galaxies can distinguish between these competing scenarios.

A crucial yet challenging task in galaxy evolution studies is the identification of distant merging galaxies, a task that suffers from a variety of issues ranging from telescope sensitivities and limitations to the inherently chaotic morphologies of young galaxies. In this paper, we use random forests and convolutional neural networks to identify high-redshift JWST Cosmic Evolution Early Release Science Survey (CEERS) galaxy mergers. We train these algorithms on simulated 3 < z < 5 CEERS galaxies created from the IllustrisTNG subhalo morphologies and the Santa Cruz SAM light cone. We apply our models to observed CEERS galaxies at 3 < z < 5. We find that our models correctly classify ∼60%–70% of simulated merging and nonmerging galaxies; better performance on the merger class comes at the expense of misclassifying more nonmergers. We could achieve more accurate classifications, as well as test for a dependency on physical parameters such as gas fraction, mass ratio, and relative orbits, by curating larger training sets. When applied to real CEERS galaxies using visual classifications as ground truth, the random forests correctly classified 40%–60% of mergers and nonmergers at 3 < z < 4 but tended to classify most objects as nonmergers at 4 < z < 5 (misclassifying ∼70% of visually classified mergers). On the other hand, the CNNs tended to classify most objects as mergers across all redshifts (misclassifying 80%–90% of visually classified nonmergers). We investigate what features the models find most useful, as well as the characteristics of false positives and false negatives, and also calculate merger rates derived from the identifications made by the models.

The exponential growth of astronomical literature poses significant challenges for researchers navigating and synthesizing general insights or even domain-specific knowledge. We present pathfinder , a machine learning framework designed to enable literature review and knowledge discovery in astronomy, focusing on semantic searching with natural language instead of syntactic searches with keywords. Utilizing state-of-the-art large language models (LLMs) and a corpus of 385,166 peer-reviewed papers from the Astrophysics Data System, pathfinder offers an innovative approach to scientific inquiry and literature exploration. Our framework couples advanced retrieval techniques with LLM-based synthesis to search astronomical literature by semantic context as a complement to currently existing methods that use keywords or citation graphs. It addresses complexities of jargon, named entities, and temporal aspects through time-based and citation-based weighting schemes. We demonstrate the tool’s versatility through case studies, showcasing its application in various research scenarios. The system’s performance is evaluated using custom benchmarks, including single-paper and multipaper tasks. Beyond literature review, pathfinder offers unique capabilities for reformatting answers in ways that are accessible to various audiences (e.g., in a different language or as simplified text), visualizing research landscapes, and tracking the impact of observatories and methodologies. This tool represents a significant advancement in applying artificial intelligence to astronomical research, aiding researchers at all career stages in navigating modern astronomy literature.

Galaxy formation and evolution involve a variety of effectively stochastic processes that operate over different timescales. The extended regulator model provides an analytic framework for the resulting variability (or “burstiness”) in galaxy-wide star formation due to these processes. It does this by relating the variability in Fourier space to the effective timescales of stochastic gas inflow, equilibrium, and dynamical processes influencing giant molecular clouds' creation and destruction using the power spectral density (PSD) formalism. We use the connection between the PSD and autocovariance function for general stochastic processes to reformulate this model as an autocovariance function, which we use to model variability in galaxy star formation histories (SFHs) using physically motivated Gaussian processes in log star formation rate (SFR) space. Using stellar population synthesis models, we then explore how changes in model stochasticity can affect spectral signatures across galaxy populations with properties similar to the Milky Way and present-day dwarfs, as well as at higher redshifts. We find that, even at fixed scatter, perturbations to the stochasticity model (changing timescales vs. overall variability) leave unique spectral signatures across both idealized and more realistic galaxy populations. Distributions of spectral features including H α and UV-based SFR indicators, H δ and Ca H and K absorption-line strengths, D _n (4000), and broadband colors provide testable predictions for galaxy populations from present and upcoming surveys with the Hubble Space Telescope, James Webb Space Telescope, and Nancy Grace Roman Space Telescope. The Gaussian process SFH framework provides a fast, flexible implementation of physical covariance models for the next generation of spectral energy distribution modeling tools. Code to reproduce our results can be found at https://github.com/kartheikiyer/GP-SFH .

We demonstrate that the inference of galaxy stellar masses via spectral energy distribution (SED) fitting techniques for galaxies formed in the first billion years after the Big Bang carries fundamental uncertainties owing to the loss of star formation history (SFH) information from the very first episodes of star formation in the integrated spectra of galaxies. While this early star formation can contribute substantially to the total stellar mass of high-redshift systems, ongoing star formation at the time of detection outshines the residual light from earlier bursts, hampering the determination of accurate stellar masses. As a result, order-of-magnitude uncertainties in stellar masses can be expected. We demonstrate this potential problem via direct numerical simulation of galaxy formation in a cosmological context. In detail, we carry out two cosmological simulations with significantly different stellar feedback models, which span a significant range in SFH burstiness. We compute the mock SEDs for these model galaxies at z = 7 via calculations of 3D dust radiative transfer, and then backward fit these SEDs with prospector SED fitting software. The uncertainties in derived stellar masses that we find for z > 7 galaxies motivate the development of new techniques and/or priors for SFH to model star formation in the early Universe.

We combine deep imaging data from the CEERS early release JWST survey and Hubble Space Telescope imaging from CANDELS to examine the size–mass relation of star-forming galaxies and the morphology–quenching relation at stellar masses M _⋆ ≥ 10 ^9.5 M _⊙ over the redshift range 0.5 < z < 5.5. In this study with a sample of 2450 galaxies, we separate star-forming and quiescent galaxies based on their star formation activity and confirm that star-forming and quiescent galaxies have different morphologies out to z = 5.5, extending the results of earlier studies out to higher redshifts. We find that star-forming and quiescent galaxies have typical Sérsic indices of n ∼ 1.3 and n ∼ 4.3, respectively. Focusing on star-forming galaxies, we find that the slope of the size–mass relation is nearly constant with redshift, as was found previously, but shows a modest increase at z ∼ 4.2. The intercept in the size–mass relation declines out to z = 5.5 at rates that are similar to what earlier studies found. The intrinsic scatter in the size–mass relation is relatively constant out to z = 5.5.

The 3D geometries of high-redshift galaxies remain 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 James Webb Space Telescope Cosmic Evolution Early Release Science observations with $\mathrm{log}{M}_{* }/{M}_{\odot }=9.0\mbox{--}10.5$ at z = 0.5–8.0. We reproduce previous results from the Hubble Space Telescope Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey in a fraction of the computing time and constrain the mean ellipticity, triaxiality, size, and covariances with samples as small as ∼50 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 ∼1 for $\mathrm{log}{M}_{* }/{M}_{\odot }=9.0\mbox{--}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\mbox{--}\mathrm{log}a$ diagram with an excess of low- b / a , large- $\mathrm{log}a$ galaxies. The dwarf prolate fraction rises from ∼25% at z = 0.5–1.0 to ∼50%–80% at z = 3–8. Our results imply a second kind of disk settling from oval (triaxial) to more circular (axisymmetric) shapes with time. We simultaneously constrain the 3D size–mass relation and its dependence on 3D geometry. High-probability prolate and oblate candidates show remarkably similar Sérsic indices ( n ∼ 1), nonparametric 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.

We present a new sample of 158 galaxies at redshift z > 7.5 selected from deep James Webb Space Telescope (JWST) NIRCam imaging of five widely separated sight lines in the CANUCS survey. Two-thirds of the pointings and 80% of the galaxies are covered by 12–14 NIRCam filters, including seven to nine medium bands, providing accurate photometric redshifts and robustness against low-redshift interlopers. A sample of 28 galaxies at z > 7.5 with spectroscopic redshifts shows a low systematic offset and scatter in the difference between photometric and spectroscopic redshifts. We derive the galaxy UV luminosity function at redshifts 8–12, finding a slightly higher normalization than previously seen with the Hubble Space Telescope at redshifts 8–10. We observe a steeper decline in the galaxy space density from z = 8 to 12 than found by most JWST Cycle 1 studies. In particular, we find only eight galaxies at z > 10 and none at z > 12.5, with no z > 10 galaxies brighter than F277W AB = 28 or M _UV = −20 in our unmasked, delensed survey area of 53.4 arcmin ^2 . We attribute the lack of bright z > 10 galaxies in CANUCS compared to GLASS and CEERS to intrinsic variance in the galaxy density along different sight lines. The evolution in the CANUCS luminosity function between z = 8 and 12 is comparable to that predicted by simulations that assume a standard star formation efficiency without invoking any special adjustments.

The physical processes responsible for shaping how galaxies form and quench over time leave imprints on both the spatial (galaxy morphology) and temporal (star formation history; SFH) tracers that we use to study galaxies. While the morphology–SFR connection is well studied, the correlation with past star formation activity is not as well understood. To quantify this, we present Katachi (形), an interpretable convolutional neural network framework that learns the connection between the factors regulating star formation in galaxies on different spatial and temporal scales. Katachi is trained on 9904 galaxies at 0.02 < z < 0.1 in the SDSS-IV MaNGA DR17 sample to predict stellar mass ( M _* ; root mean square error (RSME) 0.22 dex), current star formation rate (SFR; RMSE 0.31 dex), and half-mass time ( t _50 ; RMSE 0.23 dex). This information allows us to reconstruct nonparametric SFHs for each galaxy from gri imaging alone. To quantify the morphological features informing the SFH predictions, we use SHAP (SHapley Additive exPlanations). We recover the expected trends of M _* governed by the growth of galaxy bulges, as well as those of SFR correlating with spiral arms and other star-forming regions. We also find that the SHAP maps of D4000 are more complex than those of M _* and SFR, and that morphology is correlated with t _50 even at fixed mass and SFR. Katachi serves as a scalable public framework to predict galaxy properties from large imaging surveys including Rubin, Roman, and Euclid, with large data sets of high signal-to-noise ratio imaging across limited photometric bands.

We analyze the evolution of massive (log _10 [ M _⋆ / M _⊙ ] > 10) galaxies at z ∼ 1–4 selected from JWST Cosmic Evolution Early Release Survey (CEERS). We infer the physical properties of all galaxies in the CEERS NIRCam imaging through spectral energy distribution (SED) fitting with dense basis to select a sample of high-redshift massive galaxies. Where available we include constraints from additional CEERS observing modes, including 18 sources with MIRI photometric coverage, and 28 sources with spectroscopic confirmations from NIRSpec or NIRCam WFSS. We sample the recovered posteriors in stellar mass from SED fitting to infer the volume densities of massive galaxies across cosmic time, taking into consideration the potential for sample contamination by active galactic nuclei. We find that the evolving abundance of massive galaxies tracks expectations based on a constant baryon conversion efficiency in dark matter halos for z ∼ 1–4. At higher redshifts, we observe an excess abundance of massive galaxies relative to this simple model, resulting in a shallower decline of observed volume densities of massive galaxies. These higher abundances can be explained by modest changes to star formation physics and/or the efficiencies with which star formation occurs in massive dark matter halos, and are not in tension with modern cosmology.

The UltraViolet imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey Fields (UVCANDELS) program provides deep Hubble Space Telescope (HST) F275W and F435W imaging over four CANDELS fields (GOODS-N, GOODS-S, COSMOS, and Extended Groth Strip). We combine this newly acquired UV imaging with existing HST imaging from CANDELS as well as existing ancillary data to obtain robust photometric redshifts and reliable estimates for galaxy physical properties for over 150,000 galaxies in the ∼430 arcmin ^2 UVCANDELS area. Here, we leverage the power of the new UV photometry to not only improve the photometric redshift measurements in these fields, but also constrain the full redshift probability distribution combining multiple redshift-fitting tools. Furthermore, using the full UV-to-IR photometric data set, we measure the galaxy physical properties by fitting templates from population synthesis models with two different parameterizations (flexible and fixed form) of the star formation histories (SFHs). Compared to the flexible SFH parameterization, we find that the fixed-form SFHs systematically underestimate the galaxy stellar masses, both at the low-mass (≲10 ^9 M _⊙ ) and high-mass (≳10 ^10 M _⊙ ) end, by as much as ∼0.5 dex. This underestimation is primarily due the limited ability of fixed-form SFH parameterization to simultaneously capture the chaotic nature of star formation in these galaxies.

We report the discovery of a low-mass z = 5.200 ± 0.002 galaxy that is in the process of ceasing its star formation. The galaxy, MACS0417-z5BBG, is multiply imaged with magnification factors ∼40 by the galaxy cluster MACS J0417.5-1154, observed as part of the CAnadian NIRISS Unbiased Cluster Survey (CANUCS). Using observations of MACS0417-z5BBG with a JWST/NIRSpec Prism spectrum and NIRCam imaging, we investigate the mechanism responsible for the cessation of star formation of the galaxy and speculate about possibilities for its future. Using spectrophotometric fitting, we find a remarkably low stellar mass of ${M}_{* }=4.3{\pm }_{0.8}^{0.9}\times {10}^{7}{M}_{\odot }$ , less than 1% of the characteristic stellar mass at z ∼ 5. We measure a delensed rest-UV half-light radius in the source plane of $30{\pm }_{5}^{7}$ pc and measure a star formation rate from H α of $0.14{\pm }_{0.12}^{0.17}$ M _⊙ yr ^−1 . We find that under the assumption of a double power-law star formation history, MACS0417-z5BBG has seen a recent rise in star formation, peaking ∼10–30 Myr ago and declining precipitously since then. Together, these measurements reveal a low-mass, extremely compact galaxy which is in the process of ceasing star formation. We investigate the possibilities of mechanisms that have led to the cessation of star formation in MACS0417-z5BBG, considering stellar and active galactic nuclei (AGN) feedback and environmental processes. We can likely rule out most environmental processes but leave open the possibility of a low-mass AGN that does not leave a broad-line imprint on the spectrum or that MACS0417-z5BBG could be a star-forming galaxy in the lull of a bursty star formation history.

Lyman-break galaxy (LBG) candidates at z ≳ 10 are rapidly being identified in James Webb Space Telescope (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.6 σ SCUBA-2 detection at 850 μ m 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.

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 James Webb Space Telescope (JWST) CEERS NIRCam observations. Our sample consists of 850 galaxies at z > 3 detected in both Hubble Space Telescope (HST)/WFC3 and CEERS JWST/NIRCam images, enabling a comparison of HST and JWST morphologies. We conduct a set of visual classifications, with each galaxy in the sample classified three times. We also measure quantitative morphologies across all NIRCam filters. We find that galaxies at z > 3 have a wide diversity of morphologies. Galaxies with disks make up 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 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 observational effects impact the visibility of morphological features at high redshift. On average, Spheroid-only galaxies have a higher Sérsic index, smaller size, and higher axis ratio than disk or 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.

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. Each galaxy has been imaged with the JWST Near-Infrared Camera, which provides rest-frame optical morphologies at these redshifts. We also derive stellar masses and star formation rates for each host by fitting its spectral energy distribution using a combination of galaxy and AGN templates. We find that three of the AGN hosts have spheroidal morphologies, one is a bulge-dominated disk, and one is dominated by pointlike emission. None are found to show strong morphological disturbances that might indicate a recent interaction or merger event. When compared to a sample of mass-matched inactive galaxies, we find that the AGN hosts have morphologies that are less disturbed and more bulge-dominated. Notably, all four of the resolved hosts have rest-frame optical colors consistent with a quenched or poststarburst stellar population. The presence of AGN in passively evolving galaxies at z > 3 is significant because a rapid feedback mechanism is required in most semianalytic models and cosmological simulations to explain the growing population of massive quiescent galaxies observed at these redshifts. Our findings show that AGN can continue to inject energy into these systems after their star formation is curtailed, potentially heating their halos and preventing renewed star formation. Additional observations will be needed to determine what role this feedback may play in helping to quench these systems and/or maintain their quiescent state.

We present JWST NIRSpec spectroscopy for 11 galaxy candidates with photometric redshifts of z ≃ 9 − 13 and M _UV ∈ [ −21, −18] newly identified in NIRCam images in the Cosmic Evolution Early Release Science Survey. We confirm emission line redshifts for 7 galaxies at z = 7.762–8.998 using spectra at ∼1–5 μ m either with the NIRSpec prism or its three medium-resolution ( R ∼ 1000) gratings. For z ≃ 9 photometric candidates, we achieve a high confirmation rate of ≃90%, which validates the classical dropout selection from NIRCam photometry. No robust emission lines are identified in three galaxy candidates at z > 10, where the strong [O iii ] and H β lines would be redshifted beyond the wavelength range observed by NIRSpec, and the Ly α continuum break is not detected with the sensitivity of the current data. Compared with Hubble Space Telescope-selected bright galaxies ( M _UV ≃ −22) that are similarly spectroscopically confirmed at z ≃ 8 − 9, these NIRCam-selected galaxies are characterized by lower star formation rates (SFRs; SFR ≃ 4 M _⊙ yr ^−1 ) and lower stellar masses (≃10 ^8 M _⊙ ), but with higher specific SFR (≃40 Gyr ^−1 ), higher [O iii ]+H β equivalent widths (≃1100 Å), and elevated production efficiency of ionizing photons ( $\mathrm{log}({\xi }_{\mathrm{ion}}/\mathrm{Hz}\,{\mathrm{erg}}^{-1})\simeq 25.8$ ) induced by young stellar populations (

We report the discovery of a candidate galaxy with a photo- z of z ∼ 12 in the first epoch of the James Webb Space Telescope (JWST) Cosmic Evolution Early Release Science Survey. Following conservative selection criteria, we identify a source with a robust z _phot = ${11.8}_{-0.2}^{+0.3}$ (1 σ uncertainty) with m _F200W = 27.3 and ≳7 σ detections in five filters. The source is not detected at λ < 1.4 μ m in deep imaging from both Hubble Space Telescope (HST) and JWST and has faint ∼3 σ detections in JWST F150W and HST F160W, which signal a Ly α break near the red edge of both filters, implying z ∼ 12. This object (Maisie’s Galaxy) exhibits F115W − F200W > 1.9 mag (2 σ lower limit) with a blue continuum slope, resulting in 99.6% of the photo- z probability distribution function 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 _* / M _⊙ ∼ 8.5 and is highly star-forming (log sSFR ∼ −8.2 yr ^−1 ), with a blue rest-UV color ( β ∼ −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 that smoothly decline with increasing redshift. Should follow-up spectroscopy validate this redshift, our universe was already aglow with galaxies less than 400 Myr after the Big Bang.

We investigate the connection between the regulation of star formation and the cycling of baryons both within and in and out of galaxies. We use idealized numerical simulations of Milky Way–mass galaxies, in which we vary the galaxy morphology and stellar feedback strength. By following individual gas parcels through the disk, spiral arms, and massive star-forming clumps, we quantify how gas moves through the different phases of the interstellar medium (ISM) and forms stars. We show that the residence time of gas in the dense ISM phase ( τ _SF ), the nature of spiral arms, and the clump properties depend on both the galaxy morphology and stellar feedback. We quantify signatures of the baryon cycle within galaxies using the temporal and spatial power spectrum density (PSD) of the star formation rate (SFR). Stronger stellar feedback leads to more bursty star formation while the correlation timescale of the SFH is longer, because stronger feedback dissolves the dense ISM phase, leading to a more homogeneous ISM and a decrease in τ _SF . The bulge strength has a similar effect: the deep gravitational potential in a bulge-dominant galaxy imposes a strong shear force that breaks apart gas clumps in the ISM; this subsequently inhibits the fragmentation of gas and therefore the star formation in the disk, leading to a decrease in the spatial power on scales of ∼1 kpc. We conclude that measurements of the temporal and spatial PSD of the SFR can provide constraints on the baryon cycle and the star formation process.

The study of galaxy evolution hinges on our ability to interpret multiwavelength galaxy observations in terms of their physical properties. To do this, we rely on spectral energy distribution (SED) models, which allow us to infer physical parameters from spectrophotometric data. In recent years, thanks to wide and deep multiwave band galaxy surveys, the volume of high-quality data have significantly increased. Alongside the increased data, algorithms performing SED fitting have improved, including better modeling prescriptions, newer templates, and more extensive sampling in wavelength space. We present a comprehensive analysis of different SED-fitting codes including their methods and output with the aim of measuring the uncertainties caused by the modeling assumptions. We apply 14 of the most commonly used SED-fitting codes on samples from the CANDELS photometric catalogs at z ∼ 1 and z ∼ 3. We find agreement on the stellar mass, while we observe some discrepancies in the star formation rate (SFR) and dust-attenuation results. To explore the differences and biases among the codes, we explore the impact of the various modeling assumptions as they are set in the codes (e.g., star formation histories, nebular, dust and active galactic nucleus models) on the derived stellar masses, SFRs, and A _V values. We then assess the difference among the codes on the SFR–stellar mass relation and we measure the contribution to the uncertainties by the modeling choices (i.e., the modeling uncertainties) in stellar mass (∼0.1 dex), SFR (∼0.3 dex), and dust attenuation (∼0.3 mag). Finally, we present some resources summarizing best practices in SED fitting.

We study the spatially resolved stellar populations of 444 galaxies at 0.3 < z < 6.0 in two clusters (WHL 0137–08 and MACS 0647+70) and a blank field, combining imaging data from the Hubble Space Telescope and JWST to perform spatially resolved spectral energy distribution (SED) modeling using piXedfit . The high spatial resolution of the imaging data combined with magnification from gravitational lensing in the cluster fields allows us to resolve a large fraction of our galaxies (109) to subkiloparsec scales. At redshifts around cosmic noon and higher (2.5 ≲ z ≲ 6.0), we find mass-doubling times to be independent of radius, inferred from flat specific star formation rate (sSFR) radial profiles and similarities between the half-mass and half-SFR radii. At lower redshifts (1.5 ≲ z ≲ 2.5), a significant fraction of our star-forming galaxies shows evidence for nuclear starbursts, inferred from a centrally elevated sSFR and a much smaller half-SFR radius compared to the half-mass radius. At later epochs, we find more galaxies suppress star formation in their centers but are still actively forming stars in the disk. Overall, these trends point toward a picture of inside-out galaxy growth consistent with theoretical models and simulations. We also observe a tight relationship between the central mass surface density and global stellar mass with ∼0.38 dex scatter. Our analysis demonstrates the potential of spatially resolved SED analysis with JWST data. Future analysis with larger samples will be able to further explore the assembly of galaxy mass and the growth of their structures.

Despite the importance of feedback from active galactic nuclei (AGNs) in models of galaxy evolution, observational constraints on the influence of AGN feedback on star formation remain weak. To this end, we have compared the star formation trends of 279 low-redshift AGN galaxies with 558 inactive control galaxies using integral field unit spectroscopy from the Sloan Digital Sky Survey-IV Mapping Nearby Galaxies at Apache Point Observatory survey. With a Gaussian-process-based methodology, we reconstruct nonparametric star formation histories in spatially resolved spaxels covering the face of each galaxy. Based on the galaxy-wide star formation rates (SFRs) alone, we find no obvious signatures of AGN feedback. However, the AGN galaxies have significantly suppressed central (kiloparsec-scale) SFRs, lying up to a factor of 2 below those of the control galaxies, providing direct observational evidence of AGN feedback suppressing star formation. The suppression of central SFRs in the AGN galaxies began in the central regions ∼6 Gyr ago (redshift z ∼ 0.7), taking place over a few gigayears. A small subset of the AGN galaxies were rapidly driven to quiescence shortly before being observed (in the last 500 Myr), potentially indicating instances of AGN-driven feedback. More frequently, however, star formation continues in the AGN galaxies, with suppression primarily in the central regions. This is suggestive of a picture in which integrated (gigayear-timescale) AGN feedback can significantly affect central star formation, but may be inefficient in driving galaxy-wide quenching in low-redshift galaxies, instead leaving them in the green valley.

The UltraViolet imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey Fields (UVCANDELS) program provides Hubble Space Telescope (HST)/UVIS F275W imaging for four CANDELS fields. We combine this UV imaging with existing HST/near-IR grism spectroscopy from 3D-HST+AGHAST to directly compare the resolved rest-frame UV and H α emission for a sample of 979 galaxies at 0.7 < z < 1.5, spanning a range in stellar mass of 10 ^8−11.5 M _⊙ . Using a stacking analysis, we perform a resolved comparison between homogenized maps of rest-UV and H α to compute the average UV-to-H α luminosity ratio (an indicator of burstiness in star formation) as a function of galactocentric radius. We find that galaxies below stellar mass of ∼10 ^9.5 M _⊙ , at all radii, have a UV-to-H α ratio higher than the equilibrium value expected from constant star formation, indicating a significant contribution from bursty star formation. Even for galaxies with stellar mass ≳10 ^9.5 M _⊙ , the UV-to-H α ratio is elevated toward their outskirts ( R / R _eff > 1.5), suggesting that bursty star formation is likely prevalent in the outskirts of even the most massive galaxies, but is likely overshadowed by their brighter cores. Furthermore, we present the UV-to-H α ratio as a function of galaxy surface brightness, a proxy for stellar mass surface density, and find that regions below ∼10 ^7.5 M _⊙ kpc ^−2 are consistent with bursty star formation, regardless of their galaxy stellar mass, potentially suggesting that local star formation is independent of global galaxy properties at the smallest scales. Last, we find galaxies at z > 1.1 to have bursty star formation, regardless of radius or surface brightness.