求助高斯优化离子液体与硝酸铕结合后期出现波动不收敛的问题

soda123 · 发表于 Post on 2022-9-27 09:32:05

小弟最近刚接触计算离子液体方面，在结构优化一个离子液体与硝酸铕结合时高斯总是不收敛，已经尝试更换了溶剂模型，但仍然在后期优化时上下波动巨大，望论坛内老师指导下面附上我的输入文件和高斯优化能量图

%NprocShared=40
%mem=160000mb
%chk=sucesseun1888dehp.chk
# opt=(cart,maxstep=5,notrust) freq genecp int=fine scrf=(solvent=n-dodecane,read) pbe1pbe
scf=(maxcycle=520,vshift,novaracc,xqc)
n1888
0 7
N -4.94050400 -0.40444400 -1.56212800
O -1.24281700 -1.25458500 -1.39469900
C -1.08410800 -2.22837000 -2.20938600
O -1.96287300 -2.80669500 -2.84890500
C 0.35647500 -2.69598800 -2.42007700
O 1.22127700 -1.93083900 -1.59783700
C 2.58590600 -2.08912800 -1.87464400
C 3.32325400 -1.11200500 -0.97223800
O 2.65880900 -0.35264300 -0.24511100
N 4.65755300 -1.07814900 -0.99179400
C -4.92096100 -1.33610900 -2.72572600
C -4.53875200 -1.14358800 -0.29025800
C -5.19740500 -2.50920000 -0.04182600
C -4.31260600 -3.72572800 -0.33521200
C -3.16252600 -3.96387500 0.64590700
C -3.59080700 -4.18470700 2.10763400
C -3.38213300 -2.97036300 3.01154000
C -3.84688800 -3.19430100 4.44689900
C -3.86901000 0.67028000 -1.76731000
C -4.03834000 1.64359200 -2.93502700
C -3.58872900 -1.99111700 5.34611800
C -4.84836400 2.93193400 -2.72551400
C -4.27764200 3.91887000 -1.70258400
C -2.95824700 4.56126400 -2.13848600
C -2.53906200 5.82023300 -1.37174200
C -2.03165600 5.63832600 0.06287100
C -3.09983600 5.39589900 1.12596900
C -6.27486200 0.28252200 -1.42768100
C -7.48024100 -0.53011200 -0.98603100
C -8.73219500 0.35264300 -1.04220500
C -9.99700300 -0.32668300 -0.49689300
C -10.31850900 -0.12167400 0.98801800
C -9.32628200 -0.68300400 2.01734000
C -8.22753900 0.28459700 2.45860200
C -7.28247500 -0.32495900 3.48828400
C 5.46963900 -1.94487500 -1.84429500
C 5.84874200 -3.28624800 -1.21155400
C 6.72196000 -3.21208200 0.04308700
C 8.12963800 -2.63974600 -0.15187700
C 9.04824000 -3.51593500 -1.00130900
C 10.44097400 -2.93855800 -1.26862700
C 11.30254400 -2.67434200 -0.02283500
C 11.21948400 -1.24908200 0.51870300
C 5.32889500 -0.03998200 -0.20273300
C 5.44701200 1.27449400 -0.97000100
C 6.12517100 2.37937700 -0.15973000
C 7.59386400 2.10940200 0.17004500
C 8.32025800 3.30206500 0.79385700
C 7.79329000 3.71905700 2.16491200
C 8.59310300 4.85789600 2.79175700
C 8.06415400 5.27623300 4.15852600
H 0.61922200 -2.56140300 -3.47802400
H 0.43797100 -3.76266700 -2.17540800
H 2.79207400 -1.85181100 -2.92962700
H 2.90588000 -3.12392800 -1.68979800
H -3.92634200 -1.79057800 -2.79751200
H -5.14969200 -0.77420100 -3.63116900
H -5.67568200 -2.10910600 -2.58769700
H -3.45343000 -1.25668200 -0.36271500
H -4.74651600 -0.43839900 0.51942100
H -6.12470100 -2.60739000 -0.61015200
H -5.50995100 -2.53607600 1.00913500
H -3.89763300 -3.64955000 -1.34768600
H -4.96225600 -4.61225700 -0.32809400
H -2.44079500 -3.13933100 0.59991600
H -2.61620400 -4.84463100 0.28538300
H -3.02077600 -5.02150600 2.53293100
H -4.64573500 -4.49993900 2.14519800
H -3.90301500 -2.09264700 2.59827800
H -2.31520600 -2.71224200 3.01030900
H -4.91966800 -3.43851000 4.45245600
H -3.33266000 -4.07568800 4.85655800
H -3.80004300 1.20245000 -0.81361900
H -2.94195700 0.10621400 -1.88574200
H -3.01080400 1.93038400 -3.18837100
H -4.40904100 1.11890600 -3.82428700
H -3.91640300 -2.17769000 6.37535500
H -4.12170700 -1.10377600 4.98164800
H -2.52067500 -1.74559700 5.37160500
H -5.88925500 2.70799800 -2.46351000
H -4.89499100 3.43901200 -3.69976300
H -4.16097400 3.42640700 -0.72823200
H -5.02471300 4.71056700 -1.54702300
H -3.05794500 4.83830400 -3.19958800
H -2.14507300 3.82651700 -2.09562400
H -3.37191800 6.54136800 -1.37217000
H -1.73387200 6.29338300 -1.95027500
H -1.48004500 6.54806800 0.33890600
H -1.29530400 4.82356100 0.09099500
H -2.65732500 5.42465100 2.12761200
H -3.88454200 6.16376500 1.07969500
H -3.57922700 4.41756800 1.01889300
H -6.11123300 1.10262300 -0.72085500
H -6.47423100 0.72376600 -2.40753100
H -7.33716100 -0.90984000 0.02940500
H -7.64026700 -1.39687700 -1.63943000
H -8.55585300 1.29245200 -0.49955300
H -8.89849400 0.63728600 -2.09025800
H -9.94773900 -1.40468000 -0.71091500
H -10.85781300 0.04740400 -1.06555400
H -11.29903800 -0.58355100 1.16205800
H -10.46078600 0.95272900 1.17752200
H -8.87801900 -1.61284700 1.63543100
H -9.88224400 -0.97828300 2.91776600
H -8.69963400 1.18294100 2.88053500
H -7.64935000 0.63489400 1.59314900
H -6.51539000 0.39027700 3.80544600
H -7.82938900 -0.64539100 4.38310500
H -6.76926300 -1.20859700 3.08665800
H 6.36627100 -1.37516600 -2.10883700
H 4.94125800 -2.12026000 -2.78723400
H 4.93015200 -3.83434500 -0.96501100
H 6.35162800 -3.87680400 -1.98836700
H 6.19326900 -2.63458100 0.81295300
H 6.81314300 -4.22813800 0.45016000
H 8.57670700 -2.50607000 0.84130100
H 8.08146300 -1.63333200 -0.59427300
H 8.57276300 -3.71794600 -1.97137200
H 9.15572500 -4.49376600 -0.50832000
H 10.34537400 -2.00737000 -1.84690900
H 10.96680100 -3.64392900 -1.92520700
H 12.35220900 -2.88494900 -0.26455200
H 11.03116700 -3.38985300 0.76677100
H 11.55235900 -0.52795200 -0.23781800
H 11.85887100 -1.12336600 1.40000100
H 10.19984200 -0.97267200 0.80719000
H 6.31180700 -0.42789500 0.07259600
H 4.75493500 0.10779600 0.71551800
H 6.00629300 1.10516900 -1.90188000
H 4.43855700 1.59632300 -1.25583800
H 5.55317200 2.54720800 0.76246900
H 6.05848700 3.31449100 -0.73237400
H 8.11572300 1.82089600 -0.75473000
H 7.67883800 1.24982100 0.85091300
H 8.26160700 4.15936500 0.10668200
H 9.38829800 3.05761400 0.88498700
H 7.80951500 2.84990800 2.84006300
H 6.74050500 4.02406900 2.08732700
H 8.58244800 5.72229900 2.11274400
H 9.64619700 4.55452300 2.88045100
H 8.09296800 4.43898800 4.86640900
H 7.02331800 5.61533700 4.09199200
H 8.65496700 6.09441700 4.58552300
N 0.03968300 1.79181500 -2.40599400
O 1.11027200 1.15545600 -2.15526300
O -0.06977900 2.49774000 -3.38722800
O -0.90426400 1.64009500 -1.57704300
N 0.32869500 -1.58185000 2.28781000
O -0.35517400 -1.90566100 1.26076900
O 0.29749000 -2.24208300 3.30391700
O 1.03901500 -0.54495200 2.17129100
N -0.68835600 2.06950900 1.54544400
O 0.45775700 2.13401400 1.01068200
O -1.10283200 2.93354100 2.28550700
O -1.38460700 1.04070700 1.26470300
Eu 0.27911000 -0.03000200 -0.14958300
C H O N 0
6-31g(d)
****
Eu 0
S 5 1.00
70059.4200000 0.000097
10776.2350000 0.000707
2482.4901000 0.003362
702.1526000 0.010206
216.7926000 0.016597
S 1 1.00
41.4950000 1.0
S 1 1.00
29.4969000 1.0
S 1 1.00
15.0309000 1.0
S 1 1.00
3.7772000 1.0
S 1 1.00
1.9219000 1.0
S 1 1.00
0.7131000 1.0
S 1 1.00
0.3097000 1.0
S 1 1.00
0.0549000 1.0
S 1 1.00
0.0227000 1.0
P 6 1.00
4058.0207000 0.000048
1023.3259000 0.000357
339.3852000 0.001687
122.3485000 0.005693
29.0558000 0.072398
20.6269000 -0.046807
P 1 1.00
14.6179000 1.0
P 1 1.00
5.4444000 1.0
P 1 1.00
2.7414000 1.0
P 1 1.00
1.2283000 1.0
P 1 1.00
0.5769000 1.0
P 1 1.00
0.2489000 1.0
P 1 1.00
0.0800000 1.0
D 6 1.00
389.0841000 0.000534
117.7602000 0.004258
44.9869000 0.017160
22.8599000 0.012427
7.8331000 0.270416
4.1134000 0.445272
D 1 1.00
2.0999000 1.0
D 1 1.00
1.0043000 1.0
D 1 1.00
0.3422000 1.0
D 1 1.00
0.1051000 1.0
F 5 1.00
128.2341000 0.003424
46.2685000 0.034313
20.7366000 0.123194
9.4479000 0.251824
4.3149000 0.351292
F 1 1.00
1.9121000 1.0
F 1 1.00
0.7897000 1.0
F 1 1.00
0.2859000 1.0
G 4 1.00
20.7366000 0.026516
9.4479000 0.099066
4.3149000 0.213378
1.9121000 0.453713
G 1 1.00
0.7897000 1.0
G 1 1.00
0.2859000 1.0
****
Eu 0
ECP28MWB 5 28
H-Komponente
1
2 1.000000 0.000000
S-H
1
2 23.471384 607.659331
P-H
1
2 16.772479 264.385476
D-H
1
2 13.981343 115.381375
F-H
1
2 23.962888 -49.400794
G-H
1
2 21.232458 -26.748273
surface=ses addsph

复制代码

wzkchem5 · 发表于 Post on 2022-9-27 15:44:10

本帖最后由 wzkchem5 于 2022-9-27 08:47 编辑

从能量最低的那个结构开始继续优化，此外考虑不要用opt=cart。
另外，这个计算必须加色散校正，因为你的体系明显是色散作用非常重要的体系，如果不加色散的话，配体的烷基“尾巴”摆的角度肯定不对，不加色散的结果是不能发表的。这个可能也和不收敛有一定关系，参见https://mattermodeling.stackexch ... structure/9668#9668

soda123 · 发表于 Post on 2022-9-27 22:12:02

wzkchem5 发表于 2022-9-27 15:44
从能量最低的那个结构开始继续优化，此外考虑不要用opt=cart。
另外，这个计算必须加色散校正，因为你的体 ...

谢谢您的回答！我再计算试一下

soda123 · 发表于 Post on 2022-9-29 14:59:05

wzkchem5 发表于 2022-9-27 15:44
从能量最低的那个结构开始继续优化，此外考虑不要用opt=cart。
另外，这个计算必须加色散校正，因为你的体 ...

老师您好，我按照您说的方法计算后在后期又出现了波动的问题，但波动的幅度明显减小，请问我现在是等待任务再算下去看下一步趋势还是再做一些修改再计算

soda123 · 发表于 Post on 2022-9-29 14:59:49

这是现在计算的图片

wzkchem5 · 发表于 Post on 2022-9-29 15:45:30

soda123 发表于 2022-9-29 07:59
老师您好，我按照您说的方法计算后在后期又出现了波动的问题，但波动的幅度明显减小，请问我现在是等待任 ...

检查一下能量波动是因为结构突变还是波函数突变

soda123 · 发表于 Post on 2022-9-29 21:54:25

结构突变老师在两个结构优化点之间我观察到高斯优化的结构出现了大的波动

soda123 · 发表于 Post on 2022-9-29 21:55:31

wzkchem5 发表于 2022-9-29 15:45
检查一下能量波动是因为结构突变还是波函数突变

结构突变老师，我观察高斯优化的两个突变点的优化结构出现了大的变化

wzkchem5 · 发表于 Post on 2022-9-29 22:34:13

soda123 发表于 2022-9-29 14:55
结构突变老师，我观察高斯优化的两个突变点的优化结构出现了大的变化

那可能就是高斯优化这种体系本来就比较低效。如果仍然没有收敛趋势的话，可以把任务kill掉，然后从能量最低的结构继续优化

soda123 · 发表于 Post on 2022-9-29 23:44:17

wzkchem5 发表于 2022-9-29 22:34
那可能就是高斯优化这种体系本来就比较低效。如果仍然没有收敛趋势的话，可以把任务kill掉，然后从能量最 ...

好的老师，谢谢您的回答！

soda123 · 发表于 Post on 2022-10-14 08:46:51

wzkchem5 发表于 2022-9-29 22:34
那可能就是高斯优化这种体系本来就比较低效。如果仍然没有收敛趋势的话，可以把任务kill掉，然后从能量最 ...

老师您好，我从能量最低开始算后出现了几何优化不收敛的问题，我也尝试了sob老师文章中的方法，但仍然出现后期震荡的情况，麻烦您解答一下下面是我的优化步数图和输出文件

soda123 · 发表于 Post on 2022-10-14 08:47:42

图式如下

soda123 · 发表于 Post on 2022-10-14 08:48:40

这个是log文件

wzkchem5 · 发表于 Post on 2022-10-14 14:23:50

soda123 发表于 2022-10-14 01:46
老师您好，我从能量最低开始算后出现了几何优化不收敛的问题，我也尝试了sob老师文章中的方法，但仍然出 ...

这种体系本来就难收敛，因为太柔性了。几个可以尝试的点
（1）输出文件明确说了不建议用GDIIS，因为从这个版本的高斯开始，GDIIS基本不会是最好的选择。试试关掉GDIIS
（2）可能初始Hessian质量不好，试试计算一下初始Hessian
（3）可能格点不够大（考虑到Eu比较重，需要比较大的径向格点，而把f电子积准可能对角向格点也有一定要求），试试增加格点

soda123 · 发表于 Post on 2022-10-18 16:02:18

wzkchem5 发表于 2022-10-14 14:23
这种体系本来就难收敛，因为太柔性了。几个可以尝试的点
（1）输出文件明确说了不建议用GDIIS，因为从这 ...

好的，谢谢老师！

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 求助高斯优化离子液体与硝酸铕结合后期出现波动不收敛的问题

浏览过的版块