tecktom.com

Lanford の定理と不可逆性

2026年10月9日8分で読めますAI生成教科書『虎の巻』第8章 Boltzmann 方程式から流体方程式へ
目次
  1. 主張
  2. 証明の骨格と難所
  3. 不可逆性の起源
  4. 木下の層

AI 生成・木下が未勉強

近似の数理 第8章 §8.3 ★★ Boltzmann 方程式が B 型の極限として導けること、ただし短時間に限ることを主張として押さえる

主張

剛体球系の BBGKY 階層は、kk 粒子周辺分布 fN(k)f_N^{(k)} の方程式が、第 k+1k+1 粒子との衝突を表す項 Ck,k+1fN(k+1)C_{k,k+1}f_N^{(k+1)} を通して fN(k+1)f_N^{(k+1)} を含む連鎖である。Boltzmann–Grad 極限で衝突項の係数 (N−k)d2(N-k)d^2 は一定値に収束し、形式的には極限で Boltzmann 階層と呼ばれる連鎖が残る。Boltzmann 階層は、各 kk について f(k)=f⊗kf^{(k)}=f^{\otimes k}(ff の kk 重の積)と置くと、ff が Boltzmann 方程式を満たすことと同値になる。

定理 8.6(Lanford の定理(1975) ★★)

f0(x,v)f_0(\bm{x},\bm{v}) を周期的な箱(または R3\mathbb{R}^3)上の連続な確率密度で、ある C,β>0C,\beta>0 について f0(x,v)≤Ce−β∣v∣2f_0(\bm{x},\bm{v})\le C\mathrm{e}^{-\beta\lvert \bm{v}\rvert ^2} を満たすものとする。直径 dd の剛体球 NN 個の初期分布を f0⊗Nf_0^{\otimes N}(球どうしが重ならない配置に制限して規格化したもの)とする。このとき C,βC,\beta だけで決まる t∗>0t_*>0 があって、0≤t<t∗0\le t<t_* では、Boltzmann–Grad 極限(N→∞N\to\infty、d→0d\to0、Nd2Nd^2 一定)で

fN(1)(t)→f(t),fN(k)(t)→f(t)⊗k(k≥2)f_N^{(1)}(t)\to f(t),\qquad f_N^{(k)}(t)\to f(t)^{\otimes k}\quad(k\ge2)

が成り立つ(収束は相空間の測度 00 の集合を除いた各点で)。ここで ff は f0f_0 を初期値とする Boltzmann 方程式の解である。

2 番目の主張をカオスの伝播という。初期時刻に積の形(独立)だった分布が、時刻 tt でも極限では積の形にとどまる。分子的混沌の仮定は、この意味で仮定でなく定理の帰結になる。t∗t_* の大きさは平均自由時間の数分の 1 程度であり、気体の中の分子がおおむね 1 回も衝突しない程度の時間しかカバーしない。

証明の骨格と難所

証明は次の 3 段からなる。

  1. 級数展開。BBGKY 階層を Duhamel の公式(線型方程式の定数変化法の公式。f˙=Af+g\dot f=Af+g の解を f(t)=etAf(0)+∫0te(t−s)Ag(s) dsf(t)=\mathrm{e}^{tA}f(0)+\int_0^t\mathrm{e}^{(t-s)A}g(s)\,\mathrm{d} s と書く)で積分形に直し、反復代入すると、fN(1)(t)f_N^{(1)}(t) は「時刻 00 から tt までの間に jj 回の衝突で新しい粒子が加わる」項の和 ∑j≥0(⋯ )\sum_{j\ge0}(\cdots) に展開される。Boltzmann 階層についても同じ形の級数がある。
  2. 級数の一様収束。第 jj 項を絶対値で評価すると、衝突演算子 Ck,k+1C_{k,k+1} が速度 ∣v∣\lvert \bm{v}\rvert の重みを 1 つ増やす(失う)ため、毎回 Gauss の重み e−β∣v∣2\mathrm{e}^{-\beta\lvert \bm{v}\rvert ^2} の β\beta を少しずつ削りながら評価することになり、第 jj 項はおよそ (Ct)j(Ct)^j で抑えられる。t<t∗=1/Ct<t_*=1/C なら級数は NN によらず一様に収束する。
  3. 項ごとの収束。各項は、新しい粒子が加わっていく「衝突の木」に沿った軌道の積分である。粒子系の軌道と Boltzmann 階層の軌道が食い違うのは、すでに衝突した粒子どうしが再び衝突する(再衝突)か、加わった粒子が初期時刻に既存の粒子と重なる場合に限られ、そのような配置の測度は d→0d\to0 で 00 になることを示す。

難所は第 2 段である。第 2 段の評価は獲得項と損失項を別々に絶対値で抑えており、両者の打ち消し合いや、保存則・H 定理のような構造をまったく使っていない。このため評価は時間とともに幾何級数的に悪化し、短時間を越えられない。Boltzmann 方程式自体の解は長時間存在する場合でも、この方法では導出の時間を延ばせない。第 3 段の再衝突の評価は技術的に長いが、短時間に限れば測度論的な評価で処理できる。

長時間への延長

2024 年から 2025 年にかけて、Deng・Hani・Ma が arXiv に公開した一連のプレプリントは、級数の項を組合せ論的に組み替えて打ち消し合いを取り込むことで、Boltzmann 方程式の解が存在する限りの長時間にわたって剛体球系から Boltzmann 方程式を導出し、さらに既知の流体力学極限と接続して流体方程式まで到達したと主張している。本書の執筆時点で、この結果が査読を経て確定した状況を本書は確認していない。結果が確定すれば、図 8.1 の Newton 力学から流体方程式までの矢印が(適切な状況設定のもとで)一本につながることになる。

不可逆性の起源

Newton の運動方程式は時間反転について対称である。一方、Boltzmann 方程式は H 定理によって時間の向きを持つ。可逆な力学から不可逆な方程式が出てくることは矛盾に見え、Boltzmann の時代から議論の的になった(Loschmidt の反論)。

Lanford の定理はこの点について次のように答える。定理が保証するのは、初期分布が積の形(独立)であるときに、時刻 t>0t>0 で極限が Boltzmann 方程式に従うことである。時刻 tt の状態で全粒子の速度を反転したものを新たな初期状態にすると、Newton 力学に従えば時間を逆にたどって元に戻る(HH が増える)。しかし速度を反転した状態は、衝突直後の粒子対の間に相関を持っており、積の形をしていない。定理の仮定を満たさないので、Boltzmann 方程式に従う必要はない。つまり Boltzmann 方程式は「典型的な」(積の形の確率分布から見てほとんどすべての)初期条件に対して成り立つ極限であり、不可逆性は方程式でなく初期条件の選び方に由来する。

極限は時間反転と可換でない

Boltzmann–Grad 極限の操作と時間反転の操作は可換でない。極限を取ってから時間を反転すると、Boltzmann 方程式の時間反転(HH が増える方程式)になるが、これは粒子系の極限として現れない。物理の近似が元の方程式の対称性を壊すのは珍しくなく、壊れ方は「どの初期条件の集合で極限を取ったか」で決まる。第6章の平均場極限が時間反転対称性を保つ(Vlasov 方程式は可逆)のと対照的である。

例題 8.4 級数の評価は解の存在より先に壊れる ★★

オリジナルLanford の証明の第 2 段の事情を、数の列の階層で見る。ak(t)a_k(t)(k=1,2,…k=1,2,\dots)が

dakdt=−k ak+1,ak(0)=1\frac{\mathrm{d} a_k}{\mathrm{d} t}=-k\,a_{k+1},\qquad a_k(0)=1

を満たすとする(kk 番目の量の発展に k+1k+1 番目の量が現れる点が BBGKY 階層と同じ形である)。

  1. ak(t)=a(t)ka_k(t)=a(t)^k の形を仮定すると aa は a˙=−a2\dot a=-a^2 を満たすことを示し、ak(t)=(1+t)−ka_k(t)=(1+t)^{-k} がすべての t≥0t\ge0 で解であることを確かめよ。
  2. 積分形 ak(t)=1−k∫0tak+1(s) dsa_k(t)=1-k\int_0^ta_{k+1}(s)\,\mathrm{d} s を反復代入して、a1(t)a_1(t) を tt の冪級数として求めよ。
  3. 各項を絶対値で評価したときの級数が収束する範囲を求め、t=12t=\frac12 で真の値と評価の値を比べよ。
解答

(1) ak=aka_k=a^k を代入すると左辺は kak−1a˙ka^{k-1}\dot a、右辺は −kak+1-ka^{k+1} なので、a˙=−a2\dot a=-a^2 と同値である。a(0)=1a(0)=1 の解は a=1/(1+t)a=1/(1+t) で、すべての t≥0t\ge0 で有限である。このとき ak=(1+t)−ka_k=(1+t)^{-k} で、直接微分しても ddt(1+t)−k=−k(1+t)−k−1=−kak+1\frac{\mathrm{d} }{\mathrm{d} t}(1+t)^{-k}=-k(1+t)^{-k-1}=-ka_{k+1} となる。積の形 ak=aka_k=a^k が保たれることがカオスの伝播に、a˙=−a2\dot a=-a^2 が Boltzmann 方程式に当たる。

(2) jj 回代入すると、a1a_1 に現れる第 jj 項は係数 (−1)j⋅1⋅2⋯j=(−1)jj!(-1)^j\cdot1\cdot2\cdots j=(-1)^jj! と、時間の jj 重積分 ∫0t∫0s1⋯∫0sj−1dsj⋯ds1=tj/j!\int_0^t\int_0^{s_1}\cdots\int_0^{s_{j-1}}\mathrm{d} s_j\cdots\mathrm{d} s_1=t^j/j! の積になる。したがって

a1(t)=∑j≥0(−1)j j! tjj!=∑j≥0(−t)j.a_1(t)=\sum_{j\ge0}(-1)^j\,j!\,\frac{t^j}{j!}=\sum_{j\ge0}(-t)^j.

kk 番目の量が k+1k+1 番目の量を呼ぶたびに係数 kk が掛かり、それが 1/j!1/j! を打ち消している。

(3) 各項の絶対値は tjt^j なので、評価の級数 ∑tj=1/(1−t)\sum t^j=1/(1-t) は t<1t<1 でしか収束しない。t=12t=\frac12 では真の値が a1=2/3a_1=2/3、絶対値の評価が 22 である。真の解 1/(1+t)1/(1+t) はすべての時刻で存在するのに、符号の打ち消し合いを捨てた評価は t=1t=1 で発散する。級数そのもの ∑(−t)j\sum(-t)^j も t≥1t\ge1 では収束しないので、級数展開という方法自体がこの時刻を越えられない。Lanford の定理の t∗t_* が Boltzmann 方程式の解の存在時間でなく、この種の評価の収束半径で決まっていることが、この小さな例で見える(検算スクリプト)。

ほかの巻では

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『虎の巻』の記事一覧へ