Skip to content

Commit 8e84729

Browse files
committed
Add LiveCoMS cite for Grossfield et al.
1 parent e47cbe8 commit 8e84729

2 files changed

Lines changed: 21 additions & 5 deletions

File tree

paper/basic_training.bib

Lines changed: 16 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1155,3 +1155,19 @@ @article{Hirel:2015:ComputerPhysicsCommunications
11551155
keywords = {Atomistic simulations,Dislocation,File conversion,Nye tensor,Polycrystal},
11561156
pages = {212-219},
11571157
}
1158+
1159+
1160+
@article{Grossfield:2019:LivingJ.Comput.Mol.Sci.,
1161+
title = {Best {{Practices}} for {{Quantification}} of {{Uncertainty}} and {{Sampling Quality}} in {{Molecular Simulations}} [{{Article}} v1.0]},
1162+
volume = {1},
1163+
doi = {10.33011/livecoms.1.1.5067},
1164+
language = {en},
1165+
number = {1},
1166+
journal = {Living Journal of Computational Molecular Science},
1167+
author = {Grossfield, Alan and Patrone, Paul N and Roe, Daniel R and Schultz, Andrew J and Siderius, Daniel W and Zuckerman, Daniel M},
1168+
year = {2019},
1169+
pages = {24},
1170+
}
1171+
1172+
1173+

paper/basic_training.tex

Lines changed: 5 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -700,12 +700,12 @@ \subsubsection{Production}
700700
This is closely related to the above discussion of equilibration.
701701
Depending on the relaxation timescales involved, one may realize only after analysis of a ``production'' trajectory that the system was still equilibrating in some sense.
702702

703-
A separate Best Practices document addresses the critical issues of convergence and error analysis; we refer the reader there for more details (\url{https://github.com/dmzuckerman/Sampling-Uncertainty}).
703+
A separate Best Practices document addresses the critical issues of convergence and error analysis; we refer the reader there for more details~\cite{Grossfield:2019:LivingJ.Comput.Mol.Sci.} (\url{https://github.com/dmzuckerman/Sampling-Uncertainty}).
704704
For more specific details on procedures and parameters used in production simulations, see the appropriate best practices document for the system of interest.
705705

706706
One other key consideration in production is what data to store, and how often.
707707
Storing data especially frequently can be tempting, but utilizes a great deal of storage space and does not actually provide significant value in most situations.
708-
Particularly, observations made in MD simulations are correlated in time (e.g. see \url{https://github.com/dmzuckerman/Sampling-Uncertainty}) so storing data more frequently than the autocorrelation time results in storage of essentially redundant data.
708+
Particularly, observations made in MD simulations are correlated in time (e.g. see \url{https://github.com/dmzuckerman/Sampling-Uncertainty}~\cite{Grossfield:2019:LivingJ.Comput.Mol.Sci.}) so storing data more frequently than the autocorrelation time results in storage of essentially redundant data.
709709
Thus, storing data more frequently than intervals of the autocorrelation time is generally unnecessary.
710710
Of course, the autocorrelation time is not known \emph{a priori} which can make it necessary to store \emph{some} redundant data.
711711
Disk space may also be a limiting factor that dictates the frequency of storing data, and should at least be considered.
@@ -1208,7 +1208,7 @@ \section{Should you run MD?}
12081208
A critical question \emph{before} preparing an MD simulation of your system is whether you even \emph{should} use MD for your system in view of the resources you have and what information you hope to obtain.
12091209
MD is a tool, but it may not be the right tool for your problem.
12101210
Before beginning any study, it is critical to sort out what questions you want to answer, what resources (computational and otherwise) you have at your disposal, and whether you have any information about your system(s) of interest that indicate you can realistically expect to answer those questions given a set of MD simulations.
1211-
Try to understand basic concepts of statistical uncertainty (\cite{Grossfield:2009:AnnuRepComputChem} and \url{https://github.com/dmzuckerman/Sampling-Uncertainty}) and use these to make an educated guess regarding your chances of extracting pertinent and reliable information from your simulation.
1211+
Try to understand basic concepts of statistical uncertainty (\cite{Grossfield:2009:AnnuRepComputChem} and \url{https://github.com/dmzuckerman/Sampling-Uncertainty}~\cite{Grossfield:2019:LivingJ.Comput.Mol.Sci.}) and use these to make an educated guess regarding your chances of extracting pertinent and reliable information from your simulation.
12121212

12131213
As noted above, the frequency of the fastest vibrational motions in a system of interest sets a fundamental limit on the timestep which, given fixed computational resources, sets a limit on how much simulation time can be covered with any reasonable amount of computer time.
12141214
Thus, as noted in Section~\ref{sec:intro}, the longest all-atom MD simulations are on the microsecond to millisecond timescale.
@@ -1240,7 +1240,7 @@ \section{Use your MD simulations and interpret the results with care and caution
12401240
For example, even a very short, unequilibrated MD simulation can produce movies which appear interesting and, by virtue of the fact that they result from MD, reveal the positions of all the atoms in a system as a function of time.
12411241
It's easy to run several short MD simulations where (for example) the composition of the system is varied, and conclude that any observed differences are a result of variations in composition.
12421242
But as noted in Section~\ref{sec:velocities}, even simulations started from the \emph{same} structure but slightly different initial positions or velocities will diverge over time yielding different results, so perhaps any differences are simply a result of this divergence rather than due to the change in conditions.
1243-
Thus, analysis will require great care and caution to avoid overinterpreting data, and error analysis becomes particularly critical (as discussed in \url{https://github.com/dmzuckerman/Sampling-Uncertainty}).
1243+
Thus, analysis will require great care and caution to avoid overinterpreting data, and error analysis becomes particularly critical (as discussed in \url{https://github.com/dmzuckerman/Sampling-Uncertainty}~\cite{Grossfield:2019:LivingJ.Comput.Mol.Sci.}).
12441244

12451245
In summary, then, do not use MD simulations simply to make movies and inspect these. Considerable care must be exercised to avoid overinterpeting the full atomistic detail they provide.
12461246
While movies in some cases can be useful, proper error analysis is always essential.
@@ -1260,7 +1260,7 @@ \section{Conclusions}
12601260

12611261
Our focus here has been on the basics --- focusing on things you need to understand before beginning to prepare simulations for yourself.
12621262
Additionally, we have primarily focused on issues relating to how simulations are conducted, and leave data analysis for a separate treatment.
1263-
As a starting point relating to data analysis, readers should probably review the Best Practices document on sampling and uncertainty estimation (\url{https://github.com/dmzuckerman/Sampling-Uncertainty}).
1263+
As a starting point relating to data analysis, readers should probably review the Best Practices document on sampling and uncertainty estimation (\url{https://github.com/dmzuckerman/Sampling-Uncertainty}~\cite{Grossfield:2019:LivingJ.Comput.Mol.Sci.}).
12641264

12651265
Please remember that this is an updatable work, so we welcome contributions and suggestions via our GitHub issue tracker at \url{https://github.com/MobleyLab/basic_simulation_training}.
12661266

0 commit comments

Comments
 (0)