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: README.md
+5Lines changed: 5 additions & 0 deletions
Display the source diff
Display the rich diff
Original file line number
Diff line number
Diff line change
@@ -27,10 +27,15 @@ This paper is being developed as a living document, open to changes from the com
27
27
## List of Released Versions
28
28
<!-- update this when you decide to release a version either by preprint or when submitted to LiveCoMS-->
29
29
30
+
## Citing this work
31
+
32
+
This work has not yet been formally published at LiveCoMS; however, you can cite v1.0 of it using the DOI [10.33011/livecoms.1.1.5957](http://dx.doi.org/10.33011/livecoms.1.1.5957).
33
+
30
34
## Changelog
31
35
<!-- Here, record summaries of important changes. A granular discussion of changes will be kept in GitHub by issue tracking.-->
32
36
- Sept. 8, 2017: D. Mobley created repo, put in initial files from LiveCoMS.
33
37
- Spring/early summer 2018: Finalize who will be involved and write/edit first version of the paper
34
38
- Sept. 2, 2018: Final draft of version 1 submitted to LiveCoMS
35
39
- Nov. 5-6, 2018: Make editorial revisions suggested by peer reviewers and Victoria Lim.
36
40
- Nov. 23, 2018: Check references using [`fixbibtex`](https://github.com/jaimergp/fixbibtex), incorporate fixes for problems it caught; addresses a number of typos/missing references caught by Michael Shirts.
41
+
- Nov. 28, 2018: Add DOI, include citation to Grossfield et al. LiveCoMS article with DOI.
This is closely related to the above discussion of equilibration.
701
701
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.
702
702
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}).
704
704
For more specific details on procedures and parameters used in production simulations, see the appropriate best practices document for the system of interest.
705
705
706
706
One other key consideration in production is what data to store, and how often.
707
707
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.
709
709
Thus, storing data more frequently than intervals of the autocorrelation time is generally unnecessary.
710
710
Of course, the autocorrelation time is not known \emph{a priori} which can make it necessary to store \emph{some} redundant data.
711
711
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?}
1208
1208
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.
1209
1209
MD is a tool, but it may not be the right tool for your problem.
1210
1210
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.
1212
1212
1213
1213
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.
1214
1214
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
1240
1240
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.
1241
1241
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.
1242
1242
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.}).
1244
1244
1245
1245
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.
1246
1246
While movies in some cases can be useful, proper error analysis is always essential.
@@ -1260,7 +1260,7 @@ \section{Conclusions}
1260
1260
1261
1261
Our focus here has been on the basics --- focusing on things you need to understand before beginning to prepare simulations for yourself.
1262
1262
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.}).
1264
1264
1265
1265
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}.
0 commit comments