Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “N-body cosmology”

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.

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Numerical techniques for large cosmological N-body simulations

Techniques for carrying out large N-body simulations of the gravitational evolution of clustering in the fundamental cube of an infinite periodic universe are described and compared. The accuracy of the forces derived from several commonly used particle mesh schemes is examined, showing how submesh resolution can be achieved by including short-range forces between particles by direct summation techniques. The time integration of the equations of motion is discussed, and the accuracy of the codes for various choices of 'time' variable and time step is tested by considering energy conservation as well as by direct analysis of particle trajectories. Methods for generating initial particle positions and velocities corresponding to a growing mode representation of a specified power spectrum of linear density fluctuations are described. The effects of force resolution are studied and different simulation schemes are compared. An algorithm is implemented for generating initial conditions by varying the number of particles, the initial amplitude of density fluctuations, and the initial peculiar velocity field.

Efstathiou, G.↗
NASA NTRS1985

Cluster-cluster clustering

The cluster correlation function xi sub c(r) is compared with the particle correlation function, xi(r) in cosmological N-body simulations with a wide range of initial conditions. The experiments include scale-free initial conditions, pancake models with a coherence length in the initial density field, and hybrid models. Three N-body techniques and two cluster-finding algorithms are used. In scale-free models with white noise initial conditions, xi sub c and xi are essentially identical. In scale-free models with more power on large scales, it is found that the amplitude of xi sub c increases with cluster richness; in this case the clusters give a biased estimate of the particle correlations. In the pancake and hybrid models (with n = 0 or 1), xi sub c is steeper than xi, but the cluster correlation length exceeds that of the points by less than a factor of 2, independent of cluster richness. Thus the high amplitude of xi sub c found in studies of rich clusters of galaxies is inconsistent with white noise and pancake models and may indicate a primordial fluctuation spectrum with substantial power on large scales.

Barnes, J.↗
NASA NTRS1989

Recovering the full velocity and density fields from large-scale redshift-distance samples

A new method for extracting the large-scale three-dimensional velocity and mass density fields from measurements of the radial peculiar velocities is presented. Galaxies are assumed to trace the velocity field rather than the mass. The key assumption made is that the Lagrangian velocity field has negligible vorticity, as might be expected from perturbations that grew by gravitational instability. By applying the method to cosmological N-body simulations, it is demonstrated that it accurately reconstructs the velocity field. This technique promises a direct determination of the mass density field and the initial conditions for the formation of large-scale structure from galaxy peculiar velocity surveys.

Bertschinger, Edmund↗
NASA NTRS1990

Potential, velocity, and density fields from sparse and noisy redshift-distance samples - Method

A method for recovering the three-dimensional potential, velocity, and density fields from large-scale redshift-distance samples is described. Galaxies are taken as tracers of the velocity field, not of the mass. The density field and the initial conditions are calculated using an iterative procedure that applies the no-vorticity assumption at an initial time and uses the Zel'dovich approximation to relate initial and final positions of particles on a grid. The method is tested using a cosmological N-body simulation 'observed' at the positions of real galaxies in a redshift-distance sample, taking into account their distance measurement errors. Malmquist bias and other systematic and statistical errors are extensively explored using both analytical techniques and Monte Carlo simulations.

Dekel, Avishai↗
NASA NTRS1991

Primordial fluctuations and non-linear structure

Several aspects of nonlinear gravitational instability are explored using two numerical experiments, each of which employs a series of 3D cosmological N-body simulations. In the first experiment, all of the initial high-frequency power is truncated above some critical wavenumber, in the second experiment, initial Fourier components are replaced with wavenumbers greater than critical by waves taken from an independent realization with the same power spectrum but with unrelated phases. These initial conditions are evolved for different values of critical wavenumbers to see how the progressive elimination or substitution of initial high-frequency components affects the final, nonlinear structure.

Little, Blane↗
NASA NTRS1994

Power spectrum, correlation function, and tests for luminosity bias in the CfA redshift survey

We describe and apply a method for directly computing the power spectrum for the galaxy distribution in the extension of the Center for Astrophysics Redshift Survey. Tests show that our technique accurately reproduces the true power spectrum for k greater than 0.03 h Mpc(exp -1). The dense sampling and large spatial coverage of this survey allow accurate measurement of the redshift-space power spectrum on scales from 5 to approximately 200 h(exp -1) Mpc. The power spectrum has slope n approximately equal -2.1 on small scales (lambda less than or equal 25 h(exp -1) Mpc) and n approximately -1.1 on scales 30 less than lambda less than 120 h(exp -1) Mpc. On larger scales the power spectrum flattens somewhat, but we do not detect a turnover. Comparison with N-body simulations of cosmological models shows that an unbiased, open universe CDM model (OMEGA h = 0.2) and a nonzero cosmological constant (CDM) model (OMEGA h = 0.24, lambda(sub zero) = 0.6, b = 1.3) match the CfA power spectrum over the wavelength range we explore. The standard biased CDM model (OMEGA h = 0.5, b = 1.5) fails (99% significance level) because it has insufficient power on scales lambda greater than 30 h(exp -1) Mpc. Biased CDM with a normalization that matches the Cosmic Microwave Background (CMB) anisotropy (OMEGA h = 0.5, b = 1.4, sigma(sub 8) (mass) = 1) has too much power on small scales to match the observed galaxy power spectrum. This model with b = 1 matches both Cosmic Background Explorer Satellite (COBE) and the small-scale power spect rum but has insufficient power on scales lambda approximately 100 h(exp -1) Mpc. We derive a formula for the effect of small-scale peculiar velocities on the power spectrum and combine this formula with the linear-regime amplification described by Kaiser to compute an estimate of the real-space power spectrum. Two tests reveal luminosity bias in the galaxy distribution: First, the amplitude of the pwer spectrum is approximately 40% larger for the brightest 50% of galaxies in volume-limited samples that have M(sub lim) greater than M*. This bias in the power spectrum is independent of scale, consistent with the peaks-bias paradigm for galaxy formation. Second, the distribution of local density around galaxies shows that regions of moderate and high density contain both very bright (M less than M* = -19.2 + 5 log h) and fainter galaxies, but that voids preferentially harbor fainter galaxies (approximately 2 sigma significance level).

Park, Changbom↗
NASA NTRS1986

Correlation of QSO absorption lines in universes dominated by cold dark matter

Theoretical predictions for the redshift correlations between QSO absorption-line systems are investigated in the context of 'cold dark matter' cosmological models. Particles in 'particle-mesh' N-body simulations are interpreted as absorbing clouds at epochs corresponding to mean redshifts, z, of 0.0, 1.25, and 3.0. The velocity correlation function for absorbing clouds is found by passing lines-of-sight through the systems and computing velocity differences for those particles which lie close to the lines. It depends strongly on z and Omega but only weakly, if at all, on the number density, diameter or mass of the clouds. Two interpretations are possible: (1) the heavy element absorption systems are associated with galaxies which are an unbiased sample of the mass distribution in an Omega(0) = 0.2 universe or (2) the Lyman-alpha absorbers are an unbiased sample of the mass in an Omega(0) = 1 universe and the heavy-element absorption systems, like galaxies, are more strongly clustered than the mass.

Salmon, J.↗
NASA NTRS1993

The effect of tidal fields on the shapes and kinematics of dark halos

We have carried out a series of N-body simulations to investigate the effect of tidal shear on the structure and kinematics of dark halos. We simulate the collapse of density perturbations using a tree code as described in Dubinski & Carlberg (1991). Density peaks are selected from a random realization of a CDM density field and used as the initial conditions for N-body simulations. We use an experimental approach to examine the effects of tidal shear on collapse. The cosmological tidal field is treated as an external time dependent potential whose strength and orientation can be varied freely. We examine the effects of the tidal field with two experiments. In the first experiment, we simulate a sample of 14 dark halos from the collapse of density peaks in the presence of a 1(sigma) tidal field. In the second experiment, we use the same initial conditions though the tidal field is turned off allowing an experimental control for comparison to highlight the influence of tidal shear on the development of the structure and kinematics of the dark halos.

Dubinski, John↗
NASA NTRS1993

Testing the frozen flow approximation

We investigate the accuracy of the frozen-flow approximation (FFA), recently proposed by Matarrese, et al. (1992), for following the nonlinear evolution of cosmological density fluctuations under gravitational instability. We compare a number of statistics between results of the FFA and n-body simulations, including those used by Melott, Pellman & Shandarin (1993) to test the Zel'dovich approximation. The FFA performs reasonably well in a statistical sense, e.g. in reproducing the counts-in-cell distribution, at small scales, but it does poorly in the crosscorrelation with n-body which means it is generally not moving mass to the right place, especially in models with high small-scale power.

Lucchin, Francesco↗
NASA NTRS1989

On the generation of a bubbly universe - A quantitative assessment of the CfA slice

A first attempt is made to calculate the properties of the matter distribution in a universe filled with overlapping bubbles produced by multiple explosions. Each spherical shell follows the cosmological Sedov-Taylor solution until it encounters another shell. Thereafter, mergers are allowed to occur in pairs on the basis of N-body results. At the final epoch, the matrix of overlapping shells is populated with 'galaxies' and the properties of slices through the numerically constructed cube compare well with CfA survey results for specified initial conditions. A statistic is found which measures the distance distribution from uniformly distributed points to the nearest galaxies on the projected plane which appears to provide a good measure of the bubbly character of the galaxy distribution. In a quantitative analysis of the CfA 'slice of the universe', a very good match is found between simulation and the real data for final average bubble radii of (13.5 + or - 1.5)/h Mpc with formal filling factor 1.0-1.5 or actual filling factor of 65-80 percent.

Ostriker, J. P.↗
NASA NTRS1985

On percolation as a cosmological test

Difficulties in the use of percolation as a complementary statistic for the galaxy clustering pattern are pointed out by studying simple toy models and dynamical N-body models that represent the competing clustering scenarios. The percolation properties are found not to be very sensitive to the presence of pancakes and strings once they are clumpy, and hence they do not distinguish properly between models that are very different. In the case of very smooth pancakes, the ability to percolate depends on sampling parameters, such as the mean number density and the volume, in a way which is unknown a priori because it depends on the same properties that the test ought to measure. This problem could, in principle, be eased by using volume-limited samples of high mean number density (an order of magnitude denser than the CfA redshift survey volume limited at 4000 km/s) and by comparing to models of identical number density and volume. An alternative approach, based on the sampling effects themselves, may provide a qualitative test for pancakes in samples of lower densities.

Dekel, A.↗
NASA NTRS1994

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.↗
NASA NTRS1991

Constraints on cosmological models from the IRAS QDOT survey results

The results of the recent IRAS survey are used to explore constraints on the nature of the Galaxy density field. Realistic simulations of the survey made from N-body simulations of the n = -1 power-law power spectrum and the Omega = 1 CDM models show that the power-law model is basically consistent with all QDOT variance results, while the CDM model has a smaller variance than the observation at the 20/h Mpc Gaussian smoothing scale. It is concluded that the QDOT results are best explained if the power spectrum of the Galaxy density field has a slope of about -1 near a wavelength of about 100/h Mpc, while the CDM model has a slope of about 0 at that scale. The possibility of marginal non-Gaussian behavior at the 20/h Mpc scale reported by Saunders et al. (1990) is examined for gravitational instability models with Gaussian initial conditions. It is argued that the observed (delta exp 3) versus (delta exp 2) in redshift space is not what these models would predict.

Park, Changbom↗
NASA NTRS1994

Optical signatures of high-redshift galaxy clusters

We combine an N-body and gasdynamic simulation of structure formation with an updated population synthesis code to explore the expected optical characteristics of a high-redshift cluster of galaxies. We examine a poor (2 keV) cluster formed in a biased, cold dark matter cosmology and employ simple, but plausible, threshold criteria to convert gas into stars. At z = 2, the forming cluster appears as a linear chain of very blue (g-r approximately equals 0) galaxies, with 15 objects brighter than r = 25 within a 1 square arcmin field of view. After 2 Gyr of evolution, the cluster viewed at z = 1 displays both freshly infalling blue galaxies and red galaxies robbed of recent accretion by interaction with the hot intracluster medium. The range in G-R colors is approximately 3 mag at z = 1, with the reddest objects lying at sites of highest galaxy density. We suggest that red, high-redshift galaxies lie in the cores of forming clusters and that their existence indicates the presence of a hot intracluster medium at redshifts z approximately equals 2. The simulated cluster viewed at z = 2 has several characteristics similar to the collection of faint, blue objects identified by Dressler et al. in a deep Hubble Space Telescope observation. The similarities provide some support for the interpretation of this collection as a high-redshift cluster of galaxies.

Evrard, August E.↗
NASA NTRS1991

Voids in a neutrino-dominated universe

In a neutrino-dominated universe, galaxies are expected to form only in large-scale sheets and filaments. Most of space should be filled by low-density regions devoid of galaxies. In this paper, N-body simulations are used to estimate the size distribution for these regions for quantitative comparison with the observed voids in recent red-shift surveys. The theoretical distribution depends very weakly on the mode or epoch of galaxy formation. With very conservative assumptions, at best marginal consistency is found even for cosmological parameters as extreme as Omega(v) = 1 and H(0) = 100 km/s per/Mpc. Any significant reduction in either H(0) or Omega(v) leads to predicted void sizes much larger than those observed. This difficulty arises because the observed voids are rarely completely empty.

Zeng, Ning↗
NASA NTRS1993

The peculiar velocities of rich clusters in the hot and cold dark matter scenarios

We present the results of a study of the peculiar velocities of rich clusters of galaxies. The peculiar motion of rich clusters in various cosmological scenarios is of interest for a number of reasons. Observationally, one can measure the peculiar motion of clusters to greater distances than galaxies because cluster peculiar motions can be determined to greater accuracy. One can also test the slope of distance indicator relations using clusters to see if galaxy properties vary with environment. We have used N-body simulations to measure the amplitude and rms cluster peculiar velocity as a function of bias parameter in the hot and cold dark matter scenarios. In addition to measuring the mean and rms peculiar velocity of clusters in the two models, we determined whether the peculiar velocity vector of a given cluster is well aligned with the gravity vector due to all the particles in the simulation and the gravity vector due to the particles present only in the clusters. We have investigated the peculiar velocities of rich clusters of galaxies in the cold dark matter and hot dark matter galaxy formation scenarios. We have derived peculiar velocities and associated errors for the scenarios using four values of the bias parameter ranging from b = 1 to b = 2.5. The growth of the mean peculiar velocity with scale factor has been determined and compared to that predicted by linear theory. In addition, we have compared the orientation of force and velocity in these simulations to see if a program such as that proposed by Bertschinger and Dekel (1989) for elliptical galaxy peculiar motions can be applied to clusters. The method they describe enables one to recover the density field from large scale redshift distance samples. The method makes it possible to do this when only radial velocities are known by assuming that the velocity field is curl free. Our analysis suggests that this program if applied to clusters is only realizable for models with a low value of the bias parameter, i.e., models in which the peculiar velocities of clusters are large enough that the errors do not render the analysis impracticable.

Rhee, George F.↗
NASA NTRS1993

An improved reconstruction method for cosmological density fields

This paper proposes some improvements to existing reconstruction methods for recovering the initial linear density and velocity fields of the universe from the present large-scale density distribution. We derive the Eulerian continuity equation in the Zel'dovich approximation and show that, by applying this equation, we can trace the evolution of the gravitational potential of the universe more exactly than is possible with previous approaches based on the Zel'dovich-Bernoulli equation. The improved reconstruction method is tested using N-body simulations. When the Zel'dovich-Bernoulli equation describes the formation of filaments, then the Zel'dovich continuity equation also follows the clustering of clumps inside the filaments. Our reconstruction method recovers the true initial gravitational potential with an rms error about 3 times smaller than previous methods. We examine the recovery of the initial distribution of Fourier components and find the scale at which the recovered phases are scrambled with respect their true initial values. Integrating the Zel'dovich continuity equation back in time, we can improve the spatial resolution of the reconstruction by a factor of about 2.

Gramann, Mirt↗
NASA NTRS1989

Cosmological velocity correlations - Observations and model predictions

By applying the present simple statistics for two-point cosmological peculiar velocity-correlation measurements to the actual data sets of the Local Supercluster spiral galaxy of Aaronson et al. (1982) and the elliptical galaxy sample of Burstein et al. (1987), as well as to the velocity field predicted by the distribution of IRAS galaxies, a coherence length of 1100-1600 km/sec is obtained. Coherence length is defined as that separation at which the correlations drop to half their zero-lag value. These results are compared with predictions from two models of large-scale structure formation: that of cold dark matter and that of baryon isocurvature proposed by Peebles (1980). N-body simulations of these models are performed to check the linear theory predictions and measure sampling fluctuations.

Gorski, Krzysztof M.↗