The terahertz frequency spectrum of pentaerythritol tetranitrate (PETN) is calculated using Discover with the COMPASS force field, CASTEP and PWscf. The calculations are compared to each other and to terahertz spectra (0.3-3 THz) of crystalline PETN recorded at 4 K. A number of analysis methods are used to characterise the calculated normal modes.