Galactic Stellar Halo Research Articles

Stellar Streams Discovered in the Dark Energy Survey

Abstract We perform a search for stellar streams around the Milky Way using the first 3 yr of multiband optical imaging data from the Dark Energy Survey (DES). We use DES data covering ∼5000 deg2 to a depth of g > 23.5 with a relative photometric calibration uncertainty of <1%. This data set yields unprecedented sensitivity to the stellar density field in the southern celestial hemisphere, enabling the detection of faint stellar streams to a heliocentric distance of ∼50 kpc. We search for stellar streams using a matched filter in color–magnitude space derived from a synthetic isochrone of an old, metal-poor stellar population. Our detection technique recovers four previously known thin stellar streams: Phoenix, ATLAS, Tucana III, and a possible extension of Molonglo. In addition, we report the discovery of 11 new stellar streams. In general, the new streams detected by DES are fainter, more distant, and lower surface brightness than streams detected by similar techniques in previous photometric surveys. As a by-product of our stellar stream search, we find evidence for extratidal stellar structure associated with four globular clusters: NGC 288, NGC 1261, NGC 1851, and NGC 1904. The ever-growing sample of stellar streams will provide insight into the formation of the Galactic stellar halo, the Milky Way gravitational potential, and the large- and small-scale distribution of dark matter around the Milky Way.

The Profile of the Galactic Halo from Pan-STARRS1 3π RR Lyrae

Abstract We characterize the spatial density of the Pan-STARRS1 (PS1) sample of Rrab stars to study the properties of the old Galactic stellar halo. This sample, containing 44,403 sources, spans galactocentric radii of 0.55 kpc ≤ R gc ≤ 141 kpc with a distance precision of 3% and thus is able to trace the halo out to larger distances than most previous studies. After excising stars that are attributed to dense regions such as stellar streams, the Galactic disk and bulge, and halo globular clusters, the sample contains ∼11,000 sources within 20 kpc ≤ R gc ≤ 131 kpc. We then apply forward modeling using Galactic halo profile models with a sample selection function. Specifically, we use ellipsoidal stellar density models ρ(l, b, R gc) with a constant and a radius-dependent halo flattening q(R gc). Assuming constant flattening q, the distribution of the sources is reasonably well fit by a single power law with n = 4.40 − 0.04 + 0.05 and q = 0.918 − 0.014 + 0.016 and comparably well fit by an Einasto profile with n = 9.53 − 0.28 + 0.27 , an effective radius r eff = 1.07 ± 0.10 kpc, and a halo flattening of q = 0.923 ± 0.007. If we allow for a radius-dependent flattening q(R gc), we find evidence for a distinct flattening of q ∼ 0.8 of the inner halo at ∼25 kpc. Additionally, we find that the south Galactic hemisphere is more flattened than the north Galactic hemisphere. The results of our work are largely consistent with many earlier results (e.g., Watkins et al.; Iorio et al.). We find that the stellar halo, as traced in RR Lyrae stars, exhibits a substantial number of further significant over- and underdensities, even after masking all known overdensities.

Unmixing the Galactic halo with RR Lyrae tagging

We show that tagging RR Lyrae stars according to their location in the period-amplitude diagram can be used to shed light on the genesis of the Galactic stellar halo. The mixture of RR Lyrae of ab type, separated into classes along the lines suggested by Oosterhoff, displays a strong and coherent evolution with Galactocentric radius. The change in the RR Lyrae composition appears to coincide with the break in the halo's radial density profile at ~25 kpc. Using simple models of the stellar halo, we establish that at least three different types of accretion events are necessary to explain the observed RRab behavior. Given that there exists a correlation between the RRab class fraction and the total stellar content of a dwarf satellite, we hypothesize that the field halo RRab composition is controlled by the mass of the progenitor contributing the bulk of the stellar debris at the given radius. This idea is tested against a suite of cosmological zoom-in simulations of Milky Way-like stellar halo formation. Finally, we study some of the most prominent stellar streams in the Milky Way halo and demonstrate that their RRab class fractions follow the trends established previously.

Disentangling the Galactic Halo with APOGEE. II. Chemical and Star Formation Histories for the Two Distinct Populations

Abstract The formation processes that led to the current Galactic stellar halo are still under debate. Previous studies have provided evidence for different stellar populations in terms of elemental abundances and kinematics, pointing to different chemical and star formation histories (SFHs). In the present work, we explore, over a broader range in metallicity ( ), the two stellar populations detected in the first paper of this series from metal-poor stars in DR13 of the Apache Point Observatory Galactic Evolution Experiment (APOGEE). We aim to infer signatures of the initial mass function (IMF) and the SFH from the two α-to-iron versus iron abundance chemical trends for the most APOGEE-reliable α-elements (O, Mg, Si, and Ca). Using simple chemical-evolution models, we infer the upper mass limit (M up) for the IMF and the star formation rate, and its duration for each population. Compared with the low-α population, we obtain a more intense and longer-lived SFH, and a top-heavier IMF for the high-α population.

The first all-sky view of the Milky Way stellar halo with Gaia+2MASS RR Lyrae

We exploit the first \gaia data release to study the properties of the Galactic stellar halo as traced by RR Lyrae. We demonstrate that it is possible to select a pure sample of RR Lyrae using only photometric information available in the Gaia+2MASS catalogue. The final sample contains about 21600 RR Lyrae covering an unprecedented fraction ($\sim60\%$) of the volume of the Galactic inner halo ($\text{R}<28$ kpc). We study the morphology of the stellar halo by analysing the RR Lyrae distribution with parametric and non-parametric techniques. Taking advantage of the uniform all-sky coverage, we test halo models more sophisticated than usually considered in the literature, such as those with varying flattening, tilt and/or offset of the halo with respect to the Galactic disc. A consistent picture emerges: the inner halo is well reproduced by a smooth distribution of stars settled on triaxial ellipsoids. The minor axis is perpendicular to the Milky Ways disc, while the major axis is misaligned by $\sim20^{\circ}$ from the Galactic Y axis. The elongation along the major axis is mild ($\text{p}=1.27$), and the vertical flattening is shown to evolve from a squashed state with $\text{q}\approx0.57$ in the centre to a more spherical $\text{q}\approx0.75$ at the outer edge of our dataset. The density slope is well approximated by a single power-law with exponent $\alpha=-2.96$. Within the range probed, we see no significant evidence for a change of the radial density slope, out of the plane tilt or an offset of the halo with respect to the Galaxy's centre.

An Upper Limit on the Milky Way Mass from the Orbit of the Sagittarius Dwarf Satellite

Abstract As one of the most massive Milky Way (MW) satellites, the Sagittarius dwarf galaxy played an important role in shaping the Galactic disk and stellar halo morphologies. The disruption of Sagittarius over several close-in passages has populated the halo of our Galaxy with large-scale tidal streams and offers a unique diagnostic tool for measuring its gravitational potential. Here, we test different progenitor mass models for the MW and Sagittarius by modeling the full infall of the satellite. We constrain the mass of the Galaxy based on the kinematics of the satellite remnant and multiple tidal streams of Sagittarius. Our semianalytic modeling of the orbital dynamics agrees with full N-body simulations, and favors low values for the MW mass, . This conclusion eases the tension between ΛCDM and the observed parameters of the MW satellites.

Unveiling the stellar halo with TGAS

AbstractThe detailed study of the Galactic stellar halo may hold the key to unlocking the assembly history of the Milky Way. Here, we present a machine learning model for selecting metal poor stars from the TGAS catalogue using 5 dimensional phase-space information, coupled with optical and near-IR photometry. We characterise the degree of substructure in our halo sample in the Solar neighbourhood by measuring the velocity correlation function.

Particle tagging and its implications for stellar population dynamics

We establish a controlled comparison between the properties of galactic stellar halos obtained with hydrodynamical simulations and with `particle tagging'. Tagging is a fast way to obtain stellar population dynamics: instead of tracking gas and star formation, it `paints' stars directly onto a suitably defined subset of dark matter particles in a collisionless, dark-matter-only simulation.Our study shows that there are conditions under which particle tagging generates good fits to the hydrodynamical stellar density profiles of a central Milky-Way-like galaxy and its most prominent substructure. Phase-space diffusion processes are crucial to reshaping the distribution of stars in infalling spheroidal systems and hence the final stellar halo. We conclude that the success of any particular tagging scheme hinges on this diffusion being taken into account, at a minimum by making use of `live' tagging schemes, in which particles are regularly tagged throughout the evolution of a galaxy.

Stellar halo hierarchical density structure identification using (F)OPTICS

Context. The stellar halo holds some of the best preserved fossils of Galactic formation history that can be detected as overdensities. The detection and analysis of merger by-products within the halo enables the reconstruction of the accretion history of the Milky Way. Upcoming large-scale all-sky surveys such as Gaia and The Large Synoptic Survey Telescope (LSST) will provide a huge and rich data set, which at the same time poses challenges for automated halo debris detection. Aims. We investigate the overdensity detection algorithm Ordering Points To Identify the Clustering Structure (OPTICS) as a method to identify tidal debris in the Galactic halo with large-scale surveys, as well as the variant F OPTICS which is capable of handling data sets with multi-dimensional uncertainty ellipsoids. Methods. We applied OPTICS to the a simulated Galactic stellar Halo to assess the detection performance. Additionally, we tested the performance of F OPTICS is tested by introducing uncertainty ellipsoids to the 6D phase space of two test cases. We present the Jaccard index as an alternative way to test the stability of halo debris overdensity detections without the need for a local background density estimate. Results. We optimized the OPTICS overdensity detection algorithm so that it has a slightly superlinear run-time complexity, making the method suitable for large-scale surveys. Our test on a mock galactic halo in 6D phase space shows an excellent capability to not only detect the compact dense clusters, but also the larger streams that cover a significant part of the sky. The output of OPTICS, the so-called 2D reachability diagram, proved to be a very useful tool to grasp the size, density, and substructure of the overdensities without needing to resort to complex projections of the 6D phase space. Using F OPTICS, we show the effects of introducing uncertainty ellipsoids in the 6D phase space on the retrieved tidal streams, and how the detectability of a cluster depends on whether its size and density is sufficiently large to overcome the effects of the uncertainties on the attributes.

The stellar halo in the inner Milky Way: predicted shape and kinematics

Abstract We have used N-body simulations for the Milky Way to investigate the kinematic and structural properties of the old metal-poor stellar halo in the barred inner region of the Galaxy. We find that the extrapolation of the density distribution for bulge RR Lyrae stars, ρ ∼ r−3, approximately matches the number density of RR Lyrae in the nearby stellar halo. We follow the evolution of such a tracer population through the formation and evolution of the bar and box/peanut bulge in the N-body model. We find that its density distribution changes from oblate to triaxial, and that it acquires slow rotation in agreement with recent measurements. The maximum radial velocity is ∼15–25 km s−1 at |l| = 10°–30°, and the velocity dispersion is ∼120 km s−1. Even though the simulated metal-poor halo in the bulge has a barred shape, just 12 per cent of the orbits follow the bar, and it does not trace the peanut/X structure. With these properties, the RR Lyrae population in the Galactic bulge is consistent with being the inward extension of the Galactic metal-poor stellar halo.

Carbon‐enhanced metal‐poor stars in different environments

AbstractThe origin of carbon‐enhanced metal‐poor (CEMP) stars and their possible connections with the chemical elements produced by the first stellar generations is still highly debated. We briefly review observations of CEMP stars in different environments (Galactic stellar halo, ultra‐faint, and classical dwarf galaxies) and interpret their properties using cosmological chemical‐evolution models for the formation of the Local Group. We discuss the implications of current observations for the properties of the first stars, clarify why the fraction of carbon‐enhanced to carbon‐normal stars varies in dwarf galaxies with different luminosity, and discuss the origin of the first CEMP(‐no) star found in the Sculptor dwarf galaxy.

The age of the halo as determined from halo field stars

The age of the Galactic halo is a critical parameter that can constrain the origin of the stellar halo. In general, the Galactic stellar halo is believed to be very old. However, different independent measurements and techniques based on various types of stars are required so that the age estimates of the Galactic halo are accurate, robust, and reliable. In this work, we provide a novel approach to determine the age of the halo with turn-off stars. We first carefully select 63 field halo turn-off stars from the literature. Then, we compare them with the GARSTEC model, which takes the process of atomic diffusion into account in the B − V vs. metallicity plane. Finally, we run Monte Carlo simulations to consider the uncertainty of the color index and obtain the age of 10.5 ± 1.5 Gyr. This result is in agreement with previous studies. Future works are needed to collect more turn-off samples with more accurate photometry to reduce the uncertainty of the age.

Fading features found in the kinematics of the far-reaching Milky Way stellar halo

We test the long-term kinematical stability of a Galactic stellar halo model, due to Kafle, et al. (2012), who study the kinematics of approximately 5000 blue horizontal branch (BHB) stars in the Sloan Digital Sky Survey (SDSS). The velocity dispersion $\sigma$ and anisotropy parameter $\beta$ of the stars have been determined as functions of Galactocentric radius, over the range $6<R_\mathrm{GC} < 25$ kpc, and show a strong dip in the anisotropy profile at $R_\mathrm{GC}\sim17$ kpc. By directly integrating orbits of particles in a 3-D model of the Galactic potential with these characteristics, we show that the $\sigma$ and $\beta$ profiles quickly evolve on a time scale of a $\mathrm{few}\times10$ Myr whereas the density $\rho$ profile remains largely unaffected. We suggest that the feature is therefore transient. The origin of such features in the Galactic halo remains unclear.

Strong RR Lyrae excess in the Hercules-Aquila Cloud

We map the large-scale sub-structure in the Galactic stellar halo using accurate 3D positions of ~14,000 RR Lyrae reported by the Catalina Sky Survey. In the heliocentric distance range of 10-25 kpc, in the region of the sky approximately bounded by 30{\deg} < l < 55{\deg} and -45{\deg} < b < -25{\deg}, there appears to be a strong excess of RRab stars. This overdensity, peaking at 18 kpc, is most likely associated with the so-called Hercules-Aquila Cloud, previously detected using Main Sequence tracers at similar distances in the Sloan Digital Sky Survey data. Our analysis of the period-amplitude distribution of RR Lyrae in this region indicates that the HAC is dominated by the Oosterhoff I type population. By comparing the measured RR Lyrae number density to models of a smooth stellar halo, we estimate the significance of the observed excess and provide an updated estimate of the total luminosity of the Cloud's progenitor.

Tick, Tock, Tick, Tock...

The Gaia dataset will require a huge leap forward in terms of modelling of the Milky Way. Two problems are highlighted here. First, models of the Galactic Bar remain primitive as compared to the Galactic Disk and Stellar Halo. Although Schwarzschild and N-body methods are useful, the future belongs to Made-to-Measure (M2M) models which have significant advantages in terms of storage and flexibility. Second, the Milky Way potential will need much better representation than hitherto. Most models still use very simple building blocks (Miyamoto-Nagai disks or Hernquist bulges) and these will not be fit for purpose in the Gaia Era. Expansions in terms of basis functions offer the possibility of incorporating cosmological information as priors, as well as mych greater adaptability.

Galactic stellar populations: current and new perspectives

Present studies on the evolution of the Milky Way are driven and shaped by how we conceive its stellar populations, in an on going process started by W. Baade seventy years ago. Despite much progress and advances in our understanding of these populations, inspection of their main properties is however hardly indicative of the path the Milky Way has followed to build up its mass. This is not only a matter of (stellar) age measurement, but more so the consequence of how we interprete the structures that we see in our Galaxy, often through the filter of our definitions of stellar populations. The panorama presented in the following pages opens the possibility that the present “filter” is not fully adequate. I start these Lectures with a summary of the main properties of the disks, bulge, and halo, and then present some of the new directions in the interpretation of the structure and evolution of the disk(s), with emphasis on chemical evolution. I discuss recent results in our understanding of the bulge, its stellar components and chemical evolution. Finally, I present the ideas currently proposed to explain the formation of the Galactic stellar halo. I conclude by examining how deeply all these new results question our present definition of stellar populations.

STRUCTURE AND POPULATION OF THE NGC 55 STELLAR HALO FROM A SUBARU/SUPRIME-CAM SURVEY

As part of our survey of galactic stellar halos, we investigate the structure and stellar populations of the northern outer part of the stellar halo in NGC55, a member galaxy of the Sculptor Group, using deep and wide-field V- and I-band images taken with Subaru/Suprime-Cam. Based on the analysis of the color-magnitude diagrams (CMDs) for red-giant-branch (RGB) stars, we derive a tip of RGB (TRGB)-based distance modulus to the galaxy of (m-M)_0 = 26.58 +/- 0.11 (d = 2.1 +/- 0.1 Mpc). From the stellar density maps, we detect the asymmetrically disturbed, thick disk structure and two metal-poor overdense substructures in the north region of NGC55, which may correspond to merger remnants associated with hierarchical formation of NGC55's halo. In addition, we identify a diffuse metal-poor halo extended out to at least z ~ 16 kpc from the galactic plane. The surface-brightness profiles toward the z-direction perpendicular to the galactic plane suggest that the stellar density distribution in the northern outer part of NGC55 is described by a locally isothermal disk at z <~ 6 kpc and a likely diffuse metal-poor halo with V-band surface brightness of mu_V >~ 32 mag arcsec^{-2}, where old RGB stars dominate. We derive the metallicity distributions (MDs) of these structures on the basis of the photometric comparison of RGB stars with the theoretical stellar evolutionary models. The MDs of the thick disk structures show the peak and mean metallicity of [Fe/H]peak ~ -1.4 and [Fe/H]mean ~ -1.7, respectively, while the outer substructures show more metal-poor features than the thick disk structure. Combined with the current results with our previous study for M31's halo, we discuss the possible difference in the formation process of stellar halos among different Hubble types.

When was the Large Magellanic Cloud accreted on to the Galaxy?

Using fully self-consistent N-body models for the dynamical evolution of the Large Magellanic Cloud (LMC) in the Galaxy, we show that if the LMC initially has an extended old stellar halo before its commencement of tidal interaction with the Galaxy, physical properties of the stars stripped from the LMC stellar halo can have fossil information as to when and where the LMC was accreted onto the Galaxy for the first time. If the epoch of the first LMC accretion onto the Galaxy from outside its viriral radius is more than ~4 Gyr ago (i.e., at least two pericenter passages), the stars stripped from the stellar halo of the LMC can form an irregular polar ring or a thick disk with a size of ~100 kpc and rotational kinematics. On the other hand, if the LMC was first accreted onto the Galaxy quite recently (~ 2 Gyr ago), the stripped stars form shorter leading and trailing stellar stream at R=50-120 kpc. Also distributions of the stripped stars in phase space between the two cases can be significantly different. The derived differences in structure and kinematics of the stripped stars therefore suggest that if we compare the observed three-dimensional (3D) distribution and kinematics of the outer Galactic stellar halo along the polar-axis, then we can give strong constraints on the past orbit of the LMC. We find that the orbital properties of the LMC in the successful formation models of the Magellanic Stream (MS) are consistent with those predicted from recent LambdaCDM cosmological simulations. We conclude that the LMC was accreted onto the Galaxy more than ~4 Gyr ago so that interaction between the LMC, the Small Magellanic Cloud (SMC), and the Galaxy could form the MS and its Leading Arms (LAs).

On the identification of merger debris in the Gaia era

We model the formation of the Galactic stellar halo via the accretion of satellite galaxies onto a time-dependent semi-cosmological galactic potential. Our goal is to characterize the substructure left by these accretion events in a close manner to what may be possible with the {\it Gaia} mission. We have created a synthetic {\it Gaia} Solar Neighbourhood catalogue by convolving the 6D phase-space coordinates of stellar particles from our disrupted satellites with the latest estimates of the {\it Gaia} measurement errors, and included realistic background contamination due to the Galactic disc(s) and bulge. We find that, even after accounting for the expected observational errors, the resulting phase-space is full of substructure. We are able to successfully isolate roughly 50% of the different satellites contributing to the `Solar Neighbourhood' by applying the Mean-Shift clustering algorithm in energy-angular momentum space. Furthermore, a Fourier analysis of the space of orbital frequencies allows us to obtain accurate estimates of time since accretion for approximately 30% of the recovered satellites.

THE FRACTION OF GLOBULAR CLUSTER SECOND-GENERATION STARS IN THE GALACTIC HALO

Many observational studies have revealed the presence of multiple stellar generations in Galactic globular clusters. These studies suggest that second-generation stars make up a significant fraction of the current mass of globular clusters, with the second-generation mass fraction ranging from ~50 to 80 per cent in individual clusters. In this Letter we carry out hydrodynamical simulations to explore the dependence of the mass of second-generation stars on the initial mass and structural parameters and stellar initial mass function (IMF) of the parent cluster. We then use the results of these simulations to estimate the fraction, f_{SG,H}, of the mass of the Galactic stellar halo composed of second-generation stars that originated in globular clusters. We study the dependence of f_{SG,H} on the parameters of the initial mass function of the Galactic globular cluster system. For a broad range of initial conditions, we find that the fraction of mass of the Galactic stellar halo in second-generation stars is always small, f_{SG,H}<4-6 per cent for a Kroupa-1993 IMF and f_{SG,H}<7-9 per cent for a Kroupa-2001 IMF.