-
The Clustering of High Redshift Galaxies in the Cold Dark Matter Scenario
Authors:
Neal Katz,
Lars Hernquist,
David H. Weinberg
Abstract:
We investigate the clustering of high redshift galaxies in five variants of the cold dark matter (CDM) scenario, using hydrodynamic cosmological simulations that resolve the formation of systems with circular velocities $v_c \geq 100 \kms$ ($Ω=1$) or $v_c \geq 70 \kms$ ($Ω=0.4$). Although the five models differ in their cosmological parameters and in the shapes and amplitudes of their mass power…
▽ More
We investigate the clustering of high redshift galaxies in five variants of the cold dark matter (CDM) scenario, using hydrodynamic cosmological simulations that resolve the formation of systems with circular velocities $v_c \geq 100 \kms$ ($Ω=1$) or $v_c \geq 70 \kms$ ($Ω=0.4$). Although the five models differ in their cosmological parameters and in the shapes and amplitudes of their mass power spectra, they predict remarkably similar galaxy clustering at $z=2$, 3, and 4. The galaxy correlation functions show almost no evolution over this redshift range, even though the mass correlation functions grow steadily in time. Despite the fairly low circular velocity threshold of the simulations, the high redshift galaxies are usually highly biased tracers of the underlying mass distribution; the bias factor evolves with redshift and varies from model to model. Predicted correlation lengths for the resolved galaxy population are $2-3\hmpc$ (comoving) at $z=3$. More massive galaxies tend to be more strongly clustered. These CDM models have no difficulty in explaining the strong observed clustering of Lyman-break galaxies, and some may even predict excessive clustering. Because the effects of bias obscure differences in mass clustering, it appears that Lyman-break galaxy clustering will not be a good test of cosmological models but will instead provide a tool for constraining the physics of galaxy formation.
△ Less
Submitted 18 June, 1998;
originally announced June 1998.
-
Reconstruction Analysis of Galaxy Redshift Surveys: A Hybrid Reconstruction Method
Authors:
Vijay K. Narayanan,
David H. Weinberg
Abstract:
In reconstruction analysis of galaxy redshift surveys, one works backwards from the observed galaxy distribution to the primordial density field in the same region, then evolves the primordial fluctuations forward in time with an N-body code. This incorporates assumptions about the cosmological parameters, the properties of primordial fluctuations, and the biasing relation between galaxies and m…
▽ More
In reconstruction analysis of galaxy redshift surveys, one works backwards from the observed galaxy distribution to the primordial density field in the same region, then evolves the primordial fluctuations forward in time with an N-body code. This incorporates assumptions about the cosmological parameters, the properties of primordial fluctuations, and the biasing relation between galaxies and mass. These can be tested by comparing the reconstruction to the observed galaxy distribution, and to peculiar velocity data. This paper presents a hybrid reconstruction method that combines the `Gaussianization'' technique of Weinberg(1992) with the dynamical schemes of Nusser & Dekel(1992) and Gramann(1993). We test the method on N-body simulations and on N-body mock catalogs that mimic the depth and geometry of the Point Source Catalog Redshift Survey and the Optical Redshift Survey. This method is more accurate than Gaussianization or dynamical reconstruction alone. Matching the observed morphology of clustering can limit the bias factor b, independent of Omega. Matching the cluster velocity dispersions and z-space distortions of the correlation function xi(s,mu) constrains the parameter beta=Omega^{0.6}/b. Relative to linear or quasi-linear approximations, a fully non-linear reconstruction makes more accurate predictions of xi(s,mu) for a given beta, thus reducing the systematic biases of beta measurements and offering further scope for breaking the degeneracy between Omega and b. It also circumvents the cosmic variance noise that limits conventional analyses of xi(s,mu). It can also improve the determination of Omega and b from joint analyses of redshift
& peculiar velocity surveys as it predicts the fully non-linear peculiar velocity distribution at each point in z-space.
△ Less
Submitted 17 June, 1998;
originally announced June 1998.
-
Cosmological Simulations with Scale-Free Initial Conditions II: Radiative Cooling
Authors:
J. Michael Owen,
David H. Weinberg,
Jens V. Villumsen
Abstract:
The growth of structure from scale-free initial conditions is one of the most important tests of cosmological simulation methods, providing a realistically complex problem in which numerical results can be compared to rigorous analytic scaling laws. Previous studies of this problem have incorporated gravitational dynamics and adiabatic gas dynamics, but radiative cooling, an essential element of…
▽ More
The growth of structure from scale-free initial conditions is one of the most important tests of cosmological simulation methods, providing a realistically complex problem in which numerical results can be compared to rigorous analytic scaling laws. Previous studies of this problem have incorporated gravitational dynamics and adiabatic gas dynamics, but radiative cooling, an essential element of the physics of galaxy formation, normally violates the conditions necessary for self-similar evolution. We show that for any specified value of the initial power spectrum index n [where P(k)~k^n], there is a family of power-law cooling functions that preserves self-similarity by ensuring that the cooling time in an object of the characteristic mass M_* is a fixed fraction t_C of the Hubble time. We perform hydrodynamic numerical simulations with an Einstein-de Sitter cosmology, a baryon fraction of 5%, Gaussian initial conditions, two different power spectrum indices, and four values of t_C for each index, ranging from no cooling to strong cooling. We restrict the numerical simulations to two dimensions in order to allow exploration of a wide parameter space with adequate dynamic range. In all cases, the simulations are remarkably successful at reproducing the analytically predicted scalings of the mass function of dissipated objects and the gas temperature distributions and cooled gas fractions in collapsed systems. While similar success with 3-D simulations must still be demonstrated, our results have encouraging implications for numerical studies of galaxy formation, indicating that simulations with resolution comparable to that in many current studies can accurately follow the collapse and dissipation of baryons into the dense, cold systems where star formation is likely to occur.
△ Less
Submitted 7 May, 1998;
originally announced May 1998.
-
Mock 2dF and SDSS galaxy redshift surveys
Authors:
Shaun Cole,
Steve Hatton,
David H. Weinberg,
Carlos S. Frenk
Abstract:
We present a comprehensive set of mock 2dF and SDSS galaxy redshift surveys constructed from a set of large, high-resolution cosmological N-body simulations. The radial selection functions and geometrical limits of the catalogues mimic those of the genuine surveys. The catalogues span a wide range of cosmologies, including both open and flat universes. In all the models the galaxy distributions…
▽ More
We present a comprehensive set of mock 2dF and SDSS galaxy redshift surveys constructed from a set of large, high-resolution cosmological N-body simulations. The radial selection functions and geometrical limits of the catalogues mimic those of the genuine surveys. The catalogues span a wide range of cosmologies, including both open and flat universes. In all the models the galaxy distributions are biased so as to approximately reproduce the observed galaxy correlation function on scales of 1--10 Mpc/h In some cases models with a variety of different biasing prescriptions are included. All the mock catalogues are publically available at http://star-www.dur.ac.uk/~cole/mocks/main.html . We expect these catalogues to be a valuable aid in the development of the new algorithms and statistics that will be used to analyse the 2dF and SDSS surveys when they are completed in the next few years. Mock catalogues of the PSCZ survey of IRAS galaxies are also available at the same WWW location.
△ Less
Submitted 26 June, 1998; v1 submitted 26 January, 1998;
originally announced January 1998.
-
Simulating Gas at High Redshift
Authors:
Neal Katz,
Lars Hernquist,
David H. Weinberg
Abstract:
We discuss simulations of gas at high redshift. We briefly review the methods used and the results for quasar absorption lines. We present gas mass functions and galaxy correlation functions for 5 different cosmological models. Galaxies should be detectable at redshifts greater than 2 by SKAI, and measurements of the gas mass functions and galaxy correlation functions could be used to discrimina…
▽ More
We discuss simulations of gas at high redshift. We briefly review the methods used and the results for quasar absorption lines. We present gas mass functions and galaxy correlation functions for 5 different cosmological models. Galaxies should be detectable at redshifts greater than 2 by SKAI, and measurements of the gas mass functions and galaxy correlation functions could be used to discriminate between different cosmological models.
△ Less
Submitted 9 January, 1998;
originally announced January 1998.
-
Constraints on the Effects of Locally-biased Galaxy Formation
Authors:
Robert J. Scherrer,
David H. Weinberg
Abstract:
While it is well-known that "biased galaxy formation" can increase the strength of galaxy clustering, it is less clear whether straightforward biasing schemes can change the shape of the galaxy correlation function on large scales. Here we consider "local" biasing models, in which the galaxy density field $δ_g$ at a point is a function of the matter density field $δ$ at that point: $δ_g = f(δ)$.…
▽ More
While it is well-known that "biased galaxy formation" can increase the strength of galaxy clustering, it is less clear whether straightforward biasing schemes can change the shape of the galaxy correlation function on large scales. Here we consider "local" biasing models, in which the galaxy density field $δ_g$ at a point is a function of the matter density field $δ$ at that point: $δ_g = f(δ)$. We consider both deterministic biasing, where $f$ is simply a function, and stochastic biasing, in which the galaxy density is a random variable whose distribution depends on the matter density: $δ_g=X(δ)$. We show that even when this mapping is performed on a highly nonlinear density field with a hierarchical correlation structure, the correlation function $ξ$ is simply scaled up by a constant, as long as $ξ<< 1$. In stochastic biasing models, the galaxy autocorrelation function behaves exactly as in deterministic models, with the mean value of $X$ for a given value of $δ$ taking the role of the deterministic bias function. We extend our results to the power spectrum P(k), showing that for sufficiently small k, the effect of local biasing is equivalent to the multiplication of P(k) by a constant, with the addition of a constant term. If a cosmological model predicts a large-scale mass correlation function in conflict with the shape of the observed galaxy correlation function, then the model cannot be rescued by appealing to a complicated but local relation between galaxies and mass.
△ Less
Submitted 15 December, 1997;
originally announced December 1997.
-
New Statistical Measures of the Lya Forest Spectra for Accurate Comparison to Theoretical Models
Authors:
J. Miralda-Escude,
M. Rauch,
W. L. W. Sargent,
T. A. Barlow,
D. H. Weinberg,
L. Hernquist,
N. Katz,
R. Cen,
J. P. Ostriker
Abstract:
We propose a new method of analysis for the \lya forest, namely to measure the 1-point and 2-point joint probability distribution of the transmitted flux. The results for a sample of seven observed quasars and from two simulations of structure formation are shown and compared. Statistically significant differences in the 2-point function between the results of the numerical simulations and the o…
▽ More
We propose a new method of analysis for the \lya forest, namely to measure the 1-point and 2-point joint probability distribution of the transmitted flux. The results for a sample of seven observed quasars and from two simulations of structure formation are shown and compared. Statistically significant differences in the 2-point function between the results of the numerical simulations and the observations are easily found. The analysis we suggest is very simple to apply to observed data sets, and we discuss its superiority over the traditional Voigt-profile fitting algorithms for accurate comparison to the predictions of theoretical models.
△ Less
Submitted 21 October, 1997;
originally announced October 1997.
-
Hubble Flow Broadening of the Lyman-alpha Forest and its Implications
Authors:
David H. Weinberg,
Lars Hernquist,
Neal Katz,
Rupert Croft,
Jordi Miralda-Escudé
Abstract:
Lyman-alpha forest lines in QSO spectra have typical widths of 20-50 km/s. Low column density absorbers in cosmological simulations are large, diffuse structures, and the Hubble flow across the spatially extended absorber is usually the dominant contribution to the width of its associated absorption line. Thermal broadening is unimportant over most of the spectrum, and peculiar velocities tend t…
▽ More
Lyman-alpha forest lines in QSO spectra have typical widths of 20-50 km/s. Low column density absorbers in cosmological simulations are large, diffuse structures, and the Hubble flow across the spatially extended absorber is usually the dominant contribution to the width of its associated absorption line. Thermal broadening is unimportant over most of the spectrum, and peculiar velocities tend to make absorption features narrower rather than broader. As a consequence of Hubble flow broadening, there is a close relation between local Lyman-alpha optical depth and local neutral hydrogen density, which is well approximated by the Gunn-Peterson formula. The physics that governs the unshocked intergalactic medium leads to a tight correlation between the neutral hydrogen density and the underlying gas and dark matter overdensity. For many purposes, it is simpler to regard a Lyman-alpha forest spectrum as a continuous, non-linear map of the density field rather than a collection of discrete lines. This continuous field view of the Lyman-alpha forest can be applied to measurement of the baryon density parameter, testing of cosmological models, and robust determination of the shape and amplitude of the primordial mass power spectrum.
△ Less
Submitted 30 September, 1997;
originally announced September 1997.
-
A Study of Quasar Clustering at z>2.7 from the Palomar Transit Grism Survey
Authors:
Andrew W. Stephens,
Donald P. Schneider,
Maarten Schmidt,
James E. Gunn,
David H. Weinberg
Abstract:
The quest for structure indicators at earlier and earlier times in the evolution of the universe has led to the search for objects with ever higher redshifts. The Palomar Transit Grism Survey has produced a large sample of high redshift quasars ($z>2.7$), allowing statistical analysis of correlation between quasar positions. In this study, clustering is identified through comparison with…
▽ More
The quest for structure indicators at earlier and earlier times in the evolution of the universe has led to the search for objects with ever higher redshifts. The Palomar Transit Grism Survey has produced a large sample of high redshift quasars ($z>2.7$), allowing statistical analysis of correlation between quasar positions. In this study, clustering is identified through comparison with $100 000$ Monte Carlo generated, randomly populated volumes, which are identical to the observed region in spatial coordinates, redshift distribution, and number of quasars. Three pairs have been observed with comoving separations of 11.34, 12.97, and 24.13 h_{50}^{-1}Mpc (assuming $q_0=0.5$), smaller separations than would be expected to arise by chance in an unclustered distribution. Selection effects are ruled out as a false source of clustering by scrambling the observed quasar coordinates and redshifts, which gives a pair separation distribution nearly identical to that of the Monte Carlo distribution. Tests using the distribution of pair separations and nearest neighbor distances show that the observed pairs have a probability less than 0.1% of arising in an unclustered distribution. Using a maximum likelihood technique to estimate the correlation length $r_0$, assuming $ξ(r) = (r/r_0)^{-1.8}$, we find $r_0 = 35\pm15 h_{50}^{-1}Mpc$ (comoving, $q_0=0.5$, 1$σ$ errors), a value much larger than the correlation length of present-day galaxies.
△ Less
Submitted 26 August, 1997;
originally announced August 1997.
-
Simulating Cosmic Structure Formation
Authors:
David H. Weinberg,
Neal Katz,
Lars Hernquist
Abstract:
We describe cosmological simulation techniques and their application to studies of cosmic structure formation, with particular attention to recent hydrodynamic simulations of structure in the high redshift universe. Collisionless N-body simulations with Gaussian initial conditions produce a pattern of sheets, filaments, tunnels, and voids that resembles the observed large scale galaxy distributi…
▽ More
We describe cosmological simulation techniques and their application to studies of cosmic structure formation, with particular attention to recent hydrodynamic simulations of structure in the high redshift universe. Collisionless N-body simulations with Gaussian initial conditions produce a pattern of sheets, filaments, tunnels, and voids that resembles the observed large scale galaxy distribution. Simulations that incorporate gas dynamics and dissipation form dense clumps of cold gas with sizes and masses similar to the luminous parts of galaxies. Models based on inflation and cold dark matter predict a healthy population of high redshift galaxies, including systems with star formation rates of 20 M_{\sun}/year at z=6. At z~3, most of the baryons in these models reside in the low density intergalactic medium, which produces fluctuating Lyman-alpha absorption in the spectra of background quasars. The physical description of this ``Lyman-alpha forest'' is particularly simple if the absorption spectrum is viewed as a 1-dimensional map of a continuous medium instead of a collection of lines. The combination of superb observational data and robust numerical predictions makes the Lyman-alpha forest a promising tool for testing cosmological models.
△ Less
Submitted 24 August, 1997; v1 submitted 22 August, 1997;
originally announced August 1997.
-
Recovery of the Power Spectrum of Mass Fluctuations from Observations of the Lyman-alpha Forest
Authors:
Rupert A. C. Croft,
David H. Weinberg,
Neal Katz,
Lars Hernquist
Abstract:
We present a method to recover the shape and amplitude of the power spectrum of mass fluctuations, P(k), from observations of the high redshift \lya forest. The method is motivated by the physical picture of the \lya forest that has emerged from hydrodynamic cosmological simulations and related semi-analytic models, which predicts a tight correlation between the \lya optical depth and the underl…
▽ More
We present a method to recover the shape and amplitude of the power spectrum of mass fluctuations, P(k), from observations of the high redshift \lya forest. The method is motivated by the physical picture of the \lya forest that has emerged from hydrodynamic cosmological simulations and related semi-analytic models, which predicts a tight correlation between the \lya optical depth and the underlying matter density. We monotonically map the QSO spectrum to a Gaussian density field, measure its 3-d P(k), and normalize by evolving cosmological simulations with this P(k) until they reproduce the observed power spectrum of the QSO flux. Imposing the observed mean \lya opacity as a constraint makes the derived P(k) normalization insensitive to the choice of cosmological parameters, ionizing background spectrum, or reionization history. Thus, in contrast to estimates of P(k) from galaxy clustering, there are no uncertain "bias parameters" in the recovery of the mass power spectrum. We test the full procedure on SPH simulations of 3 cosmological models and show that it recovers their true mass power spectra on comoving scales ~1-10/h Mpc, the upper scale being set by the size of the simulation boxes. The procedure works even for noisy (S/N ~ 10), moderate resolution (~40 km/s pixels) spectra. We present an illustrative application to Q1422+231; the recovered P(k) is consistent with an Ω=1, σ_8=0.5 CDM model. Application to large QSO samples should yield the power spectrum of mass fluctuations on small scales at z ~ 2-4. (Compressed)
△ Less
Submitted 3 August, 1997;
originally announced August 1997.
-
Quasar Clustering and Spacetime Geometry
Authors:
Piotr A. Popowski,
David H. Weinberg,
Barbara S. Ryden,
Patrick S. Osmer
Abstract:
The non-Euclidean geometry of spacetime induces an anisotropy in the apparent correlation function of high-redshift quasars. This effect can constrain the cosmological constant Λindependent of any assumptions about evolution of luminosities, sizes, or clustering. We examine the prospects for distinguishing between low-density (Ω_0=0.1-0.4) cosmological models with flat and open space geometry us…
▽ More
The non-Euclidean geometry of spacetime induces an anisotropy in the apparent correlation function of high-redshift quasars. This effect can constrain the cosmological constant Λindependent of any assumptions about evolution of luminosities, sizes, or clustering. We examine the prospects for distinguishing between low-density (Ω_0=0.1-0.4) cosmological models with flat and open space geometry using the quasar samples anticipated from the 2dF and Sloan redshift surveys. We show that even these large quasar surveys are likely to reside in the "sparse sampling" regime, so that measurement errors in the correlation function obey Poisson statistics. As a result: (a) one can devise a simple maximum-likelihood scheme for estimating clustering parameters, (b) one can generate Monte Carlo realizations of correlation function measurements without creating artificial quasar distributions, and (c) for fixed quasar number, a deeper survey over a smaller area has greater statistical power than a shallow, large-area survey. Adopting recent (quite uncertain) estimates of the quasar correlation length, we find that the 2dF and Sloan samples can provide clear discrimination between flat and open geometries for Ω<= 0.2 but only marginal discrimination for Ω= 0.4. Clear discrimination is possible for Ω= 0.4 if the true quasar correlation length is a factor of two larger, and a high-density survey of 30,000 quasars in 200 square degrees would provide clear discrimination even for the lower correlation length.
△ Less
Submitted 15 July, 1997;
originally announced July 1997.
-
Testing Cosmological Models Against the Abundance of Damped Lyman-Alpha Absorbers
Authors:
Jeffrey P. Gardner,
Neal Katz,
David H. Weinberg,
Lars Hernquist
Abstract:
We calculate the number of damped Lyman-alpha absorbers expected in various popular cosmological models as a function of redshift and compare our predictions with observed abundances. The Press-Schechter formalism is used to obtain the distribution of halos with circular velocity in different cosmologies, and we calibrate the relation between circular velocity and absorption cross-section using…
▽ More
We calculate the number of damped Lyman-alpha absorbers expected in various popular cosmological models as a function of redshift and compare our predictions with observed abundances. The Press-Schechter formalism is used to obtain the distribution of halos with circular velocity in different cosmologies, and we calibrate the relation between circular velocity and absorption cross-section using detailed gas dynamical simulations of a ``standard'' cold dark matter (CDM) model. Because of this calibration, our approach makes more realistic assumptions about the absorption properties of collapsed objects than previous, analytic calculations of the damped Lyman-alpha abundance. CDM models with Omega_0=1, H_0=50, baryon density Omega_b=0.05, and scale-invariant primeval fluctuations reproduce the observed incidence and redshift evolution of damped Lyman-alpha absorption to within observational uncertainty, for both COBE normalization (sigma_8=1.2) and a lower normalization (sigma_8=0.7) that better matches the observed cluster abundance at z=0. A tilted (n=0.8, sigma_8=0.7) CDM model tends to underproduce absorption, especially at z=4. With COBE normalization, a CDM model with Omega_0=0.4, Omega_{Lambda}=0.6 gives an acceptable fit to the observed absorption; an open CDM model is marginally acceptable if Omega_0 is at least 0.4 and strongly inconsistent with the z=4 data if Omega_0=0.3. Mixed dark matter models tend not to produce sufficient absorption, being roughly comparable to tilted CDM models if Omega_{nu} = 0.2 and failing drastically if Omega_{nu} = 0.3.
△ Less
Submitted 15 May, 1997;
originally announced May 1997.
-
Cosmological Constraints from High-Redshift Damped Lyman-Alpha Systems
Authors:
Chung-Pei Ma,
Edmund Bertschinger,
Lars Hernquist,
David H. Weinberg,
Neal Katz
Abstract:
Any viable cosmological model must produce enough structure at early epochs to explain the amount of gas associated with high-redshift damped Ly$α$ systems. We study the evolution of damped Ly$α$ systems at redshifts $z\ge 2$ in cold dark matter (CDM) and cold+hot dark matter (CDM+HDM) models using both N-body and hydrodynamic simulations. Our approach incorporates the effects of gas dynamics, a…
▽ More
Any viable cosmological model must produce enough structure at early epochs to explain the amount of gas associated with high-redshift damped Ly$α$ systems. We study the evolution of damped Ly$α$ systems at redshifts $z\ge 2$ in cold dark matter (CDM) and cold+hot dark matter (CDM+HDM) models using both N-body and hydrodynamic simulations. Our approach incorporates the effects of gas dynamics, and we find that all earlier estimates which assumed that all the baryons in dark matter halos would contribute to damped Ly$α$ absorption have overestimated the column density distribution $f(N)$ and the fraction of neutral dense gas $Ω_g$ in damped Ly$α$ systems. The differences are driven by ionization of hydrogen in the outskirts of galactic halos and by gaseous dissipation near the halo centers, and they tend to exacerbate the problem of late galaxy formation in CDM+HDM models. We only include systems up to the highest observed column density $N\sim 10^{21.8}$ cm$^{-2}$ in the estimation of $Ω_g$ for a fair comparison with data. If the observed $f(N)$ and $Ω_g$ inferred from a small number of confirmed and candidate absorbers are robust, the amount of gas in damped Ly$α$ systems at high redshifts in the $Ω_ν=0.2$ CDM+HDM model falls well below the observations.
△ Less
Submitted 15 May, 1997;
originally announced May 1997.
-
Cosmological Simulations with Scale-Free Initial Conditions I: Adiabatic Hydrodynamics
Authors:
J. Michael Owen,
David H. Weinberg,
August E. Evrard,
Lars Hernquist,
Neal Katz
Abstract:
We analyze hierarchical structure formation based on scale-free initial conditions in an Einstein-de Sitter universe, including a baryonic component. We present three independent, smoothed particle hydrodynamics (SPH) simulations, performed with two different SPH codes (TreeSPH and P3MSPH) at two resolutions. Each simulation is based upon identical initial conditions, which consist of Gaussian d…
▽ More
We analyze hierarchical structure formation based on scale-free initial conditions in an Einstein-de Sitter universe, including a baryonic component. We present three independent, smoothed particle hydrodynamics (SPH) simulations, performed with two different SPH codes (TreeSPH and P3MSPH) at two resolutions. Each simulation is based upon identical initial conditions, which consist of Gaussian distributed initial density fluctuations that have an n=-1 power spectrum. The baryonic material is modeled as an ideal gas subject only to shock heating and adiabatic heating and cooling. The evolution is expected to be self-similar in time, and under certain restrictions we identify the expected scalings for many properties of the distribution of collapsed objects in all three realizations. The distributions of dark matter masses, baryon masses, and mass and emission weighted temperatures scale quite reliably. However, the density estimates in the central regions of these structures are determined by the degree of numerical resolution. As a result, mean gas densities and luminosities obey the expected scalings only when calculated within a limited dynamic range in density contrast. The temperatures and luminosities of the groups show tight correlations with the baryon masses, which can be well-represented by power-laws. The Press-Schechter (PS) approximation predicts the distribution of group dark matter and baryon masses fairly well, though it tends to overestimate the baryon masses. Combining the PS mass distribution with the measured relations for T(M) and L(M) predicts the temperature and luminosity distributions reasonably, though there are some discrepancies at high temperatures/luminosities. The three simulations agree well for the properties of groups that are resolved by 32 or more particles.
△ Less
Submitted 24 April, 1998; v1 submitted 14 May, 1997;
originally announced May 1997.
-
Large-scale structure in COBE-normalized cold dark matter cosmogonies
Authors:
Shaun Cole,
David H. Weinberg,
Carlos S. Frenk,
Bharat Ratra
Abstract:
We study the clustering of the mass distribution in COBE-normalized open and flat CDM models using large N-body simulations. With an age of the universe of 14 Gyr (12 Gyr) for the flat (open) models and a baryon density fixed by nucleosynthesis constraints, the observed abundance of rich galaxy clusters leads to tight constraints on the density parameter; 0.25<Omega_0<0.4 (0.4<Omega_0<0.5) for f…
▽ More
We study the clustering of the mass distribution in COBE-normalized open and flat CDM models using large N-body simulations. With an age of the universe of 14 Gyr (12 Gyr) for the flat (open) models and a baryon density fixed by nucleosynthesis constraints, the observed abundance of rich galaxy clusters leads to tight constraints on the density parameter; 0.25<Omega_0<0.4 (0.4<Omega_0<0.5) for flat (open) models. The upper limits on Omega_0 can be relaxed if one lowers the H_0 and increases the age of the universe, but h<0.25 is required to allow Omega_0=1. The constraints also change if one allows tilted primordial power spectra. An acceptable Omega_0=1 model with h=0.5 can be constructed, with a tilt of n=0.8. We compare the evolved mass correlation functions of the most promising models with those of APM galaxies. The flat models have steep correlation functions at small scales and require galaxies to be antibiased on scales r<8 Mpc/h. The open models require little or no antibias on small scales and a positive bias on large scales. The lack of a positive bias on small scales is difficult to reconcile with the virial analysis of galaxy clusters which, if Omega_0>0.2, implies that galaxies are overabundant in clusters relative to the field. The tilted Omega_0=1 model, on the other hand, does require that galaxies be positively biased on all scales. We also compute the topology of isodensity contours in these models, obtaining theoretical predictions that are less sensitive to galaxy bias.
△ Less
Submitted 10 February, 1997;
originally announced February 1997.
-
Characterization of Lyman Alpha Spectra and Predictions of Structure Formation Models: A Flux Statistics Approach
Authors:
Rupert A. C. Croft,
David H. Weinberg,
Lars Hernquist,
Neal Katz
Abstract:
In gravitational instability models, \lya absorption arises from a continuous fluctuating medium, so that spectra provide a non-linear one-dimensional ``map'' of the underlying density field. We characterise this continuous absorption using statistical measures applied to the distribution of absorbed flux. We describe two simple members of a family of statistics which we apply to simulated spect…
▽ More
In gravitational instability models, \lya absorption arises from a continuous fluctuating medium, so that spectra provide a non-linear one-dimensional ``map'' of the underlying density field. We characterise this continuous absorption using statistical measures applied to the distribution of absorbed flux. We describe two simple members of a family of statistics which we apply to simulated spectra in order to show their sensitivity as probes of cosmological parameters (H$_{0}$, $Ω$, the initial power spectrum of matter fluctuations) and the physical state of the IGM. We make use of SPH simulation results to test the flux statistics, as well as presenting a preliminary application to Keck HIRES data.
△ Less
Submitted 22 January, 1997;
originally announced January 1997.
-
The Topology of Large Scale Structure in the 1.2 Jy IRAS Redshift Survey
Authors:
Zacharias A. M. Protogeros,
David H. Weinberg
Abstract:
We measure the topology (genus) of isodensity contour surfaces in volume limited subsets of the 1.2 Jy IRAS redshift survey, for smoothing scales $λ=4\hmpc$, $7\hmpc$, and $12\hmpc$. At $12\hmpc$, the observed genus curve has a symmetric form similar to that predicted for a Gaussian random field. At the shorter smoothing lengths, the observed genus curve shows a modest shift in the direction of…
▽ More
We measure the topology (genus) of isodensity contour surfaces in volume limited subsets of the 1.2 Jy IRAS redshift survey, for smoothing scales $λ=4\hmpc$, $7\hmpc$, and $12\hmpc$. At $12\hmpc$, the observed genus curve has a symmetric form similar to that predicted for a Gaussian random field. At the shorter smoothing lengths, the observed genus curve shows a modest shift in the direction of an isolated cluster or ``meatball'' topology. We use mock catalogs drawn from cosmological N-body simulations to investigate the systematic biases that affect topology measurements in samples of this size and to determine the full covariance matrix of the expected random errors. We incorporate the error correlations into our evaluations of theoretical models, obtaining both frequentist assessments of absolute goodness-of-fit and Bayesian assessments of models' relative likelihoods. We compare the observed topology of the 1.2 Jy survey to the predictions of dynamically evolved, unbiased, gravitational instability models that have Gaussian initial conditions. The model with an $n=-1$, power-law initial power spectrum achieves the best overall agreement with the data, though models with a low-density cold dark matter power spectrum and an $n=0$ power-law spectrum are also consistent. The observed topology is inconsistent with an initially Gaussian model that has $n=-2$, and it is strongly inconsistent with a Voronoi foam model, which has a non-Gaussian, bubble topology.
△ Less
Submitted 20 January, 1997;
originally announced January 1997.
-
Metal Lines Associated with Lyman Alpha Absorbers: A Comparison of Theory and Observations
Authors:
Uffe Hellsten,
Romeel Dave,
Lars Hernquist,
David H. Weinberg,
Neal Katz
Abstract:
We study metal line absorption of CIV, CII, SiIV, and NV at redshifts z=3.5 to z=2 within the framework of a cosmological model for the Lyman alpha forest, comparing the results of numerical simulations to recent observations by Songaila & Cowie (1996, SC). In agreement with Rauch, Haehnelt & Steinmetz (1996), we find that the observed mean value of the CIV/HI ratio around z=3 is reproduced if a…
▽ More
We study metal line absorption of CIV, CII, SiIV, and NV at redshifts z=3.5 to z=2 within the framework of a cosmological model for the Lyman alpha forest, comparing the results of numerical simulations to recent observations by Songaila & Cowie (1996, SC). In agreement with Rauch, Haehnelt & Steinmetz (1996), we find that the observed mean value of the CIV/HI ratio around z=3 is reproduced if a uniform metallicity of [C/H] = -2.5 is assumed in our model, but that the observed scatter in this ratio is larger than predicted, implying a scatter in the metallicity of the absorbing systems of roughly an order of magnitude. We further argue that absorbers with relatively low column densities (log(N_HI) < 15) likely have a mean metallicity [C/H] less than -2.5, a result which is basically independent of the model considered, and similar to what is predicted in the simulations by Gnedin & Ostriker (1997). Our model predicts no significant evolution in the mean values of metal line column densities between z=3.5 and z=2. Comparison of the predicted and observed numbers of SiIV and NV systems suggests that the photoionizing background radiation field around z=3 is somewhat softer than that proposed by Haardt & Madau (1996). Our model does not account for the increase in the SiIV/CIV ratio at z=3.2 found by SC. While SC suggested that the increase could be explained by a softening of the radiation field due to HeII absorption at z > 3, such a modification does not raise the mean value of SiIV/CIV in our simulation because it shifts numerous weak SiIV systems to just above the detection limit, thus keeping the mean column density of observable SiIV systems low.
△ Less
Submitted 9 January, 1997;
originally announced January 1997.
-
A Lower Bound on the Cosmic Baryon Density
Authors:
David H. Weinberg,
Jordi Miralda-Escude,
Lars Hernquist,
Neal Katz
Abstract:
We derive analytic lower bounds on the cosmic baryon density by requiring that the high-redshift IGM contain enough neutral hydrogen to produce the observed \lya absorption in quasar spectra. The key theoretical assumption is that absorbing structures are no more extended in redshift space than in real space. This assumption might not hold if \lya clouds are highly overdense and thermally broade…
▽ More
We derive analytic lower bounds on the cosmic baryon density by requiring that the high-redshift IGM contain enough neutral hydrogen to produce the observed \lya absorption in quasar spectra. The key theoretical assumption is that absorbing structures are no more extended in redshift space than in real space. This assumption might not hold if \lya clouds are highly overdense and thermally broadened, but it is likely to hold in the gravitational instability picture for the \lya forest suggested by cosmological simulations, independently of the details of the cosmological model. The other ingredients in these bounds are an estimate of the UV background from quasars, a temperature $T~10^4K$ for the "warm" photoionized IGM that produces most of the \lya absorption, a value of the Hubble constant, and observational estimates of the mean \lya flux decrement $\Dbar$ or, for a more restrictive bound, the distribution function $P(τ)$ of \lya optical depths. With plausible parameter values, the mean decrement bound implies a baryon density parameter $\Omb \ga 0.0125/h^2$. With conservative values, the bound weakens to $\Omb \ga 0.005/h^2$, but the required clustering of the IGM is then incompatible with other properties of quasar spectra. A recent observational determination of $P(τ)$ implies $\Omb \ga 0.0125/h^2$ even for a conservative estimate of the UV background, and $\Omb \ga 0.018/h^2$ for a more reasonable estimate. These bounds are consistent with recent low estimates of the primordial deuterium-to-hydrogen ratio $\dtoh$, which imply $\Omb \approx 0.025/h^2$, but the $P(τ)$ bound can only be reconciled with high $\dtoh$ estimates by abandoning standard big bang nucleosynthesis or the gravitational instability picture for the origin of the \lya forest. (Shortened abstract.)
△ Less
Submitted 5 January, 1997;
originally announced January 1997.
-
The Opacity of the Lyman Alpha Forest and Implications for Omega_{baryon} and the Ionizing Background
Authors:
M. Rauch,
J. Miralda-Escude,
W. L. W. Sargent,
T. A. Barlow,
D. H. Weinberg,
L. Hernquist,
N. Katz,
R. Cen,
J. P. Ostriker
Abstract:
We have measured the distribution function of the flux decrement caused by Lyman alpha forest absorption in a new sample of high resolution QSO spectra. The observations are compared to the results from two simulations of the Lya forest: an Eulerian Lambda-CDM model, and an SPH standard CDM model. Good agreement between the shapes of simulated and observed distributions is achieved by globally s…
▽ More
We have measured the distribution function of the flux decrement caused by Lyman alpha forest absorption in a new sample of high resolution QSO spectra. The observations are compared to the results from two simulations of the Lya forest: an Eulerian Lambda-CDM model, and an SPH standard CDM model. Good agreement between the shapes of simulated and observed distributions is achieved by globally scaling the optical depth to match the mean flux decrements. This procedure amounts to a measurement of the parameter Omega_b^2 h^3 / Gamma (where Omega_b is the baryonic matter density and Gamma is the HI ionization rate). Estimating a lower limit Gamma > 7 10^{-13} s^{-1} from the abundance of known QSOs, we derive a lower limit to the baryon density, Omega_b h^2>0.021(0.017) for the Lambda-CDM (SCDM) model. In both cases the large values are inconsistent with some recent D/H determinations (Rugers & Hogan 1996a,b), favoring a low deuterium abundance as reported by Tytler, Fan & Burles (1996). Adopting a fixed Omega_b, we can determine the evolution of the ionizing radiation field. Our models predict the intensity to be approximately constant with redshift, consistent with the assumption that the ionizing background is produced by known quasars for z < 3. However, additional sources of ionizing photons are required at higher redshift.
△ Less
Submitted 31 December, 1996;
originally announced December 1996.
-
Intergalactic Helium Absorption in Cold Dark Matter Models
Authors:
Rupert A. C. Croft,
David H. Weinberg,
Neal Katz,
Lars Hernquist
Abstract:
Observations from the HUT and the HST have recently detected HeII absorption along the lines of sight to two high redshift quasars. We use cosmological simulations with gas dynamics to investigate HeII absorption in the cold dark matter (CDM) theory of structure formation. We consider two Omega=1 CDM models with different normalizations and one Omega_0=0.4 CDM model, all incorporating the photoi…
▽ More
Observations from the HUT and the HST have recently detected HeII absorption along the lines of sight to two high redshift quasars. We use cosmological simulations with gas dynamics to investigate HeII absorption in the cold dark matter (CDM) theory of structure formation. We consider two Omega=1 CDM models with different normalizations and one Omega_0=0.4 CDM model, all incorporating the photoionizing UV background spectrum computed by Haardt & Madau (1996). The simulated gas distribution, combined with the H&M spectral shape, accounts for the relative observed values of taubar_HI and taubar_HeII, the effective mean optical depths for HI and HeII absorption. If the background intensity is as high as H&M predict, then matching the absolute values of taubar_HI and taubar_HeII requires a baryon abundance larger (by factors between 1.5 and 3 for the various CDM models) than our assumed value of Omega_b h^2=0.0125. The simulations reproduce the evolution of taubar_heII over the observed redshift range, 2.2 < z < 3.3, if the HeII photoionization rate remains roughly constant. HeII absorption in the CDM simulations is produced by a diffuse, fluctuating, intergalactic medium, which also gives rise to the HI ly-alpha forest. Much of the HeII opacity arises in underdense regions where the HI optical depth is very low. We compute statistical properties of the HeII and HI absorption that can be used to test the CDM models and distinguish them from an alternative scenario in which the HeII absorption is caused by discrete, compact clouds. The CDM scenario predicts that a substantial amount of baryonic material resides in underdense regions at high redshift. HeII absorption is the only sensitive probe of such extremely diffuse, intergalactic gas, so it can provide a vital test of this fundamental prediction.
△ Less
Submitted 26 May, 1997; v1 submitted 6 November, 1996;
originally announced November 1996.
-
Galaxy Structure, Dark Matter, and Galaxy Formation
Authors:
David H. Weinberg
Abstract:
The structure of galaxies, the nature of dark matter, and the physics of galaxy formation were the interlocking themes of DM 1996: Dark and Visible Matter in Galaxies and Cosmological Implications. In this conference summary report, I review recent observational and theoretical advances in these areas, then describe highlights of the meeting and discuss their implications. I include as an append…
▽ More
The structure of galaxies, the nature of dark matter, and the physics of galaxy formation were the interlocking themes of DM 1996: Dark and Visible Matter in Galaxies and Cosmological Implications. In this conference summary report, I review recent observational and theoretical advances in these areas, then describe highlights of the meeting and discuss their implications. I include as an appendix the lyrics of The Dark Matter Rap: A Cosmological History for the MTV Generation.
△ Less
Submitted 1 October, 1996;
originally announced October 1996.
-
The Population of Damped Lyman-alpha and Lyman Limit Systems in the Cold Dark Matter Model
Authors:
Jeffrey P. Gardner,
Neal Katz,
Lars Hernquist,
David H. Weinberg
Abstract:
Lyman limit and damped Lyman-alpha absorption systems probe the distribution of collapsed, cold gas at high redshift. Numerical simulations that incorporate gravity and gas dynamics can predict the abundance of such absorbers in cosmological models. We develop a semi-analytical method to correct the numerical predictions for the contribution of unresolved low mass halos, and we apply this method…
▽ More
Lyman limit and damped Lyman-alpha absorption systems probe the distribution of collapsed, cold gas at high redshift. Numerical simulations that incorporate gravity and gas dynamics can predict the abundance of such absorbers in cosmological models. We develop a semi-analytical method to correct the numerical predictions for the contribution of unresolved low mass halos, and we apply this method to the Katz et al. (1996) simulation of the standard cold dark matter model ($Ω=1$, $h=0.5$, $Ω_b=0.05$, $σ_8=0.7$). Using this simulation and higher resolution simulations of individual low mass systems, we determine the relation between a halo's circular velocity $v_c$ and its cross section for producing Lyman limit or damped absorption. We combine this relation with the Press-Schechter formula for the abundance of halos to compute the number of absorbers per unit redshift. The resolution correction increases the predicted abundances by about a factor of two at z=2, 3, and 4, bringing the predicted number of damped absorbers into quite good agreement with observations. Roughly half of the systems reside in halos with circular velocities $v_c>100\kms$ and half in halos with $35\kms<v_c<100\kms$. Halos with $v_c>150\kms$ typically harbor two or more systems capable of producing damped absorption. Even with the resolution correction, the predicted abundance of Lyman limit systems is a factor of three below observational estimates, signifying either a failure of standard CDM or a failure of these simulations to resolve the systems responsible for most Lyman limit absorption. By comparing simulations with and without star formation, we find that depletion of the gas supply by star formation affects absorption line statistics at $z>=2$ only for column densities exceeding $N_{HI}=10^{22} cm^{-2}$.
△ Less
Submitted 9 September, 1996;
originally announced September 1996.
-
Analysis of the Lyman-alpha Forest in Cosmological Simulations Using Voigt-Profile Decomposition
Authors:
Romeel Davé,
Lars Hernquist,
David H. Weinberg,
Neal Katz
Abstract:
We use an automated Voigt-profile fitting procedure to extract statistical properties of the Ly$α$ forest in a numerical simulation of an $Ω=1$, cold dark matter (CDM) universe. Our analysis method is similar to that used in most observational studies of the forest, and we compare the simulations to recently published results derived from Keck HIRES spectra. With the Voigt-profile decomposition…
▽ More
We use an automated Voigt-profile fitting procedure to extract statistical properties of the Ly$α$ forest in a numerical simulation of an $Ω=1$, cold dark matter (CDM) universe. Our analysis method is similar to that used in most observational studies of the forest, and we compare the simulations to recently published results derived from Keck HIRES spectra. With the Voigt-profile decomposition analysis, the simulation reproduces the large number of weak lines ($N_{\rm HI}\la 10^{14} cm^{-2}$) found in the HIRES spectra. At $z=3$, the $b$-parameter distribution has a median of $35\kms$ and a dispersion of $20\kms$, in reasonable agreement with the observed values. The comparison between our new analysis and recent data strengthens earlier claims that the Ly$alpha$ forest arises naturally in hierarchical structure formation as photoionized gas falls into dark matter potential wells.
△ Less
Submitted 6 September, 1996;
originally announced September 1996.
-
Photoionization, Numerical Resolution, and Galaxy Formation
Authors:
David H. Weinberg,
Lars Hernquist,
Neal Katz
Abstract:
Using cosmological simulations that incorporate gas dynamics and gravitational forces, we investigate the influence of photoionization by a UV radiation background on the formation of galaxies. In our highest resolution simulations, we find that photoionization has essentially no effect on the baryonic mass function of galaxies at $z=2$, down to our resolution limit of $5e9 M_\sun$. We do, howev…
▽ More
Using cosmological simulations that incorporate gas dynamics and gravitational forces, we investigate the influence of photoionization by a UV radiation background on the formation of galaxies. In our highest resolution simulations, we find that photoionization has essentially no effect on the baryonic mass function of galaxies at $z=2$, down to our resolution limit of $5e9 M_\sun$. We do, however, find a strong interplay between the mass resolution of a simulation and the microphysics included in the computation of heating and cooling rates. At low resolution, a photoionizing background can appear to suppress the formation of even relatively massive galaxies. However, when the same initial conditions are evolved with a factor of eight better mass resolution, this effect disappears. Our results demonstrate the need for care in interpreting the results of cosmological simulations that incorporate hydrodynamics and radiation physics. For example, we conclude that a simulation with limited resolution may yield more realistic results if it ignores some relevant physical processes, such as photoionization. At higher resolution, the simulated population of massive galaxies is insensitive to the treatment of photoionization and star formation, but it does depend significantly on the amplitude of the initial density fluctuations. By $z=2$, an $Ω=1$ cold dark matter model normalized to produce the observed masses of present-day clusters has already formed galaxies with baryon masses exceeding $1e11 M_\sun$.
△ Less
Submitted 29 April, 1996;
originally announced April 1996.
-
Imaging the Lyman-alpha Forest
Authors:
Andrew Gould,
David H. Weinberg
Abstract:
We show that it is now possible to image optically thick $\lya$ clouds in fluorescent $\lya$ emission with a relatively long ($\sim 20 $hr) integration on a large ($\sim 10 $m) telescope. For a broad range of column densities ($N\gsim 10^{18.5} \cm^{-2}$), the flux of $\lya$ photons from recombination cascades is equal to $\sim 0.6$ times the flux of ionizing photons, independent of the geometry…
▽ More
We show that it is now possible to image optically thick $\lya$ clouds in fluorescent $\lya$ emission with a relatively long ($\sim 20 $hr) integration on a large ($\sim 10 $m) telescope. For a broad range of column densities ($N\gsim 10^{18.5} \cm^{-2}$), the flux of $\lya$ photons from recombination cascades is equal to $\sim 0.6$ times the flux of ionizing photons, independent of the geometry of the cloud. Additional $\lya$ photons are produced by collisional excitations when these are the cloud's primary cooling mechanism. For typical physical conditions expected in optically thick clouds, these mechanisms together lead to a $\lya$ emission flux that is $\sim (2/3) \VEVν/ν_0$ times the flux of ionizing photons, where $\VEVν$ is the mean frequency of ionizing background photons and $ν_0$ is the Lyman limit frequency. Hence, measurement of the surface brightness from an optically thick cloud (known to exist, e.g., from a quasar absorption line) gives a direct measure of the energy in the ionizing radiation background. Moreover, in the same long slit spectrum one could hope to detect emission from $\sim 200$ other $\lya$ systems. Such detections would allow one to make a 2-dimensional map of the $\lya$ forest. By taking a series of such spectra, one could map the forest in three dimensions, revealing structure in the high-redshift universe.
△ Less
Submitted 21 December, 1995; v1 submitted 20 December, 1995;
originally announced December 1995.
-
Small Scale Structure and High Redshift HI
Authors:
David H. Weinberg,
Lars Hernquist,
Neal Katz,
Jordi Miralda-Escude
Abstract:
Cosmological simulations with gas dynamics suggest that the Lyman-alpha forest is produced mainly by "small scale structure" --- filaments and sheets that are the high redshift analog of today's galaxy superclusters. There is no sharp distinction between Lyman-alpha clouds and "Gunn-Peterson" absorption produced by the fluctuating IGM -- the Lyman-alpha forest {\it is} the Gunn-Peterson effect.…
▽ More
Cosmological simulations with gas dynamics suggest that the Lyman-alpha forest is produced mainly by "small scale structure" --- filaments and sheets that are the high redshift analog of today's galaxy superclusters. There is no sharp distinction between Lyman-alpha clouds and "Gunn-Peterson" absorption produced by the fluctuating IGM -- the Lyman-alpha forest {\it is} the Gunn-Peterson effect. Lyman limit and damped Lyman-alpha absorption arises in the radiatively cooled gas of forming galaxies. At $z~2-3$, most of the gas is in the photoionized, diffuse medium associated with the Lyman-alpha forest, but most of the {\it neutral} gas is in damped Lyman-alpha systems. We discuss generic evolution of cosmic gas in a hierarchical scenario of structure formation, with particular attention to the prospects for detecting 21cm emission from high redshift HI. A scaling argument based on the present-day cluster mass function suggests that objects with $M_{HI} >~ 5e11 h^{-1} \msun$ should be extremely rare at $z~3$, so detections with existing instruments will be difficult. An instrument like the proposed Square Kilometer Array could detect individual damped Lyman-alpha systems at high redshift, making it possible to map structure in the high redshift universe in much the same way that today's galaxy redshift surveys map the local large scale structure.
△ Less
Submitted 4 December, 1995;
originally announced December 1995.
-
Hydrodynamic Simulations of Galaxy Formation. II. Photoionization and the Formation of Low Mass Galaxies
Authors:
Anne A. Thoul,
David H. Weinberg
Abstract:
Photoionization by the high-redshift ultraviolet radiation background heats low density gas before it falls into dark matter potential wells, and it eliminates the neutral hydrogen and singly ionized helium that dominate cooling of primordial gas at temperatures of $10^4-10^5$K. We investigate the influence of photoionization on galaxy formation using high-resolution simulations with a 1-dimensi…
▽ More
Photoionization by the high-redshift ultraviolet radiation background heats low density gas before it falls into dark matter potential wells, and it eliminates the neutral hydrogen and singly ionized helium that dominate cooling of primordial gas at temperatures of $10^4-10^5$K. We investigate the influence of photoionization on galaxy formation using high-resolution simulations with a 1-dimensional, spherically symmetric, Lagrangian hydrodynamics/gravity code. We find that the presence of a photoionizing background suppresses the formation of galaxies with circular velocities $v_{circ} \simlt 30\,$\kms\ and substantially reduces the mass of cooled baryons in systems with circular velocities up to $v_{circ} \sim 50\,$\kms. Above $v_{circ} \sim 75\,$\kms, photoionization has no significant effect. Photoionization exerts its influence primarily by heating gas before collapse; the elimination of line cooling processes is less important. We discuss the implications of these results for hierarchical theories of galaxy formation.
△ Less
Submitted 31 October, 1995;
originally announced October 1995.
-
Cosmological Simulations with TreeSPH
Authors:
Neal Katz,
David H. Weinberg,
Lars Hernquist
Abstract:
We describe numerical methods for incorporating gas dynamics into cosmological simulations and present illustrative applications to the cold dark matter (CDM) scenario. Our evolution code, a version of TreeSPH (Hernquist \& Katz 1989) generalized to handle comoving coordinates and periodic boundary conditions, combines smoothed--particle hydrodynamics (SPH) with the hierarchical tree method for…
▽ More
We describe numerical methods for incorporating gas dynamics into cosmological simulations and present illustrative applications to the cold dark matter (CDM) scenario. Our evolution code, a version of TreeSPH (Hernquist \& Katz 1989) generalized to handle comoving coordinates and periodic boundary conditions, combines smoothed--particle hydrodynamics (SPH) with the hierarchical tree method for computing gravitational forces. The Lagrangian hydrodynamics approach and individual time steps for gas particles give the algorithm a large dynamic range, which is essential for studies of galaxy formation in a cosmological context. The code incorporates radiative cooling for an optically thin, primordial composition gas in ionization equilibrium with a user-specified ultraviolet background. We adopt a phenomenological prescription for star formation that gradually turns cold, dense, Jeans-unstable gas into collisionless stars, returning supernova feedback energy to the surrounding medium. In CDM simulations, some of the baryons that fall into dark matter potential wells dissipate their acquired thermal energy and condense into clumps with roughly galactic masses. The resulting galaxy population is insensitive to assumptions about star formation; we obtain similar baryonic mass functions and galaxy correlation functions from simulations with star formation and from simulations without star formation in which we identify galaxies directly from the cold, dense gas.
△ Less
Submitted 20 September, 1995;
originally announced September 1995.
-
Damped Lyman-alpha and Lyman Limit Absorbers in the Cold Dark Matter Model
Authors:
Neal Katz,
David H. Weinberg,
Lars Hernquist,
Jordi Miralda-Escude
Abstract:
We study the formation of damped \lya and Lyman limit absorbers in a hierarchical clustering scenario using a gas dynamical simulation of an $Ω= 1$, cold dark matter universe. In the simulation, these high column density systems are associated with forming galaxies. Damped \lya absorption, $N_{HI} \simgt 10^{20.2}\cm^{-2}$, arises along lines of sight that pass near the centers of relatively mas…
▽ More
We study the formation of damped \lya and Lyman limit absorbers in a hierarchical clustering scenario using a gas dynamical simulation of an $Ω= 1$, cold dark matter universe. In the simulation, these high column density systems are associated with forming galaxies. Damped \lya absorption, $N_{HI} \simgt 10^{20.2}\cm^{-2}$, arises along lines of sight that pass near the centers of relatively massive, dense protogalaxies. Lyman limit absorption, $10^{17}\cm^{-2} \simlt N_{HI} \simlt 10^{20.2}\cm^{-2}$, develops on lines of sight that pass through the outer parts of such objects or near the centers of smaller protogalaxies. The number of Lyman limit systems is less than observed, while the number of damped \lya systems is quite close to the observed abundance. Damped absorbers are typically $\sim 10$ kpc in radius, but the population has a large total cross section because the systems are much more numerous than present day $L_*$ galaxies. Our results demonstrate that high column density systems like those observed arise naturally in a hierarchical theory of galaxy formation and that it is now possible to study these absorbers directly from numerical simulations.
△ Less
Submitted 20 September, 1995;
originally announced September 1995.
-
The Lyman-alpha Forest in the Cold Dark Matter Model
Authors:
Lars Hernquist,
Neal Katz,
David H. Weinberg,
Jordi Miralda-Escude
Abstract:
Cosmological simulations with gas provide a detailed description of the intergalactic medium, making possible predictions of neutral hydrogen absorption in the spectra of background QSOs. We present results from a high-resolution calculation of an $Ω=1$ cold dark matter model. Our simulation reproduces many of the observed properties of the \lya forest surprisingly well. The distribution of HI c…
▽ More
Cosmological simulations with gas provide a detailed description of the intergalactic medium, making possible predictions of neutral hydrogen absorption in the spectra of background QSOs. We present results from a high-resolution calculation of an $Ω=1$ cold dark matter model. Our simulation reproduces many of the observed properties of the \lya forest surprisingly well. The distribution of HI column densities agrees with existing data to within a factor of $\sim$ two over most of the range from $10^{14}\cm^{-2}$ to $10^{22}\cm^{-2}$; i.e., from unsaturated \lya forest lines to damped \lya systems. The equivalent width distribution matches the observed exponential form with a characteristic width $W_* \approx 0.3$ angstroms. The distribution of b-parameters appears consistent with that derived from QSO spectra. Most of the low column density absorption arises in large, flattened structures of moderate or even relatively low overdensity, so there is no sharp distinction between the \lya forest and the ``Gunn-Peterson'' absorption produced by the smooth intergalactic medium. Our results demonstrate that a \lya forest like that observed develops naturally in a hierarchical clustering scenario with a photoionizing background. Comparison between simulations and high-resolution QSO spectra should open a new regime for testing theories of cosmic structure formation.
△ Less
Submitted 20 September, 1995;
originally announced September 1995.
-
The Sloan Digital Sky Survey
Authors:
James E. Gunn,
David H. Weinberg
Abstract:
We summarize the plans for and the current status of the Sloan Digital Sky Survey, a digital imaging and spectroscopic survey of $π$ steradians in the northern Galactic cap. The CCD photometric survey will produce images in five bands to limiting magnitudes of order 23. The spectroscopic survey will obtain redshifts of $10^6$ galaxies (a complete sample to a limiting magnitude $r' \sim 18$) and…
▽ More
We summarize the plans for and the current status of the Sloan Digital Sky Survey, a digital imaging and spectroscopic survey of $π$ steradians in the northern Galactic cap. The CCD photometric survey will produce images in five bands to limiting magnitudes of order 23. The spectroscopic survey will obtain redshifts of $10^6$ galaxies (a complete sample to a limiting magnitude $r' \sim 18$) and $10^5$ quasars ($g' \sim 19$). Repeated imaging of a 200 deg$^2$ strip in the southern Galactic cap will yield information about variable objects and a co-added photometric catalog roughly two magnitudes deeper than the northern survey. A dedicated 2.5-meter telescope, a large multi-CCD camera, and two fiber-fed double spectrographs are under construction and should be operational by fall of 1995. The main galaxy redshift sample will have a median redshift $\langle z \rangle \approx 0.1$.
△ Less
Submitted 20 December, 1994;
originally announced December 1994.
-
Constraints onOmega from the IRAS Redshift Surveys
Authors:
Shaun Cole,
Karl B. Fisher,
David H. Weinberg
Abstract:
We measure the anisotropy of the redshift-space power spectrum in the 1.2-Jy and QDOT redshift surveys of IRAS-selected galaxies. On large scales, this anisotropy is caused by coherent peculiar motions, and gravitational instability theory predicts a distortion of the power spectrum that depends only on the ratio $β\equiv f(Ω)/b \approx Ω^{0.6}/b$, where Omega is the cosmological density paramet…
▽ More
We measure the anisotropy of the redshift-space power spectrum in the 1.2-Jy and QDOT redshift surveys of IRAS-selected galaxies. On large scales, this anisotropy is caused by coherent peculiar motions, and gravitational instability theory predicts a distortion of the power spectrum that depends only on the ratio $β\equiv f(Ω)/b \approx Ω^{0.6}/b$, where Omega is the cosmological density parameter and $b$ is the bias parameter. On small scales, the distortion is dominated by the random velocity dispersion in non-linear structures. We fit the observed anisotropy with an analytic model that incorporates two parameters, beta, and a small-scale velocity dispersion sigma_v. Tests on N-body simulations show that this model recovers beta quite accurately on the scales accessible to the existing IRAS redshift surveys.
Applying our procedure to the 1.2-Jy and QDOT surveys, we find beta=0.52 +/- 0.13 and beta=0.54 +/- 0.3, respectively. These results imply Omega approximately 0.35 if galaxies trace mass, or a bias factor of about 2 if Omega=1.
△ Less
Submitted 17 December, 1994;
originally announced December 1994.
-
Hydrodynamic Simulations of Galaxy Formation. I. Dissipation and the Maximum Mass of Galaxies
Authors:
A. A. Thoul,
D. H. Weinberg
Abstract:
We describe an accurate, one-dimensional, spherically symmetric, Lagrangian hydrodynamics/gravity code, designed to study the effects of radiative cooling and photo-ionization on the formation of protogalaxies. The code can treat an arbitrary number of fluid shells (representing baryons) and collisionless shells (representing cold dark matter). As a test of the code, we reproduce analytic soluti…
▽ More
We describe an accurate, one-dimensional, spherically symmetric, Lagrangian hydrodynamics/gravity code, designed to study the effects of radiative cooling and photo-ionization on the formation of protogalaxies. The code can treat an arbitrary number of fluid shells (representing baryons) and collisionless shells (representing cold dark matter). As a test of the code, we reproduce analytic solutions for the pulsation behavior of a polytrope and for the self-similar collapse of a spherically symmetric, cosmological perturbation. In this paper, we concentrate on the effects of radiative cooling, examining the ability of collapsing perturbations to cool within the age of the universe. In contrast to some studies based on order-of- magnitude estimates, we find that cooling arguments alone cannot explain the sharp upper cutoff observed in the galaxy luminosity function.
△ Less
Submitted 9 December, 1994; v1 submitted 3 October, 1994;
originally announced October 1994.
-
Studying Structure Formation With Large Galaxy Redshift Surveys
Authors:
David H. Weinberg
Abstract:
I outline the connections between some of the most widely used statistical measures of galaxy clustering and the fundamental issues in the theory of structure formation. I devote particular attention to the problem of biasing, i.e. to a possible difference between the distribution of galaxies and the distribution of mass. Using numerical experiments, I show that a local morphology-density relati…
▽ More
I outline the connections between some of the most widely used statistical measures of galaxy clustering and the fundamental issues in the theory of structure formation. I devote particular attention to the problem of biasing, i.e. to a possible difference between the distribution of galaxies and the distribution of mass. Using numerical experiments, I show that a local morphology-density relation leads to different slopes and amplitudes for the correlation functions of different galaxy types on small scales, but that the correlation functions on large scales all have the same shape, differing only by constant factors. I also examine a variety of biasing models in which the efficiency of galaxy formation depends on local properties of the mass distribution. While these purely local influences do lead to a large-scale bias, in all cases the bias factor becomes constant on scales larger than the galaxy correlation length.
△ Less
Submitted 29 September, 1994;
originally announced September 1994.
-
A Test of the Adhesion Approximation for Gravitational Clustering
Authors:
A. L. Melott,
S. F. Shandarin,
D. H. Weinberg
Abstract:
We quantitatively compare a particle implementation of the adhesion approximation to fully non--linear, numerical nbody simulations. Our primary tool, cross--correlation of nbody simulations with the adhesion approximation, indicates good agreement, better than that found by the same test performed with the Zel'dovich approximation (hereafter ZA). However, the cross--correlation is not as good a…
▽ More
We quantitatively compare a particle implementation of the adhesion approximation to fully non--linear, numerical nbody simulations. Our primary tool, cross--correlation of nbody simulations with the adhesion approximation, indicates good agreement, better than that found by the same test performed with the Zel'dovich approximation (hereafter ZA). However, the cross--correlation is not as good as that of the truncated Zel'dovich approximation (TZA), obtained by applying the Zel'dovich approximation after smoothing the initial density field with a Gaussian filter. We confirm that the adhesion approximation produces an excessively filamentary distribution. Relative to the nbody results, we also find that: (a) the power spectrum obtained from the adhesion approximation is more accurate than that from ZA or TZA, (b) the error in the phase angle of Fourier components is worse than that from TZA, and (c) the mass distribution function is more accurate than that from ZA or TZA. It appears that adhesion performs well statistically, but that TZA is more accurate dynamically, in the sense of moving mass to the right place.
Subject Heading: Galaxies, formation, clustering--large--scale structure of the Universe
△ Less
Submitted 30 November, 1993;
originally announced November 1993.
-
Fourier Analysis of Redshift Space Distortions and the Determination of Omega
Authors:
Shaun Cole,
Karl B. Fisher,
David H. Weinberg
Abstract:
The peculiar velocities of galaxies distort the pattern of galaxy clustering in redshift space, making the redshift space power spectrum anisotropic. In the linear regime, the strength of this distortion depends only on the ratio $β\equiv f(Ω)/b \approx Ω^{0.6}/b$, where $Ω$ is the cosmological density parameter and $b$ is the bias parameter. We derive a linear theory estimator for $β$ based on…
▽ More
The peculiar velocities of galaxies distort the pattern of galaxy clustering in redshift space, making the redshift space power spectrum anisotropic. In the linear regime, the strength of this distortion depends only on the ratio $β\equiv f(Ω)/b \approx Ω^{0.6}/b$, where $Ω$ is the cosmological density parameter and $b$ is the bias parameter. We derive a linear theory estimator for $β$ based on the harmonic moments of the redshift space power spectrum. Using N-body simulations, we examine the impact of non-linear gravitational clustering on the power spectrum anisotropy and on our $β$-estimator. Non-linear effects can be important out to wavelengths $λ\sim 50$Mpc/h or larger; in most cases, they lower the quadrupole moment of the power spectrum and thereby depress the estimate of $β$ below the true value. With a sufficiently large redshift survey, the scaling of non-linear effects may allow separate determinations of $Ω$ and $b$.
We describe a practical technique for measuring the anisotropy of the power spectrum from galaxy redshift surveys, and we test the technique on mock catalogues drawn from the N-body simulations. Preliminary application of our methods to the 1.2 Jy IRAS galaxy survey yields $β_{est} \sim 0.3-0.4 $ at wavelengths $λ\sim 30-40$Mpc/h . Non-linear effects remain important at these scales, so this estimate of $β$ is probably lower than the true value.
△ Less
Submitted 4 August, 1993;
originally announced August 1993.
-
Cosmic Voids and Biased Galaxy Formation
Authors:
Blane Little,
David H. Weinberg
Abstract:
Using cosmological $N$-body simulations and the void probability function (VPF), we investigate the statistical properties of voids within a wide range of initially Gaussian models for the origin of large-scale structure. We pay particular attention to the ability of the VPF to diagnose `biased galaxy formation': the preferential formation of galaxies in regions of high background density and co…
▽ More
Using cosmological $N$-body simulations and the void probability function (VPF), we investigate the statistical properties of voids within a wide range of initially Gaussian models for the origin of large-scale structure. We pay particular attention to the ability of the VPF to diagnose `biased galaxy formation': the preferential formation of galaxies in regions of high background density and corresponding suppression of galaxy formation in regions of low background density. Biased models have systematically higher VPFs than unbiased models, but for a given biasing scheme the VPF is relatively insensitive to the value of the bias factor $b$. Thus, while the VPF can distinguish unbiased models from some biased models, it is probably not a useful way to constrain the bias factor; uncertainties in the appropriate choice of biasing prescription overwhelm the mild dependence on $b$. We compare the predictions of our models to the most extensive VPF observations published to date. These data do {\it not} require strong biasing; Gaussian models in which galaxies trace mass can reproduce the VPF data to within the errors expected from the current finite volume fluctuations. Models with the moderate biasing predicted by cosmological simulations with gas dynamics yield a slightly better match to the data. Models in which galaxy formation is strongly suppressed in low density regions produce an excess of large, empty voids.
△ Less
Submitted 11 June, 1993;
originally announced June 1993.