Source author record

V. R. Eke

V. R. Eke appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

14 published item(s)

preprint2020arXiv

Atmospheric Erosion by Giant Impacts onto Terrestrial Planets

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.

preprint2020arXiv

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

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.

preprint2015arXiv

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

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.

preprint2015arXiv

The effect of craters on the lunar neutron flux

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

preprint2015arXiv

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

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.

preprint2013arXiv

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

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.

preprint2012arXiv

A quantitative comparison of lunar orbital neutron data

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.

preprint2009arXiv

Galaxy morphology in the LambdaCDM cosmology

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.

preprint2004arXiv

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

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.

preprint2004arXiv

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

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

preprint2004arXiv

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

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.

preprint2004arXiv

Where are the stars?

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.

preprint1998arXiv

Formation of Disc Galaxies

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

preprint1997arXiv

The Edinburgh-Durham Southern Galaxy Catalogue - VIII: The Cluster Galaxy Luminosity Function

We have re-examined the nature of the cluster galaxy luminosity function using the data from the Edinburgh-Durham Southern Galaxy Catalogue and the Edinburgh-Milano Redshift Survey. We derive a best fit luminosity function over the range M(bj)=-18 to -21, for a composite sample of 22 of the richer clusters that has M*=-20.16+/-0.02 and alpha=-1.22+/-0.04. The dominant error in these values results from the choice of background subtraction method. From extensive simulations we can show that when the LF is fitted over this narrow range, it is difficult to discriminate against bright values of M*in the single cluster fits, but that faint values provide a strong test of the universality of the luminosity function. We find that all the individual cluster data are well fit by a Schechter function with alpha fixed at -1.25, and that <=10% of these have fitted values of M* that disagree from the average at the 99% confidence level. We further show that fitting only a single parameter Schechter function to composite subsets of the data can give erroneous results for the derived M*. By considering two parameter fits, the results of Monte-Carlo simulations and direct two-sample chi-squared tests we conclude that there is only weak evidence for differences between the data when broken down into subsets based on physical properties (Bautz-Morgan class, richness, velocity dispersion): from our simulations, only the evidence for a difference between subsets based on velocity dispersion may in fact be significant. We find no evidence at all that a Schechter function is not a good model for the intrinsic cluster luminosity function over this absolute magnitude range. Models that invoke strong evolution of galaxy luminosity of all galaxies within clusters are inconsistent with our results.