高斯收敛问题

高斯收敛问题首先，我们必须理解收敛是什么意思。在自洽场（SCF）计算中，自洽循环中，首先产生一个轨道占据的初始猜测，1） 然后根据此轨道占据构造电荷密度和哈密顿量。2） 对角化哈密顿量，得到新的轨道能级和占据。3） 产生新的电荷分布和哈密顿量，重复步骤2）经过一定次数的循环后，某次循环前和循环后的电荷密度差别小于一定的标准，我们称之为收敛。如果以上过程不能收敛，则gaussian给出convergencefailure的警告。如果SCF计算收敛失败，你首先会采取哪些技巧呢？这里是我们强烈推荐的首选方法。1考虑使用更小的基组由于一定的基组对应于一定精度和速度，所以更换基组并不在所有的情况下都适用。方法是首先用小基组进行计算，由前一个波函得到用于大基组计算的初始猜测（Guess=Read自动进行）。2增加最大循环步数Gaussian默认的最大循环步数为64（SCF=DM或SCF=QC方法则为512），如果循环次数超过这个数目则会汇报convergencefailure。在一定的情况下，不收敛的原因仅仅是因为最大循环步数不够。可以通过设置maxcyc来增大最大循环步数。更多的SCF迭代（SCF（MaxCycle=N），其中N是迭代数）。这很少有帮助，但值得一试。3放宽收敛标准如果接近SCF但未达到，收敛标准就会放松或者忽略收敛标准。这通常用于不是在初始猜测而是在平衡结构收敛的几何优化。SCF=Sleazy放松收敛标准，Conver选项给出更多的控制。在Gaussian中设置收敛标准的命令行是Conv=8表示10-8，你可以将这个标准降低，比如改为10-5等等。微调结构是比较有效的方法，能量收敛标准不宜太低，最好不低于10-7.4尝试改变初始构型首先略微减小键长，接下来略微增加键长，接下来再对结构作一点改变。5尝试能级移动Levelshifting(SCF=Vshift)如果不收敛的原因是波函数的震荡行为，通常是因为在相近的能量上的泰的混合。对于这种情况，我们可以采用levelshifting的方法。Levelshifting的含义是人工的升高非占据轨道的能级，以防止和最高占据轨道之间的混合，以达到收敛的目的。在Gaussian中此方法的关键词为SCF=Vshift6使用强制的收敛方法SCF=QCSCF=QC通常最佳，但在极少数情况下SCF=DM更快。此关键字将大大增加计算时间，但是收敛的机会更大。不要忘记给计算额外增加一千个左右的迭代。应当测试这个方法获得的波函，保证它最小，并且正好不是稳定点(使用Stable关键字)。7对开壳层体系，尝试收敛到同一分子的闭壳层离子，接下来用作开壳层计算的初始猜测。添加电子可以给出更合理的虚轨道，但是作为普遍的经验规则，阳离子比阴离子更容易收敛。选项Guess=Read定义初始猜测从Gaussian计算生成的checkpoint文件中读取。8一些程序通过减小积分精度加速SCF。对于使用弥散函数，长程作用或者低能量激发态的体系，必须使用高积分精度：SCF=NoVarAcc。9改变模型或方法可以考虑改变模型方法。比较常见的方法有HF,GVB,MCSCF,CASSCF,MPn等。改变模型方法通常也会收敛性质。通常，精度更高的方法更难收敛。精度比较低的方法产生的计算结果可以作为高精度计算的初始猜测。考虑使用不同理论级别的计算。这并不总是实用的，但除此之外，增加迭代数量总是使得计算时间和使用更高理论级别差不多。10关闭DIIS外推(SCF=NoDIIS)。同时进行更多的迭代(SCF=(MaxCycle=N)。11试着改用DIIS之外其它方法(SCF=SD或SCF=SSD)。Gaussian不收敛的可能原因及对策：1由于体系有很多能量相近的能级，导致计算不收敛。如果计算中采用的是Hartree-Fock方法或者其他的混合形式的交换相关势(如B3LYP)，则可以尝试一下的方法进行改进。scf=(noincfock,conver=11,maxcyc=1025)iop(5/22=20)用以上的命令行将强制Gaussian采用EDIIS的算法。此算法计算量更大，但是更加稳定。在命令行里也增加了最大循环数，以增加收敛的到基态的可能性。检查是否有初始文件错误常见初级错误：自旋多重度错误变量赋值为整数变量没有赋值或多重赋值键角小于等于0度，大于等于180度分子描述后面没有空行二面角判断错误，造成两个原子距离过近分子描述一行内两次参考同一原子，或参考原子共线3.SCF（自洽场）不收敛则一般是L502错误，省却情况做64个cycle迭代（G03缺省128cycles）修改坐标，使之合理改变初始猜Guess=Huckel或其他的，看Guess关键词。增加叠代次数SCFCYC=N（对小分子作计算时最好不要增加，很可能结构不合理）d.iop（5/13=1）这样忽略不收敛，继续往下做。分子对称性改变修改坐标，强制高对称性或放松对称性给出精确的、对称性确定的角度和二面角。如CH4的角度给到109.47122放松对称性判据Symm=loose不做对称性检查iop（2/16=1）（最好加这个选项）iop（2/16=2）则保持新的对称性来计算Opt时收敛的问题修改坐标，使之合理增加叠代次数optcyc=N优化过渡态，若势能面太平缓，则不好找到。到草窣IdOd由也。髭苗玄暮囹甲目会好与暮香至耆苗誓果吊IdOd-0T8IAI000乙宜召也’果吊号此乙宜场爰扁身0内¥¥¥¥'此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区宜［M（DdM44±5±）M44asod由瘤44凿密qg宜却曾（iuepunpaypouj/si）=ldojo-由累旦会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□$中累旦（翠应）D丑W（aiepdnisruiou）ido叩匿〔回'习争阳\「/［®由查腥等一母凿□$°0■［御旎也’（V&0暗鬓阳欧旱裕〔也欧弱£一"）0£*¥兽（OT=8/T）doi用OPT11.优化过渡态只能做一个STEP原因是负本征数目不对添加iop（1/11）=1或者noeigentest（eigentest是表示优化过渡态检测分子振动的本征值，过渡态只有一个负值，但优化的时候往往出现两个或者更多的情况，默认的是如果出现多于一个的情况就停止优化计算，这往往时不必要的，noeigentest表示优化时不坐此检测，实际上优化过渡态的时候出现两个三个很正常的，只是第一个一般负的很大，接着的比较小，不会影响过渡态的搜寻，当然有时候需要自己判断是不是接近真实的过渡态。其实经验上也是过渡态的本征虚频越小，相对越难找。）对于（L502,L508,L9999）出错的对策对于一个优化计算，它的过程是先做一个SCF计算，得到这个构型下的能量，然后优化构型，再做SCF，然后再优化构型。。。因此，会有两种不收敛的情况：一是在某一步的SCF不收敛（L502错误），或者构型优化没有找到最后结果（L9999错误）。预备知识:计算时保存chk文件，可以在后续计算中使用guess=read读初始猜测.对于SCF不收敛，通常有以下的解决方法：使用小基组，或低级算法计算，得到scf收敛的波函数，用guess=read读初始波函数。使用scf=qc，这个计算会慢，而且需要用stable关键字来测试结果是否波函数稳定。如果这个还不收敛，会提示L508错误。改变键长，一般是缩小一点，有时会有用。计算相同体系的其他电子态，比如相应的阴离子、阳离子体系或单重态体系，得到的收敛波函数作为初始猜测进行计算。对于优化不收敛，即L9999错误，实际上是在规定的步数内没有完成优化，即还没有找到极小值点。(或者对于过渡态优化，还没有找到过渡态)这有几种可能性：看一下能量的收敛的情况，可能正在单调减小，眼看有收敛的趋势，这样的情况下，只要加大循环的步数(opt(maxcycle=200)),可能就可以解决问题了。加大循环步数还不能解决的(循环步数有人说超过200再不收敛，再加也不会有用了，这虽然不一定绝对正确，但200步应该也差不多了)，有两种可能。一是查看能量，发现能量在振荡了，且变化已经很小了，这时可能重新算一下，或者构型稍微变一下，继续优化，就可以得到收敛的结果(当然也有麻烦的，看运气和经验了)；二是构型变化太大，和你预计的差别过大，这很可能是你的初始构型太差了，优化不知道到哪里去了，这时最好检查一下初始构型，再从头优化。对于L9999快达到收敛时，考虑减小优化步长有时对于能量振荡的情况也是有用的，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学习笔记(一)默认分类2008-11-1717:41:53阅读124评论0字号：大中小订阅-- 用Gaussian计算不收敛的问题问题一：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)网格仍是出现以下错误信息：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请问各位大侠有没有什么解决方法？问题二：>>>>>>>>>>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一般而言，解决收敛的问题有以下几种，第一，波函数初猜，不稳定的波函数，积分难以收敛是肯定的，解决方案是作stable=opt，做完以后，在此基础上从check文件中读入波函数。——link1——%chk#phfoptguess=readgeom=check第二，过渡金属原子会出现d轨道分裂问题，如果分裂能不够大，电子显然会在高能轨道和低能轨道之间跳动，能量会上下波动，难以收敛，解决方案是加入vshift关键词第三，氢键，过渡金属，默认的129cycle所难以达到收敛，所以加大scf(cyc=500)可能可以解决问题。第四，初始结构一直容易被人所忽略，所以很多方法没有解决问题的时候，还是希望能仔细检查一遍结构第五，scf的时候，检查一下力常数也是好的，通常的解决方案作一下calcfc对于第二个收敛的问题，我正在用SCF(vshift,QC)尝试，得过几天才知道行不行。至少我已尝试了将内坐标改为直角坐标；加大循环次数(800仍不收敛)从check文件中读入，这些都不行。对于第一个问题，我是用massage考虑分子间相互作用，可是好多都出现第一个问题，根本不让你算，用了integral(grid=ultrafine)还是不行，输入文件如下：(这是一个不含金属原子的，有好几个含金属的输入，也是出现这个情况。芝麻582:1.在Guess关键字中使用Core,Huckel或Mix选项，试验不同的初始猜测。2.对开壳层体系，尝试收敛到同一分子的闭壳层离子，接下来用作开壳层计算的初始猜测。添加电子可以给出更合理的虚轨道，但是作为普遍的经验规则，阳离子比阴离子更容易收敛。选项Guess=Read定义初始猜测从Gaussian计算生成的checkpoint文件中读取。3.另一个初始猜测方法是首先用小基组进行计算，由前一个波函得到用于大基组计算的初始猜测(Guess=Read自动进行)。4.尝试能级移动(SCF=Vshift)。如果接近SCF但未达到，收敛标准就会放松或者忽略收敛标准。这通常用于不是在初始猜测而是在平衡结构收敛的几何优化。SCF=Sleazy放