Search papers, labs, and topics across Lattice.
This paper introduces a refined force-matching protocol for developing accurate and reproducible water models for molecular simulation. The key improvements are the use of PBE0 hybrid-functional reference forces, incorporating exact-exchange effects, and a separable nonlinear least-squares formulation that analytically eliminates linear force-field amplitudes. Applying this protocol to CP2K reference forces, the resulting water model accurately reproduces target force distributions and radial distribution functions, demonstrating the effectiveness of the approach.
A new force-matching protocol slashes the complexity of water model development, enabling more accurate and reproducible simulations by analytically eliminating linear parameters and using hybrid-functional forces.
Force-matched water models provide a practical route from first-principles reference data to long classical and path-integral molecular simulations. Previous flexible four-site potentials in the spirit of q-TIP4P/F showed that fitting analytic models to density-functional-theory forces can reproduce key structural features of liquid water while retaining the efficiency required for quantum-nuclear sampling. Here we introduce two refinements aimed at making this strategy more accurate and more reproducible for molecular simulation. First, the production fit is based on PBE0 hybrid-functional reference forces and therefore includes the Hartree--Fock exact-exchange contribution in the electronic-structure target. Second, the parametrization is formulated as a separable nonlinear least-squares problem in which all linear force-field amplitudes are eliminated analytically for every trial set of nonlinear shape parameters. The resulting force-matching protocol lowers the dimension of the nonlinear search, reduces compensation between heterogeneous parameters, and enables a controlled comparison of Lennard-Jones and Buckingham oxygen--oxygen repulsion terms. Applied to CP2K reference forces for liquid water, the projected optimization reproduces target force distributions and yields radial distribution functions close to first-principles and neutron-scattering benchmarks. The Buckingham representation gives a more flexible short-range repulsive wall than a purely Lennard-Jones form, and the final flexible model remains stable in path-integral simulations. The method provides a transparent workflow for deriving simulation-ready water potentials from accurate hybrid density-functional reference forces.
