We present integral field spectroscopy of galaxy cluster Abell 3827, using ALMA and VLT/MUSE. It reveals an unusual configuration of strong gravitational lensing in the cluster core, with at least seven lensed images of a single background spiral galaxy. Lens modelling based on HST imaging had suggested that the dark matter associated with one of the cluster's central galaxies may be offset. The new spectroscopic data enable better subtraction of foreground light, and better identification of multiple background images. The inferred distribution of dark matter is consistent with being centered on the galaxies, as expected by LCDM. Each galaxy's dark matter also appears to be symmetric. Whilst we do not find an offset between mass and light (suggestive of self-interacting dark matter) as previously reported, the numerical simulations that have been performed to calibrate Abell 3827 indicate that offsets and asymmetry are still worth looking for in collisions with particular geometries. Meanwhile, ALMA proves exceptionally useful for strong lens image identifications.

We compare predictions of large-scale cosmological hydrodynamical simulations for neutral hydrogen absorption signatures in the vicinity of 1e11 - 1e12.5 MSun haloes with observational measurements. Two different hydrodynamical techniques and a variety of prescriptions for gas removal in high density regions are examined. Star formation and wind feedback play only secondary roles in the HI absorption signatures outside the virial radius, but play important roles within. Accordingly, we identify three distinct gaseous regions around a halo: the virialized region, the mesogalactic medium outside the virial radius arising from the extended haloes of galaxies out to about two turnaround radii, and the intergalactic medium beyond. Predictions for the amount of absorption from the mesogalactic and intergalactic media are robust across different methodologies, and the predictions agree with the amount of absorption observed around star-forming galaxies and QSO host galaxies. Recovering the measured amount of absorption within the virialized region, however, requires either a higher dynamic range in the simulations, additional physics, or both.

Galaxy cluster Abell 3827 hosts the stellar remnants of four almost equally bright elliptical galaxies within a core of radius 10kpc. Such corrugation of the stellar distribution is very rare, and suggests recent formation by several simultaneous mergers. We map the distribution of associated dark matter, using new Hubble Space Telescope imaging and VLT/MUSE integral field spectroscopy of a gravitationally lensed system threaded through the cluster core. We find that each of the central galaxies retains a dark matter halo, but that (at least) one of these is spatially offset from its stars. The best-constrained offset is 1.62+/-0.48kpc, where the 68% confidence limit includes both statistical error and systematic biases in mass modelling. Such offsets are not seen in field galaxies, but are predicted during the long infall to a cluster, if dark matter self-interactions generate an extra drag force. With such a small physical separation, it is difficult to definitively rule out astrophysical effects operating exclusively in dense cluster core environments - but if interpreted solely as evidence for self-interacting dark matter, this offset implies a cross-section sigma/m=(1.7+/-0.7)x10^{-4}cm^2/g x (t/10^9yrs)^{-2}, where t is the infall duration.

Collisions between galaxy clusters provide a test of the non-gravitational forces acting on dark matter. Dark matter's lack of deceleration in the `bullet cluster collision' constrained its self-interaction cross-section \sigma_DM/m < 1.25cm2/g (68% confidence limit) for long-ranged forces. Using the Chandra and Hubble Space Telescopes we have now observed 72 collisions, including both `major' and `minor' mergers. Combining these measurements statistically, we detect the existence of dark mass at 7.6\sigma significance. The position of the dark mass has remained closely aligned within 5.8+/-8.2 kpc of associated stars: implying a self-interaction cross-section \sigma_DM/m < 0.47 cm2/g (95% CL) and disfavoring some proposed extensions to the standard model.

We perform cosmological simulations of the intergalactic medium (IGM) at redshift z ~ 3 using the numerical gravity-hydrodynamics codes GADGET-3 and Enzo for the purpose of modelling the gaseous environments of galaxies. We identify haloes in the simulations using three different algorithms. Different rank orderings of the haloes by mass result, introducing a limiting factor in identifying haloes with observed galaxies. We also compare the physical properties of the gas between the two codes, focussing primarily on the gas outside the virial radius, motivated by recent HI absorption measurements of the gas around z ~ 2 - 3 galaxies. The internal dispersion velocities of the gas in the haloes have converged for a box size of 30 comoving Mpc, but the centre-of-mass peculiar velocities of the haloes have not up to a box size of 60 comoving Mpc. The density and temperature of the gas within the instantaneous turn-around radii of the haloes are adequately captured for box sizes 30 Mpc on a side, but the results are highly sensitive to the treatment of unresolved, rapidly cooling gas, with the gas mass fraction within the virial radius severely depleted by star formation in the GADGET-3 simulations. Convergence of the gas peculiar velocity field on large scales requires a box size of at least 60 Mpc. Outside the turn-around radius, the physical state of the gas agrees to 30 percent or better both with box size and between simulation methods. We conclude that generic IGM simulations make accurate predictions for the intergalactic gas properties beyond the halo turn-around radii, but the gas properties on smaller scales are highly dependent on star formation and feedback implementations.

We develop a statistical method to measure the interaction cross-section of Dark Matter, exploiting the continuous minor merger events in which small substructures fall into galaxy clusters. We find that by taking the ratio of the distances between the galaxies and Dark Matter, and galaxies and gas in accreting sub-halos, we form a quantity that can be statistically averaged over a large sample of systems whilst removing any inherent line-of-sight projections. In order to interpret this ratio as a cross-section of Dark Matter we derive an analytical description of sub-halo infall which encompasses; the force of the main cluster potential, the drag on a gas sub-halo, a model for Dark Matter self-interactions and the resulting sub-halo drag, the force on the gas and galaxies due to the Dark Matter sub-halo potential, and finally the buoyancy on the gas and Dark Matter. We create mock observations from cosmological simulations of structure formation and find that collisionless Dark Matter becomes physically separated from X-ray gas by up to 20h^-1 kpc. Adding realistic levels of noise, we are able to predict achievable constraints from observational data. Current archival data should be able to detect a difference in the dynamical behaviour of Dark Matter and standard model particles at 6 sigma, and measure the total interaction cross-section sigma/m with 68% confidence limits of +/- 1cm2g^-1. We note that this method is not restricted by the limited number of major merging events and is easily extended to large samples of clusters from future surveys which could potentially push statistical errors to 0.1cm^2g^-1.

Direct evidence for the existence of dark matter and measurements of its interaction cross-section have been provided by the physical offset between dark matter and intra- cluster gas in merging systems like the Bullet Cluster. Although a smaller signal, this effect is more abundant in minor mergers where infalling substructure dark matter and gas are segregated. In such low-mass systems the gravitational lensing signal comes primarily from weak lensing. A fundamental step in determining such an offset in sub- structure is the ability to accurately measure the positions of dark matter sub-peaks. Using simulated Hubble Space Telescope observations, we make a first assessment of the precision and accuracy with which we can measure infalling groups using weak gravitational lensing. We demonstrate that using an existing and well-used mass re- construction algorithm can measure the positions of 1.5x1013 M substructures that have parent halos ten times more massive with a bias of less than 0.3 . In this regime, our analysis suggests the precision is sufficient to detect (at 3 {\sigma} statistical significance) the expected mean offset between dark matter and baryonic gas in infalling groups from a sample of 50 massive clusters.

The first major star-forming galaxies and Active Galactic Nuclei will produce Balmer and higher order extended haloes during the Epoch of Reionization through the scattering of Lyman resonance line photons off the surrounding neutral intergalactic gas. The optical depth dependence of the scattering rates will produce a signal sensitive to both the density and velocity fluctuations of the gas, offering the possibility of probing the ionization region and flow field surrounding young star-forming galaxies. The requirements for detecting the haloes in the infra-red using a space-based telescope are discussed, along with an assessment of the possibility of detecting the haloes using the Tunable Filter Imager on the James Webb Space Telescope.

We examine the impact of helium reionization on the structure of the intergalactic medium (IGM). We model the reionization using a radiative transfer (RT) code coupled to the combined gravity hydrodynamics code Enzo. Neutral hydrogen and helium are initially ionized by a starburst spectrum, which is allowed to gradually evolve into a power law spectrum over the redshift interval 3 < z < 4. The temperature-density relation of the gas is found to fan out and flatten following HeII reionization, with an inversion for overdensities above 5. Peculiar velocities of up to 10 km/s are induced by the increased pressure, with the gas density field distorted over large coherent regions by 10-20%, and the dark matter by levels of 1%. The photoionization-induced flows may thus distort the matter power spectrum at comoving wavenumbers k > 0.5 h/Mpc by a few percent by z = 2. Absorption spectra for HI and HeII are drawn from the simulations, and absorption lines are fit to the spectra. A median Doppler parameter of 35 km/s is obtained for the HI absorption systems at z = 3. Dividing into subsamples optically thick and optically thin at line centre reveals that the optically thick systems undergo only mild evolution while the optically thin systems evolve rapidly following HeII reionization. A comparison between HeII and HI absorption features shows a broad distribution in the HeII and HI column density ratio, peaking near the measured value and only slightly narrower than measured. A comparison with approximate simulation methods shows moderately good agreement in the absorption line properties, but not to the precision to which they may be measured.

We examine the temperature structure of the intergalactic medium IGM) surounding a hard radiation source, such as a Quasi-Stellar Object (QSO), as it responds to the onset of helium reionization by the source. We model the reionization using a radiative transfer (RT) code coupled to a particle-mesh (PM) N-body code. Neutral hydrogen and helium are initially ionized by a starburst spectrum, which is allowed to gradually evolve into a power law spectrum (fnu ~ nu^(-0.5)). Multiple simulations were performed with different times for the onset and dominance of the hard spectrum, with onset redshifts ranging from z = 3.5 to 5.5. The source is placed in a high-density region to mimic the expected local environment of a QSO. Simulations with the source placed in a low-density environment were also performed as control cases to explore the role of the environment on the properties of the surrounding IGM. We find in both cases that the IGM temperature within the HeIII region produced exceeds the IGM temperature before full helium reionization, resulting in a "thermal proximity effect", but that the temperature in the HeIII region increases systematically with distance from the source. With time the temperature relaxes with a reduced spread as a function of impact parameter along neighbouring lines of sight, although the trend continues to persist until z = 2. Such a trend could be detected using the widths of intervening metal absorption systems using high resolution, high signal-to-noise ratio spectra.

We examine the temperature structure of the IGM due to the passage of individual ionisation fronts using a radiative transfer (RT) code coupled to a particle-mesh (PM) N-body code. Multiple simulations were performed with different spectra of ionising radiation: a power law (goes as nu^{-0.5}), miniquasar, starburst, and a time-varying spectrum that evolves from a starburst spectrum to a power law. The RT is sufficiently resolved in time and space to correctly model both the ionisation state and the temperature across the ionisation front. We find the post-ionisation temperature of the reionised intergalactic medium (IGM) is sensitive to the spectrum of the source of ionising radiation, which may be used to place strong constraints on the nature of the sources of reionisation. Radiative transfer effects also produce large fluctuations in the HeII to HI number density ratio eta. The spread in values is smaller than measured, except for the time-varying spectrum. For this case, the spread evolves as the spectral nature of the ionising background changes. Large values for eta are found in partially ionised HeII as the power-law spectrum begins to dominate the starburst, suggesting that the large eta values measured may be indicating the onset of the HeII reionisation epoch.

We use numerical simulations with hydrodynamics to demonstrate that a class of cold fronts in galaxy clusters can be attributed to oscillations of the dark matter distribution. The oscillations are initiated by the off-axis passage of a low-mass substructure. From the simulations, we derive three observable morphological features indicative of oscillations: 1) The existence of compressed isophotes; 2) The regions of compression must be alternate (opposite and staggered) and lie on an axis passing through the center of the cluster; 3) The gradient of each compression region must pass through the center of the cluster. Four of six clusters reported in the literature to have cold fronts have morphologies consistent with the presence of oscillations. The clusters with oscillations are A496, A1795, A2142, and RX J1720.1+2638. Galaxy clusters A2256 and A3667 are not consistent so the cold fronts are interpreted as group remnants. The oscillations may be able to provide sufficient energy to solve the cooling-flow problem and, importantly, provide it over an extended duration.

An XMM-Newton imaging spectroscopy analysis of the galaxy cluster A1644 is presented. A1644 is a complex merging system consisting of a main and a sub cluster. A trail of cool, metal-rich gas has been discovered close to the sub cluster. The combination of results from X-ray, optical, and radio data, and a comparison to a hydrodynamical simulation suggest that the sub cluster has passed by the main cluster off-axis and a fraction of its gas has been stripped off during this process. Furthermore, for this merging system, simple effects are illustrated which can affect the use of clusters as cosmological probes. Specifically, double clusters may affect estimates of the cluster number density when treated as a single system. Mergers, as well as cool cores, can alter the X-ray luminosity and temperature measured for clusters, causing these values to differ from those expected in equilibrium.

We have analyzed the Chandra, BeppoSax, and ROSAT observations of Abell 754 and report evidence of a soft, diffuse X-ray component. The emission is peaked in the cluster center and is detected out to 8' from the X-ray center. Fitting a thermal model to the combined BeppoSax and PSPC spectra show excess emission below 1 keV in the PSPC and above 100 keV in the BeppoSax PDS. The source 26W20 is in the field of view of the PDS. The addition of a powerlaw with the spectral parameters measured by Silverman et al. (1998) for 26W20 successfully models the hard component in the PDS. The remaining excess soft emission can be modeled by either a low temperature, 0.75 - 1.03 keV component, or by a powerlaw with a steep spectral index, 2.3. Addition of a second thermal component model provides a much better fit to the data than does the addition of a non-thermal component. The Chandra temperature map does not show any region cooler than 6.9 keV within the region where the cool component was detected. Simulations of the emission from embedded groups were performed and compared with the Chandra temperature map which show groups are a plausible source of ~1 keV emission. The cool component is centrally peaked in the cluster and the gas density and temperature are relatively high arguing against the WHIM as the source of the X-ray emission. X-ray emission from elliptical galaxies is not high enough to provide the total cool component luminosity, 7.0x10^43 ergs s^-1. The peak of the cool component is located between the low frequency radio halos arguing against a non-thermal interpretation for the emission. We conclude that emission from embedded groups is the most likely origin of the cool component in Abell 754.

Analysis of a 30,000 s X-ray observation of the Abell 3266 galaxy cluster with the ACIS on board the Chandra Observatory has produced several new insights into the cluster merger. The intracluster medium has a non-monotonically decreasing radial abundance profile. We argue that the most plausible origin for the abundance enhancement is unmixed, high abundance subcluster gas from the merger. The enrichment consists of two stages: off-center deposition of a higher abundance material during a subcluster merger followed by a strong, localized intracluster wind that acts to drive out the light elements, producing the observed abundance enhancement. The wind is needed to account for both an increase in the heavy element abundance and the lack of an enhancement in the gas density. Dynamical evidence for the wind includes: (1) a large scale, low surface brightness feature perpendicular to the merger axis that appears to be an asymmetric pattern of gas flow to the northwest, away from the center of the main cluster, (2) compressed gas in the opposite direction (toward the cluster center), and (3), the hottest regions visible in the temperature map coincide with the proposed merger geometry and the resultant gas flow. The Chandra data for the central region of the main cluster shows a slightly cooler, filamentary region that is centered on the central cD galaxy and is aligned with the merger axis directly linking the dynamical state of the cD to the merger. Overall, the high spectral/spatial resolution Chandra observations support our earlier hypothesis (Henriksen, Donnelly, & Davis 1999) that we are viewing a minor merger in the plane of the sky.

We examine the necessity of resolving the Jeans mass in hydrodynamical simulations of the hierarchical formation of cosmological structures. We consider a standard two-component fluid and use smoothed particle hydrodynamics to model the baryonic component. It is found that resolution of the Jeans mass is not necessary to extract bulk properties of bound objects that are independent of resolution, provided that the objects have undergone a number of merger events. The degree of merging must be sufficient to populate the objects with substructure of similar total mass. After this criterion is met, neither the density profiles of structures nor their density-temperature distributions depend appreciably on the resolution of the simulation. Resolution of the Jeans mass by the imposition of a minimum temperature is not found to adversely affect the state of the matter in the final clusters. The baryon concentration profiles still indicate a deficit of gas with respect to the dark matter in the interiors of structures.

The baryon fraction of galaxy clusters in numerical simulations is found to be dependant on the cluster formation method. In all cases, the gas is anti-biased compared with the dark matter. However, clusters formed hierarchically are found to be more depleted in baryons than clusters formed non-hierarchically. There is a depletion of 10 to 15% for hierarchically formed clusters while the depletion is less than 10% for those clusters formed non-hierarchically. This difference is dependent on the mass of the clusters. The mean baryon enrichment profile for the hierarchically formed clusters shows an appreciable baryon enhancement around the virial radius not seen in the clusters formed without substructure. If this phenomenon applies to real clusters, it implies that determinations of the baryon fractions in clusters of galaxies require data extending beyond the virial radius of the clusters in order to achieve an unbiased value.

The significance of hierarchical clustering on the density profile and mass-temperature scaling relation for galaxy clusters is examined using hydrodynamic N-body simulations. Clusters formed hierarchically are compared with clusters formed with the initial density fluctuations on sub-cluster scales removed via smoothing. The universal profile, as described by Navarro, Frenk, and White, is not a by-product of hierarchical clustering. It is found to fit the mean profiles of clusters formed both hierarchically and otherwise. The Hernquist profile is also found to fit the data well. The characteristic radius, r_s, moves outward from 0.1 R_200 to 0.2 R_200 when the initial substructure is eliminated. Interior to r_s, rho_DM is proportional to r^-1.8, regardless of initial smoothing. Exterior to this radius, the profile is marginally shallower in the non-hierarchical case, with rho_DM proportional to r^-2.4 compared with rho_DM proportional to r^-2.7. The mass-temperature scaling relation maintains the form T proportional to M^2/3, regardless of cluster formation method. The normalisation varies at the 20% level, which is at the level of the intrinsic scatter, with the non-hierarchical simulations producing the cooler clusters.

We present a study of hydrodynamic drag forces in smoothed particle simulations. In particular, the deceleration of a resolution-limited cold clump of gas moving through a hot medium is examined. It is found that the drag for subsonic velocities exceeds that predicted by simple physical approximations. This is shown to be a result of the hydrodynamical method which encourages the accretion of particles from the hot medium onto a shell around the cold clump, effectively increasing the radius of the clump. For sonic and supersonic velocities, the drag is shown to be dependent on the effective cross section of the clump. The consequences for cosmological simulations are discussed.

We analyse the performance of twelve different implementations of Smoothed Particle Hydrodynamics (SPH) using seven tests designed to isolate key hydrodynamic elements of cosmological simulations which are known to cause the SPH algorithm problems. In order, we consider a shock tube, spherical adiabatic collapse, cooling flow model, drag, a cosmological simulation, rotating cloud-collapse and disc stability. In the implementations special attention is given to the way in which force symmetry is enforced in the equations of motion. We study in detail how the hydrodynamics are affected by different implementations of the artificial viscosity including those with a shear-correction modification. We present an improved first-order smoothing-length update algorithm that is designed to remove instabilities that are present in the Hernquist and Katz (1989) algorithm. For all tests we find that the artificial viscosity is the most important factor distinguishing the results from the various implementations. The second most important factor is the way force symmetry is achieved in the equation of motion. Most results favour a kernel symmetrization approach. The exact method by which SPH pressure forces are included has comparatively little effect on the results. Combining the equation of motion presented in Thomas and Couchman (1992) with a modification of the Monaghan and Gingold (1983) artificial viscosity leads to an SPH scheme that is both fast and reliable.