Primordial Gravitational Waves Spectrum in the Coupled-Scalar-Tachyon Bounce Universe

We derive in detail the equations of motion for the tensorial modes of primordial metric perturbations in the Coupled-Scalar-Tachyon Bounce Universe. We solve for the gravitational wave equations in the pre-bounce contraction and the post-bounce expansion epochs. To match the solutions of the tensor perturbations, we idealise the bounce process yet retaining the essential physical properties of the bounce universe. We put forward two matching conditions: one ensures the continuity of the gravitational wave functions and the other respects the symmetric nature of the bounce dynamics. The matching conditions connect the two independent modes of gravitational waves solutions before and after the bounce. We further analyse the scale dependence and time dependence of the gravitational waves spectra in the bounce universe and compare them with the primordial spectrum in the single field inflation scenario. We discuss the implications to early universe physics and present model independent observational signatures extracted from the bounce universe.


Introduction
Encouraged by the direct detections of the gravitational waves in the spiralling blackhole and neutron star systems [1][2][3][4][5][6][7] exactly one century after the prediction of their existence by Einstein, a number of observational effects are being set up to detect primordial gravitational waves originated from the "birth" of our currently observed universe, hence called the "primordial gravitational waves" or "relic gravitational waves" [8][9][10][11][12][13][14][15][16]. This shall open up a new window in the study of early universe cosmology and high energy physics, which, in turn, have profound impacts on our understanding of the physics of the earliest epoch of our universe: Gravitational waves decouple from the hot plasma at thermal equilibrium at T decouple ∼ T P lanck , much earlier than all other particles and propagate freely afterwards. Gravitational waves thus capture the snapshots of the stochastic background of the universe which encode invaluable quantum-gravity information of this earliest epoch of our universe. They encode information of physical processes at extremely high energy otherwise unattainable by particle accelerators on Earth and will perhaps shed light on the origin of matter. In short gravitational waves physics from the earliest epoch of the universe thus nicely compliment and greatly enhance the existent discriminating power of cosmology data brought about by the current era of precision cosmology.
Bounce universe -in which the early universe is proposed to have a phase of matter-dominated contraction prior to the usual "Big Bang" and the ensuing (but not necessarily exponential) expansion -has emerged over the past decade or so to be a viable alternative to inflation. The discovery that the primordial spectrum of density perturbations generated in the contraction phase being (albeit naively) scale invariant in a seminal paper by D.Wands [17] (See also [18,19].) has since then encouraged a continual stream of research effort. Reviews have been written over the years to timely document this line of research activities [20][21][22][23][24][25][26].
In this project we are going to obtain a complete solution to the gravitational waves equations in a simplified bounce universe model, taking careful account of two possible matching conditions [27,28]. We also compute the gravitational waves spectrum generated from this bounce universe and compare it with the gravitational waves spectrum in a single field inflation model, recently obtained in [29]. The bounce universe model, the gravitational waves spectrum of which we will analyse in detail, was derived from low energy effective action of string theory. It built on the idea of non-BPS D-branes and anti-D-branes annihilation [30][31][32]. The tachyon effective potential had been used, independently, by Gibbons [33] and Sen [34] to model inflation.
By introducing a low energy coupling between the tachyon and Higgs fields living on the D-branes we were able to make a model of a bounce universe, called the CST bounce universe. Not only can the CSTB model generate enough inflation to solve the "flatness," "homogeneity," and "horizon" problems of the Big Bang model [35,36], it can also address the "cosmic singularity" problem because the universe could be dynamically stabilised at a non-zero minimal scale. In [37] we obtained the primordial spectrum of the scalar perturbations and established its scale invariance and its stability under time evolution. In the latter series of work [25,[38][39][40][41] we undertook a model independent study of dark matter creation and evolution in a bounce universe. We were able to extract a characteristic curve relating dark matter coupling and dark matter mass should the production of the dark matter take place via an out-of-thermalequilibrium route in a bounce universe. This paper is organised as follows. In Section 2 we discuss the dynamics of a bounce universe model built out of the tachyon and Higgs fields, and their low energy effective coupling as much as it is needed to study the dynamics of the tensorial modes of metric perturbations in the bounce universe. Even though the equation of curvature perturbations is simple and elegant in the CSTB universe, the equations of tensor perturbations become more complex due to a frictional term. In Section 2 we derive the equation of tensor modes by varying the Lagrangian of coupled scale and tachyon fields in a perturbed background metric. In Section 3 we analyse the behaviour of the tensor modes as they cross the effective horizon. After a qualitative discussion, we solve the tensor perturbations equations in the contracting and the expanding phases of CSTB universe and put forward two possible matching conditions. The different power spectra arising from the continuous matching and the symmetric matching condition are then presented. In Section 4 various observational signatures are computed for the tensor modes from the bounce universe and their implications discussed. A side-byside comparison with the single-field inflationary scenario is presented in Section 4.2; and afterwards we conclude with a summary of results and an outlook. Pertinent resulted of the GW spectrum from single-field inflation recently obtained by Alba and Maldacena [29] is reproduced in Appendix A for easy reference and comparison.
2 The dynamics of tensor perturbations in the CSTB universe A string cosmology model, called the Coupled Scale-Tachyon Bounce (CSTB) universe, starts its life from the phase of tachyon-matter-dominated contraction, onward to a bounce at H = 0 and then into a phase of post-bounce expansion from which our observed universe is produced. A scale invariant and stable spectrum of primordial matter density perturbations is produced during the contraction phase [37]. The spectrum of the scalar perturbations is then evolved through the bounce; and its scale invariance and stability are shown to be unaffected by the dynamics of the potentially strong gravity during the bounce process via an AdS/CFT analysis [42]. The rest of this paper is devoted to the computation of the tensor modes of metric perturbations. The spectrum will be evolved to the time of proton-electron capture and obtain the B-mode correlations in the CMB polarisations. The spectral index and its running as well as the tensor-to-scalar ratio of the metric perturbations will also be presented.
The primordial density perturbations generated during the contraction phase were first found to be scale invariance by Wands [17] (See also [18,19].). The horizon exit and reentry of these primordial density perturbations are indeed similar to the analogous process in the inflationary scenario. However a closer look of the detailed dynamics reveals that each wavelength of the perturbations in a bounce universe exits the Hubble radius at different times. This makes each k-mode of the perturbations pick up an implicit time dependence which adds to the difficulty of model building [19]. Care must be exercised when computing the spectrum of primordial perturbations. After taking careful account of this implicit time dependence the spectrum of CST bounce model is proven to be truly time independent and scale independent [37].

Dynamics of the CST bounce universe
In this section the dynamics of the background universe will be recalled as it is needed when we discuss the dynamics of the tensor modes inside the bounce universe. Just as the name coupled scalar-tachyon suggests, the CSTB model utilises a tachyon field T and its interaction with a scalar (Higgs) field φ 1 . The two-field potential creates a dynamical false vacuum when the tachyon T is being pulled up its potential hill and locked at the peak by the fast oscillations of the scaler field φ around the potential minimum of the latter. The minimal coupling of these two fields introduced in [35] makes a bounce universe possible. CST bounce model utilises and extends the tachyon inflation models first introduced by Gibbons [33] and by Sen [34].
Analytically, the coupled scalar-tachyon bounce model can be described by an effective field theory, the Langrangian density of which is comprised of (2.1) In (Eq. 2.1), L (T ) stands for the open string tachyon Lagrangian, L (φ) is the Lagrangian of the scalar and −λφ 2 T 2 a low-energy effective coupling between the tachyon and the scalar, with λ denoting a dimensionless coupling constant. The tachyon Lagrangian is written as (with the metric convention, diag{−, +, +, +}): with V (T ) corresponding to the weakly attractive potential between a stack of D-branes and anti-D-branes [30][31][32][33][34]: . (2. 3) The tachyon Lagrangian describes the onset of an annihilation process of a pair of D4-branes and anti-D4-branes in a closed universe 2 .
In the string language, the scalar field φ can be simply viewed as the distance between the two stacks of D-branes and anti-D-branes.
The Friedmann equation is thus, In the CST bounce universe model, the energy density of the isotropic and homogeneous background universe is given by, where all the fields are functions of time only. 8πG N ≡ m p −2 will be used interchangeably in the rest of the paper. The cosmic evolution of the universe can be obtained, once the initial boundary conditions of the universe are prescribed. As we shall see the prescription of boundary conditions in the bounce universe is an integral part of the ensuing discussion on gravitational waves spectrum.
Cosmic evolution before the bounce -a matter dominated contraction phase: In a bounce universe model the universe is postulated to undergo a period of matter-dominated contraction in which the universe comes into thermal contact to restore causality, and hence addresses the Horizon Problem. Casted in a closed FRLW background, the CST bounce universe undergoes a contraction phase dominated by the tachyon matter during a reversal of tachyon condensation [44]. As the universe contracts, each mode inside the universe gains energy. In particular the tachyon is being pulled up its potential hill as a result of its coupling with the Higgs field; and eventually, being dynamically locked at its false vacuum by the fast oscillations of the Higgs. The stack of D-branes and anti-D-branes approach each other and this instability at the onset of co-annihilation is reflected in the tachyon potential as the tachyon is approaching its potential maximum. The stack of D-branes and anti-Dbranes approach each other due to the weakly attractive potential between the stacks of D-branes and anti-D-branes [30,33,34,45]. The energy of the universe during this process is progressively dominated by this vacuum energy. During the reversal of tachyon condensation the Equation of State changes from w = 0 (tachyon away from its peak) to w = −1 (as the tachyon climbs up its potential peak), at the same time the underlying universe goes from a contraction phase to a deflation.
The Bounce Point: At the bounce point of the CSTB universe, it is the interplay of the curvature term and this vacuum energy that brings about a bounce during which the underlying universe evolves from H < 0 (contraction) to H = 0 (bounce), and then H > 0 (expansion) . Furthermore the universe is dynamically stabilised at a minimal radius, 1 Figure 1. A two-field potential can create a dynamical false vacuum during the bounce phase. The blue solid line is the potential hill of tachyon field T while the red dashed line corresponds to the potential valley governing the Higgs field φ while at the peak of the tachyon potential. When the Higgs field is performing fast oscillation in its potential valley, the tachyon field is dynamically locked at its zero point. The locked tachyon field behaves like "dark energy" sourcing an exponential expansion (the locked inflation), which resembles an inflationary universe. As the universe expands, the amplitude of the scaler field is redshifted and D3-anti-D3-brane pairs annihilate to trigger the tachyon condensation. The tachyon field finally rolls down the blue potential hill to "reheat" the universe.
when H = 0. g s is the perturbative string coupling constant and m s is the string mass scale. As the tachyon reaches its potential peak it is being dynamically locked at its potential peak by the fast oscillations of the Higgs field around its minimum at φ = 0. The Friedmann equation is approximately the same throughout with the Hubble parameter goes from H < 0 to H = 0 and then to H > 0 corresponding to the background universe going through a deflation, a bounce and then a period of inflation. Fig. 1 depicts the two-field potential and how the dynamical false vacuum takes shape. A locked inflation ensues in the false vacuum.
The cosmic evolution after the bounce: At the peak of the tachyon potential the universe undergoes a "locked inflation" driven by the tachyon's vacuum energy. As the locked inflation proceeds, the amplitude of φ oscillations are red-shifted. As the Higgs oscillation energy ceases to dominate the energy density, the tachyon is released from its potential peak and it undergoes a tachyon condensation while the background universe is expanding in power law. The tachyon behaves like a cold matter [44,[46][47][48]. In other words, the equation of state (EoS) changes from w = −1 to w = 0 during the tachyon condensation 3 . After the tachyon condenses, the universe goes onto a decelerating expansion stage driven by the tachyon matter (w = 0), which evolves eventually into our present universe.

Tensorial Perturbations
Different from the inflationary scenario, the Coupled-Scalar-Tachyon bounce (CSTB) universe model generates perturbations during the pre-bounce contracting phase. In the cosmological perturbations theory, the metric perturbations are sourced by the perturbations of the background matter fields. To start the analysis, the classical Lagrangian of CSTB model is expanded up to second order in perturbations, δT and δφ, where R and g, being the Ricci scalar and the determinant of the metric respectively, should also be expanded to second order in metric perturbations: µν , T (0) , φ (0) and a are governed by the classical equations of motion. To study the perturbations we take background spacetime to be [49,50], There are two independent spin-2 tensorial modes of perturbations in thek =ẑ direction: With the perturbed action (Eq. 2.7) one obtains the equations of motion for the two independent modes of metric tensors perturbations, h + and h × : (2.10) and The equations of motion (Eq. 2.10) and (Eq. 2.11) show that the spatial derivatives of the perturbations always couple to the spatial derivatives of the background fields to first order. Under the assumption of the Cosmological Principle that the Universe is homogeneous and isotropic at large scales, the spatial derivative of the background fields vanishes. The equations of motion (Eq. 2.10) and (Eq. 2.11) thus reduce tö Applying Fourier transformation to (Eq. 2.12) and going to conformal time, η, we obtain, an equation obeyed by each Fourier mode, h k . denotes a derivative with respect to the conformal time and h k = h k (η, k).
Discussion: One can indeed check that the perturbation terms of the tachyon and scalar fields sourcing the tensor modes originate from g µν (0) ∂ µ T ∂ ν T and g µν (0) ∂ µ φ∂ ν φ in the perturbed Lagrangian (Eq. 2.7). Only those terms that are multiplied by the tensor perturbations h + and h × in the perturbed Lagrangian remain, upon variations with respect to h + and h × . As a consequence, the spatial derivatives of the perturbations coupled to the spatial derivatives of the background fields can eventually remain in the equations of motion for the tensor modes. The time derivatives of h and those of δT and δφ never have a chance to couple with the tensor perturbations of h to the quadric order. This is simply because the tensor perturbations are the metric perturbations in the purely spatial part. This confirms an important observation: tensorial perturbations of a homogeneous and isotropic universe cannot be sourced by the background fields and their fluctuations [27,49,50].

The primordial tensor perturbations in a simplified bounce
The evolution of the background universe is governed by the classical equations of motion. In the context of cosmic perturbations theory the known array of "matter" content is seeded by the small/quantum fluctuations of fields above this classical background. Care should be exercised to ensure that these primordial "matter" density perturbations remain small perturbations at all time; otherwise these matter constituents would have contributed significantly in the cosmic budget and influence the evolution of the underlying cosmic background. In [37,51] the spectrum of primordial density perturbations of the scalar modes are obtained. Care has been taken to compute the implicit time dependence in the k-modes as they exit their effective horizon at different times. The final spectrum of primordial density perturbations is shown to be scale invariant as well as stable against time evolution [37,42].
Likewise the dynamics of the tensorial perturbations in the metric can be studied in an analogous manner. For the current analysis of the tensorial modes in metric perturbations, we write down the analytic expressions of the scale factor using the EoS at each cosmic epoch: with the subscript 'ec' corresponding to the end of the reverse tachyon condensation [47]. This corresponds to the end of the matter-dominated contraction epoch in the CSTB universe, after which the tachyon will get up to the peak of its potential hill and trigger a period of locked inflation. We shall henceforth denote the Hubble parameter at this point by with the corresponding scale factor denoted by a ec and the physical time denoted by t ec .

The horizon crossing of tensor modes
From the dynamics of the scale factor in each epoch one can introduce the effective horizon as seen by each Fourier mode by taking the inverse of the Hubble parameter, Note here we have made no prior assumption of the Hubble parameter being a constant, as it is usually assumed in the case of inflation scenario. The effective horizon is simply the Hubble radius in conformal time, as shown in Fig. 2.
The horizon crossing conditions require that the wavelength of a given k-mode be equal to its effective Hubble horizon: The scale factor evolves as a power law in conformal time during the tachyon matter domination, leading to k|η|= |ν| .
Effective Horizon Figure 2. The evolution of the effective horizon in conformal time, in which the k-modes are the "co-moving" k-modes and do not change with the conformal time. It is thus obviously that one k-mode crosses the effective horizon four times (the black dots on the diagram) in a bounce universe. We will henceforth simplify the bounce phase (the light blue region) by a point; only the exit and reentry during the tachyon dominated contracting and expansion phases (the dark blue region) are considered when solving perturbations equations (Eq. 3.8).
Analysis of the scalar perturbations in [37,51] lends the theoretical support to the said idealisation.
One can see from Fig. 2 that long wavelength modes exit the horizon earlier and reenter later than short wavelength ones. As the conformal time approaching zero in (Eq. 3.6), the inequality k|η|< |ν| implies that the wavelength of the k-mode is larger than the effective horizon. On the other hand, the large time limit η → ∞ corresponds to the epochs when the k-modes are well inside the horizon. In summary, kη → 0, outside the horizon; kη → +∞, in the far future.

Idealisation of the bounce process
Since the coupled tachyon condensation and the reversal of tachyon condensation [46] persist a short time relative to the contraction and expansion phases of CSTB universe, we could have set them to two points on the time axis when studying the cosmic evolution [35,36]. This is especially the case when studying dynamics of the primordial perturbations of matter [37]. Furthermore we shall ignore the evolution at the bounce point by treating the bounce phase as a point. This is justified by the fact that we are interested in the large scale structures resulting from the gravitational waves: the long-wavelength modes are not sensitive to the small scale fluctuations in the cosmic background as they have long decoupled from the dynamics (out of the horizon at a much earlier time).
The evolution of the CSTB universe is then turned into a two-phase evolution: the pre-bounce contraction and the post-bounce expansion, with the deflation, the locked inflation, and the inflation be collectively represented by "the bounce point." As a result, the scale factor of a bounce universe is given by where a ec is the scale factor at which the end of the tachyon-matter-dominated contration phase. Note that a min is the smallest scale the CSTB universe can possibly attain at the "bounce point" given by which is lower than the Planck scale, where g s is the perturbative string coupling, l s is the string length and l pl is the Planck length.

Solving the tensor mode equations
With the equations of motion for the Fourier modes of the tensor perturbations (Eq. 2.13) is given in Section 2 above, one applies the usual change of variables, u = ah, and arrives at the familiar equation, In the matter-dominated phases of the CSTB model, the scale factor evolves by power law with ν = 2. The solution of (Eq. 3.11) takes a general form (3.12) Whereas the coefficients C k and D k are constants in η they can be functions in k, to be determined by initial conditions and matching conditions. Thus (3.13) are the most general solutions to (Eq. 3.10), with J j being the Bessel functions in the j-th order. The two branches of solutions in (Eq. 3.13) correspond to the contraction (η < 0) and the expansion (η > 0) phases. They are to be connected by the matching conditions at the bounce point when η = 0. Upon obtaining the complete wave function of gravitational waves in the bounce universe one can then study the properties of the primordial gravitational waves spectrum and their observational signatures.
Bunch-Davies vacuum as an initial condition: In the far past we expect the universe would have relaxed into its ground state and a natural choice is the Bunch-Davies vacuum as an initial condition for the primordial gravitational waves at a time far prior to the bounce when the modes are well inside the horizon. Matching the first branch of the solution(Eq. 3.14) to the Bunch-Davies vacuum(Eq. 3.18), we find the coefficients during the contracting phase: Given the functional forms of tensorial perturbations, once we find the right combinations of coefficients (C, D) to capture the physics of interest at the bounce point, the entire evolution of gravitational waves in the CST bounce universe is then fully determined. In particular we need to find the rules to connect coefficients (C exp , D exp ) to (C ctr , D ctr ) as the expansion phase meets with the contraction phase at η = 0. In the following we suggest two reasonable matching conditions to capture the essential physics of our interest.

Matching conditions
To state the problems we are facing, we need to match the two solutions (Eq.3.13) to the gravitational waves equation at the bounce point, where the solution (Eq.3.15) in the contracting phase meets with the solution (Eq.3.16) in the expansion phase, when the conformal time approaches zero. This is visualised in Fig 3 and Fig. 4 numerically with H ec = 20, a ec = 1 and k = 1.   Fig. 3(a) and Fig. 3(c) on the left; whereas the gravitational waves functions after the bounce are depicted in Fig. 3(b) and Fig. 3(d) on the right. Our task is to match the two wave functions on the right to a combination of wave functions on the left by relating (C ctr , D ctr ) to (C exp , D exp ) at η = 0.

I. The continuous matching condition
The first matching condition we propose is the conventional one: to ensure the continuity of the wave functions across the bounce by demanding the wave functions and their derivatives be continuous throughout the bounce phase following [27]. (3.20) The D ctr ( Hecaec 2k ) 3 term in the matching condition (Eq. 3.20) is much larger than C ctr and thus dominating C exp .

Remarks on matching conditions
The seemingly complicated physical process of "bounce" does allow a simple handling when gravitational waves spectrum calculations are concerned, as done in Section 3. The spectrum of primordial gravitational waves we are interested in lies in the low frequency (low energy) range. These modes of gravitational waves decouple much earlier long before the scale of matter genesis sets in. Furthermore the minimal radius of the Coupled-Scalar-Tachyon bounce universe is at/below the string scale (3.9). Even though the physics may be complicated at the bounce point, it does not affect the low-energy modes of gravitational waves which we are interested in. Thus we can simplify the bounce phase as we did in Section 3. For our purpose we note that the overall continuity of GW wave-functions at the idealised bounce point may break down. Fortunately there is another way to fix the coefficients across the bounce point.

II. A symmetric matching condition
If we examine the scale factor before the bounce as well as after the bounce, we find that the derivative of the scale factor gets a minus sign through a bounce point with the same expression. We thus expect the derivative of tensor modes changes by a minus sign after the bounce point. This allows us to fix the coefficients after the bounce by (3.21) Combining the initial condition provided by the Bunch-Davis vacuum (Eq. 3.19) and the symmetric matching condition (Eq. 3.21), we obtain the expressions of (C exp , D exp ) In Fig. 4 the symmetrically matched wavefunctions are shown in green.
The power spectra resulting from the two matching conditions The full analytical expressions for the tensorial metric perturbations throughout the bounce can thus be obtained by combining the matched coefficients after the bounce given by (Eq. 3.20), or (Eq. 3.22) as the case may be, with the coefficients before the  In subfigure 4(a), the blue curve corresponds to the GW wavefunction resulting from the continuous matching while the green curve corresponds to the symmetrical matching condition. Subfigures 4(b) amplifies the details of the wavefunctions at the vicinity of the bounce point. One sees that a smooth curve is obtained with the the continuous condition; while a symmetric wavefunction is obtained with the symmetric matching. After the bounce point, both matching conditions produce oscillatory modes. Subfigures 4(c) and 4(d) show the scales of these oscillations. One also notes that the continuous matching condition produces a significantly larger amplitude after the bounce than the symmetrical matching condition.
bounce as provided by the Bunch-Davies vacuum (Eq. 3.19). The full wave-funtions are summarised in Table 1 for the continuous gluing condition and the symmetric gluing condition.

The time dependence in the power spectra
The η-dependence of the gravitational wave functions resulting from both gluing conditions can be found in Table 1. We see that the wave function is indeed symmetric in the case of the symmetric matching. While being outside of the horizon, the physical

Continuous Matching
Symmetrical Matching Analytical Solution η ν Analytical Solution η ν Before horizon exit 1 √ 2k From bounce to reentry  Table 2. Scale dependence followed from two different matching methods in the simplified bounce, and from the inflationary universe modes are not frozen in either cases. This is a remarkable difference from inflation universe, as already noted in [17]. The primordial perturbations with the symmetrical matching have the same form as the BD vacuum while the continuous matching condition produces much larger oscillations. Both matching conditions yield the same η-dependence after the horizon reentry as the primordial perturbations decays in the same way as the universe expands.

The k-dependence in the power spectra
We compared the k-dependence followed from two different matching conditions in the simplified CST bounce universe with that of the single-field inflationary scenario [29]. The results are presented in Table 2.

From power spectra to observations
Two complete power spectra -one by imposing the continuity of wave-functions and their derivatives at the bounce point, and the other one by insisting on the symmetry of the wave-functions about the bounce point -are obtained in Section 3 above. A few salient points of the power spectra are discussed here. The continuous matching condition generates a power spectrum identical to the one from single-field inflation, as recently obtained by Alba and Maldacena [29], because h is proportional to k − 7 2 (as shown in Table 1). A scale-invariant primordial power spectrum is generated upon horizon exit; a red-tilt present-time spectrum results after horizon reentry. One can compare (Eq. 4.1) in the case of bounce with (Eq. A. 27) in the case of inflation. On the other hand, the symmetric matching condition produces another scale invariant primordial power spectrum, thanks also to the horizon exit mechanism. The present time power spectrum after horizon re-entry is, however, blue tilt, To compare the analytical results with cosmological observations, we calculate the ratio of the tensor perturbations to the curvature perturbations and further compute the B-mode spectrum of Cosmic Microwave Background (CMB).

The Tensor-to-Scalar Ratio
To compare the amplitude of tensor modes to the curvature perturbations one computes the "tensor-to-scalar ratio" as defined by, with the tensor modes power spectrum P h given by (Eq. 4.1) and the spectrum of the curvature perturbation defined as The power spectrum of curvature perturbations in the CSTB universe has been obtained in [37], where the spectrum of δT (P δT ≡ k 3 |δT k | 2 2π 2 ) has been shown to be scale independent and stable against time evolution. As Hubble parameter has time dependence different k-mode exits the horizon at different times. Only when this implicit time dependence is properly taken in account that one can establish the overall time dependence and scale dependence in the final power spectrum of the scalar perturbations in the metric [37,42,51]. Let us recall that in the CST bounce universe perturbations of the matter field δT obey the equation at super-Hubble scales, whose solution is In the discussion on the tensor modes above, we have set the BD vacuum as the initial condition in the far past. The same BD vacuum is set as the initial condition for δT , giving C cur = 1 m pl . Taking kη → +∞ in (Eq. 4.7) we obtain curvature perturbations after the horizon reentry, This leads to the tensor-to-scalar ratio, where κ is a dimensionful quantity defined byṪ c = κH [35,37]. Reinstating the suppressed mass scales we have, with m s being the string mass scale.
The continuously matched power spectrum: The power spectrum resulting from the continuous matching condition can be read off in Table. 1, The symmetrically matched power spectrum: In the symmetrically matched gravitational waves, the power spectrum of tensor modes after horizon re-entry is, given in Table. 1, And the ratio of "tensor to scalar" in the symmetrical case becomes, (4.14) Let us note here the quantity Hecaec k of (Eq. 4.12) is quite big, so is the ratio of "tensor to scalar" (and the amplitude in Fig. 4(a)) in the continuous case, as compared to analogous quantities in the symmetric case. As κ is obtained from numerical solutions of the CSTB model and it is determined by the linear relation between the vacuum expectation value of the tachyon field, < T c >, and Hubble parameter during the phase of tachyon matter dominated contraction, < T c >= κH. The numerical value of κ is pretty stable, 1 < κ < 10, over a wide range of parameter values [37]. It is then unlikely to fiddle with the parametric values of the CSTB model to suppress the high tensor to scalar ratio. Observations of tensor modes, or the lack of it, favours the small tensor to scalar ratio and hence the symmetric gluing condition in a bounce universe. Once the tensor to scalar ratio is measured we can use it to put an upper bound on the ratio of string mass to the Planck mass. This small tensor to scalar ratio is still the Occam's razor for alternatives to the inflationary paradigm.

Tensor modes of CSTB universe versus inflation
In this section we attempt to make a side by side comparison on the similarities and the differences of the gravitational waves spectra obtained from the CST bounce universe model and the single-field inflation model [29]. The complete GW wavefunctions from the bounce model and from the single-field inflation model are plotted together in Fig. 5. The tensor modes of both matching conditions in the bounce case are time-varying while the primordial tensor perturbations of the inflationary universe are time-invariant [17]. The k-dependence in these three cases are summarised in Table 2. From Fig.5 one sees that the continuously matched bounce model shares a similar curve with the inflationary model after the horizon reentry, which explains the same k-dependence and the similar functional forms of these two GW spectra. Detailed calculations relevant to the GW spectrum in the single-field inflation model [29] is reproduced in the Appendix A for easy reference. Bounce-continuous Bounce-symmetrical Inflation Figure 5. The gravitational wave functions in the three cases discussed in the paper.

The BB power spectrum of CMB
The B-modes of CMB, partly generated by the primordial gravitational waves, encodes the history of the universe earlier than the last scattering surface [52][53][54]. In this paper we use the Code for Anisotropies in the Microwave Background (CAMB) [55] developed from CMBFAST [56] to compute the BB power spectrum generated by the primordial GW in the simplified CST bounce universe 4 .
The BB power spectrum is presented in Fig.6, which is plotted by setting the tensor mode index to zero, and the equation of state before horizon reentry to zero corresponding to the tachyon matter dominated contraction era. These are done in line with the inflation model for a fair comparison. We have also computed and plotted the B-modes in the single-field inflationary universe and the BB power spectrum of the lensed-ΛCDM [56]. We have set the tensor-to-scalar ratio to be 0.2 when computing the BB power spectrum generated by primordial GW. From Fig. 6 one finds that the BB power spectrum in the range (10 < l < 50) from the simplified CST bounce universe tracks closely in shape with the spectrum produced by the inflation model and is consistently lower in amplitude.
Considering the data of BB power spectrum of CMB and the lensed-ΛCDM curve,  Figure 6. In both subfigures, the red curves, and blue curves, correspond to the BB power spectrum generated by the primordial gravitational waves of the inflationary universe, and the simplified CSTB universe, respectively. The black dashed curve is the BB spectrum of the lensed-ΛCDM. In Fig. 6(a), the orange dots, with error bars, are BICEP2 data released in 2014 [10] while the green dots are data released in 2015 [57]. The dots inside Fig. 6(b) are the data after subtracting the lense-ΛCDM BB power spectrum part. This is done so to compare the "pure" gravitational waves contributions to the BB spectrum from both models.
the primordial power spectrum may be responsible for the B-modes of CMB at large angle when the ratio is around 0.2. However, with further detection and joint data analysis done by BICEP2 with Planck, it was proposed that part of the B-mode signals was attributed to the Galactic dust and hence the ratio was further compressed [16,[58][59][60][61]. A "lensed-ΛCDM + tensor + dust" framework has since been worked out, and the tensor-to-scalar ratio has thus been suppressed to be lower than 0.07. We hope future detections of the primordial gravitational waves, the like of LISA, Ali, and TianQin, can measure with high precision the primordial tensor perturbations. Only then one can hope to discern the fine details in the gravitational waves spectra as predicted by the large repertoire of early universe models.

Conclusion
In this paper we have derived the equations of motion for the tensorial modes of the metric perturbations in a bounce universe whose dynamics is governed by the Coupled Scalar and Tachyon fields (CST) arising in two stacks of colliding D-branes and anti-D-branes [35,36]. We obtain two sets of independent solutions in each of the following epochs: when the modes are inside the horizon well before the bounce, and well outside the horizon during the matter dominated contraction phase characteristic of the bounce Items Contiuous-matching Symmetrical-matching Inflation Analytical form of tensor modes Table 1  Table 1  Table 4 Primordial  Table 3. Tensor modes comparison of simplified CSTB universe and inflationary cosmos universe, and then finally when the gravitational waves re-enter the horizon while the universe is in matter domination again. As gravitational waves are quite immune to the potentially complicated physics due to matter generation and subsequent developments close to the bounce point due to its extremely weak coupling to matter fields, one can henceforth simplify the bounce process. To obtain the full wave functions of the gravitational waves throughout the whole cosmic evolution and yet capture the interesting physics pertinent to low-l modes of gravitational waves, we propose two gluing conditions for the gravitational waves inside a bounce universe: One insisting on the continuity of the wave functions and their derivatives at the idealised bounce point; the other one respecting the symmetry of the bounce dynamics about the bounce point. Both matching conditions predict scale invariance of the primordial spectrum when the modes are outside of the effective horizon, and oscillation modes of the primordial tensor perturbations after the bounce. However, the continuous matching condition leads to a much larger amplitude of the primordial tensor modes than the symmetric matching condition. The power spectrum of the primordial tensor perturbations of continuous matching is red tilt while that of symmetric matching becomes blue in late time. A detailed comparison is done with the primordial spectrum obtained from the single-field inflation model [29]. We summary all features derived in this paper and their comparison with inflation in Table 5.
Higher precision measurements of BB power spectrum in the Cosmic Microwaves Background can place better constraints on the inflationary paradigm and its alternatives. Direct detections of primordial gravitational waves from LISA, ALI, or TianQin and the future array of low-frequency GW detectors can eventually tell whether our universe went through a bounce or an inflation.
Meanwhile particle physics experiments aiming at direct detections of dark matter and determination of their coupling with Standard Model particles are in full swing. As bounce universe allows an out of thermal equilibrium route of dark matter creation in the contraction phase, dark matter mass and their coupling are constrained by the present-day relic abundance [25,40,41]. One may be in for surprise earlier than expected, as it provides a way to distinguish inflation [62] from a bounce universe.

A The primordial gravitational waves in the Standard Inflationary Cosmology
The primordial power spectrum of the gravitational waves arisen in the inflation model is obtained recently by Alba and Maldacena [29]. We review their calculations in detail here for pedagogical purpose and for quick comparison with our results presented in the main body of our paper.

A.1 The single field slow-rolling inflation and the evolution of the scale factor
Based on Einstein's theory of relativity, one can construct a model to realise inflation with a simple scalar field Lagrangian, L (φ) = −∂ µ φ∂ µ φ − V (φ). The energy density and the pressure in the stress-energy tensor are will be driving an exponential expansion of the underlying universe provided that the EoS (Equation of State) can stay ω = −1 for long enough time. The easiest is to impose the slow-roll conditions [63]: In other words, the kinetic energy of the inflaton is negligible in comparison to its potential energy and its potential hill being ultra-flat.
With an almost constant Hubble parameter during an inflation, one can obtain an analytical expression of the scale factor by solving the Friedmann equations: where the zero time t = 0 is set at the time when inflation ends and a Λ denotes the scale factor of the universe attained at the end of inflation. The minus and the positive branches of the scale factor evolution correspond to the inflationary stage and the postinflation expansion respectively. The index n in the power law expansion after inflation is determined by the EoS of the cosmos, n = 2 3(1+ω) . For the study of gravitational waves in an inflationary universe we recall a few key details of the single field inflation model below.
The scale factor of the inflationary universe (Eq. A.4) and its time derivative remains continuous during the whole cosmological evolution. The physics at the end of the inflation is however complicated due to matter production: "reheating" is a process in which matter is purported to be produced from the decay of the inflaton. For a recent discussion on the interesting physics that can happen during reheating and its implications on reheating temperature, the readers are referred to [64]. In fact one can turn this around and ask how reheating and parametric resonance can have imprints on CMB; and whether CMB can tell us anything about the highest temperature the universe has ever reached in the earliest epoch of evolution. For the present discussion, however, "reheating" and the interactions of matter occurred thereafter are completely ignored, as they play no role in the cosmic evolution around the time when tensor modes exit the Hubble horizon.
When the tensor modes cross the horizon during the inflation and in the postinflation expansion, the analytical form of the scale factor is given by (Eq. A.4). In conformal time η = dt a(t) , the scale factor takes the form: with ν = 2 1+3ω .

A.2 The equation of linearised tensorial perturbations of the background metric in the single-field inflation model
The isotropic homogeneous background universe is described by the FLRW metric with the time variable t being the physical time. The tensor perturbations can be expressed as g µν = g (0) µν + δg µν . The metric tensor perturbations are traceless and divergence-less, There are thus two independent spin-2 tensor modes of metric perturbations, describing the propagation of gravitational waves at the speed of light in thek =ẑ direction.
The equations of tensor perturbations can be obtained by expanding the action, The determinant of the metric, √ −g, and the Ricci scalar, R, can be expanded up to the quadratic order in h ij as follows, and where the dot denotes derivative with respect to the physical time. In the derivation of (Eq. A.10) and (Eq. A.11), the condition (Eq. A.7) has been used. Roman indices i, j, k run over 1, 2, 3.
Expanding the action (Eq. A.9) to the quadratic order and by variational method we obtain the equation of tensor perturbations of the metric, Since we are interested in the full spectrum of gravitational waves, we need the equation of motion for h in Fourier space, where we have changed the time variable to conformal time, η, and prime denotes a derivative with respect to η. Solving (Eq. A.13), we obtain the evolution of gravitational waves spectrum.

A.3 Solving the equation of tensor modes in a single-field inflation model
When the scale factor expands with a power law a ∝ η ν , with the substitution u = ah, (Eq. A.13) can always be transformed to the form, 14) The solution of u(η) is a linear combination of η Applying the scale factor of inflationary universe (Eq. A.5) into (Eq. A.14), the solution to the tensor modes can be expressed as, (A. 15) The two branches of the solutions in (Eq.A.15) correspond to inflation (η < 0) and the post-inflationary expansion (η > 0). They are connected at the end of inflation (η = 0). Typically, the coefficients (C inf , D inf ) of the wave functions during inflation are determined by initial conditions; and the coefficients (C exp , D exp ) of the wave functions during the post-inflation expansion are then related to (C inf , D inf ) at the end of inflation.
The horizon crossing conditions of the tensor modes: The horizon crossing behaviour of the tensor modes during and after the inflation is similar to the analogous horizon crossing process in the bounce universe, the latter has been discussed in Sec. 3.1 and we refer to D. Wands's original paper [17] for more details. Thus applying the Equation of State at horizon crossing (Eq. 3.6), one obtains the condition for a given k-mode of gravitational wave to exit the horizon, k aH = kη = −1 (A. 16) and reenter the horizon, kη = ν. (A.17) The horizon crossing conditions hence divide the gravitational wave function into three stages. The inequality k|η|< |ν| (η → 0) corresponds to physical modes being well outside the horizon, when at large time limit (η → ∞) all k-modes being well inside the horizon. In summary, kη → −∞, in the far past; kη → 0, outside the horizon; kη → +∞, in the far future.
(A. 18) which is the same as in the bounce case (Eq. 3.7).
The short wavelength and long wavelength limits: Applying the limits (Eq. A.18) to the general solutions (Eq. A.15), one obtains the expression for each k-mode in each of the following situation: Bunch-Davies vacuum: The universe is expected to start from a time-independent ground state in the far past. A natural choice of the initial condition of the universe is the Bunch-Davies vacuum, where m 2 pl ≡ 1 8πG N ≡ 1 is the reduced Plank mass. The 1 √ 2k term arises due to the normalisation of the wave function, |u| 2 (2k) = 1, with u k (η) = a(η)h k m pl e −ikη being regarded as a quantized field.
Combining the solution in far past (Eq. A.19) with the initial condition (Eq. A.23), the coefficients during inflation are therefore determined, Once this is done the coefficients of the tensor perturbations during the post-inflationary epoch are also determined using the matching conditions on the super-Hubble scale.

A.4 Matching conditions
It is a natural to assume that the gravitational waves and their derivatives be continuous during the cosmological evolution. The continuity of h(η) and h (η) thus relate the post-inflation coefficients (C exp , D exp ) to (C inf , D inf ) at η → 0 by, (A.25) Regardless of the constant terms, the matching condition acts to transfer k H λ a 0 between coefficients. Since the tensor modes are stretched outside of the horizon near the end of the inflation, the long wavelength limit k H Λ a 0 → 0 makes it possible to compare the contributions of (C inf , D inf ) to (C exp , D exp ). Obviously, the C inf term is ( k H λ a 0 ) 3 smaller than the D inf term in C exp . The matched wave functions are presented numerically, with H Λ = 20, a 0 = 1, k = 1 and ν = 2, in Fig. 7 and Fig. 8.

A.5 The power spectrum
Combining the coefficients during the inflation and post-inflationary expansion and the solutions of every cosmic stage, the analytical form of each epoch is summarized in Table. 4. The primordial power spectrum becomes P h = k 3 |h| 2 2π 2 ∝ k 0 , (A. 26) while the power spectrum after the horizon reentry is 27)