Molecular simulation from modern statistics: Continuous-time, continuous-space, exact
arXiv:2305.02979 · doi:10.1038/s41598-024-66172-0
Abstract
In a world made of atoms, the computer simulation of molecular systems, such as proteins in water, plays an enormous role in science. Software packages that perform these computations have been developed for decades. In molecular simulation, Newton's equations of motion are discretized and long-range potentials are treated through cutoffs or spacial discretization, which all introduce approximations and artifacts that must be controlled algorithmically. Here, we introduce a paradigm for molecular simulation that is based on modern concepts in statistics and is rigorously free of discretizations, approximations, and cutoffs. Our demonstration software reaches a break-even point with traditional molecular simulation at high precision. We stress the promise of our paradigm as a gold standard for critical applications and as a future competitive approach to molecular simulation.
arXiv comment19 pages, 4 figures; 18 pages supplementary materials, 1 supplementary figure
References in corpus (20)
- Two-step melting in two dimensions: First-order liquid-hexatic transition
- How to mesh up Ewald sums (II): An accurate error estimate for the P3M algorithm
- Event-chain Monte Carlo algorithms for hard-sphere systems
- Generalized event-chain Monte Carlo: Constructing rejection-free global-balance algorithms from infinitesimal steps
- Rejection-free Monte-Carlo sampling for general potentials
- The good, the bad and the user in soft matter simulations
- Piecewise Deterministic Markov Processes for Continuous-Time Monte Carlo
- Event-chain Monte Carlo algorithms for three- and many-particle interactions
- Cell-veto Monte Carlo algorithm for long-range systems
- All-atom computations with irreversible Markov chains
- Efficient Equilibration of Hard Spheres with Newtonian Event Chains
- JeLLyFysh-Version1.0 -- a Python application for all-atom event-chain Monte Carlo
- Event-chain Monte Carlo with factor fields
- Multithreaded event-chain Monte Carlo with local times
- Hard-disk dipoles and non-reversible Markov chains
- Concepts in Monte Carlo sampling
- A Comparison of Popular Point Configurations on
- Sparse hard-disk packings and local Markov chains
- Direction-sweep Markov chains
- A new algorithm for electrostatic interactions in Monte Carlo simulations of charged particles