Using Perturbative Least Action to Recover Cosmological Initial Conditions

We introduce a new method for generating initial conditions consistent with highly nonlinear observations of density and velocity fields. Using a variant of the Least Action method, called Perturbative Least Action (PLA), we show that it is possible to generate several different sets of initial conditions, each of which will satisfy a set of highly nonlinear observational constraints at the present day. We then discuss a code written to test and apply this method and present the results of several simulations.Comment: 24 pages, 6 postscript figures. Accepted for publication in Astrophysical Journa

Ionizing radiation fluctuations and large-scale structure in the Lyman-alpha forest

We investigate the large-scale inhomogeneities of the hydrogen ionizing radiation field in the Universe at redshift z=3. Using a raytracing algorithm, we simulate a model in which quasars are the dominant sources of radiation. We make use of large scale N-body simulations of a LambdaCDM universe, and include such effects as finite quasar lifetimes and output on the lightcone, which affects the shape of quasar light echoes. We create Lya forest spectra that would be generated in the presence of such a fluctuating radiation field, finding that the power spectrum of the Lya forest can be suppressed by as much as 15 % for modes with k=0.05-1 Mpc/h. This relatively small effect may have consequences for high precision measurements of the Lya power spectrum on larger scales than have yet been published. We also investigate another radiation field probe, the cross-correlation of quasar positions and the Lya forest. For both quasar lifetimes which we simulate (10^7 yr and 10^8 yr), we expect to see a strong decrease in the Lya absorption close to other quasars (the ``foreground'' proximity effect). We then use data from the Sloan Digital Sky Survey First Data Release to make an observational determination of this statistic. We find no sign of our predicted lack of absorption, but instead increased absorption close to quasars. If the bursts of radiation from quasars last on average < 10^6 yr, then we would not expect to be able to see the foreground effect. However, the strength of the absorption itself seems to be indicative of rare objects, and hence much longer total times of emission per quasar. Variability of quasars in bursts with timescales > 10^4yr and < 10^6 yr could reconcile these two facts.Comment: Submitted to ApJ, 21 pages, 17 postscript figures, emulateapj.st

Constraining quasar host halo masses with the strength of nearby Lyman-alpha forest absorption

Using cosmological hydrodynamic simulations we measure the mean transmitted flux in the Lyman alpha forest for quasar sightlines that pass near a foreground quasar. We find that the trend of absorption with pixel-quasar separation distance can be fitted using a simple power law form including the usual correlation function parameters r_{0} and \gamma so that ( = \sum exp(-tau_eff*(1+(r/r_{0})^(-\gamma)))). From the simulations we find the relation between r_{0} and quasar mass and formulate this as a way to estimate quasar host dark matter halo masses, quantifying uncertainties due to cosmological and IGM parameters, and redshift errors. With this method, we examine data for ~3000 quasars from the Sloan Digital Sky Survey (SDSS) Data Release 3, assuming that the effect of ionizing radiation from quasars (the so-called transverse proximity effect) is unimportant (no evidence for it is seen in the data.) We find that the best fit host halo mass for SDSS quasars with mean redshift z=3 and absolute G band magnitude -27.5 is log10(M/M_sun) = 12.48^{+0.53}_{-0.89}. We also use the Lyman-Break Galaxy (LBG) and Lyman alpha forest data of Adelberger et al in a similar fashion to constrain the halo mass of LBGs to be log10(M/M_sun) = 11.13^{+0.39}_{-0.55}, a factor of ~20 lower than the bright quasars. In addition, we study the redshift distortions of the Lyman alpha forest around quasars, using the simulations. We use the quadrupole to monopole ratio of the quasar-Lyman alpha forest correlation function as a measure of the squashing effect. We find that this does not have a measurable dependence on halo mass, but may be useful for constraining cosmic geometry.Comment: 10 pages, 11 figures, submitted to MNRA

Global simulations of aerosol processing in clouds

An explicit and detailed representation of in-droplet and in-crystal aerosol particles in stratiform clouds has been introduced in the global aerosol-climate model ECHAM5-HAM. The new scheme allows an evaluation of the cloud cycling of aerosols and an estimation of the relative contributions of nucleation and collision scavenging, as opposed to evaporation of hydrometeors in the global aerosol processing by clouds. On average an aerosol particle is cycled through stratiform clouds 0.5 times. The new scheme leads to important changes in the simulated fraction of aerosol scavenged in clouds, and consequently in the aerosol wet deposition. In general, less aerosol is scavenged into clouds with the new prognostic treatment than what is prescribed in standard ECHAM5-HAM. Aerosol concentrations, size distributions, scavenged fractions and cloud droplet concentrations are evaluated and compared to different observations. While the scavenged fraction and the aerosol number concentrations in the marine boundary layer are well represented in the new model, aerosol optical thickness, cloud droplet number concentrations in the marine boundary layer and the aerosol volume in the accumulation and coarse modes over the oceans are overestimated. Sensitivity studies suggest that a better representation of below-cloud scavenging, higher in-cloud collision coefficients, or a reduced water uptake by seasalt aerosols could reduce these biases

Recovering the Primordial Density Fluctuations: A comparison of methods

We present a comparative study of six different methods for reversing the gravitational evolution of a cosmological density field to recover the primordial fluctuations: linear theory, the Gaussianization mapping scheme, two different quasi-linear dynamical schemes based on the Zel'dovich approximation, a Hybrid dynamical-Gaussianization method and the Path Interchange Zel'dovich Approximation (PIZA). The final evolved density field from an N-body simulation constitutes our test case. We use a variety of statistical measures to compare the initial density field recovered from it to the true initial density field, using each of the six different schemes. These include point-by-point comparisons of the density fields in real space, the individual modes in Fourier space, as well as global statistical properties such as the genus, the PDF of the density, and the distribution of peak heights and their shapes. We find linear theory to be the most inaccurate of all the schemes. The Gaussianization scheme is the least accurate after linear theory. The two quasi-linear dynamical schemes are more accurate than Gaussianization, although they break down quite drastically when used outside their range of validity - the quasi-linear regime. The complementary beneficial aspects of the dynamical and the Gaussianization schemes are combined in the Hybrid method. We find this Hybrid scheme to be more accurate and robust than either Gaussianization or the dynamical method alone. The PIZA scheme performs substantially better than the others in all point-by-point comparisons. However, it produces an oversmoothed initial density field, with a smaller number of peaks than expected, but recovers the PDF of the initial density with impressive accuracy on scales as small as 3Mpc/h.Comment: 39 pages, including 13 Figures, submitted to Ap

Generation of Vorticity and Velocity Dispersion by Orbit Crossing

We study the generation of vorticity and velocity dispersion by orbit crossing using cosmological numerical simulations, and calculate the backreaction of these effects on the evolution of large-scale density and velocity divergence power spectra. We use Delaunay tessellations to define the velocity field, showing that the power spectra of velocity divergence and vorticity measured in this way are unbiased and have better noise properties than for standard interpolation methods that deal with mass weighted velocities. We show that high resolution simulations are required to recover the correct large-scale vorticity power spectrum, while poor resolution can spuriously amplify its amplitude by more than one order of magnitude. We measure the scalar and vector modes of the stress tensor induced by orbit crossing using an adaptive technique, showing that its vector modes lead, when input into the vorticity evolution equation, to the same vorticity power spectrum obtained from the Delaunay method. We incorporate orbit crossing corrections to the evolution of large scale density and velocity fields in perturbation theory by using the measured stress tensor modes. We find that at large scales (k~0.1 h/Mpc) vector modes have very little effect in the density power spectrum, while scalar modes (velocity dispersion) can induce percent level corrections at z=0, particularly in the velocity divergence power spectrum. In addition, we show that the velocity power spectrum is smaller than predicted by linear theory until well into the nonlinear regime, with little contribution from virial velocities.Comment: 27 pages, 14 figures. v2: reorganization of the material, new appendix. Accepted by PR

Evolution of the Cluster Mass and Correlation Functions in LCDM Cosmology

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.Comment: 20 pages, 9 figures, accepted for publication in Ap

Spatial Correlation Function of X-ray Selected AGN

We present a detailed description of the first direct measurement of the spatial correlation function of X-ray selected AGN. This result is based on an X-ray flux-limited sample of 219 AGN discovered in the contiguous 80.7 deg^2 region of the ROSAT North Ecliptic Pole (NEP) Survey. Clustering is detected at the 4 sigma level at comoving scales in the interval r = 5-60 h^-1 Mpc. Fitting the data with a power law of slope gamma=1.8, we find a correlation length of r_0 = 7.4 (+1.8, -1.9) h^-1 Mpc (Omega_M=0.3, Omega_Lambda=0.7). The median redshift of the AGN contributing to the signal is z_xi=0.22. This clustering amplitude implies that X-ray selected AGN are spatially distributed in a manner similar to that of optically selected AGN. Furthermore, the ROSAT NEP determination establishes the local behavior of AGN clustering, a regime which is poorly sampled in general. Combined with high-redshift measures from optical studies, the ROSAT NEP results argue that the AGN correlation strength essentially does not evolve with redshift, at least out to z~2.2. In the local Universe, X-ray selected AGN appear to be unbiased relative to galaxies and the inferred X-ray bias parameter is near unity, b_X~1. Hence X-ray selected AGN closely trace the underlying mass distribution. The ROSAT NEP AGN catalog, presented here, features complete optical identifications and spectroscopic redshifts. The median redshift, X-ray flux, and X-ray luminosity are z=0.41, f_X=1.1*10^-13 cgs, and L_X=9.2*10^43 h_70^-2 cgs (0.5-2.0 keV), respectively. Unobscured, type 1 AGN are the dominant constituents (90%) of this soft X-ray selected sample of AGN.Comment: 17 pages, 8 figures, accepted for publication in ApJ, a version with high-resolution figures is available at http://www.eso.org/~cmullis/papers/Mullis_et_al_2004b.ps.gz, a machine-readable version of the ROSAT NEP AGN catalog is available at http://www.eso.org/~cmullis/research/nep-catalog.htm

The richness dependence of galaxy cluster correlations: Results from a redshift survey of rich APM clusters

We analyse the spatial clustering properties of a new catalogue of very rich galaxy clusters selected from the APM Galaxy Survey. These clusters are of comparable richness and space density to Abell Richness Class $\geq 1$ clusters, but selected using an objective algorithm from a catalogue demonstrably free of artificial inhomogeneities. Evaluation of the two-point correlation function $\xi_{cc}(r)$ for the full sample and for richer subsamples reveals that the correlation amplitude is consistent with that measured for lower richness APM clusters and X-ray selected clusters. We apply a maxmimum likelihood estimator to find the best fitting slope and amplitude of a power law fit to $\xi_{cc}(r)$, and to estimate the correlation length $r_{0}$ (the value of $r$ at which $\xi_{cc}(r)$ is equal to unity). For clusters with a mean space density of 1.6\times 10^{-6}\hmpccc (equivalent to the space density of Abell Richness $\geq 2$ clusters), we find r_{0}=21.3^{+11.1}_{-9.3} \hmpc (95% confidence limits). This is consistent with the weak richness dependence of $\xi_{cc}(r)$ expected in Gaussian models of structure formation. In particular, the amplitude of $\xi_{cc}(r)$ at all richnesses matches that of $\xi_{cc}(r)$ for clusters selected in N-Body simulations of a low density Cold Dark Matter model.Comment: MNRAS submitted, 9 pages, LaTeX (mn), 7 figures. Also available at http://www-astronomy.mps.ohio-state.edu/~racc