Gas phase standard state (298.15 K, 1 atm) isomerization energies (∆ isom E) were calculated using the M062X functional with the QZVP, 6-311++G(d,p), 6-311++G(2d,2p), and cc-pVTZ basis sets against the 24 reactions in the ISOL set of benchmark isomerization energies for large organic molecules. The M062X functional appears to offer comparable isomerization energy prediction performance to the best performing currently available dispersion corrected functionals against this benchmark dataset.The performance of various levels of quantum chemical methods for estimating the isomerization energies (∆ isom E) and enthalpies (∆ isom H) of organic compounds has been the focus of a number of investigations over the past several years [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15]. This attention is primarily driven by the poor performance of many popular density functionals -most notably the B3LYP functional -with regard to these types of intramolecular rearrangements. Several new types of functionals and corrections have been developed in response to these reported theoretical shortcomings. However, prior ∆ isom H/∆ isom E benchmarking efforts have concentrated on relatively small molecules (<8 heavy atoms).Recently, Huenerbein et al. [15] have proposed the ISOL set of 24 isomerization reactions for large organic compounds (Figure 1). In their benchmarking study, the authors obtained reference ∆ isom E at the SCS-MP3/CBS level on B97D/TZVP optimized structures. Based on previous efforts showing the M062X density functional performs well for computing ∆ isom H/∆ isom E of smaller organic molecules [6,[10][11][12], in the current work we assess the performance of this functional for calculating ∆ isom E for the ISOL benchmark set.Single point calculations were performed using Gaussian 09 [16] and the M062X density functional [17] with the QZVP [18], 6-311++G(d,p) and 6-311++G(2d,2p) [19][20][21][22][23], and cc-pVTZ [24,25] basis sets on the B97D/TZVP optimized geometries from the ISOL set in Huenerbein et al. [15], and corresponding ∆ isom E were calculated along with associated method specific error metrics against the reference SCS-MP3/CBS ∆ isom E (Table 1) (QZVP,p), and 6-311++G(2d,2p)) or two (cc-pVTZ) outliers present with maximum absolute deviations ranging from 13.8 to 17.8 kcal mol −1 depending on the basis set, also similar to that of the PBE0 (MAXD=9.1 kcal mol −1 ; 0 outliers), B2GP-PLYP-D (MAXD=11.5 kcal mol −1 ; 1 outlier), BMK-D (MAXD=17.1 kcal mol −1 ; 1 outlier), and PW6B95-D (MAXD=11.9 kcal mol −1 ;