Normalized to: Teodoro, L.

**37** article(s) in total. **81** co-authors, from **1** to **14** common article(s). Median position in authors list is **3,0**.

[1] oai:arXiv.org:1901.09934 [pdf] - 2076543

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

Comments: Fixed typo in Appx B equation and updated urls

Submitted: **2019-01-28**, last modified: **2020-04-07**

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.

[2] oai:arXiv.org:2003.14406 [pdf] - 2072865

Physical Properties of 299 NEOs Manually Recovered in Over Five Years of
NEOWISE Survey Data

Comments: Accepted for publication to the AAS Planetary Science Journal

Submitted: **2020-03-31**

Thermal infrared measurements of near-Earth objects provide critical data for
constraining their physical properties such as size. The NEOWISE mission has
been conducting an all-sky infrared survey to gather such data and improve our
understanding of this population. While automated routines are employed to
identify the majority of moving objects detected by NEOWISE, a subset of
objects will have dynamical properties that fall outside the window detectable
to these routines. Using the population of known near-Earth objects, we have
conducted a manual search for detections of these objects that were previously
unreported. We report 303 new epochs of observations for 299 unique near-Earth
objects of which 239 have no previous physical property characterization from
the NEOWISE Reactivation mission. As these objects are drawn from a list with
inherent optical selection biases, the distribution of measured albedos is
skewed to higher values than is seen for the diameter-selected population
detected by the automated routines. These results demonstrate the importance
and benefit of periodic searches of the archival NEOWISE data.

[3] oai:arXiv.org:2002.02977 [pdf] - 2045531

Atmospheric Erosion by Giant Impacts onto Terrestrial Planets

Comments: 17 pages, 13 figures, submitted to ApJ

Submitted: **2020-02-07**

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 for the first time. 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.

[4] oai:arXiv.org:1803.07083 [pdf] - 1710628

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

Comments: Published in ApJ. 14 pages, 11 figures

Submitted: **2018-03-19**, last modified: **2018-07-03**

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.

[5] oai:arXiv.org:1802.09508 [pdf] - 1739881

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

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

Submitted: **2018-02-26**

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.

[6] oai:arXiv.org:1612.02414 [pdf] - 1580830

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

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

Submitted: **2016-12-07**, last modified: **2017-10-31**

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.

[7] oai:arXiv.org:1708.00518 [pdf] - 1586653

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

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

Submitted: **2017-08-01**

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.

[8] oai:arXiv.org:1611.05395 [pdf] - 1532665

How thick are Mercury's polar water ice deposits?

Comments: 10 pages, 6 figures, submitted to Icarus

Submitted: **2016-11-16**

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.

[9] oai:arXiv.org:1505.05403 [pdf] - 1263029

The effect of craters on the lunar neutron flux

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

Submitted: **2015-05-20**, last modified: **2015-08-17**

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

[10] oai:arXiv.org:1501.05712 [pdf] - 925671

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

Comments: 8 pages, 7 figures and 3 tables

Submitted: **2015-01-22**

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.

[11] oai:arXiv.org:1409.1209 [pdf] - 916752

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

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

Submitted: **2014-09-03**, last modified: **2015-01-05**

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.

[12] oai:arXiv.org:1210.1203 [pdf] - 900738

Reconstructing the velocity field beyond the local universe

Comments: Accepted for publication in General Relativity and Gravitation (GRG)

Submitted: **2012-10-03**, last modified: **2014-11-17**

We present a maximum probability approach to reconstructing spatial maps of
the peculiar velocity field at redshifts $z\sim0.1$, where the velocities have
been measured from distance indicators (DI) such as $D_n-\sigma$ relations or
Tully-Fisher. With the large statistical uncertainties associated with DIs, our
reconstruction method aims to recover the underlying true peculiar velocity
field by reducing these errors with the use of two physically motivated
filtering prior terms. The first constructs an estimate of the velocity field
derived from the galaxy over-density {\deltag} and the second makes use of the
matter linear density power spectrum {\pk}. Using $N$-body simulations we find,
with an SDSS-like sample ($N_{gal}\simeq33$ per deg$^2$), an average
correlation coefficient value of $r=0.55\pm{0.02}$ between our reconstructed
velocity field and that of the true velocity field from the simulation.
However, with a suitably high number density of galaxies from the next
generation surveys (e.g. $N_{gal}\simeq140$ per deg$^2$) we can achieve an
average $r=0.70\pm{0.02}$ out to moderate redshifts $z\sim0.1$. This will prove
useful for future tests of gravity, as these relatively deep maps are
complementary to weak lensing maps at the same redshift.

[13] oai:arXiv.org:1312.4749 [pdf] - 862457

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

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

Submitted: **2013-12-17**, last modified: **2014-08-25**

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.

[14] oai:arXiv.org:1304.8123 [pdf] - 1166327

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

Comments: 12 pages, 13 figures

Submitted: **2013-04-30**

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.

[15] oai:arXiv.org:1108.2048 [pdf] - 474715

A quantitative comparison of lunar orbital neutron data

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

Submitted: **2011-08-09**, last modified: **2012-02-09**

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.

[16] oai:arXiv.org:1107.0858 [pdf] - 1077764

Completeness III: identifying characteristic systematics and evolution
in galaxy redshift surveys

Comments: 18 pages, 12 figures, 1 table, Accepted version for publication in
MNRAS

Submitted: **2011-07-05**, last modified: **2011-12-11**

This paper continues our development of non-parametric tests for analysing
the completeness in apparent magnitude of magnitude-redshift surveys. The
purpose of this third and final paper in our completeness series is two-fold:
firstly we explore how certain forms of incompleteness for a given flux-limited
galaxy redshift survey would manifest themselves in the ROBUST Tc and Tv
completeness estimators introduced in our earlier papers; secondly we provide a
comprehensive error propagation for these estimators. By using both real
surveys and Monte Carlo mock survey data, we have found distinct,
characteristic behaviour of the Tc and Tv estimators which identify
incompleteness in the form of e.g. missing objects within a particular
magnitude range. Conversely we have identified signatures of `over'
completeness, in cases where a survey contains a small region in apparent
magnitude that may have too many objects relative to the rest of the data set.
We also demonstrate how incompleteness resulting from luminosity evolution can
be identified and provide a framework for using our estimators as a robust tool
for constraining models of luminosity evolution.
Finally we explore the error propagation for Tc andTv. This builds on
Completeness II by allowing the definition of these estimators, and their
errors, via an adaptive procedure that accounts for the effects of sampling
error on the observed distribution of apparent magnitude and redshift in a
survey.

[17] oai:arXiv.org:1005.3888 [pdf] - 169766

Completeness II: A signal-to-noise approach for completeness estimators
applied to galaxy magnitude-redshift surveys

Comments: 9 pages, 6 figures, MNRAS accepted

Submitted: **2010-05-21**

This is the second paper in our completeness series which addresses some of
the issues raised in the previous article by Johnston, Teodoro & Hendry (2007)
in which we developed statistical tests for assessing the completeness in
apparent magnitude of magnitude-redshift surveys defined by two flux limits.
The statistics, Tc and Tv, associated with these tests are non-parametric and
defined in terms of the observed cumulative distribution function of sources;
they represent powerful tools for identifying the true flux limit and/or
characterising systematic errors in magnitude-redshift data. In this paper we
present a new approach to constructing these estimators that resembles an
"adaptive smoothing" procedure - i.e. by seeking to maintain the same amount
the information, as measured by the signal-to-noise ratio, allocated to each
galaxy. For consistency with our previous work, we apply our improved
estimators to the Millennium Galaxy Catalogue (MGC) and the Two Degree Field
Galaxy Redshift Survey (2dFGRS) data, and demonstrate that one needs to use a
signal-to-noise appropriately tailored for each individual catalogue to
optimise the performance of the completeness estimators. Furthermore, unless
such an adaptive procedure is employed, the assessment of completeness may
result in a spurious outcome if one uses other estimators present in the
literature which have not been designed taking into account "shot noise" due to
sampling.

[18] oai:arXiv.org:0909.3581 [pdf] - 1017604

On the surface extraction of electrons in a pulsar

Comments: 8 pages, 1 figure, accepted for publication in MNRAS

Submitted: **2009-09-19**

We present a novel description of how energetic electrons may be ejected from
the pulsar interior into the atmosphere, based on the collective electrostatic
oscillations of interior electrons confined to move parallel to the magnetic
field. The size of the interior magnetic field influences the interior plasma
frequency, via the associated matter density compression. The plasma
oscillations occur close to the regions of maximum magnetic field curvature,
that is, close to the magnetic poles where the majority of magnetic flux
emerges. Given that these oscillations have a density-dependent maximum
amplitude before wave-breaking occurs, such waves can eject energetic electrons
using only the self-field of the electron population in the interior. Moreover,
photons emitted by electrons in the bulk of the oscillation can escape along
the field lines by virtue of the lower opacity there (and the fact that they
are emitted predominantly in this direction), leading to features in the
spectra of pulsars.

[19] oai:arXiv.org:0802.1935 [pdf] - 10063

Peculiar Velocities into the Next Generation: Cosmological Parameters
From Large Surveys without Bias from Nonlinear Structure

Comments: Accepted for publication in MNRAS, 12 pages, 5 figures. V2 Discussion
clarified, 1 figure added, improvements to the text and figures; V3 Figure 5
Plotting error corrected, SN1a contours smaller

Submitted: **2008-02-13**, last modified: **2008-11-03**

We investigate methods to best estimate the normalisation of the mass density
fluctuation power spectrum (sigma_8) using peculiar velocity data from a survey
like the Six degree Field Galaxy Velocity Survey (6dFGSv). We focus on two
potential problems (i) biases from nonlinear growth of structure and (ii) the
large number of velocities in the survey. Simulations of LambdaCDM-like models
are used to test the methods. We calculate the likelihood from a full
covariance matrix of velocities averaged in grid cells. This simultaneously
reduces the number of data points and smooths out nonlinearities which tend to
dominate on small scales. We show how the averaging can be taken into account
in the predictions in a practical way, and show the effect of the choice of
cell size. We find that a cell size can be chosen that significantly reduces
the nonlinearities without significantly increasing the error bars on
cosmological parameters. We compare our results with those from a principal
components analysis following Watkins et al (2002) and Feldman et al (2003) to
select a set of optimal moments constructed from linear combinations of the
peculiar velocities that are least sensitive to the nonlinear scales. We
conclude that averaging in grid cells performs equally well. We find that for a
survey such as 6dFGSv we can estimate sigma_8 with less than 3% bias from
nonlinearities. The expected error on sigma_8 after marginalising over Omega_m
is approximately 16 percent.

[20] oai:arXiv.org:0810.2478 [pdf] - 17418

The spatial distribution of polar hydrogen deposits on the Moon

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

Submitted: **2008-10-14**

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.

[21] oai:arXiv.org:0810.1762 [pdf] - 17291

Drag-Gravity torques on galaxies in clusters: radial small-scale
alignment effects

Comments: 6 pages

Submitted: **2008-10-09**

We calculate the torque on galaxies in clusters due to gravity and to
dynamical friction forces in order to study the possible origin of small-scale
alignment effects as the result of interactions with their environment. The
equation of motion for the position angle of a galaxy is derived by using a ple
model. We find that weak radial alignment effects can be produced by this
mechanism involving only the most massive galaxies. We also introduce a
dependence on the cluster eccentricity to our equations in order to explore the
alignment of galaxies with the cluster's major axis. We find that in the inner
regions of high eccentricity clusters, alignments of massive galaxies with the
cluster's major axis dominate over the radial ones. This mechanism could
account for the observed alignment effects of the most massive galaxies with
the major axis of their host cluster. Our results suggest that dynamical
friction is a viable generator of alignment only for the most massive cluster
galaxies. For the observed alignments of normal galaxies a primordial origin
has to be explored.

[22] oai:arXiv.org:0805.0842 [pdf] - 12379

Effects of correlation between merging steps on the global halo
formation

Comments: 8 pages, submitted to MNRAS

Submitted: **2008-05-07**

The excursion set theory of halo formation is modified by adopting the
fractional Brownian motion, to account for possible correlation between merging
steps. We worked out analytically the conditional mass function, halo merging
rate and formation time distribution in the spherical collapse model. We also
developed an approximation for the ellipsoidal collapse model and applied it to
the calculation of the conditional mass function and the halo formation time
distribution. For models in which the steps are positively correlated, the halo
merger rate is enhanced when the accreted mass is less than $\sim 25M^*$, while
for the negatively correlated case this rate is reduced. Compared with the
standard model in which the steps are uncorrelated, the models with positively
correlated steps produce more aged population in small mass halos and more
younger population in large mass halos, while for the models with negatively
correlated steps the opposite is true. An examination of simulation results
shows that a weakly positive correlation between successive merging steps
appears to fit best. We have also found a systematic effect in the measured
mass function due to the finite volume of simulations. In future work, this
will be included in the halo model to accurately predict the three point
correlation function estimated from simulations.

[23] oai:arXiv.org:0803.3212 [pdf] - 260177

Resolving the universe with multifractals

Comments: To appear in the refereed Proceedings `A Century of Cosmology: Past,
Present and Future', San Servolo (Venezia, Italy), August 2007, to be
published by `Il Nuovo Cimento'. Extended version - 7 pages, 4 Figures

Submitted: **2008-03-21**

We present a new method for dealing with geometrical selection effects in
galaxy surveys while using a multifractal framework. The power of multifractal
analysis lies in its connection to higher order moments, in that it not only
probes clustering on different scales but also different densities. Therefore
any incompleteness issues must be correctly addressed before blindly applying
this technique to real survey data.

[24] oai:arXiv.org:0705.4632 [pdf] - 1773

Unplugging the Universe: the neglected electromagnetic consequence of
decoupling

Comments: 4 pages, no figures

Submitted: **2007-05-31**, last modified: **2007-06-04**

This letter concentrates on the non-equilibrium evolution of magnetic field
structures at the onset of recombination, when the charged particle current
densities decay as neutrals are formed.
We consider the effect that a decaying magnetic flux has on the acceleration
of particles via the transient induced electric field. Since the residual
charged-particle number density is small as a result of decoupling, we shall
consider the magnetic and electric fields essentially to be imposed, neglecting
the feedback from any minority accelerated population.
We find that the electromagnetic treatment of this phase transition can
produce energetic electrons scattered throughout the Universe. Such particles
could have a significant effect on cosmic evolution in several ways: (i) their
presence could delay the effective end of the recombination era; (ii) they
could give rise to plasma concentrations that could enhance early gravitational
collapse of matter by opposing cosmic expansion to a greater degree than
neutral matter could; (iii) they could continue to be accelerated, and become
the seed for reionisation at the later epoch $z \approx 10$.

[25] oai:arXiv.org:astro-ph/0703040 [pdf] - 89849

Completeness - I. Revised, reviewed and revived

Comments: 12 pages, 7 figures. MNRAS accepted

Submitted: **2007-03-02**

We have extended and improved the statistical test recently developed by
Rauzy for assessing the completeness in apparent magnitude of
magnitude-redshift surveys. Our improved test statistic retains the robust
properties -- specifically independence of the spatial distribution of galaxies
within a survey -- of the $T_c$ statistic introduced in Rauzy's seminal paper,
but now accounts for the presence of both a faint and bright apparent magnitude
limit. We demonstrate that a failure to include a bright magnitude limit can
significantly affect the performance of Rauzy's $T_c$ statistic. Moreover, we
have also introduced a new test statistic, $T_v$, defined in terms of the
cumulative distance distribution of galaxies within a redshift survey. These
test statistics represent powerful tools for identifying and characterising
systematic errors in magnitude-redshift data. We apply our completeness test to
three major redshift surveys: The Millennium Galaxy Catalogue (MGC), The Two
Degree Field Galaxy Redshift Survey (2dFGRS), and the Sloan Digital Sky Survey
(SDSS). We confirm that MGC and SDSS are complete up to the published (faint)
apparent magnitude limit Furthermore, we show that, unless a bright limit is
included for 2dFGRS, the data-set displays significant incompleteness at
magnitudes brighter than the published limit of $m_{b_{j}}=19.45$ mag.

[26] oai:arXiv.org:astro-ph/0702729 [pdf] - 89783

A Cautionary Note on Cosmological Magnetic Fields

Comments: 5 pages

Submitted: **2007-02-27**

This note is concerned with potentially misleading concepts in the treatment
of cosmological magnetic fields by magnetohydrodynamical (MHD) modelling. It is
not a criticism of MHD itself but rather a cautionary comment on the validity
of its use in cosmology. Now that cosmological data are greatly improved
compared with a few decades ago, and even better data are imminent, it makes
sense to revisit original modelling assumptions and examine critically their
shortcomings in respect of modern science. Specifically this article argues
that ideal MHD is a poor approximation around recombination, since it
inherently restricts evolutionary timescales, and is often misapplied in the
existing literature.

[27] oai:arXiv.org:astro-ph/0506395 [pdf] - 316293

Precision Determination of the Mass Function of Dark Matter Halos

Comments: 4 pages, 3 figures, submitted to ApJL

Submitted: **2005-06-16**

The predicted mass function of dark matter halos is essential in connecting
observed galaxy cluster counts and models of galaxy clustering to the
properties of the primordial density field. We determine the mass function in
the concordance $\Lambda$CDM cosmology, as well as its uncertainty, using
sixteen $1024^3$-particle nested-volume dark-matter simulations, spanning a
mass range of over five orders of magnitude. Using the nested volumes and
single-halo tests, we find and correct for a systematic error in the
friends-of-friends halo-finding algorithm. We find a fitting form and full
error covariance for the mass function that successfully describes the
simulations' mass function and is well-behaved outside the simulations'
resolutions. Estimated forecasts of uncertainty in cosmological parameters from
future cluster count surveys have negligible contribution from remaining
statistical uncertainties in the central cosmology multiplicity function. There
exists a potentially non-negligible cosmological dependence (non-universality)
of the halo multiplicity function.

[28] oai:arXiv.org:astro-ph/0406174 [pdf] - 65331

The Distribution of Dark Matter in the Universe on Scales of 10^10 M_sun
to 10^15 M_sun

Comments: 4 pages, 5 figures, talk delivered at the Rencontres de Moriond -
Cosmology: Exploring the Universe (La Thuile, 28 March - 4 April, 2004), to
appear in the proceedings

Submitted: **2004-06-07**

The use of parallel computers and increasingly sophisticated software has
allowed us to perform a large suite of N-body simulations using from $10^8$,to
$10^9$ particles. We will report on our recent convergence tests of the halo
mass function, N-point correlation functions, power spectrum and pairwise
velocity from very large high resolution treecode N-body simulations. Rather
than basing results on just one or two large simulations, now one can
investigate the role of different numerical and physical effects on the
statistics used to characterize the mass distribution of the Universe.

[29] oai:arXiv.org:astro-ph/0308027 [pdf] - 58306

IRAS PSCz v.s. 1.2-Jy Velocity and Density Fields: A Spherical Harmonics
Comparison

Comments: 12 pages, 15 figures, submitted to ApJ

Submitted: **2003-08-02**

We have used the two IRAS redshift surveys, 1.2-Jy (Fisher et al 1995) and
PSCz (Saunders et al 2000), to model the linear velocity fields within a
redshift of 8,000 km/s and have compared them in redshift space. The two
velocity fields only differ significantly in their monopole components. The
monopole discrepancy cannot be solely ascribed to shot-noise errors and
incomplete sky coverage. The mismatch seems to arise from incompleteness of the
PSCz catalog at fluxes $\le$ 1.2 Jy. The 1.2-Jy and PSCz higher order velocity
multipoles, particularly the dipole and quadrupole components, appear to be
consistent, suggesting that the dipole residuals found by Davis, Nusser and
Willick (1996) when comparing 1.2-Jy and MarkIII velocity fields probably
originates from the MarkIII velocities calibration procedure rather than from
uncertainties in the model velocity field. Our results illustrate the
efficiency of the spherical harmonics decomposition techniques in detecting
possible differences between real and model velocity fields. Similar analyses
shall prove to be very useful in the future to test the reliability of next
generation model velocity fields derived from new redshift catalogs like 2dFGRS
(Colless et al 2001), SDSS (York et al 2000), 6dF and 2MRS.

[30] oai:arXiv.org:astro-ph/0101349 [pdf] - 40487

Self-Interacting Scalar Dark Matter and Higgs Decay

Comments: Talk given by R. Rosenfeld at Cosmology 2000, Lisbon, July 12-15th.
To be published in the electronic Proceedings. One reference added

Submitted: **2001-01-19**, last modified: **2001-01-26**

Self-interacting dark matter has been suggested in order to overcome the
difficulties of the Cold Dark Matter model on galactic scales. We argue that a
scalar gauge singlet coupled to the Higgs boson, leading to an invisibly
decaying Higgs, is an interesting candidate for this self-interacting dark
matter particle.

[31] oai:arXiv.org:astro-ph/0011474 [pdf] - 1591854

PSCz vs. 1.2 Jy velocity fields: a Spherical Harmonics Comparison

Comments: 6 pages, 3 figures, uses World Scientific style. To appear in the
proceedings of X ENAA, 27-28 Jul. 2000, Lisbon, Portugal

Submitted: **2000-11-25**

We perform a detailed comparison of the IRAS PSCz and 1.2-Jy spherical
harmonic coefficients of the predicted velocity fields in redshift space. The
monopole terms predicted from the two surveys show some differences. Faint
galaxies are responsible for this mismatch that disappears when extracting a
PSCz subsample of galaxies with fluxes larger than 1.2 Jy. The analysis of PSCz
dipole components confirms the same inconsistencies found by Davis, Nusser and
Willick (1996) between the IRAS 1.2-Jy gravity field and MARKIII peculiar
velocities.
Shot-noise, which is greatly reduced in our PSCz gravity field, cannot be
responsible for the observed mismatch.

[32] oai:arXiv.org:astro-ph/0011171 [pdf] - 1507841

Large Scale Structure and Cosmic Rays revisited

Comments: 6 pages, 2 figures, uses World Scientific style. To appear in the
proceedings of "Third International Workshop on New Worlds in Astro-Particle
Physics", 1-3 Set. 2000, Faro, Portugal and X ENAA, 27-28 Jul. 2000, Lisbon,
Portugal

Submitted: **2000-11-08**

We investigate the possibility that ultra high energy cosmic rays (E > 10^19
eV) are related to the distribution of matter on large scales. The large scale
structure (LSS) data stems from the recent IRAS PSCz redshift survey. We
present preliminary predictions drawn from an anisotropic distribution of
sources which follows the galaxy distribution.

[33] oai:arXiv.org:astro-ph/0003350 [pdf] - 35224

Self-interacting Dark Matter and Invisibly Decaying Higgs

Comments: 4 pages, Revtex

Submitted: **2000-03-23**

Self-interacting dark matter has been suggested in order to overcome the
difficulties of the Cold Dark Matter model on galactic scales. We argue that a
scalar gauge singlet coupled to the Higgs boson, which could lead to an
invisibly decaying Higgs, is an interesting candidate for this self-interacting
dark matter particle. We also present estimates on the abundance of these
particles today as well as consequences to non-Newtonian forces.

[34] oai:arXiv.org:astro-ph/9909190 [pdf] - 108226

Density and Velocity Fields from the PSCz Survey

Comments: 8 pages, 8 figures. To appear in 'Towards an Understanding of Cosmic
Flows', Victoria, July 1999, eds Courteau,S., Strauss,M., Willick,J. PASP

Submitted: **1999-09-10**

We present the results for the predicted density and peculiar velocity fields
and the dipole from the PSCz survey of 15,000 IRAS galaxies over 84% of the
sky. We find a significant component to the dipole arising between 6000 and
15,000 km/s, but no significant component from greater distances. The
misalignment with the CMB is 20 degrees. The most remarkable feature of the
PSCz model velocity field is a coherent large-scale flow along the baseline
connecting Perseus-Pisces, the Local Supercluster, Great Attractor and the
Shapley Concentration. We have measured the parameter beta using the amplitude
of the dipole, bulk flow and point by point comparisons between the individual
velocities of galaxies in the MarkIII and SFI datasets, and the large-scale
clustering distortion in redshift space.All our results are consistent with
beta = 0.6 +- 0.1.

[35] oai:arXiv.org:astro-ph/9908358 [pdf] - 108030

PSCz-1.2 Jy Comparison: A Spherical Harmonics Approach

Comments: 4 pages Latex, 2 figures. Cosmic Flows Workshop, Victoria, Canada,
July 1999, ed. S. Courteau, M. Strauss & J. Willick, ASP series

Submitted: **1999-08-31**

We perform a detailed comparison of the IRAS PSCz and 1.2-Jy spherical
harmonic coefficients of the density and velocity fields in redshift space. The
monopole terms predicted from the two surveys show some differences. The
mismatch between the velocity monopoles arises from faint sources and
disappears when extracting a PSCz subsample of galaxies with fluxes larger than
1.2 Jy. The analysis of PSCz dipole moments confirms the same inconsistencies
found by Davis, Nusser and Willick (1996) between the IRAS 1.2-Jy gravity field
and MARK III peculiar velocities. We conclude that shot-noise, which is greatly
reduced in our PSCz gravity field, cannot be responsible for the observed
mismatch.

[36] oai:arXiv.org:astro-ph/9901366 [pdf] - 104960

A Non Parametric Model for the Cosmic Velocity Field

Comments: 29 pages, Latex, 25 figures (23 ps-files + 2 gif file). Accepted by
MNRAS. The correct text now replaces an old draft version, erroneously placed
on the net

Submitted: **1999-01-26**, last modified: **1999-02-03**

We present a self consistent nonparametric model of the local cosmic velocity
field based on the density distribution in the PSCz redshift survey of IRAS
galaxies. The error analysis, carried out on mock PSCz catalogues constructed
from N-body simulations, reveals uncertainties of ~70 km/sec. The denser
sampling provided by the PSCz survey compared to previous IRAS galaxy surveys
allows us to reconstruct the velocity field out to larger distances. The most
striking feature of the model velocity field is a coherent large-scale
streaming motion along a baseline connecting Perseus-Pisces, the Local
Supercluster, the Great Attractor, and the Shapley Concentration. We find no
evidence for back-infall onto the Great Attractor. Instead, material behind and
around the Great Attractor in inferred to be streaming towards the Shapley
Concentration The PSCz model velocities compare well with those predicted from
the 1.2Jy redshift survey of IRAS galaxies and with those predicted from the
distribution of Abell/ACO clusters, out to 140 Mpc/h. Comparison of the
real-space density fields (or, alternatively, the peculiar velocity fields)
inferred from the PSC$z$ and cluster catalogues gives a relative (linear) bias
parameter between clusters and IRAS galaxies of 4.4 +/- 0.6. Finally, we
compare the cumulative bulk flows predicted from the PSCz gravity field with
those measured from the Mark III and SFI catalogues of peculiar velocities. A
conservative estimate of beta=Omega**0.6/b, where b is the bias parameter for
IRAS galaxies, gives beta =0.6 + 0.22 -0.15, in agreement with other recent
determinations.

[37] oai:arXiv.org:astro-ph/9901087 [pdf] - 1943617

Likelihood Analysis of the Local Group Acceleration

Comments: replace earlier version, references added, errorbars corrected

Submitted: **1999-01-08**, last modified: **1999-01-12**

We compute the acceleration on the Local Group using 11206 IRAS galaxies from
the recently completed all-sky PSCz redshift survey. Measuring the acceleration
vector in redshift space generates systematic uncertainties due to the redshift
space distortions in the density field. We therefore assign galaxies to their
real space positions by adopting a non-parametric model for the velocity field
that solely relies on the linear gravitational instability and linear biasing
hypotheses. Remaining systematic contributions to the measured acceleration
vector are corrected for by using PSCz mock catalogues from N-body experiments.
The resulting acceleration vector points approx. 15 degrees away from the CMB
dipole apex, with a remarkable alignment between small and large scale
contributions. A considerable fraction of the measured acceleration is
generated within 40 h-1 Mpc with a non-negligible contribution from scales
between 90 and 140 h-1 Mpc after which the acceleration amplitude seems to have
converged. The local group acceleration from PSCz appears to be consistent with
the one determined from the IRAS 1.2 Jy galaxy catalogue once the different
contributions from shot noise have been taken into account. The results are
consistent with the gravitational instability hypothesis and do not indicate
any strong deviations from the linear biasing relation on large scales. A
maximum-likelihood analysis of the comulative PSCz dipole is performed within a
radius of 150 h-1 Mpc in which we account for nonlinear effects, shot noise and
finite sample size. We obtain \beta = 0.70 (+ 0.35)(-0.2) at 1 \sigma
confidence level.