We demonstrate that the confinement of half-quantized vortices (HQVs) in coherently coupled Bose-Einstein condensates simulates certain aspects of the confinement in a theory like
SU(2) quantum chromodynamics (QCD) in 2+1 space-time dimensions. By identifying the circulation of superfluid velocity as the baryon number and the relative phase between two components as a dual gluon, we identify HQVs in a single component as electrically charged particles with a half baryon number. Further, we show that only singlet states of the relative phase of two components can stably exist as bound states of vortices, that is, a pair of vortices in each component (a baryon) and a pair of a vortex and an antivortex in the same component (a meson). We then study the dynamics of a baryon and meson; a baryon is static at equilibrium and rotates once it deviates from the equilibrium, while a meson moves with constant velocity. For both baryons and mesons we verify a linear confinement and determine that they are broken, thus creating other baryons or mesons in the middle when two constituent vortices are separated by more than some critical distance, resembling QCD.
Vortex reconnections and interactions play a fundamental role in the dynamics of fluids and turbulent flows, both in the classical and quantum regime. Studying vortices in a clean system like ultracold gases can therefore help, as a bottom-up approach, to understand the physics in a wider context, including superfluid helium, polariton condensates, fluid dynamics and turbulence, neutron stars and cosmological models. So far vortex-vortex interaction was studied in Bose-Einstein condensates either in rotating systems with the observation of regular Abrikosov lattices or in flat condensates across the Berezinskii-Kosterlitz-Thouless transition. In both cases the geometrical constraints allowed to study just a planar interaction among aligned or anti-aligned vortices. The present study instead is carried out in an axisymmetric cigar-shaped Bose-Einstein condensate. Vortices in prolate structures are also known as solitonic vortices. This geometry is especially suitable for investigating vortex interactions. Indeed vortices are oriented perpendicularly to the condensate axis to minimize their length, hence energy, and, because of the cylindrical symmetry, the orientation of a vortex in the radial plane has no constraints. These facts permit interactions to occur with different incoming relative angles between the vortex lines and with different relative velocities. The strong confinement, acting along the radial axis, enhances also interesting effects due to the boundaries.