Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2008.09854 · Quantum Science and Technology

A variational quantum algorithm for Hamiltonian diagonalization

Jinfeng Zeng, Chenfeng Cao, Chao Zhang, Pengxiang Xu, Bei Zeng

Published Jul 6, 2021• 17 citations
arXivHTMLPDF

Abstract

Abstract Hamiltonian diagonalization is at the heart of understanding physical properties and practical applications of quantum systems. It is highly desired to design quantum algorithms that can speedup Hamiltonian diagonalization, especially those can be implemented on near-term quantum devices. In this work, we propose a variational quantum algorithm for Hamiltonian diagonalization (VQHD) of quantum systems, which explores the important physical properties, such as temperature, locality, and correlation, of the system. The key idea is that the thermal states of the system encode the information of eigenvalues and eigenstates of the system Hamiltonian. To obtain the full spectrum of the Hamiltonian, we use a quantum imaginary time evolution algorithm with high temperature, which prepares a thermal state with a small correlation length. With Trotterization, this then allows us to implement each step of imaginary time evolution by a local unitary transformation on only a small number of sites. Diagonalizing these thermal states hence leads to a full knowledge of the Hamiltonian eigensystem. We apply our algorithm to diagonalize local Hamiltonians and return results with high precision. Our VQHD algorithm sheds new light on the applications of near-term quantum computers.

DOI
10.1088/2058-9565/ac11a7

Related papers

Variational Hamiltonian Diagonalization for Dynamical Quantum Simulation

Benjamin Commeau, Marco Cerezo, Zoë Holmes +3 more

Dynamical quantum simulation may be one of the first applications to see quantum advantage. However, the circuit depth of standard Trotterization methods can rapidly exceed the coherence time of noisy quantum computers. This has led to recent proposals for variational approaches to dynamical simulation. In this work, we aim to make variational dynamical simulation even more practical and near-term. We propose a new algorithm called Variational Hamiltonian Diagonalization (VHD), which approximately transforms a given Hamiltonian into a diagonal form that can be easily exponentiated. VHD allows for fast forwarding, i.e., simulation beyond the coherence time of the quantum computer with a fixed-depth quantum circuit. It also removes Trotterization error and allows simulation of the entire Hilbert space. We prove an operational meaning for the VHD cost function in terms of the average simulation fidelity. Moreover, we prove that the VHD cost function does not exhibit a shallow-depth barren plateau, i.e., its gradient does not vanish exponentially. Our proof relies on locality of the Hamiltonian, and hence we connect locality to trainability. Our numerical simulations verify that VHD can be used for fast-forwarding dynamics.

arXiv:2009.02559•Bulletin of the American Physical Society•Sep 5, 2020•42 citations

Variational quantum state diagonalization

Ryan LaRose, Arkin Tikku, Étude O’Neel-Judy +2 more

Abstract Variational hybrid quantum-classical algorithms are promising candidates for near-term implementation on quantum computers. In these algorithms, a quantum computer evaluates the cost of a gate sequence (with speedup over classical cost evaluation), and a classical computer uses this information to adjust the parameters of the gate sequence. Here we present such an algorithm for quantum state diagonalization. State diagonalization has applications in condensed matter physics (e.g., entanglement spectroscopy) as well as in machine learning (e.g., principal component analysis). For a quantum state ρ and gate sequence U, our cost function quantifies how far $$U\rho U^\dagger$$ U ρ U † is from being diagonal. We introduce short-depth quantum circuits to quantify our cost. Minimizing this cost returns a gate sequence that approximately diagonalizes ρ. One can then read out approximations of the largest eigenvalues, and the associated eigenvectors, of ρ. As a proof-of-principle, we implement our algorithm on Rigetti’s quantum computer to diagonalize one-qubit states and on a simulator to find the entanglement spectrum of the Heisenberg model ground state.

arXiv:1810.10506•npj Quantum Information•Jun 26, 2019•195 citations

Variational quantum algorithms for discovering Hamiltonian spectra

Tyson Jones, Suguru Endo, Sam McArdle +2 more

Calculating the energy spectrum of a quantum system is an important task, for example to analyze reaction rates in drug discovery and catalysis. There has been significant progress in developing algorithms to calculate the ground state energy of molecules on near-term quantum computers. However, calculating excited state energies has attracted comparatively less attention, and it is currently unclear what the optimal method is. We introduce a low depth, variational quantum algorithm to sequentially calculate the excited states of general Hamiltonians. Incorporating a recently proposed technique [O. Higgott, D. Wang, and S. Brierley, arXiv:1805.08138], we employ the low depth swap test to energetically penalize the ground state, and transform excited states into ground states of modified Hamiltonians. We use variational imaginary time evolution as a subroutine, which deterministically propagates toward the target eigenstate. We discuss how symmetry measurements can mitigate errors in the swap test step. We numerically test our algorithm on Hamiltonians which encode 3SAT optimization problems of up to 18 qubits, and the electronic structure of the lithium hydride molecule. As our algorithm uses only low depth circuits and variational algorithms, it is suitable for use on near-term quantum hardware.

arXiv:1806.05707•Physical Review A•Jun 6, 2019•255 citations

Discovering Hamiltonian spectra with variational quantum imaginary time simulation

Suguru Endo, Tyson Jones, Sam McArdle +2 more

Calculating the energy spectrum of a quantum system is an important task, for example to analyse reaction rates in drug discovery and catalysis. There has been significant progress in developing algorithms to calculate the ground state energy of molecules on near-term quantum computers. However, calculating excited state energies has attracted comparatively less attention, and it is currently unclear what the optimal method is. We introduce a low depth, variational quantum algorithm to sequentially calculate the excited states of general Hamiltonians. Incorporating a recently proposed technique, we employ the low depth swap test to energetically penalise the ground state, and transform excited states into ground states of modified Hamiltonians. We use variational imaginary time evolution as a subroutine, which deterministically propagates towards the target eigenstate. We discuss how symmetry measurements can mitigate errors in the swap test step. We numerically test our algorithm on Hamiltonians which encode 3SAT optimisation problems of up to 18 qubits, and the electronic structure of the Lithium Hydride molecule. As our algorithm uses only low depth circuits and variational algorithms, it is suitable for use on near-term quantum hardware.

•Jun 14, 2018•1 citation

Quantum computing of phonon spectra and thermal properties of crystalline solids

Naman Khandelwal, Bikash K. Behera, Ashok Kumar +1 more

Abstract Variational quantum algorithms offer a promising framework for solving eigenvalue problems on near-term quantum hardware, yet their applicability beyond electronic structure calculations remains relatively unexplored. In this work, we investigate the quantum computing of lattice vibrational and thermodynamical properties by applying the variational quantum eigensolver and variational quantum deflation to phonon Hamiltonians derived from first-principles force constants obtained using density functional theory. The mass-weighted dynamical matrix is mapped onto a qubit-encoded Hermitian operator, enabling computation of the full set of acoustic and optical phonon branches of crystalline silicon and graphene using a reduced qubit register and direct benchmarking against classical diagonalization. The quantum-computed phonon spectrum is further used to evaluate vibrational entropy, constant-volume specific heat, and thermal expansion coefficient, reproducing expected low-temperature quantum behavior and the high-temperature Dulong–Petit limit. We further demonstrate that combined error mitigation strategies help recover phonon dispersions and thermodynamic behavior consistent with expected trends on near-term quantum hardware. Although classical phonon methods remain computationally superior, our results establish phonon-based thermodynamics as a stringent and physically transparent benchmark for assessing variational quantum algorithms on near-term quantum devices.

Physica Scripta•Jun 16, 2026•0 citations

Understanding Quantum Imaginary Time Evolution and its Variational form

Andreu Anglés-Castillo, Luca Ion, Tanmoy Pandit +3 more

Many computationally hard problems can be encoded in quantum Hamiltonians. The solution to these problems is given by the ground states of these Hamiltonians. A state-of-the-art algorithm for finding the ground state of a Hamiltonian is the so-called Quantum Imaginary Time Evolution (QITE) which approximates imaginary time evolution by a unitary evolution that can be implemented in quantum hardware. In this paper, we review the original algorithm together with a comprehensive computer program, as well as, the variational version of it.

arXiv:2510.02015•Oct 2, 2025•0 citations
Abstract Hamiltonian diagonalization is at the heart of understanding physical properties and practical applications of quantum systems. It is highly desired to design quantum algorithms that can speedup Hamiltonian diagonalization, especially those can be implemented on near-term quantum devices. In this work, we propose a variational quantum algorithm for Hamiltonian diagonalization (VQHD) of quantum systems, which explores the important physical properties, such as temperature, locality, and correlation, of the system. The key idea is that the thermal states of the system encode the information of eigenvalues and eigenstates of the system Hamiltonian. To obtain the full spectrum of the Hamiltonian, we use a quantum imaginary time evolution algorithm with high temperature, which prepares a thermal state with a small correlation length. With Trotterization, this then allows us to implement each step of imaginary time evolution by a local unitary transformation on only a small number of sites. Diagonalizing these thermal states hence leads to a full knowledge of the Hamiltonian eigensystem. We apply our algorithm to diagonalize local Hamiltonians and return results with high precision. Our VQHD algorithm sheds new light on the applications of near-term quantum computers.
Dynamical quantum simulation may be one of the first applications to see quantum advantage. However, the circuit depth of standard Trotterization methods can rapidly exceed the coherence time of noisy quantum computers. This has led to recent proposals for variational approaches to dynamical simulation. In this work, we aim to make variational dynamical simulation even more practical and near-term. We propose a new algorithm called Variational Hamiltonian Diagonalization (VHD), which approximately transforms a given Hamiltonian into a diagonal form that can be easily exponentiated. VHD allows for fast forwarding, i.e., simulation beyond the coherence time of the quantum computer with a fixed-depth quantum circuit. It also removes Trotterization error and allows simulation of the entire Hilbert space. We prove an operational meaning for the VHD cost function in terms of the average simulation fidelity. Moreover, we prove that the VHD cost function does not exhibit a shallow-depth barren plateau, i.e., its gradient does not vanish exponentially. Our proof relies on locality of the Hamiltonian, and hence we connect locality to trainability. Our numerical simulations verify that VHD can be used for fast-forwarding dynamics.
Abstract Variational hybrid quantum-classical algorithms are promising candidates for near-term implementation on quantum computers. In these algorithms, a quantum computer evaluates the cost of a gate sequence (with speedup over classical cost evaluation), and a classical computer uses this information to adjust the parameters of the gate sequence. Here we present such an algorithm for quantum state diagonalization. State diagonalization has applications in condensed matter physics (e.g., entanglement spectroscopy) as well as in machine learning (e.g., principal component analysis). For a quantum state ρ and gate sequence U, our cost function quantifies how far UρU†U\rho U^\daggerUρU† U ρ U † is from being diagonal. We introduce short-depth quantum circuits to quantify our cost. Minimizing this cost returns a gate sequence that approximately diagonalizes ρ. One can then read out approximations of the largest eigenvalues, and the associated eigenvectors, of ρ. As a proof-of-principle, we implement our algorithm on Rigetti’s quantum computer to diagonalize one-qubit states and on a simulator to find the entanglement spectrum of the Heisenberg model ground state.
Calculating the energy spectrum of a quantum system is an important task, for example to analyze reaction rates in drug discovery and catalysis. There has been significant progress in developing algorithms to calculate the ground state energy of molecules on near-term quantum computers. However, calculating excited state energies has attracted comparatively less attention, and it is currently unclear what the optimal method is. We introduce a low depth, variational quantum algorithm to sequentially calculate the excited states of general Hamiltonians. Incorporating a recently proposed technique [O. Higgott, D. Wang, and S. Brierley, arXiv:1805.08138], we employ the low depth swap test to energetically penalize the ground state, and transform excited states into ground states of modified Hamiltonians. We use variational imaginary time evolution as a subroutine, which deterministically propagates toward the target eigenstate. We discuss how symmetry measurements can mitigate errors in the swap test step. We numerically test our algorithm on Hamiltonians which encode 3SAT optimization problems of up to 18 qubits, and the electronic structure of the lithium hydride molecule. As our algorithm uses only low depth circuits and variational algorithms, it is suitable for use on near-term quantum hardware.
Calculating the energy spectrum of a quantum system is an important task, for example to analyse reaction rates in drug discovery and catalysis. There has been significant progress in developing algorithms to calculate the ground state energy of molecules on near-term quantum computers. However, calculating excited state energies has attracted comparatively less attention, and it is currently unclear what the optimal method is. We introduce a low depth, variational quantum algorithm to sequentially calculate the excited states of general Hamiltonians. Incorporating a recently proposed technique, we employ the low depth swap test to energetically penalise the ground state, and transform excited states into ground states of modified Hamiltonians. We use variational imaginary time evolution as a subroutine, which deterministically propagates towards the target eigenstate. We discuss how symmetry measurements can mitigate errors in the swap test step. We numerically test our algorithm on Hamiltonians which encode 3SAT optimisation problems of up to 18 qubits, and the electronic structure of the Lithium Hydride molecule. As our algorithm uses only low depth circuits and variational algorithms, it is suitable for use on near-term quantum hardware.
Abstract Variational quantum algorithms offer a promising framework for solving eigenvalue problems on near-term quantum hardware, yet their applicability beyond electronic structure calculations remains relatively unexplored. In this work, we investigate the quantum computing of lattice vibrational and thermodynamical properties by applying the variational quantum eigensolver and variational quantum deflation to phonon Hamiltonians derived from first-principles force constants obtained using density functional theory. The mass-weighted dynamical matrix is mapped onto a qubit-encoded Hermitian operator, enabling computation of the full set of acoustic and optical phonon branches of crystalline silicon and graphene using a reduced qubit register and direct benchmarking against classical diagonalization. The quantum-computed phonon spectrum is further used to evaluate vibrational entropy, constant-volume specific heat, and thermal expansion coefficient, reproducing expected low-temperature quantum behavior and the high-temperature Dulong–Petit limit. We further demonstrate that combined error mitigation strategies help recover phonon dispersions and thermodynamic behavior consistent with expected trends on near-term quantum hardware. Although classical phonon methods remain computationally superior, our results establish phonon-based thermodynamics as a stringent and physically transparent benchmark for assessing variational quantum algorithms on near-term quantum devices.
Many computationally hard problems can be encoded in quantum Hamiltonians. The solution to these problems is given by the ground states of these Hamiltonians. A state-of-the-art algorithm for finding the ground state of a Hamiltonian is the so-called Quantum Imaginary Time Evolution (QITE) which approximates imaginary time evolution by a unitary evolution that can be implemented in quantum hardware. In this paper, we review the original algorithm together with a comprehensive computer program, as well as, the variational version of it.