The Stochastic Gravitational-Wave Background (SGWB) is expected to be a key observable for Gravitational-Wave (GW) interferometry. Its detection will open a new window on early Universe cosmology, on the astrophysics of compact objects and, as shown in this Letter, on the particle physics content of the Universe. In this Letter we show that, besides their effects on the Cosmic Microwave Background (CMB) and on Large Scale Structure (LSS), relativistic particles in the early Universe leave a clear imprint on the anisotropies of the SGWB. In particular we show that a change in the number of decoupled relativistic particles shifts the angular power spectrum of the SGWB, as both the Sachs-Wolfe (SW) and the Integrated Sachs-Wolfe (ISW) terms are affected. Being very large-angle effects, these lead to new testable predictions for future GW interferometers.

Given the recent advances in gravitational-wave detection technologies, the detection and characterisation of gravitational-wave backgrounds (GWBs) with the Laser Interferometer Space Antenna (LISA) is a real possibility. To assess the abilities of the LISA satellite network to reconstruct anisotropies of different angular scales and in different directions on the sky, we develop a map-maker based on an optimal quadratic estimator. The resulting maps are maximum likelihood representations of the GWB intensity on the sky integrated over a broad range of frequencies. We test the algorithm by reconstructing known input maps with different input distributions and over different frequency ranges. We find that, in an optimal scenario of well understood noise and high frequency, high SNR signals, the maximum scales LISA may probe are $\ell_{\rm max} \lesssim 15$. The map-maker also allows to test the directional dependence of LISA noise, providing insight on the directional sky sensitivity we may expect.

The detection and characterization of the Stochastic Gravitational Wave Background (SGWB) is one of the main goals of Gravitational Wave (GW) experiments. The observed SGWB will be the combination of GWs from cosmological (as predicted by many models describing the physics of the early Universe) and astrophysical origins, which will arise from the superposition of GWs from unresolved sources whose signal is too faint to be detected. Therefore, it is important to have a proper modeling of the astrophysical SGWB (ASGWB) in order to disentangle the two signals; moreover, this will provide additional information on astrophysical properties of compact objects. Applying the Cosmic Rulers formalism, we compute the observed ASGWB angular power spectrum, hence using gauge invariant quantities, accounting for all effects intervening between the source and the observer. These are the so-called projection effects, which include Kaiser, Doppler and gravitational potentials effect. Our results show that these projection effects are the most important at the largest scales, and they contribute to up to tens of percent of the angular power spectrum amplitude, with the Kaiser term being the largest at all scales. While the exact impact of these results will depend on instrumental and astrophysical details, a precise theoretical modeling of the ASGWB will necessarily need to include all these projection effects.

The Einstein Telescope (ET), a proposed European ground-based gravitational-wave detector of third-generation, is an evolution of second-generation detectors such as Advanced LIGO, Advanced Virgo, and KAGRA which could be operating in the mid 2030s. ET will explore the universe with gravitational waves up to cosmological distances. We discuss its main scientific objectives and its potential for discoveries in astrophysics, cosmology and fundamental physics.

An observable stochastic background of gravitational waves is generated whenever primordial black holes are created in the early universe thanks to a small-scale enhancement of the curvature perturbation. We calculate the anisotropies and non-Gaussianity of such stochastic gravitational waves background which receive two contributions, the first at formation time and the second due to propagation effects. The former contribution can be generated if the distribution of the curvature perturbation is characterized by a local and scale-invariant shape of non-Gaussianity. Under such an assumption, we conclude that a sizeable magnitude of anisotropy and non-Gaussianity in the gravitational waves would suggest that primordial black holes may not comply the totality of the dark matter.

The stochastic gravitational wave background (SGWB) contains a wealth of information on astrophysical and cosmological processes. A major challenge of upcoming years will be to extract the information contained in this background and to disentangle the contributions of different sources. In this paper we provide the formalism to extract, from the correlation of three signals in the Laser Interferometer Space Antenna (LISA), information about the tensor three-point function, which characterizes the non-Gaussian properties of the SGWB. This observable can be crucial to discriminate whether a SGWB has a primordial or astrophysical origin. Compared to the two-point function, the SGWB three-point function has a richer dependence on the gravitational wave momenta and chiralities. It can be used therefore as a powerful discriminator between different models. For the first time we provide the response functions of LISA to a general SGWB three-point function. As examples, we study in full detail the cases of an equilateral and squeezed SGWB bispectra, and provide the explicit form of the response functions, ready to be convoluted with any theoretical prediction of the bispectrum to obtain the observable signal. We further derive the optimal estimator to compute the signal-to-noise ratio. Our formalism covers general shapes of non-Gaussianity, and can be extended straightaway to other detector geometries. Finally, we provide a short overview of models of the early universe that can give rise to a non-Gaussian SGWB.

A future detection of the Stochastic Gravitational Wave Background (SGWB) with GW experiments is expected to open a new window on early universe cosmology and on the astrophysics of compact objects. In this paper we study SGWB anisotropies, that can offer new tools to discriminate between different sources of GWs. In particular, the cosmological SGWB inherits its anisotropies both (i) at its production and (ii) during its propagation through our perturbed universe. Concerning (i), we show that it typically leads to anisotropies with order one dependence on frequency. We then compute the effect of (ii) through a Boltzmann approach, including contributions of both large-scale scalar and tensor linearized perturbations. We also compute for the first time the three-point function of the SGWB energy density, which can allow one to extract information on GW non-Gaussianity with interferometers. Finally, we include non-linear effects associated with long wavelength scalar fluctuations, and compute the squeezed limit of the 3-point function for the SGWB density contrast. Such limit satisfies a consistency relation, conceptually similar to what found in the literature for the case of CMB perturbations.

Parity violating interactions in the early Universe can source a stochastic gravitational wave background (SGWB) with a net circular polarization. In this paper, we study possible ways to search for circular polarization of the SGWB with interferometers. Planar detectors are unable to measure the net circular polarization of an isotropic SGWB. We discuss the possibility of using the dipolar anisotropy kinematically induced by the motion of the solar system with respect to the cosmic reference frame to measure the net circular polarization of the SGWB with planar detectors. We apply this approach to LISA, re-assessing previous analyses by means of a more detailed computation and using the most recent instrument specifications, and to the Einstein Telescope (ET), estimating for the first time its sensitivity to circular polarization. We find that both LISA and ET, despite operating at different frequencies, could detect net circular polarization with a signal-to-noise ratio of order one in a SGWB with amplitude $h^2 \Omega_\text{GW} \simeq 10^{-11}$. We also investigate the case of a network of ground based detectors. We present fully analytical, covariant formulas for the detector overlap functions in the presence of circular polarization. Our formulas do not rely on particular choices of reference frame, and can be applied to interferometers with arbitrary angles among their arms.

Why is the Universe so homogeneous and isotropic? We summarize a general study of a $\gamma$-law perfect fluid alongside an inhomogeneous, massless scalar gauge field (with homogeneous gradient) in anisotropic spaces with General Relativity. The anisotropic matter sector is implemented as a $j$-form (field-strength level), where $j\,\in\,\{1,3\}$, and the spaces studied are Bianchi space-times of solvable type. Wald's no-hair theorem is extended to include the $j$-form case. We highlight three new self-similar space-times: the Edge, the Rope and Wonderland. The latter solution is so far found to exist in the physical state space of types I,II, IV, VI$_0$, VI$_h$, VII$_0$ and VII$_h$, and is a global attractor in I and V. The stability analysis of the other types has not yet been performed. This paper is a summary of ~[1], with some remarks towards new results which will be further laid out in upcoming work.

Since the very beginning of astronomy the location of objects on the sky has been a fundamental observational quantity that has been taken for granted. While precise two dimensional positional information is easy to obtain for observations in the electromagnetic spectrum, the positional accuracy of current and near future gravitational wave detectors is limited to between tens and hundreds of square degrees, which makes it extremely challenging to identify the host galaxies of gravitational wave events or to confidently detect any electromagnetic counterparts. Gravitational wave observations provide information on source properties and distances that is complementary to the information in any associated electromagnetic emission and that is very hard to obtain in any other way. Observing systems with multiple messengers thus has scientific potential much greater than the sum of its parts. A gravitational wave detector with higher angular resolution would significantly increase the prospects for finding the hosts of gravitational wave sources and triggering a multi-messenger follow-up campaign. An observatory with arcminute precision or better could be realised within the Voyage 2050 programme by creating a large baseline interferometer array in space and would have transformative scientific potential. Precise positional information of standard sirens would enable precision measurements of cosmological parameters and offer new insights on structure formation; a high angular resolution gravitational wave observatory would allow the detection of a stochastic background and resolution of the anisotropies within it; it would also allow the study of accretion processes around black holes; and it would have tremendous potential for tests of modified gravity and the discovery of physics beyond the Standard Model.

The Stochastic Gravitational Wave Background (SGWB) is expected to be a key observable for Gravitational Wave (GW) interferometry. Its detection will open a new window on early universe cosmology and on the astrophysics of compact objects. Using a Boltzmann approach, we study the angular anisotropies of the GW energy density, which is an important tool to disentangle the different cosmological and astrophysical contributions to the SGWB. Anisotropies in the cosmological background are imprinted both at its production, and by GW propagation through the large-scale scalar and tensor perturbations of the universe. The first contribution is not present in the Cosmic Microwave Background (CMB) radiation (as the universe is not transparent to photons before recombination), causing an order one dependence of the anisotropies on frequency. Moreover, we provide a new method to characterize the cosmological SGWB through its possible deviation from a Gaussian statistics. In particular, the SGWB will become a new probe of the primordial non-Gaussianity of the large-scale cosmological perturbations.

We present a set of tools to assess the capabilities of LISA to detect and reconstruct the spectral shape and amplitude of a stochastic gravitational wave background (SGWB). We first provide the LISA power-law sensitivity curve and binned power-law sensitivity curves, based on the latest updates on the LISA design. These curves are useful to make a qualitative assessment of the detection and reconstruction prospects of a SGWB. For a quantitative reconstruction of a SGWB with arbitrary power spectrum shape, we propose a novel data analysis technique: by means of an automatized adaptive procedure, we conveniently split the LISA sensitivity band into frequency bins, and fit the data inside each bin with a power law signal plus a model of the instrumental noise. We apply the procedure to SGWB signals with a variety of representative frequency profiles, and prove that LISA can reconstruct their spectral shape. Our procedure, implemented in the code SGWBinner, is suitable for homogeneous and isotropic SGWBs detectable at LISA, and it is also expected to work for other gravitational wave observatories.

Classifying inflationary scenarios according to their scaling properties is a powerful way to connect theory with observations. A useful tool to make such a classification is the beta-function formalism. By describing inflation in terms of renormalization group equations, within this framework, it is possible to define universality classes, which can be considered as sets of theories that share a common scale invariant limit. In this paper we apply the formalism to define such classes of universality for models of inflation where the inflaton is coupled to gauge fields. We show that the formalism may consistently be extended to capture the peculiar features of these models such as statistical anisotropy. We also obtain some consistency conditions which serve as useful guidelines for model building.

We forecast the prospective of detection for a stochastic gravitational wave background sourced by cosmological first-order phase transitions. We focus on first-order phase transitions with negligible plasma effects, and consider the experimental infrastructures built by the end of the LISA mission. We make manifest the synergy among LISA, pulsar time array experiments, and ground-based interferometers. For phase transitions above the TeV scale or below the electroweak scale, LISA can detect the corresponding gravitational wave signal together with Einstein Telescope, SKA or even aLIGO-aVIRGO-KAGRA. For phase transitions at the electroweak scale, instead, LISA can be the only experiment observing the gravitational wave signal. In case of detection, by using a parameter reconstruction method that we anticipate in this work, we show that LISA on its own has the potential to determine when the phase transition occurs and, consequently, the energy scale above which the standard model of particle physics needs to be modified. The result may likely guide the collider community in the post-LHC era.

In this paper the dynamics of free gauge fields in Bianchi type I-VII$_{h}$ space-times is investigated. The general equations for a matter sector consisting of a $p$-form field strength ($p\,\in\,\{1,3\}$), a cosmological constant ($4$-form) and perfect fluid in Bianchi type I-VII$_{h}$ space-times are computed using the orthonormal frame method. The number of independent components of a $p$-form in all Bianchi types I-IX are derived and, by means of the dynamical systems approach, the behaviour of such fields in Bianchi type I and V are studied. Both a local and a global analysis are performed and strong global results regarding the general behaviour are obtained. New self-similar cosmological solutions appear both in Bianchi type I and Bianchi type V, in particular, a one-parameter family of self-similar solutions,"Wonderland ($\lambda$)" appears generally in type V and in type I for $\lambda=0$. Depending on the value of the equation of state parameter other new stable solutions are also found ("The Rope" and "The Edge") containing a purely spatial field strength that rotates relative to the co-moving inertial tetrad. Using monotone functions, global results are given and the conditions under which exact solutions are (global) attractors are found.

We develop a scenario of inflation with spontaneously broken time and space diffeomorphisms, with distinctive features for the primordial tensor modes. Inflationary tensor fluctuations are non adiabatic, and can acquire a mass during the inflationary epoch. They can evade the Higuchi bound around de Sitter space, thanks to interactions with the fields driving expansion. Correspondingly, the primordial stochastic gravitational wave background (SGWB) is characterised by a tuneable scale dependence, and can be detectable at interferometer scales. In this set-up, tensor non-Gaussianity can be parametrically enhanced in the squeezed limit. This induces a coupling between long and short tensor modes, leading to a specific quadrupolar anisotropy in the primordial SGWB spectrum, which can be used to build estimators for tensor non-Gaussianity. We analyse how our inflationary system can be tested with interferometers, also discussing how an interferometer can be sensitive to a primordial anisotropic SGWB.

In this paper we investigate the past evolution of an anisotropic Bianchi I universe in $R+R^2$ gravity. Using the dynamical system approach we show that there exists a new two-parameters set of solutions that includes both an isotropic "false radiation" solution and an anisotropic generalized Kasner solution, which is stable. We derive the analytic behaviour of the shear from a specific property of $f(R)$ gravity and the analytic asymptotic form of the Ricci scalar when approaching the initial singularity. Finally we numerically checked our results.

Supersolid inflation is a class of inflationary theories that simultaneously breaks time and space reparameterization invariance during inflation, with distinctive features for the dynamics of cosmological fluctuations. We investigate concrete realizations of such a scenario, including non-minimal couplings between gravity and the fields driving inflation. We focus in particular on the dynamics of primordial gravitational waves and discuss how their properties depend on the pattern of symmetry breaking that we consider. Tensor modes can have a blue spectrum, and for the first time we build models in which the squeezed limit of primordial tensor bispectra can be parametrically enhanced with respect to standard single-field scenarios. At leading order in a perturbative expansion, the tensor-to-scalar ratio depends only on the parameter controlling the breaking of space-reparameterization. It is independent from the quantities controlling the breaking of time-reparameterization, and this represents a difference with respect to standard single-field inflationary models.

We investigate the potential for the LISA space-based interferometer to detect the stochastic gravitational wave background produced from different mechanisms during inflation. Focusing on well-motivated scenarios, we study the resulting contributions from particle production during inflation, inflationary spectator fields with varying speed of sound, effective field theories of inflation with specific patterns of symmetry breaking and models leading to the formation of primordial black holes. The projected sensitivities of LISA are used in a model-independent way for various detector designs and configurations. We demonstrate that LISA is able to probe these well-motivated inflationary scenarios beyond the irreducible vacuum tensor modes expected from any inflationary background.

Primordial Gravitational Waves are the next target of modern cosmology. They represent a window on the early Universe and the only probe of the physics and microphysics of the inflationary period. When the production of GWs happens in scenarios richer than the standard single-field slow-roll, the GW signal becomes potentially detectable also on scales smaller than the Cosmic Microwave Background. LISA will be extremely complementary to CMB experiments to extract information about primordial inflationary models and in particular to probe phases of the inflationary period for which we have very poor knowledges.

Despite its continued observational successes, there is a persistent (and growing) interest in extending cosmology beyond the standard model, $\Lambda$CDM. This is motivated by a range of apparently serious theoretical issues, involving such questions as the cosmological constant problem, the particle nature of dark matter, the validity of general relativity on large scales, the existence of anomalies in the CMB and on small scales, and the predictivity and testability of the inflationary paradigm. In this paper, we summarize the current status of $\Lambda$CDM as a physical theory, and review investigations into possible alternatives along a number of different lines, with a particular focus on highlighting the most promising directions. While the fundamental problems are proving reluctant to yield, the study of alternative cosmologies has led to considerable progress, with much more to come if hopes about forthcoming high-precision observations and new theoretical ideas are fulfilled.

The effective field theory of inflation is a powerful tool for obtaining model independent predictions common to large classes of inflationary models. It requires only information about the symmetries broken during the inflationary era, and on the number and nature of fields that drive inflation. In this paper, we consider the case for scenarios that simultaneously break time reparameterization and spatial diffeomorphisms during inflation. We examine how to analyse such systems using an effective field theory approach, and we discuss several observational consequences for the statistics of scalar and tensor modes. For example, examining the three point functions, we show that this symmetry breaking pattern can lead to an enhanced amplitude for the squeezed bispectra, and to a distinctive angular dependence between their three wavevectors. We also discuss how our results indicate prospects for constraining the level of spatial diffeomorphism breaking during inflation.

Solid inflation is an effective field theory of inflation in which isotropy and homogeneity are accomplished via a specific combination of anisotropic sources (three scalar fields that individually break isotropy). This results in specific observational signatures that are not found in standard models of inflation: a non-trivial angular dependence for the squeezed bispectrum, and a possibly long period of anisotropic inflation (to drive inflation, the "solid" must be very insensitive to any deformation, and thus background anisotropies are very slowly erased). In this paper we compute the expected level of statistical anisotropy in the power spectrum of the curvature perturbations of this model. To do so, we account for the classical background values of the three scalar fields that are generated on large (superhorizon) scales during inflation via a random walk sum, as the perturbation modes leave the horizon. Such an anisotropy is unavoidably generated, even starting from perfectly isotropic classical initial conditions. The expected level of anisotropy is related to the duration of inflation and to the amplitude of the squeezed bispectrum. If this amplitude is close to its current observational limit (so that one of the most interesting predictions of the model can be observed in the near future), we find that a level of statistical anisotropy $\gtrsim 3\%$ in the power spectrum is to be expected, if inflation lasted $\gtrsim 20-30$ e-folds more than the final $50-60$ efolds required to generare the CMB modes. We also comment and point out various similarities between solid inflation and models of inflation where a suitable coupling of the inflaton to a vector kinetic term $F^{2}$ gives frozen and scale invariant vector perturbations on superhorizon scales.

Primordial statistical anisotropy is a key indicator to investigate early Universe models and has been probed by the cosmic microwave background (CMB) anisotropies. In this paper, we examine tensor-mode CMB fluctuations generated from anisotropic gravitational waves, parametrised by $P_h({\bf k}) = P_h^{(0)}(k) [ 1 + \sum_{LM} f_L(k) g_{LM} Y_{LM} (\hat{\bf k}) ]$, where $P_h^{(0)}(k)$ is the usual scale-invariant power spectrum. Such anisotropic tensor fluctuations may arise from an inflationary model with noncommutativity of fields. It is verified that in this model, an isotropic component and a quadrupole asymmetry with $f_0(k) = f_2(k) \propto k^{-2}$ are created and hence highly red-tilted off-diagonal components arise in the CMB power spectra, namely $\ell_2 = \ell_1 \pm 2$ in $TT$, $TE$, $EE$ and $BB$, and $\ell_2 = \ell_1 \pm 1$ in $TB$ and $EB$. We find that B-mode polarisation is more sensitive to such signals than temperature and E-mode polarisation due to the smallness of large-scale cosmic variance and we can potentially measure $g_{00} = 30$ and $g_{2M} = 58$ at 68% CL in a cosmic-variance-limited experiment. Such a level of signal may be measured in a PRISM like experiment, while the instrumental noise contaminates it in the $Planck$ experiment. These results imply that it is impossible to measure the noncommutative parameter if it is small enough for the perturbative treatment to be valid. Our formalism and methodology for dealing with the CMB tensor statistical anisotropy are general and straightforwardly applicable to other early Universe models.

We investigate parity-violating signatures of temperature and polarization bispectra of the cosmic microwave background (CMB) in an inflationary model where a rolling pseudoscalar produces large equilateral tensor non-Gaussianity. By a concrete computation based on full-sky formalism, it is shown that resultant CMB bispectra have nonzero signals in both parity-even $(\ell_1 + \ell_2 + \ell_3 = {\rm even})$ and parity-odd $(\ell_1 + \ell_2 + \ell_3 = {\rm odd})$ spaces, and are almost uncorrelated with usual scalar-mode equilateral bispectra. These characteristic signatures and polarization information help to detect such tensor non-Gaussianity. Use of both temperature and E-mode bispectra potentially improves of $400\%$ the detectability with respect to an analysis with temperature bispectrum alone. Considering B-mode bispectrum, the signal-to-noise ratio may be able to increase by 3 orders of magnitude. We present the $1\sigma$ uncertainties of a parameter depending on a coupling constant and a rolling condition for the pseudoscalar expected in the ${\it Planck}$ and the proposed PRISM experiments.

In the model of solid / elastic inflation, inflation is driven by a source that has the field theoretical description of a solid. To allow for prolonged slow roll inflation, the solid needs to be extremely insensitive to the spatial expansion. We point out that, because of this property, the solid is also rather inefficient in erasing anisotropic deformations of the geometry. This allows for a prolonged inflationary anisotropic solution, providing the first example with standard gravity and scalar fields only which evades the conditions of the so called cosmic no-hair conjecture. We compute the curvature perturbations on the anisotropic solution, and the corresponding phenomenological bound on the anisotropy. Finally, we discuss the analogy between this model and the f (phi) F^2 model, which also allows for anisotropic inflation thanks to a suitable coupling between the inflaton phi and a vector field. We remark that the bispectrum of the curvature perturbations in solid inflation is enhanced in the squeezed limit and presents a nontrivial angular dependence, as had previously been found for the f (phi) F^2 model.

A suitable coupling of the inflaton phi to a vector kinetic term F^2 gives frozen and scale invariant vector perturbations. We compute the cosmological perturbations zeta that result from such coupling by taking into account the classical vector field that unavoidably gets generated at large scales during inflation. This generically results in a too anisotropic power spectrum of zeta. Specifically, the anisotropy exceeds the 1% level (10% level) if inflation lasted ~5 e-folds (~50 e-folds) more than the minimal amount required to produce the CMB modes. This conclusion applies, among others, to the application of this mechanism for magnetogenesis, for anisotropic inflation, and for the generation of anisotropic perturbations at the end of inflation through a waterfall field coupled to the vector (in this case, the unavoidable contribution that we obtain is effective all throughout inflation, and it is independent of the waterfall field). For a tuned duration of inflation, a 1% (10%) anisotropy in the power spectrum corresponds to an anisotropic bispectrum which is enhanced like the local one in the squeezed limit, and with an effective local f_{NL} ~3 (~30). More in general, a significant anisotropy of the perturbations may be a natural outcome of all models that sustain higher than 0 spin fields during inflation.