物理化学学报 >> 2011, Vol. 27 >> Issue (11): 2547-2552.doi: 10.3866/PKU.WHXB20111107

理论与计算化学 上一篇    下一篇

水溶液中氢键寿命的定义和弛豫机理

张霞1, 张强1, 赵东霞2   

  1. 1. 渤海大学化学化工与食品安全学院, 辽宁锦州 121000;
    2. 辽宁师范大学化学化工学院, 辽宁大连 116029
  • 收稿日期:2011-07-08 修回日期:2011-08-22 发布日期:2011-10-27
  • 通讯作者: 张强 E-mail:zhangqiang@bhu.edu.cn
  • 基金资助:

    国家自然科学基金(20873055)资助项目

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).

摘要: 氢键弛豫过程决定水溶液中分子的动力学行为, 氢键寿命作为理论和实验结果中一个重要的参数通常用以描述溶液中的氢键动力学性质. 本文采用SPC/E-P2和SPC/E-OPLS两种力场方法对二甲基亚砜(DMSO)水溶液体系进行分子动力学模拟, 计算了四种不同定义下的氢键寿命, 并进行了对比阐述. 连续性氢键寿命τC和基于动力学平衡方法下的氢键寿命τR较短, 它们忽略了不成功的氢键交换过程, 稳定态模型下的氢键寿命τPR能够真实体现氢键转换过程, 可作为标度衡量其它量值. 间歇性氢键寿命τI相对最长, 重复统计了氢键成功交换后的恢复概率. 随着二甲基亚砜浓度增大, τCτIτRτPR不断增大, 分子平动扩散系数在中等浓度达到极小, 说明氢键寿命与分子活动性无决定性关联, 同类型氢键在不同浓度下的寿命差异表明氢键寿命具有分子环境依赖性; 水分子和二甲基亚砜分子氢键个数降低, 氢键受激角度扭曲和拉伸概率下降, τCτR不断接近, 氢键交换受体密度降低使氢键交换速率下降, τPR不断增大, 氢键寿命与其周围氢键密度密切相关. τIτPR的比值体现了氢键交换的局域性, 其变化趋势与分子平动活动性趋势相同. 两种力场下氢键寿命的差异也说明氢键的寿命具有明显的理论模型依赖性.

关键词: 氢键寿命, 分子动力学模拟, 氢键交换, 相关函数, 分子力场

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