sobereva 发表于 2023-4-12 18:15 感谢sob老师 |
Tiwmdh 发表于 2023-4-12 15:06 决定怎么算静电作用 不同Poission solver各有不同特点,必须选择合适的。几句话说不清楚,北京科音CP2K第一性原理培训班(http://www.keinsci.com/workshop/KFP_content.html)里很详细讲这些 SOLVER analytic强烈不建议使用,随真空区收敛巨慢 |
sobereva 发表于 2023-4-11 13:03 sob老师,我想再问一下CP2K中poisson的作用是什么,为啥同样的体系,poisson选项里写成PERIODIC NONE, SOLVER ANALYTIC整个构型就会散掉,而写成PERIODIC XYZ , SOLVER ANALYTIC就很稳定呢 |
本帖最后由 Tiwmdh 于 2023-4-11 14:52 编辑 sobereva 发表于 2023-4-11 13:03 图四给出的是一个Rh的台阶结构,在MD模拟中周期性边界条件,不需要固定原子都很稳定,原子总共有三层。现在想要从原始的团簇中切出来一个台阶结构,如图八中所示,然后发现跑MD过程中不稳定,图一是额外加了一层底层,并把加了的这层固定之后运行的结果,仍然不稳定,初始构型是图七,图二是认为可能是因为三层太少了所以切成了一个台运行的结果,初始构型是图三。图五和图六是另外一个构型的MD初末态图,同样出现了原子迁移构型不稳定的现象,现在知道原因了,感谢SOB老师 (我不知道咋弄的把图的顺序调乱了) |
|
1 看不明白问题和给出的图片的逻辑 保证有4层,并且计算时要把底端和四周边缘的原子都冻住 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-21 09:31 , Processed in 0.181656 second(s), 25 queries , Gzip On.