On a remarkable electromagnetic field in the Einstein Universe

TL;DR

This paper constructs a time-dependent electromagnetic field solution in the Einstein universe, demonstrating a knotted, finite energy, radiating configuration that can be translated into flat spacetime, with implications for cosmological models.

Contribution

It introduces a novel, explicit electromagnetic solution in the Einstein universe and relates it to flat spacetime via conformal equivalence, expanding understanding of electromagnetic fields in curved and cosmological settings.

Findings

Derived a knotted electromagnetic solution in Einstein universe

Mapped the solution to flat spacetime using conformal equivalence

Analyzed electromagnetic fields in expanding Friedmann models

Abstract

We present a time-dependent solution of the Maxwell equations in the Einstein universe, whose electric and magnetic fields, as seen by the stationary observers, are aligned with the Clifford parallels of the $3$-sphere $S^3$. The conformal equivalence between Minkowski's spacetime and (a region of) the Einstein cylinder is then exploited in order to obtain a knotted, finite energy, radiating solution of the Maxwell equations in flat spacetime. We also discuss similar electromagnetic fields in expanding closed Friedmann models, and compute the matter content of such configurations.

Figures

Equations

$$\mathrm{d}{\theta }^i=-{\epsilon }_{ijk}{\theta }^j\wedge {\theta }^k,$$

$$X_0={\partial }_t,[X_i,X_j]=2{\epsilon }_{ijk}{X}_k.$$

$$dF=d\star F=0,$$

$$F=E^i{\theta }^i\wedge \mathrm{d}t+\frac{1}{2}{B}^i{\epsilon }_{ijk}{\theta }^j\wedge {\theta }^k.$$

$$\{\begin{array}{l}{X}_i(E^i)=X_i(B^i)=0,\\ \dot{B^i}-2E^i+{\epsilon }_{ijk}{X}_j(E^k)=0,\\ \dot{E^i}+2B^i-{\epsilon }_{ijk}{X}_j(B^k)=0.\end{array}$$

$$\{\begin{array}{l}\dot{B^i}=2E^i,\\ \dot{E^i}=-2B^i.\end{array}$$

$$E=E_0\cos (2t)X_1,B=E_0\sin (2t)X_1,$$

$${\tilde{g}}_{ab}={\Omega }^{-2}{g}_{ab},\Omega =2\cos \left(\frac{t-\psi }{2}\right)\cos \left(\frac{t+\psi }{2}\right),$$

$$\tilde{H}(t)=\underset{{\tilde{\Sigma }}_t}{\int }\tilde{T}({\tilde{X}}_0{\tilde{X}}_0)\mathrm{d}(\widetilde{vol}),$$

$${\tilde{T}}_{ab}={\tilde{F}}_{ac}{\tilde{F}}_b^c-\frac{1}{4}{\tilde{g}}_{ab}{\tilde{F}}_{cd}{\tilde{F}}^{cd}$$

$$\tilde{H}(t)={\int }_{{\Sigma }_t}\Omega T(X_0,X_0)\mathrm{d}(vol).$$

$$\tilde{H}(t)=\frac{{E_0}^2}{12}\left(9\pi \sin t+\pi \sin 3t+12{\pi }^2\cos t-12\pi t\cos t\right).$$

$$E=-\sin (2t)X_1,B=\cos (2t)X_1$$

$$n_{t=0}=-\frac{1}{{\pi }^2}\underset{\{t=0\}}{\int }A\wedge F=\frac{1}{2{\pi }^2}\underset{S^3}{\int }{\theta }^1\wedge {\theta }^2\wedge {\theta }^3=1,$$

$$g=-\mathrm{d}{t}^2+a(t){\gamma }_{S^3},$$

$${\bar{\theta }}^0={\theta }^0=\mathrm{d}t,{\bar{\theta }}^i=a(t){\theta }^i,$$

$$F=E^1{\bar{\theta }}^1\wedge {\bar{\theta }}^0+B^1{\bar{\theta }}^2\wedge {\bar{\theta }}^3,$$

$$\{\begin{array}{l}\frac{d}{dt}\left(a^2{B}^1\right)=2a{E}^1,\\ \frac{d}{dt}\left(a^2{E}^1\right)=-2a{B}^1,\end{array}$$

$$E^1=\frac{E_0}{a^2}\cos \left(2\int \frac{dt}{a}\right),B^1=\frac{E_0}{a^2}\sin \left(2\int \frac{dt}{a}\right).$$

$$\frac{1}{2}\left({\left(E^1\right)}^2+{\left(B^1\right)}^2\right)=\frac{{E_0}^2}{2a^4},$$

$$G_{ab}+\Lambda g_{ab}=8\pi T_{ab}^{\text{tot}},$$

$$T_{ab}^{\text{tot}}=T_{ab}+T_{ab}^{\text{mat}}.$$

$$T_{00}^{\text{mat}}=\rho ,T_{11}^{\text{mat}}=p.$$

$$\left[\begin{array}{cccc}\mu U^0{U}^0& \mu U^0{U}^1& 0& 0\\ \mu U^0{U}^1& \mu U^1{U}^1& 0& 0\\ 0& 0& 0& 0\\ 0& 0& 0& 0\end{array}\right]+\left[\begin{array}{cccc}\mu U^0{U}^0& -\mu U^0{U}^1& 0& 0\\ -\mu U^0{U}^1& \mu U^1{U}^1& 0& 0\\ 0& 0& 0& 0\\ 0& 0& 0& 0\end{array}\right]$$

$$=\left[\begin{array}{cccc}2\mu U^0{U}^0& 0& 0& 0\\ 0& 2\mu U^1{U}^1& 0& 0\\ 0& 0& 0& 0\\ 0& 0& 0& 0\end{array}\right]=\left[\begin{array}{cccc}\rho & 0& 0& 0\\ 0& p& 0& 0\\ 0& 0& 0& 0\\ 0& 0& 0& 0\end{array}\right].$$

$$\{\begin{array}{l}3\left(\frac{1+{\dot{a}}^2}{a^2}\right)-\Lambda =4\pi \frac{{E_0}^2}{a^4}+8\pi \rho ,\\ -2\frac{\ddot{a}}{a}-\left(\frac{1+{\dot{a}}^2}{a^2}\right)+\Lambda =-4\pi \frac{{E_0}^2}{a^4}+8\pi p,\\ -2\frac{\ddot{a}}{a}-\left(\frac{1+{\dot{a}}^2}{a^2}\right)+\Lambda =4\pi \frac{{E_0}^2}{a^4}.\end{array}$$

$$\frac{a}{2}+\frac{a{\dot{a}}^2}{2}-\frac{\Lambda a^3}{6}-\frac{2\pi }{3}\frac{{E_0}^2}{a}=\frac{4\pi }{3}\rho a^3.$$

$$\rho =\frac{3M}{4\pi a^3}+\frac{{E_0}^2}{a^4},$$

$$p=\frac{{E_0}^2}{a^4}.$$

$$\frac{1}{2}{\dot{a}}^2+\left(-\frac{\Lambda a^2}{6}-\frac{2\pi {E_0}^2}{a^2}-\frac{M}{a}\right)=-\frac{1}{2}.$$

$$\frac{1}{2}{\dot{a}}^2+V(a)=h.$$

Full text

Partially funded by FCT/Portugal through projects UID/MAT/04459/2013 and PTDC/MAT-ANA/1275/2014.

Introduction

In the early days of general relativity there was a common belief that the Universe had to be eternal and unchanging. Following this concept, Albert Einstein introduced the first cosmological model, subsequently named after him. It consisted of a compact spatial hypersurface with positive curvature, the $3$-sphere $S^{3}$, which did not change under the flow of time. In order to overcome the attraction force of matter and obtain a static configuration, a new term had to be introduced in the field equations – the famous cosmological constant. Although numerous observations have since excluded this kind of static models in favor of expanding ones, the Einstein universe is still a very important explicit solution of the Einstein equations, partly because of the conformal equivalence between (a region of) this model and Minkowski’s spacetime (the standard background for quantum field theories). For example, since Maxwell’s equations are conformally invariant, one can think of any configuration of electromagnetic fields in the Einstein universe as a configuration in flat spacetime.

The Einstein universe is a particular case of, and conformal to, the closed Friedmann-Lemaître-Robertson-Walker (FLRW) models, where the assumption of staticity is relaxed by introducing a scale factor that allows the universe to expand or contract. Using this relation, one can extend solutions of Maxwell’s equations in the Einstein universe to the closed FLRW models. In general, when one allows the spatial volume to change in time, the energy of the electromagnetic configuration will also change according to an appropriate power of the scale factor.

The common denominator of the models described above is the $3$-sphere $S^{3}$. This manifold is an important example of a non-trivial fibre bundle, given by the Hopf fibration, whose base is the usual $2$-sphere $S^{2}$, and whose fibres are circles $S^{1}$ (the Clifford parallels). It is natural to look for solutions of the Maxwell equations whose electric and magnetic fields, as measured by the static observers, are aligned with these fibres. Such ansatz will be used in the first section to find a remarkable solution of Maxwell’s equations in the Einstein universe, which will then be interpreted as a knotted, finite energy, radiating electromagnetic field in Minkowski’s spacetime. The extension of this solution to closed FLRW models will be carried out in the second section. Finally, the matter distribution that must be added to obtain a self-consistent solution of the Einstein-Maxwell equations will be determined in the third section.

We follow the conventions of [1, 2]; in particular, we use the Einstein summation convention and a system of units in which $c=G=1$.

1. Electromagnetic field in the Einstein universe

We are interested in finding solutions of the Maxwell equations in the cylinder $\mathcal{E}=\mathbb{R}\times S^{3}$ with the standard Lorentzian metric. This manifold represents a static universe with positive cosmological constant $\Lambda$, called the Einstein universe. By choosing the radius of this universe as our unit length we can assume, without any loss of generality, that $\mathcal{E}$ is the product of $\mathbb{R}$ by the unit sphere.

Since $S^{3}\cong SU(2)$ is a Lie group, so is $\mathcal{E}$. We can introduce the following left-invariant orthonormal tetrad of one-forms on this manifold: Take $\theta^{0}=\mathrm{d}t$ to be the standard element of the holonomic basis on the cylinder, defining the cross section foliation. On each leaf, an orthonormal triad $\{\theta^{1},\theta^{2},\theta^{3}\}$ is chosen in such a way that111Because we have chosen an orthonormal triad, we do not have to worry about the position of the spatial indices $i,j,k,\ldots$ being subscript or superscript.

[TABLE]

thus reproducing the $\mathfrak{su}(2)$ Lie algebra structure. In terms of the dual orthonormal vector tetrad, the choice of $\theta^{0}$ and relation (1) can be written as

[TABLE]

We look for a solution of the Maxwell equations in vacuum,

[TABLE]

with the Faraday two-form

[TABLE]

Written in this basis, they are

[TABLE]

The simplest nontrivial solution can be obtained from the assumption that the components $E^{i}$ and $B^{i}$ of the electric and magnetic fields $E$ and $B$ measured by the stationary observers depend only on time. Using that ansatz, we arrive at the relations

[TABLE]

A solution of system (6) will include trigonometric functions with different initial phases and field strength factors. We choose the simplest one for the conformal analysis,

[TABLE]

because the results will not change qualitatively under more general assumptions222The property that all field lines (as seen by the stationary observers) are closed is unstable: any generic electromagnetic perturbation will destroy this feature. Moreover, such small perturbations will not decay away in time, as they would in Minkowski’s spacetime, since $S^{3}$ is compact..

Recall that Minkowski’s spacetime is conformal to an open region of the Einstein universe [2, 3]. More precisely, the Minkowski metric $\widetilde{g}$ can be written as

[TABLE]

where $(\psi,\theta,\varphi)$ are the standard hyperspherical coordinates on the 3-sphere. In these coordinates, Minkowski’s spacetime corresponds to the region $0\leq|t|+\psi<\pi$ (Fig. 1).

Since Maxwell’s equations are conformally invariant, we can interpret (7) as an electromagnetic field in flat spacetime333See [4, 5] for a similar contruction exploiting the conformal relation between Minkowski’s and de Sitter’s spacetimes.. The energy will of course change under the conformal transformation. Its value will be given by the integral expression

[TABLE]

where

[TABLE]

is the energy-momentum tensor of $\widetilde{F}$, $\widetilde{X}_{0}=\Omega X_{0}$ is $X_{0}$ normalized by the Minkowski metric $\widetilde{g}$, and $\widetilde{\Sigma}_{t}$ is the spacelike hypersurface corresponding to the hyperspherical cap $\Sigma_{t}$ of constant $t$ (Fig. 1). Note that $\widetilde{\Sigma}_{t}$ is not a Cauchy surface for Minkowski’s spacetime unless $t=0$, as it approaches $\mathscr{I^{-}}$ for $t<0$ and $\mathscr{I^{+}}$ for $t>0$.

Using $\widetilde{T}_{ab}=\Omega^{2}{T}_{ab}$ and ${\mathrm{d}(\widetilde{vol})}=\Omega^{-3}\mathrm{d}(vol)$, we can reexpress the energy as an integral over the hyperspherical cap $\Sigma_{t}$,

[TABLE]

It is straightforward to see that this quantity is finite. Moreover, it is a decreasing function of the parameter $t>0$, which indicates that the solution under consideration describes radiation fields. This can be seen by looking at the time dependence of the energy (Fig. 2.), given, for $t>0$, by the formula

[TABLE]

The electromagnetic field (7) is related to the simplest initial conditions given by Rañada in [6, Section 2] (see also [7, 8]), where only the initial magnetic field is nonzero. In fact, the phase-shifted solution (with $E_{0}=1$)

[TABLE]

is the Cauchy development of initial conditions of that type. Due to the fact that it is aligned along the Clifford parallels, the Hopf index $n$ is equal to one in our case. This fact can be confirmed using an integral expression for this quantity, namely

[TABLE]

where $F=\mathrm{d}A$.

2. Expanding models

In this section we turn our attention to the expanding models. Consider a metric of the form

[TABLE]

where $\gamma_{S^{3}}$ is the standard metric of $S^{3}$ with unit radius; this corresponds to a closed FLRW model, widely used in modern cosmology [1].

In analogy with the previous case, we introduce the tetrad $\{\bar{\theta}^{0},\bar{\theta}^{1},\bar{\theta}^{2},\bar{\theta}^{3}\}$,

[TABLE]

and the dual vector basis $\{\bar{X}_{0},\bar{X}_{1},\bar{X}_{2},\bar{X}_{3}\}$.

We are looking for a solution of the Maxwell equations similar to (7), namely with a Faraday tensor of the form

[TABLE]

where $E^{1}$ and $B^{1}$ depend only on time. Substituting into the vacuum Maxwell equations, (3), we have

[TABLE]

yielding

[TABLE]

The Einstein universe case corresponds to $a\equiv 1$. Note that the $T_{00}$ component of the electromagnetic stress-energy tensor is

[TABLE]

indicating that the energy density of the solution is decaying with the fourth power of the radius of the universe.

3. Matter content

The solutions of Maxwell’s equations discussed in the last sections were derived under the assumption of a fixed background metric, satisfying the Einstein field equations

[TABLE]

where $G_{ab}$ is the Einstein tensor, $g_{ab}$ is the metric and $T^{\text{tot}}_{ab}$ is the total stress-energy tensor. One can ask what kind of matter has to be present for this configuration to be self-consistent. In order to answer this question, we consider the decomposition of the full stress-energy tensor into the parts coming from the electromagnetic field and matter,

[TABLE]

We already used the expression for the time-time component of $T_{ab}$ when computing the energy of the solution (7); the full expression for this tensor is given by (10). It is straightforward to see that it has diagonal form for any solution of the system (18), $T_{ab}=\operatorname{diag}(D,-D,D,D)$, where in general $D=D(t)$.

Let us start from the assumption that matter can be described by a density $\rho$ and a principal pressure $p$ aligned with $X_{1}$,

[TABLE]

A possible interpretation of this anisotropic pressure can be given by superimposing the stress-energy tensors of two pressureless fluids with the same mass density $\mu$, moving in opposite directions along $X_{1}$, with the same velocity with respect to the fundamental observers:

[TABLE]

Another possible interpretation can be given by considering cosmic strings aligned with the Clifford parallels, expanding and contracting with space while preserving the structure of the Hopf fibration.

Because the metric in this model is similar to that of the closed FLRW universe, we end up with a system analogous to the Friedmann equations with the additional terms from electromagnetic stress-energy tensor,

[TABLE]

As it turns out, we can deal with system (23) using a similar procedure as in the standard case [1]. The first equation is equivalent to

[TABLE]

After taking a time derivative of both sides and using the third equation from (23), we find

[TABLE]

where $M$ is another integration constant. We can use now the second equation from (23) to obtain the pressure $p$ as

[TABLE]

Note that the matter satisfies the weak, strong and dominant energy conditions if $M\geq 0$, as in this case $0<p\leq\rho$. Nevertheless, $p/\rho$ always approaches $1$ as $a$ tends to zero, indicating that the matter particles approach the speed of light near the Big Bang or the Big Crunch. In the limit case $M=0$ we have $p=\rho$, meaning that the particles are actually moving at the speed of light.

The behavior of the scale factor $a$ can be derived from an effective potential formulation. Using (24) and (25) we have

[TABLE]

We can think about each term as a kinetic part, an effective potential $V(a)$, and a conserved energy $h=-1/2$, respectively,

[TABLE]

The potential is a function of $a$, but also depends on the three parameters $\Lambda$, $M$ and $E_{0}$. Clearly we need a positive cosmological constant to have an (unstable) equilibrium point, so that we can reproduce the previous results concerning the Einstein universe. In this case, $V(a)$ will have the shape depicted on Fig. 3. Depending on the values of the parameters, we can have expanding solutions, recollapsing solutions, de Sitter-like bouncing solutions, and of course the unstable static solution and its asymptotes. Note that in units for which the radius of the static solution is $a=1$, we have

[TABLE]

(hence we must have ${E_{0}}^{2}\leq\frac{1}{8\pi}$).

Finally, we note that except in the case of the Einstein universe and its two asymptotes, the electromagnetic field performs a finite number $n$ of oscillations during the history of the universe, proportional to its lifetime as measured in conformal time:

[TABLE]

Therefore, at any given epoch these fields will appear mostly constant and may perhaps serve as models of primordial magnetic fields [9].