ラベル 物理 の投稿を表示しています。 すべての投稿を表示
ラベル 物理 の投稿を表示しています。 すべての投稿を表示

2026年10月12日月曜日

ナビエ−ストークス方程式(3)

ナビエ−ストークス方程式(2)からの続き


図:ナビエ−ストークス方程式を巡って(ChatGPTによる)


この10日間の動きを調べます。この10日間で大きく動きました。要点を先に:クレイ研究所が公式声明を出して「どうやら決着したらしい」と認めた一方で、フィールズ賞受賞者25人が連名でAI企業のやり方を批判する宣言を出した——この2つが同じ9月11日に起きています。

 ① クレイ数学研究所の公式声明(9月11日)── 最大の変化

前回お伝えした時点では「受理せず、未解決のまま掲載」でしたが、9月11日にCMIが正式な声明「Navier-Stokes Announcement」を出しました。要点は3つ。

- CMIは「ナビエ・ストークス問題がどうやら決着したという発表」に世界の数学界とともに興奮を共有すると述べ、この仕事の背後にある革新が分析され問いただされることで「人間の新しい理解の波」が解き放たれることを期待する、としています。
- 賞の規則が「何が達成されたかを評価し、功績を割り当てる」手続きを定めており、そのプロセスは「意図的に急がない」が、随時アップデートを出すとしています。
- また、近年この分野周辺の進展(一部はクレイ研究賞で顕彰)によって問題が間もなく解決されるかもしれないという期待が高まっており、新技術が数学研究を加速する能力の向上がその期待をさらに強めた、とも書いています。

**重要な読み方**:この声明は**OpenAIの名を一度も挙げていません**。「apparently been settled(どうやら決着した)」という受動態で、主語は伏せられている。しかも規則が定めるのは「評価」と「**功績の割り当て**(assigning credit)」の両方だと明示しており、これは優先権争いを正面から視野に入れた書き方です。そして問題ページ自体は今も「**Active**」のタグのままです。

つまり **「解かれたらしい」とは認めたが、「誰が解いたか」も「本当に正しいか」もまだ何も決めていない**、というのが9月11日の到達点です。規則上、査読誌掲載+2年の一般受容期間が必要なので、公式決着は早くても2028年以降になります。


 ② フィールズ賞受賞者25名の共同宣言(9月11日)── 最大の余波

同じ日に、mathandai.org で「**A Severe Misalignment of AI in Mathematics**(数学におけるAIの深刻なミスアライメント)」という宣言が公開されました。Terence Taoが自身のMastodonとブログで公表し、Leiden宣言と同様に追加署名を受け付けています。Taoによれば、この声明は受賞者たちの前週の議論から生まれ、緊急性が高いと判断して急いで公開されたとのことです。

**署名者(初期25名)**にはTao、Scholze、Villani、Deligne、Kontsevich、Bhargava、Viazovska、Hairer、Figalli、Maynard、Huh、Ngô Bảo Châu、そして**森重文**(1990年)が並びます。注目すべきは **Pierre-Louis Lions(1994年)** の名前があること——彼自身が圧縮性ナビエ・ストークスの弱解理論の構築者であり、まさに当事者領域の人です。最古参のDeligne(1978年)から今年受賞のYu Deng(2026年)までが並んでいます。

**宣言の論旨**(要約すると):

1. LLMの数学能力は劇的に向上し、多くの分野で主要な未解決問題を解けるところまで来た。しかし**AI企業が数学の問題をベンチマークとして解こうとする動きは、数学という科学と数学コミュニティにとって有害**であり、両者の目標は深刻にずれている。
2. 有名問題は**灯台や標識**として機能してきた。解けたこと自体より、そこで生まれた新しい洞察と手法が、講演・議論・簡略化という長く骨の折れる過程を経てコミュニティに吸収され、最終的に教科書になることに意味がある。
3. **問題を解くことは、概念的理解という本来の目標の道具・代理にすぎない**。「真/偽」の言明をどんどん速く大量生産することは、新しいアイデアに息を吹き込むどころか、肥沃な土壌を破壊しかねない。
4. 急いで発表されるため、きちんとした書き上げ、新手法の抽出、**他者の先行研究の引用**の時間がない。これは他の創造的職業と同様、**深刻な帰属・剽窃の問題**を引き起こす。
5. AIが生んだアイデアの世話をし数学の正典に統合する人間がいなければ、そのアイデアは本当には生きず、**数学者どうしの人間的な継承の連鎖**が失われる。
6. ただし**AI反対の宣言ではない**。AIには本物の数学的研究と理解を高め加速する潜在力がある、と明記しています。決定的なのは、この技術を制御する人間の選択だ、と。

**押さえるべきニュアンス**:宣言はOpenAIの証明が間違っているとは一言も言っておらず、企業名も挙げていません。争点は**正しさではなく、プロセス・帰属・人間の理解**です。この区別を外すと議論が噛み合わなくなります。

背景記事として、Taoは宣言と併せてEconomist誌の記事「Top mathematicians are outraged by OpenAI's methods」を挙げています。


③ 引用漏れ問題 ── 宣言の主張を裏書きする実例

これは小さいようで象徴的な出来事です。OpenAIの当初の論文は、CórdobaとMartínez-Zoroaの論文も、AlpögeとBuckmasterの仕事も引用していませんでした。EPFLのGonzalo Cao-Laboraが参考文献の乏しさを指摘した後、9月8日中に更新版でCórdobaらの論文が追加されました。改訂後のPDFは166ページになり、歴史的経緯の記述が拡張されて、Córdoba–Martínez-Zoroaの先行構成、Zhengとの関連研究、古典的な波の不安定性研究への言及が加わりました。主定理と周期版の系は変わっていません。

つまり宣言が言う「急いだ発表による引用の欠落」が、まさに引き金となった事例そのものに起きていた。数時間で直ったとはいえ、**外部の数学者に指摘されるまで気づかれなかった**という事実が残ります。


 ④ 検証の現状 ── 今どこまで確かめられているか

前回の時点で「誰も証明を見ていない」という報道がありましたが、これは現在は誤りです。公開されているものを整理します。

|公開物―――|内容―――――――――――|
|NS論文  |166ページ(改訂後)|
|Euler論文 |別立ての57ページの原稿。外力なし・滑らかな初期データからの全空間での爆発を主張|
|Lean形式化|github.com/openai/NavierStokesAndEuler。Lean 4.34.0-rc2 + Mathlib。リポジトリの指示はキャッシュ取得と `lake build` の実行。外部チェッカーで再検証する手順の文書も別途付属|

誰でもクローンして数分でビルドでき、`sorry`(未証明の穴)の有無も自分で監査できる状態です。

しかし、これで終わりではないという点が繰り返し強調されています。Lean証明書が保証するのは「Leanとその読み込むライブラリの中で、形式的な定義から形式的な言明が導かれる」ことだけで、それらの定義が数学者の読むClay問題そのものであることは保証しないし、専門家が166ページの本文を読む作業の代わりにもならない。

現時点での判定:
・誤りの指摘は公には出ていません。
・Clay問題との対応については、問題文の原著者(Fefferman)から励ましになる予備的評価は得られているが、コミュニティの独立した受容はまだ。
・Palasekが指摘した障害は、あくまで「この構成から外力を取り除くルート」に対するものであり、Clayの仮定を満たす外力付きNSの結果を反証するものでも、あらゆる外力なしのアプローチを封じるものでもない。


 ⑤ 論争のその後 ── NYT報道などで出た新事実

ニューヨーク・タイムズ(Kenneth Chang記者)の取材で、いくつか新しい情報が出ています。

・OpenAI側にPDEの専門家がいなかった。Bubeckは、数学者のチームはいるが誰も偏微分方程式の専門家ではないと認めています。「問題文を与えただけで数日で到達する方向ではない」というBuckmasterの論点に対しては、これは有利にも不利にも読めます(内部に土地勘がなかったからこそ噂の方向づけが効いた、とも、専門家抜きでこれを出したのはモデルの力だ、とも)。
・提示された条件の全容。Bubeckが提示したのは、Buckmasterが先に結果を出してよい、OpenAI論文の筆頭著者になってよい、潤沢な計算資源を提供する、さらにBuckmasterらこそ100万ドルに値すると表明する、という内容だったと報じられています。決裂の核心は「Leventを著者に入れたくない」という一点だった。
・数学界の反応。プリンストン高等研究所のPeter Sarnakは、この種の駆け引きは聞いたことがない、彼が数学者らしく振る舞ったことを称賛すると述べています。
・Wikipedia に独立項目が立った。「Navier–Stokes priority controversy」という専用記事ができ、ミレニアム問題やナビエ・ストークスの各記事も「クレイ研究所にも独立した数学コミュニティにも未検証で、優先権争いの対象」と記述されています。歴史的事件として整理され始めた、ということです。


⑥ まだ起きていないこと(=今後の注目点)

10日経って、動いていないことの方がむしろ重要かもしれません。

(1) 証明の誤りは見つかっていない。同時に、専門家による本格的な読解結果もまだ公表されていません。166ページ+57ページ+Leanコードを読み切るには数ヶ月単位がかかります。
(2) データアクセス問題は未決着。OpenAIは9月10日の調査結果追記でBuckmasterのCodexプロンプトの影響を否定しましたが、非識別化データの寄与は排除できないという留保はそのままです。第三者による検証の道はありません。
(3) 外力を消す問題(A/B)は完全に未着手。ここが本丸であることは変わっていません。
(4) 2つ目のミレニアム問題の噂が流れていますが(Hodge予想やBSD予想の名が挙がる)、確認された情報ではありません。現時点では噂として扱うべきです。
(5) CMI(クレイ数学研究所)の次のアップデートがいつ出るか。声明は「随時アップデートする」とだけ書いています。


⑦ブログにまとめるなら

前回申し上げた3つの軸(C/D vs A/B、アイデアの出所は人間、答えと理解の脱結合)は、この10日でむしろ公式に裏書きされた形です。特に:

・「答えと理解の脱結合」というTaoの言葉は、25人のフィールズ賞受賞者の連名宣言という形で制度化されました。宣言の「有名問題は灯台であって、ベンチマークではない」「問題を解くことは理解という目的の代理にすぎない」という論理は、日本語圏でもそのまま通用する強い主張です。
・「アイデアの出所は人間」は、引用漏れとその即日修正という具体的な出来事として現れました。
・一方で「AIは何も生み出していない」式の論も同様に不正確で、Clayが「どうやら決着した」と書いた事実、Lean形式化が誰でもビルド可能な状態で公開されている事実は動きません。

(注:Claudeが)個人的に面白いと思うのは、CMIの声明とフィールズ賞宣言が同じ日に出て、正反対の温度を持っていることです。制度の側(CMI)は「規則に従って粛々と、功績の割り当ても含めて評価する」と落ち着いて構え、研究者の側は「このままでは職業としての数学が壊れる」と警鐘を鳴らした。この温度差そのものが、今回の事件の本質を映しているように思います。

出典:
[1]Navier-Stokes Announcement — Clay Mathematics Institute
[9]The Navier-Stokes controversy — Singham NYT報道の引用

2026年10月11日日曜日

ナビエ−ストークス方程式(2)

ナビエ−ストークス方程式(1)からの続き


図:ナビエ−ストークス方程式(ncu_inoue_labから引用)

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

A(Claude):

1.整理の枠組み:条件は「2値」ではなく「連続的な階段」

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

|軸――――――|パラメータ――――――|効果――――――――――――――――|
|① 空間次元  |2D / 3D       |最大の分岐。渦伸長項の有無     |
|② 散逸の強さ|$(-\Delta)^\alpha$ の $\alpha$(0=Euler, 1=NS)|大きいほど易しい|
|③ 外力   |$f$ の正則性クラス(超臨界→…→$C^\infty$急減少→$f\equiv0$)|滑らかにするほど難しい|
|④ 圧縮性  |非圧縮 / 圧縮性    |困難の種類が変わる         |
|⑤ 境界   |アリ / ナシ      |境界アリの方が易しい        |
|⑥ 解のクラス|古典解・強解 / Leray–Hopf弱解 / 一般の弱解|弱くするほど「易しい(が非一意)」|
|⑦ 対称性   |軸対称(swirlなし) など |対称性が強いほど易しい       |

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


2.支配原理:超臨界性がすべてを決める

分数階散逸 $(-\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万ドルが埋まっているわけです。


3.軸① 次元: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も自動的に大域正則です。

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


4.軸② 粘性:Eulerの方が難しいのか易しいのか(重要な逆説)

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

・爆発を証明したい → Eulerの方が易しい。粘性は特異点を潰そうと働くので、爆発を作るには粘性という敵と戦わねばならない。実際、爆発の証明はほぼすべてEuler側(ないし劣散逸側)で先に達成されています。
・大域正則性を証明したい → NSの方が望みがある。粘性のおかげでエネルギー散逸、Leray弱解の存在、Caffarelli–Kohn–Nirenbergの部分正則性といった構造が使える。Eulerにはこれらが一切ない。

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

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


5.軸③ 外力:ここが「階段」であることが最大の要点

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

|段|外力 $f$ のクラス―――――――|難易度――――――|状況――――――|
|0|超臨界($u$ から逆算した $f$)|ほぼ自明|爆発する $u$ を勝手に決めて $f:=\partial_t u+(u\cdot\nabla)u-\nu\Delta u+\nabla p$ と置くだけ。$f$ も一緒に爆発するので無意味|
|1|$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クラス内での爆発=方程式の力学が生む爆発であり、外力が押し潰しているのではない|
|2|空間的に滑らか(時間方向は粗い)|さらに難しい | Córdoba–Martínez-Zoroa の IPM(2024)|
|3|時空一様に $C^\infty$ かつ全導関数が急減少(Clayの条件(5))|非常に難しい|Alpöge–Buckmaster–Coiculescu(IPM)→ 2026年8月のBoussinesq・3D Euler → OpenAIのNS主張|
|4|$f \equiv 0$           |本丸|3D Euler:OpenAIが主張(未検証)/3D NS(A/B):完全に未解決|

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


6.軸④ 圧縮性:難しさの「種類」が変わる(同じ物差しに乗らない)

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

|――――――――――|特異点の種類―――――――――|状態――――――――――|
|圧縮性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)$ で非局所的に決まり、情報が無限速度で伝わる。圧縮性では有限伝播速度で双曲型。解析の道具立てが根本的に異なります。


7.軸⑤ 境界:あると易しい

境界の役割は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が主張している地点です。


8.軸⑥ 解のクラス:弱くすると「易しくなるが無意味になる」

これも独立した軸です。

・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$)を要求しているのはこのためです。


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

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

|モデル――――――――|性格――――――――|状態――――――――|
|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を入れた瞬間に未解決になる——この一点にも渦伸長の効き方が凝縮しています。


10.全体マップ:何がどこにあるか

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

|設定―――――――――――――――|状態―――――――――――――――|
|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の正則性            |❌ 未解決|


11.まとめると

(1) 難易度の主軸は「超臨界性の度合い」であり、次元・散逸指数・外力の正則性はすべてその一つの物差しの上に並ぶ。
(2) 3D NSは臨界点から $1/4$ だけ外れている。この微小なギャップが90年の壁。
(3) 「外力アリ/ナシ」は2値ではなく5段の階段で、Clayが要求する段3(時空一様に滑らか+急減少)と段0(逆算した外力)は天と地ほど違う。2026年9月のニュースは、段3に到達したという主張。
(4) 粘性ナシ(Euler)は「爆発を示すには易しく、正則性を示すには難しい」という逆説的な位置にある。
(5) 圧縮性は難易度の軸ではなく、問題の種類を変える軸。圧縮性の特異点(衝撃波)は証明済みかつ理論的に処理済みで、ミレニアム問題とは別世界の話。
(6) 残る最難関は依然として 3D NS、外力なし、境界なし、滑らかデータ、古典解(A/B)。そしてそこへ至る最大の技術的障壁は、外力を消したときに粘性散逸が成長機構を食い潰してしまう点にある。

2026年10月10日土曜日

ナビエ−ストークス方程式(1)


図:これがそれなのか?(ChatGPTによる)

Q(koshix):
ミレニアム懸賞問題のひとつである,ナビエ−ストークス方程式をAIが解いたという話が聞こえてきました。出来るだけ分かりやすく解説してほしいのです。
(1) 流体力学におけるナビエ−ストークス方程式の位置づけについて
(2) ナビエ−ストークス方程式はどこまでわかってきたか,何が難しいのか
(3) ミレニアム検証問題におけるナビエ−ストークス方程式に関するの正確な条件
(4) AIが何を解いて,それはどこまで画期的なことなのかそうではないの
(5) Buckmaster らを巡る混乱の争点を正確に解きほぐす
(6) 今後どうなるのか

A(ChatGPT + Claude の第一報):

0.全体像

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

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

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.どこまで分かっていて、何が難しいのか

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は「無限カスケード」を組み合わせて特異点を作る解析的手法を発展させた。ただし従来の構成では、合成した外力が十分に滑らかではなかった。残された難所は、特異形成のカスケードを維持したまま、外力を最後まで滑らかに保つことであった。

3.ミレニアム問題の正確な条件

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のほうが深い問題だと受け取られてきた。

4.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時代には、この二つが初めて大規模に分離し始めている。

5.Buckmasterをめぐる混乱

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

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

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

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

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

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

6. 今後どうなるか

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

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

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

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

今回の出来事を「AIがミレニアム問題を解いたか」という二値問題だけで見ると、本質を取り逃がす。重要なのは次の三点である。
・公式ルールと数学者の実感の乖離、
・人間が生んだアイデアとAIによる高速探索の関係、
・そして「答えを得ること」と「理解すること」の分離。

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


2026年10月9日金曜日

ニュートリノ(2)

ニュートリノ(1)からの続き


図:T2Kによるニュートリノ振動の研究(Geminiによる)


8.人工ニュートリノで振動を精密に測る: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年)、加速器で作った人工のニュートリノで「どう振動するか」を測る時代に入り、いまや「ニュートリノと反ニュートリノで振動が違うか」という宇宙の物質の起源につながる問いに挑んでいる、という流れになります。

2026年10月8日木曜日

自然界の左右対称性

ホモキラリティ(1)からの続き


図:宇宙から生命へキラリティの進化(ChatGPTによる)

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

A(ChatGPT):

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

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

ただし、「関連する研究をすべて」とすると、核物理学の対称性、結晶構造、生体分子の立体構造、さらには一般的な対称性の自発的破れまで対象が広がります。
そこで以下では、関連度を三段階に分類します。
◎ 直接関連:左右非対称性、鏡像異性、パリティの破れが研究の中心。
 ○  強い関連:左右対称性が研究の重要な要素。
△ 周辺関連:より一般的な対称性や立体構造の研究。左右性そのものが受賞理由ではない。

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

1.年代順一覧

受賞年・賞

受賞者と研究内容

関連度

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年の南部の研究は、対称性の破れの種類を区別する必要があります。


2.とくに重要な研究の流れ

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

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

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

2001年・化学賞:人工的に左右を作り分ける
ノールズ、野依、シャープレスは、キラルな触媒を使って一方の鏡像異性体を優先的に合成する方法を開発しました。特に医薬品製造で大きな意味を持つ成果です。

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


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

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

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

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


2026年10月7日水曜日

ニュートリノ(1)


図:ニュートリノ探究の歴史(Geminiによる)

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

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

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


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

A(Claude):

0.主役:ニュートリノとは

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

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

最大の発生源は太陽です。太陽の中心で水素がヘリウムに変わる核融合のたびに、電子ニュートリノが生まれます(太陽ニュートリノ)。
・地球に届く量:1 cm²あたり毎秒約650億個
・手のひら10cm四方なら、毎秒約6〜7兆個
ほとんど反応しないので、昼は空から、夜は地球を貫いて足元から、同じように体を素通りしていきます。人体で実際に反応するのは、一生に1回あるかないかです。
このほかに、地球内部の放射性元素から出るもの、宇宙線が作る大気ニュートリノ、ビッグバンの名残などがありますが、数は太陽のものよりずっと少数です。

2.最初の謎:太陽ニュートリノが足りない(1960年代末〜)

デービスは、地下の金鉱に塩素を含む約600トンの液体を置き、電子ニュートリノが塩素をアルゴンに変える反応を利用して、生じたアルゴン原子を数えました。結果は理論予想の約3分の1でした。
太陽のモデルが間違っているのか、実験が間違っているのか、それともニュートリノ自体に何かあるのか。これが30年以上続く「太陽ニュートリノ問題」の始まりです。

3.カミオカンデ:ニュートリノを「見る」装置(1983年〜)

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

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

数の多い太陽ニュートリノに、他の信号はかき消されないのか?
かき消されません。三つの物差しで区別できるからです。

(1) エネルギー:太陽ニュートリノの大半は0.4 MeV以下と低エネルギーで、約7 MeV以上しか見えない検出器にはそもそも映りません。見えるのは全体の約1万分の1だけです。しかも反応確率はエネルギーが高いほど大きいので、GeV級の大気ニュートリノと実際の事象数は同程度になり、エネルギー域も重なりません。
(2) 方向:太陽ニュートリノに弾かれた電子は、太陽と反対の方向へ飛びます。
(3) 時間:定常的に降る太陽ニュートリノに対して、超新星のニュートリノは短時間に集中して来ます。
低エネルギー側で本当に邪魔になるのは、太陽ニュートリノではなく、水や岩盤に含まれるラドンなどの放射性物質です。

二つの成果
・1987年 超新星1987A:16万光年先の超新星爆発から、13秒間に11個のニュートリノを捉えました。ニュートリノ天文学の誕生で、小柴昌俊さんの2002年ノーベル賞の理由になりました。
・1989年 太陽ニュートリノ:事象が太陽の方向から来ていることを示し、本当に太陽から来ていることをリアルタイムで確かめました。数はやはり予想の約半分で、デービスの不足が実験の誤りではないことも裏付けました。

4.第二の謎:大気ニュートリノの異常(1980年代後半)

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

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

異常の発見
梶田隆章さんは、カミオカンデの大気ニュートリノの事象を解析し、リングの形で2種類を見分けました。
・ミューオン(ミューニュートリノ由来)→ くっきりしたリング
・電子(電子ニュートリノ由来)→ ぼやけたリング
その結果、ミューニュートリノの割合が予想の2対1より明らかに小さいことが分かりました。ミューニュートリノがどこかで消えているように見えるこの現象が、「大気ニュートリノ異常」です。

5.決着①:スーパーカミオカンデ(1998年)

5万トンの水を使う後継機で、梶田さんたちは飛来方向別にニュートリノを数えました。
ここで「全方向から一様に来る」ことが決定的に効いてきます。
・真上から来るもの:飛行距離は約15 km
・真下から来るもの:地球を貫いて約1万3000 km
もとの量は上下でほぼ対称なはずなのに、結果は次のとおりでした。

|種類―――――――|上から(短距離)―|下から(長距離)|
|ミューニュートリノ|ほぼ予想どおり  |約半分に減少  |
|電子ニュートリノ |予想どおり    |予想どおり   |

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

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

6.決着②:SNO(2001〜2002年)と太陽ニュートリノ問題の解決

カナダのSNOは重水を使い、
・電子ニュートリノだけに反応する反応
・3種類すべてに等しく反応する反応
を同時に測りました。結果は、電子ニュートリノだけなら予想の約3分の1、3種類の合計なら予想どおりでした。

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

整理すると
・大気ニュートリノ:減ったのはミューニュートリノ(主にタウへ)
・太陽ニュートリノ:減ったのは電子ニュートリノ(ミューとタウへ)

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

7.次の舞台:アイスキューブと宇宙ニュートリノ(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 | ノーベル賞(ハルツェン)    |高エネルギーニュートリノの発見 |

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

2026年8月17日月曜日

ヴィルヘルム・レンツ

 イジング模型からの続き



写真:ヴィルヘルム・レンツ(ロストック大学から引用)

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

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

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


2026年8月16日日曜日

イジング模型


図:1次元イジング模型(ChatGPTによる)

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


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

A(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次元イジングモデルは低温ではスピンがそろいやすいが,有限温度では真の強磁性相転移を起こさない。

2026年8月14日金曜日

定数がわからない

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


図:エントロピーがわからない(ChatGPTによる)

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


A(Claude):

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

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


1.そこまでの道筋

始まりは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$ を含むから示量変数になっていることは確かめた。しかし定数が一致するのかどうかがわからない。


2.四年後の答え合わせ

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


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


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


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


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


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


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