# AI and quantum mechanics in drug discovery: how ML force fields speed it up

> AI speeds up drug discovery by learning quantum chemistry: ML force fields trained on DFT data simulate how drugs bind, near DFT accuracy and far faster.

- URL: https://computese.com/accelerating-drug-discovery-with-ai/
- Author: Duong Quan Nguyen, CEO, Computese
- Published: 2024-10-16
- Updated: 2026-09-25
- Topics: AI & automation

## In short
- How tightly a drug binds is set by electrons, so accurate energies need quantum mechanics. Exact methods cannot keep up: gold-standard coupled cluster is often impractical beyond a dozen atoms, and DFT cost grows roughly with the cube of the electron count.
- Machine-learned force fields such as ANI and MACE are trained on quantum chemistry data, then predict energies and forces close to that accuracy far faster: ANI-2x runs about a million times faster than the DFT it learned from.
- Open datasets drive the field: SPICE (2023) holds over 1.1 million drug-like conformations, and Open Molecules 2025 (OMol25) more than 100 million DFT calculations on systems of up to 350 atoms.
- In practice ML potentials join docking, free energy calculations and AlphaFold-style structures, often as hybrid ML/MM models. Quantum computers are still research: a 2022 estimate for one drug-metabolizing enzyme needed about 4.6 million physical qubits.
- Every prediction is limited by its training data and reference method. Charged, reactive or metal-containing molecules need extra checks, and only compounds made and measured in the lab count as results.

AI accelerates drug discovery with quantum mechanics by learning from it. Neural networks trained on quantum chemistry calculations, such as density functional theory (DFT), predict the energies and forces between atoms close to DFT accuracy but orders of magnitude faster, so simulations of how a drug binds its target can run at a scale exact calculations never could.

That matters because the numbers that decide whether a molecule becomes a drug candidate, such as how tightly it binds its target protein, are set by electrons, and computing electrons exactly is slow. This guide explains why drug design needs quantum mechanics, why DFT runs out of road, how neural network potentials such as ANI and MACE learn from open datasets like SPICE and Open Molecules 2025 (OMol25), where they sit next to docking, free energy calculations and AlphaFold-style structure prediction, what quantum computers may add later, and where all of it fails. How AI is used across the rest of the pipeline, from choosing targets to generative chemistry and clinical trials, is covered in [the role of AI in accelerating drug discovery](https://computese.com/the-role-of-ai-in-accelerating-drug-discovery/). This article is health information about research methods, not medical advice, and it recommends no product or treatment.

## Why drug binding is a quantum-mechanical problem

A small-molecule drug works by fitting into a pocket on a target protein and staying there. Whether it stays is a question of energy: the binding free energy, which sums everything that changes when the molecule leaves water and settles into the pocket. Those contributions are electronic. Hydrogen bonds, the way a molecule's electron cloud shifts near a charged side chain (polarization), the energy needed to twist a bond into the shape the pocket demands (a torsion), and any chemistry involving metal ions or bonds that form and break all follow from how electrons arrange themselves. Quantum mechanics is the physics of that arrangement. Solving the electronic Schrödinger equation for a fixed set of atomic positions gives the energy, and the force on each atom follows from how that energy changes as the atoms move.

Routine simulations in drug discovery usually leave the electrons out. A classical force field, also called molecular mechanics (MM), treats atoms as balls with fixed partial charges, joined by spring-like bonds and angles whose parameters are fitted in advance. The [MACE-OFF authors](https://arxiv.org/abs/2312.15211) note that such empirical force fields have dominated biomolecular simulation for over 50 years and are widely used in drug discovery, but generally lack the accuracy and transferability of first-principles methods. The weak spots are specific: a [2020 preprint on free energy calculations](https://www.biorxiv.org/content/10.1101/2020.07.29.227959v1) names the modelling of a ligand's internal energy, especially its torsions, as a major limitation of MM force fields.

The stakes are small numbers. Chemists usually treat an error of about 1 kcal/mol as the bar for chemical accuracy, and the [SPICE dataset paper](https://www.nature.com/articles/s41597-022-01882-6) uses the same standard. Experiments are not much tighter: when different groups measured the same protein-ligand affinities, their results differed by 0.77 to 0.95 kcal/mol (root-mean-square, depending on how the data was curated), according to a ChEMBL survey summarized in a [2023 benchmark paper](https://www.nature.com/articles/s42004-023-01019-9). A method that is off by 2 or 3 kcal/mol cannot rank candidates that measurements themselves barely separate.

Some chemistry cannot be approximated cheaply at all. Cytochrome P450 enzymes carry out the oxidation that is the most common way the human body metabolizes drugs, and their iron-containing active site has what chemists call multireference character: no single simple picture of the electrons is adequate. A [2022 PNAS study](https://pmc.ncbi.nlm.nih.gov/articles/PMC9499570/) chose it as a test case for quantum computers. We come back to it below.

## Why exact quantum chemistry is too slow for drug design

Quantum chemistry offers a ladder of methods that trade accuracy for speed. At the top, coupled-cluster theory with single, double and perturbative triple excitations, CCSD(T), extrapolated to the complete basis set limit (CBS), systematically approaches the exact solution and is considered a gold standard. It is also, according to the [ANI-1ccx authors](https://www.nature.com/articles/s41467-019-10827-4), often impractical for systems of more than a dozen atoms.

Density functional theory is the workhorse one rung down. It is much faster than coupled cluster, but it needs an empirically chosen density functional, so its results depend on that choice. Its cost still grows roughly with the cube of the number of electrons, which, the [OMol25 authors](https://arxiv.org/abs/2505.08762) write, stands in the way of routine large-scale screening and severely limits long simulations and structures of more than a few hundred atoms.

Now compare that with what a binding simulation needs. A small protein in water is already large: the MACE-OFF team's test of the protein crambin in explicit water had about 18,000 atoms. Molecular dynamics advances time in steps of one to a few femtoseconds (MACE-OFF's liquid simulations used 1 fs, the OpenMM 8 benchmarks 4 fs), so one nanosecond of motion takes hundreds of thousands to a million steps, each needing energies and forces for every atom. Free energy calculations need many such simulations per compound. Calling DFT at every step is out of reach, which is why the field relied on force fields in the first place.

| Method                     | What it models                               | Cost                                                          | Typical role in drug discovery                   |
| -------------------------- | -------------------------------------------- | ------------------------------------------------------------- | ------------------------------------------------ |
| Coupled cluster, CCSD(T)   | Electrons, close to exactly                  | Often impractical beyond about a dozen atoms at the CBS limit | Reference energies for small fragments           |
| Density functional theory  | Electrons, through an approximate functional | Grows roughly with the cube of the electron count             | Training data, fragment and reaction energies    |
| Machine-learned potential  | Energies learned from DFT or coupled cluster | Grows linearly with atoms; far below DFT, above force fields  | Ligand energetics and torsions, ML/MM simulation |
| Classical force field (MM) | Fixed charges, springs and fitted parameters | The cheapest physics-based option                             | Dynamics and free energy of whole proteins       |
| Docking scoring function   | A fast empirical score for one pose          | Cheap enough to screen 170 million compounds in a 2019 study  | Virtual screening of very large libraries        |

## How machine-learned force fields learn quantum chemistry

A machine-learned interatomic potential (MLIP), also called a neural network potential or machine learning force field, replaces the quantum calculation with a trained function. It takes the elements and 3D positions of the atoms and returns the energy. Forces come from the derivative of that energy with respect to the positions, so the model drops into molecular dynamics code wherever a force field would. What makes it scale is locality: the total energy is written as a sum of atomic contributions, each computed from the atoms within a cutoff radius of a few ångströms. Cost then grows linearly with the number of atoms, where the most efficient quantum methods grow with its cube, as the ANI-1ccx paper points out. For machine-learned potentials in crystals and other materials, see [AI in materials science](https://computese.com/revolutionizing-material-science-with-ai/).

Training is the expensive, one-time part. Researchers run DFT, or a more accurate method, on a large set of small molecules and molecular fragments in many distorted shapes, and the network learns to reproduce their energies and forces. Because each atom only sees its neighbourhood, a model trained on small fragments can then be applied to much larger systems, as long as the local environments look like the ones it was trained on.

![Small molecule fragments pass through a computer rack that calculates their energies. The results train an orange neural network, which then drives a simulation of a molecule inside a protein pocket.](https://computese.com/images/blog/accelerating-drug-discovery-with-ai/training.34cab8b40b-1536.webp)

*The quantum calculations are paid for once, on small fragments; the trained network reuses them at every step of a large simulation.*

Three lines of work stand out:

- **ANI.** [ANI-1](https://arxiv.org/abs/1610.08935), first posted in October 2016, built on a modified version of the symmetry functions of Behler and Parrinello. It was trained on molecules of up to 8 heavy atoms from the GDB databases, using hydrogen, carbon, nitrogen and oxygen, and matched reference DFT on molecules of up to 54 atoms with root-mean-square errors as low as 0.56 kcal/mol. ANI-1ccx (2019) was trained on DFT data first, then retrained by transfer learning on a set of CCSD(T)/CBS calculations chosen to span chemical space; it approaches coupled-cluster accuracy on benchmarks that include drug-like torsions and runs billions of times faster than CCSD(T)/CBS. [ANI-2x](https://europepmc.org/article/MED/32543858) (2020) added sulfur, fluorine and chlorine. Its authors note that those seven elements make up about 90% of drug-like molecules, and that it predicts DFT energies with a speed-up of about a million times.
- **MACE.** [MACE](https://arxiv.org/abs/2206.07697) (2022) is an equivariant message-passing network: atoms exchange information with their neighbours in a way that respects rotation, and its higher-order, four-body messages need only two message-passing layers, which keeps it fast. MACE-OFF, first posted in December 2023, applies it to organic molecules. It was trained on a neutral subset of the SPICE dataset covering ten elements, reproduces torsion scans of molecules it has not seen, and describes organic liquids and molecular crystals. Its authors also ran nanosecond simulations of the solvated protein crambin, at about 300,000 steps per day on one NVIDIA A100 GPU.
- **Models trained on OMol25.** The OMol25 release in 2025 came with baseline models trained on the new data, including eSEN, GemNet-OC, MACE and UMA, and published results for each on the same tests.

The speed-up is relative. The [OpenMM 8 paper](https://pmc.ncbi.nlm.nih.gov/articles/PMC10846090/) notes that ML potentials, although much faster than high-level quantum chemistry, are still orders of magnitude slower than most classical force fields, which often makes them impractical for whole macromolecules on their own. That gap is why the hybrid ML/MM simulations described below matter.

## The datasets behind the models: SPICE and OMol25

An ML potential is only as good as its reference data, and for years that data was small or narrow. Two open datasets now set the baseline for drug-like chemistry.

| Dataset | Published                     | Size                                                            | Reference method              | Coverage                                                                                              |
| ------- | ----------------------------- | --------------------------------------------------------------- | ----------------------------- | ----------------------------------------------------------------------------------------------------- |
| SPICE   | January 2023, Scientific Data | Over 1.1 million conformations                                  | DFT, ωB97M-D3(BJ)/def2-TZVPPD | Drug-like small molecules, dimers, dipeptides, solvated amino acids; 15 elements, charged and neutral |
| OMol25  | May 2025, arXiv preprint      | More than 100 million DFT calculations, systems up to 350 atoms | DFT, ωB97M-V/def2-TZVPD       | 83 elements: small molecules, biomolecules, metal complexes, electrolytes; varied charge and spin     |

SPICE was designed for one job: training potentials that simulate drug-like small molecules interacting with proteins. It stores forces as well as energies, plus multipole moments and bond orders, and its authors showed that models trained on it reach chemical accuracy across a broad region of chemical space.

Open Molecules 2025 is roughly a hundred times larger. Posted on 13 May 2025 by Meta's FAIR research lab with academic and national-laboratory collaborators, it holds more than 100 million DFT calculations (over 140 million in the paper's March 2026 revision) at one consistent, high level of theory, and took about 6.6 billion CPU core-hours by the authors' count. Its biomolecular part is aimed squarely at binding. Protein pockets with their bound ligands were cut out of experimental structures, and random drug-like molecules from the ZINC20, ChEMBL and GEOM databases were docked into pockets to add more. The authors also recomputed existing datasets, such as ANI-2x and SPICE2, at the same level of theory, so models are not trained on a patchwork of methods. The evaluations that ship with it include protein-ligand interaction energies and ligand strain: the energy difference between a ligand's bound shape and its most stable shape.

## Where ML potentials fit: docking, free energy and ML/MM

In a structure-based project, computational methods work as a funnel, and each stage is slower and more accurate than the one before.

**Docking** places each molecule in the pocket and scores the pose with a fast function. It is cheap enough for enormous libraries. In a [2019 Nature study](https://www.nature.com/articles/s41586-019-0917-9), a team docked 170 million make-on-demand compounds against two targets, the enzyme AmpC β-lactamase and the D4 dopamine receptor, then synthesized and tested 44 and 549 of the top-ranked molecules. It worked: it found new families of inhibitors, and crystal structures confirmed the predicted poses for AmpC. The same study also shows what a docking score is: for the D4 receptor, hit rates fell almost steadily as scores got worse, so a score shifts the odds rather than measuring an affinity.

**Free energy calculations** do the precise part. Free energy perturbation (FEP) and related alchemical methods simulate one ligand being gradually transformed into another, in the pocket and in water, to compute their difference in binding free energy. The 2023 benchmark paper, by scientists at Schrödinger testing FEP+ from the company's software suite, calls FEP the most consistently accurate way to predict relative binding affinities. Across 13,732 predictions, 57.5% were within 1 kcal/mol of experiment and 12.9% were off by more than 2 kcal/mol. The authors conclude that with careful preparation of protein and ligand structures, FEP can match the reproducibility of experiments. With standard MM force fields, these methods deliver 1 to 2 kcal/mol accuracy in well-behaved protein-ligand systems, according to the 2020 preprint.

**ML/MM** is where machine-learned potentials enter the workflow. The idea mirrors QM/MM: model the ligand, where force fields are weakest, with an ML potential, and the protein and water with a classical force field. OpenMM 8, the open-source molecular dynamics toolkit (published in December 2023), made this a few lines of code. Any [PyTorch](https://computese.com/latest-ai-tools-and-frameworks-a-comparative-analysis/) model can supply forces, and a mixed system is built from an ordinary force field plus a list of the atoms to treat with ML. At publication, it offered ANI-1ccx and ANI-2x as ready-made potentials.

![A small molecule sits in the cleft of a grey protein surrounded by water molecules. An orange dashed ring encloses only the small molecule, the part handled by the machine-learned potential.](https://computese.com/images/blog/accelerating-drug-discovery-with-ai/mlmm.e86891cae7-1536.webp)

*ML/MM spends accuracy where force fields are weakest, on the ligand, and keeps the rest of the system cheap.*

The payoff shows up in the numbers. In the 2020 preprint, researchers corrected standard MM free energy calculations with an ANI-2x/MM potential in a post-processing step. On a benchmark series of inhibitors of the kinase TYK2, the error in absolute binding free energies fell from 0.97 to 0.47 kcal/mol, for a small additional computational cost. It is one target and one ligand series, but it shows where the gain comes from: better energetics for the ligand itself.

## Where AlphaFold-style structure prediction fits

Every method above starts from a 3D structure of the target and a plausible pose for the ligand. When no experimental structure exists, deep learning now supplies one. [AlphaFold 3](https://www.nature.com/articles/s41586-024-07487-w), published in Nature on 8 May 2024, predicts the joint structure of complexes of proteins, nucleic acids, small molecules, ions and modified residues. On the PoseBusters benchmark, its authors report far greater accuracy for protein-ligand interactions than state-of-the-art docking tools. AlphaFold's wider place in AI-driven science is covered in [AI-enabled scientific discoveries](https://computese.com/ai-enabled-scientific-discoveries/).

Two limits keep physics in the loop. AlphaFold 3 predicts static structures like those in the Protein Data Bank, not how a complex moves in solution, and like other generative models it can invent plausible-looking structure in disordered regions. Newer co-folding models reach toward energies too. The open-source [Boltz-2](https://github.com/jwohlwend/boltz) (2025, MIT licence) predicts binding affinity alongside structure, and its developers describe it as the first deep learning model to approach the accuracy of FEP while running 1,000 times faster. That is the developers' own claim, on their own benchmarks.

In practice the tools stack rather than compete. An experimental or predicted structure provides the pocket and a starting pose; docking or co-folding narrows the list; FEP or ML/MM simulation tests the energetics and motion of the short list before anything is made.

## Will quantum computers speed up drug discovery?

Not yet. First, a distinction: everything above runs on ordinary computers. Quantum mechanics in drug discovery today means quantum chemistry software on CPUs and GPUs. Quantum computers are a separate technology that may one day run some of those calculations.

Their promise is the part classical methods handle worst: molecules whose electrons are strongly correlated, such as metal-containing active sites. The 2022 PNAS study, by Google Quantum AI with Boehringer Ingelheim and others, took cytochrome P450 as its test case. It estimated that computing the ground-state energy of the largest model of a key P450 intermediate would need about 4.6 million physical qubits and 73 hours at a physical error rate of 0.1%, or about 500,000 qubits and 25 hours if error rates fell to 0.001%. Its conclusion was careful: simulating these enzymes at a large enough scale has the potential to be a problem where quantum computers beat classical ones.

Set that against the hardware as announced:

- **December 2024:** Google's [Willow chip](https://blog.google/technology/research/google-willow-quantum-chip/), with 105 qubits, showed errors falling as more qubits were used for error correction.
- **June 2025:** IBM [set out a plan](https://www.ibm.com/quantum/blog/large-scale-ftqc) to deliver Starling, a fault-tolerant machine running 100 million gates on 200 logical qubits, by 2029.
- **October 2025:** Google reported [Quantum Echoes](https://blog.google/technology/research/quantum-echoes-willow-verifiable-quantum-advantage/) on Willow, which it calls the first verifiable quantum advantage, running 13,000 times faster than the best classical algorithm on a leading supercomputer. In a separate proof of principle with UC Berkeley, it studied two molecules of 15 and 28 atoms, and the results matched traditional NMR. Google suggests quantum-enhanced NMR could one day help show how drug candidates bind their targets.

> [!NOTE]
> Treat quantum computing for drug discovery as research. The P450 estimate needs millions of physical qubits; Willow has 105. None of these sources reports a drug candidate found or optimized on a quantum computer.

## The limits: training data, reference methods and the lab

Every stage in this chain can be wrong in a way the next stage does not see, and the failure modes are well documented.

A local model is blind beyond its cutoff. Many ML potentials let atoms interact directly only within a cutoff radius, typically around 6 Å, while the interaction between two charges fades slowly with distance, as 1/r. OMol25 tests this explicitly, with separate short-range and long-range evaluations, and the MACE-OFF authors state that without explicit long-range interactions their model applies only to neutral, non-radical and non-reactive systems.

![Atoms close to a central atom sit inside a dashed circle and are linked to it. A charged group farther away, drawn in orange, lies outside the circle, and its faint line of influence stops at the edge.](https://computese.com/images/blog/accelerating-drug-discovery-with-ai/cutoff.0cc1138953-1536.webp)

*A local model sees only what falls inside its cutoff, so charged groups and long-range effects need their own checks.*

- **A model is capped by its reference.** An ML potential trained on DFT inherits DFT's errors. In the MACE-OFF paper, DFT torsion barriers differed from coupled-cluster values by 0.2 kcal/mol on average, which the authors describe as the best their model could possibly reach on that test.
- **Accuracy drops outside the training data.** ANI-2x covers seven elements and ANI-1ccx four. On OMol25's out-of-distribution test, models did better on neutral organic molecules and biomolecules than on electrolytes and metal complexes.
- **Preparation matters.** In the 2023 FEP benchmark, adding both protonation states of the ligands for one target (TNKS2) cut the error from 2.10 to 1.60 kcal/mol, with no change to the method.
- **Structures are predictions.** AlphaFold 3's poses are static, and it can hallucinate structure in disordered regions.
- **Experiments disagree too.** With independent measurements differing by 0.77 to 0.95 kcal/mol, no method can be judged more precisely than the data it is compared with.

> [!IMPORTANT]
> Only the lab decides. Every result in the 2019 docking study is a compound that was synthesized and tested, and for AmpC, crystal structures confirmed the predicted poses. Everything upstream narrows what gets made; it does not replace making it.

Before trusting an ML potential on a project, work through these checks:

1. **Check the domain.** Compare the elements, charge states, metals and any reactive chemistry in your molecules with what the model was trained on.
2. **Spot-check against the reference.** Recompute a sample of your own structures at the training level of theory and compare energies and forces.
3. **Benchmark retrospectively.** Run the full workflow on a compound series you have already measured before using it to choose what to make.
4. **Record provenance.** Store the model version, training data and reference method with every prediction, so a result can be traced and repeated.

## The data engineering around the models

Much of the work in this field is data work. OMol25 alone took about 6.6 billion CPU core-hours, and its value depends on bookkeeping: one level of theory throughout, quality checks on every calculation, and test splits held out by composition so a model's error on unfamiliar chemistry is measured rather than assumed. Any team building on data like this faces the problems of every large machine learning dataset: versioned data, recorded provenance and evaluations that run before a new model is trusted.

Those are general data and machine learning engineering problems, not chemistry. We do not work in drug discovery or molecular modelling. Our [data platform](https://computese.com/services/data-platform/) service builds pipelines with tests, lineage and access control, and serves governed datasets for ML and AI. For teams putting language models into their own workflows, our [AI and automation](https://computese.com/services/ai-automation/) service tests every model, prompt or source change against an evaluation set of real cases before it goes live. More on the subject is under the [AI topic](https://computese.com/category/ai/).

## Key terms
- **Density functional theory (DFT)**: A quantum-mechanical method that computes a molecule's energy from its electron density through an approximate functional. The usual source of training data for ML potentials; its cost grows roughly with the cube of the number of electrons.
- **CCSD(T)**: Coupled-cluster theory with single, double and perturbative triple excitations, a gold-standard quantum chemistry method. Very accurate and very expensive: at the complete basis set limit, often impractical beyond about a dozen atoms.
- **Force field (molecular mechanics)**: A classical model that treats atoms as charged balls joined by springs, with parameters fitted in advance. Fast enough for whole proteins in water, but it has no electrons.
- **Machine-learned interatomic potential (MLIP)**: A neural network trained on quantum chemistry data that predicts the energy and forces of a set of atoms from their elements and positions. Also called a neural network potential or an ML force field.
- **Chemical accuracy**: An error of about 1 kcal/mol or less, the level chemists commonly use as the standard for computed energies.
- **Free energy perturbation (FEP)**: A simulation method that computes the difference in binding free energy between two related ligands by gradually transforming one into the other, in the protein and in water.
- **ML/MM**: A hybrid simulation in which part of the system, usually the ligand, is modelled with an ML potential and the rest, protein and water, with a classical force field.
- **Docking**: Placing a molecule in a protein's binding pocket and scoring the pose with a fast function, used to screen libraries of millions of compounds.
- **Out-of-distribution**: Inputs unlike a model's training data, such as new elements, charge states or reactions, where its errors are usually larger and harder to predict.

## Common questions

### How is quantum mechanics used in drug discovery?

Mostly as quantum chemistry on ordinary computers. Methods such as DFT compute the energies of drug-like molecules, pieces of binding pockets and reactions, where the behaviour of electrons decides the answer. The results are used directly for small systems and, increasingly, as training data for machine-learned force fields that carry near-quantum accuracy into simulations of whole proteins.

### What is a machine learning force field?

A neural network that predicts the energy of a group of atoms from their elements and positions, trained to reproduce quantum chemistry calculations. Its forces drive molecular dynamics the way a classical force field does. Examples include ANI-2x, which predicts DFT energies about a million times faster than DFT, and MACE-OFF for organic molecules.

### Is DFT accurate enough for drug design?

DFT is the standard reference for training data and fragment energies, but it is an approximation whose results depend on the chosen functional. It is also too slow for whole proteins over time, because its cost grows roughly with the cube of the electron count. Coupled cluster, CCSD(T), at the complete basis set limit is the more reliable reference, but is often impractical beyond about a dozen atoms.

### Can AI replace free energy perturbation (FEP)?

Not as a rule, yet. Hybrid ML/MM potentials already improve FEP: in a 2020 preprint, an ANI-2x/MM correction roughly halved the error on a TYK2 inhibitor benchmark, from 0.97 to 0.47 kcal/mol. Co-folding models such as Boltz-2 claim FEP-like affinity accuracy at far lower cost, but that is the developers' own benchmark, and any model needs testing on measured data for the target at hand.

### Will quantum computers speed up drug discovery?

Possibly, for molecules with strongly correlated electrons such as metal-containing enzymes, but it is still research. A 2022 estimate for a cytochrome P450 model needed about 4.6 million physical qubits at a 0.1% error rate. Google's Willow chip, announced in December 2024, has 105 qubits.

### What is the OMol25 dataset?

Open Molecules 2025, posted by Meta's FAIR lab and collaborators in May 2025: more than 100 million DFT calculations at the ωB97M-V/def2-TZVPD level, covering 83 elements and systems of up to 350 atoms, including protein-ligand pockets. It was released publicly with baseline models and evaluations for training and testing ML potentials.

## Sources
1. [Approaching coupled cluster accuracy with a general-purpose neural network potential through transfer learning](https://www.nature.com/articles/s41467-019-10827-4), Nature Communications
2. [The Open Molecules 2025 (OMol25) Dataset, Evaluations, and Models](https://arxiv.org/abs/2505.08762), arXiv (Meta FAIR and collaborators)
3. [MACE-OFF: Transferable Short Range Machine Learning Force Fields for Organic Molecules](https://arxiv.org/abs/2312.15211), arXiv
4. [SPICE, A Dataset of Drug-like Molecules and Peptides for Training Machine Learning Potentials](https://www.nature.com/articles/s41597-022-01882-6), Scientific Data
5. [The maximal and current accuracy of rigorous protein-ligand binding free energy calculations](https://www.nature.com/articles/s42004-023-01019-9), Communications Chemistry
6. [Towards chemical accuracy for alchemical free energy calculations with hybrid physics-based machine learning / molecular mechanics potentials](https://www.biorxiv.org/content/10.1101/2020.07.29.227959v1), bioRxiv
7. [Reliably assessing the electronic structure of cytochrome P450 on today's classical computers and tomorrow's quantum computers](https://pmc.ncbi.nlm.nih.gov/articles/PMC9499570/), PNAS (via PubMed Central)
8. [ANI-1: An extensible neural network potential with DFT accuracy at force field computational cost](https://arxiv.org/abs/1610.08935), arXiv
9. [Extending the Applicability of the ANI Deep Learning Molecular Potential to Sulfur and Halogens](https://europepmc.org/article/MED/32543858), Journal of Chemical Theory and Computation (via Europe PMC)
10. [MACE: Higher Order Equivariant Message Passing Neural Networks for Fast and Accurate Force Fields](https://arxiv.org/abs/2206.07697), arXiv
11. [OpenMM 8: Molecular Dynamics Simulation with Machine Learning Potentials](https://pmc.ncbi.nlm.nih.gov/articles/PMC10846090/), The Journal of Physical Chemistry B (via PubMed Central)
12. [Ultra-large library docking for discovering new chemotypes](https://www.nature.com/articles/s41586-019-0917-9), Nature
13. [Accurate structure prediction of biomolecular interactions with AlphaFold 3](https://www.nature.com/articles/s41586-024-07487-w), Nature
14. [Boltz: official repository for the Boltz biomolecular interaction models](https://github.com/jwohlwend/boltz), GitHub (Boltz developers)
15. [Meet Willow, our state-of-the-art quantum chip](https://blog.google/technology/research/google-willow-quantum-chip/), Google
16. [Our Quantum Echoes algorithm is a big step toward real-world applications for quantum computing](https://blog.google/technology/research/quantum-echoes-willow-verifiable-quantum-advantage/), Google
17. [IBM lays out clear path to fault-tolerant quantum computing](https://www.ibm.com/quantum/blog/large-scale-ftqc), IBM