Potential functions and harmonic (omega(i)) and anharmonic (nu(i)) fundamental frequencies have been calculated for 1,1-difluorocyclopropane (DFCP) and its d(4) and d(2) isotopomers using the program Gaussian 03. B3LYP and MP2 models were employed, each with the bases 6-311++G** and cc-pVTZ. Anharmonicity corrections Delta(i) = omega(i) - nu(i) are listed and shown to be different for symmetric and antisymmetric CH stretching modes in situations where Fermi resonance appears to be absent. The same effect is missing in C2H4, for which similar calculations were made. The quadratic force fields for DFCP have been scaled in symmetry coordinate space with 15 scale factors both to observed frequencies nu(obsd) and also to omega(obsd), where omega(obsd) = nu(obsd) + Delta. With nu(obsd) especially, different scale factors are needed for the symmetric and antisymmetric CH stretching force constants due to their differing anharmonicities. The source of the latter in the quartic and cubic force field is explored. MP2 calculations of valence interaction force constants involving the stretching of bonds on a common carbon atom are preferred to those from a B3LYP model. In either model, scaling to omega(obsd) rather than to nu(obsd) does not remove the necessity of varying scale factors for differing types of motion in the same group. Theoretical values of the five quartic centrifugal distortion constants are listed for the normal species and compared with new experimental data. The predictions are sufficiently good to be useful in fitting pure rotational transitions. A weakness is identified in the current Gaussian 03 code for the calculation of vibration-rotation quantities, and limitations are noted in the manner in which Fermi resonance is handled.