PARPHOM: PARallel PHOnon calculator for Moiré systems

The introduction of a twist between two layers of two-dimensional materials has opened up a new and exciting field of research known as twistronics. In these systems, the phonon dispersions show significant renormalization and enhanced electron-phonon interactions as a function of the twist angle. However, the large system size of the resulting moiré patterns in these systems makes phonon calculations computationally challenging. In this paper, we present PARPHOM, a powerful code package designed to address these challenges. PARPHOM enables the generation of force constants, computation of phononic band structures, and determination of density of states in twisted 2D material systems. Moreover, PARPHOM provides essential routines to investigate the finite temperature dynamics in these systems and analyze the chirality of the phonon bands. This paper serves as an introduction to PARPHOM, highlighting its capabilities and demonstrating its utility in unraveling the intricate phononic properties of twisted 2D materials.