pith. sign in

arxiv: 2509.20951 · v2 · submitted 2025-09-25 · ❄️ cond-mat.mtrl-sci · cond-mat.soft

Computing finite--temperature elastic constants with noise cancellation

Debashish Mukherji , Marcus M\"uller , Martin H. M\"user This is my paper

Pith reviewed 2026-05-18 14:15 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci cond-mat.soft
keywords elastic constantsfinite temperaturenoise cancellationmolecular dynamics simulationstress fluctuationsthermostatting schemesamorphous materials
0
0 comments X

The pith

Elastic constants at finite temperature can be computed with greatly reduced thermal noise by canceling fluctuations between strained and reference systems using identical thermostats.

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

The authors generalize a noise-cancellation technique to elastic constants by applying a small strain to an equilibrated system and simulating both the strained and unstrained versions with matching thermostatting schemes. This causes thermal fluctuations in the stress to largely cancel, leaving a cleaner signal for the elastic response. The idea is proven in theory and tested on one-dimensional models before being applied to real materials including crystals and amorphous solids. The result is a practical way to get accurate finite-temperature elastic constants without excessive simulation lengths or anharmonic interference.

Core claim

Stress differences evaluated from simulations of a slightly strained solid and its reference system, both run with identical thermostatting, permit the determination of elastic constants with substantially reduced thermal noise in thermal ordered and disordered systems.

What carries the argument

Noise cancellation via stress subtraction between identically thermostatted strained and unstrained (or oppositely strained) configurations.

If this is right

  • Elastic constants become accessible in simulations where direct methods suffer from poor signal-to-noise.
  • The approach works for both crystalline and amorphous materials.
  • Identical thermostatting ensures cancellation without introducing extra bias.
  • Generalization from piezoelectric methods indicates it can apply to other coupling coefficients.

Where Pith is reading between the lines

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

  • Similar cancellation strategies might apply to computing other thermodynamic response functions like specific heat or thermal expansion.
  • Integration with machine learning potentials could further accelerate such calculations for large systems.
  • Testing on systems with strong anharmonicity would reveal the limits of the noise reduction.

Load-bearing premise

Applying a small strain and using identical thermostatting schemes in the strained and reference systems produces stress differences whose thermal fluctuations cancel without residual bias or anharmonic contamination.

What would settle it

A direct comparison of elastic constants computed with this noise-cancellation method against known analytical values in a one-dimensional harmonic chain or against long-time averaged results in a three-dimensional crystal would confirm or refute the noise reduction and accuracy.

Figures

Figures reproduced from arXiv: 2509.20951 by Debashish Mukherji, Marcus M\"uller, Martin H. M\"user.

Figure 1
Figure 1. Figure 1: compares the time–dependent estimate for the modulus of linear, ideally harmonic chains. To this end, a heterogeneous and a homogeneous chain are consid￾ered. The heterogeneous chain has short bonds (as = 0.9, ks = 10) and long bonds (al = 1.1, kl = 0.52631 . . .), where the exact numerical value for kl is chosen such that the static elastic modulus is E∞ = 1. For the homoge￾neous chain, as = al = 1 and ks… view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Relaxation of the elastic–tensor estimator for a lin [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Instantaneous standard deviation for data like that [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Monomer structures of poly(methyl methacrylate) [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Time evolution of the tetragonal shear–modulus esti [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Same as Figure 5, however, for silicon at a tempera [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: The temporal relaxation of the elastic modulus [PITH_FULL_IMAGE:figures/full_fig_p009_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: A summary of the components of the elastic modulus ten [PITH_FULL_IMAGE:figures/full_fig_p010_8.png] view at source ↗
read the original abstract

Elastic constants are central material properties, frequently reported in experimental and theoretical studies. While their computation is straightforward in the absence of thermal fluctuations, finite--temperature methods often suffer from poor signal--to--noise ratios or the presence of strong anharmonic effects. Here, we show how to compute elastic constants in thermal ordered and disordered systems by generalizing a noise--cancellation method originally developed for piezoelectric coupling coefficients. A slight strain is applied to an equilibrated solid. Simulations of both the strained and unstrained (or oppositely strained) reference systems are performed using identical thermostatting schemes. As demonstrated theoretically and with generic one--dimensional models, this allows stress differences to be evaluated and elastic constants to be determined with much reduced thermal noise. We then apply this approach across a diverse set of systems, spanning crystalline argon, ordered silicon as well as amorphous silicon, poly(methyl methacrylate), and cellulose derivatives.

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 / 1 minor

Summary. The paper claims that finite-temperature elastic constants can be computed with reduced thermal noise by applying a small strain to an equilibrated solid and evaluating stress differences between the strained system and an unstrained (or oppositely strained) reference system, both run with identical thermostatting schemes. Theoretical arguments and generic one-dimensional model tests are used to demonstrate the noise cancellation, followed by applications to crystalline argon, ordered silicon, amorphous silicon, PMMA, and cellulose derivatives.

Significance. If the noise cancellation remains unbiased, the approach offers a practical improvement over standard fluctuation-based methods for elastic constants at finite temperature, especially in noisy or disordered systems. The generalization from prior piezoelectric work, combined with explicit 1D tests and checks across five material classes, provides a solid foundation for broader use in materials simulations.

major comments (1)
  1. [Abstract, paragraph on method generalization] Abstract, paragraph on method generalization: the central assumption that identical thermostatting schemes ensure the stress difference Δσ = σ(ε) − σ(0) has both greatly reduced variance and an unbiased expectation value equal to the linear elastic response is load-bearing for the claim. Finite strain alters the Hamiltonian, so anharmonic terms can produce non-identical equilibrium measures and fluctuation spectra; this may leave residual variance (Cov term strictly less than sum of variances) and possible systematic bias in the finite-difference derivative. The 1D model demonstrations do not automatically guarantee clean cancellation for 3-D tensorial or disordered systems such as amorphous silicon.
minor comments (1)
  1. [Applications section] Applications section: full error-bar reporting and explicit convergence checks with respect to strain amplitude and simulation length would make the empirical validation across the five materials more quantitative and reproducible.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for highlighting the importance of rigorously justifying the noise-cancellation assumption. We address the major comment below and have revised the manuscript to strengthen the presentation of the theoretical limits and empirical validations.

read point-by-point responses

  1. Referee: [Abstract, paragraph on method generalization] Abstract, paragraph on method generalization: the central assumption that identical thermostatting schemes ensure the stress difference Δσ = σ(ε) − σ(0) has both greatly reduced variance and an unbiased expectation value equal to the linear elastic response is load-bearing for the claim. Finite strain alters the Hamiltonian, so anharmonic terms can produce non-identical equilibrium measures and fluctuation spectra; this may leave residual variance (Cov term strictly less than sum of variances) and possible systematic bias in the finite-difference derivative. The 1D model demonstrations do not automatically guarantee clean cancellation for 3-D tensorial or disordered systems such as amorphous silicon.

    Authors: We agree that the similarity of equilibrium measures under identical thermostatting is central and that finite strain modifies the Hamiltonian. Our theoretical section derives the cancellation explicitly in the small-strain limit, where the difference in measures contributes only at O(ε²) and does not bias the first-order response; the covariance reduction follows directly from the shared thermostat. The 1D models are presented as illustrative, but the manuscript already contains direct applications to three-dimensional disordered systems (amorphous silicon, PMMA, cellulose derivatives) that demonstrate practical noise reduction and consistency with independent calculations. To address the referee’s concern, we have added a dedicated paragraph in the discussion clarifying the small-strain regime, the order of neglected terms, and the empirical checks performed across material classes. We have also included a brief note on the method’s expected limitations for large strains or extreme anharmonicity. revision: yes

Circularity Check

0 steps flagged

Noise-cancellation procedure is a direct stress-difference computation with independent theoretical support

full rationale

The derivation begins from the standard definition of elastic constants via stress-strain response and introduces a computational protocol that correlates two trajectories through identical thermostatting after a small strain is applied. This correlation is shown to reduce variance in the difference Δσ by explicit calculation in the paper's theoretical section and in generic 1-D models; the reduction follows from the shared random forces and is not obtained by fitting any parameter to the target elastic constants themselves. Subsequent applications to argon, silicon, and polymers are presented as numerical illustrations rather than inputs that define the method. No self-citation chain is required to close the central argument, and the cited prior piezoelectric work supplies only the original noise-cancellation idea, not a uniqueness theorem or ansatz that forces the present result. The procedure therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The method rests on standard molecular-dynamics assumptions plus one tunable parameter (strain amplitude) chosen to remain in the linear regime; no new particles or forces are postulated.

free parameters (1)
  • strain amplitude
    Must be small enough for linear response yet large enough for measurable stress difference; value is chosen by the user.
axioms (1)
  • domain assumption The reference and strained systems remain in thermal equilibrium under identical thermostatting
    Invoked when claiming that thermal fluctuations cancel upon subtraction.

pith-pipeline@v0.9.0 · 5689 in / 1291 out tokens · 34526 ms · 2026-05-18T14:15:35.739861+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

42 extracted references · 42 canonical work pages

  1. [1]

    The value of k1d is 10 -1 10 0 10 1 10 2 1.0 2.0 3.0 4.0 low damping high damping all inst

    This expression corresponds to the curvature of the 6–12 LJ potential at its minimum, located at 2 1/6r0. The value of k1d is 10 -1 10 0 10 1 10 2 1.0 2.0 3.0 4.0 low damping high damping all inst. std's reduced by factor of 100 FIG. 2: Relaxation of the elastic–tensor estimator for a lin ear Lennard–Jones chain using the noise–cancellation techniq ue wit...

    work page

  2. [2]

    As discussed, the re- laxation function C(t) in Eq

    Interpretation In Figures 2 and 3, we demonstrated the influence of damping on the elastic modulus estimator for a one– dimensional Lennard–Jones chain. As discussed, the re- laxation function C(t) in Eq. (10) decays most rapidly when a mode is critically damped. In real systems, how- ever, the presence of a distribution of eigenfrequencies complicates thi...

    work page

  3. [3]

    At small damping times– where stiff vibrational modes are critically damped– they equilibrate rapidly

    Application Figure 7 presents results for two different samples, each simulated using a range of Langevin damping times τT. At small damping times– where stiff vibrational modes are critically damped– they equilibrate rapidly. How- ever, this also leads to a slower overall relaxation of C11. In contrast, larger damping times result in slower equili- bration...

    work page

  4. [4]

    Born and K

    M. Born and K. Huang, Dynamical Theory of Crystal Lattices. Oxford, UK: Oxford University Press, 1954

    work page 1954

  5. [5]

    Strain fluctuations and elastic constants,

    M. Parrinello and A. Rahman, “Strain fluctuations and elastic constants,” The Journal of Chemical Physics, vol. 76, no. 5, pp. 2662–2666, 1982

    work page 1982

  6. [6]

    Statistical ensembles and molecular dynamics studies of anisotropic solids,

    J. R. Ray and A. Rahman, “Statistical ensembles and molecular dynamics studies of anisotropic solids,” The Journal of Chemical Physics, vol. 80, no. 9, p. 4423–4428, 1984

    work page 1984

  7. [7]

    On the atomic theory of elasticity,

    K. Huang, “On the atomic theory of elasticity,” Proc. R. Soc. London, Ser. A, vol. 203, no. 1073, pp. 178–194, 1950

    work page 1950

  8. [8]

    Second-order elas- tic constants of a solid under stress,

    T. H. K. Barron and M. L. Klein, “Second-order elas- tic constants of a solid under stress,” Proc. Phys. Soc., vol. 85, no. 3, pp. 523–532, 1965

    work page 1965

  9. [9]

    E. B. Tadmor and R. E. Miller, Modeling Materials: Continuum, Atomistic and Multiscale Techniques. Cam- bridge University Press, 2011. Available online: https: //doi.org/10.1017/CBO9781139003582

    work page doi:10.1017/cbo9781139003582 2011

  10. [10]

    Analytic elastic coefficients in molec- ular calculations: Finite strain, nonaffine displacements, 13 and many-body interatomic potentials,

    J. Grießer, L. Fr´ erot, J. A. Oldenstaedt, M. H. M¨ user, and L. Pastewka, “Analytic elastic coefficients in molec- ular calculations: Finite strain, nonaffine displacements, 13 and many-body interatomic potentials,” Phys. Rev. Mater., vol. 7, p. 073603, 2023

    work page 2023

  11. [11]

    Isothermal elas- tic constants for argon. theory and monte carlo calcula- tions,

    D. Squire, A. Holt, and W. Hoover, “Isothermal elas- tic constants for argon. theory and monte carlo calcula- tions,” Physica, vol. 42, p. 388–397, May 1969

    work page 1969

  12. [12]

    Elastic constants and statistical ensembles in molecular dynamics,

    J. R. Ray, “Elastic constants and statistical ensembles in molecular dynamics,” Computer Physics Reports, vol. 8, p. 109–151, Aug. 1988

    work page 1988

  13. [13]

    Elastic constants of quantum solids by path integral simulations,

    P. Sch¨ offel and M. H. M¨ user, “Elastic constants of quantum solids by path integral simulations,” Physical Review B, vol. 63, May 2001

    work page 2001

  14. [14]

    Elastic moduli and thermal conductivity of quantum materials at finite temperature,

    D. A. Folkner, Z. Chen, G. Barbalinardo, F. Knoop, and D. Donadio, “Elastic moduli and thermal conductivity of quantum materials at finite temperature,” Journal of Applied Physics, vol. 136, Dec. 2024

    work page 2024

  15. [15]

    Piezoelectric coefficient s by molecular dynamics simulations in the constant stress ensemble: A case study of quartz,

    D. Herzbach and M. H. M¨ user, “Piezoelectric coefficient s by molecular dynamics simulations in the constant stress ensemble: A case study of quartz,” Computer Physics Communications, vol. 174, no. 1, pp. 17–23, 2006

    work page 2006

  16. [16]

    Generalized monte carlo perturbation algo- rithms for correlated sampling and a second-order tay- lor series approach,

    H. Rief, “Generalized monte carlo perturbation algo- rithms for correlated sampling and a second-order tay- lor series approach,” Annals of Nuclear Energy, vol. 11, p. 455–476, Jan. 1984

    work page 1984

  17. [17]

    Interatomic potentials: achievements and challenges,

    M. H. M¨ user, S. V. Sukhomlinov, and L. Pastewka, “Interatomic potentials: achievements and challenges,” Advances in Physics: X, vol. 8, Nov. 2022

    work page 2022

  18. [18]

    Polymorphic transitions in single crystals: A new molecular dynamics method,

    M. Parrinello and A. Rahman, “Polymorphic transitions in single crystals: A new molecular dynamics method,” J. Appl. Phys., vol. 52, no. 12, p. 7182, 1981

    work page 1981

  19. [19]

    Advanced langevin thermostats: Properties, ex- tensions to rheology, and a lean momentum-conserving approach,

    S. Agarwal, S. V. Sukhomlinov, M. Honecker, and M. H. M¨ user, “Advanced langevin thermostats: Properties, ex- tensions to rheology, and a lean momentum-conserving approach,” The Journal of Chemical Physics, vol. 163, p. 124110, 09 2025

    work page 2025

  20. [20]

    Lammps-a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales,

    A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolin- tineanu, 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, vol. 271, p. 108171, 2022

    work page 2022

  21. [21]

    Gromacs: High perfor- mance molecular simulations through multi-level paral- lelism from laptops to supercomputers,

    M. J. Abraham, T. Murtola, R. Schulz, S. P´ all, J. C. Smith, B. Hess, and E. Lindahl, “Gromacs: High perfor- mance molecular simulations through multi-level paral- lelism from laptops to supercomputers,” vol. 1,, pp. 19– 25, 2015

    work page 2015

  22. [22]

    Correlations in the motion of atoms in liq- uid argon,

    A. Rahman, “Correlations in the motion of atoms in liq- uid argon,” Phys. Rev., vol. 136, pp. A405–A411, 1964

    work page 1964

  23. [23]

    X-ray diffraction study of solid argon,

    C. S. Barrett and L. Meyer, “X-ray diffraction study of solid argon,” The Journal of Chemical Physics, vol. 41, no. 4, pp. 1078–1081, 1964

    work page 1964

  24. [24]

    A simple and ef- fective verlet-type algorithm for simulating langevin dy- namics,

    N. Grønbech-Jensen and O. Farago, “A simple and ef- fective verlet-type algorithm for simulating langevin dy- namics,” Molecular Physics, vol. 111, no. 8, pp. 983–991, 2013

    work page 2013

  25. [25]

    Canoni- cal sampling through velocity rescaling,

    G. Bussi, D. Donadio, and M. Parrinello, “Canoni- cal sampling through velocity rescaling,” Journal of Chemical Physics, vol. 126, p. 014101, 2007

    work page 2007

  26. [26]

    New empirical approach for the structure and energy of covalent systems,

    J. Tersoff, “New empirical approach for the structure and energy of covalent systems,” Phys. Rev. B, vol. 37, pp. 6991–7000, 1988

    work page 1988

  27. [27]

    De- velopment of bond-order potentials that can reproduce the elastic constants and melting point of silicon for classical molecular dynamics simulation,

    T. Kumagai, S. Izumi, S. Hara, and S. Sakai, “De- velopment of bond-order potentials that can reproduce the elastic constants and melting point of silicon for classical molecular dynamics simulation,” Computational Materials Science, vol. 39, no. 2, pp. 457–464, 2007

    work page 2007

  28. [28]

    Depleted depletion drives polymer swelling in poor sol- vent mixtures,

    D. Mukherji, C. M. Marques, T. St¨ uhn, and K. Kremer, “Depleted depletion drives polymer swelling in poor sol- vent mixtures,” Nature Communications, vol. 8, p. 1374, 2017

    work page 2017

  29. [29]

    Comparing simulated specific heat of liquid polymers and oligomers to experiments,

    H. Gao, T. P. W. Menzel, M. H. M¨ user, and D. Mukherji, “Comparing simulated specific heat of liquid polymers and oligomers to experiments,” Phys. Rev. Mater., vol. 5, p. 065605, 2021

    work page 2021

  30. [30]

    Computing the thermal transport coeffi- cient of neutral amorphous polymers using exact vibra- tional density of states: Comparison with experiments,

    D. Mukherji, “Computing the thermal transport coeffi- cient of neutral amorphous polymers using exact vibra- tional density of states: Comparison with experiments,” Physical Review Materials, vol. 8, p. 085601, 2024

    work page 2024

  31. [31]

    Particle mesh ewald: An N-log(N) method for ewald sums in large sys- tems,

    T. Darden, D. York, and L. Pedersen, “Particle mesh ewald: An N-log(N) method for ewald sums in large sys- tems,” The Journal of Chemical Physics, vol. 98, no. 12, pp. 10089–10092, 1993

    work page 1993

  32. [32]

    Heat flow in solvent-free, dense amorphous and semi-crystalline cellulose derivatives,

    D. Mukherji, T. E. de Oliveira, D. G. Cahill, and M. M¨ uller, “Heat flow in solvent-free, dense amorphous and semi-crystalline cellulose derivatives,” submitted on 18th August, 2025

    work page 2025

  33. [33]

    Glycam06: A generalizable biomolecu- lar force field. carbohydrates,

    K. N. Kirschner, A. B. Yongye, S. M. Tschampel, J. Gonz´ alez-Outeiri˜ no, C. R. Daniels, B. L. Foley, and R. J. Woods, “Glycam06: A generalizable biomolecu- lar force field. carbohydrates,” Journal of Computational Chemistry, vol. 29, no. 4, pp. 622–655, 2008

    work page 2008

  34. [34]

    Charmm-gui gly- can modeler for modeling and simulation of carbohy- drates and glycoconjugates,

    S.-J. Park, J. Lee, Y. Qi, N. R. Kern, H. S. Lee, S. Jo, I. Joung, K. Joo, J. Lee, and W. Im, “Charmm-gui gly- can modeler for modeling and simulation of carbohy- drates and glycoconjugates,” Glycobiology, vol. 29, no. 4, pp. 320–331, 2019

    work page 2019

  35. [35]

    Electronic supplementary information– document num - ber to be added by the editors,

    “Electronic supplementary information– document num - ber to be added by the editors,”

    work page

  36. [36]

    The dynamics of supercooled silica: acoustic modes and boson peak,

    J. Horbach, W. Kob, and K. Binder, “The dynamics of supercooled silica: acoustic modes and boson peak,” Journal of Non-Crystalline Solids, vol. 235-237, pp. 320– 324, 1998

    work page 1998

  37. [37]

    Elastic constants of quan- tum solids by path integral simulations,

    P. Sch¨ offel and M. H. M¨ user, “Elastic constants of quan- tum solids by path integral simulations,” Phys. Rev. B, vol. 63, p. 224108, 2001

    work page 2001

  38. [38]

    The elastic constants of solid argon determined by stim- ulated brillouin scattering,

    H. Meixner, P. Leiderer, P. Berberich, and E. L¨ uscher, “The elastic constants of solid argon determined by stim- ulated brillouin scattering,” Physics Letters A, vol. 40, no. 3, pp. 257–258, 1972

    work page 1972

  39. [39]

    Electronic effects in the elastic constants of n-type silicon,

    J. J. Hall, “Electronic effects in the elastic constants of n-type silicon,” Phys. Rev., vol. 161, pp. 756–761, 1967

    work page 1967

  40. [40]

    Elastic constants of defected and amorphous silicon with the environment-dependent interatomic potential,

    C. L. Allred, X. Yuan, M. Z. Bazant, and L. W. Hobbs, “Elastic constants of defected and amorphous silicon with the environment-dependent interatomic potential,” Phys. Rev. B, vol. 70, p. 134113, 2004

    work page 2004

  41. [41]

    Thermal conductivity, heat capacity, and elas- tic constants of water-soluble polymers and polymer blends,

    X. Xie, D. Li, T. Tsai, J. Liu, P. V. Braun, and D. G. Cahill, “Thermal conductivity, heat capacity, and elas- tic constants of water-soluble polymers and polymer blends,” Macromolecules, vol. 49, pp. 972–978, 2016

    work page 2016

  42. [42]

    Thermal conductivity of commodity polymers under high pressures: Macroscopic heat flow, monomer–level energy transfer, and quantum effects,

    O. H. M. de Alencar, J. Wu, M. M¨ uller, and D. Mukherji, “Thermal conductivity of commodity polymers under high pressures: Macroscopic heat flow, monomer–level energy transfer, and quantum effects,” to be submitted, 2025

    work page 2025