Rabi-Coupled Spinor BECs

• Rabi-coupled spinor Bose–Einstein condensates are quantum fluids where an external Rabi field coherently mixes internal atomic states, locking phases and inducing magnetic transitions.
• They reshape miscibility, collective excitations, and nonlinear dynamics, leading to tunable phase boundaries, soliton formation, and altered equations of state.
• Mean-field models and Bogoliubov analyses reveal that Rabi coupling introduces spectral gaps in the spin channel while stabilizing density modes across diverse trapping configurations.

Rabi-coupled spinor Bose–Einstein condensates are multicomponent quantum fluids in which two or more internal atomic states are coherently mixed by an external field, typically with negligible momentum transfer in the pure Rabi limit. In the pseudo-spin-$1/2$ case, the coupling acts as a transverse field in spin space, locks the relative phase between components, and converts an otherwise gapless spin sector into a gapped one; in lattices, traps, and spin–orbit-coupled settings it also reshapes miscibility, collective excitations, topological defects, and superfluid transport (Recati et al., 2021). The subject spans uniform mixtures, state-dependent optical lattices, harmonically confined condensates, magnetic and bright solitons, and recent nonlinear-regime experiments showing that coherent coupling can even alter the equation of state itself (Eid et al., 29 Jul 2025).

1. Microscopic formulation and mean-field description

For a two-component condensate with hyperfine states $\sigma=1,2$, the standard second-quantized Hamiltonian in the Rabi-coupled case is

$$\hat H = \int d^3 r\, \Bigg\{ \sum_{\sigma=1,2}\hat \psi_\sigma^\dagger\Big(-\frac{\hbar^2\nabla^2}{2m}+V(\mathbf r)\Big)\hat \psi_\sigma + \frac{g_{11}}{2} \hat n_1^2 + \frac{g_{22}}{2} \hat n_2^2 + g_{12} \hat n_1 \hat n_2 - \frac{\hbar \Omega}{2}\big(\hat \psi_1^\dagger \hat \psi_2 + \hat \psi_2^\dagger \hat \psi_1\big) + \frac{\hbar \delta}{2}(\hat n_1 - \hat n_2) \Bigg\},$$

with $g_{\sigma\sigma'}=4\pi \hbar^2 a_{\sigma\sigma'}/m$, Rabi frequency $\Omega$, and detuning $\delta$ (Recati et al., 2021). At mean-field level, the order parameters $\psi_1,\psi_2$ obey coupled Gross–Pitaevskii equations,

$$i\hbar \partial_t \psi_1 = \left(-\frac{\hbar^2\nabla^2}{2m}+V+g_{11} n_1+g_{12} n_2+\frac{\hbar\delta}{2}\right)\psi_1 - \frac{\hbar\Omega}{2}\psi_2,$$

$$i\hbar \partial_t \psi_2 = \left(-\frac{\hbar^2\nabla^2}{2m}+V+g_{22} n_2+g_{12} n_1-\frac{\hbar\delta}{2}\right)\psi_2 - \frac{\hbar\Omega}{2}\psi_1,$$

where $n_\sigma=|\psi_\sigma|^2$ (Recati et al., 2021).

A useful spinor rewriting emphasizes the competition between interaction-induced imbalance and coherent locking. In the dimensionless lattice formulation used for two-component optical-lattice problems, the interaction-plus-Rabi contribution can be written as

$\sigma=1,2$0

with $\sigma=1,2$1, $\sigma=1,2$2, and $\sigma=1,2$3 (He et al., 2021). In this form, the Rabi term explicitly penalizes relative-phase and relative-density fluctuations.

External potentials enter in several distinct ways. Spin-independent lattices impose $\sigma=1,2$4, matter gratings correspond to $\sigma=1,2$5, $\sigma=1,2$6, and Zeeman lattices use opposite-sign state-dependent potentials, $\sigma=1,2$7 and $\sigma=1,2$8 (He et al., 2021). The same coherent-coupling structure extends beyond pseudo-spin-$\sigma=1,2$9: in spin-1 condensates, the Rabi coupling is implemented as a term proportional to $$\hat H = \int d^3 r\, \Bigg\{ \sum_{\sigma=1,2}\hat \psi_\sigma^\dagger\Big(-\frac{\hbar^2\nabla^2}{2m}+V(\mathbf r)\Big)\hat \psi_\sigma + \frac{g_{11}}{2} \hat n_1^2 + \frac{g_{22}}{2} \hat n_2^2 + g_{12} \hat n_1 \hat n_2 - \frac{\hbar \Omega}{2}\big(\hat \psi_1^\dagger \hat \psi_2 + \hat \psi_2^\dagger \hat \psi_1\big) + \frac{\hbar \delta}{2}(\hat n_1 - \hat n_2) \Bigg\},$$0, directly mixing the $$\hat H = \int d^3 r\, \Bigg\{ \sum_{\sigma=1,2}\hat \psi_\sigma^\dagger\Big(-\frac{\hbar^2\nabla^2}{2m}+V(\mathbf r)\Big)\hat \psi_\sigma + \frac{g_{11}}{2} \hat n_1^2 + \frac{g_{22}}{2} \hat n_2^2 + g_{12} \hat n_1 \hat n_2 - \frac{\hbar \Omega}{2}\big(\hat \psi_1^\dagger \hat \psi_2 + \hat \psi_2^\dagger \hat \psi_1\big) + \frac{\hbar \delta}{2}(\hat n_1 - \hat n_2) \Bigg\},$$1 components (Muruganandam et al., 2021).

2. Uniform phases, miscibility, and the equation of state

In the homogeneous pseudo-spin-$$\hat H = \int d^3 r\, \Bigg\{ \sum_{\sigma=1,2}\hat \psi_\sigma^\dagger\Big(-\frac{\hbar^2\nabla^2}{2m}+V(\mathbf r)\Big)\hat \psi_\sigma + \frac{g_{11}}{2} \hat n_1^2 + \frac{g_{22}}{2} \hat n_2^2 + g_{12} \hat n_1 \hat n_2 - \frac{\hbar \Omega}{2}\big(\hat \psi_1^\dagger \hat \psi_2 + \hat \psi_2^\dagger \hat \psi_1\big) + \frac{\hbar \delta}{2}(\hat n_1 - \hat n_2) \Bigg\},$$2 problem with $$\hat H = \int d^3 r\, \Bigg\{ \sum_{\sigma=1,2}\hat \psi_\sigma^\dagger\Big(-\frac{\hbar^2\nabla^2}{2m}+V(\mathbf r)\Big)\hat \psi_\sigma + \frac{g_{11}}{2} \hat n_1^2 + \frac{g_{22}}{2} \hat n_2^2 + g_{12} \hat n_1 \hat n_2 - \frac{\hbar \Omega}{2}\big(\hat \psi_1^\dagger \hat \psi_2 + \hat \psi_2^\dagger \hat \psi_1\big) + \frac{\hbar \delta}{2}(\hat n_1 - \hat n_2) \Bigg\},$$3, it is convenient to introduce

$$\hat H = \int d^3 r\, \Bigg\{ \sum_{\sigma=1,2}\hat \psi_\sigma^\dagger\Big(-\frac{\hbar^2\nabla^2}{2m}+V(\mathbf r)\Big)\hat \psi_\sigma + \frac{g_{11}}{2} \hat n_1^2 + \frac{g_{22}}{2} \hat n_2^2 + g_{12} \hat n_1 \hat n_2 - \frac{\hbar \Omega}{2}\big(\hat \psi_1^\dagger \hat \psi_2 + \hat \psi_2^\dagger \hat \psi_1\big) + \frac{\hbar \delta}{2}(\hat n_1 - \hat n_2) \Bigg\},$$4

At zero detuning, the coherent coupling produces a competition between an unpolarized paramagnetic state and a polarized ferromagnetic state. The paramagnetic solution $$\hat H = \int d^3 r\, \Bigg\{ \sum_{\sigma=1,2}\hat \psi_\sigma^\dagger\Big(-\frac{\hbar^2\nabla^2}{2m}+V(\mathbf r)\Big)\hat \psi_\sigma + \frac{g_{11}}{2} \hat n_1^2 + \frac{g_{22}}{2} \hat n_2^2 + g_{12} \hat n_1 \hat n_2 - \frac{\hbar \Omega}{2}\big(\hat \psi_1^\dagger \hat \psi_2 + \hat \psi_2^\dagger \hat \psi_1\big) + \frac{\hbar \delta}{2}(\hat n_1 - \hat n_2) \Bigg\},$$5 is stable when

$$\hat H = \int d^3 r\, \Bigg\{ \sum_{\sigma=1,2}\hat \psi_\sigma^\dagger\Big(-\frac{\hbar^2\nabla^2}{2m}+V(\mathbf r)\Big)\hat \psi_\sigma + \frac{g_{11}}{2} \hat n_1^2 + \frac{g_{22}}{2} \hat n_2^2 + g_{12} \hat n_1 \hat n_2 - \frac{\hbar \Omega}{2}\big(\hat \psi_1^\dagger \hat \psi_2 + \hat \psi_2^\dagger \hat \psi_1\big) + \frac{\hbar \delta}{2}(\hat n_1 - \hat n_2) \Bigg\},$$6

whereas the ferromagnetic phase occurs for

$$\hat H = \int d^3 r\, \Bigg\{ \sum_{\sigma=1,2}\hat \psi_\sigma^\dagger\Big(-\frac{\hbar^2\nabla^2}{2m}+V(\mathbf r)\Big)\hat \psi_\sigma + \frac{g_{11}}{2} \hat n_1^2 + \frac{g_{22}}{2} \hat n_2^2 + g_{12} \hat n_1 \hat n_2 - \frac{\hbar \Omega}{2}\big(\hat \psi_1^\dagger \hat \psi_2 + \hat \psi_2^\dagger \hat \psi_1\big) + \frac{\hbar \delta}{2}(\hat n_1 - \hat n_2) \Bigg\},$$7

with polarization

$$\hat H = \int d^3 r\, \Bigg\{ \sum_{\sigma=1,2}\hat \psi_\sigma^\dagger\Big(-\frac{\hbar^2\nabla^2}{2m}+V(\mathbf r)\Big)\hat \psi_\sigma + \frac{g_{11}}{2} \hat n_1^2 + \frac{g_{22}}{2} \hat n_2^2 + g_{12} \hat n_1 \hat n_2 - \frac{\hbar \Omega}{2}\big(\hat \psi_1^\dagger \hat \psi_2 + \hat \psi_2^\dagger \hat \psi_1\big) + \frac{\hbar \delta}{2}(\hat n_1 - \hat n_2) \Bigg\},$$8

The corresponding critical coupling is

$$\hat H = \int d^3 r\, \Bigg\{ \sum_{\sigma=1,2}\hat \psi_\sigma^\dagger\Big(-\frac{\hbar^2\nabla^2}{2m}+V(\mathbf r)\Big)\hat \psi_\sigma + \frac{g_{11}}{2} \hat n_1^2 + \frac{g_{22}}{2} \hat n_2^2 + g_{12} \hat n_1 \hat n_2 - \frac{\hbar \Omega}{2}\big(\hat \psi_1^\dagger \hat \psi_2 + \hat \psi_2^\dagger \hat \psi_1\big) + \frac{\hbar \delta}{2}(\hat n_1 - \hat n_2) \Bigg\},$$9

Accordingly, coherent coupling suppresses polarization and replaces the usual immiscibility threshold by a magnetic transition controlled by $g_{\sigma\sigma'}=4\pi \hbar^2 a_{\sigma\sigma'}/m$0 (Recati et al., 2021).

The same stabilization appears in the homogeneous limit discussed for optical-lattice superflows. There, phase separation in a Rabi-coupled two-component condensate occurs only if

$g_{\sigma\sigma'}=4\pi \hbar^2 a_{\sigma\sigma'}/m$1

with $g_{\sigma\sigma'}=4\pi \hbar^2 a_{\sigma\sigma'}/m$2. The shift relative to the uncoupled criterion $g_{\sigma\sigma'}=4\pi \hbar^2 a_{\sigma\sigma'}/m$3 expresses the energetic cost imposed by coherent locking on relative-density imbalance (He et al., 2021).

With Rashba-type spin–orbit coupling and harmonic confinement, the phase diagram becomes richer but the Rabi term remains a miscibility-promoting control parameter. In quasi-2D, two interaction regimes were studied: weak–weak ($g_{\sigma\sigma'}=4\pi \hbar^2 a_{\sigma\sigma'}/m$4) and weak–strong ($g_{\sigma\sigma'}=4\pi \hbar^2 a_{\sigma\sigma'}/m$5). The reported ground states include plane-wave, half-quantum-vortex, elongated plane-wave, semi-vortex, mixed-mode, shell-like immiscible, and several stripe phases. For weak–strong interactions, an immiscible-to-miscible transition occurs as $g_{\sigma\sigma'}=4\pi \hbar^2 a_{\sigma\sigma'}/m$6 is increased at fixed Rashba coupling $g_{\sigma\sigma'}=4\pi \hbar^2 a_{\sigma\sigma'}/m$7, and the critical Rabi coupling decreases with increasing $g_{\sigma\sigma'}=4\pi \hbar^2 a_{\sigma\sigma'}/m$8 (Ravisankar et al., 2021).

A recent extension shows that coherent coupling can modify not only phase boundaries but the equation of state itself. For a homogeneous dressed state,

$g_{\sigma\sigma'}=4\pi \hbar^2 a_{\sigma\sigma'}/m$9

with $\Omega$0 and $\Omega$1. In the low-density regime one obtains

$\Omega$2

where $\Omega$3 and $\Omega$4 depend on $\Omega$5 and $\Omega$6; when $\Omega$7, the interaction energy saturates and tends toward an effective coupling

$\Omega$8

This saturation phenomenon was experimentally demonstrated in coherently coupled $\Omega$9 mixtures (Eid et al., 29 Jul 2025).

3. Collective excitations, spectral gaps, and instability channels

The most basic spectral consequence of Rabi coupling is the separation of density and spin channels. In the paramagnetic uniform phase, with free-particle energy $\delta$0, the two Bogoliubov branches are

$\delta$1

$\delta$2

The density branch is gapless with sound speed $\delta$3, whereas the spin branch is gapped,

$\delta$4

The spin gap closes at the magnetic critical point and reopens across it (Recati et al., 2021).

This structure underlies a common physical interpretation: coherent coupling acts as a uniform transverse field that locks the relative phase and suppresses long-wavelength spin supercurrents, while leaving the density phonon intact. In long-wavelength hydrodynamics, the spin sector acquires an explicit mass term proportional to $\delta$5, with $\delta$6 the relative phase (Recati et al., 2021).

When spin–orbit coupling is added, the excitation spectrum acquires additional structure. In quasi-2D Rashba-coupled mixtures, Bogoliubov analysis revealed phonon, maxon, and roton features. Increasing the SO coupling deepens the roton minimum and can generate imaginary frequencies, while increasing the Rabi coupling suppresses the imaginary part of the spectrum and eventually removes rotonization. In the stable regime the relevant eigenvectors are density-like; in unstable regimes they become spin-like, which directly connects instability to out-of-phase fluctuations (Ravisankar et al., 2021).

This separation between energetic and dynamical criteria becomes especially sharp in lattice settings. There, Landau instability is diagnosed by negative eigenvalues of the Hermitian matrix $\delta$7, whereas dynamical instability corresponds to complex eigenfrequencies of the non-Hermitian BdG operator itself. Dynamical instability always implies Landau instability, but not conversely; Rabi-coupled Zeeman lattices provide a particularly clear realization of this distinction (He et al., 2021).

4. Optical lattices, harmonic confinement, and superfluid transport

State-dependent lattices turn Rabi-coupled mixtures into a precise setting for studying superfluid breakdown. In a one-dimensional optical lattice, stationary flows are Bloch states

$\delta$8

and stability is determined by the Bogoliubov spectrum around these nonlinear bands (He et al., 2021).

For spin-independent lattices, the first Brillouin zone is $\delta$9. In the miscible case $\psi_1,\psi_2$0, the critical quasimomenta for dynamical and Landau instabilities coincide at approximately $\psi_1,\psi_2$1 and are essentially independent of $\psi_1,\psi_2$2. In the immiscible case $\psi_1,\psi_2$3, all lowest-band Bloch states are unstable at $\psi_1,\psi_2$4, but once $\psi_1,\psi_2$5 exceeds a threshold, the region near the Brillouin-zone center becomes stable and $\psi_1,\psi_2$6 jumps to approximately $\psi_1,\psi_2$7 (He et al., 2021).

The Zeeman lattice is more distinctive. Because

$\psi_1,\psi_2$8

the linear Hamiltonian with Rabi coupling obeys the half-period symmetry

$\psi_1,\psi_2$9

which extends the first Brillouin zone to $$i\hbar \partial_t \psi_1 = \left(-\frac{\hbar^2\nabla^2}{2m}+V+g_{11} n_1+g_{12} n_2+\frac{\hbar\delta}{2}\right)\psi_1 - \frac{\hbar\Omega}{2}\psi_2,$$0. In this setting the dynamical-instability region in the $$i\hbar \partial_t \psi_1 = \left(-\frac{\hbar^2\nabla^2}{2m}+V+g_{11} n_1+g_{12} n_2+\frac{\hbar\delta}{2}\right)\psi_1 - \frac{\hbar\Omega}{2}\psi_2,$$1 plane assumes a $$i\hbar \partial_t \psi_1 = \left(-\frac{\hbar^2\nabla^2}{2m}+V+g_{11} n_1+g_{12} n_2+\frac{\hbar\delta}{2}\right)\psi_1 - \frac{\hbar\Omega}{2}\psi_2,$$2-like shape, and the critical thresholds separate strongly: $$i\hbar \partial_t \psi_1 = \left(-\frac{\hbar^2\nabla^2}{2m}+V+g_{11} n_1+g_{12} n_2+\frac{\hbar\delta}{2}\right)\psi_1 - \frac{\hbar\Omega}{2}\psi_2,$$3 increases with $$i\hbar \partial_t \psi_1 = \left(-\frac{\hbar^2\nabla^2}{2m}+V+g_{11} n_1+g_{12} n_2+\frac{\hbar\delta}{2}\right)\psi_1 - \frac{\hbar\Omega}{2}\psi_2,$$4 and saturates near $$i\hbar \partial_t \psi_1 = \left(-\frac{\hbar^2\nabla^2}{2m}+V+g_{11} n_1+g_{12} n_2+\frac{\hbar\delta}{2}\right)\psi_1 - \frac{\hbar\Omega}{2}\psi_2,$$5, whereas $$i\hbar \partial_t \psi_1 = \left(-\frac{\hbar^2\nabla^2}{2m}+V+g_{11} n_1+g_{12} n_2+\frac{\hbar\delta}{2}\right)\psi_1 - \frac{\hbar\Omega}{2}\psi_2,$$6 remains near $$i\hbar \partial_t \psi_1 = \left(-\frac{\hbar^2\nabla^2}{2m}+V+g_{11} n_1+g_{12} n_2+\frac{\hbar\delta}{2}\right)\psi_1 - \frac{\hbar\Omega}{2}\psi_2,$$7 and becomes approximately $$i\hbar \partial_t \psi_1 = \left(-\frac{\hbar^2\nabla^2}{2m}+V+g_{11} n_1+g_{12} n_2+\frac{\hbar\delta}{2}\right)\psi_1 - \frac{\hbar\Omega}{2}\psi_2,$$8-independent at large coupling. This is the “significant separation” between Landau and dynamical instabilities emphasized in the lattice study (He et al., 2021).

Matter gratings, defined by $$i\hbar \partial_t \psi_1 = \left(-\frac{\hbar^2\nabla^2}{2m}+V+g_{11} n_1+g_{12} n_2+\frac{\hbar\delta}{2}\right)\psi_1 - \frac{\hbar\Omega}{2}\psi_2,$$9 and $$i\hbar \partial_t \psi_2 = \left(-\frac{\hbar^2\nabla^2}{2m}+V+g_{22} n_2+g_{12} n_1-\frac{\hbar\delta}{2}\right)\psi_2 - \frac{\hbar\Omega}{2}\psi_1,$$0, do not inherit the half-period symmetry. Their instability diagrams resemble the spin-independent case, but $$i\hbar \partial_t \psi_2 = \left(-\frac{\hbar^2\nabla^2}{2m}+V+g_{22} n_2+g_{12} n_1-\frac{\hbar\delta}{2}\right)\psi_2 - \frac{\hbar\Omega}{2}\psi_1,$$1 is slightly larger than $$i\hbar \partial_t \psi_2 = \left(-\frac{\hbar^2\nabla^2}{2m}+V+g_{22} n_2+g_{12} n_1-\frac{\hbar\delta}{2}\right)\psi_2 - \frac{\hbar\Omega}{2}\psi_1,$$2, and both increase with $$i\hbar \partial_t \psi_2 = \left(-\frac{\hbar^2\nabla^2}{2m}+V+g_{22} n_2+g_{12} n_1-\frac{\hbar\delta}{2}\right)\psi_2 - \frac{\hbar\Omega}{2}\psi_1,$$3 before saturating. The effect is stronger for $$i\hbar \partial_t \psi_2 = \left(-\frac{\hbar^2\nabla^2}{2m}+V+g_{22} n_2+g_{12} n_1-\frac{\hbar\delta}{2}\right)\psi_2 - \frac{\hbar\Omega}{2}\psi_1,$$4, because the Rabi term competes directly with the interaction term favoring population imbalance (He et al., 2021).

Harmonic traps add a different phenomenology. In effective one-dimensional SO- and Rabi-coupled condensates, “the harmonic trap causes a strong reduction of the multi-peak nature of the condensate and it increases its density.” For repulsive interactions, increasing the SO coupling produces a less dense condensate with more multi-peak structure, whereas increasing positive Rabi coupling increases density with an almost constant number of peaks. Negative Rabi coupling generates a central notch and produces profiles resembling dark-in-bright solitons for both repulsive and attractive interactions (Chiquillo, 2018).

5. Magnetic solitons, bright solitons, and driven nonlinear dynamics

The relative-phase locking imposed by Rabi coupling supports nonlinear spin excitations that are absent in uncoupled mixtures. In weakly polarized, weakly Rabi-coupled binary condensates with $$i\hbar \partial_t \psi_2 = \left(-\frac{\hbar^2\nabla^2}{2m}+V+g_{22} n_2+g_{12} n_1-\frac{\hbar\delta}{2}\right)\psi_2 - \frac{\hbar\Omega}{2}\psi_1,$$5, two classes of magnetic solitons were identified. The first are $$i\hbar \partial_t \psi_2 = \left(-\frac{\hbar^2\nabla^2}{2m}+V+g_{22} n_2+g_{12} n_1-\frac{\hbar\delta}{2}\right)\psi_2 - \frac{\hbar\Omega}{2}\psi_1,$$6 solitons, characterized by a $$i\hbar \partial_t \psi_2 = \left(-\frac{\hbar^2\nabla^2}{2m}+V+g_{22} n_2+g_{12} n_1-\frac{\hbar\delta}{2}\right)\psi_2 - \frac{\hbar\Omega}{2}\psi_1,$$7 jump of the relative phase; the static Son–Stephanov domain wall is the zero-velocity, zero-magnetization limit of this family. The second are $$i\hbar \partial_t \psi_2 = \left(-\frac{\hbar^2\nabla^2}{2m}+V+g_{22} n_2+g_{12} n_1-\frac{\hbar\delta}{2}\right)\psi_2 - \frac{\hbar\Omega}{2}\psi_1,$$8 solitons, which do not exhibit any asymptotic phase jump. In harmonic traps, a moving $$i\hbar \partial_t \psi_2 = \left(-\frac{\hbar^2\nabla^2}{2m}+V+g_{22} n_2+g_{12} n_1-\frac{\hbar\delta}{2}\right)\psi_2 - \frac{\hbar\Omega}{2}\psi_1,$$9 soliton can evolve into a $n_\sigma=|\psi_\sigma|^2$0 soliton and back again while oscillating about the trap center (Qu et al., 2016).

Within the reduced one-dimensional theory for traveling waves, the static balanced limit gives a sine–Gordon equation for the relative phase,

$n_\sigma=|\psi_\sigma|^2$1

with kink solution

$n_\sigma=|\psi_\sigma|^2$2

and characteristic width

$n_\sigma=|\psi_\sigma|^2$3

For the $n_\sigma=|\psi_\sigma|^2$4 family, the termination of the branch occurs at the Landau critical velocity

$n_\sigma=|\psi_\sigma|^2$5

The effective mass can be positive or negative depending on the branch and on $n_\sigma=|\psi_\sigma|^2$6 (Qu et al., 2016).

Attractive and nonautonomous settings support a different class of coherent structures. In the Manakov limit $n_\sigma=|\psi_\sigma|^2$7, a quasi-1D Rabi-coupled system with time-dependent trap can be reduced by a unitary rotation and similarity transformation to the autonomous Manakov model. Exact bright vector solitons follow, and the integrability condition is the Riccati-type relation

$n_\sigma=|\psi_\sigma|^2$8

For the resulting non-autonomous bright solitons, the normalized momentum equals the soliton velocity, $n_\sigma=|\psi_\sigma|^2$9, so the Quintero criterion gives $\sigma=1,2$00, implying stability for the cases analyzed (Kanna et al., 2017).

Spin–orbit-coupled bright solitons acquire further structure under Rabi mixing. In the integrable SOC-only case, exact bright-soliton envelopes exist under a trap–nonlinearity constraint; when Rabi coupling is added, integrability is lost but numerical simulations show that the reinforcement of Rabi coupling generates a stripe phase of the bright solitons and does not destabilize them in the parameter regimes studied (Vinayagam et al., 2017). A related mechanism appears for self-attractive pseudo-spin-$\sigma=1,2$01 condensates in static random potentials: a uniform Rabi field redistributes density between components, changes the spatial profile, and thereby produces a net force through the Ehrenfest equation

$\sigma=1,2$02

so that spin dynamics is converted into center-of-mass motion when interactions are not SU(2)-symmetric (Mardonov et al., 2018).

6. Experimental realizations, numerical frameworks, and extensions

The canonical experimental platforms are alkali gases with two hyperfine states serving as the pseudo-spin-$\sigma=1,2$03 degree of freedom. In $\sigma=1,2$04, spin-dependent lattices can be engineered with far-detuned, linearly polarized laser beams, while Rabi coupling is implemented by microwave or radio-frequency transitions. The principal observables include critical velocities in moving lattices, quench growth of unstable modes, momentum-space spectroscopy of Bogoliubov branches, current decay, and direct observation of spatial demixing (He et al., 2021).

A recent experiment on $\sigma=1,2$05 provided a direct demonstration of coherent-coupling-induced saturation of the interaction energy. Using the states $\sigma=1,2$06 and $\sigma=1,2$07 at $\sigma=1,2$08, with scattering lengths $\sigma=1,2$09, $\sigma=1,2$10, and $\sigma=1,2$11, the system was probed over $\sigma=1,2$12. The longitudinal expansion displayed a pronounced minimum near $\sigma=1,2$13, and the width of the low-interaction region broadened as $\sigma=1,2$14 decreased, in agreement with the predicted crossover into the saturating regime (Eid et al., 29 Jul 2025).

Theoretical work in this area has relied heavily on split-step Crank–Nicolson schemes, imaginary-time propagation for stationary states, real-time propagation for dynamics, and Bogoliubov–de Gennes diagonalization for excitation spectra. Representative studies used such methods for harmonic SO- and Rabi-coupled NPSEs and GPEs, and for quasi-2D Rashba-coupled excitation spectra (Chiquillo, 2018, Ravisankar et al., 2021). For spin-1 condensates, an OpenMP Fortran solver was developed for rotating two-dimensional systems with Rashba or Dresselhaus SO coupling and Rabi mixing, using split-step Crank–Nicolson propagation in both imaginary and real time (Muruganandam et al., 2021).

Rabi coupling is not confined to pseudo-spin-$\sigma=1,2$15. In spin-1 condensates, the coupling acts through $\sigma=1,2$16, mixes the $\sigma=1,2$17 manifolds, and participates in vortex formation, rotation physics, and spin-texture dynamics (Muruganandam et al., 2021). A further extension uses co-propagating Laguerre–Gaussian beams to generate spin–orbit–angular-momentum coupling in a spin-1 condensate; in that setting, a quench of the quadratic Zeeman shift induces many-body Rabi oscillations between distinct quantum phases in the three-minima regime of the single-particle spectrum (Chen et al., 2015).

The broader theoretical program remains active. Open questions explicitly identified for coherently coupled mixtures include critical dynamics across the magnetic transition, dissipative spin hydrodynamics, beyond-mean-field droplet physics, and the controlled dynamics of topological excitations such as solitons, vortices, and confined vortex molecules (Recati et al., 2021). Taken together, these developments establish Rabi-coupled spinor Bose–Einstein condensates as a unifying arena in which coherent internal-state mixing, interaction asymmetry, and external structuring fields can be tuned independently, yielding a wide spectrum of magnetic, hydrodynamic, and nonlinear phenomena.