| 查看: 697 | 回复: 5 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】用VASP做MD，能将每步的速度信息都给出吗？
|
|||
| 请教大家了，我现在由于要看一个MD进行到不同程度，然后冷却的效果，所以想做一个较长的MD，然后取他们在不同程度的构型和速度接着做冷却，但是VASP，我找说明找了一遍，包括NWRITE的书写，和VASPRUN.XML,和XDATACAR的说明都没有找到怎么把每步的速度信息给出，谢谢大家了。 |
回复此楼» 猜你喜欢
求调剂
已经有23人回复
-
301求调剂
已经有15人回复
-
304求调剂（085602，过四级，一志愿985）
已经有17人回复
-
302分求调剂 一志愿安徽大学085601
已经有12人回复
-
288环境专硕,求调材料方向
已经有23人回复
-
环境专硕调剂
已经有6人回复
-
22408 调剂材料
已经有6人回复
-
285求调剂
已经有12人回复
-
求调剂
已经有6人回复
-
085600材料与化工301分求调剂院校
已经有19人回复
2楼2009-12-21 13:05:09
3楼2009-12-21 13:20:28
★ ★ ★ ★ ★ ★
ice_rain(金币+1,VIP+0):谢谢交流 12-21 17:50
limaojlu(金币+5,VIP+0): 12-21 18:19
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,每一次的速度都会写出. 方法不唯一, 不过都类似. |
4楼2009-12-21 14:46:43
5楼2009-12-21 18:19:40
6楼2009-12-21 18:22:13
