State preparation and symmetries

We demonstrate the importance of symmetries in Variational Quantum Eigensolver (VQE) algorithms to prepare the ground or specific low-lying states of quantum Hamiltonians. We examine two spin problems, one with random all-to-all couplings inspired by neutrino flavor evolution in supernovae, and the standard Heisenberg spin Hamiltonian on a $4 \times 3$ lattice. The neutrino Hamiltonian has the total spin $J$ and third component $J_{\rm{z}}$ as its only symmetries. The Heisenberg model has these symmetries plus translational invariance and reflection symmetry. We demonstrate that the convergence of variational methods is dramatically improved by keeping all symmetries. In both cases a nearly exact solution can be obtained in cases where standard unconstrained variational algorithms fail. Since variational algorithms can use standard Trotter steps as part of the optimization, allowing additional correlations that obey all the symmetries of the Hamiltonian will speed convergence of variational algorithms. This will lead to faster convergence than standard projection algorithms.