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 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.

The LISA telescope will provide the first opportunity to probe the scenario of a first-order phase transition happening close to the electroweak scale. By now, it is evident that the main contribution to the GW spectrum comes from the sound waves propagating through the plasma. Current estimates of the GW spectrum are based on numerical simulations of a scalar field interacting with the plasma or on analytical approximations -the so-called sound shell model. In this work we present a novel setup to calculate the GW spectra from sound waves. We use a hybrid method that uses a 1d simulation (with spherical symmetry) to evolve the velocity and enthalpy profiles of a single bubble after collision and embed it in a 3d realization of multiple bubble collisions, assuming linear superposition of the velocity and enthalpy. The main advantage of our method compared to 3d hydrodynamic simulations is that it does not require to resolve the scale of bubble wall thickness. This makes our simulations more economical and the only two relevant physical length scales that enter are the bubble size and the shell thickness (that are in turn enclosed by the box size and the grid spacing). The reduced costs allow for extensive parameter studies and we provide a parametrization of the final GW spectrum as a function of the wall velocity and the fluid kinetic energy.

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 ∼ 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 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 σ 2 d lnk ∆ 2 R , where ∆ 2 R (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 ∆ 2 R from primordial black hole and gravitational wave searches, we demonstrate that the resulting O(σ) 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 O(σ 2 ).