Published January 23, 2025 | Version v1
Dataset Open

Transport coefficients from equilibrium molecular dynamics

  • 1. COSMO—Laboratory of Computational Science and Modelling, IMX, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
  • 2. SISSA—Scuola Internazionale Superiore di Studi Avanzati, 34136 Trieste, Italy
  • 3. Department of Physics, Informatics and Mathematics, Università degli Studi di Modena e Reggio Emilia, 41125 Modena, Italy
  • 4. CNR-Nano S3—Istituto Nanoscienze, 41125 Modena, Italy

* Contact person

Description

The determination of transport coefficients through the time-honoured Green-Kubo theory of linear response and equilibrium molecular dynamics requires significantly longer simulation times than those of equilibrium properties, while being further hindered by the lack of well-established data-analysis techniques to evaluate the statistical accuracy of the results. Leveraging recent advances in the spectral analysis of the current time series associated to molecular trajectories, we introduce a new method to estimate the full (diagonal as well as off-diagonal) Onsager matrix of transport coefficients from a single statistical model. This approach, based on the knowledge of the statistical distribution of the Onsager-matrix samples in the frequency domain, unifies the evaluation of diagonal (conductivities and viscosities) and off-diagonal (e.g., thermoelectric) transport coefficients within a comprehensive framework, significantly improving the reliability of transport coefficient estimation for materials ranging from molten salts to solid-state electrolytes. We validate the accuracy of this method against existing approaches using benchmark data on molten cesium fluoride and liquid water, and conclude our presentation with the computation of various transport coefficients of the Li₃PS₄ solid-state electrolyte.

Files

File preview

All files

Files (405.0 MiB)

Name Size
md5:6cfe9681bf170b917e6e30a5e471cf21
1.1 KiB Preview Download
md5:b46639efbb0fc08cab67a8aa2c58ed38
1.6 MiB Preview Download
md5:82bfa8a603d9fe586d50a9c9761ad2d9
343 Bytes Download
md5:66ab73f02db7f98648ebc9320aef58be
1.3 MiB Preview Download
md5:e2cd7273b19f383788a52ef7f7210b8a
127.1 MiB Preview Download
md5:1a73c8dbb538e03919d3479a04427833
7.4 KiB Preview Download
md5:bc7e84c614b4334e0f6b480ad524e6f7
3.4 KiB Preview Download
md5:79004101334e942eec0551ad6966fd12
13.0 KiB Preview Download
md5:57abb5afba3f5995f26cb9ded689ecbe
8.7 KiB Preview Download
md5:50e803f60f78fdd440bccd7eee7fa3d5
48.5 KiB Preview Download
md5:0b84c469df29cf9f143492897ca00168
275.0 MiB Preview Download

References

Journal reference (Paper in which the method is described)
P. Pegolo, E. Drigo, F. Grasselli, S. Baroni, accepted at the Journal of Chemical Physics (2025)