| 16 | 
  | 
 | 
| 17 | 
  | 
\frontmatter | 
| 18 | 
  | 
 | 
| 19 | 
< | 
\title{APPLICATION AND DEVELOPMENT OF MOLECULAR DYNAMICS TECHNIQUES FOR THE  | 
| 20 | 
< | 
STUDY OF WATER}      | 
| 19 | 
> | 
\title{DEVELOPMENT OF MOLECULAR DYNAMICS TECHNIQUES FOR THE  | 
| 20 | 
> | 
STUDY OF WATER AND BIOCHEMICAL SYSTEMS}      | 
| 21 | 
  | 
\author{Christopher Joseph Fennell}  | 
| 22 | 
  | 
\work{Dissertation} | 
| 23 | 
  | 
\degprior{B.Sc.} | 
| 28 | 
  | 
\maketitle | 
| 29 | 
  | 
 | 
| 30 | 
  | 
\begin{abstract} | 
| 31 | 
+ | 
 | 
| 32 | 
+ | 
This dissertation comprises a body of research in the field of | 
| 33 | 
+ | 
classical molecular simulations, with particular emphasis placed on | 
| 34 | 
+ | 
the proper depiction of water. This work is arranged such that the | 
| 35 | 
+ | 
techniques and models used within are first developed and tested | 
| 36 | 
+ | 
before being applied and compared with experimental results. With this | 
| 37 | 
+ | 
organization in mind, it is appropriate that the first chapter deals | 
| 38 | 
+ | 
primarily the technique of molecular dynamics and technical | 
| 39 | 
+ | 
considerations needed to correctly perform molecular simulations. | 
| 40 | 
+ | 
 | 
| 41 | 
+ | 
Building on this framework, the second chapter discusses correction | 
| 42 | 
+ | 
techniques for handling the long-ranged electrostatic interactions | 
| 43 | 
+ | 
common in molecular simulations. Particular focus is placed on a | 
| 44 | 
+ | 
shifted-force ({\sc sf}) modification of the damped shifted Coulombic | 
| 45 | 
+ | 
summation method. In this work, {\sc sf} is shown to be nearly | 
| 46 | 
+ | 
equivalent to the more commonly utilized Ewald summation in | 
| 47 | 
+ | 
simulations of condensed phases. Since the {\sc sf} technique is | 
| 48 | 
+ | 
pairwise, it scales as $\mathcal{O}(N)$ and lacks periodicity | 
| 49 | 
+ | 
artifacts introduced through heavy reliance on the reciprocal-space | 
| 50 | 
+ | 
portion of the Ewald sum. The electrostatic damping technique used | 
| 51 | 
+ | 
with {\sc sf} is then extended beyond simple charge-charge | 
| 52 | 
+ | 
interactions to include point-multipoles. Optimal damping parameter | 
| 53 | 
+ | 
settings are also determined to ensure proper depiction of the | 
| 54 | 
+ | 
dielectric behavior of molecular systems. Presenting this technique | 
| 55 | 
+ | 
early enables its application in the systems discussed in the later | 
| 56 | 
+ | 
chapters and shows how it can improve the quality of various molecular | 
| 57 | 
+ | 
simulations. | 
| 58 | 
+ | 
 | 
| 59 | 
+ | 
The third chapter applies the above techniques and focuses on water | 
| 60 | 
+ | 
model development, specifically the single-point soft sticky dipole | 
| 61 | 
+ | 
(SSD) model. In order to better depict water with SSD in computer | 
| 62 | 
+ | 
simulations, it needed to be reparametrized. This work results in the | 
| 63 | 
+ | 
development of SSD/RF and SSD/E, new variants of the SSD model | 
| 64 | 
+ | 
optimized for simulations with and without a reaction field | 
| 65 | 
+ | 
correction. These new single-point models are more efficient than the | 
| 66 | 
+ | 
common multi-point partial charge models and better capture the | 
| 67 | 
+ | 
dynamic properties of water. SSD/RF can be successfully used with | 
| 68 | 
+ | 
damped {\sc sf} through the multipolar extension of the technique | 
| 69 | 
+ | 
described in the previous chapter. Discussion on the development of | 
| 70 | 
+ | 
the two-point tetrahedrally restructured elongated dipole (TRED) water | 
| 71 | 
+ | 
model is also presented, and this model is optimized for use with the | 
| 72 | 
+ | 
damped {\sc sf} technique. Though there remain some algorithmic | 
| 73 | 
+ | 
complexities that need to be addressed (logic for neglecting | 
| 74 | 
+ | 
charge-quadrupole interactions between other TRED molecules) to use | 
| 75 | 
+ | 
this model in general simulations, it is approximately twice as | 
| 76 | 
+ | 
efficient as the commonly used three-point water models (i.e. TIP3P | 
| 77 | 
+ | 
and SPC/E). | 
| 78 | 
+ | 
 | 
| 79 | 
+ | 
Continuing in the direction of model applications, the final chapter | 
| 80 | 
+ | 
deals with a unique polymorph of ice that was discovered while | 
| 81 | 
+ | 
performing water simulations with the fast simple water models | 
| 82 | 
+ | 
discussed in the previous chapter. This form of ice, called | 
| 83 | 
+ | 
``imaginary ice'' (Ice-$i$), has a low-density structure which is | 
| 84 | 
+ | 
different from any known polymorph observed in either experiment or | 
| 85 | 
+ | 
computer simulation studies. The free energy analysis discussed here | 
| 86 | 
+ | 
shows that this structure is in fact the thermodynamically preferred | 
| 87 | 
+ | 
form of ice for both the single-point and commonly used multi-point | 
| 88 | 
+ | 
water models under the chosen simulation conditions. It is shown that | 
| 89 | 
+ | 
inclusion of electrostatic corrections is necessary to obtain more | 
| 90 | 
+ | 
realistic results; however, the free energies of the various | 
| 91 | 
+ | 
polymorphs (both imaginary and real) in many of these models is shown | 
| 92 | 
+ | 
to be so similar that choice of system properties, like the volume in | 
| 93 | 
+ | 
$NVT$ simulations, can directly influence the ice polymorph expressed. | 
| 94 | 
+ | 
 | 
| 95 | 
  | 
\end{abstract} | 
| 96 | 
  | 
 | 
| 97 | 
  | 
\begin{dedication} | 
| 113 | 
  | 
thank Christopher Harrison and Dr. Steven Corcelli for additional | 
| 114 | 
  | 
discussions and comments. Finally, I would like to thank my parents, | 
| 115 | 
  | 
Edward P.~Fennell and Rosalie M.~Fennell, for providing the | 
| 116 | 
< | 
opprotunities and encouragement that allowed me to pursue my | 
| 116 | 
> | 
opportunities and encouragement that allowed me to pursue my | 
| 117 | 
  | 
interests, and I would like to thank my wife, Kelley, for her | 
| 118 | 
< | 
unwaivering support. | 
| 118 | 
> | 
unwavering support. | 
| 119 | 
  | 
\end{acknowledge} | 
| 120 | 
  | 
 | 
| 121 | 
  | 
\mainmatter |