计算出错，显示Inv3 failed in PCMMkU，怎么办？

guyifan

想询问各位大佬一个问题，最近在做铁矿物吸附金属酸根离子溶液构型，该构型加入了显式溶剂模型（水分子）、隐式溶剂模型（iefpcm）以及D3色散校正，但是能量迟迟不收敛还出了Inv3 failed in PCMMkU (L502/L508 报错)以下报错：并根据大佬的帖子做出了以下尝试，但是依旧不能正常运行，想请各位大佬帮忙查看一下。输入文件附下：

http://bbs.keinsci.com/thread-6682-1-1.html（参考网址）

**********************************************************
输入文件：

%nprocs=48
%chk=ccgoethite_monodentate_b3lyp_6-311+G2dfp_6w_D3+iefpcm_2.chk
%mem=100000MB
# opt=(gdiis,z-matrix,maxcycle=1000) b3lyp/genecp nosymm charge scf=xqc empiricaldispersion=gd3bj scrf=(smd,solvent=iefpcm,read)

Si-OH4_6w

-59 96
O        0        26.5370270        -10.3979490        5.4412930
O        0        26.6160530        -10.2375900        8.3863150
O        0        27.6699860        -7.9590870        3.7858480
O        0        27.7925680        -7.9029830        6.7700960
O        0        27.7402330        -7.9277910        9.7778870
O        0        29.1641690        -10.1585080        3.7719780
O        0        29.1417380        -10.2356080        6.7071360
O        0        29.2297680        -10.2039550        9.7748710
O        0        30.2172670        -8.2260150        5.3289190
O        0        30.1969230        -8.3333650        8.2154440
O        0        26.0142520        -13.0837560        6.4838770
O        0        28.4878960        -12.5285990        5.1533230
O        0        25.9822430        -12.6002790        3.6260490
H        0        29.1619500        -11.0753770        4.0633980
H        0        29.0941470        -11.1681280        6.4478490
H        0        29.0493780        -11.1332880        9.6264050
H        0        30.2308400        -7.3646470        5.7583940
H        0        30.1664030        -7.4207620        7.9203190
H        0        26.4727550        -10.5087550        7.4698450
O        0        26.7140490        -10.2424820        1.9664660
O        0        26.8630470        -10.3385760        11.3351380
O        0        22.9020170        -8.3855390        6.7749190
O        0        32.8161850        -12.9681650        6.7646320
O        0        25.2950740        -7.8672190        5.2363530
O        0        25.3092240        -7.8822570        8.3012460
O        0        24.2160010        -5.9404550        6.7879460
H        0        25.6042810        -8.7986070        5.2377520
H        0        25.6718360        -8.8101130        8.3458490
H        0        24.1372150        -4.9768370        6.8901870
H        0        23.6785160        -8.9820570        6.7773450
H        0        26.5577610        -10.9987430        2.5816630
H        0        26.6636110        -10.5150420        10.3872120
H        0        31.9499290        -13.4225340        6.7136400
O        0        31.6324540        -10.5502450        5.2321290
O        0        31.6848250        -10.5417310        8.2943060
O        0        30.1832430        -8.3075510        2.3040070
O        0        30.2315130        -8.3157680        11.2272650
O        -1        26.6882910        -5.4905150        5.2655910
O        -1        26.7026040        -5.5084490        8.2759010
O        -1        32.8401940        -8.2970820        6.7716350
O        -1        29.1982970        -5.5546540        6.7722080
H        -1        28.7551600        -6.4544800        6.7982800
O        -1        25.3252780        -7.7502710        2.2576200
O        -1        25.3253060        -7.7502460        11.2881330
O        -1        35.2450000        -7.7502650        5.2678180
O        -1        35.2450280        -7.7502750        8.2780530
O        -1        35.2450330        -12.3302880        5.2678460
O        -1        35.2450430        -12.3302870        8.2780710
O        -1        34.1935420        -5.9897440        3.7627580
O        -1        34.1935420        -5.9897530        9.7831270
O        -1        34.1935490        -10.5697500        3.7627470
O        -1        34.1935440        -10.5697500        9.7831390
O        -1        29.2335920        -5.4602840        3.7628410
O        -1        29.2336870        -5.4601460        9.7835040
O        -1        29.2336160        -10.0402600        0.7525450
O        -1        29.2336480        -10.0402330        12.7933440
O        -1        26.7586120        -5.5395940        2.2577050
O        -1        26.7585720        -5.5395930        11.2880230
O        -1        32.7599580        -3.6205880        6.7726650
O        -1        27.8002910        -3.2496160        3.7629800
O        -1        27.8002920        -3.2494010        9.7834160
O        -1        27.8004220        -7.8296120        0.7527300
O        -1        27.8003260        -7.8296040        12.7932470
O        -1        21.7988390        -5.9105630        5.2678430
O        -1        21.7988840        -5.9104890        8.2781320
O        -1        31.7184230        -10.4903260        2.2577910
O        -1        31.7184390        -10.4904010        11.2881600
O        -1        22.8403110        -8.2004830        3.7627510
O        -1        22.8403640        -8.2004780        9.7830220
O        -1        32.7601390        -12.7805090        3.7627950
O        -1        32.7601480        -12.7805020        9.7831410
O        -1        24.2130800        -5.9533900        3.7477070
O        -1        24.2625140        -5.9437340        9.8289570
O        -1        34.2394510        -5.9688280        6.7471830
O        -1        34.2577610        -10.5449480        6.7708340
O        -1        30.3199540        -3.6567690        5.2472020
O        -1        30.3199830        -3.6552340        8.2957910
O        -1        32.7699770        -8.2101400        3.7574530
O        -1        32.8142540        -8.1923120        9.7825710
O        -1        27.8159620        -3.1789700        6.7712730
O        -1        31.7366880        -5.8889910        5.2437780
O        -1        31.7157420        -5.8953820        8.3199740
Fe        0        23.2736060        -7.1068030        5.1486930
Fe        0        23.2905810        -7.1132520        8.4084380
Fe        -1        33.1569350        -7.1605250        5.2513450
Fe        -1        33.1613070        -7.1692570        8.2994110
Fe        0        33.2264640        -11.7319040        5.1914290
Fe        0        33.2385560        -11.7506940        8.3499190
Fe        0        26.1658090        -6.6262390        3.7282430
Fe        0        26.1515780        -6.6492710        6.7799010
Fe        0        26.1988410        -6.6341420        9.8286650
Fe        -1        31.2385880        -4.8538140        6.7750910
Fe        0        31.1592790        -9.3926730        3.6892050
Fe        0        31.1802960        -9.5213270        6.7750770
Fe        0        31.1903190        -9.4049180        9.8476330
Fe        -1        28.2119540        -4.3386200        5.2673080
Fe        -1        28.2118640        -4.3338130        8.2715210
Fe        0        28.2373330        -8.9722980        2.2046950
Fe        0        28.3369750        -8.9254450        5.2578770
Fe        0        28.2181520        -9.0446690        8.3139950
Fe        0        28.2976510        -8.9615810        11.3392870
Mo        0        26.7500140        -12.1470690        5.1744640
O        0        26.9151420        -13.6248880        1.3564280
H        0        26.7374990        -14.5531170        1.5593240
H        0        26.5952440        -13.1589810        2.1707040
O        0        30.1023300        -14.3760080        6.1464280
H        0        29.5935450        -13.6092210        5.7877820
H        0        29.3935020        -14.9193250        6.5311380
O        0        27.4085010        -15.4521930        6.5274910
H        0        26.8908930        -14.6161120        6.6323930
H        0        27.0724550        -16.0539290        7.2034980
O        0        29.5956800        -13.4328910        1.6488190
H        0        30.0427800        -13.3665920        0.8006040
H        0        28.6338020        -13.4679780        1.4484660
O        0        30.1932680        -15.5116860        3.4981960
H        0        30.2111510        -15.1256910        4.3939430
H        0        30.0064170        -14.7775880        2.8838980
O        0        26.0956430        -15.5030030        3.8202520
H        0        26.5964330        -15.6237890        4.6424320
H        0        25.9284900        -14.5426800        3.7929470

22.5496426        -8.8159460        2.9242268        0.1        0.0
21.9828525        -8.5515047        4.3176472        0.1        0.0
23.4290757        -9.0156207        4.1569181        0.1        0.0
23.6978307        -7.8494963        3.2078572        0.1        0.0
22.2516059        -7.3853802        3.3685776        0.1        0.0
22.5496414        -8.8159503        10.6216797        0.1        0.0
21.9830787        -7.5656589        9.9519198        0.1        0.0
23.4294048        -7.5803407        10.4425575        0.1        0.0
23.6976040        -8.8353372        9.6143995        0.1        0.0
22.2512769        -8.8206634        9.1237648        0.1        0.0
32.4694417        -13.3959429        2.9242248        0.1        0.0
31.9030888        -13.1321098        4.3179381        0.1        0.0
33.3495107        -13.5953622        4.1565035        0.1        0.0
33.6171940        -12.4288906        3.2075661        0.1        0.0
32.1707712        -11.9656385        3.3689921        0.1        0.0
32.4694464        -13.3959506        10.6216812        0.1        0.0
33.0794809        -12.0047013        10.4632759        0.1        0.0
31.7702320        -12.3494986        9.7563278        0.1        0.0
33.7500480        -13.2115010        9.8099970        0.1        0.0
32.4407978        -13.5562952        9.1030406        0.1        0.0
27.0493588        -4.9240761        1.4191257        0.2        0.0
25.7337274        -5.5183750        1.9178470        0.2        0.0
26.9358812        -6.4267527        1.6678951        0.2        0.0
27.7835930        -5.5606603        2.5974520        0.2        0.0
26.5814345        -4.6522892        2.8474065        0.2        0.0
27.0493584        -4.9240714        12.1267810        0.2        0.0
25.7335948        -5.2223070        11.4108051        0.2        0.0
26.9351140        -4.7106660        10.6187340        0.2        0.0
27.7837255        -5.8567338        11.1657197        0.2        0.0
26.5822022        -6.3683752        11.9577835        0.2        0.0
33.3412579        -4.3869524        6.2818008        0.2        0.0
32.4633970        -3.2291045        5.8111071        0.2        0.0
32.1790218        -2.8540309        7.2641239        0.2        0.0
33.6207454        -2.9679884        6.7729560        0.2        0.0
33.0568864        -4.0118747        7.7348122        0.2        0.0
27.5095489        -2.6340581        2.9242289        0.2        0.0
28.1276072        -4.0217833        3.0824379        0.2        0.0
26.8129568        -3.6867947        3.7840574        0.2        0.0
28.7875234        -2.8122261        3.7414378        0.2        0.0
27.4728710        -2.4772396        4.4430647        0.2        0.0
27.5095479        -2.6340628        10.6216842        0.2        0.0
28.1276465        -2.3689729        9.2503754        0.2        0.0
26.8129705        -3.1387087        9.3595776        0.2        0.0
28.7875100        -3.3603065        10.2067425        0.2        0.0
27.4728304        -4.1300481        10.3159410        0.2        0.0
27.5095413        -7.2140700        -0.0859698        0.2        0.0
26.9430504        -8.4644366        0.5837100        0.2        0.0
28.3894093        -8.4496022        0.0931734        0.2        0.0
28.6574321        -7.1946071        0.9213908        0.2        0.0
27.2110724        -7.2094338        1.4119244        0.2        0.0
27.5095424        -7.2140732        13.6318926        0.2        0.0
26.9429449        -7.4782497        12.2383436        0.2        0.0
28.3892561        -7.0145144        12.3993829        0.2        0.0
28.6575378        -8.1807904        13.3483921        0.2        0.0
27.2112256        -8.6445247        13.1873617        0.2        0.0
20.9381509        -5.2580033        5.2678632        0.2        0.0
21.7358286        -5.9934983        4.1928945        0.2        0.0
21.1493076        -6.7670998        5.3720154        0.2        0.0
22.4482148        -5.0539024        5.1637010        0.2        0.0
21.8616868        -5.8275025        6.3428259        0.2        0.0
20.9381508        -5.2580034        8.2780632        0.2        0.0
21.7358407        -5.9934820        7.2030925        0.2        0.0
21.1493065        -6.7671015        8.3821952        0.2        0.0
22.4482160        -5.0539009        8.1739213        0.2        0.0
21.8616744        -5.8275186        9.3530278        0.2        0.0
32.0092505        -11.1059479        1.4191258        0.2        0.0
30.6957111        -10.8177507        2.1432507        0.2        0.0
31.9074492        -11.3154428        2.9286122        0.2        0.0
32.7414086        -10.1632502        2.3720540        0.2        0.0
31.5296745        -9.6655586        1.5866850        0.2        0.0
32.0092543        -11.1059548        12.1267779        0.2        0.0
32.4383948        -9.7429891        11.5873675        0.2        0.0
30.9677388        -10.0121200        11.8996682        0.2        0.0
32.4693828        -10.9688756        10.6768502        0.2        0.0
30.9987236        -11.2380112        10.9891570        0.2        0.0
24.8559233        -8.3197020        1.4828009        0.4        0.0
24.6082216        -8.1175412        2.9618442        0.4        0.0
25.9659573        -8.5783963        2.4781436        0.4        0.0
26.0422801        -7.3830233        1.5534552        0.4        0.0
24.6845464        -6.9221612        2.0371547        0.4        0.0
34.6628680        -10.0003211        2.9878935        0.4        0.0
33.2396014        -10.4024469        3.3078929        0.4        0.0
34.3144656        -11.4600638        3.1817139        0.4        0.0
35.1474996        -10.7370284        4.2176068        0.4        0.0
34.0726363        -9.6794184        4.3437899        0.4        0.0
35.0225510        -10.3724542        9.1360622        0.4        0.0
33.7063153        -9.6855833        9.4285203        0.4        0.0
33.3645520        -10.7670267        10.4302615        0.4        0.0
34.6628683        -10.0003146        10.5580116        0.4        0.0
34.6807859        -11.4538978        10.1377952        0.4        0.0
29.7029683        -6.0297052        2.9878982        0.4        0.0
30.1876148        -9.6562624        1.0481410        0.4        0.0
29.1128609        -9.2198142        0.0764364        0.4        0.0
30.0653392        -9.3671550        12.7842948        0.4        0.0
28.7510205        -9.4340322        13.5312184        0.4        0.0
28.4019609        -10.7134252        12.8024401        0.4        0.0
23.8454702        -6.1469304        10.6504918        0.75         0.0
35.1734173        -6.0058554        6.7729600        0.75         0.0
33.9906050        -5.0309739        6.7729600        0.75         0.0
35.1734173        -10.5858654        6.7729600        0.75         0.0
33.9906050        -9.6109839        6.7729600        0.75         0.0
30.7134298        -3.8569304        4.4005282        0.75         0.0
30.3404390        -2.7381346        5.4486288        0.75         0.0
30.7134298        -3.8569304        9.1453918        0.75         0.0
30.3404390        -2.7381346        8.0972912        0.75         0.0
30.7134312        -8.4369273        1.3903183        0.75         0.0
30.3404357        -7.3181345        2.4384195        0.75         0.0
30.7134312        -8.4369273        12.1555917        0.75         0.0
30.3404357        -7.3181345        11.1074905        0.75         0.0

20-82 0
6-31G
****
1-19 103-120 0
6-311+G(2df,p)
****
Fe Mo 0
LanL2DZ
****

Fe Mo 0
LanL2DZ

surface=sas

xsc6

你最好输入文件输出文件上传一下，你这个基组感觉有点怪怪，opt也别用smd 然后就是你这电荷和自选多重度确定对么

guyifan

电荷和自选多重度倒是没有问题，因为添加水分子做显式溶剂模型没出现问题，但是加入iefpcm反而出现问题，smd的话我再去修改下，这就上传下输出文件。谢谢您的指导。

guyifan

收敛相也迟迟不收敛，输出文件20M,受到上传限制。

Freeman

救命。输入文件太逆天了。就挑几个问题说吧。
1
这个问题真的需要设置额外点电荷吗？大概不用吧。删掉charge关键词以及坐标和基组之间的一堆数。
2
电荷和自旋多重度太离谱了。问一下老师和师兄师姐，确认一下是不是这样。
3
溶剂化模型写错了。iefpcm不是溶剂，而是模型，和smd是一个等级的概念。要么写smd，要么写iefpcm，不能两个都写。高斯不认识iefpcm这个溶剂，就当成默认的水了。倒是歪打正着。
4
输出文件写的cavity type仍然是范德华，说明高斯没有正确读出surface=sas这句话。可能这句话要写在基组前面。

guyifan

Freeman 发表于 2024-4-23 14:23
救命。输入文件太逆天了。就挑几个问题说吧。
1
这个问题真的需要设置额外点电荷吗？大概不用吧。删掉cha ...

好的好的，太谢谢啦