58 research outputs found
A rigorous sequential update strategy for parallel kinetic Monte Carlo simulation
The kinetic Monte Carlo (kMC) method is used in many scientific fields in
applications involving rare-event transitions. Due to its discrete stochastic
nature, efforts to parallelize kMC approaches often produce unbalanced time
evolutions requiring complex implementations to ensure correct statistics. In
the context of parallel kMC, the sequential update technique has shown promise
by generating high quality distributions with high relative efficiencies for
short-range systems. In this work, we provide an extension of the sequential
update method in a parallel context that rigorously obeys detailed balance,
which guarantees exact equilibrium statistics for all parallelization settings.
Our approach also preserves nonequilibrium dynamics with minimal error for many
parallelization settings, and can be used to achieve highly precise sampling
CLAIMED -- the open source framework for building coarse-grained operators for accelerated discovery in science
In modern data-driven science, reproducibility and reusability are key
challenges. Scientists are well skilled in the process from data to
publication. Although some publication channels require source code and data to
be made accessible, rerunning and verifying experiments is usually hard due to
a lack of standards. Therefore, reusing existing scientific data processing
code from state-of-the-art research is hard as well. This is why we introduce
CLAIMED, which has a proven track record in scientific research for addressing
the repeatability and reusability issues in modern data-driven science. CLAIMED
is a framework to build reusable operators and scalable scientific workflows by
supporting the scientist to draw from previous work by re-composing workflows
from existing libraries of coarse-grained scientific operators. Although
various implementations exist, CLAIMED is programming language, scientific
library, and execution environment agnostic.Comment: Received IEEE OSS Award 2023 -
https://conferences.computer.org/services/2023/symposia/oss.htm
Nonequilibrium candidate Monte Carlo: A new tool for efficient equilibrium simulation
Metropolis Monte Carlo simulation is a powerful tool for studying the
equilibrium properties of matter. In complex condensed-phase systems, however,
it is difficult to design Monte Carlo moves with high acceptance probabilities
that also rapidly sample uncorrelated configurations. Here, we introduce a new
class of moves based on nonequilibrium dynamics: candidate configurations are
generated through a finite-time process in which a system is actively driven
out of equilibrium, and accepted with criteria that preserve the equilibrium
distribution. The acceptance rule is similar to the Metropolis acceptance
probability, but related to the nonequilibrium work rather than the
instantaneous energy difference. Our method is applicable to sampling from both
a single thermodynamic state or a mixture of thermodynamic states, and allows
both coordinates and thermodynamic parameters to be driven in nonequilibrium
proposals. While generating finite-time switching trajectories incurs an
additional cost, driving some degrees of freedom while allowing others to
evolve naturally can lead to large enhancements in acceptance probabilities,
greatly reducing structural correlation times. Using nonequilibrium driven
processes vastly expands the repertoire of useful Monte Carlo proposals in
simulations of dense solvated systems
Monte Carlo Techniques for Drug Design: The Success Case of PELE
This chapter summarizes the most representative software packages that readily allow running Monte Carlo (MC) simulations in relevant scenarios for drug design. It explores in detail the Protein Energy Landscape Exploration (PELE) program, providing first the main characteristics of the technique, followed by an analysis of the different application studies in mapping protein‐ligand interactions. The ligand, formed by a rigid core and a set of rotatable side chains, is perturbed by translating and rotating it. PELE creates a list of perturbation poses, and then chooses the one with the lowest system energy. PELE was originally designed to map ligand migration pathways: its first applications aimed at finding exit pathways starting from ligand‐bound crystallographic structures. Additional applied studies have centered on modeling enzymatic mechanisms and engineering; the same techniques applied in mapping protein‐drug interactions can be used in the study of substrate recognition by enzymes.Along the development of PELE in the last years, we gratefully acknowledge financial support from the European Union (in particular from the ERC program) and from the Catalan and Spanish Governments. In addition we want to thank all present and past members from the EAPM lab. at BSC for their dedication and work.Peer ReviewedPostprint (author's final draft
Equilibrium molecular thermodynamics from Kirkwood sampling.
We present two methods for barrierless equilibrium sampling of molecular systems based on the recently proposed Kirkwood method (J. Chem. Phys. 2009, 130, 134102). Kirkwood sampling employs low-order correlations among internal coordinates of a molecule for random (or non-Markovian) sampling of the high dimensional conformational space. This is a geometrical sampling method independent of the potential energy surface. The first method is a variant of biased Monte Carlo, where Kirkwood sampling is used for generating trial Monte Carlo moves. Using this method, equilibrium distributions corresponding to different temperatures and potential energy functions can be generated from a given set of low-order correlations. Since Kirkwood samples are generated independently, this method is ideally suited for massively parallel distributed computing. The second approach is a variant of reservoir replica exchange, where Kirkwood sampling is used to construct a reservoir of conformations, which exchanges conformations with the replicas performing equilibrium sampling corresponding to different thermodynamic states. Coupling with the Kirkwood reservoir enhances sampling by facilitating global jumps in the conformational space. The efficiency of both methods depends on the overlap of the Kirkwood distribution with the target equilibrium distribution. We present proof-of-concept results for a model nine-atom linear molecule and alanine dipeptide.This research was funded by the European Research Council
and EPSRC grant EP/I001352/1. Y.O. was supported, in part,
by the JSPS Grant-in-Aid for Scientific Research on Innovative
Areas (“Dynamical Ordering and Integrated Functions”).This is the final published version. It first appeared at http://pubs.acs.org/doi/abs/10.1021/acs.jpcb.5b01800
Computational Identification of Protein Catalytic Sites: Tests, Validation
This project is one element of the analysis “pipeline” to characterize an organism that previously has not been well-studied. Once a protein of unknown structure has been computationally modeled (based on its sequence similarity to proteins with solved structures), then catalytic sites are identified on the model by comparison to a library of known sites. This work tested the identification algorithms with a set of proteins that have known structures and catalytic sites
Recommended from our members
Nonequilibrium candidate Monte Carlo: A new tool for efficient equilibrium simulation
Metropolis Monte Carlo simulation is a powerful tool for studying the equilibrium properties of matter. In complex condensed-phase systems, however, it is difficult to design Monte Carlo moves with high acceptance probabilities that also rapidly sample uncorrelated configurations. Here, we introduce a new class of moves based on nonequilibrium dynamics: candidate configurations are generated through a finite-time process in which a system is actively driven out of equilibrium, and accepted with criteria that preserve the equilibrium distribution. The acceptance rule is similar to the Metropolis acceptance probability, but related to the nonequilibrium work rather than the instantaneous energy difference. Our method is applicable to sampling from both a single thermodynamic state or a mixture of thermodynamic states, and allows both coordinates and thermodynamic parameters to be driven in nonequilibrium proposals. While generating finite-time switching trajectories incurs an additional cost, driving some degrees of freedom while allowing others to evolve naturally can lead to large enhancements in acceptance probabilities, greatly reducing structural correlation times. Using nonequilibrium driven processes vastly expands the repertoire of useful Monte Carlo proposals in simulations of dense solvated systems
Recommended from our members
Nonequilibrium candidate Monte Carlo: A new tool for efficient equilibrium simulation
Metropolis Monte Carlo simulation is a powerful tool for studying the equilibrium properties of matter. In complex condensed-phase systems, however, it is difficult to design Monte Carlo moves with high acceptance probabilities that also rapidly sample uncorrelated configurations. Here, we introduce a new class of moves based on nonequilibrium dynamics: candidate configurations are generated through a finite-time process in which a system is actively driven out of equilibrium, and accepted with criteria that preserve the equilibrium distribution. The acceptance rule is similar to the Metropolis acceptance probability, but related to the nonequilibrium work rather than the instantaneous energy difference. Our method is applicable to sampling from both a single thermodynamic state or a mixture of thermodynamic states, and allows both coordinates and thermodynamic parameters to be driven in nonequilibrium proposals. While generating finite-time switching trajectories incurs an additional cost, driving some degrees of freedom while allowing others to evolve naturally can lead to large enhancements in acceptance probabilities, greatly reducing structural correlation times. Using nonequilibrium driven processes vastly expands the repertoire of useful Monte Carlo proposals in simulations of dense solvated systems
- …