Simulated Stellar Halos Research Articles

Possibilities and Limitations of Kinematically Identifying Stars from Accreted Ultra-faint Dwarf Galaxies

The Milky Way has accreted many ultra-faint dwarf galaxies (UFDs), and stars from these galaxies can be found throughout our Galaxy today. Studying these stars provides insight into galaxy formation and early chemical enrichment, but identifying them is difficult. Clustering stellar dynamics in 4D phase space (E, L z , J r , J z ) is one method of identifying accreted structure that is currently being utilized in the search for accreted UFDs. We produce 32 simulated stellar halos using particle tagging with the Caterpillar simulation suite and thoroughly test the abilities of different clustering algorithms to recover tidally disrupted UFD remnants. We perform over 10,000 clustering runs, testing seven clustering algorithms, roughly twenty hyperparameter choices per algorithm, and six different types of data sets each with up to 32 simulated samples. Of the seven algorithms, HDBSCAN most consistently balances UFD recovery rates and cluster realness rates. We find that, even in highly idealized cases, the vast majority of clusters found by clustering algorithms do not correspond to real accreted UFD remnants and we can generally only recover 6% of UFDs remnants at best. These results focus exclusively on groups of stars from UFDs, which have weak dynamic signatures compared to the background of other stars. The recoverable UFD remnants are those that accreted recently, z accretion ≲ 0.5. Based on these results, we make recommendations to help guide the search for dynamically linked clusters of UFD stars in observational data. We find that real clusters generally have higher median energy and J r , providing a way to help identify real versus fake clusters. We also recommend incorporating chemical tagging as a way to improve clustering results.

A Tilt in the Dark Matter Halo of the Galaxy

Recent observations of the stellar halo have uncovered the debris of an ancient merger, Gaia–Sausage–Enceladus (GSE), estimated to have occurred ≳8 Gyr ago. Follow-up studies have associated GSE with a large-scale tilt in the stellar halo that links two well-known stellar overdensities in diagonally opposing octants of the Galaxy (the Hercules–Aquila Cloud and Virgo Overdensity; HAC and VOD). In this paper, we study the plausibility of such unmixed merger debris persisting over several gigayears in the Galactic halo. We employ the simulated stellar halo from Naidu et al., which reproduces several key properties of the merger remnant, including the large-scale tilt. By integrating the orbits of these simulated stellar halo particles, we show that adoption of a spherical halo potential results in rapid phase mixing of the asymmetry. However, adopting a tilted halo potential preserves the initial asymmetry in the stellar halo for many gigayears. The asymmetry is preserved even when a realistic growing disk is added to the potential. These results suggest that HAC and VOD are long-lived structures that are associated with GSE and that the dark matter halo of the Galaxy is tilted with respect to the disk and aligned in the direction of HAC–VOD. Such halo–disk misalignment is common in modern cosmological simulations. Lastly, we study the relationship between the local and global stellar halo in light of a tilted global halo comprised of highly radial orbits. We find that the local halo offers a dynamically biased view of the global halo due to its displacement from the Galactic center.

The stellar halos of ETGs in the IllustrisTNG simulations: The photometric and kinematic diversity of galaxies at large radii

Context. Early-type galaxies (ETGs) are found to follow a wide variety of merger and accretion histories in cosmological simulations. Aims. We characterize the photometric and kinematic properties of simulated ETG stellar halos, and compare them to the observations. Methods. We selected a sample of 1114 ETGs in the TNG100 simulation and 80 in the higher-resolution TNG50. These ETGs span a stellar mass range of 1010.3 − 1012 M⊙ and they were selected within the range of g − r colour and λ-ellipticity diagram populated by observed ETGs. We determined photometric parameters, intrinsic shapes, and kinematic observables in their extended stellar halos. We compared the results with central IFU kinematics and ePN.S planetary nebula velocity fields at large radii, studying the variation in kinematics from center to halo, and connecting it to a change in the intrinsic shape of the galaxies. Results. We find that the simulated galaxy sample reproduces the diversity of kinematic properties observed in ETG halos. Simulated fast rotators (FRs) divide almost evenly in one third having flat λ profiles and high halo rotational support, a third with gently decreasing profiles, and another third with low halo rotation. However, the peak of rotation occurs at larger R than in observed ETG samples. Slow rotators (SRs) tend to have increased rotation in the outskirts, with half of them exceeding λ = 0.2. For M* > 1011.5 M⊙ halo rotation is unimportant. A similar variety of properties is found for the stellar halo intrinsic shapes. Rotational support and shape are deeply related: the kinematic transition to lower rotational support is accompanied by a change towards rounder intrinsic shape. Triaxiality in the halos of FRs increases outwards and with stellar mass. Simulated SRs have relatively constant triaxiality profiles. Conclusions. Simulated stellar halos show a large variety of structural properties, with quantitative but no clear qualitative differences between FRs and SRs. At the same stellar mass, stellar halo properties show a more gradual transition and significant overlap between the two families, despite the clear bimodality in the central regions. This is in agreement with observations of extended photometry and kinematics.

Quantifying the Stellar Halo's Response to the LMC's Infall with Spherical Harmonics

Abstract The vast majority of the mass in the Milky Way (MW) is in dark matter (DM); we therefore cannot directly observe the MW mass distribution and have to use tracer populations in order to infer properties of the MW DM halo. However, MW halo tracers do not only feel the gravitational influence of the MW itself. Tracers can also be affected by MW satellites; Garavito-Camargo et al. (2109) demonstrate that the Large Magellanic Cloud (LMC) induces a density wake in the MW DM, resulting in large-scale kinematic patterns in the MW stellar halo. In this work, we use spherical harmonic expansion (SHE) of the velocity fields of simulated stellar halos in an effort to disentangle perturbations on large scales (e.g., due to the LMC itself, as well as the LMC-induced DM wake) and small scales (due to substructure). Using the Garavito-Camargo et al. simulations, we demonstrate how the different terms in the SHE of the stellar velocity field reflect the different wake components and show that these signatures are a strong function of the LMC mass. An exploration of model halos built from accreted dwarfs suggests that stellar debris from massive, recent accretion events can produce much more power in the velocity angular power spectra than the perturbation from the LMC-induced wake. We therefore consider two models for the Sagittarius (Sgr) stream—the most recent, massive accretion event in the MW apart from the LMC—and find that the angular power on large scales is generally dominated by the LMC-induced wake, even when Sgr is included. We conclude that SHE of the MW stellar halo velocity field may therefore be a useful tool in quantifying the response of the MW DM halo to the LMC’s infall.

SIGNATURES OF KINEMATIC SUBSTRUCTURE IN THE GALACTIC STELLAR HALO

Tidal debris from infalling satellites can leave observable structure in the phase-space distribution of the Galactic halo. Such substructure can be manifest in the spatial and/or velocity distributions of the stars in the halo. This paper focuses on a class of substructure that is purely kinematic in nature, with no accompanying spatial features. To study its properties, we use a simulated stellar halo created by dynamically populating the Via Lactea II high-resolution N-body simulation with stars. A significant fraction of the stars in the inner halo of Via Lactea share a common speed and metallicity, despite the fact that they are spatially diffuse. We argue that this kinematic substructure is a generic feature of tidal debris from older mergers and may explain the detection of radial-velocity substructure in the inner halo made by the Sloan Extension for Galactic Understanding and Exploration. The GAIA satellite, which will provide the proper motions of an unprecedented number of stars, should further characterize the kinematic substructure in the inner halo. Our study of the Via Lactea simulation suggests that the stellar halo can be used to map the speed distribution of the local dark matter (DM) halo, which has important consequences for DM direct-detection experiments.

KINEMATICS OF THE STELLAR HALO AND THE MASS DISTRIBUTION OF THE MILKY WAY USING BLUE HORIZONTAL BRANCH STARS

Here we present a kinematic study of the Galactic halo out to a radius of $\sim$ 60 kpc, using 4664 blue horizontal branch (BHB) stars selected from the SDSS/SEGUE survey, to determine key dynamical properties. Using a maximum likelihood analysis, we determine the velocity dispersion profiles in spherical coordinates ($\sigma_{r}$, $\sigma_{\theta}$, $\sigma_{\phi}$) and the anisotropy profile ($\beta$). The radial velocity dispersion profile ($\sigma_{r}$) is measured out to a galactocentric radius of $r \sim 60$ kpc, but due to the lack of proper-motion information, $\sigma_{\theta}$, $\sigma_{\phi}$ and $\beta$ could only be derived directly out to $r \sim25$ kpc. From a starting value of $\beta\approx 0.5$ in the inner parts ($9<r/\kpc<12$), the profile falls sharply in the range $r \approx 13-18$ kpc, with a minimum value of $\beta=-1.2$ at $r=17$ kpc, rising sharply at larger radius. In the outer parts, in the range $25<r/\kpc<56$, we predict the profile to be roughly constant with a value of $\beta\approx 0.5$. The newly discovered kinematic anomalies are shown not to arise from halo substructures. We also studied the anisotropy profile of simulated stellar halos formed purely by accretion and found that they cannot reproduce the sharp dip seen in the data. From the Jeans equation, we compute the stellar rotation curve ($v_{\rm circ}$) of the Galaxy out to $r \sim 25$ kpc. The mass of the Galaxy within $r \lesssim 25$ kpc is determined to be $2.1 \times 10^{11}$ $M_{\sun}$, and with a 3-component fit to $v_{\rm circ}(r)$, we determine the virial mass of the Milky Way dark matter halo to be $M_{\rm vir} = 0.9 ^{+0.4}_{-0.3} \times 10^{12}$ $M_{\sun}$ ($R_{\rm vir} = 249^{+34}_{-31}$ kpc).

On the origin of the stellar halo and multiple stellar populations in the globular cluster NGC 1851

We propose that the observed stellar halo around the globular cluster (GC) NGC 1851 is evidence for its formation in the central region of its defunct host dwarf galaxy. We numerically investigate the long-term dynamical evolution of a nucleated dwarf galaxy embedded in a massive dark matter halo under the strong tidal field of the Galaxy. The dwarf galaxy is assumed to have a stellar nucleus (or a nuclear star cluster) that could be the progenitor for NGC 1851. We find that although the dark matter halo and the stellar envelope of the host dwarf of NGC 1851 can be almost completely stripped during its orbital evolution around the Galaxy, a minor fraction of stars in the dwarf can remain trapped by the gravitational field of the nucleus. The stripped nucleus can be observed as NGC 1851 with no/little dark matter whereas stars around the nucleus can be observed as a diffuse stellar halo around NGC 1851. The simulated stellar halo has a symmetric distribution with a power-law density slope of ~ -2 and shows no tidal tails within ~200pc from NGC 1851. We show that two GCs can merge with each other to form a new nuclear GC embedded in field stars owing to the low stellar velocity dispersion of the host dwarf. This result makes no assumption on the ages and/or chemical abundances of the two merging GCs. Thus the observed stellar halo and characteristic multiple stellar populations in NGC 1851 suggest that NGC 1851 could have formed initially in the central region of an ancient dwarf galaxy.

SUBSTRUCTURE IN THE STELLAR HALOS OF THE AQUARIUS SIMULATIONS

We characterize substructure in the simulated stellar halos of Cooper et al. (2010) which were formed by the disruption of satellite galaxies within the cosmological N-body simulations of galactic halos of the Aquarius Project. These stellar halos exhibit a wealth of tidal features: broad overdensities and very narrow faint streams akin to those observed around the Milky Way. The substructures are distributed anisotropically on the sky, a characteristic that should become apparent in the next generation of photometric surveys. The normalized RMS of the density of stars on the sky appears to be systematically larger for our halos compared to the value estimated for the Milky Way from main sequence turn-off stars in the Sloan Digital Sky Survey. We show that this is likely to be due in part to contamination by faint QSOs and redder main sequence stars, and might suggest that ~10% of the Milky Way halo stars have formed in-situ.

GALAXIA: A CODE TO GENERATE A SYNTHETIC SURVEY OF THE MILKY WAY

We present here a fast code for creating a synthetic survey of the Milky Way. Given one or more color-magnitude bounds, a survey size and geometry, the code returns a catalog of stars in accordance with a given model of the Milky Way. The model can be specified by a set of density distributions or as an N-body realization. We provide fast and efficient algorithms for sampling both types of models. As compared to earlier sampling schemes which generate stars at specified locations along a line of sight, our scheme can generate a continuous and smooth distribution of stars over any given volume. The code is quite general and flexible and can accept input in the form of a star formation rate, age metallicity relation, age velocity dispersion relation and analytic density distribution functions. Theoretical isochrones are then used to generate a catalog of stars and support is available for a wide range of photometric bands. As a concrete example we implement the Besancon Milky Way model for the disc. For the stellar halo we employ the simulated stellar halo N-body models of Bullock & Johnston (2005). In order to sample N-body models, we present a scheme that disperses the stars spawned by an N-body particle, in such a way that the phase space density of the spawned stars is consistent with that of the N-body particles. The code is ideally suited to generating synthetic data sets that mimic near future wide area surveys such as GAIA, LSST and HERMES. As an application we study the prospect of identifying structures in the stellar halo with a simulated GAIA survey. We plan to make the code publicly available at http://galaxia.sourceforge.net.

THE DUAL ORIGIN OF STELLAR HALOS. II. CHEMICAL ABUNDANCES AS TRACERS OF FORMATION HISTORY

Fully cosmological, high resolution N-Body + SPH simulations are used to investigate the chemical abundance trends of stars in simulated stellar halos as a function of their origin. These simulations employ a physically motivated supernova feedback recipe, as well as metal enrichment, metal cooling and metal diffusion. As presented in an earlier paper, the simulated galaxies in this study are surrounded by stellar halos whose inner regions contain both stars accreted from satellite galaxies and stars formed in situ in the central regions of the main galaxies and later displaced by mergers into their inner halos. The abundance patterns ([Fe/H] and [O/Fe]) of halo stars located within 10 kpc of a solar-like observer are analyzed. We find that for galaxies which have not experienced a recent major merger, in situ stars at the high [Fe/H] end of the metallicity distribution function are more [alpha/Fe]-rich than accreted stars at similar [Fe/H]. This dichotomy in the [O/Fe] of halo stars at a given [Fe/H] results from the different potential wells within which in situ and accreted halo stars form. These results qualitatively match recent observations of local Milky Way halo stars. It may thus be possible for observers to uncover the relative contribution of different physical processes to the formation of stellar halos by observing such trends in the halo populations of the Milky Way, and other local L* galaxies.

A GROUP FINDING ALGORITHM FOR MULTIDIMENSIONAL DATA SETS

We describe a density-based hierarchical group finding algorithm capable of identifying structures and substructures of any shape and density in multidimensional data sets where each dimension can be a numeric attribute with arbitrary measurement scale. This has applications in a wide variety of fields from finding structures in galaxy redshift surveys, to identifying halos and subhalos in N-body simulations and group finding in Local Group chemodynamical data sets. In general, clustering schemes require an a priori definition of a metric (a non-negative function that gives the distance between two points in a space) and the quality of clustering depends upon this choice. The general practice is to use a constant global metric which is optimal only if the clusters in the data are self-similar. For complex data configurations even the most finely tuned constant global metric turns out to be suboptimal. Moreover, the correct choice of metric also becomes increasingly important as the number of dimensions increase. To address these problems, we present an entropy-based binary space partitioning algorithm which uses a locally adaptive metric for each data point. The metric is employed to calculate the density at each point and a list of its nearest neighbors, and this information is then used to form a hierarchy of groups. Finally, the ratio of maximum to minimum density of points in a group is used to estimate the significance of the groups. Setting a threshold on this significance can effectively screen out groups arising due to Poisson noise and helps organize the groups into meaningful clusters. For a data set of N points, the algorithm requires only O(N) space and O(N(log N)3) time which makes it ideally suitable for analyzing large data sets. As an example, we apply the algorithm to identify structures in a simulated stellar halo using the full six-dimensional phase space coordinates.

The Accretion Origin of the Milky Way’s Stellar Halo

We have used data from the Sloan Digital Sky Survey (SDSS) Data Release 5 to explore the overall structure and substructure of the stellar halo of the Milky Way using ~4 million color-selected main-sequence turnoff stars with -->0.2 18.5 ? r 0.5 3.7 ? 1.2 ? 108 M?. The density profile of the stellar halo is approximately -->? r??, where ? -->2 > ? > ? 4. Yet, we found that all smooth and symmetric models were very poor fits to the distribution of stellar halo stars because the data exhibit a great deal of spatial substructure. We quantified deviations from a smooth oblate/triaxial model using the rms of the data around the model profile on scales 100 pc, after accounting for the (known) contribution of Poisson uncertainties. Within the DR5 area of the SDSS, the fractional rms deviation ?/total of the actual stellar distribution from any smooth, parameterized halo model is 40%: hence, the stellar halo is highly structured. We compared the observations with simulations of galactic stellar halos formed entirely from the accretion of satellites in a cosmological context by analyzing the simulations in the same way as the SDSS data. While the masses, overall profiles, and degree of substructure in the simulated stellar halos show considerable scatter, the properties and degree of substructure in the Milky Way's halo match well the properties of a typical stellar halo built exclusively out of the debris from disrupted satellite galaxies. Our results therefore point toward a picture in which an important fraction of the stellar halo of the Milky Way has been accreted from satellite galaxies.

Phase‐Space Distributions of Chemical Abundances in Milky Way–Type Galaxy Halos

Motivated by upcoming data from astrometric and spectroscopic surveys of the Galaxy, we explore the chemical abundance properties and phase-space distributions in hierarchically formed stellar halo simulations set in a ΛCDM universe. Our stellar halo metallicities increase with stellar halo mass. The slope of the [Fe/H]-M* trend mimics that of the satellite galaxies that were destroyed to build the halos, implying that the relation propagates hierarchically. All simulated halos contain a significant fraction of old stellar populations accreted more than 10 Gyr ago, and in a few cases some intermediate-age populations exist. In contrast with the Milky Way, many of our simulated stellar halos contain old stellar populations that are metal-rich, originating in the early accretion of massive satellites (M* ~ 109 M☉). We suggest that the (metal-rich) stellar halo of M31 falls into this category, while the more metal-poor halo of the Milky Way is lacking in early massive accretion events. Interestingly, our hierarchically formed stellar halos often have nonnegligible metallicity gradients in both [Fe/H] and [α/Fe]. These gradients extend a few tens of kiloparsecs, and can be as large as 0.5 dex in [Fe/H] and 0.2 dex in [α/Fe], with the most metal-poor halo stars typically buried within the central ~5 kpc of the galaxy. Finally, we find that chemical abundances can act as a rough substitute for time of accretion of satellite galaxies, and, based on this finding, we propose a criterion for identifying tidal streams spatially by selecting stars with [α/Fe] ratios below solar.