統計熱力学入門 第1章 §1.1 ★★★ 統計・熱力学を貫く数学的背骨。共役変数の付け替えという発想が熱力学と統計力学のあらゆる場面で使える。
動機:導関数で関数を符号化する
★★ 傾き p = f ′ ( x ) p=f'(x) p = f ′ ( x ) で関数を符号化する動機。S ↔ T S\leftrightarrow T S ↔ T 付け替えの原型。
ある凸関数 f ( x ) f(x) f ( x ) を考える。これを「点ごとの値 ( x , f ( x ) ) (x,f(x)) ( x , f ( x )) の集まり」として記述する代わりに、「各点での接線 の集まり」として記述したい、というのがルジャンドル変換の根本動機である。接線の傾きを p = f ′ ( x ) p=f'(x) p = f ′ ( x ) とすると、f f f が(狭義)凸なら p p p と x x x は一対一に対応する。したがって x x x の代わりに p p p を独立変数に選び替えても情報は失われない。
なぜこんな付け替えが嬉しいのか。物理では、しばしば「自然な独立変数」と「制御しやすい独立変数」が一致しない。たとえばエネルギー U ( S , … ) U(S,\dots) U ( S , … ) はエントロピー S S S の関数として書くのが自然だが、実験では S S S より温度 T = ∂ U / ∂ S T=\partial U/\partial S T = ∂ U / ∂ S の方が制御しやすい。S ↔ T S\leftrightarrow T S ↔ T の付け替えこそが、まさにこの傾きと独立変数の交換であり、ルジャンドル変換が熱力学の言語そのものになる理由である。
なぜ「ただ代入する」ではいけないのか
独立変数を x x x から p p p に替えたいだけなら、p = f ′ ( x ) p=f'(x) p = f ′ ( x ) を x x x について解いて x = x ( p ) x=x(p) x = x ( p ) を f f f に代入し、f ~ ( p ) ≡ f ( x ( p ) ) \tilde f(p)\equiv f(x(p)) f ~ ( p ) ≡ f ( x ( p )) とすればよさそうに見える。ルジャンドル変換がわざわざ p x px p x を引くのは何のためか、という疑問はここで必ず生じる。答えは「素朴な代入では情報が落ちるから」であり、例を1つ計算すれば即座に分かる。
f ( x ) = x 2 f(x)=x^2 f ( x ) = x 2 と、それを a a a だけ平行移動した g ( x ) = ( x − a ) 2 g(x)=(x-a)^2 g ( x ) = ( x − a ) 2 を比べる。前者は p = 2 x p=2x p = 2 x より x = p / 2 x=p/2 x = p /2 、後者は p = 2 ( x − a ) p=2(x-a) p = 2 ( x − a ) より x = a + p / 2 x=a+p/2 x = a + p /2 である。素朴な代入では
f ~ ( p ) = p 2 4 , g ~ ( p ) = p 2 4 \tilde f(p) = \frac{p^2}{4},\qquad \tilde g(p) = \frac{p^2}{4} f ~ ( p ) = 4 p 2 , g ~ ( p ) = 4 p 2
となり、両者は完全に一致してしまう 。f ~ \tilde f f ~ だけを渡されても、元の放物線がどこにあったか(a a a が何だったか)はもはや復元できない。傾きの情報だけでは、グラフを x x x 方向に平行移動する自由度が失われるのである。一方ルジャンドル変換では
f ∗ ( p ) = p 2 4 , g ∗ ( p ) = a p + p 2 4 f^\ast(p) = \frac{p^2}{4},\qquad g^\ast(p) = a p + \frac{p^2}{4} f ∗ ( p ) = 4 p 2 , g ∗ ( p ) = a p + 4 p 2
となって a a a がちゃんと残る。p x px p x を引くという一手は、装飾でも形式的な約束でもなく、変数の付け替えで情報を落とさないための最小限の補正 である。ここがルジャンドル変換の存在理由であり、f ∗ ∗ = f f^{\ast\ast}=f f ∗∗ = f (2回変換すると元に戻る)が成り立つのも、情報が落ちていないことの言い換えにほかならない。
接線の切片による定義
★★★ f ∗ ( p ) = sup x [ p x − f ( x ) ] f^\ast(p)=\sup_x[px-f(x)] f ∗ ( p ) = sup x [ p x − f ( x )] の定義。ルジャンドル変換の核。
傾き p p p の接線が y y y 軸を切る切片を考える。点 x x x における接線は
y = f ( x ) + p ( x ′ − x ) , p = f ′ ( x ) y = f(x) + p\,(x'-x), \qquad p=f'(x) y = f ( x ) + p ( x ′ − x ) , p = f ′ ( x )
であり(x ′ x' x ′ が直線の変数)、x ′ = 0 x'=0 x ′ = 0 とおくと切片は f ( x ) − p x = − ( p x − f ( x ) ) f(x)-p x = -(px-f(x)) f ( x ) − p x = − ( p x − f ( x )) である。そこで
g ( p ) ≡ p x − f ( x ) ∣ x = x ( p ) , p = f ′ ( x ) g(p) \equiv p x - f(x)\Big|_{x=x(p)}, \qquad p=f'(x) g ( p ) ≡ p x − f ( x ) x = x ( p ) , p = f ′ ( x )
と定義する(切片の符号を反転したもの)。これがルジャンドル変換 f ∗ = g f^\ast=g f ∗ = g である。
定義 (title=★★★ 定義:ルジャンドル変換(凸関数の場合))
微分可能な凸関数 f : R → R f:\mathbb{R}\to\mathbb{R} f : R → R に対し、そのルジャンドル変換 f ∗ f^\ast f ∗ を
f ∗ ( p ) = sup x [ p x − f ( x ) ] \boxed{\,f^\ast(p) \;=\; \sup_{x}\;\bigl[\,p\,x - f(x)\,\bigr]\,} f ∗ ( p ) = x sup [ p x − f ( x ) ]
で定義する。上限を与える x = x ( p ) x=x(p) x = x ( p ) は停留条件
∂ ∂ x [ p x − f ( x ) ] = 0 ⟺ p = f ′ ( x ) \frac{\partial }{\partial x}\bigl[px-f(x)\bigr]=0 \;\Longleftrightarrow\; p=f'(x) ∂ x ∂ [ p x − f ( x ) ] = 0 ⟺ p = f ′ ( x )
で決まる(f f f が凸ゆえ停留点は最大点)。このとき
f ∗ ( p ) = p x ( p ) − f ( x ( p ) ) . f^\ast(p) = p\,x(p) - f\bigl(x(p)\bigr). f ∗ ( p ) = p x ( p ) − f ( x ( p ) ) .
sup \sup sup を用いた定義の利点は、f f f が微分可能でなくても、凸でありさえすれば(折れ点や直線部分があっても)意味を持つことである。この形は大偏差理論・統計力学で本質的になる。
幾何学的意味:接線族の包絡線
★ 点と直線の双対性・包絡線像。理解の補強。
ルジャンドル変換は「点と直線の双対性」と理解するのが最も見通しがよい。凸関数 f f f のグラフは、それを下から支える接線たちの包絡線(envelope) として完全に再現できる。各傾き p p p に対し、対応する接線の切片が − f ∗ ( p ) -f^\ast(p) − f ∗ ( p ) で与えられる。つまり f f f の情報は「点集合 { ( x , f ( x ) ) } \{(x,f(x))\} {( x , f ( x ))} 」と「接線集合 { ( p , − f ∗ ( p ) ) } \{(p,-f^\ast(p))\} {( p , − f ∗ ( p ))} 」のどちらでも等価に表現でき、両者は互いの双対である。
基本的性質
★★ 共役対応・対合性・凸性保存・( f ∗ ) ′ ′ = 1 / f ′ ′ (f^\ast)''=1/f'' ( f ∗ ) ′′ = 1/ f ′′ 。応答関数の関係に直結。
title=★★ ルジャンドル変換の主要な性質
共役変数の対応 :p = f ′ ( x ) p=f'(x) p = f ′ ( x ) かつ x = ( f ∗ ) ′ ( p ) x=(f^\ast)'(p) x = ( f ∗ ) ′ ( p ) 。すなわち変換は独立変数 x x x とその共役量 p p p (傾き)を入れ替える。
対合性(involution) :f f f が凸ならば ( f ∗ ) ∗ = f (f^\ast)^\ast=f ( f ∗ ) ∗ = f 。二回変換すると元に戻る。
凸性の保存 :f ∗ f^\ast f ∗ もまた(下半連続な)凸関数である。
Young の不等式 :任意の x , p x,p x , p に対し f ( x ) + f ∗ ( p ) ≥ p x f(x)+f^\ast(p)\ge px f ( x ) + f ∗ ( p ) ≥ p x 。等号は p = f ′ ( x ) p=f'(x) p = f ′ ( x ) のとき。
スケーリング・平行移動 :g ( x ) = f ( x ) + c g(x)=f(x)+c g ( x ) = f ( x ) + c なら g ∗ ( p ) = f ∗ ( p ) − c g^\ast(p)=f^\ast(p)-c g ∗ ( p ) = f ∗ ( p ) − c ;g ( x ) = f ( x − a ) g(x)=f(x-a) g ( x ) = f ( x − a ) なら g ∗ ( p ) = f ∗ ( p ) + p a g^\ast(p)=f^\ast(p)+pa g ∗ ( p ) = f ∗ ( p ) + p a 。
性質1の証明(x = ( f ∗ ) ′ ( p ) x=(f^\ast)'(p) x = ( f ∗ ) ′ ( p ) )
f ∗ ( p ) = p x ( p ) − f ( x ( p ) ) f^\ast(p)=p\,x(p)-f(x(p)) f ∗ ( p ) = p x ( p ) − f ( x ( p )) を p p p で微分すると、連鎖律で
d f ∗ d p = x ( p ) + p d x d p − f ′ ( x ) d x d p = x ( p ) + ( p − f ′ ( x ) ) ⏟ = 0 d x d p = x ( p ) . \frac{\mathrm{d} f^\ast}{\mathrm{d} p}
= x(p) + p\,\frac{\mathrm{d} x}{\mathrm{d} p} - f'(x)\,\frac{\mathrm{d} x}{\mathrm{d} p}
= x(p) + \underbrace{\bigl(p-f'(x)\bigr)}_{=0}\frac{\mathrm{d} x}{\mathrm{d} p}
= x(p). d p d f ∗ = x ( p ) + p d p d x − f ′ ( x ) d p d x = x ( p ) + = 0 ( p − f ′ ( x ) ) d p d x = x ( p ) .
括弧内は停留条件 p = f ′ ( x ) p=f'(x) p = f ′ ( x ) により消える(これは「包絡線定理」の最も簡単な例である)。したがって f f f の傾きが f ∗ f^\ast f ∗ の独立変数に、f ∗ f^\ast f ∗ の傾きが f f f の独立変数になる、という美しい双対性が成り立つ。□ \square □
凸性と二階微分の関係
p = f ′ ( x ) p=f'(x) p = f ′ ( x ) より d p / d x = f ′ ′ ( x ) \mathrm{d} p/\mathrm{d} x=f''(x) d p / d x = f ′′ ( x ) 、かつ x = ( f ∗ ) ′ ( p ) x=(f^\ast)'(p) x = ( f ∗ ) ′ ( p ) より d x / d p = ( f ∗ ) ′ ′ ( p ) \mathrm{d} x/\mathrm{d} p=(f^\ast)''(p) d x / d p = ( f ∗ ) ′′ ( p ) 。両者は逆数の関係にあるから
( f ∗ ) ′ ′ ( p ) = 1 f ′ ′ ( x ) . (f^\ast)''(p) = \frac{1}{f''(x)}. ( f ∗ ) ′′ ( p ) = f ′′ ( x ) 1 .
f f f が狭義凸(f ′ ′ > 0 f''>0 f ′′ > 0 )なら f ∗ f^\ast f ∗ も狭義凸(( f ∗ ) ′ ′ > 0 (f^\ast)''>0 ( f ∗ ) ′′ > 0 )である。この「曲率の逆数」関係は、後で熱容量や圧縮率といった応答関数(二階微分量)どうしの関係に直結する。
多変数の場合と部分ルジャンドル変換
★★ 一部の変数だけ付け替える部分変換。熱力学ポテンシャル体系の直接の道具。
多変数 f ( x 1 , … , x n ) f(x_1,\dots,x_n) f ( x 1 , … , x n ) では、変数の一部だけ を変換する「部分ルジャンドル変換」が重要になる。たとえば x 1 x_1 x 1 のみを共役変数 p 1 = ∂ f / ∂ x 1 p_1=\partial f/\partial x_1 p 1 = ∂ f / ∂ x 1 に付け替えるとき、
g ( p 1 , x 2 , … , x n ) = sup x 1 [ p 1 x 1 − f ( x 1 , x 2 , … , x n ) ] g(p_1,x_2,\dots,x_n) = \sup_{x_1}\bigl[\,p_1 x_1 - f(x_1,x_2,\dots,x_n)\,\bigr] g ( p 1 , x 2 , … , x n ) = x 1 sup [ p 1 x 1 − f ( x 1 , x 2 , … , x n ) ]
とする。変換しなかった変数 x 2 , … x_2,\dots x 2 , … はそのまま残る。このとき微分の規則は
d g = x 1 d p 1 − ∑ k ≥ 2 ∂ f ∂ x k d x k \mathrm{d} g = x_1\,\mathrm{d} p_1 - \sum_{k\ge 2}\frac{\partial f}{\partial x_k}\,\mathrm{d} x_k d g = x 1 d p 1 − k ≥ 2 ∑ ∂ x k ∂ f d x k
となり、x 1 x_1 x 1 と p 1 p_1 p 1 の役割が入れ替わる一方、他の変数の共役量は符号付きで引き継がれる。熱力学ポテンシャルの体系は、まさにこの部分ルジャンドル変換を S , V , N S,V,N S , V , N の各方向に対して選択的に適用したものに他ならない。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com