如何在分子模拟中实现 “reflecting particle membrane”效果？

含光君

读到一篇文献介绍了一种 “reflecting particle membrane” （RPM）方法，其主要功能是在分子模拟中设置一个虚拟的膜，沿指定方向到达该膜平面的分子可以自由通过，逆着指定方向到达该膜的分子有概率会被反弹。示意图如下：



文献链接：https://journals.aps.org/pre/pdf/10.1103/PhysRevE.57.7259

该方法可用于在纳米流体中模拟一个压力梯度（Nano Lett. 2008, 8, 9, 2788–2793）：


如果要在GMX或LAMMPS中实现上述RPM效果，请问大家有没有好的思路或者替代的方法？