tecktom.com

有限の変換:括弧を重ねた級数

2026年10月15日8分で読めますAI生成教科書『虎の巻』第1部 第11章 Poisson 括弧(生成子で動かす)
目次
  1. 流れと Lie 級数
  2. 回転と調和振動子
  3. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第1部 第11章 §1.11.4 ★★ 無限小の変換を積み重ねた有限の変換は、A+ϵ{A,G}+ϵ22!{{A,G},G}+⋯A+\epsilon\{A,G\}+\frac{\epsilon^2}{2!}\{\{A,G\},G\}+\cdots の級数(Lie 級数)で書ける

流れと Lie 級数

無限小変換 (1.11.1) を、ϵ\epsilon を小さく刻んで何回も繰り返すと、有限の大きさの変換になる。これは連立微分方程式

dqids=∂G∂pi,dpids=−∂G∂qi(1.11.14)\frac{\mathrm{d} q_i}{\mathrm{d} s}=\frac{\partial G}{\partial p_i},\qquad \frac{\mathrm{d} p_i}{\mathrm{d} s}=-\frac{\partial G}{\partial q_i}\tag{1.11.14}

をパラメータ ss について 00 から ϵ\epsilon まで解くことと同じである。G=HG=H, s=ts=t ならこれは Hamilton 方程式そのものになる。初期点 x=(q,p)x=(q,p) から出た解の s=ϵs=\epsilon での点を ΦϵG(x)\Phi^G_\epsilon(x) と書き、GG の流れと呼ぶ。図 1.11.1 の線は、この流れの軌道である。

関数 AA の値が流れに沿ってどう変わるかを、ϵ\epsilon の級数で表す。

命題 1.11.5(Lie 級数 ★★)

AA と GG が解析的(収束する冪級数で表せる)なら、∣ϵ∣\lvert \epsilon\rvert が十分小さい範囲で

A(ΦϵG(x))=∑n=0∞ϵnn!(DG nA)(x)=A+ϵ{A,G}+ϵ22!{{A,G},G}+ϵ33!{{{A,G},G},G}+⋯(1.11.15)A\bigl(\Phi^G_\epsilon(x)\bigr)=\sum_{n=0}^\infty\frac{\epsilon^n}{n!}\bigl(\mathsf D_G^{\,n}A\bigr)(x) =A+\epsilon\{A,G\}+\frac{\epsilon^2}{2!}\{\{A,G\},G\}+\frac{\epsilon^3}{3!}\{\{\{A,G\},G\},G\}+\cdots\tag{1.11.15}

が成り立つ。右辺は xx で評価し、DGA={A,G}\mathsf D_GA=\{A,G\} である。

この命題は、第1部第9章の式 (1.9.9)(配置空間の向き K\mathcal K の Lie 級数)を、相空間の流れの微分演算子 DG\mathsf D_G に当てたものである。

証明

a(s)≔A(ΦsG(x))a(s)\coloneqq A(\Phi^G_s(x)) とおく。式 (1.11.14) と連鎖律から、a′(s)a'(s) は {A,G}\{A,G\} を点 ΦsG(x)\Phi^G_s(x) で評価したものになる。{A,G}\{A,G\} も相空間の関数なので、同じ議論をくり返すと a(n)(s)=(DG nA)(ΦsG(x))a^{(n)}(s)=(\mathsf D_G^{\,n}A)(\Phi^G_s(x)) を得る。解析的な右辺を持つ常微分方程式の解は ss について解析的なので、a(s)a(s) は s=0s=0 のまわりで Taylor 級数に展開でき、s=0s=0 での係数が a(n)(0)=(DG nA)(x)a^{(n)}(0)=(\mathsf D_G^{\,n}A)(x) である。

形式的には、式 (1.11.15) は A(ΦϵG)=eϵDGAA(\Phi^G_\epsilon)=\mathrm{e}^{\epsilon\mathsf D_G}A と書くことができる。無限小変換 δA=ϵDGA\delta A=\epsilon\mathsf D_GA を指数関数にしたものが有限の変換になる。この形は 1.11.6節で、量子力学のユニタリ演算子 e−iϵG^/ℏ\mathrm{e}^{-\mathrm{i}\epsilon\hat G/\hbar} と並べる。

回転と調和振動子

回転の生成子 LzL_z で xx を動かす。DLzx={x,Lz}=−y\mathsf D_{L_z}x=\{x,L_z\}=-y, DLzy={y,Lz}=x\mathsf D_{L_z}y=\{y,L_z\}=x なので、DLz\mathsf D_{L_z} を当てるたびにx→−y→−x→y→xx\to-y\to-x\to y\to x と4回で元に戻る。式 (1.11.15) の偶数番目の項と奇数番目の項を分けて足すと

x(ΦθLz)=x(1−θ22!+θ44!−⋯ )−y(θ−θ33!+⋯ )=xcos⁡θ−ysin⁡θx\bigl(\Phi^{L_z}_\theta\bigr)=x\Bigl(1-\frac{\theta^2}{2!}+\frac{\theta^4}{4!}-\cdots\Bigr)-y\Bigl(\theta-\frac{\theta^3}{3!}+\cdots\Bigr)=x\cos\theta-y\sin\theta

となる。同じく y(ΦθLz)=xsin⁡θ+ycos⁡θy(\Phi^{L_z}_\theta)=x\sin\theta+y\cos\theta で、反時計回りの角度 θ\theta の回転が出る。無限小の回転の生成子を級数にすると、有限の角度の回転の行列が出る。

調和振動子 H=p22m+12mω2q2H=\frac{p^2}{2m}+\frac12m\omega^2q^2 の時間発展も同じ形になる。DHq={q,H}=p/m\mathsf D_Hq=\{q,H\}=p/m, DHp={p,H}=−mω2q\mathsf D_Hp=\{p,H\}=-m\omega^2q からDH2q=−ω2q\mathsf D_H^2q=-\omega^2q なので、偶数番目と奇数番目を分けて足すと

q(t)=qcos⁡ωt+pmωsin⁡ωt(1.11.16)q(t)=q\cos\omega t+\frac{p}{m\omega}\sin\omega t\tag{1.11.16}

を得る(右辺の q,pq,p は時刻 00 の値)。運動方程式を解かずに、括弧を重ねるだけで解が出た。

級数を有限の項で打ち切ると、短い時間では解をよく近似するが、時間が長くなると外れる。図 1.11.2 は q=1q=1, p=0p=0 のときの打ち切りを比べたもので、打ち切る次数を上げるほど cos⁡ωt\cos\omega t に重なる範囲が延びる(ωt=π/2\omega t=\pi/2 での誤差は、2次まで・4次まで・8次までで 0.230.23, 0.0200.020, 2.5×10−52.5\times10^{-5})。冪級数としては、この例は任意の tt で収束する。

図

図 1.11.2 調和振動子の時間発展の Lie 級数 (1.11.15)(q=1q=1, p=0p=0)を (ωt)n(\omega t)^n の項までで打ち切ったもの(実線)と、厳密解 cos⁡ωt\cos\omega t(点線)。打ち切る次数を上げるほど、厳密解に重なる時間の範囲が延びる

級数が終わる場合と収束しない場合

DG\mathsf D_G を何回か当てると 00 になる関数では、Lie 級数は有限項で終わり、有限の変換が多項式で正確に書ける(下の例題)。一方、AA や GG が解析的でなければ命題 1.11.5 は成り立つとは限らない。解析的でも、級数の収束半径が有限になることはある(流れが有限の時間で無限遠へ飛ぶ H=p2/2−q4H=p^2/2-q^4 などでは、解そのものが有限の時刻までしか存在しない)。その場合でも、十分小さい ϵ\epsilon では式 (1.11.15) がそのまま成り立つ。

例題 1.11.4 落下運動の Lie 級数は有限項で終わる ★★

オリジナルH=p22m+mgqH=\dfrac{p^2}{2m}+mgq について次に答えよ。

  1. (1) DHq\mathsf D_Hq, DH2q\mathsf D_H^2q, DH3q\mathsf D_H^3q を求めよ。
  2. (2) Lie 級数 (1.11.15) で q(t)q(t) と p(t)p(t) を求めよ(右辺の q,pq,p は時刻 00 の値)。
  3. (3) 生成子 G=pG=p による有限の変換で HH がどう変わるかを、Lie 級数で求めよ。
解答

(1) DHq={q,H}=p/m\mathsf D_Hq=\{q,H\}=p/m。DH2q={p/m,H}=1m{p,H}=1m(−mg)=−g\mathsf D_H^2q=\{p/m,H\}=\frac1m\{p,H\}=\frac1m(-mg)=-g。定数の括弧は 00 なので DH3q=0\mathsf D_H^3q=0。

(2) 級数は3項で終わり、q(t)=q+pmt−12gt2q(t)=q+\frac{p}{m}t-\frac12gt^2。同様に DHp=−mg\mathsf D_Hp=-mg, DH2p=0\mathsf D_H^2p=0 から p(t)=p−mgtp(t)=p-mgt。等加速度運動の公式が、括弧を2回重ねただけで正確に出る。

(3) DpH={H,p}=∂H∂q=mg\mathsf D_pH=\{H,p\}=\frac{\partial H}{\partial q}=mg、Dp2H={mg,p}=0\mathsf D_p^2H=\{mg,p\}=0 なので H(Φϵp)=H+ϵmgH(\Phi^p_\epsilon)=H+\epsilon mg。pp は qq を ϵ\epsilon だけ上へずらす並進なので、位置エネルギーが mgϵmg\epsilon 増える。{H,p}≠0\{H,p\}\neq0 なので、定理 1.11.3 から pp は保存しない(重力のポテンシャル mgqmgq は並進で値を変える)。

木下の勉強ステータス

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

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

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