Showing 1–50 of 76 results for author: Eke, V

A recent impact origin of Saturn's rings and mid-sized moons

Authors: Luís F. A. Teodoro, Jacob A. Kegerreis, Paul R. Estrada, Matija Ćuk, Vincent R. Eke, Jeffrey N. Cuzzi, Richard J. Massey, Thomas D. Sandnes

Abstract: We simulate the collision of precursor icy moons analogous to Dione and Rhea as a possible origin for Saturn's remarkably young rings. Such an event could have been triggered a few hundred million years ago by resonant instabilities in a previous satellite system. Using high-resolution smoothed particle hydrodynamics simulations, we find that this kind of impact can produce a wide distribution of… ▽ More We simulate the collision of precursor icy moons analogous to Dione and Rhea as a possible origin for Saturn's remarkably young rings. Such an event could have been triggered a few hundred million years ago by resonant instabilities in a previous satellite system. Using high-resolution smoothed particle hydrodynamics simulations, we find that this kind of impact can produce a wide distribution of massive objects and scatter material throughout the system. This includes the direct placement of pure-ice ejecta onto orbits that enter Saturn's Roche limit, which could form or rejuvenate rings. In addition, fragments and debris of rock and ice totalling more than the mass of Enceladus can be placed onto highly eccentric orbits that would intersect with any precursor moons orbiting in the vicinity of Mimas, Enceladus, or Tethys. This could prompt further disruption and facilitate a collisional cascade to distribute more debris for potential ring formation, the re-formation of the present-day moons, and evolution into an eventual cratering population of planeto-centric impactors. △ Less

Submitted 28 September, 2023; v1 submitted 26 September, 2023; originally announced September 2023.

Comments: 22 pages, 13 figures, published in ApJ. Animations available at https://www.youtube.com/@jkeger_et_al

Journal ref: ApJ, 955, 2, 137 (2023)

Immediate origin of the Moon as a post-impact satellite

Authors: Jacob A. Kegerreis, Sergio Ruiz-Bonilla, Vincent R. Eke, Richard J. Massey, Thomas D. Sandnes, Luís F. A. Teodoro

Abstract: The Moon is traditionally thought to have coalesced from the debris ejected by a giant impact onto the early Earth. However, such models struggle to explain the similar isotopic compositions of Earth and lunar rocks at the same time as the system's angular momentum, and the details of potential impact scenarios are hotly debated. Above a high resolution threshold for simulations, we find that gian… ▽ More The Moon is traditionally thought to have coalesced from the debris ejected by a giant impact onto the early Earth. However, such models struggle to explain the similar isotopic compositions of Earth and lunar rocks at the same time as the system's angular momentum, and the details of potential impact scenarios are hotly debated. Above a high resolution threshold for simulations, we find that giant impacts can immediately place a satellite with similar mass and iron content to the Moon into orbit far outside the Earth's Roche limit. Even satellites that initially pass within the Roche limit can reliably and predictably survive, by being partially stripped then torqued onto wider, stable orbits. Furthermore, the outer layers of these directly formed satellites are molten over cooler interiors and are composed of around 60% proto-Earth material. This could alleviate the tension between the Moon's Earth-like isotopic composition and the different signature expected for the impactor. Immediate formation opens up new options for the Moon's early orbit and evolution, including the possibility of a highly tilted orbit to explain the lunar inclination, and offers a simpler, single-stage scenario for the origin of the Moon. △ Less

Submitted 4 October, 2022; originally announced October 2022.

Comments: 15 pages, 12 figures, published in ApJL. Animations available at https://www.youtube.com/channel/UCGnGqCiYLWvYN992mBkBnhg/

Journal ref: ApJL 937 L40 (2022)

Dealing with density discontinuities in planetary SPH simulations

Authors: Sergio Ruiz-Bonilla, Josh Borrow, Vincent R. Eke, Jacob A. Kegerreis, Richard J. Massey, Thomas D. Sandnes, Luis F. A. Teodoro

Abstract: Density discontinuities cannot be precisely modelled in standard formulations of smoothed particles hydrodynamics (SPH) because the density field is defined smoothly as a kernel-weighted sum of neighbouring particle masses. This is a problem when performing simulations of giant impacts between proto-planets, for example, because planets typically do have density discontinuities both at their surfa… ▽ More Density discontinuities cannot be precisely modelled in standard formulations of smoothed particles hydrodynamics (SPH) because the density field is defined smoothly as a kernel-weighted sum of neighbouring particle masses. This is a problem when performing simulations of giant impacts between proto-planets, for example, because planets typically do have density discontinuities both at their surfaces and at any internal boundaries between different materials. The inappropriate densities in these regions create artificial forces that effectively suppress mixing between particles of different material and, as a consequence, this problem introduces a key unknown systematic error into studies that rely on SPH simulations. In this work we present a novel, computationally cheap method that deals simultaneously with both of these types of density discontinuity in SPH simulations. We perform standard hydrodynamical tests and several example giant impact simulations, and compare the results with standard SPH. In a simulated Moon-forming impact using $10^7$ particles, the improved treatment at boundaries affects at least 30% of the particles at some point during the simulation. △ Less

Submitted 7 February, 2022; v1 submitted 1 February, 2022; originally announced February 2022.

Comments: 9 pages, 8 figures. Submitted to MNRAS

The surprising accuracy of isothermal Jeans modelling of self-interacting dark matter density profiles

Authors: Andrew Robertson, Richard Massey, Vincent Eke, Joop Schaye, Tom Theuns

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

Submitted 29 January, 2021; v1 submitted 16 September, 2020; originally announced September 2020.

Comments: 27 pages, 19 figures, accepted in MNRAS, updated to match accepted version

Atmospheric Erosion by Giant Impacts onto Terrestrial Planets: A Scaling Law for any Speed, Angle, Mass, and Density

Authors: Jacob A. Kegerreis, Vincent R. Eke, David C. Catling, Richard J. Massey, Luis F. A. Teodoro, Kevin J. Zahnle

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

Submitted 30 September, 2020; v1 submitted 8 July, 2020; originally announced July 2020.

Comments: Published in ApJL. 12 pages, 6 figures

The effect of pre-impact spin on the Moon-forming collision

Authors: Sergio Ruiz-Bonilla, Vincent R. Eke, Jacob A. Kegerreis, Richard J. Massey, Luis F. A. Teodoro

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

Submitted 6 July, 2020; originally announced July 2020.

Comments: Submitted to MNRAS. 11 pages, 9 figures

Space-Based Measurement of the Neutron Lifetime using Data from the Neutron Spectrometer on NASA's MESSENGER Mission

Authors: Jack T. Wilson, David J. Lawrence, Patrick N. Peplowski, Vincent R. Eke, Jacob A. Kegerreis

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

Submitted 17 June, 2020; originally announced June 2020.

Journal ref: Phys. Rev. Research 2, 023316 (2020)

Understanding the large inferred Einstein radii of observed low-mass galaxy clusters

Authors: Andrew Robertson, Richard Massey, Vincent Eke

Abstract: We assess a claim that observed galaxy clusters with mass ~10^14 Msun are more centrally concentrated than predicted in LCDM. We generate mock strong gravitational lensing observations, taking the lenses from a cosmological hydrodynamical simulation, and analyse them in the same way as the real Universe. The observed and simulated lensing arcs are consistent with one another, with three main effec… ▽ More We assess a claim that observed galaxy clusters with mass ~10^14 Msun are more centrally concentrated than predicted in LCDM. We generate mock strong gravitational lensing observations, taking the lenses from a cosmological hydrodynamical simulation, and analyse them in the same way as the real Universe. The observed and simulated lensing arcs are consistent with one another, with three main effects responsible for the previously claimed inconsistency. First, galaxy clusters containing baryonic matter have higher central densities than their counterparts simulated with only dark matter. Second, a sample of clusters selected because of the presence of pronounced gravitational lensing arcs preferentially finds centrally concentrated clusters with large Einstein radii. Third, lensed arcs are usually straighter than critical curves, and the chosen image analysis method (fitting circles through the arcs) overestimates the Einstein radii. After accounting for these three effects, LCDM predicts that galaxy clusters should produce giant lensing arcs that match those in the observed Universe. △ Less

Submitted 1 June, 2020; v1 submitted 21 February, 2020; originally announced February 2020.

Comments: 8 pages, 3 figures, updated to match MNRAS version

Atmospheric Erosion by Giant Impacts onto Terrestrial Planets

Authors: J. A. Kegerreis, V. R. Eke, R. J. Massey, L. F. A. Teodoro

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

Submitted 20 July, 2020; v1 submitted 7 February, 2020; originally announced February 2020.

Comments: 19 pages, 13 figures, published in ApJ

What does strong gravitational lensing? The mass and redshift distribution of high-magnification lenses

Authors: Andrew Robertson, Graham P. Smith, Richard Massey, Vincent Eke, Mathilde Jauzac, Matteo Bianconi, Dan Ryczanowski

Abstract: 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 redshif… ▽ More 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 <10^12 solar masses for lensing. We also calculate the probability that point sources at different redshifts are strongly lensed. At low redshift, high magnifications are extremely unlikely. Each z=0.5 source produces, on average, 5x10^-7 images with magnification greater than ten; for z =2 this increases to about 2x10^-5. Our results imply that searches for strongly lensed optical transients, including the optical counterparts to strongly lensed gravitational waves, can be optimized by monitoring massive galaxies, groups and clusters rather than concentrating on an individual population of lenses. △ Less

Submitted 1 June, 2020; v1 submitted 4 February, 2020; originally announced February 2020.

Comments: 15 pages, 5 figures, updated to match MNRAS version

Planetary Giant Impacts: Convergence of High-Resolution Simulations using Efficient Spherical Initial Conditions and SWIFT

Authors: J. A. Kegerreis, V. R. Eke, P. G. Gonnet, D. G. Korycansky, R. J. Massey, M. Schaller, L. F. A. Teodoro

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

Submitted 7 April, 2020; v1 submitted 28 January, 2019; originally announced January 2019.

Comments: Fixed typo in Appx B equation and updated urls

Observable tests of self-interacting dark matter in galaxy clusters: cosmological simulations with SIDM and baryons

Authors: Andrew Robertson, David Harvey, Richard Massey, Vincent Eke, Ian G. McCarthy, Mathilde Jauzac, Baojiu Li, Joop Schaye

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

Submitted 1 June, 2020; v1 submitted 12 October, 2018; originally announced October 2018.

Comments: 19 pages, 13 figures, updated to match MNRAS version

Consequences of Giant Impacts on Early Uranus for Rotation, Internal Structure, Debris, and Atmospheric Erosion

Authors: J. A. Kegerreis, L. F. A. Teodoro, V. R. Eke, R. J. Massey, D. C. Catling, C. L. Fryer, D. G. Korycansky, M. S. Warren, K. J. Zahnle

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

Submitted 3 July, 2018; v1 submitted 19 March, 2018; originally announced March 2018.

Comments: Published in ApJ. 14 pages, 11 figures

Report number: LA-UR-18-22206

Journal ref: ApJ 861 (2018) 52

Image Reconstruction Techniques in Neutron and Gamma-Ray Spectroscopy: Improving Lunar Prospector Data

Authors: Jack T. Wilson, David J. Lawrence, Patrick N. Peplowski, Joshua T. S. Cahill, Vincent R. Eke, Richard J. Massey, Luis F. A. Teodoro

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

Submitted 26 February, 2018; originally announced February 2018.

Comments: 17 pages, 12 figures, submitted to JGR-Planets

The diverse density profiles of galaxy clusters with self-interacting dark matter plus baryons

Authors: Andrew Robertson, Richard Massey, Vincent Eke, Sean Tulin, Hai-Bo Yu, Yannick Bahé, David J. Barnes, Richard G. Bower, Robert A. Crain, Claudio Dalla Vecchia, Scott T. Kay, Matthieu Schaller, Joop Schaye

Abstract: We present the first simulated galaxy clusters (M_200 > 10^14 Msun) with both self-interacting dark matter (SIDM) and baryonic physics. They exhibit a greater diversity in both dark matter and stellar density profiles than their counterparts in simulations with collisionless dark matter (CDM), which is generated by the complex interplay between dark matter self-interactions and baryonic physics. D… ▽ More We present the first simulated galaxy clusters (M_200 > 10^14 Msun) with both self-interacting dark matter (SIDM) and baryonic physics. They exhibit a greater diversity in both dark matter and stellar density profiles than their counterparts in simulations with collisionless dark matter (CDM), which is generated by the complex interplay between dark matter self-interactions and baryonic physics. Despite variations in formation history, we demonstrate that analytical Jeans modelling predicts the SIDM density profiles remarkably well, and the diverse properties of the haloes can be understood in terms of their different final baryon distributions. △ Less

Submitted 22 February, 2018; v1 submitted 24 November, 2017; originally announced November 2017.

Comments: 6 pages, 1 figure, v2 matches accepted MNRAS: Letters version

Equatorial locations of water on Mars: Improved resolution maps based on Mars Odyssey Neutron Spectrometer data

Authors: Jack T. Wilson, Vincent R. Eke, Richard J. Massey, Richard C. Elphic, William C. Feldman, Sylvestre Maurice, Luis F. A. Teodoro

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

Submitted 1 August, 2017; originally announced August 2017.

Comments: Accepted for publication in Icarus, 16 pages, 12 figures

Cosmic particle colliders: simulations of self-interacting dark matter with anisotropic scattering

Authors: Andrew Robertson, Richard Massey, Vincent Eke

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

Submitted 18 March, 2017; v1 submitted 12 December, 2016; originally announced December 2016.

Comments: 12 pages, 8 figures, v2 matches accepted MNRAS version

Evidence for a Localised Source of the Argon in the Lunar Exosphere

Authors: Jacob A. Kegerreis, Vincent R. Eke, Richard J. Massey, Simon K. Beaumont, Rick C. Elphic, Luis F. Teodoro

Abstract: We perform the first tests of various proposed explanations for observed features of the Moon's argon exosphere, including models of: 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 Ex… ▽ More We perform the first tests of various proposed explanations for observed features of the Moon's argon exosphere, including models of: 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. A localised source can explain the observations, with a trade-off between an unexpectedly localised 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\times10^{21}$ atoms/s. A more diffuse source, weighted by the near-surface potassium, requires much higher rates of ~$1.1\times10^{22}$ atoms/s, corresponding to a mean lifetime of just 1.4 lunar days. We do not address the mechanism for producing a localised 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. △ Less

Submitted 31 October, 2017; v1 submitted 7 December, 2016; originally announced December 2016.

Comments: Published in JGR: Planets. 19 pages, 10 figures

Journal ref: Journal of Geophysical Research: Planets, Early View 26 October 2017

How thick are Mercury's polar water ice deposits?

Authors: Vincent R. Eke, David J. Lawrence, Luıs F. A. Teodoro

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

Submitted 16 November, 2016; originally announced November 2016.

Comments: 10 pages, 6 figures, submitted to Icarus

What does the Bullet Cluster tell us about self-interacting dark matter?

Authors: Andrew Robertson, Richard Massey, Vincent Eke

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

Submitted 2 February, 2017; v1 submitted 13 May, 2016; originally announced May 2016.

Comments: 21 pages, 16 figures, v2 matches accepted MNRAS version

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

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

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

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

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

The offsets between galaxies and their dark matter in Lambda-CDM

Authors: Matthieu Schaller, Andrew Robertson, Richard Massey, Richard G. Bower, Vincent R. Eke

Abstract: 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%) of the simulated galaxies display an offset smaller than the gravitational softening length of the simulations (Plummer-equivalent $ε= 700$ pc… ▽ More 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%) 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 $σ\leq 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. △ Less

Submitted 22 July, 2015; v1 submitted 20 May, 2015; originally announced May 2015.

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

Journal ref: MNRASL, 2015, Volume 453, Issue 1, pp. L58-L62

The effect of craters on the lunar neutron flux

Authors: V. R. Eke, K. E. Bower, S. Diserens, M. Ryder, P. E. L. Yeomans, L. F. A. Teodoro, R. C. Elphic, W. C. Feldman, B. Hermalyn, C. M. Lavelle, D. J. Lawrence

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

Submitted 17 August, 2015; v1 submitted 20 May, 2015; originally announced May 2015.

Comments: 14 pages, 13 figures, minor changes to match published version

Journal ref: 2015, J. Geophys. Res. Planets, 120

Self-Interacting Dark Matter Scattering Rates Through Cosmic Time

Authors: Andrew Robertson, Richard Massey, Vincent Eke, Richard Bower

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

Submitted 16 January, 2017; v1 submitted 8 May, 2015; originally announced May 2015.

Comments: v2 matches accepted MNRAS version

The Local-time variations of Lunar Prospector epithermal-neutron data

Authors: L. F. A. Teodoro, D. J. Lawrence, V. R. Eke, R. E. Elphic, W. C. Feldman, S. Maurice, M. A. Siegler, D. A. Paige

Abstract: We assess local-time variations of epithermal-neutron count rates measured by the Lunar Prospector Neutron Spectrometer. We investigate the nature of these variations and find no evidence to support the idea that such variations are caused by diurnal variations of hydrogen concentration across the lunar surface. Rather we find an anticorrelation between instrumental temperature and epithermal-neut… ▽ More We assess local-time variations of epithermal-neutron count rates measured by the Lunar Prospector Neutron Spectrometer. We investigate the nature of these variations and find no evidence to support the idea that such variations are caused by diurnal variations of hydrogen concentration across the lunar surface. Rather we find an anticorrelation between instrumental temperature and epithermal-neutron count rate. We have also found that the measured counts are dependent on the temperatures of the top decimeters of the lunar subsurface as constrained by the Lunar Reconnaissance Orbiter Diviner Lunar Radiometer Experiment temperature measurements. Finally, we have made the first measurement of the effective leakage depth for epithermal-neutrons of ~20 cm. △ Less

Submitted 22 January, 2015; originally announced January 2015.

Comments: 8 pages, 7 figures and 3 tables

Evidence for explosive silicic volcanism on the Moon from the extended distribution of thorium near the Compton-Belkovich Volcanic Complex

Authors: J. T. Wilson, V. R. Eke, R. J. Massey, R. C. Elphic, B. L. Jolliff, D. J. Lawrence, E. W. Llewellin, J. N. McElwaine, L. F. A. Teodoro

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

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

Comments: v2: 13 pages, 8 figures, accepted version, to be published in JGR-Planets; v1: 11 pages, 6 figures, submitted to JGR-Planets

Galaxy and Mass Assembly (GAMA): The halo mass of galaxy groups from maximum-likelihood weak lensing

Authors: Jiaxin Han, Vincent R. Eke, Carlos S. Frenk, Rachel Mandelbaum, Peder Norberg, Michael D. Schneider, John A. Peacock, Yipeng Jing, Ivan Baldry, Joss Bland-Hawthorn, Sarah Brough, Michael J. I. Brown, Jochen Liske, Jon Loveday, Aaron S. G. Robotham

Abstract: We present a maximum-likelihood weak lensing analysis of the mass distribution in optically selected spectroscopic Galaxy Groups (G$^3$Cv5) in the Galaxy And Mass Assembly (GAMA) survey, using background Sloan Digital Sky Survey (SDSS) photometric galaxies. The scaling of halo mass, $M_h$, with various group observables is investigated. Our main results are: 1) the measured relations of halo mass… ▽ More We present a maximum-likelihood weak lensing analysis of the mass distribution in optically selected spectroscopic Galaxy Groups (G$^3$Cv5) in the Galaxy And Mass Assembly (GAMA) survey, using background Sloan Digital Sky Survey (SDSS) photometric galaxies. The scaling of halo mass, $M_h$, with various group observables is investigated. Our main results are: 1) the measured relations of halo mass with group luminosity, virial volume and central galaxy stellar mass, $M_\star$, agree very well with predictions from mock group catalogues constructed from a GALFORM semi-analytical galaxy formation model implemented in the Millennium $Λ$CDM N-body simulation; 2) the measured relations of halo mass with velocity dispersion and projected half-abundance radius show weak tension with mock predictions, hinting at problems in the mock galaxy dynamics and their small scale distribution; 3) the median $M_h|M_\star$ measured from weak lensing depends more sensitively on the lognormal dispersion in $M_\star$ at fixed $M_h$ than it does on the median $M_\star|M_h$. Our measurements suggest an intrinsic dispersion of $σ_{\log(M_\star)}\sim 0.15$; 4) Comparing our mass estimates with those in the catalogue, we find that the G$^3$Cv5 mass can give biased results when used to select subsets of the group sample. Of the various new halo mass estimators that we calibrate using our weak lensing measurements, group luminosity is the best single-proxy estimator of group mass. △ Less

Submitted 4 November, 2014; v1 submitted 27 April, 2014; originally announced April 2014.

Comments: MNRAS accepted

Lunar polar craters -- icy, rough or just sloping?

Authors: Vincent R. Eke, Sarah A. Bartram, David A. Lane, David Smith, Luis F. A. Teodoro

Abstract: Circular Polarisation Ratio (CPR) mosaics from Mini-SAR on Chandrayaan-1 and Mini-RF on LRO are used to study craters near to the lunar north pole. The look direction of the detectors strongly affects the appearance of the crater CPR maps. Rectifying the mosaics to account for parallax also significantly changes the CPR maps of the crater interiors. It is shown that the CPRs of crater interiors in… ▽ More Circular Polarisation Ratio (CPR) mosaics from Mini-SAR on Chandrayaan-1 and Mini-RF on LRO are used to study craters near to the lunar north pole. The look direction of the detectors strongly affects the appearance of the crater CPR maps. Rectifying the mosaics to account for parallax also significantly changes the CPR maps of the crater interiors. It is shown that the CPRs of crater interiors in unrectified maps are biased to larger values than crater exteriors, because of a combination of the effects of parallax and incidence angle. Using the LOLA Digital Elevation Map (DEM), the variation of CPR with angle of incidence has been studied. For fresh craters, CPR~0.7 with only a weak dependence on angle of incidence or position interior or just exterior to the crater, consistent with dihedral scattering from blocky surface roughness. For anomalous craters, the CPR interior to the crater increases with both incidence angle and distance from the crater centre. Central crater CPRs are similar to those in the crater exteriors. CPR does not appear to correlate with temperature within craters. Furthermore, the anomalous polar craters have diameter-to-depth ratios that are lower than those of typical polar craters. These results strongly suggest that the high CPR values in anomalous polar craters are not providing evidence of significant volumes of water ice. Rather, anomalous craters are of intermediate age, and maintain sufficiently steep walls that deep regolith does not cover all rough surfaces. △ Less

Submitted 25 August, 2014; v1 submitted 17 December, 2013; originally announced December 2013.

Comments: 16 pages, 17 figures, submitted to Icarus (2014, 241, 66)

The properties of warm dark matter haloes

Authors: Mark R. Lovell, Carlos S. Frenk, Vincent R. Eke, Adrian Jenkins, Liang Gao, Tom Theuns

Abstract: 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 cosmolo… ▽ More 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 (<10^7Msun in our case) outnumber genuine haloes. We have developed a robust algorithm to identify these spurious objects and remove them from our halo catalogues. We find that the WDM subhalo mass function is suppressed by well over an order magnitude relative to the CDM case for masses <10^9Msun. Requiring that there should be at least as many subhaloes as there are observed satellites in the Milky Way leads to a conservative lower limit to the (thermal equivalent) WDM particle mass of ~1.5\rmn{keV}. WDM haloes and subhaloes have cuspy density distributions that are well described by NFW or Einasto profiles. Their central densities are lower for lower WDM particle masses and none of the models we have considered suffer from the "too big to fail" problem recently highlighted by Boylan-Kolchin et al. △ Less

Submitted 19 May, 2014; v1 submitted 6 August, 2013; originally announced August 2013.

Comments: 20 pages, 21 Figures, Published in MNRAS, minor updates to particle mass calculations: 1.4keV -> 1.5keV

How well do we know the polar hydrogen distribution on the Moon?

Authors: L. F. A. Teodoro, V. R. Eke, R. C. Elphic, W. C. Feldman, D. J Lawrence

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

Submitted 30 April, 2013; originally announced April 2013.

Comments: 12 pages, 13 figures

Spatial and luminosity distributions of galactic satellites

Authors: Quan Guo, Shaun Cole, Vincent Eke, Carlos Frenk, John Helly

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

Submitted 14 January, 2013; originally announced January 2013.

Comments: 11 pages, 10 figures, submitted to MNRAS

Constraining Extended Gamma-ray Emission from Galaxy Clusters

Authors: Jiaxin Han, Carlos S. Frenk, Vincent R. Eke, Liang Gao, Simon D. M. White, Alexey Boyarsky, Denys Malyshev, Oleg Ruchayskiy

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

Submitted 6 December, 2012; v1 submitted 29 July, 2012; originally announced July 2012.

Comments: updated to match published version, with more references added. arXiv admin note: substantial text overlap with arXiv:1201.1003

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 427, Issue 2, pp. 1651-1665 (2012)

Satellite Galaxy Number Density Profiles in the Sloan Digital Sky Survey

Authors: Quan Guo, Shaun Cole, Vincent Eke, Carlos Frenk

Abstract: 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… ▽ More 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 magnitudes fainter than their primaries. We find 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. △ Less

Submitted 5 September, 2012; v1 submitted 5 January, 2012; originally announced January 2012.

Comments: 15 pages, 13 figures, accepted by MNRAS

Evidence for extended gamma-ray emission from galaxy clusters

Authors: Jiaxin Han, Carlos S. Frenk, Vincent R. Eke, Liang Gao, Simon D. M. White

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

Submitted 1 August, 2012; v1 submitted 4 January, 2012; originally announced January 2012.

Comments: This paper is now superseded by arxiv:1207.6749

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 427, Issue 2, pp. 1651-1665 (2012)

A quantitative comparison of lunar orbital neutron data

Authors: V. R. Eke, L. F. A. Teodoro, D. J. Lawrence, R. C. Elphic, W. C. Feldman

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

Submitted 9 February, 2012; v1 submitted 9 August, 2011; originally announced August 2011.

Comments: 19 pages, 14 figures, updated version accepted by ApJ

Journal ref: Eke et al, 2012, ApJ, 747, 6

The Baryons in the Milky Way Satellites

Authors: Owen H. Parry, Vincent R. Eke, Carlos S. Frenk, Takashi Okamoto

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

Submitted 22 June, 2011; v1 submitted 17 May, 2011; originally announced May 2011.

Comments: 16 pages, 11 figures, submitted to MNRAS. Correction to calculation of central satellite masses (projected values were used for simulated satellites in original version)

The haloes of bright satellite galaxies in a warm dark matter universe

Authors: Mark Lovell, Vincent Eke, Carlos Frenk, Liang Gao, Adrian Jenkins, Tom Theuns, Jie Wang, Simon White, Alexey Boyarsky, Oleg Ruchayskiy

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

Submitted 16 January, 2012; v1 submitted 14 April, 2011; originally announced April 2011.

Comments: 8 pages, 7 figures, accepted by MNRAS. Text and figures updated

Report number: CERN-PH-TH-2011-172

Journal ref: Mon. Not. Roy. Astr. Soc., 420 (2012) 2318-2324

The Satellite Luminosity Functions of Galaxies in SDSS

Authors: Quan Guo, Shaun Cole, Vincent Eke, Carlos Frenk

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

Submitted 2 June, 2011; v1 submitted 13 January, 2011; originally announced January 2011.

Comments: 13 pages, 12 figures

The stellar and hot gas content of low-mass galaxy clusters

Authors: Michael L. Balogh, Pasquale Mazzotta, Richard G. Bower, Vince Eke, Herve Bourdin, Ting Lu, Tom Theuns

Abstract: We analyse the stellar and hot gas content of 18 nearby, low-mass galaxy clusters, detected in redshift space and selected to have a dynamical mass 3E14<M/Msun<6E14, as measured from the 2dF Galaxy Redshift Survey. We combine X-ray measurements from both Chandra and XMM with ground-based near-infrared observations from CTIO, AAT and CFHT to compare the mass in hot gas and stars to the dynamical ma… ▽ More We analyse the stellar and hot gas content of 18 nearby, low-mass galaxy clusters, detected in redshift space and selected to have a dynamical mass 3E14<M/Msun<6E14, as measured from the 2dF Galaxy Redshift Survey. We combine X-ray measurements from both Chandra and XMM with ground-based near-infrared observations from CTIO, AAT and CFHT to compare the mass in hot gas and stars to the dynamical mass and state of the clusters. Only 13 of the clusters are detected in X-ray emission, and for these systems we find that a range of 7-20 per cent of their baryonic mass, and <3 per cent of their dynamical mass, is detected in starlight, similar to what is observed in more massive clusters. In contrast, the five undetected clusters are underluminous in X-ray emission, by up to a factor 10, given their stellar mass. Although the velocity distribution of cluster members in these systems is indistinguishable from a Gaussian, all show subtle signs of being unrelaxed: either they lack a central, dominant galaxy, or the bright galaxy distribution is less concentrated and/or more elongated than the rest of the sample. Thus we conclude that low-mass clusters and groups selected from the velocity distribution of their galaxies exhibit a dichotomy in their hot gas properties. Either they are detected in X-ray, in which case they generally lie on the usual scaling relations, or they are completely undetected in X-ray emission. The non-detections may be partly related to the apparently young dynamical state of the clusters, but it remains a distinct possibility that some of these systems are exceptionally devoid of hot emitting gas as the result of its expulsion or rarefaction. △ Less

Submitted 13 December, 2010; v1 submitted 2 November, 2010; originally announced November 2010.

Comments: MNRAS, in press. Replacement to correct error caught in proofs. Two of the undetected clusters do not have Chandra or XMM data and the limits are based on ROSAT data

A Halo Expansion (HEX) Technique for Approximating Simulated Dark Matter Haloes

Authors: Ben Lowing, Adrian Jenkins, Vincent Eke, Carlos Frenk

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

Submitted 8 July, 2011; v1 submitted 29 October, 2010; originally announced October 2010.

Comments: 16 pages, 12 figures. Revised version with new title. Accepted for publication in MNRAS

Connected structure in the 2dFGRS

Authors: David N. A. Murphy, Vincent R. Eke, Carlos S. Frenk

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

Submitted 19 January, 2011; v1 submitted 11 October, 2010; originally announced October 2010.

Comments: 10 pages, 11 figures. Accepted for publication in MNRAS

The Link between Galactic Satellite Orbits and Subhalo Accretion

Authors: Mark Lovell, Vincent Eke, Carlos Frenk, Adrian Jenkins

Abstract: We calculate the orbital angular momentum of dark matter subhaloes in the Aquarius simulations of cold dark matter galactic haloes. We calculate the orientation of their angular momentum relative to that of the spin vector of their host halo and find a variety of different configurations. All six Aquarius haloes contain statistically significant populations of subhalo orbits that are aligned with… ▽ More We calculate the orbital angular momentum of dark matter subhaloes in the Aquarius simulations of cold dark matter galactic haloes. We calculate the orientation of their angular momentum relative to that of the spin vector of their host halo and find a variety of different configurations. All six Aquarius haloes contain statistically significant populations of subhalo orbits that are aligned with the main halo spin. All haloes posses a population of subhaloes that rotates in the same direction as the main halo and three of them possess, in addition, a population that rotates in the opposite direction. These configurations arise from the filamentary accretion of subhaloes. Quasi-planar distributions of coherently rotating satellites, such as those inferred in the Milky Way and other galaxies, arise naturally in simulations of a $Λ$CDM universe. △ Less

Submitted 18 April, 2011; v1 submitted 3 August, 2010; originally announced August 2010.

Comments: 10 pages, 8 figures, higher quality versions of some figures available at http://astro.dur.ac.uk/~d50wse/. Accepted by MNRAS

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

Authors: Robert A. Crain, Tom Theuns, Claudio Dalla Vecchia, Vincent R. Eke, Carlos S. Frenk, Adrian Jenkins, Scott T. Kay, John A. Peacock, Frazer R. Pearce, Joop Schaye, Volker Springel, Peter A. Thomas, Simon D. M. White, Robert P. C. Wiersma

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

Submitted 23 June, 2009; originally announced June 2009.

Comments: 23 pages, 14 figures, 2 tables, submitted to MNRAS. Comments welcome

The angular momentum of cold dark matter haloes with and without baryons

Authors: Philip Bett, Vincent Eke, Carlos S. Frenk, Adrian Jenkins, Takashi Okamoto

Abstract: 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 medi… ▽ More 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(<r) proportional to r. The direction of the vector varies considerably with radius: the median angle between the inner (< 0.25 Rvir) and total (< Rvir) angular momentum vectors is about 25degr. To investigate how baryons affect halo spin, we use another high resolution simulation, which includes gas cooling, star formation and feedback. This simulation produces a sample of galaxies with a realistic distribution of disc-to-total ratios. The formation of the galaxy spins up the dark matter within 0.1 Rvir such that the specific halo angular momentum increases by approx 50% in the median. The dark matter angular momentum becomes better aligned, but there remains a broad distribution of (mis-)alignments between the halo and the central galaxy, with a median angle between their angular momenta of ~ 30degr. Galaxies have a range of orientations relative to the shape of the halo: half of them have their minor axes misaligned by more than 45degr. Finally, we align a sample of haloes according to the orientation of their galaxies and stack the projected mass distributions. Although the individual haloes are aspherical, galaxy--halo misalignments produce a stacked mass distribution that is indistinguishable from circular. If the misalignment found in our simulations is realistic, it will be extremely difficult for weak lensing studies to measure the shape of dark haloes using this technique. △ Less

Submitted 18 January, 2010; v1 submitted 16 June, 2009; originally announced June 2009.

Comments: 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

Journal ref: MNRAS 2010, 404, p1137-1156

The spatial distribution of polar hydrogen deposits on the Moon

Authors: V. R. Eke, L. F. A. Teodoro, R. C. Elphic

Abstract: A new analysis of the Lunar Prospector epithermal neutron data is presented, providing an improved map of the distribution of hydrogen near to the lunar poles. This is achieved using a specially developed pixon image reconstruction algorithm to deconvolve the instrumental response of the Lunar Prospector's neutron spectrometer from the observed data, while simultaneously suppressing the statisti… ▽ More A new analysis of the Lunar Prospector epithermal neutron data is presented, providing an improved map of the distribution of hydrogen near to the lunar poles. This is achieved using a specially developed pixon image reconstruction algorithm to deconvolve the instrumental response of the Lunar Prospector's neutron spectrometer from the observed data, while simultaneously suppressing the statistical noise. The results show that these data alone require the hydrogen to be concentrated into the cold traps at up to 1 wt% water-equivalent hydrogen. This combination of localisation and high concentration suggests that the hydrogen is present either in the form of a volatile compound or as solar wind protons implanted into small regolith grains. △ Less

Submitted 14 October, 2008; originally announced October 2008.

Comments: 19 pages. 5 colour figures. Submitted to Icarus

Galaxy morphology in the LambdaCDM cosmology

Authors: O. H. Parry, V. R. Eke, C. S. Frenk

Abstract: We investigate the origins of galaxy morphology (defined by bulge-to-total K-band luminosity) in the LambdaCDM cosmology using two galaxy formation models ,based on the Millennium simulation, by Bower et al. (the Durham model) and De Lucia & Blaizot (the MPA model). There are many similarities, but also fundamental disagreements in the predictions of the two models for galaxy morphology. For exa… ▽ More We investigate the origins of galaxy morphology (defined by bulge-to-total K-band luminosity) in the LambdaCDM cosmology using two galaxy formation models ,based on the Millennium simulation, by Bower et al. (the Durham model) and De Lucia & Blaizot (the MPA model). There are many similarities, but also fundamental disagreements in the predictions of the two models for galaxy morphology. For example, taking into account uncertainties in the available observational data, both produce a realistic morphological mix today, but its evolution is very different. A main cause of this and other differences is the treatment of disk instabilities which play a more prominent role in the Durham model. Our analysis confirms previous theoretical predictions that elliptical galaxies form most of their stars before the bulk of the galaxy is assembled. Spirals tend to have later `assembly' times as a consequence of in-situ star formation. With the exception of the brightest ellipticals (stellar mass > 2.5e11 M_sun/h), we find that major mergers are not the primary mechanism by which most spheroids (ellipticals and spiral bulges) assemble their mass. In fact, the majority of ellipticals (and the overwhelming majority of spirals) never experience a major merger (above the resolution limit of our simulation.) Most ellipticals and spiral bulges acquire their stellar mass through minor mergers or disk instabilities. These conclusions are common to both models. The rotation properties of spheroids may help to constrain the importance of disk instabilities in these models. △ Less

Submitted 16 April, 2009; v1 submitted 26 June, 2008; originally announced June 2008.

Comments: Accepted for publication in MNRAS 15 pages, 8 figures Substantial revisions to original following referee's report

Testing Cold Dark Matter with the hierarchical buildup of stellar light

Authors: Michael L. Balogh, Ian G. McCarthy, Richard G. Bower, Vincent R. Eke

Abstract: (Abridged) We demonstrate that the tenet of hierarchical structure growth leads directly to a robust, falsifiable prediction for the correlation between stellar fraction (fstar) and total system mass (M500) of galaxy groups and clusters. This prediction is relatively insensitive to the details of baryonic physics or cosmological parameters. In particular, if the fstar-M500 relation is fixed and… ▽ More (Abridged) We demonstrate that the tenet of hierarchical structure growth leads directly to a robust, falsifiable prediction for the correlation between stellar fraction (fstar) and total system mass (M500) of galaxy groups and clusters. This prediction is relatively insensitive to the details of baryonic physics or cosmological parameters. In particular, if the fstar-M500 relation is fixed and does not evolve with redshift, CDM models predict the logarithmic slope of this relation to be b>-0.3. This constraint can be weakened if the fstar-M500 relation evolves strongly, but this implies more stars must be formed in situ in groups at low redshift. Conservatively requiring that at least half the stars in groups were formed by z=1, the constraint from evolution models is b>-0.35. Since the most massive clusters (M500=1E15 Msun) are observed to have fstar=0.01, this means that groups with M500=5E13 Msun must have fstar<0.03. Recent observations by Gonzalez et al. (2007) indicate a much steeper relation, with fstar>0.04 in groups, leading to b=-0.64. If confirmed, this would rule out hierarchical structure formation models: today's clusters could not have been built from today's groups, or even from the higher-redshift progenitors of those groups. We perform a careful analysis of these and other data to identify the most important systematic uncertainties in their measurements. Although correlated uncertainties on stellar and total masses might explain the steep observed relation, the data are only consistent with theory if the observed group masses are systematically underestimated. △ Less

Submitted 7 January, 2008; originally announced January 2008.

Comments: MNRAS, in press. 13 typeset pages

The baryon fraction of LambdaCDM haloes

Authors: Robert A. Crain, Vincent R. Eke, Carlos S. Frenk, Adrian Jenkins, Ian G. McCarthy, Julio F. Navarro, Frazer R. Pearce

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

Submitted 19 October, 2006; originally announced October 2006.

Comments: 9 pages, 5 figures (Figs 1. and 2 reduced in quality), 1 table, submitted to MNRAS. Version with high-resolution figures can be obtained from http://star-www.dur.ac.uk/~rcrain/baryonfractions/

Journal ref: Mon.Not.Roy.Astron.Soc.377:41-49,2007

The spin and shape of dark matter haloes in the Millennium simulation of a LambdaCDM universe

Authors: Philip Bett, Vincent Eke, Carlos S. Frenk, Adrian Jenkins, John Helly, Julio Navarro

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

Submitted 20 December, 2006; v1 submitted 28 August, 2006; originally announced August 2006.

Comments: 24 pages, 25 figures. Final version accepted for publication in MNRAS. Significant re-arranging and clarification of text from previous versions, science content unchanged. Final five figures at reduced resolution; full resolution version available at http://star-www.dur.ac.uk/~dph3peb1/spindistro_web

Journal ref: Mon.Not.Roy.Astron.Soc.376:215-232,2007

Galaxy Groups in the 2dFGRS: the number density of groups

Authors: V. R. Eke, C. M. Baugh, S. Cole, C. S. Frenk, J. F. Navarro

Abstract: The abundance of galaxy clusters as a function of mass is determined using the 2dFGRS Percolation-Inferred Galaxy Group (2PIGG) catalogue. This is used to estimate the amplitude of the matter fluctuation spectrum, parametrised by the linear theory rms density fluctuations in spheres of 8Mpc/h, sigma_8. The best-fitting value for this parameter is highly correlated with the mean matter density in… ▽ More The abundance of galaxy clusters as a function of mass is determined using the 2dFGRS Percolation-Inferred Galaxy Group (2PIGG) catalogue. This is used to estimate the amplitude of the matter fluctuation spectrum, parametrised by the linear theory rms density fluctuations in spheres of 8Mpc/h, sigma_8. The best-fitting value for this parameter is highly correlated with the mean matter density in the Universe, Omega_m, and is found to satisfy sigma_8=0.25Omega_m^{-0.92-4.5(Omega_m-0.22)^2}+/-10%(statistical)+/-20% (systematic) for 0.18<Omega_m<0.50, assuming that Omega_m+Omega_Lambda=1. This gives sigma_8=0.89 when Omega_m=0.25. A ~20% correction has been applied to undo the systematic bias inherent in the measurement procedure. Mock catalogues, constructed from large cosmological N-body simulations, are used to help understand and model these systematic errors. The abundance of galaxy groups as a function of group bj band luminosity is also determined. This is used in conjunction with the halo mass function, determined from simulations, to infer the variation of halo mass-to-light ratio over four orders of magnitude in halo mass. The mass-to-light ratio shows a minimum value of 100hM_sol/L_sol in the bj band at a total group luminosity of L_bJ ~ 5.10^9Lsol/h^2. Together with the observed Tully-Fisher relation, this implies that the observed rotation speed of Tully-Fisher galaxies is within ~10% of the typical circular speed of haloes hosting brightest galaxies of the same luminosity. △ Less

Submitted 21 October, 2005; originally announced October 2005.

Comments: Submitted to MNRAS, 12 pages, 17 figures (1 in colour)

Journal ref: Mon.Not.Roy.Astron.Soc.370:1147-1158,2006