Hyperfine Coupling Tensors of the Benzosemiquinone Radical Anion from Car–Parrinello Molecular Dynamics

Gauging the importance of structural parameters for hyperfine coupling constants in organic radicals

The identification of fundamental relationships between atomic configuration and electronic structure typically requires experimental empiricism or systematic theoretical studies. Here, we provide an alternative statistical approach to gauge the importance of structure parameters, i.e., bond lengths, bond angles, and dihedral angles, for hyperfine coupling constants in organic radicals. Hyperfine coupling constants describe electron-nuclear interactions defined by the electronic structure and are experimentally measurable, for example, by electron paramagnetic resonance spectroscopy. Importance quantifiers are computed with the machine learning algorithm neighborhood components analysis using molecular dynamics trajectory snapshots. Atomic-electronic structure relationships are visualized in matrices correlating structure parameters with coupling constants of all magnetic nuclei. Qualitatively, the results reproduce common hyperfine coupling models. Tools to use the presented procedure for other radicals/paramagnetic species or other atomic structure-dependent parameters are provided.

How accurate is density functional theory in predicting spin density? An insight from the prediction of hyperfine coupling constants

Electron paramagnetic resonance (EPR) spectroscopy has been proven to be an important technique for studying paramagnetic systems. Probably, the most accessible EPR parameter and the one that provides a significant amount of information about molecular structure and spin density is the hyperfine coupling constant (HFCC). Hence, accurate quantum-chemical modeling of HFCCs is frequently essential to the adequate interpretation of EPR spectra. It requires the precise spin density, which is the difference between the densities of α- and β-electrons, and thus, its quality is expected to reflect the quality of the total electron density. The question of which approximate exchange-correlation density functional yields sufficiently accurate HFCCs, and thus, the spin density remains open. To assess the performance of well-established density functionals for calculating HFCCs, we used a series of 26 small paramagnetic species and compared the obtained results to the CCSD reference values. The performance of DFT was also tested on EPR-studied o-semiquinone radical interacting with water molecules and Mg²⁺ cation. The HFCCs were additionally calculated by the DLPNO-CCSD method, and this wave function-based technique was found superior to all functionals we tested. Although some functionals were found, on average, to be fairly efficient, we found that the most accurate functional is system-dependent, and therefore, the DLPNO-CCSD method should be preferred for theoretical investigations of the HFCCs and spin density.

Toward an Understanding of the Ambiguous Electron Paramagnetic Resonance Spectra of the Iminoxy Radical from o-Fluorobenzaldehyde Oxime: Density Functional Theory and ab Initio Studies

Iminoxy radicals (R1R2C═N—O•) possess an inherent ability to exist as E and Z isomers. Although isotropic hyperfine couplings for the species with R1 = H allow one to distinguish between E and Z, unequivocal assignment of the parameters observed in the EPR spectra of the radicals without the hydrogen atom at the azomethine carbon to the right isomer is not a simple task. The iminoxyl derived from o-fluoroacetophenone oxime (R1 = CH3 and R2 = o-FC6H5) appears to be a case in point. Moreover, for its two isomers the rotation of the o-FC6H5 group brings into existence the syn and anti conformers, depending on the mutual orientation of the F atom and C═N—O• group, making a description of hyperfine couplings to structure even more challenging. To accomplish this, a vast array of theoretical methods (DFT, OO-SCS-MP2, QCISD) was used to calculate the isotropic hyperfine couplings. The comparison between experimental and theoretical values revealed that the E isomer is the dominant radical form, for which a fast interconversion between anti and syn conformers is expected. In addition, the origin of the significant AF increase with solvent polarity was analyzed.

Theoretical Characterisation of Phosphinyl Radicals and Their Magnetic Properties: g Matrix

The g matrices (g tensors) of various phosphinyl radicals (R2 P(.) ) were calculated using the DFT and multireference configuration interaction (MRCI) methods. The g matrices were distinctly dependent on the molecular structure of the radical. To thoroughly examine this dependence, the contributions from individual atoms and excited states were calculated. The former revealed the gain from the phosphorus atom to be preeminent unless PO or PS bonds are present in the radical molecule. The contributions owing to excited states arising from electronic transitions between doubly occupied molecular orbitals and the SOMO were clearly positive, as in the case of semiquinone and niroxide radicals. The transitions from the phosphorus lone pair were of paramount importance. Surprisingly, unlike for semiquinones and nitroxides, a significant negative contribution was observed from excitations from the SOMO to unoccupied molecular orbitals. For radicals with PO bonds, this contribution to the g2 component was dominant.

First principles calculation of inhomogeneous broadening in solid-state cw-EPR spectroscopy

We present a scheme for the first-principles calculation of EPR lineshapes for continuous-wave-EPR spectroscopy (cw-EPR) of spin centers in complex chemical environments. We specifically focus on poorly characterized systems, e.g. powders and frozen glasses with variable microsolvation structures. Our approach is based on ab initio molecular dynamics simulations and ab initio calculations of the ensemble of g- and A-tensors along the trajectory. The method incorporates temperature effects as well as the full anharmonicity of the intra- and intermolecular degrees of freedom of the system. We apply this scheme to compute the lineshape of a prototypical spin probe, the nitrosodisulfonate dianionic radical (Fremy's salt), dissolved in a 50 : 50 mixture of water and methanol. We are able to determine the specific effect of variations of local solvent composition and microsolvation structure on the cw-EPR lineshape. Our molecular dynamics reveal a highly anisotropic solvation structure with distinct spatial preferences for water and methanol around Fremy's salt that can be traced back to a combination of steric and polar influences. The overall solvation structure and conformational preferences of Fremy's salt as found in our MD simulations agree very well with the results obtained from EPR and orientation-selective ENDOR spectroscopy performed on the frozen glass. The simulated EPR lineshapes show good agreement with the experimental spectra. When combined with our MD results, they characterize the lineshape dependence on local morphological fluctuations.

Predicting the ¹H and ¹³C NMR spectra of paramagnetic Ru(III) complexes by DFT

Nuclear shieldings, including the Fermi-contact and pseudocontact terms, have been calculated with density functional theory (DFT) (nonrelativistic and relativistic) methods in several Ru(III) complexes, thereby predicting (1)H and (13)C paramagnetic shifts. A fair agreement with experimental values is observed. Structural, magnetic and dynamic parameters have also been input to the Solomon-Bloembergen equation in order to predict signal lineshapes. It is shown that DFT-predicted paramagnetic shifts can greatly aid in obtaining and understanding NMR spectra of paramagnetic Ru(III) complexes.