corner
corner

Phys. Rev. B 67, 064422 (2003) [10 pages]

Monte Carlo technique with a quantified time step: Application to the motion of magnetic moments

Download: PDF (306 kB) Buy this article Export: BibTeX or EndNote (RIS)

O. Chubykalo
Instituto de Ciencia de Materiales de Madrid, CSIC, Cantoblanco, E-28049 Madrid, Spain

U. Nowak
Theoretische Tieftemperaturphysik, Gerhard-Mercator-Universität-Duisburg, D-47048 Duisburg, Germany

R. Smirnov-Rueda
Departamento de Matematica Aplicada, Facultad de Ciencias Matematicas, Universidad Complutense de Madrid, 28040 Madrid, Spain

M. A. Wongsam
Department of Physics, University of Durham, South Road, Durham, D1 3LE, United Kingdom

R. W. Chantrell
Seagate Research, River Parks Commons, 2403 Sydney Street, Pittsburgh, Pennsylvania 15203-2116

J. M. Gonzalez
Instituto de Ciencia de Materiales de Madrid, CSIC, Cantoblanco, E-28049 Madrid, Spain
Instituto de Magnetismo Aplicado, RENFE/UCM, Las Rozas, Madrid, Spain

Received 21 March 2002; published 28 February 2003

The viability of the time quantified Metropolis Monte Carlo technique to describe the dynamics of magnetic systems is discussed. Similar to standard Brownian motion, the method is introduced basing on the comparison between the Monte Carlo trial step and the mean squared deviation of the direction of the magnetic moment. The Brownian dynamics approach to the time evolution of a magnetic moment is investigated and expressions for the mean square deviations are obtained. However, the principle difference between the standard Brownian motion and the magnetic moments dynamics is the presence of the spin precession which constitutes the reversible part of the dynamics. Although some part of the precession contributes to the diffusion coefficient, it also gives rise to athermal, energy conserving motion which cannot be taken into account by Monte Carlo methods. It is found that the stochastic motion of a magnetic moment falls into one of two possible regimes: (i) precession dominated motion, (ii) nonprecessional motion, according to the value of the damping constant and anisotropy strength and orientation. Simple expressions for the diffusion coefficient can be obtained in both cases for diffusion dominated motion, i.e., where the athermal precessional contribution can be neglected. These simple expressions are used to convert the Monte Carlo steps to real time units. The switching time for magnetic particles obtained by the Monte Carlo with time quantification is compared with the numerical integration of the Landau-Lifshitz-Gilbert equations with thermal field contribution and with some well known asymptotic formulas.

© 2003 The American Physical Society

URL:
http://link.aps.org/doi/10.1103/PhysRevB.67.064422
DOI:
10.1103/PhysRevB.67.064422
PACS:
75.40.Gb, 75.40.Mg, 75.50.Tt