Chaotic behavior of the lattice Yang-Mills on CUDA

The Yang-Mills fields plays important role in the strong interaction, which describes the quark gluon plasma. The non-Abelian gauge theory provides the theoretical background understanding of this topic. The real time evolution of the classical fields is derived by the Hamiltonian for SU(2) gauge field tensor. The microcanonical equations of motion is solved on 3 dimensional lattice and chaotic dynamics was searched by the monodromy matrix. The entropy-energy relation was presented by Kolmogorov-Sinai entropy. We used block Hessenberg reduction to compute the eigenvalues of the current matrix. While the purely CPU based algorithm can handle effectively only a small amount of values, the GPUs provide enough performance to give more computing power to solve the problem.