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.     

On the long‐range relativistic effects in the 15N NMR chemical shifts of halogenated azines

Long‐range β‐ and γ‐relativistic effects of halogens in ¹⁵N NMR chemical shifts of 20 halogenated azines (pyridines, pyrimidines, pyrazines, and 1,3,5‐triazines) are shown to be unessential for fluoro‐, chloro‐, and bromo‐derivatives (1–2 ppm in average). However, for iodocontaining compounds, β‐ and γ‐relativistic effects are important contributors to the accuracy of the ¹⁵N calculation. Taking into account long‐range relativistic effects slightly improves the agreement of calculation with experiment. Thus, mean average errors (MAE) of ¹⁵N NMR chemical shifts of the title compounds calculated at the non‐relativistic and full 4‐component relativistic levels in gas phase are accordingly 7.8 and 5.5 ppm for the range of about 150 ppm. Taking into account solvent effects within the polarizable continuum model scheme marginally improves agreement of computational results with experiment decreasing MAEs from 7.8 to 7.4 ppm and from 5.5 to 5.3 ppm at the non‐relativistic and relativistic levels, respectively. The best result (MAE: 5.3 ppm) is achieved at the 4‐component relativistic level using Keal and Tozer's KT3 functional used in combination with Dyall's relativistic basis set dyall.av3z with taking into account solvent effects within the polarizable continuum solvation model. The long‐range relativistic effects play a major role (of up to dozen of parts per million) in ¹⁵N NMR chemical shifts of halogenated nitrogen‐containing heterocycles, which is especially crucial for iodine derivatives. This effect should apparently be taken into account for practical purposes.     

Halogen effect on structure and 13C NMR chemical shift of 3,6‐disubstituted‐N‐alkyl carbazoles

Structures of selected 3,6‐dihalogeno‐N‐alkyl carbazole derivatives were calculated at the B3LYP/6‐311++G(3df,2pd) level of theory, and their ¹³C nuclear magnetic resonance (NMR) isotropic shieldings were predicted using density functional theory (DFT). The model compounds contained 9H, N‐methyl and N‐ethyl derivatives. The relativistic effect of Br and I atoms on nuclear shieldings was modeled using the spin–orbit zeroth‐order regular approximation (ZORA) method. Significant heavy atom shielding effects for the carbon atom directly bonded with Br and I were observed (~?10 and ~?30 ppm while the other carbon shifts were practically unaffected). The decreasing electronegativity of the halogen substituent (F, Cl, Br, and I) was reflected in both nonrelativistic and relativistic NMR results as decreased values of chemical shifts of carbon atoms attached to halogen (C3 and C6) leading to a strong sensitivity to halogen atom type at 3 and 6 positions of the carbazole ring. The predicted NMR data correctly reproduce the available experimental data for unsubstituted N‐alkylcarbazoles. Copyright © 2013 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.     

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.     

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.     

Full four‐component relativistic calculations of the one‐bond 77Se–13C spin‐spin coupling constants in the series of selenium heterocycles and their parent open‐chain selenides

Four‐component relativistic calculations of ⁷⁷Se–¹³C spin–spin coupling constants have been performed in the series of selenium heterocycles and their parent open‐chain selenides. It has been found that relativistic effects play an essential role in the selenium–carbon coupling mechanism and could result in a contribution of as much as 15–25% of the total values of the one‐bond selenium–carbon spin‐spin coupling constants. In the overall contribution of the relativistic effects to the total values of ¹J(Se,C), the scalar relativistic corrections (negative in sign) by far dominate over the spin‐orbit ones (positive in sign), the latter being of less than 5%, as compared to the former (ca 20%). A combination of nonrelativistic second‐order polarization propagator approach (CC2) with the four‐component relativistic density functional theory scheme is recommended as a versatile tool for the calculation of ¹J(Se,C). Solvent effects in the values of ¹J(Se,C) calculated within the polarizable continuum model for the solvents with different dielectric constants (ε 2.2–78.4) are next to negligible decreasing negative ¹J(Se,C) in absolute value by only about 1 Hz. The use of the locally dense basis set approach applied herewith for the calculation of ⁷⁷Se–¹³C spin‐spin coupling constants is fully justified resulting in a dramatic decrease in computational cost with only 0.1–0.2‐Hz loss of accuracy. Copyright © 2014 John Wiley & Sons, Ltd.     

One‐bond 29Si‐1H spin‐spin coupling constants in the series of halosilanes: benchmark SOPPA and DFT calculations,relativistic effects,and vibrational corrections

A number of most representative second order polarization propagator approach (SOPPA) based wavefunction methods, SOPPA, SOPPA(CC2) and SOPPA(CCSD), and density functional theory (DFT) based methods, B3LYP, PBE0, KT2, and KT3, have been benchmarked in the calculation of the one‐bond ²⁹Si‐¹H spin‐spin coupling constants in the series of halosilanes SiH_nX_4?n (X = F, Cl, Br, I), both at the non‐relativistic and full four‐parameter Dirac's relativistic levels taking into account vibrational corrections. At the non‐relativistic level, the wavefunction methods showed much better results as compared with those of DFT. At the DFT level, out of four tested functionals, the Perdew, Burke, and Ernzerhof's PBE0 showed best performance. Taking into account, relativistic effects and vibrational corrections noticeably improves wavefunction methods results, but generally worsens DFT results. Copyright © 2013 John Wiley & Sons, Ltd.     

Quantum‐chemical analyses of aromaticity,UV spectra,and NMR chemical shifts in plumbacyclopentadienylidenes stabilized by Lewis bases

We carried out a series of zeroth‐order regular approximation (ZORA)‐density functional theory (DFT) and ZORA‐time‐dependent (TD)‐DFT calculations for molecular geometries, NMR chemical shifts, nucleus‐independent chemical shifts (NICS), and electronic transition energies of plumbacyclopentadienylidenes stabilized by several Lewis bases, (Ph)₂(^tBuMe₂Si)₂C₄PbL₁L₂ (L₁, L₂ = tetrahydrofuran, Pyridine, N‐heterocyclic carbene), and their model molecules. We mainly discussed the Lewis‐base effect on the aromaticity of these complexes. The NICS was used to examine the aromaticity. The NICS values showed that the aromaticity of these complexes increases when the donation from the Lewis bases to Pb becomes large. This trend seems to be reasonable when the 4n‐Huckel rule is applied to the fractional π‐electron number. The calculated ¹³C‐ and ²⁰⁷Pb‐NMR chemical shifts and the calculated UV transition energies reasonably reproduced the experimental trends. We found a specific relationship between the ¹³C‐NMR chemical shifts and the transition energies. As we expected, the relativistic effect was essential to reproduce a trend not only in the ²⁰⁷Pb‐NMR chemical shifts and J[Pb‐C] but also in the ¹³C‐NMR chemical shifts of carbons adjacent to the lead atom. © 2014 Wiley Periodicals, Inc.     

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.     

New implementation of molecular double point‐group symmetry in four‐component relativistic Gaussian‐type spinors

A new, practical implementation of double‐group symmetry to relativistic Gaussian spinors is presented for four‐component relativistic molecular calculations. We show that the systematic adaptability to irreducible representations under arbitrary point‐group symmetry, as well as Kramers (time‐reversal) symmetry, is inherent in the present basis spinors, which possess the analytic structure of Dirac atomic spinors. The implementation of double‐group symmetry entails significant computational efficiencies in the relativistic second‐order Møller–Plesset perturbation calculation on Au₂ and the density functional theory (DFT) calculation with the B3LYP functional on octahedral UF₆, in which the highest symmetries used are, respectively, C and D. The four‐component B3LYP equilibrium geometry of UF₆ is reported. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

An Investigation of Lanthanum Coordination Compounds by Using Solid‐State 139La NMR Spectroscopy and Relativistic Density Functional Theory

Lanthanum‐139 NMR spectra of stationary samples of several solid La^III coordination compounds have been obtained at applied magnetic fields of 11.75 and 17.60 T. The breadth and shape of the ¹³⁹La NMR spectra of the central transition are dominated by the interaction between the ¹³⁹La nuclear quadrupole moment and the electric field gradient (EFG) at that nucleus; however, the influence of chemical‐shift anisotropy on the NMR spectra is non‐negligible for the majority of the compounds investigated. Analysis of the experimental NMR spectra reveals that the ¹³⁹La quadrupolar coupling constants (C_Q) range from 10.0 to 35.6 MHz, the spans of the chemical‐shift tensor (Ω) range from 50 to 260 ppm, and the isotropic chemical shifts (δ_iso) range from ?80 to 178 ppm. In general, there is a correlation between the magnitudes of C_Q and Ω, and δ_iso is shown to depend on the La coordination number. Magnetic‐shielding tensors, calculated by using relativistic zeroth‐order regular approximation density functional theory (ZORA‐DFT) and incorporating scalar only or scalar plus spin–orbit relativistic effects, qualitatively reproduce the experimental chemical‐shift tensors. In general, the inclusion of spin–orbit coupling yields results that are in better agreement with those from the experiment. The magnetic‐shielding calculations and experimentally determined Euler angles can be used to predict the orientation of the chemical‐shift and EFG tensors in the molecular frame. This study demonstrates that solid‐state ¹³⁹La NMR spectroscopy is a useful characterization method and can provide insight into the molecular structure of lanthanum coordination compounds.     

First example of the correlated calculation of the one‐bond tellurium–carbon spin–spin coupling constants: Relativistic effects,vibrational corrections,and solvent effects

This work reports on the comprehensive calculation of the NMR one‐bond spin–spin coupling constants (SSCCs) involving carbon and tellurium, ¹J(¹²⁵Te,¹³C), in four representative compounds: Te(CH₃)₂, Te(CF₃)₂, Te(C?CH)₂, and tellurophene. A high‐level computational treatment of ¹J(¹²⁵Te,¹³C) included calculations at the SOPPA level taking into account relativistic effects evaluated at the 4‐component RPA and DFT levels of theory, vibrational corrections, and solvent effects. The consistency of different computational approaches including the level of theory of the geometry optimization of tellurium‐containing compounds, basis sets, and methods used for obtainig spin–spin coupling values have also been discussed in view of reproducing the experimental values of the tellurium–carbon SSCCs. Relativistic corrections were found to play a major role in the calculation of ¹J(¹²⁵Te,¹³C) reaching as much as almost 50% of the total value of ¹J(¹²⁵Te,¹³C) while relativistic geometrical effects are of minor importance. The vibrational and solvent corrections account for accordingly about 3–6% and 0–4% of the total value. It is shown that taking into account relativistic corrections, vibrational corrections and solvent effects at the DFT level essentially improves the agreement of the non‐relativistic theoretical SOPPA results with experiment. © 2016 Wiley Periodicals, Inc.     

1H NMR analysis of O‐methyl‐inositol isomers: a joint experimental and theoretical study

Density functional theory (DFT) calculations of ¹H NMR chemical shifts for l ‐quebrachitol isomers were performed using the B3LYP functional employing the 6‐31G(d,p) and 6‐311 + G(2d,p) basis sets. The effect of the solvent on the B3LYP‐calculated NMR spectrum was accounted for using the polarizable continuum model. Comparison is made with experimental ¹H NMR spectroscopic data, which shed light on the average uncertainty present in DFT calculations of chemical shifts and showed that the best match between experimental and theoretical B3LYP ¹H NMR profiles is a good strategy to assign the molecular structure present in the sample handled in the experimental measurements. Among four plausible O‐methyl‐inositol isomers, the l ‐quebrachitol 2a structure was unambiguously assigned based only on the comparative analysis of experimental and theoretical ¹H NMR chemical shift data. The B3LYP infrared (IR) spectrum was also calculated for the four isomers and compared with the experimental data, with analysis of the theoretical IR profiles corroborating assignment of the 2a structure. Therefore, it is confirmed in this study that a combined experimental/DFT spectroscopic investigation is a powerful tool in structural/conformational analysis studies. Copyright © 2012 John Wiley & Sons, Ltd.     

The Rearrangement of 2,2‐Diphenyl‐1‐[(E)‐2‐phenylethenyl]cyclopropane to 3,4,4‐Triphenylcyclopent‐1‐ene: a DFT Analysis

A computational study on the rearrangement of 2,2‐diphenyl‐1‐[(E)‐2‐phenylethenyl]cyclopropane ( 1 ) is presented, using density functional theory (DFT), (U)B3LYP with the 6‐31G* basis set (DFT1) and (U)M05‐2X with the 6‐311+G** basis set (DFT2). In agreement with a biradical character of the transition structure (TS) or intermediate, the potential‐energy hypersurface is lowered by the influence of three conjugated Ph groups. Surprisingly, two conformations of the geminal diphenyl group (different twist angles) induce two different minimum‐energy pathways for the rearrangement. Independent of the functional used, the first hypersurface harbors true biradical intermediates, whereas the second energy surface is a flat, slightly ascending slope from the starting material to the TS. The functional (U)M05‐2X with the basis set 6‐311+G** provides realistic energies which seem to be close to experiment. The activation energy for racemization of enantiomers of 1 is lower than that of rearrangement by 2.5 kcal mol^?1, in agreement with experiment.     

Density functional calculations of 19F and 235U NMR chemical shifts in uranium (VI) chloride fluorides UF6−nCln: Influence of the relativistic approximation and role of the exchange‐correlation functional

Relativistic density functional theory (DFT) has been applied to the calculation of the ¹⁹F nuclear magnetic resonance (NMR) chemical shifts of the title compounds. It is shown that, while large‐core effective core potentials (ECP) fail completely for the calculation of ligand NMR chemical shifts in uranium compounds, small‐core ECPs are a valid relativistic method for this purpose. In an earlier study of the same systems, certain differences between theory and experiment had been observed, for instance, in the relative chemical shift of the A₄ and X sites in UF₅Cl. The reason for these deviations has been investigated further in the current paper. By comparing different relativistic methods, it is shown that the relativistic approximation is not responsible for these deviations. The role of the approximation to the exchange‐correlation (XC) functional of DFT has been probed, and generalized gradient approximations (GGA) as well as hybrid DFT methods have been investigated. None of these methods corrects the mentioned errors. It is argued that the neglect of environmental factors (solvent effects) remains as a possible error source, although the approximate XC functional appears to be the more likely cause of the problem. ²³⁵U NMR shieldings and chemical shifts have been calculated, and the trends predicted earlier have been confirmed. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

All‐electron relativistic computations on the low‐lying electronic states,bond length,and vibrational frequency of CeF diatomic molecule with spin–orbit coupling effects

Ab initio all‐electron computations have been carried out for Ce⁺ and CeF, including the electron correlation, scalar relativistic, and spin–orbit coupling effects in a quantitative manner. First, the n‐electron valence state second‐order multireference perturbation theory (NEVPT2) and spin–orbit configuration interaction (SOCI) based on the state‐averaged restricted active space multiconfigurational self‐consistent field (SA‐RASSCF) and state‐averaged complete active space multiconfigurational self‐consistent field (SA‐CASSCF) wavefunctions have been applied to evaluations of the low‐lying energy levels of Ce⁺ with [Xe]4f¹5d¹6s¹ and [Xe]4f¹5d² configurations, to test the accuracy of several all‐electron relativistic basis sets. It is shown that the mixing of quartet and doublet states is essential to reproduce the excitation energies. Then, SA‐RASSCF(CASSCF)/NEVPT2 + SOCI computations with the Sapporo(‐DKH3)‐2012‐QZP basis set were carried out to determine the energy levels of the low‐lying electronic states of CeF. The calculated excitation energies, bond length, and vibrational frequency are shown to be in good agreement with the available experimental data. © 2018 Wiley Periodicals, Inc.     

Theoretical and experimental 1H, 13C and 15N NMR studies of N‐alkylation of substituted tetrazolo[1,5‐a]pyridines

The ¹⁵N as well as ¹H and ¹³C chemical shifts of nine substituted tetrazolopyridines and their corresponding tetrazolopyridinium salts have been determined by using NMR spectroscopy at the natural abundance level of all nuclei in CD₃CN. In this paper, we report, for the first time, the N‐alkylation reaction of electron deficient tetrazolopyridines. The treatment of tetrazolopyridines 5–13 with one equivalent of trialkyloxonium tetrafluoroborate leads to a mixture of two isomers, i.e. N3‐ and N2‐alkyl tetrazolo[1,5‐a]pyridinium salts. It has been observed that the N3‐isomer is always the major isomer, except in the case of the CF₃ substituent, where the two isomers are obtained in the same amount. The quaternary tetrazolopyridinium nitrogen N3 is shielded by around 100 ppm (parts per million) with respect to the parent tetrazolopyridine. Experimental data are interpreted by means of density functional theory (DFT) calculations, including solvent‐induced effects, within the conductor‐like polarizable continuum model (CPCM). Good agreements between theoretical and experimental ¹H, ¹³C and ¹⁵N NMR were found. The combination of multinuclear magnetic resonance spectroscopy with gauge including atomic orbital (GIAO) DFT calculations is a powerful tool in the structural elucidation for both neutral and cationic heterocycles and in the determination of the orientation of N‐alkylation of tetrazolopyridines. Copyright © 2009 John Wiley & Sons, Ltd.