Showing 1–50 of 56 results for author: Ludlow, A

Milky Way-mass disc galaxies with low-mass stellar haloes have diverse merger histories

Authors: Katy L. Proctor, Aaron D. Ludlow, Claudia del P. Lagos, Aaron S. G. Robotham

Abstract: We use the Eagle simulation to study the relationship between the stellar haloes of Milky-Way-mass galaxies and their recent merger histories. The stellar mass ratio of the most massive merger that occurred since $z=1$ is strongly correlated with the $z=0$ fraction of ex situ stars in the galaxy, $f_{\rm ex\ situ}$, but is weakly correlated with stellar halo mass fraction, $f_{\rm SH}$, particular… ▽ More We use the Eagle simulation to study the relationship between the stellar haloes of Milky-Way-mass galaxies and their recent merger histories. The stellar mass ratio of the most massive merger that occurred since $z=1$ is strongly correlated with the $z=0$ fraction of ex situ stars in the galaxy, $f_{\rm ex\ situ}$, but is weakly correlated with stellar halo mass fraction, $f_{\rm SH}$, particularly for disc galaxies. Contrary to common belief, our results suggest that disc galaxies with low mass stellar haloes do not necessarily have quiescent merger histories; in fact, roughly one quarter have experienced a merger with stellar mass ratio $> 0.1$ since $z=1$. We demonstrate that the population of disc galaxies with low $f_{\rm SH}$ and active merger histories undergo mergers with satellites whose orbits are more circular than average and are approximately co-planar with the disc; instead of contributing significantly to the stellar halo, these mergers lead to discs that contain a substantial fraction of ex situ stellar mass and are thicker and more extended than those of quiescent galaxies. Such mergers also supply fuel that often incites a significant episode of in situ star formation in the disc. Our results suggest that seemingly quiescent disc galaxies with low-mass stellar haloes actually have diverse merger histories, which limits the extent to which stellar haloes alone can act as observable tracers of the hierarchical assembly of galaxies. △ Less

Submitted 16 July, 2024; originally announced July 2024.

Comments: 11 pages, 8 figures. Submitted to MNRAS

The distribution of stellar orbits in Eagle galaxies -- the effect of mergers, gas accretion, and secular evolution

Authors: Giulia Santucci, Claudia Del P. Lagos, Katherine E. Harborne, Aaron Ludlow, Caro Foster, Richard McDermid, Adriano Poci, Katy L. Proctor, Sabine Thater, Glenn van de Ven, Ling Zhu, Daniel Walo Martin

Abstract: The merger history of a galaxy is thought to be one of the major factors determining its internal dynamics, with galaxies having undergone different types or mergers (e.g. dry, minor or major mergers) predicted to show different dynamical properties. We study the instantaneous orbital distribution of galaxies in the Eagle simulation, colouring the orbits of the stellar particles by their stellar a… ▽ More The merger history of a galaxy is thought to be one of the major factors determining its internal dynamics, with galaxies having undergone different types or mergers (e.g. dry, minor or major mergers) predicted to show different dynamical properties. We study the instantaneous orbital distribution of galaxies in the Eagle simulation, colouring the orbits of the stellar particles by their stellar age, in order to understand whether stars form in particular orbits (e.g. in a thin or thick disc). We first show that Eagle reproduces well the observed stellar mass fractions in different stellar orbital families as a function of stellar mass and spin parameter at z = 0. We find that the youngest stars reside in a thin disc component that can extend to the very inner regions of galaxies, and that older stars have warmer orbits, with the oldest ones showing orbits consistent with both hot and counter-rotating classifications, which is consistent with the trend found in the Milky-Way and other disc galaxies. We also show that counter-rotating orbits trace galaxy mergers - in particular dry mergers, and that in the absence of mergers, counter-rotating orbits can also be born from highly misaligned gas accretion that leads to star formation. △ Less

Submitted 10 October, 2023; originally announced October 2023.

Comments: Submitted to MNRAS, 21 pages

Identifying the discs, bulges, and intra-halo light of simulated galaxies through structural decomposition

Authors: Katy L. Proctor, Claudia del P. Lagos, Aaron D. Ludlow, Aaron S. G. Robotham

Abstract: We perform a structural decomposition of galaxies identified in three cosmological hydrodynamical simulations by applying Gaussian Mixture Models (GMMs) to the kinematics of their stellar particles. We study the resulting disc, bulge, and intra-halo light (IHL) components of galaxies whose host dark matter haloes have virial masses in the range $M_{200}=10^{11}$-- $10^{15}\,{\rm M_\odot}$. Our dec… ▽ More We perform a structural decomposition of galaxies identified in three cosmological hydrodynamical simulations by applying Gaussian Mixture Models (GMMs) to the kinematics of their stellar particles. We study the resulting disc, bulge, and intra-halo light (IHL) components of galaxies whose host dark matter haloes have virial masses in the range $M_{200}=10^{11}$-- $10^{15}\,{\rm M_\odot}$. Our decomposition technique isolates galactic discs whose mass fractions, $f_{\rm disc}$, correlate strongly with common alternative morphology indicators; for example, $f_{\rm disc}$ is approximately equal to $κ_{\rm co}$, the fraction of stellar kinetic energy in co-rotation. The primary aim of our study, however, is to characterise the IHL of galaxies in a consistent manner and over a broad mass range, and to analyse its properties from the scale of galactic stellar haloes up to the intra-cluster light. Our results imply that the IHL fraction, $f_{\rm IHL}$, has appreciable scatter and is strongly correlated with galaxy morphology: at fixed stellar mass, the IHL of disc galaxies is typically older and less massive than that of spheroids. Above $M_{200}\approx 10^{13}\,{\rm M_\odot}$, we find, on average, $f_{\rm IHL}\approx 0.45$, albeit with considerable scatter. The transition radius beyond which the IHL dominates the stellar mass of a galaxy is roughly $30\,{\rm kpc}$ for $M_{200}\lesssim 10^{12.8}\,{\rm M_\odot}$, but increases strongly towards higher masses. However, we find that no alternative IHL definitions -- whether based on the ex-situ stellar mass, or the stellar mass outside a spherical aperture -- reproduce our dynamically-defined IHL masses. △ Less

Submitted 6 December, 2023; v1 submitted 24 July, 2023; originally announced July 2023.

Comments: 16 pages, 11 figures. MNRAS published

Spurious heating of stellar motions by dark matter particles in cosmological simulations of galaxy formation

Authors: Aaron D. Ludlow, S. Michael Fall, Matthew J. Wilkinson, Joop Schaye, Danail Obreschkow

Abstract: We use two cosmological simulations to study the impact of spurious heating of stellar motions within simulated galaxies by dark matter (DM) particles. The simulations share the same numerical and subgrid parameters, but one used a factor of 7 more DM particles. Many galaxy properties are unaffected by spurious heating, including their masses, star formation histories, and the spatial distribution… ▽ More We use two cosmological simulations to study the impact of spurious heating of stellar motions within simulated galaxies by dark matter (DM) particles. The simulations share the same numerical and subgrid parameters, but one used a factor of 7 more DM particles. Many galaxy properties are unaffected by spurious heating, including their masses, star formation histories, and the spatial distribution of their gaseous baryons. The distribution and kinematics of stellar and DM particles, however, are affected. Below a resolution-dependent virial mass, $M_{200}^{\rm spur}$, galaxies have higher characteristic velocities, larger sizes, and more angular momentum in the simulation with lower DM mass resolution; haloes have higher central densities and lower velocity dispersions. Above $M_{200}^{\rm spur}$, galaxies and haloes have similar properties in both runs. The differences arise due to spurious heating, which transfers energy from DM to stellar particles, causing galaxies to heat up and haloes to cool down. The value of $M_{200}^{\rm spur}$ can be derived from an empirical disc heating model, and coincides with the mass below which the predicted {\em spurious} velocity dispersion exceeds the {\em measured} velocity dispersion of simulated galaxies. We predict that galaxies in the $100^3\, {\rm Mpc}^3$ \eagle\, run and IllustrisTNG-100 are robust to spurious collisional effects at their half-mass radii provided $M_{200}^{\rm spur}\approx 10^{11.7}{\rm M_\odot}$; for the $25^3\, {\rm Mpc}^3$ \eagle\, run and IllustrisTNG-50, we predict $M_{200}^{\rm spur}\approx 10^{11}{\rm M_\odot}$. Suppressing spurious heating at smaller/larger radii, or for older/younger stellar populations, requires haloes to be resolved with more/fewer DM particles. △ Less

Submitted 4 September, 2023; v1 submitted 9 June, 2023; originally announced June 2023.

Comments: 17 pages, 15 figures. MNRAS published. A basic python script for calculating spurious disc heating rates is available at https://github.com/AaronDLudlow/hot-disc.git

The impact of stochastic modeling on the predictive power of galaxy formation simulations

Authors: Josh Borrow, Matthieu Schaller, Yannick M. Bahe, Joop Schaye, Aaron D. Ludlow, Sylvia Ploeckinger, Folkert S. J. Nobels, Edoardo Altamura

Abstract: All modern galaxy formation models employ stochastic elements in their sub-grid prescriptions to discretise continuous equations across the time domain. In this paper, we investigate how the stochastic nature of these models, notably star formation, black hole accretion, and their associated feedback, that act on small ($<$ kpc) scales, can back-react on macroscopic galaxy properties (e.g. stellar… ▽ More All modern galaxy formation models employ stochastic elements in their sub-grid prescriptions to discretise continuous equations across the time domain. In this paper, we investigate how the stochastic nature of these models, notably star formation, black hole accretion, and their associated feedback, that act on small ($<$ kpc) scales, can back-react on macroscopic galaxy properties (e.g. stellar mass and size) across long ($>$ Gyr) timescales. We find that the scatter in scaling relations predicted by the EAGLE model implemented in the SWIFT code can be significantly impacted by random variability between re-simulations of the same object, even when galaxies are resolved by tens of thousands of particles. We then illustrate how re-simulations of the same object can be used to better understand the underlying model, by showing how correlations between galaxy stellar mass and black hole mass disappear at the highest black hole masses ($M_{\rm BH} > 10^8$ M$_\odot$), indicating that the feedback cycle may be interrupted by external processes. We find that although properties that are collected cumulatively over many objects are relatively robust against random variability (e.g. the median of a scaling relation), the properties of individual galaxies (such as galaxy stellar mass) can vary by up to 25\%, even far into the well-resolved regime, driven by bursty physics (black hole feedback) and mergers between galaxies. We suggest that studies of individual objects within cosmological simulations be treated with caution, and that any studies aiming to closely investigate such objects must account for random variability within their results. △ Less

Submitted 22 September, 2023; v1 submitted 15 November, 2022; originally announced November 2022.

Comments: Accepted for publication in MNRAS

The physical drivers of gas turbulence in simulated disc galaxies

Authors: Esteban Jiménez, Claudia del P. Lagos, Aaron D. Ludlow, Emily Wisnioski

Abstract: We use the EAGLE cosmological simulations to study the evolution of the vertical velocity dispersion of cold gas, $σ_{z}$, in central disc galaxies and its connection to stellar feedback, gravitational instabilities, cosmological gas accretion and galaxy mergers. To isolate the impact of feedback, we analyse runs that turn off stellar and (or) AGN feedback in addition to a run that includes both.… ▽ More We use the EAGLE cosmological simulations to study the evolution of the vertical velocity dispersion of cold gas, $σ_{z}$, in central disc galaxies and its connection to stellar feedback, gravitational instabilities, cosmological gas accretion and galaxy mergers. To isolate the impact of feedback, we analyse runs that turn off stellar and (or) AGN feedback in addition to a run that includes both. The evolution of $σ_z$ and its dependence on stellar mass and star formation rate in EAGLE are in good agreement with observations. Galaxies hosted by haloes of similar virial mass, $\rm M_{200}$, have similar $σ_z$ values even in runs where feedback is absent. The prevalence of local instabilities in discs is uncorrelated with $σ_z$ at low redshift and becomes only weakly correlated at high redshifts and in galaxies hosted by massive haloes. $σ_z$ correlates most strongly with the specific gas accretion rate onto the disc as well as with the degree of misalignment between the inflowing gas and the disc's rotation axis. These correlations are significant across all redshifts and halo masses, with misaligned accretion being the primary driver of high gas turbulence at redshifts $z \lesssim 1$ and for halo masses $\rm M_{200} \lesssim 10^{11.5} M_{\odot}$. Galaxy mergers increase $σ_z$, but because they are rare in our sample, they play only a minor role in its evolution. Our results suggest that the turbulence of cold gas in EAGLE discs results from a complex interplay of different physical processes whose relative importance depends on halo mass and redshift. △ Less

Submitted 13 July, 2023; v1 submitted 18 October, 2022; originally announced October 2022.

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

The impact of spurious collisional heating on the morphological evolution of simulated galactic discs

Authors: Matthew J. Wilkinson, Aaron D. Ludlow, Claudia del P. Lagos, S. Michael Fall, Joop Schaye, Danail Obreschkow

Abstract: We use a suite of idealised N-body simulations to study the impact of spurious heating of star particles by dark matter particles on the kinematics and morphology of simulated galactic discs. We find that spurious collisional heating leads to a systematic increase of the azimuthal velocity dispersion ($σ_φ$) of stellar particles and a corresponding decrease in their mean azimuthal velocities (… ▽ More We use a suite of idealised N-body simulations to study the impact of spurious heating of star particles by dark matter particles on the kinematics and morphology of simulated galactic discs. We find that spurious collisional heating leads to a systematic increase of the azimuthal velocity dispersion ($σ_φ$) of stellar particles and a corresponding decrease in their mean azimuthal velocities ($\overline{v}_φ$). The rate of heating is dictated primarily by the number of dark matter halo particles (or equivalently, by the dark matter particle mass at fixed halo mass) and by radial gradients in the local dark matter density along the disc; it is largely insensitive to the stellar particle mass. Galaxies within haloes resolved with fewer than $\approx 10^6$ dark matter particles are particularly susceptible to spurious morphological evolution, irrespective of the total halo mass (with even more particles required to prevent heating of the galactic centre). Collisional heating transforms galactic discs from flattened structures into rounder spheroidal systems, causing them to lose rotational support in the process. It also affects the locations of galaxies in standard scaling relations that link their various properties: at fixed stellar mass, it increases the sizes of galaxies, and reduces their mean stellar rotation velocities and specific angular momenta. Our results urge caution when extrapolating simulated galaxy scaling relations to low masses where spurious collisional effects can bias their normalisation, slope and scatter. △ Less

Submitted 30 January, 2023; v1 submitted 16 August, 2022; originally announced August 2022.

Comments: Minor changes to match MNRAS published version. 20 pages, 13 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 519, Issue 4, March 2023, Pages 5942-5961

The cosmology dependence of the concentration-mass-redshift relation

Authors: Daniel López-Cano, Raúl E. Angulo, Aaron D. Ludlow, M. Zennaro, S. Contreras, Jonás Chaves-Montero, G. Aricò

Abstract: The concentrations of dark matter haloes provide crucial information about their internal structure and how it depends on mass and redshift -- the so-called concentration-mass-redshift relation, denoted $c(M,z)$. We present here an extensive study of the cosmology-dependence of $c(M,z)$ that is based on a suite of 72 gravity-only, full N-body simulations in which the following cosmological paramet… ▽ More The concentrations of dark matter haloes provide crucial information about their internal structure and how it depends on mass and redshift -- the so-called concentration-mass-redshift relation, denoted $c(M,z)$. We present here an extensive study of the cosmology-dependence of $c(M,z)$ that is based on a suite of 72 gravity-only, full N-body simulations in which the following cosmological parameters were varied: $σ_{8}$, $Ω_{\mathrm{M}}$, $Ω_{\mathrm{b}}$, $n_{\mathrm{s}}$, $h$, $M_ν$, $w_{0}$ and $w_{\mathrm{a}}$. We characterize the impact of these parameters on concentrations for different halo masses and redshifts. In agreement with previous works, and for all cosmologies studied, we find that there exists a tight correlation between the characteristic densities of dark matter haloes within their scale radii, $r_{-2}$, and the critical density of the Universe at a suitably defined formation time. This finding, when combined with excursion set modelling of halo formation histories, allows us to accurately predict the concentrations of dark matter haloes as a function of mass, redshift, and cosmology. We use our simulations to test the reliability of a number of published models for predicting halo concentration and highlight when they succeed or fail to reproduce the cosmological $c(M,z)$ relation. △ Less

Submitted 27 July, 2022; originally announced July 2022.

Comments: 11 pages, 9 figures

Drivers of asymmetry in synthetic H I emission-line profiles of galaxies in the EAGLE simulation

Authors: Aditya Manuwal, Aaron D. Ludlow, Adam R. H. Stevens, Ruby J. Wright, Aaron S. G. Robotham

Abstract: We study the shapes of spatially integrated H I emission-line profiles of galaxies in the EAGLE simulation using three separate measures of the profile's asymmetry. We show that the subset of EAGLE galaxies whose gas fractions and stellar masses are consistent with those in the xGASS survey also have similar H I line asymmetries. Central galaxies with symmetric H I line profiles typically correspo… ▽ More We study the shapes of spatially integrated H I emission-line profiles of galaxies in the EAGLE simulation using three separate measures of the profile's asymmetry. We show that the subset of EAGLE galaxies whose gas fractions and stellar masses are consistent with those in the xGASS survey also have similar H I line asymmetries. Central galaxies with symmetric H I line profiles typically correspond to rotationally supported H I and stellar disks, but those with asymmetric line profiles may or may not correspond to dispersion-dominated systems. Galaxies with symmetric H I emission lines are, on average, more gas rich than those with asymmetric lines, and also exhibit systematic differences in their specific star formation rates, suggesting that turbulence generated by stellar or AGN feedback may be one factor contributing to H I line asymmetry. The line asymmetry also correlates strongly with the dynamical state of a galaxy's host dark matter halo: older, more relaxed haloes host more-symmetric galaxies than those hosted by unrelaxed ones. At fixed halo mass, asymmetric centrals tend to be surrounded by a larger number of massive subhaloes than their symmetric counterparts, and also experience higher rates of gas accretion and outflow. At fixed stellar mass, central galaxies have, on average, more symmetric H I emission lines than satellites; for the latter, ram pressure and tidal stripping are significant sources of asymmetry. △ Less

Submitted 29 December, 2021; v1 submitted 23 September, 2021; originally announced September 2021.

Comments: Accepted for publication in MNRAS

Spurious heating of stellar motions in simulated galactic disks by dark matter halo particles

Authors: Aaron D. Ludlow, S. Michael Fall, Joop Schaye, Danail Obreschkow

Abstract: We use idealized N-body simulations of equilibrium stellar disks embedded within course-grained dark matter haloes to study the effects of spurious collisional heating on disk structure and kinematics. Collisional heating artificially increases the vertical and radial velocity dispersions of disk stars, as well as the thickness and size of disks; the effects are felt at all galacto-centric radii.… ▽ More We use idealized N-body simulations of equilibrium stellar disks embedded within course-grained dark matter haloes to study the effects of spurious collisional heating on disk structure and kinematics. Collisional heating artificially increases the vertical and radial velocity dispersions of disk stars, as well as the thickness and size of disks; the effects are felt at all galacto-centric radii. The integrated effects of collisional heating are determined by the mass of dark matter halo particles (or equivalently, by the number of particles at fixed halo mass), their local density and characteristic velocity dispersion, but are largely insensitive to the stellar particle mass. The effects can therefore be reduced by increasing the mass resolution of dark matter in cosmological simulations, with limited benefits from increasing the baryonic (or stellar) mass resolution. We provide a simple empirical model that accurately captures the effects of spurious collisional heating on the structure and kinematics of simulated disks, and use it to assess the importance of disk heating for simulations of galaxy formation. We find that the majority of state-of-the-art zoom simulations, and a few of the highest-resolution, smallest-volume cosmological runs, are in principle able to resolve thin stellar disks in Milky Way-mass haloes, but most large-volume cosmological simulations cannot. For example, dark matter haloes resolved with fewer than $\approx 10^6$ particles will collisionally heat stars near the stellar half-mass radius such that their vertical velocity dispersion increases by $\gtrsim 10$ per cent of the halo's virial velocity in approximately one Hubble time. △ Less

Submitted 23 September, 2021; v1 submitted 7 May, 2021; originally announced May 2021.

Comments: Minor changes to match MNRAS accepted version. Main text: 19 pages, 10 Figure (Appendix: 5 pages, 6 figures)

The orbital evolution of UFDs and GCs in an evolving Galactic potential

Authors: Benjamin M. Armstrong, Kenji Bekki, Aaron D. Ludlow

Abstract: We use the second Gaia data release to investigate the kinematics of 17 ultra-faint dwarf galaxies (UFDs) and 154 globular clusters (GCs) in the Milky Way, focusing on the differences between static and evolving models of the Galactic potential. An evolving potential modifies a satellite's orbit relative to its static equivalent, though the difference is small compared to existing uncertainties on… ▽ More We use the second Gaia data release to investigate the kinematics of 17 ultra-faint dwarf galaxies (UFDs) and 154 globular clusters (GCs) in the Milky Way, focusing on the differences between static and evolving models of the Galactic potential. An evolving potential modifies a satellite's orbit relative to its static equivalent, though the difference is small compared to existing uncertainties on orbital parameters. We find that the UFD Boötes II is likely on its first passage around the Milky Way. Depending on the assumed mass of the Milky Way, the UFDs Triangulum II, Hydrus I, Coma Berenices, Draco II, and Ursa Major II, as well as the GC Pyxis, may also be on first infall so may be useful for constraining the mass of the Galaxy. We identify a clear kinematic distinction between metal-rich (${\rm [Fe/H]}>-1.1$) and metal-poor GCs (${\rm [Fe/H]}\leq-1.1$). Although most metal-rich clusters occupy predominately prograde orbits, with low eccentricities ($e\approx 0.35$) and similar specific angular momenta and orbital planes as the Galactic disc, 7 show potentially retrograde orbits, the origin of which is unclear. Metal-poor clusters have more diverse orbits, higher eccentricities ($e\approx 0.65$), and half have orbital planes offset from the disc by 60 to 120 degrees. The UFDs have similar $θ$ and $φ$ to the metal-poor GCs, suggesting a similar origin. We provide a catalogue of orbital parameters for UFDs and GCs for two different Galaxy masses and their observational uncertainties. △ Less

Submitted 6 December, 2020; v1 submitted 25 November, 2020; originally announced November 2020.

Comments: 23 pages, 12 figures, accepted for publication in MNRAS

Pre-processing, group accretion and the orbital trajectories of associated subhaloes

Authors: Lucie Bakels, Aaron D. Ludlow, Chris Power

Abstract: We use a high-resolution cosmological dark matter-only simulation to study the orbital trajectories of haloes and subhaloes in the environs of isolated hosts. We carefully tally all apsis points and use them to distinguish haloes that are infalling for the first time from those that occupy more evolved orbits. We find that roughly 21 per cent of subhaloes within a host's virial radius are currentl… ▽ More We use a high-resolution cosmological dark matter-only simulation to study the orbital trajectories of haloes and subhaloes in the environs of isolated hosts. We carefully tally all apsis points and use them to distinguish haloes that are infalling for the first time from those that occupy more evolved orbits. We find that roughly 21 per cent of subhaloes within a host's virial radius are currently on first infall, and have not yet reached their first orbital pericentre; roughly 44 per cent are still approaching their first apocentre after infall. For the range of host masses studied, roughly half of all accreted systems were pre-processed prior to infall, and about 20 per cent were accreted in groups. We confirm that the entire population of accreted subhaloes -- often referred to as "associated" subhaloes -- extend far beyond the virial radii of their hosts, with roughly half currently residing at distances that exceed $\approx 1.2\times r_{200}$. Many of these backsplash haloes have gained orbital energy since infall, and occupy extreme orbits that carry them well past their initial turnaround radii. Such extreme orbits are created during the initial accretion and dissolution of loosely bound groups, but also through penetrating encounters between subhaloes on subsequent orbits. The same processes may also give rise to unexpectedly abrupt losses of orbital energy. These effects combine, giving rise to a large variation in the ratio of sequent apocentres for accreted systems. We find that, within 2 virial radii from host centres, the concentrations of first-infall halos are remarkably similar those of isolated field halos, whereas backsplash haloes, as well as systems that were pre-processed, are considerably more concentrated. △ Less

Submitted 7 January, 2021; v1 submitted 12 August, 2020; originally announced August 2020.

Comments: 17 pages, 12 figures, submitted to MNRAS

Observational constraints on the slope of the radial acceleration relation at low accelerations

Authors: Kyle A. Oman, Margot M. Brouwer, Aaron D. Ludlow, Julio F. Navarro

Abstract: The radial acceleration relation (RAR) locally relates the `observed' acceleration inferred from the dynamics of a system to the acceleration implied by its baryonic matter distribution. The relation as traced by galaxy rotation curves is one-to-one with remarkably little scatter, implying that the dynamics of a system can be predicted simply by measuring its density profile as traced by e.g. stel… ▽ More The radial acceleration relation (RAR) locally relates the `observed' acceleration inferred from the dynamics of a system to the acceleration implied by its baryonic matter distribution. The relation as traced by galaxy rotation curves is one-to-one with remarkably little scatter, implying that the dynamics of a system can be predicted simply by measuring its density profile as traced by e.g. stellar light or gas emission lines. Extending the relation to accelerations below those usually probed by practically observable kinematic tracers is challenging, especially once accounting for faintly emitting baryons, such as the putative warm-hot intergalactic medium, becomes important. We show that in the low-acceleration regime, the (inverted) RAR predicts an unphysical, declining enclosed baryonic mass profile for systems with `observed' acceleration profiles steeper than $g_{\rm obs}\propto r^{-1}$ (corresponding to density profiles steeper than isothermal - $ρ(r)\propto r^{-2}$). If the RAR is tantamount to a natural law, such acceleration profiles cannot exist. We apply this argument to test the compatibility of an extrapolation of the rotation curve-derived RAR to low accelerations with data from galaxy-galaxy weak lensing, dwarf spheroidal galaxy stellar kinematic, and outer Milky~Way dynamical measurements, fully independent of the uncertainties inherent in direct measurements of the baryonic matter distribution. In all cases we find that the data weakly favour a break to a steeper low-acceleration slope. Improvements in measurements and modelling of the outer Milky~Way, and weak lensing, seem like the most promising path toward stronger constraints on the low-acceleration behaviour of the RAR. △ Less

Submitted 16 December, 2020; v1 submitted 11 June, 2020; originally announced June 2020.

Comments: 11 pages, 8 figures (+appendices), v2 corrects a factual error in Sec. 2

A comparison of $\text{H}_2$ formation models at high redshift

Authors: Alexander Schäbe, Emilio Romano-Díaz, Cristiano Porciani, Aaron D. Ludlow, Matteo Tomassetti

Abstract: Modelling the molecular gas that is routinely detected through CO observations of high-redshift galaxies constitutes a major challenge for ab initio simulations of galaxy formation. We carry out a suite of cosmological hydrodynamic simulations to compare three approximate methods that have been used in the literature to track the formation and evolution of the simplest and most abundant molecule,… ▽ More Modelling the molecular gas that is routinely detected through CO observations of high-redshift galaxies constitutes a major challenge for ab initio simulations of galaxy formation. We carry out a suite of cosmological hydrodynamic simulations to compare three approximate methods that have been used in the literature to track the formation and evolution of the simplest and most abundant molecule, H$_2$. Namely, we consider: i) a semi-empirical procedure that associates H$_2$ to dark-matter haloes based on a series of scaling relations inferred from observations, ii) a model that assumes chemical equilibrium between the H$_2$ formation and destruction rates, and iii) a model that fully solves the out-of-equilibrium rate equations and accounts for the unresolved structure of molecular clouds. We study the impact of finite spatial resolution and show that robust H$_2$ masses at redshift $z\approx 4$ can only be obtained for galaxies that are sufficiently metal enriched in which H$_2$ formation is fast. This corresponds to H$_2$ reservoirs with masses $M_{\mathrm{H_2}}\gtrsim 6\times 10^9 \mathrm{M}_\odot$. In this range, equilibrium and non-equilibrium models predict similar molecular masses (but different galaxy morphologies) while the semi-empirical method produces less H$_2$. The star formation rates as well as the stellar and H$_2$ masses of the simulated galaxies are in line with those observed in actual galaxies at similar redshifts that are not massive starbursts. The H$_2$ mass functions extracted from the simulations at $z\approx 4$ agree well with recent observations that only sample the high-mass end. However, our results indicate that most molecular material at high $z$ lies yet undetected in reservoirs with $10^9<M_{\mathrm H_2}<10^{10} \mathrm{M}_\odot$. △ Less

Submitted 3 September, 2020; v1 submitted 9 March, 2020; originally announced March 2020.

Comments: 19 pages, 12 figures, 1 table, accepted for publication in MNRAS, comments are welcome

Local Group star formation in warm and self-interacting dark matter cosmologies

Authors: Mark R. Lovell, Wojciech Hellwing, Aaron Ludlow, Jesús Zavala, Andrew Robertson, Azadeh Fattahi, Carlos S. Frenk, Jennifer Hardwick

Abstract: The nature of the dark matter can affect the collapse time of dark matter haloes, and can therefore be imprinted in observables such as the stellar population ages and star formation histories of dwarf galaxies. In this paper we use high resolution hydrodynamical simulations of Local Group-analogue (LG) volumes in cold dark matter (CDM), sterile neutrino warm dark matter (WDM) and self-interacting… ▽ More The nature of the dark matter can affect the collapse time of dark matter haloes, and can therefore be imprinted in observables such as the stellar population ages and star formation histories of dwarf galaxies. In this paper we use high resolution hydrodynamical simulations of Local Group-analogue (LG) volumes in cold dark matter (CDM), sterile neutrino warm dark matter (WDM) and self-interacting dark matter (SIDM) models with the EAGLE galaxy formation code to study how galaxy formation times change with dark matter model. We are able to identify the same haloes in different simulations, since they share the same initial density field phases. We find that the stellar mass of galaxies depends systematically on resolution, and can differ by as much as a factor of two in haloes of a given dark matter mass. The evolution of the stellar populations in SIDM is largely identical to that of CDM, but in WDM early star formation is instead suppressed. The time at which LG haloes can begin to form stars through atomic cooling is delayed by $\sim$200~Myr in WDM models compared to CDM. It will be necessary to measure stellar ages of old populations to a precision of better than 100~Myr, and to address degeneracies with the redshift of reionization -- and potentially other baryonic processes -- in order to use these observables to distinguish between dark matter models. △ Less

Submitted 25 September, 2020; v1 submitted 25 February, 2020; originally announced February 2020.

Comments: 17 pages, 13 figures. Version accepted by MNRAS. Contact: lovell@hi.is

Characterising the Structure of Halo Merger Trees Using a Single Parameter: The Tree Entropy

Authors: Danail Obreschkow, Pascal J. Elahi, Claudia del P. Lagos, Rhys J. J. Poulton, Aaron D. Ludlow

Abstract: Linking the properties of galaxies to the assembly history of their dark matter haloes is a central aim of galaxy evolution theory. This paper introduces a dimensionless parameter $s\in[0,1]$, the "tree entropy", to parametrise the geometry of a halo's entire mass assembly hierarchy, building on a generalisation of Shannon's information entropy. By construction, the minimum entropy ($s=0$) corresp… ▽ More Linking the properties of galaxies to the assembly history of their dark matter haloes is a central aim of galaxy evolution theory. This paper introduces a dimensionless parameter $s\in[0,1]$, the "tree entropy", to parametrise the geometry of a halo's entire mass assembly hierarchy, building on a generalisation of Shannon's information entropy. By construction, the minimum entropy ($s=0$) corresponds to smoothly assembled haloes without any mergers. In contrast, the highest entropy ($s=1$) represents haloes grown purely by equal-mass binary mergers. Using simulated merger trees extracted from the cosmological $N$-body simulation SURFS, we compute the natural distribution of $s$, a skewed bell curve peaking near $s=0.4$. This distribution exhibits weak dependences on halo mass $M$ and redshift $z$, which can be reduced to a single dependence on the relative peak height $δ_{\rm c}/σ(M,z)$ in the matter perturbation field. By exploring the correlations between $s$ and global galaxy properties generated by the SHARK semi-analytic model, we find that $s$ contains a significant amount of information on the morphology of galaxies $-$ in fact more information than the spin, concentration and assembly time of the halo. Therefore, the tree entropy provides an information-rich link between galaxies and their dark matter haloes. △ Less

Submitted 11 February, 2020; v1 submitted 27 November, 2019; originally announced November 2019.

Comments: 20 pages, 18 figures

To beta or not to beta: can higher-order Jeans analysis break the mass-anisotropy degeneracy in simulated dwarfs?

Authors: Anna Genina, Justin I. Read, Carlos S. Frenk, Shaun Cole, Alejandro Benitez-Llambay, Aaron D. Ludlow, Julio F. Navarro, Kyle A. Oman, Andrew Robertson

Abstract: We test a non-parametric higher-order Jeans analysis method, GravSphere, on 32 simulated dwarf galaxies comparable to classical Local Group dwarfs like Fornax. The galaxies are selected from the APOSTLE suite of cosmological hydrodynamics simulations with Cold Dark Matter (CDM) and Self-Interacting Dark Matter (SIDM) models, allowing us to investigate cusps and cores in density distributions. We f… ▽ More We test a non-parametric higher-order Jeans analysis method, GravSphere, on 32 simulated dwarf galaxies comparable to classical Local Group dwarfs like Fornax. The galaxies are selected from the APOSTLE suite of cosmological hydrodynamics simulations with Cold Dark Matter (CDM) and Self-Interacting Dark Matter (SIDM) models, allowing us to investigate cusps and cores in density distributions. We find that, for CDM dwarfs, the recovered enclosed mass profiles have a bias of no more than 10 per cent, with a 50 per cent scatter in the inner regions and a 20 per cent scatter near the half-light radius, consistent with standard mass estimators. The density profiles are also recovered with a bias of no more than 10 per cent and a scatter of 30 per cent in the inner regions. For SIDM dwarfs, the mass and density profiles are recovered within our 95 per cent confidence intervals, but are biased towards cuspy dark matter distributions. This is mainly due to a lack of sufficient constraints from the data. We explore the sources of scatter in the accuracy of the recovered profiles and suggest a $χ^2$ statistic to separate successful models from biased ones. Finally, we show that the uncertainties on the mass profiles obtained with GravSphere are smaller than those for comparable Jeans methods, and that they can be further improved if stronger priors, motivated by cosmological simulations, are placed on the velocity anisotropy. We conclude that GravSphere is a promising Jeans-based approach for modelling dark matter distributions in dwarf galaxies. △ Less

Submitted 5 August, 2020; v1 submitted 20 November, 2019; originally announced November 2019.

Comments: 22 pages, 13 figures. Accepted in MNRAS

Baryonic clues to the puzzling diversity of dwarf galaxy rotation curves

Authors: Isabel M. E. Santos-Santos, Julio F. Navarro, Andrew Robertson, Alejandro Benítez-Llambay, Kyle A. Oman, Mark R. Lovell, Carlos S. Frenk, Aaron D. Ludlow, Azadeh Fattahi, Adam Ritz

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

Submitted 28 April, 2020; v1 submitted 20 November, 2019; originally announced November 2019.

Comments: 20 pages, 16 figures, 1 table. Accepted for publication in MNRAS. Updated acknowledgements

Numerical convergence of hydrodynamical simulations of galaxy formation: the abundance and internal structure of galaxies and their cold dark matter haloes

Authors: Aaron D. Ludlow, Joop Schaye, Matthieu Schaller, Richard Bower

Abstract: We address the issue of numerical convergence in cosmological smoothed particle hydrodynamics simulations using a suite of runs drawn from the EAGLE project. Our simulations adopt subgrid models that produce realistic galaxy populations at a fiducial mass and force resolution, but systematically vary the latter in order to study their impact on galaxy properties. We provide several analytic criter… ▽ More We address the issue of numerical convergence in cosmological smoothed particle hydrodynamics simulations using a suite of runs drawn from the EAGLE project. Our simulations adopt subgrid models that produce realistic galaxy populations at a fiducial mass and force resolution, but systematically vary the latter in order to study their impact on galaxy properties. We provide several analytic criteria that help guide the selection of gravitational softening for hydrodynamical simulations, and present results from runs that both adhere to and deviate from them. Unlike dark matter-only simulations, hydrodynamical simulations exhibit a strong sensitivity to gravitational softening, and care must be taken when selecting numerical parameters. Our results--which focus mainly on star formation histories, galaxy stellar mass functions and sizes--illuminate three main considerations. First, softening imposes a minimum resolved escape speed, $v_ε$, due to the binding energy between gas particles. Runs that adopt such small softening lengths that $v_ε\gt 10\,{\rm km s^{-1}}$ (the sound speed in ionised $\sim 10^4\,{\rm K}$ gas) suffer from reduced effects of photo-heating. Second, feedback from stars or active galactic nuclei may suffer from numerical over-cooling if the gravitational softening length is chosen below a critical value, $ε_{\rm eFB}$. Third, we note that small softening lengths exacerbate the segregation of stars and dark matter particles in halo centres, often leading to the counter-intuitive result that galaxy sizes {\em increase} as softening is reduced. The structure of dark matter haloes in hydrodynamical runs respond to softening in a way that reflects the sensitivity of their galaxy populations to numerical parameters. △ Less

Submitted 7 May, 2020; v1 submitted 14 August, 2019; originally announced August 2019.

Comments: 28 pages, 15 figures; Accepted 2020 January 29 (Figures with reduced resolution; hi-res figures available on request)

Energy equipartition between stellar and dark matter particles in cosmological simulations results in spurious growth of galaxy sizes

Authors: Aaron D. Ludlow, Joop Schaye, Matthieu Schaller, Jack Richings

Abstract: The impact of 2-body scattering on the innermost density profiles of dark matter haloes is well established. We use a suite of cosmological simulations and idealised numerical experiments to show that 2-body scattering is exacerbated in situations where there are two species of unequal mass. This is a consequence of mass segregation and reflects a flow of kinetic energy from the more to less massi… ▽ More The impact of 2-body scattering on the innermost density profiles of dark matter haloes is well established. We use a suite of cosmological simulations and idealised numerical experiments to show that 2-body scattering is exacerbated in situations where there are two species of unequal mass. This is a consequence of mass segregation and reflects a flow of kinetic energy from the more to less massive particles. This has important implications for the interpretation of galaxy sizes in cosmological hydrodynamic simulations, which nearly always model stars with less massive particles than are used for the dark matter. We compare idealised models as well as simulations from the EAGLE project that differ only in the mass resolution of the dark matter component, but keep sub-grid physics, baryonic mass resolution and gravitational force softening fixed. If the dark matter particle mass exceeds the mass of stellar particles, then galaxy sizes--quantified by their projected half-mass radii, ${\rm R_{50}}$--increase systematically with time until ${\rm R_{50}}$ exceeds a small fraction of the redshift-dependent mean inter-particle separation, $l$ (${\rm R_{50}}\geq 0.05\times l$). Our conclusions should also apply to simulations that adopt different hydrodynamic solvers, subgrid physics or adaptive softening, but in that case may need quantitative revision. Any simulation employing a stellar-to-dark matter particle mass ratio greater than unity will escalate spurious energy transfer from dark matter to baryons on small scales. △ Less

Submitted 10 July, 2019; v1 submitted 24 March, 2019; originally announced March 2019.

Comments: 6 pages, 3 figures; replaced by MNRAS accepted version

The connection between halo concentrations and assembly histories: a probe of gravity?

Authors: Piotr Oleśkiewicz, Carlton M. Baugh, Aaron Ludlow

Abstract: We use two high resolution N-body simulations, one assuming general relativity and the other the Hu-Sawicki form of $f\left(R\right)$ gravity with $\vert \bar{f}_{R} \vert = 10^{-6}$, to investigate the concentration--formation time relation of dark matter haloes. We stack haloes in logarithmically spaced mass bins to fit median density profiles and extract median formation times. At fixed mass, h… ▽ More We use two high resolution N-body simulations, one assuming general relativity and the other the Hu-Sawicki form of $f\left(R\right)$ gravity with $\vert \bar{f}_{R} \vert = 10^{-6}$, to investigate the concentration--formation time relation of dark matter haloes. We stack haloes in logarithmically spaced mass bins to fit median density profiles and extract median formation times. At fixed mass, haloes in modified gravity are more concentrated than those in GR, especially at low masses and at low redshift, and do not follow the concentration--formation time relation seen in GR. We assess the sensitivity of the relation to how concentration and formation time are defined, as well as to the segregation of the halo population by the amount of gravitational screening. We find a clear difference between halo concentrations and assembly histories displayed in modified gravity and those in GR. Existing models for the mass--concentration--redshift relation that have gained success in cold and warm dark matter models require revision in $f\left(R\right)$ gravity. △ Less

Submitted 18 March, 2019; originally announced March 2019.

Comments: 12 pages, 10 figures

Galaxy And Mass Assembly (GAMA): Environmental Quenching of Centrals and Satellites in Groups

Authors: L. J. M. Davies, A. S. G. Robotham, C. del P. Lagos, S. P. Driver, A. R. H. Stevens, Y. M. Bahé, M. Alpaslan, M. N. Bremer, M. J. I. Brown, S. Brough, J. Bland-Hawthorn, L. Cortese, P. Elahi, M. W. Grootes, B. W. Holwerda, A. D. Ludlow, S. McGee, M. Owers, S. Phillipps

Abstract: Recently a number of studies have found a similarity between the passive fraction of central and satellite galaxies when controlled for both stellar and halo mass. These results suggest that the quenching processes that affect galaxies are largely agnostic to central/satellite status, which contradicts the traditional picture of increased satellite quenching via environmental processes such as str… ▽ More Recently a number of studies have found a similarity between the passive fraction of central and satellite galaxies when controlled for both stellar and halo mass. These results suggest that the quenching processes that affect galaxies are largely agnostic to central/satellite status, which contradicts the traditional picture of increased satellite quenching via environmental processes such as stripping, strangulation and starvation. Here we explore this further using the Galaxy And Mass Assembly (GAMA) survey which extends to ~2dex lower in stellar mass than SDSS, is more complete for closely-separated galaxies (>95% compared to >70%), and identifies lower-halo-mass groups outside of the very local Universe (M$_{\mathrm{halo}}\sim10^{12}$M$_{\odot}$ at $0.1<z<0.2$). As far as possible we aim to replicate the selections, completeness corrections and central/satellite division of one of the previous studies but find clear differences between passive fractions of centrals and satellites. We also find that our passive fractions increase with both halo-to-satellite mass ratio and central-to-second rank mass ratio. This suggests that quenching is more efficient in satellites that are low-mass for their halo ($i.e$ at high halo-to-satellite mass ratio in comparison to low halo-to-satellite mass ratio) and are more likely to be passive in older groups - forming a consistent picture of environmental quenching of satellites. We then discuss potential explanations for the previously observed similarity, such as dependence on the group-finding method. △ Less

Submitted 6 January, 2019; originally announced January 2019.

Comments: 15 pages, 8 figures, MNRAS accepted

Halo concentrations from extended Press-Schechter merger histories

Authors: Andrew J. Benson, Aaron Ludlow, Shaun Cole

Abstract: We apply the model relating halo concentration to formation history proposed by Ludlow et al. to merger trees generated using an algorithm based on excursion set theory. We find that while the model correctly predicts the median relation between halo concentration and mass, it underpredicts the scatter in concentration at fixed mass. Since the same model applied to N-body merger trees predicts the… ▽ More We apply the model relating halo concentration to formation history proposed by Ludlow et al. to merger trees generated using an algorithm based on excursion set theory. We find that while the model correctly predicts the median relation between halo concentration and mass, it underpredicts the scatter in concentration at fixed mass. Since the same model applied to N-body merger trees predicts the correct scatter, we postulate that the missing scatter is due to the lack of any environmental dependence in merger trees derived from excursion set theory. We show that a simple modification to the merger tree construction algorithm, which makes merger rates dependent on environment, can increase the scatter by the required amount, and simultaneously provide a qualitatively correct correlation between environment and formation epoch in the excursion set merger trees. △ Less

Submitted 14 December, 2018; originally announced December 2018.

Comments: 10 pages, 7 figures, submitted to MNRAS

Numerical convergence of simulations of galaxy formation: the abundance and internal structure of cold dark matter haloes

Authors: Aaron D. Ludlow, Joop Schaye, Richard Bower

Abstract: We study the impact of numerical parameters on the properties of cold dark matter haloes formed in collisionless cosmological simulations. We quantify convergence in the median spherically-averaged circular velocity profiles for haloes of widely varying particle number, as well as in the statistics of their structural scaling relations and mass functions. In agreement with prior work focused on si… ▽ More We study the impact of numerical parameters on the properties of cold dark matter haloes formed in collisionless cosmological simulations. We quantify convergence in the median spherically-averaged circular velocity profiles for haloes of widely varying particle number, as well as in the statistics of their structural scaling relations and mass functions. In agreement with prior work focused on single haloes, our results suggest that cosmological simulations yield robust halo properties for a wide range of gravitational softening parameters, $ε$, provided: 1) $ε$ is not larger than a "convergence radius", $r_{\rm conv}$, which is dictated by 2-body relaxation and determined by particle number, and 2) a sufficient number of timesteps are taken to accurately resolve particle orbits with short dynamical times. Provided these conditions are met, median circular velocity profiles converge to within $\approx 10$ per cent for radii beyond which the local 2-body relaxation timescale exceeds the Hubble time by a factor $κ\equiv t_{\rm relax}/t_{\rm H}\gt 0.177$, with better convergence attained for higher $κ$. We provide analytic estimates of $r_{\rm conv}$ that build on previous attempts in two ways: first, by highlighting its explicit (but weak) softening-dependence and, second, by providing a simpler criterion in which $r_{\rm conv}$ is determined entirely by the mean inter-particle spacing, $l$; for example, better than $10$ per cent convergence in circular velocity for $r\gt 0.05\,l$. We show how these analytic criteria can be used to assess convergence in structural scaling relations for dark matter haloes as a function of their mass or maximum circular speed. △ Less

Submitted 10 July, 2019; v1 submitted 13 December, 2018; originally announced December 2018.

Comments: 24 pages, 15 figures; MNRAS-accepted version

The signal of decaying dark matter with hydrodynamical simulations

Authors: Mark R. Lovell, David Barnes, Yannick Bahé, Joop Schaye, Matthieu Schaller, Tom Theuns, Sownak Bose, Robert A. Crain, Claudio dalla Vecchia, Carlos S. Frenk, Wojciech Hellwing, Scott T. Kay, Aaron D. Ludlow, Richard G. Bower

Abstract: Dark matter particles may decay, emitting photons. Drawing on the EAGLE family of hydrodynamic simulations of galaxy formation -- including the APOSTLE and C-EAGLE simulations -- we assess the systematic uncertainties and scatter on the decay flux from different galaxy classes, from Milky Way satellites to galaxy clusters, and compare our results to studies of the 3.55~keV line. We demonstrate tha… ▽ More Dark matter particles may decay, emitting photons. Drawing on the EAGLE family of hydrodynamic simulations of galaxy formation -- including the APOSTLE and C-EAGLE simulations -- we assess the systematic uncertainties and scatter on the decay flux from different galaxy classes, from Milky Way satellites to galaxy clusters, and compare our results to studies of the 3.55~keV line. We demonstrate that previous detections and non-detections of this line are consistent with a dark matter interpretation. For example, in our simulations the width of the the dark matter decay line for Perseus-analogue galaxy clusters lies in the range 1300-1700~\kms. Therefore, the non-detection of the 3.55~keV line in the centre of the Perseus cluster by the {\it Hitomi} collaboration is consistent with detections by other instruments. We also consider trends with stellar and halo mass and evaluate the scatter in the expected fluxes arising from the anisotropic halo mass distribution and from object-to-object variations. We provide specific predictions for observations with {\it XMM-Newton} and with the planned X-ray telescopes {\it XRISM} and {\it ATHENA}. If future detections of unexplained X-ray lines match our predictions, including line widths, we will have strong evidence that we have discovered the dark matter. △ Less

Submitted 11 October, 2018; originally announced October 2018.

Comments: 20 pages, 16 Figures. Highlights: Figs 4, 13, 15. To be submitted to MNRAS, comments welcome. Contact: lovell@hi.is

Baryon-induced dark matter cores in the EAGLE simulations

Authors: Alejandro Benitez-Llambay, Carlos S. Frenk, Aaron D. Ludlow, Julio F. Navarro

Abstract: We examine the formation of dark matter (DM) cores in dwarf galaxies simulated with the EAGLE model of galaxy formation. As in earlier work, we find that the star formation (SF) gas density threshold ($ρ_{\rm th}$) plays a critical role. At low thresholds (LT), gas is unable to reach densities high enough to dominate the gravitational potential before being dispersed by feedback from supernovae. L… ▽ More We examine the formation of dark matter (DM) cores in dwarf galaxies simulated with the EAGLE model of galaxy formation. As in earlier work, we find that the star formation (SF) gas density threshold ($ρ_{\rm th}$) plays a critical role. At low thresholds (LT), gas is unable to reach densities high enough to dominate the gravitational potential before being dispersed by feedback from supernovae. LT runs show little effect on the inner DM profile, even in systems with extended and bursty SF, two ingredients often cited as critical for core formation. For higher thresholds, gas is able to dominate the gravitational potential before being ejected by feedback. This can lead to a substantial reduction in the inner DM content, but only if the gas is gravitationally important over an extended period of time, allowing the halo to contract before gas removal. Rapid assembly and removal of gas in short SF bursts is less effective at altering the inner DM content. Subsequent gas accretion may draw DM back in and reform a cusp, unless SF is bursty enough to prevent it, preserving the core. Thus, for the EAGLE SF+feedback model, there is no simple relation between core formation and SF history, contrary to recent claims. The dependence of the inner DM content of dwarfs on $ρ_{\rm th}$ hinders robust predictions and the interpretation of observations. A simulation of a $(12 \rm \ Mpc)^3$ volume with high $ρ_{\rm th}$ results in dwarfs with sizeable cores over a limited halo mass range, but with insufficient variety in mass profiles to explain the observed diversity of dwarf galaxy rotation curves. △ Less

Submitted 23 August, 2019; v1 submitted 9 October, 2018; originally announced October 2018.

Comments: Updated to match the published version

First observation with global network of optical atomic clocks aimed for a dark matter detection

Authors: P. Wcisło, P. Ablewski, K. Beloy, S. Bilicki, M. Bober, R. Brown, R. Fasano, R. Ciuryło, H. Hachisu, T. Ido, J. Lodewyck, A. Ludlow, W. McGrew, P. Morzyński, D. Nicolodi, M. Schioppo, M. Sekido, R. Le Targat, P. Wolf, X. Zhang, B. Zjawin, M. Zawada

Abstract: We report on the first earth-scale quantum sensor network based on optical atomic clocks aimed at dark matter (DM) detection. Exploiting differences in the susceptibilities to the fine-structure constant of essential parts of an optical atomic clock, i.e. the cold atoms and the optical reference cavity, we can perform sensitive searches for dark matter signatures without the need of real-time comp… ▽ More We report on the first earth-scale quantum sensor network based on optical atomic clocks aimed at dark matter (DM) detection. Exploiting differences in the susceptibilities to the fine-structure constant of essential parts of an optical atomic clock, i.e. the cold atoms and the optical reference cavity, we can perform sensitive searches for dark matter signatures without the need of real-time comparisons of the clocks. We report a two orders of magnitude improvement in constraints on transient variations of the fine-structure constant, which considerably improves the detection limit for the standard model (SM) - DM coupling. We use Yb and Sr optical atomic clocks at four laboratories on three continents to search for both topological defect (TD) and massive scalar field candidates. No signal consistent with a dark-matter coupling is identified, leading to significantly improved constraints on the DM-SM couplings. △ Less

Submitted 12 June, 2018; originally announced June 2018.

Journal ref: Science Advances, 4, eaau4869 (2018)

The vertical structure of gaseous galaxy discs in cold dark matter halos

Authors: Alejandro Benitez-Llambay, Julio F. Navarro, Carlos S. Frenk, Aaron D. Ludlow

Abstract: We study the vertical structure of polytropic, $P\propto ρ^Γ$, centrifugally-supported gaseous discs embedded in cold dark matter (CDM) halos. At fixed radius $R$, the shape of the vertical density profile depends only weakly on whether the disc is self-gravitating (SG) or not (NSG). The disc thickness, set by the midplane sound speed and circular velocity, $(c_s/V_c)R$, in the NSG case, and by th… ▽ More We study the vertical structure of polytropic, $P\propto ρ^Γ$, centrifugally-supported gaseous discs embedded in cold dark matter (CDM) halos. At fixed radius $R$, the shape of the vertical density profile depends only weakly on whether the disc is self-gravitating (SG) or not (NSG). The disc thickness, set by the midplane sound speed and circular velocity, $(c_s/V_c)R$, in the NSG case, and by the sound speed and surface density, $c_s^2/GΣ$, in SG discs, is smaller than either of these scales. SG discs are typically Toomre unstable, NSG discs are stable. Exponential discs in CDM halos with roughly flat circular velocity curves generally "flare" outwards. For the polytropic equation of state of the EAGLE simulations, discs whose mass and size match observational constraints are stable (NSG) for $M_d< 3\times 10^9\, M_\odot$ and unstable (SG) at higher masses, if fully gaseous. We test these analytic results using a set of idealized SPH simulations and find excellent agreement. Our results clarify the role of the gravitational softening on the thickness of simulated discs, and on the onset of radial instabilities. EAGLE low-mass discs are non-self-gravitating so the softening plays no role in their vertical structure. High-mass discs, on the other hand, are expected to be self-gravitating and unstable, and may be artificially thickened and stabilized unless gravity is well resolved. Simulations with spatial resolution high enough to not compromise the vertical structure of a disc also resolve the onset of their instabilities, but the converse is not true: resolving instabilities does not guarantee that the vertical structure is resolved. △ Less

Submitted 25 July, 2017; originally announced July 2017.

Comments: Submitted to MNRAS

The origin of the mass discrepancy-acceleration relation in $Λ$CDM

Authors: Julio F. Navarro, Alejandro Benítez-Llambay, Azadeh Fattahi, Carlos S. Frenk, Aaron D. Ludlow, Kyle A. Oman, Matthieu Schaller, Tom Theuns

Abstract: We examine the origin of the mass discrepancy--radial acceleration relation (MDAR) of disk galaxies. This is a tight empirical correlation between the disk centripetal acceleration and that expected from the baryonic component. The MDAR holds for most radii probed by disk kinematic tracers, regardless of galaxy mass or surface brightness. The relation has two characteristic accelerations; $a_0$, a… ▽ More We examine the origin of the mass discrepancy--radial acceleration relation (MDAR) of disk galaxies. This is a tight empirical correlation between the disk centripetal acceleration and that expected from the baryonic component. The MDAR holds for most radii probed by disk kinematic tracers, regardless of galaxy mass or surface brightness. The relation has two characteristic accelerations; $a_0$, above which all galaxies are baryon-dominated; and $a_{\rm min}$, an effective minimum aceleration probed by kinematic tracers in isolated galaxies. We use a simple model to show that these trends arise naturally in $Λ$CDM. This is because: (i) disk galaxies in $Λ$CDM form at the centre of dark matter haloes spanning a relatively narrow range of virial mass; (ii) cold dark matter halo acceleration profiles are self-similar and have a broad maximum at the centre, reaching values bracketed precisely by $a_{\rm min}$ and $a_0$ in that mass range; and (iii) halo mass and galaxy size scale relatively tightly with the baryonic mass of a galaxy in any successful $Λ$CDM galaxy formation model. Explaining the MDAR in $Λ$CDM does not require modifications to the cuspy inner mass profiles of dark haloes, although these may help to understand the detailed rotation curves of some dwarf galaxies and the origin of extreme outliers from the main relation. The MDAR is just a reflection of the self-similar nature of cold dark matter haloes and of the physical scales introduced by the galaxy formation process. △ Less

Submitted 26 June, 2017; v1 submitted 19 December, 2016; originally announced December 2016.

Comments: 8 pages, 3 figures, MNRAS submitted

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 Mass-Discrepancy Acceleration Relation: a Natural Outcome of Galaxy Formation in Cold Dark Matter halos

Authors: Aaron D. Ludlow, Alejandro Benitez-Llambay, Matthieu Schaller, Tom Theuns, Carlos S. Frenk, Richard Bower, Joop Schaye, Robert A. Crain, Julio F. Navarro, Azadeh Fattahi, Kyle A. Oman

Abstract: We analyze the total and baryonic acceleration profiles of a set of well-resolved galaxies identified in the EAGLE suite of hydrodynamic simulations. Our runs start from the same initial conditions but adopt different prescriptions for unresolved stellar and AGN feedback, resulting in diverse populations of galaxies by the present day. Some of them reproduce observed galaxy scaling relations, whil… ▽ More We analyze the total and baryonic acceleration profiles of a set of well-resolved galaxies identified in the EAGLE suite of hydrodynamic simulations. Our runs start from the same initial conditions but adopt different prescriptions for unresolved stellar and AGN feedback, resulting in diverse populations of galaxies by the present day. Some of them reproduce observed galaxy scaling relations, while others do not. However, regardless of the feedback implementation, all of our galaxies follow closely a simple relationship between the total and baryonic acceleration profiles, consistent with recent observations of rotationally supported galaxies. The relation has small scatter: different feedback implementations -- which produce different galaxy populations -- mainly shift galaxies along the relation, rather than perpendicular to it. Furthermore, galaxies exhibit a characteristic acceleration, $g_{\dagger}$, above which baryons dominate the mass budget, as observed. These observations, consistent with simple modified Newtonian dynamics, can be accommodated within the standard cold dark matter paradigm. △ Less

Submitted 11 April, 2017; v1 submitted 24 October, 2016; originally announced October 2016.

Comments: 6 pages, 3 figures, PRL in press

Journal ref: Phys. Rev. Lett. 118, 161103 (2017)

Einasto Profiles and the Dark Matter Power Spectrum

Authors: Aaron D. Ludlow, Raul E. Angulo

Abstract: We study the mass accretion histories (MAHs) and density profiles of dark matter halos using N-body simulations of self-similar gravitational clustering from scale-free power spectra, $P(k)\propto k^n$. We pay particular attention to the density profile curvature, which we characterize using the shape parameter, $α$, of an Einasto profile. In agreement with previous findings our results suggest th… ▽ More We study the mass accretion histories (MAHs) and density profiles of dark matter halos using N-body simulations of self-similar gravitational clustering from scale-free power spectra, $P(k)\propto k^n$. We pay particular attention to the density profile curvature, which we characterize using the shape parameter, $α$, of an Einasto profile. In agreement with previous findings our results suggest that, despite vast differences in their MAHs, the density profiles of virialized halos are remarkably alike. Nonetheless, clear departures from self-similarity are evident: for a given spectral index, $α$ increases slightly but systematically with "peak height", $ν\equivδ_{sc}/σ(M,z)$, regardless of mass or redshift. More importantly, however, the "$α-ν$" relation depends on $n$: the steeper the initial power spectrum, the more gradual the curvature of both the mean MAHs and mean density profiles. These results are consistent with previous findings connecting the shapes of halo mass profiles and MAHs and imply that dark matter halos are not structurally self-similar but, through the merger history, retain a memory of the linear density field from which they form. △ Less

Submitted 14 October, 2016; originally announced October 2016.

Comments: MNRAS letters, accepted

Wide binaries in ultra-faint galaxies: a window onto dark matter on the smallest scales

Authors: Jorge Peñarrubia, Aaron D. Ludlow, Julio Chanamé, Matthew G. Walker

Abstract: We carry out controlled $N$-body simulations that follow the dynamical evolution of binary stars in the dark matter (DM) haloes of ultra-faint dwarf spheroidals (dSphs). We find that wide binaries with semi-major axes $a\gtrsim a_t$ tend to be quickly disrupted by the tidal field of the halo. In smooth potentials the truncation scale, $a_t$, is mainly governed by (i) the mass enclosed within the d… ▽ More We carry out controlled $N$-body simulations that follow the dynamical evolution of binary stars in the dark matter (DM) haloes of ultra-faint dwarf spheroidals (dSphs). We find that wide binaries with semi-major axes $a\gtrsim a_t$ tend to be quickly disrupted by the tidal field of the halo. In smooth potentials the truncation scale, $a_t$, is mainly governed by (i) the mass enclosed within the dwarf half-light radius ($R_h$) and (ii) the slope of the DM halo profile at $R\approx R_h$, and is largely independent of the initial eccentricity distribution of the binary systems and the anisotropy of the stellar orbits about the galactic potential. For the reported velocity dispersion and half-light radius of Segue I, the closest ultra-faint, our models predict $a_t$ values that are a factor 2--3 smaller in cuspy haloes than in any of the cored models considered here. Using mock observations of Segue I we show that measuring the projected two-point correlation function of stellar pairs with sub-arcsecond resolution may provide a useful tool to constrain the amount and distribution of DM in the smallest and most DM-dominated galaxies. △ Less

Submitted 30 May, 2016; originally announced May 2016.

Comments: Accepted to MNRAS Letters

Earth-mass haloes and the emergence of NFW density profiles

Authors: Raul E. Angulo, Oliver Hahn, Aaron Ludlow, Silvia Bonoli

Abstract: We simulate neutralino dark matter ($χ$DM) haloes from their initial collapse, at $\sim$ earth mass, up to a few percent solar. Our results confirm that the density profiles of the first haloes are described by a $\sim r^{-1.5}$ power-law. As haloes grow in mass, their density profiles evolve significantly. In the central regions, they become shallower and reach on average $\sim r^{-1}$, the asymp… ▽ More We simulate neutralino dark matter ($χ$DM) haloes from their initial collapse, at $\sim$ earth mass, up to a few percent solar. Our results confirm that the density profiles of the first haloes are described by a $\sim r^{-1.5}$ power-law. As haloes grow in mass, their density profiles evolve significantly. In the central regions, they become shallower and reach on average $\sim r^{-1}$, the asymptotic form of an NFW profile. Using non-cosmological controlled simulations, we observe that temporal variations in the gravitational potential caused by major mergers lead to a shallowing of the inner profile. This transformation is more significant for shallower initial profiles and for a higher number of merging systems. Depending on the merger details, the resulting profiles can be shallower or steeper than NFW in their inner regions. Interestingly, mergers have a much weaker effect when the profile is given by a broken power-law with an inner slope of $-1$ (such as NFW or Hernquist profiles). This offers an explanation for the emergence of NFW-like profiles: after their initial collapse, $r^{-1.5}$ $χ$DM haloes suffer copious major mergers, which progressively shallows the profile. Once an NFW-like profile is established, subsequent merging do not change the profile anymore. This suggests that halo profiles are not universal but rather a combination of (1) the physics of the formation of the microhaloes and (2) their early merger history -- both set by the properties of the dark matter particle -- as well as (3) the resilience of NFW-like profiles to perturbations. △ Less

Submitted 23 June, 2017; v1 submitted 11 April, 2016; originally announced April 2016.

Comments: Matches version accepted by MNRAS. 16 pages, 14 figures

The Mass-Concentration-Redshift Relation of Cold and Warm Dark Matter Halos

Authors: Aaron D. Ludlow, Sownak Bose, Raúl E. Angulo, Lan Wang, Wojciech A. Hellwing, Julio F. Navarro, Shaun Cole, Carlos S. Frenk

Abstract: We use a suite of cosmological simulations to study the mass-concentration-redshift relation, $c({\rm M},z)$, of dark matter halos. Our simulations include standard $Λ$-cold dark matter (CDM) models, and additional runs with truncated power spectra, consistent with a thermal warm dark matter (WDM) scenario. We find that the mass profiles of CDM and WDM halos are self-similar and well approximated… ▽ More We use a suite of cosmological simulations to study the mass-concentration-redshift relation, $c({\rm M},z)$, of dark matter halos. Our simulations include standard $Λ$-cold dark matter (CDM) models, and additional runs with truncated power spectra, consistent with a thermal warm dark matter (WDM) scenario. We find that the mass profiles of CDM and WDM halos are self-similar and well approximated by the Einasto profile. The $c({\rm M},z)$ relation of CDM halos is monotonic: concentrations decrease with increasing virial mass at fixed redshift, and decrease with increasing redshift at fixed mass. The mass accretion histories (MAHs) of CDM halos are also scale-free, and can be used to infer concentrations directly. These results do not apply to WDM halos: their MAHs are not scale-free because of the characteristic scale imposed by the power-spectrum suppression. Further, the WDM $c({\rm M},z)$ relation is non-monotonic: concentrations peak at a mass scale dictated by the truncation scale, and decrease at higher and lower masses. We show that the assembly history of a halo can still be used to infer its concentration, provided that the total mass of its progenitors is considered (the "collapsed mass history"; CMH), rather than just that of its main ancestor. This exploits the scale-free nature of CMHs to derive a simple scaling that reproduces the mass-concentration-redshift relation of both CDM and WDM halos over a vast range of halo masses and redshifts. Our model therefore provides a robust account of the mass, redshift, cosmology and power spectrum dependence of dark matter halo concentrations. △ Less

Submitted 10 May, 2016; v1 submitted 11 January, 2016; originally announced January 2016.

Comments: 21 pages, 16 figures

Dark-ages Reionization & Galaxy Formation Simulation II: Spin and concentration parameters for dark matter haloes during the Epoch of Reionization

Authors: Paul W. Angel, Gregory B. Poole, Aaron D. Ludlow, Alan R. Duffy, Paul M. Geil, Simon J. Mutch, Andrei Mesinger, J. Stuart B. Wyithe

Abstract: We use high resolution N-Body simulations to study the concentration and spin parameters of dark matter haloes in the mass range $10^8\, {\rm M}_{\odot}\, h^{-1} < {\rm M} < 10^{11}\, {\rm M}_{\odot}\, h^{-1}$ and redshifts $5{<}z{<}10$, corresponding to the haloes of galaxies thought to be responsible for reionization. We build a sub-sample of equilibrium haloes and contrast their properties to… ▽ More We use high resolution N-Body simulations to study the concentration and spin parameters of dark matter haloes in the mass range $10^8\, {\rm M}_{\odot}\, h^{-1} < {\rm M} < 10^{11}\, {\rm M}_{\odot}\, h^{-1}$ and redshifts $5{<}z{<}10$, corresponding to the haloes of galaxies thought to be responsible for reionization. We build a sub-sample of equilibrium haloes and contrast their properties to the full population that also includes unrelaxed systems. Concentrations are calculated by fitting both NFW and Einasto profiles to the spherically-averaged density profiles of individual haloes. After removing haloes that are out-of-equilibrium, we find a $z{>}5$ concentration$-$mass ($c(M)$) relation that is almost flat and well described by a simple power-law for both NFW and Einasto fits. The intrinsic scatter around the mean relation is $Δc_{\rm{vir}}{\sim1}$ (or 20 per cent) at $z=5$. We also find that the analytic model proposed by Ludlow et al. reproduces the mass and redshift-dependence of halo concentrations. Our best-fit Einasto shape parameter, $α$, depends on peak height, $ν$, in a manner that is accurately described by $α{=}0.0070ν^2{+}0.1839$. The distribution of the spin parameter, $λ$, has a weak dependence on equilibrium state; $λ$ peaks at roughly ${\sim}0.033$ for our relaxed sample, and at ${\sim}0.04$ for the full population. The spin--virial mass relation has a mild negative correlation at high redshift. △ Less

Submitted 1 December, 2015; originally announced December 2015.

The formation of CDM haloes II: collapse time and tides

Authors: Mikolaj Borzyszkowski, Aaron D. Ludlow, Cristiano Porciani

Abstract: We use two cosmological simulations of structure formation in the LambdaCDM scenario to study the evolutionary histories of dark-matter haloes and to characterize the Lagrangian regions from which they form. We focus on haloes identified at redshift z_id=0 and show that the classic ellipsoidal collapse model systematically overestimates their collapse times. If one imposes that halo collapse takes… ▽ More We use two cosmological simulations of structure formation in the LambdaCDM scenario to study the evolutionary histories of dark-matter haloes and to characterize the Lagrangian regions from which they form. We focus on haloes identified at redshift z_id=0 and show that the classic ellipsoidal collapse model systematically overestimates their collapse times. If one imposes that halo collapse takes place at z_id, this model requires starting from a significantly lower linear density contrast than what is measured in the simulations at the locations of halo formation. We attempt to explain this discrepancy by testing two key assumptions of the model. First, we show that the tides felt by collapsing haloes due to the surrounding large-scale structure evolve non-linearly. Although this effect becomes increasingly important for low-mass haloes, accounting for it in the ellipsoidal collapse model only marginally improves the agreement with N-body simulations. Second, we track the time evolution of the physical volume occupied by forming haloes and show that, after turnaround, it generally stabilizes at a well-defined redshift, z_c>z_id, contrary to the basic assumption of extended Press-Schechter theory based on excursion sets. We discuss the implications of this result for understanding the origin of the mass-dependence and scatter in the linear threshold for halo formation. Finally, we show that, when tuned for collapse at z_c, a modified version of the ellipsoidal collapse model that also accounts for the triaxial nature of protohaloes predicts their linear density contrast in an unbiased way. △ Less

Submitted 7 November, 2014; v1 submitted 28 May, 2014; originally announced May 2014.

Comments: 15 pages, 11 figures, MNRAS in press

Journal ref: Monthly Notices of the Royal Astronomical Society 2014 445 (3): 4124-4136

Atomic carbon as a powerful tracer of molecular gas in the high-redshift Universe: perspectives for ALMA

Authors: Matteo Tomassetti, Cristiano Porciani, Emilio Romano-Diaz, Aaron D. Ludlow, Padelis P. Papadopoulos

Abstract: We use a high-resolution hydrodynamic simulation that tracks the non-equilibrium abundance of molecular hydrogen within a massive high-redshift galaxy to produce mock Atacama Large Millimeter Array (ALMA) maps of the fine-structure lines of atomic carbon, CI 1-0 and CI 2-1. Inspired by recent observational and theoretical work, we assume that CI is thoroughly mixed within giant molecular clouds an… ▽ More We use a high-resolution hydrodynamic simulation that tracks the non-equilibrium abundance of molecular hydrogen within a massive high-redshift galaxy to produce mock Atacama Large Millimeter Array (ALMA) maps of the fine-structure lines of atomic carbon, CI 1-0 and CI 2-1. Inspired by recent observational and theoretical work, we assume that CI is thoroughly mixed within giant molecular clouds and demonstrate that its emission is an excellent proxy for H2. Nearly all of the H2 associated with the galaxy can be detected at redshifts z<4 using a compact interferometric configuration with a large synthesized beam (that does not resolve the target galaxy) in less than 4 h of integration time. Low-resolution imaging of the \CI lines (in which the target galaxy is resolved into three to four beams) will detect ~80 per cent of the H2 in less than 12 h of aperture synthesis. In this case, the resulting data cube also provides the crucial information necessary for determining the dynamical state of the galaxy. We conclude that ALMA observations of the CI 1-0 and 2-1 emission are well-suited for extending the interval of cosmic look-back time over which the H2 distributions, the dynamical masses, and the Tully-Fisher relation of galaxies can be robustly probed. △ Less

Submitted 9 September, 2014; v1 submitted 2 May, 2014; originally announced May 2014.

Comments: 5 pages, 4 figures, accepted for publication in MNRAS Letters

Simulating the H2 content of high-redshift galaxies

Authors: Matteo Tomassetti, Cristiano Porciani, Emilio Romano-Diaz, Aaron D. Ludlow

Abstract: We introduce a sub-grid model for the non-equilibrium abundance of molecular hydrogen in cosmological simulations of galaxy formation. We improve upon previous work by accounting for the unresolved structure of molecular clouds in a phenomenological way which combines both observational and numerical results on the properties of the turbulent interstellar medium. We apply the model to a cosmologic… ▽ More We introduce a sub-grid model for the non-equilibrium abundance of molecular hydrogen in cosmological simulations of galaxy formation. We improve upon previous work by accounting for the unresolved structure of molecular clouds in a phenomenological way which combines both observational and numerical results on the properties of the turbulent interstellar medium. We apply the model to a cosmological simulation of the formation of a Milky-Way-sized galaxy at z=2, and compare the results to those obtained using other popular prescriptions that compute the equilibrium abundance of H2. In these runs we introduce an explicit link between star formation and the local H2 abundance, and perform an additional simulation in which star formation is linked directly to the density of cold gas. In better agreement with observations, we find that the simulated galaxy produces less stars and harbors a larger gas reservoir when star formation is regulated by molecular hydrogen. In this case, the galaxy is composed of a younger stellar population as early star formation is inhibited in small, metal poor dark-matter haloes which cannot efficiently produce H2. The number of luminous satellites orbiting within the virial radius of the galaxy at z=2 is reduced by 10-30 per cent in models with H2-regulated star formation. △ Less

Submitted 8 December, 2014; v1 submitted 27 March, 2014; originally announced March 2014.

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

The Mass-Concentration-Redshift Relation of Cold Dark Matter Halos

Authors: Aaron D. Ludlow, Julio F. Navarro, Raul E. Angulo, Michael Boylan-Kolchin, Volker Springel, Carlos Frenk, Simon D. M. White

Abstract: We use the Millennium Simulation series to investigate the mass and redshift dependence of the concentration of equilibrium cold dark matter (CDM) halos. We extend earlier work on the relation between halo mass profiles and assembly histories to show how the latter may be used to predict concentrations for halos of all masses and at any redshift. Our results clarify the link between concentration… ▽ More We use the Millennium Simulation series to investigate the mass and redshift dependence of the concentration of equilibrium cold dark matter (CDM) halos. We extend earlier work on the relation between halo mass profiles and assembly histories to show how the latter may be used to predict concentrations for halos of all masses and at any redshift. Our results clarify the link between concentration and the ``collapse redshift'' of a halo as well as why concentration depends on mass and redshift solely through the dimensionless ``peak height'' mass parameter, $ν(M,z)=δ_{\rm crit}(z)/σ(M,z)$. We combine these results with analytic mass accretion histories to extrapolate the $c(M,z)$ relations to mass regimes difficult to reach through direct simulation. Our model predicts that, at given $z$, $c(M)$ should deviate systematically from a simple power law at high masses, where concentrations approach a constant value, and at low masses, where concentrations are substantially lower than expected from extrapolating published empirical fits. This correction may reduce the expected self-annihilation boost factor from substructure by about one order of magnitude. The model also reproduces the $c(M,z)$ dependence on cosmological parameters reported in earlier work, and thus provides a simple and robust account of the relation between cosmology and the mass-concentration-redshift relation of CDM halos. △ Less

Submitted 3 December, 2013; originally announced December 2013.

Comments: 13 pages, 10 figures, submitted to MNRAS

The Mass Profile and Accretion History of Cold Dark Matter Halos

Authors: Aaron D. Ludlow, Julio F. Navarro, Michael Boylan-Kolchin, Philip E. Bett, Raul E. Angulo, Ming Li, Simon D. M. White, Carlos Frenk, Volker Springel

Abstract: We use the Millennium Simulation series to study the relation between the accretion history (MAH) and mass profile of cold dark matter halos. We find that the mean density within the scale radius, r_{-2} (where the halo density profile has isothermal slope), is directly proportional to the critical density of the Universe at the time when the main progenitor's virial mass equals the mass enclosed… ▽ More We use the Millennium Simulation series to study the relation between the accretion history (MAH) and mass profile of cold dark matter halos. We find that the mean density within the scale radius, r_{-2} (where the halo density profile has isothermal slope), is directly proportional to the critical density of the Universe at the time when the main progenitor's virial mass equals the mass enclosed within r_{-2}. Scaled to these characteristic values of mass and density, the mean MAH, expressed in terms of the critical density of the Universe, M(ρ_{crit}(z)), resembles that of the enclosed density profile, M(<ρ>), at z=0. Both follow closely the NFW profile, suggesting that the similarity of halo mass profiles originates from the mass-independence of halo MAHs. Support for this interpretation is provided by outlier halos whose accretion histories deviate from the NFW shape; their mass profiles show correlated deviations from NFW and are better approximated by Einasto profiles. Fitting both M(<ρ>) and M(ρ_{crit}) with either NFW or Einasto profiles yield concentration and shape parameters that are correlated, confirming and extending earlier work linking the concentration of a halo with its accretion history. These correlations also confirm that halo structure is insensitive to initial conditions: only halos whose accretion histories differ greatly from the NFW shape show noticeable deviations from NFW in their mass profiles. As a result, the NFW profile provides acceptable fits to hot dark matter halos, which do not form hierarchically, and for fluctuation power spectra other than CDM. Our findings, however, predict a subtle but systematic dependence of mass profile shape on accretion history which, if confirmed, would provide strong support for the link between accretion history and halo structure we propose here. △ Less

Submitted 3 December, 2013; v1 submitted 1 February, 2013; originally announced February 2013.

Comments: 12 pages, 8 figures, MNRAS 432 1103L (2013)

Journal ref: MNRAS 432 1103L (2013)

The Dynamical State and Mass-Concentration Relation of Galaxy Clusters

Authors: Aaron D. Ludlow, Julio F. Navarro, Ming Li, Raul E. Angulo, Michael Boylan-Kolchin, Philip E. Bett

Abstract: We use the Millennium Simulation series to study how the dynamical state of dark matter halos affects the relation between mass and concentration. We find that a large fraction of massive systems are identified when they are substantially out of equilibrium and in a particular phase of their dynamical evolution: the more massive the halo, the more likely it is found at a transient stage of high co… ▽ More We use the Millennium Simulation series to study how the dynamical state of dark matter halos affects the relation between mass and concentration. We find that a large fraction of massive systems are identified when they are substantially out of equilibrium and in a particular phase of their dynamical evolution: the more massive the halo, the more likely it is found at a transient stage of high concentration. This state reflects the recent assembly of massive halos and corresponds to the first pericentric passage of recently-accreted material when, before virialization, the kinetic and potential energies reach maximum and minimum values, respectively. This result explains the puzzling upturn in the mass-concentration relation reported in recent work for massive halos; indeed, the upturn disappears when only dynamically-relaxed systems are considered in the analysis. Our results warn against applying simple equilibrium models to describe the structure of rare, massive galaxy clusters and urges caution when extrapolating scaling laws calibrated on lower-mass systems, where such deviations from equilibrium are less common. The evolving dynamical state of galaxy clusters ought to be carefully taken into account if cluster studies are to provide precise cosmological constraints. △ Less

Submitted 1 February, 2013; v1 submitted 5 June, 2012; originally announced June 2012.

Comments: 8 Pages. Minor changes to match published version

Journal ref: MNRAS 2012, 427, p1322-1328

The spatial and velocity bias of linear density peaks and proto-haloes in the Lambda cold dark matter cosmology

Authors: Anna Elia, Aaron D. Ludlow, Cristiano Porciani

Abstract: We use high resolution N-body simulations to investigate the Lagrangian bias of cold dark matter haloes within the LCDM cosmology. Our analysis focuses on "proto-haloes", which we identify in the simulation initial conditions with the subsets of particles belonging to individual redshift-zero haloes. We then calculate the number-density and velocity-divergence fields of proto-haloes and estimate t… ▽ More We use high resolution N-body simulations to investigate the Lagrangian bias of cold dark matter haloes within the LCDM cosmology. Our analysis focuses on "proto-haloes", which we identify in the simulation initial conditions with the subsets of particles belonging to individual redshift-zero haloes. We then calculate the number-density and velocity-divergence fields of proto-haloes and estimate their auto spectral densities. We also measure the corresponding cross spectral densities with the linear matter distribution. We use our results to test a Lagrangian-bias model presented by Desjacques and Sheth which is based on the assumption that haloes form out of local density maxima of a specific height. Our comparison validates the predicted functional form for the scale-dependence of the bias for both the density and velocity fields. We also show that the bias coefficients are accurately predicted for the velocity divergence. On the contrary, the theoretical values for the density bias parameters do not accurately match the numerical results as a function of halo mass. This is likely due to the simplistic assumptions that relate virialized haloes to density peaks of a given height in the model. We also detect appreciable stochasticity for the Lagrangian density bias, even on very large scales. These are not included in the model at leading order but correspond to higher order corrections. △ Less

Submitted 15 February, 2012; v1 submitted 17 November, 2011; originally announced November 2011.

Comments: 10 pages, 4 figures. Matches version accepted for publication in MNRAS

The formation of CDM haloes I: Collapse thresholds and the ellipsoidal collapse model

Authors: Aaron D. Ludlow, Mikolaj Borzyszkowski, Cristiano Porciani

Abstract: In the excursion set approach to structure formation initially spherical regions of the linear density field collapse to form haloes of mass $M$ at redshift $z_{\rm id}$ if their linearly extrapolated density contrast, averaged on that scale, exceeds some critical threshold, $δ_{\rm c}(z_{\rm id})$. The value of $δ_{\rm c}(z_{\rm id})$ is often calculated from the spherical or ellipsoidal collapse… ▽ More In the excursion set approach to structure formation initially spherical regions of the linear density field collapse to form haloes of mass $M$ at redshift $z_{\rm id}$ if their linearly extrapolated density contrast, averaged on that scale, exceeds some critical threshold, $δ_{\rm c}(z_{\rm id})$. The value of $δ_{\rm c}(z_{\rm id})$ is often calculated from the spherical or ellipsoidal collapse model, which provide well-defined predictions given auxiliary properties of the tidal field at a given location. We use two cosmological simulations of structure growth in a $Λ$ cold dark matter scenario to quantify $δ_{\rm c}(z_{\rm id})$, its dependence on the surrounding tidal field, as well as on the shapes of the Lagrangian regions that collapse to form haloes at $z_{\rm id}$. Our results indicate that the ellipsoidal collapse model provides an accurate description of the mean dependence of $δ_{\rm c}(z_{\rm id})$ on both the strength of the tidal field and on halo mass. However, for a given $z_{\rm id}$, $δ_{\rm c}(z_{\rm id})$ depends strongly on the halo's characteristic formation redshift: the earlier a halo forms, the higher its initial density contrast. Surprisingly, the majority of haloes forming $today$ fall below the ellipsoidal collapse barrier, contradicting the model predictions. We trace the origin of this effect to the non-spherical shapes of Lagrangian haloes, which arise naturally due to the asymmetry of the linear tidal field. We show that a modified collapse model, that accounts for the triaxial shape of protohaloes, provides a more accurate description of the measured minimum overdensities of recently collapsed objects. △ Less

Submitted 6 November, 2014; v1 submitted 28 July, 2011; originally announced July 2011.

Comments: MNRAS in press

Journal ref: Monthly Notices of the Royal Astronomical Society 2014 445 (3): 4110-4123

A large-scale R-matrix calculation for electron-impact excitation of the Ne$^{2+}$ O-like ion

Authors: B. M. McLaughlin, Teck-Ghee Lee, J. A. Ludlow, E. Landi, S. D. Loch, M. S. Pindzola, C. P. Ballance

Abstract: The five J$Π$ levels within a $np^2$ or $np^4$ ground state complex provide an excellent testing ground for the comparison of theoretical line ratios with astrophysically observed values, in addition to providing valuable electron temperature and density diagnostics. The low temperature nature of the line ratios ensure that the theoretically derived values are sensitive to the underlying atomic st… ▽ More The five J$Π$ levels within a $np^2$ or $np^4$ ground state complex provide an excellent testing ground for the comparison of theoretical line ratios with astrophysically observed values, in addition to providing valuable electron temperature and density diagnostics. The low temperature nature of the line ratios ensure that the theoretically derived values are sensitive to the underlying atomic structure and electron-impact excitation rates. Previous R-matrix calculations for the Ne$^{2+}$ O-like ion exhibit large spurious structure in the cross sections at higher electron energies, which may affect Maxwellian averaged rates even at low temperatures. Furthermore, there is an absence of comprehensive excitation data between the excited states that may provide newer diagnostics to compliment the more established lines discussed in this paper. To resolve these issues, we present both a small scale 56-level Breit-Pauli (BP) calculation and a large-scale 554 levels R-matrix Intermediate Coupling Frame Transformation (ICFT) calculation that extends the scope and validity of earlier JAJOM calculations both in terms of the atomic structure and scattering cross sections. Our results provide a comprehensive electron-impact excitation data set for all transitions to higher $n$ shells. The fundamental atomic data for this O-like ion is subsequently used within a collisional radiative framework to provide the line ratios across a range of electron temperatures and densities of interest in astrophysical observations. △ Less

Submitted 21 July, 2011; v1 submitted 24 February, 2011; originally announced February 2011.

Comments: 17 pages, 8 figures

Journal ref: J.Phys.B: At. Mol. Opt. Phys. 44 (2011) 175206

The Density and Pseudo-Phase-Space Density Profiles of CDM halos

Authors: Aaron D. Ludlow, Julio F. Navarro, Michael Boylan-Kolchin, Volker Springel, Adrian Jenkins, Carlos S. Frenk, Simon D. M. White

Abstract: Cosmological N-body simulations indicate that the spherically-averaged density profiles of cold dark matter halos are accurately described by Einasto profiles, where the logarithmic slope is a power-law of adjustable exponent, γ=dlnρ/dlnr ~ r^α$. The pseudo-phase-space density (PPSD) profiles of CDM halos also show remarkable regularity, and are well approximated by simple power laws, Q(r)=ρ/σ^3 ~… ▽ More Cosmological N-body simulations indicate that the spherically-averaged density profiles of cold dark matter halos are accurately described by Einasto profiles, where the logarithmic slope is a power-law of adjustable exponent, γ=dlnρ/dlnr ~ r^α$. The pseudo-phase-space density (PPSD) profiles of CDM halos also show remarkable regularity, and are well approximated by simple power laws, Q(r)=ρ/σ^3 ~ r^-χ. We show that this is expected from dynamical equilibrium considerations, since Jeans' equations predict that the pseudo-phase-space density profiles of Einasto halos should resemble power laws over a wide range of radii. For the values of αtypical of CDM halos, the inner Q profiles of equilibrium halos deviate significantly from a power law only very close to the center, and simulations of extremely high-resolution would be needed to detect such deviations unambiguously. We use an ensemble of halos drawn from the Millennium-II simulation to study which of these two alternatives describe best the mass profile of CDM halos. Our analysis indicates that at the resolution of the best available simulations, both Einasto and power-law PPSD profiles (with adjustable exponents αand χ, respectively) provide equally acceptable fits to the simulations. A full account of the structure of CDM halos requires understanding how the shape parameters that characterize departures from self-similarity, like αor χ, are determined by evolutionary history, environment or initial conditions. △ Less

Submitted 31 January, 2011; originally announced February 2011.

Comments: 9 pages, 7 figures, submitted to MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.415:3895-3902,2011

The Peaks Formalism and the Formation of Cold Dark Matter Haloes

Authors: Aaron D. Ludlow, Cristiano Porciani

Abstract: We use two cosmological simulations of structure formation to study the conditions under which dark matter haloes emerge from the linear density field. Our analysis focuses on matching sites of halo collapse to local density maxima, or "peaks", in the initial conditions of the simulations and provides a crucial test of the central ansatz of the peaks formalism. By identifying peaks on a variety of… ▽ More We use two cosmological simulations of structure formation to study the conditions under which dark matter haloes emerge from the linear density field. Our analysis focuses on matching sites of halo collapse to local density maxima, or "peaks", in the initial conditions of the simulations and provides a crucial test of the central ansatz of the peaks formalism. By identifying peaks on a variety of smoothed, linearly extrapolated density fields we demonstrate that as many as ~70% of well-resolved dark matter haloes form preferentially near peaks whose characteristic masses are similar to that of the halo, with more massive haloes showing a stronger tendency to reside near peaks initially. We identify a small but significant fraction of haloes that appear to evolve from peaks of substantially lower mass than that of the halo itself. We refer to these as "peakless haloes" for convenience. By contrasting directly the properties of these objects with the bulk of the proto-halo population we find two clear differences: 1) their initial shapes are significantly flatter and more elongated than the predominantly triaxial majority, and 2) they are, on average, more strongly compressed by tidal forces associated with their surrounding large scale structure. Using the two-point correlation function we show that peakless haloes tend to emerge from highly clustered regions of the initial density field implying that, at fixed mass, the accretion geometry and mass accretion histories of haloes in highly clustered environments differ significantly from those in the field. This may have important implications for understanding the origin of the halo assembly bias, of galaxy properties in dense environments and how environment affects the morphological transformation of galaxies near groups and rich galaxy clusters. △ Less

Submitted 17 March, 2011; v1 submitted 10 November, 2010; originally announced November 2010.

Comments: 13 pages, 11 figures, published in MNRAS

Assembly History and Structure of Galactic Cold Dark Matter Halos

Authors: Jie Wang, J. F. Navarro, C. S. Frenk, S. D. M. White, V. Springel, A. Jenkins, A. Helmi, A. Ludlow, M. Vogelsberger

Abstract: We use the Aquarius simulation series to study the imprint of assembly history on the structure of Galaxy-mass cold dark matter halos. Our results confirm earlier work regarding the influence of mergers on the mass density profile and the inside-out growth of halos. The inner regions that contain the visible galaxies are stable since early times and are significantly affected only by major mergers… ▽ More We use the Aquarius simulation series to study the imprint of assembly history on the structure of Galaxy-mass cold dark matter halos. Our results confirm earlier work regarding the influence of mergers on the mass density profile and the inside-out growth of halos. The inner regions that contain the visible galaxies are stable since early times and are significantly affected only by major mergers. Particles accreted diffusely or in minor mergers are found predominantly in the outskirts of halos. Our analysis reveals trends that run counter to current perceptions of hierarchical halo assembly. For example, major mergers (i.e. those with progenitor mass ratios greater than 1:10) contribute little to the total mass growth of a halo, on average less than 20 per cent for our six Aquarius halos. The bulk is contributed roughly equally by minor mergers and by "diffuse" material which is not resolved into individual objects. This is consistent with modeling based on excursion-set theory which suggests that about half of this diffuse material should not be part of a halo of any scale. Interestingly, the simulations themselves suggest that a significantly fraction is not truly diffuse, since it was ejected from earlier halos by mergers prior to their joining the main system. The Aquarius simulations resolve halos to much lower mass scales than are expected to retain gas or form stars. These results thus confirm that most of the baryons from which visible galaxies form are accreted diffusely, rather than through mergers, and they suggest that only relatively rare major mergers will affect galaxy structure at later times. △ Less

Submitted 2 September, 2010; v1 submitted 30 August, 2010; originally announced August 2010.

Comments: 11 pages, 10 figs, submitted to MNRAS, a superfluous paragraph at the end of sec.2 taken out, the position of plots rearranged

Secondary Infall and the Pseudo-Phase-Space Density Profiles of Cold Dark Matter Halos

Authors: Aaron D. Ludlow, Julio F. Navarro, Volker Springel, Mark Vogelsberger, Jie Wang, Simon D. M. White, Adrian Jenkins, Carlos S. Frenk

Abstract: We use N-body simulations to investigate the radial dependence of the density and velocity dispersion in cold dark matter (CDM) halos. In particular, we explore how closely Q rho/sigma^3, a surrogate measure of the phase-space density, follows a power-law in radius. Our study extends earlier work by considering, in addition to spherically-averaged profiles, local Q-estimates for individual parti… ▽ More We use N-body simulations to investigate the radial dependence of the density and velocity dispersion in cold dark matter (CDM) halos. In particular, we explore how closely Q rho/sigma^3, a surrogate measure of the phase-space density, follows a power-law in radius. Our study extends earlier work by considering, in addition to spherically-averaged profiles, local Q-estimates for individual particles, Q_i; profiles based on the ellipsoidal radius dictated by the triaxial structure of the halo, Q_i(r'); and by carefully removing substructures in order to focus on the profile of the smooth halo, Q^s. The resulting Q_i^s(r') profiles follow closely a power law near the center, but show a clear upturn from this trend near the virial radius, r_{200}. The location and magnitude of the deviations are in excellent agreement with the predictions from Bertschinger's spherical secondary-infall similarity solution. In this model, Q \propto r^{-1.875} in the inner, virialized regions, but departures from a power-law occur near r_{200} because of the proximity of this radius to the location of the first shell crossing - the shock radius in the case of a collisional fluid. Particles there have not yet fully virialized, and so Q departs from the inner power-law profile. Our results imply that the power-law nature of $Q$ profiles only applies to the inner regions and cannot be used to predict accurately the structure of CDM halos beyond their characteristic scale radius. △ Less

Submitted 13 January, 2010; originally announced January 2010.

Comments: 10 pages, 8 figures, submitted to MNRAS

Phase-space structure in the local dark matter distribution and its signature in direct detection experiments

Authors: Mark Vogelsberger, Amina Helmi, Volker Springel, Simon D. M. White, Jie Wang, Carlos S. Frenk, Adrian Jenkins, Aaron Ludlow, Julio F. Navarro

Abstract: We study predictions for dark matter phase-space structure near the Sun based on high-resolution simulations of six galaxy halos taken from the Aquarius Project. The local DM density distribution is predicted to be remarkably smooth; the density at the Sun differs from the mean over a best-fit ellipsoidal equidensity contour by less than 15% at the 99.9% confidence level. The local velocity dist… ▽ More We study predictions for dark matter phase-space structure near the Sun based on high-resolution simulations of six galaxy halos taken from the Aquarius Project. The local DM density distribution is predicted to be remarkably smooth; the density at the Sun differs from the mean over a best-fit ellipsoidal equidensity contour by less than 15% at the 99.9% confidence level. The local velocity distribution is also very smooth, but it differs systematically from a (multivariate) Gaussian distribution. This is not due to the presence of individual clumps or streams, but to broad features in the velocity modulus and energy distributions that are stable both in space and time and reflect the detailed assembly history of each halo. These features have a significant impact on the signals predicted for WIMP and axion searches. For example, WIMP recoil rates can deviate by ~10% from those expected from the best-fit multivariate Gaussian models. The axion spectra in our simulations typically peak at lower frequencies than in the case of multivariate Gaussian velocity distributions. Also in this case, the spectra show significant imprints of the formation of the halo. This implies that once direct DM detection has become routine, features in the detector signal will allow us to study the dark matter assembly history of the Milky Way. A new field, "dark matter astronomy", will then emerge. △ Less

Submitted 2 December, 2008; originally announced December 2008.

Comments: 15 pages, 15 figures, submitted to MNRAS. Pictures, movies, and further information on the Aquarius Project may be found at http://www.mpa-garching.mpg.de/aquarius/