“…Instead, dispersion-corrected MGGA functionals are required to get the proper balance of forces necessary to accurately model water interactions (Pestana et al, 2017). In previous work (Cockayne and Nelson, 2015;Cockayne, 2017), we were able to reproduce benchmarking energies and geometries of small water clusters quite well by using a combination of a PBEsol (Perdew et al, 2008) gradient term, the RTPSS MGGA density functional (Perdew et al, 2009;Sun et al, 2011), empirical van der Waals forces at the Grimme D2 level (Grimme, 2006), and a Hubbard U value of 7.05 for the water oxygen, but the combination of parameters that we chose was admittedly ad hoc and not derived self-consistently. In this work, we seek a more transferable density functional, one that has been benchmarked across different systems in the literature.…”