The isotropic hyperfine coupling constants of several organic radicals including CH 3 , CH 2 , CH 2 Ϫ , C 2 H 5 , C 2 H 3 , H 2 CN, C 6 H 7 , and C 3 H 5 are calculated analytically using the coupled cluster ͑CC͒ ''relaxed density'' matrix approach. We employ three different commonly used basis sets with CCSD and CCSD͑T͒ in order to calibrate expected accuracy. The Chipman basis set combined with the CCSD͑T͒ method performs best for carbon isotropic hyperfine coupling constants with a mean absolute deviation within 8% compared to experiment. The corresponding mean absolute deviation for hydrogen isotropic hyperfine coupling constants from experiment is 12%. We show that the UHF, ROHF, and quasi ͑QRHF͒ reference function CCSD spin densities are effectively numerically equivalent in the notorious case of the allyl radical.