二维微正则Ising模型临界相变的Q2R方法-二维微正则Isin.PDF

2006年11月 重庆大学学报(自然科学版) NOV．2006 第29卷第1 1期 Journal of Chongqing Uniw~rsity(Natural Science Edition) Vo1．29 No．11 文章编号：1000—582X(2006)11—0115—04 二维微正则Ising模型临界相变的Q2 R方法 李立奇，王宗笠 (重庆大学数理学院，重庆 400030) 摘 要：利用元胞自动机Q2R／WQ2R方法研究了二维微正则Ising模型的临界相变，计算结果表明 此数值仿真能够显示二维晶格在孤立条件下的自发磁化，并且观察到了相变临界现象，得到了临界能 量．文中确立的方法可以为三维Ising模型~~4h-真计算打下基础． 关键词：Ising模型；元胞自动机模拟；Q2R／WQ2R方法；临界能量；平均磁化强度 中图分类号：0469 文献标识码：A Ising模型是相变统计理论的一个核心模型．它 旋相互作用能为F=一JZS s ．其中．，为自旋相互作用 是利用格点自旋来描述相变性质．Ising模型确实是一 强度系数，可令J=1．求和只对每个格点的最近邻点进 种最简单却很有效的理论模型，对一维和二维的Ising 行．Ising模型中自旋s 的能量来自它与Von Neumann 模型人们已得出了严格的解析解-l ]．三维Ising模 邻居¨ (位于其上下左右)自旋．s 的耦合．自旋对(S ， 型的严格求解是困难的．一些文献中介绍了求解该模 ．s )的能量计算是：当3个自旋方向相同时为一-，，当两 型的Bragg—Williams近似 和Bethe近似 (考虑 个自旋方向不同时为 ．考虑Ising自旋动力学模型的微 了短程有序)，但推导过程比较复杂． 正则元胞自动机模拟，也就是说，正在考察的是一个与 以上结果是前人在正则条件下，即保持温度不变 环境无能量交换的系统．Q2R规则通过在离散的时间 的条件下利用元胞自动机方法获得的．事实上，由于 演化过程中，网格中任何一个格点自旋的翻转都要保证 相变过程极为迅速，系统与外界来不及交换能量，故笔 局部能量不变来保证全局的能量守恒．Q2R规定：当且 者利用元胞自动机Q2R／WQ2R方法将系统视为封闭 仅当这个运动不引起任何能量交换时，即如果自旋．s 的能量守恒系统，微正则系统，在此条件下，Ising模型 向上的邻居数与自旋向下的邻居数相同，则．s 翻转成 应该存在临界相变，存在临界能量．为了避免配分函 一 s ．但是由于所有的自旋运动都是同时发生的，判断 数的计算，运用数值模拟方法，则能直观地显示Ising 是否翻转是基于邻居不改变的前提，如果邻居也翻转， 自旋点阵指向分布随着能量的下降而由无序向有序的 则无法保证能量守恒．为解决这个问题，可把状态更新 演化过程．Q2R／WQ2R方法正是通过将二维正方网格 分为2步：第一步，根据偶数自旋构形翻转奇数位的自 划分为两套格子来实现状态更新，吻合了实际情况中 旋；第二步，根据奇数自旋构形翻转偶数位的自旋．对 所有的自旋运动的同时性+通过研究二维微正~sing 于二维正方网格，就将原网格划分为两套格子，与国际 模型临界相变，可以找到系统平均磁化随能量的变化， 象棋棋盘黑白格子的分布类似． 尤其是在临界能量附近的行为． WQ2R规则：数值模拟过程显示出：Q2R规则并 不满足egordy条件，即丧失各态历经性，系统在演化 1 Q2R规则和WQ2R规则 过程中不能遍历初始能量 。 给定的相空间区域．这 Q2R规则：Q2R规则是Vishniac 于上世纪8O 一 点可由以下例