We use density functional tight-binding (DFTB) theory to calculate the surface energies of two energetic crystals: monoclinic β-1, 3,5,3,5, and tetragonal pentaerythritol tetranitrate (PETN). The results are then employed to determine crystal shapes using the Bravais− Friedel−Donnay−Harker, attachment energy, and surface energy models. We find that energy-based models yield predictions in good agreement with experimental observations. Additionally, we propose a simple model that reframes surface energy as a measure of the lost intermolecular interactions during the formation of a surface from the bulk. The model accurately captures the results from the DFTB calculations and enables us to explain and predict surface energies as a function of the local molecular environment.