解答 (1) [ Q k , P k ′ ] = 1 N ∑ j , l e − i k j a e i k ′ l a [ u j , p l ] = i ℏ N ∑ j e i ( k ′ − k ) j a = i ℏ δ k k ′ \left[Q_k,\,P_{k'}\right]=\frac1N\sum_{j,l}\mathrm{e}^{-\mathrm{i} kja}\mathrm{e}^{\mathrm{i} k'la}\left[u_j,\,p_l\right]=\frac{\mathrm{i}\hbar}{N}\sum_j\mathrm{e}^{\mathrm{i}(k'-k)ja}=\mathrm{i}\hbar\delta_{kk'} [ Q k , P k ′ ] = N 1 ∑ j , l e − i k j a e i k ′ l a [ u j , p l ] = N i ℏ ∑ j e i ( k ′ − k ) j a = i ℏ δ k k ′ 。u j u_j u j が実なので Q k † = Q − k Q_k^\dagger=Q_{-k} Q k † = Q − k 、P k † = P − k P_k^\dagger=P_{-k} P k † = P − k 。∑ j p j 2 = ∑ k P k P − k \sum_jp_j^2=\sum_kP_kP_{-k} ∑ j p j 2 = ∑ k P k P − k 、∑ j ( u j + 1 − u j ) 2 = ∑ k ∣ e i k a − 1 ∣ 2 Q k Q − k = ∑ k 4 sin 2 ( k a / 2 ) Q k Q − k \sum_j(u_{j+1}-u_j)^2=\sum_k\lvert \mathrm{e}^{\mathrm{i} ka}-1\rvert ^2Q_kQ_{-k}=\sum_k4\sin^2(ka/2)Q_kQ_{-k} ∑ j ( u j + 1 − u j ) 2 = ∑ k ∣ e i k a − 1 ∣ 2 Q k Q − k = ∑ k 4 sin 2 ( k a /2 ) Q k Q − k なので
H = ∑ k ( P k P − k 2 M + 1 2 M ω k 2 Q k Q − k ) . H=\sum_k\Bigl(\frac{P_kP_{-k}}{2M}+\frac12M\omega_k^2Q_kQ_{-k}\Bigr). H = k ∑ ( 2 M P k P − k + 2 1 M ω k 2 Q k Q − k ) .
これは式 (5.2.9) と同じ形で、ϕ k → M Q k \phi_{\bm{k}}\to\sqrt M\,Q_k ϕ k → M Q k 、π k → P − k / M \pi_{\bm{k}}\to P_{-k}/\sqrt M π k → P − k / M と対応する([ Q k , P − k † ] = [ Q k , P k ] = i ℏ \left[Q_k,\,P_{-k}^\dagger\right]=\left[Q_k,\,P_k\right]=\mathrm{i}\hbar [ Q k , P − k † ] = [ Q k , P k ] = i ℏ が式 (5.2.8) に当たる)。したがって例題 5.2.2 と同じ計算で [ b k , b k ′ † ] = δ k k ′ \left[b_k,\,b^\dagger_{k'}\right]=\delta_{kk'} [ b k , b k ′ † ] = δ k k ′ 、式 (5.2.13) と同じ計算で k ≠ 0 k\ne0 k = 0 の部分が ∑ ℏ ω k ( b k † b k + 1 2 ) \sum\hbar\omega_k(b^\dagger_kb_k+\frac12) ∑ ℏ ω k ( b k † b k + 2 1 ) になる。k = 0 k=0 k = 0 は ω 0 \omega_0 ω 0 の項が無い(ω k = 0 = 0 \omega_{k=0}=0 ω k = 0 = 0 )ので振動子にならず、P 0 2 / ( 2 M ) = ( ∑ j p j ) 2 / ( 2 N M ) P_0^2/(2M)=(\sum_jp_j)^2/(2NM) P 0 2 / ( 2 M ) = ( ∑ j p j ) 2 / ( 2 N M ) は列全体の重心の自由な運動である(第2部第4章の Noether の全運動量)。
(2) E 0 = ∑ k ≠ 0 ℏ ω k / 2 E_0=\sum_{k\ne0}\hbar\omega_k/2 E 0 = ∑ k = 0 ℏ ω k /2 。N → ∞ N\to\infty N → ∞ で 1 N ∑ k → a 2 π ∫ − π / a π / a d k \frac1N\sum_k\to\frac{a}{2\pi}\int_{-\pi/a}^{\pi/a}\mathrm{d} k N 1 ∑ k → 2 π a ∫ − π / a π / a d k 、θ = k a / 2 \theta=ka/2 θ = k a /2 とおいて
E 0 N = a 2 π ∫ − π / a π / a ℏ ω 0 ∣ sin k a 2 ∣ d k = ℏ ω 0 π ∫ − π / 2 π / 2 ∣ sin θ ∣ d θ = 2 ℏ ω 0 π = 0.637 ℏ ω 0 . \frac{E_0}{N}=\frac{a}{2\pi}\int_{-\pi/a}^{\pi/a}\hbar\omega_0\lvert \sin\frac{ka}{2}\rvert \mathrm{d} k=\frac{\hbar\omega_0}{\pi}\int_{-\pi/2}^{\pi/2}\lvert \sin\theta\rvert \,\mathrm{d}\theta=\frac{2\hbar\omega_0}{\pi}=0.637\,\hbar\omega_0 . N E 0 = 2 π a ∫ − π / a π / a ℏ ω 0 ∣ sin 2 k a ∣ d k = π ℏ ω 0 ∫ − π /2 π /2 ∣ sin θ ∣ d θ = π 2ℏ ω 0 = 0.637 ℏ ω 0 .
2 ℏ ω 0 2\hbar\omega_0 2ℏ ω 0 との比は 1 / π = 0.318 1/\pi=0.318 1/ π = 0.318 。格子では波数に上限 π / a \pi/a π / a があるので、5.2.3 節の場と違って原子1個あたりの零点エネルギーは有限である。場の零点エネルギーの発散は、格子間隔 a → 0 a\to0 a → 0 で波数の上限が無くなることに当たる。
(3) u j = N − 1 / 2 ∑ k Q k e i k j a u_j=N^{-1/2}\sum_kQ_k\mathrm{e}^{\mathrm{i} kja} u j = N − 1/2 ∑ k Q k e i k j a で、Q k = ℏ / ( 2 M ω k ) ( b k + b − k † ) Q_k=\sqrt{\hbar/(2M\omega_k)}\,(b_k+b^\dagger_{-k}) Q k = ℏ/ ( 2 M ω k ) ( b k + b − k † ) (k ≠ 0 k\ne0 k = 0 )。真空では ⟨ Q k Q k ′ ⟩ = ℏ 2 M ω k δ k ′ , − k \langle Q_kQ_{k'} \rangle =\frac{\hbar}{2M\omega_k}\delta_{k',-k} ⟨ Q k Q k ′ ⟩ = 2 M ω k ℏ δ k ′ , − k (第4部第4章の ⟨ q 2 ⟩ = ℏ / ( 2 m ω ) \langle q^2 \rangle =\hbar/(2m\omega) ⟨ q 2 ⟩ = ℏ/ ( 2 mω ) のモード版)なので
⟨ ( u j − u j + r ) 2 ⟩ = 1 N ∑ k ≠ 0 ∣ 1 − e i k r a ∣ 2 ℏ 2 M ω k = ℏ M N ∑ k ≠ 0 1 − cos k r a ω k . \langle (u_j-u_{j+r})^2 \rangle =\frac1N\sum_{k\ne0}\lvert 1-\mathrm{e}^{\mathrm{i} kra}\rvert ^2\frac{\hbar}{2M\omega_k}=\frac{\hbar}{MN}\sum_{k\ne0}\frac{1-\cos kra}{\omega_k}. ⟨( u j − u j + r ) 2 ⟩ = N 1 k = 0 ∑ ∣ 1 − e i k r a ∣ 2 2 M ω k ℏ = M N ℏ k = 0 ∑ ω k 1 − cos k r a .
(k = 0 k=0 k = 0 は 1 − cos 0 = 0 1-\cos0=0 1 − cos 0 = 0 なので相対的なずれに寄与しない。)N → ∞ N\to\infty N → ∞ で積分にし、θ = k a / 2 \theta=ka/2 θ = k a /2 、1 − cos 2 r θ = 2 sin 2 r θ 1-\cos2r\theta=2\sin^2r\theta 1 − cos 2 r θ = 2 sin 2 r θ を使うと
⟨ ( u j − u j + r ) 2 ⟩ = ℏ π M ω 0 ∫ − π / 2 π / 2 sin 2 r θ ∣ sin θ ∣ d θ = 2 ℏ π M ω 0 ∑ n = 1 r ∫ 0 π / 2 sin ( ( 2 n − 1 ) θ ) d θ = 2 ℏ π M ω 0 ∑ n = 1 r 1 2 n − 1 . \langle (u_j-u_{j+r})^2 \rangle =\frac{\hbar}{\pi M\omega_0}\int_{-\pi/2}^{\pi/2}\frac{\sin^2r\theta}{\lvert \sin\theta\rvert }\mathrm{d}\theta
=\frac{2\hbar}{\pi M\omega_0}\sum_{n=1}^r\int_0^{\pi/2}\sin((2n-1)\theta)\,\mathrm{d}\theta=\frac{2\hbar}{\pi M\omega_0}\sum_{n=1}^r\frac{1}{2n-1}. ⟨( u j − u j + r ) 2 ⟩ = π M ω 0 ℏ ∫ − π /2 π /2 ∣ sin θ ∣ sin 2 r θ d θ = π M ω 0 2ℏ n = 1 ∑ r ∫ 0 π /2 sin (( 2 n − 1 ) θ ) d θ = π M ω 0 2ℏ n = 1 ∑ r 2 n − 1 1 .
S r S_r S r は 1 2 ln r \frac12\ln r 2 1 ln r のように際限なく増えるので、基底状態でも遠く離れた原子の相対位置の揺らぎは r → ∞ r\to\infty r → ∞ で発散する。被積分関数 ( 1 − cos k r a ) / ω k (1-\cos kra)/\omega_k ( 1 − cos k r a ) / ω k は小さな k k k で ω k ≈ a ω 0 ∣ k ∣ \omega_k\approx a\omega_0\lvert k\rvert ω k ≈ a ω 0 ∣ k ∣ のため 1 / ∣ k ∣ 1/\lvert k\rvert 1/ ∣ k ∣ のように振る舞い、∫ d k / ∣ k ∣ \int\mathrm{d} k/\lvert k\rvert ∫ d k / ∣ k ∣ の対数の発散が r r r の対数として現れる。
(4) ⟨ ( u j − u j + r ) 2 ⟩ = ( 0.1 a ) 2 \langle (u_j-u_{j+r})^2 \rangle =(0.1a)^2 ⟨( u j − u j + r ) 2 ⟩ = ( 0.1 a ) 2 とおくと 2 π × 1.0 × 10 − 3 × S r = 0.01 \frac{2}{\pi}\times1.0\times10^{-3}\times S_r=0.01 π 2 × 1.0 × 1 0 − 3 × S r = 0.01 、すなわち S r = 5 π = 15.71 S_r=5\pi=15.71 S r = 5 π = 15.71 。1 2 ln r + ln 2 + γ 2 = 15.71 \frac12\ln r+\ln2+\frac\gamma2=15.71 2 1 ln r + ln 2 + 2 γ = 15.71 を解いて ln r = 29.45 \ln r=29.45 ln r = 29.45 、r = 6.2 × 10 12 r=6.2\times10^{12} r = 6.2 × 1 0 12 。隣の原子との相対的なずれは ( 2 / π ) × 10 − 3 a = 0.025 a \sqrt{(2/\pi)\times10^{-3}}\,a=0.025\,a ( 2/ π ) × 1 0 − 3 a = 0.025 a と小さく、それが 0.1 a 0.1a 0.1 a に育つには 10 12 10^{12} 1 0 12 個を超える原子をまたぐ必要がある。対数の増え方は非常に遅い。
連続の極限では、原子の列は第2部第4章の弦、すなわち 1+1 次元の質量 0 0 0 の場になる。例題 5.2.5 の場で同じ量を計算すると∫ d k 2 π 1 − cos k x ∣ k ∣ \int\frac{\mathrm{d} k}{2\pi}\frac{1-\cos kx}{\lvert k\rvert } ∫ 2 π d k ∣ k ∣ 1 − c o s k x が小さな k k k (長い波長)で対数発散する。質量 0 0 0 の場では ω k → 0 \omega_k\to0 ω k → 0 のモードの零点の揺らぎ ℏ / ( 2 ω k ) \hbar/(2\omega_k) ℏ/ ( 2 ω k ) が大きくなり、1次元ではそれを積分した量が赤外(長い波長)で発散する。質量があれば ω k ≥ m \omega_k\ge m ω k ≥ m なので発散しない。3次元の結晶では、同じ和が ∫ d 3 k ( 1 − cos k ⋅ r ) / ∣ k ∣ \int\mathrm{d}^3k\,(1-\cos\bm{k}\cdot\bm{r})/\lvert \bm{k}\rvert ∫ d 3 k ( 1 − cos k ⋅ r ) / ∣ k ∣ になり、d 3 k = 4 π k 2 d k \mathrm{d}^3k=4\pi k^2\mathrm{d} k d 3 k = 4 π k 2 d k の因子 k 2 k^2 k 2 のおかげで小さな k k k で収束するので、r → ∞ r\to\infty r → ∞ でも相対的なずれは有限にとどまり、結晶の長距離の秩序は量子揺らぎで壊れない。1次元だけで秩序が対数的に壊れるのは、小さな k k k のモードの数が次元で違うためである。