Numerical Simulations on Stationary States for Rotating Two-Component Bose-Einstein Condensates

We numerically study stationary states such as ground, symmetric and central vortex states as well as their energy diagrams for rotating two-component Bose-Einstein condensates (BECs), which are stationary solutions of time-independent coupled Gross-Pitaevskii equations with an angular momentum rotational term. We compute these stationary states by using normalized gradient flows with a backward Euler finite difference discretization, which is proved to be efficient and accurate. By using this discretization, we find various ground state configurations with several vortices for the two-component BECs, which are novel but not found in rotating two-component BECs; we also find that the critical angular velocity at which the ground states lose symmetry.