Correct me if I am wrong, please. Ugas would be obtained from a simulation of a single molecule in vacuum. And U liq would be the total potential energy of the system (which is given by g_energy) minus the total bonded potential energy (bonds+angles+dihedrals). Am I thinking right? Thanks. -Pallavi Banerjee