You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Copy file name to clipboardExpand all lines: paper/basic_training.tex
+12-12Lines changed: 12 additions & 12 deletions
Original file line number
Diff line number
Diff line change
@@ -244,7 +244,7 @@ \subsubsection{Books}
244
244
Depending on the background, the practitioner can choose one or more of the following books to either learn or refresh their basic knowledge of thermodynamics.
245
245
Here are some works we find particularly helpful:
246
246
\begin{itemize}
247
-
\item Atkins and De Paula's ``Physical Chemistry''\cite{AtkinsBook}, chapters 1 to 4.
247
+
\item Atkins and De Paula's ``Physical Chemistry''~\cite{AtkinsBook}, chapters 1 to 4.
248
248
\item McQuarrie and Simon's extensive work, ``Physical Chemistry: A Molecular Approach''~\cite{McQuarrie:1997:}
\item Kittel and Kroemer's ``Thermal Physics''~\cite{Kittel:1980:}
@@ -728,7 +728,7 @@ \subsubsection{Background and How They Work}
728
728
729
729
Thermostat algorithms work by altering the Newtonian equations of motion that are inherently microcanonical (constant energy).
730
730
Thus, it is preferable that a thermostat not be used if it is desired to calculate dynamical properties such as diffusion coefficients; instead, the thermostat should be turned off after equilibrating the system to the desired temperature.
731
-
However, while all thermostats give non-physical dynamics, some have been found to have little effect on the calculation of particular dynamical properties, and they are commonly used during the production simulation as well\cite{Basconi:2013:JChemTheoryComput}.
731
+
However, while all thermostats give non-physical dynamics, some have been found to have little effect on the calculation of particular dynamical properties, and they are commonly used during the production simulation as well~\cite{Basconi:2013:JChemTheoryComput}.
732
732
733
733
There are several ways to categorize the many thermostatting algorithms that have been developed.
734
734
For example, thermostats can be either deterministic or stochastic depending on whether they use random numbers to guide the dynamics, and they can be either global or local depending on whether they are coupled to the dynamics of the full system or of a small subset.
The simple velocity rescaling thermostat is one of the easiest thermostats to implement; however, this thermostat is also one of the most non-physical thermostats.
759
759
This thermostat relies on rescaling the momenta of the particles such that the simulation's instantaneous temperature exactly matches the target temperature~\cite{thermostatAlgorithms2005}.
760
760
Similarly to the Gaussian thermosat, simple velocity rescaling aims to sample the isokinetic ensemble rather than the canonical ensemble.
761
-
However, it has been shown that the simple velocity rescaling fails to properly sample the isokinetic ensemble except in the limit of extremely small timesteps\cite{Braun:2018}.
762
-
Its usage can lead to simulation artifacts, so it is not recommended\cite{Harvey:1998:JCompChem,Braun:2018}.
761
+
However, it has been shown that the simple velocity rescaling fails to properly sample the isokinetic ensemble except in the limit of extremely small timesteps~\cite{Braun:2018}.
762
+
Its usage can lead to simulation artifacts, so it is not recommended~\cite{Harvey:1998:JCompChem,Braun:2018}.
763
763
764
764
\item\textbf{Berendsen}
765
765
766
-
The Berendsen\cite{berendsen1984molecular} thermostat (also known as the weak coupling thermostat) is similar to the simple velocity rescaling thermostat, but instead of rescaling velocities completely and abruptly to the target kinetic energy, it includes a relaxation term to allow the system to more slowly approach the target.
766
+
The Berendsen~\cite{berendsen1984molecular} thermostat (also known as the weak coupling thermostat) is similar to the simple velocity rescaling thermostat, but instead of rescaling velocities completely and abruptly to the target kinetic energy, it includes a relaxation term to allow the system to more slowly approach the target.
767
767
Although the Berendsen thermostat allows for temperature fluctuations, it samples neither the canonical distribution nor the isokinetic distribution.
768
-
Its usage can lead to simulation artifacts, so it is not recommended\cite{Harvey:1998:JCompChem,Braun:2018}.
768
+
Its usage can lead to simulation artifacts, so it is not recommended~\cite{Harvey:1998:JCompChem,Braun:2018}.
769
769
770
770
\item\textbf{Bussi-Donadio-Parrinello (Canonical Sampling through Velocity Rescaling)}
771
771
772
-
The Bussi\cite{Bussi:2007:JChemPhys:Canonical} thermostat is similar to the simple velocity rescaling and Berendsen thermostats, but instead of rescaling to a single kinetic energy that corresponds to the target temperature, the rescaling is done to a kinetic energy that is stochastically chosen from the kinetic energy distribution dictated by the canonical ensemble.
772
+
The Bussi~\cite{Bussi:2007:JChemPhys:Canonical} thermostat is similar to the simple velocity rescaling and Berendsen thermostats, but instead of rescaling to a single kinetic energy that corresponds to the target temperature, the rescaling is done to a kinetic energy that is stochastically chosen from the kinetic energy distribution dictated by the canonical ensemble.
773
773
Thus, this thermostat properly samples the canonical ensemble.
774
774
Similarly to the Berendsen thermostat, a user-specified time coupling parameter can be chosen to vary how abruptly the velocity rescaling takes place
775
-
The choice of time coupling constant does not affect structural properties, and most dynamical properties are fairly independent of the coupling constant within a broad range\cite{Bussi:2007:JChemPhys:Canonical}.
775
+
The choice of time coupling constant does not affect structural properties, and most dynamical properties are fairly independent of the coupling constant within a broad range~\cite{Bussi:2007:JChemPhys:Canonical}.
776
776
777
777
\item\textbf{Andersen}
778
778
779
-
The Andersen\cite{andersen1980molecular} thermostat works by selecting particles at random and having them ``collide'' with a heat bath by giving the particle a new velocity sampled from the Maxwell-Boltzmann distribution.
779
+
The Andersen~\cite{andersen1980molecular} thermostat works by selecting particles at random and having them ``collide'' with a heat bath by giving the particle a new velocity sampled from the Maxwell-Boltzmann distribution.
780
780
The number of particles affected, the time between ``collisions'', and how often it is applied to the system are possible variations of this thermostat.
781
781
The Andersen thermostat does reproduce the canonical ensemble.
782
782
However, it should only be used to sample structural properties, as dynamical properties can be greatly affected by the abrupt collisions.
783
783
784
784
\item\textbf{Langevin}
785
785
786
-
The Langevin\cite{schneider1978molecular} thermostat supplements the microcanonical equations of motion with Brownian dynamics, thus including the viscosity and random collision effects of an implicit solvent.
786
+
The Langevin~\cite{schneider1978molecular} thermostat supplements the microcanonical equations of motion with Brownian dynamics, thus including the viscosity and random collision effects of an implicit solvent.
787
787
It uses a general equation of the form $F = F_{interaction} + F_{friction} + F_{random}$, where $F_{interaction}$ is the standard interactions calculated during the simulation, $F_{friction}$ is the damping used to tune the ``viscosity'' of the implicit bath, and $F_{random}$ effectively gives random collisions with solvent molecules.
788
788
The frictional and random forces are coupled through a user-specified friction damping parameter. Careful consideration must be taken when choosing this parameter; in the limit of a zero damping parameter, both frictional and random forces go to zero and the dynamics become microcanonical, and in the limit of an infinite damping parameter, the dynamics are purely Brownian.
The choice of ``mass'' of the fictitious particle (which in many simulation packages is instead expressed as a time damping parameter) can be important as it affects the fluctuations that will be observed.
795
795
For many reasonable choices of the mass, dynamics are well-preserved~\cite{Basconi:2013:JChemTheoryComput}.
796
796
This is one of the most widely implemented and used thermostats.
797
-
However, it should be noted that with small systems, ergodicity can be an issue\cite{martyna1992nose,thermostatAlgorithms2005}.
797
+
However, it should be noted that with small systems, ergodicity can be an issue~\cite{martyna1992nose,thermostatAlgorithms2005}.
798
798
This can become important even in systems with larger numbers of particles if a portion of the system does not interact strongly with the remainder of the system, such as in alchemical free energy calculations when a solute or ligand is non-interacting.
799
799
Martyna et al.~\cite{martyna1992nose} discovered that by chaining thermostats, ergodicity can be enhanced, and most implementations of this thermostat use Nos\'{e}-Hoover chains.
800
800
@@ -848,7 +848,7 @@ \subsubsection{Background and How They Work}
848
848
To sample from the isothermal-isobaric ensemble (NPT), a thermostating algorithm like the ones discussed earlier must also be applied.
849
849
850
850
Much of the background information on barostats is analogous to thermostats.
851
-
The pressure of a molecular dynamics simulation is commonly measured using the virial theorem (an expectation value relating to positions and forces)\cite{ShellNotes, LeachBook}.
851
+
The pressure of a molecular dynamics simulation is commonly measured using the virial theorem (an expectation value relating to positions and forces)~\cite{ShellNotes, LeachBook}.
852
852
When pairwise interactions and periodic boundary conditions are considered, different approaches are often utilized~\cite{allenTildesleyLiquids, tuckermanBook, ShellNotes}.
853
853
Regardless, these formulas give pressure as a time-averaged quantity, similar to the temperature.
854
854
If we use these formulas to calculate the pressure for a single snapshot, this quantity is referred to as the instantaneous pressure.
0 commit comments