[论文解读] Hysteresis, Avalanches, and Noise: Numerical Methods

Matthew C. Kuntz, Olga Perković+3 more|arXiv (Cornell University)|Sep 7, 1998

Theoretical and Computational Physics被引用 45

一句话总结

本文提出两种高效的数值算法——有序列表法与基于位的存储方法，用于模拟大规模零温随机场伊辛模型，实现了高达十亿自旋的模拟。这些方法的时间复杂度为 O(N log N)，内存使用最高可减少 96 倍，从而能够准确提取临界指数，并观察到临界无序度 R_c ≈ 2.16 附近的分形雪崩动力学。

ABSTRACT

In studying the avalanches and noise in a model of hysteresis loops we have developed two relatively straightforward algorithms which have allowed us to study large systems efficiently. Our model is the random-field Ising model at zero temperature, with deterministic albeit random dynamics. The first algorithm, implemented using sorted lists, scales in computer time as O(N log N), and asymptotically uses N (sizeof(double)+ sizeof(int)) bits of memory. The second algorithm, which never generates the random fields, scales in time as O(N \log N) and asymptotically needs storage of only one bit per spin, about 96 times less memory than the first algorithm. We present results for system sizes of up to a billion spins, which can be run on a workstation with 128MB of RAM in a few hours. We also show that important physical questions were resolved only with the largest of these simulations.

研究动机与目标

• 实现以往极限之外的大规模无序系统滞后与雪崩模拟。
• 为零温随机场伊辛模型开发内存与时间高效的数值算法。
• 通过模拟高达十亿自旋的系统，准确提取普适临界指数。
• 观察并表征临界无序度 R_c 附近的分形雪崩结构与幂律分布。
• 为研究人员提供开源代码与工具，以复现和扩展模拟。

提出的方法

• 有序列表算法使用有序数据结构追踪自旋翻转，时间复杂度为 O(N log N)，内存使用与 N × (sizeof(double) + sizeof(int)) 成正比。
• 基于位的算法仅用一位存储每个自旋的翻转状态，与有序列表方法相比，内存使用减少约 96 倍。
• 模型采用 D 维超立方晶格，周期性边界条件，当局部场符号变化时执行确定性自旋翻转。
• 当自旋的有效场 h_i^eff = J∑s_j + h_i + H 发生符号变化（由于邻近自旋翻转或外场 H(t) 增加）时，自旋发生翻转。
• 模拟采用 H(t) 从 -∞ 到 ∞ 的缓慢上升过程，以捕捉所有雪崩与滞后回线的动力学。
• 代码支持动态可视化、磁化率追踪、雪崩大小分布与关联函数计算，利用 C++ 标准模板库实现。

实验结果

研究问题

• RQ1 在三维中，系统从有限雪崩过渡到无限雪崩的临界无序度 R_c 是多少？
• RQ2 在临界点附近，雪崩大小与时间分布如何标度，其关联临界指数为何？
• RQ3 基于位的算法是否可在 128MB 工作站上数小时内模拟十亿自旋的系统？
• RQ4 大雪崩的分形维数与三维空间相比如何，这对其结构有何含义？
• RQ5 当 R 不在 R_c 时，关联函数 G(x,R) 与雪崩大小分布 D(S,R) 如何偏离幂律？

主要发现

• 基于位的算法可在 128MB 工作站上数小时内实现十亿自旋的模拟，时间复杂度为 O(N log N)，且内存使用极低。
• 在 R = 2.3（距 R_c ≈ 2.16 仅 6% 以内）时，大雪崩表现出复杂、分形的分支结构，内部存在空洞，且随时间逐渐填充，表明其分形维数略小于三维。
• 雪崩大小分布 D(S,R) 在 R_c 附近呈现幂律分布，指数 τ̃ ≈ 1.5，与实验观测的 Barkhausen 噪声一致。
• 关联函数 G(x,R) 在小距离处呈现幂律衰减，但在 R 不在 R_c 时于大距离处出现偏离，表明非临界区域行为非普适。
• 滞后回线在 R_c ≈ 2.16 处出现宏观跳跃，标志着无限雪崩的出现与临界标度行为的开始。
• 有序列表与基于位的算法均实现 O(N log N) 复杂度，但基于位的方法将内存使用减少约 96 倍，使大规模模拟成为可能。