tecktom.com

出てくるもの3:Hamiltonian 密度とエネルギーの流れ

2026年10月15日9分で読めますAI生成教科書『虎の巻』第2部 第5章 場の Lagrangian 密度と場の Euler–Lagrange 方程式
目次
  1. 正準運動量の密度と Hamiltonian 密度
  2. 場の Hamilton 方程式と Poisson 括弧
  3. エネルギーの保存と流れ
  4. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第2部 第5章 §2.5.5 ★★ 各点で Legendre 変換すると π(x)\pi(\bm{x}) と H\mathcal H ができる。{ϕ(x),π(y)}=δ(3)(x−y)\{\phi(\bm{x}),\pi(\bm{y})\}=\delta^{(3)}(\bm{x}-\bm{y}) が場の量子化の出発点

正準運動量の密度と Hamiltonian 密度

粒子では pi=∂L∂q˙ip_i=\frac{\partial L}{\partial \dot q_i} を定義し、速度について Legendre 変換して H=∑ipiq˙i−LH=\sum_ip_i\dot q_i-L を作った(第1部第8章)。場では一般化座標の番号 ii が空間の点 x\bm{x} に変わったので、各点で同じことをする。

定義 2.5.4(場の正準運動量と Hamiltonian 密度)

π(t,x)=∂L∂(∂tϕ),H=π ∂tϕ−L,H[ϕ,π]=∫d3x  H(2.5.22)\pi(t,\bm{x})=\frac{\partial \mathcal L}{\partial (\partial_t\phi)},\qquad \mathcal H=\pi\,\partial_t\phi-\mathcal L,\qquad H[\phi,\pi]=\int\mathrm{d}^3x\;\mathcal H\tag{2.5.22}

とし、H\mathcal H を ϕ\phi、∇ϕ\nabla\phi、π\pi の関数として書いたものをHamiltonian 密度と呼ぶ。

π\pi は、LL を「番号 x\bm{x} の速度 ∂tϕ(x)\partial_t\phi(\bm{x})」で汎関数微分したもの δLδ(∂tϕ(x))\frac{\delta L}{\delta(\partial_t\phi(\bm{x}))} に等しい。∑ipiq˙i\sum_ip_i\dot q_i の和が積分に変わって ∫π ∂tϕ d3x\int\pi\,\partial_t\phi\,\mathrm{d}^3x になるので、H=∫π ∂tϕ d3x−LH=\int\pi\,\partial_t\phi\,\mathrm{d}^3x-L は粒子の式そのものである。

弦では π=ρ ∂tϕ\pi=\rho\,\partial_t\phi で、

H=π22ρ+Y2(∂xϕ)2(2.5.23)\mathcal H=\frac{\pi^2}{2\rho}+\frac Y2(\partial_x\phi)^2\tag{2.5.23}

となる。運動エネルギーの密度と位置エネルギーの密度の和で、例題 2.5.1 のエネルギー密度 E\mathcal E と一致する。

場の Hamilton 方程式と Poisson 括弧

Hamilton 方程式は、偏微分を汎関数微分に置き換えて

∂tϕ(t,x)=δHδπ(x),∂tπ(t,x)=−δHδϕ(x)(2.5.24)\partial_t\phi(t,\bm{x})=\frac{\delta H}{\delta\pi(\bm{x})},\qquad \partial_t\pi(t,\bm{x})=-\frac{\delta H}{\delta\phi(\bm{x})}\tag{2.5.24}

になる。弦なら δHδπ=π/ρ\frac{\delta H}{\delta\pi}=\pi/\rho、δHδϕ=−Y∂x2ϕ\frac{\delta H}{\delta\phi}=-Y\partial_x^2\phi(例題 2.5.2 の計算)だから、∂tϕ=π/ρ\partial_t\phi=\pi/\rho と ∂tπ=Y∂x2ϕ\partial_t\pi=Y\partial_x^2\phi で、2つを合わせると波動方程式 (2.5.15) に戻る。2つの汎関数 F[ϕ,π]F[\phi,\pi]、G[ϕ,π]G[\phi,\pi] の Poisson 括弧は

{F,G}=∫d3x (δFδϕ(x)δGδπ(x)−δFδπ(x)δGδϕ(x))(2.5.25)\{F,G\}=\int\mathrm{d}^3x\,\Bigl(\frac{\delta F}{\delta\phi(\bm{x})}\frac{\delta G}{\delta\pi(\bm{x})}-\frac{\delta F}{\delta\pi(\bm{x})}\frac{\delta G}{\delta\phi(\bm{x})}\Bigr)\tag{2.5.25}

で定める(∑i\sum_i が ∫d3x\int\mathrm{d}^3x に変わっただけである)。F=ϕ(y)F=\phi(\bm{y})、G=π(z)G=\pi(\bm{z}) と置き、式 (2.5.9) を使うと

{ϕ(y),π(z)}=∫d3x  δ(3)(x−y) δ(3)(x−z)=δ(3)(y−z)(2.5.26)\{\phi(\bm{y}),\pi(\bm{z})\}=\int\mathrm{d}^3x\;\delta^{(3)}(\bm{x}-\bm{y})\,\delta^{(3)}(\bm{x}-\bm{z})=\delta^{(3)}(\bm{y}-\bm{z})\tag{2.5.26}

を得る。粒子の {qi,pj}=δij\{q_i,p_j\}=\delta_{ij} の場の版で、場の量子化(第5部第2章)はこの括弧を交換子 [⋅,⋅]/(iℏ)[\cdot,\cdot]/(\mathrm{i}\hbar) に置き換えるところから始まる。

Schrödinger 場 (2.5.18) では πψ=∂L∂(∂tψ)=iℏ ψ∗\pi_\psi=\frac{\partial \mathcal L}{\partial (\partial_t\psi)}=\mathrm{i}\hbar\,\psi^* で、正準運動量が新しい量でなく ψ∗\psi^* そのものになる。L\mathcal L が時間微分について1次なので、Legendre 変換で速度を消す必要がなく、ψ\psi と iℏψ∗\mathrm{i}\hbar\psi^* がはじめから正準変数の組になっている。H=πψ∂tψ−L=ℏ22m∣∇ψ∣2+V∣ψ∣2\mathcal H=\pi_\psi\partial_t\psi-\mathcal L=\frac{\hbar^2}{2m}\lvert \nabla\psi\rvert ^2+V\lvert \psi\rvert ^2 は、量子力学のエネルギーの期待値の被積分関数である。

エネルギーの保存と流れ

L\mathcal L が tt を陽に含まないとき、弦のエネルギー H=∫0ℓH dxH=\int_0^\ell\mathcal H\,\mathrm{d} x の時間変化を計算する。波動方程式を使うと

∂tH+∂xS=0,S=−Y ∂xϕ ∂tϕ(2.5.27)\partial_t\mathcal H+\partial_x\mathcal S=0,\qquad \mathcal S=-Y\,\partial_x\phi\,\partial_t\phi\tag{2.5.27}

が成り立つ(左辺を展開すると ∂tϕ (ρ ∂t2ϕ−Y∂x2ϕ)\partial_t\phi\,(\rho\,\partial_t^2\phi-Y\partial_x^2\phi) になる)。これを xx で積分すると

dHdt=−[S]0ℓ=Y∂xϕ ∂tϕ∣x=ℓ−Y∂xϕ ∂tϕ∣x=0.(2.5.28)\frac{\mathrm{d} H}{\mathrm{d} t}=-\bigl[\mathcal S\bigr]_{0}^{\ell}=Y\partial_x\phi\,\partial_t\phi\Big|_{x=\ell}-Y\partial_x\phi\,\partial_t\phi\Big|_{x=0} .\tag{2.5.28}

固定端では ∂tϕ=0\partial_t\phi=0、自由端では ∂xϕ=0\partial_x\phi=0 だから、どちらの端でも右辺は 00 で、HH は保存する。式 (2.5.27) は「エネルギー密度の増える速さ=流れ込むエネルギーの流れ」という局所的な保存則で、S\mathcal S はエネルギーの流れ(単位時間に点 xx を右へ通過するエネルギー)である。電磁場ではこれが Poynting ベクトルになる。L\mathcal L の時間並進の不変性から H\mathcal H と S\mathcal S を一般に作る手続きが場の Noether の定理で、次の次の章(第2部第7章)でエネルギー運動量テンソルとしてまとめる。

例題 2.5.5 定在波のエネルギーの出入り ★★

オリジナル両端を固定した長さ ℓ\ell の弦の基本モード ϕ=Asin⁡kxcos⁡ωt\phi=A\sin kx\cos\omega t(k=π/ℓk=\pi/\ell、ω=vk\omega=vk)について次に答えよ。

  1. (1) 運動エネルギー ∫0ℓπ22ρdx\int_0^\ell\frac{\pi^2}{2\rho}\mathrm{d} x と位置エネルギー ∫0ℓY2(∂xϕ)2dx\int_0^\ell\frac Y2(\partial_x\phi)^2\mathrm{d} x を求め、和 HH が時間によらないことを示せ。
  2. (2) ρ=0.10 kg m−1\rho=0.10\,\mathrm{kg\,m^{-1}}、Y=10 NY=10\,\mathrm N、ℓ=1.0 m\ell=1.0\,\mathrm m、A=1.0×10−3 mA=1.0\times10^{-3}\,\mathrm m のとき HH を求めよ。
  3. (3) 弦の中点 x=ℓ/2x=\ell/2 を通るエネルギーの流れ S\mathcal S を求め、左半分のエネルギーが時間によらないことを確かめよ。
解答

(1) π=ρ ∂tϕ=−ρAωsin⁡kxsin⁡ωt\pi=\rho\,\partial_t\phi=-\rho A\omega\sin kx\sin\omega t。∫0ℓsin⁡2kx dx=∫0ℓcos⁡2kx dx=ℓ/2\int_0^\ell\sin^2kx\,\mathrm{d} x=\int_0^\ell\cos^2kx\,\mathrm{d} x=\ell/2 を使うと

∫0ℓπ22ρdx=ρω2A2ℓ4sin⁡2ωt,∫0ℓY2(∂xϕ)2dx=Yk2A2ℓ4cos⁡2ωt=ρω2A2ℓ4cos⁡2ωt.\int_0^\ell\frac{\pi^2}{2\rho}\mathrm{d} x=\frac{\rho\omega^2A^2\ell}{4}\sin^2\omega t,\qquad \int_0^\ell\frac Y2(\partial_x\phi)^2\mathrm{d} x=\frac{Yk^2A^2\ell}{4}\cos^2\omega t=\frac{\rho\omega^2A^2\ell}{4}\cos^2\omega t .

和は H=14ρω2A2ℓH=\frac14\rho\omega^2A^2\ell で、時間によらない。エネルギーは運動と位置の間を行き来する。

(2) v=10 m s−1v=10\,\mathrm{m\,s^{-1}}、ω=πv/ℓ=31.4 s−1\omega=\pi v/\ell=31.4\,\mathrm s^{-1} で、H=14×0.10×(31.4)2×(1.0×10−3)2×1.0=2.47×10−5 JH=\frac14\times0.10\times(31.4)^2\times(1.0\times10^{-3})^2\times1.0=2.47\times10^{-5}\,\mathrm J。

(3) S=−Y∂xϕ ∂tϕ=YA2kωcos⁡kxsin⁡kxcos⁡ωtsin⁡ωt\mathcal S=-Y\partial_x\phi\,\partial_t\phi=YA^2k\omega\cos kx\sin kx\cos\omega t\sin\omega t。x=ℓ/2x=\ell/2 では cos⁡(π/2)=0\cos(\pi/2)=0 なので S=0\mathcal S=0 がすべての時刻で成り立つ。左半分 0≤x≤ℓ/20\le x\le\ell/2 には固定端からも中点からもエネルギーが流れ込まないので、式 (2.5.28) を左半分に当てはめれば、そのエネルギーは一定である。定在波は左右に進む2つの波の重ね合わせだが、全体としては腹と節の間でエネルギーが行き来するだけで、節を越えて流れない。

ほかの巻では

木下の勉強ステータス

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

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

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