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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 其他 » 【其他】在线答疑:VASP的问题
31525/32返回列表首页上一页232425262728下一页末页
查看: 15954  |  回复: 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下一页末页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 290调剂生物0860 +13 哇哈哈，。 2026-04-11 14/700 2026-04-11 17:28 by grlsj
[考研] 270求调剂 +10 杨乐369 2026-04-11 10/500 2026-04-11 16:18 by 明月此时有
[考研] 295求调剂 +3 桂秋二十八 2026-04-05 5/250 2026-04-11 11:36 by zhq0425
[考研] 0854调剂 +5 音像店听花鼓戏 2026-04-10 5/250 2026-04-11 10:49 by qingpingzhu
[考研] 还有化工二轮调剂的学校吗 5+14 化工人999 2026-04-09 48/2400 2026-04-11 10:27 by 89436494
[考研] 284求调剂 +12 archer.. 2026-04-10 13/650 2026-04-11 08:44 by zhq0425
[考研] 一志愿华南理工大学331分材料求调剂 +9 天下ww 2026-04-09 9/450 2026-04-10 22:58 by Ftglcn90
[考研] 284求调剂 +9 让我上岸吧阿西 2026-04-09 11/550 2026-04-10 19:18 by 靖jing
[考研] 一志愿0703化学招61最终排名62化学求调剂 +24 招61排名62 2026-04-07 28/1400 2026-04-10 16:15 by yx54321
[考研] 求调剂 +11 翩翩一书生 2026-04-09 13/650 2026-04-10 10:27 by liuhuiying09
[考研] 08工学 309分求调剂 +6 Yin DY 2026-04-08 6/300 2026-04-10 09:18 by Delta2012
[考研] 085601初试330分找调剂 +10 流心奶黄包l 2026-04-09 10/500 2026-04-10 08:14 by Sammy2
[考研] 材料专硕调剂 +16 哈哈哈吼吼吼哈 2026-04-07 17/850 2026-04-09 21:16 by wutongshun
[考研] 322求调剂，08工科 +3 今天是个小号 2026-04-08 3/150 2026-04-09 15:53 by wp06
[考研] 311求调剂 +6 surte 2026-04-08 13/650 2026-04-09 14:00 by surte
[考研] 286求调剂 +19 Faune 2026-04-08 20/1000 2026-04-09 08:36 by 哦哦123
[考研] 388求调剂 +6 四川王涛 2026-04-07 8/400 2026-04-08 00:17 by JourneyLucky
[考研] 085405软件工程301分求调剂，专硕可跨专业，四六级已过 +3 静静想想 2026-04-05 3/150 2026-04-06 15:23 by nepu_uu
[考研] 332求调剂 +17 小小孟... 2026-04-05 18/900 2026-04-06 09:51 by 蓝云思雨
[考研] 08专硕275调剂 +5 AaAa7420 2026-04-05 5/250 2026-04-05 18:01 by jkddd
信息提示
关闭
请填处理意见
关闭