使用LAMMPS内建应力计算功能输出的应力张量结果与文献不同

Heinrich

我使用LAMMPS模拟LJ流体Argon的气液两相系统，并用以下代码输出了沿z方向上各bin的应力张量。

compute cc1 all chunk/atom bin/1d z lower 1.0 units box

compute stress_atom all stress/atom NULL

fix fix_stress all ave/chunk 1 10000000 10000000 cc1 &
    density/number &
    c_stress_atom[1] c_stress_atom[2] c_stress_atom[3] &
    file profile.stress

模拟的结果如附件图片所示。可以看到在我的结果中，法向应力分布在界面处存在一正一负两个峰值（从左数第2张图中蓝线），而在以往的文献中的分子动力学模拟结果中（右侧第1张图），法向应力在界面附近恒定为0。请问我的模拟和计算方法存在什么问题？LAMMPS INPUT文件也上传到附件中了。