计算化学公社

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

[Gaussian/gview] 用CASSCF计算时,怎么确定轨道数和电子数,怎么计算CI?

[复制链接 Copy URL]

163

帖子

0

威望

1025

eV
积分
1188

Level 4 (黑子)

新手第一次用算CASSCF,1)、请问怎么确定分子active space的轨道数以及电子数?
2)、想用CASSCF算Conical Intersections,手册上写加关键词Opt=Conical,但是尝试了一下出现了一个报错

输入文件是#p CASSCF(4,6,StateAverage)/def2TZVP Opt=Conical freq
报错是  MCSCF will prepare Gradient Dif. and Deriv. Cps.
Enter MCSCF program.
NBasis=  690 NCore=   68 NVal=    6 NVirt=  616 IFrame=0 Acc=  1.00D-08.
End of file reading weights.
不知道该怎么解决

6万

帖子

99

威望

6万

eV
积分
125169

管理员

公社社长

2#
发表于 Post on 2019-6-4 21:51:48 | 只看该作者 Only view this author
你设了态平均就必须指定各个态的权重。论坛里之前讨论过,手册里也明确说了
这种任务不要同时带着freq,没有意义

活性空间选取问题论坛里专门有相关资料http://bbs.keinsci.com/thread-126-1-1.html。也可以照着类似的研究文献看他们都怎么设
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

163

帖子

0

威望

1025

eV
积分
1188

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2019-6-5 14:04:42 | 只看该作者 Only view this author
sobereva 发表于 2019-6-4 21:51
你设了态平均就必须指定各个态的权重。论坛里之前讨论过,手册里也明确说了
这种任务不要同时带着freq,没 ...

谢谢老师,请问各个态的权重在输入文件中应该怎么写呢?
我算了NRoot=3,在坐标数据下面写了

【坐标】

0.333333330.333333330.33333333

然后报错信息是 Wanted an integer as input.
Found a floating point number as input.
0.333333330.333333330.33333333                                                  
?
是我输入的各个态权重不对么?

4289

帖子

4

威望

9550

eV
积分
13919

Level 6 (一方通行)

MOKIT开发者

4#
发表于 Post on 2019-6-5 14:11:19 | 只看该作者 Only view this author
Abandon-fmt 发表于 2019-6-5 14:04
谢谢老师,请问各个态的权重在输入文件中应该怎么写呢?
我算了NRoot=3,在坐标数据下面写了

是你没把Title行、电荷行、权重行之间的空行写对。
另外,注意CASSCF一般直接做是不会收敛的,还是要先选好活性轨道然后进行计算,否则就算不报错,结果也一般没有意义。
自动做多参考态计算的程序MOKIT

163

帖子

0

威望

1025

eV
积分
1188

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2019-6-5 14:24:54 | 只看该作者 Only view this author
zjxitcc 发表于 2019-6-5 14:11
是你没把Title行、电荷行、权重行之间的空行写对。
另外,注意CASSCF一般直接做是不会收敛的,还是要先 ...

谢谢老师,请问选取活性轨道是要做一步计算还是可以直接根据分子结构进行选择?请问有没有例子可以借鉴一下?

4289

帖子

4

威望

9550

eV
积分
13919

Level 6 (一方通行)

MOKIT开发者

6#
发表于 Post on 2019-6-5 14:32:45 | 只看该作者 Only view this author
Abandon-fmt 发表于 2019-6-5 14:24
谢谢老师,请问选取活性轨道是要做一步计算还是可以直接根据分子结构进行选择?请问有没有例子可以借鉴一 ...

仔细看sob老师提供的链接(以及链接帖子里的链接)
自动做多参考态计算的程序MOKIT

163

帖子

0

威望

1025

eV
积分
1188

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2019-6-5 19:43:42 | 只看该作者 Only view this author
zjxitcc 发表于 2019-6-5 14:32
仔细看sob老师提供的链接(以及链接帖子里的链接)

好的好的,谢谢老师

本版积分规则 Credits rule

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

GMT+8, 2026-2-25 04:53 , Processed in 0.227813 second(s), 20 queries , Gzip On.

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