计算化学公社

标题: 长链脂肪酸优化不收敛 [打印本页]

作者
Author: chanyeolhgh 时间: 2020-9-8 21:10
标题: 长链脂肪酸优化不收敛
Sob老师您好，在长链脂肪酸C20H40O2的优化中，总是不收敛。看了您的解决SCF不收敛，以及DFT-D色散校正的文章后，尝试过int=ultrafine、SCF=NoVarAcc、SCF(novaracc,noincfock) int=acc2e=12、SCF=XQC、SCF=QC，但还是不收敛，求Sob老师指点，麻烦老师了，感谢老师。
(, 下载次数 Times of downloads: 39) (, 下载次数 Times of downloads: 34)

作者
Author: biogon 时间: 2020-9-8 21:39
脂肪酸阴离子为什么是开壳层的

作者
Author: chanyeolhgh 时间: 2020-9-8 22:04

biogon 发表于 2020-9-8 21:39
脂肪酸阴离子为什么是开壳层的

老师含有奇数个电子不是开壳层的吗认为C20H40O2阴离子不是含有177个电子所以是开壳层的跪求老师指点感谢老师

作者
Author: ldatea 时间: 2020-9-8 22:07

chanyeolhgh 发表于 2020-9-8 22:04
老师含有奇数个电子不是开壳层的吗认为C20H40O2阴离子不是含有177个电子所以是开壳层的跪求老师指点 ...

你把氢原子数搞错了吧，应该是C20H39O2-。中性脂肪酸分子C20H40O2电离（-H+）得到C20H39O2-

作者
Author: chanyeolhgh 时间: 2020-9-9 07:33
Sob老师，我想算的是脂肪酸的电子亲和能，在优化和能量计算上一直都是在中性结构基础上只在输入文件中更改电荷数和自旋多重度，没考虑过电离出氢离子，我这样做是错误的是吗？感谢老师

作者
Author: chanyeolhgh 时间: 2020-9-9 07:33

ldatea 发表于 2020-9-8 22:07
你把氢原子数搞错了吧，应该是C20H39O2-。中性脂肪酸分子C20H40O2电离（-H+）得到C20H39O2-

Sob老师，我想算的是脂肪酸的电子亲和能，在优化和能量计算上一直都是在中性结构基础上只在输入文件中更改电荷数和自旋多重度，没考虑过电离出氢离子，我这样做是错误的是吗？感谢老师

作者
Author: biogon 时间: 2020-9-9 09:02

chanyeolhgh 发表于 2020-9-8 22:04
老师含有奇数个电子不是开壳层的吗认为C20H40O2阴离子不是含有177个电子所以是开壳层的跪求老师指点 ...

原子数显然有问题

作者
Author: chanyeolhgh 时间: 2020-9-9 11:18

biogon 发表于 2020-9-9 09:02
原子数显然有问题

Sob老师，我想算的是脂肪酸的电子亲和能，在优化和能量计算上一直都是在中性结构基础上只在输入文件中更改电荷数和自旋多重度，没考虑过原子数的问题，一直都是中性的原子数基础上更改电荷和自旋多重度，我这样做的一直都是错误的是吗？感谢老师，跪求老师指点

作者
Author: ldatea 时间: 2020-9-9 12:12
本帖最后由 ldatea 于 2020-9-9 12:26 编辑

chanyeolhgh 发表于 2020-9-9 07:33
Sob老师，我想算的是脂肪酸的电子亲和能，在优化和能量计算上一直都是在中性结构基础上只在输入文件中更 ...

原来你算的是自由基阴离子啊。算电子亲和能确实应该算自由基阴离子的能量。
你可以考虑算一下中性分子在B3LYP/6-31+g(d)下的能量，然后回到你的自由基阴离子用guess=read读取收敛的波函数作为初猜。
我看你输出文件里，Gap=0.005我觉得是HOMO-LUMO gap太小导致收敛困难，可以考虑尝试能级移动法SCF=vshift=x，x一般为300~500。SCF=Fermi也可以尝试一下。

作者
Author: chanyeolhgh 时间: 2020-9-9 13:52

ldatea 发表于 2020-9-9 12:12
原来你算的是自由基阴离子啊。算电子亲和能确实应该算自由基阴离子的能量。
你可以考虑算一下中性分子在 ...

好的老师感谢您我这就去试一试

作者
Author: chanyeolhgh 时间: 2020-9-13 10:16

ldatea 发表于 2020-9-9 12:12
原来你算的是自由基阴离子啊。算电子亲和能确实应该算自由基阴离子的能量。
你可以考虑算一下中性分子在 ...

Sob老师，我在尝试过int=ultrafine、SCF=NoVarAcc、SCF(novaracc,noincfock) int=acc2e=12、SCF=XQC、SCF=QC SCF=Vshift=300、SCF=Fermi 之后还是不收敛
输入信息
%mem=4000MB
%nprocshared=6
#p b3lyp/6-31+g(d) opt=tight  freq int=ultrafine geom=connectivity

Title Card Required

-1 2

输出信息
Cycle 128  Pass 1  IDiag  1:
RMSU=  8.18D-09 CP:  1.00D+00  1.00D+00 -9.54D-01  3.00D+00  3.00D+00
                  CP:  1.68D+00  3.00D+00
E= -935.958813925482    Delta-E=    -0.006074134781 Rises=F Damp=F
DIIS: error= 3.58D-02 at cycle 128 NSaved=  20.
NSaved=20 IEnMin= 1 EnMin= -936.446305629932    IErMin= 1 ErrMin= 1.60D-02
ErrMax= 3.58D-02  0.00D+00 EMaxC= 1.00D-01 BMatC= 1.14D+00 BMatP= 8.44D-01
IDIUse=2 WtCom= 0.00D+00 WtEn= 1.00D+00
Rare condition: small coef for last iteration:  0.000D+00
Coeff-En: 0.502D+00 0.000D+00 0.497D+00 0.000D+00 0.000D+00 0.000D+00
Coeff-En: 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00
Coeff-En: 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00
Coeff-En: 0.000D+00 0.130D-02
RMSU=  8.18D-09 CP:  1.00D+00  1.00D+00 -9.54D-01  3.00D+00  3.00D+00
                  CP:  1.68D+00  3.00D+00
E= -935.958813925482    Delta-E=    -0.006074134781 Rises=F Damp=F
DIIS: error= 3.58D-02 at cycle 128 NSaved=  20.
NSaved=20 IEnMin= 1 EnMin= -936.446305629932    IErMin= 1 ErrMin= 1.60D-02
ErrMax= 3.58D-02  0.00D+00 EMaxC= 1.00D-01 BMatC= 1.14D+00 BMatP= 8.44D-01
IDIUse=2 WtCom= 0.00D+00 WtEn= 1.00D+00
Rare condition: small coef for last iteration:  0.000D+00
RMSDP=1.86D-03 MaxDP=9.58D-02 DE=-6.07D-03 OVMax= 5.96D-03

Cycle 127  Pass 1  IDiag  1:
RMSU=  6.12D-10 CP:  1.00D+00  4.70D-05  1.78D+00 -3.11D-01  1.34D+00
                  CP:  1.44D-02
E= -935.952739790702    Delta-E=       0.006073984388 Rises=F Damp=F
DIIS: error= 3.58D-02 at cycle 127 NSaved=  20.
NSaved=20 IEnMin= 1 EnMin= -936.446305629932    IErMin= 1 ErrMin= 1.60D-02
ErrMax= 3.58D-02  0.00D+00 EMaxC= 1.00D-01 BMatC= 1.16D+00 BMatP= 8.44D-01
IDIUse=2 WtCom= 0.00D+00 WtEn= 1.00D+00
Coeff-En: 0.502D+00 0.000D+00 0.000D+00 0.196D+00 0.130D+00 0.141D+00
Coeff-En: 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00
Coeff-En: 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00
Coeff-En: 0.000D+00 0.315D-01
Coeff:    0.502D+00 0.000D+00 0.000D+00 0.196D+00 0.130D+00 0.141D+00
Coeff:    0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00
Coeff:    0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00
Coeff:    0.000D+00 0.315D-01
Gap=    0.236 Goal= None Shift= 0.000
Gap=    0.239 Goal= None Shift= 0.000
RMSDP=1.86D-03 MaxDP=9.58D-02 DE= 6.07D-03 OVMax= 8.13D-03

Cycle 128  Pass 1  IDiag  1:
RMSU=  8.18D-09 CP:  1.00D+00  1.00D+00 -9.54D-01  3.00D+00  3.00D+00
                  CP:  1.68D+00  3.00D+00
E= -935.958813925482    Delta-E=    -0.006074134781 Rises=F Damp=F
DIIS: error= 3.58D-02 at cycle 128 NSaved=  20.
NSaved=20 IEnMin= 1 EnMin= -936.446305629932    IErMin= 1 ErrMin= 1.60D-02
ErrMax= 3.58D-02  0.00D+00 EMaxC= 1.00D-01 BMatC= 1.14D+00 BMatP= 8.44D-01
IDIUse=2 WtCom= 0.00D+00 WtEn= 1.00D+00
Rare condition: small coef for last iteration:  0.000D+00
Coeff-En: 0.502D+00 0.000D+00 0.497D+00 0.000D+00 0.000D+00 0.000D+00
Coeff-En: 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00
Coeff-En: 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00
Coeff-En: 0.000D+00 0.130D-02
Coeff:    0.502D+00 0.000D+00 0.497D+00 0.000D+00 0.000D+00 0.000D+00
Coeff:    0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00
Coeff:    0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00
Coeff:    0.000D+00 0.130D-02
Gap=    0.233 Goal= None Shift= 0.000
Gap=    0.237 Goal= None Shift= 0.000
RMSDP=1.86D-03 MaxDP=9.58D-02 DE=-6.07D-03 OVMax= 5.96D-03

>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UB3LYP) =  -935.958813925    A.U. after  129 cycles
         NFock=128  Conv=0.19D-02    -V/T= 2.0035
<Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 0.5000 <S**2>= 0.7510 S= 0.5005
<L.S>= 0.000000000000E+00
KE= 9.326675769758D+02 PE=-5.305785533323D+03 EE= 1.880342144656D+03
Annihilation of the first spin contaminant:
S**2 before annihilation    0.7510, after    0.7500
Convergence failure -- run terminated.
Error termination via Lnk1e in /public/home/users/sues001/g09/l502.exe at Sat Sep 12 22:37:10 2020.
Job cpu time:    0 days  4 hours  9 minutes  3.2 seconds.
File lengths (MBytes):  RWF= 233 Int=    0 D2E=    0 Chk=    7 Scr=    1
  老师，输出信息不是很理解这表示的是能量有减小的趋势吗  我还应该增加哪些关键词才会有帮助呢感谢老师跪求指点

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