超流動とボース＝アインシュタイン凝縮：絶対零度近傍で現れる巨視的量子現象の物理
現代物理学において、絶対零度（0 K = -273.15 ℃）近傍の極低温世界は、我々の日常的な直感を大きく裏切る驚異的な現象の宝庫である。その中でも「超流動（Superfluidity）」と「ボース＝アインシュタイン凝縮（Bose-Einstein Condensation, BEC）」は、ミクロな量子力学的性質がマクロなスケールで直接的に観測される「巨視的量子現象」の代表格として、数多くの物理学者を魅了し続けてきた。
粘性が完全にゼロになり、容器の壁を自発的に這い上がる液体ヘリウム。光を当てると細いノズルから液体が噴き出す熱機械効果（噴水効果）。そして、何兆個もの原子が全く同じ量子状態に落ち込み、単一の巨大な物質波として振る舞うBEC。本稿では、これらの現象がいかにして発見され、理論的に解明されてきたのか、その歴史的背景から高度な理論的骨格に至るまで、極めて詳細かつ学術的な視点から深く掘り下げていく。
第1章：絶対零度への挑戦とヘリウム液化の歴史
1.1 低温物理学の黎明と永久気体の液化
19世紀後半、物理学の大きなテーマの一つは「すべての気体は冷却と加圧によって液化できるか」という問題であった。マイケル・ファラデーらの努力により多くの気体が液化されたが、酸素、窒素、水素、そしてヘリウムなどは「永久気体」と呼ばれ、液化は極めて困難とされていた。
しかし、熱力学の発展とジュール＝トムソン効果を用いた冷却技術の進歩により、1877年に酸素が液化され、1898年にはジェームズ・デュワーによって水素の液化（沸点約 20 K）が達成された。残る最後の未踏の気体、それがヘリウムであった。
1.2 ヘリウムの液化とカメルリング・オネスの栄光
オランダのライデン大学に低温研究所を設立したヘイケ・カメルリング・オネス（Heike Kamerlingh Onnes）は、入念な準備と巨大な液化装置の構築を行った。1908年7月10日、彼は液体水素を予冷に用い、さらにジュール＝トムソン膨張を繰り返すことで、ついにヘリウムの液化に成功した。液体ヘリウムの沸点は 4.2 K であり、人類は絶対零度まであとわずか4度という極低温の世界への扉を開いたのである。この業績により、カメルリング・オネスは1913年にノーベル物理学賞を受賞した。
液体ヘリウムの生成は、後に1911年の超伝導の発見（水銀の電気抵抗が 4.2 K でゼロになる現象）へと直結するが、液体ヘリウムそのものの持つ異常な性質が明らかになるまでには、さらに数十年を要した。
1.3 ラムダ転移と超流動の発見
1920年代から1930年代にかけて、ウィレム・ヘンドリック・ケーソンらにより、液体ヘリウムをさらに冷却していくと、約 2.17 K において比熱が対数的な発散を示すことが見出された。この比熱の温度変化のグラフがギリシャ文字の「$\lambda$（ラムダ）」に似ていることから、この相転移は「$\lambda$（ラムダ）転移」と名付けられた。2.17 K（$T_\lambda$）より高温の相を「ヘリウムI」、低温の相を「ヘリウムII」と呼ぶ。
ヘリウムIIは、通常の液体とは全く異なる異常な振る舞いを示した。最も劇的な発見は1938年に相次いで報告された。モスクワのピョートル・カピッツァ（Pyotr Kapitza）、およびケンブリッジ大学のジョン・アレン（John F. Allen）とドン・マイズナー（Don Misener）は独立に、ヘリウムIIが極めて細い毛細管を抵抗なく（粘性ゼロで）流れ抜けることを発見したのである。
カピッツァはこの異常な流動性を、超伝導（Superconductivity）からの類推で「超流動（Superfluidity）」と命名した。カピッツァが行った実験では、ガラスの円盤を向かい合わせて微小な隙間（数ミクロン）を作り、そこを液体ヘリウムが流れる際の粘性を測定した。ヘリウムIでは通常の粘性流体としての抵抗が観測されたが、$T_\lambda$ を下回った瞬間に流出速度が劇的に増大し、見かけ上の粘性係数が少なくとも通常の $10^4$ 分の1以下に急減することが確認された。
さらに奇妙なことに、回転粘度計を用いた測定（円柱を液体中で回転させる実験）では、ヘリウムIIは依然として有限の粘性を持っているかのように振る舞うという矛盾した結果が得られた。この「細管では粘性ゼロ、回転体には粘性抵抗を及ぼす」という一見不可解な現象は、後の「二流体モデル」によって見事に解決されることとなる。
第2章：ボース＝アインシュタイン凝縮（BEC）の理論的骨格
液体ヘリウムの超流動をミクロな視点から理解する鍵は、量子統計力学にある。ヘリウム4原子は陽子2個、中性子2個、電子2個から構成され、そのスピンの総和は整数（0）となる。すなわち、ヘリウム4は「ボース粒子（ボソン）」である。ボース粒子の集団が絶対零度に向かって冷却されるとき何が起きるのか、その理論的骨格を導出する。
2.1 ボース統計とフェルミ統計
量子力学において、同種粒子は原理的に区別できない。粒子の入れ替えに対する波動関数の対称性から、粒子は2つのカテゴリーに大別される。
- ボース粒子（ボソン）：スピンが整数の粒子。粒子の入れ替えに対して波動関数は対称（符号が変わらない）。パウリの排他律に従わず、同一の量子状態に複数の粒子が同時に入ることができる。
- フェルミ粒子（フェルミオン）：スピンが半整数の粒子。粒子の入れ替えに対して波動関数は反対称（符号が変わる）。パウリの排他律に従い、同一の量子状態には1つの粒子しか入ることができない。
ここで、粒子数が非負であるためには、すべての状態 $i$ において $\epsilon_i - \mu > 0$ でなければならない。基底状態のエネルギーを $\epsilon_0 = 0$ とおくと、化学ポテンシャルは常に $\mu \leq 0$ を満たす必要がある。
2.2 熱的ド・ブロイ波長と巨視的波動関数の出現
$$ \lambda_{dB} = \sqrt{\frac{2\pi \hbar^2}{m k_B T}} $$高温では $\lambda_{dB}$ は非常に短く、粒子間の平均距離 $d = (V/N)^{1/3}$ よりもはるかに小さい。このとき、粒子は古典的なビリヤード球のように振る舞い、マクスウェル＝ボルツマン統計で十分に近似できる。 しかし、温度 $T$ を下げていくと、$\lambda_{dB}$ は次第に長くなる。そして、ある臨界温度 $T_c$ において、波長 $\lambda_{dB}$ が粒子間距離 $d$ と同程度（$\lambda_{dB} \sim d$）になる。このとき、個々の粒子の波動関数が空間的に重なり合い始め、量子力学的な干渉効果がマクロなスケールで顕在化する。
さらに温度を下げると、多数のボース粒子が最低エネルギー状態（基底状態）に雪崩を打って落ち込む。この状態では、無数の粒子が完全に同一の位相を持った単一の巨大な物質波、すなわち「巨視的波動関数（Macroscopic Wavefunction） $\Psi(\mathbf{r}) = \sqrt{n_0(\mathbf{r})} e^{i\theta(\mathbf{r})}$」として振る舞うようになる。これがボース＝アインシュタイン凝縮（BEC）である。
2.3 凝縮転移温度 $T_c$ の数学的導出
$$ N = \sum_i \frac{1}{e^{\beta(\epsilon_i - \mu)} - 1} $$$$ D(\epsilon) = \frac{V}{(2\pi)^2} \left( \frac{2m}{\hbar^2} \right)^{3/2} \sqrt{\epsilon} $$$$ N_{ex} = \int_0^\infty D(\epsilon) \frac{1}{e^{\beta(\epsilon - \mu)} - 1} d\epsilon $$$$ N_{max} = \frac{V}{(2\pi)^2} \left( \frac{2m}{\hbar^2} \right)^{3/2} \int_0^\infty \frac{\sqrt{\epsilon}}{e^{\beta \epsilon} - 1} d\epsilon $$$$ \int_0^\infty \frac{\sqrt{\epsilon}}{e^{\beta \epsilon} - 1} d\epsilon = (k_B T)^{3/2} \int_0^\infty \frac{x^{1/2}}{e^x - 1} dx = (k_B T)^{3/2} \Gamma(3/2) \zeta(3/2) $$$$ N_{max} = V \left( \frac{m k_B T}{2\pi \hbar^2} \right)^{3/2} \zeta(3/2) = \frac{V}{\lambda_{dB}^3} \zeta(3/2) $$もし全粒子数 $N$ が $N_{max}$ を超えるとどうなるか。余った粒子 $N_0 = N - N_{max}$ は、積分（$\epsilon > 0$）には含まれないエネルギー $\epsilon = 0$ の基底状態に凝縮するしかない。これがBECの数学的メカニズムである。
凝縮が始まる臨界温度 $T_c$ は、$N = N_{max}(T_c)$ となる温度として定義される。これを $T_c$ について解くと、 $$ T_c = \frac{2\pi \hbar^2}{m k_B} \left( \frac{n}{\zeta(3/2)} \right)^{2/3} $$ （ただし $n = N/V$ は数密度） この式は、相互作用のない理想気体における厳密なBEC転移温度を与える。液体ヘリウムの密度を代入すると $T_c \approx 3.1 K$ となり、実際の $\lambda$転移温度 2.17 K に近い値となる。この一致は、ヘリウムIIの超流動が本質的にBECに関連した現象であることを強く示唆している。ただし、液体ヘリウムは粒子間の相互作用が非常に強い「強相関量子液体」であるため、理想気体のモデルからずれが生じているのである。
第3章：ティサとランダウの二流体モデル
超流動ヘリウムが示す「細管では粘性ゼロで流れるが、回転体には粘性抵抗を示す」という奇妙な性質を説明するために、1938年にラースロ・ティサ（László Tisza）は画期的なフェノメノロジー（現象論）を提唱した。その後、1941年にレフ・ランダウ（Lev Landau）が量子力学的な微視的基礎に基づくより洗練された理論として完成させたのが「二流体モデル（Two-Fluid Model）」である。
3.1 常流動成分と超流動成分の共存
$$ \rho = \rho_n + \rho_s $$$$ \mathbf{j} = \rho_n \mathbf{v}_n + \rho_s \mathbf{v}_s $$ここで、
- 常流動成分（Normal fluid component, $\rho_n$）：有限の粘性を持ち、エントロピーを運ぶ成分。熱的な励起状態にある粒子集団に対応する。
- 超流動成分（Superfluid component, $\rho_s$）：粘性が完全にゼロで、エントロピーを持たない（絶対零度と同じ状態にある）成分。巨視的な量子状態（BEC状態）にある粒子集団に対応する。
絶対零度（$T = 0$）ではすべてが超流動成分（$\rho_s = \rho, \rho_n = 0$）であるが、温度が上がるにつれて常流動成分の割合が増え、$T_\lambda$ で超流動成分が消失する（$\rho_s = 0, \rho_n = \rho$）。
このモデルを用いると、先の粘性の矛盾は見事に解決される。極細管を通る実験では、粘性を持つ常流動成分は管の壁との摩擦でせき止められるが、粘性のない超流動成分のみが管をスリ抜けるため、見かけ上の粘性はゼロになる。一方、バルク液体中で円柱を回転させる実験では、円柱は常流動成分と摩擦を起こすため、有限の粘性が観測されるのである。
3.2 熱機械効果（噴水効果）とメカノカロリック効果
超流動成分が「エントロピーを持たない」という性質は、熱と物質の流れが直接的に結合する驚くべき現象を引き起こす。その代表が「熱機械効果（Thermomechanical effect）」、または「噴水効果（Fountain effect）」である。
$$ \Delta P = \rho S \Delta T $$（$S$ は単位質量あたりのエントロピー）
逆に、スーパーリークを通じてヘリウムIIを押し出すと、流れ出てくるのはエントロピーを持たない超流動成分のみであるため、残された液体の温度が上昇する。これをメカノカロリック効果と呼ぶ。
3.3 第二音波（Second Sound）の物理
通常の流体において「音波」とは密度の疎密波（圧力波）である。しかしヘリウムIIの中では、二つの独立した速度場 $\mathbf{v}_n$ と $\mathbf{v}_s$ が存在するため、2種類の波が伝播できる。
- 第一音波（First Sound）：常流動成分と超流動成分が同じ位相で振動する波。通常の密度波（圧力波）に相当する。
- 第二音波（Second Sound）：常流動成分と超流動成分が逆位相で（互いに逆向きに）振動する波。密度は一定に保たれたまま、常流動成分（エントロピーの塊）と超流動成分の割合が空間的に振動するため、これは事実上「温度の波（熱波）」として伝播する。
通常の物質では熱は拡散現象として伝わるが、ヘリウムII中では熱が減衰のない「波」として音速（第二音速）で伝播する。この現象もランダウによって予言され、後に実験的に実証されたことで、二流体モデルの正しさが決定的に証明された。
第4章：ランダウの超流動臨界速度と素励起スペクトル
ボース＝アインシュタイン凝縮の理論は、理想気体を前提としていたため、ヘリウム原子間に働く強い相互作用を説明できなかった。相互作用を持たない理想ボース気体では、実は超流動性は現れない（不純物散乱によって簡単に励起されてしまう）。レフ・ランダウは、液体ヘリウムを「絶対零度の基底状態にある背景液」とその上に乗る「素励起（Elementary Excitations）の気体」として定式化した。
4.1 フォノンとロトン：エネルギー・運動量分散曲線
ランダウは、ヘリウムIIの熱的な性質を説明するために、素励起のエネルギー $\epsilon$ と運動量 $p$ の関係（分散関係 $\epsilon(p)$）を以下のように仮定した。
- $$ \epsilon(p) = c p $$
（$c$ は第一音速）
- $$ \epsilon(p) = \Delta + \frac{(p - p_0)^2}{2\mu} $$
ここで、$\Delta$ はロトン・ギャップ・エネルギー、$p_0$ は極小値をとる運動量、$\mu$ はロトンの有効質量である。
この奇妙な分散曲線は、後に中性子非弾性散乱実験によって極めて正確に実測され、ランダウの天才的な直感が見事に証明された。
4.2 臨界速度 $v_c$ による摩擦ゼロの散乱不可能証明
ランダウの最大の功績は、この分散関係を用いて「なぜヘリウムIIは壁との摩擦なしに流れることができるのか」という超流動の根本的なメカニズムを力学的に証明したことにある（ランダウの臨界速度の定理）。
$$ \Delta E = \mathbf{v} \cdot \mathbf{p} = \epsilon(p) $$$$ v \geq \frac{\epsilon(p)}{p} $$$$ v_c = \min \left[ \frac{\epsilon(p)}{p} \right] $$もし液体ヘリウムが理想ボース気体であれば、自由粒子の分散関係 $\epsilon(p) = p^2 / 2m$ を持つため、$v_c = \min [p/2m] = 0$ となり、どんなに遅く動いても励起が起きてしまい超流動にはならない。 しかし、現実の液体ヘリウムでは相互作用により分散曲線が線形（$\epsilon = c p$）から始まるため、原点付近での傾きは $c$（音速、約 240 m/s）となる。分散曲線の形状全体から最小値を求めると、ロトン・ミニマムに向かう接線の傾きが $v_c$ を与え、その値は約 60 m/s となる。 （※実際の実験で観測される臨界速度はこれよりはるかに小さい数 cm/s 程度であるが、これは後にリチャード・ファインマンらによって「量子渦」の生成メカニズムとして解明されることになる。）
このように、系のエネルギー分散関係に「ギャップ」や「有限の傾き」が存在することが、超流動性（摩擦ゼロでの流動）を保証する絶対的な条件なのである。
第5章：量子渦（Quantized Vortex）とトポロジカル欠陥
ランダウの臨界速度理論は超流動の微視的起源を説明したが、実験的に観測される臨界速度が理論値よりはるかに低い理由は謎のままであった。この謎を解き明かし、超流動の巨視的波動関数の位相（位相幾何学、トポロジー）と結びつけたのが、ラルス・オンサーガー（Lars Onsager）とリチャード・ファインマン（Richard Feynman）による「量子渦」の理論である。
5.1 巨視的波動関数の位相と循環の量子化
$$ \mathbf{v}_s = \frac{\hbar}{m} \nabla \theta $$$$ \kappa = \oint_C \mathbf{v}_s \cdot d\mathbf{r} = \frac{\hbar}{m} \oint_C \nabla \theta \cdot d\mathbf{r} $$$$ \Delta\theta = 2\pi n \quad (n = 0, \pm 1, \pm 2, \dots) $$$$ \kappa = \frac{\hbar}{m} (2\pi n) = n \frac{h}{m} $$すなわち、超流動体における循環は連続的な値をとることができず、$h/m$（プランク定数 $h$ をヘリウム原子の質量 $m$ で割った値）を単位として完全に「量子化」されているのである。これが「循環の量子化」である。
5.2 回転する液体ヘリウム中の渦格子とトポロジカル欠陥
通常の粘性流体を入れた容器を回転させると、流体は剛体のように全体が回転し、放物面を描く。しかし、超流動ヘリウム（の超流動成分）は $\nabla \times \mathbf{v}_s = (\hbar/m) \nabla \times \nabla \theta = 0$、すなわち非回転（渦なし）の性質を持つため、容器をゆっくり回しても中身は一緒に回らない。 ところが、回転速度を徐々に上げていくと、ある臨界角速度を超えた瞬間に、超流動体の中に無数の「量子渦（Quantized Vortex）」の線（渦糸）が侵入してくる。
量子渦の中心（コア）では位相が特異点となり定義できないため、超流動密度 $n_0$ がゼロになる（つまり常流動になるか真空になる）。コアの周囲を流体が量子化された循環 $\kappa = h/m$ （通常は $n=1$ の最小状態がエネルギー的に安定）で回転している。 回転速度をさらに上げると、量子渦の数が増加し、それらが互いに反発し合って規則正しい三角格子（アブリコソフ格子に類似した渦格子）を形成する。巨視的に見れば、この無数の量子渦の集合体が、あたかも流体全体が剛体回転しているかのように振る舞うのである。
ランダウの臨界速度の不一致も、この量子渦によって説明された。毛細管の壁を流れる際、非常に小さな流速であっても、壁の微小な凹凸で量子渦が生成される。この量子渦の生成に必要なエネルギーが、ロトン励起よりもはるかに低いため、実験的な臨界速度はランダウの理論値より極端に小さくなっていたのである。量子渦は超流動を破壊する「トポロジカル欠陥」として機能する。
第6章：レーザー冷却による気体BECの実現と現代物理への展開
液体ヘリウムは粒子間相互作用が極めて強いため、凝縮率は絶対零度でも約10%程度にとどまり、純粋な理想気体のBECからは程遠いものであった。真の意味での「希薄気体におけるBEC」を実現するためには、原子間相互作用が無視できるほど希薄な気体を、絶対零度近くまで冷却する必要があった。
6.1 アルカリ原子気体のBEC達成
1980年代から1990年代にかけて、レーザー冷却（ドップラー冷却、シジフォス冷却）と磁気トラップ、さらに蒸発冷却という革新的な技術が開発され、原子気体をナノケルビン（$10^{-9}$ K）の領域まで冷却することが可能となった。 そして1995年、コロラド大学（JILA）のエリック・コーネル（Eric Cornell）とカール・ワイマン（Carl Wieman）はルビジウム（$^{87}$Rb）原子気体を用いて、さらに同年、MITのヴォルフガング・ケッテルレ（Wolfgang Ketterle）はナトリウム（$^{23}$Na）原子気体を用いて、人類史上初めて希薄アルカリ原子気体におけるボース＝アインシュタイン凝縮の直接観測に成功した。彼らはこの功績により2001年にノーベル物理学賞を受賞した。
ケッテルレのグループが行った実験では、2つの独立したBECを空間的に干渉させることで、明確な「干渉縞」が観測された。これは、巨視的なスケールの物質塊が単一の物質波として振る舞っていることを示す決定的な証拠であった。
6.2 光格子中の超流動・モット絶縁体転移
現代のBEC研究は、単なる現象の確認を越え、「量子シミュレーター」としての応用へと発展している。対向するレーザー光を干渉させて作る周期的なポテンシャル（光格子、Optical Lattice）の中にBECを導入すると、原子は固体の結晶格子中の電子のように振る舞う。 レーザーの強度を調整することで、原子間の相互作用とホッピング（隣のサイトへのトンネル効果）の比率を自在に制御できる。相互作用が弱いときは原子が系全体を位相の揃った状態で自由に動き回る「超流動状態」となるが、相互作用を強くしていくと、各格子点に原子が1つずつ局在し、マクロな位相のコヒーレンスが完全に破壊される「モット絶縁体（Mott Insulator）状態」へと量子相転移を起こす。 この超流動-モット絶縁体転移（2002年、Greinerらにより観測）は、強相関電子系（高温超伝導体など）のハバードモデルを理想的な人工系で精密に検証する道を切り開いた。
6.3 超伝導（クーパー対のBEC）との統一的理解
最後に、超流動と超伝導の深い関係に触れておこう。超伝導は、金属中の伝導電子（フェルミ粒子）が、フォノン（格子振動）を媒介とした引力によって2個ずつペア（クーパー対）を形成し、そのペアがボース粒子のように振る舞ってBECを起こす現象と解釈できる（BCS理論）。 近年では、フェルミ原子気体を用いて、原子間の相互作用を磁場（フェッシュバッハ共鳴）によって制御し、分子のBEC（BECリミット）からクーパー対のBCS状態（BCSリミット）へと連続的に状態を変化させる「BCS-BECクロスオーバー」の研究が盛んに行われている。
超流動とBECは、ミクロな量子力学の支配する法則が、日常的なマクロのスケールにまで拡張された究極の物理現象である。粘性ゼロの流動、循環の量子化、物質波の干渉。これらはすべて、自然界が絶対零度という極限状態で我々に見せる、宇宙の最も根本的で純粋な姿に他ならないのである。