Far-Ultraviolet (FUV) radiation greatly exceeds ultraviolet, supernovae and winds in the energy budget of young star clusters but is poorly modelled in galaxy simulations. We present results of the first full isolated galaxy disk simulations to include FUV radiation self-consistently. This is the first science application of the TREVR radiative transfer algorithm. We find that FUV radiation alone cannot regulate star formation. However, FUV radiation is crucial for producing warm neutral gas. FUV is also a long-range feedback and is more important in the outer disks of galaxies. We also use the super-bubble feedback model, which depends only on the supernova energy per stellar mass, is more physically realistic than common, parameter-driven alternatives and thus better constrains supernova feedback impacts. FUV and supernovae together can regulate star formation, reproducing the tight Kennicutt-Schmidt relation observed for an individual galaxy, without producing too much hot ionized medium and with less disruption to the ISM compared to supernovae alone.

We present a new method to study the characteristic scales of collapse and fragmentation in galactic disks. Clump formation is seeded in simulations via controlled perturbations with a specified wavelength and velocity. These are applied to otherwise quiet gas disks ranging from analogues of present day spirals to gas-rich, high-redshift galaxies. The results are compared to linear theory, turbulently perturbed disks and observations. The results reflect the expectations of linear, non-axisymmetric theory with a finite window for growth into a bound clump. We identify two new modes of clump formation: rotation-driven fission and fragmentation of tidal tails, though both are expected to rarely contribute to clump formation in observed disks. We find that bound clumps are generally much smaller than the so-called Toomre mass. The preferred scale for fragmentation increases with the disk gas mass but cannot produce bound objects larger than $\sim10^9$ M$_{\odot}$. The most likely bound clump mass increases from $3\times10^6$ in low mass disks up to $5\times10^8$ M$_{\odot}$. We conclude that observed massive stellar and gaseous clumps on 1 kpc scales at high redshift are most likely aggregates of many initially distinct bound clumps.

We present TREVR (Tree-based REVerse Ray Tracing), a general algorithm for computing the radiation field, including absorption, in astrophysical simulations. TREVR is designed to handle large numbers of sources and absorbers; it is based on a tree data structure and is thus suited to codes that use trees for their gravity or hydrodynamics solvers (e.g. Adaptive Mesh Refinement). It achieves computational speed while maintaining a specified accuracy via controlled lowering of the resolution of both sources and rays from each source. TREVR computes the radiation field in order ${N\log N_{source}}$ time without absorption and order ${N \log N_{source} \log{N}}$ time with absorption. These scalings arise from merging sources of radiation according to an opening angle criterion and walking the tree structure to trace a ray to a depth that gives the chosen accuracy for absorption. The absorption-depth refinement criterion is unique to TREVR. We provide a suite of tests demonstrating the algorithm's ability to accurately compute fluxes, ionization fronts and shadows.

Efficient estimators of Fourier-space statistics for large number of objects rely on Fast Fourier Transforms (FFTs), which are affected by aliasing from unresolved small scale modes due to the finite FFT grid. Aliasing takes the form of a sum over images, each of them corresponding to the Fourier content displaced by increasing multiples of the sampling frequency of the grid. These spurious contributions limit the accuracy in the estimation of Fourier-space statistics, and are typically ameliorated by simultaneously increasing grid size and discarding high-frequency modes. This results in inefficient estimates for e.g. the power spectrum when desired systematic biases are well under per-cent level. We show that using interlaced grids removes odd images, which include the dominant contribution to aliasing. In addition, we discuss the choice of interpolation kernel used to define density perturbations on the FFT grid and demonstrate that using higher-order interpolation kernels than the standard Cloud in Cell algorithm results in significant reduction of the remaining images. We show that combining fourth-order interpolation with interlacing gives very accurate Fourier amplitudes and phases of density perturbations. This results in power spectrum and bispectrum estimates that have systematic biases below 0.01% all the way to the Nyquist frequency of the grid, thus maximizing the use of unbiased Fourier coefficients for a given grid size and greatly reducing systematics for applications to large cosmological datasets.

We have explored the outskirts of dark matter haloes out to 2.5 times the virial radius using a large sample of halos drawn from Illustris, along with a set of zoom simulations (MUGS). Using these, we make a systematic exploration of the shape profile beyond R$_{vir}$. In the mean sphericity profile of Illustris halos we identify a dip close to the virial radius, which is robust across a broad range of masses and infall rates. The inner edge of this feature may be related to the virial radius and the outer edge with the splashback radius. Due to the high halo-to-halo variation this result is visible only on average. However, in four individual halos in the MUGS sample, a decrease in the sphericity and a subsequent recovery is evident close to the splashback radius. We find that this feature persists for several Gyr, growing with the halo. This feature appears at the interface between the spherical halo density distribution and the filamentary structure in the environment. The shape feature is strongest when there is a high rate of infall, implying that the effect is due to the mixing of accreting and virializing material. The filamentary velocity field becomes rapidly mixed in the halo region inside the virial radius, with the area between this and the splashback radius serving as the transition region. We also identify a long-lasting and smoothly evolving splashback region in the radial density gradient in many of the MUGS halos.

We explore the regulation of star formation in star-forming galaxies through a suite of high-resolution isolated galaxy simulations. We use the SPH code GASOLINE, including photoelectric heating and metal cooling, which produces a multi-phase interstellar medium. We show that representative star formation and feedback sub-grid models naturally lead to a weak, sub-linear dependence between the amount of star formation and changes to star formation parameters. We incorporate these sub-grid models into an equilibrium pressure-driven regulation framework. We show that the sub-linear scaling arises as a consequence of the non-linear relationship between scale height and the effective pressure generated by stellar feedback. Thus, simulated star-formation regulation is sensitive to how well vertical structure in the ISM is resolved. Full galaxy disks experience density waves which drive locally time-dependent star formation. We develop a simple time-dependent, pressure-driven model that reproduces the response extremely well.

We explore when supernovae can (and cannot) regulate the star formation and bulge growth in galaxies based on a sample of 18 simulated galaxies. The simulations include key physics such as evaporation and conduction, neglected in prior work, and required to correctly model superbubbles resulting from stellar feedback. We show that for galaxies with virial masses $>10^{12}\;M_\odot$, supernovae alone cannot prevent excessive star formation. This failure occurs due to a shutdown of galactic winds, with wind mass loadings falling from $\eta\sim10$ to $\eta<1$. In more massive systems, this transfer of baryons to the circumgalactic medium falters earlier on and the galaxies diverge significantly from observed galaxy scaling relations and morphologies. The decreasing efficiency is simply due to a deepening potential well preventing gas escape. This implies that non-supernova feedback mechanisms must become dominant for galaxies with stellar masses greater than $\sim4\times10^{10}\;M_\odot$. The runaway growth of the central stellar bulge, strongly linked to black hole growth, suggests that feedback from active galactic nuclei is the probable mechanism. Below this mass, supernovae alone are able to produce a realistic stellar mass fraction, star formation history and disc morphology.

We explore the chemical distribution of stars in a simulated galaxy. Using simulations of the same initial conditions but with two different feedback schemes (MUGS and MaGICC), we examine the features of the age-metallicity relation (AMR), and the three-dimensional age-metallicity-[O/Fe] distribution, both for the galaxy as a whole and decomposed into disc, bulge, halo, and satellites. The MUGS simulation, which uses traditional supernova feedback, is replete with chemical substructure. This sub- structure is absent from the MaGICC simulation, which includes early feedback from stellar winds, a modified IMF and more efficient feedback. The reduced amount of substructure is due to the almost complete lack of satellites in MaGICC. We identify a significant separation between the bulge and disc AMRs, where the bulge is considerably more metal-rich with a smaller spread in metallicity at any given time than the disc. Our results suggest, however, that identifying the substructure in observations will require exquisite age resolution, on the order of 0.25 Gyr. Certain satellites show exotic features in the AMR, even forming a 'sawtooth' shape of increasing metallicity followed by sharp declines which correspond to pericentric passages. This fact, along with the large spread in stellar age at a given metallicity, compromises the use of metallicity as an age indicator, although alpha abundance provides a more robust clock at early times. This may also impact algorithms that are used to reconstruct star formation histories from resolved stellar populations, which frequently assume a monotonically-increasing AMR.

We explore the effect of galactic environment on properties of molecular clouds. Using clouds formed in a large-scale galactic disc simulation, we measure the observable properties from synthetic column density maps. We confirm that a significant fraction of unbound clouds forms naturally in a galactic disc environment and that a mixed population of bound and unbound clouds can match observed scaling relations and distributions for extragalactic molecular clouds. By dividing the clouds into inner and outer disc populations, we compare their distributions of properties and test whether there are statistically significant differences between them. We find that clouds in the outer disc have lower masses, sizes, and velocity dispersions as compared to those in the inner disc for reasonable choices of the inner/outer boundary. We attribute the differences to the strong impact of galactic shear on the disc stability at large galactocentric radii. In particular, our Toomre analysis of the disc shows a narrowing envelope of unstable masses as a function of radius, resulting in the formation of smaller, lower mass fragments in the outer disc. We also show that the star formation rate is affected by the environment of the parent cloud, and is particularly influenced by the underlying surface density profile of the gas throughout the disc. Our work highlights the strengths of using galaxy-scale simulations to understand the formation and evolution of cloud properties - and the star formation within them - in the context of their environment.

We present the first cosmological galaxy evolved using the modern smoothed particle hydrodynamics (SPH) code GASOLINE2 with superbubble feedback. We show that superbubble-driven galactic outflows powered by Type II supernovae alone can produce $\rm{L^*}$ galaxies with flat rotation curves with circular velocities $\sim 200\; \rm{km/s}$, low bulge-to-disc ratios, and stellar mass fractions that match observed values from high redshift to the present. These features are made possible by the high mass loadings generated by the evaporative growth of superbubbles. Outflows are driven extremely effectively at high redshift, expelling gas at early times and preventing overproduction of stars before $z=2$. Centrally concentrated gas in previous simulations has often lead to unrealistically high bulge to total ratios and strongly peaked rotation curves. We show that supernova-powered superbubbles alone can produce galaxies that agree well with observed properties without the need for additional feedback mechanisms or increased feedback energy. We present additional results arising from properly modelled hot feedback.

We present a new analytic estimate for the energy required to create a constant density core within a dark matter halo. Our new estimate, based on more realistic assumptions, leads to a required energy that is orders of magnitude lower than is claimed in earlier work. We define a core size based on the logarithmic slope of the dark matter density profile so that it is insensitive to the functional form used to fit observed data. The energy required to form a core depends sensitively on the radial scale over which dark matter within the cusp is redistributed within the halo. Simulations indicate that within a region of comparable size to the active star forming regions of the central galaxy that inhabits the halo, dark matter particles have their orbits radially increased by a factor of 2--3 during core formation. Thus the inner properties of the dark matter halo, such as halo concentration, and final core size, set the energy requirements. As a result, the energy cost increases slowly with halo mass as M$_{\rm{h}}^{0.3-0.7}$ for core sizes $\lesssim1$ kpc. We use the expected star formation history for a given dark matter halo mass to predict dwarf galaxy core sizes. We find that supernovae alone would create well over 4 kpc cores in $10^{10}$ M$_{\odot}$ dwarf galaxies \emph{if} 100% of the energy were transferred to dark matter particle orbits. We can directly constrain the efficiency factor by studying galaxies with known stellar content and core size, such as Fornax. We find that the efficiency of coupling between stellar feedback and dark matter orbital energy need only be at the 1% level or less to explain Fornax's 1 kpc core.

We examine gas accretion and subsequent star formation in representative galaxies from the McMaster Unbiased Galaxy Simulations (Stinson et al. 2010). Accreted gas is bimodal with a natural temperature division at $10^5$ K, near the peak of the cooling curve. Cold-mode accretion dominates inflows at early times, creating a peak in total accretion at redshift z=2-4 and declining exponentially below z$\sim$2. Hot-mode accretion peaks near z=1-2 and declines gradually. Hot-mode exceeds cold-mode accretion at z$\sim$1.8 for all four galaxies rather than when the galaxy reaches a characteristic mass. Cold-mode accretion can fuel immediate star formation, while hot-mode accretion preferentially builds a large, hot gas reservoir in the halo. Late-time star formation relies on reservoir gas accreted 2-8 Gyr prior. Thus, the reservoir allows the star formation rate to surpass the current overall gas accretion rate. Stellar feedback cycles gas from the interstellar medium back into the hot reservoir. Stronger feedback results in more gas cycling, gas removal in a galactic outflow and less star formation overall, enabling simulations to match the observed star formation history. For lower mass galaxies in particular, strong feedback can delay the star formation peak to z=1-2 from the accretion peak at z=2-4.

We present a new stellar feedback model that reproduces superbubbles. Superbubbles from clustered young stars evolve quite differently to individual supernovae and are substantially more efficient at generating gas motions. The essential new components of the model are thermal conduction, sub-grid evaporation and a sub-grid multi-phase treatment for cases where the simulation mass resolution is insufficient to model the early stages of the superbubble. The multi-phase stage is short compared to superbubble lifetimes. Thermal conduction physically regulates the hot gas mass without requiring a free parameter. Accurately following the hot component naturally avoids overcooling. Prior approaches tend to heat too much mass, leaving the hot ISM below $10^6$ K and susceptible to rapid cooling unless ad-hoc fixes were used. The hot phase also allows feedback energy to correctly accumulate from multiple, clustered sources, including stellar winds and supernovae. We employ high-resolution simulations of a single star cluster to show the model is insensitive to numerical resolution, unresolved ISM structure and suppression of conduction by magnetic fields. We also simulate a Milky Way analog and a dwarf galaxy. Both galaxies show regulated star formation and produce strong outflows.

We present a framework that explains the commonly observed variation in light element abundances in globular clusters. If globular clusters form in the centres of dwarf galaxies, they will be pumped onto larger orbits as star formation progresses. The potential well will only retain the moderate velocity AGB ejecta, the expected source of enrichment, but not supernova ejecta. There is no need to increase the initial cluster mass, a requirement of self-enrichment scenarios, as all the stars within the dwarf can contribute. As the clusters move through the dwarf centre they sweep up a mix of AGB ejecta and in-falling pristine gas to form a second generation of stars. The specific mix will vary in time and is thus able to explain the spread in second generation abundances observed in different clusters. The globular clusters will survive to the present day or be stripped as part of the hierarchical merging process of larger galaxies. We illustrate how this process may operate using a high-resolution simulation of a dwarf galaxy at high redshift.

The limits of available computing power have forced models for the structure of stellar halos to adopt one or both of the following simplifying assumptions: (1) stellar mass can be "painted" onto dark matter particles in progenitor satellites; (2) pure dark matter simulations that do not form a luminous galaxy can be used. We estimate the magnitude of the systematic errors introduced by these assumptions using a controlled set of stellar halo models where we independently vary whether we look at star particles or painted dark matter particles, and whether we use a simulation in which a baryonic disk galaxy forms or a matching pure dark matter simulation that does not form a baryonic disk. We find that the "painting" simplification reduces the halo concentration and internal structure, predominantly because painted dark matter particles have different kinematics than star particles even when both are buried deep in the potential well of the satellite. The simplification of using pure dark matter simulations reduces the concentration further, but increases the internal structure, and results in a more prolate stellar halo. These differences can be a factor of 1.5-7 in concentration (as measured by the half-mass radius) and 2-7 in internal density structure. Given this level of systematic uncertainty, one should be wary of overinterpreting differences between observations and the current generation of stellar halo models based on dark matter only simulations when such differences are less than an order of magnitude.

We study the effect of warm dark matter (WDM) on hydrodynamic simulations of galaxy formation as part of the Making Galaxies in a Cosmological Context (MaGICC) project. We simulate three different galaxies using three WDM candidates of 1, 2 and 5 keV and compare results with pure cold dark matter simulations. WDM slightly reduces star formation and produces less centrally concentrated stellar profiles. These effects are most evident for the 1 keV candidate but almost disappear for $m_{\mathrm{WDM}}>2$ keV. All simulations form similar stellar discs independent of WDM particle mass. In particular, the disc scale length does not change when WDM is considered. The reduced amount of star formation in the case of 1 keV particles is due to the effects of WDM on merging satellites which are on average less concentrated and less gas rich. The altered satellites cause a reduced starburst during mergers because they trigger weaker disc instabilities in the main galaxy. Nevertheless we show that disc galaxy evolution is much more sensitive to stellar feedback than it is to WDM candidate mass. Overall we find that WDM, especially when restricted to current observational constraints ($m_{\mathrm{WDM}}>2$ keV), has a minor impact on disc galaxy formation.

We present a cosmological hydrodynamical simulation of a representative volume of the Universe, as part of the Making Galaxies in a Cosmological Context (MaGICC) project. MaGICC uses a thermal implementation for supernova and early stellar feedback. This work tests the feedback model at lower resolution across a range of galaxy masses, morphologies and merger histories. The simulated sample compares well with observations of high redshift galaxies ($z \ge 2$) including the stellar mass - halo mass ($M_\star - M_h$ ) relation, the Galaxy Stellar Mass Function (GSMF) at low masses ($M_\star \lt 5 \times 10^{10} M_\odot$ ) and the number density evolution of low mass galaxies. The poor match of $M_\star - M_h$ and the GSMF at high masses ($M_\star \ge 5 \times 10^{10} M_\odot$ ) indicates supernova feedback is insufficient to limit star formation in these haloes. At $z = 0$, our model produces too many stars in massive galaxies and slightly underpredicts the stellar mass around $L_\star$ mass galaxy. Altogether our results suggest that early stellar feedback, in conjunction with supernovae feedback, plays a major role in regulating the properties of low mass galaxies at high redshift.

We analyse the structure and chemical enrichment of a Milky Way-like galaxy with a stellar mass of 2 10^{10} M_sun, formed in a cosmological hydrodynamical simulation. It is disk-dominated with a flat rotation curve, and has a disk scale length similar to the Milky Way's, but a velocity dispersion that is ~50% higher. Examining stars in narrow [Fe/H] and [\alpha/Fe] abundance ranges, we find remarkable qualitative agreement between this simulation and observations: a) The old stars lie in a thickened distribution with a short scale length, while the young stars form a thinner disk, with scale lengths decreasing, as [Fe/H] increases. b) Consequently, there is a distinct outward metallicity gradient. c) Mono-abundance populations exist with a continuous distribution of scale heights (from thin to thick). However, the simulated galaxy has a distinct and substantive very thick disk (h_z~1.5 kpc), not seen in the Milky Way. The broad agreement between simulations and observations allows us to test the validity of observational proxies used in the literature: we find in the simulation that mono-abundance populations are good proxies for single age populations (<1 Gyr) for most abundances.

We analyze the orbits of stars and dark matter particles in the halo of a disk galaxy formed in a cosmological hydrodynamical simulation. The halo is oblate within the inner ~20 kpc and triaxial beyond this radius. About 43% of orbits are short axis tubes - the rest belong to orbit families that characterize triaxial potentials (boxes, long-axis tubes and chaotic orbits), but their shapes are close to axisymmetric. We find no evidence that the self-consistent distribution function of the nearly oblate inner halo is comprised primarily of axisymmetric short-axis tube orbits. Orbits of all families, and both types of particles are highly eccentric with mean eccentricity >0.6. We find that randomly selected samples of halo stars show no substructure in "integrals of motion" space. However individual accretion events can be clearly identified in plots of metallicity versus formation time. Dynamically young tidal debris is found primarily on a single type of orbit. However, stars associated with older satellites become chaotically mixed during the formation process (possibly due to scattering by the central bulge and disk, and baryonic processes), and appear on all four types of orbits. We find that the tidal debris in cosmological hydrodynamical simulations experiences significantly more chaotic evolution than in collisionless simulations, making it much harder to identify individual progenitors using phase space coordinates alone. However by combining information on stellar ages and chemical abundances with the orbital properties of halo stars in the underlying self-consistent potential, the identification of progenitors is likely to be possible.

The high-resolution, SPH galaxies of the McMaster Unbiased Galaxy Survey (MUGS) are used to examine the satellite systems of sixteen model host galaxies. Each galaxy has a different mass, angular momentum and merger history that yield a rich set of satellite luminosity functions. With new observations of distant satellite systems, we can compare these luminosity functions to satellite systems beyond the Local Group. We find that the luminosity functions of our simulations compare well to observations when the luminosity functions are scaled according to host mass. We use the recently-found relationship between dwarf satellites and host mass in distant satellite systems (Trentham & Tully 2009) to normalize a theoretical, complete luminosity function for the Milky Way (Koposov et al. 2008). The luminosity function of satellites, expressed as a function of the host mass, is given by dN/dM_V= 3.5 x M_host^(0.91) x 10^(0.1 x M_V-10.2), where host mass is given in M_sun. The mass of a host galaxy can be used to predict the number of dwarf satellites and even when considering spiral and elliptical hosts separately this relation holds.

We introduce the Making Galaxies in a Cosmological Context (MaGICC) program of smoothed particle hydrodynamics (SPH) simulations. We describe a parameter study of galaxy formation simulations of an L* galaxy that uses early stellar feedback combined with supernova feedback to match the stellar mass--halo mass relationship. While supernova feedback alone can reduce star formation enough to match the stellar mass--halo mass relationship, the galaxy forms too many stars before z=2 to match the evolution seen using abundance matching. Our early stellar feedback is purely thermal and thus operates like a UV ionization source as well as providing some additional pressure from the radiation of massive, young stars. The early feedback heats gas to >10^6 K before cooling to 10^4 K. The pressure from this hot gas creates a more extended disk and prevents more star formation prior to z=1 than supernovae feedback alone. The resulting disk galaxy has a flat rotation curve, an exponential surface brightness profile, and matches a wide range of disk scaling relationships. The disk forms from the inside-out with an increasing exponential scale length as the galaxy evolves. Overall, early stellar feedback helps to simulate galaxies that match observational results at low and high redshifts.

We explore the circumgalactic medium (CGM) of two simulated star-forming galaxies with luminosities L ~ 0.1 and 1 L* generated using the smooth particle hydrodynamic code GASOLINE. These simulations are part of the Making Galaxies In a Cosmological Context (MAGICC) program in which the stellar feedback is tuned to match the stellar mass-halo mass relationship. For comparison, each galaxy was also simulated using a 'lower feedback' (LF) model which has strength comparable to other implementations in the literature. The 'MAGICC feedback' (MF) model has a higher incidence of massive stars and an approximately two times higher energy input per supernova. Apart from the low-mass halo using LF, each galaxy exhibits a metal-enriched CGM that extends to approximately the virial radius. A significant fraction of this gas has been heated in supernova explosions in the disc and subsequently ejected into the CGM where it is predicted to give rise to substantial O VI absorption. The simulations do not yet address the question of what happens to the O VI when the galaxies stop forming stars. Our models also predict a reservoir of cool H I clouds that show strong Ly\alpha absorption to several hundred kpc. Comparing these models to recent surveys with the Hubble Space Telescope, we find that only the MF models have sufficient O VI and H I gas in the CGM to reproduce the observed distributions. In separate analyses, these same MF models also show better agreement with other galaxy observables (e.g. rotation curves, surface brightness profiles and H I gas distribution). We infer that the CGM is the dominant reservoir of baryons for galaxy haloes.

By means of high-resolution cosmological hydrodynamical simulations of Milky Way-like disc galaxies, we conduct an analysis of the associated stellar metallicity distribution functions (MDFs). After undertaking a kinematic decomposition of each simulation into spheroid and disc sub-components, we compare the predicted MDFs to those observed in the solar neighbourhood and the Galactic bulge. The effects of the star formation density threshold are visible in the star formation histories, which show a modulation in their behaviour driven by the threshold. The derived MDFs show median metallicities lower by 0.2-0.3 dex than the MDF observed locally in the disc and in the Galactic bulge. Possible reasons for this apparent discrepancy include the use of low stellar yields and/or centrally-concentrated star formation. The dispersions are larger than the one of the observed MDF; this could be due to simulated discs being kinematically hotter relative to the Milky Way. The fraction of low metallicity stars is largely overestimated, visible from the more negatively skewed MDF with respect to the observational sample. For our fiducial Milky Way analog, we study the metallicity distribution of the stars born "in situ" relative to those formed via accretion (from disrupted satellites), and demonstrate that this low-metallicity tail to the MDF is populated primarily by accreted stars. Enhanced supernova and stellar radiation energy feedback to the surrounding interstellar media of these pre-disrupted satellites is suggested as an important regulator of the MDF skewness.

We present a new model for the formation of stellar halos in dwarf galaxies. We demonstrate that the stars and star clusters that form naturally in the inner regions of dwarfs are expected to migrate from the gas rich, star forming centre to join the stellar spheroid. For dwarf galaxies, this process could be the dominant source of halo stars. The effect is caused by stellar feedback-driven bulk motions of dense gas which, by causing potential fluctuations in the inner regions of the halo, couple to all collisionless components. This effect has been demonstrated to generate cores in otherwise cuspy cold dark matter profiles and is particularly effective in dwarf galaxy haloes. It can build a stellar spheroid with larger ages and lower metallicities at greater radii without requiring an outside-in formation model. Globular cluster-type star clusters can be created in the galactic ISM and then migrate to the spheroid on 100\thinspace Myr timescales. Once outside the inner regions they are less susceptible to tidal disruption and are thus long lived; clusters on wider orbits may be easily unbound from the dwarf to join the halo of a larger galaxy during a merger. A simulated dwarf galaxy ($\text{M}_{vir}\simeq10^{9}\text{M}_{\odot}$ at $z=5$) is used to examine this gravitational coupling to dark matter and stars.

We examine the chemical properties of 5 cosmological hydrodynamical simulations of an M33-like disc galaxy which have been shown to be consistent with the morphological characteristics and bulk scaling relations expected of late-type spirals. These simulations are part of the Making Galaxies In a Cosmological Context (MaGICC) Project, in which stellar feedback is tuned to match the stellar mass -- halo mass relationship. Each realisation employed identical initial conditions and assembly histories, but differed from one another in their underlying baryonic physics prescriptions, including (a) the efficiency with which each supernova energy couples to the ISM, (b) the impact of feedback associated with massive star radiation pressure, (c) the role of the minimum shut-off time for radiative cooling of Type II SNe remnants, (d) the treatment of metal diffusion, and (e) varying the IMF. Our analysis focusses on the resulting stellar metallicity distribution functions (MDFs) in each simulated (analogous) `solar neighbourhood' and central `bulge' region. We compare the simulated MDFs' skewness, kurtosis, and dispersion (inter-quartile, inter-decile, inter-centile, and inter-tenth-percentile regions) with that of the empirical solar neighbourhood MDF and Local Group dwarfs. We find that the MDFs of the simulated discs are more negatively skewed, with higher kurtosis, than those observed locally. We can trace this difference to the simulations' tight and correlated age-metallicity relations (compared with that of the Milky Way), suggesting that these relations within `dwarf' discs might be steeper than in L* discs and/or the degree of stellar orbital re-distribution and migration inferred locally has not been captured in their entirety, at the resolution of our simulations. The important role of metal diffusion in ameliorating the over-production of extremely metal-poor stars is highlighted.

We compare the results of thirteen cosmological gasdynamical codes used to simulate the formation of a galaxy in the LCDM structure formation paradigm. The various runs differ in their hydrodynamical treatment (SPH, moving-mesh and AMR) but share the same initial conditions and adopt their latest published model of cooling, star formation and feedback. Despite the common halo assembly history, we find large code-to-code variations in the stellar mass, size, morphology and gas content of the galaxy at z=0, due mainly to the different implementations of feedback. Compared with observation, most codes tend to produce an overly massive galaxy, smaller and less gas-rich than typical spirals, with a massive bulge and a declining rotation curve. A stellar disk is discernible in most simulations, though its prominence varies widely from code to code. There is a well-defined trend between the effects of feedback and the severity of the disagreement with observation. Models that are more effective at limiting the baryonic mass of the galaxy come closer to matching observed galaxy scaling laws, but often to the detriment of the disk component. Our conclusions hold at two different numerical resolutions. Some differences can also be traced to the numerical techniques: more gas seems able to cool and become available for star formation in grid-based codes than in SPH. However, this effect is small compared to the variations induced by different feedback prescriptions. We conclude that state-of-the-art simulations cannot yet uniquely predict the properties of the baryonic component of a galaxy, even when the assembly history of its host halo is fully specified. Developing feedback algorithms that can effectively regulate the mass of a galaxy without hindering the formation of high-angular momentum stellar disks remains a challenge.

We examine radial and vertical metallicity gradients using a suite of disk galaxy simulations, supplemented with two classic chemical evolution approaches. We determine the rate of change of gradient and reconcile differences between extant models and observations within the `inside-out' disk growth paradigm. A sample of 25 disks is used, consisting of 19 from our RaDES (Ramses Disk Environment Study) sample, realised with the adaptive mesh refinement code RAMSES. Four disks are selected from the MUGS (McMaster Unbiased Galaxy Simulations) sample, generated with the smoothed particle hydrodynamics (SPH) code GASOLINE, alongside disks from Rahimi et al. (GCD+) and Kobayashi & Nakasato (GRAPE-SPH). Two chemical evolution models of inside-out disk growth were employed to contrast the temporal evolution of their radial gradients with those of the simulations. We find that systematic differences exist between the predicted evolution of radial abundance gradients in the RaDES and chemical evolution models, compared with the MUGS sample; specifically, the MUGS simulations are systematically steeper at high-redshift, and present much more rapid evolution in their gradients. We find that the majority of the models predict radial gradients today which are consistent with those observed in late-type disks, but they evolve to this self-similarity in different fashions, despite each adhering to classical `inside-out' growth. We find that radial dependence of the efficiency with which stars form as a function of time drives the differences seen in the gradients; systematic differences in the sub-grid physics between the various codes are responsible for setting these gradients. Recent, albeit limited, data at redshift z=1.5 are consistent with the steeper gradients seen in our SPH sample, suggesting a modest revision of the classical chemical evolution models may be required.

A clear prediction of the Cold Dark Matter model is the existence of cuspy dark matter halo density profiles on all mass scales. This is not in agreement with the observed rotation curves of spiral galaxies, challenging on small scales the otherwise successful CDM paradigm. In this work we employ high resolution cosmological hydro-dynamical simulations to study the effects of dissipative processes on the inner distribution of dark matter in Milky-Way like objects (M~1e12 Msun). Our simulations include supernova feedback, and the effects of the radiation pressure of massive stars before they explode as supernovae. The increased stellar feedback results in the expansion of the dark matter halo instead of contraction with respect to N-body simulations. Baryons are able to erase the dark matter cuspy distribution creating a flat, cored, dark matter density profile in the central several kpc of a massive Milky-Way like halo. The profile is well fit by a Burkert profile, with fitting parameters consistent with the observations. In addition, we obtain flat rotation curves as well as extended, exponential stellar disk profiles. While the stellar disk we obtain is still partially too thick to resemble the MW thin disk, this pilot study shows that there is enough energy available in the baryonic component to alter the dark matter distribution even in massive disc galaxies, providing a possible solution to the long standing problem of cusps vs. cores.

Encoded within the morphological structure of galaxies are clues related to their formation and evolutionary history. Recent advances pertaining to the statistics of galaxy morphology include sophisticated measures of concentration (C), asymmetry (A), and clumpiness (S). In this study, these three parameters (CAS) have been applied to a suite of simulated galaxies and compared with observational results inferred from a sample of nearby galaxies. The simulations span a range of late-type systems, with masses between ~1e10 Msun and ~1e12 Msun, and employ star formation density thresholds between 0.1 cm^-3 and 100 cm^-3. We have found that the simulated galaxies possess comparable concentrations to their real counterparts. However, the results of the CAS analysis revealed that the simulated galaxies are generally more asymmetric, and that the range of clumpiness values extends beyond the range of those observed. Strong correlations were obtained between the three CAS parameters and colour (B-V), consistent with observed galaxies. Furthermore, the simulated galaxies possess strong links between their CAS parameters and Hubble type, mostly in-line with their real counterparts.

We present a study of satellites in orbit around a high-resolution, smoothed particle hydrodynamics (SPH) galaxy simulated in a cosmological context. The simulated galaxy is approximately the same mass as the Milky Way. The cumulative number of luminous satellites at z = 0 is similar to the observed system of satellites orbiting the Milky Way although an analysis of the satellite mass function reveals an order of magnitude more dark satellites than luminous. Some of the dark subhalos are more massive than some of the luminous subhalos at z = 0. What separates luminous and dark subhalos is not their mass at z = 0, but the maximum mass the subhalos ever achieve. We study the effect of four mass loss mechanisms on the subhalos: ultraviolet (UV) ionising radiation, ram pressure stripping, tidal stripping, and stellar feedback, and compare the impact of each of these four mechanisms on the satellites. In the lowest mass subhalos, UV is responsible for the majority of the baryonic mass loss. Ram pressure stripping removes whatever mass remains from the low mass satellites. More massive subhalos have deeper potential wells and retain more mass during reionisation. However, as satellites pass near the centre of the main halo, tidal forces cause significant mass loss from satellites of all masses. Satellites that are tidally stripped from the outside can account for the luminous satellites that are lower mass than some of the dark satellites. Stellar feedback has the greatest impact on medium mass satellites that had formed stars, but lost all their gas by z = 0. Our results demonstrate that the missing satellite problem is not an intractable issue with the cold dark matter cosmology, but is rather a manifestation of baryonic processes.

We present the McMaster Unbiased Galaxy Simulations (MUGS), the first 9 galaxies of an unbiased selection ranging in total mass from 5$\times10^{11}$ M$_\odot$ to 2$\times10^{12}$ M$_\odot$ simulated using n-body smoothed particle hydrodynamics (SPH) at high resolution. The simulations include a treatment of low temperature metal cooling, UV background radiation, star formation, and physically motivated stellar feedback. Mock images of the simulations show that the simulations lie within the observed range of relations such as that between color and magnitude and that between brightness and circular velocity (Tully-Fisher). The greatest discrepancy between the simulated galaxies and observed galaxies is the high concentration of material at the center of the galaxies as represented by the centrally peaked rotation curves and the high bulge-to-total ratios of the simulations determined both kinematically and photometrically. This central concentration represents the excess of low angular momentum material that long has plagued morphological studies of simulated galaxies and suggests that higher resolutions and a more accurate description of feedback will be required to simulate more realistic galaxies. Even with the excess central mass concentrations, the simulations suggest the important role merger history and halo spin play in the formation of disks.

It has been hypothesized that the cosmic microwave background (CMB) provides a temperature floor for collapsing protostars that can regulate the process of star formation and result in a top-heavy initial mass function at high metallicity and high redshift. We examine whether this hypothesis has any testable observational consequences. First we determine, using a set of hydrodynamic galaxy formation simulations, that the CMB temperature floor would have influenced the majority of stars formed at redshifts between z=3 and 6, and probably even to higher redshift. Five signatures of CMB-regulated star formation are: (1) a higher supernova rate than currently predicted at high redshift; (2) a systematic discrepancy between direct and indirect measurements of the high redshift star formation rate; (3) a lack of surviving globular clusters that formed at high metallicity and high redshift; (4) a more rapid rise in the metallicity of cosmic gas than is predicted by current simulations; and (5) an enhancement in the abundances of alpha elements such as O and Mg at metallicities -2 < [Fe/H] < -0.5. Observations are not presently able to either confirm or rule out the presence of these signatures. However, if correct, the top-heavy IMF of high-redshift high-metallicity globular clusters could provide an explanation for the observed bimodality of their metallicity distribution.

We estimate the impact of weak lensing by strongly nonlinear cosmological structures on the cosmic microwave background. Accurate calculation of large $\ell$ multipoles requires N-body simulations and ray-tracing schemes with both high spatial and temporal resolution. To this end we have developed a new code that combines a gravitational Adaptive Particle-Particle, Particle-Mesh (AP3M) solver with a weak lensing evaluation routine. The lensing deviations are evaluated while structure evolves during the simulation so that all evolution steps--rather than just a few outputs--are used in the lensing computations. The new code also includes a ray-tracing procedure that avoids periodicity effects in a universe that is modeled as a 3-D torus in the standard way. Results from our new simulations are compared with previous ones based on Particle-Mesh simulations. We also systematically investigate the impact of box volume, resolution, and ray-tracing directions on the variance of the computed power spectra. We find that a box size of $512 h^{-1}$ Mpc is sufficient to provide a robust estimate of the weak lensing angular power spectrum in the $\ell$-interval (2,000--7,000). For a reaslistic cosmological model the power $[\ell(\ell+1)C_{\ell}/2\pi]^{1/2}$ takes on values of a few $\mu K$ in this interval, which suggests that a future detection is feasible and may explain the excess power at high $\ell$ in the BIMA and CBI observations.

We demonstrate a significant difference in the angular momentum transport properties of galactic disks between regions in which the interstellar medium is single phase or two phase. Our study is motivated by observations of HI in extended galactic disks which indicate velocity dispersions of nonthermal origin, suggesting that turbulence in the gas may be contributing significantly to the observed dispersion. To address this, we have implemented a shearing-box framework within the FLASH code. The new code was used to perform local simulations of galactic disks that incorporate differential rotation, self-gravity, vertical stratification, hydrodynamics and cooling. These simulations explore plausible mechanisms for driving turbulent motions via the thermal and self-gravitational instabilities coupling to differential rotation. Where a two-phase medium develops, gravitational angular momentum transporting stresses are much greater, creating a possible mechanism for transferring energy from galactic rotation to turbulence. In simulations where the disk conditions do not trigger the formation of a two-phase medium, it is found that perturbations to the flow damp without leading to a sustained mechanism for driving turbulence. The differing angular momentum transport properties of the single- and two-phase regimes of the disk suggest that a significant, dynamically motivated division can be drawn between the two, even when this division occurs far outside the star formation cutoff in a galactic disk.

We show that current clustering observations of quasars and luminous AGN can be explained by a merger model augmented by feedback from outflows. Using numerical simulations large enough to study clustering out to 25 comoving h^{-1} Mpc, we calculate correlation functions, biases, and correlation lengths as a function of AGN redshift and optical and X-ray luminosity. At optical wavelengths, our results match a wide range of current observations and generate predictions for future data sets. We reproduce the weak luminosity dependence of clustering over the currently well-measured range, and predict a much stronger dependence at higher luminosities. The increase in the amplitude of binary quasar clustering observed in the Sloan Digital Sky Survey (SDSS) is also reproduced and is predicted to occur at higher redshift, an effect that is due to the one halo term in the correlation function. On the other hand, our results do not match the rapid evolution of the correlation length observed in the SDSS at z\simeq 3, a discrepancy that is at least partially due to differences in the scales probed by our simulation versus this survey. In fact, we show that changing the distances sampled from our simulations can produce changes as large as 40% in the fitted correlation lengths. Finally, in the X-ray, our simulations produce correlation lengths similar to that observed in the Chandra Deep Field (CDF) North, but not the significantly larger correlation length observed in the CDF South.

Dwarf galaxies pose significant challenges for cosmological models. In particular, current models predict a dark matter density that is divergent at the center, in sharp contrast with observations which indicate an approximately constant central density core. Energy feedback, from supernova explosions and stellar winds, has been proposed as a major factor shaping the evolution of dwarf galaxies. We present detailed cosmological simulations with sufficient resolution both to model the relevant physical processes and to directly assess the impact of stellar feedback on observable properties of dwarf galaxies. We show that feedback drives large-scale, bulk motion of the interstellar gas resulting in significant gravitational potential fluctuations and a consequent reduction in the central matter density, bringing the theoretical predictions in agreement with observations.

One of the most important and poorly-understood issues in structure formation is the role of outflows driven by active galactic nuclei (AGN). Using large-scale cosmological simulations, we compute the impact of such outflows on the small-scale distribution of the cosmic microwave background (CMB). Like gravitationally-heated structures, AGN outflows induce CMB distortions both through thermal motions and peculiar velocities, by processes known as the thermal and kinetic Sunyaev-Zel'dovich (SZ) effects, respectively. For AGN outflows the thermal SZ effect is dominant, doubling the angular power spectrum on arcminute scales. But the most distinct imprint of AGN feedback is a substantial increase in the thermal SZ distortions around elliptical galaxies, post-starburst ellipticals, and quasars, which is linearly proportional to the outflow energy. While point source subtraction is difficult for quasars, we show that by appropriately stacking microwave measurements around early-type galaxies, the new generation of small-scale microwave telescopes will be able to directly measure AGN feedback at the level important for current theoretical models.

(Abridged) We explore the idea that the anti-hierarchical turn-off observed in the quasar population arises from self-regulating feedback, via an outflow mechanism. Using a detailed hydrodynamic simulation we calculate the luminosity function of quasars down to a redshift of z=1 in a large, cosmologically representative volume. Outflows are included explicitly by tracking halo mergers and driving shocks into the surrounding intergalactic medium. Our results are in excellent agreement with measurements of the spatial distribution of quasars, and we detect an intriguing excess of galaxy-quasar pairs at very short separations. We also reproduce the anti-hierarchical turnoff in the quasar luminosity function, however, the magnitude of the turn-off falls short of that observed as well as that predicted by analogous semi-analytic models. The difference can be traced to the treatment of gas heating within galaxies. The simulated galaxy cluster L_X-T relationship is close to that observed for z~1 clusters, but the simulated galaxy groups at z=1 are significantly perturbed by quasar outflows, suggesting that measurements of X-ray emission in high-redshift groups could well be a "smoking gun" for the AGN heating hypothesis.

The standard cosmological model, now strongly constrained by direct observation at early epochs, is very successful in describing the structure of the evolved universe on large and intermediate scales. Unfortunately, serious contradictions remain on smaller, galactic scales. Among the major small-scale problems is a significant and persistent discrepancy between observations of nearby galaxies, which imply that galactic dark matter (DM) haloes have a density profile with a flat core, and the cosmological model, which predicts that the haloes should have divergent density (a cusp) at the centre. Here we use numerical N-body simulations to show that random bulk motions of gas in small primordial galaxies, of the magnitude expected in these systems, result in a flattening of the central DM cusp on short timescales (of order 10^8 years). Gas bulk motions in early galaxies are driven by supernova explosions which result from ongoing star formation. Our mechanism is general and would have operated in all star-forming galaxies at redshifts z>~ 10. Once removed, the cusp cannot be reintroduced during the subsequent mergers involved in the build-up of larger galaxies. As a consequence, in the present universe both small and large galaxies would have flat DM core density profiles, in agreement with observations.

We discuss the design and implementation of HYDRA_OMP a parallel implementation of the Smoothed Particle Hydrodynamics-Adaptive P3M (SPH-AP3M) code HYDRA. The code is designed primarily for conducting cosmological hydrodynamic simulations and is written in Fortran77+OpenMP. A number of optimizations for RISC processors and SMP-NUMA architectures have been implemented, the most important optimization being hierarchical reordering of particles within chaining cells, which greatly improves data locality thereby removing the cache misses typically associated with linked lists. Parallel scaling is good, with a minimum parallel scaling of 73% achieved on 32 nodes for a variety of modern SMP architectures. We give performance data in terms of the number of particle updates per second, which is a more useful performance metric than raw MFlops. A basic version of the code will be made available to the community in the near future.

Higher order cumulants of point processes, such as skew and kurtosis, require significant computational effort to calculate. The traditional counts-in-cells method implicitly requires a large amount of computation since, for each sampling sphere, a count of particles is necessary. Although alternative methods based on tree algorithms can reduce execution time considerably, such methods still suffer from shot noise when measuring moments on low amplitude signals. We present a novel method for calculating higher order moments that is based upon first top-hat filtering the point process data on to a grid. After correcting for the smoothing process, we are able to sample this grid using an interpolation technique to calculate the statistics of interest. The filtering technique also suppresses noise and allows us to calculate skew and kurtosis when the point process is highly homogeneous. The algorithm can be implemented efficiently in a shared memory parallel environment provided a data-local random sampling technique is used. The local sampling technique allows us to obtain close to optimal speed-up for the sampling process on the Alphaserver GS320 NUMA architecture.

By fitting a flexible stellar anisotropy model to the observed surface brightness and line-of-sight velocity dispersion profiles of Draco we derive a sequence of cosmologically plausible two-component (stars + dark matter) models for this galaxy. The models are consistent with all the available observations and can have either cuspy Navarro-Frenk-White or flat-cored dark matter density profiles. The dark matter halos either formed relatively recently (at z~2...7) and are massive (up to ~5x10^9 M_Sun), or formed before the end of the reionization of the universe (z~7...11) and are less massive (down to ~7x10^7 M_Sun). Our results thus support either of the two popular solutions of the "missing satellites" problem of Lambda cold dark matter cosmology - that dwarf spheroidals are either very massive, or very old. We carry out high-resolution simulations of the tidal evolution of our two-component Draco models in the potential of the Milky Way. The results of our simulations suggest that the observable properties of Draco have not been appreciably affected by the Galactic tides after 10 Gyr of evolution. We rule out Draco being a "tidal dwarf" - a tidally disrupted dwarf galaxy. Almost radial Draco orbits (with the pericentric distance <15 kpc) are also ruled out by our analysis. The case of a harmonic dark matter core can be consistent with observations only for a very limited choice of Draco orbits (with the apocentric-to-pericentric distances ratio of <2.5).

Using a large set of high resolution numerical simulations incorporating non-equilibrium molecular hydrogen chemistry and a constant source of external radiation, we study gas collapse in previously photo-ionized mini-galaxies with virial temperatures less than 10^4 K in the early universe (redshifts z=10-20). We confirm that the mechanism of positive feedback of ionizing radiation on star formation in mini-galaxies proposed by Ricotti, Gnedin, & Shull (2002) can be efficient despite a significant flux of metagalactic photo-dissociating radiation. We derive critical fluxes for the Lyman-Werner background radiation sufficient to prevent the collapse of gas in mini-galaxies as a function of the virial mass of the halo and redshift. In our model, the formation of mini-galaxies in defunct HII regions is most efficient at large redshifts (z>15) and/or for large local gas overdensity delta>10. We show that non-equilibrium chemistry plays an important dynamical role not only during the initial evolutionary phase leading to the gas becoming gravitationally unstable inside the mini-halo, but also at the advanced stages of the core collapse, resulting in efficient gas accretion in the core region. We speculate on a possible connection between our objects and metal-poor globular clusters and dwarf spheroidal galaxies.

We study the clustering properties of metals in the intergalactic medium (IGM) as traced by 619 CIV, 81 SiIV, N >= 10^12 cm^-2 and 316 MgII, and 82 FeII N >= 10^11.5 cm^-2 absorption components in 19 high signal-to-noise (60-100 per pixel), high resolution (R = 45000) quasar spectra. Over the redshift range probed (1.5-3.0), CIV and SiIV trace each other closely and their line-of-sight correlation functions exhibit a steep decline at large separations and a flatter profile below ~ 150 km s^-1, with a large overall bias. These features do not depend on column depth. Carrying out a detailed SPH simulation (2 X 320^3, 57 Mpc^3 comoving), we show that this behavior can not be reproduced by models in which the IGM metallicity is constant or a local function of density. However, the CIV correlation function is consistent with a model in which metals are confined within bubbles with a typical radius Rs = 2 comoving Mpc about sources of mass >= Ms = 10^12 solar masses at z=3. Our lower redshift (0.5-2) measurements of the MgII and FeII correlation functions also uncover a steep decline at large separations and a flatter profile at small separations, but the clustering is even higher, and the turn-over is shifted to ~ 75 km s^-1. Again these features do not change with column depth. We describe an analytical bubble model for these species, which come from regions that are too compact to be simulated numerically, deriving best-fit values of R_s ~ 2.4 Mpc and M_s ~ 10^12 solar masses. Equally good fits to all four species are found in a similarly biased high-redshift enrichment model in which metals are placed within 2.4 comoving Mpc of 3 x 10^9 solar mass sources at z = 7.5.

The cold dark matter model has become the leading theoretical paradigm for the formation of structure in the Universe. Together with the theory of cosmic inflation, this model makes a clear prediction for the initial conditions for structure formation and predicts that structures grow hierarchically through gravitational instability. Testing this model requires that the precise measurements delivered by galaxy surveys can be compared to robust and equally precise theoretical calculations. Here we present a novel framework for the quantitative physical interpretation of such surveys. This combines the largest simulation of the growth of dark matter structure ever carried out with new techniques for following the formation and evolution of the visible components. We show that baryon-induced features in the initial conditions of the Universe are reflected in distorted form in the low-redshift galaxy distribution, an effect that can be used to constrain the nature of dark energy with next generation surveys.

We propose a simple model for the formation of dwarf spheroidal galaxies, in which stars are assumed to have formed from isothermal gas in hydrostatic equilibrium inside extended dark matter halos. After expelling the leftover gas, the stellar system undergoes a dynamical relaxation inside the dark matter halo. These models can adequately describe the observed properties of three (Draco, Sculptor, and Carina) out of four Galactic dwarf spheroidal satellites studied in this paper. We suggest that the fourth galaxy (Fornax), which cannot be fitted well with our model, is observed all the way to its tidal radius. Our best fitting models have virial masses of ~10^9 M_Sun, halo formation redshifts consistent with the age of oldest stars in these dwarfs, and shallow inner dark matter density profiles (with slope gamma -0.5...0). The inferred temperature of gas is ~10^4 K. In our model, the "extratidal" stars observed in the vicinity of some dwarf spheroidal galaxies are gravitationally bound to the galaxies and are a part of the extended stellar halos. The inferred virial masses make Galactic dwarf spheroidals massive enough to alleviate the "missing satellites" problem of LCDM cosmologies.

We use giga-particle N-body simulations to study galaxy cluster populations in Hubble Volumes of LCDM (Omega_m=0.3, Omega_Lambda=0.7) and tCDM (Omega_m=1) world models. Mapping past light-cones of locations in the computational space, we create mock sky surveys of dark matter structure to z~1.4 over 10,000 sq deg and to z~0.5 over two full spheres. Calibrating the Jenkins mass function at z=0 with samples of ~1.5 million clusters, we show that the fit describes the sky survey counts to <~20% acccuracy over all redshifts for systems larger than poor groups (M>5e13 Msun/h). Fitting the observed local temperature function determines the ratio beta of specific thermal energies in dark matter and intracluster gas. We derive a scaling with power spectrum normalization beta \propto sigma8^{5/3}, and measure a 4% error on sigma8 arising from cosmic variance in temperature-limited cluster samples. Considering distant clusters, the LCDM model matches EMSS and RDCS X-ray-selected survey observations under economical assumptions for intracluster gas evolution. Using transformations of mass-limited cluster samples that mimic sigma8 variation, we explore SZ search expectations for a 10 sq deg survey complete above 10^{14} Msun/h. Cluster counts are shown to be extremely sensitive to sigma8 uncertainty while redshift statistics, such as the sample median, are much more stable. For LCDM, the characteristic temperature at fixed sky surface density is a weak function of redshift, implying an abundance of hot clusters at z>1. Assuming constant beta, four kT>8 keV clusters lie at z>2 and 40 kT>5 keV clusters lie at z>3 on the whole sky. Detection of Coma-sized clusters at z>1 violate LCDM at 95% confidence if their surface density exceeds 0.003 per sq deg, or 120 on the whole sky.

We calculate the cooling times at constant density for halos with virial temperatures from 100 K to 10^5 K that originate from a 3-sigma fluctuation of a CDM power spectrum in three different cosmologies. Our intention is to determine the first objects that can cool to low temperatures, but not to follow their dynamical evolution. We identify two generations of halos: those with low virial temperatures, Tvir < 9000 K that remain largely neutral, and those with larger virial temperatures that become ionized. The lower-temperature, lower-mass halos are the first to cool to 75 percent of their virial temperature. The precise temperature and mass of the first objects are dependent upon the molecular hydrogen (H2) cooling function and the cosmological model. The higher-mass halos collapse later but, in this paradigm, cool much more efficiently once they have done so, first via electronic transitions and then via molecular cooling: in fact, a greater residual ionization once the halos cool below 9000 K results in an enhanced H2 production and hence a higher cooling rate at low temperatures than for the lower-mass halos, so that within our constant-density model it is the former that are the first to cool to really low temperatures. We discuss the possible significance of this result in the context of CDM models in which the shallow slope of the initial fluctuation spectrum on small scales leads to a wide range of halo masses (of differing overdensities) collapsing over a small redshift interval. This ``crosstalk'' is sufficiently important that both high- and low-mass halos collapse during the lifetimes of the massive stars which may be formed at these epochs. Further investigation is thus required to determine which generation of halos plays the dominant role in early structure formation.

We present a smoothed particle hydrodynamic (SPH) simulation that reproduces a galaxy that is a moderate facsimile of those observed. The primary failing point of previous simulations of disk formation, namely excessive transport of angular momentum from gas to dark matter, is ameliorated by the inclusion of a supernova feedback algorithm that allows energy to persist in the model ISM for a period corresponding to the lifetime of stellar associations. The inclusion of feedback leads to a disk at a redshift $z=0.52$, with a specific angular momentum content within 10% of the value required to fit observations. An exponential fit to the disk baryon surface density gives a scale length within 17% of the theoretical value. Runs without feedback, with or without star formation, exhibit the drastic angular momentum transport observed elsewhere.

We have performed cosmological simulations in a $\Lambda$CDM cosmology with and without radiative cooling, in order to study the effect of cooling on the cluster scaling relations. Our simulations consist of 4.1 million particles each of gas and dark matter within a box-size of 100 $h^{-1}$ Mpc and the run with cooling is the largest of its kind to have been evolved to $z=0$. Our cluster catalogues both consist of over 400 objects and are complete in mass down to $\sim 10^{13} h^{-1} {\rm M_{\odot}}$. We contrast the emission-weighted temperature-mass ($T_{\rm ew}-M$) and bolometric luminosity-temperature ($L_{\rm bol}-T_{\rm ew}$) relations for the simulations at $z=0$. We find that radiative cooling \emph{increases} the temperature of intracluster gas and \emph{decreases} its total luminosity, in agreement with the results of Pearce et al. Furthermore, the temperature dependence of these effects flattens the slope of the $T_{\rm ew}-M$ relation and steepens the slope of the $L_{\rm bol}-T_{\rm ew}$ relation. Inclusion of radiative cooling in the simulations is sufficient to reproduce the observed X-ray scaling relations without requiring excessive non-gravitational energy injection.

We create mock pencil-beam redshift surveys from very large cosmological $N$-body simulations of two Cold Dark Matter cosmogonies, an Einstein-de Sitter model ($\tau$CDM) and a flat model with $\Omega_0 =0.3$ and a cosmological constant ($\Lambda$CDM). We use these to assess the significance of the apparent periodicity discovered by Broadhurst et al. (1990). Simulation particles are tagged as `galaxies' so as to reproduce observed present-day correlations. They are then identified along the past light-cones of hypothetical observers to create mock catalogues with the geometry and the distance distribution of the Broadhurst et al. data. We produce 1936 (2625) quasi-independent catalogues from our $\tau$CDM ($\Lambda$CDM) simulation. A couple of large clumps in a catalogue can produce a high peak at low wavenumbers in the corresponding one-dimensional power spectrum, without any apparent large-scale periodicity in the original redshift histogram. Although the simulated redshift histograms frequently display regularly spaced clumps, the spacing of these clumps varies between catalogues and there is no `preferred' period over our many realisations. We find only a 0.72 (0.49) per cent chance that the highest peak in the power spectrum of a $\tau$CDM ($\Lambda$CDM) catalogue has a peak-to-noise ratio higher than that in the Broadhurst et al. data. None of the simulated catalogues with such high peaks shows coherently spaced clumps with a significance as high as that of the real data. We conclude that in CDM universes, the kind of regularity observed by Broadhurst et al. has a priori probability well below $10^{-3}$.

We present results of large N-body-hydrodynamic simulations of galaxy formation. Our simulations follow the formation of galaxies in cubic volumes of side 100Mpc, in two versions of the cold dark matter (CDM) cosmogony: the standard, Omega=1 SCDM model and the flat, Omega=0.3 LCDM model. Over 2000 galaxies form in each of these simulations. We examine the rate at which gas cools and condenses into dark matter halos. This roughly tracks the cosmic star formation rate inferred from observations at various redshifts. Galaxies in the simulations form gradually over time in the hierarchical fashion characteristic of the CDM cosmogony. In the LCDM model, substantial galaxies first appear at z=5 and the population builds up rapidly until z=1 after which the rate of galaxy formation declines as cold gas is consumed and the cooling time of hot gas increases. In both cosmologies, the present-day K-band luminosity function of the simulated galaxies resembles observations. The galaxy autocorrelation functions differ significantly from those of the dark matter. At the present epoch there is little bias in either model between galaxies and dark matter on large scales, but a significant anti-bias on scales of 1Mpc/h and a positive bias on scales of 100kpc/h. The galaxy correlation function evolves little with redshift in the range z=0-3, and depends on the luminosity of the galaxy sample. The projected pairwise velocity dispersion of the galaxies is much lower than that of the dark matter on scales less than 2Mpc/h. Applying a virial mass estimator to the largest galaxy clusters recovers the cluster virial masses in an unbiased way. (Abridged)

We present a fitting formula for the non-linear evolution of the bispectrum in CDM models, obtained from measurements in high resolution numerical simulations. The formula interpolates between the perturbative and highly non-linear regimes, and generalizes previous results obtained for scale-free initial conditions.

We combine data from a number of N-body simulations to predict the abundance of dark halos in Cold Dark Matter universes over more than 4 orders of magnitude in mass. A comparison of different simulations suggests that the dominant uncertainty in our results is systematic and is smaller than 10--30% at all masses, depending on the halo definition used. In particular, our ``Hubble Volume'' simulations of \tcdm and \lcdm cosmologies allow the abundance of massive clusters to be predicted with uncertainties well below those expected in all currently planned observational surveys. We show that for a range of CDM cosmologies and for a suitable halo definition, the simulated mass function is almost independent of epoch, of cosmological parameters, and of initial power spectrum when expressed in appropriate variables. This universality is of exactly the kind predicted by the familiar Press-Schechter model, although this model predicts a mass function shape which differs from our numerical results, overestimating the abundance of ``typical'' halos and underestimating that of massive systems.

We are making available on the WWW a selection of the archived data from N-body simulations carried out by the Virgo consortium and related groups. This currently includes: (i) time-slice, lightcone and cluster data from the two $10^9$-particle Hubble volume simulations described by Evrard 1998; (ii) time-slice data from simulations of 4 different cold dark matter cosmological models with $256^3$ particles analysed by Jenkins et al 1998; (iii) Dark halo catalogs, merger trees and galaxy catalogs from the GIF project described by Kauffmann et al 1999. Basic software is supplied to read the data. The data can be accessed from: http://www.mpa-garching.mpg.de/Virgo/data_download.html

We have extracted over 400 clusters, covering more than 2 decades in mass, from three simulations of the TCDM cosmology. This represents the largest, uniform catalogue of simulated clusters ever produced. The clusters exhibit a wide variety of density-profiles. Only a minority are well-fit in their outer regions by the widely used density profile of Navarro, Frenk & White (1977). Others have steeper outer density profiles, show sharp breaks in their density profiles, or have significant substructure. If we force a fit to the NFW profile, then the best-fit concentrations decline with increasing mass, but this is driven primarily by an increase in substructure as one moves to higher masses. The measured temperature-mass relations for properties measured within a sphere enclosing a fixed overdensity all follow the self-similar form, T\propto M^{2/3}, however the normalisation is lower than in observed clusters. The temperature-mass relations for properties measured within a fixed physical radius are significantly steeper then this. Both can be accurately predicted using the NFW model.

We use very large cosmological N--body simulations to obtain accurate predictions for the two-point correlations and power spectra of mass-limited samples of galaxy clusters. We consider two currently popular cold dark matter (CDM) cosmogonies, a critical density model ($\tau$CDM) and a flat low density model with a cosmological constant ($\Lambda$CDM). Our simulations each use $10^9$ particles to follow the mass distribution within cubes of side $2h^{-1}$Gpc ($\tau$CDM) and $3h^{-1}$Gpc ($\Lambda$CDM) with a force resolution better than $10^{-4}$ of the cube side. We investigate how the predicted cluster correlations increase for samples of increasing mass and decreasing abundance. Very similar behaviour is found in the two cases. The correlation length increases from $r_0=12$ -- 13$h^{-1}$Mpc for samples with mean separation $d_{\rm c}=30h^{-1}$Mpc to $r_0=22$-- 27$h^{-1}$Mpc for samples with $d_{\rm c}=100h^{-1}$Mpc. The lower value here corresponds to $\tau$CDM and the upper to $\Lambda$CDM. The power spectra of these cluster samples are accurately parallel to those of the mass over more than a decade in scale. Both correlation lengths and power spectrum biases can be predicted to better than 10% using the simple model of Sheth, Mo & Tormen (2000). This prediction requires only the linear mass power spectrum and has no adjustable parameters. We compare our predictions with published results for the APM cluster sample. The observed variation of correlation length with richness agrees well with the models, particularly for $\Lambda$CDM. The observed power spectrum (for a cluster sample of mean separation $d_{\rm c}=31h^{-1}$Mpc) lies significantly above the predictions of both models.

In this paper, we investigate the effect of cooling on the X-ray properties of galaxy clusters. We have performed N-body, hydrodynamical simulations both with and without the effects of radiative cooling, but neglecting the effects of star formation and feedback. We show that radiative cooling produces an inflow of high-entropy gas from the outer parts of the cluster, thus \emph{raising} the cluster temperature and \emph{decreasing} the X-ray luminosity. With radiative cooling clusters are on average three to five times less luminous in X-rays than the same cluster simulated without cooling. However, we do not produce a large constant-density core in either the gas or the dark matter distributions. Our results contradict previous work in which cooling raises the X-ray luminosity and deposits an unreasonably large amount of mass in the central cluster galaxy. We achieve this by selecting our numerical resolution in such a way that a reasonable fraction of the baryonic material cools and by decoupling the hot and cold gas in our simulations, a first step towards modelling multiphase gas. We emphasise that globally cooling a sensible amount of material is vital and the presence or absence of massive central concentrations of cold baryonic material has a dramatic effect upon the resultant X-ray properties of the clusters.

We have included the chemical rate network responsible for the formation of molecular Hydrogen in the N-body hydrodynamic code, Hydra, in order to study the formation of the first cosmological at redshifts between 10 and 50. We have tested our implementation of the chemical and cooling processes by comparing N-body top hat simulations with theoretical predictions from a semi-analytic model and found them to be in good agreement. We find that post-virialization properties are insensitive to the initial abundance of molecular hydrogen. Our main objective was to determine the minimum mass ($M_{SG}(z)$) of perturbations that could become self gravitating (a prerequisite for star formation), and the redshift at which this occurred. We have developed a robust indicator for detecting the presence of a self-gravitating cloud in our simulations and find that we can do so with a baryonic particle mass-resolution of 40 solar masses. We have performed cosmological simulations of primordial objects and find that the object's mass and redshift at which they become self gravitating agree well with the $M_{SG}(z)$ results from the top hat simulations. Once a critical molecular hydrogen fractional abundance of about 0.0005 has formed in an object, the cooling time drops below the dynamical time at the centre of the cloud and the gas free falls in the dark matter potential wells, becoming self gravitating a dynamical time later.

Numerical simulations of galaxy formation require a number of parameters. Some of these are intrinsic to the numerical integration scheme (eg the timestep), while others describe the physical model (eg the gas metallicity). In this paper, we present results of a systematic exploration of the effects of varying a subset of these parameters on simulations of galaxy formation. We use N-body and ``Smoothed Particle Hydrodynamics'' techniques to follow the evolution of cold dark matter and gas in a small volume. We compare a fiducial model to 24 different simulations, in which one parameter at a time is varied, focussing on properties such as the relative fraction of hot and cold gas, and the abundance and masses of galaxies. We find that for reasonable choices of numerical values, many parameters have relatively little effect on the galaxies, with the notable exception of the parameters that control the resolution of the simulation and the efficiency with which gas cools.

We present the results of weak gravitational lensing statistics in four different cosmological $N$-body simulations. The data has been generated using an algorithm for the three-dimensional shear, which makes use of a variable softening facility for the $N$-body particle masses, and enables a physical interpretation for the large-scale structure to be made. Working in three-dimensions also allows the correct use of the appropriate angular diameter distances. Our results are presented on the basis of the filled beam approximation in view of the variable particle softening scheme in our algorithm. The importance of the smoothness of matter in the universe for the weak lensing results is discussed in some detail. The low density cosmology with a cosmological constant appears to give the broadest distributions for all the statistics computed for sources at high redshifts. In particular, the range in magnification values for this cosmology has implications for the determination of the cosmological parameters from high-redshift Type Ia Supernov\ae. The possibility of determining the density parameter from the non-Gaussianity in the probability distribution for the convergence is discussed.

We present a detailed investigation of different approaches to modeling feedback in simulations of galaxy formation. Gas-dynamic forces are evaluated using Smoothed Particle Hydrodynamics (SPH) while star formation and supernova feedback are included using a three parameter model which determines the star formation rate normalization, feedback energy and lifetime of feedback regions. The star formation rate is calculated using a Lagrangian Schmidt Law for all gas particles which satisfy temperature, density and convergent flow criteria. Feedback is incorporated as thermal heating of the ISM. We compare the effects of distributing this energy over the smoothing scale or depositing it on a single particle. Radiative losses are prevented from heated particles by adjusting the density used in radiative cooling. We test the models on the formation of galaxies from cosmological initial conditions and also on isolated Milky Way and dwarf galaxies. Extremely violent feedback is necessary to produce a gas disk with angular momentum remotely close to that of observed disk galaxies. This is a result of the extreme central concentration of the dark halos in the sCDM model, and the pervasiveness of the core-halo angular momentum transport mechanism. We emphasize that the disks formed in hierarchical simulations are partially a numerical artifact produced by the minimum mass scale of the simulation acting as a highly efficient `support' mechanism. Disk formation is strongly affected by the treatment of dense regions in SPH, which along with the difficulty of representing the hierarchical formation process, means that realistic simulations of galaxy formation require far higher resolution than currently used.

We present results from a Smoothed Particle Hydrodynamic (SPH) simulation of galaxy formation that exceeds the minimum resolution requirement suggested by Steinmetz & Muller (1993) of 3*10^4 SPH particles per galaxy. Using the multiple mass technique an effective resolution of a little over one billion particles is attained within a 48 Mpc cube. We find that even with an SPH mass resolution of 1.5*10^6 solar masses and a plausible feedback algorithm, the cooling catastrophe continues to be a problem for Einstein-de Sitter CDM cosmologies. Increasing resolution also appears to exacerbate the core-halo angular momentum transport problem.

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.

In this paper we use state-of-the-art N-body hydrodynamic simulations of a cosmological volume of side 100Mpc to produce many galaxy clusters simultaneously in both the standard cold dark matter (SCDM) cosmology and a cosmology with a positive cosmological constant (LCDM). We have performed simulations of the same volume both with and without the effects of radiative cooling, but in all cases neglect the effects of star formation and feedback. With radiative cooling clusters are on average five times less luminous in X-rays than the same cluster simulated without cooling. The importance of the mass of the central galaxy in determining the X-ray luminosity is stressed.

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 have simulated the formation of an X-ray cluster in a cold dark matter universe using 12 different codes. The codes span the range of numerical techniques and implementations currently in use, including SPH and grid methods with fixed, deformable or multilevel meshes. The goal of this comparison is to assess the reliability of cosmological gas dynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at different epochs, global properties such as mass, temperature and X-ray luminosity, and radial profiles of various dynamical and thermodynamical quantities. On the whole, the agreement among the various simulations is gratifying although a number of discrepancies exist. Agreement is best for properties of the dark matter and worst for the total X-ray luminosity. Even in this case, simulations that adequately resolve the core radius of the gas distribution predict total X-ray luminosities that agree to within a factor of two. Other quantities are reproduced to much higher accuracy. For example, the temperature and gas mass fraction within the virial radius agree to about 10%, and the ratio of specific kinetic to thermal energies of the gas agree to about 5%. Various factors contribute to the spread in calculated cluster properties, including differences in the internal timing of the simulations. Based on the overall consistency of results, we discuss a number of general properties of the cluster we have modelled.

We discuss the first results from two successful simulations of galaxy formation within a cosmological volume. With over 2000 large objects forming in each we have sufficient numbers to reliably produce both galaxy correlation and luminosity functions. We find that the observed galaxy counts are well fitted by these models and that the galaxies display an almost un-evolving correlation function back to a redshift of 3 which closely resembles the featureless observed form and amplitude.

We have simulated the growth of structure in two 100 Mpc boxes for LCDM and SCDM universes. These N-body/SPH simulations include a gaseous component which is able to cool radiatively. A fraction of the gas cools into cold dense objects which we identify as galaxies. In this article we give a preliminary analysis of the clustering behaviour of these galaxies concentrating on the LCDM model. We find a galaxy correlation function which is very close to a power law and which evolves relatively little with redshift. The pairwise dispersions of the galaxies are significantly lower than the dark matter. The LCDM model gives a surprisingly good match to the observational determinations of the galaxy correlation function and pairwise dispersions.

We discuss early results from the first large N-body/hydrodynamical simulation to resolve the formation of galaxies in a cold dark matter universe. The simulation follows the formation of galaxies by gas cooling within dark halos of mass a few times $10^{11}\Msun$ and above, in a flat universe with a positive cosmological constant. Over 2200 galaxies form in our simulated volume of $(100 \Mpc)^3$. Assigning luminosities to the model galaxies using a spectral population synthesis model results in a K-band luminosity function in excellent agreement with observations. The two-point correlation function of galaxies in the simulation evolves very little since $z=3$ and has a shape close to a power-law over four orders of magnitude in amplitude. At the present day, the galaxy correlation function in the simulation is antibiased relative to the mass on small scales and unbiased on large scales. It provides a reasonable match to observations.

We investigate the effects of weak gravitational lensing in the standard Cold Dark Matter cosmology, using an algorithm which evaluates the shear in three dimensions. The algorithm has the advantage of variable softening for the particles, and our method allows the appropriate angular diameter distances to be applied to every evaluation location within each three-dimensional simulation box. We investigate the importance of shear in the distance-redshift relation, and find it to be very small. We also establish clearly defined values for the smoothness parameter in the relation, finding its value to be at least 0.88 at all redshifts in our simulations. From our results, obtained by linking the simulation boxes back to source redshifts of 4, we are able to observe the formation of structure in terms of the computed shear, and also note that the major contributions to the shear come from a very broad range of redshifts. We show the probability distributions for the magnification, source ellipticity and convergence, and also describe the relationships amongst these quantities for a range of source redshifts. We find a broad range of magnifications and ellipticities; for sources at a redshift of 4, 97{1/2}% of all lines of sight show magnifications up to 1.3 and ellipticities up to 0.195. There is clear evidence that the magnification is not linear in the convergence, as might be expected for weak lensing, but contains contributions from higher order terms in both the convergence and the shear.

The bispectrum provides a characteristic signature of gravitational instability that can be used to probe the Gaussianity of the initial conditions and the bias of the galaxy distribution. We study how this signature is affected by redshift distortions using perturbation theory and high-resolution numerical simulations. We obtain perturbative results for the multipole expansion of the redshift-space bispectrum which provide a natural way to break the degeneracy between bias and $\Omega$ present in measurements of the redshift-space power spectrum. We propose a phenomenological model that incorporates the perturbative results and also describes the bispectrum in the transition to the non-linear regime. We stress the importance of non-linear effects and show that inaccurate treatment of these can lead to significant discrepancies in the determination of bias from galaxy redshift surveys. At small scales we find that the bispectrum monopole exhibits a strong configuration dependence that reflects the velocity dispersion of clusters. Therefore, the hierarchical model for the three-point function does not hold in redshift-space.

We have developed a new three-dimensional algorithm, based on the standard P$^3$M method, for computing deflections due to weak gravitational lensing. We compare the results of this method with those of the two-dimensional planar approach, and rigorously outline the conditions under which the two approaches are equivalent. Our new algorithm uses a Fast Fourier Transform convolution method for speed, and has a variable softening feature to provide a realistic interpretation of the large-scale structure in a simulation. The output values of the code are compared with those from the Ewald summation method, which we describe and develop in detail. With an optimal choice of the high frequency filtering in the Fourier convolution, the maximum errors, when using only a single particle, are about 7 per cent, with an rms error less than 2 per cent. For ensembles of particles, used in typical $N$-body simulations, the rms errors are typically 0.3 per cent. We describe how the output from the algorithm can be used to generate distributions of magnification, source ellipticity, shear and convergence for large-scale structure.

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.

We report on analyses of cluster samples obtained from the Hubble Volume Simulations. These simulations, an $\Omega=1$ model named $\tau$CDM and a flat low $\Omega$ model with a cosmological constant ($\Lambda$CDM), comprise the largest computational efforts to date in numerical cosmology. We investigate the presence of massive galaxy clusters at $z\approx 0.8$. The $\tau$CDM model fails to form clusters at such a redshift. However, due to the small number of observed clusters around $z\approx 0.8$ and the uncertainties in the determinations of their masses, this conclusion still is somewhat preliminary. We produce cluster catalogs at $z=0$ for both cosmologies and investigate their two--point correlation function $\xi$. We show that the relationship between the mean density of subsamples of clusters, expressed via their mean separation $d_{\rm c}$, and the correlation length $r_0$, defined through $\xi(r_0) = 1$, is not linear but turns over gently for large $d_{\rm c}$. An analytic prediction by Mo & White (1996) overpredicts $r_0$. The results from the analysis of the APM cluster data by Croft et al. (1997) are nicely matched by the $\Lambda$CDM model.

We have developed a parallel Particle-Particle, Particle-Mesh (P3M) simulation code for the Cray T3E parallel supercomputer that is well suited to studying the time evolution of systems of particles interacting via gravity and gas forces in cosmological contexts. The parallel code is based upon the public-domain serial Adaptive P3M code of Couchman et al. The algorithm resolves gravitational forces into a long-range component computed by discretizing the mass distribution and solving Poisson's equation on a grid using an FFT convolution method, and a short-range component computed by direct force summation for sufficiently close particle pairs. The code consists primarily of a particle-particle computation parallelized by domain decomposition over blocks of neighbour-cells, a more regular mesh calculation distributed in planes along one dimension, and several transformations between the two distributions. The load balancing of the P3M code is static, since this greatly aids the ongoing implementation of parallel adaptive refinements of the particle and mesh systems. Great care was taken throughout to make optimal use of the available memory, so that a version of the current implementation has been used to simulate systems of up to 10^9 particles with a 1024^3 mesh for the long-range force computation. These are the largest Cosmological N-body simulations of which we are aware. We discuss these memory optimizations as well as those motivated by computational performance. Performance results are very encouraging, and, even without refinements, the code has been used effectively for simulations in which the particle distribution becomes highly clustered as well as for other non-uniform systems of astrophysical interest.

We investigate the peculiar velocities predicted for galaxy clusters by theories in the cold dark matter family. A widely used hypothesis identifies rich clusters with high peaks of a suitably smoothed version of the linear density fluctuation field. Their peculiar velocities are then obtained by extrapolating the similarly smoothed linear peculiar velocities at the positions of these peaks. We test these ideas using large high resolution N-body simulations carried out within the Virgo supercomputing consortium. We find that at early times the barycentre of the material which ends up in a rich cluster is generally very close to a high peak of the initial density field. Furthermore the mean peculiar velocity of this material agrees well with the linear value at the peak. The late-time growth of peculiar velocities is, however, systematically underestimated by linear theory. At the time clusters are identified we find their rms peculiar velocity to be about 40% larger than predicted. Nonlinear effects are particularly important in superclusters. These systematics must be borne in mind when using cluster peculiar velocities to estimate the parameter combination $\sigma_8\Omega^{0.6}$.

We present an analysis of the clustering evolution of dark matter in four cold dark matter (CDM) cosmologies. We use a suite of high resolution, 17-million particle, N-body simulations which sample volumes large enough to give clustering statistics with unprecedented accuracy. We investigate both a flat and an open model with Omega_0=0.3, and two models with Omega=1, one with the standard CDM power spectrum and the other with the same power spectrum as the Omega_0=0.3 models. The amplitude of primordial fluctuations is set so that the models reproduce the observed abundance of rich galaxy clusters by the present day. The mass 2-point correlation function and power spectrum of all the simulations differ significantly from those of the observed galaxy distribution, in both shape and amplitude. Thus, for any of these models to provide an acceptable representation of reality, the distribution of galaxies must be biased relative to the mass in a non-trivial, scale-dependent, fashion. In the Omega=1 models the required bias is always greater than unity, but in the Omega_0=0.3 models an "antibias" is required on scales smaller than \sim 5\hmpc. The mass correlation functions in the simulations are well fit by recently published analytic models. The velocity fields are remarkably similar in all the models, whether they be characterised as bulk flows, 1-particle or pairwise velocity dispersions. This similarity is a direct consequence of our adopted normalisation. The small-scale pairwise velocity dispersion of the dark matter is somewhat larger than recent determinations from galaxy redshift surveys, but the bulk-flows predicted by our models are broadly in agreement with most available data.

We study the topology of the Virgo N-body simulations and compare it to the 1.2-Jy redshift survey of IRAS galaxies by means of the genus statistic. Four high-resolution simulations of variants of the CDM cosmology are considered: a flat standard model (SCDM), a variant of it with more large-scale power (tCDM), and two low density universes, one open (OCDM) and one flat (LCDM). The fully sampled N-body simulations are examined down to strongly nonlinear scales, both with spatially fixed smoothing, and with an adaptive smoothing technique. While the tCDM, LCDM, and OCDM simulations have very similar genus statistics in the regime accessible to fixed smoothing, they can be separated with adaptive smoothing at small mass scales. In order to compare the N-body models with the 1.2-Jy survey, we extract large ensembles of mock catalogues from the simulations. These mock surveys are used to test for systematic effects in the genus analysis and to establish the distribution of errors of the genus curve. We find that a simple multivariate analysis of the genus measurements is compromised both by non-Gaussian distributed errors and by noise that dominates the covariance matrix. We therefore introduce a principal components analysis of the genus curve. With a likelihood ratio test we find that the 1.2-Jy data favours the LCDM, tCDM and OCDM models compared to SCDM. When genus measurements for different smoothing scales are combined, the SCDM model can be excluded at a 99% confidence level, while the other three models fit the 1.2-Jy data well. (abridged)

Cosmic structure simulations have improved enormously over the past decade, both in terms of the resolution which can be achieved, and with the addition of hydrodynamic and other techniques to formerly purely gravitational methods. This is an informal, and perhaps idiosyncratic, overview of the state and progress of cosmological simulations over this period. I will discuss the strategies that are used to understand cosmic structure formation and the requirements of a successful simulation. The computational demands of cosmological simulations in general will be highlighted as will specific features of the various different algorithms that are used. The incessant push for greater resolution has lead to an increasing use of parallel computers. The desirability - or otherwise - of this trend will be discussed together with alternative techniques which sidestep the naive drive for greater resolution.

We present results of a numerical renormalization approximation to the self- similar growth of clustering of pressureless dust out of a power-law spectrum of primeval Gaussian mass density fluctuations (index n) in an Einstein-de Sitter cosmological model. The self-similar two-point correlation function, xi, seems to be well established. The renormalization solutions for xi show a satisfying insensitivity to the parameters in the method, and at n=-1 and 0 are close to the Hamilton et al. formula for interpolation between the large-scale perturbative limit and stable small-scale clustering. The solutions are tested by comparing the mean relative peculiar velocity <v_ij> of particle pairs and the velocity derived from xi under the assumption of self-similar evolution. Both the renormalization and a comparison conventional N-body solution are in reasonable agreement with the test, although the conventional approach does slightly better at large separations and the renormalization approach slightly better at small separations. Other comparisons of renormalization and conventional solutions are more demanding and the results less satisfactory. Maps of particle positions in redshift space in the renormalization solutions show more nearly empty voids and less prominent walls than do the conventional solutions. The rms relative velocity dispersion is systematically smaller in the renormalization solution. There also are sizeable differences in the frequency distributions of clump masses in the renormalization and conventional solutions. The third moment S_3 from the distribution of mass within cells is in reasonable agreement with second-order perturbation theory on large scales, while on scales less than the clustering length S_3 is roughly consistent with hierarchical clustering but is heavily affected by shot noise.

We have developed a parallel Particle-Particle, Particle-Mesh (P^3M) simulation code for the T3E well suited to studying the time evolution of systems of particles interacting via gravity and gas forces in cosmological contexts. The parallel code is based upon the public-domain serial Adaptive P^3M code of Couchman et al(1). The algorithm resolves gravitational forces into a long range component computed by discretizing the mass distribution and solving Poisson's equation on a grid using an FFT convolution method, and a short range component computed by direct force summation for sufficiently close particle pairs. The code consists primarily of a particle-particle computation parallelized by domain decomposition over blocks of neighbor-cells, a more regular mesh calculation distributed in planes along one dimension, and several transformations between the two distributions. Great care was taken throughout to make optimal use of the available memory, so that the current implementation is capable of simulating systems approaching 10^9 particles using a 1024^3 mesh for the long range force computation. These are thus among the largest N-body simulations ever carried out. We discuss these memory optimizations as well as those motivated by computational performance. Results from production runs have been very encouraging, and even prior to the implimentation of the full adaptive scheme the code has been used effectively for simulations in which the particle distribution becomes highly clustered as well as for other non-uniform systems of astrophysical interest.

The evolution of the two-point correlation function, xi(r,z), and the pairwise velocity dispersion, sigma(r,z), for both the matter and halo population, in three different cosmological models: (Omega_M,Omega_Lambda)=(1,0), (0.2,0) and (0.2,0.8) are described. If the evolution of xi is parameterized by xi(r,z)=(1+z)^{-(3+eps)}xi(r,0), where xi(r,0)=(r/r_0)^{-gamma}, then eps(mass) ranges from 1.04 +/- 0.09 for (1,0) to 0.18 +/- 0.12 for (0.2,0), as measured by the evolution of at 1 Mpc (from z ~ 5 to the present epoch). For halos, eps depends on their mean overdensity. Halos with a mean overdensity of about 2000 were used to compute the halo two-point correlation function tested with two different group finding algorithms: the friends of friends and the spherical overdensity algorithm. It is certainly believed that the rate of growth of this xihh will give a good estimate of the evolution of the galaxy two-point correlation function, at least from z ~ 1 to the present epoch. The values we get for eps(halos) range from 1.54 for (1,0) to -0.36 for (0.2,0), as measured by the evolution of xi(halos) from z ~ 1.0 to the present epoch. These values could be used to constrain the cosmological scenario. The evolution of the pairwise velocity dispersion for the mass and halo distribution is measured and compared with the evolution predicted by the Cosmic Virial Theorem (CVT). According to the CVT, sigma(r,z)^2 ~ G Q rho(z) r^2 xi(r,z) or sigma proportional to (1+z)^{-eps/2}. The values of eps measured from our simulated velocities differ from those given by the evolution of xi and the CVT, keeping gamma and Q constant: eps(CVT) = 1.78 +/- 0.13 for (1,0) or 1.40 +/- 0.28 for (0.2,0).

We investigate the internal structure of clusters of galaxies in high-resolution N-body simulations of 4 different cosmologies. There is a higher proportion of disordered clusters in critical-density than in low-density universes, although the structure of relaxed clusters is very similar in each. Crude measures of substructure, such as the shift in the position of the centre-of-mass as the density threshold is varied, can distinguish the two in a sample of just 20 or so clusters; it is harder to differentiate between clusters in open and flat models with the same density parameter. Most clusters are in a quasi-steady state within the virial radius and are well-described by the density profile of Navarro, Frenk & White (1995).

A revised version of Hydra, an adaptive particle-particle, particle-mesh plus smoothed particle hydrodynamics code developed by; Hugh Couchman, University of Western Ontario, Frazer Pearce, University of Durham, Peter Thomas, University of Sussex, is now available as a tar file from; http://coho.astro.uwo.ca/pub/hydra/hydra.html or http://star-www.maps.susx.ac.uk/~pat/hydra/hydra.html .

We describe the first parallel implementation of an adaptive particle-particle, particle-mesh code with smoothed particle hydrodynamics. Parallelisation of the serial code, ``Hydra'', is achieved by using CRAFT, a Cray proprietary language which allows rapid implementation of a serial code on a parallel machine by allowing global addressing of distributed memory. The collisionless variant of the code has already completed several 16.8 million particle cosmological simulations on a 128 processor Cray T3D whilst the full hydrodynamic code has completed several 4.2 million particle combined gas and dark matter runs. The efficiency of the code now allows parameter-space explorations to be performed routinely using $64^3$ particles of each species. A complete run including gas cooling, from high redshift to the present epoch requires approximately 10 hours on 64 processors. In this paper we present implementation details and results of the performance and scalability of the CRAFT version of Hydra under varying degrees of particle clustering.

We report on work done by the Virgo consortium, an international collaboration set up in order to study the formation and evolution of Large Scale Structure using N-body simulations on the latest generation of parallel supercomputers. We show results of 256^3 particle simulations of the formation of clusters in four Dark Matter models with different cosmological parameters. Normalizing the models such that one obtains the correct abundance of rich clusters yields an interesting result: The peculiar velocities of the clusters are almost independent of \Omega, and depend only weakly on \Gamma, the shape parameter of the power spectrum. Thus, it is nearly impossible to distinguish between high and low \Omega models on the basis of the peculiar velocities.

We report on work in progress by the Virgo consortium, a collaboration set up to carry out large simulations of the formation of galaxies and large-scale structure exploiting the latest generation of parallel supercomputers. We show results of $256^3$ particle N-body simulations of the clustering evolution of dark matter in four cold dark matter models with different cosmological parameters. The high resolution and large volume of these simulations allows us to determine reliably the mass autocorrelation function for pair separations in the range $40\hkpc$ to $20\hmpc$. Comparison of these with the observed galaxy correlation function shows that for any of these models to be viable, the distribution of galaxies must be biased relative to the distribution of mass in a non-trivial, scale-dependent fashion. In particular, low $\Omega_0$ models require the galaxies to be more ``weakly'' clustered than the mass at small and intermediate pair separations. Simulations which include the evolution of gas show that cold gas knots form with approximately the abundance expected on theoretical grounds, although a few excessively massive objects grow near the centres of rich clusters. The locations where these cold gas knots form are, in general, biased relative to the distribution of mass in a scale-dependent way. Some of these biases have the required sign but they are, for the most part, weaker than is necessary for agreement with observations. The antibias present in our low $\Omega_0$ N-body/SPH simulation appears to be related to the merging and disruption of galaxies in rich clusters.

We intend to make sets of cosmological simulations available, with a large number of different output times that may be placed side-by-side to produce a complete history of the universe stretching back to high redshift. Currently there is only one series of runs under preparation as a trial. Gauging by the response, we will consider further releases. The first series consists of N-body simulations with CDM-like spectra (Gamma=0.25) in a 100/h Mpc box, in three different cosmologies: Einstein-de Sitter, low-density flat and low-density open. The parameters of each run are: 1) Run 501: Omega=1.0, lambda=0.0, sigma_8=0.64. N=128**3. 38 time-slices back to z=5.8 (a further 11 back to z=23.4 on request). 2) Run 502: Omega=0.3, lambda=0.7, sigma_8=1.22. N=86**3. 49 time-slices back to z=3.6. 3) Run 503: Omega=0.3, lambda=0.0, sigma_8=1.06. N=86**3. 49 time-slices back to z=6.2 Check out our Web pages for more details: http://coho.astro.uwo.ca/pub/consort.html

We present an implementation of Smoothed Particle Hydrodynamics (SPH) in an adaptive-mesh PPPM algorithm. The code evolves a mixture of purely gravitational particles and gas particles. The code retains the desirable properties of previous PPPM--SPH implementations; speed under light clustering, naturally periodic boundary conditions and accurate pairwise forces. Under heavy clustering the cycle time of the new code is only 2--3 times slower than for a uniform particle distribution, overcoming the principal disadvantage of previous implementations\dash a dramatic loss of efficiency as clustering develops. A 1000 step simulation with 65,536 particles (half dark, half gas) runs in one day on a Sun Sparc10 workstation. The choice of time integration scheme is investigated in detail. A simple single-step Predictor--Corrector type integrator is most efficient. A method for generating an initial distribution of particles by allowing a a uniform temperature gas of SPH particles to relax within a periodic box is presented. The average SPH density that results varies by $\sim\pm1.3$\%. We present a modified form of the Layzer--Irvine equation which includes the thermal contribution of the gas together with radiative cooling. Tests of sound waves, shocks, spherical infall and collapse are presented. Appropriate timestep constraints sufficient to ensure both energy and entropy conservation are discussed. A cluster simulation, repeating Thomas and

Several simple mergers between model galaxy clusters containing a mixture of gas and dark matter are examined, testing the coupling of the gas to the underlying collisionless material. The gas is shocked, irreversibly dissipating the energy fed into it by the collisionless component and forms a resolved constant-density core. For the dark matter, however, admixture of phase space vacuum is not very efficient and a constant-density core is not produced. In the final state the central gas has little residual kinetic energy, indicating that streaming motions do not help to support the gas.