English

Enhanced Diffusion Sampling: Efficient Rare Event Sampling and Free Energy Calculation with Diffusion Models

Machine Learning 2026-02-19 v1 Artificial Intelligence Machine Learning Biological Physics Chemical Physics

Abstract

The rare-event sampling problem has long been the central limiting factor in molecular dynamics (MD), especially in biomolecular simulation. Recently, diffusion models such as BioEmu have emerged as powerful equilibrium samplers that generate independent samples from complex molecular distributions, eliminating the cost of sampling rare transition events. However, a sampling problem remains when computing observables that rely on states which are rare in equilibrium, for example folding free energies. Here, we introduce enhanced diffusion sampling, enabling efficient exploration of rare-event regions while preserving unbiased thermodynamic estimators. The key idea is to perform quantitatively accurate steering protocols to generate biased ensembles and subsequently recover equilibrium statistics via exact reweighting. We instantiate our framework in three algorithms: UmbrellaDiff (umbrella sampling with diffusion models), Δ\DeltaG-Diff (free-energy differences via tilted ensembles), and MetaDiff (a batchwise analogue for metadynamics). Across toy systems, protein folding landscapes and folding free energies, our methods achieve fast, accurate, and scalable estimation of equilibrium properties within GPU-minutes to hours per system -- closing the rare-event sampling gap that remained after the advent of diffusion-model equilibrium samplers.

Keywords

Cite

@article{arxiv.2602.16634,
  title  = {Enhanced Diffusion Sampling: Efficient Rare Event Sampling and Free Energy Calculation with Diffusion Models},
  author = {Yu Xie and Ludwig Winkler and Lixin Sun and Sarah Lewis and Adam E. Foster and José Jiménez Luna and Tim Hempel and Michael Gastegger and Yaoyi Chen and Iryna Zaporozhets and Cecilia Clementi and Christopher M. Bishop and Frank Noé},
  journal= {arXiv preprint arXiv:2602.16634},
  year   = {2026}
}
R2 v1 2026-07-01T10:41:38.753Z