LJ势能曲线的绘制（二）

diaolanxinyu

接着上一个帖子：为下一步大体系的计算我测试了一个小体系：是氧负离子自由基与甲烷，先对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]


谢谢大家！

sobereva

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

yflchx

带负电荷，静电作用那么强，用LJ势能表达式合适吗？从物理上来说。

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