N2+ 激发态势能面扫描问题

TDHFjiang

老师好，我最近在做N2+ 的激发态势能扫描，我用的是EOMCCSD/aug-cc-PVTZ , 一共算了10个态，从N-N 1a.u. 到 5.a.u，扫描关键词如下图(1), 扫描结果如下图（2）
参考的势能面图为下图（3）参考文献：PHYSICAL REVIEW A 78, 033430 (2008)  ， 师兄用molpro 扫描势能面为下图（4)
问题1. 出现莫名其妙的点，如图（2），是不是关键词写错了
问题2. 没有奇怪的点，我这个扫描的结果，和别人文献中相差很大，算这个没有经验，感谢大家的帮助
问题3. 如图（4）， 用molpro算出的结果，除了图中标注的拐点，与参考文献中不一致，其他都对的很好，请问这个如何解决呢？

TDHFjiang

图中没有拐点

ZMINER

请问 您是用高斯扫的吗？

TDHFjiang

ZMINER 发表于 2018-12-26 10:22
请问 您是用高斯扫的吗？

是的

ZMINER

TDHFjiang 发表于 2018-12-26 10:33
是的

等老师回答， 我学习学习

Mikasa

1.你参考的文献并不是能量曲线的原始出处，从图的标注可以看出这些是“simplified potential energy curve”。
2.Gaussian扫描的可能是部分点没有收敛到正确的态，而且你也没指定态的对称性，原始文献中扫描考虑了对称性。另外，你师兄molpro的结果应该考虑对称行了吧，也是EOMCCSD的么？处理这种多参考的体系Molpro要比Gaussian强很多，并且EOMCCSD精度并不很好。

TDHFjiang

Mikasa 发表于 2018-12-27 02:38
1.你参考的文献并不是能量曲线的原始出处，从图的标注可以看出这些是“simplified potential energy curve ...

感谢您的回复~

1. https://doi.org/10.1088/0953-4075/39/2/006
A. Ehresmann, L. Werner, S. Klumpp, S. Lucht, H. Schmoranzer,
S. Mickat, R. Schill, K.-H. Schartner, Ph. V. Demekhin,
M. P. Lemeshko, and V. L. Sukhorukov, J. Phys. B 39, 283
(2006)
不好意思，这个是其中的参考文献
2. 他用的是MRCI, 我用Molcas 试试看

xylz6188

是不是自旋多重度弄错了

Mikasa

TDHFjiang 发表于 2018-12-27 09:41
感谢您的回复~

1. https://doi.org/10.1088/0953-4075/39/2/006 [/backcolor ...

为啥不统一方法呢，为啥要试molcas呀。
先明确一下你想要什么方法，ms-caspt2精度应该够用，molpro很合适。

TDHFjiang

Mikasa 发表于 2018-12-28 02:27
为啥不统一方法呢，为啥要试molcas呀。
先明确一下你想要什么方法，ms-caspt2精度应该够用，molpro很合 ...

好的，感谢回答，正在尝试用Molpro

SigFig

可以试试更大的基组(aug-cc-PVQZ和aug-cc-PV5Z)来解决PES断裂.

关于CI 和CC 的优劣: CI由于Size Inconsistency 的问题在长距离上的精度肯定不如CC（即使是看起来更好也是right result for wrong reason)。
Size Inconsistency的意思是E(AB) != E(A) + E(B)

关于multi-reference (MR)：如果需要和文献作比较，需要MR-EOM-CCSD(orca里就有)，不过从你师兄的结果来看不太需要。