Sensitivity Analysis of Distribution Functions of Liquid Water

Sheng‐Bai Zhu, Chung F. Wong

Research output: Contribution to journalArticlepeer-review

Abstract

<div class="line" id="line-29"> <span style="font-family: Lora, serif; font-size: 20px;"> We have carried out sensitivity&nbsp;analysis&nbsp;studies to identify the most important potential parameters of the SPC and the TIP3P flexible&nbsp;water&nbsp;models in determining three&nbsp;distribution functions&nbsp;of liquid water&mdash;the site&ndash;site radial&nbsp;distribution functions,&nbsp;the&nbsp;distribution function&nbsp;of the&nbsp;interaction&nbsp;energy&nbsp;of a&nbsp;water&nbsp;molecule with its surrounding [ </span> <i style="font-family: Lora, serif; font-size: 20px;"> P </i> <i> <span style="font-family: Lora, serif; font-size: 12px;"> u </span> </i> <span style="font-family: Lora, serif; font-size: 20px;"> ( </span> <i style="font-family: Lora, serif; font-size: 20px;"> u </i> <i> <span style="font-family: Lora, serif; font-size: 12px;"> b </span> </i> <span style="font-family: Lora, serif; font-size: 20px;"> )], and the&nbsp;distribution function&nbsp;of the local&nbsp;electric field&nbsp;at the oxygen of a&nbsp;water&nbsp;molecule projected along the permanent&nbsp;dipole moment&nbsp;vector of the&nbsp;water&nbsp;molecule [ </span> <i style="font-family: Lora, serif; font-size: 20px;"> P </i> <i> <span style="font-family: Lora, serif; font-size: 12px;"> e </span> </i> <span style="font-family: Lora, serif; font-size: 20px;"> ( </span> <b style="font-family: Lora, serif; font-size: 20px;"> E </b> <span style="font-family: Lora, serif; font-size: 20px;"> &sdot;&mu; </span> <span style="font-family: Lora, serif; font-size: 12px;"> 0 </span> <span style="font-family: Lora, serif; font-size: 20px;"> )]. The site&ndash;site radial&nbsp;distribution functions&nbsp;of each&nbsp;water&nbsp;model are most sensitive to the equilibrium O&ndash;H&nbsp;bond&nbsp;length ( </span> <i style="font-family: Lora, serif; font-size: 20px;"> r </i> <span style="font-family: Lora, serif; font-size: 12px;"> OH </span> <span style="font-family: Lora, serif; font-size: 20px;"> ) and the Lennard&hyphen;Jones radius (&sigma;) of the model. In addition to these two parameters, the cut&hyphen;off radius ( </span> <i style="font-family: Lora, serif; font-size: 20px;"> R </i> <span style="font-family: Lora, serif; font-size: 12px;"> cut </span> <span style="font-family: Lora, serif; font-size: 20px;"> ) in the&nbsp;reaction&nbsp;field geometry and the atomic partial charges are also important in affecting&nbsp; </span> <i style="font-family: Lora, serif; font-size: 20px;"> P </i> <i> <span style="font-family: Lora, serif; font-size: 12px;"> u </span> </i> <span style="font-family: Lora, serif; font-size: 20px;"> ( </span> <i style="font-family: Lora, serif; font-size: 20px;"> u </i> <i> <span style="font-family: Lora, serif; font-size: 12px;"> b </span> </i> <span style="font-family: Lora, serif; font-size: 20px;"> ) of each model, especially the wings of the distribution. The oxygen charge affects the low&nbsp;energy&nbsp;wing of&nbsp; </span> <i style="font-family: Lora, serif; font-size: 20px;"> P </i> <i> <span style="font-family: Lora, serif; font-size: 12px;"> u </span> </i> <span style="font-family: Lora, serif; font-size: 20px;"> ( </span> <i style="font-family: Lora, serif; font-size: 20px;"> u </i> <i> <span style="font-family: Lora, serif; font-size: 12px;"> b </span> </i> <span style="font-family: Lora, serif; font-size: 20px;"> ) more than the high&nbsp;energy&nbsp;wing, whereas the hydrogen charge mainly affects the high&nbsp;energy&nbsp;wing and has little effect on the low&nbsp;energy&nbsp;wing of the distribution. As for&nbsp; </span> <i style="font-family: Lora, serif; font-size: 20px;"> P </i> <i> <span style="font-family: Lora, serif; font-size: 12px;"> e </span> </i> <span style="font-family: Lora, serif; font-size: 20px;"> ( </span> <b style="font-family: Lora, serif; font-size: 20px;"> E </b> <span style="font-family: Lora, serif; font-size: 20px;"> &sdot;&mu; </span> <span style="font-family: Lora, serif; font-size: 12px;"> 0 </span> <span style="font-family: Lora, serif; font-size: 20px;"> ), which provides an indirect check of the validity of the effective charge approximation in accounting for molecular&nbsp;polarizability,&nbsp;the leading factors are&nbsp; </span> <i style="font-family: Lora, serif; font-size: 20px;"> r </i> <span style="font-family: Lora, serif; font-size: 12px;"> OH </span> <span style="font-family: Lora, serif; font-size: 20px;"> , &sigma;, and&nbsp; </span> <i style="font-family: Lora, serif; font-size: 20px;"> R </i> <span style="font-family: Lora, serif; font-size: 12px;"> cut </span> <span style="font-family: Lora, serif; font-size: 20px;"> . These parameters affect the distribution wings of&nbsp; </span> <i style="font-family: Lora, serif; font-size: 20px;"> P </i> <i> <span style="font-family: Lora, serif; font-size: 12px;"> e </span> </i> <span style="font-family: Lora, serif; font-size: 20px;"> ( </span> <b style="font-family: Lora, serif; font-size: 20px;"> E </b> <span style="font-family: Lora, serif; font-size: 20px;"> &sdot;&mu; </span> <span style="font-family: Lora, serif; font-size: 12px;"> 0 </span> <span style="font-family: Lora, serif; font-size: 20px;"> ) the most. The key advantage of the sensitivity&nbsp;analysis&nbsp;technique is that it provides a systematic and economical means for studying the role of each parameter of a&nbsp;water&nbsp;model in affecting the properties of the model. &thinsp; </span></div>
Original languageAmerican English
JournalJournal of Chemical Physics
Volume99
DOIs
StatePublished - Jan 12 1993
Externally publishedYes

Disciplines

  • Atomic, Molecular and Optical Physics
  • Analytical Chemistry

Cite this