计算化学公社

标题: 含镧元素体系SCF不收敛，是镧元素本身的原因吗？ [打印本页]

作者
Author: CCCC 时间: 2018-3-5 15:30
标题: 含镧元素体系SCF不收敛，是镧元素本身的原因吗？
%chk=La2_C72.chk
%mem=6000MB
%nproc=2
#p B3LYP/Gen pseudo=cards Freq opt SCF=XQC optcyc=100 scfcyc=128
碳原子用的是6-311G* 镧元素用的是LANL2DZ
这是输入文件关键词部分，请问是镧元素本身的原因还是其他的呢？非常感谢！

作者
Author: sobereva 时间: 2018-3-5 15:43
不要用pseudo=cards
怎么使用混合基组和赝势看
详解Gaussian中混合基组、自定义基组和赝势基组的输入
http://sobereva.com/60

垃圾关键词太多，看
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331

对镧系不要用lanl2DZ，改用SDD，结果明显更好，而且lanl2DZ和6-311G*在档次上根本不搭
谈谈赝势基组的选用
http://sobereva.com/373

以上问题都解决后，如果还不收敛到时候再说

作者
Author: 我本是个娃娃 时间: 2018-3-5 17:00
另外有过渡金属尽量选择纯泛函，比如M06L这样的能够处理静态相关强的体系，B3LYP有的时候也不是特别的好

作者
Author: sobereva 时间: 2018-3-5 18:36

我本是个娃娃发表于 2018-3-5 17:00
另外有过渡金属尽量选择纯泛函，比如M06L这样的能够处理静态相关强的体系，B3LYP有的时候也不是特别的好

但纯泛函代价是往往更难收敛。
鉴于目前SCF收敛还没解决，最好先在B3LYP层面上解决了再考虑其它泛函

作者
Author: 我本是个娃娃 时间: 2018-3-5 21:25

sobereva 发表于 2018-3-5 18:36
但纯泛函代价是往往更难收敛。
鉴于目前SCF收敛还没解决，最好先在B3LYP层面上解决了再考虑其它泛函

了解了，谢谢大博士

作者
Author: helpme 时间: 2018-3-6 08:48
算过镧系的Ce，确实难收敛。即使收敛一般也特别慢。

作者
Author: CCCC 时间: 2018-3-6 09:07

helpme 发表于 2018-3-6 08:48
算过镧系的Ce，确实难收敛。即使收敛一般也特别慢。

方便告诉我你用的哪个关键词吗？

作者
Author: helpme 时间: 2018-3-6 18:33
没有用什么特别的，就是SCF不收敛的那一套方法。（见sob的相关帖子）

作者
Author: CCCC 时间: 2018-3-7 15:32

sobereva 发表于 2018-3-5 15:43
不要用pseudo=cards
怎么使用混合基组和赝势看
详解Gaussian中混合基组、自定义基组和赝势基组的输入

SOb老师我昨天重新改了输入文件，今天收敛了，但是在L1110模块出问题了，没有遇到过，还请老师指教！以下是输出文件中的出错部分：
                           The (ff|ff) Classes
Not enough memory for any path:
MxKBra= 25 MxKKet= 25 Memory= 13326052.
After MaxDer, Avail=FFFFFFFFFFFFFFFFFFFFFFFFFFFFFTFFFFFFFFFFFFFFFFFFFF
After MemPth, Avail=FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF
After MxKxxx, Avail=FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF
After UseAlp/DoOS, Avail=FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF
No path selected in FillAv.
Error termination via Lnk1e in /export/app/g03/l1110.exe at Tue Mar  6 22:48:19 2018.

                           -------------------
                           The (ff|ff) Classes
Not enough memory for any path:
MxKBra= 25 MxKKet= 25 Memory= 13326052.
After MaxDer, Avail=FFFFFFFFFFFFFFFFFFFFFFFFFFFFFTFFFFFFFFFFFFFFFFFFFF
After MemPth, Avail=FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF
After MxKxxx, Avail=FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF
After UseAlp/DoOS, Avail=FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF
No path selected in FillAv.
Error termination via Lnk1e in /export/app/g03/l1110.exe at Tue Mar  6 22:48:20 2018.

作者
Author: sobereva 时间: 2018-3-7 17:40

CCCC 发表于 2018-3-7 15:32
SOb老师我昨天重新改了输入文件，今天收敛了，但是在L1110模块出问题了，没有遇到过，还请老师指教！以 ...

从这上看不出来。
尝试g09或g16，届时如果还有问题，上传压缩后的输入和输出文件

作者
Author: CCCC 时间: 2018-3-7 19:46

sobereva 发表于 2018-3-7 17:40
从这上看不出来。
尝试g09或g16，届时如果还有问题，上传压缩后的输入和输出文件

是这样的老师，我现在借的其他老师的服务器来算，上面有g03 g09 g16，我提交作业的时候输入的命令是 g09<*.com>*.out &,可是从输出文件来看用的是g03,这是什么原因？量化初学者，问的问题没什么深度，谢谢老师回复！

作者
Author: sobereva 时间: 2018-3-7 20:51

CCCC 发表于 2018-3-7 19:46
是这样的老师，我现在借的其他老师的服务器来算，上面有g03 g09 g16，我提交作业的时候输入的命令是 g09* ...

最好谁管服务器问谁。远程不好说清楚
自己排查的话，运行which g09，看g09命令是否确实是指向g09目录下的g09可执行文件

作者
Author: CCCC 时间: 2018-3-7 20:58

sobereva 发表于 2018-3-7 20:51
最好谁管服务器问谁。远程不好说清楚
自己排查的话，运行which g09，看g09命令是否确实是指向g09目录 ...

[*@g0112 ~/c72]$ which g09
/export/app/g09/g09
是g09……

作者
Author: CCCC 时间: 2018-3-8 09:14

sobereva 发表于 2018-3-7 20:51
最好谁管服务器问谁。远程不好说清楚
自己排查的话，运行which g09，看g09命令是否确实是指向g09目录 ...

sob老师我已经找管理员老师改好了谢谢老师

作者
Author: CCCC 时间: 2018-3-8 19:17

sobereva 发表于 2018-3-7 20:51
最好谁管服务器问谁。远程不好说清楚
自己排查的话，运行which g09，看g09命令是否确实是指向g09目录 ...

老师，我换成g09之后还是出现了L1110的错误，输入文件和输出文件在压缩文件中。

作者
Author: sobereva 时间: 2018-3-8 21:29

CCCC 发表于 2018-3-8 19:17
老师，我换成g09之后还是出现了L1110的错误，输入文件和输出文件在压缩文件中。

内存不够

要么减小基组到6-31G*，要么加大内存，要么尝试减小并行核数

作者
Author: CCCC 时间: 2018-3-10 11:36

sobereva 发表于 2018-3-8 21:29
内存不够

要么减小基组到6-31G*，要么加大内存，要么尝试减小并行核数

谢谢老师已经顺利解决！再请教一个小问题，对于基组6-311++G(3df,3pd) 和6-311++G**，
我刚刚在G09上计算单点能，用前者的时候一直提示我L302错误改成后者之后L302模块已经算过了，请问这是什么原因？

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