Exploring Inflation Options for Warm Dark Matter Coupled to the Higgs Boson ()

Bruce Hoeneisen^{}

Universidad San Francisco de Quito, Quito, Ecuador.

**DOI: **10.4236/ijaa.2023.133013
PDF
HTML XML
108
Downloads
480
Views
Citations

Universidad San Francisco de Quito, Quito, Ecuador.

We extend the Standard Model with a scalar warm dark matter field *S* with an interaction with the Higgs boson *∅*. This warm dark matter scenario is in agreement with cosmological observations if *S* and *∅* come into thermal and diffusive equilibrium before the temperature drops below the Higgs boson mass *m*_{H}. We study inflation driven by the fields *∅* or *S*, and also study preheating and reheating, in order to constrain the parameters of this extension of the Standard Model. It is remarkable that, with the current data, these models pass a closure test with no free parameters.

Keywords

Share and Cite:

Hoeneisen, B. (2023) Exploring Inflation Options for Warm Dark Matter Coupled to the Higgs Boson. *International Journal of Astronomy and Astrophysics*, **13**, 217-235. doi: 10.4236/ijaa.2023.133013.

1. Introduction

The present study is motivated by two recent observations. Firstly, the precise measurements of the Higgs boson mass *m _{H}* at the Large Hadron Collider (LHC), and of the top quark mass

842 ± 13 m/s, and 1471 ± 19 m/s, for spin 0, $\frac{1}{2}$ , and 1 dark matter particles, respectively [4] . Note that spin $\frac{1}{2}$ and spin 1 particles $S$ are disfavored by the

measurements.
$S$ has a predicted mass
${M}_{S}=150\pm 2\text{\hspace{0.17em}}\text{eV}$ [4] . If
$S$ is real, then
$S$ is the dark matter particle.^{1} If
$S$ is complex, it may decay to two massive gauge bosons *V* which are the dark matter in this case [10] . Additional observables studied, that obtain consistent results, are the redshift of first galaxies [11] , the ultra-violet luminosity of first galaxies [6] , and the related “reionization optical depth” [6] . For a summary of measurements see [6] .

For these reasons, in the present study we explore inflation options of the Standard Model extended with a scalar warm dark matter field $S$ with a contact coupling $-\frac{1}{2}{\lambda}_{hS}\left({\varphi}^{\u2020}\varphi \right){S}^{2}$ to the Higgs boson $\varphi $ .

In Section 2 we take a journey towards the past, based on current observations and on the assumption of the validity of the Standard Model up to high energies. In Sections 3 to 7 we take a journey towards the future, starting from an assumed early stage of inflation of the universe driven by either
$S$ or
$\varphi $ , and the assumption that quantum fluctuations during inflation are the seeds of galaxies and of the observed large-scale power spectrum of density perturbations. The meeting, or not, of the two journeys provides a consistency test of the preceding assumptions. It is quite amazing, and non-trivial, that current data and the listed assumptions do pass the consistency test *with no adjustable parameters*! A summary of results is presented in the concluding section.

2. The Expansion History of the Universe

A summary of the expansion history of the universe is presented in Figure 1. The horizontal axis
$x\equiv {M}_{P}r/a$ is the comoving distance in Planck units, and the vertical axis
$y=a$ is the expansion parameter
$a$ , both with natural logarithmic scales.
${M}_{P}\equiv \sqrt{1/{G}_{N}}\approx 1.2\times {10}^{19}\text{\hspace{0.17em}}\text{GeV}$ is the Planck mass (we adopt the definition of [1] ). The two vertical lines correspond to two particles fixed in space: the vertical axis is the observer, *i.e.* us, and the red line is a point fixed in space at the reference comoving distance
$1/k=\text{Mpc}/0.05=\mathrm{exp}\left(134.9\right)/{M}_{P}$ . The comoving distance to the horizon relative to us, indicated by the blue line, is
${M}_{P}/\left(aH\right)$ , where the Hubble parameter *H* is the relative expansion rate. The

Figure 1. Past and future history of the universe. Vertical axis: natural logarithm of the expansion parameter
$a$ , *i.e.* “e-folds”. Horizontal axis: natural logarithm of comoving distances
${M}_{P}r/a$ in Planck units. The red line corresponds to a fixed point in space at a reference comoving distance
$1/k=\text{Mpc}/0.05=\mathrm{exp}\left(134.9\right)/{M}_{P}$ . The blue line is the comoving distance to the horizon
${M}_{P}/\left(aH\right)$ in Planck units. The thin green box (with a width too small to be resolved on this scale) is the *prediction*
$\left({M}_{P}/\left({a}_{f}{H}_{f}\right)\mathrm{,}{a}_{f}\right)$ for the end of inflation with the observed values of
${P}_{\zeta}\left(k\right)$ and
${n}_{s}$ , at 68% confidence (from Table 2), *with no free parameter*! The Cosmic Microwave Background radiation (CMB) that we see today propagates to us from point “X”. Symbols “−” and “×” illustrate segments of “light” rays (or should we say “dark” rays?) inside and outside of the horizon (see text).

expansion parameter $a\left(t\right)$ , in a spatially flat and homogeneous universe, evolves according to the Friedmann equation [1] :

${H}^{2}\equiv {\left(\frac{1}{a}\frac{\text{d}a}{\text{d}t}\right)}^{2}=\frac{8\pi}{3{M}_{P}^{2}}\rho \equiv {H}_{0}^{2}\frac{\rho}{{\rho}_{\text{crit}}}\mathrm{.}$ (1)

where
${\rho}_{\text{crit}}\equiv 3{H}_{0}^{2}/\left(8\pi {G}_{N}\right)$ is the critical density. Today,
$H\left({t}_{0}\right)\equiv {H}_{0}$ . The density has contributions from radiation, and matter, and we have included the cosmological constant Λ (or vaccuum energy) contribution into *ρ*:
${\rho}_{\text{crit}}{\Omega}_{\Lambda}={\rho}_{\text{crit}}\Lambda /\left(3{H}_{0}^{2}\right)$ :

$\rho \equiv {\rho}_{\text{crit}}\left(\frac{{\Omega}_{r}}{{a}^{4}}+\frac{{\Omega}_{m}}{{a}^{3}}+{\Omega}_{\Lambda}\right)\mathrm{.}$ (2)

It is handy to have the solutions to (1) when one term in *ρ* dominates, see Table 1.

Table 1. *t* or
$a$ dependence of some cosmological quantities when only one type of density dominates.
$1/H$ is the proper distance to the horizon.
$1/\left(aH\right)$ is the comoving distance to the horizon. The proper distance to a “galaxy”, *i.e.* a fixed point in space, is
$\propto a$ . The comoving distance to a “galaxy” is a constant. See Figure 1.

The blue straight lines in Figure 1 can be understood by noting that the comoving distance to the horizon 1/(
$aH$ ) is proportional to
${a}^{-\text{1}}$ ,
$a$ ,
${a}^{1/2}$ , and
${a}^{-\text{1}}$ , for *ρ* dominated, respectively, by vacuum energy during inflation until
$a={a}_{f}$ , radiation until
${a}_{\text{eq}}$ , matter until
${a}_{\Lambda}$ , and the cosmological constant (extended to the future to be seen on the scale of the figure). (There are intervals of
$a$ in the range
${a}_{f}<a<{a}_{r}$ with
$\rho \propto {a}^{-3}$ if the inflaton field has a dominating quadratic potential, but this will not be the case in the present scenario.)

Note that a point in space at the reference comoving distance
$1/k$ exits the horizon at
${a}_{i}$ , and re-enters the horizon at
${a}_{r}$ , *i.e.*

$k={a}_{i}{H}_{i}={a}_{r}{H}_{r}.$ (3)

Radiation decouples from matter at
${a}_{\text{dec}}$ as indicated by the black horizontal line in Figure 1. The Cosmic Microwave Background Radiation that we observe today propagates freely since *t*_{dec} from the “CMB” point marked with an X in Figure 1.

The slope of “light” rays is

$\frac{\text{d}y}{\text{d}x}=\frac{1}{\pm {\left(rH\right)}^{-1}\pm 1}\mathrm{,}$ (4)

so $\text{d}y/\text{d}x\approx 0$ if $r$ is inside the horizon, while $\text{d}y/\text{d}x\approx \pm 1$ if $r$ is outside the horizon. Segments of light rays are shown with symbols “−” and “×” in Figure 1.

The Einstein equation of General Relativity obtains (1), and also the energy conservation equation

$\frac{\text{d}\rho}{\text{d}t}+3H\left(\rho +p\right)=\mathrm{0,}$ (5)

where *p* is the pressure:

$p={\rho}_{\text{crit}}\left(\frac{{\Omega}_{r}}{3{a}^{4}}+0-{\Omega}_{\Lambda}\right)\mathrm{.}$ (6)

Let us now obtain the temperature in the range
${m}_{t}<T\left(a\right)<{T}_{\text{reh}}$ , where *T _{reh}* is the reheating temperature when the Standard Model particles and
$S$ obtain thermal equilibrium. From conservation of entropy

${g}_{\text{reh}}{a}^{3}{T}^{3}={g}_{\gamma}{a}_{0}^{3}{T}_{0}^{3}+{g}_{\nu}{a}_{0}^{3}{T}_{\nu 0}^{3}+{g}_{S}{a}_{0}^{3}{T}_{S0}^{3}\mathrm{,}$ (7)

where ${T}_{0}=2.7255\text{\hspace{0.17em}}\text{K}$ is the present photon temperature [1] , and $g\equiv {N}_{b}+\frac{7}{8}{N}_{f}$

for each particle species, with *N _{b}* or

$\frac{a}{{a}_{0}}=0.332\frac{{T}_{0}}{T}\mathrm{.}$ (8)

The corresponding density is

$\rho =\frac{{\pi}^{2}}{30}{g}_{\text{reh}}{T}^{4}\mathrm{.}$ (9)

3. Inflation with a Quadratic Potential

In this Section we take $S$ to be the field that drives inflation. Our metric for flat space is

$\text{d}{s}^{2}={g}_{\mu \nu}\text{d}{x}^{\mu}\text{d}{x}^{\nu}=\text{d}{t}^{2}-a{\left(t\right)}^{2}\left(\text{d}{x}^{2}+\text{d}{y}^{2}+\text{d}{z}^{2}\right)\mathrm{.}$ (10)

Greek indices go from 0 to 3, while latin indices go from 1 to 3. The Lagrangian, wave equation, energy-momentum tensor, energy density, and pressure of a real scalar field $S$ , are [1]

$\mathcal{L}=\frac{1}{2}{\partial}_{\mu}S\cdot {\partial}^{\mu}S-V\left(S\right)\mathrm{,}$ (11)

$0=\frac{\partial \mathcal{L}}{\partial S}-{\partial}_{\mu}\frac{\partial \mathcal{L}}{\partial \left({\partial}_{\mu}S\right)}=-\frac{\text{d}V\left(S\right)}{\text{d}S}-{\partial}_{\mu}{\partial}^{\mu}S\mathrm{,}$ (12)

${T}^{\mu}{}_{\nu}=\frac{\partial \mathcal{L}}{\partial \left({\partial}_{\mu}S\right)}{\partial}_{\nu}S-{g}^{\mu}{}_{\nu}\mathcal{L}={\partial}^{\mu}S\cdot {\partial}_{\nu}S-{g}^{\mu}{}_{\nu}\mathcal{L}$ (13)

${T}^{0}{}_{0}=\rho =\frac{1}{2}{\stackrel{\dot{}}{S}}^{2}+V\left(S\right)\mathrm{,}$ (14)

$-{T}^{i}{}_{j}=p{\delta}_{j}^{i}=\left[\frac{1}{2}{\stackrel{\dot{}}{S}}^{2}-V\left(S\right)\right]{\delta}_{j}^{i}\mathrm{.}$ (15)

We have neglected spatial derivatives in these last two equations. *ρ* and *p* are the density and pressure of the homogeneous fluid in its rest frame. From (5), (10) and (14) we obtain

$\stackrel{\xa8}{S}+\left(3H+\Gamma \right)\stackrel{\dot{}}{S}=-\frac{\text{d}V}{\text{d}S}\mathrm{.}$ (16)

Dots are derivatives with respect to *t*. In this equation we have included the total decay rate Γ of *S* to discuss reheating in Sections 6 and 7. It is convenient to define the “slow-roll” parameters:

$\epsilon \equiv \frac{{M}_{P}^{2}}{16\pi}{\left(\frac{\text{d}V/\text{d}S}{V}\right)}^{2}\text{\hspace{1em}}\text{and}\text{\hspace{1em}}\eta \equiv \frac{{M}_{P}^{2}}{8\pi}\frac{{\text{d}}^{2}V/\text{d}{S}^{2}}{V}\mathrm{.}$ (17)

Initial conditions for
$S\left({t}_{1}\right)\equiv {S}_{1}$ and
$a\left({t}_{1}\right)\equiv {a}_{1}$ are such that
$V\left({S}_{1}\right)\gg {\stackrel{\dot{}}{S}}^{2}/2$ at
${t}_{1}$ . While
$\stackrel{\xa8}{S}$ can be neglected in (16) we are in the “slow-roll” regime. Necessary conditions for the slow-roll approximation to be valid are
$\epsilon \ll 1$ and
$\eta \ll 1$ . In this approximation
$\epsilon \approx -\stackrel{\dot{}}{H}/{H}^{2}$ . Note that inflation ends, *i.e.*
$\stackrel{\xa8}{a}=0$ , when
$-\stackrel{\dot{}}{H}/{H}^{2}=1$ .

In the slow-roll approximation, *H*, *ρ*, and *p* are approximately constant, and
$a\left(t\right)\propto {\text{e}}^{Ht}$ grows approximately exponentially with time. In the slow-roll approximation, the number of “e-folds” of growth of
$a\left(t\right)$ until
$\epsilon =1$ , is

$N\equiv \mathrm{ln}\left(\frac{{a}_{f}}{{a}_{i}}\right)={\displaystyle {\int}_{{t}_{i}}^{{t}_{f}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}H\text{d}t=-\frac{2\sqrt{\pi}}{{M}_{P}}{\displaystyle {\int}_{{S}_{i}}^{{S}_{f}}}\frac{\text{d}S}{\sqrt{\epsilon}}\mathrm{.}$ (18)

The dimensionless power spectrum^{2} of density fluctuations at horizon re-entry is [1] [13] [14]

${P}_{\zeta}\left(k\right)=\frac{{H}_{i}^{2}}{\pi {\epsilon}_{i}{M}_{P}^{2}}\approx \frac{1}{24{\pi}^{2}}\frac{{V}_{i}}{{\epsilon}_{i}}{\left(\frac{8\pi}{{M}_{P}^{2}}\right)}^{2}\mathrm{,}$ (19)

and the tensor gravitational wave power spectrum at horizon re-entry is

${P}_{t}\left(k\right)={\left(\frac{{H}_{i}}{2\pi}\right)}^{2}\frac{64\pi}{{M}_{P}^{2}}\approx \frac{2}{3{\pi}^{2}}{V}_{i}{\left(\frac{8\pi}{{M}_{P}^{2}}\right)}^{2}\mathrm{.}$ (20)

Note that ${H}_{i}$ and ${\epsilon}_{i}$ correspond to ${a}_{i}$ . The power spectrum ${P}_{\zeta}\left(k\right)$ at the reference wavenumber ${k}_{\text{ref}}=0.05\text{\hspace{0.17em}}{\text{Mpc}}^{-1}$ is measured to be $\mathrm{exp}\left(3.044\pm 0.014\right)\times {10}^{-10}\approx 2.1\times {10}^{-9}$ [1] . The tensor-to-scalar ratio for single-field slow-roll inflation is [1]

$r\equiv \frac{{P}_{t}}{{P}_{\zeta}}=16{\epsilon}_{i}$ (21)

A “pedestrian” derivation of (19) is presented in Appendix A.

Let us now consider the quadratic potential

$V\left(S\right)=\frac{1}{2}{\stackrel{\xaf}{m}}_{S}^{2}{S}^{2}\mathrm{.}$ (22)

The slow-roll parameters are

$\epsilon =\eta =\frac{1}{4\pi}{\left(\frac{{M}_{P}}{S}\right)}^{2}\mathrm{.}$ (23)

The solution of the evolution equations in the slow-roll approximation are

$S={S}_{i}-\frac{{\stackrel{\xaf}{m}}_{S}{M}_{P}}{\sqrt{12\pi}}\left(t-{t}_{i}\right)\mathrm{,}$ (24)

$a\left(t\right)={a}_{i}\mathrm{exp}\left(\sqrt{\frac{4\pi}{3}}\frac{{\stackrel{\xaf}{m}}_{S}}{{M}_{P}}{\displaystyle {\int}_{{t}_{i}}^{t}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}S\text{d}t\right)\mathrm{.}$ (25)

End of inflation, *i.e.*
$\epsilon =1$ , corresponds to
${S}_{f}={M}_{P}/\sqrt{4\pi}$ . The number of e-folds of inflation is:

$N\equiv \mathrm{ln}\left(\frac{{a}_{f}}{{a}_{i}}\right)=2\pi \frac{{S}_{i}^{2}-{S}_{f}^{2}}{{M}_{P}^{2}}.$ (26)

The power spectrum of density fluctuations, and the spectral tilt, are

${P}_{\zeta}=\frac{16\pi}{3}{\left(\frac{{\stackrel{\xaf}{m}}_{S}}{{M}_{P}}\right)}^{2}{\left(\frac{{S}_{i}}{{M}_{P}}\right)}^{4},$ (27)

${n}_{s}-1\equiv \frac{\text{d}\mathrm{ln}{P}_{\zeta}}{\text{d}\mathrm{ln}k}\approx -6{\epsilon}_{i}+2{\eta}_{i}\approx -\frac{2}{N}.$ (28)

From the measured
${n}_{s}\approx 0.965$ , we obtain
$N\approx 57.1$ and
${S}_{i}\approx 3.0{M}_{P}$ , and from the measured
${P}_{\zeta}\approx 2.1\times {10}^{-9}$ , we obtain
${\stackrel{\xaf}{m}}_{S}\approx 1.2\times {10}^{-6}{M}_{P}$ . In conclusion, a quadratic potential is appropriate for cold dark matter, not for the warm dark matter studied in this article. Also, preheating with a quadratic potential requires
${M}_{S}\gg {m}_{H}$ , *i.e.* cold dark matter. We therefore turn to a study of inflation with a quartic potential.

4. Inflation with a Quartic Potential

To the Standard Model we add a warm dark matter real scalar field
$S$ , with* Z*_{2} symmetry
$S\leftrightarrow -S$ , with a contact coupling
${\lambda}_{hS}$ to the Higgs boson field
$\varphi $ . (An alternative is complex
$S$ with a Lagrangian with *U*(1) symmetry that is, or is not, broken by the
$S$ and
$\varphi $ ground state [10] . Complex
$S$ only requires minimal changes of notation, e.g.
${S}^{2}\to {S}^{\mathrm{*}}S$ .) The action is

${\mathcal{S}}_{J}={\displaystyle \int {\text{d}}^{4}x}\sqrt{-g}\left\{{\mathcal{L}}_{\text{grav}}+{\mathcal{L}}_{S}+{\mathcal{L}}_{S\varphi}+{\mathcal{L}}_{SM}\right\}\mathrm{,}$ (29)

${\mathcal{L}}_{\text{grav}}=-\left(\frac{{M}_{P}^{2}}{2\left(8\pi \right)}+{\xi}_{h}{\varphi}^{\u2020}\varphi +{\xi}_{S}{S}^{2}\right)R\mathrm{,}$ (30)

${\mathcal{L}}_{S}=+\frac{1}{2}{\partial}_{\mu}S\cdot {\partial}^{\mu}S-\frac{1}{2}{\stackrel{\xaf}{m}}_{S}^{2}{S}^{2}-\frac{{\lambda}_{S}}{4}{S}^{4}\mathrm{,}$ (31)

${\mathcal{L}}_{S\varphi}=-\frac{1}{2}{\lambda}_{hS}\left({\varphi}^{\u2020}\varphi \right){S}^{2}\mathrm{.}$ (32)

Note that the Higgs field
$\varphi $ and dark matter field
$S$ have non-minimal couplings
${\xi}_{h}$ and
${\xi}_{S}$ to the curvature scalar *R*. These terms are needed for inflation (see Table 2 and Table 3) and also for renormalization [15] . The Higgs part of the Standard Model Lagrangian
${\mathcal{L}}_{SM}$ , in the unitary gauge with

$\varphi =\frac{1}{\sqrt{2}}\left(\begin{array}{c}0\\ {v}_{h}+h\end{array}\right)\mathrm{,}$ (33)

is

${\mathcal{L}}_{h}=+\frac{1}{2}{\partial}_{\mu}h\cdot {\partial}^{\mu}h-\frac{{\lambda}_{h}}{4}{h}^{4}-{\lambda}_{h}{v}_{h}{h}^{3}-{\lambda}_{h}{v}_{h}^{2}{h}^{2}\mathrm{.}$ (34)

Table 2. Parameters of the metric version of general relativity, obtained from
${n}_{s}=0.965\pm 0.004$ ,
${P}_{\zeta}\left({k}_{\text{ref}}\right)\equiv {\Delta}_{R}^{2}=\mathrm{exp}\left(3.044\pm 0.014\right)\times {10}^{-10}$ [1] , and
$k=0.05\text{\hspace{0.17em}}{\text{Mpc}}^{-1}$ , and Equations (42) to (49). The Higgs boson *h* is assumed to drive inflation. For
$S$ driving inflation replace
$h\to S$ . Uncertainties have 68% confidence, and are correlated.

Table 3. Parameters of the Palatini version of general relativity, obtained from
${n}_{s}=0.965\pm 0.004$ ,
${P}_{\zeta}\left({k}_{\text{ref}}\right)=\mathrm{exp}\left(3.044\pm 0.014\right)\times {10}^{-10}$ [1] , and
$k=0.05\text{\hspace{0.17em}}{\text{Mpc}}^{-1}$ , and Equations (42) to (49). Other parameters remain unconstrainted. The Higgs boson *h* is assumed to drive inflation. For
$S$ driving inflation replace
$h\to S$ . Uncertainties have 68% confidence, and are correlated.

The dark matter particle mass is ${M}_{S}=\sqrt{{\stackrel{\xaf}{m}}_{S}^{2}+\frac{1}{2}{\lambda}_{hS}{v}_{h}^{2}}=150\pm 2\text{\hspace{0.17em}}\text{eV}$ [4] , with ${\stackrel{\xaf}{m}}_{S}^{2}$

assumed positive, while ${\lambda}_{hS}$ is negative. The Higgs boson mass is ${m}_{H}=\sqrt{2{\lambda}_{h}}{v}_{h}\approx 125\text{\hspace{0.17em}}\text{GeV}$ . These masses play no significant role in inflation.

In this section we will assume that the field *h* drives inflation, and therefore set
$S=0$ initially, and so neglect the term with
${\xi}_{S}$ .

The Lagrangian can be brought to the “Einstein form”, without the non-minimal couplings to *R*, with a transformation of the metric [15]

${\stackrel{^}{g}}_{\mu \nu}={\Omega}^{2}{g}_{\mu \nu}\mathrm{,}\text{\hspace{1em}}{\Omega}^{2}=1+{\xi}_{h}{h}^{2}\frac{8\pi}{{M}_{P}^{2}}\mathrm{,}$ (35)

and a redefinition of the field. We distinguish two versions of General Relativity: the “metric” version, and the “Palatini” version. In the “Palatini” version of General Relativity the metric ${g}_{\mu \nu}$ and the affine connection ${\Gamma}_{\beta \gamma}^{\alpha}$ are treated as independent variables. The metric and Palatini versions differ if the non-minimal couplings to gravity, ${\xi}_{h}$ or ${\xi}_{S}$ , are different from zero. In the metric version, the needed field redefinition and the transformed potential, for $\chi \gg \sqrt{6/\left(8\pi \right)}{M}_{P}$ , are [1] [15] :

$h\approx \frac{{M}_{P}}{\sqrt{8\pi {\xi}_{h}}}\mathrm{exp}\left(\frac{\chi \sqrt{8\pi}}{\sqrt{6}{M}_{p}}\right)\mathrm{,}$ (36)

$U\left(\chi \right)\approx \frac{{\lambda}_{h}}{4{\xi}_{h}^{2}}{\left(\frac{{M}_{P}^{2}}{8\pi}\right)}^{2}{\left[1+\mathrm{exp}\left(-\frac{2\chi \sqrt{8\pi}}{\sqrt{6}{M}_{P}}\right)\right]}^{-2}$ (37)

In the Palatini version, the needed field redefinition and the transformed potential are [16] :

$h=\frac{{M}_{P}}{\sqrt{8\pi {\xi}_{h}}}\mathrm{sinh}\left(\frac{\sqrt{8\pi {\xi}_{h}}\chi}{{M}_{P}}\right)\mathrm{.}$ (38)

$U\left(\chi \right)=\frac{{\lambda}_{h}}{4{\xi}_{h}^{2}}{\left(\frac{{M}_{P}^{2}}{8\pi}\right)}^{2}{\left[\mathrm{tanh}\left(\frac{\sqrt{8\pi {\xi}_{h}}\chi}{{M}_{P}}\right)\right]}^{4}\mathrm{.}$ (39)

The slow-roll parameters are defined as

$\epsilon \equiv \frac{{M}_{P}^{2}}{16\pi}{\left(\frac{\text{d}U/\text{d}\chi}{U}\right)}^{2}\approx -\frac{\text{d}H/\text{d}t}{{H}^{2}}\mathrm{,}$ (40)

$\eta \equiv \frac{{M}_{P}^{2}}{8\pi}\frac{{\text{d}}^{2}U/\text{d}{\chi}^{2}}{U}\mathrm{.}$ (41)

In the slow-roll approximation, while $\epsilon \ll 1$ and $\eta \ll 1$ , in the metric [15] or Palatini [16] versions of General Relativity, we obtain

$\epsilon \approx \frac{4}{3{\xi}_{h}^{2}{h}^{4}}{\left(\frac{{M}_{P}^{2}}{8\pi}\right)}^{2}\text{\hspace{1em}}\text{or}\text{\hspace{1em}}\approx \frac{8}{{\xi}_{h}{h}^{4}}{\left(\frac{{M}_{P}^{2}}{8\pi}\right)}^{2}\mathrm{,}$ (42)

$\eta \approx -\frac{4}{3{\xi}_{h}{h}^{2}}\left(\frac{{M}_{P}^{2}}{8\pi}\right)\text{\hspace{1em}}\text{or}\text{\hspace{1em}}\approx -\frac{8}{{h}^{2}}\frac{{M}_{P}^{2}}{8\pi}\mathrm{,}$ (43)

$N\approx \frac{3}{4}{\xi}_{h}\left({h}_{i}^{2}-{h}_{f}^{2}\right)\left(\frac{8\pi}{{M}_{P}^{2}}\right)\text{\hspace{1em}}\text{or}\text{\hspace{1em}}\approx \frac{\pi}{{M}_{P}^{2}}\left({h}_{i}^{2}-{h}_{f}^{2}\right)\mathrm{,}$ (44)

${n}_{s}-1\equiv \frac{\text{d}\mathrm{ln}{P}_{\zeta}}{\text{d}\mathrm{ln}k}\approx -\frac{2}{N}\text{\hspace{1em}}\text{or}\text{\hspace{1em}}\approx -\frac{2}{N}\mathrm{,}$ (45)

$r\equiv \frac{{P}_{t}\left({k}_{\text{ref}}\right)}{{P}_{\zeta}\left({k}_{\text{ref}}\right)}\approx \frac{12}{{N}^{2}}\text{\hspace{1em}}\text{or}\text{\hspace{1em}}\approx \frac{2}{{\xi}_{h}{N}^{2}}\mathrm{,}$ (46)

${\rho}_{i}={U}_{i}=\frac{{\lambda}_{h}}{4{\xi}_{h}^{2}}{\left(\frac{{M}_{P}^{2}}{8\pi}\right)}^{2}\text{\hspace{1em}}\text{or}\text{\hspace{1em}}=\frac{{\lambda}_{h}}{4{\xi}_{h}^{2}}{\left(\frac{{M}_{P}^{2}}{8\pi}\right)}^{2}\mathrm{,}$ (47)

${\rho}_{f}\approx {U}_{f}\approx 0.134{U}_{i}\text{\hspace{1em}}\text{or}\text{\hspace{1em}}\approx {\left[\mathrm{tanh}\left(\sqrt{8{\xi}_{h}}\right)\right]}^{4}{U}_{i}\mathrm{,}$ (48)

${P}_{\zeta}\left({k}_{\text{ref}}\right)=\frac{{H}_{i}^{2}}{\pi {\epsilon}_{i}{M}_{P}^{2}}\approx \frac{{\lambda}_{h}}{2}{\left(\frac{{h}_{i}}{{M}_{P}}\right)}^{4}\text{\hspace{1em}}\text{or}\text{\hspace{1em}}\approx \frac{{\lambda}_{h}}{12{\xi}_{h}}{\left(\frac{{h}_{i}}{{M}_{P}}\right)}^{4}\mathrm{.}$ (49)

Let us consider the metric version of General Relativity. Note that the measured
${n}_{s}$ obtains *N*, *r*,
$\sqrt{{\xi}_{h}}{h}_{i}$ ,
${\epsilon}_{i}$ , and
${\eta}_{i}$ , see Table 2. The measured
${P}_{\zeta}\left(k\right)$ obtains
${\lambda}_{h}^{1/4}{h}_{i}$ . Together they determine
${\lambda}_{h}/{\xi}_{h}^{2}$ ,
${\rho}_{i}$ ,
${\rho}_{f}$ ,
${H}_{i}$ , and *H _{f}*. Finally,

The preceding results, with the interchange $S\leftrightarrow h$ , are also valid for inflation driven by $S$ instead of $\varphi $ . We note that ${v}_{h}$ and ${M}_{S}$ have no significant effect on inflation.

5. Comments on Higgs Inflation

The Higgs quartic coupling ${\lambda}_{h}$ runs from ${\lambda}_{h}\left({m}_{t}\right)\approx 0.126$ at ${m}_{t}$ to

$\begin{array}{c}{\lambda}_{h}\left({M}_{P}\right)=-0.0143-0.0066\left(\frac{{m}_{t}}{\text{GeV}}-173.34\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}+0.0018\frac{{\alpha}_{3}\left({M}_{Z}\right)-0.1184}{0.0007}+0.0029\left(\frac{{m}_{H}}{\text{GeV}}-125.15\right)\mathrm{,}\end{array}$ (50)

at
${M}_{P}$ [17] . Note that, within 3*σ* uncertainties,
${\lambda}_{h}\left({M}_{P}\right)$ may be negative, zero or positive! The relation
${\lambda}_{h}\left({M}_{P}\right)\approx 0$ is not a prediction of the Standard Model since it depends on input parameters to the model, e.g.
${\lambda}_{h}$ ,
${v}_{h}$ and
${y}_{t}$ . Either
${\lambda}_{h}\left({M}_{P}\right)\approx 0$ is a coincidence, or some unknown high energy extension of the Standard Model relates these parameters. Since
${\lambda}_{h}\left({M}_{P}\right)$ may be positive up to
${M}_{P}$ , the possible need for new physics is pushed to very high energies. Note that the Higgs mass
${m}_{H}=\sqrt{2{\lambda}_{h}}{v}_{h}$ has relatively little running [17] [18] , so
${v}_{h}^{2}$ running compensates
${\lambda}_{h}$ running.
${\beta}_{\lambda}$ also vanishes around the Planck mass. “Near-criticality holds for both the Higgs quartic and the top Yukawa coupling” [17] .

From [17] and [19] we conclude that Electro Weak Symmetry Breaking (EWSB) at ${T}_{\text{EWSB}}=159.5\pm 1.5\text{\hspace{0.17em}}\text{GeV}$ [19] is not due to ${v}_{h}^{2}$ running to zero at this temperature, but to a fraction of the Higgs bosons excited to above the potential

barrier ${\left(\frac{1}{4}{\lambda}_{h}\right)}^{1/4}{v}_{h}={m}_{H}/{\left(16{\lambda}_{h}\right)}^{1/4}$ . We note that this potential barrier grows

with the renormalization energy scale $\mu $ slower than the temperature $T=\mu $ , which allows the restoration of the electroweak symmetry $U{\left(1\right)}_{Y}\times SU{\left(2\right)}_{L}$ at $T>{T}_{\text{EWSB}}$ .

The renormalization energy scale at which ${\lambda}_{h}$ crosses zero depends critically on the Yukawa coupling ${y}_{t}$ , on the SU(3) gauge coupling ${g}_{3}$ , and, to a lesser extent, on ${g}_{2}$ . During inflation the Higgs boson does not have time to become “dressed”, so we should use ${\lambda}_{h}\left({M}_{P}\right)$ (not ${\lambda}_{h}\left({m}_{t}\right)$ ) which depends critically on ${m}_{t}$ , ${g}_{3}$ , and ${m}_{H}$ , and may even be zero or negative, and so, Higgs inflation may break down.

The renormalization group equations for the Standard Model extended with a scalar field coupled to the Higgs boson, have been studied in [20] [21] and [22] . Our conclusions, for the warm dark matter extension of the Standard Model studied in this article, with the expected
${\lambda}_{hS}$ of order 10^{−8} to 10^{−7}, is that the running of Standard Model parameters is not significantly changed, while
${\lambda}_{S}\left(\mu \right)$ increases with
$\mu $ , and so does not become zero. For example, if we set
${\lambda}_{S}\left({M}_{Z}\right)=0.126$ , we obtain
${\lambda}_{S}\left({M}_{P}\right)\approx 0.5$ .

6. Preheating after Inflation

If decays of the inflaton field can be neglected during inflation, then ${T}^{3}{a}^{3}$ is constant, and the temperature is exponentially small, or zero, at the end of this “cold” inflation. A “reheating” period is needed after the end of cold inflation at ${a}_{f}$ in order to start the radiation-dominated hot big-bang cosmology. After inflation “the field oscillates about a potential minimum and decays. How this occurs is unknown in any detail.” [23] Steven Weinberg makes a similar comment [14] . This oscillation and decay is known as “preheating”. For a potential $V\left(h\right)\propto {h}^{n}$ , the equation of state is $w\equiv p/\rho =\left(n-2\right)/\left(n+2\right)$ [24] . So, while the oscillation is dominated by the quartic potential, the equation of state is $w=1/3$ , as for radiation. If the potential becomes quadratic, $w=0$ , and the inflaton field expands as matter (but this will not be the case in the present scenario).

Assume $S$ is the inflaton field. During preheating $S$ oscillates anharmonically with angular frequency $\omega \left(t\right)$ [25] (defined as 2π divided by the period), and

$\rho ={H}^{2}\frac{3{M}_{P}^{2}}{8\pi}=\frac{1}{{\left(2t\right)}^{2}}\frac{3{M}_{P}^{2}}{8\pi}={\rho}_{f}{\left(\frac{{a}_{f}}{a}\right)}^{4}\approx \frac{{\omega}^{2}{S}_{\text{peak}}^{2}}{2},$ (51)

$\rho \equiv \frac{{\pi}^{2}}{30}{g}_{\text{reh}}{T}^{4}\approx \frac{{\pi}^{2}}{30}107.75{\left(\frac{0.332{T}_{0}}{a}\right)}^{4},$ (52)

where $T$ is a pseudo-temperature since thermal equilibrium has not yet been reached. Note that in the metric version of General Relativity ${a}_{f}$ and ${\rho}_{f}$ are obtained from measurements of ${n}_{s}$ and ${P}_{\zeta}\left(k\right)$ , and so all parameters are known during preheating. Somehow we need to create particle pairs with opposite momenta $\omega $ and $-\omega $ . To this end we need density and/or curvature perturbations that couple to the oscillating inflaton field.

Local thermal and chemical equilibrium needs to be achieved before the temperature
$T$ drops below
${T}_{\text{BBN}}\approx 1$ MeV to not upset the predictions of Big Bang Nucleosynthesis (BBN) [1] . Thermal and chemical equilibrium between
$S$ and *h* needs to be achieved before
$T$ drops below *m _{H}*, so that dark matter agrees with observations [4] .

6.1. Preheating after Higgs Inflation

This subsection is presented to illustrate issues with preheating. The evolution of the coherent inflaton field is given by (16) (where in the case of Higgs inflation we need to replace
$S$ by *h*). When the expansion rate
$H\propto \sqrt{\rho}\propto {a}^{-2}\propto {t}^{-1}$ drops below the decay rate Γ, *h* decays effectively and, with (9), defines a preheating pseudo-temperature
${T}_{\text{pre}}$ (“pseudo” because thermal equilibrium has not yet been reached) [1] :

${T}_{\text{pre}}\approx 0.2{\left(\frac{100}{{g}_{\text{reh}}}\right)}^{1/4}\sqrt{{M}_{P}{\Gamma}_{\text{tot}}}\mathrm{.}$ (53)

The theoretical total width of the Higgs boson is
${\Gamma}_{\text{tot}}=4.07\times {10}^{3}\text{\hspace{0.17em}}\text{GeV}$ [1] (neglecting decays to *SS*), so a very rough first estimate is

${T}_{\text{pre}}\approx {10}^{8}\text{\hspace{0.17em}}\text{GeV}\mathrm{.}$ (54)

However, it is not a Higgs boson *h* with mass *m _{H}* (defined by a quadratic potential) that decays, it is a coherent Higgs boson field (

$Y\stackrel{\u02dc}{L}hR$ ,
$g{g}^{\prime}h{v}_{h}V{V}^{\prime}$ ,
$g{g}^{\prime}{h}^{2}V{V}^{\prime}$ ,
$\frac{1}{2}{\lambda}_{hS}h{v}_{h}{S}^{2}$ ,
$\frac{1}{4}{\lambda}_{hS}{h}^{2}{S}^{2}$ , and
$\frac{1}{4}{\lambda}_{h}{h}^{4}$ , where *L*

(*R*) are left-handed (right-handed) quarks and leptons, and *V* and *V'* are gauge bosons
${W}_{\mu}$ ,
${B}_{\mu}$ . At
$T>{T}_{\text{EWSB}}$ the electroweak symmetry
$SU{\left(2\right)}_{L}\times U{\left(1\right)}_{Y}$ is restored, so
${v}_{h}$ becomes effectively zero, and several decay channels are closed. The terms
$Y\stackrel{\u02dc}{L}hR$ can be neglected because the Yukawa couplings
$Y\ll {T}_{\text{pre}}/{v}_{h}$ . The surviving terms are quadratic in the oscillating inflaton field (*h*), and need a special treatment.

There are several other issues [26] . For example, if there are already $n$ bosons in some final state, then the decay rate to this state needs to be multiplied by $n+1$ . Also, since ${\lambda}_{h}\left({M}_{P}\right)$ at Planck energies is very small (or even negative!) Higgs inflation may break down, and $S$ -inflation may need to come to the rescue.

6.2. Preheating by Parametric Resonance

Consider inflation driven by the field
$S$ . Inflation ends, *i.e.*
$\stackrel{\xa8}{a}=0$ , when

$\epsilon \approx -\stackrel{\dot{}}{H}/{H}^{2}=1$ . Thereafter the homogeneous classical field $S$ oscillates coherently around the minimum of the potential $V\left(S\right)$ . We assume that the quartic term dominates the potential $V\left(S\right)$ . The problem of preheating is how to convert the homogeneous classically oscillating coherent field into quantum mechanical particles. The classical equations to be solved, while $\Gamma <H$ , are (1) with

$\rho =\frac{1}{2}{\stackrel{\dot{}}{S}}^{2}+\frac{1}{4}{\lambda}_{S}{S}^{4}$ , and (16). The term in the Lagrangian (30) with ${\xi}_{S}$ , that is

important early on during inflation, can be neglected when inflation ends and $S$ oscillates. In our case of interest of warm dark matter, we also neglect the term

with ${\stackrel{\xaf}{m}}_{S}^{2}$ . Then the potential of interest is $\frac{1}{4}{\lambda}_{S}{S}^{4}+\frac{1}{4}{\lambda}_{hS}{h}^{2}{S}^{2}$ . This “conformally

invariant” potential has been studied in detail in [25] . The analytical solution obtains $\rho \propto {a}^{-4}$ , as in a radiation-dominated universe. The oscillation has an amplitude that decays as ${S}_{\mathrm{max}}\propto {a}^{-1}\propto {t}^{-1/2}$ due to the expansion of the universe, and the angular frequency of the non-sinusoidal oscillation $\omega \equiv 2\pi /\tau $ , where $\tau $ is the period, becomes redshifted:

$\omega =\frac{2\pi}{7.416}{\left(\frac{\pi}{6{\lambda}_{S}{M}_{P}^{2}}\right)}^{1/4}{t}^{-1/2}\sqrt{{\lambda}_{S}}{M}_{P}\mathrm{.}$ (55)

The classically oscillating field (
$S$ ) decays into
$S$ and *h* particles that are treated quantum-mechanically. Let
${S}_{k}$ and *h _{k}* be the quantum mechanical amplitudes of the particles

${\stackrel{\xa8}{h}}_{k}+3\frac{\stackrel{\dot{}}{a}}{a}{\stackrel{\dot{}}{h}}_{k}+\left(\frac{{k}^{2}}{{a}^{2}}+\frac{1}{2}{\lambda}_{hS}{S}^{2}\right){h}_{k}=\mathrm{0,}$ (56)

${\stackrel{\xa8}{S}}_{k}+3\frac{\stackrel{\dot{}}{a}}{a}{\stackrel{\dot{}}{S}}_{k}+\left(\frac{{k}^{2}}{{a}^{2}}+3{\lambda}_{S}{S}^{2}\right){S}_{k}=\mathrm{0,}$ (57)

where
$k$ is the comoving momentum, and
$k/a$ is the proper (or physical) momentum, of the mode. Dots indicate derivatives with respect to time *t*. The terms proportional to
${S}^{\text{2}}$ in (56) and (57) come from
${\text{d}}^{2}V/\text{d}{h}^{2}$ or
${\text{d}}^{2}V/\text{d}{S}^{2}$ , respectively. These terms are the source of the parametric resonance. (Parametric resonance is understood by every child on a swing: the child stands up every time the swing passes its lowest point!) Note that both terms inside the parentheses are proportional to
${a}^{-\text{2}}$ , which allows an elimination of
$a$ by changes of variables. The solutions of these equations depend on the parameter combination
${\lambda}_{hS}/{\lambda}_{S}$ , and on the time independent comoving momentum
$k$ of the modes. The oscillations of
$S$ cause parametric resonances with numbers of particles *n _{S}* and

${n}_{h}\approx \frac{1}{2}{\text{e}}^{2{\mu}_{h}x}\text{\hspace{1em}}\text{with}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}x={\left(\frac{6{\lambda}_{S}{M}_{P}^{2}}{\pi}\right)}^{1/4}{t}^{1/2}\mathrm{,}$ (58)

and
$0<{\mu}_{h}<0.2377$ . Note that the exponential growth of *n _{h}* starts with the vacuum initial quantum field amplitude for the field

This is only half of the story: we have reviewed the increase of *n _{h}*, but not the decrease of the coherent inflaton field (
$S$ ). This “backreaction” is studied in [25] .

7. Reheating

For more general presentations see [27] and [26] . Here we will study two channels that need to bring warm dark matter
$S$ into thermal and diffusive equilibrium with the Higgs boson before the temperature drops below *m _{H}*. This is a requirement of cosmological observations [4] .

7.1. Thermalization of $S$ by Higgs Boson Annihilations $hh\to SS$

Consider Higgs inflation followed by Standard Model reheating to a temperature
${T}_{\text{reh}}\equiv {T}_{1}$ . Here we study thermalization of
$S$ by Higgs boson annihilation
$hh\to SS$ , which must occur before the temperature drops below
${T}_{2}={m}_{H}$ to satisfy the “no freeze-in and no freeze-out” warm dark matter scenario. The Higgs boson is assumed to remain in thermal equilibrium. The universe evolves as (1) with density (9). The expansion parameter
$a\left(t\right)$ , before
$T$ drops below *m _{t}*, is related to temperature
$T$ by (8). The increase of the number of dark matter particles
$S$ per unit time and unit comoving volume is

$\frac{\text{d}\left({n}_{S}{a}^{3}\right)}{\text{d}t}=2\sigma {n}_{h}^{2}{a}^{3}$ , with (59)

$\sigma \left(hh\to SS\right)\approx \frac{{\lambda}_{hS}^{2}}{256\pi {\left(2\cdot 2.7T\right)}^{2}}$ , and (60)

${n}_{h}=0.1218\cdot {T}^{3}$ . (61)

From the above equations, integrating from
${n}_{S1}{a}_{1}^{3}$ and
${T}_{\text{1}}$ , to
${n}_{S2}{a}_{2}^{3}$ and
${T}_{\text{2}}$ , neglected
${n}_{S1}{a}_{1}^{3}$ (to be conservative) and
$\text{1}/{T}_{\text{1}}$ , and setting
${n}_{S2}={n}_{h2}$ for thermal and diffusive equilibrium, we obtain the coupling constant
${\lambda}_{hS}$ needed to obtain thermal and diffusive equilibrium between
$S$ and *h* at temperature
${T}_{\text{2}}$ with this annihilation channel only:

${\lambda}_{hS}=-1.3\times {10}^{3}{\left(\frac{{T}_{2}}{{M}_{P}}\right)}^{1/2}\mathrm{.}$ (62)

For example, for ${T}_{2}=246\text{\hspace{0.17em}}\text{GeV}$ , we need ${\lambda}_{hS}\approx -5.8\times {10}^{-6}$ . For ${T}_{2}=125\text{\hspace{0.17em}}\text{GeV}$ , we need ${\lambda}_{hS}\approx -4.1\times {10}^{-6}$ . A similar result is obtained for the channel $\sigma \left({W}^{+}{W}^{-}\to {h}^{\mathrm{*}}\to SS\right)$ [10] .

7.2. Thermalization of $S$ by Higgs Boson Decays $h\to SS$

For Higgs boson decay to $SS$ , the decay probability per unit time below ${T}_{\text{EWSB}}=160\text{\hspace{0.17em}}\text{GeV}$ , is

$\Gamma \left(h\to SS\right)=\frac{{\lambda}_{hS}^{2}{v}_{h}^{2}}{64\pi {m}_{H}}\mathrm{.}$ (63)

We set
$\Gamma \approx H$ to estimate the temperature of equilibrium. The coupling
${\lambda}_{hS}$ required to obtain thermal and diffusive equilibrium of
$S$ with *h* at
${T}_{2}={m}_{H}$ , assuming this is the only channel, is
${\lambda}_{hS}\approx -1\times {10}^{-7}$ .

7.3. Self-Interacting Dark Matter $S$

There is mounting evidence that dark matter is self-interacting: see [28] and the extensive list of references therein. From observations of galaxy NGC5044, the following estimate of the dark matter-dark matter self-interaction cross-section per unit mass is obtained: ${\sigma}_{\text{DM-DM}}/{m}_{\text{DM}}=0.165\pm 0.025\text{\hspace{0.17em}}{\text{cm}}^{\text{2}}/\text{g}$ (for a dark matter velocity dispersion of about 300 km/s) [28] . We consider two interaction channels, $SS\to SS$ (that depends on ${\lambda}_{S}$ ), and $SS\to {h}^{\mathrm{*}}\to SS$ (that depends on ${\lambda}_{hS}$ ), and assume one dominates the other, so we neglect their interference. If $SS\to SS$ dominates, we obtain ${\lambda}_{S}\approx 4.5\times {10}^{-8}$ . If $SS\to {h}^{\mathrm{*}}\to SS$ dominates, we obtain ${\lambda}_{hS}\approx -8.8\times {10}^{-5}$ . The coupling $\left|{\lambda}_{hS}\right|$ needs to be less than 0.03 to not exceed the limit on the Higgs boson invisible decay width [10] . Calculating $\sigma \left(SS\to {h}^{\mathrm{*}}\to \stackrel{\u02dc}{\nu}\nu \right)$ , and $n=0.1218\cdot {m}_{\text{DM}}^{3}$ , we verify that $t=1/\left(\sigma cn\right)$ is much greater than the age of the universe.

8. Conclusions

We have extended the Standard Model with a scalar warm dark matter field
$S$ that couples to the Higgs boson
$\varphi $ . This scenario is in agreement with recent cosmological observations if
$S$ and
$\varphi $ come into thermal and diffusive equilibrium before the temperature drops below *m _{H}* [4] . The extended model has two scalar fields,
$S$ and
$\varphi $ . Either of these fields may drive inflation. We study several inflation options, and arrive at the following conclusions:

1) $S$ -inflation with quadratic potential: If the field that drives inflation is a

scalar $S$ with a quadratic potential $\frac{1}{2}{\stackrel{\xaf}{m}}_{S}^{2}{S}^{2}$ , we can satisfy the observed inflation

constraints ${n}_{s}$ and ${P}_{\zeta}\left(k\right)$ with a mass ${\stackrel{\xaf}{m}}_{S}\approx 1.2\times {10}^{-6}{M}_{P}$ . This mass corresponds to cold dark matter, and so $S$ would not be the warm dark matter favored by observations.

2) $S$ -inflation with quartic potential: If the field that drives inflation is a scalar

$S$ with a dominating quartic potential $\frac{1}{4}{\lambda}_{S}{S}^{4}$ , and a non-minimal coupling

$-{\xi}_{S}{S}^{2}R$ to the curvature scalar
$R$ , we can satisfy the observed inflation constraints
${n}_{s}$ and
${P}_{\zeta}\left(k\right)$ if
${\lambda}_{S}/{\xi}_{S}^{2}=\left(4.6\pm 1.1\right)\times {10}^{-10}$ for the metric version of General Relativity (see Table 2), or
${\lambda}_{S}/{\xi}_{S}=\left(7.6\pm 1.8\right)\times {10}^{-11}$ for the Palatini version (see Table 3). Thermal and diffusive equilibrium between
$S$ and *h* is achieved at
$T={m}_{H}$ if
$\left|{\lambda}_{hS}\right|\gtrsim {10}^{-7}$ .
$S$ -inflation does not suffer from critical
${\lambda}_{S}$ running at Planck energies. For the metric version of General Relativity, the predicted tensor-to-scalar ratio is
$r=0.0037\pm 0.0009$ , while it is unconstrained in the Palatini version. In summary, this scenario passes current observational tests.

3) Complex field
$S$ : An alternative to the real dark matter field
$S$ with *Z*_{2} symmetry, is a complex
$S$ with a Lagrangian with
$U{\left(1\right)}_{S}$ symmetry, with the ground state of
$S$ and
$\varphi $ breaking the
$U{\left(1\right)}_{S}$ and
$U{\left(1\right)}_{Y}\times SU{\left(2\right)}_{L}$ symmetries [10] . In this scenario,
$S$ decays to two massive gauge bosons *V* that are the dark matter. The present analysis is also valid for this complex
$S$ field with minor notational modifications, e.g.
${S}^{2}\to {S}^{\mathrm{*}}S$ .

4) Higgs inflation: Higgs inflation is similar to
$S$ -inflation with quartic potential. The main difference is an uncertain quartic self-coupling
${\lambda}_{h}\left({M}_{P}\right)$ at the Planck energy scale, that, within 3*σ* uncertainties, may even be zero or negative, in which case Higgs inflation breaks down!

5) The mass *M _{S}* of
$S$ , and the Higgs vacuum expectation value

The closure test presented in Figure 1, obtained from observations of *n _{s}* and
${P}_{\zeta}\left(k\right)$

Appendix

A Pedestrian Derivation of ${P}_{\zeta}(\; k\; )$

The purpose of this appendix is to try to understand the quantum origin of the density perturbations that are observed in the Cosmic Microwave Background radiation, which are the seeds of galaxy formation and of the large-scale structure of the universe. This derivation is outlined in [23] . For rigorous derivations see [1] [13] , or [14] .

Consider the comoving length scale $\text{1}/k$ that exits the horizon at ${a}_{i}$ , and re-enters the horizon at ${a}_{r}$ : $k={a}_{i}{H}_{i}={a}_{r}{H}_{r}$ . See Figure 1. We assume that during inflation, and well within the horizon, the field $S$ has a homogeneous component, plus random quantum fluctuations $\delta S$ . We apply periodic boundary conditions to a cube of volume ${L}^{3}$ , and consider the lowest non-zero mode with $\omega =2\pi /L$ , and assume that the energy in the fluctuating field in this mode

is $\frac{1}{2}\omega $ , as in the quantum harmonic oscillator. We neglect masses, so take $\omega =k$ .

This energy is
${L}^{3}{\langle {\stackrel{\dot{}}{S}}^{2}\rangle}_{k}={L}^{3}{\omega}^{2}{\langle \delta {S}^{2}\rangle}_{k}$ , so
${\langle \delta {S}^{2}\rangle}_{k}=1/\left(4\pi {L}^{2}\right)$ . As the universe expands, the box expands as
$L\propto a$ ,
$\omega =k\propto {a}^{-1}$ and
$\langle \delta {S}^{2}\rangle \propto {a}^{-2}$ . When the length scale
$\text{1}/k$ exits the horizon,
$L=2\pi /k$ , and
${H}_{i}\approx 1/L$ , so, within a factor of *O*(1), we obtain the result of [1] for the power spectrum of vacuum field fluctuations at horizon exit
${a}_{i}$ :

${\langle \delta {S}_{i}^{2}\rangle}_{k}={\left(\frac{{H}_{i}}{2\pi}\right)}_{k}^{2}\mathrm{.}$ (64)

We now consider a “small” volume well within the horizon. The density and field evolve as (5) and (16), which neglect spatial derivatives. Inflation ends at ${\epsilon}_{f}=1$ , so the density ${\rho}_{f}$ is independent of $\delta S$ . The fluctuation ${\langle \delta {S}_{i}^{2}\rangle}_{k}$ , at ${a}_{i}$ during slow-roll inflation, causes a spread ${\langle \delta {t}_{f}^{2}\rangle}_{k}$ of the time ${t}_{f}$ of the end of inflation:

${\langle \delta {t}_{f}^{2}\rangle}_{k}=\frac{{\langle \delta {S}_{i}^{2}\rangle}_{k}}{{\left(\text{d}S/\text{d}t\right)}_{i}^{2}}={\langle \delta {S}_{i}^{2}\rangle}_{k}{\left(\frac{3{H}_{i}}{-\text{d}V/\text{d}S}\right)}^{2}$ (65)

$={\left(\frac{{H}_{i}}{2\pi}\right)}_{k}^{2}{\left(\frac{3{H}_{i}}{\text{d}V/\text{d}S}\right)}^{2}=\frac{1}{\pi {\epsilon}_{i}{M}_{P}^{2}}\mathrm{.}$ (66)

(Drawing a graph $V\left(t\right)$ helps.)

We must now pass from slow roll inflation to expansion with $\rho \propto {t}^{-2}$ characteristic of radiation (or matter), or of oscillations of $S$ in a quartic potential at the end of inflation. The time fluctuation ${\langle \delta {t}_{f}^{2}\rangle}_{k}$ causes a density perturbation at a fixed time ${t}_{{f}^{\prime}}$ :

$\frac{\delta {\rho}_{{f}^{\prime}}}{{\rho}_{{f}^{\prime}}}=\frac{\text{d}{\rho}_{{f}^{\prime}}/\text{d}t}{{\rho}_{{f}^{\prime}}}\sqrt{{\langle \delta {t}_{f}^{2}\rangle}_{k}}=-\frac{2}{{t}_{{f}^{\prime}}}\sqrt{{\langle \delta {t}_{f}^{2}\rangle}_{k}}=-4{H}_{f}\sqrt{{\langle \delta {t}_{f}^{2}\rangle}_{k}}.$ (67)

For
$\rho \propto {t}^{-2}$ , and
$\delta \rho \propto {t}^{-2}$ ,
$\delta \rho /\rho $ is independent of *t*, so
$\delta {\rho}_{f}/{\rho}_{f}=\delta {\rho}_{r}/{\rho}_{r}$ . Also
$3{H}_{f}\approx {H}_{i}$ , see Table 2, so, within a factor *O*(1), we obtain the result of [1] :

${P}_{\zeta}\left(k\right)\equiv {\left(\frac{\langle \delta {\rho}_{r}^{2}\rangle}{{\rho}_{r}^{2}}\right)}_{k}={\left(\frac{{H}_{i}^{2}}{\pi {\epsilon}_{i}{M}_{P}^{2}}\right)}_{k}\mathrm{.}$ (68)

NOTES

^{1}“Thermal relic” warm dark matter with mass
$\gg 150\text{\hspace{0.17em}}\text{eV}$ freezes-out and decouples from the Standard Model sector while non-relativistic (to not exceed the measured dark matter density). Dark matter
$S$ is colder than this “thermal relic” by a factor
$<{T}_{S}/{T}_{\gamma}=0.345$ . Hence, to compare *M _{S}* with limits that can be found in the literature on the “thermal relic” dark matter mass it is necessary to multiply the latter by a factor <0.345

^{2}The dimensionless power specturm
$P\left(k\right)$ is defined as follows [1] .
$\delta \left(x\right)\equiv \left(\rho \left(x\right)-\langle \rho \rangle \right)/\langle \rho \rangle ={\displaystyle \sum {\delta}_{k}\mathrm{exp}\left(-ik\cdot x\right)}$ (with periodic boundary conditions in a box of volume
${L}^{3}$ ). The density variance is then
$\langle {\delta}^{2}\rangle ={\displaystyle \sum}{\left|{\delta}_{k}\right|}^{2}={\displaystyle \int}\frac{4\pi {k}^{2}\text{d}k}{{\left(2\pi /L\right)}^{3}}{\left|{\delta}_{k}\right|}^{2}\equiv {\displaystyle \int \text{d}\mathrm{ln}\left(k\right)P\left(k\right)}$ , so the dimensionless power spectrum is
$P\left(k\right)\equiv \frac{{k}^{3}}{2{\pi}^{2}}{L}^{3}{\left|{\delta}_{k}\right|}^{2}$ . In [13] ,
$P\left(k\right)$ is written as
${\Delta}_{R}^{2}\left(k\right)$ .

Conflicts of Interest

The author declares no conflicts of interest regarding the publication of this paper.

[1] |
Particle Data Group, Zyla, P.A., et al. (2020) Review of Particle Physics. Progress of Theoretical and Experimental Physics, 2020, 083C01. https://doi.org/10.1093/ptep/ptaa104 |

[2] |
Ema, J., Mukaida, K. and van de Vis, J. (2021) Renormalization Group Equations of Higgs-R2 Inflation. Journal of High Energy Physics, 2021, Article No. 109. https://arxiv.org/pdf/2008.01096.pdf https://doi.org/10.1007/JHEP02(2021)109 |

[3] |
Hamada, Y., Kawai, H., Oda, K. and Park, S.C. (2014) Higgs Inflation from Standard Model Criticality. https://arxiv.org/pdf/1408.4864.pdf |

[4] | Hoeneisen, B. (2023) A Data Driven Solution to the Dark Matter Problem. European Journal of Applied Sciences, 11, 473-481. |

[5] |
Hoeneisen, B. (2022) Measurement of the Dark Matter Velocity Dispersion with Dwarf Galaxy Rotation Curves. International Journal of Astronomy and Astrophysics, 12, 363-381. https://doi.org/10.4236/ijaa.2022.124021 |

[6] |
Hoeneisen, B. (2022) Measurement of the Dark Matter Velocity Dispersion with Galaxy Stellar Masses, UV Luminosities, and Reionization. International Journal of Astronomy and Astrophysics, 12, 258-272. https://doi.org/10.4236/ijaa.2022.123015 |

[7] |
Hoeneisen, B. (2022) Comments on Warm Dark Matter Measurements and Limits. International Journal of Astronomy and Astrophysics, 12, 94-109. https://doi.org/10.4236/ijaa.2022.121006 |

[8] |
Hoeneisen, B. (2023) A Study of Warm Dark Matter, the Missing Satellites Problem, and the UV Luminosity Cut-Off. International Journal of Astronomy and Astrophysics, 13, 25-38. https://doi.org/10.4236/ijaa.2023.131002 |

[9] |
Hoeneisen, B. (2020) What Is Dark Matter Made of? Presented at the 3rd World Summit on Exploring the Dark Side of the Universe. Guadeloupe Islands. https://inspirehep.net/files/7cfb2bf406baf315315e389e6eff3809 |

[10] |
Hoeneisen, B. (2021) Adding Dark Matter to the Standard Model. International Journal of Astronomy and Astrophysics, 11, 59-72. https://doi.org/10.4236/ijaa.2021.111004 |

[11] |
Hoeneisen, B. (2022) Warm Dark Matter and the Formation of First Galaxies. Journal of Modern Physics, 13, 932-948. https://doi.org/10.4236/jmp.2022.136053 |

[12] |
Hoeneisen, B. (2021) Active-Sterile Neutrino Oscillations and Leptogenesis. Journal of Modern Physics, 12, 1248-1266. https://doi.org/10.4236/jmp.2021.129077 |

[13] |
Baumann, D. (2012) TASI Lectures on Inflation. https://arxiv.org/pdf/0907.5424.pdf |

[14] | Weinberg, S. (2008) Cosmology. Oxford University Press, Oxford. |

[15] |
Bezrukov, F. and Shaposhnikov, M. (2008) The Standard Model Higgs Boson as the Inflaton. https://arxiv.org/pdf/0710.3755.pdf |

[16] |
Shaposhnikov, M., Shkerin, A. and Zell, S. (2021) Standard Model Meets Gravity: Electroweak Symmetry Breaking and Inflation. Physical Review D, 103, Article ID: 033006. https://arxiv.org/pdf/2001.09088.pdf https://doi.org/10.1103/PhysRevD.103.033006 |

[17] |
Buttazzo, D., et al. (2013) Investigating the Near-Criticality of the Higgs Boson. Journal of High Energy Physics, 12, 89. https://doi.org/10.1007/JHEP12(2013)089 |

[18] |
Bian, L. (2013) RGE of the Higgs Mass in the Context of the SM. https://arxiv.org/pdf/1303.2402.pdf |

[19] |
D’Onofrio, M. and Rummukainen, K. (2018) The Standard Model Cross-Over on the Lattice. https://arxiv.org/pdf/1508.07161.pdf |

[20] | Kawana, K. (2019) Multiple Point Principle of the Standard Model with Scalar Singlet Dark Matter and Right Handed Neutrinos. arXiv: 1411.2097. |

[21] | Haba, N., Kaneta, K. and Takahashi, R. (2013) Planck Scale Boundary Conditions in the Standard Model with Singlet Scalar Dark Matter. arXiv: 1312.2089. |

[22] |
Costa, R., Morais, A.P., Sampaio, M.O.P. and Santos, R. (2015) Two-Loop Stability of a Complex Singlet Extended Standard Model. Physical Review D, 92, Article ID: 025024. https://doi.org/10.1103/PhysRevD.92.025024 |

[23] |
Green, D. (2014) Inflation and the Higgs Scalar (Lecture Notes). https://arxiv.org/pdf/1412.2107.pdf |

[24] |
Maity, D. and Saha, P. (2019) CMB Constraints on Dark Matter Phenomenology via Reheating in Minimal Plateau Inflation. Physics of the Dark Universe, 25, Article ID: 100317. https://arxiv.org/pdf/1804.10115.pdf https://doi.org/10.1016/j.dark.2019.100317 |

[25] |
Greene, P.B., Kofman, L., Linde, A. and Starobinsky, A.A. (1997) Structure of Resonance in Preheating after Inflation. Physical Review D, 56, 6175. https://arxiv.org/pdf/hep-ph/9705347.pdf https://doi.org/10.1103/PhysRevD.56.6175 |

[26] |
Lozanov, K. (2019) Lectures on Reheating after Inflation. https://arxiv.org/pdf/1907.04402.pdf |

[27] |
Cai, R.G., Guo, Z.K. and Wang, S.J. (2018) Reheating Phase Diagram for Single-Field Slow-Roll Inflationary Models. https://arxiv.org/pdf/1501.07743.pdf |

[28] |
Gopica, K. and Desai, S. (2023) Constraints on Self-Interacting Dark Matter from Relaxed Galaxy Groups. Physics of the Dark Universe, 42, Article ID: 101291. https://arxiv.org/pdf/2307.05880.pdf |

Journals Menu

Contact us

+1 323-425-8868 | |

customer@scirp.org | |

+86 18163351462(WhatsApp) | |

1655362766 | |

Paper Publishing WeChat |

Copyright © 2024 by authors and Scientific Research Publishing Inc.

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.