Linearisation with Cosmological Perturbation Theory

We propose a new method to linearise cosmological mass density fields using higher order Lagrangian perturbation theory (LPT). We demonstrate that a given density field can be expressed as the sum of a linear and a nonlinear component which are tightly coupled to each other by the tidal field tensor within the LPT framework. The linear component corresponds to the initial density field in Eulerian coordinates, and its mean relation with the total field can be approximated by a logarithm (giving theoretical support to recent attempts to find such component). We also propose to use a combination of the linearisation method and the continuity equation to find the mapping between Eulerian and Lagrangian coordinates. In addition, we note that this method opens the possibility of use directly higher order LPT on nonlinear fields. We test our linearization scheme by applying it to the z~0.5 density field from an N-body simulation. We find that the linearised version of the full density field can be successfully recovered on >~5 h^{-1}Mpc, reducing the skewness and kurtosis of the distribution by about one and two orders of magnitude, respectively. This component can also be successfully traced back in time, converging towards the initial unevolved density field at z~100. We anticipate a number of applications of our results, from predicting velocity fields to estimates of the initial conditions of the universe, passing by improved constraints on cosmological parameters derived from galaxy clustering via reconstruction methods.Comment: 14 pages, 8 figure

The Initial Conditions of the Universe from Constrained Simulations

I present a new approach to recover the primordial density fluctuations and the cosmic web structure underlying a galaxy distribution. The method is based on sampling Gaussian fields which are compatible with a galaxy distribution and a structure formation model. This is achieved by splitting the inversion problem into two Gibbs-sampling steps: the first being a Gaussianisation step transforming a distribution of point sources at Lagrangian positions -which are not a priori given- into a linear alias-free Gaussian field. This step is based on Hamiltonian sampling with a Gaussian-Poisson model. The second step consists on a likelihood comparison in which the set of matter tracers at the initial conditions is constrained on the galaxy distribution and the assumed structure formation model. For computational reasons second order Lagrangian Perturbation Theory is used. However, the presented approach is flexible to adopt any structure formation model. A semi-analytic halo-model based galaxy mock catalog is taken to demonstrate that the recovered initial conditions are closely unbiased with respect to the actual ones from the corresponding N-body simulation down to scales of a ~ 5 Mpc/h. The cross-correlation between them shows a substantial gain of information, being at k ~ 0.3 h/Mpc more than doubled. In addition the initial conditions are extremely well Gaussian distributed and the power-spectra follow the shape of the linear power-spectrum being very close to the actual one from the simulation down to scales of k ~ 1 h/Mpc.Comment: 6 pages, 4 figures; accepted at MNRAS after minor correction

Itinerant-localized dual character of a strongly-correlated superfluid Bose gas in an optical lattice

We investigate a strongly-correlated Bose gas in an optical lattice. Extending the standard-basis operator method developed by Haley and Erdos to a boson Hubbard model, we calculate excitation spectra in the superfluid phase, as well as in the Mott insulating phase, at T=0. In the Mott phase, the excitation spectrum has a finite energy gap, reflecting the localized character of atoms. In the superfluid phase, the excitation spectrum is shown to have an itinerant-localized dual structure, where the gapless Bogoliubov mode (which describes the itinerant character of superfluid atoms) and a band with a finite energy gap coexist. We also show that the rf-tunneling current measurement would give a useful information about the duality of a strongly-correlated superfluid Bose gas near the superfluid-insulator transition.Comment: 10 pages, 4 figure

Neutrino-driven explosions twenty years after SN1987A

The neutrino-heating mechanism remains a viable possibility for the cause of the explosion in a wide mass range of supernova progenitors. This is demonstrated by recent two-dimensional hydrodynamic simulations with detailed, energy-dependent neutrino transport. Neutrino-driven explosions were not only found for stars in the range of 8-10 solar masses with ONeMg cores and in case of the iron core collapse of a progenitor with 11 solar masses, but also for a ``typical'' progenitor model of 15 solar masses. For such more massive stars, however, the explosion occurs significantly later than so far thought, and is crucially supported by large-amplitude bipolar oscillations due to the nonradial standing accretion shock instability (SASI), whose low (dipole and quadrupole) modes can develop large growth rates in conditions where convective instability is damped or even suppressed. The dominance of low-mode deformation at the time of shock revival has been recognized as a possible explanation of large pulsar kicks and of large-scale mixing phenomena observed in supernovae like SN 1987A.Comment: 11 pages, 6 figures; review proceeding for "Supernova 1987A: 20 Years After: Supernovae and Gamma-Ray Bursters" AIP, New York, eds. S. Immler, K.W. Weiler, and R. McCra

Cosmic flows and the expansion of the Local Universe from nonlinear phase-space reconstructions

This work investigates the impact of cosmic flows and density perturbations on Hubble constant $H_0$ measurements using nonlinear phase-space reconstructions of the Local Universe (LU). In particular, we rely on 25 constrained N-body simulations of the LU using the 2MRS galaxy sample within distances of 90 Mpc/h. These have been randomly extended up to volumes enclosing 360 Mpc/h with augmented Lagrangian perturbation theory (750 simulations in total), accounting in this way for gravitational mode coupling from larger scales, correcting for periodic boundary effects, and estimating systematics of missing attractors ($\sigma_{\rm large}=134$ km/s). We report on Local Group speed reconstructions, which are compatible with CMB-dipole measurements: $|v_{\rm LG}|=685\pm137$ km/s. The direction $(l,b)=(260.5\pm 13.3,39.1\pm 10.4)^\circ$ is compatible with observations after considering the variance of large scales. Accounting for large scales, our local bulk flow estimations assuming a $\Lambda$CDM model are compatible with estimates based on velocity data derived from the Tully-Fisher relation. We focus on low redshift supernova measurements ($0.01<z<0.025$), which have been found to disagree with probes at larger distances. Our analysis indicates that there are two effects contributing to this tension. First, the anisotropic distribution of supernovae aligns with the velocity dipole and induces a systematic boost in $H_0$. Second, a divergent region surrounding the Virgo Supercluster is responsible for an additional positive bias in $H_0$. Taking these effects into account yields a correction of $\Delta H_0=-1.76 \pm 0.21$ km/s/Mpc, thereby reducing the tension between local probes and more distant probes. Effectively $H_0$ is lower by about $2\%$.Comment: accepted by MNRAS 10 pages, 7 figure

Cosmological Structure Formation with Augmented Lagrangian Perturbation Theory

We present a new fast and efficient approach to model structure formation with Augmented Lagrangian Perturbation Theory (ALPT). Our method is based on splitting the displacement field into a long and a short-range component. The long-range component is computed by second order LPT (2LPT). This approximation contains a tidal nonlocal and nonlinear term. Unfortunately, 2LPT fails on small scales due to severe shell crossing and a crude quadratic behaviour in the low density regime. The spherical collapse (SC) approximation has been recently reported to correct for both effects by adding an ideal collapse truncation. However, this approach fails to reproduce the structures on large scales where it is significantly less correlated with the N-body result than 2LPT or linear LPT (the Zeldovich approximation). We propose to combine both approximations using for the short-range displacement field the SC solution. A Gaussian filter with a smoothing radius r_S is used to separate between both regimes. We use the result of 25 dark matter only N-body simulations to benchmark at z=0 the different approximations: 1st, 2nd, 3rd order LPT, SC and our novel combined ALPT model. This comparison demonstrates that our method improves previous approximations at all scales showing ~25% and ~75% higher correlation than 2LPT with the N-body solution at k = 1 and 2 h Mpc^-1, respectively. We conduct a parameter study to determine the optimal range of smoothing radii and find that the maximum correlation is achieved with r_S = 4 - 5 h^-1 Mpc. This structure formation approach could be used for various purposes, such as setting-up initial conditions for N-body simulations, generating mock galaxy catalogues, cosmic web analysis or for reconstructions of the primordial density fluctuations.Comment: 6 pages and 4 figure

Modelling Baryon Acoustic Oscillations with Perturbation Theory and Stochastic Halo Biasing

In this work we investigate the generation of mock halo catalogues based on perturbation theory and nonlinear stochastic biasing with the novel PATCHY-code. In particular, we use Augmented Lagrangian Perturbation Theory (ALPT) to generate a dark matter density field on a mesh starting from Gaussian fluctuations and to compute the peculiar velocity field. ALPT is based on a combination of second order LPT (2LPT) on large scales and the spherical collapse model on smaller scales. We account for the systematic deviation of perturbative approaches from N-body simulations together with halo biasing adopting an exponential bias model. We then account for stochastic biasing by defining three regimes: a low, an intermediate and a high density regime, using a Poisson distribution in the intermediate regime and the negative binomial distribution to model over-dispersion in the high density regime. Since we focus in this study on massive halos, we suppress the generation of halos in the low density regime. The various nonlinear and stochastic biasing parameters, and density thresholds (five) are calibrated with the large BigMultiDark N-body simulation to match the power spectrum of the corresponding halo population. Our mock catalogues show power spectra, both in real- and redshift-space, which are compatible with N-body simulations within about 2% up to k ~ 1 h Mpc^-1 at z = 0.577 for a sample of halos with the typical BOSS CMASS galaxy number density. The corresponding correlation functions are compatible down to a few Mpc. We also find that neglecting over-dispersion in high density regions produces power spectra with deviations of 10% at k ~ 0.4 h Mpc^-1. These results indicate the need to account for an accurate statistical description of the galaxy clustering for precise studies of large-scale surveys.Comment: 5 pages, 4 figure