1 Introduction

The aim of this work is to show that there exists a correspondence between the seemingly different physical models of a cosmological model using dark energy and Maxwell equations. The link between these two formalisms appears when one considers the quantum version of the cosmological model [1] using the Breit prescription [2] for spin particles.

The representation of dark energy using a quintessence field with a potential, allows us to describe the cosmological dynamics in fashion which is analogous to the description of the dynamics of a relativistic particle. In the case of quintessence described by a free massless scalar field in a Universe with cosmological constant, the resultant theory may be quantized by using a Klein–Gordon scheme, giving rise to the Wheeler–DeWitt equation [3,4,5,6,7,8,9,10,11,12,13]. However, using the Breit prescription [2], the same model can be quantized as a spinorial theory. This procedure yields a Majorana version for the quantum cosmology which happens to be supersymmetric [1].

In addition of the cosmological implications of such theory, the aim of this article is to show its direct correspondence with the description of propagating electromagnetic fields in a medium using Maxwell equations. We can identify the relative permittivity and permeability of the medium with parameters of the quantum cosmological model. As we show, this implies that the supersymmetric Majorana quantum cosmology can be studied in an analogue electromagnetic system using either normal materials or negative-index metamaterials (NIMs) [14, 15].

2 Quantum cosmology with dark energy

Consider an isotropic and homogeneous Friedmann–Robertson–Walker (FRW) spacetime with a line element [16]

$$\begin{aligned} ds^2=dt^2- a(t)^2\left[ \frac{dr^2}{1-kr^2} +r^2(d\theta ^2+ \sin ^2\theta d\phi ^2)\right] , \end{aligned}$$
(1)

where a(t) is the scale factor, and the curvature constant \(k=\pm 1,0\). The evolution of a FRW cosmology with cosmological constant \(\Lambda \), interacting with a quintessence (massless scalar) field \(\phi (x^\beta )\) characterized by a potential \({{\mathcal {V}}}(\phi )\), can be found by using Einstein equations [1] (\(8\pi G/c^4=1\), where G is the gravitational constant and c is the speed of light)

$$\begin{aligned} 2\frac{\ddot{a}}{a}+\left( \frac{\dot{a}}{a}\right) ^2 +\ \frac{k}{a^2} +\frac{1}{2}\ {{\dot{\phi }}}^2-V(\phi )= & {} 0\, , \end{aligned}$$
(2)
$$\begin{aligned} 3\left( \frac{\dot{a}}{a}\right) ^2 +3\frac{k}{a^2}-\left( \frac{1}{2}{{\dot{\phi }}}^2+ V(\phi ) \right)= & {} 0 , \end{aligned}$$
(3)

where we introduced \(V(\phi )={{\mathcal {V}}}(\phi ) - \Lambda \). Also, the Klein–Gordon equation for the quintessence field is

$$\begin{aligned} \ddot{\phi }+ 3\frac{\dot{a}}{a} {{\dot{\phi }}} + \frac{d V(\phi )}{d\phi }=0 . \end{aligned}$$
(4)

The above system describes the evolution of a Friedmann–Robertson–Walker–Quintessence (FRWQ) Universe.

In Ref. [1] was shown that the Lagrangian

$$\begin{aligned} L_F= \sqrt{{\bar{V}}(\xi , \theta )\, e^{2\xi } ( {{\dot{\theta }}}^2- {{\dot{\xi }}}^2 )} , \end{aligned}$$
(5)

gives rise to all three equations (2)–(4), with variables \(\xi =\ln ({2\sqrt{6}} a^{3/2}/3)\), \(\theta ={3\phi }/({2\sqrt{6}})\) [here \({{\dot{\xi }}}={d\xi }/{d\lambda }\) and \({{\dot{\theta }}}={d\theta }/{d\lambda }\), with \(\lambda \) is an arbitrary parameter], and the general potential

$$\begin{aligned} {\bar{V}}(\xi ,\theta )=3 \left( \frac{3}{8}\right) ^{{1}/{3}}k\, e^{2\xi /3}-\left( \frac{3}{8}\right) e^{2\xi } V(\theta ) . \end{aligned}$$
(6)

This Lagrangian (5) shows that the FRWQ cosmology evolves as a relativistic particle moving in two dimensional spacetime under the influence of potential \({\bar{V}}\). It can be proved that the Lagrangian (5) gives rise to the FRWQ equations by recalling that the Jacobi–Maupertuis and Fermat principles [17] yield identical equations of motion in classical mechanics and geometrical (ray) optics except for the fact that Fermat principle also produces a constraint equation. Notice that Lagrangian (5) can be written as the one for a relativistic particle in a two-dimensional conformally flat spacetime

$$\begin{aligned} L_F= \sqrt{g_{\mu \nu } \left( \frac{dx^\mu }{d\lambda }\right) \left( \frac{dx^\nu }{d\lambda }\right) } , \end{aligned}$$
(7)

with the metric \(g_{\mu \nu }=\Omega ^2\, \eta _{\mu \nu }\) (where \(\eta _{\mu \nu }\) is the flat spacetime metric, and \(x^0=\theta \) and \(x^1=\xi \)), and the conformal factor

$$\begin{aligned} \Omega \equiv \sqrt{{\bar{V}}}e^\xi =\left[ 3 \left( \frac{3}{8}\right) ^{\frac{1}{3}}k\, e^{8\xi /3} -\frac{3}{8} V(\theta ) e^{4\xi }\right] ^{1/2} , \end{aligned}$$
(8)

Thus, the FRWQ system is equivalent to a relativistic particle moving in a two-dimensional conformally flat spacetime, where the quintessence field plays the role of an effective time.

The above results allow to find the quantum version of the FRWQ cosmology, through quantization of the Lagrangian (5). We restrict ourselves to cases where \({\bar{V}} >0\), in order to avoid problems in the quantization procedure. We also consider static manifolds [18] such that there exists a family of spacelike surfaces orthogonal to a timelike Killing vector. Therefore, Lagrangian (5) is always real, which can be proved by using the definitions of \(\xi \), \(\theta \) and \({\bar{V}}\). All of this implies that \({\partial _\theta g_{\mu \nu }}=0\), or \({\partial _\theta }{\bar{V}}=0\), describing a constant potential \(V(\theta )\) (becoming essentially equal to a cosmological constant). Thus, the associated Noether conservation law implies that \({{\dot{\theta }}}\) does not change sign, and \(\theta \) can be used as the evolution time variable. Using this \(\theta \)–time, it can be shown at classical level [19] that the Hamiltonian for Lagrangian (5) is \( H=\sqrt{g_{00}}\sqrt{1-g^{11}\pi ^2}=\sqrt{g_{00}}\sqrt{1+{\pi ^2}/{\Omega ^{2}}}\), where \(\pi \) is the canonical momentum, and \(\sqrt{g_{00}}=\Omega \). This classical Hamiltonian is used to construct the quantum Hamiltonian operator \({\mathcal {H}}\) as [1]

$$\begin{aligned} {\mathcal {H}}=\Omega ^{1/2}\sqrt{1+{{\hat{p}}}^2}\, \Omega ^{1/2} , \quad {{\hat{p}}}=\sqrt{-g^{11}}{{\hat{\pi }}}=-\frac{i}{\Omega }\frac{\partial }{\partial \xi }\, , \end{aligned}$$
(9)

where \({{\hat{p}}}\) is the momentum operator and \({{\hat{\pi }}}=-i{\partial }_\xi \). Therefore, the quantum equation that describes the quantization of the FRWQ system is

$$\begin{aligned} i\hbar \frac{\partial \Psi }{\partial \theta }={\mathcal {H}}\Psi , \end{aligned}$$
(10)

where \(\Psi \) is the wavefunction for the quantum FRWQ cosmology. The quantum equation (10) emerges because of the direct correspondence between the FRW geometry and the quintessence scalar field at a Fermat-like Lagrangian level. To find a solution, we need a quantization procedure to solve the square-root of the Hamiltonian operator (9). It is costumary to solve the square-root using a spinless particle approach to get a Klein–Gordon equation [20], which gives origin to the Wheeler–DeWitt super-Hamiltonian formalism. However, notice that there is no restriction for the quantization scheme possible to be used. In principle, we can solve the square-root using matrices, obtaining the quantization of a relativistic particle which leads to the Dirac equation. This procedure was introduced by Breit [2], that shows that there is a correspondence between the Dirac and the relativistic pointlike particle Hamiltonians. Breit’s interpretation [2] identifies the Dirac matrices as \({\pi }/{H}\rightarrow \alpha \), and \(\sqrt{1-({\pi }/{H})^2}\rightarrow \beta \). These identifications are consistent with the postulates of Dirac electron’s theory [2, 21, 22]. By the Breit’s prescription, the Hamiltonian operator (9) becomes [1]

$$\begin{aligned} {{\mathcal {H}}}=\Omega ^{1/2} (\alpha \cdot {{\hat{p}}}+\beta ) \Omega ^{1/2} , \end{aligned}$$
(11)

where \(\alpha \) and \(\beta \) are the two-dimensional flat spacetime Dirac matrices, as the effective curvature is already taken into account in \(\Omega \) (with \(\hbar =1\)). Moreover, now the wavefunction \(\Psi \) [in Eq. (10)] is a two-dimensional spinor.

With this Hamiltonian, and defining the wavefunction \(\Phi =\sqrt{\Omega }\Psi \), we finally find from Eq. (10) the spinor quantum equation [1, 23]

$$\begin{aligned} i\gamma ^0\frac{\partial \Phi }{\partial \theta }+ i\gamma ^1\frac{\partial \Phi }{\partial \xi }=\Omega \Phi . \end{aligned}$$
(12)

where \(\gamma ^0=\beta \) and \(\gamma ^1=\gamma ^0\alpha \). The above equation corresponds to a Quantum Cosmology theory for the FRWQ system, modelling now the Universe as a spin particle in a two-dimensional conformally flat spacetime [1].

In order to obtain real wavefunctions \(\Phi \), the matrices in Eq. (12) should correspond to the two-dimensional Majorana representation

$$\begin{aligned} \gamma ^0=\left( \begin{array}{c@{\quad }c} 0 &{} -i\\ i &{}0\\ \end{array} \right) , \quad \gamma ^1=\left( \begin{array}{c@{\quad }c} i &{} 0\\ 0 &{}-i\\ \end{array} \right) . \end{aligned}$$
(13)

In this form, Eq. (12) becomes a set of supersymmetric equations of quantum mechanics [1, 24,25,26,27,28], which can only be obtained in the Majorana picture. The implications of this system were thoroughly studied in Ref. [1], showing that Eq. (12) with matrices (13) represents a Supersymmetric Majorana quantum cosmology. This can be seen by defining

$$\begin{aligned} \Phi (\theta ,\xi )=\left( \begin{array}{c} \varphi _+(\xi ) e^{E\theta }\\ \varphi _-(\xi ) e^{-E\theta } \\ \end{array} \right) , \end{aligned}$$
(14)

to find from (12) the set of supersymmetric equations of quantum mechanics [1, 23,24,25,26,27,28]

$$\begin{aligned} Q_\pm \varphi _\pm =E\varphi _\mp . \end{aligned}$$
(15)

with the operators \(Q_\pm =\pm {d_\xi }+\Omega \). These equations are supersymmetric with the two spinor components being super–partners of each other. Each wavefunction satisfies \(H_\pm \varphi _\pm ={E^2}\varphi _\pm \), with the Hamiltonians operators \(H_\pm =-{d_\xi ^2}+{{\mathcal {W}}}_\pm \), and potentials \({{\mathcal {W}}}_\pm =\mp {d_\xi \Omega }+\Omega ^2\) [1]. Also, \(Q_\pm \) correspond to supercharge operators. This theory can produce new versions of quantum cosmological models. For example, it can be shown that in flat curvature case, the Universe behaves as a diatomic molecule subject to the Morse potential [1].

3 Correspondence to Maxwell equations

In this section we focus the correspondence between the theory (12) and electromagnetism. In order to make this correspondence manifest, let us consider Maxwell equations in a medium with permittivity \(\epsilon \) and permeability \(\mu \) in the absence of charges [29]

$$\begin{aligned} \frac{\partial \mathbf{D}}{\partial t}= & {} \nabla \times \mathbf{H} ,\nonumber \\ \frac{\partial \mathbf{B}}{\partial t}= & {} -\nabla \times \mathbf{E} , \end{aligned}$$
(16)

with the electric field \(\mathbf{E}\), the magnetic field \(\mathbf{B}\), the displacement field \(\mathbf{D}=\epsilon \, \mathbf{E}\), the magnetization field \(\mathbf{H}=\mathbf{B}/\mu \) (chosing the speed of light \(c=1\)). Furthermore, \(\nabla \cdot \mathbf{D}=0\) and \(\nabla \cdot \mathbf{B}=0\).

In order to show the correspondence with the Quantum Cosmology model, we assume that the permittivity and the permeability are not constant. Also, let us consider a two-dimensional spacetime system (one temporal and one spatial dimension), with spatial variations in, let us say, the \({{\hat{e}}}_z\)-direction. We choose transverse fields \(\mathbf{B}(t,z)=B(t,z) {{\hat{e}}}_x\), and \(\mathbf{D}(t,z)=D(t,z) {{\hat{e}}}_y\), such that \(\mathbf{B}\cdot {{\hat{e}}}_z=0=\mathbf{D}\cdot {{\hat{e}}}_z\), and \(\mathbf{B}\cdot \mathbf{D} =0\). Besides, the time-independent permittivity and permeability have the same spatial dependence, \(\epsilon =\epsilon (z)\) and \(\mu =\mu (z)\). Then, Maxwell equations (16) acquire the form

$$\begin{aligned} \frac{\partial {{\mathcal {D}}}}{\partial t}= & {} -\frac{B}{{\hat{\mu }}^{2}}\frac{\partial {{\hat{\mu }}}}{\partial z}+\frac{1}{{{\hat{\mu }}}}\frac{\partial B}{\partial z}\, ,\nonumber \\ \frac{\partial B}{\partial t}= & {} -\frac{{{\mathcal {D}}}}{{\hat{\epsilon }}^{2}}\frac{\partial {{\hat{\epsilon }}}}{\partial z}+\frac{1}{{{\hat{\epsilon }}}}\frac{\partial {{\mathcal {D}}}}{\partial z} . \end{aligned}$$
(17)

where \({{\mathcal {D}}}=\mu _0 D\), the relative permittivity is \({{\hat{\epsilon }}}=\epsilon /\epsilon _0\), the relative permeability is \({{\hat{\mu }}}=\mu /\mu _0\), and \(\epsilon _0\) and \(\mu _0\) are the free-space permittivity and permeability respectively (with \(\epsilon _0\mu _0=1\)).

By introducing the spinor

$$\begin{aligned} \Phi =\left( \begin{array}{c} B \\ {{\mathcal {D}}} \end{array} \right) , \end{aligned}$$
(18)

we can rewrite system (17) in a simple way as

$$\begin{aligned} i\gamma ^0\frac{\partial \Phi }{\partial t}+i\Gamma _1\frac{\partial \Phi }{\partial z}=\Gamma _2\Phi , \end{aligned}$$
(19)

with \(\gamma ^0\) given in (13), and the matrices

$$\begin{aligned} \Gamma _1=\left( \begin{array}{c@{\quad }c} i/{{\hat{\mu }}} &{} 0 \\ 0 &{} -i/{{\hat{\epsilon }}} \end{array} \right) ,\quad \Gamma _2=\left( \begin{array}{c@{\quad }c} -{{\hat{\mu }}}'/{{\hat{\mu }}^{2}} &{} 0 \\ 0 &{} {{\hat{\epsilon }}}'/{\hat{\epsilon }}^{2} \end{array} \right) ,\quad \end{aligned}$$
(20)

with \({{\hat{\mu }}}'={\partial {{\hat{\mu }}}}/{\partial z}\) and \({{\hat{\epsilon }}}'={\partial {{\hat{\epsilon }}}}/{\partial z}\). The spinor form (19) of Maxwell equations is general for a two-dimensional spacetime. For any general permittivity and permeability, Eq. (19) does not coincide with the Quantum cosmology equation (12). For example, for vacuum, \(\Gamma _2=0\).

However, we can show that there exists a regime in which Maxwell equations and the Supersymmetric Majorana Quantum Cosmology coincide. Let us consider the following form for the relative permittivity and relative permeability

$$\begin{aligned} {{\hat{\epsilon }}}(z)\approx \pm 1+\lambda (z) ,\quad {{\hat{\mu }}}(z)\approx \pm 1-\lambda (z) , \end{aligned}$$
(21)

in terms of a function \(\lambda \) to be determined. We focus our attention in the regime when \(|\lambda |\ll 1\). The relative permeability and permittivity of this material are almost equal, \(\sqrt{{{\hat{\mu }}}/{{\hat{\epsilon }}}} \approx 1-\lambda \). With the relative permittivity and permeability given by (21), the matrices (20) becomes

$$\begin{aligned} \Gamma _1\approx & {} \pm ({1+\lambda ^2} ) \gamma ^1\pm {i\lambda }\, \text{1 } ,\nonumber \\ \Gamma _2\approx & {} \frac{d \lambda }{d z}\, \text{1 } \mp i\frac{d\lambda ^2}{d z}\gamma ^1 , \end{aligned}$$
(22)

with the matrix \(\gamma ^1\) given by (13), and the two-dimensional unit matrix \(\mathbf{1}\). In general, for any \({{\hat{\epsilon }}}\) and \({{\hat{\mu }}}\), Eq. (19) does not have a conserved current \({\bar{\Phi }}\gamma ^\mu \Phi =\Phi ^\dag \Phi +\Phi ^\dag \gamma ^0\gamma ^1\Phi \). Nonetheless, for materials satisfying (21) and (22), Eq. (19) can have a conserved current if \(|\lambda |\ll 1\) and that the variation ranges of B and D are much larger than the variation range of \(\lambda \), i.e., \(\lambda (\partial \Phi /\partial z)\ll (d\lambda /dz)\Phi \). Using these approximations, Maxwell equation (19) can be written in an approximated form as

$$\begin{aligned} i\gamma ^0\frac{\partial \Phi }{\partial t}+i\gamma ^1\frac{\partial \Phi }{\partial z}=\frac{d \lambda }{d z}\Phi , \end{aligned}$$
(23)

which has the usual conservartion law for a Dirac equation. With all of the above, the correspondence between Maxwell equations (23) and the supersymmetric quantum Majorana cosmologies equation (12) is now evident by redefining

$$\begin{aligned} t= & {} \theta ,\nonumber \\ \xi= & {} \pm z ,\nonumber \\ \lambda= & {} \lambda (\xi )= \int \Omega \, d\xi =-\frac{4}{3V}e^{-4\xi }\, \Omega ^3 . \end{aligned}$$
(24)

This last equation implies that this correspondence is only valid in the regime \({\Omega ^3}e^{-4\xi }\ll {|V|}\), which is the limit needed in order to keep current conservation.

Results (24) establish the complete correspondence between a Supersymmetric Majorana Quantum cosmological model and Maxwel equations. Materials satisfying (21) must have almost equal relative permeability and permittivity, and they can be both positive or both negative. In the former case, we are describing a normal material with such a property. In case that both permeability and permittivity are negative, we are describing a NIM [14, 15, 30]. In general, the refraction index is

$$\begin{aligned} n(\xi )=\sqrt{{{\hat{\epsilon }}}(\xi ){{\hat{\mu }}}(\xi )}\approx \pm \left( 1-\frac{8}{9V^2}e^{-8\xi } \Omega ^6\right) , \end{aligned}$$
(25)

where the positive sign describes a normal material, and the negative refractive index represents a NIM in which a wave propagates backwards.

There are materials with relative permittivity and relative permeability that can have the form (21). Composite ferrites [31,32,33,34] can, under appropriated conditions, achieve almost matching permeability and permittivity values by shining radiation of different frequencies on the material. This implies that the above results can be tested in an analogue fashion using those materials. For example, for a spatially flat cosmology, the Maxwell equations and the Supersymmetric Majorana quantum cosmology coincide for \(\lambda (\xi )=\sqrt{-3 {V}/32} \exp (2\xi )\), with constant \(V<0\), and refraction index \(n(\xi )\approx \pm [1+3 V\exp (4\xi )/64]\). For this case, Maxwell equations can only describe a quantum cosmology for \(\xi \ll \ln (11/|V|)\). Furthermore, the boundary conditions for \(\Phi \) given in Eq. (18) (which depend on \(\Omega \)) should be suitable for simulating a quantum cosmology. For \(k=0\), the boundary conditions for \(\xi \rightarrow \infty \) (\(a\rightarrow \infty \)) can be established to obtain a vanishing magnetic field at infinity in one spatial dimension (see Ref. [1]).

Our proposal is in the same spirit than similar ones for analogue optical systems for quantum cosmologies [35, 36], and for gravity in general (see for example Refs. [37,38,39,40]). However, our result establishes an “optical” analogue for a new kind of spinor quantized cosmological model. The proposed analogue electromagnetic media that correspond to the quantum cosmology is time-independent but space-dependent, which is an approach opposite to previous attempts [35]. Relations (21) are satisfied by certain tunable metamaterials [41, 42] and composite ferrites [31,32,33,34] used to operate at a wider range of frequencies. All of the above makes of this Majorana Supersymetric quantum cosmological model a system worth to be studied by studying wave propagation in Maxwell equations in the appropriated media.