计算材料学Ising模型实验报告.docVIP

Monte Carlo实验报告 项目名称：Ising 模型 二、项目内容概要 1、编译和运行 进入实验的文件夹：cd□~/sourcecode/2D_Ising 文件夹里有源代码mc2d.f和输入文件in.2d 阅读理解并编辑输入文件：gedit□in.2d 之后编译mc2d.f f95 mc2d.f -o mc2d.exe 运行可执行文件 ./mc2d.exe 查看刚刚生成的四个输出文件，四个文件的内容如下： file1.out：温度；时间；单位原子能量；单位原子磁化强度 file2.out：温度；单位原子能量；能量变化；单位原子磁化强度；磁化强度变化；单位原子热容 file3.out：温度；自旋构型 file 4.out：温度；能量升高而被接受的数目；能量下降而被接受的数目；被拒绝的数目 2、gnuplot 作图 作温度与能量图：p “file2.out” u 1:2 w p ps 3 pt 5 作出file2.out 中第1 列与第2 列数据； 作温度与磁化强度图：p “file2.out” u 1:4 w p ps 3 pt 5 作出file2.out 中第1 列与第4 列数据 作温度与热容图：p “file2.out” u 1:6 w p ps 3 pt 5 作出file2.out 中第1 列与第6 列数据 三、项目实施方法/原理 1925 年，伊辛提出描写铁磁体的简化模型：设有N 个自旋组成的d 维晶格 (d=1,2,3)，第i 格点自旋为Si=±1(i=1,2,…N; ±代表上下)。只考虑最近邻作用， 相互作用能为±J(J0 为铁磁性, J0 为反铁磁性)，平行为-J，反平行为J。 伊辛模型的蒙特卡洛模拟基本步骤如下： 项目实施结果： 各种情况下能量温度曲线 铁磁正方形点阵温度和能量曲线 铁磁三角形点阵能量与温度曲线 反铁磁性正方形点阵能量温度曲线 反铁磁性正方形点阵外场为1时能量温度曲线 反铁磁性正方形点阵外场为0.5时能量温度曲线 各种情况下磁化强度和温度的关系曲线 铁磁正方形点阵磁化强度能量曲线 铁磁三角形点阵磁化强度温度曲线 反铁磁性正方形点阵磁化强度温度曲线 反铁磁性正方形点阵磁化强度温度曲线（外场为0.5） 反铁磁性正方形点阵磁化强度温度曲线（外场为1） 铁磁性正方形点阵磁化强度温度曲线（外场为0.5） 铁磁性正方形点阵磁化强度温度曲线（外场为0.5） 各种情况下热容和温度的关系图 铁磁正方形点阵热容能量曲线 铁磁三角形点阵热容能量曲线 反铁磁正方形点阵热容能量曲线 反铁磁正方形点阵热容能量曲线（外场为1） 反铁磁正方形点阵热容能量曲线（外场为0.5） 铁磁正方形点阵热容能量曲线（外场为0.5） 铁磁正方形点阵热容能量曲线（外场为1） 五、项目小结： 1.在保持原参数不变的情况下，可以得出，温度越高，原子热运动越剧烈，因此单个原子的能量也就越高。且在3K左右的温度范围变化率为最大。 通过比较可以发现，在没有加入磁场的情况下，能量随着温度的改变发生的变化与材料是否是铁磁性或反铁磁的没有太大的关系。但是和温度是有很大的关系的。因为反铁磁材料的磁矩为0，不会和外加磁场相互作用，故其几乎没有变化。 当材料为三角形点阵的时候其能量的上升要比正方形点整的要缓慢一些，但是变化规律类似，都会有先缓慢上升在快速上升最后又变为缓慢上升的过程。猜想是因为三角形结构的晶体最近邻有六个原子与中心原子相互作用，而正方形结构只有四个，相互作用力增大因此升温对其影响相对较小。 加入磁场以后，通过对比可以发现铁磁性材料中，加入磁场之后能量的变化与不加磁场的变化不同，这是因为铁磁性材料原子温度升高热震动加剧，但是由于有磁矩，外加磁场会和其相互作用，减弱了加热使震动加剧的效果。故其上升的速度相对没有外加磁场要慢。 2.对于铁磁性材料来说，不同点阵结构会使得材料的磁化性能不同，三角形点阵结构磁化强度为零的温度要比正