24小时热门版块排行榜    

查看: 15859  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料080500调剂求收留 +7 一颗meteor 2026-03-13 7/350 2026-03-20 17:47 by 无懈可击111
[考研] 一志愿中南化学337求调剂 +4 niko- 2026-03-19 5/250 2026-03-20 13:49 by 促天成
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 14/700 2026-03-20 10:46 by EBSD
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 材料与化工求调剂 +7 为学666 2026-03-16 7/350 2026-03-19 14:48 by 尽舜尧1
[考研] 0817调剂 +3 没有答案_ 2026-03-14 3/150 2026-03-19 09:51 by Xu de nuo
[考研] 材料专业求调剂 +5 hanamiko 2026-03-18 5/250 2026-03-18 20:19 by 星空星月
[考研] 一志愿武理材料305分求调剂 +5 想上岸的鲤鱼 2026-03-18 6/300 2026-03-18 17:53 by 无际的草原
[考研] 297求调剂 +8 戏精丹丹丹 2026-03-17 8/400 2026-03-18 14:30 by laoshidan
[考研] 303求调剂 +4 睿08 2026-03-17 6/300 2026-03-18 11:01 by Iveryant
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 278求调剂 +5 烟火先于春 2026-03-17 5/250 2026-03-18 08:43 by 星空星月
[基金申请] 被我言中:新模板不强调格式了,假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 311求调剂 +3 冬十三 2026-03-13 3/150 2026-03-13 20:41 by JourneyLucky
信息提示
请填处理意见