Structure of irradiated polytetrafluoroethylene according to DFT calculations of NMR chemical shifts

Methods of density functional theory were used to calculate ¹H, ¹³C, and ¹⁹F magnetic shielding tensors for C₇F_nH_mO_l model molecules, which can arise as fragments from radiation exposure of polytetrafluoroethylene.     

Enthalpies of formation from B3LYP calculations

We have calculated the geometries, energies, and normal vibrations of 845 compounds containing the elements H, C, N, O, F, Al, Si, P, S, and Cl using hybrid density functional theory in order to investigate the accuracy of atom-additive schemes for predicting enthalpies of formation at 298 K. The results give a more realistic estimate of the accuracy of density functional calculations than some overoptimistic earlier correlations. We have also calculated atom-additive schemes for the zero-point energies and enthalpic corrections to the energies. Remarkably, it is not important to include the vibrational or rotational contributions, which can be estimated well within a purely Born-Oppenheimer regression model.     

Ab initio calculations of NMR chemical shifts

The nuclear magnetic resonance chemical shift is one of the most powerful properties available for structure determination at the molecular level. A review of advances made in the ab initio calculation of chemical shielding during the past five years is presented. Specifically, progress in the areas including the effects of an unpaired electron, electron correlation, and relativistic effects into ab initio chemical shielding calculations, the tensor nature of the chemical shift, and intramolecular and intermolecular effects on the chemical shift will be covered.     

Quantum-chemical calculations of NMR chemical shifts of organic molecules: VI. Accuracy of DFT calculations of 29Si chemical shifts of four-coordinate silicon compounds

Systematic analysis of factors affecting the accuracy of DFT calculations of ²⁹Si NMR chemical shifts in four-coordinate silicon compounds showed that the best agreement with the experimental values is attained using B1PW91 and PBE0 functionals in combination with the TZP basis set. In calculations of ²⁹Si chemical shifts by quantum-chemical methods particular attention should be given to the contribution of relativistic spin-orbit interaction and conformational equilibrium.     

DFT calculations of 1H and 13C NMR chemical shifts in transition metal hydrides

Transition metal hydrides are of great interest in chemistry because of their reactivity and their potential use as catalysts for hydrogenation. Among other available techniques, structural properties in transition metal (TM) complexes are often probed by NMR spectroscopy. In this paper we will show that it is possible to establish a viable methodological strategy in the context of density functional theory, that allows the determination of 1H NMR chemical shifts of hydride ligands attached to transition metal atoms in mononuclear systems and clusters with good accuracy with respect to experiment. 13C chemical shifts have also been considered in some cases. We have studied mononuclear ruthenium complexes such as Ru(L)(H)(dppm)2 with L = H or Cl, cationic complex [Ru(H)(H2O)(dppm)2]+ and Ru(H)2(dppm)(PPh3)2, in which hydride ligands are characterized by a negative 1H NMR chemical shift. For these complexes all calculations are in relatively good agreement compared to experimental data with errors not exceeding 20% except for the hydrogen atom in Ru(H)2(dppm)(PPh3)2. For this last complex, the relative error increases to 30%, probably owing to the necessity to take into account dynamical effects of phenyl groups. Carbonyl ligands are often encountered in coordination chemistry. Specific issues arise when calculating 1H or 13C NMR chemical shifts in TM carbonyl complexes. Indeed, while errors of 10 to 20% with respect to experiment are often considered good in the framework of density functional theory, this difference in the case of mononuclear carbonyl complexes culminates to 80%: results obtained with all-electron calculations are overall in very satisfactory agreement with experiment, the error in this case does not exceed 11% contrary to effective core potentials (ECPs) calculations which yield errors always larger than 20%. We conclude that for carbonyl groups the use of ECPs is not recommended, although their use could save time for very large systems, for instance in cluster chemistry. The reliance of NMR chemical shielding on dynamical effects, such as intramolecular rearrangements or trigonal twists, is also examined for H2Fe(CO)4, K+[HFe(CO)](-), HMn(CO)5 and HRe(CO)5. The accuracy of the theory is also examined for complexes with two dihydrogen ligands (Tp*RuH(H2)2 and [FeH(H2)(DMPE)2]+) and a ruthenium cluster, [H3Ru4(C6H6)4(CO)]+. It is shown that for all complexes studied in this work, the effect of the ligands on the chemical shielding of hydrogen coordinated to metal is suitably calculated, thus yielding a very good correlation between experimental chemical shifts and theoretical chemical shielding.     

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.     

Quantum-chemical calculations of NMR chemical shifts of organic molecules: XII. Calculation of the 13C NMR chemical shifts of fluoromethanes at the DFT level

Calculation was carried out of chemical shifts in ¹³C NMR spectra for a series of fluoromethanes CH _n F_4?n (n = 0–4) by the methods of the electron density functional theory GIAO-DFT taking in consideration the solvent effect in the framework of the polarizable continuum model Tomasi IEF-PCM. The best results were obtained at the use of Keal-Tozer KT3 functional combined with Pople standard basis sets 6-311G(d,p) and 6-311++G(d,p), and also with Jensen special set pcS-2 containing tight p-functions. The optimum reference in the calculation of chemical shifts in ¹³C NMR spectra for the fluoromethanes series is TMS.     

Complete basis set B3LYP NMR calculations of CDCl3 solvent's water fine spectral details

The assignment of singlet at 1.55 ppm and the 1:1:1 triplet at 1.519 ppm to H(2)O and HOD in the 400 MHz (1)H NMR spectrum of CDCl(3) solvent were supported by complete basis set (CBS) GIAO-B3LYP calculated chemical shift and the CBS B3LYP estimated (2)J(D,H) spin-spin coupling constant (SSCC). The CBS fitting of B3LYP/cc-pCVxZ and B3LYP/pcJ-n predicted SSCC values, the accurate value of (2)J(D,H) = -1.082 +/- 0.030 Hz of HOD in chloroform-d(1) and the H/D isotopic shift of 0.0307(1) ppm were reported for the first time. The agreement between CBS B3LYP predicted chemical shift, spin-spin values and experiment was good.     

DFT-GIAO calculations of 19F NMR chemical shifts for perfluoro compounds

The (19)F NMR shieldings for 53 kinds of perfluoro compounds were calculated by the B3LYP-GIAO method using the 6-31G(d), 6-31+G(d), 6-31G(d,p), 6-31++G(d,p), 6-311G(d,p), 6-311++G(d,p), 6-311G(2d,2p), 6-311++G(2d,2p), 6-311++G(2df,2p), 6-311++G(3d,2p), and 6-311++G(3df,2p) basis sets. The diffuse functions markedly reduce the difference between the calculated and experimental chemical shifts. The calculations using the 6-31++G(d,p) basis set give the chemical shifts within 10 ppm deviations from experimental values except for the fluorine nuclei attached to an oxygen atom, a four- and a six-coordinated sulfur atom, and FC(CF(3))(2) attached to a sulfur atom.     

DFT(B3LYP) calculations of silatranes and their radical cations. First ionization potentials

DFT(B3LYP) and 2 quantum chemical calculations have been performed for 1-substituted silatranes XSi(OCH₂CH₂)N (X = H, CH₃, CH₂Cl, F), their radical cations, and first ionization potentials (IP₁) of these silatranes. The calculated values of IP₁ agree well with the experiment and make it possible to assign the first band to IP₁ in the photoelectron spectra. Analysis of spin density distribution and electronic charges in the radical cations suggests that ionization occurs mainly due to the lone electron pair of nitrogen, participating in intramolecular coordination. The N → Si interaction is broken, and the N...Si distance increases to 335–340 pm.     

1H NMR spectra of alcohols in hydrogen bonding solvents: DFT/GIAO calculations of chemical shifts

Proton nuclear magnetic resonance (NMR) shifts of aliphatic alcohols in hydrogen bonding solvents have been computed on the basis of density functional theory by applying the gauge‐including atomic orbital method to geometry‐optimized alcohol/solvent complexes. The OH proton shifts and hydrogen bond distances for methanol or ethanol complexed with pyridine depend very much on the functional employed and very little on the basis set, provided it is sufficiently large to give the correct quasi‐linear hydrogen bond geometry. The CH proton shifts are insensitive to both the functional and the basis set. NMR shifts for all protons in several alcohol/pyridine complexes are calculated at the Perdew, Burke and Ernzerhof PBE0/cc‐pVTZ//PBE0/6‐311 + G(d,p) level in the gas phase. The results correlate with the shifts for the pyridine‐complexed alcohols, determined by analysing data from the NMR titration of alcohols against pyridine. More pragmatically, computed shifts for a wider range of alcohols correlate with experimental shifts in neat pyridine. Shifts for alcohols in dimethylsulfoxide, based on the corresponding complexes in the gas phase, correlate well with the experimental values, but the overall root mean square difference is high (0.23 ppm), shifts for the OH, CH OH and other CH protons being systematically overestimated, by averages of 0.42, 0.21 and 0.06 ppm, respectively. If the computed shifts are corrected accordingly, a very good correlation is obtained with a gradient of 1.00 ± 0.01, an intercept of 0.00 ± 0.02 ppm and a root mean square difference of 0.09 ppm. This is a modest improvement on the result of applying the CHARGE programme to a slightly different set of alcohols. Some alcohol complexes with acetone and acetonitrile were investigated both in the gas phase and in a continuum of the relevant solvent. Copyright © 2015 John Wiley & Sons, Ltd.     

Relativistic four-component DFT calculations of 1H NMR chemical shifts in transition-metal hydride complexes: unusual high-field shifts beyond the Buckingham-Stephens model

State-of-the-art relativistic four-component DFT-GIAO-based calculations of (1)H NMR chemical shifts of a series of 3d, 4d, and 5d transition-metal hydrides have revealed significant spin-orbit-induced heavy atom effects on the hydride shifts, in particular for several 4d and 5d complexes. The spin-orbit (SO) effects provide substantial, in some cases even the dominant, contributions to the well-known characteristic high-field hydride shifts of complexes with a partially filled d-shell, and thereby augment the Buckingham-Stephens model of off-center paramagnetic ring currents. In contrast, complexes with a 4d(10) and 5d(10) configuration exhibit large deshielding SO effects on their hydride (1)H NMR shifts. The differences between the two classes of complexes are attributed to the dominance of π-type d-orbitals for the true transition-metal systems compared to σ-type orbitals for the d(10) systems.     

Density functional calculations of NMR chemical shifts and ESR g-tensors

An overview is given on recent advances of density functional theory (DFT) as applied to the calculation of nuclear magnetic resonance (NMR) chemical shifts and electron spin resonance (ESR) g-tensors. This is a new research area that has seen tremendous progress and success recently; we try to present some of these developments. DFT accounts for correlation effects efficiently. Therefore, it is the only first-principle method that can handle NMR calculations on large systems like transition-metal complexes. Relativistic effects become important for heavier element compounds; here we show how they can be accounted for. The ESR g-tensor is related conceptually to the NMR shielding, and results of g-tensor calculations are presented. DFT has been very successful in its application to magnetic properties, for metal complexes in particular. However, there are still certain shortcomings and limitations, e.g., in the exchange-correlation functional, that are discussed as well. Received: 24 October 1997 / Accepted: 19 December 1997     

B3LYP calculations on bishomoaromaticity in substituted semibullvalenes*

B3LYP/6‐31G* calculations on the degenerate rearrangements of substituted semibullvalenes spuriously predict the relative enthalpies of the bishomoaromatic TSs to be lower than the experimental values. However, the calculations do make the useful and experimentally testable prediction that the two cyano and two phenyl substituents in 2,6‐dicyano‐4,8‐diphenylsemibullvalene ( 9d ) are more likely than four cyano substituents in 2,4,6,8‐tetracyanosemibullvalene ( 9f ) or the four phenyl substituents in 2,4,6,8‐tetraphenylsemibullvalene ( 9g ) to produce a semibullvalene that has a bishomoaromatic equilibrium geometry in the gas phase. The major reason for the surprising finding that 9d is more likely to be bishomoaromatic than 9g is shown to be steric interactions between the phenyl groups at C‐2 and C‐8 and at C‐4 and C‐6 in bishomoaromatic structure 10g . These interactions inhibit the conjugative stabilization of 10g ; but they are absent in bishomoaromatic structure 10d , where cyano groups replace the phenyl groups at C‐2 and C‐6 in 10g . © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1565–1573, 2001     

MINDO/3 calculations of ESCA chemical shifts

Calculations of ESCA chemical shifts, using Jolly's equivalent core approximation and the MINDO/3 semi-empirical SCF MO method, have given results in reasonable agreement with experiment.     

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.     

Density functional calculations of the 13C NMR chemical shifts in (9,0) single-walled carbon nanotubes

The electronic structure and (13)C NMR chemical shift of (9,0) single-walled carbon nanotubes (SWNTs) are investigated theoretically. Shielding tensor components are also reported. Density functional calculations were carried out for C(30)-capped and H-capped fragments which serve as model systems for the infinite (9,0) SWNT. Based on the vanishing HOMO-LUMO gap, H-capped nanotube fragments are predicted to exhibit "metallic" behavior. The (13)C chemical shift approaches a value of approximately 133 ppm for the longest fragment studied here. The C(30)-capped SWNT fragments of D(3d)/D(3h) symmetry, on the other hand, are predicted to be small-gap semiconductors just like the infinite (9,0) SWNT. The differences in successive HOMO-LUMO gaps and HOMO and LUMO energies, as well as the (13)C NMR chemical shifts, converge slightly faster with the fragment's length than for the H-capped tubes. The difference between the H-capped and C(30)-capped fragments is analyzed in some detail. The results indicate that (at least at lengths currently accessible to quantum chemical computations) the H-capped systems represent less suitable models for the (9,0) SWNT because of pronounced artifacts due to their finite length. From our calculations for the C(30)-capped fragments, the chemical shift of a carbon atom in the (9,0) SWNT is predicted to be about 130 ppm. This value is in reasonably good agreement with experimental estimates for the (13)C chemical shift in SWNTs.     

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.     

Quantum-chemical calculations of NMR chemical shifts of organic molecules: XIV. Solvation effects in calculations of chemical shifts in 13C NMR spectra of chlorine-containing compounds

Analysis of precision factors in calculations of ¹³C NMR chemical shifts in the series of saturated and unsaturated organochlorine compounds was performed in the framework of the method of electron density functional theory GIAO-DFT-KT3/pcS-2 in the gas phase and with accounting for solvent effect by the polarized continuum model IEF-PCM. The accounting for solvation effects in calculations of ¹³C NMR chemical shifts within the framework of the IEF-PCM model is not fundamental for organochlorine compounds, yet it considerably improves the precision of calculations up to 2.5 ppm.     

Application of quantum chemical calculations of C NMR chemical shifts to quinoxaline structure determination

Comparison of experimental and theoretical (GIAO DFT) ¹³C NMR chemical shifts allows the reliable assignment of isomeric structures of heteroaromatic compounds. This methodology was applied to establish the structures of isomeric quinoxalines. A modern 1D NOE technique permitted independent proof of the proposed structures.