Phase Separation and Dynamics of Two-component Bose-Einstein Condensates

We study the interactions between two atomic species in a binary Bose-Einstein condensate to revisit the conditions for miscibility, oscillatory dynamics between the species, steady state solutions and their stability. By employing a variational approach for a quasi one-dimensional, two-atomic species, condensate we obtain equations of motion for the parameters of each species: amplitude, width, position and phase. A further simplification leads to a reduction of the dynamics into a simple classical Newtonian system where components oscillate in an effective potential with a frequency that depends on the harmonic trap strength and the interspecies coupling parameter. We develop explicit conditions for miscibility that can be interpreted as a phase diagram that depends on the harmonic trap’s strength and the interspecies species coupling parameter. We numerically illustrate the bifurcation scenario whereby non-topological, phase-separated states of increasing complexity emerge out of a symmetric state, as the interspecies coupling is increased. The symmetry-breaking dynamical evolution of some of these states is numerically monitored and the associated asymmetric states are also explored.

💡 Research Summary

The paper investigates the interaction of two atomic species in a binary Bose‑Einstein condensate (BEC) confined in a quasi‑one‑dimensional harmonic trap. Starting from the coupled Gross‑Pitaevskii (GP) equations, the authors assume equal intra‑species interaction strengths (g₁₁ = g₂₂) and identical trapping potentials for both components, which reduces the problem to two coupled nonlinear Schrödinger equations with a single inter‑species coupling parameter g and trap frequency Ω.

A variational approach is employed: each component is approximated by a Gaussian wave packet centered at ±B, with time‑dependent amplitude A, width W, position B, overall phase C, wave‑number D, and chirp E. Substituting this ansatz into the Lagrangian and performing spatial integration yields an effective Lagrangian from which six coupled ordinary differential equations (ODEs) for the parameters are derived (Eqs. 8‑13). These ODEs capture the essential dynamics: amplitude decay (dA/dt = –AE), motion of the centers (dB/dt = D + 2BE), width evolution (dW/dt = 2EW), and the evolution of phase‑related variables.

Steady‑state solutions are obtained by setting all time derivatives to zero. For zero separation (B = 0) the mixed state amplitudes and widths are given analytically (Eqs. 14‑16). For non‑zero separation, transcendental relations (Eqs. 17‑19) determine A, B, and W. The analysis reveals a pitchfork bifurcation in the equilibrium separation B as the inter‑species coupling g is varied. At small g the mixed (B = 0) state is stable; beyond a critical g_cr it loses stability and two symmetry‑broken, phase‑separated states (±B ≠ 0) become stable.

Comparisons between the variational ODE predictions and full numerical integration of the GP equations show qualitative agreement. However, the variational model identifies a subcritical pitchfork (point A) and a nearby saddle‑node (point B), whereas the full GP system exhibits a supercritical pitchfork (point D). The discrepancy is attributed to the limited flexibility of the Gaussian ansatz, which cannot capture fine details near the bifurcation.

An explicit formula for the critical coupling, \