ORCA 5.0.2如何手动添加基组和ECP信息

neocc

请问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

hebrewsnabla

存在gtoaux=autobuild这个选项吗？

chrinide

要把默认的ECP删除，用delECP

neocc

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

neocc

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

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

chrinide

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

对头

neocc

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

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

chrinide

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

是的

neocc

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

wzkchem5

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

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

neocc

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

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

wzkchem5

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

在orca手册里搜索PAtom

neocc

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

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

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