Acta Phys. -Chim. Sin. ›› 2011, Vol. 27 ›› Issue (11): 2547-2552.doi: 10.3866/PKU.WHXB20111107

• THEORETICAL AND COMPUTATIONAL CHEMISTRY • Previous Articles     Next Articles

Hydrogen Bond Lifetime Definitions and the Relaxation Mechanism in Water Solutions

ZHANG Xia1, ZHANG Qiang1, ZHAO Dong-Xia2   

  1. 1. Institute of Chemistry, Chemical Engineering and Food Safety, Bohai University, Jinzhou 121000, Liaoning Province, P. R. China;
    2. Institute of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, Liaoning Proveince, P. R. China
  • Received:2011-07-08 Revised:2011-08-22 Published:2011-10-27
  • Contact: ZHANG Qiang E-mail:zhangqiang@bhu.edu.cn
  • Supported by:

    The project was supported by the National Natural Science Foundation of China (20873055).

Abstract: The molecular dynamics behaviors in water solutions are determined by the hydrogen bond (H-bond) relaxations. The H-bond lifetime, as an important experimental and theoretical parameter, is often used to explore the general kinetics of H-bond dynamics. In this work, four different H-bond lifetimes were defined and calculated in dimethyl sulfoxide (DMSO)-water mixtures with two widely-used combined force fields, SPC/E-P2 and SPC/E-OPLS. The continuous and kinetic based H-bond lifetimes, τC and τR, are always shorter than the τPR of stable states due to neglecting of unsuccessful H-bond exchanges. The intermittent H-bond lifetime τI was found to be the longest because of a recount of the reforming events after the successful switching event. The H-bond lifetimes, τC, τI, τR, and τPR increase with the mole fraction of DMSO (xD). This trend is not consistent with that of the molecular diffuse constants. This shows that the molecular mobility is not a decisive factor to the H-bond lifetime. The environment-dependent H-bond lifetimes suggest that the stronger H-bonds should not always remain longer time. The H-bond coordination numbers of water and DMSO decrease with xD. The distortion and elongation probability of the H-bond that was induced by surrounding molecules decreases and, therefore, so the τC and τR approach each other at the limiting concentrations in this work. The facts above show that the labeled H-bond lifetime is closely related to the H-bond density around it. One H-bond switching event only takes place on one new available acceptor there. The localized character of H-bond relaxation is consistent with the trend of the molecular mobility trend. The H-bond lifetimes also rely on the theoretical model used in the simulations.

Key words: Hydrogen bond lifetime, Molecular dynamics simulation, Hydrogen bond exchange, Correlation function, Force field

MSC2000: 

  • O643