The calculation of indirect nuclear spin-spin coupling constants in large molecules

We present calculations of indirect nuclear spin-spin coupling constants in large molecular systems, performed using density functional theory. Such calculations, which have become possible because of the use of linear-scaling techniques in the evaluation of the Coulomb and exchange-correlation contributions to the electronic energy, allow us to study indirect spin-spin couplings in molecules of biological interest, without having to construct artificial model systems. In addition to presenting a statistical analysis of the large number of short-range coupling constants in large molecular systems, we analyse the asymptotic dependence of the indirect nuclear spin-spin coupling constants on the internuclear separation. In particular, we demonstrate that, in a sufficiently large one-electron basis set, the indirect spin-spin coupling constants become proportional to the inverse cube of the internuclear separation, even though the diamagnetic and paramagnetic spin-orbit contributions to the spin-spin coupling constants separately decay as the inverse square of this separation. By contrast, the triplet Fermi contact and spin-dipole contributions to the indirect spin-spin coupling constants decay exponentially and as the inverse cube of the internuclear separation, respectively. Thus, whereas short-range indirect spin-spin coupling constants are usually dominated by the Fermi contact contribution, long-range coupling constants are always dominated by the negative diamagnetic spin-orbit contribution and by the positive paramagnetic spin-orbit contribution, with small spin-dipole and negligible Fermi contact contributions.     

Algebraic-diagrammatic construction polarization propagator approach to indirect nuclear spin-spin coupling constants

A new polarization propagator approach to indirect nuclear spin-spin coupling constantans is formulated within the framework of the algebraic-diagrammatic construction (ADC) approximation and implemented at the level of the strict second-order approximation scheme, ADC(2). The ADC approach possesses transparent computational procedure operating with Hermitian matrix quantities defined with respect to physical excitations. It is size-consistent and easily extendable to higher orders via the hierarchy of available ADC approximation schemes. The ADC(2) method is tested in the first applications to HF, N(2), CO, H(2)O, HCN, NH(3), CH(4), C(2)H(2), PH(3), SiH(4), CH(3)F, and C(2)H(4). The calculated indirect nuclear spin-spin coupling constants are in good agreement with the experimental data and results of the second-order polarization propagator approximation method. The computational effort of the ADC(2) scheme scales as n(5) with respect to the number of molecular orbitals n, which makes this method promising for applications to larger molecules.     

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.     

Basis-set dependence of nuclear spin-spin coupling constants

The convergence of NMR indirect spin-spin coupling constants with the extension of the basis set is analyzed, based on calculations carried out at the multiconfigurational self-consistent-field level for the HF and H₂O systems. For the dominant and difficult Fermi-contact contribution, the standard correlation-consistent basis sets of electronic-structure theory are not suitable, lacking flexibility in the core region. Improved but not satisfactory convergence of the couplings is observed when decontracting the s functions of the correlation-consistent cc-pVXZ basis sets for 2≤X≤6. Next, by systematically extending these basis sets with tight s functions, new sets are obtained that are sufficiently flexible for accurate calculations of indirect nuclear spin-spin couplings, without sacrificing the smooth convergence behavior of the correlation-consistent basis sets. Received: 22 September 1997 / Accepted: 30 December 1997     

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.     

NDDO MO calculations: isotropic hyperfine coupling constants and nuclear spin-spin coupling constants

The nonempirical NDDO MO method in its unrestricted form has been used to evaluate isotropic hyperfine coupling constants and nuclear spin-spin coupling constants. Satisfactory agreement with INDO and experimental results is obtained.     

Calculation of indirect nuclear spin-spin coupling constants within the regular approximation for relativistic effects

A new method for calculating the indirect nuclear spin-spin coupling constant within the regular approximation to the exact relativistic Hamiltonian is presented. The method is completely analytic in the sense that it does not employ numeric integration for the evaluation of relativistic corrections to the molecular Hamiltonian. It can be applied at the level of conventional wave function theory or density functional theory. In the latter case, both pure and hybrid density functionals can be used for the calculation of the quasirelativistic spin-spin coupling constants. The new method is used in connection with the infinite-order regular approximation with modified metric (IORAmm) to calculate the spin-spin coupling constants for molecules containing heavy elements. The importance of including exact exchange into the density functional calculations is demonstrated.     

Stereochemical dependence of NMR geminal spin-spin coupling constants

In this work it was sought to explore the versatility of geminal spin-spin coupling constants, (2)J(XY) SSCCs, as probes for stereochemical studies. A set of compounds, where their experimental (2)J(XY) SSCCs through the X-C-Y molecular fragment are predicted to be sensitive to hyperconjugative interactions involving either bonding or antibonding orbitals containing the C carbon atom ('coupling pathway'), were analyzed. SSCC calculations were performed for some selected examples using the second order polarization propagator approximation (SOPPA) method or within the DFT-B3LYP framework. Hyperconjugative interactions were calculated within the Natural Bond Orbital (NBO) approach. Results are condensed in two qualitative rules: Rule I(M)-hyperconjugative interactions transferring charge into the coupling pathway yield a positive increase to the Fermi contact (FC), contribution to (2)K(XY) reduced spin-spin coupling constants (RSSCC), and Rule II(M)-hyperconjugative interactions transferring charge from the coupling pathway yield a negative increase to the FC contribution to (2)K(XY) RSSCC.     

MBOHO calculations of phosphorus-carbon nuclear spin-spin coupling constants

Summary The novel generalized correlation of the nuclear spin-spin coupling constants with the atomic hybrids and net charges is employed to give a new relationship for calculating the directly bonded phosphorus-carbon coupling constants by use of the maximum bond order hybrid orbital procedure together with the extended Hückel molecular orbital calculation. The calculated coupling constants of phosphorus-carbon are all in good agreement with the experimental data, which shows that the new relationship obtained in the present paper is quite satisfactory for calculation of the phosphorus-carbon coupling constants.The project was supported by the National Natural Science Foundation of China and the Excellent University Teacher's Foundation of State Education Commission of China     

Thermal and solvent effects on NMR indirect spin-spin coupling constants of a prototypical Chagas disease drug

NMR J-couplings across hydrogen bonds reflect the static and dynamic character of hydrogen bonding. They are affected by thermal and solvent effects and can therefore be used to probe such effects. We have applied density functional theory (DFT) to compute the NMR (n)J(N,H) scalar couplings of a prototypical Chagas disease drug (metronidazole). The calculations were done for the molecule in vacuo, in microsolvated cluster models with one or few water molecules, in snapshots obtained from molecular dynamics simulations with explicit water solvent, and in a polarizable dielectric continuum. Hyperconjugative and electrostatic effects on spin-spin coupling constants were assessed through DFT calculations using natural bond orbital (NBO) analysis and atoms in molecules (AIM) theory. In the calculations with explicit solvent molecules, special attention was given to the nature of the hydrogen bonds formed with the solvent molecules. The results highlight the importance of properly incorporating thermal and solvent effects into NMR calculations in the condensed phase.     

Illumination of the effect of the overlap of lone-pairs on indirect nuclear spin-spin coupling constants

The effect of electron lone-pairs on the Fermi-contact (FC) contribution to indirect nuclear spin-spin coupling constants is analyzed using new tools for their interpretation. In particular, visualization of spin-spin coupling pathways using the coupling deformation density (CDD) has been employed. Furthermore, the recently developed perturbation-stable localization procedure has been applied for decomposition of CDD and the calculated value of couplings into contributions from localized molecular orbitals (LMOs). Correlation between the overlap of densities of LMOs representing lone-pairs and the Fermi-contact contribution to spin-spin coupling constants has been demonstrated. A new way for analyzing spin-spin couplings using the expansion of CDD as a linear combination of the products of molecular orbitals has been suggested. The considered examples include two- and three-bond phosphor-phosphor couplings. Significance of the obtained insight is not restricted to spin-spin couplings of nuclei possessing lone-pairs, as demonstrated in the example of vicinal hydrogen-hydrogen coupling in ethane.     

NMR spin-spin coupling constants in polymethine dyes as polarity indicators

Herein, we explore the use of spin-spin coupling constants (SSCCs) in merocyanine (MCYNE) dyes as indicators of polarity. For this purpose, we use Car-Parrinello hybrid quantum mechanics/molecular mechanics (QM/MM) to determine the structures of MCYNE in solvents of different polarity, followed by computations of the SSCCs by using QM/MM linear-response theory. The molecular geometry of MCYNE switches between neutral, cyanine-like, and zwitterionic depending on the polarity of the solvent. This structural variation is clearly reflected in the proton SSCCs in the polymethine backbone, which are highly sensitive to the dielectric nature of the environment; this mechanism can be used as a "polarity indicator" for different microenvironments. This result is highlighted by computing the SSCCs of the MCYNE probe in the cavity of the beta-lactoglobulin protein. The computed SSCCs clearly indicate a non-polar hydrophobic dielectric nature of this cavity.     

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.     

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.     

A review with comprehensive data on experimental indirect scalar NMR spin-spin coupling constants across hydrogen bonds

Scalar NMR spin-spin coupling constants across hydrogen bonds are fundamental in structural studies and as test grounds for theoretical calculations. Since they are scattered among many articles of different kinds, it seems useful to collect them in the most comprehensive way.     

Ab initio study of the NMR shielding constants and spin-spin coupling constants in cyclopropene

Summary Ab initio calculations of parameters which characterize the NMR spectrum are presented for the cyclopropene molecule. The London orbitals CHF (or GIAO-CHF, Gauge-Independent Atomic Orbital Coupled Hartree-Fock) results for the shielding constants are in good agreement with the experimental data, accurately determined, and with otherab initio values. The calculations of the NMR spin-spin coupling constants have been performed using the Multiconfiguration Time-Dependent Hartree-Fock (MC TDHF) approach. Different basis sets and MC SCF wavefunctions were used to estimate the accuracy of the results. Good agreement is obtained with the coupling constants estimated using the available experimental data.Dedicated to Professor Werner Kutzelnigg on the occasion of his 60th birthday     

Coupled hartree-fock and second order polarization propagator calculations of indirect nuclear spin-spin coupling constants for diatomic molecules

We report ab initio calculation of indirect nuclear spin-spin coupling constants for the HD, FH, CO and CH⁺ molecules using the first (coupled Hart     

Density functional theory calculation of indirect nuclear magnetic resonance spin-spin coupling constants in C(70)

We calculate NMR spin-spin coupling constants in the C70 fullerene by means of density functional theory. We show that using a hybrid density functional (B3LYP) and an adequate basis set (cc-pCVDZ-sd), excellent agreement with experimental values can be achieved for one-bond couplings. These benchmark calculations suggest that theoretical predictions of NMR spin-spin couplings can be extremely valuable for discerning structural information of fullerenes.     

Ab initio study of the J(CC) indirect nuclear spin-spin coupling constants in cyclobutane and related systems

Non-empirical calculations using the equations-of-motion approach, which incorporates the main portion of the electron correlation effects, are reported for the carbon-carbon nuclear spin-spin coupling constants in cyclobutane, bicyclobutane, tricyclobutane, cyclobutene, cyclobutyne, cyclobutadiene, bicyclobutene, methylenecyclopropane, and methylenecyclopropene. The results provide an overall picture of the influences exerted on sign and magnitude of the J(CC) by progressive condensation, unsaturation, and branching rearrangement of the cyclobutane frame.