Search

The James Webb Space Telescope is revolutionising our ability to understand the host galaxies and local environments of high-z quasars. Here we obtain a comprehensive understanding of the host galaxy of the z=7.08 quasar J1120+0641 by combining NIRSpec integral field spectroscopy with NIRCam photometry of the host continuum emission. Our emission line maps reveal that this quasar host is undergoing a merger with a bright companion galaxy. The quasar host and the companion have similar dynamical masses of $\sim10^{10}M_\odot$, suggesting that this is a major galaxy interaction. Through detailed quasar subtraction and SED fitting using the NIRCam data, we obtain an estimate of the host stellar mass of $M_{\ast}\simeq2.6\times10^9M_\odot$, with $M_{*}\simeq5.0\times10^9M_\odot$ for the companion galaxy. Using the H$\beta$ Balmer line we estimate a virial black hole mass of $M_{\rm{BH}}\simeq1.4\times10^9 M_\odot$. Thus, J1120+0641 has an extreme black hole - stellar mass ratio of $M_{\rm{BH}}/M_\ast\simeq0.54$, which is ~3 dex larger than expected by the local scaling relations between black hole and stellar mass. J1120+0641 is powered by an overmassive black hole with the highest reported black hole - stellar mass ratio, in a quasar host that is currently undergoing a major merger -- these new insights highlight the power of JWST for measuring and understanding these extreme first quasars., Comment: Submitted to A&A, comments welcome

We present a sample of 1956 individual stellar clumps at redshift 0.73, reflecting the rapid increase towards the peak of the cosmic star formation history. We conclude that the results achieved over the studied redshift range are in agreement with expectation of in-situ clump formation scenario from large-scale disk fragmentation., Comment: Submitted to MNRAS, comments welcome

Dwarf galaxies hold the key to crucial frontiers of astrophysics, however, their faintness renders spectroscopy challenging. Here we present the JWST Cycle 2 survey, All the Little Things (ALT, PID 3516), which is designed to seek late-forming Pop III stars and the drivers of reionization at $z\sim6-7$. ALT has acquired the deepest NIRCam grism spectroscopy yet (7-27 hr), at JWST's most sensitive wavelengths (3-4 $\mu$m), covering the powerful lensing cluster Abell 2744. Over the same 30 arcmin$^2$, ALT's ultra-deep F070W+F090W imaging ($\sim$30 mag) enables selection of very faint sources at $z>6$. We demonstrate the success of ALT's novel ``butterfly" mosaic to solve spectral confusion and contamination, and introduce the ``Allegro" method for emission line identification. By collecting spectra for every source in the field of view, ALT has measured precise ($R\sim1600$) redshifts for 1630 sources at $z=0.2-8.5$. This includes one of the largest samples of distant dwarf galaxies: [1015, 475, 50] sources less massive than the SMC, Fornax, and Sculptor with $\log(M_{*}/M_{\odot})<$[8.5, 7.5, 6.5]. We showcase ALT's discovery space with: (i) spatially resolved spectra of lensed clumps in galaxies as faint as $M_{\rm{UV}}\sim-15$; (ii) large-scale clustering -- overdensities at $z$=[2.50, 2.58, 3.97, 4.30, 5.66, 5.77, 6.33] hosting massive galaxies with striking Balmer breaks; (iii) small-scale clustering -- a system of satellites around a Milky Way analog at $z\sim6$; (iv) spectroscopically confirmed multiple images that help constrain the lensing model underlying all science in this legacy field; (v) sensitive star-formation maps based on dust-insensitive tracers such as Pa$\alpha$; (vi) direct spectroscopic discovery of rare sources such as AGN with ionized outflows. These results provide a powerful proof of concept for how grism surveys maximize the potential of strong lensing fields., Comment: Submitted to the Open Journal of Astrophysics. Figs. 6 and 10 illustrate the quality of the spectra and imaging, while Fig. 12 summarizes the yield of the survey. Comments warmly welcomed and greatly appreciated. The ALT DR1 catalog is available at https://zenodo.org/records/13871850

The H{\alpha} nebular emission line is an optimal tracer for recent star formation in galaxies. With the advent of JWST, this line has recently become observable at z>3 for the first time. We present a catalog of 1013 H{\alpha} emitters at 3.73 obtained based purely on spectroscopic data, robustly tracing galaxy star formation rates (SFRs) beyond the peak of the cosmic star formation history. We compare our results with theoretical predictions from three different simulations and find good agreement at z~4-6. The UV LFs of this spectroscopically-confirmed sample are in good agreement with pre-JWST measurements obtained with photometrically-selected objects. Finally, we derive SFR functions and integrate these to compute the evolution of the cosmic star-formation rate densities across z~4-6, finding values in good agreement with recent UV estimates from Lyman-break galaxies, which imply a continuous decrease in SFR density by a factor of 3x over z~4 to z~6. Our work shows the power of NIRCam grism observations to efficiently provide new tests for early galaxy formation models based on emission line statistics., Comment: 17 pages, 14 figures

We present the Red Unknowns: Bright Infrared Extragalactic Survey (RUBIES), providing JWST/NIRSpec spectroscopy of red sources selected across ~150 arcmin$^2$ from public JWST/NIRCam imaging in the UDS and EGS fields. RUBIES novel observing strategy offers a well-quantified selection function: the survey is optimised to reach high (>70%) completeness for bright and red (F150W-F444W>2) sources that are very rare. To place these rare sources in context, we simultaneously observe a reference sample of the 23$ using only the G395M disperser. The RUBIES data reveal a highly diverse population of red sources that span a broad redshift range ($z_{spec}\sim1-9$), with photometric redshift scatter and outlier fraction that are 3 times higher than for similarly bright sources that are less red. This diversity is not apparent from the photometric SEDs. Only spectroscopy reveals that the SEDs encompass a mixture of galaxies with dust-obscured star formation, extreme line emission, a lack of star formation indicating early quenching, and luminous active galactic nuclei. As a first demonstration of our broader selection function we compare the stellar masses and rest-frame U-V colours of the red sources and our reference sample: red sources are typically more massive ($M_*\sim10^{10-11.5} M_\odot$) across all redshifts. However, we find that the most massive systems span a wide range in U-V colour. We describe our data reduction procedure and data quality, and publicly release the reduced RUBIES data and vetted spectroscopic redshifts of the first half of the survey through the DJA., Comment: 21 pages, 16 figures; submitted to A&A

Our understanding of the assembly timeline of the Milky Way has been transforming along with the dramatic increase in astrometric and spectroscopic data available over the past several years. Many substructures in chemo-dynamical space have been discovered and identified as the remnants of various galactic mergers. To investigate the timeline of these mergers we select main sequence turn off & subgiant stars (MSTOs) from the H3 survey, finding members in seven metal poor components of the halo: GSE, the Helmi Streams, Thamnos, Sequoia, Wukong/LMS-1, Arjuna, and I'itoi. We also select out the metal poor in situ disk to facilitate comparison to the evolution of the Milky Way itself at these early epochs. We fit individual isochrone ages to the MSTOs in each of these substructures and use the resulting age distributions to infer simple star formation histories. For GSE we resolve an extended star formation history that truncates $\approx10$ Gyr ago, as well as a clear age -- metallicity relation. From this age distribution and measured star formation history we infer that GSE merged with the Milky Way at a time $9.5-10.2$ Gyr ago, in agreement with previous estimates. We infer that the other mergers occurred at various times ranging from $9-13$ Gyr ago, and that the metal poor component of the disk built up within only a few billion years. These results reinforce the emerging picture that both the disk and halo of the Milky Way experienced a rapid assembly., Comment: 27 pages, 17 figures

We report the spectroscopic discovery of a massive quiescent galaxy at $z_{\rm spec}=7.29\pm0.01$, just $\sim700\,$Myr after the Big Bang. RUBIES-UDS-QG-z7 was selected from public JWST/NIRCam and MIRI imaging from the PRIMER survey and observed with JWST/NIRSpec as part of RUBIES. The NIRSpec/PRISM spectrum reveals one of the strongest Balmer breaks observed thus far at $z>6$, no emission lines, but tentative Balmer and Ca absorption features, as well as a Lyman break. Simultaneous modeling of the NIRSpec/PRISM spectrum and NIRCam and MIRI photometry (spanning $0.9-18\,\mu$m) shows that the galaxy formed a stellar mass of log$(M_*/M_\odot)=10.23^{+0.04}_{-0.04}$ in a rapid $\sim 100-200\,$Myr burst of star formation at $z\sim8-9$, and ceased forming stars by $z\sim8$ resulting in $\log \rm{sSFR/yr}^{-1}<-10$. We measure a small physical size of $209_{-24}^{+33}\,{\rm pc}$, which implies a high stellar mass surface density within the effective radius of $\log(\Sigma_{*,\rm e}/{\rm M_\odot\,kpc}^{-2})=10.85_{-0.12}^{+0.11}$ comparable to the densities measured in quiescent galaxies at $z\sim2-5$. The 3D stellar mass density profile of RUBIES-UDS-QG-z7 is remarkably similar to the central densities of local massive ellipticals, suggesting that at least some of their cores may have already been in place at $z>7$. The discovery of RUBIES-UDS-QG-z7 has strong implications for galaxy formation models: the estimated number density of quiescent galaxies at $z\sim7$ is $>100\times$ larger than predicted from any model to date, indicating that quiescent galaxies have formed earlier than previously expected., Comment: 21 pages, 7 figures, submitted to ApJ

Galaxy clusters are the most massive, gravitationally-bound structures in the Universe, emerging through hierarchical structure formation of large-scale dark matter and baryon overdensities. Early galaxy ``proto-clusters'' are believed to be important physical drivers of the overall cosmic star-formation rate density and serve as ``hotspots'' for the reionization of the intergalactic medium. Our understanding of the formation of these structures at the earliest cosmic epochs is, however, limited to sparse observations of their galaxy members, or based on phenomenological models and cosmological simulations. Here we report the detection of a massive neutral, atomic hydrogen (HI) gas reservoir permeating a galaxy proto-cluster at redshift $z=5.4$, observed one billion years after the Big Bang. The presence of this cold gas is revealed by strong damped Lyman-$\alpha$ absorption features observed in several background galaxy spectra taken with JWST/NIRSpec in close on-sky projection. While overall the sightlines probe a large range in HI column densities, $N_{\rm HI} = 10^{21.7}-10^{23.5}$ cm$^{-2}$, they are similar across nearby sightlines, demonstrating that they probe the same dense, neutral gas. This observation of a massive, large-scale overdensity of cold neutral gas challenges current large-scale cosmological simulations and has strong implications for the reionization topology of the Universe., Comment: Submitted

JWST has revealed a large population of ultra-violet (UV)-bright galaxies at $z\gtrsim 10$ and possibly overly massive galaxies at $z\gtrsim 7$, challenging standard galaxy formation models in the $\Lambda$CDM cosmology. We use an empirical galaxy formation model to explore the potential of alleviating these tensions through an Early Dark Energy (EDE) model, originally proposed to solve the Hubble tension. Our benchmark model demonstrates excellent agreement with the UV luminosity functions (UVLFs) at $4\lesssim z \lesssim10$ in both $\Lambda$CDM and EDE cosmologies. In the EDE cosmology, the UVLF measurements at $z\simeq 12$ based on spectroscopically confirmed galaxies exhibit no tension with the benchmark model. Photometric constraints at $12 \lesssim z\lesssim 16$ can be fully explained within EDE via either moderately increased star formation efficiencies ($\epsilon_{\ast}\sim 3-10\%$ at $M_{\rm halo}\sim 10^{10.5}\,{\rm M}_\odot$) or enhanced UV variabilities ($\sigma_{\rm UV}\sim 0.8-1.3$ mag at $M_{\rm halo}\sim 10^{10.5}\,{\rm M}_\odot$) that are within the scatter of hydrodynamical simulation predictions. A similar agreement is difficult to achieve in $\Lambda$CDM, especially at $z\gtrsim 14$, where the required $\sigma_{\rm UV}$ exceeds the maximum value seen in simulations. Furthermore, the implausibly large cosmic stellar mass densities inferred from some JWST observations are no longer in tension with cosmology when the EDE is considered. Our findings highlight EDE as an intriguing unified solution to a fundamental problem in cosmology and the recent tensions raised by JWST observations. Data at the highest redshifts reached by JWST ($z \sim 14-16$) will be crucial for differentiating modified galaxy formation physics from new cosmological physics., Comment: Minor changes to references and adding a few new lines to compare in plots. Submitted to MNRAS

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

The infall of the Large Magellanic Cloud (LMC) is predicted to displace the inner Milky Way (MW), imprinting an apparent 'reflex motion' on the observed velocities of distant halo stars. We construct the largest all-sky spectroscopic dataset of luminous red giant stars from $50-160$ kpc, including a new survey of the southern celestial hemisphere. We fit the full 6D kinematics of our data to measure the amplitude and direction of the inner MW's motion towards the outer halo. The observed velocity grows with distance such that, relative to halo stars at $100$ kpc, the inner MW is lurching at $\approx 40$ km s$^{-1}$ towards a recent location along the LMC's past orbit. Our measurements align with N-body simulations of the halo's response to a $1.8 \times 10^{11} M_\odot$ LMC on first infall, suggesting that the LMC is at least 15% as massive as the MW. Our findings highlight the dramatic disequilibrium of the MW outskirts, and will enable more accurate measurements of the total mass of our Galaxy., Comment: 25 pages, 15 figures. Submitted to ApJ

With stunning clarity, JWST has revealed the Universe's first billion years. The scientific community is analyzing a wealth of JWST imaging and spectroscopic data from that era, and is in the process of rewriting the astronomy textbooks. Here, 1.5 years into the JWST science mission, we provide a snapshot of the great progress made towards understanding the initial chapters of our cosmic history. We highlight discoveries and breakthroughs, topics and issues that are not yet understood, and questions that will be addressed in the coming years, as JWST continues its revolutionary observations of the Early Universe. While this compendium is written by a small number of authors, invited to ISSI Bern in March 2024 as part of the 2024 ISSI Breakthrough Workshop, we acknowledge the work of a large community that is advancing our collective understanding of the evolution of the Early Universe., Comment: review article written by the attendees of the 2024 ISSI breakthrough workshop "The first billion year of the Universe", submitted. Comments welcome

The identification of red, apparently massive galaxies at $z>7$ in early JWST photometry suggests a strongly accelerated timeline compared to standard models of galaxy growth. A major uncertainty in the interpretation is whether the red colors are caused by evolved stellar populations, dust, or other effects such as emission lines or AGN. Here we show that three of the massive galaxy candidates at $z=6.7-8.4$ have prominent Balmer breaks in JWST/NIRSpec spectroscopy from the RUBIES program. The Balmer breaks demonstrate unambiguously that stellar emission dominates at $\lambda_{\rm rest} = 0.4\,\mu$m, and require formation histories extending hundreds of Myr into the past in galaxies only 600--800 Myr after the Big Bang. Two of the three galaxies also show broad Balmer lines, with H$\beta$ FWHM $>2500~{\rm km\,s^{-1}}$, suggesting that dust-reddened AGN contribute to, or even dominate, the SEDs of these galaxies at $\lambda_{\rm rest}\gtrsim 0.6\,\mu$m. All three galaxies have relatively narrow [O III] lines, seemingly ruling out a high-mass interpretation if the lines arise in dynamically-relaxed, inclined disks. Yet, the inferred masses also remain highly uncertain. We model the high-quality spectra using Prospector to decompose the continuum into stellar and AGN components, and explore limiting cases in stellar/AGN contribution. This produces a wide range of possible stellar masses, spanning $M_\star \sim 10^9 - 10^{11}\,{\rm M_{\odot}}$. Nevertheless, all fits suggest a very early and rapid formation, most of which follow with a truncation in star formation. Potential origins and evolutionary tracks for these objects are discussed, from the cores of massive galaxies to low-mass galaxies with over-massive black holes. Intriguingly, we find all of these explanations to be incomplete; deeper and redder data are needed to understand the physics of these systems., Comment: Accepted for publication in ApJL; 20 pages, 10 figures, 1 table

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

We present the first results on the spatial distribution of star formation in 454 star-forming galaxies at $4.8

In this paper, we describe the "Medium Bands, Mega Science" JWST Cycle 2 survey (JWST-GO-4111) and demonstrate the power of these data to reveal both the spatially-integrated and spatially-resolved properties of galaxies from the local universe to the era of cosmic dawn. Executed in November 2023, MegaScience obtained ~30 arcmin^2 of deep multiband NIRCam imaging centered on the z~0.3 Abell 2744 cluster, including eleven medium-band filters and the two shortest-wavelength broad-band filters, F070W and F090W. Together, MegaScience and the UNCOVER Cycle 1 treasury program provide a complete set of deep (~28-30 mag) images in all NIRCam medium- and broad-band filters. This unique dataset allows us to precisely constrain photometric redshifts, map stellar populations and dust attenuation for large samples of distant galaxies, and examine the connection between galaxy structures and formation histories. MegaScience also includes ~17 arcmin^2 of NIRISS parallel imaging in two broad-band and four medium-band filters from 0.9-4.8um, expanding the footprint where robust spectral energy distribution (SED) fitting is possible. We provide example SEDs and multi-band cutouts at a variety of redshifts, and use a catalog of JWST spectroscopic redshifts to show that MegaScience improves both the scatter and catastrophic outlier rate of photometric redshifts by factors of 2-3. Additionally, we demonstrate the spatially-resolved science enabled by MegaScience by presenting maps of the [OIII] line emission and continuum emission in three spectroscopically-confirmed z>6 galaxies. We show that line emission in reionization-era galaxies can be clumpy, extended, and spatially offset from continuum emission, implying that galaxy assembly histories are complex even at these early epochs. We publicly release fully reduced mosaics and photometric catalogs for both the NIRCam primary and NIRISS parallel fields., Comment: 21 pages, 10 figures. Fully reduced imaging, photometric catalogs, and photometric redshift fits publicly available at https://jwst-uncover.github.io/megascience/

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

Neutron star mergers (NSMs) produce copious amounts of heavy r-process elements after a time delayed inspiral process. Once NSMs are present in a galaxy, r-process elements, such as Eu, are expected to significantly increase with time. Yet, there has been limited observational data in support of Eu increasing within Local Group galaxies. We have obtained high-resolution Magellan/MIKE observations of 43 metal-poor stars in the Gaia-Sausage/Enceladus tidally disrupted galaxy with $-2.5 < \rm{[Fe/H]} < -1$. For the first time, we find a clear rise in [Eu/Mg] with increasing [Mg/H] within one galaxy. We use a chemical evolution model to study how such a rise can result from the interplay of prompt and delayed r-process enrichment events. Delayed r-process sources are required to explain the rise and subsequent leveling off of [Eu/Mg] in this disrupted galaxy. However, the rise may be explained by delayed r-process sources with either short ($\sim 10$ Myr) or long ($\sim 500$ Myr) minimum delay times. Future studies on the nature of r-process sources and their enrichment processes in the GSE will require additional stars in the GSE at even lower metallicities than the present study., Comment: 8 pages, 3 figures. Submitted to ApJ. Comments are welcomed!

The increasingly neutral intergalactic gas at $z>6$ impacts the Lyman-$\alpha$ flux observed from galaxies. One luminous galaxy, COLA1, stands out because of its unique double-peaked Ly$\alpha$ line at $z=6.6$, unseen in any simulation of reionization. Here we present JWST/NIRCam wide-field slitless spectroscopy in a 21 arcmin$^2$ field centered on COLA1. We find 141 galaxies spectroscopically-selected through the [OIII]($\lambda4969,5008$) doublet at $5.35M_{\rm UV}>-19.5$)., Comment: 15 pages, 13 figures. Accepted for publication in A&A. Main figs. 1, 2, 10 (SFR surface density), 11 (spectra and positions of galaxies in the bubble). Corrected error in escape fraction estimate from SFR surface density

Within the established framework of structure formation, galaxies start as systems of low stellar mass and gradually grow into far more massive galaxies. The existence of massive galaxies in the first billion years of the Universe, suggested by recent observations, appears to challenge this model, as such galaxies would require highly efficient conversion of baryons into stars. An even greater challenge in this epoch is the existence of massive galaxies that have already ceased forming stars. However, robust detections of early massive quiescent galaxies have been challenging due to the coarse wavelength sampling of photometric surveys. Here we report the spectroscopic confirmation with the James Webb Space Telescope of the quiescent galaxy RUBIES-EGS-QG-1 at redshift $z=4.90$, 1.2 billion years after the Big Bang. Deep stellar absorption features in the spectrum reveal that the galaxy's stellar mass of $10^{11}\,M_\odot$, corroborated by the mass implied by its gas kinematics, formed in a short $200\,$Myr burst of star formation, after which star formation activity dropped rapidly and persistently. According to current galaxy formation models, systems with such rapid stellar mass growth and early quenching are too rare to plausibly occur in the small area probed spectroscopically with JWST. Instead, the discovery of RUBIES-EGS-QG-1 implies that early massive quiescent galaxies can be quenched earlier or exhaust gas available for star formation more efficiently than currently assumed., Comment: 9 figures

JWST is producing high-quality rest-frame optical and UV spectra of faint galaxies at $z>4$ for the first time, challenging models of galaxy and stellar populations. One galaxy recently observed at $z=5.943$, GS9422, has nebular line and UV continuum emission that appears to require a high ionizing photon production efficiency. This has been explained with an exotic stellar initial mass function (IMF), 10-30x more top-heavy than a Salpeter IMF (Cameron et al. 2023). Here we suggest an alternate explanation to this exotic IMF. We use a new flexible neural net emulator for CLOUDY, Cue, to infer the shape of the ionizing spectrum directly from the observed emission line fluxes. By describing the ionizing spectrum with a piece-wise power-law, Cue is agnostic to the source of the ionizing photons. Cue finds that the ionizing radiation from GS9422 can be approximated by a double power law characterized by $\frac{Q_\mathrm{HeII}}{Q_\mathrm{H}} = -1.5$, which can be interpreted as a combination of young, metal-poor stars and a low-luminosity active galactic nucleus (AGN) with $F_{\nu} \propto \lambda ^ {2}$ in a 65%/35% ratio. This suggests a significantly lower nebular continuum contribution to the observed UV flux (24%) than a top-heavy IMF ($\gtrsim80$%), and hence, necessitates a damped Lyman-$\alpha$ absorber (DLA) to explain the continuum turnover bluewards of $\sim1400$ Angstrom. While current data cannot rule out either scenario, given the immense impact the proposed top-heavy IMF would have on models of galaxy formation, it is important to propose viable alternative explanations and to further investigate the nature of peculiar high-z nebular emitters., Comment: 14 pages, 5 figures, accepted for publication in ApJL

Recent observations from the EIGER JWST program have measured for the first time the quasar-galaxy cross-correlation function at $z\approx6$. The auto-correlation function of faint $z\approx6$ quasars was also recently estimated. These measurements provide key insights into the properties of quasars and galaxies at high redshift and their relation with the host dark matter halos. In this work, we interpret these data building upon an empirical quasar population model that has been applied successfully to quasar clustering and demographic measurements at $z\approx2-4$. We make use of a new, large-volume N-body simulation with more than a trillion particles, FLAMINGO-10k, to model quasars and galaxies simultaneously. We successfully reproduce observations of $z\approx6$ quasars and galaxies (i.e., their clustering properties and luminosity functions), and infer key quantities such as their luminosity-halo mass relation, the mass function of their host halos, and their duty cycle/occupation fraction. Our key findings are: (i) quasars reside on average in $\approx10^{12.5}\,{\rm M}_\odot$ halos (corresponding to $\approx5\sigma$ fluctuations in the initial conditions of the linear density field), but the distribution of host halo masses is quite broad; (ii) the duty cycle of (UV-bright) quasar activity is relatively low ($\approx1\%$); (iii) galaxies (that are bright in [OIII]) live in much smaller halos ($\approx10^{10.9}\,{\rm M}_\odot$) and have a larger duty cycle (occupation fraction) of $\approx13\%$. Finally, we focus on the inferred properties of quasars and present a homogeneous analysis of their evolution with redshift. The picture that emerges reveals a strong evolution of the host halo mass and duty cycle of quasars at $z\approx2-6$, and calls for new investigations of the role of quasar activity across cosmic time., Comment: Accepted for publication in MNRAS; 15 pages, 4 figures

Combining the public JWST/NIRCam imaging programs CEERS, PRIMER and JADES, spanning a total area of $\sim500\,{\rm arcmin}^2$, we obtain a sample of $>$30,000 galaxies at $z_{\rm phot}\sim4-9$ that allows us to perform a complete, rest-optical selected census of the galaxy population at $z>3$. Comparing the stellar mass $M_*$ and the UV-slope $\beta$ distributions between JWST- and HST-selected samples, we generally find very good agreement and no significant biases. Nevertheless, JWST enables us to probe a new population of UV-red galaxies that was missing from previous HST-based Lyman Break Galaxy (LBG) samples. We measure galaxy stellar mass functions (SMFs) at $z\sim4-9$ down to limiting masses of $10^{7.5}-10^{8.5}\,{\rm M_\odot}$, finding steep low mass slopes over the entire redshift range, reaching values of $\alpha\approx-2$ at $z\gtrsim6$. At the high-mass end, UV-red galaxies dominate at least out to $z\sim6$. The implied redshift evolution of the SMF suggests a rapid build-up of massive dust-obscured or quiescent galaxies from $z\sim6$ to $z\sim4$ as well as an enhanced efficiency of star formation towards earlier times ($z\gtrsim6$). Finally, we show that the galaxy mass density grows by a factor $\sim20\times$ from $z\sim9$ to $z\sim4$. Our results emphasize the importance of rest-frame optically-selected samples in inferring accurate distributions of physical properties and studying the mass build-up of galaxies in the first 1.5 Gyr of cosmic history., Comment: 31 pages, 17 figures, published in MNRAS

We expect luminous ($M_{1450}\lesssim-26.5$) high-redshift quasars to trace the highest density peaks in the early universe. Here, we present observations of four $z\gtrsim6$ quasar fields using JWST/NIRCam in imaging and widefield slitless spectroscopy mode and report a wide range in the number of detected [OIII]-emitting galaxies in the quasars' environments, ranging between a density enhancement of $\delta\approx65$ within a $2$ cMpc radius - one of the largest proto-clusters during the Epoch of Reionization discovered to date - to a density contrast consistent with zero, indicating the presence of a UV-luminous quasar in a region comparable to the average density of the universe. By measuring the two-point cross-correlation function of quasars and their surrounding galaxies, as well as the galaxy auto-correlation function, we infer a correlation length of quasars at $\langle z\rangle=6.25$ of $r_0^{\rm QQ}=22.0^{+3.0}_{-2.9}~{\rm cMpc}\,h^{-1}$, while we obtain a correlation length of the [OIII]-emitting galaxies of $r_0^{\rm GG}=4.1\pm0.3~{\rm cMpc}\,h^{-1}$. By comparing the correlation functions to dark-matter-only simulations we estimate the minimum mass of the quasars' host dark matter halos to be $\log_{10}(M_{\rm halo, min}/M_\odot)=12.43^{+0.13}_{-0.15}$ (and $\log_{10}(M_{\rm halo, min}^{\rm [OIII]}/M_\odot) = 10.56^{+0.05}_{-0.03}$ for the [OIII]-emitters), indicating that (a) luminous quasars do not necessarily reside within the most overdense regions in the early universe, and that (b) the UV-luminous duty cycle of quasar activity at these redshifts is $f_{\rm duty}\ll1$. Such short quasar activity timescales challenge our understanding of early supermassive black hole growth and provide evidence for highly dust-obscured growth phases or episodic, radiatively inefficient accretion rates., Comment: ApJ accepted

The JWST discovery of ``little red dots'' (LRDs) is reshaping our picture of the early Universe, yet the physical mechanisms driving their compact size and UV-optical colors remain elusive. Here we report an unusually bright LRD ($z=3.1$) observed as part of the RUBIES program. This LRD exhibits broad emission lines (FWHM $\sim4000$km/s), a blue UV continuum, a clear Balmer break and a red continuum sampled out to rest 4 $\mu$m with MIRI. We develop a new joint galaxy and AGN model within the Prospector Bayesian inference framework and perform spectrophotometric modeling using NIRCam, MIRI, and NIRSpec/Prism observations. Our fiducial model reveals a $M_*\sim 10^9M_\odot$ galaxy alongside a dust-reddened AGN driving the optical emission. Explaining the rest-frame optical color as a reddened AGN requires $A_{\rm v}\gtrsim4$, suggesting that a great majority of the accretion disk energy is re-radiated as dust emission. Yet despite clear AGN signatures, we find a surprising lack of hot torus emission, which implies that either the dust emission in this object must be cold, or the red continuum must instead be driven by a massive, evolved stellar population of the host galaxy -- seemingly inconsistent with the high EW broad lines (H$\alpha$ EW $\sim800$\AA). The widths and luminosities of Pa$\beta$, Pa$\delta$, Pa$\gamma$, and H$\alpha$ imply a modest black hole mass of $M_{\rm BH}\sim10^8M_\odot$. Additionally, we identify a narrow blue-shifted HeI absorption in G395M spectra, signaling an ionized outflow with kinetic energy up to $\sim1$\% the luminosity of the AGN. The low redshift of RUBIES-BLAGN-1 combined with the depth and richness of the JWST imaging and spectroscopic observations provide a unique opportunity to build a physical model for these so-far mysterious LRDs, which may prove to be a crucial phase in the early formation of massive galaxies and their supermassive black holes., Comment: 20 pages, 8 figures, 4 tables

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

We illustrate the formation and evolution of the Milky Way over cosmic time, utilizing a sample of 10 million red giant stars with full chemodynamical information, including metallicities and $\alpha$-abundances from low-resolution Gaia XP spectra. The evolution of angular momentum as a function of metallicity - a rough proxy for stellar age, particularly for high-[$\alpha$/Fe] stars - displays three distinct phases: the disordered and chaotic protogalaxy, the kinematically-hot old disk, and the kinematically-cold young disk. The old high-$\alpha$ disk starts at [Fe/H] $\approx -1.0$, 'spinning up' from the nascent protogalaxy, and then exhibits a smooth 'cooldown' toward more ordered and circular orbits at higher metallicities. The young low-$\alpha$ disk is kinematically cold throughout its metallicity range, with its observed properties modulated by a strong radial gradient. We interpret these trends using Milky Way analogs from the TNG50 cosmological simulation, identifying one that closely matches the kinematic evolution of our Galaxy. This halo's protogalaxy spins up into a relatively thin and misaligned high-$\alpha$ disk at early times, which is subsequently heated and torqued by a major gas-rich merger. The merger contributes a large amount of low-metallicity gas and angular momentum, from which the kinematically cold low-$\alpha$ stellar disk is subsequently born. This simulated history parallels several observed features of the Milky Way, particularly the decisive 'GSE' merger that likely occurred at $z \approx 2$. Our results provide an all-sky perspective on the emerging picture of our Galaxy's three-phase formation, impelled by the three physical mechanisms of spinup, merger, and cooldown., Comment: 21 pages, 12 figures, submitted to ApJ

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

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

JWST is revealing a new population of dust-reddened broad-line active galactic nuclei (AGN) at redshifts $z\gtrsim5$. Here we present deep NIRSpec/Prism spectroscopy from the Cycle 1 Treasury program UNCOVER of 15 AGN candidates selected to be compact, with red continua in the rest-frame optical but with blue slopes in the UV. From NIRCam photometry alone, they could have been dominated by dusty star formation or AGN. Here we show that the majority of the compact red sources in UNCOVER are dust-reddened AGN: $60\%$ show definitive evidence for broad-line H$\alpha$ with FWHM$\, >2000$ km/s, for $20\%$ current data are inconclusive, and $20\%$ are brown dwarf stars. We propose an updated photometric criterion to select red $z>5$ AGN that excludes brown dwarfs and is expected to yield $>80\%$ AGN. Remarkably, among all $z_{\rm phot}>5$ galaxies with F277W$-$F444W$>1$ in UNCOVER at least $33\%$ are AGN regardless of compactness, climbing to at least $80\%$ AGN for sources with F277W$-$F444W$>1.6$. The confirmed AGN have black hole masses of $10^7-10^9$ M$_{\odot}$. While their UV-luminosities ($-16>M_{\rm UV}>-20$ AB mag) are low compared to UV-selected AGN at these epochs, consistent with percent-level scattered AGN light or low levels of unobscured star formation, the inferred bolometric luminosities are typical of $10^7-10^9$ M$_{\odot}$ black holes radiating at $\sim 10-40\%$ of Eddington. The number densities are surprisingly high at $\sim10^{-5}$ Mpc$^{-3}$ mag$^{-1}$, 100 times more common than the faintest UV-selected quasars, while accounting for $\sim1\%$ of the UV-selected galaxies. While their UV-faintness suggest they may not contribute strongly to reionization, their ubiquity poses challenges to models of black hole growth., Comment: 23 pages, 9 figures, 5 tables, submitted to ApJ

We report {\em JWST}/NIRCam measurements of quasar host galaxy emissions and supermassive black hole (SMBH) masses for six quasars at $5.9

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 present the first detailed chemical-abundance analysis of stars from the dwarf-galaxy stellar stream Wukong/LMS-1 covering a wide metallicity range ($-3.5 < \rm[Fe/H] \lesssim -1.3$). We find abundance patterns that are effectively indistinguishable from the bulk of Indus and Jhelum, a pair of smaller stellar streams proposed to be dynamically associated with Wukong/LMS-1. We confirmed a carbon-enhanced metal-poor star ($\rm[C/Fe] > +0.7$ and $\rm[Fe/H] \sim -2.9$) in Wukong/LMS-1 with strong enhancements in Sr, Y, and Zr, which is peculiar given its solar-level [Ba/Fe]. Wukong/LMS-1 stars have high abundances of $\alpha$ elements up to $\rm[Fe/H] \gtrsim -2$, which is expected for relatively massive dwarfs. Towards the high-metallicity end, Wukong/LMS-1 becomes $\alpha$-poor, revealing that it probably experienced fairly standard chemical evolution. We identified a pair of N- and Na-rich stars in Wukong/LMS-1, reminiscent of multiple populations in globular clusters. This indicates that this dwarf galaxy contained at least one globular cluster that was completely disrupted in addition to two intact ones previously known to be associated with Wukong/LMS-1, which is possibly connected to similar evidence found in Indus. From these $\geq$3 globular clusters, we estimate the total mass of Wukong/LMS-1 to be ${\approx}10^{10} M_\odot$, representing ${\sim}1$% of the present-day Milky Way. Finally, the [Eu/Mg] ratio in Wukong/LMS-1 continuously increases with metallicity, making this the first example of a dwarf galaxy where the production of $r$-process elements is clearly dominated by delayed sources, presumably neutron-star mergers., Comment: Accepted to MNRAS. New version fixes abundance uncertainties, which significantly affects elements like Al and N. Use new version of abundance files for correct error bars

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

We report new observations of the cool diffuse gas around 29, $2.3$9. MgII absorption kinematics are shifted from the systemic redshift of host galaxies with a median absolute velocity of 135 km/s and standard deviation of 85 km/s. The high kinematic offset and large radial separation ($R> 1.3 R_{200}$), suggest that five out of the seven MgII absorption systems are gravitationally not bound to the galaxies. In contrast, most cool circumgalactic media at $z<1$ are gravitationally bound. The high incidence of unbound MgII gas in this work suggests that towards the end of reionization, galaxy halos are in a state of remarkable disequilibrium, and are highly efficient in enriching the intergalactic medium. Two strongest MgII absorption systems are detected at $z\sim$ 4.22 and 4.5, the former associated with a merging galaxy system and the latter associated with three kinematically close galaxies. Both these galaxies reside in local galaxy over-densities, indicating the presence of cool MgII absorption in two "proto-groups" at $z>4$., Comment: 23 pages, 12 figures, 4 tables, Accepted for publication in ApJ

The Magellanic Stream (MS) - an enormous ribbon of gas spanning $140^\circ$ of the southern sky trailing the Magellanic Clouds - has been exquisitely mapped in the five decades since its discovery. However, despite concerted efforts, no stellar counterpart to the MS has been conclusively identified. This stellar stream would reveal the distance and 6D kinematics of the MS, constraining its formation and the past orbital history of the Clouds. We have been conducting a spectroscopic survey of the most distant and luminous red giant stars in the Galactic outskirts. From this dataset, we have discovered a prominent population of 13 stars matching the extreme angular momentum of the Clouds, spanning up to $100^\circ$ along the MS at distances of $60-120$ kpc. Furthermore, these kinemetically-selected stars lie along a [$\alpha$/Fe]-deficient track in chemical space from $-2.5 < \mathrm{[Fe/H]} < -0.5$, consistent with their formation in the Clouds themselves. We identify these stars as high-confidence members of the Magellanic Stellar Stream. Half of these stars are metal-rich and closely follow the gaseous MS, whereas the other half are more scattered and metal-poor. We argue that the metal-rich stream is the recently-formed tidal counterpart to the MS, and speculate that the metal-poor population was thrown out of the SMC outskirts during an earlier interaction between the Clouds. The Magellanic Stellar Stream provides a strong set of constraints - distances, 6D kinematics, and birth locations - that will guide future simulations towards unveiling the detailed history of the Clouds., Comment: 21 pages, 12 figures. Submitted to ApJ

The James Webb Space Telescope (JWST) is revolutionizing our knowledge of $z>5$ galaxies and their actively accreting black holes. Using the JWST Cycle 1 Treasury program Ultradeep NIRSpec and NIRCam ObserVations before the Epoch of Reionization (UNCOVER) in the lensing field Abell 2744, we report the identification of a sample of little red dots at $3 < z_{\rm{phot}} < 7$ that likely contain highly-reddened accreting supermassive black holes. Using a NIRCam-only selection to F444W$<27.7$ mag, we find 26 sources over the $\sim45$ arcmin$^{2}$ field that are blue in F115W$-$F200W$\sim0$ (or $\beta_{\rm UV}\sim-2.0$ for $f_{\lambda} \propto \lambda^\beta$), red in F200W$-$F444W = $1-4$ ($\beta_{\rm opt} \sim +2.0$), and are dominated by a point-source like central component. Of the 20 sources with deep ALMA 1.2-mm coverage, none are detected individually or in a stack. For the majority of the sample, SED fits to the JWST+ALMA observations prefer models with hot dust rather than obscured star-formation to reproduce the red NIRCam colors and ALMA 1.2-mm non-detections. While compact dusty star formation can not be ruled out, the combination of extremely small sizes ($\langle r_e \rangle\approx50$ pc after correction for magnification), red rest-frame optical slopes, and hot dust can by explained by reddened broad-line active galactic nuclei (AGNs). Our targets have faint $M_{\rm 1450} \approx -14\ \, {\rm to} -18$ mag but inferred bolometric luminosities of $L_{\rm bol} = 10^{43}-10^{46}$ erg/s, reflecting their obscured nature. If the candidates are confirmed as AGNs with upcoming UNCOVER spectroscopy, then we have found an abundant population of reddened luminous AGN that are at least ten times more numerous than UV-luminous AGN at the same intrinsic bolometric luminosity., Comment: submitted to ApJ

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)

We report a spectroscopic search for hypervelocity white dwarfs (WDs) that are runaways from Type Ia supernovae (SNe Ia) and related thermonuclear explosions. Candidates are selected from Gaia data with high tangential velocities and blue colors. We find six new runaways, including four stars with radial velocities (RVs) $>1000\,\rm km\,s^{-1}$ and total space velocities $\gtrsim 1300\,\rm km\,s^{-1}$. These are most likely the surviving donors from double-degenerate binaries in which the other WD exploded. The other two objects have lower minimum velocities, $\gtrsim 600\,\rm km\,s^{-1}$, and may have formed through a different mechanism, such as pure deflagration of a WD in a Type Iax supernova. The four fastest stars are hotter and smaller than the previously known "D$^6$ stars," with effective temperatures ranging from $\sim$20,000 to $\sim$130,000 K and radii of $\sim 0.02-0.10\,R_{\odot}$. Three of these have carbon-dominated atmospheres, and one has a helium-dominated atmosphere. Two stars have RVs of $-1694$ and $-2285\rm \,km\,s^{-1}$ -- the fastest systemic stellar RVs ever measured. Their inferred birth velocities, $\sim 2200-2500\,\rm km\,s^{-1}$, imply that both WDs in the progenitor binary had masses $>1.0\,M_{\odot}$. The high observed velocities suggest that a dominant fraction of the observed hypervelocity WD population comes from double-degenerate binaries whose total mass significantly exceeds the Chandrasekhar limit. However, the two nearest and faintest D$^6$ stars have the lowest velocities and masses, suggesting that observational selection effects favor rarer, higher-mass stars. A significant population of fainter low-mass runaways may still await discovery. We infer a birth rate of D$^6$ stars that is consistent with the SN Ia rate. The birth rate is poorly constrained, however, because the luminosities and lifetimes of $\rm D^6$ stars are uncertain., Comment: 26 pages, 17 figures. Accepted to OJA

The onset of galaxy formation is thought to be initiated by the infall of neutral, pristine gas onto the first protogalactic halos. However, direct constraints on the abundance of neutral atomic hydrogen (HI) in galaxies have been difficult to obtain at early cosmic times. Here we present spectroscopic observations with JWST of three galaxies at redshifts $z=8.8 - 11.4$, about $400-600$ Myr after the Big Bang, that show strong damped Lyman-$\alpha$ absorption ($N_{\rm HI} > 10^{22}$ cm$^{-2}$) from HI in their local surroundings, an order of magnitude in excess of the Lyman-$\alpha$ absorption caused by the neutral intergalactic medium at these redshifts. Consequently, these early galaxies cannot be contributing significantly to reionization, at least at their current evolutionary stages. Simulations of galaxy formation show that such massive gas reservoirs surrounding young galaxies so early in the history of the universe is a signature of galaxy formation in progress., Comment: Submitted

We have derived luminosity functions, and set constraints on the UV luminosity and SFR density from z~17 to z~8, using the three most-studied JWST/NIRCam data sets, the SMACS0723, GLASS Parallel, and CEERS fields. We first used our own selections on two independent reductions of these datasets using the latest calibrations. 18 z~8, 12 z~10, 5 z~13, and 1 z~17 candidate galaxies are identified over these fields in our primary reductions, with a similar number of candidates in our secondary reductions. We then use these two reductions, applying a quantitative discriminator, to segregate the full set of z>~8 candidates reported over these fields from the literature, into three different samples, ``robust,'' ``solid,'' and ``possible''. Using all of these samples we then derive UV LF and luminosity density results at $z\geq8$, finding substantial differences. For example, including the full set of ``solid'' and ``possible'' z>~12 candidates from the literature, we find UV luminosity densities which are ~7x and ~20x higher than relying on the ``robust'' candidates alone. These results indicate the evolution of the UV LF and luminosity densities at z>~8 is still extremely uncertain, emphasizing the need for spectroscopy and deeper NIRCam+optical imaging to obtain reliable results. Nonetheless, even with the very conservative ``robust'' approach to selections, both from our own and those of other studies, we find the luminosity density from luminous (M(UV)<-19) galaxies to be ~2x larger than is easily achievable using constant star-formation efficiency models, similar to what other early JWST results have suggested., Comment: 26 pages, 10 figures, 15 tables, accepted for publication in MNRAS, figure 7 shows the key results

The majority of the Milky Way's stellar halo consists of debris from our Galaxy's last major merger, the Gaia-Sausage-Enceladus (GSE). In the past few years, stars from GSE have been kinematically and chemically studied in the inner $30$ kpc of our Galaxy. However, simulations predict that accreted debris could lie at greater distances, forming substructures in the outer halo. Here we derive metallicities and distances using Gaia DR3 XP spectra for an all-sky sample of luminous red giant stars, and map the outer halo with kinematics and metallicities out to $100$ kpc. We obtain follow-up spectra of stars in two strong overdensities - including the previously identified Outer Virgo Overdensity - and find them to be relatively metal-rich and on predominantly retrograde orbits, matching predictions from simulations of the GSE merger. We argue that these are apocentric shells of GSE debris, forming $60-90$ kpc counterparts to the $15-20$ kpc shells that are known to dominate the inner stellar halo. Extending our search across the sky with literature radial velocities, we find evidence for a coherent stream of retrograde stars encircling the Milky Way from $50-100$ kpc, in the same plane as the Sagittarius stream but moving in the opposite direction. These are the first discoveries of distant and structured imprints from the GSE merger, cementing the picture of an inclined and retrograde collision that built up our Galaxy's stellar halo., Comment: 22 pages, 10 figures. Accepted to ApJ

The Milky Way's stellar halo, which extends to $>100$ kpc, encodes the evolutionary history of our Galaxy. However, most studies of the halo to date have been limited to within a few kpc of the Sun. Here, we characterize differences between this local halo and the stellar halo in its entirety. We construct a composite stellar halo model by combining observationally motivated N-body simulations of the Milky Way's nine most massive disrupted dwarf galaxies that account for almost all of the mass in the halo. We find that (1) the representation by mass of different dwarf galaxies in the local halo compared to the whole halo can be significantly overestimated (e.g., the Helmi Streams) or underestimated (e.g., Cetus) and (2) properties of the overall halo (e.g., net rotation) inferred via orbit integration of local halo stars are significantly biased, because e.g., highly retrograde debris from Gaia-Sausage-Enceladus is missing from the local halo. Therefore, extrapolations from the local to the global halo should be treated with caution. From analysis of a sample of 11 MW-like simulated halos, we identify a population of recently accreted ($\lesssim5$ Gyrs) and disrupted galaxies on high angular momenta orbits that are entirely missing from local samples, and awaiting discovery in the outer halo. Our results motivate the need for surveys of halo stars extending to the Galaxy's virial radius., Comment: 12 pages, 6 figures. Submitted to ApJ. Comments warmly welcomed

We present the first constraints on the prevalence of z>10 galaxies in the Hubble Ultra Deep Field (HUDF) leveraging new NIRCam observations from JEMS (JWST Extragalactic Medium-band Survey). These NIRCam observations probe redward of 1.6$\mu$m, beyond the wavelength limit of HST, allowing us to search for galaxies to z>10. These observations indicate that the highest redshift candidate identified in the HUDF09 data with HST, UDFj-39546284, has a redshift of z>11.5, as had been suggested in analyses of the HUDF12/XDF data. This has now been confirmed with JWST NIRSpec. This source is thus the most distant galaxy discovered by HST in its >30 years of operation. Additionally, we identify nine other z~8-13 candidate galaxies over the HUDF, two of which are new discoveries that appear to lie at z~11-12. We use these results to characterize the evolution of the UV luminosity function (LF) from z~15 to z~8.7. While our LF results at z~8.7 and z~10.5 are consistent with previous findings over the HUDF, our new LF estimates at z~12.6 are higher than other results in the literature, potentially pointing to a milder evolution in the UV luminosity density from z~12.6. We emphasize that our LF results are uncertain given the small number of z~12.6 sources and limited volume probed. The new NIRCam data also indicate that the faint z~8-13 galaxies in the HUDF/XDF show blue UV-continuum slopes beta~-2.7, high specific star formation rates ~24.5 Gyr**-1, and high EW (~1300A) [OIII]+H$\beta$ emission, with two z~8.5 sources showing [OIII]+H$\beta$ EWs of ~2300 A., Comment: 21 pages, 14 figures, 7 tables, resubmitted to MNRAS, after replying to the referee report

We model the stellar abundances and ages of two disrupted dwarf galaxies in the Milky Way stellar halo: Gaia-Sausage Enceladus (GSE) and Wukong/LMS-1. Using a statistically robust likelihood function, we fit one-zone models of galactic chemical evolution with exponential infall histories to both systems, deriving e-folding timescales of $\tau_\text{in} = 1.01 \pm 0.13$ Gyr for GSE and $\tau_\text{in} = 3.08^{+3.19}_{-1.16}$ Gyr for Wukong/LMS-1. GSE formed stars for $\tau_\text{tot} = 5.40^{+0.32}_{-0.31}$ Gyr, sustaining star formation for $\sim$$1.5 - 2$ Gyr after its first infall into the Milky Way $\sim$10 Gyr ago. Our fit suggests that star formation lasted for $\tau_\text{tot} = 3.36^{+0.55}_{-0.47}$ Gyr in Wukong/LMS-1, though our sample does not contain any age measurements. The differences in evolutionary parameters between the two are qualitatively consistent with trends with stellar mass $M_\star$ predicted by simulations and semi-analytic models of galaxy formation. Our fitting method is based only on poisson sampling from an evolutionary track and requires no binning of the data. We demonstrate its accuracy by testing against mock data, showing that it accurately recovers the input model across a broad range of sample sizes ($20 \leq N \leq 2000$) and measurement uncertainties ($0.01 \leq \sigma_\text{[$\alpha$/Fe]}, \sigma_\text{[Fe/H]} \leq 0.5$; $0.02 \leq \sigma_{\log_{10}(\text{age})} \leq 1$). Our inferred values of the outflow mass-loading factor reasonably match $\eta \propto M_\star^{-1/3}$ as predicted by galactic wind models. Due to the generic nature of our derivation, this likelihood function should be applicable to one-zone models of any parametrization and easily extensible to other astrophysical models which predict tracks in some observed space., Comment: 25 pages; 13 figures; submitted to MNRAS; comments welcome

Massive disk galaxies like our Milky Way should host an ancient, metal-poor, and centrally concentrated stellar population. This population reflects the star formation and enrichment in the few most massive progenitor components that coalesced at high redshift to form the proto-Galaxy. While metal-poor stars are known to reside in the inner few kiloparsecs of our Galaxy, current data do not yet provide a comprehensive picture of such a metal-poor "heart" of the Milky Way. We use information from Gaia DR3, especially the XP spectra, to construct a sample of 2 million bright (BP $<15.5$ mag) giant stars within $30^\circ$ of the Galactic Center with robust [M/H] estimates, $\delta$ [M/H] $\lesssim 0.1$. For most sample members we can calculate orbits based on Gaia RVS velocities and astrometry. This sample reveals an extensive, ancient, and metal-poor population that includes $\sim 18,000$ stars with $-2.7<$ [M/H] $<-1.5$, representing a stellar mass of $\gtrsim 5\times 10^7$ M$_\odot$. The spatial distribution of these [M/H] $<-1.5$ stars has a Gaussian extent of only $\sigma_{\mathrm{R_{GC}}} \sim 2.7$ kpc around the Galactic center, with most of these orbits being confined to the inner Galaxy. At high orbital eccentricities, there is clear evidence for accreted halo stars in their pericentral orbit phase. Stars with [M/H] $< -2$ show no net rotation, whereas those with [M/H] $\sim -1$ are rotation dominated. Most of the tightly bound stars show $[\alpha/\text{Fe}]$-enhancement and [Al/Fe]-[Mn/Fe] abundance patterns expected for an origin in the more massive portions of the proto-Galaxy. These central, metal-poor stars most likely predate the oldest part of the disk ($\tau_{\text{age}}\approx 12.5$ Gyrs), which implies that they formed at $z\gtrsim 5$, forging the proto-Milky Way., Comment: 21 pages, 9 figures, submitted to ApJ

Modern Galactic surveys have revealed an ancient merger that dominates the stellar halo of our Galaxy (\textit{Gaia}-Sausage-Enceladus, GSE). Using chemical abundances and kinematics from the H3 Survey, we identify 5559 halo stars from this merger in the radial range $r_{\text{Gal}}=6-60\text{ kpc}$. We forward model the full selection function of H3 to infer the density profile of this accreted component of the stellar halo. We consider a general ellipsoid with principal axes allowed to rotate with respect to the Galactocentric axes, coupled with a multiply-broken power law. The best-fit model is a triaxial ellipsoid (axes ratios 10:8:7) tilted $25^\circ$ above the Galactic plane towards the Sun and a doubly-broken power law with breaking radii at 12 kpc and 28 kpc. This result resolves the long-standing dichotomy in literature values of the halo breaking radius, being at either $\sim15\text{ kpc}$ or $\sim30\text{ kpc}$ assuming a singly-broken power law. N-body simulations suggest that the breaking radii are connected to apocenter pile-ups of stellar orbits, and so the observed double-break provides new insight into the initial conditions and evolution of the GSE merger. Furthermore, the tilt and triaxiality of the stellar halo could imply that a fraction of the underlying dark matter halo is also tilted and triaxial. This has important implications for dynamical mass modeling of the Galaxy as well as direct dark matter detection experiments., Comment: Submitted to ApJ; comments warmly welcome