-
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 12 July, 2024; v1 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.
-
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.