|
发布时间: 2025-12-12 10:37
正文摘要:本帖最后由 MVK 于 2025-12-13 21:02 编辑 各位老师好，我目前在做一个沸石负载 Pt–Sn 金属团簇催化丙烷脱氢（PDH）的计算课题。载体为沸石骨架，目前使用的模型包含 576 个 Si/O 原子，13 个Pt/S ... |
MVK 发表于 2025-12-14 09:56 其实，你可以对现在最高值做一次频率计算了，如果显示只有一个虚频，而且把虚频放在JMOL里面观察是跟预期断键方向一致就说明成功了。CINEB控制精度是K_SPRING，我是习惯一开始0.08跑一次，然后0.05、0.02再跑一次。粗收敛的时候multiwfn默认的收敛限往上调一些 |
小冷~715 发表于 2025-12-14 00:03 不好意思老师昨晚看错日期了 实际只有不到1/10 今早检查发现有一个跑满300帧达到步数上限结束了但是看趋势是对的 请问您之前说的“粗收敛取最高点附近的三个结构再插点”里的粗收敛是指使用更低的收敛限吗 如果是的话调整到多少合适呀 我看目前的收敛限为MAX_DR 2E-4 MAX_FORCE 4.5E-4 RMS_DR 1E-4 RMS_FORCE 3E-4 |
小冷~715 发表于 2025-12-14 00:03 不会吧 我用的就是并行超算云 这个应该用的人很多吧 而且上面贴了参考文献的机理图 我参考着他的机理来算的（10.1038/s41586-025-09168-8）还是篇nature |
MVK 发表于 2025-12-13 22:49 1天两千多也太夸张了，这是什么黑心超算。我之前是指每个副本安排64核，八个副本就是512核，4个128核节点，这样可能会帮你赶在ddl前弄完。过渡态这东西里面经验以及化学直觉占比很大，你如果初末态相差很离谱，跑很久都收敛不了的。你跟着社长那个cp2k班学一下，收获会很大。而且我印象中PDH不是单纯一个丙烷脱氢吧，它还有积碳啥的，说明这东西并不是简单的一级解离反应。多参考一下实验动力学和跟一些经典量化文献，后面摆结构会更合理一些。 |
sobereva 发表于 2025-12-13 21:38 老师我感觉自己的设置也没有问题呀 这个小体系怎么会算的这么慢呢 |
本帖最后由 MVK 于 2025-12-13 22:57 编辑 小冷~715 发表于 2025-12-13 21:43 唉是的 我现在就在用超算算 目前是128个核去算8个结构跑的NEB 今天只跑了两个就烧了快2000RMB了 而且还没有结果 心态快炸了 ddl只有一周了12个过渡态一个都没找到 感觉只能放弃了  |
MVK 发表于 2025-12-13 20:52 我估计你是返稿吧，有ddl。一般这种时间紧的过渡态都是用超算，一次调用几百核。你之前是用64核算一个结构优化，然后现在inp里面写的是用64核算8个结构，每个结构是分到8个核，效率肯定不一样。还好你结构小，不然就爆内存了。MFI在十元环沸石里面算小的了，近几年的十元环沸石单胞就有700到1000个原子。相同原子数条件下，我觉得沸石比过渡金属氧化物这种算的快多了，沸石也容易收敛。 |
MVK 发表于 2025-12-13 20:58 区区303原子，当然是对CP2K而言的小体系 用的方式不对自然会慢。不仅是计算级别，还注意监控计算过程/收敛情况、image数和计算用的副本数的设置。所有计算要点在北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）都做了很充分的讲解，没学过的话必须系统性学一遍，否则踩了坑很多坑都浑然不知，若已经学过的话好好复习一遍 最快的CP2K都算不动的就没有任何其它程序能算得动了 |
本帖最后由 MVK 于 2025-12-13 21:03 编辑 sobereva 发表于 2025-12-13 04:28 老师这个体系真的小吗 我用切了一半的沸石团簇来计算 replica设置成8跑NEB仍然奇慢 今天128核算了快十个小时了还没算完 虽然大家都说cp2k跑的快但是自己也没有用过其他程序...这样的速度在第一性原理程序里是算快的吗 还是说我哪里还是设置的不对QAQ(上面又补充了下目前的输入文件) |
本帖最后由 MVK 于 2025-12-13 21:02 编辑 小冷~715 发表于 2025-12-13 14:38 老师我用一半大小的沸石团簇128核replica设置成8 泛函基组设置为PBE-D3(BJ)/DZVP-MOLOPT-SR-GTH开启OT计算仍然是奇慢无比 今天提交后算了快九个小时还没算完 请问是我哪里设置的不对吗(我在上面又贴了下现在的输入文件麻烦老师您帮忙看看吧我实在是没招了QAQ) 而且请问这里每个插点分配一个节点或者半个节点怎么做到呀 |
| 可以用机器学习势试一试。 |
| 也是做沸石+烷烃脱氢/CO2加氢方向。我觉得算的慢最大问题在于你插点太多，而且只用了一个双路节点。我之前计算是每个插点分配一个节点或者半个节点。CINEB的目的就是减少插点的数量，你这种前后变化不大的没必要插那么多。第一次插点6-8个，然后粗收敛取最高点附近的三个结构再插点，提高精度收敛到TS。IS和FS用第一次插点收敛的两个前后最低点做进一步结构优化。cp2k算沸石很快的 |
|
2 可以混用不同基组，包括全电子和赝势基组。优化目的DZVP-MOLOPT-SR-GTH性价比最高，精度也够，不需要考虑别的 CP2K支持不同理论方法描述的区域拼接。但对你这种很小的体系没必要折腾 |
我先自问自答一波 用multiwfn主功能300选项7的几何处理功能把结构调整了下沸石切掉了一半 第一个问题解决了 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-1-24 04:22 , Processed in 0.183923 second(s), 26 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites