下载文档原格式
大家好：我根据自己在计算过程中遇到的收敛性的问题，做了一点关于CASTEP计算参数调试的总结。
当然，这些只是我的个人心得体会，大家如果觉得有什么问题请提出指正，或是如果有自己独到的见解，请跟帖加上来，供大家一起学习提高。
如果计算在Max. cycle全部用完之后仍然不能收敛，要先查看一下.castep文件，如果能量比刚开始计算时有明显下降，而且震动幅度也明显减小，只是还没有降到收敛标准以下，那说明Max. cycle还不够，要增加这个数值。
但多数情况下，我们会发现，能量并没有比刚开始计算时的能量有明显下降，而是围绕着某个平均值在反复地震荡，而且震动幅度也没有减小。
那说明是参数设定的问题，此时就算再怎么增加Max. cycle的数值，计算最终还是不会收敛的。
这个时候，最好是增加energy cutoff的数值，降低K point set的数值，增加smearing width数值（当然，如果通过调试前两个数值就可以使计算收敛的话，就最好不要动这个数值）。
再有，如果计算是使用的density mixing的话，那要将mix_charge_amp，mix_spin_amp的数值降低，将mix_energy_cutoff的数值增加到 energy cutoff 数值的3到4倍，增加mix_history_length的数值。
当然，如果调试上面的参数仍然发现收敛很辛苦的话，就将mixing scheme转成All Band/EDFT,这会解决收敛问题，但是计算时间将会是使用density mixing 计算所用时间的3倍以上。
1. 参数调试过于经验化，如果楼主能结合对程序原理的认识，或许能够把每个参数调节的意义给阐述的更为清晰些；从优化算法，体系特性，物理过程等方面着手来进行分析，或许可以得到更加全面的认识（如果楼主能够结合上述内容对CASTEP收敛进行详细的分析和总结的话，我会考虑接受楼主的文章，并且发表在小木虫计算与模拟研究RCS期刊上）。
1. Fluent计算中对网格质量的要求1）网格质量参数：Skewness （不能高于0.95，最好在0.90以下；越小越好）Change in Cell-Size （也是Growth Rate，最好在1.20以内，最高不能超过1.40）Aspect Ratio （一般控制在5：1以内，边界层网格可以适当放宽）Alignment with the Flow（就是估计一下网格线与流动方向是否一致，要求尽量一致，以减少假扩散）2）网格质量对于计算收敛的影响：高Skewness的单元对计算收敛影响很大，很多时候计算发散的原因就是网格中的仅仅几个高Skewness的单元。
举个例子：共有112,000个单元，仅有7个单元的Skewness超过了0.95，在进行到73步迭代时计算就发散了！高长宽比的单元使离散方程刚性增加，使迭代收敛减慢，甚至困难。
也就是说，Aspect Ratio 尽量控制在推荐值之内。
3）网格质量对精度的影响：相邻网格单元尺寸变化较大，会大大降低计算精度，这也是为什么高连续方程残差的原因。
网格线与流动是否一致也会影响计算精度。
4）网格单元形状的影响：非结构网格比结构网格的截断误差大，因此，为提高计算精度计，请大家尽量使用结构网格，对于复杂几何，在近壁这些对流动影响较大的地方尽量使用结构网格，在其他次要区域使用非结构网格。
2. 关于fluent模拟的论文不好发发CFD专业期刊估计不行，但发一些应用型的期刊应该可以，见过不少。
数值与实验结果对比是必不可少的，如果没有实验对比，就难了。
好不好发看你模拟的什么？模拟的深度！简单的两维的常温常压的空气流场模拟应该不好发！你要是模拟的3维的模拟国际前沿问题，考虑了别人没有考虑的东西，好发的很！Fluent做出的东西,一来看,有没有相关的实验数据作对比,或者是验证你做的模型.二来,你自己做理论模型的研究,加入到Fluent中去;还是有希望的.一般的比较好发，我们实验室全是做模拟的，就发有个叫计算机与应用化学的。
【分享】高斯对不收敛问题的对策首先，我们必须理解收敛是什么意思。
在自洽场(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是迭代数)。
这很少有帮助，但值得一试。
3 放宽收敛标准如果接近SCF但未达到，收敛标准就会放松或者忽略收敛标准。
这通常用于不是在初始猜测而是在平衡结构收敛的几何优化。
SCF=Sleazy放松收敛标准，Conver选项给出更多的控制。
4 尝试改变初始构型首先略微减小键长，接下来略微增加键长，接下来再对结构作一点改变。
5 尝试能级移动Level shifting (SCF=Vshift)如果不收敛的原因是波函数的震荡行为，通常是因为在相近的能量上的泰的混合。
对于这种情况，我们可以采用level shifting的方法。
Level shifting的含义是人工的升高非占据轨道的能级，以防止和最高占据轨道之间的混合，以达到收敛的目的。
化学计算中帮助几何优化收敛的常用方法<来自小木虫>文/Sobereva First release: 2012-Oct-13几何优化，也就是寻找势能面极小点结构的过程。
量子化学计算中几何优化不收敛是个老生常谈的问题，在各种论坛里、群里都已经反复讨论过很多遍了，但是还是时常看到有人问，而且现有的讨论也都不怎么全面，所以觉得有必要撰文谈一下。
所谓几何优化不收敛，也就是始终，或者很难达到收敛要求。
通常会伴随着震荡行为，即受力、几何结构变化随优化步数呈现周期性趋势。
解决这种问题必须在结合经验和理论知识的前提下，通过考察实际收敛的趋势，尝试各种可能奏效处理办法。
本文列举一些常用的解决不收敛，也包括加速收敛的办法。
其中很多方法可以相互结合使用以达到更好的效果。
这里假定用户是用Gaussian，很多方法在其它程序中也可以类似地使用。
先说一下收敛标准。
Gaussian中判断几何优化收敛有四个标准，在默认收敛设定下，这四个标准是：最大受力<0.00045；方均根受力<0.00030；最大位移<0.00180；方均根位移<0.00120当这四个标准都满足了，达成四个YES，就宣告收敛。
另外，优化过程中只要受力小于预定的收敛限100倍，哪怕位移还没低于收敛限，则也算作已收敛。
这主要考虑到势能面非常非常缓的大的柔性分子，相对于这样尺度的分子，几何结构收敛到那么精确意义不大，放宽位移收敛限避免了收敛太慢。
有时候优化出错，不是因为几何收敛问题，而是因为每一步优化中连能量计算都没能完成。
优化也可能朝着明显错误的方向进行而导致难以收敛，这极有可能是理论方法、基组、电子态及其它诸多选项的设定不合理。
这些方面和优化不收敛问题本身没关系，所以不会在本文提到。
1 尝试不同的优化方法优化几何结构的方法有很多，以前我在《过渡态、反应路径的计算方法及相关问题》（）当中详细介绍过的很多搜索过渡态的方法其实和搜索势能面极小点（即几何优化）的方法本质是一致的。
SCF不收敛怎么办？几何优化不收敛怎么办？虽然二者有联系，但实际上是两个相对独立的问题，要分开回答。
这里提供的只是解决方案，所有的方法都用上了还是不收敛得情况也是有的。
首先要分清SCF不收敛和几何构型优化不收敛：SCF不收敛指的是自洽场叠代不收敛，可以认为是对指定结构的波函数不断优化的过程，是为了找到这个某个指定结构下能量最低的波函数，而几何构型优化是对结构的优化的过程，是为了找到某个指定的组分下能量极小结构（注意，不一定是能量最小结构）。
在量子化学计算的几何构型优化中，每一步的几何构型优化都包含的很多次的SCF计算。
1、SCF不收敛的解决方案。
(1) 可以加大SCF的循环次数，默认的循环次数是128次，通过SCF=(maxcycle=n)来设置最大循环次数n。
建议不要超过512，更多的循换没有必要。
(2) 如果加大循环次数不管用，在分子有对称性的情况下，使用SCF=dsymm(在每一SCF 迭代中，进行密度矩阵的对称化)关键词来强制密度对称，有时可以收敛。
另外，此关键词很多时候对"SCF is confused”这种错误很管用。
(3) 使用SCF=symm关键词，使用的前提同上，有时可以收敛。
(4) 如果(2)(3)两步都不行，可以将对称的分子中的某几个原子的位置微调，使分子丧失对称性。
这等效于nosymm关键词，但个人经验，这种方式比nosymm好用的多。
(5) 如果还不行，只能拿出杀手锏了，就是使用QC(要求使用二次收敛SCF程序)，但不建议直接使用，而是使用XQC(对于一级SCF 尚未收敛的情况，再加入额外的SCF=QC 步骤。
)关键词，比如SCF=(maxcycle=80,XQC)，意思是如果SCF正常计算在80个循环之内不收敛才进行昂贵的QC计算，因为SCF不收敛多数在几个优化的过程中出现，无法判断哪一步优化的时候会出现SCF不收敛，所以用XQC比纯粹使用qc要省时的多。
(6) 中级用户可以在输入文件的井号“#”开头那一行井号后面加上字母"p"来输出更多的信息，其中就有自洽场叠代的信息，分析原因可能会对采用什么方法提供指导。
虽然二者有联系，但实际上是两个相对独立的问题，要分开回答。
这里提供的只是解决方案，所有的方法都用上了还是不收敛得情况也是有的。
首先要分清scf不收敛和几何构型优化不收敛：scf不收敛指的是自洽场叠代不收敛（什么？没听说过什么叫自洽场？那还是回去学习些量化基础知识再开展计算吧！），可以认为是对指定结构的波函数不断优化的过程，是为了找到这个某个指定结构下能量最低的波函数，而几何构型优化是对结构的优化的过程，是为了找到某个指定的组分下能量极小结构（注意，不一定是能量最小结构）。
在量子化学计算的几何构型优化中，每一步的几何构型优化都包含的很多次的scf计算。
1、scf不收敛的解决方案。
(1) 可以加大scf的循环次数，默认的循环次数是128次，通过scf=(maxcycle=n)来设置最大循环次数n。
建议不要超过512，更多的循换没有必要。
(2) 如果加大循环次数不管用，在分子有对称性的情况下，使用scf=dsymm关键词来强制密度对称，有时可以收敛。
另外，此关键词很多时候对"scf is confused”这种错误很管用。
(3) 使用scf=symm关键词，使用的前提同上，有时可以收敛。
(4) 如果(2)(3)两步都不行，可以将对称的分子中的某几个原子的位置微调，使分子丧失对称性。
这等效于nosymm关键词，但个人经验，这种方式比nosymm好用的多。
(5) 如果还不行，只能拿出杀手锏了，就是使用qc，但不建议直接使用，而是使用xqc关键词，比如scf=(maxcycle=80,xqc)，意思是如果scf正常计算(dc)在80个循环之内不收敛才进行昂贵的qc计算，因为scf不收敛多数在几个优化的过程中出现，无法判断哪一步优化的时候会出现scf不收敛，所以用xqc比纯粹使用qc要省时的多。
(6) 中级用户可以在输入文件的井号“#”开头那一行井号后面加上字母"p"来输出更多的信息，其中就有自洽场叠代的信息，分析原因可能会对采用什么方法提供指导。
量子化学计算中帮助几何优化收敛的常用方法几何优化，也就是寻找势能面极小点结构的过程。
量子化学计算中几何优化不收敛是个老生常谈的问题，在各种论坛里、群里都已经反复讨论过很多遍了，但是还是时常看到有人问，而且现有的讨论也都不怎么全面，所以觉得有必要撰文谈一下。
所谓几何优化不收敛，也就是始终，或者很难达到收敛要求。
通常会伴随着震荡行为，即受力、几何结构变化随优化步数呈现周期性趋势。
解决这种问题必须在结合经验和理论知识的前提下，通过考察实际收敛的趋势，尝试各种可能奏效处理办法。
本文列举一些常用的解决不收敛，也包括加速收敛的办法。
其中很多方法可以相互结合使用以达到更好的效果。
这里假定用户是用Gaussian，很多方法在其它程序中也可以类似地使用。
先说一下收敛标准。
Gaussian中判断几何优化收敛有四个标准，在默认收敛设定下，这四个标准是：最大受力<0.00045；方均根受力<0.00030；最大位移<0.00180；方均根位移<0.00120 当这四个标准都满足了，达成四个YES，就宣告收敛。
另外，优化过程中只要受力小于预定的收敛限100倍，哪怕位移还没低于收敛限，则也算作已收敛。
这主要考虑到势能面非常非常缓的大的柔性分子，相对于这样尺度的分子，几何结构收敛到那么精确意义不大，放宽位移收敛限避免了收敛太慢。
有时候优化出错，不是因为几何收敛问题，而是因为每一步优化中连能量计算都没能完成。
优化也可能朝着明显错误的方向进行而导致难以收敛，这极有可能是理论方法、基组、电子态及其它诸多选项的设定不合理。
这些方面和优化不收敛问题本身没关系，所以不会在本文提到。
1 尝试不同的优化方法优化几何结构的方法有很多，以前我在《过渡态、反应路径的计算方法及相关问题》（/sobereva/item/c5a8de0b88de1c10acdc70d7）当中详细介绍过的很多搜索过渡态的方法其实和搜索势能面极小点（即几何优化）的方法本质是一致的。
GaussianFAQ（不断更新中）新问题可以到/bbs/viewthread.php?tid=1734563跟帖GaussianFAQ改进历程2009.09.25 V 0.01 GaussianFAQ(Beta版)诞生2009.09.30 V 1.00 为这个FAQ增加了插图11张，GaussianFAQ第一版正式发布2009.10.02 V 1.10 修改了第一章第四节的回答2009.10.04 V 2.00 增加了第二章第七、八、九节三个问题及其答案2009.10.06 V 2.10 修改了第一章第三节的回答2009.10.07 V 3.00 增加了第二章第十、十一节两个问题及其答案2009.10.07 V 4.00 更换了第一章第一节这个问题的答案2009.10.08 V 4.10 修改了第二章第九节的回答2009.10.08 V 5.00 增加了第二章第十二节这个问题及其答案2009.10.10 V 6.00 增加了第一章第五节这个问题及其答案2009.10.14 V 6.01 改了几个错别字2009.10.15 V 6.02 又改了几个错别字2009.10.18 V 7.00 增加第一章第六节这个问题及其答案2009.10.19 V 8.00 增加第一章第七节这个问题及其答案2009.10.20 V 9.00 增加第二章第十三节这个问题及其答案2009.11.27 V 9.10 修改第一章第四节的答案2009.12.08 V10.00 增加了第二章第十四、十五节两个问题及其答案2009.12.09 V11.00 增加了第四章第六节这个问题及其答案2009.12.12 V12.00 增加了第四章第七节这个问题及其答案2009.12.15 V13.00 增加了第四章第八节这个问题及其答案2009.12.15 V13.10 修改第一章第二节的答案2009.12.21 V13.20 修改第一章第四节的答案2009.12.28 V14.00 增加了第一章第八节这个问题及其答案2009.12.31 V14.10 修改了第四章第三节的答案2010.01.05 V14.20 修改了第一章第四节的答案2010.01.06 V14.21 修改了第一章第四节的答案中的一个bug2010.01.07 V15.00 增加了第四章第九节这个问题及其答案2010.01.08 V15.01 将原第一章第七节调整为第二章第十六节2010.01.09 V16.00 更换了第二章第七节的答案2010.01.11 V17.00 增加了第二章第十七节这个问题及其答案2010.01.11 V18.00 更换了第一章第五节的答案2010.01.12 V18.10 修改了第二章第九节的答案2010.01.13 V19.00 增加了第二章第十八节这个问题及其答案2010.01.13 V19.10 修改了第二章第十七节的答案2010.01.13 V20.00 增加了第二章第十九节这个问题及答案2010.01.14 V21.00 增加了第二章第二十节这个问题及答案2010.01.20 V21.01 修改了第四章第九节的一个bug2010.01.22 V21.10 修改了第四章第九节的答案2010.01.27 V22.00 增加了第四章第十节这个问题及答案2010.02.18 V23.00 增加了第二章第二十一节这个问题及答案2010.02.19 V24.00 增加了第二章第二十二节这个问题及答案2010.02.19 V25.00 增加了了目录2010.04.02 V26.00 增加了第三章第五节这个问题及答案2010.04.06 V26.10 修改了第四章第九节的答案2010.05.26 V26.11 略微修改了第四章第九节的答案2010.05.27 V26.12 略微修改了第一章第四节的答案2010.07.25 V27.00 增加了第二章第二十三节这个问题及答案2010.08.21 V27.01 略微修改了第二章第二十三节的答案请大家继续提意见版本号升级说明：GaussianFAQ的版本号采用x.xx的方式彻底更换某个问题的答案，或者增加新问题及其解答，大版本号加1修改对某个问题的回答，小版本号加1改点儿bug，小小版本号加1该帖子主要讨论Gaussian操作问题。
GaussianFAQ实践版（V25.00）【兼在线答疑】作者: yjcmwgk(站内联系TA)发布: 2009-09-25 GaussianFAQ（不断更新中）新问题可以到/bbs/viewthread.php?tid=1734563跟帖GaussianFAQ改进历程2009.09.25 V 0.01 GaussianFAQ(Beta版)诞生2009.09.30 V 1.00 为这个FAQ增加了插图11张，GaussianFAQ第一版正式发布2009.10.02 V 1.10 修改了第一章第四节的回答2009.10.04 V 2.00 增加了第二章第七、八、九节三个问题及其答案2009.10.06 V 2.10 修改了第一章第三节的回答2009.10.07 V 3.00 增加了第二章第十、十一节两个问题及其答案2009.10.07 V 4.00 更换了第一章第一节这个问题的答案2009.10.08 V 4.10 修改了第二章第九节的回答2009.10.08 V 5.00 增加了第二章第十二节这个问题及其答案2009.10.10 V 6.00 增加了第一章第五节这个问题及其答案2009.10.14 V 6.01 改了几个错别字2009.10.15 V 6.02 又改了几个错别字2009.10.18 V 7.00 增加第一章第六节这个问题及其答案2009.10.19 V 8.00 增加第一章第七节这个问题及其答案2009.10.20 V 9.00 增加第二章第十三节这个问题及其答案2009.11.27 V 9.10 修改第一章第四节的答案2009.12.08 V10.00 增加了第二章第十四、十五节两个问题及其答案2009.12.09 V11.00 增加了第四章第六节这个问题及其答案2009.12.12 V12.00 增加了第四章第七节这个问题及其答案2009.12.15 V13.00 增加了第四章第八节这个问题及其答案2009.12.15 V13.10 修改第一章第二节的答案2009.12.21 V13.20 修改第一章第四节的答案2009.12.28 V14.00 增加了第一章第八节这个问题及其答案2009.12.31 V14.10 修改了第四章第三节的答案2010.01.05 V14.20 修改了第一章第四节的答案2010.01.06 V14.21 修改了第一章第四节的答案中的一个bug 2010.01.07 V15.00 增加了第四章第九节这个问题及其答案2010.01.08 V15.01 将原第一章第七节调整为第二章第十六节2010.01.09 V16.00 更换了第二章第七节的答案2010.01.11 V17.00 增加了第二章第十七节这个问题及其答案2010.01.11 V18.00 更换了第一章第五节的答案2010.01.12 V18.10 修改了第二章第九节的答案2010.01.13 V19.00 增加了第二章第十八节这个问题及其答案2010.01.13 V19.10 修改了第二章第十七节的答案2010.01.13 V20.00 增加了第二章第十九节这个问题及答案2010.01.14 V21.00 增加了第二章第二十节这个问题及答案2010.01.20 V21.01 修改了第四章第九节的一个bug2010.01.22 V21.10 修改了第四章第九节的答案2010.01.27 V22.00 增加了第四章第十节这个问题及答案2010.02.18 V23.00 增加了第二章第二十一节这个问题及答案2010.02.19 V24.00 增加了第二章第二十二节这个问题及答案2010.02.19 V25.00 增加了了目录2010.04.02 V26.00 增加了第三章第五节这个问题及答案2010.04.06 V26.10 修改了第四章第九节的答案2010.05.26 V26.11 略微修改了第四章第九节的答案2010.05.27 V26.12 略微修改了第一章第四节的答案2010.07.25 V27.00 增加了第二章第二十三节这个问题及答案2010.08.21 V27.01 略微修改了第二章第二十三节的答案请大家继续提意见版本号升级说明：GaussianFAQ的版本号采用x.xx的方式彻底更换某个问题的答案，或者增加新问题及其解答，大版本号加1 修改对某个问题的回答，小版本号加1改点儿bug，小小版本号加1该帖子主要讨论Gaussian操作问题。
首先要分清scf不收敛和几何构型优化不收敛：scf不收敛指的是自洽场叠代不收敛（什么？没听说过什么叫自洽场？那还是回去学习些量化基础知识再开展计算吧），可以认为是对指定结构的波函数不断优化的过程，是为了找到这个某个指定结构下能量最低的波函数，而几何构型优化是对结构的优化的过程，是为了找到某个指定的组分下能量极小结构（注意，不一定是能量最小结构）。
在量子化学计算的几何构型优化中，每一步的几何构型优化都包含的很多次的scf计算。
1、scf不收敛的解决方案。
(1) 可以加大scf的循环次数，默认的循环次数是128次，通过scf=(maxcycle=n)来设置最大循环次数n。
建议不要超过512，更多的循换没有必要。
(2) 如果加大循环次数不管用，在分子有对称性的情况下，使用scf=dsymm关键词来强制密度对称，有时可以收敛。
另外，此关键词很多时候对"scf is confused”这种错误很管用。
(3) 使用scf=symm关键词，使用的前提同上，有时可以收敛。
(4) 如果(2)(3)两步都不行，可以将对称的分子中的某几个原子的位置微调，使分子丧失对称性。
这等效于nosymm关键词，但个人经验，这种方式比nosymm好用的多。
(5) 如果还不行，只能拿出杀手锏了，就是使用qc，但不建议直接使用，而是使用xqc关键词，比如scf=(maxcycle=80,xqc)，意思是如果scf正常计算(dc)在80个循环之内不收敛才进行昂贵的qc计算，因为scf不收敛多数在几个优化的过程中出现，无法判断哪一步优化的时候会出现scf不收敛，所以用xqc比纯粹使用qc要省时的多。
(6) 中级用户可以在输入文件的井号“#”开头那一行井号后面加上字母"p"来输出更多的信息，其中就有自洽场叠代的信息，分析原因可能会对采用什么方法提供指导。
(7) 前面有虫子提到一个有用的方案但没说清楚，我这里补充一下：如果用用小基组计算，scf可以收敛，那么保存好检查点文件，换成大基组的时候从检查点文件中读取初始猜测(使用guess=read关键词)，有时可以算过去。
(8) 如果你计算中使用了含有弥散函数的基组（比如6-31+G*基组中的“+”，Aug-cc-pVTZ基组中的Aug等），可以先去掉弥散函数（比如例子中的两个基组变为6-31G*和cc-pVTZ），这样很容易收敛，需要的话，换成带有弥散函数的基组时，从检查点文件的读初始猜测。
原理与(7)类似，有时候管用。
注意：(7)和(8)不适合很大的体系，此时，使用小基组或者不带弥散函数的基组与后来使用的大基组在基组数量上差得太多（如果你本来就遇到了scf收敛问题，而使用这两个方案时，机组数量差距超过512个，多数时候要出问题），反倒容易引起scf不收敛，因为读检查点文件时，对那些多出来的基组的猜测是空的，而直接计算时最起码还有个半经验的计算提供了相对合理的初始猜测。
(9) 上述方法有时可以组合使用，如果经过各种组合处理都不收敛，那么放弃吧，你的分子的电子结构太差了。
2、几何构型优化不收敛的解决方案。
引起这个问题的原因比较多，可能也说不全，尽量做到全面。
(1) 如果是很小的分子（10原子以内），初始结构可能离平衡结构比较远，又在输出文件优化的最后一步判断是否收敛的位置（可以通过查找Threshold字段找到，它下面有四个判断项，都是YES才代表优化收敛）看到了“-- Number of steps exceeded, NStep=xxx”，可以通过加大优化的循环次数来解决问题。
使用关键词opt=(maxcycle=n).
(2) 优化到如果四个收敛标准中前两项早就收敛了，而后两项尤其是第三项不收敛，这需要判断原因。
如果每一步优化过程中的能量始终在下降，那么可以继续让他算，超过了最大步数停掉了的话把结构拿出来，重新提交就行了。
如果能量忽大忽小出现跳跃，说明遇到了比较平坦的势能面，或者优化的算法不好。
此时建议按如下顺序处理：
a. 使用opt=maxstep=n来缩小最大步长，即原子移动的距离，n的默认值是30，只能设置整
数。
b. 使用关键词opt=gdiis，很多时候可以很快收敛。
c. 使用opt=calcfc关键词，在几何优化的第一步计算力常数，为优化定一个大方向，类似在第一个三岔路口选一条路。
d. 使用opt=calcall关键词，在几何优化计算的每一步都计算力常数，类似于在每一个三岔路口都来选一条路。
这个关键词不到万不得已不要使用，非常耗时。
(3) 四个判断标准的后两项早就收敛了，但前两项似乎固定在某个数值上了。
遇到这种情况，首先去检查你优化第一步的时候给出的每个轨道得本征值部分，也就是轨道的能量，如果HOMO和LUMO的能量相同或者非常接近，那么你的设置有问题，即自旋多重度不对。
如果HOMO-LUMO能量没问题，那么将目前的计算停掉，提出结构，使用opt=gdiis关键词，很多时候管用。
(4) 如果是在几何构型优化的过程中出现因为scf不收敛导致无法优化，可以先用上面scf 不收敛的解决方案来处理，如果都不行，可以通过在#后面添加p输出更多信息，如果看到每次scf能量都能收敛到10的-6次方以下然后开始出现能量的跳跃，就是收敛不到10的-8次方（这种情况还比较常见），可以使用scf=(conver=6)来让几何优化继续下去，等优化好了，再将此关键词去掉，我的经验，优化的结构一般不会有太大的差异。
如果不到10的-6能量就开始跳越，不适合此方法。
(5) 如果几何优化过程中发现键变得异常短，而且此时几何优化总是不收敛。
遇到这种情况，问题出在你的基组使用，肯定是使用了赝势基组却没有读入赝势。
由于赝势基组将内层电子和原子核用一个有效核心势来表示，你没有读入赝势就说明没有读入这个有效核心势，那么相当于把原子的内层电子都给除掉了，巨大的核引力把原子之间的距离变得异常短也不足为奇。
总结
赠送以下资料
总结会讲话稿
在这金秋十月，丹桂飘香的季节，当我们国人还沉浸在“两个奥运，同样精彩”，“神七成功问天”的喜悦中，我们又一次迎来了九月份“6s”总结暨颁奖大会。
首先我代表“6s”小组感谢各位管理和全体同仁在“6s”运动中付出的汗水和努力！正是有了你们的积极参与，“6s”才得如火如荼，生机盎然，正是有了你们的不懈支持，我们对“6s”运动充满了信心！
9月份参加“6s”考评的有20个单位，以车间为单位的成型一车间在9月份“6s”考评中荣获“团体第一名”，主任陈新华；以线（组）为单位的成型一车间包边线荣获线（组）第一名，线长罗润生；成型一车间c线荣获线（组）第二名，线长龙勇；成型二车间e线荣获线（组）第三名，线长郑权。
我们用最热烈的掌声祝贺以上获奖单位！最值得一提的是成型一车间和包边线成功卫冕“团体第一名”“线（组）第一名”的殊荣，都说“江山易攻不易守”，在考评分数悬殊不
是很大的情况下，能再次登上这光荣的宝座，我们再次把掌声献给他们，谢谢他们对“6s”运动的大力支持。
在9月份，公司进入棉鞋生产高峰期，各位管理人员和员工忙于新品试样与生产，不免对“6s”运动产生疏忽，造成考评成绩整体不是很高，在这里我希望在9月份“6s”运动考评中成绩不是很理想的线（组）、车间，不要气馁，要打起精神，鼓起勇气，不能被一时的落后消磨了你们的锐气，在“6s”面前我们戈美其人是不会低头的！
借这个机会，我觉得还是有必要把“6s”运动的意义跟在座的各位讲讲：“6s”运动是一项全方位、全过程、全员参与的系统工程，是一项长期的工作，是与生产相辅相成的活动，它既能服务于生产，也能监督生产，想必这一点，各位已体会至深，我们要把“整理、整顿、清扫、清洁、素养、安全”融入到我们工作的每一个细节中。
为了改善和增加作业面积，做到现场无杂物，安全通道畅通，从而提高工作效率，消除管理上的混放、混料等差错事故。
我们就应该彻底地将要与不要的东西归类，并将不要的东西加以处理，需对“留之无用，弃之可惜”的观念予以突破。
为了最终能实现目视管理，把经过整理出来的需要的人、事、物要加以定量、定位、定品。
为了让我们有一个舒适、安全的工作环境，就得要求全
员参与，排除隐患，重视预防，彻底地将自己的工作环境四周打扫干净，设备异常时马上维修，以维护人身与财产安全，实现一个零故障，无意外事故发生的工作场所。
千万不要因小失大，对操作人员的操作技能进行训练，要持证上岗，让“勿以善小而不为，勿以恶小而为之”的理念深入到每一位员工的思想中。
为了使现场保持完美和最佳状态，我们需要秉持三个观念：（1）只有在整齐、干净的工作场所才能产生出高效率、高品质的产品；（2）6s是一种用心的行为，千万不要只在表面下功夫；（3）6s是一种随时随地的工作，而不是上下班前后的工作。
最后只有不断努力提高人员的素养，养成严格遵守规章制度的习惯和作风，我们的“6s”活动就能顺利开展，就能坚持下来。
谢谢大家！。
