REVIEW 3 major objections 5 minor 28 references
Computing Reaction and Activation Energies of Pericyclic Reactions using a Symmetry-Adapted VQE Algorithm
T0 review · 3 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read Choosing each molecule's active space by a symmetry-matched fraction rule lets VQE reproduce CCSD reaction energies within 0.8 kcal/mol and activation energies within about 6 kcal/mol, even when absolute VQE energies are off by hundreds…
desk verdict Honest benchmark data, but the symmetry mechanism collapses for the C1 species where SMV just picks the largest active space. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The symmetry-matched fraction (SMF) is the percentage of UCCSD excitation operators inside a candidate active space that transform as the same irreducible representation as the Hartree-Fock reference; the symmetry-matched value (SMV) is the raw count, used when all excitations belong to one irreducible representation (C1 species). For each reactant, product, and transition state, the authors scan active spaces from (2e,3o) to (8e,8o) drawn from orbital pools (four occupied plus four virtual for C2v/Cs species, four plus two for C1 products, four plus five for transition states), then select the active space with the largest SMF/SMV. Energies come from VQE with a UCCSD ansatz, parity mapping with tapering, and SLSQP optimization on a noiseless simulator. The machinery ensures that the truncation of the wavefunction is symmetry-consistent across the reaction coordinate, which is what allows the large absolute errors to cancel in reaction and activation energy differences.
What would settle it
Compute reaction and activation energies for at least one additional pericyclic reaction (or perform a permutation test on the energy tables already supplied in the Supporting Information), comparing the maximum-SMF/SMV combination against randomly selected or chemically intuitive active-space combinations; if the symmetry-selected combination is not consistently closer to CCSD than typical arbitrary combinations, the claimed systematic cancellation is falsified.
Extended reading notes
Core claim
The central claim is that the cancellation of large absolute errors in energy differences is systematic when every species along the reaction coordinate is treated with the active space of highest SMF (or SMV for C1 species). Quantitatively: for the Diels-Alder reaction the SMF-VQE reaction energy is -34.31457 kcal/mol versus CCSD's -34.23219 (difference 0.08238 kcal/mol), and the activation energy is 19.36314 versus 24.34098 kcal/mol (difference 4.97784 kcal/mol); for the Alder-ene reaction the differences are 0.79230 and 6.24199 kcal/mol. These numbers rest on absolute VQE energies with average deviations of roughly 140-415 kcal/mol from CCSD. The same symmetry-guided criterion collapses 1008 (reaction) and 2592 (activation) possible active-space combinations for the Diels-Alder case and 558 and 1512 for the Alder-ene case into a single symmetry-consistent selection per energy. The paper further argues that active-space size alone is not a reliable guide: deviations are non-monotonic, and the transition states require the largest spaces (8e,8o) with 432 symmetry-matched excitations and 360 tapered variational parameters.
Load-bearing premise
The central assumption is that the hundreds-of-kcal/mol errors in absolute VQE energies cancel systematically in energy differences because active spaces are picked by the same SMF/SMV rule for every species; this cancellation is demonstrated for only two reactions, without error bars or a comparison against arbitrary active-space choices, so if it is a coincidence the reported accuracy would not generalize.
Editorial extensions
If this is right
- For the two reactions tested, reaction energies are reproduced within 1 kcal/mol of CCSD (0.08 and 0.79 kcal/mol) and activation energies within about 5-6 kcal/mol, despite absolute errors of 140-415 kcal/mol.
- The SMF/SMV criterion reduces the combinatorial search over active spaces to a single choice per energy, eliminating the need for chemical intuition to pick active spaces.
- Transition states are the resource bottleneck: both require an (8e,8o) active space with 432 symmetry-matched excitations and 360 tapered variational parameters, consistent with their near-degeneracy.
- Larger active spaces do not always improve VQE energetics within the truncated framework; the symmetry character of excitations, not active-space size, is the relevant selection principle.
- The protocol directly transfers from the earlier ring-strain studies to pericyclic reactions, indicating it may be a general strategy for reaction energetics in resource-constrained VQE.
Reading between the lines
- If the cancellation observed here is systematic, the SMF/SMV rule could serve as a parameter-free error-mitigation strategy for NISQ-era reaction energetics, complementing hardware error mitigation.
- A natural testable extension is a permutation analysis over the already-tabulated active-space energies: checking whether the maximum-SMF combination ranks better than typical random active-space combinations would quantify how often the criterion is genuinely predictive.
- Because C1 species have no symmetry distinction among excitations, the SMV criterion reduces to preferring larger excitation manifolds; testing more low-symmetry, substituted pericyclic reactions would clarify whether the rule is symmetry physics or simply a richness preference.
- The same logic may apply to other truncated correlated methods, such as selected CI or density matrix embedding, where absolute errors are large but relative energies along a reaction coordinate matter, suggesting SMF-style selection as a general truncation principle.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper applies a symmetry-based active-space selection criterion (SMF/SMV) to VQE/UCCSD calculations of the Diels–Alder and Alder–ene reactions in the STO-3G basis, and compares the resulting reaction and activation energies with CCSD reference values. The authors report that, although absolute VQE energies deviate from CCSD by roughly 140–415 kcal/mol per species, the SMF/SMV-selected active spaces yield reaction energies within 0.08 and 0.79 kcal/mol and activation energies within about 5–6 kcal/mol of CCSD. The paper further claims that this systematic cancellation of errors arises from the symmetry consistency of the selected active spaces across reactants, products, and transition states, and that the protocol reduces the large combinatorial space of active-space choices to a single selection per reaction.
Significance. If the reported cancellation is robust, the approach could be a useful practical heuristic for resource-constrained VQE simulations of correlated molecular systems, because it would make chemically meaningful relative energies accessible despite large absolute-energy errors in a minimal-basis, noiseless-simulator setting. The paper's strengths include the explicit enumeration of all candidate active spaces, the definition of the selection rule independently of the target energetics, and the quantification of the combinatorial reduction. However, the significance is substantially limited by the small number of reactions (two), the absence of uncertainty estimates, and the fact that for C1-symmetric species the criterion degenerates to a size-maximization rule, so the central symmetry-based mechanism is not actually operative for most of the species that determine the reported energy differences.
major comments (3)
- [Tables III and IV] For the five C1-symmetric species in this study (Diels–Alder product and TS; Alder–ene ene, product, and TS), all UCCSD excitation operators transform as the single irreducible representation of the C1 point group, so SMF = 100% for every candidate active space. Consequently SMV equals the total number of excitations, and maximizing SMV is mathematically identical to maximizing the number of excitations; within the finite orbital pools of Tables I and II this always selects the largest available active space, namely (8e,6o) in the (4o,2v) pool and (8e,8o) in the (4o,5v) pool, exactly as Table V reports. The paper acknowledges that SMF is identical for C1 species but does not state that the criterion reduces to a size-maximization heuristic. The central interpretation—that symmetry-consistent active-space selection produces systematic error cancellation—is therefore not operative for the product and transition states that determine the reported reaction and activation energies. Please provide control calculations (e.g., second-largest, smallest, or randomly selected active spaces, or a comparison against the distribution of all combinatorial energy differences) to determine whether the observed agreement is a symmetry effect or simply a size effect.
- [Section IV.C] The claim of systematic cancellation rests on only two reactions and two activation-energy comparisons, with no uncertainty estimates. The absolute per-species VQE errors are 140–415 kcal/mol, so the reported agreement requires that errors for species along the same reaction coordinate are highly correlated. Without a distribution of reaction and activation energies over the combinatorial space—which the authors already enumerate (e.g., 1008 combinations for the Diels–Alder reaction energy)—or repeated-run statistics to account for optimizer variability, the reported agreement could be accidental. Please quantify the spread of energy differences over all active-space combinations and locate the SMF/SMV-selected values within that distribution to demonstrate that the selection is not merely cherry-picking a favorable outlier.
- [Section IV.C] The statement that "the symmetry-guided framework adapts non-uniformly across the reaction coordinate and does not simply favor the largest active spaces for all molecular species" is misleading. It is true that the Cs and C2v species in this study can select compact active spaces, but for every C1 species the criterion selects the largest active space in the pool. The apparent non-uniform behavior is entirely driven by the two non-C1 species (diene and dienophile in Diels–Alder; enophile in Alder–ene), and the claim overstates the role of symmetry. Please rephrase this passage to explicitly acknowledge that for C1 species the criterion reduces to a maximum-excitation-count rule, and discuss how this affects the mechanistic interpretation of the error cancellation.
minor comments (5)
- [Supporting Information] The abstract contains the typo "Diel-Alder" (should be "Diels-Alder"); the same typo appears in the Supporting Information figure caption "Diene and Dieneophile" (should be "dienophile").
- [Section III.A] The main text states that absolute energies are provided in the Supporting Information, but the SI contains only comparison figures (S5 and S6) without numeric tables. Please include a table of absolute VQE and CCSD energies for all species and all active spaces considered, to allow readers to reproduce the reported differences and the cancellation analysis.
- [Section III.C] The paper specifies that CCSD/STO-3G single-point energies are computed on M05-2X/cc-pVDZ geometries, but it does not report the point-group symmetry enforcement or the irrep labels used in the SMF/SMV analysis. Please state explicitly which Abelian point group was used for each species and how the orbital irreps were assigned (e.g., from NWChem output) to make the selection reproducible.
- [Section III.C] Reference [14] is cited as "The Journal of Physical Chemistry A (2026)" without volume, article number, or DOI, and reference [28] is an arXiv preprint; please update these citations and briefly describe how the symmetry-matched fraction criterion was validated in those earlier works, since the present paper relies heavily on that prior framework.
- [Section III.C] The sentence "Active spaces with the largest SMF values were selected" is correct for the general case but could be misleading for C1 species, where all SMF values are identical; the subsequent text clarifies that SMV is used, but the initial wording might confuse readers. Please rephrase to state upfront that for C1 species the selection uses SMV (or, equivalently, the total excitation count).
Circularity Check
No significant circularity; the C1 SMV criterion degenerates to excitation-count maximization by definition, but the energetic comparison to CCSD is independent and unfitted.
-
self definitional
[Section III.C, Eq. (1), and the paragraph beginning 'As organic molecules generally possess lower symmetry'; Table V]
"SMF = Nsame-irrep/Ntotal-excitations × 100 (1) ... The numerator in the above expression will be referred to as the symmetry-matched value (SMV). ... In C1 symmetry, all excitations belong to the same irreducible representation, resulting in identical SMF values for all active spaces. In such cases, the SMV was used to guide active-space selection instead of the SMF."
By Eq. (1), SMV is Nsame-irrep. For a C1 species every UCCSD excitation transforms as the sole irreducible representation, so Nsame-irrep = Ntotal-excitations identically, making SMV = Ntotal-excitations by construction. 'Maximize SMV' is therefore exactly 'maximize the total number of excitations.' Within the finite pools used (4o2v for products and the ene; 4o5v for transition states), this selects AS(8e,6o) and AS(8e,8o), the largest active spaces, exactly as Table V reports for the product, TS, and ene species. Thus the 'symmetry-consistent' mechanism adds no independent constraint for five of the eight species; the selection reduces to a size-maximization heuristic.
full rationale
The central energetic claims are not circular. SMF/SMV values are computed from symmetry labels of UCCSD excitation operators using Eq. (1), before any comparison with CCSD, and the CCSD reference is external rather than used to tune the active-space choice. For the C2v and Cs species (diene, dienophile, enophile), the criterion selects non-maximal active spaces such as AS(2e,3o) for the enophile with SMF = 100, so the selection has independent content there. The self-citations to the authors' prior SMF work (refs 14 and 28) are not load-bearing because the criterion is defined in the present paper and tested against new CCSD data. The one definitional degeneracy is the C1 case documented above: for those species SMV collapses to the total excitation count, so the selected active spaces are the largest by construction. This affects the interpretation of the selection mechanism but not the derivation of the reported energy differences. The small sample size and lack of error bars are robustness concerns, not evidence of circularity. Overall, the paper's main numerical result is self-contained and not forced by its inputs.
Assumptions & free parameters
free parameters (2)
- Orbital pool composition per point group =
4 occupied + 4 virtual (C2v/Cs), 4 + 2 (C1 reactants/products), 4 + 5 (transition states)
- Active-space scan range =
(2e,3o) to (8e,8o); min and max orbitals in Tables I and II
assumptions (5)
- domain assumption UCCSD ansatz within the chosen active space is flexible enough to describe the correlated ground states and transition states.
- domain assumption CCSD/STO-3G single-point energies at DFT geometries are a valid reference for reaction and activation energies.
- domain assumption M05-2X/cc-pVDZ optimized geometries are accurate enough that single-point energies at these geometries represent the reaction path.
- ad hoc to paper Excitations transforming as the same irreducible representation as the HF reference are the ones that should define the active space.
- standard math Parity mapping with tapering preserves the relevant low-lying spectrum.
Cite this review
Pith. "Pith review of Computing Reaction and Activation Energies of Pericyclic Reactions using a Symmetry-Adapted VQE Algorithm." pith.science (2026). https://pith.science/paper/KX4BVPUQ
@misc{pith2026260804751,
author = {Pith},
title = {Pith review of: Computing Reaction and Activation Energies of Pericyclic Reactions using a Symmetry-Adapted VQE Algorithm},
year = {2026},
howpublished = {\url{https://pith.science/paper/KX4BVPUQ}},
note = {Machine review of arXiv:2608.04751}
}
read the original abstract
Pericyclic reactions provide stringent tests for quantum simulations because their mechanisms are governed by orbital symmetry and involve correlated transition states. In this work, we employ the variational quantum eigensolver (VQE) combined with a previously established symmetry-guided active-space selection protocol based on symmetry-matched fractions (SMF-VQE) to simulate Diel-Alder and Alder-ene reactions in complex systems involving extended pi-conjugation and multiple bonding. Although absolute electronic energies obtained from the current protocol exhibit significant deviations from the values computed using CCSD method, the symmetry-guided active spaces yield substantial cancellation of deviations in the energy differences. As a result, reaction energies are predicted with error (relative to CCSD) less than one kcal per mol, while activation energies are reproduced within about five to six kcal per mol. The symmetry-guided protocol also reduces the large combinatorial space of active-space choices to a single symmetry-consistent selection for each reaction.
Reference graph
Works this paper leans on
-
[1]
Y. Cao, J. Romero, J. P. Olson, M. Degroote, P. D. Johnson, M. Kieferová, I. D. Kivlichan, T. Menke, B. Peropadre, N. P. Sawaya,et al., Quantum chemistry in the age of quantum computing, Chem. Rev.119, 10856 (2019)
work page 2019
-
[2]
S. McArdle, S. Endo, A. Aspuru-Guzik, S. C. Benjamin, and X. Yuan, Quantum computational chemistry , Rev. Mod. Phys.92, 015003 (2020)
work page 2020
-
[3]
B. P. Lanyon, J. D. Whitfield, G. G. Gillett, M. E. Goggin, M. P. Almeida, I. Kassal, J. D. Biamonte, M. Mohseni, B. J. Powell, M. Barbieri,et al., Towards quantum chemistry on a quantum computer, Nat. Chem.2, 106 (2010)
work page 2010
-
[4]
Claudino, The basics of quantum computing for chemists, Int
D. Claudino, The basics of quantum computing for chemists, Int. J. Quantum Chem.122, e26990 (2022)
work page 2022
-
[5]
Kais, Introduction to quantum information and com- putation for chemistry , Quantum Inf
S. Kais, Introduction to quantum information and com- putation for chemistry , Quantum Inf. Comp. Chem. , 1 (2014)
work page 2014
-
[6]
G. Fano and S. Blinder, Quantum chemistry on a quantum computer, inMathematical Physics in Theoretical Chem- istry(Elsevier, 2019) pp. 377–400
work page 2019
- [7]
-
[8]
P. A. M. Dirac, Quantum mechanics of many-electron sys- tems, Proc. R. Soc. Lond. A.123, 714 (1929)
work page 1929
Show all 28 references
-
[9]
Preskill, Quantum computing in the nisq era and be- yond, Quantum2, 79 (2018)
J. Preskill, Quantum computing in the nisq era and be- yond, Quantum2, 79 (2018)
2018
-
[10]
Armaos, D
V. Armaos, D. A. Badounas, P. Deligiannis, and K. Lianos, Computational chemistry on quantum comput- ers: Ground state estimation, Appl. Phys. A.126, 625 9 (2020)
2020
-
[11]
R. P. Feynman, Simulating physics with computers, Int. J. Theor. Phys.21, 467 (1982)
1982
-
[12]
Aspuru-Guzik, A
A. Aspuru-Guzik, A. D. Dutoi, P. J. Love, and M. Head- Gordon, Simulated quantum computation of molecular energies, Science309, 1704 (2005)
2005
-
[13]
Peruzzo, J
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 pro- cessor, Nat. Commun.5, 4213 (2014)
2014
-
[14]
Sarkar, L
M. Sarkar, L. Roy , A. P. Gutal, A. Kumar, and M. Paran- jothy , Symmetry-match fraction for active-space selection: A simple route to accurate reaction energies in vqe/uccsd, The Journal of Physical Chemistry A (2026)
2026
-
[15]
Liepuoniute, M
I. Liepuoniute, M. Motta, T. Pellegrini, J. E. Rice, T. P. Gu- jarati, S. Gil, and G. O. Jones, Simulation of a diels–alder reaction on a quantum computer, Phys. Chem. Chem. Phys.26, 25181 (2024)
2024
-
[17]
D. A. McQuarrie and J. D. Simon,Physical chemistry: a molecular approach, Vol. 1 (University science books Sausalito, CA, 1997)
1997
-
[18]
Alteg, B
M. Alteg, B. Chevalier, O. Mestoudjian, and J.-L. Rossi, Study of adaptative derivative-assemble pseudo- trotter ansatzes in vqe through qiskit api (2022), arXiv:2210.15438 [quant-ph]
2022 arXiv
-
[19]
Singh, S
H. Singh, S. Majumder, and S. Mishra, Benchmarking of different optimizers in the variational quantum algo- rithms for applications in quantum chemistry , J. Chem. Phys.159(2023)
2023
-
[20]
Fleming,Pericyclic reactions(Oxford University Press, 2015)
I. Fleming,Pericyclic reactions(Oxford University Press, 2015)
2015
-
[21]
N. D. Epiotis, The theory of pericyclic reactions, Ange- wandte Chemie International Edition in English13, 751 (1974)
1974
-
[22]
Valiev, E
M. Valiev, E. J. Bylaska, N. Govind, K. Kowalski, T. P. Straatsma, H. J. J. Van Dam, D. Wang, J. Nieplocha, E. Aprà, T. L. Windus,et al., Nwchem: A comprehensive and scalable open-source solution for large scale molecu- lar simulations, Computer Physics Communications181, 1477 (2010)
2010
-
[23]
Electronic Hamiltonians were mapped to qubit operators using parity mapping with ta- pering to reduce the number of qubits
and Qiskit Nature 0.7.2 [24] with the STO-3G ba- sis and a UCCSD ansatz. Electronic Hamiltonians were mapped to qubit operators using parity mapping with ta- pering to reduce the number of qubits. The present systems involve extendedπinteractions, multiple bonding, and transit...
-
[24]
Javadi-Abhari, M
A. Javadi-Abhari, M. Treinish, K. Krsulich, C. J. Wood, J. Lishman, J. Gacon, S. Martiel, P. D. Nation, L. S. Bishop, A. W. Cross, B. R. Johnson, and J. M. Gambetta, Quantum computing with Qiskit (2024), arXiv:2405.08810 [quant- ph]
2024 arXiv
-
[25]
Q. N. developers and contributors, Qiskit nature 0.6.0 (2023)
2023
-
[26]
F. A. Cotton,Chemical applications of group theory(John Wiley & Sons, 1991)
1991
-
[27]
Mirman,Group theory: An intuitive approach(World Scientific, 1995)
R. Mirman,Group theory: An intuitive approach(World Scientific, 1995)
1995
-
[28]
P. H. Butler,Point group symmetry applications: methods and tables(Springer Science & Business Media, 2012)
2012
-
[29]
L. Roy , M. Sarkar, M. Tewari, A. Kumar, and M. Paran- jothy , Reaction-level consistency within the variational quantum eigensolver: Homodesmotic ring strain energies of cyclic hydrocarbons, arXiv preprint arXiv:2603.12737 (2026). Computing Reaction and Activation Energies of...
2026
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.