Tsallis dynamics using the Nosé-Hoover approach

On the basis of the Nosé-Hoover method, we developed a deterministic algorithm that produces an arbitrary probability density. An ordinary differential equation in the algorithm can realize the Tsallis distribution density. The Tsallis distribution has been considered a candidate of a distribution that represents a physical system in a heat bath. The Tsallis distribution density employed in this algorithm is defined using a full energy function form $E(x,p),$ along with the Tsallis index $q>~1.$ Using the current equation, numerical simulations were performed for simple systems and the Tsallis distributions were observed.