需要把能量温度曲线做出来之后用麦克斯韦等面积法则算吧,毕竟直接模拟会有热滞之类的现象,不管热滞的话做长时间的模拟退火?用Maxwell等面积法则算的话大概要用广义正则系综(GCE),我这有篇文献你可以看看,但是GCE怎么在GROMACS里实现我还不知道,只知道LAMMPS里有,我在想办法实现。
文献:Shao, M. Z., Wang, Y. T., & Zhou, X. (2020). Fast and accurate determination of phase transition temperature via individual generalized canonical ensemble simulation. Chinese Physics B, 29(8), 2–7. 作者Author: sobereva 时间: 2021-4-8 01:09
不同温度下分别做模拟,看什么时候不再是液体(通过扩散系数可判断)