Bottom-up Atomistic Descriptions of Top-Down Macroscopic Measurements: Computational Benchmarks for Hammett Electronic Parameters

The ability to relate substituent electronic effects to chemical reactivity is a cornerstone of physical organic chemistry and Linear Free Energy Relationships. The computation of electronic parameters is increasingly attractive since they can be obtained rapidly for structures and substituents without available experimental data and can be applied beyond aromatic substituents. Nevertheless, the description of “top-down” macroscopic observables, such as Hammett parameters using a “bottom-up” computational approach, poses several challenges for the practitioner. We have examined and benchmarked the performance of various computational charge schemes and atomic properties, the locations of the atoms used to obtain these descriptors, and their correlation with empirical Hammett parameters and rate differences resulting from electronic effects. These seemingly small choices have a much more significant impact than previously imagined, which outweighs the level of theory or basis set used. We observe a wide range of performance across the different computational protocols and observe stark and surprising differences in the ability of computational parameters to capture para vs. meta-electronic effects. In general σ(meta) predictions fare much worse than σ(para). As a result, the choice of where to compute these descriptors – for the ring carbons or the attached H or other substituent atoms – affects their ability to capture experimental electronic differences. Hirshfeld charges outperform all other computational descriptors, including several commonly used schemes such as Natural Population Analysis, while using the attached atoms also improves the statistical correlations. We have obtained linear relationships for the global prediction of experimental Hammett parameters from computed descriptors.