Recognition: 2 theorem links
· Lean TheoremExponential Scaling Barriers for Variational Quantum Eigensolvers
Pith reviewed 2026-05-15 11:54 UTC · model grok-4.3
The pith
Adaptive VQE requires an exponentially growing number of iterations as the number of molecular orbitals increases.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We demonstrate that the Rényi entropy derived from classical simulations predicts the required number of adaptive iterations of VQE with high accuracy (R² ≈ 0.99). We validate this on a benchmarking set of more than 20 different molecules with active spaces ranging from four to ten orbitals. For these molecules, we find an exponential scaling of the number of adaptive iterations, and in turn, of the circuit depth with the system size. We therefore conclude that it is unlikely that VQE in its current form is able to simulate large molecular systems with high fidelity without exponential resource requirements.
What carries the argument
The Rényi entropy computed on the classical wave function, which serves as a predictor of the number of adaptive VQE iterations needed to reach a given fidelity.
If this is right
- Circuit depth in adaptive VQE grows exponentially with the number of orbitals.
- The total quantum resources required scale exponentially rather than polynomially for larger molecular targets.
- High-fidelity ground-state calculations on large molecules remain out of reach for VQE without exponential overhead.
- Classical pre-computation of Rényi entropy can forecast the quantum cost before any quantum runs occur.
Where Pith is reading between the lines
- Similar entropy-based diagnostics could be applied to other variational quantum algorithms to forecast their resource scaling.
- The correlation suggests that any ansatz whose expressivity must match the target's entanglement will inherit comparable scaling limits.
- Hybrid classical-quantum workflows might use the entropy predictor to decide when to switch from VQE to a different method.
Load-bearing premise
The exponential growth in iteration count seen for four-to-ten orbital active spaces will continue unchanged once classical simulation of the full system becomes impossible.
What would settle it
Performing adaptive VQE on a molecule whose active space has twenty or more orbitals and measuring whether the iteration count remains polynomial rather than exponential in orbital number.
Figures
read the original abstract
The Variational Quantum Eigensolver (VQE) is widely regarded as a promising algorithm for calculating ground states of quantum systems that are intractable for classical computers. This promise is typically motivated by the hope of mitigating the exponential growth of Hilbert space with system size. Here we scrutinize how the computational cost of adaptive VQE scales with the size of the target system. We demonstrate that the R\'enyi entropy derived from classical simulations predicts the required number of adaptive iterations of VQE with high accuracy ($R^2 \approx 0.99$). We validate this on a benchmarking set of more than 20 different molecules with active spaces ranging from four to ten orbitals. For these molecules, we find an exponential scaling of the number of adaptive iterations, and in turn, of the circuit depth with the system size. We therefore conclude that it is unlikely that VQE in its current form is able to simulate large molecular systems with high fidelity without exponential resource requirements.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports numerical benchmarks of adaptive Variational Quantum Eigensolver (VQE) on more than 20 molecules with active spaces from 4 to 10 orbitals. It finds that the number of adaptive iterations scales exponentially with system size and is highly correlated (R² ≈ 0.99) with the Rényi entropy obtained from classical simulations. Based on this, the authors conclude that VQE in its current form is unlikely to simulate large molecular systems with high fidelity without requiring exponential resources.
Significance. If the observed exponential scaling and its correlation with Rényi entropy persist to larger systems, this would represent a significant finding that challenges the practical utility of VQE for quantum chemistry problems beyond small active spaces. The strength lies in the extensive benchmarking set and the independent computation of the predictor, providing a potential diagnostic tool. However, the significance is tempered by the limited range of system sizes tested.
major comments (3)
- [Results section (around the exponential fit and correlation analysis)] The exponential scaling is fitted to data from active spaces of 4-10 orbitals, but the manuscript does not provide error bars on the fit parameters or on the reported R² ≈ 0.99, making it difficult to assess the statistical robustness of the extrapolation claim.
- [Abstract and conclusion] The strong conclusion that VQE requires exponential resources for large systems rests on the assumption that the Rényi entropy-iteration correlation continues beyond the classically simulable regime; no scaling argument or additional evidence is given to support this continuity.
- [Methods or benchmarking set description] The criteria for selecting the more than 20 molecules and their active spaces are not detailed, which is necessary to evaluate whether the exponential scaling is general or specific to the chosen set.
minor comments (1)
- [Abstract] The abstract mentions 'high accuracy (R² ≈ 0.99)' but does not specify if this is for the correlation or the fit; clarify the exact quantity.
Simulated Author's Rebuttal
We appreciate the referee's detailed review and constructive suggestions. Below we provide point-by-point responses to the major comments, indicating where revisions will be made to the manuscript.
read point-by-point responses
-
Referee: The exponential scaling is fitted to data from active spaces of 4-10 orbitals, but the manuscript does not provide error bars on the fit parameters or on the reported R² ≈ 0.99, making it difficult to assess the statistical robustness of the extrapolation claim.
Authors: We agree with this observation. In the revised version, we will compute and report error bars on the exponential fit parameters using bootstrap resampling and also provide the uncertainty on the R² value to better assess the robustness of our claims. revision: yes
-
Referee: The strong conclusion that VQE requires exponential resources for large systems rests on the assumption that the Rényi entropy-iteration correlation continues beyond the classically simulable regime; no scaling argument or additional evidence is given to support this continuity.
Authors: We acknowledge that our conclusions involve an extrapolation. While we cannot provide direct evidence for systems beyond classical simulability, the strong correlation with Rényi entropy, which is computable classically, provides a basis for the prediction. We will revise the abstract and conclusion to explicitly note this extrapolation and its basis in the observed data. revision: partial
-
Referee: The criteria for selecting the more than 20 molecules and their active spaces are not detailed, which is necessary to evaluate whether the exponential scaling is general or specific to the chosen set.
Authors: We will expand the Methods section to detail the criteria used for selecting the molecules and their active spaces, including considerations for diversity in molecular types and active space sizes to support the generality of the results. revision: yes
Circularity Check
No significant circularity; empirical correlation and explicit extrapolation remain independent of inputs.
full rationale
The paper computes Rényi entropy via separate classical simulations on the same small active spaces (4-10 orbitals) and reports an observed high-accuracy correlation (R² ≈ 0.99) with the measured number of adaptive VQE iterations across >20 molecules. This correlation is presented as an empirical validation rather than a definitional identity or a parameter fitted directly to the iteration counts and then relabeled as a prediction. The exponential scaling is fitted to the observed iteration counts versus system size, and the final claim about large-molecule intractability is framed as an extrapolation resting on the untested persistence of the entropy-iteration link outside the classically simulable regime. No load-bearing self-citation, uniqueness theorem, or ansatz smuggling is invoked; the derivation chain consists of direct computation, correlation measurement, and stated extrapolation without reducing any step to its own inputs by construction.
Axiom & Free-Parameter Ledger
axioms (1)
- standard math Standard assumptions of quantum mechanics, variational principles, and the definition of Rényi entropy
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We demonstrate that the Rényi entropy derived from classical simulations predicts the required number of adaptive iterations of VQE with high accuracy (R² ≈ 0.99). ... exponential scaling of the number of adaptive iterations
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
For these molecules, we find an exponential scaling of the number of adaptive iterations, and in turn, of the circuit depth with the system size.
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
-
[1]
Y. Cao, J. Romero, J. P. Olson, M. Degroote, P. D. John- son, M. Kieferová, I. D. Kivlichan, T. Menke, B. Per- opadre, N. P. D. Sawaya, S. Sim, L. Veis, and A. Aspuru- Guzik, Quantum chemistry in the age of quantum com- puting, Chemical reviews119, 10856 (2019)
work page 2019
-
[2]
Burke, Perspective on density functional theory, The Journal of chemical physics136, 150901 (2012)
K. Burke, Perspective on density functional theory, The Journal of chemical physics136, 150901 (2012)
work page 2012
-
[3]
R. J. Bartlett and M. Musiał, Coupled-cluster theory in quantum chemistry, Reviews of Modern Physics79, 291 (2007)
work page 2007
-
[4]
S. Lee, J. Lee, H. Zhai, Y. Tong, A. M. Dalzell, A. Ku- mar, P. Helms, J. Gray, Z.-H. Cui, W. Liu, M. Kasto- ryano, R. Babbush, J. Preskill, D. R. Reichman, E. T. Campbell, E. F. Valeev, L. Lin, and G. K.-L. Chan, Eval- uating the evidence for exponential quantum advantage in ground-state quantum chemistry, Nature communica- tions14, 1952 (2023). 12
work page 1952
- [5]
-
[6]
K. D. Vogiatzis, D. Ma, J. Olsen, L. Gagliardi, and W. A. de Jong, Pushing configuration-interaction to the limit: Towards massively parallel mcscf calculations, The Jour- nal of chemical physics147, 184111 (2017)
work page 2017
-
[7]
S. R. White, Density matrix formulation for quantum renormalization groups, Physical review letters69, 2863 (1992)
work page 1992
-
[8]
G. H. Booth, A. J. W. Thom, and A. Alavi, Fermion monte carlo without fixed nodes: a game of life, death, andannihilationinslaterdeterminantspace,TheJournal of chemical physics131, 054106 (2009)
work page 2009
-
[9]
S. Sharma and G. K.-L. Chan, Spin-adapted density matrix renormalization group algorithms for quantum chemistry, The Journal of chemical physics136, 124121 (2012)
work page 2012
-
[10]
A. Aspuru-Guzik, A. D. Dutoi, P. J. Love, and M. Head- Gordon, Simulated quantum computation of molecular energies, Science (New York, N.Y.)309, 1704 (2005)
work page 2005
-
[11]
I. M. Georgescu, S. Ashhab, and F. Nori, Quantum sim- ulation, Reviews of Modern Physics86, 153 (2014)
work page 2014
-
[12]
R. Babbush, D. W. Berry, J. R. McClean, and H. Neven, Quantum simulation of chemistry with sub- linear scaling in basis size, npj Quantum Information5, 10.1038/s41534-019-0199-y (2019)
-
[13]
A. Peruzzo, J. McClean, P. Shadbolt, M.-H. Yung, X.-Q. Zhou, P. J. Love, A. Aspuru-Guzik, and J. L. O’Brien, A variational eigenvalue solver on a photonic quantum processor, Nature communications5, 4213 (2014)
work page 2014
-
[14]
H. R. Grimsley, S. E. Economou, E. Barnes, and N. J. Mayhall, An adaptive variational algorithm for exact molecular simulations on a quantum computer, Nature communications10, 3007 (2019)
work page 2019
-
[15]
A. Y. Kitaev, Quantum measurements and the abelian stabilizer problem (1995), arXiv:quant-ph/9511026 [quant-ph]
work page internal anchor Pith review Pith/arXiv arXiv 1995
- [16]
-
[17]
M. Ramôa, P. G. Anastasiou, L. P. Santos, N. J. Mayhall, E. Barnes, and S. E. Economou, Reducing the resources required by adapt-vqe using coupled exchange operators and improved subroutines, npj Quantum Information11, 10.1038/s41534-025-01039-4 (2025)
-
[18]
A. Fitzpatrick, A. Nykänen, N. W. Talarico, A. Lunghi, S. Maniscalco, G. García-Pérez, and S. Knecht, Self- consistent field approach for the variational quantum eigensolver: Orbital optimization goes adaptive, The journal of physical chemistry. A128, 2843 (2024)
work page 2024
-
[19]
Z. W. Windom, L. Bertels, D. Claudino, and R. J. Bartlett, Toward the "platinum standard" of quan- tum chemistry on quantum computers: Perturbative quadruple corrections in unitary coupled cluster theory, The Journal of chemical physics164, 10.1063/5.0296472 (2026)
-
[21]
Y. S. Yordanov, V. Armaos, C. H. W. Barnes, and D. R. M. Arvidsson-Shukur, Qubit-excitation-based adaptive variational quantum eigensolver, Communica- tions Physics4, 10.1038/s42005-021-00730-0 (2021)
-
[22]
H. L. Tang, V. O. Shkolnikov, G. S. Barron, H. R. Grim- sley, N. J. Mayhall, E. Barnes, and S. E. Economou, Qubit-adapt-vqe: An adaptive algorithm for construct- ing hardware-efficient ansätze on a quantum proces- sor, PRX Quantum2, 10.1103/PRXQuantum.2.020310 (2021)
-
[23]
A. Kandala, A. Mezzacapo, K. Temme, M. Takita, M. Brink, J. M. Chow, and J. M. Gambetta, Hardware- efficient variational quantum eigensolver for small molecules and quantum magnets, Nature549, 242 (2017)
work page 2017
- [24]
-
[25]
G.-L. R. Anselmetti, D. Wierichs, C. Gogolin, and R. M. Parrish, Local, expressive, quantum-number-preserving vqe ansätze for fermionic systems, New Journal of Physics23, 113010 (2021)
work page 2021
-
[26]
J. Lee, W. J. Huggins, M. Head-Gordon, and K. B. Wha- ley, Generalized unitary coupled cluster wave functions for quantum computation, Journal of chemical theory and computation15, 311 (2019)
work page 2019
-
[27]
H. G. A. Burton, Accurate and gate-efficient quantum ansätze for electronic states without adaptive optimiza- tion, Physical Review Research6, 10.1103/PhysRevRe- search.6.023300 (2024)
-
[28]
J. R. McClean, S. Boixo, V. N. Smelyanskiy, R. Bab- bush, and H. Neven, Barren plateaus in quantum neural network training landscapes, Nature communications9, 4812 (2018)
work page 2018
-
[29]
H.R.Grimsley, G.S.Barron, E.Barnes, S.E.Economou, and N. J. Mayhall, Adaptive, problem-tailored varia- tional quantum eigensolver mitigates rough parameter landscapes and barren plateaus, npj Quantum Informa- tion9, 10.1038/s41534-023-00681-0 (2023)
-
[30]
P. Virtanen, R. Gommers, T. E. Oliphant, M. Haber- land, T. Reddy, D. Cournapeau, E. Burovski, P. Pe- terson, W. Weckesser, J. Bright, S. J. van der Walt, M. Brett, J. Wilson, K. J. Millman, N. Mayorov, A. R. J. Nelson, E. Jones, R. Kern, E. Larson, C. J. Carey, İ. Po- lat, Y. Feng, E. W. Moore, J. VanderPlas, D. Laxalde, J. Perktold, R. Cimrman, I. Henri...
work page 2020
- [31]
-
[32]
Y. Suzuki, Y. Kawase, Y. Masumura, Y. Hiraga, M. Nakadai, J. Chen, K. M. Nakanishi, K. Mitarai, R.Imai, S.Tamiya, T.Yamamoto, T.Yan, T.Kawakubo, Y. O. Nakagawa, Y. Ibe, Y. Zhang, H. Yamashita, H. Yoshimura, A. Hayashi, and K. Fujii, Qulacs: a fast and versatile quantum circuit simulator for research pur- pose, Quantum5, 559 (2021). 13
work page 2021
- [33]
-
[34]
D. A. Kreplin,Multiconfiguration Self-Consistent Field Methods for Large Molecules, Doctoral dissertation, Uni- versity of Stuttgart, Stuttgart, Germany (2020), accessed 2025-11-28
work page 2020
-
[35]
H. G. A. Burton, D. Marti-Dafcik, D. P. Tew, and D. J. Wales, Exact electronic states with shallow quantum cir- cuitsfromglobaloptimisation,npjQuantumInformation 9, 10.1038/s41534-023-00744-2 (2023)
-
[36]
S. Carter and N. C. Handy, The geometry of formalde- hyde, Journal of Molecular Spectroscopy179, 65 (1996)
work page 1996
-
[37]
E. R. Sayfutyarova, Q. Sun, G. K.-L. Chan, and G. Knizia, Automated construction of molecular active spaces from atomic valence orbitals, Journal of chemical theory and computation13, 4063 (2017)
work page 2017
-
[38]
K. Boguslawski, P. Tecmer, Ö. Legeza, and M. Reiher, Entanglement measures for single- and multireference correlation effects, The journal of physical chemistry let- ters3, 3129 (2012)
work page 2012
-
[39]
K. Boguslawski, K. H. Marti, O. Legeza, and M. Reiher, Accurate ab initio spin densities, Journal of chemical the- ory and computation8, 1970 (2012)
work page 1970
-
[40]
C. N. Sakellaris, E. Miliordos, and A. Mavridis, First principles study of the ground and excited states of feo, feo+, and feo(-), The Journal of chemical physics134, 234308 (2011)
work page 2011
-
[41]
S. R. Miller, N. E. Schultz, D. G. Truhlar, and D. G. Leopold, A study of the ground and excited states of al3 and al3(-). ii. computational analysis of the 488 nm anion photoelectron spectrum and a reconsideration of the al3 bonddissociationenergy,TheJournalofchemicalphysics 130, 024304 (2009)
work page 2009
-
[42]
S. Vancoillie, P. Å. Malmqvist, and V. Veryazov, Poten- tial energy surface of the chromium dimer re-re-revisited with multiconfigurational perturbation theory, Journal of chemical theory and computation12, 1647 (2016)
work page 2016
-
[43]
D. A. Kreplin, P. J. Knowles, and H.-J. Werner, Second- order mcscf optimization revisited. i. improved algo- rithms for fast and robust second-order casscf conver- gence, The Journal of chemical physics150, 194106 (2019)
work page 2019
-
[44]
C. J. Cramer, M. Włoch, P. Piecuch, C. Puzzarini, and L. Gagliardi, Theoretical models on the Cu2O2 torture track: Mechanistic implications for oxytyrosinase and small-moleculeanalogues,TheJournalofPhysicalChem- istry A110, 1991 (2006)
work page 1991
-
[45]
Q. Sun, X. Zhang, S. Banerjee, P. Bao, M. Barbry, N. S. Blunt, N. A. Bogdanov, G. H. Booth, J. Chen, Z.-H. Cui, J. J. Eriksen, Y. Gao, S. Guo, J. Hermann, M. R. Hermes, K. Koh, P. Koval, S. Lehtola, Z. Li, J. Liu, N. Mardirossian, J. D. McClain, M. Motta, B. Mussard, H. Q. Pham, A. Pulkin, W. Purwanto, P. J. Robin- son, E. Ronca, E. R. Sayfutyarova, M. Sc...
work page 2020
-
[46]
P. Jordan and E. P. Wigner, Über das paulische äquiv- alenzverbot, inThe Collected Works of Eugene Paul Wigner, The Scientific Papers, edited by A. S. Wight- man (Springer Berlin Heidelberg and Imprint Springer, Berlin, Heidelberg and s.l., 1993) pp. 109–129
work page 1993
-
[47]
J. Zhao, B.-B. Chen, Y.-C. Wang, Z. Yan, M. Cheng, and Z. Y. Meng, Measuring rényi entanglement entropy withhighefficiencyandprecisioninquantummontecarlo simulations, npj Quantum Materials7, 10.1038/s41535- 022-00476-0 (2022)
-
[48]
T. Grover, Entanglement of interacting fermions in quan- tum monte carlo calculations, Physical review letters 111, 130402 (2013)
work page 2013
- [49]
- [50]
-
[51]
J. P. Coe and M. J. Paterson, Investigating multiref- erence character and correlation in quantum chemistry, Journal of chemical theory and computation11, 4189 (2015)
work page 2015
-
[52]
Exponential Scaling Barriers for Variational Quantum Eigensolvers
M. Hagelueken, D. A. Kreplin, F. Wieland, M. F. Huber, and M. Roth, Research data supporting "Exponential Scaling Barriers for Variational Quantum Eigensolvers" (2026)
work page 2026
- [53]
-
[54]
V. O. Shkolnikov, N. J. Mayhall, S. E. Economou, and E. Barnes, Avoiding symmetry roadblocks and minimiz- ing the measurement overhead of adaptive variational quantum eigensolvers, Quantum7, 1040 (2023)
work page 2023
- [55]
-
[56]
A. Ikhtiarudin, G. K. Sunnardianto, F. Fathurrahman, M. K. Agusta, and H. K. Dipojono, Shot-efficient adapt- vqe via reused pauli measurements and variance-based shot allocation (2025), arXiv:2507.16879 [quant-ph]
-
[57]
A. Nykänen, M. A. C. Rossi, E.-M. Borrelli, S. Man- iscalco, G. García-Pérez, and A. Glos, Mitigating the measurement overhead of adapt-vqe with optimized in- formationally complete generalized measurements, Phys- ical Review Research7, 10.1103/t1hr-y7c8 (2025)
-
[58]
N. Vaquero-Sabater, A. Carreras, and D. Casanova, Pruned-adapt-vqe: Compacting molecular ansätze by re- moving irrelevant operators, Journal of chemical theory and computation21, 8720 (2025)
work page 2025
-
[59]
C. Feniou, M. Hassan, D. Traoré, E. Giner, Y. Ma- day, and J.-P. Piquemal, Overlap-adapt-vqe: practical quantum chemistry on quantum computers via overlap- guided compact ansätze, Communications Physics6, 10.1038/s42005-023-01312-y (2023)
-
[60]
J. W. Mullinax, P. G. Anastasiou, J. Larson, S. E. Economou, and N. M. Tubman, Classical preoptimiza- tion approach for adapt-vqe: Maximizing the potential of high-performance computing resources to improve quan- tumsimulationofchemicalapplications,Journalofchem- ical theory and computation21, 4006 (2025)
work page 2025
-
[61]
M. S. Rudolph, J. Miller, D. Motlagh, J. Chen, A. Acharya, and A. Perdomo-Ortiz, Synergistic pretrain- ingofparametrizedquantumcircuitsviatensornetworks, Nature communications14, 8367 (2023)
work page 2023
-
[62]
M. Hagelueken, M. F. Huber, and M. Roth, Data efficient prediction of excited-state properties using quantum neu- ral networks, New Journal of Physics27, 054508 (2025). 14
work page 2025
- [63]
-
[64]
J. Hermann, Pyberny: Geometry optimization in inter- nal coordinates,https://github.com/jhrmnn/pyberny, accessed: 2026-01-06
work page 2026
-
[65]
R. Guareschi and C. Angeli, The lowest singlet states of hexatriene revisited, Theoretical Chemistry Accounts 142, 10.1007/s00214-023-03064-y (2023)
-
[66]
S. Sakai, Theoretical study on the aromaticity of o- , m-, and p-benzyne, Journal of Molecular Structure: THEOCHEM715, 101 (2005). Appendix A: Comparison of VQE Variants The results presented in this study primarily em- ploy CEO-ADAPT-VQE with the TETRIS extension, with the exception of the six-orbital calibration in Section IIIB, which uses QEB-ADAPT-VQE...
work page 2005
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.