Gromacs非平衡模拟求助

sun666

在文献看到用gromacs做的非平衡模拟，对离子通过孔道的过程施加加速度，从而计算扩散系数转化为菲克定律中的D。
想请教各位如何在gromacs对离子施加加速度进行非平衡模拟，在施加加速度后又是如何计算扩散系数与平衡模拟有何不同呢。
因为之前没有见过做非平衡模拟的例子，还请各位推荐这方面的学习资料以及例子。