Ordering of solvent molecules in the vicinity of a dipolar free radical affects its hyperfine coupling constants (hfcs). Specifically, it is demonstrated how the variation of the experimental methylene proton and muon hfcs of the muoniated cyclohexadienyl radical in several solvents and solvent mixtures of varying polarity can be accounted for by a dipole-dipole reaction field model that is based on the model of Reddoch and Konishi (J. Chem. Phys. < BO > 1979 , 70, 2121) which was developed to explain the solvent dependence of the N-14 hfc in the di-tert-butyl-nitroxide radical. Ab initio calculations were carried out with the cyclohexadienyl radical in an electric field to model the electric field arising from the electric dipole moments of the surrounding solvent molecules. An extension of the model that includes the dipole-quadrupole interaction can account for the larger hfc in benzene compared with that in octadecane, and it is predicted that the hfc will be proportional to the concentration of quadrupole moments to the 4/3 power. The influence of hydrogen bonding between the radicals' pi electrons and the OH groups of the solvent on the hfcs is also discussed. Comparison with gas-phase data permits a separation of vibrational effects and reveals that approximately 28% of the temperature dependence in water is due to increasing solvent disorder.