FLASH放射線治療を支える新モデル：電離箱における電荷収集効率の謎を解明

10月 29, 2026 · 論文解説

🎧 この記事を音声で聴く（AI生成Podcast）

解説対象論文: Charge collection efficiency in air-vented plane-parallel ionisation chambers at ultra-high dose rates: a self-consistent Garfield++ Monte Carlo model including space-charge effects and ion recombination
Physics in Medicine and Biology｜被引用数: 0（2026年10月28日時点）

FLASH放射線治療は、短時間での超高線量率照射により正常組織へのダメージを抑える画期的な治療法ですが、その正確な線量測定には大きな課題があります。特に、平行平板型電離箱（PPIC）では、イオン再結合や空間電荷による電場歪みが電荷収集効率（CCE）を低下させます。この研究では、Monte Carloシミュレーションツールキット「Garfield++」を拡張し、これらの複雑な現象を自己無撞着にモデル化するフレームワークを開発しました。この新しいモデルは、既存の分析モデルや実験データと高い精度で一致し、CCE低下の主要因が自由電子割合（FEF）の変化であることを明らかにしました。これは、FLASH線量測定の精度向上とリアルタイム補正法の開発に向けた重要な一歩となります。

研究サマリー（FINER & PICO）

Feasible実施可能性	既存のオープンソースフレームワークGarfield++を拡張し、文献との比較で妥当性を確認しており、実現可能性が高い。
Ethical倫理面	本研究は既存技術の向上と新しい医療技術の安全かつ正確な導入に貢献し、倫理的な課題は提起されていない。
Interesting面白さ	FLASH放射線治療の主要な課題である線量測定精度向上に直結し、医療画像インフォマティクス分野で高い関心を集める。
Novel新規性	空間電荷効果とイオン再結合を自己無撞着に結合した3次元モンテカルロモデルは、従来の1次元モデルを超える新規性を持つ。
Relevant切実さ	FLASH放射線治療の臨床応用を促進するために不可欠な、高精度な線量測定モデルの開発に直接的に貢献する。
Measurable測定可能性	電荷収集効率、自由電子割合、電場分布、誘導電流といった物理量を定量的にシミュレーションし、既存モデルや実験データと比較可能である。
Modifiable派生・改善	開発されたGarfield++フレームワークは、ガス混合物、電離箱の形状、照射条件などを柔軟に変更してシミュレーション可能である。
Structured文章構造	理論的背景、モデル開発、検証（既存モデル・数値データとの比較）、結果分析、不確実性検討、結論・展望と明確な構成を持つ。
PICO / PECO対象・介入・比較・結果	P(対象): 平行平板型電離箱(PPIC)。 I(介入): 超高線量率(UHDR)照射。 C(比較): 既存の分析・数値モデル。 O(結果): 電荷収集効率(CCE)、自由電子割合(FEF)、空間電荷による電場歪み。

はじめに：最先端FLASH放射線治療と線量測定の課題

どうも、Beyond the Pixelです。近年、放射線治療の分野で注目を集めている「FLASH放射線治療」をご存知でしょうか？これは、通常の放射線治療よりもはるかに高い線量率（UHDR: Ultra-High Dose Rate）で短時間に放射線を照射することで、腫瘍には同等以上の効果を維持しつつ、正常組織への副作用を劇的に低減できるという画期的な治療法です。2014年の発見以来、FLASH効果に関する研究は急速に進展し、様々な前臨床モデルや動物実験でその有効性が確認されています。例えば、マウスの腹部照射後の生存率向上、全脳照射後の認知機能保持、脳毒性低減などが報告されており、放射線治療の治療域を広げる可能性を秘めています。

しかし、この有望なFLASH放射線治療の実用化には、乗り越えるべき重要な課題があります。それは、「正確な線量測定」です。国際原子力機関（IAEA）が推奨する基準線量計である電離箱、特に平行平板型電離箱（PPIC: Plane-Parallel Ionisation Chamber）は、放射線治療において基準線量測定や治療中のビームモニタリングに広く使われています。ところが、超高線量率であるUHDR条件下では、電離箱の「電荷収集効率」（CCE: Charge Collection Efficiency）が大きく低下することが実験的に示されており、正確な線量測定を困難にしています。この効率低下の主な原因は、イオン再結合の増加と、高い電荷密度による「空間電荷」によって引き起こされる電場歪みにあります。これらの複雑な現象を正確にモデル化し、予測するためのフレームワークが求められていました。

既存モデルの限界とGarfield++の登場

従来の放射線治療では、イオン再結合による電荷収集の不完全さを補正するために、Boagモデルに基づいた補正係数が用いられてきました。しかし、FLASH放射線治療のようなUHDR条件下では、これらの分析モデルはもはや有効ではありません。特に、空間電荷によって引き起こされる電場歪みは無視できないレベルに達し、局所的に電場が数倍に増加したり、ほぼゼロにまで減少したりすることが複数の研究で示されています。

この課題に対処するため、この研究では、Monte Carlo（モンテカルロ）粒子輸送フレームワークである「Garfield++」を拡張しました。この拡張により、イオンとイオンの再結合プロセスと、空間電荷による電場の自己無撞着（self-consistent）な計算を実装し、粒子輸送、電子付着、再結合、そして動的な電場歪みを同時に考慮できるモデルを開発しました。Garfield++の柔軟性により、任意のガス混合物、大気条件、電離箱の形状、入射粒子の種類、ビームの空間的および時間的構造など、広範な照射条件下での再結合現象の包括的な調査が可能になります。

**Figure 3.** Figure 3: Overview of the main classes in Garfield++ and their interplay (from [37]).解説: Garfield++の主要なクラスとその相互作用の概要図です。Mediumクラスで媒体特性を定義し、ComponentとSensorクラスで検出器の形状と電場を構築、TrackやDriftクラスで粒子の輸送と電離をシミュレーションする流れが示されています。

Garfield++シミュレーションの典型的な構造は、まずMediumクラスで検出器の媒体（ガス混合物など）を定義します。次に、Componentクラスで検出器の形状と電場・磁場を設定し、Sensorクラスでこれらのコンポーネントを組み合わせて活性検出体積を定義し、信号計算を管理します。その後、TrackHeedクラスで一次電離をシミュレーションし、AvalancheMicroscopicやAvalancheMCといった専用の輸送クラスを用いて電子やイオンの運動をシミュレーションします。

電荷収集のメカニズム：電子付着とイオン再結合

空気中の電離箱では、入射した放射線によって電子と陽イオンのペアが生成されます。自由電子は、酸素分子などの電気陰性分子に付着して陰イオンを形成します。このプロセスは、二体付着と三体付着の二つのメカニズムで起こります。

二体付着は、酸素分子が電子を捕獲して解離する過程で起こり、最小限の解離エネルギーが必要です。水分子も同様に陰イオンを形成しますが、より高い解離エネルギーが必要です。

一方、低エネルギーの電子では、三体付着が支配的になります。このメカニズムでは、通常O2やH2Oである第三の分子が必要となり、準安定な酸素陰イオンを脱励起させることで安定な陰イオンが形成されます。

**Figure 1.** Figure 1: Electron collision energy at 20 Td in humid air with 2 % H2O.解説: 湿潤空気中、20 Tdの電場における電子衝突エネルギーの分布を示すグラフです。平均衝突エネルギーは1.135 eVであり、低エネルギーでの衝突が支配的であることがわかります。

図1は、湿潤空気中における電子の衝突エネルギー分布を示しており、この条件下では平均電子衝突エネルギーが1.14 eVと低く、三体付着プロセスが主に起こることが示唆されています。

イオン再結合には、「初期再結合」と「体積再結合」の二種類があります。初期再結合は、同じ電離トラック内で生成された反対の電荷を持つイオン間で起こり、線量率には依存しません。一方、体積再結合は異なる粒子トラックに由来するイオン間で起こり、高線量率条件下で支配的になります。この研究では、UHDR条件に焦点を当てているため、体積再結合のみを考慮しました。体積再結合は、反対電荷のイオン密度の積に比例する速度で起こります。

**Figure 2.** Figure 2: Comparison of the electron attachment coeﬀicient and electron mobility in dry and humid air obtained from Magboltz simulations with experimental measurements reported in the literature. Experimental data correspond to dry air and humid air with 2 % H2O [21, 25], while data from Hochhäuser et al. [26] correspond to humid air with 1.2 % H2O.解説: 乾燥空気および湿潤空気中での電子モビリティと電子付着係数について、Magboltzシミュレーションの結果と文献で報告されている実験データを比較したグラフです。特に低電場域を除き、両者で良好な一致が確認できます。

電子モビリティ（動きやすさ）と電子付着係数については、Magboltzを用いたシミュレーション結果が文献の実験値とよく一致することが確認されています（図2）。特に、電離箱の一般的な動作電場範囲外でのわずかな差異を除けば、良好な合致が得られました。しかし、イオンのモビリティやイオン-イオン再結合係数については、文献値に大きなばらつきがあり、これらがシミュレーションにおける不確実性の主要な原因の一つとなっています。

Garfield++におけるイオン再結合と空間電荷の実装

この研究では、Garfield++のAvalancheMCクラスにイオン-イオン再結合プロセスを実装しました。これは電子付着と同様のアプローチで、各モンテカルロ時間ステップで再結合確率を評価し、確率的にイオンをシミュレーションから除去するものです。陽イオンと陰イオンの密度を評価するために、ComponentGridクラスを拡張し、パーティクル・イン・セル（PIC）法を用いて粒子密度マップを計算できるようにしました。

**Figure 4.** Figure 4: Schematic illustration of the Particle-in-Cell (PIC) deposition scheme in a two-dimensional grid.解説: 2次元グリッドにおけるパーティクル・イン・セル（PIC）法による粒子配置スキームの概念図です。グリッドセル内の灰色粒子が、その相対位置に応じて周囲の4つのノードに重みを分散させる様子が示されています。

図4は、PIC法における粒子配置スキームの概念図です。シミュレートされた各粒子は、周囲のグリッドノードにその重みを分散させ、粒子密度を計算します。これにより、離散的なデータが平滑化され、グリッド点間の電場がより正確に、空間的に連続的に推定されます。

高電荷密度によって生じる空間電荷効果を考慮するため、ComponentGridクラス内にPoisson（ポアソン）方程式ソルバーを実装しました。このソルバーは、空間電荷密度マップから静電ポテンシャルを計算し、それに基づいて電場を更新します。電極の電位を固定するディリクレ境界条件と、電気的な流れがないことを示すノイマン境界条件を適用してPoisson方程式を解きました。

**Figure 5.** Figure 5: Diagram of the main components of the sim- ulation code. Multithreaded sequences are shown in green, while single-threaded sequences are shown in blue and red.解説: シミュレーションコードの主要なシーケンスを示す図です。粒子ペアの生成、空間電荷と密度マップの計算、電子と陽・陰イオンのドリフトという一連のプロセスが示されており、緑色の部分はマルチスレッド処理、青色と赤色の部分はシングルスレッド処理が施されています。

シミュレーションコードは、図5に示されるようなフローで実行されます。まず電子とイオンのペアを生成し、空間電荷と密度マップを計算します。その後、電子、陰イオン、陽イオンのドリフトをシミュレーションし、このプロセスを繰り返します。計算コストを削減するため、各シミュレートされた電子またはイオンが複数の物理粒子を表す「粒子重み」を使用しました。また、並列処理とグローバルな時間ステップの導入により、共有データ構造へのアクセスを調整し、効率的なシミュレーションを実現しています。

検証結果：電荷収集効率と空間電荷のダイナミクス

開発されたモンテカルロモデルは、既存の分析モデルや数値モデルと比較することで、その正確性が検証されました。

1. 空間電荷がない場合の電荷収集効率（CCE）

まず、空間電荷効果を無効にした状態で、FenwickとKumarが2023年に提案したモデルと比較しました。このモデルは、負イオンの空間分布に関するBoagモデルの単純化された仮定を改善したものです。電極ギャップサイズが1.0 mmと2.0 mmのPPICについて、本研究のモンテカルロシミュレーションはFenwickモデルと非常に良い一致を示しました。

**Figure 6.** Figure 6: Charge collection eﬀiciency as a function of ionisation density for an instantaneous pulse in PPICs with 1 mm and 2 mm gaps operated at 500 V. Garfield++ results with space-charge effects disabled (markers) are compared with the Fenwick and Boag models (lines).解説: 電極ギャップ1 mmと2 mmのPPICにおいて、瞬間的なパルス照射時の電離密度に対する電荷収集効率（CCE）を示したグラフです。Garfield++シミュレーション結果は、FenwickモデルおよびBoagモデルと、特に高電離密度でのCCE低下傾向において比較されています。

図6に示すように、電離密度が増加するにつれてCCEが低下する傾向が、両ギャップサイズで正確に再現されました。これは、本モデルが再結合効果を正確に記述できることを示しています。

2. 空間電荷の影響

次に、空間電荷効果を有効にした場合の電場を評価しました。均一に帯電した円筒状分布の解析解と比較することで、Poissonソルバーの精度が確認されました。

**Figure 7.** Figure 7: Comparison between the electric field com- puted using the Poisson solver (red) and the corre- sponding analytical solution (black) for a uniformly charged cylindrical distribution. The dashed lines are the limit of the cylinder.解説: 一様に帯電した円筒状分布に対して、Poissonソルバーによって計算された電場（赤線）と対応する解析解（黒線）を比較したグラフです。大部分の領域で両者が非常によく一致していることが示されています。

図7に示す通り、数値解による電場は解析解と大部分の領域でよく一致しました。

空間電荷が電場に与える影響を調査するため、2.5 µsのパルスで2 Gyを照射した際の電場（バイアス電圧300 V、ギャップ1 mm）の時間的変化をシミュレーションしました。

**Figure 8.** Figure 8: Temporal evolution of the electric field magnitude along the beam isocentre in a 1 mm gap PPIC operated at 300 V, for a pulse delivering 2 Gy over 2.5 µs.解説: 1 mmギャップPPIC（300 V動作）における、2.5 µsで2 Gyを照射したパルス後のビーム中心線に沿った電場強度の時間的変化を示したグラフです。空間電荷により電場が大きく歪められ、局所的に電場が数倍に増加したり、ほぼゼロにまで低下したりする様子が描かれています。

図8aは、電場が照射パルス中だけでなく、その後の電荷ドリフト段階でも空間電荷によって大きく歪められることを示しています。電場は局所的に名目電場の4倍を超える場合がある一方、電離箱の他の領域ではほぼゼロにまで低下することもあります。これはKranzerらによる結果（図8b）とも非常によく一致しています。

3. 自由電子割合（FEF）の挙動

さらに、UHDR条件下での自由電子割合（FEF: Free Electron Fraction）の挙動を調査しました。Boagの理論ではFEFは一定とみなされますが、空間電荷による強い電場歪みは、高パルス線量（DPP: Dose Per Pulse）でこの仮定を無効にする可能性があります。

**Figure 9.** Figure 9: Free electron fraction as a function of DPP for different gap sizes of PPICs at 300 V and a pulse duration of 2.5 µs. Results from the present simulation are compared with those reported by Paz-Martín et al. [48] and the constant analytical model.解説: 電極ギャップが異なるPPIC（300 V、2.5 µsパルス）において、パルス線量（DPP）に対する自由電子割合（FEF）の挙動を示したグラフです。DPPが増加するにつれてFEFが低下する傾向が、Paz-Martínらの結果および分析モデルと比較されています。

図9は、DPPに対するFEFを示しています。低DPPではFEFは一定で分析予測と一致しますが、DPPが増加するにつれてFEFは徐々に低下します。これは、空間電荷によって電場が低下する領域が広がるため、電場に強く依存する電子付着係数が増加し、電極に到達する自由電子の割合が減少するためと考えられます。Paz-Martínらの報告とも良好な一致が見られました。

4. 電荷収集効率（CCE）と誘導電流

空間電荷効果を考慮したCCEについても調査を行いました。

**Figure 10.** Figure 10: Charge collection eﬀiciency as a function of DPP for a PPIC with a 2 mm electrode gap at four applied voltages and a pulse duration of 2.5 µs. Results from the present simulation are compared with those reported by Paz-Martín et al. [48].解説: 2 mm電極ギャップのPPICにおいて、4つの異なる印加電圧（100 V, 200 V, 300 V, 400 V）でのパルス線量（DPP）に対する電荷収集効率（CCE）を示したグラフです。本シミュレーションの結果は、Paz-Martínらの報告と非常に良い一致を示しています。

図10は、異なるバイアス電圧におけるDPPに対するCCEを示しており、Paz-Martínらの結果と非常によく一致しています。これにより、本モデルがUHDR照射条件下での強い空間電荷による電場歪みの存在下でも再結合効果を正確に再現できることが実証されました。

最後に、電子とイオンのドリフトによって生成される誘導電流を計算しました。

**Figure 11.** Figure 11: Induced current generated by the drift of electrons and ions. Results from the present simulation compared with those reported by Paz-Martín et al. [48]. For a 2.5 µs pulse in a PPIC with a 2 mm electrode gap and a radius of 7.8 mm.解説: 電子とイオンのドリフトによって生成される誘導電流のグラフです。短い電子信号の後に、長いイオン誘導電流のテールが続いており、本シミュレーションの結果はPaz-Martínらの報告とよく比較されています。特にイオン誘導電流で良好な一致が見られます。

図11は、短い電子信号に続いて、より長いイオン誘導電流のテールを示しています。パルス開始時の立ち上がり部分でわずかな差異が見られるものの、イオン誘導電流はPaz-Martínらの結果と非常に良い一致を示しました。このわずかな差異は、本研究で仮定した矩形波状のパルスプロファイルと、参照研究で用いられた実験的に測定されたパルス形状の違いに起因すると考えられます。この結果はまた、DPPが電荷収集時間に有意な影響を与えないというPaz-Martínらの驚くべき観察を裏付けています。

これは、UHDR条件下でのCCE低下が、正負イオンの共存時間が長くなることによるものではなく、空間電荷による電場歪みが電子付着係数を変化させ、FEFを減少させることによって主に引き起こされることを示唆しています。FEFの減少は負イオン密度を増加させ、結果としてイオン-イオン再結合確率を高めます。

研究の限界と今後の展望

このシミュレーションにおける主要な物理的不確実性は、イオンモビリティとイオン-イオン再結合係数に起因します。これらのパラメータは、電離箱内でのイオン輸送と再結合率を直接支配しますが、文献値に大きなばらつきがあるためです。また、モンテカルロシミュレーションの数値的実装にも不確実性があります。計算コストを抑えるために、各シミュレートされた粒子が多数の物理イオンまたは電子を代表する「重み付き粒子」が用いられており、これにより統計的な変動が生じます。

さらに、粒子や電荷密度が離散的なボクセルグリッド上に表現されるため、結果は選択された空間分解能に依存します。より細かいボクセル化は電荷分布の空間的記述を改善しますが、各ボクセルに含まれる粒子の数が減少するため、統計的変動が増加します。これらの不確実性をすべて評価するには膨大な数のシミュレーションが必要となるため、本研究の範囲外とされています。

しかし、本研究で開発されたGarfield++フレームワークは、イオン-イオン再結合と空間電荷による電場歪みという、UHDR条件下でのPPICシミュレーションに不可欠な二つの主要な物理メカニズムを実装しました。これにより、UHDR条件下での電離箱における再結合現象を正確に再現できる、完全で柔軟かつ自己無撞着なモンテカルロフレームワークが提供されました。特に、横方向の均一性を仮定せず、3次元の粒子ベースのモンテカルロアプローチで電荷輸送を追跡できる点は、従来の1次元モデルとは異なる大きな利点です。これにより、ビームの空間的・時間的構造の影響をより柔軟に調査できます。

本研究の最も重要な発見は、UHDR条件下でのCCE低下が、主に自由電子割合（FEF）の変化によって支配されているという示唆です。空間電荷による電場歪みが電離箱内の広い領域で電場強度を低下させ、その結果、電子付着確率が増加し、FEFが減少します。この挙動が将来の研究によってさらに確認されれば、既存の再結合理論を改善する新しい分析モデルや補正方法の開発につながり、UHDR照射条件下でのPPIC線量測定におけるリアルタイム補正の可能性が開かれるでしょう。これは、FLASH放射線治療の線量測定モデルを改善し、その臨床応用を加速するための強固な基盤となります。

用語集

FLASH放射線治療: 超高線量率で放射線を短時間照射することで、正常組織の副作用を抑えつつ腫瘍に高い効果を期待する新しい治療法。
電離箱 (Ionisation Chamber, IC): 放射線がガスを電離して生じる電荷を測定することで、放射線量を測る装置。医療現場で広く使われる。
平行平板型電離箱 (Plane-Parallel Ionisation Chamber, PPIC): 電極が平行な平板状の電離箱で、特に電子線治療などでよく用いられるタイプ。
電荷収集効率 (Charge Collection Efficiency, CCE): 電離箱で生成された全電荷のうち、電極に収集された電荷の割合。高いほど正確な線量測定が可能。
イオン再結合 (Ion Recombination): 電離箱内で生成された陽イオンと陰イオンが、電極に収集される前に互いに結合して中性に戻る現象。線量測定の誤差につながる。
空間電荷 (Space Charge): 電離箱内のイオンや電子が蓄積されることで生じる電荷密度が、内部の電場分布を歪ませる現象。特に高線量率で顕著。
自由電子割合 (Free Electron Fraction, FEF): 電離箱内で生成された電子のうち、陰イオンに付着せず自由な状態で電極に到達する電子の割合。
モンテカルロシミュレーション (Monte Carlo Simulation): 乱数を用いて現象をシミュレートし、その統計的結果から問題を解決する計算手法。粒子輸送のモデリングに用いられる。
Garfield++: CERNで開発・保守されているオープンソースのシミュレーションツールキット。ガス中の荷電粒子の輸送と相互作用を研究するために使用される。