24小时热门版块排行榜

返回列表

当前主题已经存档。

limaojlu

金虫 (正式写手)

应助: 8 (幼儿园)
金币: 618
散金: 690
红花: 9
帖子: 426
在线: 65.2小时
虫号: 554763
注册: 2008-05-08
性别: GG
专业: 半导体物理

[交流] 【求助】用VASP做MD，能将每步的速度信息都给出吗？

请教大家了，我现在由于要看一个MD进行到不同程度，然后冷却的效果，所以想做一个较长的MD，然后取他们在不同程度的构型和速度接着做冷却，但是VASP，我找说明找了一遍，包括NWRITE的书写，和VASPRUN.XML,和XDATACAR的说明都没有找到怎么把每步的速度信息给出，谢谢大家了。

回复此楼

» 猜你喜欢

申博求助已经有4人回复
27年申博自荐-能源/环境方向已经有3人回复
闲聊已经有3人回复
亲密关系已经有13人回复
27秋申博自荐已经有4人回复

1楼 2009-12-21 12:37:16

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

xbelinda

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 391.7
红花: 1
帖子: 160
在线: 6.4小时
虫号: 587329
注册: 2008-08-22
性别: MM
专业: 无机合成和制备化学

★ ★ ★ ★ ★
limaojlu(金币+3,VIP+0): 12-21 13:20
qasd(金币+2,VIP+0):xiexie~ 12-21 19:43

呵呵，是这样。
用vasp做md的时候，每一步的每个离子速度都写在CONTCAR里。
在CONTCAR里，先是晶格坐标，然后是离子坐标，再下面就是速度了。
OUTCAR里也应该有。回头我再看看！

赞一下(13人)

回复此楼

2楼2009-12-21 13:05:09

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

limaojlu

金虫 (正式写手)

应助: 8 (幼儿园)
金币: 618
散金: 690
红花: 9
帖子: 426
在线: 65.2小时
虫号: 554763
注册: 2008-05-08
性别: GG
专业: 半导体物理

CONCAR只实时的给出某一步的信息，所以计算完后就只有最后一步了，不过这倒是一个解决的方案，我用
tail -f CONTCAR 来看，OUTCAR好像是没有速度的信息。谢谢答复

赞一下

回复此楼

3楼2009-12-21 13:20:28

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

xbelinda

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 391.7
红花: 1
帖子: 160
在线: 6.4小时
虫号: 587329
注册: 2008-08-22
性别: MM
专业: 无机合成和制备化学

★ ★ ★ ★ ★ ★
ice_rain(金币+1,VIP+0):谢谢交流 12-21 17:50
limaojlu(金币+5,VIP+0): 12-21 18:19

恩OUTCAR里的确没有速度。
我能想到的解决办法有两个，不改vasp代码的情况下，你可以写一个脚本。
你在INCAR中把NSW=50 或100 (多少步你想看一次速度) 每一步都看速度是没有什么太大意义的,因为分子动力学是用系综的,它正确的前提是各态历经,热力学量对时间的平均才是它的宏观物理量.
最好把INCAR中的LWAVE=T 打开.
这样,没做50步,你就把CONTCAR copy一下.然后继续跑,由于有波函数和离子位型和速度.不会影响你的结果.
第二个方法,我觉得比较好,就是你直接改源代码. 例如:你在main.F的开始写
open(1002,file="v.dat"

在main.F的结尾写
close(1002)
然后你找到poscar.F
在OUTPOS_TRAIL的子程序里
VTMP就是速度.
你可以在WRITE(IU,480)VTMP(1),VTMP(2),VTMP(3)
下面加上
write(1002,480)VTMP(1),VTMP(2),VTMP(3)

然后,编译vasp.
运行后,会产生v.dat,每一次的速度都会写出.
方法不唯一, 不过都类似.

赞一下(8人)

回复此楼

4楼2009-12-21 14:46:43

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖