and WM21-D3(BJ), we optimize a double hybrid density functional of the form ωB97X-MP2-D3(BJ). For this we allow the coefficients of opposite-spin and same-spin MP2 to be free parameters. In addition we let the s 6 coefficient in D3(BJ) be free to account for the long-range behavior of MP2. As with ωB97X21-D3(BJ), we use the same non-linear parameters as WM21-D3(BJ) and use ωB97X orbitals to calculate the energies. Following the procedure used previously and taking the combination which minimizes the weighted validation error with the minimum number of parameters without skipping any B97 terms.