Das Prinzip der Neutronenradiographie (nR) beruht auf der Abschwächung von Neutronen durch die Materie, die sie durchqueren. Im Gegensatz zu Röntgenstrahlen, die von der Elektronenwolke eines Atoms gestreut werden, können Neutronen von seinem Atomkern absorbiert oder gestreut werden. Neutronen sind empfindlich gegenüber leichten Elementen wie Wasserstoff (H) und können daher zur Röntgenaufnahme biologischer Anwendungen wie Tier 1,2,3,4,5,6,7 oder menschlichem Gewebe 8,9 und unterirdischen Boden-/Wurzelsystemen 10,11,12,13,14 verwendet werden ,15. Die Neutronenbildgebung ist eine komplementäre Technik zur Röntgenbildgebung, die in der Lage ist, schwere Elemente zu detektieren16,17,18. Die dämpfungsbasierte nR wird durch die linearen Dämpfungskoeffizienten der Materialien innerhalb der Probe und durch die Dicke der Probe bestimmt, wie sie durch das Beer-Lambert-Gesetz beschrieben wird, das besagt, dass der transmittierte Strahl direkt proportional zur Materialmenge und der Weglänge durch das Material ist. Somit kann der Transmissionsgrad T wie folgt berechnet werden:
(1)
wobei I 0 bzw. I die einfallenden bzw. transmittierten Strahlintensitäten sind; μ und x sind der lineare Abschwächungskoeffizient bzw. die Dicke einer homogenen Probe. Der Dämpfungskoeffizient μ ist gegeben durch:
(2)
wobei σ der Neutronendämpfungsquerschnitt der Probe (sowohl Streuung als auch Absorption), ρ ihre Dichte, NA die Avogadro-Zahl und M ihre Molmasse ist.
Der Kontrast in der Radiographie biologischer Proben mit niederenergetischen Neutronen (d. h. Energien unter 0,5 eV) ist hauptsächlich auf eine Änderung der Dichte von H (für eine feste Probendicke) zurückzuführen. Dies liegt an der Wahrscheinlichkeit der Wechselwirkung eines Neutrons mit dem H-Kern, die größer ist als bei anderen Kernen, die in biologischen Proben vorhanden sind, und an der Tatsache, dass die Dichte des H-Atoms von größter Bedeutung ist, da es das am häufigsten vorkommende Atom in biologischen Proben ist.
Seit ihren Anfängen werden die nR- und Neutronen-Computertomographie (nCT) in großem Umfang für Materialien und technische Anwendungen eingesetzt 19,20,21,22,23. Die ersten Demonstrationsexperimente zur Neutronenempfindlichkeit gegenüber H in biologischen Proben begannen Mitte der 1950er Jahre24 mit den Messungen von Pflanzenproben. Die Arbeit wurde bis in die 1960er Jahre fortgesetzt, z. B. mit der Röntgenaufnahme einer menschlichen Brust25 oder einer Ratte26, in der der Einsatz von Kontrastmitteln wie Gadoliniumoxid (Gd2O3) erforscht wurde. Darüber hinaus wurde die Hypothese aufgestellt, dass der Kontrast zwischen menschlichem Tumorgewebe und normalem Gewebe auf einen lokalen Anstieg des H-Gehalts zurückzuführen ist. Während dieser ersten Versuche kam man zu dem Schluss, dass ein erhöhter Neutronenfluss und eine höhere räumliche Auflösung die Qualität von nR verbessern und wahrscheinlich seine Popularität als ergänzende Technik für industrielle oder biomedizinische Anwendungen erhöhen würden. Die jüngsten Studien umfassen nR- und nCT-Messungen an Krebsgewebeproben1 und Schnitten tierischer Organe 2,3,27 für biomedizinische und forensische Anwendungen.
Der High Flux Isotope Reactor (HFIR) befindet sich am Oak Ridge National Laboratory, Oak Ridge, Tennessee, und ist eine leistungsstarke Neutronenquelle, die Neutronen durch Spaltreaktion erzeugt. Diese Neutronen haben Energien in der Größenordnung von 2 MeV und werden im Reaktorbecken durch kinetische Reaktionen mit schwerem Wasser "gekühlt", um Energien in der Größenordnung von 100-300 eV zu erreichen. Die Optimierung eines Neutronenexperiments, ob Streuung oder Bildgebung, beginnt mit dem Verständnis der Eigenschaften der Neutronenquelle und der Strahllinie wie ihrer Strahlintensität, Energieverteilung und der Wirkung des Hintergrunds (schnelle Neutronen, verzögerte Neutronen, Gammastrahlen). In der HFIR-Kälteleithalle, in der sich die abbildende Beamline befindet, werden die Neutronen durch kinetische Wechselwirkungen mit einem flüssigen H-Moderator weiter "gekühlt". Sie werden dann in einem gekrümmten Leitsystem von der Sichtlinie der Quelle weg transportiert, wodurch schnelle Neutronen und Gamma-Verschmutzung eliminiert werden. Wie in Abbildung 1 dargestellt, ist die CG-1D-Neutronenabbildungsstrahllinie28,29 auf einer kalten Führung platziert, was bedeutet, dass der Neutronenenergiebereich von einigen meV bis zu einigen Dutzend eV variiert (in diesem Fall reicht die entsprechende nutzbare Neutronenwellenlänge von 0,8 bis 10 Å) mit einem Fluss im Bereich von 107 n/(cm2∙s) an der Probenposition. Ein motorisiertes Blenden-/Diffusorsystem definiert die Lochblendengeometrie des Bildgebungsinstruments. Neutronen legen eine Strecke von 6,59 m in einem mit Helium (He) gefüllten Flugrohr mit Aluminiumfenstern (Al) an jedem Ende zurück. Flugröhren werden verwendet, um Neutronen zu transportieren und gleichzeitig die Luftstreuung so zu begrenzen, dass der Verlust der Strahlintensität minimal ist. Für die in diesem Manuskript beschriebenen Messungen besteht der Diffusor aus einem 1 mm dicken 50 nm Aluminiumoxid (Al2O3) Nanopulver, das in einem Al-Behälter eingeschlossen ist. Der Diffusor reduziert die vom Neutronenleiter ausgehenden Strahlartefakte (die durch die Lochblendengeometrie einer abbildenden Strahllinie vergrößert werden), da sonst scharfe horizontale und vertikale Intensitätsschwankungen im Röntgenbild sichtbar sind und die Normalisierung der Daten schwierig wird. Für die hier dargestellten Experimente werden Neutronen mit einem 25 μm dicken Lithium-6-Fluorid/Zinksulfid-Phosphor (6LiF/ZnS:Ag) in Licht umgewandelt.
Die Kollimationsoptimierung hängt von der Position der Probe zum Detektor, der erforderlichen räumlichen Auflösung und der Erfassungszeit ab. Wenn die Probe einige Zentimeter vom Szintillator entfernt sitzt, führen hohe Kollimationen (L/D über 800, wobei L der Abstand von der Lochblende des Durchmessers D und des Detektors ist) zu einer besseren räumlichen Auflösung auf Kosten des Neutronenflusses. Eine niedrige Kollimation (L/D unter 800) ist für dynamische In-situ-Studien vorzuziehen, bei denen die Zeitauflösung die räumliche Auflösung überwiegt. Für die in diesem Manuskript beschriebenen Messungen betrug L/D und räumliche Auflösung etwa 355 bzw. 75 μm. Die zeitliche Auflösung variierte je nach Signal-Rausch-Verhältnis (SNR). Die Probe wurde so nah wie möglich am Szintillator positioniert, um geometrische Verzerrungen wie Unschärfen zu reduzieren. Translations- und Rotationstische stehen zur Verfügung, um die Probe in die Nähe der Detektoren zu bringen und eine Computertomographie (CT) durchzuführen. CG-1D bietet drei Arten von Detektoren: ein ladungsgekoppeltes Gerät (CCD) mit 2048 Pixeln x 2048 Pixeln mit einem Pixelabstand von 13,5 μm, einen wissenschaftlichen komplementären Metalloxid-Halbleiter (sCMOS)-Detektor mit 2560 Pixeln x 2160 Pixeln mit einem Pixelabstand von 6,5 μm und einen Micro-Channel-Plate-Detektor (MCP)30,31 mit 512 Pixeln x 512 Pixeln mit einer Pixelgröße von 55 μm. Gestreute Neutronen werden mit ~5 mm dickem Borkautschuk absorbiert, um den Detektorchip vor Neutronen zu schützen. Diese Absorption erzeugt Gammastrahlen, die durch Blei (Pb) gestoppt werden können, das zwischen dem Borkautschuk und dem Detektor platziert wird. Jeder Detektor ist für ein anderes Sichtfeld (FOV) sowie für räumliche und zeitliche Auflösungen optimiert. Für die Messungen des Oberschenkelknochens der Ratte und der Lunge der Maus wurde der CCD-Detektor aufgrund seines großen Sichtfelds (~ 7 cm x 7 cm) und seiner angemessenen räumlichen Auflösung von ca. 75 μm verwendet. Die nCT des Wurzel-Boden-Systems der Pflanze wurde mit dem sCMOS durchgeführt, da das Ziel darin bestand, nCTs so schnell wie möglich auf Kosten des FOV (das auf ~ 5 cm x 4,2 cm begrenzt war) zu erhalten. Die räumliche Auflösung hat also offensichtlich gelitten. In diesen Detektoren werden Neutronen zu Detektionszwecken entweder in Licht oder ein Alpha-Teilchen umgewandelt. Die Drehung der Probe um ihre vertikale Achse und die Aufnahme von Röntgenaufnahmen bei aufeinanderfolgenden Drehwinkeln ermöglicht die Aufnahme von nCT. Das 3-dimensionale volumetrisch gerenderte Modell der untersuchten Probe wird durch die Verwendung des hauseigenen iMARS3D-Python-basierten Jupyter-FBP-Notebooks (Filtered-Back-Projection), pyMBIR oder einer kommerziellen Software erhalten, die im Folgenden beschrieben wird.
Schließlich werden Neutronen, die nicht mit der Probe oder dem Detektor interagiert haben, in einer Strahlstoppposition ca. 1 m stromabwärts des Detektorsystems gesammelt, um das Hintergrundrauschen zu minimieren. Der Strahlanschlag CG-1D ist 0,75 m breit, 0,5 m hoch und 35 mm dick und besteht aus B4C in Epoxidharz. Der Strahlanschlag ist mit 10 mm 95% angereichertem Lithiumcarbonat (6 Li2CO3) in einem feuerfesten Epoxidharz verstärkt, wo der Neutronenstrahl auftrifft, mit einem Hohlraum, der mit 6Li, Blei (Pb) und Stahl ausgekleidet ist, um die hohe Rate der sekundären Gammastrahlen einzudämmen. Der Trägeranschlag wird direkt an der Stahlabschirmwand des Strahlrohrs befestigt. Ein Foto des CG-1D-Strahlrohrs ist in Abbildung 2 dargestellt.
Es wurden jeweils drei Rekonstruktionsprogramme verwendet, um die drei experimentellen Daten in 3D zu rekonstruieren. Die Rekonstruktion der Lungenprobe der Maus wurde mit Octopus32 durchgeführt, einer kommerziellen Rekonstruktionssoftware, die FBP verwendet. Die Octopus-Software befindet sich auf einem Server-PC und kann verwendet werden, um die an der Beamline gesammelten Daten zu rekonstruieren. Eine Rekonstruktionssoftware namens iMARS3D ist unter CG-1D verfügbar. Es basiert auf dem Open-Source-Code TomoPY33 mit zusätzlichen Funktionen wie automatischer Neigungskorrektur, Nachbearbeitungsfiltern usw. iMARS3D umfasst die Vorverarbeitung der Daten (Subtraktion des Hintergrunds und des Rauschens), das Zuschneiden, die Medianfilterung (zur Korrektur von Gammaschlägen und toten Pixeln), die automatische Korrektur von Schwankungen der Strahlintensität und die Korrektur der Probenneigung. Sobald Sinogramme erstellt sind, ist eine weitere Datenverarbeitung wie das Entfernen von Ringartefakten und das Glätten möglich. Die verschiedenen Schritte der Rekonstruktion werden auf dem Analyseserver gespeichert (und später in den freigegebenen Vorschlagsordner verschoben), während die endgültigen 2D-Schichten sofort im freigegebenen Vorschlagsordner gespeichert werden. Der Oberschenkelknochen der Ratte wurde mit iMARS3D rekonstruiert. Die Pflanzenwurzel-/Bodenprobe wurde durch Medianfilterung der Daten mit TomoPY vorverarbeitet, gefolgt von einer Neigungsachsenkorrektur mit der SciPy-Bibliothek von Python. Die Rekonstruktion wurde unter Verwendung eines intern entwickelten Python-Pakets mit dem Namen pyMBIR (erstellt unter Verwendung von Kernen aus der ASTRA-Toolbox34) durchgeführt, das eine Reihe von tomographischen Algorithmen vom Basis-FBP bis hin zu fortgeschrittenen modellbasierten iterativen Rekonstruktionstechniken35 implementiert, die qualitativ hochwertige Rekonstruktionen aus extrem spärlichen und verrauschten Neutronendatensätzen erhalten können. Alle gerenderten Volumina, die auf den oben genannten Rekonstruktionswerkzeugen basieren, werden im Dämpfungskontrast dargestellt. Die gesamte Visualisierung erfolgte mit dem kommerziellen Visualisierungs-, Segmentierungs- und Datenanalyse-Softwarepaket AMIRA36.
Ziel dieses Manuskripts ist es, das Verfahren der Neutronenbildgebung (nR und nCT) an der HFIR CG-1D-Beamline zu demonstrieren. Diese Studie veranschaulicht auch den aktuellen Stand der Technik der nR- und nCT-Fähigkeiten für biologische Proben, insbesondere für eine Mauslunge, einen Rattenknochen und pflanzliche Wurzel-/Bodensysteme. Die Mauslunge wurde ausgewählt, um die Komplementarität von Neutronen zur Messung des Lungengewebes zu veranschaulichen, während Röntgenstrahlen hauptsächlich empfindlich für Knochen sind. Die Knochenprobe, ein Oberschenkelknochen einer Ratte, war mit einem Titan (Ti)-Implantat versehen, was den Kontrast zwischen Knochen und Metall veranschaulichte und die Möglichkeit bot, die Grenzfläche zwischen Knochen und Metall zu sehen (die mit Röntgenstrahlen schwer zu messen ist, da Metalle sie stark abschwächen4). Schließlich veranschaulicht das Pflanzen-Wurzel-Wassersystem die dreidimensionale (3D) Fähigkeit von nCT, Wurzel-/Bodensysteme in situ zu messen. Darüber hinaus werden die Vor- und Nachteile der Verwendung von nR für biologische Proben aufgezeigt. Offensichtlich kann diese Methode sicher verwendet werden, um die Wasserdynamik in einem Pflanzen-Wurzel-System zu messen, kann aber aufgrund der mit der Strahlenbelastung verbundenen Risiken nicht als lebendes Tier- oder menschliches Bildgebungsverfahren betrachtet werden, so dass Studien entweder auf (tote) Mäuse oder pathologieähnliche Messungen beschränkt sind, bei denen beispielsweise eine Gewebeprobe von einem Patienten (Tier oder Mensch) entnommen und durch Fixierung präpariert wird, bevor sie in einem Neutronenstrahl gemessen wird.