For this study, we parameterize this simple molecular model to optimally fit the thermodynamic properties of liquid water over a temperature range of 273.15-373.15 K and pressure range of 1-5000 atm. The pressure range in particular is significantly beyond what has been explored in previous fixed-charge classical water models. One could, of course, make the choice to find models that span the entire phase equilibrium of water, as was attempted by Wang et al. However, these models use only a few T and P points, rather than the entire jointly varying surface in temperature and pressure. Within the context of demonstrating this parameterization approach, we note it is difficult to sample near phase boundaries because of density fluctuations and sampling is especially hard in solid state where the kinetics is slow. Additionally, the large phase space separation of liquid and solid phase spaces prohibits reweighting techniques between these states and the subsequent calculation of observables using molecular dynamics simulations of one phase to predict the other phase. We leave extensions to such approaches for further work, as one can argue that models for the liquid phase alone are still useful.