This file contains some components of the running thermal conductivity (RTC) tensor from the homogeneous nonequilibrium molecular dynamics (HNEMD) method. It is generated when invoking the compute_hnemd keyword.

File format

If the driving force is in the \(\mu\) (\(\mu\) can be \(x\), \(y\), or \(z\)) direction, this file reads:

  • column 1: \(\kappa_{\mu x}^{\text{in}}(t)\) (in units of W/mK)

  • column 2: \(\kappa_{\mu x}^{\text{out}}(t)\) (in units of W/mK)

  • column 3: \(\kappa_{\mu y}^{\text{in}}(t)\) (in units of W/mK)

  • column 4: \(\kappa_{\mu y}^{\text{out}}(t)\) (in units of W/mK)

  • column 5: \(\kappa_{\mu z}^{\text{tot}}(t)\) (in units of W/mK)

Both \(\kappa_{\mu x}(t)\) and \(\kappa_{\mu y}(t)\) have been decomposed into comtributions from in-plane (hence superscript in) and out-of-plane (hence superscript out) vibrational modes, as described in [Fan2019]. This decomposition is useful for 2D (or layered) materials but is not necessary for 3D materials. For 3D materials, one can sum up some columns to get the conventional data. That is:

\[\begin{split}\kappa_{\mu x}^{\text{tot}}(t) &= \kappa_{\mu x}^{\text{in}}(t) + \kappa_{\mu x}^{\text{out}}(t) \\ \kappa_{\mu y}^{\text{tot}}(t) &= \kappa_{\mu y}^{\text{in}}(t) + \kappa_{\mu y}^{\text{out}}(t).\end{split}\]

Only the potential part of the heat current has been considered. To simulation systems in which the convective heat current is important, one would have to modify the source code.