计算化学公社

 找回密码 Forget password
 注册 Register
Views: 1074|回复 Reply: 16
打印 Print 上一主题 Last thread 下一主题 Next thread

[综合交流] sobEDA计算时静电力,交换力以及DFT关联能都为0是为什么?

[复制链接 Copy URL]

40

帖子

0

威望

279

eV
积分
319

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 Wenlin 于 2024-12-2 08:09 编辑

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


fragment.txt

22 Bytes, 下载次数 Times of downloads: 1

result.txt

2.3 KB, 下载次数 Times of downloads: 1

sobEDA.sh

14.56 KB, 下载次数 Times of downloads: 1

system.xyz

6.33 KB, 下载次数 Times of downloads: 2

template.gjf

249 Bytes, 下载次数 Times of downloads: 9

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

2#
发表于 Post on 2024-11-29 11:57:19 | 只看该作者 Only view this author
光凭这些没法判断,按照下文FAQ部分的红字判断
使用sobEDA和sobEDAw方法做非常准确、快速、方便、普适的能量分解分析
http://sobereva.com/685http://bbs.keinsci.com/thread-39446-1-1.html

若无绝对必要别带弥散,本来锕系体系就SCF难收敛,带了aug-雪上加霜。而如果是弱相互作用体系,也没理由B3LYP不带D3校正
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

40

帖子

0

威望

279

eV
积分
319

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2024-11-29 16:10:58 | 只看该作者 Only view this author
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

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

4#
发表于 Post on 2024-12-2 04:43:29 | 只看该作者 Only view this author
Wenlin 发表于 2024-11-29 16:10
老师,我是看了您上面说的那个帖子里常见的问题,应该是没有问题的,不加色散校正是因为我计算的体系里面 ...

没具体文件没法说
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

40

帖子

0

威望

279

eV
积分
319

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2024-12-2 08:10:55 | 只看该作者 Only view this author
sobereva 发表于 2024-12-2 04:43
没具体文件没法说

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

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

6#
发表于 Post on 2024-12-2 15:49:00 | 只看该作者 Only view this author
你的template.gjf末尾没写对。改完后跑成功了,相关文件: Am_coord.rar (131.26 KB, 下载次数 Times of downloads: 14)
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

1万

帖子

0

威望

9868

eV
积分
22108

Level 6 (一方通行)

7#
发表于 Post on 2024-12-2 15:53:40 | 只看该作者 Only view this author
Wenlin 发表于 2024-11-29 09:10
老师,我是看了您上面说的那个帖子里常见的问题,应该是没有问题的,不加色散校正是因为我计算的体系里面 ...

ORCA支持超重元素的DFT-D4校正,最重支持到Lr,可以用ORCA来算
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

40

帖子

0

威望

279

eV
积分
319

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2024-12-2 17:20:38 | 只看该作者 Only view this author
sobereva 发表于 2024-12-2 15:49
你的template.gjf末尾没写对。改完后跑成功了,相关文件:

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

40

帖子

0

威望

279

eV
积分
319

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2024-12-2 17:21:03 | 只看该作者 Only view this author
wzkchem5 发表于 2024-12-2 15:53
ORCA支持超重元素的DFT-D4校正,最重支持到Lr,可以用ORCA来算

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

40

帖子

0

威望

279

eV
积分
319

Level 3 能力者

10#
 楼主 Author| 发表于 Post on 2024-12-2 20:41:10 | 只看该作者 Only view this author
本帖最后由 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
*

1万

帖子

0

威望

9868

eV
积分
22108

Level 6 (一方通行)

11#
发表于 Post on 2024-12-2 22:15:48 | 只看该作者 Only view this author
Wenlin 发表于 2024-12-2 13:41
老师,我用orca5.0.3计算了dftd4但是报错了是因为版本不支持吗?
报错提示:WARNING: No D3 parameters  ...

从6.0.1(也可能是6.0.0)版开始才有
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

12#
发表于 Post on 2024-12-2 22:20:31 | 只看该作者 Only view this author
Wenlin 发表于 2024-12-2 20:41
老师,我用orca5.0.3计算了dftd4但是报错了是因为版本不支持吗?
报错提示:WARNING: No D3 parameters  ...

ORCA 6.0.0的D4已经支持到最高103号元素
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

105

帖子

0

威望

1822

eV
积分
1927

Level 5 (御坂)

13#
发表于 Post on 2024-12-2 23:13:28 | 只看该作者 Only view this author
d3bj和d4都有的情况下,现在推荐用哪一种更好?

40

帖子

0

威望

279

eV
积分
319

Level 3 能力者

14#
 楼主 Author| 发表于 Post on 2024-12-3 08:36:34 | 只看该作者 Only view this author
wzkchem5 发表于 2024-12-2 22:15
从6.0.1(也可能是6.0.0)版开始才有

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

40

帖子

0

威望

279

eV
积分
319

Level 3 能力者

15#
 楼主 Author| 发表于 Post on 2024-12-3 08:37:07 | 只看该作者 Only view this author
sobereva 发表于 2024-12-2 22:20
ORCA 6.0.0的D4已经支持到最高103号元素

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

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-21 23:48 , Processed in 0.220697 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list