Scattering of charged fermion to two-dimensional wormhole with constant axial magnetic flux

Scattering of charged fermion with $(1+2)$-dimensional wormhole in the presence of constant axial magnetic flux is explored. By extending the class of fermionic solutions of the Dirac equation in the curved space of wormhole surface to include normal modes with real energy and momentum, we found a quantum selection rule for the scattering of fermion waves to the wormhole. The newly found {\it momentum-angular momentum relation} implies that only fermion with the quantized momentum $k=m'/a\sqrt{q}$ can be transmitted through the hole. The allowed momentum is proportional to an effective angular momentum quantum number $m'$ and inversely proportional to the radius of the throat of the wormhole $a$. Flux dependence of the effective angular momentum quantum number permits us to select fermions that can pass through according to their momenta. A conservation law is also naturally enforced in terms of the unitarity condition among the incident, reflected, and transmitted waves. The scattering involving quasinormal modes~(QNMs) of fermionic states in the wormhole is subsequently explored. It is found that the transmitted waves through the wormhole for all scenarios involving QNMs are mostly suppressed and decaying in time. In the case of QNMs scattering, the unitarity condition is violated but a more generic relation of the scattering coefficients is established. When the magnetic flux $\phi=mhc/e$, i.e., quantized in units of the magnetic flux quantum $hc/e$, the fermion will tunnel through the wormhole with zero reflection.


I. INTRODUCTION
The physics of fermion on 2-dimensional surface has been one of the major research topics in recent years. Interesting boundary phenomena that have no bulk analogs emerge, e.g. quantum (spin) hall effects, topological matters, and more recently graphene physics [1]. Two dimensional surface can be manipulated in the laboratory so that it can be curved, strained, and twisted. A cage structure of graphene wormholes, the schwarzite, can even be produced with promising properties [2,3]. The gauge fields such as electric [4] and magnetic fields [5] can be applied to the system, leading to unique intriguing 2D phenomena, notably the well-known Landau quantization of fermionic states on a plane. It has been shown that effects of strain and gauge fields could be similar in 2D [6,7]. Behaviour of quantum particle on the curved surface in the presence of gauge fields can be considerably different from the flat situation [8,9]. Similar to the strain, curvature effects [10][11][12] can mimic gauge fields [13,14], specifically curvature and gauge connection appears in the equation of motion with equal role. For example, fermions on 2-dimensional sphere and wormhole experience spin-orbit coupling induced from the surface curvature [15][16][17][18] even in the absence of the gauge fields. Addition of axial magnetic field generates Landau quantization distinctive from the planar case. For the wormhole, the fermionic states can be in quasinormal modes with complex energies, its quantum statistics can be altered by the magnetic flux through the hole [18,19].
Dirac fermion in graphene wormhole without gauge field has been discussed in Ref. [20,21]. In Ref. [18], we investigate the effects of axial magnetic field on a charged fermion in a (1 + 2)-dimensional wormhole. This 2-dimensional wormhole is fundamentally different from the (1 + 3)-dimensional wormhole in General Relativity (see e.g. [22]), there is no time dilation in the 2-dimensional wormhole under consideration. For the constant magnetic flux scenario, the system can be solved analytically and exact solutions are found to contain "normal" (real energy but complex momentum) and quasinormal modes (QNMs).
In this work the scattering of fermions to the wormhole is explored. The more generic solutions to the equation of motion of the fermion in the magnetized wormhole are constructed in terms of hypergeometric functions. The normal modes are actually found when the (effective) angular quantum number m (see definition below) is related to the momentum of the fermion by k = m a √ q ≡ k m . Such quantized momentumangular momentum relation is unique to the 2D wormhole under consideration. In Section II, the mathematical formulation of the wormhole and fermion in curved space is established. In Section III as a review of the main results of Ref. [18], the Dirac equation in the magnetized wormhole is written and solved analytically, then the general solutions in the upper and lower plane connected to the wormhole is discussed. Matching conditions of the scattering of fermionic waves to the wormhole is considered in Section IV. Section V discusses the use of Wronskian to derive a general relation between scattering coefficients. The scattering scenarios are categorized into normal-modes and QNMs scattering and some of the results are presented in Section VI. Section VII concludes our work.

II. GEOMETRIC AND GAUGE SETUP OF THE WORMHOLE
FIG. 1: Geometric structure of a wormhole surface where a is a radius at the radius function R(u = u0 = ln q/2). And r is the radius of curvature of the wormhole surface along u direction.
In the (1+2)-dimensional curved spacetime, the line element on the surface of wormhole can be cast in the following form (2.1) The fermion will experience the effective curvature that can be addressed by considering the Dirac equation in curved spacetime where Ψ = Ψ(t, u, v) represents the Dirac spinor field on the wormhole and M represents the rest mass of the particle, c is the speed of light, e is electric charge, and A µ is the electromagnetic four-potential. The γ a are the Dirac matrices given by where σ k are the Pauli matrices defined by They obey the Clifford algebra The Pauli matrices have a useful identity that we will use later where ijk is Levi-Civita symbol. The covariant derivative of the spinor interaction with gauge field in the curved space is given by where the spin connection Γ µ [23] is where β, µ, ν ∈ {t, u, v} and the Christoffel symbols Γ βµν are defined by Then for the metric (2.1), and zero otherwise. It is then straightforward to show that the spin connections are We will apply an external magnetic field such that the z-component B z = B(z) is uniform with respect to the plane (x, y) and the magnetic flux through the circular area enclosed by the wormhole at a fixed z is constant, namely B z ∼ 1/R 2 . Due to the axial symmetry, the electromagnetic four-potential can be expressed in the axial gauge as and in the wormhole coordinates as x, y, z) = 0, 0, The magnetic field is then given by (2.10) We will now consider the Dirac equation of fermion in the wormhole in the presence of constant magnetic flux.

III. THE DIRAC EQUATION IN MAGNETIZED WORMHOLE
Utilizing the results from above equations, the Dirac equation Eq.(2.2) can be written in the form then where D is a differential operator We can define the pseudo-vector potential as Note that the effective gauge potential in the u direction, Aũ, is generated by the curvature along the v direction, Γ v . In this sense, wormhole "gravity" or curvature connection manifests itself in the form of gauge connection in the perpendicular direction. Consider a stationary state of the Dirac spinor needs to be single-valued at every point in spacetime, Ψ(t, u, v) must be a periodic function in v with period v ∈ [0, 2π], in the form From Eq.(3.2), the first term is equivalent to the Dirac operator with the pseudo gauge potential Aũ(u) in the u direction. This term is generated from the change in the radius of the hole. The second term contains gauge potential that generates a spin-orbit coupling. A similar setup has been used to study nanotubes under a sinusoidal potential [20]. Here we consider the dispersion relation for the two-dimensional fermions described by the Dirac equation in the presence of the effective potential arising from the wormhole geometrical structure. In the presence of external magnetic field B = × A along the z direction, the charged fermion moving in v direction is expected to form a stationary state with quantized angular momentum and energy, i.e. the Landau levels in the curved space with hole. To show this, we need to solve for the stationary states of the system. We can start by considering −iD×Eq.(3.5)-E + M c 2 / c ×Eq.(3.6) to obtain For constant magnetic flux, the operator D 2 in the equation of motion now takes the form where m = m − φ φ 0 , and the magnetic flux quantum φ 0 ≡ hc/e. We have used the momentum parameter To be specific, we parametrise the shape of wormhole by setting R(u) = a cosh q (u/r). They are based on a q-deformation of the usual hyperbolic functions defined by [24] Note the relevant properties . (3.11) The deformed functions reduce to hyperbolic functions when q = 1. For this choice, the wormhole will possess two Hilbert horizons at u p,m = r ln q + r 2 a 2 ± r a where R (u p,m ) = ±1. Then define a variable X(u) ≡ rR (u)/a = sinh q (u/r) to obtain , the equation of motion can be rewritten as where we assume 14) The coefficients are defined as the following Depending on the sign choices of κ 1 , κ 2 , the resulting equation of motion and the corresponding energy levels will be dependent or independent of the spin-orbit coupling term ∼ σmr/a √ q.
Define X = −i √ qY , the equation then takes the form Eq.(3.16) is the Jacobi Differential Equation, the energy levels become The solutions to Eq.(3.16) are the Jacobi polynomials [25] Φ n (Y ) = P (α0,β0) for integer n and α 0 = 2α, β 0 = 2β. Finally, the solutions of Eq.(3.9) is where X(u) = sinh q (u/r) and κ 1 , κ 2 = ±1. These are the solutions with QNMs found in Ref. [18]. The momentum k for these solutions is generally a complex quantity, so they are not exactly normal modes even when the energy E is real. As will be shown subsequently in VI A and VI B 2, we can analytically continue the solutions to the general n < 0 cases where the normal modes with real energy and momentum can be found as a special case with n = −1/2. Other quasinormal n < 0 solutions are also relevant in the scattering processes.

A. Solutions in the Upper and Lower plane
In the flat upper (R (u) = 1) and lower (R (u) = −1) plane region outside the wormhole, the equation of motion, Eq.(3.9), takes the form for u ≥ u p in the upper plane, and for u ≤ u m in the lower plane. The Hankel function of the first (second) kind corresponds to the waves propagating in +(−)û direction respectively.

IV. MATCHING CONDITIONS IN THE SCATTERING
The matching conditions of the wave functions are the equality of the complex energy (E) and angular momentum (m ), and the smooth continuity (i.e., C 1 continuity) of the wave functions between the inner and outer region of the wormhole. The momentum k will also be equal due to the relation (3.17).
Consider the incoming waves propagating in the upper plane region scatter with the wormhole. At the upper Hilbert horizon u p , the waves will be partially reflected back and partially transmitted into the inner region of the wormhole. The transmitted waves will be again partially reflected and partially transmitted into the lower plane region at the lower Hilbert horizon u m . In this scenario, there is no incoming waves at the lower Hilbert horizon, i.e., R (−) m,σ = 0. We can normalize I m,σ ≡ 1 and R (+) m,σ ≡ R m,σ where the Hankel function of the first and second kind correspond to outgoing and incoming waves, respectively. In the wormhole region u m < u < u p , the general solution from Eq. (3.19) can be expressed as where X = sinh q (u/r), α = κ 1 α 0 , β = κ 2 α * 0 , and α 0 ≡ 1 4 + i √ q σm r 2a and the condition E n,σ = E n ,σ is required. Note that the parameter α 0 also depends on m, σ.

The matching condition : At upper surface
The boundary between the inner and outer region of the wormhole is at R(u p ) ≡ R p , see Fig. 2. The wave function and the first derivative of the wave function must be continuous at the boundary. From conservation law of energy-momentum, E (in) = E (out) and k (in) = k (out) . The first boundary condition at the upper plane is where sinh q (u p /r) = r/a, and R p = qa 2 + r 2 . The derivative of the wave function at the boundary u = u p The matching condition : At lower surface Even for general q ≥ 0, the wormhole is symmetric with respect to the midpoint u 0 = ln q/2. The Hilbert horizon at the lower surface is at u m where R(u m ) ≡ R m = R p . The second boundary condition takes the form (4.4) The derivative of the wave function at the second boundary u = u m , R = R m is From the matching conditions (4.2)-(4.5), the scattering coefficients A m,n,σ , B m,n,σ , T m,σ , R m,n can be solved for which here and henceforth the subscripts will be suppressed.

V. ANALYTIC RELATION BETWEEN REFLECTION AND TRANSMISSION COEFFICIENTS
In this section a very general analytic relation between the reflection and transmission coefficients are derived from the equation of motions. The relation can be expressed in terms of Wronskian of the solutions in each region. The resulting relation is then verified numerically for each scattering scenario.
The generic form of the equation of motion can be expressed as and 3) respectively. Note that R(u) in the planes is given by (3.21). For any two solutions ϕ 1 , ϕ 2 satisfying (5.1), by multiplying the equation of motion of one solution with the other solution then subtracting the two equations we obtain 4) or in terms of the Wronskian W 12 ≡ ϕ 1 ϕ 2 − ϕ 2 ϕ 1 , By integrating this equation in the region covering the wormhole, is the resulting general relation and the subscript of the Wronskian has been omitted since the relation is valid for any pairs of solutions. The RHS of (5.6) can be calculated explicitly, Since the wormhole is symmetric with R(u m ) = R(u p ) = r 2 + qa 2 , the Wronskian at the two Hilbert horizons are always equal, W (u m ) = W (u p ).
The general relation (5.6) can be applied to the two solutions ϕ and ϕ * , specifically at u m , u p where the wave functions and their first derivatives in the two connecting regions are equal. In terms of the solutions in the outer regions, we have the relation (with the subscripts suppressed) (5.8)

VI. SCATTERING FOR NORMAL AND QUASI-NORMAL MODES
In this section we consider scattering of fermion in 2-dimensional planar surface with (1 + 2)-dimensional wormhole. Even this "wormhole" is not the actual spacetime wormhole as in GR and other gravity theories where time dilatation exists, its scattering still reveals a number of interesting properties. As discussed in details below, scattering of physical fermion with real momentum and energy has quantized behaviour that relates momentum and angular momentum of the scattered fermion. For scattering involving QNMs of the fermion, enhancement (and attenuation) of scattered fermion is a possibility. The wormhole and spin parameters are set to r = a = q = 1 = σ for all of the numerical results in this section.
A. scattering for real momentum k Naturally, incoming waves from outside region are generated with real momentum and energy. When scattered with the wormhole, only real momentum and energy states (i.e., normal modes) will be allowed to propagate through the hole and transmitted to the other outside region. The scattering waves will be partially reflected back and partially transmitted through the hole to the other side. From Eqn.
the unitarity condition. Unitarity is the consequence of reality of momentum in the scattering. Table II shows scattering coefficients for m = 1, 2, 3, 4 cases, all of which the unitarity relation (6.2) is numerically verified. Wave function profile of the matching is shown in Fig. 6(a). The m -dependence of the transmission coefficient T and R are shown in Fig. 3. Interestingly, T (m ) and R(m ) converge to oscillating functions for large m . The unitarity condition (6.2) is always obeyed.

B. Scattering involving QNMs
Scattering involving QNMs usually violates unitarity since the waves are decaying in time. It will satisfy a more general relation given by (5.8). There are three possible scenarios for such scattering. There are three scenarios for scattering involving QNMs of the wormhole. The matching condition requires the energy and inevitably the momentum k to be the same value for inner and outer regions of the wormhole. Generically since the momentum given by (3.17) is always a complex quantity with an exception of k m discussed in Sect. VI A, the momentum of wave functions in the outer regions will also need to be complex as well. Consequently, there are spatial attenuation and enhancement of waves in the outer regions as the waves are concurrently decaying in time. The physical energy to be measured in experiments is the real part of the (complex) energy used in the matching conditions.

scattering with m = 0
The scattering in this scenario corresponds to the states in the wormhole with quantized flux φ = mφ 0 , m = 1, 2, 3, ... and the orbital angular momentum of the states is m . The quantized magnetic flux is required specifically for the scattering to occur. In this case the energy could be real for massive fermion and small n. On the other hand, the momentum becomes k = i n + 1 2 /r, purely imaginary. The waves in the outer regions to be matched with the ones inside the wormhole thus need to be attenuating or enhancing along the u direction. Since the source of the fermion can locate at finite distance from the hole in experimental situation, this scattering scenario is still physically relavant. The fermionic states inside the wormhole consist of two waves with n = n , ϕ(+, −, n) and ϕ(−, +, n) where (κ 1 , κ 2 ) = (+, −), (−, +) represents +û, −û going modes respectively. The scattering coefficients for n = 0, 1, 2, 3 are presented in Table III. The transmission and reflection coefficients obey general relation (5.8) but not the unitarity condition (6.2) due to the complexity of momentum k. Remarkably, the waves simply tunnel through the wormhole with no reflection, i.e., R = 0. The scattering coefficients for n = 0, 1, 2, 3 are given in Table III. Wave function profile of the matching is shown in Fig. 6(b).

scattering with nonzero m
For m = 0 the waves inside the wormhole consist of ϕ(+, −, n) and ϕ(−, +, n ). According to the energy formula (3.17), the −û-travelling wave ϕ(−, +, n ) must have the wave function associated with n = −n − 1 in order to have the same energy, E n = E n . Since the Jacobi polynomials P n are zero for non-positive integer n, we need a more general form of the solutions. From the equation of motion (3.16) with (3.17) substituted, we can analytically continue the solution to the case where n < 0 by solving this equation as the hypergeometric equation. There are two solutions, where γ = −4n(n + 2α + 2β + 1) + (2α + 2β + 1) 2 . (6.6) The first solution has kinks or sharp turns in the inner region of wormhole while the second solution is smooth throughout. Therefore we use only the second solution (6.5) in the scattering analysis. The incoming waves are scattered by the wormhole, thus excite the fermionic modes in the process. The reflected outcomes occurred at both Hilbert horizons where u p and u m . At the same time, some of the waves in the inner region leaks out through both horizons into the outer planes. The leaking waves are attenuated or enhanced depending on the sign of Im(k). For the −û (+û) moving H ν ), the wave function is spatially enhanced (attenuated) with respect to the increase of u respectively. Wave function profile of the matching is shown in Fig. 6(c). The scattering coefficients for certain parameters are given in Table IV. The m -dependence of the transmission coefficient T and R are shown in Fig. 4. Remarkably, both T and R show a peak at m = ±6.
Scattering of waves to the wormhole with QNMs excitation is analogous to interaction of oscillating free springs with damped spring. The springs in all regions will eventually decay to zero due to the imaginary part of QNMs in the factor e −iEt for Im(E) < 0. For QNMs with Im(E) > 0, the system will be unstable, the backreaction to the curve surface from the fermion will be large.

scattering with spin flip σ ↔ −σ in the inner region
An interesting scenario is the spin-flip scattering while maintaining the same energy. Instead of the reflected waves with different n's as in Section VI B 2, the reflected waves in the inner region of the wormhole could get spin-flipped σ → −σ with (α, β) → −(α, β) and still have the same energy and momentum. From Eqn. (4.1), this corresponds to n = n, σ = −σ. The scattering coefficients for m = 0, 1 are given in Table V,VI respectively. Note the exact values of the scattering coefficients in m = 0 case to the values in Table III. Wave function profile of the matching is shown in Fig. 6(d). The m -dependence of the transmission and reflection coefficients are shown in Fig. 5.   Another possibility for scattering that results in the outgoing waves in the lower plane is to replace H (1) (kR(u)) with H (2) (−kR(u)) in the lower plane region of Eqn. (4.2) and (4.3). All of the scattering coefficients numerically solved turn out to be the same except the change in the transmission coefficients T → −iT for each corresponding case. The wave function profile in the three regions are also identical in every case.

VII. DISCUSSIONS AND CONCLUSIONS
The scattering of fermion to the 2-dimensional wormhole is calculated by considering both normal and quasinormal modes of the fermionic states in the wormhole. For incoming waves with real momentum and energy, we found the quantum selection of the states allowed to propagate through the wormhole, i.e., the momentum-spin relation given by Eqn. (6.1). Only the states with the quantized momentum related to its angular momentum can be transmitted through the hole via the normal modes and thus be maintained stationary for a long period of time. The waves will also be partially reflected with the same momentum and energy. The unitarity condition (6.2) is valid in this case. Notably, the main result, Eqn. (6.1), is also valid for uncharged or zero-flux fermion scattering with a replacement m → m.
For scattering involving QNMs of the wormhole, the process will be decaying or growing in time due to the imaginary part of the frequency. For QNMs with zero imaginary part, e.g. when m = 0 (this is not exactly a normal mode since the momentum is not real), the energy can be real for massive fermion and sufficiently small n. Even in such cases, the momentum will be pure imaginary resulting in spatial attenuation and enhancement of the wave functions in the upper and lower plane regions. The scattering with QNMs of the wormhole violates unitarity but satisfies more generic relation (5.8). Only the incoming waves with the right complex energy and momentum E = E n,m , k = k n,m will undergo a resonant scattering with the hole resulting in partially transmitted and reflected waves. The physical energy of the particle, however, is given by the real part of E. Remarkably, there will be only transmitted (tunneling) waves with no reflection for m = 0 fermion scattering, i.e., when the magnetic flux is quantized in integer multiples of the magnetic flux quantum φ = mhc/e. For general nonzero m scenario, the transmission and reflection coefficients have curious peaks at m = ±6 Lastly following the argument in Ref. [18], a few comments on fermion scattering in the graphene wormhole can be made. With c → v F 10 6 m/s and M = 0 replacement, the normal-modes energy of the fermion in the graphene wormhole is E = v F m a √ q = 0.658 nm a m √ q eV. Only quasi-fermion with this quantized energy can be transmitted through the graphene wormhole as a stationary state, and it will be accompanied by the reflected waves obeying the unitarity condition. The effects of stitching the graphene wormhole to the plane include the exchange of inequivalent Dirac points that can be taken into account by the effective flux of gauge field [20,26]. This effective flux can simply be added to the magnetic flux in our work resulting in the change of the effective angular momentum number m = m − φ total φ 0 and otherwise the same results.