The use of locally dense basis sets in the calculation of indirect nuclear spin–spin coupling constants: The vicinal coupling constants in H3C–CH2X (X=H, F, Cl, Br, I)

On the geometry dependence of the nuclear magnetic resonance chemical shift of mercury in thiolate complexes: A relativistic density functional theory study

Thiolate containing mercury(II) complexes of the general formula [Hg(SR)n ]2 - n have been of great interest since the toxicity of mercury was recognized. 199Hg nuclear magnetic resonance spectroscopy (NMR) is a powerful tool for characterization of mercury complexes. In this work, the Hg shielding constants in a series of [Hg(SR)n ]2 - n complexes are therefore investigated computationally with particular emphasis on their geometry dependence. Geometry optimizations and NMR chemical shift calculations are performed at the density functional theory (DFT) level with both the zeroth-order regular approximation (ZORA) and four-component relativistic methods. The four exchange-correlation (XC) functionals PBE0, PBE, B3LYP, and BLYP are used in combination with either Dyall's Gaussian-type (GTO) or Slater-type orbitals (STOs) basis sets. Comparing ZORA and four-component calculations, one observes that the calculated shielding constants for a given molecular geometry have a constant difference of ∼ 1070 ppm. This confirms that ZORA is an acceptable relativistic method to compute NMR chemical shifts. The combinations of four-component/PBE0/v3z and ZORA/PBE0/QZ4P are applied to explore the geometry dependence of the isotropic shielding. For a given coordination number, the distance between mercury and sulfur is the key factor affecting the shielding constant, while changes in bond and dihedral angles and even different side groups have relatively little impact.

Exploring Alternate Methods for the Calculation of High-Level Vibrational Corrections of NMR Spin-Spin Coupling Constants

Traditional nuclear magnetic resonance (NMR) calculations typically treat systems with a Born-Oppenheimer-derived electronic wave function that is solved for a fixed nuclear geometry. One can numerically account for this neglected nuclear motion by averaging over property values for all nuclear geometries with a vibrational wave function and adding this expectation value as a correction to an equilibrium geometry property value. Presented are benchmark coupled-cluster singles and doubles (CCSD) vibrational corrections to spin-spin coupling constants (SSCCs) computed at the level of vibrational second-order perturbation theory (VPT2) using the vibrational averaging driver of the CFOUR program. As CCSD calculations of vibrational corrections are very costly, cheaper electronic structure methods are explored via a newly developed Python vibrational averaging program within the Dalton Project. Namely, results obtained with the second-order polarization propagator approximation (SOPPA) and density functional theory (DFT) with the B3LYP and PBE0 exchange-correlation functionals are compared to the benchmark CCSD//CCSD(T) and experimental values. CCSD//CCSD(T) corrections are also combined with literature CC3 equilibrium geometry values to form the highest-order vibrationally corrected values available (i.e., CC3//CCSD(T) + CCSD//CCSD(T)). CCSD//CCSD(T) statistics showed favorable statistics in comparison to experimental values, albeit at an unfavorably high computational cost. A cheaper CCSD//CCSD(T) + B3LYP method showed quite similar mean absolute deviation (MAD) values as CCSD//CCSD(T), concluding that CCSD//CCSD(T) + B3LYP is optimal in terms of cost and accuracy. With reference to experimental values, a vibrational correction was not worth the cost for all of the other methods tested. Finally, deviation statistics showed that CC3//CCSD(T) + CCSD//CCSD(T) vibrational-corrected equilibrium values deteriorated in comparison to CCSD//CCSD(T) attributed to the use of a smaller basis set or lack of solvation effects for the CC3 equilibrium calculations.

Computational 19 F NMR of trifluoromethyl derivatives of alkenes, pyrimidines, and indenes

The ¹⁹ F NMR chemical shifts of 13 trifluoromethyl derivatives of alkenes, pyrimidines, and indenes were calculated at the DFT level using the BhandHLYP, BHandH, PBE, PBE0, O3LYP, B3LYP, KT2, and KT3 functionals in combination with the pcS-2 basis set. Best result was documented for the BHandHLYP functional: The mean absolute error (MAE) of 0.66 ppm for the scaled values was achieved for the range of about 20 ppm. Solvent, vibrational, and relativistic corrections were found to be rather small, especially when taken in combination, generally demonstrating a slight decrease in the difference between calculated and experimental fluorine chemical shifts. As a measure of the practical importance of these compounds, one should recall that the growing number of life science products that contain trifluoromethyl groups provides a continuing driving force for the development of an effective methodology that enables both regio- and stereoselective introduction of trifluoromethyl groups into both aliphatic and aromatic systems.

Electron Correlation or Basis Set Quality: How to Obtain Converged and Accurate NMR Shieldings for the Third-Row Elements?

The quality of theoretical NMR shieldings calculated at the quantum-chemical level depends on various theoretical aspects, of which the basis set type and size are among the most important factors. Nevertheless, not much information is available on the basis set effect on theoretical shieldings of the NMR-active nuclei of the third row. Here, we report on the importance of proper basis set selection to obtain accurate and reliable NMR shielding parameters for nuclei from the third row of the periodic table. All calculations were performed on a set of eleven compounds containing the elements Na, Mg, Al, Si, P, S, or Cl. NMR shielding tensors were calculated using the SCF-HF, DFT-B3LYP, and CCSD(T) methods, combined with the Dunning valence aug-cc-pVXZ, core-valence aug-cc-pCVXZ, Jensen polarized-convergent aug-pcSseg-n and Karlsruhe x2c-Def2 basis set families. We also estimated the complete basis set limit (CBS) values of the NMR parameters. Widely scattered nuclear shieldings were observed for the Dunning polarized-valence basis set, which provides irregular convergence. We show that the use of Dunning core-valence or Jensen basis sets effectively reduces the scatter of theoretical NMR results and leads to their exponential-like convergence to CBS. We also assessed the effect of vibrational, temperature, and relativistic corrections on the predicted shieldings. For systems with single bonds, all corrections are relatively small, amounting to less than 4% of the CCSD(T)/CBS value. Vibrational and temperature corrections were less reliable for H₃PO and HSiCH due to the high anharmonicity of the molecules. An abnormally high relativistic correction was observed for phosphorus in PN, reaching ~20% of the CCSD(T)/CBS value, while the correction was less than 7% for other tested molecules.

Computational 1 H and 13 C NMR in structural and stereochemical studies

Present review outlines the advances and perspectives of computational ¹ H and ¹³ C NMR applied to the stereochemical studies of inorganic, organic, and bioorganic compounds, involving in particular natural products, carbohydrates, and carbonium ions. The first part of the review briefly outlines theoretical background of the modern computational methods applied to the calculation of chemical shifts and spin-spin coupling constants at the DFT and the non-empirical levels. The second part of the review deals with the achievements of the computational ¹ H and ¹³ C NMR in the stereochemical investigation of a variety of inorganic, organic, and bioorganic compounds, providing in an abridged form the material partly discussed by the author in a series of parent reviews. Major attention is focused herewith on the publications of the recent years, which were not reviewed elsewhere.

Quantum Chemical Approaches to the Calculation of NMR Parameters: From Fundamentals to Recent Advances

Quantum chemical methods for the calculation of indirect NMR spin–spin coupling constants and chemical shifts are always in progress. They never stay the same due to permanently developing computational facilities, which open new perspectives and create new challenges every now and then. This review starts from the fundamentals of the nonrelativistic and relativistic theory of nuclear magnetic resonance parameters, and gradually moves towards the discussion of the most popular common and newly developed methodologies for quantum chemical modeling of NMR spectra.

An efficient method for generating property-energy consistent basis sets. New pecJ-n (n = 1, 2) basis sets for high-quality calculations of indirect nuclear spin-spin coupling constants involving 1H, 13C, 15N, and 19F nuclei

This paper presents a new method of generating property-energy consistent (PEC) basis sets that can be applied to any arbitrary molecular property. The PEC method generates a basis set that is optimized for the molecular property under interest, providing the least possible total molecular energy. The main algorithm of the PEC approach involves Monte Carlo simulations to generate random exponents in the predetermined range. In this work, the PEC method is introduced in the example of generation of new pecJ-n (n = 1, 2) basis sets suited for high-quality correlated calculations of indirect nuclear spin-spin coupling constants involving the most popular NMR-active nuclei: 1H, 13C, 15N, and 19F.

Computational NMR of Carbohydrates: Theoretical Background, Applications, and Perspectives

This review is written amid a marked progress in the calculation of NMR parameters of carbohydrates substantiated by a vast amount of experimental data coming from several laboratories worldwide. By no means are we trying to cover in the present compilation a huge amount of all available data. The main idea of the present review was only to outline general trends and perspectives in this dynamically developing area on the background of a marked progress in theoretical and computational NMR. Presented material is arranged in three basic sections: (1)-a brief theoretical introduction; (2)-applications and perspectives in computational NMR of monosaccharides; and (3)-calculation of NMR chemical shifts and spin-spin coupling constants of di- and polysaccharides.

Estimating the accuracy of calculated electron paramagnetic resonance hyperfine couplings for a lytic polysaccharide monooxygenase

Lytic polysaccharide monooxygenases (LPMOs) are enzymes that bind polysaccharides followed by an (oxidative) disruption of the polysaccharide surface, thereby boosting depolymerization. The binding process between the LPMO catalytic domain and polysaccharide is key to the mechanism and establishing structure-function relationships for this binding is therefore crucial. The hyperfine coupling constants (HFCs) from EPR spectroscopy have proven useful for this purpose. Unfortunately, EPR does not provide direct structural data and therefore the experimental EPR parameters have to be supported with parameters calculated with density functional theory. Yet, calculated HFCs are extremely sensitive to the employed computational setup. Using the LPMO Ls(AA9)A catalytic domain, we here quantify the importance of several choices in the computational setup, ranging from the use of specialized basis, the underlying structures, and the employed exchange-correlation functional. We show that specialized basis sets are an absolute necessity, and also that care has to be taken in the optimization of the underlying structure: only by allowing large parts of the protein around the active site to structurally relax could we obtain results that uniformly reproduced experimental trends. We compare our results to previously published X-ray structures and experimental HFCs for Ls(AA9)A as well as to recent experimental/theoretical results for another (AA10) family of LPMOs.

The 1 H and 13 C NMR chemical shifts of Strychnos alkaloids revisited at the DFT level

The density functional theory calculation of ¹ H and ¹³ C NMR chemical shifts in a series of ten 10 classically known Strychnos alkaloids with a strychnine skeleton was performed at the PBE0/pcSseg-2//pcseg-2 level. It was found that calculated ¹ H and ¹³ C NMR chemical shifts provided a markedly good correlation with experiment characterized by a mean absolute error of 0.08 ppm in the range of 7 ppm for protons and 1.67 ppm in the range of 150 ppm for carbons, so that a mean absolute percentage error was as small as ~1% in both cases.

Computational 1 H NMR: Part 1. Theoretical background

This is the first one of the three closely interrelated reviews to be published in Magnetic Resonance in Chemistry dealing with accordingly theoretical background, chemical applications, and biochemical studies of and by means of computational ¹ H NMR. Presented in the first part of the review is a general outline of the modern theoretical methods and accuracy factors of computational ¹ H NMR involving locally dense basis set schemes, solvent effects, vibrational corrections, and relativistic effects performed at the density functional theory and/or nonempirical levels. This review is dedicated to Prof. Stephan Sauer in view of his invaluable contribution to the field of computational nuclear magnetic resonance.

Computational protocols for calculating 13C NMR chemical shifts

The most recent results dealing with the computation of ¹³C NMR chemical shifts in chemistry (small molecules, saturated, unsaturated and aromatic compounds, heterocycles, functional derivatives, coordination complexes, carbocations, and natural products) are reviewed, paying special attention to theoretical background and accuracy, the latter involving solvent effects, vibrational corrections, and relativistic effects.

Exploring Through-Space Spin-Spin Couplings for Quantum Information Processing: Facing the Challenge of Coherence Time and Control Quantum States

Nuclear magnetic resonance (NMR) is a powerful tool for studying quantum information processing (QIP). Recently quantum technologies have been proposed to overcome the challenges in large-scale NMR QIP. Furthermore, computational chemistry can promote its improvement. Nuclear spins-¹/₂ are natural qubits and have been used in most NMR quantum computation experiments. However, molecules that enable many qubits NMR QIP implementations should meet some requirements regarding their spectroscopic properties. Exceptionally large through-space (TS) P-P spin-spin coupling constants (SSCC or J) observed in 1,8-diphosphanaphthalenes (PPN) and in naphtho[1,8- cd]-1,2-dithiole phenylphosphines (NTP) were proposed and investigated to provide more accurate control within large-scale NMR QIP. Spectroscopic properties of PPN and NTP derivatives were explored by theoretical strategies using locally dense basis sets (LDBS). ³¹P chemical shifts (δ) calculated at the B3LYP/aug-cc-pVTZ-J level and TS P-P SSCCs at the PBE1PBE/pcJ-2 (LDBS-1) level are very close to the experimental data for the PPN molecule. Differently, for the NTP dimer, PBE1PBE/pcJ-2 (LDBS-2) predicts more accurate ³¹P δ, whereas PBE1PBE/Def2-TZVP (LDBS-1) forecasts more accurate TS P-P SSCCs. From our results, PPN_o-F, PPN_o-ethyl, and PPN_o-NH₂ were the best candidates for NMR QIP, in which the large TS SSCCS could face the need of long-time quantum gates implementations. Therefore, it could overcome natural limitations concerning the development of large-scale NMR.

Relativistic heavy atom effect on the 31 P NMR parameters of phosphine chalcogenides. Part 1. Chemical shifts

Four-component density functional theory calculations of ³¹ P NMR chemical shifts have been performed for the representative series of 56 phosphine chalcogenides in order to investigate an influence of different functional groups on the heavy atom relativistic effect on the NMR chemical shifts of light phosphorous atoms (Heavy Atom on Light Atom effect). The validity of the 4-component density functional theory approach used for the wide-scale calculations of the phosphorous chemical shifts in a wide series of phosphine chalcogenides has been confirmed on a small series of 5 representative compounds with the aid of high-quality coupled cluster singles and doubles calculations taking into account solvent, vibrational, and the relativistic corrections in comparison with the experiment.

Theoretical calculations of carbon-hydrogen spin-spin coupling constants

Structural applications of theoretical calculations of carbon-hydrogen spin-spin coupling constants are reviewed covering papers published mainly during the last 10-15 years with a special emphasis on the most notable studies of hybridization, substitution and stereoelectronic effects together with the investigation of hydrogen bonding and intermolecular interactions. The wide scope of different applications of calculated carbon-hydrogen couplings in the structural elucidation of particular classes of organic and bioorganic molecules is reviewed, concentrating mainly on saturated, unsaturated, aromatic and heteroaromatic compounds and their functional derivatives, as well as on natural compounds and carbohydrates. The review is dedicated to Professor Emeritus Michael Barfield in view of his invaluable pioneering contribution to this field.

GIAO-DFT calculation of 15 N NMR chemical shifts of Schiff bases: Accuracy factors and protonation effects

¹⁵ N NMR chemical shifts in the representative series of Schiff bases together with their protonated forms have been calculated at the density functional theory level in comparison with available experiment. A number of functionals and basis sets have been tested in terms of a better agreement with experiment. Complimentary to gas phase results, 2 solvation models, namely, a classical Tomasi's polarizable continuum model (PCM) and that in combination with an explicit inclusion of one molecule of solvent into calculation space to form supermolecule 1:1 (SM + PCM), were examined. Best results are achieved with PCM and SM + PCM models resulting in mean absolute errors of calculated ¹⁵ N NMR chemical shifts in the whole series of neutral and protonated Schiff bases of accordingly 5.2 and 5.8 ppm as compared with 15.2 ppm in gas phase for the range of about 200 ppm. Noticeable protonation effects (exceeding 100 ppm) in protonated Schiff bases are rationalized in terms of a general natural bond orbital approach.

Substitution effects in the 15 N NMR chemical shifts of heterocyclic azines evaluated at the GIAO-DFT level

A systematic study of the accuracy factors for the computation of ¹⁵ N NMR chemical shifts in comparison with available experiment in the series of 72 diverse heterocyclic azines substituted with a classical series of substituents (CH₃ , F, Cl, Br, NH₂ , OCH₃ , SCH₃ , COCH₃ , CONH₂ , COOH, and CN) providing marked electronic σ- and π-electronic effects and strongly affecting ¹⁵ N NMR chemical shifts is performed. The best computational scheme for heterocyclic azines at the DFT level was found to be KT3/pcS-3//pc-2 (IEF-PCM). A vast amount of unknown ¹⁵ N NMR chemical shifts was predicted using the best computational protocol for substituted heterocyclic azines, especially for trizine, tetrazine, and pentazine where experimental ¹⁵ N NMR chemical shifts are almost totally unknown throughout the series. It was found that substitution effects in the classical series of substituents providing typical σ- and π-electronic effects followed the expected trends, as derived from the correlations of experimental and calculated ¹⁵ N NMR chemical shifts with Swain-Lupton's F and R constants.

Carbon-carbon spin-spin coupling constants: Practical applications of theoretical calculations

Practical applications of theoretical calculations of carbon-carbon spin-spin coupling constants in particular classes of organic and bioorganic molecules are reviewed, concentrating mainly on saturated, unsaturated, aromatic and heteroaromatic compounds and their functional derivatives as well as on carbohydrates and natural compounds.

On the accuracy factors and computational cost of the GIAO-DFT calculation of 15 N NMR chemical shifts of amides

The main factors affecting the accuracy and computational cost of Gauge-independent Atomic Orbital-density functional theory (GIAO-DFT) calculation of ¹⁵ N NMR chemical shifts in the benchmark series of 16 amides are considered. Among those are the choice of the DFT functional and basis set, solvent effects, internal reference conversion factor and applicability of the locally dense basis set (LDBS) scheme. Solvent effects are treated within the polarizable continuum model (PCM) scheme as well as at supermolecular level with solvent molecules considered in explicit way. The best result is found for Keal and Tozer's KT3 functional used in combination with Jensen's pcS-3 basis set with taking into account solvent effects within the polarizable continuum model. The proposed LDBS scheme implies pcS-3 on nitrogen and pc-2 elsewhere in the molecule. The resulting mean average error for the calculated ¹⁵ N NMR chemical shifts is about 6 ppm. The application of the LDBS approach tested in a series of 16 amides results in a dramatic decrease in computational cost (more than an order of magnitude in time scale) with insignificant loss of accuracy.

Calculation of 15N NMR chemical shifts: Recent advances and perspectives

Recent advances in computation of ¹⁵N NMR chemical shifts are reviewed, concentrating mainly on practical aspects of computational protocols and accuracy factors. The review includes the discussion of the level of theory, the choice of density functionals and basis sets together with taking into account solvent effects, rovibrational corrections and relativistic effects. Computational aspects of ¹⁵N NMR are illustrated for the series of neutral and protonated open-chain nitrogen-containing compounds and nitrogen heterocycles, coordination and intermolecular complexes.

A Systematic Comparison of Second-Order Polarization Propagator Approximation (SOPPA) and Equation-of-Motion Coupled Cluster Singles and Doubles (EOM-CCSD) Spin-Spin Coupling Constants for Selected Singly Bonded Molecules, and the Hydrides NH3, H2O, and HF and Their Protonated and Deprotonated Ions and Hydrogen-Bonded Complexes

Second-order polarization propagator approximation (SOPPA) and equation-of-motion coupled cluster singles and doubles (EOM-CCSD) methods have been employed for the calculation of one-bond spin-spin coupling constants in series of small molecules and ions, and of one- and two-bond coupling constants across X-H···Y hydrogen bonds. For isolated molecules, one-bond SOPPA coupling constants (1)J(X-Y) involving (13)C, (15)N, (17)O, and (19)F have larger absolute values than corresponding EOM-CCSD coupling constants, with the EOM-CCSD values being in significantly better agreement with available experimental data. The difference between SOPPA and EOM-CCSD tends to increase as the number of nonbonding electrons on the coupled atoms increases, and the SOPPA values for O-F coupling are significantly in error. Similarly, the absolute values of SOPPA one-bond coupling constants (1)J(X-H) for the hydrides NH3, H2O, and FH and their protonated and deprotonated ions are greater than EOM-CCSD values, with the largest differences occurring for F-H coupling. One- and two-bond coupling constants (1)J(X-H), (1h)J(H-Y), and (2h)J(X-Y) across X-H···Y hydrogen bonds in neutral, protonated, and deprotonated complexes formed from the hydrides are similar at SOPPA and EOM-CCSD, with the largest differences again found for (1)J(F-H) in complexes with F-H as the proton donor, and (2h)J(F-F) for (FHF)(-). The signs of (1)J(X-H), (1h)J(H-Y), and (2h)J(X-Y) are the same at both levels of theory, as is their variation across the proton-transfer coordinate in F-H···NH3. SOPPA would appear to provide a reliable and more cost-effective alternative approach for computing coupling constants across hydrogen bonds, although couplings involving F may be problematic.

Theoretical and experimental (15)N NMR study of enamine-imine tautomerism of 4-trifluoromethyl[b]benzo-1,4-diazepine system

The tautomeric structure of 4-trifluoromethyl[b]benzo-1,4-diazepine system in solution has been evaluated by means of the calculation of (15)N NMR chemical shifts of individual tautomers in comparison with the averaged experimental shifts to show that the enamine-imine equilibrium is entirely shifted toward the imine form. The adequacy of the theoretical level used for the computation of (15)N NMR chemical shifts in this case has been verified based on the benchmark calculations in the series of the push-pull and captodative enamines together with related azomethynes, which demonstrated a good to excellent agreement with experiment.

MP2 calculation of (77) Se NMR chemical shifts taking into account relativistic corrections

The main factors affecting the accuracy and computational cost of the Second-order Möller-Plesset perturbation theory (MP2) calculation of (77) Se NMR chemical shifts (methods and basis sets, relativistic corrections, and solvent effects) are addressed with a special emphasis on relativistic effects. For the latter, paramagnetic contribution (390-466 ppm) dominates over diamagnetic term (192-198 ppm) resulting in a total shielding relativistic correction of about 230-260 ppm (some 15% of the total values of selenium absolute shielding constants). Diamagnetic term is practically constant, while paramagnetic contribution spans over 70-80 ppm. In the (77) Se NMR chemical shifts scale, relativistic corrections are about 20-30 ppm (some 5% of the total values of selenium chemical shifts). Solvent effects evaluated within the polarizable continuum solvation model are of the same order of magnitude as relativistic corrections (about 5%). For the practical calculations of (77) Se NMR chemical shifts of the medium-sized organoselenium compounds, the most efficient computational protocols employing relativistic Dyall's basis sets and taking into account relativistic and solvent corrections are suggested. The best result is characterized by a mean absolute error of 17 ppm for the span of (77) Se NMR chemical shifts reaching 2500 ppm resulting in a mean absolute percentage error of 0.7%.

Towards the versatile DFT and MP2 computational schemes for 31P NMR chemical shifts taking into account relativistic corrections

The main factors affecting the accuracy and computational cost of the calculation of (31)P NMR chemical shifts in the representative series of organophosphorous compounds are examined at the density functional theory (DFT) and second-order Møller-Plesset perturbation theory (MP2) levels. At the DFT level, the best functionals for the calculation of (31)P NMR chemical shifts are those of Keal and Tozer, KT2 and KT3. Both at the DFT and MP2 levels, the most reliable basis sets are those of Jensen, pcS-2 or larger, and those of Pople, 6-311G(d,p) or larger. The reliable basis sets of Dunning's family are those of at least penta-zeta quality that precludes their practical consideration. An encouraging finding is that basically, the locally dense basis set approach resulting in a dramatic decrease in computational cost is justified in the calculation of (31)P NMR chemical shifts within the 1-2-ppm error. Relativistic corrections to (31)P NMR absolute shielding constants are of major importance reaching about 20-30 ppm (ca 7%) improving (not worsening!) the agreement of calculation with experiment. Further better agreement with the experiment by 1-2 ppm can be obtained by taking into account solvent effects within the integral equation formalism polarizable continuum model solvation scheme. We recommend the GIAO-DFT-KT2/pcS-3//pcS-2 scheme with relativistic corrections and solvent effects taken into account as the most versatile computational scheme for the calculation of (31)P NMR chemical shifts characterized by a mean absolute error of ca 9 ppm in the range of 550 ppm.

Relativistic effects in the one-bond spin-spin coupling constants involving selenium

One-bond spin-spin coupling constants involving selenium of seven different types, (1)  J(Se,X), X = (1) H, (13) C, (15)  N, (19)  F, (29) Si, (31) P, and (77) 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 (1)  J(Se,X) reaching as much as almost 170% of the total value of (1)  J(Se,Se) and up to 60-70% for the rest of (1)  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 (1)  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' (1)  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.

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

The computational study of the one-bond (29)Si-(13)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 (1)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 (1)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 (1)J(Si,C).

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 (77)Se-(13)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 (1)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 (1)J(Se,C). Solvent effects in the values of (1)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 (1)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 (77)Se-(13)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.

On the accuracy of the GIAO-DFT calculation of 15N NMR chemical shifts of the nitrogen-containing heterocycles--a gateway to better agreement with experiment at lower computational cost

The main factors affecting the accuracy and computational cost of the gauge-independent atomic orbital density functional theory (GIAO-DFT) calculation of (15)N NMR chemical shifts in the representative series of key nitrogen-containing heterocycles--azoles and azines--have been systematically analyzed. In the calculation of (15)N NMR chemical shifts, the best result has been achieved with the KT3 functional used in combination with Jensen's pcS-3 basis set (GIAO-DFT-KT3/pcS-3) resulting in the value of mean absolute error as small as 5 ppm for a range exceeding 270 ppm in a benchmark series of 23 compounds with an overall number of 41 different (15)N NMR chemical shifts. Another essential finding is that basically, the application of the locally dense basis set approach is justified in the calculation of (15)N NMR chemical shifts within the 3-4 ppm error that results in a dramatic decrease in computational cost. Based on the present data, we recommend GIAO-DFT-KT3/pcS-3//pc-2 as one of the most effective locally dense basis set schemes for the calculation of (15)N NMR chemical shifts.

First example of a high-level correlated calculation of the indirect spin-spin coupling constants involving tellurium: tellurophene and divinyl telluride

This paper documents the very first example of a high-level correlated calculation of spin-spin coupling constants involving tellurium taking into account relativistic effects, vibrational corrections and solvent effects for medium sized organotellurium molecules. The (125)Te-(1)H spin-spin coupling constants of tellurophene and divinyl telluride were calculated at the SOPPA and DFT levels, in good agreement with experimental data. A new full-electron basis set, av3z-J, for tellurium derived from the "relativistic" Dyall's basis set, dyall.av3z, and specifically optimized for the correlated calculations of spin-spin coupling constants involving tellurium was developed. The SOPPA method shows a much better performance compared to DFT, if relativistic effects calculated within the ZORA scheme are taken into account. Vibrational and solvent corrections are next to negligible, while conformational averaging is of prime importance in the calculation of (125)Te-(1)H spin-spin couplings. Based on the performed calculations at the SOPPA(CCSD) level, a marked stereospecificity of geminal and vicinal (125)Te-(1)H spin-spin coupling constants originating in the orientational lone pair effect of tellurium has been established, which opens a new guideline in organotellurium stereochemistry.

Toward the Quantum Chemical Calculation of NMR Chemical Shifts of Proteins. 2. Level of Theory, Basis Set, and Solvents Model Dependence

It has been demonstrated that the fragmentation scheme of our adjustable density matrix assembler (ADMA) approach for the quantum chemical calculations of very large systems is well-suited to calculate NMR chemical shifts of proteins [ Frank et al. Proteins2011, 79, 2189-2202 ]. The systematic investigation performed here on the influences of the level of theory, basis set size, inclusion or exclusion of an implicit solvent model, and the use of partial charges to describe additional parts of the macromolecule on the accuracy of NMR chemical shifts demonstrates that using a valence triple-ζ basis set leads to large improvement compared to the results given in the previous publication. Additionally, moving from the B3LYP to the mPW1PW91 density functional and including partial charges and implicit solvents gave the best results with mean absolute errors of 0.44 ppm for hydrogen atoms excluding H(N) atoms and between 1.53 and 3.44 ppm for carbon atoms depending on the size and also on the accuracy of the protein structure. Polar hydrogen and nitrogen atoms are more difficult to predict. For the first, explicit hydrogen bonds to the solvents need to be included and, for the latter, going beyond DFT to post-Hartree-Fock methods like MP2 is probably required. Even if empirical methods like SHIFTX+ show similar performance, our calculations give for the first time very reliable chemical shifts that can also be used for complexes of proteins with small-molecule ligands or DNA/RNA. Therefore, taking advantage of its ab initio nature, our approach opens new fields of application that would otherwise be largely inaccessible due to insufficient availability of data for empirical parametrization.

Optimized basis sets for the calculation of indirect nuclear spin-spin coupling constants involving the atoms B, Al, Si, P, and Cl

The aug-cc-pVTZ-J series of basis sets for indirect nuclear spin-spin coupling constants has been extended to the atoms B, Al, Si, P, and Cl. The basis sets were obtained according to the scheme previously described by Provasi et al. [J. Chem. Phys. 115, 1324 (2001)]. First, the completely uncontracted correlation consistent aug-cc-pVTZ basis sets were extended with four tight s and three tight d functions. Second, the s and p basis functions were contracted with the molecular orbital coefficients of self-consistent-field calculations performed with the uncontracted basis sets on the simplest hydrides of each atom. As a first illustration, we have calculated the one-bond indirect spin-spin coupling constants in BH(4)(-), BF, AlH, AlF, SiH(4), SiF(4), PH(3), PF(3), H(2)S, SF(6), HCl, and ClF at the level of density functional theory using the Becke three parameter Lee-Yang-Parr and the second order polarization propagator approximation with coupled cluster singles and doubles amplitudes.

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 (77)Se-(1)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 (77)Se-(1)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 (77)Se-(1)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.