数理解析研究所講究録 1413 巻 2005 年 231-238

231

Gross-Pitaevskii 方程式による超流動乱流のエネ
ルギースペクトル

大阪市立大学理学部物理学科 小林 未知数 (Michikazu Kobayashi)
Faculty of Science,

Osaka City University

自然界における乱流現象を要素還元的に理解すると い う視点におて、超流動乱流が近年非常に注目を浴び、活発に議論されて る。超流中の量子渦によって引き起こされる現象で動乱流とは超流動液体$\mathrm{A}\backslash$$\mathrm{A}\backslash$$4\mathrm{H}\mathrm{e}$あり、低温物理学において最も重要なトピツクの 1 つである [1]。液体原子力 Bose-Einstein
以下において、 構威要素であるは$\check{\backslash }$$4\mathrm{H}\mathrm{e}$$4\mathrm{H}\mathrm{e}$$2.17\mathrm{K}$

凝縮を起こし、粘性が消失する超流動状態, となる。 この超流動現象 1 ま、流体全体が粘性のある常流体成分と、粘性のな ‘超流体成分力\supset ら成る$\mathrm{A}$という二流体モデルにおいて記述され、 これを最も特徴づして、$\#\mathrm{e}$る現象と両或分がお互いに反対方向に流れることにより非常に大きな熱伝導度を引き起こす熱対向流がある。 両或分の相対速度力 ある値を超えると、 この熱対向流に散逸機構が生じ、 超流動乱流状態となる [1]。
$\backslash \theta$を持つ量子渦力 タングノレとなっ$\backslash \theta$

これは量子化された循環$\kappa=\oint v\mathrm{x}\mathrm{d}s$た状態であり ( は流体の速度場、$v$$\oint \mathrm{d}s$は渦回りの線積分) 、量子渦と常流体威分との間の相互摩擦力 (mutual ffiction) によって散逸力 \S |$\check{\backslash }$き起こされる [2]。

超流動乱流は発見されて以来すつとこの熱対向流を用 V‘て研究されてきたが [3]、熱対向流は超流動に固有の現象てあるため、通常の古典

232流体との対応を持たす、 それゆえに超流動乱流と古典乱流との関係は長い間謎のままであった。 このような背景のもとで近年、 Maurer 達が2 枚の回転円盤中で$[4]_{\text{、}}$

Stalp 達が振動格子中で [5] 超流動乱流を実現することに成功した。彼らの方法はまさに古典流体において乱流を作る方法と同じであり、 古典乱流との対応が可能になったという点において超流動乱流の研究は新しい局面に入ったといえよう。 彼らは IK以上という比較的高温領域の超流動乱流において、 そのエネルギースペクトルが古典乱流の最も重要な統計則てある Kolmogorov 則に従い
[6]、 超流動乱流が古典乱流との類似性を持つことを見出した。

Kolmogorov 則は発達した一様等方定常な非圧縮性古典乱流において成り立つ統計則である [6]。大きなスケールであるエネルギー注入領域からエネルギーが系に注入されると、 そのエネルギーのスケールは慣性領域にて小さくなってゆく。 慣性領域ではエネルギーが散逸されることなく、 系の詳細に依存しないスケール普遍性を持ち、エネルギースペクトノレ $E(k)$ が Kolmogorov 貝リ(1)$E(k)=C\epsilon^{2/3}k^{-5/3}$に従う。 ここで $E(k)$ は $E= \int \mathrm{d}kE(k)$ で定義されるエネノレギースペ
クトル、 $E$ は単位質量あたりの運動エネルギー、$k$は Fourier 変換の波数てある。粘性が有効になるエネルギー散逸領域においてエネルギーは散逸率 て散逸するが、 これは慣性領域におけるエネルギー流量に$\epsilon$等しい。 Kolmogorov 定数 $C$ はオーダー 1 の無次元量てある。

Maurer 達や Stalp 達が発見した超流動乱流と古典乱流の類似性は、古典流体のように振る舞う常流体成分と超流動或分との相互摩擦力を介してのカップリングによって両者が一緒になって古典乱流のように振る舞う、 という考えによって理解されている [7]。 そうすると常流233
体成分がほとんど存在しない 1K 以下といったような極低温におも超流動乱流は古典乱流との類似性を持つのであろう力\supset 、 とて$\mathrm{A}\backslash$

う大き$\mathrm{A}\backslash$

な疑問が生じる。 この疑問に対して、 以下のように考える。 古典舌 1 流における Kolmogorov 則は通常、 大きな渦が\Delta ゝさな渦ゞと分裂してゆくという Richardson カスケードという描像で理解される [6]。

し力 し$>$

古典流体では渦そのものの定義が、粘性拡散のために不明瞭であるため、 Richardson カスケードが本当にあるかどうか{は分力 1らな$\mathrm{A}\backslash _{\mathrm{Q}}-$方、 超流動乱流における量子渦は明確に定義できる位相欠陥である。

Bose-Einstein 凝縮によって形成される巨視的な量子力学的波動関数$\Phi(x, t)=\sqrt{\rho(x,t)}\mathrm{e}^{\mathrm{i}\phi(X,t)}$のダイナミクスは Gross-Pitaevskii 方程式$[8, 9]$(2)$. \frac{\partial}{\partial t}\Phi(x, t)=[\nabla^{2}-\mu+g|\Phi(x, t)|^{2}]\Phi(oe, t)$によって記述され、超流動速度場 $v(oe, t)$ は $v(oe, t)=2\nabla\phi(x, t)$ で与えられる。 ここで は化学ポテンシャル、 は粒子間相互作用の結合定$g$$\mu$数、$\rho$

は凝縮体密度、$\phi$結領域では存在せす、は凝縮体の位相である。 渦度 rot$\Phi$$v$

はの位相欠陥でのみ値を持つ。 つまり$\Phi$の単連$\Phi$の位相欠陥が量子渦の定義そのものになる。量子渦の回りてま循環力
子化され、芯のサイズは回復長$\xi=1/\sqrt{g\rho}$$\backslash 4*\pi$に量で与えられる。 高温で{は量

子渦は相互摩擦力を介して減衰するが、極低温では渦の再結合時に起こる圧縮性の素励起放出や回復長 程度まで小さくなった渦輪の消失
よを通してのみ減衰する [10]。いすれも場合にお ても減衰{ま回復長$\xi$$\mathrm{A}\backslash$$\xi$

りも短いスケールで起こり、長いスケールでは量子渦 t は安定に存在することができる。 いわば量子渦は古典流体中の渦に付随する粘性散逸といったようなよけいな自由度を取り除いた形となっており、 これを構成要素とする超流動乱流は、 Kolmogorov 則と Richardson カスケー

234

ドの関係を明らかにすることのできるプロトタイプであるかもしれない。 もしこれが真実であるなら、 極低温の超流動乱流も Kolmogorov則を示すはずである。 本研究では Gross-Pitaevskii 方程式を数値的に解くことにより、 この問題を理論的に調べた。

本研究以前にも Nore 達が Gross-Pitaevskii 方程式の数値シミュレーションを用いて、超流動乱流のエネルギースペクトルを計算して$\mathrm{A}\backslash$る[11]。 Gross-Pitaevskii は圧縮性流体の方程式であるので、彼らは運動エネルギーを圧縮成分と非圧縮戒分に分けることによって、非圧縮成分のスペクトルが一時的に Kolmogorov 則に従うことを発見した。 しかし渦の再結合時に渦のエネルギーは回復長 上りも短い圧縮性の短波長$\xi$素励起へと転化し、渦のカスケードに影響を与える [10, 12, 13, 14, 15]。
よって Kolmogorov 則とは長時間は一致しない。

本研究では、 回復長$\xi$よりも短いスケールでのみ働く散逸項を導入し、このような素励起を散逸させることで素励起から渦へのエネルギーの逆流を阻止する、 ということを行った。 Gross-Pitaevskii 方程式を高い精度で解くために、周期境界条件下でのスペクトル法を用いる [16]。
そのために Fourier 変換された Gross-Pitaevskii 方程式$\mathrm{i}\frac{\partial}{\partial t}\Phi(k, t)=[k^{2}-\mu]\Phi(k, t)$$+ \frac{g}{V^{2}}\sum\Phi(k_{1}, t)\Phi^{*}(k_{2}, t)\Phi(k-k_{1}+k_{2}, t)$(3)$k_{1},k_{2}$を解くことになるが、 ここで (3) 式の左辺の虚数単位 を$\mathrm{i}$で置き換えることにより散逸項$\gamma_{0}\theta(k-2\pi/\xi)$の形をしており$\gamma(k)$($\mathrm{i}arrow[\mathrm{i}-\gamma(k)]$を導入する。 この散逸項は $\gamma(k)=$$\theta(x)$は階段関数) 、回復長いスケールてのみ値を持つ。 よってこの散逸項は、 回復長$\xi$よりも短$\xi$よりも短い構造を持たない量子渦を散逸させることなく、 渦の再結合によって生じた短波長素励起のみを散逸させる。 このようにして量子渦のみに

235よって作られる超流動乱流を理論的に議論することが可能 [こなる。
一様で等方な乱流状態を得るために、 初期状態としては一様な凝縮体密度$\rho_{0}$と空間的にランダムな位相$\phi_{0}(x)$から出発する。 ランダムな位相は図 1 に示されるように、 空間中に等間隔$\lambda$

毎にー

力 ら$\backslash$$\pi$$\pi$間での乱数を与え、 これを滑らかにつなぐことで得られる。 初期の速度場

(b)

(a)
図 1: ランダムな位相を作る方法 (a) と$\phi_{0}(x,y,0)$$\mathrm{x}\mathrm{y}$平面におけるランダムな位相の断面図の一例 (b)$v(x, t=0)=2\nabla\phi_{0}(x)$ はランダムなので、系{はすぐに多くの渦輪を作って一様等方な乱流となる。

数値計算のパラメーターとして、結合定数 g=l、回復長 \mbox{\boldmath $\xi$}=1、初期の密度 \rho 0=1、散逸の大きさ$\gamma_{0}=1$を用いる。周期境界条件にお積 $V=32^{3}$ の立方体を $256^{3}$ の格子点に分け、空間解像度をとし、 空間微分に関してはスペクトル法を用$\iota\backslash$$\#\mathrm{e}$る体$\triangle x=0.125$

る。 また波数の解像度

は $\triangle k=2\pi/32$ である。 時間微分に関しては解像度$\triangle t=1\cross 10^{-4}$において Runge-Kutta-Verner 法を用いる [16]。

得られた超流動乱流状態に関して、 運動エネルギーの非圧縮成分EU。におけるスペクトノレ$E_{\mathrm{k}\mathrm{i}\mathrm{n}}^{\mathrm{i}}(k)$を計算する。 ここで$E_{\mathrm{k}\mathrm{i}\mathrm{n}}^{\mathrm{i}}(k)1\mathrm{h}E_{\mathrm{k}\mathrm{i}\mathrm{n}}^{\mathrm{i}}=$

238$\int \mathrm{d}kE_{\mathrm{k}\mathrm{i}\mathrm{n}}^{\mathrm{i}}(k)$は$(|\Phi|\nabla\phi)$

によって得られ、において$E_{\mathrm{k}\mathrm{i}\mathrm{n}}^{\mathrm{i}}= \int \mathrm{d}x[(|\Phi|\nabla\phi)^{\mathrm{i}}]^{2}$$\mathrm{d}\mathrm{i}\mathrm{v}(|\Phi|\nabla\phi)^{\mathrm{i}}=0$が進むにつれスペクトル$E_{\mathrm{k}\mathrm{i}\mathrm{n}}^{\mathrm{i}}(k)$$k<2\pi/\xi$ において波数 $k$ のべき$\eta$の時間変化を図 2(a) に示す。

$(|\Phi|\nabla\phi)^{\mathrm{i}}$を満たす成分である。 時間発展

は、散逸項の効かない慣性領域$E_{\mathrm{k}\mathrm{i}\mathrm{n}}^{\mathrm{i}}(k)\propto k^{-\eta}$$t_{\sim}>4$$\Delta k<$を示すようになる。指数において超流動乱流が Kolmogorov

(b)

(a)
図 2: エネルギースペクトルであり、$E_{\mathrm{k}\mathrm{i}\mathrm{n}}^{\mathrm{i}}(k)$の指数$\eta$の時間変化。 (a) は散逸あり、 (b) は散逸なし。

則$E_{\mathrm{k}\mathrm{i}\mathrm{n}}^{\mathrm{i}}(k)\propto k^{-5/3}$を入れない場合のを満たしていることが分かろ。 比較のために散逸項$\eta$の時間変化を図$2(\mathrm{b})$に示す。 この場合、$t_{\sim}>7$において再ひ Kolmogorov から外れていくことが分かる。 これは短波長の圧縮性素励起が量子渦のダイナミクスに影響を与えているということを意味している。

次に系のエネルギー散逸率 を$\epsilon$$\epsilon=\partial E_{\mathrm{k}\mathrm{i}\mathrm{n}}^{\mathrm{i}}(k)/\partial t$から求めることにより、 エネルギースペクトルと Kolmogorov 則との定量的な比較を行った。 $t=6$ における量子渦のスナップショットを図
スペクトルを図$3(\mathrm{b})$$3(\mathrm{a})$に、エネルギーに示す。 量子渦の分布はほぼ一様て等方であり、発達した一様等方な乱流ができていることを支持している。 またその

237(b)(a)

図 3: $t=6$ における量子渦のスナップショット (a) と、 エネルギースペクトル (b)。

ときのエネルギースペクトルと Kolmogorov 則との定量的な一致は非常に良い (本研究では Kolmogorov 定数$\mathrm{C}$に関して、 およそ $C\simeq 0.3$が得られた) 。このようにして本研究において、 極低温での超流動乱流における量子渦のダイナミクスは、 短波長の圧縮性素励起の影響がないときに古典乱流と類似の統計を示すことが証明された。

これまて量子力学の物性への応用はほとんど固体に限定されてきた。

しかし本研究を通して、超流動流体力学が通常の古典流体の、特に乱流現象を量子渦を用いて要素還元的に理解できる可能性を持っていることが示された。 つまり超流動という概念が、流体に有する非常に多くの自由度を減らす可能性を持っているということである。 今後、 この量子渦という要素還元的な見方が、 流体のより深い理解に貢献するてあろうことが期待される。

参考文献[1] R. J. Donnelly, Quantized vortices in helium $II$ (Cambridge University Press,Cambridge, 1991).

238[2] W. F. Vinen, Proc.

$\mathrm{R}\mathrm{o}\mathrm{y}.$Soc. A 240, 114 (1957); ibid. A 240, 128 (1957); ibid.

A240, 493 (1957).
[3] J. T. Tough, in Progress in Low Temperature Physics Vol. VIII, edited by C. J.

Gorter (North-Holland, Ameterdam, 1955), P. 133.
[4] J. Maurer and P. Tabeling, Europhys. Lett. 43 (1), 29 (1998).
[5] S. R. Stalp, L. Skrbek, and R. J. Donnelly, Phys. Rev. Lett. 82, 4831 (1999).
[6] U. Risch, Turbulence (Cambridge University Press, Cambridge, 1995).

[7] W. F. Vinen, Phys. Rev. $\mathrm{B}61,1410(2000)$ .

[8] E. P. Gross, J. Math. Phys. 4, 195 (1963).
[9] L. P. Pitaevskii, Soviet Phys.-JETP 13, 451 (1961).
[10] M. Leadbeater, T. Winie&i, D. C. Samuels, C. F. Barenghi and C. S. Adams,Phys. Rev. Lett. 86, 1410 (2001).
[11] C. Nore, M. Abid, and M. E. Brachet, Phys. Rev. Lett. 78, 3896 (1997); Phys.

Fluids 9, 2644 (1997).
[12] S. Ogawa, M. Tsubota and Y. Hattori, J. Phys. Soc.

$\mathrm{J}\mathrm{p}\mathrm{n}$. $71,813(2002)$ .

[13] N. G. Berloff, Phys. Rev. A 69, 053601 (2004).
[14] N. G. Berloff and B. V. Svistunov, Phys. Rev. A 66, 013603 (2002).
[15] J. Koplik and H. Levine, Phys. Rev. Lett. 71, 1375 (1993) ;ibid. 76, 4745 (1996).

[16] W. H. Press, S. A. Teukolsky, W. T. Vetterling and B. P. Flamery, Numerical

Recipes in $Fo\hslash mn77$ (Cambridge University Press, Cambr.idge, 1985).