Calculations Of Electron Paramagnetic Resonance Parameters Research Articles

A Joint Venture of Ab Initio Molecular Dynamics, Coupled Cluster Electronic Structure Methods, and Liquid-State Theory to Compute Accurate Isotropic Hyperfine Constants of Nitroxide Probes in Water.

The isotropic hyperfine coupling constant (HFCC, Aiso) of a pH-sensitive spin probe in a solution, HMI (2,2,3,4,5,5-hexamethylimidazolidin-1-oxyl, C9H19N2O) in water, is computed using an ensemble of state-of-the-art computational techniques and is gauged against X-band continuous wave electron paramagnetic resonance (EPR) measurement spectra at room temperature. Fundamentally, the investigation aims to delineate the cutting edge of current first-principles-based calculations of EPR parameters in aqueous solutions based on using rigorous statistical mechanics combined with correlated electronic structure techniques. In particular, the impact of solvation is described by exploiting fully atomistic, RISM integral equation, and implicit solvation approaches as offered by ab initio molecular dynamics (AIMD) of the periodic bulk solution (using the spin-polarized revPBE0-D3 hybrid functional), embedded cluster reference interaction site model integral equation theory (EC-RISM), and polarizable continuum embedding (using CPCM) of microsolvated complexes, respectively. HFCCs are obtained from efficient coupled cluster calculations (using open-shell DLPNO-CCSD theory) as well as from hybrid density functional theory (using revPBE0-D3). Re-solvation of “vertically desolvated” spin probe configuration snapshots by EC-RISM embedding is shown to provide significantly improved results compared to CPCM since only the former captures the inherent structural heterogeneity of the solvent close to the spin probe. The average values of the Aiso parameter obtained based on configurational statistics using explicit water within AIMD and from EC-RISM solvation are found to be satisfactorily close. Using either such explicit or RISM solvation in conjunction with DLPNO-CCSD calculations of the HFCCs provides an average Aiso parameter for HMI in aqueous solution at 300 K and 1 bar that is in good agreement with the experimentally determined one. The developed computational strategy is general in the sense that it can be readily applied to other spin probes of similar molecular complexity, to aqueous solutions beyond ambient conditions, as well as to other solvents in the longer run.

Controlling One-Electron vs Two-Electron Pathways in the Multi-Electron Redox Cycle of Nickel Diethyldithiocarbamate.

The unique redox cycle of NiII(dtc)2, where dtc- is N,N-diethyldithiocarbamate, in acetonitrile displays 2e- redox chemistry upon oxidation from NiII(dtc)2 → [NiIV(dtc)3]+ but 1e- redox chemistry upon reduction from [NiIV(dtc)3]+ → NiIII(dtc)3 → NiII(dtc)2. The underlying reasons for this cycle lie in the structural changes that occur between four-coordinate NiII(dtc)2 and six-coordinate [NiIV(dtc)3]+. Cyclic voltammetry (CV) experiments show that these 1e- and 2e- pathways can be controlled by the addition of pyridine-based ligands (L) to the electrolyte solution. Specifically, the addition of these ligands resulted in a 1e- ligand-coupled electron transfer (LCET) redox wave, which produced a mixture of pyridine-bound Ni(III) complexes, [NiIII(dtc)2(L)]+, and [NiIII(dtc)2(L)2]+. Although the complexes could not be isolated, electron paramagnetic resonance (EPR) measurements using a chemical oxidant in the presence of 4-methoxypyridine confirmed the formation of trans-[NiIII(dtc)2(L)2]+. Density functional theory calculations were also used to support the formation of pyridine coordinated Ni(III) complexes through structural optimization and calculation of EPR parameters. The reversibility of the LCET process was found to be dependent on both the basicity of the pyridine ligand and the scan rate of the CV experiment. For strongly basic pyridines (e.g., 4-methoxypyridine) and/or fast scan rates, high reversibility was achieved, allowing [NiIII(dtc)2(L)x]+ to be reduced directly back to NiII(dtc)2 + xL. For weakly basic pyridines (e.g., 3-bromopyridine) and/or slow scan rates, [NiIII(dtc)2(L)x]+ decayed irreversibly to form [NiIV(dtc)3]+. Detailed kinetics studies using CV reveal that [NiIII(dtc)2(L)]+ and [NiIII(dtc)2(L)2]+ decay by parallel pathways due to a small equilibrium between the two species. The rate constants for ligand dissociation ([NiIII(dtc)2(L)2]+ → [NiIII(dtc)2(L)]+ + L) along with decomposition of [NiIII(dtc)2(L)]+ and [NiIII(dtc)2(L)2]+ species were found to increase with the electron-withdrawing character of the pyridine ligand, indicating pyridine dissociation is likely the rate-limiting step for decomposition of these complexes. These studies establish a general trend for kinetically trapping 1e- intermediates along a 2e- oxidation path.