计算化学公社

 找回密码 Forget password
 注册 Register
Views: 12165|回复 Reply: 3
打印 Print 上一主题 Last thread 下一主题 Next thread

[ORCA] 体系中含有重元素,用ECP后锕系原子SCF不收敛,致计算终止

[复制链接 Copy URL]

19

帖子

0

威望

1190

eV
积分
1209

Level 4 (黑子)

各位老师,大家好!最近尝试用ORCA计算一个含锕系元素的分子,遇到由于该锕系原子的SCF不收敛导致计算终止问题,请教解决方法,非常感谢。
ORCA 版本号:Version 4.0.1
附上inp文件及out文件,以及自己尝试过修改scf收敛标准的参数

######inp文件如下:
%pal
nprocs 4
end
%MaxCore 500
%scf print[p_mos] 1 end
! B3LYP Opt Freq
%basis  oldecp "SDD"
        newGTO O "6-31G(d,p)" end
end
* xyz 0 3
U       0.000000000      1.206538000     -0.070896000
U       0.000000000     -1.206538000     -0.070896000
O       0.000000000      0.000000000      1.630608000
*


#######out文件的尾部:
------------------------------
INITIAL GUESS: MODEL POTENTIAL
------------------------------
Loading Hartree-Fock densities                     ... done
  calling /home/lip/software/orca401/orca 3632_atom92.inp > 3632_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 4694, Process 0]:
Error encountered when trying to calculate the atomic fitting density!
[64:11286] *** Process received signal ***
[64:11286] Signal: Aborted (6)
[64:11286] Signal code:  (-6)
[file orca_scf/scfguess.cpp, line 4694, Process 0]:
Error encountered when trying to calculate the atomic fitting density!
[64:11286] [ 0] /lib64/libpthread.so.0[0x32bda0f710]
[64:11286] [ 1] /lib64/libc.so.6(gsignal+0x35)[0x32bd232925]
[64:11286] [ 2] /lib64/libc.so.6(abort+0x175)[0x32bd234105]
[64:11286] [ 3] /home/lip/software/orca401/orca_scf_mpi[0x717812]
[64:11286] [ 4] /home/lip/software/orca401/orca_scf_mpi[0x4cafba]
[64:11286] [ 5] /home/lip/software/orca401/orca_scf_mpi[0x423710]
[64:11286] [ 6] /home/lip/software/orca401/orca_scf_mpi[0x42e250]
[64:11286] [ 7] /home/lip/software/orca401/orca_scf_mpi[0x438068]
[64:11286] [ 8] /lib64/libc.so.6(__libc_start_main+0xfd)[0x32bd21ed1d]
[64:11286] [ 9] /home/lip/software/orca401/orca_scf_mpi(realloc+0x261)[0x422a19]
[64:11286] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 0 on node 64 exited on signal 6 (Aborted).
--------------------------------------------------------------------------

ORCA finished by error termination in SCF
Calling Command: mpirun -np 4  /home/lip/software/orca401/orca_scf_mpi 3632.gbw b
[file orca_tools/qcmsg.cpp, line 432]:
  .... aborting the run



########atom 92.out的尾部
---------------
SCF CONVERGENCE
---------------

  Last Energy change         ...   -4.8995e-04  Tolerance :   1.0000e-06
  Last MAX-Density change    ...    2.1600e-03  Tolerance :   1.0000e-05
  Last RMS-Density change    ...    5.5371e-05  Tolerance :   1.0000e-06
  Last DIIS Error            ...    8.5035e-03  Tolerance :   1.0000e-06

             **** DENSITY FILE WAS UPDATED (3632_atom92.scfp.tmp) ****
             **** ENERGY FILE WAS UPDATED (3632_atom92.en.tmp) ****

     --------------------------------------------------------------------
                                      WARNING
     The wavefunction IS NOT YET CONVERGED! It shows however signs of
     convergence. Therefore the wavefunction will be stored and can be
     used as input for another calculation.
     DO NOT USE THIS WAVEFUNCTION  FOR ANYHTING ELSE. It is NOT RELIABLE
     --------------------------------------------------------------------
-------
TIMINGS
-------
Total SCF time: 0 days 0 hours 0 min 26 sec

Total time                  ....      26.251 sec
Sum of individual times     ....      25.739 sec  ( 98.0%)

Fock matrix formation       ....      22.417 sec  ( 85.4%)
Diagonalization             ....       2.828 sec  ( 10.8%)
Density matrix formation    ....       0.148 sec  (  0.6%)
Population analysis         ....       0.000 sec  (  0.0%)
Initial guess               ....       0.003 sec  (  0.0%)
Orbital Transformation      ....       0.000 sec  (  0.0%)
Orbital Orthonormalization  ....       0.000 sec  (  0.0%)
DIIS solution               ....       0.342 sec  (  1.3%)

-------------------------   --------------------
FINAL SINGLE POINT ENERGY      -474.841553027447   (SCF not fully converged!)
-------------------------   --------------------

------------------------------------------------------------------------------
                                      ERROR
     This SCF-wavefunction IS NOT FULLY CONVERGED!
     You can't use it for properties or numerical calculations !
     Aborting the run ...
     Please restart calculation (with larger maxiter/different convergence flags)
------------------------------------------------------------------------------


#######自己尝试过用scf关键字修改收敛标准,但输出的结果和上面的一样
! Slowconv
! NoDIIS
! SloppySCF
%scf
  TolE 1e-8
  TolRMSP 1e-8
  TolMaxP 1e-6
  MaxIter 1500
  LShift  0.1  
  ShiftErr 0.1  
end


5

帖子

0

威望

81

eV
积分
86

Level 2 能力者

2#
发表于 Post on 2017-11-25 10:51:05 | 只看该作者 Only view this author
%basis  oldecp "SDD"
        newGTO O "6-31G(d,p)" end
你好 你写的基组的意思是第一行对应U  下面加入的是O

5

帖子

0

威望

81

eV
积分
86

Level 2 能力者

3#
发表于 Post on 2017-12-4 22:21:57 | 只看该作者 Only view this author
你加个 Grid4 NoFinalGrid  试一下

140

帖子

0

威望

823

eV
积分
963

Level 4 (黑子)

4#
发表于 Post on 2018-5-8 10:43:35 | 只看该作者 Only view this author
你好,请问你解决了吗

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-21 18:58 , Processed in 0.157360 second(s), 20 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list