-
The Coyote Universe I: Precision Determination of the Nonlinear Matter Power Spectrum
Authors:
Katrin Heitmann,
Martin White,
Christian Wagner,
Salman Habib,
David Higdon
Abstract:
Near-future cosmological observations targeted at investigations of dark energy pose stringent requirements on the accuracy of theoretical predictions for the clustering of matter. Currently, N-body simulations comprise the only viable approach to this problem. In this paper we demonstrate that N-body simulations can indeed be sufficiently controlled to fulfill these requirements for the needs of…
▽ More
Near-future cosmological observations targeted at investigations of dark energy pose stringent requirements on the accuracy of theoretical predictions for the clustering of matter. Currently, N-body simulations comprise the only viable approach to this problem. In this paper we demonstrate that N-body simulations can indeed be sufficiently controlled to fulfill these requirements for the needs of ongoing and near-future weak lensing surveys. By performing a large suite of cosmological simulation comparison and convergence tests we show that results for the nonlinear matter power spectrum can be obtained at 1% accuracy out to k~1 h/Mpc. The key components of these high accuracy simulations are: precise initial conditions, very large simulation volumes, sufficient mass resolution, and accurate time stepping. This paper is the first in a series of three, with the final aim to provide a high-accuracy prediction scheme for the nonlinear matter power spectrum.
△ Less
Submitted 24 August, 2010; v1 submitted 5 December, 2008;
originally announced December 2008.
-
Simulations and cosmological inference: A statistical model for power spectra means and covariances
Authors:
Michael D. Schneider,
Lloyd Knox,
Salman Habib,
Katrin Heitmann,
David Higdon,
Charles Nakhleh
Abstract:
We describe an approximate statistical model for the sample variance distribution of the non-linear matter power spectrum that can be calibrated from limited numbers of simulations. Our model retains the common assumption of a multivariate Normal distribution for the power spectrum band powers, but takes full account of the (parameter dependent) power spectrum covariance. The model is calibrated…
▽ More
We describe an approximate statistical model for the sample variance distribution of the non-linear matter power spectrum that can be calibrated from limited numbers of simulations. Our model retains the common assumption of a multivariate Normal distribution for the power spectrum band powers, but takes full account of the (parameter dependent) power spectrum covariance. The model is calibrated using an extension of the framework in Habib et al. (2007) to train Gaussian processes for the power spectrum mean and covariance given a set of simulation runs over a hypercube in parameter space. We demonstrate the performance of this machinery by estimating the parameters of a power-law model for the power spectrum. Within this framework, our calibrated sample variance distribution is robust to errors in the estimated covariance and shows rapid convergence of the posterior parameter constraints with the number of training simulations.
△ Less
Submitted 21 September, 2008; v1 submitted 9 June, 2008;
originally announced June 2008.
-
X-ray Observations of Optically Selected Giant Elliptical-Dominated Galaxy Groups
Authors:
Alexey Voevodkin,
Christopher J. Miller,
Konstantin Borozdin,
Katrin Heitmann,
Salman Habib,
Paul Ricker,
Robert C. Nichol
Abstract:
We present a combined optical and X-ray analysis of three optically selected X-ray bright groups with giant elliptical galaxies in the center. These massive ellipticals were targeted for XMM-Newton X-ray observations based on their large velocity dispersions and their proximity to a nearby ROSAT X-ray source. Additionally, these targets are significantly brighter in the optical than their neares…
▽ More
We present a combined optical and X-ray analysis of three optically selected X-ray bright groups with giant elliptical galaxies in the center. These massive ellipticals were targeted for XMM-Newton X-ray observations based on their large velocity dispersions and their proximity to a nearby ROSAT X-ray source. Additionally, these targets are significantly brighter in the optical than their nearest neighbors. We show that one of these systems meets the standard criteria for a fossil group. While the other two systems have a prominent magnitude gap in the E/S0 ridgeline, they do not appear to have reached the fossil-like final stage of group evolution.
△ Less
Submitted 15 May, 2008; v1 submitted 26 March, 2008;
originally announced March 2008.
-
The Structure of Halos: Implications for Group and Cluster Cosmology
Authors:
Zarija Lukić,
Darren Reed,
Salman Habib,
Katrin Heitmann
Abstract:
The dark matter halo mass function is a key repository of cosmological information over a wide range of mass scales, from individual galaxies to galaxy clusters. N-body simulations have established that the friends-of-friends (FOF) mass function has a universal form to a surprising level of accuracy (< 10%). The high-mass tail of the mass function is exponentially sensitive to the amplitude of t…
▽ More
The dark matter halo mass function is a key repository of cosmological information over a wide range of mass scales, from individual galaxies to galaxy clusters. N-body simulations have established that the friends-of-friends (FOF) mass function has a universal form to a surprising level of accuracy (< 10%). The high-mass tail of the mass function is exponentially sensitive to the amplitude of the initial density perturbations, the mean matter density parameter, Omega_m, and to the dark energy controlled late-time evolution of the density field. Observed group and cluster masses, however, are usually stated in terms of a spherical overdensity (SO) mass which does not map simply to the FOF mass. Additionally, the widely used halo models of structure formation -- and halo occupancy distribution descriptions of galaxies within halos -- are often constructed exploiting the universal form of the FOF mass function. This again raises the question of whether FOF halos can be simply related to the notion of a spherical overdensity mass. By employing results from Monte Carlo realizations of ideal Navarro-Frenk-White (NFW) halos and N-body simulations, we study the relationship between the two definitions of halo mass. We find that the vast majority of halos (80-85%) in the mass-range 10^{12.5}-10^{15.5} M_sun/h indeed allow for an accurate mapping between the two definitions (~ 5%), but only if the halo concentrations are known. Nonisolated halos fall into two broad classes: those with complex substructure that are poor fits to NFW profiles and those ``bridged'' by the (isodensity-based)FOF algorithm. A closer investigation of the bridged halos reveals that the fraction of these halos and their satellite mass distribution is cosmology dependent. (abridged)
△ Less
Submitted 25 March, 2008;
originally announced March 2008.
-
Multiple Uncertainties in Time-Variant Cosmological Particle Data
Authors:
Steve Haroz,
Kwan-Liu Ma,
Katrin Heitmann
Abstract:
Though the mediums for visualization are limited, the potential dimensions of a dataset are not. In many areas of scientific study, understanding the correlations between those dimensions and their uncertainties is pivotal to mining useful information from a dataset. Obtaining this insight can necessitate visualizing the many relationships among temporal, spatial, and other dimensionalities of d…
▽ More
Though the mediums for visualization are limited, the potential dimensions of a dataset are not. In many areas of scientific study, understanding the correlations between those dimensions and their uncertainties is pivotal to mining useful information from a dataset. Obtaining this insight can necessitate visualizing the many relationships among temporal, spatial, and other dimensionalities of data and its uncertainties. We utilize multiple views for interactive dataset exploration and selection of important features, and we apply those techniques to the unique challenges of cosmological particle datasets. We show how interactivity and incorporation of multiple visualization techniques help overcome the problem of limited visualization dimensions and allow many types of uncertainty to be seen in correlation with other variables.
△ Less
Submitted 25 February, 2009; v1 submitted 15 January, 2008;
originally announced January 2008.
-
The Cosmic Code Comparison Project
Authors:
Katrin Heitmann,
Zarija Lukic,
Patricia Fasel,
Salman Habib,
Michael S. Warren,
Martin White,
James Ahrens,
Lee Ankeny,
Ryan Armstrong,
Brian O'Shea,
Paul M. Ricker,
Volker Springel,
Joachim Stadel,
Hy Trac
Abstract:
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…
▽ More
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.
△ Less
Submitted 8 June, 2007;
originally announced June 2007.
-
The Halo Mass Function: High-Redshift Evolution and Universality
Authors:
Zarija Lukic,
Katrin Heitmann,
Salman Habib,
Sergei Bashinsky,
Paul M. Ricker
Abstract:
We study the formation of dark matter halos in the concordance LCDM model over a wide range of redshifts, from z=20 to the present. Our primary focus is the halo mass function, a key probe of cosmology. By performing a large suite of nested-box N-body simulations with careful convergence and error controls (60 simulations with box sizes from 4 to 256 Mpc/h, we determine the mass function and its…
▽ More
We study the formation of dark matter halos in the concordance LCDM model over a wide range of redshifts, from z=20 to the present. Our primary focus is the halo mass function, a key probe of cosmology. By performing a large suite of nested-box N-body simulations with careful convergence and error controls (60 simulations with box sizes from 4 to 256 Mpc/h, we determine the mass function and its evolution with excellent statistical and systematic errors, reaching a few percent over most of the considered redshift and mass range. Across the studied redshifts, the halo mass is probed over 6 orders of magnitude (10^7 - 10^13.5 M_sun/h). Historically, there has been considerable variation in the high redshift mass function as obtained by different groups. We have made a concerted effort to identify and correct possible systematic errors in computing the mass function at high redshift and to explain the discrepancies between some of the previous results. We discuss convergence criteria for the required force resolution, simulation box size, halo mass range, initial and final redshift, and time stepping. Because of conservative cuts on the mass range probed by individual boxes, our results are relatively insensitive to simulation volume, the remaining sensitivity being consistent with extended Press-Schechter theory. Previously obtained mass function fits near z=0, when scaled by linear theory, are in good agreement with our results at all redshifts, although a mild redshift dependence consistent with that found by Reed and collaborators exists at low redshifts.
△ Less
Submitted 13 January, 2008; v1 submitted 14 February, 2007;
originally announced February 2007.
-
Cosmic Calibration: Constraints from the Matter Power Spectrum and the Cosmic Microwave Background
Authors:
Salman Habib,
Katrin Heitmann,
David Higdon,
Charles Nakhleh,
Brian Williams
Abstract:
Several cosmological measurements have attained significant levels of maturity and accuracy over the last decade. Continuing this trend, future observations promise measurements of the statistics of the cosmic mass distribution at an accuracy level of one percent out to spatial scales with k~10 h/Mpc and even smaller, entering highly nonlinear regimes of gravitational instability. In order to in…
▽ More
Several cosmological measurements have attained significant levels of maturity and accuracy over the last decade. Continuing this trend, future observations promise measurements of the statistics of the cosmic mass distribution at an accuracy level of one percent out to spatial scales with k~10 h/Mpc and even smaller, entering highly nonlinear regimes of gravitational instability. In order to interpret these observations and extract useful cosmological information from them, such as the equation of state of dark energy, very costly high precision, multi-physics simulations must be performed. We have recently implemented a new statistical framework with the aim of obtaining accurate parameter constraints from combining observations with a limited number of simulations. The key idea is the replacement of the full simulator by a fast emulator with controlled error bounds. In this paper, we provide a detailed description of the methodology and extend the framework to include joint analysis of cosmic microwave background and large scale structure measurements. Our framework is especially well-suited for upcoming large scale structure probes of dark energy such as baryon acoustic oscillations and, especially, weak lensing, where percent level accuracy on nonlinear scales is needed.
△ Less
Submitted 13 February, 2007;
originally announced February 2007.
-
Virial Scaling of Massive Dark Matter Halos: Why Clusters Prefer a High Normalization Cosmology
Authors:
A. E. Evrard,
J. Bialek,
M. Busha,
M. White,
S. Habib,
K. Heitmann,
M. Warren,
E. Rasia,
G. Tormen,
L. Moscardini,
C. Power,
A. R. Jenkins,
L. Gao,
C. S. Frenk,
V. Springel,
S. D. M. White,
J. Diemand
Abstract:
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…
▽ More
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/Ω_m, a value consistent with recent Sunyaev-Zel'dovich observations.
△ Less
Submitted 29 August, 2007; v1 submitted 8 February, 2007;
originally announced February 2007.
-
Cosmic Calibration
Authors:
Katrin Heitmann,
David Higdon,
Charles Nakhleh,
Salman Habib
Abstract:
The complexity and accuracy of current and future precision cosmology observational campaigns has made it essential to develop an efficient technique for directly combining simulation and observational datasets to determine cosmological and model parameters; a procedure we term calibration. Once a satisfactory calibration of the underlying cosmological model is achieved, independent predictions…
▽ More
The complexity and accuracy of current and future precision cosmology observational campaigns has made it essential to develop an efficient technique for directly combining simulation and observational datasets to determine cosmological and model parameters; a procedure we term calibration. Once a satisfactory calibration of the underlying cosmological model is achieved, independent predictions for new observations become possible. For this procedure to be effective, robust characterization of the uncertainty in the calibration process is highly desirable. In this Letter, we describe a statistical methodology which can achieve both of these goals. An application example based around dark matter structure formation simulations and a synthetic mass power spectrum dataset is used to demonstrate the approach.
△ Less
Submitted 7 June, 2006;
originally announced June 2006.
-
Nonequilibrium and Nonlinear Dynamics in Geomaterials I : The Low Strain Regime
Authors:
Donatella Pasqualini,
Katrin Heitmann,
James A. TenCate,
Salman Habib,
David Higdon,
Paul A. Johnson
Abstract:
Members of a wide class of geomaterials are known to display complex and fascinating nonlinear and nonequilibrium dynamical behaviors over a wide range of bulk strains, down to surprisingly low values, e.g., 10^{-7}. In this paper we investigate two sandstones, Berea and Fontainebleau, and characterize their behavior under the influence of very small external forces via carefully controlled reso…
▽ More
Members of a wide class of geomaterials are known to display complex and fascinating nonlinear and nonequilibrium dynamical behaviors over a wide range of bulk strains, down to surprisingly low values, e.g., 10^{-7}. In this paper we investigate two sandstones, Berea and Fontainebleau, and characterize their behavior under the influence of very small external forces via carefully controlled resonant bar experiments. By reducing environmental effects due to temperature and humidity variations, we are able to systematically and reproducibly study dynamical behavior at strains as low as 10^{-9}. Our study establishes the existence of two strain thresholds, the first, epsilon_L, below which the material is essentially linear, and the second, epsilon_M, below which the material is nonlinear but where quasiequilibrium thermodynamics still applies as evidenced by the success of Landau theory and a simple macroscopic description based on the Duffing oscillator. At strains above epsilon_M the behavior becomes truly nonequilibrium -- as demonstrated by the existence of material conditioning -- and Landau theory no longer applies. The main focus of this paper is the study of the region below the second threshold, but we also comment on how our work clarifies and resolves previous experimental conflicts, as well as suggest new directions of research.
△ Less
Submitted 6 January, 2006;
originally announced January 2006.
-
Capturing Halos at High Redshifts
Authors:
Katrin Heitmann,
Zarija Lukic,
Salman Habib,
Paul M. Ricker
Abstract:
We study the evolution of the mass function of dark matter halos in the concordance LCDM model at high redshift. We employ overlapping (multiple-realization) numerical simulations to cover a wide range of halo masses, 10^7-10^{15}h^{-1}M_\odot, with redshift coverage beginning at z=20. The Press-Schechter mass function is significantly discrepant from the simulation results at high redshifts. Of…
▽ More
We study the evolution of the mass function of dark matter halos in the concordance LCDM model at high redshift. We employ overlapping (multiple-realization) numerical simulations to cover a wide range of halo masses, 10^7-10^{15}h^{-1}M_\odot, with redshift coverage beginning at z=20. The Press-Schechter mass function is significantly discrepant from the simulation results at high redshifts. Of the more recently proposed mass functions, our results are in best agreement with Warren et al. (2005). The statistics of the simulations -- along with good control over systematics -- allow for fits accurate to the level of 20% at all redshifts. We provide a concise discussion of various issues in defining and computing the halo mass function, and how these are addressed in our simulations.
△ Less
Submitted 11 January, 2006;
originally announced January 2006.
-
Shapes and Sizes of Voids in the LCDM Universe: Excursion Set Approach
Authors:
Sergei Shandarin,
Hume A. Feldman,
Katrin Heitmann,
Salman Habib
Abstract:
We study the global distribution and morphology of dark matter voids in a LCDM universe using density fields generated by N-body simulations. Voids are defined as isolated regions of the low-density excursion set specified via density thresholds, the density thresholds being quantified by the corresponding filling factors, i.e., the fraction of the total volume in the excursion set. Our work enc…
▽ More
We study the global distribution and morphology of dark matter voids in a LCDM universe using density fields generated by N-body simulations. Voids are defined as isolated regions of the low-density excursion set specified via density thresholds, the density thresholds being quantified by the corresponding filling factors, i.e., the fraction of the total volume in the excursion set. Our work encompasses a systematic investigation of the void volume function, the volume fraction in voids, and the fitting of voids to corresponding ellipsoids and spheres. We emphasize the relevance of the percolation threshold to the void volume statistics of the density field both in the high redshift, Gaussian random field regime, as well as in the present epoch. By using measures such as the Inverse Porosity, we characterize the quality of ellipsoidal fits to voids, finding that such fits are a poor representation of the larger voids that dominate the volume of the void excursion set.
△ Less
Submitted 28 September, 2005;
originally announced September 2005.
-
Joint Efficient Dark-energy Investigation (JEDI): a Candidate Implementation of the NASA-DOE Joint Dark Energy Mission (JDEM)
Authors:
Arlin Crotts,
Peter Garnavich,
William Priedhorsky,
Salman Habib,
Katrin Heitmann,
Yun Wang,
Eddie Baron,
David Branch,
Harvey Moseley,
Alexander Kutyrev,
Chris Blake,
Edward Cheng,
Ian Dell'Antonio,
John Mackenty,
Gordon Squires,
Max Tegmark,
Craig Wheeler,
Ned Wright
Abstract:
JEDI will probe dark energy in three independent ways by measuring the expansion history of the universe: (1) using type Ia supernovae as cosmological standard candles over a range of distances, (2) using baryon oscillations as a cosmological standard ruler over a range of cosmic epochs, (3) mapping the weak gravitational lensing distortion by foreground galaxies of the images of background gala…
▽ More
JEDI will probe dark energy in three independent ways by measuring the expansion history of the universe: (1) using type Ia supernovae as cosmological standard candles over a range of distances, (2) using baryon oscillations as a cosmological standard ruler over a range of cosmic epochs, (3) mapping the weak gravitational lensing distortion by foreground galaxies of the images of background galaxies at different distances. JEDI will unravel the nature of dark energy with accuracy and precision.
JEDI is a 2m-class space telescope with the unique ability of simultaneous wide-field imaging (0.8-4.2 micron in five bands) and multi-slit spectroscopy (0.8-3.2 micron) with a field of view of 1 square degree. What makes JEDI efficient is its ability to simultaneously obtain high signal-to-noise ratio, moderate resolution slit spectra for all supernovae and ~ 5000 galaxies in its wide field of view, and to combine imaging and spectroscopy so that the appropriate balance of time is devoted to each.
JEDI will measure the cosmic expansion history H(z) as a free function to < 2% accuracy in redshift bins of 0.2-0.3. Assuming a flat universe and σ(Ω_m)=0.01 (0.03), JEDI could measure a departure from a vanilla ΛCDM model (w_0=-1, w'=0) with σ(w_0)=0.013 (0.031) and σ(w')=0.046 (0.063). JEDI will obtain the well-sampled lightcurves in Z, J, H, K, L bands and spectra of ~ 14,000 type Ia supernovae with redshifts ranging from 0 to 1.7; the redshifts of ~ 10-100 million galaxies to H ~ 23 and z ~ 4 over 1000-10,000 square degrees; and measurements of the shapes of galaxies over 1000-10,000 square degrees in Z,J,H,K,L for > 10^9 galaxies to H \~ 25.
△ Less
Submitted 14 July, 2005; v1 submitted 1 July, 2005;
originally announced July 2005.
-
Tightening Constraints from the Lyman Alpha Forest with the Flux Probability Distribution Function
Authors:
Adam Lidz,
Katrin Heitmann,
Lam Hui,
Salman Habib,
Michael Rauch,
Wallace L. W. Sargent
Abstract:
The analysis of the Lyman-alpha forest of absorption lines in quasar spectra has emerged as a potentially powerful technique to constrain the linear matter power spectrum. In most previous work, the amplitude of the ionizing background was fixed by calibrating simulations to match the observed mean transmitted flux in the Lyman-alpha forest. This procedure is undesirable in principle as it requi…
▽ More
The analysis of the Lyman-alpha forest of absorption lines in quasar spectra has emerged as a potentially powerful technique to constrain the linear matter power spectrum. In most previous work, the amplitude of the ionizing background was fixed by calibrating simulations to match the observed mean transmitted flux in the Lyman-alpha forest. This procedure is undesirable in principle as it requires the estimation of the unabsorbed quasar continuum level, a difficult undertaking subject to various sources of systematic error and bias. We suggest an alternative approach based on measuring the one-point probability distribution function (pdf) of the fluctuations in the flux about the mean, relative to the mean, i.e. the pdf of delta_f = (f - <f>)/<f>. This statistic, while sensitive to the amplitude of the ionizing background, has the virtue that its measurement does not require an estimate of the unabsorbed continuum level. We present a measurement of the pdf of delta_f from seven Keck HIRES spectra, spanning a redshift range of z = 2.2 - 4.4. To illustrate that our method is useful, we compare our measurements of the pdf of delta_f, and measurements of the flux power spectrum from Croft et al. (2002) at z = 2.72, with cosmological simulations. From this comparison, we obtain constraints on the mean transmission in the Lyman-alpha forest, the slope of the temperature-density relation, as well as the amplitude and slope of the mass power spectrum. Our methodology will be useful for obtaining more precise constraints with larger data samples from the Sloan Digital Sky Survey (SDSS).
△ Less
Submitted 6 May, 2005;
originally announced May 2005.
-
Inflationary Perturbations and Precision Cosmology
Authors:
Salman Habib,
Andreas Heinen,
Katrin Heitmann,
Gerard Jungman
Abstract:
Inflationary cosmology provides a natural mechanism for the generation of primordial perturbations which seed the formation of observed cosmic structure and lead to specific signals of anisotropy in the cosmic microwave background radiation. In order to test the broad inflationary paradigm as well as particular models against precision observations, it is crucial to be able to make accurate pred…
▽ More
Inflationary cosmology provides a natural mechanism for the generation of primordial perturbations which seed the formation of observed cosmic structure and lead to specific signals of anisotropy in the cosmic microwave background radiation. In order to test the broad inflationary paradigm as well as particular models against precision observations, it is crucial to be able to make accurate predictions for the power spectrum of both scalar and tensor fluctuations. We present detailed calculations of these quantities utilizing direct numerical approaches as well as error-controlled uniform approximations, comparing with the (uncontrolled) traditional slow-roll approach. A simple extension of the leading-order uniform approximation yields results for the power spectra amplitudes, the spectral indices, and the running of spectral indices, with accuracy of the order of 0.1% - approximately the same level at which the transfer functions are known. Several representative examples are used to demonstrate these results.
△ Less
Submitted 9 January, 2005;
originally announced January 2005.
-
Robustness of Cosmological Simulations I: Large Scale Structure
Authors:
Katrin Heitmann,
Paul M. Ricker,
Michael S. Warren,
Salman Habib
Abstract:
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 observa…
▽ More
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.
△ Less
Submitted 10 June, 2005; v1 submitted 30 November, 2004;
originally announced November 2004.
-
The Nonlinear Cosmological Matter Power Spectrum with Massive Neutrinos I: The Halo Model
Authors:
Kevork Abazajian,
Eric R. Switzer,
Scott Dodelson,
Katrin Heitmann,
Salman Habib
Abstract:
Measurements of the linear power spectrum of galaxies have placed tight constraints on neutrino masses. We extend the framework of the halo model of cosmological nonlinear matter clustering to include the effect of massive neutrino infall into cold dark matter (CDM) halos. The magnitude of the effect of neutrino clustering for three degenerate mass neutrinos with m_nu=0.9 eV is of order ~1%, wit…
▽ More
Measurements of the linear power spectrum of galaxies have placed tight constraints on neutrino masses. We extend the framework of the halo model of cosmological nonlinear matter clustering to include the effect of massive neutrino infall into cold dark matter (CDM) halos. The magnitude of the effect of neutrino clustering for three degenerate mass neutrinos with m_nu=0.9 eV is of order ~1%, within the potential sensitivity of upcoming weak lensing surveys. In order to use these measurements to further constrain--or eventually detect--neutrino masses, accurate theoretical predictions of the nonlinear power spectrum in the presence of massive neutrinos will be needed, likely only possible through high-resolution multiple particle (neutrino, CDM and baryon) simulations.
△ Less
Submitted 8 February, 2005; v1 submitted 19 November, 2004;
originally announced November 2004.
-
Inverse-Scattering Theory and the Density Perturbations from Inflation
Authors:
Salman Habib,
Katrin Heitmann,
Gerard Jungman
Abstract:
Inflationary cosmology provides a successful paradigm for solving several problems, notably the generation of density perturbations which seed the formation of observed cosmic structure. We show how to use inverse scattering theory as the basis for the inflationary reconstruction program, the goal of which is to gain information about the physics which drives inflation. Inverse scattering theory…
▽ More
Inflationary cosmology provides a successful paradigm for solving several problems, notably the generation of density perturbations which seed the formation of observed cosmic structure. We show how to use inverse scattering theory as the basis for the inflationary reconstruction program, the goal of which is to gain information about the physics which drives inflation. Inverse scattering theory provides an effective and well-motivated procedure, having a sound mathematical basis and being of sufficient generality that it can be considered the foundation for a non-parametric reconstruction program. We show how simple properties of the power spectrum translate directly into statements about the evolution of the background geometry during inflation.
△ Less
Submitted 24 September, 2004;
originally announced September 2004.
-
Characterizing Inflationary Perturbations: The Uniform Approximation
Authors:
Salman Habib,
Andreas Heinen,
Katrin Heitmann,
Gerard Jungman,
Carmen Molina-Paris
Abstract:
The spectrum of primordial fluctuations from inflation can be obtained using a mathematically controlled, and systematically extendable, uniform approximation. Closed-form expressions for power spectra and spectral indices may be found without making explicit slow-roll assumptions. Here we provide details of our previous calculations, extend the results beyond leading order in the approximation,…
▽ More
The spectrum of primordial fluctuations from inflation can be obtained using a mathematically controlled, and systematically extendable, uniform approximation. Closed-form expressions for power spectra and spectral indices may be found without making explicit slow-roll assumptions. Here we provide details of our previous calculations, extend the results beyond leading order in the approximation, and derive general error bounds for power spectra and spectral indices. Already at next-to-leading order, the errors in calculating the power spectrum are less than a per cent. This meets the accuracy requirement for interpreting next-generation CMB observations.
△ Less
Submitted 4 June, 2004;
originally announced June 2004.
-
Nonlinear and Nonequilibrium Dynamics in Geomaterials
Authors:
James A. TenCate,
Donatella Pasqualini,
Salman Habib,
Katrin Heitmann,
David Higdon,
Paul A. Johnson
Abstract:
The transition from linear to nonlinear dynamical elasticity in rocks is of considerable interest in seismic wave propagation as well as in understanding the basic dynamical processes in consolidated granular materials. We have carried out a careful experimental investigation of this transition for Berea and Fontainebleau sandstones. Below a well-characterized strain, the materials behave linear…
▽ More
The transition from linear to nonlinear dynamical elasticity in rocks is of considerable interest in seismic wave propagation as well as in understanding the basic dynamical processes in consolidated granular materials. We have carried out a careful experimental investigation of this transition for Berea and Fontainebleau sandstones. Below a well-characterized strain, the materials behave linearly, transitioning beyond that point to a nonlinear behavior which can be accurately captured by a simple macroscopic dynamical model. At even higher strains, effects due to a driven nonequilibrium state, and relaxation from it, complicate the characterization of the nonlinear behavior.
△ Less
Submitted 27 February, 2004;
originally announced March 2004.
-
The Inflationary Perturbation Spectrum
Authors:
Salman Habib,
Katrin Heitmann,
Gerard Jungman,
Carmen Molina-Paris
Abstract:
Motivated by the prospect of testing inflation from precision cosmic microwave background observations, we present analytic results for scalar and tensor perturbations in single-field inflation models based on the application of uniform approximations. This technique is systematically improvable, possesses controlled error bounds, and does not rely on assuming the slow-roll parameters to be cons…
▽ More
Motivated by the prospect of testing inflation from precision cosmic microwave background observations, we present analytic results for scalar and tensor perturbations in single-field inflation models based on the application of uniform approximations. This technique is systematically improvable, possesses controlled error bounds, and does not rely on assuming the slow-roll parameters to be constant. We provide closed-form expressions for the power spectra and the corresponding scalar and tensor spectral indices.
△ Less
Submitted 23 August, 2002;
originally announced August 2002.
-
Out of Equilibrium Dynamics of Supersymmetry at High Energy Density
Authors:
J. Baacke,
D. Cormier,
H. J. de Vega,
K. Heitmann
Abstract:
We investigate the out of equilibrium dynamics of global chiral supersymmetry at finite energy density. We concentrate on two specific models. The first is the massive Wess-Zumino model which we study in a selfconsistent one-loop approximation. We find that for energy densities above a certain threshold, the fields are driven dynamically to a point in field space at which the fermionic component…
▽ More
We investigate the out of equilibrium dynamics of global chiral supersymmetry at finite energy density. We concentrate on two specific models. The first is the massive Wess-Zumino model which we study in a selfconsistent one-loop approximation. We find that for energy densities above a certain threshold, the fields are driven dynamically to a point in field space at which the fermionic component of the superfield is massless. The state, however is found to be unstable, indicating a breakdown of the one-loop approximation. To investigate further, we consider an O(N) massive chiral model which is solved exactly in the large $N$ limit. For sufficiently high energy densities, we find that for late times the fields reach a nonperturbative minimum of the effective potential degenerate with the perturbative minimum. This minimum is a true attractor for O(N) invariant states at high energy densities, and this provides a mechanism for determining which of the otherwise degenerate vacua is chosen by the dynamics. The final state for large energy density is a cloud of massless particles (both bosons and fermions) around this new nonperturbative supersymmetric minimum. By introducing boson masses which softly break the supersymmetry, we demonstrate a see-saw mechanism for generating small fermion masses. We discuss some of the cosmological implications of our results.
△ Less
Submitted 15 October, 2001;
originally announced October 2001.
-
Dynamics of coupled bosonic systems with applications to preheating
Authors:
Daniel Cormier,
Katrin Heitmann,
Anupam Mazumdar
Abstract:
Coupled, multi-field models of inflation can provide several attractive features unavailable in the case of a single inflaton field. These models have a rich dynamical structure resulting from the interaction of the fields and their associated fluctuations. We present a formalism to study the nonequilibrium dynamics of coupled scalar fields. This formalism solves the problem of renormalizing int…
▽ More
Coupled, multi-field models of inflation can provide several attractive features unavailable in the case of a single inflaton field. These models have a rich dynamical structure resulting from the interaction of the fields and their associated fluctuations. We present a formalism to study the nonequilibrium dynamics of coupled scalar fields. This formalism solves the problem of renormalizing interacting models in a transparent way using dimensional regularization. The evolution is generated by a renormalized effective Lagrangian which incorporates the dynamics of the mean fields and their associated fluctuations at one-loop order. We apply our method to two problems of physical interest: (i) a simple two-field model which exemplifies applications to reheating in inflation, and (ii) a supersymmetric hybrid inflation model. This second case is interesting because inflation terminates via a smooth phase transition which gives rise to a spinodal instability in one of the fields. We study the evolution of the zero mode of the fields and the energy density transfer to the fluctuations from the mean fields. We conclude that back reaction effects can be significant over a wide parameter range. In particular for the supersymmetric hybrid model we find that particle production can be suppressed due to these effects.
△ Less
Submitted 23 March, 2002; v1 submitted 22 May, 2001;
originally announced May 2001.
-
Gauge Fields Out-Of-Equilibrium: A Gauge Invariant Formulation and the Coulomb Gauge
Authors:
Katrin Heitmann
Abstract:
We study the abelian Higgs model out-of-equilibrium in two different approaches, a gauge invariant formulation, proposed by Boyanovsky et al. \cite{Boyanovsky:1996dc} and in the Coulomb gauge. We show that both approaches become equivalent in a consistent one loop approximation. Furthermore, we carry out a proper renormalization for the model in order to prepare the equations for a numerical imp…
▽ More
We study the abelian Higgs model out-of-equilibrium in two different approaches, a gauge invariant formulation, proposed by Boyanovsky et al. \cite{Boyanovsky:1996dc} and in the Coulomb gauge. We show that both approaches become equivalent in a consistent one loop approximation. Furthermore, we carry out a proper renormalization for the model in order to prepare the equations for a numerical implementation. The additional degrees of freedom, which arise in gauge theories, influence the behavior of the system dramatically. A comparison with results in the 't Hooft-Feynman background gauge found by us recently, shows very good agreement.
△ Less
Submitted 24 January, 2001;
originally announced January 2001.
-
Dynamics of O(N) chiral supersymmetry at finite energy density
Authors:
J. Baacke,
D. Cormier,
H. J. de Vega,
K. Heitmann
Abstract:
We consider an O(N) version of a massive, interacting, chiral supersymmetry model solved exactly in the large N limit. We demonstrate that the system approaches a stable attractor at high energy densities, corresponding to a non-perturbative state for which the relevant field quanta are massless. The state is one of spontaneously broken O(N), which, due to the influence of supersymmetry, does no…
▽ More
We consider an O(N) version of a massive, interacting, chiral supersymmetry model solved exactly in the large N limit. We demonstrate that the system approaches a stable attractor at high energy densities, corresponding to a non-perturbative state for which the relevant field quanta are massless. The state is one of spontaneously broken O(N), which, due to the influence of supersymmetry, does not become restored at high energies. Introducing soft supersymmetry breaking to the Lagrangian results in scalar masses at the soft breaking scale m_s independent of the mass scale of supersymmetry mu, with even smaller masses for the fermions.
△ Less
Submitted 10 July, 2001; v1 submitted 30 November, 2000;
originally announced November 2000.
-
Nonequilibrium evolution and symmetry structure of the large-N $Φ^4$ model at finite temperature
Authors:
J. Baacke,
K. Heitmann
Abstract:
We consider the large-N $Φ^4$ theory with spontaneously broken symmetry at finite temperature. We study, in the large-N limit, quantum states which are characterized by a time dependent, spatially homogenous expectation value of one of the field components, $φ_N(t)$, and by quantum fluctuations of the other $N-1$ components, that evolve in the background of the classical field. Investigating suc…
▽ More
We consider the large-N $Φ^4$ theory with spontaneously broken symmetry at finite temperature. We study, in the large-N limit, quantum states which are characterized by a time dependent, spatially homogenous expectation value of one of the field components, $φ_N(t)$, and by quantum fluctuations of the other $N-1$ components, that evolve in the background of the classical field. Investigating such systems out of equilibrium has recently been shown to display several interesting features. We extend here this type of investigations to finite temperature systems. Essentially the novel features observed at T=0 carry over to finite temperature. This is not unexpected, as the main mechanisms that determine the late-time behavior remain the same. We extend two empirical - presumably exact - relations for the late-time behavior to finite temperature and use them to define the boundaries between the region of different asymptotic regimes. This results in a phase diagram with the temperature and the initial value of the classical field as parameters, the phases being characterized by spontaneous symmetry breaking resp. symmetry restoration. The time evolution is computed numerically and agrees very well with the expectations.
△ Less
Submitted 18 April, 2000; v1 submitted 31 March, 2000;
originally announced March 2000.
-
Gauge invariance of the one-loop effective action of the Higgs field in the SU(2) Higgs model
Authors:
J. Baacke,
K. Heitmann
Abstract:
The one-loop effective action of the abelian and nonabelian Higgs models has been studied in various gauges, in the context of instanton and sphaleron transition, bubble nucleation and most recently in nonequilibrium dynamics. Gauge invariance is expected on account of Nielsen' s theorem, if the classical background field is an extremum of the classical action, i.e., a solution of the classical…
▽ More
The one-loop effective action of the abelian and nonabelian Higgs models has been studied in various gauges, in the context of instanton and sphaleron transition, bubble nucleation and most recently in nonequilibrium dynamics. Gauge invariance is expected on account of Nielsen' s theorem, if the classical background field is an extremum of the classical action, i.e., a solution of the classical equation of motion. We substantiate this general statement for the one-loop effective action, as computed using mode functions. We show that in the gauge-Higgs sector there are two types of modes that satisfy the same equation of motion as the Faddeev-Popov modes. We apply the general analysis to the computation of the fluctuation determinant for bubble nucleation in the SU(2) Higgs model in the 't Hooft gauge with general gauge parameter $ξ$. We show that due to the cancellation of the modes mentioned above the fluctuation determinant is independent of $ξ$.
△ Less
Submitted 29 July, 1999; v1 submitted 27 May, 1999;
originally announced May 1999.
-
Nonequilibrium dynamics of fermions in a spatially homogeneous scalar background field
Authors:
Juergen Baacke,
Katrin Heitmann,
Carsten Patzold
Abstract:
We consider the time evolution of systems in which a spatially homogeneous scalar field is coupled to fermions. The quantum back-reaction is taken into account in one-loop approximation. We set up the basic equations and their renormalization in a form suitable for numerical computations. The initial singularities appearing in the renormalized equations are removed by a Bogoliubov transformation…
▽ More
We consider the time evolution of systems in which a spatially homogeneous scalar field is coupled to fermions. The quantum back-reaction is taken into account in one-loop approximation. We set up the basic equations and their renormalization in a form suitable for numerical computations. The initial singularities appearing in the renormalized equations are removed by a Bogoliubov transformation. The equations are then generalized to those in a spatially flat Friedmann-Robertson-Walker universe.
We have implemented the Minkowski space equations numerically and present results for the time evolution with various parameter sets. We find that fermion fluctuations are not in general as ineffective as assumed previously but show interesting features which should be studied further. In an especially interesting example we find that fermionic fluctuations can ``catalyze'' the evolution of bosonic fluctuations.
△ Less
Submitted 30 June, 1998; v1 submitted 30 May, 1998;
originally announced June 1998.
-
Renormalization of nonequilibrium dynamics at large N and finite temperature
Authors:
Juergen Baacke,
Katrin Heitmann,
Carsten Patzold
Abstract:
We generalize a previously proposed renormalization and computation scheme for nonequilibrium dynamics to include finite temperature and one-loop selfconsistency as arising in the large-N limit. Since such a scheme amounts essentially to tadpole summation, it also includes, at high temperature, the hard mass corrections proportional to T^2. We present some numerical examples at T=0 and at finite…
▽ More
We generalize a previously proposed renormalization and computation scheme for nonequilibrium dynamics to include finite temperature and one-loop selfconsistency as arising in the large-N limit. Since such a scheme amounts essentially to tadpole summation, it also includes, at high temperature, the hard mass corrections proportional to T^2. We present some numerical examples at T=0 and at finite temperature; the results reproduce the essential features of other groups. Especially, we can confirm a recently discovered sum rule for the late time behaviour.
△ Less
Submitted 23 December, 1997;
originally announced December 1997.
-
On the choice of initial states in nonequilibrium dynamics
Authors:
Juergen Baacke,
Katrin Heitmann,
Carsten Patzold
Abstract:
Imposing initial conditions to nonequilibrium systems at some time t_0 leads, in renormalized quantum field theory, to the appearance of singularities in the variable t-t_0 in relevant physical quantities, such as energy density and pressure. The origin of these `initial singularities' can be traced back to the choice of initial state. We construct here, by a Bogoliubov transformation, initial s…
▽ More
Imposing initial conditions to nonequilibrium systems at some time t_0 leads, in renormalized quantum field theory, to the appearance of singularities in the variable t-t_0 in relevant physical quantities, such as energy density and pressure. The origin of these `initial singularities' can be traced back to the choice of initial state. We construct here, by a Bogoliubov transformation, initial states such that these singularities are eliminated. While the construction is not unique it can be considered a minimal way of taking into account the nonequilibrium evolution of the system prior to t_0.
△ Less
Submitted 19 November, 1997;
originally announced November 1997.
-
Renormalization of nonequilibrium dynamics in FRW cosmology
Authors:
Jurgen Baacke,
Katrin Heitmann,
Carsten Patzold
Abstract:
We derive the renormalized nonequilibrium equations of motion for a scalar field and its quantum back reaction in a conformally flat Friedmann-Robertson-Walker universe. We use a fully covariant formalism proposed by us recently for handling numerically and analytically nonequilibrium dynamics in one-loop approximation. The system is assumed to be in a conformal vacuum state initially. We use di…
▽ More
We derive the renormalized nonequilibrium equations of motion for a scalar field and its quantum back reaction in a conformally flat Friedmann-Robertson-Walker universe. We use a fully covariant formalism proposed by us recently for handling numerically and analytically nonequilibrium dynamics in one-loop approximation. The system is assumed to be in a conformal vacuum state initially. We use dimensional regularization; we find that the counter terms can be chosen independent of the initial conditions though the divergent leading order graphs do depend on them.
△ Less
Submitted 28 July, 1997; v1 submitted 6 June, 1997;
originally announced June 1997.
-
Nonequilibrium dynamics: preheating in the SU(2) Higgs model
Authors:
Juergen Baacke,
Katrin Heitmann,
Carsten Paetzold
Abstract:
The term `preheating' has been introduced recently to denote the process in which energy is transferred from a classical inflaton field into fluctuating field (particle) degrees of freedom without generating yet a real thermal ensemble. The models considered up to now include, besides the inflaton field, scalar or fermionic fluctuations. On the other hand the typical ingredient of an inflationar…
▽ More
The term `preheating' has been introduced recently to denote the process in which energy is transferred from a classical inflaton field into fluctuating field (particle) degrees of freedom without generating yet a real thermal ensemble. The models considered up to now include, besides the inflaton field, scalar or fermionic fluctuations. On the other hand the typical ingredient of an inflationary scenario is a nonabelian spontaneously broken gauge theory. So the formalism should also be developed to include gauge field fluctuations excited by the inflaton or Higgs field. We have chosen here, as the simplest nonabelian example, the SU(2) Higgs model. We consider the model at temperature zero. From the technical point of view we generalize an analytical and numerical renormalized formalism developed by us recently to coupled channnel systems. We use the 't Hooft-Feynman gauge and dimensional regularization. We present some numerical results but reserve a more exhaustive discussion of solutions within the paramter space of two couplings and the initial value of the Higgs field to a future publication.
△ Less
Submitted 6 December, 1996; v1 submitted 5 December, 1996;
originally announced December 1996.
-
Nonequilibrium dynamics: a renormalized computation scheme
Authors:
Juergen Baacke,
Katrin Heitmann,
Carsten Patzold
Abstract:
We present a regularized and renormalized version of the one-loop nonlinear relaxation equations that determine the non-equilibrium time evolution of a classical (constant) field coupled to its quantum fluctuations. We obtain a computational method in which the evaluation of divergent fluctuation integrals and the evaluation of the exact finite parts are cleanly separated so as to allow for a wi…
▽ More
We present a regularized and renormalized version of the one-loop nonlinear relaxation equations that determine the non-equilibrium time evolution of a classical (constant) field coupled to its quantum fluctuations. We obtain a computational method in which the evaluation of divergent fluctuation integrals and the evaluation of the exact finite parts are cleanly separated so as to allow for a wide freedom in the choice of regularization and renormalization schemes. We use dimensional regularization here. Within the same formalism we analyze also the regularization and renormalization of the energy-momentum tensor. The energy density serves to monitor the reliability of our numerical computation. The method is applied to the simple case of a scalar phi^4 theory; the results are similar to the ones found previously by other groups.
△ Less
Submitted 28 May, 1998; v1 submitted 1 August, 1996;
originally announced August 1996.