高斯不收斂問題

1、首先，我們必須理解收斂是什么意思。在自洽場(chǎng)（SCF）計(jì)算中，自洽循環(huán)中，首先產(chǎn)生一個(gè)軌道占據(jù)的初始猜測(cè)，1）然后根據(jù)此軌道占據(jù)構(gòu)造電荷密度和哈密頓量。2）對(duì)角化哈密頓量，得到新的軌道能級(jí)和占據(jù)。3）產(chǎn)生新的電荷分布和哈密頓量，重復(fù)步驟2）經(jīng)過一定次數(shù)的循環(huán)后，某次循環(huán)前和循環(huán)后的電荷密度差別小于一定的標(biāo)準(zhǔn)，我們稱之為收斂。如果以上過程不能收斂，貝Ugaussian給出convergence failure的警告。如果SCF計(jì)算收斂失敗，你首先會(huì)采取哪些技巧呢？這里是我們強(qiáng)烈推薦的首選方法。 1考慮使用更小的基組 由于一定的基組對(duì)應(yīng)于一定精度和速度，所以更換基組并不在所有的情況下都適用。方法是首

2、先用小基組進(jìn)行計(jì)算，由前一個(gè)波函得到用于大基組計(jì)算的初始猜測(cè)（Guess=Read自動(dòng)進(jìn)行）。2增加最大循環(huán)步數(shù) Gaussian默認(rèn)的最大循環(huán)步數(shù)為64 （SCF=DM或SCF=QC方法則為512）,如果循環(huán)次數(shù)超過這個(gè)數(shù)目則會(huì)匯報(bào)convergence failure o在一定的情況下，不收斂的原因僅僅是因?yàn)樽畲笱h(huán)步數(shù)不夠?？梢酝ㄟ^設(shè)置 maxcyc來增大最大循環(huán)步數(shù)。更多的SCF迭代 （SCF（MaxCycle=N），其中N是迭代數(shù)）。這很少有幫助，但值得一試。3 放寬收斂標(biāo)準(zhǔn) 如果接近SCF但未達(dá)到，收斂標(biāo)準(zhǔn)就會(huì)放松或者忽略收斂標(biāo)準(zhǔn)。這通常用于不是在初始猜 測(cè)而是在平衡結(jié)構(gòu)收斂的幾何

3、優(yōu)化。SCF=Sleazy放松收斂標(biāo)準(zhǔn)，Conver選項(xiàng)給出更多的控制。在Gaussian中設(shè)置收斂標(biāo)準(zhǔn)的命令行是Conv=8表示10-8,你可以將這個(gè)標(biāo)準(zhǔn)降低，比如改為10-5等等。微調(diào)結(jié)構(gòu)是比較有效的方法，能量收斂標(biāo)準(zhǔn)不宜太低，最好不低于10-7.4 嘗試改變初始構(gòu)型首先略微減小鍵長(zhǎng)，接下來略微增加鍵長(zhǎng)，接下來再對(duì)結(jié)構(gòu)作一點(diǎn)改變。 5 嘗試能級(jí)移動(dòng) Level shifting （SCF=Vshift） 如果不收斂的原因是波函數(shù)的震蕩行為，通常是因?yàn)樵谙嘟哪芰可系奶┑幕旌稀?duì)于這種情況，我們可以采用level shifting的方法。Level shifting的含義是人工的升高非占據(jù)軌

4、道的 能級(jí)，以防止和最高占據(jù)軌道之間的混合，以達(dá)到收斂的目的。在 Gaussian中此方法的關(guān) 鍵詞為 SCF=Vshift 6 使用強(qiáng)制的收斂方法 SCF=QC SCF=QC通常最佳，但在極少數(shù)情況下SCF=DM更快。此關(guān)鍵字將大大增加計(jì)算時(shí)間，但是收斂的機(jī)會(huì)更大。不要忘記給計(jì)算額外增加一千個(gè)左右的迭代。應(yīng)當(dāng)測(cè)試這個(gè)方法獲得的波函，保證它最小，并且正好不是穩(wěn)定點(diǎn)（使用Stable關(guān)鍵字）。7對(duì)開殼層體系，嘗試收斂到同一分子的閉殼層離子，接下來用作開殼層計(jì)算的初始猜測(cè)。添加電子可以給出更合理的虛軌道，但是作為普遍的經(jīng)驗(yàn)規(guī)則，陽(yáng)離子比陰離子更容易收斂。 選項(xiàng)Guess=Read定義初始猜測(cè)從 G

5、aussian計(jì)算生成的checkpoint文件中讀取。8 一些程序通過減小積分精度加速SCF。對(duì)于使用彌散函數(shù)，長(zhǎng)程作用或者低能量激發(fā)態(tài)的體系，必須使用高積分精度：SCF=NoVarAcc。9改變模型或方法可以考慮改變模型方法。比較常見的方法有 HF,GVB,MCSCF,CASSCF,MPn等。改變模型方 法通常也會(huì)收斂性質(zhì)。通常，精度更高的方法更難收斂。精度比較低的方法產(chǎn)生的計(jì)算結(jié)果 可以作為高精度計(jì)算的初始猜測(cè)?？紤]使用不同理論級(jí)別的計(jì)算。這并不總是實(shí)用的，但除此之外，增加迭代數(shù)量總是使得計(jì)算時(shí)間和使用更高理論級(jí)別差不多。10關(guān)閉DIIS外推(SCF=NoDIIS)。同時(shí)進(jìn)行更多的迭代

6、(SCF=(MaxCycle=N)。 11試著改用 DIIS之外其它方法(SCF=SD或SCF=SSD)。Gaussian不收斂的可能原因及對(duì)策：1由于體系有很多能量相近的能級(jí)，導(dǎo)致計(jì)算不收斂。如果計(jì)算中采用的是 Hartree-Fock方法或者其他的混合形式的交換相關(guān)勢(shì)(如 B3LYP)，則 可以嘗試一下的方法進(jìn)行改進(jìn)。scf=(noincfock,conver=11,maxcyc=1025) iop(5/22=20)用以上的命令行將強(qiáng)制Gaussian采用EDIIS的算法。此算法計(jì)算量更大，但是更加穩(wěn)定。在命令行里也增加了最大循環(huán)數(shù)，以增加收斂的到基態(tài)的可能性。2. 檢查是否有初始文件錯(cuò)誤

7、常見初級(jí)錯(cuò)誤：a. 自旋多重度錯(cuò)誤b. 變量賦值為整數(shù)c. 變量沒有賦值或多重賦值d. 鍵角小于等于 0度，大于等于180度e. 分子描述后面沒有空行f. 二面角判斷錯(cuò)誤，造成兩個(gè)原子距離過近g. 分子描述一行內(nèi)兩次參考同一原子，或參考原子共線3. SCF(自洽場(chǎng))不收斂則一般是L502錯(cuò)誤，省卻情況做 64個(gè)cycle迭代(G03缺省128 cycles)a. 修改坐標(biāo)，使之合理b. 改變初始猜Guess=Huckel或其他的，看 Guess關(guān)鍵詞。c. 增加疊代次數(shù) SCFCYC=N (對(duì)小分子作計(jì)算時(shí)最好不要增加，很可能結(jié)構(gòu)不合理)d. iop (5/13=1)這樣忽略不收斂，繼續(xù)往下做

8、。4. 分子對(duì)稱性改變a. 修改坐標(biāo)，強(qiáng)制高對(duì)稱性或放松對(duì)稱性b. 給出精確的、對(duì)稱性確定的角度和二面角。如 CH4的角度給到109.47122c. 放松對(duì)稱性判據(jù)Symm=loosed. 不做對(duì)稱性檢查iop(2/16=1)(最好加這個(gè)選項(xiàng))iop(2/16=2)則保持新的對(duì)稱性來計(jì)算5. Opt時(shí)收斂的問題a. 修改坐標(biāo)，使之合理b. 增加疊代次數(shù) optcyc=N6. 優(yōu)化過渡態(tài)，若勢(shì)能面太平緩，則不好找到。iop(1/8=10)默認(rèn)30 (下一個(gè)結(jié)構(gòu)和該結(jié)構(gòu)的差別0.3?)，可改成10。如果每一步都要用到小的步長(zhǎng)，應(yīng)該加 opt(notrustupdate)7. 在CI (組態(tài))方法中

9、如QCISD(T),CCSD(T),CID 方法中，省卻最大循環(huán)50,若出錯(cuò)(L913錯(cuò)誤)解決方法：#P QCISD(maxcyc=N) 注：N < 5128. 優(yōu)化過渡態(tài)opt=TS (給出過渡態(tài))opt=qst2 (給出反應(yīng)物和產(chǎn)物)opt=qst3 (給出反應(yīng)物和產(chǎn)物和過渡態(tài) )a.用 G03 時(shí)的出錯(cuò) opt=ts 必須加 FC (force constant)寫法：opt=(TS, calcFc)or opt=(TS,calchffc)計(jì)算 HF 力常數(shù)，對(duì) QCISD , CCSD 等方法用； or opt=(TS,modRedundant)(最好寫這個(gè))b. 如果計(jì)算采用

10、 QCISD計(jì)算(不好計(jì)算 FC)則寫為 QCISD opt=(TS, calcHFFC)(用 HF 計(jì)算 FC)9. 無法寫大的 Scratch文件RWFa. 劈裂 RWF 文件 rwf=loc1,size1,loc2,size2, .,locN,-1b. 改變計(jì)算方法 MP2=Direct可以少占硬盤空間c. 限制最大硬盤 maxdisk=N GB , *MB,有些系統(tǒng)寫 2GB會(huì)出錯(cuò)，可以寫 2000MB10. FOPT出錯(cuò)原因是變量數(shù)與分子自由度數(shù)不相等?？捎?POPT或直接用OPT11. 優(yōu)化過渡態(tài)只能做一個(gè)STEP原因是負(fù)本征數(shù)目不對(duì)添加iop (1/11) =1或者noeigen

11、test(eigentest是表示優(yōu)化過渡態(tài)檢測(cè)分子振動(dòng)的本征值，過渡態(tài)只有一個(gè)負(fù)值，但優(yōu)化的時(shí)候往往出現(xiàn)兩個(gè)或者更多的情況，默認(rèn)的是如果出現(xiàn)多于一個(gè)的情況就停止優(yōu)化計(jì)算，這往往時(shí)不必要的，noeigentest表示優(yōu)化時(shí)不坐此檢測(cè)，實(shí)際上優(yōu)化過渡態(tài)的時(shí)候出現(xiàn)兩個(gè)三個(gè)很 正常的，只是第一個(gè)一般負(fù)的很大，接著的比較小，不會(huì)影響過渡態(tài)的搜尋，當(dāng)然有時(shí)候需要自己判斷是不是接近真實(shí)的過渡態(tài)。其實(shí)經(jīng)驗(yàn)上也是過渡態(tài)的本征虛頻越小，相對(duì)越難 找。)對(duì)于(L502, L508, L9999)出錯(cuò)的對(duì)策對(duì)于一個(gè)優(yōu)化計(jì)算，它的過程是先做一個(gè)SCF計(jì)算，得到這個(gè)構(gòu)型下的能量，然后優(yōu)化構(gòu)型，再做SCF,然后再優(yōu)化構(gòu)型

12、。因此，會(huì)有兩種不收斂的情況：一是在某一步的SCF不收斂(L502錯(cuò)誤)，或者構(gòu)型優(yōu)化沒有找到最后結(jié)果(L9999錯(cuò)誤)。預(yù)備知識(shí)：計(jì)算時(shí)保存chk文件，可以在后續(xù)計(jì)算中使用guess=read讀初始猜測(cè).對(duì)于SCF不收斂，通常有以下的解決方法：1. 使用小基組，或低級(jí)算法計(jì)算，得到 scf收斂的波函數(shù)，用 guess=read讀初始波函數(shù)。2. 使用scf=qc,這個(gè)計(jì)算會(huì)慢，而且需要用stable關(guān)鍵字來測(cè)試結(jié)果是否波函數(shù)穩(wěn)定。如果這個(gè)還不收斂，會(huì)提示 L508錯(cuò)誤。3. 改變鍵長(zhǎng)，一般是縮小一點(diǎn)，有時(shí)會(huì)有用。4. 計(jì)算相同體系的其他電子態(tài)，比如相應(yīng)的陰離子、陽(yáng)離子體系或單重態(tài)體系，得到的

13、收 斂波函數(shù)作為初始猜測(cè)進(jìn)行計(jì)算。對(duì)于優(yōu)化不收斂，即L9999錯(cuò)誤，實(shí)際上是在規(guī)定的步數(shù)內(nèi)沒有完成優(yōu)化，即還沒有找到極 小值點(diǎn)。(或者對(duì)于過渡態(tài)優(yōu)化，還沒有找到過渡態(tài)) 這有幾種可能性：1. 看一下能量的收斂的情況，可能正在單調(diào)減小，眼看有收斂的趨勢(shì)，這樣的情況下，只 要加大循環(huán)的步數(shù)(opt(maxcycle=200),可能就可以解決問題了。2. 加大循環(huán)步數(shù)還不能解決的(循環(huán)步數(shù)有人說超過200再不收斂，再加也不會(huì)有用了，這雖然不一定絕對(duì)正確，但 200步應(yīng)該也差不多了)，有兩種可能。一是查看能量，發(fā)現(xiàn)能 量在振蕩了，且變化已經(jīng)很小了，這時(shí)可能重新算一下，或者構(gòu)型稍微變一下，繼續(xù)優(yōu)化，就可

14、以得到收斂的結(jié)果(當(dāng)然也有麻煩的，看運(yùn)氣和經(jīng)驗(yàn)了)；二是構(gòu)型變化太大，和你預(yù)計(jì)的差別過大，這很可能是你的初始構(gòu)型太差了，優(yōu)化不知道到哪里去了，這時(shí)最好檢查一下初始構(gòu)型，再?gòu)念^優(yōu)化。3. 對(duì)于L9999快達(dá)到收斂時(shí)，考慮減小優(yōu)化步長(zhǎng)有時(shí)對(duì)于能量振蕩的情況也是有用的， opt(maxstep=1).(flyingheart )Hello,here is the summary of all the answers I got for my problem. Thanks (in alphabetically order) to Laurence Cuffe, Kerwin Dobbs, Paul

15、Kiprof, Liviu Mirica, Cory Pye and Dave Shobe.The problem was solved by increasing the number of SCF cycles, as suggested by many. This seems very often necessary with transition metals. The two other suggestions (using first a small basis set before switching to the larger, or turning off "vsh

16、ift", did not work, unfortunately. Thanks again, R.Herrmann.The original question:Hello,I recently tried to calculate the NMR chemical shifts of VOCl3 with Gaussian 98:#B3LYP/6-31G* SP NMRVOCl30 1clcl1 clcl2cl2 clcl31 clclcl3o3 ocl42 oclcl41 dih4v3 vcl52 vclcl51 dih5clcl33.506243clclcl359.999oc

17、l42.938735oclcl453.376dih4-64.586vcl52.118170vclcl534.144dih5-31.640The result (even with several variations in the input) is always an error termination when trying the "incremental Fock formation" (convergence criterion not met). The structure was optimized with GAMESS.Any suggestion wha

18、t I could do is wellcome - thank you!108 basis functions41 alpha electrons nuclear repulsion energy288 primitive gaussians41 beta electrons561.9817715949 Hartrees.One-electron integrals computed using PRISM.NBasis= 108 RedAO= T NBF= 108NBsUse= 108 1.00D-04 NBFU= 108Projected Huckel Guess.Requested c

19、onvergence on RMS density matrix=1.00D-08 within 64 cycles.Requested convergence on MAX density matrix=1.00D-06.Virtual orbitals will be shifted by 0.200 hartree.Restarting incremental Fock formation.Restarting incremental Fock formation.Restarting incremental Fock formation.>>>>>>

20、>>>> Convergence criterion not met.SCF Done: E(RB+HF-LYP) =-2359.59549049 A.U. after 65 cyclesConvg =0.5851D-03-V/T =1.9657S*2=0.0000Convergence failure - run terminated.Error termination via Lnk1e in /usr/local/g98_a7/l502.exe.Job cpu time: 0 days 0 hours 16 minutes 33.0 seconds.File le

21、ngths (MBytes):RWF= 14 Int= 0 D2E=0 Chk= 1 Scr=The answers:1. From: Kerwin D Dobbs <Kerwin.D.Dobbs-: at :-USA> To: "Rudolf Herrmann <herrmann" <herrmann-: at :-itqb.unl.pt> Subject: Re: CCL:G98: VOCl3 problem: convergence criterion not metDr. Herrmann,Try turning off the &qu

22、ot;vshift" by changing the line,#B3LYP/6-31G* SP NMRto# B3LYP/6-31G* SP NMR IOP(5/10=-1)When there is a TM in the molecule, G98 automatically turns on "VSHIFT" which usually makes SCF convergence worse. The above addition to your line turns off "VSHIFT". Give it a try and se

23、e what happens.Kerwin2. From: Cory Pye<cpye-: at :-crux.stmarys.ca>To: Rudolf Herrmann <herrmann-: at :-itqb.unl.pt>Subject: Re: CCL:G98: VOCl3 problem: convergence criterion not met Hello,your SCF is not converging. Try running, at first, HF/STO-3G. If it works, thenThen try B3L YP/6-31

24、G* afterwards with guess=read in the inp ut. It usually works like a charm for me, I always use something like this.-Cory3. From: LiviuMihailMirica<liviu-: at :-Stanford.EDU>To: Rudolf Herrmann <herrmann-: at :-itqb.unl.pt>Subject: Re: CCL:G98: VOCl3 problem: convergence criterion not me

25、tTry adding SCFcyc=1000 in your command lines. This will increase the # of SCF cycles up to the desired number (1000 in this case). that should be enough for your system.4. From: Laurence Cuffe<Laurence.Cuffe-: at :-ucd.ie>To: Rudolf Herrmann <herrmann-: at :-itqb.unl.pt>Subject: Re: CCL

26、:G98: VOCl3 problem: convergence criterion not met And it should run these keywords change the default SCF method and then increase the number of SCF cycles allowed. I got the job to complete in about 50 minutes running it in core on a linux box with a memory limit of 600MB but it this doesn't m

27、ake sense to you don't worry about it. Try using the modified line above and if there are any problems don't hesitate to come back to me.All the bestLaurence Cuffe5. From: Paul Kiprof<pkiprof-: at :->To: herrmann-: at :-itqb.unl.ptSubject: VOCL3Dear Dr. Herrmann,I would add to

28、 the inputSCF(direct,qc,maxcycle=1000)This will increase the number of convergence cycles and hopefully converge the calculation. I always do that when I have transition metal compounds.Best wishes,Paul Kiprof6. From: "Shobe, Dave" <dshobe-: at :-sud->To: 'Rudolf Herrmann' &l

29、t;herrmann-: at :-itqb.unl.pt> Subject: RE: G98: VOCl3 problem: convergence criterion not met The standard 64 SCF "cycles" (or iterations of the variational procedure) is rarely enough for transition metal compounds, even closed shell ones like VOCl3. Read the section in the Gaussian ma

30、nual describing the "scf" keyword, and you will find various options which may improve the convergence. A lot of it is trial and error, with these different scf options. The simplest is scf(maxcyc=999) which increases the number of cycles/iterations to 999 before the program "gives up

31、" and prints the error message. Scf=qc and scf(vshift=#) have also been useful in my experience.More advice: once you get a good wavefunction, make a backup of the checkpoint file, so if the wave function becomes corrupted in subsequent calculations you can always go back to it.Dr. Rudolf Herrm

32、annInstituto de TecnologiaQuimica e Biologica (ITQB)Av. da Republica, EAN, Apt. 127P-2781-901 Oeirase-mail: herrmann-: at :-itqb.unl.pt gaussian學(xué)習(xí)筆記(一)默認(rèn)分類 2008-11-17 17:41:53閱讀124評(píng)論0 字號(hào)：大中小訂閱用Gaussian計(jì)算不收斂的問題問題一：Problem detected with inexpensive integrals.Switching to full accuracy and repeating la

33、st cycle.Cycle 2 Pass 1 IDiag 1:Warning! Spurious integrated density:NE= 188 NElCor= 0 Integral= -76.06821 Tolerance=1.00D-03Consistency failure #2 in CalDSu.Error termination via Lnk1e in /home/gaussian/g98/l502.exe.Job cpu time: 1 days 3 hours 21 minutes 30.2 seconds.File lengths (MBytes): RWF= 11

34、8 Int= 0 D2E= 0 Chk= 17 Scr= 1 解決：(from ccl)As far as I know it's a dft-related error. Gaussian checks the charge resulting from the integration of the density to the number of electrons, and in case of a large difference prints this eror message. Try a denser grid with the "integral"

35、keyword.可是我采用integral (Grid=ultrafine )網(wǎng)格仍是出現(xiàn)以下錯(cuò)誤信息： Raffenetti 2 integral format.Two-electron integral symmetry is turned off.306 basis functions 778 primitive gaussians51 alpha electrons 51 beta electronsnuclear repulsion energy 865.0704012775 Hartrees.One-electron integrals computed using PRISM.1

36、 Symmetry operations used in ECPInt.ECPInt: NShTT= 17391 NPrTT= 123061 LenC2= 14761 LenP2D= 57336.LDataN: DoStor=F MaxTD1= 3 Len= 28LDataN: DoStor=T MaxTD1= 3 Len= 28NBasis= 306 RedAO= T NBF= 306NBsUse= 306 1.00D-04 NBFU= 306Projected Huckel Guess.of initial guess= 0.0000Requested convergence on RMS

37、 density matrix=1.00D-06 within 800 cycles.Requested convergence on MAX density matrix=1.00D-04.Integral accuracy reduced to 1.0D-05 until final iterations.Warning! Spurious integrated density:NE= 102 NElCor= 0 Integral= 102.20958 Tolerance=1.00D-03 Consistency failure #2 in CalDSu.Error termination

38、 via Lnk1e in ZusrZpeopleZgaussianZg98Zl502.exe.Job cpu time: 0 days 0 hours 13 minutes 24.7 seconds.File lengths (MBytes): RWF= 90 Int= 0 D2E= 0 Chk= 2 Scr= 1請(qǐng)問各位大俠有沒有什么解決方法？問題二：>>>>>>>>>> Convergence criterion not met.SCF Done: E(UB-PW91) = -1306.37542885 A.U. after 8

39、01 cyclesConvg = 0.1221D+01 -VZT = 2.0190S*2 = 0.0000KE= 1.281978320748D+03 PE=-8.302694580633D+03 EE= 3.210822249139D+03Annihilation of the first spin contaminant:S*2 before annihilation 0.0000, after 0.0000Convergence failure - run terminated.Error termination via Lnkle in ZusrZpeQpleZgaussianZg98

40、Zl502.exe.Job cpu time: 2 days 2 hours 6 minutes 40.2 seconds.File lengths (MBytes): RWF= 97 Int= 0 D2E= 0 Chk= 19 Scr= 1一般而言，解決收斂的問題有以下幾種，第一，波函數(shù)初猜，不穩(wěn)定的波函數(shù)，積分難以收斂是肯定的，解決方案是作stable=opt，做完以后，在此基礎(chǔ)上從check文件中讀入波函數(shù)。linkl%chk#p hf opt guess = read geom= check第二，過渡金屬原子會(huì)出現(xiàn) d軌道分裂問題，如果分裂能不夠大，電子顯然會(huì)在高能軌道 和低能軌道之間

41、跳動(dòng)，能量會(huì)上下波動(dòng)，難以收斂，解決方案是加入vshift關(guān)鍵詞第三，氫鍵，過渡金屬，默認(rèn)的129cycle所難以達(dá)到收斂，所以加大 scf(cyc=500)可能可以解決問題。第四，初始結(jié)構(gòu)一直容易被人所忽略，所以很多方法沒有解決問題的時(shí)候，還是希望能 仔細(xì)檢查一遍結(jié)構(gòu)第五，scf的時(shí)候，檢查一下力常數(shù)也是好的，通常的解決方案作一下calcfc對(duì)于第二個(gè)收斂的問題，我正在用SCF(vshift,QC)嘗試，得過幾天才知道行不行。至少我已嘗試了將內(nèi)坐標(biāo)改為直角坐標(biāo)；加大循環(huán)次數(shù)(800仍不收斂)從check文件中讀入，這些都不行。對(duì)于第一個(gè)問題，我是用massage考慮分子間相互作用，可是好多都

42、出現(xiàn)第一個(gè)問題，根本不讓你算，用了 integral(grid=ultrafine)還是不行，輸入文件如下：(這是一個(gè)不含金屬原子 的，有好幾個(gè)含金屬的輸入，也是出現(xiàn)這個(gè)情況。芝麻582：1.在Guess關(guān)鍵字中使用 Core, Huckel或Mix選項(xiàng)，試驗(yàn)不同的初始猜測(cè)。2. 對(duì)開殼層體系，嘗試收斂到同一分子的閉殼層離子，接下來用作開殼層計(jì)算的初始猜測(cè)。添加電子可以給出更合理的虛軌道，但是作為普遍的經(jīng)驗(yàn)規(guī)則，陽(yáng)離子比陰離子更容易收斂。選項(xiàng)Guess=Read定義初始猜測(cè)從 Gaussian計(jì)算生成的checkpoint文件中讀取。3. 另一個(gè)初始猜測(cè)方法是首先用小基組進(jìn)行計(jì)算，由前一個(gè)波函得到用于大基組計(jì)算的初始猜測(cè) (Guess=Read自動(dòng)進(jìn)行)。4. 嘗試能級(jí)移