Evaluating Computational Shortcuts in Supercell-Based Phonon Calculations of Molecular Crystals: The Instructive Case of Naphthalene

Tomas Kamencek,Sandro Wieser,Hirotaka Kojima,Natalia Bedoya-Martínez,Johannes P. Dürholt,Rochus Schmid,Egbert Zojer
DOI: https://doi.org/10.1021/acs.jctc.0c00119
2020-04-06
Abstract:Phonons crucially impact a variety of properties of organic semiconductor materials. For instance, charge- and heat transport depend on low-frequency phonons, while for other properties, such as the free energy, especially high-frequency phonons count. For all these quantities one needs to know the entire phonon band structure, whose simulation becomes exceedingly expensive for more complex systems when using methods like dispersion-corrected density functional theory (DFT). Therefore, in the present contribution we evaluate the performance of more approximate methodologies, including density functional tight binding (DFTB) and a pool of force fields (FF) of varying complexity and sophistication. Beyond merely comparing phonon band structures, we also critically evaluate to what extent derived quantities, like temperature-dependent heat capacities, mean squared thermal displacements and temperature-dependent free energies are impacted by shortcomings in the description of the phonon bands. As a benchmark system, we choose (deuterated) naphthalene, as the only organic semiconductor material for which to date experimental phonon band structures are available in the literature. Overall, the best performance amongst the approximate methodologies is observed for a system-specifically parametrized second-generation force field. Interestingly, in the low-frequency regime also force fields with a rather simplistic model for the bonding interactions (like the General Amber Force Field) perform rather well. As far as the tested DFTB parametrization is concerned, we obtain a significant underestimation of the unit cell volume resulting in a pronounced overestimation of the phonon energies in the low frequency region. This cannot be mended by relying on the DFT-calculated unit cell, since with this unit cell the DFTB phonon frequencies significantly underestimate the experiments.
Materials Science
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to evaluate the effectiveness of computational shortcuts in calculating phonons in molecular crystals in supercell groups. Specifically, the paper focuses on the influence of phonons on various properties in organic semiconductor materials, such as charge transport, heat transport, and free energy. These properties depend on the knowledge of the entire phonon band structure, and when using methods such as dispersion - corrected density functional theory (DFT) to simulate complex systems, the computational cost is very high. Therefore, this article evaluates more approximate methods, including density functional tight - binding (DFTB) and a series of classical force fields (FF) with different complexities, to determine how these methods perform when calculating the phonon band structure and its derived properties (such as temperature - dependent heat capacity, root - mean - square thermal displacement, and temperature - dependent free energy). The paper selects deuterated naphthalene as the benchmark system because it is currently the only organic semiconductor material for which experimental phonon band structure data are available in the literature. Through this research, the authors aim to explore the reliability and accuracy of these approximate methods in describing the phonon bands in the low - frequency and high - frequency regions.