The global climate crisis poses new risks to humanity, and with them, new challenges to the practices of professional astronomy. Avoiding the more catastrophic consequences of global warming by more than 1.5 degrees requires an immediate reduction of greenhouse gas emissions. According to the 2018 United Nations Intergovernmental Panel report, this will necessitate a 45% reduction of emissions by 2030 and net-zero emissions by 2050. Efforts are required at all levels, from the individual to the governmental, and every discipline must find ways to achieve these goals. This will be especially difficult for astronomy with its significant reliance on conference and research travel, among other impacts. However, our long-range planning exercises provide the means to coordinate our response on a variety of levels. We have the opportunity to lead by example, rising to the challenge rather than reacting to external constraints. We explore how astronomy can meet the challenge of a changing climate in clear and responsible ways, such as how we set expectations (for ourselves, our institutions, and our granting agencies) around scientific travel, the organization of conferences, and the design of our infrastructure. We also emphasize our role as reliable communicators of scientific information on a problem that is both human and planetary in scale.

We perform adaptive mesh refinement (AMR) and smoothed particle hydrodynamics (SPH) cosmological zoom simulations of a region around a forming galaxy cluster, comparing the ability of the methods to handle successively more complex baryonic physics. In the simplest, non-radiative case, the two methods are in good agreement with each other, but the SPH simulations generate central cores with slightly lower entropies and virial shocks at slightly larger radii, consistent with what has been seen in previous studies. The inclusion of radiative cooling, star formation, and stellar feedback leads to much larger differences between the two methods. Most dramatically, at z=5, rapid cooling in the AMR case moves the accretion shock well within the virial radius, while this shock remains near the virial radius in the SPH case, due to excess heating, coupled with poorer capturing of the shock width. On the other hand, the addition of feedback from active galactic nuclei (AGN) to the simulations results in much better agreement between the methods. In this case both simulations display halo gas entropies of 100 keV cm^2, similar decrements in the star-formation rate, and a drop in the halo baryon content of roughly 30%. This is consistent with AGN growth being self-regulated, regardless of the numerical method. However, the simulations with AGN feedback continue to differ in aspects that are not self-regulated, such that in SPH a larger volume of gas is impacted by feedback, and the cluster still has a lower entropy central core.

Building on the initial results of the nIFTy simulated galaxy cluster comparison, we compare and contrast the impact of baryonic physics with a single massive galaxy cluster, run with 11 state-of-the-art codes, spanning adaptive mesh, moving mesh, classic and modern SPH approaches. For each code represented we have a dark matter only (DM) and non-radiative (NR) version of the cluster, as well as a full physics (FP) version for a subset of the codes. We compare both radial mass and kinematic profiles, as well as global measures of the cluster (e.g. concentration, spin, shape), in the NR and FP runs with that in the DM runs. Our analysis reveals good consistency (<= 20%) between global properties of the cluster predicted by different codes when integrated quantities are measured within the virial radius R200. However, we see larger differences for quantities within R2500, especially in the FP runs. The radial profiles reveal a diversity, especially in the cluster centre, between the NR runs, which can be understood straightforwardly from the division of codes into classic SPH and non-classic SPH (including the modern SPH, adaptive and moving mesh codes); and between the FP runs, which can also be understood broadly from the division of codes into those that include AGN feedback and those that do not. The variation with respect to the median is much larger in the FP runs with different baryonic physics prescriptions than in the NR runs with different hydrodynamics solvers.

We investigate the high-redshift evolution of the restframe UV-luminosity function (LF) of galaxies via hydrodynamical cosmological simulations, coupled with an emulated observational astronomy pipeline that provides a direct comparison with observations. We do this by creating mock images and synthetic galaxy catalogs of approximately 100 square arcminute fields from the numerical model at redshifts ~ 4.5 to 10.4. We include the effects of dust extinction and the point spread function (PSF) for the Hubble WFC3 camera for comparison with space observations. We also include the expected zodiacal background to predict its effect on space observations, including future missions such as the James Webb Space Telescope (JWST). When our model catalogs are fitted to Schechter function parameters, we predict that the faint-end slope alpha of the LF evolves as alpha = -1.16 - 0.12 z over the redshift range z ~ 4.5 to 7.7, in excellent agreement with observations from e.g., Hathi et al. (2010). However, for redshifts z ~ 6 to 10.4, alpha(z) appears to display a shallower evolution, alpha = -1.79 - 0.03 z. Augmenting the simulations with more detailed physics - specifically stellar winds and supernovae (SN) - produces similar results. The model shows an overproduction of galaxies, especially at faint magnitudes, compared with the observations, although the discrepancy is reduced when dust extinction is taken into account.

We have simulated the formation of a galaxy cluster in a $\Lambda$CDM universe using twelve different codes modeling only gravity and non-radiative hydrodynamics (\art, \arepo, \hydra\ and 9 incarnations of GADGET). This range of codes includes particle based, moving and fixed mesh codes as well as both Eulerian and Lagrangian fluid schemes. The various GADGET implementations span traditional and advanced smoothed-particle hydrodynamics (SPH) schemes. The goal of this comparison is to assess the reliability of cosmological hydrodynamical 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 $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. The underlying gravitational framework can be aligned very accurately for all the codes allowing a detailed investigation of the differences that develop due to the various gas physics implementations employed. As expected, the mesh-based codes ART and AREPO form extended entropy cores in the gas with rising central gas temperatures. Those codes employing traditional SPH schemes show falling entropy profiles all the way into the very centre with correspondingly rising density profiles and central temperature inversions. We show that methods with modern SPH schemes that allow entropy mixing span the range between these two extremes and the latest SPH variants produce gas entropy profiles that are essentially indistinguishable from those obtained with grid based methods.

We examine the correlation between the star formation rate (SFR) and black hole accretion rate (BHAR) across a suite of different AGN feedback models, using the time evolution of a merger simulation. By considering three different stages of evolution, and a distinction between the nuclear and outer regions of star formation, we consider 63 different cases. Despite many of the feedback models fitting the M-\sigma\ relationship well, there are often distinct differences in the SFR-BHAR correlations, with close to linear trends only being present after the merger. Some of the models also show evolution in the SFR-BHAR parameter space that is at times directly across the long-term averaged SFR-BHAR correlation. This suggests that the observational SFR-BHAR correlation found for ensembles of galaxies is an approximate statistical trend, as suggested by Hickox et al. Decomposing the SFR into nuclear and outer components also highlights notable differences between models and there is only modest agreement with observational studies examining this in Seyfert galaxies. For the fraction of the black hole mass growth from the merger event relative to the final black hole mass, we find as much as a factor of three variation among models. This also translates into a similar variation in the post-starburst black hole mass growth. Overall, we find that while qualitative features are often similar amongst models, precise quantitative analysis shows there can be quite distinct differences.

We examine simulations of isolated galaxies to analyse the effects of localised feedback on the formation and evolution of molecular clouds. Feedback contributes to turbulence and the destruction of clouds, leading to a population of clouds that is younger, less massive, and with more retrograde rotation. We investigate the evolution of clouds as they interact with each other and the diffuse ISM, and determine that the role of cloud interactions differs strongly with the presence of feedback: in models without feedback, scattering events dramatically increase the retrograde fraction, but in models with feedback, mergers between clouds may slightly increase the prograde fraction. We also produce an estimate of the viscous time-scale due to cloud-cloud collisions, which increases with increasing strength of feedback (~20 Gyr vs ~10 Gyr), but is still much smaller than previous estimates (~1000 Gyr); although collisions become more frequent with feedback, less energy is lost in each collision than in the models without feedback.

Observations have revealed cold gas with large velocity dispersions (~300 km/s) within the hot outflows of ultra-luminous infrared galaxies (ULIRGs). This gas may trace its origin to the Rayleigh-Taylor (RT) fragmentation of a super-bubble or may arise on smaller scales. We model a ULIRG outflow at two scales to recreate this gas in three-dimensional hydrodynamic simulations using FLASH. Although resolution is limited, these models successfully produce cold gas in outflows with large velocity dispersions. Our small-scale models produce this cold gas through RT fragmentation of the super-bubble wall, but the large-scale models produce the cold gas after hot bubbles fragment the disc's gas into cold clouds which are then accelerated by thermal pressure, and supplemented by cooling within the outflow. We produce simple mock spectra to compare these simulations to observed absorption spectra and find line-widths of ~250 km/s, agreeing with the lower end of observations.

In the early universe, substantial relative "stream" velocities between the gas and dark matter arise due to radiation pressure and persist after recombination. To asses the impact of these velocities on high-redshift structure formation, we carry out a suite of high-resolution Adaptive Mesh Refinement (AMR) cosmological simulations, which use Smoothed Particle Hydrodynamic datasets as initial conditions, converted using a new tool developed for this work. These simulations resolve structures with masses as small as a few 100 M$_\odot$, and we focus on the $10^6$ M$_\odot$ "mini-halos" in which the first stars formed. At $z \approx 17,$ the presence of stream velocities has only a minor effect on the number density of halos below $10^6$ M$_\odot$, but it greatly suppresses gas accretion onto all halos and the dark matter structures around them. Stream velocities lead to significantly lower halo gas fractions, especially for $\approx 10^5$ M$_\odot$ objects, an effect that is likely to depend on the orientation of a halo's accretion lanes. This reduction in gas density leads to colder, more compact radial profiles, and it substantially delays the redshift of collapse of the largest halos, leading to delayed star formation and possibly delayed reionization. These many differences suggest that future simulations of early cosmological structure formation should include stream velocities to properly predict gas evolution, star-formation, and the epoch of reionization.

Modelling AGN feedback in numerical simulations is both technically and theoretically challenging, with numerous approaches having been published in the literature. We present a study of five distinct approaches to modelling AGN feedback within gravitohydrodynamic simulations of major mergers of Milky Way-sized galaxies. To constrain differences to only be between AGN feedback models, all simulations start from the same initial conditions and use the same star formation algorithm. Most AGN feedback algorithms have five key aspects: black hole accretion rate, energy feedback rate and method, particle accretion algorithm, black hole advection algorithm, and black hole merger algorithm. All models follow different accretion histories, with accretion rates that differ by up to three orders of magnitude at any given time. We consider models with either thermal or kinetic feedback, with the associated energy deposited locally around the black hole. Each feedback algorithm modifies the gas properties near the black hole to different extents. The particle accretion algorithms usually (but not always) maintain good agreement between the mass accreted by \dot{M} dt and the mass of gas particles removed from the simulation. The black hole advection algorithms dampen inappropriate dragging of the black holes by two-body interactions. Advecting the black hole a limited distance based upon local mass distributions has many desirably properties. The black holes merge when given criteria are met, and we find a range of merger times for different criteria. Using the M_{BH}-\sigma relation as a diagnostic of the remnants yields three models that lie within the one-sigma scatter of the observed relation and two that fall below it. The wide variation in accretion behaviours of the models reinforces the fact that there remains much to be learnt about the evolution of galactic nuclei. (abridged)

A recent approach to simulating localized feedback from active galactic nuclei by Power et al. (2011) uses an accretion disc particle to represent both the black hole and its accretion disc. We have extrapolated and adapted this approach to simulations of Milky Way-sized galaxy mergers containing black holes and explored the impact of the various parameters in this model as well as its resolution dependence. The two key parameters in the model are an effective accretion radius, which determines the radius within which gas particles are added to the accretion disc, and a viscous time-scale which determines how long it takes for material in the accretion disc to accrete on to the black hole itself. We find that there is a limited range of permitted accretion radii and viscous time-scales, with unphysical results produced outside this range. For permitted model parameters, the nuclear regions of simulations with the same resolution follow similar evolutionary paths, producing final black hole masses that are consistent within a factor of two. When comparing the resolution dependence of the model, there is a trend towards higher resolution producing slightly lower mass black holes, but values for the two resolutions studied again agree within a factor of two. We also compare these results to two other AGN feedback algorithms found in the literature. While the evolution of the systems vary, most notably the intermediate total black hole mass, the final black hole masses differ by less than a factor of five amongst all of our models, and the remnants exhibit similar structural parameters. The implication of this accretion model is that, unlike most accretion algorithms, a decoupling of the accretion rate on to the black hole and the local gas properties is permitted and obtained; this allows for black hole growth even after feedback has prevented additional accretion events on to the disc.

We examine the properties and evolution of a simulated polar disc galaxy. This galaxy is comprised of two orthogonal discs, one of which contains old stars (old stellar disc), and the other, containing both younger stars and the cold gas (polar disc) of the galaxy. By exploring the shape of the inner region of the dark matter halo, we are able to confirm that the halo shape is a oblate ellipsoid flattened in the direction of the polar disc. We also note that there is a twist in the shape profile, where the innermost 3 kpc of the halo flattens in the direction perpendicular to the old disc, and then aligns with the polar disc out until the virial radius. This result is then compared to the halo shape inferred from the circular velocities of the two discs. We also use the temporal information of the simulation to track the system's evolution, and identify the processes which give rise to this unusual galaxy type. We confirm the proposal that the polar disc galaxy is the result of the last major merger, where the angular moment of the interaction is orthogonal to the angle of the infalling gas. This merger is followed by the resumption of coherent gas infall. We emphasise that the disc is rapidly restored after the major merger and that after this event the galaxy begins to tilt. A significant proportion of the infalling gas comes from filaments. This infalling gas from the filament gives the gas its angular momentum, and, in the case of the polar disc galaxy, the direction of the gas filament does not change before or after the last major merger.

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

Analytic estimates of the viscous time-scale due to cloud-cloud collisions have been as high as thousands of Gyr. Consequently, cloud collisions are widely ignored as a source of viscosity in galactic disks. However, capturing the hydrodynamics of discs in simple analytic models is a challenge, both because of the wide dynamic range and importance of 2D and 3D effects. To test the validity of analytic models we present estimates for the viscous time-scale that are measured from three dimensional SPH simulations of disc formation and evolution. We have deliberately removed uncertainties associated with star-formation and feedback thereby enabling us to place lower bounds on the time-scale for this process. We also contrast collapse simulations with results from simulations of initially stable discs and examine the impact of numerical parameters and assumptions on our work, to constrain possible systematics in our estimates. We find that cloud-collision viscous time-scales are in the range of 0.6-16 Gyr, considerably shorter than previously estimated. This large discrepency can be understood in terms of how the efficiency of collisions is included in the analytical estimates. We find that the viscous time-scale only depends weakly on the number of clouds formed, and so while the viscous time-scale will increase with increasing resolution, this effect is too weak to alter our conclusions.

We describe a new algorithm for finding substructures within dark matter haloes from N-body simulations. The algorithm relies upon the fact that dynamically distinct substructures in a halo will have a {\em local} velocity distribution that differs significantly from the mean, i.e. smooth background halo. We characterize the large-scale mean field using a coarsely grained cell-based approach, while a kernel smoothing process is used to determined the local velocity distribution. Comparing the ratio of these two estimates allows us to identify particles which are strongly cluster in velocity space relative to the background and thus resident in substructure. From this population of outliers, groups are identified using a Friends-of-Friends-like approach. False positives are rejected using Poisson noise arguments. This approach does not require a search of the full phase-space structure of a halo, a non-trivial task, and is thus computationally advantageous. We apply our algorithm to several test cases and show that it identifies not only subhaloes, bound overdensities in phase-space, but can recover tidal streams with a high purity. Our method can even find streams which do not appear significantly overdense in either physical or phase-space.

We present an analysis of the neutral hydrogen (HI) properties of a fully cosmological hydrodynamical dwarf galaxy, run with varying simulation parameters. As reported by Governato et al. (2010), the high resolution, high star formation density threshold version of this galaxy is the first simulation to result in the successful reproduction of a (dwarf) spiral galaxy without any associated stellar bulge. We have set out to compare in detail the HI distribution and kinematics of this simulated bulgeless disk with what is observed in a sample of nearby dwarfs. To do so, we extracted the radial gas density profiles, velocity dispersion (e.g., velocity ellipsoid, turbulence), and the power spectrum of structure within the cold interstellar medium from the simulations. The highest resolution dwarf, when using a high density star formation threshold comparable to densities of giant molecular clouds, possesses bulk characteristics consistent with those observed in nature, though the cold gas is not as radially extended as that observed in nearby dwarfs, resulting in somewhat excessive surface densities. The lines-of-sight velocity dispersion radial profiles have values that are in good agreement with observed dwarf galaxies, but due to the fact that only the streaming velocities of particles are tracked, a correction to include the thermal velocities can lead to profiles that are quite flat. The ISM power spectra of the simulations appear to possess more power on smaller spatial scales than that of the SMC. We conclude that unavoidable limitations remain due to the unresolved physics of star formation and feedback within pc-scale molecular clouds.

We present a simple physical model for populating dark matter halos with Lyman Alpha Emiiters(LAEs) and predict the physical properties of LAEs at z~3-7. The central tenet of this model is that the Ly-alpha luminosity is proportional to the star formation rate (SFR) which is directly related to the halo mass accretion rate. The only free parameter in our model is then the star-formation efficiency (SFE). An efficiency of 2.5% provides the best-fit to the Ly-alpha luminosity function (LF) at redshift z=3.1, and we use this SFE to construct Ly-alpha LFs at other redshifts. Our model reproduce the Ly-alpha LFs, stellar ages, SFR ~1-10; Msun/yr, stellar masses ~ 10^7-10^8 Msun and the clustering properties of LAEs at z~3-7. We find the spatial correlation lengths ro ~ 3-6 Mpc/h, in agreement with the observations. Finally, we estimate the field-to-field variation ~ 30% for current volume and flux limited surveys, again consistent with observations. Our results suggest that the star formation, and hence Ly-alpha emission in LAEs is powered by the accretion of new material, and that the physical properties of LAEs do not evolve significantly over a wide range of redshifts. Relating the accreted mass, rather than the total mass of halos, to the Ly-alpha luminosity of LAEs naturally gives rise to the duty cycle of LAEs.

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.

Recently, ATIC and PAMELA measured an anomalously large flux of leptonic cosmic rays which may arise from dark matter self-annihilation. While the annihilation signal predicted for a smooth halo is 10^2-10^3 times smaller than the measured excess, the signal can be boosted by the presence of subhalos. We investigate the feasibility of large boost factors using a new Monte Carlo calculation technique that is constrained by previous simulation work on halo substructure. The model accounts for the observed decrease in the amount of substructure with decreasing halo mass and the scatter in halo structural parameters such as the density concentration parameter. Our results suggest that boost factors of the order of 10^2 are ruled out at more than 14 sigma. We conclude that substructure alone, at least with commonly assumed annihilation cross-sections, cannot explain the anomalous flux measured by ATIC and PAMELA.

The first objects to arise in a cold dark matter universe present a daunting challenge for models of structure formation. In the ultra small-scale limit, CDM structures form nearly simultaneously across a wide range of scales. Hierarchical clustering no longer provides a guiding principle for theoretical analyses and the computation time required to carry out credible simulations becomes prohibitively high. To gain insight into this problem, we perform high-resolution (N=720^3 - 1584^3) simulations of an Einstein-de Sitter cosmology where the initial power spectrum is P(k) propto k^n, with -2.5 < n < -1. Self-similar scaling is established for n=-1 and n=-2 more convincingly than in previous, lower-resolution simulations and for the first time, self-similar scaling is established for an n=-2.25 simulation. However, finite box-size effects induce departures from self-similar scaling in our n=-2.5 simulation. We compare our results with the predictions for the power spectrum from (one-loop) perturbation theory and demonstrate that the renormalization group approach suggested by McDonald improves perturbation theory's ability to predict the power spectrum in the quasilinear regime. In the nonlinear regime, our power spectra differ significantly from the widely used fitting formulae of Peacock & Dodds and Smith et al. and a new fitting formula is presented. Implications of our results for the stable clustering hypothesis vs. halo model debate are discussed. Our power spectra are inconsistent with predictions of the stable clustering hypothesis in the high-k limit and lend credence to the halo model. Nevertheless, the fitting formula advocated in this paper is purely empirical and not derived from a specific formulation of the halo model.

We explore the dependence of the subhalo mass function on the spectral index n of the linear matter power spectrum using scale-free Einstein-de Sitter simulations with n=-1 and n=-2.5. We carefully consider finite volume effects that may call into question previous simulations of n<-2 power spectra. Subhaloes are found using a 6D friends-of-friends algorithm in all haloes originating from high-sigma peaks. For n=-1, we find that the cumulative subhalo mass function is independent of the parameters used in the subhalo finding algorithm and is consistent with the subhalo mass function found in LCDM simulations. In particular, the subhalo mass function is well fit by a power-law with an index of alpha=-0.9, that is the mass function has roughly equal mass in subhaloes per logarithmic interval in subhalo mass. Conversely, for n=-2.5, the algorithm parameters affect the subhalo mass function since subhaloes are more triaxial with less well defined boundaries. We find that the index alpha is generally larger with alpha>=-0.75. We infer that although the subhalo mass function appears to be independent of n so long as n>=-2, it begins to flatten as n->-3. Thus, the common practice of using alpha=-1.0 may greatly overestimate the number of subhaloes at the smallest scales in the CDM hierarchy.

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.

We investigate the evolution of mass density profiles in secular disk galaxy models, paying special attention to the development of a two-component profile from a single initial exponential disk free of cosmological evolution (i.e., no accretion or interactions). As the source of density profile variations, we examine the parameter space of the spin parameter, halo concentration, virial mass, disk mass and bulge mass, for a total of 162 simulations in the context of a plausible model of star formation and feedback (GADGET-2). The evolution of the galaxy mass density profile, including the development of a two-component profile with an inner and outer segment, is controlled by the ratio of the disk mass fraction, $m_{d}$, to the halo spin parameter, $\lambda$. The location of the break between the two components and speed at which it develops is directly proportional to $m_{d}/\lambda$; the amplitude of the transition between the inner and outer regions is however controlled by the ratio of halo concentration to virial velocity. The location of the divide between the inner and outer profile does not change with time. (Abridged)

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.

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.

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 carry out a detailed numerical investigation of the spatial correlation function of the initial positions of cosmological dark matter halos. In this Lagrangian coordinate system, which is especially useful for analytic studies of cosmological feedback, we are able to construct cross-correlation functions of objects with varying masses and formation redshifts and compare them with a variety of analytical approaches. For the case in which both formation redshifts are equal, we find good agreement between our numerical results and the bivariate model of Scannapieco & Barkana (2002; SB02) at all masses, redshifts, and separations, while the model of Porciani et al. (1998) does well for all parameters except for objects with different masses at small separations. We find that the standard mapping between Lagrangian and Eulerian bias performs well for rare objects at all separations, but fails if the objects are highly-nonlinear (low-sigma) peaks. In the Lagrangian case in which the formation redshifts differ, the SB02 model does well for all separations and combinations of masses, apart from a discrepancy at small separations in situations in which the smaller object is formed earlier and the difference between redshifts or masses is large. As this same limitation arises in the standard approach to the single-point progenitor distribution developed by Lacey & Cole (1993), we conclude that a more complete understanding of the progenitor distribution is the most important outstanding issue in the analytic modeling of Lagrangian bias.

It is a well-established fact that massive cosmological objects exhibit a ``geometrical bias'' that boosts their spatial correlations with respect to the underlying mass distribution. Although this geometrical bias is a simple function of mass, this is only half of the story. We show using numerical simulations that objects that are in the midst of accreting material also exhibit a ``temporal bias,'' which further boosts their clustering far above geometrical bias levels. These results may help to resolve a discrepancy between spectroscopic and clustering mass estimates of Lyman Break Galaxies, a population of high-redshift galaxies that are caught in the act of forming large numbers of new stars.

Motivated by the observed discrepancy between the strong spatial correlations of Lyman break galaxies (LBGs) and their velocity dispersions, we consider a theoretical model in which these starbursting galaxies are associated with dark matter halos that experience appreciable infall of material. We show using numerical simulation that selecting halos that substantially increase in mass within a fixed time interval introduces a ``temporal bias'' which boosts their clustering above that of the underlying population. If time intervals consistent with the observed LBGs star formation rates of 50 solar masses per year are chosen, then spatial correlations are enhanced by up to a factor of two. These values roughly correspond to the geometrical bias of objects three times as massive. Thus, it is clear that temporal biasing must be taken into account when interpreting the properties of Lyman break galaxies.

A wide range of observational and theoretical arguments suggest that the universe experienced a period of heating and metal enrichment, most likely from starbursting dwarf galaxies. Using a hydrodynamic simulation we have conducted a uniquely detailed theoretical investigation of this epoch at the end of the cosmological ``dark ages''. Outflows strip baryons from pre-viralized halos with total masses $\lesssim10{}^{10}$ M${}_\odot$, reducing their number density and the overall star formation rate, while pushing these quantities toward their observed values. We show that the metallicity of $\lesssim10{}^{10}$ M${}_\odot$ objects increases with size, but with a large scatter, reproducing the metallicity-luminosity relation of dwarf galaxies. Galaxies $\gtrsim10{}^{10}$ M${}_\odot$ form with a roughly constant initial metallicity of 10% solar, explaining the observed lack of metal-poor disk stars in these objects. Outflows enrich roughly 20% of the simulation volume, yielding a mean metallicity of 0.3% solar, in agreement with observations of CIV in QSO absorption-line systems.

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 examine the effects of galaxy outflows on the formation of dwarf galaxies in numerical simulations of the high-redshift Universe. Using a Smoothed Particle Hydrodynamic code, we conduct two detailed simulations of a (5.2 Mpc/h)^3 comoving volume of the Universe. In both simulations we implement simple, well-motivated models of galaxy identification and star formation, while our second simulation also includes a simple ``blow-out'' model of galaxy outflows in which supernova driven winds from newly formed disk galaxies punch-out and shock the intergalactic medium while leaving the host galaxies intact. A direct comparison between these simulations suggests that there are two major mechanisms by which outflows affect dwarf formation. Firstly, the formation of an outflow slows down the further accretion of gas onto a galaxy, causing an overall decrease of approximately 50% in the total gas mass accreted by the objects in our simulations. Additionally, our simulations uncover a significant population of approximately 10^9 solar mass objects whose formation is suppressed by the mechanism of ``baryonic stripping,'' in which outflows from early galaxies strip the gas out of nearby overdense regions that would have otherwise later formed into dwarf galaxies. This mechanism may be important in explaining the observed discrepancy between the number of dwarf galaxies predicted and observed in the local group and provide a natural explanation for the formation of empty halos which may be required by the existence of the extremely gas-poor extra-galactic High-Velocity Clouds.

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