计算化学公社

标题: 请教molpro里caspt2这部分的input [打印本页]

作者
Author: banxia 时间: 2023-1-31 18:19
标题: 请教molpro里caspt2这部分的input
各位老师好。我想利用molpro的caspt2进行几何优化，挑选的活性空间是homo和lumo这两个轨道，我的输入文件这部分内容如下，结果运行报错显示：too many active orbitals。
我有在论坛里搜过，看到过有老师解答说：caspt2的活性轨道个数指的是occ数减去core数。
但是我看我之前算成功的文件，看到我以前这样输入是没有报错的，下面是算成功的输出与输入（上面部分为输出里rs2这部分显示的内容，下面为输入文件的部分），这部分程序是可以成功识别出我的active轨道和core轨道的。
这部分不是很清楚究竟我的输入文件哪里有问题，还请老师指点，谢谢各位老师~

作者
Author: pwzhou 时间: 2023-1-31 19:37
这里的Active Orbitals应该指的是occ-core，47-14，可以把core增加1，然后试试

作者
Author: banxia 时间: 2023-1-31 19:46

pwzhou 发表于 2023-1-31 19:37
这里的Active Orbitals应该指的是occ-core，47-14，可以把core增加1，然后试试

core增加1我试过，是可以的。但是我不太明白提示为什么说增加1，报错信息为啥这样子？

还有就是之前这样子写去算其他东西没有报错，看输出文件里打印出来的内容，active的个数也就是在上一步casscf里面里的occ-closed的值。这回的报错我不懂了..或许您的话是什么意思，可以稍微具体说说吗

作者
Author: banxia 时间: 2023-1-31 20:07

pwzhou 发表于 2023-1-31 19:37
这里的Active Orbitals应该指的是occ-core，47-14，可以把core增加1，然后试试

而且，core设置成15的时候，我看输出识别出来的active orbitals个数也是2

作者
Author: pwzhou 时间: 2023-1-31 20:25

banxia 发表于 2023-1-31 19:46
core增加1我试过，是可以的。但是我不太明白提示为什么说增加1，报错信息为啥这样子？还有就是之前这 ...

没那么复杂，无非就是这里报错信息里面的Active Orbitals（occ-core）和通常CASSCF方法里面的Active Orbitals所指代的不一样而已，不用纠结这么多。这个限制应该是针对结构优化的，对于能量计算应该没有这个限制，所以你之前算其他东西报错。这个报错就是程序自身的限制所导致的，搞清楚限制的是哪一部分即可。

作者
Author: banxia 时间: 2023-2-1 10:00

pwzhou 发表于 2023-1-31 20:25
没那么复杂，无非就是这里报错信息里面的Active Orbitals（occ-core）和通常CASSCF方法里面的Active Orbi ...

好的，谢谢啊

作者
Author: banxia 时间: 2023-2-2 22:43

pwzhou 发表于 2023-1-31 20:25
没那么复杂，无非就是这里报错信息里面的Active Orbitals（occ-core）和通常CASSCF方法里面的Active Orbi ...

那个我想再问一下，我这样优化出来的几何结构，去做能量计算的时候，core的个数是原来的，还是报了一样的错，好像不仅仅针对结构优化..请问这是怎么回事啊

作者
Author: 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=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）

作者
Author: banxia 时间: 2023-2-5 12:01

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呀？不知道我的问题表达清楚了没有，非常感谢您！

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