分子动力学模拟中计算最耗时的就是非键相互作用,假设你就N个原子 那么每个2之间就需要计算一次相互作用力,那么就是CN2=N*(N-1)/2 如果你是计算多体势,那么需要计算的数量会更多,所以我们需要考虑处理 最近邻像约定,尤其是针对于PBC体系 这个是不是最近邻像约定 只是个人突然想到盒子不能太小 同时也用这个图告诉大家如果是 长程作用,如库仑作用 就不能用此方法进行去做约束 原则上A会受到C的影响,但是C只是镜像 对应的是B,如果用B作为计算那么距离会边长 如果周期的C都计算,那么也是错的 这个时候我们只计算最近这个周期的位置 |
Powered by Discuz! X3.5
© 2001-2024 Discuz! Team.