Source author record

Simon D. M. White

Simon D. M. White appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.CO astro-ph.GA astro-ph astro-ph.HE hep-ph hep-th physics.comp-ph

Catalog footprint

What is connected

104works

7topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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

Building this map preview

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

preprint2022arXiv

Dark matter annihilation and the Galactic Centre Excess

We compare the surface brightness profile and morphology of the Galactic Centre Excess (GCE) identified in wide-angle $γ$-ray maps from the Fermi-Large Area Telescope to dark matter annihilation predictions derived from high-resolution $Λ$CDM magnetohydrodynamic simulations of galaxy formation. These simulations produce isolated, disc-dominated galaxies with structure, stellar populations, gas content, and stellar and halo masses comparable to those of the Milky Way. For a specific choice of annihilation cross-section, they agree well with the Fermi-LAT data over the full observed angular range, $1^{\circ}$ to $15^{\circ}$, whereas their dark-matter only counterparts, lacking any compression of the inner halo by the gravitational effects of the baryons, fail to predict emission as centrally concentrated as observed. These results provide additional support to the hypothesis that the GCE is produced by annihilating dark matter. If, however, it is produced by a different mechanism, they imply a strong upper limit on annihilation rates which can be translated into upper limits on the expected $γ$-ray flux not only from the inner Galaxy but also from any substructure, with or without stars, in the Galactic halo.

preprint2022arXiv

Disc instability and bar formation: view from the IllustrisTNG simulations

We make use of z = 0 samples of strongly barred and unbarred disc galaxies from the TNG100 and TNG50 cosmological hydrodynamical simulations to assess the performance of the simple disc instability criterion proposed by Efstathiou, Lake & Negroponte (1982) (ELN-criterion). We find that strongly barred galaxies generally assemble earlier, are more star-dominated in their central regions, and have more massive and more compact discs than unbarred galaxies. The ELN-criterion successfully identifies ~75% and ~80% of the strongly barred and the unbarred galaxies, respectively. Strongly barred galaxies that the criterion fails to identify tend to have more extended discs, higher spin values and bars that assembled later than is typical for the bulk of the barred population. The bars in many of these cases appear to be produced by an interaction with a close neighbour (i.e. to be externally triggered) rather than to result from secular growth in the disc. On the other hand, we find that unbarred galaxies misclassified as barred by the ELN-criterion typically have stellar discs similar to those of barred galaxies, although more extended in the vertical direction and less star-dominated in their central regions, possibly reflecting later formation times. In addition, the bulge component of these galaxies is significantly more prominent at early times than in the strongly barred sample. Thus, the ELN-criterion robustly identifies secular bar instabilities in most simulated disc galaxies, but additional environmental criteria are needed to account for interaction-induced bar formation.

preprint2020arXiv

Comparing galaxy formation in the L-GALAXIES semi-analytical model and the IllustrisTNG simulations

We perform a comparison, object-by-object and statistically, between the Munich semi-analytical model, L-Galaxies, and the IllustrisTNG hydrodynamical simulations. By running L-Galaxies on the IllustrisTNG dark matter-only merger trees, we identify the same galaxies in the two models. This allows us to compare the stellar mass, star formation rate and gas content of galaxies, as well as the baryonic content of subhaloes and haloes in the two models. We find that both the stellar mass functions and the stellar masses of individual galaxies agree to better than $\sim0.2\,$dex. On the other hand, specific star formation rates and gas contents can differ more substantially. At $z=0$ the transition between low-mass star-forming galaxies and high-mass, quenched galaxies occurs at a stellar mass scale $\sim0.5\,$dex lower in IllustrisTNG than in L-Galaxies. IllustrisTNG also produces substantially more quenched galaxies at higher redshifts. Both models predict a halo baryon fraction close to the cosmic value for clusters, but IllustrisTNG predicts lower baryon fractions in group environments. These differences are due primarily to differences in modelling feedback from stars and supermassive black holes. The gas content and star formation rates of galaxies in and around clusters and groups differ substantially, with IllustrisTNG satellites less star-forming and less gas-rich. We show that environmental processes such as ram-pressure stripping are stronger and operate to larger distances and for a broader host mass range in IllustrisTNG. We suggest that the treatment of galaxy evolution in the semi-analytic model needs to be improved by prescriptions which capture local environmental effects more accurately.

preprint2020arXiv

Measuring the Tidal Response of Structure Formation: Anisotropic Separate Universe Simulations using TreePM

We present anisotropic "separate universe" simulations which modify the N-body code Gadget-4 in order to represent a large-scale tidal field through an anisotropic expansion factor. These simulations are used to measure the linear, quasi-linear and nonlinear response of the matter power spectrum to a spatially uniform trace-free tidal field up to wavenumber $k = 7 h \text{Mpc}^{-1}$. Together with the response to a large-scale overdensity measured in previous work, this completely describes the nonlinear matter bispectrum in the squeezed limit. We find that the response amplitude does not approach zero on small scales in physical coordinates, but rather a constant value at $z=0$, $R_K\approx 0.5$ for $k \geq 3 h\text{Mpc}^{-1}$ up to the scale where we consider our simulations reliable, $k \leq 7 h\text{Mpc}^{-1}$ at $z=0$. This shows that even the inner regions of haloes are affected by the large-scale tidal field. We also measure directly the alignment of halo shapes with the tidal field, finding a clear signal which increases with halo mass.

preprint2020arXiv

Simulating the Complexity of the Dark Matter Sheet I: Numerical Algorithms

At early times dark matter has a thermal velocity dispersion of unknown amplitude which, for warm dark matter models, can influence the formation of nonlinear structure on observable scales. We propose a new scheme to simulate cosmologies with a small-scale suppression of perturbations that combines two previous methods in a way that avoids the numerical artefacts which have so far prevented either from producing fully reliable results. At low densities and throughout most of the cosmological volume, we represent the dark matter phase-sheet directly using high-accuracy interpolation, thereby avoiding the artificial fragmentation which afflicts particle-based methods in this regime. Such phase-sheet methods are, however, unable to follow the rapidly increasing complexity of the denser regions of dark matter haloes, so for these we switch to an N-body scheme which uses the geodesic deviation equation to track phase-sheet properties local to each particle. In addition, we present a novel high-resolution force calculation scheme based on an oct-tree of cubic force resolution elements which is well suited to approximate the force-field of our combined sheet+particle distribution. Our hybrid simulation scheme enables the first reliable simulations of the internal structure of low-mass haloes in a warm dark matter cosmology.

preprint2020arXiv

The globular cluster system of the Auriga simulations

We investigate whether the galaxy and star formation model used for the Auriga simulations can produce a realistic globular cluster (GC) population. We compare statistics of GC candidate star particles in the Auriga haloes with catalogues of the Milky Way (MW) and Andromeda (M31) GC populations. We find that the Auriga simulations do produce sufficient stellar mass for GC candidates at radii and metallicities that are typical for the MW GC system (GCS). We also find varying mass-ratios of the simulated GC candidates relative to the observed mass in the MW and M31 GC systems for different bins of galactocentric radius-metallicity (r$_{\text{gal}}$ -[Fe/H]). Overall, the Auriga simulations produce GC candidates with higher metallicities than the MW and M31 GCS and they are found at larger radii than observed. The Auriga simulations would require bound cluster formation efficiencies higher than ten percent for the metal-poor GC candidates, and those within the Solar radius should experience negligible destruction rates to be consistent with observations. GC candidates in the outer halo, on the other hand, should either have low formation efficiencies, or experience high mass loss for the Auriga simulations to produce a GCS that is consistent with that of the MW or M31. Finally, the scatter in the metallicity as well as in the radial distribution between different Auriga runs is considerably smaller than the differences between that of the MW and M31 GCSs. The Auriga model is unlikely to give rise to a GCS that can be consistent with both galaxies.

preprint2020arXiv

The Tessellation-Level-Tree: characterising the nested hierarchy of density peaks and their spatial distribution in cosmological N-body simulations

We use the Millennium and Millennium-II simulations to illustrate the Tessellation-Level-Tree (TLT), a hierarchical tree structure linking density peaks in a field constructed by voronoi tessellation of the particles in a cosmological N-body simulation. The TLT uniquely partitions the simulation particles into disjoint subsets, each associated with a local density peak. Each peak is a subpeak of a unique higher peak. The TLT can be persistence filtered to suppress peaks produced by discreteness noise. Thresholding a peak's particle list at $\sim 80\left<ρ\right>$ results in a structure similar to a standard friend-of-friends halo and its subhaloes. For thresholds below $\sim 7\left<ρ\right>$, the largest structure percolates and is much more massive than other objects. It may be considered as defining the cosmic web. For a threshold of $5\left<ρ\right>$, it contains about half of all cosmic mass and occupies $\sim 1\%$ of all cosmic volume; a typical external point is then $\sim 7h^{-1}\mathrm{Mpc}$ from the web. We investigate the internal structure and clustering of TLT peaks. Defining the saddle point density $ρ_{\mathrm{lim}}$ as the density at which a peak joins its parent peak, we show the median value of $ρ_{\mathrm{lim}}$ for FoF-like peaks to be similar to the density threshold at percolation. Assembly bias as a function of $ρ_{\mathrm{lim}}$ is stronger than for any known internal halo property. For peaks of group mass and below, the lowest quintile in $ρ_{\mathrm{lim}}$ has $b\approx 0$, and is thus uncorrelated with the mass distribution.

preprint2020arXiv

Universal structure of dark matter haloes over a mass range of 20 orders of magnitude

Cosmological models in which dark matter consists of cold elementary particles predict that the dark halo population should extend to masses many orders of magnitude below those at which galaxies can form. Here we report a cosmological simulation of the formation of present-day haloes over the full range of observed halo masses (20 orders of magnitude) when dark matter is assumed to be in the form of weakly interacting massive particles of mass approximately 100 gigaelectronvolts. The simulation has a full dynamic range of 30 orders of magnitude in mass and resolves the internal structure of hundreds of Earth-mass haloes in as much detail as it does for hundreds of rich galaxy clusters. We find that halo density profiles are universal over the entire mass range and are well described by simple two-parameter fitting formulae. Halo mass and concentration are tightly related in a way that depends on cosmology and on the nature of the dark matter. For a fixed mass, the concentration is independent of the local environment for haloes less massive than those of typical galaxies. Haloes over the mass range of 10^3 to 10^11 solar masses contribute about equally (per logarithmic interval) to the luminosity produced by dark matter annihilation, which we find to be smaller than all previous estimates by factors ranging up to one thousand.

preprint2016arXiv

A chronicle of galaxy mass assembly in the EAGLE simulation

We analyse the mass assembly of central galaxies in the EAGLE hydrodynamical simulations. We build merger trees to connect galaxies to their progenitors at different redshifts and characterize their assembly histories by focusing on the time when half of the galaxy stellar mass was assembled into the main progenitor. We show that galaxies with stellar mass $M_*<10^{10.5}M_{\odot}$ assemble most of their stellar mass through star formation in the main progenitor (`in-situ' star formation). This can be understood as a consequence of the steep rise in star formation efficiency with halo mass for these galaxies. For more massive galaxies, however, an increasing fraction of their stellar mass is formed outside the main progenitor and subsequently accreted. Consequently, while for low-mass galaxies the assembly time is close to the stellar formation time, the stars in high-mass galaxies typically formed long before half of the present-day stellar mass was assembled into a single object, giving rise to the observed anti-hierarchical downsizing trend. In a typical present-day $M_*\geq10^{11}M_{\odot}$ galaxy, around $20\%$ of the stellar mass has an external origin. This fraction decreases with increasing redshift. Bearing in mind that mergers only make an important contribution to the stellar mass growth of massive galaxies, we find that the dominant contribution comes from mergers with galaxies of mass greater than one tenth of the main progenitor's mass. The galaxy merger fraction derived from our simulations agrees with recent observational estimates.

preprint2016arXiv

A fully cosmological model of a Monoceros-like ring

We study the vertical structure of a stellar disk obtained from a fully cosmological high-resolution hydrodynamical simulation of the formation of a Milky Way-like galaxy. At the present day, the disk's mean vertical height shows a well-defined and strong pattern, with amplitudes as large as 3 kpc in its outer regions. This pattern is the result of a satellite - host halo - disk interaction and reproduces, qualitatively, many of the observable properties of the Monoceros Ring. In particular we find disk material at the distance of Monoceros ($R \sim$ 12-16 kpc, galactocentric) extending far above the mid plane ( 30$^{\circ}$, $\langle Z \rangle \sim$ 1-2 kpc) in both hemispheres, as well as well-defined arcs of disk material at heliocentric distances $\gtrsim 5$ kpc. The pattern was first excited $\approx 3$ Gyr ago as an $m=1$ mode that later winds up into a leading spiral pattern. Interestingly, the main driver behind this perturbation is a low-mass low-velocity fly-by encounter. The satellite has total mass, pericentre distance and pericentric velocity of $\sim 5\%$ of the host, $\sim 80$ kpc, and 215 km/s, respectively. The satellite is not massive enough to directly perturb the galactic disk but we show that the density field of the host dark matter halo responds to this interaction resulting in a strong amplification of the perturbative effects. This subsequently causes the onset and development of the Monoceros-like feature.

preprint2016arXiv

Gaps in globular cluster streams: giant molecular clouds can cause them too

As a result of their internal dynamical coherence, thin stellar streams formed by disrupting globular clusters (GCs) can act as detectors of dark matter (DM) substructure in the Galactic halo. Perturbations induced by close flybys amplify into detectable density gaps, providing a probe both of the abundance and of the masses of DM subhaloes. Here, we use N-body simulations to show that the Galactic population of giant molecular clouds (GMCs) can also produce gaps (and clumps) in GC streams, and so may confuse the detection of DM subhaloes. We explore the cases of streams analogous to the observed Palomar 5 and GD1 systems, quantifying the expected incidence of structure caused by GMC perturbations. Deep observations should detect such disturbances regardless of the substructure content of the Milky Way's halo. Detailed modelling will be needed to demonstrate that any detected gaps or clumps were produced by DM subhaloes rather than by molecular clouds.

preprint2016arXiv

Quiescence correlates strongly with directly-measured black hole mass in central galaxies

Roughly half of all stars reside in galaxies without significant ongoing star formation. However, galaxy formation models indicate that it is energetically challenging to suppress the cooling of gas and the formation of stars in galaxies that lie at the centers of their dark matter halos. In this Letter, we show that the dependence of quiescence on black hole and stellar mass is a powerful discriminant between differing models for the mechanisms that suppress star formation. Using observations of 91 star-forming and quiescent central galaxies with directly-measured black hole masses, we find that quiescent galaxies host more massive black holes than star-forming galaxies with similar stellar masses. This observational result is in qualitative agreement with models that assume that effective, more-or-less continuous AGN feedback suppresses star formation, strongly suggesting the importance of the black hole in producing quiescence in central galaxies.

preprint2016arXiv

The Diversity of Growth Histories of Milky Way-mass Galaxies

We use the semi-analytic model developed by Henriques et al. (2015) to explore the origin of star formation history diversity for galaxies that lie at the centre of their dark matter haloes and have present-day stellar masses in the range 5-8 $\times$ 10$^{10}$ M$_{\odot}$, similar to that of the Milky Way. In this model, quenching is the dominant physical mechanism for introducing scatter in the growth histories of these galaxies. We find that present-day quiescent galaxies have a larger variety of growth histories than star-formers since they underwent 'staggered quenching' - a term describing the correlation between the time of quenching and present-day halo mass. While halo mass correlates broadly with quiescence, we find that quiescence is primarily a function of black hole mass, where galaxies quench when heating from their active galactic nuclei becomes sufficient to offset the redshift-dependent cooling rate. In this model, the emergence of a prominent quiescent population is the main process that flattens the stellar mass-halo mass relation at mass scales at or above that of the Milky Way.

preprint2016arXiv

The EAGLE simulations of galaxy formation: public release of halo and galaxy catalogues

We present the public data release of halo and galaxy catalogues extracted from the EAGLE suite of cosmological hydrodynamical simulations of galaxy formation. These simulations were performed with an enhanced version of the GADGET code that includes a modified hydrodynamics solver, time-step limiter and subgrid treatments of baryonic physics, such as stellar mass loss, element-by-element radiative cooling, star formation and feedback from star formation and black hole accretion. The simulation suite includes runs performed in volumes ranging from 25 to 100 comoving megaparsecs per side, with numerical resolution chosen to marginally resolve the Jeans mass of the gas at the star formation threshold. The free parameters of the subgrid models for feedback are calibrated to the redshift z=0 galaxy stellar mass function, galaxy sizes and black hole mass - stellar mass relation. The simulations have been shown to match a wide range of observations for present-day and higher-redshift galaxies. The raw particle data have been used to link galaxies across redshifts by creating merger trees. The indexing of the tree produces a simple way to connect a galaxy at one redshift to its progenitors at higher redshift and to identify its descendants at lower redshift. In this paper we present a relational database which we are making available for general use. A large number of properties of haloes and galaxies and their merger trees are stored in the database, including stellar masses, star formation rates, metallicities, photometric measurements and mock gri images. Complex queries can be created to explore the evolution of more than 10^5 galaxies, examples of which are provided in appendix. (abridged)

preprint2016arXiv

The galaxy correlation function as a constraint on galaxy formation physics

We introduce methods which allow observed galaxy clustering to be used together with observed luminosity or stellar mass functions to constrain the physics of galaxy formation. We show how the projected two-point correlation function of galaxies in a large semi-analytic simulation can be estimated to better than ~10% using only a very small subsample of the subhalo merger trees. This allows measured correlations to be used as constraints in a Monte Carlo Markov Chain exploration of the astrophysical and cosmological parameter space. An important part of our scheme is an analytic profile which captures the simulated satellite distribution extremely well out to several halo virial radii. This is essential to reproduce the correlation properties of the full simulation at intermediate separations. As a first application, we use low-redshift clustering and abundance measurements to constrain a recent version of the Munich semi-analytic model. The preferred values of most parameters are consistent with those found previously, with significantly improved constraints and somewhat shifted "best" values for parameters that primarily affect spatial distributions. Our methods allow multi-epoch data on galaxy clustering and abundance to be used as joint constraints on galaxy formation. This may lead to significant constraints on cosmological parameters even after marginalising over galaxy formation physics.

preprint2016arXiv

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

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

preprint2015arXiv

Surface photometry of BCGs and intracluster stars in Lambda-CDM

We simulate the phase-space distribution of stellar mass in 9 massive Lambda-CDM galaxy clusters by applying the semi-analytic particle tagging method of Cooper et al. to the Phoenix suite of high-resolution N-body simulations (M200 = 7.5 to 33 x 10^14 Msol). The resulting surface brightness (SB) profiles of brightest cluster galaxies (BCGs) match well to observations. On average, stars formed in galaxies accreted by the BCG account for ~90 per cent of its total mass (the remainder is formed in situ). In circular BCG-centred apertures, the superposition of multiple debris clouds (each ~10 per cent of the total BCG mass) from different progenitors can result in an extensive outer diffuse component, qualitatively similar to a 'cD envelope'. These clouds typically originate from tidal stripping at z < 1 and comprise both streams and the extended envelopes of other massive galaxies in the cluster. Stars at very low SB contribute a significant fraction of the total cluster stellar mass budget: in the central 1 Mpc^2 of a z ~ 0.15 cluster imaged at SDSS-like resolution, our fiducial model predicts 80-95 per cent of stellar mass below a SB of mu_V = 26.5 mag arcsec^2 is associated with accreted stars in the envelope of the BCG. The ratio of BCG stellar mass (including this diffuse component) to total cluster stellar mass is ~30 per cent.

preprint2015arXiv

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

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

preprint2015arXiv

The EAGLE simulations of galaxy formation: calibration of subgrid physics and model variations

We present results from thirteen cosmological simulations that explore the parameter space of the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) simulation project. Four of the simulations follow the evolution of a periodic cube L = 50 cMpc on a side, and each employs a different subgrid model of the energetic feedback associated with star formation. The relevant parameters were adjusted so that the simulations each reproduce the observed galaxy stellar mass function at z = 0.1. Three of the simulations fail to form disc galaxies as extended as observed, and we show analytically that this is a consequence of numerical radiative losses that reduce the efficiency of stellar feedback in high-density gas. Such losses are greatly reduced in the fourth simulation - the EAGLE reference model - by injecting more energy in higher density gas. This model produces galaxies with the observed size distribution, and also reproduces many galaxy scaling relations. In the remaining nine simulations, a single parameter or process of the reference model was varied at a time. We find that the properties of galaxies with stellar mass <~ M* (the "knee" of the galaxy stellar mass function) are largely governed by feedback associated with star formation, while those of more massive galaxies are also controlled by feedback from accretion onto their central black holes. Both processes must be efficient in order to reproduce the observed galaxy population. In general, simulations that have been calibrated to reproduce the low-redshift galaxy stellar mass function will still not form realistic galaxies, but the additional requirement that galaxy sizes be acceptable leads to agreement with a large range of observables.

preprint2015arXiv

The unexpected diversity of dwarf galaxy rotation curves

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

preprint2015arXiv

Unifying X-ray Scaling Relations from Galaxies to Clusters

We examine a sample of $\sim 250 000$ "locally brightest galaxies" selected from the Sloan Digital Sky Survey to be central galaxies within their dark matter halos. These were previously stacked by the Planck Collaboration to measure the Sunyaev-Zel'dovich signal as a function of central galaxy stellar mass. Here, we stack the X-ray emission from these halos using data from the ROSAT All-Sky Survey. We detect emission across almost our entire sample, including emission which we attribute to hot gas around galaxies spanning a range of 1.2 dex in stellar mass (corresponding to nearly two orders of magnitude in halo mass) down to $M* = 10^{10.8} M_{\odot}$ ($M_{500} \approx 10^{12.6} M_{\odot}$). Over this range, the X-ray luminosity can be fit by a power-law, either of stellar mass or of halo mass. A single unified scaling relation between mass and $L_X$ applies for galaxies, groups, and clusters. This relation has a steeper slope than expected for self-similarity, in contrast to the $Y_{SZ}$-$M_{500}$ relation, showing the importance of non-gravitational heating. If this heating is predominantly due to AGN feedback, the lack of a break in our relation suggests that AGN feedback is tightly self-regulated and fairly gentle, in agreement with recent high-resolution simulations. Our results are consistent with earlier measurements of the $L_X$-$L_K$ relation for elliptical galaxies and of the $L_X$-$M_{500}$ relation for optically-selected galaxy clusters. However, our $L_X$-$M_{500}$ relation lies more than a factor of two below most previous relations based on X-ray-selected cluster samples. We argue that optical selection offers a less biased view of the $L_X$-$M_{500}$ relation.

preprint2014arXiv

Elemental Abundances in Milky Way-like Galaxies from a Hierarchical Galaxy Formation Model

We develop a new method to account for the finite lifetimes of stars and trace individual abundances within a semi-analytic model of galaxy formation. At variance with previous methods, based on the storage of the (binned) past star formation history of model galaxies, our method projects the information about the metals produced by each simple stellar population (SSP) in the future. Using this approach, an accurate accounting of the timings and properties of the individual SSPs composing model galaxies is possible. We analyse the dependence of our chemical model on various ingredients, and apply it to six simulated haloes of roughly Milky Way mass and with no massive close neighbour at z=0. For all models considered, the [Fe/H] distributions of the stars in the disc component are in good agreement with Milky Way data, while for the spheroid component (whose formation we model only through mergers) these are offset low with respect to observational measurements for the Milky Way bulge. This is a consequence of narrow star formation histories, with relatively low rates of star formation. The slow recycling of gas and energy from supernovae in our chemical model has important consequences on the predicted star formation rates, which are systematically lower than the corresponding rates in the same physical model but with an instantaneous recycling approximation. The halo that resembles most our Galaxy in terms of its global properties also reproduces the observed relation between the average metallicity and luminosity of the Milky Way satellites, albeit with a slightly steeper slope.

preprint2014arXiv

Lensing Model of MACS J1149.5+2223 - I. Cluster mass reconstruction

Measurements of the total logarithmic central slope of the mass profile in galaxy clusters constrain their evolution and assembly history and that of their brightest cluster galaxies. We report the first full surface brightness distribution modelling of the inner region of the galaxy cluster MACS J1149.5+2223. We compare these results with a position-based modelling approach for which we employ more than twice the previously known positional constraints. This is the first time that the detailed lensed image configuration of two non-central cluster galaxies with Einstein rings has been mapped. Due to the extended radial coverage provided by the multiple images in this system, we are able to determine the slope $\partial \log{ κ}/\partial \log{R} = -0.33$ of the total projected mass distribution from $8$ to $80~\mathrm{kpc}$. This is within the cluster-to-cluster scatter estimates from previous cluster measurements. Our reconstruction of the image surface brightness distribution of the large central spiral galaxy has a root mean square residual for all image pixels of $1.14~σ$, where $σ$ is the observational background noise. This corresponds to a reconstruction of the positions of bright clumps in the central galaxy with an rms of $0.063~\mathrm{arcsec}$.

preprint2014arXiv

Local Group galaxies emerge from the dark

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

preprint2014arXiv

Matched filter optimization of kSZ measurements with a reconstructed cosmological flow field

We develop and test a new statistical method to measure the kinematic Sunyaev-Zel'dovich (kSZ) effect. A sample of independently detected clusters is combined with the cosmic flow field predicted from a galaxy redshift survey in order to derive a matched filter that optimally weights the kSZ signal for the sample as a whole given the noise involved in the problem. We apply this formalism to realistic mock microwave skies based on cosmological $N$-body simulations, and demonstrate its robustness and performance. In particular, we carefully assess the various sources of uncertainty, cosmic microwave background primary fluctuations, instrumental noise, uncertainties in the determination of the velocity field, and effects introduced by miscentring of clusters and by uncertainties of the mass-observable relation (normalization and scatter). We show that available data (\plk\ maps and the MaxBCG catalogue) should deliver a $7.7σ$ detection of the kSZ. A similar cluster catalogue with broader sky coverage should increase the detection significance to $\sim 13σ$. We point out that such measurements could be binned in order to study the properties of the cosmic gas and velocity fields, or combined into a single measurement to constrain cosmological parameters or deviations of the law of gravity from General Relativity.

preprint2014arXiv

Satellite abundances around bright isolated galaxies II: radial distribution and environmental effects

We use the SDSS/DR8 galaxy sample to study the radial distribution of satellite galaxies around isolated primaries, comparing to semi-analytic models of galaxy formation based on the Millennium and Millennium-II simulations. SDSS satellites behave differently around high- and low-mass primaries: those orbiting objects with $M_*>10^{11}M_\odot$ are mostly red and are less concentrated towards their host than the inferred dark matter halo, an effect that is very pronounced for the few blue satellites. On the other hand, less massive primaries have steeper satellite profiles that agree quite well with the expected dark matter distribution and are dominated by blue satellites, even in the inner regions where strong environmental effects are expected. In fact, such effects appear to be strong only for primaries with $M_* > 10^{11}M_\odot$. This behaviour is not reproduced by current semi-analytic simulations, where satellite profiles always parallel those of the dark matter and satellite populations are predominantly red for primaries of all masses. The disagreement with SDSS suggests that environmental effects are too efficient in the models. Modifying the treatment of environmental and star formation processes can substantially increase the fraction of blue satellites, but their radial distribution remains significantly shallower than observed. It seems that most satellites of low-mass primaries can continue to form stars even after orbiting within their joint halo for 5 Gyr or more.

preprint2014arXiv

The diverse formation histories of simulated disc galaxies

We analyze the formation histories of 19 galaxies from cosmological smoothed particle hydrodynamics zoom-in resimulations. We construct mock three-colour images and show that the models reproduce observed trends in the evolution of galaxy colours and morphologies. However, only a small fraction of galaxies contains bars. Many galaxies go through phases of central mass growth by in-situ star formation driven by gas-rich mergers or misaligned gas infall. These events lead to accretion of low-angular momentum gas to the centres and leave imprints on the distributions of z=0 stellar circularities, radii and metallicities as functions of age. Observations of the evolution of structural properties of samples of disc galaxies at z=2.5-0.0 infer continuous mass assembly at all radii. Our simulations can only explain this if there is a significant contribution from mergers or misaligned infall, as expected in a LambdaCDM universe. Quiescent merger histories lead to high kinematic disc fractions and inside-out growth, but show little central growth after the last `destructive' merger at z>1.5. For sufficiently strong feedback, as assumed in our models, a moderate amount of merging does not seem to be a problem for the z=0 disc galaxy population, but may rather be a requirement. The average profiles of simulated disc galaxies agree with observations at z>=1.5. At z<=1, there is too much growth in size and too little growth in central mass, possibly due to the under-abundance of bars. The discrepancies may partly be caused by differences between the star formation histories of the simulations and those assumed for observations.

preprint2014arXiv

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

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

preprint2014arXiv

The re-distribution of matter in the cores of galaxy clusters

We present cosmological N-body resimulations of the assembly of the Brightest Cluster Galaxies (BCGs) in rich clusters. At $z=2$ we populate dark matter subhalos with self-gravitating stellar systems whose abundance and structure match observed high-redshift galaxies. By $z=0$, mergers have built much larger galaxies at cluster centre. Their dark matter density profiles are shallower than in corresponding dark-matter-only simulations, but their total mass density profiles (stars + dark matter) are quite similar. Differences are found only at radii where the effects of central black holes may be significant. Dark matter density slopes shallower than $γ=1.0$ occur for $r/r_{200} < 0.015$, close to the half-light radii of the BCGs. Our experiments support earlier suggestions that NFW-like profiles are an attractor for the hierarchical growth of structure in collisionless systems -- total mass density profiles asymptote to the solution found in dark-matter-only simulations over the radial range where mergers produce significant mixing between stars and dark matter. Simulated dark matter fractions are substantially higher in BCGs than in field ellipticals, reaching 80\% within the half-light radius. We also estimate that supermassive black hole mergers should create BCG cores as large as $r_{c}\sim 3 \,\mathrm{kpc}$. The good agreement of all these properties with recent observational studies of BCG structure suggests that dissipational processes have not played a dominant role in the assembly of the observed systems.

preprint2014arXiv

The Tenth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-III Apache Point Observatory Galactic Evolution Experiment

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

preprint2014arXiv

Why stellar feedback promotes disc formation in simulated galaxies

We study how feedback influences baryon infall onto galaxies using cosmological, zoom-in simulations of haloes with present mass $M_{vir}=6.9\times10^{11} M_{\odot}$ to $1.7\times10^{12} M_{\odot}$. Starting at z=4 from identical initial conditions, implementations of weak and strong stellar feedback produce bulge- and disc-dominated galaxies, respectively. Strong feedback favours disc formation: (1) because conversion of gas into stars is suppressed at early times, as required by abundance matching arguments, resulting in flat star formation histories and higher gas fractions; (2) because 50% of the stars form in situ from recycled disc gas with angular momentum only weakly related to that of the z=0 dark halo; (3) because late-time gas accretion is typically an order of magnitude stronger and has higher specific angular momentum, with recycled gas dominating over primordial infall; (4) because 25-30% of the total accreted gas is ejected entirely before z~1, removing primarily low angular momentum material which enriches the nearby inter-galactic medium. Most recycled gas roughly conserves its angular momentum, but material ejected for long times and to large radii can gain significant angular momentum before re-accretion. These processes lower galaxy formation efficiency in addition to promoting disc formation.

preprint2013arXiv

Estimating gas masses and dust-to-gas ratios from optical spectroscopy

We present a method to estimate the total gas column density, dust-to-gas and dust-to-metal ratios of distant galaxies from rest-frame optical spectra. The technique exploits the sensitivity of certain optical lines to changes in depletion of metals onto dust grains and uses photo-ionization models to constrain these physical ratios along with the metallicity (Z) and dust column density. We compare our gas column density estimates with HI and CO gas mass estimates in nearby galaxies to show that we recover their total gas mass surface density to within a factor of 2 up to a total surface gas mass density of ~75 Mo/pc^2. Our technique is independent of the conversion factor of CO to H2 and we show that a Z-dependent XCO is required to achieve good agreement between our measurements and that provided by CO and HI. However we also show that our method can not be reliably aperture corrected to total gas mass. We calculate dust-to-gas ratios for all star-forming galaxies in the SDSS DR7 and show the resulting dependence on Z agrees well with the trend inferred from modelling of the dust emission of nearby galaxies using FIR data. We also present estimates of the variation of the dust-to-metal ratio with Z and show that this is poorly constrained at Z<~0.5Z_sun. We end with a study of the inventory of gas in the central regions, defined both in terms of a fixed physical radius and as a fixed fraction of the half-light radius, of ~70,000 star-forming galaxies from the Sloan Digital Sky Survey. We show that their central gas content and gas depletion times are not accurately predicted by a single parameter, but in agreement with recent studies we find that a combination of the M* and some measure of central concentration provides a good predictor of gas content in galaxies. We also identify a population of galaxies with low surface densities of stars and very long gas depletion times.

preprint2013arXiv

Exploring the non-linear density field in the Millennium simulations with tessellations - I. The probability distribution function

We use the Delaunay Tessellation Field Estimator (DTFE) to study the one-point density distribution functions of the Millennium (MS) and Millennium-II (MS-II) simulations. The DTFE technique is based directly on the particle positions, without requiring any type of smoothing or analysis grid, thereby providing high sensitivity to all non-linear structures resolved by the simulations. In order to identify the detailed origin of the shape of the one-point density probability distribution function (PDF), we decompose the simulation particles according to the mass of their host FoF halos, and examine the contributions of different halo mass ranges to the global density PDF. We model the one-point distribution of the FoF halos in each halo mass bin with a set of Monte Carlo realizations of idealized NFW dark matter halos, finding that this reproduces the measurements from the N-body simulations reasonably well, except for a small excess present in simulation results. This excess increases with increasing halo mass. We show that its origin lies in substructure, which becomes progressively more abundant and better resolved in more massive dark matter halos. We demonstrate that the high density tail of the one-point distribution function in less massive halos is severely affected by the gravitational softening length and the mass resolution. In particular, we find these two parameters to be more important for an accurate measurement of the density PDF than the simulated volume. Combining our results from individual halo mass bins we find that the part of the one-point density PDF originating from collapsed halos can nevertheless be quite well described by a simple superposition of a set of NFW halos with the expected cosmological abundance over the resolved mass range. The transition region to the low-density unbound material is however not well captured by such an analytic halo model.

preprint2013arXiv

Galactic accretion and the outer structure of galaxies in the CDM model

We have combined the semi-analytic galaxy formation model of Guo et al. (2011) with the particle-tagging technique of Cooper et al. (2010) to predict galaxy surface brightness profiles in a representative sample of ~1900 massive dark matter haloes (10^12--10^14 M_sol) from the Millennium II Lambda-CDM N-body simulation. Here we present our method and basic results focusing on the outer regions of galaxies, consisting of stars accreted in mergers. These simulations cover scales from the stellar haloes of Milky Way-like galaxies to the 'cD envelopes' of groups and clusters, and resolve low surface brightness substructure such as tidal streams. We find that the surface density of accreted stellar mass around the central galaxies of dark matter haloes is well described by a Sersic profile, the radial scale and amplitude of which vary systematically with halo mass (M_200). The total stellar mass surface density profile breaks at the radius where accreted stars start to dominate over stars formed in the galaxy itself. This break disappears with increasing M_200 because accreted stars contribute more of the total mass of galaxies, and is less distinct when the same galaxies are averaged in bins of stellar mass, because of scatter in the relation between M_star and M_200. To test our model we have derived average stellar mass surface density profiles for massive galaxies at z~0.08 by stacking SDSS images. Our model agrees well with these stacked profiles and with other data from the literature, and makes predictions that can be more rigorously tested by future surveys that extend the analysis of the outer structure of galaxies to fainter isophotes. We conclude that it is likely that the outer structure of the spheroidal components of galaxies is largely determined by collisionless merging during their hierarchical assembly

preprint2013arXiv

Galaxy formation on the largest scales: The impact of astrophysics on the BAO peak

We investigate the effects of galaxy formation on the baryonic acoustic oscillations (BAO) peak by applying semi-analytic modelling techniques to the Millennium-XXL, a $3 \times 10^{11}$ particle N-body simulation of similar volume to the future EUCLID survey. Our approach explicitly incorporates the effects of tidal fields and stochasticity on halo formation, as well as the presence of velocity bias, spatially correlated merger histories, and the connection of all these with the observable and physical properties of galaxies. We measure significant deviations in the shape of the BAO peak from the expectations of a linear bias model built on top of the nonlinear dark matter distribution. We find that the galaxy correlation function shows an excess close to the maximum of the BAO peak ($r\sim110 Mpc/h$) and a deficit at $r\sim90 Mpc/h$. Depending on the redshift, selection criteria and number density of the galaxy samples, these bias distortions can be up to 5% in amplitude. They are, however, largely absorbed by marginalization over nuisance parameters in current analytical modelling of the BAO peak in configuration space, in particular into the parameter that controls the broadening due to nonlinear evolution. As a result, the galaxy formation effects detected here are unlikely to bias the high-precision measurements planned by the upcoming generation of wide-field galaxy surveys.

preprint2013arXiv

Modelling Element Abundances in Semi-analytic Models of Galaxy Formation

We update the treatment of chemical evolution in the Munich semi-analytic model, L-GALAXIES. Our new implementation includes delayed enrichment from stellar winds, supernovae type II (SNe-II) and supernovae type Ia (SNe-Ia), as well as metallicity-dependent yields and a reformulation of the associated supernova feedback. Two different sets of SN-II yields and three different SN-Ia delay-time distributions (DTDs) are considered, and eleven heavy elements (including O, Mg and Fe) are self-consistently tracked. We compare the results of this new implementation with data on a) local, star-forming galaxies, b) Milky Way disc G dwarfs, and c) local, elliptical galaxies. We find that the z=0 gas-phase mass-metallicity relation is very well reproduced for all forms of DTD considered, as is the [Fe/H] distribution in the Milky Way disc. The [O/Fe] distribution in the Milky Way disc is best reproduced when using a DTD with less than or equal to 50 per cent of SNe-Ia exploding within ~400 Myrs. Positive slopes in the mass-[alpha/Fe] relations of local ellipticals are also obtained when using a DTD with such a minor `prompt' component. Alternatively, metal-rich winds that drive light alpha elements directly out into the circumgalactic medium also produce positive slopes for all forms of DTD and SN-II yields considered. Overall, we find that the best model for matching the wide range of observational data considered here should include a power-law SN-Ia DTD, SN-II yields that take account of prior mass loss through stellar winds, and some direct ejection of light alpha elements out of galaxies.

preprint2013arXiv

Streams in the Aquarius stellar haloes

We use the very high resolution, fully cosmological simulations from the Aquarius project, coupled to a semi-analytical model of galaxy formation, to study the phase-space distribution of halo stars in "solar neighbourhood"-like volumes. We find that this distribution is very rich in substructure in the form of stellar streams for all five stellar haloes we have analysed. These streams can be easily identified in velocity space, as well as in spaces of pseudo-conserved quantities such as E vs. Lz. In our best-resolved local volumes, the number of identified streams ranges from ~ 300 to 600, in very good agreement with previous analytical predictions, even in the presence of chaotic mixing. The fraction of particles linked to (massive) stellar streams in these volumes can be as large as 84%. The number of identified streams is found to decrease as a power-law with galactocentric radius. We show that the strongest limitation to the quantification of substructure in our poorest-resolved local volumes is particle resolution rather than strong diffusion due to chaotic mixing.

preprint2013arXiv

The Growth in Size and Mass of Cluster Galaxies since z=2

We study the formation and evolution of Brightest Cluster Galaxies starting from a $z=2$ population of quiescent ellipticals and following them to $z=0$. To this end, we use a suite of nine high-resolution dark matter-only simulations of galaxy clusters in a $Λ$CDM universe. We develop a scheme in which simulation particles are weighted to generate realistic and dynamically stable stellar density profiles at $z=2$. Our initial conditions assign a stellar mass to every identified dark halo as expected from abundance matching; assuming there exists a one-to-one relation between the visible properties of galaxies and their host haloes. We set the sizes of the luminous components according to the observed relations for $z\sim2$ massive quiescent galaxies. We study the evolution of the mass-size relation, the fate of satellite galaxies and the mass aggregation of the cluster central. From $z=2$, these galaxies grow on average in size by a factor 5 to 10 of and in mass by 2 to 3. The stellar mass growth rate of the simulated BCGs in our sample is of 1.9 in the range $0.3<z<1.0$ consistent with observations, and of 1.5 in the range $0.0<z<0.3$. Furthermore the satellite galaxies evolve to the present day mass-size relation by $z=0$. Assuming passively evolving stellar populations, we present surface brightness profiles for our cluster centrals which resemble those observed for the cDs in similar mass clusters both at $z=0$ and at $z=1$. This demonstrates that the $Λ$CDM cosmology does indeed predict minor and major mergers to occur in galaxy clusters with the frequency and mass ratio distribution required to explain the observed growth in size of passive galaxies since $z=2$. Our experiment shows that Brightest Cluster Galaxies can form through dissipationless mergers of quiescent massive $z=2$ galaxies, without substantial additional star formation.

preprint2013arXiv

The Mass Profile and Accretion History of Cold Dark Matter Halos

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.

preprint2013arXiv

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

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.

preprint2013arXiv

The satellites of the Milky Way - Insights from semi-analytic modelling in a LambdaCDM cosmology

We combine the six high-resolution Aquarius dark matter simulations with a semi-analytic galaxy formation model to investigate the properties of the satellites of Milky Way-like galaxies. We find good correspondence with the observed luminosity function, luminosity-metallicity relation and radial distribution of the Milky Way satellites. The star formation histories of the dwarf galaxies in our model vary widely, in accordance with what is seen observationally. Ram-pressure stripping of hot gas from the satellites leaves a clear imprint of the environment on the characteristics of a dwarf galaxy. We find that the fraction of satellites dominated by old populations of stars matches observations well. However, the internal metallicity distributions of the model satellites appear to be narrower than observed. This may indicate limitations in our treatment of chemical enrichment, which is based on the instantaneous recycling approximation. Our model works best if the dark matter halo of the Milky Way has a mass of ~8 x 10^11 Msun, in agreement with the lower estimates from observations. The galaxy that resembles the Milky Way the most also has the best matching satellite luminosity function, although it does not contain an object as bright as the SMC or LMC. Compared to other semi-analytic models and abundance matching relations we find that central galaxies reside in less massive haloes, but the halo mass-stellar mass relation in our model is consistent both with hydrodynamical simulations and with recent observations.

preprint2012arXiv

Dark matter and cosmic structure

We review the current standard model for the evolution of cosmic structure, tracing its development over the last forty years and focusing specifically on the role played by numerical simulations and on aspects related to the nature of dark matter.

preprint2012arXiv

Evidence for extended gamma-ray emission from galaxy clusters

We report evidence for extended gamma-ray emission from the Virgo, Fornax and Coma clusters based on a maximum-likelihood analysis of the 3-year Fermi-LAT data. For all three clusters, excess emission is observed within three degrees of the center, peaking at the GeV scale. This emission cannot be accounted for by known Fermi sources or by the galactic and extragalactic backgrounds. If interpreted as annihilation emission from supersymmetric dark matter (DM) particles, the data prefer models with a particle mass in the range 20-60 GeV annihilating into the b-bbar channel, or 2-10 GeV and >1 TeV annihilating into mu-mu final states. Our results are consistent with those obtained by Hooper and Linden from a recent analysis of Fermi-LAT data in the region of the Galactic Centre. An extended DM annihilation profile dominated by emission from substructures is preferred over a simple point source model. The significance of DM detection is 4.4 sigma in Virgo and lower in the other two clusters. We also consider the possibility that the excess emission arises from cosmic ray (CR) induced gamma-rays, and infer a CR level within a factor of three of that expected from analytical models. However, the significance of a CR component is lower than the significance of a DM component, and there is no need for such a CR component in the presence of a DM component in the preferred DM mass range. We also set flux and cross-section upper limits for DM annihilation into the b-bbar and mu-mu channels in all three clusters.

preprint2012arXiv

Galactic star formation and accretion histories from matching galaxies to dark matter haloes

We present a new statistical method to determine the relationship between the stellar masses of galaxies and the masses of their host dark matter haloes over the entire cosmic history from z~4 to the present. This multi-epoch abundance matching (MEAM) model self-consistently takes into account that satellite galaxies first become satellites at times earlier than they are observed. We employ a redshift-dependent parameterization of the stellar-to-halo mass relation to populate haloes and subhaloes in the Millennium simulations with galaxies, requiring that the observed stellar mass functions at several redshifts be reproduced simultaneously. Using merger trees extracted from the dark matter simulations in combination with MEAM, we predict the average assembly histories of galaxies, separating into star formation within the galaxies (in-situ) and accretion of stars (ex-situ). The peak star formation efficiency decreases with redshift from 23% at z=0 to 9% at z=4 while the corresponding halo mass increases from 10^11.8M\odot to 10^12.5M\odot. The star formation rate of central galaxies peaks at a redshift which depends on halo mass; for massive haloes this peak is at early cosmic times while for low-mass galaxies the peak has not been reached yet. In haloes similar to that of the Milky-Way about half of the central stellar mass is assembled after z=0.7. In low-mass haloes, the accretion of satellites contributes little to the assembly of their central galaxies, while in massive haloes more than half of the central stellar mass is formed ex-situ with significant accretion of satellites at z<2. We find that our method implies a cosmic star formation history and an evolution of specific star formation rates which are consistent with those inferred directly. We present convenient fitting functions for stellar masses, star formation rates, and accretion rates as functions of halo mass and redshift.

preprint2012arXiv

Idealized models for galactic disk formation and evolution in 'realistic' LambdaCDM haloes

We study the dynamics of galactic disk formation and evolution in 'realistic' LambdaCDM haloes with idealized baryonic initial conditions. We add rotating spheres of hot gas at z=1.3 to two fully cosmological dark-matter-only halo (re)simulations. The gas cools according to an artificial and adjustable cooling function to form a rotationally supported galaxy. The simulations evolve in the full cosmological context until z=0. We vary the angular momentum and density profiles of the initial gas sphere, the cooling time and the orientation of the angular momentum vector to study the effects on the evolution of the disk. The final disks show realistic structural and kinematic properties. The slower the cooling/accretion processes, the higher the kinematic disk-to-bulge ratio D/B of the resulting system. We find that the initial orientation of the gas angular momentum with respect to the halo has a major effect on the resulting D/B. The most stable systems result from orientations parallel to the halo minor axis, but the sign of the spin can have a strong effect. Despite the spherical and coherently rotating initial gas distribution, the orientation of the central disk and of the outer gas components and the relative angle between the components can all change by more than 90 degrees over several billion years. Disks can form from initial conditions oriented parallel to the major axis, but there is often strong misalignment between inner and outer material. The more the orientation of the baryonic angular momentum changes during the evolution, the lower the final D/B. The behaviour varies strongly from halo to halo. Even our simple initial conditions can lead to strong bars, dominant bulges, massive, misaligned rings and counter-rotating components. We discuss how our results may relate to the failure or success of fully cosmological disk formation simulations. (abridged)

preprint2012arXiv

On the Tidal Dependence of Galaxy Properties

Using volume-limited samples drawn from The Sloan Digital Sky Survey Data Release 7 (SDSS DR7), we measure the tidal environment of galaxies, which we characterize by the ellipticity e of the potential field calculated from the smoothed spatial number density 1 + δ of galaxies. We analyze if galaxy properties, including color, Dn4000, concentration and size correlate with e, in addition to depending on 1 + δ. We find that there exists a transition smoothing scale at which correlations/anti-correlations with e reverse. This transition scale is well represented by the distance to the 3rd- nearest-neighbor of a galaxy in a volume limited sample with Mr < -20 which has a distribution peaked at ~ 2 h-1Mpc. We further demonstrate that this scale corresponds to that where the correlation between the color of galaxies and environmental density 1+δ is the strongest. For this optimal smoothing R0 no additional correlations with e are observed. The apparent dependence on tidal ellipticity e at other smoothing scales Rs can be viewed as a geometric effect, arising from the cross correlation be- tween (1+δo) and e(Rs). We perform the same analysis on numerical simulations with semi-analytical modeling (SAM) of galaxy formation. The e dependence of the galaxy properties shows similar behavior to that in the SDSS, although the color-density cor- relation is significantly stronger in the SAM. The 'optimal adaptive smoothing scale' in the SAM is also closely related to the distance to the 3rd-nearest-neighbor of a galaxy, and its characteristic value is consistent with, albeit slightly smaller than for SDSS.

preprint2012arXiv

Satellite abundances around bright isolated galaxies

We study satellite galaxy abundances in SDSS by counting photometric galaxies around isolated bright primaries. We present results as a function of the luminosity, stellar mass and colour of the satellites, and of the stellar mass and colour of the primaries. For massive primaries the luminosity and stellar mass functions of satellites are similar in shape to those of field galaxies, but for lower mass primaries they are significantly steeper. The steepening is particularly marked for the stellar mass function. Satellite abundance increases strongly with primary stellar mass, approximately in proportion to expected dark halo mass. Massive red primaries have up to a factor of 2 more satellites than blue ones of the same stellar mass. Satellite galaxies are systematically redder than field galaxies of the same stellar mass. Satellites are also systematically redder around more massive primaries. At fixed primary mass, they are redder around red primaries. We select similarly isolated galaxies from mock catalogues based on the simulations of Guo et al.(2011) and analyze them in parallel with the SDSS data. The simulation reproduces all the above trends qualitatively, except for the steepening of the satellite luminosity and stellar mass functions. Model satellites, however, are systematically redder than in the SDSS, particularly at low mass and around low-mass primaries. Simulated haloes of a given mass have satellite abundances that are independent of central galaxy colour, but red centrals tend to have lower stellar masses, reflecting earlier quenching of their star formation by feedback. This explains the correlation between satellite abundance and primary colour in the simulation. The correlation between satellite colour and primary colour arises because red centrals live in haloes which are more massive, older and more gas-rich, so that satellite quenching is more efficient.

preprint2012arXiv

Satellites and Haloes of Dwarf Galaxies

We study the abundance of satellite galaxies as a function of primary stellar mass using the SDSS/DR7 spectroscopic catalogue. In contrast with previous studies, which focussed mainly on bright primaries, our central galaxies span a wide range of stellar mass, 10^7.5 < M_*^pri/M_sun < 10^11, from dwarfs to central cluster galaxies. Our analysis confirms that the average number of satellites around bright primaries, when expressed in terms of satellite-to-primary stellar mass ratio (m_*^sat/M_*^pri), is a strong function of M_*^pri. On the other hand, satellite abundance is largely independent of primary mass for dwarf primaries (M_*^pri<10^10 M_sun). These results are consistent with galaxy formation models in the LCDM scenario. We find excellent agreement between SDSS data and semi-analytic mock galaxy catalogues constructed from the Millennium-II Simulation. Satellite galaxies trace dark matter substructure in LCDM, so satellite abundance reflects the dependence on halo mass, M_200, of both substructure and galaxy stellar mass (M_*). Since dark matter substructure is almost scale-free, the dependence of satellite abundance on primary mass results solely from the well-defined characteristic mass in the galaxy mass-halo mass relation. On dwarf galaxy scales, where models predict a power-law scaling, M_* \propto M_200^2.5, similarity is preserved and satellite abundance is independent of primary mass. For primaries brighter than the characteristic mass of the M_*-M_200 relation, satellite abundance increases strongly with primary mass. Our results provide strong support for the steep, approximately power-law dependence of dwarf galaxy mass on halo mass envisioned in LCDM galaxy formation models.

preprint2012arXiv

Shallow Dark Matter Cusps in Galaxy Clusters

We study the evolution of the stellar and dark matter components in a galaxy cluster of $10^{15} \, \rm{M_{\odot}}$ from $z=3$ to the present epoch using the high-resolution collisionless simulations of Ruszkowski & Springel (2009). At $z=3$ the dominant progenitor halos were populated with spherical model galaxies with and without accounting for adiabatic contraction. We apply a weighting scheme which allows us to change the relative amount of dark and stellar material assigned to each simulation particle in order to produce luminous properties which agree better with abundance matching arguments and observed bulge sizes at $z=3$. This permits the study of the effect of initial compactness on the evolution of the mass-size relation. We find that for more compact initial stellar distributions the size of the final Brightest Cluster Galaxy grows with mass according to $r\propto M^{2}$, whereas for more extended initial distributions, $r\propto M$. Our results show that collisionless mergers in a cosmological context can reduce the strength of inner dark matter cusps with changes in logarithmic slope of 0.3 to 0.5 at fixed radius. Shallow cusps such as those found recently in several strong lensing clusters thus do not necessarily conflict with CDM, but may rather reflect on the initial structure of the progenitor galaxies, which was shaped at high redshift by their formation process.

preprint2012arXiv

Stellar Disks in Aquarius Dark Matter Haloes

We investigate the gravitational interactions between live stellar disks and their dark matter halos, using LCDM haloes similar in mass to that of the Milky Way taken from the Aquarius Project. We introduce the stellar disks by first allowing the haloes to respond to the influence of a growing rigid disk potential from z = 1.3 to z = 1.0. The rigid potential is then replaced with star particles which evolve self-consistently with the dark matter particles until z = 0.0. Regardless of the initial orientation of the disk, the inner parts of the haloes contract and change from prolate to oblate as the disk grows to its full size. When the disk normal is initially aligned with the major axis of the halo at z=1.3, the length of the major axis contracts and becomes the minor axis by z=1.0. Six out of the eight disks in our main set of simulations form bars, and five of the six bars experience a buckling instability that results in a sudden jump in the vertical stellar velocity dispersion and an accompanying drop in the m=2 Fourier amplitude of the disk surface density. The bars are not destroyed by the buckling but continue to grow until the present day. Bars are largely absent when the disk mass is reduced by a factor of two or more; the relative disk-to-halo mass is therefore a primary factor in bar formation and evolution. A subset of the disks is warped at the outskirts and contains prominent non-coplanar material with a ring-like structure. Many disks reorient by large angles between z=1 and z=0, following a coherent reorientation of their inner haloes. Larger reorientations produce more strongly warped disks, suggesting a tight link between the two phenomena. The origins of bars and warps appear independent: some disks with strong bars show no disturbances at the outskirts, while the disks with the weakest bars show severe warps.

preprint2012arXiv

The Effect of Particle Noise in N-body Simulations of Gravitational Lensing

High resolution dark matter only simulations provide a realistic and fully general means to study the theoretical predictions of cosmological structure formation models for gravitational lensing. Due to the finite number of particles, the density field only becomes smooth on scales beyond a few times the local mean interparticle separation. This introduces noise on the gravitational lensing properties such as the surface mass density, the deflection angles and the magnification. At some small-scale mass limit, the noise due to the discreteness of the N-body simulation becomes comparable to the effects of physical substructures. We present analytic expressions to quantify the Poisson noise and study its scaling with the particle number of the simulation and the Lagrangian smoothing size. We use the Phoenix set of simulations, currently the largest available dark matter simulations of clusters to study the effect of limited numerical resolution and the gravitational strong lensing effects of substructure. We quantify the smallest resolved substructure, in the sense that the effect of the substructure on any strong lensing property is significant compared to the noise, and we find that the result is roughly independent of the strong lensing property. A simple scaling relates the smallest resolved substructures in a simulation with the resolution of the N-body simulation.

preprint2012arXiv

The effects of alignment and ellipticity on the clustering of galaxies

We investigate the effects of halo ellipticity and alignment with larger-scale structure on the galaxy correlation function. We base our analysis on the galaxy formation models of Guo et al. (2011), run on the Millennium Simulations. We quantify the importance of these properties of the galaxy distribution by randomizing the angular positions of satellite galaxies within haloes, either coherently or individually, while keeping the distance to their respective central galaxies fixed. We find that the effect of disrupting the alignment with larger-scale structure is a ~2 per cent decrease in the galaxy correlation function around r=1.8 Mpc/h. Sphericalizing the ellipsoidal distributions of galaxies within haloes decreases the correlation function by up to 20 per cent for r<1 Mpc/h. Similar results apply to power spectra and redshift-space correlation functions. Models such as those based on the Halo Occupation Distribution, which adopt a spherically averaged profile for the galaxy distributions within haloes, will therefore significantly underestimate the clustering on sub-Mpc scales.

preprint2012arXiv

The effects of halo alignment and shape on the clustering of galaxies

We investigate the effects of halo shape and its alignment with larger scale structure on the galaxy correlation function. We base our analysis on the galaxy formation models of Guo et al., run on the Millennium Simulations. We quantify the importance of these effects by randomizing the angular positions of satellite galaxies within haloes, either coherently or individually, while keeping the distance to their respective central galaxies fixed. We find that the effect of disrupting the alignment with larger scale structure is a ~2 per cent decrease in the galaxy correlation function around r=1.8 Mpc/h. We find that sphericalizing the ellipsoidal distributions of galaxies within haloes decreases the correlation function by up to 20 per cent for r<1 Mpc/h and increases it slightly at somewhat larger radii. Similar results apply to power spectra and redshift-space correlation functions. Models based on the Halo Occupation Distribution, which place galaxies spherically within haloes according to a mean radial profile, will therefore significantly underestimate the clustering on sub-Mpc scales. In addition, we find that halo assembly bias, in particular the dependence of clustering on halo shape, propagates to the clustering of galaxies. We predict that this aspect of assembly bias should be observable through the use of extensive group catalogues.

preprint2012arXiv

The Origin of Disks and Spheroids in Simulated Galaxies

In the simplest scenario, disk galaxies form predominantly in halos with high angular momentum and quiet recent assembly history, whereas spheroids are the slowly-rotating remnants of repeated merging events. We explore these assumptions using one hundred systems with halo masses similar to that of the Milky Way, identified in a series of cosmological gasdynamical simulations GIMIC. At z=0, the simulated galaxies exhibit a wide variety of morphologies, from dispersion-dominated spheroids to pure disk galaxies. Surprisingly, these morphological features are very poorly correlated with their halo properties: disks form in halos with high and low net spin, and mergers play a negligible role in the formation of spheroid stars, most of which form in-situ. More important to morphology is the coherent alignment of the angular momentum of baryons that accrete over time to form a galaxy. Spheroids tend to form when the spin of newly-accreted gas is misaligned with that of the extant galaxy, leading to the episodic formation of stars with different kinematics that cancel out the net rotation of the system. Disks, on the other hand, form out of gas that flows in with similar angular momentum to that of earlier-accreted material. Gas accretion from a hot corona thus favours disk formation, whereas gas that flows "cold", often along separate, misaligned filaments, favours the formation of spheroids. In this scenario, most spheroids consist of superpositions of stellar components with distinct kinematics, age, and metallicity, an arrangement that might survive to the present day given the paucity of major mergers. Since angular momentum is acquired largely at turnaround, morphology is imprinted early by the interplay of the tidal field and the shape of the material destined to form the galaxy.

preprint2011arXiv

Autocorrelations of stellar light and mass at z~0 and ~1: From SDSS to DEEP2

We present measurements of projected autocorrelation functions w_p(r_p) for the stellar mass of galaxies and for their light in the U, B and V bands, using data from the third data release of the DEEP2 Galaxy Redshift Survey and the final data release of the Sloan Digital Sky Survey (SDSS). We investigate the clustering bias of stellar mass and light by comparing these to projected autocorrelations of dark matter estimated from the Millennium Simulations (MS) at z=1 and 0.07, the median redshifts of our galaxy samples. All of the autocorrelation and bias functions show systematic trends with spatial scale and waveband which are impressively similar at the two redshifts. This shows that the well-established environmental dependence of stellar populations in the local Universe is already in place at z=1. The recent MS-based galaxy formation simulation of Guo et al. (2011) reproduces the scale-dependent clustering of luminosity to an accuracy better than 30% in all bands and at both redshifts, but substantially overpredicts mass autocorrelations at separations below about 2 Mpc. Further comparison of the shapes of our stellar mass bias functions with those predicted by the model suggests that both the SDSS and DEEP2 data prefer a fluctuation amplitude of sigma_8 0.8 rather than the sigma_8=0.9 assumed by the MS.

preprint2011arXiv

Chemical signatures of formation processes in the stellar populations of simulated galaxies

We study the chemical properties of the stellar populations in eight simulations of the formation of Milky-Way mass galaxies in a LCDM Universe. Our simulations include metal-dependent cooling and an explicitly multiphase treatment of the effects on the gas of cooling, enrichment and supernova feedback. We search for correlations between formation history and chemical abundance patterns. Differing contributions to spheroids and discs from in situ star formation and from accreted populations are reflected in differing chemical properties. Discs have younger stellar populations, with most stars forming in situ and with low alpha-enhancement from gas which never participated in a galactic outflow. Up to 15 per cent of disc stars can come from accreted satellites. These tend to be alpha-enhanced, older and to have larger velocity dispersions than the in situ population. Inner spheroids have old, metal-rich and alpha-enhanced stars which formed primarily in situ, more than 40 per cent from material recycled through earlier galactic winds. Few accreted stars are found in the inner spheroid unless a major merger occurred recently. Such stars are older, more metal-poor and more alpha-enhanced than the in situ population. Stellar haloes tend to have low metallicity and high alpha-enhancement. The outer haloes are made primarily of accreted stars. Their mean metallicity and alpha-enhancement reflect the masses of the disrupted satellites where they formed: more massive satellites typically have higher [Fe/H] and lower [alpha/Fe]. Surviving satellites have distinctive chemical patterns which reflect their extended, bursty star formation histories. These produce lower alpha-enhancement at given metallicity than in the main galaxy, in agreement with observed trends in the Milky Way.

preprint2011arXiv

Clues to the "Magellanic Galaxy" from Cosmological Simulations

We use cosmological simulations from the Aquarius Project to study the orbital history of the Large Magellanic Cloud (LMC) and its potential association with other satellites of the Milky Way (MW). We search for dynamical analogs to the LMC and find a subhalo that matches the LMC position and velocity at either of its two most recent pericentric passages. This suggests that the LMC is not necessarily on its first approach to the MW, provided that the virial mass of the Milky Way is M_200 ~1.8e12 M_sun. The simulation results yield specific predictions for the position and velocity of systems associated with the LMC prior to infall. If on first approach, most should lie close to the LMC because the Galactic tidal field has not yet had enough time to disperse them. If on second approach, the list of potential associates increases substantially. Interestingly, our analysis rules out an LMC association for Draco and Ursa Minor, two of the dwarf spheroidals suggested by Lynden-Bell & Lynden-Bell to form part of the "Magellanic Ghostly Stream". Our results also indicate that the direction of the orbital angular momentum is a powerful test of LMC association. This test, however, requires precise proper motions, which are unavailable for most MW satellites. Of the 4 satellites with published proper motions, only the Small Magellanic Cloud is clearly associated with the LMC. Taken at face value, the proper motions of Carina, Fornax and Sculptor rule them out as potential associates, but this conclusion should be revisited when better data become available. The dearth of satellites clearly associated with the Clouds might be solved by wide-field imaging surveys that target its surroundings, a region that may prove a fertile hunting ground for faint, previously unnoticed MW satellites.

preprint2011arXiv

SDSS-III: Massive Spectroscopic Surveys of the Distant Universe, the Milky Way Galaxy, and Extra-Solar Planetary Systems

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

preprint2011arXiv

Streams and caustics: the fine-grained structure of LCDM haloes

We present the first and so far the only simulations to follow the fine-grained phase-space structure of galaxy haloes formed from generic LCDM initial conditions. We integrate the geodesic deviation equation in tandem with the N-body equations of motion, demonstrating that this can produce numerically converged results for the properties of fine-grained phase-space streams and their associated caustics, even in the inner regions of haloes. Our effective resolution for such structures is many orders of magnitude better than achieved by conventional techniques on even the largest simulations. We apply these methods to the six Milky Way-mass haloes of the Aquarius Project. At 8 kpc from halo centre a typical point intersects about 10^14 streams with a very broad range of individual densities; the ~10^6 most massive streams contribute about half of the local dark matter density. As a result, the velocity distribution of dark matter particles should be very smooth with the most massive fine-grained stream contributing about 0.1% of the total signal. Dark matter particles at this radius have typically passed 200 caustics since the Big Bang. The peak densities on present-day caustics in the inner halo almost all lie well below the mean local dark matter density. As a result caustics provide a negligible boost (<0.1%) to the predicted local dark matter annihilation rate. The effective boost is larger in the outer halo but never exceeds about 10%. Thus fine-grained streams and their associated caustics have no effect on the detectability of dark matter, either directly in Earth-bound laboratories, or indirectly through annihilation radiation, with the exception that resonant cavity experiments searching for axions may see the most massive local fine-grained streams because of their extreme localisation in energy/momentum space. (abridged)

preprint2011arXiv

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

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.

preprint2011arXiv

The Eighth Data Release of the Sloan Digital Sky Survey: First Data from SDSS-III

The Sloan Digital Sky Survey (SDSS) started a new phase in August 2008, with new instrumentation and new surveys focused on Galactic structure and chemical evolution, measurements of the baryon oscillation feature in the clustering of galaxies and the quasar Ly alpha forest, and a radial velocity search for planets around ~8000 stars. This paper describes the first data release of SDSS-III (and the eighth counting from the beginning of the SDSS). The release includes five-band imaging of roughly 5200 deg^2 in the Southern Galactic Cap, bringing the total footprint of the SDSS imaging to 14,555 deg^2, or over a third of the Celestial Sphere. All the imaging data have been reprocessed with an improved sky-subtraction algorithm and a final, self-consistent photometric recalibration and flat-field determination. This release also includes all data from the second phase of the Sloan Extension for Galactic Understanding and Evolution (SEGUE-2), consisting of spectroscopy of approximately 118,000 stars at both high and low Galactic latitudes. All the more than half a million stellar spectra obtained with the SDSS spectrograph have been reprocessed through an improved stellar parameters pipeline, which has better determination of metallicity for high metallicity stars.

preprint2011arXiv

The Shape of Dark Matter Haloes in the Aquarius Simulations: Evolution and Memory

We use the high resolution cosmological N-body simulations from the Aquarius project to investigate in detail the mechanisms that determine the shape of Milky Way-type dark matter haloes. We find that, when measured at the instantaneous virial radius, the shape of individual haloes changes with time, evolving from a typically prolate configuration at early stages to a more triaxial/oblate geometry at the present day. This evolution in halo shape correlates well with the distribution of the infalling material: prolate configurations arise when haloes are fed through narrow filaments, which characterizes the early epochs of halo assembly, whereas triaxial/oblate configurations result as the accretion turns more isotropic at later times. Interestingly, at redshift z=0, clear imprints of the past history of each halo are recorded in their shapes at different radii, which also exhibit a variation from prolate in the inner regions to triaxial/oblate in the outskirts. Provided that the Aquarius haloes are fair representatives of Milky Way-like 10^12 Msun objects, we conclude that the shape of such dark matter haloes is a complex, time-dependent property, with each radial shell retaining memory of the conditions at the time of collapse.

preprint2010arXiv

An observer's view of simulated galaxies: disc-to-total ratios, bars, and (pseudo-)bulges

We use cosmological hydrodynamical simulations of the formation of Milky Way-mass galaxies to study the relative importance of the main stellar components, i.e., discs, bulges, and bars, at redshift zero. The main aim of this work is to understand if estimates of the structural parameters of these components determined from kinematics (as is usually done in simulations) agree well with those obtained using a photometric bulge/disc/bar decomposition (as done in observations). To perform such a comparison, we have produced synthetic observations of the simulation outputs with the Monte-Carlo radiative transfer code SUNRISE and used the BUDDA code to make 2D photometric decompositions of the resulting images (in the i and g bands). We find that the kinematic disc-to-total ratio (D/T) estimates are systematically and significantly lower than the photometric ones. While the maximum D/T ratios obtained with the former method are of the order of 0.2, they are typically >0.4, and can be as high as 0.7, according to the latter. The photometric decomposition shows that many of the simulated galaxies have bars, with Bar/T ratios in the range 0.2-0.4, and that bulges have in all cases low Sersic indices, resembling observed pseudo-bulges instead of classical ones. Simulated discs, bulges and bars generally have similar (g-i) colours, which are in the blue tail of the distribution of observed colours. This is not due to the presence of young stars, but rather to low metallicities and poor gas content in the simulated galaxies, which makes dust extinction low. Photometric decompositions thus match the component ratios usually quoted for spiral galaxies better than kinematic decompositions, but the shift is insufficient to make the simulations consistent with observed late-type systems.

preprint2010arXiv

Autocorrelations of stellar light and mass in the low-redshift Universe

The final data release of the Sloan Digital Sky Survey (SDSS) provides reliable photometry and spectroscopy for about half a million galaxies with median redshift 0.09. Here we use these data to estimate projected autocorrelation functions w_p(r_p) for the light of galaxies in the five SDSS photometric bands. Comparison with the analogous stellar mass autocorrelation, estimated in a previous paper, shows that stellar luminosity is less strongly clustered than stellar mass in all bands and on all scales. Over the full nonlinear range 10 kpc/h < r_p < 10 Mpc/h our autocorrelation estimates are extremely well represented by power laws. The parameters of the corresponding spatial functions ξ(r) = (r/r_0)^γvary systematically from r_0=4.5 Mpc/h and γ=-1.74 for the bluest band (the u band) to r_0=5.8 Mpc/h and γ=-1.83 for the reddest one (the z band). These may be compared with r_0=6.1 Mpc/h and γ=-1.84 for the stellar mass. Ratios of w_p(r_p) between two given wavebands are proportional to the mean colour of correlated stars at projected distance r_p from a randomly chosen star. The ratio of the stellar mass and luminosity autocorrelations measures an analogous mean stellar mass-to-light ratio (M*/L). All colours get redder and all mass-to-light ratios get larger with decreasing r_p, with the amplitude of the effects decreasing strongly to redder passbands. Even for the u-band the effects are quite modest, with maximum shifts of about 0.1 in u-g and about 25% in M*/L_u. These trends provide a precise characterisation of the well-known dependence of stellar populations on environment.

preprint2010arXiv

Bound and unbound substructures in Galaxy-scale Dark Matter haloes

We analyse the coarse-grained phase-space structure of the six Galaxy-scale dark matter haloes of the Aquarius Project using a state-of-the-art 6D substructure finder. Within r_50, we find that about 35% of the mass is in identifiable substructures, predominantly tidal streams, but including about 14% in self-bound subhaloes. The slope of the differential substructure mass function is close to -2, which should be compared to around -1.9 for the population of self-bound subhaloes. Near r_50 about 60% of the mass is in substructures, with about 30% in self-bound subhaloes. The inner 35 kpc of the highest resolution simulation has only 0.5% of its mass in self-bound subhaloes, but 3.3% in detected substructure, again primarily tidal streams. The densest tidal streams near the solar position have a 3-D mass density about 1% of the local mean, and populate the high velocity tail of the velocity distribution.

preprint2010arXiv

Dark matter response to galaxy formation

We have resimulated the six galaxy-sized haloes of the Aquarius Project including metal-dependent cooling, star formation and supernova feedback. This allows us to study not only how dark matter haloes respond to galaxy formation, but also how this response is affected by details of halo assembly history. In agreement with previous work, we find baryon condensation to lead to increased dark matter concentration. Dark matter density profiles differ substantially in shape from halo to halo when baryons are included, but in all cases the velocity dispersion decreases monotonically with radius. Some haloes show an approximately constant dark matter velocity anisotropy with $ β\approx 0.1-02$, while others retain the anisotropy structure of their baryon-free versions. Most of our haloes become approximately oblate in their inner regions, although a few retain the shape of their dissipationless counterparts. Pseudo-phase-space densities are described by a power law in radius of altered slope when baryons are included. The shape and concentration of the dark matter density profiles are not well reproduced by published adiabatic contraction models. The significant spread we find in the density and kinematic structure of our haloes appears related to differences in their formation histories. Such differences already affect the final structure in baryon-free simulations, but they are reinforced by the inclusion of baryons, and new features are produced. The details of galaxy formation need to be better understood before the inner dark matter structure of galaxies can be used to constrain cosmological models or the nature of dark matter.

preprint2010arXiv

Kinematics of Milky Way Satellites in a Lambda Cold Dark Matter Universe

We investigate whether the subhalos of Lambda-CDM galaxy halos have potentials consistent with the observed properties of Milky Way satellites, particularly those with high-quality photometric and kinematic data: Fornax, Leo I, Sculptor, Sextans, and Carina. We compare spherical models with isotropic velocity dispersion tensors to the observed, circularly averaged star counts, line-of-sight velocity dispersion profiles and line-of-sight velocity distributions. We identify subhalos within the six high-resolution dark matter halos of the Aquarius Project for which the spherically averaged potentials result in excellent fits to each of the five galaxies. In particular, our simple one-integral models reproduce the observations in the inner regions, proving that these data are fully consistent with Lambda-CDM expectations and do not require cored dark matter distributions. For four of the five satellites the fits require moderately cusped {\it stellar} density profiles. The star count data for Leo I, however, do require a cored distribution of star counts. Current data suggest that these five satellites may be hosted by Lambda-CDM subhalos with maximum circular velocities in the range 10 to 30 km/s.

preprint2010arXiv

Nonspherical similarity solutions for dark halo formation

We carry out fully 3-dimensional simulations of evolution from self-similar, spherically symmetric linear perturbations of a Cold Dark Matter dominated Einstein-de Sitter universe. As a result of the radial orbit instability, the haloes which grow from such initial conditions are triaxial with major-to-minor axis ratios of order 3:1. They nevertheless grow approximately self-similarly in time. In all cases they have power-law density profiles and near-constant velocity anisotropy in their inner regions. Both the power-law index and the value of the velocity anisotropy depend on the similarity index of the initial conditions, the former as expected from simple scaling arguments. Halo structure is thus not "universal" but remembers the initial conditions. On larger scales the density and anisotropy profiles show two characteristic scales, corresponding to particles at first pericentre and at first apocentre after infall. They are well approximated by the NFW model only for one value of the similarity index. In contrast, at all radii within the outer caustic the pseudo phase-space density can be fit by a single power law with an index which depends only very weakly on the similarity index of the initial conditions. This behaviour is very similar to that found for haloes formed from LCDM initial conditions and so can be considered approximately universal.

preprint2010arXiv

Relic density and CMB constraints on dark matter annihilation with Sommerfeld enhancement

We calculate how the relic density of dark matter particles is altered when their annihilation is enhanced by the Sommerfeld mechanism due to a Yukawa interaction between the annihilating particles. Maintaining a dark matter abundance consistent with current observational bounds requires the normalization of the s-wave annihilation cross section to be decreased compared to a model without enhancement. The level of suppression depends on the specific parameters of the particle model, with the kinetic decoupling temperature having the most effect. We find that the cross section can be reduced by as much as an order of magnitude for extreme cases. We also compute the mu-type distortion of the CMB energy spectrum caused by energy injection from such Sommerfeld-enhanced annihilation. Our results indicate that in the vicinity of resonances, associated with bound states, distortions can be large enough to be excluded by the upper limit |mu|<9.0x10^(-5) found by the COBE/FIRAS experiment.

preprint2010arXiv

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

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.

preprint2010arXiv

The evolution of the density of galaxy clusters and groups: denser environments at higher redshifts

We show that, observationally, the projected local density distribution in high-z clusters is shifted towards higher values compared to clusters at lower redshift. To search for the origin of this evolution, we analyze a sample of haloes selected from the Millennium Simulation and populated using semi-analytic models, investigating the relation between observed projected density and physical 3D density, using densities computed from the 10 and 3 closest neighbours. Both observationally and in the simulations, we study the relation between number of cluster members and cluster mass, and number of members per unit of cluster mass. We find that the observed evolution of projected densities reflects a shift to higher values of the physical 3D density distribution. In turn, this must be related with the globally higher number of galaxies per unit of cluster volume N/V in the past. We show that the evolution of N/V is due to a combination of two effects: a) distant clusters were denser in dark matter (DM) simply because the DM density within R_{200} (~the cluster virial radius) is defined to be a fixed multiple of the critical density of the Universe, and b) the number of galaxies per unit of cluster DM mass is remarkably constant both with redshift and cluster mass if counting galaxies brighter than a passively evolving magnitude limit. Our results highlight that distant clusters were much denser environments than today's clusters, both in galaxy number and mass, and that the density conditions felt by galaxies in virialized systems do not depend on the system mass.

preprint2010arXiv

There's no place like home? Statistics of Milky Way-mass dark matter halos

We present an analysis of the distribution of structural properties for Milky Way-mass halos in the Millennium-II Simulation (MS-II). This simulation of structure formation within the standard LCDM cosmology contains thousands of Milky Way-mass halos and has sufficient resolution to properly resolve many subhalos per host. It thus provides a major improvement in the statistical power available to explore the distribution of internal structure for halos of this mass. In addition, the MS-II contains lower resolution versions of the Aquarius Project halos, allowing us to compare our results to simulations of six halos at a much higher resolution. We study the distributions of mass assembly histories, of subhalo mass functions and accretion times, and of merger and stripping histories for subhalos capable of impacting disks at the centers of halos. We show that subhalo abundances are not well-described by Poisson statistics at low mass, but rather are dominated by intrinsic scatter. Using the masses of subhalos at infall and the abundance-matching assumption, there is less than a 10% chance that a Milky Way halo with M_vir =10^12 M_sun will host two galaxies as bright as the Magellanic Clouds. This probability rises to ~25% for a halo with M_vir=2.5 x 10^12 M_sun. The statistics relevant for disk heating are very sensitive to the mass range that is considered relevant. Mergers with infall mass : redshift zero virial mass greater than 1:30 could well impact a central galactic disk and are a near inevitability since z=2, whereas only half of all halos have had a merger with infall mass : redshift zero virial mass greater than 1:10 over this same period.

preprint2010arXiv

What is the (Dark) Matter with Dwarf Galaxies?

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

preprint2009arXiv

Abundances, masses, and weak-lensing mass profiles of galaxy clusters as a function of richness and luminosity in LambdaCDM cosmologies

We test the concordance LCDM cosmology by comparing predictions for the mean properties of galaxy clusters to observations. We use high-resolution N-body simulations of cosmic structure formation and semi-analytic models of galaxy formation to compute the abundance, mean density profile, and mass of galaxy clusters as a function of richness and luminosity, and we compare these predictions to observations of clusters in the Sloan Digital Sky Survey (SDSS) maxBCG catalogue. We discuss the scatter in the mass-richness relation, the reconstruction of the cluster mass function from the mass-richness relation, and fits to the weak-lensing cluster mass profiles. The impact of cosmological parameters on the predictions is investigated by comparing results from galaxy models based on the Millennium Simulation (MS) and another WMAP1 simulation to those from a WMAP3 simulation. We find that the simulated weak-lensing mass profiles and the observed profiles of the SDSS maxBCG clusters agree well in shape and amplitude. The mass-richness relations in the simulations are close to the observed relation, with differences lesssim 30%. The MS and WMAP1 simulations yield cluster abundances similar to those observed, whereas abundances in the WMAP3 simulation are 2-3 times lower. The differences in cluster abundance, mass, and density amplitude between the simulations and the observations can be attributed to differences in the underlying cosmological parameters, in particular the power spectrum normalisation sigma_8. Better agreement between predictions and observations should be reached with a normalisation $0.722<sigma_8<0.9$ (probably closer to the upper value), i.e. between the values underlying the two simulation sets.

preprint2009arXiv

Assembly bias and the dynamical structure of dark matter halos

Based on the Millennium Simulation we examine assembly bias for the halo properties: shape, triaxiality, concentration, spin, shape of the velocity ellipsoid and velocity anisotropy. For consistency we determine all these properties using the same set of particles, namely all gravitationally self-bound particles belonging to the most massive sub-structure of a given friends-of-friends halo. We confirm that near-spherical and high-spin halos show enhanced clustering. The opposite is true for strongly aspherical and low-spin halos. Further, below the typical collapse mass, M*, more concentrated halos show stronger clustering whereas less concentrated halos are less clustered which is reversed for masses above M*. Going beyond earlier work we show that: (1) oblate halos are more strongly clustered than prolate ones; (2) the dependence of clustering on the shape of the velocity ellipsoid coincides with that of the real-space shape, although the signal is stronger; (3) halos with weak velocity anisotropy are more clustered, whereas radially anisotropic halos are more weakly clustered; (4) for all highly clustered subsets we find systematically less radially biased velocity anisotropy profiles. These findings indicate that the velocity structure of halos is tightly correlated with environment.

preprint2009arXiv

Caustics in growing Cold Dark Matter Haloes

We simulate the growth of isolated dark matter haloes from self-similar and spherically symmetric initial conditions. Our N-body code integrates the geodesic deviation equation in order to track the streams and caustics associated with individual simulation particles. The radial orbit instability causes our haloes to develop major-to-minor axis ratios approaching 10 to 1 in their inner regions. They grow similarly in time and have similar density profiles to the spherical similarity solution, but their detailed structure is very different. The higher dimensionality of the orbits causes their stream and caustic densities to drop much more rapidly than in the similarity solution. This results in a corresponding increase in the number of streams at each point. At 1% of the turnaround radius (corresponding roughly to the Sun's position in the Milky Way) we find of order 10^6 streams in our simulations, as compared to 10^2 in the similarity solution. The number of caustics in the inner halo increases by a factor of several, because a typical orbit has six turning points rather than one, but caustic densities drop by a much larger factor. This reduces the caustic contribution to the annihilation radiation. For the region between 1% and 50% of the turnaround radius, this is 4% of the total in our simulated haloes, as compared to 6.5% in the similarity solution. Caustics contribute much less at smaller radii. These numbers assume a 100 GeV c^-2 neutralino with present-day velocity dispersion 0.03 cm s^-1, but reducing the dispersion by ten orders of magnitude only doubles the caustic luminosity. We conclude that caustics will be unobservable in the inner parts of haloes. Only the outermost caustic might potentially be detectable.

preprint2009arXiv

Cosmic Cartography of the Large-Scale Structure with Sloan Digital Sky Survey Data Release 6

We present the largest Wiener reconstruction of the cosmic density field made to date. The reconstruction is based on the Sloan Digital Sky Survey data release 6 covering the northern Galactic cap. We use a novel supersampling algorithm to suppress aliasing effects and a Krylov-space inversion method to enable high performance with high resolution. These techniques are implemented in the ARGO computer code. We reconstruct the field over a 500 Mpc cube with Mpc grid-resolution while accounting both for the angular and radial selection functions of the SDSS, and the shot noise giving an effective resolution of the order of ~10 Mpc. In addition, we correct for the redshift distortions in the linear and nonlinear regimes in an approximate way. We show that the commonly used method of inverse weighting the galaxies by the corresponding selection function heads to excess noise in regions where the density of the observed galaxies is small. It is more accurate and conservative to adopt a Bayesian framework in which we model the galaxy selection/detection process to be Poisson-binomial. This results in heavier smoothing in regions of reduced sampling density. Our results show a complex cosmic web structure with huge void regions indicating that the recovered matter distribution is highly non-Gaussian. Filamentary structures are clearly visible on scales up to ~20 Mpc. We also calculate the statistical distribution of density after smoothing the reconstruction with Gaussian kernels of different radii r_S and find good agreement with a log-normal distribution for ~10 Mpc < r_S < ~30 Mpc.

preprint2009arXiv

Evolution of the Early-Type Galaxy Fraction in Clusters since z = 0.8

We study the morphological content of a large sample of high-redshift clusters to determine its dependence on cluster mass and redshift. Quantitative morphologies are based on bulge+disk decompositions of cluster and field galaxies on deep VLT/FORS2 images of 18 optically-selected clusters at 0.45 < z < 0.80 from the ESO Distant Cluster Survey (EDisCS). Morphological content is given by the early-type galaxy fraction f_et, and early-type galaxies are selected based on their bulge fraction and image smoothness. A set of 158 SDSS clusters is analyzed exactly as the EDisCS sample to provide a robust local comparison. Our main results are: (1) f_et values for the SDSS and EDisCS clusters exhibit no clear trend as a function of sigma. (2) Mid-z EDisCS clusters around sigma = 500 km/s have f_et ~= 0.5 whereas high-z EDisCS clusters have f_et ~= 0.4 (~25% increase over 2 Gyrs). (3) There is a marked difference in the morphological content of EDisCS and SDSS clusters. None of the EDisCS clusters have f_et greater than 0.6 whereas half of the SDSS clusters lie above this value. This difference is seen in clusters of all velocity dispersions. (4) There is a strong correlation between morphology and star formation in SDSS and EDisCS clusters. This correlation holds independent of sigma and z even though the fraction of [OII] emitters decreases from z~0.8 to z~0.06 in all environments. Our results pose an interesting challenge to structural transformation and star formation quenching processes that strongly depend on the global cluster environment and suggest that cluster membership may be of lesser importance than other variables in determining galaxy properties. (ABRIDGED)

preprint2009arXiv

Formation of Isolated Dwarf Galaxies with Feedback

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

preprint2009arXiv

Galaxies-Intergalactic Medium Interaction Calculation --I. Galaxy formation as a function of large-scale environment

[Abridged] We present the first results of hydrodynamical simulations that follow the formation of galaxies to z=0 in spherical regions of radius ~20 Mpc/h drawn from the Millennium Simulation. The regions have overdensities that deviate by (-2, -1, 0, +1, +2)sigma from the cosmic mean, where sigma is the rms mass fluctuation on a scale of ~20Mpc/h at z=1.5. The simulations have mass resolution of up to 10^6 Msun/h, cover the entire range of large-scale environments and allow extrapolation of statistics to the entire 500 (Mpc/h)^3 Millennium volume. They include gas cooling, photoheating from an ionising background, SNe feedback and winds, but no AGN. We find that the specific SFR density at z <~ 10 varies systematically from region to region by up to an order of magnitude, but the global value, averaged over all volumes, reproduces observational data. Massive, compact galaxies, similar to those observed in the GOODS fields, form in the overdense regions as early as z=6, but do not appear in the underdense regions until z~3. These environmental variations are not caused by a dependence of the star formation properties on environment, but rather by a strong variation of the halo mass function from one environment to another, with more massive haloes forming preferentially in the denser regions. At all epochs, stars form most efficiently in haloes of circular velocity ~ 250 km/s. However, the star formation history exhibits a form of "downsizing" (even in the absence of AGN): the stars comprising massive galaxies at z=0 have mostly formed by z=1-2, whilst those comprising smaller galaxies typically form at later times. However, additional feedback is required to limit star formation in massive galaxies at late times.

preprint2009arXiv

Modeling the cosmological co-evolution of supermassive black holes and galaxies: II. The clustering of quasars and their dark environment

We use semi-analytic modeling on top of the Millennium simulation to study the joint formation of galaxies and their embedded supermassive black holes. Our goal is to test scenarios in which black hole accretion and quasar activity are triggered by galaxy mergers, and to constrain different models for the lightcurves associated with individual quasar events. In the present work we focus on studying the spatial distribution of simulated quasars. At all luminosities, we find that the simulated quasar two-point correlation function is fit well by a single power-law in the range 0.5 < r < 20 h^{-1} Mpc, but its normalization is a strong function of redshift. When we select only quasars with luminosities within the range typically accessible by today's quasar surveys, their clustering strength depends only weakly on luminosity, in agreement with observations. This holds independently of the assumed lightcurve model, since bright quasars are black holes accreting close to the Eddington limit, and are hosted by dark matter haloes with a narrow mass range of a few 10^12 h^{-1} M_sun. Therefore the clustering of bright quasars cannot be used to disentangle lightcurve models, but such a discrimination would become possible if the observational samples can be pushed to significantly fainter limits. Overall, our clustering results for the simulated quasar population agree rather well with observations, lending support to the conjecture that galaxy mergers could be the main physical process responsible for triggering black hole accretion and quasar activity.

preprint2009arXiv

One simulation to fit them all - changing the background parameters of a cosmological N-body simulation

We demonstrate that the output of a cosmological N-body simulation can, to remarkable accuracy, be scaled to represent the growth of large-scale structure in a cosmology with parameters similar to but different from those originally assumed. Our algorithm involves three steps: a reassignment of length, mass and velocity units, a relabelling of the time axis, and a rescaling of the amplitudes of individual large-scale fluctuation modes. We test it using two matched pairs of simulations. Within each pair, one simulation assumes parameters consistent with analyses of the first-year WMAP data. The other has lower matter and baryon densities and a 15% lower fluctuation amplitude, consistent with analyses of the three-year WMAP data. The pairs differ by a factor of a thousand in mass resolution, enabling performance tests on both linear and nonlinear scales. Our scaling reproduces the mass power spectra of the target cosmology to better than 0.5% on large scales (k < 0.1 h/Mpc) both in real and in redshift space. In particular, the BAO features of the original cosmology are removed and are correctly replaced by those of the target cosmology. Errors are still below 3% for k < 1 h/Mpc. Power spectra of the dark halo distribution are even more precisely reproduced, with errors below 1% on all scales tested. A halo-by-halo comparison shows that centre-of-mass positions and velocities are reproduced to better than 90 kpc/h and 5%, respectively. Halo masses, concentrations and spins are also reproduced at about the 10% level, although with small biases. Halo assembly histories are accurately reproduced, leading to central galaxy magnitudes with errors of about 0.25 magnitudes and a bias of about 0.13 magnitudes for a representative semi-analytic model.

preprint2009arXiv

Star formation and AGN activity in SDSS cluster galaxies

We investigate the recent and current star formation activity of galaxies as function of distance from the cluster center in a sample of 521 SDSS clusters at z<0.1. We show that when the BCGs are excluded from the galaxy sample, there is no evidence for mass segregation in the clusters, so that differences in cluster and field populations cannot simply be attributed to different mass functions. We find a marked star formation-radius relation in that almost all galaxies in the cluster core are quiescent, i.e. have terminated star formation a few Gyr ago. This star formation-radius relation is most pronounced for low-mass galaxies and is very weak or absent beyond the virial radius. The typical star formation rate of non-quiescent galaxies declines by approximately a factor of two towards the cluster center. However, the fraction of galaxies with young stellar populations indicating a recently completed starburst or a truncation of star formation does not vary significantly with radius. The fraction of star-forming galaxies which host a powerful optical AGN is also independent of clustercentric radius, indicating that the link between star formation and AGN in these galaxies operates independent of environment. The fraction of red galaxies which host a weak optical AGN decreases, however, towards the cluster center, with a similar timescale as the decline of star forming galaxies. Our results favor a scenario in which star formation is quenched slowly, on timescales similar to the cluster crossing time, and rule out significant contributions from more violent processes at least beyond cluster radii >0.1 R_200.

preprint2009arXiv

The distribution of stellar mass in the low-redshift Universe

We use a complete and uniform sample of almost half a million galaxies from the Sloan Digital Sky Survey to characterise the distribution of stellar mass in the low-redshift Universe. Galaxy abundances are well determined over almost four orders of magnitude in stellar mass, and are reasonably but not perfectly fit by a Schechter function with characteristic stellar mass m* = 6.7 x 10^10 M_sun and with faint-end slope α= -1.155. For a standard cosmology and a standard stellar Initial Mass Function, only 3.5% of the baryons in the low-redshift Universe are locked up in stars. The projected autocorrelation function of stellar mass is robustly and precisely determined for r_p < 30 Mpc/h. Over the range 10 kpc/kpc < r_p < 10 Mpc/h it is extremely well represented by a power law. The corresponding three-dimensional autocorrelation function is ξ*(r) = (r/6.1 Mpc/h)^{-1.84}. Relative to the dark matter, the bias of the stellar mass distribution is approximately constant on large scales, but varies by a factor of five for r_p < 1 Mpc/h. This behaviour is approximately but not perfectly reproduced by current models for galaxy formation in the concordance LCDM cosmology. Detailed comparison suggests that a fluctuation amplitude σ_8 ~ 0.8 is preferred to the somewhat larger value adopted in the Millennium Simulation models with which we compare our data. This comparison also suggests that observations of stellar mass autocorrelations as a function of redshift might provide a powerful test for the nature of Dark Energy.

preprint2009arXiv

The Diversity and Similarity of Simulated Cold Dark Matter Halos

We study the mass, velocity dispersion, and anisotropy profiles of $Λ$CDM halos using a suite of N-body simulations of unprecedented numerical resolution (the {\it Aquarius Project}). Our analysis confirms a number of results claimed by earlier work, and clarifies a few issues where conflicting claims may be found in the recent literature. The spherically-averaged density profile becomes progressively shallower inwards and, at the innermost resolved radius, the logarithmic slope is $γ\equiv -$d$\lnρ/$d$\ln r \simlt 1$. Asymptotic inner slopes as steep as the recently claimed $ρ\propto r^{-1.2}$ are clearly ruled out. The radial dependence of $γ$ is well approximated by a power-law, $γ\propto r^α$ (the Einasto profile). The shape parameter, $α$, varies slightly but significantly from halo to halo, implying that the mass profiles of $Λ$CDM halos are not strictly universal: different halos cannot, in general, be rescaled to look identical. Departures from similarity are also seen in velocity dispersion profiles and correlate with those in density profiles so as to preserve a power-law form for the spherically averaged pseudo-phase-space density, $ρ/σ^3\propto r^{-1.875}$. Our conclusions are reliable down to radii below 0.4% of the virial radius, providing well-defined predictions for halo structure when baryonic effects are neglected, and thus an instructive theoretical template against which the modifications induced by the baryonic components of real galaxies can be judged.

preprint2009arXiv

The Formation and Survival of Discs in a Lambda-CDM Universe

We study the formation of galaxies in a Lambda-CDM Universe using high resolution hydrodynamical simulations with a multiphase treatment of gas, cooling and feedback, focusing on the formation of discs. Our simulations follow eight haloes similar in mass to the Milky Way and extracted from a large cosmological simulation without restriction on spin parameter or merger history. This allows us to investigate how the final properties of the simulated galaxies correlate with the formation histories of their haloes. We find that, at z = 0, none of our galaxies contain a disc with more than 20 per cent of its total stellar mass. Four of the eight galaxies nevertheless have well-formed disc components, three have dominant spheroids and very small discs, and one is a spheroidal galaxy with no disc at all. The z = 0 spheroids are made of old stars, while discs are younger and formed from the inside-out. Neither the existence of a disc at z = 0 nor the final disc-to-total mass ratio seems to depend on the spin parameter of the halo. Discs are formed in haloes with spin parameters as low as 0.01 and as high as 0.05; galaxies with little or no disc component span the same range in spin parameter. Except for one of the simulated galaxies, all have significant discs at z > ~2, regardless of their z = 0 morphologies. Major mergers and instabilities which arise when accreting cold gas is misaligned with the stellar disc trigger a transfer of mass from the discs to the spheroids. In some cases, discs are destroyed, while in others, they survive or reform. This suggests that the survival probability of discs depends on the particular formation history of each galaxy. A realistic Lambda-CDM model will clearly require weaker star formation at high redshift and later disc assembly than occurs in our models.

preprint2008arXiv

A census of metals and baryons in stars in the local Universe

We combine stellar metallicity and stellar mass estimates for a large sample of galaxies drawn from the SDSS DR2 spanning wide ranges in physical properties, in order to derive an inventory of the total mass of metals and baryons locked up in stars today. Physical parameter estimates are derived from galaxy spectra with high S/N (>20). Coadded spectra of galaxies with similar velocity dispersions, absolute r-band magnitudes and 4000Å-break values are used for those regions of parameter space where individual spectra have lower S/N. We estimate the total density of metals and of baryons in stars and, from these two quantities, we obtain a mass- and volume-averaged stellar metallicity of <Z_star>=1.04+-0.14 Z_sun, i.e. consistent with solar. We also study how metals are distributed in galaxies according to their mass, morphology and age, and we then compare these distributions with the corresponding distributions of stellar mass. We find that the bulk of metals locked up in stars in the local Universe reside in massive, bulge-dominated galaxies, with red colours and high 4000Å-break values corresponding to old stellar populations. Bulge-dominated and disc-dominated galaxies contribute similar amounts to the total stellar mass density, but have different fractional contributions to the mass density of metals in stars, in agreement with the mass-metallicity relation. Bulge-dominated galaxies contain roughly 40% of the total amount of metals in stars, while disc-dominated galaxies less than 25%. Finally, at a given galaxy stellar mass, we define two characteristic ages as the median of the distributions of mass and metals as a function of age. These characteristic ages decrease progressively from high-mass to low-mass galaxies, consistent with the high formation epochs of stars in massive galaxies.

preprint2008arXiv

Are mergers responsible for universal halo properties?

N-body simulations of Cold Dark Matter (CDM) have shown that, in this hierarchical structure formation model, dark matter halo properties, such as the density profile, the phase-space density profile, the distribution of axial ratio, the distribution of spin parameter, and the distribution of internal specific angular momentum follow `universal' laws or distributions. Here we study the properties of the first generation of haloes in a Hot Dark Matter (HDM) dominated universe, as an example of halo formation through monolithic collapse. We find all these universalities to be present in this case also. Halo density profiles are very well fit by the Navarro et al (1997) profile over two orders of magnitude in mass. The concentration parameter depends on mass as $c \propto M^{0.2}$,reversing the dependence found in a hierarchical CDM universe. However, the concentration-formation time relation is similar in the two cases: earlier forming haloes tend to be more concentrated than their later forming counterparts. Halo formation histories are also characterized by two phases in the HDM case: an early phase of rapid accretion followed by slower growth. Furthermore, there is no significant difference between the HDM and CDM cases concerning the statistics of other halo properties: the phase-space density profile; the velocity anisotropy profile; the distribution of shape parameters; the distribution of spin parameter, and the distribution of internal specific angular momentum are all similar in the two cases. Only substructure content differs dramatically. These results indicate that mergers do not play a pivotal role in establishing the universalities, thus contradicting models which explain them as consequences of mergers.

preprint2008arXiv

Effects of Supernova Feedback on the Formation of Galaxies

We study the effects of Supernova (SN) feedback on the formation of galaxies using hydrodynamical simulations in a Lambda-CDM cosmology. We use an extended version of the code GADGET-2 which includes chemical enrichment and energy feedback by Type II and Type Ia SN, metal-dependent cooling and a multiphase model for the gas component. We focus on the effects of SN feedback on the star formation process, galaxy morphology, evolution of the specific angular momentum and chemical properties. We find that SN feedback plays a fundamental role in galaxy evolution, producing a self-regulated cycle for star formation, preventing the early consumption of gas and allowing disks to form at late times. The SN feedback model is able to reproduce the expected dependence on virial mass, with less massive systems being more strongly affected.

preprint2008arXiv

High Redshift Galaxy Populations and their Descendants

We study model predictions for three high-redshift galaxy populations: Lyman break galaxies at z~3 (LBGs), optically selected star-forming galaxies at z~2 (BXs), and distant red galaxies at z~2 (DRGs).Our galaxy formation model simultaneously reproduces the abundances, redshift distributions and clustering of all three observed populations. The star formation rates (SFRs) of model LBGs and BXs are lower than those quoted for real samples, reflecting different initial mass functions and scatter in model dust properties. About 85% of model galaxies selected as DRGs are star-forming, with SFRs ranging up to 100 M_sun/yr. Model LBGs, BXs and DRGs together account for less than half of all star formation over the range 1.5<z<3.2. Model BXs have metallicities which agree roughly with observation, but model LBGs are only slightly more metal-poor, in disagreement with recent observational results. Model galaxies are predominantly disk-dominated. About 30% of model galaxies with M>10^{11}M_sun are classified as LBGs or BXs at the relevant redshifts, while 65% are classified as DRGs. Almost all model LBGs and BXs are central galaxies, but about a quarter of DRGs are satellites. Half of all LBG descendants at z=2 would be identified as BX's, but very few as DRGs. Clustering increases with decreasing redshift for descendants of all three populations, becoming stronger than that of L^* galaxies by z=0, when many have become satellite galaxies. Their stellar mass growth is dominated by star formation until z~1 and thereafter by mergers. Most LBGs and DRGs end up as red ellipticals, while BXs have a more varied fate. 99% of local galaxies with M>10^{11.5}M_sun are predicted to have at least one LBG/BX/DRG progenitor, and over 70% above 10^{11}M_sun. (abbreviated)

preprint2008arXiv

The Aquarius Project: the subhalos of galactic halos

We have performed the largest ever particle simulation of a Milky Way-sized dark matter halo, and present the most comprehensive convergence study for an individual dark matter halo carried out thus far. We have also simulated a sample of 6 ultra-highly resolved Milky-way sized halos, allowing us to estimate the halo-to-halo scatter in substructure statistics. In our largest simulation, we resolve nearly 300,000 gravitationally bound subhalos within the virialized region of the halo. Simulations of the same object differing in mass resolution by factors up to 1800 accurately reproduce the largest subhalos with the same mass, maximum circular velocity and position, and yield good convergence for the abundance and internal properties of dark matter substructures. We detect up to four generations of subhalos within subhalos, but contrary to recent claims, we find less substructure in subhalos than in the main halo when regions of equal mean overdensity are compared. The overall substructure mass fraction is much lower in subhalos than in the main halo. Extrapolating the main halo's subhalo mass spectrum down to an Earth mass, we predict the mass fraction in substructure to be well below 3% within 100 kpc, and to be below 0.1% within the Solar Circle. The inner density profiles of subhalos show no sign of converging to a fixed asymptotic slope and are well fit by gently curving profiles of Einasto form. The mean concentrations of isolated halos are accurately described by the fitting formula of Neto et al. down to maximum circular velocities of 1.5 km/s, an extrapolation over some 5 orders of magnitude in mass. However, at equal maximum circular velocity, subhalos are more concentrated than field halos, with a characteristic density that is typically ~2.6 times larger and increases towards the halo centre.

preprint2007arXiv

Masses for the Local Group and the Milky Way

We use the very large Millennium Simulation of the concordance $Λ$CDM cosmogony to calibrate the bias and error distribution of Timing Argument estimators of the masses of the Local Group and of the Milky Way. From a large number of isolated spiral-spiral pairs similar to the Milky Way/Andromeda system, we find the interquartile range of the ratio of timing mass to true mass to be a factor of 1.8, while the 5% and 95% points of the distribution of this ratio are separated by a factor of 5.7. Here we define true mass as the sum of the ``virial'' masses $M_{200}$ of the two dominant galaxies. For current best values of the distance and approach velocity of Andromeda this leads to a median likelihood estimate of the true mass of the Local Group of $5.27\times 10^{12}\msun$, or $\log M_{LG}/M_\odot = 12.72$, with an interquartile range of $[12.58, 12.83]$ and a 5% to 95% range of $[12.26, 13.01]$. Thus a 95% lower confidence limit on the true mass of the Local Group is $1.81\times 10^{12}\msun$. A timing estimate of the Milky Way's mass based on the large recession velocity observed for the distant satellite Leo I works equally well, although with larger systematic uncertainties. It gives an estimated virial mass for the Milky Way of $2.43 \times 10^{12}\msun$ with a 95% lower confidence limit of $0.80 \times 10^{12}\msun$.

preprint2007arXiv

The statistics of LCDM Halo Concentrations

We use the Millennium Simulation (MS) to study the statistics of LCDM halo concentrations at z = 0. Our results confirm that the average halo concentration declines monotonically with mass; a power-law fits well the concentration-mass relation for over 3 decades in mass, up to the most massive objects to form in a LCDM universe (~ 10^15 h^-1 Msol). This is in clear disagreement with the predictions of the model proposed by Bullock et al. for these rare objects, and agrees better with the original predictions of Navarro, Frenk, & White. The large volume surveyed, together with the unprecedented numerical resolution of the MS, allow us to estimate with confidence the distribution of concentrations and, consequently, the abundance of systems with unusual properties. About one in a hundred cluster haloes (M200 >~ 3x10^14 h^-1 Msol) have concentrations exceeding c200 = 7.5, a result that may be used to interpret the likelihood of unusually strong massive gravitational lenses, such as Abell 1689, in the LCDM cosmogony. A similar fraction (1 in 100) of galaxy-sized haloes (M200 ~ 10^12 h^-1 Msol) have c200 < 4.5, an important constraint on models that attempt to reconcile the rotation curves of low surface-brightness galaxies by appealing to haloes of unexpectedly low concentration. We find that halo concentrations are independent of spin once haloes manifestly out of equilibrium are removed from the sample. Compared to their relaxed brethren, the concentrations of out-of-equilibrium haloes tend to be lower and to have more scatter, while their spins tend to be higher... (continue)

preprint2006arXiv

The large-scale structure of the Universe

Research over the past 25 years has led to the view that the rich tapestry of present-day cosmic structure arose during the first instants of creation, where weak ripples were imposed on the otherwise uniform and rapidly expanding primordial soup. Over 14 billion years of evolution, these ripples have been amplified to enormous proportions by gravitational forces, producing ever-growing concentrations of dark matter in which ordinary gases cool, condense and fragment to make galaxies. This process can be faithfully mimicked in large computer simulations, and tested by a variety of observations that probe the history of the Universe starting from just 400,000 yr after the Big Bang.

preprint2005arXiv

Quantitative Morphology of Galaxies from the SDSS I: Luminosity in Bulges and Disks

In the first paper of this series we use the publicly available code Gim2D to model the r- and i-band images of all galaxies in a magnitude-limited sample of roughly 1800 morphologically classified galaxies taken from the Sloan Digital Sky Survey. The model is a concentric superposition of two components, each with elliptical isophotes with constant flattening and position angle. The disk luminosity profile is assumed exponential, while the bulge is assumed to have a de Vaucouleurs or a Sersic profile. We find that the parameters returned by Gim2D depend little on the waveband or bulge profile used; their formal uncertainties are usually small. Nevertheless, for bright galaxies the measured distribution of b/a, the apparent disk flattening, deviates strongly from the expected uniform distribution, showing that the `disk' identified by the code frequently corresponds to an intrinsically 3-dimensional structure rather than to a true thin disk. We correct approximately for this systematic problem using the observed statistics of the b/a distribution and estimate, as a function of absolute magnitude, the mean fractions of galaxy light in disks and in `pure bulge' systems (those with no detectable disk). For the brightest galaxies the disk light fraction is about 10% and about 80% are `pure bulge' systems. For faint galaxies most of the light is in disks and we do not detect a `pure bulge' population. Averaging over the galaxy population as a whole, we find that 54 \pm 2% of the local cosmic luminosity density at both r and i comes from disks and 32 \pm 2% from `pure bulge' systems. The remaining 14 \pm 2% comes from bulges in galaxies with detectable disks.

preprint2005arXiv

The age dependence of halo clustering

We use a very large simulation of the concordance LCDM cosmogony to study the clustering of dark matter haloes. For haloes less massive than about 1e13Msun/h, the amplitude of the two-point correlation function on large scales depends strongly on halo formation time. Haloes that assembled at high redshift are substantially more clustered than those that assembled more recently. The effect is a smooth function of halo formation time and its amplitude increases with decreasing halo mass. At 1e11 Msun/h the ``oldest'' 10% of haloes are more than 5 times more strongly correlated than the ``youngest'' 10%. This unexpected result is incompatible with the standard excursion set theory for structure growth, and it contradicts a fundamental assumption of the halo occupation distribution models often used to study galaxy clustering, namely that the galaxy content of a halo of given mass is statistically independent of its larger scale environment.

preprint2000arXiv

Weakly Self-Interacting Dark Matter and the Structure of Dark Halos

We study the formation of dark halos in a $Λ$CDM universe under the assumption that Cold Dark Matter particles have a finite cross-section for elastic collisions. We compare evolution when CDM mean free paths are comparable to halo sizes with the collisionless and fluid limits. We show that a few collisions per particle per Hubble time at halo centre can substantially affect the central density profile. Cross-sections an order of magnitude larger produce sufficient relaxation for rich clusters to develop core radii in the range 100-200 $h^{-1}$kpc. The structural evolution of halos is a competition between collisional relaxation caused by individual particle interactions and violent relaxation resulting from the infall and merging processes by which clusters grow. Although our simulations concentrate on systems of cluster size, we can scale our results to address the halo structure expected for dwarf galaxies. We find that collision cross-sections sufficiently large to significantly modify the cores of such galaxies produce cluster cores which are too large and/or too round to be consistent with observation. Thus the simplest model for self-interacting dark matter is unable to improve fits to published dwarf galaxy rotation curves without violating other observational constraints.

preprint1999arXiv

Debris streams in the solar neighbourhood as relicts from the formation of the Milky Way

It is now generally believed that galaxies were built up through gravitational amplification of primordial fluctuations and the subsequent merging of smaller precursor structures. The stars of the structures that assembled to form the Milky Way should now make up much or all of its bulge and halo, in which case one hopes to find "fossil" evidence for those precursor structures in the present distribution of halo stars. Confirmation that this process is continuing came with the discovery of the Sagittarius dwarf galaxy, which is being disrupted by the Milky Way, but direct evidence that this process provided the bulk of the Milky Way's population of old stars has so far been lacking. Here we show that about ten per cent of the metal--poor stars in the halo of the Milky Way, outside the radius of the Sun's orbit, come from a single coherent structure that was disrupted during or soon after the Galaxy's formation. This object had a highly inclined orbit about the Milky Way at a maximum distance of $\sim$ 16 kpc, and it probably resembled the Fornax and Sagittarius dwarf spheroidal galaxies.

preprint1998arXiv

Tidal tails in CDM cosmologies

We study the formation of tidal tails in pairs of merging disk galaxies with structural properties motivated by current theories of cold dark matter (CDM) cosmologies. In a recent study, Dubinski, Mihos & Hernquist (1996) showed that the formation of prominent tidal tails can be strongly suppressed by massive and extended dark haloes. For the large halo-to-disk mass ratio expected in CDM cosmologies their sequence of models failed to produce strong tails like those observed in many well-known pairs of interacting galaxies. In order to test whether this effect can constrain the viability of CDM cosmologies, we construct N-body models of disk galaxies with structural properties derived in analogy to the analytical work of Mo, Mao & White (1998). With a series of self-consistent collisionless simulations of galaxy-galaxy mergers we demonstrate that even the disks of very massive dark haloes have no problems developing long tidal tails, provided the halo spin parameter is large enough. We show that the halo-to-disk mass ratio is a poor indicator for the ability to produce tails. Instead, the relative size of disk and halo, or alternatively, the ratio of circular velocity to local escape speed at the half mass radius of the disk are more useful criteria. This result holds in all CDM cosmologies. The length of tidal tails is thus unlikely to provide useful constraints on such models.

preprint1997arXiv

The Structure and Dynamical Evolution of Dark Matter Halos

(Shortened) We use N-body simulations to investigate the structure and dynamical evolution of dark matter halos in galaxy clusters. Our sample consists of nine massive halos from an EdS universe with scale free power spectrum and n = -1. Halos are resolved by ~20000 particles each, with a dynamical resolution of 20-25 kpc. Large scale tidal fields are included up to L=150 Mpc using background particles. The halo formation process can be characterized by the alternation of two dynamical configurations: a merging phase and a relaxation phase, defined by their signature on the evolution of the total mass and rms velocity. Halos spend on average one 1/3 of their evolution in the merging phase and 2/3 in the relaxation phase. Using this definition, we study the density profiles and their change during the halo history. The average density profiles are fitted by the NFW analytical model with an rms residual of 17% between the virial radius Rv and 0.01 Rv. The Hernquist (1990) profiles fits the same halos with an rms residual of 26%. The trend with mass of the scale radius of these fits is marginally consistent with that found by Cole & Lacey (1996): in comparison our halos are more centrally concentrated, and the relation between scale radius and halo mass is slightly steeper. We find a moderately large scatter in this relation, due both to dynamical evolution within halos and to fluctuations in the halo population. We analyze the dynamical equilibrium of our halos using the Jeans' equation, and find that on average they are approximately in equilibrium within their virial radius. Finally, we find that the projected mass profiles of our simulated halos are in very good agreement with the profiles of three rich galaxy clusters derived from strong and weak gravitational lensing observations.

preprint1996arXiv

Detection of Strong Evolution in the Population of Early-Type Galaxies

The standard picture holds that giant elliptical galaxies formed in a single burst at high redshift. Aging of their stellar populations subsequently caused them to fade and become redder. The Canada-France Redshift Survey provides a sample of about 125 galaxies with the luminosities and colours of passively evolving giant ellipticals and with 0.1 < z < 1. This sample is inconsistent with the standard evolutionary picture with better than 99.9% confidence. The standard Schmidt test gives <V/V_{max}> = 0.398 when restricted to objects with no detected star formation, and <V/V_{max}> = 0.410 when objects with emission lines are also included. A smaller sample of early-type galaxies selected from tha Hawaii Deep Survey gives equally significant results. With increasing redshift a larger and larger fraction of the nearby elliptical and S0 population must drop out of of the sample, either because the galaxies are no longer single units or because star formation alters their colours. If the remaining fraction is modelled as F=(1+z)^{-γ}, the data imply gamma = 1.5 +- 0.3. At z=1 only about one third of bright E and S0 galaxies were already assembled and had the colours of old passively evolving systems. We discuss the sensitivity of these results to the incompleteness corrections and stellar population models we have adopted. We conclude that neither is uncertain enough to reconcile the observations with the standard picture. Hierarchical galaxy formation models suggest that both merging and recent star formation play a role in the strong evolution we have detected.

preprint1996arXiv

Two--body heating in numerical galaxy formation experiments

We show that discreteness effects related to classical two-body relaxation produce spurious heating of the gaseous component in numerical simulations of galaxy formation. A simple analytic model demonstrates that this artificial heating will dominate radiative cooling in any simulation where the mass of an individual dark matter particle exceeds a certain critical value. This maximum mass depends only on the cooling function of the gas, on the fraction of the material in gaseous form, and (weakly) on typical temperatures in the gas. It is comparable to, or smaller than, the dark matter particle masses employed in most published simulations of cosmological hydrodynamics and galaxy formation. Any simulation which violates this constraint will be unable to follow cooling flows, although catastrophic cooling of gas may still occur in regions with very short cooling times. We use a series of N--body/smoothed particle hydrodynamics simulations to explore this effect. In simulations which neglect radiative cooling, two--body heating causes a gradual expansion of the gas component. When radiative effects are included, we find that gas cooling is almost completely suppressed for dark matter particle masses above our limit. Although our test simulations use smoothed particle hydrodynamics, similar effects, and a similar critical mass, are expected in any simulation where the dark matter is represented by discrete particles.

preprint1994arXiv

Simulations of X-ray Clusters

We present simulations of the formation and evolution of galaxy clusters in the Cold Dark Matter cosmogony. Clusters with a wide range of mass were selected from previous N-body models, and were resimulated at higher resolution using a combined N-body/Smooth Particle Hydrodynamics code. The effects of radiative cooling on the gas are neglected. While many present-day clusters are predicted to be undergoing mergers, the density profiles of those that are approximately in equilibrium are all very similar, both for the gas and for the dark matter. These profiles show no sign of a uniform density core and steepen gradually from the centre outwards. The standard $β$-model is a reasonable fit over most of the radius range observable in real clusters. However, the value obtained for the slope parameter $β_f$ increases with the outermost radius of the fit. Temperature profiles of different simulated clusters are also similar. Typically the temperature is almost uniform in the regions which emit most of the X-ray flux but drops at larger radii. The gas temperature and dark matter velocity dispersion in equilibrium clusters give values of $β_T\equiv μm_pσ_{DM}^2/kT$ which are consistent with unity provided an X-ray emission-weighted temperature is used. Larger values of $β_T$ are found in merging objects where there is a transient boost in the velocity dispersion of the system. Thus $β_T >1$ may be an observational indicator of merging in real clusters. The similar structure of clusters of differing mass results in scaling relations between the X-ray and dynamical properties of clusters identified at

preprint1994arXiv

The Assembly of Galaxies in a Hierarchically Clustering Universe

We study the formation of galaxies by using $N$-body/hydrodynamics simulations to investigate how baryons collect at the centre of dark matter halos. We treat the dark matter as a collisionless fluid and the baryons as an ideal gas. We include the effects of gravity, pressure gradients, hydrodynamical shocks, and radiative energy losses, but we neglect star formation. Our initial conditions assume a flat universe dominated by cold dark matter with a mean baryon abundance of 10\% by mass. Typical halos form through the merging of a few smaller systems which had themselves formed in a similar manner at higher redshift. The gas collects at the bottom of dark matter potential wells as soon as these are properly resolved by our simulations. There it settles into cold, tightly bound disks, and it remains cold during subsequent evolution. As their halos coalesce, these disks merge on a timescale that is consistent with dynamical friction estimates based on their {\it total} (gas + surrounding dark matter) mass. Both the merger rates of the disks and their mass spectrum are in remarkably good agreement with recent analytic models that describe the evolution of {\it dark halos} in a hierarchical universe. This very simple model of galaxy formation suffers from serious shortcomings. It predicts that most baryons should be locked up in galaxies, whereas in the real universe most baryons are thought to lie outside visible galaxies. In addition, it predicts the specific angular momentum of a disk to be only about 20\% that of its surrounding halo, corresponding to a radius smaller than that of observed spiral galaxy disks.

Simon D. M. White

What is connected

Connect this record

See the researcher in context

Building this map preview

104 published item(s)

Dark matter annihilation and the Galactic Centre Excess

Disc instability and bar formation: view from the IllustrisTNG simulations

Comparing galaxy formation in the L-GALAXIES semi-analytical model and the IllustrisTNG simulations

Measuring the Tidal Response of Structure Formation: Anisotropic Separate Universe Simulations using TreePM

Simulating the Complexity of the Dark Matter Sheet I: Numerical Algorithms

The globular cluster system of the Auriga simulations

The Tessellation-Level-Tree: characterising the nested hierarchy of density peaks and their spatial distribution in cosmological N-body simulations

Universal structure of dark matter haloes over a mass range of 20 orders of magnitude

A chronicle of galaxy mass assembly in the EAGLE simulation

A fully cosmological model of a Monoceros-like ring

Gaps in globular cluster streams: giant molecular clouds can cause them too

Quiescence correlates strongly with directly-measured black hole mass in central galaxies

The Diversity of Growth Histories of Milky Way-mass Galaxies

The EAGLE simulations of galaxy formation: public release of halo and galaxy catalogues

The galaxy correlation function as a constraint on galaxy formation physics

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

Surface photometry of BCGs and intracluster stars in Lambda-CDM

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

The EAGLE simulations of galaxy formation: calibration of subgrid physics and model variations

The unexpected diversity of dwarf galaxy rotation curves

Unifying X-ray Scaling Relations from Galaxies to Clusters

Elemental Abundances in Milky Way-like Galaxies from a Hierarchical Galaxy Formation Model

Lensing Model of MACS J1149.5+2223 - I. Cluster mass reconstruction

Local Group galaxies emerge from the dark

Matched filter optimization of kSZ measurements with a reconstructed cosmological flow field

Satellite abundances around bright isolated galaxies II: radial distribution and environmental effects

The diverse formation histories of simulated disc galaxies

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

The re-distribution of matter in the cores of galaxy clusters

The Tenth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-III Apache Point Observatory Galactic Evolution Experiment

Why stellar feedback promotes disc formation in simulated galaxies

Estimating gas masses and dust-to-gas ratios from optical spectroscopy

Exploring the non-linear density field in the Millennium simulations with tessellations - I. The probability distribution function

Galactic accretion and the outer structure of galaxies in the CDM model

Galaxy formation on the largest scales: The impact of astrophysics on the BAO peak

Modelling Element Abundances in Semi-analytic Models of Galaxy Formation

Streams in the Aquarius stellar haloes

The Growth in Size and Mass of Cluster Galaxies since z=2

The Mass Profile and Accretion History of Cold Dark Matter Halos

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

The satellites of the Milky Way - Insights from semi-analytic modelling in a LambdaCDM cosmology

Dark matter and cosmic structure

Evidence for extended gamma-ray emission from galaxy clusters

Galactic star formation and accretion histories from matching galaxies to dark matter haloes

Idealized models for galactic disk formation and evolution in 'realistic' LambdaCDM haloes

On the Tidal Dependence of Galaxy Properties

Satellite abundances around bright isolated galaxies

Satellites and Haloes of Dwarf Galaxies

Shallow Dark Matter Cusps in Galaxy Clusters

Stellar Disks in Aquarius Dark Matter Haloes

The Effect of Particle Noise in N-body Simulations of Gravitational Lensing

The effects of alignment and ellipticity on the clustering of galaxies

The effects of halo alignment and shape on the clustering of galaxies

The Origin of Disks and Spheroids in Simulated Galaxies

Autocorrelations of stellar light and mass at z~0 and ~1: From SDSS to DEEP2

Chemical signatures of formation processes in the stellar populations of simulated galaxies

Clues to the "Magellanic Galaxy" from Cosmological Simulations

SDSS-III: Massive Spectroscopic Surveys of the Distant Universe, the Milky Way Galaxy, and Extra-Solar Planetary Systems

Streams and caustics: the fine-grained structure of LCDM haloes

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

The Eighth Data Release of the Sloan Digital Sky Survey: First Data from SDSS-III

The Shape of Dark Matter Haloes in the Aquarius Simulations: Evolution and Memory

An observer's view of simulated galaxies: disc-to-total ratios, bars, and (pseudo-)bulges

Autocorrelations of stellar light and mass in the low-redshift Universe

Bound and unbound substructures in Galaxy-scale Dark Matter haloes

Dark matter response to galaxy formation

Kinematics of Milky Way Satellites in a Lambda Cold Dark Matter Universe

Nonspherical similarity solutions for dark halo formation

Relic density and CMB constraints on dark matter annihilation with Sommerfeld enhancement

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

The evolution of the density of galaxy clusters and groups: denser environments at higher redshifts

There's no place like home? Statistics of Milky Way-mass dark matter halos

What is the (Dark) Matter with Dwarf Galaxies?

Abundances, masses, and weak-lensing mass profiles of galaxy clusters as a function of richness and luminosity in LambdaCDM cosmologies