FourPhonon: An extension module to ShengBTE for computing four-phonon scattering rates and thermal conductivity

Abstract: FourPhonon is a computational package that can calculate four-phonon scattering rates in crystals. It is built within ShengBTE framework, which is a well-recognized lattice thermal conductivity solver based on Boltzmann transport equation. An adaptive energy broadening scheme is implemented for the calculation of four-phonon scattering rates. In analogy with $thirdorder.py$ in ShengBTE, we also provide a separate python script, $Fourthorder.py$, to calculate fourth-order interatomic force-constants. The extension module preserves all the nice features of the well-recognized lattice thermal conductivity solver ShengBTE, including good parallelism and straightforward workflow. In this paper, we discuss the general theory, program design, and example calculations on Si, BAs and $\mathrm{LiCoO_2}$.

• The paper introduces FourPhonon as an extension to ShengBTE to compute critical four-phonon scattering rates for precise thermal conductivity predictions.
• It leverages an adaptive energy broadening scheme and symmetry operations to efficiently calculate fourth-order interatomic force constants via first-principles methods.
• Case studies on Si, BAs, and LiCoO2 demonstrate the module’s impact on refining thermal models by capturing significant higher-order scattering effects.

Insights into the FourPhonon Computational Module for Thermal Conductivity Calculations

This essay provides a comprehensive overview of the paper "FourPhonon: An extension module to ShengBTE for computing four-phonon scattering rates and thermal conductivity." This work introduces the FourPhonon computational package, designed as an extension to the well-regarded ShengBTE framework. FourPhonon facilitates the calculation of four-phonon scattering rates in crystals, an aspect critical for accurate lattice thermal conductivity predictions, especially in materials where higher-order anharmonic interactions are non-negligible.

Context and Motivation

The understanding of phonon-phonon interactions is essential, as phonons serve as the primary heat carriers in various materials such as insulators, semiconductors, and some semimetals. Traditional models primarily focused on three-phonon scattering events. However, recent advances suggest that higher-order interactions, specifically four-phonon scattering, significantly influence thermal transport properties in certain materials, such as BAs (boron arsenide). These insights spurred the development of FourPhonon, addressing the computational inaccessibility of these higher-order interactions due to their complexity and intensive computational demands.

Methodology

FourPhonon is seamlessly integrated within the ShengBTE framework, exploiting its parallel computing capabilities. The package employs the Boltzmann transport equation (BTE) for phonons, incorporating an adaptive energy broadening scheme for four-phonon scattering rate calculations. This methodology hinges on the calculation of fourth-order interatomic force constants (IFCs) using first-principles approaches, facilitated by a newly introduced Python script, Fourthorder.py. This script leverages symmetry operations to reduce computational load, a necessary step given the high demands of four-phonon calculations.

Computational Insights

The paper illustrates the utility of FourPhonon through case studies in silicon (Si), boron arsenide (BAs), and lithium cobalt oxide (LiCoO2). Each example highlights the distinct aspects of four-phonon interactions:

1. Silicon: The study confirms that while three-phonon processes dominate, four-phonon interactions are not insignificant, particularly at higher temperatures. The computational demands are alleviated with a coarser q-point mesh for four-phonon calculations compared to three-phonon ones.
2. Boron Arsenide: BAs showcases the prominence of four-phonon scattering, resulting in substantial reductions in thermal conductivity estimates compared to predictions considering solely three-phonon events. The phase space of four-phonon interactions in BAs is notably larger, underlying their significant role in thermal transport.
3. Lithium Cobalt Oxide: The work on LiCoO2, a complex structured material, displays the comprehensive capabilities of FourPhonon in dealing with intricate systems. The study outlines the decomposition of scattering rates into different channels, reinforcing the nuanced understanding of anharmonic interactions.

Theoretical and Practical Implications

The findings elucidate the necessity of incorporating four-phonon scattering in thermal conductivity predictions when dealing with materials exhibiting strong anharmonicity. This advancement represents an essential tool for theoretical exploration and practical applications, particularly in thermal management and materials science, impacting fields such as semiconductor technology and energy storage solutions.

Future Directions

While the current implementation of FourPhonon operates effectively under the Single Mode Relaxation Time Approximation (SMRTA), future developments aim to fully integrate four-phonon interactions within an iterative BTE scheme. This progression is contingent on resolving substantial memory demands, a technical challenge highlighted by the authors.

In conclusion, the FourPhonon module extends the capability of phonon transport simulations, offering researchers a valuable resource for exploring complex phononic interactions. The inclusion of higher-order phonon scattering calculations aligns with the continuous evolution of thermal conductivity prediction methodologies, reinforcing the significance of advanced computational frameworks in material science.