使用gromacs计算液滴体系表面张力的相关问题

秋心

各位老师好，我最近使用gromacs在NVT系综下模拟了一个含有水分子、无机盐离子、有机物分子三元组成的中性液滴体系（模拟之后这个液滴并不是规则的球形），我想计算这个液滴体系的表面张力，目前我调研发现有两种方法可以实现计算表面张力，但是我不太确定那种方法对于我这个体系来说是比较严谨的。
第一种：借助压力张量的概念，使用Irving-Kirkwood方法计算得到相应的压力向量法向分量（gromacs是否也可以得到类似的压力向量法向分量？），然后带入公式计算表面张力，具体的细节请见球形纳米液滴表面张力的计算 - 知乎 (zhihu.com)。
第二种：使用gromacs的energy命令，选择SurfTen，直接得到表面张力的数值，但是sob老师在一个帖子中提到这样得到的表面张力是表面积和表面张力相乘，我不太理解，也不确定这个表面张力是否严谨；还有就是使用gromacs的energy命令提取压力，Pxx,Pyy,Pzz，然后使用gromacs手册给出的表面张力计算公式去计算相应的表面张力，但是这个方法好像只适用于界面平行于XY平面的体系。
希望各位老师前辈能指点一下迷津，万分感谢！！

wch53

请问我这种只有上层有表面活性剂的界面计算问题？就是我这种情况是否适合使用呢？我感觉两处 油水界面的界面张力是不一样的

Jus

前辈您好，请问您解决了吗，我也想计算球形界面的界面张力