The soliton nature of the super-Klein tunneling effect

We establish a relationship between the Davey–Stewartson II (DS II) integrable system in $(2{+}1)$ dimensions and quasi-exactly solvable planar interacting Dirac Hamiltonians that exhibit the super-Klein tunneling (SKT) effect. The Dirac interactions are constructed from the real and imaginary parts of breather solutions of the DS II system. In this framework, the SKT effect arises when the energy is tuned to match the constant background of the soliton, while the resulting Dirac Hamiltonians simultaneously support bound states embedded in the continuum. By imposing the SKT boundary conditions, we employ Darboux transformations to construct a general three-parameter family of DS II breather solutions that can be mapped to Dirac Hamiltonians. At the initial soliton time, the corresponding Dirac systems form a massless two-parameter family of Hermitian models with nontrivial electrostatic potentials. As the soliton time evolves, the systems become $\mathcal{PT}$-symmetric and develop a nontrivial imaginary mass term. Finally, when the soliton time is taken to be imaginary, the construction yields Hermitian Dirac systems that lack time-reversal symmetry. In all cases, we identify the emergence of quasi-symmetry transformations that preserve the SKT subspace of states while not commuting with the full Hamiltonian.

💡 Research Summary

The paper establishes a deep connection between the (2+1)-dimensional integrable Davey‑Stewartson II (DS II) system and planar Dirac Hamiltonians that display the super‑Klein tunneling (SKT) effect. By exploiting the Lax‑pair formulation of DS II, the authors identify a Dirac operator hidden in the spatial part of the Lax pair. The complex field u of DS II provides, through its real part, an electrostatic potential V(x,y,τ), while its imaginary part yields a PT‑symmetric mass term m(x,y,τ)σ₃.

A key technical step is the use of a first‑order Darboux transformation. Starting from a constant background solution (U₀ = −iλσ₂, W₀ = 0), a seed matrix Φ built from two independent spinor solutions is introduced. The intertwining operator L = ∂ₓ − Σ, with Σ = (∂ₓΦ)Φ⁻¹, generates new matrix potentials U₁ and W₁. By enforcing invariance under the discrete T‑symmetry (complex conjugation combined with σ₂), Σ is made T‑invariant, guaranteeing that the transformed potentials still solve DS II.

The resulting potential takes the compact form
U_γ(x,y,τ) = −iσ₂