- Research article
- Open Access
- Published:

# Dilaton and off-shell (non-critical string) effects in Boltzmann equation for species abundances

*PMC Physics A***volume 1**, Article number: 2 (2007)

## Abstract

In this work we derive the modifications to the Boltzmann equation governing the cosmic evolution of relic abundances induced by dilaton dissipative-source and non-critical-string terms in dilaton-driven non-equilibrium string Cosmologies. We also discuss briefly the most important phenomenological consequences, including modifications of the constraints on the available parameter space of cosmologically appealing particle physics models, imposed by recent precision data of astrophysical measurements.

**PACS codes:** 11.25.Wx, 98.80.Ft, 95.35.+d

## I Introduction

In a previous work [1] we have discussed a case study of dissipative Liouville-string cosmology, involving non critical string cosmological backgrounds, with the identification [2, 3] of target time with the world-sheet zero mode of the Liouville field [4–7]. Such cosmologies were found to asymptote (in cosmic time) with the conformal backgrounds of [8–10], which are thus viewed as equilibrium (relaxation) configurations of the non-equilibrium cosmologies.

In terms of microscopic considerations, a model for such departure from equilibrium could be considered the collision of two brane worlds, which results [11, 12] in departure from conformal invariance of the effective string theory on both the bulk and the brane world, and thus in need for Liouville dressing to restore this symmetry [4–7]. Dynamical arguments [13, 14], then, stemming from minimization of effective potentials in the low-energy string-inspired effective field theory on the brane, imply the (eventual) identification of the zero mode of the Liouville mode with (a function) of target time [2, 3]. There is an inherent *time irreversibility* in the process, which is associated with basic properties of the Liouville mode, viewed as a local (dynamical) renormalization group (RG) scale on the world-sheet of the string. This implies *relaxation* of the associated dark energy of such cosmologies, and a gravitational friction, associated with the conformal theory central charge deficit. For a recent review we refer the reader to [11, 12], where concepts and methods are outlined in some detail.

For our purposes in this letter we would like to concentrate on one interesting aspect of the non-critical string cosmologies, associated with the off-equilibrium effects on the Boltzmann equation describing relic abundances and the associated particle-physics phenomenology. Indeed, in conventional cosmologies, a study of relic abundances by means of Boltzmann equation that governs their cosmic time evolution yields important phenomenological constraints on the parameters of particle physics (supersymmetric) models using recent (WMAP [15] and other) astrophysical Cosmic Microwave Background (CMB) data. Essentially, the astrophysical data constraint severely in some cases the available phase-space distributions of favorite supersymmetric dark matter candidates such as neutralinos [16–19]. When off-equilibrium, non-critical string cosmologies are considered [1], which notably are consistent with the current astrophysical data from supernovae, as demonstrated recently [20], then there are significant modifications to the Boltzmann equation, stemming from extra sources (dilaton) and off-shell (non-critical, non equilibrium) terms, which affect the time evolution of the phase-space density of the species under consideration.

It is the purpose of this paper to derive such modifications in detail, and then use them in order to discuss briefly particle-physics models constraints, especially from the point of view of supersymmetry. With regards to this last issue, in this article we shall present an approximate analytical treatment, which will only provide hints to what may actually happen. A complete analysis requires numerical studies which are postponed for a future publication.

## 2 Modified Boltzmann Equation in Non-Critical Strings

Consider the phase space density of a species X, which is assumed coupled to the off-shell (non-critical string) background terms:

where quantities refer to the Einstein frame, and the Einstein metric is assumed to be of Robertson-Walker (RW) type. Throughout this work we follow the normalization and conventions of [8–10].

For completeness we state the main relationships between the Einstein and *σ*-model frames [8]:

where Φ is the dilaton field, and the superscript *σ* denotes quantities evaluated in the *σ*-model frame.

Since in the RW Einstein frame *g*_{00} = -1, and in the cosmic Einstein co-moving frame we have for a generic massive species with mass *m* (that we shall be interested in this work) ${p}_{\sigma}^{\mu}$ = *mdx*^{μ}/*dτ* = *m*(*dx*^{μ}/*dt*) (*dt*/*dt*_{
σ
}) = *mp*^{μ}*e*^{-Φ}, it can be readily seen that

We shall be interested in the action of the relativistic Liouville operator $\widehat{L}$, acting on a phase-space density of species X, *f*(|$\overrightarrow{p}$|, *t*, *g*^{i}), which in general depends on the off-shell backgrounds *g*^{i}= (*g*_{
ii
}, Φ), where *g*_{
ii
}are the spatial components of the RW space-time metric in Einstein frame, and Φ is the dilaton. The off-shell dependence comes about due to the interpretation of target time as a local world-sheet renormalization-group scale, the Liouville field.

We commence our analysis by recalling the relativistic form of the Liouville operator in conventional general relativity:

The second term on the r.h.s. is the relativistic form of the force, following from the geodesic equation. For a RW Universe, only the time-energy part survives from the first term, i.e.

Moreover, the connection part receives non trivial contributions only from the terms $\sum}_{i}{\Gamma}_{ii}^{0}{p}^{i}{p}^{i}\frac{\partial}{\partial E$, *i* = 1, 2, 3 a spatial index (assuming for concreteness an already compactified string theory, or a theory on a three-brane world), with ${\Gamma}_{ii}^{0}=-a\dot{a}$, where the overdot denotes derivative w.r.t. cosmic RW time *t*, identified in our string theory with the Einstein-frame RW time (2.1).

Hence, the conventional part of the Liouville operator in a RW Universe would read [21]:

The presence of off-shell, non-critical (Liouville) string [4–7] backgrounds, upon the (dynamical) identification [2, 3, 11–14] of the Liouville mode with the target time, one will receive extra contributions in the expression for the associated Liouville operator, stemming from the fact that the latter is nothing but a total time derivative.

In the context of our string discussion it is important to specify that the initial frame, from which we commence our discussion, is the so-called Einstein frame. It is in this frame that the usual RW cosmology is obtained in string theory [9], for which the expression (2.5) is valid.

To discuss the non-critical-string (off-shell) corrections to Boltzmann equation, it will be necessary to consider time derivatives in the *σ*-model frame. This is due to the fact that it is in this frame that the target time *X*^{0} ≡ *t*_{
σ
}is related simply to the Liouville mode ϕ in non-critical string theories [2, 3],

*ϕ*+ t_{σ}= 0.

This relation is obtained dynamically from minimization arguments of the low-energy effective potential of some physically interesting cosmological models, for instance those involving colliding brane worlds [11, 13]. To be precise, the initial relation (2.6) derived in the specific model of [13] reads: *ϕ*/$\sqrt{2}$ + *t*_{
σ
}= 0, with the factor of $\sqrt{2}$ arising from (logarithmic) conformal-field-theory considerations, and being crucial for yielding a Minkowski space-time in the scenario of [13], where the coordinate *X*^{0} assumed an initial Euclidean signature, appropriate for a world-sheet path-integral quantization. Nevertheless, in our approach below we shall absorb the $\sqrt{2}$ into the normalization of the units of the Regge slope *α'*, for convenience.

In such scenarios the target-space dimensionality of the string is extended to *D* + 1 initially, with two time-like coordinates, *t* and *ϕ*. Eventually, these two coordinates are identified (c.f. (2.6)) in the solution of the generalized conformal invariance conditions that express the restoration of the conformal invariance by the Liouville mode, as we discuss below. This implies that, initially, the Liouville operator (2.3) acquire extra Liouville components (*ϕ*, *p*^{ϕ}), with *p*^{ϕ}≡ *E*_{
ϕ
}. Hence, there are additional structures on the right hand side of (2.3), of the form ${E}_{\varphi}\frac{\partial}{\partial \varphi}$, together with the corresponding amendments in the Christoffel-symbol dependent ∂/∂*p*^{ϕ}terms, stemming from the extra Liouville coordinate on the target space time, which implies appropriate extra components of the extended-target-space-time metric. The (dynamical) implementation of (2.6) restricts oneself to a *D*-dimensional hypersurface in the extended target-space, with the identification of *E*_{
ϕ
}→ *E*, the ordinary (physical) energy. From now on we restrict ourselves on this hypersurface, bearing in mind however that the identification of the Liouville mode *ϕ*with the target time should only be implemented at the end, and thus any dependence on *ϕ*should be kept explicit at intermediate steps of the pertinent calculations.

The connection between time derivatives in the Einstein and stringy-*σ*-model frames is provided by the chain rule of differentiation $\frac{\partial}{\partial {t}_{\sigma}}=\frac{\partial t}{\partial {t}_{\sigma}}\frac{\partial}{\partial t}$, using (2.1). This will result in multiplicative factors of *e*^{-Φ} in front of the appropriate non-critical-string modifications of the Liouville operator (2.5).

Since, as mentioned above, the Liouville operator is essentially a total time derivative operator, and in our case time is related to a world-sheet renormalization group scale, *ρ*(*σ*, *τ*), taken to be local on the two-dimensional surface as a result of world-sheet general covariance [22, 23], the sought-for non-critical-string modifications of the Boltzmann equation emerge from the implicit dependence of the *g*^{i}background fields on *ρ*(*σ*, *τ*), that is the corresponding *β*-functions, which in a critical-string theory would vanish.

We now remark that in the approach of [2, 3], the local RG scale *ρ*(*ξ*) was identified with the dynamical Liouville mode *ϕ*, which implied that in this approach the local conformal invariance of the world-sheet of the string was restored [4–7],

*ρ* ≡ *ϕ*.

With these in mind we then modify the relativistic form (2.3), (2.4) by replacing

where *ξ* = *σ*, *τ*, denotes the world-sheet coordinates, ∫_{Σ} is a world-sheet integration, and the index "*i*" runs on both, a (discrete) background field space, {*g*_{
ii
}(*y*), Φ(*y*)}, *i* = 1, 2, 3 and a continuous *D*-dimensional (target) space-time *y*. Hence, summation over "*i*″ includes integration $\int {d}^{D}}y\sqrt{-g}\left(\mathrm{...}\right)$. Note that since the (physical) target time in the *σ*-model frame, *t*_{
σ
}, is only the world-sheet zero mode of the corresponding *σ*-model field *X*^{0}(*ξ*), eventually the world-sheet integration ∫_{Σ} in (2.8) disappears, since only the world-sheet zero mode *ρ* of the local RG scale *ρ*(*ξ*) (Liouville mode) will yield a non zero contribution. In the last equality of the right-hand side of (2.8), the quantity *η* is defined as: $\eta \equiv \frac{\partial \rho}{\partial {t}_{\sigma}}$. If we use the identification *ρ*(*ξ*) = *ϕ*(*ξ*) and Eq. (2.6), then *η* = -1, but at this stage we keep it general to demonstrate explicitly the renormalization-scheme dependence (i.e. choice of the local scale *ρ*(*ξ*)).

The non-critical string contributions $\frac{\partial {g}^{i}}{\partial \rho}={\tilde{\beta}}^{i}$ are the Weyl anomaly coefficients of the string, but with the renormalized *σ*-model couplings *g*^{i}being replaced by the Liouville-dressed [4–7] quantities. This is a particular feature of the approach of [2, 3] viewing the Liouville mode as a local world-sheet RG scale, as mentioned previously. In turn the latter is also identified (c.f. (2.6)) with (an appropriate function of) the target time *t*.

The dynamics of this latter identification is encoded in the solution of the generalized conformal invariance conditions, after Liouville dressing, which read in the *σ*-model frame [2–7]:

where the prime denotes differentiation with respect to the Liouville zero mode *ρ*, and the overall minus sign on the left-hand side of the above equation pertains to supercritical strings [8, 9], with a time like signature of the Liouville mode, for which the central charge deficit *Q*^{2} > 0 by convention.

Notice the dissipation, proportional to the (square root) of the central charge deficit *Q*, on the right-hand side of (2.9), which heralds the adjective *Dissipative* to the associate non-critical-string-inspired Cosmological model. Moreover, the Weyl anomaly coefficients $\tilde{\beta}$^{i}, *i* = {Φ, *g*_{
μν
}}, whose vanishing would guarantee local conformal invariance of the string-cosmology background, are associated with off-shell variations of a low-energy effective string-inspired action, $\mathcal{S}$[*g*^{j}]

with *z*, $\overline{z}$ (complex) world-sheet coordinates, $\mathcal{G}$_{
ij
}the Zamolodchikov metric in theory space of strings, and *V*_{
i
}the *σ*-model vertex operators associated with the *σ*-model background field *g*^{i}. It is this off-shell relation that characterizes the entire non-critical (*Q*) Cosmology framework, associated physically with a non-equilibrium situation as a result of an initial cosmically catastrophic event, at the beginning of the (irreversible) Liouville/cosmic-time flow [11, 12].

The detailed dynamics of (2.9) are encoded in the solution for the scale factor *a*(*t*) and the dilaton Φ in the simplified model considered in [1], after the identification of the Liouville mode with the target time (2.6). In fact, upon the inclusion of matter backgrounds, including dark matter species, the associated equations, after compactification to four target-space dimensions, read in the Einstein frame [1]:

where ${\tilde{\varrho}}_{m}({\tilde{p}}_{m})$ denotes the matter energy density (pressure), including dark matter contributions, and $\varrho $_{Φ} (*p*_{Φ}) the corresponding quantities for the dilaton dark-energy fluid. All derivatives in (2.11) are with respect the Einstein time *t* which is related to the Robertson-Walker cosmic time *t*_{
RW
}by *t* = *ωt*_{
RW
}. Without loss of generality we have taken *ω* = $\sqrt{3}$*H*_{0} where *H*_{0} is the Hubble constant. With this choice for *ω* the densities appearing in (2.11) are given in units of the critical density. Then one can see for instance that if the time *t*_{
RW
}is used the first of equations (2.11) above receives its familiar form in the RW geometry when the contributions of the dilaton and non-critical terms are absent.

The overdots in the above equations denote derivatives with respect to the Einstein time. Their right-hand side contain the non-critical string *off-shell* terms:

In the above equations *H* = ($\dot{a}$)/*a* is the Hubble parameter and *q* is the deceleration parameter of the Universe $q\equiv -\ddot{a}a/{\dot{a}}^{2}$, and are both functions of (Einstein frame) cosmic time. The potential $\widehat{V}$_{
all
}appearing above is defined by $\widehat{V}$_{
all
}= 2*Q*^{2} exp(2Φ) + exp(4Φ)*V* where, in order to cover more general cases, we have also allowed for a potential term in the four-dimensional action $-{\displaystyle \int {d}^{4}y\sqrt{-G}}V$ in addition to that dependent on the central charge deficit term *Q*. Although we have assumed a (spatially) flat Universe, the terms on the r.h.s., which manifest departure from the criticality, act in a sense like curvature terms as being non-zero at certain epochs. The dilaton energy density and pressure are given in this class of models by:

Notice that the dilaton field is not canonically normalized in this convention and its dimension has been set to zero.

For completeness, we mention at this point that the dependence of the central charge deficit *Q*(*t*) on the cosmic time stems from the running of the latter with the world-sheet RG scale [2, 3, 11, 12], and is provided by the Curci-Paffuti equation [24] expressing the renormalizability of the world-sheet theory. To leading order in an *α'* expansion, which we restrict ourselves in [1] and here, this equation in the Einstein frame reads:

For future use we also state here the corresponding continuity of the matter stress tensor, which is not an independent equation, but can be obtained from (2.11) by appropriate algebraic manipulations:

This expresses the non-conservation equation of matter as a result of its coupling to *both*, the dilaton source terms and the off-shell, non-equilibrium (non-critical-string) backgrounds.

A consistent solution of *a*(*t*), Φ(*t*), and the various densities, including back reaction of matter onto the space-time geometry, has been discussed in [1], where we refer the interested reader for further study. Also note that a preliminary comparison of such non-critical strings theories with astrophysical data, demonstrating consistency at present, is given in [20]. It should be remarked at this point that the solutions obtained in this study tend asymptotically to those presented in [8–10]. In such an approach, the dilaton rolls down asymptotically in (Einstein) target time as Φ = const. - ln *t*, consistent with a vanishing string coupling *g*_{
s
}= *e*^{Φ} → 0 as *t* → +∞. This situation is opposite to other studies in critical (equilibrium) string cosmologies, where the dilaton gets a non-trivial vacuum expectation value (v.e.v.) through an appropriate minimisation of its potential, being thus stabilised at early times. In such on-shell situations, as opposed to our non-equilibrium (off-shell) model studied in this and in previous articles [1, 25, 26], matter particles may acquire dilaton-v.e.v.-dependent masses, which may be heavy, and as such they decouple from a thermal bath. In contrast, as we shall discuss in the next section, our non-equilibrium rolling dilaton and running string coupling, implies that at late eras of the Universe, one may still evaluate the relic abundances of matter particles as in conventional cosmology, using an appropriately modified Boltzmann equation.

A final, but important comment should be made at this juncture. It is understood that, in order not to disturb the delicate balance between particle interactions and the expansion of the Universe during the nucleosynthesis era, it is necessary to consider models in which the dilaton may be almost constant during that era, in such a way that the associated string coupling and particle decay and/or recombination rates are not affected much by the dilaton and non-critical string terms, as compared with the conventional cosmology situation. This important and by far not complete issue is highly model dependent, and it relies on specific properties of the underlying microscopic string theory, involving detailed scenaria of compactification, moduli field stabilization *etc*. In the literatute there are toy brane-inspired, non-critical-string cosmological models, where such approximately constant dilatons during early epochs of the Universe could be realised explicitly [13, 14, 25]. We also remark that similar assumptions are made in generic dilaton quintessence scenaria of critical strings [27].

After this necessary digression we now come back to discussing the derivation of the dilaton-source and noncritical-string induced modifications to the Boltzmann equation, governing the cosmic evolution of the various species densities. It is important for the reader to bear in mind already at this stage that the Boltzmann equation does not contain any independent information from the dynamical equations (2.11), but it should be rather viewed as an effective way of describing the cosmic evolution of the density of a given species, consistent with the continuity equation (2.15) for the total matter energy density. We shall come back to this important point later on.

At the moment, let us concentrate first on the dilaton-source and non-critical-string background contributions to the Liouville operator (2.8). The presence of (time-dependent) dilaton source terms implies an explicit dependence of the phase space density of a species *f* on Φ, while the non-conformal on the world-sheet) nature of the metric and dilaton background, induce a Liouville mode *ρ* dependence through the corresponding backgrounds:

The non-critical string terms can be expressed, as we have seen above, in terms of the corresponding Weyl anomaly coefficients, which are non zero as a result of departure from conformal invariance of the pertinent string background. Despite their non-linear looking appearance when expressed in terms of the Weyl anomaly coefficients (which depend on the Ricci tensor and second covariant derivatives of the Dilaton field), such terms acquire a particularly simple linear form once the identification of time with the Liouville mode (2.6) is implemented, which in effect implies that the explicit solution of (2.9), (2.11) must be taken into account when discussing the Boltzmann equation.

With the above in mind, the form of the non-critical-string (off-shell) and dilaton-source corrections to the Liouville operator reads (expressing quantities in the Einstein frame, taking into account the fact that in the *σ*-model frame ${g}_{\mu \nu}^{\sigma}$ and Φ are treated as independent field-theory variables, and using (2.1)):

where the overdot denotes the derivative with respect to the (Einstein frame) cosmic time *t*, and the prime the derivative with respect to *ρ*, and we used the fact that:

with *g*_{00} = -1, *g*_{
ii
}= *a*^{2}, and $\frac{\partial {g}_{ii}^{\sigma}}{\partial \rho}={\tilde{\beta}}_{ii}^{\text{Grav}}$ the corresponding Weyl anomaly coefficient, in order to express ${\tilde{\beta}}_{ii}^{\text{Grav}}$ in terms of quantities in the Einstein frame. We have also been careful to separate the dilaton dependence on *ρ*(*t*) from that on the ordinary time *t* (prior to the identification (2.6)). We remind the reader, at this point, that the identification (2.6) should only be applied at the end of the calculations, and should be viewed as a specific choice of "renormalization" scheme. The "renormalization" scheme (2.6), i.e. *η* = -1, which is obtained by a physical principle (dynamics) in certain models of interest to us here [13, 14], may thus be viewed in this context as the "physical scheme", where observational cosmology can be analyzed. In such a scheme there is a simplification of the complicated non-linear expressions encoded in the Weyl anomaly coefficients.

Taking into account the implicit dependence of *f*(|$\overrightarrow{p}$|, *t*, Φ) on *g*_{
ii
}through $\left|\overrightarrow{p}\right|={\left({\displaystyle {\sum}_{i=1}^{3}{p}^{i}{p}^{i}{g}_{ii}}\right)}^{1/2}$, it is evident that $\partial f/\partial {g}_{ii}=(\partial \left|\overrightarrow{p}\right|/\partial {g}_{ii})\partial f/\partial \left|\overrightarrow{p}\right|=\frac{{p}_{i}{p}_{i}}{2\left|\overrightarrow{p}\right|}\frac{\partial f}{\partial \left|\overrightarrow{p}\right|}$, from which we can easily see that the final form of the action of the Liouville operator, including dilaton/source terms and non-critical string corrections, on the phase-space density *f* (2.16), is such that:

Upon considering the action of the above operator on the density of a given species *X*, *n* ≡ ∫ *d*^{3}*pf* we then arrive, after some straightforward momentum integration by parts, at the *modified* Boltzmann equation for a four-dimensional effective field theory after string compactification (or restriction on three-brane worlds), in the presence of non-critical (off-shell) string backgrounds and dilaton source terms:

The collision term *C*[*f*] assumes the usual form in conventional particle cosmology [21]. The reader is invited to compare the final equation for the cosmic time evolution of the density *n* in the second line of (2.20) to the continuity equation for the total energy density of matter (2.15). As discussed previously, the Boltzmann equation (2.20) should be compatible in the sense of leading to no extra information) with the conservation equation (2.15), as well as the (modified) Einstein equations (2.11).

In the above scenario, where the non-critical string contributions have been obtained through the identification of time with the (world-sheet zero-mode of the) Liouville field *ρ* (2.6), there is some *universality* in the coupling of *all* matter species, including dark matter, to these off-shell background terms, which may be traced to the equivalence of all species coupled to gravity. This is to be contrasted with the more general phenomenological case studied in [1], where scenarios involving only the coupling of the non-critical string backgrounds with the exotic dark matter species have been considered as well.

The final issue to be discussed in this section pertains to the form of the dependence of *f* on the dilaton source terms, which would survive a dilaton-driven critical-string cosmology case, such as the one considered in [27]. We constrain this form by requiring that in the Einstein frame there are two types of dependence on Φ: (i) *explicit*, of the form *e*^{-4Φ}, arising from the fact that in our approach, the phase space density is constructed as a quantity in the *σ*-model frame of the string, which is then expressed in terms of quantities in the Einstein frame. As such, it is by definition (as a density) *inversely* proportional to the proper *σ*-model frame volume ${V}^{\sigma}={\displaystyle \int {d}^{4}x}\sqrt{-{g}^{\sigma}}\propto {e}^{4\Phi}$ on account of (2.1); (ii) *implicit*, corresponding to a dependence on Φ through the Einstein-frame metric *g*_{
ii
}(2.1). Hence the general structure of *f* is of the form:

This implies that:

where in the last step we have performed appropriate partial (momentum-space) integrations.

The final form of the Liouville operation (2.20), then, reads:

We now notice that non-critical terms can be expressed in terms of the Weyl anomaly coeffcients for the (*σ*-model) graviton and dilaton backgrounds as:

where we used the Einstein frame metric to contract indices, with $\tilde{\beta}$^{Grav} denoting the graviton Weyl anomaly coefficient. In (2.24) we have taken into account that in the class of models of [1, 25], we are concentrating in this work, ${\tilde{\beta}}_{00}^{\text{Grav}}$ = 0. Thus the Boltzmann equation finally becomes

It is important to recall once more that, upon Liouville dressing, which restores the conformal invariance of the model, the graviton and dilaton Weyl anomaly coefficients satisfy (2.9), which upon the identification of the Liouville mode with (a function of the cosmic) time yield (2.11).

We next remark that in the case with non-critical string backgrounds, consistency requirements between the two equations determine the back reaction effects of matter onto space time. For instance, in the absence of matter, in the specific (but generic enough) non-critical string model of [25], we have that the present-era dilaton and scale factor of the Universe are such that $\dot{Q}\simeq 0,\dot{\Phi}+\frac{\dot{a}}{a}=0$. In this particular case the off-shell terms at late eras are practically zero, and the system has reached equilibrium. The corresponding Boltzmann equation for the dilaton species densities is then satisfied in agreement with the corresponding continuity equation (2.15) in absence of matter.

However, in the presence of matter, we know from the analysis of [1] that $\dot{\Phi}+\frac{\dot{a}}{a}\ne 0$, and the right-hand-side of the Boltzmann equation (2.20), with *n* now representing matter species, yields terms proportional to *n*. Consistency with the continuity equation (2.15), then, implies that the off-shell (non-critical string background) terms in the latter should be proportional to the total number density of matter species (back reaction effects). This and other similar consistency checks should be made when considering solving such non-equilibrium cosmologies.

## 3 Phenomenology of Particle Physics Models and Modified Boltzmann Equation

In this section we consider solutions of the modified Boltzmann equation (2.23), or equivalently (2.25), for a particle species density *n* in the physically interesting case of supersymmetric dark matter species, such as neutralinos, viewed as the lightest supersymmetric particles (LSP). Such Cold Dark Matter candidates lead to a rich phenomenology of supersymmetric particle physics models. In the context of conventional Cosmology [16–19], some of these models can be constrained significantly by the recently available astrophysical data on Cosmic Microwave Background temperature fluctuations [15]. The calculation of relic abundances will be done in some detail, in order for the reader to appreciate better the rôle of the non conventional terms in (2.25).

It is convenient to write the Boltzmann equation for the density of species *n* in a compact form that represents collectively the dilaton-dissipative-source and non-critical-string contributions as external-source Γ(*t*)*n* terms:

where we work in the physical scheme (2.6) from now on, for which *η* = -1. Depending on the sign of Γ(*t*) one has different effects on the relic abundance of the species *X* with density *n*, which we now proceed to analyze. To find an explicit expression for Γ(*t*) in our case we should substitute the solution of (2.9), more specifically (2.11), analyzed in [1]. Regarding the form of (3.1) it is nice to see that the extra terms can be cast in a simple-looking form of a source term Γ(*t*)*n* including both the dilaton dissipation and the non-critical-string terms. Of course Γ(*t*) is complicated and requires the full solution outlined in [1].

In a more familiar form, the interaction term *C*[*f*] of the above modified Boltzmann equation can be expressed in terms of the thermal average of the cross section *σ* times the Moeller velocity *v* of the annihilated particles [21]

Before the decoupling time *t*_{
f
}, *t* <*t*_{
f
}, equilibrium is maintained and thus *n* = *n*_{
eq
}for such an era. However, it is crucial to observe that, as a result of the presence of the source Γ terms, *n*_{
eq
}no longer scales with the inverse of the cubic power of the expansion radius *a*, which was the case in conventional (on-shell) cosmological models.

To understand this, let us assume that *n* = ${n}_{eq}^{(0)}$ at a very early epoch *t*_{0}. Then the solution of the modified Boltzmann equations at all times *t* <*t*_{
f
}is given by

The time *t*_{0} characterizes a very early time, which is not unreasonable to assume that it signals the exit from the inflationary period. Soon after the exit from inflation, all particles are in thermal equilibrium, for all times *t* <*t*_{
f
}, with the source term modifying the usual Boltzmann distributions in the way indicated in Eq. (3.3) above. It has been tacitly assumed that the entropy is conserved despite the presence of the source and the non-critical-string contributions. In our approach this is an *approximation*, since we know that non-critical strings lead to entropy production. However, as argued in our previous works on the subject [2, 3], the entropy increase is most significant during the inflationary era, and hence it is not inconsistent to assume that, for all practical purposes, sufficient for our phenomenological analysis in this work, there is no significant entropy production after the exit from inflation. This is a necessary ingredient for our approach, since without such an assumption no predictions can be made, even in the conventional cosmological scenarios. Thus, the picture we envisage is that at *t*_{0} the Universe entered an equilibrium phase, the entropy is conserved to a good approximation, and hence all particle species find themselves in thermal equilibrium, despite the presence of the Γ source, which *slowly* pumps in or sucks out energy, without, however, disturbing the particles' thermal equilibrium.

From the above discussion it becomes evident that it is of paramount importance to know the behaviour of the source term at all times, in order to extract information for the relic abundances, especially those concerning Dark Matter, and how these are modified from those of the standard Cosmology. Before embarking on such an enterprize and study the phenomenological consequences of particular models predicting the existence of Dark Matter, especially Supersymmetry-based ones, we must first proceed in a general way to set up the stage and discuss how the relic density is affected by the presence of the non conventional source terms discussed above.

For the sake of brevity, we shall not deploy all the details of the derivation of the relic density, but instead demonstrate the most important features and results of our approach, paying particular attention to exhibiting the differences from the conventional case. Generalizing the standard techniques [21], we assume that above the freeze-out point the density is the equilibrium density as provided by Eq. (3.3), while below this the interaction terms starts becoming unimportant. It is customary to define *x* ≡ *T*/${m}_{\tilde{\chi}}$ and restrict the discussion on a particular species $\tilde{\chi}$ of mass ${m}_{\tilde{\chi}}$, which eventually may play the role of the dominant Dark Matter candidate. It also proves convenient to trade the number density *n* for the quantity *Y* ≡ *n*/*s*, that is the number per entropy density [21]. The equation for *Y* is derived from (3.2) and is given by

where ${G}_{N}=1/{M}_{Planck}^{2}$ is the four-dimensional gravitational constant, the quantity *H* is the Hubble expansion rate, *h* denote the entropy degrees of freedom, and $\u27e8$*vσ*$\u27e9$is the thermal average of the relative velocity times the annihilation cross section and ${\tilde{g}}_{eff}$ is simply defined by the relation [21]

The reader should notice at this point that Δ$\varrho $incorporates the effects of the additional contributions due to the non-critical (off-shell) terms and the dilaton dissipative source, which are not accounted for in the *g*_{
eff
}of conventional Cosmology [21], hence the notation ${\tilde{g}}_{eff}$. We next remark that *ρ*, as well as Δ*ρ*, as functions of time are known, once one solves the cosmological equations. However, only the degrees of freedom involved in *ρ* are thermal, the rest, like the comological-constant term if present in a model, are included in Δ*ρ*. Therefore, the relation between temperature and time is provided by

while *ρ* + Δ*ρ* are involved in the evolution through (c.f. (2.11))

Thus, it is important for the reader to bear in mind that Δ*ρ* contributes to the dynamical expansion, through Eq. (3.7), but not to the thermal evolution of the Universe. The quantity ${\tilde{g}}_{eff}$, defined in (3.5), is therefore given by

The meaning of the above expression is that time has been replaced by temperature, through Eq. (3.6), after solving the dynamical equations. In terms of ${\tilde{g}}_{eff}$ the expansion rate *H* is written as

This is used in the Boltzmann equation for *Y* and the conversion from the time variable *t* to temperature or, equivalently, the variable *x*.

For *x* above the freezing point *x*_{
f
}, *Y* ≈ *Y*_{
eq
}and, upon omitting the contributions of the derivative terms *dh*/*dx*, an approximation which is also adopted in the standard cosmological treatments [21], we obtain for the solution of (3.4)

Here, ${Y}_{eq}^{(0)}$ corresponds to ${n}_{eq}^{(0)}$ and in the non-relativistic limit is given by

where *g*_{
s
}counts the particle's spin degrees of freedom.

In the regime *x* <*x*_{
f
}, *Y* >> ${Y}_{eq}^{(0)}$ the equation (3.4) can be written as

Applying (3.12) at the freezing point *x*_{
f
}and using (3.10) and (3.11), leads, after a straightforward calculation, to the determination of *x*_{
f
}= *T*_{
f
}/${m}_{\tilde{\chi}}$ through

As usual, all quantities are expressed in terms of the dimensionless *x* ≡ *T*/${m}_{\tilde{\chi}}$ and *x*_{
in
}corresponds to the time *t*_{0} discussed previously, taken to represent the exit from the inflationary period of the Universe.

The first term on the right-hand-side of (3.13) is that of a conventional Cosmology for, say, an LSP carrying *g*_{
s
}spin degrees of freedom, playing the rôle of the dominant Cold Dark Matter species in a concrete and physically promising example [16–19], which we use in this work. The quantity $\u27e8$*vσ*$\u27e9$_{
f
}is the thermal average of *vσ* at *x*_{
f
}and *g*_{*} is *g*_{
eff
}of conventional Cosmology at the freeze-out point. The same notation holds for ${\tilde{g}}_{\ast}$. In our treatment above, we chose in (3.13) to present *x*_{
f
}in such a way so as to separate the conventional contributions, which reside in the first term, from the contributions of the dilaton and the non-critical-string dynamics, which are contained within the last two terms. The latter induce a shift in the freeze-out temperature. The penultimate term on the right hand side of (3.13), due to its logarithmic nature, does not affect much the freeze-out temperature. The last term, on the other hand, is more important and, depending on its sign, may shift the freeze-out point to earlier or later times. To quantify the amount of the shift one must solve the equations (2.11) of ref. [1]. We shall do this in a forthcoming publication [28], where we shall also present a detailed analysis of the effects of the non-critical string and dilaton-source terms on the constraints on supersymmetric particle physics models, extending conventional cosmology works [16–19].

For our purposes here we note that, in order to calculate the relic abundance, we must solve (3.12) from *x*_{
f
}to today's value *x*_{0}, corresponding to a temperature *T*_{0} ≈ 2.7°*K*. Following the usual approximations we arrive at the result:

In conventional Cosmology [21]${\tilde{g}}_{\ast}$ is replaced by *g*_{*} and the last term in (3.14) is absent. The quantity *J* is $J\equiv {\displaystyle {\int}_{{x}_{0}}^{{x}_{f}}\u3008v\sigma \u3009}dx$. By replacing *Y*(*x*_{
f
}) by its equilibrium value (3.10) the ratio of the first term on the r.h.s. of (3.14) to the second is found to be exactly the same as in the no-dilaton case. Therefore, by the same token as in conventional Cosmology, the first term can be safely omitted, as long as *x*_{
f
}is of order of 1/10 or less. Furthermore, the integral on the r.h.s. of (3.14) can be simplified if one uses the fact that $\u27e8$*vσ*$\u27e9$*n* is small as compared with the expansion rate $\dot{a}$/*a* after decoupling. For the purposes of the evaluation of this integral, therefore, this term can be omitted in (3.12), as long as we stay within the decoupling regime, and one obtains:

By integration this yields *Y*(*x*) = *Y*(*x*_{0}) $\mathrm{exp}(-{\displaystyle {\int}_{{x}_{0}}^{x}\Gamma {H}^{-1}dx/x})$. Using this inside the integral in (3.14) we get

where the function *ψ*(*x*) is given by $\psi (x)\equiv x\mathrm{exp}(-{\displaystyle {\int}_{{x}_{0}}^{x}\Gamma {H}^{-1}}dx/x)$. With the exception of the prefactor on the r.h.s. of (3.16), this is identical in form to the result derived in standard treatments, if ${\tilde{g}}_{\ast}$ is replaced by *g*_{*} and the value of *x*_{
f
}, implicitly involved in the integral *J*, is replaced by its value found in ordinary treatments in which the dilaton-dynamics and non-critical-string effects are absent.

The matter density of species $\tilde{\chi}$ is then given by

where the prefactor *f* is:

It is important to recall that the thermal degrees of freedom are counted by *g*_{
eff
}(c.f. (3.6)), and not ${\tilde{g}}_{eff}$, the latter being merely a convenient device connecting the total energy, thermal and non-thermal, to the temperature *T* (c.f. (3.5)). Hence,

In deriving (3.18) only the thermal content of the Universe is used, while the dilaton and the non-critical terms do not participate. Therefore the $\tilde{\chi}$'s matter density is given by

This formula tacitly assumes that the $\tilde{\chi}$s decoupled before neutrinos. For the relic abundance, then, we derive the following approximate result

where we used the fact that: $1+{\displaystyle {\int}_{{x}_{0}}^{{x}_{f}}\frac{\Gamma {H}^{-1}}{\psi (x)}}dx=\mathrm{exp}\left({\displaystyle {\int}_{{x}_{0}}^{{x}_{f}}\frac{\Gamma {H}^{-1}}{x}dx}\right)$. In (3.20), the quantity referred to as *no-source* is the well known no-source expression

where $J\equiv {\displaystyle {\int}_{{x}_{0}}^{{x}_{f}}\u3008v\sigma \u3009}dx$. However, as already remarked, the end point *x*_{
f
}in the integration is the shifted freeze-out point as determined by Eq. (3.13). The merit of casting the relic density in such a form is that it clearly exhibits the effect of the presence of the source. Certainly, if an accurate result is required, one can proceed without approximations and handle the problem numerically as in the standard treatments. We shall present such a more complete analysis in a forthcoming publication [28].

## 4 Discussion: Source Effects on Particle Phenomenology

From the expression (3.20) above, it becomes evident that the effects of the integral involving the source Γ on the relic density evolution may be quite important. Indeed, if Γ is kept negative at all times, this results in reduction of the relic density with time, contrary to what happens in the case where Γ is positive. In the former case, predictions for supersymmetric models [16–19] can be drastically altered, since the parameter space is enlarged, leaving more room for supersymmetry, probably beyond the reach of LHC, even in the case of constrained minimal supersymmetric standard models with compact parameter spaces of the embedding minimal supergravity theory. The opposite happens in the case of positive Γ, where the parameter space is shrunk and predictions can be very restrictive to almost excluding supersymmetry, especially if the prefactor turns out to be a large number.

In order to get a rough picture of the importance of such changes in the calculation of the relic density, let us assume that $\dot{a}$/*a* + $\dot{\Phi}$ = 0, as would be the case in the present era if matter were absent [25]. The non-critical terms would also contribute little in this case and from (3.1) one would have Γ ≈ $\dot{\Phi}$. Then, as follows from (3.1), the density of the dark matter species *n* would scale (with the scale factor) like radiation, *a*^{-4} (which in our approach retains its conventional cosmology scaling behaviour, as a result of the coupling of the gauge fields to the dilaton in the low-energy string inspired effective (slightly off-shell) action, which assumes the standard form). Therefore, Γ*H*^{-1} ≈ -1 and the function *ψ*(*x*) would be *ψ*(*x*) ≈ *x*^{2}/*x*_{0}. As a result, the prefactor in (3.20) becomes *x*_{0}/*x*_{
f
}which is an enormously small number, of order ~10^{-13} or so. Such a small number would result practically to no cosmological constraints on supersymmetric models! However, if this were the case, an even more dramatic situation would have taken place: the above unconventional scaling of dark matter, which is the dominant species, could not be disentangled from radiation, and the universe would never become matter dominated, in the standard cosmology sense, thus affecting significantly structure formation.

However this situation is not realized in nature, since at the present era neither matter nor the non-critical term contributions are negligible [1]. A decent approach is to solve the cosmological equations (2.11) and thus obtain the function Γ(*t*) at all times, *t*, from today to the remote past. This would allow for precise information to be obtained on the value of the freezing point *x*_{
f
}and of the prefactor appearing in the relic density (3.20). Thus, a complete numerical treatment, along the lines presented in [1], is needed in order to tackle important phenomenological questions regarding the relic abundance of Dark Matter in this framework. Nevertheless, the smallness of the prefactor in the crude approximation we have employed above, gives us a sneak preview of the drastic changes one may be faced with, as a result of the non-critical string dynamics.

A detailed phenomenology programme along these lines has started becoming available in [28], where the effects of matter in the above framework have been properly taken into account. In what follows we shall describe briefly the most important of these results for completeness. For details we refer the reader to [28]. When ordinary matter ("dust"), with scaling *a*^{-3} is present, then, solving the dynamical equations (2.11), as in [1], results in a complicated time dependence of the dilaton field, which is such that $\dot{\Phi}$ + $\dot{a}$/*a* ≠ 0. This affects the scaling of the dark matter species, in such a way that, although dark matter components with "exotic" scaling are present in past eras, nevertheless the density of the latter can be made negligible during the present epoch and is such that the nucleosynthesis (occurring at MeV scales as usual) is also not affected. In such scenarios [28], the neutralino dark matter density can be diluted only by a factor of order $\mathcal{O}$(10), while the baryon energy density is diluted only by a factor of order $\mathcal{O}$(1), upon adjusting only one parameter, namely the equation of state of the neutralino [28], assuming it to be the dominant dark matter species. This implies that dark matter enhances its odds, and future colliders, such as the LHC, have more available supersymmetry parameter space to exploit!

It must be noted at this stage, however, that the above conclusions are rather phenomenological, in the sense of being derived by exploiting the effective dynamical equations of the low energy four-dimensional world, assuming generic matter and radiation densities, in the spirit of standard cosmology analyses. In realistic string models, the results depend crucially on many details of the underlying microscopic model, such as the compactification procedure, the details of supersymmetry breaking, the geometry of the bulk space in brane models, string loop corrections *etc*. Regarding the latter effects, we also mention that at late eras of the universe, the various scaling components of the energy density, such as *a*^{-4}, *a*^{-3}, *etc*, involve mixed contributions from *both* matter and gravitational/dilaton and radiation parts, in such a way that one may even have negative dust contributions during the present era in some models [20], e.g. as a result of Kaluza-Klein graviton modes in some brane cosmology scenarios. Such unconventional energy budget can still be made compatible at least with supernovae data [20]. It remains to be seen of course what constraints the CMB and structure formation data imply on such cosmologies, but it is evident from the above that string cosmology, and in particular supercritical string one, may lead to fascinating new results which may hide further surprises around the corner.

Moreover, in such a framework one could also tackle other important issues, such as a possible resolution of the gravitino overproduction problem in effective supergravity inflationary models [29–31]. In the non-critical string framework discussed above, the gravitino, as a member of the gravitational supermultiplet of the string, will feel the non-critical and dilaton dissipation terms, while matter effects are unimportant for this issue. Thus, its density may be substantially reduced, in a way similar to that discussed above, as compared with conventional supergravity scenarios [28]. *affair à suivre...*

## References

- 1.
Diamandis GA, Georgalas BC, Lahanas AB, Mavromatos NE, Nanopoulos DV: Phys Lett B. 2006, 642: 179-10.1016/j.physletb.2006.09.035. [arXiv:hep-th/0605181]

- 2.
Ellis JR, Mavromatos NE, Nanopoulos DV: Phys Lett B. 1992, 293: 37-10.1016/0370-2693(92)91478-R. [arXiv:hep-th/9207103]

- 3.
Ellis JR, Mavromatos NE, Nanopoulos DV: A microscopic Liouville arrow of time, Invited review for the special Issue of J. Chaos Solitons Fractals. Edited by: Castro C, El Naschie MS. 1999, Elsevier Science, Pergamon, 10: 345-363. [arXiv:hep-th/9805120]

- 4.
David F: Mod Phys Lett A. 1988, 3: 1651-10.1142/S0217732388001975.

- 5.
Distler J, Kawai H: Nucl Phys B. 1989, 321: 509-10.1016/0550-3213(89)90354-4.

- 6.
Mavromatos NE, Miramontes JL: Mod Phys Lett A. 1989, 4: 1847-10.1142/S0217732389002082.

- 7.
E D'Hoker, Kurzepa PS: Mod Phys Lett A. 1990, 5: 1411-10.1142/S0217732390001608.

- 8.
Antoniadis I, Bachas C, Ellis JR, Nanopoulos DV: Phys Lett B. 1988, 211: 393-10.1016/0370-2693(88)91882-5.

- 9.
Antoniadis I, Bachas C, Ellis JR, Nanopoulos DV: Nucl Phys B. 1989, 328: 117-10.1016/0550-3213(89)90095-3.

- 10.
Antoniadis I, Bachas C, Ellis JR, Nanopoulos DV: Phys Lett B. 1991, 257: 278-10.1016/0370-2693(91)91893-Z.

- 11.
Ellis JR, Mavromatos NE, Nanopoulos DV, Sakharov A: New J Phys. 2004, 6: 171-10.1088/1367-2630/6/1/171. [arXiv:gr-qc/0407089]

- 12.
Ellis JR, Mavromatos NE, Nanopoulos DV, Westmuckett M: Int J Mod Phys A. 2006, 21: 1379-10.1142/S0217751X06028990. [arXiv:gr-qc/0508105], and references therein

- 13.
Gravanis E, Mavromatos NE: Phys Lett B. 2002, 547: 117-10.1016/S0370-2693(02)02778-8. [arXiv:hep-th/0205298]

- 14.
Mavromatos NE: arXiv:hep-ph/0309221

- 15.
Spergel DN, et al: Wilkinson Microwave Anisotropy Probe (WMAP) three year results: Implications for cosmology. arXiv:astro-ph/0603449

- 16.
Lahanas AB, Nanopoulos DV, Spanos VC: Phys Lett B. 2001, 518: 94-10.1016/S0370-2693(01)01034-6. [arXiv:hep-ph/0107151]

- 17.
Ellis JR, Olive KA, Santoso Y, Spanos VC: Phys Lett B. 2003, 565: 176-10.1016/S0370-2693(03)00765-2. [arXiv:hep-ph/0303043]

- 18.
Chattopadhyay U, Corsetti A, Nath P: Phys Rev D. 2003, 68: 035005-10.1103/PhysRevD.68.035005. [arXiv:hep-ph/0303201]

- 19.
Lahanas AB, Mavromatos NE, Nanopoulos DV: Int J Mod Phys D. 2003, 12: 1529-10.1142/S0218271803004286. [arXiv:hep-ph/0308251], and references therein

- 20.
Ellis JR, Mavromatos NE, Mitsou VA, Nanopoulos DV: Astropart Phys. 2007, 27: 185-10.1016/j.astropartphys.2006.10.007. [arXiv:astro-ph/0604272]

- 21.
Kolb EW, Turner MS: The Early Universe. Front Phys. 1990, Redwood City, USA: Addison-Wesley, 69:

- 22.
Shore GM: Nucl Phys B. 1987, 286: 349-10.1016/0550-3213(87)90445-7.

- 23.
Osborn H: Nucl Phys B. 1991, 363: 486-10.1016/0550-3213(91)80030-P.

- 24.
Curci G, Paffuti G: Nucl Phys B. 1987, 286: 399-10.1016/0550-3213(87)90447-0.

- 25.
Diamandis GA, Georgalas BC, Mavromatos NE, Papantonopoulos E, Pappa I: Int J Mod Phys A. 2002, 17: 2241-10.1142/S0217751X02010534. [arXiv:hep-th/0107124]

- 26.
Diamandis GA, Georgalas BC, Mavromatos NE, Papantonopoulos E: Int J Mod Phys A. 2002, 17: 4567-10.1142/S0217751X02010893. [arXiv:hep-th/0203241]

- 27.
Gasperini M, Piazza F, Veneziano G: Phys Rev D. 2002, 65: 023508-10.1103/PhysRevD.65.023508. [arXiv:gr-qc/0108016]

- 28.
Lahanas AB, Mavromatos NE, Nanopoulos DV: Phys Lett B. 2007, 649: 83-10.1016/j.physletb.2007.03.058. [arXiv:hep-ph/0612152] and in progress

- 29.
Ellis JR, Linde AD, Nanopoulos DV: Phys Lett B. 1982, 118: 59-10.1016/0370-2693(82)90601-3.

- 30.
Ellis JR, Kim JE, Nanopoulos DV: Phys Lett B. 1984, 145: 181-10.1016/0370-2693(84)90334-4.

- 31.
Ellis JR, Nanopoulos DV, Sarkar S: Nucl Phys B. 1985, 259: 175-10.1016/0550-3213(85)90306-2.

## Acknowledgements

The work of A.B.L. and N.E.M. is partially co-funded by the European Social Fund and National (Greek) Resources (EPEAEK II) PYTHAGORAS and by the European Union through the Marie Curie Research and Training Network *UniverseNet* (MRTN-CT-2006-035863). The work of D.V.N. is supported by D.O.E. grant DE-FG03-95-ER-40917.

## Author information

### Affiliations

## Rights and permissions

## About this article

#### Received

#### Accepted

#### Published

#### DOI

### Keywords

- Dark Matter
- Boltzmann Equation
- Light Supersymmetric Particle
- Einstein Frame
- Liouville Operator