PARPHOM: PARallel PHOnon calculator for Moiré systems

Shinjan Mandal,Indrajit Maity,H R Krishnamurthy,Manish Jain

2024-10-28

Abstract: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.

Mesoscale and Nanoscale Physics,Materials Science,Computational Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: in distorted two - dimensional material systems (i.e., moiré systems), phonon calculations are extremely challenging due to the large number of atoms. Existing computational methods and tools have difficulty dealing with properties such as phonon dispersion and density of states in these large - scale systems. Especially when considering the atomic relaxation effect, the computational cost is extremely high. Specifically:

1. **Computational problems brought by large system sizes**: The number of atoms in moiré patterns can reach hundreds to thousands, making traditional first - principles calculations based on density functional theory (DFT) impractical. Such calculations require extremely high memory resources, and as the system size increases, the computational complexity grows exponentially.

2. **Limitations of existing tools**: Currently available codes (such as PHONOPY or PhonoLAMMPS) can only run on a single core when generating force constants, which greatly limits their ability to handle large - scale moiré systems.

To solve these problems, the author has developed a parallel computing package named PARPHOM. The main contributions of PARPHOM include:

- **Parallelizing force constant generation**: By using a parallel algorithm to generate the force constant matrix, the calculation time is significantly reduced.
- **Efficient large - scale system handling ability**: It can handle moiré systems containing thousands of atoms and calculate their phonon band structures and density of states.
- **Temperature - dependence and chirality analysis**: It provides tools for studying finite - temperature dynamics and the chirality of phonon modes.

Therefore, PARPHOM aims to provide an efficient solution to meet the high computational requirements for phonon calculations in moiré systems, thereby promoting the research progress in this field.