We study the effect of primordial scalar curvature perturbations on the propagation of gravitational waves over cosmic distances. We point out that such curvature perturbations deform the isotropic spectrum of any stochastic background of gravitational waves of primordial origin through the (integrated) Sachs-Wolfe effect. Computing the changes in the amplitude and frequency of the propagating gravitational wave induced at linear order by scalar curvature perturbations, we show that the resulting deformation of each frequency bin of the gravitational wave spectrum is described by a linearly biased Gaussian with the variance $\sigma^2 \simeq \int d\ln k \Delta_{\mathcal R}^2$, where $\Delta_{\mathcal R}^2(k)$ denotes the amplitude of the primordial curvature perturbations. The linear bias encodes the correlations between the changes induced in the frequency and amplitude of the gravitational waves. Taking into account the latest bounds on $\Delta_{\mathcal R}^2$ from primordial black hole and gravitational wave searches, we demonstrate that the resulting ${\mathcal O}(\sigma)$ deformation can be significant for extremely peaked gravitational wave spectra. We further provide an order of magnitude estimate for broad spectra, for which the net distortion is ${\mathcal O}(\sigma^2)$.

We study the Schr\"odinger-Poisson (SP) method in the context of cosmological large-scale structure formation in an expanding background. In the limit $\hbar \to 0$, the SP technique can be viewed as an effective method to sample the phase space distribution of cold dark matter that remains valid on non-linear scales. We present results for the 2D and 3D matter correlation function and power spectrum at length scales corresponding to the baryon acoustic oscillation (BAO) peak. We discuss systematic effects of the SP method applied to cold dark matter and explore how they depend on the simulation parameters. In particular, we identify a combination of simulation parameters that controls the scale-independent loss of power observed at low redshifts, and discuss the scale relevant to this effect.

Within the Halo Model of large scale structure, all matter is contained in dark matter halos. This simple yet powerful framework has been broadly applied to multiple data sets and enriched our comprehension of how matter is distributed in the Universe. In this work we extend this assumption by allowing for matter to rest not only inside halos but also within cosmic voids and in between halos and voids (which we call 'dust'). This assumption leads to additional contributions (1Void, 2Void, Halo-Void, etc.) to the predictions of correlation functions, spectra and profiles for both halos and voids. Whereas the Halo Model can only make predictions for halo quantities, the Halo Void Model extends those for void statistics and halo-void cross-correlations. We provide recipes for all new ingredients of the Halo Void (Dust) Model, such as the void abundance, linear bias and density profile and test their validity in a N-body simulation. Including voids and dust into the calculations improves the transition between the 1Halo and the 2Halo terms by up to $\sim 6\%$. It also eliminates the need to include low-mass structures on the normalization of large-scale terms, suggesting that halos and voids are complementary cosmic structures to effectively describe matter distribution on large scales of the Universe.

We study the power spectrum of dark matter density fluctuations in the framework of the Effective Field Theory of Large Scale Structures (EFTofLSS) up to three loop orders. In principle, several counter-terms may be needed to handle the short-distance sensitivity in perturbation theory. However, we show that a small number of extra coefficients are sufficient to match numerical simulations with percent accuracy when a generic renormalization prescription is implemented (allowing for running of the individual counter-terms). We show that the level of accuracy increases with respect to the two loop results, up to $k \simeq 0.4\, h$Mpc$^{-1}$ at redshift $z=0$, although the overall improvement is somewhat marginal. At the same time, we argue there is evidence that the behavior of the loop expansion in the EFTofLSS is typical of an asymptotic series, already on the brink of its maximum predictive power (at $z=0$). Hence, the inclusion of higher orders will likely deteriorate the matching to data, even at moderate values of $k$. Part of the reason for this behavior is due to large contributions to the (renormalized) power spectrum at three loop order from mildly non-linear scales, even after the UV counter-terms are included. In conclusion, the EFTofLSS to three loop orders provides the best approximation to the (deterministic part of the) power spectrum in the weakly non-linear regime at $z=0$, and higher loops are not expected to improve our level of accuracy.