Source author record

Azadeh Fattahi

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

28works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

28 published item(s)

preprint2026arXiv

Local Group analogues in a cosmological context -- I. Relating velocity structure to the cosmic web

Our Local Group, dominated in mass by the Milky Way (MW) and M31, provides a unique laboratory for testing $Λ$CDM cosmology on small scales owing to its proximity. However, its connection to the surrounding large-scale environment, which is essential for interpreting its properties, is inadequately understood. In this work, we explore the connection between Local Group analogues (LGAs) and their surrounding large-scale environments using the ABACUSSUMMIT simulation suite, highlighting the key role of the coupling energy of the MW-M31 orbit, $E_{\rm coupling}$. We find that LGAs with high $E_{\rm coupling}$ preferentially reside in denser regions, whereas those with low $E_{\rm coupling}$ tend to occupy low-density environments. Furthermore, LGAs with low $E_{\rm coupling}$ exhibit strong alignment with cosmic filaments, manifested as a pronounced polar anisotropy in the distribution of tracer haloes. By contrast, LGAs with high $E_{\rm coupling}$ show a weaker polar anisotropy but an enhanced azimuthal anisotropy, with large-scale tracer haloes preferentially lying in the plane spanned by the halo pair and the orbital spin vector. Within this framework, our Local Group is characterised by typical $E_{\rm coupling}$ residing in a relatively under-dense environment, yet it remains consistent with the 95\% range of analogue systems identified in the simulation.

preprint2022arXiv

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

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

preprint2022arXiv

Galactic tides and the Crater II dwarf spheroidal: a challenge to LCDM?

The unusually low velocity dispersion and large size of Crater II pose a challenge to our understanding of dwarf galaxies in the Lambda Cold Dark Matter (LCDM) cosmogony. The low velocity dispersion suggests either a dark halo mass much lower than the minimum expected from hydrogen cooling limit arguments, or one that is in the late stages of extreme tidal stripping. The tidal interpretation has been favoured in recent work and is supported by the small pericentric distances consistent with available kinematic estimates. We use N-body simulations to examine this interpretation in detail, assuming a Navarro-Frenk-White (NFW) profile for Crater II's progenitor halo. Our main finding is that, although the low velocity dispersion can indeed result from the effect of tides, the large size of Crater II is inconsistent with this hypothesis. This is because galaxies stripped to match the observed velocity dispersion are also reduced to sizes much smaller than the observed half-light radius of Crater II. Unless its size has been substantially overestimated, reconciling this system with LCDM requires that either (i) it is not bound and near equilibrium (unlikely, given its crossing time is shorter than the time elapsed since pericentre), or that (ii) its progenitor halo deviates from the assumed NFW profile. The latter alternative may signal that baryons can affect the inner halo cusp even in extremely faint dwarfs or, more intriguingly, may signal effects associated with the intimate nature of the dark matter, such as finite self-interactions, or other such deviations from the canonical LCDM paradigm.

preprint2022arXiv

Satellite mass functions and the faint end of the galaxy mass-halo mass relation in LCDM

The abundance of the faintest galaxies provides insight into the nature of dark matter and the process of dwarf galaxy formation. In the LCDM scenario, low mass halos are so numerous that the efficiency of dwarf formation must decline sharply with decreasing halo mass in order to accommodate the relative scarcity of observed dwarfs and satellites in the Local Group. The nature of this decline contains important clues to the mechanisms regulating the onset of galaxy formation in the faintest systems. We explore here two possible models for the stellar mass ($M_*$)-halo mass ($M_{200}$) relation at the faint end, motivated by some of the latest LCDM cosmological hydrodynamical simulations. One model includes a sharp mass threshold below which no luminous galaxies form, as expected if galaxy formation proceeds only in systems above the Hydrogen-cooling limit. In the second model, $M_*$ scales as a steep power-law of $M_{200}$ with no explicit cutoff, as suggested by recent semianalytic work. Although both models predict satellite numbers around Milky Way-like galaxies consistent with current observations, they predict vastly different numbers of ultra-faint dwarfs and of satellites around isolated dwarf galaxies. Our results illustrate how the satellite mass function around dwarfs may be used to probe the $M_*$-$M_{200}$ relation at the faint end and to elucidate the mechanisms that determine which low-mass halos "light up" or remain dark in the LCDM scenario.

preprint2021arXiv

EAGLE-Auriga: effects of different subgrid models on the baryon cycle around Milky Way-mass galaxies

Modern hydrodynamical simulations reproduce many properties of the real universe. These simulations model various physical processes, but many of these are included using `subgrid models' due to resolution limits. Although different subgrid models have been successful in modelling the effects of supernovae (SNe) feedback on galactic properties, it remains unclear if, and by how much, these differing implementations affect observable halo gas properties. In this work, we use `zoom-in' cosmological initial conditions of two volumes selected to resemble the Local Group (LG) evolved with both the Auriga and EAGLE galaxy formation models. While the subgrid physics models in both simulations reproduce realistic stellar components of $L^\star$ galaxies, they exhibit different gas properties. Namely, Auriga predicts that the Milky Way (MW) is almost baryonically closed, whereas EAGLE suggests that only half of the expected baryons reside within the halo. Furthermore, EAGLE predicts that this baryon deficiency extends to the LG, ($r \leq 1 \mathrm{~Mpc}$). The baryon deficiency in EAGLE is likely due to SNe feedback at high redshift, which generates halo-wide outflows, with high covering fractions and radial velocities, which both eject baryons and significantly impede cosmic gas accretion. Conversely, in Auriga, gas accretion is almost unaffected by feedback. These differences appear to be the result of the different energy injection methods from SNe to gas. Our results suggest that both quasar absorption lines and fast radio burst dispersion measures could constrain these two regimes with future observations.

preprint2021arXiv

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

We use a distribution function analysis to estimate the mass of the Milky Way out to 100 kpc using a large sample of halo stars. These stars are compiled from the literature, and the vast majority (~98%) have 6D phase-space information. We pay particular attention to systematic effects, such as the dynamical influence of the Large Magellanic Cloud (LMC), and the effect of unrelaxed substructure. The LMC biases the (pre-LMC infall) halo mass estimates towards higher values, while realistic stellar halos from cosmological simulations tend to underestimate the true halo mass. After applying our method to the Milky Way data we find a mass within 100 kpc of M(< 100 kpc) = 6.07 +/- 0.29 (stat.) +/- 1.21 (sys.) x 10^11 M_Sun. For this estimate, we have approximately corrected for the reflex motion induced by the LMC using the Erkal et al. model, which assumes a rigid potential for the LMC and MW. Furthermore, stars that likely belong to the Sagittarius stream are removed, and we include a 5% systematic bias, and a 20% systematic uncertainty based on our tests with cosmological simulations. Assuming the mass-concentration relation for Navarro-Frenk-White haloes, our mass estimate favours a total (pre-LMC infall) Milky Way mass of M_200c = 1.01 +/- 0.24 x 10^12 M_Sun, or (post-LMC infall) mass of M_200c = 1.16 +/- 0.24 x 10^12 M_Sun when a 1.5 x 10^11 M_Sun mass of a rigid LMC is included.

preprint2020arXiv

Baryonic clues to the puzzling diversity of dwarf galaxy rotation curves

We use a compilation of disc galaxy rotation curves to assess the role of the luminous component ("baryons") in the rotation curve diversity problem. As in earlier work, we find that rotation curve shape correlates with baryonic surface density: high surface density galaxies have rapidly-rising rotation curves consistent with cuspy cold dark matter halos; slowly-rising rotation curves (characteristic of galaxies with inner mass deficits or "cores") occur only in low surface density galaxies. The correlation, however, seems too weak to be the main driver of the diversity. In addition, dwarf galaxies exhibit a clear trend, from "cuspy" systems where baryons are unimportant in the inner mass budget to "cored" galaxies where baryons actually dominate. This trend constrains the various scenarios proposed to explain the diversity, such as (i) baryonic inflows and outflows during galaxy formation; (ii) dark matter self-interactions; (iii) variations in the baryonic mass structure coupled to rotation velocities through the "mass discrepancy-acceleration relation" (MDAR); or (iv) non-circular motions in gaseous discs. Together with analytical modeling and cosmological hydrodynamical simulations, our analysis shows that each of these scenarios has promising features, but none seems to fully account for the observed diversity. The MDAR, in particular, is inconsistent with the observed trend between rotation curve shape and baryonic importance; either the trend is caused by systematic errors in the data or the MDAR does not apply. The origin of the dwarf galaxy rotation curve diversity and its relation to the structure of cold dark matter halos remains an open issue.

preprint2020arXiv

How Unusual is the Milky Way's Assembly History?

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

preprint2020arXiv

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

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

preprint2020arXiv

The dark matter component of the Gaia radially anisotropic substructure

We study the properties of the dark matter component of the radially anisotropic stellar population recently identified in the Gaia data, using magneto-hydrodynamical simulations of Milky Way-like halos from the Auriga project. We identify 10 simulated galaxies that approximately match the rotation curve and stellar mass of the Milky Way. Four of these have an anisotropic stellar population reminiscent of the Gaia structure. We find an anti-correlation between the dark matter mass fraction of this population in the Solar neighbourhood and its orbital anisotropy. We estimate the local dark matter density and velocity distribution for halos with and without the anisotropic stellar population, and use them to simulate the signals expected in future xenon and germanium direct detection experiments. We find that a generalized Maxwellian distribution fits the dark matter halo integrals of the Milky Way-like halos containing the radially anisotropic stellar population. For dark matter particle masses below approximately 10 GeV, direct detection exclusion limits for the simulated halos with the anisotropic stellar population show a mild shift towards smaller masses compared to the commonly adopted Standard Halo Model.

preprint2020arXiv

The Edge of the Galaxy

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

preprint2020arXiv

The Milky Way total mass profile as inferred from Gaia DR2

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

preprint2019arXiv

Subhalo destruction in the Apostle and Auriga simulations

N-body simulations make unambiguous predictions for the abundance of substructures within dark matter halos. However, the inclusion of baryons in the simulations changes the picture because processes associated with the presence of a large galaxy in the halo can destroy subhalos and substantially alter the mass function and velocity distribution of subhalos. We compare the effect of galaxy formation on subhalo populations in two state-of-the-art sets of hydrodynamical CDM simulations of Milky Way mass halos, APOSTLE and AURIGA. We introduce a new method for tracking the orbits of subhalos between simulation snapshots that gives accurate results down to a few kiloparsecs from the centre of the halo. Relative to a dark matter-only simulation, the abundance of subhalos in APOSTLE is reduced by 50% near the centre and by 10% within r200. In AURIGA the corresponding numbers are 80% and 40%. The velocity distributions of subhalos are also affected by the presence of the galaxy, much more so in AURIGA than in APOSTLE . The differences on subhalo properties in the two simulations can be traced back to the mass of the central galaxies, which in AURIGA are typically twice as massive as those in APOSTLE . We show that some of the results from previous studies are inaccurate due to systematic errors in the modelling of subhalo orbits near the centre of halos.

preprint2019arXiv

The Biggest Splash

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

preprint2019arXiv

The Ophiuchus stream progenitor: a new type of globular cluster and its possible Sagittarius connection

The Ophiuchus stream is a short arc-like stellar feature of uncertain origin located $\sim 5$ kpc North of the Galactic centre. New proper motions from the second $Gaia$ data release reconcile the direction of motion of stream members with the stream arc, resolving a puzzling mismatch reported in earlier work. We use N-body simulations to show that the stream is likely only on its second pericentric passage, and thus was formed recently. The simulations suggest that the entire disrupted progenitor is visible in the observed stream today, and that little further tidal debris lies beyond the ends of the stream. The luminosity, length, width, and velocity dispersion of the stream suggest a globular cluster (GC) progenitor substantially fainter and of lower surface brightness than estimated in previous work, and unlike any other known globulars in the Galaxy. This result suggests the existence of clusters that would extend the known GC population to fainter and more weakly bound systems than hitherto known. How such a weakly-bound cluster of old stars survived until it was disrupted so recently, however, remains a mystery. Integrating backwards in time, we find that the orbits of Sagittarius and Ophiuchus passed within $\sim 5$ kpc of each other about $\sim 100$ Myrs ago, an interaction that might help resolve this puzzle.

preprint2016arXiv

The cold dark matter content of Galactic dwarf spheroidals: no cores, no failures, no problem

We examine the dark matter content of satellite galaxies in Lambda-CDM cosmological hydrodynamical simulations of the Local Group from the APOSTLE project. We find excellent agreement between simulation results and estimates for the 9 brightest Galactic dwarf spheroidals (dSphs) derived from their stellar velocity dispersions and half-light radii. Tidal stripping plays an important role by gradually removing dark matter from the outside in, affecting in particular fainter satellites and systems of larger-than-average size for their luminosity. Our models suggest that tides have significantly reduced the dark matter content of Can Ven I, Sextans, Carina, and Fornax, a prediction that may be tested by comparing them with field galaxies of matching luminosity and size. Uncertainties in observational estimates of the dark matter content of individual dwarfs have been underestimated in the past, at times substantially. We use our improved estimates to revisit the `too-big-to-fail' problem highlighted in earlier N-body work. We reinforce and extend our previous conclusion that the APOSTLE simulations show no sign of this problem. The resolution does not require `cores' in the dark mass profiles, but, rather, relies on revising assumptions and uncertainties in the interpretation of observational data and accounting for `baryon effects' in the theoretical modelling.

preprint2016arXiv

The low abundance and insignificance of dark discs in simulated Milky Way galaxies

We investigate the presence and importance of dark matter discs in a sample of 24 simulated Milky Way galaxies in the APOSTLE project, part of the EAGLE programme of hydrodynamic simulations in Lambda-CDM cosmology. It has been suggested that a dark disc in the Milky Way may boost the dark matter density and modify the velocity modulus relative to a smooth halo at the position of the Sun, with ramifications for direct detection experiments. From a kinematic decomposition of the dark matter and a real space analysis of all 24 halos, we find that only one of the simulated Milky Way analogues has a detectable dark disc component. This unique event was caused by a merger at late time with an LMC-mass satellite at very low grazing angle. Considering that even this rare scenario only enhances the dark matter density at the solar radius by 35% and affects the high energy tail of the dark matter velocity distribution by less than 1%, we conclude that the presence of a dark disc in the Milky Way is unlikely, and is very unlikely to have a significant effect on direct detection experiments.

preprint2016arXiv

The low-mass end of the baryonic Tully-Fisher relation

The scaling of disk galaxy rotation velocity with baryonic mass (the "Baryonic Tully-Fisher" relation, BTF) has long confounded galaxy formation models. It is steeper than the M ~ V^3 scaling relating halo virial masses and circular velocities and its zero point implies that galaxies comprise a very small fraction of available baryons. Such low galaxy formation efficiencies may in principle be explained by winds driven by evolving stars, but the tightness of the BTF relation argues against the substantial scatter expected from such vigorous feedback mechanism. We use the APOSTLE/EAGLE simulations to show that the BTF relation is well reproduced in LCDM simulations that match the size and number of galaxies as a function of stellar mass. In such models, galaxy rotation velocities are proportional to halo virial velocity and the steep velocity-mass dependence results from the decline in galaxy formation efficiency with decreasing halo mass needed to reconcile the CDM halo mass function with the galaxy luminosity function. Despite the strong feedback, the scatter in the simulated BTF is smaller than observed, even when considering all simulated galaxies and not just rotationally-supported ones. The simulations predict that the BTF should become increasingly steep at the faint end, although the velocity scatter at fixed mass should remain small. Observed galaxies with rotation speeds below ~40 km/s seem to deviate from this prediction. We discuss observational biases and modeling uncertainties that may help to explain this disagreement in the context of LCDM models of dwarf galaxy formation.

preprint2016arXiv

The oldest and most metal poor stars in the APOSTLE Local Group simulations

We examine the spatial distribution of the oldest and most metal poor stellar populations of Milky Way-sized galaxies using the APOSTLE cosmological hydrodynamical simulations of the Local Group. In agreement with earlier work, we find strong radial gradients in the fraction of the oldest (tform < 0.8 Gyr) and most metal poor ([Fe/H]< -2.5) stars, both of which increase outwards. The most metal poor stars form over an extended period of time; half of them form after z = 5.3, and the last 10% after z = 2.8. The age of the metal poor stellar population also shows significant variation with environment; a high fraction of them are old in the galaxy's central regions and an even higher fraction in some individual dwarf galaxies, with substantial scatter from dwarf to dwarf. Overall, over half of the stars that belong to both the oldest and most metal-poor population are found outside the solar circle. Somewhat counter-intuitively, we find that dwarf galaxies with a large fraction of metal poor stars that are very old are systems where metal poor stars are relatively rare, but where a substantial old population is present. Our results provide guidance for interpreting the results of surveys designed to hunt for the earliest and most pristine stellar component of our Milky Way.

preprint2016arXiv

Using the Multi-Object Adaptive Optics demonstrator RAVEN to observe metal-poor stars in and towards the Galactic Centre

The chemical abundances for five metal-poor stars in and towards the Galactic bulge have been determined from H-band infrared spectroscopy taken with the RAVEN multi-object adaptive optics science demonstrator and the IRCS spectrograph at the Subaru 8.2-m telescope. Three of these stars are in the Galactic bulge and have metallicities between -2.1 < [Fe/H] < -1.5, and high [alpha/Fe] ~+0.3, typical of Galactic disk and bulge stars in this metallicity range; [Al/Fe] and [N/Fe] are also high, whereas [C/Fe] < +0.3. An examination of their orbits suggests that two of these stars may be confined to the Galactic bulge and one is a halo trespasser, though proper motion values used to calculate orbits are quite uncertain. An additional two stars in the globular cluster M22 show [Fe/H] values consistent to within 1 sigma, although one of these two stars has [Fe/H] = -2.01 +/- 0.09, which is on the low end for this cluster. The [alpha/Fe] and [Ni/Fe] values differ by 2 sigma, with the most metal-poor star showing significantly higher values for these elements. M22 is known to show element abundance variations, consistent with a multi-population scenario (i.e. Marino et al. 2009, 2011; Alves-Brito et al. 2012) though our results cannot discriminate this clearly given our abundance uncertainties. This is the first science demonstration of multi-object adaptive optics with high resolution infrared spectroscopy, and we also discuss the feasibility of this technique for use in the upcoming era of 30-m class telescope facilities.

preprint2015arXiv

Dark matter annihilation radiation in hydrodynamic simulations of Milky Way haloes

We obtain predictions for the properties of cold dark matter annihilation radiation using high resolution hydrodynamic zoom-in cosmological simulations of Milky Way-like galaxies (APOSTLE project) carried out as part of the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) programme. Galactic halos in the simulation have significantly different properties from those assumed in the "standard halo model" often used in dark matter detection studies. The formation of the galaxy causes a contraction of the dark matter halo, whose density profile develops a steeper slope than the Navarro-Frenk-White (NFW) profile between $r\approx1.5$ kpc and $r\approx10$ kpc. At smaller radii, $r\lesssim1.5$ kpc, the halos develop a flatter than NFW slope. This unexpected feature may be specific to our particular choice of subgrid physics model but nevertheless the dark matter density profiles agree within 30% as the mass resolution is increased by a factor 150. The inner regions of the halos are almost perfectly spherical (axis ratios $b/a > 0.97$ within $r=1$ kpc) and there is no offset larger than 45 pc between the centre of the stellar distribution and the centre of the dark halo. The morphology of the predicted dark matter annihilation radiation signal is in broad agreement with $γ$-ray observations at large Galactic latitudes ($b\gtrsim3^\circ$). At smaller angles, the inferred signal in one of our four galaxies is similar to that which is observed but it is significantly weaker in the other three.

preprint2015arXiv

The APOSTLE simulations: solutions to the Local Group's cosmic puzzles

The Local Group of galaxies offer some of the most discriminating tests of models of cosmic structure formation. For example, observations of the Milky Way (MW) and Andromeda satellite populations appear to be in disagreement with N-body simulations of the "Lambda Cold Dark Matter" (ΛCDM) model: there are far fewer satellite galaxies than substructures in cold dark matter halos (the "missing satellites" problem); dwarf galaxies seem to avoid the most massive substructures (the "too-big-to-fail" problem); and the brightest satellites appear to orbit their host galaxies on a thin plane (the "planes of satellites" problem). Here we present results from APOSTLE (A Project Of Simulating The Local Environment), a suite of cosmological hydrodynamic simulations of twelve volumes selected to match the kinematics of the Local Group (LG) members. Applying the Eagle code to the LG environment, we find that our simulations match the observed abundance of LG galaxies, including the satellite galaxies of the MW and Andromeda. Due to changes to the structure of halos and the evolution in the LG environment, the simulations reproduce the observed relation between stellar mass and velocity dispersion of individual dwarf spheroidal galaxies without necessitating the formation of cores in their dark matter profiles. Satellite systems form with a range of spatial anisotropies, including one similar to that of the MW, confirming that such a configuration is not unexpected in ΛCDM. Finally, based on the observed velocity dispersion, size, and stellar mass, we provide new estimates of the maximum circular velocity for the halos of nine MW dwarf spheroidals.

preprint2015arXiv

The unexpected diversity of dwarf galaxy rotation curves

We examine the circular velocity profiles of galaxies in ΛCDM cosmological hydrodynamical simulations from the EAGLE and LOCAL GROUPS projects and compare them with a compilation of observed rotation curves of galaxies spanning a wide range in mass. The shape of the circular velocity profiles of simulated galaxies varies systematically as a function of galaxy mass, but shows remarkably little variation at fixed maximum circular velocity. This is especially true for low-mass dark matter-dominated systems, reflecting the expected similarity of the underlying cold dark matter haloes. This is at odds with observed dwarf galaxies, which show a large diversity of rotation curve shapes, even at fixed maximum rotation speed. Some dwarfs have rotation curves that agree well with simulations, others do not. The latter are systems where the inferred mass enclosed in the inner regions is much lower than expected for cold dark matter haloes and include many galaxies where previous work claims the presence of a constant density "core". The "cusp vs core" issue is thus better characterized as an "inner mass deficit" problem than as a density slope mismatch. For several galaxies the magnitude of this inner mass deficit is well in excess of that reported in recent simulations where cores result from baryon-induced fluctuations in the gravitational potential. We conclude that one or more of the following statements must be true: (i) the dark matter is more complex than envisaged by any current model; (ii) current simulations fail to reproduce the effects of baryons on the inner regions of dwarf galaxies; and/or (iii) the mass profiles of "inner mass deficit" galaxies inferred from kinematic data are incorrect.

preprint2014arXiv

Bent by baryons: the low mass galaxy-halo relation

The relation between galaxies and dark matter halos is of vital importance for evaluating theoretical predictions of structure formation and galaxy formation physics. We show that the widely used method of abundance matching based on dark matter only simulations fails at the low mass end because two of its underlying assumptions are broken: only a small fraction of low mass (below 10^9.5 solar masses) halos host a visible galaxy, and halos grow at a lower rate due to the effect of baryons. In this regime, reliance on dark matter only simulations for abundance matching is neither accurate nor self-consistent. We find that the reported discrepancy between observational estimates of the halo masses of dwarf galaxies and the values predicted by abundance matching does not point to a failure of LCDM, but simply to a failure to account for baryonic effects. Our results also imply that the Local Group contains only a few hundred observable galaxies in contrast with the thousands of faint dwarfs that abundance matching would suggest. We show how relations derived from abundance matching can be corrected, so that they can be used self-consistently to calibrate models of galaxy formation.

preprint2014arXiv

Local Group galaxies emerge from the dark

The "Lambda Cold Dark Matter" (LCDM) model of cosmic structure formation is eminently falsifiable: once its parameters are fixed on large scales, it becomes testable in the nearby Universe. Observations within our Local Group of galaxies, including the satellite populations of the Milky Way and Andromeda, appear to contradict LCDM predictions: there are far fewer satellite galaxies than dark matter halos (the "missing satellites" problem), galaxies seem to avoid the largest substructures (the "too big to fail" problem), and the brightest satellites appear to orbit their host galaxies on a thin plane (the "planes of satellites" problem). We present results from the first hydrodynamic simulations of the Local Group that match the observed abundance of galaxies. We find that when baryonic and dark matter are followed simultaneously in the context of a realistic galaxy formation model, all three "problems" are resolved within the LCDM paradigm.

preprint2014arXiv

The chosen few: the low mass halos that host faint galaxies

Since reionization prevents star formation in most halos below 3 x 10^9 solar masses, dwarf galaxies only populate a fraction of existing dark matter halos. We use hydrodynamic cosmological simulations of the Local Group to study the discriminating factors for galaxy formation in the early Universe and connect them to the present-day properties of galaxies and halos. A combination of selection effects related to reionization, and the subsequent evolution of halos in different environments, introduces strong biases between the population of halos that host dwarf galaxies, and the total halo population. Halos that host galaxies formed earlier and are more concentrated. In addition, halos more affected by tidal stripping are more likely to host a galaxy for a given mass or maximum circular velocity, vmax, today. Consequently, satellite halos are populated more frequently than field halos, and satellite halos of 10^8 - 10^9 solar masses or vmax of 12 - 20 km/s, similar to the Local Group dwarf spheroidals, have experienced a greater than average reduction in both mass and vmax after infall. They are on closer, more radial orbits with higher infall velocities and earlier infall times. Together, these effects make dwarf galaxies highly biased tracers of the underlying dark matter distribution.

preprint2013arXiv

Galaxy Pairs in the Local Group

Current models of galaxy formation predict that galaxy pairs of comparable magnitudes should become increasingly rare with decreasing luminosity. This seems at odds with the relatively high frequency of pairings among dwarf galaxies in the Local Group. We use literature data to show that ~30% of all satellites of the Milky Way and Andromeda galaxies brighter than M_V=-8 are found in likely physical pairs of comparable luminosity. Besides the previously recognised pairings of the Magellanic Clouds and of NGC 147/NGC 185, other candidate pairs include the Ursa Minor and Draco dwarf spheroidals, as well as the And I/And III satellites of M31. These pairs are much closer than expected by chance if the radial and angular distributions of satellites were uncorrelated; in addition, they have very similar line-of-sight velocities and luminosities that differ by less than three magnitudes. In contrast, the same criteria pair fewer than 4% of satellites in N-body/semi-analytic models that match the radial distribution and luminosity function of Local Group satellites. If confirmed in studies of larger samples, the high frequency of dwarf galaxy pairings may provide interesting clues to the formation of faint galaxies in the current cosmological paradigm.

preprint2013arXiv

The Orbital Ellipticity of Satellite Galaxies and the Mass of the Milky Way

We use simulations of Milky Way-sized dark matter haloes from the Aquarius Project to investigate the orbits of substructure haloes likely, according to a semi-analytic galaxy formation model, to host luminous satellites. These tend to populate the most massive subhaloes and are on more radial orbits than the majority of subhaloes found within the halo virial radius. One reason for this (mild) kinematic bias is that many low-mass subhaloes have apocentres that exceed the virial radius of the main host; they are thus excluded from subhalo samples identified within the virial boundary, reducing the number of subhalos on radial orbits. Two other factors contributing to the difference in orbital shape between dark and luminous subhaloes are their dynamical evolution after infall, which affects more markedly low-mass (dark) subhaloes, and a weak dependence of ellipticity on the redshift of first infall. The ellipticity distribution of luminous satellites exhibits little halo-to-halo scatter and it may therefore be compared fruitfully with that of Milky Way satellites. Since the latter depends sensitively on the total mass of the Milky Way we can use the predicted distribution of satellite ellipticities to place constraints on this important parameter. Using the latest estimates of position and velocity of dwarfs compiled from the literature, we find that the most likely Milky Way mass lies in the range $6 \times 10^{11} M_{\odot} < M_{200} < 3.1 \times 10^{12} M_{\odot}$, with a best fit value of $M_{200} = 1.1 \times 10^{12} M_{\odot}$. This value is consistent with Milky Way mass estimates based on dynamical tracers or the timing argument.