计算化学公社

标题: sobEDA计算时静电力，交换力以及DFT关联能都为0是为什么？ [打印本页]

作者
Author: Wenlin 时间: 2024-11-29 11:53
标题: sobEDA计算时静电力，交换力以及DFT关联能都为0是为什么？
本帖最后由 Wenlin 于 2024-12-2 08:09 编辑

使用高斯进行sobEDA计算后所得片段的部分能量不显示，且最终的静电力交换力以及DFT关联能都为0，求助各位这是什么原因造成的啊？以下分别是我的输入文件还有result文本文件
(, 下载次数 Times of downloads: 3) (, 下载次数 Times of downloads: 3)

作者
Author: sobereva 时间: 2024-11-29 11:57
光凭这些没法判断，按照下文FAQ部分的红字判断
使用sobEDA和sobEDAw方法做非常准确、快速、方便、普适的能量分解分析
http://sobereva.com/685（http://bbs.keinsci.com/thread-39446-1-1.html）

若无绝对必要别带弥散，本来锕系体系就SCF难收敛，带了aug-雪上加霜。而如果是弱相互作用体系，也没理由B3LYP不带D3校正

作者
Author: Wenlin 时间: 2024-11-29 16:10

sobereva 发表于 2024-11-29 11:57
光凭这些没法判断，按照下文FAQ部分的红字判断
使用sobEDA和sobEDAw方法做非常准确、快速、方便、普适的能 ...

老师，我是看了您上面说的那个帖子里常见的问题，应该是没有问题的，不加色散校正是因为我计算的体系里面有Am元素为94以后的元素  准备不加色散校正后面再手动计算DFT+D4. 然后我把输入文件中的基组改为了6-311G**,不过计算得到的结果还是一样，静电力等能量为0 而且每个片段部分的静电力等都是不显示数值的  连“0”都没有  就像这样

Energy components of fragment 2:
E_tot = -745.156195497 Hartree
E_T = 350.693200 Hartree
E_els =  Hartree
E_x =  Hartree
E_c =  Hartree
E_disp = 0 Hartree

作者
Author: sobereva 时间: 2024-12-2 04:43

Wenlin 发表于 2024-11-29 16:10
老师，我是看了您上面说的那个帖子里常见的问题，应该是没有问题的，不加色散校正是因为我计算的体系里面 ...

没具体文件没法说

作者
Author: Wenlin 时间: 2024-12-2 08:10

sobereva 发表于 2024-12-2 04:43
没具体文件没法说

老师，我把相关的文件都添加在帖子里了，麻烦您再看一下

作者
Author: sobereva 时间: 2024-12-2 15:49
你的template.gjf末尾没写对。改完后跑成功了，相关文件： (, 下载次数 Times of downloads: 14)

作者
Author: wzkchem5 时间: 2024-12-2 15:53

Wenlin 发表于 2024-11-29 09:10
老师，我是看了您上面说的那个帖子里常见的问题，应该是没有问题的，不加色散校正是因为我计算的体系里面 ...

ORCA支持超重元素的DFT-D4校正，最重支持到Lr，可以用ORCA来算

作者
Author: Wenlin 时间: 2024-12-2 17:20

sobereva 发表于 2024-12-2 15:49
你的template.gjf末尾没写对。改完后跑成功了，相关文件：

好的好的，后面好好注意格式，感谢老师

作者
Author: Wenlin 时间: 2024-12-2 17:21

wzkchem5 发表于 2024-12-2 15:53
ORCA支持超重元素的DFT-D4校正，最重支持到Lr，可以用ORCA来算

好的老师，我再用orca计算一下

作者
Author: Wenlin 时间: 2024-12-2 20:41
本帖最后由 Wenlin 于 2024-12-2 20:58 编辑

wzkchem5 发表于 2024-12-2 15:53
ORCA支持超重元素的DFT-D4校正，最重支持到Lr，可以用ORCA来算

老师，我用orca5.0.3计算了dftd4但是报错了是因为版本不支持吗？
报错提示：WARNING: No D3 parameters available for elements with Z > 94!
  ===> : Skipping actual calculation
[file orca_main/mainchk.cpp, line 4930]: Error (ORCA_MAIN): ... aborting the run

[file orca_main/mainchk.cpp, line 4930]: Error (ORCA_MAIN): ... aborting the run

输入文件：
%maxcore  5000
%pal nprocs  16 end
!B3LYP D4
! def2-TZVP def2/J SARC/J  RIJCOSX tightSCF noautostart miniprint nopop UKS
! KDIIS
!NOTRAH
%basis
Newgto AM
S 7
1    1000.1312000          -0.0001300
2       62.3077550          -0.0067790
省略
end
* xyz 2 4
C    15.59040197 10.99742361 3.02638091
省略
Am 14.32106037 20.78602400 -0.29022423
O    13.10715664 20.66283236 -1.48807777
O    15.55559399 20.88893348 0.89374186
*

作者
Author: wzkchem5 时间: 2024-12-2 22:15

Wenlin 发表于 2024-12-2 13:41
老师，我用orca5.0.3计算了dftd4但是报错了是因为版本不支持吗？
报错提示：WARNING: No D3 parameters ...

从6.0.1（也可能是6.0.0）版开始才有

作者
Author: sobereva 时间: 2024-12-2 22:20

Wenlin 发表于 2024-12-2 20:41
老师，我用orca5.0.3计算了dftd4但是报错了是因为版本不支持吗？
报错提示：WARNING: No D3 parameters ...

ORCA 6.0.0的D4已经支持到最高103号元素

作者
Author: michaelm 时间: 2024-12-2 23:13
d3bj和d4都有的情况下，现在推荐用哪一种更好？

作者
Author: Wenlin 时间: 2024-12-3 08:36

wzkchem5 发表于 2024-12-2 22:15
从6.0.1（也可能是6.0.0）版开始才有

好的感谢老师，更新到601版本后已经正常计算了

作者
Author: Wenlin 时间: 2024-12-3 08:37

sobereva 发表于 2024-12-2 22:20
ORCA 6.0.0的D4已经支持到最高103号元素

好的感谢老师，已经用新版本正常计算了

作者
Author: wzkchem5 时间: 2024-12-3 16:38

michaelm 发表于 2024-12-2 16:13
d3bj和d4都有的情况下，现在推荐用哪一种更好？

对于原子带电不多的体系，D4和D3BJ差不多好；对于原子带电数多的体系，D4比D3BJ好

作者
Author: sobereva 时间: 2024-12-4 16:35

michaelm 发表于 2024-12-2 23:13
d3bj和d4都有的情况下，现在推荐用哪一种更好？

可以默认用D4
D4比D3主要改进在过渡金属配合物、存在大小较大原子电荷（如离子晶体）等情况。对有机体系没什么区别

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