计算化学公社

标题: molcas进行CASSCF 优化轨道道读取问题 [打印本页]

作者
Author: lwang 时间: 2025-4-11 22:27
标题: molcas进行CASSCF 优化轨道道读取问题
本帖最后由 lwang 于 2025-4-11 22:29 编辑

我使用molcas进行 PCM 下鸟嘌呤 casscf 结构优化，我casscf使用外部MP2 轨道道读取，所以我的第一步 casscf 应该读取 MP2 轨道，那么在下一个循环中，casscf 是继续会读取 MP2 磁道吗？如果未读取，它应该是什么轨道？
我的输入文件如下所示：
&GATEWAY
coord =g.XYZ
basis =Aug-cc-pvdz
group = C1
RF-input
PCM-model; solvent= water
End of RF-inputRICD
>>> EXPORT MOLCAS_MAXITER = 500
>>> DO WHILE
&SEWARD
DoAnalytical
&RASSCF
FILEORB =g.1.MP2Orb
Charge = 0
spin=1
nActEl=12 0 0
Ras2=12
ciroot =1 1 1
Rlxroot=1
RFROot=1
&SLAPAF
>>> END DO

作者
Author: Freeman 时间: 2025-4-13 14:33
本帖最后由 Freeman 于 2025-4-13 14:36 编辑

fileorb读轨道不该写在while循环里。因为几何优化的每一步casscf应该读前一步的收敛的轨道作初猜，而不是读mp2的轨道作初猜。所以应该在while前先独立地做一个读mp2轨道的casscf单点，而删掉while里的fileorb关键词。在独立的casscf单点任务中，记得挑选正确的活性空间；后续几何优化则不用再挑选一遍了。

运行时，独立的casscf单点做完后，当前路径就会出现一个和输入文件同名的rasorb文件，保存了casscf的轨道。后续优化，每一步不写fileorb的casscf就会自动读这个rasorb文件做初猜，并且用该步收敛的轨道覆盖掉上一步的rasorb。

另外，用mp2的自然轨道当casscf的初猜，应该没有什么特别的好处吧。不如用hf或者dft的轨道得了。

作者
Author: lwang 时间: 2025-4-14 13:02
本帖最后由 lwang 于 2025-4-14 13:03 编辑

Freeman 发表于 2025-4-13 14:33
fileorb读轨道不该写在while循环里。因为几何优化的每一步casscf应该读前一步的收敛的轨道作初猜，而不是读 ...

好的，谢谢您。这是我的输入文件，我做了两个rasscf，这样子是不是也可以的呢

&GATEWAY
coord =g.XYZ
basis =Aug-cc-pvdz
group = C1
RF-input
PCM-model; solvent= water
End of RF-input
RICD
>>> EXPORT MOLCAS_MAXITER = 500
>>> DO WHILE
&SEWARD
DoAnalytical
>>> If (iter = 1)
&RASSCF
FILEORB =g.1.MP2Orb
Charge = 0
spin=1
nActEl=12 0 0
Ras2=12
ciroot =1 1 1
Rlxroot=1
RFROot=1
>>> EndIf
>> COPY $Project.JobIph $Project.JobOld
&RASSCF
Charge = 0
spin=1
nActEl=12 0 0
Ras2=12
ciroot =1 1 1
Rlxroot=1
RFROot=1
&SLAPAF
>>> END DO

作者
Author: gravitylead 时间: 2025-4-15 11:45
本帖最后由 gravitylead 于 2025-4-17 17:31 编辑

可以考虑在循环前使用 EMIL 命令将输入轨道文件复制为 INPORB，随后使用 lumorb 读入轨道。
输入文件的大致格式如下：
&GATEWAY
***
>> COPY /path to your orb file/***.RasOrb INPORB

>> do while
&SEWARD
***
&RASSCF
***
lumorb
&SLAPAF
>> end do

作者
Author: lwang 时间: 2025-4-16 16:42

gravitylead 发表于 2025-4-15 11:45
可以考虑在循环前使用 EMIL 命令将输入轨道文件复制为 INPORB，随后使用 lumorb 读入轨道。
输入文件的大 ...

好的，谢谢。我尝试读入的轨道为mp2轨道，然后输入之后是报错的。
部分输入如下：
basis =Aug-cc-pvdz
group = C1
RICD
>> COPY a.4.MP2Orb INPORB
>>> EXPORT MOLCAS_MAXITER = 500
>>> DO WHILE
&SEWARD
DoAnalytical
&RASSCF
Charge = 0
spin=1
nActEl=12 0 0
Ras2=12
ciroot =1 1 1
Rlxroot=1
lumorb
&SLAPAF
>>> END DO
他显示的输出文件错误是：
>>> COPY a.4.MP2Orb INPORB

.####################.
.# EMIL input error #.

作者
Author: gravitylead 时间: 2025-4-17 17:24
本帖最后由 gravitylead 于 2025-4-17 17:31 编辑

lwang 发表于 2025-4-16 16:42
好的，谢谢。我尝试读入的轨道为mp2轨道，然后输入之后是报错的。
部分输入如下：
basis =Aug-cc-pvdz ...

需要输入轨道文件的绝对路径，让 molcas 知道这个文件在哪里。
我之前的回复可能没说清楚，现在已经做了修改。

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