Journal of Astrophysics

Volume 2013 (2013), Article ID 543717, 14 pages

http://dx.doi.org/10.1155/2013/543717

## Dark Energy from the Gas of Wormholes

Dubna International University of Nature, Society and Man, Universitetskaya Street 19, Dubna 141980, Russia

Received 21 March 2013; Accepted 29 May 2013

Academic Editors: C. Ghag, M. Jamil, and S. Profumo

Copyright © 2013 A. A. Kirillov and E. P. Savelova. 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.

#### Abstract

We assume the space-time foam picture in which the vacuum is filled with a gas of virtual wormholes. It is shown that virtual wormholes form a finite (of the Planckian order) value of the energy density of zero-point fluctuations. However such a huge value is compensated by the contribution of virtual wormholes to the mean curvature and the observed value of the cosmological constant is close to zero. A nonvanishing value appears due to the polarization of vacuum in external classical fields. In the early Universe some virtual wormholes may form actual ones. We show that in the case of actual wormholes vacuum polarization effects are negligible while their contribution to the mean curvature is apt to form the observed dark energy phenomenon. Using the contribution of wormholes to dark matter and dark energy we find estimates for characteristic parameters of the gas of wormholes.

#### 1. Introduction

As is well known modern astrophysics (and, even more generally, theoretical physics) faces two key problems. Those are the nature of dark matter and dark energy. Recall that more than 90% of matter of the Universe has a nonbaryonic dark (to say, mysterious) form, while lab experiments still show no evidence for the existence of such matter. Both dark components are intrinsically incorporated in the most successful CDM (Lambda cold dark matter) model which reproduces correctly properties of the Universe at very large scales (e.g., see [1] and references therein). We point out that CDM predicts also the presence of cusps in centers of galaxies [2] and a too large number of galaxy satellites. Therefore other models are proposed, for example, like axions [3], which may avoid these. To be successful such models should involve a periodic self-interaction and therefore require a fine tuning, while in general the presence of standard nonbaryonic particles cannot solve the problem of cusps. Indeed, if we admit the existence of a self-interaction in the dark matter component, or some coupling to baryonic matter (which should be sufficiently strong to remove cusps), then we completely change properties of the dark matter component at the moment of recombination and destroy all successful predictions at very large scales. Recall that both warm and self-interacting dark matter candidates are rejected by the observing spectrum [1]. In other words, the two key observational phenomena (cores of dark matter in centers of galaxies [4–6] and spectrum) give a very narrow gap for dark matter particles which seems to require attracting some exotic objects in addition to standard nonbaryonic particles.

As it was demonstrated recently [7] the problem of cusps can be cured, if some part of nonbaryon particles is replaced by wormholes. Wormholes represent extremely heavy (in comparison to particles) objects which at very large scales behave exactly like nonbaryon cold particles, while at smaller scales (in galaxies) they strongly interact with baryons and form the observed [4–6] cored distribution. We note that stable wormholes violate necessarily the averaged null energy conditions which gives the basic argument against the existence of such objects. Without exotic matter stable wormholes may however exist in modified theories, for example, see [8] and references therein. In the case when the energy conditions hold a wormhole collapses into a couple of conjugated (of equal masses) blackholes which almost impossible to distinguish from standard primordial blackholes. However, the topological nontriviality of such objects retains and gravitational effects of a gas of wormholes considered in [9] and some results of [7] still remain valid which means that nontraversable wormholes can be used to smooth cusps in centers of galaxies. Thus, it worth expecting that wormholes may play an important role in the explanation of the dark matter phenomenon.

Saving the dark matter component CDM requires the presence (~70%) of dark energy (of the cosmological constant). Moreover, there is evidence for the start of an acceleration phase in the evolution of the Universe [10–14]. In the present paper we use virtual wormholes to estimate the contribution of zero-point fluctuations in the value of the cosmological constant. The idea to relate virtual wormholes (or baby universes) and the cosmological constant is not new, in somewhat different context it was used by Coleman in [15] and developed in [16]. Our basic aim is to demonstrate that virtual wormholes form a finite value of the energy density of zero-point fluctuations.

It is necessary to point here out to the principle difference between actual and virtual wormholes. The principle difference is that a virtual wormhole exists only for a very small period of time and at very small scales and does not necessarily obey to the Einstein equations. It represents tunnelling event and therefore, the averaged null energy condition (ANEC) cannot forbid the origin of such an object. For the future we also note that a set of virtual wormholes may work as an actual wormhole opening thus the way for an artificial construction of wormhole-type objects in lab experiments.

In the present paper we describe a virtual wormhole as follows. From the very beginning we use the Euclidean approach (e.g., see [17] and the standard textbooks [18]). Then the simplest virtual wormhole is described by the metric where and is the step function. Such a wormhole has vanishing throat length, while the step function at the junction may cause a problem in Einstein’s equation or when a topological Euler term is involved. A more careful analysis needs to consider distributional curvature and so forth, see [19]. To avoid these difficulties we may consider from the very beginning a wormhole of a finite throat length ~ where the step function is replaced with a smooth function (e.g., ). Then where it is necessary one may consider the limit only in final expressions. This insures that the Bianchi identity holds and that the above metric remains inside of the domain of usual gravity. (Unexpectedly our approach (the use of step function) seems to irritate an essential part of physicists working with wormholes (e.g., PRD reviewers, moreover there is a claim that we study some other exotic objects which are not wormholes [20]). Therefore, it is necessary to clarify our position here. We are quite aware that at present state the physics of wormholes is certainly on the most speculative side. Therefore, we do not see the difference between specific forms of a smooth metric that one may use. For example, the simplest choice gives the well-known metric which in 3-dimensions is called as the Bronnikov-Ellis metric, or the metric which includes an additional parameter (an arbitrary length of the handle) which is not called by any name but is not less trivial. In the case of actual wormholes the exact and correct form of metric may be established only upon understanding the nature and properties of exotic matter which may support such a metric as a stable solution to the Einstein equations. For this time has not come yet. In the case of virtual wormholes even this is not important; for in the complete theory one has to sum over all possible metrics which formally may be included in in (5). Our basic aim is to present sufficiently clear and simple model (let it be far from realistic) which retains basic qualitative features related to a nontrivial topology.)

In the region , and the metric (2) is flat, while the region , with the obvious transformation , is also flat for . Therefore, the regions and represent two Euclidean spaces glued at the surface of a sphere with the centre at the origin and radius . Such a space can be described with the ordinary double-valued flat metric in the region by where the coordinates describe two different sheets of space. Now, identifying the inner and outer regions of the sphere allows the construction of a wormhole which connects regions in the same space (instead of two independent spaces). This is achieved by gluing the two spaces in (3) by motions of the Euclidean space (the Poincare motions). If is the position of the sphere in coordinates , then the gluing is the rule where , which represents the composition of a translation and a rotation of the Euclidean space (Lorentz transformation). In terms of common coordinates such a wormhole represents the standard flat space in which the two spheres (with centers at positions ) are glued by the rule (4). We point out that the physical region is the outer region of the two spheres. Thus, in general, the wormhole is described by a set of parameters: the throat radius , positions of throats , and rotation matrix .

In the present paper we assume the space-time foam picture in which the vacuum is filled with a gas of virtual wormholes. We show that virtual wormholes form a finite (of the Planckian order) value of the energy density of zero-point fluctuations. However such a huge value is compensated by the contribution of virtual wormholes to the mean curvature and the observed value of the cosmological constant should be close to zero.

To achieve our aim we, in Section 2, present the construction of the generating functional in quantum field theory. The main idea is that the partition function includes the sum over field configurations and the sum over topologies. Where the sum over topologies is the sum over virtual wormholes described above. Such an approach gives a rather good leading approximation for calculation of the partition function and corresponds to the standard methods (e.g., Ritza method, etc.). In Section 3 we investigate properties of the two-point Green function. We show that the presence of the gas of virtual wormholes can be described by the topological bias exactly as it happens in the presence of actual wormholes [7, 9]. For limiting topologies when the density of virtual wormholes becomes infinite the Green function shows a good ultraviolet behavior which means that there exists a class of such systems when quantum field theories are free of divergencies. We demonstrate how the sum over topologies defines the mean value for the bias which takes the sense of a cutoff function in the space of modes. In Section 4 we explicitly demonstrate that for a particular set of virtual wormholes the bias defines not more than the projection operator on the subspace of functions obeying to the proper boundary conditions at wormhole throats. The projective nature of the bias means that wormholes merely cut some portion of degrees of freedom (modes). Phenomenologically it means that wormholes can be described by the presence of ghost fields which compensate the extra (cut by wormholes) modes. In Section 5 we show how the cutoff expresses via some dynamic parameters of wormholes. The exact definition of such parameters we leave it for the future investigation. In Section 6 we consider the origin of the cosmological constant. We demonstrate that the cosmological constant is determined by the contribution of the energy density of zero-point fluctuations and by the contribution of virtual wormholes to the mean curvature. We estimate contribution of virtual wormholes to the mean curvature and show also that wormholes lead to a finite (of the planckian order) value of which requires considering the contribution from the smaller and smaller wormholes with divergent density . We also present arguments of why in the absence of external classical fields the total value of the cosmological constant is exactly zero, while it acquires a nonvanishing value due to vacuum polarization effects (i.e., due to an additional distribution of virtual wormholes) in external fields. We also speculate the possibility of the formation of actual wormholes and in Section 7 we estimate their contribution to the dark energy. Finally in Section 8 we repeat basic results an discuss some perspectives.

#### 2. Generating Function

The basic aim of this section is to construct the generating functional which can be used to get all possible correlation functions. Consider the partition function which includes the sum over topologies and the sum over field configurations For the sake of simplicity we use from the very beginning the Euclidean approach. The action has the form and we use the notions . If we fix the topology of space by placing a set of wormholes with parameters , then the sum over field configurations gives the well-known result where is the standard expression and is the Green function for a fixed topology, that is, for a fixed set of wormholes .

Consider now the sum over topologies . To this end we restrict the sum over the number of wormholes and integrals over parameters of wormholes: where and is the density of wormholes in the configuration space . We also point out that in general the integration over parameters is not free (e.g., it obeys the obvious restriction ). This defines the generating function as The sum over topologies assumes an additional averaging out for all mean values with the measure , where which obey the obvious normalization condition . The averaging out over topologies assumes the two stages. First we fix the total number of wormholes and average over the parameters of wormholes (i.e., over parameters of a static gas of wormholes in ). Then we sum over the number of wormholes (the so-called big canonical ensemble).

The basic difficulty of the standard field theory is that the perturbation scheme based on (7) leads to divergent expressions. This remains true for every particular topology of space (i.e., for any particular finite set of wormholes), since there always exists a scale below which the space looks like the ordinary Euclidean space. What we expect is that the sum over all possible topologies will remove such a difficulty.

And indeed, the above measure (11) has the structure where is the cosmological constant related to the energy density of zero-point fluctuations calculated for a particular distribution of wormholes. (We recall that the total cosmological constant should include also the contribution from the mean curvature (55).) Any finite distribution of wormholes leads to the divergent expression and is suppressed (i.e., ). However, the sum over all possible topologies assumes also the limiting topologies , where is the density of wormholes. In this limit wormhole throats degenerate into points and the minimal scale below which the space looks like the Euclidean space is merely absent. We point out that from the rigorous mathematical standpoint such limiting topologies cannot be described in terms of smooth manifolds, since they are not locally Euclidean and does not possess a finite set of maps. In mathematics similar objects are well known, for example, fractal sets. However, if a fractal set is obtained by cutting (by means of a specific rule or iterations) portions of space, our limiting topologies are obtained by gluing (identifying) some portions (or in the limit couples of points) of the Euclidean space. The basic feature of such topologies is that QFT becomes finite on such a set. Indeed, as we shall see a particular infinite distribution of wormholes can always be chosen in such a way that the energy of zero-point fluctuations becomes a finite (e.g., see the next section or the second term in (67)). In the sum over topologies only such limiting topologies do survive (i.e., ).

#### 3. The Two-Point Green Function

From (7) we see that the very basic role in QFT plays the two point Green function. Such a Green function can be found from the equation with proper boundary conditions at wormholes, which gives . Now let us introduce the bias function as where is the ballistic (or the standard Euclidean Green function) and the bias can be presented as where are fictitious sources at positions which should be added to obey the proper boundary conditions. We point out that the bias can be explicitly expressed via parameters of wormholes; that is, . For the sake of illustration we consider first a particular example.

##### 3.1. The Bias for a Particular Distribution of Wormholes (Rarefied Gas Approximation)

Consider now the bias for a particular set of wormholes. For the sake of simplicity we consider the case when . The Green function obeys the Laplace equation with proper boundary conditions at throats (we require and to be continual at throats). The Green function for the Euclidean space is merely (and for the Fourier transform). In the presence of a single wormhole which connects two Euclidean spaces this equation admits the exact solution. For outer region of the throat the source generates a set of multipoles placed in the center of sphere which gives the corrections to the Green function in the form (we suppose the center of the sphere at the origin) where and are four-dimensional spherical harmonics, for example, see [21]. In the present section we shall consider a dilute gas approximation and, therefore, it is sufficient to retain the lowest (monopole) term only. A single wormhole which connects two regions in the same space is a couple of conjugated spheres of the radius with a distance between centers of spheres. So the parameters of the wormhole are (The additional parameter (rotation matrix ) is important only for multipoles of higher orders.) . The interior of the spheres is removed and surfaces are glued together. Then the proper boundary conditions (the actual topology) can be accounted for by adding the bias of the source In the approximation (e.g., see for details [9]) the bias for a single wormhole takes the form This form for the bias is convenient when constructing the true Green function and considering the long-wave limit; however it is not acceptable in considering the short-wave behavior and vacuum polarization effects. Indeed, the positions of additional sources are in the physically nonadmissible region of space (the interior of spheres ). To account for the finite value of the throat size we should replace in (19) the point-like source with the surface density (induced on the throat); that is, Such a replacement does not change the value of the true Green function; however, now all extra sources are in the physically admissible region of space.

In the rarefied gas approximation the total bias is additive; that is, where is given by (9). For a homogeneous and isotropic distribution , and then for the bias we find Consider the Fourier transform and using the integral we find for the expression

*(1) Example of a Finite Density of Wormholes*. Consider now a particular (of a finite density) distribution of wormholes , for example,
where is the density of wormholes. In the case this function corresponds to a single wormhole with the throat size and the distance between throats . We recall that the action (6) remains invariant under translations and rotations which straightforwardly leads to the above function. Then reduces to . Thus from (23) we find
And for the true Green function we get

In the short-wave limit (, ) and therefore . This means that at very small scales the space filled with a finite density of wormholes looks like the ordinary Euclidean space. In the long-wave limit we get which gives , while in a more general case we find , where is the density of wormholes with a particular values of and , and for the bias function (15) we get

In other words, in the long-wave limit (, ) the presence of a particular set of virtual wormholes diminishes merely the value of the charge values.

*(2) Limiting Topologies or Infinite Densities of Wormholes.* Consider now the limiting distribution when the density of wormholes . Since every throat cuts the finite portion of the volume , this case requires considering the limit . We assume that in this limit , where is a finite specific distribution. Then (23) reduces to where and the bias (15), (18) becomes
It is important that this limit still agrees with the rarefied gas approximation, for the basic gas parameter (i.e., the portion of volume cut by wormholes) tends to . The above expression is obtained in the linear approximation only. Taking into account next orders (e.g., see [22] where the case of a dense gas is also considered) we find and the true Green function in a gas of wormholes becomes
In the long-wave limit the function can be expanded as which also defines a renormalization of charge values which coincides with (27). However in this limiting case we have some freedom in the choice of , which we can use to assign an arbitrary function of . In other words we may get here a class of limiting topologies where Green functions have a good ultraviolet behavior and quantum field theories in such spaces turn out to be finite. In particular, this will certainly be the case when
This possibility however requires the further and more deep investigation.

##### 3.2. Green Function, General Consideration

The action (6) remains invariant under translations with an arbitrary which means that the measure (11) does not actually depend on the position of the center of mass of the gas of wormholes and, therefore, we may restrict ourself with homogeneous distributions of wormholes in space only. Indeed, we may define , while the integration over gives the volume of ; that is, which disappears from (11) due to the denominator. (Technically, we may first restrict a portion of in (6) to a finite volume and then in final expressions consider the limit (which represents the standard tool in thermodynamics and QFT).) In what follows we shall omit the prime from .

Let us consider the Fourier representation which in the case of a homogeneous distribution of wormholes gives ; then we find and the Green function can be taken as Then for the total partition function we find where and comes from the integration measure (i.e., from the Jacobian of transformation from to ) We point out that can be changed by means of adding to the action (6) of an arbitrary “nondynamical” constant term which depends only on topology (e.g., a topological Euler term). The multiplier defines the simplest measure for topologies. Now by means of using the expression (32) and (33) we find the two-point Green function in the form where is the cutoff function (the mean bias) which is given by (We recall that in this integral contribute only limiting topologies in which density of wormholes diverges ())

At the present stage we still cannot evaluate the exact form for the cutoff function in virtue of the ambiguity of pointed out. Such a term may include two parts. First part describes the proper dynamics of wormholes and should be considered separately. Indeed, in general wormholes are dynamical self-gravitating objects which require considering the gravitational contribution to the action. Some part of such a contribution (mean curvature induced by wormholes) is discussed in Section 6. However, since a wormhole represents an extended nonlocal object, it possesses a rather complex dynamics and this problem requires the further investigation. The second part may describe “external conditions” (e.g., an external classical field in (33)) for the mean topology. Actually the last term can be used to prescribe an arbitrary particular value for the cutoff function . Indeed, the “external conditions” can be accounted for by adding the term , where plays the role of a specific chemical potential which implicitly depends on through the equation From (33) we see that the role of such a chemical potential may play the external current or equivalently an external classical field . In quantum field theory such a term leads merely to a renormalization of the cosmological constant. By other words the mean topology (i.e., the cutoff function or mean distribution of wormholes) is driven by the cosmological constant and vice versa.

#### 4. Topological Bias as a Projection Operator

By the construction the topological bias plays the role of a projection operator onto the space of functions (a subspace of functions on ) which obey the proper boundary conditions at throats of wormholes. This means that for any particular topology (for a set of wormholes) there exists the basis in which it takes the diagonal form with eigenvalues (since ). In this section we illustrate this simple fact (which is probably not obvious for readers) by the explicit construction of the reference system for a single wormhole when physical functions become (due to the boundary conditions) periodic functions of one of coordinates.

Indeed, consider a single wormhole with parameters (i.e., , where is the throat radius and are positions of throats in space. (In general, there exists an additional parameter which defines a rotation of one of throats before gluing. However, it does not change the subsequent construction. There always exists a diffeomorphic map of coordinates which sets such a matrix to unity).) Consider now a particular solution to the equation (harmonic function) for in the presence of the wormhole, which corresponds to the situation when throats possess a unit charge/mass but those have the opposite signs. Now define the family of lines of force which obey the equation with initial conditions . Physically, such lines correspond to lines of force for a two charged particles in positions with charges ±1. We note that all points which lay on the trajectory may be taken as initial conditions and they define the same line of force with the obvious redefinition . By other words we may take as a new coordinates the parameter and portion of the coordinates orthogonal to the family of lines . Coordinates can be taken as laying in the hyperplane which is orthogonal to the vector and goes through the point . (Instead of the construction used here one may use also another method. Indeed, consider two point charges, and then the function can be taken as a new coordinate. Wormhole appears when we identify (glue) surfaces . We point out that such surfaces are not spheres, though they reduce to spheres in the limit or .)

Let be the values of the parameter at which the line intersects the throats . Then instead of we may consider a new parameter as , so that when the parameter takes the values , respectively. The gluing procedure at throats means merely that we identify points at and and all physical functions in the space with a single wormhole become periodic functions of . Thus, the coordinate transformation gives the map of the above space onto the cylinder with a specific metric (where ) whose components are also periodic in terms of . Now we can continue the coordinates to the whole space (to construct a cover of the fundamental region ) simply admitting all values this, however, requires to introduce the bias since every point and every source in the fundamental region acquires a countable set of images in the nonphysical region (inside of wormhole throats). Considering now the Fourier transforms for we find We point out that the above bias gives the unit operator in the space of periodic functions of . From the standpoint of all possible functions on it represents the projection operator (taking an arbitrary function we find that upon the projection becomes a periodic function of ; that is, only periodic functions survive).

The above construction can be easily generalized to the presence of a set of wormholes. In the approximation of a dilute gas of wormholes we may neglect the influence of wormholes on each other (at least there always exists a sufficiently smooth map which transforms the family of lines of force for “independent” wormholes onto the actual lines). Then the total bias (projection) may be considered as the product where is the bias for a single wormhole with parameters . Every such a particular bias realizes projection on a subspace of functions which are periodic with respect to a particular coordinate , while the total bias gives the projection onto the intersection of such particular subspaces (functions which are periodic with respect to every parameter ).

#### 5. Cutoff

The projective nature of the bias operator allows us to express the cutoff function via dynamic parameters of wormholes. Indeed, consider a box in and periodic boundary conditions which gives (in final expressions we consider the limit , which gives ). And let us consider the decomposition for the integration measure in (33) as where includes also the contribution from . We point out that this measure plays the role of the action for the bias . Indeed, the variation of the above expression gives the equation of motions for the bias in the form which can be found by considering the proper dynamics of wormholes. We however do not consider the problem of the dynamic description of wormholes here and leave this for the future research. Moreover, we may expect that in the first approximation one may retain the linear term only. Indeed, this takes place when is not a dynamic variable, or if we take into account that is a collective variable. Then the projective nature of the bias means that it can be phenomenologically expressed via some Fermionic ghost field (e.g., ) where the negative and positive frequency parts of the operator obey the anticommutation relations . In the absence of ghost particles we get ; that is, and wormholes are absent. In the term of the ghost field the action becomes . Therefore in the leading approximation equations of motion take the linear form .

Thus taken into account that () we find

The simplest choice gives merely , where the sum is taken over the number of fields and is given by . In the case of a set of massless fields we find where and is the effective number of degrees of freedom (the number of boson minus fermion fields). To ensure the absence of divergencies one has to consider the number of fields [23]. However, such a choice gives the simplest estimate which, in general, cannot be correct. Indeed, while its behavior at very small scales (i.e., when exceeding the Planckian scales and ) may be physically accepted, since it produces some kind of a cutoff, on the mass-shell it gives the behavior which is merely incorrect (e.g., from (27) we see that the true behavior should be ).

One may expect that the true cutoff function has a much more complex behavior. Indeed, some theoretical models in particle physics (e.g., string theory) have the property to be lower-dimensional at very small scales. The mean cutoff gives the natural tool to describe a scale-dependant dimensional reduction [24, 25]. In fact, this function defines the spectral number of modes in the interval between and as Hence we can define the effective spectral dimension of space as follows: From the empirical standpoint the dimension is verified at laboratory scales only, while the rigorous tool to define the spectral density of states (or the mean cutoff) can give the lattice quantum gravity, for example, see [26, 27] and references therein. And indeed, the spectral dimension for nonperturbative quantum gravity defined via Euclidean dynamical triangulations was calculated recently in [28]. It turns out that it runs from a value of at short distance to at large distance scales. We also point out that all observed dark matter phenomena can be explained by the fractal dimension starting from scales , for example, [29–33].

#### 6. Cosmological Constant

Let us consider the total Euclidean action [17] The variation of the above action leads to the Einstein equations where is the stress energy tensor and is a naked cosmological constant. In cosmology such equations are considered from the classical standpoint, which means that they involve characteristic scales . However, the presence of virtual wormholes at Planckian scales defines some additional contribution in both parts of these equations which can be adsorbed into the cosmological constant. Therefore the total cosmological constant can be defined as where is the energy of zero-point fluctuations. (It includes also the contribution of zero-point fluctuations of gravitons.) That is, the mean vacuum value (we recall that in the standard QFT is infinite, while wormholes form a finite value) and is a contribution of wormholes into the mean curvature due to gluing (1).

##### 6.1. Contribution of Virtual Wormholes into Mean Curvature

Consider a single wormhole whose metric is given by (1) ; then the components of the Ricci tensor are (We are much obliged to the Referee who pointed out to the subtleties when working with a step function in the metric.) where is the unite normal vector to the throat surface, , , and is the number of dimensions. The curvature scalar is Then substituting in the above equation we find where is a smooth function which only in a limit becomes a step function, , and . In the limit we find and the last two terms in (51) are negligible as compared to the first two terms, while in four dimensions the last two terms vanish. Then the curvature is concentrated on the throat of the wormhole where differs from the step function. In the limit of a vanishing throat size we have at the throat and we get In general all the terms together in (51) and analogous terms in the Ricci tensor provide that the Bianchi identity holds and the energy is conserved [19]. Therefore none of them can be dropped out. However it is easy to see the leading contribution to the integral over space comes from the first term only and in the limit of a vanishing throat size we find for a single wormhole

In the case of a set of wormholes (24) we find where is the density of wormhole throats with a fixed value of the throat size . This defines the contribution to the cosmological constant from the mean curvature as We see that this quantity is always negative.

##### 6.2. Stress Energy Tensor

In this section we consider the contribution from matter fields. In the case of a scalar field the stress energy tensor has the form Then the mean vacuum value of the stress energy tensor can be obtained directly from the two-point green function (32), (35) as By means of using the Fourier transform and the expressions (32), (35) we arrive at where the property has been used and is the cutoff function (36).

For the sake of simplicity we consider the massless case. Then by the use of the cutoff from the previous section we get the finite estimate ( is the effective number of the field helicity states) where the sum is taken over the number of fields. Since the leading contribution comes here from very small scales, we may hope that this value will not essentially change if the true cutoff function changes the behavior on the mass-shell as (e.g., if we take with as ).

To understand how wormholes remove divergencies, it will be convenient to split the bias function into two parts , where corresponds to the standard Euclidean contribution, while is the contribution of wormholes. The first part gives the well-known divergent contribution of vacuum field fluctuations with , while the second part remains finite for any finite number of wormholes and, due to the projective nature of the bias described in the previous section, it partially compensates (reduces) the value of the cosmological constant; that is, , where and is a negative finite contribution of a finite set of wormholes.

Consider now the particular distribution of virtual wormholes (24) and evaluate their contribution to the cosmological constant which is given by . Then from the expressions (22) and (24) we get which gives where For (we recall that by the construction ) this function has the value . Thus, for the contribution of wormholes we find

##### 6.3. Vacuum Value of the Cosmological Constant

From the above expression we see that to get the finite value of the cosmological constant one should consider the limit (infinite density of virtual wormholes) which requires considering the smaller and smaller wormholes. From the other hand we have the obvious restriction , where is the volume of one throat (wormholes cannot cut more than the total volume of space). (We also point out that in removing divergencies the leading role plays the zero-point energy. Indeed , while the mean curvature has the order and for we have . Moreover in the limit , one gets and therefore .) Therefore, in the leading order it seems to be sufficient to retain point-like wormholes only (i.e., consider the limit ). Then instead of (24) we may assume the vacuum distribution of virtual wormholes in the form where and has the meaning of the infinite density of point-like wormholes, while remains a finite. In this case the volume cut by wormholes vanishes and the rarefied gas approximation works well. This defines the bias and the mean cutoff (here we define the Fourier transform ) as The contribution to the mean curvature (55) can be expressed via the same function as Thus, for the total cosmological constant we get the expression We stress that all these terms should be finite. Indeed all distributions of virtual wormholes which lead to an infinite value of are suppressed in (5) by the factor ~, while the minimal value is reached when wormholes cut all of the volume of space and the action is merely . (Frankly speaking this statement is not rigorous. At first look the two last terms in (67) are independent and one may try to take an arbitrary big. If this were the case then the action would not possess the minimum at all. However cannot be arbitrary big, since it will violate the rarefied gas approximation and the linear expression (67) brakes down. Moreover, fermions give here a contribution of the opposite sign. The rigorous investigation of this problem requires the further studying and we present it elsewhere.) We also point out that here we considered the real scalar field as the matter source, while in the general case the stress energy tensor should include all existing Bose and Fermi fields (Fermi fields give a negative contribution to ).

The value of looks like a free parameter, which in quantum gravity runs with scales [27]. However at large scales its asymptotic value may be uniquely fixed by the simple arguments as follows. Indeed in quantum field theory properties of the ground state (vacuum) change when we imply an external classical fields. The same is true for the distribution of virtual wormholes , for example, see (33) and (37) and, therefore, which we describe in the next subsection. We recall that in gravity the role of the external current plays the stress energy tensor of matter fields . However one believes that in the absence of all classical fields the vacuum state should represent the most symmetric (Lorentz invariant) state which in our case corresponds to the Euclidean space. In order to be consistent with the Einstein equations this requires which uniquely fixes the value of in (67). We point out that from somewhat different considerations such a choice was advocated earlier in [15, 34, 35].

##### 6.4. Vacuum Polarization in an External Field

Consider now topology fluctuations in the presence of an external current. In the presence of an external current the distribution of virtual wormholes changes . Indeed in (33) for the case of a weak external field the contribution of the external current into the action can be expanded as , where Then using (36) we find , where is the bias related to an additional distribution of virtual wormholes and defines the dispersion of vacuum topology fluctuations (here ). The exact definition of requires the further development of a fundamental theory. In particular, it can be numerically calculated in lattice quantum gravity [28]. However it can be shown that at scales it reduces to and therefore Now comparing this function with (23) we relate the additional distribution of virtual wormholes and the external classical field as where . We point out that the above expression does not define the value which requires an additional consideration. Moreover, in general the external field does not possess a symmetry and therefore the correction does not reduce to a single cosmological constant. However, such corrections always violate the averaged null energy condition [36, 37] and may be considered as some kind of dark energy or, by other words, it represents an exotic matter. Some portion of dark energy still has the form of the cosmological constant which defines a nonvanishing present day value (we recall that fermions give a contribution of the opposite sign) where denotes an averaging over rotations.

The only unknown parameter in (72) is the dispersion which defines the intensity of topology fluctuations in the vacuum. It has also the sense of the efficiency coefficient which defines the portion of the energy of the external field spent on the formation of additional wormholes. Though the evaluation of requires the further development of a fundamental theory, one may expect that , where is the mean cutoff. It is expected that as and . This means that as and as , while it takes the maximum value at Planckian scales . By other words, the most efficient transmission of the energy into wormholes takes place for wormholes of the Planckian size. In the case when external classical fields have characteristic scales in (72) the efficiency coefficient and the cutoff become constant , , while their ratio may be estimated as , where is the effective number of fundamental fields which contribute to and , are dark energy and baryon energy densities, respectively. According to the modern picture this ratio gives , while (as ) and therefore this defines the estimate for the effective number of fundamental fields (helicity states) as .

##### 6.5. Speculations on the Formation of Actual Wormholes

As we already pointed out the additional distribution of virtual wormholes (72) reflects the symmetry of external classical fields and therefore it forms a homogeneous and isotropic background and perturbations. We recall that virtual wormholes represent an exotic form of matter. In the early Universe such perturbations start to develop and may form actual wormholes. The rigorous description of such a process represents an extremely complex and interesting problem which requires the further study. Some aspects of the behavior of the exotic density perturbations were considered in [7], while the simplest example of the formation of a wormhole-type object was discussed recently by us in [38]. Therefore we may expect that some portion of such a form of dark energy is reserved now in actual wormholes which we consider in the next section.

#### 7. Dark Energy from Actual Wormholes

Consider now the contribution to the dark energy from the gas of actual wormholes. Unlike the virtual wormholes, actual wormholes do exist at all times and, therefore, a single wormhole can be viewed as a couple of conjugated cylinders . So that the number of parameters of an actual wormhole is less , where is the radius of and is a spatial part of .

Actual wormholes also produce two kinds of contribution to the dark energy. One comes from their contribution to the mean curvature which corresponds to an exotic stress energy momentum tensor. Such a stress energy momentum tensor reflects the dark energy reserved by additional virtual wormholes discussed in the previous section. Such energy is necessary to support actual wormholes as a solution to the Einstein equations. The second part comes from vacuum polarization effects by actual wormholes. The consideration in the previous section shows that for macroscopic wormholes the second part has the order and is negligible as compared to the curvature (since macroscopic wormholes have throats ). However, for the sake of completeness and for methodological aims we describe it as well.

For rigorous evaluation of dark energy of the second type we, first, have to find the bias analogous to (19) for the topology . There are many papers treating different wormholes in this respect (e.g., see [36, 37] and references therein). However, in the present paper for an estimation we shall use a more simple trick.

##### 7.1. Beads of Virtual Wormholes (Quantum Wormhole)

Indeed, instead of the cylinders we consider a couple of chains (beads of virtual wormholes ). Such an idea was first suggested in [39] and one may call such an object as quantum wormhole. Then the bias can be written straightforwardly where with and is the step. We may expect that upon averaging over the position the bias for the beads will reproduce the bias for cylinders (at least it looks like a very good approximation). We point out that the averaging out and the sum reduces to a single integral of the zero term in (74). And moreover, the resulting total bias corresponds merely to a specific choice of the distribution function in (21). Namely, we may take where and is the distribution of cylinders, which can be taken as ( is 3-dimensional density) Using the normalization condition we find the relation , where is a 4-dimensional density of wormholes and is the effective number of wormholes on the unit length of the cylinder (i.e., the frequency with which the virtual wormhole appears at the positions ). This frequency is uniquely fixed by the requirement that the volume which cuts the bead is equal to that which cuts the cylinder (i.e., and ). Thus, we can use directly expression (23) and find (compare to (25)) where . Here the first term merely coincides with that in (25) and, therefore, it gives the contribution to the cosmological constant , while the second term describes a correction which does not reduce to the cosmological constant and requires a separate consideration.

##### 7.2. Stress Energy Tensor

From (57) we find that the stress energy tensor reduces to the two functions where and and these functions are By means of the use of the spherical coordinates , , and we get where , , and is given by For we find Thus, finally we find which upon the continuation to the Minkowsky space gives the equation of state in the form (An arbitrary gas of wormholes splits in fractions with a fixed and .) which in the case when behaves like a cosmological constant. However when such a constant is extremely small and can be neglected, while the leading contribution comes from the mean curvature.

##### 7.3. Mean Curvature

In this subsection we consider the Minkowsky space. Then the simplest actual wormhole can be described by the metric analogous to (1), for example, see [7] where . To avoid problems with the Bianchi identity and the conservation of energy the step function should be also smoothed as in (1). The stress energy tensor which produces such a wormhole can be found from the Einstein equation . Both regions and represent portions of the ordinary flat Minkowsky space and therefore the curvature is . However on the boundary it has the singularity. Since the metric (86) does not depend on time we find where is the outer normal to the throat , and are additional terms (e.g., see (49) and (51)) which in the leading order are negligible upon averaging over some portion of space . In the case of a set of wormholes this gives in the leading order where is the radius of a throat and is the position of the center of the throat in space and . In the case of a homogeneous and isotropic distribution of such throats we find (averaging over spatial directions gives ) where where stands for the trace of the stress energy momentum tensor which one has to add to the Einstein equations to support such a wormhole. It is clear that such a source violates the weak energy condition and, therefore, it reproduces the form of dark energy (i.e., ). If the density of such sources (and, resp., the density of wormholes) is sufficiently high, then this results in the observed [10–14] acceleration of the scale factor for the Friedmann space as with , for example, see also [40–42]. In terms of the 4-dimensional density of wormholes we get as and, therefore, the leading contribution indeed comes from the mean curvature.

#### 8. Estimates and Concluding Remarks

Now consider the simplest estimates. Actual wormholes seem to be responsible for the dark matter [7, 9]. Therefore, to get the estimate to the number density of wormholes is rather straightforward. First wormholes appear at scales when dark matter effects start to display themselves; that is, at scales of the order , which gives in that range the number density The characteristic size of throats can be estimated from (89) . Since the density of dark energy is , where is the critical density, then we find the estimate where is the Solar radius, /(sec Mpc)) and is the Hubble constant. We also recall that the background density of baryons generates a nonvanishing wormhole rest mass (where is the scale factor of the Universe and therefore remains constant), for example, see [7]. It produces the dark matter density related to the wormholes as . The typical mass of a wormhole is estimated as where is the Solar mass. We point out that this mass has not the direct relation to the parameters of the gas of wormholes. However it defines the moment when wormhole throats separated from the cosmological expansion. The above estimate shows that if wormholes form due to the development of perturbations in the exotic matter, then this process should start much earlier than the formation of galaxies.

Thus, we see that virtual wormholes should indeed lead to the regularization of all divergencies in QFT which agrees with recent results [28]. Therefore, they form the local finite value of the cosmological constant. In the absence of external classical fields such a value should be exactly zero at macroscopic scales. A some nonvanishing value for the cosmological constant appears as the result of vacuum polarization effects in external fields. Indeed, external fields form an additional distribution of virtual wormholes which possess an exotic stress energy tensor (some kind of dark energy). Only some part of it forms the cosmological constant, while the rest reflects the symmetry of external fields and possesses inhomogeneities. We assume that during the evolution of our Universe inhomogeneities in the exotic matter develop and may form actual wormholes. Although this problem requires the further and more deep investigation we refer to [38] where the formation of a simplest wormhole-like object has been considered. In other words, such polarization energy is reserved now in a gas of actual wormholes. We estimated parameters of such a gas and believe that such a gas may indeed be responsible for both, dark matter and dark energy phenomena.

#### References

- J. R. Primack, “The nature of dark matter,” http://arxiv.org/abs/astro-ph/0112255.
- J. Diemand, M. Zemp, B. Moore, J. Stadel, and M. Carollo, “Cusps in cold dark matter haloes,”
*Monthly Notices of the Royal Astronomical Society*, vol. 364, no. 2, pp. 665–673, 2005. View at Publisher · View at Google Scholar - E. W. Mielke and J. A. V. Vélez Pérez, “Axion condensate as a model for dark matter halos,”
*Physics Letters B*, vol. 671, no. 1, pp. 174–178, 2009. View at Google Scholar - G. Gentile, P. Salucci, U. Klein, D. Vergani, and P. Kalberla, “The cored distribution of dark matter in spiral galaxies,”
*Monthly Notices of the Royal Astronomical Society*, vol. 351, no. 3, pp. 903–922, 2004. View at Publisher · View at Google Scholar - D. T. F. Weldrake, W. J. G. de Blok, and F. Walter, “A high-resolution rotation curve of NGC 6822: a test-case for cold dark matter,”
*Monthly Notices of the Royal Astronomical Society*, vol. 340, no. 1, pp. 12–28, 2003. View at Publisher · View at Google Scholar - W. J. G. de Blok and A. Bosma, “High-resolution rotation curves of low surface brightness galaxies,”
*Astronomy and Astrophysics*, vol. 385, pp. 816–846, 2002. View at Google Scholar - A. Kirillov and E. P. Savelova, “Density perturbations in a gas of wormholes,”
*Monthly Notices of the Royal Astronomical Society*, vol. 412, no. 3, pp. 1710–1720, 2011. View at Publisher · View at Google Scholar · View at Scopus - P. Kanti, B. Kleihaus, and J. Kunz, “Stable Lorentzian wormholes in dilatonic Einstein-Gauss-Bonnet theory,”
*Physical Review D*, vol. 85, no. 4, Article ID 044007, 2012. View at Publisher · View at Google Scholar · View at Scopus - A. A. Kirillov and E. P. Savelova, “Dark matter from a gas of wormholes,”
*Physics Letters B*, vol. 660, no. 3, pp. 93–99, 2008. View at Publisher · View at Google Scholar · View at Scopus - A. G. Riess, A. V. Filippenko, P. Challis et al., “Observational evidence from supernovae for an accelerating universe and a cosmological constant,”
*The Astronomical Journal*, vol. 116, no. 3, p. 1009, 1998. View at Publisher · View at Google Scholar - S. Perlmutter, G. Aldering, G. Goldhaber et al., “Measurements of Ω and Λ from 42 high-redshift supernovae,”
*The Astrophysical Journal*, vol. 517, no. 2, p. 565, 1999. View at Publisher · View at Google Scholar - J. L. Tonry, B. P. Schmidt, B. Barris et al., “Cosmological results from high-z supernovae,”
*The Astrophysical Journal*, vol. 594, no. 1, p. 594, 2003. View at Publisher · View at Google Scholar - N. W. Halverson, E. M. Leitch, C. Pryke et al., “Degree angular scale interferometer first results: a measurement of the cosmic microwave background angular power spectrum,”
*The Astrophysical Journal*, vol. 568, no. 1, p. 38, 2002. View at Publisher · View at Google Scholar - C. B. Netterfield, P. A. R. Ade, J. J. Bock et al., “A measurement by BOOMERANG of multiple peaks in the angular power spectrum of the cosmic microwave background,”
*The Astrophysical Journal*, vol. 571, no. 2, p. 604, 2002. View at Publisher · View at Google Scholar - S. Coleman, “Why there is nothing rather than something: a theory of the cosmological constant,”
*Nuclear Physics B*, vol. 310, no. 3-4, pp. 643–668, 1988. View at Google Scholar · View at Scopus - I. Klebanov, L. Susskind, and T. Banks, “Wormholes and the cosmological constant,”
*Nuclear Physics B*, vol. 317, no. 3, pp. 665–692, 1989. View at Google Scholar · View at Scopus - S. W. Hawking, “Spacetime foam,”
*Nuclear Physics B*, vol. 114, p. 349, 1978. View at Google Scholar - T. P. Cheng and L. F. Li,
*Gauge Theory of Elementary Particles*, Clarendon Press, Oxford, UK, 1984. - A. H. Taub, “Space-times with distribution valued curvature tensors,”
*Journal of Mathematical Physics*, vol. 21, no. 6, pp. 1423–1431, 1980. View at Google Scholar · View at Scopus - A. Shatsky, privite communications.
- V. A. Fock, “Zur Theorie des Wasserstoffatoms,”
*Zeitschrift für Physik*, vol. 98, no. 3-4, pp. 145–154, 1935. View at Publisher · View at Google Scholar - E. P. Savelova, “Gas of wormholes in Euclidean quantum field theory,” http://arxiv.org/abs/1211.5106.
- A. A. Kirillov and E. P. Savelova, “On the behavior of bosonic systems in the presence of topology fluctuations,” http://arxiv.org/abs/0808.2478.
- A. A. Kirillov, “Violation of the Pauli principle and dimension of the Universe at very large distances,”
*Physics Letters B*, vol. 555, no. 1-2, pp. 13–21, 2003. View at Google Scholar - A. A. Kirillov, “Dark matter, dark charge, and the fractal structure of the Universe,”
*Physics Letters B*, vol. 535, no. 1–4, pp. 22–24, 2002. View at Publisher · View at Google Scholar - J. Ambjorn, J. Jurkiewicz, and R. Loll, “The spectral dimension of the universe is scale dependent,”
*Physical Review Letters*, vol. 95, no. 17, Article ID 171301, 4 pages, 2005. View at Publisher · View at Google Scholar - E. Manrique, S. Rechenberger, and F. Saueressig, “Asymptotically safe Lorentzian gravity,”
*Physical Review Letters*, vol. 106, no. 25, Article ID 251302, 4 pages, 2011. View at Publisher · View at Google Scholar · View at Scopus - J. Laiho and D. Coumbe, “Evidence for asymptotic safety from lattice quantum gravity,”
*Physical Review Letters*, vol. 107, no. 16, Article ID 161301, 2011. View at Publisher · View at Google Scholar · View at Scopus - A. A. Kirillov and D. Turaev, “On modification of the Newton's law of gravity at very large distances,”
*Physics Letters B*, vol. 532, no. 3-4, pp. 185–192, 2002. View at Publisher · View at Google Scholar - A. A. Kirillov and D. Turaev, “The universal rotation curve of spiral galaxies,”
*Monthly Notices of the Royal Astronomical Society*, vol. 371, no. 1, pp. L31–L35, 2006. View at Publisher · View at Google Scholar - A. A. Kirillov and D. Turaev, “Foam-like structure of the universe,”
*Physics Letters B*, vol. 656, pp. 1–8, 2007. View at Publisher · View at Google Scholar - A. A. Kirillov, “The nature of dark matter,”
*Physics Letters B*, vol. 632, no. 4, pp. 453–462, 2006. View at Publisher · View at Google Scholar - A. A. Kirillov and E. P. Savelova, “Astrophysical effects of space-time foam,”
*Gravitation and Cosmology*, vol. 14, no. 3, pp. 256–261, 2008. View at Google Scholar - S. W. Hawking and N. Turok, “Open inflation without false vacua,”
*Physics Letters B*, vol. 425, no. 1-2, pp. 25–32, 1998. View at Publisher · View at Google Scholar - Z. C. Wu, “The cosmological constant is probably zero, and a proof is possibly right,”
*Physics Letters B*, vol. 659, no. 5, pp. 891–893, 2008. View at Publisher · View at Google Scholar - R. Garattini, “Entropy and the cosmological constant: a spacetime-foam approach,”
*Nuclear Physics B*, vol. 88, no. 1–3, Article ID 991003, pp. 297–300, 2000. View at Google Scholar - M. Visser,
*Lorentzian Wormholes*, Springer, New York, NY, USA, 1996. - A. A. Kirillov and E. P. Savelova, “Artificial wormhole,” http://arxiv.org/abs/1204.0351.
- E. W. Mielke, “Knot wormholes in geometrodynamics?”
*General Relativity and Gravitation*, vol. 8, no. 3, pp. 175–196, 1977. View at Publisher · View at Google Scholar · View at Scopus - E. B. Gliner, “Algebraic properties of the energy-momentum tensor and vacuum-like states of matter,”
*Journal of Experimental and Theoretical Physics*, vol. 22, p. 378, 1966. View at Google Scholar - A. H. Guth, “Inflationary universe: a possible solution to the horizon and flatness problems,”
*Physical Review D*, vol. 23, no. 2, pp. 347–356, 1981. View at Publisher · View at Google Scholar - A. A. Linde, “A new inflationary universe scenario: a possible solution of the horizon, flatness, homogeneity, isotropy and primordial monopole problems,”
*Physics Letters B*, vol. 108, no. 6, pp. 389–393, 1982. View at Publisher · View at Google Scholar