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

Résumé Technique : Efficacité de la Collecte de Charge dans les Chambres d'Ionisation à Plan Parallèle à Éventilation d'Air sous Débits de Dose Ultra-Élevés

Énoncé du Problème
L'irradiation à débit de dose ultra-élevé (UHDR), utilisée en radiothérapie FLASH (≥40 Gy s⁻¹), induit des effets de charge d'espace significatifs au sein des chambres d'ionisation à plan parallèle (PPIC). Ces effets entraînent des réductions substantielles de l'efficacité de collecte de charge (CCE). Bien que les modèles analytiques standards (par exemple, le modèle de Boag) tiennent compte de la recombinaison ionique dans la radiothérapie conventionnelle, ils échouent dans les conditions UHDR car ils négligent les distorsions de champ électrique causées par les hautes densités de charge. La littérature existante indique que la charge d'espace peut localement distordre le champ électrique par des facteurs dépassant quatre, pourtant un cadre auto-cohérent pour modéliser l'interaction entre le transport des particules, l'attachement électronique, la recombinaison ion-ion et l'évolution dynamique du champ faisait défaut.

Méthodologie
Les auteurs ont étendu le cadre de simulation Monte Carlo Garfield++ pour combler ces lacunes. Le développement a impliqué trois implémentations techniques primaires :

1. Recombinaison Ion-Ion : La classe AvalancheMC a été étendue pour simuler la recombinaison volumique (recombinaison entre ions provenant de différentes traces) en tant que processus poissonnien. La probabilité de recombinaison est calculée à chaque pas de temps en fonction des densités ioniques locales et d'un coefficient de recombinaison ($\alpha$) spécifié.
2. Calcul du Champ Électrique de Charge d'Espace : Un solveur de Poisson a été implémenté au sein de la classe ComponentGrid. En utilisant une approche de type Particule-en-Cellule (PIC), le code dépose la charge des particules simulées sur une grille cartésienne 3D. La densité de charge d'espace résultante ($\rho$) est utilisée pour résoudre l'équation de Poisson ($\nabla^2\phi = -\rho/\varepsilon$) via un solveur de Gauss–Seidel avec relaxation successive (SOR). Cela met à jour le champ électrique de manière dynamique, tenant compte des distorsions causées par la charge accumulée.
3. Attachement et Transport des Électrons : Le modèle couple le transport des électrons aux processus d'attachement. En utilisant Magboltz, les paramètres de nuage d'électrons (mobilité et coefficients d'attachement) sont déterminés pour des mélanges gazeux spécifiques (air avec humidité variable). La simulation suit la conversion des électrons libres en ions négatifs, ce qui est fortement dépendant de l'intensité locale du champ électrique.

La simulation gère des mélanges gazeux arbitraires, des conditions atmosphériques et des structures de faisceau. Pour gérer les coûts computationnels liés aux doses UHDR (impliquant $\sim 10^9$ particules physiques), un facteur de multiplicité de particule est employé. La parallélisation est réalisée en définissant des pas de temps globaux où les dérives de particules sont calculées en parallèle, suivis d'une étape sérielle pour mettre à jour les cartes de densité et le champ électrique.

Contributions Clés

• Cadre Auto-Cohérent : L'article présente GarFLIC (Garfield++ pour les chambres d'ionisation FLASH), un outil Monte Carlo flexible qui couple de manière auto-cohérente le transport des particules, l'attachement, la recombinaison et les distorsions dynamiques du champ électrique.
• Validation par la Littérature : Le modèle a été rigoureusement validé par rapport à :
  • Fenwick et Kumar : Reproduction de la CCE et de la fraction d'électrons libres (FEF) en l'absence de charge d'espace.
  • Kranzer et al. : Reproduction de l'évolution temporelle des distorsions de champ électrique sous des impulsions UHDR.
  • Paz Martín et al. : Reproduction de l'évolution de la FEF, de la CCE en fonction de la dose par impulsion (DPP) et des profils de courant induit.
• Aperçu Mécaniste : L'étude isole les moteurs physiques spécifiques de la réduction de la CCE dans les régimes UHDR, en distinguant le temps de recombinaison des dynamiques d'attachement.

Résultats

• Distorsion du Champ Électrique : Les simulations confirment que sous des conditions UHDR (par exemple, une impulsion de 2 Gy sur 2,5 µs), la charge d'espace peut localement augmenter le champ électrique de plus d'un facteur quatre ou le réduire à presque zéro. Ces distorsions persistent tout au long de la phase de dérive de la charge, même après la collecte des électrons.
• Fraction d'Électrons Libres (FEF) : Il est démontré que la FEF est une variable dynamique plutôt qu'une constante. À mesure que la DPP augmente, la réduction du champ électrique dans les larges régions de la chambre augmente le coefficient d'attachement électronique, provoquant une diminution progressive de la FEF.
• Efficacité de la Collecte de Charge (CCE) : Le modèle reproduit fidèlement la chute de la CCE observée dans les études de référence à travers diverses tensions de polarisation et DPP.
• Mécanisme Primaire de Perte de CCE : Les simulations indiquent que la réduction de la CCE sous UHDR n'est pas principalement pilotée par une augmentation du temps de coexistence des ions positifs et négatifs (ce qui augmenterait la recombinaison par le temps). Au contraire, la réduction est principalement pilotée par la diminution de la FEF. La distorsion du champ électrique induite par la charge d'espace abaisse l'intensité locale du champ, augmentant le coefficient d'attachement des électrons aux molécules d'oxygène. Cela génère une densité plus élevée d'ions négatifs, ce qui accentue par la suite les taux de recombinaison ion-ion.

Signification et Revendications
L'article affirme que le cadre développé constitue un outil robuste, flexible et auto-cohérent pour étudier les phénomènes de recombinaison dans les chambres d'ionisation sous conditions FLASH. En démontant que la réduction de la CCE est largement gouvernée par l'évolution de la FEF (pilotée par l'attachement dépendant du champ), les auteurs suggèrent que le problème complexe de la dosimétrie UHDR peut être simplifié. Cette découverte ouvre la voie au développement de modèles analytiques améliorés ou de méthodes de correction en temps réel pour les PPIC, permettant potentiellement de dépasser les limites des formalismes actuels basés sur Boag. Les auteurs précisent explicitement que ces mécanismes sont désormais intégrés dans le logiciel open-source Garfield++, disponible pour la recherche future et le développement de modèles dosimétriques.