アルファ粒子

核融合プラズマの文脈で「アルファ粒子」といえば、重水素–三重水素(D–T)核融合反応

D + T → ⁴He (3.5 MeV) + n (14.1 MeV)

で生まれる高エネルギーのヘリウム4原子核を指す。この 3.5 MeV の荷電粒子が磁場で閉じ込められ、背景プラズマ中で減速しながらエネルギーを渡すことで、外部加熱に頼らない「燃焼プラズマ」の自己加熱を担う。反応生成エネルギーの約 20% がアルファ粒子に与えられるため、プラズマの維持に必要な加熱電力をアルファ粒子だけで賄えるかどうかが磁気核融合の成否の核心であり、その生成・減速・軌道・損失・安定性への影響・計測はプラズマ物理学の一分野として独自に発展してきた。

核融合生成アルファ粒子と自己加熱

D–T 反応の生成物のうち中性子は磁場に閉じ込められずブランケットでエネルギー回収に回るが、アルファ粒子は荷電を持つため磁場に沿って運動し、コアプラズマ内に留まる。生成直後のエネルギーは 3.5 MeV で、背景の熱イオン(10–20 keV 程度)より二桁以上高く、このエネルギー差こそが自己加熱の原動力になる。アルファ粒子はまず電子との衝突でエネルギーを失い、電子が受け取った熱が一部はイオンに再分配される。したがって反応炉では主加熱源がアルファ粒子による電子加熱となり、現在の実験装置とは加熱の構図が根本的に異なる。

自己加熱は非線形な正帰還を含む。温度が上がると核融合反応率が増え、アルファ加熱が強まり、さらに温度が上がる。このフィードバックが輸送による熱損失と釣り合う点が燃焼点であり、温度分布・乱流輸送・アルファ加熱を自己無撞着に扱わないと出力予測ができない。JET の D–T 実験データを用いた統合シミュレーションでは、アルファ粒子を乱流計算の運動種として明示的に含めなくても、アルファ加熱を熱輸送方程式の熱源として落とすか否かで電子温度分布が約 1 keV 程度変わること、またこの条件ではアルファ粒子がイオン尺度乱流の輸送に直接与える影響は小さいことが示された。ただしこれはアルファ密度が比較的低い場合の話であり、ITER 級の密度までアルファ粒子濃度を上げると、外部加熱で作られる中速イオンで知られる乱流抑制や、逆にエネルギー粒子駆動不安定による閉じ込め劣化が現れる余地があり、抑制と不安定化のバランスはエネルギーも位相空間分布も大きく異なるアルファ粒子では単純に外挿できない。

アルファ粒子の生成はプラズマ中で一様ではない。反応率は温度に強く依存するため、生成率は中心部に尖った半径分布を持ち、断面内でも磁軸近傍に集中する。モンテカルロ法で核融合率からアルファ生成をサンプリングすると、生成率の半径分布は核融合率の分布とよく一致し、断面内では中心に偏った分布として可視化できる。

図は(a)核融合率とアルファ粒子生成率の半径分布、(b)断面内のアルファ粒子生成率の空間分布を示す。生成率が磁軸近傍に集中し、半径分布として核融合率とほぼ同じ形をとる様子が読み取れる。

減速とアルファ粒子分布

生成直後の 3.5 MeV のアルファ粒子は、背景の電子・イオンとのクーロン衝突でエネルギーを失いながら減速する。減速過程の物理は、アルファ粒子の速度が背景電子の熱速度より十分大きい間は電子への摩擦が支配的で、速度が下がって電子熱速度を下回るとイオンへのエネルギー移譲が支配的になる、というものである。両者の寄与が等しくなる速度に対応するエネルギーを臨界エネルギーと呼び、ここでアルファ粒子はエネルギーを電子に渡す役割からイオンを加熱する役割へと移る。減速はほぼ一定の減速率で進むため、定常状態での速度分布はマクスウェル分布ではなく、生成エネルギーから臨界エネルギー付近までほぼ一定の密度を持ち、それ以下で急峻に落ちる「減速分布」となる。この解析的な減速分布はシミュレーション結果とよく一致することが確かめられている。

図はアルファ粒子の速度分布関数について、シミュレーション結果、解析的な減速分布、熱平衡のマクスウェル分布を比較したものである。生成エネルギーと臨界エネルギーの位置が示され、高エネルギー側で平坦な減速分布がマクスウェル分布と大きく異なることが分かる。

速度空間での分布は等方的とは限らない。生成直後は反応核が等方的ならほぼ等方的だが、磁場に平行な運動と垂直な運動(バウンス運動)の違いに加え、プラズマ電流と同方向に進む粒子と逆方向に進む粒子ではドリフト軌道が通る磁場強度の平均が異なるため、定常分布には非対称性が残る。この同逆方向の非対称性はアルファ粒子の正味のトロイダル流を生み、新古典的な輸送係数の解析にも現れる。分布関数の近似式の構築や、非弾性衝突を含む輸送の評価など、燃焼プラズマ設計の入力として減速分布を正確に記述する試みが続けられている。

単一粒子軌道と閉じ込め

3.5 MeV のアルファ粒子のラーマー半径は熱イオンよりはるかに大きく、軌道は磁力線にほぼ沿った通過粒子と、磁鏡点で反射してバナナ型の軌道を描く捕獲粒子に分かれる。捕獲粒子の軌道幅(バナナ幅)は大きく、粒子はバウンス運動の間に半径方向に大きく移動する。この軌道幅が輸送の基本的なスケールを決める。軸対称なトカマクでは、衝突によるピッチ角散乱を含めても、減速期間中の軌道損失は純粋な D–T プラズマでは無視できる程度であることがモンテカルロ計算で示されている。一方、不純物が増えて有効イオン価数が上がると損失は増え、有効イオン価数が 4 を超えるような高不純物プラズマでは無視できなくなる。

軌道の内訳はプラズマの三次元性に敏感である。軸対称平衡では通過粒子と捕獲粒子の二種でほぼ分類できるが、プラズマ中心が螺旋状に変形した三次元平衡(ヘリカルコア)では、通過と捕獲の間の遷移軌道が現れ、エネルギー空間での分布や位相空間での局在が大きく変わる。アルファ粒子圧が高い場合、このような三次元平衡の形成自体がアルファ粒子の存在と結びついて議論される(後述)。

損失:リップル、誤差場、MHD摂動

トカマクのトロイダル磁場はコイルの離散性により周方向に波打つ(トロイダルリップル)。リップルが存在すると、バナナ軌道の粒子はリップルによるバナナドリフトを受けて半径方向に拡散し、特に減速過程で損失が生じる。リップル率 1% 程度の反応炉級トカマクでは、リップルによる追加のパワー損失はアルファ粒子の総出力の 10% 程度に達し、粒子数の損失割合はパワー損失の 1.5 から 1.8 倍に上る。損失過程はリップルによるバナナドリフトが支配的であり、非エルゴード的なバナナ粒子に対する衝突拡散係数の理論評価と軌道追跡計算はよく一致する。リップルが強い場合は軌道が確率的(ストカスティック)になり、減速中の損失がさらに増える。TFTR の D–T 実験ではリップル損失のモデル化が実際の実験条件で検証された。

もう一つの損失経路は、コイル配置の不完全さや制御コイルの残留磁場による誤差場である。平衡磁場に対する摂動磁場の比が 10⁻⁴ 程度以上になると、高エネルギーアルファ粒子の半径方向エネルギー輸送率が減速率に匹敵し得ることが、熱粒子の新古典トロイダル粘性理論を高エネルギー粒子に拡張した輸送理論から示されている。これは反応炉で許容できる誤差場の大きさに直接的な制約を課すものであり、ヘリカル装置を反応炉候補とする場合の設計指針にも関わる。

MHD 摂動も損失と再分布を引き起こす。ティアリングモードなどの磁気島はアルファ粒子の輸送と損失を増やし、非線形に発達した捕獲イオン不安定領域ではアルファ粒子の異常拡散が問題となる。また、m/n = 1/1 のキンク・準交換モードが飽和して形成されるヘリカルコアでは、アルファ粒子の空間分布が広がる。適度な広がりは加熱を均一化しエネルギー粒子勾配駆動の不安定を抑える方向に働くが、過度な広がりは加熱効率を下げ、軌道分類の上でも遷移軌道の割合が増えるなど閉じ込めの質が変化する。ITER 規模のハイブリッドシナリオではヘリカルコアの自発形成が予測されており、アルファ粒子圧が高い条件下での形成可否、形成への影響、形成後の閉じ込めという三点が数値研究の中心課題となっている。

アルフベン固有モード駆動とエネルギー粒子効果

アルファ粒子は背景プラズマのアルフベン速度に匹敵する速度を持つため、エネルギー粒子として MHD 安定性に強い影響を与える。代表例がアルフベン固有モード(AE)の駆動である。AE は磁場の幾何学配置によって許される離散的な固有モードで、連続スペクトルのギャップ内に存在するため減衰が弱く、わずかなエネルギー粒子駆動で励起されやすい。アルファ粒子の空間勾配が駆動源となり、共鳴条件(粒子の周回・バウンス・サイクロトロン運動とモードの位相速度の整合)を満たす粒子がエネルギーを波に渡す。

AE が励起されると、波と粒子の共鳴相互作用によってアルファ粒子が半径方向に再分布し、中心部の粒子密度が下がって加熱効率が低下する。イオンサイクロトロン共鳴加熱(ICRH)のような外部高周波波と同時に存在すると、サイクロトロン相互作用や位相空間不変量の変化によって粒子とモードの位相がデコリレートし、再分布がさらに強まる。複数の不安定モードが同時に存在する場合の再分布は単純に加算的ではなく、モード数を増やしても粒子再分布は飽和する。また、核融合反応による分布の補充と壁への損失が続く限り AE の駆動は持続する。

実験的にも、JET の D–³He 実験でアルファ粒子が AE を駆動した証拠が得られており、D–T 燃焼プラズマで AE が実際に問題になることを裏付けている。さらに、AE との相互作用は高周波加熱との組み合わせで現れる波の吸収の問題とも結びつく。アルファ粒子によるイオンサイクロトロン周波数帯の波の吸収や、低域混成波の準線形吸収は、電流駆動の効率に影響するため反応炉設計上の評価対象である。AE による再分布・損失を抑えつつアルファ加熱を維持できるかが、燃焼プラズマの運転設計における重要な制約となっている。

燃焼プラズマにおけるアルファ粒子の計測

アルファ粒子が自己加熱を担う以上、その密度分布・エネルギー分布・損失を直接測る診断が不可欠である。歴史的に重要なのは荷電交換再結合分光(CHERS)を応用した手法で、TFTR の D–T 実験では 0.15 から 0.6 MeV のエネルギー帯の非熱的アルファ粒子の半径分布が分光的に測定された。測定された高速アルファ密度は 10¹⁷ m⁻³ のオーダーで、信号は背景の制動放射の 1% 以下という微弱なものであったが、得られた分布は新古典拡散に小さな異常拡散(実効拡散係数で 0.1 m²/s 以下)を加えた減速分布と整合し、閉じ込めがほぼ古典的であることを示した。

空間・エネルギー分布を決める核物理的手法(アルファ粒子が誘起する核反応を利用するプローブ)も開発されてきた。次世代の診断として期待されるのが集合トムソン散乱で、ITER では核融合生成アルファ粒子の計測を目的とした装置の準備が進められている。これは高エネルギー粒子の速度分布を散乱スペクトルから直接求められる手法であり、燃焼プラズマの定常計測の主力になると見込まれる。また D–T 以外の反応系でも、Zピンチ放電での陽子–ホウ素反応によるアルファ粒子生成の初期測定が行われるなど、アルファ粒子計測は反応系ごとに広がっている。さらに、ITER 級のディスラプションではアルファ粒子がランナウェイ電子の動力学に影響を与え得ることも指摘されており、計測と制御の対象は定常運転時にとどまらない。

このようにアルファ粒子研究は、単一粒子軌道論から始まって減速分布の記述、リップル・誤差場・MHD 摂動による損失の評価、AE 駆動による再分布の理解、そして燃焼プラズマでの直接計測へと、理論・シミュレーション・実験が密接に連動して発展してきた。ITER での D–T 運転が近づく中、アルファ加熱と乱流・MHD 安定性の自己無撞着な理解が、自己点火プラズマの実現を左右する物理の中心に位置し続けている。