计算化学公社

标题: 请教如何解决文献中提供的L-J势排斥项吸引项系数不一致问题 [打印本页]

作者
Author: 芝士栗子小甜饼 时间: 2019-11-23 11:37
标题: 请教如何解决文献中提供的L-J势排斥项吸引项系数不一致问题
最近想要构建nafion立场，参考的文献中L-J势的形式如下图
因为吸引项排斥项系数不一致而不知道怎样处理。
文献中用的组合规则是Lorentz-Berthelot规则

[img]file:///C:\Users\Administrator\AppData\Roaming\Tencent\Users\735863122\QQ\WinTemp\RichOle\4[F~%P`RNRL(4AU8QBY`(PJ.png[/img]

作者
Author: 芝士栗子小甜饼 时间: 2019-11-23 11:51
不知道这样的解决办法是否合理：

作者
Author: Lacrimosa 时间: 2019-11-23 12:07
好像没这么复杂吧，把R0除以2^(1/6)再带入进去就好了

作者
Author: 芝士栗子小甜饼 时间: 2019-11-23 14:57

Lacrimosa 发表于 2019-11-23 12:07
好像没这么复杂吧，把R0除以2^(1/6)再带入进去就好了

那排斥项的系数不也改变了吗

作者
Author: Lacrimosa 时间: 2019-11-23 15:23

芝士栗子小甜饼发表于 2019-11-23 14:57
那排斥项的系数不也改变了吗

你去gmx手册里看一下lj函数的形式对比一下就知道了，直接可以提出一个常数4

作者
Author: 芝士栗子小甜饼 时间: 2019-11-23 15:35

Lacrimosa 发表于 2019-11-23 15:23
你去gmx手册里看一下lj函数的形式对比一下就知道了，直接可以提出一个常数4

因为手册上LJ的表达式中吸引项和排斥项的系数是一致的，而我困扰的就是无法将文献上的这个式子变成一致的形式

作者
Author: Lacrimosa 时间: 2019-11-23 16:37

芝士栗子小甜饼发表于 2019-11-23 15:35
因为手册上LJ的表达式中吸引项和排斥项的系数是一致的，而我困扰的就是无法将文献上的这个式子变成一致的 ...

看图

作者
Author: 芝士栗子小甜饼 时间: 2019-11-23 17:05

Lacrimosa 发表于 2019-11-23 16:37
看图

我懂了！谢谢你

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