Computational modeling of polyoxotungstates by relativistic DFT calculations of (183)W NMR chemical shifts

The (183)W nuclear shielding in a variety of tungsten polyoxometalates (POM) (Lindqvist, Anderson, decatungstates, Keggin) of different shapes and charges has been modeled by DFT calculations that take into account relativistic effects, by means of the zero-order regular approximation (ZORA), and solvent effects, by the conductor-like screening model (COSMO) continuum method. The charge/surface area ratio (q/A) is proposed as an indicator of the charge density to which the solvation energies of all POMs are correlated in a satisfactory way. Among the various theoretical levels tested (ZORA scalar or spin-orbit, frozen-core or all-electron basis set, geometry optimization in the gas phase or in the continuum solvent, etc.), the best results are obtained when both geometry optimization in solvent and spin-orbit shielding are included (mean absolute error of delta=35 ppm). The quality of the computed chemical shifts depends systematically on the charge density as expressed by q/A; thus, POMs with low q/A ratios display the best agreement with experimental data. The performance of the method is such that computed values can aid the assignment of the (183)W NMR spectra of polyoxotungstates, as shown by the case of alpha-[PW(11)TiO(40)](5-), whose six signals are ranked computationally so as to almost reproduce the experimental ordering even though the signals are spaced by as little as 5 ppm.     

DFT Calculations of 31P NMR Chemical Shifts in Palladium Complexes

In this study, comparative analysis of calculated (GIAO method, DFT level) and experimental ³¹P NMR shifts for a wide range of model palladium complexes showed that, on the whole, the theory reproduces the experimental data well. The exceptions are the complexes with the P=O phosphorus, for which there is a systematic underestimation of shielding, the value of which depends on the flexibility of the basis sets, especially at the geometry optimization stage. The use of triple-ζ quality basis sets and additional polarization functions at this stage reduces the underestimation of shielding for such phosphorus atoms. To summarize, in practice, for the rapid assessment of ³¹P NMR shifts, with the exception of the P=O type, a simple PBE0/{6-311G(2d,2p); Pd(SDD)}//PBE0/{6-31+G(d); Pd(SDD)} approximation is quite acceptable (RMSE = 8.9 ppm). Optimal, from the point of view of “price–quality” ratio, is the PBE0/{6-311G(2d,2p); Pd(SDD)}//PBE0/{6-311+G(2d); Pd(SDD)} (RMSE = 8.0 ppm) and the PBE0/{def2-TZVP; Pd(SDD)}//PBE0/{6-311+G(2d); Pd(SDD)} (RMSE = 6.9 ppm) approaches. In all cases, a linear scaling procedure is necessary to minimize systematic errors.     

Multinuclear Diffusion NMR Spectroscopy and DFT Modeling: A Powerful Combination for Unraveling the Mechanism of Phosphoester Bond Hydrolysis Catalyzed by Metal‐Substituted Polyoxometalates

A detailed reaction mechanism is proposed for the hydrolysis of the phosphoester bonds in the DNA model substrate bis(4‐nitrophenyl) phosphate (BNPP) in the presence of the Zr^IV‐substituted Keggin type polyoxometalate (Et₂NH₂)₈[{α‐PW₁₁O₃₉Zr(μ‐OH) (H₂O)}₂] ? 7 H₂O (ZrK 2:2) at pD 6.4. Low‐temperature ³¹P DOSY spectra at pD 6.4 gave the first experimental evidence for the presence of ZrK 1:1 in fast equilibrium with ZrK 2:2 in purely aqueous solution. Moreover, theoretical calculations identified the ZrK 1:1 form as the potentially active species in solution. The reaction intermediates involved in the hydrolysis were identified by means of ¹H/³¹P NMR studies, including EXSY and DOSY NMR spectroscopy, which were supported by DFT calculations. This experimental/theoretical approach enabled the determination of the structures of four intermediate species in which the starting compound BNPP, nitrophenyl phosphate (NPP), or the end product phosphate (P) is coordinated to ZrK 1:1. In the proposed reaction mechanism, BNPP initially coordinates to ZrK 1:1 in a monodentate fashion, which results in hydrolysis of the first phosphoester bond in BNPP and formation of NPP. EXSY NMR studies showed that the bidentate complex between NPP and ZrK 1:1 is in equilibrium with monobound and free NPP. Subsequently, hydrolysis of NPP results in P, which is in equilibrium with its monobound form.     

Simulation of 59Co NMR Chemical Shifts in Aqueous Solution

⁵⁹Co chemical shifts were computed at the GIAO‐B3LYP level for [Co(CN)₆]^3?, [Co(H₂O)₆]³⁺, [Co(NH₃)₆]³⁺, and [Co(CO)₄]^? in water. The aqueous solutions were modeled by Car–Parrinello molecular dynamics (CPMD) simulations, or by propagation on a hybrid quantum‐mechanical/molecular‐mechanical Born–Oppenheimer surface (QM/MM‐BOMD). Mean absolute deviations from experiment obtained with these methods are on the order of 400 and 600 ppm, respectively, over a total δ(⁵⁹Co) range of about 18 000 ppm. The effect of the solvent on δ(⁵⁹Co) is mostly indirect, resulting primarily from substantial metal–ligand bond contractions on going from the gas phase to the bulk. The simulated solvent effects on geometries and δ(⁵⁹Co) values are well reproduced by using a polarizable continuum model (PCM), based on optimization and perturbational evaluation of quantum‐mechanical zero‐point corrections.     

Towards Accurate Predictions of Proton NMR Spectroscopic Parameters in Molecular Solids

The factors contributing to the accuracy of quantum-chemical calculations for the prediction of proton NMR chemical shifts in molecular solids are systematically investigated. Proton chemical shifts of six solid amino acids with hydrogen atoms in various bonding environments (CH, CH₂, CH₃, OH, SH and NH₃) were determined experimentally using ultra-fast magic-angle spinning and proton-detected 2D NMR experiments. The standard DFT method commonly used for the calculations of NMR parameters of solids is shown to provide chemical shifts that deviate from experiment by up to 1.5 ppm. The effects of the computational level (hybrid DFT functional, coupled-cluster calculation, inclusion of relativistic spin-orbit coupling) are thoroughly discussed. The effect of molecular dynamics and nuclear quantum effects are investigated using path-integral molecular dynamics (PIMD) simulations. It is demonstrated that the accuracy of the calculated proton chemical shifts is significantly better when these effects are included in the calculations.     

The Relationship between NMR Chemical Shifts of Thermally Polarized and Hyperpolarized 89Y Complexes and Their Solution Structures

Recently developed dynamic nuclear polarization (DNP) technology offers the potential of increasing the NMR sensitivity of even rare nuclei for biological imaging applications. Hyperpolarized ⁸⁹Y is an ideal candidate because of its narrow NMR linewidth, favorable spin quantum number (I= ), and long longitudinal relaxation times (T₁). Strong NMR signals were detected in hyperpolarized ⁸⁹Y samples of a variety of yttrium complexes. A dataset of ⁸⁹Y NMR data composed of 23 complexes with polyaminocarboxylate ligands was obtained using hyperpolarized ⁸⁹Y measurements or ¹H,⁸⁹Y‐HMQC spectroscopy. These data were used to derive an empirical equation that describes the correlation between the ⁸⁹Y chemical shift and the chemical structure of the complexes. This empirical correlation serves as a guide for the design of ⁸⁹Y sensors. Relativistic (DKH2) DFT calculations were found to predict the experimental ⁸⁹Y chemical shifts to a rather good accuracy.     

Atomic Contributions from Spin‐Orbit Coupling to 29Si NMR Chemical Shifts in Metallasilatrane Complexes

New members of a novel class of metallasilatrane complexes [X‐Si‐(μ‐mt)₄‐M‐Y], with M=Ni, Pd, Pt, X=F, Cl, Y=Cl, Br, I, and mt=2‐mercapto‐1‐methylimidazolide, have been synthesized and characterized structurally by X‐ray diffraction and by ²⁹Si solid‐state NMR. Spin‐orbit (SO) effects on the ²⁹Si chemical shifts induced by the metal, by the sulfur atoms in the ligand, and by heavy halide ligands Y=Cl, Br, I were investigated with the help of relativistic density functional calculations. Operators used in the calculations were constructed such that SO coupling can selectively be switched off for certain atoms. The unexpectedly large SO effects on the ²⁹Si shielding in the Ni complex with X=Y=Cl reported recently originate directly from the Ni atom, not from other moderately heavy atoms in the complex. With respect to Pd, SO effects are amplified for Ni owing to its smaller ligand‐field splitting, despite the smaller nuclear charge. In the X=Cl, Y=Cl, Br, I series of complexes the Y ligand strongly modulates the ²⁹Si shift by amplifying or suppressing the metal SO effects. The pronounced delocalization of the partially covalent M←Y bond plays an important role in modulating the ²⁹Si shielding. We also demonstrate an influence from the X ligand on the ²⁹Si SO shielding contributions originating at Y. The NMR spectra for [X‐Si‐(μ‐mt)₄‐M‐Y] must be interpreted mainly based on electronic and relativistic effects, rather than structural differences between the complexes. The results highlight the sometimes unintuitive role of SO coupling in NMR spectra of complexes containing heavy atoms.     

Computational Prediction of 1H and 13C NMR Chemical Shifts for Protonated Alkylpyrroles: Electron Correlation and Not Solvation is the Salvation

Prediction of chemical shifts in organic cations is known to be a challenge. In this article we meet this challenge for α-protonated alkylpyrroles, a class of compounds not yet studied in this context, and present a combined experimental and theoretical study of the ¹³C and ¹H chemical shifts in three selected pyrroles. We have investigated the importance of the solvation model, basis set, and quantum chemical method with the goal of developing a simple computational protocol, which allows prediction of ¹³C and ¹H chemical shifts with sufficient accuracy for identifying such compounds in mixtures. We find that density functional theory with the B3LYP functional is not sufficient for reproducing all ¹³C chemical shifts, whereas already the simplest correlated wave function model, Møller–Plesset perturbation theory (MP2), leads to almost perfect agreement with the experimental data. Treatment of solvent effects generally improves the agreement with experiment to some extent and can in most cases be accomplished by a simple polarizable continuum model. The only exception is the NH proton, which requires inclusion of explicit solvent molecules in the calculation.     

Paramagnetic NMR Shifts for Saddle-Shaped Five-Coordinate Iron(III) Porphyrin Complexes with Intermediate-Spin Structure

Calculable results: Complex density functional calculations and spin distribution analyses have been performed for planar and saddled iron(III) porphyrin complexes. The spin populations and the extent of the interactions between the metal and the porphyrin orbitals were determined, which can explain the large change of meso-carbon atom chemical shifts observed for different porphyrin ligands.     

An Investigation of 1:1 Adducts of Gallium Trihalides with Triarylphosphines by Solid‐State 69/71Ga and 31P NMR Spectroscopy

Several 1:1 adducts of gallium trihalides with triarylphosphines, X₃Ga(PR₃) (X=Cl, Br, and I; PR₃=triarylphosphine ligand), were investigated by using solid‐state ^69/71Ga and ³¹P NMR spectroscopy at different magnetic‐field strengths. The ^69/71Ga nuclear quadrupolar coupling parameters, as well as the gallium and phosphorus magnetic shielding tensors, were determined. The magnitude of the ⁷¹Ga quadrupolar coupling constants (C_Q(⁷¹Ga)) range from approximately 0.9 to 11.0 MHz . The spans of the gallium magnetic shielding tensors for these complexes, δ₁₁?δ₃₃, range from approximately 30 to 380 ppm; those determined for phosphorus range from 10 to 40 ppm. For any given phosphine ligand, the gallium nuclei are most shielded for X=I and least shielded for X=Cl, a trend previously observed for In^III–phosphine complexes. This experimental trend, attributed to spin‐orbit effects of the halogen ligands, is reproduced by DFT calculations. The signs of C_Q(^69/71Ga) for some of the adducts were determined from the analysis of the ³¹P NMR spectra acquired with magic angle spinning (MAS). The ¹J(^69/71Ga,³¹P) and ΔJ(^69/71Ga, ³¹P) values, as well as their signs, were also determined; values of ¹J(⁷¹Ga,³¹P) range from approximately 380 to 1590 Hz. Values of ¹J(^69/71Ga,³¹P) and ΔJ(^69/71Ga,³¹P) calculated by using DFT have comparable magnitudes and generally reproduce experimental trends. Both the Fermi‐contact and spin‐dipolar Fermi‐contact mechanisms make important contributions to the ¹J(^69/71Ga,³¹P) tensors. The ³¹P NMR spectra of several adducts in solution, obtained as a function of temperature, are contrasted with those obtained in the solid state. Finally, to complement the analysis of NMR spectra for these adducts, single‐crystal X‐ray diffraction data for Br₃Ga[P(p‐Anis)₃] and I₃Ga[P(p‐Anis)₃] were obtained.     

The influence of solvent molecules on NMR spectrum of barbituric acid in the DMSO solution

This work shows the modification of barbituric acid (BA) chemical shifts by dimethylsulphoxide (DMSO) molecules. The discussed changes are caused by creation of the H-bonded associates formed by barbituric acid with DMSO in solution. Free molecule of barbituric acid, the cluster of BA with two DMSO molecules and two different clusters of BA with four DMSO units are taken into consideration. The chemical shifts of these systems have been calculated and the obtained results have been compared with experimental data. Theoretical calculations predict a significant downfield shift for imino protons of barbituric acid involved in intermolecular-N-H...DMSO hydrogen bonds. The influence of the solvent molecules on other nuclei chemical shifts, especially protons of barbituric acid methylene group, is also reported. The calculations have involved Hartree-Fock and several Density Functional Theory methods. All methods correctly describe experimental ¹H and ¹³C NMR spectra of barbituric acid. The best consistence between experiment and theory is observed for the BLYP functional. Four approximations of magnetic properties calculations embedded in the Gaussian’98 package have been tested. The results of the performed calculations indicate that from a practical point of view the GIAO method should be preferred.     

Crystallographic and Dynamic Aspects of Solid‐State NMR Calibration Compounds: Towards ab Initio NMR Crystallography

The excellent results of dispersion‐corrected density functional theory (DFT‐D) calculations for static systems have been well established over the past decade. The introduction of dynamics into DFT‐D calculations is a target, especially for the field of molecular NMR crystallography. Four ¹³C ss‐NMR calibration compounds are investigated by single‐crystal X‐ray diffraction, molecular dynamics and DFT‐D calculations. The crystal structure of 3‐methylglutaric acid is reported. The rotator phases of adamantane and hexamethylbenzene at room temperature are successfully reproduced in the molecular dynamics simulations. The calculated ¹³C chemical shifts of these compounds are in excellent agreement with experiment, with a root‐mean‐square deviation of 2.0 ppm. It is confirmed that a combination of classical molecular dynamics and DFT‐D chemical shift calculation improves the accuracy of calculated chemical shifts.     

Prediction of the 13C NMR Chemical Shifts of Fluorenone Analogues by the GIAO Method

After the geometry optimization at B3LYP/6-31+G(d,p) level,the NMR calcula-tions of a series of fluorenone analogues have been carried out by GIAO method at HF/6-31+G(d) level and B3LYP/6-311G+(2d,p) level,respectively.The 13C NMR chemical shifts calculated at HF/6-31+G(d) level show better agreement with the observed values.By a series of linear correction equations (δpred=a + bδcalc),accurate prediction of 13C chemical shifts was achieved for the new fluorenone compound.The linear correlation of δpred with δexptl is excellent,and the square of correlation coefficient,r2,is up to 0.994.The maximum absolute difference between δpred and δexptl,Δδ,is 4.6 ppm,and the root-mean-square error between δpred and δexptl is only 2.6 ppm.