Breaking the Speed–Accuracy Barrier in Quantum Chemistry

Compute reaction barriers with unprecedented accuracy at a fraction of today’s computational cost, enabling reliable simulations at scale.

Quantum chemistryDensity functional theoryComputational chemistry

Written with inputs from

  • Prof. Stefan GrimmeProfessor, University of Bonn

Predicting Reactions Before the Lab Requires Fast and Accurate Methods

Imagine being able to predict how molecules interact and react before ever stepping into a laboratory. Such a capability could transform how we discover new materials, design molecules, and develop advanced manufacturing processes. Rather than relying on costly and time-consuming trial and error, scientists could explore thousands of chemical reactions computationally, rapidly identifying promising molecules and processes for experimental validation.

Computational methods offer a path toward this vision, but a fundamental challenge remains: balancing accuracy and speed. Widely used methods like density functional theory (DFT) are computationally efficient but often struggle to predict chemical reactions accurately. More accurate approaches, such as coupled-cluster theory, come at a much greater computational cost and remain impractical at scale. This trade-off limits our ability to reliably predict and explore new chemistry at scale.

Quantum chemistry lets us predict how the electrons in molecules move, and that movement decides how a reaction happens: which bonds break, which form, and how fast. Here is a schematic of a trimethylaluminum molecule meeting water, the textbook reaction of atomic layer deposition. Their electron clouds merge, a hydrogen jumps across, and then methane drifts away.

The problem

Develop a more accurate exchange-correlation functional that improves DFT predictions without sacrificing computational efficiency.

Submissions

We evaluate submissions on a hidden benchmark of reference-quality molecular reaction energies, activation barriers, binding energies, and bond dissociation energies representative of atomic layer deposition (ALD) and chemical vapor deposition (CVD) precursor chemistry.

Submissions should include:

  • Method: A description of the proposed computational approach and its underlying scientific principles.
  • Implementation: Working code and instructions sufficient to independently run and evaluate the method.
  • Results: Evidence of predictive accuracy and computational efficiency on the public benchmark (GMTKN55).

Evaluation considers:

  • Predictive accuracy: Error against the hidden benchmark.
  • Speed and computational effort: Compute time per molecule.
  • Reproducibility: Sufficient documentation to verify the reported results.

LeaderboardProvisional

Error in kcal/mol (WTMAD-2, the weighted total mean absolute deviation of Goerigk et al., Phys. Chem. Chem. Phys. 19, 32184 (2017)). All functionals are evaluated using the def2-QZVP basis set on a single CPU core, and separately on a single A100 GPU. Each calculation is repeated three times, and the median runtime is used. Other common functionals, such as r²SCAN and Skala, are being evaluated and will be added to the leaderboard.

Error (kcal/mol)
WTMAD-2 · lower is better
1251020≤ 2× PBE cost2–5× PBE cost> 5× PBE cost2 s5 s10 s20 s50 s100 sCPU time per molecule · log scaletargetunder 1 kcal/molchemical accuracy1–3 kcal/molabove 3 kcal/molLDAPBEPBE0HSE06ωB97M-VB3LYPLDA: error 26.35 kcal/mol (WTMAD-2; 95% interval 22.24–30.66), 2.4 s per molecule (repeats 2.4–2.4); on one A100 GPU 1.33 s per molecule, 1.8 GPU-minutes for the whole suitePBE: error 12.47 kcal/mol (WTMAD-2; 95% interval 10.34–14.50), 3.0 s per molecule (repeats 3.0–3.0); on one A100 GPU 1.47 s per molecule, 2.0 GPU-minutes for the whole suitePBE0: error 8.37 kcal/mol (WTMAD-2; 95% interval 7.89–8.87), 3.9 s per molecule (repeats 3.9–3.9); on one A100 GPU 1.41 s per molecule, 1.9 GPU-minutes for the whole suiteHSE06: error 8.29 kcal/mol (WTMAD-2; 95% interval 7.72–8.85), 12.2 s per molecule (repeats 12.2–12.2); on one A100 GPU 1.94 s per molecule, 2.6 GPU-minutes for the whole suiteωB97M-V: error 3.81 kcal/mol (WTMAD-2; 95% interval 3.03–4.63), 43.5 s per molecule (repeats 43.5–43.8); on one A100 GPU 3.94 s per molecule, 5.3 GPU-minutes for the whole suiteB3LYP: error 16.16 kcal/mol (WTMAD-2; 95% interval 15.32–16.99), 4.0 s per molecule (repeats 4.0–4.0); on one A100 GPU 1.50 s per molecule, 2.0 GPU-minutes for the whole suite

Next: 04 · Tuning Materials with Light