Reconstruction of cosmological initial conditions from galaxy redshift catalogues

We present and test a new method for the reconstruction of cosmological initial conditions from a full-sky galaxy catalogue. This method, called ZTRACE, is based on a self-consistent solution of the growing mode of gravitational instabilities according to the Zel'dovich approximation and higher order in Lagrangian perturbation theory. Given the evolved redshift-space density field, smoothed on some scale, ZTRACE finds via an iterative procedure, an approximation to the initial density field for any given set of cosmological parameters; real-space densities and peculiar velocities are also reconstructed. The method is tested by applying it to N-body simulations of an Einstein-de Sitter and an open cold dark matter universe. It is shown that errors in the estimate of the density contrast dominate the noise of the reconstruction. As a consequence, the reconstruction of real space density and peculiar velocity fields using non-linear algorithms is little improved over those based on linear theory. The use of a mass-preserving adaptive smoothing, equivalent to a smoothing in Lagrangian space, allows an unbiased (although noisy) reconstruction of initial conditions, as long as the (linearly extrapolated) density contrast does not exceed unity. The probability distribution function of the initial conditions is recovered to high precision, even for Gaussian smoothing scales of ~ 5 Mpc/h, except for the tail at delta >~ 1. This result is insensitive to the assumptions of the background cosmology.Comment: 19 pages, MN style, 12 figures included, revised version. MNRAS, in pres

Anthropic versus cosmological solutions to the coincidence problem

In this paper we investigate possible solutions to the coincidence problem in flat phantom dark energy models with a constant dark energy equation of state and quintessence models with a linear scalar field potential. These models are representative of a broader class of cosmological scenarios in which the universe has a finite lifetime. We show that, in the absence of anthropic constraints, including a prior probability for the models inversely proportional to the total lifetime of the universe excludes models very close to the $\Lambda {\rm CDM}$ model. This relates a cosmological solution to the coincidence problem with a dynamical dark energy component having an equation of state parameter not too close to -1 at the present time. We further show, that anthropic constraints, if they are sufficiently stringent, may solve the coincidence problem without the need for dynamical dark energy.Comment: 7 pages, 7 figure

An Analytical Approach to Inhomogeneous Structure Formation

We develop an analytical formalism that is suitable for studying inhomogeneous structure formation, by studying the joint statistics of dark matter halos forming at two points. Extending the Bond et al. (1991) derivation of the mass function of virialized halos, based on excursion sets, we derive an approximate analytical expression for the ``bivariate'' mass function of halos forming at two redshifts and separated by a fixed comoving Lagrangian distance. Our approach also leads to a self-consistent expression for the nonlinear biasing and correlation function of halos, generalizing a number of previous results including those by Kaiser (1984) and Mo & White (1996). We compare our approximate solutions to exact numerical results within the excursion-set framework and find them to be consistent to within 2% over a wide range of parameters. Our formalism can be used to study various feedback effects during galaxy formation analytically, as well as to simply construct observable quantities dependent on the spatial distribution of objects. A code that implements our method is publicly available at http://www.arcetri.astro.it/~evan/GeminiComment: 41 Pages, 11 figures, published in ApJ, 571, 585. Reference added, Figure 2 axis relabele

Evolution of the Pairwise Peculiar Velocity Distribution Function in Lagrangian Perturbation Theory

The statistical distribution of the radial pairwise peculiar velocity of galaxies is known to have an exponential form as implied by observations and explicitly shown in N-body simulations. Here we calculate its statistical distribution function using the Zel'dovich approximation assuming that the primordial density fluctuations are Gaussian distributed. We show that the exponential distribution is realized as a transient phenomena on megaparsec scales in the standard cold-dark-matter model.Comment: 19 pages, 8 Postscript figures, AAS LaTe

Luminosity Density of Galaxies and Cosmic Star Formation Rate from Lambda-CDM Hydrodynamical Simulations

We compute the cosmic star formation rate (SFR) and the rest-frame comoving luminosity density in various pass-bands as a function of redshift using large-scale \Lambda-CDM hydrodynamical simulations with the aim of understanding their behavior as a function of redshift. To calculate the luminosity density of galaxies, we use an updated isochrone synthesis model which takes metallicity variations into account. The computed SFR and the UV-luminosity density have a steep rise from z=0 to 1, a moderate plateau between z=1 - 3, and a gradual decrease beyond z=3. The raw calculated results are significantly above the observed luminosity density, which can be explained either by dust extinction or the possibly inappropriate input parameters of the simulation. We model the dust extinction by introducing a parameter f; the fraction of the total stellar luminosity (not galaxy population) that is heavily obscured and thus only appears in the far-infrared to sub-millimeter wavelength range. When we correct our input parameters, and apply dust extinction with f=0.65, the resulting luminosity density fits various observations reasonably well, including the present stellar mass density, the local B-band galaxy luminosity density, and the FIR-to-submm extragalactic background. Our result is consistent with the picture that \sim 2/3 of the total stellar emission is heavily obscured by dust and observed only in the FIR. The rest of the emission is only moderately obscured which can be observed in the optical to near-IR wavelength range. We also argue that the steep falloff of the SFR from z=1 to 0 is partly due to the shock-heating of the universe at late times, which produces gas which is too hot to easily condense into star-forming regions.Comment: 25 pages, 6 figures. Accepted version in ApJ. Substantially revised from the previous version. More emphasis on the comparison with various observations and the hidden star formation by dust extinctio

Radiative Transfer Effects during Photoheating of the Intergalactic Medium

The thermal history of the intergalactic medium (IGM) after reionization is to a large extent determined by photoheating. Here we demonstrate that calculations of the photoheating rate which neglect radiative transfer effects substantially underestimate the energy input during and after reionization. The neglect of radiative transfer effects results in temperatures of the IGM which are too low by a factor of two after HeII reionization. We briefly discuss implications for the absorption properties of the IGM and the distribution of baryons in shallow potential wells.Comment: 4 pages, 2 figures, submitted to ApJ

The shapes, orientation, and alignment of Galactic dark matter subhalos

We present a study of the shapes, orientations, and alignments of Galactic dark matter subhalos in the ``Via Lactea'' simulation of a Milky Way-size LCDM host halo. Whereas isolated dark matter halos tend to be prolate, subhalos are predominantly triaxial. Overall subhalos are more spherical than the host halo, with minor to major and intermediate to major axis ratios of 0.68 and 0.83, respectively. Like isolated halos, subhalos tend to be less spherical in their central regions. The principal axis ratios are independent of subhalo mass, when the shapes are measured within a physical scale like r_Vmax, the radius of the peak of the circular velocity curve. Subhalos tend to be slightly more spherical closer to the host halo center. The spatial distribution of the subhalos traces the prolate shape of the host halo when they are selected by the largest V_max they ever had, i.e. before they experienced strong tidal mass loss. The subhalos' orientation is not random: the major axis tends to align with the direction towards the host halo center. This alignment disappears for halos beyond 3 r_200 and is more pronounced when the shapes are measured in the outer regions of the subhalos. The radial alignment is preserved during a subhalo's orbit and they become elongated during pericenter passage, indicating that the alignment is likely caused by the host halo's tidal forces. These tidal interactions with the host halo act to make subhalos rounder over time.Comment: 12 pages, 11 figures, submitted to ApJ, v2: corrected typo in abstract ("[...] subhalos tend be less spherical in their central regions."), added a few reference

The effect of Limber and flat-sky approximations on galaxy weak lensing

We review the effect of the commonly-used Limber and flat-sky approximations on the calculation of shear power spectra and correlation functions for galaxy weak lensing. These approximations are accurate at small scales, but it has been claimed recently that their impact on low multipoles could lead to an increase in the amplitude of the mass fluctuations inferred from surveys such as CFHTLenS, reducing the tension between galaxy weak lensing and the amplitude determined by Planck from observations of the cosmic microwave background. Here, we explore the impact of these approximations on cosmological parameters derived from weak lensing surveys, using the CFHTLenS data as a test case. We conclude that the use of small-angle approximations for cosmological parameter estimation is negligible for current data, and does not contribute to the tension between current weak lensing surveys and Planck

Formation of early-type galaxies from cosmological initial conditions

We describe high resolution Smoothed Particle Hydrodynamics (SPH) simulations of three approximately $M_*$ field galaxies starting from \LCDM initial conditions. The simulations are made intentionally simple, and include photoionization, cooling of the intergalactic medium, and star formation but not feedback from AGN or supernovae. All of the galaxies undergo an initial burst of star formation at $z \approx 5$, accompanied by the formation of a bubble of heated gas. Two out of three galaxies show early-type properties at present whereas only one of them experienced a major merger. Heating from shocks and -PdV work dominates over cooling so that for most of the gas the temperature is an increasing function of time. By $z \approx 1$ a significant fraction of the final stellar mass is in place and the spectral energy distribution resembles those of observed massive red galaxies. The galaxies have grown from $z=1 \to 0$ on average by 25% in mass and in size by gas poor (dry) stellar mergers. By the present day, the simulated galaxies are old ($\approx 10 {\rm Gyrs}$), kinematically hot stellar systems surrounded by hot gaseous haloes. Stars dominate the mass of the galaxies up to $\approx 4$ effective radii ($\approx 10$ kpc). Kinematic and most photometric properties are in good agreement with those of observed elliptical galaxies. The galaxy with a major merger develops a counter-rotating core. Our simulations show that realistic intermediate mass giant elliptical galaxies with plausible formation histories can be formed from \LCDM initial conditions even without requiring recent major mergers or feedback from supernovae or AGN.Comment: accepted for publication in Ap

Fundamental Discreteness Limitations of Cosmological N-Body Clustering Simulations

We explore some of the effects that discreteness and two-body scattering may have on N-body simulations with ``realistic'' cosmological initial conditions. We use an identical subset of particles from the initial conditions for a $128^3$ Particle-Mesh (PM) calculation as the initial conditions for a variety P$^3$M and Tree code runs. We investigate the effect of mass resolution (the mean interparticle separation) since most ``high resolution'' codes only have high resolution in gravitational force. The phase-insensitive two--point statistics, such as the power spectrum (autocorrelation) are somewhat affected by these variations, but phase-sensitive statistics show greater differences. Results converge at the mean interparticle separation scale of the lowest mass-resolution code. As more particles are added, but the force resolution is held constant, the P$^3$M and the Tree runs agree more and more strongly with each other and with the PM run which had the same initial conditions. This shows high particle density is necessary for correct time evolution, since many different results cannot all be correct. However, they do not so converge to a PM run which continued the fluctuations to small scales. Our results show that ignoring them is a major source of error on comoving scales of the missing wavelengths. This can be resolved by putting in a high particle density. Since the codes never agree well on scales below the mean comoving interparticle separation, we find little justification for quantitative predictions on this scale. Some measures vary by 50%, but others can be off by a factor of three or more. Our results suggest possible problems with the density of galaxy halos, formation of early generation objects such as QSO absorber clouds, etc.Comment: Revised version to be published in Astrophysical Journal. One figure changed; expanded discussion, more information on code parameters. Latex, 44 pages, including 19 figures. Higher resolution versions of Figures 10-15 available at: ftp://kusmos.phsx.ukans.edu/preprints/nbod