`compute_hnemdec`

This keyword is used to calculate the multicomponent system thermal conductivity using the homogeneous non-equilibrium molecular dynamics Evans-Cummings algorithm (HNEMDEC) method. The results are written to the onsager.out output file.

## Syntax

```
compute_hnemdec <driving_force> <output_interval> <Fe_x> <Fe_y> <Fe_z>
```

`driving_force`

determines which type of driving force to use. It could be zero or non zero positive integer, which means thermal driving force or diffusive driving force.
In the case of zero, heat flux is equilvalent to dissipative flux, with driving force \(F_e\) in the units of Å^{-1}.
For a non zero positive integer such as \(i\), a momentum flux of the \(ith\) element in the order of the first line of nep.txt is produced as dissipative flux, with driving force \(F_e\) in the units of eV/Å.

`output_interval`

is the interval to output the onsager coefficients.

`Fe_x`

is the \(x\) direction component of the external driving force \(F_e\).

`Fe_y`

is the \(y\) direction component of the external driving force \(F_e\).

`Fe_z`

is the \(z\) direction component of the external driving force \(F_e\).

## Examples

### Example 1

```
compute_hnemdec 0 1000 0.00001 0 0
```

This means that

you want to calculate the onsager coefficients using the HNEMDEC method with heat flux as dissipative flux;

the onsager coefficients will be averaged and output every 1000 steps (the heat current and momentum current is sampled for every step);

the external driving force is along the \(x\) direction and has a magnitude of \(10^{-5}\) Å

^{-1}.

Note that one should control the temperature when using this keyword. Otherwise, the system will be heated up by the external driving force.

**Important:**
For this purpose, the Nose-Hoover chain thermostat is recommended.
The Langevin thermostat cannot be used for this purpose because it will affect the dynamics of the system.

### Example 2

```
compute_hnemdec 2 1000 0 0.00001 0
```

The external driving force is diffusive driving force that will produced a momentum flux of the second element as dissipative flux and has a magnitude of \(10^{-5}\) eV/Å. The force is applied along the \(y\) direction.