高斯不收斂，如何優(yōu)化
#p b3lyp/lanl2dz opt pop=(full,npa)
>>>>>>>>>> Convergence criterion not met.
SCF Done: E(UB+HF-LYP) = -1714.45115901 A.U. after 129 cycles
Convg = 0.8068D-05 -V/T = 2.0529
S**2 = 0.9172
KE= 1.628346096785D+03 PE=-1.275059909267D+04 EE= 5.005782636124D+03
Annihilation of the first spin contaminant:
S**2 before annihilation 0.9172, after 0.7568
Convergence failure -- run terminated.
Error termination via Lnk1e in d:program filesg03l502.exe at Fri Jun 18 18:59:12 2021.
Job cpu time: 1 days 5 hours 55 minutes 14.0 seconds.
File lengths (MBytes): RWF= 206 Int= 0 D2E= 0 Chk= 18 Scr= 1
請(qǐng)問高手，我這個(gè)從報(bào)錯(cuò)看應(yīng)該是沒有收斂，但是應(yīng)該如何優(yōu)化或者修改？救救孩子吧。。。。@jiewei
返回小木蟲查看更多
京公網(wǎng)安備 11010802022153號(hào)
看起來是高斯的Scratch文件夾里內(nèi)容被占用了。你修改好之后保存成另外一個(gè)名字的輸入文件。并且把高斯的scratch文件夾內(nèi)容清空 再計(jì)算就可以了。
包含重金屬，也不一定非要所有原子都用LanL2DZ 可以用混合基組 給重原子用LanL2DZ(這個(gè)2021年的今天已經(jīng)很不推薦了精度太差，建議SDD） 給其他輕原子（前三周期）用普通的基組比如 6-31G(D)之類的。會(huì)合理很多。
之前那個(gè)任務(wù)直接改不太行的。要改成混合基組需要特別的寫法。見我這個(gè)教程1.1節(jié) http://m.gaoyang168.com/t-14242735-1
如果基組設(shè)置沒有問題。電荷和多重度也設(shè)置的沒有什么問題。結(jié)構(gòu)本身原子間的距離并沒有太大異常而僅僅是gaussview不顯示某些鍵那一般不會(huì)有什么問題。
見我這個(gè)教程3.3節(jié) http://m.gaoyang168.com/t-14242735-1
3.3 化學(xué)鍵的顯示問題。
諸如 為什么優(yōu)化完成后“缺了什么什么鍵“，“哪個(gè)鍵變成了單/雙/三鍵”，“那個(gè)鍵不見了”等問題。
這類問題的回答是：沒有關(guān)系。gaussview怎么顯示的，都沒什么所謂。鍵顯示的符合預(yù)期也好或者不符合預(yù)期也罷都不能說明計(jì)算結(jié)果本身是對(duì)還是錯(cuò)。
用高斯做的量化計(jì)算需要注意的是量化計(jì)算中，或者說薛定諤方程中根本沒有“鍵”這種東西。一切都是波函數(shù)。也就是一個(gè)一個(gè)的分子軌道。優(yōu)化出來的就是原子間的相對(duì)位置，因?yàn)檠Χㄖ@方程需要的只是原子的位置信息。
而gaussview在顯示分子的時(shí)候會(huì)根據(jù)一個(gè)內(nèi)在的鍵長(zhǎng)標(biāo)準(zhǔn)（gaussview內(nèi)置了一套原子的半徑信息，會(huì)拿原子間的距離相應(yīng)的兩個(gè)原子的半徑之和作比較，比如鍵長(zhǎng)超過了原子半徑之和的百分之多少就屬于單鍵，之類之類的。），完全是根據(jù)原子間的距離來顯示鍵的種類的。跟電子真正在分子中的排布一毛錢關(guān)系都沒有。所以要分析原子間的成鍵需要進(jìn)行波函數(shù)分析。算算鍵級(jí)之類的去判斷。gaussview顯示的完全就是為了好看，但是沒有任何化學(xué)意義的。并且也是可以自己手動(dòng)更改鍵的種類的。
，