In previous works we proposed the Reverse Zeldovich Approximation (RZA) method, which can be used to estimate the cosmological initial conditions underlying the galaxy distribution in the Local Universe using peculiar velocity data. In this paper, we apply the technique to run constrained cosmological simulations from the RZA-reconstructed initial conditions, designed to reproduce the large-scale structure of the Local Universe. We test the method with mock peculiar velocity catalogues extracted from a reference simulation. We first reconstruct the initial conditions of this reference simulation using the mock data, and then run the reconstructed initial conditions forward in time until z = 0. We compare the resulting constrained simulations with the original simulation at z = 0 to test the accuracy of this method. We also compare them with constrained simulations run from the mock data without the addition of RZA, i.e. using only the currently established constrained realizations (CRs) method. Our resimulations are able to correctly recover the evolution of the large-scale structure underlying the data. The results show that the addition of RZA to the CRs method significantly improves both the reconstruction of the initial conditions and the accuracy of the obtained constrained resimulations. Haloes from the original simulation are recovered in the resimulations with an average accuracy of ≈2 Mpc h-1 on their position and a factor of 2 in mass, down to haloes with a mass of ≈1014 M⊙ h-1. In comparison, without RZA the resimulations recover only the most massive haloes with masses of ≈5 × 1014 M⊙ h-1 and higher, and with a systematic shift on their position of about ≈ 10 Mpc h- 1 due to the cosmic displacement field. We show that with the additional Lagrangian reconstruction step introduced by the RZA, this shift can be removed.
The Reverse Zeldovich Approximation (RZA) is a reconstruction method which allows us to estimate the cosmic displacement field from galaxy peculiar velocity data and to constrain initial conditions for cosmological simulations of the Local Universe. In this paper, we investigate the effect of different observational errors on the reconstruction quality of this method. For this, we build a set of mock catalogues from a cosmological simulation, varying different error sources like the galaxy distance measurement error (0-20 per cent), the sparseness of the data points and the maximum catalogue radius (3000-6000 km s-1). We perform the RZA reconstruction of the initial conditions on these mock catalogues and compare with the actual initial conditions of the simulation. We also investigate the impact of the fact that only the radial part of the peculiar velocity is observationally accessible. We find that the sparseness of a data set has the highest detrimental effect on RZA reconstruction quality. Observational distance errors also have a significant influence, but it is possible to compensate this relatively well with Wiener filter reconstruction. We also investigate the effect of different object selection criteria and find that distance catalogues distributed randomly and homogeneously across the sky (such as spiral galaxies selected for the Tully-Fisher method) allow for a higher reconstruction quality than in the case when data are preferentially drawn from massive objects or dense environments (such as elliptical galaxies). We find that the error of estimating the initial conditions with RZA is always dominated by the inherent non-linearity of data observed at z = 0 rather than by the combined effect of the observational errors. Even an extremely sparse data set with high observational errors still leads to a good reconstruction of the initial conditions on a scale of ≈5 Mpc h-1.
We propose a new method to recover the cosmological initial conditions of the presently observed galaxy distribution, which can serve to run constrained simulations of the Local Universe. Our method, the Reverse Zeldovich Approximation (RZA), can be applied to radial galaxy peculiar velocity data and extends the previously used constrained realizations (CR) method by adding a Lagrangian reconstruction step. The RZA method consists of applying the Zeldovich approximation in reverse to galaxy peculiar velocities to estimate the cosmic displacement field and the initial linear matter distribution from which the present-day Local Universe evolved. We test our method with a mock survey taken from a cosmological simulation. We show that the halo peculiar velocities at z = 0 are close to the linear prediction of the Zeldovich approximation, if a grouping is applied to the data to remove virial motions. We find that the addition of RZA to the CR method significantly improves the reconstruction of the initial conditions. The RZA is able to recover the correct initial positions of the velocity tracers with a median error of only 1.36 Mpc h-1 in our test simulation. For realistic sparse and noisy data, this median increases to 5 Mpc h-1. This is a significant improvement over the previous approach of neglecting the displacement field, which introduces errors on a scale of 10 Mpc h-1 or even higher. Applying the RZA method to the upcoming high-quality observational peculiar velocity catalogues will generate much more precise constrained simulations of the Local Universe.
Benítez-Llambay, A., Navarro, J. F., Abadi, M. G., Gottlöber, S., Yepes, G., Hoffman, Y., Steinmetz, M., 2013, The Astrophysical Journal
, 763, 2 , L41 Published: February 2013
We use a cosmological simulation of the formation of the Local Group of Galaxies to identify a mechanism that enables the removal of baryons from low-mass halos without appealing to feedback or reionization. As the Local Group forms, matter bound to it develops a network of filaments and pancakes. This moving web of gas and dark matter drifts and sweeps a large volume, overtaking many halos in the process. The dark matter content of these halos is unaffected but their gas can be efficiently removed by ram pressure. The loss of gas is especially pronounced in low-mass halos due to their lower binding energy and has a dramatic effect on the star formation history of affected systems. This "cosmic web stripping" may help to explain the scarcity of dwarf galaxies compared with the numerous low-mass halos expected in ΛCDM and the large diversity of star formation histories and morphologies characteristic of faint galaxies. Although our results are based on a single high-resolution simulation, it is likely that the hydrodynamical interaction of dwarf galaxies with the cosmic web is a crucial ingredient so far missing from galaxy formation models.
Libeskind, N. I., Hoffman, Y., Forero-Romero, J., Gottlober, S., Knebe, A., Steinmetz, M., Klypin, A., 2013, Monthly Notices of the Royal Astronomical Society
, 428, 3 , 2489 Published: January 2013
The alignment of dark matter (DM) haloes and the surrounding large-scale structure (LSS) is examined in the context of the cosmic web. Halo spin, shape and the orbital angular momentum of subhaloes are investigated relative to the LSS using the eigenvectors of the velocity shear tensor evaluated on a grid with a scale of 1 Mpc h-1, deep within the non-linear regime. Knots, filaments, sheets and voids are associated with regions that are collapsing along 3, 2, 1 or 0 principal directions simultaneously. Each halo is tagged with a web classification (i.e. knot halo, filament halo, etc.) according to the nature of the collapse at the halo position. The full distribution of shear eigenvalues is found to be substantially different from that tagged to haloes, indicating that the observed velocity shear is significantly biased. We find that larger mass haloes live in regions where the shear is more isotropic, namely the expansion or collapse is more spherical. A correlation is found between the halo shape and the eigenvectors of the shear tensor, with the longest (shortest) axis of the halo shape being aligned with the slowest (fastest) collapsing eigenvector. This correlation is web independent, suggesting that the velocity shear is a fundamental tracer of the halo alignment. A similar result is found for the alignment of halo spin with the cosmic web. It has been shown that high-mass haloes exhibit a spin flip with respect to the LSS: we find that the mass at which this spin flip occurs is web dependent and not universal as suggested previously. Although weaker than haloes, subhalo orbits too exhibit an alignment with the LSS, providing a possible insight into the highly correlated corotation of the Milky Way's satellite system. The present study suggests that the velocity shear tensor constitutes the natural framework for studying the directional properties of the non-linear LSS and those of haloes and galaxies.
