计算化学公社

标题: orca片段分块波函数如何设置波电荷与自旋多重度 [打印本页]

作者
Author: BinWang 时间: 2025-12-26 21:09
标题: orca片段分块波函数如何设置波电荷与自旋多重度
本帖最后由 BinWang 于 2025-12-26 21:10 编辑

各位老师：
   大家好！我用orca计算一个体系需要用片段分块方法形成初猜波函数，高斯中可以很方便的给各个片段赋予不同的电荷与自旋多重度，组成一个波函数初猜。在orca6.1中，我在frag下定义charge以后报错，报错如下：171 ***********************************
172 * Starting time: Fri Dec 26 21:03:15 2025
173 * Host name:    node3
174 * Process ID: 64228
175 * Working dir.:  /share1/task/orca_play/ex_1
176 ***********************************
177
178 NOTE: MaxCore=6000 MB was set to SCF,MP2,MDCI,CIPSI,MRCI,RASCI and CIS
179    => If you want to overwrite this, your respective input block should be placed after the MaxCore statement
180 [file orca_tools/Tool-Scanner/qcscan1.cpp, line 106]:
181       Unknown identifier in Frag block line 16  :
182       Last token : CHARGE
说明chage和spin这个关键词定义的不对，这个该怎么解决？另外高斯中允许两个片段的alpha自旋定义为+3，beta自旋定义为-3，组合起来就是自旋多重度1见http://sobereva.com/82，这个在orca中如何实现？

这是我的输入文件，去掉charge和spin的定义，则正常结束，否则报错。
1 ! X2C X2C-TZVPall UHF TightSCF
  2 %pal
  3 nprocs 8
  4 end
  5 %maxcore 6000
  6 *xyz -2 2
  7 Cu 0.00 0.00 0.00
  8 Cl 2.25 0.00 0.00
  9 Cl -2.25 0.00 0.00
10 Cl 0.00 2.25 0.00
11 Cl 0.00 -2.25 0.00
12 *
13 %frag
14 Definition
15 1 {0} end # atom 0 for fragment 1
16 charge 2
17 Spin 1
18 2 {1:4} end # atoms 1 to 4 for fragment 2
19 charge -4
20 Spin 1
21 end
22 end

作者
Author: sobereva 时间: 2025-12-27 02:49
我不用这个功能，我只能说利用Multiwfn可以很容易地实现。Multiwfn的主功能100的子功能19可以把用户提供的片段的波函数文件组合产生整体的波函数文件，然后可以导出成各种记录波函数信息的文件格式，如mwfn/molden/fch等。参考Multiwfn手册3.100.19节的说明和4.100.19节的例子。

作者
Author: wzkchem5 时间: 2025-12-27 09:45
orca里面不是这么准备片段初猜的，需要用orca_mergefrag，参见https://www.faccts.de/docs/orca/ ... nd-list-of-keywords

作者
Author: BinWang 时间: 2025-12-27 10:00

wzkchem5 发表于 2025-12-27 09:45
orca里面不是这么准备片段初猜的，需要用orca_mergefrag，参见https://www.faccts.de/docs/orca/6.1/manual ...

非常感谢王老师回复！

作者
Author: BinWang 时间: 2025-12-27 10:04

sobereva 发表于 2025-12-27 02:49
我不用这个功能，我只能说利用Multiwfn可以很容易地实现。Multiwfn的主功能100的子功能19可以把用户提供的 ...

非常感谢Sob老师回复，事情是这样的:我要用sf-X2C处理一个很大的体系，用x2c-tzvpall，multiwfn恐怕难以实现这个功能。

作者
Author: sobereva 时间: 2025-12-28 14:22

BinWang 发表于 2025-12-27 10:04
非常感谢Sob老师回复，事情是这样的:我要用sf-X2C处理一个很大的体系，用x2c-tzvpall，multiwfn恐怕难以 ...

怎么难实现？

作者
Author: BinWang 时间: 2025-12-28 15:42
本帖最后由 BinWang 于 2025-12-29 08:48 编辑

sobereva 发表于 2025-12-28 14:22
怎么难实现？

非常感谢Sob老师回复，我尝试用Multiwfn生成片段波函数，然后送到orca里面作为初猜。不知道Multiwfn是否支持x2c-tzvpall这种基组？

作者
Author: sobereva 时间: 2025-12-30 14:41

BinWang 发表于 2025-12-28 15:42
非常感谢Sob老师回复，我尝试用Multiwfn生成片段波函数，然后送到orca里面作为初猜。不知道Multiwfn ...

Multiwfn对基组没有任何限制，只要是Gaussian型基函数就可以

作者
Author: BinWang 时间: 2025-12-30 14:45

sobereva 发表于 2025-12-30 14:41
Multiwfn对基组没有任何限制，只要是Gaussian型基函数就可以

非常感谢Sob老师

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