请选择 进入手机版 | 继续访问电脑版
第13届北京科音初级量子化学培训班将于10月5~8日于北京举办,请点击此链接查看详情。这是新人一次性正确、完整学习量子化学计算的最好、最快机会,能少走无数弯路,欢迎参加并相互转告!(已报满)

计算化学公社

 找回密码
 现在注册!
查看: 209|回复: 24

[新手求助] 长链脂肪酸优化不收敛

[复制链接]

27

帖子

0

威望

49

eV
积分
76

Level 2 能力者

发表于 7 天前 | 显示全部楼层 |阅读模式
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 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
  老师,输出信息不是很理解 这表示的是能量有减小的趋势吗  我还应该增加哪些关键词才会有帮助呢 感谢老师 跪求指点

749

帖子

5

威望

3422

eV
积分
4271

Level 6 (一方通行)

发表于 7 天前 | 显示全部楼层
解决SCF不收敛问题的方法
http://sobereva.com/61
Stand on the shoulders of giants

27

帖子

0

威望

49

eV
积分
76

Level 2 能力者

 楼主| 发表于 7 天前 | 显示全部楼层
冰释之川 发表于 2020-9-13 10:44
解决SCF不收敛问题的方法
http://sobereva.com/61

感谢老师 这个方法我看到了 也尝试了一些 量化新手 对输出信息不是很了解 能麻烦老师看一下我的输出文件中是有能量减少的趋势吗

247

帖子

0

威望

1298

eV
积分
1545

Level 5 (御坂)

发表于 7 天前 | 显示全部楼层
不要重复发帖。

你是opt的第一圈就SCF不收敛,还是优化若干圈之后才SCF不收敛的?

27

帖子

0

威望

49

eV
积分
76

Level 2 能力者

 楼主| 发表于 7 天前 | 显示全部楼层
hebrewsnabla 发表于 2020-9-13 10:56
不要重复发帖。

你是opt的第一圈就SCF不收敛,还是优化若干圈之后才SCF不收敛的?

好的 老师 我下次一定注意 从第一圈开始就没有SCF(Done)应该是第一圈就没收敛
输出信息
Cycle   1  Pass 1  IDiag  1:
FoFJK:  IHMeth= 1 ICntrl=       0 DoSepK=F KAlg= 0 I1Cent=   0 FoldK=F
IRaf= 690000000 NMat=   1 IRICut=       1 DoRegI=T DoRafI=F ISym2E= 0.
FoFCou: FMM=T IPFlag=           0 FMFlag=      100000 FMFlg1=        2001
         NFxFlg=           0 DoJE=F BraDBF=F KetDBF=F FulRan=T
         wScrn=  0.000000 ICntrl=       0 IOpCl=  1 I1Cent=           0 NGrid=           0
         NMat0=    1 NMatS0=      1 NMatT0=    0 NMatD0=    1 NMtDS0=    0 NMtDT0=    0
Symmetry not used in FoFCou.
FMM levels:  10  Number of levels for PrismC:   9
E= -935.675188807761
DIIS: error= 4.08D-02 at cycle   1 NSaved=   1.
NSaved= 1 IEnMin= 1 EnMin= -935.675188807761     IErMin= 1 ErrMin= 4.08D-02
ErrMax= 4.08D-02  0.00D+00 EMaxC= 1.00D-01 BMatC= 3.57D+00 BMatP= 3.57D+00
IDIUse=3 WtCom= 5.92D-01 WtEn= 4.08D-01
Coeff-Com:  0.100D+01
Coeff-En:   0.100D+01
Coeff:      0.100D+01
Gap=    -0.123 Goal=   None    Shift=    0.000
Gap=     0.198 Goal=   None    Shift=    0.000
GapD=   -0.123 DampG=0.250 DampE=0.500 DampFc=0.2500 IDamp=-1.
Damping current iteration by 2.50D-01
RMSDP=7.43D+00 MaxDP=7.71D+02              OVMax= 9.95D-01
老师 这样的情况有什么办法吗 感谢老师

247

帖子

0

威望

1298

eV
积分
1545

Level 5 (御坂)

发表于 7 天前 | 显示全部楼层
chanyeolhgh 发表于 2020-9-13 11:06
好的 老师 我下次一定注意 从第一圈开始就没有SCF(Done)应该是第一圈就没收敛
输出信息
Cycle   1   ...

你说你尝试了这些关键词,但是贴出的输入文件中并没有体现。到底有没有使用,使用这些关键词之后到底有没有效果(结果离收敛更近)呢?

有效果的话可以多管齐下。

27

帖子

0

威望

49

eV
积分
76

Level 2 能力者

 楼主| 发表于 7 天前 | 显示全部楼层
hebrewsnabla 发表于 2020-9-13 11:14
你说你尝试了这些关键词,但是贴出的输入文件中并没有体现。到底有没有使用,使用这些关键词之后到底有没 ...

不好意思 老师 是我的问题 这次贴出来的只用了int=ultrafine 但之前尝试过很多关键词 结果都是502错误 老师 具体怎么看离收敛更近呢(老师别生气 小的问的问题太低级了 但真心的不太懂)跪求老师指点

247

帖子

0

威望

1298

eV
积分
1545

Level 5 (御坂)

发表于 7 天前 | 显示全部楼层
chanyeolhgh 发表于 2020-9-13 11:30
不好意思 老师 是我的问题 这次贴出来的只用了int=ultrafine 但之前尝试过很多关键词 结果都是502错误 老 ...
  1. RMSDP=1.86D-03 MaxDP=9.58D-02 DE=-6.07D-03 OVMax= 5.96D-03
复制代码

这一行,看RMSDP,MaxDP,DE是否接近收敛标准。

749

帖子

5

威望

3422

eV
积分
4271

Level 6 (一方通行)

发表于 7 天前 | 显示全部楼层
hebrewsnabla 发表于 2020-9-13 12:28
这一行,看RMSDP,MaxDP,DE是否接近收敛标准。

离收敛限很远
Stand on the shoulders of giants

27

帖子

0

威望

49

eV
积分
76

Level 2 能力者

 楼主| 发表于 7 天前 | 显示全部楼层

老师 您说的收敛限是指最大受力、均方根受力、最大位移和均方根位移吗 就是下面的这个吗
Maximum Force            0.000001     0.000015     YES
RMS     Force            0.000000     0.000010     YES
Maximum Displacement     0.000030     0.000060     YES
RMS     Displacement     0.000007     0.000040     YES
还是RMSDP、MAXDP、DE的值有规定的收敛标准呢?那具体的收敛标准的值是多少呢或者我应该在哪里看呢 感谢老师赐教

749

帖子

5

威望

3422

eV
积分
4271

Level 6 (一方通行)

发表于 7 天前 | 显示全部楼层
chanyeolhgh 发表于 2020-9-13 13:48
老师 您说的收敛限是指最大受力、均方根受力、最大位移和均方根位移吗 就是下面的这个吗
Maximum Force  ...

Requested convergence on RMS density matrix=1.00D-08 within 128 cycles.
Requested convergence on MAX density matrix=1.00D-06.
Requested convergence on             energy=1.00D-06.
No special actions if energy rises.
Fock matrices will be formed incrementally for  20 cycles.
Integral accuracy reduced to 1.0D-05 until final iterations.
Stand on the shoulders of giants

131

帖子

0

威望

148

eV
积分
279

Level 3 能力者

发表于 7 天前 | 显示全部楼层
chanyeolhgh 发表于 2020-9-13 13:48
老师 您说的收敛限是指最大受力、均方根受力、最大位移和均方根位移吗 就是下面的这个吗
Maximum Force  ...

不要混淆 SCF是否收敛 和 几何优化是否收敛 ,这是两个概念

2万

帖子

25

威望

3万

eV
积分
58547

管理员

公社社长

发表于 6 天前 | 显示全部楼层
本来这种柔性体系优化就相对来说不容易收敛,还写tight完全是自找麻烦。

直接把优化任务的输出文件压缩后上传,什么情况别人一看就清楚了,省得你描述不清楚,别人猜来猜去也猜不中你的真实情况,白讨论半天。
北京科音自然科学研究中心http://www.keinsci.com  致力于计算化学的发展和传播,长期开办最高水准的各种量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训,是提升计算化学研究水平的最佳选择。欢迎加入“北京科音”公众号获取培训最新消息和计算化学资讯!培训相关信息见《北京科音办的培训班FAQ》(http://bbs.keinsci.com/thread-5098-1-1.html)。
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群:思想家公社QQ群1号:18616395,2号:466017436。合计6000人。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一定会被拒绝加入。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

27

帖子

0

威望

49

eV
积分
76

Level 2 能力者

 楼主| 发表于 6 天前 | 显示全部楼层
sobereva 发表于 2020-9-14 06:12
本来这种柔性体系优化就相对来说不容易收敛,还写tight完全是自找麻烦。

直接把优化任务的输出文件压缩 ...

感谢Sob老师 给答疑的老师们带来不便深感抱歉 下次一定注意 但老师如果不写tight的话 那关键词int=ultrafine还可以用吗 对其他解决SCF不收敛的关键词的添加有影响吗 我在优化后出现了虚频 按照Sob老师优化解决虚频的帖子上的方法 gview虚频模式下微调结构(调了0.01) 或加入opt=calcall 第一种微调后还是有虚频 第二种直接不收敛 出现502错误 请问老师该怎么解决呢 都是要把tight去掉吗
第一种输入文件:
1ce5333675ea332e3e7fb570b0e2cb2.png
输出文件 12.4.7z (133.38 KB, 下载次数: 2)

93

帖子

1

威望

480

eV
积分
593

Level 4 (黑子)

发表于 6 天前 | 显示全部楼层
chanyeolhgh 发表于 2020-9-14 13:02
感谢Sob老师 给答疑的老师们带来不便深感抱歉 下次一定注意 但老师如果不写tight的话 那关键词int=ultraf ...

你第二个输出文件没正确贴上来
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2020-9-20 22:39 , Processed in 0.176518 second(s), 27 queries .

快速回复 返回顶部 返回列表