The computation of Karplus equation coefficients and their components using self-consistent field and second-order polarization propagator methods

The Karplus equation has been investigated by ab initio computation of the spin-spin coupling constants for a series of rotated ethane geometries. The couplings have been calculated at the self-consistent field (SCF) level as well as using the second-order polarization propagator approximation (SOPP A) and the second-order polarization propagator approximation with coupled cluster singles and doubles amplitudes (SOPPA(CCSD)) and have been compared with results of previous calculations. The four principal components of th e coupling constants rather than just the Fermi-contact have been calculate d, and the common supposition that the Fermi-contact term is totally domina nt has been confirmed. The derivatives of the orbital paramagnetic and orbi tal diamagnetic terms are significant but opposite in sign for the case of this rotation in ethane. It is found that the coefficients in the Karplus e quation are largely overestimated at the SCF level, whereas the SOPPA(CCSD) results are in good agreement with coefficients derived from experimental coupling constant data or the results of multiconfigurational self-consiste nt field (MCSCF) calculations. It is further observed that extending the Fo urier series in the Karplus equation to include cos(3 theta) and cos(4 thet a) terms neither significantly improves the quality of the rt nor significa ntly changes the values of the other coefficients. In order to simulate the Abraham and Pachler equation, calculations varying the nuclear charge on h ydrogen have been performed. These will allow an abstract but flexible pred iction of the effect of electronegative substituents.