一份给同行看的工程化复盘：agent 编排与失败恢复、文献检索的查询设计与去伪、每块物理的方程/数值方法/算法、四轮机制证伪与修复、以及全程的可复现实践。物理量按 measuredderivedestimatedecision 标注。
整个研究遵循一条编排纪律：主上下文只做协调、决策、综合，不做重读写检索的重活。所有实现、检索、文件写入都交给并行子代理（Task agent），保持主上下文窗口干净、可长程推进。
WebSearch/WebFetch 是「延迟加载工具」（deferred），没先用 ToolSearch 把 schema 载入就直接调用，结果空转、不写文件。诊断后在每个子代理 prompt 顶部加了强制的「STEP 0：ToolSearch query="select:WebSearch,WebFetch"」+「文件必须写盘才算完成」两条硬约束，重新派发即恢复。这类失败全部如实记进了 CHANGELOG。下列数字均从本会话的子代理 transcript（93 个后台任务文件）中实测汇总 measured。说明：这些是子代理侧合计，且包含「子代理递归再分裂」出的 helper；主控自身 token 未计入。cache_read 巨大是因为每个 agent 启动时会重读被缓存的上下文（这是正常的，且大部分走缓存、单价低）。
| 主模型 | Claude Opus 4.8 |
| 后台任务文件 | 93 个 |
| 子代理运行 | ~84（含递归 helper） |
| 后台 bash 任务 | 9（拟合/编译等） |
| output | ~1.01M |
| input（新） | ~4.45M |
| cache 读取 | ~157M |
| cache 写入 | ~17M |
| 子代理工具调用 | ~1,861 次 |
| 其中 Web 检索/精读 | 数百次 WebSearch+WebFetch |
| 主控 git 提交 | 20+ 次（研究日志） |
| 产出代码 | 13 模块 / 2310 行 |
口径说明：output ~1M token 是子代理「真正生成」的内容量；input 的「新 token」~4.4M 与 cache 读取 ~157M 不应简单相加当成本——缓存读取单价远低于新 input。把 cache 读取计入是因为多代理并行时每个 agent 都要把共享上下文读进来，这正是「多代理换上下文窗口」策略的代价与收益所在。
novelty 是论文的生死线，所以 Phase 0 是硬门槛：动任何推导前先把先验工作查穷尽。策略是四路正交分解——把"可能撞车"的空间切成互不重叠的子域并行扫，再交叉三角验证。
WebFetch 精读，而非只看搜索摘要。WebFetch 打开摘要页核对「标题+作者+年份+期刊」四项一致才采纳；否则标 [arXiv unverified]。检索中剔除了多条张冠李戴，例如：把 2412.19371 误归给 Yang Bai（实为 Rosenlyst）、把不存在的"Mohapatra-Nussinov box-diagram"当先例、把 1503.07874（其实是对撞机 EFT）误当 loop-DD。这些假归属若进了 bib，审稿人一眼能看穿。四路报告综合成一页判定 lit/00-novelty-synthesis.md：核心新意（单旋钮统一）真新；box portal 与 σ/m(v) 是支撑、部分被预占，逐条技术区分。最后做 INSPIRE 引用树终查确认优先权：
| 最近邻论文 | INSPIRE recid | 被引 | 结论 |
|---|---|---|---|
| 2602.18753（box portal 最像；第一作者本人前作） | 3122802 | 0 | 无后续，改为自引扩展 |
| 2410.23631（confining+loop DD） | 2844260 | 7 | 7 篇全为相邻但不同主题，未预占 |
→ 单旋钮 glueball↔暗π SIDM + SM-色 box portal + ADM 的组合未被占据。decision
定模型时第一个岔路是「几种暗夸克 N_f」。最初想 N_f=1（最简），但用 Goldstone 计数自检后否决了——这是个纯对称性论证，我现场用 sympy 算了一遍。
三项现场自检全过：Goldstone 计数、谱走向、SIDM 尺度自洽（取 Λ_d≈0.5 GeV → M_B≈1 GeV、m_π≈几十 MeV，落在标准 SIDM 窗口）。
连接子 {ψ,Φ}（带 SM 色 × 暗色）在一圈 box 图里跑。M_ψ,M_Φ ≫ Λ_d 做重质量展开；box 是紫外有限的，手征性禁止质量插入，匹配到维度-6 矢量算符：
数值（code/box_calc.py，主控复跑）：M_ψ=1.5 TeV、y=1 → σ_SI ≈ 4.66×10⁻⁴⁵ cm² derived。解耦论证：SIDM 强度由 (Λ_d, r) 定，核子耦合由 (y, M_ψ) 定——参数上正交。
⚠️ 第三轮审查修正：拟合把 M_B 推到 38 GeV（氙实验灵敏度峰值）后，该基准点被 LZ（σ_SI<9.2×10⁻⁴⁸ cm² @36 GeV measured，arXiv:2207.03764）排除 ~500×；y⁸ 陡标度下 y≲0.46 或 M_ψ≳7 TeV 即恢复可行。叙事从「躲开直接探测」翻转为「LZ 正在扫描 portal 参数空间」——更可证伪。
这块是整个工作的心脏，也是前两轮证伪的发生地。任务：从介子交换势出发，解出转移截面随速度的变化，做麦克斯韦平均，χ² 拟合矮星系→星系团真实数据。
lit/02。Yukawa 势的转移截面没有统一闭式，需按 (κ, β) 分区。我们把 Tulin-Yu-Zurek (1302.3898) 的 Born/经典拟合式 + Hulthén s-波（共振/量子区）逐式移植自同行评审的 CLASSICS 代码 (2011.04679)，并对齐其匹配常数：
# code/sigma_v.py — 区制分类（节选自真实实现） # kappa = m_phi/(M_B * v_rel) ; beta = 2*alpha_d*m_phi/(M_B*v_rel^2) def _sigma_combined(kappa, beta, sign): if kappa > 1.0: # 经典 / Born 区 return _sigma_classical(kappa, beta, sign) elif kappa < 0.4: # Hulthén s-波（共振/量子） return _sigma_hulthen(kappa, beta, sign, eps=1.6) else: # 0.4≤κ≤1 线性插值过渡 return _interp(kappa, beta, sign)
用 96 点 Gauss-Legendre 求积（np.polynomial.legendre.leggauss(96)）在麦克斯韦相对速率分布上算 ⟨σ_T v⟩/⟨v⟩，与数据论文约定一致；扫 v=10–3000 km/s。快速版供参数扫描、自适应 quad 版供单点高精度。
对称冻结用标准 Boltzmann ODE，变量 Y=n/s、x=M_B/T_SM，scipy.solve_ivp（BDF）求解，Gondolo-Gelmini 约定，并用 Steigman 2012 热遗迹基准校准（在 100 GeV 复现 Ωh²=0.121，误差 ~1%）。
暗π 是赝标量，单π交换给自旋-自旋/张量势，自旋无关中心 Yukawa 被 (m_π/2M_B)² 压低。code/ope_mediator_check.py 现场算出该因子 ≈3.4×10⁻⁷ derived ⟹ 轻暗π 单独不能当中心 mediator。
改试 s-道近阈共振（CGCM 1810.04709，sigma_v_resonant.py 实现并核验 BW 峰值到机器精度）。诚实拟合结果：单共振给 σ_T∝v⁻⁴，矮星系/星系团比 ~10³，而数据斜率仅 d ln(σ/m)/d ln v ≈ −0.70（~10 倍）；最优 χ²/dof≈8、E_R 顶到网格边界。否决。
真正对的中心力是各向同性-标量道的双π交换——核力中程吸引的暗版，自旋无关、不被压低、射程 ~1/(2m_π)。我现场验证：m_π~14 MeV → 有效介子 ~28 MeV，正是 Yukawa 拟合想要的；Yukawa 形状斜率 −0.75 对上数据 −0.70。r 旋钮故事完整保留，有效 α_d 当时被假定可由暗π 耦合经 TPE 导出——这个假定在第四轮被定量证伪，见下。
成稿后以敌意审稿人标准全局复审，发现：N_f=2 精确同位旋让 3 个暗π绝对稳定，而 box 是矢量流、不能让它衰变；WZW 3→2 在 N_f=2 恰好为零 ⟹ 28 MeV 稳定遗迹过闭合宇宙 ~10⁵ 倍（ξ=0.4）derived。同轮还发现 38 GeV 基准被 LZ 排除 ~500×、论文内部数字未随基准迁移同步。
Yukawa y_i ē_R q_L^i Φ_ℓ 积掉后诱导螺旋度压低的 π_d→e⁺e⁻：y=1、M=1 TeV 时 τ≈2.5 msderived，窗口 M/y≲4.5 TeV；暗π 与电子热接触到 T~m_π 后干净退场 ⟹ ΔN_eff≈0。DD 叙事同步翻转为「LZ 正在探测」；所有数字对 38 GeV 基准重新同步；暗核束缚态（αμ/m_φ≈116）作为显式 caveat 引 Wise-Zhang。深读 Cline-Perron 后在论文中如实让位其「非对称暗重子拟合同一数据」的先行地位，区分点收窄为：无调谐 TPE 势 vs 他们百分级调谐的 ERE 极点、有色 box vs 暗光子、r-dial+glueball 端点。
以 Nature Physics 审稿人标准再审时抓出更深的漏洞：π 是 Goldstone 玻色子，耦合全是导数耦合，在 m_π/f_π=0.0075 的基准点几乎自由。用手征 EFT 的单π/双π交换势（Kaiser-Brockmann-Weise NNLO；code/mediator_strength_check.py，先在真实 QCD 上校验：1 fm 处 −224 MeV、OPE f²/4π=0.072）算暗扇区：α_eff = 5.5×10⁻¹¹ derived，比需要的 0.34 小 6×10⁹ 倍。并由此得到一个普适 no-go：晕尺度速度依赖要求 m_med/M_B~10⁻³，若 mediator 是同一禁闭扇区的手征 pNGB，必然 m_π/f_π≪1 而退耦——任何禁闭暗扇区都不能用自己的π做轻 mediator。这也解释了 Cline-Perron 为何只能靠百分级调谐的极点。
近共形（walking）规范理论的标度对称性赝 Goldstone——dilaton——耦合 ∝ 质量：L ⊃ (M_B/f_d) d B̄B ⟹ α_d = M_B²/(4πf_d²)。初版取 f_d=Λ_d 得 α=0.36、与拟合 0.34 相差 6%——这在重建阶段被证明是运气：真实格点比值（SU(3) N_f=8：M_B/F_π≈11、f_d/F_π≈3.3，LSD/LatKMI/dilaton-EFT 拟合 measured）给 α_d = 0.9（0.67–1.2）。更重要的发现：把 α 钉在 0.34→1.2 任意值，χ² 几乎不变——数据测的是 (m_d, M_B)，对 α 不敏感。所以正确表述是：α 是格点输入，理论只剩两个参数拟合五个跨尺度数据点，χ²/dof=1.9（m_d=65 MeV, M_B=38 GeV）derived。dilaton 质量 m_d²∝m_q（Golterman-Shamir；Appelquist-Ingoldby-Piai）⟹ r 旋钮完整保留；格点 SU(3) N_f=8（LSD, PRD 99 014509）证实轻标量与π近简并。拟合数值一字不改（Yukawa 形式相同，仅 m_med→m_d）。代价与诚实：模型从 N_f=2 改为近共形 N_f（参考 N_f=8）；"dilaton 介导暗重子 SIDM"的思想须归功于 Zihan Wang（2601.04340，独立作者）——我们的增量是 no-go 定理、带 glueball 端点的 dial、参数无关 α 吻合、真实多尺度拟合、有色 box+LZ、ADM。
用户要求的不是"能发表"，而是"工作本身要好"。据此做了一次不留情面的再评估："一个旋钮"是描述性的包装（两端点是教科书，中段没算过非平凡物理）；真正有深度的是三样东西，围绕它们重建了脊柱（theory/07-reconstruction.md，三路并行文献研究 lit/05–07，约 270 次核验抓取）：
晕尺度速度依赖要求 m_med/M_B~10⁻³；对手征 pNGB 这迫使 m_π/f_π~10⁻²、Goldstone 退耦——π 交换比所需低 5–9 个量级（图 pion_nogo.pdf，含真实 QCD 参照点）。假设明示：θ_d=0（否则 θ 项给π标量耦合，2607.07406）、无微调极点（Cline-Perron 2%）、无额外轻场。⟹ 观测到的速度依赖蕴含近共形动力学——INSPIRE 核验无先例。
dilaton 耦合于质量：α=M_B²/(4πf_d²)。用真实 SU(3) N_f=8 格点值 ⟹ α≈0.9(0.67–1.2)；数据对 α 不敏感、只测 (m_d, M_B)；两参数拟合 χ²/dof=1.9。Zihan Wang 2601.04340 的基准隐含 M_B/f_d=500，与 dilaton EFT 自相矛盾——我们的构造自洽。
近共形理论的禁闭相变是过冷一级相变（α≫1）——与通讯作者自己的非过冷暗禁闭引力波工作（β/H~10⁴）恰为相反角落。重建阶段的粗估（T_reh∈[Λ/30,Λ/2] ⟹ f~10⁻⁶–10⁻⁴ Hz、"到 LISA 低频边缘"）已被 7.6 节的完整相变计算取代：T_reh=201 MeV=Λ_d/90、β/H≈450，峰值 2–4 μHz、h²Ω≈5×10⁻¹²，LISA 与 PTA 均不可见。诚实条件：N_f=8 有限温相变阶数格点未定；过冷稀释不对称 D≈1200，把遗迹机制和 T_reh 绑定。
7.5 的"引力波印记"是量级估计（T_reh 取 [Λ/30, Λ/2]、β/H 取 10–100）。这一轮把它真正算了一遍：用格点锚定的 dilaton 有效势做有限温相变，解泡核化、渗流、再加热和引力波谱。结果把叙事从"LISA 低频边缘"诚实降级到 μHz，但同时把"相变能否完成"从假设变成了一个有机制的结论。theory/08-dilaton-transition.md、lit/08-dilaton-transition.md、code/dilaton_transition.py。
solve_ivp DOP853，rtol 10⁻¹¹），O(3)（S_3/T）与 O(4)（S_4）两个通道。验证：Fubini–Lipatov 解析弹跳 S_4=261.76 对 8π²/(3λ)=263.19（差 0.5%）；独立 solve_bvp 边值解在倾斜双阱上复现 S_3、S_4 到 4 位有效数字；薄壁极限趋势单调 derived。| 量 | 值 | 含义 |
|---|---|---|
| 真空能 ΔV | 0.061 GeV⁴（ΔV^{1/4}=0.50 GeV） derived | 由 m_d²f_d²/16 定，不是 Λ_d⁴——walking 理论的真空能被小 β 函数压低；这正是 7.5 粗估高估 T_reh 的根源 |
| T_c | 272 MeV = Λ_d/66 derived | Δ=3.06、自由气体、c_h×3/4 各变体：273–293 MeV |
| 纯 dilaton 势（去掉夸克质量项） | S_4≈3×10⁵ ≫ 166，永不核化 derived | Fubini 估计 S_4≈8π²/3λ_eff，λ_eff∝(m_d/f_d)²=3.3×10⁻⁵；完成需 m_d/f_d≳0.18（扫描阈值 0.166–0.187），我们的 5.8×10⁻³ 低 30 倍 ⟹ 退禁闭相永久暴胀 |
| 夸克质量项触发 | χ^y 项为负且相关，零温无势垒；热势垒 ∝T⁴、倾斜 ∝T²·T^{0.07}，spinodal 点≈4 MeV derived | von Harling–Servant QCD 触发（1711.11554）的 dEFT 版：4n→y=2.07。给 dilaton 质量的同一把旋钮，就是过冷的出口 |
| T_n | 19.1 MeV（T_n/T_c=0.07，N_e=2.7 e-折） derived | O(3) 热隧穿（χ_r=2.4T，厚壁），比 spinodal 高 5 倍，是真隧穿而非滚落；渗流 T_p=18.5 MeV，3+T dI/dT=−874 ⟹ 完成 |
| β/H_* | ≈450（细网格 496 / 粗网格 482 / R_* 定义 421） derived | 近 spinodal 的 S_3/T 塌缩：31→21→18→13 MeV 时 700→230→142→53 |
| α | 1.2×10⁴ derived | 极强过冷 |
| T_reh | 201 MeV = Λ_d/90（g_*=113） derived | T_reh/T_c=0.74；7.5 假设的下限 Λ/30=0.6 GeV 都没达到 |
| 稀释因子 D=(T_reh/T_n)³ | 1.2×10³ derived | 相变前任何不对称被稀释三个量级 |
| 引力波峰 | f=2–4 μHz，h²Ω≈5×10⁻¹²（envelope 3.9 μHz/5.1×10⁻¹²；bulk-flow 2.2 μHz/4.2×10⁻¹²） derived | 声波通道 Υ≈0.007 可忽略；LISA 在 3.9 μHz 的噪声 h²Ω≈2×10⁻⁵，1 mHz 处 f⁻¹ 尾≈10⁻¹⁴ 仍低于其 10⁻¹¹–10⁻¹³ 底噪 derived；NANOGrav pt-bubble 窗需 β/H_*<10.5（H_*R_*>0.28） measured，我们 H_*R_*=6×10⁻³ ⟹ 出窗 |
| μHz 任务可及性 | μAres 应变极限≈6×10⁻¹³（1908.11391 h_c≈10⁻¹⁸ @1 μHz 换算）；小行星测距≈2×10⁻¹² @2 μHz derived | 信号高出 μAres 约三个量级；LLR/SLR 轨道共振探针（2038 预测 ≈10⁻⁹）差 ~10³ measured |
敏感性扫描（theory/08 §6.1）：g_χ∈{2.4,4.19,7.26} 给 T_n=17–19 MeV（核化发生在 spinodal 之上 5 倍处，势垒由 T⁴ 驼峰而非热质量决定，故对 g_χ 不敏感）；自由气体方案 25 MeV；Δ=3.06 与 4 差 <1%；唯一敏感参数是 y——y=2.0→22 MeV，y=2.6→3 MeV（α~10⁷）。所有变体 β/H_* 落在 250–630、T_reh 190–204 MeV、峰 1–3 μHz、h²Ω 3×10⁻¹²–2×10⁻¹¹，没有一个进入 NANOGrav 窗或 LISA 视野。
文献代理：直接抓取 27 篇 arXiv PDF 逐页核对公式（dEFT 势的符号与系数、热自由能、核化前因子、GW 谱公式、μHz 探测器灵敏度、稀释与共生成机制），产出 lit/08（含"可直接进代码的公式"节）。计算代理：写并验证弹跳解算器，跑物理点 + 10 个敏感性变体 + m_d/f_d 七点扫描，单次全扫约 10 分钟；产出 theory/08 与 dilaton_transition.pdf。
长时间数值扫描期间两个代理都出现过空转——等待后台 bash 完成时不再推进也不汇报。主控按固定间隔检查 transcript，发现停滞即发消息"nudge"（要求先写盘已有结果、再继续），两次即恢复。教训与第 1 节相同：长跑任务必须让代理边跑边落盘，主控不能把"没消息"当"在工作"。
最终拟合（code/fit_tpe.py，Yukawa 形式对 dilaton 交换同样成立）：对有效参数 (α_d, m_d, M_B) 用 scipy.differential_evolution 全局优化，χ² 含 5 个双边数据点 + 星系团单边上限罚项。结果：
内建诚实性检验（写在 fit_resonant_fast.py 里）：
每完成一个 block/判定就提交，commit message 记录物理结论与转折（含「resonant rejected」这种）。17 个提交串起完整路线。
TODO/CHANGELOG/ARCHITECTURE + memory 文件，保证上下文压缩或换会话后不丢线索；子代理产出由主控落盘。
13 个 Python 模块全部主控复跑验证；每个带单位的物理量标 measured/derived/estimate；格点/实验输入带 arXiv。
pdflatex && bibtex && pdflatex 链在单条 shell 里会间歇导致 .bbl 没读入、出现 [?] 断引用。诊断后写了 paper/build.sh，把四步拆成独立进程调用，并改用 pdftotext | grep '[?]' 验证（而非数累积 log）。这条教训也写进了脚本注释。