Relativistic Rotating Electromagnetic Fields

Departamento de Ciencias Exactas y Tecnologia, Centro Universitario de los Lagos, Universidad de Guadalajara, Enrique Daz de León 1144, Col. Paseos de la Montaña, Lagos de Moreno, Jalisco, Mexico IPICyT, Instituto Potosino de Investigación Científica y Tecnológica, Camino a la Presa San José 2055, Col. Lomas 4a Sección, 78216 San Luis Potosi S.L.P., Mexico Centro de Investigación y de Estudios Avanzados del Instituto Politecnico Nacional, Av. Instituto Politecnico Nacional 2508, C.P. 07360 Alcaldía Gustavo A. Madero, Ciudad de Mexico, Mexico


Introduction
The Poynting vector can be used to establish the state of motion of the electromagnetic field, for example, consider a charged line distribution on the x-axis, moving along this axis, while an observer at rest measures a Poynting vector, a comoving observer measures none; hence, the velocity of the comoving observer can be identified with the velocity of the electromagnetic field. Since the electromagnetic field is a distributed object, a full family of observers measuring no Poynting vector is required to fully characterize the electromagnetic field's motion.
In relativity, a reference frame is an idealization of an observer equipped with some measuring devices in some motion state; consequently, a set of reference frames is described by a time-like congruence, the world lines of a family of observers. In its turn, this congruence is described by a family of unit time-like vectors, the observers' 4-velocity field, τ.
Covariant electromagnetic fields can be classified according to their first invariant as of magnetic, electric, and null types. The second invariant generates the pure and nonpure classes. Except for pure null fields, it is well known that it is always possible to find the reference frame where the Poynting vector vanishes [1]. Mitskievich showed how to find such a reference frame; he proved that, for any nonnull pure field, it is possible to construct a simple bivector with the electromagnetic field tensor; from it, or its Hodge dual-conjugate, a unit time-like vector can be extracted. For pure fields, this vector is used to describe the reference frame in which there is only an electric field or a magnetic field, while for nonpure fields, it describes the reference frame in which both the electric and magnetic vectors are parallel to each other [2,3]. It follows that generically electromagnetic fields can have their own dynamical state of motion that may entail a plethora of physical effects and applications.
However, rotation and relativity are hard to match with each other [4] because of some puzzling features. In astrophysics, there exists the so-called light cylinder problem, the lack of understanding about what happens in the region where the magnetic field lines of a rigidly rotating pulsar reach the speed of light [4][5][6]; some authors argue that the magnetic field lines rotate with a speed greater than that of the light beyond this cylinder [7]. The electromagnetic fields have also opened fundamental questions in terms of rotational motion [8,9], and this theory provides a strong support to model several astrophysical objects, such as pulsars [10,11] and active galactic nuclei [12].
In this paper, we consider axially symmetric electromagnetic fields in rotational motion around the symmetry axis. They are characterized not only by the rotation of the reference frame but also by the fact that such fields have an angular momentum density in any other reference frame different from the one comoving with the field. Indeed, this angular momentum has been observed long ago [13]. In particular, if the reference frame is at rest with respect to the rotation axis, the Poynting vector of these rotational fields will be tangent to the closed lines around the rotation axis. According to Jackson [14], one can always associate an angular momentum density to these vectors. An equivalent approach to determine the angular momentum density is based on Nöther's theorem [15,16].
Reference frames are introduced here by means of both a rotational and a nonrotational time-like congruence in the Minkowski spacetime. The velocity four-vector, τ (the monad field), associated to them is then used to perform the 3 + 1 splitting of any four-dimensional vector and in the splitting of the field tensor into its corresponding electric and magnetic fields. The explicit general covariance of the 3 + 1 splitting formalism allows the usage of abstract representations of the tensor quantities. It is applied in special relativity when noninertial effects are evaluated without involving any a priori assumption. Therefore, it is also a natural framework in the theory of general relativity.
Hence, two families of observers will be considered herein, namely, the inertial family, at rest with respect to the rotation axis, and the noninertial family in rotational motion with respect to the same axis. Observers in the rotating family measure no Poynting vector and are considered comoving with the electromagnetic field. A causal border known as the light surface appears whenever the electromagnetic field becomes of the pure null type. This light surface divides the spacetime in at least two different regions; in these cases, the noninertial family may split into two complementary subsets, one for each region. With pure electromagnetic fields, these subfamilies are labeled by the indexes fI, IIg, while for nonpure fields by fA, Bg. No label is used for iner-tial observers. The two families of observers can be conveniently discussed by using the theory of arbitrary reference frames, also known as the 3 + 1 decomposition, orthogonal splitting, and/or the old monad formalism [15,17,18].
The rest of the paper is organized as follows. Section 2 presents the electrodynamics in arbitrary reference frames, the classification of electromagnetic fields in terms of its invariants, and the propagation of electromagnetic fields from the point of view of the electromagnetic field tensor's invariants. Section 3 introduces stationary electromagnetic fields with angular momentum; they are discussed from the point of view of the inertial observers. Section 4 deals with pure rotating electromagnetic fields described from the point of view of both the inertial and the comoving rotating observers; the velocity fields of the last ones are determined with respect to the inertial observers. Section 5 presents some examples regarding the pure rotating fields of point magnetic dipoles, which are generically important models of magnetospheres in both astrophysics and geophysics [19][20][21][22][23]. In Section 6, the corresponding 4-velocity field of comoving observers with nonpure electromagnetic fields is presented. Section 7 shows an example of an impure field, the superposition of an electric Coulomb field, and the magnetic field of a point dipole. Section 8 presents the conclusions. There are also two Appendixes, one on the basic definition of the Cartan formalism and the theory of arbitrary reference frames and a second one which gives the expressions of the charge and current densities in a rotating reference frame.
Here, we use the signature ð+, − , − , − Þ and a system of units in which c = 1 (unless explicitly shown otherwise for convenience). Greek indices are taken to run from 0 to 3 and Latin indices from 1 to 3, and we adopt the standard convention for summation over repeated indices. Furthermore, we will indicate three vectors with bold symbols.

Electromagnetic Fields in Arbitrary Reference Frames
In general, field tensors F in arbitrary reference frames are written as 2-form: On the other hand, the electromagnetic vector potential is written as a covector, A = A μ ðxÞdx μ , and the connection between F and A is given by the relationship: With respect to a given reference frame, the field tensor splits into two terms: involving the electric and magnetic covectors [15]. The * operation denotes the Hodge dual conjugation defined in Equation (A.3).

Advances in High Energy Physics
From a given electromagnetic field tensor, the corresponding electric and magnetic covectors in the reference frame represented by the monad field τ can be calculated from the following: The Poynting vector can be written as a Hodge conjugate as well as the following: See (A. 10). An arbitrary vector can be decomposed into the sum of a vector parallel to τ, the time-like component, and a vector perpendicular to τ, which is the space-like component. For example, as well known, the four-vector of the electric current can be written in terms of the four-dimensional velocity u μ and decomposed as follows: where ρ ðτÞ and j μð3Þ are the charge and current densities measured in the reference frame represented by the monad field τ , respectively, while b μν is the four-dimensional projector; b = g − τ ⊗ τ is also used to define the three-dimensional scalar product of the differential forms denoted by • in the following (see Appendix A).
The electromagnetic splitting (Equation (3)), with Equation (4), is the consequence of the component form of the Lorentz force [15]: Here, the three-dimensional velocity v of the charged particle acted by the Lorentz force follows from the general definition: where

Classification of Electromagnetic Fields and Its
Propagation. The classification of electromagnetic fields is based on their invariants (see [1,3]): where The first invariant suggests that the fields can be classified according to their sign in the following kinds: (1) magnetic type if I 1 > 0, (2) electric type if I 1 < 0, and (3) null type if I 1 = 0. On the other hand, the second invariant allows to introduce an additional subclassification in (a) pure if I 2 = 0 and (b) nonpure if I 2 ≠ 0.
Some useful identities are as follows: We move now to the deduction of an important result [3]; connecting the propagation speed of the electromagnetic field with the electromagnetic field invariants is presented. Consider the electromagnetic energy momentum tensor as follows: in Gaussian units. When it is contracted with an arbitrary monad, it contains the electromagnetic energy density and the Poynting vector in the corresponding reference frame: and the squared expression is calculated using Equation (11): These constructions are not only scalars under coordinate transformations, but they are also independent of the choice of reference frame since the right hand side does not depend on the monad.
Landau and Lifshitz [1] have shown that the propagation speed of the electromagnetic field is given by the following: 3 Advances in High Energy Physics Using Equation (15), Mitskievich has shown that the propagation speed of electromagnetic fields satisfies the following: Hence, only null pure fields do propagate with the speed of light. Equation (16) is valid for observers with arbitrary motion on the background of curved spacetime; it coincides with Equation (15) when inertial observers in flat spacetime are considered.
Taking into account that the divergence and curl operators in an arbitrary reference frame have the forms div V = − * d * V − G•E and curl V = * ðτ ∧ dVÞ, the 3 + 1 splitted form of Maxwell's equations is as follows [15]: and £ τ is the Lie derivative with respect to τ [24,25]. These are the electromagnetic equations which should be used by the observers in motion defined by the τ congruence. The additional terms in Maxwell's equations are interpreted as the charge and current densities (both electric and magnetic ones) of kinematic or fictitious nature, the latter by analogy with the fictitious forces which appear in noninertial frames [15]. When the covariant Maxwell equations are considered in their orthogonal 3 + 1 splitting form, the following remarks can be brought out: (i) From Gauss' law (Equation (19)), it follows that for a pure magnetic field, (E = 0), an electric charge density should always exist in order to compensate the kinematic (fictitious) charge induced by the rotation of the magnetic field (ii) From Ampere's law (Equation (20)), it follows that for a pure electric field, (B = 0), there should always exist an electric current density to compensate the kinematic (fictitious) current induced by the defor-mation of the reference frame and/or the nonstationarity of the electric field (iii) For pure electric fields, the law of absence of monopoles (Equation (21)) requires that the electric field should be orthogonal to the angular velocity ω. There should be no kinematic monopole densities (iv) On the other hand, Faraday's induction law (Equation (22)) for pure magnetic fields requires the nonexistence of kinematic magnetic current densities (v) In both cases, Equations (21) and (22) as originating from the constitutive Maxwell equations will be satisfied whenever the field tensor is an exact form, F = dA, because dF = 0 and ddA ≡ 0

Electromagnetic Fields with Angular Momentum
The four-potential of a stationary axially symmetric field is as follows: The corresponding electromagnetic field has the following expression: and the electric and magnetic fields (Equation (4)), with respect to an inertial reference frame represented by τ = θ ð0Þ , are as follows: The Poynting covector, in this inertial reference frame, is as follows: It can be seen that it points in the azimuthal direction; hence, according to Jackson [14], the fields (Equation (23)) possess an angular momentum density.
From Maxwell's Equation (17), the 4-current covector is readily calculated, j = ρθ ð0Þ + Jθ ð3Þ , where the charge density measured in the inertial reference frame is as follows: Advances in High Energy Physics while the corresponding current density is as follows: The corresponding expressions for a rotating reference frame can be easily found by inserting Equations (28) and (29) into Equation (B.4).
For the fields under consideration (Equation (23)), the electromagnetic invariants are as follows: Electromagnetic fields of many different types may be considered according to the choices of the functions Mðr, ϑ Þ and Nðr, ϑÞ. Pure electromagnetic rotating fields are obtained when either M = MðNÞ or N = NðMÞ. The first case corresponds to pure magnetic rotating fields, while the second to a pure electric rotating fields, but only the first case is considered here, since the second one is quite similar. The cases M = 0 or N = 0 result in fields with vanishing Poynting vector. Impure rotating fields correspond to other choices of M and N.

Pure Rotating Electromagnetic Fields
We consider the special case when M = MðNðr, ϑÞÞ. In this case, the electromagnetic field tensor is a simple bivector: where we have introduced the notation: For the special case considered here, it is possible to see from Equations (25) and (26) that in the inertial reference frame, several electric fields, can be associated to the same magnetic field: Of course, in the alternative case, N = NðMÞ, several magnetic fields can be associated to the same electric field. The fields (Equation (33)) are unipolar-induced electric fields by the rotation of the magnetic field. As it would be seen below, the derivative dM/dN is interpreted as the angular speed of the magnetic field.
The first electromagnetic field invariant becomes the following: For |ξ | <1, it is of pure magnetic type, while for |ξ | = 1, it is of pure null type; if |ξ | >1, it is of the pure electric type. In this sense, ξ from Equation (32) divides Minkowski spacetime into different regions.
In the spacetime regions where |ξ | <1, we can choose the monad field for the vanishing Poynting observers as follows: Comparing Equation (36) with Equation (A.14), one finds the velocity field of the comoving observers with the electromagnetic field in the inertial reference frame given by the following: Consequently, as it was mentioned, the derivative dM/d N may be interpreted as the angular speed of the observers I and therefore as that of the magnetic field. It is as follows: The electromagnetic field tensor (Equation (31)) in terms of the monad field (Equation (36)) is as follows: In the frame (Equation (36)), the electric field vanishes and only the magnetic field vector remains and can be obtained from comparing Equation (39) with Equation (3): while the Poynting vector vanishes in this reference frame.
For the spacetime regions where |ξ | >1, we can choose the monad field for the vanishing Poynting observers as the following:

Advances in High Energy Physics
Comparing Equation (41) with Equation (A.14), one can find the velocity field of the comoving observers with the electromagnetic field in the inertial reference frame: The electromagnetic field tensor (Equation (31)) in terms of the monad field (Equation (41)) is as follows: In the frame represented by the monad field (Equation (41)), the magnetic field vanishes and the surviving electric field vector can be obtained by comparing Equation (43) with Equation (3): and the Poynting vector vanishes in this reference frame.
The two congruences, defined by the monad fields (Equations (36) and (41)), are separated by the light surface, which is the surface where the first electromagnetic invariant vanishes. On this surface, the electromagnetic field rotates with the speed of light; the field tensor (Equation (31)) is the exterior product of a null covector and a spatial covector. Consequently, it is not possible to find a reference frame rotating in synchrony with the electromagnetic field, as expected from Mitskievich's Equation (16).

Point Magnetic Dipoles in an Arbitrary
Rotation State Rotating point magnetic dipoles are special cases of pure electromagnetic fields. They provide common models of the magnetospheres of planets and stars [19,22]. In order to consider them, the following function is used: In the inertial reference frame at rest with the rotation axis, the electric and magnetic covectors, the Poynting covector, the charge density, and the three-dimensional current density, are found, inserting Equation (45) in Equations (25), (26), (27), (28) and (29): The magnetic field corresponds to that of a point dipole of magnitude |k | aligned with the z-axis. But for each choice of MðNÞ, there is a different electric field, though they differ up to a conformal factor. The corresponding Poynting vectors differ also in the same manner. As seen before, the derivative dM/dN is connected with the rotation of the electromagnetic field (Equation (38)). The charge density structure is now dictated by this derivative, giving rise to a vast number of possibilities; see some examples below. The function MðNÞ, or its derivative, has to be chosen from physical considerations or from observations; for example, from an observed charge density, it should be possible to reconstruct the derivative by solving the differential Equation (49).
In this case, the first invariant becomes the following: with ξ given by Equation (32). Notice that on the light surface |ξ | = 1, the first invariant vanishes. However, the electric and magnetic fields, Equations (46) and (47), respectively, are everywhere well defined. This light surface is the boundary where the electric and magnetic fields have the same strength.
In the region jξj < 1, the electromagnetic field is of the pure magnetic type; rotating observers described by the monad field (Equation (36)) detect only a dipolar-like magnetic field: In the region |ξ | >1, the electromagnetic field is of the pure electric type, and the rotating observers described by the monad field (Equation (41)) detect only an electric field given by the following: In the following sections, some illustrative examples of rotating point magnetic dipoles are presented. They show the rich structure of their magnetospheres and give an idea of how different they may be. It is interesting that the charge density shows regions of different polarities, depending on whether the angular velocity and the magnetic moment are parallel or antiparallel. These regions of definite sign charge offer a safe environment for antimatter to accumulate and in principle could be used to explain the antiprotons found in Earth's magnetosphere; see [26,27]. In addition, the electromagnetic field's rotation velocities are provided and the light surfaces are presented for each case. 6 Advances in High Energy Physics

Point Magnetic Dipoles with Rigid Rotation.
To consider a rotating point magnetic dipole in rigid rotation, we use the following: which plugged in Equations (37) and (42) leads to the following: Hence, only comoving observers I move rigidly, while comoving observers II, outside the light cylinder, move with a subluminal speed; see Figure 1. Thus, at most, the electromagnetic field reaches the speed of light on this light surface; see Figure 2(a). These results counter the idea that the electromagnetic fields do rotate with a superluminal speed outside of the light cylinder; see Fendt and Ouyed [7]. Notice that the speeds of the reference frames I and II were derived from the very electromagnetic field tensor (Equation (31)) assuming only a dipolar magnetic field nature (Equation (45)) and rigid rotation (Equation (54)).
The rigid dipole is not isolated as there are clouds of charge of different signs around it. In the inertial reference frame, this charge is as follows: which is well defined on the light cylinder and beyond it. If the angular velocity and the dipolar moment are parallel to each other, there are two lobular clouds of positive charge and one toroidal cloud of negative charge, while if they are antiparallel, the signs of the clouds become inverted. Figure 3(a) shows some revolution profiles of constant charge density for Ω 0 = 10 and k = 1, in suitable units. These clouds of charges are similar to that found in numerical studies of pulsar magnetospheres; see Kalapotharakos et al. [28] and Cerutti and Beloborodov [29]. It is also a special case of the charge density found by Hones and Bergeson [23]. Kalapotharakos et al. found numerically that the pulsar magnetosphere can be extended well beyond the light cylinder.

Point Magnetic Dipoles with Differential Rotation and a
Deformed Light Cylinder. As a toy model of a differentially rotating magnetic dipole, consider the following angular speed: where Nðr, ϑÞ is given by Equation (45), and a and Ω 0 are constants. Here, Ω I is a positive bounded function, and no other physical criteria will be taken into account. (37) and (42) leads to the following: v I = Ω 0 e −ak sin 2 ϑ/r r sin ϑ, v II = 1 Ω 0 r sin ϑ e ak sin 2 ϑ/r :

Substituting Equation (57) in Equations
ð58Þ v I < 1 is inside the light surface, which is a deformed cylinder, while v II < 1 is outside it. Figure 2(b) shows a revolution profile for the light surface corresponding to v I = 1 when the values a = 1, k = 1, and Ω 0 = 10, in some suitable units, are used.
In this case, the charge density measured in the reference frame at rest with the rotation axis is as follows: Figure 3(b) shows some revolution profiles for surface of constant density charge, for a = 1, k = 1, and Ω 0 = 10. When the dipolar moment and the angular velocity are antiparallel, the signs of the charge of these surfaces become interchanged.

Point Magnetic Dipoles with Differential Rotation and a
Light-Like Torus. Another possible choice, using the same requirement as before, is as follows: v I < 1 is outside the light surface, in this case a torus-like surface, while v II < 1 is inside it; the light surface is shown in Figure 2(c). The field is of pure magnetic type outside this torus and of pure electric type inside it.
In the inertial reference frame, at rest with the rotation axis, the charge density is as follows: Some revolution profiles of constant charge density are shown in Figure 3(c), using a = 1, k = 1, and Ω 0 = 10.

Nonpure Rotating Electromagnetic Fields
To find the reference frame comoving with the electromagnetic field, one looks for the frame in which the Poynting vector vanishes, i.e., that in which the electric and magnetic

Advances in High Energy Physics
fields are parallel to each other. To achieve this goal, we introduce an auxiliar 2-form; see Mitskievich [3]: where αðxÞ is an arbitrary function chosen conveniently. The original electromagnetic field tensor can be readily obtained from the following: Using the auxiliary 2-form, we construct the auxiliary invariants: which can be expressed in terms of the electromagnetic ones: I 1 = cos 2α ð ÞI 1 + sin 2α ð ÞI 2 ,  (56)); the cases of differential rotation are shown on the center and left, respectively, for charge densities (Equations (59) and (62)). The green curve corresponds to ρ = −1, blue is for ρ = −10, red for ρ = +10, and magenta for ρ = +1. These examples suggest that magnetospheres around astrophysical objects may be quite diverse. Figures were plotted using a = 1, k = 1, and Ω 0 = 10.

Advances in High Energy Physics
The second auxiliary invariant vanishes if we choose α function to satisfy the following: When the second invariant of a 2-form vanishes, one says that the 2-form is decomposable (meaning that it can be written as a simple bivector) [30,31]. However, the decomposition is not unique, as we can write the 2-form in three different ways: Here, Inserting Equation (63) in the last part of Equation (72), they become as follows: Substituting Equations (25) and (26), Due to the vanishing of Equation (68), and the fact that E ð3Þ = B ð3Þ = 0, both Equations (69) and (70) coincide: Since Now, the auxiliary 2-form F, Equation (75), is a simple bivector. Proceeding in the same way like for the pure electromagnetic field case, it is easy to find the corresponding monad fields.
In the case when |E ð1Þ | > | B ð2Þ | , we can choose the monad field as follows: Comparing Equation (77) with Equation (A.14), one finds the velocity field of the comoving observers with the electromagnetic field in the inertial reference frame, which is given by the following: Inserting Equation (75) into Equation (64), we have the following: where The expressions for the electric and magnetic fields in the frame (Equation (85)) are readily obtained comparing Equation (79) with Equation (3): Both covectors are parallel to each other; consequently, the Poynting vector vanishes in this reference frame.
In the opposite case, |B ð2Þ | > | E ð1Þ | , the auxiliary 2-form can be rewritten as follows: Now, inserting Equation (82) into Equation (64), the field tensor can be rewritten as follows: 10 Advances in High Energy Physics Comparing Equation (85) with Equation (A.14), one finds that the velocity field of the observers comoving with the electromagnetic field in the inertial reference frame is given by the following: The expressions for the electric and magnetic fields in the frame (Equation (77)) are readily obtained comparing Equation (83) with Equation (3): Clearly, both covectors are parallel to each other; consequently, there is no Poynting vector in this frame.

Superposition of a Coulombian Electric Field and the Field of a Point Magnetic Dipole
As an example of a nonpure field, the superposition of a Coulombian electric field and the magnetic field is studied here, considering the following functions: for which the electromagnetic field tensor is as follows: From Equation (89), and using Equation (4), the electric and magnetic field covectors in the inertial reference frame represented by the monad field τ = θ ð0Þ are as follows: The electromagnetic field invariants are as follows: From the second invariant, one can see that the electromagnetic field is of the nonpure type everywhere, except in the equatorial plane. The electromagnetic field is of the null type on the surface: Outside this peanut-like surface, the electromagnetic field is of the electric type, while in its interior is of magnetic type; there is no light surface, but a light ring (or a light circle) on the peanut waist, with radius r = |k/Q | , and only there, the field is of the pure null type; see Figure 4.
From a classical point of view, the electromagnetic field of an electron is like the one discussed here, using SI units: where M e~MBohr = 9:27 × 10 −24 J/T is the electron's magnetic moment, e = 1:60 × 10 −19 C is the electron's charge, and c = 2:99 × 10 8 m/s is the speed of light. In the following two subsections, we consider the velocity field of this superposition in the equatorial plane where the field is of the pure type, and outside this plane where the field is nonpure.

The Velocity Field on the Superposition's Pure Field
Region. On the equatorial plane, where the field is of the pure type, the field tensor (Equation (89)) is already a simple bivector. Using the ideas of Section 4, it can be rewritten as follows: for r < |k/Q | , the monad field is as follows: 11 Advances in High Energy Physics which allows to write Equation (99) as follows: Comparing Equation (97) with Equation (3), the magnetic field in the reference frame described by the monad field (Equation (96)) is as follows: The field tensor (Equation (89)) can also be rewritten as follows: In addition, for r > |k/Q | , we can use as a monad field: which allows to write Equation (99) as follows: Comparing Equation (101) with Equation (3), the electric field in the reference frame represented by Equation (100) is as follows: From Equation (96), inside the light surface, on the equatorial plane, the electromagnetic field is rotating rigidly in the inertial reference frame with velocity: while, outside of the light surface, also on the equatorial plane, it is rotating differentially with velocity.
7.2. The Velocity Field on the Superposition's Nonpure Field Region. The auxiliary covectors are as follows: Inserting the corresponding components of E and B into Equations (78) and (86), we find the velocity of the comoving observers described by the monad field τ B as measured by the inertial observers at rest with the rotation axis: Qr cos α − 2k cos ϑ sin α k sin ϑ cos α = k sin α sin ϑ Qr sin α + 2k cos α cos ϑ : The two equalities come from Equation (76). One can see that on the equatorial plane, the above expression reduces to Equation (103); hence, v B = v I . Hence, v A and v B are the velocities of the rotating observers, respectively, inside and outside the surface where the first invariant vanishes (Equation (92)).
The function αðr, ϑÞ is found by plugging Equation (91) in Equation (67): The speed field corresponding to Equation (106) is shown below on Figure 5; it can be seen that the speed of the electromagnetic field vanishes on the z-axis, growing to reach the speed of the light on a ring on the waist of the peanut-like surface. The white dotted line corresponds to the vanishing of the first electromagnetic invariant; inside it, the speed field for v B is plotted, while outside for v A . The field is everywhere of the  (92)). Inside this surface, the electromagnetic field is of the magnetic type, while outside it, of the electric type, the field is of the null type on this surface. The field is of the nonpure type everywhere, except on the equatorial plane, where it is pure (I 2 = 0). Only in a circle on the peanut's waist, the electromagnetic field reaches the speed of light and there both invariants vanish, in agreement with Equation (16). Hence, there is no light surface, but a light circle. 12 Advances in High Energy Physics nonpure type, except on the equatorial plane. Only on the equatorial plane, inside the surface I 1 = 0, the electromagnetic field rotates rigidly. Notice that the speed of light is reached only on the light ring (the intersection of the surface I 1 = 0 and the equatorial plane ϑ = π/2, where I 2 = 0); only there the electromagnetic field is of the pure null type; notice that this is in complete agreement with (Equation (16)).

Conclusions
In this work, a covariant 3 + 1 splitting formalism was used to study rotating axially symmetric electromagnetic fields in special relativity. Two sets of observers have been introduced: a family of observers at rest with the symmetry axis and a family of comoving observers with the electromagnetic field in rotation around the symmetry axis. These comoving observers do not detect any Poynting vector in their corresponding reference frame. Due to the presence of light surfaces (or light circles), i.e., the regions where the electromagnetic field rotates with the speed of light, it was not possible to fill the entire space with one single family of rotating observers, instead two subfamilies were required, one at each side of the surface. From a methodological point of view, to introduce the rotating observers, it was necessary to consider two cases: (1) For pure nonnull electromagnetic fields, I 1 ≠ 0 and I 2 = 0, two subsets of rotating observers fI, IIg have been introduced by looking for the reference frame where either the electric vector or the magnetic vector is absent (2) For nonpure electromagnetic fields, I 2 ≠ 0, other two subsets of rotating observers fA, Bg have been introduced, this time by looking for the reference frame where the electric and magnetic vectors are parallel to each other In each case, the exact expressions for these observers' velocities have been given with respect to the reference frame at rest with the symmetry axis.
In the case of pure fields, it was found that the electromagnetic field can rotate in a vast number of ways, from rigid rotation to different kinds of differential rotation. It also became apparent that for each rotational state of a magnetic field, there was a different associated electric field, although they were all similar up to a conformal factor; a different associated charge density and a different light surface. This light surface is in the regions where the electromagnetic field is of the pure null type. Several examples have been given; all of them are related to rotating magnetic dipoles, to show these latter features.
In the case of nonpure fields, the calculations have been more cumbersome. However, it has been also clear that they can also rotate in a vast number of ways. As an example, a superposition of a Coulombian electric field and dipolar magnetic field has been presented. The resulting field can display different kinds of rotation, from the rigid one to the differential one. The field's rotation speed reaches that of the light only in a circle on the equatorial plane.
As expected from Equation (16) The examples employed can be used not only for pedagogical purposes, but they are also of relevance in astrophysics and cosmology. In particular, these examples suggest that the magnetospheres of stars and other astrophysical objects may be quite diverse depending on the rotational state of the electromagnetic field. In the case of numerical studies of With respect to the z-axis, the left half is approaching while the right half is receding. The white dotted line corresponds to the vanishing of the first electromagnetic invariant; inside it, the speed field for v B is plotted, while outside for v A . The field is everywhere of the nonpure type, except on the equatorial plane. Only on the equatorial plane, inside the surface I 1 = 0, the electromagnetic field rotates rigidly. Notice that the speed of light is reached only on the light ring (the intersection of the surface I 1 = 0 and the equatorial plane ϑ = π/ 2), where I 2 = 0; only there the electromagnetic field is of the pure null type; notice that this is in complete agreement with Equation (16). 13 Advances in High Energy Physics neutron stars and other hard objects, it is necessary to take into account that corotating frames do not describe correctly the situation on the light surface, instead it should be better to consider a reference frame at rest with the rotation axis.
As suggested by other authors [10,11] and the examples presented here, the charge density in the magnetosphere of neutron stars is separated into regions of positive and negative charges. These regions offer a safe environment for antimatter to accumulate in time. Perhaps from the annihilation of this antimatter with cosmic rays, the charge density could be reconstructed. If this were the case, it would be possible to determine the rotation state of the electromagnetic field around a neutron star, in particular for those stars possessing a magnetic dipolar moment aligned with the rotation axis. As an application of this theory to the study of neutron stars and pulsars, it would be possible to get the function dM/dN, from Equation (49), a key ingredient of function ξ (Equation (32)). This function determines the form of the light surface, ξ = 1, Equation (35); the 4-velocity field of the electromagnetic field, Equations (36) and (41); and the structure of the electric field, Equation (46). These latter features are known so far from numerical studies of neutron stars with a rigidly corotating electromagnetic fields around them [21,28,29], but of course the corresponding features of the real neutron stars could be more diverse than these models suggest.