Quantum nuclear ab Initio molecular dynamics study of water wires

Hsiao S. Mei, Mark Tuckerman, Diane E. Sagnella, Michael L. Klein

Research output: Contribution to journalArticle

Abstract

The structure of linear water wires with an excess proton was studied at room temperature using ab initio path integral molecular dynamics. The ab initio Car - Parrinello (CP) methodology employed the density functional theory (DFT) description of the electronic structure, and the Feynman path integral approach allowed for quantization of the nuclear degrees of freedom. Thus, the influence of proton tunneling and zero point nuclear vibrations were automatically included. Four or five water molecules were linearly arranged, with an excess proton (H* ). to form tetramer and pentamer complexes, respectively. In classical studies of the tetramer complex, the excess proton H*, centered within the wire, formed H3O+ and H5O2+ ions with the two inner water molecules. In the pentamer complex, the H* was found attached to the inner water molecule, forming a stable H3O+ ion with two covalent, hyperextended bonds that were hydrogen bonded to neighboring water molecules on both opposite sides. Although the addition of nuclear quantization via path integrals broadened the calculated distribution functions for both complexes, the overall features were unaltered, which suggests that nuclear quantum effects are minimal in these small, linear clusters. However, instantaneous path integral configurations revealed the formation of an extended H7O3+ complex predominantly in the pentamer wire, where the excess proton H* was delocalized over three adjacent water molecules simultaneously. Since the computational demands of CP make long simulations cost-prohibitive, angular distribution functions, requiring much longer simulation times, were obtained using an MP2-based empirical valence bond (EVB) model [Sagella, D. E.; Tuckerman, M. E. J. Chem. Phys. 1998, 108, 2073]. Additional classical CP calculations. where the water wire ends were solvated with additional capping waters, were also performed. In these studies, the proton was observed to be much more mobile; proton transfer occurred along the full water wire and occasionally into the waler solvation caps.

Original languageEnglish (US)
Pages (from-to)10446-10458
Number of pages13
JournalJournal of Physical Chemistry B
Volume102
Issue number50
StatePublished - Dec 10 1998

Fingerprint

Molecular dynamics
wire
Wire
molecular dynamics
Protons
Water
water
protons
Hydrogen
Molecules
Railroad cars
molecules
Distribution functions
distribution functions
Ions
Proton transfer
Covalent bonds
Angular distribution
Solvation
covalent bonds

ASJC Scopus subject areas

  • Physical and Theoretical Chemistry

Cite this

Mei, H. S., Tuckerman, M., Sagnella, D. E., & Klein, M. L. (1998). Quantum nuclear ab Initio molecular dynamics study of water wires. Journal of Physical Chemistry B, 102(50), 10446-10458.

Quantum nuclear ab Initio molecular dynamics study of water wires. / Mei, Hsiao S.; Tuckerman, Mark; Sagnella, Diane E.; Klein, Michael L.

In: Journal of Physical Chemistry B, Vol. 102, No. 50, 10.12.1998, p. 10446-10458.

Research output: Contribution to journalArticle

Mei, HS, Tuckerman, M, Sagnella, DE & Klein, ML 1998, 'Quantum nuclear ab Initio molecular dynamics study of water wires', Journal of Physical Chemistry B, vol. 102, no. 50, pp. 10446-10458.
Mei, Hsiao S. ; Tuckerman, Mark ; Sagnella, Diane E. ; Klein, Michael L. / Quantum nuclear ab Initio molecular dynamics study of water wires. In: Journal of Physical Chemistry B. 1998 ; Vol. 102, No. 50. pp. 10446-10458.
@article{ad7e7722b8a1418ea6e48b0d9efd2225,
title = "Quantum nuclear ab Initio molecular dynamics study of water wires",
abstract = "The structure of linear water wires with an excess proton was studied at room temperature using ab initio path integral molecular dynamics. The ab initio Car - Parrinello (CP) methodology employed the density functional theory (DFT) description of the electronic structure, and the Feynman path integral approach allowed for quantization of the nuclear degrees of freedom. Thus, the influence of proton tunneling and zero point nuclear vibrations were automatically included. Four or five water molecules were linearly arranged, with an excess proton (H* ). to form tetramer and pentamer complexes, respectively. In classical studies of the tetramer complex, the excess proton H*, centered within the wire, formed H3O+ and H5O2+ ions with the two inner water molecules. In the pentamer complex, the H* was found attached to the inner water molecule, forming a stable H3O+ ion with two covalent, hyperextended bonds that were hydrogen bonded to neighboring water molecules on both opposite sides. Although the addition of nuclear quantization via path integrals broadened the calculated distribution functions for both complexes, the overall features were unaltered, which suggests that nuclear quantum effects are minimal in these small, linear clusters. However, instantaneous path integral configurations revealed the formation of an extended H7O3+ complex predominantly in the pentamer wire, where the excess proton H* was delocalized over three adjacent water molecules simultaneously. Since the computational demands of CP make long simulations cost-prohibitive, angular distribution functions, requiring much longer simulation times, were obtained using an MP2-based empirical valence bond (EVB) model [Sagella, D. E.; Tuckerman, M. E. J. Chem. Phys. 1998, 108, 2073]. Additional classical CP calculations. where the water wire ends were solvated with additional capping waters, were also performed. In these studies, the proton was observed to be much more mobile; proton transfer occurred along the full water wire and occasionally into the waler solvation caps.",
author = "Mei, {Hsiao S.} and Mark Tuckerman and Sagnella, {Diane E.} and Klein, {Michael L.}",
year = "1998",
month = "12",
day = "10",
language = "English (US)",
volume = "102",
pages = "10446--10458",
journal = "Journal of Physical Chemistry B Materials",
issn = "1520-6106",
publisher = "American Chemical Society",
number = "50",

}

TY - JOUR

T1 - Quantum nuclear ab Initio molecular dynamics study of water wires

AU - Mei, Hsiao S.

AU - Tuckerman, Mark

AU - Sagnella, Diane E.

AU - Klein, Michael L.

PY - 1998/12/10

Y1 - 1998/12/10

N2 - The structure of linear water wires with an excess proton was studied at room temperature using ab initio path integral molecular dynamics. The ab initio Car - Parrinello (CP) methodology employed the density functional theory (DFT) description of the electronic structure, and the Feynman path integral approach allowed for quantization of the nuclear degrees of freedom. Thus, the influence of proton tunneling and zero point nuclear vibrations were automatically included. Four or five water molecules were linearly arranged, with an excess proton (H* ). to form tetramer and pentamer complexes, respectively. In classical studies of the tetramer complex, the excess proton H*, centered within the wire, formed H3O+ and H5O2+ ions with the two inner water molecules. In the pentamer complex, the H* was found attached to the inner water molecule, forming a stable H3O+ ion with two covalent, hyperextended bonds that were hydrogen bonded to neighboring water molecules on both opposite sides. Although the addition of nuclear quantization via path integrals broadened the calculated distribution functions for both complexes, the overall features were unaltered, which suggests that nuclear quantum effects are minimal in these small, linear clusters. However, instantaneous path integral configurations revealed the formation of an extended H7O3+ complex predominantly in the pentamer wire, where the excess proton H* was delocalized over three adjacent water molecules simultaneously. Since the computational demands of CP make long simulations cost-prohibitive, angular distribution functions, requiring much longer simulation times, were obtained using an MP2-based empirical valence bond (EVB) model [Sagella, D. E.; Tuckerman, M. E. J. Chem. Phys. 1998, 108, 2073]. Additional classical CP calculations. where the water wire ends were solvated with additional capping waters, were also performed. In these studies, the proton was observed to be much more mobile; proton transfer occurred along the full water wire and occasionally into the waler solvation caps.

AB - The structure of linear water wires with an excess proton was studied at room temperature using ab initio path integral molecular dynamics. The ab initio Car - Parrinello (CP) methodology employed the density functional theory (DFT) description of the electronic structure, and the Feynman path integral approach allowed for quantization of the nuclear degrees of freedom. Thus, the influence of proton tunneling and zero point nuclear vibrations were automatically included. Four or five water molecules were linearly arranged, with an excess proton (H* ). to form tetramer and pentamer complexes, respectively. In classical studies of the tetramer complex, the excess proton H*, centered within the wire, formed H3O+ and H5O2+ ions with the two inner water molecules. In the pentamer complex, the H* was found attached to the inner water molecule, forming a stable H3O+ ion with two covalent, hyperextended bonds that were hydrogen bonded to neighboring water molecules on both opposite sides. Although the addition of nuclear quantization via path integrals broadened the calculated distribution functions for both complexes, the overall features were unaltered, which suggests that nuclear quantum effects are minimal in these small, linear clusters. However, instantaneous path integral configurations revealed the formation of an extended H7O3+ complex predominantly in the pentamer wire, where the excess proton H* was delocalized over three adjacent water molecules simultaneously. Since the computational demands of CP make long simulations cost-prohibitive, angular distribution functions, requiring much longer simulation times, were obtained using an MP2-based empirical valence bond (EVB) model [Sagella, D. E.; Tuckerman, M. E. J. Chem. Phys. 1998, 108, 2073]. Additional classical CP calculations. where the water wire ends were solvated with additional capping waters, were also performed. In these studies, the proton was observed to be much more mobile; proton transfer occurred along the full water wire and occasionally into the waler solvation caps.

UR - http://www.scopus.com/inward/record.url?scp=0001557630&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0001557630&partnerID=8YFLogxK

M3 - Article

VL - 102

SP - 10446

EP - 10458

JO - Journal of Physical Chemistry B Materials

JF - Journal of Physical Chemistry B Materials

SN - 1520-6106

IS - 50

ER -