计算化学公社

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

[Gaussian/gview] 磁盘空间不足报错

[复制链接 Copy URL]

27

帖子

0

威望

355

eV
积分
382

Level 3 能力者

我用B3LYP结构优化以及频率优化后想做高精度的能量计算,但是出现了磁盘空间不足报错的问题。大概是两百多个原子CHONS元素。服务器系统存储有30多TB,磁盘空间肯定是足够的。g09这个临时文件夹tmp已经清空了,还是出现相同的报错,请问老师有什么办法可以解决这个问题吗?

202110242310544199..png (7.06 KB, 下载次数 Times of downloads: 61)

202110242310544199..png

202110242310145192..png (58.9 KB, 下载次数 Times of downloads: 63)

202110242310145192..png

1187

帖子

5

威望

2876

eV
积分
4163

Level 6 (一方通行)

2#
发表于 Post on 2021-10-24 23:41:07 | 只看该作者 Only view this author
本帖最后由 snljty 于 2021-10-25 09:03 编辑

肯定算不动,不用想了,100个都算不动。换ORCA。

227

帖子

0

威望

2800

eV
积分
3027

Level 5 (御坂)

3#
发表于 Post on 2021-10-25 08:33:54 | 只看该作者 Only view this author
用orca的ri-做双杂化,你的硬盘应该够了,内存不够。放弃吧

1245

帖子

1

威望

3860

eV
积分
5125

Level 6 (一方通行)

4#
发表于 Post on 2021-10-25 09:13:29 | 只看该作者 Only view this author
甭拿高斯算这玩意了
你得确定tmp文件夹下单个文件大小的上限是多少,这一算就会产生几TB一个的文件,你那文件系统不一定支持这么大的单个文件

194

帖子

0

威望

3694

eV
积分
3888

Level 5 (御坂)

5#
发表于 Post on 2021-10-25 09:36:10 | 只看该作者 Only view this author
mp2不用density fitting/ri局限性太大了

334

帖子

1

威望

6153

eV
积分
6507

Level 6 (一方通行)

6#
发表于 Post on 2021-10-25 10:02:14 | 只看该作者 Only view this author
以前在G16上试过 B2PLYPD3/def2tzvp 有机分子 100个左右的原子 2000个基函数左右,需要大概8TB的硬盘空间。(计算给了48核心180GB内存,单点能算了24小时)因为内存不够基本只达到了最低限(2OVN), 硬盘就会占的很多。高斯官方给的是MP2 Direct算法下硬盘空间跟基函数的平方成正比。那么我觉得你的体系30TB不够也属于正常现象。

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

7#
发表于 Post on 2021-10-25 11:28:09 | 只看该作者 Only view this author
Gaussian做MP2/双杂化结合def2-TZVPP这种档次的基组计算,50原子以上就根本甭考虑。再多就用ORCA
北京科音自然科学研究中心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

27

帖子

0

威望

355

eV
积分
382

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2021-10-26 20:40:11 | 只看该作者 Only view this author
那请问我这个三百个原子的分子如果想算溶质的自由能要采用什么方法呢?

202110262039273292..png (74.25 KB, 下载次数 Times of downloads: 45)

202110262039273292..png

27

帖子

0

威望

355

eV
积分
382

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2021-10-26 20:40:53 | 只看该作者 Only view this author
只能用第一种无法高精度计算那个吗

410

帖子

0

威望

2519

eV
积分
2929

Level 5 (御坂)

10#
发表于 Post on 2021-10-26 20:52:24 | 只看该作者 Only view this author
sobereva 发表于 2021-10-25 11:28
Gaussian做MP2/双杂化结合def2-TZVPP这种档次的基组计算,50原子以上就根本甭考虑。再多就用ORCA

请问一下高斯和ORCA算出来的能量差别大吗?

4289

帖子

4

威望

9541

eV
积分
13910

Level 6 (一方通行)

MOKIT开发者

11#
发表于 Post on 2021-10-26 21:23:13 | 只看该作者 Only view this author
本帖最后由 zjxitcc 于 2021-10-26 21:25 编辑
mfdsrax2 发表于 2021-10-26 20:52
请问一下高斯和ORCA算出来的能量差别大吗?

取决于你的关键词如何写;以及你想问绝对能量还是相对能量。

绝对能量:在合理关键词的情况下,例如气相下的HF, MP2, CCSD等可以严格几乎一样。各种泛函那影响因素就多了,积分格点、泛函定义的细节、SCF收敛限、积分精度、隐式溶剂模型的差异等等都有影响。所以这种问题意义不大。

相对能量(例如反应能垒):在合理关键词的情况下,所有软件算出来都趋势一致。
自动做多参考态计算的程序MOKIT

1万

帖子

0

威望

9859

eV
积分
22096

Level 6 (一方通行)

12#
发表于 Post on 2021-10-26 21:38:31 | 只看该作者 Only view this author
mfdsrax2 发表于 2021-10-26 13:52
请问一下高斯和ORCA算出来的能量差别大吗?

补充一下,HF, MP2, CCSD的差别做到10^-6 Hartree以内问题不大,但是从orca 5开始会默认开RIJCOSX,RIJCOSX会加速计算,但是和高斯的差别会扩大到10^-4 Hartree量级。
DFT的差别一般只能保证做到10^-4 Hartree量级,如果是不用RI、RIJCOSX的计算,差别可以缩小到10^-5 Hartree左右。
以上是假设用户已经知道了诸如B3LYP用VWN3和VWN5的区别、6-31G*用球谐基函数还是笛卡尔基函数的区别等等,把这些因素都控制成一样的情况下,才有以上的结论。否则差多少都有可能。具体参见http://sobereva.com/573
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
本团队长期招收研究生,有意者可私信联系

1万

帖子

0

威望

9859

eV
积分
22096

Level 6 (一方通行)

13#
发表于 Post on 2021-10-26 21:54:46 | 只看该作者 Only view this author
zzk 发表于 2021-10-26 13:40
那请问我这个三百个原子的分子如果想算溶质的自由能要采用什么方法呢?

如果一定要完全用高斯的话:用“讲究的人”的算法,但是把单点能的理论级别降低到比如说M06-2X/def2TZVP。注意有一点是你截图里没有明确写出来的,就是不管是结构优化还是单点能,如果用B3LYP必须加色散校正,如果用M06-2X则可加可不加。“讲究的人”的算法相比“本科生”算法只慢在单点能计算这一步,Shermo只费几分钟人力加几秒钟机时,没理由不用Shermo。假设你算反应自由能的话,计算结果误差在5kcal/mol量级。
但是建议下载一个orca,反正对学术界免费的,已经用高斯做的opt freq不需要重算,只用orca在高斯优化的结构的基础上算单点能,用PWPB95-D3/def2-TZVPP,开DLPNO加速,估计比你高斯做opt freq的时间要短,计算结果误差可以做到1.5kcal/mol量级。

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
zzk + 1 我很赞同

查看全部评分 View all ratings

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

管理员

公社社长

14#
发表于 Post on 2021-10-27 12:59:25 | 只看该作者 Only view this author
mfdsrax2 发表于 2021-10-26 20:52
请问一下高斯和ORCA算出来的能量差别大吗?

谈谈不同量子化学程序计算结果的差异问题
http://sobereva.com/573http://bbs.keinsci.com/thread-20122-1-1.html
北京科音自然科学研究中心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

27

帖子

0

威望

355

eV
积分
382

Level 3 能力者

15#
 楼主 Author| 发表于 Post on 2021-10-28 16:37:30 | 只看该作者 Only view this author
wzkchem5 发表于 2021-10-26 21:54
如果一定要完全用高斯的话:用“讲究的人”的算法,但是把单点能的理论级别降低到比如说M06-2X/def2TZVP ...

您好,请问DLPNO加速怎么开啊

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 06:08 , Processed in 0.217730 second(s), 26 queries , Gzip On.

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