First-principles-based calculations of vibrational normal modes in polyatomic materials with translational symmetry: application to PETN molecular crystal.

Numerical studies of vibrational energy transport and associated (non)linear infrared and Raman response in polyatomic materials require knowledge of the multidimensional vibrational potential-energy surface and the ability to perform normal-mode analysis on that potential. The presence of translational symmetry, as in crystals, leads to the observed dispersion of the unit cell normal modes and has to be accounted for in calculations of energy transfer rates and other spectroscopic quantities. Here we report on the implementation of a computational approach that combines the generalized supercell method and density functional theory electronic structure calculations to investigate the vibrational structure in translationally symmetric materials containing relatively large numbers of atoms in the unit cell (58 atoms in the present study). The method is applied to calculate the phonon and vibron dispersion relations and the vibrational density of states in pentaerythritol tetranitrate (PETN) molecular crystal which is an important energetic material. The results set the stage for future investigations of vibrational energy transport and associated nonlinear spectroscopic signatures in this class of materials.