-
Galactic Halo Cusp versus Core: Tidal Effects in Mergers
Authors:
Avishai Dekel,
Jonathan Devor,
Itai Arad
Abstract:
We show how the buildup of halos by merging satellites forces an inner cusp, with a density profile ρ\propto r^{-α} where α\to α_a \gsim 1$. Our analysis is based on a new prescription for tidal stripping as a function of α(r), using a simple toy model which matches N-body simulations. In a core of α< 1 there is tidal compression rather than stripping and the satellites sink towards the halo cen…
▽ More
We show how the buildup of halos by merging satellites forces an inner cusp, with a density profile ρ\propto r^{-α} where α\to α_a \gsim 1$. Our analysis is based on a new prescription for tidal stripping as a function of α(r), using a simple toy model which matches N-body simulations. In a core of α< 1 there is tidal compression rather than stripping and the satellites sink towards the halo center, causing a rapid steepening of the profile to α> 1. Where α> 1, the stripping of each satellite shell is preceded by gradual puffing up, which makes the stripping more efficient at larger α, causing flattening where αis large enough. Therefore, we can show using linear perturbation analysis that a sequence of mergers slowly leads to a fixed point α(r)=α_a. This result implies that a cusp is enforced as long as enough satellite material makes it into the inner halo and is deposited there. We conclude that in order to maintain a flat core, as indicated by observations, satellites must be disrupted outside the core, e.g., because of puffing up due to baryonic feedback effects.
△ Less
Submitted 14 January, 2002;
originally announced January 2002.
-
Towards a Resolution of the Galactic Spin Crisis: Mergers, Feedback, and Spin Segregation
Authors:
Ariyeh H. Maller,
Avishai Dekel
Abstract:
We model in simple terms the angular-momentum problems of galaxy formation in CDM cosmologies, and identify the key elements of a scenario that may solve them. The buildup of angular momentum is modeled via dynamical friction and tidal stripping in a sequence of mergers. We demonstrate how over-cooling in incoming halos leads to a transfer of angular momentum from the baryons to the dark matter,…
▽ More
We model in simple terms the angular-momentum problems of galaxy formation in CDM cosmologies, and identify the key elements of a scenario that may solve them. The buildup of angular momentum is modeled via dynamical friction and tidal stripping in a sequence of mergers. We demonstrate how over-cooling in incoming halos leads to a transfer of angular momentum from the baryons to the dark matter, in conflict with observations. By incorporating a simple recipe of supernova feedback, we are able to solve the problems of angular momentum in disk formation. Gas removal from the numerous small incoming halos, which merge to become the low specific angular momentum (j) component of the product, eliminates the low-j baryons. Heating and puffing-up of the gas in larger incoming halos, combined with efficient tidal stripping, reduces the angular momentum loss of baryons due to dynamical friction. Dependence of the feedback effects on the progenitor halo mass implies that the spin of baryons is typically higher for lower mass halos. The observed low baryonic fraction in dwarf galaxies is used to calibrate the characteristic velocity associated with supernova feedback, yielding ~100 km/s, within the range of theoretical expectations. We then find that the model naturally produces the observed distribution of the spin parameter among dwarf and bright disk galaxies, as well as the j profile inside these galaxies. This suggests that the model indeed captures the main features of a full scenario for resolving the spin crisis.
△ Less
Submitted 11 January, 2002;
originally announced January 2002.
-
Simulating the Formation of the Local Galaxy Population
Authors:
H. Mathis,
G. Lemson,
V. Springel,
G. Kauffmann,
S. D. M. White,
A. Eldar,
A. Dekel
Abstract:
We simulate the formation and evolution of the local galaxy population starting from initial conditions with a smoothed linear density field which matches that derived from the IRAS 1.2 Jy galaxy survey. Our simulations track the formation and evolution of all dark matter haloes more massive than 10e+11 solar masses out to a distance of 8000 km/s from the Milky Way. We implement prescriptions si…
▽ More
We simulate the formation and evolution of the local galaxy population starting from initial conditions with a smoothed linear density field which matches that derived from the IRAS 1.2 Jy galaxy survey. Our simulations track the formation and evolution of all dark matter haloes more massive than 10e+11 solar masses out to a distance of 8000 km/s from the Milky Way. We implement prescriptions similar to those of Kauffmann et al. (1999) to follow the assembly and evolution of the galaxies within these haloes. We focus on two variants of the CDM cosmology: an LCDM and a tCDM model. Galaxy formation in each is adjusted to reproduce the I-band Tully-Fisher relation of Giovanelli et al. (1997). We compare the present-day luminosity functions, colours, morphology and spatial distribution of our simulated galaxies with those of the real local population, in particular with the Updated Zwicky Catalog, with the IRAS PSCz redshift survey, and with individual local clusters such as Coma, Virgo and Perseus. We also use the simulations to study the clustering bias between the dark matter and galaxies of differing type. Although some significant discrepancies remain, our simulations recover the observed intrinsic properties and the observed spatial distribution of local galaxies reasonably well. They can thus be used to calibrate methods which use the observed local galaxy population to estimate the cosmic density parameter or to draw conclusions about the mechanisms of galaxy formation. To facilitate such work, we publically release our z=0 galaxy catalogues, together with the underlying mass distribution.
△ Less
Submitted 5 November, 2001;
originally announced November 2001.
-
Concentrations and Assembly Histories of Dark Matter Halos
Authors:
Risa H. Wechsler,
James S. Bullock,
J. R. Primack,
Andrey V. Kravtsov,
Avishai Dekel
Abstract:
We study the relation between the density profiles of dark matter halos and their mass assembly histories, using a statistical sample of halos in a high-resolution N-body simulation of the LCDM cosmology. For each halo at z=0, we identify its merger-history tree, and determine concentration parameters c_vir for all progenitors, thus providing a structural merger tree for each halo. We fit the ma…
▽ More
We study the relation between the density profiles of dark matter halos and their mass assembly histories, using a statistical sample of halos in a high-resolution N-body simulation of the LCDM cosmology. For each halo at z=0, we identify its merger-history tree, and determine concentration parameters c_vir for all progenitors, thus providing a structural merger tree for each halo. We fit the mass accretion histories by a universal function with one parameter, the formation epoch a_c, defined when the log mass accretion rate dlog M/doga falls below a critical value. We find that late forming galaxies tend to be less concentrated, such that c_vir ``observed'' at any epoch a_obs is strongly correlated with a_c via c_vir=c_1 a_obs/a_c. Scatter about this relation is mostly due to measurement errors in c_vir and a_c, implying that the actual spread in c_vir for halos of a given mass can be mostly attributed to scatter in a_c. Because of the direct connection between halo concentration and velocity rotation curves, and because of probable connections between halo mass assembly history and star formation history, the tight correlation between these properties provides an essential new ingredient for galaxy formation modeling.
△ Less
Submitted 4 November, 2001;
originally announced November 2001.
-
Concentrations of Dark Halos from their Assembly Histories
Authors:
Risa H. Wechsler,
James S. Bullock,
Joel R. Primack,
Andrey V. Kravtsov,
Avishai Dekel
Abstract:
(abridged) We study the relation between the density profiles of dark matter halos and their mass assembly histories, using a statistical sample of halos in a high-resolution N-body simulation of the LCDM cosmology. For each halo at z=0, we identify its merger-history tree, and determine concentration parameters c_vir for all progenitors, thus providing a structural merger tree for each halo. We…
▽ More
(abridged) We study the relation between the density profiles of dark matter halos and their mass assembly histories, using a statistical sample of halos in a high-resolution N-body simulation of the LCDM cosmology. For each halo at z=0, we identify its merger-history tree, and determine concentration parameters c_vir for all progenitors, thus providing a structural merger tree for each halo. We fit the mass accretion histories by a universal function with one parameter, the formation epoch a_c, defined when the log mass accretion rate dlogM/dloga falls below a critical value S. We find that late forming galaxies tend to be less concentrated, such that c_vir ``observed'' at any epoch a_o is strongly correlated with a_c via c_vir=c_1*a_o/a_c. Scatter about this relation is mostly due to measurement errors in c_v and a_c, implying that the actual spread in c_vir for halos of a given mass can be mostly attributed to scatter in a_c. We demonstrate that this relation can also be used to predict the mass and redshift dependence of c_v, and the scatter about the median c_vir(M,z), using accretion histories derived from the Extended Press-Schechter (EPS) formalism, after adjusting for a constant offset between the formation times as predicted by EPS and as measured in the simulations;this new ingredient can thus be easily incorporated into semi-analytic models of galaxy formation. The correlation found between halo concentration and mass accretion rate suggests a physical interpretation: for high mass infall rates the central density is related to the background density; when the mass infall rate slows, the central density stays approximately constant and the halo concentration just grows as R_vir. The tight correlation demonstrated here provides an essential new ingredient for galaxy formation modeling.
△ Less
Submitted 4 November, 2001; v1 submitted 9 August, 2001;
originally announced August 2001.
-
Cosmic Flows: A Status Report
Authors:
Stephane Courteau,
Avishai Dekel
Abstract:
We give a brief review of recent developments in the study of the large-scale velocity field of galaxies since the international workshop on Cosmic Flows held in July 1999 in Victoria, B.C. Peculiar velocities (PVs) yield a tight and unique constraint on cosmological characteristics, independent of Lambda and biasing, such as the cosmological matter density parameter Omega_m and the convergence…
▽ More
We give a brief review of recent developments in the study of the large-scale velocity field of galaxies since the international workshop on Cosmic Flows held in July 1999 in Victoria, B.C. Peculiar velocities (PVs) yield a tight and unique constraint on cosmological characteristics, independent of Lambda and biasing, such as the cosmological matter density parameter Omega_m and the convergence of bulk flows on large scales. Significant progress towards incorporating non-linear dynamics and improvements of velocity field reconstruction techniques have led to a rigorous control of errors and much refined cosmic flow analyses. Current investigations favor low-amplitude (< 250 km/s) bulk flows on the largest scales (< ~100 Mpc/h) probed reliably by existing redshift-distance surveys, consistent with favored LambdaCDM cosmogonies. Tidal field analyses also suggest that the Shapley Concentration, located behind the Great Attractor (GA), might play an important dynamical role, even at the Local Group. Low-amplitude density fluctuations on very large scales generate the overall large-scale streaming motions while massive attractors like the GA, and Perseus-Pisces account for smaller scale motions which are superposed on the large-scale flow. Likelihood analyses of galaxy PVs, in the framework of flat CDM cosmology, now provide tight constraints of Omega_m = 0.35 +/- 0.05. A four-fold size increase of our data base is expected in the next 4-5 years with the completion of next generation FP/TF surveys and automated supernovae searches within 20,000 km/s.
△ Less
Submitted 3 June, 2001; v1 submitted 27 May, 2001;
originally announced May 2001.
-
Modeling Angular-Momentum History in Dark-Matter Halos
Authors:
Ariyeh H. Maller,
Avishai Dekel,
Rachel S. Somerville
Abstract:
We model the acquisition of spin by dark-matter halos in semi-analytic merger trees. We explore two different algorithms; one in which halo spin is acquired from the orbital angular momentum of merging satellites, and another in which halo spin is gained via tidal torquing on shells of material while still in the linear regime. We find that both scenarios produce the characteristic spin distribu…
▽ More
We model the acquisition of spin by dark-matter halos in semi-analytic merger trees. We explore two different algorithms; one in which halo spin is acquired from the orbital angular momentum of merging satellites, and another in which halo spin is gained via tidal torquing on shells of material while still in the linear regime. We find that both scenarios produce the characteristic spin distribution of halos found in N-body simulations, namely, a log-normal distribution with mean ~0.04 and standard deviation ~0.5 in the log. A perfect match requires fine-tuning of two free parameters. Both algorithms also reproduce the general insensitivity of the spin distribution to halo mass, redshift and cosmology seen in N-body simulations. The spin distribution can be made strictly constant by physically motivated scalings of the free parameters. In addition, both schemes predict that halos which have had recent major mergers have systematically larger spin values. These algorithms can be implemented within semi-analytic models of galaxy formation based on merger trees. They yield detailed predictions of galaxy properties that strongly depend on angular momentum (such as size and surface brightness) as a function of merger history and environment.
△ Less
Submitted 30 May, 2001; v1 submitted 10 May, 2001;
originally announced May 2001.
-
Testing Tidal-Torque Theory: II. Alignment of Inertia and Shear and the Characteristics of Proto-haloes
Authors:
Cristiano Porciani,
Avishai Dekel,
Yehuda Hoffman
Abstract:
We investigate the cross-talk between the two key components of tidal-torque theory, the inertia (I) and shear (T) tensors, using a cosmological N-body simulation with thousands of well-resolved haloes. We find that the principal axes of I and T are strongly aligned, even though I characterizes the proto-halo locally while T is determined by the large-scale structure. Thus, the resultant galacti…
▽ More
We investigate the cross-talk between the two key components of tidal-torque theory, the inertia (I) and shear (T) tensors, using a cosmological N-body simulation with thousands of well-resolved haloes. We find that the principal axes of I and T are strongly aligned, even though I characterizes the proto-halo locally while T is determined by the large-scale structure. Thus, the resultant galactic spin, which plays a key role in galaxy formation, is only a residual due to ~10 per cent deviations from perfect alignment of T and I. The T-I correlation induces a weak tendency for the proto-halo spin to be perpendicular to the major axes of T and I, but this correlation is erased by non-linear effects at late times, making the observed spins poor indicators of the initial shear field. However, the T-I correlation implies that the shear tensor can be used for identifying the positions and boundaries of proto-haloes in cosmological initial conditions -- a missing piece in galaxy formation theory. The typical configuration is of a prolate proto-halo lying perpendicular to a large-scale high-density ridge, with the surrounding voids inducing compression along the major and intermediate inertia axes of the proto-halo. This leads to a transient sub-halo filament along the large-scale ridge, whose sub-clumps then flow along the filament and merge into the final halo. The centres of proto-haloes tend to lie in ~1 sigma over-density regions, but their association with linear density maxima smoothed on galactic scales is vague: only ~40 per cent of the proto-haloes contain peaks within them. Several other characteristics distinguish proto-haloes from density peaks, e.g., they tend to compress along two principal axes while many peaks compress along three axes.
△ Less
Submitted 12 January, 2002; v1 submitted 10 May, 2001;
originally announced May 2001.
-
Testing Tidal-Torque Theory: I. Spin Amplitude and Direction
Authors:
Cristiano Porciani,
Avishai Dekel,
Yehuda Hoffman
Abstract:
We evaluate the success of linear tidal-torque theory (TTT) in predicting galactic-halo spin using a cosmological N-body simulation with thousands of well-resolved haloes. The proto-haloes are identified by tracing today's haloes back to the initial conditions. The TTT predictions for the proto-haloes match, on average, the spin amplitudes of today's virialized haloes if linear growth is assumed…
▽ More
We evaluate the success of linear tidal-torque theory (TTT) in predicting galactic-halo spin using a cosmological N-body simulation with thousands of well-resolved haloes. The proto-haloes are identified by tracing today's haloes back to the initial conditions. The TTT predictions for the proto-haloes match, on average, the spin amplitudes of today's virialized haloes if linear growth is assumed until ~t_0/3, or 55-70 per cent of the halo effective turn-around time. This makes it a useful qualitative tool for understanding certain average properties of galaxies, such as total spin and angular-momentum distribution within haloes, but with a random scatter of the order of the signal itself. Non-linear changes in spin direction cause a mean error of ~50 degrees in the TTT prediction at t_0, such that the linear spatial correlations of spins on scales >~ 1 Mpc/h are significantly weakened by non-linear effects. This questions the usefulness of TTT for predicting intrinsic alignments in the context of gravitational lensing. We find that the standard approximations made in TTT, including a second order expansion of the Zel'dovich potential and a smoothing of the tidal field, provide close-to-optimal results.
△ Less
Submitted 12 January, 2002; v1 submitted 8 May, 2001;
originally announced May 2001.
-
The Large-Scale Tidal Velocity Field
Authors:
Y. Hoffman,
A. Eldar,
S. Zaroubi,
A. Dekel
Abstract:
We present a method for decomposing the cosmological velocity field in a given volume into its divergent component due to the density fluctuations inside the volume, and its tidal component due to the matter distribution outside the volume. The input consists of the density and velocity fields that are reconstructed either by POTENT or by Wiener Filter from a survey of peculiar velocities. The t…
▽ More
We present a method for decomposing the cosmological velocity field in a given volume into its divergent component due to the density fluctuations inside the volume, and its tidal component due to the matter distribution outside the volume. The input consists of the density and velocity fields that are reconstructed either by POTENT or by Wiener Filter from a survey of peculiar velocities. The tidal field is further decomposed into a bulk velocity and a shear field. The method is applied here to the Mark III data within a sphere of radius 60 Mpc/h about the Local Group, and to the SFI data for comparison. We find that the tidal field contributes about half of the Local-Group velocity with respect to the CMB, with the tidal bulk velocity pointing to within ~ 30 degrees of the CMB dipole. The eigenvector with the largest eigenvalue of the shear tensor is aligned with the tidal bulk velocity to within ~ 40 degrees. The tidal field thus indicates the important dynamical role of a super attractor of mass (2-5) x 10^17 M_sun/h Omega^0.4 at ~ 150 Mpc/h, coinciding with the Shapley Concentration. There is also a hint for the dynamical role of two big voids in the Supergalactic Plane. The results are consistent for the two data sets and the two methods of reconstruction.
△ Less
Submitted 11 February, 2001;
originally announced February 2001.
-
Nonlinear Peculiar-Velocity Analysis and PCA
Authors:
A. Dekel,
A. Eldar,
L. Silberman,
I. Zehavi
Abstract:
We allow for nonlinear effects in the likelihood analysis of peculiar velocities, and obtain ~35%-lower values for the cosmological density parameter and for the amplitude of mass-density fluctuations. The power spectrum in the linear regime is assumed to be of the flat LCDM model (h=0.65, n=1) with only Om_m free. Since the likelihood is driven by the nonlinear regime, we "break" the power spec…
▽ More
We allow for nonlinear effects in the likelihood analysis of peculiar velocities, and obtain ~35%-lower values for the cosmological density parameter and for the amplitude of mass-density fluctuations. The power spectrum in the linear regime is assumed to be of the flat LCDM model (h=0.65, n=1) with only Om_m free. Since the likelihood is driven by the nonlinear regime, we "break" the power spectrum at k_b=0.2 h/Mpc and fit a two-parameter power-law at k>k_b. This allows for an unbiased fit in the linear regime. Tests using improved mock catalogs demonstrate a reduced bias and a better fit. We find for the Mark III and SFI data Om_m=0.35+-0.09$ with sigma_8*Om_m^0.6=0.55+-0.10 (90% errors). When allowing deviations from \lcdm, we find an indication for a wiggle in the power spectrum in the form of an excess near k~0.05 and a deficiency at k~0.1 h/Mpc --- a "cold flow" which may be related to a feature indicated from redshift surveys and the second peak in the CMB anisotropy. A chi^2 test applied to principal modes demonstrates that the nonlinear procedure improves the goodness of fit. The Principal Component Analysis (PCA) helps identifying spatial features of the data and fine-tuning the theoretical and error models. We address the potential for optimal data compression using PCA.
△ Less
Submitted 28 January, 2001;
originally announced January 2001.
-
Cosmological Density and Power Spectrum from Peculiar Velocities: Nonlinear Corrections and PCA
Authors:
L. Silberman,
A. Dekel,
A. Eldar,
I. Zehavi
Abstract:
We allow for nonlinear effects in the likelihood analysis of galaxy peculiar velocities, and obtain ~35%-lower values for the cosmological density parameter Om and the amplitude of mass-density fluctuations. The power spectrum in the linear regime is assumed to be a flat LCDM model (h=0.65, n=1, COBE) with only Om as a free parameter. Since the likelihood is driven by the nonlinear regime, we "b…
▽ More
We allow for nonlinear effects in the likelihood analysis of galaxy peculiar velocities, and obtain ~35%-lower values for the cosmological density parameter Om and the amplitude of mass-density fluctuations. The power spectrum in the linear regime is assumed to be a flat LCDM model (h=0.65, n=1, COBE) with only Om as a free parameter. Since the likelihood is driven by the nonlinear regime, we "break" the power spectrum at k_b=0.2 h/Mpc and fit a power law at k>k_b. This allows for independent matching of the nonlinear behavior and an unbiased fit in the linear regime. The analysis assumes Gaussian fluctuations and errors, and a linear relation between velocity and density. Tests using proper mock catalogs demonstrate a reduced bias and a better fit. We find for the Mark3 and SFI data Om_m=0.32+-0.06 and 0.37+-0.09 respectively, with sigma_8*Om^0.6 = 0.49+-0.06 and 0.63+-0.08, in agreement with constraints from other data. The quoted 90% errors include cosmic variance. The improvement in likelihood due to the nonlinear correction is very significant for Mark3 and moderately so for SFI. When allowing deviations from LCDM, we find an indication for a wiggle in the power spectrum: an excess near k=0.05 and a deficiency at k=0.1 (cold flow). This may be related to the wiggle seen in the power spectrum from redshift surveys and the second peak in the CMB anisotropy. A chi^2 test applied to modes of a Principal Component Analysis (PCA) shows that the nonlinear procedure improves the goodness of fit and reduces a spatial gradient of concern in the linear analysis. The PCA allows addressing spatial features of the data and fine-tuning the theoretical and error models. It shows that the models used are appropriate for the cosmological parameter estimation performed. We address the potential for optimal data compression using PCA.
△ Less
Submitted 24 April, 2001; v1 submitted 21 January, 2001;
originally announced January 2001.
-
Galaxy Formation at z~3: Constraints from Spatial Clustering
Authors:
Risa H. Wechsler,
Rachel S. Somerville,
James S. Bullock,
Tsafrir S. Kolatt,
Joel R. Primack,
George R. Blumenthal,
Avishai Dekel
Abstract:
We use N-body simulations combined with semi-analytic models to compute the clustering properties of modeled galaxies at z~3, and confront these predictions with the clustering properties of the observed population of Lyman-break galaxies (LBGs). Several scenarios for the nature of LBGs are explored, which may be broadly categorized into models in which high-redshift star formation is driven by…
▽ More
We use N-body simulations combined with semi-analytic models to compute the clustering properties of modeled galaxies at z~3, and confront these predictions with the clustering properties of the observed population of Lyman-break galaxies (LBGs). Several scenarios for the nature of LBGs are explored, which may be broadly categorized into models in which high-redshift star formation is driven by collisional starbursts and those in which quiescent star formation dominates. For each model, we make predictions for the LBG overdensity distribution, the variance of counts-in-cells, the correlation length, and close pair statistics. Models which assume a one-to-one relationship between massive dark-matter halos and galaxies are disfavored by close pair statistics, as are models in which colliding halos are associated with galaxies in a simplified way. However, when modeling of gas consumption and star formation is included using a semi-analytic treatment, the quiescent and collisional starburst models predict similar clustering properties and none of these models can be ruled out based on the available clustering data. None of the ``realistic'' models predict a strong dependence of clustering amplitude on the luminosity threshold of the sample, in apparent conflict with some observational results.
△ Less
Submitted 25 January, 2001; v1 submitted 14 November, 2000;
originally announced November 2000.
-
Spin Profile of Galactic Halos and Disk Formation
Authors:
A. Dekel,
J. S. Bullock,
C. Porciani,
A. V. Kravtsov,
T. S. Kolatt,
A. A. Klypin,
J. R. Primack
Abstract:
We summarize recent developments in the study of the origin of halo spin profiles and preliminary implications on disk formation. The specific angular-momentum distributions within halos in N-body simulations match a universal shape, M(<j) \propto j/(j_0+j). It is characterized by a power law over most of the mass, and one shape parameter in addition to the spin parameter λ. The angular momentum…
▽ More
We summarize recent developments in the study of the origin of halo spin profiles and preliminary implications on disk formation. The specific angular-momentum distributions within halos in N-body simulations match a universal shape, M(<j) \propto j/(j_0+j). It is characterized by a power law over most of the mass, and one shape parameter in addition to the spin parameter λ. The angular momentum tends to be aligned throughout the halo and of cylindrical symmetry. Even if angular momentum is conserved during baryonic infall, the resultant disk density profile is predicted to deviate from exponential, with a denser core and an extended tail. A slightly corrected version of the scaling relation due to linear tidal-torque theory is used to explain the origin of a typical power-law profile in shells, j(M) \propto M^s with s \gsim 1. While linear theory crudely predicts the amplitudes of halo spins, it is not a good predictor of their directions. Independently, mergers of halos are found to produce a similar profile due to j transfer from the orbit to the product halo via dynamical friction and tidal stripping. The halo spin is correlated with having a recent major merger, though this correlation is weakened by mass loss. These two effects are due to a correlation between the spins of neighboring halos and their orbit, leading to prograde mergers.
△ Less
Submitted 31 October, 2000;
originally announced November 2000.
-
A universal angular momentum profile for galactic halos
Authors:
James S. Bullock,
Avishai Dekel,
Tsafrir S. Kolatt,
Andrey V. Kravtsov,
Anatoly A. Klypin,
Cristiano Porciani,
Joel R. Primack
Abstract:
[Abridged] We study the angular-momentum profiles of a statistical sample of halos drawn from a high-resolution N-body simulation of the LCDM cosmology. We find that the cumulative mass distribution of specific angular momentum, j, in a halo of mass Mv is well fit by a universal function, M(<j) = Mv μj/(j_0+j). This profile is defined by one shape parameter (μor j_0) in addition to the global sp…
▽ More
[Abridged] We study the angular-momentum profiles of a statistical sample of halos drawn from a high-resolution N-body simulation of the LCDM cosmology. We find that the cumulative mass distribution of specific angular momentum, j, in a halo of mass Mv is well fit by a universal function, M(<j) = Mv μj/(j_0+j). This profile is defined by one shape parameter (μor j_0) in addition to the global spin parameter λ. It follows a power-law over most of the mass, and flattens at large j, with the flattening more pronounced for small values of μ. Compared to a uniform sphere in solid-body rotation, most halos have a higher fraction of their mass in the low- and high-j tails of the distribution. The spatial distribution of angular momentum in halos tends to be cylindrical and is well-aligned within each halo for ~80% of the halos. We investigate two ideas for the origin of this profile. The first is based on a revised version of linear tidal-torque theory combined with extended Press-Schechter mass accretion, and the second focuses on j transport in minor mergers. Finally, we briefly explore implications of the M(<j) profile on the formation of galactic disks assuming that j is conserved during an adiabatic baryonic infall. The implied gas density profile deviates from an exponential disk, with a higher density at small radii and a tail extending to large radii. The steep central density profiles may imply disk scale lengths that are smaller than observed. This is reminiscent of the "angular-momentum problem" seen in hydrodynamic simulations, even though we have assumed perfect j conservation. A possible solution is to associate the central excesses with bulge components and the outer regions with extended gaseous disks.
△ Less
Submitted 2 November, 2000; v1 submitted 31 October, 2000;
originally announced November 2000.
-
Origin and Destiny of Dark Matter Halos: Cosmological Matter Exchange and Metal Enrichment
Authors:
Tsafrir S. Kolatt,
James S. Bullock,
Avishai Dekel,
Joel R. Primack,
Yair Sigad,
Andrey V. Kravtsov,
Anatoly A. Klypin
Abstract:
[Abridged] We analyze the exchange of dark matter between halos, subhalos, and their environments in a high-resolution cosmological N-body simulation of a Lambda CDM cosmology. At each analyzed redshift z we divide the dark matter particles into 4 components: (i) isolated galactic halos, (ii) subhalos, (iii) the diffuse medium of group and cluster halos, and (iv) the background outside of virial…
▽ More
[Abridged] We analyze the exchange of dark matter between halos, subhalos, and their environments in a high-resolution cosmological N-body simulation of a Lambda CDM cosmology. At each analyzed redshift z we divide the dark matter particles into 4 components: (i) isolated galactic halos, (ii) subhalos, (iii) the diffuse medium of group and cluster halos, and (iv) the background outside of virialized halos. We follow the time evolution of the mass distribution and flows between these components and provide fitting functions for the exchange rates. We use our derived exchange rates to gauge the importance of metal redistribution in the universe due solely to gravity-induced interactions. The diffuse metallicity in clusters is predicted to be ~40% that in isolated galaxies (~55% of groups) at z=0, and should be lower only slightly by z=1, consistent with observations. The metallicity of the diffuse media in poor groups is expected to be lower by a factor of 5 by z~2, in agreement with the observed metallicity of damped Ly$α$ systems. The metallicity of the background IGM is predicted to be (1-3)x10^{-4} that of z=0 clusters, also consistent with observations. The agreement of predicted and observed trends indicates that gravitational interaction alone may play an important role in metal enrichment of the intra-cluster and intergalactic media.
△ Less
Submitted 11 October, 2000;
originally announced October 2000.
-
Interaction rates of dark-matter halos and subhalos
Authors:
Tsafrir S. Kolatt,
James S. Bullock,
Yair Sigad,
Andrey V. Kravtsov,
Anatoly A. Klypin,
Joel R. Primack,
Avishai Dekel
Abstract:
[Abridged] The interaction rates of dark-matter halos and subhalos, are computed using high-resolution cosmological N-body simulations of the Lambda-CDM model. Although the number fraction of subhalos of mass $>2x10^{11}\hsolmass$ is only ~10%, we find that the interaction rate of such subhalos is relatively high because they reside in high density environments. At low redshift, the subhalo coll…
▽ More
[Abridged] The interaction rates of dark-matter halos and subhalos, are computed using high-resolution cosmological N-body simulations of the Lambda-CDM model. Although the number fraction of subhalos of mass $>2x10^{11}\hsolmass$ is only ~10%, we find that the interaction rate of such subhalos is relatively high because they reside in high density environments. At low redshift, the subhalo collisions dominate the total collision rate, and even at z=3 they are involved in more than 30% of all collisions. About 40% of the "major" mergers are between subhalos. Therefore subhalo interactions must be incorporated in models of structure formation. We find that the collision rate between halos in physical density units, is $\propto (1+z)^δ$, with $δ= 3-4$, in agreement with earlier simulations and most observational data. We test previous analytic estimates of the interaction rates of subhalos based on statistical models, which could be very inaccurate because of the small number of subhalos and the variation of conditions within small host halos. We find that, while such statistical estimates may severely overestimate the rate within hosts $< 10^{13} \hsolmass$, typical of high redshifts, they are valid for larger hosts regardless of the number of subhalos in them. We find the Makino & Hut (1997) estimate of the subhalo merger rate to be valid for hosts $ > 6x10^{11}\hsolmass$ at all redshifts. The collision rate between subhalos and the central object of their host halo is approximated relatively well using the timescale for dynamical friction in circular orbits. This approximation fails in ~40% of the cases, partly because of deviations from the assumption of circular orbits and partly because of the invalidity of the assumption that the subhalo mass is negligible.
△ Less
Submitted 11 October, 2000;
originally announced October 2000.
-
Cosmological Parameters from Velocities, CMB and Supernovae
Authors:
S. L. Bridle,
I. Zehavi,
A. Dekel,
O. Lahav,
M. P. Hobson,
A. N. Lasenby
Abstract:
We compare and combine likelihood functions of the cosmological parameters Omega_m, h and sigma_8, from peculiar velocities, CMB and type Ia supernovae. These three data sets directly probe the mass in the Universe, without the need to relate the galaxy distribution to the underlying mass via a "biasing" relation. We include the recent results from the CMB experiments BOOMERANG and MAXIMA-1. Our…
▽ More
We compare and combine likelihood functions of the cosmological parameters Omega_m, h and sigma_8, from peculiar velocities, CMB and type Ia supernovae. These three data sets directly probe the mass in the Universe, without the need to relate the galaxy distribution to the underlying mass via a "biasing" relation. We include the recent results from the CMB experiments BOOMERANG and MAXIMA-1. Our analysis assumes a flat Lambda CDM cosmology with a scale-invariant adiabatic initial power spectrum and baryonic fraction as inferred from big-bang nucleosynthesis. We find that all three data sets agree well, overlapping significantly at the 2 sigma level. This therefore justifies a joint analysis, in which we find a joint best fit point and 95 per cent confidence limits of Omega_m=0.28 (0.17,0.39), h=0.74 (0.64,0.86), and sigma_8=1.17 (0.98,1.37). In terms of the natural parameter combinations for these data sigma_8 Omega_m^0.6 = 0.54 (0.40,0.73), Omega_m h = 0.21 (0.16,0.27). Also for the best fit point, Q_rms-ps = 19.7 muK and the age of the universe is 13.2 Gyr.
△ Less
Submitted 12 June, 2000;
originally announced June 2000.
-
Strong Evolution in the Luminosity-Velocity Relation at z>1?
Authors:
James S. Bullock,
Avishai Dekel,
Joel R. Primack,
Rachel S. Somerville
Abstract:
We present a method for constraining the evolution of the galaxy luminosity-velocity (LV) relation in hierarchical scenarios of structure formation. The comoving number density of dark-matter halos with circular velocity of 200 km/s is predicted in favored CDM cosmologies to be nearly constant over the redshift range 0<z<5. Any observed evolution in the density of bright galaxies implies in turn…
▽ More
We present a method for constraining the evolution of the galaxy luminosity-velocity (LV) relation in hierarchical scenarios of structure formation. The comoving number density of dark-matter halos with circular velocity of 200 km/s is predicted in favored CDM cosmologies to be nearly constant over the redshift range 0<z<5. Any observed evolution in the density of bright galaxies implies in turn a corresponding evolution in the LV relation. We consider several possible forms of evolution for the zero-point of the LV relation and predict the corresponding evolution in galaxy number density. The Hubble Deep Field suggests a large deficit of bright (M_V < -19) galaxies at 1.4 < z < 2. If taken at face value, this implies a dimming of the LV zero-point by roughly 2 magnitudes. Deep, wide-field, near-IR selected surveys will provide more secure measurements to compare with our predictions.
△ Less
Submitted 16 May, 2000;
originally announced May 2000.
-
Velocity and Mass Functions of Galactic Halos: Evolution and Environmental Dependence
Authors:
Yair Sigad,
Tsafrir S. Kolatt,
James S. Bullock,
Andrey V. Kravtsov,
Anatoly A. Klypin,
Joel R. Primack,
Avishai Dekel
Abstract:
We study the distribution functions of mass and circular velocity for dark matter halos in N-body simulations of the $Λ$CDM cosmology, addressing redshift and environmental dependence. The dynamical range enables us to resolve subhalos and distinguish them from "distinct" halos. The mass function is compared to analytic models, and is used to derive the more observationally relevant circular vel…
▽ More
We study the distribution functions of mass and circular velocity for dark matter halos in N-body simulations of the $Λ$CDM cosmology, addressing redshift and environmental dependence. The dynamical range enables us to resolve subhalos and distinguish them from "distinct" halos. The mass function is compared to analytic models, and is used to derive the more observationally relevant circular velocity function. The distribution functions in the velocity range 100--500 km/s are well fit by a power-law with two parameters, slope and amplitude. We present the parameter dependence on redshift and provide useful fitting formulae. The amplitudes of the mass functions decrease with z, but, contrary to naive expectation, the comoving density of halos of a fixed velocity ~200 km/s actually increases out to z=5. This is because high-z halos are denser, so a fixed velocity corresponds to a smaller mass. The slope of the velocity function at z=0 is as steep as ~ -4, and the mass and velocity functions of distinct halos steepen with increasing z, while the functions of subhalos do not steepen with $z$, and become even flatter at z>2. A simple observable prediction is that the slope of the velocity function of isolated galaxies is steeper than that of galaxies in groups by as much as unity, reflecting the density biasing of high-velocity halos. We confirm that the Press-Schechter approximation typically overestimates the halo mass function by a factor of ~2, while modified approximations provide improved predictions.
△ Less
Submitted 16 May, 2000;
originally announced May 2000.
-
Measuring the Nonlinear Biasing Function from a Galaxy Redshift Survey
Authors:
Yair Sigad,
Enzo Branchini,
Avishai Dekel
Abstract:
We present a simple method for evaluating the nonlinear biasing function of galaxies from a redshift survey. The nonlinear biasing is characterized by the conditional mean of the galaxy density fluctuation given the underlying mass density fluctuation, or by the associated parameters of mean biasing and nonlinearity (following Dekel & Lahav 1999). Using the distribution of galaxies in cosmologic…
▽ More
We present a simple method for evaluating the nonlinear biasing function of galaxies from a redshift survey. The nonlinear biasing is characterized by the conditional mean of the galaxy density fluctuation given the underlying mass density fluctuation, or by the associated parameters of mean biasing and nonlinearity (following Dekel & Lahav 1999). Using the distribution of galaxies in cosmological simulations, at smoothing of a few Mpc, we find that the mean biasing can be recovered to a good accuracy from the cumulative distribution functions (CDFs) of galaxies and mass, despite the biasing scatter. Then, using a suite of simulations of different cosmological models, we demonstrate that the matter CDF is robust compared to the difference between it and the galaxy CDF, and can be approximated for our purpose by a cumulative log-normal distribution of 1+δwith a single parameter σ. Finally, we show how the nonlinear biasing function can be obtained with adequate accuracy directly from the observed galaxy CDF in redshift space. Thus, the biasing function can be obtained from counts in cells once the rms mass fluctuation at the appropriate scale is assumed a priori. The relative biasing function between different galaxy types is measurable in a similar way. The main source of error is sparse sampling, which requires that the mean galaxy separation be smaller than the smoothing scale. Once applied to redshift surveys such as PSCz, 2dF, SDSS, or DEEP, the biasing function can provide valuable constraints on galaxy formation and structure evolution.
△ Less
Submitted 12 April, 2000; v1 submitted 8 February, 2000;
originally announced February 2000.
-
Non-linear Stochastic Galaxy Biasing in Cosmological Simulations
Authors:
R. S. Somerville,
G. Lemson,
Y. Sigad,
A. Dekel,
G. Kauffmann,
S. D. M. White
Abstract:
We study the biasing relation between dark-matter halos or galaxies and the underlying mass distribution, using cosmological $N$-body simulations in which galaxies are modelled via semi-analytic recipes. The nonlinear, stochastic biasing is quantified in terms of the mean biasing function and the scatter about it as a function of time, scale and object properties. The biasing of galaxies and hal…
▽ More
We study the biasing relation between dark-matter halos or galaxies and the underlying mass distribution, using cosmological $N$-body simulations in which galaxies are modelled via semi-analytic recipes. The nonlinear, stochastic biasing is quantified in terms of the mean biasing function and the scatter about it as a function of time, scale and object properties. The biasing of galaxies and halos shows a general similarity and a characteristic shape, with no galaxies in deep voids and a steep slope in moderately underdense regions. At $\sim 8\hmpc$, the nonlinearity is typically $\lsim 10$ percent and the stochasticity is a few tens of percent, corresponding to $\sim 30$ percent variations in the cosmological parameter $β=Ω^{0.6}/b$. Biasing depends weakly on halo mass, galaxy luminosity, and scale. The time evolution is rapid, with the mean biasing larger by a factor of a few at $z\sim 3$ compared to $z=0$, and with a minimum for the nonlinearity and stochasticity at an intermediate redshift. Biasing today is a weak function of the cosmological model, reflecting the weak dependence on the power-spectrum shape, but the time evolution is more cosmology-dependent, relecting the effect of the growth rate. We provide predictions for the relative biasing of galaxies of different type and color, to be compared with upcoming large redshift surveys. Analytic models in which the number of objects is conserved underestimate the evolution of biasing, while models that explicitly account for merging provide a good description of the biasing of halos and its evolution, suggesting that merging is a crucial element in the evolution of biasing.
△ Less
Submitted 3 December, 1999;
originally announced December 1999.
-
Cosmic Flows 99: Conference Summary
Authors:
Avishai Dekel
Abstract:
I address the following issues: All bulk velocity measurements (but one) are consistent with our standard gravitational instability theory. New accurate data and reconstruction methods allow high-resolution dynamical analysis nearby, revealing Virgo, Ursa Major and Fornax as attractors. Large peculiar-velocity surveys enable robust reconstruction of the dynamical fields on the Great-Attractor sc…
▽ More
I address the following issues: All bulk velocity measurements (but one) are consistent with our standard gravitational instability theory. New accurate data and reconstruction methods allow high-resolution dynamical analysis nearby, revealing Virgo, Ursa Major and Fornax as attractors. Large peculiar-velocity surveys enable robust reconstruction of the dynamical fields on the Great-Attractor scale. A decomposition of the velocity field into its local and tidal components indicates the presence of big perturbations further away. Cluster velocities start exploring very large scales, revealing Coma, Shapely and other mass enhancements, and constraining a possible local Hubble bubble. Supernovae type Ia (SNIa) are very promising for cosmic flow analysis. Peculiar velocities do provide unique valuable constraints on cosmological parameters, e.g., 0.3<Omega_m<1 (95% confidence) independent of biasing. Jointly with other data they can confine other parameters such as Omega_Lambda, h, sigma_8, n, and the biasing. Nontrivial features of the biasing scheme can explain much of the span of estimates for beta. Quantitative error analysis is essential in our maturing field; every method ought to be calibrated with suitable mock catalogs, that are offered as benchmarks.
△ Less
Submitted 29 November, 1999;
originally announced November 1999.
-
Cluster versus POTENT Density and Velocity Fields: Cluster Biasing and Omega
Authors:
E. Branchini,
I. Zehavi,
M. Plionis,
A. Dekel
Abstract:
The density and velocity fields as extracted from the Abell/ACO clusters are compared to the corresponding fields recovered by the POTENT method from the Mark~III peculiar velocities of galaxies. In order to minimize non-linear effects and to deal with ill-sampled regions we smooth both fields using a Gaussian window with radii ranging between $12 - 20\hmpc$. The density and velocity fields with…
▽ More
The density and velocity fields as extracted from the Abell/ACO clusters are compared to the corresponding fields recovered by the POTENT method from the Mark~III peculiar velocities of galaxies. In order to minimize non-linear effects and to deal with ill-sampled regions we smooth both fields using a Gaussian window with radii ranging between $12 - 20\hmpc$. The density and velocity fields within $70\hmpc$ exhibit similarities, qualitatively consistent with gravitational instability theory and a linear biasing relation between clusters and mass. The random and systematic errors are evaluated with the help of mock catalogs. Quantitative comparisons within a volume containing $\sim\!12$ independent samples yield $\betac\equivΩ^{0.6}/b_c=0.22\pm0.08$, where $b_c$ is the cluster biasing parameter at $15\hmpc$. If $b_c \sim 4.5$, as indicated by the cluster correlation function, our result is consistent with $Ω\sim 1$.
△ Less
Submitted 3 November, 1999;
originally announced November 1999.
-
Clustering of high-redshift galaxies: relating LBGs to dark matter halos
Authors:
Risa H. Wechsler,
James S. Bullock,
Rachel S. Somerville,
Tsafrir Kolatt,
Joel R. Primack,
George R. Blumenthal,
Avishai Dekel
Abstract:
We investigate the clustering properties of high-redshift galaxies within three competing scenarios for assigning luminous galaxies to dark matter halos from N-body simulations: a one galaxy per massive halo model, a quiescent star formation model, and a collisional starburst model. We compare these models to observations of Lyman-Break galaxies at z~3$ With current data and the simple statistic…
▽ More
We investigate the clustering properties of high-redshift galaxies within three competing scenarios for assigning luminous galaxies to dark matter halos from N-body simulations: a one galaxy per massive halo model, a quiescent star formation model, and a collisional starburst model. We compare these models to observations of Lyman-Break galaxies at z~3$ With current data and the simple statistic used here, one cannot rule out any of these models, but we see potential for finding distinguishing features using statistics that are sensitive to the tails of the distribution, and statistics based on the number of multiple galaxies per halo, which we explore in an ongoing study.
△ Less
Submitted 20 October, 1999;
originally announced October 1999.
-
Galaxy Interactions and Starbursts at High Redshift
Authors:
R. S. Somerville,
G. Rosenfeld,
T. S. Kolatt,
A. Dekel,
J. C. Mihos,
J. R. Primack
Abstract:
Using high resolution N-body simulations with hydrodynamics and star formation, we investigate interactions and the resulting starbursts in galaxies with properties typical of $z\sim 3$. We apply spectral population models to produce mock-HST images, and discuss the observed magnitude, color, and morphological appearance of our simulated galaxies in both the rest-UV and rest-visual bands.
Using high resolution N-body simulations with hydrodynamics and star formation, we investigate interactions and the resulting starbursts in galaxies with properties typical of $z\sim 3$. We apply spectral population models to produce mock-HST images, and discuss the observed magnitude, color, and morphological appearance of our simulated galaxies in both the rest-UV and rest-visual bands.
△ Less
Submitted 19 October, 1999;
originally announced October 1999.
-
Cosmological Parameters and Power Spectrum from Peculiar Velocities
Authors:
Idit Zehavi,
Avishai Dekel
Abstract:
The power spectrum of mass density fluctuations is evaluated from the Mark III and the SFI catalogs of peculiar velocities by a maximum likelihood analysis, using parametric models for the power spectrum and for the errors. The applications to the two different data sets, using generalized CDM models with and without COBE normalization, give consistent results. The general result is a relatively…
▽ More
The power spectrum of mass density fluctuations is evaluated from the Mark III and the SFI catalogs of peculiar velocities by a maximum likelihood analysis, using parametric models for the power spectrum and for the errors. The applications to the two different data sets, using generalized CDM models with and without COBE normalization, give consistent results. The general result is a relatively high amplitude of the power spectrum, with σ_8 Ω_m^{0.6} = 0.8+/-0.2 at 90% confidence. Casting the results in the Ω_m - Ω_Λplane, yields complementary constraints to those of the high-redshift supernovae, together favoring a nearly flat, unbound and accelerating universe with comparable contributions from Ω_m and Ω_Λ. Further implications on the cosmological parameters, arising from a joint analysis of the velocities together with small-scale CMB anisotropies and the high-redshift supernovae, are also briefly described.
△ Less
Submitted 29 September, 1999;
originally announced September 1999.
-
Profiles of dark haloes: evolution, scatter, and environment
Authors:
James S. Bullock,
Tsafrir S. Kolatt,
Yair Sigad,
Rachel S. Somerville,
Andrey V. Kravtsov,
Anatoly A. Klypin,
Joel R. Primack,
Avishai Dekel
Abstract:
We study dark-matter halo density profiles in a high-resolution N-body simulation of an LCDM cosmology. Our statistical sample contains ~5000 haloes in the range 10^{11}-10^{14} M_sun. The profiles are parameterized by an NFW form with two parameters, an inner radius r_s and a virial radius r_v, and we define the halo concentration c_v = r_v/r_s. First, we find that, for a given halo mass, the r…
▽ More
We study dark-matter halo density profiles in a high-resolution N-body simulation of an LCDM cosmology. Our statistical sample contains ~5000 haloes in the range 10^{11}-10^{14} M_sun. The profiles are parameterized by an NFW form with two parameters, an inner radius r_s and a virial radius r_v, and we define the halo concentration c_v = r_v/r_s. First, we find that, for a given halo mass, the redshift dependence of the median concentration is c_v ~ 1/(1+z), corresponding to a roughly constant r_s with redshift. We present an improved analytic treatment of halo formation that fits the measured relations between halo parameters and their redshift dependence. The implications are that high-redshift galaxies are predicted to be more extended and dimmer than expected before. Second, we find that the scatter in log(c_v) is ~0.18, corresponding to a scatter in maximum rotation velocities of dV/V ~ 0.12. We discuss implications for modelling the Tully-Fisher relation, which has a smaller reported intrinsic scatter. Third, haloes in dense environments tend to be more concentrated than isolated haloes. These results suggest that c_v is an essential parameter for the theory of galaxy modelling, and we briefly discuss implications for the universality of the Tully-Fisher relation, the formation of low surface brightness galaxies, and the origin of the Hubble sequence.
△ Less
Submitted 10 October, 2000; v1 submitted 15 August, 1999;
originally announced August 1999.
-
Young galaxies: what turns them on?
Authors:
Tsafrir S. Kolatt,
James S. Bullock,
Rachel S. Somerville,
Yair Sigad,
Patrik Jonsson,
Andrey V. Kravtsov,
Anatoly A. Klypin,
Joel R. Primack,
Sandra M. Faber,
Avishai Dekel,
.
Abstract:
Lyman break galaxies (LBGs) at $z \sim 3$ exhibit number densities and clustering similar to local $L_*$ galaxies; however, their star formation rates (SFRs) are much higher. We explore the scenario in which LBGs are starburst galaxies triggered by collisions, and thus provide an explanation for these key properties. The relative paucity of starburst galaxies at low redshift can be attributed to…
▽ More
Lyman break galaxies (LBGs) at $z \sim 3$ exhibit number densities and clustering similar to local $L_*$ galaxies; however, their star formation rates (SFRs) are much higher. We explore the scenario in which LBGs are starburst galaxies triggered by collisions, and thus provide an explanation for these key properties. The relative paucity of starburst galaxies at low redshift can be attributed to a much higher collision rate in the past. We use high-resolution cosmological $N$-body simulations and a hierarchical halo finder to estimate the galaxy collision rate as a function of time in the popular $Λ$CDM cosmological model. We find that bright collisional starbursts are frequent enough to account for most of the high-$z$ ($2.5-4.5$) LBGs. Moreover, many of the objects are of relatively small mass, but they cluster about large-mass halos. They therefore exhibit strong clustering, similar to that observed and stronger than that of the relevant massive halos.
△ Less
Submitted 7 June, 1999;
originally announced June 1999.
-
Evidence for a Positive Cosmological Constant from Flows of Galaxies and Distant Supernovae
Authors:
Idit Zehavi,
Avishai Dekel
Abstract:
Recent observations of high-redshift supernovae seem to suggest that the global geometry of the Universe may be affected by a `cosmological constant', which acts to accelerate the expansion rate with time. But these data by themselves still permit an open universe of low mass density and no cosmological constant. Here we derive an independent constraint on the lower bound to the mass density, ba…
▽ More
Recent observations of high-redshift supernovae seem to suggest that the global geometry of the Universe may be affected by a `cosmological constant', which acts to accelerate the expansion rate with time. But these data by themselves still permit an open universe of low mass density and no cosmological constant. Here we derive an independent constraint on the lower bound to the mass density, based on deviations of galaxy velocities from a smooth universal expansion. This constraint rules out a low-density open universe with a vanishing cosmological constant, and together the two favour a nearly flat universe in which the contributions from mass density and the cosmological constant are comparable. This type of universe, however, seems to require a degree of fine tuning of the initial conditions that is in apparent conflict with `common wisdom'.
△ Less
Submitted 29 September, 1999; v1 submitted 16 April, 1999;
originally announced April 1999.
-
Large-Scale Power Spectrum and Cosmological Parameters from SFI Peculiar Velocities
Authors:
Wolfram Freudling,
Idit Zehavi,
Luiz N. da Costa,
Avishai Dekel,
Amiram Eldar,
Riccardo Giovanelli,
Martha P. Haynes,
John J. Salzer,
Gary Wegner,
Saleem Zaroubi
Abstract:
We estimate the power spectrum of mass density fluctuations from peculiar velocities of galaxies by applying an improved maximum-likelihood technique to the new all-sky SFI catalog. Parametric models are used for the power spectrum and the errors, and the free parameters are determined by assuming Gaussian velocity fields and errors and maximizing the probability of the data given the model. It…
▽ More
We estimate the power spectrum of mass density fluctuations from peculiar velocities of galaxies by applying an improved maximum-likelihood technique to the new all-sky SFI catalog. Parametric models are used for the power spectrum and the errors, and the free parameters are determined by assuming Gaussian velocity fields and errors and maximizing the probability of the data given the model. It has been applied to generalized CDM models with and without COBE normalization. The method has been carefully tested using artificial SFI catalogs. The most likely distance errors are found to be similar to the original error estimates in the SFI data. The general result that is not very sensitive to the prior model used is a relatively high amplitude of the power spectrum. For example, at k=0.1 h/Mpc we find P(k)Ω^{1.2}=(4.4+/-1.7)X10^3 (Mpc/h)^3. An integral over the power spectrum yields σ_8Ω^{0.6}=0.82+/-0.12. Model-dependent constraints on the cosmological parameters are obtained for families of CDM models. For example, for COBE-normalized ΛCDM models (scalar fluctuations only), the maximum-likelihood result can be approximated by Ωn^{2} h_{60}^{1.3} =0.58+/-0.11. The formal random errors quoted correspond to the 90% confidence level. The total uncertainty, including systematic errors associated with nonlinear effects, may be larger by a factor of ~2. These results are in agreement with an application of a similar method to other data (Mark III).
△ Less
Submitted 9 April, 1999;
originally announced April 1999.
-
Cosmological Implications of Lyman-Break Galaxy Clustering
Authors:
Risa H. Wechsler,
Michael A. K. Gross,
Joel R. Primack,
George R. Blumenthal,
Avishai Dekel,
Rachel S. Somerville
Abstract:
We review our analysis of the clustering properties of ``Lyman-break'' galaxies (LBGs) at redshift z~3, previously discussed in Wechsler et al (1998). We examine the likelihood of spikes found by Steidel et al (1998) in the redshift distribution of LBGs, within a suite of models for the evolution of structure in the Universe. Using high-resolution dissipationless N-body simulations, we analyze d…
▽ More
We review our analysis of the clustering properties of ``Lyman-break'' galaxies (LBGs) at redshift z~3, previously discussed in Wechsler et al (1998). We examine the likelihood of spikes found by Steidel et al (1998) in the redshift distribution of LBGs, within a suite of models for the evolution of structure in the Universe. Using high-resolution dissipationless N-body simulations, we analyze deep pencil-beam surveys from these models in the same way that they are actually observed, identifying LBGs with the most massive dark matter halos. We find that all the models (with SCDM as a marginal exception) have a substantial probability of producing spikes similar to those observed, because the massive halos are much more clumped than the underlying matter -- i.e., they are biased. Therefore, the likelihood of such a spike is not a good discriminator among these models. The LBG correlation functions are less steep than galaxies today (gamma~1.4), but show similar or slightly longer correlation lengths. We have extened this analysis and include a preliminary comparison to the new data presented in Adelberger et al (1998). We also discuss work in progress, in which we use semi-analytic models to identify Lyman-break galaxies within dark-matter halos.
△ Less
Submitted 20 December, 1998;
originally announced December 1998.
-
POTENT Reconstruction from Mark III Velocities
Authors:
A. Dekel,
A. Eldar,
T. Kolatt,
A. Yahil,
J. A. Willick,
S. M. Faber,
S. Courteau,
D. Burstein
Abstract:
We present an improved POTENT method for reconstructing the velocity and mass density fields from radial peculiar velocities, test it with mock catalogs, and apply it to the Mark III Catalog. Method improvments: (a) inhomogeneous Malmquist bias is reduced by grouping and corrected in forward or inverse analyses of inferred distances, (b) the smoothing into a radial velocity field is optimized to…
▽ More
We present an improved POTENT method for reconstructing the velocity and mass density fields from radial peculiar velocities, test it with mock catalogs, and apply it to the Mark III Catalog. Method improvments: (a) inhomogeneous Malmquist bias is reduced by grouping and corrected in forward or inverse analyses of inferred distances, (b) the smoothing into a radial velocity field is optimized to reduce window and sampling biases, (c) the density is derived from the velocity using an improved nonlinear approximation, and (d) the computational errors are made negligible. The method is tested and optimized using mock catalogs based on an N-body simulation that mimics our cosmological neighborhood, and the remaining errors are evaluated quantitatively. The Mark III catalog, with ~3300 grouped galaxies, allows a reliable reconstruction with fixed Gaussian smoothing of 10-12 Mpc/h out to ~60 Mpc/h. We present maps of the 3D velocity and mass-density fields and the corresponding errors. The typical systematic and random errors in the density fluctuations inside 40 Mpc/h are \pm 0.13 and \pm 0.18. The recovered mass distribution resembles in its gross features the galaxy distribution in redshift surveys and the mass distribution in a similar POTENT analysis of a complementary velocity catalog (SFI), including the Great Attractor, Perseus-Pisces, and the void in between. The reconstruction inside ~40 Mpc/h is not affected much by a revised calibration of the distance indicators (VM2, tailored to match the velocities from the IRAS 1.2Jy redshift survey). The bulk velocity within the sphere of radius 50 Mpc/h about the Local Group is V_50=370 \pm 110 km/s (including systematic errors), and is shown to be mostly generated by external mass fluctuations. With the VM2 calibration, V_50 is reduced to 305 \pm 110 km/s.
△ Less
Submitted 9 December, 1998;
originally announced December 1998.
-
Wiener Reconstruction of Large-Scale Structure from Peculiar Velocities
Authors:
Saleem Zaroubi,
Yehuda Hoffman,
Avishai Dekel
Abstract:
We present an alternative, Bayesian method for large-scale reconstruction from observed peculiar velocity data. The method stresses a rigorous treatment of the random errors and it allows extrapolation into poorly sampled regions in real space or in k-space. A likelihood analysis is used to determine the fluctuation power spectrum, followed by a Wiener Filter (WF) analysis to obtain the minimum-…
▽ More
We present an alternative, Bayesian method for large-scale reconstruction from observed peculiar velocity data. The method stresses a rigorous treatment of the random errors and it allows extrapolation into poorly sampled regions in real space or in k-space. A likelihood analysis is used to determine the fluctuation power spectrum, followed by a Wiener Filter (WF) analysis to obtain the minimum-variance mean fields of velocity and mass density. Constrained Realizations (CR) are then used to sample the statistical scatter about the WF mean field. The WF/CR method is applied as a demonstration to the Mark III data with 1200 km/s, 900 km/s, and 500 km/s resolutions. The main reconstructed structures are consistent with those extracted by the POTENT method. A comparison with the structures in the distribution of IRAS 1.2Jy galaxies yields a general agreement. The reconstructed velocity field is decomposed into its divergent and tidal components relative to a cube of +/-8000 km/s centered on the Local Group. The divergent component is very similar to the velocity field predicted from the distribution of IRAS galaxies. The tidal component is dominated by a bulk flow of 194 +/- 32 km/s towards the general direction of the Shapley concentration, and it also indicates a significant quadrupole.
△ Less
Submitted 19 October, 1998;
originally announced October 1998.
-
Galaxy Biasing: Nonlinear, Stochastic and Measurable
Authors:
Avishai Dekel
Abstract:
I describe a general formalism for galaxy biasing (Dekel & Lahav 1998) and its application to measurements of beta (=Omega^0.6/b), e.g. via direct comparisons of light and mass and via redshift distortions. The linear and deterministic relation g=b*d between the density fluctuation fields of galaxies g and mass d is replaced by the conditional distribution P(g|d) of these as random fields, smoot…
▽ More
I describe a general formalism for galaxy biasing (Dekel & Lahav 1998) and its application to measurements of beta (=Omega^0.6/b), e.g. via direct comparisons of light and mass and via redshift distortions. The linear and deterministic relation g=b*d between the density fluctuation fields of galaxies g and mass d is replaced by the conditional distribution P(g|d) of these as random fields, smoothed at a given scale and at a given time. The mean biasing and its nonlinearity are characterized by the conditional mean <g|d>=b(d)*d and the local scatter by the conditional variance s_b^2(d). This scatter arises from hidden effects on galaxy formation and from shot noise. For applications involving second-order local moments, the biasing is defined by three natural parameters: the slope b_h of the regression of g on d (replacing b), a nonlinearity parameter b_t, and a scatter parameter s_b. The ratio of variances b_v^2 and the correlation coefficient r mix these parameters. The nonlinearity and scatter lead to underestimates of order b_t^2/b_h^2 and s_b^2/b_h^2 in the different estimators of beta, which may partly explain the range of estimates. Local stochasticity affects the redshift-distortion analysis only by limiting the useful range of scales. In this range, for linear stochastic biasing, the analysis reduces to Kaiser's formula for b_h (not b_v) independent of the scatter. The distortion analysis is affected by nonlinearity but in a weak way. Estimates of the nontrivial features of the biasing scheme are made based on simulations and toy models, and a new method for measuring them via distribution functions is proposed.
△ Less
Submitted 23 September, 1998;
originally announced September 1998.
-
Evaluating Semi-Analytic Halo Merging Histories
Authors:
Rachel S. Somerville,
Gerard Lemson,
Tsafrir S. Kolatt,
Avishai Dekel
Abstract:
We evaluate the accuracy of semi-analytic merger-trees by comparing them with the merging histories of dark-matter halos in N-body simulations, focusing on the joint distribution of the number of progenitors and their masses. We first confirm that the halo mass function as predicted directly by the Press-Schechter (PS) model deviates from the simulations by up to 50% depending on the mass scale…
▽ More
We evaluate the accuracy of semi-analytic merger-trees by comparing them with the merging histories of dark-matter halos in N-body simulations, focusing on the joint distribution of the number of progenitors and their masses. We first confirm that the halo mass function as predicted directly by the Press-Schechter (PS) model deviates from the simulations by up to 50% depending on the mass scale and redshift, while the means of the projected distributions of progenitor number and mass for a halo of a given mass are more accurately predicted by the Extended PS model. We then use the full merger trees to study the joint distribution as a function of redshift and parent-halo mass. We find that while the deviation of the mean quantities due to the inaccuracy of the Extended PS model partly propagates into the higher moments of the distribution, the merger-tree procedure does not introduce a significant additional source of error. In particular, certain properties of the merging history such as the mass ratio of the progenitors and the total accretion rate are reproduced quite accurately for galaxy sized halos ($\sim 10^{12}\msun$), and less so for larger masses. We conclude that although there could be $\sim 50%$ deviations in the absolute numbers and masses of progenitors and in the higher order moment of these distributions, the relative properties of progenitors for a given halo are reproduced fairly well by the merger trees. They can thus provide a useful framework for modelling galaxy formation once the above-mentioned limitations are taken into account.
△ Less
Submitted 27 July, 1998;
originally announced July 1998.
-
Stochastic Nonlinear Galaxy Biasing
Authors:
Avishai Dekel,
Ofer Lahav
Abstract:
We propose a general formalism for galaxy biasing and apply it to methods for measuring cosmological parameters, such as regression of light versus mass, the analysis of redshift distortions, measures involving skewness and the cosmic virial theorem. The common linear and deterministic relation g=b*d between the density fluctuation fields of galaxies g and mass d is replaced by the conditional d…
▽ More
We propose a general formalism for galaxy biasing and apply it to methods for measuring cosmological parameters, such as regression of light versus mass, the analysis of redshift distortions, measures involving skewness and the cosmic virial theorem. The common linear and deterministic relation g=b*d between the density fluctuation fields of galaxies g and mass d is replaced by the conditional distribution P(g|d) of these random fields, smoothed at a given scale at a given time. The nonlinearity is characterized by the conditional mean <g|d>=b(d)*d, while the local scatter is represented by the conditional variance s_b^2(d) and higher moments. The scatter arises from hidden factors affecting galaxy formation and from shot noise unless it has been properly removed. For applications involving second-order local moments, the biasing is defined by three natural parameters: the slope b_h of the regression of g on d, a nonlinearity b_t, and a scatter s_b. The ratio of variances b_v^2 and the correlation coefficient r mix these parameters. The nonlinearity and the scatter lead to underestimates of order b_t^2/b_h^2 and s_b^2/b_h^2 in the different estimators of beta (=Omega^0.6/b_h). The nonlinear effects are typically smaller. Local stochasticity affects the redshift-distortion analysis only by limiting the useful range of scales, especially for power spectra. In this range, for linear stochastic biasing, the analysis reduces to Kaiser's formula for b_h (not b_v), independent of the scatter. The distortion analysis is affected by nonlinear properties of biasing but in a weak way. Estimates of the nontrivial features of the biasing scheme are made based on simulations and toy models, and strategies for measuring them are discussed. They may partly explain the range of estimates for beta.
△ Less
Submitted 30 November, 1998; v1 submitted 14 June, 1998;
originally announced June 1998.
-
A Local Hubble Bubble from SNe Ia?
Authors:
Idit Zehavi,
Adam G. Riess,
Robert P. Kirshner,
Avishai Dekel
Abstract:
We analyze the monopole in the peculiar velocities of 44 Type Ia supernovae (SNe Ia) to test for a local void. The sample extends from 20 to 300 Mpc/h, with distances, deduced from light-curve shapes, accurate to ~6%. Assuming Omega_m=1 and Omega_lambda=0, the most significant deviation we find from the Hubble law is an outwards flow of (6.6+/-2.2)% inside a sphere of radius 70 Mpc/h as would be…
▽ More
We analyze the monopole in the peculiar velocities of 44 Type Ia supernovae (SNe Ia) to test for a local void. The sample extends from 20 to 300 Mpc/h, with distances, deduced from light-curve shapes, accurate to ~6%. Assuming Omega_m=1 and Omega_lambda=0, the most significant deviation we find from the Hubble law is an outwards flow of (6.6+/-2.2)% inside a sphere of radius 70 Mpc/h as would be produced by a void of ~20% underdensity surrounded by a dense shell. This shell roughly coincides with the local Great Walls. Monte Carlo analyses, using Gaussian errors or bootstrap resampling, show the probability for chance occurrence of this result out of a pure Hubble flow to be ~2%. The monopole could be contaminated by higher moments of the velocity field, especially a quadrupole, which are not properly probed by the current limited sky coverage. The void would be less significant if Omega_m is low and Omega_lambda is high. It would be more significant if one outlier is removed from the sample, or if the size of the void is constrained a-priori. This putative void is not in significant conflict with any of the standard cosmological scenarios. It suggests that the Hubble constant as determined within 70 Mpc/h could be overestimated by ~6% and the local value of Omega may be underestimated by ~20%. While the present evidence for a local void is marginal in this data set, the analysis shows that the accumulation of SNe Ia distances will soon provide useful constraints on elusive and important aspects of regional cosmic dynamics.
△ Less
Submitted 7 April, 1998; v1 submitted 19 February, 1998;
originally announced February 1998.
-
Implications of Spikes in the Redshift Distribution of $z\sim3$ Galaxies
Authors:
Risa H. Wechsler,
Michael A. K. Gross,
Joel R. Primack,
George R. Blumenthal,
Avishai Dekel
Abstract:
We address the high peaks found by Steidel et al (1997) in the redshift distribution of ``Lyman-break'' objects (LBOs) at redshift z~3. The highest spike represents a relative overdensity of 2.6 in the distribution of LBOs in pixels of comoving size ~10Mpc/h. We examine the likelihood of such a spike in the redshift distribution within a suite of models for the evolution of structure in the Univ…
▽ More
We address the high peaks found by Steidel et al (1997) in the redshift distribution of ``Lyman-break'' objects (LBOs) at redshift z~3. The highest spike represents a relative overdensity of 2.6 in the distribution of LBOs in pixels of comoving size ~10Mpc/h. We examine the likelihood of such a spike in the redshift distribution within a suite of models for the evolution of structure in the Universe, including models with Omega=1 (SCDM and CHDM) and with Omega=0.4-0.5 (LCDM and OCDM). Using high-resolution dissipationless N-body simulations, we analyze deep pencil-beam surveys from these models in the same way that they are actually observed, identifying LBOs with the most massive dark matter halos. We find that all the models (with SCDM as a marginal exception) have a substantial probability of producing spikes similar to those observed, because the massive halos are much more clumped than the underlying matter -- i.e., they are biased. Therefore, the likelihood of such a spike is not a good discriminator among these models. We find in these models that the mean biasing parameter b of LBOs with respect to dark matter varies within a range b ~2-5 on a scale of ~10Mpc/h. We also compute the two-body correlation functions of LBOs predicted in these models. The LBO correlation functions are less steep than galaxies today (gamma ~1.4), but show similar or slightly longer correlation lengths.
△ Less
Submitted 30 April, 1998; v1 submitted 10 December, 1997;
originally announced December 1997.
-
IRAS versus POTENT Density Fields on Large Scales: Biasing and Omega
Authors:
Y. Sigad,
A. Eldar,
A. Dekel,
M. A. Strauss,
A. Yahil
Abstract:
The galaxy density field as extracted from the IRAS 1.2 Jy redshift survey is compared to the mass density field as reconstructed by the POTENT method from the Mark III catalog of peculiar velocities. The reconstruction is done with Gaussian smoothing of radius 12 h^{-1}Mpc, and the comparison is carried out within volumes of effective radii 31-46 h^{-1}Mpc, containing approximately 10-26 indepe…
▽ More
The galaxy density field as extracted from the IRAS 1.2 Jy redshift survey is compared to the mass density field as reconstructed by the POTENT method from the Mark III catalog of peculiar velocities. The reconstruction is done with Gaussian smoothing of radius 12 h^{-1}Mpc, and the comparison is carried out within volumes of effective radii 31-46 h^{-1}Mpc, containing approximately 10-26 independent samples. Random and systematic errors are estimated from multiple realizations of mock catalogs drawn from a simulation that mimics the observed density field in the local universe. The relationship between the two density fields is found to be consistent with gravitational instability theory in the mildly nonlinear regime and a linear biasing relation between galaxies and mass. We measure beta = Omega^{0.6}/b_I = 0.89 \pm 0.12 within a volume of effective radius 40 h^{-1}Mpc, where b_I is the IRAS galaxy biasing parameter at 12 h^{-1}Mpc. This result is only weakly dependent on the comparison volume, suggesting that cosmic scatter is no greater than \pm 0.1. These data are thus consistent with Omega=1 and b_I\approx 1. If b_I>0.75, as theoretical models of biasing indicate, then Omega>0.33 at 95% confidence. A comparison with other estimates of beta suggests scale-dependence in the biasing relation for IRAS galaxies.
△ Less
Submitted 13 August, 1997;
originally announced August 1997.
-
Cosmological Implications of Large-Scale Flows
Authors:
Avishai Dekel
Abstract:
Cosmological implications of the observed large-scale peculiar velocities are reviewed, alone or combined with redshift surveys and CMB data. The latest version of the POTENT method for reconstructing the underlying three-dimensional velocity and mass-density fields is described. The initial fluctuations and the nature of the dark matter are addressed via statistics such as bulk flow and mass po…
▽ More
Cosmological implications of the observed large-scale peculiar velocities are reviewed, alone or combined with redshift surveys and CMB data. The latest version of the POTENT method for reconstructing the underlying three-dimensional velocity and mass-density fields is described. The initial fluctuations and the nature of the dark matter are addressed via statistics such as bulk flow and mass power spectrum. The focus is on constraining the mass density parameter $Ω$, directly or via the parameter $β$ which involves the unknown relation between galaxies and mass. The acceptable range for $Ω$ is found to be $0.4-1.0$. The range of $β$ estimates is likely to reflect non-trivial features in the galaxy biasing scheme, such as scale dependence. Similar constraints on $Ω$ and $Λ$ from global measures are summarized.
△ Less
Submitted 6 May, 1997;
originally announced May 1997.
-
Maximum-Likelihood Comparisons of Tully-Fisher and Redshift Data: Constraints on Omega and Biasing
Authors:
Jeffrey A. Willick,
Michael A. Strauss,
Avishai Dekel,
Tsafrir Kolatt
Abstract:
We compare Tully-Fisher (TF) data for 838 galaxies within cz=3000 km/sec from the Mark III catalog to the peculiar velocity and density fields predicted from the 1.2 Jy IRAS redshift survey. Our goal is to test the relation between the galaxy density and velocity fields predicted by gravitational instability theory and linear biasing, and thereby to estimate $β_I = Ω^{0.6}/b_I,$ where $b_I$ is t…
▽ More
We compare Tully-Fisher (TF) data for 838 galaxies within cz=3000 km/sec from the Mark III catalog to the peculiar velocity and density fields predicted from the 1.2 Jy IRAS redshift survey. Our goal is to test the relation between the galaxy density and velocity fields predicted by gravitational instability theory and linear biasing, and thereby to estimate $β_I = Ω^{0.6}/b_I,$ where $b_I$ is the linear bias parameter for IRAS galaxies. Adopting the IRAS velocity and density fields as a prior model, we maximize the likelihood of the raw TF observables, taking into account the full range of selection effects and properly treating triple-valued zones in the redshift-distance relation. Extensive tests with realistic simulated galaxy catalogs demonstrate that the method produces unbiased estimates of $β_I$ and its error. When we apply the method to the real data, we model the presence of a small but significant velocity quadrupole residual (~3.3% of Hubble flow), which we argue is due to density fluctuations incompletely sampled by IRAS. The method then yields a maximum likelihood estimate $β_I=0.49\pm 0.07$ (1-sigma error). We discuss the constraints on $Ω$ and biasing that follow if we assume a COBE-normalized CDM power spectrum. Our model also yields the 1-D noise noise in the velocity field, including IRAS prediction errors, which we find to be be 125 +/- 20 km/sec.
△ Less
Submitted 31 December, 1996; v1 submitted 30 December, 1996;
originally announced December 1996.
-
Measuring Omega
Authors:
Avishai Dekel,
David Burstein,
Simon D. M. White
Abstract:
We were asked to debate the value of the cosmological mass-density parameter Omega. Is Omega_m=1 in accordance with the simplest model? Is Omega_m much smaller as indicated by some observations? There is conflicting evidence. We lay out the various methods for measuring Omega_m, mention new developments and current estimates, and focus on prospects versus the associated difficulties. We try to s…
▽ More
We were asked to debate the value of the cosmological mass-density parameter Omega. Is Omega_m=1 in accordance with the simplest model? Is Omega_m much smaller as indicated by some observations? There is conflicting evidence. We lay out the various methods for measuring Omega_m, mention new developments and current estimates, and focus on prospects versus the associated difficulties. We try to shed light on the uncertainties that are responsible for the span of estimates for Omega_m.
We divide the methods into the following four classes:
a. Global measures of the properties of space-time that constrain combinations of Omega_m and the other cosmological parameters (Lambda, Ho, to).
b. Virialized Systems. Nonlinear dynamics within galaxies and clusters on comoving scales 1-10 Mpc/h.
c. Large-scale structure. Mildly-nonlinear gravitational dynamics of fluctuations on scales 10-100 Mpc/h, in particular cosmic flows.
d. Growth rate of fluctuations from the last scattering of the CMB or from high redshift galaxies to the present.
Methods (b) typically yield low values of Omega_m=0.2-0.3. Methods (a) and (c) typically indicate higher values of Omega_m=0.4-1.
△ Less
Submitted 13 November, 1996;
originally announced November 1996.
-
Large Scale Power Spectrum from Peculiar Velocities Via Likelihood Analysis
Authors:
S. Zaroubi,
I. Zehavi,
A. Dekel,
Y. Hoffman,
T. Kolatt
Abstract:
The power spectrum (PS) of mass density fluctuations, independent of `biasing', is estimated from the Mark III catalog of peculiar velocities using Bayesian statistics. A parametric model is assumed for the PS, and the free parameters are determined by maximizing the probability of the model given the data. The method has been tested using detailed mock catalogs. It has been applied to generaliz…
▽ More
The power spectrum (PS) of mass density fluctuations, independent of `biasing', is estimated from the Mark III catalog of peculiar velocities using Bayesian statistics. A parametric model is assumed for the PS, and the free parameters are determined by maximizing the probability of the model given the data. The method has been tested using detailed mock catalogs. It has been applied to generalized CDM models with and without COBE normalization.
The robust result for all the models is a relatively high PS, with $P(k) Ω^{1.2} = (4.8 \pm 1.5) \times 10^3 (Mpc/h)^3$ at $k=0.1 h/Mpc$. An extrapolation to smaller scales using the different CDM models yields $σ_8 Ω^{0.6} = 0.88 \pm 0.15$. The peak is weakly constrained to the range $0.02 \leq k \leq 0.06 h/Mpc$. These results are consistent with a direct computation of the PS (Kolatt & Dekel 1996). When compared to galaxy-density surveys, the implied values for $β$ ($\equiv Ω^{0.6}/b$) are of order unity to within 25%.
The parameters of the COBE-normalized, flat CDM model are confined by a 90% likelihood contour of the sort $Ωh_{50}^μn^ν= 0.8 \pm 0.2$, where $μ= 1.3$ and $ν= 3.4, 2.0$ for models with and without tensor fluctuations respectively. For open CDM the powers are $μ= 0.95$ and $ν= 1.4$ (no tensor fluctuations). A $Γ$-shape model free of COBE normalization yields only a weak constraint: $Γ= 0.4 \pm 0.2$.
△ Less
Submitted 3 February, 1997; v1 submitted 28 October, 1996;
originally announced October 1996.
-
Homogeneous Velocity-Distance Data for Peculiar Velocity Analysis. III. The Mark III Catalog of Galaxy Peculiar Velocities
Authors:
Jeffrey A. Willick,
Stephane Courteau,
S. M. Faber,
David Burstein,
Avishai Dekel,
Michael A. Strauss
Abstract:
This is the third in a series of papers in which we assemble and analyze a homogeneous catalog of peculiar velocity data. In Papers I and II, we described the Tully-Fisher (TF) redshift-distance samples that constitute the bulk of the catalog, and our methodology for obtaining mutually consistent TF calibrations for these samples. In this paper, we supply further technical details of the treatme…
▽ More
This is the third in a series of papers in which we assemble and analyze a homogeneous catalog of peculiar velocity data. In Papers I and II, we described the Tully-Fisher (TF) redshift-distance samples that constitute the bulk of the catalog, and our methodology for obtaining mutually consistent TF calibrations for these samples. In this paper, we supply further technical details of the treatment of the data, and present a subset of the catalog in tabular form. The full catalog, known as the Mark III Catalog of Galaxy Peculiar Velocities, is available in accessible on-line databases, as described herein. The electronic catalog incorporates not only the TF samples discussed in Papers I and II, but also elliptical galaxy Dn-sigma samples originally presented elsewhere. The relative zeropointing of the elliptical and spiral data sets is discussed here.
The basic elements of the Mark III Catalog are the observables for each object (redshift, magnitude, velocity width, etc.) and inferred distances derived from the TF or Dn-sigma relations. Malmquist bias-corrected distances are computed for each catalog object using density fields obtained from the IRAS 1.2 Jy redshift survey. Distances for both individual objects and groups are provided. A variety of auxiliary data, including distances and local densities predicted from the IRAS redshift survey reconstruction method, are tabulated as well. We study the distributions of TF residuals for three of our samples and conclude that they are well-approximated as Gaussian. For the Mathewson et al. sample we demonstrate a significant decrease in TF scatter with increasing velocity width. We test for, but find no evidence of, a correlation between TF residuals and galaxy morphology.
△ Less
Submitted 25 October, 1996;
originally announced October 1996.
-
From Local Velocities to Microwave Background
Authors:
S. Zaroubi,
N. Sugiyama,
J. Silk,
Y. Hoffman,
A. Dekel
Abstract:
The mass density field as extracted from peculiar velocities in our cosmological neighborhood is mapped back in time to the CMB in two ways. First, the density power spectrum ($P_k$) is translated into a temperature angular power spectrum of sub-degree resolution ($C_l$) and compared to observations. Second, the local density field is translated into a temperature map in a patch on the last-scat…
▽ More
The mass density field as extracted from peculiar velocities in our cosmological neighborhood is mapped back in time to the CMB in two ways. First, the density power spectrum ($P_k$) is translated into a temperature angular power spectrum of sub-degree resolution ($C_l$) and compared to observations. Second, the local density field is translated into a temperature map in a patch on the last-scattering surface of a distant observer. A likelihood analysis of the Mark III peculiar velocity data have constrained the range of parameters for $P_k$ within the family of COBE-normalized CDM models (Zaroubi et al 1996), favoring a slight tilt in the initial spectrum, $n<1$. The corresponding range of $C_l$'s is plotted against current observations, indicating that the CMB data can tighten the constraints further: only models with ``small'' tilt ($n\sim 0.9$) and ``high'' baryonic content ($Ω_b \sim 0.1$) could survive the two data sets simultaneously. The local mass density field that has been recovered from the velocities via a Wiener method is convolved with a Boltzmann calculation to recover $10'$ resolution temperature maps as viewed from different directions. The extent of the CMB patch and the amplitude of fluctuations depend on the choice of cosmological parameters, e.g., the local $100\hmpc$ sphere corresponds to $90'$ to $30'$ at the CMB for $Ω$ between 1 and 0 respectively. The phases of the temperature map are correlated with those of the density field, contrary to the contribution of the Sachs-Wolfe effect alone. This correlation suggests the possibility of an inverse reconstruction of the underlying density field from CMB data with interesting theoretical implications.
△ Less
Submitted 16 October, 1996;
originally announced October 1996.
-
`Comparing the Predicted Abell/ACO Cluster & the MarkIII Galaxy Density and Velocity Fields', Contribution to the Wide Field Spectroscopy Meeting held in Athens, 20-24 May 1996
Authors:
M. Plionis,
E. Branchini,
I. Zehavi,
A. Dekel
Abstract:
We compare the predicted Abell/ACO cluster density and velocity fields, using linear theory, with that based on the MarkIII catalog and the POTENT machinery within 7000 km/sec. Both density & velocity fields are in good qualitative agreement indicating that apart from a linear biasing factor the optical Abell/ACO clusters do trace the underlying mass distribution quite well.
We compare the predicted Abell/ACO cluster density and velocity fields, using linear theory, with that based on the MarkIII catalog and the POTENT machinery within 7000 km/sec. Both density & velocity fields are in good qualitative agreement indicating that apart from a linear biasing factor the optical Abell/ACO clusters do trace the underlying mass distribution quite well.
△ Less
Submitted 20 September, 1996;
originally announced September 1996.
-
Large Scale Power Spectrum from Peculiar Velocities Via Likelihood Analysis
Authors:
S. Zaroubi,
A. Dekel,
Y. Hoffman,
T. Kolatt
Abstract:
The power spectrum (PS) of {\it mass} density fluctuations, independent of ``biasing", is estimated from the Mark3 Catalog of Peculiar Velocities of galaxies using Bayesian statistics. A parametric model is assumed for the PS, and the free parameters are determined by maximizing the posterior probability of the model given the data. The method has been tested using detailed mock catalogs. We use…
▽ More
The power spectrum (PS) of {\it mass} density fluctuations, independent of ``biasing", is estimated from the Mark3 Catalog of Peculiar Velocities of galaxies using Bayesian statistics. A parametric model is assumed for the PS, and the free parameters are determined by maximizing the posterior probability of the model given the data. The method has been tested using detailed mock catalogs. We use generalized CDM models with and without COBE normalization.
The robust result for all the models tested is $P(k) Ω^{1.2} = (4.1 \pm 0.7)x10^3 (Mpc h^{-1})^3$ at $k=0.1 h Mpc^{-1}$, with the peak constrained to the range $0.03 < k < 0.06 h Mpc^{-1}$. It is consistent with a direct computation of the PS (Kolatt \& Dekel 1996). When compared to galaxy-density surveys, the implied values for $β$ ($\equiv Ω^{0.6}/b$) are of order unity to within 25\%.
A $Γ$-shape model, free of COBE normalization, is constrained by the velocity data to $Γ=0.5\pm 0.15$ and $σ_8 Ω^{0.6}=0.85\pm 0.1$. Within the family of COBE-normalized CDM models, the best tilted model ($Ω=1$, $n \approx 0.84 h_{50}^{-0.65}$) and the best open model ($n=1$, $Ω\approx 0.75 h_{50}^{-1.0}$) are more likely than the best $Λ$ model ($n=1$, $Λ=1-Ω$, $Ω\approx 0.70 h_{50}^{-1.2}$). The most likely CDM model with $Ω\leq 1$ is found to be of $Ω=1$, $h\approx 0.75$, and a {\it tilted} spectrum of $n= 0.8\pm 0.02$ with tensor fluctuations.The tightest constraint for the tilted-$Λ$ models is of the sort $Ωh_{50}^{1.2} n^ν= 0.7 \pm 0.08$, with $ν=3.8$ and $1.85$ with and without tensor fluctuations respectively.
△ Less
Submitted 14 March, 1996;
originally announced March 1996.
-
Large-Scale Power-Spectrum from Peculiar Velocities
Authors:
Tsafrir Kolatt,
Avishai Dekel
Abstract:
The power spectrum (PS) of mass density fluctuations, in the range $0.05 \leq k \leq 0.2 \ihmpc$, is derived from the Mark III catalog of peculiar velocities of galaxies, independent of ``biasing". It is computed from the density field as recovered by POTENT with Gaussian smoothing of $12\hmpc$, within a sphere of radius $\sim 60\hmpc$ about the Local Group. The density is weighted inversely by…
▽ More
The power spectrum (PS) of mass density fluctuations, in the range $0.05 \leq k \leq 0.2 \ihmpc$, is derived from the Mark III catalog of peculiar velocities of galaxies, independent of ``biasing". It is computed from the density field as recovered by POTENT with Gaussian smoothing of $12\hmpc$, within a sphere of radius $\sim 60\hmpc$ about the Local Group. The density is weighted inversely by the square of the errors. The PS is corrected for the effects of smoothing, random errors, sparse sampling and finite volume using mock catalogs that mimic in detail the Mark III catalog and the dynamics of our cosmological neighborhood. The mock catalogs are also used for error analysis. The PS at $k = 0.1 \ihmpc$ (for example) is $(4.6 \pm 1.4) \times 10^3 Ω^{-1.2} \3hmpc$, with a local logarithmic slope of $-1.45 \pm 0.5$. An itegration yields $σ_8 Ω^{0.6} \simeq 0.7-0.8$, depending on where the PS peak is. Direct comparisons of the mass PS with the galaxy PS derived from sky and redshift surveys show a similarity in shape and yield for $β\equiv Ω^{0.6}/b$ values in the range $0.77-1.21$, with a typical error of $\pm 0.1$ per galaxy sample. A comparison of the mass PS at $\sim 100\hmpc$ with the large-angle CMB fluctuations by COBE provides constraints on cosmological parameters and on the slope of the initial PS ($n$). The ``standard" CDM model is marginally rejected at the $\sim 2σ$ level, while each of the following modifications leads to a good fit: $n \lsim 1$, $Ω_ν\sim 0.3$, or $Ω\lsim 1$. Values of $Ω$ as low as $\sim 0.2$ are ruled out with high confidence (independent of $Λ$).
△ Less
Submitted 8 November, 1996; v1 submitted 19 December, 1995;
originally announced December 1995.
-
Simulating Our Cosmological Neighborhood: Mock Catalogs for Velocity Analysis
Authors:
Tsafrir Kolatt,
Avishai Dekel,
Galit Ganon,
Jeffrey Willick
Abstract:
We describe the construction of an N-body simulation that mimics the true velocity and mass-density fields in a box of side 256\hmpc about the Local Group , and the production of mock catalogs that mimic in detail current catalogs of redshifts and peculiar velocities. Our main purpose is to provide a tool for developing and testing reconstruction methods, but the different components of the meth…
▽ More
We describe the construction of an N-body simulation that mimics the true velocity and mass-density fields in a box of side 256\hmpc about the Local Group , and the production of mock catalogs that mimic in detail current catalogs of redshifts and peculiar velocities. Our main purpose is to provide a tool for developing and testing reconstruction methods, but the different components of the method can be used on their own in other applications. The initial conditions are based on the IRAS 1.2Jy redshift survey, assuming that galaxies trace mass and Ω=1. A density field smoothed is recovered from the redshift survey, using quasi-linear theory and a power-preserving filter. The corresponding potential field is traced back to the linear regime using the Zel'dovich-Bernoulli equation. Small-scale power is added by means of constrained realization to mimic fluctuations on galactic scales. The gravitating system is evolved forward in time with a PM code of 2\hmpc resolution. The result reproduces the real dynamical structures on large scales and the statistical properties of the structure down to galactic scales. "Galaxies" are identified via a linear biasing scheme (b=1.35) and they are divided into "spirals" and "ellipticals" to obey Dressler's morphology-density relation. The galaxies are assigned internal-velocity parameters and absolute magnitudes scattered about an assumed mean Tully-Fisher relation. They are then "observed" as magnitude-limited samples, trying to mimic the selection criteria of the data sets constituting the Mark III catalog of peculiar velocities. Artificial IRAS 1.2Jy redshift surveys are also compiled. The simulations and mock catalogs will be made available electronically as bench marks for testing reconstruction methods.
△ Less
Submitted 12 September, 1995;
originally announced September 1995.