请教关于Linux系统下使用optDFTw进行w值的计算过程

yezhonghua · 发表于 Post on 2017-1-4 17:49:23

本帖最后由 yezhonghua 于 2017-1-4 18:03 编辑

请教一下linux下算过RS w优化的大神一下：
sober老师在http://sobereva.com/346提供了一个optDFTw linux程序；由于接触linux的比较浅，所以不太明白如何进行计算。请教一下做过的大神计算步骤：
templet.gjf是单点能计算文件（先用b3lyp/6-31g*进行opt freq得到比较好的构型,然后用LC-wPBE/6-31g*进行单点能）。
按照老师提供的做法我将gjf文件命名为templet.gjf 将optDFTw权限改为775，房子同一个目录下，执行./optDFTw。然后出现图1和2所示情况，出现报错。
请教一下是否是我的g09的问题或者是我的计算过程有问题，谢谢！
附件为templateg.gjf。
谢谢老师

sobereva · 发表于 Post on 2017-1-4 18:10:48

你先搞个最简单的水体系，如果能成功，说明9成是当前体系SCF不收敛所致，观看高斯输出文件以确认。
解决方法：
解决SCF不收敛问题的方法
http://sobereva.com/61

yezhonghua · 发表于 Post on 2017-1-4 18:14:48

sobereva 发表于 2017-1-4 18:10
你先搞个最简单的水体系，如果能成功，说明9成是当前体系SCF不收敛所致，观看高斯输出文件以确认。
解决方 ...

谢谢sober老师，我的计算过程应该没有问题吧?

sobereva · 发表于 Post on 2017-1-4 18:15:54

没问题

yezhonghua · 发表于 Post on 2017-1-4 19:27:36

本帖最后由 yezhonghua 于 2017-1-4 19:46 编辑

sobereva 发表于 2017-1-4 18:10
你先搞个最简单的水体系，如果能成功，说明9成是当前体系SCF不收敛所致，观看高斯输出文件以确认。
解决方 ...

谢谢sober老师。
我采用高斯提供的联苯分子进行计算也出现上述问题，提示找不到阿尔法轨道本征值。
我用的是09.A.02版本，我按照附件的提示网上说可能是我的高斯版本不支持wPBE的IOP？

yezhonghua · 发表于 Post on 2017-1-4 19:39:58

sobereva 发表于 2017-1-4 18:15
没问题

谢谢sober老师，我又用联苯分子试了一下，也是一样的错误，是否可能是我的g9A版本不支持wPBE 的IOP的原因吧？

yezhonghua · 发表于 Post on 2017-1-4 19:58:15

参照高斯手册应该是支持LC-wPBE的
输入文件如下：
%chk=template.chk
%mem=32GB
%nprocshared=24
#p LC-wPBE/6-31g(d)

但是计算输出却出现了报错（如附件），网上查了一下所A版本不支持LC-wPBE的IOP，需要升级到D.01版本。
请教一下大家是否遇到过这个问题，谢谢！

ggdh · 发表于 Post on 2017-1-4 20:23:04

yezhonghua 发表于 2017-1-4 19:39
谢谢sober老师，我又用联苯分子试了一下，也是一样的错误，是否可能是我的g9A版本不支持wPBE 的IOP的原因 ...

要做调控起码要B版本的gaussian

sobereva · 发表于 Post on 2017-1-4 20:48:47

A.02假若不支持LC-wPBE，一开始就报错了。你用E或者D.01跑一次，如果正常，说明是A.02的bug

yezhonghua · 发表于 Post on 2017-1-5 08:31:00

sobereva 发表于 2017-1-4 20:48
A.02假若不支持LC-wPBE，一开始就报错了。你用E或者D.01跑一次，如果正常，说明是A.02的bug

多谢sober老师

yezhonghua · 发表于 Post on 2017-1-5 11:10:42

本帖最后由 yezhonghua 于 2017-1-5 11:22 编辑

sobereva 发表于 2017-1-4 20:48
A.02假若不支持LC-wPBE，一开始就报错了。你用E或者D.01跑一次，如果正常，说明是A.02的bug

sober老师，我换用09D版本的就没有问题了，可能还是版本的问题。
另外还要麻烦sober老师一下：
1、对于这个优化w值，我看您是直接计算单点能的。是否应该先用B3LYP/6-31g*进行基态(N.gjf)的opt+frq，得到一个较好的构型，然后用您提供的程序进行（LC-wPBE泛函下）w优化单点能计算。
2、LC-wPBE的w优化直接采用单点能优化是不是很妥，对于80个原子体系的有机分子，是否在template.gjf文件中加入 opt freq 进行构型优化（至少应该加入opt吧），但是这样的运算10几圈的任务计算量肯定很大。所以，想问一下老师这个直接通过单点能计算优化得到的w值有没有必要进行opt freq计算？
如果有冒犯之处，还请老师见谅，主要是我觉得初始构型可能对于w的值存在一定的影响。
谢谢老师!

冰释之川 · 发表于 Post on 2017-1-5 17:02:50

本帖最后由冰释之川于 2017-1-5 17:10 编辑

yezhonghua 发表于 2017-1-5 11:10
sober老师，我换用09D版本的就没有问题了，可能还是版本的问题。
另外还要麻烦sober老师一下：
1、对于 ...

1. w值的优化肯定是基于优化好的基态构型坐标来找的。换句话说，你要在该结构下做一系列不同w值的单点能。2. 不用在template.gjf模版里改成opt freq，这么做不符合IP-tuning的原理！

总之，泛函调控都是基于优化好的分子构型下来实施的！，后续寻找最佳w值则不需要画蛇添足地再优化什么结构了！

yezhonghua · 发表于 Post on 2017-1-5 20:29:42

冰释之川发表于 2017-1-5 17:02
1. w值的优化肯定是基于优化好的基态构型坐标来找的。换句话说，你要在该结构下做一系列不同w值的单点能 ...

多谢冰释大神

yezhonghua · 发表于 Post on 2017-1-6 21:31:19

sobereva 发表于 2017-1-4 20:48
A.02假若不支持LC-wPBE，一开始就报错了。你用E或者D.01跑一次，如果正常，说明是A.02的bug

sober 老师，我优化文献上一个结构，文献值w=0.183；但是我优化了14圈了，J=0.002，依然没有收敛；但是w出现下降又上升又下降的趋势（如附件1所示），是不是可以确定w在0.1828-0.1830之间。
谢谢！

冰释之川 · 发表于 Post on 2017-1-6 22:03:23

yezhonghua 发表于 2017-1-6 21:31
sober 老师，我优化文献上一个结构，文献值w=0.183；但是我优化了14圈了，J=0.002，依然没有收敛；但是w ...

由于Sob老师的程序采用的是基于Brent算法来最小化J^2的，一次次迭代有高低起伏是正常的，你等结果跑完即可。默认的收敛限下，最多跑个十七、八圈w值就收敛了

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

[Gaussian/gview] 请教关于Linux系统下使用optDFTw进行w值的计算过程

g09A .02 是否支持LC-wPBE的IOP

评分 Rate

浏览过的版块