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
+9-5Lines changed: 9 additions & 5 deletions
Original file line number
Diff line number
Diff line change
@@ -310,6 +310,10 @@ \subsubsection{Key concepts}
310
310
In this type of non-equilibrium situation, the distribution of configurations will not follow a Boltzmann factor distribution.
311
311
Specialized simulation approaches are available to study such systems~\cite{Chong:2017:CurrentOpinioninStructuralBiology, Zuckerman:2017:AnnuRevBiophys} but they are not beginner-friendly.
312
312
Non-equilibrium molecular concepts pertinent to cell biology have been discussed at an introductory level (e.g. \url{http://www.physicallensonthecell.org/}).
313
+
Notably, many experiments are conducted at non-equilibrium conditions; for example, membrane diffusion coefficients are commonly measured by setting up a concentration gradient across the membrane and measuring the flux.
314
+
It can be tempting to the beginner to setup an MD simulation in the same manner as such an experiment.
315
+
However, maintaining non-equilibrium conditions is typically more complicated in an MD simulation than in an experiment as large reservoirs are commonly required.
316
+
Frequently there exist equilibrium methods that can be used to obtain the same or similar information as a non-equilibrium experiment; users should seek to obtain familiarity with such methods before choosing to conduct a non-equilibrium MD simulation.
313
317
314
318
315
319
\subsubsection{Books}
@@ -1006,7 +1010,7 @@ \subsubsection{Choosing an appropriate timestep}
1006
1010
For all-atom simulations with constraints on the high-frequency bonds, timesteps can be commonly increased to 2 fs; coarse-grained simulations with particles of higher mass and smaller force constants can have much larger timesteps.
1007
1011
After choosing a timestep, a test simulation should be run in the microcanonical ensemble to ensure that the choice of timestep yields dynamics that conserve energy.
1008
1012
Methods also exist to increase the timestep beyond the limit imposed by the system's highest-frequency motion.
1009
-
Many such enhanced timestepping algorithms exist such as multiple-timestep methods which separately integrate high-frequency and low-frequency motion and schemes which repartition atomic masses to decrease the highest-frequency motion seen in the system\cite{Berne:1999:Molecular,Hopkins:2015:JCTC:Long}.
1013
+
Some examples of these enhanced timestepping algorithms include multiple-timestep methods which separately integrate high-frequency motion from low-frequency motion and schemes which repartition atomic masses to decrease the highest-frequency motion seen in the system\cite{Berne:1999:Molecular,Hopkins:2015:JCTC:Long}.
However, here, because we use point charge electrostatics, $\rho$ is a set of delta functions.
1048
1052
1049
1053
The Ewald method is based on (temporarily) replacing the point charge distributions by smooth charge distributions in order to apply existing numerical techniques to solve this partial differential equation (PDE).
1050
-
The most common smooth function used in Ewald method is the Gaussian distribution, although other distributions have been used as well.
1054
+
The most common smooth function used in the Ewald method is the Gaussian distribution, although other distributions have been used as well.
1051
1055
Thus the overall charge distribution is divided into a short-range or ``direct space'' component ($\rho^{sr}$) involving the original point charges screened by the Gaussian-distributed charge of the same magnitude (Figure~\ref{fig:screening}) but opposite sign, and a long-range component involving Gaussian-distributed charges of the original sign ($\rho^{lr}$).
The potential due to long-range charge interactions does not decay rapidly, and thus requires consideration of all periodic copies.
1073
-
This would pose severe problems if calculated via direct summation, but the smoothness of the charge $\rho^{lr}$ (and hence potential ($\phi^{lr}$)) allows the use of fast PDE solvers.
1074
-
Specifically, while the sum is long-ranged in real space, taking the Fourier transform converts it into a sum in reciprocal space which is short-ranged in reciprocal space, damped by a factor $\exp{-k^2 \sigma^2/2}$ where $k$ is the reciprocal space vector and $\sigma$ is the width of the Gaussian.
1077
+
This would pose severe problems if calculated via direct summation, but the smoothness of the charge $\rho^{lr}$ (and hence potential ($\phi^{lr}$) allows the use of fast PDE solvers.
1078
+
Specifically, while the sum is long-ranged in real space, taking the Fourier transform converts it into a sum in reciprocal space which is short-ranged in reciprocal space, damped by a factor $\exp\left(-k^2\sigma^2/2\right)$ where $k$ is the reciprocal space vector and $\sigma$ is the width of the Gaussian.
1075
1079
1076
1080
The final term in Ewald summation is a so-called self term which gets subtracted out of the overall sum; it is calculated only once at the beginning of the simulation as it depends only on the charge magnitudes and not their positions.
1077
1081
It also does not contribute to the force.
@@ -1141,7 +1145,7 @@ \subsubsection{Grid based Ewald summation}
1141
1145
\begin{checklist}{Determine handling of cutoffs}
1142
1146
\begin{itemize}
1143
1147
\item As a general rule, electrostatics are long-range enough that either the cutoff needs to be larger than the system size (for finite systems) or
1144
-
periodicity is needed along with full treatment of long-range electrostatics (Section~\ref{sec:classical_electrostatics} % changed section ref since the folloing does not currently exist (Section~\ref{lr_electrostatics})
1148
+
periodicity is needed along with full treatment of long-range electrostatics (Section~\ref{sec:classical_electrostatics})% changed section ref since the folloing does not currently exist (Section~\ref{lr_electrostatics})
1145
1149
\item Nonpolar interactions can often be safely treated with cutoffs of 1-1.5 nm as long as the system size is at least twice that, but long-range dispersion corrections may be needed (Section~\ref{sec:force_fields})
0 commit comments