我记得培训时学习,md程序是给定初始坐标和速度,根据分子力场计算能量,受力,更新位置和速度,如此循环。那么比如体系就一个离子,给它个初速度,没有成键也没有非键,理想nve情况应该是匀速直线运动下去吧。我在NAMD2.12 CPU版本里试验,通过设定初始温度的方式是可以正常run的,不过速度似乎并不是稳定的。但是如果不给温度而是使用velocities参数提供pdb速度文件就会吐核退出,没有报错信息。pdb速度文件即使是前面通过初始温度模拟产生的可读vel也不行。而当用binvelocities提供restart.vel文件就可以了。
问题似乎出在了namd读取pdb速度文件上。请教这样问题如何解决?所用文件如下
(, 下载次数 Times of downloads: 2)