|
本帖最后由 乘风万里 于 2024-4-18 21:39 编辑 可以先用 tblite+ase 对分子晶体快速优化，https://github.com/tblite/tblite，最后用VASP优化。但如果对ase-api不熟悉，还是算了。 |
本帖最后由 2532251334 于 2024-3-9 10:11 编辑 logzzz 发表于 2024-3-5 10:23 老师我想问下您，如果先固定晶胞优化，那晶胞参数要怎么优化，先优化好结构在用isif=3优化晶胞和结构吗？我查了下6是只优化晶胞，这两个的优化顺序有啥建议吗？ |
喝杯热水 发表于 2024-3-8 14:03 好的多谢 |
2532251334 发表于 2024-3-8 10:09 我自己测试过，更少的NPAR反而具有更高的计算效率。 |
本帖最后由 2532251334 于 2024-3-8 10:12 编辑 喝杯热水 发表于 2024-3-6 11:07 您好，这个为什么要设置成2？能详细说下吗？ 我在超算上64个核，设置NCORE=16, KPAR=4，一个电子步90s;改成12,2，一个电子步只用40s，这是为啥？ |
| 设置一下NPAR=2，会加速几倍。 |
乘风万里 发表于 2024-3-5 20:31 谢谢，找到了，这样看起来边界地方还是会有好多不完整的 |
|
本帖最后由 乘风万里 于 2024-3-5 20:34 编辑 VESTA周期性显示问题点：Edit -> Bonds -> Search additional atoms if A1 is included in the boundary 给每个键选前两个都可以 
|
akakcolin 发表于 2024-3-5 15:36 谢谢，明白了。 |
卡开发发 发表于 2024-3-5 14:29 谢谢您，学到了很多，非常感谢！ |
|
楼上说的很好 注意到用的2018的intel编译器，可以无脑换最新intel编译器重新build，会快一些 另外一定要考虑色散 用dft算分子晶体一般情况下对计算精度是有要求的，用vasp大部分情况下是合适的 |
2532251334 发表于 2024-3-5 11:16 1、速度上的问题显然可以还有优化余地，ALGO用F(Fast亦即RMM-DIIS)替代N(Normal亦即Davidson)，我这边24核心（7R32）/48G内存大概最长的一个电子步是80s。另外楼上的建议是有效的，可以进行k点测试，假如仅用Gamma点就足够，那么用vasp_gam配合仅Gamma点的k点还能快出来不少。 2、初始结构也许不够好，这个时候可以适当用一些低等级的方法，如DFTB+或者xTB甚至是神经网络势程序粗略优化一下，另外POTIM对IBRION=2（亦即CG+LineSearch）在第一步有很大影响，一般来说可以取个稍微小一些的值确保结构不要炸掉。 3、一个是你的结构本身优化后不是很好，一个是vesta的显示模式并不是很好，例如边缘原子和成键显示的处理方式不是特别得当。 另： 1、可以使用一些LCAO程序，但不是非要用哪种程序就行。除了上述的CP2K，你其实也可以试试看OpenMX，但是功能不是很多，但正因为如此一共也没有几个参数，如果对计算原理比较熟悉的情况基本上一两天就能掌握操作。 2、需要特别指出的是，平面波确实一般情况下计算要略慢于LCAO，但在精度稍高的赝势支持下只要动能截断收敛就能拿到足够高精度的解，衡量计算精度不少人提出了不少测试的方案，其中一种知名的方案叫做delta测试，上述推荐OpenMX也是因为至少在Standard这个级别下OpenMX的精度在delta测试中表现甚至不会比平面波精度差很多，而速度上我们这边测试24核心情况下基本上计算时间大概在20s每电子步。 3、具体要使用什么程序取决于你想要通过计算解决什么问题，如果你只是优化一下晶体结构，很可能力场就做的不错，大可不必费时费力用一个还要看DFT和各个基组、赝势色散修正方式的精度脸色这么个方案。 |
| 参与人数Participants 1 | eV +2 | 收起 理由Reason |
|---|---|---|
|
| + 2 |
sobereva 发表于 2024-3-5 11:43 明白了，谢谢sob老师 |
2532251334 发表于 2024-3-5 10:59 分清楚内存和硬盘 固态还是机械无所谓。CP2K计算速度跟硬盘速度无关 哪个分区有地方客户机就装哪 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2024-11-24 23:00 , Processed in 0.177509 second(s), 26 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites