Search

Many distant objects can only be detected, or become more scientifically valuable, if they have been highly magnified by strong gravitational lensing. We use EAGLE and BAHAMAS, two recent cosmological hydrodynamical simulations, to predict the probability distribution for both the lens mass and lens redshift when point sources are highly magnified by gravitational lensing. For sources at a redshift of two, we find the distribution of lens redshifts to be broad, peaking at z=0.6. The contribution of different lens masses is also fairly broad, with most high-magnification lensing due to lenses with halo masses between 10^12 and 10^14 solar masses. Lower mass haloes are inefficient lenses, while more massive haloes are rare. We find that a simple model in which all haloes have singular isothermal sphere density profiles can approximately reproduce the simulation predictions, although such a model over-predicts the importance of haloes with mass, 15 pages, 5 figures, updated to match MNRAS version

We perform simulations of giant impacts onto the young Uranus using smoothed particle hydrodynamics (SPH) with over 100 million particles. This 100--1000$\times$ improvement in particle number reveals that simulations with below 10^7 particles fail to converge on even bulk properties like the post-impact rotation period, or on the detailed erosion of the atmosphere. Higher resolutions appear to determine these large-scale results reliably, but even 10^8 particles may not be sufficient to study the detailed composition of the debris -- finding that almost an order of magnitude more rock is ejected beyond the Roche radius than with 10^5 particles. We present two software developments that enable this increase in the feasible number of particles. First, we present an algorithm to place any number of particles in a spherical shell such that they all have an SPH density within 1% of the desired value. Particles in model planets built from these nested shells have a root-mean-squared velocity below 1% of the escape speed, which avoids the need for long precursor simulations to produce relaxed initial conditions. Second, we develop the hydrodynamics code SWIFT for planetary simulations. SWIFT uses task-based parallelism and other modern algorithmic approaches to take full advantage of contemporary supercomputer architectures. Both the particle placement code and SWIFT are publicly released., Fixed typo in Appx B equation and updated urls

We establish the feasibility of measuring the neutron lifetime via an alternative, space-based class of methods, which use neutrons generated by galactic cosmic ray spallation of planets' surfaces and atmospheres. Free neutrons decay via the weak interaction with a mean lifetime of around 880 s. This lifetime constrains the unitarity of the CKM matrix and is a key parameter for studies of Big-Bang nucleosynthesis. However, current laboratory measurements, using two independent approaches, differ by over 4$\sigma$. Using data acquired in 2007 and 2008 during flybys of Venus and Mercury by NASA's MESSENGER spacecraft, which was not designed to make this measurement, we estimate the neutron lifetime to be $780\pm60_\textrm{stat}\pm70_\textrm{syst}$ s, thereby demonstrating the viability of this new approach.

A set of key observations over the Cassini spacecraft's tenure has constrained Saturn's rings' age to be less than a few 100 Myr effectively ruling out currently accepted ring origin scenarios, all of which require that the rings are ancient or primordial. We propose a new scenario motivated from evidence of a comparably recent dynamical instability ~100 Myr ago which would have led to collisions between Saturn's pre-existing mid-size icy moons, opening the door to possible ring formation during that epoch. Successfully testing this scenario requires better understanding of collisional outcomes. Toward that end, we introduce a new suite of simulations modeling impacts between Saturn's icy moons using the next generation smoothed hydrodynamical and gravity code SWIFT. The unprecedented spatial resolution achieved in these simulations (108.5 particles within the simulation box) allows us to depict the myriad of gravitationally bound objects formed during icy moon collisions which may afterwards evolve both thermally and dynamically to re-accrete or collide with other bodies. Our unprecedented high resolution further allows us to determine a size distribution of fragments which can be used to inform crater impact distributions on newly accreted or remaining moons.

Recent claims of observational evidence for self-interacting dark matter (SIDM) have relied on a semi-analytic method for predicting the density profiles of galaxies and galaxy clusters containing SIDM. We present a thorough description of this method, known as isothermal Jeans modelling, and then test it with a large ensemble of haloes taken from cosmological simulations. Our simulations were run with cold and collisionless dark matter (CDM) as well as two different SIDM models, all with dark matter only variants as well as versions including baryons and relevant galaxy formation physics. Using a mix of different box sizes and resolutions, we study haloes with masses ranging from 3e10 to 3e15 Msun. Overall, we find that the isothermal Jeans model provides as accurate a description of simulated SIDM density profiles as the Navarro-Frenk-White profile does of CDM halos. We can use the model predictions, compared with the simulated density profiles, to determine the input DM-DM scattering cross-sections used to run the simulations. This works especially well for large cross-sections, while with CDM our results tend to favour non-zero (albeit fairly small) cross-sections, driven by a bias against small cross-sections inherent to our adopted method of sampling the model parameter space. The model works across the whole halo mass range we study, although including baryons leads to DM profiles of intermediate-mass (10^12 - 10^13 Msun) haloes that do not depend strongly on the SIDM cross-section. The tightest constraints will therefore come from lower and higher mass haloes: dwarf galaxies and galaxy clusters., Comment: 27 pages, 19 figures, accepted in MNRAS, updated to match accepted version

We simulate the hypothesised collision between the proto-Earth and a Mars-sized impactor that created the Moon. Amongst the resulting debris disk in some impacts, we find a self-gravitating clump of material. It is roughly the mass of the Moon, contains $\sim1\%$ iron like the Moon, and has its internal composition resolved for the first time. The clump contains mainly impactor material near its core but becomes increasingly enriched in proto-Earth material near its surface. A graduated composition has recently been measured in the oxygen isotope ratios of Apollo samples, suggesting incomplete mixing between proto-Earth and impactor material that formed the Moon. However, the formation of the Moon-sized clump depends sensitively on the spin of the impactor. To explore this, we develop a fast method to construct models of multi-layered, rotating bodies and their conversion into initial conditions for smoothed particle hydrodynamical (SPH) simulations. We use our publicly available code to calculate density and pressure profiles in hydrostatic equilibrium, then generate configurations of over a billion particles with SPH densities within $1\%$ of the desired values. This algorithm runs in a few minutes on a desktop computer, for $10^7$ particles, and allows direct control over the properties of the spinning body. In comparison, relaxation or spin-up techniques that take hours on a supercomputer before the structure of the rotating body is even known. Collisions that differ only in the impactor's initial spin reveal a wide variety of outcomes: a merger, a grazing hit-and-run, or the creation of an orbiting proto-Moon., Comment: Submitted to MNRAS. 11 pages, 9 figures

We use data from the Lunar Prospector Neutron Spectrometer to make the second space-based measurement of the free neutron lifetime finding $\tau_n=887 \pm 14_\text{stat}{\:^{+7}_{-3\:\text{syst}}}$ s, which is within 1$\sigma$ of the accepted value. This measurement expands the range of planetary bodies where the neutron lifetime has been quantified from space, and by extending the modeling to account for non-uniform elemental composition, we mitigated a significant source of systematic uncertainty on the previous space-based lifetime measurement. This modeling moves space-based neutron lifetime measurement towards the ultimate goal of reducing the magnitude of the systematics on a future space-measurement to the level of those seen in laboratory-based experiments., Comment: 10 pages, 9 figures, 2 tables, accepted in PRC; revised modeling and systematic error analysis following referee's suggestion, reduced uncertainties

We examine the mechanisms by which atmosphere can be eroded by giant impacts onto Earth-like planets with thin atmospheres, using 3D smoothed particle hydrodynamics simulations with sufficient resolution to directly model the fate of low-mass atmospheres. We present a simple scaling law to estimate the fraction lost for any impact angle and speed in this regime. In the canonical Moon-forming impact, only around 10% of the atmosphere would have been lost from the immediate effects of the collision. There is a gradual transition from removing almost none to almost all of the atmosphere for a grazing impact as it becomes more head-on or increases in speed, including complex, non-monotonic behaviour at low impact angles. In contrast, for head-on impacts, a slightly greater speed can suddenly remove much more atmosphere. Our results broadly agree with the application of 1D models of local atmosphere loss to the ground speeds measured directly from our simulations. However, previous analytical models of shock-wave propagation from an idealised point-mass impact significantly underestimate the ground speeds and hence the total erosion. The strong dependence on impact angle and the interplay of multiple non-linear and asymmetrical loss mechanisms highlight the need for 3D simulations in order to make realistic predictions., Comment: 19 pages, 13 figures, published in ApJ

We present improved resolution maps of the Lunar Prospector Neutron Spectrometer thermal, epithermal and fast neutron data and Gamma-Ray Spectrometer Th-line fluxes via global application of pixon image reconstruction techniques. With the use of mock data sets, we show that the pixon image reconstruction method compares favorably with other methods that have been used in planetary neutron and gamma-ray spectroscopy. The improved thermal neutron maps are able to clearly distinguish variations in composition across the lunar surface, including within the lunar basins of Hertzsprung and Schrodinger. The improvement in resolution reveals a correlation between albedo and thermal neutron flux within the craters. The consequent increase in dynamic range confirms that Hertzsprung basin contains one of the most anorthositic parts of the lunar crust, including nearly pure anorthite over a region tens of km in diameter. At Orientale, the improvement in spatial resolution of the epithermal neutron data show that there is a mismatch between measures of regolith maturity that sample the surface and those that probe the near-subsurface, which suggests a complex layering scenario., Comment: 17 pages, 12 figures, submitted to JGR-Planets

We present a map of the near subsurface hydrogen distribution on Mars, based on epithermal neutron data from the Mars Odyssey Neutron Spectrometer. The map's spatial resolution is approximately improved two-fold via a new form of the pixon image reconstruction technique. We discover hydrogen-rich mineralogy far from the poles, including ~10 wt. % water equivalent hydrogen (WEH) on the flanks of the Tharsis Montes and greater than 40 wt. % WEH at the Medusae Fossae Formation (MFF). The high WEH abundance at the MFF implies the presence of bulk water ice. This supports the hypothesis of recent periods of high orbital obliquity during which water ice was stable on the surface. We find the young undivided channel system material in southern Elysium Planitia to be distinct from its surroundings and exceptionally dry; there is no evidence of hydration at the location in Elysium Planitia suggested to contain a buried water ice sea. Finally, we find that the sites of recurring slope lineae (RSL) do not correlate with subsurface hydration. This implies that RSL are not fed by large, near-subsurface aquifers, but are instead the result of either small (less than 120 km diameter) aquifers, deliquescence of perchlorate and chlorate salts or dry, granular flows., Accepted for publication in Icarus, 16 pages, 12 figures

We perform the first tests of various proposed explanations for observed features of the Moon's argon exosphere, including models of the following: spatially varying surface interactions; a source that reflects the lunar near-surface potassium distribution; and temporally varying cold trap areas. Measurements from the Lunar Atmosphere and Dust Environment Explorer (LADEE) and the Lunar Atmosphere Composition Experiment (LACE) are used to test whether these models can reproduce the data. The spatially varying surface interactions hypothesized in previous work cannot reproduce the persistent argon enhancement observed over the western maria. They also fail to match the observed local time of the near-sunrise peak in argon density, which is the same for the highland and mare regions and is well reproduced by simple surface interactions with a ubiquitous desorption energy of 28 kJ mol−1. A localized source can explain the observations, with a trade-off between an unexpectedly localized source or an unexpectedly brief lifetime of argon atoms in the exosphere. To match the observations, a point-like source requires source and loss rates of ∼1.9 × 1021 atoms s−1. A more diffuse source, weighted by the near-surface potassium, requires much higher rates of ∼1.1 × 1022 atoms s−1, corresponding to a mean lifetime of just 1.4 lunar days. We do not address the mechanism for producing a localized source, but demonstrate that this appears to be the only model that can reproduce the observations. Large, seasonally varying cold traps could explain the long-term fluctuation in the global argon density observed by LADEE, but not that by LACE.

An estimate is made of the thickness of the radar-bright deposits in craters near to the north pole of Mercury. To construct an objective set of craters for this measurement, an automated crater finding algorithm is developed and applied to a digital elevation model based on data from the Mercury Laser Altimeter on board the MESSENGER spacecraft. This produces a catalogue of 663 craters with diameters exceeding 4 km, northwards of latitude +55 degrees. A subset of 12 larger, well-sampled and fresh polar craters are selected to search for correlations between topography and radar same-sense backscatter cross-section. It is found that the typical excess height associated with the radar-bright regions within these fresh polar craters is (50+/-35)m. This puts an approximate upper limit on the total polar water ice deposits on Mercury of 3e15 kg., Comment: 10 pages, 6 figures, submitted to Icarus

We present BAHAMAS-SIDM, the first large-volume, (400/h Mpc)^3, cosmological simulations including both self-interacting dark matter (SIDM) and baryonic physics. These simulations are important for two primary reasons: 1) they include the effects of baryons on the dark matter distribution 2) the baryon particles can be used to make mock observables that can be compared directly with observations. As is well known, SIDM haloes are systematically less dense in their centres, and rounder, than CDM haloes. Here we find that that these changes are not reflected in the distribution of gas or stars within galaxy clusters, or in their X-ray luminosities. However, gravitational lensing observables can discriminate between DM models, and we present a menu of tests that future surveys could use to measure the SIDM interaction strength. We ray-trace our simulated galaxy clusters to produce strong lensing maps. Including baryons boosts the lensing strength of clusters that produce no critical curves in SIDM-only simulations. Comparing the Einstein radii of our simulated clusters with those observed in the CLASH survey, we find that at velocities around 1000 km/s an SIDM cross-section of sigma/m > 1 cm^2/g is likely incompatible with observed cluster lensing., 19 pages, 13 figures, updated to match MNRAS version

We present a new scaling law to predict the loss of atmosphere from planetary collisions for any speed, angle, impactor mass, target mass, and body compositions, in the regime of giant impacts onto broadly terrestrial planets with relatively thin atmospheres. To this end, we examine the erosion caused by a wide range of impacts, using 3D smoothed particle hydrodynamics simulations with sufficiently high resolution to directly model the fate of low-mass atmospheres around 1% of the target's mass. Different collision scenarios lead to extremely different behaviours and consequences for the planets. In spite of this complexity, the fraction of lost atmosphere is fitted well by a power law. Scaling is independent of the system mass for a constant impactor mass ratio. Slow atmosphere-hosting impactors can also deliver a significant mass of atmosphere, but always accompanied by larger proportions of their mantle and core. Different Moon-forming impact hypotheses suggest that around 10 to 60% of a primordial atmosphere could have been removed directly, depending on the scenario. We find no evident departure from the scaling trends at the extremes of the parameters explored. The scaling law can be incorporated readily into models of planet formation., Published in ApJL. 12 pages, 6 figures

In Miniature Radio Frequency (Mini‐RF) radar images, anomalous craters are those having a high circular polarization ratio (CPR) in their interior but not exterior to their rims. Previous studies found that most CPR‐anomalous craters contain permanently shadowed regions (PSRs) and that their population is overabundant in the polar regions. However, there is considerable controversy in the interpretation of these signals: both water ice deposits and rocks/surface roughness have been proposed as the source of the elevated CPR values. To resolve this controversy, we have systematically analyzed >4000 impact craters with diameters between 2.5 and 24 km in the Mini‐RF radar image and Diviner rock abundance (RA) map. We first constructed two controlled orthorectified global mosaics using 6818 tracks of Mini‐RF raw data, and then analyzed the correlations between radar CPR and surface slope, rock abundance, and depth/diameter ratios of impact craters. Our results show that CPR‐anomalous craters are distributed relatively uniformly across the lunar surface, with no apparent difference in CPR between the polar, potentially icy, and non‐polar, not icy, craters. Most CPR‐anomalous craters are relatively young with a large depth/diameter ratio, and they actually represent an intermediate stage of crater evolution. Comparison with a two‐component radar scattering model suggests that rocks and surface roughness are major contributors to the observed CPR values. Using craters of 4.7–22 km in diameter with known ages, we find that craters spend up to 120 Ma with a high exterior RA, and ∼3 Ga in the CPR‐anomalous phase.

We perform a suite of smoothed particle hydrodynamics simulations to investigate in detail the results of a giant impact on the young Uranus. We study the internal structure, rotation rate, and atmospheric retention of the post-impact planet, as well as the composition of material ejected into orbit. Most of the material from the impactor's rocky core falls in to the core of the target. However, for higher angular momentum impacts, significant amounts become embedded anisotropically as lumps in the ice layer. Furthermore, most of the impactor's ice and energy is deposited in a hot, high-entropy shell at a radius of ~3 Earth radii. This could explain Uranus' observed lack of heat flow from the interior and be relevant for understanding its asymmetric magnetic field. We verify the results from the single previous study of lower resolution simulations that an impactor with a mass of at least 2 Earth masses can produce sufficiently rapid rotation in the post-impact Uranus for a range of angular momenta. At least 90% of the atmosphere remains bound to the final planet after the collision, but over half can be ejected beyond the Roche radius by a 2 or 3 Earth mass impactor. This atmospheric erosion peaks for intermediate impactor angular momenta (~3*10^36 kg m^2 s^-1). Rock is more efficiently placed into orbit and made available for satellite formation by 2 Earth mass impactors than 3 Earth mass ones, because it requires tidal disruption that is suppressed by the more massive impactors., Comment: Published in ApJ. 14 pages, 11 figures

The variation of remotely sensed neutron count rates is measured as a function of cratercentric distance using data from the Lunar Prospector Neutron Spectrometer. The count rate, stacked over many craters, peaks over the crater centre, has a minimum near the crater rim, and at larger distances, it increases to a mean value that is up to 1% lower than the mean count rate observed over the crater. A simple model is presented, based upon an analytical topographical profile for the stacked craters fitted to data from the Lunar Orbiter Laser Altimeter (LOLA). The effect of topography coupled with neutron beaming from the surface largely reproduces the observed count rate profiles. However, a model that better fits the observations can be found by including the additional freedom to increase the neutron emissivity of the crater area by ~0.35% relative to the unperturbed surface. It is unclear what might give rise to this effect, but it may relate to additional surface roughness in the vicinities of craters. The amplitude of the crater-related signal in the neutron count rate is small, but not too small to demand consideration when inferring water-equivalent hydrogen (WEH) weight percentages in polar permanently shaded regions (PSRs). If the crater-wide count rate excess is concentrated into a much smaller PSR, then it can lead to a significantly biased inferred WEH weight percentage. For instance, it may increase the inferred WEH for Cabeus crater at the Moon's South Pole from ~1% to ~4%.

We reconstruct the abundance of thorium near the Compton-Belkovich Volcanic Complex on the Moon, using data from the Lunar Prospector Gamma Ray Spectrometer. We enhance the resolution via a pixon image reconstruction technique, and find that the thorium is distributed over a larger ($40 \mathrm{km}\times 75$ km) area than the ($25 \mathrm{km}\times 35$ km) high albedo region normally associated with Compton-Belkovich. Our reconstructions show that inside this region, the thorium concentration is $14\!-\!26$ ppm. We also find additional thorium, spread up to $300$ km eastward of the complex at $\sim\!2$ ppm. The thorium must have been deposited during the formation of the volcanic complex, because subsequent lateral transport mechanisms, such as small impacts, are unable to move sufficient material. The morphology of the feature is consistent with pyroclastic dispersal and we conclude that the present distribution of thorium was likely created by the explosive eruption of silicic magma., v2: 13 pages, 8 figures, accepted version, to be published in JGR-Planets; v1: 11 pages, 6 figures, submitted to JGR-Planets

We perform numerical simulations of the merging galaxy cluster 1E 0657-56 (the Bullet Cluster), including the effects of elastic dark matter scattering. In a similar manner to the stripping of gas by ram pressure, dark matter self-interactions would transfer momentum between the two galaxy cluster dark matter haloes, causing them to lag behind the collisionless galaxies. The absence of an observed separation between the dark matter and stellar components in the Bullet Cluster has been used to place upper limits on the cross-section for dark matter scattering. We emphasise the importance of analysing simulations in an observationally-motivated manner, finding that the way in which the positions of the various components are measured can have a larger impact on derived constraints on dark matter's self-interaction cross-section than reasonable changes to the initial conditions for the merger. In particular, we find that the methods used in previous studies to place some of the tightest constraints on this cross-section do not reflect what is done observationally, and overstate the Bullet Cluster's ability to constrain the particle properties of dark matter. We introduce the first simulations of the Bullet Cluster including both self-interacting dark matter and gas. We find that as the gas is stripped it introduces radially-dependent asymmetries into the stellar and dark matter distributions. As the techniques used to determine the positions of the dark matter and galaxies are sensitive to different radial scales, these asymmetries can lead to erroneously measured offsets between dark matter and galaxies even when they are spatially coincident., 21 pages, 16 figures, v2 matches accepted MNRAS version

A detailed comparison is made of results from the Lunar Prospector Neutron Spectrometer (LPNS) and the Lunar Exploration Neutron Detector Collimated Sensors for EpiThermal Neutrons (LEND CSETN). Using the autocorrelation function and power spectrum of the polar count rate maps produced by these experiments, it is shown that the LEND CSETN has a footprint that is at least as big as would be expected for an omni-directional detector at an orbital altitude of 50 km. The collimated flux into the field of view of the collimator is negligible. Arguments put forward asserting otherwise are considered and found wanting for various reasons. The maps of lunar polar hydrogen with the highest contrast, i.e. spatial resolution, are those resulting from pixon image reconstructions of the LPNS data. These typically provide weight percentages of water equivalent hydrogen that are accurate to 30% within the polar craters., 12 pages, 13 figures

Well-motivated elementary particle candidates for the dark matter, such as the sterile neutrino, behave as warm dark matter (WDM).For particle masses of order a keV, free streaming produces a cutoff in the linear fluctuation power spectrum at a scale corresponding to dwarf galaxies. We investigate the abundance and structure of WDM haloes and subhaloes on these scales using high resolution cosmological N-body simulations of galactic haloes of mass similar to the Milky Way's. On scales larger than the free-streaming cutoff, the initial conditions have the same power spectrum and phases as one of the cold dark matter (CDM) haloes previously simulated by Springel et al as part of the Virgo consortium Aquarius project. We have simulated four haloes with WDM particle masses in the range 1.4-2.3keV and, for one case, we have carried out further simulations at varying resolution. N-body simulations in which the power spectrum cutoff is resolved are known to undergo artificial fragmentation in filaments producing spurious clumps which, for small masses (, 20 pages, 21 Figures, Published in MNRAS, minor updates to particle mass calculations: 1.4keV -> 1.5keV

We investigate the luminosity functions (LFs) and projected number density profiles of galactic satellites around isolated primaries of different luminosities. We measure these quantities for model satellites placed into the Millennium and Millennium II dark matter simulations by the GALFORM semi-analytic galaxy formation model for different bins of primary galaxy magnitude and we investigate their dependence on satellite luminosity. We compare our model predictions to the data of Guo et al. from the Sloan Digital Sky Survey Data Release 8 (SDSS DR8). First, we use a mock light-cone catalogue to verify that the method we used to count satellites in the SDSS DR8 is unbiased. We find that the radial distributions of model satellites are similar to those around comparable primary galaxies in the SDSS DR8, with only slight differences at low luminosities and small projected radii. However, when splitting the satellites by colour, the model and SDSS satellite systems no longer resemble one another, with many red model satellites, in contrast to the dominant blue fraction at similar luminosity in SDSS. The few model blue satellites are also significantly less centrally concentrated in the halo of their stacked primary than their SDSS counterparts. The implications of this result for the GALFORM model are discussed., 11 pages, 10 figures, submitted to MNRAS

Monthly Notices of the Royal Astronomical Society, 352 (1), ISSN:0035-8711, ISSN:1365-2966, ISSN:1365-8711

We investigate how self-interacting dark matter (SIDM) with anisotropic scattering affects the evolution of isolated dark matter haloes as well as systems with two colliding haloes. For isolated haloes, we find that the evolution can be adequately captured by treating the scattering as isotropic, as long as the isotropic cross-section is appropriately matched to the underlying anisotropic model. We find that this matching should not be done using the momentum transfer cross-section, as has been done previously. Matching should instead be performed via a modified momentum transfer cross-section that takes into account that dark matter particles can be relabelled after they scatter, without altering the dynamics. However, using cross-sections that are matched to give the same behaviour in isolated haloes, we find that treating dark matter scattering as isotropic under-predicts the effects of anisotropic dark matter scattering when haloes collide. In particular, the DM-galaxy offset induced by SIDM in colliding galaxy clusters is larger when we simulate the underlying particle model, than if we use a matched isotropic model. On the other hand, well motivated particle models with anisotropic scattering typically have cross-sections with a strong velocity dependence, and we discover a previously unrecognised effect that suppresses DM-galaxy offsets in colliding clusters making it hard for these systems to provide competitive constraints on such particle models., Comment: 12 pages, 8 figures, v2 matches accepted MNRAS version

Cold dark matter models predict the existence of a large number of substructures within dark matter halos. If the cold dark matter consists of weakly interacting massive particles, their annihilation within these substructures could lead to diffuse GeV emission that would dominate over the annihilation signal of the host halo. In this work we search for GeV emission from three nearby galaxy clusters: Coma, Virgo and Fornax. We first remove known extragalactic and galactic diffuse gamma-ray backgrounds and point sources from the Fermi 2-year catalog and find a significant residual diffuse emission in all three clusters. We then investigate whether this emission is due to (i) unresolved point sources; (ii) dark matter annihilation; or (iii) cosmic rays (CR). Using 45 months of Fermi-LAT data we detect several new point sources (not present in the Fermi 2-year point source catalogue) which contaminate the signal previously analyzed by Han et al.(arXiv:1201.1003). Including these and accounting for the effects of undetected point sources, we find no significant detection of extended emission from the three clusters studied. Instead, we determine upper limits on emission due to dark matter annihilation and cosmic rays. For Fornax and Virgo the limits on CR emission are consistent with theoretical models, but for Coma the upper limit is a factor of 2 below the theoretical expectation. Allowing for systematic uncertainties associated with the treatment of CR, the upper limits on the cross section for dark matter annihilation from our clusters are more stringent than those from analyses of dwarf galaxies in the Milky Way. We rule out the thermal cross section for supersymmetric dark matter particles for masses as large as 100 GeV (depending on the annihilation channel).

We study the spatial distribution of satellite galaxies around isolated primaries using the Sloan Digital Sky Survey (SDSS) spectroscopic and photometric galaxy catalogues. We select isolated primaries from the spectroscopic sample and search for potential satellites in the much deeper photometric sample. For specific luminosity primaries we obtain robust statistical results by stacking as many as ∼50 000 galaxy systems. We find no evidence for any anistropy in the satellite galaxy distribution relative to the major axes of the primaries. We derive accurate projected number density profiles of satellites down to 4 mag fainter than their primaries. We find that the normalized satellite profiles generally have a universal form and can be well fitted by projected NFW profiles. The NFW concentration parameter increases with decreasing satellite luminosity while being independent of the luminosity of the primary except for very bright primaries. The profiles of the faintest satellites show deviations from the NFW form with an excess at small galactocentric projected distances. In addition, we quantify how the radial distribution of satellites depends on the colour of the satellites and on the colour and concentration of their primaries.

High resolution N-body simulations of galactic cold dark matter haloes indicate that we should expect to find a few satellite galaxies around the Milky Way whose haloes have a maximum circular velocity in excess of 40 kms. Yet, with the exception of the Magellanic Clouds and the Sagittarius dwarf, which likely reside in subhaloes with significantly larger velocities than this, the bright satellites of the Milky Way all appear to reside in subhaloes with maximum circular velocities below 40 kms. As recently highlighted by Boylan-Kolchin et al., this discrepancy implies that the majority of the most massive subhaloes within a cold dark matter galactic halo are much too concentrated to be consistent with the kinematic data for the bright Milky Way satellites. Here we show that no such discrepancy exists if haloes are made of warm, rather than cold dark matter because these haloes are less concentrated on account of their typically later formation epochs. Warm dark matter is one of several possible explanations for the observed kinematics of the satellites.

We investigate the formation and evolution of satellite galaxies using smoothed particle hydrodynamics (SPH) simulations of a Milky Way(MW)-like system, focussing on the best resolved examples, analogous to the classical MW satellites. Comparing with a pure dark matter simulation, we find that the condensation of baryons has had a relatively minor effect on the structure of the satellites' dark matter halos. The stellar mass that forms in each satellite agrees relatively well over three levels of resolution (a factor of ~64 in particle mass) and scales with (sub)halo mass in a similar way in an independent semi-analytical model. Our model provides a relatively good match to the average luminosity function of the MW and M31. To establish whether the potential wells of our satellites are realistic, we measure their masses within observationally determined half-light radii, finding that the most massive examples have somewhat higher mass-to-light ratios than those derived for the MW dSphs from stellar kinematic data. A statistical test yields a ~9 percent probability that the simulated and observationally derived distributions of masses are consistent. Our results may suggest that either the MW halo is less massive than assumed in our simulations (~1.4e12 M_sun) or that there is substantial scatter in the satellite luminosity function or distribution of mass-to-light ratios at fixed host halo mass. Alternatively, feedback processes not properly captured by our simulations may have reduced the central densities of (sub)halos, or the subhalos may have initially formed with lower concentrations as would be the case, for example, if the dark matter were made of warm, rather than cold particles.

We study the luminosity function of satellite galaxies around isolated primaries using the Sloan Digital Sky Survey (SDSS) spectroscopic and photometric galaxy samples. We select isolated primaries from the spectroscopic sample and search for potential satellites in the much deeper photometric sample. For primaries of similar luminosity to the Milky Way and M31, we are able to stack as many as $\sim\negthinspace 20,000$ galaxy systems to obtain robust statistical results. We derive the satellite luminosity function extending almost 8 magnitudes fainter than the primary galaxy. We also determine how the satellite luminosity function varies with the luminosity, colour and concentration of the primary. We find that, in the mean, isolated primaries of comparable luminosity to the Milky Way and M31 contain about a factor of two fewer satellites brighter than $M_V=-14$ than the average of the Milky Way and M31.

We apply a basis function expansion method to create a time-evolving density/potential approximation of the late growth of simulated N-body dark matter haloes. We demonstrate how the potential of a halo from the Aquarius Project can be accurately represented by a small number of basis functions, and show that the halo expansion (HEX) method provides a way to replay simulations. We explore the level of accuracy of the technique as well as some of its limitations. We find that the number of terms included in the expansion must be large enough to resolve the large-scale distribution and shape of the halo but, beyond this, additional terms result in little further improvement. Particle and subhalo orbits can be integrated in this realistic, time-varying halo potential approximation, at much lower cost than the original simulation, with high fidelity for many individual orbits, and a good match to the distributions of orbital energy and angular momentum. Statistically, the evolution of structural subhalo properties, such as mass, half-mass radius and characteristic circular velocity, are very well reproduced in the halo expansion approximation over several gigayears. We demonstrate an application of the technique by following the evolution of an orbiting subhalo at much higher resolution than can be achieved in the original simulation. Our method represents a significant improvement over commonly used techniques based on static analytical descriptions of the halo potential.

We describe and apply a simple prescription for defining connected structures in galaxy redshift surveys. The method is based upon two passes with a friends-of-friends groupfinder. The first pass uses a cylindrical linking volume to find galaxy groups and clusters, in order to suppress the line-of-sight smearing introduced by the large random velocities of galaxies within these deep potential wells. The second pass, performed with a spherical linking volume, identifies the connected components. This algorithm has been applied to the 2dFGRS, within which it picks out a total of 7,603 systems containing at least two galaxies and having a mean redshift less than 0.12. Connected systems with many members appear filamentary in nature, and the algorithm recovers two particularly large filaments within the 2dFGRS. For comparison, the algorithm is has also been applied to LambdaCDM mock galaxy surveys. While the model population of such systems is broadly similar to that in the 2dFGRS, it does not generally contain such extremely large structures.

[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, Comment: 23 pages, 14 figures, 2 tables, submitted to MNRAS. Comments welcome

We investigate the baryon fraction in dark matter haloes formed in non-radiative gas-dynamical simulations of the LambdaCDM cosmogony. By combining a realisation of the Millennium Simulation (Springel et al.) with a simulation of a smaller volume focussing on dwarf haloes, our study spans five decades in halo mass, from 10^10 Msun/h to 10^15 Msun/h. We find that the baryon fraction within the halo virial radius is typically 90% of the cosmic mean, with an rms scatter of 6%, independently of redshift and of halo mass down to the smallest resolved haloes. Our results show that, contrary to the proposal of Mo et al. (2005), pre-virialisation gravitational heating is unable to prevent the collapse of gas within galactic and proto-galactic haloes, and confirm the need for non-gravitational feedback in order to reduce the efficiency of gas cooling and star formation in dwarf galaxy haloes. Simulations including a simple photoheating model (where a gas temperature floor of T_{floor} = 2x10^4 K is imposed from z=11) confirm earlier suggestions that photoheating can only prevent the collapse of baryons in systems with virial temperatures T_{200} < ~2.2 T_{floor} ~ 4.4x10^4 K (corresponding to a virial mass of M_{200} ~ 10^10 Msun/h and a circular velocity of V_{200} ~ 35 km/s). Photoheating may thus help regulate the formation of dwarf spheroidals and other galaxies at the extreme faint-end of the luminosity function, but it cannot, on its own, reconcile the abundance of sub-L* galaxies with the vast number of dwarf haloes expected in the LambdaCDM cosmogony. The lack of evolution or mass dependence seen in the baryon fraction augurs well for X-ray cluster studies that assume a universal and non-evolving baryon fraction to place constraints on cosmological parameters.

We investigate the spins and shapes of over a million dark matter haloes identified at z=0 in the Millennium simulation. Our sample spans halo masses ranging from dwarf galaxies to rich galaxy clusters. The very large dynamic range of this simulation enables the distribution of spins and shapes and their variation with halo mass and environment to be characterised with unprecedented precision. We compare results for haloes identified using different algorithms, including a novel one based on the branches of the halo merger trees. We investigate (and remove) biases in the estimate of angular momentum introduced by both the algorithms themselves and by numerical effects. We find that for this many objects, the traditional lognormal function is no longer an adequate description of the distribution of the dimensionless spin parameter lambda, and we provide a different function that gives a better fit. The variation of spin with halo mass is weak but detectable, although the trend depends strongly on the halo definition used. The haloes exhibit a range of shapes, with a preference for prolateness over oblateness. More massive haloes tend to be less spherical and more prolate. We find that the more spherical haloes have spin in the median, and those closest to spherical have a spin independent of mass. The most massive have a spin independent of shape. The majority of haloes have their angular momentum vector aligned with their minor axis and perpendicular to their major axis. We find that higher tend to be more clustered, with a stronger effect for more massive haloes. (abridged)

We use the ‘Evolution and Assembly of GaLaxies and their Environments’ (EAGLE) suite of hydrodynamical cosmological simulations to measure offsets between the centres of stellar and dark matter components of galaxies. We find that the vast majority (>95 per cent) of the simulated galaxies display an offset smaller than the gravitational softening length of the simulations (Plummer-equivalent ϵ = 700 pc), both for field galaxies and satellites in clusters and groups. We also find no systematic trailing or leading of the dark matter along a galaxy's direction of motion. The offsets are consistent with being randomly drawn from a Maxwellian distribution with σ ≤ 196 pc. Since astrophysical effects produce no feasible analogues for the 1.62+0.47−0.49 kpc offset recently observed in Abell 3827, the observational result is in tension with the collisionless cold dark matter model assumed in our simulations.

We estimate the rate of dark matter scattering in collapsed structures throughout the history of the Universe. If the scattering cross-section is velocity-independent, then the canonical picture is correct that scatterings occur mainly at late times. The scattering rate peaks slightly at redshift z~6, and remains significant today. Half the scatterings occur after z~1, in structures more massive than 10^12 M_sun. Within a factor of two, these numbers are robust to changes in the assumed astrophysics, and the scatterings would be captured in cosmological simulations. However, for particle physics models with a velocity-dependent cross-section (as for Yukawa potential interactions via a massive mediator), the scattering rate peaks before z~20, in objects with mass less than 10^4 M_sun. These precise values are sensitive to the redshift-dependent mass-concentration relation and the small-scale cutoff in the matter power spectrum. In extreme cases, the qualitative effect of early interactions may be reminiscent of warm dark matter and strongly affect the subsequent growth of structure. However, these scatterings are being missed in existing cosmological simulations with limited mass resolution., v2 matches accepted MNRAS version

The 2dFGRS is used in conjunction with the 2MASSXSC to study the near-IR light and stellar mass content of the local Universe. Mock galaxy catalogues, constructed from cosmological N-body simulations and semi-analytical galaxy formation models, are used to gauge the accuracy with which quantities can be recovered. The mean luminosity densities of the Universe are found to be rho_J=(3.57+/-0.11)*10^8 h Lsol/Mpc^3 and rho_KS=(7.04+/-0.23)*10^8 h Lsol/Mpc^3 (statistical uncertainty). Using the 2PIGG catalogue, the group dynamical mass-to-light ratio in the K_S band is found to increase by a factor of ~3 when going from groups with total bJ-band luminosities of 3*10^10 h^-2 Lsol to rich clusters, which have typical values of Upsilon_K~80 h Upsilon_sol. Taking into account the bias introduced by uncertainties in estimating galaxy stellar masses from luminosities, a value of Omega_* h=(0.99+/-0.03)*10^-3 is measured, assuming that a Kennicutt stellar IMF is applicable to all galaxies. The 2PIGGs are then used to study the distribution of the stellar content of the local Universe. The three main conclusions are: (1) a slowly rising stellar M/L_KS is found with the clusters having the largest value of ~0.6 Upsilon_sol, (2) in contrast, the fraction of mass in stars decreases with increasing group size, reaching ~5*10^-3 h for the rich clusters, and (3) in answer to the question posed in the title, most stellar mass is contained in Local Group-sized objects (M~2*10^12 h^-1 Msol) with only ~2% in clusters with M>5*10^14 h^-1 Msol.

We present a detailed analysis of the disk component of a simulated galaxy formed in the ΛCDM cosmogony. At redshift z = 0, two distinct dynamical components are easily identified solely on the basis of the orbital parameters of stars in the galaxy: a slowly rotating, centrally concentrated spheroid and a disklike component largely supported by rotation. The disk may be further decomposed into a thin, dynamically cold component with stars on nearly circular orbits and a hotter, thicker component with orbital parameters transitional between the thin disk and the spheroid. Supporting evidence for the presence of distinct thick- and thin-disk components is provided, as in the Milky Way, by the double-exponential vertical structure of the disk and in abrupt changes in the vertical velocity distribution as a function of stellar age. The dynamical origin of these components offers intriguing clues to the assembly of spheroids and disks in the Milky Way and other spiral galaxies. The spheroid is old and has essentially no stars younger than the time elapsed since the last major accretion event, ~8 Gyr ago for the system we consider here. The majority of thin-disk stars, on the other hand, form after the merging activity is over, although a significant fraction (~15%) of thin-disk stars are old enough to predate the last major merger event. This unexpected population of old-disk stars consists mainly of the tidal debris of satellites whose orbital plane was coincident with the disk and whose orbits were circularized by dynamical friction prior to full disruption. More than half of the stars in the thick disk share this origin, part of a trend that becomes more pronounced with age: 9 out of 10 stars presently in the old (age of 10 Gyr) disk component were actually brought into the disk by satellites. By contrast, only one in two stars belonging to the old spheroid are tidal debris; the rest may be traced to a major merger event that dispersed the luminous progenitor at z ~ 1.5 and seeded the formation of the spheroid. Our results highlight the role of satellite accretion events in shaping the disk, as well as the spheroidal, component and reveal some of the clues to the assembly process of a galaxy preserved in the detailed dynamics of old stellar populations.

We present a detailed analysis of the dynamical and photometric properties of a disk galaxy simulated in the cold dark matter (CDM) cosmogony. The galaxy is assembled through a number of high-redshift mergers followed by a period of quiescent accretion after z1 that lead to the formation of two distinct dynamical components: a spheroid of mostly old stars and a rotationally supported disk of younger stars. The surface brightness profile is very well approximated by the superposition of an R1/4 spheroid and an exponential disk. Each photometric component contributes a similar fraction of the total luminosity of the system, although less than a quarter of the stars form after the last merger episode at z1. In the optical bands the surface brightness profile is remarkably similar to that of Sab galaxy UGC 615, but the simulated galaxy rotates significantly faster and has a declining rotation curve dominated by the spheroid near the center. The decline in circular velocity is at odds with observation and results from the high concentration of the dark matter and baryonic components, as well as from the relatively high mass-to-light ratio of the stars in the simulation. The simulated galaxy lies 1 mag off the I-band Tully-Fisher relation of late-type spirals but seems to be in reasonable agreement with Tully-Fisher data on S0 galaxies. In agreement with previous simulation work, the angular momentum of the luminous component is an order of magnitude lower than that of late-type spirals of similar rotation speed. This again reflects the dominance of the slowly rotating, dense spheroidal component, to which most discrepancies with observation may be traced. On its own, the disk component has properties rather similar to those of late-type spirals: its luminosity, its exponential scale length, and its colors are all comparable to those of galaxy disks of similar rotation speed. This suggests that a different form of feedback than adopted here is required to inhibit the efficient collapse and cooling of gas at high redshift that leads to the formation of the spheroid. Reconciling, without fine-tuning, the properties of disk galaxies with the early collapse and high merging rates characteristic of hierarchical scenarios such as CDM remains a challenging, yet so far elusive, proposition.

We examine in detail the dynamical structure of an elliptical galaxy simulated in the Lambda CDM scenario. The morphology of the galaxy evolves dramatically over time in response to the mode and timing of mass accretion; smooth deposition of cooled gas leads to the formation of centrifugally supported disks, whilst major mergers disperse stellar disks into spheroids. These two modes of accretion alternate successively until z~0.6, when the galaxy undergoes one last major (1:2) merger that consumes much of the remaining gas into stars. Little gas cools and accretes subsequently and, as a result, most stars at z=0 are in a spheroidal component that resembles present-day elliptical galaxies. Dynamically, the galaxy is well approximated by an E4 oblate rotator. Boxy isophotes are obtained when the galaxy is seen face-on and Vrot/sigma<, 17 pages, 20 color figures, replaced with published version ApJ

We use gravitational lens models and X-ray spectral analysis of ten X-ray luminous galaxy clusters at z~0.2 to study the impact of cluster substructure on attempts to normalize the matter power spectrum. We estimate that unrelaxed clusters are 30% hotter than relaxed clusters causing sigma8 to be over-estimated by 20% if the cluster selection function is not accounted for correctly. This helps to explain the wide range in sigma8 derived from different techniques, sigma8~0.6-1, and offers a physically motivated explanation for some of the discrepancy. We also identify two further systematics in our analysis: (i) extrapolation of small field-of-view mass measurements to the cluster virial radius and (ii) projection of 3-dimensional masses contained in numerical simulations to the 2-dimensional information that is available from observations. We combine quantitative estimates of these two effects with our model fitting to estimate from the current data that sigma8=0.75+/-0.05(statistical)+/-0.15(systematic), where the systematic error reflects the extrapolation and projection uncertainties. All three systematics (substructure, extrapolation and projection) are fundamental to future cluster-based measurements of sigma8 regardless of the techniques employed. However, we identify gravitational lensing as the tool of choice for such studies, because a combination of strong- and weak-lensing offers the most direct route to control the systematics and thus achieve an unbiased comparison between observation and theory., Submitted to Astrophysical Journal Letters, 4 pages, 2 figures

We investigate the influence of the cooling epoch on the formation of galaxies in a cold dark matter dominated universe. Isolated haloes, with circular speeds typical of spiral galaxies, have been selected from a low resolution numerical simulation for re-simulation at higher resolution with dark matter and gas components. The initial conditions are evolved with two smoothed particle hydrodynamics codes, TREESPH and GRAPESPH. In previous SPH simulations, strong outward transport of angular momentum has led to the formation of disc-like systems with much smaller angular momenta than observed in real disc galaxies. Here we investigate whether this problem can be circumvented if feedback processes prevent disc formation until late epochs. The results of varying the cooling epoch for each of five different haloes are analysed. When cooling and star formation occur at early times, stellar discs are destroyed during merger events and we observe similar catastrophic transport of angular momentum as seen in previous work. With cooling suppressed until z=1, discs can form by the present day with angular momenta comparable to those of observed disc galaxies. We conclude that feedback processes, which prevent gas from collapsing until late epochs, are an essential ingredient in disc galaxy formation.

The evolution of the abundance of galaxy clusters depends sensitively on the value of the cosmological density parameter, Ω0. Recent ASCA data are used to quantify this evolution as measured by the cluster X-ray temperature function. A χ2 minimization fit to the cumulative temperature function, as well as a maximum-likelihood estimate (which requires additional assumptions about cluster luminosities), leads to the estimate Ω0 ≈ 0.45 ± 0.25 (1σ statistical error). Various systematic uncertainties are considered, none of which significantly enhances the probability that Ω0 = 1. These conclusions hold for models with or without a cosmological constant, i.e., with Λ0 = 0 or Λ0 = 1 − Ω0. The statistical uncertainties are at least as large as any of the individual systematic errors that have been considered here, suggesting that additional temperature measurements of distant clusters will allow an improvement in this estimate. An alternative method that uses the highest redshift clusters to place an upper limit on Ω0 is also presented and tentatively applied, with the result that Ω0 1 can be ruled out at the 98 per cent confidence level. Whilst this method does not require a well-defined statistical sample of distant clusters, there are still modelling uncertainties that preclude a firmer conclusion at this time.

We have re-examined the nature of the cluster galaxy luminosity function using the data from the Edinburgh-Durham Southern Galaxy Catalogue and the Edinburgh-Milano Redshift Survey. We derive a best fit luminosity function over the range M(bj)=-18 to -21, for a composite sample of 22 of the richer clusters that has M*=-20.16+/-0.02 and alpha=-1.22+/-0.04. The dominant error in these values results from the choice of background subtraction method. From extensive simulations we can show that when the LF is fitted over this narrow range, it is difficult to discriminate against bright values of M*in the single cluster fits, but that faint values provide a strong test of the universality of the luminosity function. We find that all the individual cluster data are well fit by a Schechter function with alpha fixed at -1.25, and that, To be published in MNRAS 15 pages, 9 postscript figures

The population of rich galaxy clusters evolves much more rapidly in a universe with critical density than one with low density, thus offering the possibility of determining the cosmological density parameter, Omega_0. We quantify this evolution using the Press-Schechter formalism which we extend to flat models with a cosmological constant. Using new large N-body simulations, we verify that this formalism accurately predicts the abundance of rich clusters as a function of redshift in various cosmologies. We normalise the models by comparing them to the local abundance of clusters as a function of their X-ray temperature which we rederive from data compiled by Henry & Arnaud. This gives values of the rms density fluctuation in spheres of radius 8 Mpc/h of sigma_8 = (0.50+/- 0.04) Omega_0^{-0.47+0.10 Omega_0} if Lambda_0=0 and sigma_8 = (0.50 +/- 0.04) Omega_0^{-0.53+0.13 Omega_0} if Lambda_0=1-Omega_0. These values depend very weakly on the shape of the power spectrum. We then examine how the distributions of mass, X-ray temperature and Sunyaev-Zel'dovich decrement evolve as a function of Omega_0. We present the expected distributions at z=0.33 and z=0.5 and the predicted number counts of the largest clusters. We find that even at z=0.33, these distributions depend very strongly on Omega_0 and only weakly on Lambda_0. For example, at this redshift, we expect 20 times as many clusters per comoving volume with M>3.5 10^{14} Msol/h and 5 times as many clusters with kT>5 keV if Omega_0=0.3 than if Omega_0=1. The splitting in the integrated counts is enhanced by the larger volume element in low Omega_0 models. There is therefore a real prospect of estimating Omega_0 from forthcoming surveys of intermediate redshift clusters that will determine their masses, X-ray temperatures or SZ decrements.

Data from the Lunar Exploration Neutron Detector (LEND) Collimated Sensors for Epithermal Neutrons (CSETN) are used in conjunction with a model based on results from the Lunar Prospector mission to quantify the extent of the background in the LEND CSETN. A simple likelihood analysis implies that at least 90% of the lunar component of the LEND CSETN flux results from high energy epithermal neutrons passing through the walls of the collimator. Thus, the effective full-width at half-maximum of the LEND CSETN is comparable with that of the omni-directional Lunar Prospector Neutron Spectrometer. The resulting map of high energy epithermal neutrons offers the opportunity to probe the hydrogen abundance at low latitudes, and provide constraints on the distribution of lunar water., 19 pages, 14 figures, updated version accepted by ApJ

[1] Although controversial in its physical form, there is mounting evidence of hydrogen enhancements at the lunar poles. The permanently shadowed locales are potential sites for significant concentrations of cold-trapped volatiles, including water ice. We derive maps of the lunar hydrogen distribution near the poles by applying a pixon image reconstruction algorithm to the Lunar Prospector epithermal neutron data coupled with a new map of cold trap locations derived from the KAGUYA (SELENE) altimetry measurements. The results presented in this article require the hydrogen to be concentrated into “cold traps.”

We investigate the magnitude and internal alignment of the angular momentum of cold dark matter haloes in simulations with and without baryons. We analyse the cumulative angular momentum profiles of hundreds of thousands of haloes in the Millennium simulation and in a smaller, but higher resolution, simulation, in total spanning 5 orders of magnitude in mass. For haloes of a given mass, the median specific angular momentum increases as j(, Comment: 21 pages, 27 figures. Accepted for publication in MNRAS (abstract shown here is slightly abridged). Replaced with accepted version, which adds significant clarifications and expansion of discussions, including additional figures. All science content remains unchanged