tecktom.com

出てくるもの2:基準座標で LL を分ける

2026年10月15日7分で読めますAI生成教科書『虎の巻』第2部 第4章 ばねでつないだ原子の列から場へ
目次
  1. 基準座標
  2. Hamiltonian と各モードのエネルギー
  3. のモードと Noether の保存量
  4. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第2部 第4章 §2.4.4 ★★★ 離散 Fourier 変換で LL は NN 個の独立な調和振動子の和になる。k=0k=0 の振動子が全運動量で、各モードのエネルギーは別々に保存する

基準座標

前の節では運動方程式に平面波を入れた。ここでは同じことを LL の段階で行い、LL そのものを独立な部分に分ける。基準座標 QkQ_k を、ずれ uju_j の離散 Fourier 変換

Qk=1N∑j=0N−1uj e−ikja,uj=1N∑kQk eikja(2.4.10)Q_k=\frac1{\sqrt N}\sum_{j=0}^{N-1}u_j\,\mathrm{e}^{-\mathrm{i} kja},\qquad u_j=\frac1{\sqrt N}\sum_kQ_k\,\mathrm{e}^{\mathrm{i} kja}\tag{2.4.10}

で定める(kk は Brillouin 域の中の NN 個の値)。前に 1/N1/\sqrt N を振り分けたので、この変換はユニタリで、逆変換で元に戻る。指数の符号は記法の約束(空間の順変換は e−ikx\mathrm{e}^{-\mathrm{i} kx})に合わせた。uju_j は実数なので Q−k=Qk∗Q_{-k}=Q_k^* で、NN 個の複素数 QkQ_k が持つ独立な実数の個数は NN のままである。

式 (2.4.10) を LL に代入する。直交関係 1N∑jei(k−k′)ja=δkk′\frac1N\sum_j\mathrm{e}^{\mathrm{i}(k-k')ja}=\delta_{kk'} を使うと、運動エネルギーは∑j12mu˙j2=∑k12mQ˙kQ˙−k\sum_j\frac12m\dot u_j^2=\sum_k\frac12m\dot Q_k\dot Q_{-k} になる。ばねの伸びは uj+1−uj=1N∑kQk(eika−1)eikjau_{j+1}-u_j=\frac1{\sqrt N}\sum_kQ_k(\mathrm{e}^{\mathrm{i} ka}-1)\mathrm{e}^{\mathrm{i} kja} で、∣eika−1∣2=4sin⁡2(ka/2)\lvert \mathrm{e}^{\mathrm{i} ka}-1\rvert ^2=4\sin^2(ka/2) から

L=∑k12m(Q˙kQ˙−k−ω(k)2QkQ−k)=∑k12m(∣Q˙k∣2−ω(k)2∣Qk∣2)(2.4.11)L=\sum_k\frac12m\bigl(\dot Q_k\dot Q_{-k}-\omega(k)^2Q_kQ_{-k}\bigr)=\sum_k\frac12m\bigl(\lvert \dot Q_k\rvert ^2-\omega(k)^2\lvert Q_k\rvert ^2\bigr)\tag{2.4.11}

になる。ω(k)\omega(k) は式 (2.4.7) である。原子の座標では隣どうしが結びついていた LL が、基準座標では kk ごとの和に分かれた。QkQ_k と Q−kQ_{-k} だけは互いに結びつくが、これは1つの実数の振動子を複素数で書いた結果で、QkQ_k の実部と虚部に分ければ実数の振動子2個になる。

Hamiltonian と各モードのエネルギー

QkQ_k に共役な運動量は、式 (2.4.11) で Q˙k\dot Q_k が kk の項と −k-k の項の2か所に現れることに注意して

Pk=∂L∂Q˙k=mQ˙−kP_k=\frac{\partial L}{\partial \dot Q_k}=m\dot Q_{-k}

である。第1部第8章の Legendre 変換 H=∑kPkQ˙k−LH=\sum_kP_k\dot Q_k-L から

H=∑k(PkP−k2m+12m ω(k)2QkQ−k)=∑kEk,Ek=12m(∣Q˙k∣2+ω(k)2∣Qk∣2)(2.4.12)H=\sum_k\Bigl(\frac{P_kP_{-k}}{2m}+\frac12m\,\omega(k)^2Q_kQ_{-k}\Bigr)=\sum_kE_k,\qquad E_k=\frac12m\bigl(\lvert \dot Q_k\rvert ^2+\omega(k)^2\lvert Q_k\rvert ^2\bigr)\tag{2.4.12}

を得る。Hamilton 方程式 Q˙k=∂H/∂Pk=P−k/m\dot Q_k=\partial H/\partial P_k=P_{-k}/m、P˙k=−∂H/∂Qk=−mω(k)2Q−k\dot P_k=-\partial H/\partial Q_k=-m\omega(k)^2Q_{-k} から Q¨k=−ω(k)2Qk\ddot Q_k=-\omega(k)^2Q_k で、EkE_k は kk ごとに保存する。列のエネルギーは NN 個の振動子に分かれて入り、調和近似の範囲では、ある振動子から別の振動子へは移らない。

k=0k=0 のモードと Noether の保存量

k=0k=0 では ω(0)=0\omega(0)=0 で、式 (2.4.11) の k=0k=0 の項は 12mQ˙02\frac12m\dot Q_0^2 だけになる。Q0Q_0 は LL に現れない循環座標で、共役な運動量 P0=mQ˙0P_0=m\dot Q_0 が保存する。式 (2.4.10) から Q0=1N∑jujQ_0=\frac1{\sqrt N}\sum_ju_j なので、

N P0=∑jmu˙j\sqrt N\,P_0=\sum_jm\dot u_j

は列の全運動量である。第2節の要点の箱で、一様な並進 uj→uj+ϵu_j\to u_j+\epsilon に Noether の定理を当てて得た保存量と同じものが、基準座標では「角振動数 00 のモードの運動量」として現れた。連続的な対称性があると、その向きに動かしても位置エネルギーが変わらないので、角振動数 00 のモードが必ず1つある。

例題 2.4.4 1個だけずらして放したときのエネルギーの分かれ方 ★★

オリジナルNN 個の原子の輪で、時刻 00 に原子 j=0j=0 だけを dd ずらし、ほかの原子は動かさず、全部を静かに放す。

  1. (1) 各モードのエネルギー EkE_k を求めよ。
  2. (2) ∑kEk\sum_kE_k が初めの位置エネルギーに等しいことを確かめよ。
  3. (3) いちばん多くのエネルギーを受け取るモードはどれか。その後の時刻で、そのモードのエネルギーはどうなるか。
解答

(1) 式 (2.4.10) で uj=d δj0u_j=d\,\delta_{j0} なので Qk(0)=d/NQ_k(0)=d/\sqrt N。静かに放したので Q˙k(0)=0\dot Q_k(0)=0。式 (2.4.12) から

Ek=12m ω(k)2d2N=2Kd2Nsin⁡2ka2.E_k=\frac12m\,\omega(k)^2\frac{d^2}{N}=\frac{2Kd^2}{N}\sin^2\frac{ka}{2}.

(2) k=2πl/(Na)k=2\pi l/(Na) について ∑l=0N−1sin⁡2(πl/N)=N/2\sum_{l=0}^{N-1}\sin^2(\pi l/N)=N/2(sin⁡2=12(1−cos⁡)\sin^2=\frac12(1-\cos) と、cos⁡(2πl/N)\cos(2\pi l/N) の和が 00 であることから)。よって ∑kEk=Kd2\sum_kE_k=Kd^2。初めの位置エネルギーは、原子 00 の左右の2本のばねがそれぞれ dd だけ伸び縮みしているので 2×12Kd2=Kd22\times\frac12Kd^2=Kd^2 で、一致する。

(3) sin⁡2(ka/2)\sin^2(ka/2) が最大の k=π/ak=\pi/a(NN が偶数のとき。隣どうしが逆向きに動くモード)が最大の 2Kd2/N2Kd^2/N を受け取る。1個だけをずらした変位は短い波長の成分を多く含むからである。各 EkE_k は保存するので、後の時刻でもこのモードのエネルギーは 2Kd2/N2Kd^2/N のままで、原子の運動が複雑に見えても、エネルギーのモードへの分かれ方は最初に決まったまま変わらない。N=6N=6 で運動方程式を時刻 20m/K20\sqrt{m/K} まで数値積分しても、各 EkE_k は初期値から変わらないことを確かめた(検算 scripts/verify_ch2-04.py)。

エネルギーがモードの間で混ざらないこと

統計力学では、十分長く待つとエネルギーがすべての自由度に等しく分かれる(等分配)と考える。ところが調和近似の列では、各 EkE_k が別々に保存するので、エネルギーは最初に入ったモードにとどまり、等分配は起こらない。混ぜる役は、捨てた3次以上の項(非調和項)にある。Fermi・Pasta・Ulam は、弱い非調和項を足した列を計算機で解き、予想に反してエネルギーがなかなか等分配されないことを見つけた(E. Fermi, J. Pasta, S. Ulam (1955), Los Alamos Scientific Laboratory Report LA-1940)。等分配の問題は 第3部第6章で扱う。

ほかの巻では

木下の勉強ステータス

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

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

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