Transport coefficients from equilibrium molecular dynamics
Creators
- 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
files_description.md
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)