zzk 发表于 2021-10-28 09:37 %mp2 dlpno true end 具体看orca说明书的双杂化泛函章节 |
zzk 发表于 2021-10-28 16:51 看卢老师博文《详谈Multiwfn产生ORCA量子化学程序的输入文件的功能》http://bbs.keinsci.com/thread-13560-1-1.html,在里面搜DLPNO。 |
zzk 发表于 2021-10-28 16:37 只需要在!后加DLPNO-CCSD(T) 就可以吗 |
wzkchem5 发表于 2021-10-26 21:54 您好,请问DLPNO加速怎么开啊 |
mfdsrax2 发表于 2021-10-26 20:52 谈谈不同量子化学程序计算结果的差异问题 http://sobereva.com/573(http://bbs.keinsci.com/thread-20122-1-1.html) |
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量级。 |
| 参与人数Participants 1 | eV +1 | 收起 理由Reason |
|---|---|---|
|
| + 1 | 我很赞同 |
mfdsrax2 发表于 2021-10-26 13:52 补充一下,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 |
本帖最后由 zjxitcc 于 2021-10-26 21:25 编辑 mfdsrax2 发表于 2021-10-26 20:52 取决于你的关键词如何写;以及你想问绝对能量还是相对能量。 绝对能量:在合理关键词的情况下,例如气相下的HF, MP2, CCSD等可以严格几乎一样。各种泛函那影响因素就多了,积分格点、泛函定义的细节、SCF收敛限、积分精度、隐式溶剂模型的差异等等都有影响。所以这种问题意义不大。 相对能量(例如反应能垒):在合理关键词的情况下,所有软件算出来都趋势一致。 |
sobereva 发表于 2021-10-25 11:28 请问一下高斯和ORCA算出来的能量差别大吗? |
| 只能用第一种无法高精度计算那个吗 |
| Gaussian做MP2/双杂化结合def2-TZVPP这种档次的基组计算,50原子以上就根本甭考虑。再多就用ORCA |
| 以前在G16上试过 B2PLYPD3/def2tzvp 有机分子 100个左右的原子 2000个基函数左右,需要大概8TB的硬盘空间。(计算给了48核心180GB内存,单点能算了24小时)因为内存不够基本只达到了最低限(2OVN), 硬盘就会占的很多。高斯官方给的是MP2 Direct算法下硬盘空间跟基函数的平方成正比。那么我觉得你的体系30TB不够也属于正常现象。 |
| mp2不用density fitting/ri局限性太大了 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-21 07:37 , Processed in 0.177745 second(s), 26 queries , Gzip On.