Pith. sign in

REVIEW 3 major objections 6 minor 85 references

Embedded quantum computing recovers experimental surface-reaction barriers and site preferences on copper by solving environment-aware active-space Hamiltonians.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.5

2026-07-30 14:09 UTC pith:TXOJRLC3

load-bearing objection Solid hybrid embedding + hardware QSCI on real Cu(111) benchmarks; frozen V_emb is a real but already-flagged accuracy risk, not a broken claim. the 3 major comments →

arxiv 2607.27009 v1 pith:TXOJRLC3 submitted 2026-07-29 quant-ph physics.chem-ph

Embedded quantum computing for many-body surface reaction

Dedong Wan , Xiaopeng Li , Yi Fan , Jie Liu , Xiongzhi Zeng , Zhenyu Li This is my paper
classification quant-ph physics.chem-ph PACS 03.67.Lx71.15.Mb82.65.+r
keywords quantum computingdensity-functional embeddingsurface chemistryQSCIheterogeneous catalysisCu(111)active spaceNEVPT2
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

Predictive catalysis needs many-body accuracy at the reactive site without discarding the extended metal surface. This paper introduces QC-DFET, which carves out a surface cluster, embeds it in a density-functional potential from the full slab, and reduces the problem to a compact active-space qubit Hamiltonian that stays chemically continuous along a reaction path. Hardware-sampled quantum-selected configuration interaction, with a classical second-order perturbation correction, then supplies the correlated energy. On Cu(111), the method matches bidirectional experimental H2 barriers to within a few hundredths of an electron-volt, restores the observed top-site preference and adsorption strength of CO, and brings the reverse barrier of the H2COO* branch into line with experiment while predicting a lower forward barrier to HCOOH*. The claim is that near-term quantum hardware can already act as a targeted many-body component inside realistic surface chemistry, not only as a solver for isolated molecules.

Core claim

QC-DFET shows that DFET-derived, reaction-consistent active-space Hamiltonians solved by hardware-sampled quantum-selected configuration interaction (up to 28 qubits) and corrected by strongly contracted NEVPT2 reproduce a hierarchy of experimentally constrained Cu(111) energetics: bidirectional H2 barriers within about 0.02 eV, CO top-site preference and adsorption strength inside the experimental window, and the H2COO* reverse barrier within 0.03 eV of experiment while identifying HCOOH* formation as the kinetically favored forward branch on the ideal terrace.

What carries the argument

QC-DFET: a density-functional embedding potential freezes the metallic environment into a cluster Hamiltonian; MBECAS-SR selects a compact, orbital-continuous active space along the reaction coordinate; hardware QSCI samples the important determinants; SC-NEVPT2 recovers residual dynamic correlation outside the active space.

Load-bearing premise

The embedding potential is optimized at the DFT level and then held fixed, so the metal environment never readjusts to the correlated electron density at the reaction site.

What would settle it

Recompute the H2 TS–IS and TS–FS barriers with an embedding potential updated self-consistently from the quantum-computing density; if the near-chemical-accuracy match to the experimental 0.77 eV and 0.59 eV barriers collapses, the fixed-potential premise fails.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • Near-term quantum processors can be used for determinant sampling inside embedded surface Hamiltonians rather than full variational optimization of the entire surface problem.
  • Keeping active orbitals continuous from reactant through transition state to product is as decisive as the quantum solver for balanced forward and reverse barriers.
  • Competing catalytic branches with different charge-redistribution patterns become separable at correlated accuracy on ideal metal terraces.
  • Larger embedded active spaces and non-ideal interfaces (steps, oxides, solvent, potential) are the natural next targets as hardware and selection protocols improve.

Where Pith is reading between the lines

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

  • On more ionic or oxide-supported catalysts, a frozen DFT embedding potential may systematically misplace charge flow at polarized transition states even if the active-space solver is exact.
  • Without automated active-space selection beyond chemically guided screening, the workflow will struggle to scale to full reaction networks with many intermediates.
  • Because the quantum role is subspace sampling plus classical diagonalization, classical selected-CI improvements remain competitive until embedded active spaces outgrow exact classical diagonalization.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. The manuscript introduces QC-DFET, a density-functional embedding framework that maps surface-reaction active spaces onto compact qubit Hamiltonians retaining a metallic environment via a DFET embedding potential. Reaction-consistent active spaces are built with MBECAS-SR; static correlation is treated by hardware-sampled QSCI on the Zuchongzhi processor (up to 28 qubits) and residual dynamic correlation by SC-NEVPT2. On Cu(111) the method is tested on three experimentally constrained problems: bidirectional H2 dissociation/desorption barriers, CO adsorption site preference and strength, and competing HCOO* hydrogenation branches. Reported QC-DFET energetics match experimental H2 barriers to ~0.01–0.02 eV, recover the top-site CO preference and adsorption energy within the experimental window, and bring the H2COO* reverse barrier to within 0.03 eV of experiment while predicting a lower forward barrier to HCOOH*.

Significance. If the reported accuracy holds under the stated approximations, the work supplies a concrete, experimentally anchored route for incorporating near-term quantum hardware into correlated surface chemistry rather than isolated molecular models. Strengths include a hierarchy of external experimental benchmarks (not only internal consistency), direct comparison to independent high-level methods (DMC, XYG3, prior ECASPT2/ENEVPT2), orbital-continuity and Mayer bond-order analyses that give chemically interpretable corrections, and an explicit non-claim of quantum advantage at ≤28 qubits. The quantum-hardware-in-the-loop design (sampling determinants, classical subspace diagonalization plus NEVPT2) is well matched to present devices and is a useful architectural contribution for the field.

major comments (3)
  1. [Methods, Density-functional embedding Theory; Discussion; Figs. 2b, 6c,d] Methods (Density-functional embedding Theory, Eqs. 1–4) and Discussion state that V_emb is DFT-optimized and held fixed during the correlated/QC calculation. Figs. 2b and 6c,d already show that emb-QC charge redistribution differs from emb-DFT at the H2 TS and along both formate branches—precisely where a frozen potential can bias state-balanced barriers or branching. The central claim of matched IS/TS/FS and site energetics under one environment-aware Hamiltonian therefore rests on an unquantified frozen-V_emb error. Please add either (i) a quantitative estimate of the barrier/site-energy sensitivity to V_emb (e.g., a one-shot update from the QC density, or a comparison against a self-consistent or alternative embedding potential) or (ii) a clear bound arguing why the residual error is smaller than the reported ~0.01–0.03 eV agreement with experiment. Without this, the state-balance cla
  2. [Results, Fig. 3; Methods, Dynamic-correlation correction, Eq. (12)] Fig. 3b shows that intermediate ADAPT-UCCGSD states with raw active-space errors of ~0.08–0.29 eV are brought within chemical accuracy only after SC-NEVPT2. Hardware QSCI (step-10 ansatz) is likewise not a fully converged active-space solution. The manuscript should state more explicitly, for each benchmark, the partition of the final QC-DFET energy correction into (a) the QSCI/selected-CI active-space contribution relative to emb-DFT and (b) the NEVPT2 increment. Without that decomposition it is hard to judge how much of the experimental agreement is carried by the quantum-sampled subspace versus the classical perturbative correction—information that is load-bearing for the claim that hardware-sampled QSCI is central to the workflow.
  3. [Supplementary Section 2; Results, Correlated reassessment of HCOO* hydrogenation] MBECAS-SR (Supplementary Section 2) still requires chemically chosen initial core orbitals and a system-dependent screening threshold r_ε (reported windows differ across H2, CO, and the two formate channels). The H2 comparison to prior ENEVPT2 active spaces (Fig. S4) is persuasive for that case, but for formate hydrogenation—where orbital character changes strongly along the path—please document sensitivity of the forward/reverse barriers to r_ε and to the initial core set (e.g., one tighter and one looser threshold, or omission of one core pair). A short sensitivity table would show that the branching conclusion (HCOOH* kinetically preferred; H2COO* reverse barrier near experiment) is not an artifact of a single active-space choice.
minor comments (6)
  1. [Abstract; Results] Abstract and opening claim “up to 28 qubits”; clarify in the main text which benchmark uses CAS(16e,14o)/28 qubits versus 24-qubit cases so readers can map hardware resources to each chemistry problem without hunting the SI.
  2. [Methods, Quantum-computing active-space solvers] Eq. numbering in the QSCI subsection restarts at (1)–(3) after earlier Eqs. (1)–(7); renumber continuously to avoid citation ambiguity.
  3. [Fig. 2a] Fig. 2a caption cites “XYG3:PBE-D3BJ (inferred from Fig. 3d of Ref. 14)”; state in the main text or SI how the inferred barriers were extracted so the comparison is reproducible.
  4. [Acknowledgements] Acknowledgements contain a duplicated phrase (“discussions and suggestions discussions and valuable guidance”); clean up.
  5. [Supplementary Figures] Several SI figure captions use H₂ / H2 inconsistently; standardize notation for H2, H2COO*, etc.
  6. [Methods, Dynamic-correlation correction and total-energy evaluation] In Methods, total-energy evaluation, “E_cluster^DFT” is introduced with slightly inconsistent subscript notation relative to E_cl^{emb-DFT}; align symbols with Eq. (12).

Circularity Check

0 steps flagged

No load-bearing circularity: QC-DFET energetics are computed from fixed DFT embedding and selected active spaces, then checked against external experiment and independent high-level methods.

full rationale

The central claims are bidirectional H2 barriers, CO site preference/adsorption strength, and formate-branch barriers on Cu(111). These are obtained from a standard embedding energy formula E_QC-DFET = E_tot^PW-DFT + (E_cl^emb-QC-NEVPT2 − E_cl^emb-DFT), with V_emb optimized once at DFT and held fixed, MBECAS-SR active spaces chosen by chemically guided correlation-energy increments and orbital continuity (not by fitting target barriers), and QSCI/SC-NEVPT2 solving that Hamiltonian. Reported numbers are then compared to external experimental anchors (H2 0.77/0.59 eV, CO ~−0.55 eV, H2COO* reverse barrier) and to independent methods (DMC, XYG3, ECASPT2, FEMION, prior ENEVPT2). Active-space thresholds and ADAPT step counts are set for compactness and hardware depth, not tuned to those benchmarks. Self-citations (MBECAS, Q2Chemistry, authors’ prior arXiv) supply tools and prior molecular workflow, not uniqueness theorems or fitted parameters that force the surface energetics. The frozen-V_emb approximation is an accuracy risk, not a by-construction identity between inputs and claimed predictions. Score 1 only for ordinary non-load-bearing methodological self-reference.

Axiom & Free-Parameter Ledger

3 free parameters · 5 axioms · 2 invented entities

The central claim rests on standard electronic-structure and quantum-computing machinery plus three methodological choices that are not independently forced: frozen DFT embedding, chemically seeded MBECAS-SR thresholds, and early-ADAPT + NEVPT2 error recovery. No new physical entities are postulated.

free parameters (3)
  • MBECAS-SR screening threshold r_ε = ~0.0001–0.0035 (system-dependent)
    Absolute correlation-energy-increment cutoff that decides which substrate and adsorbate orbitals enter the active space; system-dependent windows (e.g. 0.0001–0.0002 for H2, 0.0008–0.0013 for CO) are chosen by hand for compactness.
  • ADAPT step cutoff for hardware circuits (step-5 / step-10) = step-10 (production)
    Early ADAPT-UCCGSD operators are retained to keep circuit depth hardware-compatible; the cutoff is a pragmatic resource choice, not derived from a convergence theorem.
  • QSCI determinant retention (N_det or p_cut) and recovery iterations = ≤5 recovery iterations, 1e-6 Ha
    Subspace size and configuration-recovery loop (≤5 iterations, 1e-6 Ha) control the variational energy; values are set for robustness against device noise.
axioms (5)
  • domain assumption Wu–Yang density-matching DFET yields an embedding potential that correctly transfers the metallic environment into the cluster one-electron Hamiltonian.
    Invoked in Methods (Density-functional embedding Theory) and kept fixed thereafter; standard in the Carter-group DFET literature but not re-proven here.
  • domain assumption A compact active space selected by MBECAS-SR plus creeping captures the static correlation that differentiates IS/TS/FS and adsorption sites.
    Core premise of the reaction-consistent protocol (Results; SI §2); chemically guided initial cores are required.
  • domain assumption SC-NEVPT2 on an approximate quantum-selected wave function recovers residual dynamic correlation to chemical accuracy.
    Used to bring raw QC errors inside 0.043 eV (Fig. 3b); assumes the selected subspace already contains the dominant configurations.
  • domain assumption Jordan–Wigner mapped fermionic excitation circuits plus configuration recovery produce a sufficiently unbiased determinant sample on present superconducting hardware.
    Hardware QSCI workflow (Methods; SI §4); noise is mitigated but not eliminated.
  • domain assumption PBE-D3BJ geometries and ZPE corrections are adequate references for the correlated single-point energetics.
    All structures and vibrational corrections come from periodic PBE-D3BJ (Computational details).
invented entities (2)
  • QC-DFET workflow no independent evidence
    purpose: Name the combined DFET + MBECAS-SR + hardware QSCI + SC-NEVPT2 pipeline for surface reactions.
    Organizational label for existing components plus the new active-space protocol; not a new physical object.
  • MBECAS-SR no independent evidence
    purpose: Reaction-consistent active-space selection with orbital creeping for surfaces.
    Extension of the authors’ prior molecular MBECAS; selection rules and thresholds are defined in this work.

pith-pipeline@v1.2.0-daily-grok45 · 30369 in / 3372 out tokens · 58879 ms · 2026-07-30T14:09:28.208278+00:00 · methodology

0 comments
read the original abstract

Predictive simulations of catalytic interfaces require correlated electronic-structure treatments that describe localized chemical transformations while retaining the influence of the extended metallic environment. We introduce QC-DFET, a quantum-computing density-functional embedding framework that maps surface-reaction active spaces to compact, environment-aware qubit Hamiltonians. A reaction-consistent active-space protocol preserves orbital continuity along reaction coordinates, while quantum-selected configuration interaction based on measurements from the Zuchongzhi superconducting quantum processor and strongly contracted perturbation theory capture static and dynamic correlation. On Cu(111), QC-DFET treats active spaces up to 28 qubits and is validated through a hierarchy of experimentally constrained surface-chemistry challenges. H2 dissociation/desorption tests balanced bond breaking and recombination barriers, CO adsorption tests site selectivity and metal-adsorbate bonding, and formate hydrogenation tests competing hydrogenation branches with different kinetic and thermodynamic signatures. Across these cases, QC-DFET reproduces bidirectional H2 barriers, recovers the observed top-site preference and adsorption strength of CO, and reconciles the experimentally benchmarked H2COO* reverse barrier with the lower forward barrier to HCOOH*. These results establish embedded quantum computing as a practical route to correlated surface-reaction energetics.

discussion (0)

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

Reference graph

Works this paper leans on

85 extracted references · 5 linked inside Pith

  1. [1]

    K., Abild-Pedersen, F., Studt, F

    Nørskov, J. K., Abild-Pedersen, F., Studt, F. & Bligaard, T. Density functional theory in surface chemistry and catalysis. Proceedings of the National Academy of Sciences 108, 937–943 (2011)

    2011

  2. [2]

    Seh, Z. W. et al. Combining theory and experiment in electrocatalysis: insights into materials design. Science 355, eaad4998 (2017)

    2017

  3. [3]

    Birdja, Y . Y . et al. Advances and challenges in understanding the electrocatalytic conversion of carbon dioxide to fuels. Nature Energy 4, 732–745 (2019)

    2019

  4. [4]

    & Meyer, J

    Kroes, G.-J. & Meyer, J. Best-of-both-worlds computational approaches to difficult-to-model dissociation reactions on metal surfaces. Chemical Science 16, 480–506 (2025)

    2025

  5. [5]

    J., Mori-Sánchez, P

    Cohen, A. J., Mori-Sánchez, P. & Yang, W. Challenges for density functional theory. Chemical Reviews 112, 289–320 (2012)

    2012

  6. [6]

    Christensen, R., Hansen, H. A. & Vegge, T. Identifying systematic DFT errors in catalytic reactions. Catalysis Science & Technology 5, 4946–4949 (2015)

    2015

  7. [7]

    & Truhlar, D

    Verma, P. & Truhlar, D. G. Status and challenges of density functional theory. Trends in Chemistry 2, 302– 318 (2020)

    2020

  8. [8]

    Cao, K., Füchsel, G., Kleyn, A. W. & Juurlink, L. B. F. Hydrogen adsorption and desorption from Cu (111) and Cu (211). Physical Chemistry Chemical Physics 20, 22477–22488 (2018)

    2018

  9. [9]

    Díaz, C. et al. Chemically accurate simulation of a prototypical surface reaction: H₂ dissociation on Cu (111). Science 326, 832–834 (2009)

    2009

  10. [10]

    Feibelman, P. J. et al. The CO/Pt (111) puzzle. The Journal of Physical Chemistry B 105, 4018 –4025 (2001)

    2001

  11. [11]

    & Perdew, J

    Patra, A., Peng, H., Sun, J. & Perdew, J. P. Rethinking CO adsorption on transition-metal surfaces: effect of density-driven self-interaction errors. Physical Review B 100, 035442 (2019)

    2019

  12. [12]

    & Carter, E

    Zhao, Q. & Carter, E. A. Revisiting competing paths in electrochemical CO₂ reduction on copper via embedded correlated wavefunction theory. Journal of Chemical Theory and Computation 16, 6528–6538 (2020)

    2020

  13. [13]

    & Wei, C.-M

    Hsing, C.-R., Chang, C.-M., Cheng, C. & Wei, C.-M. Quantum Monte Carlo studies of CO adsorption on transition metal surfaces. The Journal of Physical Chemistry C 123, 15659–15664 (2019)

    2019

  14. [14]

    Chen, Z., Liu, Z. & Xu, X. Accurate descriptions of molecule-surface interactions in electrocatalytic CO₂ reduction on the copper surfaces. Nature Communications 14, 936 (2023)

    2023

  15. [15]

    Molecular orbital view of chemisorbed carbon monoxide

    Blyholder, G. Molecular orbital view of chemisorbed carbon monoxide. The Journal of Physical Chemistry 68, 2772–2777 (1964)

    1964

  16. [16]

    Grabow, L. C. & Mavrikakis, M. Mechanism of methanol synthesis on Cu through CO₂ and CO hydrogenation. ACS Catal. 1, 365–384 (2011)

    2011

  17. [17]

    Takeyasu, K. et al. Hydrogenation of formate species using atomic hydrogen on a Cu (111) model catalyst. Journal of the American Chemical Society 144, 12158–12166 (2022)

    2022

  18. [18]

    E., Gluth, O

    Ruehl, G., Harman, S. E., Gluth, O. M., LaV oy, D. H. & Campbell, C. T. Energetics of adsorbed formate and formic acid on Cu (111) by calorimetry. ACS Catalysis 12, 10950–10960 (2022)

    2022

  19. [19]

    Chen, Z., Liu, Z. & Xu, X. Clarifying the methanol synthesis mechanism via CO₂ hydrogenation on the Cu (111) surface: insights from accurate doubly hybrid density functionals. ACS Catalysis 15, 5039–5045 (2025)

    2025

  20. [20]

    Shi, B. X. et al. Many -body methods for surface chemistry come of age: achieving consensus with experiments. Journal of the American Chemical Society 145, 25372–25381 (2023)

    2023

  21. [21]

    Huang, Z. et al. Advancing surface chemistry with large-scale ab-initio quantum many-body simulations. arXiv preprint arXiv:2412.18553 (2024)

    Pith/arXiv arXiv 2024

  22. [23]

    Bartlett, R. J. & Musiał, M. Coupled-cluster theory in quantum chemistry. Reviews of Modern Physics 79, 291–352 (2007)

    2007

  23. [24]

    Chan, G. K. -L. & Sharma, S. The density matrix renormalization group in quantum chemistry. Annual Review of Physical Chemistry 62, 465–481 (2011)

    2011

  24. [25]

    Levine, D. S. et al. CASSCF with extremely large active spaces using the adaptive sampling configuration interaction method. Journal of Chemical Theory and Computation 16, 2340–2354 (2020)

    2020

  25. [26]

    O., Mosquera, M

    Jones, L. O., Mosquera, M. A., Schatz, G. C. & Ratner, M. A. Embedding methods for quantum chemistry: applications from materials to life sciences. Journal of the American Chemical Society 142, 3281 –3295 (2020)

    2020

  26. [27]

    Govind, N., Wang, Y . A. & Carter, E. A. Electronic-structure calculations by first-principles density-based embedding of explicitly correlated systems. The Journal of Chemical Physics 110, 7677–7688 (1999)

    1999

  27. [28]

    Klüner, T., Govind, N., Wang, Y . A. & Carter, E. A. Periodic density functional embedding theory for complete active space self-consistent field and configuration interaction calculations: ground and excited states. The Journal of Chemical Physics 116, 42–54 (2002)

    2002

  28. [29]

    & Carter, E

    Libisch, F., Huang, C. & Carter, E. A. Embedded correlated wavefunction schemes: theory and applications. Accounts of Chemical Research 47, 2768–2775 (2014)

    2014

  29. [30]

    & Carter, E

    Yu, K., Libisch, F. & Carter, E. A. Implementation of density functional embedding theory within the projector-augmented-wave method and applications to semiconductor defect states. The Journal of Chemical Physics 143 (2015)

    2015

  30. [32]

    Zhao, Q., Martirez, J. M. P. & Carter, E. A. Revisiting understanding of electrochemical CO₂ reduction on Cu (111): competing proton -coupled electron transfer reaction mechanisms revealed by embedded correlated wavefunction theory. Journal of the American Chemical Society 143, 6152–6164 (2021)

    2021

  31. [33]

    Cao, Y . et al. Quantum chemistry in the age of quantum computing. Chemical Reviews 119, 10856–10915 (2019)

    2019

  32. [34]

    & Yang, J

    Liu, J., Fan, Y ., Li, Z. & Yang, J. Quantum algorithms for electronic structures: basis sets and boundary conditions. Chemical Society Reviews 51, 3263–3279 (2022)

    2022

  33. [35]

    & Yang, J

    Zeng, X., Fan, Y ., Liu, J., Li, Z. & Yang, J. Quantum neural network inspired hardware adaptable ansatz for efficient quantum simulation of chemical systems. Journal of Chemical Theory and Computation 19, 8587–8597 (2023)

    2023

  34. [37]

    R., Economou, S

    Grimsley, H. R., Economou, S. E., Barnes, E. & Mayhall, N. J. An adaptive variational algorithm for exact molecular simulations on a quantum computer. Nature Communications 10, 3007 (2019)

    2019

  35. [38]

    & Wigner, E

    Jordan, P. & Wigner, E. Über das Paulische Äquivalenzverbot. Zeitschrift für Physik 47, 631–651 (1928)

    1928

  36. [39]

    Romero, J. et al. Strategies for quantum computing molecular energies using the unitary coupled cluster ansatz. Quantum Science and Technology 4, 014008 (2019)

    2019

  37. [40]

    J., Head -Gordon, M

    Lee, J., Huggins, W. J., Head -Gordon, M. & Whaley, K. B. Generalized unitary coupled cluster wave functions for quantum computation. Journal of Chemical Theory and Computation 15, 311–324 (2018)

    2018

  38. [41]

    S., Arvidsson-Shukur, D

    Yordanov, Y . S., Arvidsson-Shukur, D. R. M. & Barnes, C. H. W. Efficient quantum circuits for quantum computational chemistry. Physical Review A 102, 062612 (2020)

    2020

  39. [42]

    Zeng, X. et al. Transformer refined quantum sampling for strongly correlated electronic structure. arXiv preprint arXiv:2605.24617 (2026)

    Pith/arXiv arXiv 2026

  40. [43]

    Kanno, K. et al. Quantum-selected configuration interaction: classical diagonalization of Hamiltonians in subspaces selected by quantum computers. arXiv preprint arXiv:2302.11320 (2023)

    Pith/arXiv arXiv 2023

  41. [44]

    Robledo-Moreno, J. et al. Chemistry beyond the scale of exact diagonalization on a quantum -centric supercomputer. Science Advances 11, eadu9991 (2025)

    2025

  42. [45]

    Gujarati, T. P. et al. Quantum computation of reactions on surfaces using local embedding. npj Quantum Information 9, 88 (2023)

    2023

  43. [46]

    A., Neumann Barros Ferreira, R

    Barroca, M. A., Neumann Barros Ferreira, R. & Steiner, M. Exploring qubit -ADAPT-VQE for materials discovery in direct air capture. APL Quantum 1 (2024)

    2024

  44. [47]

    J., Chen, J

    Kattel, S., Ramírez, P. J., Chen, J. G., Rodriguez, J. A. & Liu, P. Active sites for CO₂ hydrogenation to methanol on Cu/ZnO catalysts. Science 355, 1296–1299 (2017)

    2017

  45. [48]

    Kordus, D. et al. Shape -dependent CO₂ hydrogenation to methanol over Cu₂O nanocubes supported on ZnO. Journal of the American Chemical Society 145, 3016–3030 (2023)

    2023

  46. [49]

    & Yang, W

    Wu, Q. & Yang, W. A direct optimization method for calculating density functionals and exchange – correlation potentials from electron densities. The Journal of Chemical Physics 118, 2498–2509 (2003)

    2003

  47. [50]

    & Malrieu, J

    Angeli, C., Cimiraglia, R. & Malrieu, J. -P. N -electron valence state perturbation theory: a fast implementation of the strongly contracted variant. Chemical Physics Letters 350, 297–305 (2001). Supplementary Information for Embedded quantum computing for many-body surface reaction Dedong Wan1, Xiaopeng Li1, Yi Fan1, Jie Liu2, Xiongzhi Zeng1*, Zhenyu Li1...

    2001

  48. [51]

    Grimme’s D3 dispersion correction 5,6 with the Becke –Johnson damping7 function was included throughout

    Computational details 1.1 Periodic DFT Calculations Periodic Kohn–Sham DFT calculations were performed using the Vienna Ab initio Simulation Package1,2 (V ASP, version 6.4.2) with the all-electron frozen-core projector augmented -wave (PAW)3 method and the PBE exchange –correlation functional 4. Grimme’s D3 dispersion correction 5,6 with the Becke –Johnso...

  49. [52]

    Therefore, active-space construction is essential and should be as automated and physics -based as possible

    MBECAS-SR active-space construction For certain chemical problems, quantitatively reliable results can be obtained by correlating only a limited number of orbitals and electrons directly involved in the reaction process. Therefore, active-space construction is essential and should be as automated and physics -based as possible. MBECAS was recently propose...

  50. [53]

    Definition and verification of the initial core orbitals. First, a mean-field calculation, such as HF or DFT, is performed to identify the core bonding orbitals of the adsorbed molecule that are directly involved in bond formation or cleavage. These orbitals can usually be recognized from chemically intuitive bo nding patterns. CASCI(2e,2o) calculations a...

  51. [54]

    MBECAS-SR(𝒜0) analysis is performed for the structure in which the adsorbed molecule is fully dissociated, referred to here as the IS

    MBECAS-SR screening at the reference structure. MBECAS-SR(𝒜0) analysis is performed for the structure in which the adsorbed molecule is fully dissociated, referred to here as the IS. We choose the IS as the reference because it exhibits the strongest wavefunction complexity along the reaction pathway and therefore requires the most complete orbital descri...

  52. [55]

    creeping

    Propagation along the potential -energy surface. To ensure robust and reaction -consistent convergence along the potential -energy surface (PES), a “creeping” strategy is employed. The converged CASSCF orbitals from the preceding geometry are used as the initial guess for the next geometry. The propagation begins from the fully dissociated IS, so that the...

  53. [56]

    For each electronic -structure level 𝑋 , where 𝑋 = emb-DFT or 𝑋 = emb-QC , we used the spin - summed AO density matrix 𝐷𝑋 and the AO overlap matrix 𝑆

    Wavefunction analysis 3.1 Orbital-resolved Mayer bond-order analysis for CO–Cu bonding To analyze how the correlated embedded quantum treatment modifies the CO –Cu bonding pattern, we computed orbital-resolved Mayer-type bond-order descriptors from the one-particle density matrices obtained at the emb-DFT and emb-QC levels. For each electronic -structure ...

  54. [57]

    S10a shows the convergence of the ADAPT-GUCCSD energy error and the corresponding circuit resources at representative ADAPT steps

    Hardware circuit construction and optimization 4.1 Implementation details Fig. S10a shows the convergence of the ADAPT-GUCCSD energy error and the corresponding circuit resources at representative ADAPT steps. The ADAPT -GUCCSD ansatz systematically lowers the energy error during the iterative operator selection and eventually reaches the chemical-accurac...

  55. [58]

    S1 Embedded Cu cluster models and representative optimized embedding potentials

    Supplementary Figures Fig. S1 Embedded Cu cluster models and representative optimized embedding potentials. First Cu 12 clusters were constructed from a √21 × √21 × 4 Cu (111) slab and used for the H 2 dissociation/desorption The second Cu12 cluster for cluster size robustness analysis and the Cu14 cluster used for the HCOO* + H* hydrogenation pathways to...

  56. [59]

    S1 CO adsorption energetics on the Cu 12 cluster

    Supplementary Tables Table. S1 CO adsorption energetics on the Cu 12 cluster. Energies are reported relative to the top adsorption site. Adsorption site Energy relative to top / eV top 0.00 bridge 0.51 fcc 0.46 hcp 0.73 abv 0.54 References:

  57. [60]

    & Furthmüller, J

    Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Computational Materials Science 6, 15–50 (1996)

    1996

  58. [61]

    & Furthmüller, J

    Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane- wave basis set. Physical Review B 54, 11169 (1996)

    1996

  59. [62]

    Blöchl, P. E. Projector augmented-wave method. Physical Review B 50, 17953 (1994)

    1994

  60. [63]

    P., Burke, K

    Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Physical Review Letters 77, 3865 (1996)

    1996

  61. [64]

    & Krieg, H

    Grimme, S., Antony, J., Ehrlich, S. & Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H –Pu. The Journal of Chemical Physics 132, 154104 (2010)

    2010

  62. [65]

    & Goerigk, L

    Grimme, S., Ehrlich, S. & Goerigk, L. Effect of the damping function in dispersion corrected density functional theory. Journal of Computational Chemistry 32, 1456–1465 (2011)

    2011

  63. [66]

    Becke, A. D. & Johnson, E. R. A density -functional model of the dispersion interaction. The Journal of Chemical Physics 123, 154101 (2005)

    2005

  64. [67]

    Monkhorst, H. J. & Pack, J. D. Special points for Brillouin-zone integrations. Physical Review B 13, 5188 (1976)

    1976

  65. [68]

    & Paxton, A

    Methfessel, M. & Paxton, A. T. High-precision sampling for Brillouin-zone integration in metals. Physical Review B 40, 3616 (1989)

    1989

  66. [69]

    Henkelman, G., Uberuaga, B. P. & Jónsson, H. A climbing image nudged elastic band method for finding saddle points and minimum energy paths. The Journal of Chemical Physics 113, 9901–9904 (2000)

    2000

  67. [70]

    & Jónsson, H

    Henkelman, G. & Jónsson, H. A dimer method for finding saddle points on high dimensional potential surfaces using only first derivatives. The Journal of Chemical Physics 111, 7010–7022 (1999)

    1999

  68. [71]

    Heyden, A., Bell, A. T. & Keil, F. J. Efficient methods for finding transition states in chemical reactions: comparison of improved dimer method and partitioned rational function optimization method. The Journal of Chemical Physics 123, 224101 (2005)

    2005

  69. [72]

    Zhao, Q., Zhang, X., Martirez, J. M. P. & Carter, E. A. Benchmarking an embedded adaptive sampling configuration interaction method for surface reactions: H₂ desorption from and CH₄ dissociation on Cu (111). Journal of Chemical Theory and Computation 16, 7078–7088 (2020)

    2020

  70. [73]

    & Carter, E

    Yu, K. & Carter, E. A. V ASPEmbedding. GitHub, https://github.com/EACcodes/V ASPEmbedding

  71. [74]

    Krauter, C. M. & Carter, E. A. EmbeddingIntegralGenerator. GitHub, https://github.com/EACcodes/EmbeddingIntegralGenerator

  72. [75]

    Sun, Q. et al. PySCF: the Python -based simulations of chemistry framework. Wiley Interdisciplinary Reviews: Computational Molecular Science 8, e1340 (2018)

    2018

  73. [76]

    Sun, Q. et al. Recent developments in the PySCF program package. The Journal of Chemical Physics 153, 024109 (2020)

    2020

  74. [77]

    Dunning, T. H. Jr. Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen. The Journal of Chemical Physics 90, 1007–1023 (1989)

    1989

  75. [78]

    Peterson, K. A. & Puzzarini, C. Systematically convergent basis sets for transition metals. II. Pseudopotential-based correlation consistent basis sets for the group 11 (Cu, Ag, Au) and 12 (Zn, Cd, Hg) elements. Theoretical Chemistry Accounts 114, 283–296 (2005)

    2005

  76. [79]

    & Preuss, H

    Dolg, M., Wedig, U., Stoll, H. & Preuss, H. Energy -adjusted ab initio pseudopotentials for the first row transition elements. The Journal of Chemical Physics 86, 866–872 (1987)

    1987

  77. [80]

    Fan, Y . et al. Q² Chemistry: a quantum computation platform for quantum chemistry. arXiv preprint arXiv:2208.10978 (2022)

    Pith/arXiv arXiv 2022

  78. [81]

    Saki, A. A. et al. Qiskit addon: sample -based quantum diagonalization. GitHub, https://github.com/Qiskit/qiskit-addon-sqd (2024)

    2024

  79. [82]

    Zeng, X. et al. Accurate chemical reaction modeling on a noisy intermediate-scale quantum computer with an active space-based workflow. The Journal of Physical Chemistry Letters 16, 11709–11716 (2025)

    2025

  80. [83]

    G., Pedersen, T

    Delcey, M. G., Pedersen, T. B., Aquilante, F. & Lindh, R. Analytical gradients of the state -average complete active space self-consistent field method with density fitting. The Journal of Chemical Physics 143, 044110 (2015)

    2015

Showing first 80 references.