Home Electrodynamics in Euclidean Space Time Geometries
Article Open Access

Electrodynamics in Euclidean Space Time Geometries

  • Jörn Schliewe EMAIL logo
Published/Copyright: December 30, 2019

Abstract

In this article it is proven that Maxwell’s field equations are invariant for a real orthogonal Cartesian space time coordinate transformation if polarization and magnetization are assumed to be possible in empty space. Furthermore, it is shown that this approach allows wave propagation with finite field energy transport. To consider the presence of polarization and magnetization an alternative Poynting vector has been defined for which the divergence gives the correct change in field energy density.

1 Introduction

The foundation of the electromagnetic theory is the work of many people and culminated in the publications of J. C. Maxwell in 1865 [1] and 1873 [2]. The field equations in the modern form, known as Maxwell’s equations, are validated in countless applications, like antennas, radar and electronic components. Around 1900 H. A. Lorentz [3] and others [4] tried to explain the behavior of electromagnetic fields in moving mater. Background was the idea of an ether which transports the field energy of light. In 1905 A. Einstein published his famous theory of special relativity [5]. This was the start of a new way of thinking. Electromagnetic energy moves without any carrier charges through empty space. Einstein’s theory which is based on the geometrical properties of the Lorentz transformation is enforced by Maxwell’s electromagnetic theory of empty space, e.g. [6, p. 660], [7, p. 645] and [8, p. 595]. The criticism on Einstein’s initial postulate of constant velocity of light has been refuted by the fact that Maxwell’s equations of empty space are invariant under Lorentz transformation. W. Pauli said: "It is clear that one would then have to abandon not only the idea of the existence of an aether but also Maxwell’s equations for the vacuum, so that the whole of electrodynamics would have to be constructed anew" [4, p. 6]. Attempts to find new electromagnetic theories like the one of Ritz [9] failed to explain all experiments. Nevertheless, the theory of Ritz has been considered to be noteworthy as Einstein’s initial postulate of constant velocity of light, even though supported by experiments (e.g. [10, 11]), is quite contrary to our experience with mechanical or acoustical waves [12, p. 75]. However, from a theoretical point of view, it would be preferable if this constant velocity would be a result and not an initial postulate of the theory.

Starting only from the principle of relativity and basic geometrical considerations [13, 14] or [15] showed that a set of transformation groups satisfies the principle of relativity without the initial assumption of a constant velocity of light. Different transformations can be distinguished by a universal space time constant. These transformations include the Galilean transformation as a trivial solution where the velocity of light tends towards infinity [16], the Lorentz transformation, but also a transformation based on a simple orthogonal Cartesian rotation in space time. As the Newtonian mechanics are invariant under a Galilean transformation rigorous discussions to validate the Lorentz transformation based theory of relativity against the Galilean transformation can be found in literature, e.g. [11, 16]. However, a rigorous comparison of the Lorentz transformation with a simple orthogonal Cartesian rotation in space time cannot be found.

Beginning around 1949 the emptiness of the space was somehow relativized by the possibility of vacuum fluctuations which allow the creation of particle-antiparticle pairs of virtual particles in empty space [17]. Quantum electrodynamics which successfully predicted the Lande factor of magnetic moment of electrons relies on the concept of virtual particles [18, 19].

The main aim of this article is to give a theoretical justification for the necessity of a comparison of the Lorentz transformation based theory of relativity with the simple orthogonal Cartesian rotation in space time. The two main theoretical statements presented in this article support this necessity. Firstly, it will be proven that Maxwell’s field equations are invariant under a simple orthogonal Cartesian rotation in space time if electric and magnetic polarization are considered to happen in space. Secondly, it will be shown that this approach leads to bounded wave solutions which contains finite field energy transport in contrast to the plane wave solution.

2 Definitions

In this section the most important definitions are discussed on the basis of Maxwell’s field equations [6, p. 34], [7, p. 31], [20, p. 425] or [8, p. 3]. Ampers law 1 and the law of induction 2 can be written in the following form:

(1) rot(H)=J+ϵ(P+E)t
(2) rot(E)=μ(M+H)t

rot () is the curl, H is magnetic field strength, J is current density, P is polarization, E is electric field strength, M is magnetization, ε is the permittivity of empty space which gives the ability for polarization and finally μ is the permeability of empty space which gives the ability for magnetization.

Using these definitions the electric flux density results in

(3) D=ϵ(P+E)

and the magnetic flux density results in

(4) B=μ(M+H)

These equations imply the basic definitions of the polarization and magnetization which are used throughout this paper. The permittivity ε and permeability μ describe the basic properties of empty space. Using both properties the velocity of light c and the impedance of empty space Z0 can be derived:

(5) c=1ϵμ
(6) Z0=μϵ.

Both quantities are properties of propagating electromagnetic field waves and show directely that the properties of empty space define the properties of wave propagation.

3 Orthogonal Cartesian Space Time Transformation

As a most basic approach this transformation can be derived from general group-theoretical considerations like in [13], [14] or [15]. Considering two systems moving relative to each other in z direction this leads to [4, p. 11]:

(7) x=x
(8) y=y
(9) z=11αv2c2(zvt)
(10) t=11αv2c2(tαvc2z)

For α = 0 this represents the Galilean transformation, for α = 1 this gives the Lorentz transformation and finaly with α = −1 this results in an orthogonal Cartesian space time transformation:

(11) x=x
(12) y=y
(13) z=11+v2c2(zvt)
(14) t=11+v2c2(t+vc2z)

In the following γ will be defined as γ=11+v2c2.

Figure 1 shows the geometrical representation of an orthogonal Cartesian space time transformation relation between two systems moving relative to each other in z direction. It can be seen that the orthogonal relation between time and z is ensured which represents a simple rotation in an Euclidean space. In reflection of this simple rotation neither the Lorentz transformation nor the Galilean transformation is orthognal Cartesian in four dimensional space.

Figure 1 Geometrical representation of an orthogonal space time transformation relation between two relatively in z direction moving systems
Figure 1

Geometrical representation of an orthogonal space time transformation relation between two relatively in z direction moving systems

For the proof of invariance the transformation of the differentials is essential. As space and time depend on each other the differentials have to be derived by use of the chain rule:

(15) x=x
(16) y=y
(17) z=zzz+tzt=γ(z=vc2t)
(18) t=ttt+ztz=γ(tvz)

4 Proof of Invariance

4.1 Current density, polarization and magnetization

To prove Maxwell’s equations to be invariant under the Euclidean space time transformation the role of polarization and magnetization has to be changed from an amplification of the field quantities by εr and μr to quantities which influence even the empty space. Polarization is the separation of virtual particles and magnetization is the rotation or spinning of virtual particles. With this definition the polarization and magnetization can exhibit a different distribution in space compared to the field quantities.

As explained in [6, p. 29] or [7, p. 27] the introduction of the time derivative of the electric flux density in ampers law 1 is necessary to fulfill the equation of continuity. Figure 2 illustrates the relation with a simple example. A wire connects two electrodes. If a current reduces the charges on the electrodes, the corresponding electric field should be reduced in the same way.

Figure 2 Example for the continuity equation
Figure 2

Example for the continuity equation

The divergence of the right hand side of Ampers law 1 should be equal to zero. This gives:

(19) 0=div(J+ϵPt+ϵEt)

As dicussed in [6, p. 88] or [7, p. 87] the time derivative of the polarisation can be understood as densities of polarisation currents and consequentely the so called bound charge densities can be introduced in contrast to the free charge densities ρ. As "bound" sounds too restricting they will be decribed as virtual charge densities ρv throughout this article.

In the following the electric field will be defined by div (ϵE)=ρ+ρv and the electric flux density by div (D) = ρ. With the latter the equation of continuity is defined.

(20) 0=div(J)+ρt

Proving the equation of continuity for invariance gives:

0=div(J)+ρt0=Jxx+Jyy+Jzz+ρt0=Jxx+Jyy+γ(z+vc2t)Jz+γ(tvz)ρ0=Jxx+Jyy+zγ(Jzvρ)+tγ(ρ+vc2Jz)

According to this derivation the following transformation rules should apply:

(21) Jx=Jx
(22) Jy=Jy
(23) Jz=11+v2c2(Jzvρ)
(24) ρ=11+v2c2(ρ+vc2Jz)

A similar discussion is necessary for the magnetic flux density. The divergence of the right hand side of the law of induction 2 should be equal to zero. This gives:

(25) 0=div(B)=μdiv(M)+μdiv(H)

In contrast to the usual approach div (M) and div is (H) not necessarily equal to zero. In fact this has been proven to exist in spin ice materials [21, 22, 23]. The physical manifestation of these divergences are not magnetic monopoles but the spin of virtual particles. That the magnetic field and the magnetization are not source free, is nothing unexpected for calculations of magnetostatic fields on boundaries of magnetizable materials [6, p. 303] or [7, p. 297].

4.2 Invariance of Maxwell’s equations

The same method used for the equation of continuity is now used to prove the invariance of Amperes law.

(26) rot(H)=J+ϵPt+ϵEt
(27) 0=rot(H)JϵPtϵEt

This can be separated to three equations:

(28) 0=HzyHyzJxϵPxtϵExt
(29) 0=HxzHzxJyϵPytϵEyt
(30) 0=HyxHxyJzϵPztϵEzt

Firstly, transforming equation 28:

(31) 0=Hzyγ(z+vc2t)HyJxγ(tvz)ϵPxγ(tvz)ϵEx=Hzyzγ(HyvϵPxvϵEx)JxγϵPxtϵtγ(Ex+vμHy)

To get the same field transformation rule for the electric field like Einstein, the equation has to be expanded to:

(32) 0=Hzyzγ(HyvϵPxϵEx)JxγϵPxtγvc2t(My+2Hy)ϵtγ(Exvμ(My+Hy))=Hzyzγ(Hyvϵ(Px+Ex))Jxϵtγ(Px+vμ(My+2Hy))ϵtγ(Exvμ(My+Hy))

Then equation 29:

(33) 0=γ(z+vc2t)HxHzxJyγ(tvz)ϵPyγ(tvz)ϵEy=zγ(Hx+vϵPy+vϵEy)HzxJyγϵPytϵtγ(EyvμHx)zγ(Hx+vϵPy+vϵEy)HzxJyγϵPyt+γvϵμt(Mx+2Hx)ϵtγ(Ey+vμ(Mx+Hx))=zγ(Hx+vϵ(Py+Ey))HzxJyϵtγ(Pyvμ(Mx+2Hx))ϵtγ(Ey+vμ(Mx+Hx))

To ensure invariance according to the real orthogonal spacetime for both equations, the following field transformation equations have to be used:

(34) Ex=γ(Exvμ(Mx+Hy))Ey=γ(Ey+vμ(Mx+Hx))Hx=γ(Hx+vϵ(Py+Ey))Hy=γ(Hyvϵ(Px+Ex))Hz=HzPx=γ(Px+vμ(My+2Hy))Py=γ(Pyvμ(Mx+2Hx))

Last but not least the relations in 34 can be applied on equation 30 and it can be seen that Amperes law is indeed invariant to the real orthogonal space time transformation.

(35) 0=x(1γHy+vϵPx+vϵEx)y(1γHxvϵPyvϵEy)(1γJz+vρ)ϵγ(tvz)Pzϵγ(tvz)Ez=1γHxx+vϵPxx+vϵExx1γHxy+vϵPyy+vϵEyy1γJzvρϵγPzt+ϵγvPzzϵγEzt+vϵγEzz=HyxHxxJzγvρ+vϵγPxx+vϵγPyy+vϵγExx+vϵγEyyϵγ2Pzx+ϵγ2vPzzϵγ2Ezt+vϵγ2Ezz=HyxHxyJzγvρ+vϵγ(Px+Ex)x+vϵγ(Py+Ey)yϵγ2Pztϵγ2Ezt+μϵγ2vz(Pz+Ez)

Using the definition for the electric flux density D = ϵ(P+E) and introducing the div (D) we get:

(36) 0=HyxHxyJzγvρ+vγdiv(D)vϵγ(Pz+Ez)zϵγ2Pztϵγ2Ezt+ϵγ2vz(Pz+Ez)=HyxHxyJz+γv(div(D)ρ)vϵγ2(z+vc2t)(Pz+Ez)ϵγ2Pztϵγ2Ezt+ϵγ2vz(Pz+Ez)

Now realizing that div (D)=ρ and canceling the z terms leads to:

(37) 0=HyxHxyJzϵγ2v2c2t(Pz+Ez)ϵγ2t(Pz+Ez)
=HyxHxyJzϵγ2(1+v2c2)t(Pz+Ez)

And finaly using γ2=11+v2c2 results in:

(38) 0=HyxHxyJzϵt(Pz+Ez)=HyxHxyJzϵPztϵEzt

From this derivation the next field transformation rules are derived:

(39) Ez=Ez
(40) Pz=Pz

Now the same method is applied on the law of induction.

(41) rot(E)=μMtμHt
(42) 0=rot(E)+μMt+μHt

The three seperated equations are as follows:

(43) 0=EzyEyz+μMxt+μHxt
(44) 0=ExzEzx+μMyt+μHyt
(45) 0=EyxExy+μMzt+μHzt

The transformation of the first equation 43 leads to:

(46) 0=Ezyγ(z+vc2t)Ey+μγ(tvz)(Mx+Hx)=Ezyzγ(Ey+vμ(Mx+Hx))+μγHxt+μtγ(MxvϵEy)

Transforming this with the corresponding rules of 34 results in:

(47) 0=EzyEyz+t(HxγvϵPyγvϵEy)+μtγ(MxvϵEy)
=EzyEyz+μtγ(Mxvϵ(Py+Ey))+Hxt

And the second equation 44 gives:

(48) 0=γ(z+vc2t)ExEzx+μγ(tvz)(My+Hy)=zγ(Exvμ(My+Hy))Ezx+μγHyt+μtγ(My+vϵEx)

Transforming this with the corresponding rules of 34 results in:

(49) 0=ExzEzx+μt(Hy+γvϵPx+γvϵEx)+μtγ(My+vϵEx)=ExzEzx+μtγ(My+vϵPx+2vϵEx)+μHyt

From these results it can be seen that the field transformation equations, derived from Amperes law, transform the first two equations of the law of induction, if the following additional transformation rules are introduced:

(50) Mx=γ(Mxvϵ(Py+2Ey))
(51) My=γ(Myvϵ(Px+2Ex))

Last but not least the transformation rules can be applied on equation 45.

(52) 0=x(1γEyvμ(Mx+Hx))y(1γEx+vμ(My+Hy))+γ(tvμz)(Mz+Hz)=EyxExyμγvx(Mx+Hx)μγvy(My+Hy)+γ2(tμvz)(Mz+Hz)

Using B=μ(M+H) and div (B)=0 leads to:

(53) 0=EyxExyγvdiv(B)+μγvz(Mz+Hz)+μγ2(tvz)(Mz+Hz)=EyxExy+μγ2v(z+vc2t)(Mz+Hz)
+μγ2(tvz)(Mz+Bz)

And finaly using γ2=11+v2c2 result in:

(54) 0=EyxExy+μγ2(1+v2c2)t(Mz+Hz)=EyxExy+μMzt+μHzt

This proves the invariance of the field equations to the real orthogonal spacetime. From this derivation the last field transformation rule has been found:

(55) Mz=Mz

In this proof it is shown that the following field transformation rules have to be applied:

(56) Ex=γ(Exvμ(My+Hy))Ey=γ(Eyvμ(Mx+Hx))Ez=Ez
(57) Hx=γ(Hxvϵ(Py+Ey))Hy=γ(Hyvϵ(Px+Ex))Hz=Hz

Both rules are identical to the transformation rules derived by A. Einstein [5], except that γ has changed. This is a very important result, as it shows that our well established electromagnetic theory is in fact invariant under orthogonal Cartesian space time transformation. However, in addition to the field transformation rules 56 and 57 additional transformation rules for the polarization and magnetization have to be considered.

(58) Px=γ(Px+vμ(My+2Hy))Py=γ(Pyvμ(Mx+2Hx))Pz=Pz
(59) Mx=γ(Mxvϵ(Py+2Ey))My=γ(My+vϵ(Px+2Ex))Mz=Mz

This is a new and very useful insight.

For arbitrary directions the transformations 56 to 59 have the form of:

(60) E=γ(E+v×μ(M+H))+(1γ)vE| v |2v
(61) H=γ(Hv×ϵ(P+E))+(1γ)vH v2v
(62) P=γ(Pv×μ(M+2H))+(1γ)vP| v |2v
(63) M=γ(M+v×ϵ(P+2E))+(1γ)vM v2v

To get the transformation rules for moved bulky bodies like Einstein defined in [24], we can add equations 61 and 63, or 62 and 60 to get the transformation rules for the flux densities.

(64) E=γ(E+v×B)+(1γ)vE| v |2v
(65) H=γ(Hv×D)+(1γ)vH v2v
(66) D=γ(Dvc2×H)+(1γ)vD| v |2v
(67) B=γ(B+vc2×E)+(1γ)vB v 2v

These equations are the same as for the Lorentz transformation based field transformation [24], except the γ' and the signs in front of the vc2 terms in the flux density transformation rules.

For velocities much smaller than c this simplifies to:

(68) EE+v×B
(69) HHv×D
(70) DD
(71) BB

which is in perfect agreement with our experience and valid for all technical applications.

4.3 The full set of wave equations

When polarization and magnetization can happen in empty space, wave equations for all quantities which are invariant under orthogonal Cartesian transformations are needed. This invariant wave operator is:

(72) 2x2+2y2+2z2+1c22t2

Applying the transformation rules for differentials 18 results in:

(73) 2x2+2y2+2z2+1c22t2=2x2+2y2+γ2(z+vc2z)(z+vc2t)+γ21c2(tvz)(tvz)=2x2+2y2+γ2(2z2+2vc22tz+v2c22t2)+γ2(1c22t22vc22zt+v2c22z2)=2x2+2y2+γ2(1+v2c2)2z2+γ2(1+v2c2)1c22t2=2x2+2y2+2z2+1c22t2

This proves the invariance of the new wave operator 72.

Applying the standard derivation for the wave equation like in [6, S. 414], with an additional curl on the field equations, leads to:

(74) Δ E 1 c 2 2 E t 2 = μ J t + 1 c 2 2 P t 2 + μ r o t M t + g r a d d i v E Δ H 1 c 2 2 H t 2 = 1 c 2 2 M t 2 r o t J ϵ r o t P t + g r a d d i v H

Introducing wave equations which are invariant to the orthogonal Cartesian space time transformation 74 gives:

(75) Δ E + 1 c 2 2 E t 2 = 0
(76) Δ H + 1 c 2 2 H t 2 = 0

To get the identities of 75 and 76 the following equations should be satisfied:

(77) 0=21c22Et2+μJt+1c22Pt2
+μrot(M)t+grad(div(E))0=21c22Ht2+1c22Mt2rot(J)ϵrot(P)t+grad(div(H))

which is

(78) rot(M)t=Jtϵ2Pt21μgrad(div(E))2ϵ2Et2
(79) rot(P)t=1ϵrot(J)+μ2Mt2+1ϵgrad(div(H))+2μ2Ht2

These equations give very interesting relations. Please note that they are only valid for oscillating fields where the field quantities fulfill the new wave equations. Using them for static cases will possibly lead to wrong results. Applying the curl on the first equation and replacing rot (P) and rot (J) by the second equation gives, with an integration with respect to time, the following non invariant wave equation for the magnetic moments:

(80) Δ M 1 c 2 2 M t 2 = g r a d d i v M + g r a d d i v H + 2 1 c 2 2 H t 2 + r o t 1 μ g r a d d i v E + 2 ϵ 2 E t 2 t

Applying the curl on the second equation and replacing rot(M)t by the first gives wave equations for the polarization and current densities:

(81) t Δ P 1 c 2 2 P t 2 + 1 ϵ Δ J 1 c 2 2 J t 2 = g r a d d i v 1 ϵ J + P t + g r a d d i v E t + 2 c 2 3 E t 3 r o t 1 ϵ r o t d i v H + 2 μ 2 H t 2

To get the corresponding invariant wave equations

(82) Δ M + 1 c 2 2 M t 2 = 0
(83) Δ P + 1 c 2 2 P t 2 = 0
(84) Δ J + 1 c 2 2 J t 2 = 0

the following equations should be satisfied:

(85) 0 = 2 1 c 2 2 M t 2 + g r a d d i v M + g r a d d i v H + 2 1 c 2 2 H t 2 + r o t 1 μ g r a d d i v E + 2 ϵ 2 E t 2 t 0 = 2 1 ϵ c 2 2 J t 2 + 2 1 c 2 3 P t 3 + g r a d d i v 1 ϵ J + P t + g r a d d i v E t + 2 c 2 3 E t 3 r o t 1 ϵ g r a d d i v H + 2 μ 2 H t 2

From div (B)=0 follows div (M)+div(H)=0. . Further-more:

(86) div(J)=ρt=div(D)t=ϵ(div(P)t+div(E)t)

Applying rot (div (. . .)) = 0 in addition, the remaining parts give:

(87) 0=21c22Mt2+21c22Ht2+rot(2ϵ2Et2)t0=21ϵc22Jt2+21c23Et3+2c23Et3rot(2μ2Ht2)

or

(88) 0=2c2t(rot(E)+μMt+μHt)0=2c22t2(rot(H)JϵPtϵEt)

Both identities are satisfied if Maxwell’s field equations 1 and 2 are satisfied. Now it has been proven that a set of solutions for the field equations satisfies wave equations for E,B,M,PandJ. This does not mean that every solution of a wave equation fulfills the field equations. However, the wave equations have to be fulfilled to find a solution for the field equations.

For these wave equations the simple plane wave solution is not possible. Two solutions for these equations have been found which are worth to mention. Assuming propagation in z direction, the first solution shows what happens if a plane of electromagnetic fields is applied:

(89) f(x,y,z,t)=sin(ωt)ekz

This solution explains the effect of tunneling. An oscillating field (ω = k · c) leads to a penetration of energy into empty space.

The second solution uses the perpendicular dimensions to compensate for the propagation direction and the time dimensions:

(90) f(x,y,z,t)=sin(k(zct))ek(| x |+| y |)

This solution is a traveling wave with finite field energy transport due to the exponential decay in the perpendicular directions.

4.4 Propagation of Rays of Light

Assume the following vector potential Av of a linearly polarized ray of light propagating in z direction:

(91) Av=(E0cksin(k(zct))ek(| x |+| y |)00)

This vector potential is not divergence free, resulting in the following non zero virtual charge density distribution:

(92) ρv=2E0kZ0cx| x |cos(k(zct))ek(| x |+| y |)

Figure 3 shows the virtual charge density distribution necessary to fulfill the conservation of charge rule. The virtual charges are divided on both sides of the polarization axes x. The movement of the charges is therefore perpendicular to the propagation direction and the wave is of transversal nature accordingly.

Figure 3 Illustration of the virtual charge density distribution of a linear polarized ray of light in the plane of propagation z and the polarisation axes x. Black illustrates maximal negative, gray zero and white maximal positive virtual charge density.
Figure 3

Illustration of the virtual charge density distribution of a linear polarized ray of light in the plane of propagation z and the polarisation axes x. Black illustrates maximal negative, gray zero and white maximal positive virtual charge density.

The corresponding fields still have some degrees of freedom which can be found by using inverse vector operators [25]. Field configurations have been chosen which lead to a sensible energy transport.

(93) E=E0(2cos(k(zct))xy| x || y |cos(k(zct))x| x |sin(k(zct)))ek(| x |+| y |)
(94) H=E0Z0(x| x |y| y |cos(k(zct))2cos(k(zct))y| y |sin(k(zct)))ek(| x |+| y |)

Figure 4 shows the norm of the field distribution in the plane perpendicular to direction of propagation, which in this case is z. It can be seen that there are two orthogonal preferred axes which are in alignment with the two field

Figure 4 Illustration of the field distribution in the plane perpendicular to direction of propagation z (norm of 

 E→ or H→ ).
$\left. \vec{E}\,\text{or}\,\vec{H} \right).$White is maximal intensity and black is zero intensity.
Figure 4

Illustration of the field distribution in the plane perpendicular to direction of propagation z (norm of EorH ). White is maximal intensity and black is zero intensity.

directions EandB in the linearly polarized example. It can be seen that this field configuration has finite energy in the x-y-plane.

From the equation of continuity, which ensures the conservation of charge, the density of polarization currents can be derived with the assumption that the free charge density is zero. This gives the solution corresponding to Amperes law if the field quantities 93 and 94 are used.

(95) Jp=E0kZ0(1xy| x || y |0)sin(k(zct))ek(| x |+| y |)

The density of polarization currents shows 90 phase shift corresponding to the charge density of equation 92. This seems correct as the maximum velocity of the virtual charges corresponds to the minimum displacement and vice versa. Figure 5 shows the densities of polarization currents in the yx-plane at t = T/4 and z = 0. Keeping in mind that positive and negative virtual charges are involved and therefore negative virtual charges change the direction of the density of polarization current vectors with respect to the velocity of the moving virtual charge, it can be seen that the positive charges in the lower half space will move upwards (+x) and the negative charges in the upper half space will move downwards (−x) with 90 phase shift. At the same time the virtual charges will concentrate in the center of the y axes. The whole movement will be a continuous alternation of spreading the virtual charges in space and concentrate them back into the center.

Figure 5 Illustration of the density of polarization currents in y-x plane at t = T/4 and z = 0.
Figure 5

Illustration of the density of polarization currents in y-x plane at t = T/4 and z = 0.

The solution corresponding to 93 and 94 for M is:

(96) M=E0c(x| x |y| y |cos(k(zct))cos(k(zct))0)ek(| x |+| y |)

Checking this set of solutions for conservation laws it has been found that div (P)+div(E)=0 and div (M) + div (H) = 0 holds for every point in space.

Figure 6 shows the magnetic momenta caused by overlapping spin momenta of the involved virtual charge densities. Please note that alignment within the fields is possible, more or less independent for currents and spins.

Figure 6 Illustration of the magnetic momenta necessary to satisfy the ray equations
Figure 6

Illustration of the magnetic momenta necessary to satisfy the ray equations

This derivation just shows what the ray of light solution could look like. This is no rigorously derived proof that it will indeed happen in this way. This proof would need a clear concept about virtual charges and how they behave. This would definetely go beyond the aim of this paper.

In contrast to the plane wave solution the proposed ray of light solution needs field components in direction of propagation to realize a finite but divergence free configuration. It is clear that finite fields can only be pure vortex fields if two space dimensions are involved. Figure 7 shows the vertices of the flux densities. It is easy to prove that these fields are divergence free and fulfill the new ray equations.

Figure 7 Illustration of the flux density vortices in direction of propagation z. The 

D→ and B→
$\vec{D}\ \text{and}\ \vec{B}$fields look the same, however the shown planes are x-z for 

D→
$\vec{D}$and y-z for 

B→
$\vec{B}$
Figure 7

Illustration of the flux density vortices in direction of propagation z. The DandB fields look the same, however the shown planes are x-z for D and y-z for B

The flux densities are:

(97) D=ϵE0(cos(k(zct))0x| x |sin(k(zct)))ek(| x |+| y |)
(98) B=E0c(0cos(k(zct))y| y |sin(k(zct)))ek(| x |+| y |)

It can be seen that the flux densities are in phase but on orthogonal planes. Therefore a ray of light is a TEM wave. Finally, the ray of light looks like expected for electromagnetic waves, e.g. [6, p. 35, fig. 1.26] or [7, p. 33, fig. 1.26].

4.5 Poynting vector and energy transport

Defining the Poynting vector as S=E×H and appling the divergence div () and Maxwell’s equations leads to [6, p. 109], [7, p. 108]:

(99) div(S)=div(E×H)=Hrot(E)Erot(H)=HBtEDtEJ

This Poynting vector is usually interpreted as the energy flow of the electromagnetic field as equation 99 shows that the source of the energy flowing out of a volume is a change of energy within the volume or an ohmic source (or loss) [6, p. 110]. However, in the presence of virtual charges 99 is not equal to a change of field energy. To come up with a better solution a new definition of the Poynting vector is needed.

Lets start with div (E×M). Replacing rot (E) by 2 and rot (M) by 79 simplifying by setting J=0 and grad (div(E))=0 gives:

(100) div(E×M)=Mrot(E)Erot(M)=M(μM+μH)t+E(ϵP+2ϵE)t

Now div (P×H) is calculated. Replacing rot (H) by 1 and rot (P) by 80 simplifying by setting J=0 and grad (div (H)) = 0 gives:

(101) div(P×H)
=Hrot(P)Prot(H)=H(μM+2μH)tP(ϵP+ϵE)t

Finally div (P×M) is calculated. Replacing rot (P) by 80 and rot (M) by 79 with same simplifications like above gives:

(102) div(P×M)=Mrot(P)Prot(M)=M(μM+2μH)t+P(ϵP+2ϵE)t

Redefining the pointing vector in the following form

(103) S=12P×H12E×M12P×M=12E×H12c2D×B

the divergence of the vector relates directly to the change of field energy:

(104) div(12E×H12c2D×B)=12μ(M+H)Ht12μ(M+H)tH12ϵ(P+E)Et12ϵ(P+E)tE=12Hμ(M+H)t12Eϵ(P+E)t

The sign and the 1/2 are chosen to ensure the same inherent energy definition like in 99, but 104 enables arbitrary orientations of the field and polarization vectors.

The newly defined Poynting vector for the ray of light results in:

(105) S=E02Z0(00cos(k(zct))2)e2k(| x |+| y |)

Please note that the time average of the divergence of this vector is zero because on every point the energy flows in and out.

This new Poynting vector can be used as an energy flow of electromagnetic field waves as it makes sense for many applications. However, the Poynting vector is not the only possible way to transfer energy. This is possible for magnetic or electric coupling without any Poynting vector in space.

5 Conclusion

In this article the invariance of Maxwell’s field equations to an orthogonal Cartesian space time transformation has been proven if polarization and magnetization are allowed to happen in space. The ray of light solution acts as a prove of usefulness of this new approach. This solution gives a detailed description in terms of electromagnetic fields of what a ray of light could look like. This is in fact not possible with the plane wave solution which shows infinite field energy transport due to the unbounded fields perpendicular to the propagation direction. It is theoretically unsatisfactory that the theory of relativity gives local field solutions a prominent meaning as carrier of energy density and momentum but the plane wave solution carries infinite total field energy.

Possible applications of the bounded solution are diffraction theories and ray tracing. Even the effect of tunneling appears to be a field solution in the Euclidean space. The new role of the magnetization and polarization, used in this derivation, can help to find new description models of material properties on a microscopic scale.

Future work should be dedicated to compare the orthogonal Cartesian relativity with the Lorentz transformation based relativity as the Lorentz transformation acts as a standard reference to prove new theories. In this context it was important to show that both geometrical spaces give valid solutions to Maxwell’s field equations and the principle of relativity. Finally, deducing the constance of velocity of light within the Euclidean space would be of great theoretical value.

References

[1] Maxwell J. C., A dynamical theory of the electromagnetic field, Philos. Trans. Royal Soc., 1865, 155, 459–512.10.1098/rstl.1865.0008Search in Google Scholar

[2] Maxwell J. C., A Treatise on Electricity and Magnetism, Macmillan, New York, 1873Search in Google Scholar

[3] Lorentz H. A., Versuch einer Theorie der elektrischen und magnetischen Erscheinungen in bewegten Körpern., E. J. Brill, Leiden, 1895Search in Google Scholar

[4] Pauli W., Theory of Relativity, Pergamon Press, Oxford England, 1958Search in Google Scholar

[5] Einstein A., Zur Elektrodynamik bewegter Körper, Annalen der Physik, 1905, 17, 891–921.10.1002/andp.19053221004Search in Google Scholar

[6] Lehner G., Elektromagnetische Feldtheorie, 6th ed., Springer-Verlag, Berlin Heidelberg, 2008Search in Google Scholar

[7] Lehner G, Electromagnetic field theory for engineers and physicists, Springer-Verlag, Berlin Heidelberg, 201010.1007/978-3-540-76306-2Search in Google Scholar

[8] Jackson J. D., Klassische Elektrodynamik, 4th ed. ,Walter de Gruyter, Berlin, 200610.1515/9783110200034Search in Google Scholar

[9] W. Ritz, Recherges critiques sur l’electrodynamique generale, Ann. Chim. Phys., 1908, 13, 145–275.Search in Google Scholar

[10] Lämmerzahl C., Special relativity and Lorentz invariance, Annalen der Physik, 2005, 14, 71–102.10.1002/andp.200410127Search in Google Scholar

[11] Heintze J., Bock P., Lehrbuch zur Experimentalphysik, Springer Spektrum, Berlin Heidelberg, 201410.1007/978-3-642-41210-3Search in Google Scholar

[12] Stratton J. A., Electromagnetic Theory, reissue ed., John Wiley and Sons Inc., Hoboken New Jersey, 2007Search in Google Scholar

[13] Ignatowsky W., Das Relativitätsprinzip, Archiv der Mathematik und Physik, 1910, 17, 1-24.Search in Google Scholar

[14] Ignatowsky W., Das Relativitätsprinzip, Archiv der Mathematik und Physik, 1911, 18, 17-40.Search in Google Scholar

[15] Frank W., Rothe H., Über die Transformation der Raumzeitkoordinaten von ruhenden auf bewegte Systeme, Annalen der Physik, 1911, 34, 825-855.10.1002/andp.19113390502Search in Google Scholar

[16] Bartelmann M., Feuerbacher B., Krüger T., Lüst D., Rebhan A., Wipf A., Theoretische Physik, Springer, Berlin Heidelberg, 201510.1007/978-3-642-54618-1Search in Google Scholar

[17] Dyson F. J., The radiation theories of Tomonaga, Schwinger and Feynman, Phys. Rev., 1949, 75, 486-502.10.1103/PhysRev.75.486Search in Google Scholar

[18] Feynman R. P., QED Die seltsame Theorie des Lichts und der Materie, Piper, München Zürich, 2006Search in Google Scholar

[19] Berestetzki W. B., Lifschitz E. M., Pitajewski L. P., Quantenelektrodynamik, 7th ed., Verlag Harri Deutsch, Frankfurt am Main, 2009Search in Google Scholar

[20] Griflths D. J., Elektrodynamik, 3rd ed., Pearson, München, 2011Search in Google Scholar

[21] Tchernyshov O.,Magnetism: Freedom for the poles, Nature, 2008, 451, 22–23.10.1038/451022bSearch in Google Scholar PubMed

[22] Castelnovo C., Moessner R., Sondhi S. L., Magnetic monopoles in spin ice, Nature, 2008, 451, 42–45.10.1038/nature06433Search in Google Scholar PubMed

[23] Gingras M. J. P., Observing monopoles in magnetic analog of ice, Science, 2009, 326, 375–376.10.1126/science.1181510Search in Google Scholar PubMed

[24] Einstein A., Laub J., Über die elektromagnetischen Grundgleichungen für bewegte Körper, Annalen der Physik, 1908, 26, 532–540.10.1002/andp.19083310806Search in Google Scholar

[25] Sahoo S., Inverse Vector Operations, arXiv:0804.2239 math-ph, 2010Search in Google Scholar

Received: 2019-05-07
Accepted: 2019-10-29
Published Online: 2019-12-30

© 2019 J. Schliewe, published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Regular Articles
  2. Non-equilibrium Phase Transitions in 2D Small-World Networks: Competing Dynamics
  3. Harmonic waves solution in dual-phase-lag magneto-thermoelasticity
  4. Multiplicative topological indices of honeycomb derived networks
  5. Zagreb Polynomials and redefined Zagreb indices of nanostar dendrimers
  6. Solar concentrators manufacture and automation
  7. Idea of multi cohesive areas - foundation, current status and perspective
  8. Derivation method of numerous dynamics in the Special Theory of Relativity
  9. An application of Nwogu’s Boussinesq model to analyze the head-on collision process between hydroelastic solitary waves
  10. Competing Risks Model with Partially Step-Stress Accelerate Life Tests in Analyses Lifetime Chen Data under Type-II Censoring Scheme
  11. Group velocity mismatch at ultrashort electromagnetic pulse propagation in nonlinear metamaterials
  12. Investigating the impact of dissolved natural gas on the flow characteristics of multicomponent fluid in pipelines
  13. Analysis of impact load on tubing and shock absorption during perforating
  14. Energy characteristics of a nonlinear layer at resonant frequencies of wave scattering and generation
  15. Ion charge separation with new generation of nuclear emulsion films
  16. On the influence of water on fragmentation of the amino acid L-threonine
  17. Formulation of heat conduction and thermal conductivity of metals
  18. Displacement Reliability Analysis of Submerged Multi-body Structure’s Floating Body for Connection Gaps
  19. Deposits of iron oxides in the human globus pallidus
  20. Integrability, exact solutions and nonlinear dynamics of a nonisospectral integral-differential system
  21. Bounds for partition dimension of M-wheels
  22. Visual Analysis of Cylindrically Polarized Light Beams’ Focal Characteristics by Path Integral
  23. Analysis of repulsive central universal force field on solar and galactic dynamics
  24. Solitary Wave Solution of Nonlinear PDEs Arising in Mathematical Physics
  25. Understanding quantum mechanics: a review and synthesis in precise language
  26. Plane Wave Reflection in a Compressible Half Space with Initial Stress
  27. Evaluation of the realism of a full-color reflection H2 analog hologram recorded on ultra-fine-grain silver-halide material
  28. Graph cutting and its application to biological data
  29. Time fractional modified KdV-type equations: Lie symmetries, exact solutions and conservation laws
  30. Exact solutions of equal-width equation and its conservation laws
  31. MHD and Slip Effect on Two-immiscible Third Grade Fluid on Thin Film Flow over a Vertical Moving Belt
  32. Vibration Analysis of a Three-Layered FGM Cylindrical Shell Including the Effect Of Ring Support
  33. Hybrid censoring samples in assessment the lifetime performance index of Chen distributed products
  34. Study on the law of coal resistivity variation in the process of gas adsorption/desorption
  35. Mapping of Lineament Structures from Aeromagnetic and Landsat Data Over Ankpa Area of Lower Benue Trough, Nigeria
  36. Beta Generalized Exponentiated Frechet Distribution with Applications
  37. INS/gravity gradient aided navigation based on gravitation field particle filter
  38. Electrodynamics in Euclidean Space Time Geometries
  39. Dynamics and Wear Analysis of Hydraulic Turbines in Solid-liquid Two-phase Flow
  40. On Numerical Solution Of The Time Fractional Advection-Diffusion Equation Involving Atangana-Baleanu-Caputo Derivative
  41. New Complex Solutions to the Nonlinear Electrical Transmission Line Model
  42. The effects of quantum spectrum of 4 + n-dimensional water around a DNA on pure water in four dimensional universe
  43. Quantum Phase Estimation Algorithm for Finding Polynomial Roots
  44. Vibration Equation of Fractional Order Describing Viscoelasticity and Viscous Inertia
  45. The Errors Recognition and Compensation for the Numerical Control Machine Tools Based on Laser Testing Technology
  46. Evaluation and Decision Making of Organization Quality Specific Immunity Based on MGDM-IPLAO Method
  47. Key Frame Extraction of Multi-Resolution Remote Sensing Images Under Quality Constraint
  48. Influences of Contact Force towards Dressing Contiguous Sense of Linen Clothing
  49. Modeling and optimization of urban rail transit scheduling with adaptive fruit fly optimization algorithm
  50. The pseudo-limit problem existing in electromagnetic radiation transmission and its mathematical physics principle analysis
  51. Chaos synchronization of fractional–order discrete–time systems with different dimensions using two scaling matrices
  52. Stress Characteristics and Overload Failure Analysis of Cemented Sand and Gravel Dam in Naheng Reservoir
  53. A Big Data Analysis Method Based on Modified Collaborative Filtering Recommendation Algorithms
  54. Semi-supervised Classification Based Mixed Sampling for Imbalanced Data
  55. The Influence of Trading Volume, Market Trend, and Monetary Policy on Characteristics of the Chinese Stock Exchange: An Econophysics Perspective
  56. Estimation of sand water content using GPR combined time-frequency analysis in the Ordos Basin, China
  57. Special Issue Applications of Nonlinear Dynamics
  58. Discrete approximate iterative method for fuzzy investment portfolio based on transaction cost threshold constraint
  59. Multi-objective performance optimization of ORC cycle based on improved ant colony algorithm
  60. Information retrieval algorithm of industrial cluster based on vector space
  61. Parametric model updating with frequency and MAC combined objective function of port crane structure based on operational modal analysis
  62. Evacuation simulation of different flow ratios in low-density state
  63. A pointer location algorithm for computer visionbased automatic reading recognition of pointer gauges
  64. A cloud computing separation model based on information flow
  65. Optimizing model and algorithm for railway freight loading problem
  66. Denoising data acquisition algorithm for array pixelated CdZnTe nuclear detector
  67. Radiation effects of nuclear physics rays on hepatoma cells
  68. Special issue: XXVth Symposium on Electromagnetic Phenomena in Nonlinear Circuits (EPNC2018)
  69. A study on numerical integration methods for rendering atmospheric scattering phenomenon
  70. Wave propagation time optimization for geodesic distances calculation using the Heat Method
  71. Analysis of electricity generation efficiency in photovoltaic building systems made of HIT-IBC cells for multi-family residential buildings
  72. A structural quality evaluation model for three-dimensional simulations
  73. WiFi Electromagnetic Field Modelling for Indoor Localization
  74. Modeling Human Pupil Dilation to Decouple the Pupillary Light Reflex
  75. Principal Component Analysis based on data characteristics for dimensionality reduction of ECG recordings in arrhythmia classification
  76. Blinking Extraction in Eye gaze System for Stereoscopy Movies
  77. Optimization of screen-space directional occlusion algorithms
  78. Heuristic based real-time hybrid rendering with the use of rasterization and ray tracing method
  79. Review of muscle modelling methods from the point of view of motion biomechanics with particular emphasis on the shoulder
  80. The use of segmented-shifted grain-oriented sheets in magnetic circuits of small AC motors
  81. High Temperature Permanent Magnet Synchronous Machine Analysis of Thermal Field
  82. Inverse approach for concentrated winding surface permanent magnet synchronous machines noiseless design
  83. An enameled wire with a semi-conductive layer: A solution for a better distibution of the voltage stresses in motor windings
  84. High temperature machines: topologies and preliminary design
  85. Aging monitoring of electrical machines using winding high frequency equivalent circuits
  86. Design of inorganic coils for high temperature electrical machines
  87. A New Concept for Deeper Integration of Converters and Drives in Electrical Machines: Simulation and Experimental Investigations
  88. Special Issue on Energetic Materials and Processes
  89. Investigations into the mechanisms of electrohydrodynamic instability in free surface electrospinning
  90. Effect of Pressure Distribution on the Energy Dissipation of Lap Joints under Equal Pre-tension Force
  91. Research on microstructure and forming mechanism of TiC/1Cr12Ni3Mo2V composite based on laser solid forming
  92. Crystallization of Nano-TiO2 Films based on Glass Fiber Fabric Substrate and Its Impact on Catalytic Performance
  93. Effect of Adding Rare Earth Elements Er and Gd on the Corrosion Residual Strength of Magnesium Alloy
  94. Closed-die Forging Technology and Numerical Simulation of Aluminum Alloy Connecting Rod
  95. Numerical Simulation and Experimental Research on Material Parameters Solution and Shape Control of Sandwich Panels with Aluminum Honeycomb
  96. Research and Analysis of the Effect of Heat Treatment on Damping Properties of Ductile Iron
  97. Effect of austenitising heat treatment on microstructure and properties of a nitrogen bearing martensitic stainless steel
  98. Special Issue on Fundamental Physics of Thermal Transports and Energy Conversions
  99. Numerical simulation of welding distortions in large structures with a simplified engineering approach
  100. Investigation on the effect of electrode tip on formation of metal droplets and temperature profile in a vibrating electrode electroslag remelting process
  101. Effect of North Wall Materials on the Thermal Environment in Chinese Solar Greenhouse (Part A: Experimental Researches)
  102. Three-dimensional optimal design of a cooled turbine considering the coolant-requirement change
  103. Theoretical analysis of particle size re-distribution due to Ostwald ripening in the fuel cell catalyst layer
  104. Effect of phase change materials on heat dissipation of a multiple heat source system
  105. Wetting properties and performance of modified composite collectors in a membrane-based wet electrostatic precipitator
  106. Implementation of the Semi Empirical Kinetic Soot Model Within Chemistry Tabulation Framework for Efficient Emissions Predictions in Diesel Engines
  107. Comparison and analyses of two thermal performance evaluation models for a public building
  108. A Novel Evaluation Method For Particle Deposition Measurement
  109. Effect of the two-phase hybrid mode of effervescent atomizer on the atomization characteristics
  110. Erratum
  111. Integrability analysis of the partial differential equation describing the classical bond-pricing model of mathematical finance
  112. Erratum to: Energy converting layers for thin-film flexible photovoltaic structures
Downloaded on 21.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/phys-2019-0077/html
Scroll to top button