Rather than constraining this expansion to fit a fixed function, the coefficients could be fitted directly to data, and this forms an integral part of the development described in this paper. In this way the opposite and same spin coefficients could be optimized separately, however, in this work we continue to use the opposite-spin only UW12 (osUW12) approach used in B-LYP-osUW12-D3(BJ).