Fast recall of state-history in kinetic Monte Carlo simulations utilizing the Zobrist key

被引:11
作者
Mason, DR
Hudson, TS
Sutton, AP
机构
[1] Univ Oxford, Dept Mat, Oxford OX1 3PH, England
[2] Aalto Univ, Lab Computat Engn, FIN-02015 Espoo, Finland
基金
欧盟地平线“2020”; 英国工程与自然科学研究理事会;
关键词
vacancies; vacancy diffusion; kinetic Monte Carlo; residence time; hash function; on-the-fly kMC;
D O I
10.1016/j.cpc.2004.09.007
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
We present a novel application of the Zobrist hashing method, known in the computer chess literature, to simulation of diffusional phase transformations in metal alloys. A history of previously visited states can be easily maintained, allowing very fast lookup of energies and transition rates calculated earlier in the simulation. The method has been applied to the simulation of a Fe-1at.%Cu system, with simple potentials and a transition rate for diffusional events approximated from the difference in internal energy between trial states. In this simple model at temperatures of 1073 K we find that 61.2% of the states considered during the simulation have been seen previously, and that this proportion rises to 85.1% at 773 K and even to 99.9% at 373 K. Rapid recall of these states reduces the computational time taken for the same sequence of atom-vacancy exchange moves by a factor of 6.3 at 773 K rising to over 100 at 373 K. We suggest that a similar speedup factor will be found using more sophisticated models of diffusion and that the method can, with small modifications, be applied to a wide range of kinetic Monte Carlo simulations of atomistic diffusion processes. (C) 2004 Elsevier B.V. All rights reserved.
引用
收藏
页码:37 / 48
页数:12
相关论文
共 19 条
[1]   Computer simulation of point defect properties in dilute Fe-Cu alloy using a many-body interatomic potential [J].
Ackland, GJ ;
Bacon, DJ ;
Calder, AF ;
Harry, T .
PHILOSOPHICAL MAGAZINE A-PHYSICS OF CONDENSED MATTER STRUCTURE DEFECTS AND MECHANICAL PROPERTIES, 1997, 75 (03) :713-732
[2]  
[Anonymous], COMPUT AUTOM
[3]   NEW ALGORITHM FOR MONTE-CARLO SIMULATION OF ISING SPIN SYSTEMS [J].
BORTZ, AB ;
KALOS, MH ;
LEBOWITZ, JL .
JOURNAL OF COMPUTATIONAL PHYSICS, 1975, 17 (01) :10-18
[4]   Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points [J].
Henkelman, G ;
Jónsson, H .
JOURNAL OF CHEMICAL PHYSICS, 2000, 113 (22) :9978-9985
[5]   A dimer method for finding saddle points on high dimensional potential surfaces using only first derivatives [J].
Henkelman, G ;
Jónsson, H .
JOURNAL OF CHEMICAL PHYSICS, 1999, 111 (15) :7010-7022
[6]   Long time scale kinetic Monte Carlo simulations without lattice approximation and predefined event table [J].
Henkelman, G ;
Jónsson, H .
JOURNAL OF CHEMICAL PHYSICS, 2001, 115 (21) :9657-9666
[7]   Computer simulation of phase decomposition process generating precipitates harder than matrix [J].
Ikeda, H ;
Matsuda, H .
MATERIALS TRANSACTIONS JIM, 1996, 37 (08) :1413-1421
[8]   Theoretical studies of atomic-scale processes relevant to crystal growth [J].
Jónsson, H .
ANNUAL REVIEW OF PHYSICAL CHEMISTRY, 2000, 51 :623-653
[9]  
Knuth D. E., 1973, The Art of Computer Programming Volume 3, Sorting and Searching, VIII
[10]  
MARSLAND TA, 1992, ENCY ARTIFICIAL INTE, P224