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.

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. 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. We find no evident departure from the trend at the extremes of the parameters explored. The scaling law can readily be incorporated into models of planet formation.

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.

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.

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.

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.

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.

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.

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.

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.

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

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

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.