Quantum Work Network — 有限条件での数値観察記録

1. 要約

何を調べたか

小さな量子ネットワークにパルスでエネルギーを流し込んだとき、受け手（load、負荷）にどれだけ「取り出せる仕事」が届くか、そしてそれが位相ノイズ（dephasing）によってどう変わるかを、数値的に計算しました。実装はRustで、Milestone 1からMilestone 15Aまで、前の段階のコード・CSV・レポートを書き換えずに積み上げる形で進めています。

どの範囲まで調べたか

• 鎖の長さは N=3 から N=7 まで（N=8 は正式な鎖長比較に含めていません）。
• 時間は基本的に t=0 から t=10 まで。M13Aの分岐実験と M15A の抽出は t=8.50 まで。
• ノイズの与え方は fixed-per-site（1サイトあたりを固定）と fixed-total（合計を固定）を明確に区別しています。
• M15A のあと、補足としてC2プロトコルを t=10 まで延長し、長時間窓のC0基準と比較しました。これは新しいMilestone番号を与えていません。

何が観察されたか

• equal-input（投入エネルギーをそろえた）条件で、N=7 は N=3 より高く・遅く・細いピークを持ち、一方で 0→10 の W 時間面積は小さくなりました（M11g）。
• matched N=7 条件では、loadのergotropyの対角（population）成分がすべての保存時刻で 0 であり、W = W_coh と分解されました（M12A）。
• 同じ状態から出発し、分岐後に鎖-load結合を半分にすると、有限窓のなかで W_coh の減少が弱まりました（両方の分岐開始時刻で）（M13A）。
• t=8.50 での1回の理想的な局所抽出では、C2 の gross work が C0 より大きく、その差は抽出前の W_coh 差と数値精度の範囲で一致しました（M15A）。

何はまだ分からないか

• これらの観察の背後にある物理機構。因果関係は一切確認していません。
• 量子的優位性、古典手法に対する優位性、一般的なスケーリング則。
• 実際の制御装置のコスト、有限時間のスイッチング、繰り返し動作、t=10 より先の挙動。
• 「鎖が長いほど良い」という読み方。M10の指標順位は指標ごとに向きが異なります。

2. モデルと測定量

設定はMilestone 1で凍結され、以降維持されています（hbar=1）。

系の構成

2準位のサイト（qubit：|0>＝空、|1>＝励起の2状態を持つ最小単位）を N 個つなげた鎖と、その一方の端に付けた1個の3準位 load（負荷。エネルギーの受け手）からできています。既定値は、サイト間結合 J=1、サイト-load結合 g=0.25、角振動数はすべて 1。鎖のオンサイトハミルトニアンは omega * sum_i |1><1|_i。駆動するのはサイト0（一方の端）で、loadが付くのはサイト N-1（反対の端）です。

drive（駆動）

driveは、外部から当てるパルスのことです。エネルギーを鎖に注入する入口の役割をします。

包絡線は f(t)=sin^2(pi t / tau)（0<=t<=tau、それ以外は0）、既定値 tau=3.2、Omega=0.2 です。t > tau = 3.2 ではdriveは厳密に0になります。

dephasing（位相ノイズ）

dephasingは、量子的な位相の関係（コヒーレンス）を壊す種類のノイズです。各鎖サイトに L_phi,j = sqrt(gamma_phi/2) sigma_z,j として与えます。loadには直接のノイズを入れていません。時間発展はLindblad型のマスター方程式です。

local ergotropy（取り出せる仕事）

ergotropy（エルゴトロピー）は、ある量子状態からユニタリ操作だけで取り出せる仕事の最大値です。ここではloadの縮約状態から計算します。

L = -i (I kron H - H^T kron I)
    + sum_k [ L_k* kron L_k
            - 1/2 I kron (L_k^dagger L_k)
            - 1/2 (L_k^dagger L_k)^T kron I ]

量の読み方（重要）

• W_time_area / E_time_area は状態量の時間積分であり、累積抽出仕事ではありません。
• usable_fraction は ergotropy を load エネルギーで割った記述的な比であり、独立した性能利得ではありません。
• W/Ein は全体効率ではありません。
• XGamma は dephasingカーネルで重み付けしたコヒーレンス露出の診断量 x_gamma(t) = sum_ab Gamma[a,b] |rho[a,b](t)|^2 の時間積分です。損失仕事、散逸エネルギー、熱、エントロピー生成、効率、損傷のいずれでもありません。
• W_coh / W_pop は load ergotropy のエネルギー基底での記述的な分割であり、独立に保存される量ではありません。
• 相互情報量（MI）は全相関の観測量（単位はnat）であり、仕事でもergotropyでもentanglementでもありません。
• switching work は理想的な瞬間ハミルトニアン跳びの帳簿であり、アクチュエータや制御装置のコストではありません。
• state / switching / extraction の3つの帳簿は別々に報告し、足し引きして正味利益にすることはありません。

3. Milestoneの流れ

各段階の数値の詳細は個別の MILESTONE_*.md レポートとCSVにあります。ここでは進行の流れと代表的な結果だけを並べます。

4. 現在の主な観察

以下は3種類に分けて示します。ラベルの文字（観察記述的サポート未確認）でも区別できます。色だけに頼らずお読みください。

5. M12〜M15Aの観察の流れ

以下は観察が得られた順序に沿って並べています。この並びを因果の機構として読まないでください。どの段階でも、原因と結果の関係は確認されていません。

(1) W = W_coh という有限条件での分解結果（M12A）

M12はM11Mのmatched N=7条件を固定したまま（TOTAL_GAMMA=1.5、Omega = 0.18748395731510084、dt=0.0025、t=0..10、1001保存点、真空初期状態）、loadの取り出せる仕事のピーク後の減少について3つの記述的な問いを立てます：どの成分が変わるか（12A）、どんな時間順序か（12B）、その周りでload-鎖の相関はどう振る舞うか（12C）。マッチングのやり直し、再最適化、根の探索、細dtの実行、他のN/gamma/Omegaはいずれも行っていません。

12Aでは、同じ固定条件で新しい軌道を1本走らせ、既存のエネルギー基底の規約（H_L = diag(0,1,2)）で load ergotropy を W = W_pop + W_coh に分けます（W_pop = W(Delta_H(rho_L)) が対角＝population部分）。すべての保存時刻で W_pop = 0、したがって W = W_coh。t = 7.70 → 10.00 で Delta W = Delta W_coh = -1.0587822145716654e-3、Delta W_pop = 0、Delta E_load = -1.8634507496082930e-4、Delta E_passive = 8.7243713961083613e-4。観察された符号付きのWの減少は、まるごと W_coh の減少として記述されます。

ただし事前登録していた比のルールは両方の成分が正の値で減少することを要求しており、L_pop = 0 により比が定義できなくなったため、正式な分類は component_change_inconclusive です。この判定は事前確約したルールが適用できなかったことを表すもので、分解が失敗したという意味ではありません。分解自体は厳密です（max |W - (W_pop + W_coh)| = 0）。数値品質20/20 PASS、M11M集約8/8 PASS（最大絶対差0）、正式入力SHA-256 7/7不変。population部分とcoherent部分は独立に保存される量でも、別々の物理エネルギーでもありません。

(2) 時間的な前後関係（M12B・M12C）

12Bは新しい時間発展なしで、12Aの保存1001点の系列だけを、事前登録した5区間の持続ルール（M12B_PRECOMMIT.md。平滑化、フィット、FFT、極値の補間なし）で読みます。保存格子上の順序は次のとおりです。

正式分類 coherence_and_ergotropy_decline_before_negative_current（3/5/7区間の感度はすべて一致）。持続的な負の電流は W_coh の減少開始より 0.46 遅く始まるので、「電流が負に転じたことが W_coh の減少を始めさせた」という単純な説明は時間順序から支持されません。これは順序についての記述だけであり、電流と W_coh のあいだにどちらの向きの因果関係も確立されていません。チェック15/15 PASS、W = W_coh の最大差0、入力SHA-256 3/3不変。電流の積分は状態由来量の時間積分であり、仕事、効率、損失、散逸、熱ではありません。

12Cは同じ固定条件で新しい完全軌道を1本走らせ、load-鎖の相互情報量 I(C:L) = S_chain + S_load - S_total（自然対数、nat）、個々のエントロピー、purityを全1001保存時刻で評価します。MIの大域ピークと5区間の持続的減少はどちらも t = 7.51 にあり、ピーク値は 4.5985207125390781e-3 nats。これは t=7.70 の W_coh ピーク/減少より 0.19 早い。正式分類 mutual_information_peaks_before_Wcoh_decline（3/5/7区間の感度格子で不変）。数値品質29/29 PASS、既存の正式成果物のSHA-256 30/30不変、M12Aとの共通量は最大絶対差0で再現。相互情報量は古典相関と量子相関を含む全相関の観測量であり、仕事でもergotropyでもentanglementでもありません。MIの減少が W_coh の減少を引き起こしたとか、W がMIに変換されたといった主張はしていません。

(3) couplingを半分にした短時間branch（M13A・M13B）

13Aは、M12Cと同じ保存済みの全系状態から出発し、分岐時刻以降の生成子だけを変えた短い分岐を、0.01の保存格子で t=8.50 まで伝播させました。分岐開始は t=7.50 と t=7.70 の2つ（正式保存行は合計728行）。

• B0 基準の継続（M12Cを厳密に再現、最大絶対差0）
• B1 分岐後に鎖のdephasingを0にする
• B3 分岐後に鎖-load結合を半分にする（g: 0.25 -> 0.125）
• B4 結合を0にする境界コントロール（最大不変量ドリフト 1.9984014443252818e-15）

B3は両方の分岐開始で同じ正式分類 weaker_coupling_weakens_Wcoh_decline を与えました。

B1は開始時刻で結果が異なり、平均していません（7.50では dephasing_counterfactual_mixed、7.70では future_dephasing_change_has_small_effect）。B4は両方の開始で zero_coupling_boundary_control_pass。数値品質の行はすべてPASS（最大トレース誤差 1.776e-15、エルミート性誤差0、フォールバック不使用、分岐後のdriveは0）。

許される要約は正確に次のとおりです：同じ状態から出発し、分岐後に鎖-load結合を半分にすると、通常の結合を保つ基準に対して、有限窓のなかで W_coh の減少が一貫して弱くなった（両方の分岐開始で）。明示的に主張しないこと：結合が W_coh 減少の唯一の原因であること、弱い結合が常に良いこと、最適な結合が見つかったこと、結合を半分にすることが何らかの機構で W_coh を守ると示されたこと。B1とB3のパーセントは、どちらの操作が「より効く」かを比べる形では対照していません。W_coh の時間面積は状態量の時間積分であり、累積抽出仕事ではありません。

ここまでの位置：M12はピーク後の W_coh 減少について、成分（W = W_coh、12A）、時間順序（12B）、load-鎖相関（12C）を記述しました。M13Aは同一状態からの生成子変更に対する有限時間の応答を調べ、両方の分岐開始で結合を半分にすると W_coh の減少が弱まることを見いだしました。M13Bはその分岐の内部応答順序を、利用できる保存解像度では一意に決められませんでした。結合への感度は確認された観察であり、具体的な物理機構と最適条件は未決定のままです。

(4) switching-work accounting（M14）

M14は、M13Aの結合変更を行うのにモデルのハミルトニアン帳簿の内側でどれだけの仕事がかかるかを問い、その帳簿をstateの帳簿と厳密に分けます。以降で使う固定のプロトコル名：

g=0、追加の結合値、スイッチバック、ramp、フィードバック、複数回スイッチはいずれも認めていません。

14A（設計のみ）は設計と監査のMilestoneで、正式な制御軌道も、結合の掃引も、最適化も、rampも実行していません。実装された相互作用項（H_CL(g) = g V_CL、V_CL = sigma_plus_N b + sigma_minus_N b_dagger、g に厳密に線形）を監査し、driveの包絡線が t > tau=3.2 で厳密に0であること（したがってスイッチ後のプロトコル差はdrive投入の差ではないこと）を確認し、理想quench規約 W_switch_on_system = Tr[rho(t_s)(H_after - H_before)]（正＝外部制御器からモデル系へ供給。quenchを挟んで状態は不変）を固定し、3つの分離した帳簿（状態量、switching work、スイッチ後エネルギー）を固定しました。これらは合算も減算もされません。伝播なしの単位チェックは 11/11 PASS。再利用の監査では、M13Aの成果物が完全な密度行列も Tr[rho V_CL] も保存していないことが分かり、判定は single_prefix_recomputation_required（基準prefixを1本だけ再実行）。最終の設計判定は coupling_control_experiment_ready、推奨は Plan B。

14 Plan B（実行）では、基準prefixを t=0 から t=7.70 までちょうど1本伝播させ、t=7.50 と t=7.70 で完全な密度行列を保存し、理想的な同一状態スイッチを両方で評価しました。新しいスイッチ後の分岐、結合値、掃引、ramp、抽出、最適化、細dt実行、t>8.50 の伝播はすべて 0。M13AのC0/C1/C2のスイッチ後状態軌道はそのまま再利用しています。

どちらの理想quenchのエネルギー跳びも固定の 1e-12 監査許容差を下回るので、正しい記述はこの許容差のもとで数値的にゼロです（状態 numerically_zero_within_1e-12_audit_tolerance）。ハミルトニアン跳びの仕事が数値的にゼロであることは、実際の制御器が無料であることを意味しません：このモデルにはアクチュエータのハミルトニアン、有限のramp、帯域幅、装置コストの帳簿がありません。M13Aの状態値を再利用した正式な複合分類は、C1・C2の両方で Wcoh_improves_with_load_energy_tradeoff + switching_accounting_available。端点 W_coh と W_coh 面積がそれぞれ1%以上改善する一方、端点のloadエネルギー保持率は固定の99%閾値を下回ります（C1：W_coh 端点 +6.142%、面積 +1.467%、保持率 97.159%。C2：W_coh 端点 +6.147%、面積 +2.103%、保持率 98.961%）。数値チェック 22/22 PASS、M13AのB0量は両スイッチ時刻で最大絶対差 3.839e-13 で再現、ソース/正式入力のハッシュ15/15不変。報告されたswitching workの比はスケールの診断量であり効率ではありません。switchingの帳簿は、利益や「正味で守られた W_coh」や制御効率の主張には変換されません。

(5) t=8.50での単一回のideal extraction（M15A）

M15AはC0とC2だけを 7.70 から 8.50 まで伝播させ、単一の固定評価時刻 t = 8.50 で1回の理想的な局所抽出を行います。正式な比較は C2 − C0 のみ。判定 state_improvement_preserved_after_extraction、必須の数値失敗は 0。

gross-workの差は抽出前の W_coh の差を残差 9.1072982488782372e-18（比 1.000000000000）で再現します。相互作用項の跳びの差が 2.5412084893412970e-14 と小さいため、bare-loadの会計と全ハミルトニアンの会計は比較について同じ符号を与えます。数値と来歴のチェックは 40/40 PASS（トレース、エルミート性、positivity、抽出後のergotropyゼロ、M13A端点の再利用、M14状態のSHA-256、switching会計の再利用を含む）。

許される結論は正確に次のとおりです：t=8.50 において、C2に対する1回の理想的な局所抽出は基準C0より大きなgross extracted workを与えた。C2−C0のgross-work差は抽出前の W_coh 差と数値精度の範囲で一致した。相互作用項を含む全ハミルトニアン会計でも比較の符号は変わらなかった。つまり、状態レベルの改善はこの1回の理想的な局所抽出を通じて保たれました。

帳簿は分けたままです。W_coh、loadエネルギー、switching work、W_load_gross、W_operation_out_total を単一の「利益」に足し合わせることはしません。M14 Plan B のC2のswitching work（8.8624029002963674e-16）は参照のために併記しているだけで、どちらの抽出仕事の列とも組み合わせません。15Aで行っていないこと：新しい結合条件、掃引、ramp、最適化、フィードバック、測定、複数回の抽出、他のN/gamma/Omega、t=8.50 を超える伝播。ここにあるのは、最適制御でも、仕事のただの増加でも、効率の改善でも、確認された正味利益でも、制御器コストがゼロという主張でも、実験ハードウェアの優位でも、量子的優位性でもありません。

(6) C2のt=10延長（補足。新Milestoneではありません）

この補足比較は、t=8.50 で観察された C2 − C0 の bare-load ergotropy 差が、その窓の端点による見かけのものかどうかを問います。新しいMilestoneではありません。M14 Plan B で正式保存された t=7.70 の全系状態を、C2（g=0.125）で 1回・連続に t=10.00 まで伝播させました。t=8.50 での抽出なし、スイッチバックなし、再マッチングなし、正式軌道の分割なし。固定条件：N=7、TOTAL_GAMMA=1.5、サイトあたり gamma 1.5/7、Omega=0.18748395731510084、再開時刻 7.70、終了時刻 10.00、内部 dt=0.0025、保存間隔 0.01、内部920ステップ、231保存点、driveは厳密に0。C0の長時間窓基準は M12A のCSVで、明示的な成果物監査で固定しました。主分類 C2_advantage_positive_at_t10。

t=8.50 以降、差は150個の保存区間すべてで増加し（減少0、ほぼ平坦0）、ゼロ近傍に達することも負に転じることもありませんでした。差の大域最小は t=8.50、大域最大は t=10.00 で、内部に局所極値はありません。これはC2−C0の差についての記述であって、C2そのものについてではありません。C2は t=8.50 以降に新しいWピークを作っていません。窓内のWの最大は t=7.70 のまま（3.3852501212993710e-3。C0の最大と同じ値）で、差が広がったのはC0の減少の方が速いためです（t=10 の保持率はC2が 87.807172%、C0が 68.723664%）。7.70..10.00 の ergotropy 時間面積は A_C0 = 6.5655097929267292e-3、A_C2 = 7.2916111984482872e-3、A_Delta = 7.2610140552155800e-4、相対差 +11.059330%。これらは瞬時の bare-load ergotropy の時間面積であり、累積抽出仕事ではありません。t=10 での付随するエネルギー差は Delta E_load = +2.106309466093e-5、Delta E_passive = -6.249614096404e-4、エネルギー恒等式の残差は0。site7→loadの電流は既存の定義のもとでの観察としてのみ記録され、ergotropy差の説明には使っていません。

監査：必須の数値チェック 28/28 PASS、既存成果物の再現チェック 15/15 PASS、SHA監査 PASS。C0の短窓/長窓の観測量の残差は0、C2の t=7.70..8.50 の再現残差は0、M15A の t=8.50 の値は最大残差 9.107e-18 で再現します。固有値ソルバのフォールバック数0、正式ランタイム 682.479 秒。

範囲：これは1本の保存状態の1回の伝播を、1つの固定基準に対して比べた、有限条件・保存格子上の比較です。この延長では抽出を行っていません。プロジェクト内の正式な抽出は、いまも M15A の t=8.50 における1回の理想的な局所抽出だけです。これは、一般の弱結合やC2の優越、最適な結合・スイッチ時刻・抽出時刻、いかなる因果機構、passiveエネルギーのergotropyへの変換、switchingと抽出の仕事を合わせた正味利益、繰り返しや周期動作、定常出力、t=10 より先の挙動、他のN/gamma/Omega/g、効率の改善、量子的優位性のいずれも確立しません。ここで advantage という語が使われている場合、それはこの固定された有限条件のもとでの bare-load ergotropy 差の符号だけを意味します。

6. 検証と再現

観察検証済みの部分

各段階のレポートは cargo fmt のPASS、解析バイナリのPASS、そしてリリーステスト数を記録しています（例：M2は26テスト、M5aは 47 passed / 0 failed / 1 ignored、M10bは101、M10cは104、M10 Finalは107、M11cは110、M11dは113、M11eは116、M11f と M11h–M11k は119、いずれも 0 failed / 1 ignored）。9c validationの実行時チェックは 47/47 PASS、頑健固有値診断の必須チェックは 10/10 PASS。11j と 11k は解析の前後で入力をSHA-256で検証しています。

後半の段階もそれぞれのランタイムと監査を報告します。M14Aの設計チェックは時間発展を必要とせず（ビルド後 3.455 秒、11/11 PASS）。M14 Plan Bは3080ステップの基準prefixを t=7.70 まで 2177.4 秒（36.29分。構築 5.2 秒、prefix 2169.8 秒、2状態の診断/書き出し 2.4 秒）で伝播し、22/22チェックPASS、15/15ハッシュ不変。M15Aは 40/40 の数値・来歴チェックPASSを報告します。Plan Bで保存された2つの 384×384 状態は、行優先リトルエンディアンの complex-f64 バイナリ（QWNRHO1 メタデータ付き、各 2,359,320 バイト）として保存され、m14_plan_b_state_files.csv で個別にハッシュされています。補足のC2延長は正式ランタイム 682.479 秒（11.37分）で、必須数値チェック 28/28 PASS、既存成果物の再現チェック 15/15 PASS、SHA監査PASS、固有値ソルバのフォールバック数0。抽出していないC2状態2つ（c2_pre_extraction_state_t850.bin、c2_pre_extraction_state_t1000.bin）を同じ QWNRHO1 形式で、それぞれSHA-256を記録して保存します。

再現ZIPの役割

完全な再現パッケージは、リポジトリの zip/ ディレクトリにMilestoneごとのZIPアーカイブとして置かれています。ソースコード、レポート、CSV出力、監査ファイル、その段階に必要な保存状態は、対応するアーカイブの中に入っています。GitHubリポジトリのルートに個別ファイルとしてすべてが展開されているわけではありません。

該当するアーカイブをダウンロードして展開し、その Cargo.toml があるディレクトリから以下を実行します。Rustツールチェーンが必要です（Cargo.toml：edition 2021、依存 nalgebra、num-complex、thiserror、dev依存 approx）。アーカイブの中身は段階によって異なるため、下に挙げたバイナリやデータファイルは、その段階のパッケージにしか無いことがあります。

cargo fmt --all -- --check
cargo test --release --offline
cargo build --release --offline
# opt-in 24-dimensional smoke test:
cargo test --release full_24d_short_time_smoke_test -- --ignored --nocapture

cargo run --release --offline --bin time_dependent_sanity          # M5a
cargo run --release --offline --bin dephasing_kernel_benchmark     # M8c
cargo run --release --offline --bin n7_noise_free_full             # M9a
cargo run --release --offline --bin n7_all_site_noisy_full         # M9b
cargo run --release --offline --bin fixed_total_noise_comparison   # M9c
cargo run --release --offline --bin n7_t002_eigen_diagnostic       # M9c diagnostic
cargo run --release --offline --bin n7_fixed_total_validation      # M9c validation (source of truth)

未確認検証で分かっていない部分

上記のチェックはいずれも、この実装が指定の条件で指定の値を出したことの確認です。物理的な機構、因果関係、他条件への一般化、実デバイス上での挙動については、何も確認していません。詳しくは次節をご覧ください。

7. 限界と未確認事項

このリポジトリが提供しないもの

ここにあるのは「この物理モデル、これらの初期条件、この有限時間、このRust実装について直接観察された値」です。とくに、このプロジェクトは以下を行いません。

• いかなる量子的優位性、あるいは量子が古典手法に勝つことの証明。
• 普遍的なスケーリング則、指数則/べき則の減衰、熱力学極限、N→無限大の挙動の確立。
• 長い鎖が一般に有利であるという主張、ノイズ位置や保護強度の普遍則。
• 曲線変換のフィットや減衰フィットを物理機構の証明として扱うこと（これらは記述的な形の診断にすぎません）。
• W_time_area を累積抽出仕事として、usable_fraction を独立の性能利得として、XGamma を損失仕事・散逸エネルギー・効率・損傷として扱うこと。
• レポートに書かれていない因果的説明の主張。

未確認What has not been confirmed

• 量子的優位性/量子超越性、あるいは古典手法に対する優越。
• 普遍的スケーリング則、指数則/べき則の減衰、熱力学極限、N→無限大の挙動。
• 長い鎖の一般的な有利さ、ノイズ位置/保護強度の普遍則。
• 観察された形の背後にある物理機構（群速度、境界反射、モードのうなり、entanglement/相関のフロント、いかなる因果関係も）。また変換/減衰フィットが特定の物理過程に対応すること。
• 抽出の制御コスト、有限時間のスイッチオフ/抽出、繰り返しサイクル、連続供給、定常出力、長時間安定性。
• 11Lで確認した1本のmatched N=7軌道を超える dt 半分化の収束（任意のdt、他のmatched条件、細dt固有のマッチング根）。マッチング根の一意性/単調性。
• TOTAL_GAMMA=3.0 でのequal-input結果、他の Omega 根、N>7、t>10。N=8は正式な鎖長比較に含まれていません。他の外部ソルバcrate（9c validationの範囲外）。（equal-inputのN=4/5/6は、いまのところ有限鎖のイベント構造診断としてのみ存在し、スケーリング則ではありません。）
• いかなるparity機構、普遍的な奇偶効果、parity秩序変数、統計的に有意なparity差。N=3..7 の odd_even_alternation_candidate は記述的な有限サイズの候補にすぎず、N=7はカテゴリ的な奇数群の共通性を破ります。
• 11Mのサイト分解診断の因果的な読み方（電流の変化がWの減少を引き起こす、電流のゼロ交差/反転が境界反射を証明する、passiveエネルギーの増加が熱化やデコヒーレンスを証明する、のいずれも）。
• 弱結合分岐の応答の内部時間順序。3/5/7区間 × epsilon の感度格子で B3−B0 の持続的差の順序は安定せず、現在の保存間隔と持続ルールでは、電流、loadエネルギー、passiveエネルギー、Cl1、MI、W_coh が応答する順序を一意に決められませんでした（M13B、temporal_response_inconclusive）。分岐時刻での電流の即時差は構造的（電流の式に変更した結合係数が含まれる）であり、応答順序の証拠ではありません。
• M12/M13 の観察の背後にある物理機構。結合への感度は観察されています（M13A）が、機構は特定されていません。結合が W_coh 減少の唯一の原因であること、弱い結合が一般に良いこと、最適な結合が存在する/見つかったこと、passiveエネルギーがcoherent ergotropyに変換されたこと、相関の抑制が仕事を守ったこと、相互情報量の減少が W_coh の減少を引き起こした（または変換された）ことは、いずれも示されていません。
• M12/M13 の結果の、この単一の固定matched N=7条件、これら2つの分岐開始時刻、有限窓（分岐 t<=8.50、軌道 t<=10）、この保存解像度を超える一般化。M12/M13 の中では抽出、回収、制御、繰り返しサイクル動作はいずれも行っていません（のちのC2のt=10延長は別の補足チェックで、同じく抽出を行っていません）。
• M14/M15A で使った理想的・瞬間的な理想化を超えるもの。とくに未確認：有限時間のスイッチやramp（モデル化したquenchは数学的に持続時間ゼロで、制御器ハミルトニアンも帯域幅制限もありません）。アクチュエータや制御器の実装コスト（ハミルトニアン跳びの仕事が数値的にゼロであることは、装置コストがゼロという意味ではありません）。抽出ユニタリの実装コスト。switching workと抽出仕事を合わせた正味利益（帳簿は意図的に分けたままで、決して合算しません）。繰り返し抽出やサイクル動作。一般的な長時間安定性、定常状態、t=10 より先の挙動（M13AとM15Aの正式な窓と単一の抽出時刻はいまも t=8.50 のままで、別の補足チェックがC2を t=10 まで伝播させたため、t=9.68 と t=10 での保存格子上の比較は行われています。抽出は伴いません）。パラメータの一般化（他の結合値、掃引、N、gamma、Omega、分岐時刻、ここで固定した以外の条件）。いかなる最適制御（最適値は探索も発見もしておらず、弱い結合が一般に良いとも示されていません）。実ハードウェア上の実験的優位。そして量子的優位性。
• 古典的な波動/確率モデルに対する公平なベースライン。文献に対する新規性/優先権（文献調査は行っていません）。

8. 確認を要する項目（Items requiring confirmation）

2つの元READMEを統合した際、数値や判断の矛盾は見つかりませんでした（Milestone 10は両方に現れ、統合済みです）。次の各点は提供資料からは解決できなかったため、推測せずに未解決のまま残しています。

• スイート全体のCargoテスト数「90/90」。9c validationの実行時チェックは 47/47 PASS、頑健固有値診断の必須チェックは 10/10 PASS ですが、Cargoテストスイート全体について別途言及されている「90/90」は資料から確認できないため、ここでは採用しません。
• Milestone 11g のリリーステスト数。11gのレポートは整形と解析実行のPASSを記録していますが、数値的なリリーステスト数を明記していません。
• 2つの異なる N=7 W_max は矛盾ではありません。3.8081e-3（M9c fixed-total、Omega=0.2）と 3.3853e-3（M11g equal-input、Omega~0.18748）は別の軌道・条件に属します。混同しないでください。誤読を避けるためここに明記します。
• M14/M15A のバイナリ名とリポジトリパス。14A、14 Plan B、15A の各レポートは著者の作業マシン上のソース/出力パスを示していますが、リポジトリ相対のバイナリ名は書かれていないため、Repository layout の bin/ 対応表には推測でエントリを追加していません。
• M14Aの単位チェックCSVと、M15Aの numerical_checks.csv は別の成果物です。14Aは伝播なしの単位チェック11行を記録し、M15Aの numerical_checks.csv は抽出段階の40/40チェックを保持します。同じ一般的な名前のファイルが N=3..7 イベント比較にも存在します。CSVを集める際に3つを混同しないでください。
• 数値的にゼロのswitching workはコストの記述ではありません。W_switch_on_system の 1e-15/1e-16 桁の値は、モデル化したハミルトニアン跳びだけについて1e-12の監査許容差を下回るという意味です。プロジェクトにはアクチュエータ、ramp、帯域幅、装置コストのモデルが存在しないため、これらから制御器コストの結論は導けません。

9. 詳細資料

以下のファイル名は、該当するMilestoneの再現ZIPの中のファイルを指します。一部のレポートは、ブラウザで読めるように reports/ にも置かれている場合があります。各レポートは自分のCSV一式を列挙しています。リポジトリ全体は GitHub にあります。

ここではファイル名のみを示します（個々のファイルへの直接URLは推測していません）。

• 基礎と中心的な比較MILESTONE_4_RESULT.md ／ MILESTONE_5A/5B/5C ／ MILESTONE_6A_REPORT.md ／ Milestone_4-6a_研究結果ノート.pdf（4–6aの要約）
• ノイズ位置と保護MILESTONE_7A/7B/7C/7D_REPORT.md
• 鎖長と実行可能性MILESTONE_8A/8B/8C_REPORT.md
• N=7 完全実行MILESTONE_9A/9B_REPORT.md
• 9c（正本と記録）MILESTONE_9C_VALIDATION.md（最終・正本） ／ MILESTONE_9C_DIAGNOSTIC.md ／ MILESTONE_9C_REPORT.md ／ 主要CSV：n7_fixed_total_validation_summary.csv, ..._checks.csv, ..._trajectory_comparison.csv, fixed_total_noise_final_comparison.csv, robust_eigen_diagnostic_unit_checks.csv
• fixed-total / XGammaMILESTONE_10A/10B/10C_REPORT.md ／ MILESTONE_10_FINAL_REPORT.md
• equal-input N=3 vs N=7MILESTONE_11C_PRECHECK_REPORT.md ／ MILESTONE_11D/11E_REPORT.md ／ MILESTONE_11F/11G/11H_REPORT.md ／ MILESTONE_11I_MOUNTAIN_MECHANISM_DESIGN.md ／ MILESTONE_11J_ASYMMETRIC_TRANSFORM_REPORT.md ／ MILESTONE_11K_POST_PEAK_DECAY_REPORT.md ／ 主要CSV：input_matching_interpolated_trial_summary.csv, equal_input_timeseries_shape_summary.csv, equal_input_peak_widths.csv, equal_input_curve_transform_models.csv, equal_input_asymmetric_transform_models.csv, equal_input_post_peak_decay_models.csv
• matched N=7 dt半分化（11L）MILESTONE_11L_MATCHED_N7_DT_HALVING_REPORT.md
• matched N=7 サイト分解輸送（11M）MILESTONE_11M_SITE_RESOLVED_TRANSPORT_REPORT.md ／ 主要CSV：matched_n7_site_resolved_events.csv
• N=3 Wピーク分岐チェックMILESTONE_N3_W_PEAK_BRANCH_REPORT.md ／ 主要CSV：n3_W_peak_branch_comparison.csv
• equal-input N=3..7 イベント比較MILESTONE_N3_7_EVENT_COMPARISON_REPORT.md（FINAL_REPORT.md として提供） ／ 主要CSV：event_summary.csv, odd_even_temporal_comparison.csv, numerical_checks.csv
• ピーク後のergotropy分解（12A）MILESTONE_12A_LOCAL_ERGOTROPY_DECOMPOSITION_REPORT.md
• 時間的関係（12B）MILESTONE_12B_TEMPORAL_RELATION_REPORT.md ／ M12B_PRECOMMIT.md
• load-鎖相関（12C）MILESTONE_12C_LOAD_CHAIN_CORRELATION_REPORT.md ／ 主要CSV：m12c_interval_changes.csv
• 反実仮想の分岐実験（13A）MILESTONE_13A_FORMAL_SHORT_HORIZON_COUNTERFACTUAL_BRANCH_REPORT.md ／ 主要CSV：m13a_formal_branch_timeseries.csv（正式728行）
• 分岐の時間応答（13B）MILESTONE_13B_WEAKER_COUPLING_TEMPORAL_RESPONSE_REPORT.md ／ M13B_PRECOMMIT.md ／ 主要CSV：m13b_input_audit.csv, m13b_maximum_differences.csv, m13b_interval_summary.csv, m13b_current_event_comparison.csv, m13b_energy_decomposition.csv, m13b_normalized_delta_timeseries.csv, m13b_start_time_comparison.csv
• 結合制御の設計（14A、設計のみ）MILESTONE_14A_COUPLING_CONTROL_AND_SWITCHING_WORK_DESIGN.md
• switching-work会計（14 Plan B）MILESTONE_14_PLAN_B_SWITCHING_ACCOUNTING_REPORT.md ／ 主要CSV：m14_plan_b_switching_work.csv, m14_plan_b_classification.csv, m14_plan_b_scale_diagnostics.csv, m14_plan_b_state_files.csv
• 端点抽出会計（15A）MILESTONE_15A_ENDPOINT_EXTRACTION_ACCOUNTING.md ／ 主要CSV：endpoint_extraction_comparison.csv, gross_work.csv, energy_ledgers.csv, numerical_checks.csv
• C2のt=10延長（補足、新Milestone番号なし）MILESTONE_C2_EXTENSION_TO_T10_COMPARISON.md ／ C2_EXTENSION_INPUT_AUDIT.md ／ 主要CSV：c0_c2_extended_summary.csv, c0_c2_extended_comparison_timeseries.csv, c0_c2_extended_time_area_summary.csv, c0_c2_extended_difference_events.csv, c2_extension_numerical_checks.csv, c2_extension_input_manifest.csv

README.md                 # integrated project overview
LICENSE                   # MIT license
zip/                      # Milestone-specific complete reproduction packages
  *.zip

引用とライセンス

License: MIT（Cargo.toml / LICENSE を参照）。引用の際は、このリポジトリと、使った結果の根拠となる個別の MILESTONE_*.md レポートを参照してください。各レポートが、その数値の正確な条件と留保を定義しています。