Higgs-like spectator field as the origin of structure

We show that the observed primordial perturbations can be entirely sourced by a light spectator scalar field with a quartic potential, akin to the Higgs boson, provided that the field is sufficiently displaced from vacuum during inflation. The framework relies on the indirect modulation of reheating, which is implemented without any direct coupling between the spectator field and the inflaton and does not require non-renormalisable interactions. The scenario gives rise to local non-Gaussianity with $f_{\rm NL}\simeq 5$ as the typical signal. As an example model where the indirect modulation mechanism is realised for the Higgs boson, we study the Standard Model extended with right-handed neutrinos. For the Standard Model running we find, however, that the scenario analysed does not seem to produce the observed perturbation.

Introduction. The cosmological perturbations observed in the Cosmic Microwave Background (CMB) are most commonly interpreted as a product of the inflaton field quantum fluctuations, which are stretched to macroscopic scales during the inflationary expansion in the early Universe.
Alternatively, the observed perturbation could be sourced by fluctuations of spectator fields which have no dynamical effect during inflation but affect the dynamics later on. Well-known examples include the curvaton model [1][2][3][4][5] and the modulated reheating scenario [6,7]. In the latter, the inflaton field decay width is modulated by spectator field fluctuations which sources the curvature perturbation and typically gives rise to significant non-Gaussianities [8,9]. Modulated reheating has been extensively studied in the past, most usually in setups which rely on the direct coupling of the inflaton field to light spectator scalars -for instance, the Higgs boson. We refer the reader to Ref. [8] for the general formalism of modulated reheating, and to Refs. [9][10][11][12][13][14][15][16][17] for models with a specific focus on the Higgs field. Further related work can be found in Refs. [18][19][20][21][22]. Spectator field fluctuations can also indirectly modulate the reheating, by introducing spatial dependence in the inflaton decay width through kinematic blocking and thereby source the primordial perturbation [11-13, 15, 16, 23]. In our previous work [23], we analysed the consequences of the Standard Model (SM) Higgs field in this framework, assuming that * alexandros.karam@kbfi.ee † tommi.markkanen@kbfi.ee ‡ luca.marzola@cern.ch § sami.t.nurmi@jyu.fi ¶ martti.raidal@cern.ch * * a.rajantie@imperial.ac.uk prior to the reheating process the field was in the vacuum state determined by the equilibrium configuration on a de Sitter background [24,25]. Although the Higgs boson does give rise to a significant scalar perturbation spectrum, such equilibrium state generically leads to power spectra characterised by a blue tilt [26,27], incompatible with the observed red-tilted CMB spectrum. In spite of that, quite naturally, spectator fields not in their vacuum state can generate the observed spectrum. Both in the modulated reheating and curvaton models, primordial perturbations with a red-tilted spectrum can be generated by spectator fields that possess non-vanishing mean values and slowly roll towards their respective vacuum states. In the curvaton setup, a crucial requirement is that after the end of inflation the spectator energy density needs to grow comparable to the dominant energy component in order for the mechanism to source significant perturbations. In the case of the Higgs boson, or of any spectator field with a quartic potential, this is problematic because the energy density of the field never dilutes slower than the dominant radiation component during the Hot Big Bang epoch [13,28]. On the other hand, in the modulated reheating scenario, such constraint can be avoided as the efficiency of the mechanism does not necessarily depend on the spectator energy density.
In this work we study a modulated reheating mechanism that allows a Higgs-like spectator field h -that is, driven by a quartic potential -to source the observed primordial perturbation spectrum despite having a subdominant energy density throughout the cosmic evolution. The setup is similar to that of Ref. [26], but in this work we concentrate on the limit of large spectator field values which allows to generate the observed red spectral tilt. We discuss how the mechanism could be implemented in extensions of the SM, with the purpose of identifying the spectator field with the Higgs boson. The main conclusion of our work is that a set-up with a quartic spectator with no coupling to the inflaton and without non-renormalizable operators can solely generate a curvature perturbation that is consistent with current observational bounds. The model. Following Ref. [23], we consider a Lagrangian given by where h is the spectator scalar field, the scalar singlet φ is the inflaton and Ψ is a Dirac fermion. The form of the inflaton potential considered in Eq. (1) is not a prerequisite of our setup and is only used to provide a simple and concrete template. It is crucial to our purposes that the decay rate of the inflaton is subject to the indirect modulation from a spectator field, such that the decay channel is blocked by kinematics when the spectator field has the value Since the Lagrangian also allows for decay of the spectator field, for simplicity we restrict here to parameters such that the effective mass of the spectator field, √ 3λh, is much smaller than that of the fermion, y h h, and thus the process can be safely ignored. It is worth pointing out that if we identify h with the SM Higgs boson, the remaining couplings to SM fields, such as gauge bosons, can be safely neglected as they do not lead to any significant depletion of the Higgs field value over time scales relevant for the modulation effect, corresponding to O(1) Higgs oscillations after the inflationary expansion [29]. Note also that while the tree-level action (1) does not directly couple the inflaton to the Higgs boson, such a coupling is generated by fermions at the loop level. However, because of the small values of the Yukawa coupling y h that we are considering in this work, we neglect these radiative effects.
Denoting by h * the field value at the horizon crossing of a mode k * = a * H * during inflation, we now focus on the case h * > h kin , indicating that the inflaton decay is kinematically blocked by the spectator field via the indirect modulation mechanism. We assume the h-field to be light if compared to the inflation scale, hence it fluctuates locally around its VEV and the kinematic blocking of the inflaton decay is thus lifted at different times in different locations in the Universe. This implements the modulated reheating mechanism [7] in our scheme and allows for the production of significant curvature perturbations [8,13,23]. Specifically, we require that the spectator field effective mass satisfies V (h * )/(3H 2 * ) < 0.01, which for our model translates to λ h h 2 * < 0.01H 2 * . For h * < M P , the condition implies Ω h * 10 −3 , so the spectator field energy density is necessarily subdominant with respect to the inflaton contribution.
As for the spectrum, at the equilibrium, spectator field fluctuations in a de Sitter background result only in a blue tilted spectra [26,27] (n s > 1) and therefore cannot source the observed primordial perturbation. The situation is different for field configurations far from equilibrium, i.e., in the mean field limit h * > h 2 eq , when fluctuations around the mean field value h * yield either red or blue tilts. In the setup specified by Eq. (1), the out-of-equilibrium configuration h * > h 2 eq corresponds to an atypical field configuration, or initial condition, which we suppose to be accidentally realised in our observable patch. This is similar to what is often assumed in the context of the curvaton scenario, for example see [1]. Alternatively, one might think of modifying the spectator field potential such that the initial condition could be dynamically realised. Further study of this question is however beyond the scope of this work. Here we simply assume the mean field limit h * > h 2 eq due to the phenomenological reason that it is required to get the observed spectral tilt. Using the de Sitter equilibrium result h 2 eq 0.36H * /λ 1/4 h [25], the validity of the mean field limit requires h * > 0.36H * /λ 1/4 h . Combining this with the mass bound λ h h 2 * < 0.01H 2 * , our constraints are specified by which further implies the bound λ h < 0.006. Analytical estimates. Simple analytical estimates of the spectrum and non-Gaussianity of curvature perturbations produced by the indirect modulation mechanism can be obtained in the limit where the inflaton decays rapidly after the kinematic blocking is lifted. This occurs when h(t kin ) = h kin for the first time in the evolution of the spectator field value. In the discussion, we neglect the slow evolution of h during the inflationary expansion and approximate its value at the end of this process by h end h * . The spectator field starts to oscillate when H osc ≈ √ 3λ h h * , and the field behaviour prior to the first zero crossing can be approximated quite accurately by [23] h where H ∝ a −3/2 for a Universe dominated by inflaton oscillations in a quadratic potential.
In order to ensure the prompt decay of the inflaton field during the first available window allowed by kinematics, we impose that ρ φ /ρ tot falls below 10 −5 as |h| < h kin during the first half of the first h oscillation -i.e. the first period when h moves from h kin through zero to −h kin . For ρ φ /ρ tot < 10 −5 , the inflaton contribution to the number of e-folds from later times, O(ρ φ /ρ tot ) Hdt, can be neglected when studying perturbations of order ∆N ∼ 10 −5 . As a first approximation, we can then approximate that the Universe changes from matter to radiation domination at t = t kin , and from Eq. (5) it follows that By using ρ φ ∝ exp(− Γdt), and estimating the integral over the first window allowed by kinematics as Γdt ∼ 2Γ 0 t k t0 dt, where t 0 and t k are determined by solving Eq. (5) for h(t k ) = h kin and h(t 0 ) = 0, the fast decay condition − Γdt > ln10 −5 translates into the bound Here we have also assumed the limit h kin h * . The numerical analysis presented below indeed confirms that this inequality well determines the regime where our analytical approximations are fully applicable.
The curvature perturbation in Fourier space can be computed in terms of the δN expression [30] as where N is the number of e-folds from the horizon crossing k * = a * H * to a final constant energy surface after the inflaton decay, a prime denotes differentiation with respect to the field value h * at horizon crossing, and δh n * are convolutions. Throughout this work we concentrate on the limit where the inflaton contribution to the curvature is negligible, ζ φ ≈ (1/ √ 2 * )H * /(2πM P ) 10 −5 , and therefore omit it in all expressions.
When the condition in Eq. (7) is satisfied, we can well approximate the evolution of the Universe using a sudden transition from matter to radiation domination at t kin , which yields [23] N ≈ −(1/6)H kin /H kin , with H kin given by Eq. (6). Substituting this into Eq. (8), the power spectrum of curvature perturbation becomes where in the last step we expanded Eq. (6) in the limit h * h kin . The spectral index is given by [31] where ≡ −Ḣ/H 2 is computed at the horizon crossing. Comparing our result with the observations, P ζ (k * ) = (2.100 ± 0.030) × 10 −9 and n s (k * ) = 0.965 ± 0.004 at the pivot scale k * = 0.05 Mpc −1 [32], we find that the analytical estimates agree with the measurements for h * /H * 580 and * 0.018 (assuming λ h 10 −8 so that the positive contribution from the second term in Eq. (10) is negligible). For h * 580H * , the assumed mean field condition (4) implies λ h 1.5×10 −13 . Note also that the derivation of the estimate (9) assumes the fast decay condition in Eq. (7), which constrains from above the range of λ h values for which the analytical estimate can be used. For instance, taking y φ = 1, y h = 10 −3 , m φ /H * = 0.1 and setting h * = 580H * , Eq. (7) yields λ h 1.7 × 10 −12 and, for this particular example, cases with larger values of λ h need to be studied numerically -see the results further below.
The second order term in Eq. (8) contributes to the local bispectrum with a non-Gaussianity amplitude f NL = (5/6)N /(N 2 ) given by where, in the last step, we have again expanded Eq. (6) in the limit h * h kin . The prediction f NL ∼ 5 is a characteristic signal of our mechanism in the limit of fast decay (i.e. when Eq. (7) holds), and when h * h kin so that O(h kin /h * ) terms in Eq. (6) can be neglected. This analytical result is confirmed by the numerical analysis presented below. We remark that this level of non-Gaussianity is still compatible with the current bound f local NL = −0.9 ± 5.1 [33] at the 1σ confidence level. Importantly, the next generation surveys will reach the sensitivity required to either detect the signal or rule out the mechanism as the source of the primordial perturbation in the aforementioned limits. In particular, upcoming large-scale structure probes such as DESI and Euclid are expected to constrain the primordial non-Gaussianity with uncertainties of O(1) [34][35][36]. Numerical results. The analytical approximations we proposed apply to a regime where the inflaton decays efficiently in the kinematically allowed window prior to the first zero crossing of the spectator field value. To go beyond this limit, we numerically solve the following system of equationṡ h + 3Hḣ + λ h h 3 = 0, where Γ denotes the real part of Eq. (2). The initial conditions are determined by matching the solution of the system composed by Eq. (12) with the first line replaced byφ + 3Hφ + V (φ) = 0, at a matching time chosen well after the end of inflation (so that w φ = 0) and well before t kin , so that Γ = 0. The initial conditions for this second system of equations are set by the inflationary slow roll solution, together withḣ * = 0, ρ r * = 0 and the initial spectator field value h * . Note that modelling the inflaton decay by inserting the diffusion term Γφ in its equation of motion would not be justified asφ/Γ φ when the decay takes place. We have checked that our results do not depend on the choice of the matching time within the given boundaries.
In our computation we set the inflaton parameters in Eq. (1) to m φ = 1.55 × 10 9 GeV, λ φ = 3.75 × 10 −21 and y φ = 1, which yields H * ≈ 1.6 × 10 10 GeV and * ≈ 0.017 . . . 0.018 at the horizon crossing of the pivot scale k * = 0.05 Mpc −1 . The precise values vary depending on the reheating process controlled by the spectator field couplings. For this choice of inflaton parameters, the inflaton contribution to the curvature perturbation is negligible, ζ 2 φ ∼ 10 −17 , and will therefore be omitted henceforth. Note also that our results are not limited to the precise form of the inflaton potential used in Eq. (1). Any other potential leading to same H * and * , and reducing to a quadratic form at the end of inflation, would give essentially the same results for the curvature perturbation induced by the modulation mechanism.
To compute the curvature perturbation, we construct a grid of initial values h * (the results shown in the figures are computed using a grid of 100 points and a step size ∆χ * = 0.78H * ). For each grid point we numerically evolve the system Eq. (12) up to a common final ρ chosen such that ρ φ /ρ < 10 −10 and determine the corresponding number of e-folds N (h * ). From this data we numerically evaluate the first and second derivatives of N (h * ) with respect to the initial field value h * , corresponding to the coefficients of the δh * (k * ) powers in Eq. (8). For the spectator field perturbations at the horizon crossing δh * (k * ) in Eq. (8), we use the linear perturbation theory result for the two point function amplitude, P δh * (k * ) = (H * /(2π)) 2 , and neglect any non-Gaussianity in δh * (k * ). The spectrum of the curvature perturbation and the non-Gaussianity parameter f NL are then given by and the spectral index by Eq. (10). All these quantities are readily evaluated from the numerical data. The results for P ζ , n s and f NL at the pivot scale k * = 0.05 Mpc −1 are shown in Figs. 1, 2 and 3, respectively. We have checked that decreasing the final value of ρ does not affect the results, confirming thereby that the curvature perturbation has relaxed to a constant well within  [32]. The circles mark the smallest field value which satisfies the mean-field condition assumed in our analysis.
and, correspondingly, the upper limit of the fast decay region decreases with larger values of λ h . As expected, in the fast decay regime the numerical results agree relatively well with the analytical approximations in Eqs. (9), (10) and (11). For the cases shown in Figs. 1 and 3, the analytical estimates for P ζ and f NL (obtained by using the full form of Eq. (6) in Eqs. (9) and (11)) respectively deviate less than 10% and 15% from the corresponding numerical results in the fast decay region h * /H * 0.6λ Increasing h * further delays the inflaton decay successively to the second, third, or following zero crossing of the spectator field value, leading to a drastic amplification of the non-Gaussianity and increasingly complicated non-monotonous forms for both P ζ (h * ) and f NL (h * ), as seen in Figs. 1 and 3, respectively. This appears to be a generic feature of the setup, indicating that configurations not in the fast decay regime (7) tend to generate unacceptably large non-Gaussianity, possibly apart from tuned configurations around points where f NL accidentally crosses zero.
As for the Yukawa coupling y h , larger values of this quantity decrease h kin in Eq. (3). This narrows the window where the inflaton decay is kinematically allowed around the spectator zero crossings and makes the decay process less efficient. Correspondingly, the region where the spectator field value crosses zero more than once before the decay, i.e. the region where the fast decay condition (7) does not hold, is pushed towards smaller values of h * as y h grows. The dependence of P ζ and f NL as functions of y h is illustrated in Figs. 4 and 5, respectively, for a sample of y h values. Implications for inflationary model building. In our analysis, the inflaton potential affects directly only the value of the spectral index Eq. (10). This is not specific to our setup and applies also to generic modulated reheating and curvaton models when the curvature perturbation is dominated by the spectator field. We remark that a generic prediction of such setups is a negligible tensor-to-scalar ratio r, which follows from the low inflationary scale required to keep the perturbations sourced by the inflaton subdominant. In particular, this means that potentials leading to too large tensor perturbations in the inflaton dominated limit can use the mechanism proposed here, as well as other modulated reheating or curvaton models, to achieve agreement with observations. An extensive analysis of inflaton potentials compatible with the curvaton framework was presented in Encyclopaedia Curvatonis [37] and qualitatively similar conclusions apply to our setup.
Embedding into the SM. The above analyses relied on the interplay between the inflaton, a generic spectator field characterised by a quartic potential and a fermion that interacts with both the scalar fields. By regarding the parameters in the Lagrangian (1) as free and working at the level of classical potential, the presented analytic and numerical computations have highlighted the limits where the framework can match the CMB observations. In Ref. [23] we have shown that the required modulation mechanism is already implemented in a well-known extension of the SM, where the particle content considers a singlet inflaton and right-handed Majorana neutrinos. It is therefore of interest to investigate whether the same setup allows for the generation of the CMB through the dynamics discussed in the present work.
If the spectator field is to be identified with the Higgs boson, our analysis applies provided that the effective Higgs potential can be approximated by the tree level form: V (h) = λ h (µ)h 4 /4. The renormalization scale is set to µ ∼ √ y t h, with y t being the top-quark Yukawa coupling, corresponding to the largest effective mass scale when h * H * . The effective potential has a maximum where λ h (h max ) = −β h (h max )/4 and the presented analysis, which used V = λ h h 4 /4 with a constant coupling, is applicable only for |h| |h max |. Using SM two-loop beta functions in the R ξ gauge and MS scheme, for a (world average) top-quark mass of m t = 172.9 GeV we find λ h (h max ) O(10 −5 ). However, to produce the observed spectral tilt in the quartic setup, we require λ h (h * ) < 10 −8 -see the discussion below Eq. (10). Moreover, even smaller values in the ballpark of λ h (h * ) 10 −11 are necessary to avoid large non-Gaussianities, as shown in Fig. 3. The specific points where f NL (h * ) in Fig. 3 accidentally crosses zero might represent a possible caveat. However, although we have not studied this in detail, it is expected that even such configurations lead to too large non-Gaussianity in the form of the trispectrum amplitude g NL ∝ N /N 3 . We thus conclude that, assuming the SM running of couplings, the Higgs boson seemingly fails to source the observed primordial perturbation through the mechanism studied here.
One may of course ask if it is possible to circumvent these problems in extended setups with new physics coupled to the Higgs boson. Here we do not address this important question in detail, but we briefly comment on possibly viable phenomenological structures. First, a negative threshold correction λ h (h * ) → λ h (h * ) + δλ h 10 −10 could make the scenario work, at the cost of finetuning δλ h and h * close to the threshold scale. Also new scalar degrees coupled to the Higgs would act in favour of the scenario as they contribute positively to β h and make it possible to decrease λ h (h max ) = −β h (h max )/4. An interacting fixed point λ h (µ) 10 −11 , β h (µ) = 0 at, or below, the scale h * would also work, and it would be interesting to investigate if such fixed point could arise, for instance, from anomalously broken scale invariance. More generally, the tree-level Higgs potential could differ from the quartic form, for example due to the nonminimal curvature coupling ξRh 2 or couplings to other fields with non-vanishing VEVs. This would change the reheating dynamics compared to our analysis, leading to potentially different conclusions. One could also think of scenarios where h * is located beyond the vacuum present at large field values, contrary to what we have assumed. Suppose, for instance, that the corresponding minimum h min is brought to sub-Planckian values by the interaction with an additional scalar field used to generate the right-handed neutrino Majorana masses via symmetry breaking. One could then arrange h * > h kin > h min , so that the reheating would commence at h ∼ h kin and, if thermal corrections rapidly lift the minimum h min , the Higgs could relax into the electroweak vacuum. A more careful assessment of these possibilities, however, requires dedicated analyses that go beyond the scope of this work. Conclusions. We have shown that a spectator field with a quartic potential can alone source the observed primordial perturbation through the modulated reheating mechanism, realised without direct couplings between the spectator and the inflaton field and by using only renormalisable operators. To obtain a red spectrum , the spectator field needs to be sufficiently displaced from vacuum during inflation, its quartic coupling λ h 10 −8 and, as in generic spectator models, the inflationary dynamics needs to yield a suitable value for the slow roll parameter H at the horizon crossing of observable modes. The setup gives rise to primordial non-Gaussianities of the local type. We find that when the inflaton decay completes sufficiently fast (during the first spectator oscillation that lifts the kinematical blocking associated to the modulation channel), the bispectrum amplitude is set by f NL ∼ 5, well in agreement with the current bounds and a testable signature by upcoming surveys. For λ h 10 −11 , this is the case for a relatively broad range of initial spectator field values h * . For larger values of λ h , the inflaton decay is slower which gives rise to a growing tension between obtaining the observed amplitude of perturbations and maintaining the non-Gaussianity within the observational bounds.
Because the primordial perturbations are not directly related to the inflaton potential, the proposed framework, as well as modulated reheating and curvaton models in general, allows a wide range of inflationary models to come in agreement with present data, including those built on quadratic and quartic potentials. The predicted scalar-to-tensor ratio is generally negligible.
The scenario we discussed allows, in principle, for a straightforward identification of the spectator field with the Higgs boson in a popular extension of the Standard Model with Majorana right-handed neutrinos. Although the required indirect modulation mechanism can be effectively implemented, matching the observed perturbation needs a modification of the Standard Model renormalisation group equations. In fact, using the Standard Model running and approximating the Higgs effective potential with the tree-level term, we find that the required value of the quartic coupling, λ h 10 −11 , cannot be obtained in the regime where a quartic form well approximates the full potential. Whereas more precise studies that use the full form of the Higgs effective potential are needed in order to fully assess the possibility, we have briefly discussed promising phenomenological extensions of the Standard Model which could allow the Higgs boson to source the observed perturbations in the considered scenario based on a quartic potential. It would be of interest to study how these proposals could be realised in concrete particle physics setups.