molcas进行CASSCF 优化轨道道读取问题

lwang

我使用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

Freeman

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


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

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

lwang

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

gravitylead

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

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

lwang

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 #.

gravitylead

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

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