Page 15 - PR104

Especificaciones de la simulación de Monte Carlo Numerosos parámetros de entrada son esenciales para una correcta simulación con MCNP6, tales como, la definición detallada del transporte físico, tanto de fotones como de electrones incluyendo procesos de dispersion y de absorción, la especificación de los materiales que atraviesan las partícu- las y el empleo de técnicas de reducción de la varianza. Para seguir la trayectoria tanto de fotones como de elec- trones hay que activar tanto MODE P como E. En ambos casos, con el objetivo de que la incertidumbre en los pun- tos evaluados sea menor del 3 %, se simulan al menos 108 partículas iniciales. Los electrones cuya energía está por debajo de 1 keV no se rastrean, sin embargo, se considera que depositan la energía localmente. Para acelerar los cálculos sin comprometer la precision de los resultados, el código MCNP6 (version 6.1.1) se ejecuta en un cluster Quasar del grupo de investigación SENUBIO ISIRYM usando el protocolo MPI con 32 procesadores en paralelo. Los detalles definidos en el archivo de entrada del sof- tware MCNP6 se recogen en la Tabla I, siguiendo las reco- mendaciones de AAPM TG-268 Report \[9\]. Modelado del haz de electrones La simulación del modelo con MCNP6 incluye el transporte del haz de electrones dirigido desde el cabezal del acelerador lineal móvil, atravesando el aplicador hasta llegar al maniquí. El parámetro principal que determina el espectro ener- gético de la fuente de electrones, considerando que el perfil de intensidad sigue una distribución gaussiana, es la energía media de los electrones y FWHM, el ancho del perfil a la mitad de la altura máxima, 12 MeV y 10 %, respectivamente. Registros de resultados (tallies) El código de Monte Carlo se basa en un proceso estadís- tico iterativo para estimar el transporte de rutas aleatorias y las interacciones de fotones y electrones que son emi- tidos desde una fuente tridimensional. Los resultados de dichas interacciones entre partículas pueden ser recogidos mediante tarjetas de registro denominadas tallies. En este modelo, para obtener el flujo de partículas por cm2, se incluye en el archivo de entrada de MCNP6 una tarjeta específica de registro, llamada FMESH4. Este tally proporciona el flujo de unidades de #/ cm2·s para cada partícula emitida. Para pasar los resultados de flujo a dosis media se introduce un factor de conversion, de forma que los resultados se obtienen en Gy/s para cada punto del mallado. El coeficiente de conversión utilizado es el poder de frenado para electrones y el coeficiente de absorción másico (men/r) para fotones (datos obtenidos del NIST, National Institute of Standards and Technology) \[10\]. RESULTADOS Con el objetivo de validar los resultados obtenidos, se lleva a cabo un análisis de los estudios disponibles sobre dosimetría de aceleradores lineales para radioterapia intraoperatoria. Sin embargo, dado que estos aceleradores presentan unas caracteríasticas particulares y diferentes a las de los aceleradores convencionales, la información dis- ponible es limitada. Por esta razón, los datos experimenta- les se extraen de un estudio de dosimetría realizado por la Universidad Católica de Murcia \[5\] sobre el Mobetron, dada la similitud en la metodología utilizada. En este estudio \[5\], se calibra el acelerador mediante dosis absorbida en agua usando una cámara de ionización pla- no-paralela de tipo ROOS indicada en el protocolo TRS-398 de la IAEA \[11\]. Este protocolo recomienda para aceleradores móviles realizar las medidas en las condiciones más desfa- vorables posibles, es decir, a la máxima energía disponible y con el aplicador que tenga mayor diámetro, en este caso 12 MeV y 10 cm, respectivamente. Teniendo en cuenta estas especificaciones, la dosis absorbida en agua de refe- rencia para el Mobetron es de 0.966 cGy/MU. Colaboraciones Código MCNP6 Versión 1.1 compilado con Intel v12, con indicador de optimización O1 Validación Previamente validado Duración 10 h (tiempo total de la CPU teniendo en cuenta los 32 procesadores). Todas las simulaciones se llevan a cabo con 108 historias Descripción de la fuente Haz de electrones: distribución de electrones gaussiana (media:12 MeV, FWHM:10 %) Sección transversal MCNPLIB84 EL03 Parámetros de transporte Energía límite de fotones = 1 keV Energía límite de electrones = 10 keV Herramientas de reducción de varianza (a) Sesgo de fuente (en dirección) (b) Ponderación de celdas (por importancias) Magnitudes medidas Dosis absorbida en aire Incertidumbre estadística ≤3 % (k = 2, incertidumbre máxima de todas las magnitudes medidas) Tabla I. Resumen de las principales características del método de Monte Carlo utilizado. DOSIS DISPERSA EN RADIOTERAPIA INTRAOPERATORIA 15