| 查看: 595 | 回复: 0 | ||
[资源]
Amber MD模拟中的 restart file 坐标值过大的问题
|
|
今天在用Amber跑MD的时候出现了以下的错误 forrtl: severe (64): input conversion error, unit 9, file /misc/home/mabiao/amberMD/casestudy/ham/md_0921/02/testequil/heat.rst.007 Image PC Routine Line Source sander.MPI 0000000000CE8F8F Unknown Unknown Unknown sander.MPI 0000000000CE746E Unknown Unknown Unknown sander.MPI 0000000000CBABD4 Unknown Unknown Unknown sander.MPI 0000000000C7A5A6 Unknown Unknown Unknown sander.MPI 0000000000C7A207 Unknown Unknown Unknown sander.MPI 0000000000C951AF Unknown Unknown Unknown sander.MPI 00000000004F184C Unknown Unknown Unknown google了之后看到有别的朋友也遇到过类似问题，具体请参考如下连接：http://archive.ambermd.org/200907/0008.html 有网友解释说是因为跑长时间MD的时候有时会出现rst文件中的原子坐标值超出预定值而出现＊＊＊＊＊字符，造成后续计算在读取这个rst文件的时候出错。如果是做水环境周期边界计算的话，在input文件中加入iwrap ＝ 1就能防止出现原子坐标超大的问题。 我的问题不适合上面的做法，因为我没有加水并且不是周期边界条件计算，加了iwrap＝1，反而出错：Error: IWRAP=1 cannot be used without a periodic box. 最后仔细分析出错的rst文件，发现坐标变大，出现＊＊＊＊符号的原子是我系统的倒数第八个原子，最后八个原子的坐标也比别的大。看了一下初始结构的PDB文件发现最后八个原子是我在leap中为了中和带电量而加的八个Cl原子。 因为做的是真空环境模拟不用考虑电中性，而我画蛇添足多加了几个原子后计算过程中他们的坐标被放大很多，最后出现rst文件出错。 因为是新手，所以很多时候照着教程死板硬套就容易出现意向不到的错误，在这里借小木虫论坛把我的学习经验写出来，以给后来的朋友做下参考！ 有不妥之处希望高手们能给指出！ |
回复此楼» 猜你喜欢
294求调剂
已经有8人回复
-
324求调剂
已经有14人回复
-
285求调剂
已经有6人回复
-
一志愿上海大学生物学346
已经有3人回复
-
322求调剂：一志愿湖南大学 材料与化工（085600），已过六级。
已经有16人回复
-
调剂
已经有10人回复
-
283求调剂
已经有7人回复
-
0703调剂
已经有6人回复
-
288环境专硕,求调材料方向
已经有6人回复
-
求调剂
已经有5人回复
找到一些相关的精华帖子，希望有用哦~
北航版MATLABr2011a教程(张志涌)电子版[这个是真的2011a]之[doc版]
已经有451人回复
- 后处理中display时如何指定坐标面？
已经有5人回复
- opls下L-J参数转换成amber下L-J参数的问题
已经有9人回复
- 生测数据的标准差大于平均值，怎么办？
已经有7人回复
- Amber tools 中的antechamber 必须要在安装了amber11的前体下才能使用么?
已经有4人回复
- 求Gromacs中Glycam力场数据包。。非常感谢
已经有5人回复
- DPD 关于粒子数密度的取值
已经有3人回复
- 有关MS 的建模 优化及计算的一些问题
已经有6人回复
- amber中硅原子键参数如何获取？
已经有8人回复
- AMBER中qm/mm的问题
已经有5人回复
- [求助]AMBER数据分析问题
已经有9人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
