解析力学から様々な理論分野を組み立てる 第2部 第4章 §2.4.1 ★★★ 原子の列の L L L を書き、ばね定数 K K K が原子間ポテンシャルの2階微分であることと、そこで調和近似を使ったことを確かめる
この部のここまでの章で扱ったのは、粒子が1個から数個の系だった。電磁場や重力場のような「場」は、空間の各点に1つずつ自由度を持つので、自由度が無限にある。自由度が無限にある系に Lagrangian をどう書けばよいかは、いきなり考えると分かりにくい。この章では、すでに書き方を知っている有限個の自由度の系から出発し、自由度を増やしていく極限として場の L L L を作る。
選ぶのは、質量 m m m の原子を N N N 個、ばね定数 K K K のばねで一直線につないだ列である。この列の L L L は
L = ∑ j 1 2 m u ˙ j 2 − ∑ j 1 2 K ( u j + 1 − u j ) 2 L=\sum_{j}\frac12m\dot u_j^2-\sum_j\frac12K(u_{j+1}-u_j)^2 L = j ∑ 2 1 m u ˙ j 2 − j ∑ 2 1 K ( u j + 1 − u j ) 2
の1行で書ける(u j u_j u j は j j j 番目の原子のつり合いの位置からのずれ)。これは結晶の格子振動のいちばん簡単な模型で、固体物理でフォノンと呼ぶものの古典版でもある。当てる道具は、第2部第1章 の小振動の方法(L L L を2次形式にして基準座標に分ける)と、第1部第10章 の Noether の定理である。列には「原子の番号を1つずらしても同じ」という対称性があり、この対称性によって基準座標の形が平面波に決まる。
出てくるものは3つある。1つめは、角振動数 ω \omega ω と波数 k k k の関係(分散関係)ω = 2 K / m ∣ sin ( k a / 2 ) ∣ \omega=2\sqrt{K/m}\,\lvert \sin(ka/2)\rvert ω = 2 K / m ∣ sin ( k a /2 )∣ である。2つめは、L L L が N N N 個の独立な調和振動子の和に分かれること、つまり基準座標である。3つめは、格子間隔 a a a を 0 0 0 に送ったとき、列の L L L が連続な弦の L = ∫ ( 1 2 ρ ∂ t u 2 − 1 2 Y ∂ x u 2 ) d x L=\int\bigl(\frac12\rho\,\partial_tu^2-\frac12Y\,\partial_xu^2\bigr)\mathrm{d} x L = ∫ ( 2 1 ρ ∂ t u 2 − 2 1 Y ∂ x u 2 ) d x になることである。最後の式の被積分関数がLagrangian 密度 で、これが 第2部第5章 以降で場を扱うときの出発点になる。連続の極限で何を捨てたか(短い波長のモードと、分散関係の曲がり)も式で確かめ、場の理論の「紫外の切断」への伏線として残す。
列の L L L
直線上に質量 m m m の原子が N N N 個並び、つり合いのとき j j j 番目の原子は x j = j a x_j=ja x j = j a にあるとする(a a a は格子定数、j = 0 , 1 , … , N − 1 j=0,1,\dots,N-1 j = 0 , 1 , … , N − 1 )。原子は列の向きにだけ動き、時刻 t t t の位置を j a + u j ( t ) ja+u_j(t) j a + u j ( t ) と書く。一般化座標は N N N 個のずれ u 0 , … , u N − 1 u_0,\dots,u_{N-1} u 0 , … , u N − 1 である。両端の扱いは、この節から第4節までは周期境界条件 u j + N = u j u_{j+N}=u_j u j + N = u j にする。列を輪につないだと考えてよい。輪にすると端が無くなり、どの原子も同じ立場になるので、対称性が見やすくなる。壁や自由端のある列は第5節と総合問題で扱う。
隣り合う原子は、ばね定数 K K K のばねでつながっているとする。j j j 番目と j + 1 j+1 j + 1 番目の間のばねの伸びは ( j a + a + u j + 1 ) − ( j a + u j ) − a = u j + 1 − u j (ja+a+u_{j+1})-(ja+u_j)-a=u_{j+1}-u_j ( j a + a + u j + 1 ) − ( j a + u j ) − a = u j + 1 − u j なので、
定義 2.4.1(原子の列の Lagrangian)
L = ∑ j = 0 N − 1 1 2 m u ˙ j 2 − ∑ j = 0 N − 1 1 2 K ( u j + 1 − u j ) 2 , u N ≡ u 0 (2.4.1) L=\sum_{j=0}^{N-1}\frac12m\dot u_j^2-\sum_{j=0}^{N-1}\frac12K(u_{j+1}-u_j)^2,\qquad u_N\equiv u_0\tag{2.4.1} L = j = 0 ∑ N − 1 2 1 m u ˙ j 2 − j = 0 ∑ N − 1 2 1 K ( u j + 1 − u j ) 2 , u N ≡ u 0 ( 2.4.1 )
を原子の列の Lagrangian とする。
図 2.4.1 ばねでつないだ原子の列。破線がつり合いの位置 j a ja j a 、黒丸が時刻 t t t の原子で、ずれ u j u_j u j が一般化座標になる。ばねの伸びは隣どうしのずれの差 u j + 1 − u j u_{j+1}-u_j u j + 1 − u j だけで決まる
運動エネルギーは原子ごとの和、位置エネルギーはばねごとの和である。位置エネルギーが「ずれの差」だけで書けていることが、この章の全部の結果の元になる。列全体を同じだけずらす u j → u j + ϵ u_j\to u_j+\epsilon u j → u j + ϵ では、どのばねの伸びも変わらない。
ばね定数はポテンシャルの2階微分
実際の結晶の原子の間にばねがあるわけではない。原子の間には距離 r r r だけによる相互作用のポテンシャル ϕ ( r ) \phi(r) ϕ ( r ) があり、つり合いの間隔 a a a は ϕ ′ ( a ) = 0 \phi'(a)=0 ϕ ′ ( a ) = 0 で決まる。隣どうしの距離は a + ( u j + 1 − u j ) a+(u_{j+1}-u_j) a + ( u j + 1 − u j ) なので、ϕ \phi ϕ を a a a のまわりで Taylor 展開すると
ϕ ( a + δ ) = ϕ ( a ) + 1 2 ϕ ′ ′ ( a ) δ 2 + 1 6 ϕ ′ ′ ′ ( a ) δ 3 + O ( δ 4 ) , δ = u j + 1 − u j (2.4.2) \phi(a+\delta)=\phi(a)+\frac12\phi''(a)\,\delta^2+\frac16\phi'''(a)\,\delta^3+\mathcal{O}(\delta^4),\qquad \delta=u_{j+1}-u_j\tag{2.4.2} ϕ ( a + δ ) = ϕ ( a ) + 2 1 ϕ ′′ ( a ) δ 2 + 6 1 ϕ ′′′ ( a ) δ 3 + O ( δ 4 ) , δ = u j + 1 − u j ( 2.4.2 )
になる。1次の項は ϕ ′ ( a ) = 0 \phi'(a)=0 ϕ ′ ( a ) = 0 で消える。定数 ϕ ( a ) \phi(a) ϕ ( a ) は運動方程式を変えないので落とし、3次以上の項を捨てると、K = ϕ ′ ′ ( a ) K=\phi''(a) K = ϕ ′′ ( a ) のばねと同じ式 (2.4.1) になる。3次以上を捨てる近似を調和近似 と呼ぶ。ここではポテンシャルを ϕ \phi ϕ と書き、運動エネルギーとの差を取る前の全体の位置エネルギー V V V と区別した。
調和近似が使えるかどうかは、ずれの差 δ \delta δ が a a a に比べて十分小さいかで決まる。3次の項と2次の項の比は ϕ ′ ′ ′ ( a ) δ / ( 3 ϕ ′ ′ ( a ) ) \phi'''(a)\delta/(3\phi''(a)) ϕ ′′′ ( a ) δ / ( 3 ϕ ′′ ( a )) で、例題 2.4.1 で見るように、典型的な原子間ポテンシャルではこの比が δ / a \delta/a δ / a の数倍になる。
原子を質点とみなすこと
この L L L は、原子を「位置だけを持つ質点」とし、原子どうしの相互作用を位置だけの関数 ϕ ( r ) \phi(r) ϕ ( r ) で書いている。実際の原子は原子核と電子からなり、ϕ ( r ) \phi(r) ϕ ( r ) は電子の状態を原子核の位置ごとに解いて得るエネルギーである。原子核が電子よりずっと重く遅いので、電子は原子核の位置に即座に追随するとみなす(Born–Oppenheimer 近似)。この近似の構造は 近似の数理 第7章 で、そもそも原子が周期的に並ぶことの理由(結晶化問題)は、数学としては未解決の問題で、その状況は 近似の数理 第13章 で述べる。
例題 2.4.1 Lennard-Jones 型のポテンシャルからばね定数を出す ★★
オリジナル 原子間のポテンシャルを ϕ ( r ) = V 0 [ ( r m / r ) 12 − 2 ( r m / r ) 6 ] \phi(r)=V_0\bigl[(r_{\mathrm m}/r)^{12}-2(r_{\mathrm m}/r)^6\bigr] ϕ ( r ) = V 0 [ ( r m / r ) 12 − 2 ( r m / r ) 6 ] (V 0 > 0 V_0>0 V 0 > 0 )とする。
(1) つり合いの間隔が a = r m a=r_{\mathrm m} a = r m であることを示し、ばね定数 K = ϕ ′ ′ ( r m ) K=\phi''(r_{\mathrm m}) K = ϕ ′′ ( r m ) を求めよ。
(2) ϕ ′ ′ ′ ( r m ) \phi'''(r_{\mathrm m}) ϕ ′′′ ( r m ) を求め、式 (2.4.2) の3次の項と2次の項の比を δ / r m \delta/r_{\mathrm m} δ / r m で表せ。
(3) 比が 0.1 0.1 0.1 になるずれ δ \delta δ は格子定数の何% か。
解答 (1) s = r m / r s=r_{\mathrm m}/r s = r m / r とおくと ϕ = V 0 ( s 12 − 2 s 6 ) \phi=V_0(s^{12}-2s^6) ϕ = V 0 ( s 12 − 2 s 6 ) 、d s d r = − s / r \frac{\mathrm{d} s}{\mathrm{d} r}=-s/r d r d s = − s / r なのでϕ ′ ( r ) = − 12 V 0 r ( s 12 − s 6 ) \phi'(r)=-\dfrac{12V_0}{r}(s^{12}-s^6) ϕ ′ ( r ) = − r 12 V 0 ( s 12 − s 6 ) 。r = r m r=r_{\mathrm m} r = r m (s = 1 s=1 s = 1 )で 0 0 0 になるので、つり合いの間隔は a = r m a=r_{\mathrm m} a = r m である。もう1回微分すると、r − n r^{-n} r − n の2階微分が n ( n + 1 ) r − n − 2 n(n+1)r^{-n-2} n ( n + 1 ) r − n − 2 であることから
ϕ ′ ′ ( r m ) = V 0 r m 2 ( 12 ⋅ 13 − 2 ⋅ 6 ⋅ 7 ) = 72 V 0 r m 2 = K . \phi''(r_{\mathrm m})=\frac{V_0}{r_{\mathrm m}^2}\bigl(12\cdot13-2\cdot6\cdot7\bigr)=\frac{72V_0}{r_{\mathrm m}^2}=K . ϕ ′′ ( r m ) = r m 2 V 0 ( 12 ⋅ 13 − 2 ⋅ 6 ⋅ 7 ) = r m 2 72 V 0 = K .
(2) r − n r^{-n} r − n の3階微分は − n ( n + 1 ) ( n + 2 ) r − n − 3 -n(n+1)(n+2)r^{-n-3} − n ( n + 1 ) ( n + 2 ) r − n − 3 なので
ϕ ′ ′ ′ ( r m ) = V 0 r m 3 ( − 12 ⋅ 13 ⋅ 14 + 2 ⋅ 6 ⋅ 7 ⋅ 8 ) = V 0 r m 3 ( − 2184 + 672 ) = − 1512 V 0 r m 3 . \phi'''(r_{\mathrm m})=\frac{V_0}{r_{\mathrm m}^3}\bigl(-12\cdot13\cdot14+2\cdot6\cdot7\cdot8\bigr)=\frac{V_0}{r_{\mathrm m}^3}(-2184+672)=-\frac{1512V_0}{r_{\mathrm m}^3}. ϕ ′′′ ( r m ) = r m 3 V 0 ( − 12 ⋅ 13 ⋅ 14 + 2 ⋅ 6 ⋅ 7 ⋅ 8 ) = r m 3 V 0 ( − 2184 + 672 ) = − r m 3 1512 V 0 .
3次の項と2次の項の比の大きさは
∣ ϕ ′ ′ ′ ∣ ∣ δ ∣ / 6 ϕ ′ ′ / 2 = 1512 3 ⋅ 72 ∣ δ ∣ r m = 7 ∣ δ ∣ r m . \frac{\lvert \phi'''\rvert \,\lvert \delta\rvert /6}{\phi''/2}=\frac{1512}{3\cdot72}\,\frac{\lvert \delta\rvert }{r_{\mathrm m}}=7\,\frac{\lvert \delta\rvert }{r_{\mathrm m}} . ϕ ′′ /2 ∣ ϕ ′′′ ∣ ∣ δ ∣ /6 = 3 ⋅ 72 1512 r m ∣ δ ∣ = 7 r m ∣ δ ∣ .
(3) 7 ∣ δ ∣ / r m = 0.1 7\lvert \delta\rvert /r_{\mathrm m}=0.1 7 ∣ δ ∣ / r m = 0.1 から ∣ δ ∣ / r m = 1 / 70 ≈ 0.014 \lvert \delta\rvert /r_{\mathrm m}=1/70\approx0.014 ∣ δ ∣ / r m = 1/70 ≈ 0.014 、つまり格子定数の約 1.4 % 1.4\% 1.4% である。比の係数が 1 1 1 でなく 7 7 7 になるのは、このポテンシャルが近づく側(r < r m r<r_{\mathrm m} r < r m )で急に立ち上がり、つり合いのまわりで強く非対称だからである。調和近似の「ずれが小さい」は、δ / a ≪ 1 \delta/a\ll1 δ / a ≪ 1 よりも一桁厳しい δ / a ≪ 1 / 7 \delta/a\ll1/7 δ / a ≪ 1/7 を意味する(検算 scripts/verify_ch2-04.py)。
近似の裏付け
型 B(極限) | 小さいもの ずれと格子定数の比 δ / a → 0 \delta/a\to0 δ / a → 0 (調和近似) | 証明 条件付き(振幅を 0 0 0 に送る極限で、有限の時間区間に限る。長い時間で非調和項がモードの間にエネルギーを移すかどうかは別の問題) | 行き先 近似の数理 第2章
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com