Observable Gravitational Waves from Hyperkination in Palatini Gravity and Beyond

.


Introduction
The most compelling solution to the fine-tuning of the initial conditions of the Big Bang model is the theory of Cosmic Inflation [1,2,3,4,5,6,7].Inflation manages in a single shot to explain away the horizon and flatness problems and also to provide the primordial density perturbations necessary for the formation of the large-scale structure we observe in the Universe [8,9,10,11,12,13], that is the distribution of galaxy clusters and superclusters.The primordial density perturbations reflect themselves onto the Cosmic Microwave Background (CMB) radiation, through the Sachs-Wolfe effect [14].Precision observations [15,16] of the acoustic peaks in the CMB primordial temperature anisotropy have verified in spectacular detail the predictions of Cosmic Inflation, such that the rival paradigm for structure formation (that of cosmic strings) has collapsed [17].Consequently, Cosmic Inflation is considered a necessary addition to the concordance model, ΛCDM, towards a standard model of cosmology.
Another generic prediction of Cosmic Inflation is coming within reach of observability in the near future.Indeed, soon after its proposal, it was realised that Cosmic Inflation gives rise to a stochastic background of primordial gravitational waves [1,18,19,20,21,22].These gravitational waves (GWs) are tensor perturbations of the spacetime metric, generated in much the same way as the scalar curvature perturbations behind the primordial density perturbations, for which there is overwhelming evidence in the CMB, as mentioned above.Because of this, great interest has been developed in recent years for the observability of the inflation-produced GWs either indirectly, through the B-mode polarization of the CMB [18], or directly from interferometers [19].
Gravitational waves were predicted by Einstein's general relativity at the beginning of the twentieth century.Almost exactly a hundred years afterwards, GWs were directly observed by LIGO (Laser Interferometer Gravitational-Wave Observatory) and Virgo in 2015 [23,24].This seminal observation heralded the birth of gravitational wave astronomy, which enables the study of compact objects, such as astrophysical black holes, which are typically shrouded by opaque accretion disks.It also allows, in principle, a glimpse of the very early Universe, well beyond the last-scattering surface, where the CMB was emitted.As such, there is hope to detect the stochastic primordial GW background from Cosmic Inflation.Such observations will allow the study of inflation at scales much different than the ones which correspond to the CMB primordial anisotropy, opening up a new window in the understanding of fundamental physics at extremely high energies (comparable to the energy of grand unification), which is behind the process of Cosmic Inflation and remains a mystery to this day.
This has, in part, motivated a number of future GW detection missions.In the near future, Advanced LIGO (plus Virgo and KAGRA) [25,26,27,28,29] (LVK) and the space interferometer LISA (Laser Interferometer Space Antenna) [30,31,32] are coming online; the launch date of LISA is in 2037.Another space interferometer DECIGO (DECi-hertz Interferometer Gravitational wave Observatory) [33,34,35] is also planned to be launched in the 2030s.More are to follow, such as BBO (Big Bang Observer) [36], a proposed successor to LISA.It seems an ideal time to investigate GW production by inflation and its potential observational signatures.
However, there is a challenge in the study of the inflation-produced primordial GW background.The background signal is too weak for any currently operational GW detector to observe, and it may be decades before an observation can be made.Indeed, were the early Universe dominated by radiation, as assumed by the concordance model, the primordial GW spectrum would be flat, i.e. like white noise, where the GW density parameter per logarithmic frequency interval Ω GW (f ) is constant over the range of frequencies f corresponding to the GW modes that re-enter the horizon during the radiation dominated period (they have been pushed out of the horizon during inflation).The constant value of the flat spectrum is very low, and the hope of detecting in the near future such inflation-generated primordial GWs is little [20].
Fortunately, this is not the end of our hopes for detecting primordial GWs.While there is observational evidence of the early Universe being radiation dominated, provided by the delicate process of Big Bang Nucleosynthesis (BBN) taking place a mere few seconds after the Big Bang itself, what the state of affairs was before BBN is still unknown.If the Universe's history before BBN was not dominated by radiation, then the primordial GW spectrum does not need to be flat.This opens up the possibility of a boosted GW spectrum, possible to detect even in the near future.
An early realisation of this possibility was provided by modelling quintessential inflation [37] (see Refs. [38,39] for recent reviews).Quintessential inflation aims to explain in a unified way both Cosmic Inflation in the early Universe and Dark Energy at present.Most quintessential inflation models consider non-oscillatory inflation [40,41] driven by a scalar field (the inflaton) with a runaway potential, which can play the role of quintessence at late times and explain the accelerated expansion of the Universe at present [42,43,44,45,46,47,48,49,50,51,52,53].In such models, there is a period after the end of inflation but before reheating (i.e. the onset of the radiation era) when the kinetic energy density of the inflaton field dominates the Universe.This period is called kination [54] (see also [55,56,57]), characterised by a stiff equation of state with a barotropic parameter w = p/ρ = 1.For GW modes that re-enter the horizon during kination, the spectrum is peaked with Ω GW (f ) ∝ f [58,59,60,61,62,63,64,65,66].Unfortunately, this peak corresponds to very high frequencies, which will be unobservable in the near future.Extending the period of kination does extend the peak to lower, possibly observable frequencies, but then the peak becomes too large and the resulting primordial GWs cannot but affect and destabilise the BBN process [58,59,67,42].
After the direct detection of GWs, there has been much interest in considering modifications of the history of the Universe, safely before BBN, to boost the primordial GW signal at observable frequencies.In Ref. [68], it was shown that Ω GW (f ) ∝ f −2( 1−3w 1+3w ) , where w is the barotropic parameter of the Universe (w = 1/3 for radiation domination).In Refs.[69] and [70] models were considered where there is a period of matter domination followed by kination, which would create a mountain-like peak in Ω GW (see also Ref. [68]).Another possibility is to consider a stiff period after inflation that is not kination with w = 1, but has a milder value of w ≈ 1/2 and can be extended down to observable frequencies without destabilising the BBN because the peak is not so steep as in usual kination [71].A realisation of this in hybrid inflation with a non-canonical waterfall field was investigated in Refs.[72,73].
In this paper, we consider a different possibility, motivated by Palatini modified gravity [74,75].The cosmological consequences of Palatini modified gravity with L ∝ R + R2 and a non-minimally coupled scalar field were first considered in [76,77] in the context of inflation and subsequently in [78,79,80,81,82] in the context of quintessential inflation (see also [83,84] for reviews).When switching to the Einstein frame, the scalar field obtains an additional quartic kinetic term 1 .In most cases considered, this term plays a negligible role in the dynamics of the scalar field.However, there are models for which this is not the case.We investigate in detail what happens when the scalar field dominating the Universe is governed by the quartic kinetic term in a period we call hyperkination.We show that the barotropic parameter of the Universe during hyperkination is the same as that of radiation domination, w = 1/3.As a result, in a realistic model of non-oscillatory inflation with a runaway inflaton potential, we consider a post-inflationary period of hyperkination, followed by a period of regular kination, when the kinetic energy of the inflaton is quadratic as usual.Kination is followed by radiation domination after reheating.This evolution results in a truncated peak in the GW spectrum, which can be safely extended down to observable frequencies without destabilising BBN.We calculate analytically the GW spectrum during all phases of hyperkination, kination and radiation and we verify our findings numerically.We explore the parameter space and show that we can obtain a boosted primordial GW signal with unique characteristics, which will be well-detectable by forthcoming observations.If such a signal is indeed detected, it will be a strong hint of non-canonical kinetic terms for the inflaton field from Palatini modified gravity or some other appropriate kinflation or k-essence model.
The paper is organized as follows.In section 2, we discuss the Palatini R 2 models, introduce hyperkination, and embed it into the full expansion history of the Universe.In section 3, we consider the primordial GWs, including their initial conditions as fluctuations of the quantum vacuum.Section 4 details our analytical computation of the GW evolution.We compare our GW spectra to observational bounds in Section 5 and conclude in Section 6.Throughout the paper, we use natural units with c = ℏ = 1 and 8πG = m −2 P , where m P = 2.43 × 10 18 GeV is the reduced Planck mass.The signature of our metric is (−1, +1, +1, +1). 1 In Ref. [85] it was shown that the addition of the Holst and Holst 2 terms in the usual Palatini quadratic action can generate a modification of the higher-order kinetic term.

Quartic kinetic terms from Palatini R 2 inflation
We begin by considering a Jordan frame action in the Palatini formulation of the form where φ is the inflation field and h(φ) is its non-minimal coupling function, which usually assumes the form 2 h(φ) = m 2 P + ξφ 2 .The parameter α is assumed to be positive definite and we leave the potential V (φ) unspecified.The symbol ψ describes other matter components.This action was first considered in [76,77] in the context of inflation and then in [78,79,80,81] in the context of quintessential inflation (see also [83,84] for reviews).
In the Palatini formulation of gravity, the connection Γ and the metric g µν are independent variables.The Ricci tensor R µν (Γ ) only depends on the connection and the Ricci scalar is defined as R ≡ g µν R µν (Γ ).The connection Γ can be determined by varying the action (1) but, due to the non-minimal coupling function h(φ) and the αR 2 term, it will differ from the standard Levi-Civita form.
Following [76,77], we can eliminate the αR 2 term by introducing an auxiliary scalar field χ ≡ 2αR.Then, by performing a Weyl transformation of the form ḡµν = Ω 2 g µν = [χ+h(φ)]g µν we bring the action to the canonical form with a minimally coupled scalar field.The resulting action will depend on two scalar fields: φ and χ.However, in contrast to the usual metric formalism, the auxiliary field χ is non-dynamical in the Palatini formalism.This means that one can vary the action with respect to χ, solve the resulting constraint equation, and then eliminate χ altogether from the action.
After this procedure, the resulting action in the Einstein frame reads [76,77] in order to render the quadratic kinetic term canonical, where the bars indicate quantities in the Einstein frame.Note that the process of transforming from the Jordan to the Einstein frame has generated a quartic kinetic term3 and a modified potential U which will in general display a plateau for growing V , approaching the asymptotic value m4 P /(4α) [76].Also, importantly, in the present work, we concentrate on the early era when the other matter components ψ are a perfect fluid of radiation.In this limit, the coupling between the inflaton and the matter action in the last term of Eq. ( 2) disappears [80].
Neglecting the last term for the moment, we can rewrite the action as with where The action in Eq. ( 5) belongs to the general class of k-inflation [89] (where inflation is kinetically driven) or k-essence [90,91,92] (where the non-canonical kinetic terms can behave as quintessence). 4  Varying the action in Eq. ( 2) we can obtain the equation of motion for ϕ, which reads [76] 1 Then, from the non-zero components of the energymomentum tensor we can obtain the energy density and pressure of the field, which read [94] To complete the equations of motion, the Hubble parameter can be written as Again, the above equations differ from those of a standard canonical scalar field due to the higher-order kinetic terms.In the limit α → 0 they reduce to the minimal case.The bars are dropped in what follows to avoid clutter.Unless otherwise stated we always work in the Einstein frame.
The plateau in U mentioned above is ideal for slowroll inflation, and can easily produce CMB observables compatible with observations for simple forms of the potential V [76,77].However, it restricts the inflationaryand thus post-inflationary-energy density to values lower than m 4 P /(4α).Unfortunately, this severely restricts the parameter space considered in the following sections.One way to overcome this problem is to consider an α that experiences a drastic change at the end of inflation but remains constant afterwards.This is possible if α depends on a degree of freedom that changes its value when inflation ends.A toy model discussing this possibility is presented in Appendix A. Another example of a model describing the full inflationary history may be the one studied in Ref. [79], as long as it is enhanded with the hybrid mechanism discused in Appendix A. Moreover, we point out that the Palatini R 2 models considered here act as an inspiration for the extra quartic kinetic terms in the action, but our analysis is more general, and we do not specify the details of the inflationary part of the model.

Kinetic domination
While the quartic kinetic terms in Eq. ( 2) are negligible during slow-roll inflation [76], they may play an impor-tant role in the post-inflationary Universe.We consider next such a scenario; a period of kinetic domination, where the potential V is negligible and the field rolls forward freely.In this limit, Eqs. ( 8) and (10) become 1 + 3α φ2 and It is instructive to change the time variable to the number of elapsing e-folds N = ln a, with dN = Hdt, and eliminate H.We can assume φ > 0 without loss of generality.The field time derivatives are related as Note that, due to the scaling with the heavily φ-dependent H, the limit φ → 0 corresponds to ϕ ′ → √ 6m P , and φ → ∞ corresponds to ϕ ′ → 0. Eqs.(11) and ( 12) become where w ϕ ≡ p ϕ /ρ ϕ is the barotropic parameter of the field.Note that α dropped out of the equation of motion: changing α rescales the time and energy density but leaves quantities like ϕ, N , and w ϕ untouched.If φ is small-that is, 3 2 α φ2 ≪ m 4 P and ϕ ′ ≈ √ 6m Pthe quartic extra kinetic terms are small, and Eq. ( 14) give where c is an integration constant and we are concerned with the large N limit.We see that ϕ ′ = √ 6m P is an attractor.It corresponds to standard kination [95,96,54,43,97,98] with a quickly diluting energy density and w ϕ ≈ 1.
In the opposite limit of 3 2 α φ2 ≫ m 4 P and ϕ ′ ≈ 0, the quartic kinetic terms dominate, and Eq. ( 14) gives We name this phase hyperkination.The extra kinetic terms modify the dynamics so that the energy density dilutes only as fast as radiation with w ϕ ≈ 1/3.Hyperkination only lasts for a limited time.As spatial expansion dilutes the field's kinetically dominated energy density, φ decreases and ϕ ′ grows.The quartic kinetic terms are diluted faster than the quadratic ones, and eventually the latter take over.Consequently, the field transitions into standard kination.We can use Eqs.( 15) and ( 16) to approximate the time evolution of ϕ ′ as it approaches the kination attractor as where ϕ ′ 0 is the initial value of ϕ ′ at N = 0, taken below to be the end of inflation.Tuning ϕ ′ 0 lets us modify the length of hyperkination, which we define as 6 Figure 1 compares Eq. ( 17) to a numerical solution of Eq. ( 14) in an example case.Due to the exponential growth of ϕ ′ , the transition from hyperkination to kination is fast.Let us define the beginning of standard kination N kin as the moment when both addends inside the parenthesis in the energy density in Eq. ( 12) become equal.Using Eqs. ( 13) and (17), this condition reads Thus, N hyp ≃ N kin and we conclude that it is a good approximation to assume an instantaneous transition between hyperkination and kination.We end this section with a relation between α, the energy density at the start of hyperkination (end of 6 With the restriction ρ ϕ < m 4 P /(4α) discussed at the end of section 2.1, we would have ϕ ′ 0 > 2m P at N = 0, and Eq. ( 17) restricts hyperkination to last less than 0.20 e-folds, a negligible amount.As mentioned, we omit this restriction in this paper.H=10 13 GeV N kin Fig. 1: N -derivative of the field obtained from the numerical simulation (full blue line) and its initial approximation given in Eq. (17) (dashed orange line) as functions of N .The dashed vertical line, labelled N kin , corresponds to the time at which kination starts in the numerical simulation, defined here as the moment at which both addends inside the parenthesis in the energy density in Eq. ( 12) become equal, while the dashed horizontal line corresponds to ϕ ′ = √ 6m P .In the legend, H denotes the Hubble parameter at the end of inflation H end .
inflation) ρ end , and N hyp .Equation ( 14) together with the definition of N hyp gives where we have assumed a non-negligible duration for hyperkination O(N hyp ) ∼ 1 in the last step.Note that N hyp = 0 corresponds to α = 0, as it should.

Full cosmic evolution
Let us now embed a period of kination into a full history of the Universe.Initially, during cosmic inflation, the field energy density is dominated by potential energy.
Once inflation ends, the potential drops to zero and the field's velocity increases as the potential energy is transformed into kinetic energy.In typical models, the field is trapped into a potential minimum, oscillating there and decaying into a thermal bath of particles, reheating the Universe.In our models of interest, the postinflationary potential is of the runaway type-that is, flat and low-and the field keeps rolling onward under kinetic domination.If the quartic kinetic terms dominate, this phase starts with hyperkination, transitioning into standard kination later, as described above.
To reheat the Universe, we assume a small amount of radiation is produced at the end of inflation e.g. through Ricci reheating [99,100,101].During hyperkination, the radiation energy density is diluted as fast as that of the field, ρ r,ϕ ∝ a −4 , so radiation stays subdominant.However, when standard kination starts, the field energy density dilutes faster, ρ ϕ ∝ a −6 , and the radiation fraction grows until it overtakes the field.The Universe is reheated and radiation domination starts.We assume this to take place at high energies, above the BBN temperature T BBN ≈ 1 MeV; afterwards, the Universe follows the standard ΛCDM expansion history.
The behaviour of the system can be solved from the Friedmann equations combined with the first equations from Eqs. ( 11) and (12).Here Ω end r is the radiation energy density fraction (parameter) at the end of inflation and a end and H end are the scale factor and Hubble parameter at the end of inflation.Figure 2 shows the behaviour of the energy densities solved numerically; for details on the numerical implementation, see Appendix B. It also shows the corresponding evolution of the barotropic parameter w, defined as the ratio between the total pressure and energy density of the Universe, taking values from w = 1/3 (hyperkination) to w = 1 (standard kination) back to w = 1/3 (radiation domination).
In summary, we assume a cosmological evolution where inflation is followed by two phases: hyperkination and kination, in this order.Reheating, which takes places at temperatures larger than T BBN , signals the end of these phases.After reheating, the conventional cosmic evolution with radiation and matter dominated eras follows 7 .
The non-standard expansion history opens the door for new phenomenology.For one, it changes the matching between scales in the early and late Universe.Indeed, when inflation is followed by a stiff cosmological era with barotropic parameter w, the number of inflationary e-folds is increased by [48,57,78,103] It follows that hyperkination, for which w = 1/3, has a vanishing contribution.This is not the case for kination, with w = 1.Thus, in our scenario we have where ρ kin ≪ V end is the energy density at the end of hyperkination and the onset of kination proper.Typically, this increases the remaining number of inflationary efolds after the cosmological scales exit the horizon to at most N * ≃ 65, which implies ∆N ≲ 5, something that must be taken into account when calculating the inflationary observables. 8ll in all, the CMB scales exit the Hubble radius approximately 60-65 e-folds before the end of inflation instead of the standard 50-60, see e.g.[78,107].This affects inflationary model building, although the effects are mitigated with respect to the standard kination scenario.In addition, the spectrum of primordial GWs is altered in ways that are sensitive to the duration of hyperkination.

Tensor perturbations
To study the behaviour of GWs, we write the metric tensor as g µν = a 2 (η µν + h µν ), where η µν is the Minkowski metric so that a 2 η µν ≡ ḡµν is the unperturbed FLRW metric, and h µν is a small perturbation.We expand the action in Eq. ( 2) to second order in h µν , keeping only the tensor modes 9 , which evolve independently of other perturbations in linear perturbation inflation [1] or Higgs inflation [104] (or α-attractors [105]), the scalar spectral index is . With N * = 60 this results in n s = 0.966.If we have N * = 65 instead, then n s = 0.969, which is still within the 1-σ contour of the Planck satellite observations [106]. 9The tensor perturbations obey ∂ µ h µν = 0 and h µ µ = 0.
theory.The result is (compare to, e.g., [108]) where η is the conformal time related to the cosmic time t by dt = a dη, and a prime denotes a derivative with respect to η.Here s indexes the two gravitational wave polarisations, and the polarization amplitudes h s are defined through the Fourier decompositions so that h s ⃗ k describes oscillations of a given polarization in directions perpendicular to the wave vector ⃗ k.
The amplitudes h s behave as massless scalar fields, up to normalization, following the Klein-Gordon equation with wave solutions.Here H ≡ a ′ /a and ∇ 2 ≡ ∂ i ∂ i where i is summed over the spatial indices.The corresponding energy-momentum tensor is so that the GW energy density reads

Quantization
The primordial GWs originate as quantum vacuum fluctuations during inflation.To quantize them, we first go to the canonically normalized variables v s = m P ah s / √ 2, so that (after integration by parts) the action in Eq. ( 25) becomes This is the Minkowski space action for a free field with mass a ′′ /a, quantized the standard way by writing where âs ⃗ k , âs † ⃗ k are the ladder operators following the canonical commutation relations Time evolution is delegated to the mode functions v s k , which follow the Mukhanov-Sasaki equations derived from Eq. (32), Note that, due to the ladder operators, the mode functions v s k differ in normalization from the classical Fourier modes v s ⃗ k .Abusing the notation slightly, we differentiate these by writing k instead of ⃗ k as the mode function index-in an FLRW background, the quantum mode functions only depend on the magnitude of the wave vector and not its direction.Analogously, we define ĥs = √ 2v s /(am P ), h s k = √ 2v s k /(am P ).Deep inside the Hubble radius, k ≫ H, the GWs do not feel the expansion of space, the mass term a ′′ /a is negligible, and Eq. ( 35) has the standard vacuum solution When the mode functions follow Eq. ( 36), the state annihilated by âs k is the Bunch-Davies vacuum [109]; we take the perturbations to start in this vacuum state during inflation.Over their cosmic evolution, the modes stretch and exit the Hubble radius, evolving beyond Eq.(36).After inflation, they re-enter the Hubble radius, this time following the general sub-Hubble form We will solve the coefficients λ ± (k) for a given cosmic history in section 4.2; since the Mukhanov-Sasaki equation conserves the Wronskian of its solutions, we have , set by the initial vacuum in Eq. (36).The coefficient λ − contains the GW excitations, the part beyond the vacuum solution in Eq. (36).
Let us next consider the energy density of the GWs induced by the above process.The late-time GW energy density is dominated by high-k, sub-Hubble modes, for which Eq. ( 37) applies.Using this result, we replace h s by ĥs in the energy-momentum tensor in Eq. ( 29) and compute its expectation value.The result is where we used the Wrosnkian condition, and the fact that the integration limit k > H restricts us to sub-Hubble modes.In the last line, we have taken the polarization sum (starting from the Bunch-Davies vacuum, λ ± are identical for both polarizations).Note that, regardless of λ ± , the final term of 1/2 makes Eq. ( 38) diverge for large k-this is the usual energy density vacuum divergence of quantum field theory.One can regularize the result by normal ordering the ladder operators in ρGW .However, this has to be done with the late-time ladder operators which annihilate the late-time Bunch-Davies vacuum.These are related to the original ladder operators âs k by a Bogoliubov transformation; for a detailed discussion, see e.g.Ref. [109].The regularized energy density becomes In practice, all of our modes of interest are highly excited with |λ − | ≫ 1, so that Eqs.(38) and ( 39) are approximately equal.In this limit, the vacuum contribution is negligible and the GWs are essentially classical.

Energy density scaling and the problem with kination
From Eq. ( 39), we see that the sub-Hubble GWs scale as radiation, with ρ GW ∝ a −4 , as expected for massless degrees of freedom.In cosmology with a standard expansion history, only a small amount of GWs are generated during inflation, and they always stay subdominant compared to the background radiation energy density.However, during kination, the background dilutes faster than radiation, and the gravitational wave fraction grows.The resulting gravitational wave spectrum is peaked and tends to either clash with bounds on the number of relativistic degrees of freedom during BBN or be hard to observe in gravitational wave experiments [58,59,60,61,62,63].In the following sections, we will demonstrate that adding a period of hyperkination helps with this issue, opening a wider parameter space for allowed GW spectra.
4 Analytical solution

Solving the background
Let us move on to solve the GW spectrum analytically.The first step is to solve the background dynamics, in particular the scale factor a, in the presence of radiation, as a function of the conformal time η.This provides us with a ′′ /a, allowing us to later solve the Mukhanov-Sasaki equation for the GW mode functions.
The scale factor evolves through different epochs during the cosmic history: inflation, hyperkination, kination, and radiation domination.The transitions between the epochs, assumed to be instantaneous, happen at conformal times η end (end of inflation and start of hyperkination), η kin (end of hyperkination and start of kination), and η reh (end of kination and start of radiation domination, i.e., reheating), which we will also solve in terms of the model parameters below.We use the same indices to refer to various variables evaluated at these times.We require the continuity of a(η) and its derivative at the transition times; between them, we solve a(η) from If we know how the Universe's energy density ρ scales in a, we can integrate and invert Eq. ( 40) to obtain a(η) epoch by epoch.We will normalize the scale factor so that and write a = e N , so that N counts the e-folds since the end of inflation.
For inflation, we assume a generic slow-roll inflationary phase, with the end of inflation η end < 0 determined by the usual condition where ϵ is the first Hubble slow-roll parameter.For the reader's benefit, we will express our GW mode functions as a function of ϵ, approximated to be constant.In the example spectra we consider in section 5, we work in the pure de Sitter limit ϵ = 0. To avoid clutter (and slightly abusing the notation), we will use H, as for pure de Sitter, to refer to the Hubble parameter at the end of inflation H end .
For hyperkination, we get ρ(N ) from Eq. ( 14), where ϕ ′ follows the first branch of Eq. ( 17) and we write the initial field velocity ϕ ′ 0 in terms of N hyp as explained below the equation.For kination and radiation domina-tion, we use the standard results ρ ∝ a −6 and ρ ∝ a −4 .With these, the full behaviour of the scale factor becomes For the hyperkination expression, we used which follows from Eq. ( 41) and the first line in Eq. ( 43).
We also used Eq. ( 20) with 3H 2 m 2 P = ρ end to eliminate α.For a long hyperkination period with N hyp ≳ 1, we can approximate the expression as where the right-hand-side is exactly the scale factor for a radiation-dominated universe, compare to the last line in Eq. (43).Note that the last approximation stops being valid at large times η ∼ e N hyp /H and one needs to use the middle expression instead.This is the case below, when we obtain an analytical estimate for η kin .
For kination and radiation domination, the constants in Eq. ( 43) are to be read off from the end values during the previous phase.Using Eq. ( 45), we have and , In practice, it is a good approximation to use Let us next estimate the conformal times for the rest of the transition points.We do this by solving an equation where a is expressed in two different ways, through Eq. ( 43) and through a condition related to our model parameters.
As a reminder, we define the beginning of kination as the time at which both addends inside the parenthesis in the energy density in Eq. ( 12) become equal.Since this happens at large times η ∼ e N hyp /H, we use the middle expression in Eq. ( 45) together with Eq. ( 19) to obtain so that The time of reheating η reh can be estimated by noting that the total energy density during kination scales as ρ ∝ a −6 , while that of the radiation scales as ρ r ∝ a −4 .Thus, the density parameter of radiation during kination scales as Ω r ∝ a 2 .By reheating, radiation is the dominant component, that is, so that where we used Ω kin r ≈ Ω end r , since the field and radiation redshift similarly during hyperkination, together with the approximation |η end | ≪ η kin ≪ η reh yielding a reh ≈ H √ 2η kin η reh from Eqs. ( 46) and (47).We also used Eq. ( 49) for a kin and Eq. ( 50) for η kin .

The gravitational wave mode functions
The next step is to obtain expressions for the GW mode functions.We proceed by matching the solutions and their derivatives at the transitions between epochs.To simplify the expressions, we do the matching in the super-Hubble limit, which gives an excellent approximation except for modes entering the horizon around the transitions.Our goal is to obtain the coefficients λ − (k) from Eq. (37) for each mode so that we can read off their asymptotic, sub-Hubble behaviour.We report the details of the somewhat technical calculations in Appendix C, while in the present section we simply give the main results, as well as a comparison between the analytical and numerical solutions in Fig. 3 (for details on the numerics see Appendix B).
We then get The constants η c can be read from the previous section, giving where we defined for kination Note that a ′′ = 0 during hyperkination.This feature is shared with the period of radiation domination, during which the spectrum is flat, a result that was originally derived in Ref. [110].Therefore we expect the peak from kination to be truncated by a secondary plateau.
With this, we can proceed to solve the Mukhanov-Sasaki equation (35).Making the change of variables x = kη (x = −kη during inflation when η < 0) and redefining the mode functions as v = √ xg, it can be recast as a Bessel equation the most general solution of which is given by where H (1) ν and H (2) ν are Hankel functions of the first and second kind respectively.Using the values of ν from above, the solutions during inflation, hyperkination, kination, and radiation domination become where we fixed the coefficients c 1,2 during inflation so that in the initial sub-Hubble regime, −kη ≫ 1, the mode functions obey the Bunch-Davies vacuum conditions in Eq. (36).The constants and phases in the other branches have been chosen so that the coefficients α ± , β ± , and γ ± correspond to the λ ± of Eq. ( 37) in the late sub-Hubble limit kη ≫ 1.Their values are fixed by requiring the continuity of v s k and its derivative at the transition times η end , η kin , and η reh .Matching the branches in the super-Hubble limit yields where and z reh ≃ η reh is z from (56) evaluated at η reh .For the scale-invariant case with ϵ → 0, f (ϵ) → 1, the moduli squared of the coefficients take the simplified forms Note that since we did the matchings at the super-Hubble limit, the expressions in Eqs. ( 60)-( 62) and ( 64) only apply for modes that are super-Hubble during the corresponding transition.To find the final behaviour of a mode, we take the last transition where this applies, track the following mode function from Eq. ( 59) to the sub-Hubble limit, where it takes the form in Eq. (37), and equate the α − , β − , or γ − with the coefficient λ − .Indeed, after a mode has settled to its asymptotic sub-Hubble behaviour, its evolution is trivial-redshfiting gently like radiation-and it won't be sensitive to further changes in the equation of state of the Universe.
From the Mukhanov-Sasaki solutions in Eq. ( 59) we can also deduce the metric perturbations h s k .Using the scale factor expressions, a ≃ Hη, a ≃ H √ 2η kin η, and a ≃ H η kin /(2η reh )η during hyperkination, kination, and radiation domination, respectively, and using Eqs.( 60)-( 62), we get where j 0 (kη) = π/(2kη)J 1/2 (kη) = sin kη/(kη) is a spherical Bessel function of the first kind and J 0 is a Bessel function of the first kind.For a comparison with the numerical solutions in the scale-invariant case, see Fig. 3.We do not include the inflationary metric perturbations in Eq. ( 65) as they do not simplify as nicely as the others.
Note that in the super-Hubble limit, all the expressions in Eq. ( 65) freeze to where the last one is the standard scale-invariant result.
Note that this result holds also for inflation.In principle, one can use this as an initial condition and solve the Klein-Gordon equation ( 28) to obtain Eq. ( 65) separately in each phase without the matching procedure described above 10 .One can then use Eq. ( 38) to obtain the unregularized GW energy density.We use this method in our numerical solutions.The expressions for α ± , β ± , and γ ± are still needed to regularize the integral in Eq. ( 38), and they are the conventional way to express the GW excitations in the literature.
5 Gravitational wave observations

Gravitational wave spectrum
We are finally in a position to calculate the spectral energy density of the primordial GW background.It is defined as where ρ is the total energy density of the Universe and ρ GW (k, η) is the contribution to the GW energy density from modes around k, given by Eq. ( 39) for the dominant, sub-Hubble modes.Here λ − is to be matched to α − , β − , or γ − as explained above.
To evaluate Eq. ( 67) at a specific time, we note that the radiation energy density can be written as (remember our normalization a end = 1) 10 In particular, GWs at the CMB scales stay frozen throughout the kination and hyperkination periods and are thus not affected by the non-standard background evolution.The same is true for the curvature perturbation R-see [111] for a linear treatment of R in a model with a non-standard kinetic sector, Appendix B of [94] for an application to Palatini R 2 models, and [112] for a general proof that R freezes at super-horizon scales. 11We neglect the change in the effective number of relativistic species contributing to the entropy g * S (T ) and to the energy density g * (T ).This introduces an additional mild scale dependence into the spectrum.For further details, we refer the reader to Ref. [113], and in particular to Fig. 4 therein.
In particular, using the current radiation temperature and total energy density, T 0 = 2.7 K = 0.23×10 −9 MeV and ρ 0 = 1.05 × 10 −120 m 4 P [114], we obtain ρ 0 r = 8.79 × 10 −125 m 4 P and Ω 0 r = 8.37 × 10 −5 .We use the index '0' to refer to quantities today.With this and the de Sitter limit results in Eq. ( 64) together with Eqs. ( 50) and ( 52), the GW spectrum today becomes Below, we will refer to the different branches as Ω rad GW , Ω kin GW , and Ω hyp GW .The boundary values are given by k = H at the end of inflation and at the transition times.Using Eqs. ( 46) and ( 47), we get where we approximated |η end | ≪ η kin ≪ η reh .
In our figures, we show the spectrum as a function of f , the GW frequency today.To relate f to our wavenumber 12 k, we use Eq. ( 68) and ρ = 3H 2 m 2 P , yielding An important frequency is the one that corresponds to BBN.It does not depend on the early expansion history, and we can solve it explicitly as where we used ρ BBN ≃ 3 × 10 −86 m 4 P .We present f BBN as a vertical dotted line in our graphs.
We show a comparison between the numerical and analytical spectra, for an example set of parameters, in Fig. 4. We see that the analytical expressions for the spectrum are very accurate.In Fig. 5 we present some example analytical spectra superimposed with the sensitivity curves for future GW experiments.
From Eqs. ( 70) and ( 71), we can straightforwardly understand the shape of the spectrum.The height of the first plateau, corresponding to hyperkination, is given by the combination H 2 /(Ω end r m 2 P ), i.e., the larger the energy density at the end of inflation and the smaller the reheating efficiency, the larger the energy density spectrum amplitude will be.The third free parameter of our theory, the number of e-folds of hyperkination N hyp , controls the length of the boosted spectrum; the longer the hyperkination period lasts, the more stretched the boosted spectrum is.In contrast, the height of the second plateau depends on H 2 /m 2 P , i.e., it depends on the energy scale at the end of inflation only, the standard result from a scenario with no period of kinetic domination, originally derived in Ref. [110].Both plateaus are connected via a region growing linearly with the frequency f , corresponding to the kination period.At large frequencies, the spectrum is cut off at the last mode to be excited by inflation.At small frequencies, there is no cutoff; the first line in Eq. ( 70) applies to all modes that re-enter during radiation domination.
Although it is easier to understand the shape of the spectrum in terms of N hyp , the free parameter in the action in Eq. ( 2) is α.For this reason, we present below our results regarding the parameter space of the theory in terms of α and not N hyp .The two are related by Eq. (20).For completeness, we present here the spectrum in terms of α, and with k replaced with f : where H=10 13 GeV 4: Analytical spectral energy density of the primordial GWs (dashed orange) and its numerical counterpart (full green).For details on the numerical solution, we refer the reader to Appendix B. The vertical dotted lines represent the frequencies associated with the start of kination, reheating and BBN, while the horizontal dashed line represents the BBN bound on the spectrum.The numerical spectral energy density is not well resolved at the largest frequencies because the modes re-entering the horizon right after inflation are never frozen as assumed in the code.This leads to the unphysical upslope around 10 11 Hz.The parameters used are N hyp = 15, Ω end r = 10 −10 and H = 10 13 GeV.
Note that the frequencies of the modes that cause the truncated peak, corresponding to hyperkination and kination, are always between f reh and f end , given by Eqs. ( 76) and (75), respectively.The specific values depend on the Hubble parameter at the end of inflation H, the density parameter of radiation at the end of inflation Ω end r and α.In order to give some indicative values, let us assume GUT scale inflation H ≃ 10 13 GeV and electroweak-scale reheating ρ(η reh ) ≃ (200 GeV) 4 , corresponding to Ω end r = 10 −10 .Changing α obviously leaves f end unchanged.In Table 1, we show f reh and f end for a few different α.Note that they are larger than f BBN , as they should be.Table 1: Values of the frequencies corresponding to reheating f reh and the end of inflation f end for different values of α, given that H = 10 13 GeV and Ω end r = 10 −10 .

Parameter space and detectability
In the present section, we put our model to the test and analyse the detectability of the generated spectrum of primordial GWs in the presence of a period of hyperkination after inflation.Since our analytical expression for the spectrum approximates very well its numerical counterpart, as can be seen from Fig. 4, we use it in order to compare with the sensitivity curves of various detectors, namely LISA [30,31,32], ET [115,116], LVK observing runs O3 and O5 [25,26,27,28,29], SKA [117], DECIGO [33,34,35] and BBO [36].For each of them, we run a scan over the parameter space {α, Ω end r , H}.The successful parameter space can be found in Fig. 7.
Before we describe how the parameter space scan is performed, we comment on some bounds that need to be imposed.First, BBN should happen during the period of radiation domination.In other words, at (and below) the frequency associated with BBN, the spectrum needs to be in its lower plateau, i.e., f reh > f BBN , where f reh is given by Eq. ( 76).This imposes a bound on the maximum value α can take.Solving for α in Eq. (76) gives GeV, 10 −12 } for the blue curve, {25, 7.9 × 10 11 GeV, 10 −9 } for the orange curve, {20, 7.9 × 10 10 GeV, 10 −5 } for the green curve and {29.5, 1.7 × 10 13 GeV, 10 −8 } for the red curve.We also show lines parallel (dashed gray) to the kination part of the spectrum.If not for the hyperkination period the spectra would violate the BBN bound.
Importantly, we note here that the specific value we use for f BBN in Eq. ( 73) comes from T BBN = 1 MeV.However, recent studies [69,102] have shown that the stiff era is restricted to occur at temperatures T > 2.5 MeV.This means that the value in Eq. ( 73) would become a factor of 2.5 larger, and the bound in Eq. ( 77) a factor of 0.026 smaller.However, given that the available parameter space for α spans many order of magnitude (see Fig. 7), this change does not affect our results appreciably.Nevertheless, the reader should keep in mind that our bound T > 1 MeV is an approximate one.
In addition, the GW energy density at BBN must be low enough not to disturb the standard results.Eqs.(67) and (68) give Ω BBN GW = Ω 0 GW /Ω 0 r , allowing us to translate the bound to into the GW energy density today, yielding [118] where h ≈ 0.7 is the dimensionless Hubble constant.In practice, however, for all detectors except LVK O5 and ET, this bound is irrelevant.Indeed, it is sufficient to impose that the hyperkination plateau be below the minimum of LVK O3, the region excluded by now by LVK, which is below the BBN bound.Note that for LVK O5 and ET there exists some parameter space where the hyperkination plateau is between both limits.We take this into account in the scans by showing the excluded region from LVK O3 in Fig. 6.There, for each value of H and Ω end r , we show the maximum value of α, labelled α max , below which the signal is not observationally excluded.
We can also impose an upper bound on the energy scale at the end of inflation.Using the slow-roll expression for the amplitude of the scalar power spectrum, we can write the Hubble parameter at CMB scales as where A s = 2.1×10 −9 [114] and r is the tensor-to-scalar ratio.The latest constraint on r is r < 0.036 [16].The energy scale at the end of inflation is always lower than at CMB scales, so Eq. ( 79) provides an upper bound on H at the end of inflation, Further, the plateau corresponding to radiation domination should be below the one corresponding to hyperkination, but this is not strictly guaranteed by our approximative spectrum if the kination period is short.To ensure this condition is satisfied, we impose see Eq. ( 74).
The logic for the parameter scan is as follows.We consider a grid in the (H, Ω end r ) plane, with the values of H lying in the interval [10 6 , 4.7 × 10 13 ] GeV and those of Ω end r lying in the interval [10 −20 , 0.81], both in steps of 0.5 in logarithmic units.Then, for each point in the grid, we find the minimum value α min , such that our spectrum is detectable by the specific experiment we are considering.Since the effect of increasing α (or, analogously, N hyp ) is to stretch the flat region corresponding to hyperkination, if a signal is detectable for α min , it will also be detectable for every α > α min .Note that for LVK O5 and ET, for a certain region in the (H, Ω end r ) plane, there is also a maximum value that α can take, see Fig. 6.This limitation exists only for values where the height of the hyperkination plateau is above the minimum of the LVK O3 sensitivity curve.
In order to determine whether a signal can be detected, we compute the power-law integrated curves (PLIC) [119] for each experiment.Then, for each set of parameters, we find the minimum α min such that the energy density spectrum is at least as large as the PLIC under consideration.An easy way to picture this procedure is to realise that the spectra with α = α min are tangent to the PLICs.Increasing α increases the length of the hyperkination plateau, so if the spectrum is tangent to a PLIC, it will be above it for some frequency range if α > α min .
In Fig. 5, we show some example spectra with a large enough SNR, superimposed with the sensitivity curves for all considered experiments.In the same figure, we also show a grid of lines with the same slope as Ω kin GW (f, η 0 ) to showcase how in a setup with inflation being followed by usual kination most of the signals would violate the BBN bound.Hyperkination fixes this by truncating the spectrum and introducing a new plateau at high frequencies.
We report the results of parameter space scans as contour plots in Fig. 7. There, for each pair (H, Ω end r ), we give the minimum α min such that the signal is detectable, for each experiment.We emphasize that the totality of the successful parameter space is contained in these figures.Besides the maximum value of H from Eq. ( 80), the parameter space is bounded at small H by the BBN timing condition f reh > f BBN , at small Ω end r by the BBN energy density condition in Eq. ( 78) and the LVK O3 exclusion bound, and at large Ω end r by the requirement that the higher hyperkination plateau must reach the lower end of the sensitivity band for the given experiment.
We conclude that there is ample parameter space to accommodate detectability by all experiments.Indeed, as can be seen from Fig. 7, for a Hubble parameter H ≲ 10 13 GeV, somewhat below the GUT scale, and a reheating efficiency in the range of 10 −15 ≲ Ω end r ≲ 10 −2 , which can be easily accommodated by a variety of reheating mechanisms [40,120,99,100,101,121,122], For each value of H and Ω end r , there is a maximum value for α, labelled α max , above which the signal is observationally excluded.
we can always find a detectable signal.We emphasize that the size of the parameter space is large, and there is no need for fine-tuning to obtain a detectable signal.Indeed, in Fig. 7 we report the minimum value α has to take in order for the signal to be detectable.However, any value of α larger than α min also leads to a detectable signal.
The value of α min is quite large for most experiments.This can be understood from Eq. (20).Indeed, we can find a lower bound on α min by taking the limit N hyp ≪ 1.It gives Using a GUT energy scale ρ end ∼ 10 −10 m 4 P , considering an almost non-existent period of hyperkination with N hyp = 0.1, we obtain a rough lower bound α min ≳ 10 10 .As soon as we have a larger N hyp , α min grows exponentially with it.This is in line with our latest work [80], where we study quintessential inflation with an action of the form in Eq. ( 1).There, we find α ∼ 10 10 for successful quintessential inflation, without considerable hyperkination.For each value of H and Ω end r , there is a minimum value for α, labelled α min , above which the signal is always detectable (minus the excluded region in Fig. 6 for LVK O5 and ET).

Discussion and conclusions
We have investigated the spectrum of primordial gravitational waves (GWs) generated by cosmic inflation in a model where after inflation but before reheating we have a period when the Universe is dominated by the kinetic energy density of the inflaton scalar field ϕ, when the field is characterised by both the usual quadratic kinetic term and also by a higher-order quartic kinetic term.This is natural in theories of quadratic R + αR 2 gravity in the Palatini formalism, where in the Einstein frame the quartic kinetic term is proportional to α, the coefficient of quadratic gravity.However, we can equally well envisage a k-inflation scenario where the kinetic term of the scalar field includes a term ∝ αX 2 , where X = 1 2 φ2 .This kinetically dominated period is divided into two parts.In the first part, the inflaton kinetic energy density is dominated by the higher-order kinetic term; a period which we call hyperkination.In the second part, the higher-order kinetic term becomes negligible and the inflaton kinetic energy density is dominated by the usual quadratic term; a period called kination.We have shown that, while kination is a stiff phase with barotropic parameter w = p/ρ = 1, as is well known, hyperkination is not; the barotropic parameter during hyperkination is that of radiation w = 1/3.As a result, the modes of inflation-generated primordial GWs which re-enter the horizon during hyperkination form a flat spectrum, in the same way as the modes which reenter the horizon after reheating, in the usual radiation era.However, during usual kination, the GW spectrum is not flat but the GW density parameter per frequency logarithmic interval is Ω GW (f ) ∝ f .This means that, for modes re-entering the horizon after inflation and before reheating, the GW signal is boosted.This boost corresponds to a truncated peak in the GW spectrum; truncated because the spectrum corresponding to hyperkination is flat but it can be of much larger amplitude than that corresponding to the eventual radiation era.Consequently, the period of kinetic domination (kination + hyperkination) can be made to last longer and the boosted spectrum to extend to lower frequencies without the danger of the production of excessive primordial GWs.In particular, the truncated spectrum can avoid the upper bound imposed by the requirement that Big Bang Nucleosynthesis (BBN) remains undisturbed.Thus, primordial GWs in all observable frequencies can be enhanced without a problem.
We have analytically and numerically studied thoroughly the inflationary production and the subsequent evolution of GW modes and obtained the resulting GW spectrum, linking it with the model parameters.The characteristic shape of the spectrum will be testable in the near future by forthcoming experiments, such as advanced LIGO-Virgo-KAGRA, LISA, DECIGO, BBO and ET, as depicted in Fig. 5.If observed, such a spectrum can provide insight into the underlying theory, such as the energy scale of inflation, the reheating efficiency and the coefficient α.The latter is directly related to the duration of the hyperkination phase.Indeed, when hyperkination lasts N hyp , then Eq. ( 20) suggests where ρ end = 3H 2 m 2 P is the energy density at the end of inflation, and H is the corresponding Hubble scale.Typically, inflation is at the scale of grand unification, which implies H 2 ∼ 10 −10 m 2 P .In this case, the above suggests that e N hyp ∼ 10 −3 α 1/4 , which means that Note that, in the usual Starobinsky R 2 inflation we have α = 1.1 × 10 9 .Such large values of alpha are nonperturbative, but this is no more a problem in our setup than it is in Starobinsky gravity.Important information can also be deduced by the amplitude of the truncated peak corresponding to hyperkination.Indeed, Eq. ( 70) suggests that the value of the GW spectrum on the hyperkination plateau is given by where Ω 0 r ≃ 10 −4 is the density parameter of radiation at present and Ω end r is the density parameter of radiation at the end of inflation, also called reheating efficiency, because the larger it is the sooner reheating takes place.As discussed, in order not to destabilise BBN, we need Ω hyp GW < 10 −6 .Thus, we obtain a lower bound on the reheating efficiency as Ω end r > (H/m P ) 2 .Typically for inflation we have H 2 ∼ 10 −10 m 2 P , which implies Ω end r > 10 −10 .In an effort to stay generic, we have not considered a specific mechanism for producing the radiation which eventually reheats the Universe.We note however, that a number of such mechanisms exist, such as instant preheating [40,120], curvaton reheating [121,122] or Ricci reheating [99,100,101] to name but some.It is even possible to avoid introducing additional degrees of freedom and consider that reheating occurs due to the dissipating properties of the inflaton field itself, as discussed in Ref. [103], where such processes become negligible after inflation.
Additional important information can be obtained by the observation of the frequency of the knee in the GW spectrum, shown in Figs. 4 and 5, which is given by f kin in Eq. ( 75).Combining this with Eq. ( 85), in the large N hyp limit, we obtain Where ρ 0 = 3H 2 0 m 2 P is the energy density of the Universe at present.Putting the numbers in the above, we find Observations might provide the values of the left-handside of the above, which means that α could be estimated provided H is known (e.g.H 2 ∼ 10 −10 m 2 P for inflation at the grand unified energy scale).
In Fig. 7 we display our findings with respect to observability by different missions, such as LVK 05, ET, BBO, LISA DECIGO and SKA.There, we show the minimum value α has to take in order for the spectrum to be detectable.Above this value, which we label α min , the spectrum is always detectable.We see that observability requires that the reheating efficiency is smaller the lower the inflation energy scale is (the lower H is). Also, the values of α min are larger for large inflationary energy scales.For LVK 05 and LISA we find that observability requires α min ∼ 10 30−60 , while for ET, BBO and DECIGO the numbers are smaller α min ∼ 10 10−50 .For the reheating efficiency, we find that observability requires that the density parameter of radiation at the end of inflation is Ω end GW ≳ 10 −16 , a value which may increase up to unity or so in the case of ET, BBO or DECIGO.Such a high reheating efficiency implies that the kinetic regime is very small or even non-existent (prompt reheating).This is possible because, the ET, BBO and DECIGO might be able to detect very faint signals at frequencies higher than LISA, which means that they could even marginally observe the flat GW spectrum generated by the usual radiation era (no kinetic epoch).This is why there is a region (for ET, BBO and DECIGO) when H is large (H ∼ 10 13 GeV) where suddenly α can be very small (or even zero).The parameter space for this is very small though.
We conclude that, with our mechanism, the observability of primordial GWs is much enhanced compared to traditional models.We obtained concrete predictions involving H, α and the reheating efficiency in the case the characteristic form of the GW spectrum-a truncated peak-is indeed observed.Observation of the primordial GW signal would not only confirm another prediction of cosmic inflation but would also be a tantalising hint towards the quantum nature of gravity, which is behind the assumption of the Bunch-Davies vacuum in Eq. (36).Forthcoming GW observations may reveal new and surprising details about the physics of inflation and fundamental physics in general.Our work serves to explore such a possibility.correct magnitude to satisfy the coincidence requirement.Inflation, however, as described above would not work.Indeed, the original hybrid inflation model of Ref. [123], which is characterised by the inflationary potential in Eq. (A.2), produces a blue spectral index for the scalar curvature perturbation.
As shown in Ref. [78], things change when we embed the above model in R + R 2 Palatini modified gravity.We assume that λ is small enough, such that the potential in Eq. (A.2) during inflation is V ≃ 1 2 m 2 φ 2 .Then, the inflationary plateau is due to the quadratic gravity term, which flattens the potential and creates the inflationary plateau with U inf ≃ m 4 P /4α as discussed in Sec.2.1.As mentioned, the scenario with λ = 0 = g was investigated in Ref. [78], which found that successful quintessential inflation in achieved if m ∼ 10 13 GeV and ( 12 g 2 M 2 µ 6 ) 1/8 ∼ 10 GeV, which means µ ∼ g −1/3 10 −4 GeV.The assumption g > m/M ∼ 10 −3 suggests that µ ≲ O(MeV).
In Ref. [78] it was shown that for successful quintessential inflation with this model we need α ∼ 10 8 , so that U 1/4 inf ∼ 10 16 GeV.The canonical inflaton field rolls down the Palatini inflationary plateau U inf until it triggers the GUT phase transition and sends the GUT Higgs field to its VEV.Then, the potential V is reduced drastically so that the system exits the Palatini plateau and U ≃ V .
The change of the expectation value of the GUT Higgs field χ at the phase transition not only terminates inflation but may also affect the value of α provided the latter depends on χ.Indeed, suppose that α = α(χ) = α 0 e κχ/M , (A.4) where κ = O(10) is a coefficient and α 0 ∼ 10 8 .Before the phase transition, χ = 0 and α = α 0 ∼ 10 8 .After the phase transition, χ = M ∼ 10 16 GeV and κχ/M ≲ 10 2 .As a consequence, α becomes huge.Indeed, for the range κ = 5 − 166 we find α ∼ 10 10-80 , which comfortably includes the values considered in Fig. 7.Note that α should not depend on the inflaton field, α ̸ = α(φ), because the latter changes substantially during kination and hyperkination, while α is taken to be constant.Finally, it must be pointed out that the period of hyperkination in the post-inflationary history would modify the treatment of Ref. [78] somewhat.As a result, the value of µ for successful coincidence might change, but this is beyond the scope of the present work.

Appendix B: Numerical solutions
To check our analytical results, we solve numerically the time evolution of the background composed of the field and radiation and the GW mode functions.The full set of equations reads 1 + 3α φ2 Many of the variables vary by orders of magnitude during cosmic evolution.To make numerics easier, we define new, rescaled variables x, y, and Z, a new time variable s, and a constant s 0 through where H 0 is the initial Hubble parameter.Definitions in Eq. (B.6) are chosen to ensure the new numerical quantities remain of order one throughout the computation.The equations of motion become where a circle over a variable indicates a derivative with respect to the new time variable s.
The initial conditions for the field velocity and fluid energy density are set as described in the text, engineered to match a desired end-of-inflation Hubble parameter H end , duration of hyperkination N hyp , and initial radiation energy density fraction Ω end r .We then follow their evolution from the end of inflation until the BBN temperature is reached, see Fig. 2. The gravitational wave modes are evolved from their frozen super-Hubble state in Eq. ( 66) starting somewhat before they re-enter the Hubble radius, until somewhat after the re-entry, after which they are taken to behave as radiation.To get the mode energy density, we use the first equation in Eq. ( 38)-as explained in the text, the error related to regularization is negligible for all relevant modes.Iterated over a number of modes, this produces the spectra in Fig. 4.
with the transition from inflation to hyperkination, which takes place at η end .During the hyperkination, the Mukhanov Sasaki equation reads (see Eqs. (35) and ( 55 We now take the super-Hubble (small argument) limit x end ≪ 1.Noting that the leading contributions come from the terms proportional to H Making the change of variables y ≡ k (η − η kin /2 + 1/H) (where y = kz in the notation of Eq. ( 56)) and redefining the mode functions as g = √ yv, this equation can be recast as a Bessel equation with ν = 0 (see Eq. ( 57)).Thus, the solution reads √ y e −iπ/4 β + (k)H where the overall constant and phase has been chosen such that the mode functions have a simple sub-Hubble (y ≫ 1) limit, as discussed below Eq. ( 59).We match this equation (and its derivative) with Eq. (C.9) (and its derivative) at time η kin , i.e., at 1 (y kin ) sin (kη kin ) . (C.25) We can now take the super-Hubble limit kη kin ≪ 1.Using kη kin = 2y kin , the term in brackets multiplying H The

Fig. 2 :
Fig.2: Left: Logarithm of the energy density of the Universe (full black), the field (dashed orange) and the background radiation fluid (dashed blue) as a function of the number of e-folds calculated from the end of inflation, obtained by numerically solving the system.Right: Barotropic parameter of the Universe from the same computation.The vertical dashed lines correspond to the start of kination, reheating, and the BBN.The parameters for both panels are N hyp = 15, Ω end r = 10 −10 and H = 10 13 GeV.

Fig. 3 :
Fig.3: Comparison between the analytical solution (solid blue lines) and its numerical counterpart (dashed orange) of the imaginary part of the mode functions h s k as a function of the elapsing number of e-folds when the mode enters the horizon during the hyperkination (top left), kination (top right) and radiation domination (bottom left) periods.The match is excellent, except when the wavenumber of the mode is comparable to the horizon size at a transition (bottom right).The vertical dashed lines represent the time of horizon crossing k = aH and the times at which kination starts N kin and reheating happens N reh .The parameters for all panels are N hyp = 15, Ω end r = 10 −10 and H = 10 13 GeV.

Fig. 6 :
Fig.6: Parameter space of the theory excluded by LVK O3.For each value of H and Ω end r , there is a maximum value for α, labelled α max , above which the signal is observationally excluded.

Fig. 7 :
Fig. 7: Parameter space of the theory for the minimum α such that the signal is detectable by LVK O5 (top left), ET (top middle), DECIGO (top right, BBO (bottom left), and LISA (bottom middle) and SKA (bottom right).For each value of H and Ω end r , there is a minimum value for α, labelled α min , above which the signal is always detectable (minus the excluded region in Fig.6for LVK O5 and ET).