The values of were extracted using gmx energy, and the standard errors were obtained from block averaging performed by gmx analyze. The 0 values were taken from the simulations of the pure air-water interface with the corresponding water model. The error of ⇧ was estimated as the sum of the standard errors of the corresponding 0 and values.