Gromacs里关于截断半径的一个问题,谢谢!
在读Gromacs手册的时候,有2段话描述了截断半径。对这2段和rlist,rcoulomb,rvdw比较疑惑,想请教各位大虾一下。1. GROMACS employs a pair list that contains those particle pairs for which non-bonded forces must be calculated. The pair list contains atoms i, a displacement vector for atom i, and all particles j that are within rlist of this particular image of atom i. The list is updated every nstlist steps, where nstlist is typically 10.
2. In the neighbor list, all interaction pairs that fall within rlist are stored. Furthermore, the interactions between pairs that do not fall within rlist but do fall within max(rcoulomb,rvdw) are computed during NS. The forces and energy are stored separately and added to short-range forces at every time step between successive NS. If rlist = max(rcoulomb,rvdw), no forces are evaluated during neighbor list generation.
我对上面部分的理解是,两个原子间的距离小于rlist的时候才计算这些原子之间的短程相互作用。相互作用的原子对,如果距离在rlist之内的就被存储下来,没有落到rlist内但是落到max(rcoulomb,rvdw)内的,会在邻近搜索期间被计算。
那么rlist和rcoulomb, rvdw具体是什么关系呢??不知道能不能解释的详细点,谢谢!
还有最后一句话,“If rlist = max(rcoulomb,rvdw), no forces are evaluated during neighbor list generation”,不怎么理解,请大家指教,谢谢!!
页:
[1]