scholarly journals The messy merger of a large satellite and a Milky Way-like galaxy

2020 ◽  
Vol 642 ◽  
pp. L18 ◽  
Author(s):  
Helmer H. Koppelman ◽  
Roy O. Y. Bos ◽  
Amina Helmi

Aims. About 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 a suite of simulations that includes an experiment that produces a good match to the kinematics of nearby halo stars inferred from Gaia data. Methods. 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. Results. We confirm that Gaia-Enceladus probably fell in on a retrograde, 30° inclination orbit. We find that while 75% of the debris in our preferred simulation has high eccentricity (> 0.8), roughly 9% has eccentricity lower 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. Conclusions. When a massive disc galaxy undergoes a merger event, it leaves behind debris with a complex phase-space structure, a wide 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.

2020 ◽  
Vol 500 (4) ◽  
pp. 5214-5228
Author(s):  
Louise Welsh ◽  
Ryan Cooke ◽  
Michele Fumagalli

ABSTRACT We investigate the intrinsic scatter in the chemical abundances of a sample of metal-poor ([Fe/H] < −2.5) Milky Way halo stars. We draw our sample from four historic surveys and focus our attention on the stellar Mg, Ca, Ni, and Fe abundances. Using these elements, we investigate the chemical enrichment of these metal-poor stars using a model of stochastic chemical enrichment. Assuming that these stars have been enriched by the first generation of massive metal-free stars, we consider the mass distribution of the enriching population alongside the stellar mixing and explosion energy of their supernovae. For our choice of stellar yields, our model suggests that the most metal-poor stars were enriched, on average, by $\hat{N}_{\star }=5^{+13}_{-3}~(1\sigma)$ Population III stars. This is comparable to the number of enriching stars inferred for the most metal-poor DLAs. Our analysis therefore suggests that some of the lowest mass structures at z ∼ 3 contain the chemical products from < 13(2σ) Population III enriched minihaloes. The inferred IMF is consistent with that of a Salpeter distribution and there is a preference towards ejecta from minimally mixed hypernovae. However, the estimated enrichment model is sensitive to small changes in the stellar sample. An offset of ∼ 0.1 dex in the [Mg/Ca] abundance is shown to be sensitive to the inferred number of enriching stars. We suggest that this method has the potential to constrain the multiplicity of the first generation of stars, but this will require: (1) a stellar sample whose systematic errors are well understood; and, (2) documented uncertainties associated with nucleosynthetic yields.


2019 ◽  
Vol 14 (S351) ◽  
pp. 24-33
Author(s):  
Miho N. Ishigaki

AbstractI would like to review recent efforts of detailed chemical abundance measurements for field Milky Way halo stars. Thanks to the advent of wide-field spectroscopic surveys up to a several kpc from the Sun, large samples of field halo stars with detailed chemical measurements are continuously expanding. Combination of the chemical information and full six dimensional phase-space information is now recognized as a powerful tool to identify cosmological accretion events that have built a sizable fraction of the present-day stellar halo. Future observational prospects with wide-field spectroscopic surveys and theoretical prospects with supernova nucleosynthetic yields are also discussed.


2015 ◽  
Vol 11 (S317) ◽  
pp. 241-246
Author(s):  
Allyson A. Sheffield ◽  
Kathryn V. Johnston ◽  
Katia Cunha ◽  
Verne V. Smith ◽  
Steven R. Majewski

AbstractWe report chemical abundances for a sample of 66 M giants with high S/N high-resolution spectroscopy in the inner halo of the Milky Way. The program giant stars have radial velocities that vary significantly from those expected for stars moving on uniform circular orbits in the Galactic disk. Thus, based on kinematics, we expect a sample dominated by halo stars. Abundances are derived for α-elements and neutron capture elements. By analyzing the multi-dimensional abundance space, the formation site of the halo giants – in-situ or accreted – can be assessed. Of particular interest are a class of stars that form in-situ, deep in the Milky Way's gravitational potential well, but are “kicked out” of the disk into the halo due to a perturbation event. We find: (1) our sample is dominated by accreted stars and (2) tentative evidence of a small kicked-out population in our Milky Way halo sample.


2010 ◽  
Vol 62 (5) ◽  
pp. 1369-1369 ◽  
Author(s):  
Miho Ishigaki ◽  
Masashi Chiba ◽  
Wako Aoki

2010 ◽  
Vol 62 (1) ◽  
pp. 143-178 ◽  
Author(s):  
Miho Ishigaki ◽  
Masashi Chiba ◽  
Wako Aoki

2019 ◽  
Vol 14 (S353) ◽  
pp. 109-112
Author(s):  
Jiaxin Han ◽  
Wenting Wang ◽  
Zhaozhou Li

AbstractThe inference of the Milky Way halo mass requires modelling the phase space structure of dynamical tracers, with different tracers following different models and having different levels of sensitivity to the halo mass. For steady-state models, deviations from steady-state in the tracer distribution lead to an irreducible stochastic bias. This bias is small for satellite galaxies and dark matter particles, but as large as a factor of 2 for halo stars. This is consistent with the picture that satellite galaxies closely trace the underlying phase space distribution of dark matter particles, while halo stars are less phase-mixed. As a result, the use of only ~100 satellite galaxies can achieve a significantly higher accuracy than that achievable with a much larger sample of halo stars.


2020 ◽  
Vol 500 (3) ◽  
pp. 3462-3480
Author(s):  
I Marini ◽  
A Saro ◽  
S Borgani ◽  
G Murante ◽  
E Rasia ◽  
...  

ABSTRACT Cosmological N-body simulations represent an excellent tool to study the formation and evolution of dark matter (DM) haloes and the mechanisms that have originated the universal profile at the largest mass scales in the Universe. In particular, the combination of the velocity dispersion σv with the density ρ can be used to define the pseudo-entropy $S(r)=\sigma _\mathrm{v}^2/\rho ^{\, 2/3}$, whose profile is well described by a simple power law $S\propto \, r^{\, \alpha }$. We analyse a set of cosmological hydrodynamical re-simulations of massive galaxy clusters and study the pseudo-entropy profiles as traced by different collisionless components in simulated galaxy clusters: DM, stars, and substructures. We analyse four sets of simulations, exploring different resolution and physics (N-body and full hydrodynamical simulations) to investigate convergence and the impact of baryons. We find that baryons significantly affect the inner region of pseudo-entropy profiles as traced by substructures, while DM particles profiles are characterized by an almost universal behaviour, thus suggesting that the level of pseudo-entropy could represent a potential low-scatter mass-proxy. We compare observed and simulated pseudo-entropy profiles and find good agreement in both normalization and slope. We demonstrate, however, that the method used to derive observed pseudo-entropy profiles could introduce biases and underestimate the impact of mergers. Finally, we investigate the pseudo-entropy traced by the stars focusing our interest in the dynamical distinction between intracluster light and the stars bound to the brightest cluster galaxy: the combination of these two pseudo-entropy profiles is well described by a single power law out to almost the entire cluster virial radius.


2020 ◽  
Vol 493 (2) ◽  
pp. 1913-1925
Author(s):  
David M Hernandez ◽  
Sam Hadden ◽  
Junichiro Makino

ABSTRACT N-body integrations are used to model a wide range of astrophysical dynamics, but they suffer from errors which make their orbits diverge exponentially in time from the correct orbits. Over long time-scales, their reliability needs to be established. We address this reliability by running a three-body planetary system over about 200 e-folding times. Using nearby initial conditions, we can construct statistics of the long-term phase-space structure and compare to rough estimates of resonant widths of the system. We compared statistics for a wide range of numerical methods, including a Runge–Kutta method, Wisdom–Holman method, symplectic corrector methods, and a method by Laskar and Robutel. ‘Improving’ an integrator did not increase the phase-space accuracy, but simply increasing the number of initial conditions did. In fact, the statistics of a higher order symplectic corrector method were inconsistent with the other methods in one test.


2020 ◽  
Vol 501 (1) ◽  
pp. 1511-1519
Author(s):  
Junjie Luo ◽  
Weipeng Lin ◽  
Lili Yang

ABSTRACT Symplectic algorithms are widely used for long-term integration of astrophysical problems. However, this technique can only be easily constructed for separable Hamiltonian, as preserving the phase-space structure. Recently, for inseparable Hamiltonian, the fourth-order extended phase-space explicit symplectic-like methods have been developed by using the Yoshida’s triple product with a mid-point map, where the algorithm is more effective, stable and also more accurate, compared with the sequent permutations of momenta and position coordinates, especially for some chaotic case. However, it has been found that, for the cases such as with chaotic orbits of spinning compact binary or circular restricted three-body system, it may cause secular drift in energy error and even more the computation break down. To solve this problem, we have made further improvement on the mid-point map with a momentum-scaling correction, which turns out to behave more stably in long-term evolution and have smaller energy error than previous methods. In particular, it could obtain a comparable phase-space distance as computing from the eighth-order Runge–Kutta method with the same time-step.


Author(s):  
Alis J Deason ◽  
Denis Erkal ◽  
Vasily Belokurov ◽  
Azadeh Fattahi ◽  
Facundo A Gómez ◽  
...  

Abstract We use a distribution function analysis to estimate the mass of the Milky Way out to 100 kpc using a large sample of halo stars. These stars are compiled from the literature, and the vast majority ($\sim \! 98\%$) have 6D phase-space information. We pay particular attention to systematic effects, such as the dynamical influence of the Large Magellanic Cloud (LMC), and the effect of unrelaxed substructure. The LMC biases the (pre-LMC infall) halo mass estimates towards higher values, while realistic stellar halos from cosmological simulations tend to underestimate the true halo mass. After applying our method to the Milky Way data we find a mass within 100 kpc of M( < 100kpc) = 6.07 ± 0.29(stat.) ± 1.21(sys.) × 1011M⊙. For this estimate, we have approximately corrected for the reflex motion induced by the LMC using the Erkal et al. model, which assumes a rigid potential for the LMC and MW. Furthermore, stars that likely belong to the Sagittarius stream are removed, and we include a 5% systematic bias, and a 20% systematic uncertainty based on our tests with cosmological simulations. Assuming the mass-concentration relation for Navarro-Frenk-White haloes, our mass estimate favours a total (pre-LMC infall) Milky Way mass of M200c = 1.01 ± 0.24 × 1012M⊙, or (post-LMC infall) mass of M200c = 1.16 ± 0.24 × 1012 M⊙ when a 1.5 × 1011M⊙ mass of a rigid LMC is included.


Sign in / Sign up

Export Citation Format

Share Document