| 查看: 655 | 回复: 5 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】用VASP做MD，能将每步的速度信息都给出吗？
|
|||
| 请教大家了，我现在由于要看一个MD进行到不同程度，然后冷却的效果，所以想做一个较长的MD，然后取他们在不同程度的构型和速度接着做冷却，但是VASP，我找说明找了一遍，包括NWRITE的书写，和VASPRUN.XML,和XDATACAR的说明都没有找到怎么把每步的速度信息给出，谢谢大家了。 |
回复此楼» 猜你喜欢
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有7人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有8人回复
-
申请2026年博士
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有7人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
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