-
GD-1 Stellar Stream and Cocoon in the DESI Early Data Release
Authors:
Monica Valluri,
Parker Fagrelius,
Sergey. E. Koposov,
Ting S. Li,
Oleg Y. Gnedin,
Eric F. Bell,
Raymond G. Carlberg,
Andrew P. Cooper,
Jessia N. Aguilar,
Carlos Allende Prieto,
Vasily Belokurov,
Leandro Beraldo e Silva,
David Brooks,
Amanda Byström,
Todd Claybaugh,
Kyle Dawson,
Arjun Dey,
Peter Doel,
Jaime E. Forero-Romero,
Enrique Gaztañaga,
Satya Gontcho A Gontcho,
Klaus Honscheid,
T . Kisner,
Anthony Kremin,
A. Lambert
, et al. (27 additional authors not shown)
Abstract:
We present ~ 126 new spectroscopically identified members of the GD-1 tidal stream obtained with the 5000-fiber Dark Energy Spectroscopic Instrument (DESI). We confirm the existence of a ``cocoon'' which is broad (FWHM~2.932deg~460pc) and kinematically hot (velocity dispersion, sigma~5-8km/s) component that surrounds a narrower (FWHM~0.353deg~55pc) and colder (sigma~ 2.2-2.6km/s) thin stream compo…
▽ More
We present ~ 126 new spectroscopically identified members of the GD-1 tidal stream obtained with the 5000-fiber Dark Energy Spectroscopic Instrument (DESI). We confirm the existence of a ``cocoon'' which is broad (FWHM~2.932deg~460pc) and kinematically hot (velocity dispersion, sigma~5-8km/s) component that surrounds a narrower (FWHM~0.353deg~55pc) and colder (sigma~ 2.2-2.6km/s) thin stream component (based on a median per star velocity precision of 2.7km/s). The cocoon extends over at least a ~ 20deg segment of the stream observed by DESI. The thin and cocoon components have similar mean values of [Fe/H]: -2.54+/- 0.04dex and -2.45+/-0.06dex suggestive of a common origin. The data are consistent with the following scenarios for the origin of the cocoon. The progenitor of the GD-1 stream was an accreted globular cluster (GC) and: (a) the cocoon was produced by pre-accretion tidal stripping of the GC while it was still inside its parent dwarf galaxy; (b) the cocoon is debris from the parent dwarf galaxy; (c) an initially thin GC tidal stream was heated by impacts from dark subhalos in the Milky Way; (d) an initially thin GC stream was heated by a massive Sagittarius dwarf galaxy; or a combination of some these. In the first two cases the velocity dispersion and mean metallicity are consistent with the parent dwarf galaxy having a halo mass of ~0^9\msun. Future DESI spectroscopy and detailed modeling may enable us to distinguish between these possible origins.
△ Less
Submitted 8 July, 2024;
originally announced July 2024.
-
Kinematics and Dynamics of the Galactic Bar revealed by Gaia Long Period Variables
Authors:
Han-Yuan Zhang,
Vasily Belokurov,
N. Wyn Evans,
Sarah G. Kane,
Jason L. Sanders
Abstract:
We take low-amplitude, long period variable (LA-LPV) candidates in Gaia DR3 as tracers of the kinematics and dynamics of the Milky Way bar. LA-LPVs, like other LPVs, have high luminosities and follow a tight period-luminosity relation, but unlike e.g. Mira variables, their radial velocity measurements are reliable due to their smaller pulsation amplitudes. We supplement the Gaia astrometric and ra…
▽ More
We take low-amplitude, long period variable (LA-LPV) candidates in Gaia DR3 as tracers of the kinematics and dynamics of the Milky Way bar. LA-LPVs, like other LPVs, have high luminosities and follow a tight period-luminosity relation, but unlike e.g. Mira variables, their radial velocity measurements are reliable due to their smaller pulsation amplitudes. We supplement the Gaia astrometric and radial velocity measurements with distance moduli assigned using a period-luminosity relation to acquire full 6D phase space information. The assigned distances are validated by comparing to geometric distances and StarHorse distances, which shows biases less than $\sim5\%$. Our sample provides an unprecedented panoramic picture of the inner Galaxy with minimal selection function effects. We map the kinematics of the inner Milky Way and find a significant kinematic signature corresponding to the Galactic bar. We measure the pattern speed of the Galactic bar using the continuity equation and find $Ω_{\rm b}=34.1\pm2.4$ km s$^{-1}$ kpc$^{-1}$. We develop a simple, robust and model-independent method to measure the dynamical length of the bar using only kinematics and find $R_{\rm b}\sim4.0$ kpc. We validate both measurements using N-body simulations. Assuming knowledge of the gravitational potential of the inner Milky Way, we analyse the orbital structure of the Galactic bar using orbital frequency ratios. The $x_1$ orbits are the dominant bar-supporting orbital family in our sample. Amongst the selected bar stars, the $x_1 v_1$ or "banana" orbits constitute a larger fraction ($\sim 15\%$) than other orbital families in the bar, implying that they are the dominant family contributing to the Galactic X-shape, although contributions from other orbital families are present.
△ Less
Submitted 10 June, 2024;
originally announced June 2024.
-
Trojan globular clusters: radial migration via trapping in bar resonances
Authors:
Adam M. Dillamore,
Stephanie Monty,
Vasily Belokurov,
N. Wyn Evans
Abstract:
We search for globular clusters (GCs) trapped in resonances with the bar of the Milky Way. By integrating their orbits in a potential with a decelerating bar, we select 10 whose orbits are significantly changed by its presence. Most of these are trapped in the corotation resonance (CR), including M22 and 47 Tuc. The decelerating bar is capable of transporting these GCs to their current positions f…
▽ More
We search for globular clusters (GCs) trapped in resonances with the bar of the Milky Way. By integrating their orbits in a potential with a decelerating bar, we select 10 whose orbits are significantly changed by its presence. Most of these are trapped in the corotation resonance (CR), including M22 and 47 Tuc. The decelerating bar is capable of transporting these GCs to their current positions from much lower energies, angular momenta, and radii. Our results indicate that the bar is likely to have reshaped the Milky Way's globular cluster system via its resonances. We also discuss implications for the origins of specific GCs, including the possible nuclear star cluster M22. Finally, we consider the effects of the bar on the tidal tails of a trapped GC, by running simulations of stars stripped from 47 Tuc. Instead of forming narrow tails, the stripped stars make up a diffuse extended halo around the cluster, consistent with observations of 47 Tuc.
△ Less
Submitted 23 May, 2024;
originally announced May 2024.
-
Euclid: Early Release Observations -- Deep anatomy of nearby galaxies
Authors:
L. K. Hunt,
F. Annibali,
J. -C. Cuillandre,
A. M. N. Ferguson,
P. Jablonka,
S. S. Larsen,
F. R. Marleau,
E. Schinnerer,
M. Schirmer,
C. Stone,
C. Tortora,
T. Saifollahi,
A. Lançon,
M. Bolzonella,
S. Gwyn,
M. Kluge,
R. Laureijs,
D. Carollo,
M. L. M. Collins,
P. Dimauro,
P. -A. Duc,
D. Erkal,
J. M. Howell,
C. Nally,
E. Saremi
, et al. (174 additional authors not shown)
Abstract:
Euclid is poised to make significant advances in the study of nearby galaxies in the local Universe. Here we present a first look at 6 galaxies observed for the Nearby Galaxy Showcase as part of the Euclid Early Release Observations acquired between August and November, 2023. These targets, 3 dwarf galaxies (HolmbergII, IC10, NGC6822) and 3 spirals (IC342, NGC2403, NGC6744), range in distance from…
▽ More
Euclid is poised to make significant advances in the study of nearby galaxies in the local Universe. Here we present a first look at 6 galaxies observed for the Nearby Galaxy Showcase as part of the Euclid Early Release Observations acquired between August and November, 2023. These targets, 3 dwarf galaxies (HolmbergII, IC10, NGC6822) and 3 spirals (IC342, NGC2403, NGC6744), range in distance from about 0.5 Mpc to 8.8 Mpc. Our assessment of the surface brightness depths in the stacked Euclid images confirms previous estimates in 100 arcsec^2 regions of 1sigma=30.5 mag/arcsec^2 for VIS, but slightly deeper than previous estimates for NISP with 1sigma=29.2-29.4 mag/arcsec^2. By combining Euclid HE, YE, and IE into RGB images, we illustrate the large field-of-view covered by a single Reference Observing Sequence, together with exquisite detail on parsec scales in these nearby galaxies. Radial surface brightness and color profiles demonstrate galaxy colors in agreement with stellar population synthesis models. Standard stellar photometry selection techniques find approximately 1.3 million stars across the 6 galaxy fields. Euclid's resolved stellar photometry allows us to constrain the star-formation histories of these galaxies, by disentangling the distributions of young stars, as well as asymptotic giant branch and red giant branch stellar populations. We finally examine 2 galaxies individually for surrounding satellite systems. Our analysis of the ensemble of dwarf satellites around NGC6744 reveals a new galaxy, EDwC1, a nucleated dwarf spheroidal at the end of a spiral arm. Our new census of the globular clusters around NGC2403 yields 9 new star-cluster candidates, 8 of which with colors indicative of evolved stellar populations. In summary, our investigation of the 6 Showcase galaxies demonstrates that Euclid is a powerful probe of the anatomy of nearby galaxies [abridged].
△ Less
Submitted 22 May, 2024;
originally announced May 2024.
-
Euclid: Early Release Observations -- Programme overview and pipeline for compact- and diffuse-emission photometry
Authors:
J. -C. Cuillandre,
E. Bertin,
M. Bolzonella,
H. Bouy,
S. Gwyn,
S. Isani,
M. Kluge,
O. Lai,
A. Lançon,
D. A. Lang,
R. Laureijs,
T. Saifollahi,
M. Schirmer,
C. Stone,
Abdurro'uf,
N. Aghanim,
B. Altieri,
F. Annibali,
H. Atek,
P. Awad,
M. Baes,
E. Bañados,
D. Barrado,
S. Belladitta,
V. Belokurov
, et al. (240 additional authors not shown)
Abstract:
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 t…
▽ More
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.
△ Less
Submitted 22 May, 2024;
originally announced May 2024.
-
Close Binary Fractions in accreted and in-situ Halo Stars
Authors:
Dolev Bashi,
Vasily Belokurov,
Simon Hodgkin
Abstract:
The study of binary stars in the Galactic halo provides crucial insights into the dynamical history and formation processes of the Milky Way. In this work, we aim to investigate the binary fraction in a sample of accreted and in-situ halo stars, focusing on short-period binaries. Utilising data from Gaia DR3, we analysed the radial velocity (RV) uncertainty $σ_{\mathrm{RV}}$ distribution of a samp…
▽ More
The study of binary stars in the Galactic halo provides crucial insights into the dynamical history and formation processes of the Milky Way. In this work, we aim to investigate the binary fraction in a sample of accreted and in-situ halo stars, focusing on short-period binaries. Utilising data from Gaia DR3, we analysed the radial velocity (RV) uncertainty $σ_{\mathrm{RV}}$ distribution of a sample of main-sequence stars. We used a novel Bayesian framework to model the dependence in $σ_{\mathrm{RV}}$ of single and binary systems allowing us to estimate binary fractions $F$ in a sample of bright ($G_{\mathrm{RVS}}$ < 12) Gaia sources. We selected the samples of in-situ and accreted halo stars based on estimating the 6D phase space information and affiliating the stars to the different samples on an action-angle vs energy ($L_{\mathrm{z}}-E$) diagram. Our results indicate a higher, though not significant, binary fraction in accreted stars compared to the in-situ halo sample. We further explore binary fractions using cuts in $E$ and $L_z$ and find a higher binary fraction in both high-energy and prograde orbits that might be explained by differences in metallicity. By cross-matching our Gaia sample with APOGEE DR17 catalogue, we confirm the results of previous studies on higher binary fractions in metal-poor stars and find the fractions of accreted and in-situ halo stars consistent with this trend. Our finding provides new insights into binary stars' formation processes and dynamical evolution in the primordial Milky Way Galaxy and its accreted dwarf Galaxies.
△ Less
Submitted 20 May, 2024;
originally announced May 2024.
-
The ratio of [Eu/$α$] differentiates accreted/in-situ Milky Way stars across metallicities, as indicated by both field stars and globular clusters
Authors:
Stephanie Monty,
Vasily Belokurov,
Jason L. Sanders,
Terese T. Hansen,
Charli M. Sakari,
Madeleine McKenzie,
GyuChul Myeong,
Ellot Y. Davies,
Anke Ardern-Arentsen,
Davide Massari
Abstract:
We combine stellar orbits with the abundances of the heavy, $r$-process element europium and the light, $α$-element, silicon to separate in-situ and accreted populations in the Milky Way across all metallicities. At high orbital energy, the accretion-dominated halo shows elevated values of [Eu/Si], while at lower energies, where many of the stars were born in-situ, the levels of [Eu/Si] are lower.…
▽ More
We combine stellar orbits with the abundances of the heavy, $r$-process element europium and the light, $α$-element, silicon to separate in-situ and accreted populations in the Milky Way across all metallicities. At high orbital energy, the accretion-dominated halo shows elevated values of [Eu/Si], while at lower energies, where many of the stars were born in-situ, the levels of [Eu/Si] are lower. These systematically different levels of [Eu/Si] in the MW and the accreted halo imply that the scatter in [Eu/$α$] within a single galaxy is smaller than previously thought. At the lowest metallicities, we find that both accreted and in-situ populations trend down in [Eu/Si], consistent with enrichment via neutron star mergers. Through compiling a large dataset of abundances for 46 globular clusters (GCs), we show that differences in [Eu/Si] extend to populations of in-situ/accreted GCs. We interpret this consistency as evidence that in $r$-process elements, GCs trace the star formation history of their hosts, motivating their use as sub-Gyr timers of galactic evolution. Furthermore, fitting the trends in [Eu/Si] using a simple galactic chemical evolution model, we find that differences in [Eu/Si] between accreted and in-situ MW field stars cannot be explained through star formation efficiency alone. Finally, we show that the use of [Eu/Si] as a chemical tag between GCs and their host galaxies extends beyond the Local Group, to the halo of M31 - potentially offering the opportunity to do Galactic Archaeology in an external galaxy.
△ Less
Submitted 24 May, 2024; v1 submitted 14 May, 2024;
originally announced May 2024.
-
Stochastic star formation and the abundance of $z>10$ UV-bright galaxies
Authors:
Andrey Kravtsov,
Vasily Belokurov
Abstract:
We use a well-motivated galaxy formation framework to predict stellar masses, star formation rates (SFR), and ultraviolet (UV) luminosities of galaxy populations at redshifts $z\in 5-16$, taking into account stochasticity of SFR in a controlled manner. We demonstrate that the model can match observational estimates of UV luminosity functions (LFs) at $5<z<10$ with a modest level of SFR stochastici…
▽ More
We use a well-motivated galaxy formation framework to predict stellar masses, star formation rates (SFR), and ultraviolet (UV) luminosities of galaxy populations at redshifts $z\in 5-16$, taking into account stochasticity of SFR in a controlled manner. We demonstrate that the model can match observational estimates of UV luminosity functions (LFs) at $5<z<10$ with a modest level of SFR stochasticity, resulting in the scatter of absolute UV luminosity at a given halo mass of $σ_{M_{\rm UV}}\approx 0.75$. To match the observed UV LFs at $z\approx 11-13$ and $z\approx 16$ the SFR stochasticity should increase so that $σ_{M_{\rm UV}}\approx 1-1.3$ and $\approx 2$, respectively. Model galaxies at $z\approx 11-13$ have stellar masses and SFRs in good agreement with existing measurements. The median fraction of the baryon budget that was converted into stars, $f_\star$, is only $f_\star\approx 0.005-0.05$, but a small fraction of galaxies at $z=16$ have $f_\star>1$ indicating that SFR stochasticity cannot be higher. We discuss several testable consequences of the increased SFR stochasticity at $z>10$. The increase of SFR stochasticity with increasing $z$, for example, prevents steepening of UV LF and even results in some flattening of UV LF at $z\gtrsim 13$. The median stellar ages of model galaxies at $z\approx 11-16$ are predicted to decrease from $\approx 20-30$ Myr for $M_{\rm UV}\gtrsim -21$ galaxies to $\approx 5-10$ Myr for brighter ones. Likewise, the scatter in median stellar age is predicted to decrease with increasing luminosity. The scatter in the ratio of star formation rates averaged over 10 and 100 Myr should increase with redshift. Fluctuations of ionizing flux should increase at $z>10$ resulting in the increasing scatter in the line fluxes and their ratios for the lines sensitive to ionization parameter.
△ Less
Submitted 7 May, 2024;
originally announced May 2024.
-
Gaia DR3 detectability of unresolved binary systems
Authors:
Alfred Castro-Ginard,
Zephyr Penoyre,
Andrew R. Casey,
Anthony G. A. Brown,
Vasily Belokurov,
Tristan Cantat-Gaudin,
Ronald Drimmel,
Morgan Fouesneau,
Shourya Khanna,
Evgeny P. Kurbatov,
Adrian M. Price-Whelan,
Hans-Walter Rix,
Richard L. Smart
Abstract:
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 po…
▽ More
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
△ Less
Submitted 22 April, 2024;
originally announced April 2024.
-
Radial halo substructure in harmony with the Galactic bar
Authors:
Adam M. Dillamore,
Vasily Belokurov,
N. Wyn Evans
Abstract:
Overdensities in the radial phase space $(r,v_r)$ of the Milky Way's halo have previously been associated with the phase-mixed debris of a highly radial merger event, such as Gaia Sausage-Enceladus. We present and test an alternative theory in which the overdense 'chevrons' are instead composed of stars trapped in resonances with the Galactic bar. We develop an analytic model of resonant orbits in…
▽ More
Overdensities in the radial phase space $(r,v_r)$ of the Milky Way's halo have previously been associated with the phase-mixed debris of a highly radial merger event, such as Gaia Sausage-Enceladus. We present and test an alternative theory in which the overdense 'chevrons' are instead composed of stars trapped in resonances with the Galactic bar. We develop an analytic model of resonant orbits in the isochrone potential, and complement this with a test particle simulation of a stellar halo in a realistic barred Milky Way potential. These models are used to predict the appearance of action space $(J_φ,J_r)$ and radial phase space in the Solar neighbourhood. They are able to reproduce almost all salient features of the observed chevrons. In particular, both the analytic model and simulation predict that the chevrons are more prominent at $v_r<0$ when viewed near the Sun, as is observed by Gaia. This is inconsistent with formation by an ancient merger event. We also associate individual chevrons with specific resonances. At a bar pattern speed of $Ω_\mathrm{b}=35$ km s$^{-1}$kpc$^{-1}$, the two most prominent prograde chevrons align very closely with the corotation and outer Lindblad resonances. The former can be viewed as a highly eccentric extension of the Hercules stream. Finally, our model predicts that the $v_r$ asymmetry changes sign as a function of Galactic radius and azimuth, and we find evidence that this is indeed the case in the Milky Way.
△ Less
Submitted 22 February, 2024;
originally announced February 2024.
-
Galactic Archaeology with Gaia
Authors:
Alis J. Deason,
Vasily Belokurov
Abstract:
The Gaia mission has revolutionized our view of the Milky Way and its satellite citizens. The field of Galactic Archaeology has been piecing together the formation and evolution of the Galaxy for decades, and we have made great strides, with often limited data, towards discovering and characterizing the subcomponents of the Galaxy and its building blocks. Now, the exquisite 6D phase-space plus che…
▽ More
The Gaia mission has revolutionized our view of the Milky Way and its satellite citizens. The field of Galactic Archaeology has been piecing together the formation and evolution of the Galaxy for decades, and we have made great strides, with often limited data, towards discovering and characterizing the subcomponents of the Galaxy and its building blocks. Now, the exquisite 6D phase-space plus chemical information from Gaia and its complementary spectroscopic surveys has handed us a plethora of data to pour over as we move towards a quantitative rather than qualitative view of the Galaxy and its progenitors. We review the state of the field in the post-Gaia era, and examine the key lessons that will dictate the future direction of Galactic halo research.
△ Less
Submitted 1 July, 2024; v1 submitted 19 February, 2024;
originally announced February 2024.
-
Uniting Gaia and APOGEE to unveil the cosmic chemistry of the Milky Way disc
Authors:
Tristan Cantat-Gaudin,
Morgan Fouesneau,
Hans-Walter Rix,
Anthony G. A. Brown,
Ronald Drimmel,
Alfred Castro-Ginard,
Shourya Khanna,
Vasily Belokurov,
Andrew R. Casey
Abstract:
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 b…
▽ More
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.
△ Less
Submitted 10 January, 2024;
originally announced January 2024.
-
Multiple Populations and a CH Star Found in the 300S Globular Cluster Stellar Stream
Authors:
Sam A. Usman,
Alexander P. Ji,
Ting S. Li,
Andrew B. Pace,
Lara R. Cullinane,
Gary S. Da Costa,
Sergey E. Koposov,
Geraint F. Lewis,
Daniel B. Zucker,
Vasily Belokurov,
Joss Bland-Hawthorn,
Peter S. Ferguson,
Terese T. Hansen,
Guilherme Limberg,
Sarah L. Martell,
Madeleine McKenzie,
Joshua D. Simon
Abstract:
Milky Way globular clusters (GCs) display chemical enrichment in a phenomenon called multiple stellar populations (MSPs). While the enrichment mechanism is not fully understood, there is a correlation between a cluster's mass and the fraction of enriched stars found therein. However, present-day GC masses are often smaller than their masses at the time of formation due to dynamical mass loss. In t…
▽ More
Milky Way globular clusters (GCs) display chemical enrichment in a phenomenon called multiple stellar populations (MSPs). While the enrichment mechanism is not fully understood, there is a correlation between a cluster's mass and the fraction of enriched stars found therein. However, present-day GC masses are often smaller than their masses at the time of formation due to dynamical mass loss. In this work, we explore the relationship between mass and MSPs using the stellar stream 300S. We present the chemical abundances of eight red giant branch member stars in 300S with high-resolution spectroscopy from Magellan/MIKE. We identify one enriched star characteristic of MSPs and no detectable metallicity dispersion, confirming that the progenitor of 300S was a globular cluster. The fraction of enriched stars (12.5\%) observed in our 300S stars is less than the 50\% of stars found enriched in Milky Way GCs of comparable present-day mass ($\sim10^{4.5}$\msun). We calculate the mass of 300S's progenitor and compare it to the initial masses of intact GCs, finding that 300S aligns well with the trend between the system mass at formation and enrichment. 300S's progenitor may straddle the critical mass threshold for the formation of MSPs and can therefore serve as a benchmark for the stellar enrichment process. Additionally, we identify a CH star, with high abundances of \textit{s}-process elements, probably accreted from a binary companion. The rarity of such binaries in intact GCs may imply stellar streams permit the survival of binaries that would otherwise be disrupted.
△ Less
Submitted 4 January, 2024;
originally announced January 2024.
-
Hints of a disrupted binary dwarf galaxy in the Sagittarius stream
Authors:
Elliot Y. Davies,
Stephanie Monty,
Vasily Belokurov,
Adam M. Dillamore
Abstract:
In this work, we look for evidence of a non-unity mass ratio binary dwarf galaxy merger in the Sagittarius stream. Simulations of such a merger show that, upon merging with a host, particles from the less-massive galaxy will often mostly be found in the extended stream and less-so in the central remnant. Motivated by these simulations, we use APOGEE DR17 chemical data from approximately 1100 stars…
▽ More
In this work, we look for evidence of a non-unity mass ratio binary dwarf galaxy merger in the Sagittarius stream. Simulations of such a merger show that, upon merging with a host, particles from the less-massive galaxy will often mostly be found in the extended stream and less-so in the central remnant. Motivated by these simulations, we use APOGEE DR17 chemical data from approximately 1100 stars in both the Sagittarius remnant and stream to look for evidence of contamination from a second dwarf galaxy. This search is initially justified by the idea that disrupted binary dwarf galaxies provide a possible explanation of the Sagittarius bifurcation, and the location of the massive, chemically peculiar globular cluster NGC 2419 found within the stream of Sagittarius. We separate the Sagittarius data into its remnant and stream and compare the [Mg/Fe] content of the two populations. In particular, we select [Mg/Fe] to search for hints of unique star formation histories among our sample stars. Comparing the stream and remnant populations, we find regions have distinct [Mg/Fe] distributions for fixed [Fe/H], in addition to distinct chemical tracks in [Mg/Fe] -- [Fe/H] abundance space. We show that there are large regions of the tracks for which the probability of the two samples being drawn from the same distribution is very low (p < 0.05). Furthermore, we show that the two tracks can be fit with unique star formation histories using simple, one zone galactic chemical evolution models. While more work must be done to discern whether the hypothesis presented here is true, our work hints at the possibility that Sagittarius may consist of two dwarf galaxy progenitors.
△ Less
Submitted 13 December, 2023;
originally announced December 2023.
-
The Pristine Inner Galaxy Survey (PIGS) VIII: Characterising the orbital properties of the ancient, very metal-poor inner Milky Way
Authors:
Anke Ardern-Arentsen,
Giacomo Monari,
Anna B. A. Queiroz,
Else Starkenburg,
Nicolas F. Martin,
Cristina Chiappini,
David S. Aguado,
Vasily Belokurov,
Ray Carlberg,
Stephanie Monty,
GyuChul Myeong,
Mathias Schultheis,
Federico Sestito,
Kim A. Venn,
Sara Vitali,
Zhen Yuan,
Hanyuan Zhang,
Sven Buder,
Geraint F. Lewis,
William H. Oliver,
Zhen Wan,
Daniel B. Zucker
Abstract:
The oldest stars in the Milky Way (born in the first few billion years) are expected to have a high density in the inner few kpc, spatially overlapping with the Galactic bulge. We use spectroscopic data from the Pristine Inner Galaxy Survey (PIGS) to study the dynamical properties of ancient, metal-poor inner Galaxy stars. We compute distances using StarHorse, and orbital properties in a barred Ga…
▽ More
The oldest stars in the Milky Way (born in the first few billion years) are expected to have a high density in the inner few kpc, spatially overlapping with the Galactic bulge. We use spectroscopic data from the Pristine Inner Galaxy Survey (PIGS) to study the dynamical properties of ancient, metal-poor inner Galaxy stars. We compute distances using StarHorse, and orbital properties in a barred Galactic potential. With this paper, we release the spectroscopic AAT/PIGS catalogue (13 235 stars). We find that most PIGS stars have orbits typical for a pressure-supported population. The fraction of stars confined to the inner Galaxy decreases with decreasing metallicity, but many very metal-poor stars (VMP, [Fe/H] < -2.0) stay confined (~ 60% stay within 5 kpc). The azimuthal velocity v$_φ$ also decreases between [Fe/H] = -1.0 and -2.0, but is constant for VMP stars (at ~ 40 km/s). The carbon-enhanced metal-poor (CEMP) stars in PIGS appear to have similar orbital properties compared to normal VMP stars. Our results suggest a possible transition between two spheroidal components - a more metal-rich, more concentrated, faster rotating component, and a more metal-poor, more extended and slower/non-rotating component. We propose that the former may be connected to pre-disc in-situ stars (or those born in large building blocks), whereas the latter may be dominated by contributions from smaller galaxies. This is an exciting era where large metal-poor samples, such as in this work (as well as upcoming surveys, e.g., 4MOST), shed light on the earliest evolution of our Galaxy.
△ Less
Submitted 16 April, 2024; v1 submitted 6 December, 2023;
originally announced December 2023.
-
On the existence of a very metal-poor disc in the Milky Way
Authors:
Hanyuan Zhang,
Anke Ardern-Arentsen,
Vasily Belokurov
Abstract:
There has been a discussion for many years on whether the disc in the Milky Way extends down to low metallicity. We aim to address the question by employing a large sample of giant stars with radial velocities and homogeneous metallicities based on the Gaia DR3 XP spectra. We study the 3D velocity distribution of stars in various metallicity ranges, including the very-metal poor regime (VMP, [M/H]…
▽ More
There has been a discussion for many years on whether the disc in the Milky Way extends down to low metallicity. We aim to address the question by employing a large sample of giant stars with radial velocities and homogeneous metallicities based on the Gaia DR3 XP spectra. We study the 3D velocity distribution of stars in various metallicity ranges, including the very-metal poor regime (VMP, [M/H] $<-2.0$). We find that a clear disc population starts to emerge only around [M/H] $\sim -1.3$, and is not visible for [M/H] $<-1.6$. Using Gaussian Mixture Modeling (GMM), we show that there are two halo populations in the VMP regime: one stationary and one with a net prograde rotation of $\sim80\,\mathrm{km/s}$. In this low-metallicity range, we are able to place constraints on the contribution of a rotation-supported disc sub-population to a maximum of $\sim 3$\%. We compare our results to previous claims of discy VMP stars in both observations and simulations and find that having a prograde halo component could explain most of these.
△ Less
Submitted 15 November, 2023;
originally announced November 2023.
-
In-situ vs accreted Milky Way globular clusters: a new classification method and implications for cluster formation
Authors:
Vasily Belokurov,
Andrey Kravtsov
Abstract:
We present a new scheme for the classification of the in-situ and accreted globular clusters (GCs). The scheme uses total energy $E$ and $z$-component of the orbital angular momentum and is calibrated using [Al/Fe] abundance ratio. We demonstrate that such classification results in the GC populations with distinct spatial, kinematic, and chemical abundance distributions. The in-situ GCs are distri…
▽ More
We present a new scheme for the classification of the in-situ and accreted globular clusters (GCs). The scheme uses total energy $E$ and $z$-component of the orbital angular momentum and is calibrated using [Al/Fe] abundance ratio. We demonstrate that such classification results in the GC populations with distinct spatial, kinematic, and chemical abundance distributions. The in-situ GCs are distributed within the central 10 kpc of the Galaxy in a flattened configuration aligned with the MW disc, while the accreted GCs have a wide distribution of distances and a spatial distribution close to spherical. In-situ and accreted GCs have different $\rm [Fe/H]$ distributions with the well-known bimodality present only in the metallicity distribution of the in-situ GCs. Furthermore, the accreted and in-situ GCs are well separated in the plane of $\rm [Al/Fe]-[Mg/Fe]$ abundance ratios and follow distinct sequences in the age--$\rm [Fe/H]$ plane. The in-situ GCs in our classification show a clear disc spin-up signature -- the increase of median $V_φ$ at metallicities $\rm [Fe/H]\approx -1.3÷-1$ similar to the spin-up in the in-situ field stars. This signature signals the MW's disc formation, which occurred $\approx 11.7-12.7$ Gyrs ago (or at $z\approx 3.1-5.3$) according to GC ages. In-situ GCs with metallicities of $\rm [Fe/H]\gtrsim -1.3$ were thus born in the Milky Way disc, while lower metallicity in-situ GCs were born during early, turbulent, pre-disc stages of the evolution of the Galaxy and are part of its Aurora stellar component.
△ Less
Submitted 27 September, 2023;
originally announced September 2023.
-
Taking the Milky Way for a spin: disc formation in the ARTEMIS simulations
Authors:
Adam M. Dillamore,
Vasily Belokurov,
Andrey Kravtsov,
Andreea S. Font
Abstract:
We investigate the formation (spin-up) of galactic discs in the ARTEMIS simulations of Milky Way-mass galaxies. In almost all galaxies discs spin up at higher [Fe/H] than the Milky Way (MW). Those that contain an analogue of the Gaia Sausage-Enceladus (GSE) spin up at a lower average metallicity than those without. We identify six galaxies with spin-up metallicity similar to that of the MW, which…
▽ More
We investigate the formation (spin-up) of galactic discs in the ARTEMIS simulations of Milky Way-mass galaxies. In almost all galaxies discs spin up at higher [Fe/H] than the Milky Way (MW). Those that contain an analogue of the Gaia Sausage-Enceladus (GSE) spin up at a lower average metallicity than those without. We identify six galaxies with spin-up metallicity similar to that of the MW, which form their discs $\sim 8-11$ Gyr ago. Five of these experience a merger similar to the GSE. The spin-up times correlate with the halo masses at early times: galaxies with early spin-up have larger virial masses at a lookback time $t_L=12$ Gyr. The fraction of stars accreted from outside the host galaxy is smaller in galaxies with earlier spin-ups. Accreted fractions small enough to be comparable to the MW are only found in galaxies with the earliest disc formation and large initial virial masses ($M_\mathrm{200c} \approx2\times10^{11}M_\odot$ at $t_L=12$ Gyr). We find that discs form when the halo's virial mass reaches a threshold of $M_\mathrm{200c}\approx(6\pm3)\times10^{11}M_\odot$, independent of the spin-up time. However, the failure to form a disc in other galaxies appears to be instead related to mergers at early times. We also find that discs form when the central potential is not particularly steep. Our results indicate that the MW assembled its mass and formed its disc earlier than the average galaxy of a similar mass.
△ Less
Submitted 18 December, 2023; v1 submitted 15 September, 2023;
originally announced September 2023.
-
Disrupted dwarf binary merger as the possible origin of NGC 2419 and Sagittarius stream substructure
Authors:
Elliot Y. Davies,
Vasily Belokurov,
Stephanie Monty,
N. Wyn Evans
Abstract:
Using $N$-body simulations, we demonstrate that satellite dwarf galaxy pairs which undergo significant mixing ($\sim 6$ Gyr) can have their respective most bound particles separated great distances upon subsequently merging with a more massive host. This may provide an explanation as to the origin of the complex globular cluster NGC 2149, which is found within the tail of the Sagittarius dwarf sph…
▽ More
Using $N$-body simulations, we demonstrate that satellite dwarf galaxy pairs which undergo significant mixing ($\sim 6$ Gyr) can have their respective most bound particles separated great distances upon subsequently merging with a more massive host. This may provide an explanation as to the origin of the complex globular cluster NGC 2149, which is found within the tail of the Sagittarius dwarf spheroidal galaxy, yet separated from its central remnant by over 100 kpc. Dynamical investigations could support the chemical evidence which already points to the NGC 2419 being a nuclear star cluster. Motivated by the distinct nature of NGC 2419, we run a suite of simulations whereby an initial pre-infall merger of two satellites is followed by a post-infall merger of the remnant into a MW-like host potential. We present a striking example from our suite in this work, in which this separation is reproduced by the most bound particles of the two pre-infall satellites. Additionally, this double merger scenario can induce unusual on-sky features in the tidal debris of the post-infall merger, such as clouds, overdensities, and potentially new arms that could resemble the bifurcation observed in Sagittarius.
△ Less
Submitted 3 August, 2023;
originally announced August 2023.
-
Stellar halo striations from assumptions of axisymmetry
Authors:
Elliot Y. Davies,
Adam M. Dillamore,
Vasily Belokurov,
N. Wyn Evans
Abstract:
Motivated by the LMC's impact on the integral of motion space of the stellar halo, we run an $N$-body merger simulation to produce a population of halo-like stars. We subsequently move to a test particle simulation, in which the LMC perturbs this debris. When an axisymmetric potential is assumed for the final snapshot of the $N$-body merger remnant, a series of vertical striations in $(L_z, E)$ sp…
▽ More
Motivated by the LMC's impact on the integral of motion space of the stellar halo, we run an $N$-body merger simulation to produce a population of halo-like stars. We subsequently move to a test particle simulation, in which the LMC perturbs this debris. When an axisymmetric potential is assumed for the final snapshot of the $N$-body merger remnant, a series of vertical striations in $(L_z, E)$ space form as the LMC approaches its pericentre. These result from the formation of overdensities in angular momentum owing to a relationship between the precession rate of near radial orbits and the torquing of these orbits by the LMC. This effect is heavily dependent on the shape of the inner potential. If a quadrupole component of the potential is included these striations become significantly less apparent due to the difference in precession rate between the two potentials. The absence of these features in data, and the dramatic change in orbital plane precession rate, discourages the use of an axisymmetric potential for highly eccentric orbits accreted from a massive GSE-like merger. Given the link between appearance of these striations and the shape of the potential, this effect may provide a new method of constraining the axisymmetry of the halo.
△ Less
Submitted 30 June, 2023;
originally announced July 2023.
-
Nitrogen enrichment and clustered star formation at the dawn of the Galaxy
Authors:
Vasily Belokurov,
Andrey Kravtsov
Abstract:
Anomalously high nitrogen-to-oxygen abundance ratios [N/O] are observed in globular clusters (GCs), among the field stars of the Milky Way (MW), and even in the gas in a $z\approx 11$ galaxy. Using data from the APOGEE Data Release 17 and the Gaia Data Release 3, we present several independent lines of evidence that most of the MW's high-[N/O] stars were born in situ in massive bound clusters duri…
▽ More
Anomalously high nitrogen-to-oxygen abundance ratios [N/O] are observed in globular clusters (GCs), among the field stars of the Milky Way (MW), and even in the gas in a $z\approx 11$ galaxy. Using data from the APOGEE Data Release 17 and the Gaia Data Release 3, we present several independent lines of evidence that most of the MW's high-[N/O] stars were born in situ in massive bound clusters during the early, pre-disk evolution of the Galaxy. Specifically, we show that distributions of metallicity [Fe/H], energy, the angular momentum $L_z$, and distance of the low-metallicity high-[N/O] stars match the corresponding distributions of stars of the Aurora population and of the in-situ GCs. We also show that the fraction of in-situ field high-[N/O] stars, $f_{\rm N/O}$, increases rapidly with decreasing metallicity. During epochs when metallicity evolves from $\rm [Fe/H]=-1.5$ to $\rm [Fe/H]=-0.9$, the Galaxy spins up and transitions from a turbulent Aurora state to a coherently rotating disk. This transformation is accompanied by many qualitative changes. In particular, we show that high N/O abundances similar to those observed in GN-z11 were common before the spin-up ($\rm [Fe/H]\lesssim -1.5$) when up to $\approx 50\%-70\%$ of the in-situ stars formed in massive bound clusters. The dramatic drop of $f_{\rm N/O}$ at $\rm [Fe/H]\gtrsim -0.9$ indicates that after the disk emerges the fraction of stars forming in massive bound clusters decreases by two orders of magnitude.
△ Less
Submitted 31 May, 2023;
originally announced June 2023.
-
The Kinematics, Metallicities, and Orbits of Six Recently Discovered Galactic Star Clusters with Magellan/M2FS Spectroscopy
Authors:
Andrew B. Pace,
Sergey E. Koposov,
Matthew G. Walker,
Nelson Caldwell,
Mario Mateo,
Edward W. Olszewski,
Ian U. Roederer,
John I. Bailey III,
Vasily Belokurov,
Kyler Kuehn,
Ting S. Li,
Daniel B. Zucker
Abstract:
We present Magellan/M2FS spectroscopy of four recently discovered Milky Way star clusters (Gran 3/Patchick~125, Gran 4, Garro 01, LP 866) and two newly discovered open clusters (Gaia 9, Gaia 10) at low Galactic latitudes. We measure line-of-sight velocities and stellar parameters ([Fe/H], $\log{g}$, $T_{\rm eff}$, [Mg/Fe]) from high resolution spectroscopy centered on the Mg triplet and identify 2…
▽ More
We present Magellan/M2FS spectroscopy of four recently discovered Milky Way star clusters (Gran 3/Patchick~125, Gran 4, Garro 01, LP 866) and two newly discovered open clusters (Gaia 9, Gaia 10) at low Galactic latitudes. We measure line-of-sight velocities and stellar parameters ([Fe/H], $\log{g}$, $T_{\rm eff}$, [Mg/Fe]) from high resolution spectroscopy centered on the Mg triplet and identify 20-80 members per star cluster. We determine the kinematics and chemical properties of each cluster and measure the systemic proper motion and orbital properties by utilizing Gaia astrometry. We find Gran 3 to be an old, metal-poor (mean metallicity of [Fe/H]=-1.84) globular cluster located in the Galactic bulge on a retrograde orbit. Gran 4 is an old, metal-poor ([Fe/H]=-1.84) globular cluster with a halo-like orbit that happens to be passing through the Galactic plane. The orbital properties of Gran 4 are consistent with the proposed LMS-1/Wukong and/or Helmi streams merger events. Garro 01 is metal-rich ([Fe/H]=-0.30) and on a near circular orbit in the outer disk but its classification as an open cluster or globular cluster is ambiguous. . Gaia 9 and Gaia 10 are among the most distant known open clusters at $R_{GC}\sim 18, 21.2~kpc$ and most metal-poor with [Fe/H]~-0.50,-0.46 for Gaia 9 and Gaia 10, respectively. LP 866 is a nearby, metal-rich open cluster ([Fe/H]$=+0.1$). The discovery and confirmation of multiple star clusters in the Galactic plane shows the power of {\it Gaia} astrometry and the star cluster census remains incomplete.
△ Less
Submitted 8 September, 2023; v1 submitted 13 April, 2023;
originally announced April 2023.
-
Estimating the selection function of Gaia DR3 sub-samples
Authors:
A. Castro-Ginard,
A. G. A. Brown,
Z. Kostrzewa-Rutkowska,
T. Cantat-Gaudin,
R. Drimmel,
S. Oh,
V. Belokurov,
A. R. Casey,
M. Fouesneau,
S. Khanna,
A. M. Price-Whelan,
H. W. Rix
Abstract:
Understanding which sources are present in an astronomical catalogue and which are not is crucial for the accurate interpretation of astronomical data. In particular, for the multidimensional Gaia data, filters and cuts on different parameters or measurements introduces a selection function that may unintentionally alter scientific conclusions in subtle ways. We aim to develop a methodology to est…
▽ More
Understanding which sources are present in an astronomical catalogue and which are not is crucial for the accurate interpretation of astronomical data. In particular, for the multidimensional Gaia data, filters and cuts on different parameters or measurements introduces a selection function that may unintentionally alter scientific conclusions in subtle ways. We aim to develop a methodology to estimate the selection function for different sub-samples of stars in the Gaia catalogue. Comparing the number of stars in a given sub-sample to those in the overall Gaia catalogue, provides an estimate of the sub-sample membership probability, as a function of sky position, magnitude and colour. This estimate must differentiate the stochastic absence of sub-sample stars from selection effects. When multiplied with the overall Gaia catalogue selection function this provides the total selection function of the sub-sample. We present the method by estimating the selection function of the sources in Gaia DR3 with heliocentric radial velocity measurements. We also compute the selection function for the stars in the Gaia-Sausage/Enceladus sample, confirming that the apparent asymmetry of its debris across the sky is merely caused by selection effects. The developed method estimates the selection function of the stars present in a sub-sample of Gaia data, given that the sub-sample is completely contained in the Gaia parent catalogue (for which the selection function is known). This tool is made available in a GaiaUnlimited Python package.
△ Less
Submitted 30 March, 2023;
originally announced March 2023.
-
Stellar halo substructure generated by bar resonances
Authors:
Adam M. Dillamore,
Vasily Belokurov,
N. Wyn Evans,
Elliot Y. Davies
Abstract:
Using data from the Gaia satellite's Radial Velocity Spectrometer Data Release 3 (RVS, DR3), we find a new and robust feature in the phase space distribution of halo stars. It is a prominent ridge at constant energy and with angular momentum $L_z>0$. We run test particle simulations of a stellar halo-like distribution of particles in a realistic Milky Way potential with a rotating bar. We observe…
▽ More
Using data from the Gaia satellite's Radial Velocity Spectrometer Data Release 3 (RVS, DR3), we find a new and robust feature in the phase space distribution of halo stars. It is a prominent ridge at constant energy and with angular momentum $L_z>0$. We run test particle simulations of a stellar halo-like distribution of particles in a realistic Milky Way potential with a rotating bar. We observe similar structures generated in the simulations from the trapping of particles in resonances with the bar, particularly at the corotation resonance. Many of the orbits trapped at the resonances are halo-like, with large vertical excursions from the disc. The location of the observed structure in energy space is consistent with a bar pattern speed in the range $Ω_\mathrm{b}\approx35-40$ km s$^{-1}$ kpc$^{-1}$. Overall, the effect of the resonances is to give the inner stellar halo a mild, net spin in the direction of the bar's rotation. As the distribution of the angular momentum becomes asymmetric, a population of stars with positive mean $L_z$ and low vertical action is created. The variation of the average rotational velocity of the simulated stellar halo with radius is similar to the behaviour of metal-poor stars in data from the APOGEE survey. Though the effects of bar resonances have long been known in the Galactic disc, this is strong evidence that the bar can drive changes even in the diffuse and extended stellar halo through its resonances.
△ Less
Submitted 20 July, 2023; v1 submitted 28 February, 2023;
originally announced March 2023.
-
Peeking beneath the precision floor -- II. Probing the chemo-dynamical histories of the potential globular cluster siblings, NGC 288 and NGC 362
Authors:
Stephanie Monty,
David Yong,
Davide Massari,
Madeleine McKenzie,
GyuChul Myeong,
Sven Buder,
Amanda I. Karakas,
Ken C. Freeman,
Anna F. Marino,
Vasily Belokurov,
N. Wyn Evans
Abstract:
The assembly history of the Milky Way (MW) is a rapidly evolving subject, with numerous small accretion events and at least one major merger proposed in the MW's history. Accreted alongside these dwarf galaxies are globular clusters (GCs), which act as spatially coherent remnants of these past events. Using high precision differential abundance measurements from our recently published study, we in…
▽ More
The assembly history of the Milky Way (MW) is a rapidly evolving subject, with numerous small accretion events and at least one major merger proposed in the MW's history. Accreted alongside these dwarf galaxies are globular clusters (GCs), which act as spatially coherent remnants of these past events. Using high precision differential abundance measurements from our recently published study, we investigate the likelihood that the MW clusters NGC 362 and NGC 288 are galactic siblings, accreted as part of the Gaia-Sausage-Enceladus (GSE) merger. To do this, we compare the two GCs at the 0.01 dex level for 20+ elements for the first time. Strong similarities are found, with the two showing chemical similarity on the same order as those seen between the three LMC GCs, NGC 1786, NGC 2210 and NGC 2257. However, when comparing GC abundances directly to GSE stars, marked differences are observed. NGC 362 shows good agreement with GSE stars in the ratio of Eu to Mg and Si, as well as a clear dominance in the r- compared to the s-process, while NGC 288 exhibits only a slight r-process dominance. When fitting the two GC abundances with a GSE-like galactic chemical evolution model, NGC 362 shows agreement with both the model predictions and GSE abundance ratios (considering Si, Ni, Ba and Eu) at the same metallicity. This is not the case for NGC 288. We propose that the two are either not galactic siblings, or GSE was chemically inhomogeneous enough to birth two similar, but not identical clusters with distinct chemistry relative to constituent stars.
△ Less
Submitted 14 April, 2023; v1 submitted 13 February, 2023;
originally announced February 2023.
-
The 3D kinematics of stellar substructures in the periphery of the Large Magellanic Cloud
Authors:
Camila Navarrete,
David S. Aguado,
Vasily Belokurov,
Denis Erkal,
Alis Deason,
Lara Cullinane,
Julio Carballo-Bello
Abstract:
We report the 3D kinematics of 27 Mira-like stars in the northern, eastern and southern periphery of the Large Magellanic Cloud (LMC), based on Gaia proper motions and a dedicated spectroscopic follow-up. Low-resolution spectra were obtained for more than 40 Mira-like candidates, selected to trace known substructures in the LMC periphery. Radial velocities and stellar parameters were derived for a…
▽ More
We report the 3D kinematics of 27 Mira-like stars in the northern, eastern and southern periphery of the Large Magellanic Cloud (LMC), based on Gaia proper motions and a dedicated spectroscopic follow-up. Low-resolution spectra were obtained for more than 40 Mira-like candidates, selected to trace known substructures in the LMC periphery. Radial velocities and stellar parameters were derived for all stars. Gaia data release 3 astrometry and photometry were used to discard outliers, derive periods for those stars with available light curves, and determine their photometric chemical types. The 3D motion of the stars in the reference frame of the LMC revealed that most of the stars, in all directions, have velocities consistent with being part of the LMC disk population, out of equilibrium in the radial and vertical directions. A suite of N-body simulations was used to constrain the most likely past interaction history between the Clouds given the phase-space distribution of our targets. Model realizations in which the Small Magellanic Cloud (SMC) had three pericentric passages around the LMC best resemble the observations. The interaction history of those model realizations has a recent SMC pericentric passage ($\sim$320 Myr ago), preceded by an SMC crossing of the LMC disk at $\sim$0.97 Gyr ago, having a radial crossing distance of only $\sim$4.5 kpc. The previous disk crossing of the SMC was found to occur at $\sim$1.78 Gyr ago, with a much larger radial crossing distance of $\sim$10 kpc.
△ Less
Submitted 9 February, 2023;
originally announced February 2023.
-
Accelerated phase-mixing in the stellar halo due to a rotating bar
Authors:
Elliot Y. Davies,
Adam M. Dillamore,
Eugene Vasiliev,
Vasily Belokurov
Abstract:
In a galaxy merger, the stars tidally stripped from the satellite and accreted onto the host galaxy undergo phase mixing and form finely-grained structures in the phase space. However, these fragile structures may be destroyed in the subsequent galaxy evolution, in particular, by a rotating bar that appears well after the merger is completed. In this work, we investigate the survivability of phase…
▽ More
In a galaxy merger, the stars tidally stripped from the satellite and accreted onto the host galaxy undergo phase mixing and form finely-grained structures in the phase space. However, these fragile structures may be destroyed in the subsequent galaxy evolution, in particular, by a rotating bar that appears well after the merger is completed. In this work, we investigate the survivability of phase-space structures in the presence of a bar. We find that a bar with amplitude and pattern speed similar to those of the Milky Way would blur and destroy a substantial amount of the substructure that consists of particles with pericentre radii comparable to the bar length. While this appears to be in tension with the recent discovery of phase-space chevrons in \textit{Gaia} DR3 data, the most prominent chevrons in our simulations can still be recovered when applying the same analysis procedure as in observations. Moreover, the smoothing effect is less pronounced in the population of stars whose angular momenta have the opposite sign to the bar pattern speed.
△ Less
Submitted 3 March, 2023; v1 submitted 10 January, 2023;
originally announced January 2023.
-
The wide-field, multiplexed, spectroscopic facility WEAVE: Survey design, overview, and simulated implementation
Authors:
Shoko Jin,
Scott C. Trager,
Gavin B. Dalton,
J. Alfonso L. Aguerri,
J. E. Drew,
Jesús Falcón-Barroso,
Boris T. Gänsicke,
Vanessa Hill,
Angela Iovino,
Matthew M. Pieri,
Bianca M. Poggianti,
D. J. B. Smith,
Antonella Vallenari,
Don Carlos Abrams,
David S. Aguado,
Teresa Antoja,
Alfonso Aragón-Salamanca,
Yago Ascasibar,
Carine Babusiaux,
Marc Balcells,
R. Barrena,
Giuseppina Battaglia,
Vasily Belokurov,
Thomas Bensby,
Piercarlo Bonifacio
, et al. (190 additional authors not shown)
Abstract:
WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, will see first light in late 2022. WEAVE comprises a new 2-degree field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable 'mini' integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrogr…
▽ More
WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, will see first light in late 2022. WEAVE comprises a new 2-degree field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable 'mini' integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrograph covering the wavelength range 366$-$959\,nm at $R\sim5000$, or two shorter ranges at $R\sim20\,000$. After summarising the design and implementation of WEAVE and its data systems, we present the organisation, science drivers and design of a five- to seven-year programme of eight individual surveys to: (i) study our Galaxy's origins by completing Gaia's phase-space information, providing metallicities to its limiting magnitude for $\sim$3 million stars and detailed abundances for $\sim1.5$ million brighter field and open-cluster stars; (ii) survey $\sim0.4$ million Galactic-plane OBA stars, young stellar objects and nearby gas to understand the evolution of young stars and their environments; (iii) perform an extensive spectral survey of white dwarfs; (iv) survey $\sim400$ neutral-hydrogen-selected galaxies with the IFUs; (v) study properties and kinematics of stellar populations and ionised gas in $z<0.5$ cluster galaxies; (vi) survey stellar populations and kinematics in $\sim25\,000$ field galaxies at $0.3\lesssim z \lesssim 0.7$; (vii) study the cosmic evolution of accretion and star formation using $>1$ million spectra of LOFAR-selected radio sources; (viii) trace structures using intergalactic/circumgalactic gas at $z>2$. Finally, we describe the WEAVE Operational Rehearsals using the WEAVE Simulator.
△ Less
Submitted 31 October, 2023; v1 submitted 7 December, 2022;
originally announced December 2022.
-
The Pandora project. I: the impact of radiation and cosmic rays on baryonic and dark matter properties of dwarf galaxies
Authors:
Sergio Martin-Alvarez,
Debora Sijacki,
Martin G. Haehnelt,
Marion Farcy,
Yohan Dubois,
Vasily Belokurov,
Joakim Rosdahl,
Enrique Lopez-Rodriguez
Abstract:
Enshrouded in several well-known controversies, dwarf galaxies have been extensively studied to learn about the underlying cosmology, notwithstanding that physical processes regulating their properties are poorly understood. To shed light on these processes, we introduce the Pandora suite of 17 high-resolution (3.5 parsec half-cell side) dwarf galaxy formation cosmological simulations. Commencing…
▽ More
Enshrouded in several well-known controversies, dwarf galaxies have been extensively studied to learn about the underlying cosmology, notwithstanding that physical processes regulating their properties are poorly understood. To shed light on these processes, we introduce the Pandora suite of 17 high-resolution (3.5 parsec half-cell side) dwarf galaxy formation cosmological simulations. Commencing with thermo-turbulent star formation and mechanical supernova feedback, we gradually increase the complexity of physics incorporated leading to full-physics models combining magnetism, on-the-fly radiative transfer and the corresponding stellar photoheating, and SN-accelerated cosmic rays. We investigate combinations of these processes, comparing them with observations to constrain what are the main mechanisms determining dwarf galaxy properties. We find hydrodynamical `SN feedback-only' simulations struggle to produce realistic dwarf galaxies, leading either to overquenched or too centrally concentrated, dispersion dominated systems when compared to observed field dwarfs. Accounting for radiation with cosmic rays results in extended and rotationally-supported systems. Spatially `distributed' feedback leads to realistic stellar and HI masses as well as kinematics. Furthermore, resolved kinematic maps of our full-physics models predict kinematically distinct clumps and kinematic misalignments of stars, HI and HII after star formation events. Episodic star formation combined with its associated feedback induces more core-like dark matter central profiles, which our `SN feedback-only' models struggle to achieve. Our results demonstrate the complexity of physical processes required to capture realistic dwarf galaxy properties, making tangible predictions for integral field unit surveys, radio synchrotron emission, and for galaxy and multi-phase interstellar medium properties that JWST will probe.
△ Less
Submitted 25 September, 2023; v1 submitted 16 November, 2022;
originally announced November 2022.
-
A correlation between accreted stellar kinematics and dark matter halo spin in the ARTEMIS simulations
Authors:
Adam M. Dillamore,
Vasily Belokurov,
N. Wyn Evans,
Andreea S. Font
Abstract:
We report a correlation between the presence of a Gaia-Sausage-Enceladus (GSE) analogue and dark matter halo spin in the ARTEMIS simulations of Milky Way-like galaxies. The haloes which contain a large population of accreted stars on highly radial orbits (like the GSE) have lower spin on average than their counterparts with more isotropic stellar velocity distributions. The median modified spin pa…
▽ More
We report a correlation between the presence of a Gaia-Sausage-Enceladus (GSE) analogue and dark matter halo spin in the ARTEMIS simulations of Milky Way-like galaxies. The haloes which contain a large population of accreted stars on highly radial orbits (like the GSE) have lower spin on average than their counterparts with more isotropic stellar velocity distributions. The median modified spin parameters $λ^\prime$ differ by a factor of $\sim1.7$ at the present-day, with a similar value when the haloes far from virial equilibrium are removed. We also show that accreted stars make up a smaller proportion of the stellar populations in haloes containing a GSE analogue, and are stripped from satellites with stellar masses typically $\sim4$ times smaller. Our findings suggest that the higher spin of DM haloes without a GSE-like feature is due to mergers with large satellites of stellar mass $\sim10^{10}M_\odot$, which do not result in prominent radially anisotropic features like the GSE.
△ Less
Submitted 5 December, 2022; v1 submitted 2 November, 2022;
originally announced November 2022.
-
Ironing the folds: The phase space chevrons of a GSE-like merger as a dark matter subhalo detector
Authors:
Elliot Y. Davies,
Eugene Vasiliev,
Vasily Belokurov,
N. Wyn Evans,
Adam M. Dillamore
Abstract:
Recent work uncovered features in the phase space of the Milky Way's stellar halo which may be attributed to the last major merger. When stellar material from a satellite is accreted onto its host, it phase mixes and appears finely substructured in phase space. For a high-eccentricity merger, this substructure most clearly manifests as numerous wrapping chevrons in $(v_r, r)$ space, corresponding…
▽ More
Recent work uncovered features in the phase space of the Milky Way's stellar halo which may be attributed to the last major merger. When stellar material from a satellite is accreted onto its host, it phase mixes and appears finely substructured in phase space. For a high-eccentricity merger, this substructure most clearly manifests as numerous wrapping chevrons in $(v_r, r)$ space, corresponding to stripes in $(E, θ_r)$ space. We introduce the idea of using this substructure as an alternative subhalo detector to cold stellar streams. We simulate an N-body merger akin to the GSE and assess the impact of subhaloes on these chevrons. We examine how their deformation depends on the mass, pericentre, and number of subhaloes. To quantify the impact of perturbers, we utilise the appearance of chevrons in $(E, θ_r)$ space to introduce a new quantity -- the ironing parameter. We show that: (1) a single flyby of a massive ($\sim 10^{10}$ M$_{\odot}$) subhalo with pericentre comparable to, or within, the shell's apocentre smooths out the substructure, (2) a single flyby of a low mass ($\lesssim 10^8$ M$_{\odot}$) has negligible effect, (3) multiple flybys of subhalos derived from a subhalo mass function between $10^7-10^{10}$ M$_{\odot}$ cause significant damage if deep within the potential, (4) the effects of known perturbers (e.g. Sagittarius) should be detectable and offer constraints on their initial mass. The sensitivity to the populations of subhaloes suggests that we should be able to place an upper limit on the Milky Way's subhalo mass function.
△ Less
Submitted 5 December, 2022; v1 submitted 11 October, 2022;
originally announced October 2022.
-
Chemical and stellar properties of early-type dwarf galaxies around the Milky Way
Authors:
Vasily Belokurov,
N. Wyn Evans
Abstract:
Early-type dwarfs (ETDs) are the end points of the evolution of low-mass galaxies whose gas supply has been extinguished. The cessation of star-formation lays bare the ancient stellar populations. A wealth of information is stored in the colours, magnitudes, metallicities and abundances of resolved stars of the dwarf spheroidal and ultra-faint galaxies around the Milky Way, allowing their chemistr…
▽ More
Early-type dwarfs (ETDs) are the end points of the evolution of low-mass galaxies whose gas supply has been extinguished. The cessation of star-formation lays bare the ancient stellar populations. A wealth of information is stored in the colours, magnitudes, metallicities and abundances of resolved stars of the dwarf spheroidal and ultra-faint galaxies around the Milky Way, allowing their chemistry and stellar populations to be studied in great detail. Here, we summarize our current understanding, which has advanced rapidly over the last decade thanks to the flourishing of large-scale astrometric, photometric and spectroscopic surveys. We emphasise that the primeval stellar populations in the ETDs provide a unique laboratory to study the physical conditions on small scales at epochs beyond z=2. We highlight the observed diversity of star-formation and chemical enrichment histories in nearby dwarfs. These data can not yet be fully deciphered to reveal the key processes in the dwarf evolution but the first successful attempts have been made to pin down the sites of heavy element production.
△ Less
Submitted 19 September, 2022;
originally announced September 2022.
-
Energy wrinkles and phase-space folds of the last major merger
Authors:
Vasily Belokurov,
Eugene Vasiliev,
Alis J. Deason,
Sergey E. Koposov,
Azadeh Fattahi,
Adam M. Dillamore,
Elliot Y. Davies,
Robert J. J. Grand
Abstract:
Relying on the dramatic increase in the number of stars with full 6D phase-space information provided by the Gaia Data Release 3, we discover unambiguous signatures of phase-mixing in the stellar halo around the Sun. We show that for the stars likely belonging to the last massive merger, the (v_r,r) distribution contains a series of long and thin chevron-like overdensities. These phase-space sub-s…
▽ More
Relying on the dramatic increase in the number of stars with full 6D phase-space information provided by the Gaia Data Release 3, we discover unambiguous signatures of phase-mixing in the stellar halo around the Sun. We show that for the stars likely belonging to the last massive merger, the (v_r,r) distribution contains a series of long and thin chevron-like overdensities. These phase-space sub-structures are predicted to emerge following the dissolution of a satellite, when its tidal debris is given time to wind up, thin out and fold. Additionally, the observed energy and angular momentum (E, L_z) distribution appears more prograde at high energies, possibly revealing the original orbital angular momentum of the in-falling galaxy. The energy distribution of the debris is strongly asymmetric with a peak at low E -- which, we surmise, may be evidence of the dwarf's rapid sinking -- and riddled with wrinkles and bumps. If these small-scale energy inhomogeneities have been seeded during or immediately after the interaction with the Milky Way, and are not due to the spatial restriction of our study, then making use of the (v_r,r) chevrons to constrain the time of the merger becomes cumbersome. Nonetheless, we demonstrate that similar phase-space and (E,L_z) sub-structures are present in numerical simulations of galaxy interactions, both in bespoke N-body runs and in cosmological hydrodynamical zoom-in suites. The remnant traces of the progenitor's disruption and the signatures of the on-going phase-mixing discovered here will not only help to constrain the properties of our Galaxy's most important interaction, but also can be used as a novel tool to map out the Milky Way's current gravitational potential and its perturbations.
△ Less
Submitted 23 August, 2022;
originally announced August 2022.
-
An empirical model of the Gaia DR3 selection function
Authors:
Tristan Cantat-Gaudin,
Morgan Fouesneau,
Hans-Walter Rix,
Anthony G. A. Brown,
Alfred Castro-Ginard,
Ronald Drimmel,
David W. Hogg,
Andrew R. Casey,
Shourya Khanna,
Semyeong Oh,
Adrian M. Price Whelan,
Vasily Belokurov,
Andrew K. Saydjari,
Gregory M. Green
Abstract:
Interpreting and modelling astronomical catalogues requires an understanding of the catalogues' completeness or selection function: objects of what properties had a chance to end up in the catalogue. Here we set out to empirically quantify the completeness of the overall Gaia DR3 catalogue. This task is not straightforward because Gaia is the all-sky optical survey with the highest angular resolut…
▽ More
Interpreting and modelling astronomical catalogues requires an understanding of the catalogues' completeness or selection function: objects of what properties had a chance to end up in the catalogue. Here we set out to empirically quantify the completeness of the overall Gaia DR3 catalogue. This task is not straightforward because Gaia is the all-sky optical survey with the highest angular resolution to date and no consistent ``ground truth'' exists to allow direct comparisons.
However, well-characterised deeper imaging enables an empirical assessment of Gaia's $G$-band completeness across parts of the sky.
On this basis, we devised a simple analytical completeness model of Gaia as a function of the observed $G$ magnitude and position over the sky, which accounts for both the effects of crowding and the complex Gaia scanning law. Our model only depends on a single quantity: the median magnitude $M_{10}$ in a patch of the sky of catalogued sources with $\texttt{astrometric_matched_transits}$ $\leq 10$. $M_{10}$ reflects elementary completeness decisions in the Gaia pipeline and is computable from the Gaia DR3 catalogue itself and therefore applicable across the whole sky. We calibrate our model using the Dark Energy Camera Plane Survey (DECaPS) and test its predictions against Hubble Space Telescope observations of globular clusters. We find that our model predicts Gaia's completeness values to a few per cent across the sky. We make the model available as a part of the $\texttt{gaiasf}$ Python package built and maintained by the GaiaUnlimited project: $\texttt{https://github.com/gaia-unlimited/gaiaunlimited}$
△ Less
Submitted 6 September, 2022; v1 submitted 19 August, 2022;
originally announced August 2022.
-
Overview of the DESI Milky Way Survey
Authors:
Andrew P. Cooper,
Sergey E. Koposov,
Carlos Allende Prieto,
Christopher J. Manser,
Namitha Kizhuprakkat,
Adam D. Myers,
Arjun Dey,
Boris T. Gaensicke,
Ting S. Li,
Constance Rockosi,
Monica Valluri,
Joan Najita,
Alis Deason,
Anand Raichoor,
Mei-Yu Wang,
Yuan-Sen Ting,
Bokyoung Kim,
Andreia Carrillo,
Wenting Wang,
Leandro Beraldo e Silva,
Jiwon Jesse Han,
Jiani Ding,
Miguel Sanchez-Conde,
Jessica N. Aguilar,
Steven Ahlen
, et al. (40 additional authors not shown)
Abstract:
We describe the Milky Way Survey (MWS) that will be undertaken with the Dark Energy Spectroscopic Instrument (DESI) on the Mayall 4m telescope at the Kitt Peak National Observatory. Over the next 5 yr DESI MWS will observe approximately seven million stars at Galactic latitudes |b|>20 degrees, with an inclusive target selection scheme focused on the thick disk and stellar halo. MWS will also inclu…
▽ More
We describe the Milky Way Survey (MWS) that will be undertaken with the Dark Energy Spectroscopic Instrument (DESI) on the Mayall 4m telescope at the Kitt Peak National Observatory. Over the next 5 yr DESI MWS will observe approximately seven million stars at Galactic latitudes |b|>20 degrees, with an inclusive target selection scheme focused on the thick disk and stellar halo. MWS will also include several high-completeness samples of rare stellar types, including white dwarfs, low-mass stars within 100pc of the Sun, and horizontal branch stars. We summarize the potential of DESI to advance understanding of Galactic structure and stellar evolution. We introduce the final definitions of the main MWS target classes and estimate the number of stars in each class that will be observed. We describe our pipelines for deriving radial velocities, atmospheric parameters, and chemical abundances. We use ~500,000 spectra of unique stellar targets from the DESI Survey Validation program (SV) to demonstrate that our pipelines can measure radial velocities to ~1 km/s and [Fe/H] accurate to ~0.2 dex for typical stars in our main sample. We find the stellar parameter distributions from ~100 sq. deg of SV observations with >90% completeness on our main sample are in good agreement with expectations from mock catalogs and previous surveys.
△ Less
Submitted 20 February, 2023; v1 submitted 17 August, 2022;
originally announced August 2022.
-
RR Lyrae From Binary Evolution: Abundant, Young and Metal-Rich
Authors:
Alexey Bobrick,
Giuliano Iorio,
Vasily Belokurov,
Joris Vos,
Maja Vuckovic,
Nicola Giacobbo
Abstract:
RR Lyrae are a well-known class of pulsating horizontal branch stars widely used as tracers of old, metal-poor stellar populations. However, mounting observational evidence shows that a significant fraction of these stars may be young and metal-rich. Here, through detailed binary stellar evolution modelling, we show that all such metal-rich RR Lyrae can be naturally produced through binary interac…
▽ More
RR Lyrae are a well-known class of pulsating horizontal branch stars widely used as tracers of old, metal-poor stellar populations. However, mounting observational evidence shows that a significant fraction of these stars may be young and metal-rich. Here, through detailed binary stellar evolution modelling, we show that all such metal-rich RR Lyrae can be naturally produced through binary interactions. Binary companions of these RR Lyrae stars formed through binary interactions partly strip their progenitor's envelopes during a preceding red giant phase. As a result, stripped horizontal branch stars become bluer than their isolated stellar evolution counterparts and thus end up in the instability strip. In contrast, in the single evolution scenario, the stars can attain such colours only at large age and low metallicity. While binary-made RR Lyrae can possess any ages and metallicities, their Galactic population is relatively young (1 to 9 Gyr) and dominated by the Thin Disc and the Bulge. We show that Galactic RR Lyrae from binary evolution are produced at rates compatible with the observed metal-rich population and have consistent G-band magnitudes, Galactic kinematics and pulsation properties. Furthermore, these systems dominate the RR Lyrae population in the Solar Neighbourhood. We predict that all metal-rich RR Lyrae have an A, F, G or K-type companion with a long orbital period (P > 1000 d). Observationally characterising the orbital periods and masses of such stellar companions will provide valuable new constraints on mass and angular momentum-loss efficiency for Sun-like accretors and the nature of RR Lyrae populations.
△ Less
Submitted 13 January, 2024; v1 submitted 8 August, 2022;
originally announced August 2022.
-
Milky Way's eccentric constituents with $Gaia$, APOGEE $\&$ GALAH
Authors:
G. C. Myeong,
Vasily Belokurov,
David S. Aguado,
N. Wyn Evans,
Nelson Caldwell,
James Bradley
Abstract:
We report the results of an unsupervised decomposition of the local stellar halo in the chemo-dynamical space spanned by the abundance measurements from APOGEE DR17 and GALAH DR3. In our Gaussian Mixture Model, only four independent components dominate the halo in the Solar neighborhood, three previously known $Aurora$, $Splash$ and Gaia-Sausage/Enceladus (GS/E) and one new, $Eos$. Only one of the…
▽ More
We report the results of an unsupervised decomposition of the local stellar halo in the chemo-dynamical space spanned by the abundance measurements from APOGEE DR17 and GALAH DR3. In our Gaussian Mixture Model, only four independent components dominate the halo in the Solar neighborhood, three previously known $Aurora$, $Splash$ and Gaia-Sausage/Enceladus (GS/E) and one new, $Eos$. Only one of these four is of accreted origin, namely the GS/E, thus supporting the earlier claims that the GS/E is the main progenitor of the Galactic stellar halo. We show that $Aurora$ is entirely consistent with the chemical properties of the so-called Heracles merger. In our analysis in which no predefined chemical selection cuts are applied, $Aurora$ spans a wide range of [Al/Fe] with a metallicity correlation indicative of a fast chemical enrichment in a massive galaxy, the young Milky Way. The new halo component dubbed $Eos$ is classified as $in situ$ given its high mean [Al/Fe]. $Eos$ shows strong evolution as a function of [Fe/H], where it changes from being the closest to GS/E at its lowest [Fe/H] to being indistinguishable from the Galactic low-$α$ population at its highest [Fe/H]. We surmise that at least some of the outer thin disk of the Galaxy started its evolution in the gas polluted by the GS/E, and $Eos$ is an evidence of this process.
△ Less
Submitted 30 January, 2023; v1 submitted 15 June, 2022;
originally announced June 2022.
-
The stellar halo in Local Group Hestia simulations III. Chemical abundance relations for accreted and in-situ stars
Authors:
Sergey Khoperskov,
Ivan Minchev,
Noam Libeskind,
Vasily Belokurov,
Matthias Steinmetz,
Facundo A. Gomez,
Robert J. J. Grand,
Yehuda Hoffman,
Alexander Knebe,
Jenny G. Sorce,
Martin Sparre,
Elmo Tempel,
Mark Vogelsberger
Abstract:
Since the chemical abundances of stars are the fossil records of the physical conditions in galaxies, they provide the key information for recovering the assembly history of galaxies. In this work, we explore the chemo-chrono-kinematics of accreted and in-situ stars, by analyzing six M31/MW analogues from the HESTIA suite of cosmological hydrodynamics zoom-in simulations of the Local Group. We fou…
▽ More
Since the chemical abundances of stars are the fossil records of the physical conditions in galaxies, they provide the key information for recovering the assembly history of galaxies. In this work, we explore the chemo-chrono-kinematics of accreted and in-situ stars, by analyzing six M31/MW analogues from the HESTIA suite of cosmological hydrodynamics zoom-in simulations of the Local Group. We found that the merger debris are chemically distinct from the survived dwarf galaxies. The mergers debris have abundances expected for stars originating from dwarfs that had their star formation activity quenched at early times. Accreted stellar haloes, including individual debris, reveal abundance gradients in the ELz, where the most metal-rich stars have formed in the inner parts of the disrupted systems before the merger and mainly contribute to the central regions of the hosts. Therefore, we suggest that abundance measurements in the inner MW will allow constraining better the parameters of building blocks of the MW stellar halo. The MDFs of the individual debris show several peaks and the majority of debris have lower metallicity than the in-situ stars for Lz>0, while non-rotating and retrograde accreted stars are similar to the in-situ. Prograde accreted stars show a prominent knee in the [Fe/H]-[Mg/Fe] plane while the retrograde stars typically deposit to a high-[Mg/Fe] sequence. We found that the metal-poor stars ([Fe/H]<-1) of the HESTIA galaxies exhibit between zero to 80 km/s net rotation which is consistent with the Aurora population. At higher metallicities, we detect a sharp transition (spin-up) from the turbulent phase to a disk-like rotation. Mergers debris are similar in the [Fe/H]-[Mg/Fe] plane. However, combining a set of abundances allows to capture chemical patterns corresponding to different debris, which are the most prominent as a function of stellar age.
△ Less
Submitted 12 September, 2023; v1 submitted 11 June, 2022;
originally announced June 2022.
-
An Approach to Quantum 2D Gravity
Authors:
Vladimir V. Belokurov,
Evgeniy T. Shavgulidze
Abstract:
We consider a model of 2D gravity with the action quadratic in curvature and represent path integrals as integrals over the SL(2, R) invariant Gaussian functional measure. We reduce these path integrals to the products of Wiener path integrals and calculate the correlation function of the metric in the first perturbative order.
We consider a model of 2D gravity with the action quadratic in curvature and represent path integrals as integrals over the SL(2, R) invariant Gaussian functional measure. We reduce these path integrals to the products of Wiener path integrals and calculate the correlation function of the metric in the first perturbative order.
△ Less
Submitted 10 June, 2022;
originally announced June 2022.
-
The stellar halo in Local Group Hestia simulations II. The accreted component
Authors:
Sergey Khoperskov,
Ivan Minchev,
Noam Libeskind,
Misha Haywood,
Paola Di Matteo,
Vasily Belokurov,
Matthias Steinmetz,
Facundo A. Gomez,
Robert J. J. Grand,
Yehuda Hoffman,
Alexander Knebe,
Jenny G. Sorce,
Martin Sparre,
Elmo Tempel,
Mark Vogelsberger
Abstract:
In the Milky Way, recent progress in the exploration of its assembly history is driven by the tremendous amount of high-quality data delivered by Gaia, which has revealed a number of substructures potentially linked to several ancient accretion events. In this work, aiming to explore the phase-space structure of accreted stars, we analyze six M31/MW analogues from the HESTIA suite of cosmological…
▽ More
In the Milky Way, recent progress in the exploration of its assembly history is driven by the tremendous amount of high-quality data delivered by Gaia, which has revealed a number of substructures potentially linked to several ancient accretion events. In this work, aiming to explore the phase-space structure of accreted stars, we analyze six M31/MW analogues from the HESTIA suite of cosmological hydrodynamics zoom-in simulations of the Local Group. We found that all the HESTIA galaxies experience a few dozen mergers but only 1-4 mergers have the stellar mass ratio >0.2 where, depending on the halo definition, the most massive merger contributes from 20% to 70% of the total stellar halo. Individual merger remnants show diverse density distributions at z=0, significantly overlapping with each other and with the in-situ stars in the ELz, UV and RVphi coordinates. The mergers debris often change their position in the ELz with time due to the galactic mass growth and the non-axisymmetry of the potential. In agreement with previous works, we show that even individual merger debris exhibit a number of distinct ELz features. In the UV plane, all HESTIA galaxies reveal radially hot, non-rotating or weakly counter-rotating, Gaia-Sausage-like features. We found an age gradient in Elz space for the individual debris, where the youngest stars, formed in the inner regions of accreting systems, deposit to the innermost regions of the host. The bulk of these stars is being formed during the last stages of accretion, making it possible to date the merger. In actions space (Jr, Jz, Jφ), the mergers debris do not appear as isolated substructures but are instead scattered over a large parameters area and overlapping with the in-situ stars. We also introduce a purely kinematic space (Jz/Jr-eccentricity), where different merger debris can be disentangled better from each other and from the in-situ stars.
△ Less
Submitted 12 September, 2023; v1 submitted 9 June, 2022;
originally announced June 2022.
-
The stellar halo in Local Group Hestia simulations I. The in-situ component and the effect of mergers
Authors:
Sergey Khoperskov,
Ivan Minchev,
Noam Libeskind,
Misha Haywood,
Paola Di Matteo,
Vasily Belokurov,
Matthias Steinmetz,
Facundo A. Gomez,
Robert J. J. Grand,
Yehuda Hoffman,
Alexander Knebe,
Jenny G. Sorce,
Martin Sparre,
Elmo Tempel,
Mark Vogelsberger
Abstract:
Theory suggests that mergers play an important role in shaping galactic discs and stellar haloes, which was observationally confirmed in the MW thanks to Gaia data. In this work, aiming to probe the contribution of mergers to the in situ stellar halo formation, we analyse six M31/MW analogues from the HESTIA suite of cosmological hydrodynamical zoom-in simulations of the LG. We found that all the…
▽ More
Theory suggests that mergers play an important role in shaping galactic discs and stellar haloes, which was observationally confirmed in the MW thanks to Gaia data. In this work, aiming to probe the contribution of mergers to the in situ stellar halo formation, we analyse six M31/MW analogues from the HESTIA suite of cosmological hydrodynamical zoom-in simulations of the LG. We found that all the HESTIA galaxies experience between one to four mergers with stellar mass ratios between 0.2 and 1 relative to the host at the time of the merger. These significant mergers, with a single exception, happened 7-11Gyr ago. The overall impact of the most massive mergers in HESTIA is clearly seen as a sharp increase in the orbital eccentricity (and a corresponding decrease in the rotational velocity Vphi of pre-existing disc stars of the main progenitor, thus nicely reproducing the Splash-, Plume-like feature that was discovered in the MW. We do find a correlation between mergers and close pericentric passages of massive satellites and bursts of star formation in the in situ component. Massive mergers sharply increase the disc velocity dispersion of the in situ stars; however, the latest significant merger often heats up the disc up to the numbers when the contribution of the previous ones is less prominent in the age-velocity dispersion relation. In HESTIA galaxies, the in situ halo is an important component of the inner stellar halo where its fraction is about 30-40%, while in the outer parts it typically does not exceed ~5% beyond 15 kpc. The simulations suggest that this component of the stellar haloes continues to grow well after mergers conclude; however, the most significant contribution comes from stars that formed recently before the merger. The orbital analysis of the HESTIA galaxies suggests that wedges in Rmax-Zmax space are mainly populated by the stars born between significant mergers.
△ Less
Submitted 12 September, 2023; v1 submitted 9 June, 2022;
originally announced June 2022.
-
Binary parameters from astrometric and spectroscopic errors--candidate hierarchical triples and massive dark companions in Gaia DR3
Authors:
Shion Andrew,
Zephyr Penoyre,
Vasily Belokurov,
N. Wyn Evans,
Semyeong Oh
Abstract:
We show how astrometric and spectroscopic errors introduced by an unresolved binary system can be combined to give estimates of the binary period and mass ratio. This can be performed analytically if we assume we see one or more full orbits over our observational baseline, or numerically for all other cases. We apply this method to Gaia DR3 data, combining the most recent astrometric and spectrosc…
▽ More
We show how astrometric and spectroscopic errors introduced by an unresolved binary system can be combined to give estimates of the binary period and mass ratio. This can be performed analytically if we assume we see one or more full orbits over our observational baseline, or numerically for all other cases. We apply this method to Gaia DR3 data, combining the most recent astrometric and spectroscopic data. We compare inferred periods and mass ratios calculated using our method with orbital parameters measured for non-single stars in Gaia DR3 and find good agreement. Finally, we use this method to search the subset of the Gaia DR3 RVS dataset with rv_method_used=1 for compact object candidates. We select sources with significant astrometric and spectroscopic errors ($RUWE_{ast}>1.25$ and $RUWE_{spec}>2$), with large inferred mass ratios, and large inferred companion masses ($q>1$ and $M_c>3 M_\odot$) giving a catalogue of 4,641 candidate hierarchical triples and Main Sequence+Compact Object pairs. We apply more stringent cuts, and impose low levels of photometric variability to remove likely triples ($RUWE_{phot}<2$) producing a gold sample of 45 candidates.
△ Less
Submitted 8 October, 2022; v1 submitted 9 June, 2022;
originally announced June 2022.
-
First semi-empirical test of the white dwarf mass-radius relationship using a single white dwarf via astrometric microlensing
Authors:
Peter McGill,
Jay Anderson,
Stefano Casertano,
Kailash C. Sahu,
Pierre Bergeron,
Simon Blouin,
Patrick Dufour,
Leigh C. Smith,
N. Wyn Evans,
Vasily Belokurov,
Richard L. Smart,
Andrea Bellini,
Annalisa Calamida,
Martin Dominik,
Noé Kains,
Jonas Klüter,
Martin Bo Nielsen,
Joachim Wambsganss
Abstract:
In November 2019, the nearby single, isolated DQ-type white dwarf LAWD 37 (WD 1142-645) aligned closely with a distant background source and caused an astrometric microlensing event. Leveraging astrometry from \Gaia{} and followup data from the \textit{Hubble Space Telescope} we measure the astrometric deflection of the background source and obtain a gravitational mass for LAWD~37. The main challe…
▽ More
In November 2019, the nearby single, isolated DQ-type white dwarf LAWD 37 (WD 1142-645) aligned closely with a distant background source and caused an astrometric microlensing event. Leveraging astrometry from \Gaia{} and followup data from the \textit{Hubble Space Telescope} we measure the astrometric deflection of the background source and obtain a gravitational mass for LAWD~37. The main challenge of this analysis is in extracting the lensing signal of the faint background source whilst it is buried in the wings of LAWD~37's point spread function. Removal of LAWD 37's point spread function induces a significant amount of correlated noise which we find can mimic the astrometric lensing signal. We find a deflection model including correlated noise caused by the removal of LAWD~37's point spread function best explains the data and yields a mass for LAWD 37 of $0.56\pm0.08 M_{\odot}$. This mass is in agreement with the theoretical mass-radius relationship and cooling tracks expected for CO core white dwarfs. Furthermore, the mass is consistent with no or trace amounts of hydrogen that is expected for objects with helium-rich atmospheres like LAWD 37. We conclude that further astrometric followup data on the source is likely to improve the inference on LAWD 37's mass at the $\approx3$ percent level and definitively rule out purely correlated noise explanations of the data. This work provides the first semi-empirical test of the white dwarf mass-radius relationship using a single, isolated white dwarf and supports current model atmospheres of DQ white dwarfs and white dwarf evolutionary theory.
△ Less
Submitted 28 January, 2023; v1 submitted 3 June, 2022;
originally announced June 2022.
-
On the Hunt for the Origins of the Orphan--Chenab Stream: Detailed Element Abundances with APOGEE and Gaia
Authors:
Keith Hawkins,
Adrian M. Price-Whelan,
Allyson A. Sheffield,
Aidan Z. Subrahimovic,
Rachael L. Beaton,
Vasily Belokurov,
Denis Erkal,
Sergey E. Koposov,
Richard R. Lane,
Chervin F. P. Laporte,
Christian Nitschelm
Abstract:
Stellar streams in the Galactic halo are useful probes of the assembly of galaxies like the Milky Way. Many tidal stellar streams that have been found in recent years are accompanied by a known progenitor globular cluster or dwarf galaxy. However, the Orphan--Chenab (OC) stream is one case where a relatively narrow stream of stars has been found without a known progenitor. In an effort to find the…
▽ More
Stellar streams in the Galactic halo are useful probes of the assembly of galaxies like the Milky Way. Many tidal stellar streams that have been found in recent years are accompanied by a known progenitor globular cluster or dwarf galaxy. However, the Orphan--Chenab (OC) stream is one case where a relatively narrow stream of stars has been found without a known progenitor. In an effort to find the parent of the OC stream, we use astrometry from the early third data release of ESA's Gaia mission (Gaia EDR3) and radial velocity information from the SDSS-IV APOGEE survey to find up to 13 stars that are likely members of the OC stream. We use the APOGEE survey to study the chemical nature (for up to 13 stars) of the OC stream in the $α$ (O, Mg, Ca, Si, Ti, S), odd-Z (Al, K, V), Fe-peak (Fe, Ni, Mn, Co, Cr) and neutron capture (Ce) elemental groups. We find that the stars that make up the OC stream are not consistent with a mono-metallic population and have a median metallicity of --1.92~dex with a dispersion of 0.28 dex. Our results also indicate that the $α$-elements are depleted compared to the known Milky Way populations and that its [Mg/Al] abundance ratio is not consistent with second generation stars from globular clusters. The detailed chemical pattern of these stars indicates that the OC stream progenitor is very likely to be a dwarf spheroidal galaxy with a mass of ~10$^6$ M$_\odot$.
△ Less
Submitted 27 May, 2022;
originally announced May 2022.
-
The impact of a massive Sagittarius dSph on GD-1-like streams
Authors:
Adam M. Dillamore,
Vasily Belokurov,
N. Wyn Evans,
Adrian M. Price-Whelan
Abstract:
We investigate the effects of a massive ($\gtrsim4\times10^{10}M_\odot$) Sagittarius dwarf spheroidal galaxy (Sgr) on stellar streams using test particle simulations in a realistic Milky Way potential. We find that Sgr can easily disrupt streams formed more than $\sim3$ Gyr ago, while stars stripped more recently are generally unaffected. In certain realizations, Sgr is able to produce asymmetry b…
▽ More
We investigate the effects of a massive ($\gtrsim4\times10^{10}M_\odot$) Sagittarius dwarf spheroidal galaxy (Sgr) on stellar streams using test particle simulations in a realistic Milky Way potential. We find that Sgr can easily disrupt streams formed more than $\sim3$ Gyr ago, while stars stripped more recently are generally unaffected. In certain realizations, Sgr is able to produce asymmetry between the leading and trailing tails of Pal 5, qualitatively similar to observations. Using data from the Gaia space telescope and elsewhere, we fit models to the GD-1 stream in the presence of a Sgr with various initial masses. While the best-fitting models do show perturbations resulting from interactions with Sgr, we find that the level of disruption is not significantly greater than in the observed stream. To investigate the general effects of Sgr on a population of streams, we generate 1000 mock streams on GD-1-like orbits with randomized orientations. Some streams show clear evidence of disruption, becoming folded on the sky or developing asymmetry betweeen their two tails. However, many survive unaffected and the peak surface brightness of stars is decreased by no more than $\sim0.3$ mag/arcsec$^2$ on average. We conclude that Sgr having an initial mass of $\gtrsim4\times10^{10}M_\odot$ is compatible with the survival and detection of streams formed more than 3 Gyr ago.
△ Less
Submitted 1 September, 2022; v1 submitted 26 May, 2022;
originally announced May 2022.
-
Information content of BP/RP spectra in Gaia DR3
Authors:
Callum E. C. Witten,
David S. Aguado,
Jason L. Sanders,
Vasily Belokurov,
N. Wyn Evans,
Sergey E. Koposov,
Carlos Allende Prieto,
Francesca De Angeli,
Mike J. Irwin
Abstract:
Gaia Data Release 3 has provided the astronomical community with the largest stellar spectroscopic survey to date ($>$ 220 million sources). The low resolution (R$\sim$50) blue photometer (BP) and red photometer (RP) spectra will allow for the estimation of stellar atmospheric parameters such as effective temperature, surface gravity and metallicity. We create mock Gaia BP/RP spectra and use Fishe…
▽ More
Gaia Data Release 3 has provided the astronomical community with the largest stellar spectroscopic survey to date ($>$ 220 million sources). The low resolution (R$\sim$50) blue photometer (BP) and red photometer (RP) spectra will allow for the estimation of stellar atmospheric parameters such as effective temperature, surface gravity and metallicity. We create mock Gaia BP/RP spectra and use Fisher information matrices to probe the resolution limit of stellar parameter measurements using BP/RP spectra. The best-case scenario uncertainties that this analysis provides are then used to produce a mock-observed stellar population in order to evaluate the false positive rate (FPR) of identifying extremely metal-poor (EMP) stars. We conclude that the community will be able to confidently identify metal-poor stars at magnitudes brighter than $G = 16$ using BP/RP spectra. At fainter magnitudes true detections will start to be overwhelmed by false positives. When adopting the commonly-used $G < 14$ limit for metal-poor star searches, we find a FPR for the low-metallicity regimes [Fe/H] $<$ -2, -2.5 and -3 of just 14$\%$, 33$\%$ and 56$\%$ respectively, offering the potential for significant improvements on previous targeting campaigns. Additionally, we explore the chemical sensitivity obtainable directly from BP/RP spectra for Carbon and $α$-elements. We find an absolute Carbon abundance uncertainty of $σ_{A(C)} < 1$ dex for Carbon-enriched metal-poor (CEMP) stars, indicating the potential to identify a CEMP stellar population for follow-up confirmation with higher resolution spectroscopy. Finally, we find that large uncertainties in $α$-element abundance measurements using BP/RP spectra means that efficiently obtaining these abundances will be challenging.
△ Less
Submitted 10 August, 2022; v1 submitted 24 May, 2022;
originally announced May 2022.
-
Charting Galactic Accelerations with Stellar Streams and Machine Learning
Authors:
Jacob Nibauer,
Vasily Belokurov,
Miles Cranmer,
Jeremy Goodman,
Shirley Ho
Abstract:
We present a data-driven method for reconstructing the galactic acceleration field from phase-space measurements of stellar streams. Our approach is based on a flexible and differentiable fit to the stream in phase-space, enabling a direct estimate of the acceleration vector along the stream. Reconstruction of the local acceleration field can be applied independently to each of several streams, al…
▽ More
We present a data-driven method for reconstructing the galactic acceleration field from phase-space measurements of stellar streams. Our approach is based on a flexible and differentiable fit to the stream in phase-space, enabling a direct estimate of the acceleration vector along the stream. Reconstruction of the local acceleration field can be applied independently to each of several streams, allowing us to sample the acceleration field due to the underlying galactic potential across a range of scales. Our approach is methodologically different from previous works, since a model for the gravitational potential does not need to be adopted beforehand. Instead, our flexible neural-network-based model treats the stream as a collection of orbits with a locally similar mixture of energies, rather than assuming that the stream delineates a single stellar orbit. Accordingly, our approach allows for distinct regions of the stream to have different mean energies, as is the case for real stellar streams. Once the acceleration vector is sampled along the stream, standard analytic models for the galactic potential can then be rapidly constrained. We find our method recovers the correct parameters for a ground-truth triaxial logarithmic halo potential when applied to simulated stellar streams. Alternatively, we demonstrate that a flexible potential can be constrained with a neural network, though standard multipole expansions can also be constrained. Our approach is applicable to simple and complicated gravitational potentials alike, and enables potential reconstruction from a fully data-driven standpoint using measurements of slowly phase-mixing tidal debris.
△ Less
Submitted 30 March, 2023; v1 submitted 24 May, 2022;
originally announced May 2022.
-
The Magellanic Edges Survey III. Kinematics of the disturbed LMC outskirts
Authors:
L. R. Cullinane,
A. D. Mackey,
G. S. Da Costa,
D. Erkal,
S. E. Koposov,
V. Belokurov
Abstract:
We explore the structural and kinematic properties of the outskirts of the Large Magellanic Cloud (LMC) using data from the Magellanic Edges Survey (MagES) and Gaia EDR3. Even at large galactocentric radii ($8^\circ<R<11^\circ$), we find the north-eastern LMC disk is relatively unperturbed: its kinematics are consistent with a disk of inclination ~$36.5^\circ$ and line-of-nodes position angle ~…
▽ More
We explore the structural and kinematic properties of the outskirts of the Large Magellanic Cloud (LMC) using data from the Magellanic Edges Survey (MagES) and Gaia EDR3. Even at large galactocentric radii ($8^\circ<R<11^\circ$), we find the north-eastern LMC disk is relatively unperturbed: its kinematics are consistent with a disk of inclination ~$36.5^\circ$ and line-of-nodes position angle ~$145^\circ$ east of north. In contrast, fields at similar radii in the southern and western disk are significantly perturbed from equilibrium, with non-zero radial and vertical velocities, and distances significantly in front of the disk plane implied by our north-eastern fields. We compare our observations to simple dynamical models of the Magellanic/Milky Way system which describe the LMC as a collection of tracer particles within a rigid potential, and the Small Magellanic Cloud (SMC) as a rigid Hernquist potential. A possible SMC crossing of the LMC disk plane ~400 Myr ago, in combination with the LMC's infall to the Milky Way potential, can qualitatively explain many of the perturbations in the outer disk. Additionally, we find the claw-like and arm-like structures south of the LMC have similar metallicities to the outer LMC disk ([Fe/H]~-1), and are likely comprised of perturbed LMC disk material. The claw-like substructure is particularly disturbed, with out-of-plane velocities >60 km s$^{-1}$ and apparent counter-rotation relative to the LMC's disk motion. More detailed N-body models are necessary to elucidate the origin of these southern features, potentially requiring repeated interactions with the SMC prior to ~1 Gyr ago.
△ Less
Submitted 10 March, 2022;
originally announced March 2022.
-
From dawn till disk: Milky Way's turbulent youth revealed by the APOGEE+Gaia data
Authors:
Vasily Belokurov,
Andrey Kravtsov
Abstract:
We use accurate estimates of aluminium abundance provided as part of the APOGEE Data Release 17 and Gaia Early Data Release 3 astrometry to select a highly pure sample of stars with metallicity $-1.5\lesssim {\rm [Fe/H]}\lesssim 0.5$ born in-situ in the Milky Way proper. We show that the low-metallicity ([Fe/H]$\lesssim -1.3$) in-situ component that we dub Aurora is kinematically hot with an appro…
▽ More
We use accurate estimates of aluminium abundance provided as part of the APOGEE Data Release 17 and Gaia Early Data Release 3 astrometry to select a highly pure sample of stars with metallicity $-1.5\lesssim {\rm [Fe/H]}\lesssim 0.5$ born in-situ in the Milky Way proper. We show that the low-metallicity ([Fe/H]$\lesssim -1.3$) in-situ component that we dub Aurora is kinematically hot with an approximately isotropic velocity ellipsoid and a modest net rotation. Aurora stars exhibit large scatter in metallicity and in a number of element abundance ratios. The median tangential velocity of the in-situ stars increases sharply with increasing metallicity between [Fe/H]$=-1.3$ and $-0.9$, the transition that we call the spin-up. The observed and theoretically expected age-metallicity correlations imply that this increase reflects a rapid formation of the Milky Way disk over $\approx 1-2$ Gyrs. The transformation of the stellar kinematics as a function of [Fe/H] is accompanied by a qualitative change in chemical abundances: the scatter drops sharply once the Galaxy builds up a disk during later epochs corresponding to [Fe/H]$>-0.9$. Results of galaxy formation models presented in this and other recent studies strongly indicate that the trends observed in the Milky Way reflect generic processes during the early evolution of progenitors of MW-sized galaxies: a period of chaotic pre-disk evolution, when gas is accreted along cold narrow filaments and when stars are born in irregular configurations, and subsequent rapid disk formation. The latter signals formation of a stable hot gaseous halo around the MW progenitor, which changes the mode of gas accretion and allows development of coherently rotating disk.
△ Less
Submitted 9 March, 2022;
originally announced March 2022.
-
A gap in the double white dwarf separation distribution caused by the common-envelope evolution: astrometric evidence from Gaia
Authors:
Valeriya Korol,
Vasily Belokurov,
Silvia Toonen
Abstract:
The trajectory of the center of light of an unresolved binary is different from that of its center of mass. Binary-induced stellar centroid wobbling can therefore be detected as an excess in the goodness-of-fit of the single-star astrometric model. We use reduced $χ^2$ of the astrometric fit in the {\it Gaia} Early Data Release 3 to detect the likely unresolved double white dwarfs (DWDs). Using pa…
▽ More
The trajectory of the center of light of an unresolved binary is different from that of its center of mass. Binary-induced stellar centroid wobbling can therefore be detected as an excess in the goodness-of-fit of the single-star astrometric model. We use reduced $χ^2$ of the astrometric fit in the {\it Gaia} Early Data Release 3 to detect the likely unresolved double white dwarfs (DWDs). Using parallax-based distances we convert the excess of reduced $χ^2$ into the amplitude of the centroid wobble $δa$, which is proportional to the binary separation $a$. The measured $δa$ distribution drops towards larger wobble amplitudes and shows a break around $δa \approx 0.2$ where it steepens. The integral of the distribution yields DWD fraction of $6.5 \pm 3.7$ per cent in the range $0.01 < a (\text{au}) < 2$. Using synthetic models of the Galactic DWDs we demonstrate that the break in the $δa$ distribution corresponds to one side of a deep gap in the DWD separation distribution at around $a\approx 1$ au. Model DWDs with separations less than several au shrink dramatically due to (al least one) common envelope phase, reshaping the original separation distribution, clearing a gap and creating a pile-up of systems with $a\approx 0.01$ au and $δa < 0.01$. Our models reproduce the overall shape of the observed $δa$ distribution and its normalisation, however the predicted drop in the numbers of DWDs beyond the break is steeper than in the data.
△ Less
Submitted 7 March, 2022;
originally announced March 2022.