Engineering PapersSearch

Engineering topics

Babul, Arif

Publications and source records attributed to Babul, Arif.

Testing the gravitational instability hypothesis?

We challenge a widely accepted assumption of observational cosmology: that successful reconstruction of observed galaxy density fields from measured galaxy velocity fields (or vice versa), using the methods of gravitational instability theory, implies that the observed large-scale structures and large-scale flows were produced by the action of gravity. This assumption is false, in that there exist nongravitational theories that pass the reconstruction tests and gravitational theories with certain forms of biased galaxy formation that fail them. Gravitational instability theory predicts specific correlations between large-scale velocity and mass density fields, but the same correlations arise in any model where (a) structures in the galaxy distribution grow from homogeneous initial conditions in a way that satisfies the continuity equation, and (b) the present-day velocity field is irrotational and proportional to the time-averaged velocity field. We demonstrate these assertions using analytical arguments and N-body simulations. If large-scale structure is formed by gravitational instability, then the ratio of the galaxy density contrast to the divergence of the velocity field yields an estimate of the density parameter Omega (or, more generally, an estimate of beta identically equal to Omega(exp 0.6)/b, where b is an assumed constant of proportionality between galaxy and mass density fluctuations. In nongravitational scenarios, the values of Omega or beta estimated in this way may fail to represent the true cosmological values. However, even if nongravitational forces initiate and shape the growth of structure, gravitationally induced accelerations can dominate the velocity field at late times, long after the action of any nongravitational impulses. The estimated beta approaches the true value in such cases, and in our numerical simulations the estimated beta values are reasonably accurate for both gravitational and nongravitational models. Reconstruction tests that show correlations between galaxy density and velocity fields can rule out some physically interesting models of large-scale structure. In particular, successful reconstructions constrain the nature of any bias between the galaxy and mass distributions, since processes that modulate the efficiency of galaxy formation on large scales in a way that violates the continuity equation also produce a mismatch between the observed galaxy density and the density inferred from the peculiar velocity field. We obtain successful reconstructions for a gravitational model with peaks biasing, but we also show examples of gravitational and nongravitational models that fail reconstruction tests because of more complicated modulations of galaxy formation.

Babul, Arif

Reionization in a cold dark matter universe: The feedback of galaxy formation on the intergalactic medium

We study the coupled evolution of the intergalactic medium (IGM) and the emerging structure in the universe in the context of the cold dark matter (CDM) model, with a special focus on the consequences of imposing reionization and the Gunn-Peterson constraint as a boundary condition on the model. We have calculated the time-varying density of the IGM by coupling our detailed, numerical calculations of the thermal and ionization balance and radiative transfer in a uniform, spatially averaged IGM of H and He, including the mean opacity of an evolving distribution of gas clumps which correspond to quasar absorption line clouds, to the linearized equations for the growth of density fluctuations in both the gaseous and dark matter components in a CDM universe. We use the linear growth equations to identify the fraction of the gas which must have collapsed out at each epoch, an approach similar in spirit to the so-called Press-Schechter formalism. We identify the IGM density with the uncollapsed baryon fraction. The collapsed fraction is postulated to be a source of energy injection into the IGM, by radiation or bulk hydrodynamical heating (e.g., via shocks) or both, at a rate which is marginally enough to satisfy the Gunn-Peterson constraint at z less than 5. Our results include the following: (1) We find that the IGM in a CDM model must have contained a substantial fraction of the total baryon density of the universe both during and after its reionization epoch. (2) As a result, our previous conclusion that the observed Quasi-Stellar Objects (QSOs) at high redshift are not sufficient to ionize the IGM enough to satisfy the Gunn-Peterson constraint is confirmed. (3) We predict a detectable He II Gunn-Peterson effect at 304(1 + z) A in the spectra of quasars at a range of redshift z greater than or approx. 3, depending on the nature of the sources of IGM reionization. (4) We find, moreover, that a CDM model with high bias parameter b (i.e., b greater than or approx. 2) cannot account for the baryon content of the universe at z approximately 3 observed in quasar absorption line gas unless Omega (sub B) significantly exceeds the maximum value allowed by big bang nucleocynthesis. (5) For a CDM model with bias parameter within the allowed range of (lower) values, the lower limit to Omega(sub B) imposed by big bang nucleosynthesis (Omega(sub B) h(sup 2) greater than or equal to 0.01) combines with our results to yield the minimum IGM density for the CDM fodel. For CDM with b = 1 (Cosmic Background Explorer (COBE) normalization), we find Omega(sub IGM)(sup min) (z approximately 4) approx. equal 0.02-0.03, and Omega(sub IGM)(sup min)(z approximately 0) approx. equal 0.005-0.03, depending upon the nature of the sources of IGM reionization. (6) In general, we find that self-consistent reionization of the IGM by the collapsed baryon fraction has a strong effect on the rate of collapse. (7) As a further example, we show that the feedback effect on the IGM of energy release by the collapsed baryon fraction may explain the slow evolution of the observed comoving QSO number density between z = 5 and z = 2, followed by the sharp decline after z = 2.

Shapiro, Paul R.

The large-scale morphology of IRAS galaxies

At present, visual inspection is the only method for comparing the large-scale morphologies in the distribution of galaxies to those in model universes generated by N-body simulations. To remedy the situation, we have developed a set of three structure functions (S1, S2, S3) that quantify the degree of large-scale prolateness, oblateness, and sphericity/uniformity of a 3-D particle distribution and have applied them to a volume-limited (less than = 4000 km/s) sample of 699 IRAS galaxies with f sub 60 greater than 1.2 Jy. To determine the structure functions, we randomly select 500 galaxies as origins of spherical windows of radius R sub w, locate the centroid of the galaxies in the window (assuming all galaxies have equal mass) and then, compute the principal moments of inertia (I sub 1, I sub 2, I sub 3) about the centroid. Each S sub i is a function of (I sub 2)/(I sub 1) and (I sub 3)/I sub 1). S1, S2, and S3 tend to unity for highly prolate, oblate, and uniform distributions, respectively and tend to zero otherwise. The resulting 500 values of S sub i at each scale R sub w are used to construct a histogram.

Babul, Arif

Faint blue counts from formation of dwarf galaxies at z approximately equals 1

The nature of faint blue objects (FBO's) has been a source of much speculation since their detection in deep CCD images of the sky. Their high surface density argues against them being progenitors of present-day bright galaxies and since they are only weakly clustered on small scales, they cannot be entities that merged together to form present-day galaxies. Babul & Rees (1992) have suggested that the observed faint blue counts may be due to dwarf elliptical galaxies undergoing their initial starburst at z is approximately equal to 1. In generic hierarchical clustering scenarios, however, dwarf galaxy halos (M is approximately 10(exp 9) solar mass) are expected to form at an earlier epoch; for example, typical 10(exp 9) solar mass halos will virialize at z is approximately equal to 2.3 if the power-spectrum for the density fluctuations is that of the standard b = 2 cold dark matter (CDM) model. Under 'ordinary conditions' the gas would rapidly cool, collect in the cores, and undergo star-formation. Conditions at high redshifts are far from 'ordinary'. The intense UV background will prevent the gas in the dwarf halos from cooling, the halos being released from their suspended state only when the UV flux has diminished sufficiently.

Babul, Arif

Testing the pressure-confined Ly alpha cloud model

The Ly-alpha absorption line forest, seen in quasar spectra, is generally interpreted as being due to cosmologically distributed 'clouds' of primordial gas. Analyses of the observations reveal that the number distribution can be described by power laws: dN/dz is proportional to (1 + z)(sup gamma) and dN/dN(sub HI) is proportional to N(sub HI)(sup -beyda), where N(sub HI) is the HI column density. The typical values for power law indices range between 2 is approximately less than gamma is approximately less than 2.6 and 1.7 is approximately less than gamma is approximately less than 1.9. One model postulates that the Ly-alpha clouds are optically thin entities, photoionized by the background UV flux, J(sub nu) is proportional to (1 + z)(sup j), and confined by an adiabatically evolving intercloud medium (ICM): P(z) is proportional to (1 + z)(sup 5). Analytic studies of this model suggest that the ensuing Ly-alpha line statistics can account for the observations (in particular, the dN/dz and the dN/dN(sub HI) distributions) if the cloud mass spectrum is a power law dN/dN is proportional to M(sup -delta), delta is approximately 1.9, and j is approximately 4. One of the simplifying assumptions incorporated into these studies is the existence of a large mass range for the clouds at all epochs, the validity of which is questionable. The pressure-confined model is investigated using a 1-D spherically symmetric hydrodynamical code to simulate cloud evolution over the epoch 1.8 less than z less than 6. This enables us to relax many of the assumptions incorporated in the analytic studies.

Babul, Arif

Modelling the spatial distribution of QSO absorption lines

A formalism has been developed for testing specific models for the three-dimensional spatial distributions of Ly-alpha clouds against the line-of-sight observations. The formalism can be used to compute the redshift density of lines and the distribution of line-interval sizes along a line-of-sight for comparison with the observations. A few simple models are considered: a Poisson distribution of Ly-alpha clouds, clouds arising from density fluctuations with a white-noise power spectrum, clouds distributed only on shells of constant comoving radii, clouds restricted to shells that are expanding in a self-similar fashion, and a Poisson distribution of clouds in the 'walls' between small-scale voids.

Babul, Arif

The evolving intergalactic medium - The uncollapsed baryon fraction in a cold dark matter universe

The time-varying density of the intergalactic medium (IGM) is calculated by coupling detailed numerical calculations of the thermal and ionization balance and radiative transfer in a uniform IGM of H and He to the linearized equations for the growth of density fluctuations in both gases and a dark component in a cold dark matter universe. The IGM density is identified with the collapsed baryon fraction. It is found that even if the IGM is never reheated, a significant fraction of the baryons remain uncollapsed at redshifts of four. If instead the collapsed fraction releases enough ionizing radiation or thermal energy to reionize the IGM by z greater than four as required by the Gunn-Peterson (GP) constraint, the uncollapsed fraction at z of four is even higher. The known quasar distribution is insufficient to supply the ionizing radiation necessary to satisfy the GP constraint in this case and, if stars are instead responsible, a substantial metallicity must have been produced by z of four.

Shapiro, Paul R.

IRAS galaxies and the large-scale structure in the CfA slice

The spatial distributions of the IRAS and the optical galaxies in the first CfA slice are compared. The IRAS galaxies are generally less clustered than optical ones, but their distribution is essentially identical to that of late-type optical galaxies. The discrepancy between the clustering properties of the IRAS and optical samples in the CfA slice region is found to be entirely due to the paucity of IRAS galaxies in the core of the Coma cluster. The spatial distributions of the IRAS and the optical galaxies, both late and early types, outside the dense core of the Coma cluster are entirely consistent with each other. This conflicts with the prediction of the linear biasing scenario.

Babul, Arif

Ly-alpha clouds - Isolated structures in cold dark matter cosmogony?

A filter designed to study isolated density fluctuations is introduced and used to identify Ly-alpha cloud candidates in a biased CDM model. It is found that the number density of cloud candidates at z = 2.5 is too small to account for the observed number of lines. If Ly-alpha clouds have their origins in primordial density fluctuations, the CDM power spectrum does not have sufficient power on subgalactic scales to account for the observed structures.

Babul, Arif

Gamma-ray bursts from superconducting cosmic strings at large redshifts

The relation between cusp events and gamma-rays bursts is investigated. The optical depth of the universe to X-rays and gamma-rays of various energies is calculated and discussed. The cosmological evolution of cosmic strings is examined, and the energetics and time-scales related to the cusp phenomena are estimated. It is noted that it is possible to have energy bursts with a duration of a few seconds or less from cusps at z = 1000; the maximum amount of energy associated with such an event is limited to 10 to the 7th ergs/sq cm.

Babul, Arif