# Wormholes and Flux Tubes in the 7D Gravity on the Principal Bundle with SU(2) Gauge Group as the Extra Dimensions

## I Introduction

The presence of the off-diagonal components of the multidimensional (MD) metric can essentially change the properties of the appropriate MD gravitational equations. The physical reason for this is evident: these components of the MD metric are connected with the physical gauge fields (electromagnetic or Yang-Mills gauge fields) and excluding the physical fields can lead to this essential change. Let us repeat the following theorem which is useful for understanding how many degrees of freedom we have. , :

Let $`G`$ be a structural group of the principal bundle. Then there is a one-to-one correspondence between the $`G`$-invariant metrics

$$ds^2=\phi (x^\alpha )\underset{a=5}{\overset{dimG}{}}\left[\sigma ^a+A_\mu ^a(x^\alpha )dx^\mu \right]^2+g_{\mu \nu }(x^\alpha )dx^\mu dx^\nu $$
(1)   
on the total space $`𝒳`$ and the triples $`(g_{\mu \nu },A_\mu ^a,\phi )`$. Here $`g_{\mu \nu }`$ is the 4D Einstein’s pseudo - Riemannian metric on the base; $`A_\mu ^a`$ are the gauge fields of the group $`G`$ ( the nondiagonal components of the multidimensional metric); $`\phi \gamma _{ab}`$ is the symmetric metric on the fibre ($`\sigma _a=\gamma _{ab}\sigma ^b,\gamma _{ab}=\delta _{ab}`$, $`a=5,\mathrm{},`$ dim $`G`$ is the index on the fibre and $`\mu =0,1,2,3`$ is the index on the base).

According to this theorem we have the following independent degrees of freedom: the scalar field $`\phi (x^\alpha )`$, the gauge fields $`A_\mu ^a(x^\alpha )`$ and the 4D metric $`g_{\mu \nu }(x^\alpha )`$. Note that all fields in this MD gravity can depend only on the spacetime points (points on the base) as the total space is $`G`$-invariant. Such kind of MD gravity can easily solve the following problem: why the physical degrees of freedom do not depend on the extra dimensions (ED). In addition to this, a topological structure of the ED is given that leads to a decrease of the numbers of equations connected with the ED in comparison with an ordinary MD gravity with non-fixed ED. Usually, the number of these equations is too large: this results in essential problems with the compactification. Therefore, it is necessary to introduce some external fields in order to have a compactified ED. In our case the vacuum MD gravitational equations is sufficient to obtain the solutions with the compactified ED. In MD gravity on the principal bundle the preferable choice of the coordinate transformations as the initial interpretation of Kaluza-Klein gravity is

$`y^a`$ $`=`$ $`y^a(y^a)+f^a\left(x^\mu \right),`$ (2)   
$`x^\mu `$ $`=`$ $`x^\mu \left(x^\mu \right).`$ (3)   
here $`y^a`$ are the coordinates on the fibre and $`x^\mu `$ are the coordinates on the base. The first term in (2) means that the choice of coordinate system on the fibre is arbitrary. The second term indicates that in addition we can move the origin of the coordinate system on each fibre on the value $`f^a(x^\mu )`$. It is well known that such a transformation law (2) leads to a local gauge transformation for the appropriate non-Abelian field (see for overview ). That means, the coordinate transformation (2) and (3) are the most natural transformations for the MD gravitation on the principal bundle. Certainly we can perform the much more general coordinate transformations:

$`y^a`$ $`=`$ $`y^a(y^a,x^\mu ),`$ (4)   
$`x^\mu `$ $`=`$ $`x^\mu (y^a,x^\mu ).`$ (5)   
But in this case we will mix the points of the fibre <sup>§</sup><sup>§</sup>§which are the elements of some group. with the points of the base which are the ordinary spacetimes points.. Then the new coordinates $`y^a`$ and $`x^\mu `$ are not the coordinates along the fibre and the base of the given bundle.

We can introduce the new coordinates and kill the $`A_\mu ^a`$. But in this case the initial 4D metric can be changed very radical. For example, the initial static spherically symmetric 4D metric can become nonstatic and nondiagonal. This situation is very clear in the initial Kaluza interpretation of 5D gravity with the constant and nonvariable $`G_{55}`$ component of the metric. In this case we have the ordinary vacuum electrogravity which is equivalent to the 5D gravity with the nondynamical $`G_{55}`$ component of MD metric. And we can choose the new coordinates $`x^B`$ so that the initial 5D metric:

$$ds^2=\left(d\chi +A_\mu dx^\mu \right)^2+g_{tt}dt^2+g_{rr}dr^2+r^2\left(d\theta ^2+\mathrm{sin}^2\theta d\varphi ^2\right)$$
(6)   
will be

$$ds^2=G_{55}d\chi _{}^{}{}_{}{}^{2}+g_{\mu \nu }dx^\mu dx^\nu $$
(7)   
Evidently metric (7) will be nonstatic, nondiagonal and depend on the 5<sup>th</sup> coordinate. Of course, the new coordinate system $`x_{}^{}{}_{}{}^{A}`$ ($`A=0,1,2,3,5`$) is worse than the initial $`x^A`$. This remark allow us to say that the MD metric with the off-diagonal components of the metric $`G_\mu ^a`$ can have unusual properties in comparison with the solutions where $`G_\mu ^a=0`$.

## II 7D ansatz and equations

Here we will consider the MD gravity on the principal bundle with SU(2) this is the gauge group of the weak interaction structural group. In this case the extra dimensions is SU(2) group <sup>\**</sup><sup>\**</sup>\**topologically it means that the ED are the $`S^3`$ sphere. and the 4D physical spacetime is the base of this bundle. We will search a solution for the following 7D metric

$$ds^2=\frac{\mathrm{\Sigma }^2(r)}{u^3(r)}dt^2dr^2a(r)\left(d\theta ^2+\mathrm{sin}^2\theta d\varphi ^2\right)r_0^2u(r)\left(\sigma ^a+A_\mu ^adx^\mu \right)^2$$
(8)   
here $`r_0`$ is some constant, $`\sigma ^a`$ ($`a=5,6,7`$) are the Maurer-Cartan form with relation $`d\sigma ^a=ϵ_{bc}^a\sigma ^b\sigma ^c`$

$`\sigma ^1`$ $`=`$ $`{\displaystyle \frac{1}{2}}(\mathrm{sin}\alpha d\beta \mathrm{sin}\beta \mathrm{cos}\alpha d\gamma ),`$ (9)   
$`\sigma ^2`$ $`=`$ $`{\displaystyle \frac{1}{2}}(\mathrm{cos}\alpha d\beta +\mathrm{sin}\beta \mathrm{sin}\alpha d\gamma ),`$ (10)   
$`\sigma ^3`$ $`=`$ $`{\displaystyle \frac{1}{2}}(d\alpha +\mathrm{cos}\beta d\gamma ),`$ (11)   
where $`0\beta \pi ,0\gamma 2\pi ,0\alpha 4\pi `$ are the Euler angles. We choose the potential $`A_\mu ^a`$ in the ordinary monopole-like form

$`A_\theta ^a`$ $`=`$ $`{\displaystyle \frac{1}{2}}(1f(r))\{\mathrm{sin}\varphi ;\mathrm{cos}\varphi ;0\},`$ (12)   
$`A_\varphi ^a`$ $`=`$ $`{\displaystyle \frac{1}{2}}(1f(r))\mathrm{sin}\theta \{\mathrm{cos}\varphi \mathrm{cos}\theta ;\mathrm{sin}\varphi \mathrm{cos}\theta ;\mathrm{sin}\theta \},`$ (13)   
$`A_t^a`$ $`=`$ $`v(r)\{\mathrm{sin}\theta \mathrm{cos}\varphi ;\mathrm{sin}\theta \mathrm{sin}\varphi ;\mathrm{cos}\theta \},`$ (14)   
Let us introduce the color electric $`E_i^a`$ and magnetic $`H_i^a`$ fields

$`E_i^a`$ $`=`$ $`F_{ti}^a,`$ (15)   
$`H_i^a`$ $`=`$ $`\sqrt{\gamma }ϵ_{ijk}\sqrt{g_{tt}}F^{ajk}`$ (16)   
here the field strength components are defined via $`F_{\mu \nu }^a=A_{\nu ,\mu }^aA_{\mu ,\nu }^a+ϵ_{bc}^aA_\mu ^bA_\nu ^c`$, $`\gamma `$ is the determinant of the 3D space matrix, ($`i,j=1,2,3`$) are the space index. In our case we have

$`E_rv^{},E_{\theta ,\varphi }vf,`$ (17)   
$`H_r{\displaystyle \frac{\mathrm{\Sigma }}{u^{3/2}}}{\displaystyle \frac{1f^2}{a}},H_{\theta ,\varphi }f^{}`$ (18)   
In order to have the wormhole-like (WH) solution we must demand that the functions $`\mathrm{\Sigma }(r),u(r),a(r),f(r)`$ are even functions and $`v(r)`$ is an odd function. This means that at origin ($`r=0`$) we have only the radial $`E_r`$ and $`H_r`$ fields that indicate the presence of a flux tube of color electric and magnetic fields across the throat of this WH-like solution. The substitution to the 7D gravitational equations <sup>††</sup><sup>††</sup>††The deduction of these gravitational equations for any group $`G`$ is given in the Appendix.

$`R_\mu ^{\overline{A}}`$ $`=`$ $`0`$ (19)   
$`R_{\overline{5}}^{\overline{5}}+R_{\overline{6}}^{\overline{6}}+R_{\overline{7}}^{\overline{7}}`$ $`=`$ $`0`$ (20)   
leads to the following system of equations

$`R_{\overline{2}\overline{2}}2R_{\overline{3}\overline{3}}+R{\displaystyle \frac{\mathrm{\Sigma }^{\prime \prime }}{\mathrm{\Sigma }}}+{\displaystyle \frac{a^{}\mathrm{\Sigma }^{}}{a\mathrm{\Sigma }}}{\displaystyle \frac{4}{r_0^2u}}{\displaystyle \frac{r_0^2u}{4a}}f^2{\displaystyle \frac{r_0^2u}{8a^2}}\left(f^21\right)^2`$ $`=`$ $`0`$ (21)   
$`R_{\overline{2}\overline{2}}{\displaystyle \frac{1}{2}}R24{\displaystyle \frac{\mathrm{\Sigma }^{}u^{}}{\mathrm{\Sigma }u}}24{\displaystyle \frac{u^2}{u^2}}+16{\displaystyle \frac{a^{}\mathrm{\Sigma }^{}}{a\mathrm{\Sigma }}}+4{\displaystyle \frac{a^2}{2a^2}}{\displaystyle \frac{16}{a}}+`$ (22)   
$`4{\displaystyle \frac{r_0^2u^4}{\mathrm{\Sigma }^2}}v^22{\displaystyle \frac{r_0^2u}{a}}f^28{\displaystyle \frac{r_0^2u^4}{a\mathrm{\Sigma }^2}}f^2v^2+{\displaystyle \frac{r_0^2u}{a^2}}\left(f^21\right)^2{\displaystyle \frac{48}{ur_0^2}}`$ $`=`$ $`0`$ (23)   
$`R_{\overline{3}\overline{3}}=R_{\overline{4}\overline{4}}{\displaystyle \frac{a^{\prime \prime }}{a}}+{\displaystyle \frac{a^{}\mathrm{\Sigma }^{}}{a\mathrm{\Sigma }}}{\displaystyle \frac{2}{a}}+{\displaystyle \frac{r_0^2u}{4a}}f^2{\displaystyle \frac{r_0^2u^4}{a\mathrm{\Sigma }^2}}f^2v^2+{\displaystyle \frac{r_0^2u}{4a^2}}\left(f^21\right)^2`$ $`=`$ $`0`$ (24)   
$`R_{\overline{a}}^{\overline{a}}{\displaystyle \frac{u^{\prime \prime }}{u}}+{\displaystyle \frac{u^{}\mathrm{\Sigma }^{}}{u\mathrm{\Sigma }}}{\displaystyle \frac{u^2}{u^2}}+{\displaystyle \frac{u^{}a^{}}{ua}}{\displaystyle \frac{4}{r_0^2u}}+`$ (25)   
$`{\displaystyle \frac{r_0^2u^4}{3\mathrm{\Sigma }^2}}v^2{\displaystyle \frac{r_0^2u}{6a}}f^2+{\displaystyle \frac{2r_0^2u^4}{3a\mathrm{\Sigma }^2}}f^2v^2{\displaystyle \frac{r_0^2u}{12a^2}}\left(f^21\right)^2`$ $`=`$ $`0`$ (26)   
$`R_{\overline{1}\overline{5}}v^{\prime \prime }+v^{}\left({\displaystyle \frac{\mathrm{\Sigma }^{}}{\mathrm{\Sigma }}}+4{\displaystyle \frac{u^{}}{u}}+{\displaystyle \frac{a^{}}{a}}\right){\displaystyle \frac{2}{a}}vf^2`$ $`=`$ $`0,`$ (27)   
$`R_{\overline{3}\overline{5}}f^{\prime \prime }+f^{}\left({\displaystyle \frac{\mathrm{\Sigma }^{}}{\mathrm{\Sigma }}}+4{\displaystyle \frac{u^{}}{u}}\right)+4{\displaystyle \frac{u^3}{\mathrm{\Sigma }^2}}fv^2{\displaystyle \frac{f}{a}}\left(f^21\right)`$ $`=`$ $`0`$ (28)   
Note that the equations (27) and (28) are the “Yang-Mills” equations. This system of ordinary differential equations is extremely difficult for the analytical investigation <sup>‡‡</sup><sup>‡‡</sup>‡‡although there exists a closed-form solution for the simplest case when $`a(r)=const`$ (flux tube solution ).. Therefore we will search for a numerical solution of this system.

## III Numerical Investigation

Now we must write the initial conditions. At the origin ($`r=0`$) we can expand all functions by this manner

$`a(x)`$ $`=`$ $`1+{\displaystyle \frac{a_2}{2}}x^2+\mathrm{},`$ (29)   
$`\mathrm{\Sigma }(x)`$ $`=`$ $`\mathrm{\Sigma }_0+{\displaystyle \frac{\mathrm{\Sigma }_2}{2}}x^2+\mathrm{},`$ (30)   
$`u(x)`$ $`=`$ $`u_0+{\displaystyle \frac{u_2}{2}}x^2+\mathrm{},`$ (31)   
$`v(x)`$ $`=`$ $`v_1x+{\displaystyle \frac{v_3}{6}}x^3+\mathrm{},`$ (32)   
$`f(x)`$ $`=`$ $`f_0+{\displaystyle \frac{f_2}{2}}x^2+\mathrm{},`$ (33)   
here we introduce the dimensionless coordinate $`x=r/\sqrt{a_0}`$ ($`a_0=a(0)`$) and redefine $`a(x)/a_0a(x)`$, $`\sqrt{a_0}v(x)v(x)`$, $`r_0^2/a_0r_0^2`$. Then we can rescale time and the constant $`r_0`$ so that $`\mathrm{\Sigma }_0=u_0=1`$. Thus we have only the following initial conditions for the numerical calculations

$`a_0=1,u_0=1,\mathrm{\Sigma }_0=1,v_0=0,f_0=f_0,`$ (34)   
$`a_0^{}=0,u_0^{}=0,\mathrm{\Sigma }_0^{}=0,v_0^{}=v_1,f_0^{}=0,`$ (35)   
We see that our system depends on two parameters only: $`f_0`$ and $`v_1`$ <sup>\**</sup><sup>\**</sup>\**of course the 7D metric depends on the three parameters: $`a_0`$, $`f_0`$ and $`v_1`$.. The constrained equation (23) for the initial data give us

$$r_0^2=\frac{1+\sqrt{1+3\left[v_1^2+\frac{1}{4}\left(f_0^21\right)^2\right]}}{\left[v_1^2+\frac{1}{4}\left(f_0^21\right)^2\right]}$$
(36)   
this equation is written in dimensionless variables. The numerical calculations of Eq.’s (21)-(28) are presented in the Figs. 1-5. The initial data for these calculations are the following: $`f_0=0.2,v_1=0.3,\mathrm{\hspace{0.33em}0.5},\mathrm{\hspace{0.33em}0.6},\mathrm{\hspace{0.33em}0.61},\mathrm{\hspace{0.33em}0.615},\mathrm{\hspace{0.33em}1.0},\mathrm{\hspace{0.33em}2.0}`$.

## IV Physical discussion

Figs. 1-5 show us that in accordance with some relation between $`f_0`$ and $`v_1`$ there are two types of solutions:

1. There exists some value of the radial coordinate $`r_1`$ that $`a(\pm r_1)=0`$, $`u(\pm r_1)=s(\pm r_1)=\mathrm{}`$. Probably this means that at $`r=\pm r_1`$ points we have a singularity. This 4D part of the MD metric we can name as gravitational flux tube. As we have some flux of color electric/magnetic field between the points $`r=r_1`$ and $`r=+r_1`$ where the color electric/magnetic charges are located.
2. There exists some value of the radial coordinate $`r_2`$ that $`a(\pm r_2),s(\pm r_2)<\mathrm{}`$, $`u(\pm r_2)=0`$. The value $`u(\pm r_2)=0`$ means that the interval $`ds^2=0`$ on the hypersurface $`r=\pm r_2,t,\theta ,\varphi =const`$. Since the value of $`a(\pm r_2)`$ is finite we can name this type of the solutions as the WH-like.

This type of the solutions is close to the 5D case investigated in . In this Ref. it was shown that the 5D vacuum Kaluza-Klein theory for the metric ($`r_0`$ is radius of the U(1) gauge group, $`Q=nr_0`$ is the magnetic charge, $`\chi `$ is the 5<sup>th</sup> coordinate)

$`ds^2`$ $`=`$ $`e^{2\nu (r)}dt^2r_0^2e^{2\psi (r)2\nu (r)}\left[d\chi +\omega (r)dt+n\mathrm{cos}\theta d\varphi \right]^2`$ (37)   
$``$ $`dr^2a(r)(d\theta ^2+\mathrm{sin}^2\theta d\varphi ^2),`$ (38)   
has the following solutions

1. E $`>`$ H <sup>\*†</sup><sup>\*†</sup>\*†$`E=q/a`$ and $`H=Q/a`$ are the electric and magnetic charges.. WH-like solution located between two $`ds^2=0`$ hypersurfaces.
2. E = H. Infinite flux tube with constant electric and magnetic fields.
3. E $`<`$ H. Finite flux tube located between two singularities at points $`r=\pm r_1`$ where are the electric and magnetic charges.

In the first 5D case the solution exists for $`|r|>r_0`$ where $`e^{2\nu }<0`$ and the metric is asymptotically flat. The whole construction can be interpreted as the Euclidean WH with the Lorentzian throat <sup>\*‡</sup><sup>\*‡</sup>\*‡The words Euclidean and Lorentzian can be exchanged. . The question is: whether this situation can be kept in the 7D case for the second type of the solutions? In Ref’s ,  a similar idea had been investigated about the changing of the signature of the 4D metric in some MD metric on the regular $`T`$-hole horizon.

Also in our case the numerical calculations show that on the $`(f_0,v_1)`$ plane there is a curve that separates regions with the different solution type. Evidently, in the first rough approximation the equation for this curve is

$$\frac{a_0^{\prime \prime }}{a_0}=\frac{2}{a_0}\frac{r_0^2u_0}{a_0^2}\left(f_0^21\right)^2=0$$
(39)   
In Fig. 6 these regions with the different type of solutions are shown.

In the $`(f_0,v_1)`$ plane we can single out a few cases allowing a more detailed analysis.

### A $`f_0=\pm 1,v_1=0`$

Immediately from the “Yang-Mills” equations (27), (28) we see that

$$v(r)=1,f(r)=\pm 1$$
(40)   
All terms with gauge fields in (21)-(26) vanish, and at the origin $`(r=0)`$ the equation for the initial data give us

$$\frac{1}{a_0}+\frac{3}{u_0r_0}=0$$
(41)   
As $`a(r),u(r)>0`$ this relation cannot be satisfied. This allows us to say that in the absence of the gauge fields (off-diagonal components of the MD metric) the gravitational flux tubes and WH-like solutions do not exist.

### B $`f=0`$

In this case (27) equation is easy to integrate

$$v^{}=\frac{q\mathrm{\Sigma }}{r_0^2au^4}$$
(42)   
And from (17), (18) we see that there are only the radial color component of the electric and magnetic fields.

It is very interesting to note that there is some very simple analytical solution in the case $`a(r)=const`$ . The solution is

$`a(r)`$ $`=`$ $`{\displaystyle \frac{2q^2}{7}}={\displaystyle \frac{r_0^2}{8}}=const`$ (43)   
$`\mathrm{\Sigma }(r)`$ $`=`$ $`\mathrm{cosh}\left({\displaystyle \frac{7r}{2\sqrt{2}q}}\right)`$ (44)   
$`v(r)`$ $`=`$ $`{\displaystyle \frac{\sqrt{2}}{r_0}}\mathrm{sinh}\left({\displaystyle \frac{7r}{2\sqrt{2}q}}\right)`$ (45)   
$`u(r)`$ $`=`$ $`1`$ (46)   
$`q`$ $`=`$ $`\sqrt{{\displaystyle \frac{7a}{2}}}`$ (47)   
We can apply the following gauge transformation to the potentials of Eqs. (12 \- 14) with $`f(r)=0`$

$$A_\mu ^{}=S^1A_\mu Si(_\mu S^1)S$$
(48)   
where

$$S=\left(\begin{array}{cc}\mathrm{cos}\frac{\theta }{2}& e^{i\varphi }\mathrm{sin}\frac{\theta }{2}\\ e^{i\varphi }\mathrm{sin}\frac{\theta }{2}& \mathrm{cos}\frac{\theta }{2}\end{array}\right)$$
(49)   
With this gauge transformation we find that the gauge potentials become

$`A_\theta ^a`$ $`=`$ $`(0;0;0)`$ (50)   
$`A_\phi ^a`$ $`=`$ $`(\mathrm{cos}\theta 1)(0;0;1)`$ (51)   
$`A_t^a`$ $`=`$ $`v(r)(0;0;1)`$ (52)   
In fact this means that the potential (50)-(52) is an Abelian one. This solution corresponds to the points $`x_1`$ and $`x_2`$ in Fig. 6 and we can call it an infinite gravitational flux tube. It is very close to the Levi-Civita-Robinson-Bertotti solution , ,  in 4D electrogravity. This remark poses a very interesting problem: how are the other solutions on the curve $`C_1`$ $`C_2`$ (Fig. 6) ? Are these solutions infinite flux tubes or something else Unfortunately, here the numerical investigation is more difficult because of a sensitive dependence on initial data. This fact becomes also clear by the following consideration: These lines are the borderlines between two different domains of attraction, and there this sensitive behaviour is quite typical. ?

For the question which other Lie groups can lead to similar results we can say that for any $`N>2`$ there is an inclusion SU$`(2)\mathrm{SU}(N)`$. This means that for any such SU$`(N)`$ we can use the ansatz (12)-(14). As one knows, also all other compact semi-simple Lie groups (besides U(1), of course), have SU(2) as one of their subgroups. Therefore, for all these Lie groups similar statements are also valid. For the other Lie groups it is complicated to find out the spherically symmetric ansatz for the off-diagonal components of the MD metric.

We would like to point out the following dimensional reduction expression for the Ricci scalar $`R(E)`$ on the total space of the principal bundle

$`{\displaystyle d^4xd^dy\sqrt{|detG_{AB}|}R(E)}=`$ (53)   
$`V_G{\displaystyle d^4x\sqrt{|g|}\phi ^{d/2}\left[R(M)+R(G)\frac{1}{4}\phi F_{\mu \nu }^aF_a^{\mu \nu }+\frac{1}{4}d(d1)_\mu \phi ^\mu \phi \right]}`$ (54)   
here $`R(M),R(G)`$ are the Ricci scalars of the base and structural group of the principal bundle respectively; $`V_G`$ is the volume of the group $`G`$, cf. , eq. (8.12). Immediately we see that for an arbitrary Lie group $`G`$ we face the problem of how to get the appropriate ansatz for the gauge field $`A_\mu ^a`$.

In this context it should be noted that also in the models discussed here, a generalized Birkhoff theorem (cf. ) is valid. This means that a spherically symmetric solution possesses a further isometry without additional assumptions. Of course, this result rests on the symmetries chosen for our model.

## V The possible physical applications

Probably, the most interesting above-mentioned solution in 7D gravity on the principal bundle with the SU(2) structural group is the WH-like solution. We remember that if the $`G_{55},G_{66},G_{77}`$ components of the MD metric are not the dynamical variables then this MD gravity is equivalent to the pure 4D gravity + SU(2) Yang-Mills theory. It can be supposed that in the Universe there exist regions where the $`G_{55},G_{66},G_{77}`$ metric components are nondynamical variables and that there exist other regions where these components are dynamical variables. The composite WH can be of such a kind. It means that we have our above-mentioned WH-like solution as a throat of composite WH and two 4D Yang-Mills black holes attached to this throat (see Fig. 7).

Such construction can polarize a space-time foam by the following way . Without presence of the 7D (5D) throat the handles of the space-time foam are located disordered in the space-time (Fig. 8a). But after the appearance of the throat of composite WH the location of these handles is ordered (Fig. 8b). Such model can be a geometrical model of the renormalization of the color/electric charge (the case for the electric charge is described in ). This composite WH with the polarization of the space-time foam is a continuation of Wheeler’s idea about “charge without charge” and “mass without mass” in the vacuum gravity.

In this connection we can recall Wheeler’s quote in Ref. : “We are therefore led to consider the view that the electron is nothing but a collective state of excitation of the foam-like medium $`\mathrm{}`$ In other words the electron is not a natural starting point for the description of nature, according to the present reinterpretation of the views of Lorentz. Instead it is a first order correction to vacuum physics. That vacuum, that zero order state of affairs, with its enormous concentrations of electromagnetic energy and multiply-connected topologies, has to be described properly before one has the starting point for a proper perturbation theoretic development.” May be all these words can be used to a neutrino with mass, i.e. neutrino is a wormhole in the polarized space-time foam. In these cases there remains to solve the problem of the geometric description of spin. A similar 6-dimensional model, called $`M`$-fluxbranes, has been discussed recently in .

It is interesting to note that the above-mentioned gravitational flux tube solutions cannot be the throat of the composite WH as they have a singularity on the place of the $`ds^2=0`$ hypersurface for the WH-like solution. It imposes some restriction on the possible relation between color electric and magnetic fields in the composite WH.

## Acknowledgements

We would like to thank A. Kirillov and B. Schmidt for useful comments. VD is grateful for financial support by a Georg Forster Research Fellowship from the Alexander von Humboldt Foundation. HJS thanks the DFG for financial support.

## A Gravitational equations for the gravity on the principal bundle

We consider the MD gravity on the principal bundle with the structural group $`G`$ . In this case the extra dimensions are the group $`G`$, and the 4D physical spacetime is the base of this bundle.

According to above-mentioned theorem the MD metric on the total space is

$$ds^2=\mathrm{\Sigma }_{\overline{A}}\mathrm{\Sigma }^{\overline{A}}$$
(A1)   
where

$`\mathrm{\Sigma }^{\overline{A}}`$ $`=`$ $`h_B^{\overline{A}}dx^B,`$ (A2)   
$`\mathrm{\Sigma }^{\overline{a}}`$ $`=`$ $`\phi (x^\alpha )\sigma ^{\overline{a}}+h_\mu ^{\overline{a}}(x^\alpha )dx^\mu ,`$ (A3)   
$`\sigma ^{\overline{a}}`$ $`=`$ $`h_b^{\overline{a}}dx^b,`$ (A4)   
$`\mathrm{\Sigma }^{\overline{\mu }}`$ $`=`$ $`h_\nu ^{\overline{\mu }}(x^\alpha )dx^\nu `$ (A5)   
here $`x^B`$ are the coordinates on the total space ($`B=0,1,2,3,5,\mathrm{},N`$, is the MD index dim $`(G)=N`$), $`x^a`$ is the coordinates on the group $`G`$ ($`a=5,\mathrm{},N`$), $`x^\mu =0,1,2,3`$ are the coordinates on the base of the bundle, $`\overline{A}`$ is the $`N`$-bein index, $`h_B^{\overline{A}}`$ is the $`N`$-bein, $`\sigma ^{\overline{a}}`$ are the 1-forms on the group $`G`$ satisfying $`d\sigma ^{\overline{a}}=f_{bc}^a\sigma ^{\overline{b}}\sigma ^{\overline{c}}`$, ($`f_{bc}^a`$) are the structural constants for the group $`G`$. We must note that the functions $`\phi ,h_\mu ^{\overline{a}},h_\nu ^{\overline{\mu }}`$ can depend only on the $`x^\mu `$ points on the base as the fibres of our bundle are locally homogeneous spaces. The matrix $`h_B^{\overline{A}}`$ has the following form

$$h_B^{\overline{A}}=\left(\begin{array}{cc}\phi h_b^{\overline{a}}& h_\mu ^{\overline{a}}\\ 0& h_\mu ^{\overline{\nu }}\end{array}\right)$$
(A6)   
The inverse matrix $`h_{\overline{A}}^B`$ is

$$h_{\overline{A}}^B=\left(\begin{array}{cc}\phi ^1h_{\overline{a}}^b& h_{\overline{\mu }}^b\\ 0& h_{\overline{\mu }}^\nu \end{array}\right)$$
(A7)   
here $`h_{\overline{\mu }}^b=\phi ^1h_{\overline{a}}^bh_\nu ^{\overline{a}}h_{\overline{\mu }}^\nu `$. Also we see that we have only the following degrees of freedom: $`\phi (x^\alpha ),h_{\overline{\mu }}^b(x^\alpha )`$ and $`h_{\overline{\mu }}^\nu (x^\alpha )`$, $`h_{\overline{b}}^a`$ is given and not varying. Varying with respect to $`h_{\overline{\mu }}^A=(h_{\overline{\mu }}^a,h_{\overline{\mu }}^\nu )`$ leads to the equations

$$R_A^{\overline{\mu }}\frac{1}{2}h_A^{\overline{\mu }}R=0$$
(A8)   
here $`\overline{A}=\overline{a}`$, $`\overline{\nu }`$, $`R_B^{\overline{A}}`$ is the MD Ricci tensor. Let $`x^a`$ be the coordinates on the group $`G`$ then

$$\phi \sigma ^{\overline{a}}=\phi h_b^{\overline{a}}dx^b$$
(A9)   
Varying with respect to $`\phi (x^\mu )`$ leads to the following result

$$\frac{\delta }{\delta \phi }\left(hR\right)=\frac{\delta \left(h_{\overline{a}}^b/\phi \right)}{\delta \phi }\frac{\delta \left(hR\right)}{\delta \left(h_{\overline{a}}^b/\phi \right)}=\frac{1}{\phi ^2}h_{\overline{a}}^b\left(R_b^{\overline{a}}\frac{1}{2}h_b^{\overline{a}}R\right)=0$$
(A10)   
here $`h=deth_B^{\overline{A}}`$, $`R`$ is the MD Ricci scalar for the metric on the total space. As $`h_{\overline{a}}^\nu =0`$ we can write

$$h_{\overline{a}}^b\left(R_b^{\overline{a}}\frac{1}{2}h_b^{\overline{a}}R\right)+h_{\overline{a}}^\nu \left(R_\nu ^{\overline{a}}\frac{1}{2}h_\nu ^{\overline{a}}R\right)=h_{\overline{a}}^A\left(R_A^{\overline{a}}\frac{1}{2}h_A^{\overline{a}}R\right)=0$$
(A11)   
From (A8) and (A11) we see that

$$h_{\overline{a}}^A\left(R_A^{\overline{a}}\frac{1}{2}h_A^{\overline{a}}R\right)+h_{\overline{\mu }}^A\left(R_A^{\overline{\mu }}\frac{1}{2}h_A^{\overline{\mu }}R\right)=h_{\overline{B}}^A\left(R_A^{\overline{B}}\frac{1}{2}h_A^{\overline{B}}R\right)=0.$$
(A12)   
This means that

$$R=0.$$
(A13)   
Hence from (A11) we can write

$$h_{\overline{a}}^AR_A^{\overline{a}}=R_{\overline{a}}^{\overline{a}}$$
(A14)   
Finally we have the following equation system for the MD gravity on the principal bundle

$`R_A^{\overline{\mu }}`$ $`=`$ $`0`$ (A15)   
$`R_{\overline{a}}^{\overline{a}}=R_{\overline{5}}^{\overline{5}}+\mathrm{}+R_{\overline{N}}^{\overline{N}}`$ $`=`$ $`0.`$ (A16)   
We note that the (A16) Eq. is an analog for the Brans-Dicke scalar gravity.

In addition we see that (A15) can be wrote as

$$R_{\overline{\mu }A}=0orh_{\overline{A}}^BR_{\overline{\mu }B}=R_{\overline{\mu }\overline{A}}=0$$
(A17)   
Figure captions

Fig. 1. Function $`a(x)`$. Initial data: $`f_0=0.2`$, for curve 1: $`v_1=0.3`$, for curve 2 : $`v_1=0.5`$, for curve 3: $`v_1=0.6`$, for curve 4: $`v_1=0.61`$, for curve 5: $`v_1=0.615`$, for curve 6: $`v_1=1.0`$, and for curve 7: $`v_1=2.0`$.

Fig. 2. Function $`u(x)`$. Initial data as for Fig. 1.

Fig. 3. Function $`\mathrm{\Sigma }(x)`$. Initial data as for Fig. 1.

Fig. 4. Function $`f(x)`$. Initial data as for Fig. 1.

Fig. 5. Function $`v(x)`$. Initial data as for Fig. 1.

Fig. 6. The curves C<sub>1</sub> and C<sub>2</sub> separate the regions with different types of the solutions. In the 2-regions with $`a_0^{\prime \prime }<0`$ we have the flux tube solutions and for the 1-regions $`a_0^{\prime \prime }>0`$ \- the wormhole-like solutions.

Fig. 7. The composite wormhole with the 1-throat as MD wormhole-like solution and two black 2-holes attached to its 3-ends.

Fig. 8. A polarization of space-time foam in the presence of MD insertion. Fig. 8a presents the unpolarized space-time foam. Fig. 8b presents the polarized space-time foam. 1 are the virtual wormholes, 2 is the multidimensional insertion.