We theoretically explore atomic Bose-Einstein condensates (BECs) subject to position-dependent spin-orbit coupling (SOC). This SOC can be produced by cyclically laser coupling four internal atomic ground (or metastable) states in an environment where the detuning from resonance depends on position. The resulting spin-orbit coupled BEC phase-separates into domains, each of which contain density modulations - stripes - aligned either along the x or y direction. In each domain, the stripe orientation is determined by the sign of the local detuning. When these stripes have mismatched spatial periods along domain boundaries, non-trivial topological spin textures form at the interface, including skyrmions-like spin vortices and anti-vortices. In contrast to vortices present in conventional rotating BECs, these spin-vortices are stable topological defects that are not present in the corresponding homogenous stripe-phase spin-orbit coupled BECs.
Equilibrium vortex formation in rotating binary Bose gases with a rotating frequency higher than the harmonic trapping frequency is investigated theoretically. We consider the system being evaporatively cooled to form condensates and a combined numerical scheme is applied to ensure the binary system being in an authentic equilibrium state. To keep the system stable against the large centrifugal force of ultrafast rotation, a quartic trapping potential is added to the existing harmonic part. Using the Thomas-Fermi approximation, a critical rotating frequency Omega_c is derived, which characterizes the structure with or without a central density hole. Vortex structures are studied in detail with rotation frequency both above and below ?Omega_c and with respect to the miscible, symmetrically separated, and asymmetrically separated phases in their nonrotating ground-state counterparts.
Formation of stationary 3D wave patterns generated by a small point-like impurity moving through a Bose-Einstein condensate with supersonic velocity is studied. Asymptotic formulae for a stationary far-field density distribution are obtained. Comparison with three-dimensional numerical simulations demonstrates that these formulae are accurate enough already at distances from the obstacle equal to a few wavelengths.
