Vol.18 No.1 郭春泰：Reverse Monte Carlo计算机模拟 .31· 通过RMC,求出3个偏RDF,这对无法测定偏结构因子的体系是十分有益的.图2 表明MgBi2具有典型的离子化合物（例NaCl,ZnCl)结构形式.Mg为正离子，Bi为 负离子.在低Q处，S+-负向峰和S,+,S的正向峰的排列形式是由于正负离子耦合以 达到最大的离子屏蔽· 1.400r RMC的另外一个优点是，除可正常计算 Mg-BI-Mg 1200 配位数以处，还可给出近邻配位的分布，因此 会1.000 如何定义配位壳层便显得不是十分重要了·图 言0.800 4为MgBi,液体中的配位分布.由图可知，对 于ne-,在Mg周围有3~4个Bi时配位的可能 慧dipthe 0.400 性最大， 0.200 由于RMC可给出三维立体结构，因此可 0000L 统计计算出配位对的取向或化学键的键角·图 -1.000-0.600 -0.2000.200 0.6001.000 Cos0(°) 5是键Mg-Bi-Mg的角度分布，由图可知 图5熔融Mg2监2g脱2和MgB风.键角的分布 Mg-Bi的角度为60°. 3结束语 尽管RMC方法已经受到了广泛的重视，但它的多维结构信息的唯一性问题还需进一步 检验证实.由统计力学和实验证明，对于偶势(Pairwise potential))Φ(z)占绝对主导作用 的液体，RMC是正确的.但对于多体势占主导作用的液体，在应用RMC时应多加注 意，此时最好将RMC和MMC相比较， 参考文献 I McGreevy R L,Pusztai L.Reverse Monte Carlo Simulation.Mol Sim,1988(1):359-367 2 McGreey R L,Howe M A.Condensed Matter,The Structure of Molten LiCl.J Phs,1989(1):9957~9962 3 Keen D A.McGreevy R L.Structural Modelling of Glass Using RMC.Nature,1990,344:423~424 4 Howe M A,McGreevy R L.Determination of Three Body Correlations in Liquids by RMC.Phys Chem Liq,1991,24:1~12 5 Chuntai Guo.Ph.D.Thesis],University of Bristol,1994 Reverse Monte Carlo Computer Simulation Guo Chuntai ng,19902:2773 Department of Physical Chemistry.USTB,Beijing 100083,PRC ABDTRACT Reverse Monte Carlo Computer(RMC),simulation is introduced in defail,and a brief comment is given.Application of RMC is also introduced by way of an example. KEY WORDS computerized simulation/solution structure,monte carlo metlodVb l . 18 N O . 1 郭春泰 : R e w 晓祀 M o n t e C妞ir o 计算机模 拟 通 过 R M C , 求 出 3 个 偏 R D F , 这 对 无 法 测 定 偏 结 构 因 子 的 体 系 是 十 分 有 益 的 . 图 2 表 明 M g 3 iB Z 具 有 典 型 的 离 子 化 合 物 (例 N a C I , Z n C 1 2 ) 结 构 形 式 . M g 为 正 离 子 , iB 为 负 离 子 . 在 低 Q处 , S + 一 负 向 峰 和 S 、 、 , S 一 达 到 最 大 的 离 子 屏 蔽 . R M C 的另外 一 个 优点 是 , 除 可正 常计 算 配位数 以处 , 还 可给 出近邻 配位 的分布 . 因此 如 何定 义 配位壳层便 显得不 是十分 重要 了 . 图 4 为 M g 尹1 2液 体 中的配位 分布 . 由图可知 , 对 于 n gM 一 。 , 在 M g 周围有 3 一 4个 iB 时配 位 的可 能 性最 大 . 由于 R M C 可给 出三维立 体结 构 , 因 此 可 统计计算 出配位 对 的取 向或化 学键 的键角 . 图 5 是 键M g 一 iB 一 M g 的角度 分 布 . 由 图 可 知 M g 一 iB 的角度 为 60 0 . 的 正 向峰 的 排 列 形 式 是 由 于 正 负离 子藕合以 { g M g 一 B I 一 州 … 尹入 沁协沙州 一 一 { \ 一从 ( ( D )屯/) 团 co · g 一 1 0 00 一 0 . 600 一 0 2 0 0 0 2 0 0 0 6 0 1 . 因 0 oC s (口/( 。 ) ) 图 5 熔融 M忿妙 . 2) 、 M g J残 和 M g 3 . I ` M 键角的分布 3 结束 语 尽管 R M C 方法 已 经受到 了 广泛 的重视 , 但 它的多 维结 构信 息的唯 一性 问题还需 进 一 步 检验 证 实 . 由统计力 学和 实验证 明 , 对 于 偶 势 (aP i“ 喊哭 p o ent 血l) 中 (r . 2 ) 占绝 对 主 导 作 用 的液 体 , R M C 是 正 确 的 . 但 对 于 多 体 势 占 主 导 作 用 的液 体 , 在 应 用 R M C 时 应 多 加 注 意 , 此 时最 好将 R M C 和 M M C 相 比较 . 参 考 文 献 1 M d 3 代笼竹 R L , P 出川旧1 L . R e记眨祀 M o n te C妞r】0 Sir n 创巨山n . M o l S而 , 19 8 ( l ) : 3 59 一 3 6 7 Z M以3 代犯 y R L , H o 忱 M A . 〔b 以北m ed M a t et r , , 11℃ S tr u Ct 切 re o f M o l te n L IC I . J P比 , 1989 ( l ) :卯5 7 一 望拓2 3 K“ 泊 D A, M困 er 理 R L . Str u tC t皿 I M以允 ha 唱 of G al s sU ign R M C . 卜抽 t uer , l 望刹〕 , 少抖 : 4 2 3 一 4 24 4 H o 讹 M .A M 以3 代先 W R L . L兄 t e r 面n a it o n of T b 代兄 佳劝y 山讹h t i o 飞 in L 峋山ds 勿 R M C . P h那 C 比m qjL , l卯l , 24 : l 一 12 5 hC un at i G ou . 【hP . D . T l l ` is 」 , U in vers ity o f B出ot l , l卯4 R e ve sr e M o n t e C a r l o C o m P u t e r S im u l a ti o n G z止0 C h u n t a i ng , l男X2) : 27 7 3 eD p a r t m e n t o f p h y s i e a l C h e m i s t 即 , U s T B , B e ij i n g 10 0 0 5 3 , P R e A B D T R A C T R e v e rs e M o n t e C a r l o C o m P u t e r ( R M C ) , s im u l a t i o n 1 5 i n t r o d u e d i n d e fa il , a n d a b r i e f co m e n t 1 5 g i v e n . A P P lica t i o n o f R M C 1 5 a l s o i n t r o d u e d b y w a y o f a n e x a m P l e . K E Y W O R D S co m P u t e r i z e d s im u l a t i o n / s o l u t i o n s t ur e t u er , m o n t e c a r l o m e t l o d