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

Resumo Técnico: Eficiência de Recolha de Carga em Câmaras de Ionização de Plano-Paralelo com Ventilação de Ar a Taxas de Dose Ultra-Altas

Definição do Problema
A irradiação com taxas de dose ultra-altas (UHDR), utilizada na radioterapia FLASH (≥40 Gy s⁻¹), induz efeitos significativos de carga espacial em câmaras de ionização de plano-paralelo (PPIC). Estes efeitos levam a reduções substanciais na eficiência de recolha de carga (CCE). Os modelos analíticos padrão (por exemplo, o modelo de Boag) que explicam a recombinação em radioterapia convencional, falham sob condições de UHDR porque negligenciam as distorções do campo elétrico causadas por altas densidades de carga. A literatura existente indica que a carga espacial pode distorcer localmente o campo elétrico em fatores superiores a quatro, contudo, carecia de um quadro auto-consistente para modelar a interação entre o transporte de partículas, a anexação de eletrões, a recombinação ião–ião e a evolução dinâmica do campo.

Metodologia
Os autores estenderam o framework de simulação Monte Carlo Garfield++ para colmatar estas lacunas. O desenvolvimento envolveu três implementações técnicas primárias:

1. Recombinação Ião–Ião: A classe AvalancheMC foi estendida para simular a recombinação volumétrica (recombinação entre iões de diferentes traços) como um processo Poissoniano. A probabilidade de recombinação é calculada em cada passo de tempo com base nas densidades locais de iões e num coeficiente de recombinação ($\alpha$) especificado.
2. Cálculo do Campo Elétrico de Carga Espacial: Um solver de Poisson foi implementado dentro da classe ComponentGrid. Utilizando uma abordagem de Partícula-em-Célula (PIC), o código deposita carga de partículas simuladas numa grelha cartesiana 3D. A densidade de carga espacial resultante ($\rho$) é utilizada para resolver a equação de Poisson ($\nabla^2\phi = -\rho/\varepsilon$) através de um solver Gauss–Seidel com Sobrerrelaxação Sucessiva (SOR). Isto atualiza o campo elétrico dinamicamente, contabilizando as distorções causadas pela carga acumulada.
3. Anexação e Transporte de Eletrões: O modelo acopla o transporte de eletrões com processos de anexação. Utilizando o Magboltz, os parâmetros do enxame de eletrões (mobilidade e coeficientes de anexação) são determinados para misturas gasosas específicas (ar com humidade variável). A simulação rastreia a conversão de eletrões livres em iões negativos, o que é altamente dependente da força do campo elétrico local.

A simulação lida com misturas gasosas arbitrárias, condições atmosféricas e estruturas de feixe. Para gerir os custos computacionais para doses de UHDR (envolvendo $\sim 10^9$ partículas físicas), é empregado um fator de multiplicidade de peso da partícula. A paralelização é alcançada definindo passos de tempo globais onde os drifts de partículas são computados em paralelo, seguidos de um passo serial para atualizar os mapas de densidade e o campo elétrico.

Principais Contribuições

• Framework Auto-consistente: O artigo apresenta o GarFLIC (Garfield++ for FLASH Ionisation Chamber), uma ferramenta Monte Carlo flexível que acopla auto-consistentemente o transporte de partículas, a anexação, a recombinação e as distorções dinâmicas do campo elétrico.
• Validação Contra a Literatura: O modelo foi rigorosamente validado contra:
  • Fenwick e Kumar: Reproduzindo a CCE e a fração de eletrões livres (FEF) na ausência de carga espacial.
  • Kranzer et al.: Reproduzindo a evolução temporal das distorções do campo elétrico sob pulsos de UHDR.
  • Paz Martín et al.: Reproduzindo a evolução da FEF, a CCE como função da dose por pulso (DPP) e perfis de corrente induzida.
• Insight Mecanístico: O estudo isola os drivers físicos específicos da redução da CCE em regimes de UHDR, distinguindo entre o tempo de recombinação e a dinâmica de anexação.

Resultos

• Distorção do Campo Elétrico: As simulações confirmam que, sob condições de UHDR (por exemplo, um pulso de 2 Gy sobre 2,5 µs), a carga espacial pode aumentar localmente o campo elétrico em mais de um fator de quatro ou reduzi-lo a quase zero. Estas distorções persistem durante toda a fase de drift de carga, mesmo após a recolha de eletrões.
• Fração de Eletrões Livres (FEF): Demonstra-se que a FEF é uma variável dinâmica e não uma constante. À medida que a DPP aumenta, a redução do campo elétrico em grandes regiões da câmara aumenta o coeficiente de anexação de eletrões, causando uma diminuição progressiva na FEF.
• Eficiência de Recolha de Carga (CCE): O modelo reproduz com precisão a queda na CCE observada em estudos de referência através de várias tensões de polarização e DPPs.
• Mecanismo Primário de Perda de CCE: As simulações indicam que a redução na CCE sob UHDR não é impulsionada primariamente por um aumento no tempo de coexistência de iões positivos e negativos (o que aumentaria a recombinação via tempo). Em vez disso, a redução é principalmente impulsionada pela diminuição da FEF. A distorção do campo elétrico induzida pela carga espacial reduz a força do campo local, aumentando o coeficiente de anexação de eletrões ao oxigénio molecular. Isto gera uma maior densidade de iões negativos, o que subsequentemente aumenta as taxas de recombinação ião–ião.

Significância e Alegações
O artigo alega que o framework desenvolvido fornece uma ferramenta robusta, flexível e auto-consistente para investigar fenómenos de recombinação em câmaras de ionização sob condições de FLASH. Ao demonstrar que a redução da CCE é largamente governada pela evolução da FEF (impulsionada pela anexação dependente do campo), os autores sugerem que o complexo problema da dosimetria de UHDR pode ser simplificado. Esta descoberta abre a possibilidade para o desenvolvimento de modelos analíticos melhorados ou métodos de correção em tempo real para PPICs, potencialmente indo além das limitações dos atuais formalismos baseados em Boag. Os autores declaram explicitamente que estes mecanismos foram integrados no software de código aberto Garfield++, disponível para futuras investigações e desenvolvimento de modelos dosimétricos.