The performance of a large variety of contemporary density functional theory (DFT), double-hybrid DFT, and high-level Gaussian-n (Gn) procedures has been evaluated for the calculation of bond dissociation energies (BDEs) and radical stabilization energies (RSEs) associated with NAX bonds (X ¼ H, Cl). The chosen set of 62 NAX systems (31 NAH and 31 NACl) span a wide range of biologically relevant species. As reference values, we used benchmark-quality W2w data that we recently obtained as part of a systematic thermochemical study of substituent effects in these systems. Of the Gn schemes, the modified G4 procedures (G4-5H and G4(MP2)-6X) perform somewhat better than the corresponding standard G4 procedures for the NAX BDEs of these systems. For the NAH RSEs, G3X, G3X(MP2), G3X(MP2)-RAD, G4-5H, and G4(MP2)-6X emerge as excellent performers, with mean absolute deviations (MADs) from the benchmark W2w values of 0.9-1.4 kJ mol -1 . However, for the NACl RSEs, G4 is the best performer, with an MAD of 1.7 kJ mol -1 . The BDEs of both NAH and NACl bonds represent a challenge for DFT procedures. In particular, only a handful of functionals (namely, B3P86, M05-2X, M06-2X, and ROB2-PLYP) perform well, with MADs 4.5 kJ mol À1 for both bond types. Nearly all of the considered DFT procedures perform signifi-