The Importance of Solvent Effects in Calculations of NMR Coupling Constants at the Doubles Corrected Higher Random-Phase Approximation

In this work, 242 NMR spin–spin coupling constants (SSCC) in 20 molecules are calculated, either with correlated wave function methods, SOPPA and HRPA(D), or with density functional theory based on the B3LYP, BHandH, or PBE0 functionals. The calculations were carried out with and without treatment of solvation via a polarizable continuum model in both the geometry optimization step and/or the SSCC calculation, and thereby, four series of calculations were considered (the full-vacuum calculation, the full-solvent calculation, and the two cross combinations). The results were compared with experimental results measured in a solvent. With the goal of reproducing experimental values, we find that the performance of the PBE0 and BHandH SSCCs improves upon including solvation effects. On the other hand, the quality of the B3LYP SSCCs worsens with the inclusion of solvation. Solvation had almost no effect on the performance of the SOPPA and HRPA(D) calculations. We find that the PBE0-based calculations of the spin–spin coupling constants have the best agreement with the experimental data.