arXiv:2011.06401v1 [math-ph] 12 Nov 2020 ntetredmninld itrsaetm AdS space-time Sitter de three-dimensional new the of set in complete a construct constru we is Also, equation method. Dirac integration the to solutions exact variab of of basis separation the admit not does develope which method space the homogeneous of structure general The var it. of supplement separation of method well-known the from differe differs partial method linear Thi of space. method homogeneous integration a equ noncommutative of be transformations to of shown group Lie is a metric on invariant an with equation Dirac The ahmtc ujc lsicto 00 54,1B8 58J70 17B08, 35Q41, 2010: Classification Subject Mathematics method homogen integration; noncommutative equation; Dirac Keywords: 31.15.xh 03.65.Pm, 02.20.Qs, numbers: PACS i method. which integration functions, elementary of terms in found are obtained ∗ † [email protected] [email protected] edvlpannomttv nerto ehdfrteDira the for method integration noncommutative a develop We ocmuaieitgaino h ia qainin equation Dirac the of integration Noncommutative os oyehi nvriy ei v. 0 os,Russi Tomsk, 30, ave., Lenin University, Polytechnic Tomsk eateto hoeia hsc,TmkSaeUniversity State Tomsk Physics, Theoretical of Department nvriyNvsbraaS.1 os,Rsi,634050 Russia, Tomsk, 1, Sq. Novosobornaya University oooony q ,Tmk usa 300and 634050 Russia, Tomsk, 1, Sq. Novosobornaya eateto hoeia hsc,TmkState Tomsk Physics, Theoretical of Department ooeeu spaces homogeneous .V Shapovalov V. A. .I Breev I. A. Abstract 1 3 sn h ehddvlpd h solutions The developed. method the using ∗ † tdepiil ytenoncommutative the by explicitly cted hrceitco h noncommutative the of characteristic s ossae;idcdrpeettos orbit representations; induced spaces; eous ta qain nLegop.This groups. Lie on equations ntial silsrtdwt neapeo a of example an with illustrated is d xc ouin oteDrcequation Dirac the to solutions exact e nteDrceuto.However, equation. Dirac the in les alsadt oeetn a often can extent some to and iables losu oeetvl pl the apply effectively to us allows s qaini ooeeu spaces. homogeneous in equation c vln oasse fequations of system a to ivalent ,634034 a, I. INTRODUCTION

Exact solutions of the relativistic wave equations in strong gravitational and electromag- netic fields are the basis for studying quantum effects in the framework of quantum field theory in curved space-time (see, e.g. [1–6]). A construction of the complete set of exact solutions to these equations in many cases is associated with the presence of integrals of motion. For example, to separate the variables in a wave equation, it is necessary to have dim M 1 commuting integrals, where M is the space of independent variables. In this − paper, by integrability of the wave equation we mean an explicit possibility of reducing the original equation to a system of ordinary differential equations, the solution of which provides a complete set of solutions to the original wave equation. The best-known technique for such a reduction is based on the method of separation of variables (SoV) (various aspects of the SoV method can be found, e.g., in [7–9]). There is a broad scope of research dealing with separation of variables in relativistic quantum wave equations, mainly for the Klein-Gordon and Dirac equations, and with classification of external fields admitting SoV in these equations (see, e.g., [10] and references therein). This motivates the development of methods for the exact integration of wave equations other than SoV that can give some new possibilities in relativistic quantum theory. In this regard, we focus on homogeneous spaces as geometric objects with high symmetry. We also note that most of the physically interesting problems and effects are associated with gravitational fields possessing symmetries. Mathematically, these symmetries indicate the presence of various groups of transformations that leave invariant the gravitational field. Representing the space-time as a homogeneous space with a group-invariant metric, we can consider a large class of gravitational fields and cosmological models [11, 12] with rich symmetries, and the corresponding relativistic equations in these fields have integrals of motion. We note that the relativistic wave equations on a homogeneous space may not allow sep- aration of variables. The matter is that in accordance with the theorem of Refs. [13, 14], for the separation of variables in the wave equation in an appropriate coordinate system the equation should admit a complete set of mutually commuting symmetry operators (integrals of motion, details can be found in [13, 14], see also [15]). Therefore, the problem arises of constructing exact solutions to the wave equation in the case when it has symmetry oper-

2 ators, but they do not form a complete set and separation of variables can not be carried out. We consider the noncommutative integration method (NCIM) based on noncommu- tative algebras of symmetry operators admitted by the equation [16–20]. This method can be thought as a generalization of the method of SoV. A reduction of the wave equation to a system of ODEs according to the NCIM (we use the term noncommutative reduction) can be carried out in a way that is substantially different from the method of separation of variables. We note that the method of noncommutative integration has shown its effectiveness in constructing bases of exact solutions to the Klein-Gordon and Dirac equations in some spaces with invariance groups. For instance, the NCIM was applied to the Klein-Gordon equation in homogeneous spaces with an invariant metric in [19, 20]. The polarization vacuum effect of a scalar field in a homogeneous space was studied using NCIM in [19–21]. The noncommutative reduction of the Dirac equation to a system of ordinary differential equations in the Riemannian and pseudo-Riemannian spaces with a nontrivial group of motions was considered in [22–27]. In Refs. [28, 29] the NCIM was applied to the Dirac equation in the four-dimensional flat space and in the de Sitter space. The Dirac equation on Lie groups that can be a special case of homogeneous spaces with a trivial isotropy subgroup, was explored in terms of the NCIM in Refs. [30, 31]. It may also be worth noting that the application of the NCIM to the Dirac equation can give a new class of its exact solutions, different from the solutions obtained by SoV. In cases where the Dirac equation does not admit separation of variables, the NCIM provides an uncontested option for constructing complete sets of solutions. The physical meaning of the solutions obtained by this method depends on the specifics of the problem being solved and requires special research in each case. In the present work we consider non-commutative symmetries of the Dirac equation in homogeneous spaces. We also develop the method of noncommutative integration of the Dirac equation in homogeneous spaces. Using the group-theoretic approach, we reduce the Dirac equation on the homogeneous space to such a system of equations on the transforma- tion group that lets us to apply the noncommutative reduction and construct exact solutions of the Dirac equation. In this paper, for the first time, we explicitly take into account the identities for generators of the transformation group in the problem of noncommutative

3 reduction for the Dirac equation. The work is organized as follows. In Section II we briefly introduce basic concepts and notations from the theory of homogeneous spaces [32–34], to be used later. A construction of invariant differential operator with matrix coefficients on a homogeneous space is introduced in Section III following Refs.[35, 36]. Also in this section, we show the connection between generators of the representation of a Lie group on a homogeneous space and the other representation induced by representation of a subgroup, whose action on a homogeneous space has a stationary point. In the next Section IV, we introduce a special irreducible representation of the Lie algebra of the Lie group of transformations of a homogeneous space using the Kirillov orbit method [37], that is necessary for noncommutative reduction. In Section V we present the Dirac equation in a homogeneous space with an invariant metric in terms of an invariant matrix operator of the first order. The spinor connection and symmetry operators of the Dirac equation are shown to define isotropy representation in a spinor space. Generators of the spinor representation are found explicitly. We also introduce a system of differential equations on the Lie group of transformations of a homogeneous space, which is equivalent to the original Dirac equation in a homogeneous space. Then, in Section VI, we present a noncommutative reduction of the Dirac equation on a homogeneous space, using the irreducible λ-representation introduced in section IV and functional relations between symmetry operators (identities) for the Dirac equation. In Section VII, we consider a homogeneous space with an invariant metric that does not admit separation of variables for the Klein–Gordon and Dirac equations. In this case a com- plete set of exact solutions of the Dirac equation is constructed using the noncommutative reduction (Section VI).

The next Section VIII is devoted to the Dirac equation in the (2+1) anti-de Sitter AdS3 - dimensional space. In this homogeneous space there are identities between the generators of the representation of the group SO(1, 3) that is taken into account when the noncommutative

reduction is applied. The Dirac equation admits separation of variables in the AdS3 space, but the separable solutions are expressed through the special functions and have a complex form. The NCIM, being applied to this problem, results in the other complete set of exact

solutions to the Dirac equation in the AdS3 space and these solutions are presented in terms

4 of elementary functions. In Section IX we give our conclusion remarks.

II. INVARIANT METRIC ON A HOMOGENEOUS SPACE

This section introduces some basic concepts and notations of the homogeneous space theory with an invariant metric. Let G be a simply connected real Lie group with a Lie algebra g, M be a homogeneous space with right action of the group G, (x, g) R x = xg M for x M, g G. For any → g ∈ ∈ ∈ x M there exists an isotropic subgroup H G. Denote by H = H a closed stabilizer of ∈ x ∈ x0 a point x M, and let h be a Lie algebra of H. The homogeneous space M is diffeomorphic 0 ∈ to a quotient manifold G/H of right cosets [Hg] of the Lie group G by H. A transformation group G can be regarded as a principal bundle (G,π,M,H) with a structure group H, a base M, and a canonical projection π : G M, π(e) = x , where → 0 e is the identity element of G. An arbitrary point g G can be represented uniquely as ∈ g = hs(x), where x = x g = π(g) M, h H, and s : M G is a local and smooth 0 ∈ ∈ → section of G, π s = id. ◦ Differential of the canonical projection π : T G T M is a surjective map that allows ∗ g → π(g) any tangent vector τ T M on a homogeneous space to be represented as π ζ, where ∈ x ∗ ζ T G is a tangent vector on G. ∈ g In turn, a linear space of the Lie algebra g T G is decomposed into a direct sum of ≃ e subspaces g = h m, where m = π (g) T M is a complement to h = ker π (e), i.e. ⊕ ∗ ≃ x0 ∗ X = Xm + Xh holds for any X g, where Xm = π X, Xh h. ∈ ∗ ∈ We introduce an invariant metric on the homogeneous space M. Let , m be a non- h· ·i degenerate Ad(H) - invariant scalar product on the subspace m,

[X,Y ]m,Zm m + Ym, [X,Z]m m =0, X h, Y,Z g. (2.1) h i h i ∈ ∈

By action of a Lie group G with right shifts on the homogeneous space M, we define the inner product throughout the space M as

u, v = (R −1 ) u, (R −1 ) v m, u,v T M, x = π(g). (2.2) h ix h g ∗ g ∗ i ∈ x

The Ad(H)-invariance (2.1) is necessary and sufficient for the inner product (2.2) to be

5 invariant with respect to the action of G on M. The inner product defines an invariant metric on the homogeneous space M [33]. On the principal bundle (G,π,M,H), we introduce local coordinates (xi, hα), (i = 1, ... , dim M, α = 1,..., dim H) of the direct product U H, where xi are local coordinates × { } in a domain U M of a trivialization covering x U with local coordinates hα in ⊂ 0 ∈ { } the subgroup H, e = 0 . So the local coordinates of an element g = hs(x) from some H { } i i α α neighborhood of the identity Ge can be represented as g = x , g = h . We choose a section s : M G so that equalities sa(x)= xa and sα(x) = 0 hold over the domain U. → The tangent vectors

ea = s∗∂xa , eα = ∂hα , |x=x0 |h=eH form a basis e = e e , (A = 1,..., dim g) of the Lie algebra g T M T H, { A} { a}∪{ α} ≃ x0 ⊕ eH where e is a basis of the Lie algebra h T H, and e is a basis of the linear space { α} ≃ eH { a} m T M. In the conjugate space g∗ T ∗G, we introduce the dual basis eA by the ≃ x0 ≃ e { } condition eA, e = δA where δA is the Kronecker symbol. The right-invariant basis vector h Bi B B fields η (g) = (R ) e and the corresponding right-invariant 1-forms σA(g) = (R )∗eA A − g ∗ A − g in local coordinates (x, h) have the form

i α β ηa(x, h)= ηa(x, h)∂xi + ηa (x, h)∂hα , ηα(x, h)= ηα(h)∂hβ ,

a a i α α i α β σ (x, h)= σi (x, h)dx , σ (x, h)= σi (x, h)dx + σβ (h)dh , α,β =1,...,H.

∗ ′ ′ Here (Rg)∗, (Rg) are differentials of the right shifts Rg(g )= gg on the Lie group G. C The right-invariant vector fields ηA satisfy the commutation relations [ηA, ηB]= CABηC , while the right-invariant 1-forms σA satisfy the Maurer-Cartan relations, 2 dσC = CC σA − AB ∧ B A C σ . Here CAB = [eA, eB] are the structure constants of g with indices A,B,C = 1,..., dim G. We consider an invariant metric on the homogeneous space M using a coordinate system in a domain U of trivialization. A symmetric non-degenerate square matrix G = e , e m ab h a bi in a subspace m g satisfies the Ad(H)-condition ⊂ a a GabCcα + GacCbα =0, (2.3) a, b =1,..., dim M, α =1,..., dim H.

The invariant metric tensor in local coordinates xi is written as [38]: { } a b gij(x)= Gabσi (x, eH )σj (x, eH ), i, j =1,..., dim M. (2.4)

6 The contravariant components of the metric tensor are

ij ab i j ab −1 g (x)= G ηa(x, eH )ηb (x, eH ), G =(Gab) .

In what follows we will need the Christoffel symbols of the Levi-Civita connection with respect to a G-invariant metric gM given by [19, 33]

Γi (x) =Γa σb(x, e )σc(x, e )ηi (x, e ) jk bc j H k H a H − σb(x, e )ηi (x, e ) − j H b,k H − Ca σb(x, e )σα(x, e )ηi (x, e ). (2.5) − bα j H k H a H a ab Here i, j, k =1,..., dim M, and Γbc are determined by G of the quadratic form G and the structure constants of the Lie algebra g, 1 1 Γa = Ca Gad [G Ce + G Ce ] . (2.6) bc −2 bc − 2 ec bd eb cd Thus, in a homogeneous space with invariant metric, the Levi-Civita connection is defined by algebraic properties of the homogeneous space.

III. INDUCED REPRESENTATIONS AND INVARIANT FIRST-ORDER DIF- FERENTIAL OPERATOR WITH MATRIX COEFFICIENTS

Consider algebraic conditions for an invariant first-order linear differential operator with matrix coefficients on a homogeneous space M. We follow Ref. [35] where a more general case of invariant linear matrix differential operator of the second-order was studied. Denote by C∞(M,V ) and C∞(G,V ) the two spaces of functions that map a homogeneous space M and a transformation group G, respectively, to a linear space V . The last one can be regarded as a representation space of the algebra gl(V ). Functions on the homogeneous space M can be considered as defined on a Lie group G, but invariant over the fibers H of the bundle G [33]. In our case, when the functions take values in a vector space V , the space C∞(M,V ) is isomorphic to a subspace of the function space ˆ = ϕ C∞(G,V ) ϕ(hg)= U(h)ϕ(g), h H , F { ∈ | ∈ } where U(h) is an exact representation of the isotropy group H in V . For any function ϕ ˆ, ∈ F we have ϕ(g)= ϕ(hs(x)) = U(h)ϕ(s(x)), g =(x, h). (3.1)

7 Then we can identify ϕ(s(x)) with a function ϕ C∞(M, V ). Equation (3.1) gives an ∈ explicit form of the isomorphism ˆ C∞(M,V ). Differentiating relation (3.1) with respect F ≃ α to h and assuming h = eH , we obtain ∂U(h) (η +Λ ) ϕ(g)=0, Λ = , α =1,..., dimH. (3.2) α α α ∂hα h=eH

Here, Λα are representation operators of the algebra h on the space V . Equation (3.2) is a consequence of the condition (3.1) in the definition of ˆ. The isotropy subgroup H is F assumed to be connected. Then the conditions (3.1) and (3.2) are equivalent.

From (3.2) we can see that a linear differential operator R = R(g, ∂g) leaves invariant the function space ˆ, if F (η +Λ )R(g, ∂ )ϕ(g) = [η +Λ , R(g, ∂ )]ϕ(g)=0, ϕ ˆ. (3.3) α α g α α g ∈ F Thus, the space L( ˆ) of linear differential operators R(g, ∂ ): ˆ ˆ consists of linear F g F → F differential operators on C∞(G,V ) provided that

[η +Λ , R(g, ∂ )] =0. (3.4) α α g |Fˆ Then given relation (3.1), the action of R(g, ∂ ) L( ˆ) on a function ϕ(g) from the space g ∈ F ˆ is written as F −1 R(g, ∂g)ϕ(g)= U(h) U (h)RU(h) ϕ(s(x)). (3.5)

Multiplying equation (3.3) by U −1(h) and given η U(h)=
Λ U(h), we obtain α − α

−1 −1 U (h)[ηα +Λα, R(g, ∂g)]U(h)ϕ(s(x)) = [ηα, U (h)R(g, ∂g)U(h)]ϕ(s(x)) =

−1 ηα U (h)R(g, ∂g)U(h)ϕ(s(x)) =0.
−1 From here it follows that the operator U (h)R(g, ∂g)U(h) is independent of h and (3.5) can be written as

R(g, ∂g)ϕ(g)= U(h)RM (x, ∂x)ϕ(s(x)), (3.6)

−1 RM (x, ∂x) U (h)R(g, ∂g)U(h) = ≡ h=eH = R(g, ∂ )U(h) .
g ||h=eH That is, for any operator R(g, ∂ ) of L( ˆ) there exists an operator R on the homogeneous g F M ∞ space M acting on functions of the space C (M,V ). We say that the operator RM (x, ∂x)

8 is the projection of the operator R(g, ∂g): RM (x, ∂x) = π∗R(g, ∂g). For example, for a first-order linear differential operator b

a α R1(g, ∂g)= B (x, h)∂xa + B (x, h)∂hα + B(x, h)

the projection acts as follows:

(1) a α RM (x, ∂x)= π∗R1(g, ∂g)= B (x, eH )∂xa + B (x, eH )Λα + B(x, eH ). (3.7)

∞ On the other hand, any linearb differential operator RM defined on C (M,V ) corresponds to an operator R(g, ∂ )= U(h)R (x, ∂ )U −1(h) L( ˆ). g M x ∈ F Thus, we have the isomorphism L( ˆ) L(C∞(M,V )) whose explicit form is given by (3.6). F ≃ Let ξX (g)=(Lg)∗X be a left-invariant vector field on the Lie group G, where (Lg)∗ :

T ′ G T ′ G is the left shift differential L (g′)= gg′ on G, X g. g → gg g ∈ Since the left-invariant vector fields commute with right-invariant ones, the condition of projectivity (3.4) is fulfilled. Using (3.7), we find the corresponding operator on the homogeneous space as

X(x)=ˆπ ξ (g)= X(x)+ ξα (x, e )Λ , X g, (3.8) ∗ X X H α ∈ a a d tX a Xe(x)= ξ (x)∂xa , ξ (x)= (xe ) , (3.9) X X dt t=0

where X(x) are the generators of the action of the group G on M, note that X(x) act in the space C∞(M). It is easy to verify the following commutation relations for operators (3.8):

−1 −1 [X, Y ]= U (h)ξX U(h), U ξY U(h) = h=eH −1 −1 ^ = U (h)[ξX ,ξY ]U(h) = U (h)ξ[X,Y ]U(h ) = [X,Y ], e e h=eH h=eH     for all X,Y g. Consequently, the operators X corresponding to the left-invariant vector ∈ fields ξ are generators of a transformation group acting on C∞(M,V ) and are a contin- X e uation of the vector fields X(x) obeying the same commutation relations. We show how the projection of left-invariant vector fields is related to representations of the group on homogeneous space. The space ˆ invariant under right shifts on the Lie group G and the F ′ ′ operators Tg, acting by the rule (Tgϕ)(g ) = ϕ(gg ), define a representation of the group G

9 on ˆ induced by the representation U(h) of the subgroup H. According to relation (3.1), F we have:

(Tgϕ)(s(x)) = ϕ(gs(x)) = ϕ(h(x, g)s(xg)) = U(h(x, g))ϕ(s(xg)), (3.10) where h(x, g) H is the factor of the homogeneous space [37], which is determined from ∈ the system of equations

s(x)g = h(x, g)s(xg), h(x, e)= eH .

In view of the isomorphism ˆ C∞(M,V ), we obtain from (3.10) a representation of the F ≃ Lie group G on the space of functions ψ C∞(M,V ), ∈

(Tgψ)(x)= U(h(x, g))ψ(xg).

This representation is called the induced representation of the group G on the homogeneous space M. Note that

α α α α ξX (g)h (x, g)= ξX (g)(h(x, g)s(xg)) = ξX (g)(s(x)g) = ξX (s(x)g),

whence immediately follows the expression for the derivative of the factor at the identity element d hα(x, etX ) = ξα (x, e ). (3.11) dt X H t=0

It is easy to see that the operators X(x), as described by (3.9) and (3.11), are differentials of the representation T on the homogeneous space M: g e d X(x)ψ(x)= T ψ (x) . dt exp(tX) t=0
e Thus, the projection of left-invariant vector fields on the group gives the infinitesimal operators of the representation of Tg induced by the representation U(h) of the subgroup H. An operator R (x, ∂ ) L(C∞(M,V )) is invariant under the action of the Lie group of M x ∈ transformations, if RM (x, ∂x) commutes with X:

[R (x, ∂ ), X] = [U −1(h)RUe(h), U −1(h)ξ U(h)] = M x X |h=eH = U −1(h)[R(g, ∂ ),ξ ]U(h) =0. e g X |h=eH

10 It follows that the operator RM (x, ∂x) is invariant with respect to the transformation group if and only if the corresponding operator R(g, ∂ ) L( ˆ) commutes with the left-invariant g ∈ F vector fields: [R(g, ∂ ),ξ ]=0, X g. (3.12) g X ∈ Let R(1)(x, ∂ ) C∞(M,V ) be a linear differential operator of the first order, invariant with M x ∈ respect to the group action. By (3.12), this operator corresponds to a first-order polynomial of right-invariant vector fields:

a α R(1)(g, ∂g)= B ηa(x, h)+ B ηα(h)+ B.˜

α α As a result of the projection, the expression B ηα(h) becomes constant B Λα, which can be (1) ˜ α eliminated in the operator RM (x, ∂x) by changing the variable B = B + B Λα. Therefore, α we can put B = 0 without loss of generality. If we substitute the operator R(1)(g, ∂g) in the projectivity condition (3.4), then we obtain

a a [ηα +Λα, R(1)(g, ∂g)] Fˆ = ([b , Λα]ηa + b [ηa, ηα] + [B, Λα])Fˆ = a b a a β [B , Λ ]+ B C η ˆ + [B, Λ ] B C Λ =0. α bα a|F α − aα β Also we have a system of algebraic equations
for the coefficients Ba and B:

a b a [B , Λα]+ B Cbα =0, (3.13) [B, Λ ] BaCβ Λ =0. (3.14) α − aα β

When equations (3.13) – (3.14) are fulfilled, the projection of R(1)(g, ∂g) on the homogeneous space results in the desired form of the invariant linear differential operator of the first order:

(1) a i a α RM (x, ∂x)=ˆπ∗R(1)(g, ∂g)= B ηa(x, eH )∂xi + B ηa (x, eH )Λα + B. (3.15)

So, any linear differential operator of the first order acting on the functions of C∞(M,V ) and being invariant with respect to the action of the transformation group has the form (3.15) where the matrix coefficients Ba and B satisfy the algebraic system of equations (3.13) —

(3.14). The matrices Λα are generators of the isotropy subgroup H in a linear space V .

IV. λ-REPRESENTATION OF A LIE ALGEBRA

In this section we describe a special representation of the Lie algebra g using the orbit method [37]. The direct and inverse Fourier transforms on the Lie group G are introduced,

11 that in what follows are necessary for the noncommutative reduction of the Dirac equation on the homogeneous space M. Here we also use some results of the previous section. First, we describe an orbit classification for the coadjoint representation of Lie groups following conventions of Refs. [39, 40]. A degenerate Poisson–Lie bracket,

C ∂φ(f) ∂ψ(f) ∞ ∗ φ, ψ (f)= f, [dφ(f), dψ(f)] = CABfC , φ,ψ C (g ), (4.1) { } h i ∂fA ∂fB ∈

∗ endows the space g with a Poisson structure. Here fA are coordinates of a linear functional f = f eA g∗ relative to the dual basis eA . The number ind g of functionally independent A ∈ Casimir functions Kµ(f) relative to the bracket (4.1) is called the index of the Lie algebra g. A coadjoint representation on g∗, Ad∗: G g∗ g∗, stratifies g∗ into orbits of the × → coadjoint representation (K-orbits). The restriction of the bracket (4.1) on orbits is non- degenerate and coincides with the Poisson bracket generated by the symplectic Kirillov form

ωλ. The orbits of maximal dimension dim (0) = dim g ind g are called non-degenerate, O − and the those of less dimension are singular. We denote by (s) the orbits of dimension Oλ dim g ind g 2s, s = 0,..., (dim g ind g)/2 passing through the functional λ g∗, and − − − ∈ (s) a number s is called the orbit singularity index. A tangent space T (s) to the orbit f Oλ Oλ at a point f is the linear span of vector fields

∂ C YA(f)= CAB(f) , CAB(f)= CABfC , ∂fB

so that the orbit dimension is given by the rank of the matrix CAB(f). The rank takes a constant value on the orbit, dim (s) = rank C (λ). The Kirillov form on tangent vectors Oλ AB to the orbit (s) is ω (Y ,Y )= λ, [Y ,Y ] . Oλ λ A B h A B i ∗ The space g can be decomposed into a sum of disjoint invariant algebraic surfaces Ms consisting of orbits of the same dimension dim g ind g 2s: − −

M = f g∗ F 1(f)=0 , 0 ∈ | ¬ dim g ind g M = f g∗ F s(f)=0,
F s+1(f)=0 , s =1,..., − 1, s ∈ | ¬ 2 − − ∗ dimg indg M  = f g F 2 (f)=0
, (dim g−ind g)/2 ∈ | n o 12 where F s(f) denotes the set of all minors of the matrix C (f) = CC f of size dim g AB AB C − ind g 2s +2; the notation F s(f) = 0 implies that all the corresponding minors at the point − f vanish, and (F s(f) = 0) means that at the point f, the corresponding minors do not ¬ vanish simultaneously. In the general case, the surface Ms is disconnected.

In what follows, by M(s) we denote a connected component of the surface Ms containing the orbit (s), and (s) will be called the (s)-type orbit. Each component of M is uniquely O O (s) (s) determined by a set of homogeneous polynomials Fα (f) satisfying the system

(s) YA(f)Fα (f) F (s)(f)=0 =0. (s) The non-constant functions Kµ (f) on M( s) are called the (s)-type Casimir functions if

they commute with any function on M(s) with respect to the bracket (4.1). In other words, (s) the functions Kµ (f) are invariants of the adjoint representation and are determined by the system (s) YA(f)Kµ (f) =0, µ =1,...,r(s). f∈M(s)

The number of functionally independent solutions of this system is determined by the di-

mension of the surface M(s):

r = dim M (dim g ind g 2s). (s) (s) − − − Denote by Ω(s) Rr(s) a set of values of the mapping K(s) : M Rr(s) and introduce a ⊂ (s) → locally invariant subset

(s) = f M K(s)(f)= ω(s), µ =1,...,r , ω(s) Ω(s) . Oω ∈ (s) | µ µ (s) ∈   (s) (s) If the Casimir functions Kµ (f) are single-valued, then the level surface ω consists of a O (s) ∗ countable set of orbits. We call ω the class of orbits. As a result, the space g consists of O a union of connected invariant disjoint algebraic surfaces M(s), which in turn is the union of (s) the orbit classes ω : O g∗ = M = (s). (4.2) (s) Oω (s) (s) [(s) [(s) ω [∈Ω Consider a quotient space B(s) = M(s)/G, dimB(s) = r(s), whose points are the orbits of one class, (s) M . We introduce a local section λ(j) of the bundle M with base B Oλ ∈ (s) (s) (s) using real parameters j =(j ,...,j ) taking their values in a domain J Rr(s) : 1 r(s) ⊂ (s) (s) (s) (s) ∂ωµ (j) F (λ(j))=0, Kµ (λ(j)) = ωµ (j), det =0. ∂jν 6

13 Let (s) be a K-orbit of (s)-type passing through a covector λ = λ(j) g∗ and belonging Oλ(j) ∈ to the same class of orbits for all j J. ∈ Using the Kirillov orbit method [37], we construct a unitary irreducible representation of the Lie group G on a given orbit. This representation can be constructed if and only if for the functional λ there exists a subalgebra p gC in the complex extension gC of the Lie ⊂ algebra g satisfying the conditions: 1 λ, [p, p] =0, dim p = dim g dim (s). (4.3) h i − 2 Oλ The subalgebra p is called the polarization of the functional λ. In (4.3), it is assumed that the functionals from the space g∗ are extended to gC by linearity. Moreover, real polarizations always exist for nilpotent and completely solvable Lie algebras, and the complex polarizations always exist for solvable Lie groups [41]. For non-degenerate orbits (0) there always exists, Oλ generally speaking, a complex polarization. In this paper, for simplicity, we restrict ourselves to the case when p is the real polarization. Denote by P a closed subgroup of the Lie group G whose Lie algebra is p. The Lie group acts on the right homogeneous space Q G/P : q′ = qg. According to the orbit method, we ≃ introduce a unitary one-dimensional irreducible representation of the Lie group P , which, in the neighborhood of the identity element of the group, has the form i U λ(eX ) = exp λ,X , X p. (4.4) ~h i ∈   The representation of the Lie group G corresponding to the orbit (s) is induced using Oλ (4.4) as

λ ∆G(p(q,g)) λ λ+i~β (Tg ψ)(q)= U (p(q,g))ψ(qg)= U (p(q,g))ψ(qg), (4.5) s∆P (p(q,g)) ∆ (p) β = d log G , s∆H (p) u=eP

−1 where ∆G(g) = det Adg is the module of the Lie group G, ∆P (p) = detAdp is the module of the subgroup P , p P , and e is the identity element of P . A function p(q,g) is the ∈ P factor of the homogeneous space Q. λ λ The functions ψ (q; g)=(Tg ψ)(q) on the group G satisfy a condition similar to (3.1):

ψλ(q; pg)= U λ+i~β(p)ψλ(q; g), u P. ∈ 14 The space of all such functions will be denoted by λ. Restriction of the left-invariant vector F fields ξX (g) to a homogeneous space Q, as follows from results of section III, is correctly defined, and the explicit form of the corresponding operator on the homogeneous space is given by (3.6):

λ+i~β −1 λ+i~β ℓX (q, ∂q,λ)= U (p) ξX (g)U (p) , (4.6) p=eP   [ℓ (q, ∂ ,λ),ℓ (q, ∂ ,λ)] =ℓ (q, ∂ ,λ),X,Y g. (4.7) X q Y q [X,Y ] q ∈

Equation (4.6) shows that ℓX (q, ∂q,λ) are infinitesimal operators of the induced repre- sentation (4.5), d ℓ (q, ∂ ,λ)ψ(q)= (T λ ψ)(q) . X q dt exp(tX) t=0

Denote by L(Q, h,λ) a space of functions on Q where representation (4.5) is defined. The representation (4.5) is unitary with respect to a scalar product of the function space

L2(Q, h,λ):

1 dim Q (ψ1, ψ2)= ψ1(q)ψ2(q)dµ(q), dµ(q)= ρ(q)dq ...dq . (4.8) ZQ The function ρ(q) is determined from the Hermitian condition for the operators iℓ (q, ∂ ,λ) − X q with respect to this scalar product (4.8). The irreducible representation of the Lie algebra g by the linear operators of the first order (4.6) dependent on dim (s)/2 variables is called λ - representation of the Lie algebra Oλ g and it was introduced in Ref. [16]. Let us describe a recipe for calculating operators of the λ-representation on the class of (s) orbits ω . We find a local section λ(j) of the bundle M with the base B for which there O (s) (s) exists a polarization p with a basis e . Let us then fix the basis e′ in the additional { α} { a} subspace m = p⊥ to the subalgebra p. In a neighborhood of the identity element of the Lie group G, we introduce local coordinates of the second kind,

− dim Q ′ dim Q−1 ′ ′ pdim pe pdim p 1e p1e q e q e q1e g(p, q)= e dim p e dim p−1 ...e 1 e dim Q e dim Q−1 ...e 1 ,     and write the left-invariant vector fields in local coordinates (p, q):

a α ξX (p, q)= ξX (q)∂qa + ξX (q,p)∂pα .

15 Next, operators of the λ representation are defined using (4.6) as

a i α ~ ℓX (q, ∂q,λ)= ξX (q)∂qa + ~ξX (q, eP )(λα + i βα) , 1 β = Tr (ad ) Tr ad . α 2 α − α|p h  i In other words, finding of these operators is reduced to calculating the left-invariant vector fields on the group G in the trivialization domain of the principal bundle of this group in the fibrations P = exp(p). Let us write the representation operators (4.5) in integral form

λ ′ λ (Tg ψ)(q)= ψ(q )Dqq′ (g)dµ(q), (4.9) ZQ

λ ∆G(p(q,g)) λ ′ Dqq′ (g)= U (p(q,g))δ(qg,q ), s∆H (p(q,g))

where δ(q, q′) is the generalized delta-function with respect to the measure dµ(q). The λ generalized kernels Dqq′ (g) of this representation satisfy the following properties:

λ λ λ ′′ Dqq′ (g1g2)= Dqq′′ (g1)Dq′′q′ (g2)dµ(q ), ZQ λ λ −1 λ ′ Dqq′ (g)= Dq′q(g ),Dqq′ (e)= δ(q, q ),

and the system of equations

λ (ηX (g)+ ℓX (q, ∂q,λ)) Dqq′ (g)=0,

′ ′ λ ξX (g)+ ℓX (q , ∂q ,λ) Dqq′ (g)=0. (4.10)   Here g ,g G. Let Gλ = g G Ad∗λ = λ be a stabilizer of the functional λ g∗ with 1 2 ∈ ∈ | g ∈ λ λ λ the Lie algebra g , and the vectors eα,α =1,..., dim g form some fixed basis for g p. { } ⊂ As P is the stabilizer of the point q = 0 in the homogeneous space Q, and the group Gλ lies in P , we get the equality

i ~ λ ℓα(0, ∂q,λ)= ~ (λα + i βα) , α =1,..., dim g .

Restricting the first equality (4.10) to the subgroup Gλ and setting q = 0, we find

i λ λ η (h)+ (λ + i~β ) D ′ (h)=0, h G . (4.11) α ~ α α 0q ∈   16 The solution of the system (4.11) up to a constant factor can be represented as

λ i λ α D ′ (h) = exp σ (h)+ σ (h)β . (4.12) 0q −~ α  Z Z  λ λ α The subalgebra g is subordinate to the covector λ, and the 1-forms σ (h) and σ (h)βα are closed in Gλ. Thus, the integral in (4.12) is well-defined. The local solution (4.12) can be extended to a global one, if the integral on the right-hand side of (4.12) over any closed curve Γ on the subgroup Gλ is a multiple of 2πi. Note that since the 1-form σλ(h) is closed, the value of this integral depends only on the homological class to which the curve Γ belongs. Therefore, for a global solution of the system (4.11), the following condition should be satisfied: 1 σλ(h)= n Z. (4.13) 2π~ Γ ∈ λ Γ∈HI1(G ) In other words, the 1-form σλ(h) should belong to an integral cohomology class from H1(Gλ, Z). In the case of a simply connected group G, the condition (4.13) is equivalent to the condition of integral orbit (s), proposed by A. A. Kirillov [37]: Oλ 1 ω = n Z. 2π~ λ Γ ∈ 2 Γ∈HZ(Oλ)

Thus, for a simply connected group the coadjoint orbit (s) is integral if the equality (4.13) Oλ is fulfilled. λ A set of generalized functions Dqq′ (g) satisfying the system (4.10) was studied in Refs. [16, 39] and the hypothesis was proposed that this set of generalized functions has the properties of completeness and orthogonality for a certain choice of the measure dµ(λ) in the parameter space J:

λ λ ′ ′ D ′ (g)D ′ (g)dµ(g) = ∆ (s(q))δ(q, q˜)δ(˜q , q )δ(λ,λ˜ ), (4.14) qeq qq G G Z ee ′ λ λ dµ(q ) D ′ (˜g)D ′ (g)dµ(q) dµ(λ)= δ(˜g,g). (4.15) qq qq ∆ (s(q)) ZQ×Q×J G Here δ(g) is the generalized Dirac delta function with respect to the left Haar measure dµ(g) on the Lie group G. For compact Lie groups, the relations (4.14)-(4.15) hold by the Peters-Weil theorem [42]. Note that although there is no rigorous proof of the relations (4.14)-(4.15), in each case it is easy to verify directly their validity.

17 Consider the space L(G,λ,dµ(g)) of functions of the form

λ ′ λ −1 ′ ψ (g)= ψ(q, q ,λ)Dqq′ g dµ(q )dµ(q). (4.16) ZQ
Here, a function ψ(q, q′,λ) of the two variables q and q′ belongs to the space L(Q, h,λ). The inverse transform reads

′ −1 ′ λ −1 ψ(q, q ,λ) = ∆G (s(q )) ψλ(g)Dqq′ (g )dµR(g), (4.17) ZQ×Q×J −1 where we have used (4.14)–(4.15), and dµR(g)= dµ(g ) is the right Haar measure on the Lie group G. λ The action of the operators ξX (g) and ηX (g) on the function ψ (g) from L2(G,λ,dµ(g)), † according to (4.16) and (4.17), corresponds to action of the operators ℓX (q, ∂q,λ) and ′ ′ ℓX (q , ∂q′ ,λ) on the function ψ(q, q ,λ) respectively,

ξ (g)ψλ(g) ℓ† (q, ∂ ,λ)ψ(q, q′,λ), X ⇐⇒ X q λ ′ ′ η (g)ψ (g) ℓ (q , ∂ ′ ,λ)ψ(q, q ,λ). X ⇐⇒ X q (s) (s) The functions (4.16) are eigenfunctions for the Casimir operators Kµ (i~ξ)= Kµ ( i~η). − Indeed, from the system (4.10) we can obtain

K(s)(i~ξ)ψλ(g)= K(s)( i~η)ψλ(g)= µ µ − (s) ′ λ −1 ′ Kµ ( i~ℓ(q, ∂ ,λ))ψ(q, q ,λ) D ′ g dµ(q )= − q qq ZQ  
(s) ′ ′ λ −1 ′ K ( i~ℓ(q , ∂ ′ ,λ))ψ(q, q ,λ) D ′ g dµ(q ). µ − q qq ZQ

(s) ′ ′ It follows that the operators Kµ ( i~ℓ(q ,λ)) are independent of q and − K(s)(i~ξ)ψλ(g) κ(s)(λ)ψ(q, q′,λ), µ ⇐⇒ µ (s) ′ (s) (s) (s) K ( i~ℓ(q , ∂ ′ ,λ)) = κ (λ), κµ (λ)= κ (λ), µ − q µ µ (s) (s) lim κµ (λ)= ωµ (λ). ~→0 Thus, as a result of the generalized Fourier transform (4.16), the left and the right fields become the operators of λ-representations, and the Casimir operators become constants. This fact is a key point for the method of noncommutative integration of linear partial differential equations on Lie groups, since it allows one to reduce the original differential equation with dim G independent variables to an equation with fewer independent variables equal to dim Q.

18 V. DIRAC EQUATION IN HOMOGENEOUS SPACE

In this section, we consider the Dirac equation in a n-dimensional homogeneous space M with an invariant metric. We shall assume that in the homogeneous space M an invariant

metric gM and the Levi-Civita connection are given. Denote by VΨ a space of spinor fields on M. We write the Dirac equation in the space M as an equation in a n-dimensional Lorentz manifold M with the metric (~ is the Planck constant) as follows [43]:

i~γi(x)[ +Γ (x)] m ψ(x)=0, i =1, . . . , n. (5.1) ∇i i − Here is the covariant derivative corresponding
to the Levi-Civita connection on M, m ∇i is mass of the field ψ C∞(M,V ), ψ(x) is a column with 2⌊n/2⌋ components, γi(x) are ∈ Ψ 2⌊n/2⌋ 2⌊n/2⌋ gamma matrices, × γ (x),γ (x) =2g (x)E, i, j =1, . . . , n, (5.2) { i j } ij where E denotes the 2⌊n/2⌋ 2⌊n/2⌋ identity matrix, Γ (x) is the spinor connection satisfying × i the conditions [ ,γ (x)]=0, Tr Γ (x) = 0. The spinor connection Γ (x) can be written as ∇i j i i follows [43]: 1 Γ (x)= ( γ (x))γk(x). (5.3) i −4 ∇i k We seek a solution of (5.2) with the decomposition

i a i a i a γ (x)=ˆγ ηa(x, eH ), γˆ = γ (x)σi (x, eH ). (5.4)

The constant matricesγ ˆa satisfy the algebraic equations

γˆa, γˆb =2GabE, a, b =1, . . . , n. (5.5) { } For the Dirac matrices with subscripts using (2.4) we have

γ (x)= g (x)γi(x)=ˆγ σa(x, e ), γˆ G γˆa. (5.6) i ij a i H b ≡ ab The spinor connection is given by the following Lemma.

i Lemma 1 The spinor connection Γ(x) = γ (x)Γi(x) on the homogeneous space M with

invariant metric gM reads 1 Γ(x)=ˆγa (Γ + ηα(x, e )Λ ) , Γ = Γd γˆbγˆ , a a H α a −4 ba d 1 Λ = G Ca [ˆγb, γˆc]. (5.7) α −8 ac αb

19 i Proof The function Γ(x)= γ (x)Γi(x) with Γi(x) given by (5.3), can be written as

1 i k l Γ(x)= γ (x)γ (x) ∂ i γ (x) Γ (x)γ (x) , (5.8) 4 x k − ki l
l where Γki(x) are the Christoffel symbols, and ∂xi is the partial derivative. Substituting (2.5), (5.4), and (5.6) in (5.8), we obtain

1 Γ(x)=Γ+ Cd γˆaγˆbγˆ σα(x, e )ηi (x, e ). 4 bα d i H a H

d b Using property (2.3) of the invariant metric, we reduce the expression Cbαγˆ γˆd to the form 1 Cd γˆbγˆ = Cd G γˆbγˆc = Cd G γˆbγˆc = Cd G [ˆγb, γˆc]=4Λ . bα d bα dc − cα db 2 bα dc α

From the chain of equalities

σα(x, e )ηi (x, e )= σα(x, e )ηA(x, e ) σα(e )ηβ(x, e )= i H a H A H a H − β H a H = δα ( δα)ηβ(x, e )= ηα(x, e ), a − − β a H a H

we obtain for the spinor connection the required expression (5.7).

Thus, the Dirac equation in the homogeneous space M with an invariant metric gM and the Dirac matrices of the form (5.4) takes the form

(x, ∂ )ψ = mψ, DM x a i α (x, ∂ )= i~γˆ η (x, e )∂ i +Γ + η (x, e )Λ . DM x a H x a a H α   A set of matrices Λα determines a spinor representation of the isotopy subgroup H in the space VΨ.

Lemma 2 The matrices Λα are generators of the isotropy subgroup H representation on the space VΨ.

Proof We prove that the matrices Λα satisfy the commutation relations

γ [Λα, Λβ]= CαβΛγ. (5.9)

The commutator of Λα and Λβ can be written as

1 1 [Λ , Λ ]= Cd [Λ , γˆbγˆ ]= Cd [Λ , γˆb]ˆγ +ˆγb[Λ , γˆ ] . (5.10) α β −4 βb α d −4 βb α d α d
20 a Using (2.3), (5.5) and (5.6), we find the commutator of Λα withγ ˆ : 1 1 [Λ , γˆa]= Cd [ˆγbγˆ , γˆa]= Cd δaγˆb Gabγˆ = α 4 bα d 2 bα d − d 1 = GbdCa GabCd γˆ = Ca γˆb.
(5.11) 2 bα − bα d bα
Similarly, for the γ-matrices with lower indices we have

b [Λα, γˆa]= Cαaγˆb. (5.12)

Substitution of (5.11)–(5.12) in (5.10) yields 1 [Λ , Λ ]= Cd Ce Ce Cd γˆbγˆ . (5.13) α β 4 βe βb − βb αe d
The expression inside the parentheses can be written in the form

Cd Ce Ce Cd = CA Cd + CA Cd + CA Cd + Cγ Cd . (5.14) βe βb − βb αe αb βA bβ αA βα bA αβ bγ   By the Jacobi identity for the structure constants, the expression inside the square brackets is equal to zero. Substituting (5.14) in (5.13), we obtain (5.9). The Dirac operator (x) is a differential operator of the first order with matrix coef- DM ∞ ficients acting on the functions from C (M,VΨ) (spinors). Let us associate this operator with the corresponding operator on the transformation group. The key theorem of our work is

Theorem 1 The Dirac operator (x, ∂ ) on the homogeneous space M with invariant DM x metric gM is a projection of an operator

(x, ∂ )=ˆπ (g, ∂ ), DM x ∗DG g (g, ∂ ) i~γˆa[η (g)+Γ ] L( ˆ ). (5.15) DG g ≡ a a ∈ FΨ Proof Comparing the Dirac operator (x, ∂ ) with invariant matrix differential operator DM x of the first order (3.15) on the homogeneous space M, we obtain

a a a B = i~γˆ , B = i~γˆ Γa. (5.16)

Projection (5.15) is determined if the coefficients Ba and B of the form (5.16) satisfy equa- a tions (3.13)-–(3.14). From (5.11) it follows that the commutator of Λα andγ ˆ satisfies the first condition in (3.13). In this case the condition (3.14) is reduced to

β a [Γ, Λα]= Caαγˆ Λβ. (5.17)

21 The commutator of Γ and Λα can be presented in terms of the commutator [Λα, Γa] as

a a [Λα, Γ] = [Λα, γˆ ]Γa +ˆγ [Λα, Γa]. (5.18)

Using (5.7) and (5.11)–(5.12), we get 1 [Λ , Γ ]= Γd [Λ , γˆb]ˆγ +ˆγb[Λ , γˆ ] = α a −4 ba α d α d 1 = Γd Cc γˆbγˆ Cb γˆcγˆ .
(5.19) −4 ba αd c − αc d
Substituting (5.19) into (5.18), we find 1 [Λ , Γ] = (Ce Γc Cc Γe + Ce Γc )ˆγaγˆbγˆ . (5.20) α 4 αa be − αc ba αb ea c From (2.6) and the Jacobi identity for structure constants of the Lie algebra g it follows that

c e e c e c β c CαeΓba = CαaΓbe + CαbΓea + CαaCβb. (5.21)

Substituting (5.21) into (5.20), we obtain (5.17). Thus, relations (3.13) and (3.14) are satisfied and the lemma 1 holds. It follows that the Dirac operator (x, ∂ ) can be obtained DM x a a by projection of the operator B ηa + B where B and B are given by (5.16), and we come to (5.15). From this statement we immediately obtain

Corollary 1 The generators

a α X(x)= ξX (x)∂xa + ξX (x, eH )Λα (5.22)

of representation of the Lie algebrae g in the space VΨ are symmetry operators of the Dirac operator on the homogeneous space M, i.e. DM [X(x), (x, ∂ )]=0, DM x ^ [X(x), X(x)] = [X,Y ](x),X,Y g. e ∈ In view of the isomorphisme ˆ eC∞(M,V ) and Theorem 1, the Dirac equation (5.1) on F ≃ Ψ M is equivalent to the following system of equations on the transformation group G:

(g, ∂ )ψ(g)= mψ(g), DG g  (5.23) (ηα +Λα) ψ(g)=0.

The function ϕ(x)= ψ(s(x)) = ψ(x, eH ) provides a solution of the original Dirac equation (5.1) on the homogeneous space M.

22 VI. NONCOMMUTATIVE INTEGRATION

We will look for a solution to the system (5.23) as a set of functions

′ λ −1 ′ ψσ(g)= ψσ(q )Dqq′ g dµ(q ), σ =(q,λ), (6.1) ZQ
′ where the function ψσ(q ) is a spinor, each component of which belongs to the function space ′ λ −1 L(Q, h,λ) with respect to the variable q , and Dqq′ (g ) is introduced by (4.9). Using (4.10) we can then reduce the system (5.23) to the equations

(q′, ∂ ′ ,λ)ψ (q′)= mψ (q′), Dℓ q σ σ  (6.2) ′ ′ (ℓα(q ,λ)+Λα) ψσ(q )=0,

where 

′ a ′ (q , ∂ ′ ,λ)= i~γˆ [ℓ (q ,λ)+Γ ]. (6.3) Dℓ q a a

We call the operator (q′, ∂ ′ ,λ) in (6.3) the Dirac operator in the λ-representation. The Dℓ q ′ number of independent variables q in (6.2) is dim Q. The set of functions ψσ(x)= ψσ(s(x)) gives us a solution of the original Dirac equation (5.1) on the homogeneous space M. It follows from the equations

ξ (g)ψ (g)= U(h)X(x)ψ (x)= ℓ (q,λ)ψ (g)= U(h) ℓ (q,λ)ψ (x) , X σ σ − X σ {− X σ } e that the solutions ψσ(x) of the Dirac equation satisfy the system

X(x)+ ℓX (q,λ) ψσ(x)=0, (6.4) n o where X is given by (3.8). The algebraice relations between operators of the λ-representation should correspond to the algebraic relations between the generators X(x) for compatibility e of the system (6.4). More precisely, the corollary of the system (6.4) is to be fulfilled for any e homogeneous function F of X(x):

F (Xe (x))ψ (x)= F ( ℓ (q, ∂ ,λ))ψ (x). σ − X q σ

This condition is obviously satisfiede for the commutator of two operators (X,Y g), ∈

[X(x), Y (x)]ψ (x)= [ℓ (q, ∂ ,λ),ℓ (q, ∂ ,λ)]ψ (x), σ − X q Y q σ e e 23 and for the Casimir functions, we have

(s) ~ −1 (s) ~ Kµ (i X)ψσ(x)= U (h)Kµ (i ξ(g))U(h) ψσ(x)= −1 (s) ~ (s)
= U (eh)Kµ (i ξ(g))ψσ(g)= κµ (λ)ψσ(x).

The homogeneous functions Γ C∞(g∗) provided that Γ(X(x)) 0, can exist on the dual ∈ ≡ space g∗ to the space M. These functions are called identities on the homogeneous space M.

The number of functionally independent identities iM is called the index of the homogeneous space M. In Ref. [40] it was shown that any homogeneous function Γ C∞(g∗) satisfying ∈ the condition

Γ(f) ⊥ =0, (6.5) |f∈h

(sM ) is an identity. In the same Ref. [40] it was shown that the functions Fα (f), where sM =

λ (sM ) dim g ind g /2, are identities, and all other identities are Casimir functions Kµ (f) − such that

(sM ) K (f) ⊥ 0. µ f∈h ≡

A description of identities for the generators X(x) of the form (3.8) is given by the following Lemma which is essential result of our work: e

(sM ) Lemma 3 The identities for the operators X(x) are generated by the functions Fα (f)

(sM ) and Kµ (f). e Proof Suppose that a homogeneous function Γ′ C∞(g∗) is an identity for the generators ∈ X(x), i.e., Γ′(X) 0. Then the symbol ≡ e e Γ′(x, p)=Γ′(X(x, p)) C∞(g∗,V ) ∈ of the operator Γ′(X) also equals zero fore all (x, p), and

e a α X(x, p)= X (x)pa + ξX (x, eH )Λα, where the constants p are coordinatese of the covector f x = p dxa T ∗M. At a given point a a ∈ x 0 (x0,p ), we have

0 0 0 Xa(x0,p )= pa, Xα(x0,p )=Λα, f x0 = p0dx0 T ∗ M h⊥. e a ∈ ex0 ≃

24 Expanding Γ′(X(x, p)) in terms of the basis B of matrices in the vector space V and putting x = x , p = f x0 , we get: 0 e

Γ′(X(x ,p0))=Γ′(p0,...p0 , Λ ,..., Λ )= BσΓ (p0,...p0 ) 0. 0 1 dim M 1 dim H σ 1 dim M ≡ e As a result, for each function Γσ(f) we come to equation (6.5). The last one shows that ′ the functions Γσ(f) are identities on a homogeneous space, and the identities Γ (f) for the operators X(x) have the following structure:

e ′ ′ α Γ (X(x))=Γ (X(x)+ ξ (x, eH )Λα)= = BσΓ (X(x)) 0. e σ ≡

From this one can see that the number of functionally independent identities between X(x) does not exceed the index i of the homogeneous space, and the functions Γ (f) depend on M σ e identities on the homogeneous space. For the compatibility of the system (6.4), we have to take into account the identities between the generators X(x), Γ′(X) 0; namely we impose the following conditions on the ≡ operators of the λ-representation: e e

Γ′( ℓ(q, ∂ ,λ)) 0. (6.6) − q ≡

A class of orbits and corresponding parameters j should be restricted by (6.6). For instance, for the case X(x)= X(x), the condition (6.6) is reduced to

e F (sM )( i~ℓ(q, ∂ ,λ)) 0, α − q ≡ K(sM )( i~ℓ(q, ∂ ,λ)) = κ(sM )(λ(j))=0. (6.7) µ − q µ

The first condition in (6.7) says that the λ-representation has to be constructed by the class of orbits (sM ), and the second one imposes a restriction on the parameters j. In Oλ(j) Ref. [44], a λ-representation satisfying (6.7) is called a λ-representation corresponding to the homogeneous space M. Thus, condition (6.7) is stronger than (6.6). One of the important results of our work is the fact that when performing a noncommutative reduction of the Dirac equation, it is necessary to use the weaker condition (6.6) for the correct application of the noncommutative integration method.

25 The second equation of the system (6.2) can be written as

∂ i αa (q′) + ξα(q′, e )(λ + i~β )+Λ ψ (q′)=0. (6.8) α ∂q′a ~ α P α α α σ   We look for a solution of (6.8) in the form

′ ′ ψσ(q )= Rσ(q )ψσ(v),

′ where Rσ(q ) is a certain function, and ψσ(v) is an arbitrary function of the characteristics v = v(q′) of the system (6.2). We carry out a one-to-one change of variables q′ = q′(v,w), where w = w(q′) are some coordinates additional to v. By V and W we denote domains of the variables v and w, respectively. The measure dµ(q′) in the new variables takes the form dµ(q′) = ρ(v,w)dµ(v)dµ(w). Then the solution of the original Dirac equation can be represented as

λ ψσ(x)= ψσ(v)Dqv (x) dµ(v), ZV λ λ −1 Dqv (x)= Rσ(v,w)Dqq′(v,w) g ρ(v,w)dµ(w). ZW
Substituting the solution ψσ(x) into the Dirac equation (6.2), we obtain a linear first-order differential equation with matrix coefficients for the function ψσ(v) with the number of independent variables v equal to dim V .

VII. THE METRIC THAT DOES NOT ADMIT SEPARATION OF VARIABLES IN THE DIRAC EQUATION

Consider a four-dimensional homogeneous space M with a transformation group G whose Lie algebra g is defined in some basis e by nonzero commutation relations { A}

[e , e ]= e , [e , e ]= e , [e , e ]= e , 1 4 − 1 1 5 2 2 3 1 [e , e ]= e , [e , e ]= 2e , [e , e ]= e , 2 4 2 3 4 − 3 3 5 4 [e , e ]= 2e . 4 5 − 5

The Lie algebra g is a semidirect product of the two-dimensional commutative ideal R2 = span e , e and the three-dimensional simple algebra sl(2) = span e , e , e . We also take { 1 2} { 3 4 5} h = e as the one-dimensional subalgebra. { 5} 26 Denote by (xa, hα) local coordinates on a trivialization domain U of the group G so that

1 4 3 2 1 g(x, h)= eh e5 ex e4 ex e3 ex e2 ex e1 , (7.1) x =(x1, x2, x3, x4) M. ∈

The group G is unimodular and ∆G(g) = 1. A symmetric non-degenerate matrix

0 00 c − 3  0 c4 c3 c2  (Gab)= , ck = const, c3 =0,  0 c 0 0  6  3     c c 0 c  − 3 2 1    defines an invariant metric on the space M,

x4 e 4 4 4 2 ds2 = 2c e2x dx1 x3dx2 dx3 c ex dx1 x3dx2 2 c e3(x ) d(x3)2+ c2 2 − − 1 − − 4 3 

+2c dx3 + x3dx4 dx2 2c dx1dx4 . (7.2) 3 − 3 
The metric (7.2) has nonzero scalar curvature R =6c1. The group generators in the canon- ical coordinates (7.1) have the form:

2 X1(x)= ∂x1 , X2(x)= ∂x2 , X3(x)= x ∂x1 + ∂x3 ,

1 2 3 X (x)= x ∂ 1 + x ∂ 2 2x ∂ 3 + ∂ 4 , 4 − x x − x x 1 3 2 3 X (x)= x ∂ 2 (x ) ∂ 3 + x ∂ 4 . 5 x − x x

The vector fields ξA(x, eH ), in turn, are determined by the expressions

−2x4 ξa(x, h)= Xa(x), ξ5(x, h)= X5(x)+ e ∂h1 .

The right-invariant vector fields in the canonical coordinates (7.1) are

−x4 3 x4 x4 −x4 3 η (x, h)= e + h x e ∂ 1 h e ∂ 2 , η (x, h)= e (x ∂ 1 + ∂ 2 ), 1 − 1 x − 1 x 2 x x  −2x4  2 η (x, h)= e ∂ 3 h ∂ 4 + h ∂ 1 , η (x, h)= ∂ 4 +2h ∂ 1 , 3 − x − 1 x 1 h 4 − x 1 h η (x, e )= ∂ 1 . 5 H − h

The gamma matricesγ ˆa can be presented in terms of the standard Dirac gamma matrices

27 γa as follows:

1 1 2 2 c2 1 2 3 γˆ = γ + γ , γˆ = γ + γ + √ c4γ , −c3 − c γˆ3 = 3 γ3 + iγ4 ,
−√ c − 4 c
2c c γˆ4 = 3 γ1 γ2 2 γ3 + iγ4 1 γ1 + γ2 . − 2 − − √ c − c − 4 3


The spin connection is independent of local coordinates and has the form

1 Γ(x)=Γ γˆ4Λ, Λ= γˆ1γˆ3, − 2c3 c1 1 1 3 c2 1 3 4 1 2 3 4 c1 1 3c3 3 4 Γ= 2 γˆ γˆ γˆ + 2 γˆ γˆ γˆ + γˆ γˆ γˆ γˆ + γˆ 2ˆγ . 4c3 4c3 c3 − 4c3 4c3 −

The Dirac operator in local coordinates is

1 3 2x4 2 x4 2 −2x4 3 4 4 = i~ γˆ + x e γˆ ∂ 1 + e γˆ ∂ 2 e γˆ ∂ 3 +ˆγ ∂ 4 +Γ γˆ Λ. DM − y x − x x −  h i  The first-order symmetry operators are defined by (5.22):

−2x4 Xa(x)= Xa(x), X5(x)= X5(x)+ e Λ.

e e The metric (7.2) generally does not admit the Yano vector field and the Yano-Killing tensor field, so the Dirac equation does not admit spin symmetry operators. As a result, the Dirac equation has only two commuting symmetry operators X (x), X (x) of the first { 1 2 } order. However, the Dirac equation admits a third-order symmetry operator e e

K(X)=X (x) X (x) X (x)+ X (x) X (x) X (x) X (x) X (x) X (x), (7.3) 5 · 1 · 1 1 · 2 · 4 − 2 · 2 · 4

where Xe Y e= (XYe + YXe)/2 is thee symmetrizede e producte of thee operatorse X and Y . · As a consequence, the metric (7.2) does not admit separation of variables for the Dirac equation. Note that the Klein-Gordon equation also admits only three commuting symmetry operators X (x),X (x),K(X) . One of them, K(X), is the third-order operator, and the { 1 2 } Klein-Gordon equation also does not admit separation of variables. We now carry out a noncommutative reduction of the Dirac equation. First, we describe orbits of the coadjoint representation of the Lie group G. The Lie algebra g admits the Casimir function K(f)= f f 2 + f f f f 2f . Expansion (4.2) takes 5 1 2 4 1 − 2 3 28 the form

M = f g∗ (f = f = 0) , 0 = f M K (f)= ω0 , 0 { ∈ | ¬ 1 2 } Oω ∈ 0 | 1 1 Ω0 = R1, dim 0 =4,  Oλ ∗ 1 (1) 1 M = f g f = f =0, (f = 0) , 1 = f M K (f)= ω , 1 { ∈ | 1 2 ¬ } Oω ∈ 0 | 1 1 n o M = 2 = f =0 , 2 O { } ∗ 0 1 g = M0 M1 M2, M0 = ω, M1 = ω1 , ∪ ∪ 0 0 O 1 1 O ω[∈Ω ω[∈Ω (1) 2 where K1 (f) = f4 +4f3f5 is the Casimir function of s = 1 type. Each non-degenerate orbit from the class 0 passes through the covector λ(j) = (1, 0, 0, 0, j) and K (λ(j)) = j, Oω 1 where j R. ∈ The covector λ(j) admits a real polarization p = e , e , e , and the λ-representation { 1 2 5} corresponding to the class of orbits 0 is given by Oλ(j) i i ℓ (q, ∂ ,λ)= q , ℓ (q, ∂ ,λ)= q , 1 q ~ 1 2 q −~ 2

ℓ (q, ∂ ,λ)= q ∂ 2 , ℓ (q, ∂ ,λ)= q ∂ 1 q ∂ 2 , 3 q 1 q 4 q 1 q − 2 q i j ℓ (q, ∂ ,λ)= q ∂ 1 + , 5 q 2 q ~ 2 q1 K ( i~ℓ)= j. (7.4) 1 − The operators i~ℓ (q, ∂ ,λ) are symmetric with respect to the measure dµ(q) = dq dq , − X q 1 2 Q = R2. Now solving equations (4.10), we get

x4 λ −1 i 2 1 j e 1 D ′ (g ) = exp q x q x + qq ~ 2 1 ′ ′ " − q2 q1 − q1 !#! ×

4 4 δ h q′ + q e−x , q′ δ ex q x3 q , q′ , × − 1 2 1 1 − 1 − 2 2   
 where δ(q, q′) is the generalized Dirac delta function. The completeness and orthogonality conditions (4.15)-(4.14) are satisfied for the measure

dµ(λ)=(2π~)−3dj, j R1. ∈ Orbits from the class 0 satisfy the integral condition (4.13). Oλ(j) The homogeneous space M has zero index, iM = 0, and does not have identities that have to be taken into account in the method of noncommutative integration. So, the λ- representation (7.4) corresponds to the homogeneous space M.

29 Integrating the equation

′ ′ (ℓ5(q , ∂q′ ,λ)+Λ) ψσ(q )=0, σ =(q1, q2, j)

yields ij q′ ψ (q′) = exp 1 Λ ψ (q′ ). (7.5) σ ~q′ q′ − q′ σ 2  1 2 2  Substituting (7.5) into the Dirac equation in the λ-representation (6.2), we obtain the

ordinary differential equation for the spinor ψσ,

~ d ′ ′ ′ i ′ ψσ(q2)= M(q2)ψσ(q2), dq2 1 j M(q′ )= γˆ4 i~Γ+ q′ γˆ2 + γˆ3 m . 2 c q′ 2 q′2 − 1 2  2  Then we obtain the solution as ′ i q2 ψ (q′ )= C exp M(u)du , (7.6) σ 2 σ −~ Z1 ! where Cσ is the normalization factor. The solution to the Dirac equation in local coordinates can be obtained by substituting ′ (7.5) into (6.1) and integrating over q for h1 = 0:

ij −2x4 ~ q1e Λ i 4 ψ (x) = exp q1 − exp q x2 q x1 ψ ex q q x3 . (7.7) σ q q x3 ~ 2 − 1 σ 2 − 1 2 − 1 !     
 The equations (7.6) and (7.7) provide a complete and orthogonal set of solutions to the Dirac equation on the homogeneous space M with the metric (7.2).

VIII. THE DIRAC EQUATION IN AdS3 SPACE

Consider a three-dimensional de Sitter space M as a homogeneous space with the de Sitter group of transformations G = SO(1, 3) and the isotropy subgroup H being the Lorentz group SO(1, 2). The de Sitter space M has constant positive curvature and M is topologically isomorphic to R1 S2. × The group SO(1, 3) is defined as a rotation group of a 4-dimensional pseudo-Euclidean space with the metric (G ) = diag(1, 1, 1, 1). The Lie algebra g = so(1, 3) of the AB − − − group SO(1, 3) in the basis E A < B is defined by the commutation relations { AB | } [E , E ]= G E G E + G E G E , AB CD AD BC − AC BD BC AD − BD AC

30 where A, B, C, D =1,..., 4. The basis EAB can be represented as

Eab = eab, (a < b), Ea4 = ea/ε,

2 [ea, eb]= ε eab, a, b =1,..., 3.

Here the basis eab forms the isotropy subalgebra h = so(1, 2), and ε is a parameter defining the curvature of the de Sitter space. We introduce the canonical coordinates of the second kind on the Lie group G so that

h3e23 h2e13 h1e12 ye3 xe2 te1 g(t, x, y, h1, h2, h3)= e e e e e e , (8.1)

The group generators in the canonical coordinates can be written as

X1(x)= ∂t,

X2(x) = sinh(εt) tan(εx)∂t + cosh(εt)∂x,

X3(x) = sinh(εt) sec(εx) tan(εy)∂t+

+ cosh(εt) sin(εx) tan(εy)∂x+

+ cosh(εt) cos(εx)∂y,

−1 X12(x)= ε (cosh(εt) tan(εx)∂t + sinh(εt)∂x) ,

−1 X13(x)= ε cosh(εt) sec(εx) tan(εy)∂t+

+ sinh( εt) sin(εx) tan(εy)∂x+

+ sinh(εt) cos(εx)∂y , X (x)= ε−1 (cos(εx) tan(εy
)∂ sin(εx)∂ ) . 23 x − y

The vector fields ξ(x, eH), in turn, are determined by

ξ1(x, eH )= X1(x),

ξ2(x, eH )= X2(x)+ ǫ sinh(tǫ) sec(xǫ)∂h1 ,

ξ3(x, eH )= X3(x)+ ǫ sinh(tǫ) tan(xǫ) tan(yǫ)∂h1

+ ǫ sinh(tǫ) sec(yǫ)∂h2 +

+ ǫ cosh(tǫ) sin(xǫ) sec(yǫ)∂h3 ,

31 ξ12(x, eH )= X12(x) + cosh(tǫ) sec(xǫ)∂h1 ,

ξ13(x, eH )= X13(x) + cosh(tǫ) tan(xǫ) tan(yǫ)∂h1

+ cosh(tǫ) sec(yǫ)∂h2 +

+ sinh(tǫ) sin(xǫ) sec(yǫ)∂h3 ,

ξ14(x, eH )= X14(x) + cos(xǫ) sec(yǫ)∂h3 .

The right-invariant vector fields in the canonical coordinates (8.1) are

η (x, e )= sec(εx) sec(εy)∂ 1 H − x− ε tan(εx) sec(εy)∂ 1 ε tan(εy)∂ 2 , − h − h η (x, e )= sec(εy)∂ ε tan(εy)∂ 3 , 2 H − y − h η (x, e )= ∂ , η (x, e )= ∂ 1 , 3 H − z 12 H − h η (x, e )= ∂ 2 , η (x, e )= ∂ 3 . 13 H − h 23 H − h The 2–form (G ) = diag(1, 1, 1) defines an invariant metric on the space M in local ab − − coordinates as dx2 dy2 ds2 = g (x)dxidxj = ρ(x, y)2 dt2 , ij − cos2(εx) − ρ(x, y)2   ρ(x, y) = cos(εx) cos(εy).

Whence the scalar curvature of the space M reads R =6ε2. The gamma matricesγ ˆa in the (2 + 1)-dimensional space in terms of the Pauli matrices σ , σ , σ are { x y z} 1 2 3 γˆ = σz, γˆ = isσx, γˆ = iσy.

Here the parameter s = Trˆγ1γˆ2γˆ3 = 1 is called a pseudospin. The matrices ± is is is Λ = γˆ3, Λ = γˆ2, Λ = γˆ1 12 − 2 13 2 23 − 2 realize a spinor representation of the isotropy subalgebra in a space of two-dimensional

spinors. In our case Γa = 0 and the spin connection takes the simple form ε Γ(x) =ˆγaηα(x, e )Λ = γˆ2 tan(εx) sec(εy)+ εγˆ3 tan(εy). a H α 2 Then the Dirac operator in local coordinates reads ε = i sec(εx) sec(εy)ˆγ1∂ i sec(εy)ˆγ2 ∂ tan(εx) iγˆ3 (∂ ε tan(εy)) . DM − t − x − 2 − y −   32 The first-order symmetry operators, as defined by (5.22), are given by

X1(x)= X1(x), isε X (x)= X (x) sinh(εt) sec(εx)ˆγ3, e2 2 − 2 isε 1 Xe3(x)= X3(x) cosh(εt) sin(εx) sec(εy)ˆγ − 2 − sinh(εt) sec(εy )ˆγ2 e − − sinh(εt) tan(εx) tan(εy)ˆγ3 , − isε X (x)= X (x) cosh(εt) sec(εx
)ˆγ3, 12 12 − 2 is 1 Xe13(x)= X13(x) sinh(εt) sin(εx) sec(εy)ˆγ − 2 − cosh(εt) sec(εy )ˆγ2+ e − + cosh(εt) tan(εx) tan(εy)ˆγ3 , isε X (x)= X (x) cos(εx) sec(εy
)ˆγ1. 14 14 − 2

Our aim is to constructe a complete set of exact solutions of the Dirac equation corresponding to the operator using the noncommutative integration method. DM The Lie algebra g admits the following two Casimir functions:

ε2 ε K (f)= pap labl , K (f)= eABCDL L , 1 a − 2 ab 2 8 AB CD

ABCD 1234 where e is the Levi-Civita symbol (e = 1); (pa,lab) are coordinates of the covector a ab a ab with respect to the basis e , e , i.e. f = pae + labe , and LAB are the coordinates with AB AB respect to the basis E  , f = LABE . The raising and lowering indices are performed

using the matrix GAB. The expansion (4.2) in our case takes the form:

M = f g∗ (f = 0) , 0 { ∈ | ¬ } 0 = f M K (f)= ω0, K (f)= ω0 , Oω ∈ 0 | 1 1 2 2 Ω0 = R2, dim 0 =4, Oλ M = , 1 = , M = 2 = f =0 , 1 {∅} O {∅} 2 O { } ∗ 0 g = M0 M2, M0 = ω. ∪ 0 0 O ω[∈Ω Each non-degenerate orbit from the class 0 passes through the covector λ(j)=(j , 0, 0, 0, 0, j ) Oω 1 2 33 characterized by two real parameters j =(j , j ) R2, and 1 2 ∈ K (λ(j)) = ω0(j)= j2 ε2j2, K (λ(j)) = ω0(j)= j j , 1 1 1 − 2 2 2 1 2 ∂ω0 (j) det µ =2 j2 + ε2j2 . ∂j 1 2 ν
If K (f)=0 for f p⊥ = f ea f g∗ , then the Casimir operator K (X) is an identity 2 ∈ { a | ∈ } 2 on the homogeneous space M:

K (X)= X X X X + X X =0. 2 3 · 12 − 2 · 13 1 · 23

Note that K2(X) is proportional to the Dirac operator: s e K (X)= X X X X + X X = . 2 3 · 12 − 2 · 13 1 · 23 2DM The covector λ(j) admitse the reale polarizatione e e e e

p = e , e + εe , e + εe , e . { 1 2 12 3 13 23} The corresponding λ-representation for the class of orbits 0 is represented in Appendix Oλ(j) A (see Eqs. (9.1)). The Casimir operators in the λ-representation are

K ( i~ℓ)= κ (λ)= j2 ε2j2 +(ε~)2, 1 − 2 1 − 2 K ( i~ℓ)= κ (λ)= j j . 2 − 2 1 2

′ ′ The equation (ℓα(q , ∂q′ ,λ)+Λα) cλ(q ) = 0 provided that j2 = s/2 has a nonzero solution

3 ij1 − − ~ (cos(εq1) cos(εq2)) 2 ε cλ(q)= cos(εq1) cos(εq2)+1 ×

ε(q + q ) ε(q q ) cos(εq1) cos(εq2)+1 cos 1 2 + is cos 1 − 2 . × 2 2 is sin(εq ) cos(εq ) sin(εq )       1 − 1 2   The Dirac equation in the λ-representation,

a ′ ′ i~γˆ ℓ (q , ∂ ′ ,λ) m c (q )=0, { a q − } λ is reduced to the algebraic equation j + m = 0, then we have j = m and j = s/2. That 1 1 − 2 is, the eigenvalue of the Casimir operator K1(i~X) is determined by the particle mass m, and the eigenvalue of the second Casimir operator, K2(i~X), depends on the parameter s: 1 1 κ (λ)= m2 ε2 +(ε~)2, κ (λ)= ms. 1 − 4 2 −2

34 The solution of the original Dirac equation in our case reads

−tℓ1(q,∂q ,λ) −xℓ2(q,∂q ,λ) −yℓ3(q,∂q.λ) ψσ(x)= e e e cλ(q), σ =(q1, q2) (8.2)

Here, the exponentials of operators of the λ-representation for the fixed j = m and 1 − j2 = s/2 act on a function according to Appendix A, Eqs. (9.2). From here one can see that the solution (8.2) depends on two quantum numbers q1 and q2, which are not eigenvalues for symmetry operators. The explicit form of the solution (8.2) is cumbersome, but it is expressed in terms of elementary functions.

IX. CONCLUSION

In this paper, we have explored the Dirac equation with an invariant metric on a ho- mogeneous space M of arbitrary dimension and developed the noncommutative integration method for this equation based on the ideas of symmetry analysis and the Lie group theory. The Dirac equation and its symmetry are convenient to study in terms of algebraic structures associated with homogeneous spaces, and the theory of Lie group representations can be effectively applied for constructing exact solutions. Using a special choice of the local frame and right-invariant vector fields on the Lie group of transformations G, we have obtained the spin connection (5.7). The Dirac equation is shown to be equivalent to a system of linear differential equations with constant matrix coefficients on the Lie group G given by (5.23) that is the starting point for noncommutative integration. In Refs. [22–27] an early version of the NCIM was used to construct exact solutions to the Dirac equation in four-dimensional space-time where the λ-representation was constructed directly by definition (4.7), and the domain of the variables q was not associated with a Lagrangian submanifold to the K-orbit. The desired solutions were constructed by means of joint integration of the system (6.4) in local coordinates together with the original Dirac equation (5.1). Differently to the above early method, the main idea here is the noncommutative re- duction of the corresponding system of equations on the Lie group G and the connection between the solutions of this system and the original Dirac equation. The noncommutative reduction is defined here using a special irreducible λ-representation

35 of the Lie algebra g of the Lie group G, which we introduce using the orbit method [37]. The key point of the method developed is based on the fact that there exist the identities connecting symmetry operators on a homogeneous space. For the Dirac equation, as follows from the lemma 3, the number of identities is either less than for the Klein–Gordon equation or they are completely absent. For the Klein–Gordon equation, the number of identities is determined by the index of the homogeneous space [19]. The problem of describing identities for the symmetry operators of the Dirac equation on the homogeneous space is constructively solved for the first time. The reduced system (6.2) in the λ-representation depending on a smaller number of independent variables q′ is obtained. What is remarkable is the fact that the solutions obtained for the Dirac equation on a homogeneous space are closely related to the two principal bundles of the transformation group. The local coordinates x on the space M are determined by means of the principal H-bundle of the group G with the isotropy subgroup H of M, and the set of quantum numbers q is connected with the principal P -bundle of the same group G and the subgroup P . The NCIM developed in the paper for the Dirac equation is illustrated by two non-trivial examples. In one of them, described in Section VII, we have found by using the NCIM a complete set of solutions to the Dirac equation (the MCIM-solutions) in the case when the metric does not admit separation of variables neither in the Klein-Gordon equation nor in the Dirac equation. The solutions obtained are eigenfunctions of the symmetry operator of the third-order (7.3) and are parameterized by three parameters (q1, q2, j).

The second example is the three-dimensional de Sitter space AdS3 with the transformation group SO(1, 3) (Section VIII). In this case, the Casimir operator K ( i~X) is proportional 2 − to the Dirac operator: e =2sK( i~X), DM − and the spectrum of the Casimir operator K ( i~X) gives the mass m of the spinor field. − 1 − e Noteworthy, the function K (f) is an identity in the homogeneous space M, but when 2 e substituting the extended operators i~X, it is no longer an identity, and this leads to − the original Dirac operator. If we consider the Klein–Gordon equation in AdS , the oper- e 3 ator K ( i~X) is proportional to the operator of the equation, and the second operator − 1 − K ( i~X) is identically zero that corresponds to the case s = 0. So the noncommutative 2 − e integration of the Dirac equation is different from the noncommutative integration of the

36 Klein-Gordon equation. The complete set of exact solutions (8.2) of the Dirac equation

found by the NCIM is parameterized by two real parameters (q1, q2) and is expressed by means of elementary functions, while the separation of variables gives the basis solutions to the Dirac equation in terms of special functions. The parameters q of solutions (7.7) and (8.2) obtained by the NCIM are in general not eigenvalues of an operator, a fact that crucially distinguishes them from solutions obtained by separation of variables. Nevertheless, the NCIM-solutions can be effectively applied in order to study quantum effects in homogeneous spaces (see, e.g., [19, 20]). The NCIM-solutions of the Dirac equation may have a wide range of applications in the theory of fermion fields [45, 46], quantum cosmology [47, 48] and other problems of field theory. The NCIM can be applied also to the Dirac-type equation for theoretical models in the condensed matter (graphene, topological insulators, etc.) [49, 50]. Note that the technique proposed in the article can be easily generalized to the case of spaces having new spatial dimensions much larger than the weak scale, as large as a millimeter for the case of two extra dimensions [51]. Finally, we note that the NCIM reveals new aspects, both related to the symmetry of the Dirac equation and its integrability, and to study the properties of new solutions constructed. One of the problems is to find out the meaning of the parameters q entering into the NCIM-solutions which, in the general case, do not have to be eigenvalues of operators representing observables. One can notice some similarity of the NCIM-solutions with well- studied coherent states [53]. In particular, the action of the group on the set of Q data of quantum numbers is defined, that can be found in the theory of coherent states [52–55]. However, the analysis of the parameters is the subject of special research.

ACKNOWLEDGEMENTS

Breev and Shapovalov were partially supported by Tomsk State University under the International Competitiveness Improvement Program; Breev was partially supported by the Russian Foundation for Basic Research (RFBR) under the project No. 18-02- 00149; Shapovalov was partially supported by Tomsk Polytechnic University under the International Competitiveness Improvement Program and by RFBR and Tomsk region according to the research project No. 19-41-700004.

37 APPENDIX A. λ-REPRESENTATION OF LIE ALGEBRA so(1, 3)

The λ -representation for the class of orbits 0 is written as Oλ(j)

ℓ1(q, ∂q,λ) = sin(εq1) cos(εq2)∂q1 + sec(εq1) sin(εq2)∂q2 + i i + ε + ~j1 cos(εq1) cos(εq2)+ ~εj2 tan(εq1) sin(εq2),   i ℓ (q, ∂ ,λ) =cos(εq )∂ 1 + tan(εq ) sin(εq )∂ 2 εj sec(εq ) sin(εq ), 2 q 2 q 1 2 q − ~ 2 1 2

ℓ3(q, ∂q,λ)=∂q2 ,

−1 i ℓ (q, ∂ ,λ)=ε cos(εq )∂ 1 + ε + j sin(εq ) , 12 q − 1 q ~ 1 1    

−1 ℓ (q, ∂ ,λ)=ε sin(εq ) sin(εq )∂ 1 sec(εq ) cos(εq )∂ 2 + 13 q 1 2 q − 1 2 q  i i + ε + ~j1 cos(εq1) sin(εq2)+ ~εj2 tan(εq1) cos(εq2) ,   

−1 i ℓ (q, ∂ ,λ)=ε sin(εq )∂ 1 tan(εq ) cos(εq )∂ 2 + εj sec(εq ) cos(εq ) . (9.1) 23 q 2 q − 1 2 q ~ 2 1 2   The exponentials of operators ℓ (q, ∂ ,λ), ℓ (q, ∂ ,λ) and ℓ (q, ∂ ,λ) for the fixed j = m 1 q 2 q 3 q 1 − and j2 = s/2 act on a function as follows:

e−yℓ3(q,∂q,λ)f(q)= f(q , q y), 1 2 −

is cot εq −xℓ (q,∂ ,λ) arccot 2 tan εq1 e 2 q f(q)= e 2  sin εq1 Φ cos εq sin εq , arctan εx , 2 1 2 cos εq −   2   2 2 2 − 1 is arctan(a tan b) −1 √1 a tan b −1 a tan b +1 Φ2(a, b)= e 2 f ε arctan − ,ε arcsec , √a2 tan2 b +1 s 1 a2     −    (sin εq )im/ε e−tℓ1(q,∂q ,λ)f(q)= 1 sin εq (sin εq cos εq is sin εq )× 1 1 2 − 2

p cos εq1 cos εq2 sinh εt + cosh εt Φ3 sin εq2 cot εq1, ,  2  × 1 (cos εq cos εq ) − 1 2  q  1/2 im/ε √ 2 2 b 1 is ab Φ3(a, b)= b√a +1 − − b2√a2 +1 (a2 + 1) b2 1! ×   − √a2 +1b p √a2 +1√b2 1 f ε−1arcsec ,ε−1arccos − . (9.2) (a2 + 1) b2 1! (a2 + 1) b2 1!  − −  p p 38 [1] Birrell, N.; Davies P. Quantum Fields in Curved Space; Cambridge Univ. Press: Cambrigde, 1986; p. 340 [2] Grib, A; Mostepanenko, V; Mamayev, S. Vacuum quantum effects in strong fields; Fridmann Lab: St. Petersburg, Russia, 1994; 361 p. [3] Ford, L. Vacuum polarization in a nonsimply connected spacetime. Phys. Rev. D 1980, 21, 933 [4] Frolov V. P.; Zel’Nikov, A. I. Vacuum polarization of massive fields near rotating black holes. Phys. Rev. D 1984, 29, 1057 [5] Wald, R. M. The back reaction effect in particle creation in curved spacetime. Comm. in Math. Phys. 1977, 54 , 1–19 [6] Kadoyoshi, T.; Sugamoto, A.; Nojiri, S. I.; Odintsov, S. D. Vacuum polarization of supersym- metric D-brane in the constant electromagnetic field. Mod. Phys. Lett. A 1998, 13, 1531–1537 [7] Kalnins, E. Separation of Variables in Riemannian Spaces of Constant Curvature; Wiley: New York, USA, 1986; 172 p. [8] Kalnins, E.G.; Miller, W. Jr.; Williams, G.C.; Recent advances in the use of separation of variables methods in general relativity. Philos. Trans. Roy. Soc. London Ser. A 1992, 340, 337—352 [9] Miller, W. Symmetry and Separation of Variables; Cambridge University Press: Cambridge, UK, 1984; p. 318. [10] Bagrov, V.; Gitman, D. Exact solutions of relativistic wave equations; Kluwer Academic Pub- lishers: Dordrecht, Netherlands, 1990; 324 p. [11] Stephani, H.; Kramer, D.; MacCallum, M.; Hoenselaers, C; Herlt, E. Exact solutions of Ein- stein’s field equations, 2th ed.; Cambridge University Press, Camridge, UK, 2003; 732 p. [12] Ryan, M.; Lawrence, C. Homogeneous relativistic cosmologies, reprint ed.; Princeton Univer- sity Press, New Jersey, USA, 2015; 338 p. [13] Shapovalov, V. N. Stackel spaces. Siberian Mathematical Journal 1979, 20, 790–800 [14] Shapovalov, V. N. Symmetry and separation of variables in Hamilton-Jacobi equations. I. Sov. Phys. J. 1978, 21, 1124–1129

39 [15] Obukhov, V. Hamilton–Jacobi Equation for a Charged Test Particle in the St¨ackel Space of Type (2.0). Symmetry 2020, 12, 1289. [16] Shapovalov, A.V.; Shirokov, I.V. Noncommutative integration of linear differential equations. Theor. Math. Phys. 1995, 104, 921-–934 [17] Shirokov, I.V. Symmetry in Nonlinear Mathematical Physics. Proceedings of Institute of Math- ematics of NAS of Ukraine 2003, 50, 246—251 [18] Baranovskii, S.P. ; Mikheev, V. V.; Shirokov, I. V. Quantum Hamiltonian systems on K-orbits: semiclassical spectrum of the asymmetric top. Theor. Math. Phys. 2001, 129, 1311–1319 [19] Breev, A.I.; Magazev, A.A.; Shirokov, I.V. Vacuum polarization of a scalar field on Lie groups and homogeneous spaces. Theor. Math. Phys. 2011, 167, 468—483 [20] Breev, A.I. Scalar field vacuum polarization on homogeneous spaces with an invariant metric. Theor. Math. Phys. 2014, 178, 59-–75 [21] Breev, A. I.; Kozlov, A. V. Vacuum Averages of the Energy-Momentum Tensor of a Scalar Field in Homogeneous Spaces with a Conformal Metric. Russ. Phys. J. 2016,58, 1248–1257 [22] Fedoseev, V. G.; Shapovalov, A. V.; Shirokov, I. V. On noncommutative solution of the Dirac equation in Riemann space with a dynamical group. Izv. Vuz. Phys. 1991, 34, 43–46 [23] Shapovalov, A. V.; Shirokov, I. V. Noncommutative integration of Klein-Gordon and Dirac equations with movement group. Izv. Vuz. Phys. 1991, 34, 33–38 [24] Varaksin, O. L.; Shirokov, I. V. Integration of the Dirac equation, which does not presume complete separation of variables, in St¨ackel spaces. Russ. Phys. J. 1996, 39, 27–32 [25] Varaksin, O. L.; Klishevich, V. V. Integration of Dirac equation in Riemannian spaces with five-dimensional group of motions. Russ. Phys. J. 1997, 40, 727–731 [26] Klishevich, V. V. Integration of the Dirac equation in Riemannian space with group of motions. I. Russ. Phys. J. 2000, 43, 1038–1043 [27] Klishevich, V. V. Exact solution of Dirac and Klein-Gordon-Fock equations in a curved space admitting a second Dirac operator. Class. and Quantum Grav. 2001, 18, 3735 [28] Tyumentsev, V. A.; Klishevich, V. V. Noncommutative Integration of the Dirac Equation in a Flat Space and in the de Sitter Space. Russ. Phys. J. 2003, 46, 891–896 [29] Klishevich, V. V.; Tyumentsev, V. A. On the solution of the Dirac equation in de Sitter space. Class. and Quantum Grav. 2005, 22, 4263

40 [30] Breev, A. I.; Shirokov, I. V. Polarization of a spinor field vacuum on manifolds of the Lie groups Russ. Phys. J. 2009, 52, 823–832 [31] Breev, A. I.; Magazev, A. A. Integration of the Dirac equation on Lie groups in an external electromagnetic field admitting a noncommutative symmetry algebra. Russ. Phys. J. 2017, 59, 2048–2058 [32] Dubrovin, B.; Fomenko, A.; Novikov, S.Modern Geometry and Applications, Part III: Intro- duction to Homology Theory; Graduate Texts in Mathematics; Springer Verlag, New York, USA, 1990; 418 p. [33] Kobayashi, S.; Nomizu, K. Foundations of Differential Geometry, Vol. I; Wiley-Interscience: New York, USA, 1996; 348 p. [34] Arvanitoge¯orgos, A. An Introduction to Lie Groups and the Geometry of Homogeneous Spaces; Student Mathematical Library, V. 22; American Mathematical Society: Rhode Island, USA, 2003; 141 p. [35] Kurnyavko, O. L.; Shirokov, I.V. Construction of invariant scalar particle wave equations on Riemannian manifolds with external gauge fields. Theor. Math. Phys. 2008, 156, 1169—1179 [36] Baranovskiii, S.P.; Shirokov, I.V. Prolongations of vector fields on Lie groups and homogeneous spaces. Theor. Math. Phys. 2003, 135, 510-–519 [37] Kirillov, A. Lectures on the Orbit Method; Graduate Studies in Mathematics V. 64; American Mathematical Society: Rhode Island, USA, 2004; 408 p. [38] Magazev, A. A.; Mikheyev, V. V.; Shirokov, I. V. Computation of composition functions and invariant vector fields in terms of structure constants of associated Lie algebras. SIGMA. 2015, 11, 066 [39] Shirokov, I. V. Darboux coordinates on K-orbits and the spectra of Casimir operators on Lie groups. Theor. Math. Phys. 2000, 123, 754—767 [40] Shirokov, I. V. Identities and invariant operators on homogeneous spaces. Theor. Math. Phys. 2001, 126, 326-–338 [41] Dixmier, J. Enveloping algebras; North-Holland Mathematical Studies Vol. 14; Elsevier Science & Technology Books, 1977; 376 p. [42] Barut, A; Raczka, R. Theory of Group Representations and Applications, 2th ed.; World Scientific Publishing Company: London, UK, 1986; 740 p.

41 [43] Bagrov, V.G.; Shapovalov, A.V.; Yevseyevich, A.A. Separation of variables in the Dirac equa- tion in Stackel spaces. II. External gauge fields. Class. Quantum Grav. 1991, 8, 163-–173 [44] Breev, A. I.; Goncharovskii, M. M.; Shirokov, I. V. Klein-Gordon equation with a special type of nonlocal nonlinearity in commutative homogeneous spaces with invariant metric. Russ. Phys. J. 2013, 56, 731–739 [45] Hack, T. Cosmological applications of algebraic quantum field theory in curved spacetimes, 1st ed.; SpringerBriefs in Mathematical Physics V. 6; Springer: New York, USA, 2016; 131 p. [46] Toms, D. J. Effective action for the Yukawa model in curved spacetime. J. High Energ. Phys. 2018, 5, 139 [47] Nojiri, S. I.; Odintsov, S. D. Effective equation of state and energy conditions in phan- tom/tachyon inflationary cosmology perturbed by quantum effects. Phys. Lett. B 2003, 571, 1–10 [48] Brevik, I.; Milton, K. A.; Odintsov, S. D.; Osetrin, K. E. Dynamical Casimir effect and quantum cosmology. Phys. Rev. D 2000, 62, 064005 [49] Vafek, O.; Vishwanath, A. Dirac Fermions in Solids - from High Tc cuprates and Graphene to Topological Insulators and Weyl Semimetals. ARCMP 2014, 5, 83–112 [50] Klimchitskaya, G. L.; Mostepanenko, V. M. Creation of quasiparticles in graphene by a time- dependent electric feld. Phys. Rev. D 2013, 87, 125011 [51] Antoniadis, I.; Arkani-Hamed, N. ;Dimopoulos, S.; Dvali, G. New dimensions at a millimeter to a Fermi and superstrings at a TeV. Phys. Lett. B 1998, 436, 257–263 [52] Malkin, I.; Manko, V. Dynamic symmetry and coherent states of quantum systems; Nauka: Moscow, Nauka, 1979; [53] Perelomov, A. Generalized coherent states and their applications; Theoretical and Mathemat- ical Physics; Springer: Netherlands, (2012); 418 p. [54] Perelomov, A. M. Coherent states for arbitrary Lie group. Commun. Math. Phys 1972, 26, 222–236 [55] Bagrov, V.; Gitman, D. The Dirac equation and its solutions; Walter de Gruyter GmbH: Brelin/Boston, 2014; 444 p.

42