Frenkel electron on an arbitrary electromagnetic background and magnetic Zitterbewegung Alexei A. Deriglazov∗ and Andrey M. Pupasov-Maksimov† Depto. de Matem´ atica, ICE, Universidade Federal de Juiz de Fora, MG, Brasil (Dated: January 31, 2014)

arXiv:1401.7641v2 [hep-th] 30 Jan 2014

We present Lagrangian which implies both necessary constraints and dynamical equations for position and spin of relativistic spin one-half particle. The model is consistent for any value of magnetic moment µ and for arbitrary electromagnetic background. Our equations coincide with those of Frenkel in the approximation in which the latter have been obtained by Frenkel. Transition from approximate to exact equations yields two structural modifications of the theory. First, Frenkel condition on spin-tensor turns into the Pirani condition. Second, canonical momentum is no more proportional to velocity. Due to this, even when µ = 1 (Frenkel case), the complete and approximate equations predict different behavior of spinning particle. The difference between momentum and velocity means extra contribution into spin-orbit interaction. To estimate the contribution, we found exact solution to complete equations for the case of uniform magnetic field. While BMT electron moves around the circle, our particle experiences magnetic Zitterbewegung, that is oscillates in the direction of magnetic field with amplitude of order of Compton wavelength for the fast particle.

I.

INTRODUCTION

Consistent and complete description of spin effects of the relativistic electron is achieved in QED on the base of Dirac equation. However, starting from the pioneer works [1–3] and up to present date, interpretation of final results in some cases is under permanent and controversial debates in various theoretical and experimental set-ups [4–13]. Understanding of the spin precession in the case of arbitrary magnetic moment in an external electromagnetic field is important in the development of experimental technics for measurements of anomalous magnetic moment [14, 15]. In accelerator physics [16] it is important to control resonances leading to depolarization of a beam. In the case of vertex electrons carrying arbitrary angular momentum, semiclassical description can also be useful [17]. So the relationship among classical and quantum descriptions remains an important step of analysis, providing the interpretation of results of QFT computations in usual terms: particles and their interactions. Hence an actual task is to develop, in a systematic form, the classical model of an electron [18–28, 30] which would be as close as possible to the Dirac equation. May be the best candidates for classical equations of relativistic electron are those of Frenkel [18, 19] and Bargmann, Michel and Telegdi (BMT) [20]. In this work we continue detailed analysis of these equations started in [31, 32], and construct Lagrangian which gives generalization of these equations on the case of arbitrary electromagnetic background. Even for non interacting theory, search for Lagrangian which yields the right number of degrees of freedom and the right non relativistic limit represents rather non trivial problem [19, 22, 24–30]. In [32] we have solved this problem, considering spin as a composed quantity (inner angular momentum) constructed from non-Grassmann vector-like variable and its conjugated momentum [27, 28, 30]. The desired properties are guaranteed by six constraints. So, in Hamiltonian formulation the model represents rather nontrivial example of a constrained system. Phase space of the model turns out to be curved manifold equipped, in a natural way, with the structure of fiber bundle. Detailed analysis of the underlying geometry has been presented in [31]. This allowed us to develop the proper quantization scheme. In [32] we have performed both canonical (in physical-time parametrization) and manifestly covariant (in arbitrary parametrization) quantization of free model, and established the relation with one-particle sector of Dirac equation as well as with quantum theory of two-component Klein-Gordon equation developed by Feynmann and Gell-Mann [33]. It has been demonstrated that various known in the literature non-covariant, covariant and manifestly-covariant operators of position and spin acquire clear meaning and interpretation in the Lagrangian model of Frenkel electron. In particular, we have found the manifestly covariant form of position and spin operators in the space of positive-energy Dirac spinors.

∗ Electronic † Electronic

address: [email protected] address: [email protected]

2 In the Hamiltonian formulation appeared two second-class constraints which, at the end, supply the Frenkel condition on spin-tensor. They depend on both position and spin-sector variables. This leads to new properties as compared with non relativistic spin [31]. The constraints must be taken into account by transition from Poisson to Dirac bracket, this leads to nonvanishing classical brackets for the position variables. In the result, the position space is endowed, in a natural way, with noncommutative structure which originates from accounting of spin degrees of freedom. Our model represents an example of situation when physically interesting noncommutative relativistic particle emerges in a natural way. For the case, the ”parameter of non-commutativity” is proportional to spin-tensor. As a consequence, operators corresponding to position of the electron are non-commutative (they can be identified [32] with Pryce (d) operators). This implies that effects of non-commutativity could be presented at the Compton wave length, in contrast to conventional expectations of non-commutativity at Planck length. There are a lot of candidates for spin and position operators of the relativistic electron [1–3, 21, 34]. Different position observables coincide when we consider standard quasi-classical limit. So, in absence of a systematically constructed classical model of an electron it is difficult to understand the difference between these operators. Our approach allows us to do this, after realizing all them at the classical level. Besides, all the candidates obey the same equations in free theory, so the question of which of them are the true position and spin is a matter of convention. The situation changes in interacting theory, considered below, where we can distinguish the variables according their classical dynamics in an external field. In the present work we return to the classical mechanics and study an interacting theory. In section II we show that our Lagrangian admits interaction with an arbitrary electromagnetic background1 . The model contains two coupling constants - charge e and the interaction constant µ of basic spin variables with Fµν . Provisionally, we call this magnetic moment. The theory is consistent for arbitrary values of µ. For the position variable we have the minimal interaction term ec Aµ x˙ µ . As for spin, when the particle has non-vanishing magnetic moment, this interacts with electromagnetic field in a highly nonlinear way. This turns out to be necessary for preservation of the number and algebraic structure of constraints in the passage from free to interacting theory. In section III we present and analyze equations of motion in Hamiltonian and Lagrangian formulations. We show that they follow from simple and expected Hamiltonian (33), when we deal with the Dirac bracket. We compare our equations with those of Frenkel [18, 19]. Frekel considered the case µ = 1, and found his equations in the quadratic approximation on spin-tensor. We show that our exact equations coincide with those of Frenkel in these limits. Hence our Lagrangian gives complete Frenkel equations for arbitrary field and magnetic moment. Frenkel tensor can be used to construct BMT-type four-vector. We write the corresponding equations of motion. While Hamiltonian equations can be rewritten in closed form in terms of BMT vector, see Eq. (60)-(62), we do not achieved this for Lagrangian equations in our theory. In the Lagrangian form, the equation for BMT vector contains Frenkel tensor, see Eq. (63). It seems that Frenkel spin in our theory represents more fundamental object as compared with BMT spin. In section IV we find exact solution to our equations in the uniform magnetic field. With respect to BMT equations our model takes into account two physical effects. First, magnetic moment interacted with a magnetic field results in additional mass of electron. Second, in the case of anomalous magnetic moment the velocity and the momentum are not collinear, this modifies the Lorentz force. Our model naturally incorporates both these effects and leads to small corrections of the trajectory and spin precession. While equations of motion have a rather complicated structure, in the case of uniform magnetic field there are a lot of symmetries and hence integral of motions providing complete analytical solution.

II.

LAGRANGIAN AND HAMILTONIAN OF INTERACTING THEORY

To start with, we shortly describe the structure of free theory [32, 32]. Configuration space of the model consist of the position xµ (τ ) = (ct, x) and the vector-like variable of spin ω µ (τ ) = (ω 0 , ω), ω = (ω1 , ω2 , ω3 ) in an arbitrary parametrization τ . pµ and π µ are conjugate momenta for xµ and ω µ . The variables in the free theory are subject to the constraints T1 = p2 + (mc)2 = 0 ,

(1) 2

T3 = π 2 − a3 = 0 ,

T4 = ω 2 − a4 = 0 ,

T5 = ωπ = 0 ,

T6 = pω = 0 ,

1

Interaction with an arbitrary curved background is presented in [35].

where

T7 = pπ = 0 ,

a3 =

3~ , 4a4

(2) (3)

3 where ωπ = −ω 0 π 0 + ωπ and so on. As the Hamiltonian action functional, we simply take LH = px˙ + π ω˙ − H, with Hamiltonian H = gi Ti in the form of linear combination the constraints Ti multiplied by auxiliary variables gi , i = 1, 3, 4, 5, 6, 7. The constraint T3 belong to first-class and is related with local spin-plane symmetry presented in the theory [31]. The basic spin-sector variables change under the symmetry, so they do not represent an observables quantities. As the observable quantity we take the Frenkel spin-tensor J µν J µν = 2(ω µ π ν − ω ν π µ ).

(4)

The constraints (2) and (3) imply the following restrictions (in the free theory the conjugated momentum is proportional to four-velocity, pµ ∼ uµ ) J µν pν = 0,

J 2 = 6~2 .

(5)

Spacial components of the Frenkel tensor can be used to construct the quantity Si =

1 ijk  Jjk , 4

(6)

which we identify with non relativistic spin of Pauli theory. In the interacting theory pµ turns into canonical momentum P µ , which for non uniform fields or/and µ 6= 1 does not proportional to four-velocity. Hence in this case the Frenkel condition J µν uν = 0 turns into the Pirani condition [36–38] J µν Pν = 0. Frenkel tensor can be used to construct four-vector 1 sµ (τ ) ≡ p µναβ pν Jαβ , 2 4 −p

then sµ pµ = 0,

s2 =

3~2 . 4

(7)

In our theory, even in the case of interaction, the condition sµ Pµ = 0 implies sµ uµ = 0, see Eq. (58) below. So we can identify sµ with BMT vector [20]. In the rest frame, spacial components si of BMT vector coincide with S i . In arbitrary frame, they are related as follows:   p0 pi pj Si = p (8) δij − 0 2 sj . (p ) −p2 In the free theory the equation (7) can be inverted, J µν = − √ 2

−p2

µναβ pα sβ , so the two quantities are mathematically

equivalent. In the interacting theory, we have p → P = P(u, J, F ), see Eqs. (28) and (29), so (7) became non linear equation. In [32] we developed Lagrangian formulation of the theory. Excluding conjugate momenta from LH , we obtained the Lagrangian action. Further, excluding the auxiliary variables, one after another, we obtained various equivalent formulations of the model. In the end, we got the ”minimal” formulation without auxiliary variables. This reads2 Z √ √ √ 1 S = dτ − mc −xN ˙ x˙ + a3 ωN ˙ ω˙ − g4 (ω 2 − a4 ), (9) 2 µ

ν

where N µν ≡ η µν − ωωω2 is projector on the plane transverse to the direction of ω µ . The action is written in a dt parametrization τ which obeys dτ > 0, this implies g1 (τ ) > 0 and p0 > 0. Interaction with external background should not spoil the number and algebraic properties of constraints (1)-(3). We do not know how to achieve this for the minimal action3 . For instance, the natural reparametrization-invariant interaction ec Aµ x˙ µ + µFµν ω µ ω˙ ν , even for vanishing magnetic moment, leads to the theory with the number and algebraic structure of constraints different from those of free theory. So we start with the equivalent Lagrangian with four auxiliary variables, g1 , g3 , g4 and g7 , this turns out to be appropriate to our aims. Of course, the auxiliary variables will be excluded from final equations of motion, see Eqs. (45)-(51). To introduce coupling of the position variable with an electro-magnetic field, we add the minimal interaction term Aµ x˙ µ . As for spin, we propose to modify derivative of ω as follows ω˙ µ → Dω µ = ω˙ µ − g1

2 3

eµ (F ω)µ . c

(10)

The last term in (9) represents kinematic (velocity-independent) constraint which is well known from classical mechanics. So, we might follow the classical-mechanics prescription to exclude g4 as well. But this would lead to lose of manifest covariance of the formalism. Hanson an Regge [22] have found highly non linear interaction for the case of their relativistic top. It would be interesting to apply their formalism to our minimal action.

4 This is the only term which we have found to be consistent with the constraints Ti . Lagrangian reads L=

1 e g4 g1 g3 [g (xN ˙ x) ˙ − 2g7 (xN ˙ Dω) + g1 (DωN Dω)] + Aµ x˙ µ − (ω 2 − a4 ) − m2 c2 + a3 , ˜ 3 c 2 2 2 2 det G

(11)

˜ = g1 g3 − g 2 . Since L contains the auxiliary variables, even for µ = 0 we have higher nonlinear interaction. where detG 7 As a consequence, motion of spin influences motion of the particle and vice versa. We first establish whether our Lagrangian gives the desired constraints. The momenta read 1 e ∂L = (g3 N x˙ µ − g7 N Dω µ ) + Aµ , µ ˜ ∂ x˙ c det G 1 ∂L µ µ = (−g N x ˙ = + g N Dω ) , 7 1 ˜ ∂ ω˙ µ det G

pµ = πµ

πg i =

∂L = 0. ∂ g˙ i

(12)

(13)

According to (13), the momenta πgi represent primary constraints, πgi = 0. Using the property N ω = 0 of the projector N , from Eqs. (12) follow more primary constraints, T5 = πω = 0 and T6 = Pω = 0. It has been denoted e P µ = pµ − Aµ . (14) c To write Hamiltonian, we solve the system (12) with respect to projected velocities N x˙ µ = g1 P µ + g7 π µ ,

N Dω µ = g3 π µ + g7 P µ .

(15)

Using these expressions as well as the identities P x˙ = PN x, ˙ π ω˙ = πN ω˙ we obtain Hamiltonian H = px+π ˙ ω−L+λ ˙ a Φa in the form  g
g4 2
g1  2 µe 3 H= P − (JF ) + m2 c2 + π 2 − a3 + ω − a4 + λ5 (ωπ) + λ6 (Pω) + g7 (Pπ) + λgi πgi , (16) 2 2c 2 2 where λ5 and λ6 appear as Lagrangian multipliers for primary constraints T5 and T6 . We denote (JF ) = J µν Fµν and so on. From (16) we conclude that T1 , T3 T4 and T7 appear as secondary constraints when we impose the compatibility conditions π˙ gi = {πgi , H} = 0. The second, third and fourth stages of the Dirac-Bergmann algorithm can be resumed as follows T1 = 0 T3 = 0 T4 = 0

⇒ ⇒ ⇒

T5 = 0

⇒

T6 = 0

⇒

T7 = 0

⇒

λ6 C + g7 D = 0 , λ5 = 0 , λ5 = 0 , a3 a3 ⇒ λg4 = λg3 , g4 = g3 , a4 a4 g1 C − g7 M 2 c2 = 0 , ⇒ λg7 = f (λg1 ) , 2 2

g1 D + λ6 M c = 0 .

(17) (18) (19) (20) (21) (22)

We have denoted M 2 = m2 −

e(2µ + 1) Fµν J µν , 4c3

(23)

e eµ e eµ (ω∂)(JF ), D = − (µ − 1)(πF P) + (π∂)(JF ). (24) C = − (µ − 1)(ωF P) + c 4c c 4c Eq. (17) turns out to be a consequence of (21) and (22), λ6 (21) + g7 (22) = g1 (17), and can be omitted. Eq. (21) determines g7 = MC2 c2 g1 while (22) gives the lagrangian multiplier λ6 = − MD2 c2 g1 . The Dirac-Bergmann algorithm stops at the fourth stage. This yields all the desired constraints Ta , a = 1, 3, 4, 5, 6, 7. Two auxiliary variables, g1 and g3 , and the corresponding Lagrange multipliers λg1 , λg3 have not been determined. It is useful to summarize the algebra of Poisson brackets between constraints in a compact form, see Table I. We note that Poisson brackets of T1 and T˜3 = T3 + aa43 T4 vanish on the constraint surface, so they form the first-class subset. The presence of two first-class constraints is in correspondence with the fact that two lagrangian multipliers remain undetermined within the Dirac procedure. Matrix of Poisson brackets of the remaining constraints, T4 , T5 , T6 and T7 , is non degenerated, so this is a set of second-class constraints. All this is in correspondence with free theory [32]. In resume, interaction introduced does not spoil the structure and algebraic properties of Hamiltonian constraints of the free theory.

5 TABLE I: Algebra of constraints T1

T3

T4

T5

T6

T7

T1 = P − µe µν F Jµν + m2 c2 2c T3 = π 2 − a3

0

0

0

0

-2C

-2D

0

0

−4T5

−2(a3 + T3 )

−2T7

0

T4 = ω 2 − a4

0

4T5

0

2(T4 + a4 )

0

2T6

T5 = ωπ

0

0

−T6

T7

T6 = Pω

2C

2T7

0

T6

0

T1 − M 2 c2

T7 = Pπ

2D

0

−2T6

−T7

−T1 + M 2 c2

0

2

III.

2(T3 + a3 ) −2(a4 + T4 )

EXACT FRENKEL EQUATIONS ON ARBITRARY BACKGROUND A.

Hamiltonian equations of motion

The Hamiltonian (16) determines evolution of the basic variables through the Poisson bracket q˙ = {q, H}. Equivalently, we can pass from Poisson to Dirac bracket constructed on the base of the second-class constraints T4 , T5 , T6 , T7 . The list of Dirac brackets is presented in the Appendix. After that, our highly nonlinear interaction turns out to be hidden in the Dirac bracket: the constraints can be used in the Hamiltonian (16), this gives the expression  g
g1  2 µe 3 H1 = P − (JF ) + m2 c2 + J 2 − 8a3 a4 . (25) 2 2c 2 Equations of motion now can be obtained with help of H1 and the Dirac bracket, q˙ = {q, H1 }DB . They read µe e (26) x˙ µ = g1 uµ , P˙ µ = g1 (F u)µ + g1 ∂ µ (JF ), c 4c eµ ω˙ µ = g1 (F ω)µ + g3 π µ + g7 P µ , c a3 eµ π˙ µ = g1 (F π)µ − g3 ω µ − g6 P µ , c a4

(27)

where ∂ µ (JF ) = J αβ ∂ µ Fαβ . According to (21) and (22), the four-velocity uµ is not proportional to canonical momentum P µ g7 λ6 (28) uµ = P µ + π µ + ω µ = T µ ν P ν + Y µ . g1 g1 We have denoted e(µ − 1) µe (JF )µν , Yµ =− J µα ∂α (JF ) . (29) 2c3 M 2 8M 2 c3 All the basic variables have ambiguous evolution. xµ and P µ have one-parametric ambiguity due to g1 (they change under reparametrizations) while ω and π have two-parametric ambiguity due to g1 and g3 (they change under reparametrizations and spin-plane symmetry). The quantities xµ , P µ and the spin-tensor J µν are spin-plane invariants. Their equations of motion form a closed system   e(µ − 1) µe µ µα x˙ µ = g1 P µ + (JF P) − J ∂ (JF ) , (30) α 2c3 M 2 8M 2 c3 µe e ˙ µ + g1 ∂ µ (JF ) , (31) P˙ µ = (F x) c  4c  e µ−1 µ J˙µν = g1 µF [µ α J αν] − 2 2 (JF P)[µ P ν] − P [µ J ν]α ∂α (JF ) . (32) c c M 4M 2 c2 T µν = η µν +

6 The last term in (25) does not contributes to the equations of motion for x, P and J, and can be omitted. Then the Hamiltonian for these variables acquires a simple and expected form  g1  2 µe P − (33) H= (JF ) + m2 c2 . 2 2c The interaction yields two essential structural modifications of the theory. Free theory implies the Frenkel condition, J µν x˙ ν = 0, and pµ ∼ x˙ µ . Interaction modifies not only dynamical equations but also the Frenkel condition, the latter necessarily turns into the Pirani condition J µν Pν = 0 ,

(34)

where, due to (28), P is not proportional to x˙ µ . Then Eqs. (30) and (31) imply that the interaction leads to modification of the Lorentz-force equation even for uniform field. Only for the non anomalous value of magnetic moment, µ = 1, and uniform electromagnetic field the equations (21) and (22) would be the same as in free theory, λ6 = g7 = 0. Then T µν = η µν , Y µ = 0, and four-velocity becomes proportional to P µ . Contribution of anomalous magnetic moment µ 6= 1 to the difference between u and P is proportional to cJ3 ∼ c~3 , while the term with a gradient 2 2 of field is proportional to Jc3 ∼ ~c3 . The remaining ambiguity due to g1 presented in the equations (30)-(32) reflects the reparametrization symmetry of the theory. Assuming that the functions xµ (τ ), pµ (τ ) and J µν (τ ) represent the physical variables xi (t), pµ (t) and J µν (t) in the parametric form, their equations read dx0 ui dxi = c 0, = c, dt u dt dP µ e µe = 0 F µν uν + 0 ∂ µ (JF ) , dt u 4u c ˙µν dJ µν = J . dt g1 u0

(35) (36) (37)

As it should be, they have unambiguous dynamics. Equations (30)- (32) are written in arbitrary parametrization of the world-line. In the next subsection we exclude P µ and g1 , and then analyze the resulting equations in the proper-time parameterizations. This allow us to compare them with original Frenkel equations.

B.

Lagrangian form of equations

Hamiltonian equations from the previous section can be rewritten in the Lagrangian form for the set x, J. Let us analyze relation between velocity and momentum given by the Hamiltonian equation x˙ µ = g1 (T µ ν P ν + Y µ ) ,

(38)

where T µν = η µν − aJ µα Fα ν ,

a=−

e(µ − 1) . 2c3 M 2

(39)

Matrix T is invertible, the inverse matrix T˜ has the same structure (we used the identity (JF J)µν = − 21 (JF )J µν which implied by (4)) T˜µν = η µν + bJ µα Fα ν ,

b=

2a . 2 + a(JF )

(40)

Let us first exclude P and g1 from equation (31). From Eq. (38), we can express P through x˙ Pµ =

1 ˜µ ν T ν x˙ − T˜µ ν Y ν . g1

We can express g1 calculating square of the following expression 1 P µ + T˜µ ν Y ν = T˜µ ν x˙ ν , g1

(41)

7 which yields 1 P 2 + (T˜Y )µ (T˜Y )µ = 2 (T˜x) ˙ µ (T˜x) ˙ µ. g1 We used that Pµ T˜µ ν = Pµ and Pµ Y µ = 0. From the last equation we find g1 v u √ u (T˜x) ˙ 2 −g x˙ x˙ t g1 = ≡ , µe(JF ) 2 2 2 ˜ m rc (T Y ) − m c + 2c

(42)

where we have introduced the symmetric matrix gµν = (T˜T T˜)µν ,

(43)

and the radiation mass m2r = m2 −

gY Y µe (JF ) − 2 . 2c3 c

(44)

√ In the natural parametrization −g x˙ x˙ = c, we have g1 = m−1 r , that is the auxiliary variable, which appeared in front of mass-shell constraint T1 = 0, is the inverse radiation mass. Due to the identity (T˜Y )µ = bY µ we also can write gY Y = b2 Y 2 . Using (41) and (42) in (31) and (32) we write closed system of equations for xµ and J µν in the form " # √ d (T˜x) ˙ µ e µe −g x˙ x˙ µ µ µ ˜ mr c √ ˙ + ∂ (JF ) , (45) − (T Y ) = (F x) dτ c 4mr c2 −g x˙ x˙ √ e −g x˙ x˙ J˙µν = µF [µ α J αν] + c m c r √       −g x˙ x˙ ˜ [µ µ−1 mr c µ e ∂ (JF ) , (T˜x) ˙ [µ − (T Y ) J ν]α 2 2 √ F T˜x˙ − F T˜Y − α c mr c c M 4M 2 c2 −g x˙ x˙ α mr c α J µν T˜να ( √ x˙ − Y α ) = 0 . −g x˙ x˙

(46)

(47)

Let us compare them with Frenkel equations. Frenkel found equations of motion consistent with the condition J µν uν = 0 up to order O3 (J, F, ∂F ). Besides, in our phe considered the case µ = 1. Taking these approximations ˙ 2 = c, we arrive at those of Frenkel (our J is 2mc of Frenkel J) equations in the proper-time parametrization −(x) e i e d h e e e µ (m − (JF ))x˙ µ + J µα ∂α (JF ) = (F x) ˙ µ+ ∂ (JF ) , 3 2 3 dτ 4mc 8m c c 4mc   e 1 [µ ν]α J˙µν = F [µ α J αν] − x ˙ J ∂ (JF ) , α mc 4mc2

J µν x˙ ν = 0 .

(48)

(49)

In general case, our equations (45)-(47) involve two types of corrections as compared with those of Frenkel. First, energy of magnetic moment in non uniform field leads to the contribution − gYc2Y into the Frenkel radiation mass, see (43). Second, when µ 6= 0, a contribution arises because the Frenkel condition which has been satisfied for the free particle, turns into Pirani condition in the interacting theory. Its Lagrangian form is written in (47). In the result, the components J 0i vanish in the frame P µ = (P 0 , ~0) instead of the rest frame. Hence our model predicts small dipole electric moment of the particle. The structure of our equations simplified significantly for the stationary homogeneous field ∂α F µν = 0. In this case (45) and (46) read " # mr (T˜x) ˙ µ e p ˙ = 2 (F x) ˙ µ, (50) c −(g x˙ x) ˙

8 J˙µν =

" # p m2r (µ − 1) e αν] [µ [µ ˜ ν] ˜ p − µ −(g x˙ x)F ˙ (JF T x) ˙ (T x) ˙ . αJ mr c2 M 2 −(g x˙ x) ˙ (J T˜x) ˙ µ = 0.

Contract (50) with (T˜x) ˙ µ . Then r.h.s. vanishes. Using the identity vµ m ˙r

q

−(T˜x) ˙ 2 = 0,

(51)

(52) h

µ √v˙ −v 2

i

˙ ≡ 0 at the l.h.s., we obtain

or m ˙ r = 0,

(53)

on-shell. Hence F µν Jµν is a conserved quantity. Then T1 = 0 implies that P 2 is a conserved quantity as well. Contracting (50) with T we can further simplify this equation   e d mr c2 mr x˙ µ ˜˙ √ = 2 (F 0 x) T T. ˙ µ, F0 = TF − √ dτ c −g x˙ x˙ e −g x˙ x˙

(54)

Let us choose a parametrization which implies gµν x˙ µ x˙ ν = −c2 .

(55)

Since g x˙ x˙ = x˙ 2 + O(J 2 ), in the linear approximation on J this is just the proper-time parametrization. Then the equations (51) and (54) read e d(T˜x) ˙ µ = (F x) ˙ µ, dτ mr c J˙µν =

or,

x ¨µ =

e (F 0 x) ˙ µ, mr c

F0 = TF −

mr c ˜˙ T T, e

  e m2r (µ − 1) [µ ˜ ν] [µ αν] ˜ (JF T x) ˙ (T x) ˙ . µF α J − mr c c2 M 2

(56)

(57)

So, when µ 6= 0, the exact equations differ from the approximate equations (48) and (49) even for uniform field. C.

BMT vector in Frenkel theory

Since J µν Pν = 0, the spin-tensor is equivalent to the four-vector (7) where we replace pµ → P µ . Then sµ Pµ = 0. Due to Eqs. (28) and (29) together with (4), sµ also obeys the condition sµ uµ = sµ x˙ µ = 0.

(58)

The physical dynamics can be described using sµ instead of J µν . Eq. (58) suggests that sµ could be candidate for BMT-vector in our model. Using the identities J µν = √

2 µναβ sα Pβ , −P 2

µναβ Jαβ = √

4 P [µ sν] , −P 2

(59)

to represent J µν through sµ in Eqs. (30)-(32), we obtain the closed system of equations for spin-plane invariant quantities  e(µ − 1) µναβ µ √ x˙ = g1 P µ +  (F P)ν sα Pβ + 3 c M 2 −P 2 i µe µαγδ ρβγ 0 δ 0 0 Pδ 0 ∂α Fρβ ,   s P s (60) γ δ γ 2M 2 c3 P 2 µe e ˙ µ + g1 √ P˙ µ = (F x) αβγρ sγ P ρ ∂ µ F αβ , (61) c 2c −P 2   eµ 1 1 ˙ µ s˙ µ = g1 (F s)µ + 2 (sF P)P µ − 2 (Ps)P . (62) c P P These equations valid for arbitrary electro-magnetic fields. Let us consider the case of uniform field discussed by Bargmann Michel and Telegdi. Then we can compare these equations with BMT equations. First, we should exclude

9 P and g1 from equations (61) and (62) using (60). In contrast to (38), where x˙ µ is a linear function of P µ and J µν , in (60) x˙ µ is a non-linear function of P µ and sµ . Inverse function which express P µ as a function of x˙ µ , sµ exists, though we can’t find its explicit form even in the case of uniform fields. Formally using (41) and (42) in the case of uniform fields, ∂α F µν = 0, we get s˙ µ =

p eµ e √ −g x˙ x˙ (F s)µ − [(µ − 1)(sF x) ˙ + µb(sF JF x)] ˙ (T˜x) ˙ µ, mr c2 mr c2 −g xx ˙

sx˙ = 0 .

(63)

Eq. (63) contains J but for weak fields the corresponding contribution can be neglected. In the uniform field and in the parametrization (55) we have Eq. (56) for x and s˙ µ =

e eµ [(µ − 1)(sF x) ˙ + µb(sF JF x)] ˙ (T˜x) ˙ µ. (F s)µ − mr c mr c3

(64)

This can be compared with BMT-equations e (F x) ˙ µ, mc

(65)

e eµ (µ − 1)(sF x) ˙ x˙ µ . (F s)µ − mc mc3

(66)

x ¨µ =

s˙ µ =

We can also introduce BMT-tensor dual to sµ µν JBM T =

2 µναβ  sα x˙ β . c

Due to (66) this obeys the equation   e (µ − 1) αν] µν [µ ν] [µ J˙BM = (J F x) ˙ x ˙ . µF J − BM T α T BM T mc c2

(67)

This can be compared with (57). Obtaining their equation (66) in uniform field, Bargmann, Michel and Telegdi supposed that the motion of particle (65) is independent from the motion of spin. Besides they looked for the equation for sµ linear on s and F . Obtaining Eqs. (56) and (64) we have not made any supposition of such a kind. Our approach is based on the variational formulation which satisfies all the necessary symmetries. The exact equations (56) and (64) involve two types of essential corrections as compared with BMT equations. First, an energy of magnetic moment in electromagnetic fields leads to the radiation mass mr . Second, anomalous magnetic moment affects trajectory of a particle. Neglecting non linear on F and s terms in our equations (56), (64), we obtain those of Bargmann, Michel and Telegdi. The same holds if we take the proper-time parametrization, x˙ µ x˙ µ = −c2 , instead of (55). IV.

EXACT SOLUTION IN UNIFORM MAGNETIC FIELD

BMT equations give important information about spin kinematics for relativistic particles. Exact solutions to the Dirac equation in a constant magnetic field can be used to obtain solutions to BMT equations [39]. Integrability of BMT equations in the case of rather general electromagnetic backgrounds were studied in [40]. In this section we would like to study the behavior of our model in more details. We take BMT vector of our model as the basic quantity, to compare dynamics of our model with those of BMT. So we take for the analysis the closed system of equations (60). BMT equations were derived as approximate equations for the description of polarization effects in accelerators. They do not include possible influence of spin to the Lorentz force, and are linear in spin, field and anomalous magnetic moment. Therefore, comparing our results with BMT we will take into account this approximations. Consider a particle with initial momentum P µ (0) and BMT spin sµ (0) moving in the uniform magnetic field directed along z-axis, B = Bez , of a laboratory Cartesian coordinate system defined by an orthonormal basis (ex , ey , ez ). We already established that P 2 and (F J) are integrals of motion for uniform fields. In the case of uniform magnetic field we have (F J) = 4γ [(BS) − (Bβ)(βS)] . Initial values P µ (0), sµ (0) should satisfy to the constraints of our model (T1 = 0, s2 = 3~2 /4, Ps = 0).

10 Here and through the rest of this section we use following notations P0 γ=√ , −P 2

β=

~ P , P0

in accordance with our construction of Lorentz invariant SO(3) spin fiber bundle [31]. Here, γ plays a role of relativistic factor which in the limit of free electron reads γ = (1 − v2 /c2 )−1/2 . Denote by β module of vector β. The quantity (39), a=−

4m2 c3

2e(µ − 1) , − e(2µ + 1)(F J)

is also a constant, which practically (for the magnetic fields smaller than Schwinger field) can be taken as a ≈ − e(µ−1) 2m2 c3 . The Hamiltonian equations of motions (60)-(62) written in the parametrization of physical time read dx dt dP dt dS dt dS 0 dt u0

dx0 u , = c, u0 dt dP 0 e [u, B] , = 0, u0  dt  1 e eµ [S, B] + (S, P, B)P − 2 0 (S, u, B)P , u0 P2 P u eP 0 [µ(S, P, B) − (S, u, B)] , u0 P 2
 P 0 1 + 2aγ β 2 (BS) − (Bβ)(βS) , 2a P [1 + 2aγ ((BS) − (Bβ)(βS))] − B(PS) , γ

= c

(68)

=

(69)

= = =

u =

(70) (71) (72) (73)

where [S, B] and (S, P, B) mean the vector and mixed product of 3-dimensional vectors. The velocity is not collinear to the momentum P µ , its projection to the 3-hyperplane orthogonal to P µ is proportional to anomalous magnetic moment. As a result, in general there is no common rest frame for velocity and momentum. From (69) follows that P 0 = const, hence P 2 = const, γ = const and β 2 = const. Another integral of motion is the projection of momentum to the magnetic field, (PB) = const. To simplify our calculations we assume without loosing generality that initial vector of momentum is orthogonal to magnetic field, (PB) = 0. Indeed, other values of (PB) can be obtained by boosts along B which do not modify electromagnetic tensor (B0 = B, E0 = E = 0). For the motion with momentum orthogonal to magnetic field we obtain the following system of equations dP = Ωp [P, ez ] , dt   dx Ωp 2aΩs = c P− ez (PS) , dt eB γµe 0 dS = Ω0s (S, β, ez ) , dt dS = Ωs [S, ez ] + Ω0s β(S, β, ez ) , dt where we use the following constants (frequencies) eB(1 + 2aγ(BS)) , P 0 (1 + 2aγβ 2 (BS)) µeB = 0 , P (1 + 2aγβ 2 (BS)) eBP 0 (µ − 1 − 2aγ(BS)) = = γ 2 (Ωs − Ωp ) . P 2 (1 + 2aγβ 2 (BS))

(74) (75) (76) (77)

Ωp =

(78)

Ωs

(79)

Ω0s

(80)

Multiplying (77) by B we find that (BS) = const as it should be. From (74) we find that the vector P rotates with a constant circular frequency Ωp in the plane, orthogonal to magnetic field. In the orthonormal basis (ex , ey , ez ) solution for P yields P = |P (0) | (ex cos (Ωp t + φp ) + ey sin (Ωp t + φp )) .

(81)

11 For simplicity we choose ex to provide φp = 0. To solve (77) we first note that Sz = const. From constraints SP = 0, S 2 = 34 ~2 one can see that the spin S for a particle with momentum P = P 0 β belong to the following ellipsoid Si (δ ij − β i β j )Sj =

3 2 ~ , 4

(82)

√ which is obtained from a sphere with radius ~ 3/2 by stretching in γ times in the direction of P. The main principal axis of this ellipsoid is always directed along P. Therefore we write Sx ex + Sy ey = S1 τ 1 + S2 τ 2 where the new basis is τ 1 = ex cos(Ωp t) + ey sin(Ωp t) ,

then

dτ 1 = Ωp τ 2 , dt

τ 2 = −ex sin(Ωp t) + ey cos(Ωp t)) ,

then

dτ 2 = −Ωp τ 1 , dt

τ 3 = ez , which is determined by two principal axes of the ellipsoid. Note, that we just choose convenient variables to solve the differential equation without transition to another reference frame. As a result we obtain simple equations dS1 = −Ω0s S2 , dt

dS2 Ω0 = 2s S1 , dt γ

(83)

which describe evolution of spin on the ellipsoid (82). Auxiliary radius vector S1 τ 1 + S2 τ 2 rotates with circular frequency Ω0s γ −1 . Solutions to Eqs. (83) read  0   0  Ωs S (0) Ωs (0) S1 (t) = S cos t+φ , S2 (t) = sin t+φ . (84) γ γ γ ~ moves on the ellipse which is obtained as intersection of ellipsoid (82) and plane Sz = const. S (0) The radius vector S √ is nothing but the semi-major axis of this ellipse, therefore it is restricted by the following interval 0 ≤ S (0) ≤ γ 3~/2. In terms of initial variables spacial components of 4-vector S µ evolves as follows   0   0   Ωs Ωs 1 S(t) = ex S (0) cos t + φ cos (Ωp t) − sin t + φ sin (Ωp t) + γ γ γ

ey S

(0)

  0   0   Ωs 1 Ωs cos t + φ sin (Ωp t) + sin t + φ cos (Ωp t) + ez Sz(0) , γ γ γ

(0)

where constants S (0) , Sz , φ are restricted by (82). Note that the angular velocity of precession of vector S 0 around B is time-dependent.   Nevertheless, the helicity (PS) changes with the constant rate, Ωs /γ. Indeed, (PS) = S (0) |P (0) | cos

Ω0s γ t

+φ .

Now we can substitute solutions S(t) and P(t) into equation for ~x   0  dx Ωp 2aΩs (0) Ωs (0) = c|P | (~ex cos(Ωp t) + ~ey sin(Ωp t)) − S ~ez cos t+φ . dt eB γµe γ Integrating the last equation we find the trajectory

c|P

(0)

 |

x(t) = ~xc +  0  1 2aΩs Ωs (~ex sin(Ωp t) − ~ey cos(Ωp t)) − 0 S (0)~ez sin t+φ . eB Ωs µe γ

(85)

Trajectory represents sum of two motions: circular motion in the plane orthogonal to B and oscillations along B. These oscillations accompany variations of the helicity. Vector xc defines the center of circle. The amplitude of oscillations along B ∆z = −

2caΩs (0) (0) |P |S . βλC , Ω0s µe

(86)

12

B

S

~blc

v S

P

v

P

FIG. 1: Momentum, velocity, spin and trajectory of a charged spinning particle in the uniform magnetic field

less than the Compton wave-length. The trajectory of the particle is shown in Figure 1. Exact solution for this particular case demonstrates some important features of our model. There are two essentially different situations with µ = 1 and µ 6= 1. In the case of usual magnetic moment µ = 1, helicity is an integral of motion, additional oscillations vanish and the particle moves along circular trajectory in the plane orthogonal to the magnetic field (dotted line in Figure 1). In the case of anomalous magnetic moment µ − 1 6= 0, helicity oscillate and affects the trajectory of particle. For a small anomalous magnetic moment helicity and z-coordinate oscillate with a very slow rate. These slow oscillations of trajectory along B with the amplitude of Compton wavelength can be called magnetic Zitterbewegung. Magnetic Zitterbewegung in our model appears due to modification of Lorentz force for the spinning particle. Correction to the angular velocity of orbital motion Ωp is given by   eB eγ(µ − 1) Ωp ≈ 1+ (SB) + o(~, µ − 1, B) , (87) γmc m2 c3 where eB/(γmc) is the angular velocity for spinless or BMT particle. Frequency of helicity variations also corrected by high-order terms  Ω0s eB  eγ = (µ − 1) 1 − 2 3 (1 + (µ − 1)β 2 )(SB) + o(~, µ − 1, B) , γ mc m c

(88)

eB eB from the value (µ − 1) mc = ( g2 − 1) mc computed by Bargmann, Michel and Telegdi [20]. The corrections are small, and for the experiments discussed by them, our equations gives practically the same results. Therefore our model is compatible with these experiments. We believe that other physical situations may be realized, where the corrections could became notable. For instance, this may be the case of quasiparticles with large magnetic moment [17].

V.

CONCLUSIONS

In this work we have presented solution to the problem which has been posed by Frenkel in 1926. He noticed that search for variational formulation which takes into account the spin-tensor constraint J µν x˙ ν = 0 represents rather non trivial problem. He found equations of motion consistent with this condition in the approximation O3 (J, F, ∂F ), and when anomalous magnetic moment is µ = 1. We have found Lagrangian action (11) for charged spinning particle which implies all the desired constraints and equations of motion without approximation. They remain consistent for any value of magnetic moment and for an arbitrary electromagnetic background. Besides, due to the constraints

13 (1)-(3), our action guarantees the right number of both spacial and spin degrees of freedom. In the above mentioned approximations, our equations coincide with those of Frenkel. In the recent work [35], we also demonstrated that the classical spinning particle has an expected behavior in arbitrary curved background. With the Lagrangian and Hamiltonian formulations at hands, we can unambiguously construct quantum mechanics of the spinning particle and establish its relation with the Dirac equation. For the free theory, this has been done in the work [32]. We showed that this gives one-particle sector of the Dirac equation. Due to the second-class constraints (3), the positions xi obey to classical brackets with nonvanishing right hand side, see (90). So, in the Dirac theory they realized by non commutative operators which we identified with Pryce (d) center-of-mass [1]. Since namely xi has an expected behavior (56) as the position of spinning particle in classical interacting theory, our model argue in favor of covariant Pryce (d) operator as the position operator of Dirac theory. In resume, we have constructed variational formulation for relativistic spin one-half particle which is self consistent and has reasonable behavior on both classical and quantum level. As we have seen, interaction necessarily modifies some basic relations of the model. In the free theory the conjugated momentum is proportional to velocity, pµ ∼ x˙ µ and the Frenkel condition holds. This is no more true in interacting theory. The Frenkel condition turns into the Pirani condition, J µν Pν = 0, where the canonical momentum is not collinear to velocity. The advantage of Lagrangian formulation is that this gives exact relation between them (see also Eqs. (38), (39) and (42) ) x˙ µ = g1 (T µ ν P ν + Y µ ) ,

T µν = η µν + O(µ − 1, J),

Y µ = O(∂F, J) .

(89)

Only when µ = 1 and ∂α Fµν = 0, the interacting and free theory have the same structure. To see the meaning of the deformation, we can compare our Hamiltonian equations for Y µ = 0: x˙ µ = g1 (T µ ν P ν + Y µ ), P˙ µ = g1ce (F P)µ with the standard expressions x˙ µ = g1 P µ , P˙ µ = g1ce (F P)µ . From this it is clear that excluding P from our equations, we obtain extra contributions to the standard expression for the Lorentz force, x ¨µ = g1ce (F x) ˙ µ + O(J). So the modification (89) mean that complete theory yields an extra spin-orbit interaction as compared with the approximate Frenkel and BMT equations. We studied possible effects of this spin-orbit interaction in the case of uniform magnetic field. The exact analytical solution was obtained. Besides oscillations of the helicity first calculated by Bargmann, Michel and Telegdi, the particle with anomalous magnetic moment experiences effect of magnetic Zitterbewegung of the trajectory. Usual circular motion in the plane orthogonal to B perturbed by slow oscillations along B with the amplitude of order of Compton wavelength. The Larmor frequency (87) and the frequency of helicity oscillations (88) are also shifted by small corrections. It would be interesting to construct an experiment which could detect these possible corrections, for instance due to resonance effects. Another possibility is an artificial simulation of a point-like system with spin and a large anomalous magnetic moment. This can be inspired by simulations of Zitterbewegung itself with a trapped ions [12]. VI.

ACKNOWLEDGMENTS

The work of AAD has been supported by the Brazilian foundation CNPq. AMPM thanks CAPES for the financial support (programm PNPD/2011).

[1] M. H. L. Pryce, The mass-centre in the restricted theory of relativity and its connexion with the quantum theory of elementary particles, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences 195 (1948) 62. [2] T. D. Newton and E. P. Wigner, Localized states for elementary systems, Rev. Mod. Phys. 21 (1949) 400. [3] L. L. Foldy and S. A. Wouthuysen, On the Dirac theory of spin 1/2 particles and its non-relativistic limit, Phys. Rev. 78 (1950) 29. [4] L. E. Kidder, C. M. Will and A. G. Wiseman, Phys. Rev. D 47 (1993) R4183. [5] I. B. Khriplovich and A. A. Pomeransky, Gravitational Interaction of Spinning Bodies, Center-of-Mass Coordinate and Radiation of Compact Binary Systems, Phys. Lett. A 216 (1996) 7. [6] I. B. Khriplovich and A. A. Pomeransky, Equations of motion of spinning relativistic particle in external fields, J. Exp. Theor. Phys. 86 (1998) 839. [7] A. A. Pomeransky and R. A. Sen’kov, Quadrupole interaction of relativistic quantum particle with external fields, Phys. Lett. B 468 (1999) 251. [8] Y. N. Obukhov and D. Puetzfeld, Dynamics of test bodies with spin in de Sitter spacetime, Phys. Rev. D 83 (2011) 044024. [9] R. L. Jaffe and A. Manobar, Nuc. Phys. B 337 (1990) 509.

14 [10] M. Wakamatsu, Eur. Phys. J. A 44 (2010) 297. [11] A. Selva, J. Magnin, and L. Masperi. Bargmann-michel-telegdi equation and one-particle relativistic approach, Il Nuovo Cimento B Series 11, 111(7):855–862, 1996. [12] F. Zahringer, G. Kirchmair, R. Gerritsma, E. Solano, R. Blatt, and C. F. Roos, Nature (London) 463 (2010) 68. [13] A. A. Deriglazov, Classical-mechanical models without observable trajectories and the Dirac electron, Phys. Lett. A 377 (2012) 13. [14] J. H. Field, E. Picasso and F. Combley, Tests of fundamental physical theories from measurements of free charged leptons, Sov. Phys. Uspekhi, 22 (1979) 199 (in Russian). [15] J. P. Miller, E. de Rafael, and B. L. Roberts. Muon (g-2): experiment and theory, Reports on Progress in Physics, 70 (2007) 795. [16] G. H. Hoffstaetter, H. S. Dumas, and J. A. Ellison. Adiabatic invariance of spin-orbit motion in accelerators, Phys. Rev. ST Accel. Beams, 9 (2006) 014001. [17] D. V. Karlovets. Electron with orbital angular momentum in a strong laser wave, Phys. Rev. A 86 (2012) 062102. [18] J. Frenkel, Die elektrodynamik des rotierenden elektrons, Zeitschrift f¨ ur Physik, 37(4-5) (1926) 243. [19] J. Frenkel, Spinning electrons, Nature 117 (1926) 653. [20] V. Bargmann, L. Michel, and V. L. Telegdi, Phys. Rev. Lett. 2, 435 (1959). [21] H. C. Corben. Classical and quantum theories of spinning particles, Holden-Day, San Francisco, 1968. [22] A. J. Hanson and T. Regge, The relativistic spherical top, Annals of Physics, 87(2) (1974) 498. [23] F. A. Berezin and M. S. Marinov, Particle spin dynamics as the Grassmann variant of classical mechanics, Ann. Phys. 104 (1977) 336. [24] A. P. Balachandran, G. Marmo, A. Stern and B. S. Skagerstam, Phys. Lett. B 89 (1980) 199. [25] A. O. Barut and W. Thacker, Phys. Rev. D 31 (1985) 1386. [26] A. O. Barut and A. J. Bracken, Zitterbewegung and the internal geometry of the electron, Phys. Rev. D 23 (1981) 2454. [27] P. Grassberger, Classical charged particles with spin, Journal of Physics A: Mathematical and General, 11(7) (1978) 1221. [28] G. Cognola, L. Vanzo, S. Zerbini, and R. Soldati, On the lagrangian formulation of a charged spinning particle in an external electromagnetic field, Physics Letters B 104 (1981) 67. [29] B. P. Kosyakov, On the inert properties of particles in classical theory, Fiz. Elem. Chast. Atom. Yadra 34 (2003) 1564 (Engl. transl.: Phys. Part. Nucl. 34 (2003) 808). [30] A. A. Deriglazov, Variational problem for the Frenkel and the Bargmann-Michel-Telegdi (BMT) equations, Mod. Phys. Lett. A 28 (2013) 1250234; arXiv:1204.2494. [31] A. A. Deriglazov and A. M. Pupasov-Maksimov, Geometric constructions underlying relativistic description of spin on the base of non-Grassmann vector-like variable, arXiv:1311.7005 [32] A. A. Deriglazov and A. M. Pupasov-Maksimov, Lagrangian for Frenkel electron and position’s non-commutativity due to spin, arXiv:1312.6247. [33] R. P. Feynman and M. Gell-Mann, Theory of the Fermi interaction, Phys. Rev. 109 (1958) 193. [34] R. P. Feynman, Quantum electrodynamics, W A Benjamin, 1961. [35] W. G. Ramirez, A. A. Deriglazov and A. M. Pupasov-Maksimov, Frenkel electron and a spinning body in a curved background, arXiv:1311.5743. [36] F. A. E. Pirani, Acta. Phys. pol. 15 (1956) 389. [37] W. G. Dixon, Nuovo. Cim. 34 (1964) 317. [38] W. Tulczyjew, Acta Phys. Pol. 18 (1959) 393. [39] V. A. Bordovitsyn and A. N. Myagkii. Verification of the semiclassical method for an electron moving in a homogeneous magnetic field, Phys. Rev. E 64 (2001) 046503. [40] A. E. Lobanov and O. S. Pavlova. Solutions of the classical equation of motion for a spin in electromagnetic fields, Theor. Math. Phys. 121 (1999) 1691.

Appendix. Dirac brackets

We construct Dirac brackets that take into account the second-class pairs T3 , T4 , T5 and T6 . We will calculate them iteratively in the case of arbitrary electromagnetic background. Then Dirac brackets of the free theory can be obtained by substitution F µν = 0. We start from the pair of second class constraints is T6 and T7 , 467 = {T6 , T7 } = P 2 +

e Fµν J µν . 4c

At the constraint surface 467 = −M 2 c2 . The Poisson brackets of initial variables with constraints T6 and T7 are given in table II. Using table II we calculate the Dirac brackets of basic variables with respect to constraints T6 and T7 {Q1 , Q2 }67 = {Q1 , Q2 } +

1 ({Q1 , T6 }{T7 , Q2 } − {Q1 , T7 }{T6 , Q2 }) . 467

15 TABLE II: Constraints vs variables xµ

{, }

Pµ

πµ

ωµ

J µν

T6 = Pω −ω µ − ec F µν ων P µ 0 2(ω µ P ν − ω ν P µ ) T7 = Pπ −π µ − ec F µν πν 0 −P µ 2(π µ P ν − π ν P µ )

The brackets read {xµ , xν }67 =

−J µν , 2467

{xµ , P ν }67 = η µν +

(90)

e µν , J µα Fα ν ≡ T(0) 2c467

−ω µ P ν , 467 −π µ P ν = , 467
1 = J µα P β − J µβ P α , 467

(91)

{xµ , ω ν }67 =

(92)

{xµ , π ν }67

(93)

{xµ , J αβ }67

{ω µ , ω ν }67 = 0 , {π µ , π ν }67 = 0 ,

(95) (96)

P µP ν ≡ Gµν , 467 = 2(ω α Gµβ − ω β Gµα ) , = 2(π α Gµβ − π β Gµα ) , = 2(Gµα J νβ − Gµβ J να − Gνα J µβ + Gνβ J µα ) ,

{ω µ , π ν }67 = η µν − {ω µ , J αβ }67 {π µ , J αβ }67 {J µν , J αβ }67

(94)

e2 e e µν αν F + (F JF )µν = F µ α T(0) , c 2467 c2 c e F µα ωα P ν , = − 467 c e = − F µα πα P ν , 467 c e F µν (P α Jν β − P β Jν α ) . = − 467 c

(97) (98) (99) (100)

{P µ , P ν }67 =

(101)

{P µ , ω ν }67

(102)

{P µ , π ν }67 {P µ , J αβ }67

(103) (104)

µν We have defined define tensor Gµν as the Dirac bracket of spin variables ω µ and π ν . Besides, T(0) = T µν (µ = 0), µν where T µν = η µν − e(µ−1) 2c467 (JF ) . On the next step we calculate Dirac brackets for the pair {T4 , T5 }67 = 2(T4 + a4 ),

{Q1 , Q2 }4567 = {Q1 , Q2 }67 +

1 ({Q1 , T4 }67 {T5 , Q2 }67 − {Q1 , T5 }67 {T4 , Q2 }67 ) . 2ω 2

The Dirac brackets {, }67 of initial variables with T4 and T5 are given in table III. From table III it is seen that variables xµ , P µ , J µν have vanishing Dirac brackets {, }67 with constraints T4 and T5 . Therefore, new Dirac brackets {, }4567 coincide with old Dirac brackets {, }67 when at least one of arguments is a function Z(xµ , P µ , J µν ) of variables xµ , P µ and J µν only, {Z, Q}4567 = {Z, Q}67 ,

Z = Z(xµ , P µ , J µν ) .

(105) (106)

16 TABLE III: Constraints vs variables {, }67

xµ P µ π µ

T4 = ω 2 − a4 0 T5 = ωπ 0

ω µ J µν

0 2ω µ 0 0 π µ −ω µ

0 0

We omit subscripts of brackets, so that { , } means { , }4567 . The only modification in the Dirac brackets comes from the basic variables in the spin sector {ω µ , ω ν } = 0 , {ω µ , π ν } = η µν {π µ , π ν } =

(107) ωµ ων P µP ν , − − 467 ω2

−J µν . 2ω 2

(108) (109) (110)

Thus the complete list of Dirac brackets {, } consist of expressions (90)-(94), (107)-(109) and (98)-(104). Now the second class constraints can be put equal to zero, therefore we can rewrite the Hamiltonian as follows H=

g1 2 eµ (P − (F J) + m2 c2 ) . 2 2c

(111)

The constraint T3 can also be excluded since it has zero Dirac brackets with T1 and with all spin-plane invariant variables of the theory. This Hamiltonian generates evolution P˙ µ = {P µ , H} ,

x˙ µ = {xµ , H} ,

J˙µν = {J µν , H} .

To check consistency of our calculations, let us obtain equations of motion using the Dirac brackets {, } (and taking into account that at the constrained surface 467 = −M 2 c2 = −m2 c2 − e(2µ+1) 4c3 (F J) ). Equation for coordinate reads x˙ µ =

g1 {xµ , P ν }Pν − g1  g1 η µν −

g1 µ 2 eµ {x , P − (F J)} = 2 2c

eµ µ αβ eµ {x , J }Fαβ − g1 {xµ , F αβ }Jαβ = 4c 4c


g1 eµ g1 eµ µρ e J µα Fα ν Pν − P α J µβ − P β J µα Fαβ − J ∂ρ (F J) = 2 3 2 3 2M c 4M c 8M 2 c3

 g1 η µν −

 e eµ eµ ν µα J F Pν + g1 J µβ Fβα P α − g1 J µρ ∂ρ (F J) = α 2M 2 c3 2M 2 c3 8M 2 c3 

g1 η

µν

 e(µ − 1) µα ν eµ J Fα Pν − g1 J µρ ∂ρ (F J) = g1 uµ , + 2M 2 c3 8M 2 c3

Equation for momentum reads g1 eµ P˙ µ = {P µ , P 2 − (F J)} = 2 2c = g1 {P µ , P ν }Pν − g1  = g1

eµ µ αβ eµ {P , J }Fαβ − g1 {P µ , xρ }∂ρ F αβ Jαβ = 4c 4c

 e µν e2 e2 µ µα νβ F + F F J P + g F µν P α J β ν Fαβ + αβ ν 1 c 2M 2 c4 2M 2 c4

17 g1

e = g1 F µα c

 ηα

ν

 eµ  µρ e η − J ρα Fα µ ∂ρ (F J) = 2 3 4c 2M c

  e(µ − 1) eµ eµ βν ρ Jαβ F Jαρ ∂ (F J) + g1 ∂ µ (F J) = + Pν − 2 3 2 3 2M c 8M c 4c e eµ = g1 F µα uα + g1 ∂ µ (F J) , c 4c

Equation for spin-tensor reads eµ g1 (F J)} = J˙αβ = {J αβ , P 2 − 2 2c = g1 {J αβ , P µ }Pµ − g1

= g1

e M 2 c3

F µν (P α J β ν − P β J α ν )Pµ + g1

g1

= −g1

eµ µα νβ (G J − Gµβ J να )Fµν + c


eµ 1 P α J µβ − P β J µα ∂µ (F J) = 2 2 4c M c

e eµ eµ P [α J β] ν F νµ Pµ + g1 (F α ν J νβ − F β ν J να ) + g1 2 3 P [α J β] ν F νµ Pµ + M 2 c3 c M c

g1

e = g1 c

eµ αβ µν eµ {J , J }Fµν − g1 {J αβ , xµ }∂µ (F J) = 4c 4c


eµ 1 P α J µβ − P β J µα ∂µ (F J) = 2 2 4c M c

  1 (µ − 1) [α β] νµ [α β] µ α νβ β να P J ν F Pµ − P J µ ∂ (F J) . µ(F ν J − F ν J ) + M 2 c2 4M 2 c2

The Hamiltonian is proportional to the first class constraint T1 . Therefore equations of motion contain arbitrary function g1 (τ ) which is related with reparameterization invariance of the model. To obtain unambiguous equations of evolution we can impose the gauge x0 = cτ . The gauge is often called canonical gauge. Constraint T1 together with this condition form a pair of second class constraints. We have {x0 − cτ, T1 } = 2u0 , and Dirac brackets in the canonical gauge read {Q1 , Q2 }τ = {Q1 , Q2 } +

1 ({Q1 , G1 }{T1 , Q2 } − {Q1 , T1 }{G1 , Q2 }) . 2u0

The Dirac brackets of constraint T1 and canonical gauge G1 with physical variables are given in table IV. There compact notations g1 P˙ µ ≡ {P µ , T1 } , 2 g1 J˙µν ≡ {J µν , T1 } , 2 uµ ≡

  e(µ − 1) µα ν eµ η µν + J F Pν − J µρ ∂ρ (F J) . α 2M 2 c3 8M 2 c3

18 TABLE IV: Constraints vs variables (canonical gauge) {, }4567 T1 = P 2 − 0

eµ (F J) 2c

G1 = x − cτ

xµ −2uµ

Pµ −2 ˙ µ P

−1 J 0µ 2467

0µ T(0)

g1

ωµ µ

2CP 467

πµ µ

− 2eµ (F ω)µ 2DP c 467 −ω 0 P µ 467

J µν −2 ˙µν J

− 2eµ (F π)µ c g1 −π 0 P µ −1 J 0[ν 467 467

P µ]

eµ e (ω∂)(JF ) , C = − (µ − 1)(ωF P) + c 4c e eµ D = − (µ − 1)(πF P) + (π∂)(JF ) , c 4c were used. The Dirac brackets which take into account the canonical gauge are as follow. Spacial sector:
−1 u0 J µν − uµ J 0ν + uν J 0µ , {xµ , xν }τ = 0 2u 467
uµ e eµ {xµ , P ν }τ = η µν − 0 η 0ν + 0 u0 J µα − uµ J 0α + uα J 0µ Fα ν − 0 J 0µ ∂ ν (F J) , u 2u c467 8u c467   c ˙ µ 0ν c ˙ µ 0ν e µ ν 0 µ αν {P , P }τ = 0 u F α T(0) − P T(0) + P T(0) . u c eg1 eg1

(112) (113) (114)

Frenkel sector:   1 {J µν , J αβ }τ = {J µν , J αβ } − J˙µν J 0[β P α] − J˙αβ J 0[ν P µ] , 0 g1 u 467   −1 1 {xµ , J αβ }τ = 0 uµ J 0[α P β] − u0 J µ[α P β] − 0 J 0µ J˙αβ , u 467 2u 467 g1   e eµ 1 0µ ˙αβ {P µ , J αβ }τ = 0 F µ ν u0 P [α J β]ν − uν P [β J α]0 − 0 T(0) J ∂ µ (F J)P [β J α]0 . u c467 u g1 4u0 c467

(115) (116) (117)

Basic spin variables: eµω 0 (F µα P ν − F να P µ )ωα , 2u0 c467  
P µP ν eµω 0 (π 0 C + ω 0 D) ωµ ων µν = η − − π 0 F µα ωα P ν − ω 0 F να πα P µ , 1+ − 467 u0 c467 ω2 u0 c467 −J µν eµπ 0 = − (F µα P ν − F να P µ )πα . 2ω 2 2u0 c467

{ω µ , ω ν }τ = −

(118)

{ω µ , π ν }τ

(119)

{π µ , π ν }τ

(120)

Other mixed brackets: µ

ν

{x , ω }τ {xµ , π ν }τ {ω µ , J αβ }τ {π µ , J αβ }τ {P µ , ω ν }τ {P µ , π ν }τ

  0µ   uµ ω 0 P ν 1 J CP ν eµ −ω µ P ν ν (F ω) + + 0 − , = 467 u 2467 467 c 467     −π µ P ν 1 J 0µ DP ν eµ uµ π 0 P ν = + 0 − (F π)ν + , 467 u 2467 467 c 467 !  ˙αβ  CP µ 1 eµ 0 µJ µ 0[β α] α µβ β µα = 2(ω G − ω G ) − 0 ω P + − (F ω) J P , u 467 g1 467 c !  ˙αβ  DP µ 1 eµ α µβ β µα 0 µJ µ 0[β α] = 2(π G − π G ) − 0 π P + − (F π) J P , u 467 g1 467 c !   e 1 CP ν eµ P˙ µ ω 0 P ν 0µ µα ν ν = − F ωα P − 0 T(0) − (F ω) − , 467 c u 467 c g1 467 !   e 1 DP ν eµ P˙ µ π 0 P ν 0µ µα ν ν = − F πα P − 0 T(0) − (F π) − . 467 c u 467 c g1 467

(121) (122) (123)

(124)

(125)

(126)

19 In the free theory the algebra of Dirac brackets simplifies significantly. In this case P µ = pµ , uµ = pµ , J˙µν = P˙ µ = 0, 467 = p2 , and in an arbitrary parametrization τ , we have the following brackets Basic variables of spin: {ω µ , ω ν } = 0,

{ω µ , π ν } = g µν −

{xµ , ω ν } = −

ω µ pν , p2

ωµ ων , ω2

{π µ , π ν } = −

{xµ , π ν } = −

1 µν J ; 2ω 2

π µ pν ; p2

(127)

(128)

Spacial sector: {xµ , xν } = −

1 µν J , 2p2

{xµ , pν } = η µν ,

{pµ , pν } = 0.

(129)

Frenkel sector: {J µν , J αβ } = 2(g µα J νβ − g µβ J να − g να J µβ + g νβ J µα ) ,

{xµ , J αβ } = BMT-sector: take sµ = √1 4

−p2

1 µα β (J p − J µβ pα ) , p2

(130)

(131)

µναβ pν Jαβ , then {sµ , sν } = − p

{xµ , sν } = −

1 −p2

µναβ pα sβ =

1 µν J , 2

(132)

sµ pν 1 pµ sν =− p µναβ Jαβ − 2 . 2 p p 4 −p2

(133)

Other Dirac brackets vanish. In the equations (127) and (141) it has been denoted gµ ν ≡ δµ ν − Together with g˜µ ν ≡

pµ pν p2 ,

pµ pν . p2

(134)

this forms a pair of projectors g + g˜ = 1,

g 2 = g,

g˜2 = g˜,

g˜ g = 0.

(135)

The free Dirac brackets which take into account the canonical gauge are as follows. Basic variables of spin: {ω µ , ω ν }τ = 0,

{ω µ , π ν }τ = g µν −

ωµ ων , ω2

{π µ , π ν }τ = −

−ω µ pν ω 0 pµ pν + 0 2 , 2 p p p µ ν 0 µ ν −π p π p p + 0 2 , = p2 p p

1 µν J , 2ω 2

(136)

{xµ , ω ν }τ =

(137)

{xµ , π ν }τ

(138)

Spacial sector:
−1 p0 J µν − pµ J 0ν + pν J 0µ , 2p0 p2 pµ = η µν − 0 η 0ν , {pµ , pν } = 0 . p

{xµ , xν }τ =

(139)

{xµ , pν }τ

(140)

20 Frenkel sector: {J µν , J αβ }τ = 2(g µα J νβ − g µβ J να − g να J µβ + g νβ J µα ) , {xµ , J αβ }τ =

 −1  µ 0[α β] p J p − p0 J µ[α p β] , 0 2 p p

(141)

(142)

BMT-sector: 1 1 {sµ , sν }τ = − p µναβ pα sβ = J µν , 2 −p2

(143)

1 pµ (s0 pµ pν − sµ pν p0 ) p µναβ Jαβ + 0ναβ Jαβ = {xµ , sν }τ = − p . p0 p2 4 −p2 4p0 −p2

(144)

Other Dirac brackets vanish. Here we define Dirac brackets for all phase space variables. After transition to the Dirac brackets the second-class constraints can be used as strong equalities, therefore it is enough to consider Dirac brackets at the constraint surface only. Then, explicit form of the Dirac brackets depends on the choice of independent variables. For instance, in the free theory considered in the gauge of physical time we can present s0BM T and p0 in terms of independent variables sBM T , p, x (s p) s0 = p , p2 + (mc)2

p0 =

p

p2 + (mc)2 .

(145)

The non vanishing Dirac brackets are {xi , xj }τ =

ijk sk , mcp0

{xi , pj }τ = δ ij ,

{pi , pj }τ = 0,

(146)

  p0 ijk (s p)pk {s , s }τ =  sk − , mc p20

(147)

pj . (mc)2

(148)

i

j

{xi , sj }D =



si −

(s p)pi p20