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 in our test simulation. For realistic sparse and noisy data, this median increases to 5 Mpc/h. This is a significant improvement over the previous approach of neglecting the displacement field, which introduces errors on a scale of 10 Mpc/h 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.

The Reverse Zeldovich Approximation (RZA) is a reconstruction method which allows 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%), the sparseness of the data points, and the maximum catalogue radius (3000 - 6000 km/s). 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 dataset 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 if when data is 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 dataset with high observational errors still leads to a good reconstruction of the initial conditions on a scale of about 5 Mpc/h.

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 (CR) method. Our re-simulations 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 CR 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 re-simulations with an average accuracy of about 2 Mpc/h on their position and a factor of 2 in mass, down to haloes with a mass of approx 10^14 M_odot/h. In comparison, without RZA the re-simulations recover only the most massive haloes with masses of about 5x10^14 M_odot/h and higher, and with a systematic shift on their position of about approx 10 Mpc/h due to the cosmic displacement field. We show that with the additional Lagrangian reconstruction step introduced by the RZA, this shift can be removed.

Using a dark matter only Constrained Local UniversE Simulation (CLUES) we examine the existence of subhaloes that change their affiliation from one of the two prominent hosts in the Local Group (i.e. the Milky Way and the Andromeda galaxy) to the other, and call these objects "renegade subhaloes". In light of recent claims that the two Magellanic Clouds (MCs) may have originated from another region (or even the outskirts) of the Local Group or that they have been spawned by a major merger in the past of the Andromeda galaxy, we investigate the nature of such events. However, we cannot confirm that renegade subhaloes enter as deep into the potential well of their present host nor that they share the most simplest properties with the MCs, namely mass and relative velocity. Our simulation rather suggests that these renegade subhaloes appear to be flying past one host before being pulled into the other. A merger is not required to trigger such an event, it is rather the distinct environment of our simulated Local Group facilitating such behavior. Since just a small fraction of the full z=0 subhalo population are renegades, our study indicates that it will be intrinsically difficult to distinguish them despite clear differences in their velocity, radial distribution, shape and spin parameter distributions.

Despite greatly improved observational methods, the presence of magnetic fields at cosmological scales and their role in the process of large-scale structure formation still remains unclear. In this paper we want to address the question how the presence of a hypothetical primordial magnetic field on large scales influences the cosmic structure formation in numerical simulations. As a tool for carrying out such simulations, we present our new numerical code AMIGA. It combines an N-body code with an Eulerian grid-based solver for the full set of MHD equations in order to conduct simulations of dark matter, baryons and magnetic fields in a self-consistent way in a fully cosmological setting. Our numerical scheme includes effective methodes to ensure proper capturing of shocks and highly supersonic flows and a divergence-free magnetic field. The high accuracy of the code is demonstrated by a number of numerical tests. We then present a series of cosmological MHD simulations and confirm that, in order to have a significant effect on the distribution of matter on large scales, the primordial magnetic field strength would have to be significantly higher than the current observational and theoretical constraints.