计算化学公社

标题: 怎样在orca中直接对轨道进行localization并且重新投影得到新的CI和compositions？ [打印本页]

作者
Author: sxh3231796 时间: 2025-8-28 22:49
标题: 怎样在orca中直接对轨道进行localization并且重新投影得到新的CI和compositions？
我是用orca计算casscf收敛成功，并输出对应的CI和compositions，但是要怎么样在orca中直接对轨道进行localization并且重新投影得到新的CI 和compositions.要是用什么样的关键词啊
! def2-tzvp SlowConv

%moinp "Fe-cation-orca-orbital.gbw"

%casscf
  nel 4
  norb  4
  maxiter 1
  PrintWF csf
  TPrintWF 0.01
end

%loc
LocMet    PM
end

我是用这个做了localization，但是怎么进一步投影啊

作者
Author: hebrewsnabla 时间: 2025-8-28 23:39
actorbs locorbs

作者
Author: sxh3231796 时间: 2025-8-29 00:15

hebrewsnabla 发表于 2025-8-28 23:39
actorbs locorbs

谢谢老师，能在问您一下吗，这个是我要再做了localization之后用这个关键词再不迭代的重新做一边casscf吗

作者
Author: sobereva 时间: 2025-8-29 05:39
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并清楚、准确反映出帖子具体内容，避免有任何歧义和含糊性，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题 “ORCA 计算casscf” 改了，以后务必注意

作者
Author: sobereva 时间: 2025-8-29 05:41

sxh3231796 发表于 2025-8-29 00:15
谢谢老师，能在问您一下吗，这个是我要再做了localization之后用这个关键词再不迭代的重新做一边casscf吗

参考>=第二届北京科音高级量子化学培训班（http://www.keinsci.com/KAQC）详细讲ORCA做CASSCF和多参考计算部分的幻灯片：
(, 下载次数 Times of downloads: 1)

完全不需要你自己写
%loc
LocMet PM
end

常规casscf计算指定了actorbs locorbs后，自动在迭代收敛后对活性轨道做定域化并存到gbw里，并且最后给你的CI系数也是基于定域化轨道的。例如培训里计算1,4-脱氢苯的双自由基态这个例子：

(, 下载次数 Times of downloads: 1)

作者
Author: sxh3231796 时间: 2025-8-29 15:10

sobereva 发表于 2025-8-29 05:41
参考>=第二届北京科音高级量子化学培训班（http://www.keinsci.com/KAQC）详细讲ORCA做CASSCF和多参考计 ...

谢谢sob老师，那如果想用PM做定域化需要什么关键词啊，我昨天用了上面那位老师给的关键词套在我的inp里做了计算，但是输出的CI和composition和NBO 没有区别呀
! def2-TZVP SlowConv

%base "Fe-casscf-local-print"    # 可选：统一输出前缀
%moinp "Fe-cation-orca-orbital.gbw" # 读取你上一次收敛的MOs/密度/几何

%casscf
  nel    4
  norb 4
  maxiter 1                # 不再优化，只做一次检查与打印（稳妥）
  ActOrbs LocOrbs
  PrintWF       csf       # 打印 CSF/配置
  ci
TPrintWF 0.005
  end
end

作者
Author: sobereva 时间: 2025-8-29 17:51

sxh3231796 发表于 2025-8-29 15:10
谢谢sob老师，那如果想用PM做定域化需要什么关键词啊，我昨天用了上面那位老师给的关键词套在我的inp里做 ...

哪来的NBO，别突然冒出来一个无关的东西

actorbs产生定域化轨道只能用FB

actorbs做轨道定域化会影响CI系数，在我讲的培训班的例子幻灯片里已经体现得没法更清楚了

作者
Author: sxh3231796 时间: 2025-8-29 18:06
本帖最后由 sxh3231796 于 2025-8-29 18:07 编辑

sobereva 发表于 2025-8-29 17:51
哪来的NBO，别突然冒出来一个无关的东西

actorbs产生定域化轨道只能用FB

就是natrural orbitals，在我不加关键词之前，直接输出的CI 和我加了关键词输出的CI 完全一致，那请问老师还有用PM localization的工具吗？或者我用multiwfn做定域化轨道，然后用别的软件做波函数映射可以吗？我之前的同事使用molpro做的，其中有PM定域化，但是我们没有这个软件了

作者
Author: sobereva 时间: 2025-8-30 05:10

sxh3231796 发表于 2025-8-29 18:06
就是natrural orbitals，在我不加关键词之前，直接输出的CI 和我加了关键词输出的CI 完全一致，那请问老 ...

NBO和natrural orbitals (NO)有天壤之别

没有具体输入输出文件，没法判断怎么“一致”

不理解为什么非要用PM定域化不可。FB一般足以达到同样目的

原理上可以把Multiwfn专门对活性轨道混合产生的LMO连同内、外空间的正则轨道给ORCA作为初猜波函数然后做CASSCF，届时用actorbs unchanged避免迭代收敛后产生自然轨道，且用ActConstraints 0（避免默认的superCI_PT算法自动在迭代过程每一轮产生自然轨道）。

作者
Author: sxh3231796 时间: 2025-8-30 16:35
本帖最后由 sxh3231796 于 2025-8-30 16:36 编辑

sobereva 发表于 2025-8-30 05:10
NBO和natrural orbitals (NO)有天壤之别

没有具体输入输出文件，没法判断怎么“一致”

               <<<<<<<<<<<<<<<<<<INITIAL CI STATE CHECK>>>>>>>>>>>>>>>>>>

BLOCK  0 MULT= 1 NROOTS= 1
ROOT 0:  E= -3565.1578898217 Eh
   0.89140 [    0]: 2200
   0.08592 [    1]: 2110
   0.00638 [ 15]: 0202
   0.00563 [    2]: 2101
   0.00514 [    7]: 1201

               <<<<<<<<<<<<<<<<<<INITIAL CI STATE CHECK>>>>>>>>>>>>>>>>>>
---------------------------------------------
CAS-SCF STATES FOR BLOCK  0 MULT= 1 NROOTS= 1
---------------------------------------------

ROOT 0:  E= -3565.1578898217 Eh
   0.89140 [    0]: 2200
   0.08592 [    1]: 2110
   0.00638 [ 15]: 0202
   0.00563 [    2]: 2101
   0.00514 [    7]: 1201
因为我之前已经casscf收敛了，但是加了关键词后，得到的composition 和CI 和initial check（没加关键词，用natural orbital 计算）没有区别啊

%base "Fe-casscf-local-print"    # 可选：统一输出前缀
%moinp "Fe-cation-orca-orbital.gbw" # 读取你上一次收敛的MOs/密度/几何

%casscf
  nel    4
  norb 4
  maxiter 1                # 不再优化，只做一次检查与打印（稳妥）
  ActOrbs LocOrbs
  PrintWF       csf       # 打印 CSF/配置
  ci
TPrintWF 0.005
  end
end
这个是我的输入文件，想问下老师有什么问题吗？

作者
Author: sobereva 时间: 2025-8-31 07:51

sxh3231796 发表于 2025-8-30 16:35
BLOCK  0 MULT= 1 NROOTS= 1
ROOT 0:  E= -3565.1578898217 Eh
   0.8914 ...

直接上传完整输出文件
看输出信息里的提示，弄清楚到底这次做了轨道定域化没有

作者
Author: sxh3231796 时间: 2025-9-1 13:43

sobereva 发表于 2025-8-31 07:51
直接上传完整输出文件
看输出信息里的提示，弄清楚到底这次做了轨道定域化没有

sob老师，这个是我的输出文件

作者
Author: sobereva 时间: 2025-9-2 09:41

sxh3231796 发表于 2025-9-1 13:43
sob老师，这个是我的输出文件

搜local并没搜到做定域化的信息，大抵因为你写了maxiter 1导致定域化没做。去掉maxiter 1
做轨道定域化的时候会有这种信息：
(, 下载次数 Times of downloads: 2)

另外，建议加上nopop避免输出轨道系数，要不然文件忒大，输出那些系数也没意义，有gbw转出来的molden文件随时都可以在Multiwfn里显示轨道系数

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