计算化学公社

标题: 计算能量分解出错，求指导 [打印本页]

作者
Author: F凌兰 时间: 2016-12-26 16:56
标题: 计算能量分解出错，求指导

请问：我在计算IF-HF能量分解时输入文件如附件所示，输出结果出错，有点不太懂怎么处理，还望指教！
我的体系所用方法是wB97XD,F\H\F原子用6-311++G**基组，I原子用赝势SDD基组。
(, 下载次数 Times of downloads: 15)

作者
Author: sobereva 时间: 2016-12-26 17:47
$ECP里面按照F,H,F,I的顺序来定义。调过后还不行的话把输出文件上传

作者
Author: mengzi 时间: 2016-12-26 19:18
DFTTYP=wB97XD这个是不是应该是DFTTYP=wb97x-D ?

作者
Author: F凌兰 时间: 2016-12-27 09:49

mengzi 发表于 2016-12-26 19:18
DFTTYP=wB97XD这个是不是应该是DFTTYP=wb97x-D ?

好像对这个没什么影响，第一次还显示加上不行，后来也没提示，我是按不加 “-” 处理的

作者
Author: F凌兰 时间: 2016-12-27 09:59

sobereva 发表于 2016-12-26 17:47
$ECP里面按照F,H,F,I的顺序来定义。调过后还不行的话把输出文件上传

谢谢老师的指导，改了顺序还是出错
这个好像问题出在需要设置DFT网格，一直有这一项，如下所述：

DFT GRIDS ARE AVAILABLE ONLY FOR INTEGER ATOMIC
NUMBERS, SO -SPARKLES- ARE NOT ALLOWED.

但是我不太懂这个该如何设置，还望多多指教，谢谢！一并上传输入及输出文件
(, 下载次数 Times of downloads: 4) (, 下载次数 Times of downloads: 3)

作者
Author: sobereva 时间: 2016-12-27 11:10

F凌兰发表于 2016-12-27 09:59
谢谢老师的指导，改了顺序还是出错
这个好像问题出在需要设置DFT网格，一直有这一项，如下所述：

原子名后头写上原子序号，例如
$DATA
Molecule specification
C1
C 6.0             0.01646410 -0.58657651 0.00000000
H 1.0             -0.89975307 -1.13924829 0.00000000
O 8.0             -0.00849072 0.71448419 0.00000000
H 1.0             0.86889421 -1.05716614 0.00000000
$END

作者
Author: F凌兰 时间: 2016-12-27 15:15

sobereva 发表于 2016-12-27 11:10
原子名后头写上原子序号，例如
$DATA
Molecule specification

非常感谢

作者
Author: F凌兰 时间: 2016-12-27 15:25

mengzi 发表于 2016-12-26 19:18
DFTTYP=wB97XD这个是不是应该是DFTTYP=wb97x-D ?

不好意思，我又试了一下，这个确实要改的，应该加 “-”，谢谢。

作者
Author: F凌兰 时间: 2016-12-27 16:56

sobereva 发表于 2016-12-27 11:10
原子名后头写上原子序号，例如
$DATA
Molecule specification

请问老师：我的体系是做EDA，方法采用wB97X-D,基组采用赝势SDD//6-311++G**，输出结果总是报错
因为wB97X-D方法自带DFT-D2矫正，这也是我希望用的，您看是不是我的$DFT 行设置不对？
见附件输入输出文件，谢谢！
(, 下载次数 Times of downloads: 4) (, 下载次数 Times of downloads: 5)

作者
Author: sobereva 时间: 2016-12-27 18:28

F凌兰发表于 2016-12-27 16:56
请问老师：我的体系是做EDA，方法采用wB97X-D,基组采用赝势SDD//6-311++G**，输出结果总是报错
因为wB97 ...

已经提示了 DFT-D METHODS ARE NOT SUPPORTED

作者
Author: F凌兰 时间: 2016-12-27 20:49

sobereva 发表于 2016-12-27 18:28
已经提示了 DFT-D METHODS ARE NOT SUPPORTED

嗯嗯，谢谢老师。多有打扰，抱歉。
提示说不支持DFT-D，那我还能用wb97xd这个方法吗？这个方法是我要用的方法
但GAMESS说明中提示可以算色散，我按照说明中的设置尝试了一下，未提到的按照默认进行，
这样提示错误，还有其他办法可以保留wb97xd这个方法吗？
谢谢

作者
Author: sobereva 时间: 2016-12-28 08:29

F凌兰发表于 2016-12-27 20:49
嗯嗯，谢谢老师。多有打扰，抱歉。
提示说不支持DFT-D，那我还能用wb97xd这个方法吗？这个方法是我要用 ...

不能。本身LMOEDA也没考虑DFT-D校正的情况。
你用个诸如B3LYP这样的做能量分解，然后色散的贡献用grimme的DFT-D3程序算出来。

作者
Author: F凌兰 时间: 2016-12-28 09:13

sobereva 发表于 2016-12-28 08:29
不能。本身LMOEDA也没考虑DFT-D校正的情况。
你用个诸如B3LYP这样的做能量分解，然后色散的贡献用grim ...

哦哦，好的，非常感谢老师。

作者
Author: F凌兰 时间: 2016-12-28 20:53

sobereva 发表于 2016-12-28 08:29
不能。本身LMOEDA也没考虑DFT-D校正的情况。
你用个诸如B3LYP这样的做能量分解，然后色散的贡献用grim ...

老师，您好！
我想请教一下：
(, 下载次数 Times of downloads: 10) (, 下载次数 Times of downloads: 2)
见附件我的输入、输出文件（正常结束），其中输出文件中 TOTAL INTERACTION ENERGY HF OR DFT E=-0.009561 -6.00
其中，ENERGY HF OR DFT 是代表一样的？还是我的$CONTRL行有错误？

使用SCFTYP=RHF 和 DFTTYP 命令有冲突吗？SCFTYP总是要加的，还是分情况？
若用了SCFTYP就需要$SCF？还是收敛命令总是要加的？

谢谢！

作者
Author: sobereva 时间: 2016-12-28 21:22

F凌兰发表于 2016-12-28 20:53
老师，您好！
我想请教一下：

你做HF那里就是HF的相互作用能，你做DFT那里就是DFT的相互作用能
SCFTYP设哪种设定的是SCF类型而已，相当于高斯里面设R、U、RO

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