计算化学公社

 找回密码 Forget password
 注册 Register

请教molpro里caspt2这部分的input

查看数: 2637 | 评论数: 8 | 收藏 Add to favorites 2
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2023-1-31 18:19

正文摘要:

各位老师好。我想利用molpro的caspt2进行几何优化,挑选的活性空间是homo和lumo这两个轨道,我的输入文件这部分内容如下,结果运行报错显示:too many active orbitals。 我有在论坛里搜过,看到过有老师解答说:ca ...

回复 Reply

banxia 发表于 Post on 2023-2-5 12:01:41
qinjiu 发表于 2023-2-5 00:16
你之前计算成功的是用了C2v点群。core设置是(4,0,2,0),occ是(11,7,7,4),active orb=occ-core= 7+7+5+4 ...

好的好的,我明白了,这里molpro的active是包括occ里的一部分的,我思维固定了,所以一直不理解。那请问您这里说的两种解决办法的的第二种是什么意思呀,是molpro说明里面geometry的Symmetry specification这一部分吗?设置这个为什么可以解决呀。我还想问一下,利用df-rs2来对分子进行几何优化的时候,因为有df,好像优化的分子不能有symmetry,但是我目标的是c2v,请问df-rs2可以把c1的初始结构优化到c2v吗?因为我目前尝试了好多遍,优化出来的结构虽然是c2v接近,但是还是c1。是我初始结构的问题,还是它就只能继续优化到c1呀?不知道我的问题表达清楚了没有,非常感谢您!
qinjiu 发表于 Post on 2023-2-5 00:16:34
你之前计算成功的是用了C2v点群。core设置是(4,0,2,0),occ是(11,7,7,4),active orb=occ-core= 7+7+5+4=23
最新算的没有设置对称性,所以是C1点群。core设置是14,occ是47, active orb= occ-core=33
molpro的CASPT2最多支持计算32个active orb ,超过就会提示 too many active orbitals
你要么调整occ和core卡,要么设置体系的对称性 (symmetry,x,y)

评分 Rate

参与人数
Participants 2
eV +10 收起 理由
Reason
banxia + 5 谢谢
王二葛 + 5 大佬学通了 molpro

查看全部评分 View all ratings

banxia 发表于 Post on 2023-2-2 22:43:59
pwzhou 发表于 2023-1-31 20:25
没那么复杂,无非就是这里报错信息里面的Active Orbitals(occ-core)和通常CASSCF方法里面的Active Orbi ...

那个我想再问一下,我这样优化出来的几何结构,去做能量计算的时候,core的个数是原来的,还是报了一样的错,好像不仅仅针对结构优化..请问这是怎么回事啊
banxia 发表于 Post on 2023-2-1 10:00:03
pwzhou 发表于 2023-1-31 20:25
没那么复杂,无非就是这里报错信息里面的Active Orbitals(occ-core)和通常CASSCF方法里面的Active Orbi ...

好的,谢谢啊
pwzhou 发表于 Post on 2023-1-31 20:25:04
banxia 发表于 2023-1-31 19:46
core增加1我试过,是可以的。但是我不太明白提示为什么说增加1,报错信息为啥这样子?还有就是之前这 ...

没那么复杂,无非就是这里报错信息里面的Active Orbitals(occ-core)和通常CASSCF方法里面的Active Orbitals所指代的不一样而已,不用纠结这么多。这个限制应该是针对结构优化的,对于能量计算应该没有这个限制,所以你之前算其他东西报错。这个报错就是程序自身的限制所导致的,搞清楚限制的是哪一部分即可。
banxia 发表于 Post on 2023-1-31 20:07:26
pwzhou 发表于 2023-1-31 19:37
这里的Active Orbitals应该指的是occ-core,47-14, 可以把core增加1,然后试试

而且,core设置成15的时候,我看输出识别出来的active orbitals个数也是2
banxia 发表于 Post on 2023-1-31 19:46:53
pwzhou 发表于 2023-1-31 19:37
这里的Active Orbitals应该指的是occ-core,47-14, 可以把core增加1,然后试试

core增加1我试过,是可以的。但是我不太明白提示为什么说增加1,报错信息为啥这样子?还有就是之前这样子写去算其他东西没有报错,看输出文件里打印出来的内容,active的个数也就是在上一步casscf里面里的occ-closed的值。这回的报错我不懂了..或许您的话是什么意思,可以稍微具体说说吗
pwzhou 发表于 Post on 2023-1-31 19:37:00
这里的Active Orbitals应该指的是occ-core,47-14, 可以把core增加1,然后试试

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

GMT+8, 2026-2-22 20:51 , Processed in 0.204675 second(s), 26 queries , Gzip On.

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