计算化学公社

标题: 请教表面张力的单位换算 [打印本页]

作者
Author: 霜晨月 时间: 2018-10-13 20:34
标题: 请教表面张力的单位换算
用g_mmpbsa计算结合自由能时，SASA model参数里面有一项是
;Gamma (Surface Tension) kJ/(mol A^2)
gamma = 0.0226778
这个应该就是水的表面张力。但我查得水的表面张力应为72.80 mN/m （或72.80 mJ/m^2），尝试着换算，怎么也换算不出 0.0226778 kJ/(mol A^2)来，总是差了很多倍。
请问mN/m该如何换算成 kJ/(mol A^2)？关键是分母上的那个mol量纲如何消掉？

（问这个问题的目的是，我的溶剂是环己烷，24.95mN/m，该如何换算成kJ/(mol A^2)？）

作者
Author: Daniel_Arndt 时间: 2018-10-13 23:54
Avogadro's constant: NA=6.02214076×10^23 mol^-1

作者
Author: 霜晨月 时间: 2018-10-14 09:43

Daniel_Arndt 发表于 2018-10-13 23:54
Avogadro's constant: NA=6.02214076×10^23 mol^-1

谢谢，这个我试过了。
72.80 mN/m = 72.80 mJ/m^2 = 72.80e-3 J/m^2 = 72.80e-23 J/A^2
乘以6.02214076e+23/mol = 438.4 J/(mol A^2) = 0.4384 KJ/(mol A^2)，与上面的数值差了19.33倍。
不知道问题出在哪里？

作者
Author: hdmd 时间: 2018-10-15 11:40
我一直对目前MD常用基于压力计算表面张力的方法不认同（注意，是液体的表面张力不是液-液的界面张力）

举个例子：我分别用1nm的水膜和2nm的水膜计算表面张力，数值是不是应该一样？如果认为一样，那么我告诉你一个事实：如果在npt系综下（或者是x,y,z分别控压），水分子足够多的情况下，水膜会由1nm变为2nm最终充满整个空间，换句话来讲，在定温定压条件下2nm的水膜比1nm的水膜更稳定！MD计算表面张力所得数值是与计算体系直接相关的，这明显与事实相反！

非常欢迎反驳，我也对人们为什么采用这种方法计算感到好奇，也许有我没考虑到的因素。

ps：个人认为，计算表面张力以水滴模型计算才比较靠谱，但是本人不主要做这方面工作所以没有深入研究

作者
Author: diaok 时间: 2018-10-17 14:36

hdmd 发表于 2018-10-15 11:40
我一直对目前MD常用基于压力计算表面张力的方法不认同（注意，是液体的表面张力不是液-液的界面张力）

...

水膜表面张力就不该npt吧？
平衡后表面还会存在？

作者
Author: hdmd 时间: 2018-10-17 15:42

diaok 发表于 2018-10-17 14:36
水膜表面张力就不该npt吧？
平衡后表面还会存在？

计算是用nvt系综算

我的意思是1nm的膜和2nm的膜计算结果不一样

npt那个只是为了分析

欢迎光临计算化学公社 (http://bbs.keinsci.com/)