计算材料学Ising模型实验报告Word文件下载.docx

1、 file1.out：温度；时间；单位原子能量；单位原子磁化强度 file2.out：能量变化；单位原子磁化强度；磁化强度变化；单位原子热容 file3.out：自旋构型 file4.out：能量升高而被接受的数目；能量下降而被接受的数目；被拒绝的数目2、gnuplot 作图 作温度与能量图：p “file2.out” u 1:2 w p ps 3 pt 5 作出file2.out 中第1 列与第2 列数据； 作温度与磁化强度图：4 w p ps 3 pt 5 作出file2.out 中第1 列与第4 列数据作温度与热容图：6 w p ps 3 pt 5 作出file2.out 中第1 列与第

2、6 列数据三、项目实施方法/原理1925 年，伊辛提出描写铁磁体的简化模型：设有N 个自旋组成的d 维晶格（d=1,2,3），第i 格点自旋为Si=1（i=1,2,N; 代表上下）。只考虑最近邻作用，相互作用能为J（J0 为铁磁性, Jexp（-dE/T） /判断是否接受翻转 end EE（1,q）=compute（s,j）/（n-2）*（n-2）; /记录能量 mag（1,q）=sum（sum（s）2）/（n-2）*（n-2）; /记录磁化强度endmagEEfunction y=bou（s1）/设定边界的条件函数n=length（s1）;s2=s1;for i=1:n s2（1,i）=s2（n-1,i）;for m=1: s2（i,1）=s2（i,n-1）;y=s2;function y=compute（s,J） /计算能量的函数n=length（s）;E=0;for i=2:n-1 for j=2: E=E+J*（s（i,j）*（s（i+1,j）+s（i,j+1）+s（i-1,j）+s（i,j-1）; end end y=E;6、参考书目 张帆，周伟敏.材料性能学M.上海交通大学出版社 张祥，陈东保，陈武鸣.二维伊辛模型蒙特卡罗模拟J.南京大学学报1997,33（1）