计算化学公社

标题: LJ势中σ、ε的含义以及关于原子半径的疑问 [打印本页]

作者
Author: jimulation 时间: 2020-12-22 23:29
标题: LJ势中σ、ε的含义以及关于原子半径的疑问
阅读了社长写的区分各种原子半径的博文http://sobereva.com/255，结合GMX的使用过程，对LJ势函数有两个疑问：
(, 下载次数 Times of downloads: 104)
(, 下载次数 Times of downloads: 104)

第一，σ是LJ势为0时的原子间距离，这一距离与博文中讲的原子半径有何关系？第二，ε(势阱深度)是原子间距离r=rmin=2^(1/6)σ时的LJ势，ε看似与原子所处的环境密切相关？

举例来说，同样是Na+，NaCl晶体中与水溶液中的σ值有无显著不同？晶体中的ε值是否明显强于溶液中(因晶体中相互作用更强)？此外，对于未知LJ参数的某种离子晶体(比如Cu2O)，有没有方法粗略估计两种原子的σ和ε值？

作者
Author: sobereva 时间: 2020-12-23 08:17
1 这得看什么原子半径。大多数原子半径的定义和sigma并无关系
UFF的范德华半径对应的是UFF力场的同种元素间LJ势的r_min的一半

2 有的力场对于同种元素的不同原子类型用不同的范德华参数，这和实际环境中原子带的电子数有密切关系。

没有简单的、普适的粗略估计方法。

作者
Author: jimulation 时间: 2020-12-23 09:57

sobereva 发表于 2020-12-23 08:17
1 这得看什么原子半径。大多数原子半径的定义和sigma并无关系
UFF的范德华半径对应的是UFF力场的同种元素 ...

谢谢社长的回复，大致明白了，我再理解下

作者
Author: ocbrother 时间: 2020-12-23 10:08
请问你这里的两个计算VDW相互作用的公式是等价的吗。
为什么第二个公式变成了-2

作者
Author: jimulation 时间: 2020-12-23 10:25

ocbrother 发表于 2020-12-23 10:08
请问你这里的两个计算VDW相互作用的公式是等价的吗。
为什么第二个公式变成了-2

是的，可以等价转化的

作者
Author: ocbrother 时间: 2020-12-23 10:28

jimulation 发表于 2020-12-23 10:25
是的，可以等价转化的

我推了一下确实是等价的，可是书上写的是没有2的。。
应该是书上错了。。。

作者
Author: sobereva 时间: 2020-12-24 12:24

ocbrother 发表于 2020-12-23 10:28
我推了一下确实是等价的，可是书上写的是没有2的。。
应该是书上错了。。。

有没有2对应不同情况
有2时分子对应势阱极小点距离，没2时对应作用势为0的距离

作者
Author: ocbrother 时间: 2020-12-24 15:05

sobereva 发表于 2020-12-24 12:24
有没有2对应不同情况
有2时分子对应势阱极小点距离，没2时对应作用势为0的距离

谢谢老师。已经确认过是书上错了

作者
Author: xiaoruoqiu 时间: 2022-8-30 15:29
求证一个基础问题，写在itp文件里的，是公式（1）中的σ，还是公式（2）中的rmin?

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