where α(β) denote up (down) spins, ζ(η) denote up (down) isospins [i.e., electrons in the rst (second) valley], and the spatial orbitals φ j (x, y) are given by the K lowest-energy solutions of the auxiliary single-particle Hamiltonian in Eq. ( ). For clarity and convenience, these solutions are sorted in ascending energy.