Massive Dirac states bound to vortices by a boson-fermion interaction

Abstract: Results are presented from numerical simulations of the flat-space nonlinear Maxwell-Klein-Gordon-Dirac equations. The introduction of a boson-fermion interaction allows a scalar vortex to act as a harmonic trap that can confine massive Dirac bound states. A parametric analysis is performed to understand the range of boson-fermion coupling strengths, Ginzburg-Landau parameters, and fermion effective masses that support the existence of bound state solutions. Solutions are time-evolved and are observed to be stable until the fermion number density becomes large enough to collapse the spontaneously broken vacuum of the condensate. Head-on scattering simulations are performed, and traditional vortex right-angle scattering is shown to break down with increased fermion number density, $n_f$. For sufficiently large $n_f$ and low velocity, the collision of two vortices results in a pseudostable bound state with winding number $m=2$ that eventually becomes unstable and decays into two $m=1$ vortices. For large $n_f$ and collision velocity, vortex scattering is observed to produce nontopological (zero winding number) scalar bound states that are ejected from the collision. The scalar bubbles contain coherent fermion bound states in their interiors and interpolate between the spontaneously broken vacuum of the bulk and the modified vacuum induced by the boson-fermion interaction.