The variational quantum-classical algorithms are the most promising approach for achieving quantum advantage on near-term quantum simulators. Among these methods, the variational quantum eigensolver has attracted a lot of attention in recent years. While it is very effective for simulating the ground state of many-body systems, its generalization to excited states becomes very resource demanding. Here, we show that this issue can significantly be improved by exploiting the symmetries of the Hamiltonian. The improvement is even more effective for higher energy eigenstates. We introduce two methods for incorporating the symmetries. In the first approach, called hardware symmetry preserving, all the symmetries are included in the design of the circuit. In the second approach, the cost function is updated to include the symmetries. The hardware symmetry preserving approach indeed outperforms the second approach. However, integrating all symmetries in the design of the circuit could be extremely challenging. Therefore, we introduce hybrid symmetry preserving method in which symmetries are divided between the circuit and the classical cost function. This allows to harness the advantage of symmetries while preventing sophisticated circuit design.
The variational quantum-classical algorithms are the most promising approach for achieving quantum advantage on near-term quantum simulators. Among these methods, the variational quantum eigensolver has attracted a lot of attention in recent years. While it is very effective for simulating the ground state of many-body systems, its generalization to excited states becomes very resource demanding. Here, we show that this issue can significantly be improved by exploiting the symmetries of the Hamiltonian. The improvement is even more effective for higher energy eigenstates. We introduce two methods for incorporating the symmetries. In the first approach, called hardware symmetry preserving, all the symmetries are included in the design of the circuit. In the second approach, the cost function is updated to include the symmetries. The hardware symmetry preserving approach indeed outperforms the second approach. However, integrating all symmetries in the design of the circuit could be extremely challenging. Therefore, we introduce hybrid symmetry preserving method in which symmetries are divided between the circuit and the classical cost function. This allows to harness the advantage of symmetries while preventing sophisticated circuit design.
The Variational Quantum Eigensolver (VQE), as a hybrid quantum-classical algorithm, is an important tool for effective quantum computing in the current noisy intermediate-scale quantum (NISQ) era. However, the traditional hardware-efficient ansatz without taking into account symmetries requires more computational resources to explore the unnecessary regions in the Hilbert space. The conventional Subspace-Search VQE (SSVQE) algorithm, which can calculate excited states, is also unable to effectively handle degenerate states since the loss function only contains the expectation value of the Hamiltonian. In this study, the energy eigenstates of the one-dimensional Fermi-Hubbard model with two lattice sites and the two-dimensional Hubbard model with four lattice sites are calculated. By incorporating symmetries into the quantum circuits and loss function, we find that both the ground state and excited state calculations are improved greatly compared to the case without symmetries. The enhancement in excited state calculations is particularly significant. This is because quantum circuits that conserve the particle number are used, and appropriate penalty terms are added to the loss function, enabling the optimization process to correctly identify degenerate states. The results are verified through repeated simulations.
Variational quantum algorithms are among the most promising approaches for simulating interacting quantum many-body systems on noisy intermediate-scale quantum (NISQ) devices. However, the practical success of variational quantum eigensolvers (VQE) critically depends on the structure of the chosen variational ansatz. In this work, we investigate the ground-state properties of the one-dimensional antiferromagnetic Heisenberg spin-1/2 chain using both generic hardware-efficient ansatz and physics-informed, symmetry-preserving variational circuits. We benchmark variational results against exact diagonalization and noiseless simulations, and subsequently validate the approach on real IQM Garnet quantum hardware. Our results demonstrate that incorporating physical symmetries into the circuit design leads to significantly improved energy estimates, enhanced robustness against hardware noise, and clearer convergence behavior when compared to hardware-efficient ansatz under identical resource constraints. These findings highlight the importance of problem specific ansatz construction for reliable quantum simulations in the NISQ era.
Determining low-energy eigenstates in electronic many-body quantum systems is a key challenge in computational chemistry and condensed-matter physics. Hybrid quantum-classical approaches, such as the Variational Quantum Eigensolver and Quantum Subspace Methods, offer practical solutions but face limitations in circuit depth and measurement overhead. In this article, we propose a variational strategy based on symmetry-preserving cost functions to iteratively construct a reduced subspace for the extraction of low-lying energy states. We show that, under certain conditions, our approach leads to a tridiagonal representation similar to that obtained with the Lanczos algorithm. The iterative process allows control over the trade-off between circuit depth, the number of variational parameters, and the number of measurements required to achieve the desired accuracy, making it suitable for current quantum hardware. As a proof of concept, we test the proposed algorithms on H4 chain and ring, targeting both the ground-state energy and the charge gap.
computations of the ground and excited electronic states. While substantial progress has been made in designing physically informed quantum algorithms to address the ground state under the constraints of near-term quantum hardware, relatively fewer developments exist for reliably computing excited state energetics, presumably due to the growing complexity of the quantum circuits for excited eigenroots. In this work, we exploit the inherent point-group and spin symmetries of molecular systems to construct symmetry-adapted sets of orthogonal reference states which are evolved with a common symmetry-scalar unitary, enabling the simultaneous determination of ground and excited eigenroots. By minimizing the energy functional over a restricted symmetry-adapted correlated manifold, the weighted subspace-search variational quantum eigensolver enables us to directly target multiple eigenroots within each spatial and spin symmetry sector. This physically motivated selection of references and the unitary confer substantial flexibility to the weighted energy landscape, largely allowing the optimization to avoid numerical traps and converge directly to the targeted eigenstates. More crucially, symmetry resolution also enables access to a larger number of eigenstates without substantially increasing the complexity of the cost function. Through benchmarking on moderate to strongly correlated molecular systems, we demonstrate the accuracy and advantages of our approach relative to several established excited state algorithms on quantum computers.
This work studies the variational quantum eigensolver (VQE) algorithm, which is designed to determine the ground state of a quantum mechanical system by combining classical and quantum hardware. Two methods of reducing the number of required qubit manipulations, prone to induce errors, for the variational quantum eigensolver are studied. First, we formally justify the multiple ℤ2 symmetry qubit reduction scheme first sketched by Bravyi et al. [arXiv:1701.08213 (2017)]. Second, we show that even in small, but non-trivial systems such as H2, LiH, and H2O, the choice of entangling methods (gate based or native) gives rise to varying rates of convergence to the ground state of the system. Through both the removal of qubits and the choice of entangler, the demands on the quantum hardware can be reduced. We find that in general, analyzing the VQE problem is complex, where the number of qubits, the method of entangling, and the depth of the search space all interact. In specific cases however, concrete results can be shown, and an entangling method can be recommended over others as it outperforms in terms of difference from the ground state energy.
We propose a scheme to restore spatial symmetry of Hamiltonian in the variational-quantum-eigensolver (VQE) algorithm for which the quantum circuit structures used usually break the Hamiltonian symmetry. The symmetry-adapted VQE scheme introduced here simply applies the projection operator, which is Hermitian but not unitary, to restore the spatial symmetry in a desired irreducible representation of the spatial group. The entanglement of a quantum state is still represented in a quantum circuit but the nonunitarity of the projection operator is treated classically as postprocessing in the VQE framework. By numerical simulations for a spin-1/2 Heisenberg model on a one-dimensional ring, we demonstrate that the symmetry-adapted VQE scheme with a shallower quantum circuit can achieve significant improvement in terms of the fidelity of the ground state and has a great advantage in terms of the ground-state energy with decent accuracy, as compared to the non-symmetry-adapted VQE scheme. We also demonstrate that the present scheme can approximate low-lying excited states that can be specified by symmetry sectors, using the same circuit structure for the ground-state calculation.