Global Thermodynamic Properties of Complex Spin Systems Calculated from Density of States and Indirectly by Thermodynamic Integration Method

Institute of Physics, Faculty of Science, Pavol Jozef Šafárik University in Košice, Park Angelinum 9, 040 01 Košice, Slovakia

^★ e-mail: marek.semjan@student.upjs.sk
^★★ e-mail: milan.zukovic@upjs.sk

Published online: 20 January 2020

Abstract

Evaluation of global thermodynamic properties such as the entropy or the free energy of complex systems featuring a high degree of frustration or disorder is often desirable. Nevertheless, they cannot be measured directly in standard Monte Carlo simulation. Therefore, they are either evaluated indirectly from the directly measured quantities, for example by the thermodynamic integration method (TIM), or by applying more sophisticated simulation methods, such as the Wang-Landau (WL) algorithm, which can directly sample density of states. In the present investigation we compare the performance of the WL and TIM methods for the calculation of the entropy of an Ising antiferromagnetic system on a Kagome lattice – a typical example of a complex spin system with high geometrical frustration resulting in a non-zero residual entropy the value of which is exactly known. It is found that the easier to implement TIM can yield results of comparable accuracy with that of the more involved WL method.