デバイ遮蔽

プラズマとは、たくさんの荷電粒子がばらばらに動きながらも、全体としては一つの集団のように振る舞う状態です。その「集団としての振る舞い」の一番の入り口が、このページで扱うデバイ遮蔽（Debye shielding）です。デバイ遮蔽を理解すると、デバイ長という長さのものさしが手に入り、そこからプラズマ振動、準中性、そして「そもそもプラズマとは何か」という厳密な定義までを一本の道すじでたどれます。

まずは直感から（高校レベル）

プラズマの中に、正の電気を帯びた小さな玉を一つ置いたと想像してください。この玉のまわりには、もともとプラズマ中を飛び回っていた電子（負の電気）とイオン（正の電気）がたくさんいます。

正の玉は電子を引き寄せ、イオンを遠ざけます。すると玉のまわりには、負の電子が少し多く集まった「電子の雲」ができます。この雲は、玉の正の電気を外から見えにくくします。ちょうど、まぶしい電球のまわりを黒い網で覆うと、遠くから見た明るさがやわらぐのと同じです。これがデバイ遮蔽です。

電子はイオンよりずっと軽いので、すばやく動いて先に雲を作ります。その結果、少し離れた場所にいる別の粒子から見ると、玉の電気の影響はほとんど消えてしまいます。この「電気の影響が届く距離」の目安が、デバイ長（Debye length）と呼ばれる長さです。

ここで大事なのは、なぜ電気の影響が完全には消えず、ある距離までは残るのかという点です。もし粒子が冷たく（温度が低く）じっとしていれば、電子はきっちり玉に張りついて電気を完全に打ち消してしまいます。しかし実際の粒子は熱を持っていて、たえず動き回っています。この熱運動が雲をぼやけさせ、遮蔽を不完全にします。つまりデバイ長は、電気で引き寄せる力と、熱でばらけようとする勢いのつり合いで決まる長さなのです。温度が高いほど雲はぼやけてデバイ長は長くなり、粒子の密度が高いほど雲は濃くなってデバイ長は短くなります。

物理として理解する（学部レベル）

デバイ遮蔽を式で表しましょう。真空中では、電荷 $q$ のまわりの電位はクーロンポテンシャル

\phi_{\text{vac}}(r) = \frac{q}{4\pi\varepsilon_0 r}

に従い、距離 $r$ に反比例してゆっくり減衰します。ここで $\varepsilon_0$ は真空の誘電率です。ところがプラズマ中では、まわりの荷電粒子が再配置して電位を打ち消すため、湯川型（遮蔽クーロン型）のポテンシャル

\phi(r) = \frac{q}{4\pi\varepsilon_0 r}\, e^{-r/\lambda_D}

に変わります。指数関数 $e^{-r/\lambda_D}$ が効いて、距離が $\lambda_D$ を超えると電位は急速にゼロへ近づきます。この $\lambda_D$ がデバイ長です。

デバイ長は、電子温度 $T_e$ 、電子密度 $n_e$ を用いて

\lambda_D = \sqrt{\frac{\varepsilon_0 k_B T_e}{n_e e^2}}

と書けます。ここで $k_B$ はボルツマン定数、 $e$ は素電荷です。分子に温度 $T_e$ 、分母に密度 $n_e$ が来ていることが、先ほどの直感（温度が高いほど長く、密度が高いほど短い）と一致しています。

導出の筋道は次の通りです。電位 $\phi$ の中では、電子の密度はボルツマン分布 $n_e \exp(e\phi/k_B T_e)$ に従います。電位が小さい（ $e\phi \ll k_B T_e$ ）として指数関数を 1 次まで展開し、これをポアソン方程式 $\nabla^2 \phi = -\rho/\varepsilon_0$ に代入すると、 $\nabla^2 \phi = \phi/\lambda_D^2$ という形の方程式が得られます。この解が上の湯川型ポテンシャルであり、 $\lambda_D$ が自然に遮蔽の長さスケールとして現れます。

核融合プラズマ（温度 10 keV 程度、密度 $10^{20}\ \text{m}^{-3}$ 程度）では、デバイ長はおよそ 70 マイクロメートルです。この長さは、装置サイズ（メートル）よりはるかに小さく、逆に粒子間の平均距離（この密度では約 0.2 マイクロメートル）よりはるかに大きい、という二つの不等式を同時に満たします。この「装置より十分小さく、粒子間隔より十分大きい」という関係が、あとで述べるプラズマの定義の核心になります。

もう一つの重要量がプラズマ振動です。一様なプラズマで電子集団が全体としてわずかに $x$ だけずれると、正のイオンとの間に電荷分離が生じ、もとに戻そうとする電場（復元力）が働きます。電子はこの復元力で単振動を始め、その角周波数は電子プラズマ周波数（electron plasma frequency）

\omega_{pe} = \sqrt{\frac{n_e e^2}{\varepsilon_0 m_e}}

で与えられます。 $m_e$ は電子質量です。核融合プラズマでは周波数 $f_{pe} = \omega_{pe}/2\pi$ はおよそ 90 GHz になります。イオンは質量が大きいぶん応答が遅く、イオンプラズマ周波数は電子の場合の $\sqrt{m_e/m_i}$ 倍と小さくなります。水素イオンでは電子の約 43 分の 1 です。

デバイ長とプラズマ周波数は、電子の熱速度 $v_{th,e} = \sqrt{k_B T_e/m_e}$ を通じて $\lambda_D \sim v_{th,e}/\omega_{pe}$ という関係で結ばれています。これは「電子がプラズマ振動を 1 回する間に、熱運動でおよそデバイ長だけ移動する」ことを意味し、遮蔽と振動が同じ物理から生まれていることを示します。

理論を深める（大学院レベル）

デバイ球（半径 $\lambda_D$ の球）の中にどれだけ粒子が入っているかを表すのが、プラズマパラメータ（plasma parameter）です。デバイ球内の粒子数は

N_D = \frac{4}{3}\pi \lambda_D^3\, n_e

で定義され、 $N_D \gg 1$ が集団的振る舞いの条件です。核融合プラズマでは $N_D$ はおよそ $10^8$ に達します。 $N_D$ が大きいほど、多数の粒子が協調して遮蔽と振動を担うため、個々の 2 体クーロン衝突よりも集団効果が支配的になります。逆に $N_D \sim 1$ の系は「強結合プラズマ」と呼ばれ、統計的な取り扱いが破れて別の理論が必要になります。核融合プラズマは $N_D \gg 1$ の「弱結合・理想プラズマ」に分類されます。

以上を踏まえると、電離した気体がプラズマとして振る舞うための厳密な 3 条件を、次のように整理できます。

第一に、デバイ長が系のサイズ $L$ より十分小さいこと（ $\lambda_D \ll L$ ）。これにより、系の大部分で電荷が遮蔽され、準中性が成立します。

第二に、デバイ球内粒子数が多いこと（ $N_D \gg 1$ ）。これにより遮蔽が統計的な意味を持ち、集団効果が 2 体衝突を上回ります。

第三に、電子プラズマ周波数が電子と中性粒子の衝突周波数 $\nu$ より十分高いこと（ $\omega_{pe}\tau > 1$ 、 $\tau = 1/\nu$ ）。これは、電子がプラズマ振動を完了する前に中性粒子との衝突で乱されず、電磁気的な集団運動が気体的な衝突運動に勝る、という条件です。

これらは互いに独立ではなく、いずれも「電磁的な集団効果が、個別粒子の運動や熱擾乱に勝る」という一つの思想を、長さ・粒子数・時間の三つの側面から表したものです。

準中性（quasineutrality）は、デバイ長より大きなスケールでは電子密度とイオン密度がほぼ等しく（ $n_e \approx Z n_i$ ）、正味の電荷密度がほぼゼロと見なせる、という近似です。MHD（電磁流体力学）をはじめとする多くのプラズマ理論は、この準中性を出発点にしています。ただし準中性は万能ではなく、デバイ長程度の小さなスケール、高周波電場中、そしてプラズマと壁の境界では破れます。

準中性が破れる代表例がシース（sheath）です。プラズマが固体壁に接すると、動きの速い電子が先に壁へ流れ込んで壁を負に帯電させ、壁の手前にデバイ長の数倍程度の厚さを持つ正の空間電荷層ができます。これがシースで、その中では準中性が成り立ちません。シースには重要な条件があり、壁へ流れ込むイオンがイオン音速 $c_s = \sqrt{k_B T_e/m_i}$ 以上に加速されていなければ安定なシースが形成されない、というボーム条件（Bohm criterion）が知られています。シースはプラズマから壁へ流れ込む粒子とエネルギーの流れを支配し、荷電粒子の運動が壁付近でどう終端するかを決めます。プラズマ内部での粒子の運動そのものについては、荷電粒子の運動で扱います。

研究の最前線（博士課程レベル）

デバイ遮蔽とシースの物理は、古典的な理論でありながら、現在も核融合研究の最前線と直結しています。

最も活発な研究対象の一つが、プラズマ壁相互作用（plasma-wall interaction）です。トカマクのダイバータ（divertor）では、シースを通って壁に到達するイオンのエネルギーと角度が、壁材料のスパッタリング（sputtering）や損耗、熱負荷を左右します。シースを通過するイオンのエネルギー分布や、磁力線が壁に対して浅い角度で当たる磁化シース（magnetized sheath）の構造は、実験と運動論シミュレーションの両面から研究が進められています。ダイバータの工学的な設計と課題については、ダイバータを参照してください。

シースのモデリングでは、流体的な取り扱いだけでは捉えきれない現象が問題になります。壁前面での二次電子放出、非マクスウェル的な速度分布、過渡的な現象（ELM と呼ばれる周辺部の不安定性がもたらすバースト的な熱負荷）などでは、PIC（particle-in-cell）法による運動論シミュレーションが重要な役割を果たしています。ボーム条件を非マクスウェル分布や磁化された状況にどう一般化するかも、継続的に研究されているテーマです。

デバイ長は、プラズマ計測の原理にも深く関わります。ラングミュアプローブ（Langmuir probe）では、プローブ表面に形成されるシースの理論が測定精度を決めます。トムソン散乱（Thomson scattering）では、観測する波長スケール（波数の逆数）とデバイ長の大小関係が、散乱が個々の電子からの「非集団的散乱」になるか、電子の集団運動を反映した「集団的散乱」になるかを分けます。この境界を表す量がサモンパラメータ（Salpeter parameter） $\alpha = 1/(k\lambda_D)$ で、計測から密度や温度を導く際の基礎になります。

論文を読むときに頻出する英語キーワードとしては、Debye length、plasma parameter、quasineutrality、sheath、presheath、Bohm criterion、plasma-wall interaction、magnetized sheath などがあります。これらはいずれも、このページで積み上げてきた遮蔽・振動・準中性という基礎概念の延長線上にあります。

理解度チェック

問 1. プラズマ中に置いた電荷のまわりで、電気の影響がある距離までしか届かないのはなぜですか。

問 2. デバイ長の式を見て、温度が上がったときと密度が上がったときにデバイ長がどう変化しますか。

問 3. デバイ長とプラズマ周波数はどのような物理的関係で結ばれていますか。

問 4. 電離気体がプラズマとして振る舞うための 3 条件として正しい組み合わせはどれですか。

問 5. 準中性が破れる代表例であるシースとは何ですか。