计算自由能的第二个lambda时能量最小化时势能无穷大如何解决

yee

大家好，我在计算一个制冷剂分子R125进入到DMAC液体时的自由能时，第一个lambda点0可以计算无错误显示，但是计算lambda第二个点最小化时出现了图示的错误，请问是什么原因。
我主要想算R125溶解到DMAC中的亨利系数，其中亨利系数用到了二者的自由能，请问是通过 lambda=0时无相互作用，lambda=1时vdw-q，计算多个点，然后用gmx bar 或者 MBAR来计算它们的自由能来计算吗？

附件是我的mdp文件。