Abstract

:
We consider a (2 + 1)-dimensional massless Dirac equation in the presence of complex vector potentials. It is shown that such vector potentials (leading to complex magnetic fields) can produce bound states, and the Dirac Hamiltonians are η-pseudo Hermitian. Some examples have been explicitly worked out.

1. Introduction

In recent years, the massless Dirac equation in (2 + 1) dimensions has drawn a lot of attention, primarily because of its similarity to the equation governing the motion of charge carriers in graphene [1,2]. In view of the fact that electrostatic fields alone cannot provide confinement of the electrons, there have been quite a number of works on exact solutions of the relevant Dirac equation with different magnetic field configurations, for example, square well magnetic barriers [3–5], non-zero magnetic fields in dots [6], decaying magnetic fields [7], solvable magnetic field configurations [8], etc. On the other hand, at the same time, there have been some investigations into the possible role of non-Hermiticity and
symmetry [9] in graphene [10–12], optical analogues of relativistic quantum mechanics [13] and relativistic non-Hermitian quantum mechanics [14], photonic honeycomb lattice [15], etc. Furthermore, the (2 + 1)-dimensional Dirac equation with non-Hermitian Rashba and scalar interaction was studied [16]. Here, our objective is to widen the scope of incorporating non-Hermitian interactions in the (2 + 1)-dimensional Dirac equation. We shall introduce η pseudo Hermitian interactions by using imaginary vector potentials. It may be noted that imaginary vector potentials have been studied previously in connection with the localization/delocalization problem [17,18], as well as
phase transition in higher dimensions [19]. Furthermore, in the case of the Dirac equation, there are the possibilities of transforming real electric fields to complex magnetic fields and vice versa by the application of a complex Lorentz boost [20]. To be more specific, we shall consider η-pseudo Hermitian interactions [21] within the framework of the (2 + 1)-dimensional massless Dirac equation. In particular, we shall examine the exact bound state solutions in the presence of imaginary magnetic fields arising out of imaginary vector potentials. We shall also obtain the η operator, and it will be shown that the Dirac Hamiltonians are η-pseudo Hermitian.

2. The Model

The (2 + 1)-dimensional massless Dirac equation is given by:

Hψ=Eψ,H=cσ⋅P=c(0P−P+0),ψ=(ψ1ψ2)

(1)

where c is the velocity of light and:

P±=(Px±iPy)=(px+Ax)±i(py+Ay)

(2)

In order to solve Equation (1), it is necessary to decouple the spinor components. Applying the operator, H, from the left in Equation (1), we find:

c2(P−P+00P+P−)ψ=E2ψ

(3)

Let us now consider the vector potential to be:

Ax=0,Ay=f(x)

(4)

so that the magnetic field is given by:

Bz(x)=f′(x)

(5)

For the above choice of vector potentials, the component wave functions can be taken of the form:

ψ1,2(x,y)=eikyyϕ1,2(x)

(6)

Then, from (3), the equations for the components are found to be (in units of ℏ = 1):

[−d2dx2+W2(x)+W′(x)]ϕ1(x)=ε2ϕ1(x)[−d2dx2+W2(x)−W′(x)]ϕ2(x)=ε2ϕ2(x)

(7)

where ϵ = (E/c), and the function, W(x), is given by:

W(x)=ky+f(x)

(8)

2.1. Complex Decaying Magnetic Field

It is now necessary to choose the function, f(x). Our first choice for this function is:

f(x)=−(A+iB)e−x,−∞<x<∞

(9)

where A > 0 and B are constants. This leads to a complex exponentially decaying magnetic field:

Bz(x)=(A+iB)e−x

(10)

For B = 0 or a purely imaginary number (such that (A + iB) > 0), the magnetic field is an exponentially decreasing one, and we recover the case considered in [7,8].

It is not difficult to recognize V2(x) in Equation (12) as the complex analogue of the Morse potential whose solutions are well known [22,23]. Using these results, we find:

E2,n=±cky2−(ky−n)2ϕ2,n=tky−ne−t/2Ln(2ky−2n)(t),n=0,1,2,….<[ky]

(13)

where t = 2(A + iB)e−x and
Ln(a)(t) denote generalized Laguerre polynomials. The first point to note here is that for the energy levels to be real, it follows from Equation (13) that the corresponding eigenfunctions are normalizable when the condition ky ≥ 0 holds. For ky < 0, the wave functions are not normalizable, i.e., no bound states are possible.

Let us now examine the upper component, ϕ1. Since ϕ2 is known, one can always use the intertwining relation:

cP−ψ2=Eψ1

(14)

to obtain ϕ1. Nevertheless, for the sake of completeness, we present the explicit results for ϕ1. In this case, the potential analogous to Equation (12) reads:

V1(x)=ky2+(A+iB)2e−2x−(2ky−1)(A+iB)e−x

(15)

Clearly, V1(x) can be obtained from V2(x) by the replacement ky → ky − 1, and so, the solutions can be obtained from Equation (13) as:

E1,n=±cky2−(ky−n−1)2ϕ1,n=tky−n−1e−t/2Ln(2ky−2n−2)(t),n=1,2,….<[ky−1]

(16)

Note that the n = 0 state is missing from the spectrum Equation (16), so that it is a singlet state. Furthermore, E2,n+1 = E1,n, so that the ground state is a singlet, while the excited ones are doubly degenerate. Similarly, the negative energy states are also paired. In this connection, we would like to note that {H, σ3} = 0, and consequently, except for the ground state, there is particle hole symmetry. The wave functions for the holes are given by σ3ψ. The precise structure of the wave functions of the original Dirac equation are as follows (we present only the positive energy solutions):

E0=0,ψ0=(0ϕ2,0)En+1=cky2−(ky−n−1)2,ψn+1=(ϕ1,nϕ2,n+1),n=0,1,2,…

(17)

It is interesting to note that the spectrum does not depend on the magnetic field. Furthermore, the dispersion relation is no longer linear, as it should be in the presence of interactions. It is also easily checked that when the magnetic field is reversed, i.e., A → −A and B → −B with the simultaneous change of ky → −ky, the two potentials V1,2(x) = W(x) ± W′(x) go one into each other, V1(x) ↔ V2(x). Therefore, the solutions are correspondingly interchanged, ϕ1,n ↔ ϕ2,n and E1,n ↔ E2,n, but retain the same functional form as in Equations (13) and (16).

Therefore, we find that it is indeed possible to create bound states with an imaginary vector potential. We shall now demonstrate the above results for a second example.

2.2. Complex Hyperbolic Magnetic Field

Here, we choose f(x), which leads to an effective potential of the complex hyperbolic Rosen–Morse type:

3. η-Pseudo Hermiticity

where η is a Hermitian operator. It is known that eigenvalues of a η-pseudo Hermitian Hamiltonian are either all real or are complex conjugate pairs [21]. In view of the fact that in the present examples, the eigenvalues are all real, one is tempted to conclude that the interactions are η pseudo Hermitian. To this end, we first consider case 1, and following [26], let us consider the Hermitian operator:

η=e−θpx,θ=arctanBA

(25)

Then, it follows that:

ηcη−1=c,ηpxη−1=px,ηV(x)η−1=V(x+iθ)

(26)

We recall that in both the cases considered here, the Hamiltonian is of the form:

Next, to demonstrate the pseudo Hermiticity of the Dirac Hamiltonian Equation (27), let us consider the operator η′ = η ·
2, where
2 is the (2 × 2) unit matrix. Then, it can be shown that:

η′Hη′−1=H†

(32)

Thus, the Dirac Hamiltonian with a complex decaying magnetic field Equation (10) is η-pseudo Hermitian.

For the magnetic field given by Equation (19), the operator, η, can be found by using relations Equation (26). After a straightforward calculation, it can be shown that the η operator is given by:

η=e−2αpx

(33)

so that, in this second example, also, the Dirac Hamiltonian is η-pseudo Hermitian.

4. Conclusions

Here, we have studied the (2 + 1)-dimensional massless Dirac equation (we note that if a massive particle of mass m is considered, the energy spectrum in the first example would become
En=cky2+m2c2−(ky−n)2. Similar changes will occur in the second example, too). in the presence of complex magnetic fields, and it has been shown that such magnetic fields can create bound states. It has also been shown that Dirac Hamiltonians in the presence of such magnetic fields are η-pseudo Hermitian. We feel it would be of interest to study the generation of bound states using other types of magnetic fields, e.g., periodic magnetic fields.

One of us (P. R.) wishes to thank INFN Sezione di Perugia for supporting a visit during which part of this work was carried out. He would also like to thank the Physics Department of the University of Perugia for its hospitality.