高斯收斂問題

高斯收斂問題首先，我們必須理解收斂是什么意思。在自洽場(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)，我們稱之為收斂。如果以上過程不能收斂，則gaussian給出convergencefailure的警告。如果SCF計(jì)算收斂失敗，你首先會(huì)采取哪些技巧呢？這里是我們強(qiáng)烈推薦的首選方法。1考慮使用更小的基組由于一定的基組對(duì)應(yīng)于一定精度和速度，所以更換基組并不在所有的情況下都適用。方法是首先用小基組進(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)convergencefailure。在一定的情況下，不收斂的原因僅僅是因?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)收斂的幾何優(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)Levelshifting(SCF=Vshift)如果不收斂的原因是波函數(shù)的震蕩行為，通常是因?yàn)樵谙嘟哪芰可系奶┑幕旌?。?duì)于這種情況，我們可以采用levelshifting的方法。Levelshifting的含義是人工的升高非占據(jù)軌道的能級(jí)，以防止和最高占據(jù)軌道之間的混合，以達(dá)到收斂的目的。在Gaussian中此方法的關(guān)鍵詞為SCF=Vshift6使用強(qiáng)制的收斂方法SCF=QCSCF=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ī)則，陽離子比陰離子更容易收斂。選項(xiàng)Guess=Read定義初始猜測(cè)從Gaussian計(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)行更多的迭代(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)的可能性。檢查是否有初始文件錯(cuò)誤常見初級(jí)錯(cuò)誤：自旋多重度錯(cuò)誤變量賦值為整數(shù)變量沒有賦值或多重賦值鍵角小于等于0度，大于等于180度分子描述后面沒有空行二面角判斷錯(cuò)誤，造成兩個(gè)原子距離過近分子描述一行內(nèi)兩次參考同一原子，或參考原子共線3.SCF（自洽場(chǎng)）不收斂則一般是L502錯(cuò)誤，省卻情況做64個(gè)cycle迭代（G03缺省128cycles）修改坐標(biāo)，使之合理改變初始猜Guess=Huckel或其他的，看Guess關(guān)鍵詞。增加疊代次數(shù)SCFCYC=N（對(duì)小分子作計(jì)算時(shí)最好不要增加，很可能結(jié)構(gòu)不合理）d.iop（5/13=1）這樣忽略不收斂，繼續(xù)往下做。分子對(duì)稱性改變修改坐標(biāo)，強(qiáng)制高對(duì)稱性或放松對(duì)稱性給出精確的、對(duì)稱性確定的角度和二面角。如CH4的角度給到109.47122放松對(duì)稱性判據(jù)Symm=loose不做對(duì)稱性檢查iop（2/16=1）（最好加這個(gè)選項(xiàng)）iop（2/16=2）則保持新的對(duì)稱性來計(jì)算Opt時(shí)收斂的問題修改坐標(biāo)，使之合理增加疊代次數(shù)optcyc=N優(yōu)化過渡態(tài)，若勢(shì)能面太平緩，則不好找到。到草窣IdOd由也。髭苗玄暮囹甲目會(huì)好與暮香至耆苗誓果吊IdOd-0T8IAI000乙宜召也’果吊號(hào)此乙宜場(chǎng)爰扁身0內(nèi)￥￥￥￥'此N=>|S!Pxe山器副￥智閽遂I回PaJ!Q=^dIAI累旦竟艮至旎-q「'l\PO|' '*z!S'Qor”zSno|TM」％朝W3MX譴籍3MX朝WM^bjds陽￥宜累堊,6（DdM+ldH困）（DddH3|B3，Sl）=ldoasiDd區(qū)宜［M（DdM44±5±）M44asod由瘤44鑿密qg宜卻曾（iuepunpaypouj/si）=ldojo-由累旦會(huì)dSZD'asodK'暮景qr卅竟44（祀ip|e】s_L）Fd。」。（3^3163^1）=^0:累宜（lueisuoDaDJOpDd叩蟲0s=id。果吊陽明￡09困宅（翠跳:K也麟旻也端四四吊核）gsb=do（麟旻也昭四四吊核）^sb=ido（翠跳國印核）Sl=ldo翠跳EPW?8^IS>N:王；（N=^3xeiu）asiDdd#:累旦散麹（淳果沮61）果吊是’os也甄￥窖。￥易’aiDXBasDDXBasiDd□$中累旦（翠應(yīng)）D丑W（aiepdnisruiou）ido叩匿〔回'習(xí)爭(zhēng)陽\「/［?由查腥等一母鑿□$°0■［御旎也’（V&0暗鬢陽歐旱裕〔也歐弱￡一"）0￡*￥獸（OT=8/T）doi用OPT11.優(yōu)化過渡態(tài)只能做一個(gè)STEP原因是負(fù)本征數(shù)目不對(duì)添加iop（1/11）=1或者noeigentest（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)型。。。因此，會(huì)有兩種不收斂的情況：一是在某一步的SCF不收斂（L502錯(cuò)誤），或者構(gòu)型優(yōu)化沒有找到最后結(jié)果（L9999錯(cuò)誤）。預(yù)備知識(shí):計(jì)算時(shí)保存chk文件，可以在后續(xù)計(jì)算中使用guess=read讀初始猜測(cè).對(duì)于SCF不收斂，通常有以下的解決方法：使用小基組，或低級(jí)算法計(jì)算，得到scf收斂的波函數(shù)，用guess=read讀初始波函數(shù)。使用scf=qc，這個(gè)計(jì)算會(huì)慢，而且需要用stable關(guān)鍵字來測(cè)試結(jié)果是否波函數(shù)穩(wěn)定。如果這個(gè)還不收斂，會(huì)提示L508錯(cuò)誤。改變鍵長(zhǎng)，一般是縮小一點(diǎn)，有時(shí)會(huì)有用。計(jì)算相同體系的其他電子態(tài)，比如相應(yīng)的陰離子、陽離子體系或單重態(tài)體系，得到的收斂波函數(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))這有幾種可能性：看一下能量的收斂的情況，可能正在單調(diào)減小，眼看有收斂的趨勢(shì)，這樣的情況下，只要加大循環(huán)的步數(shù)(opt(maxcycle=200)),可能就可以解決問題了。加大循環(huán)步數(shù)還不能解決的(循環(huán)步數(shù)有人說超過200再不收斂，再加也不會(huì)有用了，這雖然不一定絕對(duì)正確，但200步應(yīng)該也差不多了)，有兩種可能。一是查看能量，發(fā)現(xiàn)能量在振蕩了，且變化已經(jīng)很小了，這時(shí)可能重新算一下，或者構(gòu)型稍微變一下，繼續(xù)優(yōu)化，就可以得到收斂的結(jié)果(當(dāng)然也有麻煩的，看運(yùn)氣和經(jīng)驗(yàn)了)；二是構(gòu)型變化太大，和你預(yù)計(jì)的差別過大，這很可能是你的初始構(gòu)型太差了，優(yōu)化不知道到哪里去了，這時(shí)最好檢查一下初始構(gòu)型，再從頭優(yōu)化。對(duì)于L9999快達(dá)到收斂時(shí)，考慮減小優(yōu)化步長(zhǎng)有時(shí)對(duì)于能量振蕩的情況也是有用的，opt(maxstep=1).(flyingheart)/chemistry/resources/messages/2002/10/16.005-dir/index.htmlHello,hereisthesummaryofalltheanswersIgotformyproblem.Thanks(inalphabeticallyorder)toLaurenceCuffe,KerwinDobbs,PaulKiprof,LiviuMirica,CoryPyeandDaveShobe.TheproblemwassolvedbyincreasingthenumberofSCFcycles,assuggestedbymany.Thisseemsveryoftennecessarywithtransitionmetals.Thetwoothersuggestions(usingfirstasmallbasissetbeforeswitchingtothelarger,orturningoff"vshift",didnotwork,unfortunately.Thanksagain,R.Herrmann.Theoriginalquestion:Hello,IrecentlytriedtocalculatetheNMRchemicalshiftsofVOCl3withGaussian98:#B3LYP/6-31G*SPNMRVOCl301clcl1clcl2cl2clcl31clclcl3o3ocl42oclcl41dih4v3vcl52vclcl51dih5clcl33.506243clclcl359.999ocl42.938735oclcl453.376dih4-64.586vcl52.118170vclcl534.144dih5-31.640Theresult(evenwithseveralvariationsintheinput)isalwaysanerrorterminationwhentryingthe"incrementalFockformation”(convergencecriterionnotmet).ThestructurewasoptimizedwithGAMESS.AnysuggestionwhatIcoulddoiswellcome-thankyou!108basisfunctions288primitivegaussians41alphaelectrons41betaelectronsnuclearrepulsionenergy561.9817715949Hartrees.One-electronintegralscomputedusingPRISM.NBasis=108RedAO=TNBF=108NBsUse=1081.00D-04NBFU=108ProjectedHuckelGuess.RequestedconvergenceonRMSdensitymatrix=1.00D-08within64cycles.RequestedconvergenceonMAXdensitymatrix=1.00D-06.Virtualorbitalswillbeshiftedby0.200hartree.RestartingincrementalFockformation.RestartingincrementalFockformation.RestartingincrementalFockformation.>>>>>>>>>>Convergencecriterionnotmet.SCFDone:E(RB+HF-LYP)=-2359.59549049A.U.after65cyclesConvg=0.5851D-03-V/T=1.9657S**2=0.0000Convergencefailure--runterminated.ErrorterminationviaLnk1ein/usr/local/g98_a7/l502.exe.Jobcputime:0days0hours16minutes33.0seconds.Filelengths(MBytes):RWF=14Int=0D2E=0Chk=1Scr=Theanswers:From:KerwinDDobbs<Kerwin.D.Dobbs-:at:-USA.>To:"RudolfHerrmann<herrmann"<herrmann-:at:-itqb.unl.pt>Subject:Re:CCL:G98:VOCl3problem:convergencecriterionnotmetDr.Herrmann,Tryturningoffthe"vshift"bychangingtheline,#B3LYP/6-31G*SPNMRto1#B3LYP/6-31G*SPNMRIOP(5/10=-1)WhenthereisaTMinthemolecule,G98automaticallyturnson"VSHIFT"whichusuallymakesSCFconvergenceworse.Theaboveadditiontoyourlineturnsoff"VSHIFT".Giveitatryandseewhathappens.KerwinFrom:CoryPye<cpye-:at:-crux.stmarys.ca>To:RudolfHerrmann<herrmann-:at:-itqb.unl.pt>Subject:Re:CCL:G98:VOCl3problem:convergencecriterionnotmetHello,yourSCFisnotconverging.Tryrunning,atfirst,HF/STO-3G.Ifitworks,thenThentryB3LYP/6-31G*afterwardswithguess=readintheinput.Itusuallyworkslikeacharmforme,Ialwaysusesomethinglikethis.-CoryFrom:LiviuMihailMirica<liviu-:at:-Stanford.EDU>To:RudolfHerrmann<herrmann-:at:-itqb.unl.pt>Subject:Re:CCL:G98:VOCl3problem:convergencecriterionnotmetTryaddingSCFcyc=1000inyourcommandlines.Thiswillincreasethe#ofSCFcyclesuptothedesirednumber(1000inthiscase).thatshouldbeenoughforyoursystem.—From:LaurenceCuffe<Laurence.Cuffe-:at:-ucd.ie>To:RudolfHerrmann<herrmann-:at:-itqb.unl.pt>Subject:Re:CCL:G98:VOCl3problem:convergencecriterionnotmetAnditshouldrunthesekeywordschangethedefaultSCFmethodandthenincreasethenumberofSCFcyclesallowed.Igotthejobtocompleteinabout50minutesrunningitincoreonalinuxboxwithamemorylimitof600MBbutitthisdoesn'tmakesensetoyoudon'tworryaboutit.Tryusingthemodifiedlineaboveandifthereareanyproblemsdon'thesitatetocomebacktome.AllthebestLaurenceCuffeFrom:PaulKiprof<pkiprof-:at:->To:herrmann-:at:-itqb.unl.ptSubject:VOCL3DearDr.Herrmann,IwouldaddtotheinputSCF(direct,qc,maxcycle=1000)Thiswillincreasethenumberofconvergencecyclesandhopefullyconvergethecalculation.IalwaysdothatwhenIhavetransitionmetalcompounds.Bestwishes,PaulKiprofFrom:"Shobe,Dave"<dshobe-:at:->To:'RudolfHerrmann'<herrmann-:at:-itqb.unl.pt>Subject:RE:G98:VOCl3problem:convergencecriterionnotmetThestandard64SCF"cycles”(oriterationsofthevariationalprocedure)israrelyenoughfortransitionmetalcompounds,evenclosedshelloneslikeVOCl3.ReadthesectionintheGaussianmanualdescribingthe"scf”keyword,andyouwillfindvariousoptionswhichmayimprovetheconvergence.Alotofitistrialanderror,withthesedifferentscfoptions.Thesimplestisscf(maxcyc=999)whichincreasesthenumberofcycles/iterationsto999beforetheprogram"givesup"andprintstheerrormessage.Scf=qcandscf(vshift=###)havealsobeenusefulinmyexperience.Moreadvice:onceyougetagoodwavefunction,makeabackupofthecheckpointfile,soifthewavefunctionbecomescorruptedinsubsequentcalculationsyoucanalwaysgobacktoit.Dr.RudolfHerrmannInstitutodeTecnologiaQuimicaeBiologica(ITQB)Av.daRepublica,EAN,Apt.127P-2781-901Oeirase-mail:herrmann-:at:-itqb.unl.ptgaussian學(xué)習(xí)筆記(一)默認(rèn)分類2008-11-1717:41:53閱讀124評(píng)論0字號(hào)：大中小訂閱-- 用Gaussian計(jì)算不收斂的問題問題一：Problemdetectedwithinexpensiveintegrals.Switchingtofullaccuracyandrepeatinglastcycle.Cycle2Pass1IDiag1:Warning!Spuriousintegrateddensity:NE=188NElCor=0Integral=-76.06821Tolerance=1.00D-03Consistencyfailure#2inCalDSu.ErrorterminationviaLnk1ein/home/gaussian/g98/l502.exe.Jobcputime:1days3hours21minutes30.2seconds.Filelengths(MBytes):RWF=118Int=0D2E=0Chk=17Scr=1解決：(fromccl)AsfarasIknowit'sadft-relatederror.Gaussianchecksthechargeresultingfromtheintegrationofthedensitytothenumberofelectrons,andincaseofalargedifferenceprintsthiserormessage.Tryadensergridwiththe"integral"keyword.可是我采用integral(Grid=ultrafine)網(wǎng)格仍是出現(xiàn)以下錯(cuò)誤信息：Raffenetti2integralformat.Two-electronintegralsymmetryisturnedoff.306basisfunctions778primitivegaussians51alphaelectrons51betaelectronsnuclearrepulsionenergy865.0704012775Hartrees.One-electronintegralscomputedusingPRISM.1SymmetryoperationsusedinECPInt.ECPInt:NShTT=17391NPrTT=123061LenC2=14761LenP2D=57336.LDataN:DoStor=FMaxTD1=3Len=28LDataN:DoStor=TMaxTD1=3Len=28NBasis=306RedAO=TNBF=306NBsUse=3061.00D-04NBFU=306ProjectedHuckelGuess.ofinitialguess=0.0000RequestedconvergenceonRMSdensitymatrix=1.00D-06within800cycles.RequestedconvergenceonMAXdensitymatrix=1.00D-04.Integralaccuracyreducedto1.0D-05untilfinaliterations.Warning!Spuriousintegrateddensity:NE=102NElCor=0Integral=102.20958Tolerance=1.00D-03Consistencyfailure#2inCalDSu.ErrorterminationviaLnk1ein/usr/people/gaussian/g98/l502.exe.Jobcputime:0days0hours13minutes24.7seconds.Filelengths(MBytes):RWF=90Int=0D2E=0Chk=2Scr=1請(qǐng)問各位大俠有沒有什么解決方法？問題二：>>>>>>>>>>Convergencecriterionnotmet.SCFDone:E(UB-PW91)=-1306.37542885A.U.after801cyclesConvg=0.1221D+01-V/T=2.0190S**2=0.0000KE=1.281978320748D+03 PE=-8.302694580633D+03EE=3.210822249139D+03Annihilationofthefirstspincontaminant:S**2beforeannihilation0.0000,after0.0000Convergencefailure--runterminated.ErrorterminationviaLnk1ein/usr/people/gaussian/g98/l502.exe.Jobcputime:2days2hours6minutes40.2seconds.Filelengths(MBytes):RWF=97Int=0D2E=0Chk=19Scr=1一般而言，解決收斂的問題有以下幾種，第一，波函數(shù)初猜，不穩(wěn)定的波函數(shù)，積分難以收斂是肯定的，解決方案是作stable=opt，做完以后，在此基礎(chǔ)上從check文件中讀入波函數(shù)?！猯ink1——%chk#phfoptguess=readgeom=check第二，過渡金屬原子會(huì)出現(xiàn)d軌道分裂問題，如果分裂能不夠大，電子顯然會(huì)在高能軌道和低能軌道之間跳動(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考慮分子間相互作用，可是好多都出現(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ī)則，陽離子比陰離子更容易收斂。選項(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í)移動(dòng)(SCF=Vshift)。如果接近SCF但未達(dá)到，收斂標(biāo)準(zhǔn)就會(huì)放松或者忽略收斂標(biāo)準(zhǔn)。這通常用于不是在初始猜測(cè)而是在平衡結(jié)構(gòu)收斂的幾何優(yōu)化。SCF=Sleazy放