求助：双层PTFE模型摩擦模拟计算

清风与鹿

请问各位前辈做双层聚合物摩擦模拟时，宏观给的载荷参数如何换算成施加给原子或者group的力？

例如下面文献中给约620个原子的上层模型（滑块）施加1.78MPa的载荷，文中给的13.993e-4 kcal/(mol·angstrom)是如何计算的？这个力是施加给滑块group的吗？如果把滑块看成一个刚体，是先设置fix rigid，然后再fix addforce吗？



另外，在计算如下模型的摩擦力时，一般是哪部分的受力？是固定层所受的x方向力的总和？（compute reduce sum）还是刚性滑块和牛顿层之间的力？（compute group/group）


目前看到的文献并没有具体模拟计算说明，所以结果无法复制，请各位前辈不吝赐教！