计算化学公社

标题: 是否结构被认定错误导致不收敛问题？ [打印本页]

作者
Author: FFF 时间: 2021-1-6 09:00
标题: 是否结构被认定错误导致不收敛问题？
本帖最后由 FFF 于 2021-2-26 15:53 编辑

有一个问题想请教一下各位老师，最近在用高斯做这个分子的激发态计算：

尝试过不同种类和低精度的基组、用非限制性开壳层算，但都不行，总是502错误，第一步就无法收敛。我也是最近刚开始学习理论计算，很多不太懂，在想会不会是这种内盐正负电荷距离太远导致计算机认定结构有问题？想请教一下各位前辈会是什么原因导致的，要怎么解决呢？

作者
Author: TaylorWu 时间: 2021-1-6 10:21
SCF不收敛的问题可以参看社长写的博文：http://sobereva.com/61 如果实在不能解决，可能试试稍微调整一下初始构型？

作者
Author: zjxitcc 时间: 2021-1-6 11:10
本帖最后由 zjxitcc 于 2021-1-6 11:11 编辑

这种有机体系很好收敛。由于你是初学者，强烈建议贴出你的输入文件，让大家鉴定一下（初学者经常会写很多无用关键词、错误关键词，用奇葩泛函、基组，甚至电荷自旋多重度都写错，结构建模不合理等等），因此有必要进行正确的入门。初学者的错误纠正前，无需看sob老师的61号博文。

另外，由于你是初学者，应该先做基态结构优化，结构优化收敛后再做激发态计算，不要一下子就开始激发态计算。

作者
Author: Icestone 时间: 2021-1-6 11:19
不贴出具体的问题+不给计算的输入文件=让别人猜猜看=什么都没说

作者
Author: FFF 时间: 2021-1-6 14:55
不好意思，之前提供的参考信息可能不太够，现在补充一下，请各位前辈们帮忙看一看!
用过不同精度的基组：
%nprocshared=24
%mem=13GB
%chk=1.chk
#p pbe1pbe/def2svp opt td

%nprocshared=24
%mem=13GB
%chk=1.chk
#p b3lyp/6-311g(d,p) opt td

%NProcShared=24
%Mem=13GB
%Chk=1.chk
#p B3lyp/6-31G(d,p) Opt td

%NProcShared=24
%Mem=13GB
%Chk=1.chk
#p B3lyp/3-21G(d,p) Opt td

%NProcShared=24
%Mem=13GB
%Chk=1.chk
#p B3lyp/STO-3G Opt td

以最简单的基组为例，输入文件是：
%Chk=1.chk
%Mem=13GB
%NProcShared=24
#p uB3lyp/STO-3G Opt td

[No Title]

  0  1
C 0 -9.416735 -0.468678 1.777434
C 0 -8.467436 -0.253978 0.775034
C 0 -9.024036 -0.680278 3.112334
C 0 -7.674736 -0.681778 3.473634
C 0 -6.746235 -0.468878 2.475634
C 0 -7.139236 -0.258978 1.151134
C 0 -5.252436 -0.413078 2.540834
C 0 -4.884136 -0.147878 1.092434
N 0 -5.939835 -0.068578 0.355234
C 0 -5.929661 0.185722 -1.123018
C 0 -4.492000 0.332668 -1.603712
C 0 -4.481669 0.590867 -3.104630
S 0 -2.797633 0.762996 -3.667701
O 0 -2.787798 1.008819 -5.096677
O 0 -2.080162 2.037673 -2.882823
O 0 -2.065785 -0.455746 -3.382104
H 0  -10.486111 -0.472390 1.519709
H 0 -8.766051 -0.087842 -0.270542
H 0 -9.791284 -0.846993 3.882747
H 0 -7.363816 -0.846684 4.515812
H 0 -4.827063 -1.378404 2.895762
H 0 -4.861904 0.350050 3.250685
H 0 -3.854601 -0.037983 0.720975
H 0 -6.410713 -0.668927 -1.649264
H 0 -6.491418 1.121211 -1.342238
H 0 -4.010948 1.187316 -1.077466
H 0 -3.930243 -0.602822 -1.384491
H 0 -4.962721 -0.263781 -3.630876
H 0 -5.043426 1.526357 -3.323850

遇到的报错是：
Error: segmentation violation
rax 0000000000000000, rbx ffffffffffffffff, rcx 0000149a77322237
rdx 0000149776ed7000, rsp 00007ffef0eaf838, rbp 00007ffef0eafdb0
rsi 000000000000000b, rdi 0000000000002088, r8  0000000000000020
r9  0000000000000401, r10 0000000000000007, r11 0000000000000202
r12 00007ffef0eafdf8, r13 0000000000000000, r14 0000000000000000
r15 00000000000003e6
  /lib/x86_64-linux-gnu/libpthread.so.0(+0x128a0) [0x149a77a848a0]
  /lib/x86_64-linux-gnu/libc.so.6(kill+0x7) [0x149a77322237]
  /home/lushuang/g16/l502.exe() [0x4afeb9]
  /home/lushuang/g16/l502.exe() [0x4d1e0a]
  /home/lushuang/g16/l502.exe() [0x54e80c]
  /home/lushuang/g16/l502.exe() [0x52d7bf]
  /home/lushuang/g16/l502.exe() [0x4529db]
  /home/lushuang/g16/l502.exe() [0x450c2b]
  /home/lushuang/g16/l502.exe() [0x450b34]
  /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xe7) [0x149a77304b97]
  /home/lushuang/g16/l502.exe() [0x450a29]
Aborted (core dumped)

log文件的终止步骤：

SCF Done:  E(UB3LYP) =  -2413.59116366    A.U. after  129 cycles
         NFock=128  Conv=0.19D-03    -V/T= 2.0158
<Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 0.0000 <S**2>= 0.8739 S= 0.5601
<L.S>=  0.00000000000
KE= 2.376124017047D+03 PE=-1.478914716319D+04 EE= 5.427753299593D+03
Annihilation of the first spin contaminant:
S**2 before annihilation    0.8739, after    0.0057
Convergence failure -- run terminated.

作者
Author: FFF 时间: 2021-1-6 14:57

Icestone 发表于 2021-1-6 11:19
不贴出具体的问题+不给计算的输入文件=让别人猜猜看=什么都没说

不好意思，因为是初学者，想的没有这么周到，我补充了一些信息，想向前辈们再请教下！

作者
Author: FFF 时间: 2021-1-6 14:59

zjxitcc 发表于 2021-1-6 11:10
这种有机体系很好收敛。由于你是初学者，强烈建议贴出你的输入文件，让大家鉴定一下（初学者经常会写很多无 ...

感谢前辈的建议，我正在学习您推荐的内容。我又补充了一下信息，不收敛的问题还没有解决，想请您再看看有什么好的办法吗

作者
Author: elvisng 时间: 2021-1-6 15:09
本帖最后由 elvisng 于 2021-1-7 08:41 编辑

(內容已刪除)

作者
Author: zjxitcc 时间: 2021-1-6 15:25
本帖最后由 zjxitcc 于 2021-1-6 15:27 编辑

FFF 发表于 2021-1-6 14:59
感谢前辈的建议，我正在学习您推荐的内容。我又补充了一下信息，不收敛的问题还没有解决，想请您再看看有 ...

无需去尝试无意义的基组，见sob老师博文《谈谈量子化学中基组的选择》http://sobereva.com/336

另，你这是闭壳层的，不要写U。还是那句话，先算基态，再算激发态。
看你用B3LYP没加色散，再补充一点，看sob老师博文《谈谈“计算时是否需要加DFT-D3色散校正？”》http://sobereva.com/413

作者
Author: FFF 时间: 2021-1-6 17:22

TaylorWu 发表于 2021-1-6 10:21
SCF不收敛的问题可以参看社长写的博文：http://sobereva.com/61 如果实在不能解决，可能试试稍微调整一下初 ...

（不好意思刚才被限制回帖了）非常感谢您的建议。我简化过初始结构，我在想可能是因为烷基链太长，正负电荷距离太远导致的；但又不知道有什么好的解决办法

作者
Author: zjxitcc 时间: 2021-1-6 17:28
本帖最后由 zjxitcc 于 2021-1-6 17:29 编辑

FFF 发表于 2021-1-6 17:22
（不好意思刚才被限制回帖了）非常感谢您的建议。我简化过初始结构，我在想可能是因为烷基链太长，正负电 ...

“可能是因为烷基链太长，正负电荷距离太远导致的”完全没有关系。仔细、反复阅读楼上诸位的回复就行了。要是还有问题，请贴出你新的、有问题的输入文件，不要一次性贴那么多没用的基组之类的，就贴1个有用的、关键的。

作者
Author: FFF 时间: 2021-1-6 17:57
我认真学习一下各位前辈的问题和意见，再试试，非常感谢！

作者
Author: elvisng 时间: 2021-1-6 19:56
本帖最后由 elvisng 于 2021-1-7 08:41 编辑

(內容已刪除)

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