量子蒙特卡洛的随机级数展开法量子蒙特卡洛的随机级数展开法.pptVIP

基于循环更新的随机级数展开法 三种典型的自旋系统模型 经典情形 Ising模型： XY模型： Heisenberg模型： 量子情形 Ising模型： XY模型： Heisenberg模型： 经典蒙特卡洛模拟回顾 马尔科夫过程 (1) 各态历经：有限时间内遍历尽量多的状态。 (2) 细致平衡原理：为了使系统达到平衡态，需要在状态更新时满足 可令转移概率为 整个蒙特卡洛模拟可以简单地概括为：(1) 初始化构型 (2) 预热：通过以上马尔科夫过程使系统达到平衡态 (3) 达到平衡态后的抽样 随机级数展开(SSE) 首先寻找热力学函数的形式 同理 热容 均匀磁化率 磁化强度 下面考虑一般情况下1/2自旋的Heisenberg模型，Hamilton量写成 其中 对角 非对角 对于某个n，可插入L-n个单位算符以保证算符序列总长度对于不同的n相等。并且当n2为偶数时得到 其中两个邻近传播态最多相差一对相反的自旋 至此完成了模拟的数值化。相应的限制条件为 1. 反铁磁系统 2. 对于一给定算符序列非对角算符数为偶 3. 初态经过n个算符作用后得到的末态自旋构型必须和初态相同(周期性边条件) 4. 模拟进程中需要调整L保证其大于非单位算符数n 对角算符更新 对角算符更新需遍历L个算符 算符循环更新 根据n2为偶以及周期性边条件的要求，引入算符循环更新生成非对角算符 (a) “bounce” (b) “continue-straight” (c) “switch-and-reverse” (d) “switch-and-continue” 当进入点确定时，从某一点退出的概率与从该点退出后生成的vertex权重成正比 两种循环闭合路径 不连续情形 vertex连接 算符循环更新的一般流程 1. 对于当前的算符序列构建vertex连接 2. 进行nl次算符循环更新，nl的选取要使得nl次循环访问vertex的总数约为2n 3. 对于一次循环，首先在4n个点中随机选取一个作为进入点，根据相应概率选择退出点并更新该vertex类型，退出点所连接的点作为下一个vertex的进入点，直到到达最初的点 4. 完成nl次循环后，根据新生成的vertex序列更新算符序列以及有被相应位置算符作用的自旋初态。对于未被算符作用的位置的自旋，有1/2的概率翻转，它们是“自由的”自旋 另外，在一次循环更新中需要循环长度，若经过很多vertex(通常循环长度大于100n)仍不能到达初始点，则强制退出循环并且不更新该次产生的vertex序列。 特殊情形 “bounce”这一过程的权重是自身矩阵元，在任何情况下都应该存在，但是在一些特殊情形可以将其排除。 各向同性反铁磁 各向同性铁磁 XY模型，所有矩阵元均为1/2 排除“bounce”过程可以去除循环中的回溯情形，使循环尽量做到遍历。一般情况下“bounce”过程并不能排除并且该循环也应被计入总的循环次数。 总的模拟流程 do i=1,isteps call diagonalupdate() !对角更新 call loopupdate() !循环更新 call adjustcutoff(i) !如果n+n/3L,L=n+n/3 enddo !到达平衡态 do j=1,nbins do i=1,msteps call diagonalupdate() call loopupdate() call measureobservables() !给出热力学量 enddo call writeresults(msteps,j) !每隔一定步数输出一次结果 enddo