| Summary: | We use a variational order method to explore the ground state phase diagrams of spin–orbit coupled Bose atoms in optical lattices. By properly parameterizing the order with an ansatz function for each possible phase and comparing their energies which are minimized with respect to the variational parameters, we can effectively identify the lowest energy states and depict the phase diagram. While the Mott-insulator-superfluid phase boundary is computed by the usual Gutzwiller method, the spin-ordered phase structures of the deep Mott regime as well as the superfluid regime are studied by the variational order method. For the isotropic spin–orbit coupling, the phase diagram in the deep Mott insulator regime is qualitatively similar to results derived by the classical Monte-Carlo simulations or other numerical methods. The spiral phases with different spatial periods are discussed in detail. We find three new spin-ordered phases in the deep Mott insulator regime with anisotropic spin–orbit coupling. We also identify the uniform superfluid, stripe and checkerboard supersolid phases with exotic spin orders in the superfluid regime.
|
|---|