Thermal Conductivity Modeling using Machine Learning Potentials: Application to Crystalline and Amorphous Silicon
Pith reviewed 2026-05-24 18:31 UTC · model grok-4.3
The pith
Machine learning potentials derived from density functional theory calculations enable thermal conductivity predictions for both crystalline and amorphous silicon via molecular dynamics.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The machine learning based interatomic potential is derived from density functional theory calculations by stochastically sampling the potential energy surface in the configurational space. The thermal conductivities of both amorphous and crystalline silicon are then calculated using equilibrium molecular dynamics, which agree well with experimental measurements. This work documents the procedure for training the machine-learning based potentials for modeling thermal conductivity, and demonstrates that machine-learning based potential can be a promising tool for modeling thermal conductivity of both crystalline and amorphous materials with strong disorder.
What carries the argument
Machine-learning interatomic potential trained by stochastic sampling of the DFT potential energy surface, used directly in equilibrium molecular dynamics to obtain thermal conductivity.
If this is right
- Thermal conductivity modeling becomes possible for materials with strong disorder where the phonon quasiparticle model breaks down.
- Atomistic simulations of thermal transport reach length and time scales far beyond direct first-principles calculations.
- The same stochastic sampling and training procedure supplies potentials for other crystalline and amorphous solids.
- Thermal conductivity can be obtained without solving the Boltzmann transport equation explicitly.
Where Pith is reading between the lines
- The method could be applied to other elemental or compound semiconductors to test whether the same training protocol works across different bonding types.
- It opens the possibility of computing thermal conductivity in nanostructures or at interfaces that contain both ordered and disordered regions.
- Similar machine-learning potentials might allow direct simulation of heat flow in glasses or polymers at scales where explicit phonon calculations are impossible.
Load-bearing premise
The machine-learning potential obtained by stochastically sampling the DFT potential energy surface in configurational space is sufficiently accurate and transferable to capture the atomic dynamics that determine thermal transport in both ordered and disordered silicon structures.
What would settle it
Equilibrium molecular dynamics runs with the trained potential that produce thermal conductivity values for amorphous silicon differing by more than 20 percent from the accepted experimental range would falsify the claim of sufficient accuracy and transferability.
Figures
read the original abstract
First-principles based modeling on phonon dynamics and transport using density functional theory and Boltzmann transport equation has proven powerful in predicting thermal conductivity of crystalline materials, but it remains unfeasible for modeling complex crystals and disordered solids due to the prohibitive computational cost to capture the disordered structure, especially when the quasiparticle "phonon" model breaks down. Recently, machine-learning regression algorithms show great promises for building high-accuracy potential fields for atomistic modeling with length and time scales far beyond those achievable by first-principles calculations. In this work, using both crystalline and amorphous silicon as examples, we develop machine learning based potential fields for predicting thermal conductivity. The machine learning based interatomic potential is derived from density functional theory calculations by stochastically sampling the potential energy surface in the configurational space. The thermal conductivities of both amorphous and crystalline silicon are then calculated using equilibrium molecular dynamics, which agree well with experimental measurements. This work documents the procedure for training the machine-learning based potentials for modeling thermal conductivity, and demonstrates that machine-learning based potential can be a promising tool for modeling thermal conductivity of both crystalline and amorphous materials with strong disorder.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops machine learning interatomic potentials trained on DFT data through stochastic sampling of the configurational space for silicon. These potentials are employed in equilibrium molecular dynamics simulations to calculate thermal conductivities for both crystalline and amorphous silicon, reporting agreement with experimental measurements. The work outlines the training procedure and positions ML potentials as a tool for modeling thermal transport in materials with strong disorder where phonon-based methods are limited.
Significance. If the results hold, the significance lies in demonstrating a scalable method for computing thermal conductivity in disordered systems using ML potentials derived from first-principles data. The paper provides training details, force-error metrics on held-out configurations, and comparisons of computed conductivities with experiment, including system-size checks. This approach is independent of the target observable as the potential is fitted to energies and forces only, not conductivity. It offers a promising alternative to direct DFT for complex materials.
minor comments (2)
- [Abstract] Abstract: the statement that computed conductivities 'agree well with experimental measurements' lacks quantitative metrics, error bars, or specific values; adding these would strengthen the claim without altering the central result.
- [Methods] The description of the stochastic sampling procedure for the PES would benefit from explicit mention of the number of sampled configurations and the train/test split used for the force-error metrics on held-out data.
Simulated Author's Rebuttal
We thank the referee for their careful reading and positive evaluation of our work. The recommendation for minor revision is noted, and we appreciate the recognition of the approach's potential for disordered systems. Since no specific major comments were provided in the report, we have no points requiring direct response or revision at this stage.
Circularity Check
No significant circularity
full rationale
The paper trains an ML interatomic potential solely on DFT energies and forces obtained by stochastic sampling of the configurational space. Thermal conductivity is subsequently obtained from equilibrium MD trajectories driven by that potential; conductivity itself is never an input to the fit. No equations, self-citations, or uniqueness claims reduce the final result to the training data by construction. The reported agreement with experiment therefore constitutes an independent test.
Axiom & Free-Parameter Ledger
free parameters (1)
- ML regression hyperparameters
axioms (1)
- domain assumption Stochastic sampling of the potential energy surface yields a training set representative of configurations relevant to thermal transport in both crystalline and amorphous silicon.
Reference graph
Works this paper leans on
-
[1]
W. Li, N. Mingo, L. Lindsay, D. A. Broido, D. A. Stewart and N. A. Katcho, Phys. Rev. B 85 (19) (2012)
work page 2012
-
[2]
A. Ward, D. A. Broido, D. A. Stewart and G. Deinzer, Phys. Rev. B 80 (12) (2009)
work page 2009
- [3]
-
[4]
D. A. Broido, M. Malorny, G. Birner, N. Mingo and D. A. Stewart, Appl. Phys. Lett. 91 (23), 231922 (2007)
work page 2007
- [5]
-
[6]
L. Lindsay, D. A. Broido and T. L. Reinecke, Phys. Rev. Lett. 109 (9), 095901 (2012)
work page 2012
-
[7]
G. Fugallo, A. Cepellotti, L. Paulatto, M. Lazzeri, N. Marzari and F. Mauri, Nano Lett 14 (11), 6109-6114 (2014)
work page 2014
- [8]
- [9]
- [10]
-
[11]
H. Xie, M. Hu and H. Bao, Appl. Phys. Lett. 104 (13), 131906 (2014)
work page 2014
- [12]
-
[13]
G. Qin, Q. B. Yan, Z. Qin, S. Y. Yue, M. Hu and G. Su, Phys. Chem. Chem. Phys. 17 (7), 4854-4858 (2015)
work page 2015
-
[14]
B. Sun, X. Gu, Q. Zeng, X. Huang, Y. Yan, Z. Liu, R. Yang and Y. K. Koh, Adv. Mater. 29 (3), 1603297 (2017)
work page 2017
-
[15]
J. Zhu, H. Park, J.-Y. Chen, X. Gu, H. Zhang, S. Karthikeyan, N. Wendel, S. A. Campbell, M. Dawber, X. Du, M. Li, J.-P. Wang, R. Yang and X. Wang, Advanced Electronic Materials 2 (5), 1600040 (2016)
work page 2016
-
[16]
X. Gu, B. Li and R. Yang, J. Appl. Phys. 119 (8), 085106 (2016)
work page 2016
- [17]
- [18]
-
[19]
X. Qian, X. Gu and R. Yang, J. Phys. Chem. C (2015)
work page 2015
-
[20]
X. Qian, X. Gu and R. Yang, Nano Energy 41, 394-407 (2017)
work page 2017
-
[21]
J. Yang, X. Qian, W. Pan, R. Yang, Z. Li, Y. Han, M. Zhao, M. Huang and C. Wan, Adv. Mater. 31 (24), 1808222 (2019)
work page 2019
-
[22]
H. R. Seyf, L. Yates, T. L. Bougher, S. Graham, B. A. Cola, T. Detchprohm, M.-H. Ji, J. Kim, R. Dupuis, W. Lv and A. Henry, npj Computational Materials 3, 49 (2017). 16
work page 2017
-
[23]
A. Rohskopf, H. R. Seyf, K. Gordiz, T. Tadano and A. Henry, npj Computational Materials 3 (1) (2017)
work page 2017
-
[24]
P. C. Howell, The Journal of chemical physics 137 (22), 224111 (2012)
work page 2012
- [25]
- [26]
-
[27]
A. P. Bartok, M. C. Payne, R. Kondor and G. Csanyi, Phys. Rev. Lett. 104 (13), 136403 (2010)
work page 2010
-
[28]
A. P. Bartok, J. Kermode, N. Bernstein and G. Csanyi, Physical Review X 8 (4) (2018)
work page 2018
-
[29]
V. L. Deringer and G. Csányi, Phys. Rev. B 95 (9), 094203 (2017)
work page 2017
-
[30]
V. L. Deringer, N. Bernstein, A. P. Bartok, M. J. Cliffe, R. N. Kerber, L. E. Marbella, C. P. Grey, S. R. Elliott and G. Csanyi, J Phys Chem Lett 9 (11), 2879-2885 (2018)
work page 2018
- [31]
- [32]
-
[33]
R. M. Balabin and E. I. Lomakina, Phys. Chem. Chem. Phys. 13 (24), 11710-11718 (2011)
work page 2011
-
[34]
A. P. Thompson, L. P. Swiler, C. R. Trott, S. M. Foiles and G. J. Tucker, Journal of Computational Physics 285, 316-330 (2015)
work page 2015
-
[35]
A. P. Bartók and G. Csányi, Int. J. Quantum Chem 115 (16), 1051-1057 (2015)
work page 2015
- [36]
-
[37]
C. E. Rasmussen and C. K. I. Williams, Gaussian Processes for Machine Learning. (MIT Press, 2006)
work page 2006
- [38]
- [39]
- [40]
- [41]
-
[42]
J. P. Perdew, K. Burke and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996)
work page 1996
-
[43]
W. Li, J. Carrete, N. A. Katcho and N. Mingo, Comput. Phys. Commun. 185 (6), 1747-1758 (2014)
work page 2014
-
[44]
D. Sholl and J. A. Steckel, Density functional theory: a practical introduction. (John Wiley & Sons, 2011)
work page 2011
-
[45]
A. P. Bartók, R. Kondor and G. Csányi, Phys. Rev. B 87 (18), 184115 (2013)
work page 2013
-
[46]
J. Schreiter, D. Nguyen-Tuong and M. Toussaint, Neurocomputing 192, 29-37 (2016)
work page 2016
-
[47]
W. J. Szlachta, A. P. Bartók and G. Csányi, Phys. Rev. B 90 (10), 104108 (2014)
work page 2014
-
[48]
J. Zi, K. Zhang and X. Xie, Phys. Rev. B 41 (18), 12915-12918 (1990)
work page 1990
-
[49]
A. Pedersen, L. Pizzagalli and H. Jónsson, New Journal of Physics 19 (6), 063018 (2017)
work page 2017
-
[50]
Plimpton, Journal of Computational Physics 117, 1-19 (1995)
S. Plimpton, Journal of Computational Physics 117, 1-19 (1995)
work page 1995
-
[51]
M. Puligheddu, Y. Xia, M. K. Y. Chan and G. Galli, arXiv:1902.08260 (2019)
-
[52]
S. G. Volz and G. Chen, Phys. Rev. B 61 (4), 2651 (2000)
work page 2000
-
[53]
B. L. Zink, R. Pietri and F. Hellman, Phys. Rev. Lett. 96 (5), 055902 (2006)
work page 2006
-
[54]
K. T. Regner, D. P. Sellan, Z. Su, C. H. Amon, A. J. McGaughey and J. A. Malen, Nature communications 4, 1640 (2013)
work page 2013
-
[55]
D. G. Cahill, M. Katiyar and J. R. Abelson, Phys Rev B Condens Matter 50 (9), 6077-6081 (1994)
work page 1994
-
[56]
M. C. Wingert, J. Zheng, S. Kwon and R. Chen, Semicond. Sci. Technol. 31 (11), 113003 (2016)
work page 2016
-
[57]
J. M. Larkin and A. J. H. McGaughey, Phys. Rev. B 89, 144303 (2014)
work page 2014
- [58]
- [59]
-
[60]
C. J. Glassbrenner and G. A. Slack, Phys. Rev. 134 (4A), A1058-A1069 (1964)
work page 1964
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.