Solvent effects in the GIAO‐DFT calculations of the 15N NMR chemical shifts of azoles and azines

The calculation of ¹⁵N NMR chemical shifts of 27 azoles and azines in 10 different solvents each has been carried out at the gauge including atomic orbitals density functional theory level in gas phase and applying the integral equation formalism polarizable continuum model (IEF‐PCM) and supermolecule solvation models to account for solvent effects. In the calculation of ¹⁵N NMR, chemical shifts of the nitrogen‐containing heterocycles dissolved in nonpolar and polar aprotic solvents, taking into account solvent effect is sufficient within the IEF‐PCM scheme, whereas for polar protic solvents with large dielectric constants, the use of supermolecule solvation model is recommended. A good agreement between calculated 460 values of ¹⁵N NMR chemical shifts and experiment is found with the IEF‐PCM scheme characterized by MAE of 7.1 ppm in the range of more than 300 ppm (about 2%). The best result is achieved with the supermolecule solvation model performing slightly better (MAE 6.5 ppm). Copyright © 2014 John Wiley & Sons, Ltd.     

MP2 calculation of 77Se NMR chemical shifts taking into account relativistic corrections

The main factors affecting the accuracy and computational cost of the Second‐order Möller‐Plesset perturbation theory (MP2) calculation of ⁷⁷Se NMR chemical shifts (methods and basis sets, relativistic corrections, and solvent effects) are addressed with a special emphasis on relativistic effects. For the latter, paramagnetic contribution (390–466 ppm) dominates over diamagnetic term (192–198 ppm) resulting in a total shielding relativistic correction of about 230–260 ppm (some 15% of the total values of selenium absolute shielding constants). Diamagnetic term is practically constant, while paramagnetic contribution spans over 70–80 ppm. In the ⁷⁷Se NMR chemical shifts scale, relativistic corrections are about 20–30 ppm (some 5% of the total values of selenium chemical shifts). Solvent effects evaluated within the polarizable continuum solvation model are of the same order of magnitude as relativistic corrections (about 5%). For the practical calculations of ⁷⁷Se NMR chemical shifts of the medium‐sized organoselenium compounds, the most efficient computational protocols employing relativistic Dyall's basis sets and taking into account relativistic and solvent corrections are suggested. The best result is characterized by a mean absolute error of 17 ppm for the span of ⁷⁷Se NMR chemical shifts reaching 2500 ppm resulting in a mean absolute percentage error of 0.7%. Copyright © 2015 John Wiley & Sons, Ltd.     

Towards the versatile DFT and MP2 computational schemes for 31P NMR chemical shifts taking into account relativistic corrections

The main factors affecting the accuracy and computational cost of the calculation of ³¹P NMR chemical shifts in the representative series of organophosphorous compounds are examined at the density functional theory (DFT) and second‐order Møller–Plesset perturbation theory (MP2) levels. At the DFT level, the best functionals for the calculation of ³¹P NMR chemical shifts are those of Keal and Tozer, KT2 and KT3. Both at the DFT and MP2 levels, the most reliable basis sets are those of Jensen, pcS‐2 or larger, and those of Pople, 6‐311G(d,p) or larger. The reliable basis sets of Dunning's family are those of at least penta‐zeta quality that precludes their practical consideration. An encouraging finding is that basically, the locally dense basis set approach resulting in a dramatic decrease in computational cost is justified in the calculation of ³¹P NMR chemical shifts within the 1–2‐ppm error. Relativistic corrections to ³¹P NMR absolute shielding constants are of major importance reaching about 20–30 ppm (ca 7%) improving (not worsening!) the agreement of calculation with experiment. Further better agreement with the experiment by 1–2 ppm can be obtained by taking into account solvent effects within the integral equation formalism polarizable continuum model solvation scheme. We recommend the GIAO‐DFT‐KT2/pcS‐3//pcS‐2 scheme with relativistic corrections and solvent effects taken into account as the most versatile computational scheme for the calculation of ³¹P NMR chemical shifts characterized by a mean absolute error of ca 9 ppm in the range of 550 ppm. Copyright © 2014 John Wiley & Sons, Ltd.     

DFT computational schemes for 15N NMR chemical shifts of the condensed nitrogen-containing heterocycles

A systematic density functional theory (DFT) study of the accuracy factors (functionals, basis sets, and solvent effects) for the computation of ¹⁵N NMR chemical shifts has been performed in the series of condensed nitrogen-containing heterocycles. The behavior of the most representative functionals was examined based on the benchmark calculations of ¹⁵N NMR chemical shifts in the reference set of compounds. It was found that the best agreement with experiment was achieved with OLYP functional in combination with aug-pcS-3(N)//pc-2 locally dense basis set scheme providing mean absolute error of 5.2 ppm in the range of about 300 ppm. Taking into account solvent effects was performed within a general Tomasi's polarizable continuum model scheme. It was also found that computationally demanding supermolecular solvation model computations essentially improved some “difficult” cases, as was illustrated with phenanthroline dissolved in methanol. Based on the performed calculations, some 200 unknown ¹⁵N NMR chemical shifts were predicted with a high level of confidence for about 50 real-life condensed nitrogen-containing heterocycles, which could serve as a practical guide in structural elucidation of this class of compounds.     

On the accuracy factors and computational cost of the GIAO–DFT calculation of 15N NMR chemical shifts of amides

The main factors affecting the accuracy and computational cost of Gauge‐independent Atomic Orbital–density functional theory (GIAO–DFT) calculation of ¹⁵N NMR chemical shifts in the benchmark series of 16 amides are considered. Among those are the choice of the DFT functional and basis set, solvent effects, internal reference conversion factor and applicability of the locally dense basis set (LDBS) scheme. Solvent effects are treated within the polarizable continuum model (PCM) scheme as well as at supermolecular level with solvent molecules considered in explicit way. The best result is found for Keal and Tozer's KT3 functional used in combination with Jensen's pcS‐3 basis set with taking into account solvent effects within the polarizable continuum model. The proposed LDBS scheme implies pcS‐3 on nitrogen and pc‐2 elsewhere in the molecule. The resulting mean average error for the calculated ¹⁵N NMR chemical shifts is about 6 ppm. The application of the LDBS approach tested in a series of 16 amides results in a dramatic decrease in computational cost (more than an order of magnitude in time scale) with insignificant loss of accuracy.     

GIAO‐DFT calculation of 15N NMR chemical shifts of Schiff bases: Accuracy factors and protonation effects

¹⁵N NMR chemical shifts in the representative series of Schiff bases together with their protonated forms have been calculated at the density functional theory level in comparison with available experiment. A number of functionals and basis sets have been tested in terms of a better agreement with experiment. Complimentary to gas phase results, 2 solvation models, namely, a classical Tomasi's polarizable continuum model (PCM) and that in combination with an explicit inclusion of one molecule of solvent into calculation space to form supermolecule 1:1 (SM + PCM), were examined. Best results are achieved with PCM and SM + PCM models resulting in mean absolute errors of calculated ¹⁵N NMR chemical shifts in the whole series of neutral and protonated Schiff bases of accordingly 5.2 and 5.8 ppm as compared with 15.2 ppm in gas phase for the range of about 200 ppm. Noticeable protonation effects (exceeding 100 ppm) in protonated Schiff bases are rationalized in terms of a general natural bond orbital approach.     

Calculation of <Superscript>15</Superscript>N NMR chemical shifts of amines in the framework of the density functional theory

The ¹⁵N NMR magnetic shielding constants (followed by recalculation into chemical shifts) in a representative series of amines were calculated in the framework of the density functional theory. The results were compared with experiment. The accuracy factors of this calculation (functional, basis set, solvent effects, and multistandard) were studied. Taking into account all the above factors leads to a noticeable decrease in the average absolute error in the calculation of the ¹⁵N NMR chemical shifts (from 13 to 4 ppm) in a range of their changing in the studies series of compounds of ～60 ppm (which is 6—7% in relative units).     

Towards the Accurate Calculation of 183W NMR Chemical Shifts in Polyoxometalates: The Relevance of the Structure

DFT calculations were carried out to study ¹⁸³W NMR chemical shifts in the family of the Keggin anions with formula α‐[XW₁₂O₄₀]^q? (X=B, Al, Si, P, Ga, Ge, As, Zn), in the β‐ and γ‐[SiW₁₂O₄₀]^4? geometric isomers, in the derivative Dawson anion [P₂W₁₈O₆₂]^6?, and in the most symmetrical Lindqvist [W₆O₁₉]^2? anion and its derivative [W₁₀O₃₂]^4?. In this article, we show that the geometry employed in the calculation of NMR chemical shifts in polyoxotungstates is extremely important if we want to be quantitative. Using very large basis sets of QZ4P quality and taking into account the conductor‐like screening model (COSMO) to account for solvent effects (aqueous and organic solutions), good geometries were found for the polyoxoanions. From these optimal geometries the ¹⁸³W NMR chemical shifts were computed with the more standard basis sets of TZP quality and including spin–orbit corrections inside the zero‐order regular approximation (ZORA) to describe the relativistic effects of the internal electrons. With this strategy the mean absolute error between experimental and theoretical values was found to be less than 10 ppm, which is similar to the experimental error. We also discuss how the geometry of the polyoxoanion influences on the shielding.     

Computational NMR Spectroscopy of Transition‐Metal/Nitroimidazole Complexes: Theoretical Investigation of Potential Radiosensitizers

The computed chemical shifts of transition‐metal complexes with dimetridazole (= 1,2‐dimethyl‐5‐nitro‐1H‐imidazole; 1 ), a prototypical nitro‐imidazole‐based radiosensitizer, are reported at the GIAO‐BP86 and ‐B3LYP levels for BP86/ECP1‐optimized geometries. These complexes comprise [MCl₂( 1 )₂] (M = Zn, Pd, Pt), [RuCl₂(DMSO)₂( 1 )₂], and [Rh₂(O₂CMe)₄( 1 )₂]. Available δ(¹H) and δ(¹⁵N) values, and Δδ(¹H) and Δδ(¹⁵N) coordination shifts are well‐reproduced theoretically, provided solvation and relativistic effects are taken into account by means of a polarizable continuum model and suitable methods including spin–orbit (SO) coupling, respectively. These effects are particularly important for the metal‐coordinated N‐atom, where the contributions from solvation and relativity can affect δ(¹⁵N) and Δδ(¹⁵N) values up to 10–20 ppm. The ¹⁹⁵Pt chemical shifts of cis‐ and trans‐[PtCl₂( 1 )₂] are well‐reproduced using the zero‐order regular approximation including SO coupling (ZORA‐SO). Predictions are reported for ⁹⁹Ru and ¹⁰³Rh chemical shifts, which suggest that these metal centers could be used as additional, sensitive NMR probes in their complexes with nitro‐imidazoles.     

On the accuracy of the GIAO‐DFT calculation of 15N NMR chemical shifts of the nitrogen‐containing heterocycles – a gateway to better agreement with experiment at lower computational cost

The main factors affecting the accuracy and computational cost of the gauge‐independent atomic orbital density functional theory (GIAO‐DFT) calculation of ¹⁵N NMR chemical shifts in the representative series of key nitrogen‐containing heterocycles – azoles and azines – have been systematically analyzed. In the calculation of ¹⁵N NMR chemical shifts, the best result has been achieved with the KT3 functional used in combination with Jensen's pcS‐3 basis set (GIAO‐DFT‐KT3/pcS‐3) resulting in the value of mean absolute error as small as 5 ppm for a range exceeding 270 ppm in a benchmark series of 23 compounds with an overall number of 41 different ¹⁵N NMR chemical shifts. Another essential finding is that basically, the application of the locally dense basis set approach is justified in the calculation of ¹⁵N NMR chemical shifts within the 3–4 ppm error that results in a dramatic decrease in computational cost. Based on the present data, we recommend GIAO‐DFT‐KT3/pcS‐3//pc‐2 as one of the most effective locally dense basis set schemes for the calculation of ¹⁵N NMR chemical shifts. Copyright © 2014 John Wiley & Sons, Ltd.     

Four‐component relativistic DFT calculations of 77Se NMR chemical shifts: A gateway to a reliable computational scheme for the medium‐sized organoselenium molecules

A versatile high‐accuracy computational scheme for the ⁷⁷Se nuclear magnetic resonance (NMR) chemical shifts of the medium‐sized organoselenium compounds is suggested within a framework of a full four‐component relativistic density functional theory (DFT). The main accuracy factors (DFT functionals, relativistic geometry, vibrational corrections, and solvent effects) are addressed. The best result is achieved with NMR‐oriented KT2 functional of Keal–Tozer characterized by a fairly small error of only 30 ppm for the span of about 1700 ppm (<2%). © 2015 Wiley Periodicals, Inc.     

Normal halogen dependence of 13C NMR chemical shifts of halogenomethanes revisited at the four‐component relativistic level

The ‘Normal Halogen Dependence’ of ¹³C NMR chemical shifts in the series of halogenomethanes is revisited at the four‐component relativistic level. Calculations of ¹³C NMR chemical shifts of 70 halogenomethanes have been carried out at the density functional theory (DFT) and MP2 levels with taking into account relativistic effects using the four‐component relativistic theory of Dirac‐Coulomb within the different computational methods (4RPA, 4OPW91) and hybrid computational schemes (MP2 + 4RPA, MP2 + 4OPW91). The most efficient computational protocols are derived for practical purposes. Relativistic shielding effect reaches as much as several hundreds of ppm for heavy halogenomethanes, and to account for this effect in comparison with experiment at the qualitative level, relativistic Dyall's basis sets of triple‐zeta quality or higher are to be used within the framework of the four‐component relativistic theory taking into account solvent effects. Relativistic geometrical optimization (as compared with the non‐relativistic level) is essential for the molecules containing at least two iodines at one carbon atom. Copyright © 2016 John Wiley & Sons, Ltd.     

15N chemical shifts of a series of isatin oxime ethers and their corresponding nitrone isomers

In this article, we describe the characteristic ¹⁵N chemical shifts of isatin oxime ethers and their isomer nitrone. These oxime ethers and nitrones are the alkylation reaction products of isatin oximes. In our study, the ¹⁵N chemical shifts observed in these oxime ethers were in the 402–408 (or 22–28) ppm range, although those for their corresponding nitrone series were in the 280–320 (or ?100 to ?60) ppm range. This remarkable difference in ¹⁵N NMR chemical shift values could potentially be used to determine the O‐ versus N‐alkylation of oximes, even when only one isomer is available. In this paper, the differences in ¹⁵N NMR chemical shifts serve as the basis for a discussion about how to distinguish both regioisomers derived from the oximes alkylation. Copyright © 2010 John Wiley & Sons, Ltd.     

Relativistic effect of iodine in <Superscript>13</Superscript>C NMR chemical shifts of iodomethanes from quantum chemical calculations within the framework of the full four-component relativistic Dirac—Coulomb scheme

¹³C NMR shielding constants (chemical shifts) of iodomethanes were calculated within the framework of the full four-component relativistic Dirac—Coulomb scheme. As the number of iodine atoms in the molecule increases, the relativistic counterpart of the ¹³C NMR chemical shift increases from a few tens to several hundreds of ppm. Calculations of ¹³C NMR chemical shifts of organoiodine compounds should be performed at the relativistic level using relativistic Dyall’s basis sets dyall.vXz and dyall.xvXz (x = a, c, ac, ae; X = 2, 3, 4) of at least triple-zeta quality or at the correlated non-relativistic level taking into account relativistic corrections. Solvent effects are not of prime importance; however, taking into account the solvent corrections causes the mean absolute error of determination of the ¹³C NMR chemical shifts to decrease by 1—2 ppm.     

Solid‐state 13C, 15N and 29Si NMR characterization of block copolymers with CO2 capture properties

Natural abundance solid‐state multinuclear (¹³C, ¹⁵N and ²⁹Si) cross‐polarization magic‐angle‐spinning NMR was used to study structures of three block copolymers based on polyamide and dimethylsiloxane and two polyamides, one of which including ferrocene in its structure. Assignment of most of the resonance lines in ¹³C, ¹⁵N and ²⁹Si cross‐polarization magic‐angle‐spinning NMR spectra were suggested. A comparative analysis of ¹³C isotropic chemical shifts of polyamides with and without ferrocene has revealed a systematic shift towards higher δ ‐values (de‐shielding) explained as the incorporation of paramagnetic ferrocene into the polyamide backbone. In addition, the ¹³C NMR resonance lines for ferrocene‐based polyamide were significantly broadened, because of paramagnetic effects from ferrocene incorporated in the structure of this polyamide polymer. Single resonance lines with chemical shifts ranging from 88.1 to 91.5 ppm were observed for ¹⁵N sites in all of studied polyamide samples. ²⁹Si chemical shifts were found to be around ?22.4 ppm in polydimethylsiloxane samples that falls in the range of chemical shifts for alkylsiloxane compounds. The CO₂ capture performance of polyamide‐dimethylsiloxane‐based block copolymers was measured as a function of temperature and pressure. The data revealed that these polymeric materials have potential to uptake CO₂ (up to 9.6 cm³ g^?1) at ambient pressures and in the temperature interval 30–40 °C. Copyright © 2016 John Wiley & Sons, Ltd.     

Theoretical and experimental study of 15N NMR protonation shifts

A combined theoretical and experimental study revealed that the nature of the upfield (shielding) protonation effect in ¹⁵N NMR originates in the change of the contribution of the sp²‐hybridized nitrogen lone pair on protonation resulting in a marked shielding of nitrogen of about 100 ppm. On the contrary, for amine‐type nitrogen, protonation of the nitrogen lone pair results in the deshielding protonation effect of about 25 ppm, so that the total deshielding protonation effect of about 10 ppm is due to the interplay of the contributions of adjacent natural bond orbitals. A versatile computational scheme for the calculation of ¹⁵N NMR chemical shifts of protonated nitrogen species and their neutral precursors is proposed at the density functional theory level taking into account solvent effects within the supermolecule solvation model. Copyright © 2015 John Wiley & Sons, Ltd.     

15N nuclear magnetic resonance spectroscopy. Natural abundance 15N NMR of monosubstituted indoles

¹⁵N chemical shifts of twenty-four substituted indoles have been determined in natural abundance (in organic solvents) using Fourier transform NMR. The overall chemical shift range is 27 ppm, with groups in the 2-, 3- and 5-ring positions showing the largest substituent effects. Substituents capable of resonance interaction with the indole nitrogen give shifts in the expected directions but they cannot be correlated with known substituent parameters. Compounds measured in DMSO give 0·2 to 10·2 ppm downfield shifts with respect to the same compound measured in CDCl₃. ¹³C NMR data for previously unreported compounds are also reported.     

Combined experimental and computational study of W(II), Ru(II), Pt(IV) and Cu(I) amine and amido complexes using 15N NMR spectroscopy

Using 2D proton-coupled gHSQC pulse sequences in addition to 1D ¹⁵N NMR experiments of ¹⁵N labeled systems, ¹⁵N NMR chemical shifts of a range of transition metal amido and amine complexes were determined. Tungsten(II), ruthenium(II), platinum(IV) and copper(I) complexes with aniline and their anilido variants were studied and compared to free aniline, lithium anilido and anilinium tetrafluoroborate. Upon coordination of aniline to transition metals, upfield chemical shifts of 20–60 ppm were observed. Deprotonation of the amine complexes to form amido complexes resulted in downfield chemical shifts of 40–60 ppm for all of the complexes except for the tungsten d⁴ system. For the tungsten(II) complexes, the cationic aniline complex displayed a downfield shift of approximately 56 ppm relative to the neutral anilido complex. The change in chemical shift for amine to amido conversion is proposed to depend on the ability of the amido ligand to π-bond with the metal center, which influences the magnitude of the paramagnetic screening term.     

Predicting 15N amide chemical shifts in proteins. I. An additive model for the backbone contribution

Because proteins adopt unique structures, chemically identical nuclei in proteins exhibit different chemical shifts. Amide ¹⁵N chemical shifts have been shown to vary over 20 ppm. The cause of these chemical shift inequivalencies is the different intra‐ and intermolecular interactions that individual nuclei experience at different locations in the protein structure. These chemical shift inequivalencies can be described as structural shifts, the difference between the actual chemical shift and the random coil chemical shift. As a first step toward the prediction of these amide ¹⁵N structural shifts, calculations have been carried out on acetyl‐glycine‐methyl amide to examine how a neighboring peptide group influences the amide ¹⁵N structural shifts. The ϕ,ψ dihedral angle space is completely surveyed, while all other geometrical variables are held fixed, to isolate the effect of the backbone conformation. Similar calculations for a limited number of conformations of acetyl‐glycine‐glycine‐methyl amide were carried out, where the effects of the two terminal peptide groups on the central amide ¹⁵N structural shift are examined. It is shown that the effect of the two adjacent groups can be accurately modeled by combining their individual effects additively. This provides a quite simple method to predict the backbone influence on amide ¹⁵N structural shifts in proteins. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 366–372, 2001