-
On the Galactic rotation curve inferred from the Jeans equations Assessing its robustness using Gaia DR3 and cosmological simulations
Authors:
Orlin Koop,
Teresa Antoja,
Amina Helmi,
Thomas M. Callingham,
Chervin F. P. Laporte
Abstract:
Several works have recently applied Jeans modelling to Gaia-based datasets to infer the circular velocity curve for the Milky Way. Such works have consistently found evidence for a continuous decline in the rotation curve beyond $\sim$15kpc possibly indicative of a light dark matter halo. We used Gaia DR3 RVS data, supplemented with Bayesian distances to determine the radial variation of the secon…
▽ More
Several works have recently applied Jeans modelling to Gaia-based datasets to infer the circular velocity curve for the Milky Way. Such works have consistently found evidence for a continuous decline in the rotation curve beyond $\sim$15kpc possibly indicative of a light dark matter halo. We used Gaia DR3 RVS data, supplemented with Bayesian distances to determine the radial variation of the second moments of the velocity distribution for stars close to the Galactic plane. We have used these profiles to determine the rotation curve using the Jeans equations under the assumption of axisymmetry and explored how they vary with azimuth and above and below the Galactic disk plane. We have applied the same methodology to an N-body simulation of a Milky Way-like galaxy impacted by a satellite akin the Sagittarius dwarf and to the Auriga suite of cosmological simulations. We reveal evidence of disequilibrium and deviations from axisymmetry closer in. We find that the second moment of $V_R$ flattens out at $R \gtrsim 12.5$kpc, and that the second moment of $V_φ$ is different above and below the plane for $R \gtrsim 11$kpc. The simulations indicate that these features are typical of galaxies that have been perturbed by external satellites. They also suggest that the difference between the true circular velocity curve and that inferred from Jeans equations can be as high as 15$\%$, but is likely of order 10$\%$ for the Milky Way. This is of larger amplitude than the systematics inherent to Jeans equations. However, if the density of the tracer population were truncated at large radii, the erroneous conclusion of a steeply declining rotation curve can be reached. We find that steady-state axisymmetric Jeans modelling becomes less robust at large radii, indicating that particular caution is needed when interpreting the rotation curve inferred in those regions.
△ Less
Submitted 29 May, 2024;
originally announced May 2024.
-
Filling in the Blanks: A Method to Infer the Substructure Membership and Dynamics of 5D Stars
Authors:
Thomas M. Callingham,
Amina Helmi
Abstract:
We present and test a method to infer a probability density function (PDF) for the missing vlos of a star with 5D information within $2.5$ kpc. We use stars from the Gaia DR3 RVS catalogue to describe the local orbital structure in action space. This technique also allows us to infer the probability that a 5D star is associated with the Milky Way's stellar Disc or the stellar Halo, which can be fu…
▽ More
We present and test a method to infer a probability density function (PDF) for the missing vlos of a star with 5D information within $2.5$ kpc. We use stars from the Gaia DR3 RVS catalogue to describe the local orbital structure in action space. This technique also allows us to infer the probability that a 5D star is associated with the Milky Way's stellar Disc or the stellar Halo, which can be further decomposed into known stellar substructures. The method is tested on a 6D Gaia DR3 RVS sample and a 6D Gaia sample crossmatched to groundbased spectroscopic surveys, stripped of their true vlos. The stars predicted vlos, membership probabilities, and inferred structure properties are then compared to the true 6D equivalents, allowing the method's accuracy and limitations to be studied in detail. Our predicted vlos PDFs are statistically consistent with the true vlos, with accurate uncertainties. We find that the vlos of Disc stars can be well constrained, with a median uncertainty of 26 kms. Halo stars are typically less well constrained with a median uncertainty of 72 kms, but those found likely to belong to Halo substructures can be better constrained. The dynamical properties of the total sample and subgroups, such as distributions of integrals of motion and velocities, are also accurately recovered. The group membership probabilities are statistically consistent with our initial labelling, allowing high quality sets to be selected from 5D samples by choosing a trade off between higher expected purity and decreasing expected completeness.
△ Less
Submitted 16 May, 2024;
originally announced May 2024.
-
The 33 M$_\odot$ black hole Gaia BH3 is part of the disrupted ED-2 star cluster
Authors:
E. Balbinot,
E. Dodd,
T. Matsuno,
C. Lardo,
A. Helmi,
P. Panuzzo,
T. Mazeh,
B. Holl,
E. Caffau,
A. Jorissen,
C. Babusiaux,
P. Gavras,
L. Wyrzykowski,
L. Eyer,
N. Leclerc,
A. Bombrun,
N. Mowlavi,
G. M. Seabroke,
I. Cabrera-Ziri,
T. M. Callingham,
T. Ruiz-Lara,
E. Starkenburg
Abstract:
The Gaia Collaboration has recently reported the detection of a 33 M$_\odot$ black hole in a wide binary system located in the Solar neighbourhood. Here we explore the relationship between this black hole, known as Gaia BH3, and the nearby ED-2 halo stellar stream. We study the orbital characteristics of the Gaia BH3 binary and present measurements of the chemical abundances of ED-2 member stars d…
▽ More
The Gaia Collaboration has recently reported the detection of a 33 M$_\odot$ black hole in a wide binary system located in the Solar neighbourhood. Here we explore the relationship between this black hole, known as Gaia BH3, and the nearby ED-2 halo stellar stream. We study the orbital characteristics of the Gaia BH3 binary and present measurements of the chemical abundances of ED-2 member stars derived from high-resolution spectra obtained with the VLT. We find that the Galactic orbit of the Gaia BH3 system and its metallicity are entirely consistent with being part of the ED-2 stream. The characteristics of the stream, particularly its negligible spread in metallicity and in other chemical elements as well as its single stellar population, suggest that it originated from a disrupted star cluster of low mass. Its age is comparable to that of the globular cluster M92 that has been estimated to be as old as the Universe. This is the first black hole unambiguously associated with a disrupted star cluster. We infer a plausible mass range for the cluster to be relatively narrow, between $2\times 10^3M_\odot$ and $4.2\times 10^4M_\odot$. This implies that the black hole could have formed directly from the collapse of a massive very-metal-poor star, but that the alternative scenario of binary interactions inside the cluster environment also deserves to be explored.
△ Less
Submitted 23 April, 2024; v1 submitted 17 April, 2024;
originally announced April 2024.
-
ED-2: a cold but not so narrow stellar stream crossing the Solar neighbourhood
Authors:
E. Balbinot,
A. Helmi,
T. Callingham,
T. Matsuno,
E. Dodd,
T. Ruiz-Lara
Abstract:
ED-2 is a stellar stream identified as a compact group in integrals of motion space in a local sample of halo stars from the third Gaia data release. Here we investigate its nature and possible association with known halo substructures. We explore the current properties of ED-2 members in phase-space, and also analyse the expected distribution via orbit integration. In addition, we study the metal…
▽ More
ED-2 is a stellar stream identified as a compact group in integrals of motion space in a local sample of halo stars from the third Gaia data release. Here we investigate its nature and possible association with known halo substructures. We explore the current properties of ED-2 members in phase-space, and also analyse the expected distribution via orbit integration. In addition, we study the metallicity of ED-2 using APOGEE DR17 and LAMOST DR8 (and re-calibrated DR3). ED-2 forms a compact group in the $x-z$ (or $R-z$) plane, showing a pancake-like structure as it crosses the Solar neighbourhood. Dynamically it is most similar the globular clusters NGC 3201 and NGC 6101, and the stellar stream Ylgr and Phlegethon. However, its orbit is sufficiently different that none of these objects is likely to be ED-2's progenitor. We also find ED-2 to be quite metal-poor, with all of its stars $\mathrm{[Fe/H]} \leq -2.42$, with a median $\mathrm{[Fe/H]} = -2.60^{+0.20}_{-0.21}$. At this low metallicity, it is unlikely that ED-2 stems from any known globular cluster, instead, ED-2 seems to be in a similar category as the recently discovered Phoenix and C-19 stellar streams. We find that ED-2 members are scattered across the whole sky, which is due to its current orbital phase. We predict that as this object moves to its next apocentre it will acquire an on-sky morphology that is akin to cold stellar streams. Finally, since ED-2 is nearing pericentre, we predict that additional members found below the plane should have large radial velocities, close to $\sim$ 500 km/s in the present-day direction of the globular cluster NGC 6101.
△ Less
Submitted 12 September, 2023; v1 submitted 5 June, 2023;
originally announced June 2023.
-
Can we really pick and choose? Benchmarking various selections of Gaia Enceladus/Sausage stars in observations with simulations
Authors:
Andreia Carrillo,
Alis J. Deason,
Azadeh Fattahi,
Thomas M. Callingham,
Robert J. J. Grand
Abstract:
Large spectroscopic surveys plus Gaia astrometry have shown us that the inner stellar halo of the Galaxy is dominated by the debris of Gaia Enceladus/Sausage (GES). With the richness of data at hand, there are a myriad of ways these accreted stars have been selected. We investigate these GES selections and their effects on the inferred progenitor properties using data constructed from APOGEE and G…
▽ More
Large spectroscopic surveys plus Gaia astrometry have shown us that the inner stellar halo of the Galaxy is dominated by the debris of Gaia Enceladus/Sausage (GES). With the richness of data at hand, there are a myriad of ways these accreted stars have been selected. We investigate these GES selections and their effects on the inferred progenitor properties using data constructed from APOGEE and Gaia. We explore selections made in eccentricity, energy-angular momentum (E-Lz), radial action-angular momentum (Jr-Lz), action diamond, and [Mg/Mn]-[Al/Fe] in the observations, selecting between 144 and 1,279 GES stars with varying contamination from in-situ and other accreted stars. We also use the Auriga cosmological hydrodynamic simulations to benchmark the different GES dynamical selections. Applying the same observational GES cuts to nine Auriga galaxies with a GES, we find that the Jr-Lz method is best for sample purity and the eccentricity method for completeness. Given the average metallicity of GES (-1.28 < [Fe/H] < -1.18), we use the $z=0$ mass-metallicity relationship to find an average $\rm M_{\star}$ of $\sim 4 \times 10^{8}$ $\rm M_{\odot}$. We adopt a similar procedure and derive $\rm M_{\star}$ for the GES-like systems in Auriga and find that the eccentricity method overestimates the true $\rm M_{\star}$ by $\sim2.6\times$ while E-Lz underestimates by $\sim0.7\times$. Lastly, we estimate the total mass of GES to be $\rm 10^{10.5 - 11.1}~M_{\odot}$ using the relationship between the metallicity gradient and the GES-to-in-situ energy ratio. In the end, we cannot just `pick and choose' how we select GES stars, and instead should be motivated by the science question.
△ Less
Submitted 1 June, 2023;
originally announced June 2023.
-
The Gaia DR3 view of dynamical substructure in the stellar halo near the Sun
Authors:
Emma Dodd,
Thomas M. Callingham,
Amina Helmi,
Tadafumi Matsuno,
Tomás Ruiz-Lara,
Eduardo Balbinot,
Sofie Lövdal
Abstract:
The debris from past merger events is expected and, to some extent, known to populate the stellar halo near the Sun. We aim to identify and characterise such merger debris using Gaia DR3 data supplemented by metallicity and chemical abundance information from LAMOST LRS and APOGEE for halo stars within 2.5 kpc from the Sun. We utilise a single linkage-based clustering algorithm to identify over-de…
▽ More
The debris from past merger events is expected and, to some extent, known to populate the stellar halo near the Sun. We aim to identify and characterise such merger debris using Gaia DR3 data supplemented by metallicity and chemical abundance information from LAMOST LRS and APOGEE for halo stars within 2.5 kpc from the Sun. We utilise a single linkage-based clustering algorithm to identify over-densities in Integrals of Motion space that could be due to merger debris. Combined with metallicity information and chemical abundances, we characterise these statistically significant over-densities. We find that the local stellar halo contains 7 main dynamical groups, some of in-situ and some of accreted origin, most of which are already known. We report the discovery of a new substructure, which we name ED-1. In addition, we find evidence for 11 independent smaller clumps, 5 of which are new: ED-2, 3, 4, 5 and 6 are typically rather tight dynamically, depict a small range of metallicities, and their abundances when available, as well as their location in Integrals of Motion space, suggest an accreted origin. The local halo contains an important amount of substructure, of both in-situ and accreted origin.
△ Less
Submitted 19 July, 2022; v1 submitted 22 June, 2022;
originally announced June 2022.
-
The Chemo-Dynamical Groups of Galactic Globular Clusters
Authors:
Thomas M. Callingham,
Marius Cautun,
Alis J. Deason,
Carlos S. Frenk,
Robert J. J Grand,
Federico Marinacci
Abstract:
We introduce a multi-component chemo-dynamical method for splitting the Galactic population of Globular Clusters (GCs) into three distinct constituents: bulge, disc, and stellar halo. The latter is further decomposed into the individual large accretion events that built up the Galactic stellar halo: the Gaia-Enceladus-Sausage, Kraken and Sequoia structures, and the Sagittarius and Helmi streams. O…
▽ More
We introduce a multi-component chemo-dynamical method for splitting the Galactic population of Globular Clusters (GCs) into three distinct constituents: bulge, disc, and stellar halo. The latter is further decomposed into the individual large accretion events that built up the Galactic stellar halo: the Gaia-Enceladus-Sausage, Kraken and Sequoia structures, and the Sagittarius and Helmi streams. Our modelling is extensively tested using mock GC samples constructed from the AURIGA suite of hydrodynamical simulations of Milky Way (MW)-like galaxies. We find that, on average, a proportion of the accreted GCs cannot be associated with their true infall group and are left ungrouped, biasing our recovered population numbers to approximately 80 percent of their true value. Furthermore, the identified groups have a completeness and a purity of only 65 percent. This reflects the difficulty of the problem, a result of the large degree of overlap in energy-action space of the debris from past accretion events. We apply the method to the Galactic data to infer, in a statistically robust and easily quantifiable way, the GCs associated with each MW accretion event. The resulting groups' population numbers of GCs, corrected for biases, are then used to infer the halo and stellar masses of the now defunct satellites that built up the halo of the MW.
△ Less
Submitted 9 August, 2022; v1 submitted 1 February, 2022;
originally announced February 2022.
-
The orbital phase space of contracted dark matter halos
Authors:
Thomas M. Callingham,
Marius Cautun,
Alis J. Deason,
Carlos S. Frenk,
Robert J. J. Grand,
Federico Marinacci,
Ruediger Pakmor
Abstract:
We study the orbital phase-space of dark matter (DM) halos in the AURIGA suite of cosmological hydrodynamics simulations of Milky Way analogues. We characterise halos by their spherical action distribution, a function of the specific angular momentum, and the radial action, of the DM particles. By comparing DM-only and hydrodynamical simulations of the same halos, we investigate the contraction of…
▽ More
We study the orbital phase-space of dark matter (DM) halos in the AURIGA suite of cosmological hydrodynamics simulations of Milky Way analogues. We characterise halos by their spherical action distribution, a function of the specific angular momentum, and the radial action, of the DM particles. By comparing DM-only and hydrodynamical simulations of the same halos, we investigate the contraction of DM halos caused by the accumulation of baryons at the centre. We find a small systematic suppression of the radial action in the DM halos of the hydrodynamical simulations, suggesting that the commonly used adiabatic contraction approximation can result in an underestimate of the density by ~ 8%. We apply an iterative algorithm to contract the AURIGA DM halos given a baryon density profile and halo mass, recovering the true contracted DM profiles with an accuracy of ~15%, that reflects halo-to-halo variation. Using this algorithm, we infer the total mass profile of the Milky Way's contracted DM halo. We derive updated values for the key astrophysical inputs to DM direct detection experiments: the DM density and velocity distribution in the Solar neighbourhood.
△ Less
Submitted 21 January, 2020;
originally announced January 2020.
-
The mass of the Milky Way from satellite dynamics
Authors:
Thomas Callingham,
Marius Cautun,
Alis J. Deason,
Carlos S. Frenk,
Wenting Wang,
Facundo A. Gómez,
Robert J. J. Grand,
Federico Marinacci,
Rüdgier Pakmor
Abstract:
We present and apply a method to infer the mass of the Milky Way (MW) by comparing the dynamics of MW satellites to those of model satellites in the EAGLE cosmological hydrodynamics simulations. A distribution function (DF) for galactic satellites is constructed from EAGLE using specific angular momentum and specific energy, which are scaled so as to be independent of host halo mass. In this 2-dim…
▽ More
We present and apply a method to infer the mass of the Milky Way (MW) by comparing the dynamics of MW satellites to those of model satellites in the EAGLE cosmological hydrodynamics simulations. A distribution function (DF) for galactic satellites is constructed from EAGLE using specific angular momentum and specific energy, which are scaled so as to be independent of host halo mass. In this 2-dimensional space, the orbital properties of satellite galaxies vary according to the host halo mass. The halo mass can be inferred by calculating the likelihood that the observed satellite population is drawn from this DF. Our method is robustly calibrated on mock EAGLE systems. We validate it by applying it to the completely independent suite of 30 AURIGA high-resolution simulations of MW-like galaxies: the method accurately recovers their true mass and associated uncertainties. We then apply it to ten classical satellites of the MW with 6D phase-space measurements, including updated proper motions from the GAIA satellite. The mass of the MW is estimated to be $M_{200}^{\textnormal{MW}}=1.17_{-0.15}^{+0.21}\times10^{12}M_{\odot}$ (68\% confidence limits). We combine our total mass estimate with recent mass estimates in the inner regions of the Galaxy to infer an inner dark matter (DM) mass fraction $M^\textnormal{DM}(<20~\rm{kpc})/M^\textnormal{DM}_{200}=0.12$ which is typical of ${\sim}10^{12}M_{\odot}$ $Λ$CDM haloes in hydrodynamical galaxy formation simulations. Assuming an NFW profile, this is equivalent to a halo concentration of $c_{200}^{\textnormal{MW}}=10.9^{+2.6}_{-2.0}$.
△ Less
Submitted 23 January, 2019; v1 submitted 30 August, 2018;
originally announced August 2018.