Penalty methods for variational quantum eigensolver

TL;DR

This paper analyzes the effectiveness of penalty terms in the variational quantum eigensolver (VQE), providing theoretical insights and formulas to improve symmetry-resolving energy spectrum calculations on near-term quantum computers.

Contribution

It rigorously evaluates two types of penalty terms in VQE, identifying which works correctly and providing a formula for optimal penalty magnitude.

Findings

One penalty type correctly confines eigenstates to the desired symmetry sector.

The other penalty type can give incorrect results and does not reliably enforce symmetry.

A formula is provided to determine the optimal penalty magnitude for faster VQE convergence.

Abstract

The variational quantum eigensolver (VQE) is a promising algorithm to compute eigenstates and eigenenergies of a given quantum system that can be performed on a near-term quantum computer. Obtaining eigenstates and eigenenergies in a specific symmetry sector of the system is often necessary for practical applications of the VQE in various fields ranging from high energy physics to quantum chemistry. It is common to add a penalty term in the cost function of the VQE to calculate such a symmetry-resolving energy spectrum, but systematic analysis on the effect of the penalty term has been lacking, and the use of the penalty term in the VQE has not been justified rigorously. In this work, we investigate two major types of penalty terms for the VQE that were proposed in the previous studies. We show a penalty term in one of the two types works properly in that eigenstates obtained by the VQE with the penalty term reside in the desired symmetry sector. We further give a convenient formula to determine the magnitude of the penalty term, which may lead to the faster convergence of the VQE. Meanwhile, we prove that the other type of penalty terms does not work for obtaining the target state with the desired symmetry in a rigorous sense and even gives completely wrong results in some cases. We finally provide numerical simulations to validate our analysis. Our results apply to general quantum systems and lay the theoretical foundation for the use of the VQE with the penalty terms to obtain the symmetry-resolving energy spectrum of the system, which fuels the application of a near-term quantum computer.