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
output_interval
is the interval to output the onsager coefficients.
Fe_x
is the
Fe_y
is the
Fe_z
is the
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
direction and has a magnitude of Å-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