请教分子动力学中固液界面张力计算的两个问题

ymsw_zzw

大家好，我想用分子动力学模拟来计算水在固体壁面上的界面张力，于是我用lammps导出了体系中每一个原子在x、y、z方向上的受力。在写后处理脚本的过程中，我有几个问题不太确定，想请教一下各位大佬。

1.假如固液界面是垂直于z方向的，且固液界面的面积为S。我把体系沿z方向划分成若干薄层，那么界面张力计算公式中的Pzz是否就是这一薄层中所有原子在z方向的受力之和，再除以S？Pxx和Pyy分别是这一薄层中所有原子在x方向、y方向的受力之和，再除以S？

2.计算界面张力时，按照这一公式，是否就是先计算出某一薄层的Pzz-0.5(Pxx+Pyy)，在乘以这一薄层的厚度，得到一个值比如说a。再对所有薄层的a值求和，就得到了固液界面张力？



我是一个小白，还请各位赐教！

sharklmh

这个式子的出处你有吗