双层PTFE聚合物摩擦模拟的载荷换算与摩擦力计算方法

hansu

想请教一下有相关经验的前辈，在做双层聚合物的摩擦模拟时，宏观上给出的载荷参数要怎么换算成可以施加给原子或者原子组的力？比如有文献里给包含620个原子的上层滑块模型施加了1.78MPa的载荷，文中给出的13.993e-4 kcal/(mol·angstrom)这个力的数值是怎么计算出来的？这个力是直接施加给滑块对应的原子组吗？如果要把滑块设置为刚体，是不是要先通过fix rigid命令设置刚体属性，再用fix addforce命令施加外力？另外，在计算这类模型的摩擦力时，通常是取哪个部分的受力？是固定层在x方向所受的合力（用compute reduce sum计算），还是刚性滑块和牛顿层之间的相互作用力（用compute group/group计算）？目前查阅的文献都没有给出具体的模拟计算细节，导致无法复现结果，还请各位前辈多多指点！

lauthirteen

宏观压强（P, MPa）到微观力（F, LAMMPS单位）的换算，遵循压强的基本定义：

F (力) = P (压强) × A (面积) / N_atoms

其中：

P：目标是宏观压强，例如 1.78 MPa。

A：上层滑块在受力方向上的横截面积