We investigate the possibility of measuring the primordial gravitational wave (GW) signal across 23 decades in frequencies, using the cosmic microwave background (CMB), pulsar timing arrays (PTA), and direct detection with laser and atomic interferometers. For the CMB and PTA experiments we consider the LiteBIRD mission and the Square Kilometer Array (SKA), respectively. For the interferometers we consider space mission proposals including the Laser Interferometer Space Antenna (LISA), the Big Bang Observer (BBO), the Deci-hertz Interferometer Gravitational wave Observatory (DECIGO), the $\mu$Ares experiment, the Decihertz Observatory (DO), and the Atomic Experiment for Dark Matter and Gravity Exploration in Space (AEDGE), as well as the ground-based Einstein Telescope (ET) and Cosmic Explorer (CE) proposals. We implement the mathematics needed to compute sensitivities for both CMB and interferometers, and derive the response functions for the latter from the first principles. We also evaluate the effect of the astrophysical foreground contamination in each experiment. We present binned sensitivity curves and error bars on the energy density parameter, $\Omega_{GW}h^2$, as a function of frequency for two representative classes of models for the stochastic background of primordial GW: the quantum vacuum fluctuation in the metric from single-field slow-roll inflation, and the source-induced tensor perturbation from the spectator axion-SU(2) inflation models. We find excellent prospects for joint measurements of the GW spectrum by CMB and space-borne direct detection mission proposals.

We study how the shape of the spectrum of primordial gravitational waves can be constrained by future experiments looking at the B-mode of the Cosmic Microwave Background (CMB) polarization. We implement a Principal Component Analysis (PCA) including the effects of diffuse foreground residuals, following component separation, in the uncertainty of CMB angular power spectra, and taking into account the gravitational lensing by Large Scale Structure. We perform our study by considering the capabilities of future B-mode CMB experiments such as LiteBIRD, the Simons Observatory (SO) and Stage-IV (CMB-S4), in particular in detecting deviations of the primordial tensor spectrum from the scale-invariant behavior. We find that diffuse foreground residuals impact substantially both the derivation of the PCA basis and the corresponding constraining power, in all cases. In particular, depending on which experimental specifications and which value $r$ of tensor-to-scalar ratio for cosmological perturbations are considered, adding foregrounds residuals can determine an increase as large as a factor $\sim 4$ both on the uncertainty on $r$ and on the recovery of the PCA modes. We study the limitations of the methodology, including the effect of physicality priors on the PCA, which we quantify via a Monte Carlo Markov chain (MCMC) analysis of the combined cosmological and PCA power spectrum parameter space.