请教关于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。
谢谢老师

510522199508287 · 发表于 Post on 2019-3-1 14:32:03

sobereva 发表于 2019-3-1 00:21
是

解决SCF不收敛问题的方法

谢谢老师，在找更高版本的gaussian

sobereva · 发表于 Post on 2019-3-1 00:21:14

510522199508287 发表于 2019-2-28 19:19
你好，我的跟你一样L502报错，怎么看版本呢， Gaussian 09: EM64L-G09RevA.01 这个是A版吗？

是

解决SCF不收敛问题的方法
http://sobereva.com/61

510522199508287 · 发表于 Post on 2019-2-28 19:19:07

yezhonghua 发表于 2017-1-4 19:27
谢谢sober老师。
我采用高斯提供的联苯分子进行计算也出现上述问题，提示找不到阿尔法轨道本征值。
...

你好，我的跟你一样L502报错，怎么看版本呢， Gaussian 09: EM64L-G09RevA.01 这个是A版吗？

阳光总在风雨后 · 发表于 Post on 2017-11-13 08:23:40

sobereva 发表于 2017-11-9 15:43
PBS的事我不清楚，如果最终结果都一样，那计算的估计是同一个分子。如果最终结构都和初猜结果一样，那根 ...

恩恩谢谢老师！

sobereva · 发表于 Post on 2017-11-9 15:43:51

阳光总在风雨后发表于 2017-11-4 18:14
sob老师您好！我在用optDFTw进行w值的计算时（在linux系统下），计算的不同的分子的w的值都是一样的，不 ...

PBS的事我不清楚，如果最终结果都一样，那计算的估计是同一个分子。如果最终结构都和初猜结果一样，那根本没计算。

阳光总在风雨后 · 发表于 Post on 2017-11-4 18:14:04

sobereva 发表于 2017-1-6 22:44
可以确定。如果反复不收敛，可以在这个小区间里再用scan方式确定

sob老师您好！我在用optDFTw进行w值的计算时（在linux系统下），计算的不同的分子的w的值都是一样的，不知道问题怎么改。您看一下我的脚本有什么问题：
#! /bin/sh
#PBS -S /bin/bash
#PBS -N BPy-p3C-wB97X-s0-1.gjf
#PBS -j oe

# Envrionment setting for Gaussian09
#
#PBS -o ${PBS_JOBNAME}.err
#PBS -l nodes=node20:ppn=20
#PBS -V
JOBNAME=${PBS_JOBNAME}
### FOR Gauss09 #####
export g09root=/public/software/Gaussian09
DATE=`date +"%R-%F-%A"`
WHOAMI=`whoami`_${PBS_JOBID}_${DATE}
mkdir /tmp/${WHOAMI}
export GAUSS_SCRDIR=/tmp/$WHOAMI
. $g09root/g09/bsd/g09.profile
#
# End of envrionment setting for Gaussian09

cp -r $PBS_O_WORKDIR/* /tmp/${WHOAMI}
cd /tmp/${WHOAMI}

./optDFTw

cp -r /tmp/${WHOAMI}/*  $PBS_O_WORKDIR
rm -rf /tmp/${WHOAMI}

石头 · 发表于 Post on 2017-5-24 17:24:09

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

您好，请问您换到D版本都可以正常跑完了吗？我在调控W时也是出现和您类似的错误，我用的B版本。

yezhonghua · 发表于 Post on 2017-5-23 21:15:14

sobereva 发表于 2017-5-23 13:42
博文里明确说了，这三个通常只对因为弥散函数导致的不收敛管用
绝对不要以为加了这三个就万事大吉了

谢谢sober老师，事实上，这几个加上去明显提升SCF收敛效果。

yezhonghua · 发表于 Post on 2017-5-23 21:14:38

冰释之川发表于 2017-5-23 13:46
你w优化本身就是基于优化好的几何构型来执行的，什么叫“改变一下优化后的结构”，几何结构不能再动来动 ...

谢谢大神，我按照sober老师的帖子再试试

冰释之川 · 发表于 Post on 2017-5-23 13:46:42

本帖最后由冰释之川于 2017-5-23 13:48 编辑

yezhonghua 发表于 2017-5-23 13:36
谢谢！我用的表头是#p LC-wPBE/6-31g(d) SCF(novaracc,noincfock) int=acc2e=12
基本上算是考虑到SCF不 ...

你w优化本身就是基于优化好的几何构型来执行的，什么叫“改变一下优化后的结构”，几何结构不能再动来动去了？如果你使用了某种方法SCF不收敛，那么换其他方法试试。

sobereva · 发表于 Post on 2017-5-23 13:42:25

yezhonghua 发表于 2017-5-23 13:36
谢谢！我用的表头是#p LC-wPBE/6-31g(d) SCF(novaracc,noincfock) int=acc2e=12
基本上算是考虑到SCF不 ...

博文里明确说了，这三个通常只对因为弥散函数导致的不收敛管用
绝对不要以为加了这三个就万事大吉了

yezhonghua · 发表于 Post on 2017-5-23 13:36:23

冰释之川发表于 2017-5-23 11:39
结构不变，SCF不收敛参考《解决SCF不收敛问题的方法》http://sobereva.com/61

谢谢！我用的表头是#p LC-wPBE/6-31g(d) SCF(novaracc,noincfock) int=acc2e=12
基本上算是考虑到SCF不收敛的情况了。
我现在改成#p LC-wPBE/6-31g(d) nosymm SCF(novaracc,noincfock) int=acc2e=12
如果不行的话，我可否成尝试改变一下优化后的结构？

冰释之川 · 发表于 Post on 2017-5-23 11:39:45

yezhonghua 发表于 2017-5-23 09:03
sober老师，请教一下我在基态优化好的基础上进行w值搜索，开始的N.gjf和N-1.gjf SCF都能收敛，但是N+1.gj ...

结构不变，SCF不收敛参考《解决SCF不收敛问题的方法》http://sobereva.com/61

yezhonghua · 发表于 Post on 2017-5-23 09:03:21

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

sober老师，请教一下我在基态优化好的基础上进行w值搜索，开始的N.gjf和N-1.gjf SCF都能收敛，但是N+1.gjf无法收敛，导致计算中断。说明我的基态优化结果应该是没有问题的。
因为我是在优化好基态的基础上进行w优化的，所以我需要重新优化基态结构么？还是对已经优化好的基态进行稍微的结构调整？谢谢sober老师

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

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