Nonlinear conformal-degree-preserving Dirac equations in 2+1

0 downloads 0 Views 213KB Size Report
Abstract: We propose nonlinear Dirac equations where the conformal degree of the self- interaction terms are equal to that of the Dirac operator and the coupling ...
Nonlinear conformal-degree-preserving Dirac equations in 2+1 space-time A. D. Alhaidari(a), U. Al Khawaja(a,b) and Y. H. Sabbah(b)

(b)

(a) Saudi Center for Theoretical Physics, P. O. Box 32741, Jeddah 21438, Saudi Arabia Physics Department, United Arab Emirates University, P. O. Box 15551, Al-Ain, United Arab Emirates

Abstract: We propose nonlinear Dirac equations where the conformal degree of the selfinteraction terms are equal to that of the Dirac operator and the coupling parameters are dimensionless. As such, the massless equation is conformally invariant and preserves the conformal degree for both the linear and nonlinear components. In 2+1 space-time, we use these features to propose three- and four-parameter nonlinear Dirac equation models that might be useful for applications in 2D systems such as graphene sheets, ribbons and nanotubes. PACS: 05.45.-a, 11.10.Lm, 03.65.Pm, 03.65.Ca Keywords: nonlinear interaction, Dirac equation, conformal degree, graphene

1. Introduction: In the language of conformally invariant theories, the conformal degree of a given field is the extent to (or order with) which it scales with the dimension of space and time. In n  1 Minkowski space-time, the kinetic energy term in the spinor Lagrangian is i , where  



 n 0

    is the Dirac operator, {  } are the dimensionless Dirac gamma

matrices, and    † 0 . Since the corresponding part of the action, which is the spacetime integral

  i  d

n 1

x , is conformally invariant, then the conformal degree of

i is n  1 making the conformal degree of the spinor field  12 n . Recently, it was shown that a nonlinear conformally invariant massless Dirac equation in this space-time could be written, in the conventional relativistic units   c  1 , as [1,2] 1n i    k g k   k   k  , (1)   where  k  and  k  are elements from the space of bilinear forms of the Dirac gamma matrices and  g k  are dimensionless coupling parameters. The sum of nonlinear terms in Eq. (1) runs over all independent coupling modes whose type and number depend on the dimensionality of space-time. These terms are obtained by stationary variations of the spinor field in the action. The nonlinear components of this action, which produce the right-hand side of Eq. (1), are constructed by requiring that their conformal degree be equal to n  1 and by imposing relativistic invariance on the various types of coupling modes. These modes include scalar, vector, tensor, pseudo-scalar, etc. and could be written using objects such as the following S   , Scalar: (2a) A     , (2b) Vector: 1

Tensor:

T     ,

(2c)

Pseudo-Scalar:

W   5 ,

(2d)

Pseudo-Vector:

P   5  ,

(2e)

Pseudo-Tensor: Q   5  ,

(2f)

and so on, where  5  i 0 1... n , and   

i 4 [   

     ] . However, some of these

couplings may become redundant in certain space-time dimensions. For example, in 2+1 dimension,   µ     and  5 µ 1 , where   is the completely antisymmetric tensor in 2+1 space-time. Therefore, the tensor, pseudo-scalar, pseudo-vector, and pseudo-tensor couplings are redundant. Now, the conformal degree of the spinor field in 2+1 space-time is 1. Thus, the independent nonlinear terms in the interaction Lagrangian with conformal degree 3 are the following S 3 2   

32

and  A  A  

34

34

  (  )(  )  ,

(3)

Consequently, the nonlinear massive Dirac equation (1) becomes k 4 (4)  i   m    k 1 g k  k  k  , where m is the rest mass of the spinor particle. Now, for the massless case, this nonlinear Dirac equation is invariant under the action of the full conformal group [3]. However, if the fractional powers 3 2 and 3 4 in Eq. (3) are changed, which will destroy the overall conformal degree, then the symmetry of the massless equation is reduced and invariance is only under the smaller Weyl group [3].





2. The Models: The minimum dimensional representation of the Dirac gamma matrices in 2+1 space-time





 

is 22. These could be taken as:  0   3  10 01 ,  1  i  1  i 01 10 , and  2  i 2 

 01 10  . Consequently, Eq. (4) becomes the following nonlinear matrix equation in 2+1

space-time with the most general coupling m   S S  VV  x  i y  iUU  WW        i t        ,   m   S S  VV          x  i y  iUU  WW where

S , V , U , W 

(5)

are real dimensionless coupling parameters for the scalar

coupling S and the three-component vector coupling  A0 , A1 , A2   V ,U ,W  . The selfinteraction potential components are 2 2 S (t , x, y )         ,

(6a)

2 2 V (t , x, y )   0       ,

(6b)

U (t , x, y )   1  i ( *    * ) .

(6c)

W (t , x, y )   2   *     * . (6d) We may eliminate the U component of the vector potential in the Dirac equation (5) by the nonlinear local gauge transformation   e  iU  , where    U . Substituting x

2

these nonlinear terms into Eq. (5) and gauging away the U term, we obtain the following three-parameter equation     m    2   2    2       S V i t     * *     x  i y  W          

2

       2        

 x  i y  W  *     * 2

2

 m   S       V  

2

(7)

We could have also made the problem more general by adding external linear potentials to the various coupling modes. However, we limit the present treatment to selfinteraction. We could simplify Eq. (7) by proposing the following alternative fourparameter model that also preserves the conformal degree but violates Lorentz invariance     m            x  i y             i t        , (8a)     x  i y    *    * m                  and its “conjugate”     m            x  i y    *    *     i t        . (8b)     x  i y         m                  where   and   are four independent real dimensionless parameters. Violation of Lorentz invariance in nonlinear models might not be so dissolute. In fact, it was even suggested that at some fundamental level such violation could be tied to quantum nonlinearity [4]. Moreover, we may think of these as effective models that describe a given physical phenomenon rather than a fundamental interaction. On the other hand, the following three-parameter model is Lorentz invariant but does not preserve the conformal degree in 2+1 space-time [1] 2 2     m          x  i  y  W  *     *            , (9) i t  2 2      i     *    *   m            y W            x where the coupling parameters have length dimensions (i.e., their conformal degree is +1). This could have a negative impact on the renormalizability of the theory. The Dirac equation (5) is written in Cartesian coordinates, which is suitable for the treatment of problems with rectangular configuration. On the other hand, problems with circular configuration are described properly in cylindrical coordinates (r , ) . To simplify the writing of such an equation, we may follow a procedure described in Ref. [5] whereby a local 22 transformation (r , ) is applied to the Dirac equation such that the   cylindrical projection of the gamma matrices     rˆ and     ˆ are mapped into 

r

their canonical Cartesian representations  1 and  2 , respectively. That is,   r   1 

 1 and      1   2 . Moreover, we require that the transformed Dirac Hamiltonian matrix be Hermitian. One may refer to Ref. [5] for the details of the procedure where one i  obtains (r , )  re 2 3 and the nonlinear Dirac equation (7) takes the following form in cylindrical coordinates   2 2 2      m  S     V    r r    i t         r  ri   W *  *    r 

2

 r  ri   m 

3

S r

2

W

        (10)   2 2         

*  *

r 2

   

V r

where    , which has a conformal degree of  12 . Therefore, the two-component  i   i spinor wavefunction becomes  (t , r , )  1r e 2 3 (t , r , ) . That is,   re 2   and

the density of state  (t )    †  rdrd    †  drd . Additionally, the model given by Eq. (8a) and its conjugate in Eq. (8b) become   1  m        r  ri                  r r r      . (11a) i t                i   1    *    *  m      r        r         r r r     1   m        *   *       r  ri        r r r      . (11b) i t                i   1         m                  r r  r r r  

To be able to obtain stationary solutions of the form   (t , x, y)  ei t   ( x, y) , we need to make sure that with this substitution time becomes separable in the above models. It is evident that the three-parameter Lorentz invariant model of Eq. (7) and its cylindrical version of Eq. (10) are compatible with this requirement. On the other hand, the fourparameter Lorentz-symmetry-violating model given by Eqs. (8) or Eqs. (11) are not. Consequently, we modify this latter model to read as follows in Cartesian coordinates     m          x  i y               i t        . (12)     x  i y           m                  Whereas in cylindrical coordinates, this same model reads as follows   1   r  ri   m                        r r r      . (13) i t     1            i    m                       r r  

r

r

r



3. Sample Solutions: As illustration of the above findings, we consider two relevant physical configurations in which one of the space coordinates become cyclic (i.e., the interaction is independent of that coordinate): 1. In Cartesian coordinates, we consider y-independent interactions. That is, y becomes the cyclic coordinate and the two spinor components could be written as   (t , x, y)  ei ( ky  t )   ( x) , where the real parameter k is the wavenumber along the y-axis. Such solutions are relevant to physical configurations such as nano-ribbons and nano-strips. 2. In cylindrical coordinates, we consider cylindrically symmetric interactions. That is,  becomes the cyclic coordinate and the two-component spinor wave function becomes i  i   (t , r , )  1r e 2 3 ei  (t , r )  1r e 2 3 ei (  t )  (r ) , (14) 4

where  is the azimuthal wave number. These solutions are relevant to nanosheets in the form of disks and sectors. For nano-disks and nanotubes, where 0    2 , continuity dictates that  becomes the azimuthal quantum number    12 ,  23 ,  52 ,... . In Table 1, we list exact analytic solutions in cylindrical and Cartesian coordinates for special values of the model parameters. The importance of these solutions is not only in showing that the model can be solved analytically and exactly, but also in the diversity of physical configurations of the model that differ significantly according to the choice of values of those parameters. Therefore, in principle, there are many other analytic solutions beyond those obtained here.

4. Conclusion: In conclusion, we imposed the constraint that all components of the massless Dirac equation (linear and nonlinear) in 2+1 space-time have the same conformal degree of 1 and that the coupling parameters are dimensionless. We showed that these constraints together with Lorentz invariance lead to the nonlinear three-parameter Dirac equation model given by Eq. (7) where only the mass term violates conformal invariance. We simplified the model by breaking Lorentz invariance, which resulted in the fourparameter model (and its conjugate) shown in Eq. (8a) and Eq. (8b). Additionally, to make this latter model separable in time so that stationary solutions can be obtained, we modified it as shown in Eq. (12) and Eq. (13). Sample analytic solutions were obtained. The important issues of integrability, exact solvability and range of stability of solutions of these models are left for future studies.

References [1] [2] [3] [4] [5]

A. D. Alhaidari, “Nonlinear spin and pseudo-spin symmetric Dirac equations”, Int. J. Theor. Phys. 53, 685 (2014) U. Al Khawaja, “Exact localized and oscillatory solutions of the nonlinear spin and pseudospin symmetric Dirac equations”, Phys. Rev. A 90, 052105 (2014) W. I. Fushchich and W. M. Shtelen, “On some exact solutions of the nonlinear Dirac equation”, J. Phys. A 16, 271 (1983) R. Parwani, “An information-theoretic link between space-time symmetries and quantum linearity”, Ann. Phys. 315, 419 (2005) A. Jellal, A. D. Alhaidari, and H. Bahlouli, “Confined Dirac fermions in a constant magnetic field”, Phys. Rev. A 80, 012109 (2009) pp. 2-3

Table Caption: Table 1: Sample analytic solutions in cylindrical and Cartesian coordinates for the model equations shown in the first column and for various special values of the model parameters. 5

Table 1 Model 

Wavefunction Components 



0  0  0  0 

 





tan 2

 

cot 2

 

 

   

Eq. (7) 

Eq. (12) 





Amplitude Functions 

0  0  0 

  Eq. (13) 

Parameters 

 



1



1



3

0  0    0 

2√2 3√

  0  0  0 

3 1 3 1

6

     

 

2√2