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
+5-7Lines changed: 5 additions & 7 deletions
Original file line number
Diff line number
Diff line change
@@ -176,7 +176,7 @@ \subsubsection{Key concepts}
176
176
177
177
Classical molecular models typically consist of point particles carrying mass and electric charge, as well as potentially additional interactions such as van der Waals interactions and bonded interactions of various types.
178
178
Sometimes it is much more efficient to freeze the internal degrees of freedoms and treat the molecule as a rigid body where the particles do not change their relative orientation as the whole body moves; this is commonly done, for example, for rigid models of the water molecule.
179
-
Due to the high frequency of the O-H vibrations, accurately treating water classically would require a very small time step, so for computational efficiency water is often instead treated as a rigid body.
179
+
Due to the high frequency of the O-H vibrations, accurately treating water classically would require a solving the equations of motion with a very small timestep, so for computational efficiency water is often instead treated as a rigid body.
180
180
Keeping specified objects rigid in a simulation involves applying holonomic constraints, where the rigidity is defined by imposing a minimal set of fixed bond lengths and angles through iterative procedures during the numerical integration of the equation of motion (see Section~\ref{sec:integrators} for more on constraints and integrators).
181
181
It is important to understand the concept of point particles, rigid bodies and constraints.
182
182
@@ -591,11 +591,11 @@ \subsection{Main steps of a molecular dynamics simulation}
591
591
Potentially the worst possible outcome is if the prepared system is not what you intended (e.g. it contains incorrect molecules or protonation states) but is chemically valid and well described by your force field and thus proceeds without error through the remaining steps --- and in fact this is a frequent outcome of problems in system preparation.
592
592
It should not be assumed that if a system can proceed in a well-behaved manner through the other steps, it was necessarily prepared correctly; considerable care should be taken here.
593
593
594
-
The purpose of minimization, or relaxation, is to find a local energy minimum of the starting structure so that the molecular dynamics simulation does not immediately "blow up" (i.e. the forces on any one atom are not so large that the atoms move an unreasonable distance in a single time step).
594
+
The purpose of minimization, or relaxation, is to find a local energy minimum of the starting structure so that the molecular dynamics simulation does not immediately "blow up" (i.e. the forces on any one atom are not so large that the atoms move an unreasonable distance in a single timestep).
595
595
This involves standard minimization algorithms such as steepest descent.
596
596
For a more involved discussion of minimization algorithms utilized in molecular simulation, see \citet{LeachBook}, sections 5.1-5.7.
597
597
598
-
At the end of energy minimization, it is important to achieve a system configuration with small enough forces that the desired time step will allow numerical integration of the equations of motion without overly large displacements (see \citet{LeachBook}, section 7.3.4).
598
+
At the end of energy minimization, it is important to achieve a system configuration with small enough forces that the desired timestep will allow numerical integration of the equations of motion without overly large displacements (see \citet{LeachBook}, section 7.3.4).
599
599
Such a configuration is a suitable starting point for molecular dynamics techniques.
600
600
However, this only represents a static set of positions, while the propagation of dynamics also requires a set of starting velocities.
601
601
These may be assigned in a variety of ways, but are usually randomly assigned to atoms such that the correct Maxwell-Boltzmann distribution at the desired temperature is achieved.
@@ -605,7 +605,7 @@ \subsection{Main steps of a molecular dynamics simulation}
605
605
For this reason, it is advised that a thermostatted simulation is performed prior to a desired production simulation, even if the production simulation will ultimately be done in the NVE ensemble.
606
606
Once the kinetic and potential energies fluctuate around constant values, the thermostat may be removed (if an NVE simulation is desired) and a snapshot selected that is simultaneously as close to the average kinetic and potential energies as possible.
607
607
This snapshot, containing both positions and velocities may be used to then start an NVE simulation that will correspond to a temperature close to that which is desired.
608
-
This is necessary due to the fact that only the average temperature is obtained through coupling to a thermostat (see Section~\ref{sec:thermostats}), and the temperature fluctuates with the kinetic energy at each time step.
608
+
This is necessary due to the fact that only the average temperature is obtained through coupling to a thermostat (see Section~\ref{sec:thermostats}), and the temperature fluctuates with the kinetic energy at each timestep.
609
609
Similarly, equilibration in the NPT ensemble is necessary before production in the NVT if an average density consistent with a specific pressure is desired.
610
610
In this case, the system may be scaled to the desired average volume before the production simulation.
611
611
In the above example, the NPT simulation is said to have equilibrated to a specific volume when the dimensions of the simulation box fluctuate around constant values with minimal drift.
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 as well.
681
-
This thermostat relies on rescaling the momenta of the particles every $N$time steps based on their ratio of the instantaneous temperature to the target temperature\cite{thermostatAlgorithms2005}.
681
+
This thermostat relies on rescaling the momenta of the particles every $N$timesteps based on their ratio of the instantaneous temperature to the target temperature\cite{thermostatAlgorithms2005}.
682
682
This leads to multiple issues when trying to sample data using this thermostat.
683
683
First, this method is not reversible, there lacks a way for the particles to have knowledge of their previous thermal history.
684
684
This makes any dynamical value impossible to measure (diffusion for example).
@@ -890,8 +890,6 @@ \subsubsection{Summary}
890
890
\subsection{Integrators}
891
891
\label{sec:integrators}
892
892
893
-
\todo[inline, color={yellow!20}]{DLM: Need to decide on ``time step'' versus ``timestep'' and change everywhere in whole paper; right now we use both.}
894
-
895
893
For systems consisting of more than three interacting bodies with no constrained degrees of freedom, there is no analytical solution to the equations of motion.
896
894
Instead, we must approximate the dynamics in a discrete manner.
897
895
This is usually termed numerical integration of the equations of motion.
0 commit comments