Showing 1–44 of 44 results for author: Sawala, T

Distinct distributions of elliptical and disk galaxies across the Local Supercluster as a $Λ$CDM prediction

Authors: Till Sawala, Carlos Frenk, Jens Jasche, Peter H. Johansson, Guilhem Lavaux

Abstract: Galaxies of different types are not equally distributed in the Local Universe. In particular, the supergalactic plane is prominent among the brightest ellipticals, but inconspicuous among the brightest disk galaxies. This striking difference provides a unique test for our understanding of galaxy and structure formation. Here we use the SIBELIUS DARK constrained simulation to confront the predictio… ▽ More Galaxies of different types are not equally distributed in the Local Universe. In particular, the supergalactic plane is prominent among the brightest ellipticals, but inconspicuous among the brightest disk galaxies. This striking difference provides a unique test for our understanding of galaxy and structure formation. Here we use the SIBELIUS DARK constrained simulation to confront the predictions of the standard Lambda Cold Dark Matter ($Λ$CDM) model and standard galaxy formation theory with these observations. We find that SIBELIUS DARK reproduces the spatial distributions of disks and ellipticals and, in particular, the observed excess of massive ellipticals near the supergalactic equator. We show that this follows directly from the local large-scale structure and from the standard galaxy formation paradigm, wherein disk galaxies evolve mostly in isolation, while giant ellipticals congregate in the massive clusters that define the supergalactic plane. Rather than being anomalous as earlier works have suggested, the distributions of giant ellipticals and disks in the Local Universe and in relation to the supergalactic plane are key predictions of the $Λ$CDM model. △ Less

Submitted 5 December, 2023; v1 submitted 1 November, 2023; originally announced November 2023.

Comments: 21 pages, 9 figures, published in Nature Astronomy

The Timeless Timing Argument and the Mass of the Local Group

Authors: Till Sawala, Jorge Peñarrubia, Shihong Liao, Peter H. Johansson

Abstract: The Timing Argument connects the motion of a two-body system to its mass in an expanding Universe with a finite age, under the assumption that it has evolved on a self-gravitating orbit. It is commonly applied to the present-day Milky Way-M31 system in order to infer its unknown mass from the measured kinematics. We use a set of Local Group analogues from the Uchuu simulation to investigate the Ti… ▽ More The Timing Argument connects the motion of a two-body system to its mass in an expanding Universe with a finite age, under the assumption that it has evolved on a self-gravitating orbit. It is commonly applied to the present-day Milky Way-M31 system in order to infer its unknown mass from the measured kinematics. We use a set of Local Group analogues from the Uchuu simulation to investigate the Timing Argument over cosmic time. We find that the median inferred mass remains almost constant over the past 12 Gyr, even while the haloes themselves grew in mass by more than an order of magnitude. By contrast, we find a closer, and nearly time-invariant agreement between the Timing Argument value and the mass within a sphere of radius equal to the MW-M31 separation, and we identify this as the total mass of the system. We conclude that the comparatively close present-day agreement between the Timing Argument and the sum of the halo masses reflects no underlying relation, but merely echoes the fact that the MW and M31 now contain most (but not all) of the mass of the Local Group system. △ Less

Submitted 18 August, 2023; v1 submitted 25 July, 2023; originally announced July 2023.

Comments: 6 pages, 4 figures, this version accepted to MNRAS Letters

Modelling the accretion and feedback of supermassive black hole binaries in gas-rich galaxy mergers

Authors: Shihong Liao, Peter H. Johansson, Matias Mannerkoski, Dimitrios Irodotou, Francesco Paolo Rizzuto, Stuart McAlpine, Antti Rantala, Alexander Rawlings, Till Sawala

Abstract: We introduce a new model for the accretion and feedback of supermassive black hole (SMBH) binaries to the KETJU code, which enables us to resolve the evolution of SMBH binaries down to separations of tens of Schwarzschild radii in gas-rich galaxy mergers. Our subgrid binary accretion model extends the widely used Bondi--Hoyle--Lyttleton accretion into the binary phase and incorporates preferential… ▽ More We introduce a new model for the accretion and feedback of supermassive black hole (SMBH) binaries to the KETJU code, which enables us to resolve the evolution of SMBH binaries down to separations of tens of Schwarzschild radii in gas-rich galaxy mergers. Our subgrid binary accretion model extends the widely used Bondi--Hoyle--Lyttleton accretion into the binary phase and incorporates preferential mass accretion onto the secondary SMBH, which is motivated by results from small-scale hydrodynamical circumbinary disc simulations. We perform idealised gas-rich disc galaxy merger simulations using pure thermal or pure kinetic active galactic nuclei (AGN) feedback. Our binary accretion model provides more physically motivated SMBH mass ratios, which are one of the key parameters for computing gravitational wave (GW) induced recoil velocities. The merger time-scales of our simulated SMBH binaries are in the range $t_{\rm merge}{\sim} 10$--$400$ Myr. Prograde in-plane equal-mass galaxy mergers lead to the shortest merger time-scales, as they experience the strongest starbursts, with the ensuing high stellar density resulting in a rapid SMBH coalescence. Compared to the thermal AGN feedback, the kinetic AGN feedback predicts longer merger time-scales and results in more core-like stellar profiles, as it is more effective in removing gas from the galaxy centre and quenching star formation. This suggests that the AGN feedback implementation plays a critical role in modelling SMBH coalescences. Our model will be useful for improving the modelling of SMBH mergers in gas-rich galaxies, the prime targets for the upcoming LISA GW observatory. △ Less

Submitted 18 February, 2023; v1 submitted 21 November, 2022; originally announced November 2022.

Comments: 28 pages, 13 main + 6 appendix figures, accepted for publication in MNRAS

Journal ref: MNRAS 520, 4463 (2023)

The Local Group's mass: probably no more than the sum of its parts

Authors: Till Sawala, Meri Teeriaho, Peter H. Johansson

Abstract: The total mass of the Local Group (LG) and the masses of its primary constituents, the Milky Way and M31, are important anchors for several cosmological questions. In recent years, independent measurements have consistently yielded halo masses close to $10^{12} \mathrm{M_\odot}$ for the MW, and $1-2 \times 10^{12} \mathrm{M_\odot}$ for M31, while estimates derived from the pair's kinematics via th… ▽ More The total mass of the Local Group (LG) and the masses of its primary constituents, the Milky Way and M31, are important anchors for several cosmological questions. In recent years, independent measurements have consistently yielded halo masses close to $10^{12} \mathrm{M_\odot}$ for the MW, and $1-2 \times 10^{12} \mathrm{M_\odot}$ for M31, while estimates derived from the pair's kinematics via the `timing argument' have yielded a combined mass of around $5 \times 10^{12} \mathrm{M_\odot}$. Here, we analyse the extremely large Uchuu simulation to constrain the mass of the Local Group and its two most massive members. First, we demonstrate the importance of selecting LG analogues whose kinematics are dominated by mutual interactions to a similar extent as the LG. Adopting the observed separation and radial velocity, we obtain a weighted posterior of $75_{-40}^{+65}$ kms$^{-1}$ for the uncertain transverse velocity. Via Gaussian process regression, we infer a total mass of $3.2^{+1.2}_{-0.9} \times 10^{12} \mathrm{M_\odot}$, significantly below the timing argument prediction. Importantly, we show that the remaining uncertainty is not rooted in the analysis or observational errors, but in the irreducible scatter in the kinematics-mass relation. We further find a mass for the less massive halo of $0.9_{-0.3}^{+0.6} \times 10^{12} \mathrm{M_\odot}$ and for the more massive halo of $2.3_{-0.9}^{+1.0} \times 10^{12} \mathrm{M_\odot}$, consistent with independent measurements of the masses of MW and M31, respectively. Incorporating the mass of the MW as an additional prior allows us to further constrain all measurements and determine that the MW is very likely to be the lower mass object of the two. △ Less

Submitted 13 October, 2022; originally announced October 2022.

Comments: 15 pages, submitted to MNRAS. Full source code provided, comments are welcome

Velocity-dependent J-factors for Milky Way dwarf spheroidal analogues in cosmological simulations

Authors: Keagan Blanchette, Erin Piccirillo, Nassim Bozorgnia, Louis E. Strigari, Azadeh Fattahi, Carlos S. Frenk, Julio F. Navarro, Till Sawala

Abstract: We study the impact of the dark matter velocity distribution modelling on signals from velocity-dependent dark matter annihilation in Milky Way dwarf spheroidal galaxies. Using the high resolution APOSTLE simulations, we identify analogues corresponding to Milky Way dwarf spheroidal galaxies, and from these directly determine the dark matter pair-wise relative velocity distribution, and compare to… ▽ More We study the impact of the dark matter velocity distribution modelling on signals from velocity-dependent dark matter annihilation in Milky Way dwarf spheroidal galaxies. Using the high resolution APOSTLE simulations, we identify analogues corresponding to Milky Way dwarf spheroidal galaxies, and from these directly determine the dark matter pair-wise relative velocity distribution, and compare to best-fitting Maxwell-Boltzmann distribution models. For three velocity-dependent annihilation models, p-wave, d-wave, and the Sommerfeld model, we quantify the errors introduced when using the Maxwell-Boltzmann parameterization. We extract a simple power-law relation between the maximum circular velocity of the dwarf spheroidal analogue and the peak speed of the Maxwell-Boltzmann distribution. We show that this relation can be used to accurately calculate the dark matter relative velocity distribution, and find that it allows us to estimate the dark matter annihilation signal without the need to directly calculate the relative velocity distribution for each galaxy. The scatter in the J-factors calculated from the analogues dominates the uncertainty obtained when compared to the J-factor as determined from the observational data for each dwarf spheroidal, with the largest scatter from d-wave models and the smallest from Sommerfeld models. △ Less

Submitted 30 June, 2022; originally announced July 2022.

Comments: 17 pages, 7 figures, 2 tables

The Milky Way's plane of satellites: consistent with $Λ$CDM

Authors: Till Sawala, Marius Cautun, Carlos S. Frenk, John Helly, Jens Jasche, Adrian Jenkins, Peter H. Johansson, Guilhem Lavaux, Stuart McAlpine, Matthieu Schaller

Abstract: The "plane of satellites problem" describes the arrangement of the Milky Way's 11 brightest satellite galaxies in a remarkably thin plane, possibly supported by rotation. This is in apparent contradiction to the standard cosmological model, wherein the Galaxy is surrounded by a dispersion-supported dark matter halo. Here, we show that the reported exceptional anisotropy of the satellite system is… ▽ More The "plane of satellites problem" describes the arrangement of the Milky Way's 11 brightest satellite galaxies in a remarkably thin plane, possibly supported by rotation. This is in apparent contradiction to the standard cosmological model, wherein the Galaxy is surrounded by a dispersion-supported dark matter halo. Here, we show that the reported exceptional anisotropy of the satellite system is strongly contingent on a lopsided radial distribution, which earlier simulations have failed to reproduce, combined with the close but fleeting conjunction of the two most distant satellites, Leo I and Leo II. Using Gaia proper motions, we show that the orbital pole alignment is much more common than previously reported, and reveal the plane of satellites to be transient rather than rotationally supported. Comparing to new simulations, where such short-lived planes are common, we find the Milky Way satellites to be compatible with standard model expectations. △ Less

Submitted 5 May, 2022; originally announced May 2022.

SIBELIUS-DARK: a galaxy catalogue of the Local Volume from a constrained realisation simulation

Authors: Stuart McAlpine, John C. Helly, Matthieu Schaller, Till Sawala, Guilhem Lavaux, Jens Jasche, Carlos S. Frenk, Adrian Jenkins, John R. Lucey, Peter H. Johansson

Abstract: We present SIBELIUS-DARK, a constrained realisation simulation of the local volume to a distance of 200~Mpc from the Milky Way. SIBELIUS-DARK is the first study of the \textit{Simulations Beyond The Local Universe} (SIBELIUS) project, which has the goal of embedding a model Local Group-like system within the correct cosmic environment. The simulation is dark-matter-only, with the galaxy population… ▽ More We present SIBELIUS-DARK, a constrained realisation simulation of the local volume to a distance of 200~Mpc from the Milky Way. SIBELIUS-DARK is the first study of the \textit{Simulations Beyond The Local Universe} (SIBELIUS) project, which has the goal of embedding a model Local Group-like system within the correct cosmic environment. The simulation is dark-matter-only, with the galaxy population calculated using the semi-analytic model of galaxy formation, GALFORM. We demonstrate that the large-scale structure that emerges from the SIBELIUS constrained initial conditions matches well the observational data. The inferred galaxy population of SIBELIUS-DARK also match well the observational data, both statistically for the whole volume and on an object-by-object basis for the most massive clusters. For example, the $K$-band number counts across the whole sky, and when divided between the northern and southern Galactic hemispheres, are well reproduced by SIBELIUS-DARK. We find that the local volume is somewhat unusual in the wider context of $Λ$CDM: it contains an abnormally high number of supermassive clusters, as well as an overall large-scale underdensity at the level of $\approx 5$\% relative to the cosmic mean. However, whilst rare, the extent of these peculiarities does not significantly challenge the $Λ$CDM model. SIBELIUS-DARK is the most comprehensive constrained realisation simulation of the local volume to date, and with this paper we publicly release the halo and galaxy catalogues at $z=0$, which we hope will be useful to the wider astronomy community. △ Less

Submitted 8 February, 2022; originally announced February 2022.

Comments: 26 pages, 19 figures. Paper as been accepted for publication in MNRAS. Comments welcome

The SIBELIUS Project: E Pluribus Unum

Authors: Till Sawala, Stuart McAlpine, Jens Jasche, Guilhem Lavaux, Adrian Jenkins, Peter H. Johansson, Carlos S. Frenk

Abstract: We introduce "Simulations Beyond The Local Universe" (SIBELIUS) that connect the Local Group to its cosmic environment. We show that introducing hierarchical small-scale perturbations to a density field constrained on large scales by observations provides an efficient way to explore the sample space of Local Group analogues. From more than 60 000 simulations, we identify a hierarchy of Local Group… ▽ More We introduce "Simulations Beyond The Local Universe" (SIBELIUS) that connect the Local Group to its cosmic environment. We show that introducing hierarchical small-scale perturbations to a density field constrained on large scales by observations provides an efficient way to explore the sample space of Local Group analogues. From more than 60 000 simulations, we identify a hierarchy of Local Group characteristics emanating from different scales: the total mass, orientation, orbital energy and the angular momentum are largely determined by modes above $λ$ = 1.6 comoving Mpc (cMpc) in the primordial density field. Smaller scale variations are mostly manifest as perturbations to the MW-M31 orbit, and we find that the observables commonly used to describe the Local Group -- the MW-M31 separation and radial velocity -- are transient and depend on specifying scales down to 0.2 cMpc in the primordial density field. We further find that the presence of M33/LMC analogues significantly affects the MW-M31 orbit and its sensitivity to small-scale perturbations. We construct initial conditions that lead to the formation of a Local Group whose primary observables precisely match the current observations. △ Less

Submitted 22 March, 2021; originally announced March 2021.

Comments: 14 pages, 15 figures. Submitted to MNRAS, comments welcome!

Setting the Stage: Structures from Gaussian Random Fields

Authors: Till Sawala, Adrian Jenkins, Stuart McAlpine, Jens Jasche, Guilhem Lavaux, Peter H. Johansson, Carlos S. Frenk

Abstract: We study structure formation in a set of cosmological simulations to uncover the scales in the initial density field that gave rise to the formation of present-day structures. Our simulations share a common primordial power spectrum (here Lambda-CDM), but the introduction of hierarchical variations of the phase information allows us to systematically study the scales that determine the formation o… ▽ More We study structure formation in a set of cosmological simulations to uncover the scales in the initial density field that gave rise to the formation of present-day structures. Our simulations share a common primordial power spectrum (here Lambda-CDM), but the introduction of hierarchical variations of the phase information allows us to systematically study the scales that determine the formation of structure at later times. We consider the variance in z=0 statistics such as the matter power spectrum and halo mass function. We also define a criterion for the existence of individual haloes across simulations, and determine what scales in the initial density field contain sufficient information for the non-linear formation of unique haloes. We study how the characteristics of individual haloes such as the mass and concentration, as well as the position and velocity, are affected by variations on different scales, and give scaling relations for haloes of different mass. Finally, we use the example of a cluster-mass halo to show how our hierarchical parametrisation of the initial density field can be used to create variants of particular objects. With properties such as mass, concentration, kinematics and substructure of haloes set on distinct and well-determined scales, and its unique ability to introduce variations localised in real space, our method is a powerful tool to study structure formation in cosmological simulations. △ Less

Submitted 9 March, 2020; originally announced March 2020.

Comments: 18 pages, 22 figures, submitted to MNRAS

Structure Formation with Two Periods of Inflation: Beyond PLaIn $Λ$CDM

Authors: Kari Enqvist, Till Sawala, Tomo Takahashi

Abstract: We discuss structure formation in models with a spectator field in small-field inflation which accommodate a secondary period of inflation. In such models, subgalactic scale primordial fluctuations can be much suppressed in comparison to the usual power-law $Λ$CDM model while the large scale fluctuations remain consistent with current observations. We discuss how a secondary inflationary epoch may… ▽ More We discuss structure formation in models with a spectator field in small-field inflation which accommodate a secondary period of inflation. In such models, subgalactic scale primordial fluctuations can be much suppressed in comparison to the usual power-law $Λ$CDM model while the large scale fluctuations remain consistent with current observations. We discuss how a secondary inflationary epoch may give rise to observable features in the small scale power spectrum and hence be tested by the structures in the Local Universe. △ Less

Submitted 31 May, 2019; originally announced May 2019.

Comments: 20 pages, 4 figures

The core-cusp problem: a matter of perspective

Authors: Anna Genina, Alejandro Benitez-Llambay, Carlos S. Frenk, Shaun Cole, Azadeh Fattahi, Julio F. Navarro, Kyle A. Oman, Till Sawala, Tom Theuns

Abstract: The existence of two kinematically and chemically distinct stellar subpopulations in the Sculptor and Fornax dwarf galaxies offers the opportunity to constrain the density profile of their matter haloes by measuring the mass contained within the well-separated half-light radii of the two metallicity subpopulations. Walker and Penarrubia have used this approach to argue that data for these galaxies… ▽ More The existence of two kinematically and chemically distinct stellar subpopulations in the Sculptor and Fornax dwarf galaxies offers the opportunity to constrain the density profile of their matter haloes by measuring the mass contained within the well-separated half-light radii of the two metallicity subpopulations. Walker and Penarrubia have used this approach to argue that data for these galaxies are consistent with constant-density `cores' in their inner regions and rule out `cuspy' Navarro-Frenk-White (NFW) profiles with high statistical significance, particularly in the case of Sculptor. We test the validity of these claims using dwarf galaxies in the APOSTLE (A Project Of Simulating The Local Environment) Lambda cold dark matter cosmological hydrodynamics simulations of analogues of the Local Group. These galaxies all have NFW dark matter density profiles and a subset of them develop two distinct metallicity subpopulations reminiscent of Sculptor and Fornax. We apply a method analogous to that of Walker and Penarrubia to a sample of 53 simulated dwarfs and find that this procedure often leads to a statistically significant detection of a core in the profile when in reality there is a cusp. Although multiple factors contribute to these failures, the main cause is a violation of the assumption of spherical symmetry upon which the mass estimators are based. The stellar populations of the simulated dwarfs tend to be significantly elongated and, in several cases, the two metallicity populations have different asphericity and are misaligned. As a result, a wide range of slopes of the density profile are inferred depending on the angle from which the galaxy is viewed. △ Less

Submitted 13 December, 2017; v1 submitted 19 July, 2017; originally announced July 2017.

Comments: 15 pages, 13 figures, accepted for publication in MNRAS

Tidal stripping and the structure of dwarf galaxies in the Local Group

Authors: Azadeh Fattahi, Julio F. Navarro, Carlos S. Frenk, Kyle Oman, Till Sawala, Matthieu Schaller

Abstract: The shallow faint-end slope of the galaxy mass function is usually reproduced in $Λ$CDM galaxy formation models by assuming that the fraction of baryons that turns into stars drops steeply with decreasing halo mass and essentially vanishes in haloes with maximum circular velocities $V_{\rm max}<20$-$30$ km/s. Dark matter-dominated dwarfs should therefore have characteristic velocities of about tha… ▽ More The shallow faint-end slope of the galaxy mass function is usually reproduced in $Λ$CDM galaxy formation models by assuming that the fraction of baryons that turns into stars drops steeply with decreasing halo mass and essentially vanishes in haloes with maximum circular velocities $V_{\rm max}<20$-$30$ km/s. Dark matter-dominated dwarfs should therefore have characteristic velocities of about that value, unless they are small enough to probe only the rising part of the halo circular velocity curve (i.e., half-mass radii, $r_{1/2}\ll 1$ kpc). Many dwarfs have properties in disagreement with this prediction: they are large enough to probe their halo $V_{\rm max}$ but their characteristic velocities are well below $20$ km/s. These `cold faint giants' (an extreme example is the recently discovered Crater 2 Milky Way satellite) can only be reconciled with our $Λ$CDM models if they are the remnants of once massive objects heavily affected by tidal stripping. We examine this possibility using the APOSTLE cosmological hydrodynamical simulations of the Local Group. Assuming that low velocity dispersion satellites have been affected by stripping, we infer their progenitor masses, radii, and velocity dispersions, and find them in remarkable agreement with those of isolated dwarfs. Tidal stripping also explains the large scatter in the mass discrepancy-acceleration relation in the dwarf galaxy regime: tides remove preferentially dark matter from satellite galaxies, lowering their accelerations below the $a_{\rm min}\sim 10^{-11} m/s^2$ minimum expected for isolated dwarfs. In many cases, the resulting velocity dispersions are inconsistent with the predictions from Modified Newtonian Dynamics, a result that poses a possibly insurmountable challenge to that scenario. △ Less

Submitted 15 March, 2018; v1 submitted 12 July, 2017; originally announced July 2017.

Comments: replaced to match the accepted version. 2 new figures added compared to the previous version

Post-Newtonian dynamical modeling of supermassive black holes in galactic-scale simulations

Authors: Antti Rantala, Pauli Pihajoki, Peter H. Johansson, Thorsten Naab, Natalia Lahén, Till Sawala

Abstract: We present KETJU, a new extension of the widely-used smoothed particle hydrodynamics simulation code GADGET-3. The key feature of the code is the inclusion of algorithmically regularized regions around every supermassive black hole (SMBH). This allows for simultaneously following global galactic-scale dynamical and astrophysical processes, while solving the dynamics of SMBHs, SMBH binaries and sur… ▽ More We present KETJU, a new extension of the widely-used smoothed particle hydrodynamics simulation code GADGET-3. The key feature of the code is the inclusion of algorithmically regularized regions around every supermassive black hole (SMBH). This allows for simultaneously following global galactic-scale dynamical and astrophysical processes, while solving the dynamics of SMBHs, SMBH binaries and surrounding stellar systems at sub-parsec scales. The KETJU code includes Post-Newtonian terms in the equations of motions of the SMBHs which enables a new SMBH merger criterion based on the gravitational wave coalescence timescale pushing the merger separation of SMBHs down to $\sim 0.005$ pc. We test the performance of our code by comparison to NBODY7 and rVINE. We set up dynamically stable multi-component merger progenitor galaxies to study the SMBH binary evolution during galaxy mergers. In our simulation sample the SMBH binaries do not suffer from the final-parsec problem, which we attribute to the non-spherical shape of the merger remnants. For bulge-only models, the hardening rate decreases with increasing resolution, whereas for models which in addition include massive dark matter halos the SMBH binary hardening rate becomes practically independent of the mass resolution of the stellar bulge. The SMBHs coalesce on average 200 Myr after the formation of the SMBH binary. However, small differences in the initial SMBH binary eccentricities can result in large differences in the SMBH coalescence times. Finally, we discuss the future prospects of KETJU, which allows for a straightforward inclusion of gas physics in the simulations. △ Less

Submitted 19 April, 2017; v1 submitted 21 November, 2016; originally announced November 2016.

Comments: 28 pages, 18 figures. Accepted for publication in ApJ

Tidal features of classical Milky Way satellites in a $Λ$CDM universe

Authors: M. -Y. Wang, Azadeh Fattahi, Andrew P. Cooper, Till Sawala, Louis E. Strigari, Carlos S. Frenk, Julio F. Navarro, Kyle Oman, Matthieu Schaller

Abstract: We use the APOSTLE cosmological hydrodynamic simulations to examine the effects of tidal stripping on cold dark matter (CDM) sub haloes that host three of the most luminous Milky Way (MW) dwarf satellite galaxies: Fornax, Sculptor, and Leo I. We identify simulated satellites that match the observed spatial and kinematic distributions of stars in these galaxies, and track their evolution after infa… ▽ More We use the APOSTLE cosmological hydrodynamic simulations to examine the effects of tidal stripping on cold dark matter (CDM) sub haloes that host three of the most luminous Milky Way (MW) dwarf satellite galaxies: Fornax, Sculptor, and Leo I. We identify simulated satellites that match the observed spatial and kinematic distributions of stars in these galaxies, and track their evolution after infall. We find $\sim$ 30$\%$ of subhaloes hosting satellites with present-day stellar mass $10^6$-$10^8$ $M_{\odot}$ experience $>20\%$ stellar mass loss after infall. Fornax analogues have earlier infall times compared to Sculptor and Leo I analogues. Star formation in Fornax analogues continues for $\sim3$-$6$ Gyr after infall, whereas Sculptor and Leo I analogues stop forming stars $< 2$-$3$ Gyr after infall. Fornax analogues typically show more significant stellar mass loss and exhibit stellar tidal tails, whereas Sculptor and Leo I analogues, which are more deeply embedded in their host DM haloes at infall, do not show substantial mass loss due to tides. When additionally comparing the orbital motion of the host subaloes to the measured proper motion of Fornax we find the matching more difficult; host subhaloes tend to have pericentres smaller than that measured for Fornax itself. From the kinematic and orbital data, we estimate that Fornax has lost $10-20\%$ of its infall stellar mass. Our best estimate for the surface brightness of a stellar tidal stream associated with Fornax is $Σ\sim$ 32.6 mag $ {\rm arcsec^{-2}}$, which may be detectable with deep imaging surveys such as DES and LSST. △ Less

Submitted 17 November, 2016; v1 submitted 2 November, 2016; originally announced November 2016.

Comments: 15 pages, 12 figures, submitted to MNRAS

Properties of Local Group galaxies in hydrodynamical simulations of sterile neutrino dark matter cosmologies

Authors: Mark R. Lovell, Sownak Bose, Alexey Boyarsky, Robert A. Crain, Carlos S. Frenk, Wojciech A. Hellwing, Aaron D. Ludlow, Julio F. Navarro, Oleg Ruchayskiy, Till Sawala, Matthieu Schaller, Joop Schaye, Tom Theuns

Abstract: We study galaxy formation in sterile neutrino dark matter models that differ signifi- cantly from both cold and from `warm thermal relic' models. We use the EAGLE code to carry out hydrodynamic simulations of the evolution of pairs of galaxies chosen to resemble the Local Group, as part of the APOSTLE simulations project. We compare cold dark matter (CDM) with two sterile neutrino models with 7 ke… ▽ More We study galaxy formation in sterile neutrino dark matter models that differ signifi- cantly from both cold and from `warm thermal relic' models. We use the EAGLE code to carry out hydrodynamic simulations of the evolution of pairs of galaxies chosen to resemble the Local Group, as part of the APOSTLE simulations project. We compare cold dark matter (CDM) with two sterile neutrino models with 7 keV mass: one, the warmest among all models of this mass (LA120) and the other, a relatively cold case (LA10). We show that the lower concentration of sterile neutrino subhalos compared to their CDM counterparts makes the inferred inner dark matter content of galaxies like Fornax (or Magellanic Clouds) less of an outlier in the sterile neutrino cosmologies. In terms of the galaxy number counts the LA10 simulations are indistinguishable from CDM when one takes into account halo-to-halo (or `simulation-to-simulation') scatter. In order for the LA120 model to match the number of Local Group dwarf galaxies, a higher fraction of low mass haloes is required to form galaxies than is predicted by the EAGLE simulations. As the census of the Local Group galaxies nears completion, this population may provide a strong discriminant between cold and warm dark matter models. △ Less

Submitted 22 August, 2017; v1 submitted 31 October, 2016; originally announced November 2016.

Comments: 15 pages, 15 figures, revised version accepted to MNRAS with minor changes, contact: lovell@mpia.de

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

Authors: Else Starkenburg, Kyle A. Oman, Julio F. Navarro, Robert A. Crain, Azadeh Fattahi, Carlos S. Frenk, Till Sawala, Joop Schaye

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

Submitted 16 September, 2016; originally announced September 2016.

Comments: 14 pages, 9 figures. Submitted to MNRAS, revised after referee comments

Shaken and Stirred: The Milky Way's Dark Substructures

Authors: Till Sawala, Pauli Pihajoki, Peter H. Johansson, Carlos S. Frenk, Julio F. Navarro, Kyle A. Oman, Simon D. M. White

Abstract: The predicted abundance and properties of the low-mass substructures embedded inside larger dark matter haloes differ sharply among alternative dark matter models. Too small to host galaxies themselves, these subhaloes may still be detected via gravitational lensing, or via perturbations of the Milky Way's globular cluster streams and its stellar disk. Here we use the Apostle cosmological simulati… ▽ More The predicted abundance and properties of the low-mass substructures embedded inside larger dark matter haloes differ sharply among alternative dark matter models. Too small to host galaxies themselves, these subhaloes may still be detected via gravitational lensing, or via perturbations of the Milky Way's globular cluster streams and its stellar disk. Here we use the Apostle cosmological simulations to predict the abundance and the spatial and velocity distributions of subhaloes in the range 10^6.5-10^8.5 solar masses inside haloes of mass ~ 10^12 solar masses in LCDM. Although these subhaloes are themselves devoid of baryons, we find that baryonic effects are important. Compared to corresponding dark matter only simulations, the loss of baryons from subhaloes and stronger tidal disruption due to the presence of baryons near the centre of the main halo, reduce the number of subhaloes by ~ 1/4 to 1/2, independently of subhalo mass, but increasingly towards the host halo centre. We also find that subhaloes have non-Maxwellian orbital velocity distributions, with centrally rising velocity anisotropy and positive velocity bias which reduces the number of low-velocity subhaloes, particularly near the halo centre. We parameterise the predicted population of subhaloes in terms of mass, galactocentric distance, and velocities. We discuss implications of our results for the prospects of detecting dark matter substructures and for possible inferences about the nature of dark matter. △ Less

Submitted 6 September, 2016; originally announced September 2016.

Comments: 19 pages, submitted to MNRAS, comments welcome

The properties of "dark" ΛCDM halos in the Local Group

Authors: Alejandro Benítez-Llambay, Julio F. Navarro, Carlos S. Frenk, Till Sawala, Kyle Oman, Azadeh Fattahi, Matthieu Schaller, Joop Schaye, Robert A. Crain, Tom Theuns

Abstract: We examine the baryon content of low-mass ΛCDM halos $(10^8<M_{200}/{\rm M_\odot}<5\times 10^{9})$ using the APOSTLE cosmological hydrodynamical simulations. Most of these systems are free of stars and have a gaseous content set by the combined effects of cosmic reionization, which imposes a mass-dependent upper limit, and of ram pressure stripping, which reduces it further in high-density regions… ▽ More We examine the baryon content of low-mass ΛCDM halos $(10^8<M_{200}/{\rm M_\odot}<5\times 10^{9})$ using the APOSTLE cosmological hydrodynamical simulations. Most of these systems are free of stars and have a gaseous content set by the combined effects of cosmic reionization, which imposes a mass-dependent upper limit, and of ram pressure stripping, which reduces it further in high-density regions. Halos mainly affected by reionization (RELHICs; REionization-Limited HI Clouds) inhabit preferentially low-density regions and make up a population where the gas is in hydrostatic equilibrium with the dark matter potential and in thermal equilibrium with the ionizing UV background. Their thermodynamic properties are well specified, and their gas density and temperature profiles may be predicted in detail. Gas in RELHICs is nearly fully ionized but with neutral cores that span a large range of HI masses and column densities and have negligible non-thermal broadening. We present predictions for their characteristic sizes and central column densities: the massive tail of the distribution should be within reach of future blind HI surveys. Local Group RELHICs (LGRs) have some properties consistent with observed Ultra Compact High Velocity Clouds (UCHVCs) but the sheer number of the latter suggests that most UCHVCs are not RELHICs. Our results suggest that LGRs (i) should typically be beyond 500 kpc from the Milky Way or M31; (ii) have positive Galactocentric radial velocities; (iii) HI sizes not exceeding 1 kpc, and (iv) should be nearly round. The detection and characterization of RELHICs would offer a unique probe of the small-scale clustering of cold dark matter. △ Less

Submitted 17 November, 2016; v1 submitted 5 September, 2016; originally announced September 2016.

Comments: Accepted for publication in MNRAS. Referee's comments addressed in this version

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

Authors: Azadeh Fattahi, Julio F. Navarro, Till Sawala, Carlos S. Frenk, Laura V. Sales, Kyle Oman, Matthieu Schaller, Jie Wang

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

Submitted 21 July, 2016; originally announced July 2016.

Comments: 13 pages + 3 pages Appendix, 10 Figures, Submitted to MNRAS

What to expect from dynamical modelling of galactic haloes

Authors: Wenting Wang, Jiaxin Han, Shaun Cole, Carlos Frenk, Till Sawala

Abstract: Many dynamical models of the Milky Way halo require assumptions that the distribution function of a tracer population should be independent of time (i.e., a steady state distribution function) and that the underlying potential is spherical. We study the limitations of such modelling by applying a general dynamical model with minimal assumptions to a large sample of galactic haloes from cosmologica… ▽ More Many dynamical models of the Milky Way halo require assumptions that the distribution function of a tracer population should be independent of time (i.e., a steady state distribution function) and that the underlying potential is spherical. We study the limitations of such modelling by applying a general dynamical model with minimal assumptions to a large sample of galactic haloes from cosmological $N$-body and hydrodynamical simulations. Using dark matter particles as dynamical tracers, we find that the systematic uncertainties in the measured mass and concentration parameters typically have an amplitude of 25% to 40%. When stars are used as tracers, however, the systematic uncertainties can be as large as a factor of $2-3$. The systematic uncertainties are not reduced by increasing the tracer sample size and vary stochastically from halo to halo. These systematic uncertainties are mostly driven by underestimated statistical noise caused by correlated phase-space structures that violate the steady state assumption. The number of independent phase-space structures inferred from the uncertainty level sets a limiting sample size beyond which a further increase no longer significantly improves the accuracy of dynamical inferences. The systematic uncertainty level is determined by the halo merger history, the shape and environment of the halo. Our conclusions apply generally to any spherical steady-state model. △ Less

Submitted 27 July, 2017; v1 submitted 30 May, 2016; originally announced May 2016.

Comments: 16 pages, 12 figures, published by MNRAS

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

Authors: Matthieu Schaller, Carlos S. Frenk, Azadeh Fattahi, Julio F. Navarro, Kyle A. Oman, Till Sawala

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

Submitted 2 June, 2016; v1 submitted 9 May, 2016; originally announced May 2016.

Comments: 6 pages, 2 figures. Accepted for publication in MNRAS. Data available on request

Journal ref: MNRASL, 2016, Volume 461, Issue 1, pp. L56-L61

Knowing the unknowns: uncertainties in simple estimators of galactic dynamical masses

Authors: David J. R. Campbell, Carlos S. Frenk, Adrian Jenkins, Vincent R. Eke, Julio F. Navarro, Till Sawala, Matthieu Schaller, Azadeh Fattahi, Kyle A. Oman, Tom Theuns

Abstract: The observed stellar kinematics of dispersion-supported galaxies are often used to measure dynamical masses. Recently, several analytical relationships between the stellar line-of-sight velocity dispersion, the projected (2D) or deprojected (3D) half-light radius, and the total mass enclosed within the half-light radius, relying on the spherical Jeans equation, have been proposed. Here, we make us… ▽ More The observed stellar kinematics of dispersion-supported galaxies are often used to measure dynamical masses. Recently, several analytical relationships between the stellar line-of-sight velocity dispersion, the projected (2D) or deprojected (3D) half-light radius, and the total mass enclosed within the half-light radius, relying on the spherical Jeans equation, have been proposed. Here, we make use of the APOSTLE cosmological hydrodynamical simulations of the Local Group to test the validity and accuracy of such mass estimators for both dispersion and rotation-supported galaxies, for field and satellite galaxies, and for galaxies of varying masses, shapes, and velocity dispersion anisotropies. We find that the mass estimators of Walker et al. and Wolf et al. are able to recover the masses of dispersion-dominated systems with little systematic bias, but with a 1-sigma scatter of 25 and 23 percent, respectively. The error on the estimated mass is dominated by the impact of the 3D shape of the stellar mass distribution, which is difficult to constrain observationally. This intrinsic scatter becomes the dominant source of uncertainty in the masses estimated for galaxies like the dwarf spheroidal (dSph) satellites of the Milky Way, where the observational errors in their sizes and velocity dispersions are small. Such scatter may also affect the inner density slopes of dSphs derived from multiple stellar populations, relaxing the significance with which Navarro-Frenk-White profiles may be excluded, depending on the degree to which the relevant properties of the different stellar populations are correlated. Finally, we derive a new optimal mass estimator that removes the residual biases and achieves a statistically significant reduction in the scatter to 20 percent overall for dispersion-dominated galaxies, allowing more precise and accurate mass estimates. △ Less

Submitted 21 April, 2017; v1 submitted 14 March, 2016; originally announced March 2016.

Comments: Accepted for publication in MNRAS. 26 pages, 21 figures, 3 tables

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

Authors: Laura V. Sales, Julio F. Navarro, Kyle Oman, Azadeh Fattahi, Ismael Ferrero, Mario G. Abadi, Richard Bower, Robert A. Crain, Carlos S. Frenk, Till Sawala, Matthieu Schaller, Joop Schaye, Tom Theuns, Simon D. M. White

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

Submitted 5 February, 2016; originally announced February 2016.

Comments: 10 pages, 7 figures, submitted to MNRAS

Missing dark matter in dwarf galaxies?

Authors: Kyle A. Oman, Julio F. Navarro, Laura V. Sales, Azadeh Fattahi, Carlos S. Frenk, Till Sawala, Matthieu Schaller, Simon D. M. White

Abstract: We use cosmological hydrodynamical simulations of the APOSTLE project along with high-quality rotation curve observations to examine the fraction of baryons in ΛCDM haloes that collect into galaxies. This 'galaxy formation efficiency' correlates strongly and with little scatter with halo mass, dropping steadily towards dwarf galaxies. The baryonic mass of a galaxy may thus be used to place a lower… ▽ More We use cosmological hydrodynamical simulations of the APOSTLE project along with high-quality rotation curve observations to examine the fraction of baryons in ΛCDM haloes that collect into galaxies. This 'galaxy formation efficiency' correlates strongly and with little scatter with halo mass, dropping steadily towards dwarf galaxies. The baryonic mass of a galaxy may thus be used to place a lower limit on total halo mass and, consequently, on its asymptotic maximum circular velocity. A number of observed dwarfs seem to violate this constraint, having baryonic masses up to ten times higher than expected from their rotation speeds, or, alternatively, rotating at only half the speed expected for their mass. Taking the data at face value, either these systems have formed galaxies with extraordinary efficiency - highly unlikely given their shallow potential wells - or their dark matter content is much lower than expected from ΛCDM haloes. This 'missing dark matter' is reminiscent of the inner mass deficit of galaxies with slowly-rising rotation curves, but cannot be explained away by star formation-induced 'cores' in the dark mass profile, since the anomalous deficit applies to regions larger than the luminous galaxies themselves. We argue that explaining the structure of these galaxies would require either substantial modification of the standard Lambda cold dark matter paradigm or else significant revision to the uncertainties in their inferred mass profiles, which should be much larger than reported. Systematic errors in inclination may provide a simple resolution to what would otherwise be a rather intractable problem for the current paradigm. △ Less

Submitted 2 November, 2018; v1 submitted 5 January, 2016; originally announced January 2016.

Comments: 15 pages, 5 figures, 2 tables, MNRAS accepted version

Journal ref: MNRAS 460 (2016) 3610-3623

Dark matter-radiation interactions: the structure of Milky Way satellite galaxies

Authors: J. A. Schewtschenko, C. M. Baugh, R. J. Wilkinson, C. Boehm, S. Pascoli, T. Sawala

Abstract: In the thermal dark matter (DM) paradigm, primordial interactions between DM and Standard Model particles are responsible for the observed DM relic density. In Boehm et al. (2014), we showed that weak-strength interactions between DM and radiation (photons or neutrinos) can erase small-scale density fluctuations, leading to a suppression of the matter power spectrum compared to the collisionless c… ▽ More In the thermal dark matter (DM) paradigm, primordial interactions between DM and Standard Model particles are responsible for the observed DM relic density. In Boehm et al. (2014), we showed that weak-strength interactions between DM and radiation (photons or neutrinos) can erase small-scale density fluctuations, leading to a suppression of the matter power spectrum compared to the collisionless cold DM (CDM) model. This results in fewer DM subhaloes within Milky Way-like DM haloes, implying a reduction in the abundance of satellite galaxies. Here we use very high resolution N-body simulations to measure the dynamics of these subhaloes. We find that when interactions are included, the largest subhaloes are less concentrated than their counterparts in the collisionless CDM model and have rotation curves that match observational data, providing a new solution to the "too big to fail" problem. △ Less

Submitted 13 April, 2016; v1 submitted 21 December, 2015; originally announced December 2015.

Comments: 6 pages, 3 figures, 1 table. Augmented content available for figures (cf. 'Augmented Content' section). V3: Matches version in MNRAS

Report number: DURAST/2016/0001, IPPP/16/03

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

Authors: Till Sawala, Carlos S. Frenk, Azadeh Fattahi, Julio F. Navarro, Richard G. Bower, Robert A. Crain, Claudio Dalla Vecchia, Michelle Furlong, John. C. Helly, Adrian Jenkins, Kyle A. Oman, Matthieu Schaller, Joop Schaye, Tom Theuns, James Trayford, Simon D. M. White

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

Submitted 3 November, 2015; originally announced November 2015.

Comments: 15 pages, submitted to MNRAS. arXiv admin note: text overlap with arXiv:1412.2748

Dark matter annihilation radiation in hydrodynamic simulations of Milky Way haloes

Authors: Matthieu Schaller, Carlos S. Frenk, Tom Theuns, Francesca Calore, Gianfranco Bertone, Nassim Bozorgnia, Robert A. Crain, Azadeh Fattahi, Julio F. Navarro, Till Sawala, Joop Schaye

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

Submitted 25 November, 2015; v1 submitted 7 September, 2015; originally announced September 2015.

Comments: 11 pages, 5 figures, accepted for publication in MNRAS, corrected typos

Journal ref: MNRAS, 2015, Volume 455, Issue 4, pp. 4442-4451

The APOSTLE project: Local Group kinematic mass constraints and simulation candidate selection

Authors: Azadeh Fattahi, Julio F. Navarro, Till Sawala, Carlos S. Frenk, Kyle A. Oman, Robert A. Crain, Michelle Furlong, Matthieu Schaller, Joop Schaye, Tom Theuns, Adrian Jenkins

Abstract: We use a large sample of isolated dark matter halo pairs drawn from cosmological N-body simulations to identify candidate systems whose kinematics match that of the Local Group of Galaxies (LG). We find, in agreement with the "timing argument" and earlier work, that the separation and approach velocity of the Milky Way (MW) and Andromeda (M31) galaxies favour a total mass for the pair of… ▽ More We use a large sample of isolated dark matter halo pairs drawn from cosmological N-body simulations to identify candidate systems whose kinematics match that of the Local Group of Galaxies (LG). We find, in agreement with the "timing argument" and earlier work, that the separation and approach velocity of the Milky Way (MW) and Andromeda (M31) galaxies favour a total mass for the pair of $\sim 5\times 10^{12} \,M_{\odot}$. A mass this large, however, is difficult to reconcile with the small relative tangential velocity of the pair, as well as with the small deceleration from the Hubble flow observed for the most distant LG members. Halo pairs that match these three criteria have average masses a factor of $\sim 2$ times smaller than suggested by the timing argument, but with large dispersion. Guided by these results, we have selected $12$ halo pairs with total mass in the range $1.6$-$3.6 \times 10^{12}\,M_{\odot}$ for the APOSTLE project (A Project Of Simulating The Local Environment), a suite of hydrodynamical resimulations at various numerical resolution levels (reaching up to $\sim10^{4}\,M_{\odot}$ per gas particle) that use the subgrid physics developed for the EAGLE project. These simulations reproduce, by construction, the main kinematics of the MW-M31 pair, and produce satellite populations whose overall number, luminosities, and kinematics are in good agreement with observations of the MW and M31 companions. The APOSTLE candidate systems thus provide an excellent testbed to confront directly many of the predictions of the $Λ$CDM cosmology with observations of our local Universe. △ Less

Submitted 7 October, 2020; v1 submitted 13 July, 2015; originally announced July 2015.

Comments: replaced with the accepted version

Planes of satellite galaxies: when exceptions are the rule

Authors: Marius Cautun, Sownak Bose, Carlos S. Frenk, Qi Guo, Jiaxin Han, Wojciech A. Hellwing, Till Sawala, Wenting Wang

Abstract: The detection of planar structures within the satellite systems of both the Milky Way (MW) and Andromeda (M31) has been reported as being in stark contradiction to the predictions of the standard cosmological model ($Λ$CDM). Given the ambiguity in defining a planar configuration, it is unclear how to interpret the low incidence of the MW and M31 planes in $Λ$CDM. We investigate the prevalence of s… ▽ More The detection of planar structures within the satellite systems of both the Milky Way (MW) and Andromeda (M31) has been reported as being in stark contradiction to the predictions of the standard cosmological model ($Λ$CDM). Given the ambiguity in defining a planar configuration, it is unclear how to interpret the low incidence of the MW and M31 planes in $Λ$CDM. We investigate the prevalence of satellite planes around galactic mass haloes identified in high resolution cosmological simulations. We find that planar structures are very common, and that ~10% of $Λ$CDM haloes have even more prominent planes than those present in the Local Group. While ubiquitous, the planes of satellite galaxies show a large diversity in their properties. This precludes using one or two systems as small scale probes of cosmology, since a large sample of satellite systems is needed to obtain a good measure of the object-to-object variation. This very diversity has been misinterpreted as a discrepancy between the satellite planes observed in the Local Group and $Λ$CDM predictions. In fact, ~10% of $Λ$CDM galactic haloes have planes of satellites that are as infrequent as the MW and M31 planes. The look-elsewhere effect plays an important role in assessing the detection significance of satellite planes and accounting for it leads to overestimating the significance level by a factor of 30 and 100 for the MW and M31 systems, respectively. △ Less

Submitted 7 August, 2015; v1 submitted 12 June, 2015; originally announced June 2015.

Comments: 13 pages, 15 figures, +3 pages appendix, MNRAS accepted version

Journal ref: MNRAS, 2015, 452 (4), 3838-3852

The Copernicus Complexio: a high-resolution view of the small-scale Universe

Authors: Wojciech A. Hellwing, Carlos S. Frenk, Marius Cautun, Sownak Bose, John Helly, Adrian Jenkins, Till Sawala, Maciej Cytowski

Abstract: We introduce Copernicus Complexio (COCO), a high-resolution cosmological N-body simulation of structure formation in the $Λ{\rm CDM}{}$ model. COCO follows an approximately spherical region of radius $\sim 17.4h^{-1}\,{\rm Mpc}$ embedded in a much larger periodic cube that is followed at lower resolution. The high resolution volume has a particle mass of $1.135\times10^5h^{-1}{\rm M}_{\odot}$ (60… ▽ More We introduce Copernicus Complexio (COCO), a high-resolution cosmological N-body simulation of structure formation in the $Λ{\rm CDM}{}$ model. COCO follows an approximately spherical region of radius $\sim 17.4h^{-1}\,{\rm Mpc}$ embedded in a much larger periodic cube that is followed at lower resolution. The high resolution volume has a particle mass of $1.135\times10^5h^{-1}{\rm M}_{\odot}$ (60 times higher than the Millennium-II simulation). COCO gives the dark matter halo mass function over eight orders of magnitude in halo mass; it forms $\sim 60$ haloes of galactic size, each resolved with about 10 million particles. We confirm the power-law character of the subhalo mass function, $\bar{N}(>μ)\proptoμ^{-s}$, down to a reduced subhalo mass $M_{sub}/M_{200}\equivμ=10^{-6}$, with a best-fit power-law index, $s=0.94$, for hosts of mass $\langle M_{200}\rangle=10^{12}h^{-1}{\rm M}_{\odot}$. The concentration-mass relation of COCO haloes deviates from a single power law for masses $M_{200}<\textrm{a few}\times 10^{8}h^{-1}{\rm M}_{\odot}$, where it flattens, in agreement with results by Sanchez-Conde et al. The host mass invariance of the reduced maximum circular velocity function of subhaloes, $ν\equiv V_{max}/V_{200}$, hinted at in previous simulations, is clearly demonstrated over five orders of magnitude in host mass. Similarly, we find that the average, normalised radial distribution of subhaloes is approximately universal (i.e. independent of subhalo mass), as previously suggested by the Aquarius simulations of individual haloes. Finally, we find that at fixed physical subhalo size, subhaloes in lower mass hosts typically have lower central densities than those in higher mass hosts. △ Less

Submitted 11 February, 2016; v1 submitted 24 May, 2015; originally announced May 2015.

Comments: 18 pages, 15 figures, 3 tables, version accepted for publication in MNRAS, figure 7 updated with a corrected data, minor changes to the text, conclusions unchanged

The unexpected diversity of dwarf galaxy rotation curves

Authors: Kyle A. Oman, Julio F. Navarro, Azadeh Fattahi, Carlos S. Frenk, Till Sawala, Simon D. M. White, Richard Bower, Robert A. Crain, Michelle Furlong, Matthieu Schaller, Joop Schaye, Tom Theuns

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

Submitted 8 July, 2015; v1 submitted 6 April, 2015; originally announced April 2015.

Comments: 17 pages, 6 figures, 2 tables, MNRAS accepted

Journal ref: MNRAS 452 (2015) 3650-3665

Local Group galaxies emerge from the dark

Authors: Till Sawala, Carlos S. Frenk, Azadeh Fattahi, Julio F. Navarro, Richard G. Bower, Robert A. Crain, Claudio Dalla Vecchia, Michelle Furlong, John. C. Helly, Adrian Jenkins, Kyle A. Oman, Matthieu Schaller, Joop Schaye, Tom Theuns, James Trayford, Simon D. M. White

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

Submitted 8 December, 2014; originally announced December 2014.

Comments: 14 pages, 4 figures

A new spin on disks of satellite galaxies

Authors: Marius Cautun, Wenting Wang, Carlos S. Frenk, Till Sawala

Abstract: We investigate the angular and kinematic distributions of satellite galaxies around a large sample of bright isolated primaries in the spectroscopic and photometric catalogues of the Sloan Digital Sky Survey (SDSS). We detect significant anisotropy in the spatial distribution of satellites. To test whether this anisotropy could be related to the rotating disks of satellites recently found by Ibata… ▽ More We investigate the angular and kinematic distributions of satellite galaxies around a large sample of bright isolated primaries in the spectroscopic and photometric catalogues of the Sloan Digital Sky Survey (SDSS). We detect significant anisotropy in the spatial distribution of satellites. To test whether this anisotropy could be related to the rotating disks of satellites recently found by Ibata et al. in a sample of SDSS galaxies, we repeat and extend their analysis. Ibata et al. found an excess of satellites on opposite sides of their primaries having anticorrelated radial velocities. We find that this excess is sensitive to small changes in the sample selection criteria which can greatly reduce its significance. In addition, we find no evidence for correspondingly correlated velocities for satellites observed on the same side of their primaries, which would be expected for rotating disks of satellites. We conclude that the detection of coherent rotation in the satellite population in current observational samples is not robust. We compare our data to the $Λ$CDM Millennium simulations populated with galaxies according to the semi-analytic model of Guo et al. We find excellent agreement with the spatial distribution of satellites in the SDSS data and the lack of a strong signal from coherent rotation. △ Less

Submitted 13 April, 2015; v1 submitted 28 October, 2014; originally announced October 2014.

Comments: 13 pages, 12 figures, accepted version, MNRAS

Journal ref: MNRAS, 2015, Volume 449, Issue 3, p.2576-2587

Orbital parameters of infalling satellite haloes in the hierarchical $Λ$CDM model

Authors: Lilian Jiang, Shaun Cole, Till Sawala, Carlos S. Frenk

Abstract: We present distributions of orbital parameters of infalling satellites of $Λ$CDM haloes in the mass range $10^{12}-10^{14}$M$_\odot$, which represent the initial conditions for the subsequent evolution of substructures within the host halo. We use merger trees constructed in a high resolution cosmological N-body simulation to trace satellite haloes, and identify the time of infall. We find signfic… ▽ More We present distributions of orbital parameters of infalling satellites of $Λ$CDM haloes in the mass range $10^{12}-10^{14}$M$_\odot$, which represent the initial conditions for the subsequent evolution of substructures within the host halo. We use merger trees constructed in a high resolution cosmological N-body simulation to trace satellite haloes, and identify the time of infall. We find signficant trends in the distribution of orbital parameters with both the host halo mass and the ratio of satellite-to-host halo masses. For all host halo masses, satellites whose infall mass is a larger fraction of the host halo mass have more eccentric, radially biased orbits. At fixed satellite-to-host halo mass ratio, high mass haloes are biased towards accreting satellites on slightly more radial orbits. To charactise the orbital distributions fully requires fitting the correlated bivariate distribution of two chosen orbital parameters (e.g. radial and tangential velocity or energy and angular momentum). We provide simple fits to one choice of the bivariate distributions, which when transformed faithfully, captures the behaviour of any of the projected one-dimensional distributions. △ Less

Submitted 27 January, 2015; v1 submitted 3 September, 2014; originally announced September 2014.

Comments: 14 pages, 13 figures, accepted for publication in MNRAS

The EAGLE project: Simulating the evolution and assembly of galaxies and their environments

Authors: Joop Schaye, Robert A. Crain, Richard G. Bower, Michelle Furlong, Matthieu Schaller, Tom Theuns, Claudio Dalla Vecchia, Carlos S. Frenk, I. G. McCarthy, John C. Helly, Adrian Jenkins, Y. M. Rosas-Guevara, Simon D. M. White, Maarten Baes, C. M. Booth, Peter Camps, Julio F. Navarro, Yan Qu, Alireza Rahmati, Till Sawala, Peter A. Thomas, James Trayford

Abstract: We introduce the Virgo Consortium's EAGLE project, a suite of hydrodynamical simulations that follow the formation of galaxies and black holes in representative volumes. We discuss the limitations of such simulations in light of their finite resolution and poorly constrained subgrid physics, and how these affect their predictive power. One major improvement is our treatment of feedback from massiv… ▽ More We introduce the Virgo Consortium's EAGLE project, a suite of hydrodynamical simulations that follow the formation of galaxies and black holes in representative volumes. We discuss the limitations of such simulations in light of their finite resolution and poorly constrained subgrid physics, and how these affect their predictive power. One major improvement is our treatment of feedback from massive stars and AGN in which thermal energy is injected into the gas without the need to turn off cooling or hydrodynamical forces, allowing winds to develop without predetermined speed or mass loading factors. Because the feedback efficiencies cannot be predicted from first principles, we calibrate them to the z~0 galaxy stellar mass function and the amplitude of the galaxy-central black hole mass relation, also taking galaxy sizes into account. The observed galaxy mass function is reproduced to $\lesssim 0.2$ dex over the full mass range, $10^8 < M_*/M_\odot \lesssim 10^{11}$, a level of agreement close to that attained by semi-analytic models, and unprecedented for hydrodynamical simulations. We compare our results to a representative set of low-redshift observables not considered in the calibration, and find good agreement with the observed galaxy specific star formation rates, passive fractions, Tully-Fisher relation, total stellar luminosities of galaxy clusters, and column density distributions of intergalactic CIV and OVI. While the mass-metallicity relations for gas and stars are consistent with observations for $M_* \gtrsim 10^9 M_\odot$, they are insufficiently steep at lower masses. The gas fractions and temperatures are too high for clusters of galaxies, but for groups these discrepancies can be resolved by adopting a higher heating temperature in the subgrid prescription for AGN feedback. EAGLE constitutes a valuable new resource for studies of galaxy formation. △ Less

Submitted 2 October, 2014; v1 submitted 25 July, 2014; originally announced July 2014.

Comments: Accepted for publication in MNRAS. V2: Minor changes. For images and videos, see http://eagle.strw.leidenuniv.nl/ and http://icc.dur.ac.uk/Eagle/

Journal ref: 2015, MNRAS, 446, 521

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

Authors: Till Sawala, Carlos S. Frenk, Azadeh Fattahi, Julio F. Navarro, Tom Theuns, Richard G. Bower, Robert A. Crain, Michelle Furlong, Adrian Jenkins, Matthieu Schaller, Joop Schaye

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

Submitted 24 June, 2014; originally announced June 2014.

Comments: 12 pages, 12 figures, submitted to MNRAS, comments welcome

Subhalo statistics of galactic halos: beyond the resolution limit

Authors: Marius Cautun, Wojciech A. Hellwing, Rien van de Weygaert, Carlos S. Frenk, Bernard J. T. Jones, Till Sawala

Abstract: We study the substructure population of Milky Way (MW)-mass halos in the $Λ$CDM cosmology using a novel procedure to extrapolate subhalo number statistics beyond the resolution limit of N-body simulations. The technique recovers the mean and the variance of the subhalo abundance, but not its spatial distribution. It extends the dynamic range over which precise statistical predictions can be made b… ▽ More We study the substructure population of Milky Way (MW)-mass halos in the $Λ$CDM cosmology using a novel procedure to extrapolate subhalo number statistics beyond the resolution limit of N-body simulations. The technique recovers the mean and the variance of the subhalo abundance, but not its spatial distribution. It extends the dynamic range over which precise statistical predictions can be made by the equivalent of performing a simulation with 50 times higher resolution, at no additional computational cost. We apply this technique to MW-mass halos, but it can easily be applied to halos of any mass. We find up to $20\%$ more substructures in MW-mass halos than found in previous studies. Our analysis lowers the mass of the MW halo required to accommodate the observation that the MW has only three satellites with a maximum circular velocity $V_{max}\ge30 km/s$ in the $Λ$CDM cosmology. The probability of having a subhalo population similar to that in the MW is $20\%$ for a virial mass, $M_{200}=1\times10^{12} M_\odot$ and practically zero for halos more massive than $M_{200}=2\times10^{12} M_\odot$. △ Less

Submitted 23 October, 2014; v1 submitted 29 May, 2014; originally announced May 2014.

Comments: 17 pages, 14 figures, updated to match published version in MNRAS

Journal ref: 2014, MNRAS, 445, 1820

Bent by baryons: the low mass galaxy-halo relation

Authors: Till Sawala, Carlos S. Frenk, Azadeh Fattahi, Julio F. Navarro, Richard G. Bower, Robert A. Crain, Claudio Dalla Vecchia, Michelle Furlong, Adrian Jenkins, Ian G. McCarthy, Yan Qu, Matthieu Schaller, Joop Schaye, Tom Theuns

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

Submitted 14 April, 2014; originally announced April 2014.

Comments: 8 pages, submitted to MNRAS, comments welcome

The abundance of (not just) dark matter haloes

Authors: Till Sawala, Carlos S. Frenk, Robert A. Crain, Adrian Jenkins, Joop Schaye, Tom Theuns, Jesus Zavala

Abstract: We study the effect of baryons on the abundance of structures and substructures in a Lambda-CDM cosmology, using a pair of high resolution cosmological simulations from the GIMIC project. Both simulations use identical initial conditions, but while one contains only dark matter, the other also includes baryons. We find that gas pressure, reionisation, supernova feedback, stripping, and truncated a… ▽ More We study the effect of baryons on the abundance of structures and substructures in a Lambda-CDM cosmology, using a pair of high resolution cosmological simulations from the GIMIC project. Both simulations use identical initial conditions, but while one contains only dark matter, the other also includes baryons. We find that gas pressure, reionisation, supernova feedback, stripping, and truncated accretion systematically reduce the total mass and the abundance of structures below ~10^12 solar masses compared to the pure dark matter simulation. Taking this into account and adopting an appropriate detection threshold lowers the abundance of observed galaxies with maximum circular velocities below 100 km/s, significantly reducing the reported discrepancy between Lambda-CDM and the measured HI velocity function of the ALFALFA survey. We also show that the stellar-to-total mass ratios of galaxies with stellar masses of ~10^5 - 10^7 solar masses inferred from abundance matching of the (sub)halo mass function to the observed galaxy mass function increase by a factor of ~2. In addition, we find that an important fraction of low-mass subhaloes are completely devoid of stars. Accounting for the presence of dark subhaloes below 10^10 solar masses further reduces the abundance of observable objects, and leads to an additional increase in the inferred stellar-to-total mass ratio by factors of 2 - 10 for galaxies in haloes of 10^9 - 10^10 solar masses. This largely reconciles the abundance matching results with the kinematics of individual dwarf galaxies in Lambda-CDM. We propose approximate corrections to the masses of objects derived from pure dark matter calculations to account for baryonic effects. △ Less

Submitted 24 April, 2013; v1 submitted 27 June, 2012; originally announced June 2012.

Comments: 19 pages, updated version published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, 05/2013, Volume 431, Issue 2, pp. 1366-1382

The Missing Massive Satellites of the Milky Way

Authors: Jie Wang, Carlos S. Frenk, Julio F. Navarro, Liang Gao, Till Sawala

Abstract: Recent studies suggest that only three of the twelve brightest satellites of the Milky Way (MW) inhabit dark matter halos with maximum circular velocity, V_max, exceeding 30km/s. This is in apparent contradiction with the LCDM simulations of the Aquarius Project, which suggest that MW-sized halos should have at least 8 subhalos with V_max>30km/s. The absence of luminous satellites in such massive… ▽ More Recent studies suggest that only three of the twelve brightest satellites of the Milky Way (MW) inhabit dark matter halos with maximum circular velocity, V_max, exceeding 30km/s. This is in apparent contradiction with the LCDM simulations of the Aquarius Project, which suggest that MW-sized halos should have at least 8 subhalos with V_max>30km/s. The absence of luminous satellites in such massive subhalos is thus puzzling and may present a challenge to the LCDM paradigm. We note, however, that the number of massive subhalos depends sensitively on the (poorly-known) virial mass of the Milky Way, and that their scarcity makes estimates of their abundance from a small simulation set like Aquarius uncertain. We use the Millennium Simulation series and the invariance of the scaled subhalo velocity function (i.e., the number of subhalos as a function of ν, the ratio of subhalo V_max to host halo virial velocity, V_200) to secure improved estimates of the abundance of rare massive subsystems. In the range 0.1<ν<0.5, N_sub(>ν) is approximately Poisson-distributed about an average given by <N_sub>=10.2x(ν/0.15)^(-3.11). This is slightly lower than in Aquarius halos, but consistent with recent results from the Phoenix Project. The probability that a LCDM halo has 3 or fewer subhalos with V_max above some threshold value, V_th, is then straightforward to compute. It decreases steeply both with decreasing V_th and with increasing halo mass. For V_th=30km/s, ~40% of M_halo=10^12 M_sun halos pass the test; fewer than 5% do so for M_halo>= 2x10^12 M_sun; and the probability effectively vanishes for M_halo>= 3x 10^12 M_sun. Rather than a failure of LCDM, the absence of massive subhalos might simply indicate that the Milky Way is less massive than is commonly thought. △ Less

Submitted 30 May, 2012; v1 submitted 19 March, 2012; originally announced March 2012.

Comments: 5 figures, 8 pages. Accepted for publication in MNRAS. The number of massive satellites within a bigger radius r100 and in a WMAP7 cosmology model is considered now

Local Group Dwarf Galaxies: Nature And Nurture

Authors: Till Sawala, Cecilia Scannapieco, Simon White

Abstract: We investigate the formation and evolution of dwarf galaxies in a high resolution, hydrodynamical cosmological simulation of a Milky Way sized halo and its environment. Our simulation includes gas cooling, star formation, supernova feedback, metal enrichment and UV heating. In total, 90 satellites and more than 400 isolated dwarf galaxies are formed in the simulation, allowing a systematic study o… ▽ More We investigate the formation and evolution of dwarf galaxies in a high resolution, hydrodynamical cosmological simulation of a Milky Way sized halo and its environment. Our simulation includes gas cooling, star formation, supernova feedback, metal enrichment and UV heating. In total, 90 satellites and more than 400 isolated dwarf galaxies are formed in the simulation, allowing a systematic study of the internal and environmental processes that determine their evolution. We find that 95% of satellite galaxies are gas-free at z=0, and identify three mechanisms for gas loss: supernova feedback, tidal stripping, and photo-evaporation due to re-ionization. Gas-rich satellite galaxies are only found with total masses above ~ 5x10^9 solar masses. In contrast, for isolated dwarf galaxies, a total mass of ~ 10^9 solar masses constitutes a sharp transition; less massive galaxies are predominantly gas-free at z=0, more massive, isolated dwarf galaxies are often able to retain their gas. In general, we find that the total mass of a dwarf galaxy is the main factor which determines its star formation, metal enrichment, and its gas content, but that stripping may explain the observed difference in gas content between field dwarf galaxies and satellites with total masses close to 10^9 solar masses. We also find that a morphological transformation via tidal stripping of infalling, luminous dwarf galaxies whose dark matter is less concentrated than their stars, cannot explain the high total mass-light ratios of the faint dwarf spheroidal galaxies. △ Less

Submitted 23 March, 2011; originally announced March 2011.

Comments: 18 pages, 13 figures, submitted to MNRAS

What is the (Dark) Matter with Dwarf Galaxies?

Authors: Till Sawala, Qi Guo, Cecilia Scannapieco, Adrian Jenkins, Simon D. M. White

Abstract: We present cosmological hydrodynamical simulations of the formation of dwarf galaxies in a representative sample of haloes extracted from the Millennium-II Simulation. Our six haloes have a z = 0 mass of ~10^10 solar masses and show different mass assembly histories which are reflected in different star formation histories. We find final stellar masses in the range 5 x 10^7 - 10^8 solar masses, co… ▽ More We present cosmological hydrodynamical simulations of the formation of dwarf galaxies in a representative sample of haloes extracted from the Millennium-II Simulation. Our six haloes have a z = 0 mass of ~10^10 solar masses and show different mass assembly histories which are reflected in different star formation histories. We find final stellar masses in the range 5 x 10^7 - 10^8 solar masses, consistent with other published simulations of galaxy formation in similar mass haloes. Our final objects have structures and stellar populations consistent with dwarf elliptical and dwarf irregular galaxies. However, in a Lambda CDM universe, 10^10 solar mass haloes must typically contain galaxies with much lower stellar mass than our simulated objects if they are to match observed galaxy abundances. The dwarf galaxies formed in our own and all other current hydrodynamical simulations are more than an order of magnitude more luminous than expected for haloes of this mass. We discuss the significance and possible implications of this result. △ Less

Submitted 7 December, 2010; v1 submitted 3 March, 2010; originally announced March 2010.

Comments: 12 pages, 4 figures, MNRAS accepted

The Dynamical and Chemical Evolution of Dwarf Spheroidal Galaxies

Authors: Y. Revaz, P. Jablonka, T. Sawala, V. Hill, B. Letarte, M. Irwin, G. Battaglia, A. Helmi, M. D. Shetrone, E. Tolstoy, K. A. Venn

Abstract: We present a large sample of fully self-consistent hydrodynamical Nbody/Tree-SPH simulations of isolated dwarf spheroidal galaxies (dSphs). It has enabled us to identify the key physical parameters and mechanisms at the origin of the observed variety in the Local Group dSph properties. The initial total mass (gas + dark matter) of these galaxies is the main driver of their evolution. Star format… ▽ More We present a large sample of fully self-consistent hydrodynamical Nbody/Tree-SPH simulations of isolated dwarf spheroidal galaxies (dSphs). It has enabled us to identify the key physical parameters and mechanisms at the origin of the observed variety in the Local Group dSph properties. The initial total mass (gas + dark matter) of these galaxies is the main driver of their evolution. Star formation (SF) occurs in series of short bursts. In massive systems, the very short intervals between the SF peaks mimic a continuous star formation rate, while less massive systems exhibit well separated SF bursts, as identified observationally. The delay between the SF events is controlled by the gas cooling time dependence on galaxy mass. The observed global scaling relations, luminosity-mass and luminosity-metallicity, are reproduced with low scatter. We take advantage of the unprecedentedly large sample size and data homogeneity of the ESO Large Programme DART, and add to it a few independent studies, to constrain the star formation history of five Milky Way dSphs, Sextans, LeoII, Carina, Sculptor and Fornax. For the first time, [Mg/Fe] vs [Fe/H] diagrams derived from high-resolution spectroscopy of hundreds of individual stars are confronted with model predictions. We find that the diversity in dSph properties may well result from intrinsic evolution. We note, however, that the presence of gas in the final state of our simulations, of the order of what is observed in dwarf irregulars, calls for removal by external processes. △ Less

Submitted 29 April, 2009; originally announced April 2009.

Comments: 21 Pages, 19 figures ; Accepted for publication in A&A. Higher resolution version may be downloaded here : http://obswww.unige.ch/~revaz/publications/aa2009_11734

Formation of Isolated Dwarf Galaxies with Feedback

Authors: Till Sawala, Cecilia Scannapieco, Umberto Maio, Simon D. M. White

Abstract: We present results of high resolution hydrodynamical simulations of the formation and evolution of dwarf galaxies. Our simulations start from cosmological initial conditions at high redshift. They include metal-dependent cooling, star formation, feedback from type II and type Ia supernovae and UV background radiation, with physical recipes identical to those applied in a previous study of Milky… ▽ More We present results of high resolution hydrodynamical simulations of the formation and evolution of dwarf galaxies. Our simulations start from cosmological initial conditions at high redshift. They include metal-dependent cooling, star formation, feedback from type II and type Ia supernovae and UV background radiation, with physical recipes identical to those applied in a previous study of Milky Way type galaxies. We find that a combination of feedback and the cosmic UV background results in the formation of galaxies with properties similar to the Local Group dwarf spheroidals, and that their effect is strongly moderated by the depth of the gravitational potential. Taking this into account, our models naturally reproduce the observed luminosities and metallicities. The final objects have halo masses between 2.3x10^8 and 1.1x10^9 solar masses, mean velocity dispersions between 6.5 and 9.7 kms-1, stellar masses ranging from 5x10^5 to 1.2x10^7 solar masses, median metallicities between [Fe/H] = -1.8 and -1.1, and half-light radii of the order of 200 to 300 pc, all comparable with Local Group dwarf spheroidals. Our simulations also indicate that the dwarf spheroidal galaxies observed today lie near a halo mass threshold around 10^9 solar masses, in agreement with stellar kinematic data, where supernova feedback not only suffices to completely expel the interstellar medium and leave the residual gas-free, but where the combination of feedback, UV radiation and self-shielding establishes a dichotomy of age distributions similar to that observed in the Milky Way and M31 satellites. △ Less

Submitted 23 November, 2009; v1 submitted 10 February, 2009; originally announced February 2009.

Comments: 17 pages, 14 figures. MNRAS accepted