捐赠本论坛 (Donate this site)
开启辅助访问

计算化学公社

 找回密码 Forget password
 注册 Register
计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 量子化学 (Quantum Chemistry) › 长链脂肪酸优化不收敛

长链脂肪酸优化不收敛

查看数: 8247 | 评论数: 10 | 收藏 Add to favorites 0
关灯 | 提示：支持键盘翻页<-左 右->
帖子模式
    组图打开中，请稍候......
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，但还 ...

回复 Reply

chanyeolhgh 发表于 Post on 2020-9-13 10:16:05
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
  老师，输出信息不是很理解 这表示的是能量有减小的趋势吗  我还应该增加哪些关键词才会有帮助呢 感谢老师 跪求指点
chanyeolhgh 发表于 Post on 2020-9-9 13:52:08
ldatea 发表于 2020-9-9 12:12
原来你算的是自由基阴离子啊。算电子亲和能确实应该算自由基阴离子的能量。
你可以考虑算一下中性分子在 ...

好的 老师 感谢您 我这就去试一试
ldatea 发表于 Post on 2020-9-9 12:12:39
本帖最后由 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也可以尝试一下。
chanyeolhgh 发表于 Post on 2020-9-9 11:18:44
biogon 发表于 2020-9-9 09:02
原子数显然有问题

Sob老师，我想算的是脂肪酸的电子亲和能，在优化和能量计算上一直都是在中性结构基础上只在输入文件中更改电荷数和自旋多重度，没考虑过原子数的问题，一直都是中性的原子数基础上更改电荷和自旋多重度，我这样做的一直都是错误的是吗？感谢老师，跪求老师指点
biogon 发表于 Post on 2020-9-9 09:02:47
chanyeolhgh 发表于 2020-9-8 22:04
老师 含有奇数个电子不是开壳层的吗 认为C20H40O2阴离子不是含有177个电子 所以是开壳层的  跪求老师指点 ...

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

Sob老师，我想算的是脂肪酸的电子亲和能，在优化和能量计算上一直都是在中性结构基础上只在输入文件中更改电荷数和自旋多重度，没考虑过电离出氢离子，我这样做是错误的是吗？感谢老师
chanyeolhgh 发表于 Post on 2020-9-9 07:33:12
Sob老师，我想算的是脂肪酸的电子亲和能，在优化和能量计算上一直都是在中性结构基础上只在输入文件中更改电荷数和自旋多重度，没考虑过电离出氢离子，我这样做是错误的是吗？感谢老师
ldatea 发表于 Post on 2020-9-8 22:07:35
chanyeolhgh 发表于 2020-9-8 22:04
老师 含有奇数个电子不是开壳层的吗 认为C20H40O2阴离子不是含有177个电子 所以是开壳层的  跪求老师指点 ...

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

评分 Rate

参与人数
Participants 1		eV +2 收起 理由
Reason
zjxitcc + 2 我很赞同

查看全部评分 View all ratings

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

老师 含有奇数个电子不是开壳层的吗 认为C20H40O2阴离子不是含有177个电子 所以是开壳层的  跪求老师指点 感谢老师
biogon 发表于 Post on 2020-9-8 21:39:26
脂肪酸阴离子为什么是开壳层的

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

GMT+8, 2026-2-23 10:56 , Processed in 0.177034 second(s), 26 queries , Gzip On.

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