417 |
|
average will match the ensemble average, therefore two similar |
418 |
|
trajectories in phase space should give matching statistical averages. |
419 |
|
|
420 |
< |
\subsection{\label{introSec:MDfurtheeeeer}Further Considerations} |
420 |
> |
\subsection{\label{introSec:MDfurther}Further Considerations} |
421 |
|
In the simulations presented in this research, a few additional |
422 |
|
parameters are needed to describe the motions. The simulations |
423 |
|
involving water and phospholipids in Chapt.~\ref{chaptLipids} are |
428 |
|
$z$-axis, and $\theta$ is a rotation about the new $x$-axis, and |
429 |
|
$\psi$ is a final rotation about the new $z$-axis (see |
430 |
|
Fig.~\ref{introFig:euleerAngles}). This sequence of rotations can be |
431 |
< |
accumulated into a single $3\time3$ matrix $\underline{\mathbf{A}}$ |
431 |
> |
accumulated into a single $3 \times 3$ matrix $\mathbf{A}$ |
432 |
|
defined as follows: |
433 |
|
\begin{equation} |
434 |
|
eq here |
443 |
|
\end{equation} |
444 |
|
Where $\omega^s_i$ is the angular velocity in the lab space frame |
445 |
|
along cartesian coordinate $i$. However, a difficulty arises when |
446 |
< |
attempting to integrate Eq.~\ref{introEq:MDeuleerPhi} and |
446 |
> |
attempting to integrate Eq.~\ref{introEq:MDeulerPhi} and |
447 |
|
Eq.~\ref{introEq:MDeulerPsi}. The $\frac{1}{\sin \theta}$ present in |
448 |
|
both equations means there is a non-physical instability present when |
449 |
|
$\theta$ is 0 or $\pi$. |
450 |
|
|
451 |
|
To correct for this, the simulations integrate the rotation matrix, |
452 |
< |
$\underline{\mathbf{A}}$, directly, thus avoiding the instability. |
452 |
> |
$\mathbf{A}$, directly, thus avoiding the instability. |
453 |
|
This method was proposed by Dullwebber |
454 |
|
\emph{et. al.}\cite{Dullwebber:1997}, and is presented in |
455 |
|
Sec.~\ref{introSec:MDsymplecticRot}. |