Orsolya Gereben, László Pusztai
Institute for Solid State Physics and Optics, Wigner Research Centre for Physics, Hungarian Academy of Sciences, PO Box 49, H-1525 Budapest, Hungary.
Journal of computational chemistry 2012 Nov 5An approach has been devised and tested for preserving the molecular dynamics molecular geometry taking into account energetic considerations during Reverse Monte Carlo (RMC) modeling. Instead of the commonly used fixed neighbor constraints, where molecules are held together by constraining distance ranges available for the specified atom pairs, here molecules are kept together via bond, angle, and dihedral potential energies. The scaled total potential energy contributes to the measure of the goodness-of-fit, thus, the atoms can be prevented from drifting apart. In some of the calculations (Lennard-Jones and Coulombic) nonbonding potentials were also applied. The algorithm was successfully tested for the X-ray structure factor-based structure study of liquid dimethyl trisulfide, for which material now significantly more sensible results have been obtained than during previous attempts via any earlier version of RMC modeling. It is envisaged that structural modeling of a large class of materials, primarily liquids and amorphous solids containing molecules of up to about 100 atoms, will make use of the new code in the near future. Copyright © 2012 Wiley Periodicals, Inc.
Orsolya Gereben, László Pusztai. RMC_POT: a computer code for reverse Monte Carlo modeling the structure of disordered systems containing molecules of arbitrary complexity. Journal of computational chemistry. 2012 Nov 5;33(29):2285-91
PMID: 22782785
View Full Text