XYZ Quantum Heisenberg Models with p-Orbital Bosons

We demonstrate how the spin-1/2 XYZ quantum Heisenberg model can be realized with bosonic atoms loaded in the p band of an optical lattice in the Mott regime. The combination of Bose statistics and the symmetry of the p-orbital wave functions leads to a nonintegrable Heisenberg model with antiferro-magnetic couplings. Moreover, the sign and relative strength of the couplings characterizing the model are shown to be experimentally tunable. We display the rich phase diagram in the one-dimensional case and discuss finite size effects relevant for trapped systems. Finally, experimental issues related to preparation, manipulation, detection, and imperfections are considered.