Relativistic Dynamics of Point Magnetic Moment

The covariant motion of a classical point particle with magnetic moment in the presence of (external) electromagnetic fields is revisited. We are interested in understanding Lorentz force extension involving point particle magnetic moment (Stern-Gerlach force) and how the spin precession dynamics is modified for consistency. We introduce spin as a classical particle property inherent to Poincare\'e symmetry of space-time. We propose a covariant formulation of the magnetic force based on a \lq magnetic\rq\ 4-potential and show how the point particle magnetic moment relates to the Amperian (current loop) and Gilbertian (magnetic monopole) description. We show that covariant spin precession lacks a unique form and discuss connection to $g-2$ anomaly. We consider variational action principle and find that a consistent extension of Lorentz force to include magnetic spin force is not straightforward. We look at non-covariant particle dynamics, and present a short introduction to dynamics of (neutral) particles hit by a laser pulse of arbitrary shape.


Introduction
The (relativistic) dynamics of particle magnetic moment µ, i.e. the proper time dynamics of spin s µ (τ ), has not been fully described before. Our interest in this topic originates in a multitude of current research topics: i) the ongoing effort to understand the magnetic moment anomaly of the muon [1,2]; ii) questions regarding how elementary magnetic dipoles (e.g. neutrons) interact with external fields [3,4]; iii) particle dynamics in ultra strong magnetic fields created in relativistic heavy ion collisions [5,6]; iv) magnetars, stellar objects with extreme O(10 11 ) T magnetic fields [7,8]; v) the exploration of particle dynamics in laser generated strong fields [9]; vi) neutron beam guidance and neutron storage rings [10]; and vii) the finding of unusual quantum spin dynamics when gyromagnetic ratio g = 2 [11,12]. The results we present will further improve the understanding of plasma physics in presence of inhomogeneous magnetic fields, and improve formulation of radiation reaction forces, topics not further discussed in this presentation.
In the context of the electromagnetic (EM) Maxwell-Lorentz theory we learn in the classroom that 1. The magnetic moment µ has an interaction energy with a magnetic field B The corresponding Stern-Gerlach force F SG has been written in two formats The name 'Amperian' relates to the loop current generating the force. The Gilbertian model invokes a magnetic dipole made of two magnetic monopoles. These two forces written here in rest frame of a particle are related [3,4]. We will show that a internal spin based magnetic dipole appears naturally; it does not need to be made of magnetic monopoles or current loops. We find that both force expressions in Eq. (2) are equivalent, this equivalence arises from covariant dynamics we develop and requires additional terms in particle rest frame complementing those shown in Eq. (2). 2. The torque T that a magnetic field B exercises on a magnetic dipole µ in a way that tends to align the dipole with the direction of a magnetic field B The magnetic moment is defined in general in terms of the product of Bohr magneton µ B with the gyromagnetic ratio g, |µ| ≡ gµ B . In Eq. (3) we used |s| = /2 for a spin-1/2 particle, a more general expression will be introduced in subsection 3.1.1.
We used the same coefficient µ to characterize both the Stern-Gerlach force Eq. (2) and spin precession force Eq. (3). However, there is no compelling argument to do so and we will generalize this hypothesis -it is well known that Dirac quantum dynamics of spin-1/2 particles predicts both the magnitude g = 2 and identity of magnetic moments entering Eq. (2) and Eq. (3).
While the conservation of electrical charge is rooted in gauge invariance symmetry, the magnitude of electrical charge has remained a riddle; the situation is similar for the the case of the magnetic moment µ: spin properties are rooted in the Poincaré symmetry of space-time, however, the strength of spin interaction with magnetic field, Eq. (1) and Eq. (3), is arbitrary but unique for each type of (classical) particle. Introducing the gyromagnetic ratio g we in fact create an additional conserved particle quality. This becomes clearer when we realize that the appearance of 'e' does not mean that particles we study need to be electrically charged.
First principle considerations of point particle relativistic dynamics experience some difficulties in generating Eqs. (2,3), as a rich literature on the subject shows -we will cite only work that is directly relevant to our approach; for further 70+ references see the recent numerical study of spin effects and radiation reaction in a strong electromagnetic field [9].
For what follows it is important to know that the spin precession Eq. (3) is a result of spatial rotational invariance which leads to angular and spin coupling, and thus spin dynamics can be found without a new dynamical principle has been argued e.g. by Van Dam and Ruijgrok [13] and Schwinger [14]. Similar physics content is seen in the work of Skagerstam and Stern [15,16], who considered the context of fiber bundle structure focusing on Thomas precession. Covariant generalization of the spin precession Eq. (3) is often attributed to the 1959 work by Bergmann-Michel-Telegdi [17]. However we are reminded [18,19,20] that this result was discovered already 33 years earlier by L.H. Thomas [23,24] at the time when the story of the electron gyromagnetic ratio g = 2 was unfolding. Following Jackson [18] we call the corresponding equation TBMT. J. Frenkel, who published [21,22] at the same time with L.H. Thomas explored covariant form of the Stern-Gerlach force, a task we complete in this work.
There have been numerous attempts to improve the understanding of how spin motion back-reacts into the Lorentz force generating the the Stern-Gerlach force. In the 1962 review Nyborg [25] summarized efforts to formulate the covariant theory of electromagnetic forces including particle intrinsic magnetic moment. In 1972 Itzykson and Voros [26] proposed a covariant variational action principle formulation introducing the inertia of spin I, seeking consistent variational principle but they found that no new dynamical insight resulted in this formulation.
Our study relates most to the work of Van Dam and Ruijgrok [13]. This work relies on an action principle and hence there are in the Lorentz force inconsistent terms that violate the constraint that the speed of light is constant, see e.g. their Eq. 3.11 and remarks: ' The last two terms are O(e 2 ) and will be omitted in what follows.' Other authors were proposing mass modifications to compensate these terms, a step which is equally unacceptable. For this reason our approach is intuitive, without insisting on 'in principle there is an action'. Once we have secured a consistent, unique covariant extension of the Lorentz force, we explore the natural variational principle action. We find it is not consistent and we identify the origin of the variational principle difficulties.
We develop the concept of the classical point particle spin vector in the following section 2. Our discussion relates to Casimir invariants rooted in space-time symmetry transformations. Using Poincaré group generators and Casimir eigenvalues we construct the particle momentum p µ and particle space-like spin pseudo vector s µ . In section 3 we present a consistent picture of the Stern-Gerlach force (subsection 3.1) and generalize TBMT precession equation (subsection 3.2) linear in both, the EM field and EM field derivatives. We connect the Amperian form of SG force (3.1.1) with the Gilbertian force (3.1.2). We discuss non-uniqueness of spin dynamics (3.2.3) in consideration of impact on muon g − 2 experiment. We show in section 4 that the natural choice of action for the considered dynamical system does not lead to a consistent set of equations; in this finding we align with all prior studies of Stern-Gerlach extension of the Lorentz force.
In the final part of this work section 5 we show some of the physical content of the theoretical framework. In subsection 5.1 we present a more detailed discussion of dynamical equations for the case of particle in motion with a given β = v/c and E, B fields in laboratory. In subsection 5.2 we study solution of the dynamical equations for the case of an EM light wave pulse hitting a neutral particle. We have obtained exact solutions of this problem, detail will follow under seperate cover [27]. The concluding section 6 is a brief summary of our findings.

Notation
For most of notation, see Ref. [28]. Here we note that we use SI unit system and the metric: We further recognize the totally antisymmetric covariant pseudo tensor ǫ: where 'perm' is the signature of the permutation. It is important to remember when transiting to non-covariant notation in laboratory frame of reference that the analog contravariant pseudo-tensor due to odd number of spacelike dimensions is negative for even permutations and positive for odd permutations. The Appendix B of Ref. [29] presents an introduction to ǫ. We will introduce an elementary magnetic dipole charge d -the limitations of the alphabet force us to adopt the letter d otherwise used to describe the electric dipole to be the elementary magnetic dipole charge. The magnetic dipole charge of a particle we call d converts the spin vector s to magnetic dipole vector µ, The factor c is needed in SI units since in the EM-tensor F µν has as elements E/c and B. It seems natural to introduce also s µ d = µ µ -an object µ µ can confuse and we stick to the product s µ d, however we always replace sd → µ/c. Note that we place d to right of pertinent quantities to avoid confusion such as dx.
We cannot avoid appearance in the same equation of both magnetic moment µ and vacuum permeability µ 0 .

Spin Vector
A classical intrinsic covariant spin has not been clearly defined or even identified in prior work. In some work addressing covariant dynamics of particles with intrinsic spin and magnetic moment particle spin is by implication solely a quantum phenomenon. Therefore we describe the precise origin of classical spin conceptually and introduce it in explicit terms in the following.
Considering the Poincare group of space-time symmetry transformations [30,31], it has been established that elementary particles have to be in a representation that is characterized by eigenvalues of two Casimir operators (a 'bar' marks operators) All physical point 'particles' have fixed eigenvalues of C 1 , C 2 . The quantities (with a bar)p µ andw α are differential operators constructed from generators of the symmetry transformations of space-time; that is 10 generators of the Poincaré group of symmetry transformations of 4-spacetime:p µ for translations, J for rotations and K for boosts. Once we construct suitable operator valued quantities we will transition to the physics of 'c-number' valued (without bar) variables as used in classical dynamics where all quantities will be normal numbers and rely on the eigenvalues of Casimir operators C 1 , C 2 for each type of particle. In Eq. (5) the first of the space-time operators based on generators of the four space-time translations p µ guarantees that a point particle has a conserved inertial mass m (with a value specific for any particle type). The second Casimir operator C 2 is obtained from the square of the Pauli-Lubański 4-vector Here M µν is the antisymmetric tensor (operator) created from three Lorentz-boost generators K and three space rotation generators J such that These relations help us see that The generators J , K of space-time transformations are recognized by their commutation relations. They are used in a well known way to construct representations of the Lorentz group.
In terms of the generator tensor M νµ the covariant definition of the particle spin (operator) vector is According to Eq. (8), spins µ is a pseudo vector, as required for angular dynamics. The dimension ofs µ is the same as the dimension of the generator of space rotations J. We further find thats µ is orthogonal to the 4-velocity (operator)ū µ by virtue of the antisymmetry of M ⋆ evident in the definition Eq. (6). The definition of the particle spin (operator) is unique: no other space-like (space-like given the orthogonalitys ·ū = 0) pseudo vector associated with the Poincaré group describing space-time symmetry transformations can be constructed.
We now transition to c-numbered quantities (dropping the bar): an observer '(0)' co-moving with a particle measures the 4-momentum and 4-spin s µ where according to conventionẑ-axis of the coordinate system points in direction of the intrinsic spin vector s. In the particle rest frame we see that consistent with operator equation Eq. (9); more generally, any space-like vector is normal to the time-like 4-velocity vector. For the magnitude of the spin vector we obtain We keep in mind that s 2 must always be a constant of motion in any frame of reference. Its value s · s = −s 2 is always negative, appropriate for a space-like vector. Similarly must be a constant of motion in any frame of reference and the value p 2 is positive, appropriate for a time-like vector.
As long as forces are small in the sense discussed in Ref. [28] we can act as if rules of relativity apply to both inertial and (weakly) accelerated frames of reference. This allows us to explore the action of forces on particles in their rest frame where Eq. (10) defines the state of a particle. By writing the force laws in covariant fashion we can solve for dynamical evolution of p µ (τ ), s µ (τ ) classical numbered variables.

Magnetic dipole potential and Amperian force
We have gone to great lengths in section 2 to argue for the existence of particle intrinsic spin. For all massive particles this implies the existence of a particle intrinsic magnetic dipole moment, without need for magnetic monopoles to exist or current loops. Spin naturally arises in the context of symmetries of Minkowski space-time, it is not a quantum property.
In view of above it is appropriate to study classical dynamics of particles that have both, an elementary electric charge e, and an elementary magnetic dipole charge d. The covariant dynamics beyond the Lorentz force needs to incorporate the Stern-Gerlach force. Thus the extension has to contain the elementary magnetic moment of a particle contributing to this force. To achieve a suitable generalization we introduce the magnetic potential We use dual pseudo tensor since s µ is a pseudo vector; the product in Eq. (14) results in a polar 4-vector B µ . We note that the magnetic dipole potential B µ by construction in terms of antisymmetric field pseudo tensor F ⋆ µν satisfies The additional potential energy of a particle at rest placed in this magnetic dipole potential is (16) This shows Eq. (14) describes the energy content seen in Eq. (1); all factors are appropriate.
The explicit format of this new force is obtained when we use Eq. (14) to define a new antisymmetric tensor (17) Equation (17) allows us to add to the Lorentz force In the G-tensor we note appearance in the force of the derivative of EM fields, required if we are to see the Amperian model variant of the Stern-Gerlach force Eq. (2) as a part of generalized Lorentz force. The Amperian-Stern-Gerlach (ASG) force 4-vector is obtained multiplying with u ν d the G-tensor Eq. (17). Thus the total 4-force a particle of charge e and magnetic dipole charge d experiences is In the particle rest frame we have We can use Eq. (20) to read-off from Eq. (18) the particle rest frame force to be where two contributions ∂(µ · B)/∂t to F 0 cancel. Each of the three terms originates in one of the covariant terms in the sequence shown. The result is what one calls Amperian model originating in dipoles created by current loops. This is, however, not the last word in regard to the form of the force.

Gilbertian model Stern-Gerlach force
We restate the Stern-Gerlach-Lorentz force Eq. (18), showing the derivative terms explicitly, Multiplying with s µ the last term vanishes due to antisymmetry of F ⋆ and we obtain This equation suggests that we explore as the generalized Lorentz force replacing the usual field tensor eF by F in a somewhat simpler way compared to the original H µν Eq. (18) modification.
We demonstrate now that the field modification seen in Eq. (24) leads to a different and fully equivalent format of the force. We replace in the first term in Eq. (22) F → F and add the extra term from Eq. (24) to the two reminder terms. Changing the index naming these we can write symmetrically The tensor appearing in the parentheses in the 2nd line of Eq. (25) is antisymmetric under any of the three exchanges of the indices. It is therefore proportional to the totally antisymmetric tensor ǫ αµβγ which must be contracted with some 4-vector V γ containing a gradient of the EM dual field tensor, there are two such available 4-vectors ∂ κ F ⋆ κγ which vanishes by virtue of Maxwell equations, and Thus we introduce the Gilbertian form of the 4-force Note that in our formulation the Amperian and the Gilbertian 4-forces are identical they are just written differently.
In the rest frame of a particle, see Eq. (20) the Gilbertian force Eq. (27) is It is interesting to see the mechanism by which the two formats of the forces equal to each other in the particle rest frame. With we form the force difference between Eq. (21) and Eq. (28) The terms in parenthesis cancel according to Maxwell equation confirming that both the Amperian and the Gilbertian forces are equal taking as an example the instantaneous rest frame. From now one we will use Gilbertian form of the force and in later examples we will focus on particle motion in vacuum, j µ = 0.
In this discussion of forces we kept the electrical charge e and the elementary magnetic moment 'charge' d Eq. (4) as independent qualities of a point particle. As noted in the introduction it is common to set |µ| ≡ gµ B , see below Eq. (3). Hence we can have both, charged particles without magnetic moment, or neutral particles with magnetic moment, aside from particles that have both charge and magnetic moment. For particles with both charge and magnetic moment we can write, using Gilbertian format of force where a = (g − 2)/2 is the gyromagnetic ratio anomaly. The Compton wavelengthλ = /mc defines the scale at which the spatial field inhomogeneity is relevant; note that inhomogeneities of the field are boosted in size for a particle in motion, a situation which will become more explicit in section 5.1.3.

Conventional TBMT
For particles with m = 0 differentiating Eq. (11) with respect to proper time we finḋ Here F µν is the usual EM field tensor. Equation (32) has the general TBMT solution where we used the notation u · F · s ≡ u µ F µν s ν . In Eq. (33) a is an arbitrary constant considering that the additional term multiplied with u µ vanishes. On the other hand we can read off the magnetic moment entering Eq. (3): the last term is higher order in 1/c 2 . Hence in the rest frame of the particle we see that 2(1 + a) = g i.e. Eq. (33) reproduces Eq. (3) with the magnetic moment coefficient when a = a. Therefore, as introduced, a = a is the g = 2 anomaly. However, in Eq. (33) we could for example use a = (g 2 − 4)/8 = a + a 2 /2, which classical limit of quantum dynamics in certain specific conditions implies [12]. In this case a → a up to higher order corrections. This means that measurement of a as performed in experiments [1,2] depends on derivation of the relation of a with a obtained from quantum theory. These remarks apply even before we study gradient in field corrections.

Gradient corrections to TBMT
The arguments by Schwinger, see Eqs. (31,32,33), are ideally positioned to obtain in a consistent way generalization of the TBMT equations including the gradient of fields terms required for consistency. We use Eq. (24) in Eq. (33) to obtain The dominant gradient of field correction arises for an elementary particle from the 2nd term in the first line in Eq. (34), considering the coefficient of the second line a = α 2 /2π + . . . = 1.2 × 10 −3 . One should remember that given the precision of the measurement [1,2] of a which is driven by the first term in the second line in Eq. (34) we cannot in general neglect the new 2nd term in first line in Eq. (34) even if the characteristic length defining the gradient magnitude is the Compton wavelengthλ, see Eq. (30).

Non-uniqueness of gradient corrections to TBMT
It is not self evident that the form Eq. (34) is unique. To see that a family of possible extensions TBMT arises we recall the tensor Eq. (18) H µν made of the two potentials A µ and B µ . We now consider the spin dynamics in terms of the two field tensors, F and G replacing the usual EMtensor F µν in the Schwinger solution, Eq. (33). In other words, we explore the dynamics according to Two different constants a and b are introduced now since the two terms shown involving F and G tensors could be included in Schwinger solution independently with different constants. Intuition demands that a = b. However, aside from algebraic simplicity we do not find any compelling argument for this assumption.
We return now to the definition of the G tensor Eq. (17) to obtain The first term in the first line vanishes by antisymmetry of F ⋆ tensor. We also have Using Eq. (36) and Eq. (37) we can combine in Eq. (35) the first two terms in both lines, and the last terms in both lines to obtain This equation agrees with Eq. (34) only when a = b. However, this requirement is neither mathematically nor physically necessary. For example using Eq. (26) we easily check s ·u + u ·ṡ = 0 without any assumptions about a, b. As Eq. (35) shows the physical difference between factors a and b is related to the nature of the interaction: the 'magnetic' tensor G is related to b only. Thus for a neutral particle e → 0 we see in Eq. (38) that the torque depends only on b. Conversely, when the effect of magnetic potential is negligible Eq. (38) becomes the textbook spin dynamics that depends on a alone.
To make further contact with textbook physics we note that the coefficient of the first term in Eq. should reproduce in leading order the torque coefficient in Eq. (3) as is expected from study of quantum correspondence. However, quantum correspondence could mean a = a + a 2 /2, which follows comparing exact solutions of the Dirac equation with spin precession for the case we explored [12] and which is not exactly the motion of a muon in storage ring. However, this means that in order to compare the measurement of magnetic moment of the muon carried out on macroscopic scale [1,2] with quantum computations requires a further step, the establishment of quantum correspondence at the level of precision at which the anomaly is measured.

Search for variational principle action
At the beginning of earlier discussions of a covariant extension to the Lorentz force describing the Stern-Gerlach force was always a well invented covariant action. However, the Lorentz force itself is not a consistent complement of the Maxwell equations. The existence of radiation means that an accelerated particle experiences radiation friction. The radiation-reaction force has not been incorporated into a variational principle [28,32]. Thus we should not expect that the Stern-Gerlach force must originate in a simple action. We seek a path x µ (τ ) in space-time that a particle will take considering an action that is a functional of the 4-velocity u µ (τ ) = dx µ /dτ and spin s µ (τ ). Variational principle requires an action I(u, x; s). When I respects space-time symmetries the magnitudes of particle mass and spin are preserved in the presence of electromagnetic (EM) fields. We also need to assure that u 2 = c 2 which constrains the form of force and thus I that is allowed. Moreover, we want to preserve gauge invariance of the resultant dynamics.
The component in the action that produced the LHS (inertia part) of the Lorentz force remains in discussion. To generate the Lorentz force one choice of action is We note that reparametrization of τ → kτ considering u = dx/dτ has no effect on value of I Lz .
Variation with respect to path lead to where the RHS produces upon differentiation of eA µ (x(τ )) the usual Lorentz force Multiplying Eq. (41) with mcu µ / √ u 2 we establish by antisymmetry of the tensor F µν Eq. (42) that also the product with the LHS in Eq. (41) vanishes. This means that (mcu ν / √ u 2 ) 2 = m 2 c 2 ≡ p 2 = Const. Henceforth There is a problem when we supplement in Eq. (40) the usual action I Lz by a term I m based on our prior consideration of A µ → A µ + B µ , see subsection 3.1.1. The problem one encounters is that the quantity B µ contains additional dependence on s µ (τ ) which adds another term to the force. Let us look at the situation explicitly Here the dependence on s µ (τ ) is akin to a parameter dependence; some additional consideration defines the behavior, in our case this is the TBMT equations.
Varying with respect to the path the modified action Eq. (44) we find the modified covariant force with two new contributions We applied here with A → B the result seen in Eq. (41), and the additional term L µ S2 follows by remembering to take proper time derivative of s µ . The first term Eq. (46) is as we identified previously in Eq. (18). We note that another additional term arises if and when an additional power of √ u 2 to accompany u · B as was done in [13]. An unsolved problem is created by the torque-like term, Eq. (47).
If we replace in our thoughts ds ν /dτ in Eq. (47) by the TBMT equation Eq. (33) or as would be more appropriate by its extended version Eq. (35), we see that the force L µ S2 would be quadratic in the fields containing also field derivatives. However, by assumption we modified the action limiting the new term in Eq. (44) to be linear in the fields and derivatives. Finding non linear terms we learn that this assumption was not justified. However, if we add the quadratic in fields term to the action we find following the chain of arguments just presented that a cubic term is also required and so on; with derivatives of fields appearing along.
We have searched for some time for a form that avoids this circular conundrum, but akin to previous authors we did not find one. Clearly a 'more' first principle approach would be needed to create a consistent variational principle based equation system. On the other hand we have presented before a formulation of spin dynamics which does not require a variational principle in the study the particle dynamics: as is we have obtained a dynamical equation system empirically. Our failing in the search for an underlying action is not critical. A precedent situation comes to mind here: the radiation emitted by accelerated charges introduces a 'radiation friction'which must be studied [28,32] without an available action, based on empirical knowledge about the energy loss arising for accelerated charges. In most physical cases we create a particle guiding field which is at rest in the laboratory. Particle motion occurs with respect to this prescribed field and thus in nearly all situations it is practical to study particle position z µ (τ ) in the laboratory frame of reference. Employing the Lorentzcoordinate transformations from the particle rest frame to the laboratory frame we obtain where as usual γ = 1/ 1 − β 2 and one often sees the spin written with γ 2 /(γ + 1) = (γ − 1)/β 2 . One easily checks that Eq. (48) and Eq. (49) also satisfy Eq. (11): u µ s µ = 0. A classic result of TBMT reported in textbooks is that the longitudinal polarizationβ · s for g ≃ 2 and β → 1 is a constant of motion. This shows that for a relativistic particle the magnitude of both time-like and space-like components of the spin 4-vector Eq. (49) can be arbitrarily large, even if the magnitude of the 4vector is bounded s µ s µ = −s 2 . This behavior parallels the behavior of 4-velocity u µ u µ = c 2 .
We remind that to obtain in the laboratory frame the usual Lorentz force we use the 4-velocity with respect to the Laboratory frame Eq. (48), with laboratory defined tensor F , i.e. with laboratory given E, B EM-fields Sometimes it is of advantage to transform Eq. (50) to the particle rest frame. Such a transformation L with Lu| rest = u L when used on the left hand side in Eq. (50) produces proper time differentiation of the transformation operator, see also [33]. Such transformation into a co-rotating frame of reference originates the Thomas precession term in particle rest frame for the torque equation. This term is naturally present in covariant formulation when we work in the laboratory reference frame. For the full force Eq. (26) we thus have We see that in laboratory frame of reference a covariant gradient of the fields is prescribed, i.e. that some apparatus prescribes the magnitude which allows for a moving particle with u µ | L Eq. (48) and s µ | L Eq. (49) to experience the Stern-Gerlach force F µ We have gone to extraordinary length in arguing Eq. (54) to make sure that the forthcoming finding of the Lorentz boost of field inhomogeneity is not questioned.

Magnetic potential in the laboratory frame
We evaluate in the laboratory frame the form of Eq. (14). The computation is particularly simple once we first recall the laboratory format of the Lorentz force F µ The magnetic part of the action will be evaluated (see second line below) in analogy to above. We now consider where we used in 2nd line i) F ⋆ µν follows from the usual F µν upon exchange of E/c ↔ B and ii) flip β → −β to account for contravariant and not covariant 4-velocity.
For the magnetic action potential energy of a particle in lab frame we obtain Equation (57) extends the rest frame β = 0 Eq. (16) and represents covariant generalization of Eq. (1). In ultrarelativistic limit all terms in Eq. (57) have the same magnitude.

Field to particle energy transfer
We now consider the energy gain by a particle per unit of laboratory time, that is we study the 0th component of Eq. (26) A further simplification is achieved considering where the total derivative with respect to time accounts for both, the change in time of the laboratory given field B, and the change due to change of position in the field by the moving particle. We thus find two parts where the 2nd line is of particular interest as it is proportional to γ. Focusing our attention on this last term: we can use β = cp/E and γβ = p/mc. Upon multiplication with E and remembering that c 2 pdp = EdE we obtain which in qualitative terms implies an exponential response of particle momentum as it crosses a magnetic field |p| ≃ mc e ±(| B|− B0))|µ|/mc 2 .
However, even a magnetar magnetic field of up 10 11 T will not suffice to impact electron momentum decisively in view of the smallness of the electron magnetic moment 5.810 −11 MeV/T. However, in ultrarelativistic heavy ion collisions at LHC a 10,000 stronger very non-homogeneous B-fields arise.

Properties of equations
The dynamical equations developed here have a considerably more complex form compared to the Lorentz force and TBMT spin precession in constant fields [33]. We need field gradients in the Stern-Gerlach force, and in the related correction in the TBMT equations. Since the new physics appears only in the presence of a particle magnetic moment, we simplify by considering neutral particles. We now show that the external field described by a light wave (pulse) lends itself to an analytical solution effort. This context could be of practical relevance in the study of laser interaction with magnetic atoms, molecules, the neutron and maybe neutrinos. For e = 0 our equations Eq. (26) and Eq. (38) reaḋ The external light wave field is a pulse with The derivative of the dual EM tensor for linear fixed in space pulse polarization ε µ is prime '′' indicates derivative with respect to the phase ξ.
Notice that if we contract Eq. (67) with k µ or ε µ we get zero because Levi-Civita tensor ǫ µναβ is totally antisymmetric. Therefore contracting Eq. (64) with either k µ or ε µ we find We further note that the argument of the light pulse Eq. (66) satisfies where we used Eq. (68). Thus we conclude that the particle follows the pulse such that The two conservation laws Eq. (68) and Eq. (69) along with Eq. (70) make the light pulse an interesting example amenable to an analytical solution.
We now evaluate several invariants in the laboratory frame seeking understanding of their relevance. A particle moving in the laboratory frame in consideration of Eq. (48) experiences in its rest frame a plane wave with the Doppler shifted frequency which is unbounded as it grows with particle laboratory Lorentz-γ 0 . However, k·s, the projection of spin onto plane wave 4-momentum k µ is bounded. To see this we recall the constraint Eq. (11) which in the laboratory frame reads We thus obtain where we used Eq. (73) in last equality. Since β and n = k/|k| are unit-magnitude vectors we find The magnitude of the spin vector in the lab frame is constrained by Eq. (12) where we again used Eq. (73). Combining Eq. (75) and Eq. (76) we see that except when particle is moving exactly in direction of S L (sin 2 θ = 0), the magnitude of (k · s(τ )) 2 is bounded.

Invariant acceleration and spin precession
Even without knowing the explicit form for u µ (τ ), s µ (τ ) we were able to obtain [27] the invariant acceleratioṅ This result follows using the usual trick of taking a further (proper) time derivative of Eq. (64) (multiplied by a suitable factor) and on RHS eliminatingu by using Eq. (64).
Multiplying the result with u µ and eliminating u ·ü using the 2nd differential of u 2 = c 2 produces Eq. (77). We see in Eq. (77) that the magnitude of the 4-force created by a light pulse and acting on an ultrarelativistic particle is dependent on square of the product of the 2nd derivative of pulse function with respect to ξ, f ′′ (ξ), with the Doppler shifted frequency Eq. (72). The value Eq. (77) is negative since acceleration is a space-like vector.
As we discussed below Eq. (76) the spin precession factor k·s seen in Eq. (77) is bounded. We were able to obtain a soluble formulation of the spin precession dynamics described by the dimensionless variable y = k · s(τ ) bd mc C 1 (78) which satisfies the differential equation (79) obtained performing suitable manipulations of dynamical equations prior to solving for u µ (τ ), s µ (τ ). We are seeking bounded periodic solutions of nonlinear Eq. (79) no matter how large the constant C 1 determined by the initial conditions C 2 contains the initial particle Lorentz-γ factor. One can see several possible solutions of interest of Eq. (79); for example y = sin(φ(s)) satisfies all constraints. It leads to the pendulum type differential equation and we recognize that high intensity light pulses can flip particle spin. However, there are other relevant solutions, e.g. y ∝ 1/ cosh z. Upon solution of Eq. (79) k·s(τ ) is known, given Eq. (71) we also know the dependence of Eq. (67) on proper time τ . Hence Eq. (64) can be solved for u µ and Eq. (65) can be solved for s µ resulting in an analytical solution of the dynamics of a neutral magnetic dipole moment in the field of a light pulse of arbitrary shape. The full description of the dynamics exceeds in length this presentation and will follow [27].

Conclusions
The Stern-Gerlach covariant extension of the Lorentz force has seen considerable interest as there are many immediate applications listed in first paragraph. Here we have: 1) introduced in Eq. (10) the covariant classical 4-spin vector s µ in a way expected in the context of Poincare symmetry of space-time; 2) presented a unique linear in fields form of the covariant magnetic moment potential, Eq. (14), which leads to a natural generalization of the Lorentz force; 3) shown that the resultant Amperian, Eq. (19), and Gilbertian, Eq. (26), forms of the magnetic moment force are equivalent; 4) extended the TBMT torque dynamics, Eq. (35), making these consistent with the modifications of the Lorentz force; 5) demonstrated the need to connect the magnetic moment magnitude entering the Stern-Gerlach force with the one seen in the context of torque dynamics, subsection 3.2.3; 6) shown that variational principle based dynamics has systemic failings when both position and spin are addressed within present day conceptual framework, see section 4; 7) reduced the covariant dynamical equations to laboratory frame of reference uncovering important features governing the coupled dynamics, see section 5.1; 8) obtained work done by variations of magnetic field in space-time on a particle, Eq. (61); 9) shown salient features of solutions of neutral particles with non-zero magnetic moment hit by a laser pulse, see section 5.2;