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:2610.18681

Explicitly correlated ground states with a feedback-free variational hybrid algorithm

Benjamin Corbett, Akimasa Miyake

quant-phcond-mat.str-elSubmitted Oct 13, 2026• Announced Oct 14, 2026• 0 citations
arXivHTMLPDF

Abstract

Explicitly correlated methods, which represent dominant short-range correlations by a diagonal Jastrow-type correlator, are designed to accurately describe systems of strongly correlated electrons, such as the fully periodic 2D Fermi--Hubbard model. However, as the reference state is compact in momentum space whereas the correlator acts naturally in real space, there is a fundamental mismatch of bases that classical algorithms resolve only by restricting the reference or the correlator. Here we propose a hybrid algorithm in which a quantum processor prepares a low-entanglement momentum-space matrix product state (MPS), precomputed classically by the transcorrelated method at the cost of a nonvariational energy, rotates it to real space, and measures in bases fixed by the Hamiltonian. Because the measurement bases are independent of the correlator, the correlator can be variationally optimized entirely in classical post-processing with no quantum-classical feedback loop. We demonstrate our algorithm on 32 qubits of Quantinuum's H2-2 with a bond-dimension-4 MPS for the half-filled $4 \times 4$ lattice and achieve a relative ground-state energy error of $16 \pm 4\%$ after postselection, reduced to $3.5 \pm 0.5\%$ by calibration using low-order moments from the MPS. Compared with other feedback-free hybrid algorithms, even under generous assumptions, our method requires fewer samples than quantum-selected configuration interaction and fewer two-qubit gates than quantum Krylov methods.

Related papers

Basis Adaptive Algorithm for Quantum Many-Body Systems on Quantum Computers

Anutosh Biswas, Sayan Ghosh, Ritajit Majumdar +3 more

We introduce a Basis Adaptive (BA) algorithm for hybrid quantum-classical simulation of correlated quantum many-body systems. Starting from a small set of physically motivated bitstrings, the algorithm iteratively applies a single-step first-order Trotterized circuit on a quantum processor, filters the sampled configurations by enforcing $U(1)$ spin conservation and lattice reflection symmetry, and classically diagonalizes the Hamiltonian in the resulting reduced Hilbert space. This design avoids the variational optimization overhead of VQE, the deep coherent circuits required by QPE, and the symmetry-violating subspaces that arise in SKQD. The ground-state energy error is bounded analytically by $\sqrt{8}\,\|H\|\left(1-\sqrt{α_{D_T}}\right)^{1/2}$, where $α_{D_T}$ is the probability weight captured by the $D_T$ sampled basis states. This bound connects algorithm performance directly to ground-state sparsity and explains the observed accuracy hierarchy across different phases. Benchmarked on the spin-$1/2$ Heisenberg XXZ chain (up to $N=62$ qubits on the IBM Heron processor), the algorithm achieves a $3.5\%$ energy error in the gapped Neel phase ($Δ=2.0$) and below $0.5\%$ at the ferromagnetic boundary ($Δ=-1.0$). The accuracy degrades to $28.7\%$ in the strongly quasi-long-range-ordered regime ($Δ=0.5$). Spin-spin correlation functions are reproduced across all regimes, confirming that symmetry-filtered real-time sampling provides a practical and noise-resilient pathway to ground-state properties on near-term quantum hardware.

cond-mat.str-elCross-listedarXiv:2512.12753•Dec 16, 2025•0 citations

Unraveling correlated materials' properties with noisy quantum computers: Natural-orbitalized variational quantum eigensolving of extended impurity models within a slave-boson approach

Pauline Besserve, Thomas Ayral

We propose a method for computing space-resolved correlation properties of the two-dimensional Hubbard model within a quantum-classical embedding strategy that uses a Noisy, Intermediate Scale Quantum (NISQ) computer to solve the embedded model. While previous approaches were limited to purely local, one-impurity embedded models, requiring at most four qubits and relatively shallow circuits, we solve a two-impurity model requiring eight qubits with an advanced hybrid scheme on top of the Variational Quantum Eigensolver algorithm. This iterative scheme, dubbed Natural Orbitalization (NOization), gradually transforms the single-particle basis to the approximate Natural-Orbital basis, in which the ground state can be minimally expressed, at the cost of measuring the one-particle reduced density-matrix of the embedded problem. We show that this transformation tends to make the variational optimization of existing (but too deep) ansatz circuits faster and more accurate, and we propose an ansatz, the Multi-Reference Excitation Preserving (MREP) ansatz, that achieves great expressivity without requiring a prohibitive gate count. The one-impurity version of the ansatz has only one parameter, making the ground state preparation a trivial step, which supports the optimal character of our approach. Within a Rotationally Invariant Slave Boson embedding scheme that requires a minimal number of bath sites and does not require computing the full Green's function, the NOization combined with the MREP ansatz allow us to compute accurate, space-resolved quasiparticle weights and static self-energies for the Hubbard model even in the presence of noise levels representative of current NISQ processors. This paves the way to a controlled solution of the Hubbard model with larger and larger embedded problems solved by quantum computers.

•Aug 24, 2021•1 citation

Quantum simulation of electronic structure with a transcorrelated Hamiltonian: improved accuracy with a smaller footprint on the quantum computer

Mario Motta, Tanvi P. Gujarati, Julia E. Rice +6 more

Quantum simulations of electronic structure with a transformed Hamiltonian that includes some electron correlation effects are demonstrated. The transcorrelated Hamiltonian used in this work is efficiently constructed classically, at polynomial cost, by an approximate similarity transformation with an explicitly correlated two-body unitary operator. This Hamiltonian is Hermitian, includes no more than two-particle interactions, and is free of electron-electron singularities. We investigate the effect of such a transformed Hamiltonian on the accuracy and computational cost of quantum simulations by focusing on a widely used solver for the Schrödinger equation, namely the variational quantum eigensolver method, based on the unitary coupled cluster with singles and doubles (q-UCCSD) Ansatz. Nevertheless, the formalism presented here translates straightforwardly to other quantum algorithms for chemistry. Our results demonstrate that a transcorrelated Hamiltonian, paired with extremely compact bases, produces explicitly correlated energies comparable to those from much larger bases. For the chemical species studied here, explicitly correlated energies based on an underlying 6-31G basis had cc-pVTZ quality. The use of the very compact transcorrelated Hamiltonian reduces the number of CNOT gates required to achieve cc-pVTZ quality by up to two orders of magnitude, and the number of qubits by a factor of three.

arXiv:2006.02488•Physical Chemistry Chemical Physics•Jan 1, 2020•78 citations

Scaling up the transcorrelated density matrix renormalization group

Benjamin C. Corbett, Akimasa Miyake

Explicitly correlated methods, such as the transcorrelated method which shifts a Jastrow or Gutzwiller correlator from the wave function to the Hamiltonian, are designed for high-accuracy calculations of electronic structures, but their application to larger systems has been hampered by the computational cost. We develop improved techniques for the transcorrelated density-matrix renormalization group (DMRG), in which the ground state of the transcorrelated Hamiltonian is represented as a matrix product state (MPS), and demonstrate large-scale calculations of the ground-state energy of the two-dimensional Fermi-Hubbard model. Our developments stem from three technical inventions: (i) constructing matrix product operators (MPOs) of transcorrelated Hamiltonians with low bond dimension and high sparsity, (ii) exploiting the entanglement structure of the ground states to increase the accuracy of the MPS representation, and (iii) optimizing the nonlinear parameter of the Gutzwiller correlator to mitigate the nonvariational nature of the transcorrelated method. We examine systems of size up to $12\ifmmode\times\else\texttimes\fi{}12$ lattice sites, four times larger than previous transcorrelated DMRG studies, and demonstrate that transcorrelated DMRG yields significant improvements over standard nontranscorrelated DMRG for equivalent computational effort. Transcorrelated DMRG reduces the error of the ground-state energy by $2.4\ifmmode\times\else\texttimes\fi{}$--$14\ifmmode\times\else\texttimes\fi{}$, with the smallest improvement seen for a small system at half filling and the largest improvement in a dilute closed-shell system.

arXiv:2506.07441•Physical review. B./Physical review. B•Sep 15, 2025•1 citation

Imaginary-time-enhanced feedback-based quantum algorithms for universal ground-state preparation

Thanh Nguyen Van Long, Lan Nguyen Tran, Le Bin Ho

Preparing ground states of strongly correlated quantum systems is a central goal in quantum simulation and optimization. The feedback-based quantum algorithm (FALQON) provides an attractive alternative to variational methods with a fully quantum feedback rule, but it fails in the presence of spectral degeneracies, where the feedback signal collapses and the evolution cannot reach the ground state. Using the Fermi-Hubbard model on lattices up to 3x3, we show that this breakdown appears at half-filling on the 2x2 lattice and extends to both half-filled and doped configurations on the 3x3 lattice. We then introduce an imaginary-time-enhanced FALQON ($i$FALQON) scheme, which inserts short imaginary-time evolution steps into the feedback loop. This proposed method suppresses excited-state components, escapes degenerate subspaces, restores monotonic energy descent and improves ground-state convergence across the tested fillings. We further benchmark its performance against existing methods, showing improved results for iFALQON, demonstrating the potential of hybrid real-imaginary-time feedback dynamics for ground-state preparation in strongly correlated quantum systems.

quant-pharXiv:2512.13044•Dec 16, 2025•0 citations

Unraveling correlated material properties with noisy quantum computers: Natural orbitalized variational quantum eigensolving of extended impurity models within a slave-boson approach

Pauline Besserve, Thomas Ayral

We propose a method for computing space-resolved correlation properties of the two-dimensional Hubbard model within a quantum-classical embedding strategy that uses a noisy, intermediate scale quantum (NISQ) computer to solve the embedded model. While previous approaches were limited to purely local, one-impurity embedded models, requiring at most four qubits and relatively shallow circuits, we solve a two-impurity model requiring eight qubits with an advanced hybrid scheme on top of the variational quantum eigensolver algorithm. This iterative scheme, dubbed natural orbitalization (NOization), gradually transforms the single-particle basis to the approximate natural-orbital basis, in which the ground state can be minimally expressed, at the cost of measuring the one-particle reduced density matrix of the embedded problem. We show that this transformation tends to make the variational optimization of existing (but too deep) Ansatz circuits faster and more accurate, and we propose an Ansatz, the multireference excitation-preserving (MREP) Ansatz, that achieves great expressivity without requiring a prohibitive gate count, thus bridging the gap between hardware-efficient and physically motivated strategies in variational Ansatz design. The one-impurity version of the Ansatz has only one parameter, making the ground-state preparation a trivial step, which supports the optimal character of our approach. Within a rotationally invariant slave boson embedding scheme that requires a minimal number of bath sites and does not require computing the full Green's function, the NOization combined with the MREP Ansatz allow us to compute accurate, space-resolved quasiparticle weights and static self-energies for the Hubbard model even in the presence of noise levels representative of current NISQ processors. This paves the way to a controlled solution of the Hubbard model with larger and larger embedded problems solved by quantum computers.

arXiv:2108.10780•Physical review. B./Physical review. B•Mar 7, 2022•14 citations
Explicitly correlated methods, which represent dominant short-range correlations by a diagonal Jastrow-type correlator, are designed to accurately describe systems of strongly correlated electrons, such as the fully periodic 2D Fermi--Hubbard model. However, as the reference state is compact in momentum space whereas the correlator acts naturally in real space, there is a fundamental mismatch of bases that classical algorithms resolve only by restricting the reference or the correlator. Here we propose a hybrid algorithm in which a quantum processor prepares a low-entanglement momentum-space matrix product state (MPS), precomputed classically by the transcorrelated method at the cost of a nonvariational energy, rotates it to real space, and measures in bases fixed by the Hamiltonian. Because the measurement bases are independent of the correlator, the correlator can be variationally optimized entirely in classical post-processing with no quantum-classical feedback loop. We demonstrate our algorithm on 32 qubits of Quantinuum's H2-2 with a bond-dimension-4 MPS for the half-filled 4×44 \times 44×4 lattice and achieve a relative ground-state energy error of 16±4%16 \pm 4\%16±4% after postselection, reduced to 3.5±0.5%3.5 \pm 0.5\%3.5±0.5% by calibration using low-order moments from the MPS. Compared with other feedback-free hybrid algorithms, even under generous assumptions, our method requires fewer samples than quantum-selected configuration interaction and fewer two-qubit gates than quantum Krylov methods.
We introduce a Basis Adaptive (BA) algorithm for hybrid quantum-classical simulation of correlated quantum many-body systems. Starting from a small set of physically motivated bitstrings, the algorithm iteratively applies a single-step first-order Trotterized circuit on a quantum processor, filters the sampled configurations by enforcing U(1)U(1)U(1) spin conservation and lattice reflection symmetry, and classically diagonalizes the Hamiltonian in the resulting reduced Hilbert space. This design avoids the variational optimization overhead of VQE, the deep coherent circuits required by QPE, and the symmetry-violating subspaces that arise in SKQD. The ground-state energy error is bounded analytically by 8 ∥H∥(1−αDT)1/2\sqrt{8}\,\|H\|\left(1-\sqrt{α_{D_T}}\right)^{1/2}8​∥H∥(1−αDT​​​)1/2, where αDTα_{D_T}αDT​​ is the probability weight captured by the DTD_TDT​ sampled basis states. This bound connects algorithm performance directly to ground-state sparsity and explains the observed accuracy hierarchy across different phases. Benchmarked on the spin-1/21/21/2 Heisenberg XXZ chain (up to N=62N=62N=62 qubits on the IBM Heron processor), the algorithm achieves a 3.5%3.5\%3.5% energy error in the gapped Neel phase (Δ=2.0Δ=2.0Δ=2.0) and below 0.5%0.5\%0.5% at the ferromagnetic boundary (Δ=−1.0Δ=-1.0Δ=−1.0). The accuracy degrades to 28.7%28.7\%28.7% in the strongly quasi-long-range-ordered regime (Δ=0.5Δ=0.5Δ=0.5). Spin-spin correlation functions are reproduced across all regimes, confirming that symmetry-filtered real-time sampling provides a practical and noise-resilient pathway to ground-state properties on near-term quantum hardware.
We propose a method for computing space-resolved correlation properties of the two-dimensional Hubbard model within a quantum-classical embedding strategy that uses a Noisy, Intermediate Scale Quantum (NISQ) computer to solve the embedded model. While previous approaches were limited to purely local, one-impurity embedded models, requiring at most four qubits and relatively shallow circuits, we solve a two-impurity model requiring eight qubits with an advanced hybrid scheme on top of the Variational Quantum Eigensolver algorithm. This iterative scheme, dubbed Natural Orbitalization (NOization), gradually transforms the single-particle basis to the approximate Natural-Orbital basis, in which the ground state can be minimally expressed, at the cost of measuring the one-particle reduced density-matrix of the embedded problem. We show that this transformation tends to make the variational optimization of existing (but too deep) ansatz circuits faster and more accurate, and we propose an ansatz, the Multi-Reference Excitation Preserving (MREP) ansatz, that achieves great expressivity without requiring a prohibitive gate count. The one-impurity version of the ansatz has only one parameter, making the ground state preparation a trivial step, which supports the optimal character of our approach. Within a Rotationally Invariant Slave Boson embedding scheme that requires a minimal number of bath sites and does not require computing the full Green's function, the NOization combined with the MREP ansatz allow us to compute accurate, space-resolved quasiparticle weights and static self-energies for the Hubbard model even in the presence of noise levels representative of current NISQ processors. This paves the way to a controlled solution of the Hubbard model with larger and larger embedded problems solved by quantum computers.
Quantum simulations of electronic structure with a transformed Hamiltonian that includes some electron correlation effects are demonstrated. The transcorrelated Hamiltonian used in this work is efficiently constructed classically, at polynomial cost, by an approximate similarity transformation with an explicitly correlated two-body unitary operator. This Hamiltonian is Hermitian, includes no more than two-particle interactions, and is free of electron-electron singularities. We investigate the effect of such a transformed Hamiltonian on the accuracy and computational cost of quantum simulations by focusing on a widely used solver for the Schrödinger equation, namely the variational quantum eigensolver method, based on the unitary coupled cluster with singles and doubles (q-UCCSD) Ansatz. Nevertheless, the formalism presented here translates straightforwardly to other quantum algorithms for chemistry. Our results demonstrate that a transcorrelated Hamiltonian, paired with extremely compact bases, produces explicitly correlated energies comparable to those from much larger bases. For the chemical species studied here, explicitly correlated energies based on an underlying 6-31G basis had cc-pVTZ quality. The use of the very compact transcorrelated Hamiltonian reduces the number of CNOT gates required to achieve cc-pVTZ quality by up to two orders of magnitude, and the number of qubits by a factor of three.
Explicitly correlated methods, such as the transcorrelated method which shifts a Jastrow or Gutzwiller correlator from the wave function to the Hamiltonian, are designed for high-accuracy calculations of electronic structures, but their application to larger systems has been hampered by the computational cost. We develop improved techniques for the transcorrelated density-matrix renormalization group (DMRG), in which the ground state of the transcorrelated Hamiltonian is represented as a matrix product state (MPS), and demonstrate large-scale calculations of the ground-state energy of the two-dimensional Fermi-Hubbard model. Our developments stem from three technical inventions: (i) constructing matrix product operators (MPOs) of transcorrelated Hamiltonians with low bond dimension and high sparsity, (ii) exploiting the entanglement structure of the ground states to increase the accuracy of the MPS representation, and (iii) optimizing the nonlinear parameter of the Gutzwiller correlator to mitigate the nonvariational nature of the transcorrelated method. We examine systems of size up to $12\ifmmode\times\else\texttimes\fi{}12$ lattice sites, four times larger than previous transcorrelated DMRG studies, and demonstrate that transcorrelated DMRG yields significant improvements over standard nontranscorrelated DMRG for equivalent computational effort. Transcorrelated DMRG reduces the error of the ground-state energy by $2.4\ifmmode\times\else\texttimes\fi{}$--$14\ifmmode\times\else\texttimes\fi{}$, with the smallest improvement seen for a small system at half filling and the largest improvement in a dilute closed-shell system.
Preparing ground states of strongly correlated quantum systems is a central goal in quantum simulation and optimization. The feedback-based quantum algorithm (FALQON) provides an attractive alternative to variational methods with a fully quantum feedback rule, but it fails in the presence of spectral degeneracies, where the feedback signal collapses and the evolution cannot reach the ground state. Using the Fermi-Hubbard model on lattices up to 3x3, we show that this breakdown appears at half-filling on the 2x2 lattice and extends to both half-filled and doped configurations on the 3x3 lattice. We then introduce an imaginary-time-enhanced FALQON (iiiFALQON) scheme, which inserts short imaginary-time evolution steps into the feedback loop. This proposed method suppresses excited-state components, escapes degenerate subspaces, restores monotonic energy descent and improves ground-state convergence across the tested fillings. We further benchmark its performance against existing methods, showing improved results for iFALQON, demonstrating the potential of hybrid real-imaginary-time feedback dynamics for ground-state preparation in strongly correlated quantum systems.
We propose a method for computing space-resolved correlation properties of the two-dimensional Hubbard model within a quantum-classical embedding strategy that uses a noisy, intermediate scale quantum (NISQ) computer to solve the embedded model. While previous approaches were limited to purely local, one-impurity embedded models, requiring at most four qubits and relatively shallow circuits, we solve a two-impurity model requiring eight qubits with an advanced hybrid scheme on top of the variational quantum eigensolver algorithm. This iterative scheme, dubbed natural orbitalization (NOization), gradually transforms the single-particle basis to the approximate natural-orbital basis, in which the ground state can be minimally expressed, at the cost of measuring the one-particle reduced density matrix of the embedded problem. We show that this transformation tends to make the variational optimization of existing (but too deep) Ansatz circuits faster and more accurate, and we propose an Ansatz, the multireference excitation-preserving (MREP) Ansatz, that achieves great expressivity without requiring a prohibitive gate count, thus bridging the gap between hardware-efficient and physically motivated strategies in variational Ansatz design. The one-impurity version of the Ansatz has only one parameter, making the ground-state preparation a trivial step, which supports the optimal character of our approach. Within a rotationally invariant slave boson embedding scheme that requires a minimal number of bath sites and does not require computing the full Green's function, the NOization combined with the MREP Ansatz allow us to compute accurate, space-resolved quasiparticle weights and static self-energies for the Hubbard model even in the presence of noise levels representative of current NISQ processors. This paves the way to a controlled solution of the Hubbard model with larger and larger embedded problems solved by quantum computers.