Computer simulations of dynamic response of ferrofluids on an alternating magnetic field with high amplitude

The response of ferrofluids to a high-amplitude AC magnetic field is important for several applications including magnetic hyperthermia and biodetection. In computer simulations of the dynamic susceptibility of a ferrofluid outside the linear response region, there are several problems associated with the fact that an increase in the frequency of the AC field leads to the appearance of additional computational errors, which can even lead to unphysical results. In this article, we study the dependence of the computational error arising in the computer simulation of the dynamic susceptibility on the input parameters of the numerical algorithm: the length of the time step, the total number of computer simulation periods, and averaging period. Computer simulation is carried out using the Langevin dynamics method and takes Brownian rotational relaxation of magnetic particles and interparticle interactions into account. The reference theory [Yoshida T.; Enpuku K. Jap. J. Ap. Phys. 2009] is used to estimate computational error. As a result, we give practical recommendations for choosing the optimal input parameters of the numerical algorithm, which make it possible to obtain reliable results of the dynamic susceptibility of a ferrofluid in a high-amplitude AC field in a wide frequency range.