Simulated tempering with irreversible Gibbs sampling techniques
Fahim Faizi1, Pedro J Buigues2, George Deligiannidis3
1Department of Mathematics, King's College London, Strand, WC2R 2LS London, United Kingdom.
The Journal of Chemical Physics
|December 9, 2020
Summary
We developed new irreversible algorithms for simulated tempering simulations that improve sampling efficiency. These methods break detailed balance but maintain target distribution invariance, outperforming conventional techniques.
Area of Science:
- Computational Physics
- Statistical Mechanics
- Chemical Physics
Background:
- Simulated tempering (ST) is a powerful Monte Carlo method for exploring complex energy landscapes.
- Conventional ST methods often rely on detailed balance condition (DBC), which can limit efficiency.
- Breaking DBC while preserving target distribution invariance is a key challenge.
Purpose of the Study:
- To introduce two novel irreversible algorithms for simulated tempering.
- To demonstrate improved sampling efficiency compared to conventional methods.
- To validate the algorithms on diverse systems including Ising models and molecular dynamics.
Main Methods:
- Developed irreversible algorithms based on Gibbs sampling, specifically targeting temperature swap updates.
- Broke the detailed balance condition (DBC) while ensuring skewed detailed balance for target distribution invariance.
- Tested algorithms on a 1D double-well potential, the Ising model, and alanine pentapeptide (ALA5) molecular dynamics simulations.
Main Results:
- Achieved significant gains in sampling efficiency for the Ising model, evidenced by faster relaxation times for inverse temperature, magnetic susceptibility, and energy density.
- Observed improved mixing times for inverse temperature and system energy in ALA5 simulations with numerous temperature replicas compared to standard Metropolis-Hastings (MH).
- Demonstrated favorable constant scaling in Ising spin systems for large numbers of temperature ladders, outperforming both reversible and irreversible MH algorithms.
Conclusions:
- The novel irreversible algorithms offer a more efficient alternative to conventional simulated tempering methods with DBC, without additional computational overhead.
- These algorithms show particular promise for large-scale simulations with extensive temperature ranges.
- Future applications include tailoring these irreversible methods to other dynamical variables for flattening rugged free energy landscapes.
More Related Videos
Related Concept Videos
Temperature Dependence on Reaction Rate
86.6K
The Collision Theory
Atoms, molecules, or ions must collide before they can react with each other. Atoms must be close together to form chemical bonds. This premise is the basis for a theory that explains many observations regarding chemical kinetics, including factors affecting reaction rates.
The collision theory is based on the postulates that (i) the reaction rate is proportional to the rate of reactant collisions, (ii) the reacting species collide in an orientation allowing contact between...
Atoms, molecules, or ions must collide before they can react with each other. Atoms must be close together to form chemical bonds. This premise is the basis for a theory that explains many observations regarding chemical kinetics, including factors affecting reaction rates.
The collision theory is based on the postulates that (i) the reaction rate is proportional to the rate of reactant collisions, (ii) the reacting species collide in an orientation allowing contact between...
86.6K
Isothermal Processes
4.5K
A thermodynamic process that occurs at constant temperature is called an isothermal process. Heat slowly flows into the system or out of the system to maintain thermal equilibrium. Processes involving phase changes like water evaporation into steam or freezing water into ice at a constant temperature are examples of Isothermal Processes.
An ideal gas can also undergo isothermal expansion or compression.
For example, consider 1 mole of an ideal gas inside an isolated cylinder at initial volume V...
An ideal gas can also undergo isothermal expansion or compression.
For example, consider 1 mole of an ideal gas inside an isolated cylinder at initial volume V...
4.5K
Thermal Sigmatropic Reactions: Overview
2.3K
Sigmatropic rearrangements are a class of pericyclic reactions in which a σ bond migrates from one part of a π system to another. These are intramolecular rearrangements where the total number of σ and π bonds remain unchanged.
Sigmatropic shifts are classified based on an order term [i, j ], where i and j indicate the number of atoms across which each end of the σ bond migrates. Below are examples of a [3,3] sigmatropic shift in 1,5-hexadiene, referred...
Sigmatropic shifts are classified based on an order term [i, j ], where i and j indicate the number of atoms across which each end of the σ bond migrates. Below are examples of a [3,3] sigmatropic shift in 1,5-hexadiene, referred...
2.3K
Temperature and Thermal Equilibrium
8.4K
Heat and temperature are essential concepts for everyone every day. The study of heat and temperature is part of an area of physics known as thermodynamics. It is not always easy to distinguish heat and temperature.
The concept of temperature has evolved from the common concepts of hot and cold. The scientific definition of temperature explains more than just our sense of hot and cold. Temperature is operationally defined as the quantity measured with a thermometer. Furthermore, temperature is...
The concept of temperature has evolved from the common concepts of hot and cold. The scientific definition of temperature explains more than just our sense of hot and cold. Temperature is operationally defined as the quantity measured with a thermometer. Furthermore, temperature is...
8.4K
Constant Volume Calorimetry
29.5K
Calorimeters are useful to determine the heat released or absorbed by a chemical reaction. Coffee cup calorimeters are designed to operate at constant (atmospheric) pressure and are convenient to measure heat flow (or enthalpy change) accompanying processes that occur in solution at constant pressure. A different type of calorimeter that operates at constant volume, colloquially known as a bomb calorimeter, is used to measure the energy produced by reactions that yield large amounts of heat and...
29.5K
Gibbs Free Energy
36.7K
One of the challenges of using the second law of thermodynamics to determine if a process is spontaneous is that it requires measurements of the entropy change for the system and the entropy change for the surroundings. An alternative approach involving a new thermodynamic property defined in terms of system properties only was introduced in the late nineteenth century by American mathematician Josiah Willard Gibbs. This new property is called the Gibbs free energy (G) (or simply the free...
36.7K


