Radioactive isotopes produced in core-collapse supernovae (CCSNe) provide useful insights into the underlying processes driving the collapse mechanism and the origins of elemental abundances. Their study generates a confluence of major physics research, including experimental measurements of nuclear reaction rates, astrophysical modeling, and $\gamma$-ray observations. Here we identify the key nuclear reaction rates to the nucleosynthesis of observable radioactive isotopes in explosive silicon-burning during CCSNe. Using the nuclear reaction network calculator SkyNet and current REACLIB reaction rates, we evolve temperature-density-time profiles of the innermost $0.45~M_\odot$ ejecta from the core collapse and explosion of a $12~M_\odot$ star. Individually varying 3403 reaction rates by factors of 100, we identify 141 reactions which cause significant differences in the isotopes of interest, namely, $^{43}$K, $^{47}$Ca, $^{44,47}$Sc, $^{44}$Ti, $^{48,51}$Cr, $^{48,49}$V, $^{52,53}$Mn, $^{55,59}$Fe, $^{56,57}$Co, and $^{56,57,59}$Ni. For each of these reactions, we present a novel method to extract the temperature range pertinent to the nucleosynthesis of the relevant isotope; the resulting temperatures lie within the range $T = 0.47$ to $6.15~$GK. Limiting the variations to within $1\sigma$ of STARLIB reaction rate uncertainties further reduces the identified reactions to 48 key rates, which can be used to guide future experimental research. Complete results are presented in tabular form.

With the advent of modern neutrino and gravitational wave detectors, the promise of multi-messenger detections of the next galactic core-collapse supernova has become very real. Such detections will give insight into the core-collapse supernova mechanism, the structure of the progenitor star, and may resolve longstanding questions in fundamental physics. In order to properly interpret these detections, a thorough understanding of the landscape of possible core-collapse supernova events, and their multi-messenger signals, is needed. We present detailed predictions of neutrino and gravitational wave signals from 1D simulations of stellar core collapse, spanning the landscape of core-collapse progenitors from $9-120\,\mathrm{M}_{\odot}$. In order to achieve explosions in 1D, we use the STIR model, which includes the effects of turbulence and convection in 1D supernova simulations to mimic the 3D explosion mechanism. We study the gravitational wave emission from the 1D simulations using an astroseismology analysis of the proto-neutron star. We find that the neutrino and gravitational wave signals are strongly correlated with the structure of the progenitor star and remnant compact object. Using these correlations, future detections of the first few seconds of neutrino and gravitational wave emission from a galactic core-collapse supernova may be able to provide constraints on stellar evolution independent of pre-explosion imaging and the mass of the compact object remnant prior to fallback accretion.

The core-collapse supernova (CCSN) mechanism is fundamentally three-dimensional with instabilities, convection, and turbulence playing crucial roles in aiding neutrino-driven explosions. Simulations of CCNSe including accurate treatments of neutrino transport and sufficient resolution to capture key instabilities remain amongst the most expensive numerical simulations in astrophysics, prohibiting large parameter studies in 2D and 3D. Studies spanning a large swath of the incredibly varied initial conditions of CCSNe are possible in 1D, though such simulations must be artificially driven to explode. We present a new method for including the most important effects of convection and turbulence in 1D simulations of neutrino-driven CCSNe, called Supernova Turbulence In Reduced-dimensionality, or STIR. Our new approach includes crucial terms resulting from the turbulent and convective motions of the flow. We estimate the strength of convection and turbulence using a modified mixing length theory (MLT) approach introducing a few free parameters to the model which are fit to the results of 3D simulations. For sufficiently large values of the mixing length parameter, turbulence-aided neutrino-driven explosions are obtained. We compare the results of STIR to high-fidelity 3D simulations and perform a parameter study of CCSN explosion using 200 solar-metallicity progenitor models from 9 to 120 $M_\odot$. We find that STIR is a better predictor of which models will explode in multidimensional simulations than other methods of driving explosions in 1D. We also present a preliminary investigation of predicted observable characteristics of the CCSN population from STIR, such as the distributions of explosion energies and remnant masses.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

We perform a suite of smoothed particle hydrodynamics simulations to investigate in detail the results of a giant impact on the young Uranus. We study the internal structure, rotation rate, and atmospheric retention of the post-impact planet, as well as the composition of material ejected into orbit. Most of the material from the impactor's rocky core falls in to the core of the target. However, for higher angular momentum impacts, significant amounts become embedded anisotropically as lumps in the ice layer. Furthermore, most of the impactor's ice and energy is deposited in a hot, high-entropy shell at a radius of ~3 Earth radii. This could explain Uranus' observed lack of heat flow from the interior and be relevant for understanding its asymmetric magnetic field. We verify the results from the single previous study of lower resolution simulations that an impactor with a mass of at least 2 Earth masses can produce sufficiently rapid rotation in the post-impact Uranus for a range of angular momenta. At least 90% of the atmosphere remains bound to the final planet after the collision, but over half can be ejected beyond the Roche radius by a 2 or 3 Earth mass impactor. This atmospheric erosion peaks for intermediate impactor angular momenta (~3*10^36 kg m^2 s^-1). Rock is more efficiently placed into orbit and made available for satellite formation by 2 Earth mass impactors than 3 Earth mass ones, because it requires tidal disruption that is suppressed by the more massive impactors.

We present an updated general purpose nuclear equation of state (EoS) for use in simulations of core-collapse supernovae, neutron star mergers and black hole collapse. This EoS is formulated in the context of Density Functional Theory (DFT) and is generalized to include all DFT EoSs consistent with known nuclear and astrophysical constraints. This EoS also allows for the possibility of the formation of material with a net proton excess ($Y_p > 0.5$) and has an improved treatment of the nuclear statistical equilibrium and the transition to heavy nuclei as the density approaches nuclear matter density. We include the effects of pions in the regime above nuclear matter density and incorporate all of the known mesonic and baryonic states at high temperature. We analyze how a 3-body nuclear force term in the DFT at high densities stiffens the EoS to satisfy the maximum neutron star constraint, however the density dependence of the symmetry anergy and the formation of pions at high temperatures allows for a softening of the central core in supernova collapse calculations leading to a robust explosion. We also add the possibility of a transition to a QCD chiral-symmetry-restoration and deconfinement phase at densities above nuclear matter density. This paper details the physics, and constraints on, this new EoS and presents an illustration of its implementation in both neutron stars and core-collapse supernova simulations. We present the first results from core-collapse supernova simulations with this EoS.

We summarize the impact of sterile neutrino dark matter on core-collapse supernova explosions. We explore various oscillations between electron neutrinos or mixed $\mu-\tau$ neutrinos and right-handed sterile neutrinos that may occur within a core-collapse supernova. In particular, we consider sterile neutrino masses and mixing angles that are consistent with sterile neutrino dark matter candidates as indicated by recent X-ray flux measurements. We find that the interpretation of the observed 3.5 keV X-ray excess as due to a decaying 7 keV sterile neutrino that comprises 100\% of the dark matter would have almost no observable effect on supernova explosions. However, in the more realistic case in which the decaying sterile neutrino comprises only a small fraction of the total dark matter density due to the presence of other sterile neutrino flavors, WIMPs, etc., a larger mixing angle is allowed. In this case a 7 keV sterile neutrino could have a significant impact on core-collapse supernovae. We also consider mixing between $\mu-\tau$ neutrinos and sterile neutrinos. We find, however, that this mixing does not significantly alter the explosion and has no observable effect on the neutrino luminosities at early times.

We have explored the impact of sterile neutrino dark matter on core-collapse supernova explosions. We have included oscillations between electron neutrinos or mixed $\mu,\tau$ neutrinos and right-handed sterile neutrinos into a supernova model. We have chosen sterile neutrino masses and mixing angles that are consistent with sterile neutrino dark matter candidates as indicated by recent x-ray flux measurements. Using these simulations, we have explored the impact of sterile neutrinos on the core bounce and shock reheating. We find that, for ranges of sterile neutrino mass and mixing angle consistent with most dark matter constraints, the shock energy can be significantly enhanced and even a model that does not explode can be made to explode. In addition, we have found that the presence of a sterile neutrino may lead to detectable changes in the observed neutrino luminosities.

We use large N-body simulations and empirical scaling relations between dark matter halos, galaxies, and supermassive black holes to estimate the formation rates of supermassive black hole binaries and the resulting low-frequency stochastic gravitational wave background (GWB). We find this GWB to be relatively insensitive ($\lesssim10\%$) to cosmological parameters, with only slight variation between WMAP5 and Planck cosmologies. We find that uncertainty in the astrophysical scaling relations changes the amplitude of the GWB by a factor of $\sim 2$. Current observational limits are already constraining this predicted range of models. We investigate the Poisson variance in the amplitude of the GWB for randomly-generated populations of supermassive black holes, finding a scatter of order unity per frequency bin below 10 nHz, and increasing to a factor of $\sim 10$ near 100 nHz. This variance is a result of the rarity of the most massive binaries, which dominate the signal, and acts as a fundamental uncertainty on the amplitude of the underlying power law spectrum. This Poisson uncertainty dominates at $\gtrsim 20$ nHz, while at lower frequencies the dominant uncertainty is related to our poor understanding of the astrophysical scaling relations, although very low frequencies may be dominated by uncertainties related to the final parsec problem and the processes which drive binaries to the gravitational wave dominated regime. Cosmological effects are negligible at all frequencies.

We present the nonlinear 2D galaxy power spectrum, $P(k,\mu)$, in redshift space, measured from the Dark Sky simulations, using galaxy catalogs constructed with both halo occupation distribution and subhalo abundance matching methods, chosen to represent an intermediate redshift sample of luminous red galaxies. We find that the information content in individual $\mu$ (cosine of the angle to the line of sight) bins is substantially richer then multipole moments, and show that this can be used to isolate the impact of nonlinear growth and redshift space distortion (RSD) effects. Using the $\mu<0.2$ simulation data, which we show is not impacted by RSD effects, we can successfully measure the nonlinear bias to an accuracy of $\sim 5$% at $k<0.6 h$Mpc$^{-1}$. This use of individual $\mu $ bins to extract the nonlinear bias successfully removes a large parameter degeneracy when constraining the linear growth rate of structure. We carry out a joint parameter estimation, using the low $\mu$ simulation data to constrain the nonlinear bias, and $\mu\ge0.2$ to constrain the growth rate and show that $f$ can be constrained to $\sim 26\, (22)$% to a $k_{\rm max}< 0.4\, (0.6) h$Mpc$^{-1}$ from clustering alone using a simple dispersion model, for a range of galaxy models. Our analysis of individual $\mu $ bins also reveals interesting physical effects which arise simply from different methods of populating halos with galaxies. We find a prominent turnaround scale, at which RSD damping effects are greater then the nonlinear growth, which differs not only for each $\mu$ bin but also for each galaxy model. These features may provide unique signatures which could be used to shed light on the galaxy-dark matter connection.

We present a maximum probability approach to reconstructing spatial maps of the peculiar velocity field at redshifts $z\sim0.1$, where the velocities have been measured from distance indicators (DI) such as $D_n-\sigma$ relations or Tully-Fisher. With the large statistical uncertainties associated with DIs, our reconstruction method aims to recover the underlying true peculiar velocity field by reducing these errors with the use of two physically motivated filtering prior terms. The first constructs an estimate of the velocity field derived from the galaxy over-density {\deltag} and the second makes use of the matter linear density power spectrum {\pk}. Using $N$-body simulations we find, with an SDSS-like sample ($N_{gal}\simeq33$ per deg$^2$), an average correlation coefficient value of $r=0.55\pm{0.02}$ between our reconstructed velocity field and that of the true velocity field from the simulation. However, with a suitably high number density of galaxies from the next generation surveys (e.g. $N_{gal}\simeq140$ per deg$^2$) we can achieve an average $r=0.70\pm{0.02}$ out to moderate redshifts $z\sim0.1$. This will prove useful for future tests of gravity, as these relatively deep maps are complementary to weak lensing maps at the same redshift.

We have made core-collapse supernova simulations that allow oscillations between electron neutrinos (or their anti particles) with right-handed sterile neutrinos. We have considered a range of mixing angles and sterile neutrino masses including those consistent with sterile neutrinos as a dark matter candidate. We examine whether such oscillations can impact the core bounce and shock reheating in supernovae. We identify the optimum ranges of mixing angles and masses that can dramatically enhance the supernova explosion by efficiently transporting electron anti-neutrinos from the core to behind the shock where they provide additional heating leading to much larger explosion kinetic energies. We show that this effect can cause stars to explode that otherwise would have collapsed. We find that an interesting periodicity in the neutrino luminosity develops due to a cycle of depletion of the neutrino density by conversion to sterile neutrinos that shuts off the conversion, followed by a replenished neutrino density as neutrinos transport through the core.

We present VIDE, the Void IDentification and Examination toolkit, an open-source Python/C++ code for finding cosmic voids in galaxy redshift surveys and N-body simulations, characterizing their properties, and providing a platform for more detailed analysis. At its core, VIDE uses a substantially enhanced version of ZOBOV (Neyinck 2008) to calculate a Voronoi tessellation for estimating the density field and a performing a watershed transform to construct voids. Additionally, VIDE provides significant functionality for both pre- and post-processing: for example, vide can work with volume- or magnitude-limited galaxy samples with arbitrary survey geometries, or dark matter particles or halo catalogs in a variety of common formats. It can also randomly subsample inputs and includes a Halo Occupation Distribution model for constructing mock galaxy populations. VIDE uses the watershed levels to place voids in a hierarchical tree, outputs a summary of void properties in plain ASCII, and provides a Python API to perform many analysis tasks, such as loading and manipulating void catalogs and particle members, filtering, plotting, computing clustering statistics, stacking, comparing catalogs, and fitting density profiles. While centered around ZOBOV, the toolkit is designed to be as modular as possible and accommodate other void finders. VIDE has been in development for several years and has already been used to produce a wealth of results, which we summarize in this work to highlight the capabilities of the toolkit. VIDE is publicly available at http://bitbucket.org/cosmicvoids/vide public and http://www.cosmicvoids.net.

The Dark Sky Simulations are an ongoing series of cosmological N-body simulations designed to provide a quantitative and accessible model of the evolution of the large-scale Universe. Such models are essential for many aspects of the study of dark matter and dark energy, since we lack a sufficiently accurate analytic model of non-linear gravitational clustering. In July 2014, we made available to the general community our early data release, consisting of over 55 Terabytes of simulation data products, including our largest simulation to date, which used $1.07 \times 10^{12}~(10240^3)$ particles in a volume $8h^{-1}\mathrm{Gpc}$ across. Our simulations were performed with 2HOT, a purely tree-based adaptive N-body method, running on 200,000 processors of the Titan supercomputer, with data analysis enabled by yt. We provide an overview of the derived halo catalogs, mass function, power spectra and light cone data. We show self-consistency in the mass function and mass power spectrum at the 1% level over a range of more than 1000 in particle mass. We also present a novel method to distribute and access very large datasets, based on an abstraction of the World Wide Web (WWW) as a file system, remote memory-mapped file access semantics, and a space-filling curve index. This method has been implemented for our data release, and provides a means to not only query stored results such as halo catalogs, but also to design and deploy new analysis techniques on large distributed datasets.

We present and study cosmic voids identified using the watershed void finder VIDE in the Sloan Digital Sky Survey Data Release 9, compare these voids to ones identified in mock catalogs, and assess the impact of the survey mask on void statistics such as number functions, ellipticity distributions, and radial density profiles. The nearly 1,000 identified voids span three nearly volume-limited samples from redshift z = 0.43 to 0.7. For comparison we use 98 of the publicly available 2LPT-based mock galaxy catalogs of Manera et al., and also generate our own mock catalogs by applying a Halo Occupation Distribution model to an N-body simulation. We find that the mask reduces the number density of voids at all scales by a factor of three and slightly skews the relative size distributions. This engenders an increase in the mean ellipticity by roughly 30%. However, we find that radial density profiles are largely robust to the effects of the mask. We see excellent agreement between the data and both mock catalogs, and find no tension between the observed void properties and the properties derived from {\Lambda}CDM simulations. We have added the void catalogs from both data and mock galaxy populations discussed in this work to the Public Cosmic Void Catalog at http://www.cosmicvoids.net.

To study the impact of sparsity and galaxy bias on void statistics, we use a single large-volume, high-resolution N-body simulation to compare voids in multiple levels of subsampled dark matter, halo populations, and mock galaxies from a Halo Occupation Distribution model tuned to different galaxy survey densities. We focus our comparison on three key observational statistics: number functions, ellipticity distributions, and radial density profiles. We use the hierarchical tree structure of voids to interpret the impacts of sampling density and galaxy bias, and theoretical and empirical functions to describe the statistics in all our sample populations. We are able to make simple adjustments to theoretical expectations to offer prescriptions for translating from analytics to the void properties measured in realistic observations. We find that sampling density has a much larger effect on void sizes than galaxy bias. At lower tracer density, small voids disappear and the remaining voids are larger, more spherical, and have slightly steeper profiles. When a proper lower mass threshold is chosen, voids in halo distributions largely mimic those found in galaxy populations, except for ellipticities, where galaxy bias leads to higher values. We use the void density profile of Hamaus et al. (2014) to show that voids follow a self-similar and universal trend, allowing simple translations between voids studied in dark matter and voids identified in galaxy surveys. We have added the mock void catalogs used in this work to the Public Cosmic Void Catalog at http://www.cosmicvoids.net.

Galaxy bias, the unknown relationship between the clustering of galaxies and the underlying dark matter density field is a major hurdle for cosmological inference from large-scale structure. While traditional analyses focus on the absolute clustering amplitude of high-density regions mapped out by galaxy surveys, we propose a relative measurement that compares those to the underdense regions, cosmic voids. On the basis of realistic mock catalogs we demonstrate that cross correlating galaxies and voids opens up the possibility to calibrate galaxy bias and to define a static ruler thanks to the observable geometric nature of voids. We illustrate how the clustering of voids is related to mass compensation and show that volume-exclusion significantly reduces the degree of stochasticity in their spatial distribution. Extracting the spherically averaged distribution of galaxies inside voids from their cross correlations reveals a remarkable concordance with the mass-density profile of voids.

How do observed voids relate to the underlying dark matter distribution? To examine the spatial distribution of dark matter contained within voids identified in galaxy surveys, we apply Halo Occupation Distribution models representing sparsely and densely sampled galaxy surveys to a high-resolution N-body simulation. We compare these galaxy voids to voids found in the halo distribution, low-resolution dark matter, and high-resolution dark matter. We find that voids at all scales in densely sampled surveys - and medium- to large-scale voids in sparse surveys - trace the same underdensities as dark matter, but they are larger in radius by ~20%, they have somewhat shallower density profiles, and they have centers offset by ~0.4Rv rms. However, in void-to-void comparison we find that shape estimators are less robust to sampling, and the largest voids in sparsely sampled surveys suffer fragmentation at their edges. We find that voids in galaxy surveys always correspond to underdensities in the dark matter, though the centers may be offset. When this offset is taken into account, we recover almost identical radial density profiles between galaxies and dark matter. All mock catalogs used in this work are available at http://www.cosmicvoids.net.

We report on improvements made over the past two decades to our adaptive treecode N-body method (HOT). A mathematical and computational approach to the cosmological N-body problem is described, with performance and scalability measured up to 256k ($2^{18}$) processors. We present error analysis and scientific application results from a series of more than ten 69 billion ($4096^3$) particle cosmological simulations, accounting for $4 \times 10^{20}$ floating point operations. These results include the first simulations using the new constraints on the standard model of cosmology from the Planck satellite. Our simulations set a new standard for accuracy and scientific throughput, while meeting or exceeding the computational efficiency of the latest generation of hybrid TreePM N-body methods.

We present a sample of 23 spectroscopically confirmed Type Ia supernovae that were discovered in the background of galaxy clusters targeted by ROTSE-IIIb and use up to 18 of these to determine the local (z = 0.05) volumetric rate. Since our survey is flux limited and thus biased against fainter objects, the pseudo-absolute magnitude distribution (pAMD) of SNeIa in a given volume is an important concern, especially the relative frequency of high to low-luminosity SNeIa. We find that the pAMD derived from the volume limited Lick Observatory Supernova Search (LOSS) sample is incompatible with the distribution of SNeIa in a volume limited (z<0.12) sub sample of the SDSS-II. The LOSS sample requires far more low-luminosity SNeIa than the SDSS-II can accommodate. Even though LOSS and SDSS-II have sampled different SNeIa populations, their volumetric rates are surprisingly similar. Using the same model pAMD adopted in the SDSS-II SNeIa rate calculation and excluding two high-luminosity SNeIa from our sample, we derive a rate that is marginally higher than previous low-redshift determinations. With our full sample and the LOSS pAMD our rate is more than double the canonical value. We also find that 5 of our 18 SNeIa are hosted by very low-luminosity (M_B > -16) galaxies, whereas only 1 out 79 nearby SDSS-II SNeIa have such faint hosts. It is possible that previous works have under-counted either low luminosity SNeIa, SNeIa in low luminosity hosts, or peculiar SNeIa (sometimes explicitly), and the total SNeIa rate may be higher than the canonical value.

We present a galaxy catalog simulator which turns N-body simulations with subhalos into multiband photometric mocks. The simulator assigns galaxy properties to each subhalo to reproduce the observed cluster galaxy halo occupation distribution, the radial and mass dependent variation in fractions of blue galaxies, the luminosity functions in clusters and the field, and the red-sequence in clusters. Moreover, the evolution of these parameters is tuned to match existing observational constraints. Field galaxies are sampled from existing multiband photometric surveys using derived galaxy photometric redshifts. Parametrizing an ensemble of cluster galaxy properties enables us to create mock catalogs with variations in those properties, which in turn allows us to quantify the sensitivity of cluster finding to current observational uncertainties in these properties. We present an application of the catalog simulator to characterize the selection function of a galaxy cluster finder that utilizes the cluster red-sequence galaxy clustering on the sky, in terms of completeness and contamination. We estimate systematic uncertainties due to the observational uncertainties on our simulator parameters in determining the selection function using five different sets of modified catalogs. Our estimates indicate that these uncertainties are at the $\le15$% level with current observational constraints on cluster galaxy populations and their evolution. In addition, we examine the $B_{gc}$ parameter as an optical mass indicator and measure the intrinsic scatter of the $B_{gc}$--mass relation to be approximately log normal with $\sigma_{\log_{10}M}\sim0.25$. Finally, we present tests of a red sequence overdensity redshift estimator using both simulated and real data, showing that it delivers redshifts for massive clusters with $\sim$2% accuracy out to redshifts $z\sim0.5$ with SDSS-like datasets.

We measure the clustering of dark matter halos in a large set of collisionless cosmological simulations of the flat LCDM cosmology. Halos are identified using the spherical overdensity algorithm, which finds the mass around isolated peaks in the density field such that the mean density is Delta times the background. We calibrate fitting functions for the large scale bias that are adaptable to any value of Delta we examine. We find a ~6% scatter about our best fit bias relation. Our fitting functions couple to the halo mass functions of Tinker et. al. (2008) such that bias of all dark matter is normalized to unity. We demonstrate that the bias of massive, rare halos is higher than that predicted in the modified ellipsoidal collapse model of Sheth, Mo, & Tormen (2001), and approaches the predictions of the spherical collapse model for the rarest halos. Halo bias results based on friends-of-friends halos identified with linking length 0.2 are systematically lower than for halos with the canonical Delta=200 overdensity by ~10%. In contrast to our previous results on the mass function, we find that the universal bias function evolves very weakly with redshift, if at all. We use our numerical results, both for the mass function and the bias relation, to test the peak-background split model for halo bias. We find that the peak-background split achieves a reasonable agreement with the numerical results, but ~20% residuals remain, both at high and low masses.

We are developing a purely commensal survey experiment for fast (<5s) transient radio sources. Short-timescale transients are associated with the most energetic and brightest single events in the Universe. Our objective is to cover the enormous volume of transients parameter space made available by ASKAP, with an unprecedented combination of sensitivity and field of view. Fast timescale transients open new vistas on the physics of high brightness temperature emission, extreme states of matter and the physics of strong gravitational fields. In addition, the detection of extragalactic objects affords us an entirely new and extremely sensitive probe on the huge reservoir of baryons present in the IGM. We outline here our approach to the considerable challenge involved in detecting fast transients, particularly the development of hardware fast enough to dedisperse and search the ASKAP data stream at or near real-time rates. Through CRAFT, ASKAP will provide the testbed of many of the key technologies and survey modes proposed for high time resolution science with the SKA.

We investigate methods to best estimate the normalisation of the mass density fluctuation power spectrum (sigma_8) using peculiar velocity data from a survey like the Six degree Field Galaxy Velocity Survey (6dFGSv). We focus on two potential problems (i) biases from nonlinear growth of structure and (ii) the large number of velocities in the survey. Simulations of LambdaCDM-like models are used to test the methods. We calculate the likelihood from a full covariance matrix of velocities averaged in grid cells. This simultaneously reduces the number of data points and smooths out nonlinearities which tend to dominate on small scales. We show how the averaging can be taken into account in the predictions in a practical way, and show the effect of the choice of cell size. We find that a cell size can be chosen that significantly reduces the nonlinearities without significantly increasing the error bars on cosmological parameters. We compare our results with those from a principal components analysis following Watkins et al (2002) and Feldman et al (2003) to select a set of optimal moments constructed from linear combinations of the peculiar velocities that are least sensitive to the nonlinear scales. We conclude that averaging in grid cells performs equally well. We find that for a survey such as 6dFGSv we can estimate sigma_8 with less than 3% bias from nonlinearities. The expected error on sigma_8 after marginalising over Omega_m is approximately 16 percent.

The shape of the primordial matter power spectrum Plin(k) encodes critical information on cosmological parameters. At large scales, the observable galaxy power spectrum Pobs(k) is expected to follow the shape of Plin(k), but on smaller scales the effects of nonlinearity and galaxy bias make the ratio Pobs(k)/Plin(k) scale-dependent. We develop a method that can extend the dynamic range of the Plin(k) recovery by incorporating constraints on the galaxy halo occupation distribution (HOD) from the projected galaxy correlation function wp. We devise an analytic model to calculate Pobs(k) in real-space and redshift-space. Once HOD parameters are determined by matching wp for a given cosmological model, galaxy bias is completely specified, and our analytic model predicts both the shape and normalization of Pobs(k). Applying our method to SDSS main galaxy samples, we find that the real-space Pobs(k) follows the shape of the nonlinear matter power spectrum at the 1-2% level up to k=0.2 h/Mpc. When we apply our method to SDSS LRG samples, the linear bias approximation is accurate to 5% at k<0.08 h/Mpc, but the scale-dependence of LRG bias prevents the use of linear theory at k>0.08 h/Mpc. Our HOD model prediction is in good agreement with the recent SDSS LRG Pobs(k) measurements at all measured scales (k<0.2 h/Mpc), naturally explaining the shape of Pobs(k). The "Q-model" prescription is a poor description of galaxy bias for the LRG samples, and it can lead to biased cosmological parameter estimates when measurements at k>0.1 h/Mpc are included in the analysis. We quantify the potential bias and constraints on cosmological parameters that arise from applying linear theory and Q-model fitting, and we demonstrate the utility of HOD modeling of future high precision measurements of Pobs(k) on quasi-linear scales.

Galaxy cluster merger statistics are an important component in understanding the formation of large-scale structure. Unfortunately, it is difficult to study merger properties and evolution directly because the identification of cluster mergers in observations is problematic. We use large N-body simulations to study the statistical properties of massive halo mergers, specifically investigating the utility of close halo pairs as proxies for mergers. We examine the relationship between pairs and mergers for a wide range of merger timescales, halo masses, and redshifts (0<z<1). We also quantify the utility of pairs in measuring merger bias. While pairs at very small separations will reliably merge, these constitute a small fraction of the total merger population. Thus, pairs do not provide a reliable direct proxy to the total merger population. We do find an intriguing universality in the relation between close pairs and mergers, which in principle could allow for an estimate of the statistical merger rate from the pair fraction within a scaled separation, but including the effects of redshift space distortions strongly degrades this relation. We find similar behavior for galaxy-mass halos, making our results applicable to field galaxy mergers at high redshift. We investigate how the halo merger rate can be statistically described by the halo mass function via the merger kernel (coagulation), finding an interesting environmental dependence of merging: halos within the mass resolution of our simulations merge less efficiently in overdense environments. Specifically, halo pairs with separations less than a few Mpc/h are more likely to merge in underdense environments; at larger separations, pairs are more likely to merge in overdense environments.

We measure the mass function of dark matter halos in a large set of collisionless cosmological simulations of flat LCDM cosmology and investigate its evolution at z<~2. Halos are identified as isolated density peaks, and their masses are measured within a series of radii enclosing specific overdensities. We argue that these spherical overdensity masses are more directly linked to cluster observables than masses measured using the friends-of-friends algorithm (FOF), and are therefore preferable for accurate forecasts of halo abundances. Our simulation set allows us to calibrate the mass function at z=0 for virial masses in the range 10^{11} Msol/h < M < 10^{15} Msol/h, to <~ 5%. We derive fitting functions for the halo mass function in this mass range for a wide range of overdensities, both at z=0 and earlier epochs. In addition to these formulae, which improve on previous approximations by 10-20%, our main finding is that the mass function cannot be represented by a universal fitting function at this level of accuracy. The amplitude of the "universal" function decreases monotonically by ~20-50%, depending on the mass definition, from z=0 to 2.5. We also find evidence for redshift evolution in the overall shape of the mass function.

We present a precise estimate of the bulk virial scaling relation of halos formed via hierarchical clustering in an ensemble of simulated cold dark matter cosmologies. The result is insensitive to cosmological parameters, the presence of a trace, dissipationless gas component, and numerical resolution down to a limit of ~1000 particles. The dark matter velocity dispersion scales with total mass as log(sigma_{DM}(M,z)) = log(1082.9 +- 4.0 \kms) + (0.3361 +- 0.0026) log(h(z)M_{200}/10^{15} Msun), with h(z) the dimensionless Hubble parameter. At fixed mass, the velocity dispersion likelihood is nearly log-normal, with scatter sigma_{ln sigma} = 0.0426 +- 0.015, except for a tail to higher dispersions containing 10% of the population that are merger transients. We combine this relation with the halo mass function in LCDM models, and show that a low normalization condition, S_8 = sigma_8(Omega_m/0.3)^{0.35} \= 0.69, favored by recent WMAP and SDSS analysis requires that galaxy and gas specific energies in rich clusters be 50% larger than that of the underlying dark matter. Such large energetic biases are in conflict with the current generation of direct simulations of cluster formation. A higher normalization, S_8 = 0.80, alleviates this tension and implies that the hot gas fraction within r_{500} is (0.71 +- 0.09) h_{70}^{-3/2} Omega_b/\Omega_m, a value consistent with recent Sunyaev-Zel'dovich observations.

We use measurements of the projected galaxy correlation function w_p and galaxy void statistics to test whether the galaxy content of halos of fixed mass is systematically different in low density environments. We present new measurements of the void probability function (VPF) and underdensity probability function (UPF) from Data Release Four of the Sloan Digital Sky Survey, as well as new measurements of the VPF from the full data release of the Two-Degree Field Galaxy Redshift Survey. We compare these measurements to predictions calculated from models of the Halo Occupation Distribution (HOD) that are constrained to match both w_p and the space density of galaxies. The standard implementation of the HOD assumes that galaxy occupation depends on halo mass only, and is independent of local environment. For luminosity-defined samples, we find that the standard HOD prediction is a good match to the observations, and the data exclude models in which galaxy formation efficiency is reduced in low-density environments. More remarkably, we find that the void statistics of red and blue galaxies (at L ~ 0.4L_*) are perfectly predicted by standard HOD models matched to the correlation function of these samples, ruling out "assembly bias" models in which galaxy color is correlated with large-scale environment at fixed halo mass. We conclude that the luminosity and color of field galaxies are determined predominantly by the mass of the halo in which they reside and have little direct dependence on the environment in which the host halo formed. In broader terms, our results show that the sizes and emptiness of voids found in the distribution of L > 0.2L_* galaxies are in excellent agreement with the predictions of a standard cosmological model with a simple connection between galaxies and dark matter halos. (abridged)

Current and upcoming cosmological observations allow us to probe structures on smaller and smaller scales, entering highly nonlinear regimes. In order to obtain theoretical predictions in these regimes, large cosmological simulations have to be carried out. The promised high accuracy from observations make the simulation task very demanding: the simulations have to be at least as accurate as the observations. This requirement can only be fulfilled by carrying out an extensive code validation program. The first step of such a program is the comparison of different cosmology codes including gravitation interactions only. In this paper we extend a recently carried out code comparison project to include five more simulation codes. We restrict our analysis to a small cosmological volume which allows us to investigate properties of halos. For the matter power spectrum and the mass function, the previous results hold, with the codes agreeing at the 10% level over wide dynamic ranges. We extend our analysis to the comparison of halo profiles and investigate the halo count as a function of local density. We introduce and discuss ParaView as a flexible analysis tool for cosmological simulations, the use of which immensely simplifies the code comparison task.

The clustering properties of dark matter halos are a firm prediction of modern theories of structure formation. We use two large volume, high-resolution N-body simulations to study how the correlation function of massive dark matter halos depends upon their mass and formation history. We find that halos with the lowest concentrations are presently more clustered than those of higher concentration, the size of the effect increasing with halo mass; this agrees with trends found in studies of lower mass halos. The clustering dependence on other characterizations of the full mass accretion history appears weaker than the effect with concentration. Using the integrated correlation function, marked correlation functions, and a power-law fit to the correlation function, we find evidence that halos which have recently undergone a major merger or a large mass gain have slightly enhanced clustering relative to a randomly chosen population with the same mass distribution.

We identify galaxy groups and clusters in volume-limited samples of the SDSS redshift survey, using a redshift-space friends-of-friends algorithm. We optimize the friends-of-friends linking lengths to recover galaxy systems that occupy the same dark matter halos, using a set of mock catalogs created by populating halos of N-body simulations with galaxies. Extensive tests with these mock catalogs show that no combination of perpendicular and line-of-sight linking lengths is able to yield groups and clusters that simultaneously recover the true halo multiplicity function, projected size distribution, and velocity dispersion. We adopt a linking length combination that yields, for galaxy groups with ten or more members: a group multiplicity function that is unbiased with respect to the true halo multiplicity function; an unbiased median relation between the multiplicities of groups and their associated halos; a spurious group fraction of less than ~1%; a halo completeness of more than ~97%; the correct projected size distribution as a function of multiplicity; and a velocity dispersion distribution that is ~20% too low at all multiplicities. These results hold over a range of mock catalogs that use different input recipes of populating halos with galaxies. We apply our group-finding algorithm to the SDSS data and obtain three group and cluster catalogs for three volume-limited samples that cover 3495.1 square degrees on the sky. We correct for incompleteness caused by fiber collisions and survey edges, and obtain measurements of the group multiplicity function, with errors calculated from realistic mock catalogs. These multiplicity function measurements provide a key constraint on the relation between galaxy populations and dark matter halos.

The coupling of photons and baryons by Thomson scattering in the early universe imprints features in both the Cosmic Microwave Background (CMB) and matter power spectra. The former have been used to constrain a host of cosmological parameters, the latter have the potential to strongly constrain the expansion history of the universe and dark energy. Key to this program is the means to localize the primordial features in observations of galaxy spectra which necessarily involve galaxy bias, non-linear evolution and redshift space distortions. We present calculations, based on mock catalogs produced from high-resolution N-body simulations, which show the range of behaviors we might expect of galaxies in the real universe. We investigate physically motivated fitting forms which include the effects of non-linearity, galaxy bias and redshift space distortions and discuss methods for analysis of upcoming data. In agreement with earlier work, we find that a survey of several Gpc^3 would constrain the sound horizon at z~1 to about 1%.

(Abridged) We make predictions for the pairwise velocity dispersion (PVD) of galaxies with models that are constrained to match the projected correlation function and luminosity function of galaxies in the Two-Degree Field Galaxy Redshift Survey (2dFGRS). We use these data to constrain the halo occupation distribution (HOD), then calculate the PVD by populating the halos of a high resolution N-body simulation. We examine the luminosity and scale dependence of the predicted PVD. At r<1 Mpc/h and r>3 Mpc/h, we find that the PVD decreases with increasing galaxy luminosity. This result is mostly driven by the fraction of satellite galaxies f_sat, which is well-constrained by the correlation function. We find f_sat~25% for galaxies fainter than L_star, while for brighter galaxies the satellite fraction rapidly declines, creating the decrease in the PVD with luminosity. At r=1 Mpc/h, the PVD has no dependence on luminosity because satellite galaxies dominate the statistics for all objects. Recent measurements of the PVD in Fourier space using the "dispersion model" have reported a strong decline in PVD with increasing luminosity at k=1 h/Mpc. We test this method with our HOD models, finding that there is no consistent comparison between the PVD at a given k and the true dispersion at a given value of r. This results in a luminosity dependence in k-space that is stronger than in configuration space. The luminosity dependence of the HOD results in Fourier space are consistent with those measured at k=1 h/Mpc; thus the recent measurements of the PVD are fully explainable in the context of halo occupation models. The normalization of the PVD is lower than predicted by our fiducial model, and reproducing it requires a lower value of Omega_m (~0.2 instead of 0.3), a lower value of sigma_8 (~0.7 instead of 0.9), or strong velocity bias.

(Abridged) We investigate the power of void statistics to constrain galaxy bias and the amplitude of dark matter fluctuations. We use the halo occupation distribution (HOD) framework to describe the relation between galaxies and dark matter. After choosing HOD parameters that reproduce the mean space density n_gal and projected correlation function w_p measured for galaxy samples with M_r<-19 and M_r<-21 from the Sloan Digital Sky Survey (SDSS), we predict the void probability function (VPF) and underdensity probability function (UPF) of these samples by populating the halos of a large, high-resolution N-body simulation. If we make the conventional assumption that the HOD is independent of large scale environment at fixed halo mass, then models constrained to match n_gal and w_p predict nearly identical void statistics, independent of the scatter between halo mass and central galaxy luminosity or uncertainties in HOD parameters. Models with sigma_8=0.7 and sigma_8=0.9 also predict very similar void statistics. However, the VPF and UPF are sensitive to environmental variations of the HOD in a regime where these variations have little impact on w_p. For example, doubling the minimum host halo mass in regions with large scale (5 Mpc/h) density contrast delta<-0.65 has a readily detectable impact on void probabilities of M_r<-19 galaxies, and a similar change for delta<-0.2 alters the void probabilities of M_r<-21 galaxies at a detectable level. The VPF and UPF provide complementary information about the onset and magnitude of density- dependence in the HOD. By detecting or ruling out HOD changes in low density regions, void statistics can reduce systematic uncertainties in the cosmological constraints derived from HOD modeling, and, more importantly, reveal connections between halo formation history and galaxy properties.

We provide a description of the SNSPH code--a parallel 3-dimensional radiation hydrodynamics code implementing treecode gravity, smooth particle hydrodynamics, and flux-limited diffusion transport schemes. We provide descriptions of the physics and parallelization techniques for this code. We present performance results on a suite of code tests (both standard and new), showing the versatility of such a code, but focusing on what we believe are important aspects of modeling core-collapse supernovae.

We present the small-scale (0.01<r<8 h^{-1} Mpc) projected correlation function w_p(r_p) and real space correlation function xi(r) of 24520 luminous early-type galaxies from the Sloan Digital Sky Survey Luminous Red Galaxy (LRG) sample (0.16<z<0.36). ``Fiber collision'' incompleteness of the SDSS spectroscopic sample at scales smaller than 55 arcseconds prevents measurements of the correlation function for LRGs on scales smaller than ~0.3 Mpc by the usual methods. In this work, we cross-correlate the spectroscopic sample with the imaging sample, with a weighting scheme to account for the collisions, extensively tested against mock catalogs. We correct for photometric biases in the SDSS imaging of close galaxy pairs. We find that the correlation function xi(r) is surprisingly close to a r^{-2} power law over more than 4 orders of magnitude in separation r. This result is too steep at small scales to be explained in current versions of the halo model for galaxy clustering. We infer an LRG-LRG merger rate of less than or similar to 0.6 Gyr^{-1} Gpc^{-3} for this sample. This result suggests that the LRG-LRG mergers are not the main mode of mass growth for LRGs at z<0.36.

The predicted mass function of dark matter halos is essential in connecting observed galaxy cluster counts and models of galaxy clustering to the properties of the primordial density field. We determine the mass function in the concordance $\Lambda$CDM cosmology, as well as its uncertainty, using sixteen $1024^3$-particle nested-volume dark-matter simulations, spanning a mass range of over five orders of magnitude. Using the nested volumes and single-halo tests, we find and correct for a systematic error in the friends-of-friends halo-finding algorithm. We find a fitting form and full error covariance for the mass function that successfully describes the simulations' mass function and is well-behaved outside the simulations' resolutions. Estimated forecasts of uncertainty in cosmological parameters from future cluster count surveys have negligible contribution from remaining statistical uncertainties in the central cosmology multiplicity function. There exists a potentially non-negligible cosmological dependence (non-universality) of the halo multiplicity function.

The gravitationally-driven evolution of cold dark matter dominates the formation of structure in the Universe over a wide range of length scales. While the longest scales can be treated by perturbation theory, a fully quantitative understanding of nonlinear effects requires the application of large-scale particle simulation methods. Additionally, precision predictions for next-generation observations, such as weak gravitational lensing, can only be obtained from numerical simulations. In this paper, we compare results from several N-body codes using test problems and a diverse set of diagnostics, focusing on a medium resolution regime appropriate for studying many observationally relevant aspects of structure formation. Our conclusions are that -- despite the use of different algorithms and error-control methodologies -- overall, the codes yield consistent results. The agreement over a wide range of scales for the cosmological tests is test-dependent. In the best cases, it is at the 5% level or better, however, for other cases it can be significantly larger than 10%. These include the halo mass function at low masses and the mass power spectrum at small scales. While there exist explanations for most of the discrepancies, our results point to the need for significant improvement in N-body errors and their understanding to match the precision of near-future observations. The simulation results, including halo catalogs, and initial conditions used, are publicly available.

On large scales galaxies and their halos are usually assumed to trace the dark matter with a constant bias and dark matter is assumed to trace the linear density field. We test these assumption using several large N-body simulations with 384^3-1024^3 particles and box sizes between 100-1000h/Mpc, which can both resolve the small galactic size halos and sample the large scale fluctuations. We explore the average halo bias relation as a function of halo mass and show that existing fitting formulae overestimate the halo bias by up to 20% in the regime just below the nonlinear mass. We propose a new expression that fits our simulations well. We find that the halo bias is nearly constant, b~0.65-0.7, for masses below one tenth of the nonlinear mass. We explore next the relation between the initial and final dark matter in individual Fourier modes and show that there are significant fluctuations in their ratio, ranging from 10% rms at k~0.03h/Mpc to 50% rms at k~0.1h/Mpc. We argue that these large fluctuations are caused by perturbative effects beyond the linear theory, which are dominated by long wavelength modes with large random fluctuations. Similar or larger fluctuations exist between halos and dark matter and between halos of different mass. While these fluctuations are small compared to the sampling variance, they are significant for attempts to determine the bias by relating directly the maps of galaxies and dark matter or the maps of different galaxy populations, which would otherwise be immune to sampling variance.

The use of parallel computers and increasingly sophisticated software has allowed us to perform a large suite of N-body simulations using from $10^8$,to $10^9$ particles. We will report on our recent convergence tests of the halo mass function, N-point correlation functions, power spectrum and pairwise velocity from very large high resolution treecode N-body simulations. Rather than basing results on just one or two large simulations, now one can investigate the role of different numerical and physical effects on the statistics used to characterize the mass distribution of the Universe.

Recent observations with the Chandra X-ray Observatory have provided us with the capability to discriminate point sources, such as the supermassive black hole Sgr A*, from the diffuse emission within the inner 10'' of the Galaxy. The hot plasma producing the diffuse X-radiation, estimated at ~7.6*10^{31} ergs s^{-1} arcsec^{-2} in the 2-10 keV band, has a RMS electron density ~26 cm^{-3} and a temperature kT ~= 1.3 keV, with a total inferred mass of ~0.1 solar masses. At least some of this gas must be injected into the ISM via stellar winds. In the most recent census, about 25 bright, young stars have been identified as the dominant sources of the overall mass efflux from the Galactic center. In this paper, we use detailed 3-dimensional SPH simulations to study the wind-wind interactions occurring in the inner 3 parsecs of the Galaxy, with a goal of understanding what fraction, if any, of the diffuse X-ray flux measured by Chandra results from the ensuing shock heating of the ambient medium. We conclude that this process alone can account for the entire X-ray flux observed by Chandra in the inner 10'' of the Galaxy. Understanding the X-ray morphology of the environment surrounding Sgr A* will ultimately provide us with a greater precision in modeling the accretion of gas onto this object, which appears to be relatively underluminous compared to its brethren in the nuclei of other galaxies.

A scenario is put forth for the formation of supermassive black holes at the centers of galaxies. It depends upon the formation of a Mestel disk with a flat rotation curve, where the mass is proportional to r and thickness to 1/r. Such disks should form from the collapse of uniformly rotating, isolated, gaseous clouds, either proto-galactic, galaxy-mass damped Lyman alpha clouds or the gas that survives galaxy mergers. We propose that in any case the disk will be unstable to the Rossby vortex instability (RVI). This instability grows from any large, steep pressure gradient in an optically thick disk. Such pressure gradients either occur adjacent to compact objects or could be triggered by individual supernovae in and around the disk. Upon excitation, the RVI transports angular momentum far more efficiently than the turbulence of a viscosity-based Shakura-Sunyaev disk. The critical column density condition in a L* Mestel disk predicts that the black hole mass is ~ 10^8 Msolar and is also proportional to the 1/4 power of the velocity dispersion, in accordance with observations.

In this paper, we present the results of 3-dimensional collapse simulations of rotating stars for a range of stellar progenitors. We find that for the fastest spinning stars, rotation does indeed modify the convection above the proto-neutron star, but it is not fast enough to cause core fragmentation. Similarly, although strong magnetic fields can be produced once the proto-neutron star cools and contracts, the proto-neutron star is not spinning fast enough to generate strong magnetic fields quickly after collapse and, for our simulations, magnetic fields will not dominate the supernova explosion mechanism. Even so, the resulting pulsars for our fastest rotating models may emit enough energy to dominate the total explosion energy of the supernova. However, more recent stellar models predict rotation rates that are much too slow to affect the explosion, but these models are not sophisticated enough to determine whether the most recent, or past, stellar rotation rates are most likely. Thus, we must rely upon observational constraints to determine the true rotation rates of stellar cores just before collapse. We conclude with a discussion of the possible constraints on stellar rotation which we can derive from core-collapse supernovae.

We present 3-dimensional SPH simulations of supernova explosions from 100 seconds to 1 year after core-bounce. By extending our modelling efforts to a 3-dimensional hydrodynamics treatment, we are able to investigate the effects of explosion asymmetries on mixing and gamma-ray line emergence in supernovae. A series of initial explosion conditions are implemented, including jet-like and equatorial asymmetries of varying degree. For comparison, symmetric explosion models are also calculated. A series of time slices from the explosion evolution are further analyzed using a 3-dimensional Monte Carlo gamma-ray transport code. The emergent hard X- and gamma-ray spectra are calculated as a function of both viewing angle and time, including trends in the gamma-ray line profiles. We find significant differences in the velocity distribution of radioactive nickel between the symmetric and asymmetric explosion models. The effects of this spatial distribution change are reflected in the overall high energy spectrum, as well as in the individual gamma-ray line profiles.

The new generation of gravitational wave (GW) detectors have the potential to open a novel window onto the violent dynamics of core collapse. Although it is certain that core collapse events generate gravitational radiation, understanding the characteristics of the radiation -- whether it can be measured with these detectors, and the best way to go about doing so -- is a challenging problem. In this chapter we review the promise of GWs as observational probes, including a discussion of the current state of GW detectors, and discuss the status of work to understand the waves generated by stellar core collapse.

We present the first complete 3-dimensional simulations of the core-collapse of a massive star from the onset of collapse to the resultant supernova explosion. We compare the structure of the convective instabilities that occur in 3-dimensional models with those of past 2-dimensional simulations. Although the convective instabilities are clearly 3-dimensional in nature, we find that both the size-scale of the flows and the net enhancement to neutrino heating does not differ greatly between 2- and 3-dimensional models. The explosion energy, explosion timescale, and remnant mass does not differ by more than 10% between 2- and 3-dimensional simulations.

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 propose a method to determine the cosmic mass density Omega from redshift-space distortions induced by large-scale flows in the presence of nonlinear clustering. Nonlinear structures in redshift space such as fingers of God can contaminate distortions from linear flows on scales as large as several times the small-scale pairwise velocity dispersion sigma_v. Following Peacock & Dodds (1994), we work in the Fourier domain and propose a model to describe the anisotropy in the redshift-space power spectrum; tests with high-resolution numerical data demonstrate that the model is robust for both mass and biased galaxy halos on translinear scales and above. On the basis of this model, we propose an estimator of the linear growth parameter beta = Omega^0.6/b, where b measures bias, derived from sampling functions which are tuned to eliminate distortions from nonlinear clustering. The measure is tested on the numerical data and found to recover the true value of beta to within ~10%. An analysis of the IRAS 1.2Jy galaxies yields beta = 0.8+0.4/-0.3 at a scale of 1,000 km/s which is close to optimal given the shot noise and the finite survey volume. This measurement is consistent with dynamical estimates of beta derived from both real-space and redshift-space information. The importance of the method presented here is that nonlinear clustering effects are removed to enable linear correlation anisotropy measurements on scales approaching the translinear regime. We discuss implications for analyses of forthcoming optical redshift surveys in which the dispersion is more than a factor of two greater than in the IRAS data.