License: CC Zero
arXiv:2501.13110v1 [hep-ph] 07 Jan 2025
institutetext: National Research Nuclear University MEPhI,
Kashirskoye shosse, 31, 115409 Moscow,
Russia

The Open Relativistic Two-body Problem

A.V.Koshelkin and.kosh59@gmail.com
Abstract

The open relativistic two-body problem, when two interacting particles also are in external potentials, is considered in terms of the principle of the least action. Based on the consistent modification of the relativistic version Newton’s third law in external fields, the exact covariant operator equations, which govern dynamics of either a scalar particle or each of the components of the 16 component spinor of spin-1/2121/21 / 2 fermions, depending on the particle type, are derived in the center-mass and relative motion variables, beyond the consideration in the Breit frame only. The class of external fields and interaction potentials, when the two-body problem can be covariantly reformulated in (3+1) phase space of relative motion variables, uncoupled from the center-mass motion of such a system, is outlined. In the case of fermions the new γ𝛾\gammaitalic_γ-matrix basis generating the Dirac-like equation for the 16 component spinors in the (3+1) phase space is found. The chosen basis allows us to decouple the derived Dirac-like equation into four independent equations governing the dynamics of the four-component spinors. The developed approach is examined in studying the low-energy positronium states in magnetic fields.

1 Introduction

The relativistic two-body problem, which has actively studied for more than forty years, beginning from 80th Teo71 ; Teo76 ; Droz ; Kal ; Teo78 ; Kom ; Are ; Kir ; Dom78 ; Dom80 ; Sad86-1 ; Sad86-2 ; Sad94 ; Sad89 ; Cra84 ; Cra87 ; Cra90 ; Cra01 ; Cra03 ; Cra06 ; Cra10 ; Bar87 ; Bar91 ; Log , is of great interest now in connection primarily with the indirect searches of the dark matterAgu ; Perez ; Jung ; Ber ; Gel . It also concerns researches when strong magnetic fields play an important role in forming signals from both the astrophysics objectsBro ; Bus ; Dun ; Hard06 and the strong interacting matter, arising in non-central collisions of high energy particlesFuku ; Vor ; Tuc . The value of the magnetic field strength in the latter case is estimated to be about B10181019Gsimilar-to𝐵superscript1018superscript1019𝐺B\sim 10^{18}-10^{19}Gitalic_B ∼ 10 start_POSTSUPERSCRIPT 18 end_POSTSUPERSCRIPT - 10 start_POSTSUPERSCRIPT 19 end_POSTSUPERSCRIPT italic_G at RHIC and LHC energiesFuku ; Vor ; Tuc what is of the order and more than the squared pion mass. Besides that, magnetic fields play a key role in forming signals from the interior area of neutron starsDun ; Hard06 where the magnetic field strength turns out to be of the same order of magnitude.

The first attempt to describe the relativistic two-body system was made by summing two free Dirac equations, adding to them interaction potential Brei29 ; Brei30 , that resulted in the relativistic non-invariant Breit equation. On other hand, the manifestly covariant Bethe-Salpeter equationBSE , which is suitable for solving such kind problems, was found to haveNak the negative normalized solutions in the case of two one half spin particles. The lack of the approaches mentioned above can be eliminated with Dirac’s Hamiltonian constraint techniqueDir ; Ber , decreasing a number of the Hamiltonian variables by means of establishing relations between them. It is achieved due to the Poincare invariant functions containing interaction terms Φi(p1,m1;p2,m2,q12)subscriptΦ𝑖subscript𝑝1subscript𝑚1subscript𝑝2subscript𝑚2subscript𝑞12\Phi_{i}(p_{1},m_{1};p_{2},m_{2},q_{12}...)roman_Φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ; italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , italic_q start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT … )

pi2mi2=Φi(p1,m1;p2,m2,q12,),subscriptsuperscript𝑝2𝑖subscriptsuperscript𝑚2𝑖subscriptΦ𝑖subscript𝑝1subscript𝑚1subscript𝑝2subscript𝑚2subscript𝑞12\displaystyle p^{2}_{i}-m^{2}_{i}=\Phi_{i}(p_{1},m_{1};p_{2},m_{2},q_{12},...),italic_p start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT - italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = roman_Φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ; italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , italic_q start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT , … ) , (1)

where i=1,2𝑖12i=1,2italic_i = 1 , 2 numerates particles, whereas the functions ΦisubscriptΦ𝑖\Phi_{i}roman_Φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT depends on the particle momenta pisubscript𝑝𝑖p_{i}italic_p start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT and on the parameters q12subscript𝑞12q_{12}italic_q start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT characterizing interaction between particles.

On the other hand, in non-relativistic physics the III-rd Newton law restricts a set of the interaction and external fields which allow to reduce the two-body problem to the single particle dynamics. The similar restrictions dictated by the relativistic version of the III-rd Newton law to external and interaction fields should obviously take place in the relativistic case. The relativistic version of this law is given by a formulaCra84 ; Cra87 ; Cra90 ; Cra01

Pp=0,𝑃𝑝0\displaystyle Pp=0,italic_P italic_p = 0 , (2)

whereP𝑃Pitalic_P and p𝑝pitalic_p are the total and relative momenta of particles. Such a restriction turns out to b the necessary condition of the compatibility of equation in studying of the two-body problem in the constraint dynamics approachDir ; Ber .

Based on the constraint dynamics ideasDir ; Ber the classical relativistic two-particle system has been studied Droz ; Kal ; Teo78 ; Kom for the first time, where a consistent and covariant approach has been developed for the relativistic canonical Hamiltonian mechanics. For more than forty years, beginning from 80th, the constraint dynamic approach have been actively and successfully developed in studying the quantum two-body systems, and applied for describing real compound systems in QED and QCD Sad86-1 ; Sad86-2 ; Sad89 ; Cra84 ; Cra87 ; Cra90 ; Cra01 ; Cra03 ; Cra10 . Following the constraint dynamics the two interacting particle systems, involving spin-0 bosons and/or spin-(1/2) fermions, is studied Sad86-1 ; Sad86-2 . The obtained results and is applied to the description of the quarkonium states Sad86-2 . In the framework of the constraint dynamics approach the 4×4444\times 44 × 4 matrix wave function satisfying an equation of the Pauli-Schrodinger type is constructedSad94 . Later, the developed approach was generalized to the case of the N interaction particleSad89 .

Along with researchesSad86-1 ; Sad86-2 ; Sad89 , by developing previous resultsTeo76 ; Teo78 , the quantum consideration of two spinlessCra84 and spin-1/2121/21 / 2Cra84 particles under vector and scalar interaction, have been developed, using the constraint mechanics for the two-body Klein-Gordon and Dirac equations. The hyperbolic interaction structure of these equations, established earlierCra84 ; Cra87 , allowsCra90 us to include the pseudoscalar, pseudovector, and tensor interactions in the co responding motion equations. The Hamiltonian formulation of the relativistic N-body problem in a separable two-body basis, in which the particles interact pair-wise through scalar and vector interactions, is developed in Ref.Cra01 . The worked out approachesCra84 ; Cra87 ; Cra90 ; Cra01 have been further applied to researches in a meson spectroscopy Cra06 ; Cra10 , and in nuclear-nuclear collisionsCra03 . The two-body states in the case of QED interactionBar87 and when the interaction is realized due to a neutral bosonBar91 , has been studied in terms of the action integral in the coordinate representation that leads to the non-covariant Breit-like equation in the coordinate space-time.

The dynamics of the two-body systems is studied in the Lagrange formalism in the case when such a system consists of two classical relativistic particlesDom78 ; Dom80 . The singular-Lagrangian approach to the two-body problem is derived thereinDom78 ; Dom80 .

The relativistic two-particle systems in external fields which consist of scalarDroz00 ; Droz09 and fermionBij00 particles are considered in the context of the covariant separation of the relative and mutual motion such systems. Having been mainly concentrated on the mutual motion and on constraint relations concerning external fieldDroz00 ; Bij00 ; Droz09 , the reconstruction of the relative motion spectrum as well as the spin mixing effect turn out beyond the consideration mentioned aboveDroz00 ; Bij00 ; Droz09 .

Despite the active use of the relativistic version of the third Newton lawCra84 ; Cra87 ; Cra90 ; Cra01 , the requirements for the interaction potentials that allow to say about the two-body problem in the relativistic case, as well as the modification of the third Newton law in external field potentials, have still remained to be undiscussable. The same concerns arising the special phase space for the two-body system evolution as well as relations between the seed constituent and effective phenomenological potential what are beyond previous considerationsSad86-1 ; Sad86-2 ; Sad89 ; Cra84 ; Cra87 ; Cra90 ; Cra01 ; Cra03 ; Cra10 . Besides that, even in the absence of external field, the Dirac equation derived in earlierSad86-1 ; Sad86-2 ; Sad89 ; Cra84 ; Cra87 ; Cra90 ; Cra01 ; Cra03 ; Cra10 turn out to be very complicate for analyzing the dynamics of fields in real physical situations. Therefore, the question whether it is possible to find some new relevant γ𝛾\gammaitalic_γ-matrix basis which allows to simplify the Dirac equation structure in the two-body problem, has remained unsolved.

In the present paper the two-body problem in the external vector fields is studied in the Lagrange formalism beyond the constraint Hamiltonian dynamics. The exact covariant operator equations governing scalar particles and spin-1/2121/21 / 2 fermions have been derived. The obtained equation are formulated in the standard space-time variables of the relative and center mass motion of the two-body system, beyond the consideration in the Breit frame of reference onlySad86-1 ; Sad86-2 ; Sad89 ; Cra84 ; Cra87 ; Cra90 ; Cra01 ; Cra03 ; Cra10 . The class of interaction and external field potentials when the two-body problem can be covariantly reformulated in (3+1) phase space of the relative motion variables, which are uncoupled from the center mass variables, is pointed out.

Provided that Newton’s third law occurs, the Klein-Gordon equation, which governs either dynamics of a scalar particle or describes each of a single component of the 16 components spinor of the spin-1/2121/21 / 2 fermion, depending on the particle type, has been derived in the (3+1) phase space. In the case of 1/2121/21 / 2-spin particles we find a new γ𝛾\gammaitalic_γ-matrix basis which allows us to go to the Dirac-like equation which is decoupled in four independent equations for four different components of the initial 16 component spinor. In terms of the derived Dirac-like equations the Lagrangian governing the dynamics of two spin-1/2121/21 / 2 fermions in an external field has been found. This Lagrangian generates the conservative fermion current.

In the case of pure QED the developed approach is applied to studying the low energy states of a positronium in a magnetic field, taking into account spin-mixing effect. It is shown that strong magnetic fields result in essentially reconstructing the spectrum of positronium states and decreasing the life-time of the positronium ground state.

The paper is organized as follows. The system consisting of two interacting scalar particles in external fields is studied in Section 2. Two interacting particles with spin-1/2121/21 / 2 in external fields are considered in Section 3. Section 4 is devoted to studying the low energy states of a positronium in a strong uniform magnetic field. The last section is the conclusion. The details of some calculations are given in an appendix.

2 Two interacting spinless particles in external fields

We consider two scalar particles with masses m1subscript𝑚1m_{1}italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and m2subscript𝑚2m_{2}italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT which interact by means of vector Aintν(x)superscriptsubscript𝐴𝑖𝑛𝑡𝜈𝑥A_{int}^{\nu}(x)italic_A start_POSTSUBSCRIPT italic_i italic_n italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT ( italic_x ) and scalar Sint(x)subscript𝑆𝑖𝑛𝑡𝑥S_{int}(x)italic_S start_POSTSUBSCRIPT italic_i italic_n italic_t end_POSTSUBSCRIPT ( italic_x ) potentials, and which also are in an external vector Aextν(x)superscriptsubscript𝐴𝑒𝑥𝑡𝜈𝑥A_{ext}^{\nu}(x)italic_A start_POSTSUBSCRIPT italic_e italic_x italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT ( italic_x ) and scalar Sext(x)subscript𝑆𝑒𝑥𝑡𝑥S_{ext}(x)italic_S start_POSTSUBSCRIPT italic_e italic_x italic_t end_POSTSUBSCRIPT ( italic_x ) potentials. Such external potentials transform the standard confined two-body problem into the open one. The action integral of such a system of particles consists of two terms

𝔄=i=12λi2𝑑xi(|πiμΦi(xi)|2Φi(xi)(mi+𝒮i(xi))2Φi(xi))𝔄superscriptsubscript𝑖12subscript𝜆𝑖2differential-dsubscript𝑥𝑖superscriptsubscriptsuperscript𝜋𝜇𝑖subscriptΦ𝑖subscript𝑥𝑖2subscriptsuperscriptΦ𝑖subscript𝑥𝑖superscriptsubscript𝑚𝑖subscript𝒮𝑖subscript𝑥𝑖2subscriptΦ𝑖subscript𝑥𝑖\displaystyle{\mathfrak{A}}=\sum\limits_{i=1}^{2}\frac{\lambda_{i}}{2}\int dx_% {i}\Big{(}|\pi^{\mu}_{i}\Phi_{i}(x_{i})|^{2}-\Phi^{\ast}_{i}(x_{i})(m_{i}+{% \cal S}_{i}(x_{i}))^{2}\Phi_{i}(x_{i})\Big{)}fraktur_A = ∑ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT divide start_ARG italic_λ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG ∫ italic_d italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( | italic_π start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT roman_Φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - roman_Φ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) ( italic_m start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT + caligraphic_S start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_Φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) ) (3b)
πiνpiν+𝒜iν(xi)=iiν+𝒜iν(xi),subscriptsuperscript𝜋𝜈𝑖subscriptsuperscript𝑝𝜈𝑖superscriptsubscript𝒜𝑖𝜈subscript𝑥𝑖𝑖subscriptsuperscript𝜈𝑖superscriptsubscript𝒜𝑖𝜈subscript𝑥𝑖\displaystyle\pi^{\nu}_{i}\equiv p^{\nu}_{i}+{\cal A}_{i}^{\nu}(x_{i})=i% \partial^{\nu}_{i}+{\cal A}_{i}^{\nu}(x_{i}),italic_π start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ≡ italic_p start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT + caligraphic_A start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) = italic_i ∂ start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT + caligraphic_A start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) ,

where i=1,2𝑖12i=1,2italic_i = 1 , 2 enumerates particles, Φi(xi)subscriptΦ𝑖subscript𝑥𝑖\Phi_{i}(x_{i})roman_Φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) denotes a scalar field, , piν=iiν=(i0,i)subscriptsuperscript𝑝𝜈𝑖𝑖subscriptsuperscript𝜈𝑖𝑖superscript0𝑖p^{\nu}_{i}=i\partial^{\nu}_{i}=(i\partial^{0},\;-i\nabla)italic_p start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = italic_i ∂ start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = ( italic_i ∂ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT , - italic_i ∇ ) is a particle momentum; gνμsubscript𝑔𝜈𝜇g_{\nu\mu}italic_g start_POSTSUBSCRIPT italic_ν italic_μ end_POSTSUBSCRIPT is the metric tensor which diagonal is diag=(1,1,1,1)𝑑𝑖𝑎𝑔1111diag=(1,-1,-1,-1)italic_d italic_i italic_a italic_g = ( 1 , - 1 , - 1 , - 1 ), whereas x=(x0,𝒙)𝑥superscript𝑥0𝒙x=(x^{0},\hbox{\boldmath${x}$})italic_x = ( italic_x start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT , bold_italic_x ). The parameters λisubscript𝜆𝑖\lambda_{i}italic_λ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT are the Lagrange multipliers; 𝒜iν(xi)superscriptsubscript𝒜𝑖𝜈subscript𝑥𝑖{\cal A}_{i}^{\nu}(x_{i})caligraphic_A start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) and 𝒮i(xi)subscript𝒮𝑖subscript𝑥𝑖{\cal S}_{i}(x_{i})caligraphic_S start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) are the vector and scalar potentials acting on the field Φi(xi)subscriptΦ𝑖subscript𝑥𝑖\Phi_{i}(x_{i})roman_Φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ), respectively. These potentials consist of two terms which describe interaction between particles inside the two-body and an external fields acting on each particle independently. For brevity, the coupling constants are not written explicitly, and are included into the potential 𝒜iν(xi)superscriptsubscript𝒜𝑖𝜈subscript𝑥𝑖{\cal A}_{i}^{\nu}(x_{i})caligraphic_A start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) and 𝒮i(xi)subscript𝒮𝑖subscript𝑥𝑖{\cal S}_{i}(x_{i})caligraphic_S start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ).

Further consideration requires to discuss the structure of the interaction terms in Eq.(3b). The fields governing interaction inside the confined two-body problem have to obviously take into account the states of both interaction particle, whereas any external field acting on the each of the considered particles is the same, differing by the coupling constant in the interaction term. Then, we write the interaction potentials as follows

𝒜iν(xi)=Aiν(xi)+𝑑xjΦj(xj)Aijν(xixj)Φj(xj)𝑑xjΦj(xj)Φj(xj),subscriptsuperscript𝒜𝜈𝑖subscript𝑥𝑖subscriptsuperscript𝐴𝜈𝑖subscript𝑥𝑖differential-dsubscript𝑥𝑗subscriptsuperscriptΦ𝑗subscript𝑥𝑗subscriptsuperscript𝐴𝜈𝑖𝑗subscript𝑥𝑖subscript𝑥𝑗subscriptΦ𝑗subscript𝑥𝑗differential-dsubscript𝑥𝑗subscriptsuperscriptΦ𝑗subscript𝑥𝑗subscriptΦ𝑗subscript𝑥𝑗\displaystyle{\cal A}^{\nu}_{i}(x_{i})={A}^{\nu}_{i}(x_{i})+\frac{\int dx_{j}% \Phi^{\ast}_{j}(x_{j}){A}^{\nu}_{ij}(x_{i}-x_{j})\Phi_{j}(x_{j})}{\int dx_{j}% \Phi^{\ast}_{j}(x_{j})\Phi_{j}(x_{j})},caligraphic_A start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) = italic_A start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) + divide start_ARG ∫ italic_d italic_x start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT roman_Φ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) italic_A start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT - italic_x start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) roman_Φ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) end_ARG start_ARG ∫ italic_d italic_x start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT roman_Φ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) roman_Φ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) end_ARG ,
(4a)
𝒮i(xi)=Si(xi)+𝑑xjΦj(xj)Sij(xixj)Φj(x)𝑑xjΦj(xj)Φj(xj).subscript𝒮𝑖subscript𝑥𝑖subscript𝑆𝑖subscript𝑥𝑖differential-dsubscript𝑥𝑗subscriptsuperscriptΦ𝑗subscript𝑥𝑗subscript𝑆𝑖𝑗subscript𝑥𝑖subscript𝑥𝑗subscriptΦ𝑗𝑥differential-dsubscript𝑥𝑗subscriptsuperscriptΦ𝑗subscript𝑥𝑗subscriptΦ𝑗subscript𝑥𝑗\displaystyle{\cal S}_{i}(x_{i})={S}_{i}(x_{i})+\frac{\int dx_{j}\Phi^{\ast}_{% j}(x_{j}){S}_{ij}(x_{i}-x_{j})\Phi_{j}(x)}{\int dx_{j}\Phi^{\ast}_{j}(x_{j})% \Phi_{j}(x_{j})}.caligraphic_S start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) = italic_S start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) + divide start_ARG ∫ italic_d italic_x start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT roman_Φ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) italic_S start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT - italic_x start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) roman_Φ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ( italic_x ) end_ARG start_ARG ∫ italic_d italic_x start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT roman_Φ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) roman_Φ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) end_ARG . (4b)

Here Aiνsubscriptsuperscript𝐴𝜈𝑖{A}^{\nu}_{i}italic_A start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT and Sjsubscript𝑆𝑗S_{j}italic_S start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT denote the potentials affecting the i𝑖iitalic_i-th particle by an external field, whereas Aijνsubscriptsuperscript𝐴𝜈𝑖𝑗{A}^{\nu}_{ij}italic_A start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT and Sijsubscript𝑆𝑖𝑗S_{ij}italic_S start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT supplied by a double subscribe notate the potentials acting by the j𝑗jitalic_j-th particle on the i𝑖iitalic_i-the one, inside the two-body system ( i,j=1,2,ijformulae-sequence𝑖𝑗12𝑖𝑗i,j=1,2,i\neq jitalic_i , italic_j = 1 , 2 , italic_i ≠ italic_j). All of Aijνsubscriptsuperscript𝐴𝜈𝑖𝑗{A}^{\nu}_{ij}italic_A start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT depends , in general, on the fields Φj(xj)subscriptΦ𝑗subscript𝑥𝑗\Phi_{j}(x_{j})roman_Φ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ). However, we ignore such a dependence, restricting by consideration of the so-called mean field approximation for the two-particle interaction.

As for the external potentials, to eliminate interaction between external and interaction potential we demand them orthogonality which means the following in the case of the vector potentials

Aijν(Aν)k=0,i,j,k=1,2.formulae-sequencesubscriptsuperscript𝐴𝜈𝑖𝑗subscriptsubscript𝐴𝜈𝑘0𝑖𝑗𝑘12\displaystyle{A}^{\nu}_{ij}(A_{\nu})_{k}=0,~{}~{}~{}i,j,k=1,2.italic_A start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT ( italic_A start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT = 0 , italic_i , italic_j , italic_k = 1 , 2 . (5)

The restrictions to the external scalar potentials will be discussed below in Sec.2.4.

2.1 Free particles

Let us consider the two-body problem in case of free particle to clarify selection of the Lagrange multiplier λ1,2subscript𝜆12\lambda_{1,2}italic_λ start_POSTSUBSCRIPT 1 , 2 end_POSTSUBSCRIPT.

We firstly take the Lagrange multipliers to be

λ1=(𝑑x2Φ2(x2)Φ2(x2)),λ2=(𝑑x1Φ1(x1)Φ1(x1)).formulae-sequencesubscript𝜆1differential-dsubscript𝑥2subscriptsuperscriptΦ2subscript𝑥2subscriptΦ2subscript𝑥2subscript𝜆2differential-dsubscript𝑥1subscriptsuperscriptΦ1subscript𝑥1subscriptΦ1subscript𝑥1\displaystyle\lambda_{1}=\Big{(}\int dx_{2}\Phi^{\ast}_{2}(x_{2})\Phi_{2}(x_{2% })\Big{)},~{}~{}~{}\lambda_{2}=\Big{(}\int dx_{1}\Phi^{\ast}_{1}(x_{1})\Phi_{1% }(x_{1})\Big{)}.italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = ( ∫ italic_d italic_x start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT roman_Φ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) roman_Φ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) ) , italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = ( ∫ italic_d italic_x start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT roman_Φ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) roman_Φ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) ) .

Then, we introduce the canonically conjugation variables111The same momenta are used earlierTeo71 ; Teo76

P=p1+p2,p=μ1p2μ2p1,formulae-sequence𝑃subscript𝑝1subscript𝑝2𝑝subscript𝜇1subscript𝑝2subscript𝜇2subscript𝑝1\displaystyle~{}P=p_{1}+p_{2},~{}~{}~{}p=\mu_{1}p_{2}-\mu_{2}p_{1},italic_P = italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , italic_p = italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , (7a)
p1=μ1Pp,p2=μ2P+p,formulae-sequencesubscript𝑝1subscript𝜇1𝑃𝑝subscript𝑝2subscript𝜇2𝑃𝑝\displaystyle p_{1}=\mu_{1}P-p,~{}~{}~{}p_{2}=\mu_{2}P+p,italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_P - italic_p , italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_P + italic_p , (7b)
X=12(x1+x2),x=x2x1.formulae-sequence𝑋12subscript𝑥1subscript𝑥2𝑥subscript𝑥2subscript𝑥1\displaystyle X=\frac{1}{2}(x_{1}+x_{2}),~{}~{}~{}~{}x=x_{2}-x_{1}.italic_X = divide start_ARG 1 end_ARG start_ARG 2 end_ARG ( italic_x start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_x start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) , italic_x = italic_x start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - italic_x start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT . (7c)

where

μ1=12(1+m12m22M2),μ2=12(1m12m22M2),formulae-sequencesubscript𝜇1121superscriptsubscript𝑚12superscriptsubscript𝑚22superscript𝑀2subscript𝜇2121superscriptsubscript𝑚12superscriptsubscript𝑚22superscript𝑀2\displaystyle\mu_{1}=\frac{1}{2}\left(1+\frac{m_{1}^{2}-m_{2}^{2}}{M^{2}}% \right),\mu_{2}=\frac{1}{2}\left(1-\frac{m_{1}^{2}-m_{2}^{2}}{M^{2}}\right),italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG ( 1 + divide start_ARG italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_M start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) , italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG ( 1 - divide start_ARG italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_M start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) ,
(8a)
μ1+μ2=1,M=|p1+p2|=|P|.formulae-sequencesubscript𝜇1subscript𝜇21𝑀subscript𝑝1subscript𝑝2𝑃\displaystyle\mu_{1}+\mu_{2}=1,~{}~{}~{}~{}M=|p_{1}+p_{2}|=|P|.italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = 1 , italic_M = | italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT | = | italic_P | . (8b)

In Eq. (7a) P𝑃Pitalic_P and p𝑝pitalic_p are the total and relative momentum of the two-body system, respectively, while M=|P|=|p1+p2|𝑀𝑃subscript𝑝1subscript𝑝2M=|P|=|p_{1}+p_{2}|italic_M = | italic_P | = | italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT | is its invariant mass. The momenta P,p,𝑃𝑝P,p,italic_P , italic_p , and coordinates X,x𝑋𝑥X,xitalic_X , italic_x satisfy the following commutative relations

[Pj,xk]=[pj,Xk]=0,subscript𝑃𝑗subscript𝑥𝑘subscript𝑝𝑗subscript𝑋𝑘0\displaystyle[P_{j},x_{k}]=[p_{j},X_{k}]=0,[ italic_P start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT , italic_x start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ] = [ italic_p start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT , italic_X start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ] = 0 ,
[Pj,Xk]=iδjk,[pj,xk]=iδjk,j,k=1,2,3.formulae-sequencesubscript𝑃𝑗subscript𝑋𝑘𝑖subscript𝛿𝑗𝑘formulae-sequencesubscript𝑝𝑗subscript𝑥𝑘𝑖subscript𝛿𝑗𝑘𝑗𝑘123\displaystyle[P_{j},X_{k}]=-i\delta_{jk},[p_{j},x_{k}]=-i\delta_{jk},~{}j,k=1,% 2,3.[ italic_P start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT , italic_X start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ] = - italic_i italic_δ start_POSTSUBSCRIPT italic_j italic_k end_POSTSUBSCRIPT , [ italic_p start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT , italic_x start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ] = - italic_i italic_δ start_POSTSUBSCRIPT italic_j italic_k end_POSTSUBSCRIPT , italic_j , italic_k = 1 , 2 , 3 . (9)

In terms of such new variables the action integral becomes to be equal to

𝔄+=𝑑X𝑑x𝑑X𝑑xΦ(X,x)(p2+Ew2mw2m12m22M2Pp)Φ(X,x),subscript𝔄differential-d𝑋differential-d𝑥differential-d𝑋differential-d𝑥superscriptΦ𝑋𝑥superscript𝑝2superscriptsubscript𝐸𝑤2superscriptsubscript𝑚𝑤2superscriptsubscript𝑚12superscriptsubscript𝑚22superscript𝑀2𝑃𝑝Φ𝑋𝑥\displaystyle{\mathfrak{A}}_{+}=\int dXdx\int dXdx\Phi^{\ast}(X,x)\big{(}p^{2}% +E_{w}^{2}-m_{w}^{2}-\frac{m_{1}^{2}-m_{2}^{2}}{M^{2}}Pp\big{)}\Phi(X,x),fraktur_A start_POSTSUBSCRIPT + end_POSTSUBSCRIPT = ∫ italic_d italic_X italic_d italic_x ∫ italic_d italic_X italic_d italic_x roman_Φ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( italic_X , italic_x ) ( italic_p start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_m start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - divide start_ARG italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_M start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_P italic_p ) roman_Φ ( italic_X , italic_x ) , (10)

where we introduce an effective energy and an effective mass of the two-body systemTeo76 ; Teo78 ; Cra84 ; Cra87 ,

Ew=12M(M2(m12+m22)),mw=m1m2M,formulae-sequencesubscript𝐸𝑤12𝑀superscript𝑀2superscriptsubscript𝑚12superscriptsubscript𝑚22subscript𝑚𝑤subscript𝑚1subscript𝑚2𝑀\displaystyle E_{w}=\frac{1}{2M}(M^{2}-(m_{1}^{2}+m_{2}^{2})),~{}~{}~{}m_{w}=% \frac{m_{1}m_{2}}{M},italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 italic_M end_ARG ( italic_M start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - ( italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) ) , italic_m start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT = divide start_ARG italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG italic_M end_ARG , (11)

and the two-body field Φ(X,x)Φ𝑋𝑥\Phi(X,x)roman_Φ ( italic_X , italic_x )

Φ(X,x)=Φ1(Xx/2)Φ2(X+x/2).Φ𝑋𝑥subscriptΦ1𝑋𝑥2subscriptΦ2𝑋𝑥2\displaystyle\Phi(X,x)=\Phi_{1}(X-x/2)\Phi_{2}(X+x/2).roman_Φ ( italic_X , italic_x ) = roman_Φ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_X - italic_x / 2 ) roman_Φ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_X + italic_x / 2 ) . (12)

Varying with respect to the field Φ(X,x)superscriptΦ𝑋𝑥\Phi^{\ast}(X,x)roman_Φ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( italic_X , italic_x ) we obtain the motion equation

(Ew2+p2mw2(μ1μ2)Pp)Φ(X,x)=0.superscriptsubscript𝐸𝑤2superscript𝑝2superscriptsubscript𝑚𝑤2subscript𝜇1subscript𝜇2𝑃𝑝Φ𝑋𝑥0\displaystyle\bigg{(}{E}_{w}^{2}+p^{2}-{m}_{w}^{2}-\big{(}{{\mu_{1}}}-{{\mu_{2% }}}\big{)}{P}{p}\bigg{)}\Phi(X,x)=0.( italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_p start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_m start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - ( italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) italic_P italic_p ) roman_Φ ( italic_X , italic_x ) = 0 . (13)

Let us take the Lagrange multiplier as follows

λ1=(𝑑xΦ2(x)Φ2(x)),λ2=(𝑑xΦ1(x)Φ1(x)).formulae-sequencesubscript𝜆1differential-d𝑥subscriptsuperscriptΦ2𝑥subscriptΦ2𝑥subscript𝜆2differential-d𝑥subscriptsuperscriptΦ1𝑥subscriptΦ1𝑥\displaystyle\lambda_{1}=\Big{(}\int dx\Phi^{\ast}_{2}(x)\Phi_{2}(x)\Big{)},% \lambda_{2}=-\Big{(}\int dx\Phi^{\ast}_{1}(x)\Phi_{1}(x)\Big{)}.italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = ( ∫ italic_d italic_x roman_Φ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_x ) roman_Φ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_x ) ) , italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = - ( ∫ italic_d italic_x roman_Φ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) roman_Φ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) ) .

Then, we get the action integral 𝔄subscript𝔄{\mathfrak{A}}_{-}fraktur_A start_POSTSUBSCRIPT - end_POSTSUBSCRIPT

𝔄=𝑑X𝑑x𝑑X𝑑xΦ(X,x)(Pp)Φ(X,x).subscript𝔄differential-d𝑋differential-d𝑥differential-d𝑋differential-d𝑥superscriptΦ𝑋𝑥𝑃𝑝Φ𝑋𝑥\displaystyle{\mathfrak{A}}_{-}=\int dXdx\int dXdx\Phi^{\ast}(X,x)\big{(}Pp% \big{)}\Phi(X,x).fraktur_A start_POSTSUBSCRIPT - end_POSTSUBSCRIPT = ∫ italic_d italic_X italic_d italic_x ∫ italic_d italic_X italic_d italic_x roman_Φ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( italic_X , italic_x ) ( italic_P italic_p ) roman_Φ ( italic_X , italic_x ) . (15)

Such an action integral generates the motion equation

(Pp)Φ(X,x)=0.𝑃𝑝Φ𝑋𝑥0\displaystyle({P}{p})~{}\Phi(X,x)=0.( italic_P italic_p ) roman_Φ ( italic_X , italic_x ) = 0 . (16)

This motion equation corresponds to the relativistic version of the third Newton’s lawCra84 ; Cra87 and is in entire relation to Eq.(19)Cra87 . Such a results is physically understandable since the Lagrange multipliers selection dictated by Eq.(2.1) keeps the particle interaction only, eliminating their kinematics.

Provided that Eq.(16) the motion equation takes a form

(Ew2+p2mw2)Φ(X,x)=0.superscriptsubscript𝐸𝑤2superscript𝑝2superscriptsubscript𝑚𝑤2Φ𝑋𝑥0\displaystyle\bigg{(}{E}_{w}^{2}+p^{2}-{m}_{w}^{2}\bigg{)}\Phi(X,x)=0.( italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_p start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_m start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) roman_Φ ( italic_X , italic_x ) = 0 . (17)

This equation is covariant in 13(X)13(x)tensor-productsubscriptsuperscript31𝑋subscriptsuperscript31𝑥{\mathfrak{R}}^{3}_{1}(X)\otimes{\mathfrak{R}}^{3}_{1}(x)fraktur_R start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_X ) ⊗ fraktur_R start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) space corresponding to the both center-mass and relative motion. Such types of motion can be obviously decoupled looking for the solution of Eq.(17) as

Φ(X,x)=χ(X)φ(x),Φ𝑋𝑥𝜒𝑋𝜑𝑥\displaystyle\Phi(X,x)=\chi(X)\varphi(x),roman_Φ ( italic_X , italic_x ) = italic_χ ( italic_X ) italic_φ ( italic_x ) , (18)

where χ(X)𝜒𝑋\chi(X)italic_χ ( italic_X ) is the eigenfuncion of an equation

Ew2χ(X)=Ew2(P)χ(X).superscriptsubscript𝐸𝑤2𝜒𝑋superscriptsubscript𝐸𝑤2𝑃𝜒𝑋\displaystyle{E}_{w}^{2}\chi(X)={E}_{w}^{2}(P)\chi(X).italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_χ ( italic_X ) = italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_P ) italic_χ ( italic_X ) . (19)

As a result, we go to the Eq.(20) for the function ϕ(x)italic-ϕ𝑥\phi(x)italic_ϕ ( italic_x ) where Ew2(P)superscriptsubscript𝐸𝑤2𝑃{E}_{w}^{2}(P)italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_P ) has already been a C𝐶Citalic_C-function.

(Ew2(P)+p2mw2)φ(x)=0.superscriptsubscript𝐸𝑤2𝑃superscript𝑝2superscriptsubscript𝑚𝑤2𝜑𝑥0\displaystyle\bigg{(}{E}_{w}^{2}(P)+p^{2}-{m}_{w}^{2}\bigg{)}\varphi(x)=0.( italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_P ) + italic_p start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_m start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) italic_φ ( italic_x ) = 0 . (20)

Although this equation is covariant in the 13(x)subscriptsuperscript31𝑥{\mathfrak{R}}^{3}_{1}(x)fraktur_R start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) but it cannot be reduced to the standard Klein-Gordon equation.

We look for the solution of Eq.(20) as follows

φ(x)=φ(𝒙)exp(iτp0).𝜑𝑥𝜑𝒙𝑖𝜏superscript𝑝0\displaystyle\varphi(x)=\varphi(\hbox{\boldmath${x}$})\exp(-i\tau p^{0}).italic_φ ( italic_x ) = italic_φ ( bold_italic_x ) roman_exp ( - italic_i italic_τ italic_p start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT ) . (21)

We breaking the 13(x)subscriptsuperscript31𝑥{\mathfrak{R}}^{3}_{1}(x)fraktur_R start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) covariance in favor of keeping the Klein-Gordon form of the motion equation, that turns out to be reasonable under studying two interaction particles, and introduce εwsubscript𝜀𝑤\varepsilon_{w}italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT and θwsubscript𝜃𝑤\theta_{w}italic_θ start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT

Ew(P)=εwcosθwsubscript𝐸𝑤𝑃subscript𝜀𝑤subscript𝜃𝑤\displaystyle{E}_{w}(P)=\varepsilon_{w}\cos\theta_{w}italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT ( italic_P ) = italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT roman_cos italic_θ start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT (22a)
p0=εwsinθw.superscript𝑝0subscript𝜀𝑤subscript𝜃𝑤\displaystyle p^{0}=\varepsilon_{w}\sin\theta_{w}.italic_p start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT = italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT roman_sin italic_θ start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT . (22b)

As a result, we get

(εw2+𝒑2mw2)φ(𝒙)=0,superscriptsubscript𝜀𝑤2superscript𝒑2superscriptsubscript𝑚𝑤2𝜑𝒙0\displaystyle\bigg{(}{\varepsilon}_{w}^{2}+\hbox{\boldmath${p}$}^{2}-{m}_{w}^{% 2}\bigg{)}\varphi(\hbox{\boldmath${x}$})=0,( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + bold_italic_p start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_m start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) italic_φ ( bold_italic_x ) = 0 , (23)

the obtained equation is the Klein-Gordon equation which is covariant in the 13(εw,𝒙)subscriptsuperscript31subscript𝜀𝑤𝒙{\mathfrak{R}}^{3}_{1}(\varepsilon_{w},\hbox{\boldmath${x}$})fraktur_R start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_x ) phase spaceTeo71 ; Teo76 ; Droz ; Kal ; Teo78 ; Kom ; Are ; Kir ; Dom78 ; Dom80 ; Sad86-1 ; Sad86-2 ; Sad94 ; Sad89 ; Cra84 ; Cra87 ; Cra90 ; Cra01 ; Cra03 ; Cra06 ; Cra10 rather than the standard coordinate 13(x)subscriptsuperscript31𝑥{\mathfrak{R}}^{3}_{1}(x)fraktur_R start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) space-time.

2.2 Interacting particles

In the case of interacting particles we construct the action integrals 𝔄+subscript𝔄{\mathfrak{A}}_{+}fraktur_A start_POSTSUBSCRIPT + end_POSTSUBSCRIPT and 𝔄subscript𝔄{\mathfrak{A}}_{-}fraktur_A start_POSTSUBSCRIPT - end_POSTSUBSCRIPT by taking the Lagrange multipliers λ1,2subscript𝜆12\lambda_{1,2}italic_λ start_POSTSUBSCRIPT 1 , 2 end_POSTSUBSCRIPT according to Eqs.(2.1), (2.1) as it has been done above for the free particles. Such action integrals generate the motion equations after going to the center mass and relative motion variables (see Eq.(7)).

{2(p2+Ew2mw2(μ1μ2)Pp)+𝒜12(Xx/2)+2𝒜1(Xx/2)(μ2Pp)+\displaystyle\bigg{\{}2(p^{2}+E_{w}^{2}-m_{w}^{2}-(\mu_{1}-\mu_{2})Pp)+{\cal A% }_{1}^{2}(X-x/2)+2{\cal A}_{1}(X-x/2)(\mu_{2}P-p)+{ 2 ( italic_p start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_m start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - ( italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) italic_P italic_p ) + caligraphic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_X - italic_x / 2 ) + 2 caligraphic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_X - italic_x / 2 ) ( italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_P - italic_p ) +
[(μ2Pp);𝒜1(Xx/2)]+𝒜22(X+x/2)+2𝒜2(X+x/2)(μ1P+p)+subscript𝜇2𝑃𝑝subscript𝒜1𝑋𝑥2superscriptsubscript𝒜22𝑋𝑥2limit-from2subscript𝒜2𝑋𝑥2subscript𝜇1𝑃𝑝\displaystyle[(\mu_{2}P-p);{\cal A}_{1}(X-x/2)]+{\cal A}_{2}^{2}(X+x/2)+2{\cal A% }_{2}(X+x/2)(\mu_{1}P+p)+[ ( italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_P - italic_p ) ; caligraphic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_X - italic_x / 2 ) ] + caligraphic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_X + italic_x / 2 ) + 2 caligraphic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_X + italic_x / 2 ) ( italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_P + italic_p ) +
[(μ1P+p);𝒜2(X+x/2)]2m1𝒮1(Xx/2)𝒮12(Xx/2)+subscript𝜇1𝑃𝑝subscript𝒜2𝑋𝑥22subscript𝑚1subscript𝒮1𝑋𝑥2limit-fromsubscriptsuperscript𝒮21𝑋𝑥2\displaystyle\ [(\mu_{1}P+p);{\cal A}_{2}(X+x/2)]-2m_{1}{\cal S}_{1}(X-x/2)-{% \cal S}^{2}_{1}(X-x/2)+[ ( italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_P + italic_p ) ; caligraphic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_X + italic_x / 2 ) ] - 2 italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT caligraphic_S start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_X - italic_x / 2 ) - caligraphic_S start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_X - italic_x / 2 ) +
2m2𝒮2(X+x/2)𝒮22(X+x/2)}Φ(X,x)=0,\displaystyle-2m_{2}{\cal S}_{2}(X+x/2)-{\cal S}^{2}_{2}(X+x/2)\bigg{\}}\Phi(X% ,x)=0,- 2 italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT caligraphic_S start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_X + italic_x / 2 ) - caligraphic_S start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_X + italic_x / 2 ) } roman_Φ ( italic_X , italic_x ) = 0 , (24)
2PpΦ(X,x)={𝒜12(Xx/2)+2𝒜1(Xx/2)(μ2Pp)+[(μ2Pp);𝒜1(Xx/2)]\displaystyle-2Pp~{}\Phi(X,x)=\bigg{\{}{\cal A}_{1}^{2}(X-x/2)+2{\cal A}_{1}(X% -x/2)(\mu_{2}P-p)+[(\mu_{2}P-p);{\cal A}_{1}(X-x/2)]-- 2 italic_P italic_p roman_Φ ( italic_X , italic_x ) = { caligraphic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_X - italic_x / 2 ) + 2 caligraphic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_X - italic_x / 2 ) ( italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_P - italic_p ) + [ ( italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_P - italic_p ) ; caligraphic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_X - italic_x / 2 ) ] -
𝒜22(X+x/2)2𝒜2(X+x/2)(μ1P+p)+[(μ1P+p);𝒜2(X+x/2)]superscriptsubscript𝒜22𝑋𝑥22subscript𝒜2𝑋𝑥2subscript𝜇1𝑃𝑝limit-fromsubscript𝜇1𝑃𝑝subscript𝒜2𝑋𝑥2\displaystyle{\cal A}_{2}^{2}(X+x/2)-2{\cal A}_{2}(X+x/2)(\mu_{1}P+p)+[(\mu_{1% }P+p);{\cal A}_{2}(X+x/2)]-caligraphic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_X + italic_x / 2 ) - 2 caligraphic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_X + italic_x / 2 ) ( italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_P + italic_p ) + [ ( italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_P + italic_p ) ; caligraphic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_X + italic_x / 2 ) ] -
2m1𝒮1(Xx/2)𝒮12(Xx/2)2m2𝒮2(X+x/2)𝒮22(X+x/2)}Φ(X,x),\displaystyle 2m_{1}{\cal S}_{1}(X-x/2)-{\cal S}^{2}_{1}(X-x/2)-2m_{2}{\cal S}% _{2}(X+x/2)-{\cal S}^{2}_{2}(X+x/2)\bigg{\}}\Phi(X,x),2 italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT caligraphic_S start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_X - italic_x / 2 ) - caligraphic_S start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_X - italic_x / 2 ) - 2 italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT caligraphic_S start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_X + italic_x / 2 ) - caligraphic_S start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_X + italic_x / 2 ) } roman_Φ ( italic_X , italic_x ) , (25)

where the square brackets denote, as usual, a commutator.

Obtained above Eq.(2.2) governs breaking the 3-rd Newton law which naturally takes place when there is no restriction to interaction potentials.

Substitution Eq.(2.2) into Eq.(2.2) leads to

{(Ew2+p2mw2)+𝒜12(Xx/2)μ1+2𝒜1(Xx/2)(μ2Pp)μ1+\displaystyle\bigg{\{}(E_{w}^{2}+p^{2}-m_{w}^{2})+{\cal A}_{1}^{2}(X-x/2)\mu_{% 1}+2{\cal A}_{1}(X-x/2)(\mu_{2}P-p)\mu_{1}+{ ( italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_p start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_m start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) + caligraphic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_X - italic_x / 2 ) italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + 2 caligraphic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_X - italic_x / 2 ) ( italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_P - italic_p ) italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT +
[(μ2Pp);𝒜1(Xx/2)]μ1+𝒜22(X+x/2)μ2+2𝒜2(X+x/2)(μ1P+p)μ2+subscript𝜇2𝑃𝑝subscript𝒜1𝑋𝑥2subscript𝜇1superscriptsubscript𝒜22𝑋𝑥2subscript𝜇2limit-from2subscript𝒜2𝑋𝑥2subscript𝜇1𝑃𝑝subscript𝜇2\displaystyle[(\mu_{2}P-p);{\cal A}_{1}(X-x/2)]\mu_{1}+{\cal A}_{2}^{2}(X+x/2)% \mu_{2}+2{\cal A}_{2}(X+x/2)(\mu_{1}P+p)\mu_{2}+[ ( italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_P - italic_p ) ; caligraphic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_X - italic_x / 2 ) ] italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + caligraphic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_X + italic_x / 2 ) italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + 2 caligraphic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_X + italic_x / 2 ) ( italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_P + italic_p ) italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT +
[(μ1P+p);𝒜2(X+x/2)]μ22m1𝒮1(Xx/2)μ1𝒮12(Xx/2)μ1+subscript𝜇1𝑃𝑝subscript𝒜2𝑋𝑥2subscript𝜇22subscript𝑚1subscript𝒮1𝑋𝑥2subscript𝜇1limit-fromsubscriptsuperscript𝒮21𝑋𝑥2subscript𝜇1\displaystyle[(\mu_{1}P+p);{\cal A}_{2}(X+x/2)]\mu_{2}-2m_{1}{\cal S}_{1}(X-x/% 2)\mu_{1}-{\cal S}^{2}_{1}(X-x/2)\mu_{1}+[ ( italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_P + italic_p ) ; caligraphic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_X + italic_x / 2 ) ] italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - 2 italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT caligraphic_S start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_X - italic_x / 2 ) italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - caligraphic_S start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_X - italic_x / 2 ) italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT +
2m2𝒮2(X+x/2)μ2𝒮22(X+x/2)μ2}Φ(X,x)=0.\displaystyle-2m_{2}{\cal S}_{2}(X+x/2)\mu_{2}-{\cal S}^{2}_{2}(X+x/2)\mu_{2}% \bigg{\}}\Phi(X,x)=0.- 2 italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT caligraphic_S start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_X + italic_x / 2 ) italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - caligraphic_S start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_X + italic_x / 2 ) italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT } roman_Φ ( italic_X , italic_x ) = 0 . (26)

Note, that Eq.(2.2) can be also derived by summing and subtracting the motion equation for each of the particles consisting of the two-body system.

2.3 The two-body problem in the absence of external fields

The center-mass and relative motion variables turn out to be uncoupled in Eq.(2.2) due external potentials since these potentials depend on the coordinate each a particle independently (see Eqs.(4), (4)). We firstly consider the situation when there is no external potentials acting on the two-body system. We also assume that the potentials 𝒜iμ=(𝒜10,𝑨i)=𝒜iμ(𝒙)subscriptsuperscript𝒜𝜇𝑖subscriptsuperscript𝒜01subscript𝑨𝑖subscriptsuperscript𝒜𝜇𝑖𝒙{\cal A}^{\mu}_{i}=({\cal A}^{0}_{1},\hbox{\boldmath${\cal\hbox{\boldmath${A}$% }}$}_{i})={\cal A}^{\mu}_{i}(\hbox{\boldmath${x}$})caligraphic_A start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = ( caligraphic_A start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , roman_A start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) = caligraphic_A start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( bold_italic_x ) and 𝒮i=𝒮i(𝒙)subscript𝒮𝑖subscript𝒮𝑖𝒙{\cal S}_{i}={\cal S}_{i}(\hbox{\boldmath${x}$})caligraphic_S start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = caligraphic_S start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( bold_italic_x ) , where i=1,2𝑖12i=1,2italic_i = 1 , 2 numerates particles, depend only on the spatial relative coordinate 𝒙𝒙{x}bold_italic_x. Moreover, let the following symmetry relation be

𝑨1μ(𝒙)=𝑨2μ(𝒙)𝑨μ(𝒙),subscriptsuperscript𝑨𝜇1𝒙subscriptsuperscript𝑨𝜇2𝒙superscript𝑨𝜇𝒙\displaystyle{\cal\hbox{\boldmath${A}$}}^{\mu}_{1}(\hbox{\boldmath${x}$})={% \cal\hbox{\boldmath${A}$}}^{\mu}_{2}(\hbox{\boldmath${x}$})\equiv{\cal\hbox{% \boldmath${A}$}}^{\mu}(\hbox{\boldmath${x}$}),bold_italic_A start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( bold_italic_x ) = bold_italic_A start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( bold_italic_x ) ≡ bold_italic_A start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT ( bold_italic_x ) , (27a)
𝒜i0(𝒙)=𝒜i0(𝒙),𝑨(𝒙)=𝑨(𝒙)formulae-sequencesubscriptsuperscript𝒜0𝑖𝒙subscriptsuperscript𝒜0𝑖𝒙𝑨𝒙𝑨𝒙\displaystyle{\cal A}^{0}_{i}(\hbox{\boldmath${x}$})={\cal A}^{0}_{i}(-\hbox{% \boldmath${x}$}),~{}~{}~{}{\cal\hbox{\boldmath${A}$}}(\hbox{\boldmath${x}$})=-% {\cal\hbox{\boldmath${A}$}}(-\hbox{\boldmath${x}$})caligraphic_A start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( bold_italic_x ) = caligraphic_A start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( - bold_italic_x ) , bold_italic_A ( bold_italic_x ) = - bold_italic_A ( - bold_italic_x ) (27b)
𝒮i(𝒙)=𝒮i(𝒙).subscript𝒮𝑖𝒙subscript𝒮𝑖𝒙\displaystyle{\cal S}_{i}(\hbox{\boldmath${x}$})={\cal S}_{i}(-\hbox{\boldmath% ${x}$}).caligraphic_S start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( bold_italic_x ) = caligraphic_S start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( - bold_italic_x ) . (27c)

Then, we separate the center-mass and relative motion variables, following Eqs.(18), (19). As a result, in terms of εwsubscript𝜀𝑤\varepsilon_{w}italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT and θwsubscript𝜃𝑤\theta_{w}italic_θ start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT given by Eq.(22) we obtain from Eq.(2.2)

{(εw2+(𝒑+𝑨(𝒙))2mw2)+μ1(𝒜10)2(𝒙)+μ2(𝒜20)2(𝒙)+2μ1μ2P0(𝒜10(𝒙)+𝒜20(𝒙))\displaystyle\bigg{\{}(\varepsilon_{w}^{2}+(\hbox{\boldmath${p}$}+{\cal\hbox{% \boldmath${A}$}}(\hbox{\boldmath${x}$}))^{2}-m_{w}^{2})+\mu_{1}({\cal A}^{0}_{% 1})^{2}(\hbox{\boldmath${x}$})+\mu_{2}({\cal A}^{0}_{2})^{2}(\hbox{\boldmath${% x}$})+2\mu_{1}\mu_{2}P^{0}({\cal A}^{0}_{1}(\hbox{\boldmath${x}$})+{\cal A}^{0% }_{2}(\hbox{\boldmath${x}$})){ ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( bold_italic_p + bold_italic_A ( bold_italic_x ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_m start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) + italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( caligraphic_A start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( bold_italic_x ) + italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( caligraphic_A start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( bold_italic_x ) + 2 italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_P start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT ( caligraphic_A start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( bold_italic_x ) + caligraphic_A start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( bold_italic_x ) )
+2p0(μ2𝒜20(𝒙)μ1𝒜10(𝒙))μ1𝒮12μ2𝒮22(𝒙)2(m1μ1𝒮1(𝒙)+m2μ2𝒮2(𝒙))}φ(𝒙)=0.\displaystyle+2p^{0}(\mu_{2}{\cal A}^{0}_{2}(\hbox{\boldmath${x}$})-\mu_{1}{% \cal A}^{0}_{1}(\hbox{\boldmath${x}$}))-\mu_{1}{\cal S}^{2}_{1}--\mu_{2}{\cal S% }^{2}_{2}(\hbox{\boldmath${x}$})-2(m_{1}\mu_{1}{\cal S}_{1}(\hbox{\boldmath${x% }$})+m_{2}\mu_{2}{\cal S}_{2}(\hbox{\boldmath${x}$}))\bigg{\}}\varphi(\hbox{% \boldmath${x}$})=0.+ 2 italic_p start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT ( italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT caligraphic_A start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( bold_italic_x ) - italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT caligraphic_A start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( bold_italic_x ) ) - italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT caligraphic_S start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - - italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT caligraphic_S start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( bold_italic_x ) - 2 ( italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT caligraphic_S start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( bold_italic_x ) + italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT caligraphic_S start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( bold_italic_x ) ) } italic_φ ( bold_italic_x ) = 0 .

We write the potentials 𝒜i0(𝒙)subscriptsuperscript𝒜0𝑖𝒙{\cal A}^{0}_{i}(\hbox{\boldmath${x}$})caligraphic_A start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( bold_italic_x ) and 𝒮1(𝒙)subscript𝒮1𝒙{\cal S}_{1}(\hbox{\boldmath${x}$})caligraphic_S start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( bold_italic_x ) as follows

𝒜10(𝒙)=1μ1𝒜0(𝒙)cosθA,𝒜20(𝒙)=1μ2𝒜0(𝒙)sinθA,formulae-sequencesubscriptsuperscript𝒜01𝒙1subscript𝜇1superscript𝒜0𝒙subscript𝜃𝐴subscriptsuperscript𝒜02𝒙1subscript𝜇2superscript𝒜0𝒙subscript𝜃𝐴\displaystyle{\cal A}^{0}_{1}(\hbox{\boldmath${x}$})=\frac{1}{\sqrt{\mu}_{1}}{% \cal A}^{0}(\hbox{\boldmath${x}$})\cos\theta_{A},~{}~{}{\cal A}^{0}_{2}(\hbox{% \boldmath${x}$})=\frac{1}{\sqrt{\mu}_{2}}{\cal A}^{0}(\hbox{\boldmath${x}$})% \sin\theta_{A},caligraphic_A start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( bold_italic_x ) = divide start_ARG 1 end_ARG start_ARG square-root start_ARG italic_μ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG caligraphic_A start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT ( bold_italic_x ) roman_cos italic_θ start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT , caligraphic_A start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( bold_italic_x ) = divide start_ARG 1 end_ARG start_ARG square-root start_ARG italic_μ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG caligraphic_A start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT ( bold_italic_x ) roman_sin italic_θ start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT , (29a)
𝒮1(𝒙)=1μ1𝒮0(𝒙)cosθS,𝒮2(𝒙)=1μ2𝒮0(𝒙)sinθS.formulae-sequencesubscript𝒮1𝒙1subscript𝜇1superscript𝒮0𝒙subscript𝜃𝑆subscript𝒮2𝒙1subscript𝜇2superscript𝒮0𝒙subscript𝜃𝑆\displaystyle{\cal S}_{1}(\hbox{\boldmath${x}$})=\frac{1}{\sqrt{\mu}_{1}}{\cal S% }^{0}(\hbox{\boldmath${x}$})\cos\theta_{S},~{}~{}{\cal S}_{2}(\hbox{\boldmath$% {x}$})=\frac{1}{\sqrt{\mu}_{2}}{\cal S}^{0}(\hbox{\boldmath${x}$})\sin\theta_{% S}.caligraphic_S start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( bold_italic_x ) = divide start_ARG 1 end_ARG start_ARG square-root start_ARG italic_μ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG caligraphic_S start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT ( bold_italic_x ) roman_cos italic_θ start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT , caligraphic_S start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( bold_italic_x ) = divide start_ARG 1 end_ARG start_ARG square-root start_ARG italic_μ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG caligraphic_S start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT ( bold_italic_x ) roman_sin italic_θ start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT . (29b)

Substituting Eq.(29) into Eq.(2.3) we derive

{εw2+(𝒑+𝑨(𝒙))2mw2+(𝒜0)2(𝒙)+2𝒜0(𝒙)P0cos(θAψA)2𝒜0(𝒙)p0cos(θA+ψA)\displaystyle\bigg{\{}\varepsilon_{w}^{2}+(\hbox{\boldmath${p}$}+{\cal\hbox{% \boldmath${A}$}}(\hbox{\boldmath${x}$}))^{2}-m_{w}^{2}+({\cal A}^{0})^{2}(% \hbox{\boldmath${x}$})+2{\cal A}^{0}(\hbox{\boldmath${x}$})P^{0}\cos(\theta_{A% }-\psi_{A})-2{\cal A}^{0}(\hbox{\boldmath${x}$})p^{0}\cos(\theta_{A}+\psi_{A}){ italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( bold_italic_p + bold_italic_A ( bold_italic_x ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_m start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( caligraphic_A start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( bold_italic_x ) + 2 caligraphic_A start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT ( bold_italic_x ) italic_P start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT roman_cos ( italic_θ start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT - italic_ψ start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT ) - 2 caligraphic_A start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT ( bold_italic_x ) italic_p start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT roman_cos ( italic_θ start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT + italic_ψ start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT )
𝒮2(𝒙)2𝒮(𝒙)(m1μ1cosθS+m2μ2sinθS)}φ(𝒙)=0,\displaystyle-{\cal S}^{2}(\hbox{\boldmath${x}$})-2{\cal S}(\hbox{\boldmath${x% }$})(m_{1}\sqrt{\mu}_{1}\cos\theta_{S}+m_{2}\sqrt{\mu}_{2}\sin\theta_{S})\bigg% {\}}\varphi(\hbox{\boldmath${x}$})=0,- caligraphic_S start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( bold_italic_x ) - 2 caligraphic_S ( bold_italic_x ) ( italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT square-root start_ARG italic_μ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT roman_cos italic_θ start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT + italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT square-root start_ARG italic_μ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT roman_sin italic_θ start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT ) } italic_φ ( bold_italic_x ) = 0 , (30)

where tanψA=μ2/μ1subscript𝜓𝐴subscript𝜇2subscript𝜇1\tan\psi_{A}=\sqrt{\mu_{2}/\mu_{1}}roman_tan italic_ψ start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT = square-root start_ARG italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT / italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG.

We take θAsubscript𝜃𝐴\theta_{A}italic_θ start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT and θSsubscript𝜃𝑆\theta_{S}italic_θ start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT so that the following equations takes place

P0cos(θAψA)p0cos(θA+ψA)=εwsuperscript𝑃0subscript𝜃𝐴subscript𝜓𝐴superscript𝑝0subscript𝜃𝐴subscript𝜓𝐴subscript𝜀𝑤\displaystyle P^{0}\cos(\theta_{A}-\psi_{A})-p^{0}\cos(\theta_{A}+\psi_{A})=% \varepsilon_{w}italic_P start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT roman_cos ( italic_θ start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT - italic_ψ start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT ) - italic_p start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT roman_cos ( italic_θ start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT + italic_ψ start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT ) = italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT (31)
m1μ1cosθS+m2μ2sinθS=mw.subscript𝑚1subscript𝜇1subscript𝜃𝑆subscript𝑚2subscript𝜇2subscript𝜃𝑆subscript𝑚𝑤\displaystyle m_{1}\sqrt{\mu}_{1}\cos\theta_{S}+m_{2}\sqrt{\mu}_{2}\sin\theta_% {S}=m_{w}.italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT square-root start_ARG italic_μ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT roman_cos italic_θ start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT + italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT square-root start_ARG italic_μ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT roman_sin italic_θ start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT = italic_m start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT . (32)

Then, the motion equation takes a form

{(εw+𝒜0(𝒙))2+(𝒑+𝑨(𝒙))2(mw+𝒮(𝒙))2}φ(𝒙)=0.superscriptsubscript𝜀𝑤superscript𝒜0𝒙2superscript𝒑𝑨𝒙2superscriptsubscript𝑚𝑤𝒮𝒙2𝜑𝒙0\displaystyle\bigg{\{}(\varepsilon_{w}+{\cal A}^{0}(\hbox{\boldmath${x}$}))^{2% }+(\hbox{\boldmath${p}$}+{\cal\hbox{\boldmath${A}$}}(\hbox{\boldmath${x}$}))^{% 2}-(m_{w}+{\cal S}(\hbox{\boldmath${x}$}))^{2}\bigg{\}}\varphi(\hbox{\boldmath% ${x}$})=0.{ ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT + caligraphic_A start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT ( bold_italic_x ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( bold_italic_p + bold_italic_A ( bold_italic_x ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - ( italic_m start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT + caligraphic_S ( bold_italic_x ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT } italic_φ ( bold_italic_x ) = 0 . (33)

If the solution of Eq.(33) is known the general solution of the problem can be written using Eqs.(19), (21), provided that Eq.(27) also takes place.

In the Breit frame, when 𝒑1+𝒑2=0subscript𝒑1subscript𝒑20\hbox{\boldmath${p}$}_{1}+\hbox{\boldmath${p}$}_{2}=0bold_italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + bold_italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = 0 the obtained equation Eq.(33) coincides with the motion equations have derived earlierCra84 ; Cra87 .

2.4 The two-body problem in external fields

We consider the open two-body problem when the two-particle interaction system is also in external fields. The possibility to factorize the center-mass and relative motion variables, as it has been done above in the absence of external potentials, is entirely determined on the structure of the first terms in Eq.(4), (4b).

In order to the center-mass and relative motion variables could be separated the sums of external potentials A1ν(Xx/2)+A2ν(X+x/2)subscriptsuperscript𝐴𝜈1𝑋𝑥2subscriptsuperscript𝐴𝜈2𝑋𝑥2{A}^{\nu}_{1}(X-x/2)+{A}^{\nu}_{2}(X+x/2)italic_A start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_X - italic_x / 2 ) + italic_A start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_X + italic_x / 2 ) and S1(Xx/2)+S2(X+x/2)subscript𝑆1𝑋𝑥2subscript𝑆2𝑋𝑥2{S}_{1}(X-x/2)+{S}_{2}(X+x/2)italic_S start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_X - italic_x / 2 ) + italic_S start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_X + italic_x / 2 ) should depend only on either X𝑋Xitalic_X or x𝑥xitalic_x. Since such a demand is impossible for the scalar field we exclude the external scalar field from our consideration below , setting S1(Xx/2)=S2(X+x/2)=0subscript𝑆1𝑋𝑥2subscript𝑆2𝑋𝑥20{S}_{1}(X-x/2)={S}_{2}(X+x/2)=0italic_S start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_X - italic_x / 2 ) = italic_S start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_X + italic_x / 2 ) = 0. As for the external vector fields the center-mass and relative motion variables are decoupled when the vector potentials are a linear function of coordinates

A1ν(Xx/2)+A2ν(X+x/2)=eνμαβFμα(x1±x2)β,subscriptsuperscript𝐴𝜈1𝑋𝑥2subscriptsuperscript𝐴𝜈2𝑋𝑥2superscript𝑒𝜈𝜇𝛼𝛽subscript𝐹𝜇𝛼subscriptplus-or-minussubscript𝑥1subscript𝑥2𝛽\displaystyle{A}^{\nu}_{1}(X-x/2)+{A}^{\nu}_{2}(X+x/2)=e^{\nu\mu\alpha\beta}F_% {\mu\alpha}(x_{1}\pm x_{2})_{\beta},italic_A start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_X - italic_x / 2 ) + italic_A start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_X + italic_x / 2 ) = italic_e start_POSTSUPERSCRIPT italic_ν italic_μ italic_α italic_β end_POSTSUPERSCRIPT italic_F start_POSTSUBSCRIPT italic_μ italic_α end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ± italic_x start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_β end_POSTSUBSCRIPT , (34)

where Fμαsubscript𝐹𝜇𝛼F_{\mu\alpha}italic_F start_POSTSUBSCRIPT italic_μ italic_α end_POSTSUBSCRIPT is some 4-tensor being independent on coordinates. When the plus sign occurs the external potential affects the center mass motion, while in the opposite case the external fields influence on the relative motion.

We note that the demand given by Eq.(34) takes place particularly provided that Fμαsubscript𝐹𝜇𝛼F_{\mu\alpha}italic_F start_POSTSUBSCRIPT italic_μ italic_α end_POSTSUBSCRIPT is the tensor of homogeneous electromagnetic fields.

3 Two spin-1212\frac{1}{2}divide start_ARG 1 end_ARG start_ARG 2 end_ARG particles in external fields

We consider spin-1212\frac{1}{2}divide start_ARG 1 end_ARG start_ARG 2 end_ARG particles with masses m1subscript𝑚1m_{1}italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and m2subscript𝑚2m_{2}italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT which interact each other by means of vector and scalar fields , which are Aintν(x)superscriptsubscript𝐴𝑖𝑛𝑡𝜈𝑥A_{int}^{\nu}(x)italic_A start_POSTSUBSCRIPT italic_i italic_n italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT ( italic_x ) and Sint(x)subscript𝑆𝑖𝑛𝑡𝑥S_{int}(x)italic_S start_POSTSUBSCRIPT italic_i italic_n italic_t end_POSTSUBSCRIPT ( italic_x ), respectively. These particles also are in the external vector Aextν(x)superscriptsubscript𝐴𝑒𝑥𝑡𝜈𝑥A_{ext}^{\nu}(x)italic_A start_POSTSUBSCRIPT italic_e italic_x italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT ( italic_x ). Such a two-body system is described by the 16-component spinor which each component φi(εw,𝒓)subscript𝜑𝑖subscript𝜀𝑤𝒓\varphi_{i}(\varepsilon_{w},\hbox{\boldmath${r}$})italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_r ) is governed by the Klein-Gordon equation (2.2). Provided that the interaction and external potentials satisfy Eqs.(27) (34), respectively, these equation is simlified, and φi(εw,𝒓)subscript𝜑𝑖subscript𝜀𝑤𝒓\varphi_{i}(\varepsilon_{w},\hbox{\boldmath${r}$})italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_r ) are governed by Eq.(33)

{(εw+𝒜0(𝒙))2+(𝒑+𝑨(𝒙))2(mw+𝒮(𝒙))2}φi(εw,𝒙)=0.superscriptsubscript𝜀𝑤superscript𝒜0𝒙2superscript𝒑𝑨𝒙2superscriptsubscript𝑚𝑤𝒮𝒙2subscript𝜑𝑖subscript𝜀𝑤𝒙0\displaystyle\bigg{\{}(\varepsilon_{w}+{\cal A}^{0}(\hbox{\boldmath${x}$}))^{2% }+(\hbox{\boldmath${p}$}+{\cal\hbox{\boldmath${A}$}}(\hbox{\boldmath${x}$}))^{% 2}-(m_{w}+{\cal S}(\hbox{\boldmath${x}$}))^{2}\bigg{\}}\varphi_{i}(\varepsilon% _{w},\hbox{\boldmath${x}$})=0.{ ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT + caligraphic_A start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT ( bold_italic_x ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( bold_italic_p + bold_italic_A ( bold_italic_x ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - ( italic_m start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT + caligraphic_S ( bold_italic_x ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT } italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_x ) = 0 . (35)

Following the standard way we rewrite this equation in a form

p^νgνμp^μφi(𝒙)=0,subscript^𝑝𝜈superscript𝑔𝜈𝜇subscript^𝑝𝜇subscript𝜑𝑖𝒙0\displaystyle{\hat{p}}_{\nu}g^{\nu\mu}{\hat{p}}_{\mu}\varphi_{i}(\hbox{% \boldmath${x}$})=0,over^ start_ARG italic_p end_ARG start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT italic_g start_POSTSUPERSCRIPT italic_ν italic_μ end_POSTSUPERSCRIPT over^ start_ARG italic_p end_ARG start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( bold_italic_x ) = 0 , (36)

where we introduce p^ν=((εw+𝒜0(𝒙))2(mw+𝒮(𝒙))2,(𝒑+𝑨(𝒙)))subscript^𝑝𝜈superscriptsubscript𝜀𝑤superscript𝒜0𝒙2superscriptsubscript𝑚𝑤𝒮𝒙2𝒑𝑨𝒙{\hat{p}}_{\nu}=\Big{(}\sqrt{(\varepsilon_{w}+{\cal A}^{0}(\hbox{\boldmath${x}% $}))^{2}-(m_{w}+{\cal S}(\hbox{\boldmath${x}$}))^{2}},-({{\hbox{\boldmath${p}$% }}}+{\cal\hbox{\boldmath${A}$}}(\hbox{\boldmath${x}$}))\Big{)}over^ start_ARG italic_p end_ARG start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT = ( square-root start_ARG ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT + caligraphic_A start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT ( bold_italic_x ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - ( italic_m start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT + caligraphic_S ( bold_italic_x ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , - ( bold_italic_p + bold_italic_A ( bold_italic_x ) ) ); and the metric tensor gνμsuperscript𝑔𝜈𝜇g^{\nu\mu}italic_g start_POSTSUPERSCRIPT italic_ν italic_μ end_POSTSUPERSCRIPT which non-equal to zeroth diagonal componenets are diag=(1,1.1.1){\rm diag}=(1,-1.-1.-1)roman_diag = ( 1 , - 1 . - 1 . - 1 ).

To obtain the motion equation in the Dirac-like form we express the metric tensor in terms of the γ𝛾\gammaitalic_γ-matrices of the (16×16)1616(16\times 16)( 16 × 16 ) dimension. The relations (7) which introduce new variables imply going to a new basis in the γ𝛾\gammaitalic_γ-matrix space so that the new and old γ𝛾\gammaitalic_γ-matrices are related to each other by means of equations

Γ^=μ1γ(1)I(2)+μ2γ(2)I(1),\displaystyle{\hat{\Gamma}}=\mu_{1}\stackrel{{\scriptstyle(1)}}{{\gamma}}% \otimes\stackrel{{\scriptstyle(2)}}{{I}}+\mu_{2}\stackrel{{\scriptstyle(2)}}{{% \gamma}}\otimes\stackrel{{\scriptstyle(1)}}{{I}},over^ start_ARG roman_Γ end_ARG = italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_RELOP SUPERSCRIPTOP start_ARG italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP ⊗ start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 2 ) end_ARG end_RELOP + italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_RELOP SUPERSCRIPTOP start_ARG italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP ⊗ start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 ) end_ARG end_RELOP ,
γ^=γ(2)I(1)γ(1)I(2)\displaystyle{\hat{\gamma}}=\stackrel{{\scriptstyle(2)}}{{{\gamma}}}\otimes% \stackrel{{\scriptstyle(1)}}{{I}}-\stackrel{{\scriptstyle(1)}}{{\gamma}}% \otimes\stackrel{{\scriptstyle(2)}}{{I}}over^ start_ARG italic_γ end_ARG = start_RELOP SUPERSCRIPTOP start_ARG italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP ⊗ start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 ) end_ARG end_RELOP - start_RELOP SUPERSCRIPTOP start_ARG italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP ⊗ start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 2 ) end_ARG end_RELOP (37)

where the hatted γ^^𝛾\hat{\gamma}over^ start_ARG italic_γ end_ARG-matrices form a new basis, the numbers in brackets which are over a symbol means that this matrix acts on the variables of the corresponding particle, I(1,2)superscript𝐼12\stackrel{{\scriptstyle(1,2)}}{{I}}start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 , 2 ) end_ARG end_RELOP is the unit matrix in the space of the first or second particle, respectively. These matrices satisfy the following transformation relations (see Appendix B)

Γ^0Γ^0+Γ^0Γ^0=2I(1+2);superscript^Γ0superscript^Γ0superscript^Γ0superscript^Γ02superscript𝐼12absent\displaystyle{\hat{\Gamma}}^{0}{\hat{\Gamma}}^{0}+{\hat{\Gamma}}^{0}{\hat{% \Gamma}}^{0}=2\stackrel{{\scriptstyle(1+2)}}{{I}};over^ start_ARG roman_Γ end_ARG start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT over^ start_ARG roman_Γ end_ARG start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT + over^ start_ARG roman_Γ end_ARG start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT over^ start_ARG roman_Γ end_ARG start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT = 2 start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 + 2 ) end_ARG end_RELOP ; (38)
Γ^0γa^+γa^Γ^0=0,a=1,2,3;formulae-sequencesuperscript^Γ0^superscript𝛾𝑎^superscript𝛾𝑎superscript^Γ00𝑎123\displaystyle{\hat{\Gamma}}^{0}{\hat{\gamma^{a}}}+{\hat{\gamma^{a}}}{\hat{% \Gamma}}^{0}=0,~{}~{}a=1,2,3;over^ start_ARG roman_Γ end_ARG start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT over^ start_ARG italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG + over^ start_ARG italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG over^ start_ARG roman_Γ end_ARG start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT = 0 , italic_a = 1 , 2 , 3 ; (39)
γa^γb^+γb^γa^=4δabI(1+2)+^superscript𝛾𝑎^superscript𝛾𝑏^superscript𝛾𝑏^superscript𝛾𝑎4superscript𝛿𝑎𝑏superscript𝐼12\displaystyle{\hat{\gamma^{a}}}{\hat{\gamma^{b}}}+{\hat{\gamma^{b}}}{\hat{% \gamma^{a}}}=-4\delta^{ab}\stackrel{{\scriptstyle(1+2)}}{{I}}+over^ start_ARG italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG over^ start_ARG italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG + over^ start_ARG italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG over^ start_ARG italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG = - 4 italic_δ start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 + 2 ) end_ARG end_RELOP +
(σ1aσ2b+σ2bσ1a+σ2aσ1b+σ1bσ2a)I(1+2),superscript𝐼12superscriptsubscript𝜎1𝑎superscriptsubscript𝜎2𝑏superscriptsubscript𝜎2𝑏superscriptsubscript𝜎1𝑎superscriptsubscript𝜎2𝑎superscriptsubscript𝜎1𝑏superscriptsubscript𝜎1𝑏superscriptsubscript𝜎2𝑎absent\displaystyle(\sigma_{1}^{a}\sigma_{2}^{b}+\sigma_{2}^{b}\sigma_{1}^{a}+\sigma% _{2}^{a}\sigma_{1}^{b}+\sigma_{1}^{b}\sigma_{2}^{a})\stackrel{{\scriptstyle(1+% 2)}}{{I}},( italic_σ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT italic_σ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT + italic_σ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT italic_σ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT + italic_σ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT italic_σ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT + italic_σ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT italic_σ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT ) start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 + 2 ) end_ARG end_RELOP , (40)

where 𝒔1subscript𝒔1\hbox{\boldmath${s}$}_{1}bold_italic_s start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and 𝒔2subscript𝒔2\hbox{\boldmath${s}$}_{2}bold_italic_s start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT denote particle spins, δabsuperscript𝛿𝑎𝑏\delta^{ab}italic_δ start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT is the Kronecker symbol, I(1+2)superscript𝐼12\stackrel{{\scriptstyle(1+2)}}{{I}}start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 + 2 ) end_ARG end_RELOP is the 16×16161616\times 1616 × 16 unit matrix. Such transformation relations allow us to write down the metric tensor in terms of new γ𝛾\gammaitalic_γ- matrices as follows

ΓνΓμ+ΓνΓμ=2gνμ+superscriptΓ𝜈superscriptΓ𝜇superscriptΓ𝜈superscriptΓ𝜇limit-from2superscript𝑔𝜈𝜇\displaystyle\Gamma^{\nu}\Gamma^{\mu}+\Gamma^{\nu}\Gamma^{\mu}=2g^{\nu\mu}+roman_Γ start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT roman_Γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT + roman_Γ start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT roman_Γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT = 2 italic_g start_POSTSUPERSCRIPT italic_ν italic_μ end_POSTSUPERSCRIPT + (42)
(σ1aσ2b+σ2bσ1a+σ2aσ1b+σ1bσ2a+4δab)I(1+2)/2,superscript𝐼122superscriptsubscript𝜎1𝑎superscriptsubscript𝜎2𝑏superscriptsubscript𝜎2𝑏superscriptsubscript𝜎1𝑎superscriptsubscript𝜎2𝑎superscriptsubscript𝜎1𝑏superscriptsubscript𝜎1𝑏superscriptsubscript𝜎2𝑎4superscript𝛿𝑎𝑏absent\displaystyle(\sigma_{1}^{a}\sigma_{2}^{b}+\sigma_{2}^{b}\sigma_{1}^{a}+\sigma% _{2}^{a}\sigma_{1}^{b}+\sigma_{1}^{b}\sigma_{2}^{a}+4\delta^{ab})\stackrel{{% \scriptstyle(1+2)/2}}{{I}},( italic_σ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT italic_σ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT + italic_σ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT italic_σ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT + italic_σ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT italic_σ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT + italic_σ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT italic_σ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT + 4 italic_δ start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT ) start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 + 2 ) / 2 end_ARG end_RELOP ,
Γν=(Γ^0,γa^2)(Γ0,𝚪)superscriptΓ𝜈superscript^Γ0^superscript𝛾𝑎2superscriptΓ0𝚪\displaystyle\Gamma^{\nu}=\Big{(}{\hat{\Gamma}}^{0},\frac{{\hat{\gamma^{a}}}}{% 2}\Big{)}\equiv({\Gamma}^{0},{{\hbox{\boldmath${\Gamma}$}}})roman_Γ start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT = ( over^ start_ARG roman_Γ end_ARG start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT , divide start_ARG over^ start_ARG italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG end_ARG start_ARG 2 end_ARG ) ≡ ( roman_Γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT , bold_Γ )

Taking gνμsuperscript𝑔𝜈𝜇g^{\nu\mu}italic_g start_POSTSUPERSCRIPT italic_ν italic_μ end_POSTSUPERSCRIPT from Eq.(42) and substituting it into Eq.(36), we obtain the equation

(Γ0[(εw+𝒜0(𝒙))2(mw+𝒮(𝒙))2]1/2)2limit-fromsuperscriptsuperscriptΓ0superscriptdelimited-[]superscriptsubscript𝜀𝑤superscript𝒜0𝒙2superscriptsubscript𝑚𝑤𝒮𝒙2122\displaystyle\Big{(}\Gamma^{0}[(\varepsilon_{w}+{\cal A}^{0}(\hbox{\boldmath${% x}$}))^{2}-(m_{w}+{\cal S}(\hbox{\boldmath${x}$}))^{2}]^{1/2}\Big{)}^{2}-( roman_Γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT [ ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT + caligraphic_A start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT ( bold_italic_x ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - ( italic_m start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT + caligraphic_S ( bold_italic_x ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT -
(𝚪𝒑^)2)Ψ(Ew,𝒓)=[12(𝝈1+𝝈2)𝒑^]2Ψ(εw,𝒙),\displaystyle(\hbox{\boldmath${\Gamma}$}{\hat{\hbox{\boldmath${p}$}}})^{2}\Big% {)}\Psi(E_{w},\hbox{\boldmath${r}$})=[\frac{1}{2}(\hbox{\boldmath${\sigma}$}_{% 1}+\hbox{\boldmath${\sigma}$}_{2}){\hat{\hbox{\boldmath${p}$}}}]^{2}\Psi(% \varepsilon_{w},\hbox{\boldmath${x}$}),( bold_Γ over^ start_ARG bold_italic_p end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) roman_Ψ ( italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_r ) = [ divide start_ARG 1 end_ARG start_ARG 2 end_ARG ( bold_italic_σ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + bold_italic_σ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) over^ start_ARG bold_italic_p end_ARG ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_Ψ ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_x ) , (43)

where Ψ(εw,𝒓)Ψsubscript𝜀𝑤𝒓\Psi(\varepsilon_{w},\hbox{\boldmath${r}$})roman_Ψ ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_r ) is the 16-component spinor mentioned above in the current section. Eq.(3) leads to the Dirac-like form of the motion equations for this spinor

Γ0[(εw+𝒜0(𝒙))2(mw+𝒮(𝒙))2]1/2Ψ(εw,𝒙)=subscriptΓ0superscriptdelimited-[]superscriptsubscript𝜀𝑤superscript𝒜0𝒙2superscriptsubscript𝑚𝑤𝒮𝒙212Ψsubscript𝜀𝑤𝒙absent\displaystyle\Gamma_{0}[(\varepsilon_{w}+{\cal A}^{0}(\hbox{\boldmath${x}$}))^% {2}-(m_{w}+{\cal S}(\hbox{\boldmath${x}$}))^{2}]^{1/2}\Psi(\varepsilon_{w},% \hbox{\boldmath${x}$})=roman_Γ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT [ ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT + caligraphic_A start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT ( bold_italic_x ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - ( italic_m start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT + caligraphic_S ( bold_italic_x ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT roman_Ψ ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_x ) =
(𝚪+𝑺)𝒑^Ψ(εw,𝒙),𝚪𝑺^𝒑Ψsubscript𝜀𝑤𝒙\displaystyle(\hbox{\boldmath${\Gamma}$}+\hbox{\boldmath${S}$}){\hat{\hbox{% \boldmath${p}$}}}\Psi(\varepsilon_{w},\hbox{\boldmath${x}$}),( bold_Γ + bold_italic_S ) over^ start_ARG bold_italic_p end_ARG roman_Ψ ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_x ) , (44)
Γ0[(εw+𝒜0(𝒙))2(mw+𝒮(𝒙))2]1/2Ψ(εw,𝒙)=subscriptΓ0superscriptdelimited-[]superscriptsubscript𝜀𝑤superscript𝒜0𝒙2superscriptsubscript𝑚𝑤𝒮𝒙212Ψsubscript𝜀𝑤𝒙absent\displaystyle\Gamma_{0}[(\varepsilon_{w}+{\cal A}^{0}(\hbox{\boldmath${x}$}))^% {2}-(m_{w}+{\cal S}(\hbox{\boldmath${x}$}))^{2}]^{1/2}\Psi(\varepsilon_{w},% \hbox{\boldmath${x}$})=roman_Γ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT [ ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT + caligraphic_A start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT ( bold_italic_x ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - ( italic_m start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT + caligraphic_S ( bold_italic_x ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT roman_Ψ ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_x ) =
(𝚪𝑺)𝒑^Ψ(εw,𝒙),𝚪𝑺^𝒑Ψsubscript𝜀𝑤𝒙\displaystyle(\hbox{\boldmath${\Gamma}$}-\hbox{\boldmath${S}$}){\hat{\hbox{% \boldmath${p}$}}}\Psi(\varepsilon_{w},\hbox{\boldmath${x}$}),( bold_Γ - bold_italic_S ) over^ start_ARG bold_italic_p end_ARG roman_Ψ ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_x ) , (45)

where 𝒑=i𝒑𝑖\hbox{\boldmath${p}$}=-i\nablabold_italic_p = - italic_i ∇; 𝒑^=𝒑+𝑨(𝒙)^𝒑𝒑𝑨𝒙{\hat{\hbox{\boldmath${p}$}}}={{\hbox{\boldmath${p}$}}}+{\cal\hbox{\boldmath${% A}$}}(\hbox{\boldmath${x}$})over^ start_ARG bold_italic_p end_ARG = bold_italic_p + bold_italic_A ( bold_italic_x ) , 𝑺=(𝝈1+𝝈2)/2𝑺subscript𝝈1subscript𝝈22\hbox{\boldmath${S}$}=(\hbox{\boldmath${\sigma}$}_{1}+\hbox{\boldmath${\sigma}% $}_{2})/2bold_italic_S = ( bold_italic_σ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + bold_italic_σ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) / 2 in Eqs.(3), (3) is assumed to be multiplied by I(1+2)superscript𝐼12\stackrel{{\scriptstyle(1+2)}}{{I}}start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 + 2 ) end_ARG end_RELOP. It is very important to mention that such a structure of the γμsuperscript𝛾𝜇\gamma^{\mu}italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT matrices defined by Eq.(3) leads to decoupling (see, also Appendix B) the motion equation (3) in four independent equation as compared with the previous researchesSad86-1 ; Sad86-2 ; Sad89 ; Cra84 ; Cra87 ; Cra90 ; Cra01 ; Cra03 ; Cra10 .

Assuming the Eq.(3) is the leading, let us establish the physical sense of Eqs.(3). To clarify it we change the numbers of particles in the two-body system. According to Eq. (3) this procedure consists in changing 𝒑𝒑,𝚪𝚪formulae-sequence𝒑𝒑𝚪𝚪\hbox{\boldmath${p}$}\to-\hbox{\boldmath${p}$},~{}\hbox{\boldmath${\Gamma}$}% \to-\hbox{\boldmath${\Gamma}$}bold_italic_p → - bold_italic_p , bold_Γ → - bold_Γ. Making the changes mentioned above in Eq.(3), we obtain the equation which coincides with Eq.(3). Thus, choosing the leading equation among Eqs.(3) and Eq.(3) means fixation of particle numbering.

Provided that Ψ(εw,𝒙)Ψsubscript𝜀𝑤𝒙\Psi(\varepsilon_{w},\hbox{\boldmath${x}$})roman_Ψ ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_x ) is governed by Eq.(3), the direct calculations give that the Dirac conjugated field Ψ¯(εw,𝒙)=Ψ(εw,𝒙)Γ0¯Ψsubscript𝜀𝑤𝒙superscriptΨsubscript𝜀𝑤𝒙superscriptΓ0{\bar{\Psi}}(\varepsilon_{w},\hbox{\boldmath${x}$})={\Psi}^{\dagger}(% \varepsilon_{w},\hbox{\boldmath${x}$})\Gamma^{0}over¯ start_ARG roman_Ψ end_ARG ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_x ) = roman_Ψ start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_x ) roman_Γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT satisfies an equation

Ψ¯(εw,𝒙)Γ0[(εw+𝒜0(𝒙))2(mw+𝒮(𝒙))2]1/2=¯Ψsubscript𝜀𝑤𝒙subscriptΓ0superscriptdelimited-[]superscriptsubscript𝜀𝑤superscript𝒜0𝒙2superscriptsubscript𝑚𝑤𝒮𝒙212absent\displaystyle{\bar{\Psi}}(\varepsilon_{w},\hbox{\boldmath${x}$})\Gamma_{0}[(% \varepsilon_{w}+{\cal A}^{0}(\hbox{\boldmath${x}$}))^{2}-(m_{w}+{\cal S}(\hbox% {\boldmath${x}$}))^{2}]^{1/2}=over¯ start_ARG roman_Ψ end_ARG ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_x ) roman_Γ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT [ ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT + caligraphic_A start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT ( bold_italic_x ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - ( italic_m start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT + caligraphic_S ( bold_italic_x ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT =
Ψ¯(εw,𝒙)(𝚪+𝑺)𝒑^.¯Ψsubscript𝜀𝑤𝒙𝚪𝑺^𝒑\displaystyle-{\bar{\Psi}}(\varepsilon_{w},\hbox{\boldmath${x}$})(\hbox{% \boldmath${\Gamma}$}+\hbox{\boldmath${S}$}){\hat{{\overleftarrow{\hbox{% \boldmath${p}$}}}}}.- over¯ start_ARG roman_Ψ end_ARG ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_x ) ( bold_Γ + bold_italic_S ) over^ start_ARG over← start_ARG bold_italic_p end_ARG end_ARG . (46)

3.1 Lagrangian and current of fermions in the center-mass and relative motion variables

Eqs.(3), (3) allow us to derive the Lagrangian governing the dynamics of the two-body system in an external field which is

𝔏=12Ψ¯(εw,𝒙){(Γ0[(εw+𝒜0(𝒙))2(mw+𝒮(𝒙))2]1/2\displaystyle{\mathfrak{L}}=\frac{1}{2}{\bar{\Psi}}(\varepsilon_{w},\hbox{% \boldmath${x}$})\Bigg{\{}\Big{(}\Gamma_{0}[(\varepsilon_{w}+{\cal A}^{0}(\hbox% {\boldmath${x}$}))^{2}-(m_{w}+{\cal S}(\hbox{\boldmath${x}$}))^{2}]^{1/2}-fraktur_L = divide start_ARG 1 end_ARG start_ARG 2 end_ARG over¯ start_ARG roman_Ψ end_ARG ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_x ) { ( roman_Γ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT [ ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT + caligraphic_A start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT ( bold_italic_x ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - ( italic_m start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT + caligraphic_S ( bold_italic_x ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT -
(𝚪+𝑺)𝒑^)+(Γ0[(εw+𝒜0(𝒙))2(mw+𝒮(𝒙))2]1/2+\displaystyle(\hbox{\boldmath${\Gamma}$}+\hbox{\boldmath${S}$}){\hat{\hbox{% \boldmath${p}$}}}\Big{)}+\Big{(}\Gamma_{0}[(\varepsilon_{w}+{\cal A}^{0}(\hbox% {\boldmath${x}$}))^{2}-(m_{w}+{\cal S}(\hbox{\boldmath${x}$}))^{2}]^{1/2}+( bold_Γ + bold_italic_S ) over^ start_ARG bold_italic_p end_ARG ) + ( roman_Γ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT [ ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT + caligraphic_A start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT ( bold_italic_x ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - ( italic_m start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT + caligraphic_S ( bold_italic_x ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT +
(𝚪+𝑺)𝒑^)}Ψ(εw,𝒙)\displaystyle(\hbox{\boldmath${\Gamma}$}+\hbox{\boldmath${S}$}){\hat{{% \overleftarrow{\hbox{\boldmath${p}$}}}}}\Big{)}\Bigg{\}}\Psi(\varepsilon_{w},% \hbox{\boldmath${x}$})( bold_Γ + bold_italic_S ) over^ start_ARG over← start_ARG bold_italic_p end_ARG end_ARG ) } roman_Ψ ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_x ) (47)

Such a Lagrangian generates the current

𝑱(Ew,𝒓)=Ψ¯(εw,𝒙)(𝚪+𝑺)Ψ(εw,𝒙),𝑱subscript𝐸𝑤𝒓¯Ψsubscript𝜀𝑤𝒙𝚪𝑺Ψsubscript𝜀𝑤𝒙\displaystyle{\hbox{\boldmath${J}$}}(E_{w},\hbox{\boldmath${r}$})={\bar{\Psi}(% \varepsilon_{w},\hbox{\boldmath${x}$})}(\hbox{\boldmath${\Gamma}$}+\hbox{% \boldmath${S}$})\Psi(\varepsilon_{w},\hbox{\boldmath${x}$}),bold_italic_J ( italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_r ) = over¯ start_ARG roman_Ψ end_ARG ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_x ) ( bold_Γ + bold_italic_S ) roman_Ψ ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_x ) , (48)

whose conservation can be expressed in terms of the continuity equation.

𝑱(εw,𝒙)=0.𝑱subscript𝜀𝑤𝒙0\displaystyle\nabla~{}\hbox{\boldmath${J}$}(\varepsilon_{w},\hbox{\boldmath${x% }$})=0.∇ bold_italic_J ( italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_x ) = 0 . (49)

4 A positronium in a strong magnetic field

We apply the obtained result to a study of observable physical phenomena, by considering a positronium in a strong uniform magnetic field. Such a problem is important in the astrophysics when the electromagnetic signals coming from the astrophysical objects such as pulsars or neutron starsBro ; Bus , whose energy is about 0.511MeV0.511𝑀𝑒𝑉0.511MeV0.511 italic_M italic_e italic_V, likely correspond to the electromagnetic decay of the positronium ground state into two photons.

We study positronium states in a magnetic field in the Breit frame. Then, εw=Ewsubscript𝜀𝑤subscript𝐸𝑤\varepsilon_{w}=E_{w}italic_ε start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT = italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT in left-hand side of Eq.(3). After that we take the scalar potential 𝒮(𝒓)=0𝒮𝒓0{\cal S}(\hbox{\boldmath${r}$})=0caligraphic_S ( bold_italic_r ) = 0 therein, and direct the magnetic field 𝑩𝑩{B}bold_italic_B along the OZ𝑂𝑍OZitalic_O italic_Z axis. Further, we go to the second order equitation, by acting on Eq.(3) with the left-hand side operator of Eq.(3). As a result, we get

(+2(4πα)Ewr+(4πα)2r2+iπα(𝑩×𝝆)\displaystyle\Big{(}\triangle+\frac{2(4\pi\alpha)E_{w}}{r}+\frac{(4\pi\alpha)^% {2}}{r^{2}}+i\sqrt{\pi\alpha}\nabla(\hbox{\boldmath${B}$}\times\hbox{\boldmath% ${\rho}$})-( △ + divide start_ARG 2 ( 4 italic_π italic_α ) italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT end_ARG start_ARG italic_r end_ARG + divide start_ARG ( 4 italic_π italic_α ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + italic_i square-root start_ARG italic_π italic_α end_ARG ∇ ( bold_italic_B × bold_italic_ρ ) - (50)
4παB2ρ2mw2+Ew2)(1001)Ψ(Ew,𝒓)\displaystyle 4\pi\alpha B^{2}\rho^{2}-m_{w}^{2}+E_{w}^{2}\Big{)}\begin{% pmatrix}1&0\\ 0&1\end{pmatrix}\Psi(E_{w},\hbox{\boldmath${r}$})4 italic_π italic_α italic_B start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_m start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) ( start_ARG start_ROW start_CELL 1 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 1 end_CELL end_ROW end_ARG ) roman_Ψ ( italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_r )
=πα(σz,1+σz,200σz,1+σz,2)Ψ(Ew,𝒓),absent𝜋𝛼matrixsubscript𝜎𝑧1subscript𝜎𝑧200subscript𝜎𝑧1subscript𝜎𝑧2Ψsubscript𝐸𝑤𝒓\displaystyle=-\sqrt{\pi\alpha}\begin{pmatrix}\sigma_{z,1}+\sigma_{z,2}&0\\ 0&\sigma_{z,1}+\sigma_{z,2}\end{pmatrix}\Psi(E_{w},\hbox{\boldmath${r}$}),= - square-root start_ARG italic_π italic_α end_ARG ( start_ARG start_ROW start_CELL italic_σ start_POSTSUBSCRIPT italic_z , 1 end_POSTSUBSCRIPT + italic_σ start_POSTSUBSCRIPT italic_z , 2 end_POSTSUBSCRIPT end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL italic_σ start_POSTSUBSCRIPT italic_z , 1 end_POSTSUBSCRIPT + italic_σ start_POSTSUBSCRIPT italic_z , 2 end_POSTSUBSCRIPT end_CELL end_ROW end_ARG ) roman_Ψ ( italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_r ) ,

where α=e2/4π𝛼superscript𝑒24𝜋\alpha=e^{2}/4\piitalic_α = italic_e start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / 4 italic_π is the fine structure constant; σz(1,2)subscript𝜎𝑧12\sigma_{z(1,2)}italic_σ start_POSTSUBSCRIPT italic_z ( 1 , 2 ) end_POSTSUBSCRIPT are the Pauli matrices acting on the spin variables of the first (an electron) and the second (a positron) particle. In obtaining the latest equation we have taken

Aintμ=(Aint0,0)=(4παr,0);𝒜ext=12(𝑩×𝒓).formulae-sequencesubscriptsuperscript𝐴𝜇𝑖𝑛𝑡subscriptsuperscript𝐴0𝑖𝑛𝑡04𝜋𝛼𝑟0subscript𝒜𝑒𝑥𝑡12𝑩𝒓\displaystyle A^{\mu}_{int}=(A^{0}_{int},0)=(-\frac{4\pi\alpha}{r},0);~{}{\cal A% }_{ext}=\frac{1}{2}(\hbox{\boldmath${B}$}\times\hbox{\boldmath${r}$}).italic_A start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i italic_n italic_t end_POSTSUBSCRIPT = ( italic_A start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i italic_n italic_t end_POSTSUBSCRIPT , 0 ) = ( - divide start_ARG 4 italic_π italic_α end_ARG start_ARG italic_r end_ARG , 0 ) ; caligraphic_A start_POSTSUBSCRIPT italic_e italic_x italic_t end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG ( bold_italic_B × bold_italic_r ) . (51)

Eq.(50) contains two diagonal matrices that allows us to go to two-component functions Ψ(Ew,𝒓)Ψsubscript𝐸𝑤𝒓\Psi(E_{w},\hbox{\boldmath${r}$})roman_Ψ ( italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_r )

(+2(4πα)Ewr+(4πα)2r2+iπα(𝑩×𝒓)\displaystyle\Big{(}\triangle+\frac{2(4\pi\alpha)E_{w}}{r}+\frac{(4\pi\alpha)^% {2}}{r^{2}}+i\sqrt{\pi\alpha}\nabla(\hbox{\boldmath${B}$}\times\hbox{\boldmath% ${r}$})-( △ + divide start_ARG 2 ( 4 italic_π italic_α ) italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT end_ARG start_ARG italic_r end_ARG + divide start_ARG ( 4 italic_π italic_α ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + italic_i square-root start_ARG italic_π italic_α end_ARG ∇ ( bold_italic_B × bold_italic_r ) -
4παB2ρ2mw2+Ew2)Ψ(Ew,𝒓)\displaystyle 4\pi\alpha B^{2}\rho^{2}-m_{w}^{2}+E_{w}^{2}\Big{)}\Psi(E_{w},% \hbox{\boldmath${r}$})4 italic_π italic_α italic_B start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_m start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) roman_Ψ ( italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_r )
=πα(σz,1+σz,2)Ψ(Ew,𝒓).absent𝜋𝛼subscript𝜎𝑧1subscript𝜎𝑧2Ψsubscript𝐸𝑤𝒓\displaystyle=-\sqrt{\pi\alpha}(\sigma_{z,1}+\sigma_{z,2})\Psi(E_{w},\hbox{% \boldmath${r}$}).= - square-root start_ARG italic_π italic_α end_ARG ( italic_σ start_POSTSUBSCRIPT italic_z , 1 end_POSTSUBSCRIPT + italic_σ start_POSTSUBSCRIPT italic_z , 2 end_POSTSUBSCRIPT ) roman_Ψ ( italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_r ) . (52)

We study the low energy states (S=Sz=0𝑆subscript𝑆𝑧0S=S_{z}=0italic_S = italic_S start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = 0) of a positronium, which are differed by the spin number S𝑆Sitalic_S, in a strong magnetic field 𝑩𝑩{B}bold_italic_B when the magnetic length a=(2παB)1/2𝑎superscript2𝜋𝛼𝐵12a=(2\sqrt{\pi\alpha}B)^{-1/2}italic_a = ( 2 square-root start_ARG italic_π italic_α end_ARG italic_B ) start_POSTSUPERSCRIPT - 1 / 2 end_POSTSUPERSCRIPT is much smaller than the positronium Bohr radius aB=1/2παmesubscript𝑎𝐵12𝜋𝛼subscript𝑚𝑒a_{B}=1/2\pi\alpha m_{e}italic_a start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT = 1 / 2 italic_π italic_α italic_m start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT

aaB,much-less-than𝑎subscript𝑎𝐵\displaystyle~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{% }~{}a\ll a_{B},italic_a ≪ italic_a start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT , (53)

where mesubscript𝑚𝑒m_{e}italic_m start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT is the electron mass.

However, the ground state (S=Sz=0𝑆subscript𝑆𝑧0S=S_{z}=0italic_S = italic_S start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = 0) turn out to be mixedLL4 with the triplet spin state (S=1,Sz=0)formulae-sequence𝑆1subscript𝑆𝑧0(S=1,S_{z}=0)( italic_S = 1 , italic_S start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = 0 ) due to the right-hand side terms in Eq.(4). Therefore, we look for solution of this equation as follows

Ψ(Ew,𝒓)=ψm=0,n=0(𝝆)2S=01((10)1(01)2+(1)S+1(01)1(10)2)ψS(z),Ψsubscript𝐸𝑤𝒓subscript𝜓formulae-sequence𝑚0𝑛0𝝆2superscriptsubscript𝑆01subscriptmatrix101subscriptmatrix012superscript1𝑆1subscriptmatrix011subscriptmatrix102subscript𝜓𝑆𝑧\displaystyle\Psi(E_{w},\hbox{\boldmath${r}$})=\frac{\psi_{m=0,n=0}(\hbox{% \boldmath${\rho}$})}{\sqrt{2}}\sum\limits_{S=0}^{1}\Bigg{(}\begin{pmatrix}1\\ 0\end{pmatrix}_{1}\begin{pmatrix}0\\ 1\end{pmatrix}_{2}+(-1)^{S+1}\begin{pmatrix}0\\ 1\end{pmatrix}_{1}\begin{pmatrix}1\\ 0\end{pmatrix}_{2}\Bigg{)}\psi_{S}(z),roman_Ψ ( italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_r ) = divide start_ARG italic_ψ start_POSTSUBSCRIPT italic_m = 0 , italic_n = 0 end_POSTSUBSCRIPT ( bold_italic_ρ ) end_ARG start_ARG square-root start_ARG 2 end_ARG end_ARG ∑ start_POSTSUBSCRIPT italic_S = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT ( ( start_ARG start_ROW start_CELL 1 end_CELL end_ROW start_ROW start_CELL 0 end_CELL end_ROW end_ARG ) start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( start_ARG start_ROW start_CELL 0 end_CELL end_ROW start_ROW start_CELL 1 end_CELL end_ROW end_ARG ) start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + ( - 1 ) start_POSTSUPERSCRIPT italic_S + 1 end_POSTSUPERSCRIPT ( start_ARG start_ROW start_CELL 0 end_CELL end_ROW start_ROW start_CELL 1 end_CELL end_ROW end_ARG ) start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( start_ARG start_ROW start_CELL 1 end_CELL end_ROW start_ROW start_CELL 0 end_CELL end_ROW end_ARG ) start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) italic_ψ start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT ( italic_z ) , (54)

where the indexes mean the number of a particle; ψm=0,n=0(𝝆)subscript𝜓formulae-sequence𝑚0𝑛0𝝆\psi_{m=0,n=0}(\hbox{\boldmath${\rho}$})italic_ψ start_POSTSUBSCRIPT italic_m = 0 , italic_n = 0 end_POSTSUBSCRIPT ( bold_italic_ρ ) is the wave function governing the transverse motion of a particle in the ground state in a magnetic fieldLL3 , which is equal to

ψm=0,n=0(𝝆)=1a2πexp(ρ24a2).subscript𝜓formulae-sequence𝑚0𝑛0𝝆1𝑎2𝜋superscript𝜌24superscript𝑎2\displaystyle\psi_{m=0,n=0}(\hbox{\boldmath${\rho}$})=\frac{1}{a\sqrt{2\pi}}% \exp{\left(-\frac{\rho^{2}}{4a^{2}}\right)}.italic_ψ start_POSTSUBSCRIPT italic_m = 0 , italic_n = 0 end_POSTSUBSCRIPT ( bold_italic_ρ ) = divide start_ARG 1 end_ARG start_ARG italic_a square-root start_ARG 2 italic_π end_ARG end_ARG roman_exp ( - divide start_ARG italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 4 italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) . (55)

Here the cylindrical coordinates 𝒓=𝝆+𝒆zz𝒓𝝆subscript𝒆𝑧𝑧\hbox{\boldmath${r}$}=\hbox{\boldmath${\rho}$}+\hbox{\boldmath${e}$}_{z}zbold_italic_r = bold_italic_ρ + bold_italic_e start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_z are introduced, m𝑚mitalic_m is the projection of an angular momentum along the 𝑩𝑩{B}bold_italic_B direction, n𝑛nitalic_n is a radial quantum number. Substituting Ψ(Ew,𝒓)Ψsubscript𝐸𝑤𝒓\Psi(E_{w},\hbox{\boldmath${r}$})roman_Ψ ( italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT , bold_italic_r ) given by Eq.(54) into Eq.(4), we derive

(d2dξ2+21/2π(4πα)Ewaeξ2/2erfc(|ξ|/2)+\displaystyle\Big{(}\frac{d^{2}}{d\xi^{2}}+2^{1/2}\sqrt{\pi}(4\pi\alpha)E_{w}a% ~{}e^{\xi^{2}/2}erfc(|\xi|/\sqrt{2})+( divide start_ARG italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_d italic_ξ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + 2 start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT square-root start_ARG italic_π end_ARG ( 4 italic_π italic_α ) italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT italic_a italic_e start_POSTSUPERSCRIPT italic_ξ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / 2 end_POSTSUPERSCRIPT italic_e italic_r italic_f italic_c ( | italic_ξ | / square-root start_ARG 2 end_ARG ) +
(4πα)2eξ2/22E1(ξ2/2)εS2)ψS(ξ)=2παψ1S(ξ),\displaystyle\frac{(4\pi\alpha)^{2}e^{\xi^{2}/2}}{2}E_{1}(\xi^{2}/2)-% \varepsilon^{2}_{S}\Big{)}\psi_{S}(\xi)=-2\sqrt{\pi\alpha}\psi_{1-S}(\xi),divide start_ARG ( 4 italic_π italic_α ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT italic_ξ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 end_ARG italic_E start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_ξ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / 2 ) - italic_ε start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT ) italic_ψ start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT ( italic_ξ ) = - 2 square-root start_ARG italic_π italic_α end_ARG italic_ψ start_POSTSUBSCRIPT 1 - italic_S end_POSTSUBSCRIPT ( italic_ξ ) ,

where the introduced erfc(x)𝑒𝑟𝑓𝑐𝑥erfc(x)italic_e italic_r italic_f italic_c ( italic_x ) is the complementary error function, E1(x)subscript𝐸1𝑥E_{1}(x)italic_E start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x ) is the integral exponentGra . We have also denoted in Eq.(4)

z/aξ,(mw2Ew2)a2+1+(1)S+1εS2.formulae-sequence𝑧𝑎𝜉superscriptsubscript𝑚𝑤2superscriptsubscript𝐸𝑤2superscript𝑎21superscript1𝑆1superscriptsubscript𝜀𝑆2\displaystyle z/a\equiv\xi,~{}~{}(m_{w}^{2}-E_{w}^{2})a^{2}+1+(-1)^{S+1}\equiv% \varepsilon_{S}^{2}.italic_z / italic_a ≡ italic_ξ , ( italic_m start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 1 + ( - 1 ) start_POSTSUPERSCRIPT italic_S + 1 end_POSTSUPERSCRIPT ≡ italic_ε start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT . (57)

Since the magnetic field is strong the particle222The compound particle - a positronium - is meant here. is weakly bound with respect to a motion along the OZ𝑂𝑍OZitalic_O italic_Z axis as compared with the transverse motion. Besides that, the term

U(ξ)(21/2π(4πα)Ewaeξ2/2erfc(|ξ|/2)+\displaystyle-U(\xi)\equiv\Bigg{(}2^{1/2}\sqrt{\pi}(4\pi\alpha)E_{w}ae^{\xi^{2% }/2}erfc(|\xi|/\sqrt{2})+- italic_U ( italic_ξ ) ≡ ( 2 start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT square-root start_ARG italic_π end_ARG ( 4 italic_π italic_α ) italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT italic_a italic_e start_POSTSUPERSCRIPT italic_ξ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / 2 end_POSTSUPERSCRIPT italic_e italic_r italic_f italic_c ( | italic_ξ | / square-root start_ARG 2 end_ARG ) +
(4πα)2eξ2/22E1(ξ2/2))\displaystyle\frac{(4\pi\alpha)^{2}e^{\xi^{2}/2}}{2}E_{1}(\xi^{2}/2)\Bigg{)}divide start_ARG ( 4 italic_π italic_α ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT italic_ξ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 end_ARG italic_E start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_ξ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / 2 ) ) (58)

is weakly singular since U(ξ)1/ξproportional-to𝑈𝜉1𝜉U(\xi)\propto-1/\sqrt{\xi}italic_U ( italic_ξ ) ∝ - 1 / square-root start_ARG italic_ξ end_ARG at small ξ0𝜉0\xi\to 0italic_ξ → 0. Therefore, we change the exact potential U(ξ)𝑈𝜉U(\xi)italic_U ( italic_ξ ) by its approximate form

U(ξ)δ(ξ)+ψS(ξ)U(ξ)ψS(ξ)𝑑ξ,similar-to-or-equals𝑈𝜉𝛿𝜉superscriptsubscriptsubscriptsuperscript𝜓𝑆𝜉𝑈𝜉subscript𝜓𝑆𝜉differential-d𝜉\displaystyle U(\xi)\simeq\delta(\xi)\int\limits_{-\infty}^{+\infty}\psi^{\ast% }_{S}(\xi)U(\xi)\psi_{S}(\xi)d\xi,italic_U ( italic_ξ ) ≃ italic_δ ( italic_ξ ) ∫ start_POSTSUBSCRIPT - ∞ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + ∞ end_POSTSUPERSCRIPT italic_ψ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT ( italic_ξ ) italic_U ( italic_ξ ) italic_ψ start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT ( italic_ξ ) italic_d italic_ξ , (59)

where ψS(ξ)subscript𝜓𝑆𝜉\psi_{S}(\xi)italic_ψ start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT ( italic_ξ ) satisfies an equation

(d2dξ2δ(ξ)(+ψS(ξ)U(ξ)ψS(ξ)𝑑ξ)εS2)ψS(ξ)=superscript𝑑2𝑑superscript𝜉2𝛿𝜉superscriptsubscriptsubscriptsuperscript𝜓𝑆superscript𝜉𝑈superscript𝜉subscript𝜓𝑆superscript𝜉differential-dsuperscript𝜉subscriptsuperscript𝜀2𝑆subscript𝜓𝑆𝜉absent\displaystyle\Big{(}\frac{d^{2}}{d\xi^{2}}-\delta(\xi)(\int\limits_{-\infty}^{% +\infty}\psi^{\ast}_{S}(\xi^{\prime})U(\xi^{\prime})\psi_{S}(\xi^{\prime})d\xi% ^{\prime})-\varepsilon^{2}_{S}\Big{)}\psi_{S}(\xi)=( divide start_ARG italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_d italic_ξ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG - italic_δ ( italic_ξ ) ( ∫ start_POSTSUBSCRIPT - ∞ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + ∞ end_POSTSUPERSCRIPT italic_ψ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT ( italic_ξ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) italic_U ( italic_ξ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) italic_ψ start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT ( italic_ξ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) italic_d italic_ξ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) - italic_ε start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT ) italic_ψ start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT ( italic_ξ ) =
2παψ1S(ξ),2𝜋𝛼subscript𝜓1𝑆𝜉\displaystyle-2\sqrt{\pi\alpha}\psi_{1-S}(\xi),- 2 square-root start_ARG italic_π italic_α end_ARG italic_ψ start_POSTSUBSCRIPT 1 - italic_S end_POSTSUBSCRIPT ( italic_ξ ) ,

We substitute ψS(ξ)subscript𝜓𝑆𝜉\psi_{S}(\xi)italic_ψ start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT ( italic_ξ ), given by the formula

ψS(ξ)=kSexp(kS|ξ|),subscript𝜓𝑆𝜉subscript𝑘𝑆subscript𝑘𝑆𝜉\displaystyle\psi_{S}(\xi)=\sqrt{k_{S}}\exp(-k_{S}|\xi|),italic_ψ start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT ( italic_ξ ) = square-root start_ARG italic_k start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT end_ARG roman_exp ( - italic_k start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT | italic_ξ | ) , (61)

into Eq.(4), where kS=εSsubscript𝑘𝑆subscript𝜀𝑆k_{S}=\sqrt{\varepsilon}_{S}italic_k start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT = square-root start_ARG italic_ε end_ARG start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT is a positive constant. As a result, we derive the equation to obtain kSsubscript𝑘𝑆k_{S}italic_k start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT

kS=(4πα)20𝑑ξexp(ξ2/22kSξ)E1(ξ2/2)+subscript𝑘𝑆limit-fromsuperscript4𝜋𝛼2superscriptsubscript0differential-d𝜉superscript𝜉222subscript𝑘𝑆𝜉subscript𝐸1superscript𝜉22\displaystyle k_{S}={(4\pi\alpha)^{2}}\int\limits_{0}^{\infty}d\xi\exp{({\xi^{% 2}/2}-2k_{S}\xi)}E_{1}(\xi^{2}/2)+italic_k start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT = ( 4 italic_π italic_α ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT italic_d italic_ξ roman_exp ( italic_ξ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / 2 - 2 italic_k start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT italic_ξ ) italic_E start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_ξ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / 2 ) +
23/2π(4πα)Ewa0𝑑ξexp(ξ2/22kSξ)erfc(ξ/2)superscript232𝜋4𝜋𝛼subscript𝐸𝑤𝑎superscriptsubscript0differential-d𝜉superscript𝜉222subscript𝑘𝑆𝜉𝑒𝑟𝑓𝑐𝜉2\displaystyle 2^{3/2}\sqrt{\pi}(4\pi\alpha)E_{w}a\int\limits_{0}^{\infty}d\xi% \exp{({\xi^{2}/2}-2k_{S}\xi)}erfc(\xi/\sqrt{2})2 start_POSTSUPERSCRIPT 3 / 2 end_POSTSUPERSCRIPT square-root start_ARG italic_π end_ARG ( 4 italic_π italic_α ) italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT italic_a ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT italic_d italic_ξ roman_exp ( italic_ξ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / 2 - 2 italic_k start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT italic_ξ ) italic_e italic_r italic_f italic_c ( italic_ξ / square-root start_ARG 2 end_ARG )

The validity limit of the derived equation is kS1much-less-thansubscript𝑘𝑆1k_{S}\ll 1italic_k start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT ≪ 1. At such kSsubscript𝑘𝑆k_{S}italic_k start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT the second integral is approximately π3/2/2superscript𝜋322\pi^{3/2}/\sqrt{2}italic_π start_POSTSUPERSCRIPT 3 / 2 end_POSTSUPERSCRIPT / square-root start_ARG 2 end_ARG Gra , whereas the first integral can be approximately taken to be equal to ln(aB/a)/πsubscript𝑎𝐵𝑎𝜋\ln(a_{B}/a)/\sqrt{\pi}roman_ln ( italic_a start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT / italic_a ) / square-root start_ARG italic_π end_ARGLL3 after the dimensional regularization.

As a result, we find

kS=23/2(4πα)Ewaln(aB/a)+(4πα)2π3/22.subscript𝑘𝑆superscript2324𝜋𝛼subscript𝐸𝑤𝑎subscript𝑎𝐵𝑎superscript4𝜋𝛼2superscript𝜋322\displaystyle k_{S}=2^{3/2}(4\pi\alpha)E_{w}a\ln(a_{B}/a)+(4\pi\alpha)^{2}% \frac{\pi^{3/2}}{\sqrt{2}}.italic_k start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT = 2 start_POSTSUPERSCRIPT 3 / 2 end_POSTSUPERSCRIPT ( 4 italic_π italic_α ) italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT italic_a roman_ln ( italic_a start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT / italic_a ) + ( 4 italic_π italic_α ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT divide start_ARG italic_π start_POSTSUPERSCRIPT 3 / 2 end_POSTSUPERSCRIPT end_ARG start_ARG square-root start_ARG 2 end_ARG end_ARG . (63)

In the non-relativistic case meaαmuch-greater-thansubscript𝑚𝑒𝑎𝛼m_{e}a\gg\alphaitalic_m start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_a ≫ italic_α we obtain kS=23/2(4πα/3)mealn(aB/a)subscript𝑘𝑆superscript2324𝜋𝛼3subscript𝑚𝑒𝑎subscript𝑎𝐵𝑎k_{S}=2^{3/2}(4\pi\alpha/3)m_{e}a\ln(a_{B}/a)italic_k start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT = 2 start_POSTSUPERSCRIPT 3 / 2 end_POSTSUPERSCRIPT ( 4 italic_π italic_α / 3 ) italic_m start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_a roman_ln ( italic_a start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT / italic_a ). Taking into account of Eq.(63), we derive

Ew=me2(1+2(1+(1)S+1)(mea)2\displaystyle E_{w}=\frac{m_{e}}{2}\Big{(}1+\frac{2(1+(-1)^{S+1})}{(m_{e}a)^{2% }}-italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT = divide start_ARG italic_m start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG ( 1 + divide start_ARG 2 ( 1 + ( - 1 ) start_POSTSUPERSCRIPT italic_S + 1 end_POSTSUPERSCRIPT ) end_ARG start_ARG ( italic_m start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_a ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG -
2(4πα/3)2ln2(aB/a)),\displaystyle 2(4\pi\alpha/3)^{2}\ln^{2}(a_{B}/a)\Big{)},2 ( 4 italic_π italic_α / 3 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_ln start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_a start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT / italic_a ) ) , (64)

that relates to the results obtained earlierEll .

In the opposite limiting situation meaαmuch-less-thansubscript𝑚𝑒𝑎𝛼m_{e}a\ll\alphaitalic_m start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_a ≪ italic_α, but provided that the vacuum polarization effects in a strong magnetic field are absentLL4 , we get kS=9(2π)7/2α2subscript𝑘𝑆9superscript2𝜋72superscript𝛼2k_{S}=9(2\pi)^{7/2}\alpha^{2}italic_k start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT = 9 ( 2 italic_π ) start_POSTSUPERSCRIPT 7 / 2 end_POSTSUPERSCRIPT italic_α start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, that gives

Ew1a2((1+(1)S+1)(2π)7α4/162).similar-to-or-equalssubscript𝐸𝑤1superscript𝑎21superscript1𝑆1superscript2𝜋7superscript𝛼4162\displaystyle E_{w}\simeq\frac{1}{a^{2}}((1+(-1)^{S+1})-(2\pi)^{7}\alpha^{4}/1% 62).italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT ≃ divide start_ARG 1 end_ARG start_ARG italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ( ( 1 + ( - 1 ) start_POSTSUPERSCRIPT italic_S + 1 end_POSTSUPERSCRIPT ) - ( 2 italic_π ) start_POSTSUPERSCRIPT 7 end_POSTSUPERSCRIPT italic_α start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT / 162 ) . (65)

The obtained Ewsubscript𝐸𝑤E_{w}italic_E start_POSTSUBSCRIPT italic_w end_POSTSUBSCRIPT sufficiently differs from the results has been derived above in the non-relativistic case. Eqs.(4), (65) show the stronger a magnetic field , the less the positronium level depth under the bottom of the main Landau zone.

4.1 Decay of the positronium ground state in a strong magnetic field

The decay width ΓΓ\Gammaroman_Γ of a particle is proportional toGreen

Γ|ψ(0)|2M2,similar-to-or-equalsΓsuperscript𝜓02superscript𝑀2\displaystyle\Gamma\simeq\frac{|\psi(0)|^{2}}{M^{2}},roman_Γ ≃ divide start_ARG | italic_ψ ( 0 ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_M start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , (66)

where ψ(0)𝜓0\psi(0)italic_ψ ( 0 ) and M𝑀Mitalic_M are the wave function and mass of the decaying particle, respectively.

According to Eqs.(55), (61) the squared wave function module |ψ(0)|2superscript𝜓02|\psi(0)|^{2}| italic_ψ ( 0 ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, depending on the absence or presence of a magnetic field, is of the order of

|ψB=0(0)|2aB3|ψB0(0)|21aBa2.formulae-sequencesimilar-tosuperscriptsubscript𝜓𝐵002superscriptsubscript𝑎𝐵3similar-tosuperscriptsubscript𝜓𝐵0021subscript𝑎𝐵superscript𝑎2\displaystyle|\psi_{B=0}(0)|^{2}\sim a_{B}^{-3}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}|% \psi_{B\neq 0}(0)|^{2}\sim\frac{1}{a_{B}a^{2}}.| italic_ψ start_POSTSUBSCRIPT italic_B = 0 end_POSTSUBSCRIPT ( 0 ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ∼ italic_a start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT | italic_ψ start_POSTSUBSCRIPT italic_B ≠ 0 end_POSTSUBSCRIPT ( 0 ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ∼ divide start_ARG 1 end_ARG start_ARG italic_a start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG . (67)

The estimations given by Eq.(67) show that a strong magnetic field leads to sufficient increasing, in (aB/a)21much-greater-thansuperscriptsubscript𝑎𝐵𝑎21(a_{B}/a)^{2}\gg 1( italic_a start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT / italic_a ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ≫ 1 times, the decay width as compared with the case B=0𝐵0B=0italic_B = 0. Such a behavior of the decay width in a magnetic field is a result of the additional compressing the electron-positron bound state by a magnetic field in the plane which is perpendicular to the 𝐁𝐁{\bf B}bold_B-direction.

5 Conclusion

In terms of the principle of least action we study the open relativistic two-body problem, when two interacting particles also undergo effects of external potentials. Starting from the single particle action integral for a spinless particle we derive the Lagrangian governing the dynamics of both the two-body system of scalar particles and each of the 16 components of a spinor determoning the two-fermion system in external fields, beyond the consideration in the Breit frame. The motion equations, generated such a Lagrangian, take into account the modification of Newton’s third law due to in external fields, and have, in the general case, the operator structure in the variables of the relative momentum and the momentum of the center mass of the two-body system. Provided that Newton’s third law takes place the derived equations are simplified, and turn out in a good relation to the equations have been obtained earlier in the constraint dynamic approach Cra84 ; Cra87 ; Cra90 ; Cra01 ; Cra03 ; Cra10 in the absence of external fields. Following the standard way, based on using the γ𝛾\gammaitalic_γ-matrix basis, the Klein-Gordon equation governing the evolution of each of the 16 components of the relativistic spinor is reduced to the Dirac-like equations. The new γ𝛾\gammaitalic_γ-matrix basis, when the derived Dirac-like equations turn out to be decoupled in four independent coupled equations which determine the dynamics of four-component spinors, is found . Based on the derived Dirac-like equations we obtain the Lagrangian which generates the conserving fermion current. In the case of the pure QED the developed approach has been applied to studying a positronium in a strong magnetic field. The positronium low energy states are studied, taking into account spin-mixing effectLL4 . It is shown that strong magnetic fields result in essentially reconstructing the spectrum of the positronium states and decreasing the life-time of the positronium ground state.

Appendix A The transformation relation for ΓΓ\Gammaroman_Γ- matrices

According to Eqs.(3) we have

Γ^=μ1γ(1)I(2)+μ2γ(2)I(1),\displaystyle{\hat{\Gamma}}=\mu_{1}\stackrel{{\scriptstyle(1)}}{{\gamma}}% \otimes\stackrel{{\scriptstyle(2)}}{{I}}+\mu_{2}\stackrel{{\scriptstyle(2)}}{{% \gamma}}\otimes\stackrel{{\scriptstyle(1)}}{{I}},over^ start_ARG roman_Γ end_ARG = italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_RELOP SUPERSCRIPTOP start_ARG italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP ⊗ start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 2 ) end_ARG end_RELOP + italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_RELOP SUPERSCRIPTOP start_ARG italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP ⊗ start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 ) end_ARG end_RELOP ,
𝜸^=𝜸(2)I(1)𝜸(1)I(2).\displaystyle{\hat{\hbox{\boldmath${\gamma}$}}}=\stackrel{{\scriptstyle(2)}}{{% {\hbox{\boldmath${\gamma}$}}}}\otimes\stackrel{{\scriptstyle(1)}}{{I}}-% \stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}}}\otimes\stackrel{{% \scriptstyle(2)}}{{I}}.over^ start_ARG bold_italic_γ end_ARG = start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP ⊗ start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 ) end_ARG end_RELOP - start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP ⊗ start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 2 ) end_ARG end_RELOP . (68)

The matrices Γ^^Γ{\hat{\Gamma}}over^ start_ARG roman_Γ end_ARG and 𝜸^^𝜸\hat{\hbox{\boldmath${\gamma}$}}over^ start_ARG bold_italic_γ end_ARG are explicitly given by formulas

Γ^0=(μ1γ0(1)0I(2)+μ2γ0(2)0I(1))=μ1(γ0(1)00000γ0(1)00000γ0(1)00000γ0(1)0)+μ2(γ0(2)00000γ0(2)00000γ0(2)00000γ0(2)0)\displaystyle{\hat{\Gamma}}^{0}=(\mu_{1}\stackrel{{\scriptstyle(1)}}{{\gamma^{% 0}}}\otimes\stackrel{{\scriptstyle(2)}}{{I}}+\mu_{2}\stackrel{{\scriptstyle(2)% }}{{\gamma^{0}}}\otimes\stackrel{{\scriptstyle(1)}}{{I}})=\mu_{1}\begin{% pmatrix}\stackrel{{\scriptstyle(1)}}{{\gamma^{0}}}&0&0&0\\ 0&\stackrel{{\scriptstyle(1)}}{{\gamma^{0}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(1)}}{{\gamma^{0}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(1)}}{{\gamma^{0}}}\\ \end{pmatrix}+\mu_{2}\begin{pmatrix}\stackrel{{\scriptstyle(2)}}{{\gamma^{0}}}% &0&0&0\\ 0&\stackrel{{\scriptstyle(2)}}{{\gamma^{0}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(2)}}{{\gamma^{0}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(2)}}{{\gamma^{0}}}\\ \end{pmatrix}over^ start_ARG roman_Γ end_ARG start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT = ( italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP ⊗ start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 2 ) end_ARG end_RELOP + italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP ⊗ start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 ) end_ARG end_RELOP ) = italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL end_ROW end_ARG ) + italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL end_ROW end_ARG ) (69)
𝜸^=(𝜸(2)0000𝜸(2)0000𝜸(2)0000𝜸(2))(𝜸(1)0000𝜸(1)0000𝜸(1)0000𝜸(1)),^𝜸matrixsuperscript𝜸20000superscript𝜸20000superscript𝜸20000superscript𝜸2matrixsuperscript𝜸10000superscript𝜸10000superscript𝜸10000superscript𝜸1\displaystyle{\hat{\hbox{\boldmath${\gamma}$}}}=\begin{pmatrix}\stackrel{{% \scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}}}&0&0&0\\ 0&\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}}}\\ \end{pmatrix}-\begin{pmatrix}\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${% \gamma}$}}}&0&0&0\\ 0&\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}}}\\ \end{pmatrix},over^ start_ARG bold_italic_γ end_ARG = ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL end_ROW end_ARG ) - ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL end_ROW end_ARG ) , (70)

where a number in the brackets which is over γ𝛾\gammaitalic_γ means that this matrix acts on the variable of the corresponding particle, I(a)superscript𝐼𝑎\stackrel{{\scriptstyle(a)}}{{I}}start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( italic_a ) end_ARG end_RELOP is the unit matrix in the space of the a𝑎aitalic_a particle. Then, we obtain

(Γ^0Γ^0+Γ^0Γ^0)=2Γ^0Γ^0=2(μ1γ0(1)0I(2)+μ2γ0(2)0I(1))×\displaystyle({\hat{\Gamma}}^{0}{\hat{\Gamma}}_{0}+{\hat{\Gamma}}^{0}{\hat{% \Gamma}}_{0})=2{\hat{\Gamma}}^{0}{\hat{\Gamma}}_{0}=2(\mu_{1}\stackrel{{% \scriptstyle(1)}}{{\gamma^{0}}}\otimes\stackrel{{\scriptstyle(2)}}{{I}}+\mu_{2% }\stackrel{{\scriptstyle(2)}}{{\gamma^{0}}}\otimes\stackrel{{\scriptstyle(1)}}% {{I}})\times( over^ start_ARG roman_Γ end_ARG start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT over^ start_ARG roman_Γ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + over^ start_ARG roman_Γ end_ARG start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT over^ start_ARG roman_Γ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) = 2 over^ start_ARG roman_Γ end_ARG start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT over^ start_ARG roman_Γ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 2 ( italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP ⊗ start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 2 ) end_ARG end_RELOP + italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP ⊗ start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 ) end_ARG end_RELOP ) ×
(μ1γ0(1)0I(2)+μ2γ0(2)0I(1))=2(μ12I(1+2)+μ22I(1+2)+\displaystyle(\mu_{1}\stackrel{{\scriptstyle(1)}}{{\gamma^{0}}}\otimes% \stackrel{{\scriptstyle(2)}}{{I}}+\mu_{2}\stackrel{{\scriptstyle(2)}}{{\gamma^% {0}}}\otimes\stackrel{{\scriptstyle(1)}}{{I}})=2(\mu_{1}^{2}\stackrel{{% \scriptstyle(1+2)}}{{I}}+\mu_{2}^{2}\stackrel{{\scriptstyle(1+2)}}{{I}}+( italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP ⊗ start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 2 ) end_ARG end_RELOP + italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP ⊗ start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 ) end_ARG end_RELOP ) = 2 ( italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 + 2 ) end_ARG end_RELOP + italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 + 2 ) end_ARG end_RELOP +
μ1μ2(γ0(2)0I(1)γ0(1)0I(2)+γ0(1)0I(2)γ0(2)0I(1)))\displaystyle\mu_{1}\mu_{2}(\stackrel{{\scriptstyle(2)}}{{\gamma^{0}}}\otimes% \stackrel{{\scriptstyle(1)}}{{I}}\cdot\stackrel{{\scriptstyle(1)}}{{\gamma^{0}% }}\otimes\stackrel{{\scriptstyle(2)}}{{I}}+\stackrel{{\scriptstyle(1)}}{{% \gamma^{0}}}\otimes\stackrel{{\scriptstyle(2)}}{{I}}\cdot\stackrel{{% \scriptstyle(2)}}{{\gamma^{0}}}\otimes\stackrel{{\scriptstyle(1)}}{{I}}))italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP ⊗ start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 ) end_ARG end_RELOP ⋅ start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP ⊗ start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 2 ) end_ARG end_RELOP + start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP ⊗ start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 2 ) end_ARG end_RELOP ⋅ start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP ⊗ start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 ) end_ARG end_RELOP ) )
=2(μ12I(1+2)+μ22I(1+2)+2μ1μ2(γ0(1+2)0γ0(1+2)0))=2I(1+2),\displaystyle=2(\mu_{1}^{2}\stackrel{{\scriptstyle(1+2)}}{{I}}+\mu_{2}^{2}% \stackrel{{\scriptstyle(1+2)}}{{I}}+2\mu_{1}\mu_{2}(\stackrel{{\scriptstyle(1+% 2)}}{{\gamma^{0}}}\cdot\stackrel{{\scriptstyle(1+2)}}{{\gamma^{0}}}))=2% \stackrel{{\scriptstyle(1+2)}}{{I}},= 2 ( italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 + 2 ) end_ARG end_RELOP + italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 + 2 ) end_ARG end_RELOP + 2 italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 + 2 ) end_ARG end_RELOP ⋅ start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 + 2 ) end_ARG end_RELOP ) ) = 2 start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 + 2 ) end_ARG end_RELOP ,

where I(1+2)superscript𝐼12\stackrel{{\scriptstyle(1+2)}}{{I}}start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 + 2 ) end_ARG end_RELOP is the (16×16)1616(16\times 16)( 16 × 16 ) unit matrix.

Calculation of the anti-commutator [Γ^0,𝜸^]+subscriptsuperscript^Γ0^𝜸[{\hat{\Gamma}}^{0},{\hat{\hbox{\boldmath${\gamma}$}}}]_{+}[ over^ start_ARG roman_Γ end_ARG start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT , over^ start_ARG bold_italic_γ end_ARG ] start_POSTSUBSCRIPT + end_POSTSUBSCRIPT leads to the following results

Γ^0𝜸^+𝜸^Γ^0=(μ1(γ0(1)00000γ0(1)00000γ0(1)00000γ0(1)0)+μ2(γ0(2)00000γ0(2)00000γ0(2)00000γ0(2)0))×\displaystyle{\hat{\Gamma}}^{0}{\hat{\hbox{\boldmath${\gamma}$}}}+{\hat{\hbox{% \boldmath${\gamma}$}}}{\hat{\Gamma}}^{0}=(\mu_{1}\begin{pmatrix}\stackrel{{% \scriptstyle(1)}}{{\gamma^{0}}}&0&0&0\\ 0&\stackrel{{\scriptstyle(1)}}{{\gamma^{0}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(1)}}{{\gamma^{0}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(1)}}{{\gamma^{0}}}\\ \end{pmatrix}+\mu_{2}\begin{pmatrix}\stackrel{{\scriptstyle(2)}}{{\gamma^{0}}}% &0&0&0\\ 0&\stackrel{{\scriptstyle(2)}}{{\gamma^{0}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(2)}}{{\gamma^{0}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(2)}}{{\gamma^{0}}}\\ \end{pmatrix})\timesover^ start_ARG roman_Γ end_ARG start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT over^ start_ARG bold_italic_γ end_ARG + over^ start_ARG bold_italic_γ end_ARG over^ start_ARG roman_Γ end_ARG start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT = ( italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL end_ROW end_ARG ) + italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL end_ROW end_ARG ) ) ×
((𝜸(2)0000𝜸(2)0000𝜸(2)0000𝜸(2))(𝜸(1)0000𝜸(1)0000𝜸(1)0000𝜸(1)))+((𝜸(2)0000𝜸(2)0000𝜸(2)0000𝜸(2))(𝜸(1)0000𝜸(1)0000𝜸(1)0000𝜸(1)))×\displaystyle(\begin{pmatrix}\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${% \gamma}$}}}&0&0&0\\ 0&\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}}}\\ \end{pmatrix}-\begin{pmatrix}\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${% \gamma}$}}}&0&0&0\\ 0&\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}}}\\ \end{pmatrix})+(\begin{pmatrix}\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${% \gamma}$}}}&0&0&0\\ 0&\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}}}\\ \end{pmatrix}-\begin{pmatrix}\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${% \gamma}$}}}&0&0&0\\ 0&\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}}}\\ \end{pmatrix})\times( ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL end_ROW end_ARG ) - ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL end_ROW end_ARG ) ) + ( ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL end_ROW end_ARG ) - ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL end_ROW end_ARG ) ) ×
(μ1(γ0(1)00000γ0(1)00000γ0(1)00000γ0(1)0)+μ2(γ0(2)00000γ0(2)00000γ0(2)00000γ0(2)0))subscript𝜇1matrixsuperscriptsuperscript𝛾010000superscriptsuperscript𝛾010000superscriptsuperscript𝛾010000superscriptsuperscript𝛾01subscript𝜇2matrixsuperscriptsuperscript𝛾020000superscriptsuperscript𝛾020000superscriptsuperscript𝛾020000superscriptsuperscript𝛾02\displaystyle(\mu_{1}\begin{pmatrix}\stackrel{{\scriptstyle(1)}}{{\gamma^{0}}}% &0&0&0\\ 0&\stackrel{{\scriptstyle(1)}}{{\gamma^{0}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(1)}}{{\gamma^{0}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(1)}}{{\gamma^{0}}}\\ \end{pmatrix}+\mu_{2}\begin{pmatrix}\stackrel{{\scriptstyle(2)}}{{\gamma^{0}}}% &0&0&0\\ 0&\stackrel{{\scriptstyle(2)}}{{\gamma^{0}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(2)}}{{\gamma^{0}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(2)}}{{\gamma^{0}}}\\ \end{pmatrix})( italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL end_ROW end_ARG ) + italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL end_ROW end_ARG ) ) (72)
Γ^0𝜸^+𝜸^Γ^0=superscript^Γ0^𝜸^𝜸superscript^Γ0absent\displaystyle{\hat{\Gamma}}^{0}{\hat{\hbox{\boldmath${\gamma}$}}}+{\hat{\hbox{% \boldmath${\gamma}$}}}{\hat{\Gamma}}^{0}=over^ start_ARG roman_Γ end_ARG start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT over^ start_ARG bold_italic_γ end_ARG + over^ start_ARG bold_italic_γ end_ARG over^ start_ARG roman_Γ end_ARG start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT =
μ1((γ0(1)0𝜸(2)0000γ0(1)0𝜸(2)0000γ0(1)0𝜸(2)0000γ0(1)0𝜸(2))(γ0(1)0𝜸(1)0000γ0(1)0𝜸(1)0000γ0(1)0𝜸(1)0000γ0(1)0𝜸(1)+(𝜸(2)γ0(1)00000𝜸(2)γ0(1)00000𝜸(2)γ0(1)00000𝜸(2)γ0(1)0))\displaystyle\mu_{1}(\begin{pmatrix}\stackrel{{\scriptstyle(1)}}{{\gamma^{0}}}% \stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}}}&0&0&0\\ 0&\stackrel{{\scriptstyle(1)}}{{\gamma^{0}}}\stackrel{{\scriptstyle(2)}}{{% \hbox{\boldmath${\gamma}$}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(1)}}{{\gamma^{0}}}\stackrel{{\scriptstyle(2)}}{{% \hbox{\boldmath${\gamma}$}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(1)}}{{\gamma^{0}}}\stackrel{{\scriptstyle(2)}}{{% \hbox{\boldmath${\gamma}$}}}\\ \end{pmatrix}-\begin{pmatrix}\stackrel{{\scriptstyle(1)}}{{\gamma^{0}}}% \stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}}}&0&0&0\\ 0&\stackrel{{\scriptstyle(1)}}{{\gamma^{0}}}\stackrel{{\scriptstyle(1)}}{{% \hbox{\boldmath${\gamma}$}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(1)}}{{\gamma^{0}}}\stackrel{{\scriptstyle(1)}}{{% \hbox{\boldmath${\gamma}$}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(1)}}{{\gamma^{0}}}\stackrel{{\scriptstyle(1}}{{% \hbox{\boldmath${\gamma}$}}}\\ \end{pmatrix}+\begin{pmatrix}\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${% \gamma}$}}}\stackrel{{\scriptstyle(1)}}{{\gamma^{0}}}&0&0&0\\ 0&\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}}}\stackrel{{% \scriptstyle(1)}}{{\gamma^{0}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}}}\stackrel{{% \scriptstyle(1)}}{{\gamma^{0}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}}}\stackrel{{% \scriptstyle(1)}}{{\gamma^{0}}}\\ \end{pmatrix})italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL end_ROW end_ARG ) - ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 end_ARG end_RELOP end_CELL end_ROW end_ARG ) + ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL end_ROW end_ARG ) )
(𝜸(1)γ0(1)00000𝜸(1)γ0(1)00000𝜸(1)γ0(1)00000𝜸(1)γ0(1)0))+μ2((γ0(2)0𝜸(2)0000γ0(2)0𝜸(2)0000γ0(2)0𝜸(2)0000γ0(2)0𝜸(2))(γ0(2)0𝜸(1)0000γ0(2)0𝜸(1)0000γ0(2)0𝜸(1)0000γ0(2)0𝜸(1)\displaystyle-\begin{pmatrix}\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${% \gamma}$}}}\stackrel{{\scriptstyle(1)}}{{\gamma^{0}}}&0&0&0\\ 0&\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}}}\stackrel{{% \scriptstyle(1)}}{{\gamma^{0}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}}}\stackrel{{% \scriptstyle(1)}}{{\gamma^{0}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}}}\stackrel{{% \scriptstyle(1)}}{{\gamma^{0}}}\\ \end{pmatrix})+\mu_{2}(\begin{pmatrix}\stackrel{{\scriptstyle(2)}}{{\gamma^{0}% }}\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}}}&0&0&0\\ 0&\stackrel{{\scriptstyle(2)}}{{\gamma^{0}}}\stackrel{{\scriptstyle(2)}}{{% \hbox{\boldmath${\gamma}$}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(2)}}{{\gamma^{0}}}\stackrel{{\scriptstyle(2)}}{{% \hbox{\boldmath${\gamma}$}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(2)}}{{\gamma^{0}}}\stackrel{{\scriptstyle(2)}}{{% \hbox{\boldmath${\gamma}$}}}\\ \end{pmatrix}-\begin{pmatrix}\stackrel{{\scriptstyle(2)}}{{\gamma^{0}}}% \stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}}}&0&0&0\\ 0&\stackrel{{\scriptstyle(2)}}{{\gamma^{0}}}\stackrel{{\scriptstyle(1)}}{{% \hbox{\boldmath${\gamma}$}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(2)}}{{\gamma^{0}}}\stackrel{{\scriptstyle(1)}}{{% \hbox{\boldmath${\gamma}$}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(2)}}{{\gamma^{0}}}\stackrel{{\scriptstyle(1}}{{% \hbox{\boldmath${\gamma}$}}}\\ \end{pmatrix}- ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL end_ROW end_ARG ) ) + italic_μ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL end_ROW end_ARG ) - ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 end_ARG end_RELOP end_CELL end_ROW end_ARG )
+(𝜸(2)γ0(2)00000𝜸(2)γ0(2)00000𝜸(2)γ0(2)00000𝜸(2)γ0(2)0))(𝜸(1)γ0(2)00000𝜸(1)γ0(2)00000𝜸(1)γ02)00000𝜸(1)γ0(2)0)).\displaystyle+\begin{pmatrix}\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${% \gamma}$}}}\stackrel{{\scriptstyle(2)}}{{\gamma^{0}}}&0&0&0\\ 0&\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}}}\stackrel{{% \scriptstyle(2)}}{{\gamma^{0}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}}}\stackrel{{% \scriptstyle(2)}}{{\gamma^{0}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}}}\stackrel{{% \scriptstyle(2)}}{{\gamma^{0}}}\\ \end{pmatrix})-\begin{pmatrix}\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${% \gamma}$}}}\stackrel{{\scriptstyle(2)}}{{\gamma^{0}}}&0&0&0\\ 0&\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}}}\stackrel{{% \scriptstyle(2)}}{{\gamma^{0}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}}}\stackrel{{% \scriptstyle 2)}}{{\gamma^{0}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}}}\stackrel{{% \scriptstyle(2)}}{{\gamma^{0}}}\\ \end{pmatrix}).+ ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 2 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL end_ROW end_ARG ) ) - ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ end_ARG start_ARG ( 1 ) end_ARG end_RELOP start_RELOP SUPERSCRIPTOP start_ARG italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL end_ROW end_ARG ) ) .

Then, we finally derive

Γ^0𝜸^+𝜸^Γ^0=0.superscript^Γ0^𝜸^𝜸superscript^Γ00\displaystyle{\hat{\Gamma}}^{0}{\hat{\hbox{\boldmath${\gamma}$}}}+{\hat{\hbox{% \boldmath${\gamma}$}}}{\hat{\Gamma}}^{0}=0.over^ start_ARG roman_Γ end_ARG start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT over^ start_ARG bold_italic_γ end_ARG + over^ start_ARG bold_italic_γ end_ARG over^ start_ARG roman_Γ end_ARG start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT = 0 . (74)

Let us go to the anti-commutator [𝜸a^,𝜸b^]+subscript^superscript𝜸𝑎^superscript𝜸𝑏[{\hat{\hbox{\boldmath${\gamma}$}^{a}},}{\hat{\hbox{\boldmath${\gamma}$}^{b}}}% ]_{+}[ over^ start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG , over^ start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG ] start_POSTSUBSCRIPT + end_POSTSUBSCRIPT

(𝜸a^𝜸b^+𝜸b^𝜸a^)=(𝜸a(2)aI(1)𝜸a(1)aI(2))(𝜸b(2)bI(1)𝜸b(1)bI(2))\displaystyle({\hat{\hbox{\boldmath${\gamma}$}^{a}}}{\hat{\hbox{\boldmath${% \gamma}$}^{b}}}+{\hat{\hbox{\boldmath${\gamma}$}^{b}}}{\hat{\hbox{\boldmath${% \gamma}$}^{a}}})=(\stackrel{{\scriptstyle(2)}}{{{\hbox{\boldmath${\gamma}$}^{a% }}}}\otimes\stackrel{{\scriptstyle(1)}}{{I}}-\stackrel{{\scriptstyle(1)}}{{% \hbox{\boldmath${\gamma}$}^{a}}}\otimes\stackrel{{\scriptstyle(2)}}{{I}})(% \stackrel{{\scriptstyle(2)}}{{{\hbox{\boldmath${\gamma}$}^{b}}}}\otimes% \stackrel{{\scriptstyle(1)}}{{I}}-\stackrel{{\scriptstyle(1)}}{{\hbox{% \boldmath${\gamma}$}^{b}}}\otimes\stackrel{{\scriptstyle(2)}}{{I}})( over^ start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG over^ start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG + over^ start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG over^ start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG ) = ( start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP ⊗ start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 ) end_ARG end_RELOP - start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP ⊗ start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 2 ) end_ARG end_RELOP ) ( start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP ⊗ start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 ) end_ARG end_RELOP - start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP ⊗ start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 2 ) end_ARG end_RELOP )
+\displaystyle++ (𝜸b(2)bI(1)𝜸b(1)bI(2))(𝜸a(2)aI(1)𝜸a(1)aI(2)).\displaystyle(\stackrel{{\scriptstyle(2)}}{{{\hbox{\boldmath${\gamma}$}^{b}}}}% \otimes\stackrel{{\scriptstyle(1)}}{{I}}-\stackrel{{\scriptstyle(1)}}{{\hbox{% \boldmath${\gamma}$}^{b}}}\otimes\stackrel{{\scriptstyle(2)}}{{I}})(\stackrel{% {\scriptstyle(2)}}{{{\hbox{\boldmath${\gamma}$}^{a}}}}\otimes\stackrel{{% \scriptstyle(1)}}{{I}}-\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$% }^{a}}}\otimes\stackrel{{\scriptstyle(2)}}{{I}}).( start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP ⊗ start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 ) end_ARG end_RELOP - start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP ⊗ start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 2 ) end_ARG end_RELOP ) ( start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP ⊗ start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 ) end_ARG end_RELOP - start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP ⊗ start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 2 ) end_ARG end_RELOP ) . (75)

Substituting the explicit form of 𝜸a^^superscript𝜸𝑎{\hat{\hbox{\boldmath${\gamma}$}^{a}}}over^ start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG we get

(𝜸a^𝜸b^+𝜸b^𝜸a^)=((𝜸a(2)a0000𝜸a(2)a0000𝜸a(2)a0000𝜸a(2)a)(𝜸a(1)a0000𝜸a(1)a0000𝜸a(1)a0000𝜸a(1)a))×((𝜸b(2)b0000𝜸b(2)b0000𝜸b(2)b0000𝜸b(2)b)(𝜸b(1)b0000𝜸b(1)b0000𝜸b(1)b0000𝜸b(1)b))^superscript𝜸𝑎^superscript𝜸𝑏^superscript𝜸𝑏^superscript𝜸𝑎matrixsuperscriptsuperscript𝜸𝑎20000superscriptsuperscript𝜸𝑎20000superscriptsuperscript𝜸𝑎20000superscriptsuperscript𝜸𝑎2matrixsuperscriptsuperscript𝜸𝑎10000superscriptsuperscript𝜸𝑎10000superscriptsuperscript𝜸𝑎10000superscriptsuperscript𝜸𝑎1matrixsuperscriptsuperscript𝜸𝑏20000superscriptsuperscript𝜸𝑏20000superscriptsuperscript𝜸𝑏20000superscriptsuperscript𝜸𝑏2matrixsuperscriptsuperscript𝜸𝑏10000superscriptsuperscript𝜸𝑏10000superscriptsuperscript𝜸𝑏10000superscriptsuperscript𝜸𝑏1\displaystyle({\hat{\hbox{\boldmath${\gamma}$}^{a}}}{\hat{\hbox{\boldmath${% \gamma}$}^{b}}}+{\hat{\hbox{\boldmath${\gamma}$}^{b}}}{\hat{\hbox{\boldmath${% \gamma}$}^{a}}})=(\begin{pmatrix}\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath% ${\gamma}$}^{a}}}&0&0&0\\ 0&\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}^{a}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}^{a}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}^{a}}}\\ \end{pmatrix}-\begin{pmatrix}\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${% \gamma}$}^{a}}}&0&0&0\\ 0&\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}^{a}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}^{a}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}^{a}}}\\ \end{pmatrix})\times(\begin{pmatrix}\stackrel{{\scriptstyle(2)}}{{\hbox{% \boldmath${\gamma}$}^{b}}}&0&0&0\\ 0&\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}^{b}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}^{b}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}^{b}}}\\ \end{pmatrix}-\begin{pmatrix}\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${% \gamma}$}^{b}}}&0&0&0\\ 0&\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}^{b}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}^{b}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}^{b}}}\\ \end{pmatrix})( over^ start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG over^ start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG + over^ start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG over^ start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG ) = ( ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL end_ROW end_ARG ) - ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL end_ROW end_ARG ) ) × ( ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL end_ROW end_ARG ) - ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL end_ROW end_ARG ) )
+((𝜸b(2)b0000𝜸b(2)b0000𝜸b(2)b0000𝜸b(2)b)(𝜸b(1)b0000𝜸b(1)b0000𝜸b(1)b0000𝜸b(1)b))×((𝜸a(2)a0000𝜸a(2)a0000𝜸a(2)a0000𝜸a(2)a)(𝜸a(1)a0000𝜸a(1)a0000𝜸a(1)a0000𝜸a(1)a)).matrixsuperscriptsuperscript𝜸𝑏20000superscriptsuperscript𝜸𝑏20000superscriptsuperscript𝜸𝑏20000superscriptsuperscript𝜸𝑏2matrixsuperscriptsuperscript𝜸𝑏10000superscriptsuperscript𝜸𝑏10000superscriptsuperscript𝜸𝑏10000superscriptsuperscript𝜸𝑏1matrixsuperscriptsuperscript𝜸𝑎20000superscriptsuperscript𝜸𝑎20000superscriptsuperscript𝜸𝑎20000superscriptsuperscript𝜸𝑎2matrixsuperscriptsuperscript𝜸𝑎10000superscriptsuperscript𝜸𝑎10000superscriptsuperscript𝜸𝑎10000superscriptsuperscript𝜸𝑎1\displaystyle+(\begin{pmatrix}\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${% \gamma}$}^{b}}}&0&0&0\\ 0&\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}^{b}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}^{b}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}^{b}}}\\ \end{pmatrix}-\begin{pmatrix}\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${% \gamma}$}^{b}}}&0&0&0\\ 0&\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}^{b}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}^{b}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}^{b}}}\\ \end{pmatrix})\times(\begin{pmatrix}\stackrel{{\scriptstyle(2)}}{{\hbox{% \boldmath${\gamma}$}^{a}}}&0&0&0\\ 0&\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}^{a}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}^{a}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(2)}}{{\hbox{\boldmath${\gamma}$}^{a}}}\\ \end{pmatrix}-\begin{pmatrix}\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${% \gamma}$}^{a}}}&0&0&0\\ 0&\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}^{a}}}&0&0\\ 0&0&\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}^{a}}}&0\\ 0&0&0&\stackrel{{\scriptstyle(1)}}{{\hbox{\boldmath${\gamma}$}^{a}}}\\ \end{pmatrix}).+ ( ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL end_ROW end_ARG ) - ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL end_ROW end_ARG ) ) × ( ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG start_ARG ( 2 ) end_ARG end_RELOP end_CELL end_ROW end_ARG ) - ( start_ARG start_ROW start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL start_RELOP SUPERSCRIPTOP start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 ) end_ARG end_RELOP end_CELL end_ROW end_ARG ) ) .

Making obvious transformation we derive

(𝜸a^𝜸b^+𝜸b^𝜸a^)=4gabI(1+2)/2+(σ1aσ2b+σ2bσ1a+σ2aσ1b+σ1bσ2a)I(1+2)/2^superscript𝜸𝑎^superscript𝜸𝑏^superscript𝜸𝑏^superscript𝜸𝑎4superscript𝑔𝑎𝑏superscript𝐼122superscriptsubscript𝜎1𝑎superscriptsubscript𝜎2𝑏superscriptsubscript𝜎2𝑏superscriptsubscript𝜎1𝑎superscriptsubscript𝜎2𝑎superscriptsubscript𝜎1𝑏superscriptsubscript𝜎1𝑏superscriptsubscript𝜎2𝑎superscript𝐼122absent\displaystyle({\hat{\hbox{\boldmath${\gamma}$}^{a}}}{\hat{\hbox{\boldmath${% \gamma}$}^{b}}}+{\hat{\hbox{\boldmath${\gamma}$}^{b}}}{\hat{\hbox{\boldmath${% \gamma}$}^{a}}})=4g^{ab}\stackrel{{\scriptstyle(1+2)/2}}{{I}}+(\sigma_{1}^{a}% \sigma_{2}^{b}+\sigma_{2}^{b}\sigma_{1}^{a}+\sigma_{2}^{a}\sigma_{1}^{b}+% \sigma_{1}^{b}\sigma_{2}^{a})\stackrel{{\scriptstyle(1+2)/2}}{{I}}( over^ start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG over^ start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG + over^ start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT end_ARG over^ start_ARG bold_italic_γ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT end_ARG ) = 4 italic_g start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 + 2 ) / 2 end_ARG end_RELOP + ( italic_σ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT italic_σ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT + italic_σ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT italic_σ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT + italic_σ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT italic_σ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT + italic_σ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT italic_σ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT ) start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 + 2 ) / 2 end_ARG end_RELOP
=4δabI(1+2)+(σ1aσ2b+σ2bσ1a+σ2aσ1b+σ1bσ2a)I(1+2)/2.absent4superscript𝛿𝑎𝑏superscript𝐼12superscriptsubscript𝜎1𝑎superscriptsubscript𝜎2𝑏superscriptsubscript𝜎2𝑏superscriptsubscript𝜎1𝑎superscriptsubscript𝜎2𝑎superscriptsubscript𝜎1𝑏superscriptsubscript𝜎1𝑏superscriptsubscript𝜎2𝑎superscript𝐼122absent\displaystyle=-4\delta^{ab}\stackrel{{\scriptstyle(1+2)}}{{I}}+(\sigma_{1}^{a}% \sigma_{2}^{b}+\sigma_{2}^{b}\sigma_{1}^{a}+\sigma_{2}^{a}\sigma_{1}^{b}+% \sigma_{1}^{b}\sigma_{2}^{a})\stackrel{{\scriptstyle(1+2)/2}}{{I}}.= - 4 italic_δ start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 + 2 ) end_ARG end_RELOP + ( italic_σ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT italic_σ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT + italic_σ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT italic_σ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT + italic_σ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT italic_σ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT + italic_σ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT italic_σ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT ) start_RELOP SUPERSCRIPTOP start_ARG italic_I end_ARG start_ARG ( 1 + 2 ) / 2 end_ARG end_RELOP .

References

  • (1) G.Breit, The Effect of Retardation on the Interaction of Two Electrons, Phys. Rev. 34 553 (1929).
  • (2) G.Breit, The Fine Structure of HE as a Test of the Spin Interactions of Two Electrons, Phys. Rev.. 36 383 ( 1930) ; Dirac’s Equation and the Spin-Spin Interactions of Two Electrons, Phys. Rev.. 39 553 (1932).
  • (3) I.T.Todorov,Quasipotential Equation Corresponding to the Relativistic Eikonal Approximation, Phys.Rev. D3 2351 (1971) .
  • (4) I.T.Todorov, Dynamics of Relativistic Point Particles as a Problem with Constraints, Dubna Joint Institute for Nuclear Researches, Reprint, E2-10125, (1976).
  • (5) Ph. Droz-Vincent, Hamiltonian systems of relativistic particles, Reports on Math Phys. B 8 79(1975) .
  • (6) M. Kalb and P. van Alstine, Invariant Singular Actions for the two body Relativistic Problem: A Hamiltonian Formulation, Yale Reports, C00-3075-146; C00-3075-156 (1976).
  • (7) I.T.Todorov, 1978, Sur la quantification d’ un systeme mecanique avec des contraintes de deuxieme classe , Annals of the Institute of H. Poincare A28, 207.
  • (8) A.Komar, Constraint formalism of classical mechanics, Phys. Rev.D 18 1881 (1978); Interacting relativistic particles, Phys. Rev.D 18 1887 (1978).
  • (9) R. Arens, Classical Lorentz Invariant Particles., J.Math.Phys. 12, 2415 (1971).
  • (10) A.A. Kirillov. Element of the Theory of Representations, Springer-Verlag Berlin Heidelberg,1976
  • (11) D. Dominici, J. Gomis, G. Longhi, A Lagrangian for two interacting relativistic particles , Nuovo Cim. , B 48, 152 (1978); Nuovo Cim.A Lagrangian for two interacting relativistic particles: Canonical formulation , A 48, 257 (1978).
  • (12) D. Dominici, J. Gomis, G. Longhi, A Lagrangian approach to the no-interaction theorem Nuovo Cim., bf A 56, 263 (1980).
  • (13) Ph.Droz-Vincent, Two-Body Mass-Shell Constraints in a Constant Magnetic Field (Neutral Case), Int. J. Theor. Phys., 39, 389 ( 2000).
  • (14) J. Bijtebier, A cluster-separable Born approximation for the 3D reduction of the three-fermion Bethe-Salpeter equation , J. Phys., G 26, 871 (2000).
  • (15) Ph.Droz-Vincent, Relativistic Two-Boson System in Presence of Electromagnetic Plane Wave, Int.J. Theor. Phys. , 48, 2177 (2009).
  • (16) H.Sazdjian,  it Relativistic wave equations for the dynamics of two interacting particles Phys. Rev. D33 3401 (1986) .
  • (17) H.Sazdjian, Light fermion bound states in two-particle relativistic quantum mechanics, Phys. Rev. D33 3435 (1986) .
  • (18) J.Mourad, H.Sazdjian, The two-fermion relativistic wave equations of Constraint Theory in the Pauli-Schrödinger form, J.Math.Phys. 35 6379 ( 1994).
  • (19) H.Sazdjian, N-body bound state relativistic wave equations , Ann. of Phys.(Elsevier) 191 52 (1989).
  • (20) Crater, H.W., VanAlstine, P., Two-body Dirac equations. Annals of Physics. 148 57 (1983)/
  • (21) Crater, H.W., VanAlstine, P., Two-body Dirac equations for particles interacting through world scalar and vector potentials, Phys. Rev. D36, 3007 (1987).
  • (22) H.W.Crater, P. Van Alstine, Extension of two-body Dirac equations to general covariant interactions through a hyperbolic transformation , J. Math. Phys. 31 1998 ( 1990).
  • (23) C.Y.Wong, H.W.Crater, RelativisticN-body problem in a separable two-body basis, Phys. Rev. C63 044907 (2001).
  • (24) B.Liu, H.W.Crater, Two-body Dirac equations for nucleon-nucleon scattering, Phys. Rev. C67 024001 (2003).
  • (25) Crater, H.W., Wong, C.Y, van Alstine, P., Two Gamma Quarkonium and Positronium Decays with Two-Body Dirac Equations of Constraint Dynamics, Phys. Rev. D74 054028 (2006).
  • (26) Crater, H.W., Schiermeyer, Relativistic two-body calculation of bb¯𝑏¯𝑏b\bar{b}italic_b over¯ start_ARG italic_b end_ARG-mesons radiative decays, Phys. Rev. D82 094020 (2010) .
  • (27) A.O.Barut, N. Unal, A new approach to bound state quantum electrodynamics. 1. Theory, Physica 142A 467 (1987) .
  • (28) A.O.Barut , The covariant many-body problem in quantumelectrodynamics, J.Math.Phys 32 1091 (1991).
  • (29) A.A.Logunov and A. N.Tavkhelidze, Quasi-optical approach in quantum field theory, Nuovo Cimento 29 380 (1963).
  • (30) Aguilar, M.; et al. , (AMS Collaboration), First Result from the Alpha Magnetic Spectrometer on the International Space Station: Precision Measurement of the Positron Fraction in Primary Cosmic Rays of 0.5-350 GeV , Phy.l Rev. Lett. 110, 141102 (2013) .
  • (31) J.Perez-Rios, S. T.Love, Searching for light dark matter through Positronium decay, Euro. Phys. Journal D72 , 44 ( 2018).
  • (32) Jung, D-W., Lee,K-Y., Yu, C., Positronium decays with a dark Z𝑍Zitalic_Z and fermionic dark matter , Phys. Rev. D105 095023 (2022).
  • (33) G. Bertone, D. Hooper and J. Silk, Particle Dark Matter: Evidence, Candidates and Constraints, Phys. Rep. 405 2005
  • (34) M. Geller,a S. Iwamoto,b,c G. Lee,b,d,e Y. Shadmib and O. Telemd, Dark quarkonium formation in the early universe, JHEP 06 , 135 (2018).
  • (35) B.L.Brown, Astrophys. Journal, Part 2 - Letters to the Editor, A possible galactic positron annihilation medium - Neutral hydrogen, 292, L68 (1985) .
  • (36) R.W.Bussard , R.Ramaty, R. J.Drachman, The annihilation of galactic positrons, Astrophys, Journal, 228, 928 (1979).
  • (37) R. C.Duncan, C.Thompson, Formation of Very Strongly Magnetized Neutron Stars, Astrophys Journal Lett, 392, L9 (1992).
  • (38) A.K.Harding, L. Dong,Physics of Strongly Magnetized Neutron Stars, Reports on Progress in Physics 9 2631 ( 2006) .
  • (39) K.Fukushima, D. E.Kharzeev , H.Warringa , Testing the Chiral Magnetic and Chiral Vortical Effects in Heavy Ion Collisions, Phys.Rev.Lett.,104 , 212001 (2010).
  • (40) V.Voronyuk , V.Toneev, W.Cassing, E.Bratkovskaya, V.Konchakovski et al., Electromagnetic field evolution in relativistic heavy-ion collisions, Phys.Rev.C, C83,054911 (2012).
  • (41) K.Tuchin, Electromagnetic field evolution in relativistic heavy-ion collisions, Adv.High Energy Phys. 2103, 490495 (2013) .
  • (42) H.A.Bethe, E.E Salpeter , Quantum mechanics of one- and two-electron atoms (Dover ed.). Mineola, N.Y.: Dover Publications, 2008.
  • (43) N. Nakanishi, Progress of Theoretical Physics Supplement. Oxford University Press (OUP), General Survey of the Theory of the Bethe-Salpeter Equation, 43, 1 ( 1969) .
  • (44) P.A.M. Dirac, Lectures on Quantum Mechanics,Yeshiva University, New York, 1964.
  • (45) J.L.Anderson, P.G. Bergmann,Constraints in Covariant Field Theories,, Phys. Rev. 83, 1018 (1951) .
  • (46) L.D.Landau and E.M.Lifshitz, Quantum Electrodynamics. Second edition. Course of Theoretical Physics, Volume 4, p. 575. Pergamon Press, 1982.
  • (47) L.D.Landau and E.M.Lifshitz, it Quantum Mechanics. Non-relativistic Theory. Second edition, revised and enlarged. Course of Theoretical Physics, Volume 3, p. 426. Pergamon Press, 1965.
  • (48) I.S.Gradshteyn and I.M.Ryzhik. Table of Integrals, Series, and Products. Seventh Edition. Elsevier, 2007.
  • (49) R.J.Elliott, R.Loudon, Theory of the absorption edge in semiconductors in a high magnetic field, J. Phys. Chem. Solids, 15, 196 (1960).
  • (50) Alford, J., Strickland, M., Charmonia and bottomonia in a magnetic field, Phys. Rev.D88, 105017 (2013) .
  • (51) D.Green. Lectures in particle physics, p.102., World scientific, 1994.