计算化学公社

标题: AVAS用于CASSCF轨道选取的问题 [打印本页]

作者
Author: BinWang 时间: 2026-1-14 10:51
标题: AVAS用于CASSCF轨道选取的问题
各位老师：
   大家好！我在计算一个六配位的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个占据轨道我都想考虑进去。这样的话可能会算不动。
   第二个问题是：大多数占据的活性轨道很难找到对应的空轨道
各位老师有什么建议？

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

作者
Author: BinWang 时间: 2026-1-14 16:27

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

非常感谢老师

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

作者
Author: BinWang 时间: 2026-1-15 09:09
本帖最后由 BinWang 于 2026-1-15 09:35 编辑

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

非常感谢Sob老师

，我积极争取报名参加，也会推送给需要的同学。

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