Showing papers by "Justin I. Read published in 2021"

Vintergatan - i. The origins of chemically, kinematically, and structurally distinct discs in a simulated milky way-mass galaxy

Abstract: Spectroscopic surveys of the Milky Way's stars have revealed spatial, chemical and kinematical structures that encode its history. In this work, we study their origins using a cosmological zoom simulation, VINTERGATAN, of a Milky Way-mass disc galaxy. We find that in connection to the last major merger at $z\sim 1.5$, cosmological accretion leads to the rapid formation of an outer, metal-poor, low-[$\alpha$/Fe] gas disc around the inner, metal-rich galaxy containing the old high-[$\alpha$/Fe] stars. This event leads to a bimodality in [$\alpha$/Fe] over a range of [Fe/H]. A detailed analysis of how the galaxy evolves since $z\sim 1$ is presented. We demonstrate the way in which inside-out growth shapes the radial surface density and metallicity profile and how radial migration preferentially relocates stars from the inner to the outer disc. Secular disc heating is found to give rise to increasing velocity dispersions and scaleheights with stellar age, which together with disc flaring explains several trends observed in the Milky Way, including shallower radial [Fe/H]-profiles above the midplane. We show how the galaxy formation scenario imprints non-trivial mappings between structural associations (i.e. thick and thin discs), velocity dispersions, $\alpha$-enhancements, and ages of stars, e.g. the most metal-poor stars in the low-[$\alpha$/Fe] sequence are found to have a scaleheight comparable to old high-[$\alpha$/Fe] stars. Finally, we illustrate how at low spatial resolution, comparable to the thickness of the galaxy, the proposed pathway to distinct sequences in [$\alpha$/Fe]-[Fe/H] cannot be captured.

VINTERGATAN – II. The history of the Milky Way told by its mergers

Abstract: Using the VINTERGATAN cosmological zoom simulation, we explore the contributions of the in situ and accreted material, and the effect of galaxy interactions and mergers in the assembly of a Milky Way-like galaxy. We find that the initial growth phase of galaxy evolution, dominated by repeated major mergers, provides the necessary physical conditions for the assembly of a thick, kinematically hot disk populated by high-[$\alpha$/Fe] stars, formed both in situ and in accreted satellite galaxies. We find that the diversity of evolutionary tracks followed by the simulated galaxy and its progenitors leads to very little overlap of the in situ and accreted populations for any given chemical composition. At a given age, the spread in [$\alpha$/Fe] abundance ratio results from the diversity of physical conditions in VINTERGATAN and its satellites, with an enhancement in [$\alpha$/Fe] found in stars formed during starburst episodes. Later, the cessation of the merger activity promotes the in situ formation of stars in the low-[$\alpha$/Fe] regime, in a radially extended, thin and overall kinematically colder disk, thus establishing chemically bimodal thin and thick disks, in line with observations. We draw links between notable features in the [Fe/H] - [$\alpha$/Fe] plane with their physical causes, and propose a comprehensive formation scenario explaining self-consistently, in the cosmological context, the main observed properties of the Milky Way.

VINTERGATAN III: how to reset the metallicity of the Milky Way

Abstract: Using the cosmological zoom simulation VINTERGATAN, we present a new scenario for the onset of star formation at the metal-poor end of the low-[$\alpha$/Fe] sequence in a Milky Way-like galaxy. In this scenario, the galaxy is fueled by two distinct gas flows. One is enriched by outflows from massive galaxies, but not the other. While the former feeds the inner galactic region, the latter fuels an outer gas disk, inclined with respect to the main galactic plane, and with a significantly poorer chemical content. The first passage of the last major merger galaxy triggers tidal compression in the outer disk, which increases the gas density and eventually leads to star formation, at a metallicity 0.75 dex lower than the inner galaxy. This forms the first stars of the low-[$\alpha$/Fe] sequence. These in situ stars have halo-like kinematics, similarly to what is observed in the Milky Way, due to the inclination of the outer disk which eventually aligns with the inner one via gravitational torques. We show that this tilting disk scenario is likely to be common in Milky-Way like galaxies. This process implies that the low-[$\alpha$/Fe] sequence is populated in situ, simultaneously from two formation channels, in the inner and the outer galaxy, with distinct metallicities. This contrasts with purely sequential scenarios for the assembly of the Milky Way disk and could be tested observationally.

EDGE : Two routes to dark matter core formation in ultra-faint dwarfs

Abstract: In the standard Lambda cold dark matter paradigm, pure dark matter simulations predict dwarf galaxies should inhabit dark matter haloes with a centrally diverging density 'cusp'. This is in conflict with observations that typically favour a constant density 'core'. We investigate this 'cusp-core problem' in 'ultra-faint' dwarf galaxies simulated as part of the 'Engineering Dwarfs at Galaxy formation's Edge' project. We find, similarly to previous work, that gravitational potential fluctuations within the central region of the simulated dwarfs kinematically heat the dark matter particles, lowering the dwarfs' central dark matter density. However, these fluctuations are not exclusively caused by gas inflow/outflow, but also by impulsive heating from minor mergers. We use the genetic modification approach on one of our dwarf's initial conditions to show how a delayed assembly history leads to more late minor mergers and, correspondingly, more dark matter heating. This provides a mechanism by which even ultra-faint dwarfs ($M_∗ \lt 10^5\, \text{M}_{\odot }$), in which star formation was fully quenched at high redshift, can have their central dark matter density lowered over time. In contrast, we find that late major mergers can regenerate a central dark matter cusp, if the merging galaxy had sufficiently little star formation. The combination of these effects leads us to predict significant stochasticity in the central dark matter density slopes of the smallest dwarfs, driven by their unique star formation and mass assembly histories. (Less)

Andromeda XXI - a dwarf galaxy in a low-density dark matter halo

Abstract: Andromeda XXI (And XXI) has been proposed as a dwarf spheroidal galaxy with a central dark matter density that is lower than expected in the Standard $\Lambda$ Cold Dark Matter ($\Lambda$CDM) cosmology. In this work, we present dynamical observations for 77 member stars in this system, more than doubling previous studies to determine whether this galaxy is truly a low density outlier. We measure a systemic velocity of $v_r=-363.4\pm1.0\,{\rm kms}^{-1}$ and a velocity dispersion of $\sigma_v=6.1^{+1.0}_{-0.9}\,{\rm kms}^{-1}$, consistent with previous work and within $1\sigma$ of predictions made within the modified Newtonian dynamics framework. We also measure the metallicity of our member stars from their spectra, finding a mean value of ${\rm [Fe/H]}=-1.7\pm0.1$~dex. We model the dark matter density profile of And~XXI using an improved version of \GravSphere, finding a central density of $\rho_{\rm DM}({\rm 150 pc})=2.7_{-1.7}^{+2.7} \times 10^7 \,{\rm M_\odot\,kpc^{-3}}$ at 68\% confidence, and a density at two half light radii of $\rho_{\rm DM}({\rm 1.75 kpc})=0.9_{-0.2}^{+0.3} \times 10^5 \,{\rm M_\odot\,kpc^{-3}}$ at 68\% confidence. These are both a factor ${\sim}3-5$ lower than the densities expected from abundance matching in $\Lambda$CDM. We show that this cannot be explained by `dark matter heating' since And~XXI had too little star formation to significantly lower its inner dark matter density, while dark matter heating only acts on the profile inside the half light radius. However, And~XXI's low density can be accommodated within $\Lambda$CDM if it experienced extreme tidal stripping (losing $>95\%$ of its mass), or if it inhabits a low concentration halo on a plunging orbit that experienced repeated tidal shocks.

Formation of the largest galactic cores through binary scouring and gravitational wave recoil

Abstract: Massive elliptical galaxies are typically observed to have central cores in their projected radial light profiles. Such cores have long been thought to form through ‘binary scouring’ as supermassive black holes (SMBHs), brought in through mergers, form a hard binary and eject stars from the galactic centre. However, the most massive cores, like the ∼3kpc core in A2261-BCG, remain challenging to explain in this way. In this paper, we run a suite of dry galaxy merger simulations to explore three different scenarios for central core formation in massive elliptical galaxies: ‘binary scouring’, ‘tidal deposition’, and ‘gravitational wave (GW) induced recoil’. Using the GRIFFIN code, we self-consistently model the stars, dark matter, and SMBHs in our merging galaxies, following the SMBH dynamics through to the formation of a hard binary. We find that we can only explain the large surface brightness core of A2261-BCG with a combination of a major merger that produces a small ∼1kpc core through binary scouring, followed by the subsequent GW recoil of its SMBH that acts to grow the core size. Key predictions of this scenario are an offset SMBH surrounded by a compact cluster of bound stars and a non-divergent central density profile. We show that the bright ‘knots’ observed in the core region of A2261-BCG are best explained as stalled perturbers resulting from minor mergers, though the brightest may also represent ejected SMBHs surrounded by a stellar cloak of bound stars.

EDGE: The sensitivity of ultra-faint dwarfs to a metallicity-dependent initial mass function

Abstract: We study how an observationally-motivated, metallicity-dependent initial mass function (IMF) affects the feedback budget and observables of an ultra-faint dwarf galaxy. We model the evolution of a low-mass ($\approx 8 \, \times \, 10^{8} \, \rm M_{\odot}$) dark matter halo with cosmological, zoomed hydrodynamical simulations capable of resolving individual supernovae explosions. We complement the EDGE galaxy formation model from Agertz et al. (2020) with a new prescription for IMF variations according to Geha et al. (2013). At the low metallicities typical of faint dwarf galaxies, the IMF becomes top-heavy, increasing the efficiency of supernova and photo-ionization feedback in regulating star formation. This results in a 100-fold reduction of the final stellar mass of the dwarf compared to a canonical IMF, at fixed dynamical mass. The increase in the feedback budget is nonetheless met by increased metal production from more numerous massive stars, leading to nearly constant iron content at $z=0$. A metallicity-dependent IMF therefore provides a mechanism to produce low-mass ($\rm M_{\star}\sim 10^3 \rm M_{\odot}$), yet enriched ($\rm [Fe/H]\approx -2$) field dwarf galaxies, thus opening a self-consistent avenue to populate the plateau in $\rm [Fe/H]$ at the faintest end of the mass-metallicity relation.

The local dark sector: Probing gravitation’s low-acceleration frontier and dark matter in the Solar System neighborhood

Abstract: We speculate on the development and availability of new innovative propulsion techniques in the 2040s, that will allow us to fly a spacecraft outside the Solar System (at 150 AU and more) in a reasonable amount of time, in order to directly probe our (gravitational) Solar System neighborhood and answer pressing questions regarding the dark sector (dark energy and dark matter). We identify two closely related main science goals, as well as secondary objectives that could be fulfilled by a mission dedicated to probing the local dark sector: (i) begin the exploration of gravitation’s low-acceleration regime with a spacecraft and (ii) improve our knowledge of the local dark matter and baryon densities. Those questions can be answered by directly measuring the gravitational potential with an atomic clock on-board a spacecraft on an outbound Solar System orbit, and by comparing the spacecraft’s trajectory with that predicted by General Relativity through the combination of ranging data and the in-situ measurement (and correction) of non-gravitational accelerations with an on-board accelerometer. Despite a wealth of new experiments getting online in the near future, that will bring new knowledge about the dark sector, it is very unlikely that those science questions will be closed in the next two decades. More importantly, it is likely that it will be even more urgent than currently to answer them. Tracking a spacecraft carrying a clock and an accelerometer as it leaves the Solar System may well be the easiest and fastest way to directly probe our dark environment.

Constraining Ultra Light Dark Matter with the Galactic Nuclear Star Cluster

Abstract: We use the Milky Way's nuclear star cluster (NSC) to test the existence of a dark matter 'soliton core', as predicted in ultra-light dark matter (ULDM) models. Since the soliton core size is proportional to mDM^{-1}, while the core density grows as mDM^2, the NSC (dominant stellar component within about 3 pc) is sensitive to a specific window in the dark matter particle mass, mDM. We apply a spherical isotropic Jeans model to fit the NSC line-of-sight velocity dispersion data, assuming priors on the Milky Way's supermassive black hole (SMBH) mass taken from the Gravity Collaboration et al. (2020) and stellar density profile taken from Gallego-Cano et al. (2018). We find that the current observational data reject the existence of a soliton core for a single ULDM particle with mass in the range 10^{-20.0} < mDM < 10^{-18.5} eV, assuming that the soliton core structure is not affected by the Milky Way's SMBH. We test our methodology on mock data, confirming that we are sensitive to the same range in ULDM mass as for the real data. Dynamical modelling of a larger region of the Galactic centre, including the nuclear stellar disc, promises tighter constraints over a broader range of mDM. We will consider this in future work.