Nuclear Magnetic Shielding and Spin−Spin Coupling of 1,2-13C-Enriched Acetylene in Gaseous Mixtures with Xenon and Carbon Dioxide

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.

Absolute 13C Nuclear Magnetic Shielding of Simple Isolated Molecules from Gas Phase Measurements

The new experimental value of 13C absolute shielding constant in an isolated 13CO molecule was evaluated from the 13C and 3He gas phase NMR measurements performed for 3He/13CO mixtures. The...

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.

17O and 1H NMR spectral parameters in isolated water molecules

Small amounts of water enriched in oxygen-17 were studied by 17O and 1H NMR in binary gaseous mixtures with Xe, Kr, CHF3 and CH3F and CO2. The distinct linear dependences of 17O and 1H chemical shifts and 1J(17O,1H) spin-spin coupling on the density of every gas solvent were measured. After the extrapolation of experimental results to zero density the relevant parameters in the isolated H217O molecule were determined. The same procedure was applied for H216O when its proton chemical shift was analyzed but the secondary isotope effect in the 1H shielding of H217O and H216O molecules was too small for detection. As shown, all the intermolecular effects in nuclear magnetic shielding are negative and these effects are more significant for 17O nuclei than for protons. It is consistent with the appropriate gas-to-liquid shifts of water which also indicate deshielding effects for both the investigated nuclei. On the other hand, the 1J0(17O,1H) coupling constant in H217O, which is completely free from intermolecular interactions, considerably differs from the 1J(17O,1H) experimental values obtained for water in liquid solutions. The present experimental data of the isolated H217O molecule are compared with selected results of shielding and spin-spin coupling calculations available from the literature and with the recent experimental data for a water molecule encapsulated in the C60 fullerene. Additionally, on the basis of actual results the magnetic dipole moment of the 17O nucleus is revalued for greater accuracy.

Importance of Triples Contributions to NMR Spin-Spin Coupling Constants Computed at the CC3 and CCSDT Levels

We present the first analytical implementation of CC3 second derivatives using the spin-unrestricted approach. This allows, for the first time, the calculation of nuclear spin-spin coupling constants (SSCC) relevant to NMR spectroscopy at the CC3 level of theory in a fully analytical manner. CC3 results for the SSCCs of a number of small molecules and their fluorine substituted derivatives are compared with the corresponding coupled cluster singles and doubles (CCSD) results obtained using specialized basis sets. For one-bond couplings the change when going from CCSD to CC3 is typically 1-3%, but much higher corrections were found for ¹J_CN in FCN, 15.7%, and ¹J_OF in OF₂, 6.4%. The changes vary significantly in the case of multibond couplings, with differences of up to 10%, and even 13.6% for ³J_FH in fluoroacetylene. Calculations at the coupled cluster singles, doubles, and triples (CCSDT) level indicate that the most important contributions arising from connected triple excitations in the coupled cluster expansion are accounted for at the CC3 level. Thus, we believe that the CC3 method will become the standard approach for the calculation of reference values of nuclear spin-spin coupling constants.

Rovibrational and Temperature Effects in Theoretical Studies of NMR Parameters

The demand for high precision calculations of NMR shieldings (or their related values, chemical shifts δ) and spin–spin coupling constants facilitating and supporting detailed interpretations of NMR spectra increases hand in hand with the development of computational techniques and hardware resources. Highly sophisticated calculations including even relativistic effects are nowadays possible for these properties. However, NMR parameters depend not only on molecular structure and environment but also on molecular flexibility and temperature and the apparent success of theoretical predictions for molecular equilibrium geometries creates a demand for zero-point vibrational and temperature corrections. In this chapter we describe briefly the theory behind rovibrational corrections and review then some important contributions to this field.

Theoretical and experimental NMR studies on muscimol from fly agaric mushroom (Amanita muscaria)

In this article we report results of combined theoretical and experimental NMR studies on muscimol, the bioactive alkaloid from fly agaric mushroom (Amanita muscaria). The assignment of (1)H and (13)C NMR spectra of muscimol in DMSO-d6 was supported by additional two-dimensional heteronuclear correlated spectra (2D NMR) and gauge independent atomic orbital (GIAO) NMR calculations using density functional theory (DFT). The effect of solvent in theoretical calculations was included via polarized continuum model (PCM) and the hybrid three-parameter B3LYP density functional in combination with 6-311++G(3df,2pd) basis set enabled calculation of reliable structures of non-ionized (neutral) molecule and its NH and zwitterionic forms in the gas phase, chloroform, DMSO and water. GIAO NMR calculations, using equilibrium and rovibrationally averaged geometry, at B3LYP/6-31G* and B3LYP/aug-cc-pVTZ-J levels of theory provided muscimol nuclear magnetic shieldings. The theoretical proton and carbon chemical shifts were critically compared with experimental NMR spectra measured in DMSO. Our results provide useful information on its structure in solution. We believe that such data could improve the understanding of basic features of muscimol at atomistic level and provide another tool in studies related to GABA analogs.

Benchmarking density-functional theory calculations of NMR shielding constants and spin-rotation constants using accurate coupled-cluster calculations

Accurate sets of benchmark nuclear-magnetic-resonance shielding constants and spin-rotation constants are calculated using coupled-cluster singles-doubles (CCSD) theory and coupled-cluster singles-doubles-perturbative-triples [CCSD(T)] theory, in a variety of basis sets consisting of (rotational) London atomic orbitals. The accuracy of the calculated coupled-cluster constants is established by a careful comparison with experimental data, taking into account zero-point vibrational corrections. Coupled-cluster basis-set convergence is analyzed and extrapolation techniques are employed to estimate basis-set-limit quantities, thereby establishing an accurate benchmark data set. Together with the set provided for rotational g-tensors and magnetizabilities in our previous work [O. B. Lutnæs, A. M. Teale, T. Helgaker, D. J. Tozer, K. Ruud, and J. Gauss, J. Chem. Phys. 131, 144104 (2009)], it provides a substantial source of consistently calculated high-accuracy data on second-order magnetic response properties. The utility of this benchmark data set is demonstrated by examining a wide variety of Kohn-Sham exchange-correlation functionals for the calculation of these properties. None of the existing approximate functionals provide an accuracy competitive with that provided by CCSD or CCSD(T) theory. The need for a careful consideration of vibrational effects is clearly illustrated. Finally, the pure coupled-cluster results are compared with the results of Kohn-Sham calculations constrained to give the same electronic density. Routes to future improvements are discussed in light of this comparison.

Convergence of Nuclear Magnetic Shieldings in the Kohn-Sham Limit for Several Small Molecules

Convergence patterns and limiting values of isotropic nuclear magnetic shieldings were studied for several small molecules (N2, CO, CO2, NH3, CH4, C2H2, C2H4, C2H6, and C6H6) in the Kohn-Sham limit. Individual results of calculations using dedicated families of Jensen's basis sets (pcS-n and pcJ-n) were fitted toward the complete basis set limit (CBS) using a simple two-parameter formula. Several density functionals were used; calculated vibrational corrections (ZPV) applied; and, for comparison purposes, similar calculations performed using RHF, MP2, SOPPA, SOPPA(CCSD), and CCSD(T) methods and additionally, the aug-cc-pVTZ-J basis set. Finally, the CBS estimated results were critically compared with earlier reported literature data and experimental results. Among 42 studied DFT methods, the KTn and "pure" functionals produced the most accurate heavy atom isotropic nuclear shieldings.

Theoretical and experimental 1H, 13C and 15N NMR studies of N-alkylation of substituted tetrazolo[1,5-a]pyridines

The (15)N as well as (1)H and (13)C chemical shifts of nine substituted tetrazolopyridines and their corresponding tetrazolopyridinium salts have been determined by using NMR spectroscopy at the natural abundance level of all nuclei in CD(3)CN. In this paper, we report, for the first time, the N-alkylation reaction of electron deficient tetrazolopyridines. The treatment of tetrazolopyridines 5-13 with one equivalent of trialkyloxonium tetrafluoroborate leads to a mixture of two isomers, i.e. N3- and N2-alkyl tetrazolo[1,5-a]pyridinium salts. It has been observed that the N3-isomer is always the major isomer, except in the case of the CF(3) substituent, where the two isomers are obtained in the same amount. The quaternary tetrazolopyridinium nitrogen N3 is shielded by around 100 ppm (parts per million) with respect to the parent tetrazolopyridine. Experimental data are interpreted by means of density functional theory (DFT) calculations, including solvent-induced effects, within the conductor-like polarizable continuum model (CPCM). Good agreements between theoretical and experimental (1)H, (13)C and (15)N NMR were found. The combination of multinuclear magnetic resonance spectroscopy with gauge including atomic orbital (GIAO) DFT calculations is a powerful tool in the structural elucidation for both neutral and cationic heterocycles and in the determination of the orientation of N-alkylation of tetrazolopyridines.

Prediction of spin-spin coupling constants in solution based on combined density functional theory/molecular mechanics

We present theory and implementation of calculation of spin-spin coupling constants within combined quantum mechanics/molecular mechanics methods. Special attention is given to the role of explicit solvent polarization as well as the molecular consequences due to hydrogen bonding. The model is generally applicable but is here implemented for the case of density functional theory. First applications to liquid water and acetylene in aqueous solution are presented. Good agreement between theory and experiment is obtained in both cases, thereby showing the strength of our approach. Finally, spin-spin coupling constants across hydrogen bonds are discussed considering for the first time the role of an explicit solvent on this class of spin-spin couplings.

NMR frequency and magnetic dipole moment of (3)He nucleus

We present new gas-phase NMR spectra which relate the resonance frequency of (3)He nucleus to the resonance frequency of the proton in tetramethylsilane (TMS). We discuss the dependence of (3)He resonance frequency on the density of the solvent gas, and we consider in detail the absolute shielding scales of both nuclei. Finally, we analyse the accuracy of the results, using the relationship between the resonance frequencies, absolute shielding constants and magnetic dipole moments of (1)H and (3)He nuclei.

Configurational assignment of carbon, silicon and germanium containing propynal oximes by means of 13C-1H, 13C-13C and 15N-1H spin-spin coupling constants

Configurational assignment of five carbon, silicon and germanium containing propynal oximes has been carried out by means of experimental measurements and high-level ab initio calculations of their 13C-1H, 13C-13C and 15N-1H spin-spin coupling constants. The title compounds were shown to exist in the nonequilibrium mixture of E and Z isomers with the energy difference of less than 0.3 kcal/mol calculated at the MP2/6-311G** level.

Gas-Phase NMR Measurements, Absolute Shielding Scales, and Magnetic Dipole Moments of 29Si and 73Ge Nuclei

New gas-phase NMR measurements of the shielding constants of 29Si, 73Ge, and 1H nuclei in SiH4 and GeH4 are reported. The results, extrapolated to zero density, provide accurate isolated molecule values, best suited for comparison with theoretical calculations. Using the recent ab initio results for these molecules and the measured chemical shifts, we determine the absolute shielding scales for 29Si and 73Ge. This allows us to provide new values of the nuclear magnetic dipole moments for these two nuclei; in addition, we examine the dipole moments of 13C and 119Sn.

Non-empirical calculations of NMR indirect carbon-carbon coupling constants. Part 11--saturated carbocycles: a reference data set and a practical guide to structural elucidation

Carbon-carbon spin-spin coupling constants, J(C,C), calculated at the SOPPA level for 50 mono-, bi- and polycycloalkanes in 10 previous papers are put in order and discussed on unified grounds. Basic structural trends of J(C,C) established in the original publications are summarized and briefly outlined for the representative series. Many unknown couplings are predicted with high reliability, and this provides a good reference data set and a practical guide to the structural elucidation of saturated carbocycles by means of J(C,C) coupling constants.