求助，如何指定体系的密度(ws元胞的半径）

qisames

各位老师好，
我现在想要用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