使用链接单元算法,因此,计算成本与原子数呈线性关系。
可以选择激活周期性边界条件。
应该定义一个三次区域来生成邻居列表。
原子坐标只能是 3D 格式,作为行向量。
% [nebrTab, nebrLen] = md_neighlist_cube(x, nx, rc, srgdat, bctyp, frc)% % ------------------------------------------------- --------------------- % 为立方体 % 区域中的原子系统创建邻居列表% WZ。 单%------------------------------------------------- --------------- % 输入: % x : 原子坐标,[ px, py, pz] % nx : 数量。 原子数,= size( x, 1) % rc : 截止
2021-10-12 13:41:09
3KB
matlab
1