We show that Bose-Einstein condensates in a honeycomb optical lattice are described by a nonlinear Dirac equation in the long wavelength, mean field limit. Unlike nonlinear Dirac equations posited by particle theorists, which are designed to preserve the principle of relativity, i.e., Poincaré covariance, the nonlinear Dirac equation for Bose-Einstein condensates breaks this symmetry. We present a rigorous derivation of the nonlinear Dirac equation from first principles. We provide a thorough discussion of all symmetries broken and maintained.

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 87Rb 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.

We present relativistic linear stability equations (RLSE) for quasi-relativistic cold atoms in a honeycomb optical lattice. These equations are derived from first principles and provide a method for computing stabilities of arbitrary localized solutions of the nonlinear Dirac equation (NLDE), a relativistic generalization of the nonlinear Schrödinger equation. We present a variety of such localized solutions: skyrmions, solitons, vortices, and half-quantum vortices, and study their stabilities via the RLSE. When applied to a uniform background, our calculations reveal an experimentally observable effect in the form of Cherenkov radiation. Remarkably, the Berry phase from the bipartite structure of the honeycomb lattice induces a boson-fermion transmutation in the quasi-particle operator statistics.

We present a thorough analysis of soliton solutions to the quasi-one-dimensional (quasi-1D) nonlinear Dirac equation (NLDE) for a Bose-Einstein condensate in a honeycomb lattice with armchair geometry. Our NLDE corresponds to a quasi-1D reduction of the honeycomb lattice along the zigzag direction, in direct analogy to graphene nanoribbons. Excitations in the remaining large direction of the lattice correspond to the linear subbands in the armchair nanoribbon spectrum. Analytical as well as numerical soliton Dirac spinor solutions are obtained. We analyze the solution space of the quasi-1D NLDE by finding fixed points, delineating the various regions in solution space, and through an invariance relation which we obtain as a first integral of the NLDE. We obtain spatially oscillating multi-soliton solutions as well as asymptotically flat single soliton solutions using five different methods: by direct integration; an invariance relation; parametric transformation; a series expansion; and by numerical shooting. By tuning the ratio of the chemical potential to the nonlinearity for a fixed value of the energy-momentum tensor, we can obtain both bright and dark solitons over a nonzero density background.

We analyze the vortex solution space of the 2 1 + ( )-dimensional nonlinear Dirac equation for bosons in a honeycomb optical lattice at length scales much larger than the lattice spacing. Dirac point relativistic covariance combined with s-wave scattering for bosons leads to a large number of vortex solutions characterized by different functional forms for the internal spin and overall phase of the order parameter. We present a detailed derivation of these solutions which include skyrmions, halfquantum vortices, Mermin-Ho and Anderson-Toulouse vortices for vortex windingwe obtain topological as well as non-topological solutions defined by the asymptotic radial dependence. For arbitrary values of ℓ the non-topological solutions include bright ring-vortices which explicitly demonstrate the confining effects of the Dirac operator. We arrive at solutions through an asymptotic Bessel series, algebraic closed-forms, and using standard numerical shooting methods. By including a harmonic potential to simulate a finite trap we compute the discrete spectra associated with radially quantized modes. We demonstrate the continuous spectral mapping between the vortex and free particle limits for all of our solutions.