双谐波标准映射中金比例不变环面破缺阈值曲线：Greene 剩余判据的高精度测定

1 · 研究问题

在双谐波推广标准映射 f(x) = (K₁ sin 2πx + K₂ sin 4πx)/2π 中，金比例旋转数不变环面（golden-mean invariant torus）的破缺阈值在 (K₁, K₂) 参数平面上构成怎样的临界曲线？进一步：沿该曲线，Greene 剩余判据（residue criterion）的临界剩余值 R_c 与破缺标度是否保持与单谐波情形相同的普适值？

2 · 研究背景与空白

KAM 环面破缺是哈密顿混沌的核心问题。Greene 1979 用周期轨道剩余判据测得标准映射金环面破缺阈值 K_c = 0.971635…，MacKay 的重整化理论解释了临界剩余值 R_c ≈ 0.25 与自相似标度的普适性；判据的合理性有后续严格化（MacKay 1992 等）。推广方向已有大量工作：非扭转映射（nontwist map）的剪切环面破缺（Chaos 2006）、三次/四次非扭转映射（CNSNS 2011）、体积保持映射的推广（arXiv:1205.6143）——这些确认判据可移植，也表明"换映射族测阈值"是该领域的标准研究单元。

空白在参数维度：检索未见双谐波扭转映射 (K₁, K₂) 平面上金环面破缺曲线的系统发表（两谐波会产生竞争共振结构；未找到不等于不存在，Ketoja–MacKay 一类重整化工作是否覆盖此族须第 6 周专项核查）。方法完全成熟公开，算力需求是高精度算术而非大规模，正合笔记本形态。

3 · 可检验假设

• H1 沿 K₂/K₁ ∈ [−0.5, 0.5] 的 ≥ 11 条射线上，破缺阈值可用 Fibonacci 周期轨道序列（周期至 F_n ≥ 10⁴）的剩余收敛/发散二分测定至相对精度 10⁻⁴，且 K₂=0 端点复现 Greene 值 0.971635 至 10⁻⁵。
• H2 临界曲线上各点的临界剩余值 R_c 与单谐波普适值一致（偏差 ≤ 5%）——或在某些射线方向偏离，提示重整化普适类改变（两种结果都成立）。

4 · 量化验收标准

1. 方法学校验（硬门槛）：自建管线在单谐波标准映射上复现 Greene 1979 结果:K_c = 0.971635 偏差 ≤ 10⁻⁵，且临界点剩余序列的自相似行为与 MacKay 发表图形一致。不过关则推广族结果全部无效。
2. 周期轨道求解:对称线方法（symmetry line）+ Newton 迭代，收敛判据 |fⁿ(x)−x| ≤ 10⁻²⁰（mpmath ≥ 30 位精度）；每条轨道剩余值报有效位数。
3. 阈值判定规则写死:连续 ≥ 5 个 Fibonacci 层级剩余单调发散判"破"、趋 0 判"存"，边界二分至 10⁻⁴；层级敏感性（用 F_n−1 截断复判）必须报告。
4. 收敛性检查:精度位数 30/50 双档、Fibonacci 截断两档，阈值漂移 ≤ 10⁻⁵ 才采信。
5. 独立交叉校验:≥ 3 个参数点用 Chirikov 重叠判据粗估 + Sobol' 网格轨道存活法（直接迭代看输运）验证阈值方向一致。
6. 代码开源、一键重跑。

5 · 数据与工具

用途	来源 / 工具
高精度算术	Python mpmath（任意精度内置）；Julia BigFloat 备选。周期轨道 Newton 全部自建（无现成包完整覆盖对称线方法——需核实 mpmath 下自动微分缺失，导数手推）
方法源头	Greene 1979 (J. Math. Phys. 20, 1183)、MacKay 重整化文献；对称线技术描述见 nontwist 系列（如 arXiv:nlin/0601040），仅用于校验与方法，不计入贡献
单谐波参照值	Kc = 0.971635406…（Greene/MacKay 发表值），仅校验
符号推导	sympy 推导双谐波映射的切映射与对称线
算力	单条 Fn ~ 10⁴ 周期轨道 Newton（30 位）秒级；单射线二分 ~10³ 条轨道 → 小时级；11 条射线 + 收敛档总 CPU < 100 小时

6 · 方法路径

1. 实现单谐波标准映射的对称线 + Newton 周期轨道求解（mpmath）。
2. 复现 Greene K_c 与临界剩余自相似（第 1 条验收）。
3. 第 6 周专项文献核查:双谐波族是否已被系统测过（Ketoja–MacKay、重整化文献链）。
4. 推导双谐波映射对称线结构（sympy 核验；对称线可能增多，须给可复现选线规则）。
5. 沿射线族二分测阈值，绘制 (K₁, K₂) 临界曲线。
6. 各临界点测 R_c 与剩余标度，检验普适性（H2）。
7. 三点独立交叉校验（轨道存活法），成文。

7 · 新颖性边界

本课题不声称新判据或新理论；Greene 判据、MacKay 重整化、以及在 nontwist/三次/四次/体积保持映射上的移植（Chaos 2006；CNSNS 2011；arXiv:1205.6143）均已发表——这同时证明"换映射族测破缺"是被承认的研究单元。本项目的贡献是双谐波扭转族临界曲线的测绘与 R_c 普适性的延伸检验，属"换构型 + 新参数扫描"。若第 6 周核查发现该族已被做（重整化文献覆盖两谐波扰动的可能性真实存在），差异化轴切换为:已发表工作未覆盖的旋转数（白银比例 √2−1）或负 K₂ 区——旋转数改变即样本改变，框架不动。若 H2 出现偏离，必须先排除数值假象（验收 4）后才可报告。历届丘奖 2025 finalist Smale 马蹄 Hausdorff 维数为理论构造，无重合。

8 · 决策门槛（go / no-go）

• 第 5 周末:单谐波 K_c 复现至 10⁻⁵。这是全题最硬的门槛，预计 4–5 周；失败最常见原因是 Newton 初值链断裂（沿参数延拓解决）。若第 8 周仍未过，放弃本题转 M15（同族动力系统技能可全部迁移）——这是唯一"换题级"退出，故本题排最后。
• 第 6 周末:文献核查完成，差异化轴定稿。
• 第 16 周末:双谐波对称线方案可用（周期轨道找得到）。若对称线结构混乱导致漏轨，降级路径 A：改用无对称线依赖的多点打靶 Newton（计算量 ×5，规模砍半:射线 11→5 条），临界曲线主结论保留、分辨率降档。
• 第 32 周末:若 R_c 测量散布过大无法判 5%，H2 改报区间估计，H1（临界曲线）独立成文。
• 适合对象:仅推荐给已展示高精度数值耐心与动力系统阅读能力、且接受"前 8 周可能零产出"的学生；回报是全组数学辨识度最高的结果之一。