この問いに完全な答えを与えたのが、エミー・ネーターが 1918 年にゲッティンゲンで発表した論文「変分問題の不変量」です。当時、ヒルベルトとクラインは一般相対性理論におけるエネルギー保存の奇妙な振る舞い(重力場のエネルギーが局所的にうまく定義できない)に悩んでいました。ネーターはこの問題を、変分原理の不変性という一般的な枠組みで整理し、二つの定理を証明しました。第一定理は、有限個のパラメータを持つ連続対称性が保存則を与えることを述べます。第二定理は、任意関数を含む対称性(ゲージ対称性)が保存則ではなく方程式の間の恒等式を与えることを述べ、一般相対論のエネルギー問題の正体を明らかにしました。
以下、系の配位は n n n 個の一般化座標 q = ( q 1 , … , q n ) q = (q^1, \ldots, q^n) q = ( q 1 , … , q n ) で指定され、ラグランジアン L ( q , q ˙ , t ) L(q, \dot q, t) L ( q , q ˙ , t ) は C 2 C^2 C 2 級とします。時刻 t 1 t_1 t 1 から t 2 t_2 t 2 までの経路 q ( ⋅ ) q(\cdot) q ( ⋅ ) に対する作用は
S [ q ] = ∫ t 1 t 2 L ( q ( t ) , q ˙ ( t ) , t ) d t S[q] = \int_{t_1}^{t_2} L\big(q(t), \dot q(t), t\big)\, dt S [ q ] = ∫ t 1 t 2 L ( q ( t ) , q ˙ ( t ) , t ) d t で、端点を固定した変分が消える条件からオイラー–ラグランジュ方程式
d d t ∂ L ∂ q ˙ i − ∂ L ∂ q i = 0 ( i = 1 , … , n ) \frac{d}{dt}\frac{\partial L}{\partial \dot q^i} - \frac{\partial L}{\partial q^i} = 0
\qquad (i = 1, \ldots, n) d t d ∂ q ˙ i ∂ L − ∂ q i ∂ L = 0 ( i = 1 , … , n ) が出ます(オイラー・ラグランジュ方程式(Theorem 3.4)[Lagrangian Mechanics] )。導出は ラグランジュ形式の力学 を参照してください。以下、この方程式の解を単に解 と呼び、「解に沿って」という言い方をします。
一般化運動量を
p i ( q , q ˙ , t ) = ∂ L ∂ q ˙ i ( q , q ˙ , t ) p_i(q, \dot q, t) = \frac{\partial L}{\partial \dot q^i}(q, \dot q, t) p i ( q , q ˙ , t ) = ∂ q ˙ i ∂ L ( q , q ˙ , t ) で定義します。この記号を使うとオイラー–ラグランジュ方程式は p ˙ i = ∂ L / ∂ q i \dot p_i = \partial L / \partial q^i p ˙ i = ∂ L / ∂ q i と書けます。証明の中で何度もこの形を使います。
Definition 2.1 (保存量 )
関数 I ( q , q ˙ , t ) I(q, \dot q, t) I ( q , q ˙ , t ) が保存量 (第一積分)であるとは、オイラー–ラグランジュ方程式の任意の解 q ( t ) q(t) q ( t ) に対して
d d t I ( q ( t ) , q ˙ ( t ) , t ) = 0 \frac{d}{dt}\, I\big(q(t), \dot q(t), t\big) = 0 d t d I ( q ( t ) , q ˙ ( t ) , t ) = 0 が成り立つことをいいます。すなわち I I I の値は解ごとに決まる定数です。
もっとも簡単な保存量は、座標が一つラグランジアンに現れない場合に得られます。
Definition 2.2 (循環座標 )
∂ L / ∂ q k = 0 \partial L / \partial q^k = 0 ∂ L / ∂ q k = 0 がすべての ( q , q ˙ , t ) (q, \dot q, t) ( q , q ˙ , t ) で成り立つとき、q k q^k q k を循環座標 (無視できる座標)といいます。
Proof(Proposition 2.3) 解に沿ってオイラー–ラグランジュ方程式の i = k i = k i = k 成分を書くと
d p k d t = ∂ L ∂ q k = 0 \frac{d p_k}{dt} = \frac{\partial L}{\partial q^k} = 0 d t d p k = ∂ q k ∂ L = 0 です。最後の等号は Definition 2.2 の仮定によります。したがって p k p_k p k は時間によらず一定です。
∎ この命題は強力ですが、使うには「うまい座標」を見つけなければなりません。たとえば平面上の 2 粒子系で全運動量が保存することを示すには、重心座標を導入して初めて循環座標が現れます。座標の取り方に依存せず、変換そのものを主役にして述べ直したものがネーターの定理です。
「対称性がある」という言葉を、計算できる形に翻訳します。鍵は、変換を一つずつ考えるのではなく、恒等変換から連続的につながった族 として考えることです。
Definition 3.1 (一径数変換族と生成子 )
配位空間の開集合上で定義された写像の族 Q i ( s ; q , t ) Q^i(s; q, t) Q i ( s ; q , t ) (i = 1 , … , n i = 1, \ldots, n i = 1 , … , n 、s s s は 0 0 0 の近くの実数)が、次を満たすとき一径数変換族 といいます。
Q i Q^i Q i は ( s , q , t ) (s, q, t) ( s , q , t ) について C 2 C^2 C 2 級である。
Q i ( 0 ; q , t ) = q i Q^i(0; q, t) = q^i Q i ( 0 ; q , t ) = q i (s = 0 s = 0 s = 0 で恒等変換)。
このとき
ζ i ( q , t ) = ∂ Q i ∂ s ( 0 ; q , t ) \zeta^i(q, t) = \frac{\partial Q^i}{\partial s}(0; q, t) ζ i ( q , t ) = ∂ s ∂ Q i ( 0 ; q , t ) を、この族の生成子 (無限小変換)と呼びます。
生成子は「s s s を少しだけ動かしたときに点がどの向きにどれだけ動くか」を表すベクトル場です。s s s が小さいとき Q i ( s ; q , t ) ≈ q i + s ζ i ( q , t ) Q^i(s; q, t) \approx q^i + s\,\zeta^i(q,t) Q i ( s ; q , t ) ≈ q i + s ζ i ( q , t ) が成り立ちます。
次が対称性の定義です。ここであえて弱い条件 を採用します。ラグランジアンがぴったり不変である必要はなく、時間に関する完全微分だけずれることを許します。理由は後で見ます(Example 4.6 と Exercise 7.2 )。完全微分の差はオイラー–ラグランジュ方程式を変えないので、これは物理的に自然な緩め方です。
Definition 3.2 (ラグランジアンの準対称性 )
一径数変換族 Q ( s ; q , t ) Q(s; q, t) Q ( s ; q , t ) (Definition 3.1 )が L L L の準対称性 であるとは、ある C 1 C^1 C 1 級関数 F ( q , t ) F(q, t) F ( q , t ) が存在して、任意の C 2 C^2 C 2 級曲線 q ( ⋅ ) q(\cdot) q ( ⋅ ) と任意の時刻 t t t について
∂ ∂ s ∣ s = 0 L ( Q ( s ; q ( t ) , t ) , d d t Q ( s ; q ( t ) , t ) , t ) = d d t F ( q ( t ) , t ) \frac{\partial}{\partial s}\bigg|_{s=0}
L\!\left( Q\big(s; q(t), t\big),\ \frac{d}{dt}Q\big(s; q(t), t\big),\ t \right)
= \frac{d}{dt} F\big(q(t), t\big) ∂ s ∂ s = 0 L ( Q ( s ; q ( t ) , t ) , d t d Q ( s ; q ( t ) , t ) , t ) = d t d F ( q ( t ) , t ) が成り立つことをいいます。F ≡ 0 F \equiv 0 F ≡ 0 と取れるとき、単に対称性 (厳密な不変性)といいます。
まず時間を変えない変換について定理を述べます。時間も動かす一般形は次節で扱います。
Theorem 4.1 (ネーターの定理(点変換版) )
L ( q , q ˙ , t ) L(q, \dot q, t) L ( q , q ˙ , t ) を C 2 C^2 C 2 級のラグランジアンとし、一径数変換族 Q ( s ; q , t ) Q(s; q, t) Q ( s ; q , t ) が生成子 ζ i ( q , t ) \zeta^i(q,t) ζ i ( q , t ) を持ち、関数 F ( q , t ) F(q,t) F ( q , t ) とともに L L L の準対称性である(Definition 3.2 )とする。このとき
I ( q , q ˙ , t ) = ∑ i = 1 n ∂ L ∂ q ˙ i ( q , q ˙ , t ) ζ i ( q , t ) − F ( q , t ) I(q, \dot q, t) = \sum_{i=1}^{n} \frac{\partial L}{\partial \dot q^i}(q, \dot q, t)\, \zeta^i(q, t) \;-\; F(q, t) I ( q , q ˙ , t ) = i = 1 ∑ n ∂ q ˙ i ∂ L ( q , q ˙ , t ) ζ i ( q , t ) − F ( q , t ) はオイラー–ラグランジュ方程式の保存量である(Definition 2.1 )。
Proof(Theorem 4.1) C 2 C^2 C 2 級曲線 q ( ⋅ ) q(\cdot) q ( ⋅ ) を一つ固定し、q ~ i ( s , t ) = Q i ( s ; q ( t ) , t ) \tilde q^i(s,t) = Q^i(s; q(t), t) q ~ i ( s , t ) = Q i ( s ; q ( t ) , t ) と置きます。Definition 3.2 の左辺を連鎖律で展開します。L L L の第 1 引数は q ~ \tilde q q ~ 、第 2 引数は ∂ t q ~ \partial_t \tilde q ∂ t q ~ なので
∂ ∂ s ∣ s = 0 L ( q ~ , ∂ t q ~ , t ) = ∑ i ∂ L ∂ q i ∂ q ~ i ∂ s ∣ s = 0 + ∑ i ∂ L ∂ q ˙ i ∂ ∂ s ∣ s = 0 ∂ q ~ i ∂ t \frac{\partial}{\partial s}\bigg|_{s=0} L\big(\tilde q, \partial_t \tilde q, t\big)
= \sum_i \frac{\partial L}{\partial q^i}\,\frac{\partial \tilde q^i}{\partial s}\bigg|_{s=0}
+ \sum_i \frac{\partial L}{\partial \dot q^i}\,\frac{\partial}{\partial s}\bigg|_{s=0}\frac{\partial \tilde q^i}{\partial t} ∂ s ∂ s = 0 L ( q ~ , ∂ t q ~ , t ) = i ∑ ∂ q i ∂ L ∂ s ∂ q ~ i s = 0 + i ∑ ∂ q ˙ i ∂ L ∂ s ∂ s = 0 ∂ t ∂ q ~ i となります。ここで ∂ L / ∂ q i \partial L/\partial q^i ∂ L / ∂ q i などはすべて ( q ( t ) , q ˙ ( t ) , t ) (q(t), \dot q(t), t) ( q ( t ) , q ˙ ( t ) , t ) で評価されています(s = 0 s = 0 s = 0 では q ~ = q ( t ) \tilde q = q(t) q ~ = q ( t ) 、∂ t q ~ = q ˙ ( t ) \partial_t \tilde q = \dot q(t) ∂ t q ~ = q ˙ ( t ) だからです)。
第 1 項の ∂ q ~ i / ∂ s ∣ s = 0 \partial \tilde q^i/\partial s|_{s=0} ∂ q ~ i / ∂ s ∣ s = 0 は、Definition 3.1 の生成子の定義そのもので ζ i ( q ( t ) , t ) \zeta^i(q(t), t) ζ i ( q ( t ) , t ) です。第 2 項では s s s 微分と t t t 微分の順序を交換します。Q Q Q が C 2 C^2 C 2 級で q ( ⋅ ) q(\cdot) q ( ⋅ ) も C 2 C^2 C 2 級なので q ~ \tilde q q ~ は ( s , t ) (s,t) ( s , t ) について C 2 C^2 C 2 級であり、シュワルツの定理(偏微分の順序交換(Theorem 7.1)[多変数関数の微分と偏微分] )により
∂ ∂ s ∂ q ~ i ∂ t = ∂ ∂ t ∂ q ~ i ∂ s \frac{\partial}{\partial s}\frac{\partial \tilde q^i}{\partial t}
= \frac{\partial}{\partial t}\frac{\partial \tilde q^i}{\partial s} ∂ s ∂ ∂ t ∂ q ~ i = ∂ t ∂ ∂ s ∂ q ~ i が使えます。s = 0 s = 0 s = 0 と置けば右辺は d d t ζ i ( q ( t ) , t ) \dfrac{d}{dt}\zeta^i(q(t), t) d t d ζ i ( q ( t ) , t ) です。すなわち具体的には
d d t ζ i ( q ( t ) , t ) = ∑ j ∂ ζ i ∂ q j q ˙ j + ∂ ζ i ∂ t \frac{d}{dt}\zeta^i\big(q(t), t\big)
= \sum_j \frac{\partial \zeta^i}{\partial q^j}\,\dot q^j + \frac{\partial \zeta^i}{\partial t} d t d ζ i ( q ( t ) , t ) = j ∑ ∂ q j ∂ ζ i q ˙ j + ∂ t ∂ ζ i です。以上をまとめ、Definition 3.2 の等式に代入すると、任意の C 2 C^2 C 2 級曲線について
∑ i ∂ L ∂ q i ζ i + ∑ i ∂ L ∂ q ˙ i d ζ i d t = d F d t \sum_i \frac{\partial L}{\partial q^i}\,\zeta^i
+ \sum_i \frac{\partial L}{\partial \dot q^i}\,\frac{d \zeta^i}{dt}
= \frac{dF}{dt} i ∑ ∂ q i ∂ L ζ i + i ∑ ∂ q ˙ i ∂ L d t d ζ i = d t d F が成り立ちます。
ここまでは任意の曲線についての恒等式です。ここで初めて 曲線が解であると仮定します。解に沿ってはオイラー–ラグランジュ方程式より ∂ L / ∂ q i = d d t ∂ L ∂ q ˙ i \partial L/\partial q^i = \dfrac{d}{dt}\dfrac{\partial L}{\partial \dot q^i} ∂ L / ∂ q i = d t d ∂ q ˙ i ∂ L なので、左辺は
∑ i ( d d t ∂ L ∂ q ˙ i ) ζ i + ∑ i ∂ L ∂ q ˙ i d ζ i d t = d d t ( ∑ i ∂ L ∂ q ˙ i ζ i ) \sum_i \left( \frac{d}{dt}\frac{\partial L}{\partial \dot q^i}\right)\zeta^i
+ \sum_i \frac{\partial L}{\partial \dot q^i}\,\frac{d\zeta^i}{dt}
= \frac{d}{dt}\left( \sum_i \frac{\partial L}{\partial \dot q^i}\,\zeta^i \right) i ∑ ( d t d ∂ q ˙ i ∂ L ) ζ i + i ∑ ∂ q ˙ i ∂ L d t d ζ i = d t d ( i ∑ ∂ q ˙ i ∂ L ζ i ) と積の微分公式でまとめられます。したがって解に沿って
d d t ( ∑ i ∂ L ∂ q ˙ i ζ i − F ) = 0 \frac{d}{dt}\left( \sum_i \frac{\partial L}{\partial \dot q^i}\,\zeta^i - F \right) = 0 d t d ( i ∑ ∂ q ˙ i ∂ L ζ i − F ) = 0 が成り立ちます。これが主張です。
∎ flowchart LR
A["時間の並進"] --> A2["エネルギー関数 h"]
B["空間の並進"] --> B2["全運動量 P"]
C["空間の回転"] --> C2["全角運動量 J"]
D["ガリレイ変換"] --> D2["重心の等速直線運動"] 対称性と保存量の対応。左の変換に対して、右の量が解に沿って一定に保たれる。 N N N 個の質点からなる系を考えます。質点 a a a の位置を r a \boldsymbol{r}_a r a 、質量を m a m_a m a とし、ラグランジアンを
L = ∑ a = 1 N m a 2 ∣ r ˙ a ∣ 2 − U ( r 1 , … , r N ) L = \sum_{a=1}^{N} \frac{m_a}{2}\,|\dot{\boldsymbol{r}}_a|^2 - U(\boldsymbol{r}_1, \ldots, \boldsymbol{r}_N) L = a = 1 ∑ N 2 m a ∣ r ˙ a ∣ 2 − U ( r 1 , … , r N ) とします。
Example 4.3 (全運動量の保存 )
ポテンシャル U U U が、すべての質点を同じだけ平行移動しても変わらないとします。すなわち任意のベクトル c \boldsymbol{c} c について
U ( r 1 + c , … , r N + c ) = U ( r 1 , … , r N ) U(\boldsymbol{r}_1 + \boldsymbol{c}, \ldots, \boldsymbol{r}_N + \boldsymbol{c}) = U(\boldsymbol{r}_1, \ldots, \boldsymbol{r}_N) U ( r 1 + c , … , r N + c ) = U ( r 1 , … , r N ) が成り立つとします。相互作用が距離 ∣ r a − r b ∣ |\boldsymbol{r}_a - \boldsymbol{r}_b| ∣ r a − r b ∣ だけで書ける系はこの条件を満たします。
単位ベクトル n \boldsymbol{n} n を固定し、変換族 R a ( s ) = r a + s n \boldsymbol{R}_a(s) = \boldsymbol{r}_a + s\,\boldsymbol{n} R a ( s ) = r a + s n を取ります。生成子は ζ a = n \boldsymbol{\zeta}_a = \boldsymbol{n} ζ a = n (a a a によらない定ベクトル)です。速度は R ˙ a = r ˙ a \dot{\boldsymbol{R}}_a = \dot{\boldsymbol{r}}_a R ˙ a = r ˙ a で変わらないので運動エネルギーは不変、ポテンシャルも仮定より不変です。よって
∂ ∂ s ∣ s = 0 L = 0 \frac{\partial}{\partial s}\bigg|_{s=0} L = 0 ∂ s ∂ s = 0 L = 0 となり、F = 0 F = 0 F = 0 の厳密な対称性です(Definition 3.2 )。
一般化運動量は ∂ L / ∂ r ˙ a = m a r ˙ a \partial L / \partial \dot{\boldsymbol{r}}_a = m_a \dot{\boldsymbol{r}}_a ∂ L / ∂ r ˙ a = m a r ˙ a です。Theorem 4.1 の公式に代入すると
I = ∑ a = 1 N m a r ˙ a ⋅ n = ( ∑ a = 1 N m a r ˙ a ) ⋅ n = P ⋅ n I = \sum_{a=1}^{N} m_a \dot{\boldsymbol{r}}_a \cdot \boldsymbol{n}
= \left( \sum_{a=1}^{N} m_a \dot{\boldsymbol{r}}_a \right)\cdot \boldsymbol{n}
= \boldsymbol{P}\cdot\boldsymbol{n} I = a = 1 ∑ N m a r ˙ a ⋅ n = ( a = 1 ∑ N m a r ˙ a ) ⋅ n = P ⋅ n が保存します。n \boldsymbol{n} n は任意に取れるので、P ⋅ e x \boldsymbol{P}\cdot\boldsymbol{e}_x P ⋅ e x 、P ⋅ e y \boldsymbol{P}\cdot\boldsymbol{e}_y P ⋅ e y 、P ⋅ e z \boldsymbol{P}\cdot\boldsymbol{e}_z P ⋅ e z がそれぞれ保存し、結局ベクトル P = ∑ a m a r ˙ a \boldsymbol{P} = \sum_a m_a\dot{\boldsymbol{r}}_a P = ∑ a m a r ˙ a の 3 成分すべてが保存します。
Example 4.4 (全角運動量の保存 )
同じ系で、今度は U U U が原点まわりの回転で不変、すなわち任意の回転行列 R ∈ S O ( 3 ) R \in SO(3) R ∈ S O ( 3 ) について U ( R r 1 , … , R r N ) = U ( r 1 , … , r N ) U(R\boldsymbol{r}_1, \ldots, R\boldsymbol{r}_N) = U(\boldsymbol{r}_1,\ldots,\boldsymbol{r}_N) U ( R r 1 , … , R r N ) = U ( r 1 , … , r N ) とします。相互作用が ∣ r a − r b ∣ |\boldsymbol{r}_a - \boldsymbol{r}_b| ∣ r a − r b ∣ と ∣ r a ∣ |\boldsymbol{r}_a| ∣ r a ∣ だけで書ける系はこれを満たします。
単位ベクトル n \boldsymbol{n} n を軸とし、角 s s s の回転行列を R ( s n ) R(s\boldsymbol{n}) R ( s n ) として、変換族 R a ( s ) = R ( s n ) r a \boldsymbol{R}_a(s) = R(s\boldsymbol{n})\,\boldsymbol{r}_a R a ( s ) = R ( s n ) r a を取ります。回転の生成子はよく知られた外積の形
ζ a = d d s ∣ s = 0 R ( s n ) r a = n × r a \boldsymbol{\zeta}_a = \frac{d}{ds}\bigg|_{s=0} R(s\boldsymbol{n})\,\boldsymbol{r}_a = \boldsymbol{n}\times\boldsymbol{r}_a ζ a = d s d s = 0 R ( s n ) r a = n × r a です。速度は R ˙ a ( s ) = R ( s n ) r ˙ a \dot{\boldsymbol{R}}_a(s) = R(s\boldsymbol{n})\dot{\boldsymbol{r}}_a R ˙ a ( s ) = R ( s n ) r ˙ a であり、回転行列は長さを保つので ∣ R ˙ a ( s ) ∣ = ∣ r ˙ a ∣ |\dot{\boldsymbol{R}}_a(s)| = |\dot{\boldsymbol{r}}_a| ∣ R ˙ a ( s ) ∣ = ∣ r ˙ a ∣ 、したがって運動エネルギーは s s s によらず一定です。ポテンシャルも仮定より一定です。よってこれも F = 0 F = 0 F = 0 の厳密な対称性です。
Theorem 4.1 より
I = ∑ a m a r ˙ a ⋅ ( n × r a ) I = \sum_a m_a \dot{\boldsymbol{r}}_a \cdot (\boldsymbol{n}\times\boldsymbol{r}_a) I = a ∑ m a r ˙ a ⋅ ( n × r a ) が保存します。スカラー三重積の巡回置換 a ⋅ ( b × c ) = b ⋅ ( c × a ) \boldsymbol{a}\cdot(\boldsymbol{b}\times\boldsymbol{c}) = \boldsymbol{b}\cdot(\boldsymbol{c}\times\boldsymbol{a}) a ⋅ ( b × c ) = b ⋅ ( c × a ) を a = r ˙ a \boldsymbol{a} = \dot{\boldsymbol{r}}_a a = r ˙ a 、b = n \boldsymbol{b} = \boldsymbol{n} b = n 、c = r a \boldsymbol{c} = \boldsymbol{r}_a c = r a に適用すると
r ˙ a ⋅ ( n × r a ) = n ⋅ ( r a × r ˙ a ) \dot{\boldsymbol{r}}_a\cdot(\boldsymbol{n}\times\boldsymbol{r}_a) = \boldsymbol{n}\cdot(\boldsymbol{r}_a\times\dot{\boldsymbol{r}}_a) r ˙ a ⋅ ( n × r a ) = n ⋅ ( r a × r ˙ a ) なので
I = n ⋅ ∑ a r a × m a r ˙ a = n ⋅ J I = \boldsymbol{n}\cdot \sum_a \boldsymbol{r}_a \times m_a\dot{\boldsymbol{r}}_a = \boldsymbol{n}\cdot\boldsymbol{J} I = n ⋅ a ∑ r a × m a r ˙ a = n ⋅ J です。n \boldsymbol{n} n は任意なので、全角運動量 J = ∑ a r a × m a r ˙ a \boldsymbol{J} = \sum_a \boldsymbol{r}_a\times m_a\dot{\boldsymbol{r}}_a J = ∑ a r a × m a r ˙ a の 3 成分すべてが保存します。
F = 0 F = 0 F = 0 に限定すると取り逃す保存量があります。
Example 4.6 (一様重力中の鉛直並進 )
一様重力中の質点 L = m 2 ( x ˙ 2 + y ˙ 2 + z ˙ 2 ) − m g z L = \dfrac{m}{2}(\dot x^2 + \dot y^2 + \dot z^2) - mgz L = 2 m ( x ˙ 2 + y ˙ 2 + z ˙ 2 ) − m g z を考えます。鉛直方向の並進 Z ( s ) = z + s Z(s) = z + s Z ( s ) = z + s (x , y x, y x , y は動かさない)を取ると、生成子は ζ = ( 0 , 0 , 1 ) \zeta = (0,0,1) ζ = ( 0 , 0 , 1 ) で、
∂ ∂ s ∣ s = 0 L = − m g \frac{\partial}{\partial s}\bigg|_{s=0} L = -mg ∂ s ∂ s = 0 L = − m g です。0 0 0 ではないので厳密な対称性ではありません。しかし − m g -mg − m g は定数なので、F ( q , t ) = − m g t F(q,t) = -mgt F ( q , t ) = − m g t と置けば d F / d t = − m g dF/dt = -mg d F / d t = − m g となり、Definition 3.2 の準対称性の条件を満たします。
p z = m z ˙ p_z = m\dot z p z = m z ˙ なので Theorem 4.1 より
I = m z ˙ − ( − m g t ) = m z ˙ + m g t I = m\dot z - (-mgt) = m\dot z + mgt I = m z ˙ − ( − m g t ) = m z ˙ + m g t が保存量です。実際に解 z ( t ) = z 0 + v 0 t − 1 2 g t 2 z(t) = z_0 + v_0 t - \tfrac12 g t^2 z ( t ) = z 0 + v 0 t − 2 1 g t 2 で確かめると m z ˙ + m g t = m ( v 0 − g t ) + m g t = m v 0 m\dot z + mgt = m(v_0 - gt) + mgt = mv_0 m z ˙ + m g t = m ( v 0 − g t ) + m g t = m v 0 で、確かに定数です。これはエネルギーとも運動量とも異なる、独立な保存量です(内容としては「鉛直方向の加速度が − g -g − g で一定」という情報を一つの式にまとめたものです)。F = 0 F = 0 F = 0 に固執していたら、この保存量は見えませんでした。
ここまでの定理は時間を動かしませんでした。エネルギー保存則は時間の一様性から来るはずなので、時間も変換する形に拡張する必要があります。準備として、まずエネルギー関数を導入します。
Definition 5.1 (エネルギー関数(ヤコビ積分) )
ラグランジアン L ( q , q ˙ , t ) L(q,\dot q,t) L ( q , q ˙ , t ) に対し
h ( q , q ˙ , t ) = ∑ i = 1 n q ˙ i ∂ L ∂ q ˙ i ( q , q ˙ , t ) − L ( q , q ˙ , t ) h(q, \dot q, t) = \sum_{i=1}^{n} \dot q^i\, \frac{\partial L}{\partial \dot q^i}(q,\dot q,t) \;-\; L(q, \dot q, t) h ( q , q ˙ , t ) = i = 1 ∑ n q ˙ i ∂ q ˙ i ∂ L ( q , q ˙ , t ) − L ( q , q ˙ , t ) をエネルギー関数 (ヤコビ積分)といいます。
Proposition 5.2 (エネルギー関数の時間変化 )
オイラー–ラグランジュ方程式の任意の解 q ( t ) q(t) q ( t ) に沿って
d d t h ( q ( t ) , q ˙ ( t ) , t ) = − ∂ L ∂ t ( q ( t ) , q ˙ ( t ) , t ) \frac{d}{dt}\,h\big(q(t),\dot q(t),t\big) = -\frac{\partial L}{\partial t}\big(q(t),\dot q(t),t\big) d t d h ( q ( t ) , q ˙ ( t ) , t ) = − ∂ t ∂ L ( q ( t ) , q ˙ ( t ) , t ) が成り立つ。特に L L L が t t t に陽に依存しない(∂ L / ∂ t ≡ 0 \partial L/\partial t \equiv 0 ∂ L / ∂ t ≡ 0 )ならば h h h は保存量である。
Proof(Proposition 5.2) p i = ∂ L / ∂ q ˙ i p_i = \partial L/\partial \dot q^i p i = ∂ L / ∂ q ˙ i と書きます。解に沿って積の微分公式を使うと
d h d t = ∑ i ( q ¨ i p i + q ˙ i p ˙ i ) − d L d t \frac{dh}{dt} = \sum_i \big( \ddot q^i p_i + \dot q^i \dot p_i \big) - \frac{dL}{dt} d t d h = i ∑ ( q ¨ i p i + q ˙ i p ˙ i ) − d t d L です。一方 L L L の全微分は連鎖律より
d L d t = ∑ i ∂ L ∂ q i q ˙ i + ∑ i ∂ L ∂ q ˙ i q ¨ i + ∂ L ∂ t = ∑ i ∂ L ∂ q i q ˙ i + ∑ i p i q ¨ i + ∂ L ∂ t \frac{dL}{dt} = \sum_i \frac{\partial L}{\partial q^i}\dot q^i + \sum_i \frac{\partial L}{\partial \dot q^i}\ddot q^i + \frac{\partial L}{\partial t}
= \sum_i \frac{\partial L}{\partial q^i}\dot q^i + \sum_i p_i \ddot q^i + \frac{\partial L}{\partial t} d t d L = i ∑ ∂ q i ∂ L q ˙ i + i ∑ ∂ q ˙ i ∂ L q ¨ i + ∂ t ∂ L = i ∑ ∂ q i ∂ L q ˙ i + i ∑ p i q ¨ i + ∂ t ∂ L です。差を取ると ∑ i p i q ¨ i \sum_i p_i\ddot q^i ∑ i p i q ¨ i の項が相殺して
d h d t = ∑ i q ˙ i ( p ˙ i − ∂ L ∂ q i ) − ∂ L ∂ t \frac{dh}{dt} = \sum_i \dot q^i \left( \dot p_i - \frac{\partial L}{\partial q^i} \right) - \frac{\partial L}{\partial t} d t d h = i ∑ q ˙ i ( p ˙ i − ∂ q i ∂ L ) − ∂ t ∂ L となります。ここで解に沿ってはオイラー–ラグランジュ方程式より丸括弧の中が 0 0 0 なので、d h / d t = − ∂ L / ∂ t dh/dt = -\partial L/\partial t d h / d t = − ∂ L / ∂ t を得ます。
∎ この命題は単独でも使えますが、対称性の言葉で言い直しておくと全体の構造が見えます。そのために定理を拡張します。
Theorem 5.3 (ネーターの定理(時間変換を含む形) )
L ( q , q ˙ , t ) L(q,\dot q,t) L ( q , q ˙ , t ) を C 2 C^2 C 2 級とする。C 2 C^2 C 2 級の族 T ( s ; q , t ) T(s;q,t) T ( s ; q , t ) 、Q i ( s ; q , t ) Q^i(s;q,t) Q i ( s ; q , t ) が T ( 0 ; q , t ) = t T(0;q,t) = t T ( 0 ; q , t ) = t 、Q i ( 0 ; q , t ) = q i Q^i(0;q,t) = q^i Q i ( 0 ; q , t ) = q i を満たすとし、生成子を
τ ( q , t ) = ∂ T ∂ s ( 0 ; q , t ) , ζ i ( q , t ) = ∂ Q i ∂ s ( 0 ; q , t ) \tau(q,t) = \frac{\partial T}{\partial s}(0;q,t), \qquad
\zeta^i(q,t) = \frac{\partial Q^i}{\partial s}(0;q,t) τ ( q , t ) = ∂ s ∂ T ( 0 ; q , t ) , ζ i ( q , t ) = ∂ s ∂ Q i ( 0 ; q , t ) とする。さらに、ある C 1 C^1 C 1 級関数 F ( q , t ) F(q,t) F ( q , t ) が存在して、任意の C 2 C^2 C 2 級曲線 q ( ⋅ ) q(\cdot) q ( ⋅ ) について
∂ ∂ s ∣ s = 0 [ L ( Q , d Q / d t d T / d t , T ) d T d t ] = d d t F ( q ( t ) , t ) \frac{\partial}{\partial s}\bigg|_{s=0}
\left[ L\!\left( Q,\ \frac{dQ/dt}{dT/dt},\ T \right) \frac{dT}{dt} \right]
= \frac{d}{dt}F\big(q(t),t\big) ∂ s ∂ s = 0 [ L ( Q , d T / d t d Q / d t , T ) d t d T ] = d t d F ( q ( t ) , t ) が成り立つとする(Q , T Q, T Q , T の引数は ( s ; q ( t ) , t ) (s; q(t), t) ( s ; q ( t ) , t ) 、d / d t d/dt d / d t は曲線に沿う全微分)。このとき
I = ∑ i = 1 n p i ζ i − h τ − F , p i = ∂ L ∂ q ˙ i , h = ∑ i q ˙ i p i − L I = \sum_{i=1}^{n} p_i\,\zeta^i \;-\; h\,\tau \;-\; F,
\qquad p_i = \frac{\partial L}{\partial \dot q^i},\quad
h = \sum_i \dot q^i p_i - L I = i = 1 ∑ n p i ζ i − h τ − F , p i = ∂ q ˙ i ∂ L , h = i ∑ q ˙ i p i − L は保存量である。
不変性の条件が奇妙に見えるかもしれませんが、これは「変換後の時刻で測った作用が変わらない」という条件を書き下したものです。変換後の作用を t ~ = T ( s ; q ( t ) , t ) \tilde t = T(s;q(t),t) t ~ = T ( s ; q ( t ) , t ) で積分変数変換すると d t ~ = ( d T / d t ) d t d\tilde t = (dT/dt)\,dt d t ~ = ( d T / d t ) d t となり、被積分関数が上の角括弧の中身になります。速度は d q ~ / d t ~ = ( d Q / d t ) / ( d T / d t ) d\tilde q/d\tilde t = (dQ/dt)/(dT/dt) d q ~ / d t ~ = ( d Q / d t ) / ( d T / d t ) です。
変換の生成子。もとの経路(実線)上の点が、時間方向に s τ、座標方向に s ζ だけ動いて変換後の経路(破線)に乗る。 Proof(Theorem 5.3) 曲線 q ( ⋅ ) q(\cdot) q ( ⋅ ) を固定し、s = 0 s = 0 s = 0 での値を使って各因子の s s s 微分を計算します。s = 0 s=0 s = 0 では Q = q ( t ) Q = q(t) Q = q ( t ) 、T = t T = t T = t 、d T / d t = 1 dT/dt = 1 d T / d t = 1 、( d Q / d t ) / ( d T / d t ) = q ˙ ( t ) (dQ/dt)/(dT/dt) = \dot q(t) ( d Q / d t ) / ( d T / d t ) = q ˙ ( t ) です。δ \delta δ で ∂ s ∣ s = 0 \partial_s|_{s=0} ∂ s ∣ s = 0 を表すと、生成子の定義から
δ Q i = ζ i , δ T = τ , δ ( d Q i d t ) = d ζ i d t , δ ( d T d t ) = d τ d t \delta Q^i = \zeta^i, \qquad \delta T = \tau, \qquad
\delta\!\left(\frac{dQ^i}{dt}\right) = \frac{d\zeta^i}{dt}, \qquad
\delta\!\left(\frac{dT}{dt}\right) = \frac{d\tau}{dt} δ Q i = ζ i , δ T = τ , δ ( d t d Q i ) = d t d ζ i , δ ( d t d T ) = d t d τ です(Theorem 4.1 の証明と同じくシュワルツの定理で s s s 微分と t t t 微分を交換しました)。商の微分法より
δ ( d Q i / d t d T / d t ) = ( d ζ i / d t ) ⋅ 1 − q ˙ i ( d τ / d t ) 1 2 = d ζ i d t − q ˙ i d τ d t \delta\!\left(\frac{dQ^i/dt}{dT/dt}\right)
= \frac{(d\zeta^i/dt)\cdot 1 - \dot q^i (d\tau/dt)}{1^2}
= \frac{d\zeta^i}{dt} - \dot q^i \frac{d\tau}{dt} δ ( d T / d t d Q i / d t ) = 1 2 ( d ζ i / d t ) ⋅ 1 − q ˙ i ( d τ / d t ) = d t d ζ i − q ˙ i d t d τ です。これらを使って角括弧全体を微分すると、積の微分法と連鎖律から
δ [ L d T d t ] = ∑ i ∂ L ∂ q i ζ i + ∑ i ∂ L ∂ q ˙ i ( d ζ i d t − q ˙ i d τ d t ) + ∂ L ∂ t τ + L d τ d t \delta\!\left[ L \frac{dT}{dt} \right]
= \sum_i \frac{\partial L}{\partial q^i}\zeta^i
+ \sum_i \frac{\partial L}{\partial \dot q^i}\left( \frac{d\zeta^i}{dt} - \dot q^i \frac{d\tau}{dt} \right)
+ \frac{\partial L}{\partial t}\tau
+ L\,\frac{d\tau}{dt} δ [ L d t d T ] = i ∑ ∂ q i ∂ L ζ i + i ∑ ∂ q ˙ i ∂ L ( d t d ζ i − q ˙ i d t d τ ) + ∂ t ∂ L τ + L d t d τ となります。仮定よりこれは d F / d t dF/dt d F / d t に等しく、これは任意の曲線についての恒等式です。
ここで曲線が解であると仮定します。オイラー–ラグランジュ方程式 ∂ L / ∂ q i = p ˙ i \partial L/\partial q^i = \dot p_i ∂ L / ∂ q i = p ˙ i を使うと第 1 項と第 2 項の前半がまとまり、
d d t ( ∑ i p i ζ i ) − ( ∑ i p i q ˙ i − L ) d τ d t + ∂ L ∂ t τ = d F d t \frac{d}{dt}\left( \sum_i p_i \zeta^i \right)
- \left( \sum_i p_i \dot q^i - L \right)\frac{d\tau}{dt}
+ \frac{\partial L}{\partial t}\,\tau
= \frac{dF}{dt} d t d ( i ∑ p i ζ i ) − ( i ∑ p i q ˙ i − L ) d t d τ + ∂ t ∂ L τ = d t d F を得ます。丸括弧は Definition 5.1 のエネルギー関数 h h h そのものです。さらに Proposition 5.2 より解に沿って d h / d t = − ∂ L / ∂ t dh/dt = -\partial L/\partial t d h / d t = − ∂ L / ∂ t なので
d d t ( h τ ) = d h d t τ + h d τ d t = − ∂ L ∂ t τ + h d τ d t \frac{d}{dt}\big( h\,\tau \big)
= \frac{dh}{dt}\tau + h \frac{d\tau}{dt}
= -\frac{\partial L}{\partial t}\tau + h\frac{d\tau}{dt} d t d ( h τ ) = d t d h τ + h d t d τ = − ∂ t ∂ L τ + h d t d τ が成り立ちます。すなわち上の式の第 2 項と第 3 項の和は − d d t ( h τ ) -\dfrac{d}{dt}(h\tau) − d t d ( h τ ) に等しくなります。したがって
d d t ( ∑ i p i ζ i − h τ − F ) = 0 \frac{d}{dt}\left( \sum_i p_i \zeta^i - h\,\tau - F \right) = 0 d t d ( i ∑ p i ζ i − h τ − F ) = 0 となり、主張が示されました。
∎ Corollary 5.4 (時間の一様性とエネルギー保存 )
∂ L / ∂ t ≡ 0 \partial L/\partial t \equiv 0 ∂ L / ∂ t ≡ 0 ならば、時間並進 T ( s ; q , t ) = t + s T(s;q,t) = t + s T ( s ; q , t ) = t + s 、Q i ( s ; q , t ) = q i Q^i(s;q,t) = q^i Q i ( s ; q , t ) = q i は Theorem 5.3 の意味で F = 0 F = 0 F = 0 の対称性であり、対応する保存量はエネルギー関数 h h h (の符号を変えたもの)である。
Proof(Corollary 5.4) この変換では d T / d t = 1 dT/dt = 1 d T / d t = 1 、d Q i / d t = q ˙ i dQ^i/dt = \dot q^i d Q i / d t = q ˙ i なので、角括弧の中身は L ( q ( t ) , q ˙ ( t ) , t + s ) L(q(t), \dot q(t), t+s) L ( q ( t ) , q ˙ ( t ) , t + s ) です。s s s で微分して s = 0 s = 0 s = 0 と置くと ∂ L / ∂ t \partial L/\partial t ∂ L / ∂ t となり、仮定よりこれは 0 0 0 です。0 = d F / d t 0 = dF/dt 0 = d F / d t は F = 0 F = 0 F = 0 で満たされるので、準対称性の条件(Theorem 5.3 )が成り立ちます。生成子は τ = 1 \tau = 1 τ = 1 、ζ i = 0 \zeta^i = 0 ζ i = 0 なので保存量は
I = 0 − h ⋅ 1 − 0 = − h I = 0 - h\cdot 1 - 0 = -h I = 0 − h ⋅ 1 − 0 = − h です。− h -h − h が保存することと h h h が保存することは同値なので、h h h が保存量です。これは Proposition 5.2 の後半と一致します。
∎ 学部の講義では「時間並進対称性からエネルギー T + U T + U T + U が保存する」と言われますが、正確には保存するのは h h h であって T + U T + U T + U ではありません。両者が一致するには条件が要ります。
Lemma 5.5 (オイラーの同次関数定理 )
関数 f ( q ˙ 1 , … , q ˙ n ) f(\dot q^1, \ldots, \dot q^n) f ( q ˙ 1 , … , q ˙ n ) が q ˙ \dot q q ˙ について k k k 次同次、すなわち任意の λ > 0 \lambda > 0 λ > 0 について f ( λ q ˙ ) = λ k f ( q ˙ ) f(\lambda \dot q) = \lambda^k f(\dot q) f ( λ q ˙ ) = λ k f ( q ˙ ) を満たし、かつ C 1 C^1 C 1 級ならば
∑ i = 1 n q ˙ i ∂ f ∂ q ˙ i ( q ˙ ) = k f ( q ˙ ) \sum_{i=1}^{n} \dot q^i \frac{\partial f}{\partial \dot q^i}(\dot q) = k\, f(\dot q) i = 1 ∑ n q ˙ i ∂ q ˙ i ∂ f ( q ˙ ) = k f ( q ˙ ) が成り立つ。
Proof(Lemma 5.5) f ( λ q ˙ ) = λ k f ( q ˙ ) f(\lambda\dot q) = \lambda^k f(\dot q) f ( λ q ˙ ) = λ k f ( q ˙ ) の両辺を λ \lambda λ で微分します。左辺は連鎖律より ∑ i q ˙ i ∂ f ∂ q ˙ i ( λ q ˙ ) \sum_i \dot q^i \dfrac{\partial f}{\partial \dot q^i}(\lambda\dot q) ∑ i q ˙ i ∂ q ˙ i ∂ f ( λ q ˙ ) 、右辺は k λ k − 1 f ( q ˙ ) k\lambda^{k-1}f(\dot q) k λ k − 1 f ( q ˙ ) です。λ = 1 \lambda = 1 λ = 1 と置くと主張の式を得ます。
∎ Proposition 5.6 (エネルギー関数が力学的エネルギーになる条件 )
ラグランジアンが L = T − U L = T - U L = T − U の形に分解され、次の 2 条件を満たすとする。
U = U ( q , t ) U = U(q, t) U = U ( q , t ) は速度 q ˙ \dot q q ˙ に依存しない。
T = T ( q , q ˙ , t ) T = T(q, \dot q, t) T = T ( q , q ˙ , t ) は q ˙ \dot q q ˙ について 2 次同次である。すなわち T ( q , λ q ˙ , t ) = λ 2 T ( q , q ˙ , t ) T(q, \lambda\dot q, t) = \lambda^2 T(q,\dot q,t) T ( q , λ q ˙ , t ) = λ 2 T ( q , q ˙ , t ) 。
このとき h = T + U h = T + U h = T + U である。
Proof(Proposition 5.6) 条件 1 より ∂ L / ∂ q ˙ i = ∂ T / ∂ q ˙ i \partial L/\partial \dot q^i = \partial T/\partial \dot q^i ∂ L / ∂ q ˙ i = ∂ T / ∂ q ˙ i です。条件 2 と Lemma 5.5 (k = 2 k = 2 k = 2 、q q q と t t t は固定して q ˙ \dot q q ˙ の関数とみる)より
∑ i q ˙ i ∂ T ∂ q ˙ i = 2 T \sum_i \dot q^i \frac{\partial T}{\partial \dot q^i} = 2T i ∑ q ˙ i ∂ q ˙ i ∂ T = 2 T です。したがって Definition 5.1 より
h = ∑ i q ˙ i ∂ L ∂ q ˙ i − L = 2 T − ( T − U ) = T + U h = \sum_i \dot q^i \frac{\partial L}{\partial \dot q^i} - L = 2T - (T - U) = T + U h = i ∑ q ˙ i ∂ q ˙ i ∂ L − L = 2 T − ( T − U ) = T + U となります。
∎ 条件 2 は無条件に成り立つものではありません。一般化座標からデカルト座標への変換 r a = r a ( q ) \boldsymbol{r}_a = \boldsymbol{r}_a(q) r a = r a ( q ) が時間に陽に依存しない(スクレロノミック な)場合には、r ˙ a = ∑ i ∂ r a ∂ q i q ˙ i \dot{\boldsymbol{r}}_a = \sum_i \dfrac{\partial \boldsymbol{r}_a}{\partial q^i}\dot q^i r ˙ a = ∑ i ∂ q i ∂ r a q ˙ i なので
T = ∑ a m a 2 ∣ r ˙ a ∣ 2 = 1 2 ∑ i , j a i j ( q ) q ˙ i q ˙ j T = \sum_a \frac{m_a}{2}|\dot{\boldsymbol{r}}_a|^2 = \frac{1}{2}\sum_{i,j} a_{ij}(q)\,\dot q^i \dot q^j T = a ∑ 2 m a ∣ r ˙ a ∣ 2 = 2 1 i , j ∑ a ij ( q ) q ˙ i q ˙ j となり、確かに 2 次同次です。しかし変換が時間に陽に依存する(レオノミック な)場合には q ˙ \dot q q ˙ の 1 次項や 0 次項が現れ、条件 2 が壊れます。次がその典型例です。
Example 5.7 (回転する棒の上を滑るビーズ:h は保存するが T+U は保存しない )
水平面内で原点まわりに一定角速度 ω \omega ω で回転している真っ直ぐな棒を考えます。棒の上を質量 m m m のビーズが摩擦なく滑ります。棒は外部のモーターで駆動されており、ω \omega ω は何が起きても一定に保たれるものとします。
ビーズの原点からの距離を r r r とすると、ラボ系での位置は r = ( r cos ω t , r sin ω t , 0 ) \boldsymbol{r} = (r\cos\omega t,\ r\sin\omega t,\ 0) r = ( r cos ω t , r sin ω t , 0 ) です。これは t t t に陽に依存する(レオノミックな)変換です。速度を計算すると
r ˙ = ( r ˙ cos ω t − r ω sin ω t , r ˙ sin ω t + r ω cos ω t , 0 ) \dot{\boldsymbol{r}} = (\dot r\cos\omega t - r\omega\sin\omega t,\ \dot r \sin\omega t + r\omega\cos\omega t,\ 0) r ˙ = ( r ˙ cos ω t − r ω sin ω t , r ˙ sin ω t + r ω cos ω t , 0 ) なので、三角関数の基本関係を使って
∣ r ˙ ∣ 2 = r ˙ 2 + r 2 ω 2 |\dot{\boldsymbol{r}}|^2 = \dot r^2 + r^2\omega^2 ∣ r ˙ ∣ 2 = r ˙ 2 + r 2 ω 2 です。重力は水平面に垂直で仕事をしないので、ポテンシャルは 0 0 0 と取れます。よって
L = m 2 ( r ˙ 2 + r 2 ω 2 ) . L = \frac{m}{2}\big( \dot r^2 + r^2\omega^2 \big). L = 2 m ( r ˙ 2 + r 2 ω 2 ) . L L L は t t t に陽に依存しないので、Corollary 5.4 より h h h は保存します。実際に計算すると p r = m r ˙ p_r = m\dot r p r = m r ˙ なので
h = r ˙ ( m r ˙ ) − m 2 ( r ˙ 2 + r 2 ω 2 ) = m 2 r ˙ 2 − m 2 ω 2 r 2 h = \dot r\,(m\dot r) - \frac{m}{2}(\dot r^2 + r^2\omega^2) = \frac{m}{2}\dot r^2 - \frac{m}{2}\omega^2 r^2 h = r ˙ ( m r ˙ ) − 2 m ( r ˙ 2 + r 2 ω 2 ) = 2 m r ˙ 2 − 2 m ω 2 r 2 です。一方、運動エネルギーは T = m 2 ( r ˙ 2 + r 2 ω 2 ) T = \dfrac{m}{2}(\dot r^2 + r^2\omega^2) T = 2 m ( r ˙ 2 + r 2 ω 2 ) 、ポテンシャルは U = 0 U = 0 U = 0 ですから
h = T − m ω 2 r 2 ≠ T + U . h = T - m\omega^2 r^2 \ne T + U . h = T − m ω 2 r 2 = T + U . ずれの原因は、T T T が r ˙ \dot r r ˙ について 2 次同次でないこと(r 2 ω 2 r^2\omega^2 r 2 ω 2 という 0 次の項がある)です。Proposition 5.6 の条件 2 が破れています。
差が本当に効いていることを確かめます。オイラー–ラグランジュ方程式は
d d t ( m r ˙ ) = ∂ L ∂ r = m ω 2 r ⟹ r ¨ = ω 2 r \frac{d}{dt}(m\dot r) = \frac{\partial L}{\partial r} = m\omega^2 r
\quad\Longrightarrow\quad \ddot r = \omega^2 r d t d ( m r ˙ ) = ∂ r ∂ L = m ω 2 r ⟹ r ¨ = ω 2 r で、一般解は r ( t ) = A cosh ω t + B sinh ω t r(t) = A\cosh\omega t + B\sinh\omega t r ( t ) = A cosh ω t + B sinh ω t です。簡単のため A > 0 A > 0 A > 0 、B = 0 B = 0 B = 0 の解を取ると r ˙ = A ω sinh ω t \dot r = A\omega\sinh\omega t r ˙ = A ω sinh ω t なので
h = m 2 A 2 ω 2 sinh 2 ω t − m 2 ω 2 A 2 cosh 2 ω t = − m 2 A 2 ω 2 h = \frac{m}{2}A^2\omega^2\sinh^2\omega t - \frac{m}{2}\omega^2 A^2\cosh^2\omega t
= -\frac{m}{2}A^2\omega^2 h = 2 m A 2 ω 2 sinh 2 ω t − 2 m ω 2 A 2 cosh 2 ω t = − 2 m A 2 ω 2 (双曲線関数の関係 cosh 2 x − sinh 2 x = 1 \cosh^2 x - \sinh^2 x = 1 cosh 2 x − sinh 2 x = 1 を使いました)で、確かに定数です。ところが
T = m 2 A 2 ω 2 ( sinh 2 ω t + cosh 2 ω t ) = m 2 A 2 ω 2 cosh 2 ω t T = \frac{m}{2}A^2\omega^2\big(\sinh^2\omega t + \cosh^2\omega t\big) = \frac{m}{2}A^2\omega^2\cosh 2\omega t T = 2 m A 2 ω 2 ( sinh 2 ω t + cosh 2 ω t ) = 2 m A 2 ω 2 cosh 2 ω t は時間とともに指数関数的に増大します。ビーズは外へ飛ばされながら加速していきます。増えたエネルギーはモーターが供給しています。棒の回転を一定に保つために外部が仕事をしているので、系だけを見れば力学的エネルギーは保存しないのが当然です。それでも h h h という組み合わせは保存する、というのがこの例の教訓です。
ネーターの定理は古典力学の技法にとどまりません。この節では、この原理がどこまで届くかを短く見ておきます。
場の理論へ。 有限自由度の力学から連続体・場へ移ると、ネーターの定理は保存量ではなく保存流 を与えます。「保存する量が空間の各点に分布していて、局所的な連続の方程式 ∂ t ρ + ∇ ⋅ j = 0 \partial_t \rho + \nabla\cdot \boldsymbol{j} = 0 ∂ t ρ + ∇ ⋅ j = 0 に従う」という形になり、保存の内容が「総量が変わらない」から「どこかで減ったなら必ず隣へ流れた」へと強くなります。電荷保存則はこの形で現れます。導出は Appendix にまとめました。
素粒子物理の標準模型は、まさにこの原理の上に構築されています。どの粒子がどの力を感じるかは、ラグランジアンにどんな対称性を課すかでほぼ決まります。対称性を先に決め、ラグランジアンを後から作る。これは 20 世紀後半の理論物理の標準的な手順になりました。ネーターの定理は、その手順が意味を持つ理由を与えています。
Exercise 7.1 易
平面上を運動する質量 m m m の質点が、原点からの距離だけで決まるポテンシャル U ( r ) U(r) U ( r ) の中にあるとします。極座標でのラグランジアンは
L = m 2 ( r ˙ 2 + r 2 θ ˙ 2 ) − U ( r ) L = \frac{m}{2}\big(\dot r^2 + r^2\dot\theta^2\big) - U(r) L = 2 m ( r ˙ 2 + r 2 θ ˙ 2 ) − U ( r ) です。
循環座標を見つけ、対応する保存量を書き、それが角運動量に一致することを確かめてください。
h h h を計算し、この場合は h = T + U h = T + U h = T + U となることを Proposition 5.6 で確認してください。
Solution 1. L L L に θ \theta θ は現れず θ ˙ \dot\theta θ ˙ だけが現れるので ∂ L / ∂ θ = 0 \partial L/\partial\theta = 0 ∂ L / ∂ θ = 0 、すなわち θ \theta θ は循環座標です(Definition 2.2 )。Proposition 2.3 より
p θ = ∂ L ∂ θ ˙ = m r 2 θ ˙ p_\theta = \frac{\partial L}{\partial \dot\theta} = m r^2\dot\theta p θ = ∂ θ ˙ ∂ L = m r 2 θ ˙ が保存します。デカルト座標で x = r cos θ x = r\cos\theta x = r cos θ 、y = r sin θ y = r\sin\theta y = r sin θ と書くと
x y ˙ − y x ˙ = r cos θ ( r ˙ sin θ + r θ ˙ cos θ ) − r sin θ ( r ˙ cos θ − r θ ˙ sin θ ) = r 2 θ ˙ x\dot y - y\dot x = r\cos\theta\,(\dot r\sin\theta + r\dot\theta\cos\theta) - r\sin\theta\,(\dot r\cos\theta - r\dot\theta\sin\theta) = r^2\dot\theta x y ˙ − y x ˙ = r cos θ ( r ˙ sin θ + r θ ˙ cos θ ) − r sin θ ( r ˙ cos θ − r θ ˙ sin θ ) = r 2 θ ˙ なので、p θ = m ( x y ˙ − y x ˙ ) = J z p_\theta = m(x\dot y - y\dot x) = J_z p θ = m ( x y ˙ − y x ˙ ) = J z です。これは Example 4.4 で z z z 軸まわりの回転から得た角運動量に他なりません。
2. p r = m r ˙ p_r = m\dot r p r = m r ˙ 、p θ = m r 2 θ ˙ p_\theta = mr^2\dot\theta p θ = m r 2 θ ˙ なので
h = r ˙ ( m r ˙ ) + θ ˙ ( m r 2 θ ˙ ) − [ m 2 ( r ˙ 2 + r 2 θ ˙ 2 ) − U ( r ) ] = m 2 ( r ˙ 2 + r 2 θ ˙ 2 ) + U ( r ) h = \dot r (m\dot r) + \dot\theta (mr^2\dot\theta) - \left[ \frac{m}{2}(\dot r^2 + r^2\dot\theta^2) - U(r) \right]
= \frac{m}{2}\big(\dot r^2 + r^2\dot\theta^2\big) + U(r) h = r ˙ ( m r ˙ ) + θ ˙ ( m r 2 θ ˙ ) − [ 2 m ( r ˙ 2 + r 2 θ ˙ 2 ) − U ( r ) ] = 2 m ( r ˙ 2 + r 2 θ ˙ 2 ) + U ( r ) となり、h = T + U h = T + U h = T + U です。Proposition 5.6 で確認すると、U U U は速度を含まないので条件 1 が成立し、T T T は ( r ˙ , θ ˙ ) (\dot r, \dot\theta) ( r ˙ , θ ˙ ) を λ \lambda λ 倍すると λ 2 \lambda^2 λ 2 倍になるので条件 2 も成立します。極座標への変換 x = r cos θ x = r\cos\theta x = r cos θ が t t t を含まない(スクレロノミック)ことが効いています。Example 5.7 との違いはここです。
Exercise 7.2 標準
相互作用が相対位置だけで決まる N N N 質点系
L = ∑ a = 1 N m a 2 ∣ r ˙ a ∣ 2 − U ( { r a − r b } ) L = \sum_{a=1}^{N}\frac{m_a}{2}|\dot{\boldsymbol{r}}_a|^2 - U\big(\{\boldsymbol{r}_a - \boldsymbol{r}_b\}\big) L = a = 1 ∑ N 2 m a ∣ r ˙ a ∣ 2 − U ( { r a − r b } ) に、ガリレイ変換の族 R a ( s ) = r a + s t n \boldsymbol{R}_a(s) = \boldsymbol{r}_a + s\,t\,\boldsymbol{n} R a ( s ) = r a + s t n (n \boldsymbol{n} n は固定した単位ベクトル)を作用させます。
これが厳密な対称性ではなく準対称性であることを示し、F F F を求めてください。
保存量を求め、その物理的意味を述べてください。
Solution 1. ポテンシャルは相対位置 r a − r b \boldsymbol{r}_a - \boldsymbol{r}_b r a − r b だけの関数で、すべての質点が同じ s t n s\,t\,\boldsymbol{n} s t n だけずれるので相対位置は変わらず、U U U は s s s によりません。速度は R ˙ a = r ˙ a + s n \dot{\boldsymbol{R}}_a = \dot{\boldsymbol{r}}_a + s\,\boldsymbol{n} R ˙ a = r ˙ a + s n なので
∑ a m a 2 ∣ r ˙ a + s n ∣ 2 = ∑ a m a 2 ( ∣ r ˙ a ∣ 2 + 2 s r ˙ a ⋅ n + s 2 ) \sum_a \frac{m_a}{2}|\dot{\boldsymbol{r}}_a + s\boldsymbol{n}|^2
= \sum_a \frac{m_a}{2}\Big( |\dot{\boldsymbol{r}}_a|^2 + 2s\,\dot{\boldsymbol{r}}_a\cdot\boldsymbol{n} + s^2 \Big) a ∑ 2 m a ∣ r ˙ a + s n ∣ 2 = a ∑ 2 m a ( ∣ r ˙ a ∣ 2 + 2 s r ˙ a ⋅ n + s 2 ) です(∣ n ∣ = 1 |\boldsymbol{n}| = 1 ∣ n ∣ = 1 を使いました)。s s s で微分して s = 0 s = 0 s = 0 と置くと
∂ ∂ s ∣ s = 0 L = ∑ a m a r ˙ a ⋅ n ≠ 0 \frac{\partial}{\partial s}\bigg|_{s=0} L = \sum_a m_a \dot{\boldsymbol{r}}_a\cdot\boldsymbol{n} \ne 0 ∂ s ∂ s = 0 L = a ∑ m a r ˙ a ⋅ n = 0 なので厳密な対称性ではありません。しかしこれは全微分の形に書けます。全質量を M = ∑ a m a M = \sum_a m_a M = ∑ a m a 、重心を R c m = 1 M ∑ a m a r a \boldsymbol{R}_{\mathrm{cm}} = \dfrac{1}{M}\sum_a m_a\boldsymbol{r}_a R cm = M 1 ∑ a m a r a とすると
∑ a m a r ˙ a ⋅ n = d d t ( ∑ a m a r a ⋅ n ) = d d t ( M R c m ⋅ n ) \sum_a m_a\dot{\boldsymbol{r}}_a\cdot\boldsymbol{n} = \frac{d}{dt}\left( \sum_a m_a\boldsymbol{r}_a\cdot\boldsymbol{n} \right) = \frac{d}{dt}\big( M\boldsymbol{R}_{\mathrm{cm}}\cdot\boldsymbol{n} \big) a ∑ m a r ˙ a ⋅ n = d t d ( a ∑ m a r a ⋅ n ) = d t d ( M R cm ⋅ n ) なので、F = M R c m ⋅ n F = M\boldsymbol{R}_{\mathrm{cm}}\cdot\boldsymbol{n} F = M R cm ⋅ n と取れば Definition 3.2 の準対称性の条件が満たされます。
2. 生成子は ζ a = t n \boldsymbol{\zeta}_a = t\,\boldsymbol{n} ζ a = t n 、一般化運動量は m a r ˙ a m_a\dot{\boldsymbol{r}}_a m a r ˙ a なので Theorem 4.1 より
I = ∑ a m a r ˙ a ⋅ ( t n ) − M R c m ⋅ n = ( t P − M R c m ) ⋅ n I = \sum_a m_a\dot{\boldsymbol{r}}_a\cdot(t\boldsymbol{n}) - M\boldsymbol{R}_{\mathrm{cm}}\cdot\boldsymbol{n}
= \big( t\,\boldsymbol{P} - M\boldsymbol{R}_{\mathrm{cm}} \big)\cdot\boldsymbol{n} I = a ∑ m a r ˙ a ⋅ ( t n ) − M R cm ⋅ n = ( t P − M R cm ) ⋅ n が保存します。n \boldsymbol{n} n は任意なのでベクトル t P − M R c m t\boldsymbol{P} - M\boldsymbol{R}_{\mathrm{cm}} t P − M R cm の 3 成分すべてが保存します。Example 4.3 より P \boldsymbol{P} P も保存するので、R c m ( t ) = R c m ( 0 ) + ( P / M ) t \boldsymbol{R}_{\mathrm{cm}}(t) = \boldsymbol{R}_{\mathrm{cm}}(0) + (\boldsymbol{P}/M)\,t R cm ( t ) = R cm ( 0 ) + ( P / M ) t となり、これは「重心が等速直線運動する」という主張です。
なお F = 0 F = 0 F = 0 に限定していたらこの保存量は得られませんでした。Example 4.6 と同じく、準対称性まで許す Definition 3.2 の定式化が必要になる例です。
Exercise 7.3 標準
一様磁場 B = B e z \boldsymbol{B} = B\boldsymbol{e}_z B = B e z (B B B は定数)の中を運動する電荷 e e e 、質量 m m m の荷電粒子のラグランジアンは、ベクトルポテンシャルを A ( r ) = 1 2 B × r \boldsymbol{A}(\boldsymbol{r}) = \tfrac12\boldsymbol{B}\times\boldsymbol{r} A ( r ) = 2 1 B × r 、スカラーポテンシャルを 0 0 0 に取ると
L = m 2 ∣ r ˙ ∣ 2 + e r ˙ ⋅ A ( r ) L = \frac{m}{2}|\dot{\boldsymbol{r}}|^2 + e\,\dot{\boldsymbol{r}}\cdot\boldsymbol{A}(\boldsymbol{r}) L = 2 m ∣ r ˙ ∣ 2 + e r ˙ ⋅ A ( r ) と書けます。z z z 軸まわりの回転がこの L L L の対称性であることを確かめ、対応する保存量を x , y , x ˙ , y ˙ x, y, \dot x, \dot y x , y , x ˙ , y ˙ で書き下してください。
Solution z z z 軸まわりの角 s s s の回転を R ( s ) R(s) R ( s ) とします。B \boldsymbol{B} B は z z z 軸方向なので R ( s ) B = B R(s)\boldsymbol{B} = \boldsymbol{B} R ( s ) B = B です。回転行列 R R R は行列式が 1 1 1 の直交行列なので外積と可換、すなわち R ( u × v ) = ( R u ) × ( R v ) R(\boldsymbol{u}\times\boldsymbol{v}) = (R\boldsymbol{u})\times(R\boldsymbol{v}) R ( u × v ) = ( R u ) × ( R v ) が成り立ちます。よって
A ( R r ) = 1 2 B × ( R r ) = 1 2 ( R B ) × ( R r ) = R ( 1 2 B × r ) = R A ( r ) \boldsymbol{A}(R\boldsymbol{r}) = \tfrac12\,\boldsymbol{B}\times (R\boldsymbol{r}) = \tfrac12\,(R\boldsymbol{B})\times(R\boldsymbol{r}) = R\big( \tfrac12 \boldsymbol{B}\times\boldsymbol{r} \big) = R\,\boldsymbol{A}(\boldsymbol{r}) A ( R r ) = 2 1 B × ( R r ) = 2 1 ( R B ) × ( R r ) = R ( 2 1 B × r ) = R A ( r ) です。したがって変換後の相互作用項は、R R R が内積を保つことから
e ( R r ˙ ) ⋅ A ( R r ) = e ( R r ˙ ) ⋅ ( R A ( r ) ) = e r ˙ ⋅ A ( r ) e\,(R\dot{\boldsymbol{r}})\cdot\boldsymbol{A}(R\boldsymbol{r}) = e\,(R\dot{\boldsymbol{r}})\cdot(R\boldsymbol{A}(\boldsymbol{r})) = e\,\dot{\boldsymbol{r}}\cdot\boldsymbol{A}(\boldsymbol{r}) e ( R r ˙ ) ⋅ A ( R r ) = e ( R r ˙ ) ⋅ ( R A ( r )) = e r ˙ ⋅ A ( r ) で不変です。運動エネルギーも ∣ R r ˙ ∣ = ∣ r ˙ ∣ |R\dot{\boldsymbol{r}}| = |\dot{\boldsymbol{r}}| ∣ R r ˙ ∣ = ∣ r ˙ ∣ より不変なので、F = 0 F = 0 F = 0 の厳密な対称性です。
生成子は Example 4.4 と同じく ζ = e z × r \boldsymbol{\zeta} = \boldsymbol{e}_z\times\boldsymbol{r} ζ = e z × r です。一般化運動量は
p = ∂ L ∂ r ˙ = m r ˙ + e A \boldsymbol{p} = \frac{\partial L}{\partial \dot{\boldsymbol{r}}} = m\dot{\boldsymbol{r}} + e\boldsymbol{A} p = ∂ r ˙ ∂ L = m r ˙ + e A (これは m r ˙ m\dot{\boldsymbol{r}} m r ˙ ではないことに注意してください)なので、Theorem 4.1 より
I = p ⋅ ( e z × r ) = e z ⋅ ( r × p ) I = \boldsymbol{p}\cdot(\boldsymbol{e}_z\times\boldsymbol{r}) = \boldsymbol{e}_z\cdot(\boldsymbol{r}\times\boldsymbol{p}) I = p ⋅ ( e z × r ) = e z ⋅ ( r × p ) が保存します(スカラー三重積の巡回置換を使いました)。r × m r ˙ \boldsymbol{r}\times m\dot{\boldsymbol{r}} r × m r ˙ の z z z 成分は m ( x y ˙ − y x ˙ ) m(x\dot y - y\dot x) m ( x y ˙ − y x ˙ ) です。もう一方は、ベクトル三重積の公式 u × ( v × w ) = v ( u ⋅ w ) − w ( u ⋅ v ) \boldsymbol{u}\times(\boldsymbol{v}\times\boldsymbol{w}) = \boldsymbol{v}(\boldsymbol{u}\cdot\boldsymbol{w}) - \boldsymbol{w}(\boldsymbol{u}\cdot\boldsymbol{v}) u × ( v × w ) = v ( u ⋅ w ) − w ( u ⋅ v ) を使って
r × e A = e B 2 r × ( e z × r ) = e B 2 ( e z ∣ r ∣ 2 − r ( r ⋅ e z ) ) \boldsymbol{r}\times e\boldsymbol{A} = \frac{eB}{2}\,\boldsymbol{r}\times(\boldsymbol{e}_z\times\boldsymbol{r})
= \frac{eB}{2}\Big( \boldsymbol{e}_z |\boldsymbol{r}|^2 - \boldsymbol{r}\,(\boldsymbol{r}\cdot\boldsymbol{e}_z) \Big) r × e A = 2 e B r × ( e z × r ) = 2 e B ( e z ∣ r ∣ 2 − r ( r ⋅ e z ) ) なので、その z z z 成分は e B 2 ( ( x 2 + y 2 + z 2 ) − z 2 ) = e B 2 ( x 2 + y 2 ) \dfrac{eB}{2}\big( (x^2+y^2+z^2) - z^2 \big) = \dfrac{eB}{2}(x^2+y^2) 2 e B ( ( x 2 + y 2 + z 2 ) − z 2 ) = 2 e B ( x 2 + y 2 ) です。まとめると
I = m ( x y ˙ − y x ˙ ) + e B 2 ( x 2 + y 2 ) I = m(x\dot y - y\dot x) + \frac{eB}{2}\big( x^2 + y^2 \big) I = m ( x y ˙ − y x ˙ ) + 2 e B ( x 2 + y 2 ) が保存します。力学的角運動量 m ( x y ˙ − y x ˙ ) m(x\dot y - y\dot x) m ( x y ˙ − y x ˙ ) だけでは保存しない点が重要です。実際、磁場中の粒子は x y xy x y 平面内で円運動しますが、円の中心が原点になければ原点まわりの力学的角運動量は振動します。それを補正する項が第 2 項です。
Exercise 7.4 難
G ( q , t ) G(q,t) G ( q , t ) を C 2 C^2 C 2 級の任意の関数とし、L ′ = L + d G d t L' = L + \dfrac{dG}{dt} L ′ = L + d t d G と置きます(ここで d G d t = ∑ i ∂ G ∂ q i q ˙ i + ∂ G ∂ t \dfrac{dG}{dt} = \sum_i \dfrac{\partial G}{\partial q^i}\dot q^i + \dfrac{\partial G}{\partial t} d t d G = ∑ i ∂ q i ∂ G q ˙ i + ∂ t ∂ G です)。L L L と L ′ L' L ′ は同じオイラー–ラグランジュ方程式を与えることが知られています(全微分項による不定性(Proposition 4.4)[Lagrangian Mechanics] )。いま、時間を変えない一径数変換族 Q ( s ; q , t ) Q(s;q,t) Q ( s ; q , t ) が生成子 ζ \zeta ζ と関数 F F F とともに L L L の準対称性であるとします。
同じ変換族が L ′ L' L ′ の準対称性であることを示し、対応する F ′ F' F ′ を求めてください。
L L L から作ったネーター保存量と L ′ L' L ′ から作ったネーター保存量が一致することを示してください。
Solution 1. Definition 3.2 の左辺は L L L について線形なので、L ′ L' L ′ に対する ∂ s ∣ s = 0 \partial_s|_{s=0} ∂ s ∣ s = 0 は L L L の分と d G / d t dG/dt d G / d t の分の和になります。d G / d t dG/dt d G / d t の分を計算します。変換後の G G G の値は G ( Q ( s ; q ( t ) , t ) , t ) G(Q(s;q(t),t), t) G ( Q ( s ; q ( t ) , t ) , t ) であり、変換された曲線に沿ったその全微分の s s s 微分は、Theorem 4.1 の証明と同じくシュワルツの定理で s s s 微分と t t t 微分を交換して
∂ ∂ s ∣ s = 0 d d t G ( Q ( s ; q ( t ) , t ) , t ) = d d t ( ∂ ∂ s ∣ s = 0 G ( Q ( s ; q ( t ) , t ) , t ) ) = d d t ( ∑ i ∂ G ∂ q i ζ i ) \frac{\partial}{\partial s}\bigg|_{s=0}\frac{d}{dt}G\big(Q(s;q(t),t),t\big)
= \frac{d}{dt}\left( \frac{\partial}{\partial s}\bigg|_{s=0} G\big(Q(s;q(t),t),t\big) \right)
= \frac{d}{dt}\left( \sum_i \frac{\partial G}{\partial q^i}\zeta^i \right) ∂ s ∂ s = 0 d t d G ( Q ( s ; q ( t ) , t ) , t ) = d t d ( ∂ s ∂ s = 0 G ( Q ( s ; q ( t ) , t ) , t ) ) = d t d ( i ∑ ∂ q i ∂ G ζ i ) となります。したがって
∂ ∂ s ∣ s = 0 L ′ = d F d t + d d t ( ∑ i ∂ G ∂ q i ζ i ) = d d t ( F + ∑ i ∂ G ∂ q i ζ i ) \frac{\partial}{\partial s}\bigg|_{s=0} L' = \frac{dF}{dt} + \frac{d}{dt}\left( \sum_i \frac{\partial G}{\partial q^i}\zeta^i \right)
= \frac{d}{dt}\left( F + \sum_i \frac{\partial G}{\partial q^i}\zeta^i \right) ∂ s ∂ s = 0 L ′ = d t d F + d t d ( i ∑ ∂ q i ∂ G ζ i ) = d t d ( F + i ∑ ∂ q i ∂ G ζ i ) なので、L ′ L' L ′ もこの変換族について準対称性であり
F ′ = F + ∑ i ∂ G ∂ q i ζ i F' = F + \sum_i \frac{\partial G}{\partial q^i}\zeta^i F ′ = F + i ∑ ∂ q i ∂ G ζ i と取れます。
2. L ′ L' L ′ の一般化運動量は
p i ′ = ∂ L ′ ∂ q ˙ i = ∂ L ∂ q ˙ i + ∂ ∂ q ˙ i ( ∑ j ∂ G ∂ q j q ˙ j + ∂ G ∂ t ) = p i + ∂ G ∂ q i p'_i = \frac{\partial L'}{\partial \dot q^i} = \frac{\partial L}{\partial \dot q^i} + \frac{\partial}{\partial \dot q^i}\left( \sum_j \frac{\partial G}{\partial q^j}\dot q^j + \frac{\partial G}{\partial t} \right) = p_i + \frac{\partial G}{\partial q^i} p i ′ = ∂ q ˙ i ∂ L ′ = ∂ q ˙ i ∂ L + ∂ q ˙ i ∂ ( j ∑ ∂ q j ∂ G q ˙ j + ∂ t ∂ G ) = p i + ∂ q i ∂ G です(G G G は q ˙ \dot q q ˙ を含まないので ∂ G / ∂ q j \partial G/\partial q^j ∂ G / ∂ q j も ∂ G / ∂ t \partial G/\partial t ∂ G / ∂ t も q ˙ \dot q q ˙ に依存しません)。よって Theorem 4.1 の公式より
I ′ = ∑ i p i ′ ζ i − F ′ = ∑ i ( p i + ∂ G ∂ q i ) ζ i − F − ∑ i ∂ G ∂ q i ζ i = ∑ i p i ζ i − F = I I' = \sum_i p'_i \zeta^i - F'
= \sum_i \left( p_i + \frac{\partial G}{\partial q^i} \right)\zeta^i - F - \sum_i \frac{\partial G}{\partial q^i}\zeta^i
= \sum_i p_i \zeta^i - F = I I ′ = i ∑ p i ′ ζ i − F ′ = i ∑ ( p i + ∂ q i ∂ G ) ζ i − F − i ∑ ∂ q i ∂ G ζ i = i ∑ p i ζ i − F = I となり、両者は一致します。
つまりネーター保存量は、ラグランジアンの完全微分だけの不定性に依存しません。物理的に意味のある量であるためにはこの性質が必要であり、それが成り立っているわけです。なお時間変換 τ \tau τ を含む一般の場合(Theorem 5.3 )にも同じ結論が成り立ちます。その場合は h ′ = h − ∂ G / ∂ t h' = h - \partial G/\partial t h ′ = h − ∂ G / ∂ t となり、F ′ = F + ∑ i ( ∂ G / ∂ q i ) ζ i + ( ∂ G / ∂ t ) τ F' = F + \sum_i (\partial G/\partial q^i)\zeta^i + (\partial G/\partial t)\tau F ′ = F + ∑ i ( ∂ G / ∂ q i ) ζ i + ( ∂ G / ∂ t ) τ となって、やはり I ′ = I I' = I I ′ = I が示せます。
E. Noether, “Invariante Variationsprobleme”, Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen, Mathematisch-Physikalische Klasse (1918), 235–257. 英訳は M. A. Tavel によるものが arXiv:physics/0503066 で読めます。原論文。第一定理と第二定理が区別して述べられています。
L. D. ランダウ・E. M. リフシッツ『力学(増訂第 3 版)』東京図書 — 第 1 章・第 2 章。対称性から保存則を導く筋道を、最短の記述で通しています。エネルギー・運動量・角運動量の扱いはこの記事の §4・§5 に対応します。
H. Goldstein, C. Poole, J. Safko, Classical Mechanics , 3rd ed., Addison-Wesley, 2002 — 第 2 章(保存則と対称性)、第 13 章(連続体・場に対するネーターの定理)。準対称性やヤコビ積分の扱いが丁寧です。
V. I. Arnold, Mathematical Methods of Classical Mechanics , 2nd ed., Springer, 1989 — 第 4 章。ネーターの定理を多様体上の変分問題として、仮定を明示した形で定式化しています。
Y. Kosmann-Schwarzbach, The Noether Theorems: Invariance and Conservation Laws in the Twentieth Century , Springer, 2011. 原論文の英訳と、成立の歴史的背景・その後の受容の詳細な検討。
設定。 有限個の座標 q i ( t ) q^i(t) q i ( t ) の代わりに、時空の各点で値を取る場 ϕ ( x ) \phi(x) ϕ ( x ) (x = ( t , x ) x = (t, \boldsymbol{x}) x = ( t , x ) )を考えます。ラグランジアン密度 L ( ϕ , ∂ μ ϕ ) \mathcal{L}(\phi, \partial_\mu\phi) L ( ϕ , ∂ μ ϕ ) に対して作用は S = ∫ L d 4 x S = \displaystyle\int \mathcal{L}\, d^4x S = ∫ L d 4 x で、端点固定の変分から場のオイラー–ラグランジュ方程式
∂ μ ( ∂ L ∂ ( ∂ μ ϕ ) ) = ∂ L ∂ ϕ \partial_\mu\!\left( \frac{\partial\mathcal{L}}{\partial(\partial_\mu\phi)} \right) = \frac{\partial\mathcal{L}}{\partial\phi} ∂ μ ( ∂ ( ∂ μ ϕ ) ∂ L ) = ∂ ϕ ∂ L が得られます(μ \mu μ について和を取ります)。有限自由度の場合との対応は、添字 i i i が連続変数 x \boldsymbol{x} x に、d / d t d/dt d / d t が ∂ μ \partial_\mu ∂ μ に置き換わったものです。
内部対称性と保存流。 複素スカラー場 ϕ \phi ϕ について L = ∂ μ ϕ ˉ ∂ μ ϕ − m 2 ϕ ˉ ϕ \mathcal{L} = \partial_\mu\bar\phi\,\partial^\mu\phi - m^2\bar\phi\phi L = ∂ μ ϕ ˉ ∂ μ ϕ − m 2 ϕ ˉ ϕ を取ります。位相変換の族 ϕ ↦ e i s ϕ \phi \mapsto e^{is}\phi ϕ ↦ e i s ϕ 、ϕ ˉ ↦ e − i s ϕ ˉ \bar\phi\mapsto e^{-is}\bar\phi ϕ ˉ ↦ e − i s ϕ ˉ を考えると、L \mathcal{L} L の各項は ϕ ˉ \bar\phi ϕ ˉ と ϕ \phi ϕ を 1 つずつ含むので位相因子が打ち消し合い、L \mathcal{L} L は s s s によらず不変です。生成子は δ ϕ = i ϕ \delta\phi = i\phi δ ϕ = i ϕ 、δ ϕ ˉ = − i ϕ ˉ \delta\bar\phi = -i\bar\phi δ ϕ ˉ = − i ϕ ˉ です。
不変性を s s s で微分すると
0 = ∂ L ∂ ϕ δ ϕ + ∂ L ∂ ( ∂ μ ϕ ) ∂ μ ( δ ϕ ) + ∂ L ∂ ϕ ˉ δ ϕ ˉ + ∂ L ∂ ( ∂ μ ϕ ˉ ) ∂ μ ( δ ϕ ˉ ) 0 = \frac{\partial\mathcal{L}}{\partial\phi}\delta\phi
+ \frac{\partial\mathcal{L}}{\partial(\partial_\mu\phi)}\partial_\mu(\delta\phi)
+ \frac{\partial\mathcal{L}}{\partial\bar\phi}\delta\bar\phi
+ \frac{\partial\mathcal{L}}{\partial(\partial_\mu\bar\phi)}\partial_\mu(\delta\bar\phi) 0 = ∂ ϕ ∂ L δ ϕ + ∂ ( ∂ μ ϕ ) ∂ L ∂ μ ( δ ϕ ) + ∂ ϕ ˉ ∂ L δ ϕ ˉ + ∂ ( ∂ μ ϕ ˉ ) ∂ L ∂ μ ( δ ϕ ˉ ) です。ここで場の方程式を使って ∂ L / ∂ ϕ \partial\mathcal{L}/\partial\phi ∂ L / ∂ ϕ を ∂ μ ( ∂ L / ∂ ( ∂ μ ϕ ) ) \partial_\mu(\partial\mathcal{L}/\partial(\partial_\mu\phi)) ∂ μ ( ∂ L / ∂ ( ∂ μ ϕ )) で置き換えると、有限自由度の場合とまったく同じ積の微分公式により
0 = ∂ μ j μ , j μ = ∂ L ∂ ( ∂ μ ϕ ) δ ϕ + ∂ L ∂ ( ∂ μ ϕ ˉ ) δ ϕ ˉ 0 = \partial_\mu j^\mu, \qquad
j^\mu = \frac{\partial\mathcal{L}}{\partial(\partial_\mu\phi)}\,\delta\phi
+ \frac{\partial\mathcal{L}}{\partial(\partial_\mu\bar\phi)}\,\delta\bar\phi 0 = ∂ μ j μ , j μ = ∂ ( ∂ μ ϕ ) ∂ L δ ϕ + ∂ ( ∂ μ ϕ ˉ ) ∂ L δ ϕ ˉ を得ます。この L \mathcal{L} L では ∂ L / ∂ ( ∂ μ ϕ ) = ∂ μ ϕ ˉ \partial\mathcal{L}/\partial(\partial_\mu\phi) = \partial^\mu\bar\phi ∂ L / ∂ ( ∂ μ ϕ ) = ∂ μ ϕ ˉ なので
j μ = i ( ϕ ∂ μ ϕ ˉ − ϕ ˉ ∂ μ ϕ ) j^\mu = i\big( \phi\,\partial^\mu\bar\phi - \bar\phi\,\partial^\mu\phi \big) j μ = i ( ϕ ∂ μ ϕ ˉ − ϕ ˉ ∂ μ ϕ ) です。∂ μ j μ = 0 \partial_\mu j^\mu = 0 ∂ μ j μ = 0 を成分で書くと ∂ t j 0 + ∇ ⋅ j = 0 \partial_t j^0 + \nabla\cdot\boldsymbol{j} = 0 ∂ t j 0 + ∇ ⋅ j = 0 、すなわち連続の方程式です。
保存量への還元。 空間全体で積分し、場が無限遠で十分速く減衰するとすると、ガウスの発散定理より
d d t ∫ j 0 d 3 x = − ∫ ∇ ⋅ j d 3 x = 0 \frac{d}{dt}\int j^0\, d^3x = -\int \nabla\cdot\boldsymbol{j}\, d^3x = 0 d t d ∫ j 0 d 3 x = − ∫ ∇ ⋅ j d 3 x = 0 となり、Q = ∫ j 0 d 3 x Q = \displaystyle\int j^0\,d^3x Q = ∫ j 0 d 3 x が保存します。これが電荷です。有限自由度の場合の「保存量が一つ」に対応するのは、場の理論では「保存流が一本」であり、そこから積分して保存電荷が出ます。局所的な連続の方程式のほうが情報が多く、「総量は変わらないが、離れた 2 点で同時に消えたり現れたりする」ことを禁じています。
時空並進 x μ ↦ x μ + s a μ x^\mu \mapsto x^\mu + s\,a^\mu x μ ↦ x μ + s a μ に対して同じ手順を踏むと、保存流はエネルギー・運動量テンソル T μ ν T^{\mu\nu} T μν になり、∂ μ T μ ν = 0 \partial_\mu T^{\mu\nu} = 0 ∂ μ T μν = 0 が得られます。その ν = 0 \nu = 0 ν = 0 成分の空間積分が全エネルギー、ν = i \nu = i ν = i 成分の空間積分が全運動量です。Corollary 5.4 と Example 4.3 の場への一般化がこれにあたります。