Spin evolution of spin-1 Bose-Einstein condensates

An analytical formula is obtained to describe the evolution of the average populations of spin components of spin-1 atomic gases. The formula is derived from the exact time-dependent solution of the Hamiltonian $H_{S}=c mathbf{S}^{2}$ without using approximation. Therefore it goes beyond the mean field theory and provides a general, accurate, and complete description for the whole process of non-dissipative evolution starting from various initial states. The numerical results directly given by the formula coincide qualitatively well with existing experimental data, and also with other theoretical results from solving dynamic differential equations. For some special cases of initial state, instead of undergoing strong oscillation as found previously, the evolution is found to go on very steadily in a very long duration.

💡 Research Summary

The paper presents a fully analytical description of the spin dynamics of spin‑1 Bose‑Einstein condensates (BECs) that goes beyond the usual mean‑field treatment. The authors start from the simple yet exact Hamiltonian (H_{S}=c,\mathbf{S}^{2}), where (\mathbf{S}) is the total spin operator of the N‑atom system and (c) characterises the spin‑exchange interaction. Because (\mathbf{S}^{2}) commutes with all components of (\mathbf{S}), its eigenstates (|S,M\rangle) (total spin quantum number S and its projection M) are exact stationary states with energies (E_{S}=c,S(S+1)).

The key methodological step is to expand any physically relevant initial many‑body state (|\Psi(0)\rangle) as a linear combination of these eigenstates, (|\Psi(0)\rangle=\sum_{S,M}C_{S,M}|S,M\rangle). Time evolution is then trivial: each component acquires a phase factor (\exp(-iE_{S}t/\hbar)). Consequently the full time‑dependent wavefunction is known analytically, without any approximation or numerical integration of differential equations.

From this exact wavefunction the authors derive closed‑form expressions for the expectation values of the occupation numbers of the three Zeeman components, (\langle N_{m}(t)\rangle) with (m=+1,0,-1). The result is a sum of cosine terms whose frequencies are differences of the exact eigenenergies, i.e. ((E_{S}-E_{S’})/\hbar=c,