BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Alis J. Deason

Alis J. Deason appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.GAastro-ph.COastro-ph.SRhep-ph

Catalog footprint

What is connected

28works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

28 published item(s)

preprint2022arXiv

Constraining the shape of dark matter haloes with globular clusters

We explore how diffuse stellar light and globular clusters (GCs) can be used to trace the matter distribution of their host halo using an observational methodology. For this, we use 117 simulated dark matter (DM) haloes from the $(34.4~\rm cMpc)^3$ periodic volume of the E-MOSAICS project. For each halo, we compare the stellar surface brightness and GC projected number density maps to the surface densities of DM and total mass. We find that the dominant structures identified in the stellar light and in the GCs correspond closely with those from the DM and total mass. Our method is unaffected by the presence of satellites and its precision improves with fainter GC samples. We recover tight relations between the profiles of stellar surface brightness and GC number density to those of the DM, suggesting that the profile of DM can be accurately recovered from the stars and GCs ($σ\leq0.5~$dex). We quantify the projected morphology of DM, stars and GCs, and find that the stars and GCs are more flattened than the DM. Additionally, the semi-major axes of the distribution of stars and GCs are typically misaligned by $\sim 10~$degrees from that of DM. We demonstrate that deep imaging of diffuse stellar light and GCs can place constraints on the shape, profile and orientation of their host halo. These results extend down to haloes with central galaxies $M_{\star}\geq10^{10}~M_{\odot}$, and the analysis will be applicable to future data from the Euclid, Roman and the Rubin observatories.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Dwarf stellar haloes: a powerful probe of small-scale galaxy formation and the nature of dark matter

We use N-body cosmological simulations and empirical galaxy models to study the merger history of dwarf-mass galaxies (with M_halo~10^10 M_Sun). Our input galaxy models describe the stellar mass-halo mass relation, and the galaxy occupation fraction. The number of major and minor mergers depends on the type of dark matter; in particular, minor mergers are greatly suppressed in warm dark matter models. In addition, the number of mergers that bring in stars is strongly dependent on the galaxy occupation model. For example, minor mergers are negligible for stellar halo growth in models with a high mass threshold for galaxy formation (i.e. 10^9.3 M_Sun at z=0). Moreover, this threshold for galaxy formation can also determine the relative difference (if any) between the stellar haloes of satellite and field dwarfs. Using isolated simulations of dwarf-dwarf mergers, we show that the relative frequency of major and minor mergers predict very different stellar haloes: Typically, "intermediate" dark matter merger ratios (~1:5) maximise the growth of distant stellar haloes. We discuss the observability of dwarf stellar haloes and find that the surface brightness of these features are incredibly faint. However, when several dwarfs are stacked together models that form particularly rich stellar haloes could be detectable. Finally, we show that stellar streams in the Galactic halo overlapping in phase-space with known dwarf satellites are likely remnants of their stripped stellar haloes. The mere existence of dwarf stellar haloes can already put constraints on some small-scale models, and thus observational probes should be a high priority.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Fast, Slow, Early, Late: Quenching Massive Galaxies at z~0.8

We investigate the stellar populations for a sample of 161 massive, mainly quiescent galaxies at $\langle z_{\rm obs} \rangle=0.8$ with deep Keck/DEIMOS rest-frame optical spectroscopy (HALO7D survey). With the fully Bayesian framework Prospector, we simultaneously fit the spectroscopic and photometric data with an advanced physical model (including non-parametric star-formation histories, emission lines, variable dust attenuation law, and dust and AGN emission) together with an uncertainty and outlier model. We show that both spectroscopy and photometry are needed to break the dust-age-metallicity degeneracy. We find a large diversity of star-formation histories: although the most massive ($M_{\star}>2\times10^{11}~M_{\odot}$) galaxies formed the earliest (formation redshift of $z_{\rm f}\approx5-10$ with a short star-formation timescale of $τ_{\rm SF}\lesssim1~\mathrm{Gyr}$), lower-mass galaxies have a wide range of formation redshifts, leading to only a weak trend of $z_{\rm f}$ with $M_{\star}$. Interestingly, several low-mass galaxies with have formation redshifts of $z_{\rm f}\approx5-8$. Star-forming galaxies evolve about the star-forming main sequence, crossing the ridgeline several times in their past. Quiescent galaxies show a wide range and continuous distribution of quenching timescales ($τ_{\rm quench}\approx0-5~\mathrm{Gyr}$) with a median of $\langleτ_{\rm quench}\rangle=1.0_{-0.9}^{+0.8}~\mathrm{Gyr}$ and of quenching epochs of $z_{\rm quench}\approx0.8-5.0$ ($\langle z_{\rm quench}\rangle=1.3_{-0.4}^{+0.7}$). This large diversity of quenching timescales and epochs points toward a combination of internal and external quenching mechanisms. In our sample, rejuvenation and "late bloomers" are uncommon. In summary, our analysis supports the "grow & quench" framework and is consistent with a wide and continuously-populated diversity of quenching timescales.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Snowmass2021 Cosmic Frontier White Paper: Prospects for obtaining Dark Matter Constraints with DESI

Despite efforts over several decades, direct-detection experiments have not yet led to the discovery of the dark matter (DM) particle. This has led to increasing interest in alternatives to the Lambda CDM (LCDM) paradigm and alternative DM scenarios (including fuzzy DM, warm DM, self-interacting DM, etc.). In many of these scenarios, DM particles cannot be detected directly and constraints on their properties can ONLY be arrived at using astrophysical observations. The Dark Energy Spectroscopic Instrument (DESI) is currently one of the most powerful instruments for wide-field surveys. The synergy of DESI with ESA's Gaia satellite and future observing facilities will yield datasets of unprecedented size and coverage that will enable constraints on DM over a wide range of physical and mass scales and across redshifts. DESI will obtain spectra of the Lyman-alpha forest out to z~5 by detecting about 1 million QSO spectra that will put constraints on clustering of the low-density intergalactic gas and DM halos at high redshift. DESI will obtain radial velocities of 10 million stars in the Milky Way (MW) and Local Group satellites enabling us to constrain their global DM distributions, as well as the DM distribution on smaller scales. The paradigm of cosmological structure formation has been extensively tested with simulations. However, the majority of simulations to date have focused on collisionless CDM. Simulations with alternatives to CDM have recently been gaining ground but are still in their infancy. While there are numerous publicly available large-box and zoom-in simulations in the LCDM framework, there are no comparable publicly available WDM, SIDM, FDM simulations. DOE support for a public simulation suite will enable a more cohesive community effort to compare observations from DESI (and other surveys) with numerical predictions and will greatly impact DM science.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Chemo-Dynamical Groups of Galactic Globular Clusters

We introduce a multi-component chemo-dynamical method for splitting the Galactic population of Globular Clusters (GCs) into three distinct constituents: bulge, disc, and stellar halo. The latter is further decomposed into the individual large accretion events that built up the Galactic stellar halo: the Gaia-Enceladus-Sausage, Kraken and Sequoia structures, and the Sagittarius and Helmi streams. Our modelling is extensively tested using mock GC samples constructed from the AURIGA suite of hydrodynamical simulations of Milky Way (MW)-like galaxies. We find that, on average, a proportion of the accreted GCs cannot be associated with their true infall group and are left ungrouped, biasing our recovered population numbers to approximately 80 percent of their true value. Furthermore, the identified groups have a completeness and a purity of only 65 percent. This reflects the difficulty of the problem, a result of the large degree of overlap in energy-action space of the debris from past accretion events. We apply the method to the Galactic data to infer, in a statistically robust and easily quantifiable way, the GCs associated with each MW accretion event. The resulting groups' population numbers of GCs, corrected for biases, are then used to infer the halo and stellar masses of the now defunct satellites that built up the halo of the MW.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Globular clusters as tracers of the dark matter halo: insights from the E-MOSAICS simulations

Globular clusters (GCs) are bright objects that span a wide range of galactocentric distances, and are thus probes of the structure of dark matter (DM) haloes. In this work, we explore whether the projected radial profiles of GCs can be used to infer the structural properties of their host DM haloes. We use the simulated GC populations in a sample of 166 central galaxies from the $(34.4~\rm cMpc)^3$ periodic volume of the E-MOSAICS project. We find that more massive galaxies host stellar and GC populations with shallower density profiles that are more radially extended. In addition, the metal-poor GC subpopulations tend to have shallower and more extended profiles than the metal-rich subsamples, which we relate to the preferentially accreted origin of the metal-poor GCs. We find strong correlations between the slopes and effective radii of the radial profiles of the GC populations and the structural properties of the DM haloes, such as their power-law slopes, scale radii, and concentration parameters. Accounting for a dependence on the galaxy stellar mass decreases the scatter of the two-dimensional relations. This suggests that the projected number counts of GCs, combined with their galaxy mass, trace the density profile of the DM halo of their host galaxy. When applied to extragalactic GC systems, we recover the scale radii and the extent of the DM haloes of a sample of ETGs with uncertainties smaller than $0.2~\rm dex$. Thus, extragalactic GC systems provide a novel avenue to explore the structure of DM haloes beyond the Local Group.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The mass of the Milky Way out to 100 kpc using halo stars

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 (~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(< 100 kpc) = 6.07 +/- 0.29 (stat.) +/- 1.21 (sys.) x 10^11 M_Sun. 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 M_200c = 1.01 +/- 0.24 x 10^12 M_Sun, or (post-LMC infall) mass of M_200c = 1.16 +/- 0.24 x 10^12 M_Sun when a 1.5 x 10^11 M_Sun mass of a rigid LMC is included.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

How Unusual is the Milky Way's Assembly History?

In the $Λ$CDM model of structure formation galactic haloes build up by accretion of mass and mergers of smaller subunits. The most recent massive merger event experienced by the Milky Way (MW) halo was the accretion of the Large Magellanic Cloud (LMC; which has a stellar mass of $\sim \; 10^9M_\odot$). Recent analyses of galactic stellar data from the \textit{Gaia} satellite have uncovered an earlier massive accretion event, the Gaia-Enceladus-Sausage (GES), which merged with the MW around 10 Gyr ago. Here, we use the EAGLE cosmological hydrodynamics simulation to study properties of simulated MW-mass haloes constrained to have accretion histories similar to that of the MW, specifically the recent accretion of an "LMC" galaxy and a "GES'' merger, with a quiescent period between the GES merger and the infall of the LMC (the "GES \& LMC" class). We find that $\sim 16$ per cent of MW-mass haloes have an LMC; $\sim 5$ per cent have a GES event and no further merger with an equally massive object since $z=1$; and only $0.65$ per cent belong to the LMC \& GES category. The progenitors of the MWs in this last category are much less massive than average at early times but eventually catch up with the mean. The LMC \& GES class of galaxies naturally end up in the "blue cloud" in the colour-magnitude diagram at $z=0$, tend to have a disc morphology and have a larger than average number of satellite galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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. (2019) (hereafter GC19) 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 GC19 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 Bullock & Johnston (2005) 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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Sausage & Mash: The dual origin of the Galactic thick disc and halo from the gas-rich Gaia-Enceladus-Sausage merger

We analyse a set of cosmological magneto-hydrodynamic simulations of the formation of Milky Way-mass galaxies identified to have a prominent radially anisotropic stellar halo component similar to the so-called "Gaia Sausage" found in the Gaia data. We examine the effects of the progenitor of the Sausage (the Gaia-Enceladus-Sausage, GES) on the formation of major galactic components analogous to the Galactic thick disc and inner stellar halo. We find that the GES merger is likely to have been gas-rich and contribute 10-50$\%$ of gas to a merger-induced centrally concentrated starburst that results in the rapid formation of a compact, rotationally supported thick disc that occupies the typical chemical thick disc region of chemical abundance space. We find evidence that gas-rich mergers heated the proto-disc of the Galaxy, scattering stars onto less-circular orbits such that their rotation velocity and metallicity positively correlate, thus contributing an additional component that connects the Galactic thick disc to the inner stellar halo. We demonstrate that the level of kinematic heating of the proto-galaxy correlates with the kinematic state of the population before the merger, the progenitor mass and orbital eccentricity of the merger. Furthermore, we show that the mass and time of the merger can be accurately inferred from local stars on counter-rotating orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Edge of the Galaxy

We use cosmological simulations of isolated Milky Way-mass galaxies, as well as Local Group analogues, to define the "edge" -- a caustic manifested in a drop in density or radial velocity -- of Galactic-sized haloes, both in dark matter and in stars. In the dark matter, we typically identify two caustics: the outermost caustic located at ~1.4r_200m corresponding to the "splashback" radius, and a second caustic located at ~0.6r_200m which likely corresponds to the edge of the virialized material which has completed at least two pericentric passages. The splashback radius is ill defined in Local Group type environments where the halos of the two galaxies overlap. However, the second caustic is less affected by the presence of a companion, and is a more useful definition for the boundary of the Milky Way halo. Curiously, the stellar distribution also has a clearly defined caustic, which, in most cases, coincides with the second caustic of the dark matter. This can be identified in both radial density and radial velocity profiles, and should be measurable in future observational programmes. Finally, we show that the second caustic can also be identified in the phase-space distribution of dwarf galaxies in the Local Group. Using the current dwarf galaxy population, we predict the edge of the Milky Way halo to be 292 +/- 61 kpc.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Halo's Ancient Metal-Rich Progenitor Revealed with BHB Stars

Using the data from the Sloan Digital Sky Survey and the Gaia satellite, we assemble a pure sample of $\sim$3000 Blue Horizontal Branch (BHB) stars with 7-D information, including positions, velocities and metallicities. We demonstrate that, as traced with BHBs, the Milky Way's stellar halo is largely unmixed and can not be well represented with a conventional Gaussian velocity distribution. A single-component model fails because the inner portions of the halo are swamped with metal-rich tidal debris from an ancient, head-on collision, known as the "Gaia Sausage". Motivated by the data, we build a flexible mixture model which allows us to track the evolution of the halo make-up across a wide range of radii. It is built from two components, one representing the radially anisotropic Sausage stars with their lobed velocity distribution, the other representing a more metal-poor and more isotropic component built up from minor mergers. We show that inside 25 kpc the "Sausage" contributes at least 50 % of the Galactic halo. The fraction of "Sausage" stars diminishes sharply beyond 30 kpc, which is the long-established break radius of the classical stellar halo.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Milky Way total mass profile as inferred from Gaia DR2

We determine the Milky Way (MW) mass profile inferred from fitting physically motivated models to the Gaia DR2 Galactic rotation curve and other data. Using various hydrodynamical simulations of MW-mass haloes, we show that the presence of baryons induces a contraction of the dark matter (DM) distribution in the inner regions, r<20 kpc. We provide an analytic expression that relates the baryonic distribution to the change in the DM halo profile. For our galaxy, the contraction increases the enclosed DM halo mass by factors of roughly 1.3, 2 and 4 at radial distances of 20, 8 and 1 kpc, respectively compared to an uncontracted halo. Ignoring this contraction results in systematic biases in the inferred halo mass and concentration. We provide a best-fitting contracted NFW halo model to the MW rotation curve that matches the data very well. The best-fit has a DM halo mass, $M_{200}^{\rm DM}=0.97_{-0.19}^{+0.24}\times10^{12} M_\odot$, and concentration before baryon contraction of $9.4_{-2.6}^{+1.9}$, which lie close to the median halo mass--concentration relation predicted in $Λ$CDM. The inferred total mass, $M_{200}^{\rm total}=1.08_{-0.14}^{+0.20} \times 10^{12} M_\odot$, is in good agreement with recent measurements. The model gives a MW stellar mass of $5.04_{-0.52}^{+0.43}\times10^{10} M_\odot$ and infers that the DM density at the Solar position is $ρ_{\odot}^{\rm DM}=8.8_{-0.5}^{+0.5}\times10^{-3} M_\odot \rm{pc}^{-3}\equiv0.33_{-0.02}^{+0.02}~\rm{GeV}~\rm{cm}^{-3}$. The rotation curve data can also be fitted with an uncontracted NFW halo model, but with very different DM and stellar parameters. The observations prefer the physically motivated contracted NFW halo, but the measurement uncertainties are too large to rule out the uncontracted NFW halo.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The orbital phase space of contracted dark matter halos

We study the orbital phase-space of dark matter (DM) halos in the AURIGA suite of cosmological hydrodynamics simulations of Milky Way analogues. We characterise halos by their spherical action distribution, a function of the specific angular momentum, and the radial action, of the DM particles. By comparing DM-only and hydrodynamical simulations of the same halos, we investigate the contraction of DM halos caused by the accumulation of baryons at the centre. We find a small systematic suppression of the radial action in the DM halos of the hydrodynamical simulations, suggesting that the commonly used adiabatic contraction approximation can result in an underestimate of the density by ~ 8%. We apply an iterative algorithm to contract the AURIGA DM halos given a baryon density profile and halo mass, recovering the true contracted DM profiles with an accuracy of ~15%, that reflects halo-to-halo variation. Using this algorithm, we infer the total mass profile of the Milky Way's contracted DM halo. We derive updated values for the key astrophysical inputs to DM direct detection experiments: the DM density and velocity distribution in the Solar neighbourhood.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The twisted dark matter halo of the Milky Way

We analyse systems analogous to the Milky Way (MW) in the EAGLE cosmological hydrodynamics simulation in order to deduce the likely structure of the MW's dark matter halo. We identify MW-mass haloes in the simulation whose satellite galaxies have similar kinematics and spatial distribution to those of the bright satellites of the MW, specifically systems in which the majority of the satellites (8 out of 11) have nearly co-planar orbits that are also perpendicular to the central stellar disc. We find that the normal to the common orbital plane of the co-planar satellites is well aligned with the minor axis of the host dark matter halo, with a median misalignment angle of only $17.3^\circ$. Based on this result, we infer that the minor axis of the Galactic dark matter halo points towards $(l,b)=(182^\circ,-2^\circ)$, with an angular uncertainty at the 68 and 95 percentile confidence levels of 22$^\circ$ and 43$^\circ$ respectively. Thus, the inferred minor axis of the MW halo lies in the plane of the stellar disc. The halo, however, is not homologous and its flattening and orientation vary with radius. The inner parts of the halo are rounder than the outer parts and well-aligned with the stellar disc (that is the minor axis of the halo is perpendicular to the disc). Further out, the halo twists and the minor axis changes direction by $90^\circ$. This twist occurs over a very narrow radial range and reflects variations in the filamentary network along which mass was accreted into the MW.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Simulating cosmological substructure in the solar neighbourhood

We explore the predictive power of cosmological, hydrodynamical simulations for stellar phase space substructure and velocity correlations with the Auriga simulations and Aurigaia mock-Gaia catalogues. We show that at the solar circle the Auriga simulations commonly host phase space structures in the stellar component that have constant orbital energies and arise from accreted subhaloes. These structures can persist for a few Gyrs, even after coherent streams in position space have been erased. We also explore velocity two-point correlation functions and find this diagnostic is not deterministic for particular clustering patterns in phase space. Finally, we explore these structure diagnostics with the Aurigaia catalogues and show that the current catalogues have the ability to recover some structures in phase space but careful consideration is required to separate physical structures from numerical structures arising from catalogue generation methods.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Biggest Splash

Using a large sample of bright nearby stars with accurate Gaia Data Release 2 astrometry and auxiliary spectroscopy we map out the properties of the principle Galactic components such as the "thin" and "thick" discs and the halo. We show that in the Solar neighborhood, there exists a large population of metal-rich ([Fe/H]>-0.7) stars on highly eccentric orbits. By studying the evolution of elemental abundances, kinematics and stellar ages in the plane of azimuthal velocity v_phi and metallicity [Fe/H], we demonstrate that this metal-rich halo-like component, which we dub the Splash, is linked to the alpha-rich (or "thick") disc. Splash stars have little to no angular momentum and many are on retrograde orbits. They are predominantly old, but not as old as the stars deposited into the Milky Way in the last major merger. We argue, in agreement with several recent studies, that the Splash stars may have been born in the Milky Way's proto-disc prior to the massive ancient accretion event which drastically altered their orbits. We can not, however, rule out other (alternative) formation channels. Taking advantage of the causal connection between the merger and the Splash, we put constraints of the epoch of the last massive accretion event to have finished 9.5 Gyr ago. The link between the local metal-rich and metal-poor retrograde stars is confirmed using a large suite of cutting-edge numerical simulations of the Milky Way's formation.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The little things matter: relating the abundance of ultrafaint satellites to the hosts' assembly history

Ultrafaint dwarf galaxies ($M_\star\lesssim10^{5-6}\,{\rm M}_\odot$) are relics of an early phase of galaxy formation. They contain some of the oldest and most metal-poor stars in the Universe which likely formed before the epoch of hydrogen reionisation. These galaxies are so faint that they can only be detected as satellites of the Milky Way. They are so small that they are just barely resolved in current cosmological hydrodynamics simulations. Here we combine very high resolution cosmological $N$-body simulations with a semi-analytic model of galaxy formation to study the demographics and spatial distribution of ultrafaint satellites in Milky Way-mass haloes. We show that the abundance of these galaxies is correlated with the assembly history of the host halo: at fixed mass, haloes assembled earlier contain, on average, more ultrafaint satellites today than haloes assembled later. We identify simulated galactic haloes that experience an ancient Gaia-Enceladus-Sausage-like and a recent LMC-like accretion event and find that the former occurs in 33% of the sample and the latter in 9%. Only 3% experience both events and these are especially rich in ultrafaint satellites, most acquired during the ancient accretion event. Our models predict that the radial distribution of satellites is more centrally concentrated in early-forming haloes. Accounting for the depletion of satellites by tidal interactions with the central disc, we find a very good match to the observed radial distribution of satellites in the Milky Way over the entire radial range. This agreement is mainly due to the ability of our model to track 'orphan' galaxies after their subhaloes fall below the resolution limit of the simulation.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Nature of Ionized Gas in the Milky Way Galactic Fountain

We address the spatial scale, ionization structure, mass and metal content of gas at the Milky Way disk-halo interface detected as absorption in the foreground of seven closely-spaced, high-latitude halo blue horizontal branch stars (BHBs) with heights z = 3 - 14 kpc. We detect transitions that trace multiple ionization states (e.g. CaII, FeII, SiIV, CIV) with column densities that remain constant with height from the disk, indicating that the gas most likely lies within z < 3.4 kpc. The intermediate ionization state gas traced by CIV and SiIV is strongly correlated over the full range of transverse separations probed by our sightlines, indicating large, coherent structures greater than 1 kpc in size. The low ionization state material traced by CaII and FeII does not exhibit a correlation with either N$_{\rm HI}$ or transverse separation, implying cloudlets or clumpiness on scales less than 10 pc. We find that the observed ratio log(N_SiIV/ N_CIV), with a median value of -0.69+/-0.04, is sensitive to the total carbon content of the ionized gas under the assumption of either photoionization or collisional ionization. The only self-consistent solution for photoionized gas requires that Si be depleted onto dust by 0.35 dex relative to the solar Si/C ratio, similar to the level of Si depletion in DLAs and in the Milky Way ISM. The allowed range of values for the areal mass infall rate of warm, ionized gas at the disk-halo interface is 0.0003 < dM_gas / dtdA [M_sun kpc^-2 yr^-] < 0.006. Our data support a physical scenario in which the Milky Way is fed by complex, multiphase processes at its disk-halo interface that involve kpc-scale ionized envelopes or streams containing pc-scale, cool clumps.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Isotropic at the Break? 3D Kinematics of Milky Way Halo Stars in the Foreground of M31

We present the line-of-sight (LOS) velocities for 13 distant main sequence Milky Way halo stars with published proper motions. The proper motions were measured using long baseline (5-7 years) multi-epoch HST/ACS photometry, and the LOS velocities were extracted from deep (5-6 hour integrations) Keck II/DEIMOS spectra. We estimate the parameters of the velocity ellipsoid of the stellar halo using a Markov chain Monte Carlo ensembler sampler method. The velocity second moments in the directions of the Galactic $(l,b,$ LOS) coordinate system are $\langle v^2_l \rangle^{1/2} = 138^{+43}_{-26}$ km/s, $\langle v^2_b \rangle^{1/2} = 88^{+28}_{-17}$ km/s, and $\langle v^2_{\rm{LOS}} \rangle^{1/2} = 91^{+27}_{-14}$ km/s. We use these ellipsoid parameters to constrain the velocity anisotropy of the stellar halo. Ours is the first measurement of the anisotropy parameter $β$ using 3D kinematics outside of the solar neighborhood. We find $β=-0.3^{+0.4}_{-0.9}$, consistent with isotropy and lower than solar neighborhood $β$ measurements by 2$σ$ ($β_{SN} \sim 0.5-0.7$). We identify two stars in our sample that are likely members of the known TriAnd substructure, and excluding these objects from our sample increases our estimate of the anisotropy to $β=0.1^{+0.4}_{-1.0}$, which is still lower than solar neighborhood measurements by $1σ$. The potential decrease in $β$ with Galactocentric radius is inconsistent with theoretical predictions, though consistent with recent observational studies, and may indicate the presence of large, shell-type structure (or structures) at $r \sim 25$ kpc. The methods described in this paper will be applied to a much larger sample of stars with 3D kinematics observed through the ongoing HALO7D program.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Eating Habits of Milky Way Mass Halos: Destroyed Dwarf Satellites and the Metallicity Distribution of Accreted Stars

We study the mass spectrum of destroyed dwarfs that contribute to the accreted stellar mass of Milky Way (MW) mass M_vir ~ 10^12.1 M_sun) halos using a suite of 45 zoom-in, dissipationless simulations. Empirical models are employed to relate (peak) subhalo mass to dwarf stellar mass, and we use constraints from z=0 observations and hydrodynamical simulations to estimate the metallicity distribution of the accreted stellar material. The dominant contributors to the accreted stellar mass are relatively massive dwarfs with M_star ~ 10^8-10^10 M_sun. Halos with more quiescent accretion histories tend to have lower mass progenitors (10^8-10^9 M_sun), and lower overall accreted stellar masses. Ultra-faint mass (M_star < 10^5 M_sun) dwarfs contribute a negligible amount (<< 1%) to the accreted stellar mass and, despite having low average metallicities, supply a small fraction (~2-5 %) of the very metal-poor stars with [Fe/H] < -2. Dwarfs with masses 10^5 < M_star/M_sun < 10^8 provide a substantial amount of the very metal-poor stellar material (~40-80 %), and even relatively metal-rich dwarfs with M_star > 10^8 M_sun can contribute a considerable fraction (~20-60 %) of metal-poor stars if their metallicity distributions have significant metal-poor tails. Finally, we find that the generic assumption of a quiescent assembly history for the MW halo seems to be in tension with the mass spectrum of its surviving dwarfs. We suggest that the MW could be a "transient fossil"; a quiescent halo with a recent accretion event(s) that disguises the preceding formation history of the halo.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Hubble Astrometry Initiative: Laying the Foundation for the Next-Generation Proper-Motion Survey of the Local Group

High-precision astrometry throughout the Local Group is a unique capability of the Hubble Space Telescope (HST), with potential for transformative science, including constraining the nature of dark matter, probing the epoch of reionization, and understanding key physics of galaxy evolution. While Gaia will provide unparalleled astrometric precision for bright stars in the inner halo of the Milky Way, HST is the only current mission capable of measuring accurate proper motions for systems at greater distances (> 80 kpc), which represents the vast majority of galaxies in the Local Group. The next generation of proper-motion measurements will require long time baselines, spanning many years to decades and possibly multiple telescopes, combining HST with the James Webb Space Telescope (JWST) or the Wide-Field Infrared Survey Telescope (WFIRST). However, the current HST allocation process is not conducive to such multi-cycle/multi-mission science, which will bear fruit primarily over many years. We propose an HST astrometry initiative to enable long-time-baseline, multi-mission science, which we suggest could be used to provide comprehensive kinematic measurements of all dwarf galaxies and high surface-density stellar streams in the Local Group with HST's Advanced Camera for Surveys (ACS) or Wide Field Camera 3 (WFC3). Such an initiative not only would produce forefront scientific results within the next 5 years of HST's life, but also would serve as a critical anchor point for future missions to obtain unprecedented astrometric accuracy, ensuring that HST leaves a unique and lasting legacy for decades to come.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Satellite Dwarf Galaxies in a Hierarchical Universe: Infall Histories, Group Preprocessing, and Reionization

In the Local Group (LG), almost all satellite dwarf galaxies that are within the virial radius of the Milky Way (MW) and Andromeda (M31) exhibit strong environmental influence. The orbital histories of these satellites provide the key to understanding the role of the MW/M31 halo, lower-mass groups, and cosmic reionization on the evolution of dwarf galaxies. We examine the virial-infall histories of satellites with M_star = 10^{3-9} M_sun using the ELVIS suite of cosmological zoom-in dissipationless simulations of 48 MW/M31-like halos. Satellites at z = 0 fell into the MW/M31 halos typically 5 - 8 Gyr ago at z = 0.5 - 1. However, they first fell into any host halo typically 7 - 10 Gyr ago at z = 0.7 - 1.5. This difference arises because many satellites experienced "group preprocessing" in another host halo, typically of M_vir ~ 10^{10-12} M_sun, before falling into the MW/M31 halos. Satellites with lower mass and/or those closer to the MW/M31 fell in earlier and are more likely to have experienced group preprocessing; half of all satellites with M_star < 10^6 M_sun were preprocessed in a group. Infalling groups also drive most satellite-satellite mergers within the MW/M31 halos. Finally, none of the surviving satellites at z = 0 were within the virial radius of their MW/M31 halo during reionization (z > 6), and only < 4% were satellites of any other host halo during reionization. Thus, effects of cosmic reionization versus host-halo environment on the formation histories of surviving dwarf galaxies in the LG occurred at distinct epochs, separated typically by 2 - 4 Gyr, so they are separable theoretically and, in principle, observationally.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Satellites of LMC-Mass Dwarfs: Close Friendships Ruined by Milky Way Mass Halos

Motivated by the recent discovery of several dwarfs near the Large Magellanic Cloud (LMC), we study the accretion of massive satellites onto Milky Way (MW)/M31-like halos using the ELVIS suite of N-body simulations. We identify 25 surviving LMC-mass subhalos, and investigate the lower-mass satellites that were associated with these subhalos before they fell into the MW/M31 halos. Typically, 7% of the overall z=0 satellite population of MW/M31 halos were in a surviving LMC-group before falling into the MW/M31 halo. This fraction can vary between 1% and 25%, being higher for groups with higher mass and/or more recent infall times. Groups of satellites disperse rapidly in phase space after infall, and their distances and velocities relative to the group center become statistically similar to the overall satellite population after 4-8 Gyr. We quantify the likelihood that satellites were associated with an LMC-mass group as a function of both distance and velocity relative to the LMC at z=0. The close proximity in distance of the nine Dark Energy Survey candidate dwarf galaxies to the LMC suggest that ~2-4 are likely associated with the LMC. Furthermore, if several of these dwarfs are genuine members, then the LMC-group probably fell into the MW very recently, <2 Gyr ago. If the connection with the LMC is established with follow-up velocity measurements, these "satellites of satellites" represent prime candidates to study the affects of group pre-processing on lower mass dwarfs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Touching The Void: A Striking Drop in Stellar Halo Density Beyond 50 kpc

We use A-type stars selected from Sloan Digital Sky Survey data release 9 photometry to measure the outer slope of the Milky Way stellar halo density profile beyond 50 kpc. A likelihood-based analysis is employed that models the ugr photometry distribution of blue horizontal branch (BHB) and blue straggler (BS) stars. In the magnitude range, 18.5 < g < 20.5, these stellar populations span a heliocentric distance range of: 10 kpc < D_BS < 75 kpc, 40 kpc < D_BHB < 100 kpc. Contributions from contaminants, such as QSOs, and the effect of photometric uncertainties, are also included in our modeling procedure. We find evidence for a very steep outer halo profile, with power-law index alpha ~ 6 beyond Galactocentric radii r=50 kpc, and even steeper slopes favored (alpha ~ 6-10) at larger radii. This result holds true when stars belonging to known overdensities, such as the Sagittarius stream, are included or excluded. We show that, by comparison to numerical simulations, stellar halos with shallower slopes at large distances tend to have more recent accretion activity. Thus, it is likely that the Milky Way has undergone a relatively quiet accretion history over the past several Gyr. Our measurement of the outer stellar halo profile may have important implications for dynamical mass models of the Milky Way, where the tracer density profile is strongly degenerate with total mass-estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A "Light," Centrally-Concentrated Milky Way Halo?

We discuss a novel approach to "weighing" the Milky Way dark matter halo, one that combines the latest samples of halo stars selected from the Sloan Digital Sky Survey (SDSS) with state-of-the-art numerical simulations of Milky Way analogs. The fully cosmological runs employed in the present study include "Eris", one of the highest-resolution hydrodynamical simulations of the formation of a M_vir=8e11 M_sun late-type spiral, and the dark-matter only M_vir=1.7e12 M_sun "Via Lactea II" simulation. Eris provides an excellent laboratory for creating mock SDSS samples of tracer halo stars, and we successfully compare their density, velocity anisotropy, and radial velocity dispersion profiles with the observational data. Most mock SDSS realizations show the same "cold veil" recently observed in the distant stellar halo of the Milky Way, with tracers as cold as sigma_los ~ 50 km/s between 100-150 kpc. Controlled experiments based on the integration of the spherical Jeans equation as well as a particle tagging technique applied to Via Lactea II show that a "heavy" M_vir 2e12 M_sun realistic host produces a poor fit to the kinematic SDSS data. We argue that these results offer added evidence for a "light," centrally-concentrated Milky Way halo.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Stellar mass-gap as a probe of halo assembly history and concentration: youth hidden among old fossils

We investigate the use of the halo mass-gap statistic --- defined as the logarithmic difference in mass between the host halo and its most massive satellite subhalo --- as a probe of halo age and concentration. A cosmological N-body simulation is used to study N ~25, 000 group/cluster sized halos in the mass range 10^12.5 < M_halo/M_sol < 10^14.5. In agreement with previous work, we find that halo mass-gap is related to halo formation time and concentration. On average, older and more highly concentrated halos have larger halo mass-gaps, and this trend is stronger than the mass-concentration relation over a similar dynamic range. However, there is a large amount of scatter owing to the transitory nature of the satellite subhalo population, which limits the use of the halo mass-gap statistic on an object-by-object basis. For example, we find that 20% of very large halo mass-gap systems (akin to "fossil groups") are young, and have likely experienced a recent merger between a massive satellite subhalo and the central subhalo. We relate halo mass-gap to the observable stellar mass-gap via abundance matching. Using a galaxy group catalog constructed from the Sloan Digital Sky Survey Data Release 7, we find that the star formation and structural properties of galaxies at fixed mass show no trend with stellar mass-gap. This is despite a variation in halo-age of ~2.5 Gyr over ~1.2 dex in stellar mass-gap. Thus, we find no evidence to suggest that the halo formation history significantly affects galaxy properties.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Global structure and kinematics of stellar haloes in cosmological hydrodynamic simulations

We use the Galaxies-Intergalactic Medium Interaction Calculation (GIMIC) suite of cosmological hydrodynamical simulations to study the global structure and kinematics of stellar spheroids of Milky Way mass disc galaxies. Font et al. have recently demonstrated that these simulations are able to successfully reproduce the satellite luminosity functions and the metallicity and surface brightness profiles of the spheroids of the Milky Way and M31. A key to the success of the simulations is a significant contribution to the spheroid from stars that formed in situ. While the outer halo is dominated by accreted stars, stars formed in the main progenitor of the galaxy dominate at r < ~30 kpc. In the present study we show that this component was primarily formed in a proto-disc at high redshift and was subsequently liberated from the disc by dynamical heating associated with mass accretion. As a consequence of its origin, the in situ component of the spheroid has different kinematics (namely net prograde rotation with respect to the disc) than that of the spheroid component built from the disruption of satellites. In addition, the in situ component has a flattened distribution, that is due in part to its rotation. We make comparisons with measurements of the shape and kinematics of local galaxies, including the Milky Way and M31, and stacked observations of more distant galaxies. We find that the simulated disc galaxies have spheroids of the correct shape (oblate with a median axis ratio of ~0.6 at radii of < ~30 kpc, but note there is significant system-to-system scatter in this quantity) and that the kinematics show evidence for two components (due to in situ vs. accreted), as observed. Our findings therefore add considerable weight to the importance of dissipative processes in the formation of stellar haloes and to the notion of a 'dual stellar halo'.

0 citations0 reviewsNo code linkedOpen access