逆向蒙特卡罗方法 (
英语:Reverse Monte Carlo method ),是标准Metropolis-Hastings算法的变体,用于解决逆向问题,即调节模型使其参数与实验数据达到最大的一致性。在科学和数学的分支中经常遇到逆向问题,然而这一方法可能更广泛地应用于凝聚态物理学和固体化学。
在凝聚态科学中的应用 [ 编辑]
基本方法 [ 编辑]
这一方法经常出现在凝聚态物理学中,生成与实验数据相对应的原子结构模型或是某种限制条件下的物态。
初始的结构是周期性边界条件下N个原子组成的晶胞。基于这一构型可计算一个或一个以上的可测量参数。常用的可测量参数包括对分布函数及其傅里叶变换形式,后者可由中子或x射线衍射直接获得实验数值。其余还有晶体材料的布拉格衍射参数和扩展x射线吸收精细结构参数。实验值与模拟值的比较由以下函数形式量化
χ 2 = ∑ (y obs − y calc )2 / σ 2
其中
y obs 和
y calc 分别为观测值(实验值)和计算值(模拟值),
σ 是测量精度。
根据观测值,随机选取原子进行随机运动。将这一过程重复迭代,使得
χ 2 变大或变小。当
χ 2 最小时,可认为体系处于平衡状态。
应用 [ 编辑]
逆向蒙特卡罗方法由McGreevy和Pusztai于1988年提出用于凝聚态问题的研究
[1] 。多年来这一直是液态和非晶态材料研究中,获得结构模型的唯一方法。
参考资料 [ 编辑]
^ RL McGreevy and L Pusztai, Reverse Monte Carlo Simulation: A New Technique for the Determination of Disordered Structures, Molecular Simulation 1, 359–367, 1988 [1]