SpecSwap-RMC: a novel reverse Monte Carlo approach using a discrete set of local configurations and pre-computed properties.

We present a novel approach to reverse Monte Carlo (RMC) modeling, SpecSwap-RMC, specifically applicable to structure modeling based on properties that require significant computer time to evaluate. In this approach pre-computed property data from a discrete set of local configurations are used and the configuration space is expressed in this basis. Atomistic moves are replaced with swap moves of contributions to a sample set representing the state of the simulated system. We demonstrate the approach by fitting jointly and separately the EXAFS signal and x-ray absorption spectrum (XAS) of ice Ih using a SpecSwap sample set of 80 configurations from a library of 1382 local structures with associated pre-computed spectra. As an additional demonstration we compare SpecSwap and FEFFIT fits of EXAFS data on crystalline copper, finding excellent agreement. SpecSwap-RMC thus extends RMC structure modeling to any property that can be computed from a structure irrespective of computational expense, but at the cost of a reduced configuration space. The method is general enough that it can be applied to any sets of computed properties, not necessarily limited to structure determination.

[1]  T. Straatsma,et al.  THE MISSING TERM IN EFFECTIVE PAIR POTENTIALS , 1987 .

[2]  P. Wernet,et al.  Surface structure of thin ice films , 2004 .

[3]  P. Wernet,et al.  Erratum: “Nearest-neighbor oxygen distances in liquid water and ice observed by x-ray Raman based extended x-ray absorption fine structure” [J. Chem. Phys. 127, 174504 (2007)] , 2008 .

[4]  A. K. Soper An asymmetric model for water structure , 2005 .

[5]  P. Wernet,et al.  The Structure of the First Coordination Shell in Liquid Water , 2004, Science.

[6]  E. Stern,et al.  Analysis of multiple-scattering XAFS data using theoretical standards , 1995 .

[7]  T. Pennanen,et al.  Nuclear magnetic shielding and quadrupole coupling tensors in liquid water: a combined molecular dynamics simulation and quantum chemical study. , 2004, Journal of the American Chemical Society.

[8]  A. Ankudinov,et al.  REAL-SPACE MULTIPLE-SCATTERING CALCULATION AND INTERPRETATION OF X-RAY-ABSORPTION NEAR-EDGE STRUCTURE , 1998 .

[9]  M Newville,et al.  IFEFFIT: interactive XAFS analysis and FEFF fitting. , 2001, Journal of synchrotron radiation.

[10]  T. Tyson,et al.  “Ab-initio” modelling of x-ray absorption spectra , 1991 .

[11]  G. Tóth,et al.  Direct determination of two-body potentials from measured pair structures , 2000 .

[12]  H. Ågren,et al.  Detailed study of pyridine at the C 1s and N 1s ionization thresholds: The influence of the vibrational fine structure , 2001 .

[13]  Alan K. Soper,et al.  Empirical potential Monte Carlo simulation of fluid structure , 1996 .

[14]  I. Hertel,et al.  Hydrogen bonds in liquid water studied by photoelectron spectroscopy. , 2007, The Journal of chemical physics.

[15]  M. Odelius,et al.  Are recent water models obtained by fitting diffraction data consistent with infrared/Raman and x-ray absorption spectra? , 2006, The Journal of chemical physics.

[16]  P. Wernet,et al.  Nearest-neighbor oxygen distances in liquid water and ice observed by x-ray Raman based extended x-ray absorption fine structure. , 2007, The Journal of chemical physics.

[17]  H. Ohashi,et al.  High resolution X-ray emission spectroscopy of liquid water : The observation of two structural motifs , 2008 .

[18]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[19]  A. Trapananti,et al.  Reverse Monte Carlo refinement of molecular and condensed systems by x-ray absorption spectroscopy , 2005 .

[20]  A. Lyubartsev,et al.  Assessing the electric-field approximation to IR and Raman spectra of dilute HOD in D2O. , 2009, The Journal of chemical physics.

[21]  G. Duclos New York 1987 , 2000 .

[22]  D. Nordlund,et al.  Electronic structure effects in liquid water studied by photoelectron spectroscopy and density functional theory , 2008 .

[23]  A. Lyubartsev,et al.  Theoretical Approximations to X-Ray Absorption Spectroscopy of Liquid Water and Ice , 2010 .

[24]  R. Mcgreevy,et al.  Reverse Monte Carlo simulation for the analysis of EXAFS data , 1990 .

[25]  Hans Ågren,et al.  Calculations of near-edge x-ray-absorption spectra of gas-phase and chemisorbed molecules by means of density-functional and transition-potential theory , 1998 .

[26]  András Baranyai,et al.  Conceptual and technical improvement of the reverse Monte Carlo algorithm , 1997 .

[27]  R. Mcgreevy,et al.  Reverse Monte Carlo modelling , 2001 .

[28]  Ankudinov,et al.  Multiple-scattering calculations of x-ray-absorption spectra. , 1995, Physical review. B, Condensed matter.

[29]  S. Adams,et al.  Bond valence analysis of reverse Monte Carlo produced structural models; a way to understand ion conduction in glasses , 2005 .

[30]  L. Pettersson,et al.  Intra- and intermolecular effects in the Compton profile of water , 2006 .

[31]  A. Soper Partial structure factors from disordered materials diffraction data: An approach using empirical potential structure refinement , 2005 .

[32]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[33]  R. L. McGreevy,et al.  Reverse Monte Carlo Simulation: A New Technique for the Determination of Disordered Structures , 1988 .

[34]  M. Odelius,et al.  Diffraction and IR/Raman data do not prove tetrahedral water. , 2008, The Journal of chemical physics.

[35]  A. Trapananti,et al.  Is there icosahedral ordering in liquid and undercooled metals? , 2003, Physical review letters.

[36]  P. Jóvári,et al.  A new version of the RMC++ Reverse Monte Carlo programme, aimed at investigating the structure of covalent glasses , 2007 .