Nonlinear Dirac equation in Bose-Einstein condensates: Preparation and stability of relativistic vortices

We propose a detailed experimental procedure for preparing relativistic vortices, governed by the nonlinear Dirac equation, in a two-dimensional Bose-Einstein condensate (BEC) in a honeycomb optical lattice. Our setup contains Dirac points, in direct analogy to graphene. We determine a range of practical values for all relevant physical parameters needed to realize relativistic vortices in a BEC of ⁸⁷Rb atoms. Seven distinct vortex types, including Anderson-Toulouse and Mermin-Ho skyrmion textures and half-quantum vortices, are obtained, and their discrete spectra and stability properties are calculated in a weak harmonic trap. We predict that most vortices are stable, with a lifetime between 1 and 10 s.