pith. sign in

arxiv: 2505.18127 · v3 · submitted 2025-05-23 · ❄️ cond-mat.mtrl-sci · cond-mat.mes-hall

Tuning Thermal Conductivity and Electron-Phonon Interactions in Carbon and Boron Nitride Moir\'e Diamanes via Twist Angle Manipulation

Pith reviewed 2026-05-19 13:04 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci cond-mat.mes-hall
keywords moiré diamanestwist anglelattice thermal conductivitystructural disorderband gap renormalizationcarbonboron nitrideanharmonic effects
0
0 comments X

The pith

Increasing the twist angle in carbon and boron nitride moiré diamanes reduces their in-plane lattice thermal conductivity by 4.5 to 9 times due to increased structural disorder.

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

The paper investigates the impact of interlayer twist angle on lattice thermal conductivity and band gap renormalization in boron nitride and carbon moiré diamanes. Using moment tensor potentials and solving the Boltzmann transport equation along with the Green-Kubo formula, the authors calculate thermal properties. They observe that higher twist angles lead to greater structural disorder, causing a substantial drop in in-plane thermal conductivity. This disorder also results in larger band gap changes from nuclear motion, with significant renormalization from quantum effects due to high-frequency phonons involving hydrogen atoms. Such findings offer a way to tune thermal and electronic properties for applications in thermoelectrics, microelectronics, and optoelectronics.

Core claim

We have shown that increasing the twist angle in moiré diamanes of carbon and boron nitride leads to a reduction in in-plane lattice thermal conductivity by a factor of 4.5 to 9, driven by the growth of structural disorder. This disorder also enhances band gap renormalization induced by classical nuclei motion, while high phonon frequencies from surface hydrogen bonds cause notable renormalization when quantum nuclear effects are considered. The calculations reveal 20-40% differences between Green-Kubo and Boltzmann transport equation methods, highlighting the role of high-order anharmonic contributions.

What carries the argument

Twist angle manipulation that induces structural disorder in moiré diamanes, thereby scattering phonons and modifying electron-phonon interactions.

Load-bearing premise

The moment tensor potentials accurately reproduce the interatomic forces and anharmonic interactions for the range of twist angles and disorder levels in the moiré diamanes.

What would settle it

Measuring the in-plane thermal conductivity of moiré diamane samples prepared with increasing twist angles from 0 to 30 degrees using techniques like time-domain thermoreflectance would test if the conductivity drops by the predicted factors of 4.5 to 9.

Figures

Figures reproduced from arXiv: 2505.18127 by Alexander Kvashnin, Alexander Shapeev, Leonid Chernozatonskii, Mikhail Polovinkin, Nikita Rybin, Rustam Arabov, Victor Demin.

Figure 1
Figure 1. Figure 1: Top and side views of the atomic structure of considered hydrogenated BN and [PITH_FULL_IMAGE:figures/full_fig_p014_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Validation of MTP for BNnAB. Comparison of: (a) - forces obtained with DFT [PITH_FULL_IMAGE:figures/full_fig_p015_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Validation of MTP for DnAB. Comparison of: (a) - forces obtained with DFT [PITH_FULL_IMAGE:figures/full_fig_p016_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Phonon band structures for the BNnθ Moiré lattices ((a) - BNnAB, (b) - BNn21.8, (c) - BNn27.8) calculated with MTP and DFT. For convenience of visual per￾ception, the frequency ranges are shown partially. The full graphs with uncut y axes can be found in Fig. S5. and Dnθ Moiré lattices studied here, the electronic contribution is negligible, making the lattice component dominant. To obtain the lattice ther… view at source ↗
Figure 5
Figure 5. Figure 5: Phonon band structures for the Dnθ Moiré lattices ((a) - DnAB, (b) - Dn21.8, (c) - Dn27.8) calculated with MTP and DFT. For convenience of visual perception, the frequency ranges are shown partially. The full graphs with uncut y axes can be found in Fig. S6. of phonon heat capacities, group velocities, and lifetimes are required. All these values are obtained from the phonon band structures calculated with… view at source ↗
Figure 6
Figure 6. Figure 6: Phonon heat capacities in BNnθ (a) and Dnθ (b) structures. lead to lower thermal conductivity. As shown in [PITH_FULL_IMAGE:figures/full_fig_p018_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Phonon group velocities in BNnθ (a) and Dnθ (b) structures. bonds start deviating from 90◦ . The similar situation occurs in BNnθ, as shown in Fig. S3 in the Supplementary information. This fact also proves that the degree of disorder grows with the increase of the twist angle in the considered Moiré lattices. The growing degree of disorder causes the increase of phonon scattering. This leads to the decrea… view at source ↗
Figure 8
Figure 8. Figure 8: Phonon lifetimes for BNnθ (a) and Dnθ (b) structures at T = 300 K. (a) (b) (c) [PITH_FULL_IMAGE:figures/full_fig_p020_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: C-C bond length distributions for the hydrogenated graphene Moiré lattices [PITH_FULL_IMAGE:figures/full_fig_p020_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Distributions of the angle β between the planes and interplanar bonds in the Dnθ lattices ((a) - DnAB, (b) - Dn21.8, (c) - Dn27.8) [PITH_FULL_IMAGE:figures/full_fig_p021_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Temperature dependence of the lattice thermal conductivity in BNn [PITH_FULL_IMAGE:figures/full_fig_p022_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Comparison of the LTC values obtained with BTE-based (BTE) approach and [PITH_FULL_IMAGE:figures/full_fig_p023_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Band gap renormalization in BNnθ (a) and Dnθ (b) Moiré lattices (BNnθ/Dnθ(NVT) - the distorted samples for calculating band gap values were obtained with classical NVT MD, fitted with line, BNnθ/Dnθ(NpT) - the distorted samples were generated by means of classical NpT MD, fitted with line. In the case of quantum nuclei, the values of ZPR in all Moiré lattices are high. These values are comparable with the… view at source ↗
Figure 14
Figure 14. Figure 14: Band gap renormalization in BNnθ (a) and Dnθ (b) Moiré lattices (BNnθ/Dnθ(c) - the distorted samples for calculating band gap values were obtained with classical NVT MD, fitted with line, BNnθ/Dnθ(q) - the distorted samples were generated by means of quantum harmonic sampling, as implemented in [55]), fitted with modified Varshni’s expression [59]. with higher maximal vibrational frequencies (the phonon f… view at source ↗
read the original abstract

We have investigated the effect of interlayer twist angle on lattice thermal conductivity (LTC) and band gap renormalization in boron nitride and carbon Moir\'e diamanes. Moment tensor potentials were used for calculating energies and forces of interatomic interactions. The methods based on the solution of Boltzmann transport equation (BTE) for phonons and the GreenKubo (GK) formula were utilized to calculate LTC. The 20-40 % difference in LTC values obtained with GK and BTE-based methods showed the importance of high-order anharmonic contributions to LTC. Significant reduction (by 4.5 - 9 times) of the in-plane LTC with the twist angle increase caused by the growth of structural disorder was observed in the Moir\'e diamanes. This growth of disorder also leads to higher band gap renormalization (induced by classical nuclei motion) in the structures with higher twist angles. Significant band gap renormalization values obtained considering the quantum nuclear effects are caused by the high phonon frequencies related to the bonds with hydrogen atoms on the Moir\'e diamanes surfaces. Understanding of the twist angle effect on LTC and electron-phonon coupling in the Moir\'e diamanes provides a fundamental basis for manipulating their thermal and electronic properties, making these materials promising for thermoelectrics, microelectronics and optoelectronics.

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 investigates the effect of interlayer twist angle on lattice thermal conductivity (LTC) and band gap renormalization in boron nitride and carbon Moiré diamanes. Moment tensor potentials (MTPs) are employed to compute interatomic energies and forces, with LTC evaluated via both the Boltzmann transport equation (BTE) for phonons and the Green-Kubo (GK) formula. The authors report a 20-40% discrepancy between BTE and GK results, indicating significant high-order anharmonic contributions, and a 4.5-9 times reduction in in-plane LTC with increasing twist angle, attributed to growing structural disorder. Higher twist angles are also linked to increased band gap renormalization, with notable values arising from quantum nuclear effects due to high-frequency phonons associated with surface hydrogen bonds.

Significance. If the central numerical results hold after validation, the work would offer useful quantitative guidance on twist-angle engineering of thermal transport and electron-phonon coupling in these layered materials, with potential relevance to thermoelectrics and nanoelectronics. The dual use of BTE and GK methods, together with explicit consideration of quantum nuclear effects on the band gap, represents a methodological strength. The reported LTC reduction factors and anharmonicity discrepancies are concrete and falsifiable, which strengthens the paper's contribution if the underlying potentials are shown to be reliable across the explored configurations.

major comments (1)
  1. [Methods] Methods section (description of MTP training and usage): No quantitative benchmarks are provided for the accuracy of the moment tensor potentials in reproducing anharmonic forces or energies for the full range of twist angles and associated structural disorder. Force-error statistics, phonon-dispersion comparisons to DFT reference calculations, or transferability tests at the largest twist angles are absent. Because the headline claim of a 4.5–9× LTC drop is attributed to disorder-induced scattering and rests entirely on MTP-derived forces fed into both BTE and GK calculations, the lack of such validation is load-bearing for the central attribution.
minor comments (2)
  1. [Abstract] The abstract states a '20-40 % difference' between GK and BTE LTC values but does not specify which method yields the higher values or provide error bars; adding this information would improve clarity.
  2. [Figures and Results] Figure captions and text should explicitly define the in-plane versus out-of-plane LTC components and the precise definition of 'structural disorder' (e.g., atomic displacement variance or bond-length distribution) used to correlate with the LTC reduction.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their thorough review and insightful comments on our manuscript. We address the major comment point by point below and outline the revisions we plan to make.

read point-by-point responses
  1. Referee: [Methods] Methods section (description of MTP training and usage): No quantitative benchmarks are provided for the accuracy of the moment tensor potentials in reproducing anharmonic forces or energies for the full range of twist angles and associated structural disorder. Force-error statistics, phonon-dispersion comparisons to DFT reference calculations, or transferability tests at the largest twist angles are absent. Because the headline claim of a 4.5–9× LTC drop is attributed to disorder-induced scattering and rests entirely on MTP-derived forces fed into both BTE and GK calculations, the lack of such validation is load-bearing for the central attribution.

    Authors: We agree with the referee that quantitative validation of the MTPs is crucial for supporting the central claims regarding the LTC reduction. Although the MTPs were trained on DFT data for a range of configurations including twisted structures, explicit benchmarks were not reported in the original submission. In the revised manuscript, we will add detailed information on the MTP training, including root-mean-square errors for forces and energies on validation sets that cover the full range of twist angles. Additionally, we will include comparisons of phonon dispersions calculated with MTPs versus DFT for representative low-twist and high-twist configurations to demonstrate accuracy and transferability. These additions will directly address the concerns about the reliability of the potentials for disordered systems. revision: yes

Circularity Check

0 steps flagged

No significant circularity; results from direct simulation

full rationale

The paper computes LTC via MTP-derived forces and energies fed into standard BTE and Green-Kubo solvers applied to explicitly constructed moiré structures at varying twist angles. These numerical outputs (including the reported 4.5–9× reduction) are not defined in terms of the target LTC values, nor obtained by fitting parameters to subsets of the same LTC data. No self-citation chains, uniqueness theorems, or ansatzes are invoked to justify the core computational pipeline. The derivation remains independent of its own results and is therefore self-contained.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The study rests on standard assumptions about the transferability of moment tensor potentials and the validity of the Boltzmann transport and Green-Kubo frameworks for these disordered layered systems; no new entities are postulated and no free parameters are explicitly fitted beyond the choice of twist angles.

axioms (2)
  • domain assumption Moment tensor potentials trained on reference configurations accurately capture forces and anharmonic interactions in twisted diamane structures
    Invoked when using the potentials to compute energies and forces for all twist angles
  • domain assumption The Boltzmann transport equation and Green-Kubo formula remain applicable despite increasing structural disorder from twist
    Used to obtain LTC values whose 20-40% difference is interpreted as evidence of high-order anharmonicity

pith-pipeline@v0.9.0 · 5807 in / 1573 out tokens · 37894 ms · 2026-05-19T13:04:51.797113+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

59 extracted references · 59 canonical work pages

  1. [1]

    S. Carr, D. Massatt, S. Fang, P. Cazeaux, M. Luskin, E. Kaxiras, Twistronics: Manipulating the electronic properties of two-dimensional layered structures through their twist angle, Physical Review B 95 (7) (2017) 075420

  2. [2]

    Y. Cao, V. Fatemi, S. Fang, K. Watanabe, T. Taniguchi, E. Kaxi- ras, P. Jarillo-Herrero, Unconventional superconductivity in magic-angle graphene superlattices, Nature 556 (7699) (2018) 43–50

  3. [3]

    Bistritzer, A

    R. Bistritzer, A. H. MacDonald, Moiré bands in twisted double-layer graphene, Proceedings of the National Academy of Sciences 108 (30) (2011) 12233—12237

  4. [4]

    D. Wu, Y. Pan, T. Min, Twistronics in graphene, from transfer assembly to epitaxy, Applied Sciences 10 (14) (2020) 4690

  5. [5]

    Y. Hou, J. Zhou, M. Xue, M. Yu, Y. Han, Z. Zhang, Y. Lu, Strain engineering of twisted bilayer graphene: The rise of strain-twistronics, Small (2024) 2311185

  6. [6]

    A. I. Cocemasov, D. L. Nika, A. A. Balandin, Phonons in twisted bilayer graphene, Physical Review B 88 (3) (2013) 035428

  7. [7]

    Chowdhury, V

    S. Chowdhury, V. A. Demin, L. A. Chernozatonskii, A. G. Kvashnin, Ultra-low thermal conductivity of Moiré diamanes, Membranes 12 (10) (2022) 925

  8. [8]

    Eriksson, E

    F. Eriksson, E. Fransson, C. Linderalv, Z. Fan, P. Erhart, Tuning the through-plane lattice thermal conductivity in van der Waals structures throughrotational(dis)ordering, ACSnano17(24)(2023)25565–25574. 26

  9. [9]

    Manunza, R

    L. Manunza, R. Dettori, A. Cappai, C. Melis, Thermal conductivity of graphene Moiré superlattices at small twist angles: An approach-to- equilibrium molecular dynamics and Boltzmann transport study, C 11 (3) (2025) 46

  10. [10]

    K. S. Novoselov, A. Mishchenko, A. Carvalho, A. Castro Neto, 2D ma- terials and van der Waals heterostructures, Science 353 (6298) (2016) aac9439

  11. [11]

    Zacharias, M

    M. Zacharias, M. Scheffler, C. Carbogno, Fully anharmonic nonpertur- bative theory of vibronically renormalized electronic band structures, Phys. Rev. B 102 (2020) 045126

  12. [12]

    Zheng, D

    J. Zheng, D. Shi, Y. Yang, C. Lin, H. Huang, R. Guo, B. Huang, Anharmonicity-induced phonon hardening and phonon transport en- hancement in crystalline perovskite BaZrO3, Phys. Rev. B 105 (2022) 224303

  13. [13]

    Zhu, D.A

    X. Zhu, D.A. Egger, Effect of Overdamped Phonons on the Fundamental Band Gap of Perovskites, Phys. Rev. Lett. 134 (2025) 016403

  14. [14]

    Seidl, X

    S.A. Seidl, X. Zhu, G. Reuveni, S. Aharon, C. Gehrmann, S. Caicedo- Dávila, O. Yaffe, D.A. Egger, Anharmonic fluctuations govern the band gap of halide perovskites, Phys. Rev. Materials 7 (2023) L092401

  15. [15]

    X. Yuan, Y. Zhao, Y. Sun, J. Ni, Z. Dai, Influence of quartic anhar- monicity on lattice dynamics and thermal transport properties of 16 antifluorite structures, Phys. Rev. B 110 (2024) 014304

  16. [16]

    L. A. Chernozatonskii, A. I. Kochaev, K. P. Katin, M. M. Maslov, Moiré M-nitridanes (M= Al, B, Ga) analogues of carbon diamanes:ab initio investigation of atomic structures, electronic, and mechanical properties, ACS Applied Electronic Materials 5 (10) (2023) 5677–5686

  17. [17]

    L. A. Chernozatonskii, V. A. Demin, D. G. Kvashnin, Fully hydro- genated and fluorinated bigraphenes–diamanes: Theoretical and exper- imental studies, C 7 (1) (2021) 17

  18. [18]

    Mortazavi, X

    B. Mortazavi, X. Zhuang, T. Rabczuk, A. V. Shapeev, Atomistic model- ingofthemechanicalproperties: theriseofmachinelearninginteratomic potentials, Materials Horizons 10 (6) (2023) 1956–1968. 27

  19. [19]

    Ladygin, P

    V. Ladygin, P. Y. Korotaev, A. Yanilkin, A. Shapeev, Lattice dynam- ics simulation using machine learning interatomic potentials, Computa- tional Materials Science 172 (2020) 109333

  20. [20]

    Gubaev, E

    K. Gubaev, E. V. Podryabinkin, G. L. Hart, A. V. Shapeev, Accelerat- ing high-throughput searches for new alloys with active learning of inter- atomic potentials, Computational Materials Science 156 (2019) 148–156

  21. [21]

    V. L. Deringer, G. Csányi, Machine learning based interatomic potential for amorphous carbon, Physical Review B 95 (9) (2017) 094203

  22. [22]

    Y. Zuo, C. Chen, X. Li, Z. Deng, Y. Chen, J. Behler, G. Csányi, A. V. Shapeev, A. P. Thompson, M. A. Wood, et al., Performance and cost assessment of machine learning interatomic potentials, The Journal of Physical Chemistry A 124 (4) (2020) 731–745

  23. [23]

    I. S. Novikov, K. Gubaev, E. V. Podryabinkin, A. V. Shapeev, The MLIP package: moment tensor potentials with MPI and active learning, Machine Learning: Science and Technology 2 (2) (2020) 025002

  24. [24]

    Mortazavi, I

    B. Mortazavi, I. S. Novikov, E. V. Podryabinkin, S. Roche, T. Rabczuk, A. V. Shapeev, X. Zhuang, Exploring phononic properties of two- dimensional materials using machine learning interatomic potentials, Applied Materials Today 20 (2020) 100685

  25. [25]

    E. V. Podryabinkin, E. V. Tikhonov, A. V. Shapeev, A. R. Oganov, Ac- celerating crystal structure prediction by machine-learning interatomic potentials with active learning, Physical Review B 99 (6) (2019) 064114

  26. [26]

    Korotaev, I

    P. Korotaev, I. Novoselov, A. Yanilkin, A. Shapeev, Accessing thermal conductivity of complex compounds by machine learning interatomic potentials, Physical Review B 100 (14) (2019) 144308

  27. [27]

    Ghosal, S

    S. Ghosal, S. Chowdhury, D. Jana, Electronic and thermal transport in novel carbon-based bilayer with tetragonal rings: a combined study us- ing first-principles and machine learning approach, Physical Chemistry Chemical Physics 23 (27) (2021) 14608–14616

  28. [28]

    Mortazavi, Ultrahigh thermal conductivity and strength in direct- gap semiconducting graphene-like BC6N: A first-principles and classical investigation, Carbon 182 (2021) 373–383

    B. Mortazavi, Ultrahigh thermal conductivity and strength in direct- gap semiconducting graphene-like BC6N: A first-principles and classical investigation, Carbon 182 (2021) 373–383. 28

  29. [29]

    Kresse, J

    G. Kresse, J. F¨urthmuller, Efficient iterative schemes forab initiototal- energy calculations using a plane-wave basis set, Physical Review B 54 (16) (1996) 11169–11186

  30. [30]

    P. E. Bl¨ochl, Projector augmented-wave method, Physical Review B 50 (24) (1994) 17953–17979

  31. [31]

    J. P. Perdew, K. Burke, M. Ernzerhof, Generalized gradient approxima- tion made simple, Physical Review Letters 77 (18) (1996) 3865–3868

  32. [32]

    S.Grimme, J.Antony, S.Ehrlich, H.Krieg, Aconsistentandaccurateab initio parametrization of density functional dispersion correction (DFT- D) for the 94 elements H-Pu, Journal of Chemical Physics 132 (15) (2010)

  33. [33]

    L. A. Chernozatonskii, K. P. Katin, A. I. Kochaev, M. M. Maslov, Moiré and non-twisted sp3-hybridized structures based on hexagonal boron nitride bilayers:Ab initioinsight into infrared and Raman spectra, bands structures and mechanical properties, Applied Surface Science 606 (2022) 154909

  34. [34]

    Campanera, G

    J. Campanera, G. Savini, I. Suarez-Martinez, M. Heggie, Density func- tionalcalculationsontheintricaciesofMoirépatternsongraphite, Phys- ical Review B 75 (2007) 235449

  35. [35]

    S. A. Goreinov, I. V. Oseledets, D. V. Savostyanov, E. E. Tyrtyshnikov, N. L. Zamarashkin, How to find a good submatrix, in Matrix Methods: Theory, Algorithms And Applications: Dedicated to the Memory of Gene Golub, World Scientific, 2010, pp. 247–256

  36. [36]

    S.Nosé, Aunifiedformulationoftheconstanttemperaturemoleculardy- namics methods, The Journal of chemical physics 81 (1) (1984) 511–519

  37. [37]

    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, et al., LAMMPS - a flexible simulation tool for particle- based materials modeling at the atomic, meso, and continuum scales, Computer Physics Communications 271 (2022) 108171. 29

  38. [38]

    A. Togo, L. Chaput, T. Tadano, I. Tanaka, Implementation strategies in phonopy and phono3py, Journal of Physics: Condensed Matter 35 (35) (2023) 353001

  39. [39]

    A. Togo, L. Chaput, I. Tanaka, Distributions of phonon lifetimes in brillouin zones, Phys. Rev. B 91 (2015) 094306

  40. [40]

    Mizokami, A

    K. Mizokami, A. Togo, I. Tanaka, Lattice thermal conductivities of two SiO2 polymorphs by first-principles calculations and the phonon Boltz- mann transport equation, Phys. Rev. B 97 (2018) 224306

  41. [41]

    Parlinski, Z

    K. Parlinski, Z. Q. Li, Y. Kawazoe, First-principles determination of the soft mode in cubic ZrO2, Physical Review Letters 78 (1997) 4063–4066

  42. [42]

    A. Togo, I. Tanaka, First principles phonon calculations in materials science, Scripta Materialia 108 (2015) 1–5

  43. [43]

    Mortazavi, Electronic, thermal and mechanical properties of carbon and boron nitride holey graphyne monolayers, Materials 16 (20) (2023) 6642

    B. Mortazavi, Electronic, thermal and mechanical properties of carbon and boron nitride holey graphyne monolayers, Materials 16 (20) (2023) 6642

  44. [44]

    Y. Li, X. Li, B. Wei, J. Liu, F. Pan, H. Wang, P. Cheng, H. Zhang, D. Xu, W. Bao, et al., Phonon coherence in bismuth-halide perovskite Cs3Bi2Br9 withultralowthermalconductivity, AdvancedFunctionalMa- terials (2024) 2411152

  45. [45]

    Dangić, G

    D. Dangić, G. Caldarelli, R. Bianco, I. Savić, I. Errea, Lattice ther- mal conductivity in the anharmonic overdamped regime, arXiv preprint arXiv:2410.13485 (2024)

  46. [46]

    Fiorentino, S

    A. Fiorentino, S. Baroni, From Green-Kubo to the full Boltzmann ki- netic approach to heat transport in crystals and glasses, Phys. Rev. B 107 (2023) 054311

  47. [47]

    Carbogno, R

    C. Carbogno, R. Ramprasad, M. Scheffler,Ab InitioGreen-Kubo Ap- proach for the Thermal Conductivity of Solids, Phys. Rev. Lett. 118 (2017) 175901

  48. [48]

    Knoop, M

    F. Knoop, M. Scheffler, C. Carbogno,Ab initioGreen-Kubo simulations of heat transport in solids: Method and implementation, Phys. Rev. B 107 (2023) 224304. 30

  49. [49]

    Castellano, J.P.A

    A. Castellano, J.P.A. Batista, O. Hellman, M.J. Verstraete, Mode- coupling formulation of heat transport in anharmonic materials, Phys. Rev. B 111 (2025) 094306

  50. [50]

    Boone, H

    P. Boone, H. Babaei, C. E. Wilmer, Heat flux for many-body interac- tions: corrections to LAMMPS, Journal of chemical theory and compu- tation 15 (10) (2019) 5579–5587

  51. [51]

    S. T. Tai, C. Wang, R. Cheng, Y. Chen, Revisiting many-body interac- tion heat current and thermal conductivity calculations using the mo- ment tensor potential/LAMMPS interface, Journal of Chemical Theory and Computation 21 (7) (2025) 3649-3657

  52. [52]

    Zhang, G

    H. Zhang, G. Lee, K. Cho, Thermal transport in graphene and effects of vacancy defects, Phys. Rev. B 84 (2011) 115460

  53. [53]

    F. Saiz, C. Da Silva, C.H. Amon, Prediction of Thermal Conductivity of Two-Dimensional Superlattices of Graphene and Boron Nitride by Equi- librium Molecular Dynamics, Heat Transfer and Thermal Engineering, American Society of Mechanical Engineers 8B (2015) V08BT10A026

  54. [54]

    Pereira, D

    L.F.C. Pereira, D. Donadio, Divergence of the thermal conductivity in uniaxially strained graphene, Phys. Rev. B 87 (2013) 125424

  55. [55]

    Eriksson, E

    F. Eriksson, E. Fransson, P. Erhart, The hiphive package for the ex- traction of high-order force constants by machine learning, Advanced Theory and Simulations 2 (5) (2019) 1800184

  56. [56]

    V. Blum, R. Gehrke, F. Hanke, P. Havu, V. Havu, X. Ren, K. Reuter, M. Scheffler,Ab initiomolecular simulations with numeric atom-centered orbitals, Computer Physics Communications 180 (11) (2009) 2175–2196

  57. [57]

    Miglio, V

    A. Miglio, V. Brousseau-Couture, E. Godbout, G. Antonius, Y.-H. Chan, S. G. Louie, M. Côté, M. Giantomassi, X. Gonze, Predominance of non-adiabatic effects in zero-point renormalization of the electronic band gap, npj Computational Materials 6 (1) (2020) 167

  58. [58]

    Lindsay, Isotope scattering and phonon thermal conductivity in light atom compounds: LiH and LiF, Physical Review B 94 (17) (2016) 174304

    L. Lindsay, Isotope scattering and phonon thermal conductivity in light atom compounds: LiH and LiF, Physical Review B 94 (17) (2016) 174304. 31

  59. [59]

    Y. P. Varshni, Temperature dependence of the energy gap in semicon- ductors, Physica 34 (1) (1967) 149–154. 32