Magnetic Properties of Nitroxide Spin Probes: Reliable Account of Molecular Motions and Nonspecific Solvent Effects by Time-Dependent and Time-Independent Approaches

Integration of Quantum Chemistry, Statistical Mechanics, and Artificial Intelligence for Computational Spectroscopy: The UV-Vis Spectrum of TEMPO Radical in Different Solvents

The ongoing integration of quantum chemistry, statistical mechanics, and artificial intelligence is paving the route toward more effective and accurate strategies for the investigation of the spectroscopic properties of medium-to-large size chromophores in condensed phases. In this context we are developing a novel workflow aimed at improving the generality, reliability, and ease of use of the available computational tools. In this paper we report our latest developments with specific reference to unsupervised atomistic simulations employing non periodic boundary conditions (NPBC) followed by clustering of the trajectories employing optimized feature spaces. Next accurate variational computations are performed for a representative point of each cluster, whereas intracluster fluctuations are taken into account by a cheap yet reliable perturbative approach. A number of methodological improvements have been introduced including, e.g., more realistic reaction field effects at the outer boundary of the simulation sphere, automatic definition of the feature space by continuous perception of solute–solvent interactions, full account of polarization and charge transfer in the first solvation shell, and inclusion of vibronic contributions. After its validation, this new approach has been applied to the challenging case of solvatochromic effects on the UV–vis spectra of a prototypical nitroxide radical (TEMPO) in different solvents. The reliability, effectiveness, and robustness of the new platform is demonstrated by the remarkable agreement with experiment of the results obtained through an unsupervised approach characterized by a strongly reduced computational cost as compared to that of conventional quantum mechanics and molecular mechanics models without any accuracy reduction.

A radical approach to radicals

Nitroxide radicals are characterized by a long-lived spin-unpaired electronic ground state and are strongly sensitive to their chemical surroundings. Combined with electron paramagnetic resonance spectroscopy, these electronic features have led to the widespread application of nitroxide derivatives as spin labels for use in studying protein structure and dynamics. Site-directed spin labelling requires the incorporation of nitroxides into the protein structure, leading to a new protein-ligand molecular model. However, in protein crystallographic refinement nitroxides are highly unusual molecules with an atypical chemical composition. Because macromolecular crystallography is almost entirely agnostic to chemical radicals, their structural information is generally less accurate or even erroneous. In this work, proteins that contain an example of a radical compound (Chemical Component Dictionary ID MTN) from the nitroxide family were re-refined by defining its ideal structural parameters based on quantum-chemical calculations. The refinement results show that this procedure improves the MTN ligand geometries, while at the same time retaining higher agreement with experimental data.

Towards reliable references for electron paramagnetic resonance parameters based on quantum chemistry: the case of verdazyl radicals

We present an efficient and accurate computational procedure to calculate properties measurable by EPR spectroscopy. We simulate a molecular dynamics (MD) trajectory by employing the quantum mechanically derived force field (QMDFF) [S. Grimme, J. Chem. Theory Comput., 2014, 10, 4497] and sample the trajectory at different time steps. For each snapshot EPR properties are calculated with a hybrid density functional theory (DFT) method. EPR spectra are simulated based on the averaged results. We applied the strategy to a number of previously published and novel verdazyl radicals, for which we recorded EPR spectra. The resulting simulated spectra are compatible with experiment already before employing an additional fitting step, in contrast to those from single point electronic-structure calculations. After the refinement, the experimental data are excellently reproduced, and the fitted EPR parameters do not deviate much from the calculated ones. This provides confidence in ascribing a direct physical meaning to the refined data in terms of experimental EPR parameters rather than merely considering them as mathematical fit parameters. We also find that couplings to hydrogen nuclei have a significant influence on the spectra of verdazyl radicals.

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.

Multifrequency Pulsed EPR and the Characterization of Molecular Dynamics

In fluid solution, motion-dependent processes dominate electron spin-lattice relaxation for nitroxides and semiquinones at frequencies between 250 MHz and 34 GHz. For triarylmethyl radicals, motion-dependent processes dominate spin-lattice relaxation at frequencies below about 3 GHz. The frequency dependence of relaxation provides invaluable information about dynamic processes occurring with characteristic times on the order of the electron Zeeman frequency. Relaxation mechanisms, methods of measuring spin-lattice relaxation, and motional processes for nitroxide, semiquinone, and triarylmethyl radicals are discussed.

Modeling EPR parameters of nitrogen containing conjugated radical cations

DFT investigation on conjugated radical cations containing¹⁴N nucleus to obtain accurate isotropic hyperfine coupling constants.

Electron spin-lattice relaxation mechanisms of rapidly-tumbling nitroxide radicals

Electron spin relaxation times at 295 K were measured at frequencies between 250 MHz and 34 GHz for perdeuterated 2,2,6,6-tetramethyl-4-piperidone-1-oxyl (PDT) in five solvents with viscosities that result in tumbling correlation times, τR, between 4 and 50 ps and for three (14)N/(15)N pairs of nitroxides in water with τR between 9 and 19 ps. To test the impact of structure on relaxation three additional nitroxides with τR between 10 and 26 ps were studied. In this fast tumbling regime T2(-1)~T1(-1) at frequencies up to about 9 GHz. At 34 GHz T2(-1)>T1(-1) due to increased contributions to T2(-1) from incomplete motional averaging of g-anisotropy, and T2(-1)-T1(-1) is proportional to τR. The contribution to T1(-1) from spin rotation is independent of frequency and decreases as τR increases. Spin rotation dominates T1(-1) at 34 GHz for all τR studied, and at all frequencies studied for τR=4 ps. The contribution to T1(-1) from modulation of nitrogen hyperfine anisotropy increases as frequency decreases and as τR increases; it dominates at low frequencies for τR>~15 ps. The contribution from modulation of g anisotropy is significant only at 34 GHz. Inclusion of a thermally-activated process was required to account for the observation that for most of the radicals, T1(-1) was smaller at 250 MHz than at 1-2 GHz. The significant (15)N/(14)N isotope effect, the small H/D isotope effect, and the viscosity dependence of the magnitude of the contribution from the thermally-activated process suggest that it arises from intramolecular motions of the nitroxide ring that modulate the isotropic A values.

Extension of the AMBER Force Field for Nitroxide Radicals and Combined QM/MM/PCM Approach to the Accurate Determination of EPR Parameters of DMPO-H in Solution

A computational strategy that combines both time-dependent and time-independent approaches is exploited to accurately model molecular dynamics and solvent effects on the isotropic hyperfine coupling constants of the DMPO-H nitroxide. Our recent general force field for nitroxides derived from AMBER ff99SB is further extended to systems involving hydrogen atoms in β-positions with respect to NO-moiety. The resulting force-field has been employed in a series of classical molecular dynamics simulations, comparing the computed EPR parameters from selected molecular configurations to the corresponding experimental data in different solvents. The effect of vibrational averaging on the spectroscopic parameters is also taken into account, by second-order vibrational perturbation theory involving semidiagonal third energy derivatives together first and second property derivatives.

"Pump-probe" atom-centered density matrix propagation studies to gauge anharmonicity and energy repartitioning in atmospheric reactive adducts: case study of the OH + isoprene and OH + butadiene reaction intermediates

Time-resolved "pump-probe" ab initio molecular dynamics studies are constructed to probe the stability of reaction intermediates, the mechanism of energy transfer, and energy repartitioning, for moieties involved during the interaction of volatile organic compunds with hydroxyl radical. These systems are of prime importance in the atmosphere. Specifically, the stability of reaction intermediates of hydroxyl radical adducts to isoprene and butadiene molecules is used as a case study to develop novel computational techniques involving "pump-probe" ab initio molecular dynamics. Starting with the various possible hydroxyl radical adducts to isoprene and butadiene, select vibrational modes of each of the adducts are populated with excess energy to mimic the initial conditions of an experiment. The flow of energy into the remaining modes is then probed by subjecting the excited adducts to ab initio molecular dynamics simulations. It is found that the stability of the adducts arises directly due to the anhormonically driven coupling of the modes to facilitate repartitioning of the excess vibrational energy. This kind of vibrational repartitioning has a critical influence on the energy density.

Zero-point vibrational corrections to isotropic hyperfine coupling constants in polyatomic molecules

The present work addresses isotropic hyperfine coupling constants in polyatomic systems with a particular emphasis on a largely neglected, but a posteriori significant, effect, namely zero-point vibrational corrections. Using the density functional restricted-unrestricted approach, the zero-point vibrational corrections are evaluated for the allyl radical and four of its derivatives. In addition for establishing the numerical size of the zero-point vibrational corrections to the isotropic hyperfine coupling constants, we present simple guidelines useful for identifying hydrogens for which such corrections are significant. Based on our findings, we critically re-examine the computational procedures used for the determination of hyperfine coupling constants in general as well as the practice of using experimental hyperfine coupling constants as reference data when benchmarking and optimizing exchange-correlation functionals and basis sets for such calculations.