Systematic comparison of VMEC and HINT equilibrium calculations for finite-beta LHD plasmas

本論文は、Large Helical DeviceプラズマにおけるVMECとHINTの平衡計算を体系的に比較しており、両コードは低ベータでは一致するものの、高ベータ値においては、VMECの入れ子状の磁気面仮定では表現できないエッジのストカスティシティ（確率論的性質）や磁気面の崩壊をHINTが捉えるため、両者の間に乖離が生じることを明らかにしている。

要約

非常に奇妙でねじれたオーブンの中で、完璧な丸いケーキを焼こうとしている場面を想像してみてください。核融合エネルギーの世界では、科学者たちはLHD（大型ヘリカル装置）のような「ステラレータ」と呼ばれる機械を使って、超高温のプラズマを閉じ込めています。このプラズマを安定させるためには、それを保持する磁気的な「壁」がどのような形になるべきかを正確に計算する必要があります。

この論文は、これら2つの異なる「パン屋」（コンピュータプログラム）が、プラズマが非常に高温かつ高圧になったときに、磁気的な壁の形状をどのように解明しようとするかを比較したものです。

2人のパン屋：VMECとHINT

1. VMEC（厳格な建築家）： このプログラムは、ケーキの各層が完璧で滑らかな、入れ子状の玉ねぎでなければならないと主張する建築家のようです。彼は、磁気的な壁が決して壊れたり、互いに触れ合ったりすることはないと仮定しています。単純で低圧な状況においては非常に優れていますが、弱点があります。それは、壁が乱れたり壊れたりすることを決して認めないという点です。
2. HINT（現実的な観察者）： このプログラムは、実際にケーキが焼かれる様子を見守る科学者のようなものです。彼は、層が完璧であるという前提を持ちません。代わりに、物理現象が自然に起こるままにさせておきます。もし熱が上がりすぎれば、磁気的な壁がぐらついたり、バラバラになったり、あるいは混沌とした混乱状態に陥ったりすることを許容します。

実験：熱を上げる

研究者たちは、LHD装置を用いて、3種類の異なる形状の磁気的な「オーブン」（内側にシフトしたもの、外側にシフトしたものなど）を用いて、これら2つのプログラムをテストしました。彼らはプラズマの圧力（ケーキの「熱」）を0%から5%へとゆっくりと上昇させました。

低圧のときはどうなったか？
プラズマが冷たく穏やかだったとき、2人のパン屋の意見は一致しました。磁気的な壁は滑らかで、入れ子状のまま維持されていました。まさに「厳格な建築家（VMEC）」が予測した通りでした。すべては順調でした。

熱が上がるとどうなったか？
圧力が特定の「臨界点」を超えると、2人のパン屋の間で意見の相違が生じ始めました。

• VMECは、完璧で滑らかな、膨張していく玉ねぎの層を描き続けました。彼は、プラズマがただ大きくなり、丸くなっていると考えていました。
• HINTは、異なる光景を見ていました。彼は、磁気的な壁が「ストカスティック（確率論的／カオス的）」になり始めていることに気づいたのです。

「ストカスティック」な混乱：創造的な比喩

磁力線をスパゲッティの束だと考えてみてください。

• 完璧な状態（低圧）では、スパゲッティの麺はきれいに束ねられ、互いに平行に走っています。
• 圧力が上昇すると、プリフシュル・シュリッター電流（プラズマの中に自然に形成される一種の電流）が、まるでスパゲッティを混ぜ合わせる混沌とした手のように作用します。
• 最終的に、麺同士が重なり合い、絡まり合います。これは磁気島や**ストカスティシティ（カオス性）**と呼ばれます。整然とした「玉ねぎの層」は崩壊します。

HINTはこの絡まり合いを許容するため、磁気の「ケージ（籠）」が縮小していることを検知します。エッジ部分での混沌とした混合により、実効的な体積が小さくなるのです。しかし、VMECは依然として完璧に膨張する玉ねぎを描いているため、体積が増加していると考えてしまいます。

主な知見

1. 「転換点」： 整然とした玉ねぎの層が壊れる、特定の圧力レベルが存在します。この点を超えると、VMECは壊れた壁を見ることができないため、もはや正確ではなくなります。
2. 形状の影響： 「転換点」は、機械が外側にシフトしている場合、より低い圧力で発生します。
   • 比喩： 外側にシフトした機械が、ぐらつくテーブルだと想像してください。それは、頑丈な内側シフトのテーブルよりも、簡単にひっくり返ってしまう（混沌を生み出す）ものです。外側の形状は磁場により多くの「さざ波」を生み出し、スパゲッティをより早く絡ませます。
3. 体積の減少： 外側にシフトした構成や標準的な構成において、圧力が非常に高くなると、実際のプラズマの体積（現実的なHMTモデルによる）は、磁気的な壁が崩壊するために減少し始めます。VMECはこの現象を完全に見逃しており、体積が増え続けていると考えています。

結論

この論文は、高圧の核融合プラズマにおいては、「完璧な玉ねぎ」のモデル（VMEC）だけに頼ることはできないことを示しています。磁気的な壁が崩壊し、カオスへと変化する様子を見るためには、「現実的な観察者（HINT）」が必要です。これは、特に外側にシフトした構成において顕著であり、その場合、磁場はこれらの乱雑な3次元的効果に対してより敏感になります。この研究は、より高いエネルギーを目指すにつれて、「滑らかで完璧な磁気層」という仮定が、ますます妥当性を失っていくことを裏付けています。

技術概要

技術要約：有限β LHDプラズマにおけるVMECとHINTの平衡計算の系統的な比較

問題提起
ステラレータ物理学において、三次元磁気流体力学（MHD）平衡の計算は、頻繁にVMECソルバーを用いて行われる。VMECは変分原理に基づいており、完全にネストされた磁束面が存在することを厳密に仮定している。しかし、有限βプラズマにおいては、プラズマ圧の上昇が磁場トポロジーを変化させ、磁気アイランドやストカスティック（確率的）領域を生成することがある。これらの現象は、特にネット電流ゼロの装置において、Pfirsch-Schlüter（PS）電流によって引き起こされ、ネストされた磁束面の仮定の妥当性を損なう可能性がある。磁気アイランドが重なり合うと、最後の閉じた磁束面（LCFS）が崩壊し、囲まれたプラズマ体積が減少する。この制限は、シャラノフ・シフトよりも深刻であると観察されている。一方、HINTコードは、ネストされた磁束面を仮定せずに磁場と圧力の動的方程式を解くものであるが、LHD（Large Helical Device）における様々な構成およびベータ値にわたるVMECとHINTの系統的な比較が必要であった。これは、これらの仮定の影響を明確にするためである。

手法
本研究では、LHDプラズマにおけるVMECとHINTの平衡計算の系統的な比較を行った。調査は、3つの特定の真空磁軸構成、$R_{ax}^V = 3.53$ m（内側シフト）、$3.60$m（標準）、および$3.85$m（外側シフト）に焦点を当てて実施された。解析は、オン軸ベータ値（$\beta_0$）を$0.0%$から$5.0%$まで$0.5%$刻みでカバーした。

両コードとも、初期圧力プロファイルを $p = p_0(1-s)$ （ここで $s$ は正規化されたトロイダルフラックス）として定義した。本研究では、主に以下の3つの指標を比較した：

1. 磁軸の位置（$\langle R_{ax} \rangle$）。
2. 軸上の回転変成率（$\iota/2\pi$）。
3. 明確なLCFSによって囲まれた体積。

磁気アイランドやストカスティック領域を解明するために緩和法を用いるHINTコードを、ネストされた磁束面の制約を課すVMECと比較した。

主な結果

• 低β領域： 低い$\beta_0$値において、VMECとHINTは一致した平衡解を生成した。これは、この領域ではネストされた磁束面の仮定が概ね有効であり、両コード間の3D MHD応答が同等であることを示している。
• 臨界$\beta_0$と乖離： 構成に依存する臨界$\beta_0$が特定された。この閾値を超えると、両コードの解は乖離し始めた。HINTの計算では、磁軸はVMECと比較して、より大きなシャラノフ・シフトと高い回転変成率を示した。
• 磁束面の崩壊と体積： HINTでは、3つの構成すべてにおいて、高いベータ値で明確なLCFSによって囲まれた体積が減少した。この減少は、プラズマ端部付近でのストカスティック磁場の進化による磁束面の崩壊に起因する。対照的に、ネストされた磁束面の仮定に制約されているVMECは、$\beta_0$の増加に伴い囲まれた体積が単調に増加しており、磁束面の崩壊を表現できていなかった。
• 構成依存性： 臨界$\beta_0$の値は、内側シフトから外側シフトの構成に向かって減少した。外側シフト構成（$R_{ax}^V = 3.85$ m）が、3D平衡応答に対して最も敏感であることが判明した。この感度は、$R_{ax}^V$が増加するとヘリカル・リップルが増大して強いPS電流を駆動する一方で、磁気シアが減少することに起因する。強いPS電流と低いシアの組み合わせにより、より大きな磁気アイランドが生じ、より低いベータ値で早期にストカスティシティが発生する。

意義と主張
本論文は、LHDプラズマにおける3D平衡応答が、構成が内側から外側へシフトするにつれてますます重要になることを確立している。本研究は、VMECの中心であるネストされた磁束面の仮定が、有限ベータにおいて、PS電流駆動の摂動および結果としての端部ストカスティシティによって損なわれることを示している。

著者らは、VMECとHINTは低$\beta_0$構成においては同等であるが、HINTは磁束面の崩壊とそれに伴うプラズマ体積の減少を解明することで、高ベータにおいてより物理的に完全な記述を提供すると主張している。両コードの乖離は、ネストされた磁束面の仮定の限界が一様ではなく、特にトロイディティ（PS電流を駆動する）と磁気シアの相互作用という、真空磁気構成に強く依存していることを浮き彫りにしている。これらの結果は、特に臨界ベータ閾値が低い外側シフト構成において、高ベータ平衡を解析する際に、HINTのようなストカスティシティを解像できるコードを使用することの重要性を強調している。