24小时热门版块排行榜    

CyRhmU.jpeg
查看: 430  |  回复: 0
【悬赏金币】回答本帖问题,作者qisames将赠送您 50 个金币

qisames

新虫 (初入文坛)

[求助] 求助,如何指定体系的密度(ws元胞的半径)

各位老师好,
我现在想要用vasp重复一论文中的md-dft结果,
论文采用的是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里把 调整原子之间的距离?
有没有相应的参数可以做呀

比如说,我想计算如下体系:
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




以下是poscar:
h
1.0
       10.6359996796         0.0000000000         0.0000000000
       -5.3179998398         9.2110459171         0.0000000000
        0.0000000000         0.0000000000        17.3360004425
    h
  128
direct
     0.083333336         0.166666672         0.062500000
     0.916666687         0.833333313         0.937500000
     0.833333313         0.916666687         0.062500000
     0.166666672         0.083333336         0.937500000
     0.083333336         0.916666687         0.062500000
     0.916666687         0.083333336         0.937500000
     0.916666687         0.833333313         0.187500000
     0.166666672         0.083333336         0.187500000
     0.916666687         0.083333336         0.187500000
     0.083333336         0.166666672         0.312500000
     0.916666687         0.833333313         0.687500000
     0.833333313         0.916666687         0.312500000
     0.166666672         0.083333336         0.687500000
     0.083333336         0.916666687         0.312500000
     0.916666687         0.083333336         0.687500000
     0.916666687         0.833333313         0.437500000
     0.083333336         0.166666672         0.812500000
     0.166666672         0.083333336         0.437500000
     0.833333313         0.916666687         0.812500000
     0.916666687         0.083333336         0.437500000
     0.083333336         0.916666687         0.812500000
     0.083333336         0.166666672         0.562500000
     0.833333313         0.916666687         0.562500000
     0.083333336         0.916666687         0.562500000
     0.083333336         0.416666687         0.062500000
     0.916666687         0.583333313         0.937500000
     0.583333313         0.666666687         0.062500000
     0.416666687         0.333333343         0.937500000
     0.333333343         0.916666687         0.062500000
     0.666666687         0.083333336         0.937500000
     0.916666687         0.583333313         0.187500000
     0.416666687         0.333333343         0.187500000
     0.666666687         0.083333336         0.187500000
     0.416666687         0.083333336         0.937500000
     0.583333313         0.916666687         0.062500000
     0.666666687         0.583333313         0.937500000
     0.333333343         0.416666687         0.062500000
     0.916666687         0.333333343         0.937500000
     0.083333336         0.666666687         0.062500000
     0.416666687         0.083333336         0.187500000
     0.666666687         0.583333313         0.187500000
     0.916666687         0.333333343         0.187500000
     0.083333336         0.416666687         0.312500000
     0.916666687         0.583333313         0.687500000
     0.583333313         0.666666687         0.312500000
     0.416666687         0.333333343         0.687500000
     0.333333343         0.916666687         0.312500000
     0.666666687         0.083333336         0.687500000
     0.916666687         0.583333313         0.437500000
     0.083333336         0.416666687         0.812500000
     0.416666687         0.333333343         0.437500000
     0.583333313         0.666666687         0.812500000
     0.666666687         0.083333336         0.437500000
     0.333333343         0.916666687         0.812500000
     0.416666687         0.083333336         0.687500000
     0.583333313         0.916666687         0.312500000
     0.666666687         0.583333313         0.687500000
     0.333333343         0.416666687         0.312500000
     0.916666687         0.333333343         0.687500000
     0.083333336         0.666666687         0.312500000
     0.583333313         0.916666687         0.812500000
     0.416666687         0.083333336         0.437500000
     0.333333343         0.416666687         0.812500000
     0.666666687         0.583333313         0.437500000
     0.083333336         0.666666687         0.812500000
     0.916666687         0.333333343         0.437500000
     0.083333336         0.416666687         0.562500000
     0.583333313         0.666666687         0.562500000
     0.333333343         0.916666687         0.562500000
     0.583333313         0.916666687         0.562500000
     0.333333343         0.416666687         0.562500000
     0.083333336         0.666666687         0.562500000
     0.333333343         0.166666672         0.062500000
     0.666666687         0.833333313         0.937500000
     0.833333313         0.166666657         0.062500000
     0.166666672         0.833333313         0.937500000
     0.833333313         0.666666687         0.062500000
     0.166666657         0.333333343         0.937500000
     0.666666687         0.833333313         0.187500000
     0.166666672         0.833333313         0.187500000
     0.166666657         0.333333343         0.187500000
     0.333333343         0.166666672         0.312500000
     0.666666687         0.833333313         0.687500000
     0.833333313         0.166666657         0.312500000
     0.166666672         0.833333313         0.687500000
     0.833333313         0.666666687         0.312500000
     0.166666657         0.333333343         0.687500000
     0.666666687         0.833333313         0.437500000
     0.333333343         0.166666672         0.812500000
     0.166666672         0.833333313         0.437500000
     0.833333313         0.166666657         0.812500000
     0.166666657         0.333333343         0.437500000
     0.833333313         0.666666687         0.812500000
     0.333333343         0.166666672         0.562500000
     0.833333313         0.166666657         0.562500000
     0.833333313         0.666666687         0.562500000
     0.333333343         0.666666687         0.062500000
     0.666666687         0.333333313         0.937500000
     0.666666687         0.333333313         0.187500000
     0.333333343         0.666666687         0.312500000
     0.666666687         0.333333313         0.687500000
     0.666666687         0.333333313         0.437500000
     0.333333343         0.666666687         0.812500000
     0.333333343         0.666666687         0.562500000
     0.583333313         0.166666672         0.062500000
     0.416666657         0.833333313         0.937500000
     0.833333313         0.416666657         0.062500000
     0.166666672         0.583333313         0.937500000
     0.583333313         0.416666657         0.062500000
     0.416666657         0.583333313         0.937500000
     0.416666657         0.833333313         0.187500000
     0.166666672         0.583333313         0.187500000
     0.416666657         0.583333313         0.187500000
     0.583333313         0.166666672         0.312500000
     0.416666657         0.833333313         0.687500000
     0.833333313         0.416666657         0.312500000
     0.166666672         0.583333313         0.687500000
     0.583333313         0.416666657         0.312500000
     0.416666657         0.583333313         0.687500000
     0.416666657         0.833333313         0.437500000
     0.583333313         0.166666672         0.812500000
     0.166666672         0.583333313         0.437500000
     0.833333313         0.416666657         0.812500000
     0.416666657         0.583333313         0.437500000
     0.583333313         0.416666657         0.812500000
     0.583333313         0.166666672         0.562500000
     0.833333313         0.416666657         0.562500000
     0.583333313         0.416666657         0.562500000
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 qisames 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
信息提示
请填处理意见