计算化学公社

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

[Gaussian/gview] 高斯中使用双杂DSD-PBEP86优化报错

[复制链接 Copy URL]

19

帖子

0

威望

581

eV
积分
600

Level 4 (黑子)

本帖最后由 Luckynan 于 2022-11-18 15:21 编辑

在使用高斯16优化结构时,使用双杂化泛函DSD-PBEP86计算报错。关键词#popt DSDPBEP86/6-31g(d) DENSITY,想请教一下,该如何修改?是什么原因?输入文件及输出文件如下。

202211181521203003..png (90.9 KB, 下载次数 Times of downloads: 16)

202211181521203003..png

A1DSD.gjf

4.73 KB, 下载次数 Times of downloads: 15

输入文件

A1DSD.log

103.53 KB, 下载次数 Times of downloads: 3

输出文件

1万

帖子

0

威望

9859

eV
积分
22096

Level 6 (一方通行)

2#
发表于 Post on 2022-11-18 15:54:23 | 只看该作者 Only view this author
先不论报错原因,这个计算即使正常结束结果也是不能用的,因为双杂化泛函至少要搭配3-zeta基组。双杂化泛函搭配2-zeta基组,结果还不如普通泛函搭配3-zeta基组靠谱。
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
本团队长期招收研究生,有意者可私信联系

538

帖子

1

威望

5763

eV
积分
6321

Level 6 (一方通行)

3#
发表于 Post on 2022-11-18 15:55:58 | 只看该作者 Only view this author
结构优化没有必要用双杂化泛函

211

帖子

0

威望

572

eV
积分
783

Level 4 (黑子)

4#
发表于 Post on 2022-11-18 16:26:04 | 只看该作者 Only view this author
看倒数第二行,硬盘使用超过了40GW(320GB),可能是超过了限制然后进程被杀掉了。这种体系用高斯算双杂化泛函的梯度有些勉为其难,硬要算可以尝试用ORCA开RI算。

6万

帖子

99

威望

6万

eV
积分
125140

管理员

公社社长

5#
发表于 Post on 2022-11-18 18:09:15 | 只看该作者 Only view this author
写DENSITY毫无意义

你这种体系用PBE0-D3(BJ)结合6-311G*优化得妥妥的,用双杂化结合6-31G*结果只会更烂,耗时还多两个数量级,而且rwf文件还会占极大的硬盘,临时文件目录写满了任务就会崩溃
北京科音自然科学研究中心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

19

帖子

0

威望

581

eV
积分
600

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2022-11-18 20:51:26 | 只看该作者 Only view this author
sobereva 发表于 2022-11-18 18:09
写DENSITY毫无意义

你这种体系用PBE0-D3(BJ)结合6-311G*优化得妥妥的,用双杂化结合6-31G*结果只会更烂 ...

好的,谢谢老师。还想问一下,我主要计算化合物的CT激发,使用cam-b3lyp和wb97xd均无法很好的符合实验值,PBE0-D3(BJ)的吸收值也会比实验值蓝移25nm左右,我有必要使用自定义泛函嘛?

1万

帖子

0

威望

9859

eV
积分
22096

Level 6 (一方通行)

7#
发表于 Post on 2022-11-18 21:22:46 | 只看该作者 Only view this author
Luckynan 发表于 2022-11-18 13:51
好的,谢谢老师。还想问一下,我主要计算化合物的CT激发,使用cam-b3lyp和wb97xd均无法很好的符合实验值 ...

25nm已经属于TDDFT误差范围以内了(当然,如果你算的是极紫外区的吸收,那还可以指望误差做得比25nm小),一般来说没必要也很难算得更准,即使更准也是碰巧。
一定要算得更准的话,第一必须做振动分辨光谱,第二必须显式考虑溶剂,第三必须做构象平均(包括溶剂分子的构象),第四激发能必须在DLPNO-STEOM-CCSD或更高级别下计算,TDDFT是不行的。
此外,如果你只研究吸收的话,结构优化和TDDFT的理论级别可以不统一,不要因为看到双杂化泛函算TDDFT准,就也用双杂化泛函做结构优化。
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
本团队长期招收研究生,有意者可私信联系

本版积分规则 Credits rule

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

GMT+8, 2026-2-20 19:14 , Processed in 0.257356 second(s), 28 queries , Gzip On.

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