正则系综模型.pdf

第二章 重要抽样的Monte Carlo模拟 §2.3 正则系综的统计力学模型 ±2π ，相应于+2π 的为正涡旋，−2π 的为反涡旋。相变温度下正涡旋和反涡旋是 同时出现的。对于这个体系，除了计算比热等之外，还可以计算涡度（单位面积 上的涡旋数目）与温度等的关系。 2.3.2.2 Heisenberg 模型 在 Heisenberg 模型中，自旋可取3维空间中的任意方向，即 Hamilton 量为 E =−J σ ⋅σ =−J σ σ +σ σ +σσ ， (2.3.2.2-1) ∑ i j ∑( ix jx iy jy iz jz ) i , j i, j 同样可以采用经典自旋的矢量模型。因此，Ising 模型中自旋矢量只有一个分量， X Y 模型中有两个分量，在 Heisenberg 模型中有三个分量。注意模型是以自旋矢 量分量的数目区别的，每种模型中自旋占据的格子点阵则可以分别有一维、二维 和三维。除了对J 0 的铁磁性和J 0 的反铁磁性加以区别外，在各分量之间还 可进一步加上不同的权重因子以模拟各向异性的情况。除此之外，外加磁场方向 只在空间中设定了某个自旋分量的优先取向，而对其它两个分量则无影响。 在系统最低能量的状态即基态上，所有自旋按照完全有序的平行或反平行排 列，在有限的温度下由于热激发出现能量较高的状态。在 Ising 模型中这样的激 发态是自旋的反转，在 Heisenberg 模型中可以出现周期性的自旋波激发，它是系 统的一种集体行为（图2.3.2.2-1 ），其量子称为磁子（与晶格振动的声子相对应）， 由自旋波的 Bose 统计可以推导出磁化强度与温度成T 3 2 的关系，该关系已经有 实验结果验证。而在一维X Y 模型中可以出现孤立子或孤立波激发，一维链上的 自旋发生2π 的扭转，在反铁磁情况下这种扭转有3种：一种自旋子格子上的自旋 扭转π ，而另一种扭转−π ；两种子格子的自旋各自扭转π ；一种子格子上的自 旋不变，而另一种扭转 2π 。对于三维 Heisenberg 模型，其中的一个自旋分量在 外加磁场下呈现有序，而另外两个分量出现类似于 Kosterlitz-Touless 相变的束缚 拓扑态激发。随着温度升高，涡度增加但正负涡旋束缚态开始解离（图2.3.2.2-2 ）。 图2.3.2.2-1 三种自旋系统的激发态模式：自旋 图2.3.2.2-2 反铁磁自旋系统的KT 相变，在高 波、涡旋、孤立子。 温下涡旋对发生解离。 2 －42 第二章 重要抽样的Monte Carlo模拟 §2.3 正则系综的统计力学模型 图2.3.2.3-1 Potts 模型中格点的广义自旋取值形成的畴区域。 2.3.2.3 q 态Potts 模型 统计力学中的另外一个重要的模型是q 态 Potts 模型，它在介观尺度的微结 构模拟中有着广泛的应用。Potts 模型是 Ising 模型的推广，其主要差别在于广义 的自旋σ和使用不同的 Hamilton 量。σ 的取值可为q 个离散状态，以代替 Ising 模型中只能取“向上”和“向下”的布尔自旋变量，系统的 Hamilton 量是， E =−J ∑δ ， (2.3.2.3-1) σσ i j i, j 其中每个格点上的σ取值为1, 2,…,q 。只有当最近邻格点对有相同的态时，这对 格点间才形成一条连接键。显然，Ising 模型（不考虑无意义的能量最小值）相 当于q 2 的Potts 模型。 引入自旋的离散谱后，可以将有相同广义自旋取值的区域表示为广义畴，在 微结构的模拟中，畴可以理解为晶体材料中具有相似取向的区域，广义自旋取向 则描述了晶格能量、表面能、位错密度等与结构有关的基本物理量，从而将原始 研究自旋的模型推广到可以模拟材料的微结构。如图2.3.2.3-1所示，晶粒可以用 广义自旋有