Indirect Nuclear Spin Spin Coupling Tensors Research Articles

SCFit: Software for single-crystal NMR analysis. Free vs constrained fitting

The design and implementation of a software package for the analysis of single-crystal NMR data is presented. The SCFit software can treat spectra arising from various interactions: (i) chemical shift tensor only; (ii) chemical shift tensor and quadrupolar coupling tensor; (iii) dipolar and indirect nuclear spin-spin coupling tensors; (iv) all four interactions. The software is demonstrated on recently reported 17O and 31P single-crystal NMR data for triphenylphosphine oxide and for two of its halogen-bonded cocrystals. The 17O single-crystal NMR data represent a case where all four above-mentioned interactions simultaneously affect the spectra. SCFit can fit the chemical shift and quadrupolar coupling in two ways: (i) through an unconstrained fitting process where all tensor parameters are freely optimized or (ii) through a constrained fitting process where the principal components of the tensors may be fixed to values known previously with high precision via the analysis of powder samples. The second strategy is explored in an effort to reduce the number of unknowns in the fitting process; an improvement in the precision of the resulting tensor orientations is noted in some cases.

Single-Crystal NMR Characterization of Halogen Bonds.

Oxygen-17-enriched triphenylphosphine oxide and three of its halogen-bonded cocrystals featuring 1,4-diiodotetrafluorobenzene and 1,3,5-trifluoro-2,4,6-triiodobenzene as halogen bond donors have been characterized by 31P and 17O single-crystal NMR spectroscopy. Single-crystal NMR allows for the measurement of not only the magnitudes of various NMR interaction tensors, but also their orientations relative to the crystal lattice and therefore relative to the halogen bonds themselves. 31P chemical shift tensors, 17O chemical shift tensors, 17O quadrupolar coupling tensors, and 31P-17O indirect nuclear spin-spin (J) coupling tensors are reported here for P═O···I halogen bonds. The angular deviations in the directions of the pseudo-unique components of the 31P chemical shift tensors, the 17O chemical shift tensors, and the 17O quadrupolar coupling tensors from the direction of the oxygen-iodine halogen bond correlate with the deviations in linearity of the P═O···I halogen bond. There is also a clear decrease in anisotropy and an increase in asymmetry of the J(31P,17O) coupling tensors attributable to the formation of iodine-oxygen halogen bonds. The small but quantifiable changes in the tensors are consistent with the weak nature of these halogen bonds relative to the P═O motif. Overall, this work establishes single-crystal NMR as a novel probe of halogen bonds in solids. Analysis of the results has provided insights into the correlations between the magnitude and orientation of various NMR interaction tensors and the local geometry of the halogen bond. Gauge-including projector-augmented wave computations corroborate the experimental findings.

Calculation of nuclear spin-spin coupling constants using frozen density embedding

We present a method for a subsystem-based calculation of indirect nuclear spin-spin coupling tensors within the framework of current-spin-density-functional theory. Our approach is based on the frozen-density embedding scheme within density-functional theory and extends a previously reported subsystem-based approach for the calculation of nuclear magnetic resonance shielding tensors to magnetic fields which couple not only to orbital but also spin degrees of freedom. This leads to a formulation in which the electron density, the induced paramagnetic current, and the induced spin-magnetization density are calculated separately for the individual subsystems. This is particularly useful for the inclusion of environmental effects in the calculation of nuclear spin-spin coupling constants. Neglecting the induced paramagnetic current and spin-magnetization density in the environment due to the magnetic moments of the coupled nuclei leads to a very efficient method in which the computationally expensive response calculation has to be performed only for the subsystem of interest. We show that this approach leads to very good results for the calculation of solvent-induced shifts of nuclear spin-spin coupling constants in hydrogen-bonded systems. Also for systems with stronger interactions, frozen-density embedding performs remarkably well, given the approximate nature of currently available functionals for the non-additive kinetic energy. As an example we show results for methylmercury halides which exhibit an exceptionally large shift of the one-bond coupling constants between (199)Hg and (13)C upon coordination of dimethylsulfoxide solvent molecules.

A ZORA-DFT and NLMO study of the one-bond fluorine–X indirect nuclear spin-spin coupling tensors for various VSEPR geometries

Zeroth-order regular approximation (ZORA) density functional theory (DFT) calculations of one-bond X–19F indirect nuclear spin-spin coupling (J) tensors were performed on a series of fluorine-containing compounds covering several valence shell electron pair repulsion (VSEPR) theory geometries for which J, by symmetry, is not required to be axially symmetric. The calculations show that the antisymmetric components of J are only of the same order of magnitude as the principal components of the symmetric J-coupling tensor for a few geometries, and that in cases of approximate axial symmetry along the bond, J remains nearly axially symmetric with its unique component along the bond. In general, different species having the same nominal geometry tend to have similar tensor orientations, magnitudes of anisotropy of J relative to the isotropic coupling constant, as well as the same dominant contributions from the different coupling mechanisms. Structures are also systematically modified to determine how the tensor components depend on geometrical parameters. The isotropic coupling constants are subsequently interpreted using a natural localized molecular orbital (NLMO) approach. Our results could prove to be useful for future experimental characterizations of J tensors in systems having symmetry properties that do not force J to be axially symmetric or coincident with the dipolar coupling tensor.

Evaluation of the influence of anisotropic indirect nuclear spin-spin coupling tensors on effective residual dipolar couplings for model peptides.

Residual dipolar couplings (RDCs) observed between nuclear spins in molecules in partially oriented media have become a valuable source of information for NMR spectroscopists seeking to structurally characterize biological macromolecules. Examination of the form of the direct (D) and indirect (J) nuclear spin-spin coupling Hamiltonians indicates that all observed RDCs contain an unknown contribution from the anisotropic part of J (Delta J) in addition to the direct dipolar contribution, DPQ. Here, we evaluate the influence of DeltaJ on RDCs through a series of DFT calculations on model peptides. Very small corrections to one-bond RDCs measured between heavy atoms in peptides and proteins are recommended: +0.51% for N-C' spin pairs, and +0.45% for C alpha-C' spin pairs. The corrections to RDCs involving at least one proton are negligible. This latter point is likely to be equally applicable to nucleic acids and oligosaccharides in addition to peptides and proteins. Finally, the orientations of the J(N, C') and J(C alpha, C') tensors in the molecular framework are reported for glycylglycine.

Interpretation of indirect nuclear spin-spin coupling tensors for polyatomic xenon fluorides and group 17 fluorides: results from relativistic density-functional calculations.

Significant improvements have been made recently in the calculation of NMR indirect nuclear spin-spin coupling tensors (J). In particular, the relativistic zeroth-order regular approximation density-functional theory (ZORA-DFT) approach holds great promise for the calculation of spin-spin coupling constants for a variety of chemical systems containing heavy nuclei. In the present work, the ZORA-DFT method is applied to the calculation of the complete reduced coupling tensors, K, for a range of chlorine-, bromine-, iodine-, and xenon-containing species: K(Cl,F) for ClF(2)(+), ClF(3), ClF(4)(+), ClF(5), ClF(6)(-), and ClF(6)(+); K(Br,F) for BrF(3), BrF(6)(-), and BrF(6)(+); K(I,F) for IF(4)(+) and IF(6)(+); K(Xe,F) for XeF(+), XeF(2), XeF(3)(+), XeF(4), XeF(5)(-), XeF(5)(+), and XeF(7)(+). These species represent a wide variety of geometrical bonding arrangements. Agreement between the calculated coupling constants and available experimental data is excellent, and the absolute sign of the coupling constants is provided. It is shown that (1)K(iso) may be positive or negative even within the same molecule, e.g., K(Cl,F)(iso) may be of either sign, depending on the local environment. Periodic trends in (1)K(iso) for isovalent and isostructural molecules are evident. The spin-spin coupling anisotropies, Delta K, and the orientations of the K tensors are also determined. The success of the calculations is a direct result of employing reliable geometries and considering both scalar and spin-orbit relativistic effects. The dependence of K(Cl,F)(iso) and K(Xe,F)(iso) on the local molecular and electronic structure is discussed in terms of the paramagnetic spin-orbit (PSO) and combined Fermi-contact spin-dipolar (FC+SD) coupling mechanisms. The PSO term depends strongly on the number of valence shell electron lone pairs on the central heavy atom, and the FC+SD contribution increases with the Cl[bond]F or Xe[bond]F bond length for a given series of compounds. This interpretation allows for the successful rationalization of the existing experimental data.

Periodic trends in indirect nuclear spin-spin coupling tensors: relativistic density functional calculations for interhalogen diatomics.

There have been significant advances in the calculation and interpretation of indirect nuclear spin-spin coupling (J) tensors during the past few years; however, much work remains to be done, especially for molecules containing heavy atoms where relativistic effects may play an important role. Many J tensors cannot be explained based solely on a nonrelativistic Fermi-contact mechanism. In the present work, the relativistic zeroth-order regular approximation density-functional (ZORA-DFT) implementation for the calculation of J has been applied to the complete series of homonuclear and heteronuclear diatomic halogen molecules: F(2), Cl(2), Br(2), I(2), At(2), ClF, BrF, IF, ClBr, ClI, and BrI. For all of these compounds, the reduced isotropic coupling constant (K(iso)) is positive and the reduced anisotropic coupling constant (DeltaK) is negative. With the exception of molecular fluorine, the magnitudes of K(iso) and DeltaK are shown to increase linearly with the product of the atomic numbers of the coupled nuclei. ZORA-DFT calculations of J for F(2) and ClF are in excellent agreement with the results obtained from multiconfigurational self-consistent-field calculations. The relative importance of the various coupling mechanisms is approximately constant for all of the compounds, with the paramagnetic spin-orbit term being the dominant contributor to K(iso), at approximately 70-80%. Available experimental stimulated resonant Raman spectroscopy data are exploited to extract the complete J((127)I,(127)I) tensor for iodine in two rotational states. The dependence of K(iso) and DeltaK on bond length and rovibrational state is investigated by using calculated results in combination with available experimental data. In addition to providing new insights into periodic trends for J coupling tensors, this work further demonstrates the utility of the ZORA-DFT method and emphasizes the necessity of spin-orbit relativistic corrections for J calculations involving heavy nuclei.

Second-order polarization propagator calculations of indirect nuclear spin-spin coupling tensors in the water molecule

The contribution to indirect nuclear spin-spin coupling tensors provided by the Fermi contact, the spin-dipolar, the Fermi contact/spin-dipolar crossterm, and the paramagnetic spin-orbit interactions are investigated in a zeroth-, first- (the same as the coupled Hartree-Fock method), and second-order polarization propagator approach. Numerical applications to the water molecule show that the second-order results for both the HO and the HH coupling constants are in good agreement with experimental data - especially if vibrational corrections and the diamagnetic spin-orbit contributions are taken into account. We find that the correlation corrections beyond coupled Hartree-Fock are important. We also report how the second-order results are influenced by neglect of some of the most time-consuming steps in the calculation.