高斯l502.exe报错和关于中多种基组的命名

迷茫的小飞侠 · 发表于 Post on 2016-3-16 01:22:00

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原子），三种基组的和混合定义和两种基组的混合定义是一样的吗？先表示价层电子基组再表示赝势基组。（求组论坛的各位前辈）

非常谢谢

liyuanhe211 · 发表于 Post on 2016-3-16 10:32:17

本帖最后由 liyuanhe211 于 2016-3-16 10:35 编辑

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

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

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

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

sobereva · 发表于 Post on 2016-3-16 10:42:28

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

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

迷茫的小飞侠 · 发表于 Post on 2016-3-16 17:33:50

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老师

迷茫的小飞侠 · 发表于 Post on 2016-3-16 18:04:41

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

liyuanhe211 · 发表于 Post on 2016-3-16 20:01:11

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

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

bomsaude · 发表于 Post on 2016-3-16 20:14:23

你能不能先用大核试一下？我优化过100个左右原子的稀土配合物，用的是大核，都能收敛，看到别人说，小核不好收敛，大核收敛了，用这个收敛结构看看小核能否收敛。我也不太懂。

迷茫的小飞侠 · 发表于 Post on 2016-3-16 20:55:02

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

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

非常谢谢

迷茫的小飞侠 · 发表于 Post on 2016-3-16 20:57:46

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

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

liyuanhe211 · 发表于 Post on 2016-3-16 21:03:14

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

非常谢谢

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

liyuanhe211 · 发表于 Post on 2016-3-16 21:03:48

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

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

迷茫的小飞侠 · 发表于 Post on 2016-3-16 21:09:31

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

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

迷茫的小飞侠 · 发表于 Post on 2016-3-16 21:26:46

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

非常谢谢

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

sobereva · 发表于 Post on 2016-3-16 21:38:24

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

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

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

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

迷茫的小飞侠 · 发表于 Post on 2016-3-16 22:24:42

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老师提供的资料和解答

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] 高斯l502.exe报错和关于中多种基组的命名

浏览过的版块