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:2402.12273 · Physical Review Letters

Exact Ansatz of Fermion-Boson Systems for a Quantum Device

Samuel Warren, Yuchen Wang, Carlos L. Benavides-Riveros, David A. Mazziotti

Published Aug 23, 2024• 21 citations
arXivHTMLPDF

Abstract

We present an exact Ansatz for the eigenstate problem of mixed fermion-boson systems that can be implemented on quantum devices. Based on a generalization of the electronic contracted Schrödinger equation (CSE), our approach guides a trial wave function to the ground state of any arbitrary mixed Hamiltonian by directly measuring residuals of the mixed CSE on a quantum device. Unlike density functional and coupled cluster theories applied to electron-phonon or electron-photon systems, the accuracy of our approach is not limited by the unknown exchange-correlation functional or the uncontrolled form of the exponential Ansatz. To test the performance of the method, we study the Tavis-Cummings model, commonly used in polaritonic quantum chemistry. Our results demonstrate that the CSE is a powerful tool in the development of quantum algorithms for solving general fermion-boson many-body problems.

DOI
10.1103/physrevlett.133.080202

Related papers

Verifiably exact solution of the electronic Schrödinger equation on quantum devices

Scott E. Smart, David A. Mazziotti

Quantum computers have the potential for a significant speedup of molecular computations. However, existing algorithms have limitations; quantum phase estimation (QPE) is intractable on current hardware while variational quantum eigensolvers (VQE) are dependent upon approximate wave functions without guaranteed convergence. In this paper we present an algorithm that yields verifiably exact solutions of the many-electron Schr\"odinger equation. Rather than solve the Schr\"odinger equation directly, we solve its contraction over all electrons except two, known as the contracted Schr\"odinger equation (CSE). The CSE generates a wave-function Ansatz, constructed from an iterative product of nonunitary two-body transformations, whose energy gradient with respect to the two-body operator of the current iteration vanishes if and only if the CSE is satisfied. Because the CSE implies the Schr\"odinger equation, the two-electron Ansatz provides a verifiably exact Ansatz for solving the many-electron Schr\"odinger equation. The exactness property contrasts with that of Ans\"atze built from the product of unitary two-body transformations where the gradient---the residual of the anti-Hermitian part of the CSE (ACSE)---can vanish without implying a solution of the Schr\"odinger equation. We demonstrate the algorithm on both simulators and noisy quantum computers with ${\mathrm{H}}_{2}$ dissociation and the rectangle-to-square transition in ${\mathrm{H}}_{4}$.

arXiv:2303.00758•Physical Review A•Feb 5, 2024•18 citations

Quantum simulation of bosons with the contracted quantum eigensolver

Yuchen Wang, LeeAnn M. Sager-Smith, David A. Mazziotti

Abstract Quantum computers are promising tools for simulating many-body quantum systems due to their potential scaling advantage over classical computers. While significant effort has been expended on many-fermion systems, here we simulate a model entangled many-boson system with the contracted quantum eigensolver (CQE). We generalize the CQE to many-boson systems by encoding the bosonic wavefunction on qubits. The CQE provides a compact ansatz for the bosonic wave function whose gradient is proportional to the residual of a contracted Schrödinger equation. We apply the CQE to a bosonic system, where N quantum harmonic oscillators are coupled through a pairwise quadratic repulsion. The model is relevant to the study of coupled vibrations in molecular systems on quantum devices. Results demonstrate the potential efficiency of the CQE in simulating bosonic processes such as molecular vibrations with good accuracy and convergence even in the presence of noise.

New Journal of Physics•Oct 1, 2023•21 citations

Quantum Simulation of Bosons with the Contracted Quantum Eigensolver

Yuchen Wang, LeeAnn M. Sager-Smith, David A. Mazziotti

Quantum computers are promising tools for simulating many-body quantum systems due to their potential scaling advantage over classical computers. While significant effort has been expended on many-fermion systems, here we simulate a model entangled many-boson system with the contracted quantum eigensolver (CQE). We generalize the CQE to many-boson systems by encoding the bosonic wavefunction on qubits. The CQE provides a compact ansatz for the bosonic wave function whose gradient is proportional to the residual of a contracted Schrödinger equation. We apply the CQE to a bosonic system, where $N$ quantum harmonic oscillators are coupled through a pairwise quadratic repulsion. The model is relevant to the study of coupled vibrations in molecular systems on quantum devices. Results demonstrate the potential efficiency of the CQE in simulating bosonic processes such as molecular vibrations with good accuracy and convergence even in the presence of noise.

arXiv:2307.07088•Jul 13, 2023•2 citations

Many-fermion simulation from the contracted quantum eigensolver without fermionic encoding of the wave function

Scott E. Smart, David A. Mazziotti

Quantum computers potentially have an exponential advantage over classical computers for the quantum simulation of many-fermion quantum systems. Nonetheless, fermions are more expensive to simulate than bosons due to the fermionic encoding---a mapping by which the qubits are encoded with fermion statistics. Here we generalize the contracted quantum eigensolver (CQE) to avoid fermionic encoding of the wave function. In contrast to the variational quantum eigensolver, the CQE solves for a many-fermion stationary state by minimizing the contraction (projection) of the Schr\"odinger equation onto two fermions. We avoid fermionic encoding of the wave function by contracting the Schr\"odinger equation onto an unencoded pair of particles. Solution of the resulting contracted equation by a series of unencoded two-body exponential transformations generates an unencoded wave function from which the energy and two-fermion reduced density matrix (2-RDM) can be computed. We apply the unencoded and the encoded CQE algorithms to the hydrogen fluoride molecule, the dissociation of oxygen ${\mathrm{O}}_{2}$, and a series of hydrogen chains. Both algorithms show comparable convergence towards the exact ground-state energies and 2-RDMs, but the unencoded algorithm has computational advantages in terms of state preparation and tomography.

arXiv:2205.01725•Physical Review A•Jun 15, 2022•18 citations

Quantum Many-body Simulations from a Reinforcement-Learned Exponential Ansatz

Yuchen Wang, David A. Mazziotti

Solving for the many-body wavefunction represents a significant challenge on both classical and quantum devices because of the exponential scaling of the Hilbert space with system size. While the complexity of the wavefunction can be reduced through conventional ansätze (e.g., the coupled cluster ansatz), it can still grow rapidly with system size even on quantum devices. An exact, universal two-body exponential ansatz for the many-body wavefunction has been shown to be generated from the solution of the contracted Schrödinger equation (CSE), and recently, this ansatz has been implemented without classical approximation on quantum simulators and devices for the scalable simulation of many-body quantum systems. Here we combine the solution of the CSE with a form of artificial intelligence known as reinforcement learning (RL) to generate highly compact circuits that implement this ansatz without sacrificing accuracy. As a natural extension of CSE, we reformulate the wavefunction update as a Markovian decision process and train the agent to select the optimal actions at each iteration based upon only the current CSE residual. Compact circuits with high accuracy are achieved for H3 and H4 molecules over a range of molecular geometries.

arXiv:2505.01935•May 3, 2025•0 citations

Simulating Quantum Many-Body States with Neural-Network Exponential Ansatz

Wei Zeng, Jiaji Zhang, Lipeng Chen +1 more

Preparing quantum many-body states on classical or quantum devices is a very challenging task that requires accounting for exponentially large Hilbert spaces. Although this complexity can be managed with exponential ansätze (such as in the coupled-cluster method), these approaches are often tailored to specific systems, which limits their universality. Recent work has shown that the contracted Schrödinger equation enables the construction of universal, formally exact exponential ansätze for quantum many-body physics. However, while the ansatz is capable of resolving arbitrary quantum systems, it still requires a full calculation of its parameters whenever the underlying Hamiltonian changes, even slightly. Here, inspired by recent progress in operator learning, we develop a surrogate neural network solver that generates the exponential ansatz parameters using the Hamiltonian parameters as inputs, eliminating the need for repetitive computations. We illustrate the effectiveness of this approach by training neural networks of several quantum many-body systems, including the Fermi-Hubbard model.

arXiv:2411.07886•Nov 12, 2024•0 citations
We present an exact Ansatz for the eigenstate problem of mixed fermion-boson systems that can be implemented on quantum devices. Based on a generalization of the electronic contracted Schrödinger equation (CSE), our approach guides a trial wave function to the ground state of any arbitrary mixed Hamiltonian by directly measuring residuals of the mixed CSE on a quantum device. Unlike density functional and coupled cluster theories applied to electron-phonon or electron-photon systems, the accuracy of our approach is not limited by the unknown exchange-correlation functional or the uncontrolled form of the exponential Ansatz. To test the performance of the method, we study the Tavis-Cummings model, commonly used in polaritonic quantum chemistry. Our results demonstrate that the CSE is a powerful tool in the development of quantum algorithms for solving general fermion-boson many-body problems.
Quantum computers have the potential for a significant speedup of molecular computations. However, existing algorithms have limitations; quantum phase estimation (QPE) is intractable on current hardware while variational quantum eigensolvers (VQE) are dependent upon approximate wave functions without guaranteed convergence. In this paper we present an algorithm that yields verifiably exact solutions of the many-electron Schr\"odinger equation. Rather than solve the Schr\"odinger equation directly, we solve its contraction over all electrons except two, known as the contracted Schr\"odinger equation (CSE). The CSE generates a wave-function Ansatz, constructed from an iterative product of nonunitary two-body transformations, whose energy gradient with respect to the two-body operator of the current iteration vanishes if and only if the CSE is satisfied. Because the CSE implies the Schr\"odinger equation, the two-electron Ansatz provides a verifiably exact Ansatz for solving the many-electron Schr\"odinger equation. The exactness property contrasts with that of Ans\"atze built from the product of unitary two-body transformations where the gradient---the residual of the anti-Hermitian part of the CSE (ACSE)---can vanish without implying a solution of the Schr\"odinger equation. We demonstrate the algorithm on both simulators and noisy quantum computers with H2{\mathrm{H}}_{2}H2​ dissociation and the rectangle-to-square transition in H4{\mathrm{H}}_{4}H4​.
Abstract Quantum computers are promising tools for simulating many-body quantum systems due to their potential scaling advantage over classical computers. While significant effort has been expended on many-fermion systems, here we simulate a model entangled many-boson system with the contracted quantum eigensolver (CQE). We generalize the CQE to many-boson systems by encoding the bosonic wavefunction on qubits. The CQE provides a compact ansatz for the bosonic wave function whose gradient is proportional to the residual of a contracted Schrödinger equation. We apply the CQE to a bosonic system, where N quantum harmonic oscillators are coupled through a pairwise quadratic repulsion. The model is relevant to the study of coupled vibrations in molecular systems on quantum devices. Results demonstrate the potential efficiency of the CQE in simulating bosonic processes such as molecular vibrations with good accuracy and convergence even in the presence of noise.
Quantum computers are promising tools for simulating many-body quantum systems due to their potential scaling advantage over classical computers. While significant effort has been expended on many-fermion systems, here we simulate a model entangled many-boson system with the contracted quantum eigensolver (CQE). We generalize the CQE to many-boson systems by encoding the bosonic wavefunction on qubits. The CQE provides a compact ansatz for the bosonic wave function whose gradient is proportional to the residual of a contracted Schrödinger equation. We apply the CQE to a bosonic system, where NNN quantum harmonic oscillators are coupled through a pairwise quadratic repulsion. The model is relevant to the study of coupled vibrations in molecular systems on quantum devices. Results demonstrate the potential efficiency of the CQE in simulating bosonic processes such as molecular vibrations with good accuracy and convergence even in the presence of noise.
Quantum computers potentially have an exponential advantage over classical computers for the quantum simulation of many-fermion quantum systems. Nonetheless, fermions are more expensive to simulate than bosons due to the fermionic encoding---a mapping by which the qubits are encoded with fermion statistics. Here we generalize the contracted quantum eigensolver (CQE) to avoid fermionic encoding of the wave function. In contrast to the variational quantum eigensolver, the CQE solves for a many-fermion stationary state by minimizing the contraction (projection) of the Schr\"odinger equation onto two fermions. We avoid fermionic encoding of the wave function by contracting the Schr\"odinger equation onto an unencoded pair of particles. Solution of the resulting contracted equation by a series of unencoded two-body exponential transformations generates an unencoded wave function from which the energy and two-fermion reduced density matrix (2-RDM) can be computed. We apply the unencoded and the encoded CQE algorithms to the hydrogen fluoride molecule, the dissociation of oxygen O2{\mathrm{O}}_{2}O2​, and a series of hydrogen chains. Both algorithms show comparable convergence towards the exact ground-state energies and 2-RDMs, but the unencoded algorithm has computational advantages in terms of state preparation and tomography.
Solving for the many-body wavefunction represents a significant challenge on both classical and quantum devices because of the exponential scaling of the Hilbert space with system size. While the complexity of the wavefunction can be reduced through conventional ansätze (e.g., the coupled cluster ansatz), it can still grow rapidly with system size even on quantum devices. An exact, universal two-body exponential ansatz for the many-body wavefunction has been shown to be generated from the solution of the contracted Schrödinger equation (CSE), and recently, this ansatz has been implemented without classical approximation on quantum simulators and devices for the scalable simulation of many-body quantum systems. Here we combine the solution of the CSE with a form of artificial intelligence known as reinforcement learning (RL) to generate highly compact circuits that implement this ansatz without sacrificing accuracy. As a natural extension of CSE, we reformulate the wavefunction update as a Markovian decision process and train the agent to select the optimal actions at each iteration based upon only the current CSE residual. Compact circuits with high accuracy are achieved for H3 and H4 molecules over a range of molecular geometries.
Preparing quantum many-body states on classical or quantum devices is a very challenging task that requires accounting for exponentially large Hilbert spaces. Although this complexity can be managed with exponential ansätze (such as in the coupled-cluster method), these approaches are often tailored to specific systems, which limits their universality. Recent work has shown that the contracted Schrödinger equation enables the construction of universal, formally exact exponential ansätze for quantum many-body physics. However, while the ansatz is capable of resolving arbitrary quantum systems, it still requires a full calculation of its parameters whenever the underlying Hamiltonian changes, even slightly. Here, inspired by recent progress in operator learning, we develop a surrogate neural network solver that generates the exponential ansatz parameters using the Hamiltonian parameters as inputs, eliminating the need for repetitive computations. We illustrate the effectiveness of this approach by training neural networks of several quantum many-body systems, including the Fermi-Hubbard model.