We develop a framework to describe gravitational wave propagation through a stochastic distribution of weak gravitational lenses beyond the geometric optics limit. We model the lens distribution as a static random background field and formulate the problem in the language of quenched disorder, treating the disorder averaged density matrix as the fundamental object from which observables are computed. Using the Schwinger Keldysh formalism, we construct a path-integral representation of the averaged density matrix and derive its explicit form perturbatively for a suitable class of couplings. The result naturally separates into a quadratic exponential term, which governs the suppression of phase sensitive contributions in the averaged description, and a purely oscillatory contribution, which modifies coherent propagation through a disorder-induced correction to the propagation kernel. This provides a unified description of interference, diffraction, and statistical fluctuations of the lens distribution within a single framework. We also identify the physical scales controlling the onset of coherence loss and illustrate the formalism in the case of Gaussian wave packets. More generally, the derivation applies to any system described by the same class of actions, making the framework relevant beyond gravitational wave lensing to wave propagation in disordered media.

We present a complete computation of the scalar power spectrum in the inflation without inflaton (IWI) framework, where the inflationary expansion is driven solely by a de Sitter (dS) background and scalar fluctuations arise as second-order effects sourced by tensor perturbations. By explicitly deriving and numerically integrating the full second-order kernel of the Einstein equations, we obtain a scale-invariant scalar spectrum without invoking a fundamental scalar field. In this framework, the amplitude of the scalar fluctuations is directly linked to the scale of inflation. More precisely, we show that matching the observed level of scalar fluctuations Delta 2 phi (k*) approximate to 10-9 at cosmic microwave background scales fixes the inflationary energy scale Hinf as a function of the number of observed e-folds Nobs. For Nobs 30-60, we find Hinf 5 x 1013 GeV-2 x 1010 GeV, corresponding to a tensor-to-scalar ratio r 0.01-5 x 10-9. In particular, requiring consistency with instantaneous reheating, we predict a number of e-foldings O(50) and an inflationary scale Hinf 1011 GeV. We also incorporate in our framework the quantum break time of the dS state, and we show that it imposes an upper bound on the number of particle species. Specifically, using laboratory constraints on the number of species limits the duration of inflation to Nobs less than or similar to 126 e-folds. These results establish the IWI scenario as a predictive and falsifiable alternative to standard inflaton-driven models, linking the observed amplitude of primordial fluctuations directly to the quantum nature and finite lifetime of dS space.

We extend the Schrödinger approach to large-scale structure formation beyond the Newtonian regime by working at first post-Friedmann (1PF) order. The standard Schrödinger–Poisson system gives a useful reformulation of the dynamics of a self-gravitating pressureless fluid, but it corresponds to the leading post-Friedmann, or Newtonian, limit. It therefore misses the relativistic corrections that enter at next-to-leading order and become relevant on horizon scales and for high-precision cosmological surveys. Starting from the 1PF continuity and Euler equations in a flat ΛCDM background, we identify the conserved density variable associated with covariant mass conservation. In terms of this variable, the continuity equation takes a Newtonian-like conservative form. However, even for vanishing covariant vorticity, the spatial velocity field in the cosmological frame contains a transverse 1PF component. Thus the full 1PF mass flux cannot be represented solely by the gradient of a scalar phase. We show that the Schrödinger-like formulation at 1PF order requires an effective vector potential fixed by this transverse velocity component. This vector potential contains the post-Friedmann metric vector perturbation, related to relativistic frame-dragging effects, together with nonlinear scalar terms required by the zero-vorticity condition. Equivalently, when the equation is written in scalar form, these corrections appear as an imaginary contribution to the effective potential. At leading order our system reduces to the usual Schrödinger–Poisson formulation, while at 1PF order it provides a relativistic extension of the Schrödinger description of cold matter dynamics.

Over the next two decades, gravitational-wave (GW) observations are expected to evolve from a discovery-driven endeavour into a precision tool for astrophysics, cosmology, and fundamental physics. Current second-generation ground-based detectors have established the existence of compact-binary mergers and enabled GW multi-messenger astronomy, but they remain limited in sensitivity, redshift reach, frequency coverage, and duty cycle. These limitations prevent them from addressing many fundamental open questions in cosmology. By the 2040s, wide-field electromagnetic surveys will have mapped the luminous Universe with unprecedented depth and accuracy. Nevertheless, key problems including the nature of dark matter, the physical origin of cosmic acceleration, the properties of gravity on cosmological scales, and the physical conditions of the earliest moments after the Big Bang will remain only partially constrained by electromagnetic observations alone. Progress on these fronts requires access to physical processes and epochs that do not emit light. Gravitational waves provide a unique and complementary observational channel: they propagate over cosmological distances largely unaffected by intervening matter, probe extreme astrophysical environments, and respond directly to the geometry of spacetime. In this context, next-generation GW observatories such as the Einstein Telescope (ET) will be transformative for European astronomy. Operating at sensitivities and frequencies beyond existing detectors, ET will observe binary black holes and neutron stars out to previously inaccessible redshifts, enable continuous high signal-to-noise monitoring of compact sources, and detect gravitational-wave backgrounds of astrophysical and cosmological origin. Together with space-based detectors, ET will play a central role in advancing our understanding of cosmic evolution and fundamental physics.

We investigate static, spherically symmetric halo configurations within Unified Dark Matter (UDM) scalar-field models, developing a systematic mapping between standard cold dark matter (CDM) density profiles and their UDM counterparts. Exploiting the equivalence-class structure of UDM models, we show that, in principle, different Lagrangian realisations can share the same weak-field rotation curve while exhibiting distinct field properties. We reconstruct the effective energy density, radial and tangential pressures from a phenomenological circular velocity profile, ensuring the absence of ghosts and instabilities and the preservation of the Null Energy Condition (NEC). Applying our procedure to several commonly used CDM halo profiles – including Persic, Salucci & Stel, NFW, and Burkert models – we demonstrate that their phenomenological success can be retained within a relativistic UDM framework, reproducing the observed flatness of rotation curves without introducing separate dark matter and dark energy components.

We investigate primordial non-Gaussianity in the Inflation without an Inflaton (IWI) framework, where scalar perturbations are generated at second order by primordial gravitational waves in Einstein gravity on an exact de Sitter (dS) background. Since scalar modes are produced nonlinearly from tensor modes, non-Gaussianity is an intrinsic prediction of the mechanism. We compute the corresponding scalar bispectrum, derive the relevant contribution to the three-point function of the scalar potential, and evaluate its shape numerically. We find that, unlike the scalar power spectrum, the bispectrum depends logarithmically on the ultraviolet cutoff set by the end of inflation, indicating a structural difference between the two- and three-point statistics in this scenario. Its shape is enhanced toward squeezed configurations, but its amplitude becomes strongly suppressed once the scalar power spectrum is normalized to the observed value. The resulting non-Gaussianity at CMB scales is therefore negligibly small, well below present observational sensitivity.

We employ the quantum Boltzmann equation to analyze the gravitationally generated entanglement between two remote qubits by considering two explicit microscopic models. A graviton propagator is employed as the mediator of the interactions, while the qubits are considered in a spatial superposition state. Such a setup, in the case of any entanglement generation, could potentially offer experimental evidence for the quantization of gravity. By treating the qubits as spin-1/2 particles in wave packets, we establish that the entanglement arises from forward scattering processes involving graviton exchanges. In our study, we consider both static and dynamical limits of the propagator and show that only in the dynamical limit such entangled states can be generated. These effects are observed to be diminished as the wave packet size increases. Our findings sheds more light into the gravity mediated entanglement between two spin-1/2 particles.

Primordial gravitational waves (PGWs) generate scalar density perturbations at second order. Since the induced density contrast is quadratic in the tensor field, it is intrinsically non-Gaussian. We study the imprint of this tensor-induced non-Gaussianity (NG) on the large-scale clustering of dark matter halos through its correction to halo bias. Focusing on inflationary scenarios with a peaked primordial tensor spectrum, we derive the leading scale-dependent contribution sourced by the bispectrum of the induced density field. While yielding a percent-level bias correction for massive low-redshift halos, this effect can reach an 𝒪(1) modulation for rare, high-redshift halos at z=7. Notably, the resulting signature exhibits a distinct scale dependence that is not captured by standard primordial non-Gaussianity (PNG) templates. Our results establish halo bias as a novel probe of PGWs through their imprint on the large-scale structure, providing a complementary window into the inflationary epoch.

Higher derivative quantum corrections are essential components of scalar tensor effective field theories (EFTs), yet they typically reintroduce the Ostrogradsky ghost instability that the classical theory was designed to evade. This paper resolves this fundamental tension by establishing a rigorous equivalence between two distinct criteria for theoretical consistency. We analyze a general DHOST theory augmented by Gauss Bonnet and Weyl squared operators with coefficients that are arbitrary functions of the scalar field and its kinetic term. We then pursue two independent paths: first, we derive a set of differential equations for these coefficients by demanding that the full action remains invariant under the protective gauge symmetry of the classical theory. Second, we perform a first principles Hamiltonian analysis using the ADM formalism, deriving a separate set of conditions by imposing the primary and secondary constraints required to eliminate the ghost. Our central result is the proof that these two sets of conditions, one algebraic and one dynamical, are mathematically identical. This equivalence demonstrates that the gauge symmetry is the fundamental origin of Hamiltonian stability in the quantum corrected theory and establishes the symmetry principle as a powerful and practical tool for constructing consistent, ghost free gravitational EFTs without resorting to a full Hamiltonian analysis

This chapter explores theoretical and observational strategies to use the stochastic gravitational-wave background detectable by Square Kilometre Array Observatory (SKAO) as a probe of precision cosmology. We detail the critical phenomenon of scalar-induced gravitational waves, demonstrating their unique features and their sensitivity to primordial non-Gaussianity on scales much smaller than those probed by the cosmic microwave background and large-scale structure. We investigate the phenomenology of parity violation in the early Universe through the chirality imprinted in the stochastic gravitational-wave background, demonstrating that a parity-odd primordial trispectrum can generate a detectable scale-dependent helicity. On the observational side, we point out that the standard gravitational-wave angular auto-correlation analysis is significantly limited by astrophysical shot noise. We show that cross-correlating the gravitational wave signal with independent large-scale structure tracers enhances the signal-to-noise ratio and allows us to isolate the astrophysical and cosmological components in the background. These results can be achieved only thanks to the enhanced sensitivity of SKAO, extensive sky coverage, and high angular resolution, which together make such targets observationally feasible.

We investigate the non-Gaussianity of second-order matter density perturbations induced by primordial gravitational waves (GWs). These tensor-induced scalar modes arise from local fluctuations in the GWs energy density, which is quadratic in tensor perturbations. The resulting second-order density contrast follows a chi-squared distribution, naturally exhibiting significant non-Gaussianity. We compute the bispectrum of these tensor-induced scalar modes and analyze its dependence on various primordial GWs power spectra, including scale-invariant, blue-tilted, Gaussian-bump, and monochromatic sources. We find that the bispectrum shape is inherently sensitive to the underlying GWs spectrum by construction. In particular, Gaussian-bump and monochromatic sources produce a strong signal peaking in the equilateral configuration, similar to the effect of scalar-induced tensor modes. Our findings reveal a new way to probe primordial GWs via galaxy surveys and highlight a unique feature of tensor-induced density perturbations, otherwise mimicking linear ones on sub-horizon scales.

Scalar-Induced Gravitational Waves (SIGWs), second-order tensor modes sourced by first-order scalar fluctuations in General Relativity (GR), are expected to contribute to the Stochastic Gravitational Wave Background (SGWB) potentially detectable by current and future gravitational wave interferometers. In the framework of GR, this SGWB represents an unavoidable contribution to the gravitational wave spectrum. In this paper, we go beyond GR and we investigate the behavior of SIGWs in f(R) gravity. We explore how the SIGW spectrum is influenced across a broad range of frequencies, from the nano-Hz regime, where the Pulsar Timing Array (PTA) operates, through the milli-Hz band probed by the space-based LISA detector, up to the kilo-Hz frequency range, where the ground-based LIGO/Virgo/KAGRA network is currently operational. Our results indicate that the beyond-GR correction leaves an observational imprint, mainly in the low-frequency part of the spectrum, giving the possibility to use SIGW to constrain GR on scales on which we have limited information.

We propose a novel scenario in which scalar perturbations, which seed the large-scale structure of the universe, are generated without relying on a scalar field (the inflaton). In this framework, inflation is driven by a de Sitter space time, where tensor metric fluctuations (i.e., gravitational waves) naturally arise from quantum vacuum oscillations, and scalar fluctuations are generated via second-order tensor effects. We compute the power spectrum of such scalar fluctuations and show it to be consistent with near scale invariance. We derive the necessary conditions under which scalar perturbations become significant and much larger than the tensor modes, and we identify a natural mechanism to end inflation via a transition to a radiation-dominated phase. Our proposed mechanism could remove the need for a model-dependent scenario: the choice of a scalar field, as the inflaton, to drive inflation.

The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance.

A variety of mechanisms in the early Universe lead to the generation of gravitational waves (GWs). We introduce here a novel source of GWs generated by vacuum fluctuations after inflation. Being gravitons minimally coupled particles, their quantum creation takes place during inflation, but is absent in an unperturbed Universe during the radiation-dominated epoch, since they behave as conformally invariant particles. However, the presence of inhomogeneities breaks the conformal flatness of the metric, allowing scalar metric perturbations to induce the quantum production of gravitons. We compute the resulting GW spectrum from this mechanism for different models of the primordial scalar power spectrum. We find that this GW signal peaks around the GHz frequency range, distinguishing it from other astrophysical and cosmological backgrounds and underscoring the need for detectors sensitive to these high frequencies.

The Cosmological Gravitational Wave Background (CGWB) anisotropies contain valuable information about the physics of the early universe. Given that General Relativity is intrinsically non-linear, it is important to look beyond first-order contributions in cosmological perturbations. In this work, we present a non-perturbative approach for the computation of the CGWB anisotropies at large scales, providing the extension of the initial conditions and the Sachs-Wolfe effect for the CGWB, which encodes the full non-linearity of the scalar metric perturbations. We also derive the non-perturbative expression for three-point correlation of the gravitational wave energy density perturbation in the case of an inflationary CGWB with a scale-invariant power spectrum and negligible primordial non-Gaussianity. We show that, under such conditions, the gravitational wave energy density perturbations are lognormally distributed, leading to an interesting effect such as intermittency.

The detection of primordial B modes of the cosmic microwave background (CMB) could provide information about the early stages of the Universe's evolution. The faintness of this signal requires exquisite calibration accuracy and control of instrumental systematic effects which otherwise could bias the measurements. In this work, we study the impact of an imperfect relative polarisation gain calibration on the recovered value of the tensor-toscalar ratio r for the LiteBIRD experiment, through the application of the blind Needlet Internal Linear Combination (NILC) foreground-cleaning method. We derive requirements on the relative calibration accuracy of the overall polarisation gain (Delta g(v)) for each LiteBIRD frequency channel. Our results show that minimum variance techniques, as NILC, are less sensitive to systematic gain calibration uncertainties compared to a parametric approach, if the latter is not equipped with a proper modelling of these instrumental effects. In this study, the most stringent requirements are found in the channels where the CMB signal is relatively brighter, with the tightest constraints at 166GHz (Delta g(v) similar to 0.16%). This differs from the outcome of an analogous analysis performed with a parametric method, where the tightest requirements are obtained for the foreground-dominated channels. Gain calibration uncertainties, corresponding to the derived requirements, are then simultaneously propagated into all frequency channels. By doing so, we find that the overall impact on estimated r is lower than the total gain systematic budget for LiteBIRD approximately by a factor 5, due to the correlations of the impacts of gain calibration uncertainties in different frequency channels. In order to decouple the systematic effect from the specific choice of the model, we derive the requirements assuming constant spectral parameters for the foreground emission. To assess the robustness of the obtained results against more realistic scenarios, we repeat the analysis assuming sky models of intermediate and high complexity. In these further cases, we adopt an optimised NILC pipeline, called the Multi-Clustering NILC (MC-NILC). We find that the impact of gain calibration uncertainties on r is lower than the LiteBIRD gain systematics budget for the intermediate-complexity sky model. For the high-complexity case, instead, it would be necessary to tighten the requirements by a factor 1.8.

We explore the capability of measuring lensing signals in $LiteBIRD$ full-sky polarization maps. With a $30$ arcmin beam width and an impressively low polarization noise of $2.16\,\mu$K-arcmin, $LiteBIRD$ will be able to measure the full-sky polarization of the cosmic microwave background (CMB) very precisely. This unique sensitivity also enables the reconstruction of a nearly full-sky lensing map using only polarization data, even considering its limited capability to capture small-scale CMB anisotropies. In this paper, we investigate the ability to construct a full-sky lensing measurement in the presence of Galactic foregrounds, finding that several possible biases from Galactic foregrounds should be negligible after component separation by harmonic-space internal linear combination. We find that the signal-to-noise ratio of the lensing is approximately $40$ using only polarization data measured over $90\%$ of the sky. This achievement is comparable to $Planck$'s recent lensing measurement with both temperature and polarization and represents a four-fold improvement over $Planck$'s polarization-only lensing measurement. The $LiteBIRD$ lensing map will complement the $Planck$ lensing map and provide several opportunities for cross-correlation science, especially in the northern hemisphere.

We present a study of the impact of a beam far side-lobe lack of knowledge on the measurement of the Cosmic Microwave Background B-mode signal at large scale. Beam far side-lobes induce a mismatch in the transfer function of Galactic foregrounds between the dipole and higher multipoles which degrads the performances of component separation methods. This leads to foreground residuals in the CMB map. It is expected to be one of the main source of systematic effects in future CMB polarization observations. Thus, it becomes crucial for all-sky survey missions to take into account the interplays between beam systematic effects and all the data analysis steps. LiteBIRD is the ISAS/JAXA second strategic large-class satellite mission and is dedicated to target the measurement of CMB primordial B modes by reaching a sensitivity on the tensor-to-scalar ratio r of sigma (r) <= 10(-3) assuming r = 0. The primary goal of this paper is to provide the methodology and develop the framework to carry out the end-to-end study of beam far side-lobe effects for a space-borne CMB experiment. We introduce uncertainties in the beam model, and propagate the beam effects through all the steps of the analysis pipeline, most importantly including component separation, up to the cosmological results in the form of a bias delta r. As a demonstration of our framework, we derive requirements on the calibration and modeling for the LiteBIRD's beams under given assumptions on design, simulation, component separation method and allocated error budget. In particular, we assume a parametric method of component separation with no mitigation of the far side-lobes effect at any stage of the analysis pipeline. We show that delta r is mostly due to the integrated fractional power difference between the estimated beams and the true beams in the far side-lobes region, with little dependence on the actual shape of the beams, for low enough delta r. Under our set of assumptions, in particular considering the specific foreground cleaning method we used, we find that the integrated fractional power in the far side-lobes should be known at the level of similar to 10(-4), to achieve the required limit on the bias delta r < 1.9 x 10(-5). The framework and tools developed for this study can be easily adapted to provide requirements under different design, data analysis frameworks and for other future space-borne experiments, such as PICO or CMB-Bharat. We further discuss the limitations of this framework and potential extensions to circumvent them.

The non-thermal nature of the stochastic gravitational-wave background of cosmological origin (CGWB) poses a challenge in defining the initial conditions for the graviton overdensity. Specifically, the adiabatic initial condition, which holds for Cosmic Microwave Background (CMB) photons, is not guaranteed a priori for the primordial GWs. In this letter, we compute the initial conditions for the cosmological background generated by quantum fluctuations of the metric during inflation. Our analysis reveals that adiabatic initial conditions are no longer valid. The violation of adiabaticity arises from the presence of independent tensor perturbations during inflation, which behave as two extra fields that affect the standard single-clock argument. Since the energy density of the CGWB is subdominant compared to ordinary matter, gravitational radiation plays a negligible role in Einstein's equations. Therefore, the only way to compute the initial conditions is to perturb the energy-momentum tensor defined in terms of the gravitational strain. A direct consequence of our finding is that the initial conditions from inflation enhance the total CGWB angular power-spectrum by an order of magnitude compared to the standard adiabatic case.