三维Ising模型.PDF

三维 Ising 模型 祁永晖 2001012521 清华大学物理系 一．关于选题的说明 1 ．我选的题目是：统计力学三维 isingmodel 。这个题目是当作期末考试论文交的。 2 ．另外，我又做了非线性方程-- kdv 孤子的传播和非线性方程—Sine-Gordan 孤子的振 动。 二．关于程序的说明 （一）三维isingmodel 1． 算法：Metropplis 算法 简述：假想有一个随机行走者在三维空间运动，相继两步的终点产生出点子的一个序 列。行走越长，它连接的点子就越接近真实的分布。 规则：从 X 走出试探步 X , 这个新点可以用任意方法产生，例如可以在 X 点周围的一 n t n w(x ) 个边长为δ很小的多维体中均匀地随机的选取。然后按照比值，r i 决定是接受还是 w(x ) t 拒绝这一步，如果 r 1, 那么接受这一步，即取 X = X ,如果 r1,则以概率 r 接受这一 n+1 t 步，这时要再产生一个[0,1]之间的随机数η，若ηr ,就接受这一步，否则不接受，并取 X = X ,然后从新的 X 出发，产生 X 。任意一点 X 都可以作为随机行走的起始点。 n+1 n n+1 n+2 0 2 ．原程序：具体说明看M-file isingmodel3D.m 重点说明： （1）c=sign(0.5-rand(ny,nx,nz)); %建立三维自旋矩阵 S （2 ）只考虑同前一个位形仅相差一个自旋反转的试探位形 S － S ，它们只有一个自旋 t 试探点值 ω(S ) −H (S )+H (S ) t t S S 反转了方向。已知权函数 r ＝ ＝e ，求得 α i , j ,k 出发点值 ω(S ) −2Sα ( Jf +B ) r e 。 （3 ）R(1:7,2)=(exp(-2*(JJ*(2*(1:7)-8)+B))) 例如：B=1T,JJ=0.1 时，R(1:7,2)= 0 0.4493 0 0.3012 0 0.2019 0 0.1353