Dolag, K., Sorce, J. G., Pilipenko, S., Hernández-Martínez, E., Valentini, M., Gottlöber, S., Aghanim, N., Khabibullin, I., 2023, Astronomy and Astrophysics
, 677 , A169 Published: September 2023
Context. Several observations of the Local Universe point toward the existence of very prominent structures: massive galaxy clusters and local superclusters on the one hand, but also large local voids and underdensities on the other. However, it is highly nontrivial to connect such different observational selected tracers to the underlying dark matter (DM) distribution. Aims: Therefore, constructing mock catalogs of such observable tracers using cosmological hydrodynamics simulations is needed. These simulations have to follow galaxy formation physics and also have to be constrained to reproduce the Local Universe. Such constraints should be based on observables that directly probe the full underlying gravitational field, such as the observed peculiar velocity field, to provide an independent test on the robustness of these distinctive structures. Methods: We used a 500 h‒1 Mpc constrained simulation of the Local Universe to investigate the anomalies in the local density field, as found in observations. Constructing the initial conditions based on peculiar velocities derived from the CosmicFlows-2 catalog makes the predictions of the simulations completely independent from the distribution of the observed tracer population, and following galaxy formation physics directly in the hydrodynamics simulations also allows the comparison to be based directly on the stellar masses of galaxies or X-ray luminosity of clusters. We also used the 2668 h‒1 Mpc large cosmological box from the Magneticum simulations to evaluate the frequency of finding such anomalies in random patches within simulations. Results: We demonstrate that halos and galaxies in our constrained simulation trace the local dark matter density field very differently. Thus, this simulation reproduces the observed 50% underdensity of galaxy clusters and groups within the sphere of ≈100 Mpc when applying the same mass or X-ray luminosity limit used in the observed cluster sample (CLASSIX), which is consistent with a ≈1.5σ feature. At the same time, the simulation reproduces the observed overdensity of massive galaxy clusters within the same sphere, which on its own also corresponds to a ≈1.5σ feature. Interestingly, we find that only 44 out of 15 635 random realizations (i.e., 0.28%) match both anomalies, thus making the Local Universe a ≈3σ environment. We finally compared a mock galaxy catalog with the observed distribution of galaxies in the Local Universe, finding a match to the observed factor of 2 overdensity at ∼16 Mpc as well as the observed 15% underdensity at ∼40 Mpc. Conclusions: Constrained simulations of the Local Universe which reproduce the main features of the local density field open a new window for local field cosmology, where the imprint of the specific density field and the impact on the bias through the observational specific tracers can be investigated in detail.
Khoperskov, S., Minchev, I., Libeskind, N., Belokurov, V., Steinmetz, M., Gomez, F. A., Grand, R. J. J., Hoffman, Y., Knebe, A., Sorce, J. G., Spaare, M., Tempel, E., Vogelsberger, M., 2023, Astronomy and Astrophysics
, 677 , A91 Published: September 2023
Stellar chemical abundances and kinematics provide key information for recovering the assembly history of galaxies. In this work we explore the chemo-chrono-kinematics of accreted and in situ stellar populations, by analyzing six M31/Milky Way (MW) analogues from the HESTIA suite of cosmological hydrodynamics zoom-in simulations of the Local Group. We show that elemental abundances ([Fe/H], [Mg/Fe]) of merger debris in the stellar haloes are chemically distinct from the survived dwarf galaxies, in that they are [α/Fe]-enhanced and have lower metallicity in the same stellar mass range. Therefore, mergers debris have abundances expected for stars originating from dwarfs that had their star formation activity quenched at early times. Accreted stellar haloes, including individual debris, reveal [Fe/H] and [Mg/Fe] gradients in the E ‒ Lz plane, with the most metal-rich, [α/Fe]-poor stars, which have formed in the inner parts of the disrupted systems before the merger, contributing mainly to the central regions of the host galaxies. This results in negative metallicity gradients in the accreted components of stellar haloes at z = 0, seen also for the individual merger debris. We suggest, therefore, that abundance measurements of halo stars in the inner MW will allow constraining better the parameters, such as the mass and merger time, of MW's most massive merger Gaia-Sausage-Enceladus. The metallicity distribution functions (MDFs) of the individual debris show several peaks and the majority of debris have lower metallicity than the in situ stars in the prograde part of the E ‒ Lz space. At the same time, non-rotating and retrograde accreted populations are very similar to the in situ stars in terms of [Fe/H] abundance. Prograde accreted stars show a prominent knee in the [Fe/H]-[Mg/Fe] plane, reaching up to solar [Mg/Fe], while retrograde stars typically contribute to the high-[Mg/Fe] sequence only. We find that the most metal-poor stars ([Fe/H] ≲ ‒1) of the HESTIA galaxies exhibit net rotation up to 80 km s‒1, which is consistent with the Aurora population recently identified in the MW. At higher metallicities ([Fe/H] ≈ ‒0.5 ± 0.1) we detect a sharp transition (spin-up) from the turbulent phase to a regular disk-like rotation. Different merger debris appear similar in the [Fe/H]-[Mg/Fe] plane, thus making it difficult to identify individual events. However, combining a set of abundances, and especially stellar age, makes it possible to distinguish between different debris.
Khoperskov, S., Minchev, I., Libeskind, N., Haywood, M., Di Matteo, P., Belokurov, V., Steinmetz, M., Gomez, F. A., Grand, R. J. J., Hoffman, Y., Knebe, A., Sorce, J. G., Spaare, M., Tempel, E., Vogelsberger, M., 2023, Astronomy and Astrophysics
, 677 , A90 Published: September 2023
Recent progress in understanding the assembly history of the Milky Way (MW) is driven by the tremendous amount of high-quality data delivered by Gaia (ESA), revealing a number of substructures potentially linked to several ancient accretion events. In this work we aim to explore the phase-space structure of accreted stars by analysing six M31/MW analogues from the HESTIA suite of cosmological hydrodynamics zoom-in simulations of the Local Group. We find that all HESTIA galaxies experience a few dozen mergers but only between one and four of those have stellar mass ratios > 0.2, relative to the host at the time of the merger. Depending on the halo definition, the most massive merger contributes from 20% to 70% of the total stellar halo mass. Individual merger remnants show diverse density distributions at z = 0, significantly overlapping with each other and with the in situ stars in the Lz ‒ E, (VR, Vϕ) and (R, vϕ) coordinates. Moreover, merger debris often shifts position in the Lz ‒ E space with cosmic time due to the galactic mass growth and the non-axisymmetry of the potential. In agreement with previous works, we show that even individual merger debris exhibit a number of distinct Lz ‒ E features. In the (VR, Vϕ) plane, all HESTIA galaxies reveal radially hot, non-rotating or weakly counter-rotating, Gaia-Sausage-like features, which are the remnants of the most recent significant mergers. We find an age gradient in Lz ‒ E space for individual debris, where the youngest stars, formed in the inner regions of accreting systems, deposit to the innermost regions of the host galaxies. The bulk of these stars formed during the last stages of accretion, making it possible to use the stellar ages of the remnants to date the merger event. In action space (Jr, Jz, Jϕ), merger debris do not appear as isolated substructures, but are instead scattered over a large parameter area and overlap with the in situ stars. We suggest that accreted stars can be best identified using Jr > 0.2‒0.3(104 kpc km s‒1)0.5. We also introduce a new, purely kinematic space (Jz/Jr-orbital eccentricity), where different merger debris can be disentangled better from each other and from the in situ stars. Accreted stars have a broad distribution of eccentricities, peaking at ϵ ≈ 0.6 ‒ 0.9, and their mean eccentricity tends to be smaller for systems accreted more recently.
Khoperskov, S., Minchev, I., Libeskind, N., Haywood, M., Di Matteo, P., Belokurov, V., Steinmetz, M., Gomez, F. A., Grand, R. J. J., Hoffman, Y., Knebe, A., Sorce, J. G., Spaare, M., Tempel, E., Vogelsberger, M., 2023, Astronomy and Astrophysics
, 677 , A89 Published: September 2023
Theory suggests that mergers play an important role in shaping galactic discs and stellar haloes, which was observationally confirmed in the Milky Way (MW) thanks to Gaia data. In this work, aiming to probe the contribution of mergers to the in situ stellar halo formation, we analyse six M 31 and MW analogues from the HESTIA suite of cosmological hydrodynamical zoom-in simulations of the Local Group. We found that all the HESTIA galaxies experience between one to four mergers with stellar mass ratios between 0.2 and 1 relative to the host at the time of the merger. These significant mergers, with a single exception, happened 7 ‒ 11 Gyr ago. The overall impact of the most massive mergers in HESTIA is clearly seen as a sharp increase in the orbital eccentricity (and a corresponding decrease in the rotational velocity Vϕ) of pre-existing disc stars of the main progenitor, thus nicely reproducing the Splash-, Plume-like feature that was discovered in the MW. We do find a correlation between mergers and close pericentric passages of massive satellites and bursts of the star formation in the in situ component. Massive mergers sharply increase the disc velocity dispersion of the in situ stars; however, the latest significant merger often heats up the disc up to the numbers when the contribution of the previous ones is less prominent in the age-velocity dispersion relation. In HESTIA galaxies, the in situ halo is an important component of the inner stellar halo where its fraction is about 30 ‒ 40%, while in the outer parts it typically does not exceed ≈5% beyond 15 kpc from the galactic centre. The simulations suggest that this component of the stellar haloes continues to grow well after mergers conclude; however, the most significant contribution comes from stars that formed recently before the merger. The orbital analysis of the HESTIA galaxies suggests that wedges in Rmax ‒ Zmax (apocentre - maximum height from the mid-plane) space are mainly populated by the stars born in between significant mergers.
Pfeifer, S., Valade, A., Gottlöber, S., Hoffman, Y., Libeskind, N. I., Hellwing, W. A., 2023, Monthly Notices of the Royal Astronomical Society
, 523, 4 , 5985 Published: August 2023
The aim of cosmological simulations is to reproduce the properties of the observed Universe, serving as tools to test structure and galaxy formation models. Constrained simulations of our local cosmological region up to a few hundred $h^{-1}\, \mbox{Mpc}$ , the Local Universe, are designed to reproduce the actual cosmic web of structures as observed. A question that often arises is how to judge the quality of constrained simulations against the observations of the Local Universe. Here we introduce the Local Universe model (LUM), a new methodology, whereby many constrained simulations can be judged and the 'best' initial conditions can be identified. By characterizing the Local Universe as a set of rich clusters, the model identifies haloes that serve as simulated counterparts to the observed clusters. Their merit is determined against a null hypothesis, the probability that such a counterpart could be identified in a random, unconstrained simulation. This model is applied to 100 constrained simulations using the Cosmicflows-3 data. Cluster counterparts are found for all constrained simulations, their distribution of separation from the true observed cluster position and their mass distribution are investigated. Lastly, the 'best' constrained simulation is selected using the LUM and discussed in more detail.
