This paper introduces the Wavelet Phase Harmonics (WPH) statistics. They are interpretable low-dimensional statistics which describe 2D non-Gaussian density fields. These statistics are built from WPH moments, which have been recently introduced in data science and machine learning community. In this paper, we applied WPH statistics to projected matter density fields from the Quijote N-body simulations. We find by computing the Fisher information matrix, that the WPH statistics can place more stringent constraints on 5 cosmological parameters when compared to the combination of power-spectrum and bi-spectrum. We also use the WPH statistics to successfully generate from a maximum entropy model new 2D density fields that reproduce the PDF, mean, power-spectrum, bispectrum and the Minkowski functionals of the input density fields. While separate methods have been proven very efficient for parameter estimations and statistical syntheses for for large scale structure, WPH statistics are the first statistics that can both achieve a more stringent cosmological parameter constraint, and produce a sufficiently accurate simulation of the Universe while being interpretable.

The interstellar medium (ISM) is a complex non-linear system governed by gravity and magneto-hydrodynamics, as well as radiative, thermodynamical, and chemical processes. Our understanding of it mostly progresses through observations and numerical simulations, and a quantitative comparison between these two approaches requires a generic and comprehensive statistical description. The goal of this paper is to build such a description, with the purpose to permit an efficient comparison independent of any specific prior or model. We start from the Wavelet Scattering Transform (WST), a low-variance statistical description of non-Gaussian processes, developed in data science, that encodes long-range interactions through a hierarchical multiscale approach based on the Wavelet transform. We perform a reduction of the WST through a fit of its angular dependencies, allowing to gather most of the information it contains into a few components whose physical meanings are identified, and that describe, e.g., isotropic and anisotropic behaviours. The result of this paper is the Reduced Wavelet Scattering Transform (RWST), a statistical description with a small number of coefficients that characterizes complex structures arising from non-linear phenomena, free from any specific prior. The RWST coefficients encode moments of order up to four, have reduced variances, and quantify the couplings between scales. To show the efficiency and generality of this description, we apply it successfully to three kinds of processes: fractional Brownian motions, MHD simulations, and Herschel observations in a molecular cloud. With fewer than 100 coefficients when probing 6 scales and 8 angles on 256*256 maps, we were able with the RWST to perform quantitative comparisons, to infer relevant physical properties, and to produce realistic synthetic fields.

The Quijote simulations are a set of 43100 full N-body simulations spanning more than 7000 cosmological models in the $\{\Omega_{\rm m}, \Omega_{\rm b}, h, n_s, \sigma_8, M_\nu, w \}$ hyperplane. At a single redshift the simulations contain more than 8.5 trillions of particles over a combined volume of 43100 $(h^{-1}{\rm Gpc})^3$. Billions of dark matter halos and cosmic voids have been identified in the simulations, whose runs required more than 35 million core hours. The Quijote simulations have been designed for two main purposes: 1) to quantify the information content on cosmological observables, and 2) to provide enough data to train machine learning algorithms. In this paper we describe the simulations and show a few of their applications. We also release the Petabyte of data generated, comprising hundreds of thousands of simulation snapshots at multiple redshifts, halo and void catalogs, together with millions of summary statistics such as power spectra, bispectra, correlation functions, marked power spectra, and estimated probability density functions.

The current description of fundamental interactions is based on two theories with the status of standard models. The electromagnetic and nuclear interactions are described at a quantum level by the Standard Model of particle physics, using tools like gauge theories and spontaneous symmetry breaking by the Higgs mechanism. The gravitational interaction is described on the other hand by general relativity, based on a dynamical description of space-time at a classical level. Although these models are verified to high precision in the solar system experiments, they suffer from several theoretical weaknesses and a lack of predictive power at the Planck scale as well as at cosmological scales; they are thus not viewed anymore as fundamental theories. As its phenomenology involves both these extreme scales, cosmology is then a good laboratory to probe theories going beyond these standard models. The first part of this thesis focus on cosmic strings, topological defects forming during the spontaneous symmetry breaking of grand unified theories in the early universe. I show especially how to study these defects while taking into account the complete structure of the particles physics models leading to their formation, going beyond the standard descriptions in terms of simplified toy-models. The second part is devoted to the construction and the examination of different theories of modified gravity related to the Galileon model, a model which tries in particular to explain the dark energy phenomenology.

Following previous works on generalized Abelian Proca theory, also called vector Galileon, we investigate the massive extension of an SU(2) gauge theory, i.e., the generalized SU(2) Proca model, which could be dubbed non-Abelian vector Galileon. This particular symmetry group permits fruitful applications in cosmology such as inflation driven by gauge fields. Our approach consists in building, in an exhaustive way, all the Lagrangians containing up to six contracted Lorentz indices. For this purpose, and after identifying by group theoretical considerations all the independent Lagrangians which can be written at these orders, we consider the only linear combinations propagating three degrees of freedom and having healthy dynamics for their longitudinal mode, i.e., whose pure St\"uckelberg contribution turns into the SU(2) multi-Galileon dynamics. Finally, and after having considered the curved space-time expansion of these Lagrangians, we discuss the form of the theory at all subsequent orders.

We summarize previous results on the most general Proca theory in 4 dimensions containing only first-order derivatives in the vector field (second-order at most in the associated St\"uckelberg scalar) and having only three propagating degrees of freedom with dynamics controlled by second-order equations of motion. Discussing the Hessian condition used in previous works, we conjecture that, as in the scalar galileon case, the most complete action contains only a finite number of terms with second-order derivatives of the St\"uckelberg field describing the longitudinal mode, which is in agreement with the results of JCAP 1405, 015 (2014) and Phys. Lett. B 757, 405 (2016) and complements those of JCAP 1602, 004 (2016). We also correct and complete the parity violating sector, obtaining an extra term on top of the arbitrary function of the field $A_\mu$, the Faraday tensor $F_{\mu \nu}$ and its Hodge dual $\tilde{F}_{\mu \nu}$.

We revisit the most general theory for a massive vector field with derivative self-interactions, extending previous works on the subject to account for terms having trivial total derivative interactions for the longitudinal mode. In the flat spacetime (Minkowski) case, we obtain all the possible terms containing products of up to five first-order derivatives of the vector field, and provide a conjecture about higher-order terms. Rendering the metric dynamical, we covariantize the results and add all possible terms implying curvature.