Reverse Monte Carlo Simulation: A New Technique for the Determination of Disordered Structures

Abstract We have developed a new technique, based on the standard Monte Carlo simulation method with Markov chain sampling, in which a set of three dimensional particle configurations are generated that are consistent with the experimentally measured structure factor. A(Q), and radial distribution function, g(r), of a liquid or other disordered system. Consistency is determined by a standard χ2 test using the experimental errors. No input potential is required, we present initial results for liquid argon. Since the technique can work directly from the structure factor it promises to be useful for modelling the structures of glasses or amorphous materials. It also has other advantages in multicomponent systems and as a tool for experimental data analysis.

[1]  Roy Kaplow,et al.  Atomic Arrangement in Vitreous Selenium , 1968 .

[2]  I. R. Mcdonald,et al.  Theory of simple liquids , 1998 .

[3]  J. Kern,et al.  The determination of atomic interaction energies in solid solutions from short range order coefficients—an inverse monte-carlo method , 1987 .

[4]  N. H. March,et al.  Ion-ion oscillatory potentials in liquid metals , 1964, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[5]  Denis Weaire,et al.  Modeling Tetrahedrally Bonded Random Networks by Computer , 1987 .

[6]  R. Howe,et al.  The structure of molten sodium chloride , 1975 .

[7]  R. Mcgreevy,et al.  The structure of molten CsCl , 1985 .

[8]  B. Lebech Neutron scattering data analysis: edited by M.W. Johnson, Institute of Physics Conference Series No. 81 (Institute of Physics and Adam Hilger, Bristol and Boston, 1986) pp. vii + 178, £19.00, ISBN 0 85498 1721 , 1988 .

[9]  Max Born,et al.  A general kinetic theory of liquids , 1949 .

[10]  R. G. Wenzel,et al.  Structure Factor and Radial Distribution Function for Liquid Argon at 85 °K , 1973 .

[11]  N. Cowlam,et al.  Interatomic Pair Potentials for Metallic Alloy Glasses , 1987 .

[12]  Applications of the inverse Monte Carlo method in photon beam physics , 1987 .

[13]  Weis,et al.  Pair interaction from structural data for dense classical liquids. , 1985, Physical review letters.

[14]  H. Rietveld The crystal structure of some alkaline earth metal uranates of the type M3UO6 , 1966 .

[15]  A. Baranyai,et al.  Statistical geometry of molten alkali halides , 1986 .