应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 其他 » 【其他】在线答疑:VASP的问题
31525/32返回列表首页上一页232425262728下一页末页
查看: 16483  |  回复: 314
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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 的主题更新
31525/32返回列表首页上一页232425262728下一页末页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于filecode，很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:22 by oEVWOejN9taj
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 16:47 by oEVWOejN9taj
[论文投稿] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 16:47 by oEVWOejN9taj
[考博] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:27 by oEVWOejN9taj
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:22 by oEVWOejN9taj
[基金申请] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:19 by oEVWOejN9taj
[博后之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:02 by oEVWOejN9taj
[论文投稿] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 15:07 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 WQyTGMbfH7Cx 2026-08-07 4/200 2026-08-08 15:02 by oEVWOejN9taj
[基金申请] fileCode有新解读？ +9 Tide man 2026-08-08 17/850 2026-08-08 14:56 by archvillain
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:27 by oEVWOejN9taj
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 好消息？这个有何含义？？？ +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[基金申请] 面上提前没消息，有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
信息提示
关闭
请填处理意见
关闭