Hazelnut 发表于 2023-11-28 17:27 好的，感谢。 |
朱陈 发表于 2023-11-28 17:18 高斯势只能设置宽度和高度，宽度的话不能太窄不然容易被随机波动克服，高度的话就根据自己的体系参考相关文献吧。 还可以看看vaspwiki的示例： https://www.vasp.at/wiki/index.php/PENALTYPOT |
Hazelnut 发表于 2023-11-28 17:11 这个位置好判断，但是这个宽度和斜率请问楼主是如何选择的呀？ |
朱陈 发表于 2023-11-28 17:12 是的，需要把希望受到势场作用的原子全部写上。 |
Hazelnut 发表于 2023-11-28 16:59 还有个问题请教你一下，在写ICONST的时候，是把全部的原子都写上才能在设定的区域内跑吧。 |
huxu 发表于 2023-11-23 17:32 不好意思，最近没看论坛回复。 ICONST中设置限制Z方向的方法是，Z (所需限制的原子序号) status 例如我要使用高斯形式的偏置势限制1号原子的话就是 Z 1 5 。 文件写起来感觉还是有点抽象的，如果要限制多个原子的话就得一个一个写。 然后再参照https://www.vasp.at/wiki/index.p ... _molecular_dynamics在INCAR或者PENALTYPOT文件里设置每个限制势的参数（比如势的位置、宽度、斜率等），所设置的参数的数量和顺序要和先前ICONST中相同。 |
Hazelnut 发表于 2023-11-28 16:59 谢谢你了，我的是5.4版本应该只能用高斯势函数了。 |
朱陈 发表于 2023-11-28 16:45 是的，但是需要额外在INCAR中设置FBIAS_A, FBIAS_D, 与FBIAS_R0参数，具体参考https://www.vasp.at/wiki/index.p ... _molecular_dynamics 但首先你需要确认你的vasp版本是否支持status=4，不然会在第一个离子步结束时报错。 |
Hazelnut 发表于 2023-11-22 21:20 你好，status=4是不是不需要设置PENALTYPOT文件了。 |
Hazelnut 发表于 2023-11-22 21:20 你好，请问该如何设置这个flag呢，比如是在Z方向相对于lattice vector 0.3添加高斯势函数，是直接写Z 0.3 5吗，我这么写会报错“Error reading ICONST (item 1 ): invalid definition of coordinate”，还请前辈指教！ |
|
小白先前的理解有误，flag = X, Y, and Z 是分开的，可以很方便的添加一个关于Z方向的偏置势。 折腾了两天发现踩了一个小坑，本人使用的vasp版本为5.4.4，按照vaspwiki上使用status = 4结果发现一直在第一个离子步后报参数缺失或有误的错，然后计算就直接终止了。又试着编译了vasp6.3.2版本，使用status = 4可以正常运行。 查阅了旧版本的vasp手册发现5.4.4版本应该是不支持status = 4的费米型阶跃函数（？），只能使用status = 5的高斯势函数。（vaspwiki没有写不支持，我还一直以为是输入文件格式有误，醉了）。 总之问题顺利解决了，感谢楼上大佬。 |
Aletyx 发表于 2023-11-21 09:41 感谢！已经解决了！ |
| 这个例子使用了PENALTYPOT：https://www.vasp.at/wiki/index.p ... tution_CH3Cl_-_mMD2 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2024-11-25 16:02 , Processed in 0.392393 second(s), 26 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites