计算化学公社

标题: 求助Molcas程序搜索S0/T1，S0/T2等各种态的MECP时，MCLR子程序不收敛 [打印本页]

作者
Author: xexlalalan 时间: 2021-3-8 13:27
标题: 求助Molcas程序搜索S0/T1，S0/T2等各种态的MECP时，MCLR子程序不收敛
在OpenMolcas程序中利用MECP搜索对氯苯酚的基态与各个激发态的圆锥交叉点时，执行MCLR子程序时，res(CI)总是不收敛。已尝试各种初始结构，均告失败。请问大家这种问题原因大概有哪些？该如何解决这个问题？

另：在Molcas的官方论坛上找到过一个类似的问题，提问者怀疑是软件编译的问题，不过这个问题到现在也没有人解答……（https://cobalt.itc.univie.ac.at/molcasforum/viewtopic.php?id=728）

附件中是我的输入及输出文件
如能解答，不胜感激！！！

作者
Author: 函数与激情 时间: 2021-3-8 18:55
估计是PCM导致的问题，RASSCF模块应该加上RFPERT关键词吧？

也有可能Molcas不支持PCM优化CI；

GATEWAY中把PCM关键词删除可以正常计算。

作者
Author: xexlalalan 时间: 2021-3-9 10:00

函数与激情发表于 2021-3-8 18:55
估计是PCM导致的问题，RASSCF模块应该加上RFPERT关键词吧？

也有可能Molcas不支持PCM优化CI；

非常感谢！我把PCM去掉之后确实成功了。但是我这个体系是在液相中进行反应的，去掉溶剂模型结果又不敢用。我尝试了一下对溶剂部分进行修改，换成C-PCM也不行。之前放宽收敛限也不知道为什么对态平均的那一次MCLR不起作用，感觉现在真是没辙了……

另：RFPERT好像是加在CASPT2里的？

作者
Author: zjxitcc 时间: 2021-3-9 10:25
非要加上隐式溶剂的话，可以尝试高斯CASSCF（轨道初猜质量要保证，最好用GVB轨道，次好用UNO轨道；优化步长调小一点；CASSCF圈数调大一些）。

作者
Author: 函数与激情 时间: 2021-3-9 18:24

函数与激情发表于 2021-3-8 18:55
估计是PCM导致的问题，RASSCF模块应该加上RFPERT关键词吧？

也有可能Molcas不支持PCM优化CI；

RFPERT可以加在RASSCF模块的；

Module  rasscf  key  RFPErt

  This keyword will add a constant reaction field
perturbation to the  bare nuclei Hamiltonian. The
perturbation is read from the RUNOLD (if not present
defaults to RUNFILE) and is the latest self consistent
perturbation generated  by one of the programs SCF or
RASSCF.

作者
Author: xexlalalan 时间: 2021-3-12 14:40

函数与激情发表于 2021-3-9 18:24
RFPERT可以加在RASSCF模块的；

Module rasscf key RFPErt

非常感谢回复！抱歉我前几天在做别的project，所有没来的及顾上这个。我这两天又试了一下，把第一轮RASSCF抽出Do while循环，然后把RFPERT加到Do while循环里面的RASSCF可以跑，但是依然有MCLR不收敛的问题。我现在在尝试降低MCLR的收敛限，但是不知道为什么只对前几个MCLR管用，后面的又回到默认收敛限了。
其它的一些思路也尝试过，比如换成CPCM溶剂；活性空间、基组也换过，均收效甚微。不过有几次RASSCF之后软件有警告：WARNING: Large orbital rotation，我隐隐感觉可能还是活性空间选的不好的原因。我再试试

作者
Author: xexlalalan 时间: 2021-3-12 14:44

zjxitcc 发表于 2021-3-9 10:25
非要加上隐式溶剂的话，可以尝试高斯CASSCF（轨道初猜质量要保证，最好用GVB轨道，次好用UNO轨道；优化步长 ...

谢谢！我用的是UNO做初猜，一个(12,10)的活性空间。之前就是用Gaussian算这个CASSCF光报错才转投到了Molcas的怀抱。我再钻一钻Molcas，实在不行再回去

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