Search

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., Comment: 19 pages, 13 figures, published in ApJ

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., Comment: Fixed typo in Appx B equation and updated urls

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., Comment: Published in ApJ. 14 pages, 11 figures

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%., Comment: 14 pages, 13 figures, minor changes to match published version

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., Comment: 8 pages, 7 figures and 3 tables

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., Comment: v2: 13 pages, 8 figures, accepted version, to be published in JGR-Planets; v1: 11 pages, 6 figures, submitted to JGR-Planets

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., Comment: 12 pages, 13 figures

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., Comment: 19 pages, 14 figures, updated version accepted by ApJ

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., Comment: 19 pages. 5 colour figures. Submitted to Icarus

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., Comment: Accepted for publication in MNRAS 15 pages, 8 figures Substantial revisions to original following referee's report

The abundance of galaxy clusters as a function of mass is determined using the 2dFGRS Percolation-Inferred Galaxy Group (2PIGG) catalogue. This is used to estimate the amplitude of the matter fluctuation spectrum, parametrised by the linear theory rms density fluctuations in spheres of 8Mpc/h, sigma_8. The best-fitting value for this parameter is highly correlated with the mean matter density in the Universe, Omega_m, and is found to satisfy sigma_8=0.25Omega_m^{-0.92-4.5(Omega_m-0.22)^2}+/-10%(statistical)+/-20% (systematic) for 0.18

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

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

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., Comment: 15 pages, 14 figures (1 colour), submitted to MNRAS

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., Comment: 16 pages. Submitted to MNRAS

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., Comment: 15 pages, 8 figures (1 in colour), MNRAS, 348, 866 (as before)

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)., Comment: Same as before

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

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

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

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

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., Comment: 19 pages, LaTeX, MNRAS format, 14 inlined and 2 external Postscript figures. Submitted to MNRAS

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., Comment: To be published in MNRAS 15 pages, 9 postscript figures

It has been known since the Lunar Prospector mission that the poles of the Moon evidently harbor enhanced concentrations of hydrogen [1,2]. The physical and chemical form of the hydrogen has been much debated. Using imagery from Clementine it was possible to roughly estimate permanently-shadowed regions (PSRs), and to perform image reconstructions of the Lunar Prospector epithermal neutron flux maps [3,4]. The hydrogen concentrations resulting from these reconstructions were consistent with a few weight percent water ice in selected locations. With the LCROSS impact, we now know that hydrogen in the form of ice does exist in lunar polar cold traps [5]. Armed with this information, and new data from LRO/Diviner, we can examine whether the pre-sent-day distribution of hydrogen in the form of water ice is consistent with a past large impact that delivered a large mass of volatiles to the lunar surface. These volatiles, mixed with solid impact ejecta, would then be lost from locations having high mean temperatures but would otherwise remain trapped in locations with sufficiently low mean annual temperatures [6]. The time scales for loss would depend on the location-dependent temperatures as well as impact history.