Recognition: unknown
Crossing the 12,000-atom barrier with heterogeneous quantum-classical supercomputing: quantum chemistry of protein-ligand complexes
Pith reviewed 2026-05-09 18:45 UTC · model grok-4.3
The pith
A hybrid quantum-classical method simulates protein-ligand complexes with over 12,000 atoms while matching coupled-cluster accuracy.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors show that quantum embedding allows a large molecule to be split into fragments whose electronic configurations can be sampled on 156-qubit processors and whose wavefunctions can be obtained via optimized subspace diagonalization on distributed classical hardware, producing accurate energies for dispersion- and electrostatics-dominated protein-ligand complexes at 11,608 and 12,635 atoms with over 40-fold size increase and up to 210-fold accuracy gain over prior state-of-the-art while reproducing CCSD fragment energies.
What carries the argument
Quantum embedding decomposition into fragments combined with heterogeneous sampling on quantum processors and subspace diagonalization on supercomputers, which lets independent fragment calculations reconstruct the full system's properties.
If this is right
- Protein-ligand binding energies become computable at near-CCSD quality for systems previously limited to much smaller sizes.
- The approach supplies a scalable route to systematically improvable ab initio simulations of biomolecular complexes.
- Fragment energies from the hybrid method can be used to benchmark or replace less accurate classical force fields in drug-design pipelines.
- The demonstrated parallel efficiency on supercomputers indicates that further increases in atom count remain feasible with additional hardware.
Where Pith is reading between the lines
- The same decomposition could be applied to molecular dynamics trajectories to obtain time-averaged energies without simulating every snapshot at full quantum cost.
- Extending the fragment size or adding more quantum processors might reduce the embedding error further, tightening agreement with experiment.
- The workflow creates a testbed for comparing quantum hardware performance against classical methods on chemically relevant, non-trivial Hamiltonians.
Load-bearing premise
The embedding step splits the molecule such that interactions between fragments are preserved with only negligible error, so separate simulations of the pieces can be added back together accurately.
What would settle it
Running a full coupled-cluster calculation on one of the smaller fragments or measuring an experimental observable such as binding free energy for the 11,608-atom complex and comparing it directly to the reported HQC result would show whether the claimed accuracy holds.
Figures
read the original abstract
Ab initio wavefunction methods provide accurate molecular simulations but their computational scaling restricts applications to small systems. We develop a workflow combining quantum embedding to decompose a molecule into fragments with a heterogeneous quantum-classical (HQC) method to simulate fragments. We sample fragment electronic configurations on two 156-qubit quantum processors (ibm$\_$cleveland, ibm$\_$kobe), using up to 94 qubits, running 9,200 circuits for over 100 hours, collecting $1.3 \cdot 10^9$ measurement outcomes - the most resource-intensive HQC computation for quantum chemistry to date. We compute fragment wavefunctions via optimized subspace diagonalization on two supercomputers (Fugaku, Miyabi-G), achieving 72.5$\%$ parallel efficiency with scalable distributed linear algebra kernels. We simulate two protein-ligand complexes spanning dispersion- and electrostatics-dominated regimes (11,608 and 12,635 atoms), demonstrate $>40\times$ increase in system size and up to $210\times$ improvement in accuracy over the previous state-of-the-art, with HQC matching coupled-cluster (CCSD) accuracy in fragment energies, and establish a scalable pathway for systematically improvable biomolecular simulations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents a heterogeneous quantum-classical (HQC) workflow that uses quantum embedding to decompose protein-ligand complexes into fragments, simulates the fragments on IBM quantum processors (up to 94 qubits, 9200 circuits, 1.3e9 measurements) and classical supercomputers (Fugaku, Miyabi-G) via optimized subspace diagonalization, and applies this to two systems of 11608 and 12635 atoms. It reports >40x larger systems and up to 210x accuracy gains over prior SOTA, with HQC fragment energies matching CCSD accuracy, while achieving 72.5% parallel efficiency.
Significance. If the embedding reconstruction of total energies is shown to be accurate, this constitutes a major advance in scaling wavefunction-based quantum chemistry to realistic biomolecular sizes, combining current quantum hardware with classical resources at record scale. The empirical benchmarking against independent CCSD calculations and the reported parallel efficiency on distributed linear algebra kernels are concrete strengths that support the demonstration of feasibility.
major comments (2)
- Abstract: The central claim is that the HQC method enables accurate simulation of the full 11k-12k atom protein-ligand complexes (spanning dispersion- and electrostatics-dominated regimes) with substantial accuracy gains. However, accuracy is only reported for individual fragment energies matching CCSD; no section quantifies the embedding error on the reconstructed total energy, nor demonstrates convergence with fragment size or bath depth. For long-range interactions, truncation of the fragment environment or approximate bath orbitals risks non-negligible many-body errors that do not cancel upon summation, which is load-bearing for the overall accuracy and size claims.
- Results and methods sections: The abstract states matching CCSD accuracy and large size gains but provides no details on validation protocols, error bars or statistical uncertainty on the 1.3e9 measurement outcomes, or explicit handling of hardware noise in the quantum runs. Without these, the robustness of the fragment-level CCSD matching and the claimed 210x accuracy improvement cannot be assessed.
Simulated Author's Rebuttal
We thank the referee for the constructive feedback and positive assessment of the work's significance. We address the two major comments point-by-point below, committing to revisions that strengthen the manuscript without altering its core claims.
read point-by-point responses
-
Referee: Abstract: The central claim is that the HQC method enables accurate simulation of the full 11k-12k atom protein-ligand complexes (spanning dispersion- and electrostatics-dominated regimes) with substantial accuracy gains. However, accuracy is only reported for individual fragment energies matching CCSD; no section quantifies the embedding error on the reconstructed total energy, nor demonstrates convergence with fragment size or bath depth. For long-range interactions, truncation of the fragment environment or approximate bath orbitals risks non-negligible many-body errors that do not cancel upon summation, which is load-bearing for the overall accuracy and size claims.
Authors: We agree that explicit validation of the reconstructed total energy is important for fully supporting the size and accuracy claims. The original manuscript reports fragment energies matching CCSD because the HQC simulation and accuracy benchmarks are performed at the fragment level, with the total energy obtained by summing embedded fragment contributions. To address the referee's concern regarding potential non-cancellation of long-range embedding errors, we will add a new subsection in Results that (i) describes the total-energy reconstruction formula, (ii) provides convergence data with fragment size and bath depth on representative subsystems, and (iii) reports estimated embedding errors on the total energy for the two target complexes. These additions will be based on additional post-processing of existing data and will not change the reported fragment-level results. revision: yes
-
Referee: Results and methods sections: The abstract states matching CCSD accuracy and large size gains but provides no details on validation protocols, error bars or statistical uncertainty on the 1.3e9 measurement outcomes, or explicit handling of hardware noise in the quantum runs. Without these, the robustness of the fragment-level CCSD matching and the claimed 210x accuracy improvement cannot be assessed.
Authors: The referee correctly identifies that the original text is insufficiently detailed on these points. The 1.3e9 shots were collected with standard mitigation (readout-error correction and dynamical decoupling) and statistical uncertainties were computed from binomial shot noise, but these protocols were only summarized. We will expand the Methods section with a dedicated subsection titled 'Quantum Measurement Protocol and Error Analysis' that (i) specifies the validation protocol (direct comparison of each fragment energy to classical CCSD on the same active space), (ii) reports per-fragment error bars derived from the measurement statistics, and (iii) quantifies residual hardware noise after mitigation. This will allow independent assessment of the 210x accuracy claim. revision: yes
Circularity Check
No circularity: empirical hardware demonstration benchmarked externally
full rationale
The paper reports a computational workflow applying quantum embedding decomposition followed by HQC simulation on real quantum processors and supercomputers for two large protein-ligand systems. Accuracy is established by direct comparison of fragment energies to independent CCSD calculations, with reported speedups and size increases as empirical outcomes. No derivation chain reduces a claimed result to fitted parameters, self-citations, or ansatzes by construction; the embedding error bound is treated as an assumption validated by the benchmarks rather than proven via internal equivalence. The work is self-contained against external references and does not rename known results or import uniqueness from prior author work in a load-bearing way.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Quantum embedding accurately decomposes the full electronic Hamiltonian into weakly coupled fragments
Reference graph
Works this paper leans on
-
[1]
Ultra-large-scale ab initio quantum chemical computation of bio-molecular systems: the case of spike protein of SARS-CoV-2 virus,
W.-Y . Ching, P. Adhikari, B. Jawad, and R. Podgornik, “Ultra-large-scale ab initio quantum chemical computation of bio-molecular systems: the case of spike protein of SARS-CoV-2 virus,”Comput. Struct. Biotechnol. J, vol. 19, pp. 1288–1301, 2021
2021
-
[2]
Hunting for quantum-classical crossover in condensed matter prob- lems,
N. Yoshioka, T. Okubo, Y . Suzuki, Y . Koizumi, and W. Mizukami, “Hunting for quantum-classical crossover in condensed matter prob- lems,”npj Quantum Inf, vol. 10, no. 1, p. 45, 2024
2024
-
[3]
Reliably assessing the electronic structure of cytochrome P450 on today’s classical computers and tomorrow’s quantum computers,
J. J. Goings, A. White, J. Lee, C. S. Tautermann, M. Degroote, C. Gid- ney, T. Shiozaki, R. Babbush, and N. C. Rubin, “Reliably assessing the electronic structure of cytochrome P450 on today’s classical computers and tomorrow’s quantum computers,”Proc. Natl. Acad. Sci. USA, vol. 119, no. 38, p. e2203533119, 2022
2022
-
[4]
Breaking the million- electron and 1 EFLOP/s barriers: biomolecular-scale ab initio molecular dynamics using MP2 potentials,
R. Stocks, J. L. G. Vallejo, F. C. Y . Yu, C. Snowdon, E. Palethorpe, J. Kurzak, D. Bykov, and G. M. J. Barca, “Breaking the million- electron and 1 EFLOP/s barriers: biomolecular-scale ab initio molecular dynamics using MP2 potentials,” inSC 24. IEEE Press, 2024
2024
-
[5]
Scaling correlated fragment molecular orbital calculations on Summit,
G. M. J. Barca, C. Snowdon, J. L. G. Vallejo, F. Kazemian, A. P. Rendell, and M. S. Gordon, “Scaling correlated fragment molecular orbital calculations on Summit,” inSC 22. IEEE Press, 2022
2022
-
[6]
Divide-and-conquer linear-scaling quantum chemical computations,
H. Nakai, M. Kobayashi, T. Yoshikawa, J. Seino, Y . Ikabata, and Y . Nishimura, “Divide-and-conquer linear-scaling quantum chemical computations,”J. Phys. Chem. A, vol. 127, no. 3, pp. 589–618, 2023
2023
-
[7]
D. G. Fedorov,Complete guide to the fragment molecular orbital method in GAMESS. World Scientific, 2023
2023
-
[8]
Fragment: an open-source framework for multiscale quantum chemistry based on fragmentation,
D. R. Broderick, P. E. Bowling, C. Brandt, S. Childress, J. Shockey, J. Higley, H. Dickerson, S. S. Ahmed, and J. M. Herbert, “Fragment: an open-source framework for multiscale quantum chemistry based on fragmentation,”WIREs Comput. Mol. Sci, vol. 15, p. e70058, 2025
2025
-
[9]
Fragment molecular orbital-based variational quantum eigensolver for quantum chemistry in the age of quantum computing,
H. Lim, D. H. Kang, J. Kim, A. Pellow-Jarman, S. McFarthing, R. Pellow-Jarman, H.-N. Jeon, B. Oh, J.-K. K. Rhee, and K. T. No, “Fragment molecular orbital-based variational quantum eigensolver for quantum chemistry in the age of quantum computing,”Sci. Rep, vol. 14, no. 1, p. 2422, 2024
2024
-
[10]
Multiscale quantum algorithms for quantum chemistry,
H. Ma, J. Liu, H. Shang, Y . Fan, Z. Li, and J. Yang, “Multiscale quantum algorithms for quantum chemistry,”Chem. Sci, vol. 14, p. 3190, 2023
2023
-
[11]
Quantum embedding theories to simulate condensed systems on quantum com- puters,
C. V orwerk, N. Sheng, M. Govoni, B. Huang, and G. Galli, “Quantum embedding theories to simulate condensed systems on quantum com- puters,”Nature Comp. Sci, vol. 2, pp. 424–432, 2022
2022
-
[12]
Hunting for quantum advantage in electronic structure calculations is a highly non-trivial task,
¨Ors Legeza, A. Menczer, M. A. Werner, S. S. Xantheas, F. Neese, M. Ganahl, C. Brower, S. R. Bernabeu, J. Hammond, and J. Gunnels, “Hunting for quantum advantage in electronic structure calculations is a highly non-trivial task,”arXiv:2603.28648, 2026
-
[13]
Heat-bath configuration interaction: An efficient selected configuration interaction algorithm inspired by heat-bath sampling,
A. A. Holmes, N. M. Tubman, and C. Umrigar, “Heat-bath configuration interaction: An efficient selected configuration interaction algorithm inspired by heat-bath sampling,”J. Chem. Theory Comput, vol. 12, no. 8, pp. 3674–3680, 2016
2016
-
[14]
H. Zhang and M. Otten, “From Random Determinants to the Ground State,”arXiv:2511.14734, 2025
-
[15]
Chemistry beyond the scale of exact diagonalization on a quantum- centric supercomputer,
J. Robledo-Moreno, M. Motta, H. Haas, A. Javadi-Abhari, P. Jurcevic, W. Kirby, S. Martiel, K. Sharma, S. Sharma, T. Shirakawaet al., “Chemistry beyond the scale of exact diagonalization on a quantum- centric supercomputer,”Sci. Adv, vol. 11, no. 25, p. eadu9991, 2025
2025
-
[16]
Quantum-centric compu- tation of molecular excited states with extended sample-based quantum diagonalization,
S. Barison, J. Robledo Moreno, and M. Motta, “Quantum-centric compu- tation of molecular excited states with extended sample-based quantum diagonalization,”Quant. Sci. Tech, vol. 10, no. 2, p. 025034, 2025
2025
-
[17]
T. Shirakawa, J. Robledo-Moreno, T. Itoko, V . Tripathi, K. Ueda, Y . Kawashima, L. Broers, W. Kirby, H. Pathak, H. Paiket al., “Closed- loop calculations of electronic structure on a quantum processor and a classical supercomputer at full scale,”arXiv:2511.00224, 2025
-
[18]
S. Piccinelli, S. Barison, A. Baiardi, F. Tacchino, J. Repp, I. Ron ˇcevi´c, F. Albrecht, H. L. Anderson, L. Gross, A. Curioniet al., “A note on large-scale quantum chemistry on quantum computers: the case of a molecule with half-M ¨obius topology,”arXiv:2603.08696, 2026
-
[19]
Krylov diagonalization of large many-body Hamiltonians on a quantum processor,
N. Yoshioka, M. Amico, W. Kirby, P. Jurcevic, A. Dutt, B. Fuller, S. Garion, H. Haas, I. Hamamura, A. Ivriiet al., “Krylov diagonalization of large many-body Hamiltonians on a quantum processor,”Nat. Comm, vol. 16, no. 1, p. 5014, 2025
2025
-
[20]
arXiv preprint arXiv:2512.17130 , year =
A. Shajan, D. Kaliakin, F. Liang, T. Pellegrini, H. Doga, S. Bhowmik, S. Das, A. Mezzacapo, M. Motta, and K. M. J. Merz, “Molecular quantum computations on a protein,”arXiv:2512.17130, 2025
-
[21]
Systematic improvability in quantum embedding for real materials,
M. Nusspickel and G. H. Booth, “Systematic improvability in quantum embedding for real materials,”Phys. Rev. X, vol. 12, p. 011046, 2022
2022
-
[22]
A. K. Patra, M. Mukherjee, A. Shukla, R. Bhatet al., “Towards chemically accurate and scalable quantum simulations on IQM quantum hardware: a quantum-HPC hybrid approach,”arXiv:2604.01983, 2026
-
[23]
Bridging physical intuition and hardware efficiency for correlated electronic states: the local unitary cluster Jastrow ansatz for electronic structure,
M. Motta, K. J. Sung, K. B. Whaley, M. Head-Gordon, and J. Shee, “Bridging physical intuition and hardware efficiency for correlated electronic states: the local unitary cluster Jastrow ansatz for electronic structure,”Chem. Sci, vol. 14, no. 40, pp. 11 213–11 227, 2023
2023
-
[24]
Improved parameter initialization for the (local) unitary cluster Jastrow ansatz,
W.-H. Lin, F. Liang, M. Motta, H. Zhang, K. M. M. Jr., and K. J. Sung, “Improved parameter initialization for the (local) unitary cluster Jastrow ansatz,”arXiv:2511.22476, 2025
-
[25]
Selected-Basis Diagonalization,
T. Shirakawa, “Selected-Basis Diagonalization,” Mar. 2026. [Online]. Available: https://github.com/r-ccs-cms/sbd
2026
-
[26]
GPU- accelerated selected-basis diagonalization with thrust for SQD-based algorithms,
J. Doi, T. Shirakawa, Y . Kawashima, S. Yunoki, and H. Horii, “GPU- accelerated selected-basis diagonalization with thrust for SQD-based algorithms,”arXiv:2601.16637, 2026
-
[27]
Recent developments in AMBER biomolecular simulations,
D. A. Case, D. S. Cerutti, V . W. D. Cruzeiro, T. A. Darden, R. E. Duke, M. Ghazimirsaeed, G. M. Giambasu, T. J. Giese, A. W. Gotz, J. A. Harris et al., “Recent developments in AMBER biomolecular simulations,”J. Chem. Inf. Model, vol. 65, no. 15, pp. 7835–7843, 2025
2025
-
[28]
H++: a server for estimating pKas and adding missing hydrogens to macromolecules,
J. Gordon, J. Myers, T. Folta, V . Shoja, L. Heath, and A. Onufriev, “H++: a server for estimating pKas and adding missing hydrogens to macromolecules,”Nucleic Acids Res, vol. 33, pp. 368–371, 2005
2005
-
[29]
Fast, efficient generation of high-quality atomic charges. AM1-BCC model: II. Parameterization and validation,
A. Jakalian, D. B. Jack, and C. I. Bayly, “Fast, efficient generation of high-quality atomic charges. AM1-BCC model: II. Parameterization and validation,”J. Comp. Chem, vol. 23, no. 16, pp. 1623–1641, 2002
2002
-
[30]
Complete reconstruction of an enzyme-inhibitor binding process by molecular dynamics simulations,
I. Buch, T. Giorgino, and G. De Fabritiis, “Complete reconstruction of an enzyme-inhibitor binding process by molecular dynamics simulations,” Proc. Natl. Acad. Sci. USA, vol. 108, no. 25, pp. 10 184–10 189, 2011
2011
-
[31]
The geometry of the reactive site and of the peptide groups in trypsin, trypsinogen and its complexes with inhibitors,
M. Marquart, J. Walter, J. Deisenhofer, W. Bode, and R. Huber, “The geometry of the reactive site and of the peptide groups in trypsin, trypsinogen and its complexes with inhibitors,”Struct. Sci, vol. 39, no. 4, pp. 480–490, 1983
1983
-
[32]
Homologous ligands accommodated by discrete conformations of a buried cavity,
M. Merski, M. Fischer, T. E. Balius, O. Eidam, and B. K. Shoichet, “Homologous ligands accommodated by discrete conformations of a buried cavity,”Proc. Natl. Acad. Sci. USA, vol. 112, p. 5039, 2015
2015
-
[33]
Energetic origins of specificity of ligand binding in an interior nonpolar cavity of T4- lysozyme,
A. Morton, W. A. Baase, and B. W. Matthews, “Energetic origins of specificity of ligand binding in an interior nonpolar cavity of T4- lysozyme,”Biochemistry, vol. 34, no. 27, pp. 8564–8575, 1995
1995
-
[34]
Co-design for A64FX manycore processor and “Fugaku
M. Sato, Y . Ishikawa, H. Tomita, Y . Kodama, T. Odajima, M. Tsuji, H. Yashiro, M. Aoki, N. Shida, I. Miyoshi, K. Hirai, A. Furuya, A. Asato, K. Morita, and T. Shimizu, “Co-design for A64FX manycore processor and “Fugaku”,” inSC 20, 2020, pp. 1–15
2020
-
[35]
Preliminary performance evaluation of Grace-Hopper GH200,
T. Hanawa, K. Nakajima, Y . Miki, T. Shimokawabe, K. Yamazaki, S. Sumimoto, O. Tatebe, T. Boku, D. Takahashi, A. Nukada, N. Fujita, R. Kobayashi, H. Tadano, and A. Naruse, “Preliminary performance evaluation of Grace-Hopper GH200,” inIEEE CLUSTER Workshops, 2024, pp. 184–185
2024
-
[36]
Software update: the ORCA program system - version 6.0,
F. Neese, “Software update: the ORCA program system - version 6.0,” WIREs Comput. Mol. Sci, vol. 15, no. 2, p. e70019, 2025
2025
-
[37]
Recent developments in the PySCF program package,
Q. Sun, X. Zhang, S. Banerjee, P. Bao, M. Barbry, N. S. Blunt, N. A. Bogdanov, G. H. Booth, J. Chen, Z.-H. Cuiet al., “Recent developments in the PySCF program package,”J. Chem. Phys, vol. 153, no. 2, 2020
2020
-
[38]
Cheap and near-exact CASSCF with large active spaces,
J. E. Smith, B. Mussard, A. A. Holmes, and S. Sharma, “Cheap and near-exact CASSCF with large active spaces,”J. Chem. Theory Comput, vol. 13, no. 11, pp. 5468–5478, 2017
2017
-
[39]
Block2: a comprehensive open-source framework to develop and apply state-of-the-art DMRG algorithms in electronic structure and beyond,
H. Zhai, H. R. Larsson, S. Lee, Z.-H. Cui, T. Zhu, C. Sun, L. Peng, R. Peng, K. Liao, J. T ¨olleet al., “Block2: a comprehensive open-source framework to develop and apply state-of-the-art DMRG algorithms in electronic structure and beyond,”J. Chem. Phys, vol. 159, no. 23, 2023
2023
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.