[論文レビュー] Monte Carlo simulations of the electron-gas interactions in the KATRIN experiment

J. Kellerer, F. Spanier|arXiv (Cornell University)|Dec 31, 2021

Atomic and Molecular Physics参考文献 28被引用数 3

ひとこと要約

本論文では、KATRIN実験のウィンドウレスガス状トリチウムソース（WGTS）における電子-ガス相互作用をモデル化する新しいモンテカルロシミュレーションコードKARLを提示する。このコードは、エネルギースペクトル、粒子密度、電流分布を計算するために、電子およびイオンの軌道を追跡する。主な貢献は、自己整合的かつ空間的に分解能のある電子およびイオンダイナミクスのシミュレーションであり、特にプラズマポテンシャルモデル化とニュートリノ質量感度に不可欠な、顕著な半径方向のイオン電流および位置依存スペクトル特徴を明らかにしている。

ABSTRACT

At the KATRIN experiment, the electron antineutrino mass is inferred from the shape of the $\beta$-decay spectrum of tritium. Important systematic effects in the Windowless Gaseous Tritium Source (WGTS) of the experiment include the energy loss by electron scattering, and the extended starting potential. In the WGTS, primary high-energy electrons from $\beta$-decay produce an extended secondary spectrum of electrons through various atomic and molecular processes including ionization, recombination, cluster formation and scattering. In addition to providing data essential to the simulation of energy loss processes, the electron spectrum also provides information important in the simulation of plasma processes. These simulations will then provide an insight on the starting potential. Here, a Monte Carlo approach is used to model the electron spectrum in the source for a given magnetic and electric field configuration. The spectrum is evaluated at different positions within the WGTS, which allows for a direct analysis of the spectrum close to the rear wall and detector end of the experiment. Alongside electrons, also ions are tracked by the simulation, resulting in a full description of the currents in the source.

研究の動機と目的

• KATRINのWGTSにおける荷電粒子ダイナミクスの自己整合的モンテカルロシミュレーションを構築すること。ここには、電子およびイオンを含む。
• ソース内における空間的に依存する電子エネルギースペクトル、粒子密度、および電流分布を解像すること。
• プラズマポテンシャル効果のモデル化に不可欠な、特に半径方向のイオン電流を正確に推定すること。
• 完全な粒子追跡、二次プロセス、および場の結合を含めることで、先行研究のシミュレーションを改善すること。
• 全電流および半径方向電流の実験的測定値と直接比較できるようにし、プラズマポテンシャル再構築を支援すること。

提案手法

• 2.5 Tの磁場と電場配置下で、個々の電子およびイオンの軌道をシミュレートするモンテカルロスキームを実装する。
• イオン化、再結合、弾性散乱のための更新済み断面積を用いて、電子とトリチウム分子との相互作用をモデル化する。
• 二次粒子を追跡し、自己整合的密度場を通じて局所的プラズマ密度へのフィードバックを組み込む。
• 複数の半径方向および軸方向位置に仮想障壁を設け、粒子電流密度およびエネルギースペクトルを計算する。
• ローレンツ力下での運動方程式を用いて粒子運動を統合し、弾性散乱に起因するジグロ運動効果を含める。
• 既知の物理的限界および先行結果（特に電子エネルギー損失およびスペクトル形状）と照合して、シミュレーションを検証する。

実験結果

リサーチクエスチョン

• RQ1 WGTS内での電子エネルギースペクトルは、特に後方壁およびビームチューブ近傍でどのように空間的に変化するか？
• RQ2 電子およびイオンの半径方向および軸方向電流密度は何か？また、磁場強度およびソース条件にどのように依存するか？
• RQ3 二次電子プロセス（イオン化、再結合、散乱）は、ソース内での一次電子スペクトルにどの程度の影響を及ぼすか？
• RQ4 粒子密度および電流分布は、WGTS内の入口ガス密度および温度にどのように依存するか？
• RQ5 シミュレーションは、測定可能な電流（例：後方壁電流およびDPS電流）を再現できるか？また、測定不能な成分（例：半径方向イオンフラックス）を予測できるか？

主な発見

• シミュレートされた電子スペクトルは、三つの明確な領域（β崩壊領域、熱的領域、二次電子領域）を示し、その相対的寄与はソース内の位置によって変動する。
• ビームチューブ壁への半径方向のイオン電流は無視できない大きさであり、2.5 Tの磁場下で後方壁から5 cmの位置で最大8 × 10^6 cm⁻²s⁻¹に達する。
• イオン密度プロファイルは強い半径依存性を示し、ビームチューブ壁に近づくにつれて増加し、2.5 Tおよび6.3 × 10¹⁴ cm⁻³の入口密度条件下でピーク密度が10^6 cm⁻³を超える。
• 半径電流および密度分布は磁場強度に敏感であり、1 Tの低磁場では、グリロ半径の増大および散乱誘発の半径拡산により電流が減少する。
• シミュレーションは、後方壁およびビームチューブ近傍で、電子エネルギー損失および二次電子生成がスペクトルを顕著に変化させることを明らかにした。特に低エネルギー領域で顕著である。
• コードは既知のスペクトル的特徴を正確に再現し、完全なイオン追跡および空間的に分解能のある診断を含めることで、先行研究のシミュレーションを拡張した。これにより、将来的な半径方向電流測定と直接比較可能となった。