热力学与统计物理课程设计.doc

热力学与统计物理课程设计 计算机模拟 (Matlab实现) 题目: Ising 模型铁磁-顺磁相变现象的Monte Carlo 模拟 学院: 数理学院 专业: 06 应用物理 指导老师: 刘旭东 同组人员: 姓 名 学 号 彭加福 0640502112 黎 婕 0640502102 胡进军 0640502107 Ising模型铁磁—顺磁相变现象的Monte Carlo模拟 （江苏科技大大学 06应用物理） 摘要：课程设计中首先采用Ising model的思想建立一个二维的模型，然后利用重要性抽样和Monte Carlo方法及其思想模拟铁磁-顺磁相变过程。计算了顺磁物质的能量平均值Ev、热容Cv、磁化强度M及磁化率X的值，进而研究Ev、Cv、M、X与温度T的变化关系并绘制成Ev-T图、Cv-T图、M-T图、X-T图，得出顺磁物质的内能随着温度的升高先增大而后趋于稳定值；热容Cv、磁化率X随着温度的升高先增大后减小；磁化强度M在转变温度Tc处迅速减小为零，找出铁磁相变的转变温度Tc大约为2.35。 关键字：蒙特卡罗方法、伊辛模型、重要性抽样、顺磁相变、居里点磁畴1907年法国科学家外斯系统地提出了铁磁性假说其主要内容有铁磁物质内部存在很强的分子场在分子场的作用下原子磁矩趋于同向平行排列即自发磁化至饱和称为自发磁化铁磁体自发磁化分成若干个小区域这种自发磁化至饱和的小区域称为磁畴由于各个区域磁畴的磁化方向各不相同其磁性彼此相互抵消所以大块铁磁体对外不显示磁性当温度升高时原子间距加大降低了交换作用同时热运动不断破坏原子磁矩的规则取向故自发磁化强度下降直到温度高于居里点以致完全破坏了原子磁矩的规则取向自发磁矩就不存在了材料表现为强顺磁性其磁化率与温度的关系服从居里——外斯定律， 从上式可以看出：如果J0时，则当所有自旋具有相同取向时系统具有最低的能量，相应于绝对零度下的状态。在足够低的温度下，也会有很多的自旋具有相同的取向，这就是无外场时铁磁具有自发磁化的原因。 ii. 物理量的平均值 , iii.物体的磁化率Χ , ⅳ.物体的磁化强度 , ⅴ.物体的热容 , 推导过程如下 , 又: , 所以有 , 即是 , 四、建立模型 下面用Matlab语言编写了用Monte Carlo Method显示在二维伊辛模型在低温条件下发生相变的简化程序。 1、假定有400个磁矩，呈20×20晶格排列，初始状态有200个磁矩朝上，200个磁矩朝下，并在晶格中交替排列，如下图所示（其中Ising（i,j）=1表示处于第i行第j列的晶格点上磁矩朝上并用黑色三角形表示，Ising（i,j）=-1代表该磁矩朝下，并用红色三角形表示）。 2、用rand语句产生两个随机数，取整并令之等于i和j。改变Ising（i,j）的符号计算内能U，并用以下方法判断是否接受： 若，说明磁畴磁化方向，这种情形下，也是接受的。 若，说明磁畴磁化方向(1,m)= Ising (21,m) Ising (22,m)= Ising (2,m) Ising (n,1)= Ising (n,21) Ising (n,22)= Ising (n,2) 4、本实验中取相对单位,取K=1，J=1。 5、程序流程 流程说明： 首先建立形如1中提到的模型，用22×22的矩阵表示（实为20×20的矩阵，四周的元素作为边界条件），正负号代表磁化方向，用两个for-end语句产生并绘制成图。 其次产生两个随机数sai和saj作为矩阵的行标和列标如2所述，将每一组（sai,saj）对应元素翻转（改变正负号），再算内能U并用2中的方法判断是否接受，如此运行50万次。接着计算每一次的能量E及磁化强度Mo，运行50万次，统计得到每个温度对应的Ev，Cv，M，X，绘制成图。 五、结果分析 1．平均能量与温度的关系图 从平均能量与温度的