1. Introduction
Tracing the evolution of the Universe back to its origins is one of the hardest tasks in physics, involving general relativity and quantum physics in regimes that cannot be reached by any conceivable experiments. At the same time, the early Universe perhaps provides the only window offered by nature to explore such regimes and look for possible observable consequences that may shed light on the foundations of the pillars on which our description of the world is based [
1,
2,
3]. Hubble’s law is our main observational evidence for the expansion of the Universe. If space-time participates in the dynamics of the matter it contains, it must be described through a theory that treats it as a dynamical entity. General relativity provides a beautiful geometrical framework describing the evolution of space-time and its relation to the energy-momentum tensor. Taking the present expansion of the Universe as a matter of fact, the generic prediction of general relativity for a universe filled with ordinary matter and radiation is an initial singularity called the Big Bang. Such a simple model of the Universe suffers from some well-known problems [
4]:
The horizon problem: the Universe appears homogeneous on scales that only now are coming back in causal contact.
The flatness problem: the Universe has little or no spatial curvature, requiring extremely special initial conditions.
A possible solution to these problems is a phase of accelerated expansion called inflation [
5,
6,
7], which would stretch the space between particles to distances larger than the cosmological horizon, and also make the Universe spatially flat. At the end of the inflation, the Universe is essentially empty, save for quantum fluctuations that should seed the cosmological perturbations that later evolve into large-scale structures harboring galaxies and clusters of galaxies [
8,
9].
In general, inflationary models typically assume that the accelerated expansion takes place at some energy scale below the Planck scale (e.g., at the grand unification [
7]), thus representing just a special phase within a Universe with ever-decreasing energy density and space-time curvature. In this case, any physics at energy scales higher than the inflation plays little or no role, as it is washed out by the following accelerated expansion. Therefore, any quantum gravity effects at the Planck scale and even the mere existence of the initial singularity remain hidden from our observations [
10,
11].
It is widely believed that the Big Bang singularity is an accident of the extrapolation of general relativity beyond its range of validity [
12]. When the curvature radius of the Universe becomes comparable with the Planck length, quantum effects should dominate the structure of space-time [
13]. What happens at this scale depends on the unknown physics of quantum gravity. In this respect, string theory has played a major role in inspiring theoretical physicists about the possibilities opened by a quantum theory encompassing general relativity [
14,
15]. Furthermore, the opportunity of finding any possible observational signatures of string theory through relics of an early stringy phase of the Universe has stimulated the research in many different directions [
3]. One of the possibilities explored by string cosmologists in the last four decades is that the initial singularity might be cured by string theory and replaced by a finite maximum in the curvature and energy density. Then we can extend back the history of the Universe past this maximum in a new phase known as the Pre-Big Bang [
16,
17,
18]. In this scenario, proposed by Gasperini and Veneziano in the 1990s, the Universe started from an asymptotically flat space-time, where quantum fluctuations randomly grow and nucleate bubbles undergo a super-inflation phase (an accelerated expansion with growing curvature) [
19]. When the curvature of these bubbles reaches the string scale (slightly below the Planck scale), string physics stops the super-inflation and gives rise to a decelerated expansion converging to the standard cosmological picture, including radiation, a matter-dominated phase, or any late-time evolution.
Similar to Pre-Big Bang cosmology, many other string cosmologies were proposed thereafter in which the Big Bang is replaced by a cosmic bounce, i.e., a transition between a contraction era and the present expansion [
20,
21,
22] (for a review see [
2]). Considering that the Pre-Big Bang superinflation becomes a contraction after a transformation of the action from the string frame to the Jordan frame, bouncing cosmologies can be viewed as a generalization of the Pre-Big Bang scenario, which was particularly motivated by the driving concept of T-duality [
16].
The Pre-Big Bang and bouncing cosmologies solve the problems of standard cosmology similarly to standard inflation [
23]. However, important differences arise when we look at cosmological perturbations. While standard inflation predicts a nearly scale-invariant spectrum both for scalar and tensor perturbations [
4], Pre-Big Bang cosmology (and other related string cosmologies) predicts steep blue spectra [
24,
25]. As soon as the first observations of the cosmic microwave background (CMB) in the first decade of the millennium revealed a nearly scale-invariant spectrum for scalar perturbations [
26], it became clear that a revision of the mechanism for the generation of cosmological perturbations was necessary for Pre-Big Bang cosmology.
In this contribution, in honor of the 70th birthday of Maurizio Gasperini, it is my pleasure to revive the studies of cosmological perturbations in bouncing cosmologies of the early years of the millennium in which I was personally involved. The legacy of those studies is still strong and continues to inspire theories and observations looking for stringy signatures in our sky. In
Section 2 I present the Pre-Big Bang scenario and other bouncing cosmologies. In
Section 3 I discuss the primordial spectra of perturbations generated in bouncing cosmologies. In
Section 4 we face the delicate issue of how these perturbations evolve through the bounce. In
Section 5 we will see how a scale-invariant spectrum can be re-generated after the bounce and check the observational and theoretical constraints. Finally, we draw some conclusions in
Section 6.
2. The Pre-Big Bang Phase
Pre-Big Bang cosmology assumes that at sufficiently low energies space-time is fully described by the effective action of the bosonic sector of string theory [
14]:
where
is the fundamental string length,
d is the number of spatial dimensions,
is the metric tensor,
R is the Ricci scalar,
is the scalar dilaton field, coming with a non-perturbative potential
, and
is the field strength of the Kalb–Ramond antisymmetric field
.
String theory is consistently formulated with
or
(in the M-theory version). We can thus imagine that the Pre-Big Bang phase leads to an expansion of only three spatial dimensions accompanied by a contraction of the remaining six dimensions (if we adopt
), whose volume shrinks down to the string scale. Therefore, we split the metric in the following form
where
is the conformal time,
(
) are the coordinates along the three large dimensions,
(
) are the coordinates in the internal small dimensions.
is the scale factor for the large dimensions and
is the scale factor describing the contraction of the small dimensions. We have assumed homogeneity and isotropicity within each of the two blocks of coordinates. In this framework, we may introduce an effective 4-dimensional dilaton
and write the antisymmetric field strength in terms of a single pseudoscalar axion field
as
[
18].
The dilaton is non-minimally coupled to the metric in this action. However, by a conformal transformation of the metric, we can go from the physical string frame to an equivalent Einstein frame in which the dilaton is minimally coupled. The new metric is related to the old one by
In the Einstein frame, the description of the Pre-Big Bang phase is particularly simple, since it becomes an accelerated contraction driven by the kinetic energy of the dilaton field. The exact dynamics, however, depends on the evolution of the spatial dimensions. The equations for the metric, dilaton and axion are
where
and the prime denotes derivative with respect to the conformal time.
If the axion does not contribute to the cosmic background in the Pre-Big Bang phase and the dilaton has a vanishing potential, then we have the following solutions
Since the dilaton and the moduli field enter as degenerate scalar fields, the contraction rate of the internal dimensions s remains as a free parameter. The solution for the dilaton depends on this parameter and may be growing or decaying depending on the chosen sign in the exponent.
These solutions are singular for and may describe either a Pre-Big Bang () or a Post-Big Bang (). It is believed that the full string theory removes the singularity and ensures that the Pre-Big Bang contraction bounces to a standard Post-Big Bang evolution. It is also assumed that the same string phase is characterized by a non-trivial potential for the dilaton, which is frozen in the Post-Big Bang so that the effective gravitational coupling remains fixed to the observed Planck mass .
The Pre-Big Bang phase, then, appears as an accelerated contraction governed by Equation (
9) in the Einstein frame. From the point of view of the string frame, instead, we have
which corresponds to an accelerated expansion (pole inflation) in the negative sign branch if
. Such superinflation is able to solve the problems of standard cosmology and may represent a good alternative to slow-roll inflation if we only consider the background evolution. In the Einstein frame, the exponent of the scale factor is fixed to
, whatever the relative contraction rate of the internal dimensions. Alternative bouncing cosmologies have been proposed with different contraction rates, depending on the field content and the geometry of the space-time [
21,
22]. In the next section, we will leave this exponent as a free parameter so as to derive the primordial spectra of cosmological perturbations for a general class of bouncing cosmologies. We will therefore set
in the pre-bounce phase and go back to
for the Pre-Big Bang. We note that the space-time curvature vanishes at
. Therefore, all these models start from an asymptotically flat past [
19].
It has been noticed that fast contractions with
are exposed to uncontrolled growth of anisotropies that may lead to chaotic mixmaster oscillations [
27]. The Pre-Big Bang scenario is just at the divide between safe and unsafe backgrounds. A possible solution to this problem was proposed in Ref. [
28].
4. Across the Bounce
The evolution of tensor perturbations in bouncing cosmologies is quite straightforward since nothing is expected to happen to modes outside the horizon at the bounce. The spectrum generated at the onset of the bounce should pass unaffected and thus a steep blue spectrum is expected in all models, including the Pre-Big Bang scenario [
24]. This prediction marks a sharp distinction with respect to standard inflationary models, where the tensor spectrum is always nearly scale-invariant or just slightly red. So, if any contribution of tensor modes to low multipoles of the CMB is ever measured, bouncing cosmologies would be ruled out or need to be completely reformulated. On the other hand, the blue spectrum opens the possibility to measure a stochastic gravitational-wave background at high frequencies, accessible to future interferometers [
48].
Scalar perturbations are more tricky since they directly depend on the energy-momentum source dominating the Universe. In order to step out of the initial collapse and avoid a big crunch singularity, some new physics is invoked to intervene at the quantum gravity scale. In the Pre-Big Bang scenario, for example, the string mass plays the role of a cut-off for the curvature and the energy density [
18]. The string effective action should be complemented by higher-order terms that should drive the Universe toward an expansion phase. As an interesting alternative, an explicit realization of the bounce through a non-local potential for the dilaton has been presented in Refs. [
49,
50].
Non-singular bouncing models have been obtained in many ways by adding a cosmological component violating the null energy condition at the bounce [
51,
52,
53,
54,
55,
56,
57]. In the ekpyrotic/cyclic model, the collision of our brane with the hidden brane is believed to occur through a singularity of the metric, in the sense that the distance between the two branes vanishes at the bounce [
21].
In any case, one may describe the evolution of the Universe during the bounce by complementing the Einstein equations with additional terms arising in this high-energy regime. These terms can be moved to the right-hand side and act as an effective additional source. In this respect, the violation of the null energy condition (NEC), which is needed to convert a collapse phase to an expansion, can be performed by these additional effective terms. Therefore, the generation of possible ghosts and instabilities due to negative energy density can be avoided if the NEC violation is due to these higher-order terms. Explicit realizations of ghost-free bouncing scenarios have indeed been presented [
49,
50,
56].
We may guess that the evolution of scalar perturbations across the bounce may be heavily affected by these additional terms, which may appear as effective sources in the perturbation Equations (
27)–(
29). However, this is not necessarily the case. In fact, inhomogeneities in the sources dominating the bounce necessarily descend from the inhomogeneities generated all along the previous pre-bounce collapse. Leaving apart isocurvature perturbations generated by subdominant fields in the pre-bounce, if we wish to study the evolution of the adiabatic mode, we assume that the source terms in Equations (
28) and (
29) are functions of
and
with the addition of possible powers of
from Laplacian operators:
with
F,
G,
J, and
K being regular power expansions in
[
31].
At this point, the evolution of the gauge-invariant variables
and
is fully contained in the closed set of Equations (
28) and (
29). Although we do not know the details of the background evolution during the bounce, we may identify the modes surviving after the bounce and their time dependence by putting the equations in the integral form and solving them recursively, as proposed in Ref. [
31]. The only two scales in the problem are the wave number
k and the fundamental bounce scale
, which can be identified as the scale at which the new physics triggering the bounce comes into play. This scale governs the bounce duration, energy density, and curvature. So, each integral over the conformal time
in the formal solutions of Equations (
28) and (
29) will introduce a new factor of
in the expressions without changing the
k dependence.
With these rules, it is then possible to write down the post-bounce behavior of the scalar perturbations as [
31]
where
are constants possibly related to the fundamental bounce scale
. We have also parameterized the post-bounce scale-factor as
∼
, in analogy with Equation (
14). The first two modes coincide with the dominant modes generated in the pre-bounce. However, since the contraction has been converted to an expansion, the growing mode of
is now a decaying mode [
38]. The final spectrum will be thus dominated by the constant modes. Leaving apart the last term for the moment, the dominant contribution comes from the third mode, which gives rise to the spectral index
already discussed in the pre-bounce context as yielding a scale-invariant spectrum in the limit of a dust-like contraction [
22]. The last terms contains a factor
, remarking the fact that it only exists if the
function introduced in Equation (
39) is non-zero for
[
31]. Therefore, only if during the bounce there is a dependence
∼
rather than
∼
, the spectrum of the growing mode of
is transferred to a constant mode and survives after the bounce.
This condition was already stated in a different way by modeling the bounce as a thin space-like hypersurface and applying Israel junction conditions [
42]. It clarifies a long-debated issue about the viability of ekpyrotic/cyclic models as an alternative to standard inflation for the generation of the primordial spectrum of scalar perturbations. A slow contraction indeed generates a scale-invariant spectrum that matches a decaying mode in the post-bounce, unless we have some unconventional source proportional to
rather than
in the spatial Einstein equations. No explicit example of such sources has been provided up to now. Regular bounces obtained by perfect fluids or scalar fields only involve sources with
∼
, which is bound to be proportional to
by Equation (
27). Therefore, these toy models have confirmed that the original growing mode of
decays in the post-bounce expansion [
52,
53,
54,
57]. However, it might still be possible that some new physics mechanisms may replace the
factor with another scale in the problem, being the bounce scale
itself or some geometric scale, such as the size of the extra-dimensions or similar. Therefore, it would still be interesting to continue the search for fully self-consistent bouncing cosmologies in which the same fields dominate the pre-bounce background and generate the observed spectra. Otherwise, we have no other route than looking at the perturbations of sub-dominant fields.
5. Re-Generation of Cosmological Perturbations by the Axion/Curvaton
Any credible alternative to standard inflation should contain a mechanism to generate fluctuations compatible with CMB and large-scale structure observations. This requires a nearly scale-invariant adiabatic primordial spectrum of scalar perturbations with the correct amplitude.
As a result of the studies outlined in the previous sections, we have seen that the Pre-Big Bang cosmology and other related bouncing cosmologies inspired by string theory or other quantum gravity theories are not able to naturally generate a nearly scale-invariant spectrum for scalar perturbations if we confine our attention to the dominant fields. However, we have also noticed that the axion field of the Pre-Big Bang scenario, although subdominant, typically develops a perturbation spectrum with a slope spanning a relatively wide range of possibilities depending on the specific background evolution [
45]. In particular, a scale-invariant spectrum is obtained for a particularly symmetric background in which the contraction rates of internal and external dimensions are the same. However, even in the case that such fluctuations dominate in some regimes, they would only give rise to isocurvature fluctuations, which are severely constrained by current observations.
A possible mechanism to convert isocurvature perturbations to adiabatic ones has been suggested through the so-called curvaton field [
58]. The idea is that the field responsible for the standard inflation (or the collapse/inflation for bouncing cosmologies) is not the same that generates the observed fluctuations. A second scalar field takes over as the dominant field at some point in cosmic history and then imprints its fluctuations on
and
. The resulting perturbations would be adiabatic as they come from the new dominant source. At the end of its lifetime, the curvaton decays to radiation fields, and the standard expansion history proceeds without other changes.
As soon as the curvaton mechanism was proposed, it was clear that the Pre-Big Bang axion would make an ideal candidate to implement the curvaton in a practical case [
59]. A scalar field driven by its kinetic energy redshifts faster than radiation and would never dominate in an expanding Universe. For a scalar field with a potential, we may consider two different regimes: if the potential dominates and the kinetic energy is negligible, the scalar field behaves as a cosmological constant and gives rise to an exponential expansion; if the scalar field oscillates at the bottom of its potential, the average expansion rate tends to be similar to a matter-dominated phase with
a∼
.
In the Pre-Big Bang scenario, it is believed that the dilaton field is stabilized by a non-perturbative potential acquired during the string phase. Similarly, the axion may receive a non-perturbative potential and find itself displaced from the minimum. So, as the expansion of the Universe starts with an early radiation phase, the axion slowly rolls down its potential. At some point, it will dominate over radiation. Depending on its initial value, this may happen when it is still in a slow-roll phase or when it has already reached the minimum of its potential and started its oscillations. It is in this phase that the axion fluctuations dominate and become adiabatic. As the axion finally decays, its fluctuations are naturally inherited by the new radiation field, which will dominate the following era.
Although the route for the curvaton/axion mechanism is clear, many details need to be worked out explicitly to check that no physical constraints are violated [
60,
61].
Let us assume that the initial value of the axion at the onset of the Post-Big Bang expansion is
. At this time, the dilaton is already frozen at the present value and the extra-dimensions are finally stabilized, so that the space-time dynamics is effectively governed by the Friedmann equations. Moreover, with a fixed dilaton, there is no distinction between Einstein and string frames, so
a∼
up to a constant. We will thus use the symbol
a for the scale factor without the tilde in the Post-Big Bang. The energy-momentum tensor is initially dominated by a radiation field, while the axion is still subdominant. However, a non-perturbative potential
arises for the axion, so that the background evolution is described by the following set
with
and
.
The potential for the Kalb–Ramond axion arising after the string phase is believed to be periodic. However, close to a minimum, it makes sense to deal with a quadratic approximation . Then the axion starts from some initial value such that . If the axion kinetic energy is initially negligible, we have a slow roll phase for the axion. As the axion reaches the minimum of its potential, it starts to oscillate. This occurs when the expansion rate in cosmic time ∼m. During oscillations, the energy density of the axion redshifts as , similarly to pressureless matter.
Either as a slowly rolling scalar field or as pressureless matter, the axion will sooner or later dominate over the early radiation, which redshifts as . The axion dominance will occur at . Depending on the initial value of the axion , the axion phase may start either before or after the first oscillation. In particular, we have an early dominance of the axion if (in Planck units) and a late dominance otherwise. The two cases have different background evolutions and must be treated separately.
Finally, the axion will decay through its gravitational coupling to photons when the space-time curvature becomes of the order of the decay rate at
∼
. After this phase, the standard radiation era takes place as usual and will inherit the fluctuation spectrum of the axion as its adiabatic mode. The successful transfer of the initial spectrum of the axion to the re-generated radiation has been tested analytically and numerically in Ref. [
61].
The space allowed for an axion-dominated phase in the early Universe is however limited by several requirements and observations that need to be checked carefully. We are going to list them here [
61].
First of all, the decay of the axion must occur well before the primordial nucleosynthesis or even baryogenesis. Depending on the specific baryogenesis mechanism, the lower bound may become very strict or somewhat relaxed. Here we adopt the nucleosynthesis bound, which translates into a minimal mass for the axion of
Secondly, the axion should not decay before its dominance phase. This means that we need
. This constraint is automatically satisfied for the case of the early dominance
. For the late dominance, we have
∼
and the constraint requires
On the high-energy side, we do not want the axion to dominate at the beginning of the Post-Big bang phase: it was a subdominant field in the Pre-Big Bang and thus its energy density should be somewhere below the string scale. Therefore,
Figure 1 illustrates the space left by these constraints on the background evolution. The requirement of an early radiation phase before the axion dominance depends on the bounce scale
, but still, there is a very wide range of possible masses allowed for the Kalb–Ramond axion and for its initial value [
61].
Besides the requirements on the background evolution, we now have very accurate measurements from several cosmological observations: galaxy counts, weak lensing, and in particular the CMB anisotropies, which provide the most stringent limits on the amplitude and the slope of the power spectrum. We have already discussed the spectral index of the axion fluctuations, which is transferred to and via the curvaton mechanism.
The amplitude of the axion fluctuations is typically regulated by the ratio of the string scale (coinciding with
in our Post-Big Bang history) to the Planck mass. The adiabatic spectrum inherited by the curvature perturbations can be further amplified if the phase preceding the axion dominance is very long. In fact, from Equation (
28) we see that any non-adiabatic pressure leads to a growth of
∼
. In the early radiation phase, this corresponds to
∼
∼
. Therefore, the lower the initial axion value, the higher the amplification taking place before the axion dominates. However, for large values of the axion, the axion-driven inflation amplifies the fluctuations proportionally to
. The interplay between these two effects determines the final amplitude of the adiabatic scalar perturbations encoded in the Bardeen potential as [
61]
where
and the constants
,
,
can be estimated by fitting the spectra obtained by careful numerical simulations.
Comparing the amplitude of this spectrum with the current PLANCK limits [
62] of
we can put interesting constraints on the bounce scale.
Figure 2 shows that a bounce scale
generates too large fluctuations, incompatible with the observations. This value is what is naturally expected for the string mass and would be the natural choice for a bounce driven by string effects modifying general relativity at these scales. However, starting from
, an amplitude of scalar fluctuations compatible with the PLANCK observations becomes possible.
The mechanism driving the bounce may well start being effective at scales slightly below the string mass. Indeed we do not know enough about string physics to exclude or validate the possibility that the bounce occurs at scales of the order of M. Note that this scale is quite close to the GUT scale. Indeed, the observed amplitude of cosmological perturbation indicates that the mechanism for their generation should have something to do with this scale, either through standard inflation or through a cosmic bounce.
6. Conclusions
Standard inflation provides a very simple and predictive model for the early Universe. It solves the curvature and horizon problems and naturally provides a nearly scale-invariant spectrum compatible with the observations. There is currently no compelling need to replace inflation with other models. However, everybody knows that inflation cannot be the end of the story. In particular, what happens to our space-time if we go past beyond the inflationary era remains a theoretical puzzle. Indeed inflation makes the Universe strongly independent of what existed before the accelerated expansion and any possible weirdness in space-time. However, we cannot exclude that any signatures of a quantum Universe survived through inflation and are observable today. Therefore, the investigation of the initial singularity problem and the possible alternatives offered by string theory and other quantum gravity alternatives remain important now as ever. Furthermore, although the current cosmological picture is overall very robust, some minor inconsistencies and tensions might actually be serious signals of new physics that need to be addressed. It is possible that some of these details hide the sought signatures of an ancient pre-inflationary era.
Bouncing cosmologies start from the assumption that the Big Bang singularity may just be the outcome of a hazardous extrapolation of general relativity beyond its domain of validity. A complete theory should avoid the singularity and allow us to extend our past history back to a pre-bounce growing curvature phase. Among the possible bouncing cosmologies, here we have focused on the Pre-Big Bang scenario, which uses the fields available in the bosonic sector of the heterotic string theory to build a complete history of the Universe, from an asymptotically flat past to the current expansion phase.
The main issue of bouncing cosmologies is that they generally predict steep blue power spectra for scalar perturbations, in contrast with all observations of the cosmic microwave background, galaxy counts, and large-scale structures. However, we have seen at least three ways to obtain a scale-invariant spectrum.
A dust-dominated contraction would generate a scale-invariant spectrum in the pre-bounce era, but it is possibly exposed to background instabilities.
A slow contraction would generate a scale-invariant spectrum in the growing mode, which is matched to a decaying mode in the Post-Big Bang unless some unknown physics intervenes at the bounce and allows the transfer of the spectrum to the constant mode.
The curvaton mechanism may convert the initial isocurvature fluctuations of some subdominant fields to adiabatic. In the case of the Pre-Big Bang scenario, the Kalb–Ramond axion develops a scale-invariant spectrum for a particularly symmetric Pre-Big Bang evolution and is able to pass all observational and theoretical constraints, provided the bounce scale is at M or below.
In addition to these, we should also mention cosmological models in which perturbations are not the outcome of the amplification of vacuum fluctuations, but have a thermal origin such as the string gas cosmology [
20,
63]. The detailed realization of this scenario, however, is made difficult by our ignorance of string theory.
All these proposals indicate more convoluted routes to explain what standard inflation predicts without apparent difficulties. However, we already know that inflation is too simple to be the end of the story and that something is hidden behind it. Even more intriguing, the existence of some yet unsolved tensions warns us that the Universe is indeed more complicated than the standard concordance model we hoped for just a few years ago. So, it is not impossible that we need to complicate standard inflation at a similar level as bouncing cosmologies in order to account for all the details.
We may wonder whether any distinctive features exist that may rule out one scenario or the other. Indeed, scalar perturbations are sensitive to contributions from all possible fields contained in the Universe. It is not impossible to obtain an acceptable spectrum in both scenarios. On the other hand, tensor perturbations are dramatically different: a nearly scale-invariant spectrum for inflation contrasts with the steep blue spectrum obtained in all bouncing cosmologies. In this case, it is difficult to imagine any mechanisms that may bring the tensor spectrum back to scale-invariant. Therefore, the discovery of tensor modes in CMB may definitely rule out bouncing scenarios or relegate them to a theoretical UV completion to standard inflation with no observational consequences. Conversely, a missing detection of tensor modes by more and more precise probes may put inflation in a difficult position. Finally, the continual increase in sensitivity of gravitational wave interferometers may lead to the discovery of the blue end of the spectrum of stochastic gravitational waves and provide a surprising validation of bouncing cosmologies.