化学学报 ›› 2005, Vol. 63 ›› Issue (5): 358-362. 上一篇    下一篇

研究论文

电解质水溶液传递性质的布朗动力学模拟研究

史红兵,于养信,高光华*   

  1. (清华大学化学工程系 北京 100084)
  • 投稿日期:2004-08-05 修回日期:2004-11-08 发布日期:2010-12-10
  • 通讯作者: 高光华

Study on the Transport Properties of Aqueous Electrolyte Solution by Brownian Dynamics Simulation

SHI Hong-Bing, YU Yang-Xin, GAO Guang-Hua*   

  1. (Department of Chemical Engineering, Tsinghua University, Beijing 100084)
  • Received:2004-08-05 Revised:2004-11-08 Published:2010-12-10
  • Contact: GAO Guang-Hua

在传统布朗动力学的基础上, 考虑流体力学的影响, 并且引入Smart Monte Carlo方法的接受概率, 对电解质溶液进行布朗动力学模拟, 得到不同浓度和温度下KCl溶液中离子间的径向分布函数, 并且与超网链积分方程理论计算结果进行了比较, 同时, 模拟了KCl和NaCl溶液的摩尔电导率. 模拟过程基于电解质溶液的原始模型, 溶剂被看作连续介质, 溶质分子之间的相互作用采用软核加静电的势能函数模型, 长程静电力的处理采用Ewald加和方法. 结果显示, 流体力学的作用对于电解质溶液的结构性质没有明显的影响, 但是对于传递性质的影响显著; 考虑流体力学作用的布朗动力学模拟结果与实验数据吻合良好.

关键词: 布朗动力学, 电解质溶液, 径向分布函数, 摩尔电导率

Brownian dynamics simulation was carried out to study the electrolyte solution. The effect of hydrodynamic interactions was taken into account and an acceptance criterion known from the Smart Monte Carlo method was combined with the traditional method. The radial distribution functions for both like and unlike-ion interactions were calculated at different concentration and temperature and compared with those from HNC integral equation. The molar conductivity of sodium chloride and potassium chloride aqueous solution was also obtained and compared with the experiment values. In our simulation, the solvent-averaged interaction potential between the ions was modeled by pairwise repulsive soft-core interactions and Coulomb forces that were handled by Ewald summation technique. It is found that the hydrodynamic interactions have no effect on the structure properties in contrast with an obvious effect on transporting properties. Good agreement between simulation results and the corresponding values from experiment was achieved by including hydrodynamic interactions.

Key words: brownian dynamics simulation, electrolyte solution, radial distribution function, molar conductivity