AVAS用于CASSCF轨道选取的问题

BinWang

各位老师：
       大家好！我在计算一个六配位的Ir络合物的旋轨耦合常数，该络合物的一个配体是双自由基，我用AVAS投影出想要的活性轨道。投影的的目标包含两个部分，第一个部分是双自由基相关的轨道，第二部分是Ir的d轨道，我做这个计算的目的是希望对比一下MOKIT计算出的结果与orca AVAS方法输出的结果的差异。我的输入文件如下：
     ! X2C X2C-SVPall NoIter
%pal                                   
nprocs 32                             # 第3行
end                                    
%maxcore 6000                         # 第5行：单行指令（无end！）
%scf                                   # 第6行：SCF参数块
MOInp "stable.gbw"                 # 第7行
Guess MORead
avas
system  #Ir d轨道
shell 5, 5, 5, 5, 5
l 2, 2, 2, 2, 2
m_l 0, 1, -1, 2, -2
center 0, 0, 0, 0, 0
end
system   #双自由基中心1，该自由基离域在苯环
  center 6, 8, 20, 25, 67, 75, 63
  type pz, pz, pz, pz, pz, pz, pz
end
system
  shell 2, 2, 2   #双自由基中心2
  l 1, 1, 1       # 角量子数l均为1，代表p轨道
  m_l 1, -1, 0    # 磁量子数：1 -> px, -1 -> py, 0 -> pz (此为常见约定)
  center 82, 82, 82
end
end
end                                    # 第8行：闭合%scf
%cpcm smd true
      SMDsolvent "water"
end
%basis                                 # 第12行：基组设置
  NewGTO Ir "X2C-TZVPall" end          # 第13行
end                                    # 第14行
*xyzfile 1 3 geo.xyz                   # 第15行：分子坐标文件

我计算出AVOS轨道以后遇到了两个问题：
      第一个问题是：双自由基轨道，Ir的d轨道都已经明显的显示出来了，还有双自由基与Ir形成的sigma-轨道，但这些轨道的数目太大，有12个占据轨道我都想考虑进去。这样的话可能会算不动。
       第二个问题是：大多数占据的活性轨道很难找到对应的空轨道
各位老师有什么建议？

-HMDS-

第一的话建议再筛选一下活性轨道，12个确实有点多，感觉上没有什么必要，再加上orca的casscf算法确实效率不高。
第二点的话个人操作方式，可能不是很合理但是有时候能解决问题的方法，一个是用hf算的轨道或者小cas的做初猜，因为dft做多参考的时候空轨道确实比较乱，再不就是限制对称性。

BinWang

-HMDS- 发表于 2026-1-14 16:18
第一的话建议再筛选一下活性轨道，12个确实有点多，感觉上没有什么必要，再加上orca的casscf算法确实效率不 ...

非常感谢老师

sobereva

ORCA的CASSCF效率已经算非常高了，尤其是默认的superCI_PT算法，而且本身ORCA的电子积分效率又高，CASSCF时还支持RI加速。ORCA的ICE-SCF（面对大活性空间的CASSCF形式）也很成熟，CAS(30,30)在普通服务器上都能搞得动，北京科音高级量子化学培训班（http://www.keinsci.com/KAQC）给了这个方面很详细的讲解和例子。

BinWang

sobereva 发表于 2026-1-15 01:52
ORCA的CASSCF效率已经算非常高了，尤其是默认的superCI_PT算法，而且本身ORCA的电子积分效率又高，CASSCF时 ...

非常感谢Sob老师，我积极争取报名参加，也会推送给需要的同学。