“…11 This has led to the advent of several many-body PEFs for water in which the individual many-body terms in eqn (1) are tted to large sets of reference data calculated at the coupled cluster level of theory with single, double, and iterative triple excitations, i.e., CCSD(T), the current "gold standard" for molecular interactions. [12][13][14][15][16][17] Among existing many-body PEFs, it has been shown that the MB-pol PEF achieves great accuracy by integrating classical representations of many-body electrostatics, which correctly describe interactions at large molecular separations, with data-driven 2B and 3B terms represented by permutationally invariant polynomials (PIPs) that smoothly turn on as two and three water molecules, respectively, approach each other. [15][16][17] The MB-pol 2B and 3B terms, which were tted to large sets of CCSD(T) reference data calculated in the complete basis set (CBS) limit, effectively account for non-classical interactions (e.g., Pauli repulsion, charge transfer and charge penetration) arising in regions where the electron densities of individual water molecules overlap.…”