Related papers: Ensemble Optimization Techniques for the Simulatio…
We propose a Monte Carlo algorithm designed to simulate quantum as well as classical systems at equilibrium, bridging the algorithmic gap between quantum and classical thermal simulation algorithms. The method is based on a novel…
In the absence of impurities and boundary effects, first order phase transitions are initiated by the nucleation of critical bubbles. In thermally driven transitions many systems can remain metastable for an extended time, possibly tens of…
We introduce a general Monte Carlo scheme for achieving atomistic simulations with monoelectronic Hamiltonians including the thermalization of both nuclear and electronic degrees of freedom. The kinetic Monte Carlo algorithm is used to…
Parallel tempering and population annealing are both effective methods for simulating equilibrium systems with rough free energy landscapes. Parallel tempering, also known as replica exchange Monte Carlo, is a Markov chain Monte Carlo…
We apply the recently developed adaptive ensemble optimization technique to simulate dense Lennard-Jones fluids and a particle-solvent model by broad-histogram Monte Carlo techniques. Equilibration of the simulated fluid is improved by…
Calculations of topological observables in lattice gauge theories with traditional Monte Carlo algorithms have long been known to be a difficult task, owing to the effects of long autocorrelations times. Several mitigation strategies have…
We consider several multiscale-in-time kinetic Monte Carlo models, in which some variables evolve on a fast time scale, while the others evolve on a slow time scale. In the first two models we consider, a particle evolves in a…
We present a novel Exchange Monte Carlo (EMC) method designed for application in continuous-space Path Integral Monte Carlo (PIMC) simulations at finite temperature. Traditional PIMC methods for bosonic systems suffer from long…
We describe parallel Markov chain Monte Carlo methods that propagate a collective ensemble of paths, with local covariance information calculated from neighboring replicas. The use of collective dynamics eliminates multiplicative noise and…
The most efficient MC weights for the calculation of physical, canonical expectation values are not necessarily those of the canonical ensemble. The use of suitably generalized ensembles can lead to a much faster convergence of the…
The efficiency of statistical sampling in broad-histogram Monte Carlo simulations can be considerably improved by optimizing the simulated extended ensemble for fastest equilibration. Here we describe how a recently developed feedback…
In biomolecular systems (especially all-atom models) with many degrees of freedom such as proteins and nucleic acids, there exist an astronomically large number of local-minimum-energy states. Conventional simulations in the canonical…
It was recently demonstrated that a simple Monte Carlo (MC) algorithm involving the swap of particle pairs dramatically accelerates the equilibrium sampling of simulated supercooled liquids. We propose two numerical schemes integrating the…
We show how to use the multiple histogram method to combine canonical ensemble Monte Carlo simulations made at different temperatures and densities. The method can be applied to study systems of particles with arbitrary interaction…
The recently developed density matrix quantum Monte Carlo (DMQMC) algorithm stochastically samples the N -body thermal density matrix and hence provides access to exact properties of many-particle quantum systems at arbitrary temperatures.…
The importance-sampling Monte Carlo algorithm appears to be the universally optimal solution to the problem of sampling the state space of statistical mechanical systems according to the relative importance of configurations for the…
In complex systems with many degrees of freedom such as peptides and proteins there exist a huge number of local-minimum-energy states. Conventional simulations in the canonical ensemble are of little use, because they tend to get trapped…
We present a study of the parallel tempering (replica exchange) Monte Carlo method, with special focus on the feedback-optimized parallel tempering algorithm, used for generating an optimal set of simulation temperatures. This method is…
We introduce and implement a Monte Carlo scheme to study the equilibrium statistics of polymers in the globular phase. It is based on a model of "interacting elastic lattice polymers" and allows a sufficiently good sampling of long and…
Techniques for simulating molecules whose conformations satisfy constraints are presented. A method for selecting appropriate moves in Monte Carlo simulations is given. The resulting moves not only obey the constraints but also maintain…