Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userRapid Separable Analysis of Higher Order Correlators in Large Scale Structure
194
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present an efficient separable approach to the estimation and reconstruction of the bispectrum and the trispectrum from observational (or simulated) large scale structure data. This is developed from general CMB (poly-)spectra methods which exploit the fact that the bispectrum and trispectrum in the literature can be represented by a separable mode expansion which converges rapidly (with $n_textrm{max}={cal{O}}(30)$ terms). With an effective grid resolution $l_textrm{max}$ (number of particles/grid points $N=l_textrm{max}^3$), we present a bispectrum estimator which requires only ${cal O}(n_textrm{max} times l_textrm{max}^3)$ operations, along with a corresponding method for direct bispectrum reconstruction. This method is extended to the trispectrum revealing an estimator which requires only ${cal O}(n_textrm{max}^{4/3} times l_textrm{max}^3)$ operations. The complexity in calculating the trispectrum in this method is now involved in the original decomposition and orthogonalisation process which need only be performed once for each model. However, for non-diagonal trispectra these processes present little extra difficulty and may be performed in ${cal O}(l_textrm{max}^4)$ operations. A discussion of how the methodology may be applied to the quadspectrum is also given. An efficient algorithm for the generation of arbitrary nonGaussian initial conditions for use in N-body codes using this separable approach is described. This prescription allows for the production of nonGaussian initial conditions for arbitrary bispectra and trispectra. A brief outline of the key issues involved in parameter estimation, particularly in the non-linear regime, is also given.
rate research
Read More
An important aspect of large-scale structure data analysis is the presence of non-negligible theoretical uncertainties, which become increasingly important on small scales. We show how to incorporate these uncertainties in realistic power spectrum likelihoods by an appropriate change of the fitting model and the covariance matrix. The inclusion of the theoretical error has several advantages over the standard practice of using the sharp momentum cut $k_{rm max}$. First, the theoretical error covariance gradually suppresses the information from the short scales as the employed theoretical model becomes less reliable. This allows one to avoid laborious measurements of $k_{rm max}$, which is an essential part of the standard methods. Second, the theoretical error likelihood gives unbiased constrains with reliable error bars that are not artificially shrunk due to over-fitting. In realistic settings, the theoretical error likelihood yields essentially the same parameter constraints as the standard analysis with an appropriately selected $k_{rm max}$, thereby effectively optimizing the choice of $k_{rm max}$. We demonstrate these points using the large-volume N-body data for the clustering of matter and galaxies in real and redshift space. In passing, we validate the effective field theory description of the redshift space distortions and show that the use of the one-parameter phenomenological Gaussian damping model for fingers-of-God causes significant biases in parameter recovery.
Cosmological perturbations of sufficiently long wavelength admit a fluid dynamic description. We consider modes with wavevectors below a scale $k_m$ for which the dynamics is only mildly non-linear. The leading effect of modes above that scale can be accounted for by effective non-equilibrium viscosity and pressure terms. For mildly non-linear scales, these mainly arise from momentum transport within the ideal and cold but inhomogeneous fluid, while momentum transport due to more microscopic degrees of freedom is suppressed. As a consequence, concrete expressions with no free parameters, except the matching scale $k_m$, can be derived from matching evolution equations to standard cosmological perturbation theory. Two-loop calculations of the matter power spectrum in the viscous theory lead to excellent agreement with $N$-body simulations up to scales $k=0.2 , h/$Mpc. The convergence properties in the ultraviolet are better than for standard perturbation theory and the results are robust with respect to variations of the matching scale.
We study soft limits of correlation functions for the density and velocity fields in the theory of structure formation. First, we re-derive the (resummed) consistency conditions at unequal times using the eikonal approximation. These are solely based on symmetry arguments and are therefore universal. Then, we explore the existence of equal-time relations in the soft limit which, on the other hand, depend on the interplay between soft and hard modes. We scrutinize two approaches in the literature: The time-flow formalism, and a background method where the soft mode is absorbed into a locally curved cosmology. The latter has been recently used to set up (angular averaged) `equal-time consistency relations. We explicitly demonstrate that the time-flow relations and `equal-time consistency conditions are only fulfilled at the linear level, and fail at next-to-leading order for an Einstein de-Sitter universe. While applied to the velocities both proposals break down beyond leading order, we find that the `equal-time consistency conditions quantitatively approximates the perturbative results for the density contrast. Thus, we generalize the background method to properly incorporate the effect of curvature in the density and velocity fluctuations on short scales, and discuss the reasons behind this discrepancy. We conclude with a few comments on practical implementations and future directions.
We present a case study describing efforts to optimise and modernise Modal, the simulation and analysis pipeline used by the Planck satellite experiment for constraining general non-Gaussian models of the early universe via the bispectrum (or three-point correlator) of the cosmic microwave background radiation. We focus on one particular element of the code: the projection of bispectra from the end of inflation to the spherical shell at decoupling, which defines the CMB we observe today. This code involves a three-dimensional inner product between two functions, one of which requires an integral, on a non-rectangular domain containing a sparse grid. We show that by employing separable methods this calculation can be reduced to a one-dimensional summation plus two integrations, reducing the overall dimensionality from four to three. The introduction of separable functions also solves the issue of the non-rectangular sparse grid. This separable method can become unstable in certain cases and so the slower non-separable integral must be calculated instead. We present a discussion of the optimisation of both approaches. We show significant speed-ups of ~100x, arising from a combination of algorithmic improvements and architecture-aware optimisations targeted at improving thread and vectorisation behaviour. The resulting MPI/OpenMP hybrid code is capable of executing on clusters containing processors and/or coprocessors, with strong-scaling efficiency of 98.6% on up to 16 nodes. We find that a single coprocessor outperforms two processor sockets by a factor of 1.3x and that running the same code across a combination of both microarchitectures improves performance-per-node by a factor of 3.38x. By making bispectrum calculations competitive with those for the power spectrum (or two-point correlator) we are now able to consider joint analysis for cosmological science exploitation of new data.
Power suppression of the cosmic microwave background on the largest observable scales could provide valuable clues about the particle physics underlying inflation. Here we consider the prospect of power suppression in the context of the multifield landscape. Based on the assumption that our observable universe emerges from a tunnelling event and that the relevant features originate purely from inflationary dynamics, we find that the power spectrum not only contains information on single-field dynamics, but also places strong con- straints on all scalar fields present in the theory. We find that the simplest single-field models giving rise to power suppression do not generalise to multifield models in a straightforward way, as the resulting superhorizon evolution of the curvature perturbation tends to erase any power suppression present at horizon crossing. On the other hand, multifield effects do present a means of generating power suppression which to our knowledge has so far not been considered. We propose a mechanism to illustrate this, which we dub flume inflation.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
