计算化学公社

标题: 用B3LYP/6-31G*算HOMO LUMO轨道出现错误，求助大神 [打印本页]

作者
Author: 18522780016 时间: 2018-8-9 17:28
标题: 用B3LYP/6-31G*算HOMO LUMO轨道出现错误，求助大神
大家好，sob老师好，小弟是计算新手，最近刚学会了用高斯计算分子的HOMO LUMO轨道，但是最近用B3LYP/6-31G*计算的时候不能正常终止，而利用B3LYP/6-31G时就能正常终止，但是听说6-31G的精度太差，想问问大家有什么解决方法吗？输入文件以及输出文件如下
之前在别的帖子中看到，如果出现像输出文件中“warning！！”那样的错误（图中红框标出部分）如果能正常终止就没关系，如果不能正常终止就要用不含有弥散的基组，想问一下大家有哪些是不含有弥散的基组，可能问题比较低级，希望大家不要嫌弃，谢谢了
(, 下载次数 Times of downloads: 59) (, 下载次数 Times of downloads: 51)

作者
Author: wrtgcn 时间: 2018-8-9 18:18
直接看error那行，然后在http://bbs.keinsci.com/thread-4829-1-1.html找到对应的问题和解决方法。比如你这个错误对应：
“成因：已知、且已修复的Gaussian程序Bug，自G09 C.01版本之后已修复
（G09 C.01 Release Note：“A memory allocation bug for very large systems, which could cause a failure with the message “NIJ > Max2 in MMCore,” was fixed”）
解决：更换G09 C.01或之后版本。”

作者
Author: 18522780016 时间: 2018-8-9 18:25

wrtgcn 发表于 2018-8-9 18:18
直接看error那行，然后在http://bbs.keinsci.com/thread-4829-1-1.html找到对应的问题和解决方法。比如你这 ...

非常感谢，我是利用学校的超算中心计算的，想再问您一下怎么能用G09 C.01或者以后的版本计算？是不是将下列代码中的g09替换成G09 C.01就可以啊？
#!/bin/bash
#SBATCH --nodes=1
#SBATCH --tasks-per-node=1
#SBATCH --cpus-per-task=8
#SBATCH --mail-type=END
#SBATCH --mail-user=boxu.feng@iapp.de
#SBATCH --time=2400:00:00 # walltime
#SBATCH --partition=haswell
#SBATCH --mem-per-cpu=2583M # memory per CPU core
#SBATCH -A solar

module load python/2.7.5-scipy
module load gaussian/g09
module load turbomole/6.5

export OMP_NUM_THREADS=4
OUTFILE="thbf-1.out"

g09 < thbf-1.inp > "$OUTFILE"

exit 0

作者
Author: 18522780016 时间: 2018-8-9 18:27

18522780016 发表于 2018-8-9 18:25
非常感谢，我是利用学校的超算中心计算的，想再问您一下怎么能用G09 C.01或者以后的版本计算？是不是将下 ...

这段代码是在一个文档中的，和.inp文件一起上传到服务器上的

作者
Author: wrtgcn 时间: 2018-8-9 18:54

18522780016 发表于 2018-8-9 18:25
非常感谢，我是利用学校的超算中心计算的，想再问您一下怎么能用G09 C.01或者以后的版本计算？是不是将下 ...

我们组的服务器只有一个高斯版本，一般就是直接用了。
如果你们同时有多个高斯版本的话，不太清楚是否可以用改参数的方式选用其他版本的高斯（我个人感觉应该不会有了新版还默认用旧版的吧）。这个可能你得联系相关的管理者问一下。

作者
Author: 18522780016 时间: 2018-8-9 19:28

wrtgcn 发表于 2018-8-9 18:54
我们组的服务器只有一个高斯版本，一般就是直接用了。
如果你们同时有多个高斯版本的话，不太清楚是否可 ...

好的，非常感谢

作者
Author: sobereva 时间: 2018-8-10 04:56
认真阅读
谈谈量子化学中基组的选择
http://sobereva.com/336
谈谈弥散函数和“月份”基组
http://sobereva.com/119

当前问题跟基组是否带弥散根本毫无联系，别把极化和弥散搞混，不带极化的6-31G绝对不能用
那两行warning是吓唬人的，直接无视
常见问题先在论坛首页google框里搜贴，常见问题>99%论坛里都讨论过

作者
Author: alonewolfyang 时间: 2018-8-10 09:24
如题，sob老师的这两个博文已经很好的回答了你的问题，但是个人感觉，计算前线轨道能量的话，想要计算准确，目前比较好的方法是使用调节w参数泛函的方法，Bredas组目前都用的是这个方法

作者
Author: 18522780016 时间: 2018-8-10 15:57

sobereva 发表于 2018-8-10 04:56
认真阅读
谈谈量子化学中基组的选择
http://sobereva.com/336

感谢老师

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