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

Resumen Técnico: Eficiencia de Recolección de Carga en Cámaras de Ionización de Plano Paralelo con Ventilación de Aire bajo Tasas de Dosis Ultra-Altas

Planteamiento del Problema
La irradiación de tasa de dosis ultra-alta (UHDR), utilizada en la radioterapia FLASH (≥40 Gy s⁻¹), induce efectos significativos de carga espacial dentro de las cámaras de ionización de plano paralelo (PPIC). Estos efectos conducen a reducciones sustanciales en la eficiencia de recolección de carga (CCE). Los modelos analíticos estándar (por ejemplo, el modelo de Boag) que dan cuenta de la recombinación en radioterapia convencional, fallan bajo condiciones de UHDR porque omiten las distorsiones del campo eléctrico causadas por altas densidades de carga. La literatura existente indica que la carga espacial puede distorsionar localmente el campo eléctrico por factores que superan el cuatro, pero se carecía de un marco autoconsistente para modelar la interacción entre el transporte de partículas, la captura de electrones, la recombinación ion-ión y la evolución dinámica del campo.

Metodología
Los autores extendieron el marco de simulación de Monte Carlo Garfield++ para abordar estas brechas. El desarrollo involucró tres implementaciones técnicas primarias:

1. Recombinación Ion–Ion: Se extendió la clase AvalancheMC para simular la recombinación volumétrica (recombinación entre iones de diferentes trazas) como un proceso Poissoniano. La probabilidad de recombinación se calcula en cada paso de tiempo basándose en las densidades locales de iones y un coeficiente de recombinación ($\alpha$) especificado.
2. Cálculo del Campo Eléctrico de Carga Espacial: Se implementó un resolvedor de Poisson dentro de la clase ComponentGrid. Utilizando un enfoque de Partícula en Celda (PIC), el código deposita la carga de las partículas simuladas en una rejilla cartesiana 3D. La densidad de carga espacial resultante ($\rho$) se utiliza para resolver la ecuación de Poisson ($\nabla^2\phi = -\rho/\varepsilon$) mediante un resolvedor de Gauss–Seidel con técnica de Sobre-Relajación Sucesiva (SOR). Esto actualiza el campo eléctrico dinámicamente, contabilizando las distorsiones causadas por la carga acumulada.
3. Captura de Electrones y Transporte: El modelo acopla el transporte de electrones con los procesos de captura. Utilizando Magboltz, se determinan los parámetros de la nube de electrones (movilidad y coeficientes de captura) para mezclas de gases específicas (aire con humedad variable). La simulación rastrea la conversión de electrones libres en iones negativos, lo cual es altamente dependiente de la fuerza del campo eléctrico local.

La simulación gestiona mezclas de gases arbitrarias, condiciones atmosféricas y estructuras de haz. Para gestionar los costos computacionales para dosis UHDR (que involucran $\sim 10^9$ partículas físicas), se emplea un factor de multiplicidad de partícula. La paralelización se logra definiendo pasos de tiempo globales donde los desplazamientos de las partículas se computan en paralelo, seguidos de un paso serial para actualizar los mapas de densidad y el campo eléctrico.

Contribuciones Clave

• Marco Autoconsistente: El artículo presenta GarFLIC (Garfield++ para Cámara de Ionización FLASH), una herramienta de Monte Carlo flexible que acopla de manera autoconsistente el transporte de partículas, la captura, la recombinación y las distorsiones dinámicas del campo eléctrico.
• Validación Frente a la Literatura: El modelo fue rigurosamente validado contra:
  • Fenwick y Kumar: Reproduciendo la CCE y la fracción de electrones libres (FEF) en ausencia de carga espacial.
  • Kranzer et al.: Reproduciendo la evolución temporal de las distorsiones del campo eléctrico bajo pulsos de UHDR.
  • Paz Martín et al.: Reproduciendo la evolución de la FEF, la CCE como función de la dosis por pulso (DPP) y los perfiles de corriente inducida.
• Perspectiva Mecanística: El estudio aísla los impulsores físicos específicos de la reducción de la CCE en regímenes de UHDR, distinguiendo entre el tiempo de recombinación y la dinámica de captura.

Resultos

• Distorsión del Campo Eléctrico: Las simulaciones confirman que bajo condiciones de UHDR (por ejemplo, un pulso de 2 Gy sobre 2.5 µs), la carga espacial puede aumentar localmente el campo eléctrico en más de un factor de cuatro o reducirlo casi a cero. Estas distorsiones persisten durante toda la fase de deriva de la carga, incluso después de la recolección de electrones.
• Fracción de Electrones Libres (FEF): Se demuestra que la FEF es una variable dinámica en lugar de una constante. A medida que el DPP aumenta, la reducción del campo eléctrico en regiones grandes de la cámara incrementa el coeficiente de captura de electrones, causando una disminución progresiva en la FEF.
• Eficiencia de Recolección de Carga (CCE): El modelo reproduce con precisión la caída en la CCE observada en estudios de referencia a través de varios voltajes de polarización y DPPs.
• Mecanismo Primario de Pérdida de CCE: Las simulaciones indican que la reducción en la CCE bajo UHDR no es impulsada primordialmente por un aumento en el tiempo de coexistencia de iones positivos y negativos (lo que aumentaría la recombinación por tiempo). En su lugar, la reducción es impulsada principalmente por la disminución de la FEF. La distorsión del campo eléctrico inducida por la carga espacial reduce la fuerza del campo local, aumentando la probabilidad de captura de electrones por moléculas de oxígeno. Esto genera una mayor densidad de iones negativos, lo que subsecuentemente potencia las tasas de recombinación ion–ión.

Significancia y Reivindicaciones
El artículo afirma que el marco desarrollado proporciona una herramienta robusta, flexible y autoconsistente para investigar fenómenos de recombinación en cámaras de ionización bajo condiciones FLASH. Al demostrar que la reducción de la CCE está gobernada en gran medida por la evolución de la FEF (impulsada por la captura dependiente del campo), los autores sugieren que el complejo problema de la dosimetría UHDR puede simplificarse. Este hallazgo abre la posibilidad de desarrollar modelos analíticos mejorados o métodos de corrección en tiempo real para las PPIC, moviéndose potencialmente más allá de las limitaciones de los actuales formalismos basados en Boag. Los autores declaran explícitamente que estos mecanismos han sido integrados en el software de código abierto Garfield++, disponible para futuras investigaciones y desarrollo de modelos dosimétricos.