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 establish the feasibility of measuring the neutron lifetime via an alternative, space-based class of methods, which use neutrons generated by galactic cosmic ray spallation of planets' surfaces and atmospheres. Free neutrons decay via the weak interaction with a mean lifetime of around 880 s. This lifetime constrains the unitarity of the CKM matrix and is a key parameter for studies of Big-Bang nucleosynthesis. However, current laboratory measurements, using two independent approaches, differ by over 4$\sigma$. Using data acquired in 2007 and 2008 during flybys of Venus and Mercury by NASA's MESSENGER spacecraft, which was not designed to make this measurement, we estimate the neutron lifetime to be $780\pm60_\textrm{stat}\pm70_\textrm{syst}$ s, thereby demonstrating the viability of this new approach.

We assess a claim that observed galaxy clusters with mass ~10^14 Msun are more centrally concentrated than predicted in LCDM. We generate mock strong gravitational lensing observations, taking the lenses from a cosmological hydrodynamical simulation, and analyse them in the same way as the real Universe. The observed and simulated lensing arcs are consistent with one another, with three main effects responsible for the previously claimed inconsistency. First, galaxy clusters containing baryonic matter have higher central densities than their counterparts simulated with only dark matter. Second, a sample of clusters selected because of the presence of pronounced gravitational lensing arcs preferentially finds centrally concentrated clusters with large Einstein radii. Third, lensed arcs are usually straighter than critical curves, and the chosen image analysis method (fitting circles through the arcs) overestimates the Einstein radii. After accounting for these three effects, LCDM predicts that galaxy clusters should produce giant lensing arcs that match those in the observed Universe.

Many distant objects can only be detected, or become more scientifically valuable, if they have been highly magnified by strong gravitational lensing. We use EAGLE and BAHAMAS, two recent cosmological hydrodynamical simulations, to predict the probability distribution for both the lens mass and lens redshift when point sources are highly magnified by gravitational lensing. For sources at a redshift of two, we find the distribution of lens redshifts to be broad, peaking at z=0.6. The contribution of different lens masses is also fairly broad, with most high-magnification lensing due to lenses with halo masses between 10^12 and 10^14 solar masses. Lower mass haloes are inefficient lenses, while more massive haloes are rare. We find that a simple model in which all haloes have singular isothermal sphere density profiles can approximately reproduce the simulation predictions, although such a model over-predicts the importance of haloes with mass <10^12 solar masses for lensing. We also calculate the probability that point sources at different redshifts are strongly lensed. At low redshift, high magnifications are extremely unlikely. Each z=0.5 source produces, on average, 5x10^-7 images with magnification greater than ten; for z =2 this increases to about 2x10^-5. Our results imply that searches for strongly lensed optical transients, including the optical counterparts to strongly lensed gravitational waves, can be optimized by monitoring massive galaxies, groups and clusters rather than concentrating on an individual population of lenses.

We present BAHAMAS-SIDM, the first large-volume, (400/h Mpc)^3, cosmological simulations including both self-interacting dark matter (SIDM) and baryonic physics. These simulations are important for two primary reasons: 1) they include the effects of baryons on the dark matter distribution 2) the baryon particles can be used to make mock observables that can be compared directly with observations. As is well known, SIDM haloes are systematically less dense in their centres, and rounder, than CDM haloes. Here we find that that these changes are not reflected in the distribution of gas or stars within galaxy clusters, or in their X-ray luminosities. However, gravitational lensing observables can discriminate between DM models, and we present a menu of tests that future surveys could use to measure the SIDM interaction strength. We ray-trace our simulated galaxy clusters to produce strong lensing maps. Including baryons boosts the lensing strength of clusters that produce no critical curves in SIDM-only simulations. Comparing the Einstein radii of our simulated clusters with those observed in the CLASH survey, we find that at velocities around 1000 km/s an SIDM cross-section of sigma/m > 1 cm^2/g is likely incompatible with observed cluster lensing.

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.

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 present the first simulated galaxy clusters (M_200 > 10^14 Msun) with both self-interacting dark matter (SIDM) and baryonic physics. They exhibit a greater diversity in both dark matter and stellar density profiles than their counterparts in simulations with collisionless dark matter (CDM), which is generated by the complex interplay between dark matter self-interactions and baryonic physics. Despite variations in formation history, we demonstrate that analytical Jeans modelling predicts the SIDM density profiles remarkably well, and the diverse properties of the haloes can be understood in terms of their different final baryon distributions.

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.

The observed stellar kinematics of dispersion-supported galaxies are often used to measure dynamical masses. Recently, several analytical relationships between the stellar line-of-sight velocity dispersion, the projected (2D) or deprojected (3D) half-light radius, and the total mass enclosed within the half-light radius, relying on the spherical Jeans equation, have been proposed. Here, we make use of the APOSTLE cosmological hydrodynamical simulations of the Local Group to test the validity and accuracy of such mass estimators for both dispersion and rotation-supported galaxies, for field and satellite galaxies, and for galaxies of varying masses, shapes, and velocity dispersion anisotropies. We find that the mass estimators of Walker et al. and Wolf et al. are able to recover the masses of dispersion-dominated systems with little systematic bias, but with a 1-sigma scatter of 25 and 23 percent, respectively. The error on the estimated mass is dominated by the impact of the 3D shape of the stellar mass distribution, which is difficult to constrain observationally. This intrinsic scatter becomes the dominant source of uncertainty in the masses estimated for galaxies like the dwarf spheroidal (dSph) satellites of the Milky Way, where the observational errors in their sizes and velocity dispersions are small. Such scatter may also affect the inner density slopes of dSphs derived from multiple stellar populations, relaxing the significance with which Navarro-Frenk-White profiles may be excluded, depending on the degree to which the relevant properties of the different stellar populations are correlated. Finally, we derive a new optimal mass estimator that removes the residual biases and achieves a statistically significant reduction in the scatter to 20 percent overall for dispersion-dominated galaxies, allowing more precise and accurate mass estimates.

We investigate how self-interacting dark matter (SIDM) with anisotropic scattering affects the evolution of isolated dark matter haloes as well as systems with two colliding haloes. For isolated haloes, we find that the evolution can be adequately captured by treating the scattering as isotropic, as long as the isotropic cross-section is appropriately matched to the underlying anisotropic model. We find that this matching should not be done using the momentum transfer cross-section, as has been done previously. Matching should instead be performed via a modified momentum transfer cross-section that takes into account that dark matter particles can be relabelled after they scatter, without altering the dynamics. However, using cross-sections that are matched to give the same behaviour in isolated haloes, we find that treating dark matter scattering as isotropic under-predicts the effects of anisotropic dark matter scattering when haloes collide. In particular, the DM-galaxy offset induced by SIDM in colliding galaxy clusters is larger when we simulate the underlying particle model, than if we use a matched isotropic model. On the other hand, well motivated particle models with anisotropic scattering typically have cross-sections with a strong velocity dependence, and we discover a previously unrecognised effect that suppresses DM-galaxy offsets in colliding clusters making it hard for these systems to provide competitive constraints on such particle models.

We perform numerical simulations of the merging galaxy cluster 1E 0657-56 (the Bullet Cluster), including the effects of elastic dark matter scattering. In a similar manner to the stripping of gas by ram pressure, dark matter self-interactions would transfer momentum between the two galaxy cluster dark matter haloes, causing them to lag behind the collisionless galaxies. The absence of an observed separation between the dark matter and stellar components in the Bullet Cluster has been used to place upper limits on the cross-section for dark matter scattering. We emphasise the importance of analysing simulations in an observationally-motivated manner, finding that the way in which the positions of the various components are measured can have a larger impact on derived constraints on dark matter's self-interaction cross-section than reasonable changes to the initial conditions for the merger. In particular, we find that the methods used in previous studies to place some of the tightest constraints on this cross-section do not reflect what is done observationally, and overstate the Bullet Cluster's ability to constrain the particle properties of dark matter. We introduce the first simulations of the Bullet Cluster including both self-interacting dark matter and gas. We find that as the gas is stripped it introduces radially-dependent asymmetries into the stellar and dark matter distributions. As the techniques used to determine the positions of the dark matter and galaxies are sensitive to different radial scales, these asymmetries can lead to erroneously measured offsets between dark matter and galaxies even when they are spatially coincident.

We estimate the rate of dark matter scattering in collapsed structures throughout the history of the Universe. If the scattering cross-section is velocity-independent, then the canonical picture is correct that scatterings occur mainly at late times. The scattering rate peaks slightly at redshift z~6, and remains significant today. Half the scatterings occur after z~1, in structures more massive than 10^12 M_sun. Within a factor of two, these numbers are robust to changes in the assumed astrophysics, and the scatterings would be captured in cosmological simulations. However, for particle physics models with a velocity-dependent cross-section (as for Yukawa potential interactions via a massive mediator), the scattering rate peaks before z~20, in objects with mass less than 10^4 M_sun. These precise values are sensitive to the redshift-dependent mass-concentration relation and the small-scale cutoff in the matter power spectrum. In extreme cases, the qualitative effect of early interactions may be reminiscent of warm dark matter and strongly affect the subsequent growth of structure. However, these scatterings are being missed in existing cosmological simulations with limited mass resolution.

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 use the "Evolution and Assembly of GaLaxies and their Environments" ( EAGLE ) suite of hydrodynamical cosmological simulations to measure offsets between the centres of stellar and dark matter components of galaxies. We find that the vast majority (>95%) of the simulated galaxies display an offset smaller than the gravitational softening length of the simulations (Plummer-equivalent $\epsilon = 700$ pc), both for field galaxies and satellites in clusters and groups. We also find no systematic trailing or leading of the dark matter along a galaxy's direction of motion. The offsets are consistent with being randomly drawn from a Maxwellian distribution with $\sigma \leq 196$ pc. Since astrophysical effects produce no feasible analogues for the $1.62^{+0.47}_{-0.49}$ kpc offset recently observed in Abell 3827, the observational result is in tension with the collisionless cold dark matter model assumed in our simulations.

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-likelihood weak lensing analysis of the mass distribution in optically selected spectroscopic Galaxy Groups (G$^3$Cv5) in the Galaxy And Mass Assembly (GAMA) survey, using background Sloan Digital Sky Survey (SDSS) photometric galaxies. The scaling of halo mass, $M_h$, with various group observables is investigated. Our main results are: 1) the measured relations of halo mass with group luminosity, virial volume and central galaxy stellar mass, $M_\star$, agree very well with predictions from mock group catalogues constructed from a GALFORM semi-analytical galaxy formation model implemented in the Millennium $\Lambda$CDM N-body simulation; 2) the measured relations of halo mass with velocity dispersion and projected half-abundance radius show weak tension with mock predictions, hinting at problems in the mock galaxy dynamics and their small scale distribution; 3) the median $M_h|M_\star$ measured from weak lensing depends more sensitively on the lognormal dispersion in $M_\star$ at fixed $M_h$ than it does on the median $M_\star|M_h$. Our measurements suggest an intrinsic dispersion of $\sigma_{\log(M_\star)}\sim 0.15$; 4) Comparing our mass estimates with those in the catalogue, we find that the G$^3$Cv5 mass can give biased results when used to select subsets of the group sample. Of the various new halo mass estimators that we calibrate using our weak lensing measurements, group luminosity is the best single-proxy estimator of group mass.

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.

Well-motivated elementary particle candidates for the dark matter, such as the sterile neutrino, behave as warm dark matter (WDM).For particle masses of order a keV, free streaming produces a cutoff in the linear fluctuation power spectrum at a scale corresponding to dwarf galaxies. We investigate the abundance and structure of WDM haloes and subhaloes on these scales using high resolution cosmological N-body simulations of galactic haloes of mass similar to the Milky Way's. On scales larger than the free-streaming cutoff, the initial conditions have the same power spectrum and phases as one of the cold dark matter (CDM) haloes previously simulated by Springel et al as part of the Virgo consortium Aquarius project. We have simulated four haloes with WDM particle masses in the range 1.4-2.3keV and, for one case, we have carried out further simulations at varying resolution. N-body simulations in which the power spectrum cutoff is resolved are known to undergo artificial fragmentation in filaments producing spurious clumps which, for small masses (<10^7Msun in our case) outnumber genuine haloes. We have developed a robust algorithm to identify these spurious objects and remove them from our halo catalogues. We find that the WDM subhalo mass function is suppressed by well over an order magnitude relative to the CDM case for masses <10^9Msun. Requiring that there should be at least as many subhaloes as there are observed satellites in the Milky Way leads to a conservative lower limit to the (thermal equivalent) WDM particle mass of ~1.5\rmn{keV}. WDM haloes and subhaloes have cuspy density distributions that are well described by NFW or Einasto profiles. Their central densities are lower for lower WDM particle masses and none of the models we have considered suffer from the "too big to fail" problem recently highlighted by Boylan-Kolchin et al.

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.

We investigate the luminosity functions (LFs) and projected number density profiles of galactic satellites around isolated primaries of different luminosities. We measure these quantities for model satellites placed into the Millennium and Millennium II dark matter simulations by the GALFORM semi-analytic galaxy formation model for different bins of primary galaxy magnitude and we investigate their dependence on satellite luminosity. We compare our model predictions to the data of Guo et al. from the Sloan Digital Sky Survey Data Release 8 (SDSS DR8). First, we use a mock light-cone catalogue to verify that the method we used to count satellites in the SDSS DR8 is unbiased. We find that the radial distributions of model satellites are similar to those around comparable primary galaxies in the SDSS DR8, with only slight differences at low luminosities and small projected radii. However, when splitting the satellites by colour, the model and SDSS satellite systems no longer resemble one another, with many red model satellites, in contrast to the dominant blue fraction at similar luminosity in SDSS. The few model blue satellites are also significantly less centrally concentrated in the halo of their stacked primary than their SDSS counterparts. The implications of this result for the GALFORM model are discussed.

Cold dark matter models predict the existence of a large number of substructures within dark matter halos. If the cold dark matter consists of weakly interacting massive particles, their annihilation within these substructures could lead to diffuse GeV emission that would dominate over the annihilation signal of the host halo. In this work we search for GeV emission from three nearby galaxy clusters: Coma, Virgo and Fornax. We first remove known extragalactic and galactic diffuse gamma-ray backgrounds and point sources from the Fermi 2-year catalog and find a significant residual diffuse emission in all three clusters. We then investigate whether this emission is due to (i) unresolved point sources; (ii) dark matter annihilation; or (iii) cosmic rays (CR). Using 45 months of Fermi-LAT data we detect several new point sources (not present in the Fermi 2-year point source catalogue) which contaminate the signal previously analyzed by Han et al.(arxiv:1201.1003). Including these and accounting for the effects of undetected point sources, we find no significant detection of extended emission from the three clusters studied. Instead, we determine upper limits on emission due to dark matter annihilation and cosmic rays. For Fornax and Virgo the limits on CR emission are consistent with theoretical models, but for Coma the upper limit is a factor of 2 below the theoretical expectation. Allowing for systematic uncertainties associated with the treatment of CR, the upper limits on the cross section for dark matter annihilation from our clusters are more stringent than those from analyses of dwarf galaxies in the Milky Way. We rule out the thermal cross section for supersymmetric dark matter particles for masses as large as 100 GeV (depending on the annihilation channel).

We study the spatial distribution of satellite galaxies around isolated primaries using the Sloan Digital Sky Survey (SDSS) spectroscopic and photometric galaxy catalogues. We select isolated primaries from the spectroscopic sample and search for potential satellites in the much deeper photometric sample. For specific luminosity primaries we obtain robust statistical results by stacking as many as ~50, 000 galaxy systems. We find no evidence for any anistropy in the satellite galaxy distribution relative to the major axes of the primaries. We derive accurate projected number density profiles of satellites down to 4 magnitudes fainter than their primaries. We find the normalized satellite profiles generally have a universal form and can be well fitted by projected NFW profiles. The NFW concentration parameter increases with decreasing satellite luminosity while being independent of the luminosity of the primary except for very bright primaries. The profiles of the faintest satellites show deviations from the NFW form with an excess at small galactocentric projected distances. In addition, we quantify how the radial distribution of satellites depends on the colour of the satellites and on the colour and concentration of their primaries.

We report evidence for extended gamma-ray emission from the Virgo, Fornax and Coma clusters based on a maximum-likelihood analysis of the 3-year Fermi-LAT data. For all three clusters, excess emission is observed within three degrees of the center, peaking at the GeV scale. This emission cannot be accounted for by known Fermi sources or by the galactic and extragalactic backgrounds. If interpreted as annihilation emission from supersymmetric dark matter (DM) particles, the data prefer models with a particle mass in the range 20-60 GeV annihilating into the b-bbar channel, or 2-10 GeV and >1 TeV annihilating into mu-mu final states. Our results are consistent with those obtained by Hooper and Linden from a recent analysis of Fermi-LAT data in the region of the Galactic Centre. An extended DM annihilation profile dominated by emission from substructures is preferred over a simple point source model. The significance of DM detection is 4.4 sigma in Virgo and lower in the other two clusters. We also consider the possibility that the excess emission arises from cosmic ray (CR) induced gamma-rays, and infer a CR level within a factor of three of that expected from analytical models. However, the significance of a CR component is lower than the significance of a DM component, and there is no need for such a CR component in the presence of a DM component in the preferred DM mass range. We also set flux and cross-section upper limits for DM annihilation into the b-bbar and mu-mu channels in all three clusters.

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.

High resolution N-body simulations of galactic cold dark matter haloes indicate that we should expect to find a few satellite galaxies around the Milky Way whose haloes have a maximum circular velocity in excess of 40 kms. Yet, with the exception of the Magellanic Clouds and the Sagittarius dwarf, which likely reside in subhaloes with significantly larger velocities than this, the bright satellites of the Milky Way all appear to reside in subhaloes with maximum circular velocities below 40 kms. As recently highlighted by Boylan-Kolchin et al., this discrepancy implies that the majority of the most massive subhaloes within a cold dark matter galactic halo are much too concentrated to be consistent with the kinematic data for the bright Milky Way satellites. Here we show that no such discrepancy exists if haloes are made of warm, rather than cold dark matter because these haloes are less concentrated on account of their typically later formation epochs. Warm dark matter is one of several possible explanations for the observed kinematics of the satellites.

We apply a basis function expansion method to create a time-evolving density/potential approximation of the late growth of simulated N-body dark matter haloes. We demonstrate how the potential of a halo from the Aquarius Project can be accurately represented by a small number of basis functions, and show that the halo expansion (HEX) method provides a way to replay simulations. We explore the level of accuracy of the technique as well as some of its limitations. We find that the number of terms included in the expansion must be large enough to resolve the large-scale distribution and shape of the halo but, beyond this, additional terms result in little further improvement. Particle and subhalo orbits can be integrated in this realistic, time-varying halo potential approximation, at much lower cost than the original simulation, with high fidelity for many individual orbits, and a good match to the distributions of orbital energy and angular momentum. Statistically, the evolution of structural subhalo properties, such as mass, half-mass radius and characteristic circular velocity, are very well reproduced in the halo expansion approximation over several gigayears. We demonstrate an application of the technique by following the evolution of an orbiting subhalo at much higher resolution than can be achieved in the original simulation. Our method represents a significant improvement over commonly used techniques based on static analytical descriptions of the halo potential.

We investigate the formation and evolution of satellite galaxies using smoothed particle hydrodynamics (SPH) simulations of a Milky Way(MW)-like system, focussing on the best resolved examples, analogous to the classical MW satellites. Comparing with a pure dark matter simulation, we find that the condensation of baryons has had a relatively minor effect on the structure of the satellites' dark matter halos. The stellar mass that forms in each satellite agrees relatively well over three levels of resolution (a factor of ~64 in particle mass) and scales with (sub)halo mass in a similar way in an independent semi-analytical model. Our model provides a relatively good match to the average luminosity function of the MW and M31. To establish whether the potential wells of our satellites are realistic, we measure their masses within observationally determined half-light radii, finding that the most massive examples have somewhat higher mass-to-light ratios than those derived for the MW dSphs from stellar kinematic data. A statistical test yields a ~9 percent probability that the simulated and observationally derived distributions of masses are consistent. Our results may suggest that either the MW halo is less massive than assumed in our simulations (~1.4e12 M_sun) or that there is substantial scatter in the satellite luminosity function or distribution of mass-to-light ratios at fixed host halo mass. Alternatively, feedback processes not properly captured by our simulations may have reduced the central densities of (sub)halos, or the subhalos may have initially formed with lower concentrations as would be the case, for example, if the dark matter were made of warm, rather than cold particles.

We study the luminosity function of satellite galaxies around isolated primaries using the Sloan Digital Sky Survey (SDSS) spectroscopic and photometric galaxy samples. We select isolated primaries from the spectroscopic sample and search for potential satellites in the much deeper photometric sample. For primaries of similar luminosity to the Milky Way and M31, we are able to stack as many as $\sim\negthinspace 20,000$ galaxy systems to obtain robust statistical results. We derive the satellite luminosity function extending almost 8 magnitudes fainter than the primary galaxy. We also determine how the satellite luminosity function varies with the luminosity, colour and concentration of the primary. We find that, in the mean, isolated primaries of comparable luminosity to the Milky Way and M31 contain about a factor of two fewer satellites brighter than $M_V=-14$ than the average of the Milky Way and M31.

We calculate the orbital angular momentum of dark matter subhaloes in the Aquarius simulations of cold dark matter galactic haloes. We calculate the orientation of their angular momentum relative to that of the spin vector of their host halo and find a variety of different configurations. All six Aquarius haloes contain statistically significant populations of subhalo orbits that are aligned with the main halo spin. All haloes posses a population of subhaloes that rotates in the same direction as the main halo and three of them possess, in addition, a population that rotates in the opposite direction. These configurations arise from the filamentary accretion of subhaloes. Quasi-planar distributions of coherently rotating satellites, such as those inferred in the Milky Way and other galaxies, arise naturally in simulations of a $\Lambda$CDM universe.

We describe and apply a simple prescription for defining connected structures in galaxy redshift surveys. The method is based upon two passes with a friends-of-friends groupfinder. The first pass uses a cylindrical linking volume to find galaxy groups and clusters, in order to suppress the line-of-sight smearing introduced by the large random velocities of galaxies within these deep potential wells. The second pass, performed with a spherical linking volume, identifies the connected components. This algorithm has been applied to the 2dFGRS, within which it picks out a total of 7,603 systems containing at least two galaxies and having a mean redshift less than 0.12. Connected systems with many members appear filamentary in nature, and the algorithm recovers two particularly large filaments within the 2dFGRS. For comparison, the algorithm is has also been applied to LambdaCDM mock galaxy surveys. While the model population of such systems is broadly similar to that in the 2dFGRS, it does not generally contain such extremely large structures.

We analyse the stellar and hot gas content of 18 nearby, low-mass galaxy clusters, detected in redshift space and selected to have a dynamical mass 3E14<M/Msun<6E14, as measured from the 2dF Galaxy Redshift Survey. We combine X-ray measurements from both Chandra and XMM with ground-based near-infrared observations from CTIO, AAT and CFHT to compare the mass in hot gas and stars to the dynamical mass and state of the clusters. Only 13 of the clusters are detected in X-ray emission, and for these systems we find that a range of 7-20 per cent of their baryonic mass, and <3 per cent of their dynamical mass, is detected in starlight, similar to what is observed in more massive clusters. In contrast, the five undetected clusters are underluminous in X-ray emission, by up to a factor 10, given their stellar mass. Although the velocity distribution of cluster members in these systems is indistinguishable from a Gaussian, all show subtle signs of being unrelaxed: either they lack a central, dominant galaxy, or the bright galaxy distribution is less concentrated and/or more elongated than the rest of the sample. Thus we conclude that low-mass clusters and groups selected from the velocity distribution of their galaxies exhibit a dichotomy in their hot gas properties. Either they are detected in X-ray, in which case they generally lie on the usual scaling relations, or they are completely undetected in X-ray emission. The non-detections may be partly related to the apparently young dynamical state of the clusters, but it remains a distinct possibility that some of these systems are exceptionally devoid of hot emitting gas as the result of its expulsion or rarefaction.

We investigate the magnitude and internal alignment of the angular momentum of cold dark matter haloes in simulations with and without baryons. We analyse the cumulative angular momentum profiles of hundreds of thousands of haloes in the Millennium simulation and in a smaller, but higher resolution, simulation, in total spanning 5 orders of magnitude in mass. For haloes of a given mass, the median specific angular momentum increases as j(<r) proportional to r. The direction of the vector varies considerably with radius: the median angle between the inner (< 0.25 Rvir) and total (< Rvir) angular momentum vectors is about 25degr. To investigate how baryons affect halo spin, we use another high resolution simulation, which includes gas cooling, star formation and feedback. This simulation produces a sample of galaxies with a realistic distribution of disc-to-total ratios. The formation of the galaxy spins up the dark matter within 0.1 Rvir such that the specific halo angular momentum increases by approx 50% in the median. The dark matter angular momentum becomes better aligned, but there remains a broad distribution of (mis-)alignments between the halo and the central galaxy, with a median angle between their angular momenta of ~ 30degr. Galaxies have a range of orientations relative to the shape of the halo: half of them have their minor axes misaligned by more than 45degr. Finally, we align a sample of haloes according to the orientation of their galaxies and stack the projected mass distributions. Although the individual haloes are aspherical, galaxy--halo misalignments produce a stacked mass distribution that is indistinguishable from circular. If the misalignment found in our simulations is realistic, it will be extremely difficult for weak lensing studies to measure the shape of dark haloes using this technique.

[Abridged] We present the first results of hydrodynamical simulations that follow the formation of galaxies to z=0 in spherical regions of radius ~20 Mpc/h drawn from the Millennium Simulation. The regions have overdensities that deviate by (-2, -1, 0, +1, +2)sigma from the cosmic mean, where sigma is the rms mass fluctuation on a scale of ~20Mpc/h at z=1.5. The simulations have mass resolution of up to 10^6 Msun/h, cover the entire range of large-scale environments and allow extrapolation of statistics to the entire 500 (Mpc/h)^3 Millennium volume. They include gas cooling, photoheating from an ionising background, SNe feedback and winds, but no AGN. We find that the specific SFR density at z <~ 10 varies systematically from region to region by up to an order of magnitude, but the global value, averaged over all volumes, reproduces observational data. Massive, compact galaxies, similar to those observed in the GOODS fields, form in the overdense regions as early as z=6, but do not appear in the underdense regions until z~3. These environmental variations are not caused by a dependence of the star formation properties on environment, but rather by a strong variation of the halo mass function from one environment to another, with more massive haloes forming preferentially in the denser regions. At all epochs, stars form most efficiently in haloes of circular velocity ~ 250 km/s. However, the star formation history exhibits a form of "downsizing" (even in the absence of AGN): the stars comprising massive galaxies at z=0 have mostly formed by z=1-2, whilst those comprising smaller galaxies typically form at later times. However, additional feedback is required to limit star formation in massive galaxies at late times.

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.

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.

(Abridged) We demonstrate that the tenet of hierarchical structure growth leads directly to a robust, falsifiable prediction for the correlation between stellar fraction (fstar) and total system mass (M500) of galaxy groups and clusters. This prediction is relatively insensitive to the details of baryonic physics or cosmological parameters. In particular, if the fstar-M500 relation is fixed and does not evolve with redshift, CDM models predict the logarithmic slope of this relation to be b>-0.3. This constraint can be weakened if the fstar-M500 relation evolves strongly, but this implies more stars must be formed in situ in groups at low redshift. Conservatively requiring that at least half the stars in groups were formed by z=1, the constraint from evolution models is b>-0.35. Since the most massive clusters (M500=1E15 Msun) are observed to have fstar=0.01, this means that groups with M500=5E13 Msun must have fstar<0.03. Recent observations by Gonzalez et al. (2007) indicate a much steeper relation, with fstar>0.04 in groups, leading to b=-0.64. If confirmed, this would rule out hierarchical structure formation models: today's clusters could not have been built from today's groups, or even from the higher-redshift progenitors of those groups. We perform a careful analysis of these and other data to identify the most important systematic uncertainties in their measurements. Although correlated uncertainties on stellar and total masses might explain the steep observed relation, the data are only consistent with theory if the observed group masses are systematically underestimated.

We investigate the spins and shapes of over a million dark matter haloes identified at z=0 in the Millennium simulation. Our sample spans halo masses ranging from dwarf galaxies to rich galaxy clusters. The very large dynamic range of this simulation enables the distribution of spins and shapes and their variation with halo mass and environment to be characterised with unprecedented precision. We compare results for haloes identified using different algorithms, including a novel one based on the branches of the halo merger trees. We investigate (and remove) biases in the estimate of angular momentum introduced by both the algorithms themselves and by numerical effects. We find that for this many objects, the traditional lognormal function is no longer an adequate description of the distribution of the dimensionless spin parameter lambda, and we provide a different function that gives a better fit. The variation of spin with halo mass is weak but detectable, although the trend depends strongly on the halo definition used. The haloes exhibit a range of shapes, with a preference for prolateness over oblateness. More massive haloes tend to be less spherical and more prolate. We find that the more spherical haloes have spin in the median, and those closest to spherical have a spin independent of mass. The most massive have a spin independent of shape. The majority of haloes have their angular momentum vector aligned with their minor axis and perpendicular to their major axis. We find that higher tend to be more clustered, with a stronger effect for more massive haloes. (abridged)

We investigate the baryon fraction in dark matter haloes formed in non-radiative gas-dynamical simulations of the LambdaCDM cosmogony. By combining a realisation of the Millennium Simulation (Springel et al.) with a simulation of a smaller volume focussing on dwarf haloes, our study spans five decades in halo mass, from 10^10 Msun/h to 10^15 Msun/h. We find that the baryon fraction within the halo virial radius is typically 90% of the cosmic mean, with an rms scatter of 6%, independently of redshift and of halo mass down to the smallest resolved haloes. Our results show that, contrary to the proposal of Mo et al. (2005), pre-virialisation gravitational heating is unable to prevent the collapse of gas within galactic and proto-galactic haloes, and confirm the need for non-gravitational feedback in order to reduce the efficiency of gas cooling and star formation in dwarf galaxy haloes. Simulations including a simple photoheating model (where a gas temperature floor of T_{floor} = 2x10^4 K is imposed from z=11) confirm earlier suggestions that photoheating can only prevent the collapse of baryons in systems with virial temperatures T_{200} < ~2.2 T_{floor} ~ 4.4x10^4 K (corresponding to a virial mass of M_{200} ~ 10^10 Msun/h and a circular velocity of V_{200} ~ 35 km/s). Photoheating may thus help regulate the formation of dwarf spheroidals and other galaxies at the extreme faint-end of the luminosity function, but it cannot, on its own, reconcile the abundance of sub-L* galaxies with the vast number of dwarf haloes expected in the LambdaCDM cosmogony. The lack of evolution or mass dependence seen in the baryon fraction augurs well for X-ray cluster studies that assume a universal and non-evolving baryon fraction to place constraints on cosmological parameters.

The abundance of galaxy clusters as a function of mass is determined using the 2dFGRS Percolation-Inferred Galaxy Group (2PIGG) catalogue. This is used to estimate the amplitude of the matter fluctuation spectrum, parametrised by the linear theory rms density fluctuations in spheres of 8Mpc/h, sigma_8. The best-fitting value for this parameter is highly correlated with the mean matter density in the Universe, Omega_m, and is found to satisfy sigma_8=0.25Omega_m^{-0.92-4.5(Omega_m-0.22)^2}+/-10%(statistical)+/-20% (systematic) for 0.18<Omega_m<0.50, assuming that Omega_m+Omega_Lambda=1. This gives sigma_8=0.89 when Omega_m=0.25. A ~20% correction has been applied to undo the systematic bias inherent in the measurement procedure. Mock catalogues, constructed from large cosmological N-body simulations, are used to help understand and model these systematic errors. The abundance of galaxy groups as a function of group bj band luminosity is also determined. This is used in conjunction with the halo mass function, determined from simulations, to infer the variation of halo mass-to-light ratio over four orders of magnitude in halo mass. The mass-to-light ratio shows a minimum value of 100hM_sol/L_sol in the bj band at a total group luminosity of L_bJ ~ 5.10^9Lsol/h^2. Together with the observed Tully-Fisher relation, this implies that the observed rotation speed of Tully-Fisher galaxies is within ~10% of the typical circular speed of haloes hosting brightest galaxies of the same luminosity.

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

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

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.

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.

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.

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

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.

We analyse the observed correlation between galaxy environment and H-alpha emission line strength, using volume-limited samples and group catalogues of 24968 galaxies drawn from the 2dF Galaxy Redshift Survey (Mb<-19.5) and the Sloan Digital Sky Survey (Mr<-20.6). We characterise the environment by 1) Sigma_5, the surface number density of galaxies determined by the projected distance to the 5th nearest neighbour; and 2) rho1.1 and rho5.5, three-dimensional density estimates obtained by convolving the galaxy distribution with Gaussian kernels of dispersion 1.1 Mpc and 5.5 Mpc, respectively. We find that star-forming and quiescent galaxies form two distinct populations, as characterised by their H-alpha equivalent width, EW(Ha). The relative numbers of star-forming and quiescent galaxies varies strongly and continuously with local density. However, the distribution of EW(Ha) amongst the star-forming population is independent of environment. The fraction of star-forming galaxies shows strong sensitivity to the density on large scales, rho5.5, which is likely independent of the trend with local density, rho1.1. We use two differently-selected group catalogues to demonstrate that the correlation with galaxy density is approximately independent of group velocity dispersion, for sigma=200-1000 km/s. Even in the lowest density environments, no more than ~70 per cent of galaxies show significant H-alpha emission. Based on these results, we conclude that the present-day correlation between star formation rate and environment is a result of short-timescale mechanisms that take place preferentially at high redshift, such as starbursts induced by galaxy-galaxy interactions.

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

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)

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

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

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

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

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

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.

A speedy pixon algorithm for image reconstruction is described. Two applications of the method to simulated astronomical data sets are also reported. In one case, galaxy clusters are extracted from multiwavelength microwave sky maps using the spectral dependence of the Sunyaev-Zel'dovich effect to distinguish them from the microwave background fluctuations and the instrumental noise. The second example involves the recovery of a sharply peaked emission profile, such as might be produced by a galaxy cluster observed in X-rays. These simulations show the ability of the technique both to detect sources in low signal-to-noise data and to deconvolve a telescope beam in order to recover the internal structure of a source.

Hydrodynamical simulations of galaxy formation in spatially flat Cold Dark Matter (CDM) cosmologies with and without a cosmological constant (Lambda) are described. A simple star formation algorithm is employed and radiative cooling is allowed only after redshift z=1 so that enough hot gas is available to form large, rapidly rotating stellar discs if angular momentum is approximately conserved during collapse. The specific angular momenta of the final galaxies are found to be sensitive to the assumed background cosmology. This dependence arises from the different angular momenta contained in the haloes at the epoch when the gas begins to collapse and the inhomogeneity of the subsequent halo evolution. In the Lambda-dominated cosmology, the ratio of stellar specific angular momentum to that of the dark matter halo (measured at the virial radius) has a median value of ~0.24 at z=0. The corresponding quantity for the Lambda=0 cosmology is over 3 times lower. It is concluded that the observed frequency and angular momenta of disc galaxies pose significant problems for spatially flat CDM models with Lambda=0 but may be consistent with a Lambda-dominated CDM universe.

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

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.

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.

We present results of N-body/gasdynamical simulations designed to investigate the evolution of X-ray clusters in a flat, low-density, cold dark matter (CDM) cosmogony. The density profile of the dark matter component can be fitted rather accurately by the simple formula originally proposed by Navarro, Frenk & White to describe the structure of clusters in a CDM universe with $\Omega=1$. In projection, the shape of the dark matter radial density profile and the corresponding line-of-sight velocity dispersion profile are in very good agreement with the observed profiles for galaxies in the CNOC sample of clusters. The gas in our simulated clusters is less centrally concentrated than the dark matter, and its radial density profile is well described by the familiar $\beta$-model. The total mass and velocity dispersion of our clusters can be accurately inferred (with $\sim 15%$ uncertainty) from their X-ray emission-weighted temperature. We generalize Kaiser's scaling relations for scale-free universes and show that the clusters in our simulations generally follow these relations. The agreement between the simulations and the analytical results provides a convincing demonstration of the soundness of our gasdynamical numerical techniques. The slope of the luminosity-temperature relation implied by the scaling relations, and obeyed by the simulations, is in disagreement with observations. This suggests that non-gravitational effects such as preheating or cooling must have played an important role in determining the properties of the observed X-ray emission from galaxy clusters.

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.

We use N-body simulations to examine the effects of mass outflows on the density profiles of cold dark matter (CDM) halos surrounding dwarf galaxies. In particular, we investigate the consequences of supernova-driven winds that expel a large fraction of the baryonic component from a dwarf galaxy disk after a vigorous episode of star formation. We show that this sudden loss of mass leads to the formation of a core in the dark matter density profile, although the original halo is modeled by a coreless (Hernquist) profile. The core radius thus created is a sensitive function of the mass and radius of the baryonic disk being blown up. The loss of a disk with mass and size consistent with primordial nucleosynthesis constraints and angular momentum considerations imprints a core radius which is only a small fraction of the original scale-length of the halo. These small perturbations are, however, enough to reconcile the rotation curves of dwarf irregulars with the density profiles of haloes formed in the standard CDM scenario.

The correlation function of galaxy clusters has often been used as a test of cosmological models. A number of assumptions are implicit in the comparison of theoretical expectations to data. Here we use an ensemble of ten large N-body simulations of the standard cold dark matter cosmology to investigate how cluster selection criteria and other uncertain factors influence the cluster correlation function. We consider the effects of varying the definition of a cluster, the mean number density (or equivalently the richness) in a catalogue, and the assumed normalisation of the model; we also examine the importance of redshift space distortions. We implement five different group-finding algorithms and construct cluster catalogues defined by mass, velocity dispersion or a measure of X-ray luminosity. We find that different cluster catalogues yield correlation functions which can differ by substantially more than the statistical errors in any one determination. For a given cluster selection criteria, the correlation length typically varies by $\sim 20\%$ in catalogues spanning the range of intercluster separations covered by the APM and Abell (richness class $\gsim 1$) catalogues. Distortions produced by peculiar velocities in redshift space enhance the correlation function at large separations and lead to a larger clustering length in redshift space than in real space. The sensitivity of the cluster correlation function to various uncertain model assumptions substantially weakens previous conclusions based on the comparison of model predictions with real data. Detailed modelling of cluster selection procedures including the effects of selecting from projected galaxy catalogues is required before the cluster correlation function can be regarded as a high precision constraint on cosmological models.

The population of rich galaxy clusters evolves much more rapidly in a universe with critical density than one with low density, thus offering the possibility of determining the cosmological density parameter, Omega_0. We quantify this evolution using the Press-Schechter formalism which we extend to flat models with a cosmological constant. Using new large N-body simulations, we verify that this formalism accurately predicts the abundance of rich clusters as a function of redshift in various cosmologies. We normalise the models by comparing them to the local abundance of clusters as a function of their X-ray temperature which we rederive from data compiled by Henry & Arnaud. This gives values of the rms density fluctuation in spheres of radius 8 Mpc/h of sigma_8 = (0.50+/- 0.04) Omega_0^{-0.47+0.10 Omega_0} if Lambda_0=0 and sigma_8 = (0.50 +/- 0.04) Omega_0^{-0.53+0.13 Omega_0} if Lambda_0=1-Omega_0. These values depend very weakly on the shape of the power spectrum. We then examine how the distributions of mass, X-ray temperature and Sunyaev-Zel'dovich decrement evolve as a function of Omega_0. We present the expected distributions at z=0.33 and z=0.5 and the predicted number counts of the largest clusters. We find that even at z=0.33, these distributions depend very strongly on Omega_0 and only weakly on Lambda_0. For example, at this redshift, we expect 20 times as many clusters per comoving volume with M>3.5 10^{14} Msol/h and 5 times as many clusters with kT>5 keV if Omega_0=0.3 than if Omega_0=1. The splitting in the integrated counts is enhanced by the larger volume element in low Omega_0 models. There is therefore a real prospect of estimating Omega_0 from forthcoming surveys of intermediate redshift clusters that will determine their masses, X-ray temperatures or SZ decrements.