###这里nscm我设置成1 手册意思是说 每一步都让 translation and rotation reset to zero,以免在inpcrd文件中出现*****的溢出。 “The only reason to even reset the coordinates is to prevent
the molecule from diffusing so far away from the origin that its coordinates overflow the format used in restart or trajectory files.“ 但是 这里的CL- 与蛋白离的距离仍然越来越远，这样一直模拟下去的话还是会出现溢出的报错啊，希望大家交流指点下 先谢过了 作者Author: sobereva 时间: 2015-11-14 23:12
怎么有5个蛋白？
加Cl-的目的是什么？如果是作为抗衡离子，隐式溶剂模型下不用加它。
消除平动是对于整体而言的，所以依然可以分离得很远。 作者Author: tbw@zyc 时间: 2015-11-20 20:57
抱歉出差去了，现在才看到回帖
不是5个蛋白是 取的每5ns 最后一帧的图像，能请教"CL-作为抗衡离子，隐式溶剂模型下不用加"的原因吗？ 作者Author: sobereva 时间: 2015-11-21 02:08
只有PBC计算才是要求必须把体系电荷平衡掉的，否则程序无法算。非周期性的隐式溶剂模型计算时没这个必要。你可以先尝试不加Cl模拟一下看看还是否平移。
由于有了隐式溶剂模型，蛋白和Cl-的静电相互作用很弱，而且起始温度又高，很可能吸引不住Cl-，导致自发地越跑越远。虽然设了消平动，但那是对于整个体系而言的。作者Author: tbw@zyc 时间: 2015-11-21 20:04
” Born方程是指在某种介电常数ε的溶剂中，使其中半径为a的粒子带上q电荷所引起自由能的改变为0.5*q^2/ε/a“，这是您以前发过对GB概念描述的贴子，http://sobereva.com/42 “非周期性的隐式溶剂模型计算时没这个必要。" 是说这几个CL-对常数ε没多大影响所以不予考虑吗？
可是不加的话我体系不就成带4库伦的电场体系了，这和实际不符合，不和理，请sob老师详细讲解下为什么没必要？作者Author: sobereva 时间: 2015-11-21 21:21