计算化学公社

标题: 高斯CASSCF频率计算 [打印本页]

作者
Author: ter20 时间: 2017-3-15 14:24
标题: 高斯CASSCF频率计算
最近在用高斯做CASSCF计算，体系是一个17原子，67电子的小分子正离子，遇到一个问题，当优化做完后，继续做频率时，会莫名其妙地停止，也没有任何错误信息，关键字如下：

%chk=2-631G-d-opt.chk

%nproc=16

%mem=3200MW

# CAS(10,9,NRoot=2)/6-31G(d) Guess=(read) Pop=Full Nosymm geom=allcheck opt freq=noraman

频率计算开始大约半小时后，就自动停掉，临时文件Gau-XXXX.rwf等都在，log文件末尾信息如下：

ITN= 24 MaxIt= 64 E= -382.3852066022 DE= 5.94D-09 Acc= 1.00D-08 Lan= 0

... Do an extra-iteration for final printing.

Final one electron symbolic 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.000000D+00

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

6 7 8 9

6 0.200000D+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

连续算了几个分子都是这样，不知道是何原因，求教

作者
Author: pwzhou 时间: 2017-3-15 14:54
对09和之前的版本而言，当活化空间大于（8,8）的时候，优化交叉点和算频率时请用NoCPMCSCF关键词，另外，算频率还得用Freq=Numer，这样就不会出现你这个问题了，要不然你的内存根本不够用。16是否还存在这个问题未知，谨慎估计还存在。

作者
Author: ter20 时间: 2017-3-15 15:19

pwzhou 发表于 2017-3-15 14:54
对09和之前的版本而言，当活化空间大于（8,8）的时候，优化交叉点和算频率时请用NoCPMCSCF关键词，另外，算 ...

感谢！

作者
Author: ter20 时间: 2017-3-15 18:21

pwzhou 发表于 2017-3-15 14:54
对09和之前的版本而言，当活化空间大于（8,8）的时候，优化交叉点和算频率时请用NoCPMCSCF关键词，另外，算 ...

发现很奇怪的结果，加了NoCPMCSCF和Freq=Numer后，Opt一步收敛，如下：
Item                               Value    Threshold  Converged?
Maximum Force          0.000000    0.000450    YES
RMS    Force             0.000000    0.000300    YES
Maximum Displacement 0.000000    0.001800    YES
RMS    Displacement    0.000000    0.001200    YES

频率计算也是45个虚频，感觉不太对

作者
Author: 函数与激情 时间: 2017-4-6 16:27

pwzhou 发表于 2017-3-15 14:54
对09和之前的版本而言，当活化空间大于（8,8）的时候，优化交叉点和算频率时请用NoCPMCSCF关键词，另外，算 ...

我之前cas优化交叉点的确（12,10）的空间算不了，师姐让我加了这个关键词之后可以跑，请问这个关键词是什么意思呢？

作者
Author: pwzhou 时间: 2017-4-7 08:28

ter20 发表于 2017-3-15 18:21
发现很奇怪的结果，加了NoCPMCSCF和Freq=Numer后，Opt一步收敛，如下：
Item ...

我应该是记错了，算频率不需要加这个NoCPMCSCF关键词，用Freq=Numer即可。我很少用Gaussian做CASSCF，那个关键词是之前组里有学生出现问题问Gaussian公司得到的答案，然后我记忆出了差错，sorry。

作者
Author: ter20 时间: 2017-4-8 21:08

pwzhou 发表于 2017-4-7 08:28
我应该是记错了，算频率不需要加这个NoCPMCSCF关键词，用Freq=Numer即可。我很少用Gaussian做CASSCF，那 ...

多谢回复！

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