El principio de la radiografía de neutrones (nR) se basa en la atenuación de los neutrones a través de la materia que atraviesan. A diferencia de los rayos X que son dispersados por la nube de electrones de un átomo, los neutrones pueden ser absorbidos o dispersados por su núcleo. Los neutrones son sensibles a los elementos ligeros, como el hidrógeno (H), y en consecuencia pueden utilizarse para aplicaciones biológicas radiográficas como tejidos animales 1,2,3,4,5,6,7 o humanos 8,9 y sistemas subterráneos de suelo/raíces 10,11,12,13,14 ,15. La imagen de neutrones es una técnica complementaria a la imagen de rayos X, que es capaz de detectar elementos pesados16,17,18. El nR basado en la atenuación se rige por los coeficientes de atenuación lineal de los materiales dentro de la muestra y por el espesor de la muestra, como se describe en la ley de Beer-Lambert, que establece que el haz transmitido es directamente proporcional a la cantidad de material y la longitud del camino a través del material. Por lo tanto, la transmitancia, T, se puede calcular como:
(1)
donde I 0 e I son, respectivamente, las intensidades del haz incidente y transmitido; μ y x son el coeficiente de atenuación lineal y el espesor de una muestra homogénea, respectivamente. El coeficiente de atenuación μ viene dado por:
(2)
donde σ es la sección transversal de atenuación de neutrones de la muestra (tanto dispersión como absorción), ρ es su densidad, NA es el número de Avogadro y M es su masa molar.
El contraste en la radiografía de muestras biológicas que utilizan neutrones de baja energía (es decir, energías inferiores a 0,5 eV) se debe principalmente a un cambio en la densidad de H (para un espesor de muestra fijo). Esto se debe a la probabilidad de interacción de un neutrón con el núcleo H, que es mayor que con otros núcleos presentes en muestras biológicas, y al hecho de que la densidad del átomo de H es primordial ya que es el átomo más abundante en muestras biológicas.
Desde sus primeras etapas, nR y la tomografía computarizada de neutrones (nCT) se han utilizado ampliamente para aplicaciones de materiales e ingeniería 19,20,21,22,23. Los primeros experimentos de demostración de sensibilidad neutrónica al H en muestras biológicas comenzaron a mediados de la década de 195024 con las mediciones de especímenes de plantas. El trabajo continuó durante la década de 1960 con, por ejemplo, la radiografía de un tórax humano25 o ratas26, en la que se exploró el uso de agentes de contraste, como el óxido de gadolinio (Gd2O3). Además, se planteó la hipótesis de que el contraste en el tejido tumoral humano frente al tejido normal se debía a un aumento local en el contenido de H. Durante estos ensayos iniciales, se concluyó que el aumento del flujo de neutrones y la resolución espacial mejoraría la calidad de nR y probablemente aumentaría su popularidad como técnica complementaria para aplicaciones industriales o biomédicas. Los estudios más recientes comprenden mediciones de nR y nCT realizadas en muestras de tejido canceroso1 y secciones de órganos animales 2,3,27 para aplicaciones biomédicas y forenses.
Ubicado en el Laboratorio Nacional de Oak Ridge, Oak Ridge, TN, el Reactor de Isótopos de Alto Flujo (HFIR) es una poderosa fuente de neutrones que produce neutrones por reacción de fisión. Estos neutrones tienen energías del orden de 2 MeV y son "enfriados" en la piscina del reactor por reacciones cinéticas con agua pesada para alcanzar energías del orden de 100-300 eV. La optimización de un experimento de neutrones, ya sea dispersión o imagen, comienza con la comprensión de la fuente de neutrones y las propiedades de la línea de luz, como la intensidad del haz, la distribución de energía y el efecto del fondo (neutrones rápidos, neutrones retardados, rayos gamma). En la sala de guía fría HFIR donde se encuentra la línea de luz de imágenes, los neutrones se "enfrían" aún más por interacciones cinéticas con un moderador H líquido. Luego se transportan en un sistema de guía curvo lejos de la línea de visión de la fuente, eliminando así los neutrones rápidos y la contaminación gamma. Como se ilustra en la Figura 1, la línea de haz de imágenes de neutronesCG-1D 28,29 se coloca en una guía fría, lo que implica que el rango de energía de neutrones varía de unos pocos meV a unas pocas decenas de eV (en este caso, la longitud de onda de neutrones utilizable correspondiente varía de 0.8 a 10 Å) con un flujo en el rango de 107 n / (cm2∙s) en la posición de la muestra. Un sistema motorizado de apertura/difusor define la geometría estenopeica del instrumento de imagen. Los neutrones viajan una distancia de 6,59 m en un tubo de vuelo lleno de helio (He) con ventanas de aluminio (Al) en cada extremo. Los tubos de vuelo se utilizan para transportar neutrones al tiempo que limitan la dispersión del aire de tal manera que la pérdida en la intensidad del haz sea mínima. Para las mediciones descritas en este manuscrito, el difusor está hecho de un nanopolvo de óxido de aluminio de 50 nm (Al2O3) de 1 mm de espesor encerrado en un recipiente de Al. El difusor reduce los artefactos del haz provenientes de la guía de neutrones (que se magnifican por la geometría estenopeica de una línea de luz de imagen), de lo contrario, las fluctuaciones de intensidad horizontales y verticales agudas son visibles en la radiografía y la normalización de los datos se vuelve un desafío. Para los experimentos ilustrados aquí, los neutrones se convierten en luz utilizando un fluoruro de litio-6 de 25 μm de espesor / fósforo de sulfuro de zinc (6LiF / ZnS: Ag).
La optimización de la colimación depende de la posición de la muestra al detector, la resolución espacial requerida y el tiempo de adquisición. Cuando la muestra se encuentra a pocos cm del centelleador, las colimaciones altas (L / D por encima de 800, donde L es la distancia desde la apertura del agujero de alfiler del diámetro, D y el detector) producen una mejor resolución espacial a costa del flujo de neutrones. La baja colimación (L/D por debajo de 800) es preferible para estudios dinámicos in situ cuando la resolución temporal prevalece sobre la resolución espacial. Para las mediciones descritas en este manuscrito, L/D y resolución espacial fueron aproximadamente 355 y 75 μm, respectivamente. La resolución temporal varió en función de la relación señal-ruido (SNR). La muestra se colocó lo más cerca posible del centelleador para reducir la distorsión geométrica, como el desenfoque. Las etapas de traslación y rotación están disponibles para colocar la muestra cerca de los detectores y realizar la tomografía computarizada (TC). CG-1D ofrece tres tipos de detectores: un dispositivo de carga acoplada (CCD) con 2048 píxeles x 2048 píxeles con un paso de píxeles de 13,5 μm, un detector científico complementario de semiconductores de óxido metálico (sCMOS) con 2560 píxeles x 2160 píxeles con un paso de píxeles de 6,5 μm, y un detector de placa de microcanal (MCP)30,31 con 512 píxeles x 512 píxeles con un tamaño de píxel de 55 μm. Los neutrones dispersos se absorben con goma de boro de ~ 5 mm de espesor para proteger el chip detector de ver neutrones. Esta absorción genera rayos gamma que pueden ser detenidos por el plomo (Pb) colocado entre la goma de boro y el detector. Cada detector está optimizado para un campo de visión (FOV) diferente, así como para resoluciones espaciales y temporales. Para las mediciones del fémur de rata y el pulmón del ratón, se utilizó el detector CCD por su gran capacidad FOV (~ 7 cm x 7 cm) y una resolución espacial razonable de aproximadamente 75 μm. La nCT del sistema radicular/suelo de la planta se realizó con el sCMOS, ya que el objetivo era adquirir nCT lo más rápido posible a costa del campo de visión (que se limitó a ~ 5 cm x 4,2 cm); Por lo tanto, la resolución espacial evidentemente sufrió. En estos detectores, los neutrones se convierten en luz o en una partícula alfa para fines de detección. La rotación de la muestra alrededor de su eje vertical y la adquisición de radiografías en ángulos de rotación consecutivos permite la adquisición de nCT. El modelo de renderizado volumétrico tridimensional de la muestra bajo investigación se obtiene utilizando el portátil interno de retroproyección filtrada (FBP) Jupyter basado en python iMARS3D, pyMBIR o un software comercial, todo lo cual se describe a continuación.
Finalmente, los neutrones que no han interactuado con la muestra o el detector se recogen en una posición de parada del haz aproximadamente 1 m aguas abajo del sistema detector para minimizar el ruido de fondo. El tope de viga CG-1D tiene 0,75 m de ancho, 0,5 m de alto y 35 mm de grosor y está hecho de B4C en epoxi. El tope del haz está reforzado con 10 mm de carbonato de litio enriquecido al 95% (6 Li2CO3) en un epoxi resistente al fuego donde el haz de neutrones golpea, con una cavidad revestida con 6Li, plomo (Pb) y acero diseñado para contener la alta tasa de rayos gamma secundarios. El tope de la viga está directamente unido a la pared de blindaje de acero de la línea de luz. Una fotografía de la línea de luz CG-1D se da en la Figura 2.
Se utilizaron tres programas de reconstrucción para reconstruir los tres datos experimentales en 3D, respectivamente. La reconstrucción de la muestra de pulmón de ratón se realizó utilizando Octopus32, un software de reconstrucción comercial que utiliza FBP. El software Octopus se encuentra en una PC servidor y se puede utilizar para reconstruir los datos recopilados en la línea de luz. Un software de reconstrucción, llamado iMARS3D, está disponible en CG-1D. Se basa en el código fuente abierto TomoPY33 con características añadidas como corrección automática de inclinación, filtros de post-procesamiento, etc. iMARS3D incluye el preprocesamiento de los datos (resta del fondo y el ruido), el recorte, el filtrado de medianas (para corregir los rayos gamma y los píxeles muertos), la corrección automatizada de la fluctuación de la intensidad del haz y la corrección de la inclinación de la muestra. Una vez que se crean los sinogramas, el procesamiento adicional de datos, como la eliminación de artefactos de anillo y el suavizado, son una opción. Los diferentes pasos de la reconstrucción se guardan en el servidor de análisis (y luego se mueven en la carpeta compartida de la propuesta), mientras que los segmentos 2D finales se almacenan inmediatamente en la carpeta compartida de la propuesta. El fémur de la rata fue reconstruido usando iMARS3D. La muestra de raíz / suelo de la planta se procesó previamente mediante el filtrado mediano de los datos utilizando TomoPY seguido de la corrección del eje de inclinación utilizando la biblioteca SciPy de Python. La reconstrucción se llevó a cabo utilizando un paquete de Python desarrollado internamente denominado pyMBIR (construido con kernels de la caja de herramientas ASTRA34) que implementa un conjunto de algoritmos tomográficos desde la línea de base FBP hasta técnicas avanzadas de reconstrucción iterativa basadas en modelos35 que pueden obtener reconstrucciones de alta calidad a partir de conjuntos de datos de neutrones extremadamente dispersos y ruidosos. Todos los volúmenes renderizados basados en las herramientas de reconstrucción mencionadas anteriormente se representan en contraste de atenuación. Toda la visualización se realizó utilizando el paquete de software comercial de visualización, segmentación y análisis de datos AMIRA36.
Este manuscrito tiene como objetivo demostrar el procedimiento de uso de imágenes de neutrones (nR y nCT) en la línea de luz HFIR CG-1D. Este estudio también ilustra las capacidades actuales de nR y nCT de vanguardia para muestras biológicas, específicamente un pulmón de ratón, un hueso de rata y sistemas de raíces / suelo de plantas. El pulmón de ratón fue elegido para ilustrar la complementariedad de los neutrones para medir el tejido pulmonar, mientras que los rayos X son principalmente sensibles a los huesos. La muestra ósea, un fémur de rata, tenía un implante de titanio (Ti), lo que ilustra el contraste entre el hueso y el metal, y la oportunidad de ver la interfaz hueso/metal (que es difícil de medir con rayos X ya que los metales los atenúan fuertemente4). Finalmente, el sistema de agua planta-raíz ilustra la capacidad tridimensional (3D) de nCT para medir los sistemas de raíces / suelo in situ. Además, muestra las ventajas / desventajas de usar nR para muestras biológicas. Evidentemente, este método puede usarse con seguridad para medir la dinámica del agua en un sistema planta-raíz, pero no puede considerarse como una técnica de imagen animal o humana viva debido a los riesgos asociados con la exposición a la radiación, limitando así los estudios a ratones (muertos) o mediciones similares a la patología en las que, por ejemplo, una muestra de tejido se reseca de un paciente (animal o humano) y se prepara por fijación antes de medirse en un haz de neutrones.