24小时热门版块排行榜    

查看: 15586  |  回复: 314

m1520342646

金虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
各位我是新手,想利用VASP计算氢分子的键长( 0.75102 A),程序如下:
INCAR
SYSTEM = H2
ISMEAR = 0
NSW = 0
IBRION = 2
ENCUT = 350.0
KPOINTS
Auto
0
Monk
1 1 1
0 0 0
loop
#!/bin/sh
for i in 0.72 0.74 0.76 0.78 0.80
do
cat >POSCAR < SYSTEM = H2
$i
1 0   0
0 20  0
0 0   20
2
direct
0 0 0
1 0 0
!
vasp-p
E=`grep "TOTEM" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'`
#V=`grep "volume" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'`
echo "E=$E"
echo $i $E >>a-e.dat
done
rm CHG* WAVECAR
但是运行后,有错误,说元胞体积为0,
能帮我修改一下,或者写出正确的程序?谢谢~~~
241楼2013-03-04 15:38:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞

引用回帖:
19936757楼: Originally posted by m1520342646 at 2013-03-04 15:38:26
各位我是新手,想利用VASP计算氢分子的键长( 0.75102 A),程序如下:
INCAR
SYSTEM = H2
ISMEAR = 0
NSW = 0
IBRION = 2
ENCUT = 350.0
KPOINTS
Auto
0
Monk
1 1 1
0 0 0
loop
#!/bin/sh
for i i ...

#!/bin/sh
for i in 0.036 0.037 0.038 0.039 0.04
do
cat >POSCAR < SYSTEM = H2
1.0
20 0   0
0 20  0
0 0   20
2
direct
0 0 0
$i 0 0
!
vasp-p
E=`grep "TOTEM" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'`
#V=`grep "volume" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'`
echo "E=$E"
echo $i $E >>a-e.dat
done
rm CHG* WAVECAR
真理是一点点接近的!
242楼2013-03-04 19:32:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
243楼2013-03-26 11:14:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuchangzzw

铜虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
我想问下vasp计算出的DOSCAR里的没列数据代表什么。 我分别用LORBIT = 10 和11计算了一种材料的DOSCAR
结果发现LORBIT = 10 时DOSCAR后面有这些数据
12.81861797    -26.09640211  301      7.92181128      1.00000000
    -26.096  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
    -25.967  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
    -25.837  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
    -25.707  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00

而LORBIT = 11时多了许多列数据:
     13.32148466    -26.09927272  301      8.02304946      1.00000000
    -26.099  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E
+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0
000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
0.0000E+00
    -25.968  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E
+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0
000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
0.0000E+00
我想分别画出dxy dyz dzx等5个分轨道的态密度,我该如何选择数据?
244楼2013-03-27 10:26:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞

引用回帖:
19937249楼: Originally posted by liuchangzzw at 2013-03-27 10:26:57
我想问下vasp计算出的DOSCAR里的没列数据代表什么。 我分别用LORBIT = 10 和11计算了一种材料的DOSCAR
结果发现LORBIT = 10 时DOSCAR后面有这些数据
12.81861797    -26.09640211  301      7.92181128      1.00 ...

不加LORBIT是不画分立的态密度。
LORBIT=10 是 s p d f的态密度。
LORBIT=11  s px py pz   dxy ...
真理是一点点接近的!
245楼2013-03-28 09:51:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuchangzzw

铜虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
19937285楼: Originally posted by gavinliu7390 at 2013-03-28 09:51:08
不加LORBIT是不画分立的态密度。
LORBIT=10 是 s p d f的态密度。
LORBIT=11  s px py pz   dxy ......

谢谢回答,但我就是想画出 PX PY PZ DXY DYZ DZX……这些分轨道的态密度,用LORBIT=11算出的DOSCAR里面好多列,不知道哪一列才是啊。-26.099  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E
+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0
000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
0.0000E+00
光一个能量点的态密度就有这么多数据
246楼2013-04-01 11:55:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞

引用回帖:
19937359楼: Originally posted by liuchangzzw at 2013-04-01 11:55:47
谢谢回答,但我就是想画出 PX PY PZ DXY DYZ DZX……这些分轨道的态密度,用LORBIT=11算出的DOSCAR里面好多列,不知道哪一列才是啊。-26.099  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.00 ...

每一列依次是能量,s,px,py,pz,Dxy,Dyz,Dzx ...
真理是一点点接近的!
247楼2013-04-01 19:07:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunray55

银虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
有用的帖子,

[ Last edited by sunray55 on 2013-4-28 at 00:47 ]
248楼2013-04-27 00:11:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzzaqxq

银虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
你好,我想问一下,单机能安装vasp吗?win7系统的!想练习练习,谢谢
勤能补拙是良训一分辛苦一分才
249楼2013-04-30 19:08:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞

引用回帖:
249楼: Originally posted by zzzaqxq at 2013-04-30 19:08:30
你好,我想问一下,单机能安装vasp吗?win7系统的!想练习练习,谢谢

win7不能,必须是linux,unix,或者mac。
真理是一点点接近的!
250楼2013-07-30 23:00:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gavinliu7390 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见