where L z is the elongation of the z-axis, and P xx , P yy , and P zz are the diagonal components of the pressure tensor. To perform the analysis, we used the gmx energy routine of the GROMACS suite. Experimental reference data for the surface tension are obtained from Ref. Lessons learned from multi-objective automatic optimizations of classical three-site rigid water models using microscopic and macroscopic target experimental observables