# 1 fig1

FERMILAB–PUB–98/217–A   
CERN-TH/98-227   
hep-ph/9809453   
September 1998   

Production of massive particles during reheating

Daniel J. H. Chung,<sup>a,b,</sup><sup>1</sup><sup>1</sup>1Electronic mail: djchung@theory.uchicago.edu Edward W. Kolb,<sup>b,c,</sup><sup>2</sup><sup>2</sup>2Electronic mail: rocky@rigoletto.fnal.gov and Antonio Riotto<sup>d,</sup><sup>3</sup><sup>3</sup>3Electronic mail: riotto@nxth04.cern.ch<sup>,</sup><sup>4</sup><sup>4</sup>4On leave from Department of Theoretical Physics, University of Oxford, U.K.  

<sup>a</sup>Department of Physics and Enrico Fermi Institute   
The University of Chicago, Chicago, Illinois 60637-1433   
<sup>b</sup>NASA/Fermilab Astrophysics Center   
Fermilab National Accelerator Laboratory, Batavia, Illinois 60510-0500   
<sup>c</sup>Department of Astronomy and Astrophysics and Enrico Fermi Institute   
The University of Chicago, Chicago, Illinois 60637-1433   
<sup>d</sup>Theory Division, CERN, CH-1211 Geneva 23, Switzerland 

> What is commonly called the reheat temperature, $`T_{RH}`$, is not the maximum temperature obtained after inflation. The maximum temperature is, in fact, much larger than $`T_{RH}`$. As an application of this we consider the production of massive stable dark-matter particles of mass $`M_X`$ during reheating, and show that their abundance is suppressed as a power of $`T_{RH}/M_X`$ rather than $`\mathrm{exp}(M_X/T_{RH})`$. We find that particles of mass as large as $`2\times 10^3`$ times the reheat temperature may be produced in interesting abundance. In addition to dark matter, our analysis is relevant for baryogenesis if the baryon asymmetry is produced by the baryon (or lepton) number violating decays of superheavy bosons, and also for relic ultra-high energy cosmic rays if decays of superheavy particles are responsible for the highest energy cosmic rays.
> 
> PACS number(s): 98.80.Cq

I. INTRODUCTION

At the end of inflation  the energy density of the universe is locked up in a combination of kinetic energy and potential energy of the inflaton field, with the bulk of the inflaton energy density in the zero-momentum mode of the field. Thus, the universe at the end of inflation is in a cold, low-entropy state with few degrees of freedom, very much unlike the present hot, high-entropy universe. After inflation the frozen inflaton-dominated universe must somehow be defrosted and become a high-entropy radiation-dominated universe.

The process by which the inflaton energy density is converted to radiation is known as “reheating” . The possible role of nonlinear dynamics leading to explosive particle production has recently received a lot of attention. This process, known as “preheating”  may convert a fair fraction of the inflaton energy density into other degrees of freedom, with extremely interesting cosmological effects such as symmetry restoration, baryogenesis, or production of dark matter. But the efficiency of preheating is very sensitive to the model and the model parameters. In some models the process is inefficient; in some models it is not operative at all. Even if preheating is relatively efficient, it is unlikely to remove all of the energy density of in the inflaton field. It is likely that the slow decay of the inflaton field is necessary to extract the remaining inflaton field energy.

The simplest way to envision this process is if the comoving energy density in the zero mode of the inflaton decays into normal particles, which then scatter and thermalize to form a thermal background. It is usually assumed that the decay width of this process is the same as the decay width of a free inflaton field.

There are two reasons to suspect that the inflaton decay width might be small. The requisite flatness of the inflaton potential suggests a weak coupling of the inflaton field to other fields since the potential is renormalized by the inflaton coupling to other fields . However, this restriction may be evaded in supersymmetric theories where the nonrenormalization theorem ensures a cancelation between fields and their superpartners. A second reason to suspect weak coupling is that in local supersymmetric theories gravitinos are produced during reheating. Unless reheating is delayed, gravitinos will be overproduced, leading to a large undesired entropy production when they decay after big-bang nucleosynthesis .

Of particular interest is a quantity known as the reheat temperature, denoted as $`T_{RH}`$. The reheat temperature is calculated by assuming an instantaneous conversion of the energy density in the inflaton field into radiation when the decay width of the inflaton energy, $`\mathrm{\Gamma }_\varphi `$, is equal to $`H`$, the expansion rate of the universe.

The reheat temperature is calculated quite easily . After inflation the inflaton field executes coherent oscillations about the minimum of the potential. Averaged over several oscillations, the coherent oscillation energy density redshifts as matter: $`\rho _\varphi a^3`$, where $`a`$ is the Robertson–Walker scale factor. If we denote as $`\rho _I`$ and $`a_I`$ the total inflaton energy density and the scale factor at the initiation of coherent oscillations, then the Hubble expansion rate as a function of $`a`$ is ($`M_{Pl}`$ is the Planck mass)

$$H(a)=\sqrt{\frac{8\pi }{3}\frac{\rho _I}{M_{Pl}^2}\left(\frac{a_I}{a}\right)^3}.$$
(1)   
Equating $`H(a)`$ and $`\mathrm{\Gamma }_\varphi `$ leads to an expression for $`a_I/a`$. Now if we assume that all available coherent energy density is instantaneously converted into radiation at this value of $`a_I/a`$, we can define the reheat temperature by setting the coherent energy density, $`\rho _\varphi =\rho _I(a_I/a)^3`$, equal to the radiation energy density, $`\rho _R=(\pi ^2/30)g_{}T_{RH}^4`$, where $`g_{}`$ is the effective number of relativistic degrees of freedom at temperature $`T_{RH}`$. The result is

$$T_{RH}=\left(\frac{90}{8\pi ^3g_{}}\right)^{1/4}\sqrt{\mathrm{\Gamma }_\varphi M_{Pl}}=0.2\left(\frac{200}{g_{}}\right)^{1/4}\sqrt{\mathrm{\Gamma }_\varphi M_{Pl}}.$$
(2)   
The limit from gravitino overproduction is $`T_{RH}\stackrel{<}{_{}}10^9`$ to $`10^{10}`$ GeV.

The reheat temperature is best regarded as the temperature below which the universe expands as a radiation-dominated universe, with the scale factor decreasing as $`g_{}^{1/3}T^1`$. In this regard it has a limited meaning . For instance, $`T_{RH}`$ should not be used as the maximum temperature obtained by the universe during reheating. The maximum temperature is, in fact, much larger than $`T_{RH}`$. One implication of this is that it is incorrect to assume that the maximum abundance of a massive particle species produced after inflation is suppressed by a factor of $`\mathrm{exp}(M/T_{RH})`$.

In this paper we illustrate this effect by calculating the abundance of a massive particle species produced in reheating. We show that particles of mass much greater than the eventual “reheat” temperature $`T_{RH}`$ may be created by the thermalized decay products of the inflaton. As an example, we demonstrate that a stable particle species $`X`$ of mass $`M_X`$ would be produced in the reheating process in sufficient abundance that its contribution to closure density today is approximately $`M_X^2\sigma |v|(g_{}/200)^{3/2}(2000T_{RH}/M_X)^7`$, where $`g_{}`$ is the number of effective degrees of freedom of the radiation energy density and $`\sigma |v|`$ is the thermal average of the $`X`$ annihilation cross section times the Møller flux factor. Thus, particles of mass as large as $`2000`$ times the reheat temperature may be produced in interesting abundance.

Other applications of the effect include production of massive Higgs bosons which could decay and produce the baryon asymmetry, or massive particles that could decay and produce high-energy cosmic rays.

In the next section we develop a system of Boltzmann equations describing the evolution of the energy densities of the inflaton field, radiation, and a massive particle species. In Section III we find analytic approximations to the system and estimate the contribution to the present critical density from a stable, massive particle produced in reheating. Section IV contains some numerical results which illustrate several generic features. We conclude in Section V by discussing some applications of our results. The assumption of local thermodynamic equilibrium for the light degrees of freedom is addressed in an appendix.

II. BOLTZMANN EQUATIONS DESCRIBING REHEATING

Let us consider a model universe with three components: inflaton field energy, $`\rho _\varphi `$, radiation energy density, $`\rho _R`$, and the energy density of a nonrelativistic particle species, $`\rho _X`$. We will assume that the decay rate of the inflaton field energy density is $`\mathrm{\Gamma }_\varphi `$, with a branching fraction into $`X\overline{X}`$ of $`B_X`$, and a branching fraction $`1B_X`$ into light degrees of freedom, generically referred to as radiation. We will denote the decay width of the $`X`$ as $`\mathrm{\Gamma }_X`$. We will also assume that the light degrees of freedom are in local thermodynamic equilibrium. This is by no means guaranteed, and we will return to the question in the appendix.

With the above assumptions, the Boltzmann equations describing the redshift and interchange in the energy density among the different components is

$`\dot{\rho }_\varphi +3H\rho _\varphi +\mathrm{\Gamma }_\varphi \rho _\varphi =0`$   
$`\dot{\rho }_R+4H\rho _R(1B_X)\mathrm{\Gamma }_\varphi \rho _\varphi {\displaystyle \frac{\sigma |v|}{m_X}}\left[\rho _X^2\left(\rho _X^{EQ}\right)^2\right]\mathrm{\Gamma }_X\left(\rho _X\rho _X^{EQ}\right)=0`$   
$`\dot{\rho }_X+3H\rho _XB_X\mathrm{\Gamma }_\varphi \rho _\varphi +{\displaystyle \frac{\sigma |v|}{m_X}}\left[\rho _X^2\left(\rho _X^{EQ}\right)^2\right]+\mathrm{\Gamma }_X\left(\rho _X\rho _X^{EQ}\right)=0,`$ (3)   
where dot denotes time derivative As already mentioned, $`\sigma |v|`$ is the thermal average of the $`X`$ annihilation cross section times the Møller flux factor. The equilibrium energy density for the $`X`$ particles, $`\rho _X^{EQ}`$, is determined by the radiation temperature, $`T`$.

It is useful to introduce two dimensionless constants, $`\alpha _\varphi `$ and $`\alpha _X`$, defined in terms of $`\mathrm{\Gamma }_\varphi `$ and $`\sigma |v|`$ as

$$\mathrm{\Gamma }_\varphi =\alpha _\varphi M_\varphi \sigma |v|=\alpha _XM_X^2.$$
(4)   
For a reheat temperature much smaller than $`M_\varphi `$, $`\mathrm{\Gamma }_\varphi `$ must be small. From Eq. (2), the reheat temperature in terms of $`\alpha _X`$ and $`M_X`$ is $`T_{RH}\alpha _\varphi ^{1/2}\sqrt{M_\varphi M_{Pl}}`$. For $`M_\varphi =10^{13}`$GeV, $`\alpha _\varphi `$ must be smaller than of order $`10^{13}`$. On the other hand, $`\alpha _X`$ may be as large as of order unity, or it may be small also.

In what follows we will make the simplifying assumption that $`B_X=0`$. Since we are interested in a stable particle relic, we will assume that $`\mathrm{\Gamma }_X=0`$. It is also convenient to work with dimensionless quantities that can absorb the effect of expansion of the universe. This may be accomplished with the definitions

$$\mathrm{\Phi }\rho _\varphi M_\varphi ^1a^3;R\rho _Ra^4;X\rho _XM_X^1a^3.$$
(5)   
It is also convenient to use the scale factor, rather than time, for the independent variable, so we define a variable $`x=aM_\varphi `$. With this choice the system of equations can be written as (prime denotes $`d/dx`$)

$`\mathrm{\Phi }^{}`$ $`=`$ $`c_1{\displaystyle \frac{x}{\sqrt{\mathrm{\Phi }x+R}}}\mathrm{\Phi }`$   
$`R^{}`$ $`=`$ $`c_1{\displaystyle \frac{x^2}{\sqrt{\mathrm{\Phi }x+R}}}\mathrm{\Phi }+c_2{\displaystyle \frac{x^1}{\sqrt{\mathrm{\Phi }x+R}}}\left(X^2X_{EQ}^2\right)`$   
$`X^{}`$ $`=`$ $`c_3{\displaystyle \frac{x^2}{\sqrt{\mathrm{\Phi }x+R}}}\left(X^2X_{EQ}^2\right).`$ (6)   
The constants $`c_1`$, $`c_2`$, and $`c_3`$ are given by

$$c_1=\sqrt{\frac{3}{8\pi }}\frac{M_{Pl}}{M_\varphi }\alpha _\varphi c_2=c_1\frac{M_\varphi }{M_X}\frac{\alpha _X}{\alpha _\varphi }c_3=c_2\frac{M_\varphi }{M_X}.$$
(7)   
$`X_{EQ}`$ is the equilibrium value of $`X`$, given in terms of the temperature $`T`$ as (assuming a single degree of freedom for the $`X`$ species)

$$X_{EQ}=\frac{M_X^3}{M_\varphi ^3}\left(\frac{1}{2\pi }\right)^{3/2}x^3\left(\frac{T}{M_X}\right)^{3/2}\mathrm{exp}(M_X/T).$$
(8)   
The temperature depends upon $`R`$ and $`g_{}`$, the effective number of degrees of freedom in the radiation:

$$\frac{T}{M_X}=\left(\frac{30}{g_{}\pi ^2}\right)^{1/4}\frac{M_\varphi }{M_X}\frac{R^{1/4}}{x}.$$
(9)   
It is straightforward to solve the system of equations in Eq. (S0.Ex3) with initial conditions at $`x=x_I`$ of $`R(x_I)=X(x_I)=0`$ and $`\mathrm{\Phi }(x_I)=\mathrm{\Phi }_I`$. It is convenient to express $`\rho _\varphi (x=x_I)`$ in terms of the expansion rate at $`x_I`$, which leads to

$$\mathrm{\Phi }_I=\frac{3}{8\pi }\frac{M_{Pl}^2}{M_\varphi ^2}\frac{H_I^2}{M_\varphi ^2}x_I^3.$$
(10)   
The numerical value of $`x_I`$ is irrelevant.

Before numerically solving the system of equations, it is useful to consider the early-time solution for $`R`$. Here, by early time, we mean $`H\mathrm{\Gamma }_\varphi `$, i.e., before a significant fraction of the comoving coherent energy density is converted to radiation. At early times $`\mathrm{\Phi }\mathrm{\Phi }_I`$, and $`RX0`$, so the equation for $`R^{}`$ becomes $`R^{}=c_1x^{3/2}\mathrm{\Phi }_I^{1/2}`$. Thus, the early time solution for $`R`$ is simple to obtain:

$$R\frac{2}{5}c_1\left(x^{5/2}x_I^{5/2}\right)\mathrm{\Phi }_I^{1/2}(H\mathrm{\Gamma }_\varphi ).$$
(11)   
Now we may express $`T`$ in terms of $`R`$ to yield the early-time solution for $`T`$:

$$\frac{T}{M_\varphi }\left(\frac{12}{\pi ^2g_{}}\right)^{1/4}c_1^{1/4}\left(\frac{\mathrm{\Phi }_I}{x_I^3}\right)^{1/8}\left[\left(\frac{x}{x_I}\right)^{3/2}\left(\frac{x}{x_I}\right)^4\right]^{1/4}(H\mathrm{\Gamma }_\varphi ).$$
(12)   
Thus, $`T`$ has a maximum value of

$`{\displaystyle \frac{T_{MAX}}{M_\varphi }}`$ $`=`$ $`0.77\left({\displaystyle \frac{12}{\pi ^2g_{}}}\right)^{1/4}c_1^{1/4}\left({\displaystyle \frac{\mathrm{\Phi }_I}{x_I^3}}\right)^{1/8}`$ (13)   
$`=`$ $`0.77\alpha _\varphi ^{1/4}\left({\displaystyle \frac{9}{2\pi ^3g_{}}}\right)^{1/4}\left({\displaystyle \frac{M_{Pl}^2H_I}{M_\varphi ^3}}\right)^{1/4},`$   
which is obtained at $`x/x_I=(8/3)^{2/5}=1.48`$. It is also possible to express $`\alpha _\varphi `$ in terms of $`T_{RH}`$ and obtain

$$\frac{T_{MAX}}{T_{RH}}=0.77\left(\frac{9}{5\pi ^3g_{}}\right)^{1/8}\left(\frac{H_IM_{Pl}}{T_{RH}^2}\right)^{1/4}.$$
(14)   
For an illustration, in the simplest model of chaotic inflation $`H_IM_\varphi `$ with $`M_\varphi 10^{13}`$GeV, which leads to $`T_{MAX}/T_{RH}10^3(200/g_{})^{1/8}`$ for $`T_{RH}=10^9`$GeV.

We can see from Eq. (11) that for $`x/x_I>1`$, in the early-time regime $`T`$ scales as $`a^{3/8}`$, which implies that entropy is created in the early-time regime . So if one is producing a massive particle during reheating it is necessary to take into account the fact that the maximum temperature is greater than $`T_{RH}`$, and that during the early-time evolution, $`Ta^{3/8}`$.

III. PRODUCTION OF A MASSIVE, STABLE PARTICLE SPECIES

A. Freeze out of the comoving $`X`$ energy density

In this section we develop the equation for the $`X_F`$, the final value of $`X`$, which can be found from the early-time behavior.

At early times $`\mathrm{\Phi }\mathrm{\Phi }_I`$ and $`R0`$. We will here also assume that $`XX_{EQ}`$. Numerical results confirm the validity of this approximation and show that the massive particles are never in chemical equilibrium (although presumably they are in kinetic equilibrium). The early-time equation for the development of the $`X`$ energy density is

$$X^{}=c_3\mathrm{\Phi }_I^{1/2}x^{5/2}X_{EQ}^2.$$
(15)   
$`X_{EQ}`$ is given in terms of $`M_X`$ and the temperature, which may be found from the early-time solution for $`R`$.

We can integrate Eq. (15) by approximating it as a Gaussian integral. First we rearrange Eq. (15) by making appropriate redefinitions. Define the quantities $`y`$ and $`\nu `$ by $`yX/x_I^3`$ and $`\nu (x/x_I)^{3/16}`$. Now, using Eq. (12), we can rewrite Eq. (15) as

$$y(\nu )=Q_1^\nu 𝑑\nu ^{}\mathrm{exp}[H(\nu ^{})]$$
(16)   
where we have defined

$`Q`$ $``$ $`\alpha _\varphi ^{3/4}\alpha _X\left({\displaystyle \frac{3^{1/2}2^{1/4}}{\pi ^{21/4}g_{}^{3/4}}}\right){\displaystyle \frac{M_{Pl}^{3/2}M_X}{H_I^{1/4}M_\varphi ^{9/4}}}`$   
$`H(\nu )`$ $``$ $`\lambda \left(\nu ^8\nu ^{64/3}\right)^{1/4}{\displaystyle \frac{3}{4}}\mathrm{ln}\left(\nu ^8\nu ^{64/3}\right)23\mathrm{ln}\nu `$   
$`\lambda `$ $``$ $`{\displaystyle \frac{2^{5/4}\pi ^{3/4}g_{}^{1/4}}{3^{1/2}}}{\displaystyle \frac{\alpha _\varphi ^{1/4}M_X}{M_\varphi ^{1/4}H_I^{1/4}M_{Pl}^{1/2}}}.`$ (17)   
To proceed with the Gaussian integral approximation, we assume $`\nu _0^8\nu _0^{64/3}`$ where $`\nu _0`$ is the solution to $`H^{}(\nu _0)=0`$. Then, we can easily solve $`H^{}(\nu _0)=0`$, finding $`\nu _0=\sqrt{17/2\lambda }`$, which is the point about which we Taylor expand $`H(\nu )`$ to quadratic order. Since the integrand falls to 0 rapidly away from $`\nu =\nu _0`$, and since we desire the freeze out value for $`y`$, the limits of the integrand can be taken to $`\pm \mathrm{}`$. We thereby find

$$\frac{X_F}{x_I^3}y_{\mathrm{}}\frac{3^5\alpha _\varphi ^3\alpha _X}{8\pi ^{23/2}g_{}^3}\left(\frac{H_I^2M_{Pl}^6}{M_X^8}\right)\left(\frac{\sqrt{17}}{2}\right)^{17}\mathrm{exp}(17/2).$$
(18)   
Using Eq. (2), we rewrite this in a more transparent form

$$\frac{X_F}{x_I^3}\frac{4.21\times 10^6}{(g_{}/200)^{3/2}}\frac{\alpha _XH_I^2M_{Pl}^3T_{RH}^6}{M_\varphi ^3M_X^8}.$$
(19)   
Note that this approximation should be valid as long as $`\nu _0^8\nu _0^{64/3}`$ is satisfied. If the condition is not satisfied, the suppression will be exponential in $`M_X/T`$.

B. $`\mathrm{\Omega }_Xh^2`$ in terms of $`X_F`$

After freeze out of the comoving energy density of the stable particle, $`X`$ remains constant, so $`\rho _X(x>x_F)=X_Fx_I^3M_XM_\varphi ^3(x_I/x)^3`$. For delayed reheating ($`\mathrm{\Gamma }_\varphi H_I`$) freeze out will be well before reheating. After reheating, $`\rho _X(x>x_{RH})=\rho _X(x_{RH})(x_{RH}/x)^3`$. The comoving entropy density is constant after reheating, so the radiation energy density scales as $`\rho _R(x>x_{RH})=\rho _R(x_{RH})[g_{}(T_{RH})/g_{}(T)]^{1/3}(x_{RH}/x)^4`$. Using these facts, we can express the present contribution of the massive particle species to the critical density in terms of the ratio of the energy densities at freeze out:

$$\frac{\mathrm{\Omega }_Xh^2}{\mathrm{\Omega }_Rh^2}=\frac{\rho _X(T_{RH})}{\rho _R(T_{RH})}\left(\frac{g_{}(\text{today})}{g_{}(T_{RH})}\right)^{1/3}\frac{x_0}{x_{RH}}=\frac{\rho _X(T_{RH})}{\rho _R(T_{RH})}\frac{T_{RH}}{T_0}(\text{for}x>x_{RH}),$$
(20)   
where $`x_0`$ is the present value of $`x`$ and $`T_0=2.37\times 10^{13}`$GeV is the present temperature.<sup>2</sup><sup>2</sup>2In this subsection we make the heretofore criticized approximation that the inflaton energy density scales like pressureless matter until it dumps all of its energy into radiation at the instant of “reheating.” In this instance, however, it is an appropriate approximation, as borne out by analytic approximations and the numerical calculations presented in the next section. Today, $`\mathrm{\Omega }_Rh^2=4.3\times 10^5`$, and the contribution to $`\mathrm{\Omega }h^2`$ from the massive particle is

$$\mathrm{\Omega }_Xh^2=1.5\times 10^{18}\left(\frac{T_{RH}}{10^9\mathrm{GeV}}\right)\frac{X_F}{x_I^3}\frac{M_XM_\varphi ^3}{H_I^2M_{Pl}^2}.$$
(21)   
Using the expression for $`X_F/x_I^3`$ from the previous section, we arrive at the final result

$$\mathrm{\Omega }_Xh^2=M_X^2\sigma |v|\left(\frac{g_{}}{200}\right)^{3/2}\left(\frac{2000T_{RH}}{M_X}\right)^7.$$
(22)   
IV. NUMERICAL RESULTS

An example of a numerical evaluation of the complete system in Eq. (S0.Ex3) is shown in Fig. 1. The model parameters chosen were $`M_\varphi =10^{13}`$GeV, $`\alpha _\varphi =2\times 10^{13}`$, $`M_X=1.15\times 10^{12}`$GeV, $`\alpha _X=10^2`$, and $`g_{}=200`$. The expansion rate at the beginning of the coherent oscillation period was chosen to be $`H_I=M_\varphi `$. These parameters result in $`T_{RH}=10^9`$GeV and $`\mathrm{\Omega }_Xh^2=0.3`$.

Figure 1 serves to illustrate several aspects of the problem. Just as expected, the comoving energy density of $`\varphi `$ (i.e., $`a^3\rho _\varphi `$) remains roughly constant until $`\mathrm{\Gamma }_\varphi H`$, which for the chosen model parameters occurs around $`a/a_I5\times 10^8`$. But of course, that does not mean that the temperature is zero. Notice that the temperature peaks well before “reheating.” The maximum temperature, $`T_{MAX}=10^{12}`$GeV, is reached at $`a/a_I`$ slightly larger than unity (in fact at $`a/a_I=1.48`$ as expected), while the reheat temperature, $`T_{RH}=10^9`$GeV, occurs much later, around $`a/a_I10^8`$. Note that $`T_{MAX}10^3T_{RH}`$ in agreement with Eq. (14).

From the numerical results we can justify one of the assumptions in deriving the analytical approximations. From the figure it is clear that $`XX_{EQ}`$ at the epoch of freeze out of the comoving $`X`$ number density, which occurs around $`a/a_I10^2`$. The rapid rise of the ratio after freeze out is simply a reflection of the fact that $`X`$ is constant while $`X_{EQ}`$ decreases exponentially. The relevance of the ratio is the justification of the neglect of $`X_{EQ}`$ term in Eq. (15).

A close examination of the behavior of $`T`$ shows that after the sharp initial rise of the temperature, the temperature decreases as $`a^{3/8}`$ \[as follows from Eq. (12)\] until $`H\mathrm{\Gamma }_\varphi `$, and thereafter $`Ta^1`$ as expected for the radiation-dominated era.

For the choices of $`M_\varphi `$, $`\alpha _\varphi `$, $`g_{}`$, and $`\alpha _X`$ used for the model illustrated in Fig. 1, $`\mathrm{\Omega }_Xh^2=0.3`$ for $`M_X=1.15\times 10^{12}`$GeV, in excellent agreement with the mass predicted by using Eq. (22).

V. CONCLUSIONS

Let us now analyze the implications of our findings for the GUT baryogenesis scenario, where the baryon asymmetry is produced by the baryon (or lepton) number violating decays of superheavy bosons . At the end of inflation the Universe does not contain any matter and, even more important, it is perfectly baryon symmetric—there is no dominance of matter over antimatter. This means that GUT baryogenesis may be operative only if the supermassive GUT bosons are regenerated during the stage of thermalization of the decay products of the inflaton field<sup>3</sup><sup>3</sup>3In the case in which reheating is anticipated by a period of preheating, superheavy bosons may be produced by the phenomenon of parametric resonance .. A naive estimate would lead to conclude that that the maximum number density of a massive particle species $`X`$ produced after inflation is suppressed by a factor of $`(M_X/T_{RH})^{3/2}\mathrm{exp}(M_X/T_{RH})`$ with respect to the photon number density. For such a reason, it is commonly believed that GUT baryogenesis is incompatible with models of inflation where the reheating temperature is much smaller than the GUT scale and, in general, than the mass of the $`X`$ particles, $`T_{RH}M_X`$. In fact, we have seen that the reheat temperature has a limited meaning and should not be used as the maximum temperature obtained by the universe during reheating. The maximum temperature, Eq. (14), is much larger than $`T_{RH}`$, and particles of mass much greater than the eventual reheating temperature $`T_{RH}`$ may be created by the thermalized decay products of the inflaton without any exponential suppression factor. Indeed, the number density $`n_X`$ of particles $`X`$ after freeze out and reheating may be easily inferred from Eqs. (5) and (19), and reads

$$\frac{n_X}{n_\gamma }3\times 10^4\left(\frac{100}{g_{}}\right)^{3/2}\left(\frac{T_{RH}}{M_X}\right)^7\left(\frac{M_{Pl}}{M_X}\right).$$
(23)   
In theories that invoke supersymmetry to preserve the flatness of the inflaton potential, the slow decay rate of the gravitinos, the superpartners of the gravitons, is a source of the cosmological problems because the decay products of the gravitino will destroy the <sup>4</sup>He and D nuclei by photodissociation, and in the process destroy the successful nucleosynthesis predictions. The most stringent bound comes from the resulting overproduction of D $`+`$ <sup>3</sup>He, which would require that the gravitino abundance is smaller than about $`10^{10}`$ relative to the entropy density at the time of reheating after inflation. This translates into an upper bound on the reheating temperature after inflation, $`T_{RH}/M_{Pl}\stackrel{<}{_{}}10^9`$ .

It is easy to check that for such small values of $`T_{RH}`$ the ratio in Eq. (23) is always much larger than the equilibrium value, $`n_X^{EQ}/n_\gamma =(M_X/2T_{RH})^{3/2}(\pi ^{1/2}/\xi (3))\mathrm{exp}(M_X/T_{RH})`$. This result is crucial for the out-of-equilibrium decay scenarios of baryogenesis. For instance, in theories where $`BL`$ is a spontaneously broken local symmetry, as suggested by $`SO(10)`$ unification, the cosmological baryon asymmetry can be generated by the out-of-equilibrium decay of the lightest heavy Majorana right-handed neutrino $`N_1^c`$ , whose typical mass is about $`10^{10}`$ GeV. For reheat temperatures of the order of $`10^9`$ GeV, the number density of the right-handed neutrino is about $`3\times 10^2n_\gamma `$ and one can estimate the final baryon number to be of the order of $`B(n_{N_1^c}/n_\gamma )(ϵ/g_{})10^4ϵ`$, where $`ϵ`$ is the coefficient containing one-loop suppression factor and $`CP`$ violating phases. The observed value of the baryon asymmetry, $`B10^{10}`$, is then obtained without any fine tuning of parameters.

Our findings have also important implications for the conjecture that ultra-high cosmic rays, above the Greisen-Zatsepin-Kuzmin cut-off of the cosmic ray spectrum, may be produced in decays of superheavy long-living particles . In order to produce cosmic rays of energies larger than about $`10^{13}`$ GeV, the mass of the $`X`$-particles must be very large, $`M_X\stackrel{>}{_{}}10^{13}`$ GeV and their lifetime $`\tau _X`$ cannot be much smaller than the age of the Universe, $`\tau _X\stackrel{>}{_{}}10^{10}`$ yr. With the smallest value of the lifetime, the observed flux of ultra-high energy cosmic rays will be reproduced with a rather low density of $`X`$-particles, $`\mathrm{\Omega }_X10^{12}`$. It has been suggested that $`X`$-particles can be produced in the right amount by usual collisions and decay processes taking place during the reheating stage after inflation, if the reheat temperature never exceeded $`M_X`$ . Again, assuming naively that that the maximum number density of a massive particle species $`X`$ produced after inflation is suppressed by a factor of $`(M_X/T_{RH})^{3/2}\mathrm{exp}(M_X/T_{RH})`$ with respect to the photon number density, one concludes that the reheat temperature $`T_{RH}`$ should be in the range $`10^{11}`$ to $`10^{15}`$GeV . This is a rather high value and leads to the gravitino problem in generic supersymmetric models. This is one reason alternative production mechanisms of these superheavy $`X`$-particles have been proposed . However, our analysis show that the situation is much more promising. Making use of Eq. (22), the right amount of $`X`$-particles to explain the observed ultra-high energy cosmic rays is produced for

$$\left(\frac{T_{RH}}{10^{10}\mathrm{GeV}}\right)\left(\frac{g_{}}{200}\right)^{3/14}\left(\frac{M_X}{10^{15}\mathrm{GeV}}\right),$$
(24)   
where we have assumed $`\sigma |v|M_X^2`$. Therefore, we conclude that particles as massive as $`10^{15}`$ GeV may be generated during the reheating stage in abundances large enough to explain the ultra-high energy cosmic rays even if the reheat temperature satisfies the gravitino bound.

ACKNOWLEDGEMENTS

DJHC and EWK were supported by the DOE and NASA under Grant NAG5-7092.

APPENDIX A: THERMALIZATION OF LIGHT DEGREES OF FREEDOM

The form of the Boltzmann equations we use, e.g., Eq. (S0.Ex1), assumes that the light particle decay products of the inflaton field are in local thermodynamic equilibrium (LTE). In this appendix we discuss this assumption, and the implications if it is not valid.

Before discussing the validity of the assumption, it is useful to recall why the assumption was made. In the derivation of Eq. (S0.Ex1), one starts with an equation for the rate of change of the $`X`$ number density due to the process $`\gamma \gamma XX`$ with four-momentum conservation $`p_\gamma +p_\gamma ^{}=p_X+p_X^{}`$:

$`\dot{n}_X`$ $`=`$ $`{\displaystyle \frac{d^3p_\gamma }{2E_\gamma }\frac{d^3p_\gamma ^{}}{2E_\gamma ^{}}\frac{d^3p_X}{2E_X}\frac{d^3p_X^{}}{2E_X^{}}\left(2\pi \right)^8\delta ^4\left(p_\gamma +p_\gamma ^{}p_Xp_X^{}\right)}`$ (A1)   
$`\times f_\gamma (p_\gamma )f_\gamma (p_\gamma ^{})\left|\right|_{\gamma \gamma XX}^2+\mathrm{}.`$   
Here $`f_i(p)`$ is the phase-space density of particle species $`i`$ with momentum $`p_i`$, and $`\left|\right|_{\gamma \gamma XX}^2`$ is the square of the matrix element for the process $`\gamma \gamma XX`$. With the assumption that the light particles are in LTE, the product of the light particle phase-space densities is $`f_\gamma (p_\gamma )f_\gamma (p_\gamma ^{})=\mathrm{exp}(E_\gamma /T)\mathrm{exp}(E_\gamma ^{}/T)`$. This last product is, of course, simply $`f_X^{EQ}(p_X)f_X^{EQ}(p_X^{})`$, which, after some rearrangement of Eq. (A1), leads to a term for the creation of $`X`$’s proportional to $`\left(n_X^{EQ}\right)^2`$:

$$\dot{n}_X=\sigma |v|\left(n_X^{EQ}\right)^2+\mathrm{}.$$
(A2)   
(for complete details, see ).

The factor $`\left(n_X^{EQ}\right)^2`$ in Eq. (A2) is present because not every light-particle collision has sufficient center-of-mass energy to create an $`X`$ pair. If LTE is established with temperature $`T<M_X`$, the factor $`\mathrm{exp}(2M_X/T)`$ in $`\left(n_X^{EQ}\right)^2`$ represents the fraction of the collisions with center-of-mass energy above threshold, i.e., with $`\sqrt{s}>2M_X`$.

A simple indication of whether thermalization occurs on a timescale shorter than the timescale for $`X`$ production is the ratio of the cross section for the thermalization reactions to the cross section for $`X`$ production. If the ratio is larger than unity, then thermalization of the light degrees of freedom is a good assumption.

The process of $`X`$ production involves a “hard” process, and the cross section will be $`\alpha _X/M_X^2`$, where $`\alpha _X`$ was defined in Eq. (4). In order to produce an equilibrium distribution from the original decay distribution it is necessary to change the number of particles. Therefore, the relevant cross section is the one for processes like $`\gamma \gamma \gamma \gamma \gamma `$.<sup>4</sup><sup>4</sup>4Recall that $`\gamma `$ represents a light particle, not just a photon, so $`\gamma `$ may carry electric charge, color charge, etc. Although the thermalization reaction is higher order in perturbation theory, it is a “softer” process, and radiation of a soft photon has a large cross section.

Without knowing the details of the interactions of the decay products, it is impossible to say with certainty how complete thermalization will be. But if the inflaton decay products have usual gauge interactions, the thermalization cross section will be larger than the $`X`$ production cross section, and thermalization of the inflaton decay products is likely.

Now let’s explore the consequences if LTE of the light degrees of freedom is not obtained. If the light particles are not in LTE, then the factor $`n_X^{EQ}`$ in Eq. (A2) could simply be replaced by the more general factor $`n_\gamma (E>M_X)`$.<sup>5</sup><sup>5</sup>5Of course in this case $`\sigma |v|`$ would not be a thermal average, but an average over the actual phase-space density.  Now let’s make the extreme assumption that the light degrees of freedom never interact before $`X`$ production, and that they have the original (redshifted) momentum with which they were created in inflaton decay. Assume this original momentum is $`M_\varphi /\eta `$. If $`M_\varphi /\eta `$ is greater than $`M_X`$, then the $`X`$ production rate actually will be larger than the equilibrium rate since $`n_\gamma (E>M_X)>n_X^{EQ}`$, while if $`M_\varphi /\eta `$ is less than $`M_X`$, then the $`X`$ production rate will be zero! The most reasonable assumption is that even if there is a large multiplicity in $`\varphi `$ decay, a fair fraction of $`M_\varphi `$ is carried by a few leading particles. So the effective value of $`\eta `$ is probably not too large.

The above analysis leads us to the conclusion that thermalization of the light degrees of freedom is likely unless the inflaton decay products themselves are very weakly coupled to everything. Even if thermalization does not occur, production of massive particles during reheating is not much different than our simple model suggests.