二维 Ising 模型临界指数的有限尺寸标度：跨三种蒙特卡洛算法的收敛性与临界慢化对照

来源说明：本则出自独立撰写的第二批方案。它与前 20 则同样遵循八块结构与硬门槛要求，但撰写时未做英文文献检索，新颖性边界依据的是历届获奖图谱与既有知识，而非当轮查新。因此其「需核实」条目更多，启动前须自行补一轮英文检索。

1 · 研究问题

Metropolis、Wolff 团簇与热浴三种蒙特卡洛算法，在 L = 16–128 的二维 Ising 模型上给出的临界温度与指数比 β/ν、γ/ν 是否在各自误差棒内一致并符合 Onsager 精确解？三者的自关联时间随尺寸的动力学标度指数 z 相差多少、对给定精度目标各需多少墙钟时间？

2 · 研究背景与空白

二维 Ising 模型是唯一同时具备"非平庸相变 + 精确解"的多体系统:Onsager(1944)给出 T_c = 2/ln(1+√2) ≈ 2.269(以 J/k_B 为单位)与 β = 1/8，是检验一切蒙特卡洛管线的黄金标准。局部更新算法(Metropolis、热浴)在临界点遭遇临界慢化(z ≈ 2.2),Wolff(1989)团簇算法把 z 压到接近 0——这是计算物理教学的经典叙事。

各算法的指数与 z 值在文献与教材(Newman & Barkema《Monte Carlo Methods in Statistical Physics》)中均有发表；但教学实现常常只跑单一算法、不给误差棒。

空白在于可复现性维度：同一代码框架、同一随机数流水线、同一误差估计口径下三算法的"精度-成本"完整对照表，配开源代码，公开资料中意外地少。物理结论全部已知，正因如此代码对错可被无情检验——是训练科学计算规范的理想课题，且"对照表"本身即是被计算领域承认的贡献形式。

3 · 可检验假设

• H1：三算法由 Binder 累积量交叉法得到的 T_c 互相一致且与精确值偏差 ≤ 0.5%；有限尺寸标度得到的 β/ν 与精确值 1/8÷1 = 0.125 偏差 ≤ 5%。
• H2:Metropolis 的动力学指数 z 落在 2.0–2.4,Wolff 的有效 z ≤ 0.5；等目标误差下 L = 128 临界点处 Wolff 的墙钟时间优势 ≥ 30 倍。

4 · 量化验收标准

1. 方法学校验(硬门槛)：与解析解极限对照 + 分辨率(尺寸)收敛扫描——L = 4 小格子的能量/磁化与穷举精确求和(2^16 态，直接计算)逐点一致到浮点精度;T_c 与 Onsager 值偏差 ≤ 0.5%；比热峰位随 L 的移动符合 1/L 修正趋势。此条不过关，后续全部结论无效。
2. 统计口径预先写死：一切观测量误差用分块 jackknife(块长 ≥ 10 倍积分自关联时间)；自关联时间用积分法，窗口截断判据(自洽窗口 c = 6)写死；指数拟合的置信区间用对 bin 的自助法 ≥1000 次。
3. 尺寸序列 L = 16/24/32/48/64/96/128，每尺寸临界区 ≥ 10 个温度点；热化丢弃 ≥ 20 倍自关联时间并报告。
4. 随机数:PCG64(numpy 默认)，种子全部记录；同种子重跑结果逐位一致。
5. 三算法共用观测量与分析代码，只有更新核不同(保证对照公平)，代码结构在报告中说明。
6. 可复现性：代码 + 运行配置 + 原始时间序列开源，单尺寸算例可在 1 小时内重跑验证。

5 · 数据与工具

用途	来源 / 工具
实现	Python + numba(JIT；开源)自写三算法；参考实现见 Newman & Barkema 教材伪代码
运行量估算	numba 加持的 Metropolis 单自旋更新 ~10⁷–10⁸ 次/秒(常见笔记本 CPU 公开基准量级);L=128、10⁶ 全格扫描 ≈ 1.6×10¹⁰ 次更新 ≈ 数十分钟/温度点；全项目累计 ≤ 30 CPU 小时，可夜间后台完成
精确解对照	Onsager 1944 解析结果(教材级公式),仅用于校验与对比，不计入本项目的数据贡献
分析	numpy/scipy + matplotlib;jackknife/自助法自写

6 · 方法路径

1. 写 Metropolis 单核 + 观测量框架，先过 L = 4 穷举对照(验收 1 前半)。
2. 实现 Wolff 与热浴，三算法过同一穷举对照。
3. 小尺寸(L ≤ 32)扫温,Binder 交叉定 T_c，对照 Onsager(验收 1 后半)。
4. 全尺寸生产运行(夜间后台)，原始序列落盘。
5. 自关联时间与 z 的测量(H2)，等误差成本核算(墙钟实测)。
6. 有限尺寸标度提取 β/ν、γ/ν + 自助法置信区间(H1)。
7. 独立交叉校验：用重加权(single-histogram reweighting)在单温度点数据上复算峰位，与直接扫温一致 ≤ 1σ。

7 · 新颖性边界

• 本课题不声称任何新物理:T_c、临界指数、临界慢化与团簇算法优势全部为已发表结论(Onsager 1944;Wolff 1989；教材)。物理结论已发表，不得声称为本项目发现。
• 已有工作：算法原始论文 + 教材级实现；分散的 z 测量。
• 本项目贡献(主结论)：同框架、同误差口径下三算法的收敛性与精度-成本对照表 + 一条可操作的"给定精度选算法与格子大小"判据 + 开源全套。属"方法层面的可复现性贡献"——在计算领域被承认，但定位必须讲清楚，否则会被误读为课后作业。答辩时的叙事重点是误差分析与成本核算的严谨性。
• 核对获奖清单:"'Phase Transition' in a Mechanical System: Rotation-Induced Spontaneous Symmetry Breaking"(2025 金奖)是宏观力学系统的对称破缺实验，与统计力学蒙特卡洛无方法重叠，但说明"相变"叙事受评委欢迎——报告可借该题说明相变普适性话题的听众基础，并点名区分。

8 · 决策门槛(go / no-go)

⚠️ 第 4 周末:L = 4 穷举对照通过(逐位级一致)。不通过禁止前进——此关卡专为暴露更新核 bug 设计。 ⚠️ 第 8 周末：小尺寸 T_c 达标(≤ 0.5%)。若 numba 性能不足以支撑 L = 128(单温度点 > 2 小时),降级：尺寸上限降至 L = 64(文献与教材表明该范围足以提取指数，误差棒略宽)，全部分析框架不变；同时把"Python/numba 的性能边界实测"写进成本对照一节，变损失为内容。 ⚠️ 第 20 周检查点：三算法一致性若出现 > 2σ 分歧，最可能是自关联低估——加倍运行长度复核；分歧被修复的过程写入"误差分析"一节。 ⚠️ 预算裁剪顺序：热浴算法可砍(保 Metropolis vs Wolff 二元对照)；温度点可从 10 减到 7。

---