计算conical intersection时的报错该如何解决

Abandon-fmt

仿照Exploring-chemistry-with-electric-structure-methods 第二版第九章的Advanced Exercise 9.6计算conical intersection
在计算最后一步的时候遇到了很奇怪的报错
输入文件如下：
%mem=20GB
%nprocshared=32
%chk=4.chk
%oldchk=3.chk
#p CAS(10,10,NRoot=2)/def2TZVP Opt=conical Nosymm iop(1/8=5) guess=read geom=check

4

0 1

0.5 0.5

log文件最后20行信息如下：
  Final State Averaged  Density Matrix
                1             2             3             4             5
      1  0.200000D+01
      2  0.000000D+00  0.200000D+01
      3  0.000000D+00  0.000000D+00  0.200000D+01
      4  0.000000D+00  0.000000D+00  0.000000D+00  0.200000D+01
      5  0.000000D+00  0.000000D+00  0.000000D+00  0.000000D+00  0.100000D+01
      6  0.000000D+00  0.000000D+00  0.000000D+00  0.000000D+00  0.000000D+00
      7  0.000000D+00  0.000000D+00  0.000000D+00  0.000000D+00  0.000000D+00
      8  0.000000D+00  0.000000D+00  0.000000D+00  0.000000D+00  0.000000D+00
      9  0.000000D+00  0.000000D+00  0.000000D+00  0.000000D+00  0.000000D+00
     10  0.000000D+00  0.000000D+00  0.000000D+00  0.000000D+00  0.000000D+00
                6             7             8             9            10
      6  0.100000D+01
      7  0.000000D+00  0.000000D+00
      8  0.000000D+00  0.000000D+00  0.000000D+00
      9  0.000000D+00  0.000000D+00  0.000000D+00  0.000000D+00
     10  0.000000D+00  0.000000D+00  0.000000D+00  0.000000D+00  0.000000D+00
Erroneous write. Write -1 instead of 1906968.
fd = 4
orig len = 1907160 left = 1906968
g_write



而且算得异常缓慢，59k算了六个小时
请问各位大神，这个错误是什么情况？？

Abandon-fmt

pwzhou 发表于 2019-7-26 16:25
单点，不会用STO-3G去优化构型。从小基组到大基组的唯一好处就是更容易找到你需要的活化轨道，如果用大基 ...

谢谢老师，明白了，我去尝试一下。

pwzhou

Abandon-fmt 发表于 2019-7-26 12:37
不好意思，多次打扰老师了，我还有两个问题，老师有空可以麻烦帮忙解答一下么？
1）从最初的结构优化就 ...

单点，不会用STO-3G去优化构型。从小基组到大基组的唯一好处就是更容易找到你需要的活化轨道，如果用大基组找活化轨道没问题的话，不需要走这个流程。为什么STO-3G更利于选择活化轨道，因为轨道少，在非占据轨道没有一些里德堡轨道干扰你，你测试一下就知道了。你想想在10个轨道里面选5个容易还是在100个轨道里面选5个容易？

Abandon-fmt

Daniel_Arndt 发表于 2019-7-26 08:47
本论坛里面应该有关于活化空间选取的文献，我记得不错的话。

谢谢老师，我找到了那些文献

Abandon-fmt

pwzhou 发表于 2019-7-26 10:18
不好说，至少大基组的结果明显不合理。从小基组到大基组我见的比较多的是STO-3g-->4-31G-->6-31G这样的， ...

不好意思，多次打扰老师了，我还有两个问题，老师有空可以麻烦帮忙解答一下么？
1）从最初的结构优化就用STO-3G，然后用这个结果进行活性空间的选择？
2）请问“STO-3G基函数少”为什么“非常利于选择活化轨道”？

pwzhou

Abandon-fmt 发表于 2019-7-26 07:50
小基组是def2SVP，大基组def2TZVP，跨度应该不算大吧？

不好说，至少大基组的结果明显不合理。从小基组到大基组我见的比较多的是STO-3g-->4-31G-->6-31G这样的，svp到tzvp似乎挺大的，而且这样做感觉也没太大意思，STO-3G基函数少，非常利于选择活化轨道，def2-svp已经起不到这个作用了，所以这种情况下def2-tzvp就直接从头开始选轨道吧，不要用svp的结果。

Daniel_Arndt

本论坛里面应该有关于活化空间选取的文献，我记得不错的话。

Abandon-fmt

pwzhou 发表于 2019-7-25 20:31
能量没有，是不是你的基组跨度太大

小基组是def2SVP，大基组def2TZVP，跨度应该不算大吧？

pwzhou

Abandon-fmt 发表于 2019-7-25 14:31
老师，再打扰一下，请问这个totoal的能量指的是活性空间的能量还是什么能量？这个能量有合理的范围么？
...

能量没有，是不是你的基组跨度太大

pwzhou

Abandon-fmt 发表于 2019-7-25 13:25
谢谢老师，请问这本书是Exploring chemistry with electric structure methods，第二版么？

第三版

Abandon-fmt

pwzhou 发表于 2019-7-25 09:39
看上面图片的描述，第一副最下面和第二页最上面。

老师，再打扰一下，请问这个totoal的能量指的是活性空间的能量还是什么能量？这个能量有合理的范围么？
我是拿小基组的结果作为大基组的初猜，小基组算出来的占据数都在0.02-1.98之间，大基组就不是。大基组的能量是-999.502917，小基组的能量是-998.011960，这个数值合理么？

Abandon-fmt

pwzhou 发表于 2019-7-25 09:39
看上面图片的描述，第一副最下面和第二页最上面。

谢谢老师，请问这本书是Exploring chemistry with electric structure methods，第二版么？

pwzhou

Abandon-fmt 发表于 2019-7-24 08:34
老师不好意思，没大看懂“检查活化空间是否合理的一个标准就是占据轨道占据数不能接近2，非占据轨道占据 ...

看上面图片的描述，第一副最下面和第二页最上面。

Abandon-fmt

pwzhou 发表于 2019-7-22 21:34
你的活化空间肯定有问题，除了HOMO和LUMO的占据数是1，其他占据轨道的占据数是2，非占据轨道的占据数是0 ...

老师不好意思，没大看懂“检查活化空间是否合理的一个标准就是占据轨道占据数不能接近2，非占据轨道占据数不能接近0”这句话，为什么活化空间占据轨道占据数不能接近2，非占据轨道占据数不能接近0？
还有个问题，怎么判断占据轨道占据数和非占据轨道占据数？

再次感谢老师的指导！

pwzhou

Abandon-fmt 发表于 2019-7-22 13:43
请问老师，这一步我跑了三天，结果只有1M，log文件的信息是这样的：
  Final State Averaged  Density Mat ...

你的活化空间肯定有问题，除了HOMO和LUMO的占据数是1，其他占据轨道的占据数是2，非占据轨道的占据数是0，说明这些轨道都不应该包含在活化空间里面，检查活化空间是否合理的一个标准就是占据轨道占据数不能接近2，非占据轨道占据数不能接近0