pith. sign in

arxiv: 2605.05620 · v1 · submitted 2026-05-07 · ❄️ cond-mat.mtrl-sci · physics.comp-ph

Physics-Grounded Understanding of Thermal Boundary Conductance between Ga₂O₃ and SiC from a Feedforward Neural Network Potential

Pith reviewed 2026-05-08 08:57 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci physics.comp-ph
keywords thermal boundary conductanceGa2O3/SiC heterostructureneural network potentialmolecular dynamics simulationinterface phonon transportbeta-Ga2O34H-SiCnonequilibrium heat transport
0
0 comments X

The pith

A unified feedforward neural network potential trained on bulk data predicts thermal boundary conductance at Ga2O3/SiC interfaces and explains its dependence on orientation, length, and temperature.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a single feedforward neural network potential that describes both gallium oxide and silicon carbide using training data focused on their bulk properties. This potential is then used in nonequilibrium molecular dynamics simulations to compute heat flow across the interfaces at scales that are too large for direct quantum calculations. The results show that thermal boundary conductance drops as the distance over which heat travels increases, rises as temperature goes up, and is higher when the gallium oxide is oriented along its (201) face than along its (100) face. These patterns are traced to how phonons with different travel distances interact at the boundary and to differences in bonding strength between the two orientations.

Core claim

A transferable feedforward neural network potential, validated against density-functional theory data, bulk phonon dispersions, and anisotropic thermal conductivities, enables nonequilibrium molecular dynamics simulations of thermal boundary conductance. The simulations establish that conductance decreases with transport length through attenuation of long-mean-free-path carriers, increases with temperature through enhanced incoherent and anharmonic exchange within stable spectral channels, and remains higher for the Ga2O3(201)/SiC(0001) interface than for Ga2O3(100)/SiC(0001) because of stronger bonding and vibrational coupling at the former.

What carries the argument

The feedforward neural network potential that unifies oxide and carbide bonding environments, trained primarily on bulk DFT data and phonon properties, and deployed in large-scale nonequilibrium molecular dynamics to compute interfacial heat transport.

If this is right

  • Thermal boundary conductance decreases with increasing transport length because long-mean-free-path phonons are attenuated before they reach the interface.
  • Thermal boundary conductance increases with temperature because incoherent and anharmonic processes at the interface become more effective while the main frequency channels stay the same.
  • The Ga2O3(201)/SiC(0001) interface conducts heat better than the (100) interface because of stronger atomic bonding and better vibrational overlap across the boundary.
  • A single neural network potential trained on bulk data is sufficient to predict interfacial transport in systems where separate empirical potentials would be unreliable.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same potential could be used to screen other crystal orientations or to add controlled defects and see how they alter heat flow without retraining from scratch.
  • If the orientation dependence holds in real devices, engineers could select substrate cuts that maximize heat extraction from Ga2O3 power electronics.
  • The length dependence implies that thinner interfacial layers or shorter heat paths in devices would yield higher effective conductance than bulk-interface models predict.
  • Extending the approach to other oxide-carbide pairs might reveal general rules for choosing interface planes that improve thermal performance.

Load-bearing premise

The neural network potential trained mainly on bulk properties transfers accurately to the atomically sharp, defect-free Ga2O3/SiC interface without introducing significant errors in the predicted heat flow.

What would settle it

Experimental measurements of thermal boundary conductance on high-quality, defect-free Ga2O3/SiC interfaces that show no difference between the (201) and (100) orientations, or that fail to match the simulated length and temperature trends, would falsify the transferability of the potential.

Figures

Figures reproduced from arXiv: 2605.05620 by Chen Shen, Izabela Szlufarska, Jiechen Wang, Liang Peng, Masood Mortazavi, Nuohao Liu, Pingfan Wu, Song Xue, Yue Cao, Zongfang Lin.

Figure 1
Figure 1. Figure 1: Overview of the multiscale workflow used in this work. DFT and AIMD calculations first view at source ↗
Figure 2
Figure 2. Figure 2: Supervised validation of the unified NEP. (a) Training loss curve. (b) Parity plot for energies. view at source ↗
Figure 3
Figure 3. Figure 3: Bulk vibrational and thermal validation of the learned potential. Phonon dispersions of (a) view at source ↗
Figure 4
Figure 4. Figure 4: Overall TBC trends from NEMD simulation. (a) Representative steady-state temperature view at source ↗
Figure 5
Figure 5. Figure 5: Mechanistic interpretation of the length-dependent TBC. (a) Representative frequency view at source ↗
Figure 6
Figure 6. Figure 6: Mechanistic interpretation of the temperature-dependent TBC. (a) Frequency-binned view at source ↗
Figure 7
Figure 7. Figure 7: Orientation-dependent interfacial diagnostics. (a) Planar-averaged charge-density-difference view at source ↗
Figure 8
Figure 8. Figure 8: Supplementary bin-resolved cross-interface coupling analysis at view at source ↗
read the original abstract

Ga$_2$O$_3$/SiC heterointegration is attractive for ultra-wide-bandgap power electronics, but interfacial thermal boundary conductance (TBC) remains a major heat-removal bottleneck. Direct experimental access to intrinsic atomistic interfacial transport remains limited, particularly for ideally synthesized materials with defect-free interfacial contact. First-principles simulations are too expensive at relevant length and time scales, while empirical Molecular Dynamics (MD) potentials often lack transferability across oxide and carbide bonding environments. We develop a unified feedforward neural network potential and validate it against density-functional data, bulk phonon dispersions, and anisotropic thermal-conductivity trends in both $\beta$-Ga$_2$O$_3$ and SiC. Nonequilibrium simulations show that TBC decreases with transport length, increases with temperature, and is consistently higher for Ga$_2$O$_3$$(\bar{2}01)$/SiC(0001) than for Ga$_2$O$_3$(100)/SiC(0001). These trends are explained by attenuation of long-mean-free-path carriers, enhanced incoherent and anharmonic interfacial exchange within broadly unchanged spectral channels, and stronger bonding and vibrational coupling at the $(\bar{2}01)$ interface. The results show how a single transferable feedforward neural network potential can enable large-scale transport prediction and physics-grounded mechanistic understanding of thermal boundary conductance. Code for NEP training and simulation workflows is available at the project repository https://github.com/knowhow07/TBC_Ga2O3_SiC.git

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The manuscript develops a unified feedforward neural network potential (NEP) trained on DFT data for β-Ga₂O₃ and SiC. It validates the potential against bulk DFT calculations, phonon dispersions, and anisotropic thermal conductivities in each material separately. Nonequilibrium MD simulations with this potential are then used to predict thermal boundary conductance (TBC) at Ga₂O₃/SiC heterointerfaces, showing that TBC decreases with transport length, increases with temperature, and is higher for the Ga₂O₃(¯201)/SiC(0001) orientation than for Ga₂O₃(100)/SiC(0001). These trends are attributed to attenuation of long-mean-free-path carriers, enhanced incoherent/anharmonic interfacial exchange, and stronger bonding/vibrational coupling at the (¯201) interface. Open code for training and workflows is provided.

Significance. If the bulk-trained NEP transfers reliably to the atomically sharp interface, the work demonstrates how machine-learned potentials can bridge first-principles accuracy with large-scale nonequilibrium transport simulations, yielding mechanistic explanations for TBC trends relevant to ultra-wide-bandgap power electronics. The open repository strengthens reproducibility.

major comments (1)
  1. [Abstract] Abstract (and implied validation section): Validation is reported only against bulk DFT data, phonon dispersions, and anisotropic conductivity in each material separately. No interface-specific checks—such as formation energy of the Ga₂O₃/SiC heterojunction, force errors on atoms near the boundary, or local vibrational density of states—are described. This is load-bearing for the central claims, as the reported TBC values and the attribution to 'stronger bonding and vibrational coupling at the (¯201) interface' depend on accurate cross-material interactions that may not transfer from bulk-only training.
minor comments (2)
  1. [Abstract] Abstract: Notation for the (¯201) orientation is inconsistent between 'Ga₂O₃(¯201)' and 'Ga₂O₃(201)'; standardize throughout.
  2. Consider adding a brief table or figure panel comparing absolute TBC values against any available experimental reports or prior simulations to contextualize the results.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the positive assessment of the work's significance and for highlighting the importance of interface validation. We address the major comment point by point below and will revise the manuscript accordingly to strengthen the transferability claims.

read point-by-point responses
  1. Referee: [Abstract] Abstract (and implied validation section): Validation is reported only against bulk DFT data, phonon dispersions, and anisotropic thermal conductivities in each material separately. No interface-specific checks—such as formation energy of the Ga₂O₃/SiC heterojunction, force errors on atoms near the boundary, or local vibrational density of states—are described. This is load-bearing for the central claims, as the reported TBC values and the attribution to 'stronger bonding and vibrational coupling at the (¯201) interface' depend on accurate cross-material interactions that may not transfer from bulk-only training.

    Authors: We agree that explicit validation of the NEP at the interface is important to support the TBC predictions and mechanistic attributions. The training set includes DFT data on bulk Ga₂O₃ and SiC configurations as well as mixed environments to capture cross-material interactions, and the unified architecture is designed to ensure transferability. However, we did not report interface-specific metrics such as heterojunction formation energies, near-boundary force errors, or local VDOS in the original manuscript. In the revised version, we will add these checks: (i) comparison of interface formation energy to direct DFT calculations on the same supercells, (ii) force error statistics on atoms within 5 Å of the boundary, and (iii) local vibrational density of states extracted from short MD runs at both orientations. These additions will directly address the transferability concern and provide quantitative support for the stronger bonding claim at the (¯201) interface. revision: yes

Circularity Check

0 steps flagged

No circularity; TBC trends are independent MD outputs from externally trained potential

full rationale

The derivation proceeds from external DFT data used to train the feedforward NN potential, followed by validation against bulk phonon dispersions and anisotropic thermal conductivities, then nonequilibrium MD simulations that generate TBC values and trends. These simulation outputs are not algebraically equivalent to the training inputs or fitting procedure; the potential parameters are not redefined in terms of TBC, no self-citations bear the central load, and no uniqueness theorems or ansatzes are smuggled in. The chain remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim rests on the transferability of a neural-network potential fitted to bulk DFT data; no new physical entities are postulated, but the model implicitly assumes that short-range many-body interactions learned from bulk suffice at the interface.

free parameters (1)
  • neural-network weights and biases
    Hundreds of parameters in the feedforward network are optimized against DFT energies and forces for Ga2O3 and SiC configurations.
axioms (2)
  • domain assumption The neural-network potential reproduces the correct phonon spectrum and thermal conductivity of the bulk phases
    Invoked to justify use of the potential for interfacial transport.
  • domain assumption Nonequilibrium molecular dynamics with the potential yields the intrinsic, defect-free TBC
    Standard assumption in classical MD transport studies.

pith-pipeline@v0.9.0 · 5622 in / 1409 out tokens · 32131 ms · 2026-05-08T08:57:52.272936+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

34 extracted references · 34 canonical work pages

  1. [1]

    Perspective: Ga2O3 for ultra-high power rectifiers and MOSFETS.J

    S J Pearton, Fan Ren, Marko Tadjer, and Jihyun Kim. Perspective: Ga2O3 for ultra-high power rectifiers and MOSFETS.J. Appl. Phys., 124(22):220901, December 2018

  2. [2]

    β-Gallium oxide power electronics.APL Mater ., 10(2):029201, February 2022

    Andrew J Green, James Speck, Grace Xing, Peter Moens, Fredrik Allerstam, Krister Gumaelius, Thomas Neyer, Andrea Arias-Purdue, Vivek Mehrotra, Akito Kuramata, Kohei Sasaki, Shinya Watanabe, Kimiyoshi Koshi, John Blevins, Oliver Bierwagen, Sriram Krishnamoorthy, Kevin Leedy, Aaron R Arehart, Adam T Neal, Shin Mou, Steven A Ringel, Avinash Kumar, Ankit Shar...

  3. [3]

    A review of Ga 2O3 materials, processing, and devices.Appl

    S J Pearton, Jiancheng Yang, Patrick H Cary, IV , F Ren, Jihyun Kim, Marko J Tadjer, and Michael A Mastro. A review of Ga 2O3 materials, processing, and devices.Appl. Phys. Rev., 5 (1):011301, March 2018

  4. [4]

    Toward emerging gallium oxide semiconductors: A roadmap.Fundam

    Yuan Yuan, Weibing Hao, Wenxiang Mu, Zhengpeng Wang, Xuanhu Chen, Qi Liu, Guangwei Xu, Chenlu Wang, Hong Zhou, Yanni Zou, Xiaolong Zhao, Zhitai Jia, Jiandong Ye, Jincheng Zhang, Shibing Long, Xutang Tao, Rong Zhang, and Yue Hao. Toward emerging gallium oxide semiconductors: A roadmap.Fundam. Res., 1(6):697–716, November 2021

  5. [5]

    Anisotropic thermal conductivity in single crystal β-gallium oxide.Appl

    Zhi Guo, Amit Verma, Xufei Wu, Fangyuan Sun, Austin Hickman, Takekazu Masui, Akito Kuramata, Masataka Higashiwaki, Debdeep Jena, and Tengfei Luo. Anisotropic thermal conductivity in single crystal β-gallium oxide.Appl. Phys. Lett., 106(11):111909, March 2015

  6. [6]

    Heteroepitaxial growth of β-Ga2O3 films on SiC via molecular beam epitaxy.J

    Neeraj Nepal, D Scott Katzer, Brian P Downey, Virginia D Wheeler, Luke O Nyakiti, David F Storm, Matthew T Hardy, Jaime A Freitas, Eric N Jin, Diego Vaca, Luke Yates, Samuel Graham, Satish Kumar, and David J Meyer. Heteroepitaxial growth of β-Ga2O3 films on SiC via molecular beam epitaxy.J. V ac. Sci. Technol. A, 38(6), December 2020

  7. [7]

    Thermal conductivity of 4H-SiC single crystals.J

    Rusheng Wei, Sheng Song, Kun Yang, Yingxin Cui, Yan Peng, Xiufang Chen, Xiaobo Hu, and Xiangang Xu. Thermal conductivity of 4H-SiC single crystals.J. Appl. Phys., 113(5):053503, February 2013

  8. [8]

    Anisotropic thermal conductivity of 4H and 6H silicon carbide measured using time-domain thermoreflectance.Mater

    Xin Qian, Puqing Jiang, and Ronggui Yang. Anisotropic thermal conductivity of 4H and 6H silicon carbide measured using time-domain thermoreflectance.Mater . Today Phys., 3:70–75, December 2017

  9. [9]

    Improving interfacial thermal conductivity by constructing covalent bond between Ga2O3 and SiC.Nat

    Yi Shen, Xin Qi, Yuan Li, Yunduo Guo, Qimin Huang, Wen Dai, Qilong Yuan, Lin Gu, Chengxi Ding, Wen-Jie Chen, Mingyang Yang, Zhenglin Jia, Cheng-Te Lin, Nan Jiang, Fangyuan 10 Sun, Qingchun Jon Zhang, and Hong-Ping Ma. Improving interfacial thermal conductivity by constructing covalent bond between Ga2O3 and SiC.Nat. Commun., 16(1):10723, November 2025

  10. [10]

    Study of the bonding characteristics at β-Ga2O3(201)/4H-SiC(0001) interfaces from first principles and experiment

    Bei Xu, Jichao Hu, Jiaqi Meng, Xiaomin He, Xi Wang, and Hongbin Pu. Study of the bonding characteristics at β-Ga2O3(201)/4H-SiC(0001) interfaces from first principles and experiment. Crystals (Basel), 13(2):160, January 2023

  11. [11]

    Investigation of the interface electronic characteristics of β-Ga2O3 (1 0 0)/4H-SiC (0 0 0 1).J

    Bei Xu, Jichao Hu, Jiaqi Meng, Xiaomin He, Hongjuan Cheng, Jian Wang, Xi Wang, and Hongbin Pu. Investigation of the interface electronic characteristics of β-Ga2O3 (1 0 0)/4H-SiC (0 0 0 1).J. Cryst. Growth, 601(126951):126951, January 2023

  12. [12]

    Interfacial thermal transport investigation ofβ-Ga2O3/4H-SiC power devices via active-learning- driven neuroevolution potential.Int

    Guangzheng Zhang, Shilin Dong, Lin Guo, Chenghao Li, Xinyu Wang, and Gongming Xin. Interfacial thermal transport investigation ofβ-Ga2O3/4H-SiC power devices via active-learning- driven neuroevolution potential.Int. J. Heat Mass Transf., 263(128710):128710, March 2026

  13. [13]

    Neuroevolution machine learning potentials: Combining high accuracy and low cost in atomistic simulations and application to heat transport.Phys

    Zheyong Fan, Zezhu Zeng, Cunzhi Zhang, Yanzhou Wang, Keke Song, Haikuan Dong, Yue Chen, and Tapio Ala-Nissila. Neuroevolution machine learning potentials: Combining high accuracy and low cost in atomistic simulations and application to heat transport.Phys. Rev. B., 104(10), September 2021

  14. [14]

    Gpumd: A package for constructing accurate machine-learned potentials and performing highly efficient atomistic simulations.J

    Zheyong Fan, Yanzhou Wang, Penghua Ying, Keke Song, Junjie Wang, Yong Wang, Zezhu Zeng, et al. Gpumd: A package for constructing accurate machine-learned potentials and performing highly efficient atomistic simulations.J. Chem. Phys., 157(11), 2022

  15. [15]

    A neuroevolution potential for predicting the thermal conductivity ofα,β, andϵ-Ga 2O3.Appl

    Zhanpeng Sun, Zijun Qi, Kang Liang, Xiang Sun, Zhaofu Zhang, Lijie Li, Qijun Wang, Guoqing Zhang, Gai Wu, and Wei Shen. A neuroevolution potential for predicting the thermal conductivity ofα,β, andϵ-Ga 2O3.Appl. Phys. Lett., 123(19), November 2023

  16. [16]

    High dimensions and heavy tails for natural evolution strategies

    Tom Schaul, Tobias Glasmachers, and Jürgen Schmidhuber. High dimensions and heavy tails for natural evolution strategies. InProceedings of the 13th Annual Conference on Genetic and Evolutionary Computation, pages 845–852. Association for Computing Machinery, 2011

  17. [17]

    Enhanced thermal conductance at interfaces between gold and amorphous silicon and between gold and amorphous silica.Phys

    Julien El Hajj, Christophe Adessi, Michaël de San Feliciano, Gilles Ledoux, and Samy Merabia. Enhanced thermal conductance at interfaces between gold and amorphous silicon and between gold and amorphous silica.Phys. Rev. B., 110(11), September 2024

  18. [18]

    The role of coherence in time delay estimation

    G Clifford Carter. The role of coherence in time delay estimation. InAspects of Signal Processing, pages 251–256. Springer Netherlands, Dordrecht, 1977

  19. [19]

    Microscopic description of thermal-phonon coherence: From coherent transport to diffuse interface scattering in superlattices.Phys

    B Latour, S V olz, and Y Chalopin. Microscopic description of thermal-phonon coherence: From coherent transport to diffuse interface scattering in superlattices.Phys. Rev. B Condens. Matter Mater . Phys., 90(1), July 2014

  20. [20]

    Development of a neuroevolution machine learning potential of Pd-Cu-Ni-P alloys.Mater

    Rui Zhao, Shucheng Wang, Zhuangzhuang Kong, Yunlei Xu, Kuan Fu, Ping Peng, and Cuilan Wu. Development of a neuroevolution machine learning potential of Pd-Cu-Ni-P alloys.Mater . Des., 231(112012):112012, July 2023

  21. [21]

    Insight into the effect of force error on the thermal conductivity from machine-learned potentials.Mater

    Wenjiang Zhou, Nianjie Liang, Xiguang Wu, Shiyun Xiong, Zheyong Fan, and Bai Song. Insight into the effect of force error on the thermal conductivity from machine-learned potentials.Mater . Today Phys., 50(101638):101638, January 2025

  22. [22]

    Correcting force error-induced underestimation of lattice thermal conductivity in machine learning molecular dynamics.J

    Xiguang Wu, Wenjiang Zhou, Haikuan Dong, Penghua Ying, Yanzhou Wang, Bai Song, Zheyong Fan, and Shiyun Xiong. Correcting force error-induced underestimation of lattice thermal conductivity in machine learning molecular dynamics.J. Chem. Phys., 161(1), July 2024

  23. [23]

    Three-dimensional anisotropic thermal conductivity tensor of single crystalline β-Ga2O3.Appl

    Puqing Jiang, Xin Qian, Xiaobo Li, and Ronggui Yang. Three-dimensional anisotropic thermal conductivity tensor of single crystalline β-Ga2O3.Appl. Phys. Lett., 113(23):232105, December 2018

  24. [24]

    Molecular dynamics simulations of heat transport using machine-learned potentials: A mini-review and tutorial on GPUMD with neuroevolution potentials.J

    Haikuan Dong, Yongbo Shi, Penghua Ying, Ke Xu, Ting Liang, Yanzhou Wang, Zezhu Zeng, Xin Wu, Wenjiang Zhou, Shiyun Xiong, Shunda Chen, and Zheyong Fan. Molecular dynamics simulations of heat transport using machine-learned potentials: A mini-review and tutorial on GPUMD with neuroevolution potentials.J. Appl. Phys., 135(16), April 2024. 11

  25. [25]

    A universal graph deep learning interatomic potential for the periodic table.Nat

    Chi Chen and Shyue Ping Ong. A universal graph deep learning interatomic potential for the periodic table.Nat. Comput. Sci., 2(11):718–728, 2022

  26. [26]

    Ab-initio simulations of materials using vasp: density-functional theory and beyond.J

    Jürgen Hafner. Ab-initio simulations of materials using vasp: density-functional theory and beyond.J. Comput. Chem., 29(13):2044–2078, 2008

  27. [27]

    Scuseria

    Matthias Ernzerhof and Gustavo E. Scuseria. Assessment of the perdew–burke–ernzerhof exchange-correlation functional.J. Chem. Phys., 110(11):5029–5036, 1999

  28. [28]

    A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in ’t Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. Lammps - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales.Comput. Phys. Commun., 2...

  29. [29]

    Melville S. Green. Markoff random processes and the statistical mechanics of time-dependent phenomena. ii. irreversible processes in fluids.J. Chem. Phys., 22(3):398–413, 1954

  30. [30]

    Statistical-mechanical theory of irreversible processes

    Ryogo Kubo, Mario Yokota, and Sadao Nakajima. Statistical-mechanical theory of irreversible processes. ii. response to thermal disturbance.J. Phys. Soc. Jpn., 12(11):1203–1211, 1957

  31. [31]

    Visualization and analysis of atomistic simulation data with OVITO–the open visualization tool.Model

    Alexander Stukowski. Visualization and analysis of atomistic simulation data with OVITO–the open visualization tool.Model. Simul. Mat. Sci. Eng., 18(1):015012, January 2010

  32. [32]

    VESTA 3 for three-dimensional visualization of crystal, volumetric and morphology data.J

    Koichi Momma and Fujio Izumi. VESTA 3 for three-dimensional visualization of crystal, volumetric and morphology data.J. Appl. Crystallogr ., 44(6):1272–1276, December 2011

  33. [33]

    A fast and robust algorithm for bader decomposition of charge density.Comput

    Graeme Henkelman, Andri Arnaldsson, and Hannes Jónsson. A fast and robust algorithm for bader decomposition of charge density.Comput. Mater . Sci., 36(3):354–360, June 2006

  34. [34]

    1 Nstr NstrX n=1 U NEP n −U DFT n 2 #1/2 +λ F

    Atsushi Togo, Laurent Chaput, Terumasa Tadano, and Isao Tanaka. Implementation strategies in phonopy and phono3py.J. Phys. Condens. Matter, 35(35):353001, 2023. A List of Abbreviations This section summarizes the abbreviations used throughout this work, together with their full names and brief explanations. •AIMD – Ab Initio Molecular Dynamics A simulatio...