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

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

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)

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

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

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

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

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

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

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

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)

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)

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

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

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

The effects of feedback on the morphology of galaxy discs

Authors: Takashi Okamoto, Vince R. Eke, Carlos S. Frenk, Adrian Jenkins

Abstract: We have performed hydrodynamic simulations of galaxy formation in a LCDM universe. We have followed galaxy formation in a dark matter halo, chosen to have a relatively quiet recent merger history, using different models for star formation and feedback. In all cases, we have adopted a multi-phase description of the interstellar medium and modelled star formation in quiescent and burst modes. We h… ▽ More We have performed hydrodynamic simulations of galaxy formation in a LCDM universe. We have followed galaxy formation in a dark matter halo, chosen to have a relatively quiet recent merger history, using different models for star formation and feedback. In all cases, we have adopted a multi-phase description of the interstellar medium and modelled star formation in quiescent and burst modes. We have explored two triggers for starbursts: strong shocks and high gas density, allowing for the possibility that stars in the burst may form with a top-heavy initial mass function. We find that the final morphology of the galaxy is extremely sensitive to the modelling of star formation and feedback. Starting from identical initial conditions, galaxies spanning the wide range of Hubble types, with B-band disc-to-total luminosity ratios ranging from 0.25 to 0.9, can form in the same dark matter halo. Models in which starbursts are induced by high gas density (qualitatively similar to models in which feedback is produced by AGN) generate energetic winds and result in galaxies with an early-type morphology. Models in which the starbursts are induced by strong shocks lead to extended discs. In this case, the feedback associated with the bursts suppresses the collapse of baryons in small halos, helping to create a reservoir of hot gas that is available for cooling after $z = 1$, following the bulk of the dynamical activity that builds up the halo. This gas then cools to form a large, young stellar disc. △ Less

Submitted 16 August, 2005; v1 submitted 30 March, 2005; originally announced March 2005.

Comments: 17 pages, 11 figures, accepted for Publication in MNRAS. Movies are available at http://star-www.dur.ac.uk/~takashi/research/movies.html

Journal ref: Mon.Not.Roy.Astron.Soc. 363 (2005) 1299-1314

Galaxy Groups at 0.3 <= z <= 0.55. II. Evolution to z ~ 0

Authors: D. J. Wilman, M. L. Balogh, R. G. Bower, J. S. Mulchaey, A. Oemler Jnr, R. G. Carlberg, V. R. Eke, I. Lewis, S. L. Morris, R. J. Whitaker

Abstract: We compare deep Magellan spectroscopy of 26 groups at 0.3 <= z <= 0.55, selected from the Canadian Network for Observational Cosmology 2 field survey (CNOC2), with a large sample of nearby groups from the 2PIGG catalogue (Eke et al., 2004). We find that the fraction of group galaxies with significant [OII] emission (>=5Å) increases strongly with redshift, from ~29% in 2dFGRS to ~58% in CNOC2, fo… ▽ More We compare deep Magellan spectroscopy of 26 groups at 0.3 <= z <= 0.55, selected from the Canadian Network for Observational Cosmology 2 field survey (CNOC2), with a large sample of nearby groups from the 2PIGG catalogue (Eke et al., 2004). We find that the fraction of group galaxies with significant [OII] emission (>=5Å) increases strongly with redshift, from ~29% in 2dFGRS to ~58% in CNOC2, for all galaxies brighter than ~ M*+1.75. This trend is parallel to the evolution of field galaxies, where the equivalent fraction of emission line galaxies increases from ~ 53% to ~ 75%. The fraction of emission-line galaxies in groups is lower than in the field, across the full redshift range, indicating that the history of star formation in groups is influenced by their environment. We show that the evolution required to explain the data is inconsistent with a quiescent model of galaxy evolution; instead, discrete events in which galaxies cease forming stars (truncation events) are required. We constrain the probability of truncation (P_trunc) and find that a high value is required in a simple evolutionary scenario neglecting galaxy mergers (P_trunc >~ 0.3 Gyr^{-1}). However, without assuming significant density evolution, P_trunc is not required to be larger in groups than in the field, suggesting that the environmental dependence of star formation was embedded at redshifts z >~ 0.45. △ Less

Submitted 11 January, 2005; originally announced January 2005.

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

Journal ref: Mon.Not.Roy.Astron.Soc. 358 (2005) 88-100

The 2dF Galaxy Redshift Survey: Power-spectrum analysis of the final dataset and cosmological implications

Authors: S. Cole, W. J. Percival, J. A. Peacock, P. Norberg, C. M. Baugh, C. S. Frenk, I. Baldry, J. Bland-Hawthorn, T. Bridges, R. Cannon, M. Colless, C. Collins, W. Couch, N. J. G. Cross, G. Dalton, V. R. Eke, R. De Propris, S. P. Driver, G. Efstathiou, R. S. Ellis, K. Glazebrook, C. Jackson, A. Jenkins, O. Lahav, I. Lewis , et al. (6 additional authors not shown)

Abstract: We present a power spectrum analysis of the final 2dF Galaxy Redshift Survey, employing a direct Fourier method. The sample used comprises 221,414 galaxies with measured redshifts. We investigate in detail the modelling of the sample selection. A new angular mask is derived, based on revisions to the photometric calibration. The redshift selection function is determined by dividing the survey ac… ▽ More We present a power spectrum analysis of the final 2dF Galaxy Redshift Survey, employing a direct Fourier method. The sample used comprises 221,414 galaxies with measured redshifts. We investigate in detail the modelling of the sample selection. A new angular mask is derived, based on revisions to the photometric calibration. The redshift selection function is determined by dividing the survey according to rest-frame colour, and deducing a self-consistent treatment of k-corrections and evolution for each population. The covariance matrix for the power-spectrum estimates is determined using two different approaches to the construction of mock surveys which are used to demonstrate that the input cosmological model can be correctly recovered. We are confident that the 2dFGRS power spectrum can be used to infer the matter content of the universe. On large scales, our estimated power spectrum shows evidence for the `baryon oscillations' that are predicted in CDM models. Fitting to a CDM model, assuming a primordial $n_{s}=1$ spectrum, $h=0.72$ and negligible neutrino mass, the preferred parameters are $Ω_{M} h = 0.168 \pm 0.016$ and a baryon fraction $Ω_{b} /Ω_{M} = 0.185\pm0.046$ (1$σ$ errors). The value of $Ω_{M} h$ is $1σ$ lower than the $0.20 \pm 0.03$ in our 2001 analysis of the partially complete 2dFGRS. This shift is largely due to the signal from the newly-sampled regions of space, rather than the refinements in the treatment of observational selection. This analysis therefore implies a density significantly below the standard $Ω_{M} =0.3$: in combination with CMB data from WMAP, we infer $Ω_{M} =0.231\pm 0.021$. (Abridged.) △ Less

Submitted 5 August, 2005; v1 submitted 11 January, 2005; originally announced January 2005.

Comments: Updated to match the in press MNRAS version. P(k) data and covariance matrix available at http://www.mso.anu.edu.au/2dFGRS/Public/Release/PowSpec

Journal ref: Mon.Not.Roy.Astron.Soc.362:505-534,2005

Where are the stars?

Authors: V. R. Eke, C. M. Baugh, S. Cole, C. S. Frenk, H. M. King, J. A. Peacock

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

Submitted 2 December, 2004; originally announced December 2004.

Comments: 15 pages, 14 figures (1 colour), submitted to MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc. 362 (2005) 1233-1246

The 2dF Galaxy Redshift Survey: The clustering of galaxy groups

Authors: Nelson D. Padilla, C. M. Baugh, V. R. Eke, P. Norberg, S. Cole, C. S. Frenk, D. J. Croton, I. K. Baldry, J. Bland-Hawthorn, T. Bridges, R. Cannon, M. Colless, C. Collins, W. Couch, G. Dalton, R. De Propris, S. P. Driver, G. Efstathiou, R. S. Ellis, K. Glazebrook, C. Jackson, O. Lahav, I. Lewis, S. Lumsden, S. Maddox , et al. (5 additional authors not shown)

Abstract: We measure the clustering of galaxy groups in the 2dFGRS Percolation-Inferred Galaxy Group (2PIGG) catalogue. The 2PIGG sample has 29,000 groups with at least two members. The clustering amplitude of the full 2PIGG catalogue is weaker than that of 2dFGRS galaxies, in agreement with theoretical predictions. We have subdivided the 2PIGG catalogue into samples that span a factor of 25 in median tot… ▽ More We measure the clustering of galaxy groups in the 2dFGRS Percolation-Inferred Galaxy Group (2PIGG) catalogue. The 2PIGG sample has 29,000 groups with at least two members. The clustering amplitude of the full 2PIGG catalogue is weaker than that of 2dFGRS galaxies, in agreement with theoretical predictions. We have subdivided the 2PIGG catalogue into samples that span a factor of 25 in median total luminosity. Our correlation function measurements span an unprecedented range of clustering strengths, connecting the regimes probed by groups fainter than L* galaxies and rich clusters. There is a steady increase in clustering strength with group luminosity; the most luminous groups are ten times more strongly clustered than the full 2PIGG catalogue. We demonstrate that the 2PIGG results are in very good agreement with the clustering of groups expected in the LCDM model. △ Less

Submitted 12 July, 2004; v1 submitted 25 February, 2004; originally announced February 2004.

Comments: 16 pages. Submitted to MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.352:211,2004

Galaxy groups in the 2dFGRS: the group-finding algorithm and the 2PIGG catalogue

Authors: V. R. Eke, C. M. Baugh, S. Cole, C. S. Frenk, P. Norberg, J. A. Peacock, I. K. Baldry, J. Bland-Hawthorn, T. Bridges, R. Cannon, M. Colless, C. Collins, W. Couch, G. Dalton, R. De Propris, S. P. Driver, G. Efstathiou, R. S. Ellis, K. Glazebrook, C. Jackson, O. Lahav, I. Lewis, S. Lumsden, S. Maddox, D. Madgwick , et al. (3 additional authors not shown)

Abstract: The construction of a catalogue of galaxy groups from the 2-degree Field Galaxy Redshift Survey (2dFGRS) is described. Groups are identified by means of a friends-of-friends percolation algorithm which has been thoroughly tested on mock versions of the 2dFGRS generated from cosmological N-body simulations. The tests suggest that the algorithm groups all galaxies that it should be grouping, with… ▽ More The construction of a catalogue of galaxy groups from the 2-degree Field Galaxy Redshift Survey (2dFGRS) is described. Groups are identified by means of a friends-of-friends percolation algorithm which has been thoroughly tested on mock versions of the 2dFGRS generated from cosmological N-body simulations. The tests suggest that the algorithm groups all galaxies that it should be grouping, with an additional 40% of interlopers. About 55% of the ~190000 galaxies considered are placed into groups containing at least two members of which ~29000 are found. Of these, ~7000 contain at least four galaxies, and these groups have a median redshift of 0.11 and a median velocity dispersion of 260km/s. This 2dFGRS Percolation-Inferred Galaxy Group (2PIGG) catalogue represents the largest available homogeneous sample of galaxy groups. It is publicly available on the WWW. △ Less

Submitted 26 February, 2004; v1 submitted 24 February, 2004; originally announced February 2004.

Comments: 15 pages, 8 figures (1 in colour), MNRAS, 348, 866 (as before)

Journal ref: Mon.Not.Roy.Astron.Soc.348:866,2004

Galaxy groups in the 2dFGRS: the luminous content of the groups

Authors: V. R. Eke, C. S. Frenk, C. M. Baugh, S. Cole, P. Norberg, J. A. Peacock, I. K. Baldry, J. Bland-Hawthorn, T. Bridges, R. Cannon, M. Colless, C. Collins, W. Couch, G. Dalton, R. De Propris, S. P. Driver, G. Efstathiou, R. S. Ellis, K. Glazebrook, C. Jackson, O. Lahav, I. Lewis, S. Lumsden, S. Maddox, D. Madgwick , et al. (3 additional authors not shown)

Abstract: The 2dFGRS Percolation-Inferred Galaxy Group (2PIGG) catalogue of ~29000 objects is used to study the luminous content of galaxy systems of various sizes. Mock galaxy catalogues constructed from cosmological simulations are used to gauge the accuracy with which intrinsic group properties can be recovered. A Schechter function is found to be a reasonable fit to the galaxy luminosity functions in… ▽ More The 2dFGRS Percolation-Inferred Galaxy Group (2PIGG) catalogue of ~29000 objects is used to study the luminous content of galaxy systems of various sizes. Mock galaxy catalogues constructed from cosmological simulations are used to gauge the accuracy with which intrinsic group properties can be recovered. A Schechter function is found to be a reasonable fit to the galaxy luminosity functions in groups of different mass in the real data. The characteristic luminosity L* is larger for more massive groups. However, the mock data show that the shape of the recovered luminosity function is expected to differ from the true shape, and this must be allowed for when interpreting the data. The variation of halo mass-to-light ratio with group size is studied in both these wavebands. A robust trend of increasing M/L with increasing group luminosity is found in the 2PIGG data. From groups with L_bj=10^{10}Lsol to those 100 times more luminous, the typical bj-band M/L increases by a factor of 5, whereas the rf-band M/L grows by a factor of 3.5. These trends agree well with the simulations, which also predict a minimum M/L on a scale corresponding to the Local Group. Our data indicate that if such a minimum exists, then it must occur at L<~10^{10}Lsol, below the range accurately probed by the 2PIGG catalogue. According to the mock data, the bj M/Ls of the largest groups are expected to be approximately 1.1 times the global value. Assuming that this correction applies to the real data yields an estimate of Omega_m=0.26+/-0.03 (statistical). △ Less

Submitted 26 February, 2004; v1 submitted 24 February, 2004; originally announced February 2004.

Comments: Same as before

Journal ref: Mon.Not.Roy.Astron.Soc.355:769,2004

Momentum transfer across shear flows in Smoothed Particle Hydrodynamic simulations of galaxy formation

Authors: T. Okamoto, A. Jenkins, V. R. Eke, V. Quilis, C. S. Frenk

Abstract: We investigate the evolution of angular momentum in SPH simulations of galaxy formation, paying particular attention to artificial numerical effects. We find that a cold gas disc forming in an ambient hot gas halo receives a strong hydrodynamic torque from the hot gas. By splitting the hydrodynamic force into artificial viscosity and pressure gradients, we find that the angular momentum transpor… ▽ More We investigate the evolution of angular momentum in SPH simulations of galaxy formation, paying particular attention to artificial numerical effects. We find that a cold gas disc forming in an ambient hot gas halo receives a strong hydrodynamic torque from the hot gas. By splitting the hydrodynamic force into artificial viscosity and pressure gradients, we find that the angular momentum transport is caused not by the artificial viscosity but by the pressure gradients. Using simple test simulations of shear flows, we conclude that the pressure gradient-based viscosity can be divided into two components: one due to the noisiness of SPH and the other to ram pressure. The former is problematic even with very high resolution because increasing resolution does not reduce the noisiness. On the other hand, the ram pressure effect appears only when a cold gas disc or sheet does not contain enough particles. In such a case, holes form in the disc or sheet, and then ram pressure from intra-hole hot gas, causes significant deceleration. In simulations of galactic disc formation, star formation usually decreases the number of cold gas particles, and hole formation leads to the fragmentation of the disc. To circumvent these problem, we modify the SPH algorithm, decoupling the cold from the hot gas phases, i.e. inhibiting the hydrodynamic interaction between cold and hot particles. This, a crude modelling of a multi-phase fluid in SPH cosmological simulations, leads to the formation of smooth extended cold gas discs and to better numerical convergence. The decoupling is applicable in so far as the self-gravitating gas disc with negligible external pressure is a good approximation for a cold gas disc. (abridged) △ Less

Submitted 26 June, 2003; originally announced June 2003.

Comments: 20 pages, 22 figures; accepted for publication in the MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.345:429,2003

On the nature of the ring-like structure in the outer Galactic disk

Authors: Amina Helmi, Julio F. Navarro, Andres Meza, Matthias Steinmetz, Vincent R. Eke

Abstract: We examine the tidal disruption of satellite galaxies in a cosmological simulation of the formation of a disk galaxy in the $Λ$CDM scenario. We find that the disruption of satellite galaxies in orbits roughly coplanar with the disk leads naturally to the formation of ring-like stellar structures similar to that recently discovered in the outer disk of the Milky Way. Two interpretations appear pl… ▽ More We examine the tidal disruption of satellite galaxies in a cosmological simulation of the formation of a disk galaxy in the $Λ$CDM scenario. We find that the disruption of satellite galaxies in orbits roughly coplanar with the disk leads naturally to the formation of ring-like stellar structures similar to that recently discovered in the outer disk of the Milky Way. Two interpretations appear plausible in this context. One is that the ring is a transitory, localized radial density enhancement reflecting the apocenter of particles stripped from a satellite during a recent pericentric passage (a ``tidal arc'' reminiscent of the tidal arms seen in disk galaxy mergers). In the second scenario, the ring is analogous to the ``shells'' found around some elliptical galaxies, and would result from a minor merger that took place several Gyr ago. The two interpretations differ in several ways. Tidal arcs are expected to span a limited longitude range; may carry a substantial fraction of the original mass of the satellite; should exhibit a significant velocity gradient with Galactic longitude; and are in all likelihood asymmetric above and below the plane of the disk. Shells, on the other hand, may comprise at most a small fraction of the original mass of the satellite and, due to their more relaxed state, ought to be symmetric above and below the plane, with no discernible velocity gradients across the structure. If confirmed as a tidal feature, the ring discovered by SDSS in the outer Galactic disk would strengthen the view--supported by numerical simulations--that minor mergers have played a significant role building up not only the stellar halo, but also the disk components of the Galaxy. △ Less

Submitted 16 June, 2003; v1 submitted 13 March, 2003; originally announced March 2003.

Comments: 4 pages, 2 figures, ApJ Letters in press. High resolution Fig. 2 available from http://www.astro.uu.nl/~helmi/ring. More thorough discussion of each scenario, and of observability of shells. Conclusions unchanged

Journal ref: Astrophys.J. 592 (2003) L25-L28

Simulations of Galaxy Formation in a Lambda CDM Universe III: The Dissipative Formation of an Elliptical Galaxy

Authors: Andres Meza, Julio F. Navarro, Matthias Steinmetz, Vincent R. Eke

Abstract: We examine in detail the dynamical structure of an elliptical galaxy simulated in the Lambda CDM scenario. The morphology of the galaxy evolves dramatically over time in response to the mode and timing of mass accretion; smooth deposition of cooled gas leads to the formation of centrifugally supported disks, whilst major mergers disperse stellar disks into spheroids. These two modes of accretion… ▽ More We examine in detail the dynamical structure of an elliptical galaxy simulated in the Lambda CDM scenario. The morphology of the galaxy evolves dramatically over time in response to the mode and timing of mass accretion; smooth deposition of cooled gas leads to the formation of centrifugally supported disks, whilst major mergers disperse stellar disks into spheroids. These two modes of accretion alternate successively until z~0.6, when the galaxy undergoes one last major (1:2) merger that consumes much of the remaining gas into stars. Little gas cools and accretes subsequently and, as a result, most stars at z=0 are in a spheroidal component that resembles present-day elliptical galaxies. Dynamically, the galaxy is well approximated by an E4 oblate rotator. Boxy isophotes are obtained when the galaxy is seen face-on and Vrot/sigma<<1. On the other hand, disky isophotes are found for inclinations which maximize Vrot/sigma. The line-of-sight velocity distribution is nearly Gaussian at all radii. The sign of the Gauss-Hermite skewness parameter h3 is anti-correlated with the apparent sense of rotation, in agreement with observed trends. The simulated galaxy has much higher effective surface brightness than normal ellipticals of similar luminosity, in a way reminiscent of the less common type of compact ellipticals. Our simulation shows that repeated episodes of dissipational collapse, followed by merger events, lead to stellar spheroids that are only mildly triaxial and of relatively simple kinematic structure. This is in better agreement with observation than earlier models based on dissipationless mergers of stellar disks, and a positive step towards reconciling the observed structure of ellipticals with the hierarchical merging scenario. △ Less

Submitted 20 August, 2003; v1 submitted 13 January, 2003; originally announced January 2003.

Comments: 17 pages, 20 color figures, replaced with published version ApJ

Journal ref: Astrophys.J. 590 (2003) 619-635

Simulations of Galaxy Formation in a Lambda CDM Universe II: The Fine Structure of Simulated Galactic Disks

Authors: Mario G. Abadi, Julio F. Navarro, Matthias Steinmetz, Vincent R. Eke

Abstract: We present a detailed analysis of the dynamical properties of a simulated disk galaxy assembled hierarchically in the Lambda CDM cosmogony. At z=0, two distinct dynamical components are identified on the basis of the orbital parameters of stars in the galaxy: a slowly rotating, centrally concentrated spheroid and a disk-like component largely supported by rotation. These components are also reco… ▽ More We present a detailed analysis of the dynamical properties of a simulated disk galaxy assembled hierarchically in the Lambda CDM cosmogony. At z=0, two distinct dynamical components are identified on the basis of the orbital parameters of stars in the galaxy: a slowly rotating, centrally concentrated spheroid and a disk-like component largely supported by rotation. These components are also recognized in the surface brightness profile of the galaxy, which can be very well approximated by the superposition of an R^{1/4} spheroid and an exponential disk. However, neither does the dynamically-identified spheroid follow de Vaucouleurs' law nor is the disk purely exponential, a resultwhich calls for caution when estimating the importance of the disk from traditional photometric decomposition techniques. The disk may be further decomposed into a thin, dynamically cold component with stars on nearly circularorbits and a hotter, thicker component with orbital parameters transitional between the thin disk and the spheroid. The spheroid is old, and has essentiallyno stars younger than the time elapsed since the last major accretion event ~8 Gyr ago. The majority of thin disk stars, form after the merging activity is over, although a significant fraction ~15% of thin-disk stars are old enough to predate the last major merger event. This unexpected population of old disk stars consists mainly of the tidal debris of satellites whose orbital plane was coincident with the disk and whose orbits were circularized by dynamical friction prior to full disruption. Our results highlight the role of satellite accretion events in shaping the disk and the spheroidal component and reveal some of the clues to the assembly process of a galaxy preserved in the detailed dynamics of old stellar populations. △ Less

Submitted 11 December, 2002; originally announced December 2002.

Comments: 15 pages, 14 figures, submitted to ApJ

Journal ref: Astrophys.J.597:21-34,2003

Simulations of Galaxy Formation in a Lambda CDM Universe I: Dynamical and Photometric Properties of a Simulated Disk Galaxy

Authors: Mario G. Abadi, Julio F. Navarro, Matthias Steinmetz, Vincent R. Eke

Abstract: We analyze the properties of a disk galaxy simulated with unprecedented numerical resolution in the Lambda CDM cosmogony. The galaxy is assembled through a number of high-redshift mergers followed by a period of quiescent accretion after z~1 which lead to the formation of two distinct dynamical components: a spheroid and a disk. The surface brightness profile is very well approximated by the sup… ▽ More We analyze the properties of a disk galaxy simulated with unprecedented numerical resolution in the Lambda CDM cosmogony. The galaxy is assembled through a number of high-redshift mergers followed by a period of quiescent accretion after z~1 which lead to the formation of two distinct dynamical components: a spheroid and a disk. The surface brightness profile is very well approximated by the superposition of an R^{1/4} spheroid and an exponential disk. The surface brightness profile is remarkably similar to that of Sab galaxy UGC615, but the simulated galaxy rotates significantly faster and has a declining rotation curve dominated by the spheroid near the center. The decline in circular velocity is at odds with observation and results from the high concentration of the dark matter and baryonic components, as well as from the relatively high mass-to-light ratio of the stars in the simulation. The simulated galaxy lies ~ 1 mag off the I-band Tully-Fisher relation of late-type spirals, but seems to be in reasonable agreement with Tully-Fisher data on S0 galaxies. The angular momentum of the luminous component is an order of magnitude lower than that of late-type spirals of similar rotation speed. This reflects the dominance of the slowly-rotating, dense spheroidal component. The disk component, on the other hand, has properties rather similar to those of late-type spirals. This suggests that a different form of feedback than adopted in this simulation is required to inhibit the efficient collapse and cooling of gas at high redshift that leads to the formation of the spheroid. Reconciling the properties of disk galaxies with the early collapse and high merging rates characteristic of hierarchical scenarios such as Lambda CDM remains a challenging, yet so far elusive, proposition. △ Less

Submitted 14 November, 2002; originally announced November 2002.

Comments: 15 pages, 13 figures, submitted to ApJ

Journal ref: Astrophys.J.591:499-514,2003

Measuring Sigma8 with Cluster Lensing : Biases from Unrelaxed Clusters

Authors: Graham P. Smith, Alastair C. Edge, Vincent R. Eke, Robert C. Nichol, Ian Smail, Jean-Paul Kneib

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

Submitted 8 November, 2002; originally announced November 2002.

Comments: Submitted to Astrophysical Journal Letters, 4 pages, 2 figures

Journal ref: Astrophys.J. 590 (2003) L79-L82

Surveying the Sky with the Arcminute MicroKelvin Imager: Expected Constraints on Galaxy Cluster Evolution and Cosmology

Authors: Ruediger Kneissl, Michael E. Jones, Richard Saunders, Vincent R. Eke, Anthony N. Lasenby, Keith Grainge, Garret Cotter

Abstract: We discuss prospects for cluster detection via the Sunyaev-Zel'dovich (SZ) effect in a blank field survey with the proposed interferometer array, the Arcminute MicroKelvin Imager (AMI). Clusters of galaxies selected in the SZ effect probe cosmology and structure formation with little observational bias, because the effect directly measures integrated gas pressure, and does so independently of cl… ▽ More We discuss prospects for cluster detection via the Sunyaev-Zel'dovich (SZ) effect in a blank field survey with the proposed interferometer array, the Arcminute MicroKelvin Imager (AMI). Clusters of galaxies selected in the SZ effect probe cosmology and structure formation with little observational bias, because the effect directly measures integrated gas pressure, and does so independently of cluster redshift. We use hydrodynamical simulations in combination with the Press-Schechter expression to simulate SZ cluster sky maps. These are used with simulations of the observation process to gauge the expected SZ cluster counts. Even with a very conservative choice of parameters we find that AMI will discover at least several tens of clusters every year with total mass down to 10^14 solar masses; the numbers depends on factors such as the mean matter density, the density fluctuation power spectrum and cluster gas evolution. The AMI survey itself can distinguish between these to some degree, and parameter degeneracies are largely eliminated given optical and X-ray follow-up of these clusters; this will also permit direct investigation of cluster physics and what drives the evolution. △ Less

Submitted 6 March, 2001; v1 submitted 2 March, 2001; originally announced March 2001.

Comments: 13 pages, 14 figures, submitted to MNRAS. A version with higher resolution images is available at http://www.mrao.cam.ac.uk/telescopes/ami; minor typos corrected

Journal ref: Mon.Not.Roy.Astron.Soc. 328 (2001) 783

The Power Spectrum Dependence of Dark Matter Halo Concentrations

Authors: V. R. Eke, J. F. Navarro, M. Steinmetz

Abstract: High-resolution N-body simulations are used to examine the power spectrum dependence of the concentration of galaxy-sized dark matter halos. It is found that dark halo concentrations depend on the amplitude of mass fluctuations as well as on the ratio of power between small and virial mass scales. This finding is consistent with the original results of Navarro, Frenk & White (NFW), and allows th… ▽ More High-resolution N-body simulations are used to examine the power spectrum dependence of the concentration of galaxy-sized dark matter halos. It is found that dark halo concentrations depend on the amplitude of mass fluctuations as well as on the ratio of power between small and virial mass scales. This finding is consistent with the original results of Navarro, Frenk & White (NFW), and allows their model to be extended to include power spectra substantially different from Cold Dark Matter (CDM). In particular, the single-parameter model presented here fits the concentration dependence on halo mass for truncated power spectra, such as those expected in the warm dark matter scenario, and predicts a stronger redshift dependence for the concentration of CDM halos than proposed by NFW. The latter conclusion confirms recent suggestions by Bullock et al., although this new modeling differs from theirs in detail. These findings imply that observational limits on the concentration, such as those provided by estimates of the dark matter content within individual galaxies, may be used to constrain the amplitude of mass fluctuations on galactic and subgalactic scales. The constraints on $Λ$CDM models posed by the dark mass within the solar circle in the Milky Way and by the zero-point of the Tully-Fisher relation are revisited, with the result that neither dataset is clearly incompatible with the `concordance' ($Ω_0=0.3$, $Λ_0=0.7$, $σ_8=0.9$) $Λ$CDM cosmogony. This conclusion differs from that reached recently by Navarro & Steinmetz, a disagreement that can be traced to inconsistencies in the normalization of the $Λ$CDM power spectrum used in that work. △ Less

Submitted 14 December, 2000; originally announced December 2000.

Comments: 12 pages including 7 figures using emulateapj, submitted to ApJ

Journal ref: Astrophys.J.554:114-125,2001

Cosmological parameters from cluster abundances, CMB and IRAS

Authors: S. L. Bridle, V. R. Eke, O. Lahav, A. N. Lasenby, M. P. Hobson, S. Cole, C. S. Frenk, J. P. Henry

Abstract: We combine information on cosmological parameters from cluster abundances, CMB primordial anisotropies and the IRAS 1.2 Jy galaxy redshift survey. We take as free parameters the present values of the total matter density of the universe, Omega_m, the Hubble parameter, h, sigma_8, and the IRAS biasing factor, b_IRAS. We assume that the universe is spatially flat, with a cosmological constant, and… ▽ More We combine information on cosmological parameters from cluster abundances, CMB primordial anisotropies and the IRAS 1.2 Jy galaxy redshift survey. We take as free parameters the present values of the total matter density of the universe, Omega_m, the Hubble parameter, h, sigma_8, and the IRAS biasing factor, b_IRAS. We assume that the universe is spatially flat, with a cosmological constant, and that structure formed from adiabatic initial fluctuations with a Harrison-Zel'dovich power spectrum. The nucleosynthesis value Omega_b=0.019/h^2 is adopted. We use the full three- and four- dimensional likelihood functions for each data set and marginalise these to two- and one- dimensional distributions in a Bayesian way, integrating over the other parameters. It is shown that the three data sets are in excellent agreement, with a best fit point of Omega_m= 1-Omega_Lambda=0.36, h=0.54, sigma_8=0.74, and b_IRAS=1.08. This point is within one sigma of the minimum for each data set alone. Pairs of these data sets have their degeneracies in sufficiently different directions that using only two data sets at a time is sufficient to place good constraints on the cosmological parameters. We show that the results from each of the three possible pairings of the data are also in good agreement. Finally, we combine all three data sets to obtain marginalised 68 per cent confidence intervals of 0.30<Omega_m<0.43, 0.48<h<0.59, 0.69<sigma_8<0.79, and 1.01<b_IRAS<1.16. For the best fit parameters the CMB quadrupole is Q_rms-ps=18.0 uK, Gamma=0.15, Omega_b=0.066 and the age of the universe is 16.7 Gyr. △ Less

Submitted 31 March, 1999; originally announced March 1999.

Comments: 6 pages, 4 figures. Submitted to MNRAS

Measuring Omega_0 using cluster evolution

Authors: V. R. Eke, S. Cole, C. S. Frenk, J. P. Henry

Abstract: The evolution of the galaxy cluster abundance depends sensitively on the value of the cosmological density parameter, Omega_0. Recent ASCA data are used to quantify this evolution as measured by the X-ray temperature function. A chi^2 minimisation fit to the cumulative temperature function, as well as a maximum likelihood estimate (which requires additional assumptions about cluster luminosities… ▽ More The evolution of the galaxy cluster abundance depends sensitively on the value of the cosmological density parameter, Omega_0. Recent ASCA data are used to quantify this evolution as measured by the X-ray temperature function. A chi^2 minimisation fit to the cumulative temperature function, as well as a maximum likelihood estimate (which requires additional assumptions about cluster luminosities), lead to the estimate Omega_0 \approx 0.45+/-0.2 (1-sigma statistical error). Various systematic uncertainties are considered, none of which enhance significantly the probability that Omega_0=1. These conclusions hold for models with or without a cosmological constant. The statistical uncertainties are at least as large as the individual systematic errors that have been considered here, suggesting that additional temperature measurements of distant clusters will allow an improvement in this estimate. An alternative method that uses the highest redshift clusters to place an upper limit on Omega_0 is also presented and tentatively applied, with the result that Omega_0=1 can be ruled out at the 98 per cent confidence level. Whilst this method does not require a well-defined statistical sample of distant clusters, there are still modelling uncertainties that preclude a firmer conclusion at this time. △ Less

Submitted 27 February, 1998; originally announced February 1998.

Comments: 17 pages, 15 figures, submitted to MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc. 298 (1998) 1145