马尔可夫链理论和Monte Carlo 取样的实现.ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,马尔可夫链理论和,Monte Carlo,取样的实现,Monte Carlo,取样不直接使用,P(X),而是以某种方式取样,大量的样本最终符合所需的分布,.,Monte Carlo,利用转移矩阵,从当前,x,生成一个状态,y.,马尔可夫链和,Monte Carlo,马尔可夫链是一个简单的随机过程,给定转移矩阵,W,可以得到平衡分布,P.,Monte Carlo,是马尔可夫链的计算机实现,.Monte Carlo,中,P,给定,我们需要寻找,W,使得,P=P W.,各态历经,对所有,n,n,max,所有,x,和,x,细致平衡,马尔可夫链的收敛条件,(,充要,),证明,考虑很多个平行的,Markov,链,在一个给定的某一步,有,N,r,个链处于第,r,个态,N,s,个链处于第,s,个态,.,于是在下一步从,r,态到,s,态的数目为,从,s,态到,r,态的数目为,从,r,态到,s,态的净转移的数目为,若,w(,x,r,!,x,s,),满足细致平衡条件,则上式成为,这是一个十分重要的结果,上式表明,如果二个状态之间,不满足分布,P,则这一,Markov,过程的演化结果将总是,使其趋于满足,.,这样,就证明了我们的论断,.,Metropolis,算法,(1953),Metropolis,算法取,T,为一个对称的转移矩阵,The Paper(7500 citations from 1988 to 2003),THE JOURNAL OF CHEMICAL PHYSICS VOLUME 21,NUMBER 6 JUNE,1953,Equation of State Calculations by Fast Computing Machines,NICHOLAS METROPOLIS,ARIANNA W.ROSENBLUTH,MARSHALL N.ROSENBLUTH,AND AUGUSTA H.TELLER,Los Alamos Scientific Laboratory,Los Alamos,New Mexico,AND,EDWARD TELLER,*,Department of Physics,University of Chicago,Chicago,Illinois,(Received March 6,1953),A general method,suitable for fast computing machines,for investigating such properties as equations of state for substances consisting of interacting individual molecules is described.The method consists of a modified Monte Carlo integration over configuration space.Results for the two-dimensional rigid-sphere system have been obtained on the Los Alamos MANIAC and are presented here.These results are compared to the free volume equation of state and to a four-term,virial,coefficient expansion.,1087,The Calculation,Number of particles N=224,Monte Carlo sweep 60,Each sweep took 3 minutes on MANIAC,Each data point took 5 hours,MANIAC the Computer and the Man,Seated is Nick Metropolis,the background is the MANIAC vacuum tube computer,正则分布的抽样方法,:,选择一个满足细致平衡条件的转移几率,;,产生一个,Markov,链,丢掉链的前而面,M,个状态,;,用其余状态进行物理量的计算,.,考虑从,r,态到,s,态的转移,若二状态的能量差为,则,:,当年,Metropolis,选择,:,目前常用的另一种选择是,:,应当注意的是,w,的选择并不唯一,只要满足细致平衡条件的要求即可,但不同的,w,收敛速度往往差别很大,如何选择合适的,w,以达到尽可能快的收敛速度和尽可能高的计算精度仍然是当前,Monte Carlo,算法研究的前沿课题之一,.,The,Ising,Model,-,+,+,+,+,+,+,+,+,+,+,+,+,+,+,+,+,-,-,-,-,-,-,-,-,-,-,-,-,-,-,-,-,-,-,-,The energy of configuration,s,is,E(,s,)=-J,s,i,s,j,where,i,and,j,run over a lattice,denotes nearest neighbors,s,=1,s,=,s,1,s,2,s,i,例题,Ising,模型的模拟,Ising,模型,:,式中,J,称为交换积分,h,为外场,s,i,可取值,(1,-1),称为自旋变量,.,Ising,模型是最简单的非平庸统计物理模型,它是由德国物理学家,Lenz,在二十年代提出的,这一模型可用来描述单轴各向异性磁性系统,合金等物理体系,同时也是一个十分有兴趣的理论模型,.,例题,Ising,模型的模拟,Ising,最早给出了这一模型在一维情况下的严格解,证明了在一维下这一模型不存在相变,.Onsager,于,1944,年做出了零场下这一模型在二维空间的严格解并计算了它的相变温度,比热在相变点的行为等热力学量,.,杨振宁在,1952,年解出了外场很小时二维空间的,Ising,模型,求出了序参量的临界行为,.,由于对这一模型的很多形为目前了解的比较透彻,因此它经常被用来做为检验各种数值方法或解析近似方法的标准,.,1969,年，,A.E.Ferdinand,和,M.E.Fisher,求得了有限尺寸,2D,Ising,模型在周期性边界条件下的严格解，成为检验模拟结果的一个有效标准。,感兴趣的物理量,平均能量,比热，利用下面的公式,:,3.,磁化,=,磁化率,5.,Binder 4,阶累积量,自旋相关函数,时间相关函数,2D,Ising,模型的比热,From D P Landau,Phys Rev B 13(1976)2997.,一个算法,选择一个格点,i,其自旋将考虑作翻转,s,i,!,-,s,i,.,计算与此翻转相联系的能量变化,H.,计算这一翻转的转移几率,w.,产生一在,0,1,之间均匀分布的随机数,.,如果,w,则翻转该自旋,否则,保持不变,.,不论何种情况,其结果都作为一新的状态,.,分析该状态,为计算平均值收集数据,.,有限尺寸标度与相变,涨落,:,h,A,2,i,通常与某种响应函数相联系,:,涨落耗散定理,!,若,A(x,),为系统的哈密顿量,则平均每个粒子的能量,h,E,i,为,:,比热是一个强度量,由此可以推断,:,自平均效应,.,当系统的尺度趋于无限时,其涨落趋于,0!,有限尺寸标度与相变,这实际上是中心极限定理的一个结果,.,考虑把系统分成一些小的系统,小系统之间的相互作用可以忽略时,(,这总可以做到,?!),就是一系列独立,同分布的样本,从而中心极限定理成立,.,由此,能量的分布可以写为,:,P(E),N=1,5,10,20,30,h,E,i,=1,K,B,T,2,c=1,分布宽度,/,1/N,1/2,Ising,模型,:,平均磁化,:,自发磁化,:,注意求极限的顺序,!,有限尺寸标度与相变,由于,Ising,模型的哈密顿量在,h=0,时具有,S,i,$,S,i,的对称性,所以,如果先取,h=0,则得不到自发磁化,.,自发磁化是热力学极限下的产物,由各态历经破坏而得到,.,有限尺寸标度与相变,由于,:,磁化率,:,有限尺寸标度与相变,由此,:,基于在讨论能量时的同样的分析,我们得到序参量,(,磁化,),的分布,.,当,h=0,时,由于对称性,M,和,-M,是对称的,.,在临界温度以上,M=0,所以,当,T,T,c,时,为了反映这种对称性,分布函数写为,:,有限尺寸标度与相变,P(S),N=1,5,10,20,30,h,S,i,=0.8,K,B,T,=1,有限尺寸标度与相变,如何计算序参量,?,注意到前面的分析,对于有限系统,如果直接计算,h,S,i,则在任何温度都得到的是,0!,注意到自平均效应,有下述结果,:,可以使用上述任何一种进行计算,然后外推到,N,!,1,.,有限尺寸标度与相变,磁化随,MCS,的变化，有限系统的各态历经演示。,2D,Ising,模型，,L,64,有限尺寸标度与相变,相变点会发生什么,?,在二级相变点,(,现一般称为临界点,二级相变也称为连续相变,),系统的比热,磁化率等物理量发散,发散的根源在于相关长度发散,.,这样,前面的分析将不再成立,.,在临界点,热力学量具有很有趣的标度行为,:,以单轴各向异性磁体为例,(,理论上近似以,Ising,模型描述,).,在,T,!,T,c,时,引进约化温度,:,有限尺寸标度与相变,比热,:,磁化率,:,自发磁化,:,关联长度,:,当,t=0,时,磁化与外场之间有,:,有限尺寸标度与相变,称为临界指数,是反映临界点本质的量,.,这些指数的来源与关联长度的发散有关,.,还有一个指数,这里不做介绍,.,当系统有限时,不会有发散,!Monte Carlo,只能计算有限体系,如何从,Monte Carlo,计算获得相变和临界现象的信息,?,有限尺寸标度与相变,如何找临界点,?,如何算临界指数指数,?,为此,我们考虑有限体系的标度理论,.,基本出发点,:,设有限系统的线度为,L,则关联长度最多为,L,T,C,时的序参量分布,序参量的有限尺寸效应,有限尺寸标度与相变,磁化率的有限尺寸效应,T,T,C,时的序参量分布,TT,C,时的序参量分布,考虑磁化率,:,在临界点附近有一个峰,最高处的温度定义为有限系统的临界温度,T,c,(L,),峰宽为,T.,当,L,!,1,时,T,c,(L,),!,T,c,(,1,),T,!,0,假定,:,磁化率的极大值,:,磁化在,T,c,(L,),不为,0,而是,:,Fisher,假定,:,在临界点只有一个重要的特征长度,即关联长度,.,有限尺寸标度与相变,由于,磁化率峰的宽度对应于,于是,有限尺寸标度与相变,由此得到,:,Finite-Size Scaling,Shift of,T,c,From A M,Ferrenberg,and D P Landau,Phys Rev B,44,(1991)5081,Accurate Exponent Ratio,From J S Wang,R H,Swendsen,and R,Koteck,Phys Rev.B,42,(1990)2465.,有限尺寸标度与相变,Binder,累积量,:,性质,1:T,远离临界区域,有限尺寸标度与相变,性质,2:T,T,c,L,远离临界区域,有限尺寸标度与相变,性质,3:T,T,c,L,在临界区域,U,*,是一个与,L,关系很小的普适值,.,应用,:,确定,T,c,对不同的,L,计算,U,L,对各种对,(L,L),计算比值,U,L,/U,L,并对,T,作图,这些比值的曲线将交于一点,T,c,由此方法可以非常精确地定出临界温度,.,对于一级相变,U,L,也很有用,.,证明,:,在临界区域,由于,L,所以前面给出的分布失效,.,这里的关键假定是,做为,L,S,(,对应于温度,),函数的,P(S),实际上只是两个变量的函数,.,有限尺寸标度与相变,计算二阶和四阶距,:,有限尺寸标度与相变,类似的有,:,于是,:,当,T,!,T,c,时,!,1,与,L,无关,有限尺寸标度与相变,Binder,累积量的尺寸效应,确定,T,C,的方法,有限尺寸标度与相变,不同方法计算自发磁化，外推到无限大系统时，趋于同一个值。,