tecktom.com

この章のまとめ

2026年10月17日8分で読めますAI生成教科書『虎の巻』第7部 第3章 量子流体:渦の量子化と Gross–Pitaevskii 方程式の流体形

AI 生成・木下が未勉強

量子力学の様々な構成方法 第7部 第3章 §7.3.7 ★★★ GP 方程式の流体形、渦の量子化、Landau の臨界速度を振り返り、総合問題で輪の中の永久流を解く

多数のボース粒子が同じ状態に入った凝縮体は、粒子数密度に規格化した巨視的波動関数 ψ\psi で書け、ψ\psi は平均場の変分から GP 方程式に従う。GP 方程式を極形式に書くと、第7部第1章の Madelung 方程式に相互作用の圧力 Pint=gn2/2P_{\mathrm{int}}=gn^2/2 が加わった順圧流体の方程式になり、流体の言葉と GP 方程式は、第7部第1章と同じ量子化条件のもとで同値である。量子応力と相互作用の圧力の兼ね合いで癒し長 ξ=ℏ/2mgn\xi=\hbar/\sqrt{2mgn} が決まり、ξ\xi より長い距離では量子応力が落ちて Thomas–Fermi 近似になる。相互作用を 00 にすると、共通例題①の基底状態(μ=ℏω/2\mu=\hbar\omega/2)に戻る。

流体を線形化すると Bogoliubov の分散 ωk2=cs2k2+(ℏk2/2m)2\omega_k^2=c_{\mathrm s}^2k^2+(\hbar k^2/2m)^2 が出て、長波長は音波、短波長は量子応力の項から出る粒子の分散になる。Landau の判定条件から臨界速度は音速 csc_{\mathrm s} で、斥力の相互作用が無ければ超流動は起きない。循環は κ=h/m\kappa=h/m の整数倍に量子化され、その節の線が量子渦の芯になる。芯の大きさは ξ\xi の程度で、単位長さのエネルギーは (πℏ2n0/m)ln⁡(1.464R/ξ)(\pi\hbar^2n_0/m)\ln(1.464R/\xi)、巻き数の2乗に比例するので巻き数1の渦が安定である。回る容器の中では渦の密度が 2mΩ/h2m\Omega/h(Feynman の規則)になり、荷電した対の凝縮体では同じ量子化が磁束の量子 h/2eh/2e になる。

橋の側では、GP 方程式は場の演算子の Heisenberg 方程式の期待値をコヒーレント状態で因数分解したもので、その線形化が Bogoliubov の励起である。巨視的波動関数の位相を選ぶことは位相の回転の対称性を破ることで、第8部の非同値な表現の話につながる。循環と磁束の量子化は、第10部の前量子化の整数条件と同じ形をしている。

総合問題 輪の形の罠の中の永久流 ★★★

オリジナル半径 R0R_0 の細い輪の形の罠に、質量 mm の原子の凝縮体が、輪に沿って一様な密度 n0n_0 で入っている。輪の断面の大きさは ξ\xi より大きく、R0R_0 よりずっと小さいとする。凝縮体の位相を輪に沿った角度 φ\varphi で S/ℏ=ℓφS/\hbar=\ell\varphi と置く。数値は例題 7.3.1 の原子(cs=1.83 mm/sc_{\mathrm s}=1.83\,\mathrm{mm/s})と R0=20 μmR_0=20\,\mu\mathrm m を使う。

  1. (1) 流れの速さ vℓv_\ell と1粒子あたりの流れのエネルギーを ℓ\ell で書け。ℓ\ell が整数に限られる理由を述べ、v1v_1 の値を求めよ。(使う道具:この章の第4節・第7部第1章の輪の上の流れの例題)
  2. (2) 流れ ℓ\ell が、輪の壁の凹凸との摩擦で励起を作って減速するかどうかを Landau の判定条件で調べ、減速しない最大の ℓ\ell を求めよ。(使う道具:第3節)
  3. (3) 輪を角速度 Ω\Omega で回したとき、回る座標系で E−ΩLzE-\Omega L_z を最小にする ℓ\ell を求め、ℓ\ell と ℓ+1\ell+1 が入れ替わる Ω\Omega を求めよ。ℓ=0\ell=0 から 11 へ移る角速度の値を求めよ。(使う道具:第5節)
  4. (4) 同じ形の超伝導の輪を貫く磁束 Φ\Phi を変えたときの、対の凝縮体の巻き数の入れ替わりと (3) を対応させよ。(使う道具:第5節・第7部第1章の電磁場があるときの注意・解析力学から様々な理論分野を組み立てる 第5部第6章)
  5. (5) 相互作用 g→0g\to0 の極限で、(2) の結論がどう変わるかを述べ、凝縮していることと超流動であることの違いを説明せよ。(使う道具:第1節・第3節・第6節)
解答

(1) v=∇S/m\bm{v}=\nabla S/m で、輪に沿った長さは R0φR_0\varphi なので vℓ=ℏℓ/(mR0)v_\ell=\hbar\ell/(mR_0)。密度が一様なので量子ポテンシャルは 00 で、1粒子あたりの流れのエネルギーはmvℓ2/2=ℏ2ℓ2/(2mR02)mv_\ell^2/2=\hbar^2\ell^2/(2mR_0^2)。輪を1周する循環は 2πR0vℓ=ℓh/m2\pi R_0v_\ell=\ell h/m で、(7.3.17) から ℓ∈Z\ell\in\mathbb{Z} でなければならない。第7部第1章の例題で見たとおり、流体の方程式だけなら ℓ\ell は任意の実数でよく、ℓ\ell を整数に限るのは ψ∝eiℓφ\psi\propto\mathrm{e}^{\mathrm{i}\ell\varphi} が1価であるという条件である。この流れは輪の中心を囲むので、輪の穴(密度 00 の領域)が渦の芯の役をしている。数値は ℏ/(mR0)=1.0546×10−34/(1.443×10−25×2×10−5) m/s=3.65×10−5 m/s\hbar/(mR_0)=1.0546\times10^{-34}/(1.443\times10^{-25}\times2\times10^{-5})\,\mathrm{m/s}=3.65\times10^{-5}\,\mathrm{m/s}、v1=36.5 μm/sv_1=36.5\,\mu\mathrm{m/s} である。

(2) 命題 7.3.2 から、流れの速さが音速より遅ければ、壁との摩擦で励起を作ってもエネルギーは下がらず、流れは減速しない。条件 ℏℓ/(mR0)<cs\hbar\ell/(mR_0)<c_{\mathrm s} は ℓ<mcsR0/ℏ=1.83×10−3/3.65×10−5=50.1\ell<mc_{\mathrm s}R_0/\hbar=1.83\times10^{-3}/3.65\times10^{-5}=50.1 で、ℓ≤50\ell\le50 の流れは永久流として続く。ℓ\ell を1つ下げるには、輪の断面を渦が横切って位相の巻き数を1つ減らす必要があり、それにはエネルギーの山(渦の芯を作るエネルギー)を越えなければならない。Landau の条件はこの山を越えない範囲での安定性で、実際の寿命の見積もりにはこの山の高さも要る。

(3) 全粒子が ℓ\ell の状態にあるとき、全体で Lz=NℏℓL_z=N\hbar\ell、流れのエネルギーは E=Nℏ2ℓ2/(2mR02)E=N\hbar^2\ell^2/(2mR_0^2) なので

E−ΩLzN=ℏ2ℓ22mR02−ℏΩℓ=ℏ22mR02(ℓ−mR02Ωℏ)2−mR02Ω22.\frac{E-\Omega L_z}{N}=\frac{\hbar^2\ell^2}{2mR_0^2}-\hbar\Omega\ell=\frac{\hbar^2}{2mR_0^2}\Bigl(\ell-\frac{mR_0^2\Omega}{\hbar}\Bigr)^2-\frac{mR_0^2\Omega^2}{2}.

最小にする整数 ℓ\ell は mR02Ω/ℏmR_0^2\Omega/\hbar に最も近い整数で、ℓ\ell と ℓ+1\ell+1 の値が等しくなるのは ℏ2(2ℓ+1)/(2mR02)=ℏΩ\hbar^2(2\ell+1)/(2mR_0^2)=\hbar\Omega、つまり

Ωℓ=ℏmR02(ℓ+12)\Omega_\ell=\frac{\hbar}{mR_0^2}\Bigl(\ell+\frac12\Bigr)

である。ℏ/(mR02)=3.65×10−5/(2×10−5) s−1=1.83 rad/s\hbar/(mR_0^2)=3.65\times10^{-5}/(2\times10^{-5})\,\mathrm{s^{-1}}=1.83\,\mathrm{rad/s} なので、ℓ=0\ell=0 から 11 へ移るのは Ω0=0.91 rad/s\Omega_0=0.91\,\mathrm{rad/s}。回転を速くすると、凝縮体の角運動量は連続にでなく NℏN\hbar ずつ階段状に増える。Ω\Omega が Ω0\Omega_0 を越えたあと回転を止めても、(2) により ℓ=1\ell=1 の流れは残る。

(4) 超伝導の輪では (7.3.24) から、細い輪に沿った速度は m∗v=(2πℏℓ−qΦ)/(2πR0)m^*v=(2\pi\hbar\ell-q\Phi)/(2\pi R_0) で、運動エネルギーは (ℓ−qΦ/(2πℏ))2(\ell-q\Phi/(2\pi\hbar))^2 に比例する。∣q∣=2e\lvert q\rvert =2e なので、最小にする ℓ\ell は ±Φ/Φ0\pm\Phi/\Phi_0 に最も近い整数(符号は q=−2eq=-2e の向きの取り方による)で、ℓ\ell の入れ替わりは ∣Φ∣=(ℓ+12)Φ0\lvert \Phi\rvert =(\ell+\frac12)\Phi_0 で起き、エネルギーと電流は Φ\Phi について周期 Φ0=h/2e\Phi_0=h/2e で繰り返す。(3) と比べると、回転の「磁束」2πmΩR022\pi m\Omega R_0^2 を hh で割ったものが Φ/Φ0\Phi/\Phi_0 に当たり、2πmΩℓR02=(ℓ+12)h2\pi m\Omega_\ell R_0^2=(\ell+\frac12)h は Ωℓ=(ℏ/mR02)(ℓ+12)\Omega_\ell=(\hbar/mR_0^2)(\ell+\frac12) と同じ式である。回る座標系の mΩ×rm\bm{\Omega}\times\bm{r} が qAq\bm{A} の役をするので、中性の凝縮体を回すことと、荷電した凝縮体に磁束を通すことは、同じ量子化条件の2つの実現である。

(5) g→0g\to0 で cs=gn0/m→0c_{\mathrm s}=\sqrt{gn_0/m}\to0 となり、(2) の条件 vℓ<csv_\ell<c_{\mathrm s} を満たす ℓ≠0\ell\neq0 は無くなる。理想的なボース気体の凝縮体では、どんなに遅い流れも壁との摩擦で ε(p)=p2/(2m)\varepsilon (p)=p^2/(2m) の励起を作ってエネルギーを下げられるので、永久流は減衰する。第1節で見たように g=0g=0 の凝縮体は全粒子が1粒子の基底状態(罠が調和なら共通例題①の基底状態)に入っただけのもので、凝縮は起きている。超流動に要るのは、凝縮(巨視的波動関数があること)に加えて、相互作用が分散を音波に変えて臨界速度を有限にすることである。第6節の言葉では、Bogoliubov の励起の分散が原点で線形であることが超流動の条件で、これは凝縮のモードの外の量子揺らぎの性質である。

木下の勉強ステータス

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

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

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