Planar Dirac Equation with Radial Contact Potentials

1. Introduction

Much work has recently been done on the role of electric and magnetic fields for technological applications in the confinement of relativistic massive charged particles, especially in two-dimensional systems governed by the Dirac-Weyl equation [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16], although it is known that this is not suitable for massless two-dimensional particles due to the existence of a strong Klein tunneling effect [17,18]. Another situation where electric fields are able to confine particles is when we deal with a quantum dot with some kind of symmetry, typically radial, that allows for a well-defined angular momentum quantum number. Confinement of massive particles in two dimensions has also been studied to better understand new topological materials, such as graphene, silicene, germanene, stanene and phosphorene [19,20,21], where charges acquire effective mass from spin-orbit interaction and perpendicular electric fields [22,23]. An important consequence of massive particle confinement is that it allows to observe the phenomenon of atomic collapse, which is very difficult to detect in relativistic quantum mechanics but is accessible with Dirac materials [24,25,26,27] due to the existence of the Fermi velocity $v_F<<c$ [28]. The search for resonances is also an interesting topic, and in fact confinement in Dirac materials must necessarily include bound and resonance states [29].

Although it is clear that the aforementioned technological applications are extremely interesting, in this work our focus will be on the theoretical analysis of certain two-dimensional systems that are governed by the (2+1) massive Dirac equation subjected to a type of toy potentials, the singular contact interactions in two dimensions. It is well known that contact or point interactions in one-dimensional relativistic quantum mechanics have been used to model short-range potentials. They are singular interactions whose appropriate treatment needs the use of special mathematical methods in order to avoid ambiguities (see [30] and references quoted therein). The main goal of the present work is to adapt a recent distributional approach, originally developed for the one-dimensional case [31], to study, in a mathematically rigorous manner, radially symmetric two-dimensional contact interactions supported on a circumference of radius R. From the radial symmetry it is possible to reduce the problem to one dimension (the radial one), and apply in an almost straightforward way the distributional approach of [31] to a contact potential at $r=R$. In the actual one-dimensional case, the usual four parameters of the point interactions can be rewritten in terms of four physical potentials strengths, corresponding to a scalar, a pseudoscalar and the two componentes of a (1+1) Lorentz electromagnetic point potential [31]. As it will be seen here, in the two-dimensional case with radial symmetry and a singular potential concentrated on a circumference, those four parameters will correspond to the physical strengths of a scalar and the three components of a (2+1) electromagnetic Lorentz singular potential. After obtaining the general solutions for the Dirac equation with the most general contact interaction concentrated on the circumference, we will consider some special cases, namely a purely scalar potential, a purely electrostatic potential, a purely magnetic potential and two mixtures of a scalar and an electrostatic potential. In each of these cases we will systematically investigate the confining properties of the contact potential, understood broadly as the structure of the bound states and resonances, and the conditions under which the circular barrier becomes impenetrable.

The structure of the paper is as follows. In Section 2 we review the covariant form of the Dirac equation in (2+1) dimensions with the most general external regular potential with radial symmetry. In Section 3 we consider the (2+1) Dirac equation with the most general, time-independent, contact potential supported on a circumference of radius R. In Section 4 we determine the equations for the energy eigenstates of the previously derived Dirac radial equation, including bound states, scattering states, and resonances. These quantities are subsequently investigated for some special cases of physical interest in Section 5. The paper ends with a discussion and conclusions in Section 6.

2. Planar Dirac Equation with a Radial Regular Potential

In this section we briefly review the covariant form of the Dirac equation in (2+1) dimensions with the most general external regular potential1, which is formed by a scalar $B\left(x\right)$ and a vector potential $A_{\mu }\left(x\right)$, with $x=\left(x^0,x^1,x^2\right)$ and $\mu =0,1,2$ (we are using natural units, $\hslash =c=1$):

$$\left(i{\gamma }^{\mu }{\partial }_{\mu }-m\right)\Psi \left(x\right)=\mathcal{I}\left(x\right)\Psi \left(x\right),\mathrm{where}\mathcal{I}\left(x\right)=B\left(x\right)⊮+A_{\mu }\left(x\right){\gamma }^{\mu }.$$

(1)

For convenience, we will use the following representation of the Dirac matrices, in terms of the Pauli matrices:

$${\gamma }^0={\sigma }_3,{\gamma }^1={\sigma }_3{\sigma }_1\mathrm{and}{\gamma }^2={\sigma }_3{\sigma }_2.$$

(2)

We are interested in the case where in a particular reference frame (e.g. the laboratory frame) the potentials are independent of time, and the time-independent Dirac equation reads as follows:

$$\left({\gamma }^{0}E+i{\gamma }^k{\partial }_k-m\right)\psi \left(\mathbf{r}\right)=\mathcal{I}\left(\mathbf{r}\right)\psi \left(\mathbf{r}\right),k=1,2,$$

(3)

where $\mathbf{r}=(x,y)$ is the position vector of the Dirac particle in the cartesian coordinates of the laboratory frame, $\psi \left(\mathbf{r}\right)$ is the corresponding time independent Dirac spinor and

$$\mathcal{I}\left(\mathbf{r}\right)=B\left(\mathbf{r}\right)⊮+A_{\mu }\left(\mathbf{r}\right){\gamma }^{\mu },\mu =0,1,2.$$

(4)

Even if the potential in the Dirac equation (3) is $V\left(\mathbf{r}\right)={\gamma }^0\mathcal{I}\left(\mathbf{r}\right)$, we will sometimes abuse the language and will refer to $\mathcal{I}\left(\mathbf{r}\right)$ as the potential.

Since we are only interested in the case where the potentials have radial symmetry (in the laboratory reference frame), it is natural to use polar coordinates from now on. The space cartesian components $A_k\left(\mathbf{r}\right)$ ($k=1,2$) of the Lorentz vector potential transform to the polar components $A_r(r,\theta )$ and $A_{\theta }(r,\theta )$, with $r=\left|\mathbf{r}\right|$ and $\theta$ the polar angle, as:

$$A_r(r,\theta )=A_1(x,y)cos\theta +A_2(x,y)sin\theta ,A_{\theta }(r,\theta )=-A_1(x,y)sin\theta +A_2(x,y)cos\theta ,$$

(5)

where $A_1(x,y)=A_1(rcos\theta ,rsin\theta )$, $A_2(x,y)=A_2(rcos\theta ,rsin\theta )$, and, with a slight abuse of notation, $B(r,\theta )\equiv B(rcos\theta ,rsin\theta )$, $A_0(r,\theta )\equiv A_0(rcos\theta ,rsin\theta )$. By defining

$${\gamma }^r={\gamma }^{1}cos\theta +{\gamma }^{2}sin\theta =\left(\begin{array}{cc}0& e^{-i\theta }\\ -e^{i\theta }& 0\end{array}\right),{\gamma }^{\theta }=-{\gamma }^{1}sin\theta +{\gamma }^{2}cos\theta =\left(\begin{array}{cc}0& -i{e}^{-i\theta }\\ -i{e}^{i\theta }& 0\end{array}\right),$$

(6)

we can write (4) in polar coordinates as

$$\mathcal{I}(rcos\theta ,rsin\theta )=B(r,\theta )⊮+A_0(r,\theta ){\gamma }^0+A_r(r,\theta ){\gamma }^r+A_{\theta }(r,\theta ){\gamma }^{\theta }\equiv \mathcal{I}(r,\theta ).$$

(7)

The whole potential will have radial symmetry if all the above four potentials $B(r,\theta )$, $A_0(r,\theta )$, $A_r(r,\theta )$ and $A_{\theta }(r,\theta )$ depend only on r, and in this case we will denote them simply as $B\left(r\right),A_0\left(r\right),A_r\left(r\right)$ and $A_{\theta }\left(r\right)$.2 Thus, in polar coordinates the Dirac equation (3), with a radial potential (7), reads3

$$\begin{array}{ccc}\hfill \left({\gamma }^{0}E+i{\gamma }^r{\partial }_r+{\displaystyle \frac{i}{r}}{\gamma }^{\theta }{\partial }_{\theta }-m\right)\psi (r,\theta )& =& \mathcal{I}(r,\theta )\psi (r,\theta )\hfill \\ & & =\left[B\left(r\right)⊮+A_0\left(r\right){\gamma }^0+A_r\left(r\right){\gamma }^r+A_{\theta }\left(r\right){\gamma }^{\theta }\right]\psi (r,\theta ),r>0.\hfill \end{array}$$

(8)

Taking advantage of the radial symmetry, we seek for solutions of the time independent Dirac equation that are also eigensolutions of the total angular momentum operator:

$$J_z=L_z+\Sigma ,\mathrm{where}{L}_z=-i{\partial }_{\theta }\mathrm{and}\Sigma ={\displaystyle \frac{1}{2}}{\sigma }_z={\displaystyle \frac{1}{2}}{\gamma }^0.$$

(9)

By writing the time independent Dirac spinor in the separated form

$$\psi (r,\theta )=\beta \left(\theta \right)\Phi \left(r\right),$$

(10)

where $\beta \left(\theta \right)=\left(\begin{array}{cc}{\beta }_1\left(\theta \right)& 0\\ 0& {\beta }_2\left(\theta \right)\end{array}\right)$ and $\Phi \left(r\right)=\left(\begin{array}{c}{\varphi }_1\left(r\right)\\ i{\varphi }_2\left(r\right)\end{array}\right)$, the eigenvalue equation for $J_z$ is $J_z\left(\beta \Phi \right)=j\left(\beta \Phi \right)$, which implies

$$\begin{array}{c}\hfill -i{\partial }_{\theta }{\beta }_1+{\displaystyle \frac{1}{2}}{\beta }_1=j{\beta }_1-i{\partial }_{\theta }{\beta }_2-{\displaystyle \frac{1}{2}}{\beta }_2=j{\beta }_2\end{array}$$

and, thus,

$$\beta \left(\theta \right)=\left(\begin{array}{cc}{e}^{i\left(j-{\displaystyle \frac{1}{2}}\right)\theta }& 0\\ 0& e^{i\left(j+{\displaystyle \frac{1}{2}}\right)\theta }\end{array}\right)=\left(\begin{array}{cc}{e}^{il\theta }& 0\\ 0& e^{i\left(l+1\right)\theta }\end{array}\right),l=j-{\displaystyle \frac{1}{2}}\in \mathbb{Z}.$$

(11)

It is direct to see that $\beta$ is a unitary matrix, ${\beta }^{†}\beta =⊮$.

Now, turning (10), with (11), into the Dirac equation (8), and multiplying from the left by ${\beta }^{†}$, we obtain

$$\left[\left({\beta }^{†}{\gamma }^0\beta \right)E+i\left({\beta }^{†}{\gamma }^r\beta \right){\partial }_r+{\displaystyle \frac{i}{r}}\left({\beta }^{†}{\gamma }^{\theta }\beta \right)\left({\beta }^{†}{\partial }_{\theta }\beta \right)-m\right]\Phi \left(r\right)=\left[{\beta }^{†}\mathcal{I}(r,\theta )\beta \right]\Phi \left(r\right).$$

(12)

From (11), (6) and (2) we have

$$\begin{array}{ccc}\hfill {\beta }^{†}{\gamma }^0\beta & =& {\gamma }^0,{\beta }^{†}{\gamma }^r\beta ={\gamma }^1,{\beta }^{†}{\gamma }^{\theta }\beta ={\gamma }^2,\hfill \\ \hfill {\beta }^{†}\mathcal{I}(r,\theta )\beta & =& B\left(r\right)⊮+A_0\left(r\right){\gamma }^0+A_r\left(r\right){\gamma }^1+A_{\theta }\left(r\right){\gamma }^2\equiv \tilde{\mathcal{I}}\left(r\right),\hfill \\ \hfill {\beta }^{†}\left({\partial }_{\theta }\beta \right)& =& i\left(l+{\displaystyle \frac{1}{2}}\right)⊮-{\displaystyle \frac{i}{2}}{\gamma }^0.\hfill \end{array}$$

(13)

Substituting the above results into the equation (12), we obtain the radial Dirac equation

$$\left[{\gamma }^{0}E+i{\gamma }^1\left({\partial }_r+{\displaystyle \frac{1}{2r}}\right)-{\displaystyle \frac{1}{r}}\left(l+{\displaystyle \frac{1}{2}}\right){\gamma }^2-m\right]\Phi \left(r\right)=\tilde{\mathcal{I}}\left(r\right)\Phi \left(r\right),$$

(14)

which can be written in terms of an effective potential:

$$\left[{\gamma }^{0}E+i{\gamma }^1{\partial }_r-m\right]\Phi \left(r\right)={\mathcal{I}}_{\mathrm{eff}}\left(r\right)\Phi \left(r\right),r>0,$$

(15)

where

$$\begin{array}{ccc}\hfill {\mathcal{I}}_{\mathrm{eff}}\left(r\right)& =& \tilde{\mathcal{I}}\left(r\right)-{\displaystyle \frac{i}{2r}}{\gamma }^1+{\displaystyle \frac{1}{r}}\left(l+{\displaystyle \frac{1}{2}}\right){\gamma }^2\hfill \\ & =& B\left(r\right)⊮+A_0\left(r\right){\gamma }^0+A_r\left(r\right){\gamma }^1+A_{\theta }\left(r\right){\gamma }^2-{\displaystyle \frac{i}{2r}}{\gamma }^1+{\displaystyle \frac{1}{r}}\left(l+{\displaystyle \frac{1}{2}}\right){\gamma }^2.\hfill \end{array}$$

(16)

Conservation of the Dirac Current

The Dirac current in (2+1) dimensions, in cartesian components, is given by

$$j^{\mu }\left(x\right)={\Psi }^{†}\left(x\right){\gamma }^0{\gamma }^{\mu }\Psi \left(x\right),x=\left(x^{\mu }\right),\mu =0,1,2.$$

(17)

The current conservation states that ${\partial }_{\mu }{j}^{\mu }\left(x\right)=0$, which, in the cartesian coordinates of the laboratory frame, reads as

$${\partial }_k{j}^k\left(x^{\mu }\right)=-{\partial }_0{j}^0\left(x^{\mu }\right).$$

(18)

Integrating the above on a arbitrary circle $S_{\rho }$ of radius $\rho$, and using the divergence theorem, we obtain (below, $\widehat{\mathbf{r}}$ is the unit vector in the radial direction):

$${\int }_{S_{\rho }}\nabla ·\mathbf{j}\left(x^{\mu }\right)ds={\oint }_{C_{\rho }}\mathbf{j}\left(x^{\mu }\right)·\widehat{\mathbf{r}}dl=-{\displaystyle \frac{d}{dt}}{\int }_{S_{\rho }}{j}^0\left(x^{\mu }\right)ds,2\pi \rho j^r(\rho ,t)=-{\displaystyle \frac{d}{dt}}{Q}_{\rho }\left(t\right),$$

(19)

where $Q_{\rho }\left(t\right)$ is the total probability to find the Dirac particle enclosed inside the circumference $C_{\rho }$ of radius $\rho$. Above, we identified the radial component of the Dirac current, $j^r\left(\rho ,t\right)$, and used the fact that the radial symmetry implies that it does not depend on $\theta$ (see below).

Now we restrict ourselves to the stationary case. From the radial symmetry, we have

$$j^{\mu }={\Phi }^{†}{\beta }^{†}{\gamma }^0{\gamma }^{\mu }\beta \Phi ={\Phi }^{†}\left({\beta }^{†}{\gamma }^0\beta \right)\left({\beta }^{†}{\gamma }^{\mu }\beta \right)\Phi ={\Phi }^{†}{\gamma }^0\left({\beta }^{†}{\gamma }^{\mu }\beta \right)\Phi ,\mu =0,1,2.$$

(20)

The explicit form of the component $j^0$ is

$$j^0={\Phi }^{†}\Phi =\overline{\Phi }{\gamma }^0\Phi ,\overline{\Phi }\equiv {\Phi }^{†}{\gamma }^0,$$

and, for the cartesian space components, we have

$$\begin{array}{ccc}\hfill j^1=\overline{\Phi }\left({\beta }^{†}{\gamma }^1\beta \right)\Phi & =& \overline{\Phi }\left({\gamma }^{1}cos\theta -{\gamma }^{2}sin\theta \right)\Phi =j^{r}cos\theta -j^{\theta }sin\theta ,\hfill \\ \hfill j^2=\overline{\Phi }\left({\beta }^{†}{\gamma }^2\beta \right)\Phi & =& \overline{\Phi }\left({\gamma }^{1}sin\theta +{\gamma }^{2}cos\theta \right)\Phi =j^{r}sin\theta +j^{\theta }cos\theta ,\hfill \end{array}$$

where, from (6) and (20), we have identified the r and $\theta$-components of the stationary Dirac current as

$$\begin{array}{ccc}\hfill j^r& =& j^{1}cos\theta +j^{2}sin\theta =\overline{\psi }{\gamma }^r\psi =\overline{\Phi }{\gamma }^1\Phi ,\hfill \end{array}$$

(21)

$$\begin{array}{ccc}\hfill j^{\theta }& =& -j^{1}sin\theta +j^{2}cos\theta =\overline{\psi }{\gamma }^{\theta }\psi =\overline{\Phi }{\gamma }^2\Phi .\hfill \end{array}$$

(22)

From (21) and (22) it is clear that both the current’s polar components $j^r$ and $j^{\theta }$ do not depend on the angular coordinate $\theta$, as mentioned earlier. Thus, we write $j^r(r,\theta )\equiv j\left(r\right)$ and $j^{\theta }(r,\theta )\equiv j^{\theta }\left(r\right)$. From the above we have that along the angular direction both the stationary currents $j^r$ and $j^{\theta }$ are constant and, from (19), in the stationary case we have that the quantity $r{j}^r$ vanishes, which implies that $j^r\left(r\right)=0$, $\forall r>0$. In particular, for any radius $r=\rho >0$:

$$j^r\left({\rho }^{-}\right)=j^r\left({\rho }^{+}\right)=j^r\left(\rho \right)=0,\mathrm{where}{j}^r\left({\rho }^{\pm }\right)\equiv \underset{r\to {\rho }^{\pm }}{lim}{j}^r\left(r\right).$$

(23)

Summarizing the results obtained so far, we have that, for an external, regular, and radial symmetric potential, $\mathcal{I}(r,\theta )=\left[B\left(r\right)⊮+A_0\left(r\right){\gamma }^0+A_r\left(r\right){\gamma }^r+A_{\theta }\left(r\right){\gamma }^{\theta }\right]$, the radial spinor $\Phi \left(r\right)$ must solve the one-dimensional stationary Dirac equation (15) with the effective one-dimensional potential (16). Conservation of the Dirac current in the original two-dimensional problem implies that the one-dimensional (radial) stationary current $j^r\left(r\right)$, given by (21), is continuous (and vanishes) at any r, and therefore satisfies condition (23).

Equation (15) is mathematically analogous to the Dirac equation in (1+1) dimensions, with an external potential ${\mathcal{I}}_{\mathrm{eff}}\left(r\right)$ (see Equation (3) at p. 4 of [31]). However, (15) holds only for $r>0$. In the next Section we will see how to extend it as a distributional equation on the interval $r\in [0,\infty )$, and which condition the Dirac spinor $\psi$ must satisfy at $r=0$. The analogous roles of the potentials $A_1\left(x\right)$ and $W\left(x\right)$ of the regular, true one-dimensional case (see the Equation (2) of [31]) will be played in (15) by the following combinations

$$\begin{array}{c}\hfill A_1\left(x^1\right)⟶A_r\left(r\right)-{\displaystyle \frac{i}{2r}},W\left(x^1\right)⟶-A_{\theta }\left(r\right)-{\displaystyle \frac{1}{r}}\left(l+{\displaystyle \frac{1}{2}}\right).\end{array}$$

(24)

It is worth to observe that in (15) the effective one-dimensional potential in the radial direction, ${\gamma }^0{\mathcal{I}}_{\mathrm{eff}}\left(r\right)$, is not Hermitian, due to the presence of the imaginary quantity $-{\displaystyle \frac{i}{2r}}$ added to the potential $A_r\left(r\right)$ in (24). However, this is immaterial, since for the distributional method of [31] to be applicable in the present case it suffices that the current be continuous at $r=R$. Finally, it is worth to note that the original potential ${\gamma }^0\mathcal{I}(r,\theta )$ is in fact Hermitian, and the Lorentz vector $j^{\mu }$ is conserved.

3. Dirac Equation with Radial Contact Potentials at $\mathbf{r}=\mathbf{R}$

Now we consider the Dirac equation with the most general, time-independent, contact potential supported on a circumference of radius R in the laboratory frame. Since now the potential is singular, the products between the potentials and the spinors $\Phi \left(r\right)$ appearing in (15) are not well defined. Guided by the analogy to the one-dimensional case, we follow closely the steps of reference [31] and write the time independent Dirac equation as

$$\left[{\gamma }^{0}E+i{\gamma }^1{\partial }_r-m\right]\Phi \left(r\right)=D[\Phi ]\left(r\right)+{\mathcal{I}}_{\mathrm{eff}}^{\mathrm{reg}}\left(r\right)\Phi \left(r\right),r>0.$$

(25)

where

$${\mathcal{I}}_{\mathrm{eff}}^{\mathrm{reg}}\left(r\right)=-{\displaystyle \frac{i}{2r}}{\gamma }^1+{\displaystyle \frac{1}{r}}\left(l+{\displaystyle \frac{1}{2}}\right){\gamma }^2$$

(26)

is the regular part of the effective potential. Since this term defines a regular distribution on $[0,\infty )$, which includes the origin, the equation (25), understood as a distributional equation, can be extended to include the origin, which will be assumed from now on.4 The first term in the right hand side of (25) is the contact term, and is a $2\times 1$ matrix, yet unknown, whose elements are singular (not regular) distributions, related to the singular potentials concentrated at $r=R$. The distribution space we consider here is the Schwartz space ${\mathcal{S}}^{\prime }$, which consists of continuous linear functionals acting on the Schwartz space $\mathcal{S}$ of infinitely differentiable strongly decreasing test functions. For details, we refer the reader to reference [31]. The distribution $D[\Phi ]\left(x\right)$ will be determined from the following three basic requirements (adapted from [31]):

(R1) 

$D[\Phi ]\left(r\right)$ must be a distribution concentrated at $r=R$, i.e., it must equals the zero distribution in ${\mathbb{R}}_{+}\setminus \left\{R\right\}$.

(R2) 

The two components of the radial Dirac spinor $\Phi \left(r\right)$ must correspond to regular distributions on ${\mathbb{R}}_{+}=[0,\infty )$, i.e., they must be slow growth functions and locally integrable (in the Lebesgue sense) on ${\mathbb{R}}_{+}$. This is equivalent to saying that $\Phi \left(r\right)$ must have order $s_{\Phi }\le -1$ on any closed interval $K\ni R$.

(R3) 

The Dirac current must be conserved across $r=R$. In the stationary case this implies that $j^r\left(R^{-}\right)=j^r\left(R^{+}\right)=0$.

The definition we are using above for the order of a distribution on a closed, finite interval $K\subset {\mathbb{R}}_{+}$, is the same of [31], adapted to the semi-infinite axis ${\mathbb{R}}_{+}$. Briefly, with $R\in K$, the order of the Dirac delta distribution $\delta (r-R)$ on K is $s=0$. Taking the derivative of any distribution increases its order on K by $+1$, whereas taking an indefinite integral (a primitive) decreases its order by $-1$. The order of a sum of distributions on K is the largest order of its terms. A distribution has order $s=-\infty$ on K if it is infinitely differentiable on K. A distribution has order $s=-2$ on K if on this interval it coincides with a continuous but not differentiable (in the ordinary sense) function. A distribution of order $s=-1$ corresponds to a discontinuous function, and may or may not be regular. A distribution having order $s\ge 0$ on K is always singular. For more details, we refer the reader to [31].

From requirement (R2), and balancing the orders on both sides of (25), the order of the interaction distribution $D[\Phi ]$ on any interval $K\ni R$ must be $s_{D[\Phi ]}\le 0$. This same requirement also sets a maximum for the order of the spinor $\Phi \left(r\right)$ on a closed finite interval $I_0$ containing the origin, but not containing R, i.e., $I_0=[0,b]$, $b<R$. Since on $I_0$ the contact distribution $D[\Phi ]$ equals the zero distribution, and the term (26) has order $\le -1$ (since it is regular on $I_0$, see footnote 4), by balancing the orders on both sides of (25) we conclude that $s_{{\Phi }^{\prime }}\le -1$ on $I_0$ and, consequently, $s_{\Phi }\le -2$ on $I_0$. This means that the spinor $\Phi$ must be continuous on $I_0$ and, therefore, it must be bounded at the origin.5 This is another way to justify, in a mathematically rigorous way, the usual assumption of boundedness of the Dirac spinor at $r=0$.

Theorem 3.5-2 of [32] states that a distribution that is concentrated at a single point $r=R$, and that has singular order s on a closed finite interval $K\ni R$, must in this interval coincide with a finite linear combination of Dirac deltas and its derivatives, up to order s. Since on K the order of $\delta (r-R)$ is zero, we must have $s_{D[\Phi ]}\ge 0$. We have that, on K, $s_{\Phi }\le -1$ (from (R2)) and, consequently, $s_{{\Phi }^{\prime }}\le 0$. Now, by balancing the orders on both sides of (25), we have that $s_{D[\Phi ]}=s_{{\Phi }^{\prime }}\le 0$, and therefore $s_{D[\Phi ]}=0$ on K. Requirements (R1) and (R2) thus imply that

$$D[\Phi ]\left(r\right)=\Omega [\Phi ]\delta (r-R),$$

(27)

where the $2\times 1$ matrix $\Omega [\Phi ]$ does not depend on the coordinate r, and will have a functional (linear) dependence of the spinor $\Phi$. The form of $\Omega [\Phi ]$ will completely specify the contact term, and below it will be found from requirement (R3). As will be shown, the obtained form of $\Omega [\Phi ]$ will imply matching conditions that $\Phi \left(r\right)$ must satisfy at $r=R$, similar to the true one-dimensional case. Also, it will provide the contact term $\Omega [\Phi ]\delta (r-R)$ explicitly as a well defined distribution defined on the semi-axis ${\mathbb{R}}_{+}$.

To determine $\Omega [\Phi ]$ we proceed exactly in the same way as in [31]. By taking a distributional indefinite integral of equation (25) on $K\subset {\mathbb{R}}_{+}$, $K\ni R$, we obtain

$$\left[{\gamma }^{0}E-m\right]{\Phi }^{(-1)}\left(r\right)+i{\gamma }^1\Phi \left(r\right)=\Omega [\Phi ]\theta (r-R)+F\left(r\right)+c,$$

(28)

where ${\Phi }^{(-1)}$ is any primitive of $\Phi$, $F\left(r\right)$ is a primitive of the product ${\mathcal{I}}_{\mathrm{eff}}^{\mathrm{reg}}\left(r\right)\Phi \left(r\right)\in {\mathcal{L}}_{\mathrm{loc}}^1\left(K\right)$ and c is a constant and arbitrary column matrix. Similarly to the reasoning in [31], we conclude that both functions, ${\Phi }^{(-1)}$ and $F\left(r\right)$, are continuous on K.6 By isolating the term $i{\gamma }^1\Phi \left(r\right)$ in the above equation, we find that $\Phi \left(r\right)$ must have well defined lateral limits at $r=R$. By considering the limits $r\to R^{\pm }$ in (28) and taking the difference, we obtain:

$$i{\gamma }^1\left[\Phi \left(R^{+}\right)-\Phi \left(R^{-}\right)\right]=\Omega [\Phi ].$$

(29)

Now the functional coefficient $\Omega [\Phi ]$ will be determined from requirement (R3):

$$j^r\left(R^{-}\right)=j^r\left(R^{+}\right),\mathrm{with}{j}^r\left(r\right)=\overline{\Phi }{\gamma }^1\Phi .$$

(30)

The problem of determining $\Omega [\Phi ]$ from (29), by using the current condition (30), is mathematically analogous to the one-dimensional problem treated in [31]. Therefore, we just collect the results from [31], replacing the analogous quantities to the present problem:

$$\Omega [\Phi ]=\left(B⊮+A_0{\gamma }^0+A_r{\gamma }^1+A_{\theta }{\gamma }^2\right){\displaystyle \frac{\Phi \left(R^{-}\right)+\Phi \left(R^{+}\right)}{2}},$$

(31)

where the real constants $B,A_0,A_r$ and $A_{\theta }$ are identified with the strengths of the singular potentials associated to the contact term supported on $r=R$, namely, from (27)

$$B\left(r\right)=B\delta (r-R),A_0\left(r\right)=A_0\delta (r-R),A_r\left(r\right)=A_r\delta (r-R)\mathrm{and}{A}_{\theta }\left(r\right)=A_{\theta }\delta (r-R).$$

(32)

The radial Dirac equation (25), with the contact and the regular interaction terms given explicitly, is:

$$\begin{array}{ccc}\hfill \left({\gamma }^{0}E+i{\gamma }^1{\partial }_r-m\right)\Phi \left(r\right)& =& \Omega [\Phi ]\delta (r-R)+{\mathcal{I}}_{\mathrm{eff}}^{\mathrm{reg}}\left(r\right)\Phi \left(r\right)\hfill \\ & =& \left(B⊮+A_0{\gamma }^0+A_r{\gamma }^1+A_{\theta }{\gamma }^2\right)\delta (r-R){\displaystyle \frac{\Phi \left(R^{-}\right)+\Phi \left(R^{+}\right)}{2}}\hfill \\ & & +\left[-{\displaystyle \frac{i}{2r}}{\gamma }^1+{\displaystyle \frac{1}{r}}\left(l+{\displaystyle \frac{1}{2}}\right){\gamma }^2\right]\Phi \left(r\right).\hfill \end{array}$$

(33)

A permeable contact interaction is characterised by the condition (see Equation (55) of [31])

$$A_r\ne 0\mathrm{or}{B}^2-4-A_0^2+A_{\theta }^2\ne 0,$$

(34)

which means that the solutions at $r=R^{-}$ and $r=R^{+}$ are connected by the following matching conditions

$$\Phi \left(R^{+}\right)=\Lambda \Phi \left(R^{-}\right),\mathrm{with}\Lambda =e^{i\phi }\left(\begin{array}{cc}a& ib\\ -ic& d\end{array}\right),ad-bc=1,\phi \in [0,\pi ),$$

(35)

where the $\Lambda$ matrix is written in terms of the physical strengths as7

$$\Lambda =\frac{e^{iarg\left(B^2-4-A_0^2+A_r^2+A_{\theta }^2+4i{A}_r\right)}}{\sqrt{{\left(B^2-4-A_0^2+A_r^2+A_{\theta }^2\right)}^2+16A_r^2}}\left(\begin{array}{cc}-B^2+A_0^2-A_r^2-{(A_{\theta }+2)}^2& 4i(A_0-B)\\ 4i(A_0+B)& -B^2+A_0^2-A_r^2-{(A_{\theta }-2)}^2\end{array}\right).$$

(36)

Yet, to have a unique set of $\Lambda$-parameters $\phi ,a,b,c,d$, given a set of parameters $B,A_0,$ and $A_r$ and $A_{\theta }$ satisfying the permeability conditions (34), we can write the relationships:

$$\begin{array}{ccc}\hfill \phi & =& {tan}^{-1}\left(\frac{4A_r}{B^2-4-A_0^2+A_r^2+A_{\theta }^2}\right),\phi \in [0,\pi ),\hfill \end{array}$$

(37)

$$\begin{array}{ccc}\hfill a& =& \pm \frac{A_0^2-A_r^2-B^2-{(A_{\theta }+2)}^2}{\sqrt{{\left(B^2-4-A_0^2+A_r^2+A_{\theta }^2\right)}^2+16A_r^2}},b=\pm \frac{4(A_0-B)}{\sqrt{{\left(B^2-4-A_0^2+A_r^2+A_{\theta }^2\right)}^2+16A_r^2}},\hfill \end{array}$$

(38)

$$\begin{array}{ccc}\hfill c& =& \pm \frac{-4(A_0+B)}{\sqrt{{\left(B^2-4-A_0^2+A_r^2+A_{\theta }^2\right)}^2+16A_r^2}},d=\pm \frac{A_0^2-A_r^2-B^2-{(A_{\theta }-2)}^2}{\sqrt{{\left(B^2-4-A_0^2+A_r^2+A_{\theta }^2\right)}^2+16A_r^2}},\hfill \end{array}$$

(39)

where in the above expressions the plus (minus) sign must be taken if $A_r>0$ ($A_r<0$); if $A_r=0$ we must take the same sign of $B^2-4-A_0^2+A_{\theta }^2$ (see [31]).

When the strength potentials are such that the permeability condition (34) is not met, the circular wall at $r=R$ is impenetrable, and some (or all) of the $\Lambda$ matrix parameters becomes infinite. In this case, the inner and outer solutions are independent, the (anti)particle is completely confined in the inner (outer) region of the circumference of radius R, and it cannot be transmitted to the outer (inner) region. The set of allowed energies for the (anti)particle confined in the inner region will depend on the boundary conditions at $r=R$, which in turn depend on the values of the strength potentials. In the outer region all energies $\left|E\right|>m$ are allowed; the $E=\pm m$ cases may or may not be allowed in the outer region, as we will see later. In the following sections, we will use the matching conditions (35), together with the form (36) of the $\Lambda$ matrix in terms of the physical strengths, to study bound, scattering, and resonant states for various choices of the contact potentials. Cases that give rise to an impermeable wall can be obtained by considering the corresponding limits on the $\Lambda$ parameters.

4. Energy Eigenstates

In this section we will determine the energy eigenstates for the radial Dirac equation (15), for a given angular momentum $l>0$. Consider the spinorial form of the solutions, for the inner ($r<R$, $\alpha =i$) or the outer ($r>R$, $\alpha =o$) regions

$${\Phi }_{\alpha }\left(r\right)=\left(\begin{array}{c}{\varphi }_{1,\alpha }\left(r\right)\\ {\displaystyle i{\varphi }_{2,\alpha }\left(r\right)}\end{array}\right).$$

(40)

The spinors ${\Phi }_{\alpha }\left(r\right)$ are solutions of equation (14), with $\tilde{\mathcal{I}}\left(r\right)=0$ (the free equation). In terms of the spinor components ${\varphi }_{k,\alpha }\left(r\right)$, $k=1,2$, $\alpha =i,o$, that equation reads

$$-{\displaystyle \frac{d{\varphi }_{1,\alpha }}{dr}}+{\displaystyle \frac{l}{r}}{\varphi }_{1,\alpha }=(E+m){\varphi }_{2,\alpha },{\displaystyle \frac{d{\varphi }_{2,\alpha }}{dr}}+{\displaystyle \frac{l+1}{r}}{\varphi }_{2,\alpha }=(E-m){\varphi }_{1,\alpha },\forall r\ne R.$$

(41)

When $\left|E\right|\ne m$, by isolating ${\varphi }_{2,\alpha }$ from the first of equations (41) and substituting it into the second equation, we obtain a Schrödinger-like equation for ${\varphi }_{1,\alpha }$, valid $\forall r\ne R$,

$$\begin{array}{ccc}\hfill -{\displaystyle \frac{d^2{\varphi }_{1,\alpha }}{d{r}^2}}-{\displaystyle \frac{1}{r}}{\displaystyle \frac{d{\varphi }_{1,\alpha }}{dr}}+{\displaystyle \frac{l^2}{r^2}}{\varphi }_{1,\alpha }& =& \left(E^2-m^2\right){\varphi }_{1,\alpha },\hfill \\ \hfill r^2{\displaystyle \frac{d^2{\varphi }_{1,\alpha }}{d{r}^2}}+r{\displaystyle \frac{d{\varphi }_{1,\alpha }}{dr}}+\left(p^2{r}^2-l^2\right){\varphi }_{1,\alpha }& =& 0,\hfill \end{array}$$

(42)

where we defined $p=\sqrt{E^2-m^2}$. The solution of (42) is (see equation (2.13) of [33])8

$$\begin{array}{ccc}\hfill {\varphi }_{1,\alpha }\left(r\right)& =& a_{\alpha }{J}_l\left(pr\right)+b_{\alpha }{Y}_l\left(pr\right)={\tilde{a}}_{\alpha }{H}_l^{\left(1\right)}\left(pr\right)+{\tilde{b}}_{\alpha }{H}_l^{\left(2\right)}\left(pr\right),\alpha =i,o,\hfill \end{array}$$

(43)

where $J_l$($H_l^{\left(1\right)}$) and $Y_l$($H_l^{\left(2\right)}$) are, respectively, the Bessel (Hankel) functions of the first and of the second kind. Substituting this solution into the first of equations (41), we obtain ${\varphi }_{2,\alpha }\left(r\right)$ (see also [33]):

$${\varphi }_{2,\alpha }\left(r\right)={\displaystyle \frac{p}{E+m}}\left[a_{\alpha }{J}_{l+1}\left(pr\right)+b_{\alpha }{Y}_{l+1}\left(pr\right)\right]={\displaystyle \frac{p}{E+m}}\left[{\tilde{a}}_{\alpha }{H}_{l+1}^{\left(1\right)}\left(pr\right)+{\tilde{b}}_{\alpha }{H}_{l+1}^{\left(2\right)}\left(pr\right)\right],\alpha =i,o.$$

(44)

Thus, the radial part of the energy and angular momentum eingenstates (40) of the free Dirac particle in the inner and outer regions can be written in the form

$$\begin{array}{ccc}\hfill {\Phi }_{\alpha }& =& a_{\alpha }\left(\begin{array}{c}{J}_l\left(pr\right)\\ {\displaystyle \frac{ip}{E+m}}{J}_{l+1}\left(pr\right)\end{array}\right)+b_{\alpha }\left(\begin{array}{c}{Y}_l\left(pr\right)\\ {\displaystyle \frac{ip}{E+m}}{Y}_{l+1}\left(pr\right)\end{array}\right)\hfill \\ & =& {\tilde{a}}_{\alpha }\left(\begin{array}{c}{H}_l^{\left(1\right)}\left(pr\right)\\ {\displaystyle \frac{ip}{E+m}}{H}_{l+1}^{\left(1\right)}\left(pr\right)\end{array}\right)+{\tilde{b}}_{\alpha }\left(\begin{array}{c}{H}_l^{\left(2\right)}\left(pr\right)\\ {\displaystyle \frac{ip}{E+m}}{H}_{l+1}^{\left(2\right)}\left(pr\right)\end{array}\right),\alpha =i,o.\hfill \end{array}$$

(45)

Since the inner solution ${\Phi }_i\left(r\right)$ must be bounded at the origin (see the discussion regarding requirement (R2) in Section 4), and this holds only for the spinor whose elements are Bessel functions of the first kind, from now on we will consider only the first line of equations (45), with $b_i=0$, as the solution for the inner region.

4.1. Critical and Supercritical States ($\left|E\right|=m$)

The states for which $E=m$ or $E=-m$ are called, respectively, critical and supercritical states and are important since, as we will see later, by varying the strengths of the potential a resonance can transform into a bound state, and vice versa. In the former case we say that a bound state is captured (or absorbed) from the continuum, and in the latter, that a bound state was emitted into the continuum. Captures or emissions of bound states from and to the continuum always occur at a critical or supercritical energy. To solve the Dirac equation for these states we must consider the equation (41) with $E=\pm m$.

4.1.1. Critical States

The critical states ($E=m$) solve the equation (42), with $p=0$ (an Euler equation),

$$r^2{\displaystyle \frac{d^2{\varphi }_{1,\alpha }}{d{r}^2}}+r{\displaystyle \frac{d{\varphi }_{1,\alpha }}{dr}}-l^2{\varphi }_{1,\alpha }=0,\forall r\ne R,$$

(46)

and the first of equations (41), with $E=m$:

$$-{\displaystyle \frac{d{\varphi }_{1,\alpha }}{dr}}+{\displaystyle \frac{l}{r}}{\varphi }_{1,\alpha }=2m{\varphi }_{2,\alpha }.$$

(47)

The solutions of (46) have the form $r^{\pm l}$, but since we only consider $l\ge 0$ and require the Dirac spinor be bounded at $r=0$, only the solution $r^l$ is admissible in the inner region. We also require the Dirac spinor to be bounded at infinity (a physical requirement), so the only admissible solution for the outer region is $r^{-l}$. Therefore, the solutions for (46) in the inner and outer regions are

$${\varphi }_{1,i}=c_i{r}^l\mathrm{and}{\varphi }_{1,o}=c_o{r}^{-l}$$

(48)

with $c_i$ and $c_o$ being arbitrary constants. substituting these solutions into (47), we obtain ${\varphi }_{2,i}=0$ and ${\varphi }_{2,o}=c_o{\displaystyle \frac{l}{m}}{r}^{-(l+1)}$. Therefore, the solutions for critical states in the inner and outer regions are

$${\Phi }_i^{\mathrm{crit}}=c_i\left(\begin{array}{c}{r}^l\\ 0\end{array}\right)\mathrm{and}{\Phi }_o^{\mathrm{crit}}=c_o\left(\begin{array}{c}{r}^{-l}\\ {\displaystyle \frac{il}{m}}{r}^{-(l+1)}\end{array}\right).$$

(49)

These two solutions must match at $r=R$ using the matrix $\Lambda$, according to (35),

$${\Phi }_o^{\mathrm{crit}}\left(R\right)=\Lambda {\Phi }_i^{\mathrm{crit}}\left(R\right),$$

(50)

that results in

$$c_o\left(\begin{array}{c}{R}^{-l}\\ {\displaystyle \frac{il}{m}}{R}^{-(l+1)}\end{array}\right)=c_i\Lambda \left(\begin{array}{c}{R}^l\\ 0\end{array}\right).$$

(51)

This equation has non-trivial solutions for $c_{i,o}$ if, and only if, the following secular equation is satisfied by the parameters of the matrix $\Lambda$

$$c+a{\displaystyle \frac{l}{mR}}=0.$$

(52)

From (49) we observe that critical states with $l=0,1$ are not true bound states: they are quasi-bound states, since they are bounded or vanish at infinity, but are not square integrable.9 For $l\ge 2$ the critical states are true bound states.

4.1.2. Supercritical States

To obtain the supercritical states ($E=-m$) we isolate ${\varphi }_{1,\alpha }$ in the second of the equations (41) and substitute it in the first, obtaining another Euler equation

$$r^2{\displaystyle \frac{d^2}{d{r}^2}}{\varphi }_{2,\alpha }+r{\displaystyle \frac{d}{dr}}{\varphi }_{2,\alpha }-{\left(l+1\right)}^2{\varphi }_{2,\alpha }=0,$$

(53)

whose solutions are of the form $r^{\pm (l+1)}$. Similarly to the previous case, for the inner region the solution with the positive sign is the only admissible one, while for the outer region we need only consider the solution with the negative sign. The corresponding solutions for ${\varphi }_{1,\alpha }$ are obtained from the second of the equations (41). Therefore, the admissible supercritical solutions in the inner and outer regions are

$${\Phi }_i^{\mathrm{super}}=c_i\left(\begin{array}{c}-{\displaystyle \frac{(l+1)}{m}}{r}^l\\ i{r}^{l+1}\end{array}\right)\mathrm{and}{\Phi }_o^{\mathrm{super}}=c_o\left(\begin{array}{c}0\\ i{r}^{-(l+1)}\end{array}\right).$$

(54)

Again, to match solutions at $r=R$ we use (35), ${\Phi }_o^{\mathrm{super}}\left(R\right)=\Lambda {\Phi }_i^{\mathrm{super}}\left(R\right)$, and we obtain the following secular equation for the existence of non-trivial solutions:

$$b+a{\displaystyle \frac{(l+1)}{mR}}=0.$$

(55)

From (54) we observe that, with the exception of $l=0$, the supercritical states correspond to true bound states. Let us emphasize that the phase factor $e^{i\phi }$ of the matrix $\Lambda$ does not enter the conditions (52) and (55) for the existence of critical or supercritical states.

Using the equations (38)-(39) we can write the conditions for the existence of critical and supercritical states, (52) and (55), in terms of the physical parameters B and $A_{\mu }$ ($\mu =0,1,2$).

4.2. Bound States ($\left|E\right|<m$)

The bound states are characterized by $\left|E\right|<m$, and in this case $p=iq$, with $q=\sqrt{m^2-E^2}>0$. For the outer region, it is more convenient to consider the solution in terms of Hankel functions, since it is easier to explore their asymptotic behavior as $r\to \infty$. The Hankel functions behave asymptotically as

$$H_l^{\left(1\right)}\left(iqr\right)\approx \sqrt{{\displaystyle \frac{2}{i\pi qr}}}{e}^{-qr-{\displaystyle \frac{i}{4}}(2l+1)\pi },H_l^{\left(2\right)}\left(iqr\right)\approx \sqrt{{\displaystyle \frac{2}{i\pi qr}}}{e}^{qr-{\displaystyle \frac{i}{4}}(2l-3)\pi }.$$

(56)

From the above we see that the Hankel function $H_l^{\left(2\right)}\left(iqr\right)$ explodes as r goes to infinity, while $H_l^{\left(1\right)}\left(iqr\right)$ gives an evanescent wave. So, for bound states, we must have ${\tilde{b}}_o=0$ in the second line of equations (45). The inner and outer solutions must be matched at $r=R$ by using (35) :

$$\begin{array}{c}\hfill {\tilde{a}}_o\left(\begin{array}{c}{H}_l^{\left(1\right)}\left(iqR\right)\\ {\displaystyle \frac{-q}{E+m}}{H}_{l+1}^{\left(1\right)}\left(iqR\right)\end{array}\right)=a_i\Lambda \left(\begin{array}{c}{J}_l\left(iqR\right)\\ {\displaystyle \frac{-q}{E+m}}{J}_{l+1}\left(iqR\right)\end{array}\right),\end{array}$$

(57)

that gives

$$\mathcal{M}\left(iqR\right)\left(\begin{array}{c}{\tilde{a}}_o\\ a_i\end{array}\right)=\left(\begin{array}{c}0\\ 0\end{array}\right),$$

(58)

where $\mathcal{M}\left(iqR\right)=\left({\mathcal{C}}_o\left(iqR\right){\mathcal{C}}_i\left(iqR\right)\right)$ is the matrix whose columns are

$${\mathcal{C}}_o\left(iqR\right)=\left(\begin{array}{c}{H}_l^{\left(1\right)}\left(iqR\right)\\ {\displaystyle \frac{-q}{E+m}}{H}_{l+1}^{\left(1\right)}\left(iqR\right)\end{array}\right)\mathrm{and}{\mathcal{C}}_i\left(iqR\right)=-\Lambda \left(\begin{array}{c}{J}_l\left(iqR\right)\\ {\displaystyle \frac{-q}{E+m}}{J}_{l+1}\left(iqR\right)\end{array}\right).$$

(59)

The condition for (58) to have non-trivial solutions for the coefficients ${\tilde{a}}_o$ and $a_i$ is that $det\mathcal{M}\left(iqR\right)=0$. By writing $det\mathcal{M}={\mathcal{C}}_o^T{\gamma }^1{\mathcal{C}}_i$, where the superscript T denotes the matrix transpose, we obtain the secular equation for the bound state energies:

$$\left(\begin{array}{cccc}{H}_l^{\left(1\right)}\left(iqR\right)& & & {\displaystyle \frac{-q}{E+m}}{H}_{l+1}^{\left(1\right)}\left(iqR\right)\end{array}\right){\gamma }^1\Lambda \left(\begin{array}{c}{J}_l\left(iqR\right)\\ {\displaystyle \frac{-q}{E+m}}{J}_{l+1}\left(iqR\right)\end{array}\right)=0.$$

(60)

From the above equation we observe that the phase factor $e^{i\phi }$ of the $\Lambda$ matrix again does not affect the energies of the eventual bound state. In terms of the $\Lambda$ matrix parameters, the above equation may be written as as

$$\begin{array}{ccc}& & I_l\left(qR\right)\left[aq{K}_{l+1}\left(qR\right)+c(m+E)K_l\left(qR\right)\right]\hfill \\ & & +I_{l+1}\left(qR\right)\left[b(m-E)K_{l+1}\left(qR\right)+dq{K}_l\left(qR\right)\right]=0,\hfill \end{array}$$

(61)

where $I_l\left(z\right)$ and $K_l\left(z\right)$ represent the modified Bessel functions of the first and second kinds, respectively. Again, the above secular equation can be written in terms of the physical parameters B and $A_{\mu }$ ($\mu =0,1,2$) using (38)-(39).

4.3. Scattering States ($\left|E\right|>m$)

The scattering states correspond to energies $\left|E\right|>m$, such that $p=\sqrt{E^2-m^2}>0$. The outer solution, written in terms of the Bessel functions of the first and second kinds, is given by (45):

$${\Phi }_o=a_o\left(\begin{array}{c}{J}_l\left(pr\right)\\ {\displaystyle \frac{ip}{E+m}}{J}_{l+1}\left(pr\right)\end{array}\right)+b_o\left(\begin{array}{c}{Y}_l\left(pr\right)\\ {\displaystyle \frac{ip}{E+m}}{Y}_{l+1}\left(pr\right)\end{array}\right).$$

(62)

The asymptotic behavior of these Bessel functions, that is,

$$J_l\left(pr\right)\approx \sqrt{{\displaystyle \frac{2}{\pi pr}}}cos\left(pr-{\displaystyle \frac{1}{4}}(2l+1)\pi \right),Y_l\left(pr\right)\approx \sqrt{{\displaystyle \frac{2}{\pi pr}}}sin\left(pr-{\displaystyle \frac{1}{4}}(2l+1)\pi \right),$$

suggests the following way of writing the coefficients $a_o,b_o$, in terms of the phase shift ${\delta }_l$:

$$a_o=hcos{\delta }_l,b_o=-hsin{\delta }_l,$$

(63)

where h may be complex, but ${\delta }_l$ is real. Therefore,

$${\Phi }_o=a_o\left(\begin{array}{c}{J}_l\left(pr\right)-\left(tan{\delta }_l\right)Y_l\left(pr\right)\\ {\displaystyle \frac{ip}{E+m}}\left[J_{l+1}\left(pr\right)-\left(tan{\delta }_l\right)Y_{l+1}\left(pr\right)\right]\end{array}\right),$$

(64)

and the matching condition (35) at $r=R$ gives

$$\left[\left(\begin{array}{c}{J}_l\left(pR\right)-\left(tan{\delta }_l\right)Y_l\left(pR\right)\\ {\displaystyle \frac{ip}{E+m}}\left[J_{l+1}\left(pR\right)-\left(tan{\delta }_l\right)Y_{l+1}\left(pR\right)\right]\end{array}\right)-\Lambda \left(\begin{array}{c}{J}_l\left(pR\right)\\ {\displaystyle \frac{ip}{E+m}}{J}_{l+1}\left(pR\right)\end{array}\right)\right]\left(\begin{array}{c}{a}_o\\ a_i\end{array}\right)=\left(\begin{array}{c}0\\ 0\end{array}\right),$$

(65)

where we have represented the columns of the left-hand side matrix as column matrices. To obtain non-trivial solutions for $a_o,a_i$, the determinant of the matrix on the left-hand side must be zero. Expressed in terms of the $\Lambda$ matrix parameters, this requirement results in the following expression for the phase shift ${\delta }_l$:

$$tan{\delta }_l={\displaystyle \frac{p(a-d)J_{l+1}\left(pR\right)J_l\left(pR\right)-b(E-m)J_{l+1}{\left(pR\right)}^2+c(m+E)J_l{\left(pR\right)}^2}{J_l\left(pR\right)\left[ap{Y}_{l+1}\left(pR\right)+c(m+E)Y_l\left(pR\right)\right]-J_{l+1}\left(pR\right)\left[b(E-m)Y_{l+1}\left(pR\right)+dp{Y}_l\left(pR\right)\right]}}.$$

(66)

From the above expression we can calculate the Wigner time delay [34]

$${\tau }_l\left(E\right)=2{\displaystyle \frac{\partial {\delta }_l}{\partial E}}=\left({\displaystyle \frac{2}{1+{tan}^2{\delta }_l}}\right){\displaystyle \frac{\partial }{\partial E}}tan{\delta }_l,$$

(67)

which gives us a scale for the time the particle spends in the inner region $r<R$. The above equations for the phase shift ${\delta }_l$ and the Wigner time delay ${\tau }_l$ do not depend on the phase factor $e^{i\phi }$ of the matrix $\Lambda$ and again these expressions can be expressed in terms of the physical parameters using (38)-(39).

4.4. Resonances

A very practical method for investigating resonances is to look for the energies at which the scattering states become purely outgoing states. Of course, it is not physically possible to have only purely outgoing scattering states, and for this we need to search for complex energies. The real part of these complex energies will correspond to the resonances. Although not physically possible, purely outgoing scattering states contain important information about the system: in addition to giving the resonant energies, the inverse of the imaginary parts of these complex energies provides scales for the corresponding resonant time delays, as we will see later [33].

A second, more physical way to characterize the resonant energies of the scattering states is to associate them with the peaks of the Wigner time delay ${\tau }_l\left(E\right)$ at (67). This is quite intuitive, since resonant states are associated with larger amounts of time that the particle spends within the inner region $r<R$. We will later show, in the context of the special cases investigated, that these two criteria give in general the same values for the resonances.

To identify the outcoming and incoming scattering states, it is convenient to write the scattering solution in the external region in terms of the Hankel functions (the second line of equations (45)). We can use the asymptotic behavior of the Hankel functions (56), with $q=-ip$, to identify the outcoming and incoming scattering solutions as $H_l^{\left(1\right)}$ and $H_l^{\left(2\right)}$, respectively. Thus, the purely outcoming solutions will have only the Hankel function of the first kind $H_l^{\left(1\right)}$, and thus we put ${\tilde{b}}_o=0$ in (45). The secular equation will then be exactly the same as we obtained previously for the bound states, i.e., equation (61), but now with $q=-ip=-i\sqrt{{\left(E_R+i{E}_I\right)}^2-m^2}$, and we look for solutions of that equation for the complex energy $E=E_R+i{E}_I$. Obviously, all solutions for a real energy in the interval $\left|E\right|\le m$ will be included in the set of complex solutions, and these will correspond to the bound (or quasi-bound) states.

As we continuously change the potential strengths $B,A_0,A_r$ and $A_{\theta }$ (or, equivalently, the parameters of the matrix $\Lambda$), new bound states can be created (captured from the continuum) or destroyed (emitted into the continuum). A capture(emission) generally occurs when the energy of a previous resonant(bound) state continuously crosses the (super)critical values $E=\pm m$, becoming a new bound(resonant) state. The values of the potential strengths (or $\Lambda$ parameters) for which a capture or emission from or to a bound state occurs are given by the conditions (52) and (55) for critical ($E=m$) or supercritical ($E=-m$) states, respectively.

5. Special Cases

In this section, we investigate some special cases of physical interest, where we choose the set of $\Lambda$ parameters motivated by the physical potentials singular at $r=R$ that we want to describe, according to the relations (37)-(39). For each case we systematically investigate the confinement properties of the circular potential (i.e., whether it admits bound states or is impermeable), as well as the structure of the resonances and the process of capture and emission of bound states. In all cases, we proceed as follows. First we obtain the values of the potential strengths corresponding to the critical and supercritical states, as well as the secular equation relating these strengths to the bound states energies. Next, we consider complex-energy solutions for purely outgoing scattering states, for various values of the potential strengths, and graphically investigate the structure of the resonances and the discrete set of energies associated with impermeable circular walls at $r=R$. Finally, we investigate the Wigner time delay for the scattering solutions and compare its peaks with the resonant energies calculated using the complex-energy method.

Next we will consider five types of contact potentials concentrated at $r=R$, namely: (i) a purely scalar shell, (ii) a purely electrostatic shell, (iii) a purely magnetic shell, (iv) a “$\delta$ shell", ad (v) a “${\delta }^{\prime }$ shell”.

5.1. A Purely Scalar Shell Potential

First, we consider the case of a purely scalar potential, $A_0=A_r=A_{\theta }=0$ and B arbitrary which models an infinite “kick" in the (anti)particle mass at $r=R$. From (36), the parameters of the $\Lambda$ matrix are

$$\phi =0,a=d={\displaystyle \frac{4+B^2}{4-B^2}},b=c={\displaystyle \frac{4B}{4-B^2}}.$$

(68)

The above relations are invariant under the transformation $B\to {\displaystyle \frac{4}{B}}$, except for a global sign, which does not affect the analysis of the critical, supercritical, bound states, or resonances. Then, for our purposes, it is sufficient to study the behavior of the system when $B\in [-2,2]$. If $B=\pm 2$, we have an impermeable wall, according to (34); in such cases, the boundary conditions at $r=R$ become ${\varphi }_{2,i}\left(R\right)=s{\varphi }_{1,i}$ and ${\varphi }_{2,0}\left(R\right)=-s{\varphi }_{1,o}$, with $s=\mathrm{sign}\left(B\right)$ (see [31]).

5.1.1. Critical, Supercritical and Bound States

The condition (52) for critical states reads

$${\displaystyle \frac{c}{a}}={\displaystyle \frac{4B}{4+B^2}}=-{\displaystyle \frac{l}{mR}}.$$

(69)

For $l=0$ there is only the free solution $B=0$. For $l>0$ the solutions are

$$B_{\mathrm{crit}}=-{\displaystyle \frac{2}{l}}\left(mR\pm \sqrt{m^2{R}^2-l^2}\right).$$

(70)

Thus, a necessary condition for having critical states is $0<l\le mR$. These two solutions are related to each other by the transformation $B\to {\displaystyle \frac{4}{B}}$. On the other hand, the condition (55) for supercritical states is

$${\displaystyle \frac{b}{a}}={\displaystyle \frac{4B}{4+B^2}}=-{\displaystyle \frac{l+1}{mR}},$$

which is similar to the condition for critical states, changing $l\to l+1$. We will have supercritical states only for $0\le l\le mR-1$, and for the following two values of B

$$B_{\mathrm{super}}=-{\displaystyle \frac{2}{l+1}}\left(mR\pm \sqrt{m^2{R}^2-{(l+1)}^2}\right),$$

which are also related between themselves by the transformation $B\to {\displaystyle \frac{4}{B}}$, as they should be. From the previous results we observe that if for a value of B we have a critical state for the angular momentum l, for this same value of B we will have a supercritical state for the angular momentum $l-1$.

The secular equation (61) for the bound states, with (68), becomes:

$${\displaystyle \frac{a}{c}}={\displaystyle \frac{4+B^2}{4B}}=R(E-m)I_{l+1}\left(qR\right)K_{l+1}\left(qR\right)-R(E+m)I_l\left(qR\right)K_l\left(qR\right).$$

(71)

The relationship between the bound state energies and the strength of the pure scalar singular field B is shown in Figure 1, where it is observed that emission/absoprtion

of bound states at the critical and/or supercritical energies is only allowed for angular momenta in the range $0\le l\le mR$, as mentioned above. For $l>mR$, no bound states exist, regardless of the value of the scalar field strength B.

5.1.2. Resonances

The complex energy solutions for the secular equation (71) are shown in Figure 2 as colored dots, for various positive values of the scalar strength B in the interval $[-2,2]$, for $l=1$, $m=2$ and $R=1$. The blue curve represent the locus of the complex solutions to the imaginary part of equation (71), which does not depend on the strength B.

For $B=\pm 2$ the circular wall becomes impenetrable (with boundary conditions depending on the sign of B, as seen above), and all complex energy solutions become real, and correspond to the set of discrete energies admissible for an (anti)particle confined in the inner circle, except, in the case $B=-2$, for the point (in brown in the figure) with energy $0<E_b<+m$, which corresponds to a bound state confined in the outer region of the circle (this state vanishes identically in the inner region, and is decaying in the outer region, as can be seen using the boundary conditions for this case). In the figure we observe that as $B>0$ and decreases from 2 to 0 the complex energies move away from the real axis along the blue U-shaped curves. On the other hand, as $B\in [-2,0)$ increases in absolute value, the complex solutions move towards the real axis, with the absorption of the first bound state at $E=+m$ when B reaches the value $B_{crit}\in [-2,2]$ (red dots in the figure).

5.1.3. Wigner Time Delay

Equation (66), in terms of the parameter B of (68), becomes

$$\begin{array}{ccc}\hfill tan{\delta }_l& =& \left\{2\pi BR\left[(E-m)J_{l+1}^2\left(R\sqrt{E^2-m^2}\right)-(E+m)J_l^2\left(R\sqrt{E^2-m^2}\right)\right]\right\}\hfill \\ & & \times \left\{B^2+2\pi BR\left[(E-m)J_{l+1}\left(R\sqrt{E^2-m^2}\right)Y_{l+1}\left(R\sqrt{E^2-m^2}\right)\right.\right.\hfill \\ & & {\left.\left.-(E+m)J_l\left(R\sqrt{E^2-m^2}\right)Y_l\left(R\sqrt{E^2-m^2}\right)\right]+4\right\}}^{-1},\hfill \end{array}$$

from which the Wigner time delay is obtained from (67).

Figure 3 shows this time for some values of the parameter B. We observe very good agreement between the location of the sharp peaks (or valleys, in the case of negative energies) of the Wigner time delay and the real parts of the complex energies that solve (71). These plots also confirm that the negative inverse of the imaginary parts of the complex energies provides a scale for the Wigner time delay at the resonances.

5.2. A Purely Electrostatic Shell Potential

A purely electrostatic shell at $r=R$ will have $B=A_r=A_{\theta }=0$, and an arbitrary $A_0$. From (36), the parameters of the matrix $\Lambda$ in this case are

$$\phi =0,a=d={\displaystyle \frac{4-A_0^2}{4+A_0^2}},c=-b={\displaystyle \frac{4A_0}{4+A_0^2}}.$$

(72)

We note that the transformation $A_0\to -{\displaystyle \frac{4}{A_0}}$ converts $a,b,c,d\to -a,-b,-c,-d$, which is equivalent to multiplying the matrix $\Lambda$ in (35) by $-1$. As already mentioned, multiplying $\Lambda$ by a phase factor (in fact, by any constant factor) does not affect the results for bound states, resonances, or phase delay. Thus, it is sufficient to investigate those quantities for positive values of $A_0$. We also note that there is no value of $A_0$ for which the circular wall is impermeable, since the permeability condition (34) always holds for a purely electrostatic potential. In the limit $|A_0|\to \infty$ the circular wall becomes completely transparent, since $\Lambda \to -⊮$, which is a phase times the matrix $\Lambda$ of the free case.10

5.2.1. Critical, Supercritical and Bound States

The equation (52) for critical states, written in terms of the electrostatic strength $A_0$ is $-l{A}_0^2+4mR{A}_0+4l=0.$ For $l=0$, this equation has only the trivial solution $A_0=0$, which corresponds to the absence of interaction (free case). For $l>0$, we have two solutions,

$$A_0={\displaystyle \frac{2mR}{l}}\left(1\pm \sqrt{1+{\displaystyle \frac{l^2}{m^2{R}^2}}}\right),$$

(73)

which are related to each other by the transformation $A_0\to -{\displaystyle \frac{4}{A_0}}$. Similarly, the equation (55) for the supercritical states is $-(l+1)A_0^2-4mR{A}_0+4(l+1)=0,$ whose solutions, for $l\ge 0$, are

$$A_0=-{\displaystyle \frac{2mR}{l+1}}\left(1\pm \sqrt{1+{\displaystyle \frac{{(l+1)}^2}{m^2{R}^2}}}\right),$$

(74)

that also shows the symmetry $A_0\to -{\displaystyle \frac{4}{A_0}}$.

The secular equation (61) for the bound states is written as

$${\displaystyle \frac{4-A_0^2}{4A_0}}=-R(E-m)I_{l+1}\left(qR\right)K_{l+1}\left(qR\right)-R(E+m)I_l\left(qR\right)K_l\left(qR\right),$$

(75)

expression in which we observe that any energy in the interval $(-m,m)$ will be associated to a single bound state energy for exactly two values of $A_0$. Figure 4 shows the relationship between the bound state energies $E_b$ and the value of the electric potential strength $A_0$, for three values of the angular momentum l.

5.2.2. Resonances

We now look for complex energy solutions to the secular equation (75), which can be rewritten in a convenient form as:

$$(E+m)J_l\left(pR\right)H_l^{\left(1\right)}\left(pR\right)+(E-m)J_{l+1}\left(pR\right)H_{l+1}^{\left(1\right)}\left(pR\right)={\displaystyle \frac{i}{2\pi R}}{\displaystyle \frac{4-A_0^2}{A_0}}.$$

(76)

Figure 5 shows by means of colored dots the complex energies that solve (76), for various positive values of the electrostatic strength $A_0>0$ and $m=2,R=1,l=2$.

The blue curve represents the solutions for the real part of (76), which does not depend on the potential intensity $A_0$. As $A_0$ increases in the positive direction, the complex energy solutions move along the blue curve from left to right. When $A_0$ crosses the value (74) with the minus sign ($A_0\approx 1.07$ in the figure), a bound state is captured (green dots in the figure) at $E=-m$; when $A_0$ crosses the value (73) with the plus sign ($A_0\approx 4.828$ in the figure), a bound state is emitted at $E=+m$ (brown dots in the figure).

5.2.3. Wigner Time Delay

The equation (66), for the pure electrostatic case, becomes

$$tan{\delta }_l={\displaystyle \frac{2\pi A_{0}R\left[(E+m)J_l{\left(pR\right)}^2+(E-m)J_{l+1}{\left(pR\right)}^2\right]}{A_0^2+2\pi A_{0}R\left[(E+m)J_l\left(pR\right)Y_l\left(pR\right)+(E-m)J_{l+1}\left(pR\right)Y_{l+1}\left(pR\right)\right]-4}},$$

from which we can calculate the Wigner time delay using (67), shown in Figure 6 for some values of the parameter $A_0$. Again, we can observe a very good agreement between the sharp peaks (or valleys if the energies are negative) of the Wigner time delay and the real parts of the complex energies that solve (76).

5.3. A Purely Magnetic Shell Potential

In this section we will consider a general pure magnetic potential ($B=A_0=0$), where $A_r$ and $A_{\theta }$ are not simultaneously zero. From (36), the parameters of the $\Lambda$ matrix are

$$\phi ={tan}^{-1}\left({\displaystyle \frac{4A_r}{A_r^2+A_{\theta }^2-4}}\right),a={\displaystyle \frac{1}{d}}=\pm {\displaystyle \frac{-A_r^2-{(A_{\theta }+2)}^2}{\sqrt{16A_r^2+{(A_r^2+A_{\theta }^2-4)}^2}}},b=c=0,$$

(77)

As mentioned above, the value of the phase $\phi$ and the choice of the sign in (77) are irrelevant to determining the energies of the critical, supercritical, bound states, and resonances, so we will only consider the negative sign in the equation above, which implies that $a={\displaystyle \frac{1}{d}}\ge 0$. All the analysis of the magnetic potential will be carried out in terms of the parameter a, but, for clarity, let us consider some particular limits and their relationships with the physical intensities $A_r$ and $A_{\theta }$. From (77) we have that $a=0$ (resp. $a=\infty$) if, and only if, $A_r=0$ and $A_{\theta }=-2$ (resp. $A_{\theta }=+2$), which are the cases for an impermeable circular barrier, according to (34). The case $a=1$ (the free case) corresponds to the limit when the strength $A_r$ or $A_{\theta }$ (or both) tend to infinity.

5.3.1. Critical, Supercritical and Bound States

The condition (52) for critical states in this case reads $a{\displaystyle \frac{l}{mR}}=0$. For $l=0$ any value of a admits a critical state. For $l>0$ only an impenetrable wall with $a=0$ admits a critical state, which will be a state for a particle confined in the inner region of the circular wall, as can be seen by using the boundary conditions in this case, ${\varphi }_{2,i}\left(R\right)={\varphi }_{1,o}\left(R\right)=0$. On the other hand, the condition (55) for supercritical states is $a{\displaystyle \frac{l+1}{mR}}=0$ and, again, only an impermeable barrier with $a=0$ will admit a supercritical state. In this case, however, the supercritical state will not correspond to an antiparticle confined in $r<R$, but in the external region $r>R$. The secular equation (61), with parameters (77), becomes:

$$\begin{array}{c}\hfill a{I}_l\left(qR\right)K_{l+1}\left(qR\right)+{\displaystyle \frac{1}{a}}{I}_{l+1}\left(qR\right)K_l\left(qR\right)=0\iff a^2=-{\displaystyle \frac{I_{l+1}\left(qR\right)K_l\left(qR\right)}{I_l\left(qR\right)K_{l+1}\left(qR\right)}}.\end{array}$$

(78)

The above equation has no solution for bound state energies $\left|E\right|<m$, for any $a>0$, since $a^2$ is real and nonzero, and the right-hand side of the above equation is always negative. This result can be extended to the cases $a\to +\infty$. For the case $a=0$ we arrive at the same conclusion, explicitly considering the boundary conditions for the impermeable barrier in this case (see above). In summary, a purely magnetic potential does not allow solutions with $\left|E\right|<m$. However, as we will see later, an impenetrable wall with $a=0$ will admit a set of discrete energies with $E<-m$ or $E\ge m$ for an (anti)particle confined in the inner circle $r<R$. Similarly, an impenetrable wall with $a=+\infty$ will have states confined in the inner circle for a discrete set of energies $\left|E\right|>m$, for $l>0$, and $\left|E\right|\ge 0$, for $l=0$, as can be concluded using the appropriate boundary conditions in this limiting case, namely ${\varphi }_{1,i}\left(R\right)={\varphi }_{2,o}\left(R\right)=0$.

5.3.2. Resonances

Figure 7 shows the complex energies solving the equation (78), for various values of the parameter $a\ge 0$, and for $l=4$. As mentioned above, the case $a=0$ (dark blue points in the figure) corresponds to an impermeable circular barrier with $A_r=0$ and $A_{\theta }=-2$. In this case, the purely real solutions of (78) correspond to the admissible energies for an (anti)particle confined in the inner region of the circular impenetrable wall (this set includes $E=+m$, but not $E=-m$, which corresponds to the energy of a particle confined in the outer region).

The two complex solutions close to $\left|E_R\right|=m$ correspond to remnant resonances. The case $a=\infty$ (purple dots), corresponding to $A_r=0,A_{\theta }=+2$, also produces an impenetrable circular wall (with boundary conditions on the wall different from those in the case $a=0$) and, moreover, presents a set of purely real energies corresponding to the set of discrete energies admissible for an (anti)particle confined in the inner region of the impenetrable wall. The system also presents two complex solutions, corresponding to remnant resonances close to $|E_R|=m$. We note (not shown in the figure) that, for lower values of the angular momentum ($l=0$ and $l=1$, with the remaining parameters fixed), such remnant complex resonances at the limits $a=0$ and $a=\infty$ do not appear. In both cases of impenetrable walls, the solutions in the outer region admits a continuum of energies $\left|E\right|>m$, as expected. For $a=0$, $E=+m$ does not belong to this continuum (it vanishes identically in the outer region), but $E=-m$ does (the solution with $E=-m$ vanishes trivially in the inner region).

5.3.3. Wigner Time Delay

For the purely magnetic case, equation (66) in terms of the parameter a of (77) is

$$tan{\delta }_l={\displaystyle \frac{\left(a^2-1\right)J_l\left(pR\right)J_{l+1}\left(pR\right)}{a^2{J}_l\left(pR\right)Y_{l+1}\left(pR\right)-J_{l+1}\left(pR\right)Y_l\left(pR\right)}}.$$

The Wigner time delay (67) is shown in Figure 8 for the same parameters as in Figure 7. For large values of the resonant energies (real parts of the complex energies that solve the second equation in (78)), we observe a very good agreement between these energies and the locations of the sharp peaks (or valleys for negative energies) of the Wigner time delay. For these resonances, the plots also confirms that the negative inverse of the imaginary parts of the complex energies gives a scale for the Wigner time delay. In the figure, the agreement between the real and imaginary parts of the complex energies and respectively the location and magnitude of the Wigner time delay peaks is poorer for low values of $E_R$, and these relationships are absent for the two complex resonances close to $E_R=m$. The lack of such correspondence may be understood from the fact that such resonances do not evolve to a confined state for the corresponding impenetrable barrier in the limits $a\to 0,+\infty$.

5.4. A “Delta Shell" Potential

The delta shell potential corresponds to an equal mixture of the scalar and the electrostatic potentials, with the magnetic potentials vanishing (arbitrary $A_0=B$, $A_r=A_{\theta }=0$). We call it a “delta shell" because the matrix $\Lambda$ in this case reduce to that associated to a delta potential in the non-relativistic limit of one-dimensional relativistic quantum mechanics (see [31]). From (36), the $\Lambda$ parameter values in this case are $a=d=1$, $\phi =b=0$, and $c=2A_0=2B$. From now on, we will characterize the strength of the delta shell potential by the value of $A_0=B$.

5.4.1. Critical, Supercritical and Bound States

For each l, condition (52) gives a critical state ($E=m$) for $A_0=B=-{\displaystyle \frac{l}{2mR}}$. For $l=0$ this condition degenerates to the free case (no interaction). For $l=1$ we have a quasi-bound critical state and for $l>1$ we have a true critical bound state. On the other hand, for $l\ge 0$ equation (55) has no solution for supercritical states. The secular equation (61) for the bound state energies simplifies to

$$I_l\left(qR\right)\left[q{K}_{l+1}\left(qR\right)+c(m+E)K_l\left(qR\right)\right]+q{I}_{l+1}\left(qR\right)K_l\left(qR\right)=0,q=\sqrt{m^2-E^2},$$

and then

$$A_0=B=-{\displaystyle \frac{1}{2R(E+m)I_l\left(qR\right)K_l\left(qR\right)}}.$$

(79)

Figure 9 shows the relationship between the potential strength $A_0=B$ and the bound state energies $E_b$ for three values of the angular momentum l. From this figure, and also from (79), it follows that a necessary condition for having bound states is that $A_0=B<0$, since $-m<E<m$ and $I_l\left(qR\right)K_l\left(qR\right)>0$. At the value $A_0=B=-{\displaystyle \frac{l}{2mR}}$, a bound state is captured/emitted at the critical energy $E=m$. There exits a unique bound state for $A_0=B\in (-\infty ,{\displaystyle \frac{-l}{2mR}}]$.

5.4.2. Resonances

We now look for the complex energy solutions to the secular equation (79), which now becomes a complex equation that can be conveniently rewritten as

$$(E+m)J_l\left(pR\right)H_l^{\left(1\right)}\left(pR\right)={\displaystyle \frac{i}{\pi R{A}_0}}.$$

(80)

Figure 10 shows in blue the contour line corresponding to the complex energies that solve the real part of the equation (80), which does not depend on the potential strength $A_0=B$, for $m=2,R=1$ and $l=2$. The colored points on that curve correspond to the complex energies that solve the whole complex equation, for various values of the potential strength $A_0=B$. The real part of the complex energies corresponds to the resonant energies. As $A_0=B$ increases in absolute value, the complex energies on the blue curve shift toward the real axis. When $A_0=B=-{\displaystyle \frac{l}{2mR}}$ ($=-0.5$ in the figure, denoted by red dots), the system captures/emits the unique bound state from/to the continuum. The bound state energy moves leftward along the real axis with the strength $A_0=B$ increasing in absolute value, and asymptotically approaches $E=-m$ when $A_0=B\to \pm \infty$. In these limits the circular wall becomes impermeable and all the complex energies being real. These energies, except $E=-m$, correspond to the discrete energies admissible for an (anti)particle confined in the inner region of the impenetrable circular wall. In this limit, there exits a supercritical energy $E=-m$, which corresponds to an antiparticle restricted to the outer region of the impenetrable circle, as can be inferred from the associated boundary conditions in these cases, namely ${\varphi }_{1,i}\left(R\right)={\varphi }_{1,o}\left(R\right)=0$ in (40)).

5.4.3. The Wigner Time Delay

Now, the equation (66) for the phase change becomes

$$tan{\delta }_l={\displaystyle \frac{\pi R(E+m){\left[J_l\left(Rp\right)\right]}^2{A}_0}{\pi R(E+m)J_l\left(Rp\right)Y_l\left(Rp\right)A_0-1}},$$

(81)

from which we can obtain the Wigner time delay by using (67), which is shown in Figure 11 for some values of $A_0=B$.

We also plot the real part and the negative inverse of the imaginary part of the complex energies obtained as solutions of (80). We note that the real part of the complex energies (the resonant energies) are very good approximations for the location of the Wigner time delay peaks. These results show that both methods for finding resonant energies (by calculating the real part of the complex energies associated with purely outgoing scattering states or the energies related to the Wigner time delay peaks) give essentially the same results. Moreover, the negative inverse of the imaginary part of the complex energies provides a scale for the Wigner time delay at these resonances. In the absorption/emission of a bound state at $E=+m$ the Wigner time delay diverges, which is consistent with the $-{\displaystyle \frac{1}{E_I}}$ scaling, since $E_I=0$ in the emission/absorption of a bound state.

5.5. A “${\delta }^{\prime }$ Shell" Potential

As a final special case, we will analyze the “${\delta }^{\prime }$ shell", which is an inverse combination of the scalar and electrostatic potentials, i.e., $B=-A_0,A_r=A_{\theta }=0$. This choice results in the following values for the $\Lambda$ matrix parameters: $a=d=1$, $c=\phi =0$ and $b=-2A_0=2B$. It is so named because, in the non-relativistic limit of the analogous one-dimensional case, this choice of parameters results in to the so-called ${\delta }^{\prime }$ (or the non-local delta prime) interaction.

In the following, we will show that, by means of a “spin flip" transformation, the Dirac equation (33) for a “${\delta }^{\prime }$ shell" can be obtained from that for a “$\delta$ shell" potential. Therefore, the results for the present case can easily be obtained from the previous one and will not be presented here. Let us then consider the Dirac equation (33) for a delta shell potential, $A_0=B$, $A_r=A_{\theta }=0$, and a “spin flip" transformation $\Phi \to -i{\gamma }^1\Phi$. It is evident that this equation will be invariant under this transformation if also

$$E\to -E,B\to B,A_0\to -A_0,l\to -(l+1).$$

Therefore, the transformed equation will be the Dirac equation for a ${\delta }^{\prime }$ shell, with $B=-A_0$, $A_r=A_{\theta }=0$.

6. Discussion and Conclusions

In this work, we investigate the massive Dirac equation in the plane, with a singular interaction supported on a circumference of radius R. Using the radial symmetry of the interaction, we show that this problem can be reduced to a one-dimensional point interaction problem in the radial coordinate. To study this problem, we adapt the distributional approach of [31], in which a physical interpretation of the parameters of point interactions in one dimension was given in terms of the strengths of four Lorentz singular point potentials. In this way, we were able to interpret the four parameters of the singular interaction supported on the circumference as the strengths of a scalar (B) and the three polar components ($A_0,A_r$ and $A_{\theta }$) of a Lorentz vector. It is worth noting that the usual requirement of boundedness of the Dirac spinor at the origin (sometimes stated without justification) arises from the basic assumptions of the distributional approach.

After obtaining time-independent solutions for critical, supercritical, bound, and scattering states for the most general contact potential concentrated on $r=R$, we systematically addressed five special cases of contact interactions and their confining properties. In all these cases, we investigated how varying the potential strengths modify the structure of the critical, supercritical, bound, and resonant states. We also investigated the values (or corresponding limits) of the strengths for which the circular barrier is impenetrable. The resonant energies were identified using the real parts of the complex energy solutions for purely outgoing scattering conditions, and subsequently compared to the location of the Wigner time delay peaks, with excellent overall agreement. Our results also showed that, when the real parts of the complex energies coincide with the location of the Wigner time peaks, the negative inverse of the imaginary parts provides a scale for the intensities of these peaks. From the permeability condition (34), we observed that only potentials with strengths $A_r=0$ and a non-zero B or $A_{\theta }$ can produce an impenetrable barrier on the circumference. For example, we saw that a purely electrostatic contact barrier cannot completely confine an (anti)particle either in the inner or outer circle. This result is consistent with that of reference [35] for the one-dimensional case, and recalls Klein’s paradox for point potentials. In the following, we briefly summarize the main results for the five special cases considered in this article.

A purely scalar shell (arbitrary B, $A_0=A_r=A_{\theta }=0$), modeling an infinite “kick" in the (anti)particle mass at $r=R$, displays a rich structure. The ensemble of bound states and resonances has the symmetry $B\to {\displaystyle \frac{4}{B}}$, implying that an infinite strength B will be similar to the free case. This potential can produce an impermeable wall at $r=R$ if the scalar strength is $B=\pm 2$ and can admit none, one, or two bound states (for $l>mR,l=mR-1$ and $l<mR-1$, respectively), for suitable negative values of the scalar strength B. The larger the mass, the more values of angular momentum will admit bound states, which may be absorbed/emitted at both the critical and supercritical energies, depending on the angular momentum l. For any value of B, there exits a discrete set of resonances close to the discrete energies admitted for (anti)particles completely confined by an impenetrable circular wall. The energy location of the Wigner time delay peaks agrees very well with the real part of the complex resonances.

A purely electrostatic shell (arbitrary $A_0$, $B=A_r=A_{\theta }=0$) is never impermeable, for any finite or infinite value of strength $A_0$. For positive or negative $A_0$ and in suitable ranges, it admits a single bound state that can be absorbed/emitted at both critical and supercritical energies $E=\pm m$. The potential also exhibits a discrete set of resonances for any value of strength $A_0$. The real part of the complex resonances also shows a very good agreement with the Wigner time delay peaks.

A purely magnetic shell (arbitrary $A_r$ and $A_{\theta }$ and $B=A_0=0$) does not admit bound states $\left|E\right|<m$, but it can be impenetrable for suitable combinations of finite values of the electromagnetic strengths, namely $A_0=0$ and $A_{\theta }=\pm 2$. This singular potential can therefore completely confine an (anti)particle in the inner region of the circular wall at $r=R$, with a discrete set of admissible energies depending on the value of $A_{\theta }$ ($=\pm 2$). This potential also allows a discrete set of resonances for arbitrary choices of the electromagnetic strengths, with energies very close to the discrete energies of the impenetrable cases. It is an interesting result that for large values of the angular momentum there are still complex resonances in the limit of impenetrable walls, which do not correspond to the location of any peaks in the Wigner time delay. This mismatch may be due to the fact that these resonances do not converge to a confined state within this limit. For high values of the resonant energies, there is a good agreement between the real part of the complex energies and the locations of the Wigner time delay peaks.

Finally, we considered a situation with two “mixed" potentials, a “$\delta$ shell" potential (arbitrary $A_0=B$, $A_r=A_{\theta }=0$) and a “${\delta }^{\prime }$ shell" (arbitrary $A_0=-B$, $A_r=A_{\theta }=0$). These two cases are related to each other by a “spin flip" transformation, and we detail only the “delta shell" case. The impermeability condition is met only in the limits $A_0=B\to \pm \infty$, and for both cases the system admits the same discrete set of energies for a particle confined in the inner region of the impenetrable circular wall. For finite values of the strength $A_0=B$ this potential admits a single bound state for sufficiently large and negative values of this parameter. It also admits a discrete set of resonances for any value of the strength $A_0=B$, with energies close to those corresponding to the impenetrable case. The single bound state is captured/emitted only at the critical energy $E=+m$. The results obtained for Wigner time delay are similar to those analyzed in the previous cases.

Future applications of the contact models presented here could include the analysis of the Aharonov-Bohm effect with a singular magnetic field concentrated on a circumference. This investigation may be interesting due to the nontrivial characteristics of singular electromagnetic interactions, both in one and two dimensions. It also seems necessary and interesting to consider the natural generalization of these models to three dimensions, with singular potentials concentrated on a spherical surface. Finally, we believe it is worth emphasizing that our two-dimensional models presented here can describe phenomena in Dirac materials, as mentioned in the Introduction, where the Fermi velocity at the surface acts like the speed of light, and may therefore be accessible, at least in principle, to perform experiments involving ordinary velocities and short-range potentials approximately localized on a circle.