IR and Raman spectra of isoxazole and D 3 -isoxazole are reassigned with the aid of density functional theory calculations at the B3LYP/6-311+G(d,p) level. Harmonic wavenumber values of all non-CH stretch modes give an r.m.s. deviation from experimental values of 3cm −1 when they are uniformly scaled by 0.9818. Anharmonic values obtained using a second-order perturbative approach show an r.m.s. deviation of 7.6cm −1 . Results of equal quality are obtained for D 3 -isoxazole. Predicted Raman and IR band strengths match experimental spectra very closely. Detailed re-examination of rovibrational constants and inertial defect data from high resolution IR spectra, in comparison with values from the anharmonic frequency calculation, provides confirmation of the vibrational assignments. ‘Dark state’ fundamentals ν 12 (900.2cm −1 ) and ν 15 (866cm −1 ) are detected through the perturbations they cause in nearby bands.