Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cosmological perturbation theory”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

Testing higher-order Lagrangian perturbation theory against numerical simulation. 1: Pancake models

We present results showing an improvement of the accuracy of perturbation theory as applied to cosmological structure formation for a useful range of quasi-linear scales. The Lagrangian theory of gravitational instability of an Einstein-de Sitter dust cosmogony investigated and solved up to the third order is compared with numerical simulations. In this paper we study the dynamics of pancake models as a first step. In previous work the accuracy of several analytical approximations for the modeling of large-scale structure in the mildly non-linear regime was analyzed in the same way, allowing for direct comparison of the accuracy of various approximations. In particular, the Zel'dovich approximation (hereafter ZA) as a subclass of the first-order Lagrangian perturbation solutions was found to provide an excellent approximation to the density field in the mildly non-linear regime (i.e. up to a linear r.m.s. density contrast of sigma is approximately 2). The performance of ZA in hierarchical clustering models can be greatly improved by truncating the initial power spectrum (smoothing the initial data). We here explore whether this approximation can be further improved with higher-order corrections in the displacement mapping from homogeneity. We study a single pancake model (truncated power-spectrum with power-spectrum with power-index n = -1) using cross-correlation statistics employed in previous work. We found that for all statistical methods used the higher-order corrections improve the results obtained for the first-order solution up to the stage when sigma (linear theory) is approximately 1. While this improvement can be seen for all spatial scales, later stages retain this feature only above a certain scale which is increasing with time. However, third-order is not much improvement over second-order at any stage. The total breakdown of the perturbation approach is observed at the stage, where sigma (linear theory) is approximately 2, which corresponds to the onset of hierarchical clustering. This success is found at a considerable higher non-linearity than is usual for perturbation theory. Whether a truncation of the initial power-spectrum in hierarchical models retains this improvement will be analyzed in a forthcoming work.

Buchert, T.↗

Testing higher-order Lagrangian perturbation theory against numerical simulations. 2: Hierarchical models

We present results showing an improvement of the accuracy of perturbation theory as applied to cosmological structure formation for a useful range of scales. The Lagrangian theory of gravitational instability of Friedmann-Lemaitre cosmogonies is compared with numerical simulations. We study the dynamics of hierarchical models as a second step. In the first step we analyzed the performance of the Lagrangian schemes for pancake models, the difference being that in the latter models the initial power spectrum is truncated. This work probed the quasi-linear and weakly non-linear regimes. We here explore whether the results found for pancake models carry over to hierarchical models which are evolved deeply into the non-linear regime. We smooth the initial data by using a variety of filter types and filter scales in order to determine the optimal performance of the analytical models, as has been done for the 'Zel'dovich-approximation' - hereafter TZA - in previous work. We find that for spectra with negative power-index the second-order scheme performs considerably better than TZA in terms of statistics which probe the dynamics, and slightly better in terms of low-order statistics like the power-spectrum. However, in contrast to the results found for pancake models, where the higher-order schemes get worse than TZA at late non-linear stages and on small scales, we here find that the second-order model is as robust as TZA, retaining the improvement at later stages and on smaller scales. In view of these results we expect that the second-order truncated Lagrangian model is especially useful for the modelling of standard dark matter models such as Hot-, Cold-, and Mixed-Dark-Matter.

Melott, A. L.↗

A photoionization instability in the early intergalactic medium

It is argued that any fairly uniform source of ionizing photons can be the cause of an instability in the pregalactic medium on scales larger than a photon path length. Underdense regions receive more ionizing energy per atom and reach higher temperature and entropy, driving the density down still further. Fluctuations created by this instability can lead to the formation of structures resembling protogalaxies and intergalactic clouds, obviating the need for gas clouds or density perturbations of earlier cosmological provenance, as is usually assumed in theories of galaxy and structure formation. Characteristic masses for clouds produced by the instability, with log mass in solar units plotted against log radius in kpc, are illustrated.

Hogan, Craig J.↗

Microwave anisotropies in the light of the data from the COBE satellite

The recent measurement of anisotropies in the cosmic microwave background by the Cosmic Background Explorer (COBE) satellite and the recent South Pole experiment offer an excellent opportunity to probe cosmological theories. We test a class of theories in which the universe today is flat and matter dominated, and primordial perturbations are adiabatic parameterized by an index n. In this class of theories the predicted signal in the South Pole experiment depends on n, the Hubble constant, and the baryon density. For n = 1 a large region of this parameter space is ruled out, but there is still a window open which satisfies constraints from COBE, the South Pole experiment, and big bang nucleosynthesis.

Dodelson, Scott↗

Origin of density fluctuations in extended inflation

The density fluctuations (both curvature and isocurvature) that arise due to quantum fluctuations in a simple model of extended inflation based upon the Jordan-Brans-Dicke theory are calculated. Curvature fluctuations arise due to quantum fluctuations in the Brans-Dicke field, in general have a nonscale-invariant spectrum, and can have an amplitude that is cosmologically acceptable and interesting without having to tune any coupling constant to a very small value. The density perturbations that arise due to the inflation field are subdominant. If there are other massless fields in the theory, e.g., an axion or an ilion, then isocurvature fluctuations arise in these fields too. Production of gravitational waves and the massless particles associated with excitations of the Brans-Dicke field are also discussed. Several attempts at more realistic models of extended inflation are also analyzed. The importance of the Einstein conformal frame in calculating curvature fluctuations is emphasized. When viewed in this frame, extended inflation closely resembles slow-rollover inflation with an exponential potential and the usual formula for the amplitude of curvature perturbations applies.

Kolb, Edward W.↗

Experiment-specific cosmic microwave background calculations made easier - Approximation formula for smoothed delta T/T windows

Simple and easy to implement elementary function approximations are introduced to the spectral window functions needed in calculations of model predictions of the cosmic microwave backgrond (CMB) anisotropy. These approximations allow the investigator to obtain model delta T/T predictions in terms of single integrals over the power spectrum of cosmological perturbations and to avoid the necessity of performing the additional integrations. The high accuracy of these approximations is demonstrated here for the CDM theory-based calculations of the expected delta T/T signal in several experiments searching for the CMB anisotropy.

Gorski, Krzysztof M.↗

The Mach number of the cosmic flow - A critical test for current theories

A new cosmological, self-contained test using the ratio of mean velocity and the velocity dispersion in the mean flow frame of a group of test objects is presented. To allow comparison with linear theory, the velocity field must first be smoothed on a suitable scale. In the context of linear perturbation theory, the Mach number M(R) which measures the ratio of power on scales larger than to scales smaller than the patch size R, is independent of the perturbation amplitude and also of bias. An apparent inconsistency is found for standard values of power-law index n = 1 and cosmological density parameter Omega = 1, when comparing values of M(R) predicted by popular models with tentative available observations. Nonstandard models based on adiabatic perturbations with either negative n or small Omega value also fail, due to creation of unacceptably large microwave background fluctuations.

Ostriker, Jeremiah P.↗

Theoretical Astrophysics at Fermilab

The Theoretical Astrophysics Group works on a broad range of topics ranging from string theory to data analysis in the Sloan Digital Sky Survey. The group is motivated by the belief that a deep understanding of fundamental physics is necessary to explain a wide variety of phenomena in the universe. During the three years 2001-2003 of our previous NASA grant, over 120 papers were written; ten of our postdocs went on to faculty positions; and we hosted or organized many workshops and conferences. Kolb and collaborators focused on the early universe, in particular and models and ramifications of the theory of inflation. They also studied models with extra dimensions, new types of dark matter, and the second order effects of super-horizon perturbations. S tebbins, Frieman, Hui, and Dodelson worked on phenomenological cosmology, extracting cosmological constraints from surveys such as the Sloan Digital Sky Survey. They also worked on theoretical topics such as weak lensing, reionization, and dark energy. This work has proved important to a number of experimental groups [including those at Fermilab] planning future observations. In general, the work of the Theoretical Astrophysics Group has served as a catalyst for experimental projects at Fennilab. An example of this is the Joint Dark Energy Mission. Fennilab is now a member of SNAP, and much of the work done here is by people formerly working on the accelerator. We have created an environment where many of these people made transition from physics to astronomy. We also worked on many other topics related to NASA s focus: cosmic rays, dark matter, the Sunyaev-Zel dovich effect, the galaxy distribution in the universe, and the Lyman alpha forest. The group organized and hosted a number of conferences and workshop over the years covered by the grant. Among them were:

Source record↗

The galaxy velocity field and CDM models

It is generally accepted that some kind of non-baryonic dark matter accounts for most of the mass density of the universe. Considering such a component has become, in the last decade, a key ingredient in current theories of structure formation. In particular, the Cold Dark Matter (CDM) scenario has proven to be quite successful in explaining most of the observed properties of galaxies and of their large-scale distribution. The standard CDM model is characterized by a primordial Zel'dovich spectrum, of random-phase adiabatic perturbations in a universe with density parameter omega sub 0 = 1 and vanishing cosmological constant. This poster paper presents an analysis of observational data on peculiar motion of optical galaxies in comparison to the predictions of CDM models where the assumptions of the standard scenario: omega sub 0 = 1, n = 1, and bias parameter b = 1 are relaxed. In particular, CDM models with 0 less than n less than 1 and 0.4 less than omega sub 0 less than 1 are considered.

Tormen, Giuseppe↗

The minimal power spectrum: Higher order contributions

It has been an accepted belief for some time that gravity induces a minimal tail P(k) approximately k(exp 4) in the power spectrum as k approaches 0 for distributions with no initial power on large scales. In a recent numerical experiment with initial power confined to a restricted range in k, Shandarin and Melott (1990) found a k approaches 0 tail that at early stages of evolution behaves as k(exp 4) and grows with time as a(exp 4)(t), where a(t) is the cosmological expansion factor, and at late times depends on scale as k(exp 3) and grows with time as a(exp 2)(t). I compute analytically several contributions to the power spectrum of higher order than those included in earlier work, and I apply the results to the particular case of initial power restricted to a finite range of k. As expected, in the perturbative regime P(k) approximately a(exp 4)k(exp 4) from the first correction to linear perturbation theory is the dominant term as k approaches 0. Numerical investigations show that the higher order contributions go as k(exp 4) also. However, perturbation theory alone cannot tell whether the P approximately a(exp 2)k(exp 3) result is 'nonperturbative' or a numerical artifact.

Fry, J. N.↗

A Cosmic Microscope for the Preheating Era

Light fields with spatially varying backgrounds can modulate cosmic preheating, and imprint the nonlinear effects of preheating dynamics at tiny scales on large scale fluctuations. This provides us a unique probe into the preheating era which we dub the “cosmic microscope”. We identify a distinctive effect of preheating on scalar perturbations that turns the Gaussian primordial fluctuations of a light scalar field into square waves, like a diode. The effect manifests itself as local non-Gaussianity. We present a model, “modulated partial preheating”, where this nonlinear effect is consistent with current observations and can be reached by near future cosmic probes.

Cosmology of Theories beyond the SM↗

Comparison of dynamical approximation schemes for nonlinear gravitaional clustering

We have recently conducted a controlled comparison of a number of approximations for gravitational clustering against the same n-body simulations. These include ordinary linear perturbation theory (Eulerian), the lognormal approximation, the adhesion approximation, the frozen-flow approximation, the Zel'dovich approximation (describable as first-order Lagrangian perturbation theory), and its second-order generalization. In the last two cases we also created new versions of the approximation by truncation, i.e., by smoothing the initial conditions with various smoothing window shapes and varying their sizes. The primary tool for comparing simulations to approximation schemes was cross-correlation of the evolved mass density fields, testing the extent to which mass was moved to the right place. The Zel'dovich approximation, with initial convolution with a Gaussian e(exp -k(exp 2)/k(sub G(exp 2)), where k(sub G) is adjusted to be just into the nonlinear regime of the evolved model (details in text) worked extremely well. Its second-order generalization worked slightly better. We recommend either n-body simulations or our modified versions of the Zel'dovich approximation, depending upon the purpose. The theoretical implication is that pancaking is implicit in all cosmological gravitational clustering, at least from Gaussian initial conditions, even when subcondensations are present. This in turn provides a natural explanation for the presence of sheets and filaments in the observed galaxy distribution. Use of the approximation scheme can permit extremely rapid generation of large numbers of realizations of model universes with good accuracy down to galaxy group mass scales.

Melott, Adrian L.↗

Spillway Preheating

In traditional models only an order one fraction of energy is transferred from the inflaton to radiation through nonperturbative resonance production in preheating immediately after inflation, due to backreaction effects. We propose a particle production mechanism that could improve the depletion of the inflaton energy density by up to four orders of magnitude. The improvement comes from the fast perturbative decays of resonantly produced daughter particles. They act as a “spillway” to drain these daughter particles, reducing their backreaction on the inflaton and keeping the resonant production effective for a longer period. Thus we dub the scenario “spillway preheating”. We also show that the fraction of energy density remaining in the inflaton has a simple inverse power-law scaling in the scenario. In general, spillway preheating is a much more efficient energy dissipation mechanism, which may have other applications in model building for particle physics.

Cosmology of Theories beyond the SM↗

Higher order moments of the matter distribution in scale-free cosmological simulations with large dynamic range

We calculate reduced moments (xi bar)(sub q) of the matter density fluctuations, up to order q = 5, from counts in cells produced by particle-mesh numerical simulations with scale-free Gaussian initial conditions. We use power-law spectra P(k) proportional to k(exp n) with indices n = -3, -2, -1, 0, 1. Due to the supposed absence of characteristic times or scales in our models, all quantities are expected to depend on a single scaling variable. For each model, the moments at all times can be expressed in terms of the variance (xi bar)(sub 2), alone. We look for agreement with the hierarchical scaling ansatz, according to which ((xi bar)(sub q)) proportional to ((xi bar)(sub 2))(exp (q - 1)). For n less than or equal to -2 models, we find strong deviations from the hierarchy, which are mostly due to the presence of boundary problems in the simulations. A small, residual signal of deviation from the hierarchical scaling is however also found in n greater than or equal to -1 models. The wide range of spectra considered and the large dynamic range, with careful checks of scaling and shot-noise effects, allows us to reliably detect evolution away from the perturbation theory result.

Lucchin, Francesco↗

Bias and high-order galaxy correlation functions in the APM galaxy survey

On large scales, the higher order moments of the mass distribution, S(sub J) = bar-zeta(sub J)/bar-zeta(sup J-1)(sub 2), e.g., the skewness S(sub 3) and kurtosis S(sub 4), can be predicted using nonlinear perturbation theory. Comparison of these predictions with moments of the observed galaxy distribution probes the bias between galaxies and mass. Applying this method to models with initially Gaussian fluctuations and power spectra P(k) similar to that of galaxies in the Automatic Plate Measuring (APM) survey, we find that the predicted higher order moments S(sub J)(R) are in good agreement with those directly inferred from the APM survey in the absence of bias. We use this result to place limits on the linear and nonlinear bias parameters. Models in which the extra power observed on large scales (with respect to the standard cold dark matter (CDM) model) is produced by scale-dependent bias match the APM higher order amplitudes only if nonlinear bias (rather than nonlinear gravity) generates the observed higher order moments. When normalized to Cosmic Background Explorer Differential Microwave Radiometer (COBE DMR), these models are siginificantly ruled out by the S(sub 3) observations. The cold plus hot dark matter model normalized to COBE can reproduce the APM higher order correlations if one introduces nonlinear bias terms, while the low-density CDM model with a cosmological constant does not require any bias to fit the large-scale amplitudes.

Gaztanaga, Enrique↗

Comparison of dynamical approximation schemes for non-linear gravitational clustering

We have recently conducted a controlled comparison of a number of approximations for gravitational clustering against the same n-body simulations. These include ordinary linear perturbation theory (Eulerian), the adhesion approximation, the frozen-flow approximation, the Zel'dovich approximation (describable as first-order Lagrangian perturbation theory), and its second-order generalization. In the last two cases we also created new versions of approximation by truncation, i.e., smoothing the initial conditions by various smoothing window shapes and varying their sizes. The primary tool for comparing simulations to approximation schemes was crosscorrelation of the evolved mass density fields, testing the extent to which mass was moved to the right place. The Zel'dovich approximation, with initial convolution with a Gaussian e(exp -k(exp 2)/k(exp 2, sub G)) where k(sub G) is adjusted to be just into the nonlinear regime of the evolved model (details in text) worked extremely well. Its second-order generalization worked slightly better. All other schemes, including those proposed as generalizations of the Zel'dovich approximation created by adding forces, were in fact generally worse by this measure. By explicitly checking, we verified that the success of our best-choice was a result of the best treatment of the phases of nonlinear Fourier components. Of all schemes tested, the adhesion approximation produced the most accurate nonlinear power spectrum and density distribution, but its phase errors suggest mass condensations were moved to slightly the wrong location. Due to its better reproduction of the mass density distribution function and power spectrum, it might be preferred for some uses. We recommend either n-body simulations or our modified versions of the Zel'dovich approximation, depending upon the purpose. The theoretical implication is that pancaking is implicit in all cosmological gravitational clustering, at least from Gaussian initial conditions, even when subcondensations are present.

Melott, Adrian L.↗

Dark Matter in the Universe and in the Galaxy

During the past four years, Prof. Kamionkowski and collaborators have made progress in research on the nature and distribution of dark-matter in the Universe and in the Galaxy, and on related topics in astrophysics and cosmology. We have made progress on research on the cosmic microwave background, large-scale structure, issues related to particle dark matter, and the gamma-ray-burst enigma. A significant fraction of the research supported by this ATP has been on the cosmic microwave background (CMB). Prof. Kamionkowski and collaborators showed how the polarization of the CMB could be used to detect long-wavelength gravitational waves, such as those produced by inflation. With Kosowsky, Prof. Kamionkowski calculated the amplitude of a stochastic gravitational-wave background that could be detected for a satellite experiment of a given sensitivity and angular resolution. They showed that polarization should improve the sensitivity oa MAP to these gravity waves, and that the Planck Surveyor should do even better. Prof. Kamionkowski, Caldwell, and a student calculated and illustrated the CMB temperature/polarization pattern produced by a single plane-wave gravitational wave. They calculated the amplitude of such a wave that would be detectable with MAP and Planck, and compared that with the sensitivity of traditional gravitational-wave detectors like LIGO and LISA. With Lue and Wang, the PI showed how parity violation from new high-energy physics could conceivably give rise to an observable signature in the CMB polarization. With Loeb, Prof. Kamionkowski showed how measurement of the polarization of CMB photons scattered by hot gas in a cluster could be used to determine the quadrupole moment of the CMB incident on that cluster. Prof. Kamionkowski and Jaffe calculated the amplitude of secondary anisotropies produced by scattering of CMB photons from reionized regions. Research has also been carried out on probing the large-scale distribution of mass in the Universe today, and on structure-formation theories. They investigated the possibility of determining the large-scale distribution of mass in the Universe via measurement of ellipticity-ellipticity correlations in the FIRST radio survey induced by weak gravitational lensing due to mass inhomogeneities along the line of sight. Dr. Summers, Prof. Kamionkowski, and a student investigated the distribution of protogalactic masses and angular momenta in an effort to understand how the luminosity function and angular-momentum distribution of disk galaxies arises from an initial power spectrum of density perturbations. Several projects related to the distribution and possible detection of dark matter in our Galactic halo were studied. The PI was involved in several projects involving the calculation of nuclear-reaction rates needed for stellar evolution and for predictions of solar-neutrino fluxes.

Kamionkowski, Marc↗