|
1002 | 1002 | "We find that the committor is constant within the metastable sets defined above.\n", |
1003 | 1003 | "Transition regions can be identified by committor values $\\approx 0.5$.\n", |
1004 | 1004 | "\n", |
1005 | | - "## Computing experimental observables [➜ 📓](07-expectations-and-observables.ipynb)\n", |
| 1005 | + "## Computing experimental observables [➜ 📓](06-expectations-and-observables.ipynb)\n", |
1006 | 1006 | "\n", |
1007 | 1007 | "Having thoroughly constructed, validated, and analyzed our MSM,\n", |
1008 | 1008 | "we may want to take the next step and compare our model to experimental data.\n", |
1009 | 1009 | "PyEMMA enables computation of stationary as well as dynamic experimental observables;\n", |
1010 | 1010 | "below we give give some examples of this.\n", |
1011 | 1011 | "We will make use of some external library functionality provided by MDTraj <a id=\"ref-9\" href=\"#cite-mdtraj\">mcgibbon-15</a>.\n", |
1012 | 1012 | "\n", |
1013 | | - "Theoretical backgrounds can be found here:\n", |
| 1013 | + "[Notebook 06](06-expectations-and-observables.ipynb) includes a brief summary of the theory used in below.\n", |
| 1014 | + "More in-depth descriptions of the theory and their applications to various data can be found in the following references:\n", |
1014 | 1015 | "- <a id=\"ref-10\" href=\"#cite-simon-amm\">olsson-17</a>\n", |
1015 | 1016 | "- <a id=\"ref-11\" href=\"#cite-noe-fingerprints\">noe-11</a>\n", |
1016 | 1017 | "- <a id=\"ref-12\" href=\"#cite-simon-mech-mod-nmr\">olsson-16</a>\n", |
|
1182 | 1183 | "metadata": {}, |
1183 | 1184 | "source": [ |
1184 | 1185 | "As for the stationary expectation (ensemble averages) considered above,\n", |
1185 | | - "we can use our pre-computed SASA vector to compute the auto-correlation function of trypotophan flourescene using the MSM `correlation()` method:" |
| 1186 | + "we can use our pre-computed SASA vector to compute the auto-correlation function of trypotophan flourescene using the MSM `correlation()` method (See [Notebook 06](06-expectations-and-observables.ipynb#Dynamic/kinetic-experimental-observables) for details):" |
1186 | 1187 | ] |
1187 | 1188 | }, |
1188 | 1189 | { |
|
1234 | 1235 | "Note the scale on the $y$-axis: this amplitude is likely too small to be experimentally measurable considering experimental uncertainty.\n", |
1235 | 1236 | "\n", |
1236 | 1237 | "However, using more advanced experimental setups such as stopped flow, T-jump, P-jump, and others,\n", |
1237 | | - "we can prepare our ensemble in a non-equilibrium initial condition.\n", |
| 1238 | + "we can prepare our ensemble in a non-equilibrium initial condition. (See [Notebook 06](06-expectations-and-observables.ipynb#Dynamic/kinetic-experimental-observables) for details)\n", |
1238 | 1239 | "\n", |
1239 | 1240 | "Let us say we can experimentally prepare a sample to only be in metastable state $\\mathcal{S}_1$.\n", |
1240 | 1241 | "In this case, the initial condition will be given by the metastable distribution of metastable state $\\mathcal{S}_1$, $p_0$.\n", |
|
1322 | 1323 | "\n", |
1323 | 1324 | "## Hidden Markov models [➜ 📓](07-hidden-markov-state-models.ipynb)\n", |
1324 | 1325 | "\n", |
1325 | | - "Another way of approaching metastable dynamics is with hidden Markov models (HMMs) <a id=\"ref-14\" href=\"#cite-hmm-baum-welch-alg\">baum-1970</a>.\n", |
| 1326 | + "Another way of approaching metastable dynamics is with hidden Markov models (HMMs) <a id=\"ref-14\" href=\"#cite-noe-proj-hid-msm\">noe-15</a>.\n", |
1326 | 1327 | "HMMs model the dynamics between so-called hidden states which we initiate from the metastable states found by PCCA++.\n", |
1327 | 1328 | "The estimation is less prone to discretization errors as we do not assume Markovianity in the space of our cluster centers.\n", |
1328 | 1329 | "It further provides a natural coarse graining into a given number of hidden states,\n", |
|
1757 | 1758 | "\n", |
1758 | 1759 | "<a id=\"cite-noe-dy-neut-scatt\"/><sup><a href=#ref-13>[^]</a></sup>Benjamin Lindner and Zheng Yi and Jan-Hendrik Prinz and Jeremy C. Smith and Frank Noé. 2013. _Dynamic neutron scattering from conformational dynamics. I. Theory and Markov models_. [URL](https://doi.org/10.1063/1.4824070)\n", |
1759 | 1760 | "\n", |
1760 | | - "<a id=\"cite-hmm-baum-welch-alg\"/><sup><a href=#ref-14>[^]</a></sup>Leonard E. Baum and Ted Petrie and George Soules and Norman Weiss. 1970. _A Maximization Technique Occurring in the Statistical Analysis of Probabilistic Functions of Markov Chains_. [URL](http://www.jstor.org/stable/2239727)\n", |
| 1761 | + "<a id=\"cite-noe-proj-hid-msm\"/><sup><a href=#ref-14>[^]</a></sup>Frank Noé and Hao Wu and Jan-Hendrik Prinz and Nuria Plattner. 2013. _Projected and hidden Markov models for calculating kinetics and metastable states of complex molecules_. [URL](https://doi.org/10.1063/1.4828816)\n", |
1761 | 1762 | "\n" |
1762 | 1763 | ] |
1763 | 1764 | } |
|
1778 | 1779 | "name": "python", |
1779 | 1780 | "nbconvert_exporter": "python", |
1780 | 1781 | "pygments_lexer": "ipython3", |
1781 | | - "version": "3.6.5" |
| 1782 | + "version": "3.6.6" |
1782 | 1783 | }, |
1783 | 1784 | "toc": { |
1784 | 1785 | "base_numbering": 1, |
|
0 commit comments