求助：计算制冷剂和离子液体的粘度与文献误差较大如何解决？

yee

我计算了300个R1234yf分子在278K，5atm下的粘度，mdp文件根据老师网站水的粘度教程 使用GROMACS计算粘度 来用的，但老师的mdp全部都没有continuation，而且只在第一个npt的时候生成了速度，我才用这种方法时，模拟跑不了，gromacs提示段错误，后面我把mdp文件改成第一个nvt就生成速度，后面全部continuation=yes后就能跑了，单出来的结果与文献差25%左右，energy出来的粘度err.est也很大。另外模拟的bmim+tf2n（300+300）的粘度误差差了3个数量级。
但这两者的模拟密度和文献是吻合的。我也尝试了改小cos_acceleration，但误差反而变大了。
请问是力场的问题还是mdp设置的问题，还是模拟时间不够。