ナビエ−ストークス方程式（２）

ナビエ−ストークス方程式（１）からの続き

図：ナビエ−ストークス方程式（ncu_inoue_labから引用）

Ｑ（koshix）：

外力の有無，粘性の有無，圧縮性の有無，その他の条件分類の中で，流体力学の方程式の難易度はどうなっていて，ナビエ−ストークス方程式はどこに位置していて，オイラー方程式その他はどこに対応するのかが，混乱しているのでできるだけ，整理してわかりやすく説明してください。

Ａ（Claude）：

１．整理の枠組み：条件は「2値」ではなく「連続的な階段」

混乱の最大の原因は、外力・粘性・圧縮性を「アリ／ナシ」の2値だと思ってしまうことです。実際には外力の正則性クラスも散逸の強さも連続パラメータで、難易度はその上で滑らかに変化します。整理の軸は7つあります。

｜軸――――――｜パラメータ――――――｜効果――――――――――――――――｜

｜① 空間次元 　｜2D / 3D　　　　　　　｜最大の分岐。渦伸長項の有無　　　　　｜

｜② 散逸の強さ｜$(-\Delta)^\alpha$ の $\alpha$（0＝Euler, 1＝NS）｜大きいほど易しい｜

｜③ 外力　　　｜$f$ の正則性クラス（超臨界→…→$C^\infty$急減少→$f\equiv0$）｜滑らかにするほど難しい｜

｜④ 圧縮性　　｜非圧縮 / 圧縮性　　　　｜困難の種類が変わる　　　　　　　　　｜

｜⑤ 境界　　　｜アリ / ナシ　　　　　　｜境界アリの方が易しい　　　　　　　　｜

｜⑥ 解のクラス｜古典解・強解 / Leray–Hopf弱解 / 一般の弱解｜弱くするほど「易しい（が非一意）」｜

｜⑦ 対称性 　　｜軸対称(swirlなし) など ｜対称性が強いほど易しい　　　　　　　｜

そして、これら全部を貫くマスターキーが「スケーリング臨界性」です。

２．支配原理：超臨界性がすべてを決める

分数階散逸 $(-\Delta)^\alpha$ を持つ一般化NS方程式は、

$$u_\lambda(x,t) = \lambda^{2\alpha-1}u(\lambda x,\ \lambda^{2\alpha}t)$$

で不変です。唯一のア・プリオリ量であるエネルギー $\|u\|_{L^2}^2$ はこの変換で $\lambda^{2(2\alpha-1)-n}$ 倍になる。$n=3$ で指数がゼロになるのは

$$2\alpha - 1 = \tfrac{3}{2} \quad\Longrightarrow\quad \boxed{\alpha = \tfrac{5}{4}}$$

これが臨界点です。これより上なら勝ち、下なら（現在の道具では）負け。

｜$\alpha$―――――――｜名称―――――――――｜状態―――――｜

｜$\alpha \ge 5/4$　　　　｜超散逸（hyperdissipative）｜大域正則性 証明済（Ladyzhenskaya, Lions）｜

｜$5/4$ をわずかに下回る｜対数的超臨界　　　　　｜証明済（Tao 2009, Barbato–Morandin–Romito ほか）｜

｜$\alpha = 1$　　　　　｜標準のNS方程式　　　　｜未解決（ミレニアム問題）｜

｜$0 < \alpha < 1$　　　｜劣散逸（hypodissipative）　｜粗い外力付きなら爆発 証明済｜

｜$\alpha = 0$　　　　　｜Euler方程式　　　　　　｜未解決（ただし爆発が強く予想）｜

ここが決定的に重要な感覚です。標準のNS方程式は、臨界点からわずか $1/4$ だけ超臨界側に落ちている。ほんの少し粘性を強く（非局所的に）してやれば、90年来の難問は既に解けている。その $1/4$ の隙間に、ミレニアム賞100万ドルが埋まっているわけです。

３．軸① 次元：2Dはなぜ易しいのか

3次元の渦度方程式には渦伸長項 $\omega\cdot\nabla u$ があります。

$$\partial_t\omega + (u\cdot\nabla)\omega = \underbrace{(\omega\cdot\nabla)u}_{\text{3Dのみ}} + \nu\Delta\omega$$

2次元では渦度がスカラーになりこの項が消え、$\omega$ は単に流れに沿って運ばれるだけになる。したがって $\|\omega\|_{L^\infty}$ が保存され、爆発しようがない。これが2D Eulerの大域可解性（Yudovich 1963）の本質であり、2D NSも自動的に大域正則です。

渦伸長こそが、非圧縮流体の全困難の単一の源だと言ってよい。逆に言えば、渦伸長を持つ模型方程式はどれも同じ難しさを共有します。

４．軸② 粘性：Eulerの方が難しいのか易しいのか（重要な逆説）

ここが最も混乱しやすい点です。答えは「目標による」。

・爆発を証明したい → Eulerの方が易しい。粘性は特異点を潰そうと働くので、爆発を作るには粘性という敵と戦わねばならない。実際、爆発の証明はほぼすべてEuler側（ないし劣散逸側）で先に達成されています。

・大域正則性を証明したい → NSの方が望みがある。粘性のおかげでエネルギー散逸、Leray弱解の存在、Caffarelli–Kohn–Nirenbergの部分正則性といった構造が使える。Eulerにはこれらが一切ない。

だから「NSはEulerより難しい」も「EulerはNSより難しい」も、どちらも文脈次第で正しい。前回触れた「無限小の摩擦を入れただけで挙動が本質的に変わる」というのは、この非対称性のことです。

なおEuler方程式はクレイの賞金対象外ですが、Fefferman自身が公式記述の中で、これらの問題はEuler（$\nu=0$）についても未解決かつ非常に重要だと明記しています。

５．軸③ 外力：ここが「階段」であることが最大の要点

外力を一切許さないのが本丸（A/B）ですが、外力を許す場合もその正則性で難易度が全く違う。これが今回のニュースの核心なので、階段として並べます。

｜段｜外力 $f$ のクラス―――――――｜難易度――――――｜状況――――――｜

｜０｜超臨界（$u$ から逆算した $f$）｜ほぼ自明｜爆発する $u$ を勝手に決めて $f:=\partial_t u+(u\cdot\nabla)u-\nu\Delta u+\nabla p$ と置くだけ。$f$ も一緒に爆発するので無意味｜

｜１｜$L^1_t C^{1,\epsilon}_x \cap L^\infty_t L^2_x$ 程度 | 難しい｜Córdoba–Martínez-Zoroa（3D Euler, 2023）、Córdoba–Martínez-Zoroa–Zheng（劣散逸NS, ARMA 2026）。well-posedクラス内での爆発＝方程式の力学が生む爆発であり、外力が押し潰しているのではない｜

｜２｜空間的に滑らか（時間方向は粗い）｜さらに難しい | Córdoba–Martínez-Zoroa の IPM（2024）｜

｜３｜時空一様に $C^\infty$ かつ全導関数が急減少（Clayの条件(5)）｜非常に難しい｜Alpöge–Buckmaster–Coiculescu（IPM）→ 2026年8月のBoussinesq・3D Euler → OpenAIのNS主張｜

｜４｜$f \equiv 0$　　　　　　　　　　　｜本丸｜3D Euler：OpenAIが主張（未検証）／3D NS（A/B）：完全に未解決｜

段0と段3の落差を掴むのが肝心です。「外力を許せば簡単」というのは段0の話であって、Clayが要求する段3は、$u$ が無限大に発散する最中も $f$ が最後まで完璧に滑らかで急減少していなければならない。各項が同時に巨大化しながら、外力として残る部分だけが精妙に相殺して滑らかに留まる——この設計が2023年から2026年までの3年をかけて詰められた部分です。

６．軸④ 圧縮性：難しさの「種類」が変わる（同じ物差しに乗らない）

ここも混同されがちですが、圧縮性は難易度の上下というより問題の性格を変える軸です。

｜――――――――――｜特異点の種類―――――――――｜状態――――――――――｜

｜圧縮性Euler（3D）　｜衝撃波（微分の爆発、速度自体は有界）｜有限時間で必ず形成されることが証明済（Sideris 1985、Christodoulouの精密理論）　　｜

｜圧縮性Euler（1D系）｜衝撃波　　　　　　　　　　　　　｜エントロピー解の枠組みで存在・一意性理論が整備（Glimm, Bressan）｜

｜圧縮性NS　　　　　｜真空・密度の退化　　　　　　　　｜大域弱解の存在は証明済（P.-L. Lions, Feireisl）／正則性は未解決／真空を含む古典解は爆発しうる（Xin 1998）｜

｜非圧縮Euler / NS　　｜渦の自己増幅（速度そのものが無限大）｜未解決｜

つまり圧縮性流体では「特異点が起きること」は驚きでも何でもなく、衝撃波として日常的に起きる。しかもそれは弱解（エントロピー解）の理論で数学的に乗り越え済みです。ミレニアム問題が問うている非圧縮の特異点は、衝撃波とは全く別種の、渦度が自己増幅して速度が発散するタイプであり、しかもそれを乗り越える弱解理論が存在しない（後述の非一意性のため）。

もう一つの構造的な違い：非圧縮では圧力が $-\Delta p = \nabla\cdot((u\cdot\nabla)u)$ で非局所的に決まり、情報が無限速度で伝わる。圧縮性では有限伝播速度で双曲型。解析の道具立てが根本的に異なります。

７．軸⑤ 境界：あると易しい

境界の役割はFeffermanが明確に述べています。境界のあるモデルで特異点が見つかっても、それは流体が境界との相互作用のおかげで特異点を作れることを示すにすぎない。境界なしなら、無茶をしているのは流体自身だ、と。

・Hou–Luo 2013：円筒の上下半分を逆回転させるとEulerが爆発する（数値＋誤差評価）。「最初の本格的な特異点の主張」。

・Chen–Hou 2022：それを計算機援用で厳密証明（2D Boussinesq＋境界あり3D軸対称Euler、滑らかな初期データ）。

・Elgindi 2021（Annals）：境界なし $\mathbb{R}^3$ で爆発を証明。ただし初期データは $C^{1,\alpha}$ にとどまり、$C^\infty$ ではない。

つまり3D Eulerの爆発は「境界を許せば滑らかなデータで」「滑らかさを緩めれば境界なしで」既に証明済みで、両方同時（境界なし＋滑らかなデータ＋外力なし）だけが残っていた。そこがOpenAIが主張している地点です。

８．軸⑥ 解のクラス：弱くすると「易しくなるが無意味になる」

これも独立した軸です。

・Scheffer（1993）・Shnirelman（1997）：2D Eulerに、時空コンパクト台を持つ弱解が存在する。静止した流体が外部刺激なしに勝手に動き出し、また止まる。Euler弱解の一意性は劇的に偽。

・Buckmaster–Vicol（2019）：convex integrationにより3D NSの弱解も非一意（有限エネルギーだがLeray–Hopfクラス外）。

・Isett（2018）ほか：Onsager予想の解決。$C^{1/3}$ 以下ではエネルギー散逸する弱解が作れる。

・Vishik（2018）：粗い外力を許すと2D Eulerですら非一意になる。

要するに、弱解の層では「何でも作れる」。だから「NSの弱解を作ってから正則性を示す」という古典的プログラムの後半が機能しない。ミレニアム問題が古典解（$C^\infty$）を要求しているのはこのためです。

９．軸⑦ 対称性と、周辺モデル方程式の難易度階段

爆発現象の研究は、易しいモデルから順に登ってきました。

｜モデル――――――――｜性格――――――――｜状態――――――――｜

｜Constantin–Lax–Majda　　｜1D、非局所項のみ　｜明示的に爆発（解ける）｜

｜De Gregorio / gCLM　　　｜CLM＋移流　　　　｜ほぼ解決｜

｜3D軸対称NS/Euler（swirlなし）｜$\omega_\theta/r$ が輸送される｜大域正則性 証明済｜

｜IPM（多孔質媒体、Darcy則）｜2D、密度輸送　　｜外力付き爆発 証明済｜

｜2D Boussinesq　　　　　｜3D軸対称Eulerの境界近傍と同型｜境界ありで爆発 証明済（Chen–Hou）｜

｜SQG　　　　　　　　　｜$\alpha=1/2$ 臨界散逸なら大域正則（Kiselev–Nazarov–Volberg, Caffarelli–Vasseur）。非粘性は未解決　　　　　　　｜一部解決｜

｜3D Euler　　　　　　　｜本命の一つ　　　　　　｜前述のとおり条件付きで証明済｜

｜3D NS　　　　　　　　｜最終目標　　　　　　　｜A/B 未解決｜

｜平均化NS（Tao 2014）｜エネルギー恒等式は保つ改変NS｜爆発 証明済 ← 「エネルギー法では絶対に決着しない」ことの証明｜

3D軸対称でswirlがない場合が解けているのに、swirlを入れた瞬間に未解決になる——この一点にも渦伸長の効き方が凝縮しています。

１０．全体マップ：何がどこにあるか

最後に、すべてを一枚に並べます。上ほど易しい。

｜設定―――――――――――――――｜状態―――――――――――――――｜

｜2D NS／2D Euler（滑らかデータ、外力なし） ｜✅ 大域正則性 証明済｜

｜3D 超散逸NS $\alpha\ge5/4$（対数的にわずか下も可） ｜✅ 大域正則性 証明済｜

｜3D 軸対称 NS／Euler、swirlなし 　　　　｜✅ 大域正則性 証明済｜

｜3D NS、初期値が小さい／短時間　　　　｜✅ 証明済｜

｜圧縮性Euler（3D）　　　　　　　　　　｜✅ 衝撃波形成 証明済（別種の特異点）｜

｜3D NS Leray弱解の存在　　　　　　　　｜✅ 証明済（1934）／一意性は ❓｜

｜3D NS 部分正則性（特異集合の1次元放物測度＝0）｜✅ CKN 1982。これ以上進むのが極めて困難｜

｜平均化NS／dyadicモデルの爆発　　　　　｜✅ 証明済（エネルギー法の限界の証明）｜

｜3D Euler、境界あり、滑らかデータ、外力なし｜✅ Chen–Hou 2022（計算機援用）｜

｜3D Euler、境界なし、$C^{1,\alpha}$ データ、外力なし｜✅ Elgindi 2021｜

｜3D Euler／劣散逸NS、滑らかデータ、粗い外力｜✅ Córdoba–Martínez-Zoroa 系列（2023–2026）｜

｜IPM・Boussinesq・3D Euler、時空一様に滑らかな外力｜🆕 Buckmaster–Alpöge 2026（Lean検証済、査読前）｜

｜3D NS、滑らかな外力（Clay C/D）　　　　　｜⚠️ OpenAI 2026 が主張（論文・Lean公開、独立検証待ち、Clay未受理） ｜

｜3D Euler、境界なし、滑らかデータ、外力なし｜⚠️ OpenAI が主張（未検証）｜

｜3D NS、外力なし（Clay A/B）　　　　　　　｜❌ 完全に未解決。本丸｜

｜3D NS Leray–Hopf弱解の一意性　　　　　　｜❌ 未解決｜

｜圧縮性NSの正則性　　　　　　　　　　　　｜❌ 未解決｜

１１．まとめると

(1) 難易度の主軸は「超臨界性の度合い」であり、次元・散逸指数・外力の正則性はすべてその一つの物差しの上に並ぶ。

(2) 3D NSは臨界点から $1/4$ だけ外れている。この微小なギャップが90年の壁。

(3) 「外力アリ／ナシ」は2値ではなく5段の階段で、Clayが要求する段3（時空一様に滑らか＋急減少）と段0（逆算した外力）は天と地ほど違う。2026年9月のニュースは、段3に到達したという主張。

(4) 粘性ナシ（Euler）は「爆発を示すには易しく、正則性を示すには難しい」という逆説的な位置にある。

(5) 圧縮性は難易度の軸ではなく、問題の種類を変える軸。圧縮性の特異点（衝撃波）は証明済みかつ理論的に処理済みで、ミレニアム問題とは別世界の話。

(6) 残る最難関は依然として 3D NS、外力なし、境界なし、滑らかデータ、古典解（A/B）。そしてそこへ至る最大の技術的障壁は、外力を消したときに粘性散逸が成長機構を食い潰してしまう点にある。

ナビエ−ストークス方程式（１）

図：これがそれなのか？（ChatGPTによる）

Ｑ（koshix）：

ミレニアム懸賞問題のひとつである，ナビエ−ストークス方程式をAIが解いたという話が聞こえてきました。出来るだけ分かりやすく解説してほしいのです。

(1) 流体力学におけるナビエ−ストークス方程式の位置づけについて

(2) ナビエ−ストークス方程式はどこまでわかってきたか，何が難しいのか

(3) ミレニアム検証問題におけるナビエ−ストークス方程式に関するの正確な条件

(4) AIが何を解いて，それはどこまで画期的なことなのかそうではないの

(5) Buckmaster らを巡る混乱の争点を正確に解きほぐす

(6) 今後どうなるのか

Ａ（ChatGPT + Claude の第一報）：

０．全体像

結論から述べる。2026年9月8日、OpenAIは「ナビエ−ストークス（NS）方程式のミレニアム問題を解いた」と発表した。内部モデルによる証明とLean形式化が公開されている。しかしクレイ数学研究所はまだ受理しておらず、NS問題は現在も「未解決」のままである。加えて、数学とは別に、優先権と情報流出をめぐる論争が並行して起きている。

１． 流体力学におけるナビエ−ストークス方程式

NS方程式は流体に対するニュートンの第二法則である。非圧縮・一定密度なら、

$\partial_t u + (u\cdot\nabla)u = \nu\Delta u-\nabla p+f, \qquad \nabla\cdot u=0$

と書ける。左辺第2項は移流、右辺は粘性拡散、圧力勾配、外力である。

難しさは三つある。第一に、移流項は速度場が自らを運ぶ二次の非線形項である。3次元では渦が引き伸ばされ、エネルギーが小さなスケールへ送り込まれる。第二に、圧力は独立ではなく、

$-\Delta p =\nabla\cdot((u\cdot\nabla)u)$

により速度場全体から非局所的に決まる。第三に、粘性項は平滑化として働くが、その効果と非線形増幅との競争を十分に制御できない。

航空機、気象、血流、乱流など、流体現象の基礎方程式でありながら、3次元では「滑らかな初期状態から始めた解が永久に滑らかなままか」すら分かっていない。ただし実際の流体は分子から成る。数学的特異点が存在しても、それが直ちに物理的な無限大を意味するわけではなく、連続体近似の限界を示すにとどまる可能性がある。

２．どこまで分かっていて、何が難しいのか

2次元では大域正則性が知られている。3次元でも、初期値が十分小さければ大域解が存在し、一般の場合でも短時間なら滑らかな解が存在する。

Lerayは1934年に大域的弱解の存在を示したが、その一意性と滑らかさは未解決である。Caffarelli–Kohn–Nirenbergの部分正則性定理は、適当な弱解の特異集合が非常に小さいことを示すが、特異点そのものは排除しない。オイラー方程式側にはBeale–Kato–Majda条件があり、爆発が起こるなら渦度の最大値の時間積分が発散することが分かっている。

本質的な障害はスケーリングである。NS方程式は

$ u_\lambda(x,t) = \lambda u(\lambda x,\lambda^2t)$

で不変だが、エネルギーはこの変換のもとで小さなスケールを十分に制御しない。特異点が生じるとすれば小スケールであるのに、最も基本的なア・プリオリ評価であるエネルギーがそこを押さえきれない。この「超臨界性」が核心である。

Taoは2014年、非線形項を平均化した改変NS方程式について、エネルギー恒等式を保ったまま「有限時間爆発」を構成した。単純なエネルギー法だけでは大域正則性を証明できないことが、これで鮮明になった。

近年の重要な流れがThomas HouとGuo Luo、そしてDiego CórdobaとLuis Martínez-Zoroaの仕事である。Hou–Luoは境界を持つ3次元オイラーで特異形成の強い証拠を与え、Córdoba–Martínez-Zoroaは「無限カスケード」を組み合わせて特異点を作る解析的手法を発展させた。ただし従来の構成では、合成した外力が十分に滑らかではなかった。残された難所は、特異形成のカスケードを維持したまま、外力を最後まで滑らかに保つことであった。

３．ミレニアム問題の正確な条件

Feffermanの公式記述によれば、3次元NSについて次の四つのうち一つを証明すればよい。

| 条件 | 空間 | 外力 | 主張 |

| -- | -- | -- | -- |

| A | \(\mathbb{R}^3\) | \(f=0\) | 任意の滑らかな初期値から大域的に滑らかな解が存在する |

| B | 周期空間 | \(f=0\) | 同上 |

| C | \(\mathbb{R}^3\) | 滑らかな \(f\) | 条件を満たす滑らかな初期値・外力で、大域的に滑らかな解が存在しない |

| D | 周期空間 | 滑らかな \(f\) | 同上 |

初期値と外力には、すべての導関数について非常に強い滑らかさと減衰の条件が課される。爆発する速度場を勝手に作り、

$f = \partial_tu + (u\cdot\nabla)u -

\nu\Delta u + \nabla p$

と逆算すればよいわけではない。それでは通常、爆発時刻に外力 \(f\) まで発散してしまう。

したがって(C)(D)の真の難しさは、速度場が特異になる一方で、加速・圧力・輸送・粘性の巨大な項が互いに相殺し、残る外力だけは滑らかなままにすることにある。ただし、数学者の多くが本命と考えてきたのは外力なしの(A)(B)である。(C)(D)も公式には正解だが、物理的・解析的にはA/Bのほうが深い問題だと受け取られてきた。

４．AIは何を解いたのか

Buckmaster–Alpögeは9月7日深夜、非圧縮多孔質媒体方程式、ブシネスク方程式、3次元非圧縮オイラー方程式について、滑らかな外力付きの有限時間爆発を公表した。Córdoba–Martínez-Zoroaのプログラムを、粗い外力から滑らかな外力へ押し上げる仕事である。

翌9月8日、OpenAIはNSについて、静止初期状態から滑らかな外力を加え、有限エネルギーのまま有限時間で速度が非有界になる解を構成したと発表した。主張が正しければ、Feffermanの(C)、さらに周期版の(D)を満たす。

OpenAIの説明では、9月1日に「ミレニアム問題が二つ解かれた」という噂を聞いて探索を開始し、多数のエージェントを投入して数日で候補証明に到達、その後Leanで形式化した。計算量はNSだけでメッセージ約270万通、出力約1300億トークンに達したという。

画期的な点は三つある。第一に、公式問題文の一ルートを正面から満たす可能性がある。第二に、Lean形式化が付属している。第三に、人間が長期間かけて進めてきた最先端の構成に、AIエージェント群が数日単位で到達した。

ただし留保も大きい。証明は外力付きであり、外力なしNSへの道は依然として開かれていない。根本的なアイデアの源流はCórdobaとMartínez-Zoroaにあり、Fefferman自身もこの流れの主役は彼らだと評価している。

また、Leanが保証するのは「形式化された命題が形式体系内で正しい」ことだけである。その命題がFeffermanの問題設定と真に一致するかは、人間が確認しなければならない。

Taoが指摘した「答え」と「理解」の乖離も重要である。正しい証明が得られても、なぜその構造が働くのかという数学的理解が同じ速度で得られるとは限らない。AI時代には、この二つが初めて大規模に分離し始めている。

５．Buckmasterをめぐる混乱

争点は少なくとも四つに分けて考える必要がある。

第一は優先権である。外力付きオイラーについてはBuckmaster–Alpögeが先であり、OpenAIもそれを認めている。一方でOpenAIは、自社のオイラー結果は別手法によるものだとし、NSでは類似した方向を用いたことを認めている。

第二はデータ流出である。Buckmasterは未公開の草稿をCodexセッションに入力していたため、それがモデル改善や探索に使われた可能性を疑った。OpenAIは調査の結果、直前2か月の彼のCodexプロンプトは今回のシステムに影響し得なかったと説明している。ただし一般論として、非識別化された利用データがモデル改善に使われる可能性までは完全には否定していない。

第三が最も重要で、なぜOpenAIが「滑らかな外力付き」という極めて特殊なルートに短期間で到達したのかである。Buckmaster側は、この方向はCórdoba–Martínez-Zoroaから続く非常に限定された研究路線であり、問題文を与えただけで自然に選ばれる道ではないと疑問視している。OpenAI側は、公開前の彼らの研究は見ていないと主張している。つまり、直接的な情報漏洩の有無とは別に、「問題が解かれたらしい」という噂そのものが探索空間を劇的に狭めた可能性が争点となっている。

第四は交渉過程である。Buckmasterは、著者扱いや公表方法をめぐって不適切な圧力を受けたと記録している。Bubeck側はその一部を否定しつつ、発言の一部については謝罪している。当事者の説明は一致していない。

したがって現時点で、「盗用だった」「完全に独立だった」のいずれも断定は早い。

６． 今後どうなるか

短期的には、数学者がOpenAIの論文とLean形式化を読み、形式化された定理がFeffermanの(C)(D)と真に一致するかを精査することになる。クレイ数学研究所も即断せず、通常どおり厳密に審査するとしている。公式の決着には査読と一定期間の一般承認が必要であり、時間がかかる。

中期的な焦点は、外力を消せるかどうかである。しかしforcedからunforcedへの移行は自動的ではない。外力なしNSでは、粘性による散逸が特異形成のカスケードを潰す可能性があり、そこが依然として大きな障害である。

物理・工学への直接の影響は当面ほとんどない。乱流モデルや数値流体計算が明日から変わるわけではない。意味が大きいのは、きわめて単純な運動法則が、3次元流体ではこれほど制御困難な振る舞いを許しうるという事実のほうである。

そして最大の変化は、おそらく数学研究の制度面に現れる。著者性、優先権、AI利用の開示、査読、研究データの扱い、そして「AIが得た証明を誰が理解し、誰が責任を持つのか」という問題が、一気に現実のものとなった。

今回の出来事を「AIがミレニアム問題を解いたか」という二値問題だけで見ると、本質を取り逃がす。重要なのは次の三点である。

・公式ルールと数学者の実感の乖離、

・人間が生んだアイデアとAIによる高速探索の関係、

・そして「答えを得ること」と「理解すること」の分離。

今回のナビエ−ストークス騒動は、数学の問題が一つ解けたかもしれないという話である以上に、数学を「解く」とは何を意味するのかが変わり始めた事件である。

［１］ナビエ−ストーク方程式のミレニアム懸賞問題について（OpenAI）

［２］An intuitive exploration of OpenAI’s Navier–Stokes singularity construction（Wolfram Community）

ニュートリノ（２）

ニュートリノ（１）からの続き

図：T2Kによるニュートリノ振動の研究（Geminiによる）

８．人工ニュートリノで振動を精密に測る：K2KからT2Kへ

なぜ人工のニュートリノを使うのか

太陽や大気のニュートリノは自然任せで、量・種類・エネルギー・方向を自分では決められません。加速器で作れば、

種類：ほぼ純粋なミューニュートリノ（または反ミューニュートリノ）

エネルギー：狙った値にそろえられる

方向と距離：発射地点から検出器まで固定

時刻：いつ撃ったかが正確に分かる

と、条件を「実験」として管理できます。

先駆け：K2K（1999〜2004年）

茨城県つくば市の高エネルギー加速器研究機構（KEK）から、約250 km先のスーパーカミオカンデへミューニュートリノを撃ち込みました。到着したミューニュートリノは予想より少なく、大気ニュートリノで見つかった振動を、人工ニュートリノで初めて確認しました。

T2K（Tokai to Kamioka、2010年〜）

茨城県東海村の大強度陽子加速器施設J-PARCから、295 km先のスーパーカミオカンデへ撃ち込む実験です。

ニュートリノビームの作り方

30 GeVの陽子を黒鉛の標的に当てて、パイ中間子を大量に作る

電磁ホーンという磁石でパイ中間子を前方にそろえる

約100 mのトンネルの中でパイ中間子が崩壊し、ミューニュートリノになる

電磁ホーンの電流の向きを逆にすると、逆符号のパイ中間子が集まり、反ミューニュートリノのビームになります。これがあとで大きな意味を持ちます。

ニュートリノは地球を素通りするので、トンネルは要りません。ビームはわずかに下向きに発射され、地下を一直線に通ってスーパーカミオカンデに届きます（経路のいちばん深い所は地下約1.7 km）。

オフアクシスという工夫

ビームの中心軸を、スーパーカミオカンデの方向から2.5°だけずらしています。パイ中間子の崩壊の運動学のため、軸から少しずれた方向では、ニュートリノのエネルギーが約0.6 GeVに鋭くそろいます。これは、295 kmの距離で振動の効果がちょうど最大になるエネルギーです。

前置検出器と後置検出器

発射地点から280 mの所に前置検出器を置き、振動する前のビームの量と性質を測ります。それを295 km先のスーパーカミオカンデで測った結果と比べ、途中で何が変わったかを引き算で取り出します。

見ているもの：これまでと同じ「リング」

スーパーカミオカンデで見るのは、梶田さんの大気ニュートリノのときと同じです。

ミューニュートリノ → ミューオン → くっきりしたリング

電子ニュートリノ → 電子 → ぼやけたリング

大量の太陽ニュートリノや大気ニュートリノに埋もれないのかという、前に出た問いがここでも問題になりますが、T2Kでは時間という物差しが決定的に効きます。ビームは数マイクロ秒のパルスで撃たれ、発射時刻はGPSで同期されています。その時間窓に入った事象だけを数えれば、偶然紛れ込む太陽・大気ニュートリノはほぼゼロになります。それでも届くビーム事象は、数年かけて数百個程度です。

何が発見・検証されたのか

① 電子ニュートリノの「出現」（2011年に兆候、2013年に確立）

スーパーカミオカンデで見えたのは「ミューニュートリノが減る」という消失でした。T2Kは、ミューニュートリノのビームの中に、本来ないはずの電子ニュートリノが現れることを捉えました。2011年は6事象で兆候の段階でしたが（同年3月の東日本大震災で施設が被災し、一時中断）、2013年に統計的に確立しました。

これは、ミューから電子への振動を支配する混合角（θ₁₃）がゼロでないことを意味します。同じ時期に、中国の大亜湾実験などが原子炉からの反電子ニュートリノの減少から、このθ₁₃を精密に測定しました。

② ミューニュートリノ消失の精密測定

大気ニュートリノで見つかった振動（主にミューからタウへ）の大きさと周期を、制御されたビームで高精度に測りました。結果は、ミューとタウがほぼ最大限に混ざり合っていることを示しています。

③ CP対称性の破れの探索（2016年〜）

θ₁₃がゼロでないと分かったことで、次の大問題に挑めるようになりました。ニュートリノと反ニュートリノで振動の仕方が違うか、という問題です。

・ミューニュートリノ → 電子ニュートリノ

・反ミューニュートリノ → 反電子ニュートリノ

この二つの起こりやすさを比べます。T2Kの結果では、ニュートリノ側では予想より多く、反ニュートリノ側では少なく現れました。2020年には、CP対称性が保たれている可能性を約95%の信頼度で退け、破れが大きい値が好ましいと報告しています（ただし「発見」と言える統計的確実性にはまだ届いていません）。

これがなぜ重要かというと、宇宙にはなぜ反物質がほとんどなく、物質ばかりが残ったのかという謎があり、ニュートリノでの物質・反物質の非対称が、その答えの鍵になる可能性があるからです（レプトジェネシス説）。

④ 他実験との統合

その後、T2Kとスーパーカミオカンデの大気ニュートリノの合同解析が行われました。さらに、ライバルとされてきた米国のNOvA実験とT2Kが、2025年10月に初めての合同解析を発表しています。距離（NOvAは約810 km）やエネルギーの異なる実験を組み合わせることで、3種類のニュートリノの質量の並び順（質量順序）などの未解決問題に迫っています。

次の段階：ハイパーカミオカンデ

スーパーカミオカンデの後継機で、タンクの容積は26万トン、有効体積は約19万トン（スーパーカミオカンデは2.25万トン）です。2025年8月に巨大空洞の掘削が完了し、運転開始は2028年の予定とされています。J-PARCのビームも強化し、T2Kの延長線上で、CP対称性の破れを確実に「発見」することが最大の目標です。陽子崩壊探索や、銀河系内で超新星が起きたときの大量のニュートリノ観測も担います。

年表への追加

　年――――｜出来事―――――――――――｜何が分かったか――――――――――――

1999〜2004 ｜K2K（つくば→神岡 250 km）　｜人工ニュートリノで振動を初確認

2010 　　｜T2K開始（東海→神岡 295 km）

2013 　　｜T2K　　　　　　　　　　　　｜ミュー→電子ニュートリノの出現を確立（θ₁₃≠0）

2012〜 　　｜原子炉実験（大亜湾など） 　｜θ₁₃を精密測定

2016〜20 ｜T2K（反ニュートリノとの比較）｜ CP対称性の破れの兆候

2025 　　｜T2KとNOvAの合同解析 　　　｜質量順序・CP破れへの制約

2028予定 ｜ハイパーカミオカンデ稼働 　｜CP破れの発見を目指す

一言でつなぐと、自然のニュートリノで「振動がある」ことが発見され（1998年）、加速器で作った人工のニュートリノで「どう振動するか」を測る時代に入り、いまや「ニュートリノと反ニュートリノで振動が違うか」という宇宙の物質の起源につながる問いに挑んでいる、という流れになります。

自然界の左右対称性

ホモキラリティ（１）からの続き

図：宇宙から生命へキラリティの進化（ChatGPTによる）

Ｑ（koshix）：

ノーベル賞で，自然界の左右対称性に絡んだ研究は，物理学賞，化学賞でたくさんありますが，すべてを年代順にリストアップしてください。

Ａ（ChatGPT）：

ノーベル賞に見る「自然界の左右対称性」研究史

結論から言うと、自然界の左右対称性をめぐるノーベル賞には、大きく二つの系譜があります。一つは物理学のパリティ（空間反転）対称性とその破れ、もう一つは化学の分子キラリティ（鏡像異性）と不斉合成です。両者をつなぐのが、「自然界はなぜ左と右を区別するのか」という問いです。

ただし、「関連する研究をすべて」とすると、核物理学の対称性、結晶構造、生体分子の立体構造、さらには一般的な対称性の自発的破れまで対象が広がります。

そこで以下では、関連度を三段階に分類します。

◎ 直接関連：左右非対称性、鏡像異性、パリティの破れが研究の中心。

 ○  強い関連：左右対称性が研究の重要な要素。

△ 周辺関連：より一般的な対称性や立体構造の研究。左右性そのものが受賞理由ではない。

特に2026年10月7日に発表された化学賞は、この歴史の最新の重要な到達点です。

１．年代順一覧

受賞年・賞	受賞者と研究内容	関連度
1902 化学	エミール・フィッシャー：糖類などの構造研究。糖の立体化学を体系化	○
1913 化学	アルフレート・ヴェルナー：配位化学の確立。金属錯体の光学異性を実証	◎
1937 化学	ウォルター・ハース：糖類とビタミンCの構造。糖の立体配置解明に貢献	○
1954 化学	ライナス・ポーリング：化学結合と複雑な分子構造。タンパク質のらせん構造研究にも貢献	○
1957 物理	楊振寧、李政道：弱い相互作用におけるパリティ非保存を提唱	◎
1962 化学	マックス・ペルーツ、ジョン・ケンドルー：球状タンパク質の立体構造解明	○
1963 物理	ユージン・ウィグナー：原子核・素粒子物理学における基本的対称性原理	○
1964 化学	ドロシー・ホジキン：重要な生体分子のX線結晶構造解析	○
1965 化学	ロバート・ウッドワード：複雑な有機化合物の合成。立体化学を精密に制御	○
1969 化学	デレク・バートン、オッド・ハッセル：分子の立体配座とその化学的意味	○
1975 化学	ジョン・コーンフォース、ウラジミル・プレローグ：酵素反応と有機分子の立体化学	◎
1979 物理	シェルドン・グラショウ、アブドゥス・サラム、スティーヴン・ワインバーグ：電弱統一理論。弱い相互作用のカイラル構造を説明	◎
1980 物理	ジェームズ・クローニン、ヴァル・フィッチ：中性K中間子のCP対称性の破れを発見	◎
1981 化学	福井謙一、ロアルド・ホフマン：化学反応における軌道対称性。鏡映対称性とは異なるが関連する対称性研究	△
1988 物理	レオン・レーダーマン、メルヴィン・シュワルツ、ジャック・シュタインバーガー：ミューニュートリノ発見、弱い相互作用のレプトン構造	○
1999 物理	ヘーラルト・トホーフト、マルティヌス・フェルトマン：電弱理論の量子構造の解明	○
2001 化学	ウィリアム・ノールズ、野依良治、バリー・シャープレス：触媒的不斉合成	◎
2008 物理	南部陽一郎、小林誠、益川敏英：自発的対称性の破れとCP対称性の破れの理論	◎
2021 化学	ベンジャミン・リスト、デイヴィッド・マクミラン：不斉有機触媒の開発	◎
2026 化学	アンリ・カガン、硤合憲三：不斉有機合成における非線形効果と自己触媒作用の発見	◎

この一覧は、ノーベル賞の公式受賞理由に加え、関連する研究内容も考慮したものです。したがって、◎についても受賞理由と左右対称性の関係が同程度に直接的というわけではありません。特に1979年の電弱統一理論、2008年の南部の研究は、対称性の破れの種類を区別する必要があります。

２．とくに重要な研究の流れ

1957年・物理学賞：左右対称性の破れ

李政道と楊振寧は、弱い相互作用ではパリティが保存されない可能性を指摘しました。呉健雄らによるコバルト60のベータ崩壊実験で実証され、左右対称性が自然界の普遍的法則ではないことが明らかになりました。

1975年・化学賞：生命反応の立体選択性

コーンフォースは酵素が立体的に区別される原子や反応経路を選択する仕組みを、プレローグは有機分子と反応の立体化学を解明しました。生体分子が左右を区別することを理解するための基礎です。

1980年・物理学賞：CP対称性の破れ

パリティ反転Pと粒子・反粒子の交換Cを組み合わせても、自然界は完全には対称でないことを発見しました。時間反転対称性や宇宙の物質・反物質の非対称性を考える重要な出発点です。

2001年・化学賞：人工的に左右を作り分ける

ノールズ、野依、シャープレスは、キラルな触媒を使って一方の鏡像異性体を優先的に合成する方法を開発しました。特に医薬品製造で大きな意味を持つ成果です。

2026年・化学賞：微小な左右差の増幅

カガンは不斉合成における非線形効果を、硤合は不斉自己触媒反応を研究しました。ごく小さな鏡像異性体の偏りが、反応によって大きく増幅され得ることを示した点が重要です。これは生命のホモキラリティの起源に対する有力な化学的手掛かりですが、生命の起源そのものが解明されたわけではありません。

３．物理学と化学の関係についての重要な注意点

物理学のパリティ非保存と、生命におけるホモキラリティは、どちらも「左と右の非対称性」ですが、両者の因果関係は確立していません。

弱い相互作用は本質的に左右を区別します。一方、不斉自己触媒反応では、極めて小さい初期の偏りが、化学反応によって大きく増幅されます。その初期の偏りが弱い相互作用によって生じたのか、偶然の揺らぎによるのか、別の物理化学的過程によるのかは、なお研究課題です。

私が特に重要だと考えるのは、1957年の「自然は左右を区別する」から、2026年の「微小な左右差が化学反応で増幅される」へ至る約70年間の展開です。ただし、これは一直線につながった一つの研究史ではなく、物理学と化学で独立に発展した二つの研究史を比較したものです。

ニュートリノ（１）

図：ニュートリノ探究の歴史（Geminiによる）

アイスキューブ・ニュートリノ観測所が，2026年のノーベル物理学賞をもらったので，書きためていた予定記事を前倒しで公開することに。

発端は，NHKのフロンティア「ニュートリノが開く “深宇宙”の扉」である。有名な石原安野さんが登場していた。うかつにも安野は「あや」と読むことを本人の口で始めて知った。

カミオカンデだけではなく，アイスキューブの現場の映像がみられてとても興味深かった。そこで，ニュートリノが手のひらに6兆個／秒も降り注いでいるところに驚いてClaudeにいろいろと質問した結果をまとめたのが以下の記事である。

Ｑ（koshix）：

ここまでのストーリーをまとめてください。どんな現象で，何が生じていて，何が発見され／検証されたのかを歴史的にまとめて，上記の質問や疑問に答えながら分かりやすく説明してほしいです。

Ａ（Claude）：

０．主役：ニュートリノとは

ニュートリノは、電荷を持たず、他の物質とほとんど反応しない素粒子です。1930年にパウリが、ベータ崩壊でエネルギーの勘定が合わないことを説明するために存在を予言し、1956年に原子炉からのニュートリノとして初めて検出されました。電子ニュートリノ、ミューニュートリノ、タウニュートリノの3種類（フレーバー）があります。素粒子の標準模型では、長らく質量ゼロと考えられていました。

１．私たちはニュートリノの雨の中にいる

最大の発生源は太陽です。太陽の中心で水素がヘリウムに変わる核融合のたびに、電子ニュートリノが生まれます（太陽ニュートリノ）。

・地球に届く量：1 cm²あたり毎秒約650億個

・手のひら10cm四方なら、毎秒約6〜7兆個

ほとんど反応しないので、昼は空から、夜は地球を貫いて足元から、同じように体を素通りしていきます。人体で実際に反応するのは、一生に1回あるかないかです。

このほかに、地球内部の放射性元素から出るもの、宇宙線が作る大気ニュートリノ、ビッグバンの名残などがありますが、数は太陽のものよりずっと少数です。

２．最初の謎：太陽ニュートリノが足りない（1960年代末〜）

デービスは、地下の金鉱に塩素を含む約600トンの液体を置き、電子ニュートリノが塩素をアルゴンに変える反応を利用して、生じたアルゴン原子を数えました。結果は理論予想の約3分の1でした。

太陽のモデルが間違っているのか、実験が間違っているのか、それともニュートリノ自体に何かあるのか。これが30年以上続く「太陽ニュートリノ問題」の始まりです。

３．カミオカンデ：ニュートリノを「見る」装置（1983年〜）

カミオカンデは岐阜県神岡鉱山の地下1000 mに作られた、3000トンの水タンクです。

何を見ているのか

ニュートリノそのものは見えません。ニュートリノが水中の粒子と反応して生まれたミューオンや電子が光速に近い速さで走るときに出す、青白いチェレンコフ光を、壁面の光センサーで捉えます。光はリングとして映り、その位置と形から粒子の方向とエネルギーが分かります。

数の多い太陽ニュートリノに、他の信号はかき消されないのか？

かき消されません。三つの物差しで区別できるからです。

(1) エネルギー：太陽ニュートリノの大半は0.4 MeV以下と低エネルギーで、約7 MeV以上しか見えない検出器にはそもそも映りません。見えるのは全体の約1万分の1だけです。しかも反応確率はエネルギーが高いほど大きいので、GeV級の大気ニュートリノと実際の事象数は同程度になり、エネルギー域も重なりません。

(2) 方向：太陽ニュートリノに弾かれた電子は、太陽と反対の方向へ飛びます。

(3) 時間：定常的に降る太陽ニュートリノに対して、超新星のニュートリノは短時間に集中して来ます。

低エネルギー側で本当に邪魔になるのは、太陽ニュートリノではなく、水や岩盤に含まれるラドンなどの放射性物質です。

二つの成果

・1987年 超新星1987A：16万光年先の超新星爆発から、13秒間に11個のニュートリノを捉えました。ニュートリノ天文学の誕生で、小柴昌俊さんの2002年ノーベル賞の理由になりました。

・1989年 太陽ニュートリノ：事象が太陽の方向から来ていることを示し、本当に太陽から来ていることをリアルタイムで確かめました。数はやはり予想の約半分で、デービスの不足が実験の誤りではないことも裏付けました。

４．第二の謎：大気ニュートリノの異常（1980年代後半）

大気ニュートリノとは

宇宙から来た陽子などの宇宙線が上空15〜20 kmで大気の原子核に衝突すると、パイ中間子などが生まれ、その崩壊でニュートリノができます。崩壊の連鎖から、ミューニュートリノと電子ニュートリノがおよそ2対1でできると予想されます。エネルギーはGeV級と、太陽ニュートリノよりずっと高くなります。

宇宙線は太陽から来るのではないか？

違います。大気ニュートリノのもとになる宇宙線の大半は、銀河系内の超新星残骸などで加速された銀河宇宙線です。電荷を持つため、銀河の磁場で何百万年もかき混ぜられ、地球にはほぼ全方向から一様に降ってきます。むしろ太陽は宇宙線を遮っていて、太陽の方向にはわずかに少ない「太陽の影」が観測されています。また、太陽風は外からの宇宙線を押し返すので、太陽活動の盛んな時期には宇宙線が減ります。

異常の発見

梶田隆章さんは、カミオカンデの大気ニュートリノの事象を解析し、リングの形で2種類を見分けました。

・ミューオン（ミューニュートリノ由来）→ くっきりしたリング

・電子（電子ニュートリノ由来）→ ぼやけたリング

その結果、ミューニュートリノの割合が予想の2対1より明らかに小さいことが分かりました。ミューニュートリノがどこかで消えているように見えるこの現象が、「大気ニュートリノ異常」です。

５．決着①：スーパーカミオカンデ（1998年）

5万トンの水を使う後継機で、梶田さんたちは飛来方向別にニュートリノを数えました。

ここで「全方向から一様に来る」ことが決定的に効いてきます。

・真上から来るもの：飛行距離は約15 km

・真下から来るもの：地球を貫いて約1万3000 km

もとの量は上下でほぼ対称なはずなのに、結果は次のとおりでした。

｜種類―――――――｜上から（短距離）―｜下から（長距離）｜

｜ミューニュートリノ｜ほぼ予想どおり　　｜約半分に減少　　｜

｜電子ニュートリノ　｜予想どおり　　　　｜予想どおり　　　｜

長い距離を飛ぶと減るというパターンは、ミューニュートリノが飛びながら別の種類に変わる「ニュートリノ振動」の特徴そのものです。電子ニュートリノは増えていないので、変わった先は主にタウニュートリノです。このタウの出現は、のちに欧州のOPERA実験などで確認されました。

なぜ「質量がある」ことになるのか

種類が移り変わるには、3種類がそれぞれ異なる質量を持っている必要があります。質量がすべてゼロなら振動は起きません。この発見は、標準模型を超える物理の最初の確かな証拠になりました。

６．決着②：SNO（2001〜2002年）と太陽ニュートリノ問題の解決

カナダのSNOは重水を使い、

・電子ニュートリノだけに反応する反応

・３種類すべてに等しく反応する反応

を同時に測りました。結果は、電子ニュートリノだけなら予想の約3分の1、3種類の合計なら予想どおりでした。

つまり、太陽のモデルは正しく、太陽で作られた電子ニュートリノが途中でミューやタウに変わっていたのです。30年来の太陽ニュートリノ問題も、原因は振動でした。

整理すると

・大気ニュートリノ：減ったのはミューニュートリノ（主にタウへ）

・太陽ニュートリノ：減ったのは電子ニュートリノ（ミューとタウへ）

別々の現象から同じ結論が出たことで、2015年のノーベル物理学賞が、梶田さんとSNOのマクドナルドさんに贈られました。授賞理由は「ニュートリノに質量があることを示した、ニュートリノ振動の発見」です。

７．次の舞台：アイスキューブと宇宙ニュートリノ（2010年〜）

アイスキューブは，南極点の氷の地下1450〜2450 mに、5000個余りの光センサーを約1 km³にわたって埋めた巨大検出器です。狙うのは、銀河系外の天体から来るTeV〜PeV級の宇宙ニュートリノです。

何を見ているのか

原理はカミオカンデと同じで、反応で生まれた荷電粒子のチェレンコフ光を見ます。高エネルギーのニュートリノが氷の原子核の中のクォークと反応すると、核子は粉々になってシャワーになります。弾き出された陽子が光るわけではありません。事象は主に二つの形で現れます。

・トラック型：ミューニュートリノから生まれたミューオンが数kmも走り、光の円錐を出しながら進みます。方向精度は約1°で、天体の特定の主力です。

・カスケード型：電子ニュートリノや中性カレント反応によるシャワーが、ほぼ球状の光の塊になります。エネルギーはよく測れますが、方向は約10°とぼやけます。

ここでの敵は大気ニュートリノと大気ミューオン

太陽ニュートリノはエネルギーが低すぎて、まったく問題になりません。代わりに、

・大気ミューオン（宇宙線が作ったミューオンが直接届いたもの）が毎秒約3000個 → 地球をフィルターにして、下から上へ向かう事象だけを選びます。

・大気ニュートリノが年に約10万個、宇宙ニュートリノは年に数十個 → エネルギー（100 TeV以上では宇宙由来が多数派）、種類の比率（大気はほぼミュー、宇宙由来は振動でほぼ1:1:1）、方向の集中で見分けます。

ブレーザーからのニュートリノも、届くときにはほぼ1:1:1なので、どちらの型にもなります。ただし天体の同定には、方向精度の高いトラック型が主に使われます。

成果

・2013年：約1 PeVの事象から、宇宙ニュートリノの存在を確立

・2017年：トラック型事象IceCube-170922A（約290 TeV）の方向にあったブレーザーTXS 0506+056がガンマ線で増光していることが分かり、初めて発生源の候補を特定

・2022年：セイファート銀河NGC 1068から

・2023年：天の川銀河の面から（こちらはカスケード型事象を使用）

 全体の流れ

｜ 年― ―｜出来事 ―――――――――――――｜ 何が分かったか―――――――――｜

｜ 1930　｜パウリの予言　　　　　　　　　　 ｜ 見えない粒子の存在を仮定　　　　｜

｜ 1956　｜原子炉で初検出　　　　　　　　　 ｜ ニュートリノは実在する　　　　　｜

｜ 1968〜 ｜デービスの塩素実験　　　　　　　｜ 太陽ニュートリノが約1/3しかない（謎の始まり） ｜

｜ 1987　｜カミオカンデが超新星1987Aを観測 ｜ ニュートリノ天文学の誕生　　　　｜

｜ 1989　｜カミオカンデが太陽の方向から検出｜ 太陽起源を確認、不足も裏付け　　｜

｜ 1980年代後半 ｜大気ニュートリノ異常　　　｜ ミューニュートリノが足りない　　｜

｜ 1998　｜スーパーカミオカンデ　　　　　　｜ 長距離でミューニュートリノが減る → 振動、質量あり｜

｜ 2001〜02｜SNO　　　　　　　　　　　　　｜太陽の電子ニュートリノが他の種類に変化、合計は予想どおり ｜

｜ 2002　｜ ノーベル賞（デービス、小柴）　　｜宇宙ニュートリノの検出　　　　　｜

｜ 2013〜23｜アイスキューブ　　　　　　　　 ｜銀河系外・銀河系内からの高エネルギーニュートリノ｜

｜ 2015　｜ ノーベル賞（梶田、マクドナルド）｜ニュートリノ振動の発見　　　　　｜

｜ 2026　｜ ノーベル賞（ハルツェン）　　　　｜高エネルギーニュートリノの発見　｜

一言でいえば、「数が合わない」という小さな食い違いを手がかりに、ニュートリノに質量があることが突き止められ、その後ニュートリノは宇宙を見る新しい「目」になった、という物語です。

ヴィルヘルム・レンツ

 イジング模型からの続き

写真：ヴィルヘルム・レンツ（ロストック大学から引用）

イジング模型を調べてみると，これを提案したのはアルノルト・ゾンマーフェルト（1868-1951）の弟子だったヴィルヘルム・レンツ（1888-1957）であり，1次元の場合にこれを具体化して博士論文にしたのが，レンツの学生のエルンスト・イジング（1900-1998）だということだった。

物理でレンツといえば，レンツの法則が思い浮かぶわけだ。そちらの方はドイツ系ロシア人の物理学者のハインリヒ・レンツ（1804-1865）によるものだ。ヴィルヘルム・レンツの方は，イジング模型に加えて，水素原子の量子力学におけるルンゲ＝レンツベクトルでも有名だ。

混乱しやすいレンツだが，2026年4月24日段階で，Wikipedia日本語版にヴィルヘルム・レンツの項目がなかった。英語版を翻訳して日本語版をつくることにチャレンジしてみた。Wikipediaのアカウントは取得済みだ。言語欄から新しい言語として日本語への翻訳を選択する。その後は，段落ごとに翻訳して（DeepLやClaudeの助けを借りて）リンクを日本語記事に植え替えていく。いろいろ警告はでてくるもののなんとかできあがった。

イジング模型

図：1次元イジング模型（ChatGPTによる）

おなじみのイジング模型は格子点上にスピンが配置されたモデルであり，磁性体の相転移などの考察に用いられる。1920年に，ヴィルヘルム・レンツが提案し，その学生のエルンストイジングが1924年に博士論文で1次元イジング模型を研究した。

Ｑ（koshix）：

１次元イジングモデルの導出について，できるだけわかりやすく，かつ，簡潔にお願いします。LaTeXで下さい。

Ａ（ChatGPT）：

古典論的1次元イジングモデルでは，各格子点にスピン変数，$s_i=\pm 1$を置く。

隣り合うスピンだけが相互作用するとし，外部磁場を\(h\)，相互作用定数を \(J\) とすると，ハミルトニアンは$H=-J\sum_{i=1}^{N}s_i s_{i+1}-h\sum_{i=1}^{N}s_i$である。

ここで周期境界条件$s_{N+1}=s_1$を仮定する。

温度 \(T\) における分配関数は，$Z=\sum_{\{s_i\}}\exp(-\beta H)$である。ただし$\beta=\frac{1}{k_{\mathrm B}T}$。

ハミルトニアンを代入すると，$Z=\sum_{\{s_i\}}\exp\left[\beta J\sum_{i=1}^{N}s_i s_{i+1}+\beta h\sum_{i=1}^{N}s_i\right]$となる。

ここで磁場の項を隣接する結合に半分ずつ分けると，$\beta h \sum_{i=1}^{N}s_i=\sum_{i=1}^{N}\frac{\beta h}{2}(s_i+s_{i+1})$と書ける。

したがって分配関数は，$Z=\sum_{\{s_i\}}\prod_{i=1}^{N}\exp\left[\beta J s_i s_{i+1}+\frac{\beta h}{2}(s_i+s_{i+1})\right]$となる。

ここで転送行列 \(T\) を$T_{s,s'} = \exp \left[ \beta J s s' + \frac{\beta h}{2} (s+s') \right]$ と定義する。スピン \(s,s'\) はそれぞれ \(\pm 1\) をとるので，$T= \begin{pmatrix} e^{\beta J+\beta h} & e^{-\beta J} \\ e^{-\beta J} & e^{\beta J-\beta h} \end{pmatrix} $ である。

このとき分配関数は，転送行列の積のトレースとして$ Z=\operatorname{Tr}T^N $ と書ける。転送行列の固有値を \(\lambda_+\), \(\lambda_-\) とすると，$ Z=\lambda_+^N+\lambda_-^N $ である。

転送行列の固有値は $ \lambda_{\pm} = e^{\beta J}\cosh(\beta h) \pm \sqrt{ e^{2\beta J}\sinh^2(\beta h) + e^{-2\beta J}} $ となる。

熱力学極限 $ N\to\infty $ では，大きい固有値 \(\lambda_+\) だけが支配的なので，$ Z\simeq \lambda_+^N $となる。

したがって，1スピンあたりの自由エネルギーは$ f = -\frac{1}{\beta} \lim_{N\to\infty}\frac{1}{N}\log Z = -\frac{1}{\beta}\log\lambda_+$ である。

特に外部磁場がない場合，$ h=0 $ ならば，$ \lambda_+ = e^{\beta J}+e^{-\beta J} =2\cosh(\beta J) $なので，$ f =-\frac{1}{\beta} \log \left[ 2 \cosh(\beta J) \right] $となる。

この結果から，1次元イジングモデルでは有限温度で相転移がないことも分かる。
なぜなら，自由エネルギーは $ T>0 $ で解析的であり，特異性を持たないからである。つまり，1次元イジングモデルは低温ではスピンがそろいやすいが，有限温度では真の強磁性相転移を起こさない。

定数がわからない

カノニカル分布からの続き

図：エントロピーがわからない（ChatGPTによる）

４年越しに自分（koshix）の書いた記事と過去の記事を読み比べたClaudeが，私（koshix）の文体でブログ記事を書いた。かつ，当時からの疑問点を解決したらしい。言葉の上ではわかったように見えるが，自分の腑に落ちたわけではない。それでいいのか。生成AIは知ったかぶりを増殖させるが，本当に自分が納得する時間はそこからは得られない。

Ａ（Claude）：

しばらく前に四年三ヶ月ぶりに統計力学の記事を書いた。前作としてミクロカノニカル分布にリンクを張ったが，実際に途切れていたのはその十数日後，2022年5月11日の理想気体のエントロピーである。その記事の最後の一文はこうだった。理想気体の熱力学的エントロピーと統計力学的エントロピーは，ともに $N k_B \log T^{3/2}V + const.$ の形をしているが，定数部分まで含めて同じかどうかがよく理解できていない。

そこで止まっていたのだった。

１．そこまでの道筋

始まりは2021年11月である。$N$ 個の粒子を $M$ 個の箱に配り，配置の場合の数 $W$ の対数を最大にする。制約はふたつ，粒子数一定とエネルギー一定で，これをラグランジュの未定乗数 $\alpha,\ \beta$ で処理する。この手続きを三度繰り返した。箱の中の $g_i$ 個の状態が区別できて粒子が何個でも入れるならボルツマン分布，区別できず何個でも入れるならボース分布，区別できず一個までならフェルミ分布。取り替えたのは $W_i$ の数え方だけで，あとは同じ計算である。得られるのは状態の占有率 $f_i$ であって，未定乗数はあとから $\alpha=-\mu/k_BT,\ \beta = 1/k_BT$ と同定した。同定の道具は $S = k_B \log W$ と熱力学の関係式である。ボルツマンの原理は天下りに使った。

翌2022年の4月から5月にかけては，授業のための復習が続いた。ミクロカノニカル分布で等重率の原理を確認し，カルノーサイクルの計算を練習し，モル比熱の定義が教科書ごとにそろっていないことに苛立った。そして5月6日，エントロピーそのもので立ち止まった。天下りでボルツマンの原理を導入してしまえばあとは計算だけになるが，それでは熱力学との関係がうやむやになる。だから熱力学の側から入りたい。$T$ が示強変数なら，それと相補的でエネルギーの次元を与える示量変数として $S$ を置き，$d'Q = TdS$ とする。カルノーサイクルの断熱過程で計算してみれば，たしかに経路によらない。ただしここまでは準静的過程の話であって，不可逆過程と第二法則との関係がつかない。

その関係をつけるために摩擦を持ち込んだ。教科書はクラウジウスやトムソンの原理から抽象的な熱機関の組み合わせで議論を進めていくが，具体的な構成をしていない。ならば，等温過程でだけピストンに摩擦が働くカルノーサイクルを作ればよい。5月8日に効率が $\eta' \le \eta_c$ となることまでは出たものの，$\int d'Q/T \le 0$ につながらない。詰まったのは二日である。5月10日，摩擦力による仕事 $\delta w$ として書く代わりに，摩擦で生じた熱 $\delta q = -\delta w$ として書き直したら通った。

$ \frac{Q'_{\rm H}}{T_{\rm H}} + \frac{Q'_{\rm L}}{T_{\rm L}} $

$= \frac{\delta q_{\rm H}}{T_{\rm H}}+\frac{\delta q_{\rm L}}{T_{\rm L}} \le 0 \quad \Longrightarrow \quad \oint \frac{d'Q}{T_{\rm ex}} \le 0$

同じ量を仕事の側から書くか熱の側から書くかという，それだけの違いだった。

そして翌日の5月11日である。理想気体のエントロピーを二本の道で計算した。熱力学からは $dS = dU/T + pdV/T$ に状態方程式と $U = \frac{3}{2}Nk_BT$ を代入して積分し，統計力学からは $6N$ 次元位相空間の体積を $h^{3N}$ と $N!$ で割った状態数に $S = k_B \log W$ を適用する。前者は $S_0 + Nk_B \log T^{3/2}V$，後者はザックール-テトローデの式になる。$\log$ の中身が無次元で $V/N$ を含むから示量変数になっていることは確かめた。しかし定数が一致するのかどうかがわからない。

２．四年後の答え合わせ

先日この経緯をClaudeに読ませたところ，定数の問いには答えがあると言う。要点は，二つの定数が同じかどうかを問う前に，両者の身分が違うということだった。

熱力学が与えるのは $dS = d'Q/T$ の積分であって，出てくるのはつねに差である。基準点 $S_0$ は積分定数として残り，熱力学の枠内では決まらない。決まらないのは計算が足りないからではなく，理論がそこまでしか言わないからである。一方の統計力学は絶対値を与える。ただしそれは，$h^{3N}$ で割り $N!$ で割るという二つの操作を経てのことだ。

この二つはどちらも熱力学に存在しない。$N!$ は同種粒子が区別できないという要請から来ていて，これを落とすとエントロピーが示量的でなくなる——ギブズのパラドックスである。四年前に $V/N$ の形を見て示量性を確認していたのは，まさにこの $N!$ の効き目を見ていたことになる。そして $h^{3N}$ のほうは，位相空間の体積を無次元の状態数に変えるために要る。$\log$ の中身が無次元であってほしいという，あのとき当たり前のように書いた条件が，プランク定数を呼び込んでいたのだった。

つまり答えはこうなる。両者の定数が一致するかどうかは，熱力学の内側では問えない。統計力学が与える定数は，熱力学の外にある二つの量を借りて作られているからである。逆に言えば，熱力学がエントロピーの絶対値を語れるようになるためには，公理をひとつ足す必要がある。それが第三法則で，ネルンストが化学平衡の定数から辿り着き，プランクが $S_0 = 0$ という形に整えたものだ。定数がわからないという状態は，第三法則を持たない熱力学の正常な姿だったわけである。

定数の値そのものは実験で確かめられている。ザックールとテトローデが1912年に式を出したとき，その検証に使われたのは水銀蒸気の蒸気圧だった。単原子分子の気体の蒸気圧には，この定数がそのまま効いてくる。古典熱力学の測定量の中に量子定数が顔を出す，最も早い場面のひとつである。

こうして四年三ヶ月前の問いには決着がついた。ただし，そこに至る道で立てた他の問いはまだ立ったままである。ひとつは等重率の原理の根拠で，これは2022年4月に，エルゴード定理を持ち出す教科書が田崎さんの本の以降に減ったという形で書きつけたきり，今回のカノニカル分布の導出でも前提として素通りしている。もうひとつは準静的過程と可逆過程の関係で，教科書五冊を並べて定義が食い違うことを確かめただけで終わっている。

後者については，もう少し古い記憶がある。大学一年のとき，クラス担任だった国富信彦先生の授業で，応力テンソルと準静的過程につまずいた。覆水盆に返らずの話をしながら準静的過程の説明が進むので，これは可逆なのか不可逆なのかと混乱したのだった。五十年以上たって，同じ場所にまだ立っている。定数のほうは片づいたが，こちらは片づく気配がない。

カノニカル分布

ミクロカノニカル分布からの続き

図：カノニカル分布に対してChatGPTが考えたイメージ

カノニカル分布とは，孤立系を巨大な自由度をもつ熱浴と小さな部分系に分割したとき，熱浴によって一定温度 T に保たれた部分系の状態に対して定まる確率分布である。部分系と熱浴は自由にエネルギー交換ができるが，部分系の粒子数や体積は一定である。

１．カノニカル分布の導出

孤立系のエネルギーを$E_{\rm tot}$，熱浴のエネルギーを$E_{\rm B}$，状態$i$の部分系のエネルギーを$E_i$とすると，$E_{\rm tot} = E_{\rm B} + E_i$が成り立つ。全体系の状態数 $P(i)$は，ほとんど熱浴の状態数$\Omega_{\rm B}(E_{\rm B})$で決まるので，$P(i) \propto \Omega_{\rm B} (E_{\rm tot} - E_i)$ である。

全系は孤立系なので，ミクロカノニカル集団となって，各エネルギーの状態が等確率で実現される。このとき，統計力学によって熱浴のエントロピーは，$S_{\rm B}(E) = k_{\rm B} \log \Omega_{\rm B} (E) $ と定義される。

注：これは，独立な2系に対して状態数は積になる（$\Omega_{\rm AB} = \Omega_{\rm A}\cdot \Omega_{\rm B}$）がエントロピーは加法性を満たす（$S_{\rm AB} = S_{\rm A} + S_{\rm B}$）べきだという要請からきている。

これから，$\Omega_{\rm B}(E) = \exp \Bigl(\dfrac{S_{\rm B}(E)}{k_{\rm B}}\Bigr)$

したがって，$P(i) \propto \exp \Bigl(\dfrac{S_{\rm B}(E_{\rm tot}-E_i)}{k_{\rm B}}\Bigr)$

ここで，$E_{\rm tot} >> E_i$ から，エントロピーを$E_{\rm tot}$のまわりで展開すると，$S_{\rm B}(E_{\rm tot}-E_i) \simeq S_{\rm B}(E_{\rm tot}) - \Bigl( \dfrac{\partial S_{\rm B}}{\partial E} \Bigr)_{E_{\rm tot}} E_i$

また，熱力学から，$\dfrac{\partial S}{\partial E} = \dfrac{1}{T}$ なので，

$S_{\rm B}(E_{\rm tot}-E_i) \simeq S_{\rm B}(E_{\rm tot}) - \dfrac{E_i}{T}$

したがって，$P(i) \propto \exp \Bigl( \dfrac{S_{\rm B}(E_{\rm tot})}{k_{\rm B}} \Bigr) \cdot \exp \Bigl( - \dfrac{E_i}{k_{\rm B}T} \Bigr)  \propto e^{-\beta E_i}$， $\beta = \dfrac{1}{k_{\rm B} T}$

正規化された確率分布は， $P(i) = \dfrac{e^{-\beta E_i}}{Z}$となる。

ただし，$Z=\sum_i e^{-\beta E_i}$ は分配関数である。

２．量子系の場合

古典系では状態 $i$ は位相空間上の点に対応していたが，量子系では状態はヒルベルト空間上のベクトルで表され，エネルギーはハミルトニアン $\hat{H}$ の固有値として与えられる．すなわち，$\hat{H} \, |n\rangle = E_n \, |n\rangle $ を満たす固有状態 $\{|n\rangle\}$ が存在する。

部分系に注目するために，全体系の密度演算子を $\hat{\rho}_{\rm tot}$ とし，熱浴の自由度をトレースアウトすると，$ \hat{\rho} = \mathrm{Tr}_{\rm B} \, \hat{\rho}_{\rm tot} $，この $\hat{\rho}$ が部分系の統計状態（縮約密度行列）である。

カノニカル分布において，量子系における状態は次の密度演算子，$\hat{\rho}  = \Sigma_{n} P(n) | n \rangle  \langle n |$ で表される。ここで，$ P(n) = \dfrac{1}{Z} e^{-\beta E_n} $ であるから，$\hat{\rho}  = \dfrac{1}{Z}  \sum_{n} e^{-\beta E_n}  | n \rangle  \langle n | = \dfrac{e^{-\beta \hat{H}}}{Z} $ となる。分配関数は，$Z = \sum_n e^{-\beta E_n} =   \mathrm{Tr}  e^{-\beta \hat{H}}$ となる。

任意の物理量 $\hat{A}$ の熱平均（期待値）は，$ \langle \hat{A} \rangle = \mathrm{Tr} \, (\hat{\rho} \hat{A}) = \dfrac{\mathrm{Tr} \, (e^{-\beta \hat{H}} \hat{A})}{Z} $ で与えられる。

特に内部エネルギーは，$U = \langle \hat{H} \rangle = - \frac{\partial}{\partial \beta} \log Z $ となる。

さらに，エントロピーはフォン・ノイマンエントロピーとして $ S = -k_{\rm B} \, \mathrm{Tr} \, (\hat{\rho} \log \hat{\rho}) $ で定義される．カノニカル分布を代入すると，$ S = k_{\rm B} (\log Z + \beta U) $ が得られる。この関係は，ヘルムホルツ自由エネルギー $F  = -k_{\rm B} T \log Z $ を導入すると，$F = U - TS $という熱力学関係式と一致する。

社会物理学

図：未来都市のデータネットワーク（ChatGPTによる）

Ａ（Claude + ChagGPT）：

社会物理学：歴史・体系・原理

はじめに：定義と分野の性格

「社会物理学」を，物理学に由来する数理モデル・統計力学・非線形力学・ネットワーク科学・確率過程などを用いて，広義の社会現象を研究する領域と定義する。

ただし最初に確認しておくべきは，社会物理学が物性物理学のような統一された一分野ではなく，複数の研究領域が重なり合った「研究連合」に近いということである。狭義の sociophysics（スピン模型や投票者模型を直接用いる研究）と，経済物理学・都市科学・計算社会科学まで含めた広義の social physics とでは，規模が二桁以上違う。この二重性—看板としての分野は小さく，実質としての領域は巨大—は，以下の歴史・体系・原理のすべてに影を落としている。

１．研究史

(1) 前史(19世紀)：社会にも「法則」があるという構想

「社会物理学」という語自体は古い。サン＝シモンとオーギュスト・コントは，社会現象も天文学や物理学と同様に法則的に研究できるとして physique sociale という言葉を用いた。ほぼ同時期にアドルフ・ケトレーが『人間について——社会物理学試論』（1835）で，犯罪・出生・結婚・身体測定などの統計的規則性を研究し，「平均人」の概念を提唱した。ここでいう物理学は力学というより，多数の個人を集めると統計的規則性が現れるという発想である。コントが後に自らの学問を「社会学 sociologie」と呼び直した背景には，ケトレーが「社会物理学」の名称を統計学的研究に使ってしまったという事情もあった。つまりこの分野は，命名をめぐる縄張り争いとともに生まれている。

(2) 20世紀前半：統計的規則性の発見

20世紀前半には，理論というより経験則の蓄積が進んだ。ルイス・リチャードソンは戦争・武力衝突の規模分布がべき則に従うことを示し，ジップは言語の語彙頻度と都市規模の法則（1949）を定式化した。物理学者エットーレ・マヨラナが1942年に（死後）発表された「社会科学における統計法則の価値」は，量子統計の発想を社会に持ち込む先駆的エッセイとして知られる。また人口移動や都市間交流をニュートンの万有引力に似た式で表す「重力モデル」，確率過程による人口動態，交通流の流体・気体モデルもこの時期から中盤にかけて発展した。

(3) 数理社会科学の並行系譜

一方，社会科学の内部でも，物理学者によらない重要な系譜が育っていた。シェリングの居住分離モデル（1971），グラノヴェッターの閾値モデル（1978），アクセルロッドの文化伝播モデル（1997）は，いずれも単純な個人規則から予想外の集団現象が生じることを示すエージェント型モデルであり，後の社会物理学の中心的発想を先取りしている。社会物理学の歴史を物理学者の系譜だけで語ると，この半分が抜け落ちる。

(4) 現代的 sociophysics の成立（1970〜80年代）

現代的な意味での社会物理学は，ヴォルフガング・ヴァイトリッヒらが確率過程とマスター方程式による「社会力学」を展開した1970年代に始まる（世論形成のイジング模型的定式化は1971年）。1982年にはセルジュ・ガラムらが，スピン模型や相転移の考え方を集団意思決定・投票・連合形成へ応用し，sociophysics の旗を明示的に掲げた。個人をスピンに，賛否を上下スピンに，同調圧力を相互作用に，マスメディアを外場に対応させる研究様式がここで確立する。ただし当時この試みは物理学コミュニティ内で「不真面目」と見なされ，ガラム自身が長く冷遇されたと回想しているように，周縁的な営みだった。

(5) 離陸期（1990年代）

計算機シミュレーションの普及とともに状況が変わる。ナーゲル＝シュレッケンベルクのセルオートマトン交通流モデル（1992），ヘルビングの歩行者の社会力モデル（1995）が登場し，金融市場の価格変動・所得分布・企業規模・取引ネットワークを統計物理で扱う「経済物理学 econophysics」（1995年命名，マンテーニャとスタンレーら）が独立した大領域となった。自己組織化臨界（SOC）の流行も「べき則の遍在」への関心を後押しした。協力行動と進化ゲームの空間モデル研究もこの時期に急速に発展している。

(6) ネットワーク革命（1998〜2009）

1998年のワッツ＝ストロガッツのスモールワールド模型と1999年のバラバシ＝アルバートのスケールフリー模型は決定的な転機だった。社会関係を「平均的な相互作用」ではなく，誰と誰が結ばれているかというネットワーク構造として扱う流れが確立し，感染・情報・噂・技術・政治運動の拡散がネットワーク構造に大きく左右されることが示された。

2000年代には意見ダイナミクスの標準モデル群—Sznajd模型（2000），Deffuantらの限定信頼模型（2000），Hegselmann–Krause模型（2002）—が出揃い，2009年の Castellano・Fortunato・Loreto による Reviews of Modern Physics の大部の総説「Statistical physics of social dynamics」が分野の成立を象徴する画期となった。

(7) データ駆動期（2009〜）

同じ2009年，Lazer らが Science 誌で「計算社会科学 Computational Social Science」を宣言する。携帯電話・GPS・クレジットカード・検索履歴・SNSから得られる大量の行動データにより，人間移動・都市活動・コミュニケーション・オンライン世論が実証的に研究できるようになった。ペントランドの『ソーシャル物理学』（2014）はこのビッグデータ路線を一般に広めた本だが，統計物理系の sociophysics とはやや別系統であることには注意がいる（同じ social physics の看板の下に，異なる知的系譜が同居している）。2010年代後半以降は分極化・エコーチェンバー研究が盛んになり，COVID-19 が感染症モデルと社会行動モデルの融合を一気に加速させた。

(8) 第五段階：モデルと実証の接続（2010年代後半〜現在）

現在の中心課題は，単純模型を作ること自体から，模型と現実を接続することへ移っている。具体的には，観測データから相互作用規則を推定すること，オンライン実験や自然実験で因果関係を検証すること，ネットワーク自体の変化(共進化)を扱うこと，推薦アルゴリズムと人間行動のフィードバックを扱うこと，そしてAIエージェント（近年はLLMエージェント）を社会模型に組み込むことである。現代の代表的総説は都市・交通・金融・協力・犯罪・移住・感染症・環境問題・AIとの共存までを射程に収めており，たとえば2022年の包括的総説（Jusup ら，Physics Reports）は148ページ・参考文献千件超に及ぶ。もはや一つの模型群というより複合領域である。

２．分野の体系と規模

(1) 規模測定の困難：看板と実質の乖離

規模を測るときには根本的な注意が必要である。Scopus で「sociophysics」という語そのものを検索すると，2025年半ばまでで443件にすぎない。しかし social networks，human mobility，evolutionary games，urban studies，computational social science まで含めると数万から十万件級になる。つまり看板としての社会物理学は小さいが，実質的な研究領域は非常に大きい。この乖離を無視して単一の数字を出すことはできないため，規模は次の三層で把握するのが誠実である。

第一層は狭義の社会物理学。スピン模型・投票者模型・限定信頼模型などを直接用いる研究で，数百〜数千報規模。

第二層は中核的な広義の社会物理学。経済物理学・ネットワーク科学・都市物理学・人間移動・進化ゲーム・交通・群衆を含み，数万報規模。

第三層は計算社会科学まで含む最広義の領域。SNS分析・デジタル行動・因果推論・機械学習を含み，十万報を大きく超える—ただしこの全体を「物理学」と呼ぶのは，いささか領土拡張的である。

参考までに個別の数字を挙げると，狭義の sociophysics 検索は443報，「社会的影響」という広い OpenAlex キーワードでは約2万件，進化ゲーム理論のある書誌調査では約7800報，「ネットワーク上の世論」に限定した書誌調査でも約1400報，都市研究・都市計画全体では10万報以上(うち物理学的研究は一部)となる。境界の引き方で二桁動くことがよく分かる。

(2) 第一の軸：研究対象による分類

以下は厳密な排他的分類ではなく，周辺分野を含めた文献量・研究人口・専用学会や雑誌の充実度による相対評価である。

| 領域 | 主な対象 | 典型的な手法 | 規模 |

|-------|--------------|-------------------|--------|

| 社会ネットワークと拡散 | 人間関係，SNS，情報・噂・技術の普及 | グラフ理論，パーコレーション，分岐過程 | 最大級 |

| 経済物理学・金融市場 | 価格変動，企業，所得・富の分布 | 確率過程，時系列，極値統計，エージェント模型 | 最大級 |

| 都市・交通・人間移動 | 都市規模，渋滞,通勤,移住 | スケーリング，重力模型，流体・運動論 | 最大級 |

| 感染症と行動 | 感染，接触，予防行動，リスク認知 | SIR系模型，接触ネットワーク，適応行動模型 | 大(COVID以降急拡大) |

| 協力と進化ゲーム | 協力，公共財，規範 | レプリケーター方程式，空間・ネットワークゲーム | 大 |

| 世論・意思決定・文化 | 合意，分極化，投票，言語，文化 | イジング，投票者模型，限定信頼模型，ABM | 大(狭義の本丸) |

| 群衆・歩行者・避難 | 群集流，避難，安全 | 社会力模型，セルオートマトン，連続体模型 | 中〜大 |

| 科学の科学 | 引用，研究者移動，共同研究，イノベーション | 成長模型，優先的選択，計量書誌学 | 中〜大 |

| 犯罪・紛争・テロ・移民 | 犯罪集中，暴動，戦争，連合形成 | 自励過程(Hawkes),空間模型,ゲーム | 中 |

| 環境・エネルギー・気候行動 | 消費，政策受容，社会的転換(tipping) | エージェント模型，カスケード | 中 |

| 人間・アルゴリズム・AI | 推薦，フィルターバブル，LLMエージェント社会 | 適応ネットワーク，強化学習，LLM実験 | 新興・急成長 |

(3) 第二の軸：方法による分類

対象ではなく方法で切ると，社会物理学は五つの群に整理できる。

第一群は，統計的規則性の発見（べき分布・スケーリング・揺らぎ・相関）

第二群は，相互作用するエージェントの模型（スピン・投票者・ゲーム参加者・歩行者)

第三群は，ネットワーク上の動力学（感染・情報拡散・合意形成・同期・カスケード)

第四群は，流れと空間の物理（人口移動・交通・都市・群衆）

第五群は，データ駆動型社会力学（巨大データ・統計推定・機械学習・オンライン実験・因果推論)

歴史的には第一群（ケトレー〜ジップ）から第二群（ガラム以降）へ，そして現在は第二群単独から，ネットワーク・データ・因果推論を組み合わせる第三〜第五群へと重心が移っている。この重心移動こそ,「単純模型の時代」から「モデルと実証の接続の時代」への転換の内実である。

主要な発表媒体は Physica A，Physical Review E，EPJ B / EPJ Data Science，Journal of Statistical Mechanics などの物理系誌に加え，近年は Nature Human Behaviour や PNAS など学際誌の比重が増している。

３．重要で普遍的な原理と共通の視点

前置き：「社会版ニュートンの法則」はまだない

まず結論を述べれば，社会物理学には，エネルギー保存則や運動量保存則に相当する，社会全体で成立する厳密な基本法則は確認されていない。存在するのは，異なる現象を共通の形式で見るための中間レベルの原理の集合である。以下の七つがその骨格をなす。

① ミクロな相互作用からマクロな秩序が創発する

個人の行動規則が単純でも，多数が相互作用すると，合意・分極化・居住分離・バブル・暴落・渋滞・群集事故・協力規範などが自発的に生じる。これが社会物理学の最も根本的な視点であり，統計力学の中心思想の移植である。マクロな秩序を，誰かの設計や個人の性格だけで説明せず，相互作用の帰結として説明する。シェリング型の居住分離が示すように，強い差別意識を持つ個人がいなくても，わずかな近隣選好から強い空間的分離が生じうる。社会物理学が社会科学に与えた最大の貢献は，「単純な原因と単純な結果を結ぶ直線的説明」を崩したことだろう。

② 相互作用の「幾何学」はネットワークである

物理系では空間的に近い粒子が相互作用するのに対し，社会では物理的距離よりも，家族か，友人か，同じ職場か，SNSで接続されているか，同じ情報源を見ているかが重要である。したがって社会物理学ではネットワークが物理空間に相当し，同じ個人規則でも，完全混合・格子・スモールワールド・モジュール構造・中心集中型ネットワークでは結果が全く異なる。典型例として，理想的なスケールフリーネットワーク上では感染閾値が（熱力学極限で）消失する。「誰と誰がつながっているか」の構造が，個々の性質と同等以上にマクロな帰結を決める—平均場近似が破れる場面こそ社会的に重要，という視点である。

③ 非線形フィードバックが相転移・臨界現象を生む

社会には正のフィードバックが遍在する。人気があるからさらに人気が出る。売られるから不安が増し，不安だから売られる。推薦されるから閲覧され，閲覧されるからさらに推薦される。このため比例関係（線形応答）よりも，閾値・複数安定状態・履歴依存（ヒステリシス）・臨界転移が重要になる。合意⇄意見分裂，自由流⇄渋滞，流行の発生⇄鎮静は，いずれも制御パラメータ（結合強度・密度・感染率）の変化に伴う秩序−無秩序転移として記述できる。小さな変化が何も起こさない範囲と，突然大きな転換を起こす範囲がある—社会の質的変化（革命，流行の臨界点）が「量的変化の閾値現象」として理解できるというのが，相転移とのアナロジーが最も有効に働く部分である。

 ④ 平均よりも分布と異質性が重要である——ただしべき則信仰には留保つきで

社会を「平均的個人」で代表させると多くの場合に失敗する。所得・影響力・友人数・企業規模・投稿頻度・移動距離は強く歪んだ分布を持ち，影響を与えやすい人，頑固な人（zealot），橋渡しをする人が少数存在するだけで全体の挙動が変わる。ジップ則，パレート分布，都市スケーリング，戦争・地震・SNSカスケードの規模分布—特徴的スケールを持たない分布の遍在は，社会が平衡から遠い状態にあることを示唆する。

ただし，ここには重要な留保が必要である。初期ネットワーク科学では「スケールフリー」がやや乱発されたが，Broido と Clauset の大規模検証（2019）以降，純粋なべき分布に従う実ネットワークは当初考えられたほど一般的でないことが指摘されている。都市スケーリング則でさえ，都市境界の定義や都市内部の格差の扱いによって指数や解釈が変わる。「べき則の遍在」は社会物理学の看板であると同時に，最も過剰宣伝されやすい主張でもあった。分布の重い裾（heavy tail）と異質性の重要性は揺るがないが，それが厳密なべき則かどうかは個別の統計的検証の問題である。

⑤ 社会は非平衡・開放・適応系であり，状態とネットワークが共進化する

社会は熱平衡に向かう閉じた系ではない。情報・人口・資源・制度・技術が絶えず出入りし，人々は経験から学習する。詳細釣り合いは成り立たず，定常状態は平衡分布ではない。ゆらぎは単なるノイズでなく，対称性の破れ（どちらの意見に収束するか）や希少事象（バブル崩壊）を駆動する建設的役割を果たす。自己組織化臨界のように，外部調整なしに系が臨界近傍に居続けるという描像も，この非平衡性の一つの表現である。

さらに重要なのは，個人の状態だけでなく相互作用ネットワークそのものが変化することである。意見が似るから友人になり（選択・ホモフィリー），友人だから意見が似る（社会的影響）。状態がネットワークを変え，ネットワークが状態を変えるという共進化（adaptive networks）が起こる。この二方向の因果は分極化・エコーチェンバー研究の中心的難問（ホモフィリーと影響の分離不可能性）に直結しており，平衡統計力学よりも，非平衡系・適応ネットワーク・確率過程・学習力学のほうが社会には本質的である。

⑥ 普遍性のジレンマ：普遍性は単純模型を正当化するが，機構の同定を不可能にする

統計物理の臨界現象には，ミクロな詳細がマクロな振る舞いに影響しない「普遍性クラス」が存在する。これは「人間の複雑な心理を無視した単純化しすぎたモデル」という批判への，社会物理学側の最も強い理論的正当化である。人間が原子より複雑でも，集団の統計的振る舞いは少数の本質的要素（対称性，相互作用の範囲，保存量）で決まりうる。

しかし，この同じ性質が裏返ると方法論上の最大の難点になる。ミクロの詳細がマクロに効かないなら，マクロな曲線の再現からミクロ機構を同定することはできない。実際，分極化という同一のマクロ現象は，同類選好からも，反発作用からも，マスメディアからも，推薦アルゴリズムからも，経済格差からも，地域的分離からも生じうる（等結果性）。したがって「模型が分極化を再現した」ことは，その模型が正しい説明であることを意味しない。物理模型は「その現象が単純な相互作用からでも生じうる」ことを示すには強力だが，「現実にその機構で生じた」ことの証明にはならない。再現は説明ではなく，生じうるは生じたではない。

このジレンマ—普遍性は武器であると同時に制約である—こそ，社会物理学の方法論の核心である。ここを乗り越えるには，マクロ曲線へのフィッティングを超えて，実験・因果推論・個人レベルの時系列・複数の独立データによる検証が必要になる。現代の計算社会科学が単なる巨大データ分析から因果的・実験的研究へ向かっているのは，まさにこのためである。

⑦ 観測者と模型が対象を変えてしまう（反射性）

社会系には物理系にない強い反射性（reflexivity）がある。選挙予測が発表されれば投票行動が変わり，信用評価モデルが使われれば人々は評価基準に適応し，推薦アルゴリズムが導入されればコンテンツ制作者も利用者も行動を変える。つまり社会物理学の模型は，対象を外から受動的に記述するだけでなく，導入された瞬間に系の一部になりうる。これは予測市場・アルゴリズム統治・世論誘導では決定的であり，「モデルの予測自体が予測対象を変える」という点で，量子力学の観測問題より厄介かもしれない—観測者効果が物理定数ではなく，学習し適応する主体を介して働くからである。

４．総括：社会物理学の「普遍性」をどう考えるべきか

社会物理学で真に普遍的なのは，特定の数値や法則というより，次の認識枠組みである。

多数の異質な主体が，ネットワーク上で相互作用し，ノイズ・外部環境・学習・制度とのフィードバックを通じて，非線形なマクロ現象を生む。

この枠組み—社会を，多体系・ネットワーク・非線形系・非平衡系・適応する複雑系として見ること—はかなり広い範囲で有効である。一方，「都市の指標は必ず特定の指数でスケールする」「社会ネットワークは必ずスケールフリーである」「人間はイジングスピンとして扱える」といった個別の主張まで普遍法則とみなすのは危険である。

社会物理学の弱点は，物理学で成功した単純化を社会に移植しすぎると，意味・権力・制度・歴史・文化を捨象してしまうことにある。⑥で述べたとおり，物理模型は可能性の証明であって現実性の証明ではない。この区別を忘れると，社会物理学は科学というより，数式を着せた寓話になる—白衣を着せれば仮定が事実になるわけではない。

1982年のガラムの孤立から，2020年代のパンデミック政策への実装まで。「社会に法則はあるか」というケトレー以来の問いに対して，社会物理学は「厳密な基本法則はまだないが，中間レベルの原理群と，それを実データで検証する方法論は育ちつつある」という，部分的だが定量的な肯定を与えつつある。それが2020年代半ばの現在地である。