计算化学公社

 找回密码 Forget password
 注册 Register
Views: 11684|回复 Reply: 6
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] 高斯CASSCF频率计算

[复制链接 Copy URL]

177

帖子

0

威望

1567

eV
积分
1744

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
楼主
最近在用高斯做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

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

368

帖子

0

威望

4691

eV
积分
5059

Level 6 (一方通行)

2#
发表于 Post on 2017-3-15 14:54:43 | 只看该作者 Only view this author
对09和之前的版本而言,当活化空间大于(8,8)的时候,优化交叉点和算频率时请用NoCPMCSCF关键词,另外,算频率还得用Freq=Numer,这样就不会出现你这个问题了,要不然你的内存根本不够用。16是否还存在这个问题未知,谨慎估计还存在。

评分 Rate

参与人数
Participants 2
eV +9 收起 理由
Reason
函数与激情 + 4 赞!
ter20 + 5 谢谢

查看全部评分 View all ratings

177

帖子

0

威望

1567

eV
积分
1744

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2017-3-15 15:19:58 | 只看该作者 Only view this author
pwzhou 发表于 2017-3-15 14:54
对09和之前的版本而言,当活化空间大于(8,8)的时候,优化交叉点和算频率时请用NoCPMCSCF关键词,另外,算 ...

感谢!

177

帖子

0

威望

1567

eV
积分
1744

Level 5 (御坂)

4#
 楼主 Author| 发表于 Post on 2017-3-15 18:21:32 | 只看该作者 Only view this author
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个虚频,感觉不太对

273

帖子

0

威望

4430

eV
积分
4703

Level 6 (一方通行)

5#
发表于 Post on 2017-4-6 16:27:24 | 只看该作者 Only view this author
pwzhou 发表于 2017-3-15 14:54
对09和之前的版本而言,当活化空间大于(8,8)的时候,优化交叉点和算频率时请用NoCPMCSCF关键词,另外,算 ...

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

368

帖子

0

威望

4691

eV
积分
5059

Level 6 (一方通行)

6#
发表于 Post on 2017-4-7 08:28:05 | 只看该作者 Only view this author
ter20 发表于 2017-3-15 18:21
发现很奇怪的结果,加了NoCPMCSCF和Freq=Numer后,Opt一步收敛,如下:
Item                           ...

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

177

帖子

0

威望

1567

eV
积分
1744

Level 5 (御坂)

7#
 楼主 Author| 发表于 Post on 2017-4-8 21:08:59 | 只看该作者 Only view this author
pwzhou 发表于 2017-4-7 08:28
我应该是记错了,算频率不需要加这个NoCPMCSCF关键词,用Freq=Numer即可。我很少用Gaussian做CASSCF,那 ...

多谢回复!

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-20 13:08 , Processed in 0.238730 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list