We investigate the ground state of the
SU(4) Heisenberg model on a square lattice with spatial anisotropy on each plaquette bond using the tensor-network method based on infinite projected entangled pair states. We find that the
SU(4) singlet ground state appears in the strongly anisotropic limit, whereas N\'eel and valence-bond crystal orders coexist in the nearly isotropic limit. By examining the intermediate parameter region, we identify a phase transition between these phases. The nature of the phase transition is likely to be of first order, and the transition point is estimated to be around ${J}^{\ensuremath{'}}/J\ensuremath{\approx}0.85(5)$, where
J and ${J}^{\ensuremath{'}}$ are the interaction strengths of intra- and interplaquette bonds, respectively. We also calculate the anisotropy dependence of singlet correlations on a plaquette bond, which will be useful for future experiments of ultracold atoms in optical lattices.
The SU(4) Heisenberg model can serve as a low-energy model of the Mott insulating state in materials where the spins and orbitals are highly symmetric or in systems of alkaline-earth atoms on an optical lattice. Recently, it has been argued that on the honeycomb lattice, the model exhibits a unique spin-orbital liquid phase with an algebraic decay of correlations [P. Corboz et al., Phys. Rev. X 2, 041013 (2012)]. Here we study the instability of the algebraic spin-orbital liquid toward the spontaneous formation of SU(4) singlet plaquettes (tetramerization). Using a variational Monte Carlo approach to evaluate the projected wave function of fermions with $\ensuremath{\pi}$-flux state, we find that the algebraic liquid is robust and that a small finite value of the next-nearest exchange is needed to induce tetramerization. We also studied the phase diagram of a model which interpolates between the nearest-neighbor Heisenberg model and a Hamiltonian for which the singlet-plaquette-product state is an exact ground state.