Mathematical modelling by symbolic mathematical computation - A cosmological application
Cosmology application of mathematical modeling by symbolic computation
SEARCH · Search NASA
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.
Cosmology application of mathematical modeling by symbolic computation
Numerical relativity, the solution of the Einstein equations on a computer, is one of the most challenging and exciting areas of physics. Richard Matzner has played a key role in this subject from its birth, roughly 3 decades ago, to the present. This talk will present some of the highlights of Richard's work in numerical relativity.
One-dimensional model of stellar system evolution, using computer to calculate minimum energy configuration
We compute the evolution of spatially flat, mixed cold and hot dark matter models containing both baryonic matter and two kinds of dark matter. Hydrodynamics is treated with a highly developed Eulerian hydrodynamic code (see Cen 1992). A standard particle-mesh (PM) code is also used in parallel to calculate the motion of the dark matter components. We adopt the following parameters: h equivalent to (sub 0)/100 km/s Mpc(exp -1) = 0.5, OMEGA(sub C) = 0.3, and OMEGA(sub B) = 0.06, with amplitude of the perturbation spectrum fixed by the Cosmic Background Explorer Satellite (COBE) Dark Matter Radiation (DMR) measurements (Smoot et al. 1992) being sigma (sub 8) = 0.67. Four different boxes are simulated with box sizes of L = (64, 16, 4, 1) h(exp -1) Mpc, respectively, the two small boxes providing good resolution but little valid information due to absence of large-scale power. We use 128(exp 3) approximate 10(exp 6.3) baryonic cells, 128(exp .3) cold dark matter particles, and 2 x 128(exp 3) hot dark matter particles. In addition to the dark matter we follow separately six baryonic species (H, H(+), He, He(+), He(++), e(-)) with allowance for both (nonequilibrium) collisional and radiative ionization in every cell. The background radiation field is also followed in detail with allowance made for both continuum and line processes, to allow nonequilibrium heating and cooling processes to be followed in detail. The mean final Zeldovich-Sunyaev y parameter is estimated to be y Bar = (5.4 + or - 2.7) x 10(exp -7) below currently attainable observations, with a rms fluctuation of approximately delta bar y = (0.6 + or - 3.0) x 10(exp -7) on arcminute scales. The rate of galaxy formation peaks at an even later epoch (z approximate 0.3) than in the standard (OMEGA = 1, sigma sub 8 = 0.67) cold dark matter (CDM) model (z approximate 0.5) and, at a redshift of z = 4, is nearly a factor of 100 lower than for the CDM model with the same value of sigma sub 8. With regard to mass function, the smallest objects are stabilized against collapse by thermal energy: the mass-weighted mass spectrum has a broad peak in the vicinity of M(sub B) = 10(exp 9.5) solar mass with a reasonable fit to the Schechter luminosity function if the ratio of baryon mass to blue light is approximately 4. In addition, one very large PM simulation was made in a box with size (320 h(exp - 1) Mpc) containing 3 x 200(exp 3) = 10(exp 7.4) particles. Utilizing this simulation we find that the model yields a cluster mass function which is about a factor of 4 higher than observed, but a cluster-cluster correlation length marginally lower than observed, but that both are closer to observations than in the (COBE) normalized CDM model. The one-dimensional pairwise velocity dispersion is 605 + or - 8 km/s at 1/h separation, lower than that of the DCM model normalized to COBE, but still significant higher than observations (Davis & Peebles 1983). A plausible velocity bias b(sub v) = 0.8 + or - 0.1 on this scale will reduce but not remove the discrepancy. The velocity auto-correlat ion function has a coherence length of 40/h Mpc, which is somewhat lower than the observed counterpart. In all these respects the model would be improved by decreasing the cold fraction of the dark OMEGA(sub CDM)/ (OMEGA(sub CDM) + OMEGA(sub HDB). But formation of galaxies and clusters of galaxies is much later in this model than in COBE-normalized CDM, perhaps too late. To improve on these constraints a larger ratio of OMEGA(sub CDM)/ (OMEGA(sub CDM) + OMEGA(sub HDM)) is required than the value of 0.67 adopted here. It does not seem possible to find a value for this ratio which would satisfy all tests. Overall, the model is similar both on large and intermediate scales to the standard CDM model normalized to the same value of sigma(sub B), but the problem with regard to late formation of galaxies is more severe in this model than in that CDM model. Adding hot dark matter, significantly improves the ability of the COBE-normalized CDM scenario to fit existing observations, but the model is in fact not as good as the CDM model with the same sigma(sub 8) and is still probably unsatisfactory with regard to several critical tests.
The study computes the evolution of the hot dark matter (HDM) model containing both baryonic matter and dark matter for a post recombination Friedmann-Robertson-Walker universe. A locally valid Newtonian approximation is used to model a representative piece of the universe with size much less than the horizon. For the HDM model with the present chosen normalization, the hard X-ray (1-10 keV) radiation intensity is less than that in the observations (Wu et al., 1991) by a factor of 30. In agreement with other work, it is found that baryonic matter is slightly antibiased over dark matter on the cell scale, 0.5/h Mpc = 667 kpc. The HDM model with the present chosen parameters does not overproduce X-ray-luminous clusters, and there is a negative evolution in the late epochs in the sense that the number density of X-ray clusters was higher at 0.5 redshift than at 0 redshift at the brightest end.
Explore the source record for details and available documents.
We study the topology of large-scale structure in the Center for Astrophysics Redshift Survey, which now includes approximately 12,000 galaxies with limiting magnitude m(sub B) is less than or equal to 15.5. The dense sampling and large volume of this survey allow us to compute the topology on smoothing scales from 6 to 20/h Mpc; we thus examine the topology of structure in both 'nonlinear' and 'linear' regimes. On smoothing scales less than or equal to 10/h Mpc this sample has 3 times the number of resolution elements of samples examined in previous studies. Isodensity surface of the smoothed galaxy density field demonstrate that coherent high-density structures and large voids dominate the galaxy distribution. We compute the genus-threshold density relation for isodensity surfaces of the CfA survey. To quantify phase correlation in these data, we compare the CfA genus with the genus of realizations of Gaussian random fields with the power spectrum measured for the CfA survey. On scales less than or equal to 10/h Mpc the observed genus amplitude is smaller than random phase (96% confidence level). This decrement reflects the degree of phase coherence in the observed galaxy distribution. In other words the genus amplitude on these scales is not good measure of the power spectrum slope. On scales greater than 10/h Mpc, where the galaxy distribution is rougly in the 'linear' regime, the genus ampitude is consistent with the random phase amplitude. The shape of the genus curve reflects the strong coherence in the observed structure; the observed genus curve appears broader than random phase (94% confidence level for smoothing scales less than or equal to 10/h Mpc) because the topolgoy is spongelike over a very large range of density threshold. This departre from random phase consistent with a distribution like a filamentary net of 'walls with holes.' On smoothing scales approaching approximately 20/h Mpc the shape of the CfA genus curve is consistent with random phase. There is very weak evidence for a shift of the genus toward a 'bubble-like' topology. To test cosmological models, we compute the genus for mock CfA surveys drawn from large (L greater than or approximately 400/h Mpc) N-body simulations of three variants of the cold dark matter (CDM) cosmogony. The genus amplitude of the 'standard' CDM model (omega h = 0.5, b = 1.5) differs from the observations (96% confidence level) on smoothing scales is less than or approximately 10/h Mpc. An open CDM model (omega h = 0.2) and a CDM model with nonzero cosmological constant (omega h = 0.24, lambda (sub 0) = 0.6) are consistent with the observed genus amplitude over the full range of smoothing scales. All of these models fail (97% confidence level) to match the broadness of the observed genus curve on smoothing scales is less than or equal to 10/h Mpc.
We assess the effects of the stochastic attenuation produced by intervening QSO absorption systems on the broadband colors of galaxies at cosmological distances. We compute the H I opacity of a clumpy universe as a function of redshift, including scattering in resonant lines, such as Lyman alpha, Lyman beta, Lyman gamma, and higher order members, and Lyman-continuum absorption. Both the numerous, optically thin Lyman-alpha forest clouds and the rarer, optically thick Lyman limit systems are found to contribute to the obscuration of background sources. We study the mean properties of primeval galaxies at high redshift in four broad optical passbands, U(sub n), B, G, and R. Even if young galaxies radiated a significant amount of ionizing photons, the attenuation due to the accumulated photoelectric opacity along the path is so severe that sources beyond z approximately 3 will drop out of the U(sub n) image together. We also show that the observed B-R color of distant galaxies can be much redder than expected from a stellar population. At z approximately 3.5, the blanketing by discrete absorption lines in the Lyman series is so effective that background galaxies appear, on average, 1 mag fainter in B. By z approximately 4, the observed B magnitude increment due to intergalactic absorption exceeds 2 mag. By modeling the intrinsic UV spectral energy distribution of star-forming galaxies with a stellar population synthesis code, we show that the (B-R)(sub AB) approximately 0 criterion for identifying 'flat-spectrum,' metal-producing galaxies is biased against objects at z greater than 3. The continuum blanketing from the Lyman series produces a characteristic staircase profile in the transmitted power. We suggest that this cosmic Lyman decrement might be used as a tool to identify high-z galaxies.
We present the analytical formulae for computing the magnification probability caused by cosmologically distributed galaxies. The galaxies are assumed to be singular, truncated-isothermal spheres without both evolution and clustering in redshift. We find that, for a fixed total mass, extended galaxies produce a broader shape in the magnification probability distribution and hence are less efficient as gravitational lenses than compact galaxies. The high-magnification tail caused by large galaxies is well approximated by an A exp -3 form, while the tail by small galaxies is slightly shallower. The mean magnification as a function of redshift is, however, found to be independent of the size of the lensing galaxies. In terms of the flux conservation, our formulae for the isothermal galaxy model predict a mean magnification to within a few percent with the Dyer-Roeder model of a clumpy universe.
Inhomogeneities in the large-scale distribution of matter inevitably lead to the generation of large-scale anisotropy in the cosmic background radiation. The dipole, quadrupole, and higher order fluctuations expected in an Einstein-de Sitter cosmological model have been computed. The dipole and quadrupole anisotropies are comparable to the measured values, and impose important constraints on the allowable spectrum of large-scale matter density fluctuations. A significant dipole anisotropy is generated by the matter distribution on scales greater than approximately 100 Mpc. The large-scale anisotropy is insensitive to the ionization history of the universe since decoupling, and cannot easily be reconciled with a galaxy formation theory that is based on primordial adiabatic density fluctuations.
A scalar field phi with a potential V(phi) varies as phi exp -alpha(alpha is greater than 0) has an energy density, behaving like that of a time-variable cosmological 'constant', that redshifts less rapidly than the energy densities of radiation and matter, and so might contribute significantly to the present energy density. We compute, in this spatially flat cosmology, the gravitational lensing optical depth, and the expected lens redshift distribution for fixed source redshift. We find, for the values of alpha of about 4 and baryonic density parameter Omega of about 0.2 consistent with the classical cosmological tests, that the optical depth is significantly smaller than that in a constant-Lambda model with the same Omega. We also find that the redshift of the maximum of the lens distribution falls between that in the constant-Lambda model and that in the Einstein-de Sitter model.
The volume contains recommendations for space-and ground-based programs in gravitational physics, cosmology, and cosmic-ray physics. The section on gravitation examines current and planned experimental tests of general relativity; the theory behind, and search for, gravitational waves, including sensitive laser-interferometric tests and other observations; and advances in gravitation theory (for example, incorporating quantum effects). The section on cosmology deals with the big-bang model, the standard model from elementary-particle theory, the inflationary model of the Universe. Computational needs are presented for both gravitation and cosmology. Finally, cosmic-ray physics theory (nucleosynthesis, acceleration models, high-energy physics) and experiment (ground and spaceborne detectors) are discussed.
Variations in the electron pressure can produce significant distortions of the cosmic microwave background radiation (CMBR) spectrum through the Sunyaev-Zel'dovich effect. In this paper, we compute the amplitude of these distortions in various cosmological scenarios using the results of hydrodynamical simulations combined with analytic theory. For reionized cold dark matter (CDM) with h = 0.5 probed with the OVRO experiment (Theta approximately = 0.1 deg), we find that the amplitude of these fluctuations is (Delta T/T)(sub rms) = 1 x 10(exp -5) (Omega(sub i)/0.06) if the experiment looks at a random part of the sky and 6 x 10(exp -6)(Omega(sub i)/(0.06) if bright areas are avoided. Here Omega is the density of ionized baryons in units of the critical density, rho(sub c) = 2.78 x 10(exp 11)h(exp 2)/cu Mpc. Primary fluctuations would generate (Delta T/T)(sub rms) approximately = 2 x 10(exp -6) for this model. These distortions could be a dominant source of CMBR anisotropics on small and intermediate angular scales. It is important to note that these signals are generated relatively recently, when nonlinear effects are improtant, and do not require early reionization. We also examine the contribution to anisotropies from the Doppler effect, in both the linear and nonlinear regimes. In this case, the bulk of the signal is generated before a redshift of 10, provided that the intercluster medium is ionized back to the surface of last scattering. On the smallest scales thermal effects are bigger, while, on larger scales the Doppler-induced fluctuations are greater with the exact boundary determined by the scenario. However, over a wide range of scales, the two signals are equal to within a factor of a few.
A model for the photometric evolution of galaxies has been developed and has been applied to the problem of galaxy counts. The integrated colors of galaxies are calculated using the most recently computed evolutionary tracks from Maeder and collaborators complemented with evolutionary tracks derived by other authors. The asymptotic giant branch lifetime is left as a free parameter. A series of cosmological models using different values of the cosmological constant, lambda(sub 0), and the density parameter, omega(sub 0), have been computed. The universality hypothesis of the luminosity function of galaxies has been abandoned. The influence of galaxy merging on the counts has been considered in a simple manner by assuming that the number of strongly interacting galaxies in a comoving volume increases with redshift as a power law given by (1 + z)(exp 3.8). Taking a Schechter parametrization for the luminosity function of the different types of galaxies, we are able to reproduce the observations reasonably well. We have also considered models with a Gaussian distribution for the luminosity function of the brighter galaxies that provide a poorer fit to the observations. It is shown that galaxy count data are not yet able to make unambiguous cosmological statements since evolutionary assumptions are critical. In particular, an omega(sub 0) = 1, lambda(sub 0) = 0 cosmology is shown to be consistent with the data.
We present the final nine-year maps and basic results from the Wilkinson Microwave Anisotropy Probe (WMAP) mission. The full nine-year analysis of the time-ordered data provides updated characterizations and calibrations of the experiment. We also provide new nine-year full sky temperature maps that were processed to reduce the asymmetry of the effective beams. Temperature and polarization sky maps are examined to separate cosmic microwave background (CMB) anisotropy from foreground emission, and both types of signals are analyzed in detail.We provide new point source catalogs as well as new diffuse and point source foreground masks. An updated template-removal process is used for cosmological analysis; new foreground fits are performed, and new foreground reduced are presented.We nowimplement an optimal C(exp -1)1 weighting to compute the temperature angular power spectrum. The WMAP mission has resulted in a highly constrained Lambda-CDM cosmological model with precise and accurate parameters in agreement with a host of other cosmological measurements. When WMAP data are combined with finer scale CMB, baryon acoustic oscillation, and Hubble constant measurements, we find that big bang nucleosynthesis is well supported and there is no compelling evidence for a non-standard number of neutrino species (N(sub eff) = 3.84 +/- 0.40). The model fit also implies that the age of the universe is (sub 0) = 13.772 +/- 0.059 Gyr, and the fit Hubble constant is H(sub 0) = 69.32 +/- 0.80 km/s/ Mpc. Inflation is also supported: the fluctuations are adiabatic, with Gaussian random phases; the detection of a deviation of the scalar spectral index from unity, reported earlier by the WMAP team, now has high statistical significance (n(sub s) = 0.9608+/-0.0080); and the universe is close to flat/Euclidean (Omega = −0.0027+0.0039/−0.0038). Overall, the WMAP mission has resulted in a reduction of the cosmological parameter volume by a factor of 68,000 for the standard six-parameter Lambda-CDM model, based on CMB data alone. For a model including tensors, the allowed seven-parameter volume has been reduced by a factor 117,000. Other cosmological observations are in accord with the CMB predictions, and the combined data reduces the cosmological parameter volume even further.With no significant anomalies and an adequate goodness of fit, the inflationary flat Lambda-CDM model and its precise and accurate parameters rooted in WMAP data stands as the standard model of cosmology.
A cosmological hydrodynamic code is described, which includes a routine to compute cosmic ray acceleration and transport in a simplified way. The routine was designed to follow explicitly diffusive, acceleration at shocks, and second-order Fermi acceleration and adiabatic loss in smooth flows. Synchrotron cooling of the electron population can also be followed. The updated code is intended to be used to study the properties of nonthermal synchrotron emission and inverse Compton scattering from electron cosmic rays in clusters of galaxies, in addition to the properties of thermal bremsstrahlung emission from hot gas. The results of a test simulation using a grid of 128 (exp 3) cells are presented, where cosmic rays and magnetic field have been treated passively and synchrotron cooling of cosmic ray electrons has not been included.
We obtain predictions for cosmic microwave background anisotropies at angular scales near 1 deg in the context of cold dark matter models with a nonzero cosmological constant, normalized to the Cosmic Background Explorer (COBE) Differential Microwave Radiometer (DMR) detection. The results are compared to those computed in the matter-dominated models. We show that the coherence length of the Cosmic Microwave Background (CMB) anisotropy is almost insensitive to cosmological parameters, and the rms amplitude of the anisotropy increases moderately with decreasing total matter density, while being most sensitive to the baryon abundance. We apply these results in the statistical analysis of the published data from the UCSB South Pole (SP) experiment (Gaier et al. 1992; Schuster et al. 1993). We reject most of the Cold Dark Matter (CDM)-Lambda models at the 95% confidence level when both SP scans are simulated together (although the combined data set renders less stringent limits than the Gaier et al. data alone). However, the Schuster et al. data considered alone as well as the results of some other recent experiments (MAX, MSAM, Saskatoon), suggest that typical temperature fluctuations on degree scales may be larger than is indicated by the Gaier et al. scan. If so, CDM-Lambda models may indeed provide, from a point of view of CMB anisotropies, an acceptable alternative to flat CDM models.
A set of computer models for the evolution of intercluster matter (ICM) in expanding cosmologies which includes processes of electron impact ionization/excitation, radiative and dielectronic recombination, thermal bremsstrahlung, Compton scattering, and photoionization is presented. The thermal, ionization, and excitation histories of major constituents were found as functions of time, and predictions of emission spectra, the possible distortions of the cosmic 3 K microwave background, and optical depths of absorption troughs in QSO continua at various redshifts were obtained. It is concluded that the present Friedmann cosmology can be dominated by hot intercluster matter which would contain H and He in their cosmic ratio and traces of heavier elements. For a narrow range of parameters it is possible to generate cosmologies with a closure density of ICM which do not exceed the diffuse X-ray measurements. In open cosmologies distortion of the 3 K microwave background occurs, while closed cases show negligible change from the unperturbed spectrum.