shc.out
This file contains the non-equilibrium virial-velocity correlation function \(K(t)\) and the spectral heat current (SHC) \(J_q(\omega)\), in a given direction, for a group of atoms, as defined in Eq. (18) and the left part of Eq. (20) of [Fan2019]. It is generated when invoking the compute_shc keyword.
File format
For each run, there are 3 columns and 2*Nc-1 + num_omega rows.
Here, Nc is the number of correlation steps and num_omega is the number of frequency points.
In the first 2*Nc-1 rows:
column 1: correlation time \(t\) from negative to positive, in units of ps
column 2: \(K^{\rm in}(t)\) in units of Å eV/ps
column 3: \(K^{\rm out}(t)\) in units of Å eV/ps
\(K^{\rm in}(t) + K^{\rm out}(t) = K(t)\) is exactly the expression in Eq. (18) of [Fan2019]. The in-out decomposition follows the definition in [Fan2017], which is useful for 2D materials but is not necessary for 3D materials.
In the next num_omega rows:
column 1: angular frequency \(\omega\) in units of THz
column 2: \(J_q^{\rm in}(\omega)\) in units of Å eV/ps/THz
column 3: \(J_q^{\rm out}(\omega)\) in units of Å eV/ps/THz
\(J_q^{\rm in}(\omega) + J_q^{\rm out}(\omega) = J_q(\omega)\) is exactly the left expression in Eq. (20) of [Fan2019].
Only the potential part of the heat current has been included.
If :attr: ‘group_id’ is -1, then the file follows the above rules and will contain \(K(t)\) and \(J_q(\omega)\) for each group id except for group id 0. And the contents of the :attr: ‘group_id’ are arranged from smallest to largest.