长链脂肪酸优化不收敛

chanyeolhgh

Sob老师您好，在长链脂肪酸C20H40O2的优化中，总是不收敛。看了您的解决SCF不收敛，以及DFT-D色散校正的文章后，尝试过int=ultrafine、SCF=NoVarAcc、SCF(novaracc,noincfock) int=acc2e=12、SCF=XQC、SCF=QC，但还是不收敛，求Sob老师指点，麻烦老师了，感谢老师。

biogon

脂肪酸阴离子为什么是开壳层的

chanyeolhgh

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

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

ldatea

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

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

chanyeolhgh

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

chanyeolhgh

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

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

biogon

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

原子数显然有问题

chanyeolhgh

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

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

ldatea

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也可以尝试一下。

chanyeolhgh

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

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

chanyeolhgh

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
  老师，输出信息不是很理解 这表示的是能量有减小的趋势吗  我还应该增加哪些关键词才会有帮助呢 感谢老师 跪求指点