计算化学公社

标题: ORCA 5.0.2如何手动添加基组和ECP信息 [打印本页]

作者
Author: neocc 时间: 2022-2-22 21:29
标题: ORCA 5.0.2如何手动添加基组和ECP信息
本帖最后由 neocc 于 2022-2-22 21:36 编辑

请问ORCA 5.0.2里面如何手动添加https://www.basissetexchange.org/提供的基组和ECP信息?

按照 ORCA 5.0.2手册9章5节选择基组，尝试如下inp文件：

! B3LYP D3 ma-def2-TZVP autoaux RIJCOSX opt freq tightSCF noautostart miniprint nopop
%maxcore  2000
%pal nprocs  64 end
%cpcm
smd true
SMDsolvent "water"
end
%basis
  NewGTO 92
S 7
   1    1534.9336000             0.0002230
   2    227.7483800             0.0007220
   3       30.6968310          -0.0217020
   4       18.1706260             0.1497960
   5       10.8135370          -0.2761100
   6       2.7332980             0.3765790
   7       1.4314980             0.1244080
S 7
   1    1534.9336000          -0.0006950
   2    227.7483800          -0.0022930
   3       30.6968310             0.0724440
   4       18.1706260          -0.4905590
   5       10.8135370             0.8898400
   6       2.7332980          -1.1157580
   7       1.4314980          -0.3336500
S 7
   1    1534.9336000             0.0009200
   2    227.7483800             0.0036820
   3       30.6968310          -0.1502310
   4       18.1706260             0.9039950
   5       10.8135370          -1.4859220
   6       2.7332980             1.1244040
   7       1.4314980             0.3510760
S 1
   1       0.6152980             1.0000000
S 1
   1       0.2866390             1.0000000
S 1
   1       0.0711700             1.0000000
S 1
   1       0.0305390             1.0000000
S 1
   1       0.0050000             1.0000000
P 6
   1    553.3452500          -0.0010750
   2    109.2550100          -0.0080420
   3       23.4760300          -0.0279690
   4       6.7944720             0.4690490
   5       5.4323190          -0.4926390
   6       2.7021690          -0.3757660
P 6
   1    553.3452500          -0.0016100
   2    109.2550100          -0.0110870
   3       23.4760300          -0.0466190
   4       6.7944720             0.6441530
   5       5.4323190          -0.5204550
   6       2.7021690          -0.6311540
P 1
   1       1.4938570             1.0000000
P 1
   1       0.7928170             1.0000000
P 1
   1       0.3515420             1.0000000
P 1
   1       0.1439620             1.0000000
P 1
   1       0.0050000             1.0000000
D 6
   1       81.2028580          -0.0006390
   2       18.3255750          -0.0070210
   3       10.4546990             0.0266420
   4       3.6663120          -0.1571320
   5       1.9233490          -0.1575720
   6       0.9896380          -0.0011000
D 6
   1       81.2028580          -0.0016200
   2       18.3255750          -0.0251810
   3       10.4546990             0.0896320
   4       3.6663120          -0.4094290
   5       1.9233490          -0.4822600
   6       0.9896380          -0.2109030
D 1
   1       0.4953460             1.0000000
D 1
   1       0.2044550             1.0000000
D 1
   1       0.0732730             1.0000000
D 1
   1       0.0050000             1.0000000
F 5
   1       55.3345250             0.0012940
   2       16.5886490             0.0098570
   3       4.7575180             0.1303860
   4       2.3875500             0.3180830
   5       1.1301950             0.3953140
F 1
   1       0.4895350             1.0000000
F 1
   1       0.1814200             1.0000000
F 1
   1       0.0050000             1.0000000
  end
  NewECP 92
N_core 60
lmax h
s 1
   1    16.414038690  536.516627780 2
p 1
   1    9.060556060  169.544924650 2
d 1
   1    8.831831980  142.615598370 2
f 1
   1    7.018516290 60.393076020 2
g 1
   1    12.804088440  -60.129989580 2
h 1
   1    1.000000000 0.000000000 2
  end
  gtoaux=autobuild
end
* xyz 0 1
[GEOMETRY]
*

运行后报错，提示:

[file orca_main/maininp1.cpp, line 3240]: [BASIS]: New ORCA basis convention: Please use quote delimited "..." specifiers!

将基组和ECP信息用引号""括起来同样是上述报错。

  NewGTO 92
"xxxx"
  end
  NewECP 92
   "xxxx"
  end

作者
Author: hebrewsnabla 时间: 2022-2-22 22:30
存在gtoaux=autobuild这个选项吗？

作者
Author: chrinide 时间: 2022-2-22 23:19
要把默认的ECP删除，用delECP

作者
Author: neocc 时间: 2022-2-22 23:38

hebrewsnabla 发表于 2022-2-22 22:30
存在gtoaux=autobuild这个选项吗？

抱歉，gtoaux=autobuild是老版本的内容，删除后之前的报错内容不再出现了，这个问题已解决
不过再次计算后提示
ORCA finished by error termination in GTOInt

参考如下帖子，可能是超算slurm的openmpi设置有问题
http://bbs.keinsci.com/thread-19186-1-1.html

作者
Author: neocc 时间: 2022-2-22 23:43

chrinide 发表于 2022-2-22 23:19
要把默认的ECP删除，用delECP

%basis
  NewGTO 92
   "xxxx"
  end
  NewECP 92
   "xxxx"
  end
  DelECP 92 #是这样修改嚒？
end

作者
Author: chrinide 时间: 2022-2-23 08:18

neocc 发表于 2022-2-22 23:43
%basis
NewGTO 92
"xxxx"

对头

作者
Author: neocc 时间: 2022-2-23 09:12

chrinide 发表于 2022-2-23 08:18
对头

  DelECP 92  #需要按照先删后立的顺序吗
  NewECP 92
   "xxxx"
  end

作者
Author: chrinide 时间: 2022-2-23 11:18

neocc 发表于 2022-2-23 09:12
DelECP 92  #需要按照先删后立的顺序吗
  NewECP 92
   "xxxx"

是的

作者
Author: neocc 时间: 2022-2-23 16:14

chrinide 发表于 2022-2-23 11:18
是的

终于开始跑起来了,结果还没高兴30s,又报错了:

         ************************************************************
         *       Program running with 64 parallel MPI-processes *
         *             working on a common directory             *
         ************************************************************

Diagonalization of the overlap matrix:
Smallest eigenvalue                      ... 2.956e-06
Time for diagonalization                ... 0.052 sec
Threshold for overlap eigenvalues       ... 1.000e-08
Number of eigenvalues below threshold    ... 0
Time for construction of square roots    ... 0.267 sec
Total time needed                         ... 0.490 sec

Time for model grid setup = 1.009 sec

Loading Hartree-Fock densities                   ... done
  calling orca UO2-orca_atom92.inp > UO2-orca_atom92.out in order to generate an atomic fitting density for atom 0 (U  with ECP) on-the-fly...
atom 0 (U ), assumed electronic state with S=5: 1s2 2s2 2p6 3s2 3p6 4s2 3d10 4p6 5s2 4d10 5p6 6s2 4f14 5d10 6p6 7s2 5f4          ... [file orca_scf/scfguess.cpp, line 4853, Process 0]:
Error encountered when trying to calculate the atomic fitting density!

ORCA finished by error termination in SCF
Calling Command: mpirun -np 64  orca_scf_mpi UO2-orca.gbw b UO2-orca
Check for SCF-logfiles
[file orca_main/maingopt.cpp, line 1214]: ORCA finished with an error in the energy calculation

作者
Author: wzkchem5 时间: 2022-2-23 16:26

neocc 发表于 2022-2-23 09:14
终于开始跑起来了,结果还没高兴30s,又报错了:

PModel不支持U等比较重的元素，应该改用PAtom等初始猜测

作者
Author: neocc 时间: 2022-2-23 16:46

wzkchem5 发表于 2022-2-23 16:26
PModel不支持U等比较重的元素，应该改用PAtom等初始猜测

请问具体操作需要怎么设置呢?

作者
Author: wzkchem5 时间: 2022-2-23 17:00

neocc 发表于 2022-2-23 09:46
请问具体操作需要怎么设置呢?

在orca手册里搜索PAtom

作者
Author: neocc 时间: 2022-2-23 22:52

wzkchem5 发表于 2022-2-23 17:00
在orca手册里搜索PAtom

感谢，先试一下手册里面的！PAtom

结果还是SCF不收敛，看来是此贴之外的另一个困难了

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