Pith. sign in

REVIEW 3 major objections 4 minor 125 references

Heat transport in crystalline organic semiconductors: coexistence of phonon propagation and tunneling

T0 review · 3 major / 4 minor · reviewed 2026-08-11 · deepseek-v4-flash

Pith's one-line read Phonon tunneling carries heat in organic crystals.

desk verdict First WTE treatment of organic semiconductors that plausibly resolves the BTE underestimate for naphthalene and pentacene, though the full acene-series claims are shakier than the paper's tone suggests. read the letter →

arxiv 2412.05062 v1 pith:AWHEMOYZ submitted 2024-12-06 cond-mat.mtrl-sci

classification cond-mat.mtrl-sci
keywords thermalconductivityphonontunnelingWignertransportequationorganicsemiconductorsacenesmachine-learnedinteratomicpotentialsanharmoniclatticedynamicspopulationandcoherencecontributions
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper argues that the ultra-low thermal conductivity of crystalline acenes cannot be captured by the standard picture of heat carried by particle-like phonons. It shows that adding a second, wave-like channel, phonon tunneling between overlapping vibration bands as formulated in the Wigner transport equation, makes computed conductivities match measured values for naphthalene and pentacene. This matters because the tunneling contribution grows with temperature and with molecular length, explaining the weak temperature dependence of naphthalene, the almost temperature-independent conductivity of pentacene, and a predicted minimum in pentacene's conductivity along the molecular axis near 300 K.

What carries the argument

The machinery is the Wigner transport equation (WTE), which splits the lattice thermal conductivity into the population (propagation) contribution, identical to the Boltzmann term, and the coherence (tunneling) contribution, driven by inter-mode coupling between phonon bands whose linewidths overlap. The coherence term is expressed through mode frequencies, heat capacities, generalized velocity operators, and Lorentzian functions of frequency differences and linewidths. To supply the required second- and third-order force constants for crystals with up to 216 modes per cell, the authors train system-specific Moment Tensor Potentials on dispersion-corrected DFT data and evaluate the WTE with these potentials, using the tetrahedron method for Brillouin-zone integration and the relaxation-time approximation after validating it against the linearized Boltzmann solution.

What would settle it

Measure the thermal conductivity of a pentacene single crystal along the molecular long axis between 100 K and 500 K. The paper predicts a minimum near 300 K followed by a rise at higher temperatures, while a particle-only theory predicts monotonic decrease; a monotonic decrease would falsify the large tunneling contribution. A complementary test is to recompute the third-order force constants of anthracene with a more accurate method and check whether the resulting linewidths match measured Raman linewidths; if they broaden substantially, the reported conductivities would change.

Watch

Extended reading notes

Core claim

The central claim is that in all four acenes studied, heat is carried by two coexisting mechanisms: particle-like phonon propagation and wave-like phonon tunneling. The Peierls-Boltzmann equation describes only the first and systematically underestimates the measured thermal conductivity; the Wigner transport equation, whose coherence term couples modes whose broadened bands overlap, accounts for both. With the propagation and tunneling channels combined, the calculated isotropic conductivities of naphthalene and pentacene agree quantitatively with experiment across wide temperature ranges, and the tunneling channel is responsible for the observed weak or vanishing temperature dependence, the growth of conductivity with molecular length along the backbone direction, and the appearance of a conductivity minimum near 300 K for pentacene.

Load-bearing premise

The argument depends on the machine-learned potentials reproducing the true third-order force constants closely enough that the calculated phonon linewidths, and therefore the size of the tunneling term, are trustworthy; for anthracene and tetracene the force errors lie above a threshold previously associated with 2% conductivity accuracy, and the computed Raman linewidths are systematically narrower than measured ones.

Editorial extensions

If this is right

  • Predictions for naphthalene and pentacene match measured conductivities only when the tunneling channel is included, so organic-semiconductor heat transport calculations that use the Boltzmann equation alone will systematically underestimate thermal conductivity.
  • The same compensation mechanism explains why naphthalene's conductivity falls only weakly with temperature and pentacene's is nearly temperature-invariant near 300 K.
  • Along the molecular backbone direction, the tunneling channel grows with molecular length, so tetracene and pentacene should show a minimum and then a rise in conductivity with temperature above roughly 350 K and 300 K, respectively.
  • The predicted thermal-conductivity anisotropy is opposite to the electrical-conductivity anisotropy: heat prefers the molecular backbone direction, while charge transport prefers the herringbone plane.
  • The machine-learned potential workflow makes anharmonic lattice dynamics feasible for molecular crystals with large unit cells, extending first-principles-quality heat transport calculations well beyond the acenes.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The compensation mechanism is likely generic to molecular crystals with many closely spaced optical bands; one testable extension is to apply the same WTE workflow to rubrene or functionalized acenes and look for the same flattening of conductivity with temperature.
  • Because the simulated Raman linewidths are systematically narrower than measured ones, the third-order force constants likely overestimate phonon lifetimes; if so, the absolute tunneling values may shift with temperature, although the paper's linewidth-rescaling test suggests the qualitative trends would survive.
  • A direct experimental test would be direction-resolved thermal conductivity of a pentacene single crystal: the predicted minimum near 300 K along the molecular backbone direction and the subsequent increase are specific, falsifiable signatures of the tunneling channel.
  • If correct, this picture implies that thermal management in organic devices can be engineered by choosing molecular length to control optical-band overlap, rather than only by reducing disorder or grain-boundary scattering.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper computes the lattice thermal conductivities of naphthalene, anthracene, tetracene, and pentacene using the Wigner transport equation, with harmonic and third-order force constants obtained from system-specific Moment Tensor Potentials trained on dispersion-corrected DFT data. The central claim is that the Peierls-Boltzmann picture, retaining only particle-like phonon propagation, systematically underestimates measured conductivities, whereas the full WTE including wave-like phonon tunneling quantitatively reproduces experiments for naphthalene and pentacene over the measured temperature ranges. The authors further analyze mode-resolved and directional contributions, finding that acoustic modes conduct mainly by propagation, dense optical bands contribute significantly by tunneling, and the tunneling conductivity along the molecular backbone increases with molecular length, leading to predicted nonmonotonic temperature dependence in tetracene and pentacene.

Significance. If the quantitative agreement holds, the paper is significant: it extends the Wigner transport formulation to crystalline organic semiconductors, gives a reciprocal-space decomposition of heat transport into propagation and tunneling channels for the acene family, and provides falsifiable predictions such as the [001] conductivity minimum in pentacene and the crossover to tunneling-dominated transport along the molecular axis. The manuscript has notable strengths: systematic convergence tests for q-meshes, supercells, and displacement amplitudes; validation of the relaxation-time approximation against a direct LBTE solution for naphthalene; inclusion of measured grain-size boundary scattering for pentacene; and a clear data-availability statement. The principal caveats concern the lack of direct benchmarking of MTP third-order force constants against DFT and the use of a single stochastic MTP per acene for all final WTE numbers.

major comments (3)
  1. [Validating the parametrized MTPs (Table 1; Supplementary Section 2.2.1)] The WTE conductivities used for the central comparison are computed with a single retained MTP per acene, and the anharmonic linewidths entering both κP and κC are never benchmarked directly against DFT third-order force constants. The retained pentacene MTP has a force RMSD of 7.69 meV/Å (Table 1), above the 5 meV/Å threshold that the authors themselves cite from Póta et al. as sufficient for about 2% thermal-conductivity accuracy. Because the MTP training is stochastic and only the best of three potentials is retained, I ask the authors to report κtot, κP, and κC for all three independently trained MTPs, or at least for the non-retained naphthalene and pentacene MTPs, and to quantify how much of the experimental agreement is a consequence of the particular potential draw.
  2. [Supplementary Section 3 and Supplementary Fig. 18] The comparison with measured Raman linewidths shows systematically narrower simulated linewidths for anthracene and pentacene, which is consistent with overestimated phonon lifetimes. The ±25% global rescaling in Supplementary Fig. 24 does not bound the resulting error in the central comparison, because it multiplies all linewidths by a single factor and cannot capture mode-specific or MTP-specific inaccuracies; indeed, the Raman comparison shows that some modes are reproduced and others are not. Please provide a mode-resolved sensitivity estimate, or justify why the systematic underestimate of linewidths cannot change the conclusion that tunneling is required to rationalize the experiments.
  3. [Fig. 6b and Supplementary Section 2.6] The anthracene [010] propagation and total conductivities are omitted because the tetrahedron and Gaussian-smearing Brillouin-zone integrations do not converge with respect to the broadening parameter. Since the claim that the in-plane [100]/[010] conductivities are essentially material-independent is one of the main structure-property conclusions, the missing [010] data for anthracene leaves that claim incomplete. Please either provide a converged [010] result using an independent method, such as the adaptive smearing already used in the ShengBTE comparison or the collisional-broadening treatment cited in the text, or explicitly restrict the material-independence claim to the directions that are actually converged.
minor comments (4)
  1. [Results, 'Heat transport in polycrystalline acenes'] There are several typos: 'napthalene' should be 'naphthalene', 'comparion' should be 'comparison', and 'provided in in Fig. 3' contains a duplicated 'in'.
  2. [Eq. (1)] The quantities V and N_c in Eq. (1) are used without definition; please define them explicitly in the text or in the equation caption.
  3. [Fig. 6 caption and main text] The caption states that κtot and κP of anthracene are omitted in panel b, while the main text specifically refers to κP[010] and κtot[010]; please harmonize the wording and make the omission visible in the legend itself.
  4. [References and Supplementary text] The spelling of the naphthalene experimental reference is inconsistent: 'Ueberreiter' in the main text and 'Überreiter' in the Supplementary Information; please standardize the spelling.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the WTE thermal conductivities are first-principles predictions with no parameter fitted to the target experiments.

full rationale

The derivation chain is: DFT training data -> MTPs validated against DFT forces/frequencies -> second- and third-order force constants -> WTE/BTE conductivities -> comparison with literature experiments. The experimental thermal conductivities of naphthalene and pentacene enter only as comparison benchmarks, not as training or fitting data. The only externally set parameter in the comparison, the pentacene grain size of 256 nm, is taken from the morphological characterization in the experimental paper (ref. 68) and is not adjusted to reproduce the measured conductivity. The WTE formalism is an externally established theory (refs. 27,28), and although one co-author is among the cited developers, the paper does not invoke any unpublished uniqueness theorem or ansatz from the authors' prior work; it uses the published equations. Self-citations to prior MTP/phonon work (refs. 18,47,55) supply methodology and harmonic references that are themselves benchmarked against DFT, so they are not load-bearing circular premises. Accuracy limitations (force RMSDs above 5 meV/A for some acenes, Raman linewidths narrower than experiment) are correctness risks, not circularity: they do not make the prediction equivalent to its input by construction.

Assumptions & free parameters 0 free parameters · 5 assumptions · 0 invented entities

The central claim rests on the accuracy of MTP-derived third-order force constants and the sufficiency of the third-order WTE framework. No free parameters are fitted to thermal-conductivity data: the pentacene boundary-scattering mean free path (256 nm) is taken directly from the experimental morphology characterization of ref. 68, and numerical settings (displacement amplitudes, q-meshes, supercells) are chosen via convergence tests. No new physical entities are introduced.

assumptions (5)
  • domain assumption MTPs trained on PBE+D3BJ DFT data accurately reproduce the anharmonic force constants (third order) required for the WTE.
    The MTPs are validated against DFT phonon band structures and force sets, but third-order force constants are not directly validated against DFT because of cost; force RMSDs for anthracene and tetracene exceed 8 meV/A (Table 1).
  • domain assumption The Wigner transport equation truncated at third-order anharmonicity is sufficient to describe heat transport in these crystals.
    The paper relies on WTE as implemented in phono3py with only third-order force constants; the impact is tested only via linewidth rescaling (Supp. Fig. 24), and the anthracene [010] case suggests the completed-collision limit may be insufficient.
  • domain assumption Phonon lifetimes are not overdamped (Ioffe-Regel limit) so the quasiparticle picture applies.
    Checked in Supp. Section 2.5; all relevant modes are above 1/omega.
  • domain assumption DFT-relaxed unit cells remain valid across the whole temperature range (thermal expansion neglected).
    Tested for naphthalene: 8% difference at 300 K and trends unchanged (Supp. Section 4.2).
  • domain assumption The relaxation time approximation is a good approximation for the WTE and BTE solutions.
    Validated against direct LBTE for naphthalene, differences negligible (Supp. Fig. 7).

how reviews work

0 comments
Cite this review

Pith. "Pith review of Heat transport in crystalline organic semiconductors: coexistence of phonon propagation and tunneling." pith.science (2026). https://pith.science/paper/AWHEMOYZ

@misc{pith2026241205062,
  author       = {Pith},
  title        = {Pith review of: Heat transport in crystalline organic semiconductors: coexistence of phonon propagation and tunneling},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/AWHEMOYZ}},
  note         = {Machine review of arXiv:2412.05062}
}
read the original abstract

Understanding heat transport in organic semiconductors is of fundamental and practical relevance. Therefore, we study the lattice thermal conductivities of a series of (oligo)acenes, where an increasing number of rings per molecule leads to a systematic increase of the crystals' complexity. Temperature-dependent thermal conductivity experiments in these systems disagree with predictions based on the traditional Peierls-Boltzmann framework, which describes heat transport in terms of particle-like phonon propagation. We demonstrate that accounting for additional phonon-tunneling conduction mechanisms through the Wigner Transport Equation resolves this disagreement and quantitatively rationalizes experiments. The pronounced increase of tunneling transport with temperature explains several unusual experimental observations, such as a weak temperature dependence in naphthalene's conductivity and an essentially temperature-invariant conductivity in pentacene. While the anisotropic conductivities within the acene planes are essentially material-independent, the tunneling contributions (and hence the total conductivities) significantly increase with molecular length in the molecular backbone direction, which for pentacene results in a surprising minimum of the thermal conductivity at 300K.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

125 extracted references · 65 canonical work pages

  1. [10]

    & Scheffler, M

    Carbogno, C., Ramprasad, R. & Scheffler, M. Ab Initio Green-Kubo Approach for the Thermal Conductivity of Solids. Phys. Rev. Lett. 118, 175901 (2017)

  2. [1]

    Zhang, X., Dong, H. & Hu, W. Organic Semiconductor Single Crystals for Electronics and Photonics. Advanced Materials 30, 1801048 (2018)

  3. [2]

    Newman, C. R. et al. Introduction to Organic Thin Film Transistors and Design of n-Channel Organic Semiconductors. Chem. Mater. 16, 4436–4451 (2004)

  4. [3]

    Pfeiffer, M. et al. Doped organic semiconductors: Physics and application in light emitting diodes. Organic Electronics 4, 89–103 (2003)

  5. [4]

    P., Genoe, J., Heremans, P

    Rand, B. P., Genoe, J., Heremans, P. & Poortmans, J. Solar cells utilizing small molecular weight organic semiconductors. Progress in Photovoltaics: Research and Applications 15, 659–676 (2007)

  6. [5]

    Hong, G. et al. A Brief History of OLEDs—Emitter Development and Industry Milestones. Advanced Materials 33, 2005630 (2021)

  7. [6]

    Wang, X. et al. Thermal transport in organic semiconductors. Journal of Applied Physics 130, 170902 (2021)

  8. [7]

    A., Malorny, M., Birner, G., Mingo, N

    Broido, D. A., Malorny, M., Birner, G., Mingo, N. & Stewart, D. A. Intrinsic lattice thermal conductivity of semiconductors from first principles. Applied Physics Letters 91, 231922 (2007)

Show all 125 references
  1. [8]

    & Colombo, L

    Melis, C., Dettori, R., Vandermeulen, S. & Colombo, L. Calculating thermal conductivity in a transient System DIM Natom Ndisp q-mesh 2A (2,3,2) 432 140022 9×13×9 3A (2,3,2) 578 248904 10×10×8 4A (2,2,2) 480 518580 10×10×8 5A (2,2,2) 576 746712 9×9×7 17 conduction regime: theor...

  2. [9]

    S., Talaat, K

    El-Genk, M. S., Talaat, K. & Cowen, B. J. Thermal conductivity of silicon using reverse non-equilibrium molecular dynamics. Journal of Applied Physics 123, 205104 (2018)

  3. [11]

    Selezneva, E. et al. Strong Suppression of Thermal Conductivity in the Presence of Long Terminal Alkyl Chains in Low-Disorder Molecular Semiconductors. Advanced Materials 33, 2008708 (2021)

  4. [12]

    Gueye, M

    N. Gueye, M. et al. Thermal conductivity of benzothieno-benzothiophene derivatives at the nanoscale. Nanoscale 13, 3800–3807 (2021)

  5. [13]

    & Cornil, J

    Vercouter, A., Lemaur, V., Melis, C. & Cornil, J. Computing the Lattice Thermal Conductivity of Small- Molecule Organic Semiconductors: A Systematic Comparison of Molecular Dynamics Based Methods. Advanced Theory and Simulations n/a, 2200892

  6. [14]

    & Shuai, Z

    Wang, D., Tang, L., Long, M. & Shuai, Z. Anisotropic Thermal Transport in Organic Molecular Crystals from Nonequilibrium Molecular Dynamics Simulations. J. Phys. Chem. C 115, 5940–5946 (2011)

  7. [15]

    & Wang, X

    Yang, C., Wang, W., Peng, B., Ji, W. & Wang, X. Insight into the effect of side chains on thermal transport of organic semiconductors. Nanoscale (2023) doi:10.1039/D3NR04275H

  8. [16]

    & Hopkins, P

    Giri, A. & Hopkins, P. E. Spectral Contributions to the Thermal Conductivity of C60 and the Fullerene Derivative PCBM. J. Phys. Chem. Lett. 8, 2153–2157 (2017)

  9. [17]

    Zur kinetischen Theorie der Wärmeleitung in Kristallen

    Peierls, R. Zur kinetischen Theorie der Wärmeleitung in Kristallen. Annalen der Physik 395, 1055–1101 (1929)

  10. [18]

    & Zojer, E

    Kamencek, T. & Zojer, E. Discovering structure– property relationships for the phonon band structures of hydrocarbon-based organic semiconductor crystals: the instructive case of acenes. Journal of Materials Chemistry C 10, 2532–2543 (2022)

  11. [19]

    & Zojer, E

    Legenstein, L., Reicht, L., Kamencek, T. & Zojer, E. Anisotropic Phonon Bands in H-Bonded Molecular Crystals: The Instructive Case of α-Quinacridone. ACS Mater. Au 3, 371–385 (2023)

  12. [20]

    & Galli, G

    Puligheddu, M., Xia, Y., Chan, M. & Galli, G. Computational prediction of lattice thermal conductivity: A comparison of molecular dynamics and Boltzmann transport approaches. Phys. Rev. Mater. 3, 085401 (2019)

  13. [21]

    & Stokes, H

    Esfarjani, K., Chen, G. & Stokes, H. T. Heat transport in silicon from first-principles calculations. Phys. Rev. B 84, 085204 (2011)

  14. [22]

    & Mauri, F

    Fugallo, G., Lazzeri, M., Paulatto, L. & Mauri, F. Ab initio variational approach for evaluating lattice thermal conductivity. Phys. Rev. B 88, 045430 (2013)

  15. [23]

    McGaughey, A. J. H., Jain, A., Kim, H.-Y. & Fu, B. ( 傅博). Phonon properties and thermal conductivity from first principles, lattice dynamics, and the Boltzmann transport equation. Journal of Applied Physics 125, 011101 (2019)

  16. [24]

    & Ruan, X

    Luo, Y., Yang, X., Feng, T., Wang, J. & Ruan, X. Vibrational hierarchy leads to dual-phonon transport in low thermal conductivity crystals. Nat Commun 11, 2554 (2020)

  17. [25]

    Lee, W. et al. Ultralow thermal conductivity in all- inorganic halide perovskites. Proceedings of the National Academy of Sciences 114, 8693–8697 (2017)

  18. [26]

    & Marzari, N

    Di Lucente, E., Simoncelli, M. & Marzari, N. Crossover from Boltzmann to Wigner thermal transport in thermoelectric skutterudites. Phys. Rev. Res. 5, 033125 (2023)

  19. [27]

    & Mauri, F

    Simoncelli, M., Marzari, N. & Mauri, F. Wigner Formulation of Thermal Transport in Solids. Phys. Rev. X 12, 041011 (2022)

  20. [28]

    & Mauri, F

    Simoncelli, M., Marzari, N. & Mauri, F. Unified theory of thermal transport in crystals and glasses. Nat. Phys. 15, 809–813 (2019)

  21. [29]

    Li, Y. et al. Phonon Coherence in Bismuth-Halide Perovskite Cs3Bi2Br9 With Ultralow Thermal Conductivity. Advanced Functional Materials n/a, 2411152 (2024)

  22. [30]

    Zheng, J. et al. Unravelling ultralow thermal conductivity in perovskite Cs2AgBiBr6: dominant wave-like phonon tunnelling and strong anharmonicity. npj Comput Mater 10, 1–13 (2024)

  23. [31]

    Tong, Z. et al. Predicting the Lattice Thermal Conductivity in Nitride Perovskite LaWN3 from ab initio Lattice Dynamics. Advanced Science 10, 2205934 (2023)

  24. [32]

    & Ong, W.-L

    Yang, J., Jain, A. & Ong, W.-L. Inter-channel conversion between population-/coherence-channel dictates thermal transport in MAPbI3 crystals. Materials Today Physics 28, 100892 (2022)

  25. [33]

    Simoncelli, M. et al. Temperature-invariant heat conductivity from compensating crystalline and glassy transport: from the Steinbach meteorite to furnace bricks. Preprint at https://doi.org/10.21203/rs.3.rs- 4456620/v1 (2024)

  26. [34]

    & Simoncelli, M

    Pazhedath, A., Bastonero, L., Marzari, N. & Simoncelli, M. First-principles characterization of thermal conductivity in LaPO4-based alloys. Phys. Rev. Appl. 22, 024064 (2024)

  27. [35]

    & Lindsay, L

    Thébaud, S., Berlijn, T. & Lindsay, L. Perturbation theory and thermal transport in mass-disordered alloys: Insights from Green’s function methods. Phys. Rev. B 105, 134202 (2022). 18 36.Gross, R. & Marx, A. Festkörperphysik: Kapitel 6.4. (De Gruyter, 2018). doi:10.1515/9783110559187

  28. [38]

    C., Albinati, A., Mason, S

    Capelli, S. C., Albinati, A., Mason, S. A. & Willis, B. T. M. Molecular Motion in Crystalline Naphthalene: Analysis of Multi-Temperature X-Ray and Neutron Diffraction Data. J. Phys. Chem. A 110, 11695–11703 (2006)

  29. [42]

    I., Nelson, S

    Lin, Y.-Y., Gundlach, D. I., Nelson, S. F. & Jackson, T. N. Pentacene-based organic thin-film transistors. IEEE Transactions on Electron Devices 44, 1325–1331 (1997)

  30. [43]

    E., Brooks, J

    Anthony, J. E., Brooks, J. S., Eaton, D. L. & Parkin, S. R. Functionalized Pentacene: Improved Electronic Properties from Control of Solid-State Order. J. Am. Chem. Soc. 123, 9482–9483 (2001)

  31. [44]

    He, Z., Zhang, Z., Bi, S., Chen, J. & Li, D. Conjugated Polymer Controlled Morphology and Charge Transport of Small-Molecule Organic Semiconductors. Sci Rep 10, 4344 (2020)

  32. [45]

    Visualization and analysis of atomistic simulation data with OVITO–the Open Visualization Tool

    Stukowski, A. Visualization and analysis of atomistic simulation data with OVITO–the Open Visualization Tool. Modelling Simul. Mater. Sci. Eng. 18, 015012 (2009)

  33. [46]

    Kamencek, T. et al. Evaluating Computational Shortcuts in Supercell-Based Phonon Calculations of Molecular Crystals: The Instructive Case of Naphthalene. J. Chem. Theory Comput. 16, 2716– 2735 (2020)

  34. [47]

    & Zojer, E

    Wieser, S. & Zojer, E. Machine learned force-fields for an Ab-initio quality description of metal-organic frameworks. npj Comput Mater 10, 1–18 (2024)

  35. [48]

    Shapeev, A. V. Moment Tensor Potentials: A Class of Systematically Improvable Interatomic Potentials. Multiscale Model. Simul. 14, 1153–1173 (2016)

  36. [49]

    & Shapeev, A

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

  37. [50]

    Liu, H., Qian, X., Bao, H., Zhao, C. Y. & Gu, X. High- temperature phonon transport properties of SnSe from machine-learning interatomic potential. J. Phys.: Condens. Matter 33, 405401 (2021)

  38. [51]

    Mortazavi, B. et al. Accelerating first-principles estimation of thermal conductivity by machine- learning interatomic potentials: A MTP/ShengBTE solution. Computer Physics Communications 258, 107583 (2021)

  39. [52]

    Ouyang, Y. et al. Accurate description of high- order phonon anharmonicity and lattice thermal conductivity from molecular dynamics simulations with machine learning potential. Phys. Rev. B 105, 115202 (2022)

  40. [53]

    Cui, C. et al. Machine learning interatomic potentials as efficient tools for obtaining reasonable phonon dispersions and accurate thermal conductivity: A case study of typical two-dimensional materials. Applied Physics Letters 123, 152201 (2023)

  41. [54]

    & Shapeev, A

    Rybin, N. & Shapeev, A. A moment tensor potential for lattice thermal conductivity calculations of α and β phases of Ga2O3. Journal of Applied Physics 135, 205108 (2024)

  42. [55]

    & Zojer, E

    Reicht, L., Legenstein, L., Wieser, S. & Zojer, E. Designing Accurate Moment Tensor Potentials for Phonon-Related Properties of Crystalline Polymers. Molecules 29, 3724 (2024)

  43. [56]

    Ziman, J. M. Electrons and Phonons: The Theory of Transport Phenomena in Solids. (Oxford University Press, Oxford, 2001). doi:10.1093/acprof:oso/9780198507796.001.0001

  44. [59]

    D., Gardner, J

    Morrow, J. D., Gardner, J. L. A. & Deringer, V. L. How to validate machine-learned interatomic potentials. The Journal of Chemical Physics 158, 121501 (2023)

  45. [60]

    & Simoncelli, M

    Póta, B., Ahlawat, P., Csányi, G. & Simoncelli, M. Thermal Conductivity Predictions with Foundation Atomistic Models. Preprint at http://arxiv.org/abs/2408.00755 (2024)

  46. [61]

    & Chen, Y

    Wang, Q., Zeng, Z. & Chen, Y. Revisiting phonon transport in perovskite SrTiO3: Anharmonic phonon renormalization and four-phonon scattering. Phys. Rev. B 104, 235205 (2021)

  47. [62]

    Liu, Z., Yang, X., Zhang, B. & Li, W. High Thermal Conductivity of Wurtzite Boron Arsenide Predicted by Including Four-Phonon Scattering with Machine Learning Potential. ACS Appl. Mater. Interfaces 13, 53409–53415 (2021)

  48. [63]

    & Chen, Y

    Wang, C., Wang, Q., Zhang, Q., Chen, C. & Chen, Y. Intrinsic Zn Vacancies-Induced Wavelike Tunneling of 19 Phonons and Ultralow Lattice Thermal Conductivity in Zintl Phase Sr2ZnSb2. Chem. Mater. 34, 7837–7844 (2022)

  49. [64]

    & Feng, T

    Zhou, H., Tiwari, J. & Feng, T. Understanding the flat thermal conductivity of La2Zr2O7 at ultrahigh temperatures. Phys. Rev. Materials 8, 043804 (2024)

  50. [65]

    Shen, X. et al. Amorphous-Like Ultralow Thermal Transport in Crystalline Argyrodite Cu7PS6. Advanced Science 11, 2400258 (2024)

  51. [66]

    & Errea, I

    Dangić, Đ., Caldarelli, G., Bianco, R., Savić, I. & Errea, I. Lattice thermal conductivity in the anharmonic overdamped regime. Preprint at http://arxiv.org/abs/2410.13485 (2024)

  52. [67]

    & Orthmann, H.-J

    Ueberreiter, K. & Orthmann, H.-J. Spezifische Wärme, spezifisches Volumen, Temperatur-und Wärmeleitfähigkeit einiger disubstituierter Benzole und polycyclischer Systeme. Zeitschrift für Naturforschung A 5, 101–108 (1950)

  53. [68]

    Epstein, J., Ong, W.-L., Bettinger, C. J. & Malen, J. A. Temperature Dependent Thermal Conductivity and Thermal Interface Resistance of Pentacene Thin Films with Varying Morphology. ACS Appl. Mater. Interfaces 8, 19168–19174 (2016)

  54. [69]

    Bounds on the conductivity of statistically isotropic polycrystals

    Schulgasser, K. Bounds on the conductivity of statistically isotropic polycrystals. J. Phys. C: Solid State Phys. 10, 407–417 (1977)

  55. [70]

    G., Andersson, P

    Ross, R. G., Andersson, P. & Bäckström, G. Thermal conductivity and heat capacity of benzene, naphthalene and anthracene under pressure. Molecular Physics 38, 527–533 (1979)

  56. [71]

    Y., Roth, S

    Lee, J. Y., Roth, S. & Park, Y. W. Anisotropic field effect mobility in single crystal pentacene. Applied Physics Letters 88, 252106 (2006)

  57. [72]

    P., Shim, J

    Nguyen, T. P., Shim, J. H. & Lee, J. Y. Density Functional Theory Studies of Hole Mobility in Picene and Pentacene Crystals. J. Phys. Chem. C 119, 11301– 11310 (2015)

  58. [73]

    & Benfatto, L

    Caldarelli, G., Simoncelli, M., Marzari, N., Mauri, F. & Benfatto, L. Many-body Green’s function approach to lattice thermal transport. Phys. Rev. B 106, 024312 (2022)

  59. [74]

    & Berkelbach, T

    Jasrasaria, D. & Berkelbach, T. C. Strong anharmonicity dictates ultralow thermal conductivities of type-I clathrates. Preprint at http://arxiv.org/abs/2409.08242 (2024)

  60. [75]

    Evaluation of the Carrier–Quasiparticle Scattering Superoperator

    Rossi, F. Evaluation of the Carrier–Quasiparticle Scattering Superoperator. in Theory of Semiconductor Quantum Devices: Microscopic Modeling and Simulation Strategies (ed. Rossi, F.) 353–355 (Springer, Berlin, Heidelberg, 2011). doi:10.1007/978-3-642- 10556-2_14

  61. [76]

    Zeller, R. C. & Pohl, R. O. Thermal Conductivity and Specific Heat of Noncrystalline Solids. Phys. Rev. B 4, 2029–2041 (1971)

  62. [77]

    S., Gubaev, K., Podryabinkin, E

    Novikov, I. S., Gubaev, K., Podryabinkin, E. V. & Shapeev, A. V. The MLIP package: moment tensor potentials with MPI and active learning. Mach. Learn.: Sci. Technol. 2, 025002 (2020)

  63. [78]

    Lee, N.-E., Zhou, J.-J., Agapito, L. A. & Bernardi, M. Charge transport in organic molecular semiconductors from first principles: The bandlike hole mobility in a naphthalene crystal. Phys. Rev. B 97, 115203 (2018)

  64. [79]

    & Marktanner, J

    Karl, N. & Marktanner, J. Electron and Hole Mobilities in High Purity Anthracene Single Crystals. Molecular Crystals and Liquid Crystals Science and Technology. Section A. Molecular Crystals and Liquid Crystals 355, 149–173 (2001)

  65. [80]

    & Bao, H

    Zhang, H., Gu, X., Fan, Z. & Bao, H. Vibrational anharmonicity results in decreased thermal conductivity of amorphous HfO2 at high temperature. Phys. Rev. B 108, 045422 (2023)

  66. [81]

    Lee, D. W. & Kingery, W. D. Radiation Energy Transfer and Thermal Conductivity of Ceramic Oxides. Journal of the American Ceramic Society 43, 594–607 (1960)

  67. [84]

    & Furthmüller, J

    Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Computational Materials Science 6, 15–50 (1996)

  68. [85]

    & Kresse, G

    Jinnouchi, R., Karsai, F. & Kresse, G. On-the-fly machine learning force field generation: Application to melting points. Phys. Rev. B 100, 014105 (2019)

  69. [86]

    & Asahi, R

    Jinnouchi, R., Miwa, K., Karsai, F., Kresse, G. & Asahi, R. On-the-Fly Active Learning of Interatomic Potentials for Large-Scale Atomistic Simulations. J. Phys. Chem. Lett. 11, 6946–6955 (2020)

  70. [87]

    GitLab https://gitlab.com/ashapeev/interface- lammps-mlip-2 (2024)

    Alexander Shapeev / LAMMPS-MLIP interface · GitLab. GitLab https://gitlab.com/ashapeev/interface- lammps-mlip-2 (2024)

  71. [88]

    Thompson, A. P. et al. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Computer Physics Communications 271, 108171 (2022)

  72. [89]

    Allen, M. P. & Tildesley, D. J. Computer Simulation of Liquids, Chapter 3: Molecular Dynamics. (Oxford University Press, Oxford, United Kingdom, 2017)

  73. [90]

    & Kresse, G

    Verdi, C., Karsai, F., Liu, P., Jinnouchi, R. & Kresse, G. Thermal transport and phase transitions of zirconia by on-the-fly machine-learned interatomic potentials. npj Comput Mater 7, 1–9 (2021)

  74. [91]

    VASP Wiki

    Best practices for machine-learned force fields. VASP Wiki. 20

  75. [92]

    P., Burke, K

    Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 77, 3865–3868 (1996)

  76. [93]

    & Krieg, H

    Grimme, S., Antony, J., Ehrlich, S. & Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. The Journal of Chemical Physics 132, 154104 (2010)

  77. [94]

    & Goerigk, L

    Grimme, S., Ehrlich, S. & Goerigk, L. Effect of the damping function in dispersion corrected density functional theory. Journal of Computational Chemistry 32, 1456–1465 (2011)

  78. [95]

    Blöchl, P. E. Projector augmented-wave method. Phys. Rev. B 50, 17953–17979 (1994)

  79. [96]

    GitLab https://gitlab.com/ashapeev/mlip-2 (2023)

    Alexander Shapeev / MLIP version 2 · GitLab. GitLab https://gitlab.com/ashapeev/mlip-2 (2023)

  80. [97]

    Parlinski, K., Li, Z. Q. & Kawazoe, Y. First-Principles Determination of the Soft Mode in Cubic ZrO2. Phys. Rev. Lett. 78, 4063–4066 (1997)

  81. [98]

    & Tanaka, I

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

  82. [99]

    First-principles Phonon Calculations with Phonopy and Phono3py

    Togo, A. First-principles Phonon Calculations with Phonopy and Phono3py. J. Phys. Soc. Jpn. 92, 012001 (2023)

  83. [101]

    H., Vosko, S

    MacDonald, A. H., Vosko, S. H. & Coleridge, P. T. Extensions of the tetrahedron method for evaluating spectral properties of solids. J. Phys. C: Solid State Phys. 12, 2991–3002 (1979)

  84. [102]

    E., Jepsen, O

    Blöchl, P. E., Jepsen, O. & Andersen, O. K. Improved tetrahedron method for Brillouin-zone integrations. Phys. Rev. B 49, 16223–16233 (1994)

  85. [104]

    & Marzari, N

    Simoncelli, M., Mauri, F. & Marzari, N. Thermal conductivity of glasses: first-principles theory and applications. npj Comput Mater 9, 1–22 (2023)

  86. [105]

    Иоффе, А. Ф. & Regel, A. R. Non-Crystalline, Amorphous, and Liquid Electronic Semiconductors. 237–291 (1960). Acknowledgments We acknowledge Tomas Kamencek and his work in laying the foundation for this publication. We also would like to thank Lukas Hörmann for stimulating dis...

  87. [106]

    accurate

    (e-h) for naphthalene, anthracene, tetracene, and pentacene (left to right). The initial crystal structures of the stable, monoclinic and triclinic acene polymorphs studied in this work were obtained from the Cambridge Crystallographic Data Centre1 and their identifiers are li...

  88. [107]

    Ag”) and close-lying intermolecular phonon bands, but for naphthalene and anthracene the additional “Bg

    drops monotonic and at 300 K it is roughly 40% smaller than for the tetrahedron method. To check whether the integration methods would also produce diverging results, the same test was performed for naphthalene. As indicated already above, the Gaussian smearing and tetrahedron...

  89. [108]

    acoustic participation ratio

    in panel c), according to eq.(3) in the main work. Supplementary Section 6: Acoustic phonons in the acene crystals In the main text, we refer to the acoustic phonons bands several times and a lthough, they are in principle a well-understood concept, we want to provide an illus...

  90. [109]

    R., Bruno, I

    Groom, C. R., Bruno, I. J., Lightfoot, M. P. & Ward, S. C. The Cambridge Structural Database. Acta Cryst B 72, 171–179 (2016)

  91. [110]

    C., Albinati , A., Mason, S

    Capelli, S. C., Albinati , A., Mason, S. A. & Willis, B. T. M. Molecular Motion in Crystalline Naphthalene: Analysis of Multi-Temperature X-Ray and Neutron Diffraction Data. J. Phys. Chem. A 110, 11695–11703 (2006)

  92. [111]

    Brock, C. P. & Dunitz, J. D. Temperature dependence of thermal motion in crystalline anthracene. Acta Cryst B 46, 795–806 (1990)

  93. [112]

    Holmes, D., Kumaraswamy, S., Matzger, A. J. & Vollhardt, K. P. C. On the Nature of Nonplanarity in the [N]Phenylenes. Chemistry – A European Journal 5, 3399–3412 (1999)

  94. [113]

    B., Robertson, J

    Campbell, R. B., Robertson, J. M. & Trotter, J. The crystal structure of hexacene, and a revision of the crystallographic data for tetracene. Acta Cryst 15, 289–290 (1962)

  95. [114]

    & Zojer, E

    Kamencek, T. & Zojer, E. Discovering structure –property relationships for the phonon band structures of hydrocarbon-based organic semiconductor crystals: the instructive case of acenes. Journal of Materials Chemistry C 10, 2532–2543 (2022)

  96. [115]

    & Harms, F

    Niggli, P., Wien, W. & Harms, F. Handbuch der Experimentalphysik. vol. Vol. 7 (Akad. Verlag-Ges., Leipzig, 1928)

  97. [116]

    & Hafner, J

    Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. Rev. B 47, 558–561 (1993)

  98. [117]

    & Furthmüller, J

    Kresse, G. & Furthmüller, J. Efficiency of ab -initio total energy calculations for metals and semiconductors using a plane-wave basis set. Computational Materials Science 6, 15–50 (1996)

  99. [118]

    & Furthmüller, J

    Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54, 11169–11186 (1996)

  100. [119]

    VASP Wiki

    Best practices for machine-learned force fields. VASP Wiki

  101. [120]

    NOMAD dataset: Phonons Crystalline Acenes

    Kamencek, T. NOMAD dataset: Phonons Crystalline Acenes. Preprint at https://doi.org/10.17172/NOMAD/2021.09.28-1 (2021)

  102. [121]

    Kamencek, T. et al. Evaluating Computational Shortcuts in Supercell -Based Phonon Calculations of Molecular Crystals: The Instructive Case of Naphthalene. J. Chem. Theory Comput. 16, 2716– 2735 (2020)

  103. [122]

    & Zojer, E

    Legenstein, L., Reicht, L., Kamencek, T. & Zojer, E. Anisotropic Phonon Bands in H -Bonded Molecular Crystals: The Instructive Case of α-Quinacridone. ACS Mater. Au 3, 371–385 (2023)

  104. [123]

    https://vsc.ac.at//systems/vsc-5/

    VSC: VSC-5. https://vsc.ac.at//systems/vsc-5/

  105. [124]

    Direct Solution to the Linearized Phonon Boltzmann Equation

    Chaput, L. Direct Solution to the Linearized Phonon Boltzmann Equation. Phys. Rev. Lett. 110, 265506 (2013)

  106. [125]

    & Tanaka, I

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

  107. [126]

    Иоффе, А. Ф. & Regel, A. R. Non -Crystalline, Amorphous, and Liquid Electronic Semiconductors. 237–291 (1960)

  108. [127]

    Katcho, N

    Li, W., Carrete, J., A. Katcho, N. & Mingo, N. ShengBTE: A solver of the Boltzmann transport equation for phonons. Computer Physics Communications 185, 1747–1758 (2014)

  109. [128]

    & Tanaka, I

    Togo, A., Chaput, L., Tadano, T. & Tanaka, I. Implementation strategies in phonopy and phono3py. J. Phys.: Condens. Matter 35, 353001 (2023)

  110. [129]

    & Erhart, P

    Eriksson, F., Fransson, E. & Erhart, P. The Hiphive Package for the Extraction of High-Order Force Constants by Machine Learning. Advanced Theory and Simulations 2, 1800184 (2019)

  111. [130]

    Asher, M. et al. Anharmonic Lattice Vibrations in Small -Molecule Organic Semiconductors. Advanced Materials 32, 1908028 (2020)

  112. [131]

    Bellows, J. C. & Prasad, P. N. Dephasing times and linewidths of optical transitions in molecular crystals. Temperature dependence of line shapes, linewidths, and frequencies of Raman active phonons in naphthalene. The Journal of Chemical Physics 70, 1864–1871 (1979). 37 / 37

  113. [132]

    Hess, L. A. & Prasad, P. N. Vibrational dephasing in organic solids: Temperature dependence of a Raman active localized internal mode of naphthalene. The Journal of Chemical Physics 72, 573– 579 (1980)

  114. [133]

    & Orthmann, H

    Ueberreiter, K. & Orthmann, H. -J. Spezifische Wärm e, spezifisches Volumen, Temperatur -und Wärmeleitfähigkeit einiger disubstituierter Benzole und polycyclischer Systeme. Zeitschrift für Naturforschung A 5, 101–108 (1950)

  115. [134]

    Data to ‘Heat transport in crystalline organic semiconductors: c oexistence of phonon propagation and tunneling’

    Legenstein, L. Data to ‘Heat transport in crystalline organic semiconductors: c oexistence of phonon propagation and tunneling’. TU Graz Repository https://doi.org/10.3217/t9czy -wjy83 (2024)

  116. [135]

    F., Iwanowski, K., Witt, W

    Harper, A. F., Iwanowski, K., Witt, W. C., Payne, M. C. & Simoncelli, M. Vibrational and thermal properties of amorphous alumina from first principles. Phys. Rev. Mater. 8, 043601 (2024)

  117. [136]

    & Zojer, E

    Kamencek, T., Bedoya-Martínez, N. & Zojer, E. Supplementary Material: Understanding phonon properties in isoreticular metal-organic frameworks from first principles. Phys. Rev. Materials 3, 116003 (2019)

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.