pith. sign in

arxiv: 2602.10435 · v2 · submitted 2026-02-11 · 💻 cs.CE

Multiconfiguration Pair-Density Functional Theory Calculations of Low-lying States of Complex Chemical Systems with Quantum Computers

Zhanou Liu , Yuhao Chen , Yingjin Ma , Xiao He , Yuxin Deng This is my paper

Pith reviewed 2026-05-16 03:57 UTC · model grok-4.3

classification 💻 cs.CE
keywords quantum computingmulticonfiguration pair-density functional theoryvariational quantum eigensolverstrong electron correlationchromium dimerchemical accuracynear-term quantum hardware
0
0 comments X

The pith

A hybrid quantum-classical method separates static and dynamic electron correlation to compute accurate energies for strongly correlated molecules like Cr2 on noisy near-term hardware.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes that confining static correlation to a compact multireference quantum state computed by the variational quantum eigensolver, while recovering dynamic correlation via a classical on-top density functional, produces reliable results for complex systems. This decoupling enables self-consistent orbital optimization and lowers the quantum circuit depth needed. On standard test cases the approach reproduces C2 bond lengths and benzene excitation energies to chemical accuracy. For the demanding Cr2 dimer with a 48-electron 42-orbital active space it yields a bound potential-energy curve and the correct qualitative dissociation profile even when realistic hardware noise is included. A reader would care because the strategy makes predictions on molecules too strongly correlated for classical single-reference methods feasible without waiting for fault-tolerant quantum computers.

Core claim

The hybrid variational quantum eigensolver plus multiconfiguration pair-density functional theory framework decouples electron correlation by treating static correlation exactly within the quantum multireference wave function and dynamic correlation through a classical functional evaluated on the reduced density supplied by that wave function; the resulting method recovers chemical accuracy on C2 and benzene benchmarks and produces a bound potential-energy curve with correct dissociation behavior for Cr2 despite hardware noise.

What carries the argument

Multiconfiguration pair-density functional theory, in which an on-top density functional is applied to the one- and two-body reduced densities obtained from a quantum-computed multireference state to capture dynamic correlation after static correlation has been handled by the variational quantum eigensolver.

If this is right

  • C2 equilibrium bond lengths are reproduced with a mean absolute error of 0.006 Å.
  • Benzene low-lying excitation energies are obtained with a mean absolute error of 0.048 eV.
  • A 48-electron, 42-orbital active space for Cr2 becomes tractable on near-term hardware while still yielding a bound curve.
  • Quantum resource overhead is reduced because only static correlation is retained in the variational quantum eigensolver circuit.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same separation could be applied to other open-shell transition-metal complexes whose active spaces exceed classical full-configuration-interaction limits.
  • Self-consistent orbital relaxation performed classically after each quantum iteration may systematically improve convergence for active spaces larger than those demonstrated here.
  • If the on-top functional remains reliable across a broader range of molecular geometries, the method provides a concrete route to thermochemistry and spectroscopy on systems currently inaccessible to both classical multireference and pure quantum algorithms.

Load-bearing premise

The classical on-top density functional must recover dynamic correlation accurately from the reduced-density data supplied by the quantum multireference state, and self-consistent orbital optimization must not introduce large uncontrolled errors.

What would settle it

If the computed Cr2 potential-energy curve shows no binding minimum or fails to recover the expected qualitative shape at large separation when the on-top functional is applied to the quantum-computed densities under realistic noise levels, the separation-of-correlations premise is falsified.

read the original abstract

Accurately describing strong electron correlation in complex systems remains a prominent challenge in computational chemistry as near-term quantum algorithms treating total correlation often require prohibitively deep circuits. Here we present a hybrid strategy combining the Variational Quantum Eigensolver with Multiconfiguration Pair-Density Functional Theory to efficiently decouple correlation effects. This approach confines static correlation to a compact multireference quantum state while recovering dynamic correlation through a classical on-top density functional using reduced-density information. By enabling self-consistent orbital optimization, the method significantly reduces quantum resource overheads without sacrificing physical rigor. We demonstrate chemical accuracy on standard benchmarks by reproducing C$_2$ equilibrium bond lengths and benzene excitation energies with mean absolute errors of 0.006 {\AA} and 0.048 eV respectively. Most notably, for the strongly correlated Cr$_2$ dimer requiring a large complete active space (48e, 42o), the framework yields a bound potential-energy curve and recovers qualitative dissociation behavior despite realistic hardware noise. These results establish that separating correlation types provides a practical route to reliable predictions on near-term quantum hardware.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper introduces a hybrid VQE-MC-PDFT framework that confines static correlation to a quantum multireference state in a compact active space while recovering dynamic correlation classically via an on-top pair-density functional evaluated on reduced density matrices extracted from the VQE wavefunction. Self-consistent orbital optimization is used to lower quantum resource costs. Demonstrations include chemical accuracy on C2 bond lengths (MAE 0.006 Å) and benzene excitation energies (MAE 0.048 eV), plus a bound potential-energy curve with qualitative dissociation for the strongly correlated Cr2 dimer in a (48e,42o) CAS despite hardware noise.

Significance. If the central claims hold after additional validation, the work provides a concrete route to treat larger strongly correlated systems on near-term quantum hardware by cleanly separating correlation types, with the Cr2 result illustrating feasibility for active spaces that exceed typical VQE limits. The approach builds on established MC-PDFT components without introducing new fitted parameters.

major comments (2)
  1. [Abstract] Abstract: The reported MAEs for C2 equilibrium bond lengths and benzene excitations are given without error bars, number of independent VQE runs, shot statistics, or direct baseline comparisons to pure VQE or classical MC-PDFT, which are required to establish that the hybrid separation improves upon either component alone.
  2. [Results (Cr2)] Cr2 results section: The claim that the framework yields a bound PEC and recovers qualitative dissociation for the (48e,42o) CAS relies on the classical on-top functional accurately recovering dynamic correlation from noisy quantum 1-RDM and on-top 2-RDM inputs; however, no error propagation from finite shots/hardware noise, no comparison of quantum versus exact CAS densities, and no sensitivity tests on the functional input are provided, leaving this separation step unverified.
minor comments (2)
  1. The abstract and methods should explicitly name the on-top density functional employed and the form of the self-consistent orbital optimization procedure.
  2. Notation for the reduced density matrices supplied to the functional should be clarified with explicit definitions of the 1-RDM and on-top 2-RDM extraction steps.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for their constructive and detailed comments, which have helped us identify areas for improvement in clarity and validation. We address each major comment point by point below, indicating where revisions will be made to the manuscript.

read point-by-point responses

  1. Referee: [Abstract] Abstract: The reported MAEs for C2 equilibrium bond lengths and benzene excitations are given without error bars, number of independent VQE runs, shot statistics, or direct baseline comparisons to pure VQE or classical MC-PDFT, which are required to establish that the hybrid separation improves upon either component alone.

    Authors: We agree that including error bars, the number of independent VQE runs, shot statistics, and explicit baseline comparisons would strengthen the claims. In the revised manuscript we will update the abstract and main text to report these statistical details (obtained from multiple independent runs with the specified shot counts) and add direct comparisons to pure VQE and classical MC-PDFT results on the same systems to demonstrate the improvement from the hybrid separation of static and dynamic correlation. revision: yes

  2. Referee: [Results (Cr2)] Cr2 results section: The claim that the framework yields a bound PEC and recovers qualitative dissociation for the (48e,42o) CAS relies on the classical on-top functional accurately recovering dynamic correlation from noisy quantum 1-RDM and on-top 2-RDM inputs; however, no error propagation from finite shots/hardware noise, no comparison of quantum versus exact CAS densities, and no sensitivity tests on the functional input are provided, leaving this separation step unverified.

    Authors: We acknowledge the value of additional validation for the correlation separation step. In the revised manuscript we will include an error-propagation analysis from finite-shot noise on the extracted 1-RDM and on-top 2-RDM, as well as sensitivity tests varying the on-top functional parameters. However, a direct comparison of the quantum-derived densities against exact CAS densities is not feasible for the (48e,42o) active space, as the classical full-CI calculation is intractable; this limitation is inherent to the problem size and is precisely the motivation for the quantum approach. We will explicitly note this constraint while retaining the qualitative demonstration of a bound curve under realistic hardware noise. revision: partial

standing simulated objections not resolved
  • Direct comparison of quantum versus exact CAS densities for the (48e,42o) Cr2 active space

Circularity Check

0 steps flagged

No significant circularity; hybrid VQE-MC-PDFT separation is independently grounded

full rationale

The paper's core derivation applies standard VQE to obtain a multireference state for static correlation and evaluates a classical on-top pair-density functional on the resulting reduced densities to recover dynamic correlation. No equations reduce a reported prediction to a quantity defined by the paper's own fitted parameters or self-citations. The MC-PDFT component is treated as an established external method whose validity is not derived within this work, and the Cr2 demonstration is presented as an empirical test of the hybrid workflow rather than a tautological restatement of inputs. Self-citations to foundational MC-PDFT literature are not load-bearing for the central claims.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the domain assumption that static and dynamic correlation can be cleanly decoupled and that the on-top functional remains accurate when fed reduced-density data from a quantum multireference state; no free parameters or invented entities are explicitly introduced in the abstract.

axioms (1)
  • domain assumption Standard quantum mechanics and the validity of the multiconfiguration pair-density functional approximation for recovering dynamic correlation
    Invoked when the method confines static correlation to the quantum state and recovers dynamic correlation classically via the on-top functional.

pith-pipeline@v0.9.0 · 5504 in / 1254 out tokens · 89836 ms · 2026-05-16T03:57:42.148554+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

50 extracted references · 50 canonical work pages

  1. [1]

    Hartree, D. R. The wave mechanics of an atom with a non-Coulomb central field. part i. theory and methods. Math. Proc. Camb. Phil. Soc.24, 89–110 (1928)

    work page 1928

  2. [2]

    Slater, J. C. The theory of complex spectra.Phys. Rev. 34, 1293 (1929)

    work page 1929

  3. [3]

    Slater, J. C. Note on hartree’s method.Phys. Rev.35, 210 (1930)

    work page 1930

  4. [4]

    Roothaan, C. C. J. New developments in molecular orbital theory.Rev. Mod. Phys.23, 69 (1951)

    work page 1951

  5. [5]

    K¨ ummel, H. G. A biography of the coupled cluster method.Int. J. Mod. Phys. B17, 5311–5325 (2003)

    work page 2003

  6. [6]

    Bartlett, R. J. & Musia l, M. Coupled-cluster theory in quantum chemistry.Rev. Mod. Phys.79, 291–352 (2007)

    work page 2007

  7. [7]

    Møller–plesset perturbation theory: from small molecule methods to methods for thousands of atoms.WIREs Comput

    Cremer, D. Møller–plesset perturbation theory: from small molecule methods to methods for thousands of atoms.WIREs Comput. Mol. Sci.1, 509–530 (2011)

    work page 2011

  8. [8]

    & Kohn, W

    Hohenberg, P. & Kohn, W. Inhomogeneous electron gas. Phys. Rev.136, B864 (1964)

    work page 1964

  9. [9]

    Universal variational functionals of electron densities, first-order density matrices, and natural spin- orbitals and solution of the v-representability problem

    Levy, M. Universal variational functionals of electron densities, first-order density matrices, and natural spin- orbitals and solution of the v-representability problem. Proc. Natl Acad. Sci. U.S.A.76, 6062–6065 (1979)

    work page 1979

  10. [10]

    & Rasolt, M

    Vignale, G. & Rasolt, M. Density-functional theory in strong magnetic fields.Phys. Rev. Lett.59, 2360 (1987)

    work page 1987

  11. [11]

    Nielsen, M. A. & Chuang, I. L.Quantum Computation and Quantum Information: 10th Anniversary Edition (Cambridge University Press, Cambridge, 2010)

    work page 2010

  12. [12]

    Hartree-fock on a superconducting qubit quantum computer.Science 369, 1084–1089 (2020)

    Google AI Quantum and Collaborators. Hartree-fock on a superconducting qubit quantum computer.Science 369, 1084–1089 (2020)

    work page 2020

  13. [13]

    J., Head-Gordon, M

    Lee, J., Huggins, W. J., Head-Gordon, M. & Wha- ley, K. B. Generalized unitary coupled cluster wave functions for quantum computation.J. Chem. Theory Comput.15, 311–324 (2018)

    work page 2018

  14. [14]

    M., Hohenstein, E

    Parrish, R. M., Hohenstein, E. G., McMahon, P. L. & Mart´ ınez, T. J. Quantum computation of electronic transitions using a variational quantum eigensolver. Phys. Rev. Lett.122, 230401 (2019)

    work page 2019

  15. [15]

    & Mu˜ noz Ramo, D

    Greene-Diniz, G. & Mu˜ noz Ramo, D. Generalized unitary coupled cluster excitations for multireference molecular states optimized by the variational quantum eigensolver.Int. J. Quantum Chem.121, e26352 (2021)

    work page 2021

  16. [16]

    Mizukami, W.et al.Orbital optimized unitary coupled cluster theory for quantum computer.Phys. Rev. Res. 2, 033421 (2020)

    work page 2020

  17. [17]

    O.et al.Quantum orbital-optimized uni- tary coupled cluster methods in the strongly correlated regime: Can quantum algorithms outperform their clas- sical equivalents?J

    Sokolov, I. O.et al.Quantum orbital-optimized uni- tary coupled cluster methods in the strongly correlated regime: Can quantum algorithms outperform their clas- sical equivalents?J. Chem. Phys.152, 124107 (2020)

    work page 2020

  18. [18]

    Takeshita, T.et al.Increasing the representation accu- racy of quantum simulations of chemistry without extra quantum resources.Phys. Rev. X10, 011004 (2020)

    work page 2020

  19. [19]

    Technol.6, 024004 (2021)

    Yalouz, S.et al.A state-averaged orbital-optimized hybrid quantum–classical algorithm for a democratic description of ground and excited states.Quantum Sci. Technol.6, 024004 (2021)

    work page 2021

  20. [20]

    & Mochizuki, Y

    Sugisaki, K., Kato, T., Minato, Y., Okuwaki, K. & Mochizuki, Y. Variational quantum eigensolver simu- lations with the multireference unitary coupled cluster ansatz: a case study of the c 2v quasi-reaction pathway of beryllium insertion into a h 2 molecule.Phys. Chem. Chem. Phys.24, 8439–8452 (2022)

    work page 2022

  21. [21]

    Khan, I.et al.Chemically aware unitary coupled cluster with ab initio calculations on an ion trap quantum com- puter: A refrigerant chemicals’ application.J. Chem. Phys.158, 214114 (2023)

    work page 2023

  22. [22]

    Fitzpatrick, A.et al.Self-consistent field approach for the variational quantum eigensolver: orbital optimiza- tion goes adaptive.J. Phys. Chem. A128, 2843–2856 (2024)

    work page 2024

  23. [23]

    Zhao, S.et al.Quantum computation of conical inter- sections on a programmable superconducting quantum 11 processor.J. Phys. Chem. Lett.15, 7244–7253 (2024)

    work page 2024

  24. [24]

    Roos, B. O. The complete active space self-consistent field method and its applications in electronic structure calculations.Adv. Chem. Phys.69, 399–445 (1987)

    work page 1987

  25. [25]

    Bharti, K.et al.Noisy intermediate-scale quantum algorithms.Rev. Mod. Phys.94, 015004 (2022)

    work page 2022

  26. [26]

    Li Manni, G.et al.Multiconfiguration pair-density func- tional theory.J. Chem. Theory Comput.10, 3669–3680 (2014)

    work page 2014

  27. [27]

    Chemical Science (RSC Publishing) (2026).https://doi.org/10.1039/D5SC07528A

    Chen, Y.et al.Quantum-classical hybrid computa- tion of electron transfer in a cryptochrome protein via vqe-pdft and multiscale modeling.Chem. Sci.DOI: 10.1039/D5SC07528A (2026). In press

    work page doi:10.1039/d5sc07528a 2026

  28. [28]

    Jiang, J.et al.Diabatic valence-hole states in the c2 molecule: ”putting humpty dumpty together again”.J. Phys. Chem. A126, 3090–3100 (2022)

    work page 2022

  29. [29]

    A., Umrigar, C

    Holmes, A. A., Umrigar, C. & Sharma, S. Excited states using semistochastic heat-bath configuration interac- tion.J. Chem. Phys.147, 164111 (2017)

    work page 2017

  30. [30]

    V´ eril, M.et al.Questdb: A database of highly accurate excitation energies for the electronic structure commu- nity.WIREs Comput. Mol. Sci.11, e1517 (2021)

    work page 2021

  31. [31]

    C2 spectroscopy and kinetics.J

    Martin, M. C2 spectroscopy and kinetics.J. Pho- tochem. Photobiol. A: Chem.66, 263–289 (1992)

    work page 1992

  32. [32]

    R., Zhai, H., Umrigar, C

    Larsson, H. R., Zhai, H., Umrigar, C. J. & Chan, G. K.- L. The chromium dimer: closing a chapter of quantum chemistry.J. Am. Chem. Soc.144, 15932–15937 (2022)

    work page 2022

  33. [33]

    & Yanai, T

    Kurashige, Y. & Yanai, T. Second-order perturbation theory with a density matrix renormalization group self- consistent field reference function: Theory and applica- tion to the study of chromium dimer.J. Chem. Phys. 135, 094104 (2011)

    work page 2011

  34. [34]

    Bauschlicher Jr, C. W. & Partridge, H. Cr2 revisited. Chem. Phys. Lett.231, 277–282 (1994)

    work page 1994

  35. [35]

    Li, J.et al.Accurate many-body electronic structure near the basis set limit: Application to the chromium dimer.Phys. Rev. Res.2, 012015 (2020)

    work page 2020

  36. [36]

    & Alavi, A

    Sharma, S. & Alavi, A. Multireference linearized cou- pled cluster theory for strongly correlated systems using matrix product states.J. Chem. Phys.143, 102815 (2015)

    work page 2015

  37. [37]

    & Jacquemin, D

    Loos, P.-F., Lipparini, F., Boggio-Pasqua, M., Sce- mama, A. & Jacquemin, D. A mountaineering strategy to excited states: Highly accurate energies and bench- marks for medium sized molecules.J. Chem. Theory Comput.16, 1711–1741 (2020)

    work page 2020

  38. [38]

    & Spencer, J

    Pfau, D., Axelrod, S., Sutterud, H., von Glehn, I. & Spencer, J. S. Accurate computation of quantum excited states with neural networks.Science385, eadn0137 (2024)

    work page 2024

  39. [39]

    Phys.20, 1240–1246 (2024)

    Guo, S.et al.Experimental quantum computational chemistry with optimized unitary coupled cluster ansatz.Nat. Phys.20, 1240–1246 (2024)

    work page 2024

  40. [40]

    & Olsen, J.Molecu- lar Electronic-Structure Theory(John Wiley & Sons, Chichester, 2014)

    Helgaker, T., Jorgensen, P. & Olsen, J.Molecu- lar Electronic-Structure Theory(John Wiley & Sons, Chichester, 2014)

    work page 2014

  41. [41]

    & Wigner, E

    Jordan, P. & Wigner, E. ¨Uber das paulische ¨ aquivalenzverbot.Z. Phys.47, 631–651 (1928)

    work page 1928

  42. [42]

    T., Richard, M

    Seeley, J. T., Richard, M. J. & Love, P. J. The bravyi-kitaev transformation for quantum computation of electronic structure.J. Chem. Phys.137, 224109 (2012)

    work page 2012

  43. [43]

    Bravyi, S. B. & Kitaev, A. Y. Fermionic quantum computation.Ann. Phys.298, 210–226 (2002)

    work page 2002

  44. [44]

    Commun.5, 4213 (2014)

    Peruzzo, A.et al.A variational eigenvalue solver on a photonic quantum processor.Nat. Commun.5, 4213 (2014)

    work page 2014

  45. [45]

    R., Romero, J., Babbush, R

    McClean, J. R., Romero, J., Babbush, R. & Aspuru- Guzik, A. The theory of variational hybrid quantum- classical algorithms.New J. Phys.18, 023023 (2016)

    work page 2016

  46. [46]

    Generalized trotter’s formula and system- atic approximants of exponential operators and inner derivations with applications to many-body problems

    Suzuki, M. Generalized trotter’s formula and system- atic approximants of exponential operators and inner derivations with applications to many-body problems. Commun. Math. Phys.51, 183–190 (1976)

    work page 1976

  47. [47]

    B., Zimbor´ as, Z

    Maciejewski, F. B., Zimbor´ as, Z. & Oszmaniec, M. Miti- gation of readout noise in near-term quantum devices by classical post-processing based on detector tomography. Quantum4, 257 (2020)

    work page 2020

  48. [48]

    Qufem: Fast and accurate quantum read- out calibration using the finite element method, 948–963 (Association for Computing Machinery, La Jolla, CA, USA, 2024)

    Tan, S.et al. Qufem: Fast and accurate quantum read- out calibration using the finite element method, 948–963 (Association for Computing Machinery, La Jolla, CA, USA, 2024)

    work page 2024

  49. [49]

    & Zeng, W

    Giurgica-Tiron, T., Hindy, Y., LaRose, R., Mari, A. & Zeng, W. J.Digital zero noise extrapolation for quantum error mitigation, 306–316 (IEEE, Denver, CO, USA, 2020)

    work page 2020

  50. [50]

    Digital Silk Road

    Czarnik, P., Arrasmith, A., Coles, P. J. & Cincio, L. Error mitigation with clifford quantum-circuit data. Quantum5, 592 (2021). Acknowledgments Y.D. was supported by the National Key R&D Program of China (Grant No. 2023YFA1009403), the National Nat- ural Science Foundation of China (Grant No. 62472175), the “Digital Silk Road” Shanghai International Join...

    work page 2021