Advances in High Energy Physics

Volume 2017 (2017), Article ID 7650238, 24 pages

https://doi.org/10.1155/2017/7650238

## Irreversible Thermodynamic Description of Dark Matter and Radiation Creation during Inflationary Reheating

^{1}School of Physics, Sun Yat-Sen University, Guangzhou 510275, China^{2}Yat Sen School, Sun Yat-Sen University, Guangzhou 510275, China^{3}Department of Physics, Babes-Bolyai University, Kogalniceanu Street, 400084 Cluj-Napoca, Romania^{4}Department of Mathematics, University College London, Gower Street, London WC1E 6BT, UK^{5}State Key Laboratory of Optoelectronic Material and Technology and Guangdong Province Key Laboratory of Display Material and Technology, Guangzhou, China

Correspondence should be addressed to Tiberiu Harko; ku.ca.lcu@okrah.t

Received 10 July 2017; Accepted 23 August 2017; Published 18 October 2017

Academic Editor: Orlando Luongo

Copyright © 2017 Juntong Su et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The publication of this article was funded by SCOAP^{3}.

#### Abstract

We investigate matter creation processes during the reheating period of the early Universe, by using the thermodynamic of open systems. The Universe is assumed to consist of the inflationary scalar field, which, through its decay, generates relativistic matter and pressureless dark matter. The inflationary scalar field transfers its energy to the newly created matter particles, with the field energy decreasing to near zero. The equations governing the irreversible matter creation are obtained by combining the thermodynamics description of the matter creation and the gravitational field equations. The role of the different inflationary scalar field potentials is analyzed by using analytical and numerical methods. The values of the energy densities of relativistic matter and dark matter reach their maximum when the Universe is reheated up to the reheating temperature, which is obtained as a function of the scalar field decay width, the scalar field particle mass, and the cosmological parameters. Particle production leads to the acceleration of the Universe during the reheating phase, with the deceleration parameter showing complex dynamics. Once the energy density of the scalar field becomes negligible with respect to the matter densities, the expansion of the Universe decelerates, and inflation has a graceful exit.

#### 1. Introduction

The inflationary paradigm is one of the corner stones of present day cosmology. Inflation was first proposed in [1] to solve the spatial flatness, the horizon, and the monopole problems of the early Universe. The main theoretical ingredients in inflation are a scalar field with self-interaction potential , its energy density , and pressure , respectively [2–5]. In the initial approach the inflation potential arrives at a local minimum at through supercooling from a phase transition, and then the Universe expands exponentially. However, in this scenario, called “old inflation,” there is no graceful exit to the inflationary stage. To solve this problem the “new inflation” model was proposed in [6–9]. In this model the inflationary scalar field is initially at its local maximum . The potential is required to be very flat near the minimum at , so that the field rolls slowly down, without changing much before the Universe expands exponentially. The “flatness” requirement is not easy to obtain, and another problem of this scenario is that the inflation has to begin very late and last a long time, or it may never happen in some cases [10]. By extending the new inflation model, the chaotic inflation theory was proposed [11], which solved the problems above. The model does not require a thermal equilibrium state of the early Universe, and various forms of the scalar field potentials are acceptable, so that a large number potential forms can be used in this scenario. Later a subtler theory of hybrid inflation with two scalar fields in the model was proposed [12]. The potential takes the form of , with , , , and constants. One scalar field accounts for the exponential growth of the Universe, and the other one is responsible for the graceful exit of inflation.

Presently, very precise observations of the Cosmic Microwave Background (CMB) radiation give us the opportunity to test the fundamental predictions of inflation on primordial fluctuations, such as scale independence and Gaussianity [13–17]. One can calculate cosmological parameters from CMB fluctuations and use them as constraints of the inflationary models. The slow-roll parameters and can be calculated directly from the potentials of the inflationary scalar field, so that these parameters can be observed, and the corresponding inflationary models can be tested [18]. The slow-roll parameters and the inflationary parameters like scalar spectral index and tensor-to-scalar ratio of several forms of the inflationary potentials were investigated in [19]. Further investigations of the parameters and can be found in [20–22].

The prediction of the statistical isotropy, claiming that inflation removes the classical anisotropy of the Universe, has been a central prediction of inflationary scenarios, also strongly supported by the cosmic no-hair conjecture. However, recently several observations of the large scale structure of the Universe have raised the possibility that the principles of homogeneity and isotropy may not be valid at all scales and that the presence of an intrinsic large scale anisotropy in the Universe cannot be ruled out a priori [23]. Although generally the recent Planck observations tend to confirm the fundamental principles of the inflationary paradigm, the CMB data seem to show the existence of some tension between the models and the observations. For example, the combined Planck and WMAP polarization data show that the index of the power spectrum is [14, 15], at the pivot scale Mpc^{−1}, a result which definitely rules out the exact scale-invariance (), as predicted by some inflationary models [24, 25], at more than . Moreover, this case can also be ruled out based on some fundamental theoretical considerations, since it is connected with a de Sitter phase, and, more importantly, it does not lead to a consistent inflationary model, mainly due to the graceful exit problem. On the other hand the joint constraints on (tensor-to-scalar ratio) and can put strong restrictions on the inflationary scenarios. For example, inflationary models with power law potentials of the form cannot provide an acceptable number of -folds (of the order of 50–60) in the restricted space at a level [15]. In actuality, the theories of inflation have not yet reached a general consensus, because not only the overall inflationary picture, but also the details of the theoretical models are based on physics beyond the Standard Model of particle physics.

During inflation, the exponential growth of the Universe led to a homogeneous, isotropic, and empty Universe. Elementary particles are believed to be created in the period of reheating at the end of inflation, defrosting the Universe, and transferring energy from the inflationary scalar field to matter. Reheating was suggested along with the new inflationary scenario [6], and further developed in [26–28]. When the inflationary scalar field reaches its minimum after the expansion of the Universe, it oscillates around the minimum of the potential and decays into Standard Model particles. The newly created particles interact with each other and reach a thermal equilibrium state of temperature . However, quantum field theoretical investigations indicate that the reheating period is characterised by complicated nonequilibrium processes, the main characteristic of which are initial, violent particle production via parametric resonance and inflaton decay (preheating) with a highly nonequilibrium distribution of the produced particles, subsequently relaxing to an equilibrium state. Such explosive particle production is due to parametric amplification of quantum fluctuations for the nonbroken symmetry case (appropriate for chaotic inflation) or spinodal instabilities in the broken symmetry phase [29–32].

The reheating temperature is thought to be the maximum temperature of the radiation dominated Universe after reheating. However, it was argued that the reheating temperature is not necessarily the maximum temperature of the Universe after reheating [33]. In other words, reheating might be just one of the stages of the matter creation, and other substance with larger freeze-out temperature may be created after inflation through other process [34].

In the study of reheating, the phenomenological approach introduced in [26] has been widely investigated. This approach is based on the introduction of a suitable loss term in the scalar field equation, which also appears as a source term for the energy density of the newly created matter fluid. It is this source term that if chosen appropriately is responsible for the reheating process that follows adiabatic supercooling during the de Sitter phase. In this simple model the corresponding dynamics is considered within a two component model. The first component is a scalar field, with energy density and pressure and , respectively, and with energy-momentum tensor . The second component is represented by normal matter, with energy density and pressure and , and with energy-momentum tensor . In the above equations for simplicity we assume that both cosmological fluid components share the same four velocities , normalized as . Einstein’s field equations imply the relation , which in the case of a flat geometry reduces to where is the Hubble function and is the scale factor of the Universe. In order to describe the transition process between the scalar field and the matter component it is convenient to introduce phenomenologically a friction term , describing the decay of the scalar field component and acting as a source term for the matter fluid. Hence, (1) is assumed to decompose into two separate equations

In the framework of the braneworld models, in which our Universe is a 3-brane embedded in a five-dimensional bulk, the reheating era after the inflationary period was analyzed [35], by assuming the possibility of brane-bulk energy exchange. The inflaton field was assumed to decay into normal matter only, while the dark matter is injected into the brane from the bulk. The observational constraint of an approximately constant ratio of the dark and the baryonic matter requires that the dark matter must be nonrelativistic (cold). The model predicts a reheating temperature of the order of GeV, and brane tension of the order of GeV^{4}. The composition of the Universe thus obtained is consistent with the observational data. The problem of perturbative reheating and its effects on the evolution of the curvature perturbations in tachyonic inflationary models was investigated in [36]. It was found that reheating does not affect the amplitude of the curvature perturbations, and, after the transition, the relative nonadiabatic pressure perturbation decays extremely rapidly during the early stages of the radiation dominated epoch. These behaviors ensure that the amplitude of the curvature perturbations remains unaffected during reheating. It was pointed out in [37] that, among primordial magnetogenesis models, inflation is a prime candidate to explain the current existence of cosmological magnetic fields. Their energy density decreases as radiation during the decelerating eras of the Universe, and in particular during reheating. Avoiding magnetic field backreaction is always complementary to CMB and can give stronger limits on reheating for all high energy models of inflation. The reheating dynamics after the inflation induced by -corrected models were considered [38]. The inflationary and reheating dynamics was analyzed in the Einstein frame. Observational constraints on the model were also discussed.

The processes of particle production from the inflaton, their subsequent thermalization, and evolution of inflaton/plasma system by taking dissipation of the inflaton in a hot plasma into account were investigated in [39], and it was shown that the reheating temperature is significantly affected by these effects.

In modified gravity models the reheating dynamics after the inflation is significantly modified and affects the shape of the gravitational wave background spectrum [40]. In [41] it was shown that reheating considerations may provide additional constraints to some classes of inflationary models. The current Planck satellite measurements of the Cosmic Microwave Background anisotropies constrain the kinematic properties of the reheating era for most of the inflationary models [42]. This result is obtained by deriving the marginalized posterior distributions of the reheating parameter for about 200 models. It turns out that the precision of the current CMB data is such that estimating the observational performance of an inflationary model now requires incorporating information about its reheating history. It was pointed out in [43] that the inflaton coupling to the Standard Model particles is generated radiatively, even if absent at tree level. Hence, the dynamics of the Higgs field during inflation can change dramatically. The inflaton decay and reheating period after the end of inflation in the nonminimal derivative coupling gravity model with chaotic potential was investigated in [44]. This model provides an enhanced slow-roll inflation, caused by gravitationally enhanced friction. A scale-invariant model of quadratic gravity with a nonminimally coupled scalar field was studied in [45]. At the end of inflation, the Hubble parameter and the scalar field converge to a stable fixed point through damped oscillations that can reheat the Universe in various ways. For reviews on the inflationary reheating see [46] and [47], respectively.

If the importance of the nonequilibrium and dissipative aspects of particle creation during reheating has been already pointed out a long time ago [48], its irreversible and open character has not received a similar attention. Thermodynamical systems in which matter creation occurs belong to the class of open thermodynamical systems, in which the usual adiabatic conservation laws are modified by explicitly including irreversible matter creation [49]. The thermodynamics of open systems has been also applied first to cosmology in [49]. The explicit inclusion of the matter creation in the matter energy stress tensor in the Einstein field equations leads to a three-stage cosmological model, beginning with an instability of the vacuum. During the first stage the Universe is driven from an initial fluctuation of the vacuum to a de Sitter space, and particle creation occurs. The de Sitter space exists during the decay time of its constituents (second stage) and ends, after a phase transition, in the usual FRW phase. The phenomenological approach of [49] was further discussed and generalized in [50, 51] through a covariant formulation, allowing specific entropy variations as usually expected for nonequilibrium processes. The cosmological implications of the irreversible matter creation processes in the framework of the thermodynamics of open or quantum systems have been intensively investigated in [52–75].

From the analysis of some classes of scalar field potentials, including the potential, which provide a good fit to the CMB measurements, in [76] it was found that the Planck 2015 68% confidence limit upper bound on the spectral index, , implies an upper bound on the reheating temperature of GeV. The observations also exclude instantaneous reheating. Hence the low reheating temperatures allowed by some scalar field models open the possibility that dark matter could have been produced during the reheating period, instead of the radiation dominated era in the history of the Universe. Such a possibility could lead to a major change in our understanding of the early evolution stages of the Universe, leading to different predictions for the relic density and momentum distribution of WIMPs, sterile neutrinos, and axions.

It is the purpose of the present paper to apply the thermodynamics of open systems, as introduced in [49, 50], to a cosmological fluid mixture, in which particle decay and production occur, consisting of two basic components: scalar field, and matter, respectively. From a cosmological point of view this physical situation is specific to the preheating/reheating period of the inflationary cosmological models. The thermodynamics of irreversible processes as applied to inflationary cosmological models leads to a self-consistent description of the matter creation processes, which determines the whole dynamics and future evolution of the Universe. By combining the basic principles of the thermodynamics of open systems with the cosmological Einstein equations for a flat, homogeneous, and isotropic Universe we obtain a set of first-order ordinary differential equations describing matter creation due to the decay of the scalar field. As for the matter components we specifically consider a model in which ordinary matter is a mixture of two components, a relativistic one (radiation) and a pressureless fluid, assumed to correspond to dark matter. In order to simplify the analysis of the reheating equations we reformulate them in a dimensionless form, by introducing a set of appropriate dimensionless variables. As a cosmological application of the general formalism we investigate matter creation from a scalar field by adopting different mathematical forms of the self-interaction potential of the field. More exactly, we will consider power law and exponential and Higgs type potentials. For each case the evolution of the postinflationary Universe with matter creation is investigated in detail by solving numerically the set of the cosmological reheating evolution equations. The results display the process of inflationary scalar field decaying to matter and dark matter, and the effects of the expansion of the Universe on the decay. We concentrate on the evolution of the dark matter and radiation components, which increase from zero to a maximum value, as well as on the scalar field decay. A simple model, which can be fully solved analytically, is also presented. As a general result we find that the scalar field potential also plays an important role in the reheating process, and in the decay of the scalar field. Some of the cosmological parameters of each model are also presented and analyzed in detail. The models are constrained by the observational parameters called the inflationary parameters, including the scalar spectral index , the tensor-to-scalar ratio , the number of -folds , and the reheating temperature . These parameters can be calculated directly from the potentials used in the models, and they can be used together with the observational data to constrain the free parameters in the potentials. The study of the reheating models and the development of new formalisms and constraints to the theories are undoubtedly an essential part of the completion to the theory. Analysis of the current experimental data, including the calculation of the inflationary parameters, is able to test the reheating models.

The present paper is organized as follows. The basic results in the thermodynamics of open systems are briefly reviewed in Section 2. The full set of equations describing matter creation due to the decay of a scalar field are obtained in Section 3, where the observational constraints on the inflationary and reheating models are also presented. An exact solution of the system of the reheating equations, corresponding to simple form of the scalar field (coherent wave), is obtained in Section 4. Several reheating models, corresponding to different choices of the scalar field self-interaction potential, are considered in Section 5, by numerically solving the set of reheating evolution equations. The time dynamics of the cosmological scalar field and of the matter components is obtained and analyzed in detail. We conclude and discuss our results in Section 6.

#### 2. Irreversible Thermodynamics of Matter Creation

In this Section, we briefly review the thermodynamics of matter creation in open systems, which we will use for the study of the postinflationary reheating processes. In our presentation we start from the fundamental laws of thermodynamics, we proceed to the general covariant formulation of the theory, and then we apply our results to the particular case of homogeneous and isotropic cosmological models.

##### 2.1. The First Law of Thermodynamics

The main novel element that appears in the thermodynamics of open systems with energy/matter exchange is the creation pressure, related to the variation of the particle number. The expression of the creation pressure can be derived from the fundamental laws of thermodynamics [49]. We consider a thermodynamic system consisting of particles, inside a volume . For a closed system, is a constant, and the first law of thermodynamics, the energy conservation equation, is where is the internal energy, is the heat received by the system, and is the thermodynamic pressure. We can also write the energy conservation as where the energy density , the particle number density , and , respectively. Adiabatic transformations are described by

For an open system, the particle number can vary in time, so that . By including the thermodynamic parameters of the newly created particles, the energy conservation equation is modified to [49] where is the enthalpy per unit volume. Adiabatic transformations are then expressed asThe matter creation process is a source of the internal energy of the system. From (8) we obtain the relationdescribing the variation of the matter density due to the creation processes. On the other hand the above equation can also describe physical systems in the absence of matter creation, like, for example, the photon gas, in which, for a particular temperature, the particle number varies with the volume in a fixed manner, adjusting itself to have a constant density of photons. Hence, for a photon gas the particle number is a variable, and not a constant quantity, like in an ordinary gas. The physical mechanism by which thermodynamic equilibrium can be established consists in the absorption and emission of photons by matter. For a radiation like equation of state, with , from (9) we obtain , with and constants, giving the well-known energy density-particle number relation for the ultrarelativistic photon gas.

Alternatively, we can reformulate the law of the energy conservation by introducing the matter creation pressure, so that the energy conservation equation for open systems takes a form similar to (6), with the pressure written as , where the creation pressure associated describing the creation of the matter is

##### 2.2. The Second Law of Thermodynamics

In order to formulate the second law of thermodynamics for open systems, we define first the entropy flow , and the entropy creation , so that the total entropy variation of the system is written as [49] For a closed and adiabatic system, . For an open system with matter creation, in the following we assume that matter is created in a thermal equilibrium state . Therefore the entropy increases only due to the particle creation processes.

The total differential of the entropy is where the entropy density and the chemical potential . Hence the entropy production for an open system is

The matter creation processes can also be formulated in the framework of general relativity. In order to apply the thermodynamics of open systems to cosmology, in the following we consider the formulation of the thermodynamics of open systems in a general relativistic covariant form [50]. Then, by starting from this formulation, we particularize it to the case of a homogeneous and isotropic Universe, by also deriving the entropy of the system [72].

##### 2.3. General Relativistic Description of the Thermodynamic of Open Systems

In the following we denote the energy-momentum tensor of the system as , and we introduce the entropy flux vector , the particle flux vector , and the four-velocity of the relativistic fluid. These variables are defined aswhere is the matter creation pressure, is the specific entropy per particle, and is the particle number density.

The energy conservation law requireswhile the second law of thermodynamics imposes the constraint on the entropy flux four-vector. The balance equation of the particle flux vector is given bywhere represents a particle source, or a particle sink if . For an open system, the Gibbs equation is [50]In order to derive the energy balance equation we multiply both sides of (17) by , thus obtaining where we have used the relations , , and , respectively. Hence we have obtained the energy conservation equation for open systems as

In order to obtain the entropy variation we substitute the relation (22) into the Gibbs equation (20), and we use (16) and (19) to obtain first Thus we obtain the entropy balance equation aswhere is the chemical potential and is the expansion of the fluid. From (16) we obtainwhere we have constrained the specific entropy per particle to be a constant, . From (24) and (25) it follows that and thus we obtain matter creation pressure asOne can see from the above expression that if there is no matter creation, , and then the matter creation pressure vanishes, and the entropy becomes a constant. If matter is created with and , then the matter creation pressure will be negative, .

##### 2.4. Matter Creation in Homogeneous and Isotropic Universes

In the case of a homogeneous and isotropic flat Friedmann-Robertson-Walker (FRW) geometry, the line element is given bywhere is the scale factor. In the comoving frame, the four-velocity is given by . For any thermodynamic function its derivative is given by . For the covariant derivative of the four-velocity, by taking into account the fact that , we find where , the expansion of the fluid (the Hubble function), is given by . The components of the energy-momentum tensor (15) in the presence of particle creation are

From (19) we obtain the particle number balance equation as Hence, by substituting the above equation into (14), we obtain the variation of the entropy of the open system due to particle creation processes as where is an arbitrary integration constant.

#### 3. Reheating Dynamics: Irreversible Open Systems Thermodynamics Description

Reheating is the period at the end of inflation when matter particles were created through the decay of the inflationary scalar field. The particle creation is essentially determined by the energy density and pressure of the inflationary scalar field . For all inflationary models, inflation begins at the energy scale of order GeV, or, equivalently, at an age of the Universe of the order of s [3]. Hence reheating starts at the end of inflation, at a time of the order of s, and it should end before the hot big bang at s, which means that the energy scale of reheating is from around GeV up to GeV. The reheating temperature is an important parameter of the model, which we will discuss in detail.

In this Section we derive first the general equations describing the reheating process in the framework of the thermodynamics of open systems. We also review the main observational constraints and parameters of the inflationary cosmological models, prior and during the reheating.

##### 3.1. Irreversible Thermodynamical Description of Reheating

We consider the Universe in the latest stages of inflation as an open thermodynamic system with matter creation. The three constituents of the Universe are the inflationary scalar field, the ordinary matter, with energy density and pressure , and the dark matter, with energy density and pressure , respectively. Both the ordinary matter and the dark matter are assumed to be perfect fluids. The total energy density of the system is Similarly, the total pressure and the total particle number density of the Universe at the end of inflation are given by

In the presence of particle creation processes, the stress-energy tensor of the inflationary Universe is obtained as where is the total creation pressure describing matter creation, and we have assumed that all three fluid components are comoving, thus having the same four-velocity . In the following we will restrict our analysis to the case of the flat Friedmann-Robertson-Walker Universe, with line element given by (28). The basic equations describing the reheating process in the framework of the thermodynamic of open systems can be formulated as follows.

*(a) The Particle Number Balance Equation*. In the FRW geometry the balance equation of the particle flux vector, (19), can be written for each species of particle with particle numbers , , and as

We assume that* the particle source function ** is proportional to the energy density of the inflationary scalar field *, with the proportionality factor representing a decay (creation) width. Hence the balance equations of the particle number densities are given by where , , and are generally different, nonzero functions, and , , and denote the masses of the scalar field, dark matter, and ordinary matter particles, respectively. For the whole system, the particle balance equation is

*(b) The Energy Balance Equations*. From (9) we obtain the equations of the energy density balance of the cosmological fluids asrespectively.

*(c) The Matter Creation Pressures*. From the expression of the creation pressure associated with the newly created particles Eq. (27) we find where superscript denotes the different cosmological fluids. During the reheating phase, for each component, the decay and creation pressures areThe total creation pressure is *(d) The Gravitational Field Equations*. The Einstein gravitational field equations in a flat FRW Universe with inflationary scalar field, ordinary matter, and dark matter fluids in the presence of matter creation can be written as follows:In the present paper we use natural units with , in which MeV^{−2}, where is the “reduced” Planck mass.

*(e) The Entropy Generation*. The entropy generated by the newly created matter is

##### 3.2. The Models of Reheating

Presently, the accepted paradigm for reheating is based on particle production and parametric amplification, as extensively discussed in [46, 47]. According to this model, the main process of particle production is parametric amplification in combination with a slowing expansion. This mechanism entails a description in terms of Mathieu equations, either for the scalar or fermion fields. These equations feature unstable bands with the concomitant parametric instabilities. This mechanism cannot be described in terms of a simple functions, and different bands feature different Floquet exponents (in Minkowski or near Minkowski spacetime) and correspond to different wave-vector bands. In this picture, matter is created in a highly nonthermal state (determined primarily by the position and widths of the unstable parametric bands) and thermalizes (at least in principle, although never quite proven directly) via scattering. This last step requires a clear knowledge of relaxation times, which entail understanding the microscopic processes and cross sections.

In the following we will consider models of reheating involving energy transfer between the inflationary scalar field and the newly created ordinary matter and dark matter, respectively. We assume that particles are created through an irreversible thermodynamic process. In this Subsection, we present explicitly the general equations governing the matter and dark matter creation during inflationary reheating. In our reheating models, the energy density and pressure for the inflationary scalar field have the formwhere is the self-interaction potential of the scalar field. Substituting (48) into (39), we obtain the generalized Klein-Gordon equation for the scalar field in the presence of matter creation. The equation of state of the dark matter is defined as , since we assume that the pressure of the newly created dark matter particles is . The energy density of dark matter is proportional to its particle number density, , so that , where is the mass of the dark matter particles. For pressureless matter from (9) we obtain the general relation

The newly created matter is considered to be in form of radiation (an ultrarelativistic fluid), with the equation of state . The relation between the energy density and the particle number density of radiation is .

In summary, the set of equations governing the matter creation process during reheating is written below as where , , and are considered as constants. Eq. (52) is the generalized Klein-Gordon equation, describing the evolution of the scalar field. The last term on the left hand side of this equation represents an effective friction term, describing the effects of irreversible matter creation (decay) on the scalar field dynamics. If , we recover (2), introduced on a purely phenomenological basis in [26].

###### 3.2.1. Dimensionless Form of the Reheating Equations

In order to simplify the mathematical formalism we introduce a set of dimensionless variables , , , , , , , , , and , respectively, defined as where is a constant with the physical dimension of a mass.

In the new variables (50)–(54) take the dimensionless form

Eqs. (59) and (60) can be solved exactly, and we obtain the following exact representation of the scalar field particle number and matter densities: When , the particles creation processes end, and the Universe enters in the standard expansionary phase, with its chemical composition consisting of two matter components, as well as a possible background of scalar field particles. After reheating, the energy density of the dark matter and radiation varies as where and are the dark matter and radiation energy densities at the end of reheating. The evolution of the scale factor is described by the equation giving The deceleration parameter for the postreheating phase is given by In the limit , the deceleration parameter tends to the value , indicating a decelerating expansion of the Universe.

##### 3.3. Observational Tests of Inflationary Models

In order to facilitate the comparison of the theoretical models of the reheating based on the thermodynamics of open systems with the observational results, we also present and calculate various cosmological parameters, either directly from the inflationary models or from the potentials of the inflationary scalar field.

*(a) Slow-Roll Parameters*. The slow-roll parameters were introduced in 1992 and are based on the slow-roll approximation [18]. Their calculations are presented in [4]. During the period before reheating, the inflationary scalar field has not yet decayed, and the field equation and the equation of the scalar field can be written as The slow-roll approximation requires that the terms and be sufficiently small, so that the above equations can be simplified toThe slow-roll parameters are defined as and , respectively [4, 18]. The slow-roll approximation only holds when the following conditions are satisfied: and , respectively.

*(b) Inflationary Parameters*. During inflation, the Universe expands exponentially. The scale of inflation is represented by the scale factor . The ratio of the scale factor at the end of inflation and at its beginning is always large, so that the logarithm of this quantity is of physical interest. The number of -folds is defined as , where is the time when inflation ends. The number of -folds can also be obtained from the comoving Hubble length [4]. Under the slow-roll approximation, it can be calculated as where is defined as , when the slow-roll approximation is violated, and inflation ends. Typically . Sometimes the integral is multiplied by a factor of .

The scalar spectral index is defined by the slow-roll parameters . The index is usually a function of the parameters of the model. The value of can be obtained from the observations [16] and is used to constrain the inflationary models. The tensor-to-scalar ratio is similarly defined as . One can see from their definition that the two parameters and are not independent.

*(c) Cosmological Parameters*. The deceleration parameter is defined as It indicates the nature of the acceleration of the expansion of the Universe. means that the expansion of the Universe is accelerating, while a positive indicates deceleration.

The density parameters of the inflationary scalar field, matter, and dark matter are defined as , , and . The evolution of density parameters of each matter component in the early Universe gives us a hint on the matter creation process during inflationary reheating.

###### 3.3.1. Model Parameters

In studying the reheating models, we consider different scalar field potentials that influence the dynamics of the early Universe. During inflation the energy density of the scalar field is potential energy dominated.

*(a) Scalar Field Potentials.* We picked six potentials of cosmological interest that we are going to investigate and presented them in Table 1.