Some like it hot: $R^2$ heals Higgs inflation, but does not cool it

Strong coupling in Higgs inflation at high energies hinders a joint description of inflation, reheating and low-energy dynamics. The situation may be improved with a proper UV completion of the model. A well-defined self-consistent way is to introduce an $R^2$-term into the action. In this modified model the strong coupling scale returns back to the Planck scale, which justifies the use of the perturbative methods in studies of the model dynamics after inflation. We investigate the reheating of the post-inflationary Universe, which involves two highly anharmonic oscillators strongly interacting with each other: homogeneous Higgs field and scalaron. We observe that in interesting regions of model parameter space these oscillations make longitudinal components of the weak gauge bosons tachyonic, triggering instant preheating at timescales much shorter than the Hubble time. The weak gauge bosons are heavy and decay promptly into light Standard Model particles, ensuring the onset of the radiation domination era right after inflation.


Introduction
A large variety of inflationary models lead to dynamics which make the Universe spatially flat and homogeneous, simultaneously producing matter and gravitational perturbations consistent with present observations of the cosmic microwave background and large scale structure [1][2][3]. Naturally, viable inflationary models with additional signatures deserve special attention. In particular, independent tests of the inflationary dynamics are provided in models with (a part of the) inflaton sector playing some role in late-time cosmology or lowenergy particle phenomenology, see e.g. [4][5][6][7][8]. Among these models, Higgs-driven inflation [9] is in a unique position since the main component in this inflationary model has been discovered [10,11] and its properties were extensively studied at LHC [12].
While very appealing indeed, the concept of Higgs inflation has unresolved issues both on phenomenological and theoretical sides. The phenomenological tension arises due to the Standard Model (SM) parameters -the Higgs boson mass, top quark mass, strong and weak gauge couplings -whose central values measured at the electroweak scale imply that quantum corrections to the Higgs boson self-coupling make it negative for large values of the Higgs field [13][14][15][16]. Thus, its energy density becomes negative, causing problems for cosmology in general [17] and for implementation of plain Higgs inflation.
The theoretical problems are associated with the large dimensionless constant of the nonminimal coupling to gravity. First, the large coupling constant makes the theory strongly coupled much below the Planck scale [18,19]. Second, the model is non-renormalizable because of the non-minimal coupling to gravity, which decouples the parameters of the small and large field potentials at the quantum level. Hence, the cosmological observables (amplitudes and tilts of the perturbation spectra) and low-energy observables to be measured in particle physics experiments are not related to each other, which prevents direct tests of the model.
While one may argue that the true values of the relevant physical parameters are actually at about 2σ off their presently accepted central values, making the Higgs potential positive all the way to the Planck scale, the large coupling constant is the key ingredient of the model and cannot be changed at will. However, it was argued [20] that the strong coupling scale grows with the Higgs field value. Recall that the first investigations of the post-inflationary dynamics and reheating of the Universe [21,22] yield a consistent estimate of the reheating temperature, which demonstrates that from the inflation till the radiation domination stage the system dynamics never exhibit the strong coupling behaviour. It was also shown [23] that the presence of non-renormalizable operators suppressed by the field-dependent scale has no impact on the inflationary predictions for the spectral parameters of scalar and tensor perturbations.
However, the first estimates of the reheating in Refs. [21,22] were incomplete. The crucial ingredient -the evolution of longitudinal components of the weak gauge bosonswas missed there. The dependence of this component on the external Higgs field (inflaton) turns out to be spiky [24][25][26], leading to violent production of longitudinal modes of the weak gauge bosons and very rapid reheating of the Universe. This finding is dangerous due to the strong coupling problem, since the reheating dynamics happen right inside the strong coupling domain, making any analysis unreliable. At the same time, the cosmological predictions of an inflationary model depend on the evolution of the Universe after inflation and the reheating temperature, so knowledge of the proper reheating dynamics is essential for precise calculation of the inflationary observables.
All these problems ask for a modification of Higgs inflation capable of keeping the model inside the weak coupling regime from inflation, through preheating, till the onset of the hot stage. Such a modification has been recently suggested [27,28] with a key ingredient -an R 2 term 1 added to the model Lagrangian. 2 Detailed studies of the inflationary evolution and perturbative unitarity in the scalar, gravity and gauge sectors revealed a region in the model parameter space [28] where the model remains weakly coupled up to the Planck scale and its inflationary dynamics are similar to the original Higgs inflation. These conditions open a possibility to test this cosmological model directly in particle physics experiments.
Here we refine the predictions of Higgs inflation UV-completed by the R 2 -term. Namely, we study the post-inflationary dynamics of the model and reheating of the Universe. In the interesting range of the model parameter space we observe instant generation of radiation due to the tachyonic instability of the longitudinal components of the weak bosons 3 and their subsequent decays into light SM particles. This tachyonic instability is effective only for particular dynamics at the moments of the scalaron field crossing zero, which correspond to the vicinities of the bifurcation points between two branches of the solutions. Efficient production happens at the first scalaron zero crossing only for particular values of the model parameters. However, if the conditions for the tachyonic preheating are not satisfied at the first scalaron zero crossing, they will be satisfied at a later one. Therefore, if the production does not transfer all the background energy into radiation at the first crossing, this will still happen over several scalaron oscillations. As far as the scalaron oscillates with frequency much exceeding the Hubble rate, the whole reheating process happens instantly from a cosmological viewpoint. This picture is valid for the region of the cosmologically viable model parameters sufficiently close to the Higgs like limit [27,28] of the R 2 -healed Higgs inflation. Thus, while the preheating in the model proceeds through generation of the weak gauge bosons, as in Higgs inflation, the dynamics responsible for reheating are rather different. In particular, we confirm that the spike in the mass of the longitudinal gauge bosons [24][25][26] does not produce the weak bosons in the amount sufficient for reheating [30].
It is worth noting that the current description does not take into account the details of the backreaction of the produced particles on the background. In particular, it is assumed that, as far as the production of the weak bosons is sourced by the Higgs field at short timescales, there is no immediate energy transfer from the (slower) oscillations of the scalaron background. First, this observation by itself requires that reheating happens during several scalaron oscillations, unless the model parameters correspond to a very close vicinity of the tachyonic bifurcation at the first scalaron zero crossing. Second, even in the most optimistic case the backreaction must be taken into account because the intensive energy flow from the homogeneous mode to particles changes the scalar field dynamics and consequently the production process itself. Also, there is additional energy drain from the Higgs oscillations away from scalaron zero crossing due to parametric resonance, which is subdominant for direct reheating process, but changes the details of background homogeneous field evolution. All this issues require detailed investigation, which we leave for the further study. However, given the dramatic swiftness of the process we expect the instant preheating of the Universe right after inflation to be the general model prediction in the Higgs-like region of the model parameters.
The paper is organised as follows. In Sec. 2 we describe the model and recall its inflationary dynamics. We present the equations of motion for homogeneous scalar fields and reduce them for the small field limit suitable for investigation of the background evolution after inflation. Sec. 3 contains equations for inhomogeneous modes in the scalar and gauge boson sectors. Sec. 4 describes the particle production. The obtained numerical results are discussed in Sec. 5. We conclude and finally present the predictions for the cosmological parameters in Sec. 6.

Model Lagrangian and evolution of scalar homogeneous modes
Higgs inflation, augmented with a term quadratic in curvature (Ricci) scalar R, is described by the action where ξ, β and λ are dimensionless real positive couplings, h is the Higgs field in the unitary gauge (and we neglect its present non-zero vacuum expectation value irrelevant for the largefield dynamics), g is the determinant of the metric chosen as ds 2 = dt 2 − a 2 (t)dx 2 , and the reduced Planck mass M P is defined via Planck mass M P l as M 2 P = M 2 P l /(8π). Upon the Weyl transformation g µν → g µν × e √ 2 3 φ M P with real scalar function φ(x) the action (1) takes form of the Einstein-Hilbert for gravity and two coupled scalars -Higgs h and scalaron φ, This form explicitly shows the absence of any strong coupling problems in the scalar sector right up to the Planck scale if the model parameters obey the inequality [27,28] β ξ 2 4π .
To describe the model dynamics during inflation another form of the action is more suitable. Changing the variables (h, φ) → (H, Φ) according to one arrives at the following Lagrangian in the scalar sector The model exhibits effectively single-field inflation, whose trajectory in the scalar field space is defined by nullifying the second term in the potential (5), see Refs. [27,28] for details. The matter power spectrum normalization on the the Planck measurements [31] fixes the model parameters as [27] β + ξ 2 λ ≈ 2 × 10 9 .
The second term above must dominate to give the same predictions for cosmological perturbations as in the original Higgs inflation, which together with (3) constrain the viable range of model parameters as The lower end of this range corresponds to "Higgs-like" behaviour, while the upper end gets closer to R 2 inflation [24,28]. Recall that we treat λ as a positive parameter taking its natural value of 10 −3 -10 −2 . Hence the tilt of scalar power spectrum n s − 1 and the tensor-to-scalar ratio r are as for pure-Higgs inflation [9], with N e = 50-60 being a number of e-foldings until the end of inflation. This number depends slightly on the reheating temperature, and our main task in this paper is to get a reliable estimate of this quantity and hence to refine the prediction (8) by obtaining an exact value for N e . At the onset of inflation Φ exceeds H in value. The latter remains almost constant, but then starts to evolve smoothly, and at N e 60 before the end of inflation both fields crawl along the inflationary valley, where the second term in scalar potential (5) is zero. In the transverse direction the scalar potential is always very steep, preventing the production of undesirable isocurvature perturbations [27]. The equations of motion for the homogeneous Higgs H 0 (t) and scalaron Φ 0 (t) fields follow from the action with Lagrangian (4), (5): Here, V H and V Φ stand for the potential derivatives with respect to corresponding fields, and By the end of inflation both scalar homogeneous fields are sub-Planckian [28], and actively participate in the model dynamics. The Hubble scale at the end of inflation can be estimated from the analysis of Ref. [24] as First, one can check that the smallness of H 0 (12) guarantees that terms following from the non-canonical kinetic term in (4) can be neglected in equation of motions (9) (10). Then, in the small-field regime, the scalar potential (5) can be approximated as polynomial of the fourth order in fields, Hereafter we take the leading order terms in ξ 1. Right after inflation the potential (14) refers to a system of two coupled and highly anharmonic oscillators, whose initial effective frequencies are not far from the Planck scale, but rapidly decrease as the Universe expands and the field amplitudes fall down. The scalar potential is symmetric with respect to a change of the sign of the Higgs field, but contains terms odd in Φ. The last term in the first line of (14) is of the special interest, as it can give a large negative contribution (dominating with ξ 1 over the negative contribution of the last term in second line of (14)). Note in passing that the first two terms in the first line turn into the SM Higgs selfcoupling and scalaron mass terms, which dominate the potential at very small fields. One recognises that for our range of parameters (7), not λ, but mostly ξ 2 /β defines the Higgs self-coupling constant at large sub-Planckian fields. However, below the energy scale of the scalaron mass Φ must be integrated out at the quantum level. As a result, the last term in the first line of (14) provides a contribution to the Higgs self-coupling which exactly cancels the ξ 2 /β-part below this scale, making the parameter λ solely responsible for the Higgs selfcoupling at low energies, as it must be in the SM. In principle, to verify the model, one should independently measure all the parameters in particle collisions -the Higgs self-coupling (or mass, which has been already done), scalaron mass M P / √ 3β, and scalaron-Higgs interaction (which would allow to measure ξ). Then, independently, CMB observations give the relation (6), which could then be tested against the measurements in particle collisions.
Let us characterise the motion of the uniform background fields Φ 0 and H 0 after inflation. Typically one has |Φ 0 | > |H 0 |, and the timescale corresponding to the oscillations in Φ 0 is much longer than for H 0 direction. Therefore one can split the dynamics into periods with positive and negative Φ 0 .
For Φ 0 < 0 all terms in the potential (14) are positive, and Φ 0 makes half-oscillation with a typical frequency not lower than the scalaron mass, which follows from the second term of (14). At the same time the field H 0 oscillates about the origin with frequency depending on Φ 0 . The oscillation type (harmonic or anharmonic) depends on the Higgs amplitude and Φ 0 -i.e. which of the terms in (14) dominates. In the harmonic regime, one can read the frequency from (14) as Note, that since the Higgs evolution timescale is much shorter than that of the scalaron, the Higgs contribution can be averaged over its oscillations, further increasing the scalaron effective frequency. The case of positive Φ 0 is the most interesting for particle production. Here the fields evolve along the valley where the second term in potential (5) is zero (Φ 2 0 and Φ 0 H 2 0 terms in (14) cancel each other), where in the last equality we used the assumption (7). The fields oscillate around this trajectory, which departs from the origin, evolves towards large fields, then stops and returns back to the next zero crossing of the scalaron. The typical scalaron timescale of the evolution along the valley and back can be estimated by replacing the H 0 in the quartic potential (14) by its value at the bottom of the valley H min and calculating the effective quadratic part in Φ 0 . This gives for the scalaron frequency squared in the harmonic regime which is lower, see (7), than that at Φ 0 < 0, limited from below by the scalaron frequency in R 2 -inflation [32], M 2 P /(3β), and actually corresponds better to the frequency of oscillations in pure Higgs inflation, see Refs. [21,22]. At large Φ 0 the higher order terms in (14) can make the oscillations slightly anharmonic. Finally, the frequency of oscillations of H, transverse to the valley, can be found from the second derivative of the potential (14) V HH calculated in a point along the valley floor (16) Note, this frequency parametrically coincides with that at negative Φ 0 (15). One observes, that at Φ 0 > 0 there are two equivalent valleys in the potential, for positive and negative Higgs values, H 0 > 0 and H 0 < 0. When the scalaron Φ 0 returns from negative steep potential to the positive valleys, the system can end up in either of them. The choice depends on the number of oscillations that H 0 made while the scalaron field was completing its half oscillation at Φ 0 < 0. This ratio is determined by the model parameters only (values of λ, ξ, and β), but in practice is impossible to be estimated analytically, as far as it depends on the nonlinearities of the potential and on the amplitudes of the fields in the given oscillation. Indeed, while for the first scalaron zero crossing after the end of inflation the "choice" of the valley can be found analytically given the smooth ending of inflation, for further oscillations the situation is more complicated, as far as the system approaches zero not simply along the valley (16), but with some oscillations in the H 0 direction.
Note, that there is a bifurcation point, corresponding to the system coming from Φ 0 but bouncing from zero back either into the same valley, or into the valley with the opposite sign of H 0 . This situation is rather peculiar -the exact separation between these two trajectories corresponds to unstable motion with H 0 = 0 and growing Φ 0 > 0 -or the fields "climbing" up along the middle of the barrier separating the H 0 < 0 and H 0 > 0 valleys. It can be immediately seen that for such a configuration the massive term for the Higgs field in (14) (third term) is tachyonic. This is exactly the source of the tachyonic instability that we study in detail in the next section. The closer the trajectory of the background fields is to the bifurcation situation, the longer the system stays in the tachyonic region, leading to more efficient decay of the background into inhomogeneous modes. Also, bifurcation points correspond to the maximum transfer of the energy from the motion in the Φ 0 direction into the oscillations in the H 0 direction.
We conclude the study of homogeneous field dynamics with the observation that in this oscillatory picture, the Friedman equation shows that for (12) the Hubble timescale is much larger than a period of Φ 0 oscillation, which in turn is much larger than a period of Higgs oscillations. Over several Φ 0 periods, the expansion of the Universe can safely be ignored.

Equation of motions for the inhomogeneous modes
As the two coupled scalars (14) begin to oscillate, they produce particles and eventually can reheat the Universe. Since particle production in pure R 2 -inflation [32][33][34] exploits very slow gravitational dynamics, rather than the fast electroweak dynamics that cause reheating in pure Higgs inflation [21,22,26], one expects the Higgs homogeneous fields to be responsible for the early-time reheating in the mixed Higgs-R 2 case.
Particle production in the scalar sector (scalarons and Higgs bosons) is described by the following quadratic action (derived from (4) and (5)) for inhomogeneous perturbations H(t, x) andΦ(t, x), As discussed in Sec. 2, at small homogeneous fields we set cosh H 0 / √ 6M P 1, reducing the first four terms in (19) to a pair of canonical kinetic terms. The second derivatives of the scalar potential are well approximated by (see eq. (14)) In the above expressions all terms linear in fields become negative when Φ 0 gets positive value and can potentially lead to tachyonic instabilities -as we discuss in the next section.
To the leading order in background fields one can obtain from the quadratic action (19) the equations of motion for the Fourier transforms, H k and Φ k of the Higgs and scalaron inhomogeneous modes,Ḧ Using (20)- (22), these equations to leading order arë They describe two linearly coupled oscillators with time-dependent mass terms, which may lead to particle production. When Φ 0 < 0 both diagonal terms are positive, and the offdiagonal terms are not large enough to make the squared frequency negative, initiating an instability in the system. However, the mass terms rapidly oscillate with the background Higgs field, potentially producing particles. At the positive branch, Φ 0 > 0, the heavy (Higgs-like) mass eigenstate can become tachyonic. Indeed, in the potential valley with the homogeneous Higgs mode in the minimum (15), the term in the parenthesis in (25) reduces to Higgs oscillations about H min with sufficiently large amplitudes give a contribution to (27) which may turn it negative. This induces the tachyonic instability meaning instant particle production. Note that production is more efficient at small H min , where the system just left the bifurcation point and started to evolve along the potential valley. Later H min becomes large enough to prevent this term flipping sign because of the Higgs oscillations. Anyway, the larger the oscillation amplitude, the more efficient the production. The source is the homogeneous Higgs field, so the outcome is constrained not only by the Higgs amplitude (which defines the highest momentum of produced particles saturating the total density and energy of produced particles) but most severely by the amount of total energy concealed in the Higgs field. For the weak gauge bosons the mass terms are determined by the Higgs fields [28]. The quadratic Lagrangian for the W ± -bosons reads where g is the weak gauge coupling constant; in the very small field limit, H 0 → v = 246 GeV, we restore here the SM mass term of W ± -bosons. In what follows we consider the W ± -bosons only, the case of Z-bosons is similar. Defining the field-dependent variable one obtains for the Fourier modes of transverse components of W -bosons a damped Klein-Gordon equation with time-dependent mass: While the last term in the equation above, being rapidly oscillating, generically sources the particle production, in our case with functional form (28) and rather small Higgs field immediately after inflation the production is not very efficient, see Ref. [30] for details. However, the situation is different for the longitudinal modes, which also obey the Klein-Gordon equationẄ with frequency (see Ref. [24] for its conformal analog) In case of large 3-momenta k/a m T , expression (31) reduces to For the fields obeying (12) one obtains the leading contributions where the last two terms cancel along the valley, and the Higgs oscillations about its potential minimum with sufficiently large amplitude of H 0 − H min , see (16), close to a scalaron zero crossing can give a dominant tachyonic contribution. This post-inflationary feature is responsible for the violent production of longitudinal modes of the weak gauge bosons and instant preheating in the model, as we show in the next sections.

Description of particle production
In order to study the enhancement of the inhomogeneous scalar modes, they must be quantized in a canonical manner. To do so, we first remove the determinant of the metric in (19) by rescaling fields asΦ(t, x) ≡ a 3/2 Φ(t, x),H(t, t) ≡ a 3/ 2 H(t, x), what removes the Hubble friction term in (23), (24) and adds irrelevant Hubble suppressed term to the mass. To diagonalise the mass matrix we rotate (Φ,H) → (ϕ L , ϕ H ) with time-dependent mixing angle θ The resulting mass eigenvalues are Except for short moments around Φ 0 ≈ 0 one is much heavier than the other, corresponding to the oscillations in transverse directions of the potential valley and along the valley. To qualitatively understand the reheating mechanism we can focus just on the heavy inhomogeneous Higgs-like mode. To leading order, its mass equals, see eq. (22), This expression gives the estimate for the mass mode with the largest absolute value, and coincides with m 2 H for positive mass, and m 2 L for tachyonic mass. It coincides with the diagonal mass of the mode H k in (25). At Φ 0 < 0 it is positive, while can be negative for positive scalaron. The mode frequency at the minimum along the potential valley equals (27), which allows for negative values with sufficiently large amplitude of the Higgs oscillations, as explained right after eq. (27). We confirm this numerically in Sec. 5.
The diagonal modes have physical frequencies In the Hamiltonian there are terms proportional toθ, which rotate fields and momenta into one another over a timescale ∼θ −1 . This decoherence process has been studied extensively in [35]. Numerically one finds thatθ 2 |m 2 L | at all times, except during a short region of time near Φ 0 = 0, where the mode roles exchange. Outside of this region, the L and H modes are separable. Furthermore, when ω L and ω H are varying adiabatically slowly, one has WKB solutions to (25) and (26), and a particle interpretation for ϕ (L,H) . This is true for any time when the scalaron is sufficiently far away from Φ 0 = 0, so one may legally quantize ϕ (L,H) during these adiabatic time slots. The particle production can be obtained from the solutions to eqs. (23), (24) (and, in a similar way for (30)) subject to the vacuum initial conditions f k (t) = e −iωt / 2ω(k) at t → 0, where f stands for each of the ϕ H , ϕ L , W L , W T with the respective frequency. Taking these solutions for each mode at large t when the system comes back to adiabatic behaviour, one obtains the comoving number density per phase space volume d 3 k/(2π) 3 of particles of each type [35] and the physical energy density of each species is

Numerical study and results
To study the preheating process numerically we consider three particularly illustrative values of β: β ≡ (β 1 , β 2 , β 3 ) = (1.869, 2.9, 18.0) × 10 6 . The first value corresponds exactly to the system crossing the first zero in the previously-discussed bifurcation regime, where the background Higgs oscillates right about zero for an extended period, while the scalaron oscillations are nearly trapped in the vicinity of zero. Value β 2 is somewhat off the tuned value, the background fields rapidly return to one of the positive Φ 0 potential valleys with significant oscillations on present in the H 0 direction. Finally, β 3 is far away from any special points, the Higgs oscillation amplitude along the valley floor is small. 4 The evolution of the background fields near the first scalaron zero crossings is shown for these cases in Fig. 1.
Consequently, the homogeneous Higgs gets the greatest fraction of the total energy in the first scenario, and the least in the third, see Fig. 2. The corresponding masses of the  Higgs mode is less energetic, the masses of the inhomogeneous modes are not allowed to oscillate as deeply into the tachyonic regime.
To study the preheating dynamics, the mode functions for W L , and (L, H) scalar perturbations were evolved numerically according to (23), (24) and (30), subject to the vacuum initial conditions. The total energy in the inhomogeneous excitations of each type was calculated using (39). This calculation is valid only while the backreaction on the uniform background can be neglected (or, roughly, while the energy density of the produced perturbations is below the uniform background energy density).
In Fig. 4 the physical energy density is plotted for both the Higgs perturbations and W L , and compared to the background energy density 3H 2 M 2 P . In the first two plots, the background Higgs energy is drained more or less instantaneously. However, parameters further away from the tachyonic situation lead to a significantly slower rise in the corresponding energy densities (third plot).
It is insightful to investigate the spectrum for each situation. In Fig. 5 the number density of each mode is plotted immediately after the tachyonic region for β 1 and β 2 , and at a comparable point in the scalaron oscillation for β 3 . The typically produced particles are sharply peaked at low momenta, where the frequencies become tachyonic after the scalaron zero crossing. All species are non-relativistic. A warning is in order -as far as the analysis of the current paper does not account for the backreaction on the background fields, the calculation of the particle production can be trusted only while the energy is smaller than the background energy (dashed line on the plots). Even more, as far as the particle bath is produced by the oscillations in the H 0 direction, it is safe to neglect backreaction only while below the energy in H 0 , c.f. Fig. 2. However, the latter statement is not too constraining, as far as the effective tachyonic production happens only when the transfer of energy into H 0 is effective. As far as the produced gauge bosons are not relativistic, a mechanism to transfer energy into light particles (radiation) is required. It is provided by the decays of the produced gauge bosons, which are short lived, and decay within one oscillation of Φ 0 . Indeed, the time-averaged decay width [12] of longitudinal W bosons does not allow to accumulate them over more than a single scalaron oscillation. Here, we have used α W ≈ 0.025. Meanwhile, we estimate the time-averaged Higgs decay width [12] to be with Yukawa of b-quark y b ∼ 0.02, so we need not worry about Higgs decays over preheating timescales.
We therefore see that the closer one approaches the special values of β which allow the background to be significantly drained, the faster this draining occurs. What's more, this process is unambiguously driven by the tachyonic enhancement of the inhomogeneous modes.
Between any two scalaron zero crossings, the non-relativistic longitudinal weak bosons will decay into relativistic products, and the amplitude of the background Higgs mode will fall to nearly zero. Therefore, it can be expected that the system "resets" its homogeneous background to the smooth motion along the potential valley (16) upon each scalaron crossing -up to the production of long-lived inhomogeneity in the Higgs mode.
We must now ask, how many scalaron zero crossings are required to reach this finely tuned situation of efficient tachyonic reheating. In general, the homogeneous fields will not reflect exactly along the line H 0 = 0, the system will either overshoot or undershoot and end up with a significantly smaller energy in H 0 , and, therefore, in radiation. The depth of the tachyonic mass dip of the Higgs mode between the first and second scalaron zero crossing, as a function of β, is shown in Fig. 6. Also plotted is the maximum kinetic energy  in the background Higgs direction over the positive half oscillation of the scalaron -the equivalence between the two quantities is clear. Each next peak on Fig. 6 corresponds to one more oscillation of field H 0 within the negative half-oscillation of Φ 0 , c.f. Fig. 1.
Overall, we see that if the zero crossing happened near the bifurcation, energy is efficiently transferred from the background oscillations to the inhomogeneous modes. Longitudinal gauge bosons are produced more effectively (recall there are three longitudinal components for each of W + , W − and Z-bosons) and rapidly decay to light relativistic particles, leading to full preheating. Away from the bifurcation, neither the tachyonic production is effective, nor is the amount of energy in the Higgs-like oscillation is sufficient for preheating.
The narrow peaks repeat themselves with increasing β. The position of the centres of the peaks depends on the phase of the background Higgs oscillations at the moment of scalaron crossing. While it is impossible to estimate this phase at second and all later zero crossings without considering backreaction, we expect the positions of the peaks to vary chaotically, with no preference for any particular values of β. Therefore, our specially-tuned situation becomes general if one waits a few scalaron oscillations for it to be realized. This statement is more robust in the Higgs-like regime, because the tachyonic peaks are more abundant in this region of parameter space. In any case, one expects the preheating timescale in the Higgs-like limit to be much shorter than the Hubble time, and is therefore instantaneous from a cosmological perspective.

Conclusions
We have found that the dynamics of preheating in Higgs inflation regularized by the additional R 2 term are quite involved. Without the regularizing term the system had a spike like feature at zero crossings of the field in the longitudinal gauge boson mass. In the regularized version, this feature is also present, and corresponds to the motion at small values of the Higgs field and negative values of the scalaron field, however it does not lead to significant particle production. We demonstrated that immediately after the reflection from the scalaron potential and return to the positive values of the scalaron, the system may enter a strongly tachyonic regime, leading to instant preheating. Though this happens at the first zero crossing only for a model parameters in fine tuned regions, it is expected, that this happens after several zero crossings, leading to nearly immediate reheating for values of parameters, sufficiently close to the Higgs-like regime. With this conclusion we fix completely the post-inflationary evolution of the model, with reheating temperature T reh 10 15 GeV and hence the number of e-foldings N e corresponding to the scale of matter perturbations adopted as the Planck prior. Namely, from eqs. (8) we obtain (c.f. [36,37]) N e = 59, n s = 0.97, r = 0.0034.
The relevant dynamics in the system are due to the tachyonic masses appearing in the longitudinal gauge bosons (or, Goldstone bosons in the simpler model without gauge symmetry), and, to probably a lesser extent due to tachyonic mass in the Higgs boson itself. Further study with full account for backreaction on the background evolution of the fields is required for analysis at any values of the model parameters, allowing one to consistently deal with the system with equal energy distribution between the homogeneous Higgs motion and inhomogeneous modes (weak gauge bosons and Higgs particles).
The work of FB was supported in part by the STFC research grant ST/P000800/1. The part of the work devoted to production of the weak gauge bosons was supported by Russian Science Foundation grant 19-12-00393.