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:2409.12237 · Physical Review Applied

Compressing Hamiltonians with ab initio downfolding for simulating strongly-correlated materials on quantum computers

Antonios M. Alvertis, Abid Khan, Norm M. Tubman

Published Apr 11, 2025• 7 citations
arXivHTMLPDF

Abstract

The accurate first-principles description of strongly correlated materials is an important and challenging problem in condensed matter physics. Ab initio downfolding has emerged as a way of deriving compressed many-body Hamiltonians that maintain the essential physics of strongly correlated materials. The solution of these material-specific models is still exponentially difficult to generate on classical computers, but quantum algorithms allow for a significant speed-up in obtaining the ground states of these compressed Hamiltonians. Here, we demonstrate that using quantum algorithms to obtain the properties of downfolded Hamiltonians can indeed yield high-fidelity solutions. By combining ab initio downfolding and variational quantum eigensolvers, we correctly predict the antiferromagnetic state of one-dimensional cuprate ${\mathrm{Ca}}_{2}{\mathrm{Cu}\mathrm{O}}_{3}$, the excitonic ground state of monolayer ${\mathrm{W}\mathrm{Te}}_{2}$, and the charge-ordered state of correlated metal ${\mathrm{Sr}\mathrm{VO}}_{3}$. Numerical simulations using a classical tensor network implementation of variational quantum eigensolvers allow us to simulate large models with up to 54 qubits and encompassing up to four bands in the correlated subspace, which is indicative of the complexity that our framework can address. Through these methods we demonstrate the potential of classical preoptimization and downfolding techniques for enabling efficient materials simulation using quantum algorithms.

DOI
10.1103/physrevapplied.23.044028

Related papers

Scalable Ab Initio Quantum Many-Body Methods for Crystalline Materials

Junjie Yang

Achieving materials-specific predictions for large, realistic molecules and solids with strong electron correlations remains a long-standing challenge in quantum chemistry. This dissertation develops scalable ab initio quantum many-body methods for crystalline materials, demonstrating that rigorous optimization of theoretical approximations and numerical algorithms enables predictive computation even in strongly correlated systems. The work tackles three interconnected challenges, namely describing material-specific properties of strongly correlated superconductors, achieving linear scaling with k-points for large-scale periodic calculations, and establishing reliable methods for electron-phonon interactions. The following chapters address each challenge in turn. Chapter 2 demonstrates that ab initio density matrix embedding theory, combined with symmetry-breaking strategies for superconducting order parameters, reproduces experimental trends in cuprate superconductivity, including pressure and layer effects on transition temperatures. By directly solving the electronic Schrödinger equation from material structures, these calculations capture multi-orbital covalency effects that simplified model Hamiltonians neglect, revealing their essential role in the spin fluctuations driving pairing. This predictive capability establishes a foundation for efficient computational screening of superconducting materials. Extending these calculations to larger k-point meshes and broader phase diagram explorations, however, demands more efficient computational algorithms. To overcome this computational bottleneck, Chapter 3 introduces interpolative separable density fitting framework. This development enables thermodynamic limit calculations with up to 1000 k-points, as validated on diamond, carbon dioxide, nickel monoxide, and cuprate systems. Integration with density matrix embedding theory and local natural orbital methods further yields converged ground-state energies from accurate correlated wavefunction methods such as CCSD(T) for crystalline materials in the thermodynamic limit. With computational efficiency established, the remaining challenge lies in extending the theoretical framework beyond purely electronic degrees of freedom. Chapter 4 addresses this need by establishing exponential ansatz wavefunctions as competitive approaches for electron-phonon systems. Unlike previous treatments based on perturbation theory, systematic benchmarking of coupled cluster theory and variational Lang-Firsov methods on the Holstein model demonstrates, for the first time, that these approaches accurately describe polaron formation across coupling regimes when combined with appropriate reference state optimization. These results provide a foundation for incorporating phonon effects into ab initio studies of superconducting materials. In addition to these methodological advances, GPU acceleration and MPI parallelization enhance computational throughput, enabling efficient large-scale calculations. Together, these developments constitute a computational framework for predictive ab initio quantum many-body calculations in crystalline materials, opening new avenues for understanding and designing strongly correlated quantum materials.

•Mar 5, 2026•0 citations

Quantitative Description of Strongly Correlated Materials by Combining Downfolding Techniques and Tensor Networks

Daan Vrancken, Simon Ganne, Daan Verraes +5 more

High Resolution Image Download MS PowerPoint Slide We present a high-accuracy procedure for electronic structure calculations of strongly correlated materials. To address limitations in current electronic structure methods, we employ density functional theory in combination with the constrained random phase approximation to construct an effective multiband Hubbard model, which is subsequently solved using tensor networks. Our work focuses on one-dimensional and quasi-one-dimensional materials, for which we employ the machinery of matrix product states. We apply this framework to the conjugated polymers trans -polyacetylene and polythiophene, as well as the quasi-one-dimensional charge-transfer insulator Sr 2 CuO 3 . The predicted band gaps show quantitative agreement with state-of-the-art computational techniques and experimental measurements. Beyond band gaps, tensor networks provide access to a wide range of physically relevant properties, including spin magnetization and various excitation energies. Their flexibility supports the implementation of complex Hamiltonians with longer-range interactions, while the bond dimension enables systematic control over accuracy. Furthermore, the computational cost scales efficiently with system size, demonstrating the framework’s scalability.

arXiv:2502.19588•Journal of Chemical Theory and Computation•Aug 10, 2025•0 citations

Downfolding from ab initio to interacting model Hamiltonians: comprehensive analysis and benchmarking of the DFT+cRPA approach

Yueqing Chang, Erik G. C. P. van Loon, Brandon Eskridge +8 more

Abstract Model Hamiltonians are regularly derived from first principles to describe correlated matter. However, the standard methods for this contain a number of largely unexplored approximations. For a strongly correlated impurity model system, here we carefully compare a standard downfolding technique with the best possible ground-truth estimates for charge-neutral excited-state energies and wave functions using state-of-the-art first-principles many-body wave function approaches. To this end, we use the vanadocene molecule and analyze all downfolding aspects, including the Hamiltonian form, target basis, double-counting correction, and Coulomb interaction screening models. We find that the choice of target-space basis functions emerges as a key factor for the quality of the downfolded results, while orbital-dependent double-counting corrections diminish the quality. Background screening of the Coulomb interaction matrix elements primarily affects crystal-field excitations. Our benchmark uncovers the relative importance of each downfolding step and offers insights into the potential accuracy of minimal downfolded model Hamiltonians.

arXiv:2311.05987•npj Computational Materials•Jun 21, 2024•28 citations

Requirements for building effective Hamiltonians using quantum-enhanced density matrix downfolding

Shivesh Pathak, Antonio Russo, Stefan Seritan +4 more

Density matrix downfolding (DMD) is a technique for regressing low-energy effective Hamiltonians from quantum many-body Hamiltonians. One limiting factor in the accuracy of classical implementations of DMD is the presence of difficult-to-quantify systematic errors attendant to sampling the observables of quantum many-body systems on an approximate low-energy subspace. We propose a hybrid quantum-classical protocol for circumventing this limitation, relying on the prospective ability of quantum computers to efficiently prepare and sample from states in well-defined low-energy subspaces with systematically improvable accuracy. We introduce three requirements for when this is possible, including a notion of compressibility that quantifies features of Hamiltonians and low-energy subspaces thereof for which quantum DMD might be efficient. Assuming that these requirements are met, we analyze design choices for our protocol and provide resource estimates for implementing quantum-enhanced DMD on both the doped 2-D Fermi-Hubbard model and an ab initio model of a cuprate superconductor.

arXiv:2403.01043•Mar 2, 2024•1 citation

Ab initio extended Hubbard model of short polyenes for efficient quantum computing

Yuichiro Yoshida, Nayuta Takemori, Wataru Mizukami

We propose introducing an extended Hubbard Hamiltonian derived via the ab initio downfolding method, which was originally formulated for periodic materials, toward efficient quantum computing of molecular electronic structure calculations. By utilizing this method, the first-principles Hamiltonian of chemical systems can be coarse-grained by eliminating the electronic degrees of freedom in higher energy space and reducing the number of terms of electron repulsion integral from O(N4) to O(N2). Our approach is validated numerically on the vertical excitation energies and excitation characters of ethylene, butadiene, and hexatriene. The dynamical electron correlation is incorporated within the framework of the constrained random phase approximation in advance of quantum computations, and the constructed models capture the trend of experimental and high-level quantum chemical calculation results. As expected, the L1-norm of the fermion-to-qubit mapped model Hamiltonians is significantly lower than that of conventional ab initio Hamiltonians, suggesting improved scalability of quantum computing. Those numerical outcomes and the results of the simulation of excited-state sampling demonstrate that the ab initio extended Hubbard Hamiltonian may hold significant potential for quantum chemical calculations using quantum computers.

arXiv:2404.01623•The Journal of Chemical Physics•Aug 28, 2024•6 citations

Quantum Algorithms to Determine Spin-Resolved Exchange-Correlation Potential for Strongly Correlated Materials

H. Arslan Hashim, Volodymyr Turkowski, Eduardo R. Mucciolo

Accurate exchange-correlation (XC) potentials are essential for density functional theory, yet reliable approximations remain challenging for strongly correlated systems. In this work, we present a quantum algorithmic framework to determine spin-resolved XC potentials using a variational quantum eigensolver. Using the Hubbard model as a prototypical strongly correlated lattice system, we prepare ground states in fixed spin sectors through a Hamiltonian variational ansatz initialized from the non-interacting $(U=0)$ Slater-determinant ground state. From the resulting many-body ground states, we extract the XC energy and compute the corresponding spin-resolved XC potentials via finite differences. The accuracy of the approach is benchmarked against exact diagonalization for one- and two-dimensional Hubbard systems of various lattice sizes. We demonstrate that the variational ansatz reproduces the ground-state energies and densities with high fidelity, enabling accurate construction of both magnetic and non-magnetic XC potentials. We analyzed the dependence of the XC potentials on the interaction strength, charge, spin densities, and magnetization. We also present an empirical complexity scaling relation for the computational cost of the method at a fixed fidelity. These results illustrate how quantum computing can be used to construct spin-resolved XC functionals for correlated lattice models, providing a potential pathway for improving density functional approximations in strongly correlated materials.

arXiv:2603.16605•Mar 17, 2026•0 citations
The accurate first-principles description of strongly correlated materials is an important and challenging problem in condensed matter physics. Ab initio downfolding has emerged as a way of deriving compressed many-body Hamiltonians that maintain the essential physics of strongly correlated materials. The solution of these material-specific models is still exponentially difficult to generate on classical computers, but quantum algorithms allow for a significant speed-up in obtaining the ground states of these compressed Hamiltonians. Here, we demonstrate that using quantum algorithms to obtain the properties of downfolded Hamiltonians can indeed yield high-fidelity solutions. By combining ab initio downfolding and variational quantum eigensolvers, we correctly predict the antiferromagnetic state of one-dimensional cuprate Ca2CuO3{\mathrm{Ca}}_{2}{\mathrm{Cu}\mathrm{O}}_{3}Ca2​CuO3​, the excitonic ground state of monolayer WTe2{\mathrm{W}\mathrm{Te}}_{2}WTe2​, and the charge-ordered state of correlated metal SrVO3{\mathrm{Sr}\mathrm{VO}}_{3}SrVO3​. Numerical simulations using a classical tensor network implementation of variational quantum eigensolvers allow us to simulate large models with up to 54 qubits and encompassing up to four bands in the correlated subspace, which is indicative of the complexity that our framework can address. Through these methods we demonstrate the potential of classical preoptimization and downfolding techniques for enabling efficient materials simulation using quantum algorithms.
Achieving materials-specific predictions for large, realistic molecules and solids with strong electron correlations remains a long-standing challenge in quantum chemistry. This dissertation develops scalable ab initio quantum many-body methods for crystalline materials, demonstrating that rigorous optimization of theoretical approximations and numerical algorithms enables predictive computation even in strongly correlated systems. The work tackles three interconnected challenges, namely describing material-specific properties of strongly correlated superconductors, achieving linear scaling with k-points for large-scale periodic calculations, and establishing reliable methods for electron-phonon interactions. The following chapters address each challenge in turn. Chapter 2 demonstrates that ab initio density matrix embedding theory, combined with symmetry-breaking strategies for superconducting order parameters, reproduces experimental trends in cuprate superconductivity, including pressure and layer effects on transition temperatures. By directly solving the electronic Schrödinger equation from material structures, these calculations capture multi-orbital covalency effects that simplified model Hamiltonians neglect, revealing their essential role in the spin fluctuations driving pairing. This predictive capability establishes a foundation for efficient computational screening of superconducting materials. Extending these calculations to larger k-point meshes and broader phase diagram explorations, however, demands more efficient computational algorithms. To overcome this computational bottleneck, Chapter 3 introduces interpolative separable density fitting framework. This development enables thermodynamic limit calculations with up to 1000 k-points, as validated on diamond, carbon dioxide, nickel monoxide, and cuprate systems. Integration with density matrix embedding theory and local natural orbital methods further yields converged ground-state energies from accurate correlated wavefunction methods such as CCSD(T) for crystalline materials in the thermodynamic limit. With computational efficiency established, the remaining challenge lies in extending the theoretical framework beyond purely electronic degrees of freedom. Chapter 4 addresses this need by establishing exponential ansatz wavefunctions as competitive approaches for electron-phonon systems. Unlike previous treatments based on perturbation theory, systematic benchmarking of coupled cluster theory and variational Lang-Firsov methods on the Holstein model demonstrates, for the first time, that these approaches accurately describe polaron formation across coupling regimes when combined with appropriate reference state optimization. These results provide a foundation for incorporating phonon effects into ab initio studies of superconducting materials. In addition to these methodological advances, GPU acceleration and MPI parallelization enhance computational throughput, enabling efficient large-scale calculations. Together, these developments constitute a computational framework for predictive ab initio quantum many-body calculations in crystalline materials, opening new avenues for understanding and designing strongly correlated quantum materials.
High Resolution Image Download MS PowerPoint Slide We present a high-accuracy procedure for electronic structure calculations of strongly correlated materials. To address limitations in current electronic structure methods, we employ density functional theory in combination with the constrained random phase approximation to construct an effective multiband Hubbard model, which is subsequently solved using tensor networks. Our work focuses on one-dimensional and quasi-one-dimensional materials, for which we employ the machinery of matrix product states. We apply this framework to the conjugated polymers trans -polyacetylene and polythiophene, as well as the quasi-one-dimensional charge-transfer insulator Sr 2 CuO 3 . The predicted band gaps show quantitative agreement with state-of-the-art computational techniques and experimental measurements. Beyond band gaps, tensor networks provide access to a wide range of physically relevant properties, including spin magnetization and various excitation energies. Their flexibility supports the implementation of complex Hamiltonians with longer-range interactions, while the bond dimension enables systematic control over accuracy. Furthermore, the computational cost scales efficiently with system size, demonstrating the framework’s scalability.
Abstract Model Hamiltonians are regularly derived from first principles to describe correlated matter. However, the standard methods for this contain a number of largely unexplored approximations. For a strongly correlated impurity model system, here we carefully compare a standard downfolding technique with the best possible ground-truth estimates for charge-neutral excited-state energies and wave functions using state-of-the-art first-principles many-body wave function approaches. To this end, we use the vanadocene molecule and analyze all downfolding aspects, including the Hamiltonian form, target basis, double-counting correction, and Coulomb interaction screening models. We find that the choice of target-space basis functions emerges as a key factor for the quality of the downfolded results, while orbital-dependent double-counting corrections diminish the quality. Background screening of the Coulomb interaction matrix elements primarily affects crystal-field excitations. Our benchmark uncovers the relative importance of each downfolding step and offers insights into the potential accuracy of minimal downfolded model Hamiltonians.
Density matrix downfolding (DMD) is a technique for regressing low-energy effective Hamiltonians from quantum many-body Hamiltonians. One limiting factor in the accuracy of classical implementations of DMD is the presence of difficult-to-quantify systematic errors attendant to sampling the observables of quantum many-body systems on an approximate low-energy subspace. We propose a hybrid quantum-classical protocol for circumventing this limitation, relying on the prospective ability of quantum computers to efficiently prepare and sample from states in well-defined low-energy subspaces with systematically improvable accuracy. We introduce three requirements for when this is possible, including a notion of compressibility that quantifies features of Hamiltonians and low-energy subspaces thereof for which quantum DMD might be efficient. Assuming that these requirements are met, we analyze design choices for our protocol and provide resource estimates for implementing quantum-enhanced DMD on both the doped 2-D Fermi-Hubbard model and an ab initio model of a cuprate superconductor.
We propose introducing an extended Hubbard Hamiltonian derived via the ab initio downfolding method, which was originally formulated for periodic materials, toward efficient quantum computing of molecular electronic structure calculations. By utilizing this method, the first-principles Hamiltonian of chemical systems can be coarse-grained by eliminating the electronic degrees of freedom in higher energy space and reducing the number of terms of electron repulsion integral from O(N4) to O(N2). Our approach is validated numerically on the vertical excitation energies and excitation characters of ethylene, butadiene, and hexatriene. The dynamical electron correlation is incorporated within the framework of the constrained random phase approximation in advance of quantum computations, and the constructed models capture the trend of experimental and high-level quantum chemical calculation results. As expected, the L1-norm of the fermion-to-qubit mapped model Hamiltonians is significantly lower than that of conventional ab initio Hamiltonians, suggesting improved scalability of quantum computing. Those numerical outcomes and the results of the simulation of excited-state sampling demonstrate that the ab initio extended Hubbard Hamiltonian may hold significant potential for quantum chemical calculations using quantum computers.
Accurate exchange-correlation (XC) potentials are essential for density functional theory, yet reliable approximations remain challenging for strongly correlated systems. In this work, we present a quantum algorithmic framework to determine spin-resolved XC potentials using a variational quantum eigensolver. Using the Hubbard model as a prototypical strongly correlated lattice system, we prepare ground states in fixed spin sectors through a Hamiltonian variational ansatz initialized from the non-interacting (U=0)(U=0)(U=0) Slater-determinant ground state. From the resulting many-body ground states, we extract the XC energy and compute the corresponding spin-resolved XC potentials via finite differences. The accuracy of the approach is benchmarked against exact diagonalization for one- and two-dimensional Hubbard systems of various lattice sizes. We demonstrate that the variational ansatz reproduces the ground-state energies and densities with high fidelity, enabling accurate construction of both magnetic and non-magnetic XC potentials. We analyzed the dependence of the XC potentials on the interaction strength, charge, spin densities, and magnetization. We also present an empirical complexity scaling relation for the computational cost of the method at a fixed fidelity. These results illustrate how quantum computing can be used to construct spin-resolved XC functionals for correlated lattice models, providing a potential pathway for improving density functional approximations in strongly correlated materials.