tecktom.com

橋:第二量子化・第1章の流体・ほかの部との対応

2026年10月17日11分で読めますAI生成教科書『虎の巻』第7部 第3章 量子流体:渦の量子化と Gross–Pitaevskii 方程式の流体形
目次
  1. 場の演算子から GP 方程式へ
  2. 第7部第1章の流体との関係
  3. 地図の上の位置と授業の形
  4. 木下の層

AI 生成・木下が未勉強

量子力学の様々な構成方法 第7部 第3章 §7.3.6 ★★★ GP 方程式は第二量子化の場の演算子の Heisenberg 方程式の期待値を、凝縮した状態で因数分解したものである。流体はこれで多粒子でも 3 次元の空間に戻る

場の演算子から GP 方程式へ

第6部の第二量子化(第6部第1章・第6部第2章)では、ボース粒子の場の演算子 ψ^(r)\hat\psi(\bm{r}) が[ψ^(r), ψ^†(r′)]=δ(3)(r−r′)\left[\hat\psi(\bm{r}),\,\hat\psi^\dagger(\bm{r}')\right]=\delta^{(3)}(\bm{r}-\bm{r}') を満たし、(7.3.1) は

H^=∫ψ^†(−ℏ22m∇2+V)ψ^ d3r+g2∫ψ^†ψ^†ψ^ψ^ d3r(7.3.26)\hat H=\int\hat\psi^\dagger\Bigl(-\frac{\hbar^2}{2m}\nabla^2+V\Bigr)\hat\psi\,\mathrm{d}^3r+\frac g2\int\hat\psi^\dagger\hat\psi^\dagger\hat\psi\hat\psi\,\mathrm{d}^3r\tag{7.3.26}

と書ける。Heisenberg の運動方程式(第2部第3章)iℏ∂tψ^=[ψ^, H^]\mathrm{i}\hbar\partial_t\hat\psi=\left[\hat\psi,\,\hat H\right] を計算する。交換関係から [ψ^(r), ψ^†(r′)ψ^†(r′)ψ^(r′)ψ^(r′)]=2δ(3)(r−r′)ψ^†(r′)ψ^(r′)ψ^(r′)\left[\hat\psi(\bm{r}),\,\hat\psi^\dagger(\bm{r}')\hat\psi^\dagger(\bm{r}')\hat\psi(\bm{r}')\hat\psi(\bm{r}')\right]=2\delta^{(3)}(\bm{r}-\bm{r}')\hat\psi^\dagger(\bm{r}')\hat\psi(\bm{r}')\hat\psi(\bm{r}') となるので([a^, a^†a^†]=2a^†\left[\hat a,\,\hat a^\dagger\hat a^\dagger\right]=2\hat a^\dagger と同じ計算)、

iℏ ∂tψ^=−ℏ22m∇2ψ^+Vψ^+g ψ^†ψ^ψ^(7.3.27)\mathrm{i}\hbar\,\partial_t\hat\psi=-\frac{\hbar^2}{2m}\nabla^2\hat\psi+V\hat\psi+g\,\hat\psi^\dagger\hat\psi\hat\psi\tag{7.3.27}

を得る。これは演算子の方程式で、近似は入っていない。両辺の期待値をとり ψ=⟨ψ^⟩\psi=\langle \hat\psi \rangle と書くと、左辺と最初の2項は ψ\psi だけで書けるが、最後の項は ⟨ψ^†ψ^ψ^⟩\langle \hat\psi^\dagger\hat\psi\hat\psi \rangle で、一般には ∣ψ∣2ψ\lvert \psi\rvert ^2\psi にならない。

凝縮体を、凝縮のモードについてのコヒーレント状態(ψ^\hat\psi の固有状態、第6部第6章(第6部第6章)が経路積分の基底に使ったもの)とみなすと、ψ^∣Φ⟩=ψ∣Φ⟩\hat\psi| \Phi \rangle =\psi| \Phi \rangle から ⟨ψ^†ψ^ψ^⟩=∣ψ∣2ψ\langle \hat\psi^\dagger\hat\psi\hat\psi \rangle =\lvert \psi\rvert ^2\psi が厳密に成り立ち、(7.3.27) の期待値は GP 方程式 (7.3.6) そのものになる。コヒーレント状態は時間発展で厳密には保たれないので、ここが平均場の近似である。ψ^=ψ+δψ^\hat\psi=\psi+\delta\hat\psi と分けて δψ^\delta\hat\psi の1次まで残せば、第3節の線形化を演算子で行ったことになり、Bogoliubov の励起 (7.3.15) が出る。流体の振動の量子が Bogoliubov 準粒子である。

この見方では、巨視的波動関数 ψ\psi は場の演算子の期待値で、その位相 SS は凝縮体全体が選んだ1つの値である。粒子数の決まった状態 ∣N⟩| N \rangle では ⟨N∣ψ^∣N⟩=0\langle N | \hat\psi | N \rangle =0 で、位相は定まらない。ψ≠0\psi\neq0 と置くことは、ハミルトニアン (7.3.26) が持つ位相の回転 ψ^→eiαψ^\hat\psi\to\mathrm{e}^{\mathrm{i}\alpha}\hat\psi の対称性(粒子数の保存と同じもの)を状態の側で破ることに当たる。有限の系では、この破れは粒子数の決まった状態の重ね合わせで書けるが、無限の系では位相の異なる状態どうしが別の表現に属するようになる。この区別は第8部第3章(第8部第3章)の非同値な表現と、第8部第4章(第8部第4章)の自発的対称性の破れの言葉で正確になる。粒子数を固定したまま凝縮を定義するには、1体の密度行列 ⟨ψ^†(r)ψ^(r′)⟩\langle \hat\psi^\dagger(\bm{r})\hat\psi(\bm{r}') \rangle が ∣r−r′∣→∞\lvert \bm{r}-\bm{r}'\rvert \to\infty で 00 にならず ψ∗(r)ψ(r′)\psi^*(\bm{r})\psi(\bm{r}') に近づくという条件(Penrose と Onsager の非対角長距離秩序)を使う。こちらは第4部の部分トレースで作った縮約密度行列の言葉である。

例題 7.3.6 1モードで見る平均場の因数分解 ★★

オリジナル凝縮のモードを1つの消滅演算子 a^\hat a で表し、相互作用を H^int=ℏU2a^†a^†a^a^\hat H_{\mathrm{int}}=\frac{\hbar U}{2}\hat a^\dagger\hat a^\dagger\hat a\hat a とする(UU は角振動数の次元の定数)。

  1. [a^, a^†a^†a^a^]=2a^†a^a^\left[\hat a,\,\hat a^\dagger\hat a^\dagger\hat a\hat a\right]=2\hat a^\dagger\hat a\hat a を示し、a^\hat a の Heisenberg 方程式を書け。
  2. コヒーレント状態 ∣α⟩| \alpha \rangle(a^∣α⟩=α∣α⟩\hat a| \alpha \rangle =\alpha| \alpha \rangle)で ⟨a^†a^a^⟩=∣α∣2α\langle \hat a^\dagger\hat a\hat a \rangle =\lvert \alpha\rvert ^2\alpha となることを示し、α(t)\alpha(t) の平均場の方程式を解け。
  3. 粒子数の決まった状態 ∣N⟩| N \rangle で ⟨a^⟩\langle \hat a \rangle を求め、(2) の α\alpha の位相が何を表すかを述べよ。
解答

(1) [a^, a^†a^†]=a^†[a^, a^†]+[a^, a^†]a^†=2a^†\left[\hat a,\,\hat a^\dagger\hat a^\dagger\right]=\hat a^\dagger\left[\hat a,\,\hat a^\dagger\right]+\left[\hat a,\,\hat a^\dagger\right]\hat a^\dagger=2\hat a^\dagger で、a^\hat a は a^a^\hat a\hat a と交換するので[a^, a^†a^†a^a^]=[a^, a^†a^†]a^a^=2a^†a^a^\left[\hat a,\,\hat a^\dagger\hat a^\dagger\hat a\hat a\right]=\left[\hat a,\,\hat a^\dagger\hat a^\dagger\right]\hat a\hat a=2\hat a^\dagger\hat a\hat a。Heisenberg 方程式は iℏa^˙=[a^, H^int]=ℏUa^†a^a^\mathrm{i}\hbar\dot{\hat a}=\left[\hat a,\,\hat H_{\mathrm{int}}\right]=\hbar U\hat a^\dagger\hat a\hat a となる。(7.3.27) の最後の項の1モード版である。

(2) ⟨α∣a^†=α∗⟨α∣\langle \alpha |\hat a^\dagger=\alpha^*\langle \alpha | と a^a^∣α⟩=α2∣α⟩\hat a\hat a| \alpha \rangle =\alpha^2| \alpha \rangle から ⟨α∣a^†a^a^∣α⟩=α∗α2=∣α∣2α\langle \alpha | \hat a^\dagger\hat a\hat a | \alpha \rangle =\alpha^*\alpha^2=\lvert \alpha\rvert ^2\alpha。平均場の方程式は iα˙=U∣α∣2α\mathrm{i}\dot\alpha=U\lvert \alpha\rvert ^2\alpha で、∣α∣2\lvert \alpha\rvert ^2 は保存するので α(t)=α(0) e−iU∣α∣2t\alpha(t)=\alpha(0)\,\mathrm{e}^{-\mathrm{i} U\lvert \alpha\rvert ^2t}。これは一様な凝縮体 ψ=n0 e−ign0t/ℏ\psi=\sqrt{n_0}\,\mathrm{e}^{-\mathrm{i} gn_0t/\hbar}(位相が μ/ℏ\mu/\hbar で回る)の1モード版で、U∣α∣2U\lvert \alpha\rvert ^2 が μ/ℏ\mu/\hbar に当たる。厳密な時間発展では、各 Fock 成分 ∣n⟩| n \rangle が e−iUn(n−1)t/2\mathrm{e}^{-\mathrm{i} Un(n-1)t/2} で回る。厳密な期待値は ⟨a^⟩(t)=αexp⁡[−∣α∣2(1−e−iUt)]\langle \hat a \rangle (t)=\alpha\exp\bigl[-\lvert \alpha\rvert ^2(1-\mathrm{e}^{-\mathrm{i} Ut})\bigr] で、その大きさは短い時間で ∣α∣ e−∣α∣2U2t2/2\lvert \alpha\rvert \,\mathrm{e}^{-\lvert \alpha\rvert ^2U^2t^2/2} と減る。粒子数の広がり ∣α∣\lvert \alpha\rvert の分だけ回る速さがばらつくためで、因数分解は時間 1/(U∣α∣)1/(U\lvert \alpha\rvert ) 程度より短い間だけ正しい(時間 2π/U2\pi/U で元に戻る)。∣α∣2=N→∞\lvert \alpha\rvert ^2=N\to\infty で UNUN を一定に保つと、崩れる時間 1/(UN)1/(U\sqrt N) は N\sqrt N に比例して延びる。平均場の極限は、この時間が無限に延びる極限である。

(3) a^∣N⟩=N∣N−1⟩\hat a| N \rangle =\sqrt N| N-1 \rangle は ∣N⟩| N \rangle と直交するので ⟨N∣a^∣N⟩=0\langle N | \hat a | N \rangle =0。粒子数の決まった状態では ⟨a^⟩\langle \hat a \rangle がそもそも 00 で、位相は定まらない。コヒーレント状態は粒子数の異なる状態の重ね合わせで、その相対的な位相が α\alpha の偏角である。凝縮体の位相 S/ℏS/\hbar は、この相対位相を全粒子がそろって持つことを表す。2つの凝縮体の間の干渉縞のように位相の差を測る実験では、粒子数の決まった2つの凝縮体からでも、測定のたびにある位相差が選ばれて縞が現れる。

第7部第1章の流体との関係

第7部第1章の流体と、この章の流体の違いを並べる。

第1に、密度の意味が違う。第1章の ρ\rho は1粒子の位置の確率密度で、1回の測定で見えるのは1点だけだった。この章の nn は凝縮体の粒子数密度で、NN 個の粒子が同じ ψ\psi に入っているので、1回の撮影で nn の形がそのまま写る。N=1N=1、g=0g=0 と置けば2つは同じ式になる。

第2に、住む空間が違う。多粒子の波動関数を極形式に書けば流体は 3N3N 次元に住むが、凝縮体では積の試行関数 (7.3.3) によって全粒子が同じ3次元の関数に入るので、流体は3次元の空間に戻る。その代償が平均場の近似で、凝縮から外れる粒子の割合が 00 になる極限(第1節の札)でだけ正しい。

第3に、量子化条件の意味が違う。第1章では循環の量子化は「流体から波動関数を作り直すために足す要請」で、その必要は節の線があるときにだけ現れた。この章では節の線が渦の芯として実在し、その周りの循環が実験で測れる。第7部の構成の出発点の要請のうち、流体の言葉の外にあった条件が、凝縮体では最も目立つ性質になっている。

地図の上の位置と授業の形

図 7.3.4 に、この章の結果がほかの部とつながる所をまとめた。第10部第3章(第10部第3章)の前量子化では、相空間の上の線束の曲率の積分が 2πℏ2\pi\hbar の整数倍になる条件が出る。循環の量子化 (7.3.17) と磁束の量子化 (7.3.24) は、どちらも「位相 S/ℏS/\hbar が閉曲線を1周して 2π2\pi の整数倍だけ増える」という同じ条件の別の姿である。

図

図 7.3.4 この章の結果と同値性の地図。太枠がこの章。上の矢印が第二量子化からの導き方、左が第7部第1章の1粒子の流体の拡張、右の2本が第8部・第10部へのつながり、下の2本が荷電した凝縮体と平均場の正当化である。

学部の授業では、超流動は2流体模型(全体の密度を超流動成分と常流動成分に分け、超流動成分は渦なしで粘性がないとする)で、超伝導は London 方程式 js=−(nsq2/m∗)A\bm{j}_{\mathrm s}=-(n_{\mathrm s}q^2/m^*)\bm{A}(London ゲージ)で導入する。この章の形と並べると、超流動成分の速度は v=∇S/m\bm{v}=\nabla S/m で、渦なしの条件は「速度が位相の勾配である」ことから出る。London 方程式は、(7.3.24) の速度場 v=(∇S−qA)/m∗\bm{v}=(\nabla S-q\bm{A})/m^* に ∇S=0\nabla S=0 を代入して電流 js=qnsv\bm{j}_{\mathrm s}=qn_{\mathrm s}\bm{v} を書いたものである。授業で「超流動成分は渦なし」「London 方程式」と仮定として置く2つは、この章ではどちらも巨視的波動関数の位相から導かれる。2流体模型の常流動成分は、第3節の Bogoliubov の励起が熱的に詰まったもので、有限温度の凝縮体の扱いは第7部の流体の外にある。

木下の勉強ステータス

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

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

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