计算化学公社

标题: 求助，如何指定体系的密度(ws元胞的半径） [打印本页]

作者
Author: qisames 时间: 2022-7-21 15:55
标题: 求助，如何指定体系的密度(ws元胞的半径）
本帖最后由 qisames 于 2022-7-21 16:01 编辑

各位老师好，
我现在想要用vasp重复一论文中的 DFT-MD 的运算，
论文采用的是CPMD package.
NVT 系综。
cutoff enery 35-50 ha
perdew-burke-ernzerof gga approximation
ionic time step  16 a.u.( 1 a.u.=0.0242fs)
分别用不同温度和 wigner-seitz radius.  来进行计算
温度，我知道应是用  tebeg = 2000
但是  wigner-seitz radius 也就是密度，应该何约束呢，

是不是只能在POSCAR里调整原子之间的距离？
还有有什么其它的方法？or 参数？

比如说，我想计算如下体系：
128个H原子的体系。用vesta做了超晶胞
温度2000k。 tebeg = 2000
正则系综，nose-hoover_thermostat    mdalgo=2
rs=1.6波尔半径  (0.66g/立方厘米)

以下是我的INCAR ：

ismear = 0
ibrion = 0
mdalgo = 2
isif = 2
smass = 1.0
sigma = 0.1
lreal = auto
encut = 1400
isym = 0
tebeg = 2000
nsw = 600
potim = 0.3872

欢迎光临计算化学公社 (http://bbs.keinsci.com/)