Reverse Monte Carlo: Difference between revisions

From SklogWiki
Jump to navigation Jump to search
No edit summary
Line 27: Line 27:


== References ==
== References ==
#R.L.McGreevy and L. Pusztai, ''Mol. Simulation,'' '''1''' 359-367 (1988)
#[http://dx.doi.org/10.1080/08927028808080958 R.L.McGreevy and L. Pusztai, "Reverse Monte Carlo Simulation: A New Technique for the Determination of Disordered Structures", Molecular  Simulation, '''1''' pp. 359-367 (1988)]

Revision as of 18:52, 19 February 2007

Reverse Monte Carlo (RMC) [1] is a variation of the standard Metropolis Monte Carlo (MMC) method. It is used to produce a 3 dimensional atomic model that fits a set of measurements (Neutron-, X-ray-diffraction, EXAFS etc.). In addition to measured data a number of constraints based on prior knowledge of the system (like chemocal bonds etc.) can be applied. Some examples are:

  1. Closest approach between atoms (hard sphere potential)
  2. Coordination numbers.
  3. Angels in triplets of atoms.

The algorithm for RMC can be written:

1. Start with a configuration of atoms with periodic boundary conditions. This can be a random or a crystalline configuration from a different simulation or model. 2. Calculate the total radial distribution function for this old configuration. 3. Transform to the total structure factor:

where Q is the momentum transfer and the number density. 4. Calculate the difference between the measured structure factor and the one calculated from the configuration :

this sum is taken over all experimental points is the experimental error. 5. Select and move one atom at random and calculate the new distribution function, structure factor and: 6. If accept the move and let the new configuration become the old. If then the move is accepted with probability otherwiase rejected. 7. repeat from step 5. When have reached an equilibrium the configuration is saved and can be analysed.


References

  1. R.L.McGreevy and L. Pusztai, "Reverse Monte Carlo Simulation: A New Technique for the Determination of Disordered Structures", Molecular Simulation, 1 pp. 359-367 (1988)