In practice, as the external field is non-conservative (i.e. as it is not the gradient of a potential), it will result in a constant dissipation in the system at steady state.
Hence, the use of color fields must be combined with the use of a thermostat.
In what follows, we will make use of the Nosé-Hoover thermostat, as it is a global thermostat, and it conserves linear momentum [1]}.
