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: manuscript/manuscript.tex
+24-24Lines changed: 24 additions & 24 deletions
Original file line number
Diff line number
Diff line change
@@ -248,7 +248,7 @@ \subsection{Variational approach and TICA}
248
248
249
249
\subsection{Hidden Markov state models}
250
250
251
-
\begin{figure}
251
+
\begin{figure}[ht]
252
252
\includegraphics[width=0.48\textwidth]{figure_1}
253
253
\caption{The HMM transition matrix $\tilde{\mathbf{P}}(\tau)$ propagates the hidden state trajectory $\tilde{s}(t)$ (orange circles) and, at each time step $t$, the emission into the observable state $s(t)$ (cyan circles) is governed by the emission probabilities $\bm{\chi}\left( s(t) \middle| \tilde{s}(t) \right)$.}
254
254
\label{fig:hmm-scheme}
@@ -325,6 +325,15 @@ \section{PyEMMA tutorials}
325
325
326
326
\subsection{The PyEMMA workflow}
327
327
328
+
\begin{figure}[ht]
329
+
\includegraphics[width=0.48\textwidth]{figure_2}
330
+
\caption{The PyEMMA workflow: MD trajectories are processed and discretized (first row).
331
+
A Markov state model is estimated from the resulting discrete trajectories and validated (middle row).
332
+
By iterating between data processing and MSM estimation/validation,
333
+
a dynamical model is obtained that can be analyzed (last row).}
334
+
\label{fig:workflowchart}
335
+
\end{figure}
336
+
328
337
In short, the workflow (Fig.~\ref{fig:workflowchart}) for a full analysis of an MD dataset might consist of,
329
338
\begin{itemize}
330
339
\item extracting molecular features from the raw data (01),
We demonstrate how to apply TICA, suggest how to interpret the projected coordinates, and compare the results to other dimension reduction techniques in notebook~02.
402
412
403
-
\begin{figure}
404
-
\includegraphics{figure_3}
405
-
\caption{Example analysis of the conformational dynamics of a pentapeptide backbone:
406
-
(a)~The Trp-Leu-Ala-Leu-Leu pentapeptide in licorice representation~\cite{vmd}.
407
-
(b)~The VAMP-2 score indicates which of the tested featurizations contains the highest kinetic variance.
408
-
(c)~The sample free energy projected onto the first two time-lagged independent components (ICs) at lag time $\tau=0.5$~ns shows multiple minima and
409
-
(d)~the time series of the first two ICs of the first trajectory show rare jumps.}
410
-
\label{fig:io-to-tica}
411
-
\end{figure}
412
-
413
413
\subsection{Discretization}
414
414
415
415
TICA yields a representation of our molecular simulation data with a reduced dimensionality,
@@ -421,7 +421,7 @@ \subsection{Discretization}
421
421
422
422
\subsection{MSM estimation and validation}
423
423
424
-
\begin{figure}
424
+
\begin{figure}[ht]
425
425
\includegraphics{figure_4}
426
426
\caption{Example analysis of the conformational dynamics of a pentapeptide backbone:
427
427
(a)~The convergence behavior of the implied timescales associated with the four slowest processes.
@@ -456,7 +456,7 @@ \subsection{MSM estimation and validation}
456
456
457
457
\subsection{Analyzing the MSM}
458
458
459
-
\begin{figure}
459
+
\begin{figure}[ht]
460
460
\includegraphics{figure_5}
461
461
\caption{Example analysis of the conformational dynamics of a pentapeptide backbone:
462
462
(a)~The reweighted free energy surface projected onto the first two independent components exhibits five minima which
@@ -467,7 +467,7 @@ \subsection{Analyzing the MSM}
467
467
\label{fig:msm-analysis}
468
468
\end{figure}
469
469
470
-
\begin{figure}
470
+
\begin{figure}[ht]
471
471
\includegraphics{figure_6}
472
472
\caption{Example analysis of the conformational dynamics of a pentapeptide backbone:
473
473
visualization of the transition paths from $\mathcal{S}_2$ to $\mathcal{S}_4$.
@@ -537,7 +537,7 @@ \subsection{Analyzing the MSM}
537
537
538
538
\subsection{Connecting the MSM with experimental data}
539
539
540
-
\begin{figure}
540
+
\begin{figure}[ht]
541
541
\includegraphics{figure_7}
542
542
\caption{Example analysis of the conformational dynamics of a pentapeptide backbone:
543
543
(a)~the Trp-1 SASA autocorrelation function yields a weak signal which, however,
0 commit comments