NMR spin–spin coupling constants in water molecule: equilibrium and rovibrational values

Spin-Spin Coupling Constant Based on Reference Interaction Site Model Self-Consistent Field with Constrained Spatial Electron Density

A method for computing spin-spin coupling constants (SSCCs) using the reference interaction site model self-consistent field with constrained spatial electron density (RISM-SCF-cSED) is proposed for the first time. Describing solvents using integral equation theory allows us to reflect solvent effects at atomic resolution in SSCCs while accounting for thermal fluctuations at a low computational cost. Applying the method to water, 1,1-difluoroethylene, and 1-methylaminomethylene-2-naphthalenone revealed that the solvent shift was evaluated to a greater extent than in the continuum solvent model. The origin of this phenomenon was analyzed in terms of the physical mechanisms underlying SSCCs.

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.

A new radio-frequency acoustic method for remote study of liquids

In the present work, a novel study method of conductive liquids has been proposed. It is based on a discovered phenomenon of radiofrequency anisotropy of electrolyte solution, which arises in response to mechanical excitation of the solution. The phenomenon was observed during the development of a radiofrequency polarimetric contactless cardiograph. The electric field vector rotates after its transition through the pericardial region due to the acceleration changes of blood. Numerous in vitro experiments with monochromatic and impulse acoustic waves always induced the polarization rotation of the RF wave passing through an electrolyte solution. The response obtained from the solutions on acoustic excitation of the Heaviside function form demonstrates the effect of a solution “memory”. The dynamics of this process resembles the spin glasses magnetization. We hypothesized that there was a magnetic moment change within the solution, and the possible reason for it is an appearance of electromagnetic impulse caused by the same acoustic excitation. In a further experiment, we really captured a suspected electrical potential. Given that, we can declare at least three new physical effects never observed before for an electrolyte solution. The study method itself may provide broad options for remote measurement of the electrolyte solution parameters.

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.

Natural bond orbital/natural J-coupling study of vicinal couplings

NBO-NJC decomposition of vicinal (3)J HH spin-spin coupling constants into Lewis, delocalization, and repolarization contributions are presented. A deep study allows to assign the main contributions to specific orbitals or electron delocalizations between two orbitals. (3)J HH torsional dependence and the substituent effect are analyzed according to the main orbital contributions for ethane and fluoroethane molecules using different basis sets. The torsional dependence for the energies corresponding to electron delocalization is also studied.

Benchmarking NMR indirect nuclear spin-spin coupling constants: SOPPA, SOPPA(CC2), and SOPPA(CCSD) versus CCSD

Accurate calculations of NMR indirect nuclear spin-spin coupling constants require especially optimized basis sets and correlated wave function methods such as CCSD or SOPPA(CCSD). Both methods scale as N(6), where N is the number of orbitals, which prevents routine applications to molecules with more than 10-15 nonhydrogen atoms. We have therefore developed a modification of the SOPPA(CCSD) method in which the CCSD singles and doubles amplitudes are replaced by CC2 singles and doubles amplitudes. This new method, called SOPPA(CC2), scales only as N(5), like the original SOPPA-method. The performance of the SOPPA(CC2) method for the calculation of indirect nuclear spin-spin coupling constants is compared to SOPPA and SOPPA(CCSD) employing a set of benchmark molecules. We also investigate the basis set dependence by employing three different basis sets optimized for spin-spin coupling constants, namely the HuzIV-su4, ccJ-pVTZ, and ccJ-pVQZ basis sets. The results of the corresponding CCSD calculations are used as a theoretical reference.

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.

H2O, H2, HF, F2 and F2O nuclear magnetic shielding constants and indirect nuclear spin-spin coupling constants (SSCCs) in the BHandH/pcJ-n and BHandH/XZP Kohn-Sham limits

Good performance of segmented contracted basis sets XZP, where X = D, T, Q and 5, for obtaining H(2)O, H(2), HF, F(2) and F(2)O nuclear isotropic shielding constants in the BHandH Kohn-Sham basis set limit was shown. The results of two- and three-parameter complete basis set limit extrapolation schemes were compared with experimental results, earlier literature data and benchmark ab initio results. Similar convergence patterns of shieldings obtained from calculations using general purpose XZP basis sets and from polarization-consistent basis sets pcS-n and pcJ-n, where n = 0, 1, 2, 3 and 4, designed to accurately predict magnetic properties were observed. On the contrary, the SSCCs were more sensitive to the XZP basis set size and generally less accurate than those estimated using pcJ-n basis set family. The BHandH density functional markedly outperforms B3LYP method in predicting heavy atom shieldings and SSCCs values in the studied systems.

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.

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. 1J(CH) and 3J(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.

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 (1)H, (17)O isotropic shieldings and (1)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 (1)J(OH) is similar to the more advanced methods, second-order polarization propagator approximation (SOPPA) and SOPPA(CCSD), in the basis set limit.

Complete basis set B3LYP NMR calculations of CDCl3 solvent's water fine spectral details

The assignment of singlet at 1.55 ppm and the 1:1:1 triplet at 1.519 ppm to H(2)O and HOD in the 400 MHz (1)H NMR spectrum of CDCl(3) solvent were supported by complete basis set (CBS) GIAO-B3LYP calculated chemical shift and the CBS B3LYP estimated (2)J(D,H) spin-spin coupling constant (SSCC). The CBS fitting of B3LYP/cc-pCVxZ and B3LYP/pcJ-n predicted SSCC values, the accurate value of (2)J(D,H) = -1.082 +/- 0.030 Hz of HOD in chloroform-d(1) and the H/D isotopic shift of 0.0307(1) ppm were reported for the first time. The agreement between CBS B3LYP predicted chemical shift, spin-spin values and experiment was good.

Vicinal fluorine-fluorine coupling constants: Fourier analysis

Stereochemical dependences of vicinal fluorine-fluorine nuclear magnetic resonance coupling constants (3JFF) have been studied with the multiconfigurational self-consistent field in the restricted active space approach, with the second-order polarization propagator approximation (SOPPA), and with density functional theory. The SOPPA results show the best overall agreement with experimental couplings. The relationship with the dihedral angle between the coupled fluorines has been studied by Fourier analysis, the result is very different from that of proton-proton couplings. The Fourier coefficients do not resemble those of a typical Karplus equation. The four nonrelativistic contributions to the coupling constants of 1,2-difluoroethane configurations have been studied separately showing that up to six Fourier coefficients are required to reproduce the calculated values satisfactorily. Comparison with Fourier coefficients for matching hydrogen fluoride dimer configurations suggests that the higher order Fourier coefficients (Cn> or =3) originate mainly from through-space Fermi contact interaction. The through-space interaction is the main reason 3JFF do not follow the Karplus equation.

Density functional theory study of indirect nuclear spin-spin coupling constants with spin-orbit corrections

This work outlines the calculation of indirect nuclear spin-spin coupling constants with spin-orbit corrections using density functional response theory. The nonrelativistic indirect nuclear spin-spin couplings are evaluated using the linear response method, whereas the relativistic spin-orbit corrections are computed using quadratic response theory. The formalism is applied to the homologous systems H2X (X=O,S,Se,Te) and XH4 (X=C,Si,Ge,Sn,Pb) to calculate the indirect nuclear spin-spin coupling constants between the protons. The results confirm that spin-orbit corrections are important for compounds of the H2X series, for which the electronic structure allows for an efficient coupling between the nuclei mediated by the spin-orbit interaction, whereas in the case of the XH4 series the opposite situation is encountered and the spin-orbit corrections are negligible for all compounds of this series. In addition we analyze the performance of the density functional theory in the calculations of nonrelativistic indirect nuclear spin-spin coupling constants.