Benchmark calculations of 29Si–1H spin–spin coupling constants across double bond

Benchmark calculations of geminal and vicinal ²⁹Si–¹H spin–spin coupling constants across double bond in three reference alkenylsilanes have been carried out at both DFT and SOPPA levels in comparison with experiment. At the former, four density functionals, B3LYP, B3PW91, PBE0 and KT3, were tested in combination with five representative basis sets. At the latter, three main SOPPA‐based methods, SOPPA, SOPPA(CC2) and SOPPA(CCSD), were examined in combination with the same series of basis sets. On the whole, the wavefunction methods showed much better results as compared to DFT, with the most efficient combination of SOPPA/cc‐pVTZ‐su2 characterized by a mean absolute error of only 0.4 Hz calculated for a set of nine coupling constants in three compounds with a sample span of around 40 Hz. Copyright © 2012 John Wiley & Sons, Ltd.     

Prediction of water's isotropic nuclear shieldings and indirect nuclear spin–spin coupling constants (SSCCs) using correlation‐consistent and polarization‐consistent basis sets in the Kohn–Sham basis set limit

Density functional theory (DFT) was used to estimate water's isotropic nuclear shieldings and indirect nuclear spin–spin coupling constants (SSCCs) in the Kohn–Sham (KS) complete basis set (CBS) limit. Correlation‐consistent cc‐pVxZ and cc‐pCVxZ (x = D, T, Q, 5, and 6), and their modified versions (ccJ‐pVxZ, unc‐ccJ‐pVxZ, and aug‐cc‐pVTZ‐J) and polarization‐consistent pc‐n and pcJ‐n (n = 0, 1, 2, 3, and 4) basis sets were used, and the results fitted with a simple mathematical formula. The performance of over 20 studied density functionals was assessed from comparison with the experiment. The agreement between the CBS DFT‐predicted isotropic shieldings, spin–spin values, and the experimental values was good and similar for the modified correlation‐consistent and polarization‐consistent basis sets. The BHandH method predicted the most accurate ¹H, ¹⁷O isotropic shieldings and ¹J(OH) coupling constant (deviations from experiment of about ? 0.2 and ? 1 ppm and 0.6 Hz, respectively). The performance of BHandH for predicting water isotropic shieldings and ¹J(OH) is similar to the more advanced methods, second‐order polarization propagator approximation (SOPPA) and SOPPA(CCSD), in the basis set limit. Copyright © 2008 John Wiley & Sons, Ltd.     

Structural trends of 77Se?1H spin–spin coupling constants and conformational behavior of 2‐substituted selenophenes

Experimental measurements and second‐order polarization propagator approach (SOPPA) calculations of ⁷⁷Se? ¹H spin–spin coupling constants together with theoretical energy‐based conformational analysis in the series of 2‐substituted selenophenes have been carried out. A new basis set optimized for the calculation of ⁷⁷Se? ¹H spin–spin coupling constants has been introduced by extending the aug‐cc‐pVTZ‐J basis for selenium. Most of the spin–spin coupling constants under study, especially vicinal ⁷⁷Se? ¹H couplings, demonstrated a remarkable stereochemical behavior with respect to the internal rotation of the substituent in the 2‐position of the selenophene ring, which is of major importance in the stereochemical studies of the related organoselenium compounds. Copyright © 2009 John Wiley & Sons, Ltd.     

Nonempirical calculations of the one‐bond 29Si–13C spin–spin coupling constants taking into account relativistic and solvent corrections

The computational study of the one‐bond ²⁹Si–¹³C spin–spin coupling constants has been performed at the second‐order polarization propagator approximation (SOPPA) level in the series of 60 diverse silanes with a special focus on the main factors affecting the accuracy of the calculation including the level of theory, the quality of the basis set, and the contribution of solvent and relativistic effects. Among three SOPPA‐based methods, SOPPA(MP2), SOPPA(CC2), and SOPPA(CCSD), the best result was achieved with SOPPA(CCSD) when used in combination with Sauer's basis set aug‐cc‐pVTZ‐J characterized by the mean absolute error of calculated coupling constants against the experiment of ca 2 Hz in the range of ca 200 Hz. The SOPPA(CCSD)/aug‐cc‐pVTZ‐J method is recommended as the most accurate and effective computational scheme for the calculation of ¹J(Si,C). The slightly less accurate but essentially more economical SOPPA(MP2)/aug‐cc‐pVTZ‐J and/or SOPPA(CC2)/aug‐cc‐pVTZ‐J methods are recommended for larger molecular systems. It was shown that solvent and relativistic corrections do not play a major role in the computation of the total values of ¹J(Si,C); however, taking them into account noticeably improves agreement with the experiment. The rovibrational corrections are estimated to be of about 1 Hz or 1–1.5% of the total value of ¹J(Si,C). Copyright © 2014 John Wiley & Sons, Ltd.     

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.     

Non‐empirical calculations of NMR indirect carbon–carbon coupling constants: 2. Strained polycarbocycles

High‐level non‐empirical calculations of carbon–carbon spin–spin coupling constants in a series of strained polycarbocycles have been carried out, in excellent agreement with available experimental data. The utmost importance of electronic correlation effects in this case has been demonstrated and it has been shown that the Second‐Order Polarization Propagator Approach (SOPPA) is an adequate method to account for those effects. It has been demonstrated that the most reliable basis sets to calculate J(C,C) at the SOPPA level are the correlation‐consistent basis sets of Dunning and co‐workers augmented with inner core s‐functions or decontracted in their s‐parts. The nature of the unusual bridgehead–bridgehead bonds in bicyclobutane and propellane in terms of s‐characters of bonding hybrids and also the hybridization effects in spiropentane are discussed based on the arguments derived from the current calculations of J(C,C) in the title compounds. The values of the unknown J(C,C) in propellane and spiropentane are predicted with high reliability. Copyright © 2003 John Wiley & Sons, Ltd.     

Magnitude of Finite‐Nucleus‐Size Effects in Relativistic Density Functional Computations of Indirect NMR Nuclear Spin–Spin Coupling Constants

A spherical Gaussian nuclear charge distribution model has been implemented for spin‐free (scalar) and two‐component (spin–orbit) relativistic density functional calculations of indirect NMR nuclear spin–spin coupling (J‐coupling) constants. The finite nuclear volume effects on the hyperfine integrals are quite pronounced and as a consequence they noticeably alter coupling constants involving heavy NMR nuclei such as W, Pt, Hg, Tl, and Pb. Typically, the isotropic J‐couplings are reduced in magnitude by about 10 to 15 % for couplings between one of the heaviest NMR nuclei and a light atomic ligand, and even more so for couplings between two heavy atoms. For a subset of the systems studied, viz. the Hg atom, Hg₂²⁺, and Tl? X where X=Br, I, the basis set convergence of the hyperfine integrals and the coupling constants was monitored. For the Hg atom, numerical and basis set calculations of the electron density and the 1s and 6s orbital hyperfine integrals are directly compared. The coupling anisotropies of TlBr and TlI increase by about 2 % due to finite‐nucleus effects.     

Computational NMR coupling constants: Shifting and scaling factors for evaluating 1JCH

Optimized shifting and/or scaling factors for calculating one‐bond carbon–hydrogen spin–spin coupling constants have been determined for 35 combinations of representative functionals (PBE, B3LYP, B3P86, B97‐2 and M06‐L) and basis sets (TZVP, HIII‐su3, EPR‐III, aug‐cc‐pVTZ‐J, ccJ‐pVDZ, ccJ‐pVTZ, ccJ‐pVQZ, pcJ‐2 and pcJ‐3) using 68 organic molecular systems with 88 ¹J_CH couplings including different types of hybridized carbon atoms. Density functional theory assessment for the determination of ¹J_CH coupling constants is examined, comparing the computed and experimental values. The use of shifting constants for obtaining the calculated coupling improves substantially the results, and most models become qualitatively similar. Thus, for the whole set of couplings and for all approaches excluding those using the M06 functional, the root‐mean‐square deviations lie between 4.7 and 16.4 Hz and are reduced to 4–6.5 Hz when shifting constants are considered. Alternatively, when a specific rovibrational contribution of 5 Hz is subtracted from the experimental values, good results are obtained with PBE, B3P86 and B97‐2 functionals in combination with HIII‐su3, aug‐cc‐pVTZ‐J and pcJ‐2 basis sets. Copyright © 2013 John Wiley & Sons, Ltd.     

Non‐empirical calculations of NMR indirect carbon–carbon coupling constants. Part 4: Bicycloalkanes

A systematic study of the one‐bond and long‐range J(C,C), J(C,H) and J(H,H) in the series of nine bicycloalkanes was performed at the SOPPA level with special emphasis on the coupling transmission mechanisms at bridgeheads. Many unknown couplings were predicted with high reliability. Further refinement of SOPPA computational scheme adjusted for better performance was carried out using bicyclo[1.1.1]pentane as a benchmark to investigate the influence of geometry, basis set and electronic correlation. The calculations performed demonstrated that classical ab initio SOPPA applied with the locally dense Dunning's sets augmented with inner core s‐functions used for coupled carbons and Sauer's sets augmented with tight s‐functions used for coupled hydrogens performs perfectly well in reproducing experimental values of different types of coupling constants (the estimated reliability is ca 1–2 Hz) in relatively large organic molecules of up to 11 carbon atoms. Additive coupling increments were derived for J(C,C), J(C,H) and J(H,H) based on the calculated values of coupling constants within SOPPA in the model bicycloalkanes, in reasonably good agreement with the known values obtained earlier on pure empirical grounds. Most of the bridgehead couplings in all but one bicycloalkane appeared to be essentially additive within ca 2–3 Hz while bicyclo[1.1.1]pentane demonstrated dramatic non‐additivity of ?14.5 Hz for J(C,C), +16.6 Hz for J(H,H) and ?5.5 Hz for J(C,H), in line with previous findings. Non‐additivity effects in the latter compound established at the SOPPA level should be attributed to the through‐space non‐bonded interactions at bridgeheads due to the essential overlapping of the bridgehead rear lobes which provides an additional and effective non‐bonding coupling path for the bridgehead carbons and their protons in the bicyclopentane framework. Copyright © 2003 John Wiley & Sons, Ltd.     

Fluorine spin–spin coupling constants of pentafluorobenzene revisited at the ab initio correlated levels

All possible spin–spin coupling constants, ¹⁹F–¹⁹F, ¹⁹F–¹³C, and ¹⁹F–¹H, of pentafluorobenzene were calculated at five different levels of theory, HF, DFT, SOPPA (CCSD), CCSD, and the SOPPA (CCSD)-based composite scheme with taking into account solvent, vibrational, relativistic, and correlation corrections. Most corrections were next to negligible for the long-range couplings but quite essential for the one-bond carbon–fluorine coupling constants. Hartree–Fock calculations were found to be entirely unreliable, while DFT results were comparable in accuracy with the data obtained using the wave function-based methods.     

Complete basis set prediction of methanol isotropic nuclear magnetic shieldings and indirect nuclear spin–spin coupling constants (SSCC) using polarization‐consistent and XZP basis sets and B3LYP and BHandH density functionals

Efficient B3LYP and BHandH density functionals were used to estimate methanol's nuclear magnetic isotropic shieldings and spin–spin coupling constants in the basis set limit. Polarization‐consistent pcS‐n and pcJ‐n (n = 0, 1, 2, 3 and 4), and segmented contracted XZP, where X = D, T, Q and 5, basis sets were used and the results fitted with simple mathematical formulas. The performance of the methods was assessed from comparison with experiment and higher level calculations. ¹J(CH) and ³J(HH) values were determined from very diluted solutions in deuterochloroform and compared with theoretical predictions. The agreement between complete basis set (CBS) density functional theory (DFT) predicted isotropic shieldings and spin–spin values and experiment was good. The BHandH/pcS‐n methanol shieldings obtained using structures optimized at the same level of theory are approaching the accuracy of the advanced coupled‐cluster‐singles‐doubles‐approximate triples (CCSD(T)) calculations. Copyright © 2009 John Wiley & Sons, Ltd.     

Non‐empirical calculations of NMR indirect carbon–carbon coupling constants: 3. Polyhedranes

High‐level ab initio calculations of carbon–carbon coupling constants were carried out in tetrahedrane, prismane and cubane using the SOPPA (Second‐Order Polarization Propagator Approach) computational scheme, in good agreement with available experimental data. It was found that SOPPA performs perfectly well in combination with Dunning's correlation‐consistent basis sets augmented with inner core functions; however, no improvement was observed on adding tight s‐functions. The utmost importance of electronic correlation effects decreasing the total values of computed J(C,C) in the title compounds by a factor of ～2.0–2.5 was found. Unknown values of J(C,C) in the title polyhedranes were predicted with high reliability and the latter were treated in terms of s‐characters of carbon–carbon bonds based on the additive scheme of coupling pathways. All three compounds under study showed decreased s‐characters of their carbon–carbon bonds, which is the result of their remarkable steric strain. Copyright © 2003 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.     

NMR spin–spin coupling constants in hydrogen‐bonded glycine clusters

The influence of the hydrogen bond formation on the NMR spin–spin coupling constants (SSCC), including the Fermi contact (FC), the diamagnetic spin‐orbit, the paramagnetic spin‐orbit, and the spin dipole term, has been investigated systematically for the homogeneous glycine cluster, in gas phase, containing up to three monomers. The one‐bond and two‐bond SSCCs for several intramolecular (through covalent bond) and intermolecular (across the hydrogen‐bond) atomic pairs are calculated employing the density functional theory with B3LYP and KT3 functionals and different types of extended basis sets. The ab initio SOPPA(CCSD) is used as benchmark for the SSCCs of the glycine monomer. The hydrogen bonding is found to cause significant variations in the one‐bond SSCCs, mostly due to contribution from electronic interactions. However, the nature of variation depends on the type of oxygen atom (proton‐acceptor or proton‐donor) present in the interaction. Two‐bond intermolecular coupling constants vary more than the corresponding one‐bond constants when the size of the cluster increases. Among the four Ramsey terms that constitute the total SSCC, the FC term is the most dominant contributor followed by the paramagnetic spin‐orbit term in all one‐bond interaction.     

Structural studies on aryl‐substituted enaminoketones and their thio analogues. Part II. Experimental and DFT calculated carbon–carbon coupling constants,nJ(CC)'s (n = 1–3)

Spin–spin carbon–carbon coupling constants across one, two and three bonds, J(CC), have been measured for a series of aryl‐substituted Z‐s‐Z‐s‐E enaminoketones and their thio analogues. As a result, a large set, altogether 178, of J(CC)s has been obtained. It consists of 82 couplings across one bond, 31 couplings across two bonds and 65 couplings across three bonds. Independently, the DFT calculations at the B3PW91/6‐311++G(d,p)//B3PW91/6‐311++G(d,p) level yielded a set of theoretical J(CC) values. A comparison of these two sets of data gave an excellent linear correlation with parameters a and b close to ideal; a = 0.9978 which is not far from unity and b = 0.22 Hz which is close to zero. The ¹J(CC) couplings determined for the crucial fragment of the molecules, i.e. ? C?C? C?O (or ? C?C? C?S), are: ¹J(C?C) ≈ 68 Hz (67 Hz) and ¹J(C? C) = 60.5 Hz (60.0 Hz). The corresponding couplings found for the Z‐s‐Z‐s‐E isomer of the parent enaminoketone, 4‐methylamino‐but‐3‐en‐2‐one are 64.1 and 59.3 Hz, respectively. The most sensitive towards substitution of the oxygen atom by sulfur are two‐bond couplings between the α‐vinylic and aromatic C_ipso carbon atoms, which attain 12 Hz in the enaminoketone derivatives and decrease to 5 Hz in their thio analogues. Copyright © 2009 John Wiley & Sons, Ltd.     

One‐bond 13C–13C coupling constants in alkyl‐substituted cyclopropenes: experimental and theoretical studies

Measurements of one‐bond carbon–carbon coupling constants, ¹J(C, C), were performed for two series of compounds, alkyl‐substituted cyclopropenes and cyclopropanes. The experimental data were complemented by a set of DFT‐calculated J couplings for the parent cyclopropene ( 1 ), its methyl and silyl derivatives and, additionally, for 1‐methylcyclobutene ( 3 ), 1‐methylcyclopentene ( 4 ) and 1‐methylcyclohexene ( 5 ) and good agreement was observed between the experimental and the calculated data; all the trends are perfectly maintained, including a dramatic decrease in the couplings across endocyclic single bonds in cyclopropene and its derivatives, and a significant decrease in the corresponding couplings in cyclobutene. Using the data obtained, the s characters of the carbon hybrid orbitals involved in the formation of the cyclopropene were calculated. The results clearly show that the ring closure and the related strain exerted upon the cyclopropene molecule only slightly disturb the electron structure of the double bond. The s character of the corresponding carbon orbital is 0.314 in cyclopropene vs the theoretical value of 0.333 in ethene. This is at variance with the endo‐ and exocyclic single bonds, where the s characters of the orbitals forming the endocyclic single bonds are much smaller than those of the bonds in the open‐chain compounds, i.e. 0.229 (C‐1 and/or C‐2) and 0.166 (C‐3). The s values calculated for the exocyclic CH bonds are 0.334 for C‐3 and 0.456 for C‐1 and/or C‐2. Copyright © 2002 John Wiley & Sons, Ltd.     

Indirect relativistic bridge and substituent effects from the ‘heavy’ environment on the one‐bond and two‐bond 13C1H spin–spin coupling constants

Indirect relativistic bridge effect (IRBE) and indirect relativistic substituent effect (IRSE) induced by the ‘heavy’ environment of the IV‐th, V‐th and VI‐th main group elements on the one‐bond and geminal ¹³C? ¹H spin–spin coupling constants are observed, and spin‐orbit parts of these two effects were interpreted in terms of the third‐order Rayleigh–Schrödinger perturbation theory. Both effects, IRBE and IRSE, rapidly increase with the total atomic charge of the substituents at the coupled carbon. The accumulation of IRSE for geminal coupling constants is not linear with respect to the number of substituents in contrast to the one‐bond couplings where IRSE is an essentially additive quantity. Copyright © 2015 John Wiley & Sons, Ltd.     

Relativistic effects in the one‐bond spin–spin coupling constants involving selenium

One‐bond spin–spin coupling constants involving selenium of seven different types, ¹ J(Se,X), X = ¹H, ¹³C, ¹⁵ N, ¹⁹ F, ²⁹Si, ³¹P, and ⁷⁷Se, were calculated in the series of 14 representative compounds at the SOPPA(CCSD) level taking into account relativistic corrections evaluated both at the RPA and DFT levels of theory in comparison with experiment. Relativistic corrections were found to play a major role in the calculation of ¹ J(Se,X) reaching as much as almost 170% of the total value of ¹ J(Se,Se) and up to 60–70% for the rest of ¹ J(Se,X). Scalar relativistic effects (Darwin and mass‐velocity corrections) by far dominate over spin–orbit coupling in the total relativistic effects for all ¹ J(Se,X). Taking into account relativistic corrections at both random phase approximation and density functional theory levels essentially improves the agreement of theoretical results with experiment. The most ‘relativistic’ ¹ J(Se,Se) demonstrates a marked Karplus‐type dihedral angle dependence with respect to the mutual orientation of the selenium lone pairs providing a powerful tool for stereochemical analysis of selenoorganic compounds. Copyright © 2014 John Wiley & Sons, Ltd.     

Basis set error estimation for DFT calculations of electronic g‐tensors for transition metal complexes

We present a detailed study of the basis set dependence of electronic g‐tensors for transition metal complexes calculated using Kohn–Sham density functional theory. Focus is on the use of locally dense basis set schemes where the metal is treated using either the same or a more flexible basis set than used for the ligand sphere. The performance of all basis set schemes is compared to the extrapolated complete basis set limit results. Furthermore, we test the performance of the aug‐cc‐pVTZ‐J basis set developed for calculations of NMR spin‐spin and electron paramagnetic resonance hyperfine coupling constants. Our results show that reasonable results can be obtain when using small basis sets for the ligand sphere, and very accurate results are obtained when an aug‐cc‐pVTZ basis set or similar is used for all atoms in the complex. © 2014 Wiley Periodicals, Inc.     

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.