计算化学公社

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

[新手求助] 关于大体系S0优化的疑问

[复制链接 Copy URL]

430

帖子

0

威望

784

eV
积分
1214

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
最近做了一个较大体系的结构(300原子),用6-31g能跑成功,但是用了6-31g*之后,计算中的临时文件太大,导致根目录磁盘不够。
想问一下6-31g*优化S0时,我先算opt,然后进行freq计算的结果和两个一起进行计算有区别么?这样也可降低每次计算的任务量,谢谢!

3097

帖子

29

威望

1万

eV
积分
17221

Level 6 (一方通行)

2#
发表于 Post on 2017-11-23 09:27:09 | 只看该作者 Only view this author
opt freq 并不是“一起计算”,只是“自动依次计算”而已。
自己算opt再做freq不会让计算耗时更低(写的不好还能更高)
300个原子的优化一般机器会非常费劲,尽可能的想办法简化体系、或者用其他的处理方式。

430

帖子

0

威望

784

eV
积分
1214

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2017-11-23 11:25:35 | 只看该作者 Only view this author
liyuanhe211 发表于 2017-11-23 09:27
opt freq 并不是“一起计算”,只是“自动依次计算”而已。
自己算opt再做freq不会让计算耗时更低(写的不 ...

已经是简化过的体系了,应该能算的,只不过中间文件太大,超过我的磁盘了。
谢谢

221

帖子

0

威望

6218

eV
积分
6439

Level 6 (一方通行)

跳跳猪

4#
发表于 Post on 2017-11-24 00:44:58 | 只看该作者 Only view this author
1. 磁盘根目录不够?其它目录呢,再不行外接个硬盘呢。
2. 这么大体系,感觉不如用半经验。(http://sobereva.com/214
流年似水,浮生如梦。

430

帖子

0

威望

784

eV
积分
1214

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2017-12-6 12:59:50 | 只看该作者 Only view this author
Mikasa 发表于 2017-11-24 00:44
1. 磁盘根目录不够?其它目录呢,再不行外接个硬盘呢。
2. 这么大体系,感觉不如用半经验。(http://sober ...

文献也是用DFT进行计算的,半经验不具有说服力。原因是加了弥散函数后,freq计算中间文件太大,占满根目录了。

221

帖子

0

威望

6218

eV
积分
6439

Level 6 (一方通行)

跳跳猪

6#
发表于 Post on 2017-12-7 00:42:37 | 只看该作者 Only view this author
yezhonghua 发表于 2017-12-6 12:59
文献也是用DFT进行计算的,半经验不具有说服力。原因是加了弥散函数后,freq计算中间文件太大,占满根目 ...

你一直说的是“占满根目录了”
你的意思是整个硬盘没有空间了,还是只是根目录的空间不够?
流年似水,浮生如梦。

6万

帖子

99

威望

6万

eV
积分
125151

管理员

公社社长

7#
发表于 Post on 2017-12-7 17:54:02 | 只看该作者 Only view this author
yezhonghua 发表于 2017-12-6 12:59
文献也是用DFT进行计算的,半经验不具有说服力。原因是加了弥散函数后,freq计算中间文件太大,占满根目 ...

*是极化,不是弥散
如果GAUSS_SCRDIR指定的目录写满了,而其它目录还有空间,可以分割rwf。如果还不够干脆买个4T硬盘,也没多少钱。
300原子做DFT/6-31G*的振动分析是非常困难的。
北京科音自然科学研究中心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

430

帖子

0

威望

784

eV
积分
1214

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2017-12-8 08:46:33 | 只看该作者 Only view this author
sobereva 发表于 2017-12-7 17:54
*是极化,不是弥散
如果GAUSS_SCRDIR指定的目录写满了,而其它目录还有空间,可以分割rwf。如果还不够干 ...

谢谢sober老师,那如果采用G16在freq分析上是不是能够快很多相对于G09。
因为G16对于TD-DFT的freq采用了新的解析方法,比原来的G09快很多,特别是对于激发态的freq优化。请教一下G16对于基态的 DFT freq计算也能有很快的加速么?
谢谢

6万

帖子

99

威望

6万

eV
积分
125151

管理员

公社社长

9#
发表于 Post on 2017-12-8 13:02:05 | 只看该作者 Only view this author
yezhonghua 发表于 2017-12-8 08:46
谢谢sober老师,那如果采用G16在freq分析上是不是能够快很多相对于G09。
因为G16对于TD-DFT的freq采用了 ...


对于非TDDFT/TDHF,freq速度基本没有变化。而且由于默认用了int=ultrafine,速度还更慢。(但如果是比如AVX2版G16和非AVX2的G09比,在支持AVX2的CPU上,相同设定下,G16会快一点。这决不是改进算法的功劳)
北京科音自然科学研究中心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

430

帖子

0

威望

784

eV
积分
1214

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2017-12-8 16:25:37 | 只看该作者 Only view this author
sobereva 发表于 2017-12-8 13:02
对于非TDDFT/TDHF,freq速度基本没有变化。而且由于默认用了int=ultrafine,速度还更慢。(但如果是比 ...

有人做个过一个64个原子体系的TD-B3LYP/CCPVDZ S1 freq 优化,16核心CPU。G16用了2小时,G09用了51小时,速度提升很明显。说是G16的TD用了analysis 分析比G09的数值分析快了不少。不知道对于原子数量的增加,G16对于S1 freq的提升速度相比G09是否更快。
谢谢!

6万

帖子

99

威望

6万

eV
积分
125151

管理员

公社社长

11#
发表于 Post on 2017-12-8 20:23:27 | 只看该作者 Only view this author
yezhonghua 发表于 2017-12-8 16:25
有人做个过一个64个原子体系的TD-B3LYP/CCPVDZ S1 freq 优化,16核心CPU。G16用了2小时,G09用了51小时, ...

那是因为G16支持TD的解析Hessian,G09只支持半数值Hessian
体系越大提升速度越多
北京科音自然科学研究中心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

430

帖子

0

威望

784

eV
积分
1214

Level 4 (黑子)

12#
 楼主 Author| 发表于 Post on 2017-12-9 22:03:32 | 只看该作者 Only view this author
sobereva 发表于 2017-12-8 20:23
那是因为G16支持TD的解析Hessian,G09只支持半数值Hessian
体系越大提升速度越多

多谢sober老师,非常感谢您的解答!

本版积分规则 Credits rule

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

GMT+8, 2026-2-23 08:34 , Processed in 0.159077 second(s), 20 queries , Gzip On.

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