Search

Long-lived organic radicals are promising candidates for the development of high-performance energy solutions such as organic redox batteries, transistors, and light-emitting diodes. However, “stable” organic radicals that remain unreactive for an extended time and that can be stored and handled under ambient conditions are rare. A necessary but not sufficient condition for organic radical stability is the presence of thermodynamic stabilization, such as conjugation with an adjacent π-bond or lone-pair, or hyperconjugation with a σ-bond. However, thermodynamic factors alone do not result in radicals with extended lifetimes: many resonance-stabilized radicals are transient species that exist for less than a millisecond. Kinetic stabilization is also necessary for persistence, such as steric effects that inhibit radical dimerization or reaction with solvent molecules. We describe a quantitative approach to map organic radical stability, using molecular descriptors intended to capture thermodynamic and kinetic considerations. The comparison of an extensive dataset of quantum chemical calculations of organic radicals with experimentally-known stable radical species reveals a region of this feature space where long-lived radicals are located. These descriptors, based upon maximum spin density and buried volume, are combined into a single metric, the radical stability score, that outperforms thermodynamic scales based on bond dissociation enthalpies in identifying remarkably long-lived radicals. This provides an objective and accessible metric for use in future molecular design and optimization campaigns. We demonstrate this approach in identifying Pareto-optimal candidates for stable organic radicals., Molecular descriptors encoding kinetic and thermodynamic stabilization capture the difference between transient and persistent organic radicals.

The effectiveness of the optical rotation prediction (ORP) basis set for computing specific rotations at the coupled cluster (CC) level has been evaluated for a test set of 14 chiral compounds. For this purpose, the ORP basis set has been developed for the second-row atoms present in the investigated systems (that is, for sulfur, phosphorus, and chlorine). The quality of the resulting set was preliminarily evaluated for seven molecules using time-dependent density-functional theory (TD-DFT). Rotations were calculated with the coupled cluster singles and doubles method (CCSD) as well as the second-order approximate coupled cluster singles and doubles method (CC2) with the correlation-consistent aug-cc-pVDZ and aug-cc-pVTZ basis sets and extrapolated to estimate the complete basis-set (CBS) limit for comparison with the ORP basis set. In the compounds examined here, the ORP calculations on molecules containing only first-row atoms compare favorably with results from the larger aug-cc-pVTZ basis set, in some cases lying closer to the estimated CBS limit, while results for molecules containing second-row atoms indicate that larger correlation-consistent basis sets are necessary to obtain reliable estimates of the CBS limit.

The origin of enantioselectivity in asymmetric catalysis is attributed to the energy difference between lower and higher energy diastereomeric transition states, which are respectively responsible for the formation of major and minor enantiomers. Although the increase in the number of transition state models emphasizes the role of weak non-covalent interactions in asymmetric induction, the strength of such interactions is seldom quantified. Through this article, we propose a simple and effective method of quantifying the total non-covalent interaction in stereocontrolling transition states belonging to a group of three representative asymmetric catalytic reactions involving chiral phosphoric acids. Our method relies on rational partitioning of a given transition state into two (or three) sub-units, such that the complex network of intramolecular interactions can be ameliorated to a set of intermolecular interactions between two sub-units. The computed strength of interaction obtained using the counterpoise (CP) method on suitably partitioned transition states provides improved estimates of non-covalent interactions, which are also devoid of basis set superposition error (BSSE). It has been noted that catalysts decorated with larger aromatic arms provide cumulative non-covalent interactions (C-H center dot center dot center dot pi, N-H center dot center dot center dot pi and pi center dot center dot center dot pi) to the tune of 10 to 15 kcal mol(-1). Fine-tuning of the magnitude and nature of these interactions can provide valuable avenues in the design of asymmetric catalysts.