Operator splitting algorithm for isokinetic SLLOD molecular dynamics.

We apply an operator splitting method to develop a simulation algorithm that has complete analytical solutions for the Gaussian thermostated SLLOD equations of motion [D. J. Evans and G.

The conjugate-pairing rule for non-Hamiltonian systems.

In systems that satisfy the Conjugate Pairing Rule (CPR), the spectrum of Lyapunov exponents is symmetric. The sum of each conjugate pair of exponents is identical. Since in dissipative systems the sum of all the exponents is the entropy production divided by Boltzmann's constant, the calculation of transport coefficients from the Lyapunov exponents is greatly simplified in systems that satisfy CPR.