Determining Molecular Ground State with Quantum Imaginary Time Evolution using Broken-Symmetry Wave Function

Abstract: The Hartree-Fock (HF) wave function, commonly used for approximating molecular ground states, becomes nonideal in open shell systems due to the inherent multi-configurational nature of the wave function, limiting accuracy in Quantum Imaginary Time Evolution (QITE). We propose replacing the HF wave function with a spin- and spatial-symmetry broken wave function, enhancing convergence by adding a spin operator, $S^2$ as a penalty term to the original molecular Hamiltonian. We verify that this approach provides good convergence behavior towards the lowest energy eigenstate using direct matrix exponentiation for Imaginary Time Evolution (ITE). Numerical simulations were performed on the hydrogen molecule and a square tetrahydrogen cluster using measurement-assisted unitary approximation in QITE. QITE demonstrates faster convergence to the ground state with broken symmetry (BS) compared to HF, particularly after the molecule exhibits a diradical character of 0.56 for hydrogen. Prior to this point, HF remains more effective, suggesting a transition threshold of diradical character for wave function selection. Additionally, the overlap analysis with Complete Active Space Configuration Interaction (CAS-CI) wave function shows that BS has a larger initial overlap than HF in higher-spin, multi-configurational systems like triple bond dissociation in nitrogen molecule. This method provides a pathway for improved energy simulations in open shell systems, where wave function accuracy significantly impacts downstream quantum algorithms and practical applications in quantum chemistry.