Theoretical Study of 13C Chemical Shifts Structures of Some ortho-Substituted 3-Anilino-2-nitrobenzo[b]thiophenes and 2-Anilino-3-nitrobenzo[b]thiophenes and Comparison with Experimental

¹³C NMR chemical shifts have been calculated for structures of some substituted 3‐anilino‐2‐nitrobenzo‐[b]thiophenes ( 2 o) and 2‐anilino‐3‐nitrobenzo[b]thiophenes ( 3 o) derivatives containing OH, NH₂, OMe, Me, Et, H, F, Cl and Br. The molecular structures were fully optimized using B3LYP/6‐31G(d,p). The calculation of the ¹³C shielding tensors employed the GAUSSIAN 03 implementation of the gauge‐including atomic orbital (GIAO) and continuous set of gauge transformations (CSGT) by using 6‐311++G(d,p) basis set at density functional levels of theories (DFT). The isotropic and the anisotropy parameters of chemical shielding for all compounds are calculated. The predicted ¹³C chemical shifts are derived from equation δ=δ₀+δ where δ is the chemical shift, δ is the absolute shielding, and δ₀ is the absolute shielding of the standard TMS. Excellent linear relationships have been observed between experimental and calculated ¹³C NMR chemical shifts for all derivatives     

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.     

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.     

GIAO/DFT evaluation of 13C NMR chemical shifts of selected acetals based on DFT optimized geometries

DFT/B3LYP calculations of the ground-state conformation of eight cyclic and acyclic acetals are presented and compared with experimental data. Results of single-point GIAO/DFT calculations at five different levels of theory show that isotropic shieldings need to be empirically scaled to achieve agreement with experimental chemical shifts. Statistical evaluation of data indicates that the most accurate prediction of 13C chemical shifts is achieved at the MPW1PW91/6-311G** level of theory. An empirical equation describing the relationship between delta values and shielding constants is postulated. This equation has been applied to the non-chair ground-state conformation of the six-membered acetonide and to the conformationally flexible benzodioxonine derivative. The agreement observed between the experimental and predicted chemical shifts shows that calculations at the MPW1PW91/6-311G** level of theory are adequate for addressing questions of conformation.     

DFT Calculations of 1H NMR Chemical Shifts of Geometric Isomers of Conjugated Linolenic Acids,Hexadecatrienyl Pheromones,and Model Triene-Containing Compounds: Structures in Solution and Revision of NMR Assignments

A DFT study of the ¹H NMR chemical shifts, δ(¹H), of geometric isomers of 18:3 conjugated linolenic acids (CLnAs), hexadecatrienyl pheromones, and model triene-containing compounds is presented, using standard functionals (B3LYP and PBE0) as well as corrections for dispersion interactions (B3LYP-D3, APFD, M06–2X and ωB97XD). The results are compared with literature experimental δ(¹H) data in solution. The closely spaced “inside” olefinic protons are significantly more deshielded due to short-range through-space H^…H steric interactions and appear close to or even beyond δ-values of aromatic systems. Several regularities of the computational δ(¹H) of the olefinic protons of the conjugated double bonds are reproduced very accurately for the lowest-energy DFT-optimized single conformer for all functionals used and are in very good agreement with experimental δ(¹H) in solution. Examples are provided of literature studies in which experimental resonance assignments deviate significantly from DFT predictions and, thus, should be revised. We conclude that DFT calculations of ¹H chemical shifts of trienyl compounds are powerful tools (i) for the accurate prediction of δ(¹H) even with less demanding functionals and basis sets; (ii) for the unequivocal identification of geometric isomerism of conjugated trienyl systems that occur in nature; (iii) for tackling complex problems of experimental resonance assignments due to extensive signal overlap; and (iv) for structure elucidation in solution.     

(XNR)_4立方簇合物结构与稳定性的DFT研究

用密度泛函理论(DFT),在B3LYP/CEP-121G水平上,首次对(XNR)_4 [X = C, Si,Ge,Sn,Pb;R = H,CH_3,C(CH_3)_3,Si(CH_3)_3,C_6H_5]立方簇合物的 几何构型、电子结构、振动光谱和化学键性质进行了研究,并对其分子碎片(XNR) _2进行了相同方法的优化计算。结果表明,(CNR)_4对应原5种结构均不能稳定存在 。(XNH)_4对应原5种结构稳定性较低。对于同一取代基R,簇合物的稳定性大小顺 序为：(PbNR)_4>(SnNR)_4>(GeNR)_4>(SiNR)_4>(CNR)_4。     

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.     

Syntheses and DFT/GIAO/NMR Study of Polyhedral Selenaboranes: The First Selenahexaborane(5) closo‐1‐SeB5Cl5

The six‐vertex closo‐1‐SeB₅Cl₅ ( 1 ) and the 12‐vertex closo‐1‐SeB₁₁Cl₁₁ ( 2 ) selenaboranes were prepared in low yields from the co‐pyrolysis reaction of B₂Cl₄ with Se₂Cl₂ at 330 °C. Both compounds are sublimable, colorless solids and were tentatively characterized by ¹¹B NMR and high‐resolution mass spectroscopy. DFT/GIAO/NMR calculations support octahedral and icosahedral arrangements for 1 and 2 , respectively, as expected due to their closo electron counts. Compound 1 represents the first example of a polyhedral selenaborane with a cluster size smaller than 10 vertices.     

A Comparison between the Experimental and Theoretical Investigations on Carbon‐13 Isotropic Shielding Constants of Some Tripodal Tetraamine Ligands

The direct implementation of GIAO and CSGT methods for calculation of ¹³C isotropic shielding constants of fully protonated forms of six tripodal tetraamine ligands tren, pee, ppe, tpt, epb and ppb at the Hartree‐Fock level of theory are presented. The shielding constants were determined using hybrid methods (including a mixture of Hartree‐Fock exchange and DFT exchange‐correlation) and are close to the experimental data. A splayed‐like conformation was considered for fully protonated forms of all ligands, and it was confirmed that this is the most stable conformation for the latter form of such ligands. A good linear correlation between the calculated chemical shielding at gas‐phase and experimental shift values in D₂O solution was obtained.     

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.     

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.     

Influence of the H/F replacement on the homoaromaticity of homotropylium ion: a GIAO/DFT theoretical study

The problem of homoaromaticity in mono-, di- and polyfluorinated- homotropylium cations is addressed by the B3LYP/6-311++G** DFT method. The energetic, structural and magnetic criteria are used for this purpose. They convincingly show that the ground state equilibrium species are aromatic, or in other words that the homoaromaticity is preserved by the (poly)fluorination. In contrast, a considerable decrease in the aromatic stabilization is observed in the transition structures (TS). According to the NICS(0) index, they vary form strongly antiaromatic, via weakly and non-aromatic to slightly aromatic transition states. However, the hierarchy of the aromaticity in fluorinated homotropylium ions predicted by NICS(0) is completely unrelated to that obtained by using the energy criterion assuming a kinetic definition of aromaticity. On the other hand the latter is closely related to geometric parameters of the equilibrium and transition structures.     

DFT studies of the structure, vibrational and NMR spectra of 2-amino-pyridine betaine monohydrate

Four of the most stable conformers of 2-amino-pyridine betaine (1-carboxymethyl-2-amino-pyridinium inner salt) monohydrates, 2-NH₂PB·H₂O, and one anhydrous were analyzed by the B3LYP/6-31G(d,p) calculations and compared with the X-ray data. Two types of optimized conformers can be distinguished: (a) with NH₂ and COO groups and (b) an imino tautomer with NH and COOH groups. A common feature of the optimized molecules are intramolecular hydrogen bonds between the COO⁻ and H₂N or COOH and HN groups. In the crystal both NH₂ and COO groups participate in intermolecular hydrogen bonds. The probable assignments of the anharmonic experimental solid state vibrational frequencies of 2-NH₂PB·H₂O and 2-ND₂PB·D₂O (conformer 2) based on the calculated B3LYP/6-31G(d,p) harmonic frequencies have been made. Correlations between experimental chemical shifts for 2-NH₂PB, its hydrochloride and 1-carboxyethyl-2-amino-pyridinium inner salt (¹³C and ¹H in D₂O) and GIAO/B3LYP/6-31G(d,p) calculated isotropic shielding constants, δ_exp=a+bσ_calc, are reported. Good linear regression between experimental and theoretical results for ¹³C was obtained. Only in 2-NH₂PB the hydrogen at -position is outside the linear correlation.     

Solid-state NMR studies and DFT calculations of flavonoids: baicalein, baicalin and wogonoside

Three flavonoids of pharmaceutical importance-baicalein, baicalin, and wogonoside-were isolated from a Chinese medicinal plant Scutellaria baicalensis Georgi and studied by 13C NMR in solution and solid state. Two-dimensional (2D) NMR spectroscopy in the liquid phase and dipolar dephasing (DD) experiments in magic-angle spinning (MAS) spectra enabled the assignment of 13C resonances. The cross-polarization (CP) time constants T(CH) and relaxation times T(H) (1rho) were obtained from the variable-contact time experiments. The principal elements of the 13C chemical shift tensor were determined in the spectra recorded under slow sample spinning (2 kHz) using phase-adjusted spinning sideband (PASS)-2D NMR technique, and were verified by density functional theory gauge-independent atomic orbital (DFT GIAO) calculations of shielding constants. Analysis of the 13C delta(ii) and comparison with shielding parameters calculated for different conformers of compounds 1-3 enabled the selection of the most reliable geometry in the solid phase. In all three compounds, an intramolecular hydrogen bond C5--OH...=C4 is formed; the existence of baicalein and baicalin with 'anticlockwise' orientation of OH groups is more probable.     

Substituent effects on 13C chemical shifts of ketenimines: a GIAO/HF and DFT study

GIAO/HF and DFT methods were utilized to predict the ¹³C chemical shifts of substituted ketenimines. GIAO HF/6–311+G(2d,p) and B3LYP/6–311+G(2d,p) methods were applied on the optimized B3LYP/6–31G(d) geometries and ¹³C chemical shifts of C_α and C_β of substituted ketenimines were correlated with group electronegativities. HF and DFT calculations indicated that increasing substituent group electronegativity leads to increasing chemical shift of C_β of substituted ketenimines, whereas the C_α values decrease. Copyright © 2003 John Wiley & Sons, Ltd.     

Quantum Chemistry Studies on the Fe-Cu Interactions and 31p NMR in Fe（CO）3（Ph2Ppy）2（CuXn） （Xn = Cl2^2-, Cl-, Br-）

In order to study the Fe-Cu interactions and their effects on 31p NMR, the structures of mononuclear complex Fe（CO）3fPhzPpy）a 1 and binuclear complexes Fe（CO）3（PhEPpy）z（CuXn） （2： Xn = Cl2^2-, 3： Xn = Cl-, 4： Xn = Br-） are calculated by density functional theory （DFT） PBE0 method. For complexes 1, 3 and 4, the 31p NMR chemical shifts calculated by PBE0-GIAO method are in good agreement with experimental results. The 31p chemical shift is 82.10 ppm in the designed complex 2. The Fe-Cu interactions （including Fe→Cu and Fe←Cu charge transfer） mainly exhibit the indirect interactions. Moreover, the Fe-Cu（I） interactions （mostly acting as σFe-p→4Scu and aFe-C→4Scu charge transfer） in complexes 3 and 4 are stronger than Fe-Cu（Ⅱ） interactions （mostly acting as σFe-p→4Scu and σFe-p←4Sc,） in complex 2. In complex 2, the stronger Fe←Cu interac- tions, acting as σFe-p←44SCu charge transfer, increase the electron density on P nucleus, which causes the upfield 31p chemical shift compared with mononuclear complex 1. For 3 and 4, although a little deshielding for P nucleus is derived from the delocalization of σFe-p→4Scu due to the Fe→Cu interactions, the stronger σFe-c→np charge-transfer finally increases the electron density on P nucleus. As a result, an upfield 31p chemical shift is observed compared with 1. The stability follows the order of 2〉3=4, indicating that Fe（CO）3（PhzPpy）2（CuCl2） is stable and could be synthesized experimentally. The N-Cu（Ⅱ） interaction plays an important role in the stability of 2. Because the delocalization of σFe-p→4SCu and σFe-c→πc-o weakens the a bonds of Fe-C and ~r bonds of CO, it is favorable for increasing the catalytic activity of binuclear complexes. Complexes 3 and 4 are expected to show higher catalytic activity compared to 2.     

多氯二苯硫醚结构参数和热力学性质的密度泛函理论研究

在B3LYP/6-31G*水平上对209个多氯二苯硫醚(PCDS)系列化合物进行了全优化计算, 得到各分子的结构参数和热力学参数. 研究了这些参数与氯原子在苯环各位置的取代数目和相互位置(N_PCS)之间的关系发现: 分子平均极化率(α)、焓(H $)、自由能(G $)、恒容热容( )和熵(S $)与N_PCS之间有很强的相关性, 相关性系数r²＞0.988, 分子体积(V_m)、最高占据轨道能(E_HOMO)和最低未占据轨道能(E_LUMO)与N_PCS也有较好的相关性, 相关性系数r²分别为0.949, 0.894和0.915. 设计等键反应, 计算了PCDS各异构体的标准生成热(Δ_fH $)和标准生成自由能(Δ_fG $). 根据异构体自由能的相对大小, 求得异构体的相对稳定性顺序.     

A GIAO/DFT study of 1H, 13C and 15N shieldings in amines and its relevance in conformational analysis

The 1H, 13C and 15N absolute shieldings of 13 amines were calculated at the GIAO/B3LYP/6-311++G** level. For some compounds (ethylamine, piperidine and 1-methylpiperidine) two conformations were calculated. The 13C and 15N data could be correctly correlated with experimental chemical shifts, allowing the conformation of 1-methylpiperidine to be established. The 1H NMR absolute shieldings, although less well correlated with delta values, were used to account for the anisotropy effects of the N lone pair.     

An NMR, IR and theoretical investigation of (1)H chemical shifts and hydrogen bonding in phenols

The change in (1)H NMR chemical shifts upon hydrogen bonding was investigated using both experimental and theoretical methods. The (1)H NMR spectra of a number of phenols were recorded in CDCl(3) and DMSO solvents. For phenol, 2- and 4-cyanophenol and 2-nitrophenol the OH chemical shifts were measured as a function of concentration in CDCl(3). The plots were all linear with concentration, the gradients varying from 0.940 (phenol) to 7.85 (4-cyanophenol) ppm/M because of competing inter- and intramolecular hydrogen bonding. Ab initio calculations of a model acetone/phenol system showed that the OH shielding was linear with the H...O=C distance (R) for R < 2.1 A with a shielding coefficient of - 7.8 ppm/A and proportional to cos(2)phi where phi is the H...O=C--C dihedral angle. Other geometrical parameters had little effect. It was also found that the nuclear shielding profile is unrelated to the hydrogen bonding energy profile. The dependence of the OH chemical shift on the pi density on the oxygen atom was determined as ca 40 ppm/pi electron. This factor is similar to that for NH but four times the value for sp(2) hybridized carbon atoms. The introduction of these effects into the CHARGE programme allowed the calculation of the (1)H chemical shifts of the compounds studied. The CHARGE calculations were compared with those from the ACD database and from GIAO calculations. The CHARGE calculations were more accurate than other calculations both when all the shifts were considered and also when the OH shifts were excluded. The calculations from the ACD and GIAO approaches were reasonable when the OH shifts were excluded but not as good when all the shifts were considered. The poor treatment of the OH shifts in the GIAO calculations is very likely due to the lack of explicit solvent effects in these calculations.