The same procedure was previously used to optimize the TIP3P-FB three-point model. The TIP4P-ST four-point model used the TIP4P-Ew functional form and four parameters were optimized: the virtual-site position that carries the negative charge, the hydrogen charge q H , and the Lennard-Jones parameters for oxygen σ LJ and LJ . Starting values of the parameters are given in Table .