【分享】高斯对不收敛问题的对策；首先，我们必须理解收敛是什么意思；1)然后根据此轨道占据构造电荷密度和哈密顿量；2)对角化哈密顿量，得到新的轨道能级和占据；3)产生新的电荷分布和哈密顿量，重复步骤2）；经过一定次数的循环后，某次循环前和循环后的电荷密；如果SCF计算收敛失败，你首先会采取哪些技巧呢？；由于一定的基组对应于一定精度和速度，所以更换基组；2增加最
【分享】高斯对不收敛问题的对策
首先，我们必须理解收敛是什么意思。在自洽场(SCF)计算中，自洽循环中，首先产生一个轨道占据的初始猜测，
1)然后根据此轨道占据构造电荷密度和哈密顿量。
2)对角化哈密顿量，得到新的轨道能级和占据。
3)产生新的电荷分布和哈密顿量，重复步骤2）
经过一定次数的循环后，某次循环前和循环后的电荷密度差别小于一定的标准，我们称之为收敛。
如果以上过程不能收敛，则gaussian给出convergence failure的警告。
如果SCF计算收敛失败，你首先会采取哪些技巧呢？这里是我们强烈推荐的首选方法。 1 考虑使用更小的基组
由于一定的基组对应于一定精度和速度，所以更换基组并不在所有的情况下都适用。方法是首先用小基组进行计算，由前一个波函得到用于大基组计算的初始猜测(Guess=Read自动进行)。
2 增加最大循环步数
Gaussian默认的最大循环步数为64 (SCF=DM或SCF=QC方法则为512)，如果循环次数超过这个数目则会汇报convergence failure。在一定的情况下，不收敛的原因仅仅是因为最大循环步数不够。可以通过设置maxcyc来增大最大循环步数。更多的SCF迭代(SCF(MaxCycle=N)，其中N是迭代数)。这很少有帮助，但值得一试。
放宽收敛标准
如果接近SCF但未达到，收敛标准就会放松或者忽略收敛标准。这通常用于不是在初始猜测而是在平衡结构收敛的几何优化。SCF=Sleazy放松收敛标准，Conver选项给出更多的控制。
尝试改变初始构型
首先略微减小键长，接下来略微增加键长，接下来再对结构作一点改变。
尝试能级移动Level shifting (SCF=Vshift)
如果不收敛的原因是波函数的震荡行为，通常是因为在相近的能量上的泰的混合。对于这种情况，我们可以采用level shifting的方法。Level shifting的含义是人工的升高非占据轨道的能级，以防止和最高占据轨道之间的混合，以达到收敛的目的。在Gaussian中此方法的关键词为SCF=Vshift
使用强制的收敛方法SCF=QC
SCF=QC通常最佳，但在极少数情况下SCF=DM更快。此关键字将大大增加计算时间，但是收敛的机会更大。不要忘记给计算额外增加一千个左右的迭代。应当测试这个方法获得的波函，保证它最小，并且正好不是稳定点(使用Stable关键字)。
对开壳层体系，尝试收敛到同一分子的闭壳层离子，接下来用作开壳层计算的初始猜测。添加电子可以给出更合理的虚轨道，但是作为普遍的经验规则，阳离子比阴离子更容易收敛。选项Guess=Read定义初始猜测从Gaussian计算生成的checkpoint文件中读取。
一些程序通过减小积分精度加速SCF。对于使用弥散函数，长程作用或者低能量激发态的体系，必须使用高积分精度：SCF=NoVarAcc。
改变模型或方法
可以考虑改变模型方法。比较常见的方法有HF,GVB,MCSCF,CASSCF,MPn等。改变模型方法通常也会收敛性质。通常，精度更高的方法更难收敛。精度比较低的方法产生的计算结果可以作为高精度计算的初始猜测。考虑使用不同理论级别的计算。这并不总是实用的，但除此之外，增加迭代数量总是使得计算时间和使用更高理论级别差不多。
关闭DIIS外推(SCF=NoDIIS)。同时进行更多的迭代(SCF=(MaxCycle=N))。
试着改用DIIS之外其它方法(SCF=SD或SCF=SSD)。
Gaussian不收敛的可能原因及对策：
1 由于体系有很多能量相近的能级，导致计算不收敛。
如果计算中采用的是Hartree-Fock方法或者其他的混合形式的交换相关势（如B3LYP）,则可以尝试以下的方法进行改进。
scf=(noincfock,conver=11,maxcyc=1025) iop(5/22=20)
用以上的命令行将强制Gaussian采用EDIIS的算法。此算法计算量更大，但是更加稳定。在命令行里也增加了最大循环数，以增加收敛的到基态的可能性。
2． 检查是否有初始文件错误
常见初级错误：
a. 自旋多重度错误
b. 变量赋值为整数
c. 变量没有赋值或多重赋值
d. 键角小于等于0度，大于等于180度
e. 分子描述后面没有空行
f. 二面角判断错误，造成两个原子距离过近
g. 分子描述一行内两次参考同一原子,或参考原子共线
3． SCF(自洽场)不收敛
则一般是L502错误，省却情况做64个cycle迭代(G03缺省128 cycles)
a. 修改坐标，使之合理
b. 改变初始猜
Guess=Huckel 或其他的，看Guess关键词。
c. 增加叠代次数 SCFCYC=N （对小分子作计算时最好不要增加，很可能结构不合理）
d. iop（5/13=1）这样忽略不收敛，继续往下做。
4． 分子对称性改变
a. 修改坐标，强制高对称性或放松对称性
b. 给出精确的、对称性确定的角度和二面角。 如CH4的角度给到109.47122
c. 放松对称性判据
Symm=loose
d. 不做对称性检查 iop(2/16=1) （最好加这个选项）
iop(2/16=2) 则保持新的对称性来计算
5． Opt时收敛的问题
a. 修改坐标，使之合理
b. 增加叠代次数optcyc=N
6． 优化过渡态，若势能面太平缓，则不好找到。
iop(1/8=10) 默认30（下一个结构和该结构的差别0.3?），可改成10。如果每一步都要用到小的步长，应该加opt(notrustupdate)
7． 在CI（组态）方法中
如QCISD(T),CCSD(T),CID方法中，省却最大循环50，若出错（L913错误）
解决方法：#P QCISD(maxcyc=N)
注：N≤512
8． 优化过渡态
opt=TS (给出过渡态)
opt=qst2 (给出反应物和产物)
opt=qst3 (给出反应物和产物和过渡态)
用G03时的出错
opt=ts 必须加FC (force constant)
写法：opt=(TS, calcFc)
opt=(TS,calchffc)计算HF力常数，对QCISD，CCSD等方法用；
opt=(TS,modRedundant)
(最好写这个)
b. 如果计算采用QCISD计算（不好计算FC）
则写为QCISD opt=(TS, calcHFFC) (用HF计算FC)
9. 无法写大的Scratch文件RWF
a. 劈裂RWF文件 %rwf=loc1,size1,loc2,size2,……..,locN,-1
b. 改变计算方法 MP2=Direct可以少占硬盘空间
c. 限制最大硬盘 maxdisk=N GB，****MB,有些系统写2GB会出错，可以写2000MB
10. FOPT出错 原因是变量数与分子自由度数不相等。 可用POPT 或直接用OPT
11. 优化过渡态只能做一个STEP 原因是负本征数目不对 添加 iop（1/11）=1或者noeigentest (eigentest 是表示优化过渡态检测分子振动的本征值，过渡态只有一个负值，但优化的时候往往出现两个或者更多的情况，默认的是如果出现多于一个的情况就停止优化计算，这往往时不必要的，noeigentest表示优化时不坐此检测，实际上优化过渡态的时候出现两个三个很正常的，只是第一个一般负的很大，接着的比较小，不会影响过渡态的搜寻，当然有时候需要自己判断是不是接近真实的过渡态。其实经验上也是过渡态的本征虚频越小，相对越难找。)
对于 (L502, L508, L9999)出错的对策
对于一个优化计算，它的过程是先做一个SCF计算，得到这个构型下的能量，然后优化构型，再做SCF，然后再优化构型。。。因此，会有两种不收敛的情况：一是在某一步的SCF不收敛(L502错误)，或者构型优化没有找到最后结果（L9999错误）。
预备知识:计算时保存chk文件,可以在后续计算中使用guess=read读初始猜测.
对于SCF不收敛，通常有以下的解决方法：
1. 使用小基组，或低级算法计算，得到scf收敛的波函数，用guess=read读初始波函数。
2. 使用scf=qc，这个计算会慢，而且需要用stable关键字来测试结果是否波函数稳定。如果这个还不收敛，会提示L508错误。
3. 改变键长，一般是缩小一点，有时会有用。
4. 计算相同体系的其他电子态，比如相应的阴离子、阳离子体系或单重态体系，得到的收敛波函数作为初始猜测进行计算。
对于优化不收敛,即L9999错误，实际上是在规定的步数内没有完成优化，即还没有找到极小值点。（或者对于过渡态优化，还没有找到过渡态）
这有几种可能性：
1. 看一下能量的收敛的情况，可能正在单调减小，眼看有收敛的趋势，这样的情况下，只要加大循环的步数（opt(maxcycle=200)），可能就可以解决问题了。
2. 加大循环步数还不能解决的（循环步数有人说超过200再不收敛，再加也不会有用了，这虽然不一定绝对正确，但200步应该也差不多了），有两种可能。一是查看能量，发现能量在振荡了，且变化已经很小了，这时可能重新算一下，或者构型稍微变一下，继续优化，就可以得到收敛的结果（当然也有麻烦的，看运气和经验了）；二是构型变化太大，和你预计的差别过大，这很可能是你的初始构型太差了，优化不知道到哪里去了，这时最好检查一下初始构型，再从头优化。
3. 对于L9999快达到收敛时,考虑减小优化步长有时对于能量振荡的情况也是有用的，opt(maxstep=1).(flyingheart )
三亿文库包含各类专业文献、中学教育、应用写作文书、各类资格考试、文学作品欣赏、生活休闲娱乐、外语学习资料、幼儿教育、小学教育、专业论文、17高斯计算中优化不收敛对策等内容。
(每 个轨道 1 个电子) 在高斯中,以 HF 为例,闭、开、限制开壳层计算分别...不收敛(L502 错误),或者构型优化没有找到最后结果(L9999 错误) 。 预备知识:... . Central South University 结构化学之高斯计算学院:...一种可能是,优化的几何结构不够精确,还没有达到稳...第八,优 化一个局域极小点, 收敛后从力常数本征... (每个轨道 1 个电子) 在高斯中,以 HF 为例,闭、开、限制开壳层计算分别是...5. 待补充. 对于优化不收敛,即 L9999 错误,实际上是在规定的步数内没有... 分子性质研究的一般模式: 分子性质研究的一般模式: 2 高斯中所用到的一些术语...计算收敛的标准 2.3 计算收敛的标准优化计算不能无限制的进行下去, 判定是否... 高斯牛顿_计算机软件及应用_IT/计算机_专业资料。很好...即无约束最优化问题 n 是变量数目,目标函数 f(x)...对于有些使用高斯牛顿方法局部不收敛的问题, 如果没... (每个轨 道 1 个电子) 在高斯中,以 HF 为例,闭、开、限制开壳层计算分别...不收敛(L502 错误),或者构型优化没有找到最后结果(L9999 错误) 。 预备知识:... 对于含有一些较高 R/X 值支路的网络,特别是系统中有大量的电缆线路时,应用高斯-赛德 尔法收敛性较好。当出现潮流计算不收敛时,可以选择快速解耦法,并关闭自动... 高斯常见问题_军事/政治_人文社科_专业资料。一、混合...态的优化,使得过渡态不至于收敛到其他的分子结构中...其次在计算上 还可以尝试:加大循环次数;提高... 高斯常见错误_理学_高等教育_教育专区。最近总结的高斯计算中常见的错误近来...对于优化不收敛,即 L9999 错误,实际上是在规定的步数内没有完成优化,即还没...您所在位置： &
 &  & 
高斯错误修改总结重点.doc19页
本文档一共被下载：
次 ,您可全文免费在线阅读后下载本文档。
文档加载中...广告还剩秒
需要金币：350 &&
你可能关注的文档：
··········
A list of error messages and possible solutions -
Gaussian calculations can fail with various error messages. Some error messages from .out and .log files - and possible solutions - have been compiled here to facilitate problem solving.-
These are divided into:-
Syntax and similar errors- 语法类错误
Memory and similar errors-
内存类错误
Convergence problems - 不收敛错误
Errors in solvent calculations - 溶剂中的计算错误
Errors in log files- 错误文件 -
ERROR MESSAGES IN OUTPUT FILES-
Syntax and similar errors:End of file in ZSymb.-
Error termination via Lnk1e in /global/apps/gaussian/g03.e01/g03/l101.exeSolution: The blank line after the coordinate section in the .inp file is missing.
输入文件空行丢失
Unrecognized layer "X".-（不识别层X）
Error termination via Lnk1e in /global/apps/gaussian/g03.e01/g03/l101.exeSolution: Error due to syntax error s
in coordinate section
check carefully . If error is "^M", it is caused by DOS end-of-line characters
e.g. if coordinates were written under Windows . Remove ^M from line ends using e.g. emacs. To process .inp files from command line, use sed -i
's/^M//' File.inp
Important: command does not work if ^M is written as characters - generate ^M on command line using ctrl-V ctrl-M .-
QPERR --- A SYNTAX ERROR WAS DETECTED IN THE INPUT LINE.-
Solution: Check .inp carefully for syntax errors in keywords -
Unable to locate IRWF 0 Number
Error termination via Lnk1e in /global/apps/gaussian/g03.e01/g03/l401.exe or-
FileIO operation on non-existent file.-
[...] Error termination in NtrErr:-
NtrErr Called from FileIO.Solution: Operation on .chk file was specified
e.g. geom check, opt restart , but .chk was not found. Check that:-
was specifed in .inp-
.chk has the same name as .inp-
.chk is in the same directory as .inp -
run script transports .chk to temporary folder upon job start. Run scripts downloaded here should do this. -
The combination
正在加载中，请稍后...您所在位置： &
 &  & 
高斯计算常见错误和解决方案.doc3页
本文档一共被下载：
次 ,您可全文免费在线阅读后下载本文档。
文档加载中...广告还剩秒
需要金币：100 &&
你可能关注的文档：
··········
··········
GAUSSION计算常见错误及解决方案 1. 自旋多重度错误 2. 变量赋值为整数 3. 变量没有赋值 4. 键角小于等于0度，大于等于180度 5. 分子描述后面没有空行 6. 二面角判断错误，造成两个原子距离过近 7. 分子描述一行内两次参考同一原子,或参考原子共线 运行出错 1. 自洽场不收敛??SCF a. 修改坐标，使之合理 b. 改变初始猜??Guess c. 增加叠代次数 SCFCYC N d. iop（5/13 1） 2.??分子对称性改变 a. 修改坐标，强制高对称性或放松对称性 b. 给出精确的、对称性确定的角度和二面角 c. 放松对称性判据??Symm loose d. 不做对称性检查 iop 2/16 1 3. 无法写大的Scratch文件RWF a. 劈裂RWF文件 %rwf loc1,size1,loc2,size2,……..,locN,-1 b. 改变计算方法 MP2 Direct可以少占硬盘空间 c. 限制最大硬盘 maxdisk N GB 4. FOPT出错 原因是变量数与分子自由度数不相等。 可用POPT 或直接用OPT 5. 优化过渡态只能做一个STEP 原因是负本征数目不对 添加 iop（1/11） 1 6. 组态相互作用计算中相关能叠代次数不够，增加叠代次数 QCISD（Maxcyc N） Default.Rou设置 ? 在Scratch文件夹中的Default.Rou文件中设置G03程序运行的省缺参数： ? -M-??200MW ? -P- 4 ? -#- MaxDisk 10GB ? -#- SCF Conventional or Direct ? -#- MP2 NoDirect or??Direct ? -#- OPTCYC 200 ? -#- SCFCYC 200 ? -#- IOPs 设置 如iop（2/16 1） Default.Rou设置中的冲突 ? Default route:??MaxDisk 2GB SCF Direct MP2 Direct OPTCYC 200 SCFcyc 100 iop 2/16 1
iop 5/13 1 ? ------------------ ? # ccsd/6-31G** opt ? ------------------ ? L903/L905 and L906 can only do MP2. 问题在于，ＭP2 Direct ！ ??去掉这个设置，CCSD的作业就能进行了。 因此，建议在Default设置中只设置，内存，最大硬盘，等项 振动分析的关键词 振动分析的选项 拉曼光谱 Raman 非谐性 Anharmonic 势能函数不是简谐的 内转动，受阻转子 Hinderedrotor 力常数过小
正在加载中，请稍后...小木虫 --- 600万学术达人喜爱的学术科研平台
&&查看话题
紧急求助啊，跪求大神啊，高斯计算出现错误，Convergence failure--run terminated
不懂高斯，很苦恼，计算了三次了都是出现错误
路过的都看看帮帮忙~
谢谢美女:D
关于紧急求助啊，跪求大神啊，高斯计算出现错误，Convergence failure--run terminated的相关话题在小木虫APP已经有18位虫友给出了详细回复。
赶快查看回复吧！
学术必备与600万学术达人在线互动！
扫描下载送金币
北京学而思教育科技有限公司 地址：北京市海淀区北三环甲18号中鼎大厦A座1层102室 电话：010-
浏览器进程
打开微信扫一扫
随时随地聊科研2.对于对称分子，设置scf(dsymm)强制使；3.对于对称分子，设置全局nosymm，降低对称；4.使用scf(novaracc)也可能会产生或；5.加上int(untrafine)，不过这个改；6.最佳解决方法：加关键词int(NoXCTes；Errorininternalcoordinat；可在opt中加入cartesian；例如：；%chk=TS3
2. 对于对称分子，设置scf(dsymm) 强制使用密度对称。
3. 对于对称分子，设置全局nosymm，降低对称性。
4. 使用scf(novaracc)也可能会产生或消除这个错误。
5. 加上int(untrafine)，不过这个改变了积分网格，这个计算出的能量不能与同类计算直接比较。
6. 最佳解决方法：加关键词int(NoXCTest) 网格问题(G09默认的精度检测方法不适用,加上Int=NoXCTest), Skip tests of numerical accuracy of XC quadrature (g09手册)。需要检测结果波函数是否正确。可同时加guess=indo。
Error in internal coordinate system.
可在opt中加入cartesian
%chk=TS3-D-itm3-eoc-4.chk
%mem=45000MB
%NProcShared=8
#B3LYP/6-31G** iop(5/13=1) iop(2/11=1)
#scf=(maxcycle=200)
#opt=(cartesian,TS,calcfc,noeigen,maxcycle=200) freq
#geom=allcheck
Error in internal coordinates
This topic covers how to solve the &Error in internal coordinates& problem.
Some times Gaussian quits at the Berny optimization stage with the &Error in
internal coordinates& message. This happens as Berny optimization by
default uses redundant internal coordinates. Cartesian coords are easy to define and use but can be strongly coupled to one another. Interal coords describe molecular properties (bond lengths, angle...) naturally and thus involve less coupling. There are a few types of redundant internal coords which can accelerate opt process. But when a molecule has flat geometry, it could be hard to calculate the reverse matrix of redundant internal coords. Therefore, it is better to switch back to Cartesian (opt=(...,Cartesian,...)) or use &nosymm&.
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad Berny optimization. NTrRot=
29 NAtoms=
5 IsLin=F Error in internal coordinate system. Error termination via Lnk1e in /home/cast/program/g09/l103.exe at Tue Nov 23 13:44:47 2010.
本人计算时出现如上错误，哪位达人知道是怎么回事？我觉得我的结构应该没问题 加上关键词OPT=Cartesian！
高斯对不收敛问题的对策
首先，我们必须理解收敛是什么意思。在自洽场(SCF)计算中，自洽循环中，首先产生一个轨道占据的初始猜测，
1)然后根据此轨道占据构造电荷密度和哈密顿量。
2)对角化哈密顿量，得到新的轨道能级和占据。 3)产生新的电荷分布和哈密顿量，重复步骤2） 经过一定次数的循环后，某次循环前和循环后的电荷密度差别小于一定的标准，我们称之为收敛。
如果以上过程不能收敛，则gaussian给出convergence failure的警告。 如果SCF计算收敛失败，你首先会采取哪些技巧呢？这里是我们强烈推荐的首选方法。 1 考虑使用更小的基组 由于一定的基组对应于一定精度和速度，所以更换基组并不在所有的情况下都适用。方法是首先用小基组进行计算，由前一个波函得到用于大基组计算的初始猜测(Guess=Read自动进行)。
2 增加最大循环步数
Gaussian默认的最大循环步数为64 (SCF=DM或SCF=QC方法则为512)，如果循环次数超过这个数目则会汇报convergence failure。在一定的情况下，不收敛的原因仅仅是因为最大循环步数不够。可以通过设置maxcyc来增大最大循环步数。更多的SCF迭代(SCF(MaxCycle=N)，其中N是迭代数)。这很少有帮助，但值得一试。
放宽收敛标准
如果接近SCF但未达到，收敛标准就会放松或者忽略收敛标准。这通常用于不是在初始猜测而是在平衡结构收敛的几何优化。SCF=Sleazy放松收敛标准，Conver选项给出更多的控制。在Gaussian中设置收敛标准的命令行是Conv=8 表示10-8,你可以将这个标准降低，比如改为10-5等等。微调结构是比较有效的方法，能量收敛标准不宜太低，最好不低于10-7.
尝试改变初始构型
首先略微减小键长，接下来略微增加键长，接下来再对结构作一点改变。
尝试能级移动Level shifting (SCF=Vshift) 如果不收敛的原因是波函数的震荡行为，通常是因为在相近的能量上的泰的混合。对于这种情况，我们可以采用level shifting的方法。Level shifting的含义是人工的升高非占据轨道的能级，以防止和最高占据轨道之间的混合，以达到收敛的目的。在Gaussian中此方法的关键词为SCF=Vshift 6
使用强制的收敛方法SCF=QC
SCF=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的算法。此算法计算量更大，但是更加稳定。在命令行里也增加了最大循环数，以增加收敛的到基态的可能性。
2． 检查是否有初始文件错误
常见初级错误：
a. 自旋多重度错误
b. 变量赋值为整数
c. 变量没有赋值或多重赋值
d. 键角小于等于0度，大于等于180度
e. 分子描述后面没有空行
f. 二面角判断错误，造成两个原子距离过近
g. 分子描述一行内两次参考同一原子,或参考原子共线
3． SCF(自洽场)不收敛
则一般是L502错误，省却情况做64个cycle迭代(G03缺省128 cycles)
a. 修改坐标，使之合理
b. 改变初始猜
Guess=Huckel 或其他的，看Guess关键词。
c. 增加叠代次数 SCFCYC=N （对小分子作计算时最好不要增加，很可能结构不合理） d. iop（5/13=1）这样忽略不收敛，继续往下做。
4． 分子对称性改变
a. 修改坐标，强制高对称性或放松对称性
b. 给出精确的、对称性确定的角度和二面角。 如CH4的角度给到109.47122
c. 放松对称性判据 Symm=loose
d. 不做对称性检查 iop(2/16=1) （最好加这个选项）
iop(2/16=2) 则保持新的对称性来计算
5． Opt时收敛的问题
a. 修改坐标，使之合理
b. 增加叠代次数optcyc=N
6． 优化过渡态，若势能面太平缓，则不好找到。
iop(1/8=10) 默认30（下一个结构和该结构的差别0.3?），可改成10。如果每一步都要用到小的步长，应该加opt(notrustupdate)
7． 在CI（组态）方法中
如QCISD(T),CCSD(T),CID方法中，省却最大循环50，若出错（L913错误）
解决方法：#P QCISD(maxcyc=N)
注：N≤512
8． 优化过渡态
opt=TS (给出过渡态)
opt=qst2 (给出反应物和产物)
opt=qst3 (给出反应物和产物和过渡态) a.
用G03时的出错
opt=ts 必须加FC (force constant)
写法：opt=(TS, calcFc)
opt=(TS,calchffc)计算HF力常数，对QCISD，CCSD等方法用；
opt=(TS,modRedundant)
(最好写这个)
b. 如果计算采用QCISD计算（不好计算FC）
则写为QCISD opt=(TS, calcHFFC) (用HF计算FC)
9. 无法写大的Scratch文件RWF
a. 劈裂RWF文件 %rwf=loc1,size1,loc2,size2,……..,locN,-1
b. 改变计算方法 MP2=Direct可以少占硬盘空间
c. 限制最大硬盘 maxdisk=N GB，****MB,有些系统写2GB会出错，可以写2000MB
10. FOPT出错 原因是变量数与分子自由度数不相等。 可用POPT 或直接用OPT
11. 优化过渡态只能做一个STEP 原因是负本征数目不对 添加 iop（1/11）=1或者noeigentest (eigentest 是表示优化过渡态检测分子振动的本征值，过渡态只有一个负值，但优化的时候往往出现两个或者更多的情况，默认的是如果出现多于一个的情况就停止优化计算，这往往时不必要的，noeigentest表示优化时不坐此检测，实际上优化过渡态的时候出现两个三个很正常的，只是第一个一般负的很大，接着的比较小，不会影响过渡态的搜寻，当然有时候需要自己判断是不是接近真实的过渡态。其实经验上也是过渡态的本征虚频越小，相对越难找。)
对于 (L502, L508, L9999)出错的对策
对于一个优化计算，它的过程是先做一个SCF计算，得到这个构型下的能量，然后优化构型，再做SCF，然后再优化构型。。。因此，会有两种不收敛的情况：一是在某一步的SCF不收敛(L502错误)，或者构型优化没有找到最后结果（L9999错误）。
预备知识:计算时保存chk文件,可以在后续计算中使用guess=read读初始猜测. 对于SCF不收敛，通常有以下的解决方法：
1. 使用小基组，或低级算法计算，得到scf收敛的波函数，用guess=read读初始波函数。
2. 使用scf=qc，这个计算会慢，而且需要用stable关键字来测试结果是否波函数稳定。如果这个还不收敛，会提示L508错误。
3. 改变键长，一般是缩小一点，有时会有用。
4. 计算相同体系的其他电子态，比如相应的阴离子、阳离子体系或单重态体系，得到的收敛波函数作为初始猜测进行计算。
对于优化不收敛,即L9999错误，实际上是在规定的步数内没有完成优化，即还没有找到极小值点。（或者对于过渡态优化，还没有找到过渡态）
这有几种可能性：
1. 看一下能量的收敛的情况，可能正在单调减小，眼看有收敛的趋势，这样的情况下，只要加大循环的步数（opt(maxcycle=200)），可能就可以解决问题了。
三亿文库包含各类专业文献、中学教育、行业资料、应用写作文书、幼儿教育、小学教育、16高斯错误修改总结等内容。
高斯常见错误及注意地方_能源/化工_工程科技_专业资料...,可改成 10。如果每一步都要用到小的步长,应该 ...高斯中的优化 29页 免费 总结:高斯量化计算总结 8... 高斯错误问题_计算机软件及应用_IT/计算机_专业资料。Erro r in int erna l ...《建筑工程管理与实务》笔记总结50份文档 2014年注册会计师考试 ... 绕某个键转动两个基团,有时 可以很方便地用修改二面角的方法实现:OPT=...高斯常见错误及注意地方 4页 1下载券 总结:高斯量化计算总结 8页 1下载券©... 总结:高斯量化计算总结_电子/电路_工程科技_专业资料。高斯量化计算总结总结...关于收敛问题 (L502, L508, L9999) 对于一个优化计算,它的过程是先做一个 ... 高斯模糊实现小结 高斯模糊是一种图像滤波器,它使用正态分布(高斯函数)计算模糊...Smooth2D 运行时间长,难道是二维数组比不上一维指针,或者是 Debug 版本的问题?... 让学 探究任务 图描述 掌握高斯问题的算法流程图 生掌握高斯算法 一 的规律 组内交流 完成任务 师生共同总结 一 探究任务二: 探究任务 启动 E 语言开发环境 二... 优化小结_化学_自然科学_专业资料。优化第一步:确定...(如根据晶体软件获得高斯直角坐标输入文件,软件可在...再问一个问题,如果出现虚频如何修改得 到分子稳定的... ccs3.3 错误总结_计算机软件及应用_IT/计算机_专业...[], n, m) 应改为:Matrix_DEL(float A[],...8 问题 : 在高斯牛顿迭 代是 , 欲迭代十次 , ... 第2课 利用E语言解决高斯问题_其它课程_初中教育_教育专区。《第 2 课利用 E...设计意图:总结归纳本课所学,系统掌握知识体系。 (五)评价评价项目 优秀 算法...
