Within cosmological perturbation theory, the cosmic microwave background anisotropies are usually computed from a Boltzmann hierarchy coupled to the perturbed Einstein equations. In this setup, one set of multipoles describes the temperature anisotropies, while two other sets, of electric and magnetic types, describe the polarization anisotropies. In order to reduce the number of multipoles types needed for polarization, and thus to speed up the numerical resolution, an optimal hierarchy has been proposed in the literature for Einstein-Boltzmann codes. However, it has been recently shown that the separability between directional and orbital eigenfunctions employed in the optimal hierarchy is not correct in the presence of spatial curvature. We investigate how the assumption of separability affects the optimal hierarchy, and show that it introduces relative errors of order $\Omega_K$ with respect to the full hierarchy. Despite of that, we show that the optimal hierarchy still gives extremely good results for temperature and polarization angular spectra, with relative errors that are much smaller than cosmic variance even for curvatures as large as $|\Omega_K|=0.1$. Still, we find that the polarization angular spectra from tensor perturbations are significantly altered when using the optimal hierarchy, leading to errors that are typically of order $50 |\Omega_K| \%$ on that component.

In the literature different approaches have been proposed to compute the anisotropies of the astrophysical gravitational wave background. The different expressions derived, although starting from our work Cusin, Pitrou, Uzan, Phys.Rev.D96, 103019 (2017) [1], seem to differ. This article compares the various theoretical expressions proposed so far and provides a separate derivation based on a Boltzmann approach. We show that all the theoretical formula in the literature are equivalent and boil down to the one of Ref. [1] when a proper matching of terms and integration by parts are performed. The difference between the various predictions presented for anisotropies in a cosmological context can only lie in the astrophysical modeling of sources, and neither in the theory nor in the cosmological description of the large scale structures. Finally we comment on the gauge invariance of expressions.

In the primordial Universe, neutrino decoupling occurs only slightly before electron-positron annihilations, leading to an increased neutrino energy density with order $10^{-2}$ spectral distortions compared to the standard instantaneous decoupling approximation. However, there are discrepancies in the literature on the impact it has on the subsequent primordial nucleosynthesis, in terms of both the magnitude of the abundance modifications and their sign. We review how neutrino decoupling indirectly affects the various stages of nucleosynthesis, namely, the freezing out of neutron abundance, the duration of neutron beta decay, and nucleosynthesis itself. This allows to predict the sign of the abundance variations that are expected when the physics of neutrino decoupling is taken into account. For simplicity, we ignore neutrino oscillations, but we conjecture from the detailed interplay of neutrino temperature shifts and distortions that their effect on final light element abundances should be subdominant.

The spatial and temporal discreteness of gravitational wave sources leads to shot noise that may, in some regimes, swamp any attempts at measuring the anisotropy of the gravitational wave background. Cross-correlating a gravitational wave background map with a sufficiently dense galaxy survey can alleviate this issue, and potentially recover some of the underlying properties of the gravitational wave background. We quantify the shot noise level and we explicitly show that cross-correlating the gravitational wave background and a galaxy catalog improves the chances of a first detection of the background anisotropy with a gravitational wave observatory operating in the frequency range (10Hz,100Hz), given sufficient sensitivity.

We show that the anisotropies of the astrophysical stochastic gravitational wave background in the mHz band have a strong dependence on the modelling of galactic and sub-galactic physics. We explore a wide range of self-consistent astrophysical models for stellar evolution and for the distribution of orbital parameters, all calibrated such that they predict the same number of resolved mergers to fit the number of detections during LIGO/Virgo O1+O2 observations runs. We show that different physical choices for the process of black hole collapse and cut-off in the black hole mass distribution give fractional differences in the angular power spectrum of anisotropies up to 50\% on all angular scales. We also point out that the astrophysical information which can be extracted from anisotropies is complementary to the isotropic background and individual mergers. These results underline the interest in the anisotropies of the stochastic gravitational wave background as a new and potentially rich field of research, at the cross-road between astrophysics and cosmology.

A powerful result in theoretical cosmology states that a subset of anisotropic Bianchi models can be seen as the homogeneous limit of (standard) linear cosmological perturbations. Such models are precisely those leading to Friedmann spacetimes in the limit of zero anisotropy. Building on previous works, we give a comprehensive exposition of this result, and perform the detailed identification between anisotropic degrees of freedom and their corresponding scalar, vector, and tensor perturbations of standard perturbation theory. In particular, we find that anisotropic models very close to open (i.e., negatively curved) Friedmann spaces correspond to some type of super-curvature perturbations. As a consequence, provided anisotropy is mild, its effects on all types of cosmological observables can always be computed as simple extensions of the standard techniques used in relativistic perturbation theory around Friedmann models. This fact opens the possibility to consistently constrain, for all cosmological observables, the presence of large scale anisotropies on the top of the stochastic fluctuations.

Precision on primordial abundances, deduced from observations, have now reached the percent level for 4He and deuterium. Precision on big bang nucleosynthesis (BBN) predictions should, hence, reach the same level. The uncertainty on the 4He mass fraction is strongly affected by theoretical uncertainties on the weak reaction rates that interconvert neutrons with protons. All these corrections have been calculated in a self-consistent manner and implemented in a new, and public, Mathematica code PRIMAT, together with an extensive data base of reaction rates. Both can be obtained at http://www2.iap.fr/users/pitrou/primat.htm.

This article explores the properties (amplitude and shape) of the angular power spectrum of the anisotropies of the astrophysical gravitational wave background (AGWB) focusing on the signatures of the astrophysical models describing sub-galactic physics. It demonstrates that while some parameters have negligible impact others, and in particular the stellar evolution models, the metallicity and the merger time delay distribution can result in relative differences of order 40% in the angular power spectrum of anisotropies in both the LIGO/Virgo and LISA frequency bands. It is also shown that the monopole and the anisotropic components of the AGWB are complementary and sensitive to different astrophysical parameters. It follows that AGWB anisotropies are a new observable with the potential to provide new astrophysical information that can not be accessed otherwise.

We compute radiative corrections to nuclear reaction rates that determine the outcome of the Big-Bang Nucleosynthesis (BBN). Any nuclear reaction producing a photon with an energy above $2m_e$ must be supplemented by the corresponding reaction where the final state photon is replaced by an electron-positron pair. We find that pair prodution brings a typical $0.2 \%$ enhancement to photon emission rates, resulting in a similar size corrections to elemental abundances. The exception is $^4{\rm He}$ abundance, which is completely insenstive to the small changes in the nuclear reaction rates. We also investigate the effect of vacuum polarisation on the Coulomb barrier, which brings a small extra correction when reaction rates are extrapolated from the measured energies to the BBN Gamow peak energies.

We review the general construction of distribution functions for gases of fermions and bosons (photons), emphasizing the similarities and differences between both cases. The central object which describes polarization for photons is a tensor-valued distribution function, whereas for fermions it is a vector-valued one. The collision terms of Boltzmann equations for fermions and bosons also possess the same general structure and differ only in the quantum effects associated with the final state of the reactions described. In particular, neutron-proton conversions in the early universe, which set the primordial Helium abundance, enjoy many similarities with Compton scattering which shapes the cosmic microwave background and we show that both can be handled with a Fokker-Planck type expansion. For neutron-proton conversions, this allows to obtain the finite nucleon mass corrections, required for precise theoretical predictions, whereas for Compton scattering it leads to the thermal and recoil effects which enter the Kompaneets equation. We generalize the latter to the general case of anisotropic and polarized photon distribution functions. Finally we discuss a parameterization of the photon spectrum based on logarithmic moments which allows for a neat separation between temperature shifts and spectral distortions.

We investigate the discrepancy pointed out by Jenkins et al. in Ref. [1] between the predictions of anisotropies of the astrophysical gravitational wave (GW) background, derived using different methods in Cusin et al. [2] and in Jenkins et al. [3]. We show that this discrepancy is not due to our treatment of galaxy clustering, contrary to the claim made in Ref. [1] and we show that our modeling of clustering gives results in very good agreement with observations. Furthermore we show that the power law spectrum used in Refs. [1] and [3] to describe galaxy clustering is incorrect on large scales and leads to a different scaling of the multipoles $C_\ell$. Moreover, we also explain that the analytic derivation of the gravitational wave background correlation function in Refs. [1] and [3] is mathematically ill-defined and predicts an amplitude of the angular power spectrum which depends on the (arbitrary) choice of a non-physical cut-off.

Unresolved and resolved sources of gravitational waves are at the origin of a stochastic gravitational wave background. While the computation of its mean density as a function of frequency in a homogeneous and isotropic universe is standard lore, the computation of its anisotropies requires to understand the coarse graining from local systems, to galactic scales and then to cosmology. An expression of the gravitational wave energy density valid in any general spacetime is derived. It is then specialized to a perturbed Friedmann-Lema\^itre spacetime in order to determine the angular power spectrum of this stochastic background as well as its correlation with other cosmological probes, such as the galaxy number counts and weak lensing. Our result for the angular power spectrum also provides an expression for the variance of the gravitational wave background.

Primordial nucleosynthesis is one of the three historical evidences for the big bang model, together with the expansion of the universe and the cosmic microwave background. Now that the number of neutrino families and the baryonic densities have been fixed by laboratory measurements or CMB observations, the model has no free parameter and its predictions are rigid. Departure from its predictions could provide hints or constraints on new physics or astrophysics in the early universe. Precision on primordial abundances deduced from observations have recently been drastically improved and reach the percent level for both deuterium and helium-4. Accordingly, the BBN predictions should reach the same level of precision. For most isotopes, the dominant sources of uncertainty come from those on the laboratory thermonuclear reactions. This article focuses on helium-4 whose predicted primordial abundance depends essentially on weak interactions which control the neutron-proton ratio. The rates of the various weak interaction processes depend on the experimentally measured neutron lifetime, but also includes numerous corrections that we thoroughly investigate here. They are the radiative, zero-temperature, corrections, finite nucleon mass corrections, finite temperature radiative corrections, weak-magnetism, and QED plasma effects, which are for the first time all included and calculated in a self consistent way, allowing to take into account the correlations between them, and verifying that all satisfy detailed balance. The helium-4 predicted mass fraction is $0.24709\pm0.00017$. In addition, we provide a Mathematica code (PRIMAT) that incorporates, not only these corrections but also a full network of reactions, using the best available thermonuclear reaction rates, allowing the predictions of primordial abundances up to the CNO region.

This article investigates the full Boltzmann equation up to second order in the cosmological perturbations. Describing the distribution of polarized radiation by a tensor valued distribution function, we study the gauge dependence of the distribution function and summarize the construction of the gauge-invariant distribution function. The Liouville operator which describes the free streaming of electrons, and the collision term which describes the scattering of photons on free electrons are computed up to second order. Finally, the remaining dependence in the direction of the photon momentum is handled by expanding in projected symmetric trace-free multipoles and also in the more commonly used normal modes components. The results obtained remain to be used for computing numerically the contribution in the cosmic microwave background bispectrum which arises from the evolution of second order perturbations, in order to disentangle the primordial non-Gaussianity from the one generated by the subsequent non-linear evolution.

High precision astrometry now enables to measure the time drift of astrophysical observables in real time, hence providing new ways to probe different cosmological models. This article presents a general derivation of the redshift and direction drifts for general observers. It is then applied to the standard cosmological framework of a Friedmann-Lemaitre space- time including all effects at first order in the cosmological perturbations, as well as in the class of spatially anisotropic universe models of the Bianchi I family. It shows that for a general observer, the direction drift splits into a parallax and an aberration drifts and order of magnitude estimates of these two components are provided. The multipolar decomposition of the redshift and aberration drifts is also derived and shows that the observer's peculiar velocity contributes only as a dipole whereas the anisotropic shear contributes as a quadrupole.

We present the first predictions for the angular power spectrum of the astrophysical gravitational wave background constituted of the radiation emitted by all resolved and unresolved astrophysical sources. Its shape and amplitude depend on both the astrophysical properties on galactic scales and on cosmological properties. We show that the angular power spectrum behaves as $C_{\ell}\propto 1/{\ell}$ on large scales and that relative fluctuations of the signal are of order 30% at 100 Hz. We also present the correlations of the astrophysical gravitational wave background with weak-lensing and galaxy distribution. These numerical results pave the way to the study of a new observable at the crossroad between general relativity, astrophysics and cosmology.

The gravitational wave radiation emitted by all, resolved and unresolved, astrophysical sources in the observable universe generates a stochastic background. This background has a directional dependence inherited from the inhomogeneities of the matter distribution. This article proposes a new and independent derivation of the angular dependence of its energy density by focusing on the total gravitational wave signal produced by an ensemble of incoherent sources. This approach clarifies the origin of the angular correlation and the relation between the gravitational wave signal that can be measured by interferometers and the energy density of the stochastic background.

We build a statistical description of fermions, taking into account the spin degree of freedom in addition to the momentum of particles, and we detail its use in the context of the kinetic theory of gases of fermions particles. We show that the one-particle distribution function needed to write a Liouville equation is a spinor valued operator. The degrees of freedom of this function are covariantly described by an intensity function and by a polarisation vector which are parallel transported by free streaming. Collisions are described on the microscopic level and lead to a Boltzmann equation for this operator. We apply our formalism to the case of weak interactions, which at low energies can be considered as a contact interaction between fermions, allowing us to discuss the structure of the collision term for a few typical weak-interaction mediated reactions. In particular we find for massive particles that a dipolar distribution of velocities in the interacting species is necessary to generate linear polarisation, as opposed to the case of photons for which linear polarisation is generated from the quadrupolar distribution of velocities.

The properties of the cosmic microwave background (CMB) temperature and polarisation anisotropies measured by a static, off-centered observer located in a local spherically symmetric void, are described. In particular in this paper we compute, together with the standard 2- point angular correlation functions, the off-diagonal correlators, which are no more vanishing by symmetry. While the energy shift induced by the off-centered position of the observer can be suppressed by a proper choice of the observer velocity, a lensing-like effect on the CMB emission point remains. This latter effect is genuinely geometrical (e.g. non-degenerate with a boost) and reflects in the structure of the off-diagonal correlators. At lowest order in this effect, the temperature and polarisation correlation matrices have non-vanishing diagonal elements, as usual, and all the off-diagonal terms are excited. This particular signature of a local void model allows one, in principle, to disentangle geometrical effects from local kinematical ones in CMB observations.

This article details the computation of the two-point correlators of the convergence, $E$- and $B$-modes of the cosmic shear induced by the weak-lensing by large scale structure assuming that the background spacetime is spatially homogeneous and anisotropic. After detailing the perturbation equations and the general theory of weak-lensing in an anisotropic universe, it develops a weak shear approximation scheme in which one can compute analytically the evolution of the Jacobi matrix. It allows one to compute the angular power spectrum of the $E$- and $B$-modes. In the linear regime, the existence of $B$-modes is a direct tracer of a late time anisotropy and their angular power spectrum scales as the square of the shear. It is then demonstrated that there must also exist off-diagonal correlations between the $E$-modes, $B$-modes and convergence that are linear in the geometrical shear and allow one to reconstruct the eigendirections of expansion. These spectra can be measured in future large scale surveys, such as Euclid and SKA, and offer a new tool to test the isotropy of the expansion of the universe at low redshift.

We compute the angular power spectrum of the $B$-modes of the weak-lensing shear in a spatially anisotropic spacetime. We find that there must also exist off-diagonal correlations between the $E$-modes, $B$-modes, and convergence that allow one to reconstruct the eigendirections of expansion. Focusing on future surveys such as Euclid and SKA, we show that observations can constrain the geometrical shear in units of the Hubble rate at the percent level, or even better, offering a new and powerful method to probe our cosmological model.

Redshift-space distortions are generally considered in the plane parallel limit, where the angular separation between the two sources can be neglected. Given that galaxy catalogues now cover large fractions of the sky, it becomes necessary to consider them in a formalism which takes into account the wide angle separations. In this article we derive an operational formula for the matter correlators in the Newtonian limit to be used in actual data sets, both in configuration and in Fourier spaces without relying on a plane-parallel approximation. We then recover the plane-parallel limit not only in configuration space where the geometry is simpler, but also in Fourier space, and we exhibit the first corrections that should be included in large surveys as a perturbative expansion over the plane-parallel results. We finally compare our results to existing literature, and show explicitly how they are related.

A primordial inflationary phase allows one to erase any possible anisotropic expansion thanks to the cosmic no-hair theorem. If there is no global anisotropic stress, then the anisotropic expansion rate tends to decrease. What are the observational consequences of a possible early anisotropic phase? We first review the dynamics of anisotropic universes and report analytic approximations. We then discuss the structure of dynamical equations for perturbations and the statistical properties of observables, as well as the implication of a primordial anisotropy on the quantization of these perturbations during inflation. Finally we briefly review models based on primordial vector field which evade the cosmic no-hair theorem.

We investigate the generation of magnetic fields from non-linear effects around recombination. As tight-coupling is gradually lost when approaching $z\simeq 1100$, the velocity difference between photons and baryons starts to increase, leading to an increasing Compton drag of the photons on the electrons. The protons are then forced to follow the electrons due to the electric field created by the charge displacement; the same field, following Maxwell's laws, eventually induces a magnetic field on cosmological scales. Since scalar perturbations do not generate any magnetic field as they are curl-free, one has to resort to second-order perturbation theory to compute the magnetic field generated by this effect. We reinvestigate this problem numerically using the powerful second-order Boltzmann code SONG. We show that: i) all previous studies do not have a high enough angular resolution to reach a precise and consistent estimation of the magnetic field spectrum; ii) the magnetic field is generated up to $z\simeq 10$; iii) it is in practice impossible to compute the magnetic field with a Boltzmann code for scales smaller than $1\,{\rm Mpc}$. Finally we confirm that for scales of a few ${\rm Mpc}$, this magnetic field is of order $2\times 10^{-29}{\rm G}$, many orders of magnitude smaller than what is currently observed on intergalactic scales.

This article proposes a comprehensive analysis of light propagation in an anisotropic and spatially homogeneous Bianchi I universe. After recalling that null geodesics are easily determined in such a spacetime, we derive the expressions of the redshift and direction drifts of light sources; by solving analytically the Sachs equation, we then obtain an explicit expression of the Jacobi matrix describing the propagation of narrow light beams. As a byproduct, we recover the old formula by Saunders for the angular diameter distance in a Bianchi I spacetime, but our derivation goes further since it also provides the optical shear and rotation. These results pave the way to the analysis of both supernovae data and weak lensing by the large-scale structure in Bianchi universes.

We reanalyze microlensing events in the published list of anomalous events that were observed from the OGLE lensing survey conducted during 2004-2008 period. In order to check the existence of possible degenerate solutions and extract extra information, we conduct analyses based on combined data from other survey and follow-up observation and consider higher-order effects. Among the analyzed events, we present analyses of 8 events for which either new solutions are identified or additional information is obtained. We find that the previous binary-source interpretations of 5 events are better interpreted by binary-lens models. These events include OGLE-2006-BLG-238, OGLE-2007-BLG-159, OGLE-2007-BLG-491, OGLE-2008-BLG-143, and OGLE-2008-BLG-210. With additional data covering caustic crossings, we detect finite-source effects for 6 events including OGLE-2006-BLG-215, OGLE-2006-BLG-238, OGLE-2006-BLG-450, OGLE-2008-BLG-143, OGLE-2008-BLG-210, and OGLE-2008-BLG-513. Among them, we are able to measure the Einstein radii of 3 events for which multi-band data are available. These events are OGLE-2006-BLG-238, OGLE-2008-BLG-210, and OGLE-2008-BLG-513. For OGLE-2008-BLG-143, we detect higher-order effect induced by the changes of the observer's position caused by the orbital motion of the Earth around the Sun. In addition, we present degenerate solutions resulting from the known close/wide or ecliptic degeneracy. Finally, we note that the masses of the binary companions of the lenses of OGLE-2006-BLG-450 and OGLE-2008-BLG-210 are in the brown-dwarf regime.

CMB spectral distortions are induced by Compton collisions with electrons. We review the various schemes to characterize the anisotropic CMB with a non-Planckian spectrum. We advocate using logarithmically averaged temperature moments as the preferred language to describe these spectral distortions, both for theoretical modeling and observations. Numerical modeling is simpler, the moments are frame-independent, and in terms of scattering the mode truncation is exact.

Various extensions of standard inflationary models have been proposed recently by adding vector fields. Because they are generally motivated by large-scale anomalies, and the possibility of statistical anisotropy of primordial fluctuations, such models require to introduce non-standard couplings between vector fields on the one hand, and either gravity or scalar fields on the other hand. In this article, we study models involving a vector field coupled to a scalar field. We derive restrictive necessary conditions for these models to be both stable (Hamiltonian bounded by below) and causal (hyperbolic equations of motion).

We compute the spectral distortions of the Cosmic Microwave Background (CMB) polarization induced by non-linear effects in the Compton interactions between CMB photons and cold intergalactic electrons. This signal is of the $y$-type and is dominated by contributions arising from the reionized era. We stress that it is not shadowed by the thermal SZ effect which has no equivalent for polarization. We decompose its angular dependence into $E$- and $B$-modes, and we calculate the corresponding power spectra, both exactly and using a suitable Limber approximation that allows a simpler numerical evaluation. We find that $B$-modes are of the same order of magnitude as $E$-modes. Both spectra are relatively flat, peaking around $\ell=280$, and their overall amplitude is directly related to the optical depth to reionization. Moreover, we find this effect to be one order of magnitude larger than the non-linear kinetic Sunyaev-Zel'dovich effect in galaxy clusters. Finally, we discuss how to improve the detectability of our signal by cross-correlating it with other quantities sourced by the flow of intergalactic electrons.

In this paper, we develop in detail a fully geometrical method for deriving perturbation equations about a spatially homogeneous background. This method relies on the 3+1 splitting of the background space-time and does not use any particular set of coordinates: it is implemented in terms of geometrical quantities only, using the tensor algebra package xTensor in the xAct distribution along with the extension for perturbations xPert. Our algorithm allows one to obtain the perturbation equations for all types of homogeneous cosmologies, up to any order and in all possible gauges. As applications, we recover the well-known perturbed Einstein equations for Friedmann-Lemaitre-Robertson-Walker cosmologies up to second order and for Bianchi I cosmologies at first order. This work paves the way to the study of these models at higher order and to that of any other perturbed Bianchi cosmologies, by circumventing the usually too cumbersome derivation of the perturbed equations.

Gravitational microlensing events produced by lenses composed of binary masses are important because they provide a major channel to determine physical parameters of lenses. In this work, we analyze the light curves of two binary-lens events OGLE-2006-BLG-277 and OGLE-2012-BLG-0031 for which the light curves exhibit strong deviations from standard models. From modeling considering various second-order effects, we find that the deviations are mostly explained by the effect of the lens orbital motion. We also find that lens parallax effects can mimic orbital effects to some extent. This implies that modeling light curves of binary-lens events not considering orbital effects can result in lens parallaxes that are substantially different from actual values and thus wrong determinations of physical lens parameters. This demonstrates the importance of routine consideration of orbital effects in interpreting light curves of binary-lens events. It is found that the lens of OGLE-2006-BLG-277 is a binary composed of a low-mass star and a brown dwarf companion.

In the context of cosmological perturbation theory, we derive the second order Boltzmann equation describing the evolution of the distribution function of radiation without a specific gauge choice. The essential steps in deriving the Boltzmann equation are revisited and extended given this more general framework: i) the polarisation of light is incorporated in this formalism by using a tensor-valued distribution function; ii) the importance of a choice of the tetrad field to define the local inertial frame in the description of the distribution function is emphasized; iii) we perform a separation between temperature and spectral distortion, both for the intensity and for polarisation for the first time; iv) the gauge dependence of all perturbed quantities that enter the Boltzmann equation is derived, and this enables us to check the correctness of the perturbed Boltzmann equation by explicitly showing its gauge-invariance for both intensity and polarization. We finally discuss several implications of the gauge dependence for the observed temperature.

This article derives a multipolar hierarchy for the propagation of the weak-lensing shear and convergence in a general spacetime. The origin of B-modes, in particular on large angular scales, is related to the local isotropy of space. Known results assuming a Friedmann-Lema\^itre background are naturally recovered. The example of a Bianchi I spacetime illustrates our formalism and its implications for future observations are stressed.

Local non-Gaussianity, parametrized by $f_{\rm NL}$, introduces a scale-dependent bias that is strongest at large scales, precisely where General Relativistic (GR) effects also become significant. With future data, it should be possible to constrain $f_{\rm NL} = {\cal O}(1)$ with high redshift surveys. GR corrections to the power spectrum and ambiguities in the gauge used to define bias introduce effects similar to $f_{\rm NL}= {\cal O}(1)$, so it is essential to disentangle these effects. For the first time in studies of primordial non-Gaussianity, we include the consistent GR calculation of galaxy power spectra, highlighting the importance of a proper definition of bias. We present observable power spectra with and without GR corrections, showing that an incorrect definition of bias can mimic non-Gaussianity. However, these effects can be distinguished by their different redshift and scale dependence, so as to extract the true primordial non-Gaussianity.

The CMB bispectrum generated by second-order effects at recombination can be calculated analytically when one of the three modes has a wavelength much longer than the other two and is outside the horizon at recombination. This was pointed out in \cite{Creminelli:2004pv} and here we correct their results. We derive a simple formula for the bispectrum, $f_{NL}^{loc} = - (1/6+ \cos 2 \theta) \cdot (1- 1/2 \cdot d \ln (l_S^2 C_{S})/d \ln l_S)$, where $C_S$ is the short scale spectrum and $\theta$ the relative orientation between the long and the short modes. This formula is exact and takes into account all effects at recombination, including recombination-lensing, but neglects all late-time effects such as ISW-lensing. The induced bispectrum in the squeezed limit is small and will negligibly contaminate the Planck search for a local primordial signal: this will be biased only by $f_{NL}^{loc}\approx-0.4$. The above analytic formula includes the primordial non-Gaussianity of any single-field model. It also represents a consistency check for second-order Boltzmann codes: we find substantial agreement with the CMBquick code.

The non-linear evolution of the energy density of the radiation is shown to induce spectral distortions of the cosmic microwave background both at recombination and during the reionization era. This distortion has the same spectral signature as the one produced by the re-scattering of photons by non-relativistic hot electrons, the thermal Sunyaev-Zeldovich effect, whose amplitude is quantified by a Compton y parameter. A diffuse y-sky is then expected to emerge from mode couplings in the non-linear evolution of the cosmological perturbations and to superimpose to the point source contributions of galaxy clusters. The equations describing the evolution of the y field and a hierarchy governing its angular multipoles are derived from the second order Boltzmann equation. These equations are then integrated numerically to obtain the first predicted power spectrum of the diffuse y-sky. It is found to be a remarkable tracer of the reionization history of the Universe.

This article investigates the stability of a generic Kasner spacetime to linear perturbations, both at late and early times. It demonstrates that the perturbation of the Weyl tensor diverges at late time in all cases but in the particular one in which the Kasner spacetime is the product of a two-dimensional Milne spacetime and a two-dimensional Euclidean space. At early times, the perturbation of the Weyl tensor also diverges unless one imposes a condition on the perturbations so as to avoid the most divergent modes to be excited.

Nonlinear dynamics creates vortical currents when the tight-coupling approximation between photons and baryons breaks down around the time of recombination. This generates a magnetic field at second order in cosmological perturbations, whose power spectrum is fixed by standard physics, without the need for any ad hoc assumptions. We present the fully relativistic calculation of the magnetic power spectrum, including the effects of metric perturbations, second-order velocity and the photon anisotropic stress, thus generalizing and correcting previous results. We also show that significant magnetogenesis continues to occur after recombination. The power spectrum $\sqrt{k^3 P_B} $ behaves as $ \propto k^4$ on large scales, and $\propto k^{0.5}$ on small scales, down to $\sim 1\,$Mpc. On cluster scales, the created field has strength $\sim 3\times 10^{-29}$ Gauss.

This article constructs flat-sky approximations in a controlled way in the context of the cosmic microwave background observations for the computation of both spectra and bispectra. For angular spectra, it is explicitly shown that there exists a whole family of flat-sky approximations of similar accuracy for which the expression and amplitude of next to leading order terms can be explicitly computed. It is noted that in this context two limiting cases can be encountered for which the expressions can be further simplified. They correspond to cases where either the sources are localized in a narrow region (thin-shell approximation) or are slowly varying over a large distance (which leads to the so-called Limber approximation). Applying this to the calculation of the spectra it is shown that, as long as the late integrated Sachs-Wolfe contribution is neglected, the flat-sky approximation at leading order is accurate at 1% level for any multipole. Generalization of this construction scheme to the bispectra led to the introduction of an alternative description of the bispectra for which the flat-sky approximation is well controlled. This is not the case for the usual description of the bispectrum in terms of reduced bispectrum for which a flat-sky approximation is proposed but the next-to-leading order terms of which remain obscure.

The tight-coupling approximation (TCA) used to describe the early dynamics of the baryons-photons system is systematically built to higher orders in the inverse of the interaction rate. This expansion can be either used to grasp the physical effects by deriving simple analytic solutions or to obtain a form of the system which is stable numerically at early times. In linear cosmological perturbations, we estimate numerically its precision, and we discuss the implications for the baryons acoustic oscillations. The TCA can be extended to the second order cosmological perturbations, and in particular we recover that vorticity is not generated at lowest order of this expansion.

This article presents the first computation of the complete bispectrum of the cosmic microwave background temperature anisotropies arising from the evolution of all cosmic fluids up to second order, including neutrinos. Gravitational couplings, electron density fluctuations and the second order Boltzmann equation are fully taken into account. Comparison to limiting cases that appeared previously in the literature are provided. These are regimes for which analytical insights can be given. The final results are expressed in terms of equivalent fNL for different configurations. It is found that for moments up to lmax=2000, the signal generated by non-linear effects is equivalent to fNL~5 for both local-type and equilateral-type primordial non-Gaussianity.

This article investigates the full Boltzmann equation up to second order in the cosmological perturbations. Describing the distribution of polarized radiation by using a tensor valued distribution function, the second order Boltzmann equation, including polarization, is derived without relying on the Stokes parameters.

This article provides a general study of the Hamiltonian stability and the hyperbolicity of vector field models involving both a general function of the Faraday tensor and its dual, $f(F^2,F\tilde F)$, as well as a Proca potential for the vector field, $V(A^2)$. In particular it is demonstrated that theories involving only $f(F^2)$ do not satisfy the hyperbolicity conditions. It is then shown that in this class of models, the cosmological dynamics always dilutes the vector field. In the case of a nonminimal coupling to gravity, it is established that theories involving $R f(A^2)$ or $Rf(F^2)$ are generically pathologic. To finish, we exhibit a model where the vector field is not diluted during the cosmological evolution, because of a nonminimal vector field-curvature coupling which maintains second-order field equations. The relevance of such models for cosmology is discussed.

Anisotropic and isotropic slow-roll inflation supported by n-forms are sought. Canonical field strengths and their duals are taken into account, and they are allowed to have a potential and also, when necessary for slow-roll, a nonminimal curvature coupling. New isotropic solutions are found for three-forms. It is also shown that some n-form actions are equivalent to f(R) gravity and scalar field models with possible nonminimal couplings. Anisotropic solutions are found for two-forms, generalising vector inflation, however, also the two-form is unstable during inflation due to the nonminimal coupling to curvature. We also address the stability of the isotropic background solutions supported by a triad of vectors .

This article investigates the non-linear evolution of cosmological perturbations on sub-Hubble scales in order to evaluate the unavoidable deviations from Gaussianity that arise from the non-linear dynamics. It shows that the dominant contribution to modes coupling in the cosmic microwave background temperature anisotropies on small angular scales is driven by the sub-Hubble non-linear evolution of the dark matter component. The perturbation equations, involving in particular the first moments of the Boltzmann equation for the photons, are integrated up to second order in perturbations. An analytical analysis of the solutions gives a physical understanding of the result as well as an estimation of its order of magnitude. This allows to quantify the expected deviation from Gaussianity of the cosmic microwave background temperature anisotropy and, in particular, to compute its bispectrum on small angular scales. Restricting to equilateral configurations, we show that the non-linear evolution accounts for a contribution that would be equivalent to a constant primordial non-Gaussianity of order fNL~25 on scales ranging approximately from l~1000 to l~3000.

This article investigates the predictions of an inflationary phase starting from a homogeneous and anisotropic universe of the Bianchi I type. After discussing the evolution of the background spacetime, focusing on the number of e-folds and the isotropization, we solve the perturbation equations and predict the power spectra of the curvature perturbations and gravity waves at the end of inflation. The main features of the early anisotropic phase is (1) a dependence of the spectra on the direction of the modes, (2) a coupling between curvature perturbations and gravity waves, and (3) the fact that the two gravity waves polarisations do not share the same spectrum on large scales. All these effects are significant only on large scales and die out on small scales where isotropy is recovered. They depend on a characteristic scale that can, but a priori must not, be tuned to some observable scale. To fix the initial conditions, we propose a procedure that generalises the one standardly used in inflation but that takes into account the fact that the WKB regime is violated at early times when the shear dominates. We stress that there exist modes that do not satisfy the WKB condition during the shear-dominated regime and for which the amplitude at the end of inflation depends on unknown initial conditions. On such scales, inflation loses its predictability. This study paves the way to the determination of the cosmological signature of a primordial shear, whatever the Bianchi I spacetime. It thus stresses the importance of the WKB regime to draw inflationary predictions and demonstrates that when the number of e-folds is large enough, the predictions converge toward those of inflation in a Friedmann-Lemaitre spacetime but that they are less robust in the case of an inflationary era with a small number of e-folds.

This article investigates the collisionless Boltzmann equation up to second order in the cosmological perturbations. It describes the gauge dependence of the distribution function and the construction of a gauge invariant distribution function and brightness, and then derives the gauge invariant fluid limit.

This article describes the theory of cosmological perturbations around a homogeneous and anisotropic universe of the Bianchi I type. Starting from a general parameterisation of the perturbed spacetime a la Bardeen, a complete set of gauge invariant variables is constructed. Three physical degrees of freedom are identified and it is shown that, in the case where matter is described by a scalar field, they generalize the Mukhanov-Sasaki variables. In order to show that they are canonical variables, the action for the cosmological perturbations at second order is derived. Two major physical imprints of the primordial anisotropy are identified: (1) a scalar-tensor ``see-saw'' mechanism arising from the fact that scalar, vector and tensor modes do not decouple and (2) an explicit dependence of the statistical properties of the density perturbations and gravity waves on the wave-vector instead of its norm. This analysis extends, but also sheds some light on, the quantization procedure that was developed under the assumption of a Friedmann-Lemaitre background spacetime, and allows to investigate the robustness of the predictions of the standard inflationary scenario with respect to the hypothesis on the symmetries of the background spacetime. These effects of a primordial anisotropy may be related to some anomalies of the cosmic microwave background anisotropies on large angular scales.

This note derives the analogue of the Mukhanov-Sasaki variables both for scalar and tensor perturbations in the 1+3 covariant formalism. The possibility of generalizing them to non-flat Friedmann-Lemaitre universes is discussed.

The generation of gravitational waves during inflation due to the non-linear coupling of scalar and tensor modes is discussed. Two methods describing gravitational wave perturbations are used and compared: a covariant and local approach, as well as a metric-based analysis based on the Bardeen formalism. An application to slow-roll inflation is also described.