Skip to main content

La integración de la resonancia magnética con aceleradores lineales representa uno de los avances más significativos en la radioterapia guiada por imágenes de la última década. Sistemas como Elekta Unity (1,5 T) y ViewRay MRIdian (0,35 T) le permiten visualizar tumores y estructuras críticas con resolución de tejido blando en tiempo real durante la irradiación. Esta capacidad abre perspectivas concretas para la adaptación del plan diario (radioterapia adaptativa en línea) y para reducir los márgenes de planificación en sitios de alta movilidad como el pulmón, la próstata y el páncreas.

Sin embargo, la presencia del campo magnético estático introduce perturbaciones físicas que los algoritmos configurados para entornos libres de campo no representan. El campo modifica las trayectorias de los electrones secundarios, altera la distribución de dosis en las interfaces y puede producir diferencias clínicamente relevantes en las geometrías con el aire y la baja densidad. La conclusión correcta no es que una familia algorítmica siempre sea inválida, sino que el motor utilizado en MR-Linac debe incorporar explícitamente el campo y validarse para el entorno empresarial específico.

Electron trajectories under the Lorentz force in MR-Linac
Infografía técnica del cluster de algoritmos de cálculo de dosis.

Este artículo analiza los mecanismos físicos subyacentes a estos efectos, cómo los sistemas de planificación de tratamiento (TPS) los modelan computacionalmente, los requisitos de puesta en marcha y control de calidad (QA) específicos de MR-Linac y las implicaciones prácticas para los sitios anatómicos tratados con mayor frecuencia en esta plataforma. Las distinciones entre descripción física establecida, implementación comercial y evidencia de validación clínica se dejan claras a lo largo del texto.


Por qué agregar un campo magnético cambia la dosimetría

En los aceleradores convencionales, la deposición de dosis mediante un haz de fotones está mediada principalmente por electrones secundarios. Las implementaciones clínicas diseñadas únicamente para campo magnético cero no incluyen la fuerza de Lorentz en el transporte electrónico. Por lo tanto, el nombre genérico del algoritmo no es suficiente: es necesario confirmar si la versión de TPS fue desarrollada y puesta en servicio para la geometría magnética del equipo.

Cuando un campo magnético estático B se superpone al volumen irradiado, la fuerza F = q(v × B) desvía cada electrón cargado de su trayectoria lineal. Para los electrones de energía típica en la terapia de fotones (2-20 MeV), el radio de curvatura en el campo de 1,5 T de Unity es del orden de centímetros, comparable al rango adecuado de estos electrones en el tejido blando. El resultado es que la distribución de la dosis ya no es simétrica axialmente y comienza a depender de la orientación entre el campo B, el haz de fotones y las interfaces materiales presentes en la geometría del paciente. Los fotones, al ser neutros, no se desvían; sólo se perturba el transporte electrónico secundario.

Tres consecuencias prácticas surgen directamente de este hecho. Primero, incluso en un campo homogéneo, la desviación de los electrones en el plano perpendicular a B crea asimetría en los perfiles del haz, observable en mediciones con cámaras o películas de ionización. En segundo lugar, en los límites entre el tejido y el aire (pared torácica sobre el pulmón, mucosa esofágica, superficie de la piel), el cambio abrupto en la densidad electrónica amplifica la desviación, produciendo una sobredosis o una subdosificación localizada según la geometría y la orientación. En tercer lugar, surge el llamado efecto de retorno de electrones (ERE), descrito en la siguiente sección, que es físicamente distinto de la mera asimetría del perfil.

La revisión de Sim et al. (2020) sistematiza los desafíos físicos en la radioterapia guiada por RM, abarcando la perturbación de dosis, la instrumentación y el control de calidad. Distingue explícitamente entre lo que está físicamente bien establecido, lo que ha sido validado en condiciones clínicas y lo que sigue bajo investigación, una distinción que debería guiar la implementación en cualquier servicio.


Fuerza de Lorentz y efecto de retorno de electrones

ERE es el fenómeno más ampliamente documentado de la interacción entre el campo magnético y la dosimetría en MR-Linac. Ocurre cuando los electrones secundarios penetran en regiones de baja densidad electrónica (pulmón, aire intraluminal) y, en lugar de dispersarse progresivamente como lo harían sin campo, describen trayectorias curvas que los devuelven al tejido adyacente de alta densidad. El resultado es un pico de dosis en la superficie del tejido «receptor», por ejemplo, en la pared interna de un tumor adyacente al pulmón o en la mucosa esofágica adyacente al aire.

La configuración más común en Unity es un haz perpendicular a B (campo paralelo al eje cráneo-caudal, haz entrando lateralmente). En esta configuración, el ERE es más pronunciado porque los electrones emitidos hacia adelante desde el punto de interacción tienen un componente de velocidad en el plano donde la fuerza de Lorentz es máxima. Para un haz paralelo a B, los electrones tienden a formar espirales a lo largo de las líneas de campo, reduciendo el ERE pero favoreciendo la corriente de electrones (siguiente sección).

La magnitud del ERE depende de cuatro variables principales: intensidad del campo magnético (campo más grande → mayor curvatura → ERE más pronunciado), energía del haz (los haces de mayor energía producen electrones con mayor radio y, por lo tanto, mayor distancia de retorno), ángulo entre el haz y el campo, y geometría de la interfaz (plana, cóncava o convexa). Estas dependencias se caracterizaron en trabajo de banco, con cámaras de ionización y películas, en geometrías de referencia con fantasmas de sólido y agua.

Es fundamental distinguir la evidencia disponible en diferentes contextos. ERE como fenómeno físico en geometrías simples está bien establecido y se reproduce mediante simulaciones Monte Carlo validadas. La extrapolación cuantitativa a geometrías clínicas complejas (con tumores irregulares, múltiples interfaces y variación intrafracción) es objeto de investigación activa y presenta incertidumbres que deben tenerse en cuenta en la planificación clínica.


Flujo de electrones en cavidades y regiones fuera del campo

Cuando el haz es paralelo al campo magnético, los electrones secundarios tienden a seguir las líneas del campo en lugar de dispersarse lateralmente. Este fenómeno, llamado , produce haces de electrones estrechos y altamente localizados que cruzan cavidades (tráquea, esófago, intestino) y regiones periféricas al campo de tratamiento con una fluencia significativamente mayor de lo que se esperaría en ausencia de un campo.

tiene dos consecuencias prácticas distintas. El primero es el aumento de la dosis periférica, importante para los órganos en riesgo alejados del objetivo, especialmente en pacientes pediátricos o cuando las estructuras radiosensibles están cerca del campo. El segundo es el depósito de dosis en patrones filamentosos dentro de las cavidades de aire o gas intestinal, un patrón que los algoritmos deterministas sin campo no modelan. La detección experimental de requiere dosímetros con alta resolución espacial (películas radiocrómicas o conjuntos 2D de pequeños detectores) ya que el fenómeno se manifiesta en regiones estrechas.

La implicación práctica más importante es que los protocolos de control de calidad de MR-Linac deben incluir la verificación de dosis fuera de campo como un componente explícito. En los aceleradores convencionales, la dosis periférica fuera del campo es pequeña y rara vez se controla directamente. En MR-Linac, puede producir dosis periféricas no predichas por los algoritmos convencionales. Este punto es un escenario de fracaso real, no hipotético, y debería incluirse en la formación del equipo de físico médico.


Cómo Monte Carlo y GPUMCD modelan el problema

La solución natural para el transporte de partículas en un campo magnético es el método Monte Carlo (MC). Simuladores como EGSnrc, GATE/Geant4 y MCNP permiten incluir la fuerza de Lorentz explícitamente en el paso de transporte electrónico, agregando el término qv × B a la ecuación de movimiento en cada paso. Los códigos MC validados reproducen ERE, y asimetrías de perfiles con fidelidad adecuada a partir de datos experimentales en geometrías de referencia, lo que los convierte en el punto de referencia para validar algoritmos comerciales.

El desafío clínico es el tiempo de cálculo. La reducción de la incertidumbre estadística aumenta la cantidad de historias y el costo computacional, mientras que la adaptación en línea impone una ventana operativa corta. El rendimiento debe medirse en el hardware, la red y el protocolo de incertidumbre utilizados por el propio servicio.

GPUMCD surgió originalmente como una plataforma Monte Carlo orientada a GPU, no como un algoritmo creado específicamente para Unity. Las implementaciones comerciales posteriores pueden incorporar el campo magnético e integrar flujos de Monaco/Unity, pero esta capacidad depende del producto y de la versión. Lo mismo ocurre con la magnitud de la dosis: no debe asumir dose to medium o dose to water sin consultar la documentación y los datos exportados por la versión encargada.

Tabla 1 — Comparación de algoritmos de cálculo de dosis en el contexto de MR-Linac

Familia Tipo Requisito para MR-Linac Limitación si se utiliza sin adaptación de campo
Pencil Beam / kernel Analítico o semianalítico Modelo explícito validado con campo B Puede perder asimetría, ERE y streaming
solucionador LBTE Determinista Término de transporte compatible con geometría magnética El nombre LBTE por sí solo no garantiza el soporte para el campo
Monte Carlo en GPU Estocástico Transporte electrónico con fuerza de Lorentz y modelo de haz local El resultado depende de los materiales, la cuadrícula, las estadísticas y la independencia beam model
Monte Carlo Estocástico Código y geometría validados para el campo estudiado El costo y la integración pueden limitar el uso en línea

GPUMCD es una implementación comercial cuyas validaciones clínicas continúan acumulando evidencia. Las validaciones de geometrías de referencia son más sólidas que las validaciones de casos clínicos complejos. Los escenarios que pueden desafiar las suposiciones del modelo incluyen implantes metálicos voluminosos cerca del objetivo, anatomías con cavidades de aire de geometría extremadamente irregular y pacientes fuera del rango de variación cubierto por la puesta en servicio local.


Densidad electrónica, TC sintética y asignación de materiales

En la planificación convencional, la simulación de TC proporciona números de Hounsfield (HU), convertidos en densidad relativa de electrones mediante la tabla de calibración HU-ρe. En MR-Linac, la imagen disponible en tiempo real es una resonancia magnética, cuya señal no tiene relación física directa con la densidad electrónica. Se utilizan tres enfoques para solucionar este problema:

Asignación de densidad aparente: las estructuras contorneadas reciben valores fijos de densidad electrónica de tablas de referencia (tejidos blandos, huesos, pulmones, aire). Es el método más rápido y el más utilizado clínicamente, pero introduce errores cuando la anatomía real se desvía de la referencia: los pulmones hiperinflados frente a los atelectásicos, por ejemplo, pueden tener densidades muy diferentes del valor de referencia asignado.

CT sintética (sCT): los algoritmos de aprendizaje automático o bases de atlas generan una imagen pseudo-CT a partir de la resonancia magnética, con valores de HU aproximados píxel a píxel. Permite el cálculo de dosis con resolución espacial completa y captura heterogeneidades locales no cubiertas por la densidad aparente. Falla cuando la resonancia magnética tiene artefactos de movimiento, una fuerte susceptibilidad magnética (cerca de prótesis) o cuando la anatomía del paciente está fuera de la distribución de entrenamiento del modelo.

TC de referencia con registro deformable: el TC de planificación inicial se registra deformablemente en la resonancia magnética del día. Útil para las fracciones tempranas, pero acumula incertidumbre a medida que cambia la anatomía: un nuevo derrame pleural, una variación del llenado rectal o una pérdida de peso significativa pueden invalidar el registro.

Cada abordaje tiene implícitas hipótesis y escenarios de falla que deben documentarse en el protocolo clínico del servicio. La validación del método elegido debe ser una parte explícita del proceso de puesta en servicio, con pruebas en fantasmas con densidades conocidas. Los protocolos publicados por Cusumano et al. (2022) describen este tipo de evaluación como un requisito previo al uso clínico.


Puesta en marcha y control de calidad específicos de MR-Linac

La puesta en marcha de un MR-Linac tiene todas las complejidades de un acelerador convencional, además de desafíos específicos del entorno del campo magnético. El protocolo descrito por Winkel et al. (2019) para Unity cubre todo, desde mediciones de salida y perfiles de haz hasta verificación de dosis absoluta en fantasmas con un campo magnético activo, y es un punto de partida para el desarrollo de protocolos locales.

Los desafíos instrumentales clave incluyen cámaras de ionización, cuyo factor de respuesta se ve afectado por el campo magnético de una manera que depende de la orientación de la cámara en relación con B , y conjuntos de detectores, que deben calificarse específicamente para el campo antes de su uso clínico. Las películas radiocrómicas (EBT3, EBT-XD) son los dosímetros más robustos para la verificación de la distribución 2D, ya que no se ven afectados directamente por el campo magnético. Los dosímetros MOSFET y TLD requieren validación individual en el entorno de campo antes de cualquier uso clínico.

Tabla 2 — Componentes de control de calidad específicos de MR-Linac y sus particularidades

Componente de control de calidad Desafío específico de MR-Linac Nota
Salida absoluta (cGy/UM) Factor de corrección de la cámara en el campo B (dependiente de la orientación) Verifique con el factor kB determinado en la puesta en servicio
Perfiles de haz Asimetría inducida por B; datos de verificación versus datos de referencia Campo B activo durante la medición
Coincidencia del isocentro RT/MRI El isocentro mecánico debe coincidir con el isocentro de la imagen Herramienta dedicada (fantasma Winston-Lutz adaptado)
Dosis fuera de campo Transmisión de electrones sin equivalente en linac convencional Se requieren dosímetros de alta resolución espacial
Asignación de materiales (sCT/bulk) Validación en fantasma con densidades conocidas Repetir después de TPS actualizaciones de software
Verificación del plan individual Dosis en fantasma con campo B activo Formato (2D, 3D, punto) definido en el protocolo local
Geometría de la imagen de resonancia magnética La distorsión del gradiente afecta la ubicación del objetivo Fantasma de distorsión geométrica; frecuencia mensual
Control de calidad del proceso de adaptación Registro + recontorneado + reoptimización: cada paso con fuente de error Auditoría periódica de la cadena completa

Un aspecto a menudo subestimado es el control de calidad del proceso de adaptación en línea como una cadena integrada. Verificar únicamente componentes aislados (salida, perfiles, geometría) no garantiza que el producto final (el plano adaptado entregado al paciente) sea correcto. La cadena de adaptación implica el registro de imágenes, el recontorneado automático, la reoptimización del plan y el recálculo de dosis; cualquier paso puede introducir errores no detectados por los controles individuales.


Implicaciones para pulmón, esófago, mama y pelvis

La relevancia clínica de los efectos dosimétricos del campo magnético es heterogénea entre sitios anatómicos, siendo proporcional a la presencia y variabilidad de las interfaces aire-tejido.

El pulmón es un sitio particularmente sensible. La baja densidad del parénquima aumenta la relevancia de las interfaces y la falta de equilibrio electrónico. La precisión requiere de un motor que modele el campo magnético y haya sido validado en esta geometría; aplicar una implementación destinada al campo nulo no es un reemplazo aceptable.

El esófago presenta el desafío adicional de que el volumen de aire intraluminal varía entre fracciones: deglución, reflujo, posicionamiento. El ERE en la pared esofágica adyacente al aire crea picos de dosis locales, pero la posición exacta de estos picos cambia con la geometría de la cavidad. La adaptación en línea permite monitorear esta variabilidad, que sería invisible para los sistemas IGRT convencionales; sin embargo, es difícil predecir con precisión la dosis aplicada a la mucosa si el contorno de la cavidad de aire no se actualiza con cada fracción.

Mama incluye interfaces tejido-aire en la superficie de la piel donde ERE puede resultar en una dosis superficial elevada en ciertas orientaciones de campo y pórtico. En regímenes de hipofracción, donde la dosis por fracción es mayor y la tolerancia a la dosis cutánea es un factor limitante, este efecto tiene implicaciones clínicas directas. La dosimetría de superficie en MR-Linac requiere dosímetros específicos (películas extrabolus, detectores de superficie delgada) que deben formar parte del protocolo de puesta en servicio.

La pelvis (próstata, cuello uterino, recto) tiene una proporción menor de aire que los pulmones y el esófago, lo que reduce la magnitud del ERE en relación a estos sitios. Sin embargo, las interfaces con el recto y la vejiga siguen siendo relevantes, especialmente cuando hay gas rectal variable entre las fracciones. El principal beneficio clínico de MR-Linac en la pelvis es la visibilidad superior del tejido blando en la resonancia magnética, lo que permite una delimitación precisa de la próstata, las vesículas seminales y las estructuras de los ganglios linfáticos sin contraste. La evidencia de validación dosimétrica es más sólida para la pelvis que para el pulmón y el esófago, lo que refleja tanto el número de centros que tratan la próstata en MR-Linac como la menor complejidad dosimétrica del sitio.

La aplicación de MR-Linac en sitios menos estudiados, como páncreas, hígado y tumores de cabeza y cuello, requiere precaución adicional e idealmente una recopilación prospectiva de datos dosimétricos como parte de protocolos de investigación formales.


Preguntas frecuentes

¿Se puede utilizar GPUMCD para todos los casos clínicos de Unity sin restricciones?

La implementación Monte Carlo prevista para el flujo de equipos clínicos incluye el modelado del campo magnético, pero sus límites siguen definidos por producto, versión y puesta en servicio. Los implantes, las cavidades complejas y las anatomías fuera del conjunto de prueba pueden requerir una verificación adicional. El protocolo local debe registrar estos límites sin asumir que el nombre GPUMCD garantiza un comportamiento idéntico en todas las plataformas.

¿Cuál es la diferencia práctica entre dose to medium (D~m~) y dose to water (D~w~) en el contexto de MR-Linac?

La convención reportada debe ser confirmada en la versión clínica. D~m~ y D~w~ no son etiquetas intercambiables, especialmente en materiales cuya composición difiere de la del agua. El problema práctico surge cuando las restricciones históricas y la distribución calculada utilizan convenciones diferentes; esta relación debe documentarse en el proceso de implementación.

¿La TC sintética generada a partir de resonancia magnética es lo suficientemente precisa como para reemplazar la TC de planificación en todos los sitios?

No universalmente. La precisión de la sCT varía según el sitio anatómico, la calidad de la imagen y el método de generación (aprendizaje automático, atlas o densidad aparente). Para la pelvis, los resultados publicados son en general favorables. En el caso del pulmón, la variabilidad de la densidad parenquimatosa y la mayor sensibilidad a ERE hacen que la precisión de la sCT sea más crítica. La validación del método elegido debe realizarse localmente, con fantasmas de densidad conocida, y no simplemente asumirse a partir de publicaciones de otros centros con hardware y software potencialmente diferentes.

¿Qué dosímetros son adecuados para la verificación de dosis en un campo magnético?

Las películas radiocrómicas (EBT3, EBT-XD) son la opción más robusta para distribuciones 2D ya que no se ven afectadas directamente por el campo magnético estático. Las cámaras de ionización se pueden utilizar para dosis puntuales con factores de corrección específicos de la orientación de la cámara en relación con B ; estos factores deben determinarse en la puesta en servicio, no asumirse a partir de publicaciones con diferentes geometrías. Los dosímetros MOSFET y TLD requieren validación de campo individual antes de su uso clínico. Los conjuntos de detectores deben estar específicamente calificados; No todos los modelos comerciales son adecuados sin modificaciones o factores de corrección documentados por el fabricante.

¿MR-Linac elimina la necesidad de control de calidad de imágenes convencional?

No. MR-Linac agrega requisitos de control de calidad relacionados con el campo magnético y la cadena de adaptación en línea, pero no elimina los requisitos de control de calidad de imágenes. La distorsión geométrica de la resonancia magnética por los gradientes y la susceptibilidad magnética es una fuente de error relevante: si el tumor aparece en una posición desplazada por la distorsión del gradiente, el tratamiento estará mal centrado independientemente de la calidad del cálculo de la dosis. El control de calidad de la geometría de la imagen de resonancia magnética (distorsión, resolución, contraste) es un componente obligatorio del programa de control de calidad. Además, el control de calidad del proceso adaptativo como cadena integrada (registro + recontorneado + reoptimización + recálculo) debe ser auditado periódicamente, ya que pueden surgir errores en la integración de los pasos incluso cuando cada componente individual está dentro de las especificaciones.


Referencias

  1. Sí AJ, Korreman SS, Vestergaard A. Desafíos físicos en la radioterapia guiada por imágenes por resonancia magnética. Acta Oncol. 2020. https://pubmed.ncbi.nlm.nih.gov/32370788/

  2. Winkel D, Bol GH, Kroon PS, van Asselen B, Hackett SS, Werensteijn-Honingh AM, et al. Radioterapia adaptativa: implementación clínica del concepto Elekta Unity MR-linac en un consorcio MR-linac. Radiother Oncol. 2019. https://pubmed.ncbi.nlm.nih.gov/32432405/

  3. Cusumano D, Boldrini L, Dhont J, Fiorino C, Green O, Grimwood A, et al. Inteligencia artificial en radioterapia guiada por resonancia magnética: Consideraciones médicas y físicas sobre el estado del arte y perspectivas de futuro. Phys Med. 2022. https://pubmed.ncbi.nlm.nih.gov/35235188/

  4. Raaijmakers AJE, Raaymakers BW, Lagendijk JJW. Integración de un escáner de resonancia magnética con un acelerador de radioterapia de 6 MV: aumento de dosis en las interfaces tejido-aire en un campo magnético lateral debido al retorno de electrones. Phys Med Biol. 2005;50(7):1363–1376.

  5. Raaijmakers AJE, Raaymakers BW, van der Meer S, Lagendijk JJW. Integración de un escáner de resonancia magnética con un acelerador de radioterapia de 6 MV: impacto de la orientación de la superficie en la dosis de entrada y salida debido al campo magnético transversal. Phys Med Biol. 2007;52(4):929–939.

  6. Kirkby C, Stanescu T, Rathee S, Carlone M, Murray B, Fallone BG. Dosimetría del paciente para sistemas híbridos de resonancia magnética-radioterapia. Med Phys. 2008;35(3):1019–1027.

  7. Keyvanloo A, Burke B, Warkentin B, Tadic T, Rathee S, Kirkby C, et al. Dosis cutánea en sistemas MRI-linac longitudinales y transversales utilizando Monte Carlo y modelos de campo de MRI 3D realistas. Med Phys. 2012;39(10):6509–6521.

  8. O’Brien DJ, Roberts DA, Ibbott GS, Sawakuchi GO. Dosimetría de referencia en campos magnéticos: formalismo y medidas. Med Phys. 2016;43(8):4915–4927.