Gravitational Waves: the theorist's Swiss knife

Gravitational waves provide a novel and powerful way to test astrophysical models of compact objects, early universe processes, beyond the Standard Model particle physics, dark matter candidates, Einstein's theory of General Relativity and extended gravity models, and even quantum gravity candidate theories. A short introduction to the gravitational-wave background and the method we are using to detect it will be presented. Constraints on various astrophysical/cosmological models from the non-detectability of the gravitational-wave background will be discussed. Gravitational waves from transients will be highlighted and their physical implications will be summarised.

The LIGO/Virgo interferometers have so far detected several tens of phase-coherent transient signals emitted by compact binaries (black holes, neutron stars, or combination of both).In addition, there is an ongoing search for phasecoherent persistent signals, like the continuous waves emitting by spinning neutron stars [31][32][33], and phase-incoherent transient or persistent signals; the latter one being the long-sought stochastic gravitational-wave background.The techniques used to search for GW signals depend on their nature and it is important to note that even in the case of a non-detection (like for the gravitational-wave background), still important constraints can be imposed to various theoretical models and their free parameters.
In what follows, I will briefly summarise the current status of constraints imposed to various astrophysical and cosmological models, as well to fundamental theories using the novel and powerful tool of gravitational waves.

II. GRAVITATIONAL-WAVE BACKGROUND
Our universe is permeated by cosmic microwave background electromagnetic radiation, the relic photons of the Big Bang, first discovered by Penzias and Wilson in 1965.Similarly, our universe is expected to be filled with a stochastic gravitational-wave background generated by various processes that took place in its early history.In addition, a background of gravitational waves can also emerge from the incoherent superposition of a large number of astrophysical sources, which are too weak to be detected separately, and such that the number of sources that contribute to each frequency bin is much larger than one.These two components (the stochastic cosmological and the astrophysical one) form the gravitational-wave background (GWB).
The GWB [34] can be expressed in terms of a dimensionless parameter (like all cosmological density parameters), which expresses the energy density of gravitational waves ρ gw in a logarithmic interval of frequency f , normalised by the critical energy density ρ cr = 3H 2 0 c 2 /(8πG), as with ρ gw ∼ ḣ2 , where h denotes the strain of GWs and overdot stands for derivative with respect to cosmological time t.
The GWB is characterised by random non-deterministic phase evolution of GWs emitted from a large number of distant sources.On the one hand, the GWB of cosmological origin looses any phase coherence present at emission or at horizon crossing, as metric perturbations affect the phase of GWs propagating through the inhomogeneous universe.On the other hand, the GWB of astrophysical origin, generated by overlapping signals from many individual sources below confusion limit of detector, has no phase-coherence.For a standard cosmological model, it has been shown [35] that there is complete phase decoherence of the GWB above frequencies f ∼ 10 −12 Hz.Hence, only the amplitude of GWs plays a rôle in estimating Ω gw : where h denotes the strain in Fourier space and A stands for the GW amplitude.
A detection of the GWB from unresolved compact binary coalescences (CBC) could be made by Advanced LIGO and Advanced Virgo at their design sensitivities.It would appear as noise in a single GW detector i: with the noise n being much bigger than the strain h at a given frequency f .Thus, to detect the GWB, we need to take the correlation between two detector outputs i, j: The signal-to-noise (SNR) ratio grows slowly over observing time T obs : Assume that the GWB is isotropic, Gaussian, stationary and unpolarised: where H 0 stands for the Hubble constant and A, A ′ denote the +, x polarisation modes.We build the estimator Ĉij (f ; t) for the GWB measured from detectors i, j: where si (f ; t) denotes the Fourier transform of the strain time series in detector i starting at time t, Γ ij (f ) is the normalised overlap reduction function between detectors i andj accounting for all the geometric factors that come into play when cross-correlating data from different detectors, T stands for the duration over which the Fourier transform is taken, and S 0 (f ) = 3H 2 0 /(10π 2 f 3 ) is the spectral shape for a GWB that is flat in energy density.Since the GWB is considered to be isotropic, Gaussian, stationary and unpolarised, it is well-described by a single power spectral density S gw (f ): where δ T (f − f ′ ) denotes the finite-time approximation to the dirac delta function and S gw (f ) is related to Ω gw (f ) through Assuming the GW signal and the intrinsic noise to be uncorrelated and the noise in each frequency bin to be independent, we obtain In the absence of correlated noise, ñ⋆ i (f )ñ j (f ) = 0, the Ĉij (f ) is an estimator of the dimensioneless gravitational wave energy density Ω gw (f ).However, it may be that ñ⋆ i (f )ñ j (f ) = 0; a potential source of correlated magnetic noise is Schumann resonances.
Schumann resonances are global extremely low frequency peaks in the electromagnetic field of the Earth.They are eigenmodes of the conducting spherical cavity formed by the surface of the Earth and its ionosphere, and are excited by lightning discharges.These low-frequency resonances in the Earth's global electromagnetic field could couple to the mass suspension system and electronics in the interferometers, and mimic a GWB signal.One must therefore develop a method for detecting correlated magnetic noise and separating from a true GWB signal.This issue is usually addressed by performing a coherent subtraction in the raw data using Wiener filtering.However, for weakly coupled signals it can be difficult to completely subtract the noise.Recently, an alternative novel method to minimise the possibility of a false detection was proposed and tested with realistic simulations of correlated magnetic noise in interferometric GW detectors in [36].
The correlated magnetic noise in GW detectors was modelled using data collected by magnetometers placed near the detector sites, and adopting a parameterised model for the magnetic field to GW detector coupling [36].Peaks in the detectors' strain channels were related to the peaks in the magnetometer channels via a coupling function, assumed to be constant in time, well-described by a power-law, with index β.Hence, with Ω M,ij (f ) representing the magnetic contribution: κ denotes the amplitude of the coupling at 10 Hz at the site of interferometric GW detector.The factor 10 −22 assumes that the unites of magnetometer data mi (f ) are T/Hz.Bayesian model selection was then used to distinguish between models that include correlated magnetic noise and those with a true GWB.This method originally proposed [36] and used in the framework of LIGO/Virgo detectors [37], was consequently extended in the case of the 3g Einstein Telescope interferometer [38] forecasting the necessary measures to ensure that magnetic contamination will not pose a threat to the corresponding science goals.In particular, it was shown [38] that for GWB searches below ∼ 30 Hz it will be necessary for the Einstein Telescope magnetic isolation coupling to be two to four orders of magnitude better than the one measured in the current Advanced LIGO and Virgo interferometers.
In the Advanced LIGO's and Adavanced Virgo's third observing run (O3) data from all three baselines were included in the isotropic GWB search.Combining O3 data with upper limits from the earlier first and second (O1 and O2) observing runs, it was found that the results are consistent with uncorrelated noise, hence placing upper limits on the strength of the GWB [37].More precisely, performing a Bayesian analysis that allows for the presence of both a GWB and an effective magnetic background arising from Schumann resonances, no evidence of correlated noise of magnetic origin was found.
Considering a power-law GWB of the form where f ref is a fixed reference frequency chosen to be f ref = 25 Hz and α is the spectral index describing the source of the GWB, limits on the dimensionless energy density in GWs are imposed for α = 0, 2/3, 3. The α = 0 approximately characterises a GWB originated by a network of cosmic strings, or a period of slow-roll inflation in the LIGO-Virgo frequency band.The α = 2/3 power-law describes the GWB from compact binaries coalescences when contributions from the inspiral phase dominate the GWB.This can be easily seen by writing down the quadrupole expression for the energy E gw released per frequency from the inspiral phase of two compact objects of masses m 1 , m 2 : and then integrating over all frequencies f .While this consideration is a very good approximation in the LIGO-Virgo frequency band, it is not be valid for mergers of binaries arising from population III stars [8].The α = 3 corresponds to a GWB that is flat in the strain power, S h (f ) ∝ f −3 Ω gw (f ) and describes a GWB from supernovae.
Detection of a GWB is indeed hard, as one can easily observe by comparing the dimensionless GW energy parameter to the corresponding of the cosmic microwave background (Ω cmb ): Ω gw ≪ Ω cmb ≈ 10 −5 .The challenge of claiming a true detection of the GWB will be immediately succeeded by the difficulty of relating the signal to the sources (astrophysical or cosmological) that contribute to it.This task was addressed in [39], considering that compact binary coalescences and two additional cosmological sources, namely cosmic strings (CS) and early universe first order phase transition (FOPT), comprise the GWB.It was shown [39] that the network of Advanced LIGO and Advanced Virgo, operating at design sensitivity, will not be able to allow for separation of these sources.However, 3g detectors (such as Einstein Telescope and Cosmic Explorer), which are sensitive to most individual compact binary mergers, may reduce the astrophysical signal via subtraction of individual sources, and potentially reveal a cosmological background.Performing a Bayesian analysis, it was found [39] that considering an Einstein Telescope and Cosmic Explorer network and reasonable levels of individual source subtraction, will allow us to detect a cosmological signal of Ω cs (25Hz) ≃ 4.5 × 10 −13 for cosmic strings and Ω bpl (25Hz) ≃ 2.2 × 10 −13 for a broken power-law model describing GWs produced during an early universe phase transition at f ref = 25 Hz.

A. GWB from CBC: information on compact binaries
Consider an astrophysical GWB composed by compact binaries of population I and II (pop I/II) stars.The energy density of GWs from compact binaries as a function of frequency f and model parameters θ can be expressed as where p(θ) is the probability distribution of the source parameters, dE gw (f s ; θ)/df s is the energy density emitted by as single source at a redshift z with parameters θ, E(Ω m , Ω Λ , z) = Ω m (1 + z) 3 + Ω Λ accounts for the expansion history of the universe with Ω m = 0.31 and Ω Λ = 0.69 the matter and dark energy contribution to the dimensioneless cosmological energy density, respectively, f s = (1 + z)f is the emitted frequency in the source frame, and z max (θ) is the maximal redshift at which a compact binary with parameters θ can form.
The most important quantities describing binary black holes ate the masses and spins of each black hole component.The impact on uncertainties in the distribution of binary black holes on the astrophysical GWB has been investigated in [40] using population inference from GW observations.Adopting a truncated power-law black hole mass distribution and modelling the black hole spin magnitudes by a Beta distribution, about 10,000 models have been considered, showing that the astrophysical GWB monopole is sensitive to the nature of the binary black hole population, and in particular the local merger rate.Hence, the GWB can be used as a new probe of the population of compact objects throughout cosmic history, complementary to direct observations.
The theoretically postulated population III (pop III) stars are thought to have formed at high redshifts, hence they have low metallicity as compared to pop I/II stars which are more recently formed.If pop III exists, their GWB signal will be superposed with a pop I/II signal.Several models indicate that the total CBC background is dominated by pop I/II, implying that the pop III signal will remain hidden underneath the pop I/II signal.For second generation detector networks (even by including LIGO-Hanford, LIGO-Livingston, Virgo, LIGO-India, and KAGRA) the pop III signal is lost in the pop I/II foreground and hence invisible.However, third generation terrestrial detectors, such as Einstein Telescope and Cosmic Explorer, may reveal a pop III background [8].
A CBC signal is typically modelled with a 2/3 spectrum.This is appropriate for pop I/II stars which lead to a GWB in the LIGO/Virgo frequency star at the inspiral phase.However, the older pop III stars would show up at the later merger and ringdown phases, hence there will be a clear deviation from the 2/3 spectrum.Testing filters different from a 2/3 power spectrum and performing a Bayesian analysis, it was shown [8] that a broken power-law filter is the most appropriate one to reveal a pop III GWB.
Considering a detector network composed by Einstein Telescope and two Cosmic Explorers, it was shown that one could successfully subtract enough individual merger events to unravel the elusive pop III GWB.The reason being that subtraction methods are more effective for nearby sources, so the residual pop I/II signal can become sub-dominant to the pop III residual.It was also found that a detection of a pop III GWB could reveal important information, such as the mass and redshift distribution of pop III stars [8].

B. Stochastic GWB from CS: information on particle physics models beyond the Standard Model
Cosmic strings are one-dimensional topological defects, which may have been formed in the early universe as a result of a phase transitions followed by spontaneously broken symmetries, characterised by a vacuum manifold with non-contractible closed curves.These defects have been shown to be generically formed in the context of Grand Unified Theories [41].
A cosmic string network reaches an attractor scaling solution in which its energy density remains a fixed fraction of the background energy density [42][43][44].As strings collide, they always intercommute exchanging partners and reconnect after a collision.As a result, closed loops are formed when a string self-intersects or two curved strings collide.Loops smaller than the horizon decouple from the cosmological evolution and oscillate under their own tension, slowly decaying into GWs [45].
The relativistic nature of strings leads to the formation of points where the string momentarily moves at the speed of light; these points are called cusps.The intersections of strings generate discontinuities on their tangent vecto; they are called kinks.Cusps and kinks generate gravitational wave bursts [46].Cusps produce beamed GWs in the forward direction of the cusp, kinks propagate around the string creating GWs with a fan-like emission, while the collision of two kinks radiates GWs isotropically.The high-frequency gravitational-wave spectrum of an oscillating loop is dominated by bursts emitted by cusps and kinks.The waveforms of the strain are power-law functions in the frequency domain, and grow as f −4/3 , f −5/3 f −2 for cusp, kink, and kink-kink collision cases, respectively.In addition to these distinct individual bursts, the incoherent superposition of weaker GW bursts from cosmic strings produced over the history of the universe leads to a strongly non-Gaussian stochastic GWB.
The cosmic string network also contains long strings that stretch across a Hubble volume, which are either infinite or in the form of super-horizon loops; they also emit GWs [47].However, the dominant GW contribution is the one produced by the superposition of GW radiation from sub-horizon loops along each line of sight.
The spectrum of the GWB from cosmic strings is where d 2 R i /dz/dℓ is the burst rate of type i -cusp, kink or kink-kink collision -per unit loop size ℓ and per redshift z; it depends on the number of cusps and kinks.Consider Nambu-Goto strings, parametrised by a single quantity, namely the string tension Gµ (with G Newton's constant and µ the string linear mass density; note that we set c = 1), related to the energy scale η at string formation by Gµ ∼ (η/M Pl ) 2 (with M Pl the Planck mass).To estimate the GWB generated by cosmic strings one must know the loop distribution and the average number of cusps and kinks.Numerical simulations have shown that while kinks accumulate over the cosmological evolution, the number of cusps per loop remains yet undetermined.
There are two different models [48,49] for the loop distribution based on different numerical simulations.Their differences arise from different modelling of the production and cascade of loops from the infinite string network.In addition, a new agnostic model that extends and encompasses the ones based on the numerical simulations, has been recently proposed and studied [50].
A search for a stochastic GWB was conducted using the LIGO and Virgo O3 data and combining them with those from the previous O1 and O2 runs.Performing a Bayesian analysis taking into account the precise shape of the GWB from cosmic strings, upper limits were derived on the cosmic string tension [15].More precisely, for model [49] the stochastic analysis ruled out particle physics models leading to cosmic strings with tension Gµ (9.6 × 10 −9 − 10 −6 ).For this model, the strongest limit comes from pulsar timing measurements, excluding string tension Gµ 10 −10 .For model [48], the stochastic analysis ruled out particle physics models with string tension Gµ (4.0 − 6.3) × 10 −15 ; this is the strongest limit in the literature.Two different variants of model [50] were considered.One that reproduces qualitatively the loop production function of model [49] in the radiation era and the loop production of model [48] in the matter era, for which string tension Gµ (2.1 − 4.5) × 10 −15 was ruled out aside form a small region where the number of kinks is greater than 180.For this variant, the exclusion range obtained from the stochastic LIGO/Virgo analysis is the strongest one in the literature.The other variant reproduces qualitatively the loop production function of model [48] in the radiation era and the loop production of model [49] in the matter era, for which particle physics models with string tension Gµ (4.2 − 7.0) × 10 −15 were excluded.The corresponding constraint obtained from the cosmic microwave background is almost as strong as this one.We note that the average number of cusps per oscillation on a cosmic string loop has been set to 1.A high number of cusps leads qualitatively to a similar result as increasing the number of kinks: for model [49], the constraints are weakened, whereas for all other considered models the bounds are insensitive to the number of cusps [15].
Considering the cosmic string models [48,49] for the loop distribution, it was shown [51] that LISA will be able to probe cosmic strings with tension Gµ O(10 −17 ), improving by about 6 orders of magnitude current pulsar timing arrays (PTA) constraints, and potentially 3 orders of magnitude with respect to expected constraints from next generation PTA observatories.However, to obtain these limits only the LISA noise was taken into account.Considering the cosmic string GWB in the presence of a compact binary produced astrophysical background and a galactic foreground, in addition to the LISA noise, it was shown [52] that with four years of data, LISA will be able to measure a cosmic string tension Gµ ≈ 10 −16 (for the model [53]) to Gµ ≈ 10 −15 (for the model [44,48,49]) or bigger; the galactic foreground affecting the Gµ limit more than the astrophysical background.

C. Stochastic GWB from FOPT: information on particle physics beyond the Standard Model
As the universe expands and its temperature drops, a series of phase transitions might have taken place.In the case of a first-order phase transition (FOPT), once the temperature drops below a critical value, the universe, being originally in a meta-stable phase, reaches a stable one, through a sequence of bubble nucleation, growth, and merger.During this process, a stochastic GWB is thought to be generated [54,55] by mainly three sources: bubble collisions, sound waves, and magnetohydrodynamic turbulence.Several extensions of the Standard Model predict strong FOPTs; their footprint on the stochastic GWB may thus reveal important information on beyond the Standard Model particle physics constraining their free parameters.
The stochastic GWB produced by a FOPT covers a large frequency range, with the peak frequency set by the temperature T FOPT at which the FOPT takes place.For a T FOPT ∼ (10 7 − 10 10 ) GeV, which is the case for the Peccei-Quinn axion model or high-scale supersymmetry models, the generated stochastic GWB lies within the frequency range of Advanced LIGO and Advanced Virgo.Thus, current terrestrial interferometers have the ability to constrain particle physics models at energy scales which remain inaccessible by any existing terrestrial accelerator.
Using the absence of a stochastic GWB signal in the data from the last three observing runs (O1-O3) of the LIGO/Virgo, one can constrain Ω gw sourced by FOPT, and consequently test the underlying particle physics models.The O1-O3 data were analysed [16] using either an approximated broken power-law (bpl) describing its main features in terms of the shape and the peak amplitude, or a phenomenological model accounting for contributions from bubble collisions and sound waves.In the former approach, the stochastic GWB spectrum is approximated by with n 1 = 3 from causality, n 2 equal to -4 or -1 for sound wave and for bubble collisions, respectively, and ∆ equal to 2 or 4, for sound waves and approximating bubble collisions, respectively.One then runs a Bayesian search and model selection considering both a broken power law background and a CBC background approximated by a f 2/3 power law.In the second approach, it is adopted a phenomenological model with the dominant contribution coming either from bubble collisions (bc), following or from sound waves (sw), following where h is the dimensionless Hubble parameter, ∆ stands for the amplitude, H FOPT is the Hubble parameter at T FOPT , β is the inverse time duration of the transition, κ φ denotes the fraction of vacuum energy converted into gradient energy of the scalar field, α is the energy released normalised by the radiation energy density, g ⋆ is the number of relativistic degrees of freedom chosen to be 100, S env is the spectral shape, κ sw is the fraction of vacuum energy converted into the kinetic energy of the bulk flow, v w is the bubble wall velocity, f sw denotes the present peak frequency, and Υ is a suppression factor due to the finite lifetime τ SW .The analysis following the generic power-law spectrum, led to 95% confidence level upper limits simultaneously on the GW energy density at 25 Hz from unresolved CBC and strong FOPT, Ω CBC < 6.1 × 10 −9 and Ω bpl < 4.4 × 10 −9 , respectively [16] .The analysis following the phenomenological model, constrained at 25 Hz the stochastic GWB to Ω bc < 5.0 × 10 −9 and Ω sw < 5.8 × 10 −9 at 95% confidence level due to bubble collisions and sound waves, respectively, [16] D. Stochastic GWB: information on the early universe Several mechanisms in the early universe can create parity violation manifesting itself in the production of asymmetric amounts of right-and left-handed circularly polarised isotropic GWs.Since astrophysical sources of the GWB are unlikely to create circular polarisation, a detection of a circularly polarised signal may allow us to distinguish the cosmologically sourced component of the GWB from the astrophysically sourced one.Analysing polarised GWB can then be used to constrain parity violating theories.
Right-and left-handed correlators can then be written as where • represents the ensemble average and I, V are the Stokes parameters, with V characterising the asymmetry between right-and left-handed polarised waves, and I(≥ |V |) the wave's total amplitude.Note that for V = 0, the above equation would just be the correlator for unpolarised isotropic GWB.The cross-correlation estimator reads [34] where with H 0 the Hubble parameter, T the observing time, δ T (f ) = sin(πf T )/(πf ), sd1 (f ) and sd2 (f ) the Fourier transforms of the strain time series of GW detectors d 1 and d 2 , Q(f ) the optimal filter taking into account detectors' strain power spectral densities, and F A n = e A ab d ab n the contraction of the tensor modes of polarisation A to the n th detector's geometry.We denote by γ d1d2 I the standard overlap reduction function of detectors d 1 , d 2 , and by γ d1d2 V the overlap function associated with the parity violation term.The polarisation degree, Π(f ) = V (f )/I(f ), can take any value between -1 (fully left polarisation) and 1 (fully right polarisation), with Π = 0 denoting an unpolarised isotropic GWB.
A strategy based on parameter estimation and subsequent fit of two GW models to O3 data, as well as simulated GW data (for A+ LIGO-Virgo detectors including KAGRA) was proposed in [13].The proposed method, using a hybrid frequentist-Bayesian approach, was employed in the search of a generic parity violation GW signal with a power-law behaviour, as well as a broken power-law spectral shape motivated by high energy physics.The search has shown no evidence for a parity violation signal in the most recent O3 data, while an upper limit of Ω(f ref = 25hz) = 4.9 × 10 −9 at 95% confidence limit was imposed on the amplitude of a power-law GW model.Simulating GW data for future detection prospects, and considering a chiral turbulence source in the early universe it was demonstrated the importance of having a multi-detector network, while it was shown that even if we detect a turbulence signal, we may not be able to deduce its polarisation [13].
Subsequently, it was investigated [14] the ability of current or next generation detectors to detect parity violation sources from axion inflation focusing on the quadratic potential.It was shown [14] that using two 3g detectors -Einstein Telescope (ET) and Cosmic Explorer (CE) -one may be able to claim with at least 95% confidence a positive polarisation degree Π, while a network of three 3g detectors (ET and 2 CEs) is needed in order to make a confident claim about the detection of a quadratic axion inflation signature.

E. Stochastic GWB: information on dark matter (primordial black holes)
Primordial black holes in the mass range (30 − 100)M ⊙ are interesting candidates for dark matter, as they occupy a narrow window between microlensing and cosmic microwave background constraints.By revisiting the primordial black holes merger rate to account for a more adequate description of black holes in their surrounding spacetime, it was shown that the merger rate constraints on the abundance of dark matter PBHs, imposed by LIGO, are evaded entirely [56,57].
Several studies in the literature examine whether some of the LIGO-Virgo detected CBC events are the outcome of PBH coalescences (see, for instance [58,59]).One of the mechanisms leading to PBH formation in the early universe is the large curvature fluctuations generated during an inflationary era.More precisely, PBHs result from the gravitational collapse of overdensed regions seeded by large curvature perturbations on small scales once the corresponding wavelength re-enters the horizon as the curvature perturbations exceed a critical value.Since the curvature perturbations couple to the tensor perturbations at second-order, such large curvature perturbations are accompanied by the scalar induced GWs in the radiation dominated era [60].
The spectrum of the scalar induced GWB scales as the square of the curvature power spectrum P ζ (k), which at large scales is O(10 −9 ) [61], set by cosmic microwave background data, and hence the scalar induced GWB remains unobservable by current or planned experiments.However, P ζ (k) may be considerable larger at small scales, reaching P ζ (k) ≈ O(0.01), and hence resulting to a GWB within reach of gravitational-wave observatories.
A Bayesian search of scalar induces GWB in the data from Advanced LIGO and Virgo's first, second and third observing runs was performed in [12].Avoiding the choice of a particular inflationary model, the curvature power spectrum was parametrised by a well-localised peak depending on only three parameters: its position, width and integrated power, using a log-normal function [62] where A ζ is the integrated power of the peak, ∆ ζ determines its width and k * its position.The scalar induced GW spectrum is peaked around the same wavenumber as P ζ (k), corresponding to the frequency f * /Hz = 1.6×10 −15 k * /Mpc −1 .At ∆ ζ ≪ 1, for which the amplitude of the induced GWs as well as the generated PBH abundance become independent of ∆ ζ , the peak amplitude of the GWB is Ω gw = O(10 −5 )A 2 ζ .Since the LIGO-Virgo detectors, are sensitive to frequencies around f 10 Hz, they have the potential to probe peaks in P ζ (k), at scales larger than 10 15 Mpc −1 , which re-enter the horizon at temperatures 10 8 Gev.
Performing a Bayesian analysis -including also potential contributions from CBC -no evidence was found for a scalar induced GWB in the LIGO-Virgo O1, O2 and O3 data [12].However, depending on the width ∆ ζ of the peak, 95% confidence level upper limits were imposed on the integrated power of the peak A ζ , as a function of its position k ⋆ .The obtained constraints are stronger than the ones arising form big bang nucleosynthesis and cosmic microwave background observations in the range O(10 15 ) < k * /Mpc −1 < O(10 18 ), and reach A ζ ≃ 0.02 for a narrow peak at k * ≃ 10 17 Mpc −1 .It is worth noting the constraints from LIGO and Virgo, at their design sensitivity, and from the Einstein Telescope can compete with those related to the abundance of the formed primordial black holes [12].

F. Stochastic GWB: Tests of General Relativity
Einstein's theory of General Relativity allows for the existence of only two gravitational-wave polarisation modes (the tensor plus and cross modes), while general metric theories of gravity may allow for up to four additional polarisations: two vector modes and two scalar modes (the breathing and longitudinal scalar modes).The detection of any of these additional polarisation modes would imply violation of General Relativity, while their non-detection may experimentally constrain extended theories of gravity.However, the three-detector Advanced LIGO-Virgo network is generally unable to distinguish the polarisation of transient gravitational-wave signals, like those from binary black holes.The twin LIGO detectors (Livingston and Hanford) are nearly co-oriented, implying that Advanced LIGO is sensitive to only a single polarisation mode.Moreover, even if the twin LIGO detectors were more favorablyoriented, one needs generically a network of at least six detectors to uniquely determine the polarisation content of a gravitational-wave transient.
Observation of the stochastic GWB offers a way to directly measure generic gravitational-wave polarisations.It has been shown that a stochastic GWB dominated by alternative polarisation modes may be missed by current searches optimised only for tensor polarisations.A Bayesian method to detect a generically-polarised stochastic GWB was proposed in [63] and subsequently used in the LIGO/Virgo observing runs.This (generic) method relies on the construction of two odds ratios: one to determine if a stochastic GWB has been detected, another one to quantify evidence for the presence of alternative polarisations in the background.Employing this method, it was shown [63] that flat scalar-polarised backgrounds of amplitude Ω (S) (25 Hz) ≈ 2 × 10 −8 can be confidently identified as non-tensorial with Advanced LIGO.Note that Ω gw does not necessarily have the interpretation of an energy density in modified theories of gravity; one may think of it as a measure of the strain power in each polarisation.It was also shown [63] that after three years of observation at design sensitivity, Advanced LIGO will be able to limit the amplitudes of tensor, vector, and polarisations to Ω (T) gw (25 Hz) < 1.6 × 10 −10 , Ω (V) gw (25 Hz) < 2.0 × 10 −10 , and Ω (S) gw (25 Hz) < 5.0 × 10 −10 , respectively, at 95% credibility.In the case that a stochastic GWB of mixed polarisation is detected, Advanced LIGO alone will not be able to precisely determine the parameters of the tensor, vector, and/or scalar components simultaneously due to large degeneracies between modes.It is worth noting that while addition of Advanced Virgo does not particularly increase our ability to detect or identify backgrounds of alternative polarisations, it does significantly improve our parameter estimation on power-law backgrounds, hence breaking the degeneracies of the twin LIGO analysis.
So far no evidence for scalar or vector modes was found in the data from Advanced LIGO and Advanced Virgo.The non-detection of scalar and vector polarised stochastic GWB is consistent with predictions of Einstein's theory of General Relativity.Using the log-uniform prior for Ω ref and a Gaussian prior on the spectral index for each polarisation, the upper limits Ω (S) (25 Hz) ≤ 2.1 × 10 −8 , Ω (V) (25 Hz) ≤ 7.9 × 10 −9 , Ω (S) (25 Hz) ≤ 6.4 × 10 −9 , at the 95% credible level on Ω ref for scalar, vector, and tensor polarisations, respectively, were obtained [37].
The GW data analysis of the LVK Collaboration is first done in low-latency so that public alerts can be generated that subsequently trigger follow-up astronomical observations, and then in higher latency in order to create GW catalogues.To analyse GW data from the last (O3) observing run, five pipelines were used in real time: an unmodelled burst search, and four matched-filter pipelines.The matched-filter method relies on having a model of the signal, as a function of the intrinsic (two individual component masses and two dimensionless spin vectors) and extrinsic (luminosity distance, right ascension and declination, inclination between total angular momentum and line-of-sight, time of merger, reference phase, and polarisation angle) physical parameters of the source.Using the search pipelines, a template bank of GW waveforms is created, covering the appropriate intrinsic parameter space; it is consequently used to filter against the data and produce signal-to-noise ratio time series.

A. GW transients: tests of General Relativity
Gravitational waves emitted by CBCs allow us to test GR in extreme environments of strong gravitational fields, large spacetime curvature, and velocities close to the speed of light.Considering the most significant signals, namely the ones with a false-alarm rate (FAR) < 10 −3 yr −1 , from the gravitational-wave transient catalogue GWTC02 [3], observational constraints of General Relativity (GR) were reported [22] by the LVK Collaboration, following a twofold approach.First, the consistency of GR predictions with the data was tested and then, by parametrising modifications to GW waveforms, the amount of allowed deviations from GR predictions was constrained.
To investigate whether the GR waveforms are successful in describing the data, one studies the residual strain, once the best-fit template for each event has been subtracted.Such an analysis has not found any statistically significant deviations between the observed residual power and the detector noise [22].Within GR the final state of the coalescence of two black holes will be a single perturbed Kerr black hole, gradually relaxing to its Kerr stationary state by emitting GWs which correspond to a specific set of characteristic quasinormal modes.The explicit analysis of the GWTC02 events, has shown [22] that all GW transient signals were produced by merging Kerr black holes in agreement with Einstein's theory of General Relativity.
In addition, one can use the GWTC02 events to constrain parameters of phenomenological modifications to General Relativity.In the context of GR, GWs far from their source propagate along null geodesics, with energy E and momentum p related by the dispersion relation E 2 = p 2 c 2 , with c denoting the speed of light.To investigate modifications to General Relativity, a phenomenological approach was considered introducing a modified dispersion relation [64] E 2 = p 2 c 2 + A α p α c α , with α = 0, 0.5, 1.0, 1.5, 2.5, 3.0, 3.5, 4.0 , with A α and α two phenomenological parameters.Such an expression holds for a variety of proposed extensions to General Relativity.For instance, α = 0, A α > 0 represents massive gravity with graviton mass m g = A 1/2 0 c −2 ; note that General Relativity is recovered for A α , ∀α.All cases apart α = 0 correspond to a Lorentz-violating dispersion relation.
In frequency-domain, the strain in Fourier space, within General Relativity, is A non-zero A α will lead to a frequency-dependent dephasing of the GW signal, δΦ α (f ), building up as the GW propagates towards Earth.For a given model (i.e., given the values of A α , α) the dephasing δΦ α (f ) depends on the binary's luminosity distance, the binary's detector-frame chirp mass, and the effective wavelength parameter used in the sampling, defined in terms of binary's redshift, and a distance parameter for a given cosmological model.Using the gravitational-wave transient catalogue GWTC02 [3] 90% credible upper bounds on the absolute value of the modified dispersion relation parameter A α as a function of α where obtained in [22].The analysis has shown no evidence for GW dispersion, constraining the Lorentz-violating dispersion parameters.The obtained constraint on the graviton mass is [22] with 90% credibility; it is an improvement of 1.8 over Solar System bounds.Since we have only 3 active detectors, we cannot simultaneously constrain all six possible GW polarisations allowed in generic metric theories of gravity.For GW transients, one can then compare the likelihood of having purely scalar or purely vector polarisations against the pure tensor case, predicted by General Relativity.Such an approach, using multiple binary black hole signals reported in the gravitational-wave transient catalogue GWTC02 [3], has shown [22] no evidence in favor of non-tensor GWs.

B. GW transients: tests of modified gravity models
Several modified gravity phenomenological models have been proposed in the literature, in an attempt to find an explanation for the dark energy.Gravitational waves offer the means to test/constrain such modified gravity models by studying the propagation of GWs across cosmological distances [24,25].
In the context of General Relativity, gravitational waves travelling on a four-dimensional Friedmann-Lemaître-Roberson-Walker background, obey the linearised evolution equation where A = +, × stands for the two polarisation plus and cross modes, primes denote derivatives with respect to conformal time η, related to the cosmological time through dη = dt/a(t) with a(t) the scale factor, H is the Hubble parameter in conformal time η, and Π A denotes the source term related to the anisotropic stress tensor.The GW propagation equation above, gets modified in a generic modified gravity model into where three new quantities have been introduced.The function δ(η) modifies the friction term and hence affects the amplitude of a GW propagating across cosmological distances.The tensor velocity c T can be in general time and scale dependent; in General Relativity it is equal to the speed of light c.The mass of the tensor mode m T , can be non-zero in the context of a modified gravity theory.These three quantities are in principle testable with GW data.The modification in the tensor sector leads to the gravitational-wave luminosity distance d where with Ω m , Ω r and Ω DE the present matter density, radiation density and DE fractions, respectively.A simple phenomenological parametrisation which depends on the (positive) parameters Ξ 0 and n (with Ξ 0 = 1 in General Relativity) has been proposed in [65].This approach has been employed in [24] to test several modified gravity models -Horndeski or the more general degenerate higher order scalar-tensor theories, non-local infrared modifications of gravity, bi-gravity theories, as well as theories with extra or varying dimensions -using standard sirens at GW detectors such as LISA or third-generation ground based experiments.More precisely, constructing mock source catalogues and performing Markov Chain Monte Carlo studies of the likelihood obtained from LISA standard sirens combined with other cosmological datasets, it was found [24] that the parameter Ξ 0 can be measured to the percent level accuracy.
In conclusion, it has been shown that the modified GW propagation is an important observable for advanced GW detectors; supermassibe black hole mergers binaries detectable with LISA can provide a powerful probe of modified gravity and dark energy.
C. GW transients: information on dark matter (axion-like particles) Pseudo-scalar (axion-like) particles are hypothetical scalar particles that appear in many fundamental theories.String theory compactifications generally predict a number of light axions; they represent one of the most compelling candidates for dark matter.
Consider axions that couple to nuclear matter in the same way as the QCD axion, but with masses that are relatively light.Such axions can be sourced by compact objects with a high nuclear density, for instance neutron stars.If the axion decay constant is f a 10 18 , neutron stars with radius of the order of 10 km will develop an axion profile, interpolating from ±πf a near the neutron star surface to 0 at spatial infinity.For two neutron stars, with each one within the axion's Compton wavelength, λ a ≡ 1/m a , the axion field mediates an additional force, which at leading order is Q 1,2 denote the axion charge carried by each neutron star, related to its radius R 1,2 through The axion force F a is attractive or repulsive, depending on whether the axion field values have the same or opposite sign, respectively, on the surfaces of the two neutron stars.Compact objects with a scalar charge may have important implications for binary neutron stars or neutron star-black hole mergers, leading to the emission of axion radiation and an axion-mediated fifth force, which effects the orbital dynamics leaving potentially detectable fingerprints in the inspiral GW waveform [17].More precisely, if such neuron stars form a binary system, the axion field might also radiate axion waves during the compact objects coalescence.Considering circular orbits, the radiation power at leading order reads where ω denotes the orbital frequency and r the separation between the two neutron stars of masses M 1 , M 2 .The axion radiation turns on once the orbital frequency exceeds the axion mass.
Using an effective field theory approach, the first post-Newtonian corrections to the orbital dynamics, radiated power, and gravitational waveform for binary neutron star mergers in the presence of an axion have been calculated in [18].It was shown that Advanced LIGO at design sensitivity, can potentially exclude axions wuth mass m a 10 −11 eV and decay constant f a ∼ (10 14 − 10 17 ) GeV.
Subsequently, a search for such axions has been conducted [19] by analysing the gravitational waves from the binary neutron star inspiral GW170817.Let us write the frequency domain waveform as [19] h with where Ψ GR is the phase of the waveform within General Relativity, Ψ a denotes the leading order phase correction due to the axion field, and v 2 counts the post-Newtonian order.The effects of the axion field in the waveform template can be parametrised by the axion mass m a and a dimensionless parameter γ a ≡ Q 1 Q 2 /(4πGM 1 M 2 ) that characterises the relative strength of the axion and gravitational forces between the two neutron stars.Performing a Bayesian analysis of GW17081, the posterior distribution over the waveform parameters suggests no significant evidence for such axion fields [19].The null result constrainted axions with masses below 10 −11 eV by excluding the ones with decay constants ranging from 1.6 × 10 16 GeV to 10 18 GeV at a 3σ confidence level.This is the first constraint on axions imposed from neutron star inspirals [19].
It is worth noting that such constraint may be improved if the SNR of the merger event is enhanced; the constraint on f a will improve by √ N if the SNR is improved by N .The constraint could also be improved by roughly 2 orders of magnitude by considering a neutron star-black hole merger.

D. GW transients: information on quantum gravity theories
Gravitational waves can also provide means to test quantum gravity (QG) theories, challenging one's naive claim that any late-time QG imprint will be Planck-suppressed and hence undetectable.The reason being the existence of a longrange non-perturbative mechanism, called dimensional flow (scale-dependence of the dimensionality of spacetime), shared by most QG candidate theories.
Let us consider propagation of gravitational waves in the context of QG theories.Quantisation of spacetime geometry introduces two types of change relevant for the propagation of GWs: an anomalous spacetime measure and a kinetic operator.The former captures the way that volumes scale and the latter is related to the modified dispersion relations.Given a spacetime metric, a kinetic operator and a compact source, the GW amplitude, in radial coordiantes and in the local wave zone, reads [27,28] where Γ ≃ d H /2 − d k H /d s , with d H , d k H and d s the Hausdorff dimension, the Hausdorff dimension in momentum space, and the spectral dimension, respectively.Clearly, h is the product of a dimensionless function f h , depending on the source and the type of correlation function, and a power-law distance behaviour.This is a rather general feature in quantum gravity theories.
The amplitude of GWs propagating in such a background is [27,28] h ∝ 1 with ǫ = O(1), γ = 1, and ℓ ⋆ a characteristic scale of the geometry.If there is only one fundamental length scale ℓ ⋆ = O(ℓ Pl ), the above equation relating the GW luminosity distance to the electromagnetic one, is exact and γ = Γ UV , namely scales at which QG corrections are important.If however ℓ ⋆ is a mesoscopic scale Γ meso , denoting intermediate scales where the corrections are small but non-negligible, then the above equation is valid only near the infrared regime and γ = Γ meso ≈ 1.The values of Γ UV and Γ meso are known for large variety of QG candidate theories (including group field theory, spin foams, loop quantum gravity, causal dynamical triangulations, κ-Minkowski, low-energy limit of string theory, asymptotic safety).The left-hand side of Eq. ( 37) is the strain measured in a GW interferometer, while the right-hand side gives the luminosity distance measured for the optical counterpart of the standard siren.GW observations can thus be used to constraint ℓ ⋆ and γ, by constraining the ratio d (z) as a function of the redshift z of the source.Considering two examples of multimessenger standard sirens, namely the binary neutron-star merger GW170817 and a simulated supermassive black hole merger event detectable with LISA, it was found [27] that only group field theory, spin foams or loop quantum gravity could generate a signal detectable with standard sirens.

IV. CONCLUSIONS
Gravitational waves give us the means to unravel mysteries of our universe, allowing us to test theories describing its early stages, when it was opaque to electromagnetic radiation.Using direct observations of gravitational waves we can now test particle physics models at energy scales which are not accessible to any foreseeable particle accelerator, test Einstein's theory of General Relativity in the strong gravity regime, and also constrain modified/extended gravity models.Gravitational waves offer us an independent way to study the large scale structure of our universe, and to obtain a useful insight on astrophysical properties of compact objects.Finally, they allow us to constrain some of the most promising dark matter candidates, while they open the way to test quantum gravity theories.