-
Hidden deep in the halo: Selection of a reduced proper motion halo catalogue and mining retrograde streams in the velocity space
Authors:
Akshara Viswanathan,
Else Starkenburg,
Helmer H. Koppelman,
Amina Helmi,
Eduardo Balbinot,
Anna F. Esselink
Abstract:
The Milky Way halo is one of the few galactic haloes that provides a unique insight into galaxy formation by resolved stellar populations. Here, we present a catalogue of $\sim$47 million halo stars selected independent of parallax and line-of-sight velocities, using a combination of Gaia DR3 proper motion and photometry by means of their reduced proper motion. We select high tangential velocity (…
▽ More
The Milky Way halo is one of the few galactic haloes that provides a unique insight into galaxy formation by resolved stellar populations. Here, we present a catalogue of $\sim$47 million halo stars selected independent of parallax and line-of-sight velocities, using a combination of Gaia DR3 proper motion and photometry by means of their reduced proper motion. We select high tangential velocity (halo) main sequence stars and fit distances to them using their simple colour-absolute-magnitude relation. This sample reaches out to $\sim$21 kpc with a median distance of $6.6$ kpc thereby probing much further out than would be possible using reliable Gaia parallaxes. The typical uncertainty in their distances is $0.57_{-0.26}^{+0.56}$ kpc. Using the colour range $0.45<(G_0-G_\mathrm{RP,0})<0.715$ where the main sequence is narrower, gives an even better accuracy down to $0.39_{-0.12}^{+0.18}$ kpc in distance. The median velocity uncertainty for stars within this colour range is 15.5 km/s. The distribution of these sources in the sky, together with their tangential component velocities, are very well-suited to study retrograde substructures. We explore the selection of two complex retrograde streams: GD-1 and Jhelum. For these streams, we resolve the gaps, wiggles and density breaks reported in the literature more clearly. We also illustrate the effect of the kinematic selection bias towards high proper motion stars and incompleteness at larger distances due to Gaia's scanning law. These examples showcase how the full RPM catalogue made available here can help us paint a more detailed picture of the build-up of the Milky Way halo.
△ Less
Submitted 31 January, 2023;
originally announced February 2023.
-
High-precision chemical abundances of Galactic building blocks. II. Revisiting the chemical distinctness of the Helmi streams
Authors:
Tadafumi Matsuno,
Emma Dodd,
Helmer H. Koppelman,
Amina Helmi,
Miho N. Ishigaki,
Wako Aoki,
Jingkun Zhao,
Zhen Yuan,
Kohei Hattori
Abstract:
Context: The Helmi streams are a kinematic substructure whose progenitor is likely a dwarf galaxy. Although 20 years have passed since their discovery, it is still unclear whether their members are chemically distinguishable from other halo stars in the Milky Way. Aim: We aim to precisely characterize the chemical properties of the Helmi streams. Methods: We analyzed high-resolution, high signal-t…
▽ More
Context: The Helmi streams are a kinematic substructure whose progenitor is likely a dwarf galaxy. Although 20 years have passed since their discovery, it is still unclear whether their members are chemically distinguishable from other halo stars in the Milky Way. Aim: We aim to precisely characterize the chemical properties of the Helmi streams. Methods: We analyzed high-resolution, high signal-to-noise ratio spectra for 11 Helmi stream stars through a line-by-line abundance analysis. We compared the derived abundances to homogenized literature abundances of the other halo stars, including those belonging to other kinematic substructures, such as Gaia-Enceladus and Sequoia. Results: Compared to typical halo stars, the Helmi stream members clearly show low values of [X/Fe] in elements produced by massive stars, such as Na and $α$-elements. This tendency is seen down to metallicities of at least [Fe/H]$\sim -2.2$, suggesting type~Ia supernovae already started to contribute to the chemical evolution at this metallicity. We find that the [$α$/Fe] ratio does not evolve significantly with metallicity, making the Helmi stream stars less distinguishable from Gaia-Enceladus stars at [Fe/H]$\gtrsim -1.5$. The almost constant but low value of [$α$/Fe] might be indicative of quiescent star formation with low efficiency at the beginning and bursty star formation at later times. We also find extremely low values of [Y/Fe] at low metallicity, providing further support for the claim that light neutron-capture elements are deficient in Helmi streams. While Zn is deficient at low metallicity, it shows a large spread at high metallicity. The origin of the extremely low Y abundances and Zn variations remains unclear. Conclusion: The Helmi stream stars are distinguishable from the majority of the halo stars if homogeneously derived abundances are compared.
△ Less
Submitted 6 July, 2023; v1 submitted 22 March, 2022;
originally announced March 2022.
-
Substructure in the stellar halo near the Sun. II. Characterisation of independent structures
Authors:
Tomás Ruiz-Lara,
Tadafumi Matsuno,
S. Sofie Lövdal,
Amina Helmi,
Emma Dodd,
Helmer H. Koppelman
Abstract:
In Lövdal et al, we presented a data-driven method for clustering in Integrals of Motion space and applied it to a large sample of nearby halo stars with 6D phase-space information. We identified a large number of clusters, many of which could tentatively be merged into larger groups. Our goal is to establish the reality of the clusters through a combined study of their stellar populations to gain…
▽ More
In Lövdal et al, we presented a data-driven method for clustering in Integrals of Motion space and applied it to a large sample of nearby halo stars with 6D phase-space information. We identified a large number of clusters, many of which could tentatively be merged into larger groups. Our goal is to establish the reality of the clusters through a combined study of their stellar populations to gain more insights into the accretion history of the Milky Way. We develop a procedure that quantifies the similarity of clusters based on KS tests using their metallicity distribution functions, and an isochrone fitting method to determine their average age, which is also used to compare the distribution of stars in the Colour-Absolute magnitude diagram. This allows us to group clusters into substructures, and to compare substructures with one another. The clusters identified are merged into 12 extended substructures, while 8 small clusters remain as such. The large substructures include the previously known Gaia-Enceladus, Helmi streams, Sequoia, and Thamnos 1 and 2. We identify overdensities associated with the hot thick disc and hosting a metal-poor population. Especially notable is our largest substructure which, although peaking at the metallicity characteristic of the thick disk has a well populated metal-poor component, and dynamics in-between hot thick disc and halo. We identify additional debris in the region occupied by Sequoia with distinct kinematics, likely remnants of three different accretion events with progenitors of similar mass. We also identify different trends of [Mg/Fe] vs [Fe/H] for the various substructures confirming our dissection of the nearby halo. At least 20\% of the halo near the Sun is associated to substructures. When comparing their global properties, we note that those substructures on retrograde orbits are not only more metal-poor on average but also older.
△ Less
Submitted 28 April, 2022; v1 submitted 7 January, 2022;
originally announced January 2022.
-
Substructure in the stellar halo near the Sun. I. Data-driven clustering in Integrals of Motion space
Authors:
S. Sofie Lövdal,
Tomás Ruiz-Lara,
Helmer H. Koppelman,
Tadafumi Matsuno,
Emma Dodd,
Amina Helmi
Abstract:
Aims: Develop a data-driven and statistically based method for finding such clumps in Integrals of Motion space for nearby halo stars and evaluating their significance robustly. Methods: We use data from Gaia EDR3 extended with radial velocities from ground-based spectroscopic surveys to construct a sample of halo stars within 2.5 kpc from the Sun. We apply a hierarchical clustering method that us…
▽ More
Aims: Develop a data-driven and statistically based method for finding such clumps in Integrals of Motion space for nearby halo stars and evaluating their significance robustly. Methods: We use data from Gaia EDR3 extended with radial velocities from ground-based spectroscopic surveys to construct a sample of halo stars within 2.5 kpc from the Sun. We apply a hierarchical clustering method that uses the single linkage algorithm in a 3D space defined by the commonly used integrals of motion energy $E$, together with two components of the angular momentum, $L_z$ and $L_\perp$. To evaluate the statistical significance of the clusters found, we compare the density within an ellipsoidal region centered on the cluster to that of random sets with similar global dynamical properties. We pick out the signal at the location of their maximum statistical significance in the hierarchical tree. We estimate the proximity of a star to the cluster center using the Mahalanobis distance. We also apply the HDBSCAN clustering algorithm in velocity space. Results: Our procedure identifies 67 highly significant clusters ($ > 3σ$), containing 12\% of the sources in our halo set, and in total 232 subgroups or individual streams in velocity space. In total, 13.8\% of the stars in our data set can be confidently associated to a significant cluster based on their Mahalanobis distance. Inspection of our data set reveals a complex web of relationships between the significant clusters, suggesting that they can be tentatively grouped into at least 6 main structures, many of which can be associated to previously identified halo substructures, and a number of independent substructures. This preliminary conclusion is further explored in an accompanying paper by Ruiz-Lara et al., where we also characterize the substructures in terms of their stellar populations. Conclusions: We find... (abridged version)
△ Less
Submitted 17 May, 2022; v1 submitted 7 January, 2022;
originally announced January 2022.
-
High-precision chemical abundances of Galactic building blocks. The distinct chemical abundance sequence of Sequoia
Authors:
Tadafumi Matsuno,
Helmer H. Koppelman,
Amina Helmi,
Wako Aoki,
Miho N. Ishigaki,
Takuma Suda,
Zhen Yuan,
Kohei Hattori
Abstract:
Context: Sequoia is a retrograde kinematic substructure in the nearby Galactic halo, whose properties are a matter of debate. For example, previous studies do not necessarily agree on the chemical abundances of Sequoia stars, which are important for understanding its nature. Aim: We characterize the chemical properties of a sample of stars from Sequoia by determining high-precision abundances. Met…
▽ More
Context: Sequoia is a retrograde kinematic substructure in the nearby Galactic halo, whose properties are a matter of debate. For example, previous studies do not necessarily agree on the chemical abundances of Sequoia stars, which are important for understanding its nature. Aim: We characterize the chemical properties of a sample of stars from Sequoia by determining high-precision abundances. Methods: We measured abundances of Na, Mg, Si, Ca, Ti, Cr, Mn, Ni, Zn, Y, and Ba from a differential abundance analysis on high signal-to-noise ratio, high-resolution spectra from new observations and from archival data. We compared precisely measured chemical abundances of 12 Sequoia candidates with those of typical halo stars from the literature, which also includes stars from Gaia-Enceladus. This allowed us to characterize Sequoia and compare it to another Galactic building block. The comparison was made after putting all the abundances onto the same scale using standard stars. Results: There are significant differences in [Na/Fe], [Mg/Fe], [Ca/Fe], [Ti/Fe], [Zn/Fe], and [Y/Fe] between Sequoia and Gaia-Enceladus stars at $-1.8\lesssim [\mathrm{Fe/H}]\lesssim -1.4$ in the sense that these abundance ratios are lower in Sequoia. These differences are similar to those seen between Gaia-Enceladus and in situ stars at a higher metallicity, suggesting that Sequoia is affected by type~Ia supernovae at a lower metallicity than Gaia-Enceladus. We also confirm that the low [{Mg}/{Fe}] of Sequoia is seen in the literature and in surveys, namely APOGEE DR17 and GALAH DR3, if the stars are kinematically selected in the same way. Conclusion: Sequoia stars have a distinct chemical abundance pattern and can be chemically separated from in situ stars or Gaia-Enceladus stars if abundances are measured with sufficient precision, namely $σ([\mathrm{X/Fe}])\lesssim 0.07\,\mathrm{dex}$.
△ Less
Submitted 6 July, 2023; v1 submitted 30 November, 2021;
originally announced November 2021.
-
Age Determination of Galaxy Merger Remnant Stars using Asteroseismology
Authors:
Camilla C. Borre,
Víctor Aguirre Børsen-Koch,
Amina Helmi,
Helmer H. Koppelman,
Martin B. Nielsen,
Jakob L. Rørsted,
Dennis Stello,
Amalie Stokholm,
Mark L. Winther,
Guy R. Davies,
Marc Hon,
J. M. Diederik Kruijssen,
Chervin Laporte,
Claudia Reyes,
Jie Yu
Abstract:
The Milky Way was shaped by the mergers with several galaxies in the past. We search for remnant stars that were born in these foreign galaxies and assess their ages in an effort to put upper limits on the merger times and thereby better understand the evolutionary history of our Galaxy. Using 6D-phase space information from Gaia eDR3 and chemical information from APOGEE DR16, we kinematically and…
▽ More
The Milky Way was shaped by the mergers with several galaxies in the past. We search for remnant stars that were born in these foreign galaxies and assess their ages in an effort to put upper limits on the merger times and thereby better understand the evolutionary history of our Galaxy. Using 6D-phase space information from Gaia eDR3 and chemical information from APOGEE DR16, we kinematically and chemically select $23$ red giant stars belonging to former dwarf galaxies that merged with the Milky Way. With added asteroseismology from Kepler and K2, we determine the ages of the $23$ ex-situ stars and $55$ in-situ stars with great precision. We find that all the ex-situ stars are consistent with being older than $8$ Gyr. While it is not possible to associate all the stars with a specific dwarf galaxy we classify eight of them as Gaia-Enceladus/Sausage stars, which is one of the most massive mergers in our Galaxy's history. We determine their mean age to be $9.5^{+1.2}_{-1.3}$ Gyr consistent with a merger time of $8$-$10$ Gyr ago. The rest of the stars are possibly associated with Kraken, Thamnos, Sequoia, or another extragalactic progenitor. The age determination of ex-situ stars paves the way to more accurately pinning down when the merger events occurred and hence provide tight constraints useful for simulating how these events unfolded.
△ Less
Submitted 2 November, 2021;
originally announced November 2021.
-
Substructures, Resonances and debris streams
Authors:
Emma Dodd,
Amina Helmi,
Helmer H. Koppelman
Abstract:
The local stellar halo of the Milky Way contains the debris from several past accretion events. Here we study in detail the structure and properties of nearby debris associated with the Helmi streams, originally identified as an overdensity in integrals of motion space. We use 6D phase-space information from Gaia EDR3 combined with spectroscopic surveys, and we analyse the orbits and frequencies o…
▽ More
The local stellar halo of the Milky Way contains the debris from several past accretion events. Here we study in detail the structure and properties of nearby debris associated with the Helmi streams, originally identified as an overdensity in integrals of motion space. We use 6D phase-space information from Gaia EDR3 combined with spectroscopic surveys, and we analyse the orbits and frequencies of the stars in the streams using various Galactic potentials. We also explore how the Helmi streams constrain the flattening q, of the Galactic dark matter halo. We find that the streams are split into substructures in integrals of motion space, most notably into two clumps in angular momentum space. The clumps have consistent metallicity distributions and stellar populations, supporting a common progeny. In all the realistic Galactic potentials explored, the Helmi streams stars depict a diffuse distribution close to Omega_z /Omega_R ~0.7. At the same time, the reason for the substructure in angular momentum space appears to be a Omega_z : Omega_phi resonance close to the 1:1. This resonance is exactly the 1:1 in the case that the (density) flattening of the dark halo is q = 1.2. For this halo shape the substructure in angular momenta is also long lasting. Our findings suggest that the structure of the Galactic potential leaves a clear imprint on the properties of phase-mixed debris streams.
△ Less
Submitted 27 September, 2021; v1 submitted 20 May, 2021;
originally announced May 2021.
-
Time evolution of gaps in stellar streams in axisymmetric Stäckel potentials
Authors:
Helmer H. Koppelman,
Amina Helmi
Abstract:
When a subhalo interacts with a cold stellar stream it perturbs its otherwise nearly smooth distribution of stars, and this leads to the creation of a gap. The properties of such gaps depend on the parameters of the interaction. Their characterisation could thus lead to the determination of the mass spectrum of the perturbers and possibly reveal the existence of dark subhalos orbiting the Milky Wa…
▽ More
When a subhalo interacts with a cold stellar stream it perturbs its otherwise nearly smooth distribution of stars, and this leads to the creation of a gap. The properties of such gaps depend on the parameters of the interaction. Their characterisation could thus lead to the determination of the mass spectrum of the perturbers and possibly reveal the existence of dark subhalos orbiting the Milky Way. Our goal is to construct a fully analytical model of the formation and evolution of gaps embedded in streams orbiting in a realistic Milky Way potential. To this end, we extend the model of Helmi & Koppelman (2016) for spherical potentials, and predict the properties of gaps in streams evolving in axisymmetric Stäckel potentials. We make use of action-angles and their simple behaviour to calculate the divergence of initially nearby orbits slightly perturbed by the interaction with a subhalo. Our model, corroborated by N-body experiments, predicts that the size of a gap grows linearly with time. We obtain analytical expressions for the dependencies of the growth rate on the orbit of the stream, the properties of the subhalo (mass, scale-radius), and the geometry of the encounter (relative velocity, impact parameter). We find that the density at the centre of the gap decreases with time as a power-law in the same way as the density of a stream. This results in the density-contrast between a pristine and a perturbed stream on the same orbit asymptotically reaching a constant value dependent only on the encounter parameters. We find that at a fixed age, smallish gaps are sensitive mostly to the mass of the subhalo, while gaps formed by subhalo flybys with a low relative velocity, or when the stream and subhalo move parallel, are degenerate to the encounter parameters.
△ Less
Submitted 23 November, 2020;
originally announced November 2020.
-
Frequencies, chaos and resonances: a study of orbital parameters of nearby thick disc and halo stars
Authors:
Helmer H. Koppelman,
Jorrit H. J. Hagen,
Amina Helmi
Abstract:
We study the distribution of nearby thick disc and halo stars in subspaces defined by their characteristic orbital parameters. Our aim is to establish the origin of the structure reported in particular, in the $R_{\rm max}-z_{\rm max}$ space. To this end, we compute the orbital parameters and frequencies of stars for a generic and for a Stäckel Milky Way potential. We find that for both the thick…
▽ More
We study the distribution of nearby thick disc and halo stars in subspaces defined by their characteristic orbital parameters. Our aim is to establish the origin of the structure reported in particular, in the $R_{\rm max}-z_{\rm max}$ space. To this end, we compute the orbital parameters and frequencies of stars for a generic and for a Stäckel Milky Way potential. We find that for both the thick disc and halo populations very similar prominent substructures are apparent for the generic Galactic potential, while no substructure is seen for the Stäckel model. This indicates that the origin of these features is not merger-related, but due to non-integrability of the generic potential. This conclusion is strengthened by our frequency analysis of the orbits of stars, which reveals the presence of prominent resonances, with $\sim 30\%$ of the halo stars associated to resonance families. In fact, the stars in resonances define the substructures seen in the spaces of characteristic orbital parameters. Furthermore, we find that some stars in our sample and in debris streams are on the same resonance as the Sagittarius dwarf, suggesting this system has also influenced the distribution of stars in the Galactic thick disc and halo components. Our study constitutes a step towards disentangling the imprint of merger debris from substructures driven by internal dynamics. Given their prominence, these resonant-driven overdensities could potentially be useful to constrain the exact form of the Galactic potential.
△ Less
Submitted 10 September, 2020;
originally announced September 2020.
-
Galactic potential constraints from clustering in action space of combined stellar stream data
Authors:
Stella Reino,
Elena M. Rossi,
Robyn E. Sanderson,
Elena Sellentin,
Amina Helmi,
Helmer H. Koppelman,
Sanjib Sharma
Abstract:
Stream stars removed by tides from their progenitor satellite galaxy or globular cluster act as a group of test particles on neighboring orbits, probing the gravitational field of the Milky Way. While constraints from individual streams have been shown to be susceptible to biases, combining several streams from orbits with various distances reduces these biases. We fit a common gravitational poten…
▽ More
Stream stars removed by tides from their progenitor satellite galaxy or globular cluster act as a group of test particles on neighboring orbits, probing the gravitational field of the Milky Way. While constraints from individual streams have been shown to be susceptible to biases, combining several streams from orbits with various distances reduces these biases. We fit a common gravitational potential to multiple stellar streams simultaneously by maximizing the clustering of the stream stars in action space. We apply this technique to members of the GD-1, Pal 5, Orphan and Helmi streams, exploiting both the individual and combined data sets. We describe the Galactic potential with a Stäckel model, and vary up to five parameters simultaneously. We find that we can only constrain the enclosed mass, and that the strongest constraints come from the GD-1, Pal 5 and Orphan streams whose combined data set yields $M(< 20\ \mathrm{kpc}) = 2.96^{+0.25}_{-0.26} \times 10^{11} \ M_{\odot}$. When including the Helmi stream in the data set, the mass uncertainty increases to $M(< 20\ \mathrm{kpc}) = 3.12^{+3.21}_{-0.46} \times 10^{11} \ M_{\odot}$.
△ Less
Submitted 2 February, 2021; v1 submitted 1 July, 2020;
originally announced July 2020.
-
Determination of the escape velocity of the Milky Way using a proper motion selected halo sample
Authors:
Helmer H. Koppelman,
Amina Helmi
Abstract:
The {\it Gaia} mission has provided the largest catalogue ever of sources with tangential velocity information. However, using this catalogue for dynamical studies is difficult because most of the stars lack line-of-sight velocity measurements. Recently, we presented a selection of $\sim 10^7$ halo stars with accurate distances that have been selected based on their photometry and proper motions.…
▽ More
The {\it Gaia} mission has provided the largest catalogue ever of sources with tangential velocity information. However, using this catalogue for dynamical studies is difficult because most of the stars lack line-of-sight velocity measurements. Recently, we presented a selection of $\sim 10^7$ halo stars with accurate distances that have been selected based on their photometry and proper motions. Using this sample, we model the tail of the velocity distribution with a power-law distribution, a commonly used approach first established by \cite{Leonard1990THESPEED}. For the first time ever we use tangential velocities measured accurately for an unprecedented number of halo stars to estimate the escape velocity. In the solar neighbourhood, we obtain a very precise estimate of the escape velocity which is $497^{+8}_{-8}~{\rm km/s}$. This estimate is most likely biased low, our best guess is by 10\%. As a result, the true escape velocity most likely is closer to $550~{\rm km/s}$. The escape velocity directly constrains the total mass of the Milky Way. To find the best fitting halo mass and concentration parameter we adjusted the dark (spherical NFW) halo of a realistic Milky Way potential while keeping the circular velocity at the solar radius fixed at $v_c(R_\odot) = 232.8~{\rm km/s}$. The resulting halo parameters are $M_{200}^{+10\%} = 1.11^{+0.08}_{-0.07} \cdot10^{12} ~{\rm M}_\odot$ and concentration parameter $c^{+10\%} = 11.8^{+0.3}_{-0.3}$, where we use the explicit notation to indicate that these have been corrected for the 10\% bias. The slope of the escape velocity with galactocentric distance is as expected in the inner Galaxy based on Milky Way models. Curiously, we find a disagreement beyond the solar radius which is likely an effect of a change in the shape of the velocity distribution and could be related to the presence of velocity clumps.
△ Less
Submitted 8 March, 2021; v1 submitted 29 June, 2020;
originally announced June 2020.
-
A massive mess: When a large dwarf and a Milky Way-like galaxy merge
Authors:
Helmer H. Koppelman,
Roy O. Y. Bos,
Amina Helmi
Abstract:
Circa 10 billion years ago the Milky Way merged with a massive satellite, Gaia-Enceladus. To gain insight into the properties of its debris we analyse in detail the suite of simulations from Villalobos & Helmi (2008), which includes an experiment that produces a good match to the kinematics of nearby halo stars inferred from Gaia data. We compare the kinematic distributions of stellar particles in…
▽ More
Circa 10 billion years ago the Milky Way merged with a massive satellite, Gaia-Enceladus. To gain insight into the properties of its debris we analyse in detail the suite of simulations from Villalobos & Helmi (2008), which includes an experiment that produces a good match to the kinematics of nearby halo stars inferred from Gaia data. We compare the kinematic distributions of stellar particles in the simulations and study the distribution of debris in orbital angular momentum, eccentricity and energy, and its relation to the mass-loss history of the simulated satellite. We confirm that Gaia-Enceladus probably fell in on a retrograde, 30$^\circ$ inclination orbit. We find that while 75% of the debris in our preferred simulation has large eccentricity ($> 0.8$), roughly 9% has eccentricity smaller than 0.6. Star particles lost early have large retrograde motions, and a subset of these have low eccentricity. Such stars would be expected to have lower metallicities as they stem from the outskirts of the satellite, and hence naively they could be confused with debris associated with a separate system. These considerations seem to apply to some of the stars from the postulated Sequoia galaxy. When a massive discy galaxy merges, it leaves behind debris with a complex phase-space structure, a large range of orbital properties, and a range of chemical abundances. Observationally, this results in substructures with very different properties, which can be misinterpreted as implying independent progeny. Detailed chemical abundances of large samples of stars and tailored hydrodynamical simulations are critical to resolving such conundrums.
△ Less
Submitted 13 June, 2020;
originally announced June 2020.
-
The Reduced Proper Motion selected halo: methods and description of the catalogue
Authors:
Helmer H. Koppelman,
Amina Helmi
Abstract:
The Gaia mission has provided the largest ever astrometric chart of the Milky Way. Using it to map the Galactic halo is helpful for disentangling its merger history. The identification of halo stars in Gaia DR2 with reliable distance estimates requires special methods because such stars are typically farther away and scarce. We apply the reduced proper motion (RPM) method to identify halo main seq…
▽ More
The Gaia mission has provided the largest ever astrometric chart of the Milky Way. Using it to map the Galactic halo is helpful for disentangling its merger history. The identification of halo stars in Gaia DR2 with reliable distance estimates requires special methods because such stars are typically farther away and scarce. We apply the reduced proper motion (RPM) method to identify halo main sequence stars on the basis of Gaia photometry and proper motions. Using the colour-absolute-magnitude relation for this type of stars, we calculate photometric distances. Our selection results in a set of $\sim10^7$ tentative main sequence halo stars with typical distance uncertainties of $7\%$ and with median velocity errors of 20 km/s. The median distance of our sample is $\sim 4.4$ kpc, with the faintest stars located at $\sim 16$ kpc. The spatial distribution of the stars in our sample is centrally concentrated. Visual inspection of the mean velocities of stars on the sky reveals large-scale patterns as well as clear imprints of the GD-1 stream and tentative hints of the Jhelum and Leiptr streams. Incompleteness and selection effects limit our ability to interpret the patterns reliably as well as to identify new substructures. We define a pseudo-velocity space by setting to zero the line-of-sight velocities of our sample stars. In this space, we recover several known structures such as the footprint of Gaia-Enceladus (i.e. the Gaia-Sausage) as well as the Helmi streams and some other retrograde substructures (Sequoia, Thamnos). We show that the two-point velocity correlation function reveals significant clustering on scales smaller than 100 km/s, of similar amplitude as found for the 6D Gaia halo sample. This clustering indicates the presence of nearby streams that are predominantly phase-mixed.
△ Less
Submitted 24 November, 2021; v1 submitted 15 April, 2020;
originally announced April 2020.
-
Multiple retrograde substructures in the Galactic halo: A shattered view of Galactic history
Authors:
Helmer H. Koppelman,
Amina Helmi,
Davide Massari,
Adrian M. Price-Whelan,
Tjitske K. Starkenburg
Abstract:
Aims. Several kinematic and chemical substructures have been recently found amongst Milky Way halo stars with retrograde motions. It is currently unclear how these various structures are related to each other. This Letter aims to shed light on this issue. Methods. We explore the retrograde halo with an augmented version of the Gaia DR2 RVS sample, extended with data from three large spectroscopic…
▽ More
Aims. Several kinematic and chemical substructures have been recently found amongst Milky Way halo stars with retrograde motions. It is currently unclear how these various structures are related to each other. This Letter aims to shed light on this issue. Methods. We explore the retrograde halo with an augmented version of the Gaia DR2 RVS sample, extended with data from three large spectroscopic surveys, namely RAVE, APOGEE and LAMOST. In this dataset, we identify several structures using the HDBSCAN clustering algorithm. We discuss their properties and possible links using all the available chemical and dynamical information. Results. In concordance with previous work, we find that stars with [Fe/H] $<-1$ have more retrograde motions than those with [Fe/H] $>-1$. The retrograde halo contains a mixture of debris from objects like Gaia-Enceladus, Sequoia, and even the chemically defined thick-disc. We find that the Sequoia has a smaller range in orbital energies than previously suggested and is confined to high-energy. Sequoia could be a small galaxy in itself, but since it overlaps both in integrals-of-motion space and chemical abundance space with the less bound debris of Gaia-Enceladus, its nature cannot be fully settled yet. In the low-energy part of the halo we find evidence for at least one more distinct structure: Thamnos. Stars in Thamnos are on low inclination, mildly eccentric retrograde orbits, moving at $v_φ\approx-150$ km/s, and are chemically distinct from the other structures. Conclusions. Even with the excellent Gaia DR2 data it remains challenging to piece together all the fragments found in the retrograde halo. At this point, we are very much in need of large datasets with high-quality high-resolution spectra and tailored high-resolution hydrodynamical simulations of galaxy mergers.
△ Less
Submitted 19 September, 2019;
originally announced September 2019.
-
Origin of the system of globular clusters in the Milky Way
Authors:
Davide Massari,
Helmer H. Koppelman,
Amina Helmi
Abstract:
$Context$. The assembly history experienced by the Milky Way is currently being unveiled thanks to the data provided by the $Gaia$ mission. It is likely that the globular cluster system of our Galaxy has followed a similarly intricate formation path. $Aims…
▽ More
$Context$. The assembly history experienced by the Milky Way is currently being unveiled thanks to the data provided by the $Gaia$ mission. It is likely that the globular cluster system of our Galaxy has followed a similarly intricate formation path. $Aims$. To constrain this formation path, we explore the link between the globular clusters and the known merging events that the Milky Way has experienced. $Methods$. To this end, we combined the kinematic information provided by $Gaia$ for almost all Galactic clusters, with the largest sample of cluster ages available after carefully correcting for systematic errors. To identify clusters with a common origin we analysed their dynamical properties, particularly in the space of integrals of motion. $Results$. We find that about 40% of the clusters likely formed in situ. A similarly large fraction, 35%, appear to be possibly associated to known merger events, in particular to $Gaia$-Enceladus (19%), the Sagittarius dwarf galaxy (5%), the progenitor of the Helmi streams (6%), and to the Sequoia galaxy (5%), although some uncertainty remains due to the degree of overlap in their dynamical characteristics. Of the remaining clusters, 16% are tentatively associated to a group with high binding energy, while the rest are all on loosely bound orbits and likely have a more heterogeneous origin. The resulting age-metallicity relations are remarkably tight and differ in their detailed properties depending on the progenitor, providing further confidence on the associations made. $Conclusions$. We provide a table listing the likely associations. Improved kinematic data by future Gaia data releases and especially a larger, systematic error-free sample of cluster ages would help to further solidify our conclusions.
△ Less
Submitted 10 September, 2019; v1 submitted 19 June, 2019;
originally announced June 2019.
-
Characterization and history of the Helmi streams with Gaia DR2
Authors:
Helmer H. Koppelman,
Amina Helmi,
Davide Massari,
Sebastian Roelenga,
Ulrich Bastian
Abstract:
The halo of the Milky Way has long been hypothesized to harbour significant amounts of merger debris. This view has been supported over more than a decade by wide-field photometric surveys which have revealed the outer halo to be lumpy. The recent release of Gaia DR2 is allowing us to establish that mergers also have been important and possibly built up the majority of the inner halo. In this work…
▽ More
The halo of the Milky Way has long been hypothesized to harbour significant amounts of merger debris. This view has been supported over more than a decade by wide-field photometric surveys which have revealed the outer halo to be lumpy. The recent release of Gaia DR2 is allowing us to establish that mergers also have been important and possibly built up the majority of the inner halo. In this work we focus on the Helmi streams, a group of streams crossing the Solar vicinity and known for almost two decades. We characterize their properties and relevance for the build-up of the Milky Way's halo. We identify new members of the Helmi streams in an unprecedented dataset with full phase-space information combining Gaia DR2, and the APOGEE DR2, RAVE DR5 and LAMOST DR4 spectroscopic surveys. Based on the orbital properties of the stars, we find new stream members up to a distance of 5 kpc from the Sun, which we characterize using photometry and metallicity information. We also perform N-body experiments to constrain the time of accretion and properties of the progenitor of the streams. We find nearly 600 new members of the Helmi streams. Their HR diagram reveals a broad age range, from approximately 11 to 13 Gyr, while their metallicity distribution goes from $\sim$ 2.3 to $\sim$1.0, and peaks at [Fe/H] $\sim$1.5. These findings confirm that the streams originate in a dwarf galaxy. Furthermore, we find 7 globular clusters to be likely associated, and which follow a well-defined age-metallicity sequence whose properties suggest a relatively massive progenitor object. Our N-body simulations favour a system with a stellar mass of $\sim 10^8\,\mathrm{M}_\odot$ accreted $5 - 8$ Gyr ago. The debris from the Helmi streams is an important donor to the MilkyWay halo, contributing approximately 15\% of its mass in field stars and 10\% of its globular clusters.
△ Less
Submitted 3 December, 2018;
originally announced December 2018.
-
The Pristine Survey IV: Approaching the Galactic metallicity floor with the discovery of an ultra metal-poor star
Authors:
Else Starkenburg,
David S. Aguado,
Piercarlo Bonifacio,
Elisabetta Caffau,
Pascale Jablonka,
Carmela Lardo,
Nicolas Martin,
Ruben Sanchez-Janssen,
Federico Sestito,
Kim A. Venn,
Kris Youakim,
Carlos Allende Prieto,
Anke Arentsen,
Marc Gentile,
Jonay I. Gonzalez Hernandez,
Collin Kielty,
Helmer H. Koppelman,
Nicolas Longeard,
Eline Tolstoy,
Raymond G. Carlberg,
Patrick Cote,
Morgan Fouesneau,
Vanessa Hill,
Alan W. McConnachie,
Julio F. Navarro
Abstract:
The early Universe presented a star formation environment that was almost devoid of heavy elements. The lowest metallicity stars thus provide a unique window into the earliest Galactic stages, but are exceedingly rare and difficult to find. Here we present the discovery of an ultra-metal-poor star, Pristine_221.8781+9.7844, using narrow-band Ca H&K photometry from the Pristine survey. Follow-up me…
▽ More
The early Universe presented a star formation environment that was almost devoid of heavy elements. The lowest metallicity stars thus provide a unique window into the earliest Galactic stages, but are exceedingly rare and difficult to find. Here we present the discovery of an ultra-metal-poor star, Pristine_221.8781+9.7844, using narrow-band Ca H&K photometry from the Pristine survey. Follow-up medium and high-resolution spectroscopy confirms the ultra-metal-poor nature of Pristine_221.8781+9.7844 ([Fe/H] = -4.66 +/- 0.13 in 1D LTE) with an enhancement of 0.3-0.4 dex in alpha-elements relative to Fe, and an unusually low carbon abundance. We derive an upper limit of A(C) = 5.6, well below typical A(C) values for such ultra metal-poor stars. This makes Pristine_221.8781+9.7844 one of the most metal-poor stars; in fact, it is very similar to the most metal-poor star known (SDSS J102915+172927). The existence of a class of ultra metal-poor stars with low(er) carbon abundances suggest that there must have been several formation channels in the early Universe through which long-lived, low-mass stars were formed.
△ Less
Submitted 11 July, 2018;
originally announced July 2018.
-
The merger that led to the formation of the Milky Way's inner stellar halo and thick disk
Authors:
Amina Helmi,
Carine Babusiaux,
Helmer H. Koppelman,
Davide Massari,
Jovan Veljanoski,
Anthony G. A. Brown
Abstract:
The assembly process of our Galaxy can be retrieved using the motions and chemistry of individual stars. Chemo-dynamical studies of the nearby halo have long hinted at the presence of multiple components such as streams, clumps, duality and correlations between the stars' chemical abundances and orbital parameters. More recently, the analysis of two large stellar surveys have revealed the presence…
▽ More
The assembly process of our Galaxy can be retrieved using the motions and chemistry of individual stars. Chemo-dynamical studies of the nearby halo have long hinted at the presence of multiple components such as streams, clumps, duality and correlations between the stars' chemical abundances and orbital parameters. More recently, the analysis of two large stellar surveys have revealed the presence of a well-populated chemical elemental abundance sequence, of two distinct sequences in the colour-magnitude diagram, and of a prominent slightly retrograde kinematic structure all in the nearby halo, which may trace an important accretion event experienced by the Galaxy. Here report an analysis of the kinematics, chemistry, age and spatial distribution of stars in a relatively large volume around the Sun that are mainly linked to two major Galactic components, the thick disk and the stellar halo. We demonstrate that the inner halo is dominated by debris from an object which at infall was slightly more massive than the Small Magellanic Cloud, and which we refer to as Gaia-Enceladus. The stars originating in Gaia-Enceladus cover nearly the full sky, their motions reveal the presence of streams and slightly retrograde and elongated trajectories. Hundreds of RR Lyrae stars and thirteen globular clusters following a consistent age-metallicity relation can be associated to Gaia-Enceladus on the basis of their orbits. With an estimated 4:1 mass-ratio, the merger with Gaia-Enceladus must have led to the dynamical heating of the precursor of the Galactic thick disk and therefore contributed to the formation of this component approximately 10 Gyr ago. These findings are in line with simulations of galaxy formation, which predict that the inner stellar halo should be dominated by debris from just a few massive progenitors.
△ Less
Submitted 31 October, 2018; v1 submitted 15 June, 2018;
originally announced June 2018.
-
One large blob and many streams frosting the nearby stellar halo in Gaia DR2
Authors:
Helmer H. Koppelman,
Amina Helmi,
Jovan Veljanoski
Abstract:
We explore the phase-space structure of nearby halo stars identified kinematically from Gaia DR2 data. We focus on their distribution in velocity and in "integrals of motion" space as well as on their photometric properties. Our sample of stars selected to be moving at a relative velocity of at least 210 km/s with respect to the Local Standard of Rest, contains an important contribution from the l…
▽ More
We explore the phase-space structure of nearby halo stars identified kinematically from Gaia DR2 data. We focus on their distribution in velocity and in "integrals of motion" space as well as on their photometric properties. Our sample of stars selected to be moving at a relative velocity of at least 210 km/s with respect to the Local Standard of Rest, contains an important contribution from the low rotational velocity tail of the disk(s). The $V_R$-distribution of these stars depicts a small asymmetry similar to that seen for the faster rotating thin disk stars near the Sun. We also identify a prominent, slightly retrograde "blob", which traces the metal-poor halo main sequence reported by Gaia Collaboration et al. (2018d). We also find many small clumps especially noticeable in the tails of the velocity distribution of the stars in our sample. Their HR diagrams disclose narrow sequences characteristic of simple stellar populations. This stream-frosting confirms predictions from cosmological simulations, namely that substructure is most apparent amongst the fastest moving stars, typically reflecting more recent accretion events.
△ Less
Submitted 11 June, 2018; v1 submitted 30 April, 2018;
originally announced April 2018.
-
Scratching the surface with TGAS and RAVE: disk moving groups in the Solar neighborhood
Authors:
Helmer H. Koppelman,
Titania Virginiflosia,
Lorenzo Posti,
Jovan Veljanoski,
Amina Helmi
Abstract:
There is long tradition extending more than a century on the identification of moving groups in the Solar neighborhood. However, with the advent of large kinematic surveys, and especially of the upcoming Gaia data releases there is a need for more sophisticated and automated substructure finders. We analyze the TGASxRAVE dataset to identify moving groups in the Galactic disk. These groups of stars…
▽ More
There is long tradition extending more than a century on the identification of moving groups in the Solar neighborhood. However, with the advent of large kinematic surveys, and especially of the upcoming Gaia data releases there is a need for more sophisticated and automated substructure finders. We analyze the TGASxRAVE dataset to identify moving groups in the Galactic disk. These groups of stars may then be used to map dynamical and star formation processes in the vicinity of the Sun. We use the ROCKSTAR algorithm, a "friends-of-friends"-like substructure finder in 6D phase-space, and analyze the Hertzsprung-Russell diagrams of the groups identified. We find 125 moving groups within 300 pc of the Sun, containing on average 50 stars, and with 3D velocity dispersions smaller than 10 km/s. Most of these groups were previously unknown. Our photometric analysis allows us to isolate a subsample of 30 statistically significant groups likely composed of stars that were born together.
△ Less
Submitted 24 April, 2018; v1 submitted 20 April, 2018;
originally announced April 2018.
-
The time evolution of gaps in tidal streams
Authors:
Amina Helmi,
Helmer H. Koppelman
Abstract:
We model the time evolution of gaps in tidal streams caused by the impact of a dark matter subhalo, while both orbit a spherical gravitational potential. To this end, we make use of the simple behaviour of orbits in action-angle space. A gap effectively results from the divergence of two nearby orbits whose initial phase-space separation is, for very cold thin streams, largely given by the impulse…
▽ More
We model the time evolution of gaps in tidal streams caused by the impact of a dark matter subhalo, while both orbit a spherical gravitational potential. To this end, we make use of the simple behaviour of orbits in action-angle space. A gap effectively results from the divergence of two nearby orbits whose initial phase-space separation is, for very cold thin streams, largely given by the impulse induced by the subhalo. We find that in a spherical potential the size of a gap increases linearly with time for sufficiently long timescales. We have derived an analytic expression that shows how the growth rate depends on the mass of the perturbing subhalo, its scale and its relative velocity with respect to the stream. We have verified these scalings using N-body simulations and find excellent agreement. For example, a subhalo of mass 10^8 Msun directly impacting a very cold thin stream on an inclined orbit can induce a gap that may reach a size of several tens of kpc after a few Gyr. The gap size fluctuates importantly with phase on the orbit, and it is largest close to pericentre. This indicates that it may not be fully straightforward to invert the spectrum of gaps present in a stream to recover the mass spectrum of the subhalos.
△ Less
Submitted 19 August, 2016; v1 submitted 28 June, 2016;
originally announced June 2016.