Karplus relationships for all 26 NMR spin-spin coupling constants (SSCCs) J of the pseudorotating tetrahydrofuran (THF) molecule were derived by expanding J as a function J(q,φ) of the puckering amplitude q and the pseudorotational phase angle φ. For this purpose, the conformational potential V(q,φ) of THF was determined at the MBPT(2)/ccpVTZ and B3LYP/6-31G(d,p) levels of theory. THF is a slightly hindered pseudorotor (MP2 barriers ∆E and ∆H are 0.1 -0.2 kcal/mol) with a puckering amplitude q = 0.396Å and a barrier to inversion ∆H = 4.1 kcal/mol in reasonable agreement with experimental data. The SSCCs of THF were calculated both at MBPT(2)/cc-pVTZ and B3LYP/6-31G(d,p) geometries using coupled perturbed density functional theory (CPDFT) with the B3LYP functional and a (9s,5p,1d/5s,1p)[6s,4p,1d/3s,1p] basis set. All geometrical parameters and the 26 SSCCs of THF were computed as functions of the phase angle φ and averaged to give