Lagrangian Perturbation Theory Research Articles

Reignited star formation in dwarf galaxies that were quenched during reionization

Context. Irregular dwarf galaxies of the Local Group have very varied properties and star formation histories. Some of them formed the majority of their stars very late compared to others. Extreme examples of this are Leo A and Aquarius, which reached the peak of star formation at z < 1 (more than 6 Gyr after the Big Bang). This fact seemingly challenges the ΛCDM cosmological framework because the dark matter halos of these galaxies on average should assemble the majority of their masses before z ~ 2 (<3 Gyr after the Big Bang). Aims. We investigate whether the delayed star formation histories of some irregular dwarf galaxies might be explained purely by the stochasticity of their mass assembly histories coupled with the effect of cosmic reionization. Methods. We developed a semi-analytic model to follow the accretion of baryonic matter, star formation, and stellar feedback in dark matter halos with present-day virial masses 109 M⊙ < Mdm,0 < 1011 M⊙ and with different stochastic growth histories obtained using the PINOCCHIO code based on Lagrangian perturbation theory. Results. We obtain the distributions of observable parameters and the evolution histories for these galaxies. Accretion of baryonic matter is strongly suppressed after the epoch of reionization in some models, but the galaxies continue to accrete dark matter and eventually reach enough mass for accretion of baryonic matter to begin again. These “reborn” model galaxies show delayed star formation histories that are very similar to those of Leo A and Aquarius. Conclusions. We find that the stochasticity caused by mass assembly histories is enhanced in systems with virial masses ~1010 M⊙ because of their sensitivity to the photoionizing intergalactic radiation field after the epoch of reionization. This results in qualitatively different star formation histories in late- and early-forming galaxies, and it might explain the peculiar star formation histories of irregular dwarf galaxies such as Leo A and Aquarius.

Field theoretical derivation of L\xfcscher\u2019s formula and calculation of finite volume form factors

We initiate a systematic method to calculate both the finite volume energy levels and form factors from the momentum space finite volume two-point function. By expanding the two point function in the volume we extracted the leading exponential volume correction both to the energy of a moving particle state and to the simplest non-diagonal form factor. The form factor corrections are given in terms of a regularized infinite volume 3-particle form factor and terms related to the Lüsher correction of the momentum quantization. We tested these results against second order Lagrangian and Hamiltonian perturbation theory in the sinh-Gordon theory and we obtained perfect agreement.

Stability Analysis of Fluid Flows Using Lagrangian Perturbation Theory (LPT): Application to the Plane Couette Flow

We present a new application of Lagrangian Perturbation Theory (LPT): the stability analysis of fluid flows. As a test case that demonstrates the framework we focus on the plane Couette flow. The incompressible Navier-Stokes equation is recast such that the particle position is the fundamental variable, expressed as a function of Lagrangian coordinates. The displacement due to the steady state flow is taken to be the zeroth order solution and the position is formally expanded in terms of a small parameter (generally, the strength of the initial perturbation). The resulting hierarchy of equations is solved analytically at first order. We find that we recover the standard result in the Eulerian frame: the plane Couette flow is asymptotically stable for all Reynolds numbers. However, it is also well established that experiments contradict this prediction. In the Eulerian picture, one of the proposed explanations is the phenomenon of `transient growth' which is related to the non-normal nature of the linear stability operator. The first order solution in the Lagrangian frame also shows this feature, albeit qualitatively. As a first step, and for the purposes of analytic manipulation, we consider only linear stability of 2D perturbations but the framework presented is general and can be extended to higher orders, other flows and/or 3D perturbations.

The clustering of the SDSS-IV extended Baryon Oscillation Spectroscopic Survey DR14 quasar sample: measurement of the growth rate of structure from the anisotropic correlation function between redshift 0.8 and 2.2

We present the clustering measurements of quasars in configuration space based on the Data Release 14 (DR14) of the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey. This dataset includes 148,659 quasars spread over the redshift range $0.8\leq z \leq 2.2$ and spanning 2112.9 square degrees. We use the Convolution Lagrangian Perturbation Theory (CLPT) approach with a Gaussian Streaming (GS) model for the redshift space distortions of the correlation function and demonstrate its applicability for dark matter halos hosting eBOSS quasar tracers. At the effective redshift $z_{\rm eff} = 1.52$, we measure the linear growth rate of structure $f\sigma_{8}(z_{\rm eff})= 0.426 \pm 0.077$, the expansion rate $H(z_{\rm eff})= 159^{+12}_{-13}(r_{s}^{\rm fid}/r_s){\rm km.s}^{-1}.{\rm Mpc}^{-1}$, and the angular diameter distance $D_{A}(z_{\rm eff})=1850^{+90}_{-115}\,(r_s/r_{s}^{\rm fid}){\rm Mpc}$, where $r_{s}$ is the sound horizon at the end of the baryon drag epoch and $r_{s}^{\rm fid}$ is its value in the fiducial cosmology. The quoted errors include both systematic and statistical contributions. The results on the evolution of distances are consistent with the predictions of flat $\Lambda$-Cold Dark Matter ($\Lambda$-CDM) cosmology with Planck parameters, and the measurement of $f\sigma_{8}$ extends the validity of General Relativity (GR) to higher redshifts($z>1$) This paper is released with companion papers using the same sample. The results on the cosmological parameters of the studies are found to be in very good agreement, providing clear evidence of the complementarity and of the robustness of the first full-shape clustering measurements with the eBOSS DR14 quasar sample.

Symmetries, invariants and generating functions: higher-order statistics of biased tracers

Gravitationally collapsed objects are known to be biased tracers of an underlying density contrast. Using symmetry arguments, generalised biasing schemes have recently been developed to relate the halo density contrast δh with the underlying density contrast δ, divergence of velocity θ and their higher-order derivatives. This is done by constructing invariants such as s, t, ψ,η. We show how the generating function formalism in Eulerian standard perturbation theory (SPT) can be used to show that many of the additional terms based on extended Galilean and Lifshitz symmetry actually do not make any contribution to the higher-order statistics of biased tracers. Other terms can also be drastically simplified allowing us to write the vertices associated with δh in terms of the vertices of δ and θ, the higher-order derivatives and the bias coefficients. We also compute the cumulant correlators (CCs) for two different tracer populations. These perturbative results are valid for tree-level contributions but at an arbitrary order. We also take into account the stochastic nature bias in our analysis. Extending previous results of a local polynomial model of bias, we express the one-point cumulants \U0001d4aeN and their two-point counterparts, the CCs i.e. \U0001d49epq, of biased tracers in terms of that of their underlying density contrast counterparts. As a by-product of our calculation we also discuss the results using approximations based on Lagrangian perturbation theory (LPT).

Lagrangian theory of structure formation in relativistic cosmology. IV. Lagrangian approach to gravitational waves

The relativistic generalization of the Newtonian Lagrangian perturbation theory is investigated. In previous works, the perturbation and solution schemes that are generated by the spatially projected gravitoelectric part of the Weyl tensor were given to any order of the perturbations, together with extensions and applications for accessing the nonperturbative regime. We here discuss more in detail the general first-order scheme within the Cartan formalism including and concentrating on the gravitational wave propagation in matter. We provide master equations for all parts of Lagrangian-linearized perturbations propagating in the perturbed spacetime, and we outline the solution procedure that allows one to find general solutions. Particular emphasis is given to global properties of the Lagrangian perturbation fields by employing results of Hodge-de Rham theory. We here discuss how the Hodge decomposition relates to the standard scalar-vector-tensor decomposition. Finally, we demonstrate that we obtain the known linear perturbation solutions of the standard relativistic perturbation scheme by performing two steps: first, by restricting our solutions to perturbations that propagate on a flat unperturbed background spacetime and, second, by transforming to Eulerian background coordinates with truncation of nonlinear terms.

Lagrangian perturbation theory for modified gravity

We present a formalism to compute Lagrangian displacement fields for a wide range of cosmologies in the context of perturbation theory up to third order. We emphasize the case of theories with scale dependent gravitational strengths, such as chameleons, but our formalism can be accommodated to other modified gravity theories. In the non-linear regime two qualitative features arise. One, as is well known, is that nonlinearities lead to a screening of the force mediated by the scalar field. The second is a consequence of the transformation of the Klein-Gordon equation from Eulerian to Lagrangian coordinates, producing frame-lagging terms that are important especially at large scales, and if not considered, the theory does not reduce to the $\Lambda$CDM model in that limit. We apply our formalism to compute the 1-loop power spectrum and the correlation function in $f(R)$ gravity by using different resummation schemes. We further discuss the IR divergences of these formalisms.

COLA with massive neutrinos

The effect of massive neutrinos on the growth of cold dark matter perturbations acts as a scale-dependent Newton's constant and leads to scale-dependent growth factors just as we often find in models of gravity beyond General Relativity. We show how to compute growth factors for ΛCDM and general modified gravity cosmologies combined with massive neutrinos in Lagrangian perturbation theory for use in COLA and extensions thereof. We implement this together with the grid-based massive neutrino method of Brandbyge and Hannestad in MG-PICOLA and compare COLA simulations to full N-body simulations of ΛCDM and f(R) gravity with massive neutrinos. Our implementation is computationally cheap if the underlying cosmology already has scale-dependent growth factors and it is shown to be able to produce results that match N-body to percent level accuracy for both the total and CDM matter power-spectra up to k≲ 1 h/Mpc.

Accurate halo–galaxy mocks from automatic bias estimation and particle mesh gravity solvers

Reliable extraction of cosmological information from clustering measurements of galaxy surveys requires estimation of the error covariance matrices of observables. The accuracy of covariance matrices is limited by our ability to generate sufficiently large number of independent mock catalogs that can describe the physics of galaxy clustering across a wide range of scales. Furthermore, galaxy mock catalogs are required to study systematics in galaxy surveys and to test analysis tools. In this investigation, we present a fast and accurate approach for generation of mock catalogs for the upcoming galaxy surveys. Our method relies on low-resolution approximate gravity solvers to simulate the large scale dark matter field, which we then populate with halos according to a flexible nonlinear and stochastic bias model. In particular, we extend the \textsc{patchy} code with an efficient particle mesh algorithm to simulate the dark matter field (the \textsc{FastPM} code), and with a robust MCMC method relying on the \textsc{emcee} code for constraining the parameters of the bias model. Using the halos in the BigMultiDark high-resolution $N$-body simulation as a reference catalog, we demonstrate that our technique can model the bivariate probability distribution function (counts-in-cells), power spectrum, and bispectrum of halos in the reference catalog. Specifically, we show that the new ingredients permit us to reach percentage accuracy in the power spectrum up to $k\sim 0.4\; \,h\,{\rm Mpc}^{-1}$ (within 5\% up to $k\sim 0.6\; \,h\,{\rm Mpc}^{-1}$) with accurate bispectra improving previous results based on Lagrangian perturbation theory.

Time-dependent Lagrangian Perturbation Theory with The Dynamical Dark Energy

This article gives a brief overview of the time-dependent Lagrangian perturbation theory as a theoretical tool to predict the behavior of galaxy clustering in quasi-linear regimes. While most of the Lagrangian perturbation theories are based on an Einstein-De Sitter space where the growth function is time independent, more accurate theoretical predictions could be achieved by time-dependent growth functions in a more general cosmology. This is directly applicable to the dynamical dark-energy where the clustering behavior can be recognized with time-dependent growth functions.

COLA with scale-dependent growth: applications to screened modified gravity models

We present a general parallelized and easy-to-use code to perform numerical simulations of structure formation using the COLA (COmoving Lagrangian Acceleration) method for cosmological models that exhibit scale-dependent growth at the level of first and second order Lagrangian perturbation theory. For modified gravity theories we also include screening using a fast approximate method that covers all the main examples of screening mechanisms in the literature. We test the code by comparing it to full simulations of two popular modified gravity models, namely f(R) gravity and nDGP, and find good agreement in the modified gravity boost-factors relative to ΛCDM even when using a fairly small number of COLA time steps.

Testing approximate predictions of displacements of cosmological dark matter halos

We present a test to quantify how well some approximate methods, designed to reproduce the mildly non-linear evolution of perturbations, are able to reproduce the clustering of DM halos once the grouping of particles into halos is defined and kept fixed. The following methods have been considered: Lagrangian Perturbation Theory (LPT) up to third order, Truncated LPT, Augmented LPT, MUSCLE and COLA. The test runs as follows: halos are defined by applying a friends-of-friends (FoF) halo finder to the output of an N-body simulation. The approximate methods are then applied to the same initial conditions of the simulation, producing for all particles displacements from their starting position and velocities. The position and velocity of each halo are computed by averaging over the particles that belong to that halo, according to the FoF halo finder. This procedure allows us to perform a well-posed test of how clustering of the matter density and halo density fields are recovered, without asking to the approximate method an accurate reconstruction of halos. We have considered the results at z=0,0.5,1, and we have analysed power spectrum in real and redshift space, object-by-object difference in position and velocity, density Probability Distribution Function (PDF) and its moments, phase difference of Fourier modes. We find that higher LPT orders are generally able to better reproduce the clustering of halos, while little or no improvement is found for the matter density field when going to 2LPT and 3LPT. Augmentation provides some improvement when coupled with 2LPT, while its effect is limited when coupled with 3LPT. Little improvement is brought by MUSCLE with respect to Augmentation. The more expensive particle-mesh code COLA outperforms all LPT methods, and this is true even for mesh sizes as large as the inter-particle distance. This test sets an upper limit on the ability of these methods to reproduce the clustering of halos, for the cases when these objects are reconstructed at the object-by-object level.

Efficient simulations of large-scale structure in modified gravity cosmologies with comoving Lagrangian acceleration

We implement an adaptation of the COLA approach, a hybrid scheme that combines Lagrangian perturbation theory with an N-body approach, to model non-linear collapse in chameleon and symmetron modified gravity models. Gravitational screening is modeled effectively through the attachment of a suppression factor to the linearized Klein-Gordon equations. The adapted COLA approach is benchmarked, with respect to an N-body code both for the $\Lambda$CDM scenario and for the modified gravity theories. It is found to perform well in the estimation of the dark matter power spectra, with consistency of 1 % to $k\sim2.5$ h/Mpc. Redshift space distortions are shown to be effectively modeled through a Lorentzian parameterization with a velocity dispersion fit to the data. We find that COLA performs less well in predicting the halo mass functions, but has consistency, within $1\sigma$ uncertainties of our simulations, in the relative changes to the mass function induced by the modified gravity models relative to $\Lambda$CDM. The results demonstrate that COLA, proposed to enable accurate and efficient, non-linear predictions for $\Lambda$CDM, can be effectively applied to a wider set of cosmological scenarios, with intriguing properties, for which clustering behavior needs to be understood for upcoming surveys such as LSST, DESI, Euclid and WFIRST.

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: on the measurement of growth rate using galaxy correlation functions

We present a measurement of the linear growth rate of structure, \textit{f} from the Sloan Digital Sky Survey III (SDSS III) Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12 (DR12) using Convolution Lagrangian Perturbation Theory (CLPT) with Gaussian Streaming Redshift-Space Distortions (GSRSD) to model the two point statistics of BOSS galaxies in DR12. The BOSS-DR12 dataset includes 1,198,006 massive galaxies spread over the redshift range $0.2 < z < 0.75$. These galaxy samples are categorized in three redshift bins. Using CLPT-GSRSD in our analysis of the combined sample of the three redshift bins, we report measurements of $f \sigma_8$ for the three redshift bins. We find $f \sigma_8 = 0.430 \pm 0.054$ at $z_{\rm eff} = 0.38$, $f \sigma_8 = 0.452 \pm 0.057$ at $z_{\rm eff} = 0.51$ and $f \sigma_8 = 0.457 \pm 0.052$ at $z_{\rm eff} = 0.61$. Our results are consistent with the predictions of Planck $\Lambda$CDM-GR. Our constraints on the growth rates of structure in the Universe at different redshifts serve as a useful probe, which can help distinguish between a model of the Universe based on dark energy and models based on modified theories of gravity. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in Alam et al. 2016 to produce the final cosmological constraints from BOSS.

The integrated bispectrum and beyond

The position-dependent power spectrum has been recently proposed as a descriptor of gravitationally induced non-Gaussianity in galaxy clustering, as it is sensitive to the “soft limit” of the bispectrum (i.e. when one of the wave number tends to zero). We generalise this concept to higher order and clarify their relationship to other known statistics such as the skew-spectrum, the kurt-spectra and their real-space counterparts the cumulants correlators. Using the Hierarchical Ansatz (HA) as a toy model for the higher order correlation hierarchy, we show how in the soft limit, polyspectra at a given order can be identified with lower order polyspectra with the same geometrical dependence but with renormalised amplitudes expressed in terms of amplitudes of the original polyspectra. We extend the concept of position-dependent bispectrum to bispectrum of the divergence of the velocity field Θ and mixed multispectra involving δ and Θ in the 3D perturbative regime. To quantify the effects of transients in numerical simulations, we also present results for lowest order in Lagrangian perturbation theory (LPT) or the Zel'dovich approximation (ZA). Finally, we discuss how to extend the position-dependent spectrum concept to encompass cross-spectra. And finally study the application of this concept to two dimensions (2D), for projected galaxy maps, convergence κ maps from weak-lensing surveys or maps of CMB secondaries e.g. the frequency cleaned y-parameter maps of thermal Sunyaev-Zel'dovich (tSZ) effect from CMB surveys.

Gaussian streaming with the truncated Zel’dovich approximation

We calculate the dark matter halo correlation function in redshift space using the Gaussian streaming model (GSM). To determine the scale dependent functions entering the streaming model we use local Lagrangian bias together with Convolution Lagrangian perturbation theory (CLPT) which constitutes an approximation to the Post-Zel'dovich approximation. On the basis of N-body simulations we demonstrate that a smoothing of the initial conditions with the Lagrangian radius improves the Zel'dovich approximation and its ability to predict the displacement field of proto-halos. Based on this observation we implement a "truncated" CLPT by smoothing the initial power spectrum and investigate the dependence of the streaming model ingredients on the smoothing scale. We find that the real space correlation functions of halos and their mean pairwise velocity are optimised if the coarse graining scale is chosen to be 1 Mpc/h at z=0, while the pairwise velocity dispersion is optimised if the smoothing scale is chosen to be the Lagrangian size of the halo. We compare theoretical results for the halo correlation function in redshift space to measurements within the Horizon Run 2 N-body simulation halo catalog. We find that this simple two-filter smoothing procedure in the spirit of the truncated Zel'dovich approximation significantly improves the GSM+CLPT prediction of the redshift space halo correlation function over the whole mass range from large galaxy to galaxy cluster-sized halos.

The Gaussian streaming model and convolution Lagrangian effective field theory

We update the ingredients of the Gaussian streaming model (GSM) for the redshift-space clustering of biased tracers using the techniques of Lagrangian perturbation theory, effective field theory (EFT) and a generalized Lagrangian bias expansion. After relating the GSM to the cumulant expansion, we present new results for the real-space correlation function, mean pairwise velocity and pairwise velocity dispersion including counter terms from EFT and bias terms through third order in the linear density, its leading derivatives and its shear up to second order. We discuss the connection to the Gaussian peaks formalism. We compare the ingredients of the GSM to a suite of large N-body simulations, and show the performance of the theory on the low order multipoles of the redshift-space correlation function and power spectrum. We highlight the importance of a general biasing scheme, which we find to be as important as higher-order corrections due to non-linear evolution for the halos we consider on the scales of interest to us.

Improving fast generation of halo catalogues with higher order Lagrangian perturbation theory

We present the latest version of Pinocchio, a code that generates catalogues of DM haloes in an approximate but fast way with respect to an N-body simulation. This code version extends the computation of particle and halo displacements up to 3rd-order Lagrangian Perturbation Theory (LPT), in contrast with previous versions that used Zeldovich approximation (ZA). We run Pinocchio on the same initial configuration of a reference N-body simulation, so that the comparison extends to the object-by-object level. We consider haloes at redshifts 0 and 1, using different LPT orders either for halo construction - where displacements are needed to decide particle accretion onto a halo or halo merging - or to compute halo final positions. We compare the clustering properties of Pinocchio haloes with those from the simulation by computing the power spectrum and 2-point correlation function (2PCF) in real and redshift space (monopole and quadrupole), the bispectrum and the phase difference of halo distributions. We find that 2LPT and 3LPT give noticeable improvement. 3LPT provides the best agreement with N-body when it is used to displace haloes, while 2LPT gives better results for constructing haloes. At the highest orders, linear bias is typically recovered at a few per cent level. In Fourier space and using 3LPT for halo displacements, the halo power spectrum is recovered to within 10 per cent up to $k_{max}\sim0.5\ h/$Mpc. The results presented in this paper have interesting implications for the generation of large ensemble of mock surveys aimed at accurately compute covariance matrices for clustering statistics.

Spherical collapse of dark matter haloes in tidal gravitational fields

We study the spherical collapse model in the presence of external gravitational tidal shear fields for different dark energy scenarios and investigate the impact on the mass function and cluster number counts. While previous studies of the influence of shear and rotation on $\delta_\mathrm{c}$ have been performed with heuristically motivated models, we try to avoid this model dependence and sample the external tidal shear values directly from the statistics of the underlying linearly evolved density field based on first order Lagrangian perturbation theory. Within this self-consistent approach, in the sense that we restrict our treatment to scales where linear theory is still applicable, only fluctuations larger than the scale of the considered objects are included into the sampling process which naturally introduces a mass dependence of $\delta_\mathrm{c}$. We find that shear effects are predominant for smaller objects and at lower redshifts, i. e. the effect on $\delta_\mathrm{c}$ is at or below the percent level for the $\Lambda$CDM model. For dark energy models we also find small but noticeable differences, similar to $\Lambda$CDM. The virial overdensity $\Delta_\mathrm{V}$ is nearly unaffected by the external shear. The now mass dependent $\delta_c$ is used to evaluate the mass function for different dark energy scenarios and afterwards to predict cluster number counts, which indicate that ignoring the shear contribution can lead to biases of the order of $1\sigma$ in the estimation of cosmological parameters like $\Omega_\mathrm{m}$, $\sigma_8$ or $w$.

Modeling the reconstructed BAO in Fourier space

The density field reconstruction technique, which was developed to partially reverse the nonlinear degradation of the Baryon Acoustic Oscillation (BAO) feature in the galaxy redshift surveys, has been successful in substantially improving the cosmology constraints from recent galaxy surveys such as Baryon Oscillation Spectroscopic Survey (BOSS). We estimate the efficiency of the reconstruction method as a function of various reconstruction details. To directly quantify the BAO information in nonlinear density fields before and after reconstruction, we calculate the cross-correlations (i.e., propagators) of the pre(post)-reconstructed density field with the initial linear field using a mock galaxy sample that is designed to mimic the clustering of the BOSS CMASS galaxies. The results directly provide the BAO damping as a function of wavenumber that can be implemented into the Fisher matrix analysis. We focus on investigating the dependence of the propagator on a choice of smoothing filters and on two major different conventions of the redshift-space density field reconstruction that have been used in literature. By estimating the BAO signal-to-noise for each case, we predict constraints on the angular diameter distance and Hubble parameter using the Fisher matrix analysis. We thus determine an optimal Gaussian smoothing filter scale for the signal-to-noise level of the BOSS CMASS. We also present appropriate BAO fitting models for different reconstruction methods based on the first and second order Lagrangian perturbation theory in Fourier space. Using the mock data, we show that the modified BAO fitting model can substantially improve the accuracy of the BAO position in the best fits as well as the goodness of the fits.