计算化学公社

标题: ORCA能否将两个片段（金属离子和配体）的gbw文件合成一个新的gbw文件？ [打印本页]

作者
Author: CLJJ 时间: 2023-10-10 20:34
标题: ORCA能否将两个片段（金属离子和配体）的gbw文件合成一个新的gbw文件？
各位老师好，我在分析锕系金属离子（使用全电子基组）与有机配体形成的配合物的电子构型的过程中，需要将配合物分成两个片段，一个片段是金属离子，一个片段是有机配体。基于Gaussian优化的配合物结构，将金属离子以及有机配体的坐标提取出来，用orca对两个片段分别进行单点能计算，之后对金属离子的片段调整轨道占据初猜方式，得到两个波函数gbw文件。

请问能否在orca中将这两个波函数文件直接组合在一起产生一个新的波函数gbw文件呢？使得其可以用于后续的计算。

（如果将两个gbw文件分别转为molden文件，再转为fchk文件，在Multiwfn软件中可以将两个fchk文件合在一起生成一个新的fchk文件，再将新的fchk文件转成mkl文件，再在orca中转成gbw文件，再用orca计算可以实现上述目的，但是由于我计算的是锕系金属离子，电子构型比较复杂，所以担心如此繁琐的过程会产生一些问题。）

作者
Author: wzkchem5 时间: 2023-10-10 21:14
可以，orca有一个工具叫做orca_mergefrag，在orca手册里搜一下就能找到用法。
此外如果最终目的只是算配合物，不是一定需要两个片段各自的波函数的话，可以改用flipspin等关键字一步得到想要的配合物波函数，而未必需要用片段初猜方法

作者
Author: CLJJ 时间: 2023-10-10 21:23
本帖最后由 CLJJ 于 2023-10-10 23:11 编辑

wzkchem5 发表于 2023-10-10 21:14
可以，orca有一个工具叫做orca_mergefrag，在orca手册里搜一下就能找到用法。
此外如果最终目的只是算配合 ...

谢谢老师的回复，我去阅读了orca手册的相关部分，发现orca_mergefrag关键词好像不适用于将一个开壳层片段与一个闭壳层片段进行组合，但是我的金属离子是开壳层，配体是闭壳层；并且我想构造的初猜就是将金属离子的一个5f占据轨道与一个6d未占据轨道进行调换，即将金属离子的外层价电子构型从5f5变成5f46d1，配体的电子构型保持不变，这与使用对称性破缺关键词Flipspin的情况不太相同，请问老师现在这种情况该如何解决呀？

作者
Author: sobereva 时间: 2023-10-11 01:09
倘若支持把两个开壳层的波函数合并，把闭壳层的那个当开壳层计算就完了

靠Multiwfn完全可以实现你的目的。交换两个轨道在主功能6里用29 Exchange information of two orbitals。主功能6里37 Split spatial orbitals as alpha and beta spin orbitals还可以将原本的闭壳层波函数转换成alpha和beta相同的开壳层波函数。

作者
Author: CLJJ 时间: 2023-10-11 10:18

sobereva 发表于 2023-10-11 01:09
倘若支持把两个开壳层的波函数合并，把闭壳层的那个当开壳层计算就完了

靠Multiwfn完全可以实现你的目的 ...

好的，谢谢老师的回复。

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