Showing 1–50 of 82 results for author: Deason, A J

The accreted stellar haloes of Milky Way-mass galaxies as a probe of the nature of the dark matter

Authors: Victor J. Forouhar Moreno, Azadeh Fattahi, Alis J. Deason, Fergus Henstridge, Alejandro Benítez-Llambay

Abstract: Galactic stellar haloes are largely composed of the remnants of galaxies accreted during the assembly of their host galaxies, and hence their properties reflect the mass spectrum and post-accretion evolution of their satellites. As the nature of dark matter (DM) can affect both, we explore how the properties of the accreted stellar component vary across cold (CDM), warm (WDM) and self-interacting… ▽ More Galactic stellar haloes are largely composed of the remnants of galaxies accreted during the assembly of their host galaxies, and hence their properties reflect the mass spectrum and post-accretion evolution of their satellites. As the nature of dark matter (DM) can affect both, we explore how the properties of the accreted stellar component vary across cold (CDM), warm (WDM) and self-interacting (SIDM) models. We do this by studying accreted stellar populations around eight MW-mass haloes using cosmological hydrodynamical simulations based on the EAGLE galaxy formation model, in which we find that the accreted stellar mass remains similar across models. Contrary to WDM, which only presents minor differences relative to CDM, the distribution of accreted stars in SIDM changes significantly within $0.05R_{200}$ ($10\,\mathrm{kpc}$). The central density reduces to $\langle ρ^{\mathrm{SIDM}}_{\mathrm{exsitu}} / ρ^{\mathrm{CDM}}_{\mathrm{exsitu}} \rangle = 0.3$ and has a shallower radial dependence, with logarithmic density slopes of $\langle α_{\mathrm{SIDM}} \rangle = -1.4$ vs $\langle α_{\mathrm{CDM}} \rangle = -1.7$. Additionally, stars are on more tangential orbits than their CDM counterparts, with a change in the velocity anisotropy of $\langle Δβ\rangle = - 0.2$. Finally, SIDM stellar haloes have the largest number and prominence of overdensities in radius vs radial velocity space. This is due to a combination of shorter stellar halo progenitor merging timescales and shallower host potentials, with the former resulting in less time for dynamical friction and radialisation to operate. In summary, we show that the phase-space structure of Galactic stellar haloes encode key information that can be used to distinguish and rule out different DM models. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: 13 pages, 8 figures. Submitted to MNRAS

Galactic Archaeology with Gaia

Authors: Alis J. Deason, Vasily Belokurov

Abstract: The Gaia mission has revolutionized our view of the Milky Way and its satellite citizens. The field of Galactic Archaeology has been piecing together the formation and evolution of the Galaxy for decades, and we have made great strides, with often limited data, towards discovering and characterizing the subcomponents of the Galaxy and its building blocks. Now, the exquisite 6D phase-space plus che… ▽ More The Gaia mission has revolutionized our view of the Milky Way and its satellite citizens. The field of Galactic Archaeology has been piecing together the formation and evolution of the Galaxy for decades, and we have made great strides, with often limited data, towards discovering and characterizing the subcomponents of the Galaxy and its building blocks. Now, the exquisite 6D phase-space plus chemical information from Gaia and its complementary spectroscopic surveys has handed us a plethora of data to pour over as we move towards a quantitative rather than qualitative view of the Galaxy and its progenitors. We review the state of the field in the post-Gaia era, and examine the key lessons that will dictate the future direction of Galactic halo research. △ Less

Submitted 1 July, 2024; v1 submitted 19 February, 2024; originally announced February 2024.

Comments: Accepted Review article for New Astronomy Reviews (as part of a Special Issue: "Gaia, the first crop of discoveries")

Can we really pick and choose? Benchmarking various selections of Gaia Enceladus/Sausage stars in observations with simulations

Authors: Andreia Carrillo, Alis J. Deason, Azadeh Fattahi, Thomas M. Callingham, Robert J. J. Grand

Abstract: Large spectroscopic surveys plus Gaia astrometry have shown us that the inner stellar halo of the Galaxy is dominated by the debris of Gaia Enceladus/Sausage (GES). With the richness of data at hand, there are a myriad of ways these accreted stars have been selected. We investigate these GES selections and their effects on the inferred progenitor properties using data constructed from APOGEE and G… ▽ More Large spectroscopic surveys plus Gaia astrometry have shown us that the inner stellar halo of the Galaxy is dominated by the debris of Gaia Enceladus/Sausage (GES). With the richness of data at hand, there are a myriad of ways these accreted stars have been selected. We investigate these GES selections and their effects on the inferred progenitor properties using data constructed from APOGEE and Gaia. We explore selections made in eccentricity, energy-angular momentum (E-Lz), radial action-angular momentum (Jr-Lz), action diamond, and [Mg/Mn]-[Al/Fe] in the observations, selecting between 144 and 1,279 GES stars with varying contamination from in-situ and other accreted stars. We also use the Auriga cosmological hydrodynamic simulations to benchmark the different GES dynamical selections. Applying the same observational GES cuts to nine Auriga galaxies with a GES, we find that the Jr-Lz method is best for sample purity and the eccentricity method for completeness. Given the average metallicity of GES (-1.28 < [Fe/H] < -1.18), we use the $z=0$ mass-metallicity relationship to find an average $\rm M_{\star}$ of $\sim 4 \times 10^{8}$ $\rm M_{\odot}$. We adopt a similar procedure and derive $\rm M_{\star}$ for the GES-like systems in Auriga and find that the eccentricity method overestimates the true $\rm M_{\star}$ by $\sim2.6\times$ while E-Lz underestimates by $\sim0.7\times$. Lastly, we estimate the total mass of GES to be $\rm 10^{10.5 - 11.1}~M_{\odot}$ using the relationship between the metallicity gradient and the GES-to-in-situ energy ratio. In the end, we cannot just `pick and choose' how we select GES stars, and instead should be motivated by the science question. △ Less

Submitted 1 June, 2023; originally announced June 2023.

Comments: 20 pages, 14 figures, submitted to MNRAS

The Hubble Space Telescope Survey of M31 Satellite Galaxies II. The Star Formation Histories of Ultra-Faint Dwarf Galaxies

Authors: A. Savino, D. R. Weisz, E. D. Skillman, A. Dolphin, A. A. Cole, N. Kallivayalil, A. Wetzel, J. Anderson, G. Besla, M. Boylan-Kolchin, T. M. Brown, J. S. Bullock, M. L. M. Collins, M. C. Cooper, A. J. Deason, A. L. Dotter, M. Fardal, A. M. N. Ferguson, T. K. Fritz, M. C. Geha, K. M. Gilbert, P. Guhathakurta, R. Ibata, M. J. Irwin, M. Jeon , et al. (12 additional authors not shown)

Abstract: We present the lifetime star formation histories (SFHs) for six ultra-faint dwarf (UFD; $M_V>-7.0$, $ 4.9<\log_{10}({M_*(z=0)}/{M_{\odot}})<5.5$) satellite galaxies of M31 based on deep color-magnitude diagrams constructed from \textit{Hubble Space Telescope} imaging. These are the first SFHs obtained from the oldest main sequence turn-off of UFDs outside the halo of the Milky Way (MW). We find th… ▽ More We present the lifetime star formation histories (SFHs) for six ultra-faint dwarf (UFD; $M_V>-7.0$, $ 4.9<\log_{10}({M_*(z=0)}/{M_{\odot}})<5.5$) satellite galaxies of M31 based on deep color-magnitude diagrams constructed from \textit{Hubble Space Telescope} imaging. These are the first SFHs obtained from the oldest main sequence turn-off of UFDs outside the halo of the Milky Way (MW). We find that five UFDs formed at least 50\% of their stellar mass by $z=5$ (12.6~Gyr ago), similar to known UFDs around the MW, but that 10-40\% of their stellar mass formed at later times. We uncover one remarkable UFD, \A{XIII}, which formed only 10\% of its stellar mass by $z=5$, and 75\% in a rapid burst at $z\sim2-3$, a result that is robust to choices of underlying stellar model and is consistent with its predominantly red horizontal branch. This "young" UFD is the first of its kind and indicates that not all UFDs are necessarily quenched by reionization, which is consistent with predictions from several cosmological simulations of faint dwarf galaxies. SFHs of the combined MW and M31 samples suggest reionization did not homogeneously quench UFDs. We find that the least massive MW UFDs ($M_*(z=5) \lesssim 5\times10^4 M_{\odot}$) are likely quenched by reionization, whereas more massive M31 UFDs ($M_*(z=5) \gtrsim 10^5 M_{\odot}$) may only have their star formation suppressed by reionization and quench at a later time. We discuss these findings in the context of the evolution and quenching of UFDs. △ Less

Submitted 13 September, 2023; v1 submitted 22 May, 2023; originally announced May 2023.

Comments: 18 pages, 14 figures, 5 appendices, accepted for publication in ApJ

HALO7D III: Chemical Abundances of Milky Way Halo Stars from Medium Resolution Spectra

Authors: Kevin A. McKinnon, Emily C. Cunningham, Constance M. Rockosi, Puragra Guhathakurta, Ivanna Escala, Evan N. Kirby, Alis J. Deason

Abstract: The Halo Assembly in Lambda Cold Dark Matter: Observations in 7 Dimensions (HALO7D) survey measures the kinematics and chemical properties of stars in the Milky Way (MW) stellar halo to learn about the formation of our Galaxy. HALO7D consists of Keck II/DEIMOS spectroscopy and Hubble Space Telescope-measured proper motions of MW halo main sequence turn-off (MSTO) stars in the four CANDELS fields.… ▽ More The Halo Assembly in Lambda Cold Dark Matter: Observations in 7 Dimensions (HALO7D) survey measures the kinematics and chemical properties of stars in the Milky Way (MW) stellar halo to learn about the formation of our Galaxy. HALO7D consists of Keck II/DEIMOS spectroscopy and Hubble Space Telescope-measured proper motions of MW halo main sequence turn-off (MSTO) stars in the four CANDELS fields. HALO7D consists of deep pencil beams, making it complementary to other contemporary wide-field surveys. We present the [Fe/H] and [$α$/Fe] abundances for 113 HALO7D stars in the Galactocentric radial range of $\sim 10-40$ kpc. Using the full 7D chemodynamical data (3D positions, 3D velocities, and abundances) of HALO7D, we measure the velocity anisotropy, $β$, of the halo velocity ellipsoid for each field and for different metallicity-binned subsamples. We find that two of the four fields have stars on very radial orbits, while the remaining two have stars on more isotropic orbits. Separating the stars into high, mid, and low [Fe/H] bins at $-2.2$ dex and $-1.1$ dex for each field separately, we find differences in the anisotropies between the fields and between the bins; some fields appear dominated by radial orbits in all bins while other fields show variation between the [Fe/H] bins. These chemodynamical differences are evidence that the HALO7D fields have different fractional contributions from the progenitors that built up the MW stellar halo. Our results highlight the additional information that is available on smaller spatial scales when compared to results from a spherical average of the stellar halo. △ Less

Submitted 9 May, 2023; v1 submitted 14 February, 2023; originally announced February 2023.

Comments: 34 pages, 16 figures

Unravelling the mass spectrum of destroyed dwarf galaxies with the metallicity distribution function

Authors: Alis J. Deason, Sergey E. Koposov, Azadeh Fattahi, Robert J. J. Grand

Abstract: Accreted stellar populations are comprised of the remnants of destroyed galaxies, and often dominate the `stellar haloes' of galaxies such as the Milky Way (MW). This ensemble of external contributors is a key indicator of the past assembly history of a galaxy. We introduce a novel statistical method that uses the unbinned metallicity distribution function (MDF) of a stellar population to estimate… ▽ More Accreted stellar populations are comprised of the remnants of destroyed galaxies, and often dominate the `stellar haloes' of galaxies such as the Milky Way (MW). This ensemble of external contributors is a key indicator of the past assembly history of a galaxy. We introduce a novel statistical method that uses the unbinned metallicity distribution function (MDF) of a stellar population to estimate the mass spectrum of its progenitors. Our model makes use of the well-known mass-metallicity relation of galaxies and assumes Gaussian MDF distributions for individual progenitors: the overall MDF is thus a mixture of MDFs from smaller galaxies. We apply the method to the stellar halo of the MW, as well as the classical MW satellite galaxies. The stellar components of the satellite galaxies have relatively small sample sizes, but we do not find any evidence for accreted populations with L > L_host/100. We find that the MW stellar halo has N~1-3 massive progenitors (L > 10^8 L_Sun) within 10 kpc, and likely several hundred progenitors in total. We also test our method on simulations of MW-mass haloes, and find that our method is able to recover the true accreted population within a factor of two. Future datasets will provide MDFs with orders of magnitude more stars, and this method could be a powerful technique to quantify the accreted populations down to the ultra-faint dwarf mass-scale for both the MW and its satellites. △ Less

Submitted 28 February, 2023; v1 submitted 11 January, 2023; originally announced January 2023.

Comments: 13 pages, 11 figures. MNRAS in press

The progenitor galaxies of stellar haloes as "failed" Milky Ways

Authors: Sownak Bose, Alis J. Deason

Abstract: The stellar halo of the Milky Way records the history of its interactions with dwarf galaxies, whose subsequent destruction results in the formation of an extended stellar component. Recent works have suggested that galaxies with masses comparable to the Large Magellanic Cloud (LMC, $M_\star \sim 10^9\,{\rm M}_\odot$) may be the primary building blocks of the stellar halo of our Galaxy. We use cos… ▽ More The stellar halo of the Milky Way records the history of its interactions with dwarf galaxies, whose subsequent destruction results in the formation of an extended stellar component. Recent works have suggested that galaxies with masses comparable to the Large Magellanic Cloud (LMC, $M_\star \sim 10^9\,{\rm M}_\odot$) may be the primary building blocks of the stellar halo of our Galaxy. We use cosmological simulations of the $Λ$ Cold Dark Matter model to investigate LMC-mass galaxies at $z=1-2$ using a semi-analytic model of galaxy formation. We find that LMC analogues at $z=2$ evolve until the present day along three distinct pathways: (1) those that are destroyed in Milky Way-mass hosts; (2) those that are themselves the main progenitors of Milky Way-mass galaxies; and (3) those that survive until $z=0$, with stellar mass $\sim$1.0 dex lower than typical Milky Ways. We predict that the properties of these galaxies at $z=2$ (stellar metallicities, sizes, gas content etc.) are largely indistinguishable, irrespective of which of these pathways is eventually taken; a survey targeting such galaxies in this redshift range would struggle to tell apart a 'destroyed' stellar halo progenitor from a 'surviving' LMC analogue. The only factor that determines the eventual fate of these galaxies is their proximity to a neighbouring Milky Way main progenitor at $z=2$: while the mean separation to a 'surviving' galaxy is around 7 Mpc, it is only 670 kpc to a 'destroyed' galaxy. This suggests that old stellar populations in the Milky Way may share intrinsic (i.e. non-dynamical) properties that are essentially indistinguishable from progenitors of its stellar halo. △ Less

Submitted 3 March, 2023; v1 submitted 12 September, 2022; originally announced September 2022.

Comments: 9 pages, 5 figures, 1 table. Accepted for publication in MNRAS

Energy wrinkles and phase-space folds of the last major merger

Authors: Vasily Belokurov, Eugene Vasiliev, Alis J. Deason, Sergey E. Koposov, Azadeh Fattahi, Adam M. Dillamore, Elliot Y. Davies, Robert J. J. Grand

Abstract: Relying on the dramatic increase in the number of stars with full 6D phase-space information provided by the Gaia Data Release 3, we discover unambiguous signatures of phase-mixing in the stellar halo around the Sun. We show that for the stars likely belonging to the last massive merger, the (v_r,r) distribution contains a series of long and thin chevron-like overdensities. These phase-space sub-s… ▽ More Relying on the dramatic increase in the number of stars with full 6D phase-space information provided by the Gaia Data Release 3, we discover unambiguous signatures of phase-mixing in the stellar halo around the Sun. We show that for the stars likely belonging to the last massive merger, the (v_r,r) distribution contains a series of long and thin chevron-like overdensities. These phase-space sub-structures are predicted to emerge following the dissolution of a satellite, when its tidal debris is given time to wind up, thin out and fold. Additionally, the observed energy and angular momentum (E, L_z) distribution appears more prograde at high energies, possibly revealing the original orbital angular momentum of the in-falling galaxy. The energy distribution of the debris is strongly asymmetric with a peak at low E -- which, we surmise, may be evidence of the dwarf's rapid sinking -- and riddled with wrinkles and bumps. If these small-scale energy inhomogeneities have been seeded during or immediately after the interaction with the Milky Way, and are not due to the spatial restriction of our study, then making use of the (v_r,r) chevrons to constrain the time of the merger becomes cumbersome. Nonetheless, we demonstrate that similar phase-space and (E,L_z) sub-structures are present in numerical simulations of galaxy interactions, both in bespoke N-body runs and in cosmological hydrodynamical zoom-in suites. The remnant traces of the progenitor's disruption and the signatures of the on-going phase-mixing discovered here will not only help to constrain the properties of our Galaxy's most important interaction, but also can be used as a novel tool to map out the Milky Way's current gravitational potential and its perturbations. △ Less

Submitted 23 August, 2022; originally announced August 2022.

Comments: 16 pages, 14 figures, submitted to MNRAS. Comments are welcome!

The Hubble Space Telescope Survey of M31 Satellite Galaxies I. RR Lyrae-based Distances and Refined 3D Geometric Structure

Authors: Alessandro Savino, Daniel R. Weisz, Evan D. Skillman, Andrew Dolphin, Nitya Kallivayalil, Andrew Wetzel, Jay Anderson, Gurtina Besla, Michael Boylan-Kolchin, James S. Bullock, Andrew A. Cole, Michelle L. M. Collins, M. C. Cooper, Alis J. Deason, Aaron L. Dotter, Mark Fardal, Annette M. N. Ferguson, Tobias K. Fritz, Marla C. Geha, Karoline M. Gilbert, Puragra Guhathakurta, Rodrigo Ibata, Michael J. Irwin, Myoungwon Jeon, Evan Kirby , et al. (11 additional authors not shown)

Abstract: We measure homogeneous distances to M31 and 38 associated stellar systems ($-$16.8$\le M_V \le$ $-$6.0), using time-series observations of RR Lyrae stars taken as part of the Hubble Space Telescope Treasury Survey of M31 Satellites. From $>700$ orbits of new/archival ACS imaging, we identify $>4700$ RR Lyrae stars and determine their periods and mean magnitudes to a typical precision of 0.01 days… ▽ More We measure homogeneous distances to M31 and 38 associated stellar systems ($-$16.8$\le M_V \le$ $-$6.0), using time-series observations of RR Lyrae stars taken as part of the Hubble Space Telescope Treasury Survey of M31 Satellites. From $>700$ orbits of new/archival ACS imaging, we identify $>4700$ RR Lyrae stars and determine their periods and mean magnitudes to a typical precision of 0.01 days and 0.04 mag. Based on Period-Wesenheit-Metallicity relationships consistent with the Gaia eDR3 distance scale, we uniformly measure heliocentric and M31-centric distances to a typical precision of $\sim20$ kpc (3%) and $\sim10$ kpc (8%), respectively. We revise the 3D structure of the M31 galactic ecosystem and: (i) confirm a highly anisotropic spatial distribution such that $\sim80$% of M31's satellites reside on the near side of M31; this feature is not easily explained by observational effects; (ii) affirm the thin (rms $7-23$ kpc) planar "arc" of satellites that comprises roughly half (15) of the galaxies within 300 kpc from M31; (iii) reassess physical proximity of notable associations such as the NGC 147/185 pair and M33/AND XXII; and (iv) illustrate challenges in tip-of-the-red-giant branch distances for galaxies with $M_V > -9.5$, which can be biased by up to 35%. We emphasize the importance of RR Lyrae for accurate distances to faint galaxies that should be discovered by upcoming facilities (e.g., Rubin Observatory). We provide updated luminosities and sizes for our sample. Our distances will serve as the basis for future investigation of the star formation and orbital histories of the entire known M31 satellite system. △ Less

Submitted 12 September, 2022; v1 submitted 6 June, 2022; originally announced June 2022.

Comments: 32 pages, 15 figures, accepted for publication on ApJ

Constraining the shape of dark matter haloes with globular clusters

Authors: Marta Reina-Campos, Sebastian Trujillo-Gomez, Joel L. Pfeffer, Alison Sills, Alis J. Deason, Robert A. Crain, J. M. Diederik Kruijssen

Abstract: 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… ▽ More 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. △ Less

Submitted 25 April, 2022; originally announced April 2022.

Comments: 17 pages, 15 figures; submitted to MNRAS (25/04/2022) - comments and suggestions welcome!

Observing EAGLE galaxies with JWST: predictions for Milky Way progenitors and their building blocks

Authors: Tilly A Evans, Azadeh Fattahi, Alis J Deason, Carlos S Frenk

Abstract: We present predictions, derived from the EAGLE LCDM cosmological hydrodynamical simulations, for the abundance and properties of galaxies expected to be detected at high redshift by the James Webb Space Telescope (JWST). We consider the galaxy population as a whole and focus on the sub-population of progenitors of Milky Way (MW) analogues, defined to be galaxies with accretion histories similar to… ▽ More We present predictions, derived from the EAGLE LCDM cosmological hydrodynamical simulations, for the abundance and properties of galaxies expected to be detected at high redshift by the James Webb Space Telescope (JWST). We consider the galaxy population as a whole and focus on the sub-population of progenitors of Milky Way (MW) analogues, defined to be galaxies with accretion histories similar to the MW's, that is, galaxies that underwent a merger resembling the Gaia-Enceladus-Sausage (GES) event and that contain an analogue of the Large Magellanic Cloud (LMC) satellite today. We derive the luminosity function of all EAGLE galaxies in JWST NIRCam passbands, in the redshift range z=2-8, taking into account dust obscuration and different exposure times. For an exposure time of $T=10^5$s, average MW progenitors are observable as far back as z~6 in most bands, and this changes to z~5 and z~4 for the GES and LMC progenitors, respectively. The progenitors of GES and LMC analogues are, on average, ~2 and ~1 mag fainter than the MW progenitors at most redshifts. They lie, on average, within ~60 and ~30 arcsec, respectively, of their future MW host at all times, and thus will appear within the field-of-view of JWST NIRCam. We conclude that galaxies resembling the main progenitor of the MW and its major accreted components should be observable with JWST beyond redshift 2, providing a new and unique window in studying the formation history of our own galaxy. △ Less

Submitted 5 December, 2022; v1 submitted 4 April, 2022; originally announced April 2022.

Comments: 17 page, 16 figures. Published in MNRAS

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

Authors: Monica Valluri, Solene Chabanier, Vid Irsic, Eric Armengaud, Michael Walther, Connie Rockosi, Miguel A. Sanchez-Conde, Leandro Beraldo e Silva, Andrew P. Cooper, Elise Darragh-Ford, Kyle Dawson, Alis J. Deason, Simone Ferraro, Jaime E. Forero-Romero, Antonella Garzilli, Ting Li, Zarija Lukic, Christopher J. Manser, Nathalie Palanque-Delabrouille, Corentin Ravoux, Ting Tan, Wenting Wang, Risa Wechsler, Andreia Carrillo, Arjun Dey , et al. (7 additional authors not shown)

Abstract: 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 thei… ▽ More 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. △ Less

Submitted 1 July, 2022; v1 submitted 14 March, 2022; originally announced March 2022.

Comments: Contributed white paper to Snowmass 2021, CF03; minor revisions

The Chemo-Dynamical Groups of Galactic Globular Clusters

Authors: Thomas M. Callingham, Marius Cautun, Alis J. Deason, Carlos S. Frenk, Robert J. J Grand, Federico Marinacci

Abstract: 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. O… ▽ More 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. △ Less

Submitted 9 August, 2022; v1 submitted 1 February, 2022; originally announced February 2022.

Comments: 18 pages, 11 figures, + Appendices

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

Authors: Alis J. Deason, Sownak Bose, Azadeh Fattahi, Nicola C. Amorisco, Wojciech Hellwing, Carlos S. Frenk

Abstract: 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 mod… ▽ More 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. △ Less

Submitted 1 January, 2022; v1 submitted 11 October, 2021; originally announced October 2021.

Comments: 16 pages, 12 figures. Accepted by MNRAS

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

Authors: Marta Reina-Campos, Sebastian Trujillo-Gomez, Alis J. Deason, J. M. Diederik Kruijssen, Joel L. Pfeffer, Robert A. Crain, Nate Bastian, Meghan E. Hughes

Abstract: 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… ▽ More 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. △ Less

Submitted 14 June, 2021; originally announced June 2021.

Comments: 19 pages, 11 figures and 2 tables; submitted to MNRAS, comments and/or suggestions are welcomed!

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

Authors: Sandro Tacchella, Charlie Conroy, S. M. Faber, Benjamin D. Johnson, Joel Leja, Guillermo Barro, Emily C. Cunningham, Alis J. Deason, Puragra Guhathakurta, Yicheng Guo, Lars Hernquist, David C. Koo, Kevin McKinnon, Constance M. Rockosi, Joshua S. Speagle, Pieter van Dokkum, Hassen M. Yesuf

Abstract: 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, emissio… ▽ More 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. △ Less

Submitted 5 January, 2022; v1 submitted 24 February, 2021; originally announced February 2021.

Comments: 30 pages, 12 figures + appendix, accepted by ApJ

The survival of globular clusters in a cuspy Fornax

Authors: Shi Shao, Marius Cautun, Carlos S. Frenk, Marta Reina-Campos, Alis J. Deason, Robert A. Crain, J. M. Diederik Kruijssen, Joel Pfeffer

Abstract: It has long been argued that the radial distribution of globular clusters (GCs) in the Fornax dwarf galaxy requires its dark matter halo to have a core of size $\sim 1$ kpc. We revisit this argument by investigating analogues of Fornax formed in E-MOSAICS, a cosmological hydrodynamical simulation that self-consistently follows the formation and evolution of GCs in the EAGLE galaxy formation model.… ▽ More It has long been argued that the radial distribution of globular clusters (GCs) in the Fornax dwarf galaxy requires its dark matter halo to have a core of size $\sim 1$ kpc. We revisit this argument by investigating analogues of Fornax formed in E-MOSAICS, a cosmological hydrodynamical simulation that self-consistently follows the formation and evolution of GCs in the EAGLE galaxy formation model. In EAGLE, Fornax-mass haloes are cuspy and well described by the Navarro-Frenk-White profile. We post-process the E-MOSAICS to account for GC orbital decay by dynamical friction, which is not included in the original model. Dynamical friction causes 33 per cent of GCs with masses $M_{\rm GC}\geq4\times10^4 {~\rm M_\odot}$ to sink to the centre of their host where they are tidally disrupted. Fornax has a total of five GCs, an exceptionally large number compared to other galaxies of similar stellar mass. In the simulations, we find that only 3 per cent of the Fornax analogues have five or more GCs, while 30 per cent have only one and 35 per cent have none. We find that GC systems in satellites are more centrally concentrated than in field dwarfs, and that those formed in situ (45 per cent) are more concentrated than those that were accreted. The survival probability of a GC increases rapidly with the radial distance at which it formed ($r_{\rm init}$): it is 37 per cent for GCs with $r_{\rm init} \leq 1$ kpc and 92 per cent for GCs with $r_{\rm init} \geq 1$ kpc. The present-day radial distribution of GCs in E-MOSAICS turns out to be indistinguishable from that in Fornax, demonstrating that, contrary to claims in the literature, the presence of five GCs in the central kiloparsec of Fornax does not exclude a cuspy DM halo. △ Less

Submitted 14 December, 2020; originally announced December 2020.

Comments: 15 pages, 13 figures. submitted to MNRAS, comments welcome!

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

Authors: Alis J. Deason, Denis Erkal, Vasily Belokurov, Azadeh Fattahi, Facundo A. Gómez, Robert J. J. Grand, Rüdiger Pakmor, Xiang-Xiang Xue, Chao Liu, Chengqun Yang, Lan Zhang, Gang Zhao

Abstract: We use a distribution function analysis to estimate the mass of the Milky Way out to 100 kpc using a large sample of halo stars. These stars are compiled from the literature, and the vast majority (~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. T… ▽ More 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. △ Less

Submitted 18 January, 2021; v1 submitted 26 October, 2020; originally announced October 2020.

Comments: 9 pages, 6 figures. MNRAS accepted

Detection of the LMC-induced sloshing of the Galactic halo

Authors: Denis Erkal, Alis J. Deason, Vasily Belokurov, Xiang-Xiang Xue, Sergey E. Koposov, Sarah A. Bird, Chao Liu, Iulia T. Simion, Chengqun Yang, Lan Zhang, Gang Zhao

Abstract: A wealth of recent studies have shown that the LMC is likely massive, with a halo mass $>10^{11} M_\odot$. One consequence of having such a nearby and massive neighbour is that the inner Milky Way is expected to be accelerated with respect to our Galaxy's outskirts (beyond $\sim 30$ kpc). In this work we compile a sample of $\sim 500$ stars with radial velocities in the distant stellar halo,… ▽ More A wealth of recent studies have shown that the LMC is likely massive, with a halo mass $>10^{11} M_\odot$. One consequence of having such a nearby and massive neighbour is that the inner Milky Way is expected to be accelerated with respect to our Galaxy's outskirts (beyond $\sim 30$ kpc). In this work we compile a sample of $\sim 500$ stars with radial velocities in the distant stellar halo, $r_{\rm GC}> 50$ kpc, to test this hypothesis. These stars span a large fraction of the sky and thus give a global view of the stellar halo. We find that stars in the Southern hemisphere are on average blueshifted, while stars in the North are redshifted, consistent with the expected, mostly downwards acceleration of the inner halo due to the LMC. We compare these results with simulations and find the signal is consistent with the infall of a $1.5\times10^{11} M_\odot$ LMC. We cross-match our stellar sample with \textit{Gaia} DR2 and find that the mean proper motions are not yet precise enough to discern the LMC's effect. Our results show that the outer Milky Way is significantly out of equilibrium and that the LMC has a substantial effect on our Galaxy. △ Less

Submitted 26 October, 2020; originally announced October 2020.

Comments: 9 pages, 7 figures, submitted to MNRAS. Comments welcome!

Stellar splashback: the edge of the intracluster light

Authors: Alis J. Deason, Kyle A. Oman, Azadeh Fattahi, Matthieu Schaller, Mathilde Jauzac, Yuanyuan Zhang, Mireia Montes, Yannick M. Bahé, Claudio Dalla Vecchia, Scott T. Kay, Tilly A. Evans

Abstract: We examine the outskirts of galaxy clusters in the C-EAGLE simulations to quantify the `edges' of the stellar and dark matter distribution. The radius of the steepest slope in the dark matter, commonly used as a proxy for the splashback radius, is located at ~r_200m; the strength and location of this feature depends on the recent mass accretion rate, in good agreement with previous work. Interesti… ▽ More We examine the outskirts of galaxy clusters in the C-EAGLE simulations to quantify the `edges' of the stellar and dark matter distribution. The radius of the steepest slope in the dark matter, commonly used as a proxy for the splashback radius, is located at ~r_200m; the strength and location of this feature depends on the recent mass accretion rate, in good agreement with previous work. Interestingly, the stellar distribution (or intracluster light, ICL) also has a well-defined edge, which is directly related to the splashback radius of the halo. Thus, detecting the edge of the ICL can provide an independent measure of the physical boundary of the halo, and the recent mass accretion rate. We show that these caustics can also be seen in the projected density profiles, but care must be taken to account for the influence of substructures and other non-diffuse material, which can bias and/or weaken the signal of the steepest slope. This is particularly important for the stellar material, which has a higher fraction bound in subhaloes than the dark matter. Finally, we show that the `stellar splashback' feature is located beyond current observational constraints on the ICL, but these large projected distances (>> 1 Mpc) and low surface brightnesses (mu >> 32 mag/arcsec^2) can be reached with upcoming observational facilities such as the Vera C. Rubin Observatory, the Nancy Grace Roman Space Telescope, and Euclid. △ Less

Submitted 29 November, 2020; v1 submitted 6 October, 2020; originally announced October 2020.

Comments: 12 pages, 9 figures. Accepted for publication in MNRAS

Space Telescope and Optical Reverberation Mapping Project. XII. Broad-Line Region Modeling of NGC 5548

Authors: P. R. Williams, A. Pancoast, T. Treu, B. J. Brewer, B. M. Peterson, A. J. Barth, M. A. Malkan, G. De Rosa, Keith Horne, G. A. Kriss, N. Arav, M. C. Bentz, E. M. Cackett, E. Dalla Bontà, M. Dehghanian, C. Done, G. J. Ferland, C. J. Grier, J. Kaastra, E. Kara, C. S. Kochanek, S. Mathur, M. Mehdipour, R. W. Pogge, D. Proga , et al. (133 additional authors not shown)

Abstract: We present geometric and dynamical modeling of the broad line region for the multi-wavelength reverberation mapping campaign focused on NGC 5548 in 2014. The dataset includes photometric and spectroscopic monitoring in the optical and ultraviolet, covering the H$β$, C IV, and Ly$α$ broad emission lines. We find an extended disk-like H$β$ BLR with a mixture of near-circular and outflowing gas traje… ▽ More We present geometric and dynamical modeling of the broad line region for the multi-wavelength reverberation mapping campaign focused on NGC 5548 in 2014. The dataset includes photometric and spectroscopic monitoring in the optical and ultraviolet, covering the H$β$, C IV, and Ly$α$ broad emission lines. We find an extended disk-like H$β$ BLR with a mixture of near-circular and outflowing gas trajectories, while the C IV and Ly$α$ BLRs are much less extended and resemble shell-like structures. There is clear radial structure in the BLR, with C IV and Ly$α$ emission arising at smaller radii than the H$β$ emission. Using the three lines, we make three independent black hole mass measurements, all of which are consistent. Combining these results gives a joint inference of $\log_{10}(M_{\rm BH}/M_\odot) = 7.64^{+0.21}_{-0.18}$. We examine the effect of using the $V$ band instead of the UV continuum light curve on the results and find a size difference that is consistent with the measured UV-optical time lag, but the other structural and kinematic parameters remain unchanged, suggesting that the $V$ band is a suitable proxy for the ionizing continuum when exploring the BLR structure and kinematics. Finally, we compare the H$β$ results to similar models of data obtained in 2008 when the AGN was at a lower luminosity state. We find that the size of the emitting region increased during this time period, but the geometry and black hole mass remain unchanged, which confirms that the BLR kinematics suitably gauge the gravitational field of the central black hole. △ Less

Submitted 1 October, 2020; originally announced October 2020.

Comments: 26 pages, 19 figures, 1 table, accepted for publication in ApJ

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

Authors: Emily C. Cunningham, Nicolas Garavito-Camargo, Alis J. Deason, Kathryn V. Johnston, Denis Erkal, Chervin F. P. Laporte, Gurtina Besla, Rodrigo Luger, Robyn E. Sanderson

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. (2019) (hereafter G… ▽ More 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. △ Less

Submitted 15 June, 2020; originally announced June 2020.

Comments: 25 pages, 13 figures. Accepted for publication in ApJ

How Unusual is the Milky Way's Assembly History?

Authors: Tilly A. Evans, Azadeh Fattahi, Alis J. Deason, Carlos S. Frenk

Abstract: 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 earl… ▽ More 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. △ Less

Submitted 11 May, 2020; originally announced May 2020.

Comments: 12 pages, 8 figures, submitted to MNRAS, contact Tilly Evans with comments or questions

The twisted dark matter halo of the Milky Way

Authors: Shi Shao, Marius Cautun, Alis J. Deason, Carlos S. Frenk

Abstract: 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… ▽ More 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. △ Less

Submitted 15 December, 2020; v1 submitted 6 May, 2020; originally announced May 2020.

Comments: 16 pages, 13 figures. Accepted for publication in MNRAS

Space Telescope and Optical Reverberation Mapping Project. IX. Velocity-Delay Maps for Broad Emission Lines in NGC 5548

Authors: Keith Horne, G. De Rosa, B. M. Peterson, A. J. Barth, J. Ely, M. M. Fausnaugh, G. A. Kriss, L. Pei, S. M. Adams, M. D. Anderson, P. Arevalo, T G. Beatty, V. N. Bennert, M. C. Bentz, A. Bigley, S. Bisogni, G. A. Borman, T. A. Boroson, M. C. Bottorff, W. N. Brandt, A. A. Breeveld, M. Brotherton, J. E. Brown, J. S. Brown, E. M. Cackett , et al. (133 additional authors not shown)

Abstract: We report velocity-delay maps for prominent broad emission lines, Ly_alpha, CIV, HeII and H_beta, in the spectrum of NGC5548. The emission-line responses inhabit the interior of a virial envelope. The velocity-delay maps reveal stratified ionization structure. The HeII response inside 5-10 light-days has a broad single-peaked velocity profile. The Ly_alpha, CIV, and H_beta responses peak inside 10… ▽ More We report velocity-delay maps for prominent broad emission lines, Ly_alpha, CIV, HeII and H_beta, in the spectrum of NGC5548. The emission-line responses inhabit the interior of a virial envelope. The velocity-delay maps reveal stratified ionization structure. The HeII response inside 5-10 light-days has a broad single-peaked velocity profile. The Ly_alpha, CIV, and H_beta responses peak inside 10 light-days, extend outside 20 light-days, and exhibit a velocity profile with two peaks separated by 5000 km/s in the 10 to 20 light-day delay range. The velocity-delay maps show that the M-shaped lag vs velocity structure found in previous cross-correlation analysis is the signature of a Keplerian disk with a well-defined outer edge at R=20 light-days. The outer wings of the M arise from the virial envelope, and the U-shaped interior of the M is the lower half of an ellipse in the velocity-delay plane. The far-side response is weaker than that from the near side, so that we see clearly the lower half, but only faintly the upper half, of the velocity--delay ellipse. The delay tau=(R/c)(1-sin(i))=5 light-days at line center is from the near edge of the inclined ring, giving the inclination i=45 deg. A black hole mass of M=7x10^7 Msun is consistent with the velocity-delay structure. A barber-pole pattern with stripes moving from red to blue across the CIV and possibly Ly_alpha line profiles suggests the presence of azimuthal structure rotating around the far side of the broad-line region and may be the signature of precession or orbital motion of structures in the inner disk. Further HST observations of NGC 5548 over a multi-year timespan but with a cadence of perhaps 10 days rather than 1 day could help to clarify the nature of this new AGN phenomenon. △ Less

Submitted 27 November, 2020; v1 submitted 3 March, 2020; originally announced March 2020.

Comments: 19 pages, 9 figures, ApJ in press

A tale of two populations: surviving and destroyed dwarf galaxies and the build up of the Milky Way's stellar halo

Authors: Azadeh Fattahi, Alis J. Deason, Carlos S. Frenk, Christine M. Simpson, Facundo A. Gomez, Robert J. J. Grand, Antonela Monachesi, Federico Marinacci, Ruediger Pakmor

Abstract: We use magneto-hydrodynamical simulations of Milky Way-mass haloes from the Auriga project to examine the properties of surviving and destroyed dwarf galaxies that are accreted by these haloes over cosmic time. We show that the combined luminosity function of surviving and destroyed dwarfs at infall is similar in the various Auriga haloes, and is dominated by the destroyed dwarfs. There is, howeve… ▽ More We use magneto-hydrodynamical simulations of Milky Way-mass haloes from the Auriga project to examine the properties of surviving and destroyed dwarf galaxies that are accreted by these haloes over cosmic time. We show that the combined luminosity function of surviving and destroyed dwarfs at infall is similar in the various Auriga haloes, and is dominated by the destroyed dwarfs. There is, however, a strong dependence on infall time: destroyed dwarfs have typically early infall times, $t_{infall}<6$ Gyr, whereas the majority of dwarfs accreted at $t_{infall}>10$ Gyr have survived to the present day. Because of their late infall the surviving satellites today had higher metallicites at infall than their destroyed counterparts of similar infall mass; the difference is even more pronounced for the present-day metallicites of satellites, many of which continue to form stars after infall. In agreement with previous work, we find that a small number of relatively massive destroyed dwarf galaxies dominate the mass of the stellar haloes. However, there is a significant radial dependence: while 90 per cent of the mass in the inner regions ($<\,20\,$kpc) is contributed, on average, by only 3 massive progenitors, the outer regions ($>\,100\,$kpc) typically have $\sim8$ main progenitors of relatively lower mass. Finally, we show that a few massive progenitors dominate the metallicity distribution of accreted stars, even at the metal poor end. Contrary to common assumptions in the literature, dwarf galaxies of mass $M_{*}<10^7 \, M_{\odot}$ make up less than 10 per cent of the accreted, metal poor stars ([Fe/H] $<\,-3$) in the inner $50\,$kpc. △ Less

Submitted 7 October, 2020; v1 submitted 27 February, 2020; originally announced February 2020.

Comments: Replaced with the accepted version after minor modifications

The Edge of the Galaxy

Authors: Alis J. Deason, Azadeh Fattahi, Carlos S. Frenk, Robert J. J. Grand, Kyle A. Oman, Shea Garrison-Kimmel, Christine M. Simpson, Julio F. Navarro

Abstract: 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 s… ▽ More 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. △ Less

Submitted 25 June, 2020; v1 submitted 21 February, 2020; originally announced February 2020.

Comments: 14 pages, 11 figures. MNRAS in press

The orbital phase space of contracted dark matter halos

Authors: Thomas M. Callingham, Marius Cautun, Alis J. Deason, Carlos S. Frenk, Robert J. J. Grand, Federico Marinacci, Ruediger Pakmor

Abstract: 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… ▽ More 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. △ Less

Submitted 21 January, 2020; originally announced January 2020.

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

Authors: Robert J. J. Grand, Daisuke Kawata, Vasily Belokurov, Alis J. Deason, Azadeh Fattahi, Francesca Fragkoudi, Facundo A. Gómez, Federico Marinacci, Rüdiger Pakmor

Abstract: 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 th… ▽ More 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. △ Less

Submitted 13 July, 2020; v1 submitted 16 January, 2020; originally announced January 2020.

Comments: MNRAS accepted

The Milky Way total mass profile as inferred from Gaia DR2

Authors: Marius Cautun, Alejandro Benitez-Llambay, Alis J. Deason, Carlos S. Frenk, Azadeh Fattahi, Facundo A. Gómez, Robert J. J. Grand, Kyle A. Oman, Julio F. Navarro, Christine M. Simpson

Abstract: 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 d… ▽ More 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. △ Less

Submitted 16 April, 2020; v1 submitted 11 November, 2019; originally announced November 2019.

Comments: 21 pages, 14 figures, main results: figs 3, 7, 10 and 11; MNRAS accepted version, code and resources available at https://github.com/MariusCautun/Milky_Way_mass_profile

The Biggest Splash

Authors: Vasily Belokurov, Jason L. Sanders, Azadeh Fattahi, Martin C. Smith, Alis J. Deason, N. Wyn Evans, Robert J. J. Grand

Abstract: 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 element… ▽ More 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. △ Less

Submitted 10 September, 2019; originally announced September 2019.

Comments: Submitted to MNRAS. Comments are welcome!

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

Authors: Sownak Bose, Alis J. Deason, Vasily Belokurov, Carlos S. Frenk

Abstract: 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… ▽ More 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. △ Less

Submitted 9 September, 2019; originally announced September 2019.

Comments: 16 pages, 11 figures, 1 table, comments welcome

The total stellar halo mass of the Milky Way

Authors: Alis J. Deason, Vasily Belokurov, Jason L. Sanders

Abstract: We measure the total stellar halo luminosity using red giant branch (RGB) stars selected from Gaia data release 2. Using slices in magnitude, colour and location on the sky, we decompose RGB stars belonging to the disc and halo by fitting 2-dimensional Gaussians to the Galactic proper motion distributions. The number counts of RGB stars are converted to total stellar halo luminosity using a suite… ▽ More We measure the total stellar halo luminosity using red giant branch (RGB) stars selected from Gaia data release 2. Using slices in magnitude, colour and location on the sky, we decompose RGB stars belonging to the disc and halo by fitting 2-dimensional Gaussians to the Galactic proper motion distributions. The number counts of RGB stars are converted to total stellar halo luminosity using a suite of isochrones weighted by age and metallicity, and by applying a volume correction based on the stellar halo density profile. Our method is tested and calibrated using Galaxia and N-body models. We find a total luminosity (out to 100 kpc) of L_halo = 7.9 +/- 2.0 x 10^8 L_Sun excluding Sgr, and L_halo = 9.4 +/- 2.4 x 10^8 L_Sun including Sgr. These values are appropriate for our adopted stellar halo density profile and metallicity distribution, but additional systematics related to these assumptions are quantified and discussed. Assuming a stellar mass-to-light ratio appropriate for a Kroupa initial mass function (M*/L = 1.5), we estimate a stellar halo mass of M*_halo = 1.4 +/- 0.4 x 10^9 M_Sun. This mass is larger than previous estimates in the literature, but is in good agreement with the emerging picture that the (inner) stellar halo is dominated by one massive dwarf progenitor. Finally, we argue that the combination of a ~10^9 M_Sun mass and an average metallicity of <[Fe/H]> ~ -1.5 for the Galactic halo points to an ancient (~10 Gyr) merger event. △ Less

Submitted 2 October, 2019; v1 submitted 7 August, 2019; originally announced August 2019.

Comments: 14 pages, 13 figures. Accepted for publication in MNRAS

Galactic Gas Flows from Halo to Disk: Tomography and Kinematics at the Milky Way's Disk-Halo Interface

Authors: Hannah V. Bish, Jessica K. Werk, J. Xavier Prochaska, Kate H. R. Rubin, Yong Zheng, John M. O'Meara, Alis J. Deason

Abstract: We present a novel absorption line survey using 54 blue horizontal branch stars (BHBs) in the Milky Way halo as background sources for detecting gas flows at the disk-halo interface. Distance measurements to high-latitude ($b$ > 60°) background stars at 3.1-13.4 kpc, combined with unprecedented spatial sampling and spectral resolution, allow us to examine the 3-dimensional spatial distribution and… ▽ More We present a novel absorption line survey using 54 blue horizontal branch stars (BHBs) in the Milky Way halo as background sources for detecting gas flows at the disk-halo interface. Distance measurements to high-latitude ($b$ > 60°) background stars at 3.1-13.4 kpc, combined with unprecedented spatial sampling and spectral resolution, allow us to examine the 3-dimensional spatial distribution and kinematics of gas flows near the disk. We detect absorption signatures of extraplanar CaII and NaI in Keck HIRES spectra and find that their column densities exhibit no trend with distance to the background sources, indicating that these clouds lie within 3.1 kpc of the disk. We calculate covering fractions of $f_{CaII}$ = 63%, $f_{NaI}$ = 26%, and $f_{HI}$ = 52%, consistent with a picture of the CGM that includes multi-phase clouds containing small clumps of cool gas within hotter, more diffuse gas. Our measurements constrain the scale of any substructure within these cool clouds to <0.5 kpc. CaII and NaI absorption features exhibit an intermediate-velocity (IV) component inflowing at velocities of -75 km/s < v < -25 km/s relative to the local standard of rest, consistent with previously-studied HI structures in this region. We report the new detection of an inflow velocity gradient $Δv_z$ ~ 6-9 km/s/kpc across the Galactic plane. These findings place constraints on the physical and kinematic properties of CGM gas flows through the disk-halo interface, and support a galactic fountain model in which cold gas rains back onto the disk. △ Less

Submitted 22 July, 2019; originally announced July 2019.

Comments: 18 pages, 15 figures, 2 tables, and a 10-page appendix with 9 additional figures. Submitted to ApJ

Space Telescope and Optical Reverberation Mapping Project. VIII. Time Variability of Emission and Absorption in NGC 5548 Based on Modeling the Ultraviolet Spectrum

Authors: G. A. Kriss, G. De Rosa, J. Ely, B. M. Peterson, J. Kaastra, M. Mehdipour, G. J. Ferland, M. Dehghanian, S. Mathur, R. Edelson, K. T. Korista, N. Arav, A. J. Barth, M. C. Bentz, W. N. Brandt, D. M. Crenshaw, E. Dalla Bontà, K. D. Denney, C. Done, M. Eracleous, M. M. Fausnaugh, E. Gardner, M. R. Goad, C. J. Grier, Keith Horne , et al. (142 additional authors not shown)

Abstract: We model the ultraviolet spectra of the Seyfert 1 galaxy NGC~5548 obtained with the Hubble Space Telescope during the 6-month reverberation-mapping campaign in 2014. Our model of the emission from NGC 5548 corrects for overlying absorption and deblends the individual emission lines. Using the modeled spectra, we measure the response to continuum variations for the deblended and absorption-correcte… ▽ More We model the ultraviolet spectra of the Seyfert 1 galaxy NGC~5548 obtained with the Hubble Space Telescope during the 6-month reverberation-mapping campaign in 2014. Our model of the emission from NGC 5548 corrects for overlying absorption and deblends the individual emission lines. Using the modeled spectra, we measure the response to continuum variations for the deblended and absorption-corrected individual broad emission lines, the velocity-dependent profiles of Ly$α$ and C IV, and the narrow and broad intrinsic absorption features. We find that the time lags for the corrected emission lines are comparable to those for the original data. The velocity-binned lag profiles of Ly$α$ and C IV have a double-peaked structure indicative of a truncated Keplerian disk. The narrow absorption lines show delayed response to continuum variations corresponding to recombination in gas with a density of $\sim 10^5~\rm cm^{-3}$. The high-ionization narrow absorption lines decorrelate from continuum variations during the same period as the broad emission lines. Analyzing the response of these absorption lines during this period shows that the ionizing flux is diminished in strength relative to the far-ultraviolet continuum. The broad absorption lines associated with the X-ray obscurer decrease in strength during this same time interval. The appearance of X-ray obscuration in $\sim\,2012$ corresponds with an increase in the luminosity of NGC 5548 following an extended low state. We suggest that the obscurer is a disk wind triggered by the brightening of NGC 5548 following the decrease in size of the broad-line region during the preceding low-luminosity state. △ Less

Submitted 12 July, 2019; v1 submitted 8 July, 2019; originally announced July 2019.

Comments: 50 pages, 30 figures, uses aastex62.cls. Accepted for publication in ApJ, 07/06/2019. High-level products page in MAST will go live after 7/15/2019. Replaced Figure 4 on 7/12/2019 to be more red/green color-blind friendly

Simulating cosmological substructure in the solar neighbourhood

Authors: Christine M. Simpson, Ignacio Gargiulo, Facundo A. Gómez, Robert J. J. Grand, Nicolás Maffione, Andrew P. Cooper, Alis J. Deason, Carlos Frenk, John Helly, Federico Marinacci, Rüdiger Pakmor

Abstract: 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 str… ▽ More 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. △ Less

Submitted 12 September, 2019; v1 submitted 23 May, 2019; originally announced May 2019.

Comments: 6 pages, 4 figures, accepted for publication in MNRAS letters

The effects of dynamical substructure on Milky Way mass estimates from the high velocity tail of the local stellar halo

Authors: Robert J. J. Grand, Alis J. Deason, Simon D. M. White, Christine M. Simpson, Facundo A. Gómez, Federico Marinacci, Ruediger Pakmor

Abstract: We investigate the impact of dynamical streams and substructure on estimates of the local escape speed and total mass of Milky Way-mass galaxies from modelling the high velocity tail of local halo stars. We use a suite of high-resolution, magneto-hydrodynamical cosmological zoom-in simulations, which resolve phase space substructure in local volumes around solar-like positions. We show that phase… ▽ More We investigate the impact of dynamical streams and substructure on estimates of the local escape speed and total mass of Milky Way-mass galaxies from modelling the high velocity tail of local halo stars. We use a suite of high-resolution, magneto-hydrodynamical cosmological zoom-in simulations, which resolve phase space substructure in local volumes around solar-like positions. We show that phase space structure varies significantly between positions in individual galaxies and across the suite. Substructure populates the high velocity tail unevenly and leads to discrepancies in the mass estimates. We show that a combination of streams, sample noise and truncation of the high velocity tail below the escape speed leads to a distribution of mass estimates with a median that falls below the true value by $\sim 20 \%$, and a spread of a factor of 2 across the suite. Correcting for these biases, we derive a revised value for the Milky Way mass presented in Deason et al. of $1.29 ^{+0.37}_{-0.47} \times 10^{12}$ $\rm M_{\odot}$. △ Less

Submitted 3 June, 2019; v1 submitted 23 May, 2019; originally announced May 2019.

Comments: Accepted to the Monthly Notices of the Royal Astronomical Society, Letters, with minor revision

The Nature of Ionized Gas in the Milky Way Galactic Fountain

Authors: Jessica K. Werk, Kate H. R. Rubin, Hannah V. Bish, J. X. Prochaska, Yong Zheng, John M. O'Meara, Daniel Lenz, Cameron Hummels, Alis J. Deason

Abstract: 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… ▽ More 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. △ Less

Submitted 18 November, 2019; v1 submitted 24 April, 2019; originally announced April 2019.

Comments: 12 figures, 4 tables, Accepted to ApJ. Revision features areal mass accretion rate calculation

The Pisces Plume and the Magellanic Wake

Authors: Vasily Belokurov, Alis J. Deason, Denis Erkal, Sergey E. Koposov, Julio A. Carballo-Bello, Martin C Smith, Prashin Jethwa, Camila Navarrete

Abstract: Using RR Lyrae stars in the Gaia Data Release 2 and Pan-STARRS1 we study the properties of the Pisces Overdensity, a diffuse sub-structure in the outer halo of the Milky Way. We show that along the line of sight, Pisces appears as a broad and long plume of stars stretching from 40 to 110 kpc with a steep distance gradient. On the sky Pisces's elongated shape is aligned with the Magellanic Stream.… ▽ More Using RR Lyrae stars in the Gaia Data Release 2 and Pan-STARRS1 we study the properties of the Pisces Overdensity, a diffuse sub-structure in the outer halo of the Milky Way. We show that along the line of sight, Pisces appears as a broad and long plume of stars stretching from 40 to 110 kpc with a steep distance gradient. On the sky Pisces's elongated shape is aligned with the Magellanic Stream. Using follow-up VLT FORS2 spectroscopy, we have measured the velocity distribution of the Pisces candidate member stars and have shown it to be as broad as that of the Galactic halo but offset to negative velocities. Using a suite of numerical simulations, we demonstrate that the structure has many properties in common with the predicted behaviour of the Magellanic wake, i.e. the Galactic halo overdensity induced by the in-fall of the Magellanic Clouds. △ Less

Submitted 16 April, 2019; originally announced April 2019.

Comments: Submitted to MNRAS Letters. Comments welcome!

The local high velocity tail and the Galactic escape speed

Authors: Alis J. Deason, Azadeh Fattahi, Vasily Belokurov, Wyn Evans, Robert J. Grand, Federico Marinacci, Rudiger Pakmor

Abstract: We model the fastest moving (v_tot > 300 km/s) local (D < 3 kpc) halo stars using cosmological simulations and 6-dimensional Gaia data. Our approach is to use our knowledge of the assembly history and phase-space distribution of halo stars to constrain the form of the high velocity tail of the stellar halo. Using simple analytical models and cosmological simulations, we find that the shape of the… ▽ More We model the fastest moving (v_tot > 300 km/s) local (D < 3 kpc) halo stars using cosmological simulations and 6-dimensional Gaia data. Our approach is to use our knowledge of the assembly history and phase-space distribution of halo stars to constrain the form of the high velocity tail of the stellar halo. Using simple analytical models and cosmological simulations, we find that the shape of the high velocity tail is strongly dependent on the velocity anisotropy and number density profile of the halo stars --- highly eccentric orbits and/or shallow density profiles have more extended high velocity tails. The halo stars in the solar vicinity are known to have a strongly radial velocity anisotropy, and it has recently been shown the origin of these highly eccentric orbits is the early accretion of a massive (M_star ~ 10^9 M_Sun) dwarf satellite. We use this knowledge to construct a prior on the shape of the high velocity tail. Moreover, we use the simulations to define an appropriate outer boundary of 2r_200, beyond which stars can escape. After applying our methodology to the Gaia data, we find a local (r_0=8.3 kpc) escape speed of v_esc(r_0) = 528(+24,-25) km/s. We use our measurement of the escape velocity to estimate the total Milky Way mass, and dark halo concentration: M_200,tot = 1.00(+0.31,-0.24) x 10^12 M_Sun, c_200 = 10.9(+4.4,-3.3). Our estimated mass agrees with recent results in the literature that seem to be converging on a Milky Way mass of M_200,tot ~ 10^12 M_Sun. △ Less

Submitted 28 February, 2019; v1 submitted 7 January, 2019; originally announced January 2019.

Comments: 13 pages, 10 Figures. MNRAS accepted

HALO7D II: The Halo Velocity Ellipsoid and Velocity Anisotropy with Distant Main Sequence Stars

Authors: Emily C. Cunningham, Alis J. Deason, Robyn E. Sanderson, Sangmo Tony Sohn, Jay Anderson, Puragra Guhathakurta, Constance M. Rockosi, Roeland P. van der Marel, Sarah R. Loebman, Andrew Wetzel

Abstract: The Halo Assembly in Lambda-CDM: Observations in 7 Dimensions (HALO7D) dataset consists of Keck II/DEIMOS spectroscopy and Hubble Space Telescope-measured proper motions of Milky Way (MW) halo main sequence turnoff stars in the CANDELS fields. In this paper, the second in the HALO7D series, we present the proper motions for the HALO7D sample. We discuss our measurement methodology, which makes use… ▽ More The Halo Assembly in Lambda-CDM: Observations in 7 Dimensions (HALO7D) dataset consists of Keck II/DEIMOS spectroscopy and Hubble Space Telescope-measured proper motions of Milky Way (MW) halo main sequence turnoff stars in the CANDELS fields. In this paper, the second in the HALO7D series, we present the proper motions for the HALO7D sample. We discuss our measurement methodology, which makes use of a Bayesian mixture modeling approach for creating the stationary reference frame of distant galaxies. Using the 3D kinematic HALO7D sample, we estimate the parameters of the halo velocity ellipsoid, $\langle v_φ \rangle, σ_r, σ_φ, σ_θ$, and the velocity anisotropy $β$. Using the full HALO7D sample, we find $β=0.63 \pm 0.05$ at $\langle r \rangle =24$ kpc. We also estimate the ellipsoid parameters for our sample split into three apparent magnitude bins; the posterior medians for these estimates of $β$, while consistent with one another, increase as a function of mean sample distance. Finally, we estimate $β$ in each of the individual HALO7D fields. We find that the velocity anisotropy $β$ can vary from field to field, which suggests that the halo is not phase mixed at $\langle r \rangle =24$ kpc. We explore the $β$ variation across the skies of two stellar halos from the \textit{Latte} suite of FIRE-2 simulations, finding that both simulated galaxies show $β$ variation over a similar range to the variation observed across the four HALO7D fields. The accretion histories of the two simulated galaxies result in different $β$ variation patterns; spatially mapping $β$ is thus a way forward in characterizing the accretion history of the Galaxy. △ Less

Submitted 29 October, 2018; originally announced October 2018.

Comments: 23 pages, 13 figures. Submitted to ApJ; comments welcome

The origin of galactic metal-rich stellar halo components with highly eccentric orbits

Authors: Azadeh Fattahi, Vasily Belokurov, Alis J. Deason, Carlos S. Frenk, Facundo A. Gomez, Robert J. J. Grand, Federico Marinacci, Rudiger Pakmor, Volker Springel

Abstract: Using the astrometry from the ESA's Gaia mission, previous works have shown that the Milky Way stellar halo is dominated by metal-rich stars on highly eccentric orbits. To shed light on the nature of this prominent halo component, we have analysed 28 Galaxy analogues in the Auriga suite of cosmological hydrodynamics zoom-in simulations. Some three quarters of the Auriga galaxies contain significan… ▽ More Using the astrometry from the ESA's Gaia mission, previous works have shown that the Milky Way stellar halo is dominated by metal-rich stars on highly eccentric orbits. To shed light on the nature of this prominent halo component, we have analysed 28 Galaxy analogues in the Auriga suite of cosmological hydrodynamics zoom-in simulations. Some three quarters of the Auriga galaxies contain significant components with high radial velocity anisotropy, beta > 0.6. However, only in one third of the hosts do the high-beta stars contribute significantly to the accreted stellar halo overall, similar to what is observed in the Milky Way. For this particular subset we reveal the origin of the dominant stellar halo component with high metallicity, [Fe/H]~-1, and high orbital anisotropy, beta>0.8, by tracing their stars back to the epoch of accretion. It appears that, typically, these stars come from a single dwarf galaxy with a stellar mass of order of 10^9-10^10 Msol that merged around 6-10 Gyr ago, causing a sharp increase in the halo mass. Our study therefore establishes a firm link between the excess of radially anisotropic stellar debris in the Milky Way halo and an ancient head-on collision between the young Milky Way and a massive dwarf galaxy △ Less

Submitted 7 October, 2020; v1 submitted 17 October, 2018; originally announced October 2018.

Comments: replaced with the accepted version (very minor changes)

The aftermath of the Great Collision between our Galaxy and the Large Magellanic Cloud

Authors: Marius Cautun, Alis J. Deason, Carlos S. Frenk, Stuart McAlpine

Abstract: The Milky Way (MW) offers a uniquely detailed view of galactic structure and is often regarded as a prototypical spiral galaxy. But recent observations indicate that the MW is atypical: it has an undersized supermassive black hole at its centre; it is surrounded by a very low mass, excessively metal-poor stellar halo; and it has an unusually large nearby satellite galaxy, the Large Magellanic Clou… ▽ More The Milky Way (MW) offers a uniquely detailed view of galactic structure and is often regarded as a prototypical spiral galaxy. But recent observations indicate that the MW is atypical: it has an undersized supermassive black hole at its centre; it is surrounded by a very low mass, excessively metal-poor stellar halo; and it has an unusually large nearby satellite galaxy, the Large Magellanic Cloud (LMC). Here we show that the LMC is on a collision course with the MW with which it will merge in $2.4^{+1.2}_{-0.8}$ Gyrs (68% confidence level). This catastrophic and long-overdue event will restore the MW to normality. Using the EAGLE galaxy formation simulation, we show that, as a result of the merger, the central supermassive black hole will increase in mass by up to a factor of 8. The Galactic stellar halo will undergo an equally impressive transformation, becoming 5 times more massive. The additional stars will come predominantly from the disrupted LMC, but a sizeable number will be ejected onto the halo from the stellar disc. The post-merger stellar halo will have the median metallicity of the LMC, [Fe/H]=-0.5 dex, which is typical of other galaxies of similar mass to the MW. At the end of this exceptional event, the MW will become a true benchmark for spiral galaxies, at least temporarily. △ Less

Submitted 17 January, 2019; v1 submitted 24 September, 2018; originally announced September 2018.

Comments: 12 pages, 8 figures, the published version in MNRAS

Journal ref: MNRAS 2019, Volume 483, Issue 2, p.2185-2196

HALO7D I: The Line of Sight Velocities of Distant Main Sequence Stars in the Milky Way Halo

Authors: Emily C. Cunningham, Alis J. Deason, Constance M. Rockosi, Puragra Guhathakurta, Zachary G. Jennings, Evan N. Kirby, Elisa Toloba, Guillermo Barro

Abstract: The Halo Assembly in Lambda-CDM: Observations in 7 Dimensions (HALO7D) dataset consists of Keck II/DEIMOS spectroscopy and Hubble Space Telescope-measured proper motions of Milky Way halo main sequence turnoff stars in the CANDELS fields. In this paper, we present the spectroscopic component of this dataset, and discuss target selection, observing strategy, and survey properties. We present a new… ▽ More The Halo Assembly in Lambda-CDM: Observations in 7 Dimensions (HALO7D) dataset consists of Keck II/DEIMOS spectroscopy and Hubble Space Telescope-measured proper motions of Milky Way halo main sequence turnoff stars in the CANDELS fields. In this paper, we present the spectroscopic component of this dataset, and discuss target selection, observing strategy, and survey properties. We present a new method of measuring line-of-sight (LOS) velocities by combining multiple spectroscopic observations of a given star, utilizing Bayesian hierarchical modeling. We present the LOS velocity distributions of the four HALO7D fields, and estimate their means and dispersions. All of the LOS distributions are dominated by the "hot halo": none of our fields are dominated by substructure that is kinematically cold in the LOS velocity component. Our estimates of the LOS velocity dispersions are consistent across the different fields, and these estimates are consistent with studies using other types of tracers. To complement our observations, we perform mock HALO7D surveys using the synthetic survey software Galaxia to "observe'" the Bullock & Johnston (2005) accreted stellar halos. Based on these simulated datasets, the consistent LOS velocity distributions across the four HALO7D fields indicates that the HALO7D sample is dominated by stars from the same massive (or few relatively massive) accretion event(s). △ Less

Submitted 11 September, 2018; originally announced September 2018.

Comments: 25 pages, 21 figures. Submitted to ApJ; comments welcome!

The mass of the Milky Way from satellite dynamics

Authors: Thomas Callingham, Marius Cautun, Alis J. Deason, Carlos S. Frenk, Wenting Wang, Facundo A. Gómez, Robert J. J. Grand, Federico Marinacci, Rüdgier Pakmor

Abstract: We present and apply a method to infer the mass of the Milky Way (MW) by comparing the dynamics of MW satellites to those of model satellites in the EAGLE cosmological hydrodynamics simulations. A distribution function (DF) for galactic satellites is constructed from EAGLE using specific angular momentum and specific energy, which are scaled so as to be independent of host halo mass. In this 2-dim… ▽ More We present and apply a method to infer the mass of the Milky Way (MW) by comparing the dynamics of MW satellites to those of model satellites in the EAGLE cosmological hydrodynamics simulations. A distribution function (DF) for galactic satellites is constructed from EAGLE using specific angular momentum and specific energy, which are scaled so as to be independent of host halo mass. In this 2-dimensional space, the orbital properties of satellite galaxies vary according to the host halo mass. The halo mass can be inferred by calculating the likelihood that the observed satellite population is drawn from this DF. Our method is robustly calibrated on mock EAGLE systems. We validate it by applying it to the completely independent suite of 30 AURIGA high-resolution simulations of MW-like galaxies: the method accurately recovers their true mass and associated uncertainties. We then apply it to ten classical satellites of the MW with 6D phase-space measurements, including updated proper motions from the GAIA satellite. The mass of the MW is estimated to be $M_{200}^{\textnormal{MW}}=1.17_{-0.15}^{+0.21}\times10^{12}M_{\odot}$ (68\% confidence limits). We combine our total mass estimate with recent mass estimates in the inner regions of the Galaxy to infer an inner dark matter (DM) mass fraction $M^\textnormal{DM}(<20~\rm{kpc})/M^\textnormal{DM}_{200}=0.12$ which is typical of ${\sim}10^{12}M_{\odot}$ $Λ$CDM haloes in hydrodynamical galaxy formation simulations. Assuming an NFW profile, this is equivalent to a halo concentration of $c_{200}^{\textnormal{MW}}=10.9^{+2.6}_{-2.0}$. △ Less

Submitted 23 January, 2019; v1 submitted 30 August, 2018; originally announced August 2018.

Comments: 14 pages, 15 figures, submitted to MNRAS, comments are welcomed

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

Authors: Lachlan Lancaster, Sergey E. Koposov, Vasily Belokurov, N. Wyn Evans, Alis J. Deason

Abstract: 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 sin… ▽ More 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. △ Less

Submitted 14 February, 2020; v1 submitted 11 July, 2018; originally announced July 2018.

Comments: 11 pages, 8 figures

Apocenter Pile-Up: Origin of the Stellar Halo Density Break

Authors: Alis J. Deason, Vasily Belokurov, Sergey E. Koposov, Lachlan Lancaster

Abstract: We measure the orbital properties of halo stars using 7-dimensional information provided by Gaia and the Sloan Digital Sky Survey. A metal-rich population of stars, present in both local main sequence stars and more distant blue horizontal branch stars, have very radial orbits (eccentricity ~ 0.9) and apocenters that coincide with the stellar halo "break radius" at galactocentric distance r ~ 20 k… ▽ More We measure the orbital properties of halo stars using 7-dimensional information provided by Gaia and the Sloan Digital Sky Survey. A metal-rich population of stars, present in both local main sequence stars and more distant blue horizontal branch stars, have very radial orbits (eccentricity ~ 0.9) and apocenters that coincide with the stellar halo "break radius" at galactocentric distance r ~ 20 kpc. Previous work has shown that the stellar halo density falls off much more rapidly beyond this break radius. We argue that the correspondence between the apocenters of high metallicity, high eccentricity stars and the broken density profile is caused by the build-up of stars at the apocenter of a common dwarf progenitor. Although the radially biased stars are likely present down to metallicities of [Fe/H] ~ -2 the increasing dominance at higher metallicities suggests a massive dwarf progenitor, which is at least as massive as the Fornax and Sagittarius dwarf galaxies, and is likely the dominant progenitor of the inner stellar halo. △ Less

Submitted 3 July, 2018; v1 submitted 25 May, 2018; originally announced May 2018.

Comments: 5 pages, 5 figures. Accepted for publication in ApJ Letters

Aurigaia: mock Gaia DR2 stellar catalogues from the Auriga cosmological simulations

Authors: Robert J. J. Grand, John Helly, Azadeh Fattahi, Marius Cautun, Shaun Cole, Andrew P. Cooper, Alis J. Deason, Carlos Frenk, Facundo A. Gómez, Jason A. S. Hunt, Federico Marinacci, Rüdiger Pakmor, Christine M. Simpson, Volker Springel, Dandan Xu

Abstract: We present and analyse mock stellar catalogues that match the selection criteria and observables (including uncertainties) of the Gaia satellite data release 2 (DR2). The source are six cosmological high-resolution magneto-hydrodynamic $Λ$CDM zoom simulations of the formation of Milky Way analogues from the AURIGA project. Mock data are provided for stars with $V < 16$ mag, and $V < 20$ mag at… ▽ More We present and analyse mock stellar catalogues that match the selection criteria and observables (including uncertainties) of the Gaia satellite data release 2 (DR2). The source are six cosmological high-resolution magneto-hydrodynamic $Λ$CDM zoom simulations of the formation of Milky Way analogues from the AURIGA project. Mock data are provided for stars with $V < 16$ mag, and $V < 20$ mag at $|b|>20$ degrees. The mock catalogues are made using two different methods: the public SNAPDRAGONS code, and a method based on that of Lowing et al. that preserves the phase-space distribution of the model stars. These publicly available catalogues contain 5-parameter astrometry, radial velocities, multi-band photometry, stellar parameters, dust extinction values, and uncertainties in all these quantities. In addition, we provide the gravitational potential and information on the origin of each star. By way of demonstration, we apply the mock catalogues to analyses of the young stellar disc and the stellar halo. We show that: i) the young outer stellar disc exhibits a flared distribution that is detectable in the height and vertical velocity distribution of A- and B-dwarf stars up to radii of ~15 kpc; and ii) the spin of the stellar halo out to 100 kpc can be accurately measured with Gaia DR2 RR Lyrae stars. These catalogues are well suited for comparisons with observations and should help to: i) develop and test analysis methods for the Gaia DR2 data; ii) gauge the limitations and biases of the data and iii) interpret the data in the light of theoretical predictions from realistic $ab$ $initio$ simulations of galaxy formation in the $Λ$CDM cosmological model. △ Less

Submitted 5 September, 2018; v1 submitted 23 April, 2018; originally announced April 2018.

Comments: Accepted for publication in MNRAS (1st of September, 2018) after improvements to several figures and some discussion. The mock data are available for access and download at http://icc.dur.ac.uk/data and documentation and supplementary material is available at http://auriga.h-its.org/gaiamock.html

Absolute HST Proper Motion (HSTPROMO) of Distant Milky Way Globular Clusters: Galactocentric Space Velocities and the Milky Way Mass

Authors: Sangmo Tony Sohn, Laura L. Watkins, Mark A. Fardal, Roeland P. van der Marel, Alis J. Deason, Gurtina Besla, Andrea Bellini

Abstract: We present Hubble Space Telescope (HST) absolute proper motion (PM) measurements for 20 globular clusters (GCs) in the Milky Way (MW) halo at Galactocentric distances $R_{\rm GC} \approx 10-100$ kpc, with median per-coordinate PM uncertainty 0.06 mas yr$^{-1}$. Young and old halo GCs do not show systematic differences in their 3D Galactocentric velocities, derived from combination with existing li… ▽ More We present Hubble Space Telescope (HST) absolute proper motion (PM) measurements for 20 globular clusters (GCs) in the Milky Way (MW) halo at Galactocentric distances $R_{\rm GC} \approx 10-100$ kpc, with median per-coordinate PM uncertainty 0.06 mas yr$^{-1}$. Young and old halo GCs do not show systematic differences in their 3D Galactocentric velocities, derived from combination with existing line-of-sight velocities. We confirm the association of Arp 2, Pal 12, Terzan 7, and Terzan 8 with the Sagittarius (Sgr) stream. These clusters and NGC 6101 have tangential velocity $V_{\rm tan} > 290$ km s$^{-1}$, whereas all other clusters have $V_{\rm tan} < 200$ km s$^{-1}$. NGC 2419, the most distant GC in our sample, is also likely associated with the Sgr stream, whereas NGC 4147, NGC 5024, and NGC 5053 definitely are not. We use the distribution of orbital parameters derived using the 3D velocities to separate halo GCs that either formed within the MW or were accreted. We also assess the specific formation history of e.g. Pyxis and Terzan 8. We constrain the MW mass via an estimator that considers the full 6D phase-space information for 16 of the GCs from $R_{\rm GC} = 10$ to 40 kpc. The velocity dispersion anisotropy parameter $β= 0.609^{+0.130}_{-0.229}$. The enclosed mass $M (<39.5 \rm{kpc}) = 0.61^{+0.18}_{-0.12} \times 10^{12}$ M$_{\odot}$, and the virial mass $M_\rm{vir} = 2.05^{+0.97}_{-0.79} \times 10^{12}$ M$_{\odot}$, are consistent with, but on the high side among recent mass estimates in the literature. △ Less

Submitted 14 June, 2018; v1 submitted 5 April, 2018; originally announced April 2018.

Comments: 20 pages, 10 figures, accepted for publication in ApJ

Evolution of LMC/M33-mass dwarf galaxies in the EAGLE simulation

Authors: Shi Shao, Marius Cautun, Alis J. Deason, Carlos S. Frenk, Tom Theuns

Abstract: We investigate the population of dwarf galaxies with stellar masses similar to the Large Magellanic Cloud (LMC) and M33 in the EAGLE galaxy formation simulation. In the field, galaxies reside in haloes with stellar-to-halo mass ratios of $1.03^{+0.50}_{-0.31}\times10^{-2}$ (68% confidence level); systems like the LMC, which have an SMC-mass satellite, reside in haloes about 1.3 times more massive,… ▽ More We investigate the population of dwarf galaxies with stellar masses similar to the Large Magellanic Cloud (LMC) and M33 in the EAGLE galaxy formation simulation. In the field, galaxies reside in haloes with stellar-to-halo mass ratios of $1.03^{+0.50}_{-0.31}\times10^{-2}$ (68% confidence level); systems like the LMC, which have an SMC-mass satellite, reside in haloes about 1.3 times more massive, which suggests an LMC halo mass at infall, $M_{200}=3.4^{+1.8}_{-1.2}\times10^{11}M_\odot$ (68% confidence level). The colour distribution of dwarfs is bimodal, with the red galaxies ($g-r>0.6$) being mostly satellites. The fraction of red LMC-mass dwarfs is 15% for centrals, and for satellites this fraction increases rapidly with host mass: from 10% for satellites of Milky Way (MW)-mass haloes to nearly 90% for satellites of groups and clusters. The quenching timescale, defined as the time after infall when half of the satellites have acquired red colours, decreases with host mass from ${>}5$ Gyrs for MW-mass hosts to $2.5$ Gyrs for cluster mass hosts. The satellites of MW-mass haloes have higher star formation rates and bluer colours than field galaxies. This is due to enhanced star formation triggered by gas compression shortly after accretion. Both the LMC and M33 have enhanced recent star formation that could be a manifestation of this process. After infall into their MW-mass hosts, the $g-r$ colours of LMC-mass dwarfs become bluer for the first 2 Gyrs, after which they rapidly redden. LMC-mass dwarfs fell into their MW-mass hosts only relatively recently, with more than half having an infall time of less than 3.5 Gyrs. △ Less

Submitted 28 June, 2018; v1 submitted 20 March, 2018; originally announced March 2018.

Comments: 12 pages, 14 figures. MNRAS in press