计算化学公社

标题: LJ势能曲线的绘制（二） [打印本页]

作者
Author: diaolanxinyu 时间: 2015-2-3 23:02
标题: LJ势能曲线的绘制（二）
接着上一个帖子：为下一步大体系的计算我测试了一个小体系：是氧负离子自由基与甲烷，先对O和甲烷中H距离的变化进行了scan扫描，（#p scan mp2/aug-cc-pvtz ）
输出结果中有四列数据，如下：
Summary of the potential surface scan:
N       R1          SCF             MP2
---- --------- -----------       -----------
1    1.0       -114.92160 -115.38199
2    1.1       -114.94759 -115.39834
3    1.2       -114.96646 -115.40957
4    1.3       -114.98024 -115.41777
将2，4列用origin作图（我取最后一步O和H的距离差不多41埃的单点能作为参考点），如图1

问题一：纵坐标应该写势能还是结合能或者其他？

问题二：用上述结果如何得到LJ势能参数中的epsilon和sigma呢？
LJ势能表达式：U(r)=4*epsilon *[(sigma/r)^12-(sigma/r)^6]

谢谢大家！

作者
Author: sobereva 时间: 2015-2-4 04:04
1 如果直接用体系的MP2的值，就写能量。如果是体系能量减去O和甲烷的能量，就写相互作用能。
2 用origin之类的程序拟合即得到sigma和epsilon，很容易

作者
Author: yflchx 时间: 2015-2-4 09:42
带负电荷，静电作用那么强，用LJ势能表达式合适吗？从物理上来说。

当然，拟合是没啥问题的。

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