计算化学公社

标题: 高斯l502.exe报错和关于中多种基组的命名 [打印本页]

作者
Author: 迷茫的小飞侠 时间: 2016-3-16 01:22
标题: 高斯l502.exe报错和关于中多种基组的命名
1.现在在做硅晶胞晶体中掺杂稀土金属的测试（主要是f区金属元素），在优化中一直出现l502.exe的报错而且并没有出现正常的优化结构（收敛部分基本是没有出现的），为此使用了各种解决方案，依旧没有解决.
使用的关键词是p B3lyp/Genecp NoSym Int=UltraFine SCF(dsymm,MaxCyc=400) punch(Coord) Opt(Cartesian,gdiis)，基组是根据近期的文献使用了ECP28MWB for La，6-31G* for Si H
下面内容是出现的502报错

由于相关的解决措施已经尝试，现在依旧不能得到优化结果，渴求各位前辈指导！

2. 给予以上优化我想采用混合基组（两种基组的混合定义是会的），我能否对La和Si同时进行赝势基组的定义？因为适合f区元素的基组较少且文献中采用stuttgart RSC 1997 ECP，而Si则可以使用常规的SDD进行定义. 在测试中是采用同时对La和Si都适用的ECP基组（CRENBL ECP）还是分别定义两种赝势基组：ECP28MWB for La+SDD for Si. 如果是分别定义基组，在输入结构中就需要对三种混合基组进行定义（包括用于饱和的H原子），三种基组的和混合定义和两种基组的混合定义是一样的吗？先表示价层电子基组再表示赝势基组。（求组论坛的各位前辈）

非常谢谢

作者
Author: liyuanhe211 时间: 2016-3-16 10:32
本帖最后由 liyuanhe211 于 2016-3-16 10:35 编辑

"相关的解决措施"都有什么？
观察SCF收敛趋势，除非看起来有收敛趋势，SCF=maxcycle才有用，有用的时候十分罕见。

阅读 http://sobereva.com/61 了？本论坛一搜SCF就有，如果使用了，估计就不会贴出这样的关键词。

SCF不收敛为什么要调opt的选项

Si才第三行。。。SDD就是D95哪来的赝势

作者
Author: sobereva 时间: 2016-3-16 10:42
1 多余的垃圾关键词都去掉，包括
NoSym SCF(dsymm,MaxCyc=400) punch(Coord)
不收敛看具体输出文件，分析收敛趋势。有可能解决不收敛的方法无非就是这些：
解决SCF不收敛问题的方法
http://sobereva.com/61

2 对Si不建议用混合基组，6-31G*正合适。CRENBL不建议用。稀土用SDD，其它6-31G*足矣，不用考虑其它的。

作者
Author: 迷茫的小飞侠 时间: 2016-3-16 17:33

sobereva 发表于 2016-3-16 10:42
1 多余的垃圾关键词都去掉，包括
NoSym SCF(dsymm,MaxCyc=400) punch(Coord)
不收敛看具体输出文件，分析 ...

谢谢Sob老师，我对Si使用的是6-31G*，稀土用的SDD，泛函使用了常用的B3LYP.我做了一个相关的测试，把Si晶胞中掺杂的稀土金属改成了Ni，然后使用了相同的泛函基组以及关键词，其优化是能正常进行的，
相反对于喊稀土金属的优化，我使用之前小木虫上一个关于scf与几何优化不收敛帖子里面提到的方法：
1）增大scf循环次数
2）使用scf=dsymm
3）对结构进行微调
4）scf中添加QC（以前用过一次耗时太长们所以没有采用）
5）int=ultrafine

关键点：优化中当迭代圈数达到400以后优化就终止了，且没有出现迭代能量的曲线。基于这种情况，我是不是还不要继续增大迭代的圈数？
我看了Sob提供的信息，我想尝试同时采用SCF=NoVarAcc SCF=vshift=x 这两个关键词，不知道这两个使用是否会冲突（因为Sob老师在解决方法中写道“大多数可以结合在一起使用以得到更好的效果”）

谢谢Sob老师

作者
Author: 迷茫的小飞侠 时间: 2016-3-16 18:04

liyuanhe211 发表于 2016-3-16 10:32
"相关的解决措施"都有什么？
观察SCF收敛趋势，除非看起来有收敛趋势，SCF=maxcycle才有用，有用的时候十 ...

不好意思我的表达可能有些错误。
我想说的是：SDD是可以应用到除Fr和Ra之外的原子（高斯手册中有提到“SDD的使用范围除F如何Ra之外的元素”.），Stuttgart/Dresden ECP指定的基组和芯势. 这里的Stuttgart/Dresden ECP是属于与SDD关键词等价的基组，对于到Cl的元素选择使用了全电子基组，之后的元素使用ECP. 基组高斯手册中我的理解是SDD是可以应用到大部分元素的，只是对于Ar之后的元素需要采用SDD对应的赝势. （我不知道这样的理解是否有误，如果不对请您指导，谢谢）

对于我提到的502报错，我使用了以下的方法：
1）增大scf循环次数
2）使用scf=dsymm
3）对结构进行微调
4）scf中添加QC（以前用过一次耗时太长们所以没有采用）
5）int=ultrafine

作者
Author: liyuanhe211 时间: 2016-3-16 20:01

迷茫的小飞侠发表于 2016-3-16 18:04
不好意思我的表达可能有些错误。
我想说的是：SDD是可以应用到除Fr和Ra之外的原子（高斯手册中有提到“ ...

我知道SDD技术上可以用，但是不适合。不用管小木虫，换成合适的基组，看Sob和我贴的链接。

作者
Author: bomsaude 时间: 2016-3-16 20:14
你能不能先用大核试一下？我优化过100个左右原子的稀土配合物，用的是大核，都能收敛，看到别人说，小核不好收敛，大核收敛了，用这个收敛结构看看小核能否收敛。我也不太懂。

作者
Author: 迷茫的小飞侠 时间: 2016-3-16 20:55

bomsaude 发表于 2016-3-16 20:14
你能不能先用大核试一下？我优化过100个左右原子的稀土配合物，用的是大核，都能收敛，看到别人说，小核不 ...

什么叫大核？具体是指什么？之前没有听到过大核和小核的概念、

非常谢谢

作者
Author: 迷茫的小飞侠 时间: 2016-3-16 20:57

liyuanhe211 发表于 2016-3-16 20:01
我知道SDD技术上可以用，但是不适合。不用管小木虫，换成合适的基组，看Sob和我贴的链接。

恩恩，已看您和Sob老师给的链接内容，谢谢！
有的解决方法之前已经使用了，最后的结果是“convergence criterion not met”且没有迭代过程的能量曲线. 所以不是很能明白出错的原因.

作者
Author: liyuanhe211 时间: 2016-3-16 21:03

迷茫的小飞侠发表于 2016-3-16 20:55
什么叫大核？具体是指什么？之前没有听到过大核和小核的概念、

非常谢谢

去了解（查找）大核赝势和小核赝势的相关知识。

作者
Author: liyuanhe211 时间: 2016-3-16 21:03

迷茫的小飞侠发表于 2016-3-16 20:57
恩恩，已看您和Sob老师给的链接内容，谢谢！
有的解决方法之前已经使用了，最后的结果是“convergence c ...

数值在输出里搜索，需要能量曲线可以用很多工具/小脚本自己画。

作者
Author: 迷茫的小飞侠 时间: 2016-3-16 21:09

liyuanhe211 发表于 2016-3-16 21:03
去了解（查找）大核赝势和小核赝势的相关知识。

好的，已找到大小核赝势的相关信息，谢谢.

作者
Author: 迷茫的小飞侠 时间: 2016-3-16 21:26

迷茫的小飞侠发表于 2016-3-16 20:55
什么叫大核？具体是指什么？之前没有听到过大核和小核的概念、

非常谢谢

已弄明白大核赝势和小核赝势的区别谢谢

作者
Author: sobereva 时间: 2016-3-16 21:38

迷茫的小飞侠发表于 2016-3-16 17:33
谢谢Sob老师，我对Si使用的是6-31G*，稀土用的SDD，泛函使用了常用的B3LYP.我做了一个相关的测试，把Si晶 ...

小木虫上看到的解决不收敛的办法一定是坑爹的，只会误人子弟白费时间
怎么解决要看我贴的博文
不会冲突

大核和小核的区别看
在赝势下做波函数分析的一些说明
http://sobereva.com/156

用了#P就会输出SCF迭代过程每一圈的能量变化，看看收敛趋势

作者
Author: 迷茫的小飞侠 时间: 2016-3-16 22:24

sobereva 发表于 2016-3-16 21:38
小木虫上看到的解决不收敛的办法一定是坑爹的，只会误人子弟白费时间
怎么解决要看我贴的博文
不会冲 ...

谢谢sob老师。
我在模拟的时候使用了#P，但是在Gauss View中并没有得到迭代的能量曲线图。根据Sob写的博文我想要同时使用SCF=NoVarAcc SCF=vshift=x但是结果中显示convergence criterion not met“，所以我在想我是不是还需要修改收敛限或者继续增加scf圈数。但是这里我又有了一个新的问题没在Sob老师的博文中提到”做几何优化，振东分析的时候不建议降低默认的SCF收敛限，否则可能结果不准确，还可能阻碍几何优化的收敛.“
重新编辑的关键词如下：
#p B3lyp/Genecp Int=UltraFine SCF(NoVarAcc,vshift=300) Opt （这样是否可行，我不是很确定）

谢谢Sob老师提供的资料和解答

作者
Author: bomsaude 时间: 2016-3-16 22:41

迷茫的小飞侠发表于 2016-3-16 14:24
谢谢sob老师。
我在模拟的时候使用了#P，但是在Gauss View中并没有得到迭代的能量曲线图。根据Sob写的博 ...

不管是优化正常结束的还是正在优化的输出文件，gaussview在打开文件时，界面的下面有一个选项（默认是不选的），你选上那个就有能量变化曲线了。

作者
Author: sobereva 时间: 2016-3-16 22:57

迷茫的小飞侠发表于 2016-3-16 22:24
谢谢sob老师。
我在模拟的时候使用了#P，但是在Gauss View中并没有得到迭代的能量曲线图。根据Sob写的博 ...

gview是显示不了迭代过程中能量曲线图的。自行看输出的SCF每一圈的能量，或者用grep命令批量提取以便观看
scf圈数你增加到100000000000000次也照样收敛不了，迭代100年也收敛不了。
修改收敛限到SCF=conver=6对单点可以接受，对几优化不要修改
没有万能的解决不收敛的办法，必须原理上搞懂博文中提到的不同解决方式，反复尝试。

作者
Author: 迷茫的小飞侠 时间: 2016-3-16 23:01

bomsaude 发表于 2016-3-16 22:41
不管是优化正常结束的还是正在优化的输出文件，gaussview在打开文件时，界面的下面有一个选项（默认是不 ...

我知道这个，我选择了但是没有显示. 我看了输出文件，但也没有找到为什么. 一般我们做了优化以后如果选择高斯view的open到read intermediate geometries的方式打开log文件：在查看时应该是会看到迭代能量曲线以及优化的每一步对应的结果，但是现在不仅迭代能量曲线没有而且打开后的log文件也有对应的优化结果（log文件打开后座上是没有出现）优化步数的.

作者
Author: 迷茫的小飞侠 时间: 2016-3-16 23:09

sobereva 发表于 2016-3-16 22:57
gview是显示不了迭代过程中能量曲线图的。自行看输出的SCF每一圈的能量，或者用grep命令批量提取以便观看 ...

恩恩好的. 因为我对有的原理理解上总是不很确定理解也不头侧所以不是太明白Sob老师提到的有些概念. 一般只要正常优化了，我们通过Gview打开计算的log文件后同时选择窗口最下方的”read intermediate geometries“选项，在查看log文件时左上方会有相应的优化步数，而”点击result“的optimization也会出现相应的能量随优化步数增加而变化的趋势图. 此时这里的能量不是scf的迭代能量. （之前理解错了）

作者
Author: 迷茫的小飞侠 时间: 2016-3-17 00:12

sobereva 发表于 2016-3-16 22:57
gview是显示不了迭代过程中能量曲线图的。自行看输出的SCF每一圈的能量，或者用grep命令批量提取以便观看 ...

Sob老师，我查看了SCF能量，刚开始基本是在震荡，当圈数是在30-370之间，总体趋势是下降.在370-385之间出现震荡，385-400又开呈下降趋势. 根据这样的数据，我可以得到什么信息呢？通过SCF的迭代能量我能知道是什么原因引起的SCF不收敛吗？在将錒系元素换成Ni以后可以收敛并优化成功，结合Sob老师博文中信息，也就是说在这一个含有錒系元素的优化中有可能是基组不合适？或者是该结构的HOMO-LUMO能隙较小（在纯晶胞单元中得到的能隙只有0.5eV）

作者
Author: sobereva 时间: 2016-3-17 12:12

迷茫的小飞侠发表于 2016-3-17 00:12
Sob老师，我查看了SCF能量，刚开始基本是在震荡，当圈数是在30-370之间，总体趋势是下降.在370-385之间出 ...

就是震荡导致的难收敛
锕系的SDD有大核赝势和小核赝势，如果之前是小核的改成大核试试，其它轻原子也可以换基组试试，也可以换泛函试试（用HF成分高一些的）。各种方法都试试，总比增大迭代次数上限有用得多得多得多得多

作者
Author: 迷茫的小飞侠 时间: 2016-3-17 22:14

sobereva 发表于 2016-3-17 12:12
就是震荡导致的难收敛
锕系的SDD有大核赝势和小核赝势，如果之前是小核的改成大核试试，其它轻原子也可 ...

谢谢Sob老师，在不了解大核和小核赝势之前，我做过了大核和小核的测试，两个测试都没有收敛。昨天重新使用了关键词（根据Sob老师博文修改了关键词）分别换了泛函和大核赝势，还是报错了，这次报的l301.exe错误。现在正在重新修改

作者
Author: 迷茫的小飞侠 时间: 2016-4-6 00:48

sobereva 发表于 2016-3-16 21:38
小木虫上看到的解决不收敛的办法一定是坑爹的，只会误人子弟白费时间
怎么解决要看我贴的博文
不会冲 ...

Sob老师关于l502的报错，我将博文中提到的方法都使用到了自己的体系中，查看了SCF的能量我确定是能量是在在一个特别小的范围内振动，但是结果依旧出现502错误. 如果是初始结构的电子结构不是很好，会不会对优化有影响？
急需Sob老师的指导，谢谢

作者
Author: sobereva 时间: 2016-4-6 12:35

迷茫的小飞侠发表于 2016-4-6 00:48
Sob老师关于l502的报错，我将博文中提到的方法都使用到了自己的体系中，查看了SCF的能量我确定是能量是 ...

优化是优化，SCF是SCF，两个问题不要搅在一起

作者
Author: 迷茫的小飞侠 时间: 2016-4-20 01:05

sobereva 发表于 2016-4-6 12:35
优化是优化，SCF是SCF，两个问题不要搅在一起

谢谢Sob老师，我再好好看看感觉被自己弄糊涂了.

欢迎光临计算化学公社 (http://bbs.keinsci.com/)