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.     

Compatibility of correlation‐consistent basis sets with a hybrid Hartree–Fock/density functional method

The present study examines the feasibility of combining the correlation‐consistent basis sets developed by Dunning and coworkers with the hybrid Hartree–Fock/density functional method B3LYP. Furthermore, extrapolation to the complete basis set (CBS) limit minimizes errors due to the presence of an incomplete basis set and can act as a rigorous test of the limitations of the B3LYP method. Equilibrium geometries, energies, and harmonic vibrational frequencies were determined for a series of well‐studied, yet computationally challenging, small inorganics and their respective ions. The results were then extrapolated to the CBS limit, where applicable, and compared to experiment. It was found that a union between the hybrid Hartree–Fock/density functional B3LYP method and Dunning's augmented correlation‐consistent basis sets gave results that were comparable to molecular orbital methods that explicitly account for electron correlation. Furthermore, the minimum basis set necessary to attain reasonable results for the systems studied was aug‐cc‐pVTZ. Upgrading to the aug‐cc‐pVQZ level and subsequent extrapolation to the CBS limit further improved the overall agreement with the experiment. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 207–216, 1999     

Non-empirical calculations of NMR indirect carbon-carbon coupling constants. Part 7--spiroalkanes

Carbon-carbon spin-spin coupling constants were calculated at the SOPPA level for a series of seven classical spiroalkanes, spiro[2.2]pentane, spiro[2.3]hexane, spiro[2.4]heptane, spiro[2.5]octane, spiro[3.3]heptane, spiro[4.4]nonane and spiro[5.5]undecane, with special focus upon couplings involving and/or across spiro carbons. Many interesting structural trends were investigated originating in specific geometries and unusual bonding environments at the spiro carbon.     

Non-empirical calculations of NMR indirect carbon-carbon coupling constants. Part 8--monocycloalkanes

Carbon-carbon and carbon-hydrogen spin-spin coupling constants were calculated in the series of the first six monocycloalkanes using SOPPA and SOPPA(CCSD) methods, and very good agreement with the available experimental data was achieved, with the latter method showing slightly better results in most cases, at least in those involving calculations of J(C,C). Benchmark calculations of all possible 21 coupling constants J(C,C), J(C,H) and J(H,H) in chair cyclohexane revealed the importance of using the appropriate level of theory and adequate quality of the basis sets. Many unknown couplings in this series were predicted with high confidence and several interesting structural trends (hybridization effects, multipath coupling transmission mechanisms, hyperconjugative interactions) were elucidated and are discussed based on the present calculations of spin-spin couplings.     

Non-empirical calculations of NMR indirect spin-spin coupling constants. Part 13: configurational assignment of aminosulfonylamidines

Carbon-carbon coupling constants have been experimentally measured using the INADEQUATE pulse sequence in the series of N-[1,2-bis(dialkylamino)-2-arylethylidene]arylsulfonamides obtained from N-(1-aryl-2,2,2-trichloroethyl)amides of arylsulfonic acids. Comparison of the experimental J(C,C) in this series with those calculated at the SOPPA (Second-Order Polarization Propagator Approach) level in the model aminosulfonylamidine provided an unambiguous assignment at the C=N bond of the eight-title aminosulfonylamidines to E-configuration, while the unknown J(C,C) couplings in their inaccessible diverse Z isomers have been predicted with high reliability. The established marked difference between J(C,C) of the corresponding carbon-carbon bonds in cis and trans orientations to the nitrogen lone pair in aminosulfonylamidines provides a powerful tool in the configurational assignment at the C=N bond in a wide series of the related systems containing the C=N-SO2R moiety.     

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.     

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.     

Non-empirical calculations of NMR indirect carbon-carbon coupling constants. Part 10--carbocages

A comprehensive theoretical study of nine classical caged polycycloalkanes (tetrahedrane, prismane, homoprismane, quadricyclane, cubane, pentaprismane, hexaprismane, adamantane and diamantane) was carried out with special focus on the structural behavior of their J(C,C) values calculated at the SOPPA level. The structural behavior of J(C,C) in small carbocages is dominated by steric strain whereas in medium-sized polycycloalkanes the J(C,C) values show no marked peculiarities and follow several well-defined structural trends typical of other alicyclic compounds.     

Nonempirical calculations of NMR indirect spin-spin coupling constants. Part 15: pyrrolylpyridines

Conformational study of 2-(2-pyrrolyl)pyridine and 2,6-di(2-pyrrolyl)pyridine was performed on the basis of the experimental measurements and high-level ab initio calculations of the one-bond 13C-13C, 13C-1H and 15N-1H spin-spin coupling constants showing marked stereochemical behavior upon the internal rotation around the pyrrole-pyridine interheterocyclic bonds. Both compounds were established to adopt predominant s-cis conformations with no noticeable out-of-plane deviations.     

Non-empirical calculations of NMR indirect carbon-carbon coupling constants. Part 6: propellanes

A full set of carbon-carbon coupling constants have been calculated at the SOPPA level in the series of six most representative propellanes. Special attention was focused on spin-spin couplings involving both bridgehead carbons, and these data were rationalized in terms of the multipath coupling mechanism and hybridization effects. Many unknown couplings in the propellane frameworks were predicted with high reliability.     

Nonempirical calculations of NMR indirect spin-spin coupling constants. Part 14: azomethines of the alpha, beta-unsaturated aldehydes

Configurational assignment of seven azomethines obtained from the alpha-functionally substituted and nonsubstituted alpha,beta-unsaturated aldehydes has been performed on the basis of experimental measurements and the high-level ab initio calculations of their 1J(C,C) and 1J(C,H), involving the alpha-imino carbon that demonstrated the marked stereochemical dependence of both coupling constants upon the orientation of the nitrogen lone pair in the diverse isomers of the title azomethines.     

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.     

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.     

Nuclear magnetic resonance J coupling constant polarizabilities of hydrogen peroxide: A basis set and correlation study

In this article, we present the so far most extended investigation of the calculation of the coupling constant polarizability of a molecule. The components of the coupling constant polarizability are derivatives of the nuclear magnetic resonance (NMR) indirect nuclear spin–spin coupling constant with respect to an external electric field and play an important role for both chiral discrimination and solvation effects on NMR coupling constants. In this study, we illustrate the effects of one‐electron basis sets and electron correlation both at the level of density functional theory as well as second‐order polarization propagator approximation for the small molecule hydrogen peroxide, which allowed us to perform calculations with the largest available basis sets optimized for the calculation of NMR coupling constants. We find a systematic but rather slow convergence with the one‐electron basis set and that augmentation functions are required. We observe also large and nonsystematic correlation effects with significant differences between the density functional and wave function theory methods. © 2012 Wiley Periodicals, Inc.     

Calculation of the vibration frequencies of alpha-quartz: the effect of Hamiltonian and basis set

The central-zone vibrational spectrum of alpha-quartz (SiO2) is calculated by building the Hessian matrix numerically from the analytical gradients of the energy with respect to the atomic coordinates. The nonanalytical part is obtained with a finite field supercell approach for the high-frequency dielectric constant and a Wannier function scheme for the evaluation of Born charges. The results obtained with four different Hamiltonians, namely Hartree-Fock, DFT in its local (LDA) and nonlocal gradient corrected (PBE) approximation, and hybrid B3LYP, are discussed, showing that B3LYP performs far better than LDA and PBE, which in turn provide better results than HF, as the mean absolute difference from experimental frequencies is 6, 18, 21, and 44 cm(-1), respectively, when a split valence basis set containing two sets of polarization functions is used. For the LDA results, comparison is possible with previous calculations based on the Density Functional Perturbation Theory and usage of a plane-wave basis set. The effects associated with the use of basis sets of increasing size are also investigated. It turns out that a split valence plus a single set of d polarization functions provides frequencies that differ from the ones obtained with a double set of d functions and a set of f functions on all atoms by on average less than 5 cm(-1).     

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.     

Complete basis set limit of Ab initio binding energies and geometrical parameters for various typical types of complexes

Using basis‐set extrapolation schemes for a given data set, we evaluated the binding energies and geometries at the complete basis set (CBS) limit at the levels of the second order Møller–Plesset perturbation theory (MP2) and the coupled cluster theory with singles, doubles, and perturbative triples excitations [CCSD(T)]. The systems include the hydrogen bonding (water dimer), aromatic interaction (benzene dimer), π–H interaction (benzene–water), cation–water, anion–water, π–cation interaction (cation–benzene), and π–anion interaction (anion–triazine). One extrapolation method is to exploit both BSSE‐corrected and BSSE‐uncorrected binding energies for the aug‐cc‐pVNZ (N = 2, 3, 4, …) basis set in consideration that both binding energies give the same CBS limit (CBS^B). Another CBS limit (CBS^C) is to use the commonly known extrapolation approach to exploit that the electron correlation energy is proportional to N^?3. Since both methods are complementary, they are useful for estimating the errors and trend of the asymptotic values. There is no significant difference between both methods. Overall, the values of CBS^C are found to be robust because of their consistency. However, for small N (in particular, for N = 2, 3), CBS is found to be slightly better for water–water interactions and cation–water and cation–benzene interactions, whereas CBS is found to be more reliable for bezene–water and anion–water interactions. We also note that the MP2 CBS limit value based on N = 2 and 3 combined with the difference between CCSD(T) and MP2 at N = 2 would be exploited to obtain a CCSD(T)/CBS value for aromatic–aromatic interactions and anion–π interactions, but not for cationic complexes. © 2007 Wiley Periodicals, Inc. J Comput Chem, 2008     

B3LYP density functional study on spectroscopic properties of CuO− and CuS−

Spectroscopic constant and molecular properties of CuO^? and CuS^? in their ground state were studied in detail using hybrid the HF/DF B3LYP density functional method. Four basis sets of Cu, one of them fully uncontracted, were used for this calculation. The basis sets for O and S are correlation consistent basis sets. The basis sets were developed systematically in order to see the effect of basis sets on the spectroscopic properties and also to check the consistency of the calculated properties. The B3LYP spectroscopic properties of CuO^? agree very well with the only available experimental values. Most of the spectroscopic properties of CuO^? and particularly those of CuS^? are reported for the first time. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005