The details of the stellar mass assembly of brightest cluster galaxies (BCGs) remain an unresolved problem in galaxy formation. We have developed a novel approach that allows us to construct a sample of clusters that form an evolutionary sequence, and have applied it to the Spitzer IRAC Shallow Cluster Survey (ISCS) to examine the evolution of BCGs in progenitors of present-day clusters with mass of (2.5-4.5)x10^{14}Msun. We follow the cluster mass growth history extracted from a high resolution cosmological simulation, and then use an empirical method that infers the cluster mass based on the ranking of cluster luminosity to select high-z clusters of appropriate mass from ISCS to be progenitors of the given set of z=0 clusters. We find that, between z=1.5 and 0.5, the BCGs have grown in stellar mass by a factor of 2.3, which is well-matched by the predictions from a state-of-the-art semi-analytic model. Below z=0.5 we see hints of differences in behavior between the model and observation.

We present a new method for incorporating nonthermal pressure from bulk motions of gas into an analytic model of the intracluster medium in clusters of galaxies, which is based on a polytropic equation of state and hydrostatic equilibrium inside gravitational potential wells drawn from cosmological dark matter simulations. The pressure is allowed to have thermal and nonthermal components with different radial distributions; the overall level of nonthermal support is based on the dynamical state of the halo, such that it is lower in more relaxed clusters. This level is normalized by comparison to pressure profiles derived from X-ray observations, and to a high resolution hydrodynamical simulation. The nonthermal pressure fraction measured at r_500 is typically in the range 10-20%, increasing with cluster mass and with redshift. The resulting model cluster properties are in accord with Sunyaev-Zel'dovich (SZ) effect observations of clusters. Inclusion of nonthermal pressure reduces the expected angular power spectrum of SZ fluctuations in the microwave sky by 24%.

Observations of clusters of galaxies suggest that they contain significantly fewer baryons (gas plus stars) than the cosmic baryon fraction. This `missing baryon' puzzle is especially surprising for the most massive clusters which are expected to be representative of the cosmic matter content of the universe (baryons and dark matter). Here we show that the baryons may not actually be missing from clusters, but rather are extended to larger radii than typically observed. The baryon deficiency is typically observed in the central regions of clusters (~0.5 the virial radius). However, the observed gas-density profile is significantly shallower than the mass-density profile, implying that the gas is more extended than the mass and that the gas fraction increases with radius. We use the observed density profiles of gas and mass in clusters to extrapolate the measured baryon fraction as a function of radius and as a function of cluster mass. We find that the baryon fraction reaches the cosmic value near the virial radius for all groups and clusters above 5e13 solar masses. This suggests that the baryons are not missing, they are simply located in cluster outskirts. Heating processes (shock-heating of the intracluster gas, plus supernovae and AGN feedback) that cause the gas to expand are likely explanations for these results. Upcoming observations should be able to detect these baryons.

We present templates for the Sunyaev-Zel'dovich (SZ) angular power spectrum based on four models for the nonlinear gas distribution. The frequency-dependent SZ temperature fluctuations, with thermal (TSZ) and kinetic (KSZ) contributions, are calculated by tracing through a dark matter simulation, processed to include gas in dark matter halos and in the filamentary intergalactic medium. Different halo gas models are compared to study how star formation, energetic feedback, and nonthermal pressure support influence the angular power spectrum. The standard model has been calibrated to reproduce the stellar and gas fractions and X-ray scaling relations measured from low redshift clusters and groups. The other models illustrate the current theoretical and empirical uncertainties relating to properties of the intracluster medium. Relative to the standard model, their angular power spectra differ by approximately 50% (TSZ), 20% (KSZ), and 40% (SZ at 148 GHz) for l=3000, sigma_8=0.8, and homogeneous reionization at z=10. The angular power spectrum decreases in amplitude as gas mass and binding energy is removed through star formation, and as gas is pushed out to larger radii by energetic feedback. With nonthermal pressure support, less pressure is required to maintain hydrostatic equilibrium, thus reducing the thermal contribution to the SZ power. We also calculate the SZ templates as a function of sigma_8 and quantify this dependence. Assuming C_l is proprotional to (sigma_8/0.8)^alpha, the effective scaling index ranges from 7<alpha_tsz<9, 4.5<alpha_ksz<5.5, and 6.5<alpha_sz(148 GHz)<8 at l=3000 for 0.6<sigma_8<1. The template spectra are publicly available and can be used when fitting for the SZ contribution to the cosmic microwave background on arcminute scales.

We create realistic, full-sky, half-arcminute resolution simulations of the microwave sky matched to the most recent astrophysical observations. The primary purpose of these simulations is to test the data reduction pipeline for the Atacama Cosmology Telescope (ACT) experiment; however, we have widened the frequency coverage beyond the ACT bands to make these simulations applicable to other microwave background experiments. Some of the novel features of these simulations are that the radio and infrared galaxy populations are correlated with the galaxy cluster populations, the CMB is lensed by the dark matter structure in the simulation via a ray-tracing code, the contribution to the thermal and kinetic Sunyaev-Zel'dovich (SZ) signals from galaxy clusters, groups, and the IGM has been included, and the gas prescription to model the SZ signals matches the most recent X-ray observations. Regarding the contamination of cluster SZ flux by radio galaxies, we find for 148 GHz (90 GHz) only 3% (4%) of halos have their SZ decrements contaminated at a level of 20% or more. We find the contamination levels higher for infrared galaxies. However, at 90 GHz, less than 20% of clusters with M_{200} > 2.5 x 10^{14} Msun and z<1.2 have their SZ decrements filled in at a level of 20% or more. At 148 GHz, less than 20% of clusters with M_{200} > 2.5 x 10^{14} Msun and z<0.8 have their SZ decrements filled in at a level of 50% or larger. Our models also suggest that a population of very high flux infrared galaxies, which are likely lensed sources, contribute most to the SZ contamination of very massive clusters at 90 and 148 GHz. These simulations are publicly available and should serve as a useful tool for microwave surveys to cross-check SZ cluster detection, power spectrum, and cross-correlation analyses.

The state of the hot gas in clusters of galaxies is investigated with a set of model clusters, created by assuming a polytropic equation of state (Gamma=1.2) and hydrostatic equilibrium inside gravitational potential wells drawn from a dark matter simulation. Star formation, energy input, and nonthermal pressure support are included. To match the gas fractions seen in non-radiative hydrodynamical simulations, roughly 5% of the binding energy of the dark matter must be transferred to the gas during cluster formation; the presence of nonthermal pressure support increases this value. In order to match X-ray observations, scale-free behavior must be broken. This can be due to either variation of the efficiency of star formation with cluster mass M_500, or the input of additional energy proportional to the formed stellar mass M_F. These two processes have similar effects on X-ray scalings. If 9% of the gas is converted into stars, independent of cluster mass, then feedback energy input of 1.2e-5*M_Fc^2 (or ~1.0 keV per particle) is required to match observed clusters. Alternatively, if the stellar mass fraction varies as M_500^-0.26 then a lower feedback of 4e-6*M_Fc^2 is needed, and if the stellar fraction varies as steeply as M_500^-0.49 then no additional feedback is necessary. The model clusters reproduce the observed trends of gas temperature and gas mass fraction with cluster mass, as well as observed entropy and pressure profiles; thus they provide a calibrated basis with which to interpret upcoming SZ surveys. One consequence of the increased gas energy is that the baryon fraction inside the virial radius is less than roughly 90% of the cosmic mean, even for the most massive clusters.

The complex relationship between the galaxy density field and the underlying matter field limits our ability to extract cosmological constraints from galaxy redshift surveys. Our approach is to use halos rather than galaxies to trace the underlying mass distribution. We identify Fingers-of-God (FOGs) and replace multiple galaxies in each FOG with a single halo object. This removes the nonlinear contributions of satellite galaxies, the one-halo term. We test our method on a large set of high-fidelity mock SDSS Luminous Red Galaxy (LRG) catalogs. We find that the aggressive FOG compression algorithm adopted in the LRG P(k) analysis of Tegmark et al. (2006) leads to a ~10% correction to the underlying matter power spectrum at k = 0.1 h/Mpc and ~40% correction at k=0.2 h/Mpc, thereby compromising the cosmological constraints. In contrast, the power spectrum of our reconstructed halo density field deviates from the underlying matter power spectrum by less than 1% for k less than 0.1 h/Mpc and less than 4% for k less than 0.2 h/Mpc. The reconstructed halo density field also removes the bias in the measurement of the redshift space distortion parameter beta induced by the FOG smearing of the linear redshift space distortions.

Large scale structure deflects cosmic microwave background (CMB) photons. Since large angular scales in the large scale structure contribute significantly to the gravitational lensing effect, a realistic simulation of CMB lensing requires a sufficiently large sky area. We describe simulations that include these effects, and present both effective and multiple plane ray-tracing versions of the algorithm, which employs spherical harmonic space and does not use the flat sky approximation. We simulate lensed CMB maps with an angular resolution of ~0.9 arcmin. The angular power spectrum of the simulated sky agrees well with analytical predictions. Maps generated in this manner are a useful tool for the analysis and interpretation of upcoming CMB experiments such as PLANCK and ACT.

We present a detailed analysis of the intrinsic scatter in the integrated SZ effect - cluster mass (Y-M) relation, using semi-analytic and simulated cluster samples. Specifically, we investigate the impact on the Y-M relation of energy feedback, variations in the host halo concentration and substructure populations, and projection effects due to unresolved clusters along the line of sight (the SZ background). Furthermore, we investigate at what radius (or overdensity) one should measure the integrated SZE and define cluster mass so as to achieve the tightest possible scaling. We find that the measure of Y with the least scatter is always obtained within a smaller radius than that at which the mass is defined; e.g. for M_{200} (M_{500}) the scatter is least for Y_{500} (Y_{1100}). The inclusion of energy feedback in the gas model significantly increases the intrinsic scatter in the Y-M relation due to larger variations in the gas mass fraction compared to models without feedback. We also find that variations in halo concentration for clusters of a given mass may partly explain why the integrated SZE provides a better mass proxy than the central decrement. Substructure is found to account for approximately 20% of the observed scatter in the Y-M relation. Above M_{200} = 2x10^{14} h^{-1} msun, the SZ background does not significantly effect cluster mass measurements; below this mass, variations in the background signal reduce the optimal angular radius within which one should measure Y to achieve the tightest scaling with M_{200}.

The statistics of giant arcs and large separation lensed quasars provide powerful constraints for the parameters of the underlying cosmological model. So far, most investigations have been carried out using pure dark matter simulations. Here we present a recipe for including the effects of baryon cooling (i.e. large galaxy formation) in dark matter N-body simulations that is consistent with observations of massive galaxies. Then we quantitatively compare lensing with and without applying this baryon correction to the pure dark matter case. Including the baryon correction significantly increases the frequency of giant arcs and lensed quasars, particularly on scales of 10 arcsec and smaller: the overall frequency of multiple images increases by about 25% for source redshifts between z_s = 1.5 and 7.5 and splittings larger than about 3 arcsec. The baryon rearrangement also slightly increases the fraction of quadruple images over doubles.

Using the full, three-dimensional potential of galaxy cluster halos (drawn from an N-body simulation of the current, most favored cosmology), the distribution of the X-ray emitting gas is found by assuming a polytropic equation of state and hydrostatic equilibrium, with constraints from conservation of energy and pressure balance at the cluster boundary. The resulting properties of the gas for these simulated redshift zero clusters (the temperature distribution, mass-temperature and luminosity-temperature relations, and the gas fraction) are compared with observations in the X-ray of nearby clusters. The observed properties are reproduced only under the assumption that substantial energy injection from non-gravitational sources has occurred. Our model does not specify the source, but star formation and AGN may be capable of providing this energy, which amounts to 3 to 5 x10^{-5} of the rest mass in stars (assuming ten percent of the gas initially in the cluster forms stars). With the method described here it is possible to generate realistic X-ray and Sunyaev-Zel'dovich cluster maps and catalogs from N-body simulations, with the distributions of internal halo properties (and their trends with mass, location, and time) taken into account.

We study the ability of three-frequency, arcminute-resolution microwave measurements to detect galaxy clusters via their Sunyaev-Zel'dovich (SZ) distortion of the microwave background. For this purpose, we have constructed large-area simulations of the microwave sky, and we have made them publicly available to further investigations into optimal data reduction techniques for upcoming SZ cluster surveys. In these sky simulations, galaxy clusters are modeled using N-body simulated dark matter halos plus a gas prescription for the intracluster medium that allows the small scale cluster physics such as star formation and feedback to be realistically incorporated. We also model the primary microwave background, radio and infrared point sources, galactic dust emission, and the SZ flux including kinetic and relativistic contributions. We make use of these simulations to study the cluster scaling relation between integrated SZ flux and cluster mass and find our clusters fit a power-law well, with a power-law index that is steeper than that for self-similar cluster models. Some evolution of the power-law index and normalization with redshift is also observed. These simulations are also used to study cluster detection for the Atacama Cosmology Telescope (ACT). Using a multi-frequency Wiener filter to separate clusters from other microwave components, we find that ACT alone can recover a cluster sample that is ~90% complete above 3*10^{14} M_sun and ~85% pure.

We present a new definition of subhalos in dissipationless dark matter N-body simulations, based on the coherent identification of their dynamically bound constituents. Whereas previous methods of determining the energetically bound components of a subhalo ignored the contribution of all the remaining particles in the halo (those not geometrically or dynamically associated with the subhalo), our method allows for all the forces, both internal and external, exerted on the subhalo. We demonstrate, using the output of a simulation at different timesteps, that our new method is more accurate at identifying the bound mass of a subhalo. We then compare our new method to previously adopted means of identifying subhalos by applying each to a sample of 1838 virialized halos extracted from a high resolution cosmological simulation. We find that the subhalo distributions are similar in each case, and that the increase in the binding energy of a subhalo from including all the particles located within it is almost entirely balanced by the losses due to the external forces; the net increase in the mass fraction of subhalos is roughly 10%, and the extra substructures tending to reside in the inner parts of the system. Finally, we compare the subhalo populations of halos to the sub-subhalo populations of subhalos, finding the two distributions to be similar. This is a new and interesting result, suggesting a self-similarity in the hierarchy substructures within cluster mass halos.

We determine the shape, multiplicity, size, and radial structure of superclusters in the LambdaCDM concordance cosmology from z = 0 to z = 2. Superclusters are defined as clusters of clusters in our large-scale cosmological simulation. We find that superclusters are triaxial in shape; many have flattened since early times to become nearly two-dimensional structures at present, with a small fraction of filamentary systems. The size and multiplicity functions are presented at different redshifts. Supercluster sizes extend to scales of ~ 100 - 200 Mpc/h. The supercluster multiplicity (richness) increases linearly with supercluster size. The density profile in superclusters is approximately isothermal (~ R^{-2}) and steepens on larger scales. These results can be used as a new test of the current cosmology when compared with upcoming observations of large-scale surveys.

We have identified over 2000 well resolved cluster halos, and also their associated bound subhalos, from the output of a 1024^3 particle cosmological N-body simulation (of box size 320 h^-1 Mpc and softening length 3.2 h^-1 kpc). This has allowed us to measure halo quantities in a statistically meaningful way, and for the first time analyse their distribution for a large and well resolved sample. We characterize each halo in terms of its morphology, concentration, spin, circular velocity and the fraction of their mass in substructure. We also identify those halos that have not yet reached a state of dynamical equilibrium using the virial theorem with an additional correction to account for the surface pressure at the boundary. These amount to 3.4% of our initial sample. For the virialized halos, we find a median of 5.6% of halo mass is contained within substructure, with the distribution ranging between no identified subhalos to 65%. The fraction of mass in substructure increases with halo mass with logarithmic slope of 0.44 +- 0.06. Halos tend to have a prolate morphology, becoming more so with increasing mass. Subhalos have a greater orbital angular momentum per unit mass than their host halo. Furthermore, their orbital angular momentum is typically well aligned with that of their host. Overall, we find that dimensionless properties of dark matter halos do depend on their mass, thereby demonstrating a lack of self-similarity.

We investigate the dependence of ellipticities of clusters of galaxies on cosmological parameters using large-scale cosmological simulations. We determine cluster ellipticities out to redshift unity for LCDM models with different mean densities $\Omega_m$ and amplitudes of mass fluctuation $\sigma_{8,0}$. The mean ellipticity increases monotonically with redshift for all models. Larger values of $\sigma_{8,0}$, i.e., earlier cluster formation time, produce lower ellipticities. The dependence of ellipticity on $\Omega_m$ is relatively weak in the range $0.2 \leq \Omega_m \leq 0.5$ for high mass clusters. The mean ellipticity $\bar{e}(z)$ decreases linearly with the amplitude of fluctuations at the cluster redshift $z$, nearly independent of $\Omega_m$; on average, older clusters are more relaxed and are thus less elliptical. The distribution of ellipticities about the mean is approximated by a Gaussian, allowing a simple characterization of the evolution of ellipticity with redshift as a function of cosmological parameters. At $z=0$, the mean ellipticity of high mass clusters is approximated by $\bar{e}(z=0) = 0.248-0.069 \sigma_{8,0} + 0.013 \Omega_{m,0}$. This relation opens up the possibility that, when compared with future observations of large cluster samples, the mean cluster ellipticity and its evolution could be used as a new, independent tool to constrain cosmological parameters, especially the amplitude of mass fluctuations, $\sigma_{8,0}$.

Starting with the well-known NFW dark matter halo distribution, we construct a simple polytropic model for the intracluster medium which is in good agreement with high resolution numerical hydrodynamical simulations, apply this model to a very large scale concordance dark matter simulation, and compare the resulting global properties with recent observations of X-ray clusters, including the mass-temperature and luminosity-temperature relations. We make allowances for a non-negligible surface pressure, removal of low entropy (short cooling time) gas, energy injection due to feedback, and for a relativistic (non-thermal) pressure component. A polytropic index n=5 (Gamma=1.2) provides a good approximation to the internal gas structure of massive clusters (except in the very central regions where cooling becomes important), and allows one to recover the observed M_500-T, L_x-T, and T/n_e^{2/3} ~ T^0.65 relations. Using these concepts and generalizing this method so that it can be applied to fully three-dimensional N-body simulations, one can predict the global X-ray and SZE trends for any specified cosmological model. We find a good fit to observations when assuming that twelve percent of the initial baryonic mass condenses into stars, the fraction of rest mass of this condensed component transferred back to the remaining gas (feedback) is 3.9E-5, and the fraction of total pressure from a nonthermal component is near ten percent.

We present an algorithm which is designed to allow the efficient identification and preliminary dynamical analysis of thousands of structures and substructures in large N-body simulations. First we utilise a refined density gradient system (based on DENMAX) to identify the structures, and then apply an iterative approximate method to identify unbound particles, allowing fast calculation of bound substructures. After producing a catalog of separate energetically bound substructures we check to see which of these are energetically bound to adjacent substructures. For such bound complex subhalos, we combine components and check if additional free particles are also bound to the union, repeating the process iteratively until no further changes are found. Thus our subhalos can contain more than one density maximum, but the scheme is stable: starting with a small smoothing length initially produces small structures which must be combined later, and starting with a large smoothing length produces large structures within which sub-substructure is found. We apply this algorithm to three simulations. Two which are using the TPM algorithm by Bode et al. (2000) and one on a simulated halo by Diemand (2004). For all these halos we find about 5-8% of the mass in substructures.

We compute the expected number of quasars multiply imaged by cluster size dark halos for current wide field quasar surveys by carrying out a large ensemble of ray tracing simulations through clusters from a cosmological N-body simulation of the LCDM cosmology. Our calculation predicts ~ 4 quasar lenses with splittings theta > 10" in the SDSS spectroscopic quasar sample, consistent with the recent discovery of the wide separation lens SDSSJ1004+4112 which has theta=14.6". The SDSS faint photometric quasar survey will contain ~12 multiply imaged quasars with splittings theta > 10". Of these, ~ 2 will be lenses with separations theta > 30", and ~ 2 will be at high redshift (z ~ 4).

We present a detailed investigation into which properties of CDM halos make them effective strong gravitational lenses. Strong lensing cross sections of 878 clusters from an N-body simulation are measured by ray tracing through 13,594 unique projections. We measure concentrations, axis ratios, orientations, and the amount of substructure of each cluster, and compare the lensing weighted distribution of each quantity to that of the cluster population as a whole. The concentrations of lensing clusters are on average 34% larger than the typical cluster in the Universe. Despite this bias, the anomalously high concentrations (c >14) recently measured by several groups, appear to be inconsistent with the concentration distribution in our simulations, which predict < 2% of lensing clusters should have concentrations this high. No correlation is found between lensing cross section and the amount of substructure. We introduce several types of simplified dark matter halos, and use them to isolate which properties of CDM clusters make them effective lenses. Projections of halo substructure onto small radii and the large scale mass distribution of clusters do not significantly influence cross sections. The abundance of giant arcs is primarily determined by the mass distribution within an average overdensity of ~ 10,000. A multiple lens plane ray tracing algorithm is used to show that projections of large scale structure increase the giant arc abundance by a modest amount <7%. We revisit the question of whether there is an excess of giant arcs behind high redshift clusters in the RCS survey and find that the number of high redshift (z > 0.6) lenses is in good agreement with LCDM, although our simulations predict more low redshift (z < 0.6) lenses than were observed. (abridged)

The evolution of the cluster mass function and the cluster correlation function from z = 0 to z = 3 are determined using 10^6 clusters obtained from high-resolution simulations of the current best-fit LCDM cosmology (\Omega_m = 0.27, \sigma_8 = 0.84, h = 0.7). The results provide predictions for comparisons with future observations of high redshift clusters. A comparison of the predicted mass function of low redshift clusters with observations from early Sloan Digital Sky Survey data, and the predicted abundance of massive distant clusters with observational results, favor a slightly larger amplitude of mass fluctuations (\sigma_8 = 0.9) and lower density parameter (\Omega_m = 0.2); these values are consistent within 1-\sigma with the current observational and model uncertainties. The cluster correlation function strength increases with redshift for a given mass limit; the clusters were more strongly correlated in the past, due to their increasing bias with redshift - the bias reaches b = 100 at z = 2 for M > 5 x 10^13 h^-1 M_sun. The richness-dependent cluster correlation function, represented by the correlation scale versus cluster mean separation relation, R0-d, is generally consistent with observations. This relation can be approximated as R_0 = 1.7 d^0.6 h^-1 Mpc for d = 20 - 60 h^-1 Mpc. The R0-d relation exhibits surprisingly little evolution with redshift for z < 2; this can provide a new test of the current LCDM model when compared with future observations of high redshift clusters.

The ellipticities and alignments of clusters of galaxies, and their evolution with redshift, are examined in the context of a Lambda-dominated cold dark matter cosmology. We use a large-scale, high-resolution N-body simulation to model the matter distribution in a light cone containing ~10^6 clusters out to redshifts of z=3. Cluster ellipticities are determined as a function of mass, radius, and redshift, both in 3D and in projection. We find strong cluster ellipticities: the mean ellipticity increases with redshift from 0.3 at z=0 to 0.5 at z=3, for both 3D and 2D ellipticities; the evolution is well-fit by e=0.33+0.05z. The ellipticities increase with cluster mass and with cluster radius; the main cluster body is more elliptical than the cluster cores, but the increase of ellipticities with redshift is preserved. Using the fitted cluster ellipsoids, we determine the alignment of clusters as a function of their separation. We find strong alignment of clusters for separations <100 Mpc/h; the alignment increases with decreasing separation and with increasing redshift. The evolution of clusters from highly aligned and elongated systems at early times to lower alignment and elongation at present reflects the hierarchical and filamentary nature of structure formation. These measures of cluster ellipticity and alignment will provide a new test of the current cosmological model when compared with upcoming cluster surveys.

Giant arcs in strong lensing galaxy clusters can provide a purely geometric determination of cosmological parameters, such as the dark energy density and equation of state. We investigate sources of noise in cosmography with giant arcs, focusing in particular on errors induced by density fluctuations along the line-of-sight, and errors caused by modeling uncertainties. We estimate parameter errors in two independent ways, first by developing a Fisher matrix formalism for strong lensing parameters, and next by directly ray-tracing through N-body simulations using a multi-plane lensing code. We show that for reasonable power spectra, density fluctuations from large-scale structure produce > 100% errors in cosmological parameters derived from any single sightline, precluding the use of individual clusters or golden lenses to derive accurate cosmological constraints. Modeling uncertainties similarly can lead to large errors, and we show that the use of parametrized mass models in fitting strong lensing clusters can significantly bias the inferred cosmological parameters. We lastly speculate on means by which these errors may be corrected.

We explore several physical effects on the power spectrum of the Lyman-alpha forest transmitted flux. The effects we investigate here are usually not part of hydrodynamic simulations and so need to be estimated separately. The most important effect is that of high column density absorbers with damping wings, which add power on large scales. We compute their effect using the observational constraints on their abundance as a function of column density. Ignoring their effect leads to an underestimation of the slope of the linear theory power spectrum. The second effect we investigate is that of fluctuations in the ionizing radiation field. For this purpose we use a very large high resolution N-body simulation, which allows us to simulate both the fluctuations in the ionizing radiation and the small scale LyaF within the same simulation. We find an enhancement of power on large scales for quasars and a suppression for galaxies. The strength of the effect rapidly increases with increasing redshift, allowing it to be uniquely identified in cases where it is significant. We develop templates which can be used to search for this effect as a function of quasar lifetime, quasar luminosity function, and attenuation length. Finally, we explore the effects of galactic winds using hydrodynamic simulations. We find the wind effects on the LyaF power spectrum to be be degenerate with parameters related to the temperature of the gas that are already marginalized over in cosmological fits. While more work is needed to conclusively exclude all possible systematic errors, our results suggest that, in the context of data analysis procedures where parameters of the LyaF model are properly marginalized over, the flux power spectrum is a reliable tracer of cosmological information.

To date, in almost all strong gravitational lensing analyses for modeling giant arc systems and multiple quasar images, it has been assumed that all the deflecting matter is concentrated in one lens plane at a certain distance - the thin lens approximation. However, in a few observed cases, lenses at more than one redshift have been identified as contributing to the image splitting. Here we report on a quantitative investigation of the importance and frequency of significant multiple lensing agents. We use multi-lens plane simulations to evaluate how frequently two or more lens planes combined are essential for multiple imaging, as compared with the cases where a single lens plane alone provides enough focusing to be supercritical. We find that the fraction of cases for which more than one lens plane contributes significantly to a multi-image lensing situation is a strong function of source redshift. For sources at redshift unity, 95% of lenses involve only a single mass concentration, but for a more typical scenario with, e.g., a source at a redshift of z_s = 3.8, as many as 38% of the strongly lensed quasars/arcs occur because of a significant matter contribution from one or more ADDITIONAL lens planes. In the 30% to 40% of cases when additional planes make a significant contribution, the surface mass density of the primary lens will be overestimated by about 15% to 20%, if the additional contributions are not recognized.

We compute the properties of dark matter halos with mass $10^{6.5}-10^9\msun$ at redshift $z=6-11$ in the standard cold dark matter cosmological model, utilizing a very high resolution N-body simulation. We find that dark matter halos in these mass and redshift ranges are significantly biased over matter with a bias factor in the range 2-6. The dark matter halo mass function displays a slope of $2.05\pm 0.15$ at the small mass end. We do not find a universal dark matter density profile. Instead, we find a significant dependence of the central density profile of dark matter halos on halo mass and epoch with $\alpha_0=0.4-1.0$; the high-mass ($M\ge 10^8\msun$) low-redshift ($z\sim 6$) halos occupy the high end of the range and low-mass ($M\sim 10^{7}\msun$) high-redshift ($z\sim 11$) halos occupy the low end. Additionally, for fixed mass and epoch there is a significant dispersion in $\alpha_0$ due to the stochastic assembly of halos. Our results fit a relationship of the form $\alpha_0=0.75((1+z)/7.0)^{-1.25}(M/10^7\msun)^{0.11(1+z)/7.0}$ with a dispersion about this fit of $\pm 0.5$ and no systematic dependence of variance correlated with environment. The median spin parameter of dark matter halos is $0.03-0.04$ but with a large lognormal dispersion of $\sim 0.4$. Various quantities are tabulated or fitted with empirical formulae.

We use $N$-body simulations to measure mass functions in flat cosmological models with quintessence characterized by constant $w$ with $w=-1$, -2/3 and -1/2. The results are compared to the predictions of the formula proposed by Jenkins et al. at different redshifts, in terms of FOF masses as well as Abell masses appropriate for direct comparison to observations. The formula reproduces quite well the mass functions of simulated haloes in models with quintessence. We use the cluster mass function data at a number of redshifts from Carlberg et al. to constrain $\Omega_0$, $\sigma_8$ and $w$. The best fit is obtained in the limit $w \to 0$, but none of the values of $w$ in the considered range $-1 \le w < 0$ can actually be excluded. However, the adopted value of $w$ affects significantly the constraints in the $\Omega_0-\sigma_8$ plane. Taking into account the dependence on $w$ we find $\Omega_0=0.32 \pm 0.15$ and $\sigma_8=0.85_{-0.12}^{+0.38}$ (68% c.l.). Since less negative $w$ push the confidence regions toward higher $\Omega_0$ and lower $\sigma_8$ we conclude that relaxing the assumption of $w=-1$ typically made in such comparisons may resolve the discrepancy between recent cluster mass function results (yielding rather low $\Omega_0$ and high $\sigma_8$) and most other estimates. The fact that high $w$ values are preferred may however also point towards some unknown systematics in the data or the model with constant $w$ being inadequate.

We study the evolution of the cluster correlation function and its richness-dependence from z = 0 to z = 3 using large-scale cosmological simulations. A standard flat LCDM model with \Omega_m = 0.3 and, for comparison, a tilted \Omega_m = 1 model, TSCDM, are used. The evolutionary predictions are presented in a format suitable for direct comparisons with observations. We find that the cluster correlation strength increases with redshift: high redshift clusters are clustered more strongly (in comoving scale) than low redshift clusters of the same mass. The increased correlations with redshift, in spite of the decreasing mass correlation strength, is caused by the strong increase in cluster bias with redshift: clusters represent higher density peaks of the mass distribution as the redshift increases. The richness-dependent cluster correlation function, presented as the correlation-scale versus cluster mean separation relation, R_0 - d, is found to be, remarkably, independent of redshift to z <~ 2 for LCDM and z <~ 1 for TCDM (for a fixed correlation function slope and cluster mass within a fixed comoving radius). The non-evolving R_0 - d relation implies that both the comoving clustering scale and the cluster mean separation increase with redshift for the same mass clusters so that the R_0 - d relation remains essentially unchanged. The evolution of the R_0 - d relation from z ~ 0 to z ~ 3 provides an important new tool in cosmology; it can be used to break degeneracies that exist at z ~ 0 and provide precise determination of cosmological parameters.

The cluster correlation function and its richness dependence are determined from 1108 clusters of galaxies -- the largest sample of clusters studied so far -- found in 379 deg^2 of Sloan Digital Sky Survey early data. The results are compared with previous samples of optically and X-ray selected clusters. The richness-dependent correlation function increases monotonically from an average correlation scale of ~ 12 h^{-1} Mpc for poor clusters to ~ 25 h^{-1} Mpc for the richer, more massive clusters with a mean separation of ~ 90 h^{-1} Mpc. X-ray selected clusters suggest slightly stronger correlations than optically selected clusters (~ 2-\sigma). The results are compared with large-scale cosmological simulations. The observed richness-dependent cluster correlation function is well represented by the standard flat LCDM model (\Omega_m ~= 0.3, h ~= 0.7), and is inconsistent with the considerably weaker correlations predicted by \Omega_m = 1 models. An analytic relation for the correlation scale versus cluster mean separation, r_0 - d, that best describes the observations and the LCDM prediction is r_0 ~= 2.6 \sqrt{d} (for d ~= 20 - 90 h^{-1} Mpc). Data from the complete Sloan Digital Sky Survey, when available, will greatly enhance the accuracy of the results and allow a more precise determination of cosmological parameters.

The frequency of giant arcs - highly distorted and strongly gravitationally lensed background galaxies - is a powerful test for cosmological models. Previous comparisons of arc statistics for the currently favored concordance cosmological model (a flat lambda-dominated universe) with observations have shown an apparently large discrepancy. We present here new ray-shooting results, based on a high resolution (1024^3 particles in a 320 Mpc/h box) large-scale structure simulation normalized to the WMAP observations. We follow light rays through a pseudo-3D matter distribution approximated by up to 38 lens planes, and evaluate the occurrence of arcs for various source redshifts. We find that the frequency of strongly lensed background sources is a steep function of source redshift: the optical depth for giant arcs increases by a factor of five when the background sources are moved from redshift z_s = 1.0 to z_s = 1.5. This is a consequence of a moderate decrease of the critical surface mass density for lensing, combined with the very steep cluster mass function at the high mass end. Our results are consistent with those of Bartelmann et al. (1998) if we - as they did - restrict all sources to be exactly at z_s = 1. But if we allow for a more realistic distribution of source redshifts extending to or beyond z_s > 1.5, the apparent discrepancy vanishes: the frequency of arcs is increased by about a factor of ten as compared to previous estimates, and results in roughly one arc per 20 square degrees over the sky. This prediction for an LCDM model is then in good agreement with the observed frequency of arcs. Hence we consider the ``missing arc'' problem for a concordance LCDM cosmology to be solved.

We determine the linear amplitude of mass fluctuations in the universe, sigma_8, from the abundance of massive clusters at redshifts z=0.5 to 0.8. The evolution of massive clusters depends exponentially on the amplitude of mass fluctuations and thus provides a powerful measure of this important cosmological parameter. The relatively high abundance of massive clusters observed at z>0.5, and the relatively slow evolution of their abundance with time, suggest a high amplitude of mass fluctuations: sigma_8=0.9 +-10% for Omega_m=0.4, increasing slightly to sigma_8=0.95 for Omega_m=0.25 and sigma_8=1.0 for Omega_m=0.1 (flat CDM models). We use the cluster abundance observed at z=0.5 to 0.8 to derive a normalization relation from the high-redshift clusters, which is only weakly dependent on Omega_m: sigma_8*Omega_m^0.14 = 0.78 +-0.08. When combined with recent constraints from the present-day cluster mass function (sigma_8*Omega_m^0.6=0.33 +-0.03) we find sigma_8=0.98 +-0.1 and Omega_m=0.17 +-0.05. Low sigma_8 values (<0.7) are unlikely; they produce an order of magnitude fewer massive clusters than observed.

An improved implementation of an N-body code for simulating collisionless cosmological dynamics is presented. TPM (Tree-Particle-Mesh) combines the PM method on large scales with a tree code to handle particle-particle interactions at small separations. After the global PM forces are calculated, spatially distinct regions above a given density contrast are located; the tree code calculates the gravitational interactions inside these denser objects at higher spatial and temporal resolution. The new implementation includes individual particle time steps within trees, an improved treatment of tidal forces on trees, new criteria for higher force resolution and choice of time step, and parallel treatment of large trees. TPM is compared to P^3M and a tree code (GADGET) and is found to give equivalent results in significantly less time. The implementation is highly portable (requiring a Fortran compiler and MPI) and efficient on parallel machines. The source code can be found at http://astro.princeton.edu/~bode/TPM/

We present constraints on several cosmological parameters from a combined analysis of the most recent Cosmic Microwave Background anisotropy data and the Sloan Digital Sky Survey cluster mass function. We find that the combination of the two data sets breaks several degeneracies among the parameters and provides the following constraints: $\sigma_8=0.76\pm0.09$, $\Omega_m=0.26^{+0.06}_{-0.07}$, $h=0.66^{+0.05}_{-0.06}$, $n=0.96 \pm 0.05$, $\tau_c=0.07^{+0.07}_{-0.05}$.

The mass function of clusters of galaxies is determined from 400 deg^2 of early commissioning imaging data of the Sloan Digital Sky Survey; ~300 clusters in the redshift range z = 0.1 - 0.2 are used. Clusters are selected using two independent selection methods: a Matched Filter and a red-sequence color magnitude technique. The two methods yield consistent results. The cluster mass function is compared with large-scale cosmological simulations. We find a best-fit cluster normalization relation of sigma_8*omega_m^0.6 = 0.33 +- 0.03 (for 0.1 ~< omega_m ~< 0.4), or equivalently sigma_8 = (0.16/omega_m)^0.6. The amplitude of this relation is significantly lower than the previous canonical value, implying that either omega_m is lower than previously expected (omega_m = 0.16 if sigma_8 = 1) or sigma_8 is lower than expected (sigma_8 = 0.7 if omega_m = 0.3). The best-fit mass function parameters are omega_m = 0.19 (+0.08,-0.07) and sigma_8 = 0.9 (+0.3,-0.2). High values of omega_m (>= 0.4) and low sigma_8 (=< 0.6) are excluded at >~ 2 sigma.

High-resolution N-body simulations of four popular Cold Dark Matter cosmologies (LCDM, OCDM, QCDM, and tilted SCDM), each containing 10^5 clusters of galaxies in a cubic gigaparsec volume, are used to determine the evolution of the cluster mass function from z=3 to z=0. The large volume and high resolution of these simulations allow an accurate measure of the evolution of cosmologically important (but rare) massive clusters at high redshift. The simulated mass function is presented for cluster masses within several radii typically used observationally (R=0.5, 1.0, and 1.5 Mpc/h, both comoving and physical) in order to enable direct comparison with current and future observations. The simulated evolution is compared with current observations of massive clusters at redshifts 0.3<z<0.8. The Omega_m=1 tilted SCDM model, which exhibits very rapid evolution of the cluster abundance, produces too few clusters at z>0.3 and no massive clusters at z>0.5, in stark contradiction with observations. The Omega_m=0.3 models- LCDM, OCDM, and QCDM- all exhibit considerably weaker evolution and are consistent with current data. Among these low density models, OCDM evolves the least. These trends are enhanced at high redshift and can be used to discriminate between flat and open low density models. The simulated mass functions are compared with the Press-Schechter approximation. Standard Press-Schechter predicts too many low mass clusters at z=0, and too few clusters at higher redshift. We modify the approximation by a simple parameterization of the density contrast threshold for collapse, which has a redshift dependence. This modified Press-Schechter approximation provides a good fit to the simulated mass functions.

Discrepancies have emerged between the predictions of standard cold dark matter (CDM) theory and observations of clustering on sub-galactic scales. Warm dark matter (WDM) is a simple modification of CDM in which the dark matter particles have initial velocities due either to their having decoupled as thermal relics, or having been formed via non-equilibrium decay. We investigate the nonlinear gravitational clustering of WDM with a high resolution N-body code, and identify a number of distinctive observational signatures. Relative to CDM, halo concentrations and core densities are lowered, core radii are increased, and large halos emerge with far fewer low mass satellites. The number of small halos is suppressed, and those present are formed by `top down' fragmentation of caustics, as part of a `cosmic web' connecting massive halos. Few small halos form outside this web. If we identify small halos with dwarf galaxies, their number, spatial distribution, and formation epoch appear in better agreement with the observations for WDM than they are for CDM.

The Tree-Particle-Mesh (TPM) N-body algorithm couples the tree algorithm for directly computing forces on particles in an hierarchical grouping scheme with the extremely efficient mesh based PM structured approach. The combined TPM algorithm takes advantage of the fact that gravitational forces are linear functions of the density field. Thus one can use domain decomposition to break down the density field into many separate high density regions containing a significant fraction of the mass but residing in a very small fraction of the total volume. In each of these high density regions the gravitational potential is computed via the tree algorithm supplemented by tidal forces from the external density distribution. For the bulk of the volume, forces are computed via the PM algorithm; timesteps in this PM component are large compared to individually determined timesteps in the tree regions. Since each tree region can be treated independently, the algorithm lends itself to very efficient parallelization using message passing. We have tested the new TPM algorithm (a refinement of that originated by Xu 1995) by comparison with results from Ferrell & Bertschinger's P^3M code and find that, except in small clusters, the TPM results are at least as accurate as those obtained with the well-established P^3M algorithm, while taking significantly less computing time. Production runs of 10^9 particles indicate that the new code has great scientific potential when used with distributed computing resources.

We present simple analytic approximations for the linear and fully evolved nonlinear mass power spectrum for spatially flat cold dark matter (CDM) cosmological models with quintessence (Q). Quintessence is a time evolving, spatially inhomogeneous energy component with negative pressure and an equation of state w_Q < 0. It clusters gravitationally on large length scales but remains smooth like the cosmological constant on small length scales. We show that the clustering scale is determined by the Compton wavelength of the Q-field and derive a shape parameter, \Gamma_Q, to characterize the linear mass power spectrum. The growth of linear perturbations as functions of redshift, w_Q, and matter density \Omega_m is also quantified. Calibrating to N-body simulations, we construct a simple extension of the formula by Ma (1998) that closely approximates the nonlinear power spectrum for a range of plausible QCDM models.

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 have developed a code which links scalar-mode fluctuations in the early universe with those observable at the present time by integrating the coupled, linearized, Einstein, Boltzmann, and fluid equations in a perturbed flat Robertson-Walker spacetime. The results are useful both for calculations of the cosmic microwave background anisotropy and the linear power spectrum of matter fluctuations. This paper introduces the serial and parallel codes, presents timing results, and gives sample output in graphical and animated form. This is a latex version containing most of an HTML-format technical paper submitted for Supercomputing '95. The preferred method of viewing is to point your WWW client to http://arcturus.mit.edu/SC95/