Search

It remains a mystery when our Milky Way first formed a stellar disk component that survived and maintained its disk structure from subsequent galaxy mergers. We present a study of the age-dependent structure and star formation rate of the Milky Way's disk using high-alpha stars with significant orbital angular momentum that have precise age determinations. Our results show that the radial scale length is nearly independent of age, while the vertical scale height experienced dramatic evolution. A disk-like geometry presents even for populations older than 13 Gyr, with the scale height-to-length ratio dropping below 0.5 for populations younger than 12.5 Gyr. We dub the oldest population that has maintained a disk geometry - apparently formed over 13 Gyr ago - PanGu. With an estimated present-day stellar mass of $2 \times 10^9$ $M_\odot$, PanGu is presumed to be a major stellar component of our Galaxy in the earliest epoch. The total present-day stellar mass of the whole high-alpha disk is $2 \times 10^{10}$ $M_\odot$, mostly formed during a distinct star formation rate peak of 11 $M_\odot$ per year around 11 Gyrs ago. A comparison with Milky Way analogs in the TNG50 simulations implies that our Galaxy has experienced an exceptionally quiescent dynamical history, even before the Gaia-Enceladus merger., Comment: Published in Nature Astronomy on 10 October, 2024; This manuscript is the accepted version

In the third APOKASC catalog, we present data for the complete sample of 15,808 evolved stars with APOGEE spectroscopic parameters and Kepler asteroseismology. We used ten independent asteroseismic analysis techniques and anchor our system on fundamental radii derived from Gaia $L$ and spectroscopic $T_{\rm eff}$. We provide evolutionary state, asteroseismic surface gravity, mass, radius, age, and the spectroscopic and asteroseismic measurements used to derive them for 12,418 stars. This includes 10,036 exceptionally precise measurements, with median fractional uncertainties in \nmax, \dnu, mass, radius and age of 0.6\%, 0.6\%, 3.8\%, 1.8\%, and 11.1\% respectively. We provide more limited data for 1,624 additional stars which either have lower quality data or are outside of our primary calibration domain. Using lower red giant branch (RGB) stars, we find a median age for the chemical thick disk of $9.14 \pm 0.05 ({\rm ran}) \pm 0.9 ({\rm sys})$ Gyr with an age dispersion of 1.1 Gyr, consistent with our error model. We calibrate our red clump (RC) mass loss to derive an age consistent with the lower RGB and provide asymptotic GB and RGB ages for luminous stars. We also find a sharp upper age boundary in the chemical thin disk. We find that scaling relations are precise and accurate on the lower RGB and RC, but they become more model dependent for more luminous giants and break down at the tip of the RGB. We recommend the usage of multiple methods, calibration to a fundamental scale, and the usage of stellar models to interpret frequency spacings., Comment: 43 pages, 25 figures, submitted ApJSupp. Comments welcome. Data tables available on request from pinsonneault.1@osu.edu

Ly$\alpha$ is the transition to the ground state from the first excited state of hydrogen (the most common element). Resonant scattering of this line by neutral hydrogen greatly impedes its emergence from galaxies, so the fraction of galaxies which show Ly$\alpha$ is a tracer of the neutral fraction of the intergalactic medium (IGM), and thus the history of reionization. In previous works, we used early JWST/NIRSpec data from the JWST Advanced Deep Extragalactic Survey (JADES) to classify and characterise Ly$\alpha$ emitting galaxies (LAEs). This survey is now approaching completion, and the current sample is nearly an order of magnitude larger. From a sample of 784 galaxies in JADES at $4.0

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

Early JWST studies found an apparent population of massive, compact galaxies at redshifts $z\gtrsim7$. Recently three of these galaxies were shown to have prominent Balmer breaks, demonstrating that their light at $\lambda_{\rm rest} \sim 3500$ $\r{A}$ is dominated by a stellar population that is relatively old ($\sim$200 Myr). All three also have broad H$\beta$ emission with $\sigma > 1000 \,\rm km s^{-1}$, a common feature of such `little red dots'. From S\'ersic profile fits to the NIRCam images in F200W we find that the stellar light of galaxies is extremely compact: the galaxies have half-light radii of $r_{\rm e}\sim$ 100 pc, in the regime of ultra compact dwarfs in the nearby Universe. Their masses are uncertain, as they depend on the contribution of possible light from an AGN to the flux at $\lambda_{\rm rest}>5000$ $\r{A}$. If the AGN contribution is low beyond the Balmer break region, the masses are $M_* \sim 10^{10}-10^{11}\mathrm{M}_{\odot}$, and the central densities are higher than those of any other known galaxy population by an order of magnitude. Interestingly, the implied velocity dispersions of $\sim$1500 kms$^{-1}$ are in very good agreement with the measured H$\beta$ line widths. We suggest that some of the broad lines in `little red dots' are not due to AGNs but simply reflect the kinematics of the galaxies, and speculate that the galaxies are observed in a short-lived phase where the central densities are much higher than at later times. We stress, however, that the canonical interpretation of AGNs causing the broad H$\beta$ lines also remains viable.

Stars in the Galactic disk are born on cold, nearly circular orbits with small vertical excursions. After their birth, their orbits evolve, driven by small- or large-scale perturbations in the Galactic disk's gravitational potential. Here, we study the vertical motions of young stars over their first few orbital periods, using a sample of OBA stars from \textit{Gaia} E/DR3, which includes radial velocities and ages $\tau$ from LAMOST. We constructed a parametric model for the time evolution of the stellar orbits' mean vertical actions $J_z$ as a function of Galactocentric radius, $R_{\mathrm{GC}}$. Accounting for data uncertainties, we use Markov Chain Monte Carlo (MCMC) analysis in annuli of Galactocentric radius to constrain the model parameters. Our best-fit model shows a remarkably linear increase of vertical actions with age across all Galactocentric radii examined. Orbital \textit{heating} by random scattering could offer a straightforward interpretation for this trend. However, various other dynamical aspects of the Galactic disk, such as stars being born in a warped disk, might offer alternative explanations that could be tested in the future., Comment: 10 pages, 5 figures. Draft version July 9, 2024. Typeset using LATEX default style in AASTeX631

A large population of intermediate-separation binaries, consisting of a main-sequence (MS) star and a white dwarf (WD), has recently emerged from Gaia's third data release (DR3), posing challenges to current models of binary evolution. Here we examine the $s$-process element abundances in these systems using data from GALAH DR3. Following refined sample analysis with parameter estimates based on GALAH spectra, we find a distinct locus where enhanced $s$-process elements depend on both the WD mass and metallicity, consistent with loci identified in previous asymptotic giant branch (AGB) nucleosynthesis studies with higher $s$-process yields. Notably, these enhanced abundances show no correlation with the systems' orbital parameters, supporting a history of accretion in intermediate-separation MS+WD systems. Consequently, our results form a direct observational evidence of a connection between AGB masses and $s$-process yields. We conclude by showing that the GALAH DR3 survey includes numerous Ba dwarf stars, within and beyond the mass range covered in our current sample, which can further elucidate $s$-process element distributions in MS+WD binaries., Comment: Published in ApJL; 12 pages, 10 figures

We present JWST/NIRSpec Integral Field Spectroscopy in the rest-frame optical bands of the system PJ308-21, a quasar at $z=6.2342$ caught as its host galaxy interacts with companion galaxies. We detect spatially extended emission of several emission lines (H$\alpha$, H$\beta$, [OIII], [NII], [SII], HeII), which we use to study the properties of the ionized phase of the interstellar medium: the source and hardness of the photoionizing radiation field, metallicity, dust reddening, electron density and temperature, and star formation. We also marginally detect continuum starlight emission associated with the companion sources. We find that at least two independent satellite galaxies are part of the system. While the quasar host appears highly enriched and obscured, with AGN-like photoionization conditions, the western companion shows minimal dust extinction, low metallicity ($Z\sim0.4$ Z$_\odot$), and star-formation driven photoionization. The eastern companion shows higher extinction and metallicity ($Z\sim0.8$ Z$_\odot$) compared to the western companion, and it is at least partially photoionized by the nearby quasar. We do not find any indication of AGN in the companion sources. Our study shows that while the quasar host galaxy is already very massive ($M_{\rm dyn}>10^{11}$ M$_\odot$), it is still rapidly building up by accreting two relatively massive ($M_{\rm star}\sim 10^{10}$ M$_\odot$) companion sources. This dataset showcases the power of JWST in exposing the build-up of massive galaxies in the first Gyr of the Universe., Comment: 15 pages, 16 figures. Accepted for publication in A&A

We show with Gaia XP spectroscopy that extremely metal-rich stars in the Milky Way (EMR; $[M/H]_{XP} > 0.5$) - but only those - are largely confined to a tight "knot" at the center of the Galaxy. This EMR knot is round in projection, has a fairly abrupt edge near $\sim 1.5$kpc, and is a dynamically hot system. This central knot also contains very metal-rich (VMR; $+0.2\le [M/H]_{XP} \le +0.4$) stars. However, in contrast to EMR stars, the bulk of VMR stars form an extended, highly flattened distribution in the inner Galaxy ($R_{\mathrm{GC}}\lesssim 5$ kpc). We draw on TNG50 simulations of Milky Way analogs for context and find that compact, metal-rich knots confined to $<1.5$kpc are a universal feature. In typical simulated analogs, the top 5-10% most metal-rich stars are confined to a central knot; however, in our Milky Way data this fraction is only 0.1%. Dust-penetrating wide-area near-infrared spectroscopy, such as SDSS-V, will be needed for a rigorous estimate of the fraction of stars in the Galactic EMR knot. Why in our Milky Way only EMR giants are confined to such a central knot remains to be explained. Remarkably, the central few kiloparsecs of the Milky Way harbor both the highest concentration of metal-poor stars (the `poor old heart') and almost all EMR stars. This highlights the stellar population diversity at the bottom of galactic potential wells., Comment: 11 pages, 7 figures, submitted to ApJ

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

A suite of spectroscopic surveys is producing vast sets of stellar spectra with the goal of advancing stellar physics and Galactic evolution by determining their basic physical properties. A substantial fraction of these stars are in binary systems, but almost all large-survey modeling pipelines treat them as single stars. For sets of multi-epoch spectra, spectral disentangling is a powerful technique to recover or constrain the individual components' spectra of a multiple system. So far, this approach has focused on small samples or individual objects, usually with high resolution ($R \gtrsim 10.000$) spectra and many epochs ($\gtrsim 8$). Here, we present a disentangling implementation that accounts for several aspects of few-epoch spectra from large surveys: that vast sample sizes require automatic determination of starting guesses; that some of the most extensive spectroscopic surveys have a resolution of only $\approx 2,000$; that few epochs preclude unique orbit fitting; that one needs effective regularisation of the disentangled solution to ensure resulting spectra are smooth. We describe the implementation of this code and show with simulated spectra how well spectral recovery can work for hot and cool stars at $R \approx 2000$. Moreover, we verify the code on two established binary systems, the ``Unicorn'' and ``Giraffe''. This code can serve to explore new regimes in survey disentangling in search of massive stars with massive dark companions, e.g. the $\gtrsim 200,000$ hot stars of the SDSS-V survey., Comment: Accepted to MNRAS

We forecast the expected population of active galactic nuclei (AGN) observable in the Euclid Wide Survey (EWS) and Euclid Deep Survey (EDS). Starting from an X-ray luminosity function (XLF) we generate volume-limited samples of the AGN expected in the survey footprints. Each AGN is assigned an SED appropriate for its X-ray luminosity and redshift, with perturbations sampled from empirical distributions. The photometric detectability of each AGN is assessed via mock observation of the assigned SED. We estimate 40 million AGN will be detectable in at least one band in the EWS and 0.24 million in the EDS, corresponding to surface densities of 2.8$\times$10$^{3}$ deg$^{-2}$ and 4.7$\times$10$^{3}$ deg$^{-2}$. Employing colour selection criteria on our simulated data we select a sample of 4.8$\times$10$^{6}$ (331 deg$^{-2}$) AGN in the EWS and 1.7$\times$10$^{4}$ (346 deg$^{-2}$) in the EDS, amounting to 10% and 8% of the AGN detectable in the EWS and EDS. Including ancillary Rubin/LSST bands improves the completeness and purity of AGN selection. These data roughly double the total number of selected AGN to comprise 21% and 15% of the detectable AGN in the EWS and EDS. The total expected sample of colour-selected AGN contains 6.0$\times$10$^{6}$ (74%) unobscured AGN and 2.1$\times$10$^{6}$ (26%) obscured AGN, covering $0.02 \leq z \lesssim 5.2$ and $43 \leq \log_{10} (L_{bol} / erg s^{-1}) \leq 47$. With this simple colour selection, expected surface densities are already comparable to the yield of modern X-ray and mid-infrared surveys of similar area. The relative uncertainty on our expectation for detectable AGN is 6.7% for the EWS and 12.5% for the EDS, driven by the uncertainty of the XLF., Comment: 36 pages, 21 figures, submitted to A&A

As part of the Euclid Early Release Observations (ERO) programme, we analyse deep, wide-field imaging from the VIS and NISP instruments of two Milky Way globular clusters (GCs), namely NGC 6254 (M10) and NGC 6397, to look for observational evidence of their dynamical interaction with the Milky Way. We search for such an interaction in the form of structural and morphological features in the clusters' outermost regions, which are suggestive of the development of tidal tails on scales larger than those sampled by the ERO programme. Our multi-band photometric analysis results in deep and well-behaved colour-magnitude diagrams that, in turn, enable an accurate membership selection. The surface brightness profiles built from these samples of member stars are the deepest ever obtained for these two Milky Way GCs, reaching down to $\sim30.0$ mag~arcsec$^{-2}$, which is about $1.5$ mag arcsec$^{-2}$ below the current limit. The investigation of the two-dimensional density map of NGC 6254 reveals an elongated morphology of the cluster peripheries in the direction and with the amplitude predicted by $N$-body simulations of the cluster's dynamical evolution, at high statistical significance. We interpret this as strong evidence for the first detection of tidally induced morphological distortion around this cluster. The density map of NGC 6397 reveals a slightly elliptical morphology, in agreement with previous studies, which requires further investigation on larger scales to be properly interpreted. This ERO project thus demonstrates the power of Euclid in studying the outer regions of GCs at an unprecedented level of detail, thanks to the combination of large field of view, high spatial resolution, and depth enabled by the telescope. Our results highlight the future Euclid survey as the ideal data set to investigate GC tidal tails and stellar streams., Comment: 15 pages, 18 figures. Paper accepted as part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

We provide an early assessment of the imaging capabilities of the Euclid space mission to probe deeply into nearby star-forming regions and associated very young open clusters, and in particular to check to what extent it can shed light on the new-born free-floating planet population. This paper focuses on a low-reddening region observed in just one Euclid pointing where the dust and gas has been cleared out by the hot sigma Orionis star. One late-M and six known spectroscopically confirmed L-type substellar members in the sigma Orionis cluster are used as benchmarks to provide a high-purity procedure to select new candidate members with Euclid. The exquisite angular resolution and depth delivered by the Euclid instruments allow us to focus on bona-fide point sources. A cleaned sample of sigma Orionis cluster substellar members has been produced and the initial mass function (IMF) has been estimated by combining Euclid and Gaia data. Our sigma Orionis substellar IMF is consistent with a power-law distribution with no significant steepening at the planetary-mass end. No evidence of a low-mass cutoff is found down to about 4 Jupiter masses at the young age (3 Myr) of the sigma Orionis open cluster., Comment: 13 pages, 10 figures, submitted to the A&A special issue on "Euclid on the sky"

The Euclid ERO showcase Euclid's capabilities in advance of its main mission, targeting 17 astronomical objects, from galaxy clusters, nearby galaxies, globular clusters, to star-forming regions. A total of 24 hours observing time was allocated in the early months of operation, engaging the scientific community through an early public data release. We describe the development of the ERO pipeline to create visually compelling images while simultaneously meeting the scientific demands within months of launch, leveraging a pragmatic, data-driven development strategy. The pipeline's key requirements are to preserve the image quality and to provide flux calibration and photometry for compact and extended sources. The pipeline's five pillars are: removal of instrumental signatures; astrometric calibration; photometric calibration; image stacking; and the production of science-ready catalogues for both the VIS and NISP instruments. We report a PSF with a full width at half maximum of 0.16" in the optical and 0.49" in the three NIR bands. Our VIS mean absolute flux calibration is accurate to about 1%, and 10% for NISP due to a limited calibration set; both instruments have considerable colour terms. The median depth is 25.3 and 23.2 AB mag with a SNR of 10 for galaxies, and 27.1 and 24.5 AB mag at an SNR of 5 for point sources for VIS and NISP, respectively. Euclid's ability to observe diffuse emission is exceptional due to its extended PSF nearly matching a pure diffraction halo, the best ever achieved by a wide-field, high-resolution imaging telescope. Euclid offers unparalleled capabilities for exploring the LSB Universe across all scales, also opening a new observational window in the NIR. Median surface-brightness levels of 29.9 and 28.3 AB mag per square arcsec are achieved for VIS and NISP, respectively, for detecting a 10 arcsec x 10 arcsec extended feature at the 1 sigma level., Comment: Submitted to A&A, 44 pages, 36 figures - Part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

The near-infrared calibration unit (NI-CU) on board Euclid's Near-Infrared Spectrometer and Photometer (NISP) is the first astronomical calibration lamp based on light-emitting diodes (LEDs) to be operated in space. Euclid is a mission in ESA's Cosmic Vision 2015-2025 framework, to explore the dark universe and provide a next-level characterisation of the nature of gravitation, dark matter, and dark energy. Calibrating photometric and spectrometric measurements of galaxies to better than 1.5% accuracy in a survey homogeneously mapping ~14000 deg^2 of extragalactic sky requires a very detailed characterisation of near-infrared (NIR) detector properties, as well their constant monitoring in flight. To cover two of the main contributions - relative pixel-to-pixel sensitivity and non-linearity characteristics - as well as support other calibration activities, NI-CU was designed to provide spatially approximately homogeneous (<12% variations) and temporally stable illumination (0.1%-0.2% over 1200s) over the NISP detector plane, with minimal power consumption and energy dissipation. NI-CU is covers the spectral range ~[900,1900] nm - at cryo-operating temperature - at 5 fixed independent wavelengths to capture wavelength-dependent behaviour of the detectors, with fluence over a dynamic range of >=100 from ~15 ph s^-1 pixel^-1 to >1500 ph s^-1 pixel^-1. For this functionality, NI-CU is based on LEDs. We describe the rationale behind the decision and design process, describe the challenges in sourcing the right LEDs, as well as the qualification process and lessons learned. We also provide a description of the completed NI-CU, its capabilities and performance as well as its limits. NI-CU has been integrated into NISP and the Euclid satellite, and since Euclid's launch in July 2023 has started supporting survey operations., Comment: Paper accepted for publication in A&A as part of the special issue 'Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

This paper presents the specification, design, and development of the Visible Camera (VIS) on the ESA Euclid mission. VIS is a large optical-band imager with a field of view of 0.54 deg^2 sampled at 0.1" with an array of 609 Megapixels and spatial resolution of 0.18". It will be used to survey approximately 14,000 deg^2 of extragalactic sky to measure the distortion of galaxies in the redshift range z=0.1-1.5 resulting from weak gravitational lensing, one of the two principal cosmology probes of Euclid. With photometric redshifts, the distribution of dark matter can be mapped in three dimensions, and, from how this has changed with look-back time, the nature of dark energy and theories of gravity can be constrained. The entire VIS focal plane will be transmitted to provide the largest images of the Universe from space to date, reaching m_AB>24.5 with S/N >10 in a single broad I_E~(r+i+z) band over a six year survey. The particularly challenging aspects of the instrument are the control and calibration of observational biases, which lead to stringent performance requirements and calibration regimes. With its combination of spatial resolution, calibration knowledge, depth, and area covering most of the extra-Galactic sky, VIS will also provide a legacy data set for many other fields. This paper discusses the rationale behind the VIS concept and describes the instrument design and development before reporting the pre-launch performance derived from ground calibrations and brief results from the in-orbit commissioning. VIS should reach fainter than m_AB=25 with S/N>10 for galaxies of full-width half-maximum of 0.3" in a 1.3" diameter aperture over the Wide Survey, and m_AB>26.4 for a Deep Survey that will cover more than 50 deg^2. The paper also describes how VIS works with the other Euclid components of survey, telescope, and science data processing to extract the cosmological information., Comment: Paper submitted as part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

We present the Sloan Digital Sky Survey V (SDSS-V) Local Volume Mapper (LVM). The LVM is an integral-field spectroscopic survey of the Milky Way, Magellanic Clouds, and of a sample of local volume galaxies, connecting resolved pc-scale individual sources of feedback to kpc-scale ionized interstellar medium (ISM) properties. The 4-year survey covers the southern Milky Way disk at spatial resolutions of 0.05 to 1 pc, the Magellanic Clouds at 10 pc resolution, and nearby large galaxies at larger scales totaling $>4300$ square degrees of sky, and more than 55M spectra. It utilizes a new facility of alt-alt mounted siderostats feeding 16 cm refractive telescopes, lenslet-coupled fiber-optics, and spectrographs covering 3600-9800A at R ~ 4000. The ultra-wide field IFU has a diameter of 0.5 degrees with 1801 hexagonally packed fibers of 35.3 arcsec apertures. The siderostats allow for a completely stationary fiber system, avoiding instability of the line spread function seen in traditional fiber feeds. Scientifically, LVM resolves the regions where energy, momentum, and chemical elements are injected into the ISM at the scale of gas clouds, while simultaneously charting where energy is being dissipated (via cooling, shocks, turbulence, bulk flows, etc.) to global scales. This combined local and global view enables us to constrain physical processes regulating how stellar feedback operates and couples to galactic kinematics and disk-scale structures, such as the bar and spiral arms, as well as gas in- and out-flows., Comment: 29 pages, 12 figures, accepted for publication in The Astronomical Journal

We report discovery and spectroscopic follow-up of 21 astrometric binaries containing solar-type stars and dark companions with masses near 1.4 $M_{\odot}$. The simplest interpretation is that the companions are dormant neutron stars (NSs), though ultramassive white dwarfs (WDs) and tight WD+WD binaries cannot be fully excluded. We selected targets from Gaia DR3 astrometric binary solutions in which the luminous star is on the main sequence and the dynamically-implied mass of the unseen companion is (a) more than $1.25\,M_{\odot}$ and (b) too high to be any non-degenerate star or close binary. We obtained multi-epoch radial velocities (RVs) over a period of 700 days, spanning a majority of the orbits' dynamic range in RV. The RVs broadly validate the astrometric solutions and significantly tighten constraints on companion masses. Several systems have companion masses that are unambiguously above the Chandrasekhar limit, while the rest have masses between 1.25 and 1.4 $M_{\odot}$. The orbits are significantly more eccentric at fixed period than those of typical WD + MS binaries, perhaps due to natal kicks. Metal-poor stars are overrepresented in the sample: 3 out of 21 objects (14%) have [Fe/H]$\sim-1.5$ and are on halo orbits, compared to $\sim$0.5% of the parent Gaia binary sample. The metal-poor stars are all strongly enhanced in lithium. The formation history of these objects is puzzling: it is unclear both how the binaries escaped a merger or dramatic orbital shrinkage when the NS progenitors were red supergiants, and how they remained bound when the NSs formed. Gaia has now discovered 3 black holes (BHs) in astrometric binaries with masses above 9 $M_{\odot}$, and 21 NSs with masses near $1.4\,M_{\odot}$. The lack of intermediate-mass objects in this sample is striking, supporting the existence of a BH/NS mass bimodality over 4 orders of magnitude in orbital period., Comment: 32 pages, 20 figures, accepted to OJAp

Gaia can not individually resolve very close binary systems, however, the collected data can still be used to identify them. A powerful indicator of stellar multiplicity is the sources reported Renormalized Unit Weight Error (ruwe), which effectively captures the astrometric deviations from single-source solutions. We aim to characterise the imprints left on ruwe caused by binarity. By flagging potential binary systems based on ruwe, we aim to characterise which of their properties will contribute the most to their detectability. We develop a model to estimate ruwe values for observations of Gaia sources, based on the biases to the single-source astrometric track arising from the presence of an unseen companion. Then, using the recipes from previous GaiaUnlimited selection functions, we estimate the selection probability of sources with high ruwe, and discuss what binary properties contribute to increasing the sources ruwe. We compute the maximum ruwe value which is compatible with single-source solutions as a function of their location on-sky. We see that binary systems selected as sources with a ruwe higher than this sky-varying threshold have a strong detectability window in their orbital period distribution, which peaks at periods equal to the Gaia observation time baseline. We demonstrate how our sky-varying ruwe threshold provides a more complete sample of binary systems when compared to single sky-averaged values by studying the unresolved binary population in the Gaia Catalogue of Nearby Stars. We provide the code and tools used in this study, as well as the sky-varying ruwe threshold through the GaiaUnlimited Python package, Comment: 11 pages, 10 figures. Accepted for publication in A&A

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

Using $\textit{Gaia}$ DR3 astrometry and spectroscopy, we study two new substructures in the orbit-metallicity space of the inner Milky Way: $\textit{Shakti}$ and $\textit{Shiva}$. They were identified as two confined, high-contrast overdensities in the $(L_z, E)$ distribution of bright ($G<16$) and metal-poor ($-2.5<\rm{[M/H]}<-1.0$) stars. Both have stellar masses of $M_\star \gtrsim 10^7M_\odot$, and are distributed on prograde orbits inside the Solar circle in the Galaxy. Both structures have an orbit-space distribution that points towards an $\textit{accreted}$ origin, however, their abundance patterns -- from APOGEE -- are such that are conventionally attributed to an $\textit{in situ}$ population. These seemingly contradictory diagnostics could be reconciled if we interpret the abundances [Mg/Fe], [Al/Fe], [Mg/Mn] $\textit{vs.}$ [Fe/H] distribution of their member stars merely as a sign of rapid enrichment. This would then suggest one of two scenarios. Either these prograde substructures were created by some form of resonant orbit trapping of the field stars by the rotating bar; a plausible scenario proposed by Dillamore et al. (2023). Or, $\textit{Shakti}$ and $\textit{Shiva}$ were proto-galactic fragments that formed stars rapidly and coalesced early, akin to the constituents of the $\textit{Poor Old Heart}$ of the Milky Way; just less deep in the Galactic potential and still discernible in orbit space., Comment: Published in ApJ. The paper contains 18 pages, 8 figures, 2 Tables

We offer a new way to look at the Large Magellanic Cloud through stellar mono-abundance and mono-age-mono-abundance maps. These maps are based on $\gtrsim 500\,000$ member stars with photo-spectroscopic [M/H] and age estimates from Gaia DR3 data, and they are the first area-complete, metallicity- and age-differentiated stellar maps of any disk galaxy. Azimuthally averaged, these maps reveal a surprisingly simple picture of the Milky Way's largest satellite galaxy. For any [M/H] below -0.1 dex, the LMC's radial profile is well described by a simple exponential, but with a scale length that steadily shrinks towards higher metallicities, from nearly 2.3~kpc at [M/H]$=-1.8$ to only 0.75~kpc at [M/H]$=-0.25$. The prominence of the bar decreases dramatically with [M/H], making it barely discernible at [M/H]$\lesssim -1.5$. Yet, even for metal-rich populations, the bar has little impact on the azimuthally averaged profile of the mono-abundance components. Including ages, we find that the scale length is a greater function of age than of metallicity, with younger populations far more centrally concentrated. At old ages, the scale length decreases with increasing metallicity; at young ages, the scale-length is independent of metallicity. These findings provide quantitative support for a scenario where the LMC built its stellar structure effectively outside in., Comment: Submitted to ApJ; constructive comments would be appreciated

The Near-infrared Spectrograph (NIRSpec) on the James Webb Space Telescope is uniquely suited to studying galaxies in the distant Universe with its combination of multi-object capabilities and sensitivity over a large range in wavelength (0.6-5.3 microns). Here we present the NIRSpec Wide survey, part of the NIRSpec Instrument Science Team's Guaranteed Time Observations, using NIRSpec's microshutter array to obtain spectra of more than 3200 galaxies at $z>1$ at both low- and high-resolution ($R\approx100$ and 2700) for a total of 105 hours. With 31 pointings covering $\approx$320 arcmin$^2$ across the five CANDELS fields with exquisite ancillary photometry from the Hubble Space Telescope, the NIRSpec Wide survey represents a fast and efficient way of using JWST to probe galaxies in the early Universe. Pointing centers are determined to maximize the observability of the rarest, high-value sources. Subsequently, the microshutter configurations are optimized to observe the maximum number of "census" galaxies with a selection function based primarily on HST/F160W magnitude, photometric/slitless grism redshift, and predicted \ha\ flux tracing the bulk of the galaxy population at cosmic noon ($z_{\rm med}=2.0$). We present details on the survey strategy, the target selection, an outline of the motivating science cases, and discuss upcoming public data releases to the community., Comment: Published in A&A. Data for pointings in AEGIS available for download at https://archive.stsci.edu/hlsp/wide (https://archive.stsci.edu/hlsps/wide/download_scripts/hlsp_wide_jwst_nirspec_all_multi_v1.0_spec-download.sh)

The Galaxy inner parts are generally considered to be optically symmetric, as well as kinematically symmetric for most massive early-type galaxies. At the lower-mass end, many galaxies contain lots of small patches in their velocity maps, causing their kinematics to be nonsmooth in small scales and far from symmetry. These small patches can easily be mistaken for measurement uncertainties and have not been well discussed. We used the comparison of observations and numerical simulations to demonstrate the small patches existence beyond uncertainties. For the first time we have found that the fluctuation degrees have an approximate inverse loglinear relation with the galaxy stellar surface mass densities. This tight relation among galaxies that do not show obvious optical asymmetry that traces environmental perturbations indicates that stellar motion in galaxies has inherent asymmetry besides external environment influences. The degree of the kinetic asymmetry is closely related to and constrained by the intrinsic properties of the host galaxy., Comment: 13 pages,10 figures, accepted for publication in ApJL

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

Whereas well-studied galaxies at cosmic noon are found to be baryon-dominated within the effective radius, recent JWST observations of $z\sim6-7$ galaxies with stellar masses of only $M_*\sim10^{8-9}\,{\rm M_\odot}$ surprisingly indicate that they are dark matter-dominated within $r_{\rm e}\approx 1\,$kpc. Here, we place these high-redshift measurements in the context of the TNG50 galaxy formation simulation, by measuring the central (within $1\,$kpc) stellar, gas, and dark matter masses of galaxies in the simulation. The central baryon fraction varies strongly with galaxy stellar mass in TNG50, and this $M_*$-dependence is remarkably constant across $0

We report discovery and characterization of a main-sequence G star orbiting a dark object with mass $1.90\pm 0.04 M_{\odot}$. The system was discovered via Gaia astrometry and has an orbital period of 731 days. We obtained multi-epoch RV follow-up over a period of 639 days, allowing us to refine the Gaia orbital solution and precisely constrain the masses of both components. The luminous star is a $\gtrsim 12$ Gyr-old, low-metallicity halo star near the main-sequence turnoff ($T_{\rm eff}\approx 6000$ K; $\log(g/\left[{\rm cm\,s^{-2}}\right])\approx 4.0$; $\rm [Fe/H]\approx-1.25$; $M\approx0.79 M_{\odot}$) with a highly enhanced lithium abundance. The RV mass function sets a minimum companion mass for an edge-on orbit of $M_2 > 1.67 M_{\odot}$, well above the Chandrasekhar limit. The Gaia inclination constraint, $i=68.7\pm 1.4$ deg, then implies a companion mass of $M_2=1.90\pm0.04 M_{\odot}$. The companion is most likely a massive neutron star: the only viable alternative is two massive white dwarfs in a close binary, but this scenario is disfavored on evolutionary grounds. The system's low eccentricity ($e=0.122\pm 0.002$) disfavors dynamical formation channels and implies that the neutron star likely formed with little mass loss ($\lesssim1\,M_{\odot}$) and with a weak natal kick ($v_{\rm kick}\lesssim 20\,\rm km\,s^{-1}$). The current orbit is too small to have accommodated the neutron star progenitor as a red supergiant or super-AGB star. The simplest formation scenario -- isolated binary evolution -- requires the system to have survived unstable mass transfer and common envelope evolution with a donor-to-accretor mass ratio $>10$. The system, which we call Gaia NS1, is likely a progenitor of symbiotic X-ray binaries and long-period millisecond pulsars. Its discovery challenges binary evolution models and bodes well for Gaia's census of compact objects in wide binaries., Comment: 19 pages, 12 figures, accepted to OJAp

The spatial distribution of Galactic stars with different chemical abundances encodes information on the processes that drove the formation and evolution of the Milky Way. Survey selection functions are indispensable for analysing astronomical catalogues produced by large-scale surveys. The use of these selection functions in data modelling is more complex when data from different surveys are to be modelled simultaneously. We introduce a procedure for constructing the selection function of a sample of red clump stars that have parallaxes and elemental abundances from the Gaia mission. We separately constructed the selection function of the APOGEE DR17 red clump stars, which depends on very different observables and has a very different spatial coverage. We combined the two surveys and accounted for their joint selection function to provide strong constraints on the radial and vertical density distribution of mono-abundance populations, with Gaia offering a dense coverage of the solar neighbourhood, while APOGEE reaches larger distances near the Galactic plane. We confirm that the radial density profile steepens with increasing metallicity. The combined sample also indicates a metallicity-dependent flaring of the alpha-poor disc. We provide the code for constructing the Gaia selection function we used in this study through the GaiaUnlimited Python package., Comment: 21 pages, 18 figures, accepted for publication in Astronomy & Astrophysics

Stars formed with initial mass over 50 Msun are very rare today, but they are thought to be more common in the early universe. The fates of those early, metal-poor, massive stars are highly uncertain. Most are expected to directly collapse to black holes, while some may explode as a result of rotationally powered engines or the pair-creation instability. We present the chemical abundances of J0931+0038, a nearby low-mass star identified in early followup of SDSS-V Milky Way Mapper, which preserves the signature of unusual nucleosynthesis from a massive star in the early universe. J0931+0038 has relatively high metallicity ([Fe/H] = -1.76 +/- 0.13) but an extreme odd-even abundance pattern, with some of the lowest known abundance ratios of [N/Fe], [Na/Fe], [K/Fe], [Sc/Fe], and [Ba/Fe]. The implication is that a majority of its metals originated in a single extremely metal-poor nucleosynthetic source. An extensive search through nucleosynthesis predictions finds a clear preference for progenitors with initial mass > 50 Msun, making J0931+0038 one of the first observational constraints on nucleosynthesis in this mass range. However the full abundance pattern is not matched by any models in the literature. J0931+0038 thus presents a challenge for the next generation of nucleosynthesis models and motivates study of high-mass progenitor stars impacted by convection, rotation, jets, and/or binary companions. Though rare, more examples of unusual early nucleosynthesis in metal-poor stars should be found in upcoming large spectroscopic surveys., Comment: 11 pages + 22 page appendix, accepted to ApJL

We present the UV-to-NIR size evolution of a sample of 161 quiescent galaxies (QGs) with $M_*>10^{10}M_\odot$ over $0.510^{10.6}M_\odot$. To constrain the physical mechanisms driving the apparent size evolution, we study the relationship between $R_e$ and the formation redshift ($z_{form}$) of QGs. For lower-mass QGs, this relationship is broadly consistent with $R_e\sim(1+z_{form})^{-1}$, in line with the expectation of the progenitor effect. For higher-mass QGs, the relationship between $R_e$ and $z_{form}$ depends on stellar age. Older QGs have a steeper relationship between $R_e$ and $z_{form}$ than that expected from the progenitor effect alone, suggesting that mergers and/or post-quenching continuous gas accretion drive additional size growth in very massive systems. We find that the $z>3$ QGs in our sample are very compact, with mass surface densities $\Sigma_e\gtrsim10^{10} M_\odot/\rm{kpc}^2$, and their $R_e$ are possibly even smaller than anticipated from the size evolution measured for lower-redshift QGs. Finally, we take a close look at the structure of GS-9209, one of the earliest confirmed massive QGs at $z_{spec}\sim4.7$. From UV to NIR, GS-9209 becomes increasingly compact, and its light profile becomes more spheroidal, showing that the color gradient is already present in this earliest massive QG., Comment: 28 pages, 19 figures, submitted to ApJ

We present high-precision radial velocity (RV) observations of Gaia BH1, the nearest known black hole (BH). The system contains a solar-type G star orbiting a massive dark companion, which could be either a single BH or an inner BH + BH binary. A BH + BH binary is expected in some models where Gaia BH1 formed as a hierarchical triple, which are attractive because they avoid many of the difficulties associated with forming the system through isolated binary evolution. Our observations test the inner binary scenario. We have measured 115 precise RVs of the G star, including 40 from ESPRESSO with a precision of $3$-$5$ m s$^{-1}$, and 75 from other instruments with a typical precision of $30$-$100$ m s$^{-1}$. Our observations span $2.33$ orbits of the G star and are concentrated near a periastron passage, when perturbations due to an inner binary would be largest. The RVs are well-fit by a Keplerian two-body orbit and show no convincing evidence of an inner binary. Using REBOUND simulations of hierarchical triples with a range of inner periods, mass ratios, eccentricities, and orientations, we show that plausible inner binaries with periods $P_{\text{inner}} \gtrsim 1.5$ days would have produced larger deviations from a Keplerian orbit than observed. Binaries with $P_{\text{inner}} \lesssim 1.5$ days are consistent with the data, but these would merge within a Hubble time and would thus imply fine-tuning. We present updated parameters of Gaia BH1's orbit. The RVs yield a spectroscopic mass function $f\left(M_{\text{BH}}\right)=3.9358 \pm 0.0002\,M_{\odot}$ - about $7000\sigma$ above the $\sim2.5\,M_{\odot}$ maximum neutron star mass. Including the inclination constraint from Gaia astrometry, this implies a BH mass of $M_{\text{BH}} = 9.27 \pm 0.10 ~ M_{\odot}$., Comment: 26 pages, 14 figures, Accepted to PASP. Github repository at https://github.com/pranav-nagarajan/Gaia-BH1-Precision-RVs

The orbit distribution of young stars in the Galactic disk is highly structured, from well-defined clusters to streams of stars that may be widely dispersed across the sky, but are compact in orbital action-angle space. The age distribution of such groups can constrain the timescales over which co-natal groups of stars disperse into the `field'. Gaia data have proven powerful to identify such groups in action-angle space, but the resulting member samples are often too small and have too narrow a CMD coverage to allow robust age determinations. Here, we develop and illustrate a new approach that can estimate robust stellar population ages for such groups of stars. This first entails projecting the predetermined action-angle distribution into the 5D space of positions, parallaxes and proper motions, where much larger samples of likely members can be identified over a much wider range of the CMD. It then entails isochrone fitting that accounts for a) widely varying distances and reddenings; b) outliers and binaries; c) sparsely populated main sequence turn-offs, by incorporating the age information of the low-mass main sequence; and d) the possible presence of an intrinsic age spread in the stellar population. When we apply this approach to 92 nearby stellar groups identified in 6D orbit space, we find that they are predominately young ($\lesssim 1$ Gyr), mono-age populations. Many groups are established (known) localized clusters with possible tidal tails, others tend to be widely dispersed and manifestly unbound. This new age-dating tool offers a stringent approach to understanding on which orbits stars form in the solar neighborhood and how quickly they disperse into the field., Comment: Accepted for publication in ApJ

Measuring the mass-radius relation of individual white dwarfs is an empirically challenging task that has been performed for only a few dozen stars. We measure the white dwarf mass-radius relation using gravitational redshifts and radii of 137 white dwarfs in wide binaries with main sequence companions. We obtain the space velocities to these systems using the main sequence companion, and subtract these Doppler redshifts from the white dwarfs' apparent motions, isolating their gravitational redshifts. We use Gaia data to calculate the surface temperatures and radii of these white dwarfs, thereby deriving an empirical gravitational redshift-radius relation. This work demonstrates the utility of low-resolution Galactic surveys to measure the white dwarf equation of state. Our results are consistent with theoretical models, and represent the largest sample of individual white dwarf gravitational redshift measurements to date., Comment: 12 pages, 9 figures

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 summarize the properties and initial data release of the JADES Origins Field (JOF), which will soon be the deepest imaging field yet observed with the James Webb Space Telescope (JWST). This field falls within the GOODS-S region about 8' south-west of the Hubble Ultra Deep Field (HUDF), where it was formed initially in Cycle 1 as a parallel field of HUDF spectroscopic observations within the JWST Advanced Deep Extragalactic Survey (JADES). This imaging will be greatly extended in Cycle 2 program 3215, which will observe the JOF for 5 days in six medium-band filters, seeking robust candidates for z>15 galaxies. This program will also include ultra-deep parallel NIRSpec spectroscopy (up to 104 hours on-source, summing over the dispersion modes) on the HUDF. Cycle 3 observations from program 4540 will add 20 hours of NIRCam slitless spectroscopy to the JOF. With these three campaigns, the JOF will be observed for 380 open-shutter hours with NIRCam using 15 imaging filters and 2 grism bandpasses. Further, parts of the JOF have deep 43 hr MIRI observations in F770W. Taken together, the JOF will soon be one of the most compelling deep fields available with JWST and a powerful window into the early Universe. This paper presents the second data release from JADES, featuring the imaging and catalogs from the year 1 JOF observations., Comment: Submitted to ApJ Supplement. Images and catalogs are available at https://archive.stsci.edu/hlsp/jades . A FITSmap portal to view the images is at https://jades.idies.jhu.edu

We present AspGap, a new approach to infer stellar labels from low-resolution Gaia XP spectra, including precise [$\alpha$/M] estimates for the first time. AspGap is a neural-network based regression model trained on APOGEE spectra. In the training step, AspGap learns to use XP spectra not only to predict stellar labels but also the high-resolution APOGEE spectra that lead to the same stellar labels. The inclusion of this last model component -- dubbed the hallucinator -- creates a more physically motivated mapping and significantly improves the prediction of stellar labels in the validation, particularly of [$\alpha$/M]. For giant stars, we find cross-validated rms accuracies for Teff, log g, [M/H], [$\alpha$/M] of ~1%, 0.12 dex, 0.07 dex, 0.03 dex, respectively. We also validate our labels through comparison with external datasets and through a range of astrophysical tests that demonstrate that we are indeed determining [$\alpha$/M] from the XP spectra, rather than just inferring it indirectly from correlations with other labels. We publicly release the AspGap codebase, along with our stellar parameter catalog for all giants observed by Gaia XP. AspGap enables new insights into the formation and chemo-dynamics of our Galaxy by providing precise [$\alpha$/M] estimates for 23 million giant stars, including 12 million with radial velocities from Gaia., Comment: 29 pages, 19 figures, submitted to ApJS

We present a detailed analysis of the the structure of the Local Group flocculent spiral galaxy M33, as measured using the Panchromatic Hubble Andromeda Treasury Triangulum Extended Region (PHATTER) survey. Leveraging the multiwavelength coverage of PHATTER, we find that the oldest populations are dominated by a smooth exponential disk with two distinct spiral arms and a classical central bar $-$ completely distinct from what is seen in broadband optical imaging, and the first-ever confirmation of a bar in M33. We estimate a bar extent of $\sim$1 kpc. The two spiral arms are asymmetric in orientation and strength, and likely represent the innermost impact of the recent tidal interaction responsible for M33's warp at larger scales. The flocculent multi-armed morphology for which M33 is known is only visible in the young upper main sequence population, which closely tracks the morphology of the ISM. We investigate the stability of M33's disk, finding $Q{\sim}1$ over the majority of the disk. We fit multiple components to the old stellar density distribution and find that, when considering recent stellar kinematics, M33's bulk structure favors the inclusion of an accreted halo component, modeled as a broken power-law. The best-fit halo model has an outer power-law index of $-$3 and accurately describes observational evidence of M33's stellar halo from both resolved stellar spectroscopy in the disk and its stellar populations at large radius. Integrating this profile yields a total halo stellar mass of ${\sim}5{\times}10^8\ M_{\odot}$, giving a total stellar halo mass fraction of 16%, most of which resides in the innermost 2.5 kpc., Comment: 31 pages, 17 figures, 5 tables, Accepted for publication in ApJ

We explore the kinematic gas properties of six $5.51$ thus far. The cold gas masses implied by the observed star formation rates are $\sim 10\times$ larger than the stellar masses. We find that their ionised gas is spatially resolved by JWST, with evidence for broadened lines and spatial velocity gradients. Using a simple thin-disc model, we fit these data with a novel forward modelling software that accounts for the complex geometry, point spread function, and pixellation of the NIRSpec instrument. We find the sample to include both rotation- and dispersion-dominated structures, as we detect velocity gradients of $v(r_{\rm e})\approx100-150\,{\rm km\,s^{-1}}$, and find velocity dispersions of $\sigma_0\approx 30-70\,{\rm km\,s^{-1}}$ that are comparable to those at cosmic noon. The dynamical masses implied by these models ($M_{\rm dyn}\sim10^{9-10}\,{\rm M_\odot}$) are larger than the stellar masses by up to a factor 40, and larger than the total baryonic mass (gas + stars) by a factor of $\sim 3$. Qualitatively, this result is robust even if the observed velocity gradients reflect ongoing mergers rather than rotating discs. Unless the observed emission line kinematics is dominated by outflows, this implies that the centres of these galaxies are dark-matter dominated or that star formation is $3\times$ less efficient, leading to higher inferred gas masses., Comment: Accepted for publication in A&A. Software for JWST/NIRSpec MSA modelling (slit losses, 1D LSFs and 2D model fitting) publicly available at https://github.com/annadeg/jwst-msafit

There is compelling evidence that the most massive galaxies in the Universe stopped forming stars due to the time-integrated feedback from their central super-massive black holes (SMBHs). However, the exact quenching mechanism is not yet understood, because local massive galaxies were quenched billions of years ago. We present JWST/NIRSpec integral-field spectroscopy observations of GS-10578, a massive, quiescent galaxy at redshift z=3.064. From the spectrum we infer that the galaxy has a stellar mass of $M_*=1.6\pm0.2 \times 10^{11}$ MSun and a dynamical mass $M_{\rm dyn}=2.0\pm0.5 \times 10^{11}$ MSun. Half of its stellar mass formed at z=3.7-4.6, and the system is now quiescent, with the current star-formation rate SFR<9 MSun/yr. We detect ionised- and neutral-gas outflows traced by [OIII] emission and NaI absorption. Outflow velocities reach $v_{\rm out}\approx$1,000 km/s, comparable to the galaxy escape velocity and too high to be explained by star formation alone. GS-10578 hosts an Active Galactic Nucleus (AGN), evidence that these outflows are due to SMBH feedback. The outflow rates are 0.14-2.9 and 30-300 MSun/yr for the ionised and neutral phases, respectively. The neutral outflow rate is ten times higher than the SFR, hence this is direct evidence for ejective SMBH feedback, with mass-loading capable of interrupting star formation by rapidly removing its fuel. Stellar kinematics show ordered rotation, with spin parameter $\lambda_{Re}=0.62\pm0.07$, meaning GS-10578 is rotation supported. This study shows direct evidence for ejective AGN feedback in a massive, recently quenched galaxy, thus clarifying how SMBHs quench their hosts. Quenching can occur without destroying the stellar disc., Comment: 40 pages, 9 figures, submitted to Nat. Ast., comments welcome!

A wide-field near-infrared survey of the Galactic disk and bulge/bar(s) is supported by a large representation of the community of Galactic astronomers. The combination of sensitivity, angular resolution and large field of view make Roman uniquely able to study the crowded and highly extincted lines of sight in the Galactic plane. A ~1000 deg2 survey of the bulge and inner Galactic disk would yield an impressive dataset of ~120 billion sources and map the structure of our Galaxy. The effort would foster subsequent expansions in numerous dimensions (spatial, depth, wavelengths, epochs). Importantly, the survey would benefit from early defintion by the community, namely because the Galactic disk is a complex environment, and different science goals will require trade offs., Comment: Submitted to the Roman Project on October 22 2021 in response to a call for white papers on early-definition Astrophysics opportunity

We study the angular clustering of Quaia, a Gaia- and unWISE-based catalog of over a million quasars with an exceptionally well-defined selection function. With it, we derive cosmology constraints from the amplitude and growth of structure across cosmic time. We divide the sample into two redshift bins, centered at $z=1.0$ and $z=2.1$, and measure both overdensity auto-correlations and cross-correlations with maps of the Cosmic Microwave Background convergence measured by Planck. From these data, and including a prior from measurements of the baryon acoustic oscillations scale, we place constraints on the amplitude of the matter power spectrum $\sigma_8=0.766\pm 0.034$, and on the matter density parameter $\Omega_m=0.343^{+0.017}_{-0.019}$. These measurements are in reasonable agreement with \planck at the $\sim$ 1.4$\sigma$ level, and are found to be robust with respect to observational and theoretical uncertainties. We find that our slightly lower value of $\sigma_8$ is driven by the higher-redshift sample, which favours a low amplitude of matter fluctuations. We present plausible arguments showing that this could be driven by contamination of the CMB lensing map by high-redshift extragalactic foregrounds, which should also affect other cross-correlations with tracers of large-scale structure beyond $z\sim1.5$. Our constraints are competitive with those from state-of-the-art 3$\times$2-point analyses, but arise from a range of scales and redshifts that is highly complementary to those covered by cosmic shear data and most galaxy clustering samples. This, coupled with the unprecedented combination of volume and redshift precision achieved by Quaia allows us to break the usual degeneracy between $\Omega_m$ and $\sigma_8$., Comment: Submitted to JCAP, corrects links to public data release. The per-bin selection functions used in this analysis, and the complete set of power spectra, covariance, and associated metadata used in our fiducial analysis are publicly available at https://doi.org/10.5281/zenodo.8098635 . The public Quaia catalog can be found at https://doi.org/10.5281/zenodo.8060755

We present a new, all-sky quasar catalog, Quaia, that samples the largest comoving volume of any existing spectroscopic quasar sample. The catalog draws on the 6,649,162 quasar candidates identified by the Gaia mission that have redshift estimates from the space observatory's low-resolution BP/RP spectra. This initial sample is highly homogeneous and complete, but has low purity, and 18% of even the bright ($G<20.0$) confirmed quasars have discrepant redshift estimates ($|\Delta z/(1+z)|>0.2$) compared to those from the Sloan Digital Sky Survey (SDSS). In this work, we combine the Gaia candidates with unWISE infrared data (based on the Wide-field Infrared Survey Explorer survey) to construct a catalog useful for cosmological and astrophysical quasar studies. We apply cuts based on proper motions and Gaia and unWISE colors, reducing the number of contaminants by $\sim$4$\times$. We improve the redshifts by training a $k$-nearest neighbors model on SDSS redshifts, and achieve estimates on the $G<20.0$ sample with only 6% (10%) catastrophic errors with $|\Delta z/(1+z)|>0.2$ ($0.1$), a reduction of $\sim$3$\times$ ($\sim$2$\times$) compared to the Gaia redshifts. The final catalog has 1,295,502 quasars with $G<20.5$, and 755,850 candidates in an even cleaner $G<20.0$ sample, with accompanying rigorous selection function models. We compare Quaia to existing quasar catalogs, showing that its large effective volume makes it a highly competitive sample for cosmological large-scale structure analyses. The catalog is publicly available at https://zenodo.org/records/10403370., Comment: Published in the Astrophysical Journal; updated to match published version. Catalog available at https://doi.org/10.5281/zenodo.10403370. Animation available at https://cosmo.nyu.edu/ksf/quaia.mp4

We investigate the incidence and properties of ionized gas outflows in a sample of 52 galaxies with stellar mass between $10^7$ M$_{\odot}$ and $10^9$ M$_{\odot}$ observed with ultra-deep JWST/NIRSpec MSA spectroscopy as part of the JWST Advanced Deep Extragalactic Survey (JADES). The high-spectral resolution (R2700) NIRSpec observations allowed us to identify for the first time the signature of outflows in the rest-frame optical nebular lines in low-mass galaxies at $z>3$. The incidence fraction of ionized outflows, traced by broad components, is about 25-40$\%$ depending on the intensity of the emission lines. The low incidence fraction might be due to both the sensitivity limit and the fact that outflows are not isotropic but have a limited opening angle which results in a detection only when this is directed toward our line of sight. Evidence for outflows increases slightly with stellar mass and star-formation rate. The median velocity and mass loading factor (i.e., the ratio between mass outflow rate and star formation rate) of the outflowing ionized gas are 350 km s$^{-1}$ and $\eta=2.0^{+1.6}_{-1.5}$, respectively. These are 1.5 and 100 times higher, respectively than the typical values observed in local dwarf galaxies. These outflows are able to escape the gravitational potential of the galaxy and enrich the circum-galactic medium and, potentially, the inter-galactic medium. Our results indicate that outflows can significantly impact the star formation activity in low-mass galaxies within the first 2 Gyr of the Universe., Comment: 16 pages, 3 tables

The Nancy Grace Roman Space Telescope is capable of delivering an unprecedented all-sky, high-spatial resolution, multi-epoch infrared map to the astronomical community. This opportunity arises in the midst of numerous ground- and space-based surveys that will provide extensive spectroscopy and imaging together covering the entire sky (such as Rubin/LSST, Euclid, UNIONS, SPHEREx, DESI, SDSS-V, GALAH, 4MOST, WEAVE, MOONS, PFS, UVEX, NEO Surveyor, etc.). Roman can uniquely provide uniform high-spatial-resolution (~0.1 arcsec) imaging over the entire sky, vastly expanding the science reach and precision of all of these near-term and future surveys. This imaging will not only enhance other surveys, but also facilitate completely new science. By imaging the full sky over two epochs, Roman can measure the proper motions for stars across the entire Milky Way, probing 100 times fainter than Gaia out to the very edge of the Galaxy. Here, we propose NANCY: a completely public, all-sky survey that will create a high-value legacy dataset benefiting innumerable ongoing and forthcoming studies of the universe. NANCY is a pure expression of Roman's potential: it images the entire sky, at high spatial resolution, in a broad infrared bandpass that collects as many photons as possible. The majority of all ongoing astronomical surveys would benefit from incorporating observations of NANCY into their analyses, whether these surveys focus on nearby stars, the Milky Way, near-field cosmology, or the broader universe., Comment: Submitted to the call for white papers for the Roman Core Community Survey (June 16th, 2023), and to the Bulletin of the AAS

Reionization is thought to be driven by faint star-forming galaxies, but characterizing this population has long remained very challenging. Here we utilize deep nine-band NIRCam imaging from JADES to study the star-forming and ionizing properties of 756 $z\sim6-9$ galaxies, including hundreds of very UV-faint objects ($M_\mathrm{UV}>-18$). The faintest ($m\sim30$) galaxies in our sample typically have stellar masses of $M_\ast\sim(1-3)\times10^7$ $M_\odot$ and young light-weighted ages ($\sim$50 Myr), though some show strong Balmer breaks implying much older ages ($\sim$500 Myr). We find no evidence for extremely massive galaxies ($>3\times10^{10}$ $M_\odot$). We infer a strong (factor $>$2) decline in the typical [OIII]$+$H$\beta$ EWs towards very faint $z\sim6-9$ galaxies, yet a weak UV luminosity dependence on the H$\alpha$ EWs at $z\sim6$. We demonstrate that these EW trends can be explained if fainter galaxies have systematically lower metallicities as well as more recently-declining star formation histories relative to the most UV-luminous galaxies in our sample. Our data provide evidence that the brightest galaxies are frequently experiencing a recent strong upturn in SFR. We also discuss how the EW trends may be influenced by a strong correlation between $M_\mathrm{UV}$ and Lyman continuum escape fraction. This alternative explanation has dramatically different implications for the contribution of galaxies along the luminosity function to cosmic reionization. Finally, we quantify the photometric overdensities around two $z>7$ strong Ly$\alpha$ emitters. One Ly$\alpha$ emitter lies close to a strong photometric overdensity while the other shows no significant nearby overdensity, perhaps implying that not all strong $z>7$ Ly$\alpha$ emitters reside in large ionized bubbles., Comment: Accepted for publication in MNRAS. Minor changes relative to previous version. Catalog of galaxy coordinates, inferred properties, and photometry available here: https://www.ryan-endsley.com/z6to9-jades-galaxy-catalog

We present a study of the environments of 17 Lyman-$\alpha$ (Ly$\alpha$) emitting galaxies (LAEs) in the reionisation era ($5.8 < z < 8$) identified by JWST/NIRSpec as part of the JWST Advanced Deep Extragalactic Survey (JADES). Unless situated in sufficiently (re)ionised regions, Ly$\alpha$ emission from these galaxies would be strongly absorbed by neutral gas in the intergalactic medium (IGM). We conservatively estimate sizes of the ionised regions required to reconcile the relatively low Ly$\alpha$ velocity offsets ($\Delta v_\text{Ly$\alpha$}<300\,\mathrm{km\,s^{-1}}$) with moderately high Ly$\alpha$ escape fractions ($f_\mathrm{esc,\,Ly\alpha}>5\%$) observed in our sample of LAEs, suggesting the presence of ionised hydrogen along the line of sight towards at least eight out of 17 LAEs. We find minimum physical `bubble' sizes of the order of $R_\text{ion}\sim0.1$-$1\,\mathrm{pMpc}$ are required in a patchy reionisation scenario where ionised bubbles containing the LAEs are embedded in a fully neutral IGM. Around half of the LAEs in our sample are found to coincide with large-scale galaxy overdensities seen in FRESCO at $z \sim 5.8$-$5.9$ and $z\sim7.3$, suggesting Ly$\alpha$ transmission is strongly enhanced in such overdense regions, and underlining the importance of LAEs as tracers of the first large-scale ionised bubbles. Considering only spectroscopically confirmed galaxies, we find our sample of UV-faint LAEs ($M_\text{UV}\gtrsim-20\,\mathrm{mag}$) and their direct neighbours are generally not able to produce the required ionised regions based on the Ly$\alpha$ transmission properties, suggesting lower-luminosity sources likely play an important role in carving out these bubbles. These observations demonstrate the combined power of JWST multi-object and slitless spectroscopy in acquiring a unique view of the early Universe during cosmic reionisation via the most distant LAEs., Comment: 15 pages, 8 figures. Accepted for publication in A&A

We present the properties of 17 faint Ly$\alpha$ emitting galaxies (LAEs) at $z>5.8$ from the JWST Advanced Deep Extragalactic Survey (JADES) in the Hubble Ultra Deep Field/GOODS-S. These LAEs span a redshift range $z\approx5.8-8.0$ and a UV magnitude range $M_{UV}\approx-17$ to $-20.6$, with the Ly$\alpha$ equivalent width (EW) in the range $\approx 25-350$ \AA. The detection of other rest-optical emission lines in the spectra of these LAEs enables the determination of accurate systemic redshifts and Ly{\alpha} velocity offsets, as well as the physical and chemical composition of their stars and interstellar media. These faint LAEs are consistent with metal-poor systems with high ionization parameters, similar to the general galaxy population at $z>6$. We measured an average ionizing photon production efficiency, log($\xi_\rm{ion}$/erg$^{-1}$ Hz) $\approx25.57$ across our LAEs, which does not evolve strongly with redshift. We report an anti-correlation between the Ly$\alpha$ escape fraction (f_\rm{esc}) and the velocity offset from systemic redshift, consistent with model expectations. We further find that the strength and velocity offset of Ly$\alpha$ are neither correlated with galaxy spectroscopic properties nor with $\xi_\rm{ion}$. We find a decrease in $f_\rm{esc}$(Ly$\alpha$) with redshift, indicative of decreasing sizes of ionized bubbles around LAEs at high redshifts. We used a range of galaxy properties to predict Lyman continuum $f_\rm{esc}$ for our LAEs, finding that the ionizing photon output into the intergalactic medium remains roughly constant across the observed Ly$\alpha$ EW, showing a mild increase at fainter M$_{UV}$ and at higher redshifts. We derived correlations between the ionizing photon output from LAEs and $M_{UV}$, Ly$\alpha$ EW and redshifts, which can be used to constrain the ionizing photon contribution of LAEs at $z > 6$ towards cosmic reionization., Comment: Accepted for publication in A&A. 24 pages, 13 figures, spectra presented in the appendix. This updated version includes the addition of one more LAE to the sample, ID 004296 at z=6.712

We present Balmer decrement H$\alpha$/ H$\beta$ measurements for a sample of 51 galaxies at redshifts z = 4-7 observed with the JWST/NIRSpec MSA, as part of the JADES survey. Leveraging 28-hour long exposures and the efficiency of the prism/clear configuration (but also using information from the medium-resolution gratings), we are able to probe directly the low-mass end of the galaxy population, reaching stellar masses Mstar as low as 10^7 Msun . We find that the correlation between Balmer decrement and Mstar is already established at these high redshifts, indicating a rapid build up of dust in moderately massive galaxies at such early epochs. The lowest-mass galaxies in our sample (Mstar = 1-3 x 10^7 Msun ) display a remarkably low Balmer decrement of 2.88 $\pm$ 0.08, consistent with case B, suggesting very little dust content. However, we warn that such a low observed Balmer decrement may also partly be a consequence of an intrinsically lower H$\alpha$/ H$\beta$, resulting from the extreme conditions of the ionized gas in these primeval and unevolved systems. We further compare the Balmer decrement to continuum-derived star-formation rates (SFR), finding tentative evidence of a correlation, which likely traces the underlying connection between SFR and mass of cold gas. However, we note that larger samples are required to distinguish between direct and primary correlations from indirect and secondary dependencies at such high redshifts.