Density Functional Restricted–Unrestricted/Molecular Mechanics Theory for Hyperfine Coupling Constants of Molecules in Solution

Close Relationships between NMR J-Coupling Alternation (JCA) and Molecular Properties of Carbon Chains

We propose a J-coupling alternation (JCA) value that is demonstrated to be a suitable parameter to evaluate the nuclear magnetic resonance (NMR) indirect spin-spin coupling constants (SSCCs) as a function of molecular properties of chains by increasing their length. As an application, we report a theoretical study of the SSCCs for the interactions between neighbor nuclei in increasingly patterned carbon chains within density functional theory. First, we examine the J-coupling constants between ¹H and ¹³C nuclei ( ⁿ J_HC) considering the separation distance, as well as between two adjacent ¹³C nuclei (¹ J_CC) considering their relative positions in polyynes and cumulenes. Further, we define and determine JCA in terms of the differences of ¹ J_CC, which is investigated as a function of several molecular properties, e.g., cohesive energy, characteristic vibrational frequency, average polarizability, and energy gap of the systems. We also determine JCA for other types of carbon chains, such as diphenyl-capped polyynes, polyacetylene and polythiophene. The behavior of JCA as a function of the energy gap may be related to highly π-conjugated low-band-gap carbon chains. Overall, JCA correlates very well with the electronic properties of these chains, especially with their energy gap, exhibiting positive values for pristine polyyne and polythiophene and negative values for pristine cumulene and plyacetylene. These findings indicate an alternative way to establish an appropriate SSCC descriptor that characterizes the electronic nature of the system, such as the proposed JCA value averaging the whole system, instead of using only the individual J-coupling values to give insights into the properties of large and extended systems.

Decomposition of molecular properties

We review recent work on property decomposition techniques using quantum chemical methods and discuss some topical applications in terms of quantum mechanics-molecular mechanics calculations and the constructing of properties of large molecules and clusters. Starting out from the so-called LoProp decomposition scheme [Gagliardi et al., J. Chem. Phys., 2004, 121, 4994] for extracting atomic and inter-atomic contributions to molecular properties we show how this method can be generalized to localized frequency-dependent polarizabilities, to localized hyperpolarizabilities and to localized dispersion coefficients. Some applications of the generalized decomposition technique are reviewed - calculations of frequency-dependent polarizabilities, Rayleigh scattering of large clusters, and calculations of hyperpolarizabilities of proteins.

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.

Internal Dynamics of the 3-Pyrroline-N-Oxide Ring in Spin-Labeled Proteins

Site-directed spin labeling is a versatile tool to study structure as well as dynamics of proteins using EPR spectroscopy. Methanethiosulfonate (MTS) spin labels tethered through a disulfide linkage to an engineered cysteine residue were used in a large number of studies to extract structural as well as dynamic information on the protein from the rotational dynamics of the nitroxide moiety. The ring itself was always considered to be a rigid body. In this contribution, we present a combination of high-resolution X-ray crystallography and EPR spectroscopy of spin-labeled protein single crystals demonstrating that the nitroxide ring inverts fast at ambient temperature while exhibiting nonplanar conformations at low temperature. We have used quantum chemical calculations to explore the potential energy that determines the ring dynamics as well as the impact of the geometry on the magnetic parameters probed by EPR spectroscopy.

Convergence of QM/MM and Cluster Models for the Spectroscopic Properties of the Oxygen-Evolving Complex in Photosystem II

The latest crystal structure of photosystem II at 1.9 Å resolution, which resolves the topology of the Mn4CaO5 oxygen evolving complex (OEC) at atomistic detail, enables a better correlation between structural features and spectroscopic properties than ever before. Building on the refined crystallographic model of the OEC and the protein, we present combined quantum mechanical/molecular mechanical (QM/MM) studies of the spectroscopic properties of the natural catalyst embedded in the protein matrix. Focusing on the S2 state of the catalytic cycle, we examine the convergence of not only structural parameters but also of the intracluster magnetic interactions in terms of exchange coupling constants and of experimentally relevant (55)Mn, (17)O, and (14)N hyperfine coupling constants with respect to QM/MM partitioning using five QM regions of increasing size. This enables us to assess the performance of the method and to probe second sphere effects by identifying amino acid residues that principally affect the spectroscopic properties of the OEC. Comparison between QM-only and QM/MM treatments reveals that whereas QM/MM models converge quickly to stable values, the QM cluster models need to incorporate significantly larger parts of the second coordination sphere and surrounding water molecules to achieve convergence for certain properties. This is mainly due to the sensitivity of the QM-only models to fluctuations in the hydrogen bonding network and ligand acidity. Additionally, a hydrogen bond that is typically omitted in QM-only treatments is shown to determine the hyperfine coupling tensor of the unique Mn(III) ion by regulating the rotation plane of the ligated D1-His332 imidazole ring, the only N-donor ligand of the OEC.

A gauge invariant multiscale approach to magnetic spectroscopies in condensed phase: general three-layer model, computational implementation and pilot applications

Analytical equations to calculate second order electric and magnetic properties of a molecular system embedded into a polarizable environment are presented. The treatment is limited to molecules described at the self consistent field level of theory, including Hartree-Fock theory as well as Kohn-Sham density functional theory and is extended to the Gauge-Including Atomic Orbital method. The polarizable embedding is described by means of our already implemented polarizable quantum mechanical/molecular mechanical (MM) methodology, where the polarization in the MM layer is handled by means of the fluctuating charge (FQ) model. A further layer of description, i.e, the polarizable continuum model, can also be included. The FQ(/polarizable continuum model) contributions to the properties are derived, with reference to the calculation of the magnetic susceptibility, the nuclear magnetic resonance shielding tensor, electron spin resonance g-tensors, and hyperfine couplings.

π-Stacking effects on the EPR parameters of a prototypical DNA spin label

The character and value of spin labels for probing environments like double-stranded DNA depend on the degree of change in the spin Hamiltonian parameters of the spin label induced by the environment. Herein we provide a systematic theoretical investigation of this issue, based on a density functional theory method applied to a spin labeled DNA model system, focusing on the dependence of the EPR properties of the spin label on the π stacking and hydrogen bonding that occur upon incorporating the spin label into the selected base pair inside DNA. It is found that the EPR spin Hamiltonian parameters of the spin label are only negligibly affected by its incorporation into DNA, when compared to its free form. This result gives a theoretical ground for the common empirical assumption regarding the behaviour of spin Hamiltonian parameters made in EPR based measurements of the distance between spin labels incorporated into DNA.

Improving the calculation of Electron Paramagnetic Resonance hyperfine coupling tensors for d-block metals

Calculation of hyperfine coupling constants (HFCs) of Electron Paramagnetic Resonance from first principles can be a beneficial complement to experimental data in cases where the molecular structure is unknown. We have recently investigated basis set convergence of HFCs in d-block complexes and obtained a set of basis functions for the elements Sc-Zn, which were saturated with respect to both the Fermi contact and spin-dipolar components of the hyperfine coupling tensor [Hedegård et al., J. Chem. Theory Comput., 2011, 7, 4077-4087]. Furthermore, a contraction scheme was proposed leading to very accurate, yet efficient basis sets for the elements Sc-Zn. Here this scheme is tested against a larger test set of molecules and a wider range of DFT functionals. We further investigate the regular aug-cc-pVTZ and core-valence correlation aug-cc-pCVTZ basis sets as well as another core-property basis set, CP(PPP). While aug-cc-pVTZ-J provides hyperfine coupling constants that are almost identical to the converged series (aug-cc-pVTZ-Juc), we observe that not only the regular but also the core-valence correlation basis sets provide results far from the converged results. The usage of specialized core-basis sets leads to a large and highly significant improvement of the calculated hyperfine couplings in comparison with experimental data.