At those geometries, single point energy evaluations are performed by the explicitly correlated coupled cluster method, including single, double, and (perturbatively) triple excitations (CCSD(T)-F12) within the frozen-core approximation and in conjunction with the j3Z basis set. Next, CBS extrapolation, core-valence correlation (CV), and possibly other minor terms are added at different levels depending on the specific model chemistry. Finally, the experimental values of spin-orbit couplings are employed for O, OH, SH, and Cl radicals, lowering their electronic energies by 0.9, 0.8, 2.3, and 3.5 kJ mol -1 , respectively. All the DFT computations have been performed with the Gaussian code, 37 F12 calculations with the Molpro package and CCSDT or CCSDT(Q) energy evaluations with the MRCC program. Finally, diagonal Born-Oppenheimer Corrections (DBOC) and relativistic contributions have been computed by the CFOUR code. The junChS-F12 model chemistry