Density-Functional Tight-Binding Molecular Dynamics Simulations of Excess Proton Diffusion in Ice Ih, Ice Ic, Ice III, and Melted Ice VI Phases

Aditya Wibawa Sakti, Yoshifumi Nishimura, Chien Pin Chou, Hiromi Nakai*

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

9 Citations (Scopus)


The structural, dynamical, and energetic properties of the excess proton in ice were studied using density-functional tight-binding molecular dynamics simulations. The ice systems investigated herein consisted of low-density hexagonal and cubic crystalline variants (ice Ih and Ic) and high-density structures (ice III and melted ice VI). Analysis of the temperature dependence of radial distribution function and bond order parameters served to characterize the distribution and configuration of hundreds of water molecules in a unit cell. We confirmed that ice Ih and Ic possess higher hexagonal symmetries than ice III and melted ice VI. The estimated Grotthuss shuttling diffusion coefficients in ice were larger than that of liquid water, indicating a slower proton diffusion process in high-density structures than in low-density systems. The energy barriers calculated on the basis of the Arrhenius plot of diffusion coefficients were in reasonable agreement with experimental measurement for ice Ih. Furthermore, the energy barriers for high-density structures were several times larger than those of low-density systems. The simulation results were likely related to the suppression of proton transfer in disordered water configurations, in particular, ice with low hexagonal symmetry.

Original languageEnglish
Pages (from-to)33-40
Number of pages8
JournalJournal of Physical Chemistry A
Issue number1
Publication statusPublished - 2018 Jan 11

ASJC Scopus subject areas

  • Physical and Theoretical Chemistry


Dive into the research topics of 'Density-Functional Tight-Binding Molecular Dynamics Simulations of Excess Proton Diffusion in Ice I<sub>h</sub>, Ice I<sub>c</sub>, Ice III, and Melted Ice VI Phases'. Together they form a unique fingerprint.

Cite this