To carry out the study described above, we require molecular dynamics trajectories of water box simulated at 21 temperatures (T ) × 21 pressures (P ) × water-like force fields (F ). We therefore ran 21×21×5 = 2205 independent simulations. Each run consists of energy minimization, equilibration at a constant temperature, equilibration at constant pressure, and final production run to generate a trajectory for a given set of (T, P, F ). We used optimized molecular dynamics simulation parameters, nonbonded interaction cut-offs, switching distances, and so forth, as found in our previous study and given in Section 3 in the SI. As explained above, trajectory generation is done only once because of the high phase overlap, and optimization is done with reweighting alone, with verification using additional simulations afterwards.