-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
Robustness of slow contraction to cosmic initial conditions
Authors:
Anna Ijjas,
William G. Cook,
Frans Pretorius,
Paul J. Steinhardt,
Elliot Y. Davies
Abstract:
We present numerical relativity simulations of cosmological scenarios in which the universe is smoothed and flattened by undergoing a phase of slow contraction and test their sensitivity to a wide range of initial conditions. Our numerical scheme enables the variation of all freely specifiable physical quantities that characterize the initial spatial hypersurface, such as the initial shear and spa…
▽ More
We present numerical relativity simulations of cosmological scenarios in which the universe is smoothed and flattened by undergoing a phase of slow contraction and test their sensitivity to a wide range of initial conditions. Our numerical scheme enables the variation of all freely specifiable physical quantities that characterize the initial spatial hypersurface, such as the initial shear and spatial curvature contributions as well as the initial field and velocity distributions of the scalar that drives the cosmological evolution. In particular, we include initial conditions that are far outside the perturbative regime of the well-known attractor scaling solution. We complement our numerical results by analytically performing a complete dynamical systems analysis and show that the two approaches yield consistent results.
△ Less
Submitted 8 July, 2020; v1 submitted 8 June, 2020;
originally announced June 2020.
-
Fuzzy Dark Matter Soliton Cores around Supermassive Black Holes
Authors:
Elliot Yarnell Davies,
Philip Mocz
Abstract:
We explore the effect of a supermassive black hole (SMBH) on the density profile of a fuzzy dark matter (FDM) soliton core at the centre of a dark matter halo. We numerically solve the Schrödinger-Poisson equations, treating the black hole as a gravitational point mass, and demonstrate that this additional perturbing term has a `squeezing' effect on the soliton density profile, decreasing the core…
▽ More
We explore the effect of a supermassive black hole (SMBH) on the density profile of a fuzzy dark matter (FDM) soliton core at the centre of a dark matter halo. We numerically solve the Schrödinger-Poisson equations, treating the black hole as a gravitational point mass, and demonstrate that this additional perturbing term has a `squeezing' effect on the soliton density profile, decreasing the core radius and increasing the central density. In the limit of large black hole mass, the solution approaches one akin to the hydrogen atom, with radius inversely proportional to the black hole mass. By applying our analysis to two specific galaxies (M87 and the Milky Way) and pairing it with known observational limits on the amount of centrally concentrated dark matter, we obtain a constraint on the FDM particle mass, finding that the range $10^{-22.12}$ eV $\lesssim m \lesssim 10^{-22.06}$ eV should be forbidden (taking into account additional factors concerning the life-time of the soliton in the vicinity of a black hole). Improved observational mass measurements of the black hole and total enclosed masses will significantly extend the lower-bound on the excluded FDM mass region, while self-consistent theoretical modeling of the soliton--black hole system can extend the upper-bound.
△ Less
Submitted 22 January, 2020; v1 submitted 13 August, 2019;
originally announced August 2019.
