Einstein-Cartan pseudoscalaron inflation

We study a class of early universe cosmological models based on Einstein-Cartan gravity and including a higher derivative term corresponding to a power of the Holst scalar curvature. The resulting effective action is basically given by General Relativity and an additional neutral pseudoscalar field (the pseudoscalaron), unequivocally related to the corresponding components of the torsion, that necessarily acquire a dynamics. The induced pseudoscalaron potential provides a realistic inflationary phase together with a very rich postinflationary epoch, resulting by the coupling of the pseudoscalaron to ordinary matter.

The theory of cosmological inflation [1][2][3][4][5][6] (for reviews, see [7][8][9][10][11]) provides a comprehensive explanation for the origin of the special initial conditions that support the standard Hot Big Bang (HBB) theory.These conditions include the flatness of the three-dimensional constant time hypersurfaces, the homogeneity and isotropy of the cosmic microwave background (CMB), the apparent absence of heavy relic particles such as magnetic monopoles, and the large (comoving) observable entropy.Additionally, inflation generates adiabatic, gaussian, and almost scale-invariant scalar metric perturbations.These perturbations are responsible for both matter inhomogeneities, that lead to the formation of the large-scale structures, and the temperature anisotropies observed in CMB photons.Furthermore, inflation also gives rise to tensor perturbations, i.e. gravitational waves (GW), that could be detected if the inflationary energy scale is sufficiently high [12][13][14].The inflationary stage can be triggered in various ways, although the simplest one is the single-field slow-roll scenario [15,16].In this class of models, it is assumed that below the Planck scale, a single, homogeneous, neutral (spin 0), minimally coupled and canonically normalized (pseudo)scalar field ϕ called the inflaton, plays the dominant role within the stress-energy tensor T µν of the Universe.Specifically, the inflaton field is characterized by a scalar potential V (ϕ), known as the inflationary potential, which exhibits an almost plateau-like region and a global vacuum.Initially, the inflaton field is misaligned from the minima and slowly moves through the flat region.As a consequence, the potential contribution dominates over the kinetic term ( φ2 ≪ V (ϕ)), and it mimicks the presence of a false vacuum or a (transient) effective cosmological constant, whose energy density contribution, V (ϕ) ∼ M 4  inf ≪ M 4 p , leads to an almost de Sitter expansion of the Universe.Once the inflaton crosses a total distance ∆ϕ [17][18][19][20][21][22] and reaches a slow-roll-breaking value, inflation ends.As a results, the kinetic term becomes important and the inflaton rapidly falls down to the global vacuum, around which it begins to oscillate.The inflaton field should also be coupled to the degrees of freedom of the Standard Model (SM) or of its hypothetical extensions (BSM), enabling the transfer of energy density stored in the field condensate to the SM (or BSM) sectors.The reheating of the Universe then occurs, resulting in a hot plasma of SM (or BSM) particles with the required (comoving) entropy and paving the way to the initial standard radiation dominated phase of the HBB cosmology (see [23][24][25][26][27][28][29] for perturbative reheating, [30][31][32][33][34][35][36][37] for nonperturbative reheating and [38][39][40][41][42] for reviews).
The depicted single-field slow-roll inflationary scenario should naturally arise within a more fundamental scheme, where both extensions of the SM of particle physics at a scale higher than the electroweak one and of General Relativity (GR) are consistently coupled in a full quantum theory.The only available examples of such theories are (Super)strings or M-theory [43,44].They cure the non-renormalizability [45,46] of GR and typically predict effective actions containing exactly additional exotic matter and gauge fields together with modified (super)gravities.Unfortunately, a complete (realistic) top-down model is still not available starting from ten or eleven dimensions, due to the huge number of vacua and to the technical difficulties of selecting in a clean and complete way a four dimensional preferred vacuum.In this paper, we thus take a (bottom-up) effective action point of view (see e.g.[47] for a general treatment), showing an example of a class of theories (developed in [48]) where the inflaton comes unequivocally from geometry.In the framework of an Einstein-Cartan [49][50][51][52] extension of GR (see also [53][54][55][56][57][58][59][60][61][62][63] for recent early universe developments) it is indeed possible to integrate out some components of the connection.Introducing a suitable non-linear term related to a power of the Holst curvature term [64][65][66], one gets a dynamical pseudoscalar field (the pseudoscalaron) equipped with a potential able to drive an inflationary phase and, possibly, the subsequent reheating of the Universe.The pseudoscalaron would not exist in a Riemannian geometry, and can be coupled to SM (or BSM) matter.Many other examples in the same class can be studied, giving rise to an interesting zoo of inflationary models, fully compatible with the current observational evidences.Moreover, the study of the rich post-inflationary epoch could also provide interesting windows towards an understanding of the reheating phase, leptogenesis, stochastic primordial GW background and Dark Matter.We will sketch a description of some of these issues, postponing a more detailed analysis to forthcoming papers.The plan of the paper is as follows.In Section II, we briefly review the Palatini approach to metric-affine extensions of GR [67][68][69][70][71][72][73], where the affine connection is promoted to be independent of the metric.We also review how to build a class of effective field theories where the torsion of the Einstein-Cartan connection is dynamical and classically equivalent to a pseudoscalar field, related to the so called Holst curvature term [64][65][66].Finally, we introduce our simple class of modified gravities, deriving the Lagrangian of the canonically normalized pseudoscalaron field.The pseudoscalaron potential is studied in detail in Section III, where we show how it is suited for inflation in a single-field slow-roll cosmological scenario.In Section IV, we numerically solve the inflationary equation of motion for the pseudoscalaron evolution, computing the corresponding Potential Slow-Roll Parameters (PSRP) and obtaining predictions for the scalar spectral index and the tensor-to-scalar ratio.We utilize the standard COBE normalization to constrain both the inflationary scale and the parameters of the scalar potential.In Section V, we provide an overview on the postinflationary phase for the whole class of models, stressing the dependence of the reheating epoch on the powers of the introduced non-linear Holst scalar curvature term.In particular, we discuss qualitative features of the models, leaving a quantitative deeper analysis to future publications.Section VI contains our final remarks and comments, with a summary of the results and also mentioning several open perspectives.
Throughout this paper, we use natural units with ℏ = c = 1.The reduced Planck mass is denoted as , where G N indicates the gravitational Newton constant.We adopt the "mostly plus" spacetime signature (−, +, +, +) for the four-dimensional Lorentzian spacetime metric.

II. EINSTEIN-CARTAN GRAVITY AND THE PSEUDOSCALARON FROM DYNAMICAL TORSION
The original Starobinsky cosmological model demonstrated that (trace anomaly) gravitational higher-order corrections could drive a successful inflationary phase.Subsequently, such proposal was revisited and it was soon realized that (i) inflation is basically controlled just by the R 2 correction and that (ii) this version is classically equivalent to GR plus an additional scalar field, the so-called scalaron, canonically coupled to GR and equipped with a nontrivial potential [74][75][76][77][78].This insight was also recognized in the supergravity context [79] and further extended, suggesting that effective modified gravity theories with a (phenomenological) lagrangian density of the form f (R) could correspond to the well known scalar-tensor theories [80][81][82].In an effective field theory approach that aims to describe the coupling of the Standard Model of particle physics to gravity, higher order terms also emerge quite naturally by quantum corrections1 or in the realm of UV (quantum) completions of GR, like Supergravities or (Super)string/Mtheory.Moreover, such extensions of GR are better described if the connection is promoted to be an independent field with respect to the metric.Indeed, the choice of (metric and torsionless) Levi-Civita connection in GR is a convenient option, basically related to the choice of the (natural) Einstein-Hilbert action.A metric-affine connection differs from the Levi-Civita connection by a tensor that we dubb (as in [48]) distorsion.Gravity theories based on a generic (linear) affine connection are known as "metric-affine" gravities [67][68][69][70][71][72][73].In this paper, we limit ourselves to Einstein-Cartan manifolds, where the connection has torsion but is metric compatible.The reason is that spinor fields coupled to gravity must also be included in the matter sector.Following the notations and conventions in [48], the distorsion is defined as where A ρ µ σ is the generic connection, while Γ ρ µ σ is its Levi-Civita component.Obviously, the torsion T µνρ is the antisymmetric part of the distorsion The curvature associated with A ρ µ σ is defined by and can be expressed in terms of C ρ µ σ as where R ρ µν σ is the "standard" (Levi-Civita) Riemann tensor.The curvature tensor can be contracted to provide the usual Ricci scalar curvature R ≡ R µν µν (5) and a pseudoscalar called the Holst invariant (see [64][65][66]) where ε µνρσ is the totally antisymmetric Levi-Civita tensor with √ −g ε 0123 = 1.Note that R ′ vanishes for C ρ µ σ = 0, namely when the connection is the distorsionless Levi-Civita one.This is the reason why in the standard formulation of GR R ′ is always absent.However, it plays a prominent role in the class of theories we are going to consider, where the distorsion is dynamical.As shown in [48], they can be described using an action of the form where with Φ we generically denote all fields that do not depend upon the distorsion and enter the action in combinations that respect the (global and local) symmetries present in the lagrangian, together with its scalar nature.It should be noticed that the α and β functions contain (possible) non-minimal couplings to the scalar and pseudoscalar curvatures, while ∆ is an arbitrary function of the indicated fields and curvatures that brings the non-linear terms.
It must be chosen appropriately (see ref. [48] for more details) in order to guarantee that field redefinitions do not take the model back to cases where the distorsion is non-vanishing but non-dynamical.Finally, Σ(Φ, DΦ, C) contains the "matter" fields.Notice that it depends on the distorsion through the covariant derivatives built out of the whole A connection, and thus on C.
In this paper, we focus on a simple class of theories that are exactly solvable and give rise to an interesting set of inflationary models where the inflaton can be identified with a pseudoscalar field representing exactly a pseudoscalar combination of the distorsion components, thus originating directly and unequivocally from the geometry of the underlying spacetime.To stress the inflationary scenario, we take Σ = 0 and select 2α(Φ) = M2 P (thus directly the "Einstein frame") and 4γβ(Φ) = M 2 P where γ is known as the Barbero-Immirzi parameter [83,84].Moreover, we choose ∆ depending solely by the pseudoscalar curvature and of the form 2 where p > 1 is a real number and ξ is a coupling constant with the dimensions of mass [m] 4−2p .As shown in [48], one may introduce an auxiliary pseudoscalar field z, in such a way that the considered model is classically equivalent to provided ∂ 2 ∆ ∂z 2 ̸ = 0. Indeed, the equation of motion of the auxiliary fields yields z = R ′ , giving back the previous model on-shell.It is now an easy algebraic exercise to decompose the distorsion into its irreducible components and to integrate it out.Defining the quantity it happens that its derivative sources the equations of motion of the vectorial and pseudovectorial components of the distorsion.In other words, on shell the action can be written as the sum of the Einstein-Hilbert action and the lagrangian density of the pseudoscalar field z, where R is the usual Levi Civita scalar curvature, the function K(z) reads as and the potential results The action in Eq.( 11) suggests that B(z) brings about the (non-canonical) kinetic term related to the pseudoscalaron z, that in turn certainly is not a ghost, since K(z) is always positive.As usual, a field redefinition allows to rewrite the pseudoscalar action in a canonical way.Firstly, we have to introduce a pseudoscalar field ϕ whose kinetic term is the standard one.This can be done letting The expression of K(z) in terms of B(z), allows us to establish the universal relation between the pseudoscalar field ϕ and B(z), which holds for the whole considered class of models, i.e.
where ϕ 0 is an integration constant.Secondly, we need to invert the previous relation to find z as a function of ϕ, in order to find the potential V = V [z(ϕ)].This procedure involves the solution of a complicated non-linear differential equation related to ∆(z) and its first derivative.In most cases, it is not possible to find an analytic solution.Fortunately, the simple form of our choice in Eq.( 8) allows to do that, as shown in the next Section.

III. INFLATIONARY MODEL
The assumption of Eq.( 8) allows to write the pseudoscalar sector of the action (11) in terms of the canonically normalized field ϕ.Indeed, the potential can be written as and Eq.( 15) can be explicitly inverted to give Therefore the final cosmological action is Here g is the determinant of the assumed Friedmann-Lemaitre-Robertson-Walker (FLRW) metric tensor g µν , with spacetime line element given by where dl is the line element of the three-dimensional spatial subspace, t is the cosmic time and a(t) is the dimensionless cosmic scale factor, which allows to define the standard Hubble rate H(t) = ȧ/a.Finally, the scalar potential comes out to be The evolution of the inflaton as a function of the cosmic time t is described by the standard Einstein-Klein-Gordon equations properly equipped with a set of initial conditions (for the field and its derivative) of the form where the time t * is meant to be close enough to the horizon crossing era of observed cosmological scales, typically around 50/60 e-folds before the end of inflation.The vacuum expectation value (vev) v of the field at the minimum of the potential is not constrained and typically could differ from ϕ = 0, as happens in models with a spontaneous symmetry breaking.To select an inflationary model with v = 0 and to interpret the pseudoscalaron ϕ as a particle oscillating around a flat Lorentz invariant vacuum, thus avoiding a Cosmological Constant term associated to the potential, it is customary to choose appropriately the integration constant ϕ 0 .In this respect, it is useful to introduce the combination where with γ the Barbero-Immirzi parameter.The condition V (ϕ) = 0 for ϕ = 0 implies that In conclusion, the scalar potential can be written in the convenient form where the scale of inflation can be identified with while the dependence on the field is encoded in In the spirit of an effective field theory approach, it is also useful to write the coupling strenght in the form , where ξ 0 is dimensionless.The sign of the Barbero-Immirzi parameter determines the direction of the slow-roll phase.Specifically, the slow-roll phase occurs for decreasing values of the inflaton field (i.e.φ < 0) for negative values of γ, while it occurs for increasing values of ϕ (i.e.φ > 0) for positive values of γ.In Fig. 1 both the global shape of inflationary potential and the vacuum geometry are shown, for γ = −10 −3 and some chosen pairs (p, M inf ).The Barbero-Immirzi parameter models the height and shape of the inflationary potential.It is straightforward to show that smaller values of γ (e.g.∼ |10 −2 |) implies a higher inflation scale and a shorter plateau.The parameter p controls the extension of the inflationary plateau, the asymptotics of the potential for large field values as well as the vacuum geometry.Indeed, as p increases, the plateau region becomes longer and the vacuum geometry becomes more and more cuspy.It is worth noting that the shape of the inflationary potential closely resembles that of some superstring-derived inflationary scenarios such as large volume inflation [86] and in particular fibre inflation [87][88][89] (see also [90] for an overview on (super)string-inspired inflationary cosmology).For example, the expanded form of the scalar potential for p = 2 and γ < 0 can be written as where b = 2 3 and the coefficients a i are positive and read The mathematical difference with the fibre inflation model lays on the exponential growth of V (ϕ) for large field values.
In the superstring case, the divergence of the potential is just driven by a single positive exponential term, while in the present case there are two exponential contributions related to the a 3 and a 4 coefficients.Exponential potentials of the form V ∼ e k ϕ/Mp are also interesting in string theory because they can provide some clues on the very onset of inflation, due to the "climbing" phenomenon.In certain orientifold models (see [91] for a review) corresponding to k ≥ 3/2 in d = 10, a range whose d = 4 counterpart would be for k ≥ √ 6, the dilaton can only emerge from the initial singularity "climbing" the exponential potential and reaching a turning point before a descent phase that can inject slow-roll inflation [92][93][94].With an exponential well, like in our case, within the interesting range p/(p − 1) > 3 (or p < 3/2), the behaviour near the singularity is even more complicated [95], with a chaotic sequence of bounces.
In the following sections, we will study the inflationary phase related to the potential in Eq. ( 27) for negative values of γ.

IV. INFLATIONARY PHASE
To start the analysis of the inflationary scenario, it is fundamental to characterize the slow-roll parameters and thus to compute the derivatives of the potential.One gets where and with The introduced parametrization easily provides the Potential Slow-Roll Parameters [16] purely as combinations of the functions f i (ϕ).In particular, at the leading order the first two PSRP are and namely As usual, the PSRP allow to infer the inflationary predictions.The (superhorizon solutions of the) scalar and tensor power spectra in terms of the horizon crossing epoch commonly read as and become, in terms of some horizon crossing inflaton value ϕ and The related scalar spectral index and tensor-to-scalar ratio at first order in the PSRP result or and The primary focus lies in computing such inflationary predictions for inflaton values corresponding to a time interval compatible with the number of e-folds before the end of inflation between 60 and 50, as usual.These inflaton values can be roughly or qualitatively determined by examining the final portion of the slow-roll plateau region of the scalar potential V (ϕ).Nevertheless, a preferred approach is to determine the inflaton values by exactly solving the inflationary equation of motion in Eq.( 21), reformulated using as independent "clock" variable the number of e-folds instead of the cosmic time.One gets where the ϵ N slow-roll parameter in terms of N is defined as and Eq.( 47) is a second-order differential equation, nonlinear with respect to the first derivative.It includes an effective potential term V ef f determined by the first slow-roll parameter ϵ V .Unfortunately, this equation cannot be solved analytically and requires numerical integration.The result is an inflationary field trajectory ϕ(N ) that provides inflationary predictions as functions of N , the number of e-folds before the end of inflation.Before proceeding with the numerical integration, it would be interesting to explore the slow-roll limit of the equation and attempt to find a close, albeit approximate, relation between ϕ and N .In the slow-roll approximation, where the inertial term and ϵ N are subdominant the equation can be approximated with a linear first order ODE whose integration allows to relate the number of e-folds before the end of inflation to ϵ V (ϕ) and to the field values by The number N depends on the value of ϕ at the end of inflation, defined by the condition Using Eq.( 38), one gets the two solutions where Of course, given the shape of the pseudoscalaron potential in Fig. 1, ϕ end must be chosen as the smaller, positive value among the two solutions in Eq.( 54) (since γ < 0).At the same time, the integrand of N (ϕ, ϕ end ) can be computed as Hence, the slow-roll solution N (ϕ) results The first contribution just depends on the (log)ratio of two hyperbolic cosine while the second contribution is the difference between the Gudermannian function computed at ϕ and ϕ end .In general, this equation cannot be exactly inverted in terms of the scalar field.However, it remains a very useful tool to compute pairs (N, ϕ) that can be used, for example as initial conditions to trigger the dynamics described by the Eq.( 47).For example, if γ ∼ −10 −3 , then for p = 2, ϕ(60) ∼ 5.45 while for p = 4 one gets ϕ(60) ∼ 4.99.Note that for values of the Barbero-Immirzi parameter |γ| ≫ 10 −1 Eq.( 57) can be properly approximated as It is possible to invert such expression and find where ω(N ) In Fig. 2 we provide the inflationary predictions derived through numerical integrations for N in the range [50,60].
In particular, we show the standard (n s , r) plane for the three cases of the parameter p, i.e. p = 2, 3, 4 and a set of Barbero-Immirzi γ.The white region represents the Planck constraints [96] for the scalar spectral index n s , namely n s = 0.9649 ± 0.0042 at 68% of confidence level (CL).The BICEP/Keck upper limit on the tensor-to-scalar ratio [97], i.e r 0.05 < 0.036 at 95% CL, is not visible.The computation shows that as the parameter p decreases, the inflationary predictions shift towards the central region and become more and more compatible with the Planck constraints.This means that both the scalar spectral index and the tensor-to-scalar ratio tend to decrease for a fixed value of γ.At the same time, the n s -r curves tend to cluster together as γ becomes larger (in modulus).For example, if N ∼ 55 and γ ∼ −10 −3 then n s ∼ 0.9650 and r ∼ 0.0035.On the other hand, if one considers the case p = 2 the scenario γ ∼ −10 −2 is practically disfavoured by the standard observations while for γ ∼ 0.5 × 10 −2 the compatibility with the observations is almost recovered.Instead, the general predicted amplitude of GW via the r parameter is of the order of ∼ 10 −3 and it is aligned with the plethora of models typically inspired by supersymmetry, supergravity or superstring theories.Let us also comment about the scale of inflation, namely the value of M inf .By requiring that the scalar power spectrum of Eq.( 41) coincides with the COBE measured P COBE S ∼ 2 × 10 −9 , we obtain  for some selected values of the Barbero-Immirzi parameter and the power p.The simplest case p = 2 provides a very large coupling ξ while for smaller values, the coupling values is higly suppressed.The inflationary scale could be easily of the order of ∼ 10 16 GeV.
The previous expression, fixed γ, p and N (or in other words the associated value of ϕ), allows to constrain the coupling parameter ξ that can be inferred to be It is important to note that such a coupling constant depends in a non-trivial way on the ratio of the Planck mass on a certain power of the Barbero-Immirzi parameter γ.In Tab.I, we present various scenarios for the inflationary energy scale M inf and the magnitude of the dimensionless coupling constant ξ 0 = ξM 2p−4 p .Notably, it can be observed that the coupling strength is significantly suppressed for p > 2, despite a slight variation in the inflation scale.This outcome directly stems from the power (p − 1): as p > 2, the dominant contribution becomes M −4 inf , leading to a damping effect on the value of ξ.

V. POSTINFLATIONARY PHASE
The inflationary phase typically dilutes all the preexisting energy and entropy densities, providing a very cold and almost empty Universe.However, as the pseudoscalaron reaches the corresponding slow-roll ending condition value ϕ end , it rapidly relaxes in the global vacuum and starts to oscillate, paving the way for the reheat of the Universe.The properties of the scalar field oscillations depend on the local vacuum geometry and are obviously damped by the background Universe expansion.In pseudoscalaron inflation there are basically two possibile scenarios.Let us start by the first one, already partially discussed in [85] and relative to the case in which the potential in Eq.( 27) exhibits a smooth global minimum.It happens in the case with p = 2, where the pseudoscalaron potential can be expanded around the (Lorentz invariant) vacuum as The pseudoscalaron acquires both a finite mass m 2 ϕ = d 2 V (ϕ)/dϕ 2 | ϕ=0 and higher order self-coupling contributions.The mass term naturally overwhelms the higher-order interactions and results larger than the Hubble rate H (m ϕ ≫ H).As a consequence, the (homogeneous) pseudoscalaron field begins to oscillate with a frequency ω ϕ compatible with the mass m ϕ and with a period of oscillation t ϕ shorter than the Hubble time scale In other words, a Hubble expansion tipically contains a huge number of field oscillations.At this stage of the Universe evolution, the system can be thought of as a condensate of a large number of heavy pseudoscalaron quanta with mass m ϕ and zero momentum.Thus, the inflaton can be described in terms of a perfect fluid with an averaged Equation-of-State [26][27][28] In order to properly reheat the Universe and to generate the corresponding (comoving) entropy density, the pseudoscalar inflaton must couple to the particles of the SM/BSM sectors, producing them in the open decay channels.Several decay products can appear in the Σ(Φ, DΦ, C) part of the Lagrangian of Eq. (7).For instance, one may have spinor fields minimally coupled to Einstein-Cartan gravity.Additional effective couplings among matter fields, that respect the symmetries of the total lagrangian and that can be generated by quantum corrections are also allowed.Among them, quite interesting are effective non-minimal couplings of scalar fields to the scalar curvatures as well as effective couplings of (abelian) gauge fields to the vectorial and pseudovectorial components of the distorsion.The former can give rise, for instance, to decays of the inflaton directly to the Higgs SM sector [85].The latter allow mixings of abelian gauge fields to vector or pseudovector fields, paving the way to an interpretation of these new dynamical components of the distorsion as dark photons [48].The interactions with ordinary (or dark) matter are very model-dependent, being connected to the way they mix to the "visible" photon and to the form of the gravity sector.Some effective interaction terms of the expansion can be guessed, at least at the lowest order.For instance, a non-minimally coupled additional scalar field χ can be considered.It means to have in Eq.( 7) the α term no-longer constant but function of χ.At the lowest order it results where g χ is a dimensionless coupling constant.It results in an interaction term in the effective Lagrangian of the form Analogously, for a Dirac fermion ψ minimally coupled the interaction term is of the form The exact coefficients can be obtained by integrating out the non-dynamical part of the connection, as before.Notice that the coupling strenghts c ϕii are typically functions of the Barbero-Immirzi parameter γ.The inflaton energy (density) conversion process is characterized by a perturbative phase, sometimes preceded by a preliminary nonperturbative regime, although the details strongly depend on the model parameters as well as on the couplings c ϕii .
In the case of a negligible or absent nonperturbative regime, the inflaton decay and the related particle production is just a single-particle (non collective) process, basically driven by the decay rates of the inflaton to the daughter particles Γ ϕ→i1,...,in .Such decay amplitudes typically depend on the inflaton mass, while the total decay amplitude is the sum over all the possible channels.As is well known, the impact of this process on the inflaton oscillations about the vacuum can be modeled by introducing a phenomenological friction term ∼ Γ ϕ φ in the inflaton equation of motion, that represents a further source of damping in addition to the standard Hubble friction.It is important to stress that the computed total decay amplitude Γ ϕ is crucial to establish the out-of-equilibrium decay of the pressureless pseudoscalaron (average) energy density ρ ϕ and the corresponding growing of the (SM/BSM) radiation one ρ r , via the usual Einstein-Boltzmann equations [26][27][28][29] ρϕ + 3Hγ ϕ ρ ϕ = −γ ϕ Γ ϕ ρ ϕ (67) ρr + 4Hρ r = +γ ϕ Γ ϕ ρ ϕ (68) where γ ϕ = 1 + w ϕ and w ϕ = 0 in the present case.The Hubble parameter is defined by the Freedman equation while the initial conditions are This evolution is characterized by a final reheating temperature that notoriously scales as T reh ∼ M p Γ ϕ .As qualitatively discussed in [85], the decay amplitude of the pseudoscalaron to fermion is proportional to ∼ m ϕ m 2 ψ /M 2 p .Therefore a quite large fermion mass (e.g m ψ ≲ m ϕ ) is required in order to reach a reheating temperature well above the electoweak scale.The SM does not contain such heavy fermions and an interesting possibility would be to insert in the matter sector, in addition to the SM fields, sterile Right-Handed Neutrinos (RHN), that naturally possess such large masses.Moreover, RHN could also work as the crucial ingredients to get nonthermal leptogenesis in this scenario.Clearly, a decay channel of the pseudoscalaron to non-minimally coupled scalars like that in Eq.( 64) would be sufficient to reheat the Universe in a more efficient way to the needed temperature.Indeed, the decay rate is independent of m χ and proportional to ∼ m 3 ϕ /M 2 p .As noted in [85], the scalar field χ could naturally be identified with the SM Higgs.Thus, the presence of only the distorsion in the gravity sector is sufficient to get inflation and reheating, without the necessity of additional BSM degrees of freedom.
Adding a minimally coupled scalar is also a possibility.However, there is no direct interaction with the pseudoscalaron, since the distorsion does not enter the covariant kinetic term of the additional scalar.Rather, the two scalars would naturally mix giving rise, likely, to a multifield inflation.
All the details missed or just mentioned in this brief description will be presented in a separate publication, where a quantitative and complete description of perturbative reheating and Leptogenesis in the proposed scenario will be given.Let us now briefly discuss the second class of scenarios, related to the case where p > 2. As one can immediately recognize from Eq.( 27), the global minimum around ϕ = 0 is not smooth, being of the form Therefore, the scalar potential exhibits a cuspy behaviour.Similar potentials occur in many models, like in k-inflation [98] and in string theory, e.g. in axion monodromy models and in flux compactifications [99][100][101][102].What happens is the following: first of all, the field oscillates around the minimum, but clearly the oscillations cease to be harmonic and the (would be) mass term is divergent.Fig. 3 shows an example of pseudoscalaron oscillation in the usual Minkowski spacetime (upper panel) and in the realistic expanding postinflationary Universe (lower panel).In the first case, the pseudoscalaron starts from a reference plateau field value and falls in the global vacuum.Here, it exhibits an asymmetric oscillation that is quite regular in the standard (almost) quadratic scenario, p = 2.However, as p > 2, the cuspy potential leads to a nontrivial oscillation characterized by a natural damping.In the second case, the Hubble friction plays a crucial role in suppressing the oscillation structures of the previous scenario, although one can always appreciate the difference between the quadratic and cuspy cases.In the p > 2 framework, the pseudoscalaron can still be treated as a perfect fluid that, however, is equipped with a highly nontrivial EoS given by This quantity is generically different from zero when p ̸ = 2 and in the limit of p ≫ 2, the EoS parameter takes the value −1/3, mimicking the same effect of a linear potential V (ϕ) ∼ |ϕ| around the vacuum.Certainly, the inflaton is still able to transfer energy to the decay products.Indeed, as shown in [103,104], cuspy potentials naturally lead to a so-called preheating phase, where the (highly efficient) inflaton energy density conversion process is realized through an exponential growth of (the occupation number of) the pseudoscalaron decay product modes.In other words, collective phenomena lead to a very fast reheating of the Universe because of Bose condensation effects triggered by a decay product phase space that is sufficiently populated.Preheating in the form of parametric resonance can also be treated when the involved couplings are large (see, for instance, [32]).In addition, preheating in cuspy potentials can also give rise to a copious formation of oscillons that source an additional stochastic background of primordial gravitational waves (GW) [103,104].Such GW peaks could be visible in the future observing run of the LIGO-Virgo experiment [105].To summarize, the post-reheating phase of the considered models with p > 2 is quite rich and interesting and it deserves, in our opinion, a deeper and more complete analysis.Again, we will report specific and detailed predictions on preheating (and related GW generation) in separate, forthcoming publications.

VI. CONCLUSIONS AND PROSPECTS
Fundamental theories containing (extensions of) the SM coupled to gravity, like superstrings or M-theory, typically cure the non-renormalizability [45,46] of GR and predict (Einstein-Cartan) modified gravities coupled to SM or BSM fields, together with a plethora of additional sectors.Pure Einstein-Hilbert GR is typically accompanied by quantum corrections, due to string loops and higher derivative corrections, related to massive string states.In a bottom-up approach, it is possible to include corrections in an effective action consisting of a power series in curvature terms weighted by corresponding powers of the Planck mass.In this paper, we have discussed a class of models (a subclass of those introduced in [48]) where the corrections to GR are unequivocally linked to the geometry of the (non-Riemannian) spacetime.In particular, we have considered Einstein-Cartan gravities with dynamical distorsion and higher derivative terms given by powers of the so-called Holst scalar curvature that, as known, vanishes in a Riemannian geometry.The resulting effective field theory is equivalent to GR coupled to a pseudoscalar field that can naturally drive a single-field slow-roll inflationary phase fully compatible with the current observational evidences.In addition, compatible couplings to ordinary matter provide a highly non-trivial reheating phase that depends on the power p of the (non-linear) Holst term (see Eq.( 8)).Indeed for p = 2 the potential is smooth around the minimum, giving rise to a standard field oscillation that can source both a preheating phase and/or a (subsequent) perturbative reheating phase driven by the decay rates of the inflaton to SM fields or hypothetical BSM fields.On the other hand, for p > 2 the potential is cuspy around the minimum and the standard, almost quadratic oscillation phenomenon, is absent.The inflaton energy density conversion can be realized via a rich preheating phase with a possible production of GW via oscillon dynamics.Moreover, the coupling of the pseudoscalaron to heavy RHN and their resulting decay can pave the way for nonthermal leptogenesis and the generation of the observed baryon asymmetry of the Universe.We have described in detail the class of involved models and their inflationary slow-roll phase.We have also depicted, in a qualitative way, the postinflationary phase, stressing the involved and peculiar interesting properties.The richness of the postinflationary phase requires and deserves a careful and deeper analysis.We leave it to separate, forthcoming, publications [106,107].Finally, It would be challenging to discover string theory setups where similar potentials are realized and give rise to the "climbing" phenomenon, as explained at the end of Section III.It would also be important to search for a new class of Einstein-Cartan gravities that generalize those proposed in this paper to models that are scale-invariant at any coupling, along the lines of [53] as well as [108,109].

Figure 1 .
Figure 1.First panel: Pseudoscalaron potential normalized to a chosen reference scale M 4 ref ∼ 10 62 GeV 4 .The parameter p controls the extension of the inflationary plateau.In particular, smaller values of p favor a potential uphill climbing for smaller values of the scalar field while larger values of p tend to suppress the potential uphill climbing, which begins at larger field values.Second panel: Pseudoscalaron potential around the vacuum.The parameter p also governs the vacuum geometry.It is evident that when p increases, the anharmonic terms around the vacuum become significant also for |ϕ| ≪ Mp.Moreover, for p > 2 the potential is cuspy around the minimum.

Figure 2 .
Figure 2. Observational constraints on pseudoscalaron inflation for three different values of the parameter p and a set of values for the Barbero-Immirzi parameter.The case p = 2 shows that large values (in modulus) of the γ parameter are heavily disfavoured while smaller values turns to be compatible with the Planck constraints.As p gets larger, large values of the Barbero-Immirzi parameters tend to become partially compatible.

Figure 3 .
Figure 3. Prototype pseudoscalaron oscillation in a Minkowski spacetime (upper panel) and in the more realistic postinflationary Universe (lower panel).The decay of the pseudoscalaron to daughter particles is neglected.