目的是为执行 Ising 模型仿真提供一个非常基本的起点。 欢迎任何建议。
2021-11-18 18:59:35 44KB matlab
1
Ising 模型用于通过应用 Metropolis 算法-蒙特卡洛方法模拟磁系统(正、负或随机自旋)。 运行主文件,输入晶格大小(最好是 100),然后选择一个输入自旋作为初始配置。 设置了两个不同的温度(T=2.0 和 T=2.5)。 例如,在 T=2(低温)下使用选定的正自旋进行初始化,大多数自旋是黑色的,这是因为翻转自旋的机会很小,并且材料具有铁磁性。 对于 T=2.5 的情况,由于温度较高,它会产生自旋翻转的趋势。 因此,材料失去磁化。 自旋似乎没有对齐,由此产生的配置似乎是随机无序的。 那是因为顺磁行为。 模拟的下一部分是可观测值计算:平均磁化、平均能量、平均磁化率和比热。 为了计算平均能量和磁化强度,我们必须找到能量和磁化值变化很小的时间依赖性(能量和磁化强度随时间增加变化很小的时间)因此,我们选择精度 p 并检查间隔(应满足精度的时间步数)。 这些间隔应取决于初始配
2021-11-18 18:52:56 3KB matlab
1
在无外场的情况下,自由边界的 方格上,存在着某种粒子,它具备三种自旋状态-1,0,1。我们认为t_initial=200后系统达到稳定态。系统演化至t_max=1200步。 (1)分别绘制单粒子平均磁矩 、单粒子平均能量 与温度 之间的关系图。观察并讨论图中现象。 (2)根据(1)中的现象,合理选取不同的T,分别绘制系统达到稳定态后的演化过程中,系统总磁矩的分布曲线
2021-11-18 18:33:32 6KB 蒙特卡洛 晶格 磁矩
1
二维伊辛模型的高效实现。 网格分为两个子网格(棋盘格)。 在这两个子网格上,可以分别考虑自旋。 这样就可以将计算完全矢量化。 使用简单的笔记本电脑 CPU,每秒可实现高达 3000 万次的旋转翻转。 将磁化强度的温度依赖性与分析结果进行比较。
2021-11-07 18:17:56 3KB matlab
1
IsingModel-蒙特卡洛 日语中的Ising模型的蒙特卡洛(MC)仿真简介。 Ising2d_MC 用日语记录Ising模型和MC模拟(Metropolis方法和热浴方法)。 使用Python的各向同性Ising mocdel的MC模拟代码(计算磁化强度,比热,Binder累积量和相关函数)。 Python代码中的注释为英文。 Ising2d_各向异性 使用Python(计算磁化强度)的各向异性Ising模型的MC模拟代码。
2021-10-11 20:01:55 639KB JupyterNotebook
1
要安装依赖项,请运行:pip install -r requirements.txt 要绘制一个点阵实例化一个绘图对象,参数包括: N- size of lattice B- strength of magnetic field (default is zero) start- low(cold) or high(hot) start (default is low) inc- size of increments in plots (default is 0.01) x0- starting point of plots (default is 1) x1- final point of plots (default is 5) steps-number of steps (default is 50000) T - Temperature (default is 1) 一些展示程
2021-10-11 19:31:16 2.53MB Python
1
使用Fortran 90语言编写,使用蒙特卡洛模拟方法,对线度N=5,10,20,50,100的二维Ising自旋动力学模型进行模拟。该过程使用了周期性边界条件、Metropolis 准则和马尔科夫链进行二维矩阵的演变,并计算了系统演变至平衡态的能量和绝对磁化强度。压缩包中包含完整的Fortran源码和最终的数据曲线图。源码注释齐全易懂,图例清晰明白,欢迎下载交流。
二维伊辛模型phython模拟
2021-05-10 12:02:27 4KB phython Isingmodel montecarlo
1
使用matlab对二维伊辛模型进行数值模拟。通过Metropolis准则判定磁矩是否反转。以图像的方式给出磁矩的分布随时间的变化。可以明显看到铁磁性与顺磁性的差别。
2019-12-21 21:44:37 2KB 统计物理学 伊辛模型 Ising
1