Reliable Quantum Chemical Prediction of the Localized/Delocalized Character of Organic Mixed-Valence Radical Anions. From Continuum Solvent Models to Direct-COSMO-RS

A recently proposed quantum-chemical protocol for the description of the character of organic mixed-valence (MV) compounds, close from both sides to the localized/delocalized borderline, is evaluated and extended for a series of dinitroaryl radical anions 1-6. A combination of global hybrid functionals with exact-exchange admixtures of 35% (BLYP35) or 42% (BMK) with appropriate solvent modeling allows an essentially quantitative treatment of, for example, structural symmetry breaking in Robin/Day class II systems, thermal electron transfer (ET) barriers, and intervalence charge transfer (IV-CT) excitation energies, while covering also the delocalized class III cases. Global hybrid functionals with lower exact exchange admixtures (e.g., B3LYP, M05, or M06) provide a too delocalized description, while functionals with higher exact exchange admixtures (M05-2X, M06-2X) provide a too localized one The B2PLYP double hybrid gives reasonable structures but far too small barriers in class II cases. The CAM-B3LYP range hybrid gives somewhat too high ET barriers and IV-CT energies, while the range hybrids omega B97X and LC-BLYP clearly exhibit too much exact exchange. Continuum solvent models describe the situation well in most aprotic solvents studied. The transition of 1,4-dinitrobenzene anion 1 from a class HI behavior in aprotic solvents to a class II behavior in alcohols is not recovered by continuum solvent models. In contrast, it is treated faithfully by the novel direct conductor-like screening model for real solvents (D-COSMO-RS). The D-COSMO-RS approach, the TURBOMOLE implementation of which is reported, also describes accurately the increased ET barriers of class II systems 2 and 3 in alcohols as compared to aprotic solvents and can distinguish at least qualitatively between different aprotic solvents with identical or similar dielectric constants. The dominant role of the solvent environment for the ET character of these MV radical anions is emphasized, as in contrast to some previous computational suggestions essentially all of the present systems have delocalized class III character in the gas phase. The present approach allows accurate estimates from the gas phase to aprotic and protic solvent environments, without the need for explicit ab initio molecular dynamics simulations, and without artificial constraints.