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

Sintesi Tecnica: Efficienza di Raccolta della Carica nelle Camere di Ionizzazione Planari ad Apertura d'Aria a Tassi di Dose Ultra-Elevati

Problema
L'irradiazione a tasso di dose ultra-elevato (UHDR), utilizzata nella radioterapia FLASH (≥40 Gy s⁻¹), induce significativi effetti di carica spaziale all'interno delle camere di ionizzazione planari (PPIC). Questi effetti portano a riduzioni sostanziali dell'efficienza di raccolta della carica (CCE). I modelli analitici standard (ad es. il modello di Boag) tengono conto della ricombinazione ionica nelle radioterapie convenzionali, ma falliscono in condizioni UHDR perché trascurano le distorsioni del campo elettrico causate dalle alte densità di carica. La letteratura esistente indica che la carica spaziale può distorcere localmente il campo elettrico con fattori superiori a quattro, eppure mancava un quadro auto-coerente per modellare l'interazione tra trasporto di particelle, cattura elettronica (attachment), ricombinazione iono-ionica ed evoluzione dinamica del campo.

Metodologia
Gli autori hanno esteso il framework di simulazione Monte Carlo Garfield++ per colmare queste lacune. Lo sviluppo ha comportato tre implementazioni tecniche primarie:

1. Ricombinazione Iono-Ionica: La classe AvalancheMC è stata estesa per simulare la ricombinazione volumetrica (ricombinazione tra ioni provenienti da tracce diverse) come un processo Poissoniano. La probabilità di ricombinazione è calcolata ad ogni passo temporale sulla base delle densità ioniche locali e di un coefficiente di ricombinazione specificato ($\alpha$).
2. Calcolo del Campo Elettrico da Carica Spaziale: Un risolutore di Poisson è stato implementato all'interno della classe ComponentGrid. Utilizzando un approccio Particle-in-Cell (PIC), il codice deposita la carica dalle particelle simulate su una griglia cartesiana 3D. La densità di carica spaziale risultante ($\rho$) viene utilizzata per risolvere l'equazione di Poisson ($\nabla^2\phi = -\rho/\varepsilon$) tramite un risolutore Gauss–Seidel iterativo con sovrarilassamento successivo (SOR). Questo aggiorna dinamicamente il campo elettrico, tenendo conto delle distorsioni causate dalla carica accumulata.
3. Cattura Elettronica e Trasporto: Il modello accoppia il trasporto elettronico con i processi di cattura (attachment). Utilizzando Magboltz, i parametri dello sciame elettronico (mobilità e coefficienti di cattura) vengono determinati per specifiche miscele gassose (aria con diversa umidità). La simulazione traccia la conversione di elettroni liberi in ioni negativi, un processo altamente dipendente dall'intensità del campo elettrico locale.

La simulazione gestisce miscele gassose arbitrarie, condizioni atmosferiche e strutture del fascio. Per gestire i costi computazionali per dosi UHDR (che comportano $\sim 10^9$ particelle fisiche), viene impiegato un fattore di molteplicità della particella (peso della particella). La parallelizzazione è ottenuta definendo passi temporali globali in cui i drift delle particelle sono calcolati in parallelo, seguiti da un passaggio seriale per aggiornare le mappe di densità e il campo elettrico.

Contributi Chiave

• Framework Auto-Coerente: Il documento presenta GarFLIC (Garfield++ per la Camera di Ionizzazione FLASH), uno strumento Monte Carlo flessibile che accoppia auto-coerentemente trasporto di particelle, cattura, ricombinazione e distorsioni dinamiche del campo elettrico.
• Validazione rispetto alla Letteratura: Il modello è stato rigorosamente validato contro:
  • Fenwick e Kumar: Riproducendo la CCE e la frazione di elettroni liberi (FEF) in assenza di carica spaziale.
  • Kranzer et al.: Riproducendo l'evoluzione temporale delle distorsioni del campo elettrico sotto impulsi UHDR.
  • Paz Martín et al.: Riproducendo l'evoluzione della FEF, la CCE in funzione della dose per impulso (DPP) e i profili di corrente indotta.
• Approfondimento Meccanicistico: Lo studio isola i driver fisici specifici della riduzione della CCE nei regimi UHDR, distinguendo tra il tempo di ricombinazione e la dinamica di cattura.

Risultati

• Distorsione del Campo Elettrico: Le simulazioni confermano che, in condizioni UHDR (ad es. un impulso di 2 Gy in 2,5 µs), la carica spaziale può aumentare localmente il campo elettrico di oltre un fattio quattro o ridurlo quasi a zero. Queste distorsioni persistono durante tutta la fase di drift della carica, anche dopo la raccolta degli elettroni.
• Frazione di Elettroni Liberi (FEF): La FEF è dimostrata essere una variabile dinamica piuttosto che una costante. All'aumentare del DPP, la riduzione del campo elettrico nelle grandi regioni della camera aumenta il coefficiente di cattura elettronica, causando una progressiva diminuzione della FEF.
• Efficienza di Raccolta della Carica (CCE): Il modello riproduce accuratamente il calo della CCE osservato negli studi di riferimento attraverso varie tensioni di polarizzazione e DPP.
• Meccanismo Primario di Perdita della CCE: Le simulazioni indicano che la riduzione della CCE sotto UHDR non è guidata principalmente da un aumento del tempo di coesistenza di ioni positivi e negativi (il che aumenterebbe la ricombinazione tramite il tempo). Invece, la riduzione è principalmente guidata dalla diminuzione della FEF. La distorsione del campo elettrico indotta dalla carica spaziale abbassa l'intensità del campo locale, aumentando il coefficiente di cattura elettronica alle molecole di ossigeno. Ciò genera una maggiore densità di ioni negativi, che successivamente potenzia i tassi di ricombinazione iono-ionica.

Significatività e Rivendicazioni
Il documento afferma che il framework sviluppato fornisce uno strumento robusto, flessibile e auto-coerente per investigare i fenomeni di ricombinazione nelle camere di ionizzazione in condizioni FLASH. Dimostrando che la riduzione della CCE è ampiamente governata dall'evoluzione della FEF (guidata dalla cattura dipendente dal campo), gli autori suggeriscono che il complesso problema della dosimetria UHDR possa essere semplificato. Questa scoperta apre la possibilità di sviluppare modelli analitici migliorati o metodi di correzione in tempo reale per le PPIC, andando potenzialmente oltre i limiti degli attuali formalismi basati su Boag. Gli autori dichiarano esplicitamente che questi meccanismi sono ora integrati nel software open-source Garfield++, disponibile per la ricerca futura e lo sviluppo di modelli dosimetrici.