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.     

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.     

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.     

Non‐empirical calculations of NMR indirect carbon‐carbon coupling constants. Part 5: Bridged bicycloalkanes

All possible J(C,C) of the bicarbocyclic frameworks together with J(C,H) and J(H,H) at bridgeheads in the series of six bridged bicycloalkanes, bicyclo[1.1.0]butane, bicyclo[2.1.0]pentane, bicyclo[3.1.0]hexane, bicyclo[2.2.0]hexane, bicyclo[3.2.0]heptane and bicyclo[3.3.0]octane, were calculated at the SOPPA level with correlation consistent Dunning sets cc‐pVTZ‐Cs augmented with inner core s‐functions and locally dense Sauer sets aug‐cc‐pVTZ‐J augmented with tight s‐functions and rationalized in terms of the multipath coupling mechanism and hybridization effects explaining many interesting structural trends. Copyright © 2003 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.     

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.     

Non‐empirical calculations of NMR indirect carbon–carbon coupling constants: 1. Three‐membered rings

The carbon–carbon indirect nuclear spin–spin coupling constants in cyclopropane, aziridine and oxirane were investigated by means of ab initio calculations at the RPA, SOPPA and DFT/B3LYP levels. We found that the carbon–carbon couplings are by far dominated by the Fermi contact term. Our best SOPPA and DFT results are in a very good agreement with each other and with the experimental values, whereas calculations at the RPA level of theory strongly overestimate the carbon–carbon couplings. Significant differences in the basis set dependence of the calculated carbon–carbon coupling constants obtained with either wavefunction method, RPA or SOPPA, or the density functional method, DFT/B3LYP, are observed. The SOPPA results depend much more strongly on the quality of the basis set than the results of DFT/B3LYP calculations. The medium‐sized core‐valence basis sets cc‐pCVTZ and even cc‐pCVDZ were found to perform fairly well at the SOPPA level for the one‐bond carbon–carbon couplings investigated here. Copyright © 2002 John Wiley & Sons, Ltd.     

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.     

Difluorobenzenes revisited: an experimental and theoretical study of spin–spin coupling constants for 1,2‐, 1,3‐, and 1,4‐difluorobenzene

The experimental spin–spin coupling constants (SSCCs) for 1,3‐ and 1,4‐difluorobenzene have been determined anew, and found to be consistent with previously determined values. SSCCs for 1,2‐, 1,3‐, and 1,4‐difluorobenzene have been analyzed by comparing them with the coupling constants computed using the second‐order polarization propagator approximation (SOPPA) and the equation‐of‐motion coupled cluster singles and doubles method (EOM‐CCSD). Eighty experimental values have been analyzed using SOPPA calculations, and a subset of 40 values using both SOPPA and EOM‐CCSD approaches. One‐bond coupling constants ¹J(C? C) and ¹J(C? F) are better described by EOM‐CCSD, whereas one‐bond ¹J(C? H) values are better described by SOPPA. An empirical equation is presented which allows for the prediction of unknown coupling constants from computed SOPPA values. A similar approach may prove useful for predicting coupling constants in larger systems. Copyright © 2009 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.     

Open‐chain unsaturated selanyl sulfides: stereochemical structure and stereochemical behavior of their 77Se–1H spin–spin coupling constants

Stereochemical structure of nine Z‐2‐(vinylsulfanyl)ethenylselanyl organyl sulfides has been investigated by means of experimental measurements and second‐order polarization propagator approach calculations of their ¹H–¹H, ¹³C–¹H, and ⁷⁷Se–¹H spin–spin coupling constants together with a theoretical conformational analysis performed at the MP2/6‐311G** level. All nine compounds were shown to adopt the preferable skewed s‐cis conformation of their terminal vinylsulfanyl group, whereas the favorable rotational conformations with respect to the internal rotations around the C–S and C–Se bonds of the internal ethenyl group are both skewed s‐trans. Stereochemical trends of ⁷⁷Se–¹H spin–spin coupling constants originating in the geometry of their coupling pathways and the selenium lone pair effect were rationalized in terms of the natural J‐coupling analysis within the framework of the natural bond orbital approach. Copyright © 2012 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.     

Can One Assess the π Character of a C–C Bond with the Help of the NMR Spin–Spin Coupling Constants?

Measured one‐bond spin–spin coupling constants (SSCC) ¹J(CC) can be used to describe the nature of the C–C bond, provided one is able to separate the various coupling mechanisms leading to ¹J(CC). The Fermi‐contact (FC) term probes the first‐order density at the positions of the coupling nuclei, whereas the noncontact terms (the paramagnetic spin orbit (PSO) and the spin–dipole (SD) terms) probe the π character of the C–C bond (the diamagnetic spin orbit (DSO) term can mostly be neglected). A model is tested, in which the value of the FC(CC) term is estimated with the help of measured SSCCs ¹J(CH). The difference between the measured J(CC) and the estimated FC(CC) values, Δ(CC)=PSO(CC)+SD(CC)+DSO(CC), provides a semiquantitative measure of the π character of a C–C multiple bond. The applicability and limitations of this approach are discussed by partitioning the four Ramsey terms of the SSCC ¹J(CC) into one‐ and two‐orbital contributions. The FC, PSO, and SD terms of ¹J(CC) are explained and analyzed with regard to their relationship to other C–C bond properties. It is shown that empirical relationships between measured SSCCs and the s character of a bond need reconsideration.     

Conformation analysis of d‐glucaric acid in deuterium oxide by NMR based on its JHH and JCH coupling constants

d ‐Glucaric acid (GA) is an aldaric acid and consists of an asymmetric acyclic sugar backbone with a carboxyl group positioned at either end of its structure (i.e., the C1 and C6 positions). The purpose of this study was to conduct a conformation analysis of flexible GA as a solution in deuterium oxide by NMR spectroscopy, based on J‐resolved conformation analysis using proton–proton (³J_HH) and proton–carbon (²J_CH and ³J_CH) coupling constants, as well as nuclear overhauser effect spectroscopy (NOESY). The ²J_CH and ³J_CH coupling constants were measured using the J‐resolved heteronuclear multiple bond correlation (HMBC) NMR technique. NOESY correlation experiments indicated that H2 and H5 were in close proximity, despite the fact that these protons were separated by too large distance in the fully extended form of the chain structure to provide a NOESY correlation. The validities of the three possible conformers along the three different bonds (i.e., C2? C3, C3? C4, and C4? C5) were evaluated sequentially based on the J‐coupling values and the NOESY correlations. The results of these analyses suggested that there were three dominant conformers of GA, including conformer 1 , which was H2H3:gauche, H3H4:anti, and H4H5:gauche; conformer 2 , which was H2H3:gauche, H3H4:anti, and H4H5:anti; and conformer 3 , which was H2H3:gauche, H3H4: gauche, and H4H5:anti. These results also suggested that all three of these conformers exist in equilibrium with each other. Lastly, the results of the current study suggested that the conformational structures of GA in solution were ‘bent’ rather than being fully extended. Copyright © 2016 John Wiley & Sons, Ltd.     

Conformational analysis and stereochemical dependences of 31P–1H spin–spin coupling constants of bis(2‐phenethyl)vinylphosphine and related phosphine chalcogenides

Theoretical energy‐based conformational analysis of bis(2‐phenethyl)vinylphosphine and related phosphine oxide, sulfide and selenide synthesized from available secondary phosphine chalcogenides and vinyl sulfoxides is performed at the MP2/6‐311G** level to study stereochemical behavior of their ³¹P–¹H spin–spin coupling constants measured experimentally and calculated at different levels of theory. All four title compounds are shown to exist in the equilibrium mixture of two conformers: major planar s‐cis and minor orthogonal ones, while ³¹P–¹ H spin–spin coupling constants under study are found to demonstrate marked stereochemical dependences with respect to the geometry of the coupling pathways, and to the internal rotation of the vinyl group around the P(X)‐C bonds (X = LP, O, S and Se), opening a new guide in the conformational studies of unsaturated phosphines and phosphine chalcogenides. Copyright © 2009 John Wiley & Sons, Ltd.     

Precise Measurement of Long‐Range Heteronuclear Coupling Constants by a Novel Broadband Proton–Proton‐Decoupled CPMG‐HSQMBC Method

A broadband proton–proton‐decoupled CPMG‐HSQMBC method for the precise and direct measurement of long‐range heteronuclear coupling constants is presented. The Zangger–Sterk‐based homodecoupling scheme reported herein efficiently removes unwanted proton–proton splittings from the heteronuclear multiplets, so that the desired heteronuclear couplings can be determined simply by measuring frequency differences between singlet maxima in the resulting spectra. The proposed pseudo‐1D/2D pulse sequences were tested on nucleotides, a metal complex incorporating P heterocycles, and diglycosyl (di)selenides, as well as on other carbohydrate derivatives, for the extraction of ⁿJ(¹H,³¹P), ⁿJ(¹H,⁷⁷Se), and ⁿJ(¹H,¹³C) values, respectively.     

Conformational Dependence of Through‐Space Tellurium–Tellurium Spin–Spin Coupling in Peri‐Substituted Bis(Tellurides)

Three related series of peri‐substituted bis(tellurides) bearing naphthalene, acenaphthene and acenaphthylene backbones (Nap/Acenap/Aceyl(TeY)₂ (Nap=naphthalene‐1,8‐diyl N ; Acenap=acenaphthene‐5,6‐diyl A ; Aceyl=acenaphthylene‐5,6‐diyl Ay ; Y=Ph 1 ; Fp 2 ; Tol 3 ; An‐p­ 4 ; An‐o­ 5 ; Tp 6 ; Mes 7 ; Tip 8 ) have been synthesised and their solid‐state structures determined by X‐ray crystallography. Molecular conformations were classified as a function of the two C9‐C‐Te‐C(Y) dihedral angles (θ); in the solid all members adopt AB or CCt configurations, with larger Te(aryl) moieties exclusively imposing the CCt variant. Exceptionally large J(¹²⁵Te,¹²⁵Te) spin–spin coupling constants between 3289–3848 Hz were obtained for compounds substituted by bulky Te(aryl) groups, implying these species are locked in a CCt‐type conformation. In contrast, compounds incorporating smaller Te(aryl) moieties are predicted to be rather dynamic in solution and afford much smaller J values (2050–2676 Hz), characteristic of greater populations of AB conformers with lower couplings. This conformational dependence of through‐space coupling is supported by DFT calculations.     

Experimental and theoretical study of the intramolecular C–H···N and C–H···S hydrogen bonding effects in the 1H and 13C NMR spectra of the 2‐(alkylsulfanyl)‐5‐amino‐1‐vinylpyrroles: a particular state of amine nitrogen

In the ¹H NMR spectra of the 1‐vinylpyrroles with amino‐ and alkylsulfanyl groups in 5 and 2 positions, an extraordinarily large difference between resonance positions of the H_A and H_B terminal methylene protons of the vinyl group is discovered. Also, the one‐bond ¹J(C_β,H_B) coupling constant is surprisingly greater than the ¹J(C_β,H_A) coupling constant in pyrroles under investigation, while in all known cases, there was a reverse relationship between these coupling constants. These spectral anomalies are substantiated by quantum chemical calculations. The calculations show that the amine nitrogen lone pair is removed from the conjugation with the π‐system of the pyrrole ring so that it is directed toward the H_B hydrogen. These factors are favorable to the emergence of the intramolecular C–H_B???N hydrogen bonding in the s‐cis(N) conformation. On the other hand, the spatial proximity of the sulfur to the H_B hydrogen provides an opportunity of the intramolecular C–H_B???S hydrogen bonding in the s‐cis(S) conformation. Presence of the hydrogen bond critical points as well as ring critical point for corresponding chelate ring revealed by a quantum theory of atoms in molecules (QTAIM) approach confirms the existence of the weak intramolecular C–H???N and C–H???S hydrogen bonding. Therefore, an unusual high‐frequency shift of the H_B signal and the increase in the ¹J(C_β,H_B) coupling constant can be explained by the effects of hydrogen bonding. Copyright © 2013 John Wiley & Sons, Ltd.     

A theoretical structural analysis of the factors that affect 1JNH, 1hJNH and 2hJNN in N–H···N hydrogen‐bonded complexes

Calculations of ¹ J_NH, ^1h J_NH and ^2h J_NN spin–spin coupling constants of 27 complexes presenting N–H·N hydrogen bonds have allowed to analyze these through hydrogen‐bond coupling as a function of the hybridization of both nitrogen atoms and the charge (+1, 0, ? 1) of the complex. The main conclusions are that the hybridization of N atom of the hydrogen bond donor is much more important than that of the hydrogen bond acceptor. Positive and negative charges (cationic and anionic complexes) exert opposite effects while the effect of the transition states ‘proton‐in‐the‐middle’ is considerable. Copyright © 2008 John Wiley & Sons, Ltd.     

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.