Photon Scattering by an Electric Field in Noncommutative Spacetime

As is known, the existence of a small noncommutativity between coordinates would generate nonlocal self-interactions in the electromagnetic theory. To explore some consequences of this effect on the propagation of photons we consider Moyal space half-filled with a static and homogeneous electric field and analyze electromagnetic fluctuations on top of this step-like background. Both the localization of photons and the possibility of photon production by strong electric fields are addressed. Several aspects of the Klein paradox in this setup are discussed as well.


Introduction
As a rule, photons do not interact among themselves.Nevertheless, such selfinteraction is possible in noncommutative (NC) spacetimes.Moyal spacetime is a most natural NC scenario where noncommutativity between coordinates can be represented by the Moyal ⋆-product between fields [1,2].Electrodynamics in this spacetime, known as U (1) ⋆ , has been extensively studied: as a first remarkable feature one finds that noncommutativity introduces nonlinearities in Maxwell's equations through a self-interaction which is nonlocal already at the classical level.In other words, since gauge fields-as functions of NC coordinates-do not commute, U (1) ⋆ holds a close resemblance with nonabelian Yang-Mills theories (see the recent lecture notes [3]).
The goal of the present work is to expose some aspects of this self-interaction in the scattering of photons by a static electric field background.One of our original motivations was the question of whether noncommutativity allows photon creation by intense electromagnetic fields (the effects of noncommutativity on electron/positron production by strong electric fields can be found in [4]).Although the answer requires a quantum field theoretic approach, the full set of stationary solutions contained in this article lays down the ingredients for a subsequent analysis of Schwinger effect in this framework (either through the S-matrix formalism [5,6] or with heat-kernel techniques [7]).
In spite of the similarity between U (1) ⋆ and commutative Yang-Mills theories, constant chromoelectric fields do create gluons [8] whereas homogeneous fields do not produce photons in NC space (some aspects of photon propagation on homogeneous electromagnetic fields have been analyzed in [9] and [10]).For this reason, we study the propagation of photons across the interface of two regions, one with and the other without a background electric field.For simplicity, we take an electrostatic potential that only varies along a fixed spatial direction, and such that an homogeneous electric field is suddenly built up: NC space gets then split into halves by a flat interface.Next, we study the propagation of fluctuations of this background.
More specifically, our setting is the following.We start with (d+1)-dimensional Moyal spacetime, whose coordinates satisfy [x µ , xν ] = 2iθ µν , where θ µν are elements of a real antisymmetric matrix.For simplicity, we take the time coordinate x 0 as an ordinary commuting parameter, i.e., θ µ0 = 0. We then introduce a particular background, namely, a static electric field which vanishes on the left half-space x1 < 0, takes a constant value on the right half-space x1 > 0, and is perpendicular to the interface x1 = 0.The potential A 0 = −E(|x 1 | + x1 )/2 generates such a step-like field; the constant E gives the electric field.
Due to the self-interaction introduced by noncommutativity, photons are scattered by this background.However, as will be shown, in the Moyal representation the interface effectively acts as a step of finite width with (depending on the polarization) delta-functions at the boundaries.This step scatters oblique incident photons; perpendicular beams do not experience any interaction with the electric field.According to the intensity of the background field, the energy and the direction of the incident photons one recognizes different regimes.Apart from the usual scattering-where incident photons split into partially transmitted and reflected light beams-one also finds, for a background which is not particularly strong, the propagation of edge states (viz.localized along the interface).Bound states as perturbations of solitonic solutions have already been found in [11].
On the other hand, for extremely intense electric fields (roughly, |E| ∼ 1/|θ µν |) a realization of the Klein paradox arises.Our result for the scattering of vector particles is similar to the usual (i.e., commutative) Klein paradox for scalar particles [12]: the reflection coefficient R for a wave packet is greater than one.This would indicate that at the boundary of regions with strong electric fields photon beams could be produced.
There is an interesting result of our calculations that concerns two different manifestations of the Klein paradox.As is well-known, a reflection coefficient R greater than one is found either for monochromatic solutions of Dirac equation or for localized solutions of Klein-Gordon equation.In both cases this is considered as a sign of particle creation.At the same time, one finds (both for spinor and scalar particles) the paradoxical result that the transmission coefficient T does not vanish for infinitely strong barriers [13].In the problem studied in this article we find R > 1 (for a localized packet) but, on the other hand, T → 0 as the background field grows to infinity.This shows that these two aspects are not necessarily connected.It also suggests the need to perform a second quantized approach to reliably determine the rate of photon production.
Another aspect which is usually considered in the Klein regime is the divergence of the transmission coefficient for some value of the incident energy; this is usually called super-Klein tunneling (see e.g.[14]).For a bosonic field colliding against an abrupt step this occurs at incident energy equal to half the step height.In our setting there is no super-Klein tunneling, i.e., the transmission coefficient remains finite.The reason lies in the smearing introduced by noncommutativity which turns the abrupt edge into a region of finite slope (cf.[12]).
All these features of photon propagation are analyzed in the present article both for TE as well as for TM modes; sometimes called s-polarized and p-polarized waves, respectively.The former are not affected by gauge transformations; for the latter we choose the temporal gauge.Interestingly, this gauge choice introduces a singular background in the equation of motion.We show that this is a general picture in the Klein regime-independently of the background profile.However, the origin of the singularity can be clearly identified and its consequences removed.
Since its discovery [15], the Klein paradox has been thoroughly analyzed from different perspectives.Many aspects of its relation with particle creation in quantum field theory have been thoroughly discussed-we mention, e.g., [16,17,18,19,20,21].More recently, techniques based on the space-time resolved solutions have provided further insight into the behavior of such systems (see, e.g., [22,23,24,25]).We think that the example considered in this article gives interesting information on the Klein paradox in a different setup, namely, one with propagating photons under nonlocal interactions.Moreover, for a specific photon polarization, we are lead to the problem of a single scalar field colliding against a step of finite width1 .We consequently analyze some aspects of the solutions which, to the best of our knowledge, have not been studied before (the explicit solutions were already written in [27]).
Our article is organized as follows.In Section 2 we briefly review the construction of U (1) ⋆ and focus on the dynamics of fluctuations around a classical background.In Section 3 we introduce our specific background and perform the decoupling between TE and TM modes.In Section 4 we identify the different regimes according to the energy of the incident beam.Sections 5 and 6 are devoted to the study of TE and TM modes, respectively.Singularities arising in the Klein region are resolved in Section 7. In Section 8 we analyze the existence of edge states.Finally, in Section 9, we draw our conclusions.
Coordinate operators xµ generate a NC algebra of classical fields ϕ(x) that can be represented in the space of ordinary functions ϕ(x) on R d+1 through the Weyl-Wigner transform, This establishes a one-to-one map between operators ϕ(x) and functions ϕ(x).We use xµ to denote NC coordinates and unhatted x µ for the ordinary commuting coordinates of Moyal spacetime.One can check that under this correspondence the composition of NC fields is mapped into the Groenewold-Moyal product of functions, This representation can be readily checked, for example, on the original commutator [x µ , x ν ] = x µ ⋆ x ν − x ν ⋆ x µ = 2iθ µν .It can also be shown that any cyclic trace in the algebra of operators ϕ(x) gives, up to a normalization factor, the integral on R d+1 of the corresponding functions ϕ(x) [2].Therefore, any action functional S[ϕ(x)] can be written as the integral of a Lagrangian density obtained by replacing every composition of operators with the ⋆-product of ordinary commuting fields.As a consequence, the NC generalization of local field theories turns out to be nonlocal because interactions become, roughly speaking, smeared over distances of order |θ µν |.As we will see, in the case of electrodynamics invariance under nonlocal gauge transformations leads to a self-interacting theory.
As in the commutative case, the nonlocal symmetry U (1) ⋆ is implemented through the covariant derivative that introduces a gauge field A µ (x).Throughout this article we will work in the adjoint representation so the action of D µ should be understood as where, of course, the commutator is defined with the ⋆-product.The commutator [D µ , D ν ] = i[F µν , •] defines the covariant field strength which clearly shows the resemblance between electromagnetism in Moyal spacetime and ordinary Yang-Mills theories.As briefly explained above, the generalization of Maxwell action to Moyal spacetime is simply given by Upon integration by parts the ⋆-product can be removed from quadratic terms in the action but it must be retained in cubic and quartic self-interactions.To introduce in our setting a background electric field we have included a fixed external source J µ (x).The equation of motion determines a configuration of the gauge field for each external source.Still, an immediate implication of (2.7) is that the current must be covariantly constant, Let us now assume a specific solution A µ (x) and introduce small fluctuations a µ (x) around this background.If we write the full gauge field as A µ (x) + a µ (x), the (linearized) field tensor becomes (2.9) In this expression and hereafter, F µν and D µ are computed exclusively in terms of the background field A µ .Replacing (2.9) into (2.7)we obtain the linearized equation of motion for the perturbations, Note also that, since J µ is the same for A µ and A µ + a µ , condition (2.8) now implies Before considering in detail a specific background, let us examine gauge invariance on the perturbations.In the presence of a background field, gauge invariance remains as a small U (1) ⋆ symmetry on the fluctuations, that preserves their equation of motion as long as Therefore, we can conveniently gauge the perturbations with a parameter α(x) that commutes with the external source.

The background
Let us define our specific setup.As already mentioned, we are interested on the effects of spatial noncommutativity, so we take the time coordinate x 0 as a commuting variable.As for the space coordinates, note that after an appropriate choice the antisymmetric matrix θ µν can be diagonalized in 2 × 2 blocks.In this way, space becomes decomposed into several Moyal planes such that coordinates from different planes commute with each other.For our purposes it is enough to distinguish the coordinates x 1 , x 2 of one of these planes from the remaining d − 2 spatial coordinates.Noncommutativity is then characterized by a positive parameter θ through together with a (d − 2) × (d − 2) antisymmetric matrix in the rest of Moyal space.
As will become evident soon, if the background does not depend on these d − 2 coordinates then this antisymmetric matrix is irrelevant.We now introduce a background A 0 = A 0 (x 1 ) that varies only along the direction of x 1 , and set the spatial components to zero: A 1 = . . .= A d = 0. We will shortly choose a specific x 1 -dependence for A 0 .This static background represents a purely electric field E(x 1 ) = −A ′ 0 (x 1 ) in the x 1 -direction, and originates on a static distribution of charge density J 0 = −A ′′ 0 (x 1 ); spatial components of the current vanish as well, J 1 = . . .= J d = 0.
Since the static background does not change with the transverse coordinates x ⊥ = (x 2 , . . ., x d ), we analyze fluctuations that propagate with fixed energy k 0 and transverse momentum k ⊥ = (k 2 , . . ., k d ), It is now important to point out that the commutator between any background which only depends on x 1 and any field with definite momentum k 2 can be written as an ordinary multiplication between the same field and another function which, roughly, evaluates the variation of the original function at distances ∆x 1 ∼ θk 2 .For example, with This expression shows that if k ⊥ has no component in the direction of x 2 the background is transparent to photon beams.As already mentioned, since the background only depends on x 1 , the propagation of photons will only be affected by θ, and the remaining NC parameters will have no physical consequence.Before solving (2.10) for our specific background, let us establish some convenient notation.We decompose the fluctuations a µ = (a 0 , . . ., a d ) = (a 0 , a 1 ; a ⊥ ) as follows.In an asymptotic region-where the background eventually vanishes-the fluctuations propagate with momentum ⃗ k = (k 1 , ..., k d ) = (k 1 ; k ⊥ ); as expected, in the absence of a background the wave equation (2.10) indicates that the spatial part of a µ is perpendicular to ⃗ k.We separate from the transverse field a ⊥ two components: (i) the projection a P in the plane of incidence (i.e., the component parallel to k ⊥ ), and (ii) a T , which is a (d − 3)-component vector perpendicular to the direction of x 1 and to k ⊥ .This decomposition is illustrated in figure 1 for an incoming wave entering a region of non-vanishing background (represented by the flat surface).With some abuse of notation, In the asymptotic region the electric field of an incoming wave is proportional to the spatial components of the gauge field.Therefore, a 1 , a P describe a beam with electric field in the plane of incidence (TM mode or p-polarized wave).On the other hand, the components a T give an electric field perpendicular to the plane of incidence (TE mode or s-polarized wave).
Through the equation of motion (2.10) each component of a T gets completely decoupled, TM modes then satisfy Klein-Gordon equation for a charged field of mass |k ⊥ | in interaction with an electrostatic potential δA 0 (x 1 ).After choosing a specific background δA 0 we will solve this equation in section 5.The remaining components a 0 (x 1 ), a 1 (x 1 ), a P (x 1 ) satisfy the coupled equations These equations imply condition (2.11), which can now be written as δJ 0 a 0 = 0 or, equivalently, as δA ′′ 0 a 0 = 0.Under this condition any two equations among (3.7)-(3.9)imply the third one.Consequently, we will choose the temporal gauge a 0 = 0, which trivially satisfies (2.11), and will use (3.7) and (3.8) to find a 1 and a P .This will be done in section 6.Now it is time to choose a specific profile for the background A 0 (x 1 ): We take a slope that begins abruptly at x 1 = 0 (as in figure 2), The parameter E ∈ R gives the intensity of the electric field.From the viewpoint of the NC spacetime, this potential represents the boundary of a region supporting a static homogeneous electric field.This setting can be considered as a first approach to the study of the dynamics of fluctuations at the interface where an electric field is built up.Recall that, according to (3.3), the action of A 0 (x 1 ) on oscillations (3.2) turns the slope-like potential into a broadened step of small width 2L (figure 3), with  The width of the step depends on the transverse momentum k 2 ; its slope, on the electric field E.
In fact, inserting (3.10) into (3.4)one obtains In this expression and in the rest of this article we assume sign(k 2 ) < 0; if the component k 2 is positive one simply changes the sign of E. Correspondingly, the adjoint action of the electric field is given by This shows that photons effectively interact with the electric field only in a narrow region of width 2L.
Correspondingly, the homogeneous background on the half-plane is generated by a charge density J 0 = 1 2 E δ(x 1 ) together with a second parallel plate of opposite charge density and located at x 1 = +∞.Condition δJ 0 a 0 = 0 now forces a 0 to vanish at the edges of the step.

Energy regimes
Fluctuations (3.2) are characterized by their energy k 0 and their transverse momentum k ⊥ .As already mentioned, photons incident normally to the interface-or even with no k 2 component-do not interact with the background.Let us call β the angle between the plane of incidence and the x 1 x 2 -plane.One can distinguish, both for TE and TM modes, different energy regimes according to the intensity of the electric field: (i) very strong electric fields such that θ|E| cos β > 1, and (ii) ordinary electric fields for which θ|E| cos β < 1.  Case (i) is represented in figures 4 and 5 for E < 0 and E > 0, respectively.The transverse momentum k ⊥ introduces a mass gap.Accordingly, five different energy regimes can be identified: the blue lines represent photons with positive or negative energy that collide against the step-like background and generate reflected and transmitted waves; this we call the propagation region.Orange dashed lines represent photons which propagate only at one side of the interface with an energy that lies within the mass gap at the other side.Therefore, the photon does propagate across the interface and suffers a total reflection.Finally, there exists a finite interval of energies, usually called Klein zone, for which the incident photon lies above and below the surface of the Dirac sea, depending on the side of the step.This is represented by the red horizontal lines in the figures.In this regime particles (with positive kinetic energy) at one side of the step are seen as antiparticles (with negative kinetic energy) at the other side.
In case (ii) one also finds the usual propagation region but the electric field is not strong enough to separate the mass gaps and generate a Klein zone.As a consequence, there is an overlapping between the mass gaps at either side of the interface so that one could eventually find states localized close to this surface (edge states).This possibility will be analyzed in section 8.

TE modes
We begin by solving equation of motion for TE modes a T .Recall that each of the d − 2 components-which we now generically denote as ϕ(x 1 )-satisfies the Klein-Gordon equation for a charged scalar field of mass |k ⊥ | scattered by an electrostatic potential δA 0 .Taking into account the piecewise definition of δA 0 (x 1 ), we solve (3.6) by appropriately matching at x 1 = ±L the solutions in the left region x 1 < −L, the right region x 1 > L, and the intermediate region |x 1 | < L. These we denote by ϕ L , ϕ R and ϕ I , respectively.
In the regions |x 1 | > L the background δA 0 is constant so the solutions are simply incoming and outgoing plane waves at either side of the step, with (5.4) In the intermediate region |x 1 | < L the function δA 0 is linear in x 1 so the field satisfies where we have shifted the x 1 coordinate as The solutions can be written as in terms of parabolic functions where ν = −1/2 + ik2 ⊥ /2|E|.These particular combinations where chosen to satisfy After imposing continuity of ϕ and its derivative at x 1 = ±L one obtains four equations: two of them give the coefficients A, B in the intermediate region in terms of a, b, c, d; the remaining two equations determine the elements of the S-matrix, that relates the outgoing coefficients b, c with the incoming coefficients a, d through A straightforward calculation gives The upper (lower) sign corresponds to r L (r R ).We have also defined k L = k 1 and k R = κ 1 .We use F L , G L to denote F (L), G(L), and similarly for the derivatives.Relation (5.11) parametrizes the two-dimensional space of solutions.The coefficients a and d represent incoming waves from the left and from the right, respectively.Thus, if d = 0 then b is related to the reflection coefficient and c to the transmission coefficient of a wave colliding the step from the left.Correspondingly, a = 0 describes a wave colliding from the right and the roles of b and c get interchanged.To define reflection and transmission coefficients we compute the conserved current 2 which according to (3.6) must be x 1 -independent.For x 1 < −L, expression (5.1) gives The same expression holds for x 1 > L if one makes the following replacements:  For energies in the propagation and in the Klein regimes both k 1 and κ 1 are positive real numbers, so current conservation J(−∞) = J(+∞) implies ( This expression-valid for arbitrary values of the constants a, d-proves that S is a unitary matrix.Further relations between the elements of the S-matrix can be obtained by noting that the complex conjugates of (5.1) and (5.Although these relations can be readily checked in (5.12) and (5.13), their derivation did not make use of the specific form of the background.In fact, this is a well-known result from one-dimensional scattering: reflection (and transmission) coefficients for left-or right-incident waves coincide.Reflection and transmission coefficients are then defined as 6 displays R as a function of the incident energy k 0 .The reflection coefficient decreases monotonically from 1 to 0 with increasing |k 0 |.If the energy lies within one of the mass gaps then either k 1 or κ 1 are pure imaginary so J = 0 and the incoming photon is totally reflected at the interface (R = 1).
In the Klein zone R is smaller than 1 but to get an appropriate interpretation of the reflection coefficient it is crucial to recall that for these energies phase and group velocities have opposite signs.Let us consider, for example, the case E < 0 (figure 4).For incident energies |k ⊥ | < k 0 < −2EL − |k ⊥ | both momenta k 1 and κ 1 are real but the group velocity at (5.17This is the Klein paradox which is interpreted as the result of a total reflection of the incident packet together with the creation of pairs of photons at the interface at a rate T .Alternatively, the term "Klein paradox" is sometimes used to refer to the existence of a non-vanishing transmission coefficient even for a infinitely high barrier-this is the case for a Klein-Gordon or Dirac particle colliding an ordinary abrupt step.Interestingly, this does not happen in our setting.In fact, the leading behavior of (5.8) and (5.9) for large |E| reads [28] Consequently, the behavior of T for large values of the electric field is given by (5.20) The dependence of T on the electric field E is shown in figure 7 for a wave packet.
In the propagation region 0 < T < 1.If the energy of the packet falls in the mass gap then T = 0.For strong backgrounds, T may be even greater than 1-the transmitted pulse is interpreted as photons created at the step.Nevertheless, T → 0 as |E| → ∞.
Finally, let us recall that for an ordinary scalar field pulse colliding against an abrupt step (i.e., a step of zero width) the coefficients R and T diverge if the energy is half the step height-this is usually known as super-Klein tunneling.However, the reflection coefficient shown in figure 6 (for a monochromatic wave) does never vanishes, so the corresponding coefficient for the wave packet (obtained by replacing R → 1/R) does not diverge.This is a consequence of the smearing of the step in Moyal space.In the limit case θ → 0 and |E| → ∞ (with |E|θ finite) one obtains and abrupt step and the reflection coefficient for a monochromatic wave certainly vanishes if its energy is exactly half the step height and one recovers super-Klein tunneling in this specific limit.

TM modes
Components a T are unaffected by a gauge transformation (2.12) but to determine a 0 , a 1 , a P we need to fix a specific gauge.We solve (3.7)-(3.9)by setting a 0 (x 1 ) = 0 with an appropriate parameter of the form α(x) = α(x 1 ) e −ik 0 x 0 +ik ⊥ x ⊥ .Condition (2.13) requires α(x 1 ) to vanish at the edges of the step x 1 = ±L but, since any solution a 0 (x 1 ) vanishes at those points, one can consistently choose α such that a 0 + D 0 α = 0.
For a 0 = 0 equations (3.7) and (3.8) read As already mentioned, (3.9) can be derived from these two equations.From (6.1) we eliminate the component which is parallel to the transverse momentum k ⊥ , The other component can be written as where φ(x 1 ) satisfies the Schrödinger-like equation Up to now we have not used the explicit dependence of the background on x 1 .For the step-like potential (3.12) the fourth term inside the square brackets only introduces delta-functions at the edges of the step, We now solve (6.5) in the same fashion as for the TE modes: we find a piecewise solution and match the coefficients by imposing the appropriate behavior at the edges of the step.At x 1 < −L and x 1 > L the field φ (respectively denoted by φ L and φ R ) takes the form with k 1 and κ 1 once more given by (5.3) and (5.4).At the interval |x 1 | < L the field (denoted by φ I ) satisfies where the variable z is the same as (5.6).It is interesting to make a comparison with the equation for the TM modes (5.5): the dynamics of TE modes contains an additional, singular term 1/z 2 .The effect of this term will be discussed in some detail.Solutions to (6.9) can be written as combinations of the following confluent hypergeometric functions, , − 1 2 ; i|E|z 2 ) , (6.10) The field at the intermediate region is thus determined by two coefficients A, B, where now These particular solutions have been chosen such that F (−L) = Ḡ′ (−L) = 0 and Next, we match the behavior of φ L , φ R , φ I at the edges of the step to construct the full solution φ.We demand continuity of the wave function but, in accordance with (6.6), we introduce discontinuities on the derivatives.The matching conditions then read φ R (L) = φ I (L) , (6.17) Figure 8: Gauge field components a 1 (blue) and a P (orange) with energy in the propagation regime.In black, the background field (E < 0).The deltafunctions at x 1 = ±L introduce discontinuities in the derivatives.(blue) and a P (orange) with energy in the Klein regime.The background field (black) corresponds to E < 0. Due to the singularity in (6.9), the solutions diverge at z = 0.
From these four equations we compute the elements rL , rR , tL , tR of the S-matrix which, as in (5.11), relate outgoing with incoming coefficients, As before, the upper (lower) sign corresponds to r L (r R ); we define k L = k 1 and k R = κ 1 and we use FL , ḠL , F ′ L , Ḡ′ L to denote F (L), Ḡ(L), F ′ (L), Ḡ′ (L). Figure 8 shows the full solution for the gauge field components a 1 (x), a P (x 1 ) with energy in the propagation regime.As expected, the wave is continuous but its derivative jumps at the edges of the step.Figure 9 displays instead both components a 1 (x 1 ), a P (x 1 ) in the case where the energy of the gauge field lies in the Klein zone.The most remarkable feature is the singularity at some The obvious reason for this divergence is the singularity in the differential equation (6.9) at z = 0.For energies in the propagation region z never vanishes in the interval x 1 ∈ [−L, L] but in the Klein regime there is some value of x 1 in this interval where z = 0.In consistence with figure 9, the component a 1 has a double pole, whereas a P has a simple pole.We postpone to section 7 a more detailed discussion of the physical origin of this divergence and a procedure to appropriately remove it.The interval within the red vertical dashed lines is the Klein zone.The intervals where R = 1 correspond to the mass gaps.The outer region, delimited by the orange vertical dashed lines is the propagation region.We also display the reflection coefficient for the TE modes (dashed green) as in figure 6.We have chosen k 2 = 0.25, k ⊥ = 0.5, E = −8 in length units of √ θ.
Let us comment on the solutions (6.19) and (6.20).We have obtained t L = t R and, since F , Ḡ are real functions, one readily checks |r L | = |r R |, as long as k 1 , κ 1 ∈ R. The reflection coefficient for a monochromatic wave is then given by R = |r L | 2 .As before, unitarity of the S-matrix implies R + T = 1. Figure 10 shows the reflection coefficient as a function of the photon energy k 0 .The plot is similar to the one corresponding to TE modes but with a sudden decrease in the vicinity of the mass-gap regions, where R = 1.The figure shows that, in general, TE modes suffer a higher reflection than TM modes.

The Klein zone
Let us analyze the divergences that occur in the Klein zone for TM modes.The emergence of singularities stems from the covariant derivative D 0 = ∂ 0 + iA 0 (x 1 ) which acting on stationary solutions introduces the factor k 0 − δA 0 (x 1 ).In a typical Klein paradox scenario the potential attains asymptotically constant values δA 0 (±∞); for intense backgrounds their difference δA 0 (+∞) − δA 0 (−∞) exceeds the mass gap and thus the band of energies known as Klein zone is created.For k 0 within this band there is necessarily at least one point x 1 (assuming continuous backgrounds) such that k 0 = δA 0 (x 1 ), so the covariant derivative vanishes at this point.Decoupling the components of the gauge field in (3.7)-(3.9)then leads to the appearance of a singular term (k 0 − δA 0 ) −2 in the equation of motion (third term in (6.5)).
To be more concrete, let us turn back to our background (3.12), which behaves as As a consequence, the field φ(x 1 ) scatters against a singular background of the type z −2 .Close to a singularity of this type the solutions show two different behaviors: z −1 and z 2 (as shown by (6.10) and (6.11), respectively).The regular solution (6.11) leads to finite gauge field components a 1 (x 1 ), a P (x 1 ) but How do we interpret the singular gauge field that arises from the solution (6.10)?The answer is that the singular behavior comes from a singular gauge transformation.Consider a (nonsingular) solution of the set of equations (3.7)-(3.9)before any specific gauge choice.In general a 0 (x 1 ) ̸ = 0 but we can impose the temporal gauge by choosing α(x 1 ) such that D 0 α(x 1 ) = −i(k 0 − δA 0 (x 1 ))α(x 1 ) = −a 0 (x 1 ).Therefore, in the Klein zone α(x 1 ) has a singularity at x 1 = −L − k 0 /E.As a consequence, in this gauge the time component vanishes but the spatial components a 1 , a P become singular.In fact, since α has a simple pole at z = 0, the transformation a µ → a µ + D µ α introduces a simple pole in a P and double pole in a 1 .This is precisely what one obtains if one inserts the solution φ ∼ z −1 into (6.3) and (6.4).Fortunately, it is not necessary to abandon this gauge choice: one can undo the singular gauge transformation from (6.12) to obtain a regular gauge field (figure 11) or, alternatively, one can keep the singular solution but choosing the appropriate behavior at z = 0. We choose the latter approach.
In general, one-dimensional scattering against a singularity detaches the solutions at both sides of the singular point unless a certain matching condition is specified.In other words, one could choose a certain value for the coefficients A, B at one side of the singularity and a different value at the other side.If one is interested-to give an example-in solutions of definite parity, then one should choose them such that only one of the solutions (6.10) or (6.11) remains.One could even choose A = B = 0 at one side of the singularity if it fits the physical setting.In our case, we must keep in mind that after removing the simple pole from the solution φ one must get regular fields a 0 , a 1 , a P .Therefore, the values of both coefficients A, B must not change after crossing z = 0.This solves the ambiguity that arises due to the singularity and permits to match the solutions (6.7) and (6.8) to determine the relation between the constants a, b, c, d.For this reason, the result expressed in (6.19) and (6.20) still holds in the Klein zone.
Note that, as for the TE modes, R < 1 even in the Klein zone.Of course, for a Figure 11: Gauge field components a 0 (magenta), a 1 (blue) and a P (orange) after removing the singular gauge transformation for the same parameters as in figure 9. To display all components in the same figure, a 0 (magenta) has been amplified and a 1 (blue) reduced by a factor 10. of TM modes (blue).At the left of the orange dashed line, the propagation zone; at the right of the red dashed line, the Klein zone.We also display T for TE modes (dashed green) as in figure 7. We have chosen E < 0, k 2 = 0.25, k ⊥ = 0.5, k 0 = 1 in length units of √ θ.
monochromatic wave with positive group velocity at x 1 > L, one takes c = 0 and d = 1 in (6.7)-(6.8), the reflection coefficient results 1/|r L | 2 and is thus greater than 1.As for TE modes, this is interpreted as the result of photon production: the wave is totally reflected but a flux of photons created at the step runs in both directions.
Finally, one can also evaluate for TM modes the Klein paradox in the sense of a non-vanishing transmitted wave even for an infinitely high step.As for the TE modes, this does not occur for TM modes either.Using the standard asymptotic behavior of confluent hypergeometric functions [28] one obtains the transmission coefficient for a wave packet for large values of the background electric field3 , Figure 12 shows T as a function of E.

Edge states
For very intense electric fields such that θ|E| < |k ⊥ |/|k 2 | a Klein zone is created between the left-and the right-mass gaps.Propagation with energies within this zone has been discussed in the previous sections.
On the contrary, for θ|E| < |k ⊥ |/|k 2 | there is an overlapping between both mass gaps.If the energy lies within this overlapping then both k 1 and κ 1 are purely imaginary, so it could be possible to find waves localized in the vicinity of the interface x 1 = 0 (edge states).Of course, for this to happen the solution must be appropriately matched all along the x 1 axis.
Propagation of TE modes is described by the Klein-Gordon equation with a step-like electrostatic potential that generates an homogeneous electric field in a region of width 2L, so states localized in the x 1 -direction are not to be expected.
As regards the TM modes, inspection of (6.5) shows that, apart from the same step-like potential, the field is subject to extra terms which depend on the incident energy k 0 .One of these terms includes a delta-function with a negative coefficient so the equation could admit, a priori, bounded solutions.
To obtain bound states in the x 1 -direction we look for solutions (6.7) and (6.8) with k 1 and κ 1 in the negative imaginary semi-axis, such that b = c = 0 to match the appropriate behavior at x 1 → ±∞.This condition can only be implemented if the S-matrix is singular.In conclusion, the equation det S = 0 determines the energies at which edge states could occur.
Figure 13 shows that for some values of the parameters E, θ, k 2 , k ⊥ there is a finite number of edge states.

Conclusions
Electromagnetism in Moyal space is a self-interacting theory, so perturbations of a classical configuration do not propagate freely but interact with the underlying background.Since this interaction does not show up for homogeneous backgrounds, we have considered two half-spaces separated by a flat interface: one half is free from background fields, the other is filled up with an homogeneous electric field.One each side a photon beam propagates free from interactions but suffers scattering at the interface.In this article we have determined the stationary states, both for TE and TM modes.Each gauge field component which is perpendicular to the plane of incidence (TE mode) decouples and propagates as a massive, charged scalar particle.The mass is given by the momentum component parallel to the interface.Since the gauge field belongs to the adjoint representation of the ⋆-product, each TE mode only interacts with a homogeneous electric field restricted to a narrow region around the interface.As a consequence, TE modes reproduce the dynamics of a relativistic scalar field interacting with a step-like electrostatic field: R < 1 for a monochromatic wave (but R > 1 for a wave packet in the Klein regime) and, due to the finite width of the potential, super-Klein tunneling is suppressed and the transmission coefficients vanishes as the step height tends to infinity.
To solve for TM modes one needs to choose a particular gauge.We take the temporal gauge and solve the equations of motion.The reflection coefficient is computed and compared to the result for TE modes (figure 10).Also in this case, R < 1 for a monochromatic wave (as for scalar fields; unlike fermionic fields).In the Klein regime the wave interacts with a singularity, which is an artifact of the gauge choice.A careful analysis allows one to determine the appropriate matching conditions across the singularity.After removing the singularity with a gauge transformation, we determine the profiles of the propagating gauge fields.
As the electric field goes to infinity, the transmission coefficient also vanishes.The asymptotic behavior of T for TE and TM modes obeys the same power law but with a different coefficient ((5.20) and (7.1)).
For moderate electric fields, such that there is no separation between the leftand right-mass gaps, it might be possible that the electric field supports waves localized in the vicinity of the interface.Of course, this could only be the case for TM modes, for which the one-dimensional equation (6.5) contains a delta-function with a negative coefficient.Figure 13 shows that for some values of the electric field a finite number of edge states may appear.
It is important to point out that the conserved current computed in this article-both for TE and TM modes-can be obtained from the asymptotic behavior of the canonical energy-momentum tensor T µν = − 1 2 {F µρ , ∂ ν A ρ } ⋆ + 1 4 η µν F 2 ⋆ (see [29]).As in ordinary Yang-Mills theories the construction of symmetric, traceless, gauge covariant or (covariantly) conserved energy-momentum tensors in NC gauge theories has been thoroughly studied [29,30,31].Nevertheless, the conserved current considered in the present article is only used to illustrate the asymptotic behavior of the gauge fields profiles; our main purpose is to present the classical solutions which can be subsequently used to expand the quantized gauge field.This article therefore provides the spectrum of fluctuations around a specific background.In particular, fluctuation modes with energies in the Klein zone would indicate that inhomogeneities in the electric field are capable of generating photon beams.We think that an expansion of the quantum field in the stationary modes described in this article would provide the rate of photon creation.In the Smatrix formalism, this is usually connected with the integral of the transmission coefficient-as a function of the incident energy-in the Klein region.On the other hand, the spectral decomposition of the operator of quantum fluctuations can also be used to compute the heat-trace, from which the imaginary part of the effective action (and, thus, the stability of the vacuum) can be studied.It would be interesting to compare the results with the already studied heat-kernel methods for NC theories [32,33,34].Work along this two lines of research are currently under consideration.

Figure 1 :
Figure 1: Decomposition of the gauge field into TE modes a T (green) and TM modes a 1 , a P (blue).

Figure 2 :
Figure 2: Background slope from the NC spacetime perspective.In the picture E < 0 so the electric field points towards the left.

Figure 3 :
Figure 3: The background as it interacts with incoming photons (for sign(k 2 )E > 0).The width of the step depends on the transverse momentum k 2 ; its slope, on the electric field E.

Figure 4 :
Figure 4: The black solid line represents the potential δA 0 for a strong electric field E < 0. Horizontal lines correspond to five different energy regimes.

Figure 5 :
Figure 5: The black solid line represents the potential δA 0 for a strong electric field E > 0. Horizontal lines correspond to five different energy regimes.

Figure 7 :
Figure 7: T vs. |E| for a pulse of TE modes.At the left of the orange dashed line lies the propagation zone; at the right of the red dashed line, the Klein zone.In this picture E < 0, k 2 = 0.25, k ⊥ = 0.5, k 0 = 1 in length units of √ θ.
2) provide another stationary solution but where the roles of b, c and a, d are played by a * , d * and b * , c * , respectively.Therefore, S * = S −1 = S † .This implies t R = t L and |r L | = |r R |.
) Therefore, the stationary solutions which are relevant to model the scattering of a left-incident wave packet are not obtained by choosing d = 0 in (5.1) and (5.2) but instead by taking c = 0.The reflection and transmission coefficients are then given by R = 1/|r L | 2 and T = |t L | 2 /|r L | 2 and the relation R = 1 + T > 1 holds.

Figure 12 :
Figure 12: T vs. |E| for a wave packet

Figure 13 :
Figure 13: The figure shows det S as a function of k 0 for different values of the electric field (E < 0).Each curve is defined in the interval −2EL − k P < k 0 < k P , which corresponds to the overlapping of the mass gaps.The transverse momenta were chosen as k 2 = 0.25, k P = 1 (in units of √ θ.)Each intersection with the horizontal axis gives the energy of an edge state.The number of edge states decreases with increasing |E|.