As with WM21-D3(BJ) we re-optimize the coefficients after recalculating the energies in the training set self-consistently and repeating for multiple cycles until converged. Table shows the coefficients for the hybrid functional at each stage of the procedure. Four of the Table : Optimized coefficients for the ωB97X21-D3(BJ) functional at each stage of the optimization. The first set of coefficients are fitted using energies calculated using orbitals optimised with the ωB97X functional. Subsequent fittings are performed on the training set data obtained using the functional defined in the previous column. The final functional parameters are shown in bold. x all changing by more than 0.1, though these have converged after the second iteration. We note that the amount of same spin correlation contribution is greatly reduced by this re-optimization. Also we note that the non-linear parameters for both range-separation and dispersion are not re-optimized with the same ones used for both functionals.