[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

Search Results (166)

Search Parameters:
Keywords = Geant4

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
24 pages, 6808 KiB  
Article
Single-Particle Radiation Sensitivity of Ultrawide-Bandgap Semiconductors to Terrestrial Atmospheric Neutrons
by Daniela Munteanu and Jean-Luc Autran
Crystals 2025, 15(2), 186; https://doi.org/10.3390/cryst15020186 - 15 Feb 2025
Viewed by 263
Abstract
Semiconductors characterized by ultrawide bandgaps (UWBGs), exceeding the SiC bandgap of 3.2 eV and the GaN bandgap of 3.4 eV, are currently under focus for applications in high-power and radio-frequency (RF) electronics, as well as in deep-ultraviolet optoelectronics and extreme environmental conditions. These [...] Read more.
Semiconductors characterized by ultrawide bandgaps (UWBGs), exceeding the SiC bandgap of 3.2 eV and the GaN bandgap of 3.4 eV, are currently under focus for applications in high-power and radio-frequency (RF) electronics, as well as in deep-ultraviolet optoelectronics and extreme environmental conditions. These semiconductors offer numerous advantages, such as a high breakdown field, exceptional thermal stability, and minimized power losses. This study used numerical simulation to investigate, at the material level, the single-particle radiation response of various UWBG semiconductors, such as aluminum gallium nitride alloys (AlxGa1−xN), diamond, and β-phase gallium oxide (β-Ga2O3), when exposed to ground-level neutrons. Through comprehensive Geant4 simulations covering the entire spectrum of atmospheric neutrons at sea level, this study provides an accurate comparison of the neutron radiation responses of these UWBG semiconductors focusing on the interaction processes, the number and nature of secondary ionizing products, their energy distributions, and the production of electron–hole pairs at the origin of single-event effects (SEEs) in microelectronics devices. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Density and bandgap of the ternary alloy Al<sub>x</sub>Ga<sub>1−x</sub>N as a function of the Al content <math display="inline"><semantics> <mrow> <mi>x</mi> </mrow> </semantics></math>.</p>
Full article ">Figure 2
<p>Spectrum of the atmospheric neutron (lethargic representation) at sea level recorded on the roof of the IBM Watson Research Center main building [<a href="#B48-crystals-15-00186" class="html-bibr">48</a>] in New York City. Experimental data courtesy of Paul Goldhagen (U.S. Department of Homeland Security). Percentages of the total flux of neutrons related to each domain of the spectrum are also indicated.</p>
Full article ">Figure 3
<p>Total number of interactions generated by Geant4 simulations of Al<sub>x</sub>Ga<sub>1−x</sub>N, diamond, and β-Ga<sub>2</sub>O<sub>3</sub> bulk targets under the full spectrum of atmospheric neutron shown in <a href="#crystals-15-00186-f002" class="html-fig">Figure 2</a> for an equivalent exposure time of 2.5 × 10<sup>7</sup> h.</p>
Full article ">Figure 4
<p>Number of neutron interactions in Al<sub>x</sub>Ga<sub>1−x</sub>N, diamond, and β-Ga<sub>2</sub>O<sub>3</sub> bulk targets subjected to part I, II, or III of the spectrum of atmospheric neutrons (<a href="#crystals-15-00186-f002" class="html-fig">Figure 2</a>) over 2.5 × 10<sup>7</sup> h.</p>
Full article ">Figure 5
<p>Neutron cross-section versus the neutron energy for diamond obtained from the TENDL open nuclear data library.</p>
Full article ">Figure 6
<p>Number of events of elastic and inelastic scattering and nuclear reactions in Al<sub>x</sub>Ga<sub>1−x</sub>N, diamond and β-Ga<sub>2</sub>O<sub>3</sub> bulk targets under the full spectrum of atmospheric neutrons shown in <a href="#crystals-15-00186-f002" class="html-fig">Figure 2</a> for an equivalent exposure time of 2.5 × 10<sup>7</sup> h.</p>
Full article ">Figure 7
<p>Number of secondary products in Al<sub>x</sub>Ga<sub>1−x</sub>N, diamond, and β-Ga<sub>2</sub>O<sub>3</sub> bulk targets subjected to part I, II, or III of the spectrum of atmospheric neutrons (<a href="#crystals-15-00186-f002" class="html-fig">Figure 2</a>) over 2.5 × 10<sup>7</sup> h.</p>
Full article ">Figure 8
<p>Number of events versus the shower multiplicity for Al<sub>x</sub>Ga<sub>1−x</sub>N, diamond, and β-Ga<sub>2</sub>O<sub>3</sub> bulk targets subjected to the full spectrum of atmospheric neutrons (<a href="#crystals-15-00186-f002" class="html-fig">Figure 2</a>) over 2.5 × 10<sup>7</sup> h.</p>
Full article ">Figure 9
<p>Number of secondary products versus the atomic number Z for GaN, Al<sub>0.6</sub>Ga<sub>0.4</sub>N, AlN, diamond, and β-Ga<sub>2</sub>O<sub>3</sub> bulk under the full spectrum of atmospheric neutron shown in <a href="#crystals-15-00186-f002" class="html-fig">Figure 2</a> over 2.5 × 10<sup>7</sup> h.</p>
Full article ">Figure 10
<p>Energy histograms (500 bins) of the secondary products generated in the Al<sub>x</sub>Ga<sub>1−x</sub>N alloy targets (1 cm<sup>2</sup> × 20 μm) irradiated with neutrons of energies less than 1 eV (part I of the spectrum) over 2.5 × 10<sup>7</sup> h. (<b>a</b>) GaN; (<b>b</b>) Al<sub>0.6</sub>Ga<sub>0.4</sub>N; (<b>c</b>) AlN.</p>
Full article ">Figure 11
<p>Energy histograms (500 bins) of the secondary products generated in the diamond and β-Ga<sub>2</sub>O<sub>3</sub> targets (1 cm<sup>2</sup> × 20 μm) irradiated with neutrons of energies less than 1 eV (part I of the spectrum) over 2.5 × 10<sup>7</sup> h. (<b>a</b>) Diamond; (<b>b</b>) β-Ga<sub>2</sub>O<sub>3</sub>.</p>
Full article ">Figure 12
<p>Energy histograms (500 bins) of the secondary products generated in the Al<sub>x</sub>Ga<sub>1−x</sub>N alloy targets (1 cm<sup>2</sup> × 20 μm) irradiated with neutrons of energies between 1 eV and 1 MeV (part II of the spectrum) over 2.5 × 10<sup>7</sup> h. (<b>a</b>) GaN; (<b>b</b>) Al<sub>0.6</sub>Ga<sub>0.4</sub>N; (<b>c</b>) AlN.</p>
Full article ">Figure 13
<p>Energy histograms (500 bins) of the secondary products generated in the diamond and β-Ga<sub>2</sub>O<sub>3</sub> targets (1 cm<sup>2</sup> × 20 μm) irradiated with neutrons of energies between 1 eV and 1 MeV (part II of the spectrum) over 2.5 × 10<sup>7</sup> h. (<b>a</b>) Diamond; (<b>b</b>) β-Ga<sub>2</sub>O<sub>3</sub>.</p>
Full article ">Figure 14
<p>Energy histograms (500 bins) of the secondary products generated in the Al<sub>x</sub>Ga<sub>1−x</sub>N alloy targets (1 cm<sup>2</sup> × 20 μm) irradiated with neutrons of energies higher than 1 MeV (part III of the spectrum) over 2.5 × 10<sup>7</sup> h. (<b>a</b>) GaN; (<b>b</b>) Al<sub>0.6</sub>Ga<sub>0.4</sub>N; (<b>c</b>) AlN.</p>
Full article ">Figure 15
<p>Energy histograms (500 bins) of the secondary products generated in the diamond and β-Ga<sub>2</sub>O<sub>3</sub> targets (1 cm<sup>2</sup> × 20 μm) irradiated with neutrons of energies higher than 1 MeV (part III of the spectrum) over 2.5 × 10<sup>7</sup> h. (<b>a</b>) Diamond; (<b>b</b>) β-Ga<sub>2</sub>O<sub>3</sub>.</p>
Full article ">Figure 16
<p>Number of interactions capable of depositing at least 0.7 fC of charge in the GaN, Al<sub>0.6</sub>Ga<sub>0.4</sub>N, diamond, and β-Ga<sub>2</sub>O<sub>3</sub> bulk targets subjected to part I, part II, and part III of the spectrum of atmospheric neutrons (<a href="#crystals-15-00186-f002" class="html-fig">Figure 2</a>) for an equivalent exposure time of 2.5 × 10<sup>7</sup> h.</p>
Full article ">Figure 17
<p>Total number of interactions capable of depositing at least 0.7 fC of charge in the GaN, Al<sub>0.6</sub>Ga<sub>0.4</sub>N, diamond, and β-Ga<sub>2</sub>O<sub>3</sub> bulk targets under the full spectrum of atmospheric neutron shown in <a href="#crystals-15-00186-f002" class="html-fig">Figure 2</a> for an equivalent exposure time of 2.5 × 10<sup>7</sup> h.</p>
Full article ">
16 pages, 1408 KiB  
Article
Feasibility Study of a PET Detector with a Wavelength-Shifting Fiber Readout
by Anzori Sh. Georgadze
Instruments 2025, 9(1), 2; https://doi.org/10.3390/instruments9010002 - 5 Feb 2025
Viewed by 694
Abstract
We designed and evaluated the performance of a high-resolution large-area detector for positron emission tomography (PET) based on a crystal assembly readout using wavelength-shifting (WLS) fibers, offering a cost-effective alternative to the direct readout of monolithic crystals with photodetectors. The considered detector geometries [...] Read more.
We designed and evaluated the performance of a high-resolution large-area detector for positron emission tomography (PET) based on a crystal assembly readout using wavelength-shifting (WLS) fibers, offering a cost-effective alternative to the direct readout of monolithic crystals with photodetectors. The considered detector geometries were made up of 4 × 4 assemblies of LuY2SiO5:Ce (LYSO) crystal scintillators, each with surface area of 50 × 50 mm2 and thickness of 7 or 15 mm, which were optically coupled together using optical adhesive. The crystal assembly was coupled with square cross-sections of orthogonal wavelength-shifting (WLS) fibers placed on the top and bottom of the assembly. To evaluate the characteristics of the novel detector, we used GEANT4 to perform optical photon transport in the crystal assembly and WLS fibers. The simulation results show that best position resolution achieved was 1.6 ± 0.4 mm full width at half maximum (FWHM) and 4.2 ± 0.6 mm full width at tenth maximum (FWTM) for the crystal thickness of 7 mm and 1.7 ± 0.4 mm FWHM and 6.0 ± 0.6 mm FWTM for the crystal thickness of 15 mm. Compared with a direct photosensor readout, WLS fibers can drastically reduce the number of photosensors required while covering a larger sensitive detection area. In the proposed detector design, 2N photodetectors are used to cover the same image area instead of N2 with a direct readout. This design allows for the development of a compact detector with an expanded effective field of view and reduced cost. Full article
(This article belongs to the Special Issue Medical Applications of Particle Physics, 2nd Edition)
Show Figures

Figure 1

Figure 1
<p>Schematic of detector module composed of a 4 × 4 assembly of LYSO crystals. Red boxes represent SiPMS, blue cubes represent the individual LYSO scintillators, and green rectangles represent the WLS fibers.</p>
Full article ">Figure 2
<p>Photon emission spectrum of LYSO scintillator [<a href="#B34-instruments-09-00002" class="html-bibr">34</a>] (blue color), absorption (red color), and emission (green color) spectra of BCF-91A [<a href="#B35-instruments-09-00002" class="html-bibr">35</a>], SiPM PDE (dark blue color), and absorption (magenta color) and emission (dark green color) spectra of BCF-92 [<a href="#B35-instruments-09-00002" class="html-bibr">35</a>].</p>
Full article ">Figure 3
<p>Example of one event with the secondary scintillation photons simulated in GEANT4 for a LYSO assembly. The light green lines are the tracks of optical photons, and the red boxes are SiPMs.</p>
Full article ">Figure 4
<p>(<b>a</b>,<b>b</b>) The light distribution profiles along the <math display="inline"><semantics> <mi mathvariant="italic">x</mi> </semantics></math> (<b>a</b>) and <math display="inline"><semantics> <mi mathvariant="italic">y</mi> </semantics></math> (<b>b</b>) directions in a typical event displaying photoelectron absorption-like behavior. The simulated data were fitted with a Gaussian function. The blue histogram is the simulated results, and the red histogram is the Gaussian fitting, while the green histogram is a high-resolution peak search function. (<b>c</b>) The <math display="inline"><semantics> <mi mathvariant="italic">x</mi> </semantics></math>–<math display="inline"><semantics> <mi mathvariant="italic">y</mi> </semantics></math> distribution image obtained by combining signals from the <math display="inline"><semantics> <mi mathvariant="italic">x</mi> </semantics></math> and <math display="inline"><semantics> <mi mathvariant="italic">y</mi> </semantics></math> WLS fibers. (<b>d</b>) The <math display="inline"><semantics> <mi mathvariant="italic">x</mi> </semantics></math>–<math display="inline"><semantics> <mi mathvariant="italic">y</mi> </semantics></math> distribution image in 3D.</p>
Full article ">Figure 5
<p>(<b>a</b>,<b>b</b>) The light distribution profiles along the <math display="inline"><semantics> <mi mathvariant="italic">x</mi> </semantics></math> (<b>a</b>) and <math display="inline"><semantics> <mi mathvariant="italic">y</mi> </semantics></math> (<b>b</b>) directions in a typical event displaying the Compton scattering event event. The simulated data were fitted with a Gaussian function. The blue histogram is the simulated results, and the red histogram is the Gaussian fitting, while the green histogram is a high-resolution peak search function. (<b>c</b>) The <math display="inline"><semantics> <mi mathvariant="italic">x</mi> </semantics></math>–<math display="inline"><semantics> <mi mathvariant="italic">y</mi> </semantics></math> distribution image obtained by combining signals from the <math display="inline"><semantics> <mi mathvariant="italic">x</mi> </semantics></math> and <math display="inline"><semantics> <mi mathvariant="italic">y</mi> </semantics></math> WLS fibers. (<b>d</b>) The <math display="inline"><semantics> <mi mathvariant="italic">x</mi> </semantics></math>–<math display="inline"><semantics> <mi mathvariant="italic">y</mi> </semantics></math> distribution image in 3D.</p>
Full article ">Figure 6
<p>Histogram of reconstructed interaction positions (blue line) of a 511 keV <math display="inline"><semantics> <mi>γ</mi> </semantics></math> ray beam positioned at (<math display="inline"><semantics> <mi mathvariant="italic">x</mi> </semantics></math>, <math display="inline"><semantics> <mi mathvariant="italic">y</mi> </semantics></math>) = (20 mm, 15 mm) for a crystal thickness of 7 mm (<b>a</b>,<b>b</b>) and for a crystal thickness of 15 mm (<b>c</b>,<b>d</b>). A Lorentzian fit to the distribution is also shown (red line).</p>
Full article ">Figure 7
<p>XY histograms of interaction position reconstruction of the 511 keV <math display="inline"><semantics> <mi>γ</mi> </semantics></math> ray beams positioned at (<math display="inline"><semantics> <mi mathvariant="italic">x</mi> </semantics></math>, <math display="inline"><semantics> <mi mathvariant="italic">y</mi> </semantics></math>) = (15 mm, 15 mm) and (<math display="inline"><semantics> <mi mathvariant="italic">x</mi> </semantics></math>, <math display="inline"><semantics> <mi mathvariant="italic">y</mi> </semantics></math>) = (20 mm, 15 mm), shown in both 3D (<b>top</b>) and 2D (<b>bottom</b>) representations (blue dots) for a crystal thickness of 7 mm: (<b>a</b>,<b>c</b>) no rejection and (<b>b</b>,<b>d</b>) with Compton scattering events rejection applied.</p>
Full article ">Figure 8
<p>XY histograms of interaction position reconstruction of the 511 keV <math display="inline"><semantics> <mi>γ</mi> </semantics></math> ray beams positioned at (<math display="inline"><semantics> <mi mathvariant="italic">x</mi> </semantics></math>, <math display="inline"><semantics> <mi mathvariant="italic">y</mi> </semantics></math>) = (15 mm, 15 mm) and (<math display="inline"><semantics> <mi mathvariant="italic">x</mi> </semantics></math>, <math display="inline"><semantics> <mi mathvariant="italic">y</mi> </semantics></math>) = (20 mm, 15 mm), shown in both 3D (<b>top</b>) and 2D (<b>bottom</b>) representations (blue dots) for a crystal thickness of 15 mm: (<b>a</b>,<b>c</b>) no rejection and (<b>b</b>,<b>d</b>) with Compton scattering events rejection applied.</p>
Full article ">Figure 9
<p>Response time spectrum of detected photons in a single event for 511 keV energy deposition in a LYSO crystal array readout using BCF-91A fibers (<b>a</b>) and BCF-92 fibers (<b>b</b>). (<b>c</b>) The average (n = 2000) response time spectrum of detected photons for 511 keV energy deposition in a LYSO crystal array, shown for the readout with BCF-91A fibers (blue histogram fitted with an exponential function in red) and the readout with BCF-92 fibers (green histogram fitted with an exponential function in black).</p>
Full article ">Figure 10
<p>(<b>a</b>) A 2D histogram of the positioning estimations using a ROOT-based algorithm, demonstrating interaction uniformly distributed across a grid of 20 × 20 points spaced 10 mm apart. (<b>b</b>) The <math display="inline"><semantics> <mi mathvariant="italic">z</mi> </semantics></math>-coordinate of the interaction position, simulated at six depths, regarded as the depth-of-interaction (DOI) resolution.</p>
Full article ">Figure 11
<p>Three dimensional plot of reconstructed (<math display="inline"><semantics> <mrow> <mi mathvariant="italic">x</mi> <mo>,</mo> <mi mathvariant="italic">y</mi> <mo>,</mo> <mi mathvariant="italic">z</mi> </mrow> </semantics></math>) coordinates for the 511 keV <math display="inline"><semantics> <mi>γ</mi> </semantics></math>-ray beam positioned at (<math display="inline"><semantics> <mi mathvariant="italic">x</mi> </semantics></math>, <math display="inline"><semantics> <mi mathvariant="italic">y</mi> </semantics></math>) = (30 mm, 30 mm) shown in with no rejection (<b>a</b>) and with Compton scattering event rejection applied (<b>b</b>). An energy spectrum of 511 keV <math display="inline"><semantics> <mi>γ</mi> </semantics></math>-rays for the modeled WLS-PET detector is shown together with the corresponding Gaussian fit through the corresponding full energy peak (red line) (<b>c</b>).</p>
Full article ">Figure 12
<p>The 3D (<b>a</b>) and 2D (<b>b</b>) histograms of the average number of detected photons, demonstrating uniformity of energy resolution across the detector’s sensitive area. In the histogram, darker yellow bins indicate fewer detected photons, while brighter yellow bins represent a higher number of detected photons.</p>
Full article ">Figure 13
<p>Schematic of a brain PET scanner (<b>a</b>) and a total-body PET scanner (<b>b</b>), both composed of modules from a 4 × 4 assembly of optically coupled LYSO crystals (gray). The interaction of annihilation <math display="inline"><semantics> <mi>γ</mi> </semantics></math>-rays with the opposite detection modules produces scintillation photons (green), which are shared within the LYSO assembly due to optical coupling. WLS fibers are not shown.</p>
Full article ">
10 pages, 1831 KiB  
Article
B-10-Based Macrostructured Cathode for Neutron Detectors
by Alexander G. Kolesnikov, Aleksey K. Kurilkin, Viktor I. Bodnarchuk, Alexander S. Ovodov and Marat R. Gafurov
Coatings 2025, 15(2), 168; https://doi.org/10.3390/coatings15020168 - 2 Feb 2025
Viewed by 611
Abstract
This paper focuses on the search for the desired thickness of the 10B4C thin-film coating, as well as the macrostructuring of the aluminum foil substrate used as a cathode in the production of a multiwire gas detector of thermal neutrons. [...] Read more.
This paper focuses on the search for the desired thickness of the 10B4C thin-film coating, as well as the macrostructuring of the aluminum foil substrate used as a cathode in the production of a multiwire gas detector of thermal neutrons. The impact of the 10B4C film thickness from 1.0 to 2.5 μm and of the 0.05 mm thick Λ-shaped macrostructured aluminum foil substrate with an angle at the Λ-vertex from 10 to 45 degrees, with a height from 0.5 to 4 mm and a distance between the Λ-structures from 0.1 to 0.8 relative units on the neutron registration efficiency 1.8Å, was investigated. Numerical modeling of the electrostatic field was carried out using the Elcut software package. The interaction of neutrons with the 10B4C thin-film coating was modeled using the Monte Carlo method in the Geant4 program. The optimal values of the geometrical parameters for the best neutron registration efficiency were determined. Full article
(This article belongs to the Section Bioactive Coatings and Biointerfaces)
Show Figures

Figure 1

Figure 1
<p>Schematic view of macrostructured cathode.</p>
Full article ">Figure 2
<p>Electrostatic field distribution inside detector with macrostructured cathode at <span class="html-italic">H</span> = 2 mm.</p>
Full article ">Figure 3
<p>(<b>a</b>) Efficiency of detector at <math display="inline"><semantics> <mi>α</mi> </semantics></math> <math display="inline"><semantics> <msup> <mn>10</mn> <mo>∘</mo> </msup> </semantics></math> and threshold 0 keV. (<b>b</b>) Efficiency of detector at <math display="inline"><semantics> <mi>α</mi> </semantics></math> <math display="inline"><semantics> <msup> <mn>30</mn> <mo>∘</mo> </msup> </semantics></math> and threshold 0 keV. (<b>c</b>) Efficiency of detector at <math display="inline"><semantics> <mi>α</mi> </semantics></math> <math display="inline"><semantics> <msup> <mn>10</mn> <mo>∘</mo> </msup> </semantics></math> and threshold 120 keV. (<b>d</b>) Efficiency of detector at <math display="inline"><semantics> <mi>α</mi> </semantics></math><math display="inline"><semantics> <msup> <mn>30</mn> <mo>∘</mo> </msup> </semantics></math> and threshold 120 keV.</p>
Full article ">Figure 4
<p>(<b>a</b>) Efficiency of detector at <span class="html-italic">d</span>(<math display="inline"><semantics> <msub> <mi mathvariant="normal">B</mi> <mn>4</mn> </msub> </semantics></math>C) = 1000 nm and threshold 0 keV. (<b>b</b>) Efficiency of detector at <span class="html-italic">d</span>(<math display="inline"><semantics> <msub> <mi mathvariant="normal">B</mi> <mn>4</mn> </msub> </semantics></math>C) = 1000 nm and threshold 120 keV. (<b>c</b>) Neutron using factor at <span class="html-italic">d</span>(<math display="inline"><semantics> <msub> <mi mathvariant="normal">B</mi> <mn>4</mn> </msub> </semantics></math>C) = 1000 nm and threshold 0 keV. (<b>d</b>) Neutron using factor at <span class="html-italic">d</span>(<math display="inline"><semantics> <msub> <mi mathvariant="normal">B</mi> <mn>4</mn> </msub> </semantics></math>C) = 1000 nm and threshold 120 keV.</p>
Full article ">Figure 5
<p>(<b>a</b>) Efficiency of detector at <span class="html-italic">d</span>(<math display="inline"><semantics> <msub> <mi mathvariant="normal">B</mi> <mn>4</mn> </msub> </semantics></math>C) = 1500 nm and threshold 0 keV. (<b>b</b>) Efficiency of detector at <span class="html-italic">d</span>(<math display="inline"><semantics> <msub> <mi mathvariant="normal">B</mi> <mn>4</mn> </msub> </semantics></math>C) = 1500 nm and threshold 120 keV. (<b>c</b>) Neutron using factor at <span class="html-italic">d</span>(<math display="inline"><semantics> <msub> <mi mathvariant="normal">B</mi> <mn>4</mn> </msub> </semantics></math>C) = 1500 nm and threshold 0 keV. (<b>d</b>) Neutron using factor at <span class="html-italic">d</span>(<math display="inline"><semantics> <msub> <mi mathvariant="normal">B</mi> <mn>4</mn> </msub> </semantics></math>C) = 1500 nm and threshold 120 keV.</p>
Full article ">Figure 6
<p>(<b>a</b>) Efficiency of detector at <span class="html-italic">d</span>(<math display="inline"><semantics> <msub> <mi mathvariant="normal">B</mi> <mn>4</mn> </msub> </semantics></math>C) = 2000 nm and threshold 0 keV. (<b>b</b>) Efficiency of detector at <span class="html-italic">d</span>(<math display="inline"><semantics> <msub> <mi mathvariant="normal">B</mi> <mn>4</mn> </msub> </semantics></math>C) = 2000 nm and threshold 120 keV. (<b>c</b>) Neutron using factor at <span class="html-italic">d</span>(<math display="inline"><semantics> <msub> <mi mathvariant="normal">B</mi> <mn>4</mn> </msub> </semantics></math>C) = 2000 nm and threshold 0 keV. (<b>d</b>) Neutron using factor at <span class="html-italic">d</span>(<math display="inline"><semantics> <msub> <mi mathvariant="normal">B</mi> <mn>4</mn> </msub> </semantics></math>C) = 2000 nm and threshold 120 keV.</p>
Full article ">Figure 7
<p>(<b>a</b>) Efficiency of detector at <span class="html-italic">d</span>(<math display="inline"><semantics> <msub> <mi mathvariant="normal">B</mi> <mn>4</mn> </msub> </semantics></math>C) = 2500 nm and threshold 0 keV. (<b>b</b>) Efficiency of detector at <span class="html-italic">d</span>(<math display="inline"><semantics> <msub> <mi mathvariant="normal">B</mi> <mn>4</mn> </msub> </semantics></math>C) = 2500 nm and threshold 120 keV. (<b>c</b>) Neutron using factor at <span class="html-italic">d</span>(<math display="inline"><semantics> <msub> <mi mathvariant="normal">B</mi> <mn>4</mn> </msub> </semantics></math>C) = 2500 nm and threshold 0 keV. (<b>d</b>) Neutron using factor at <span class="html-italic">d</span>(<math display="inline"><semantics> <msub> <mi mathvariant="normal">B</mi> <mn>4</mn> </msub> </semantics></math>C) = 2500 nm and threshold 120 keV.</p>
Full article ">
13 pages, 1268 KiB  
Article
Simulation and Analysis of Imaging Process of Phosphor Screens for X-Ray Imaging of Streak Tube Using Geant4-Based Monte Carlo Method
by Zichen Wang, Riyi Lin, Yuxiang Liao, Lin Tang, Zhenhua Wu, Diwei Liu, Renbin Zhong and Kaichun Zhang
Sensors 2025, 25(3), 881; https://doi.org/10.3390/s25030881 - 31 Jan 2025
Viewed by 633
Abstract
Ultrafast diagnostic technology has caused breakthroughs in fields such as inertial confinement fusion, particle accelerator research, and laser-induced phenomena. As the most widely used tool for ultrafast diagnostic technology, investigating the characteristics of streak cameras in the imaging process and streak tubes’ complex [...] Read more.
Ultrafast diagnostic technology has caused breakthroughs in fields such as inertial confinement fusion, particle accelerator research, and laser-induced phenomena. As the most widely used tool for ultrafast diagnostic technology, investigating the characteristics of streak cameras in the imaging process and streak tubes’ complex physical processes is significant for its overall development. In this work, the imaging process of a streak camera is modeled and simulated using Geant4-based Monte Carlo simulations. Based on the selected phosphor screen P43 (Gd2O2S: Tb) and charged coupled device (CCD) sensor parameters, Monte Carlo simulation models of phosphor screens and CCD sensors (We refer to the sensor parameters of the US company onsemi’s KAF-50100 sensor, but some adjustments are made during the simulation), implemented with the toolkit Geant4, are used to study the electron beam to generate fluorescence on phosphor and photoelectrons on CCD sensors. The physical process of a high-energy electron beam hitting a phosphor screen and imaging on the CCD camera is studied. Meanwhile, merits such as the luminous efficiency of the selected phosphor, spatial resolution of the phosphor screen, and spatial resolution of the selected CCD sensor are analyzed. The simulation results show that the phosphor screen and CCD sensor simulation models can accurately simulate the selected components’ performance parameters with the imaging process’ simulation results precisely reflecting the distribution of output electrons in the streak image tube. References for simulation and device selection in the subsequent research on streak cameras can be provided. Full article
(This article belongs to the Section Physical Sensors)
Show Figures

Figure 1

Figure 1
<p>Geant4 architecture block diagram.</p>
Full article ">Figure 2
<p>Simulation results of phosphor screen resolution. (<b>a</b>) 60 lp/mm electron beam; (<b>b</b>) photon distribution corresponding to 60 lp/mm electron beam; (<b>c</b>) 95 lp/mm electron beam; (<b>d</b>) photon distribution corresponding to 95 lp/mm electron beam; (<b>e</b>) 120 lp/mm electron beam; (<b>f</b>) photon distribution corresponding to 120 lp/mm electron beam.</p>
Full article ">Figure 3
<p>Simulation results of CCD imaging. (<b>a</b>) photon distribution of 30 lp/mm; (<b>b</b>) CCD imaging for photon of 30 lp/mm; (<b>c</b>) CCD imaging for photon of 60 lp/mm; (<b>d</b>) CCD imaging for photon of 90 lp/mm.</p>
Full article ">Figure 4
<p>Simulation results of streak camera imaging. (<b>a</b>) distribution of output electrons from streak image tube; (<b>b</b>) distribution of incident photons on CCD; (<b>c</b>) imaging results on CCD camera.</p>
Full article ">
11 pages, 1292 KiB  
Article
Design and Simulation of a Muon Detector Using Wavelength-Shifting Fiber Readouts for Border Security
by Anzori Sh. Georgadze
Instruments 2025, 9(1), 1; https://doi.org/10.3390/instruments9010001 - 27 Jan 2025
Viewed by 905
Abstract
Cosmic ray muon tomography is a promising method for the non-invasive inspection of shipping containers and trucks. It leverages the highly penetrating cosmic muons and their interactions with various materials to generate three-dimensional images of large and dense objects, such as inter-modal shipping [...] Read more.
Cosmic ray muon tomography is a promising method for the non-invasive inspection of shipping containers and trucks. It leverages the highly penetrating cosmic muons and their interactions with various materials to generate three-dimensional images of large and dense objects, such as inter-modal shipping containers, which are typically opaque to conventional X-ray radiography techniques. One of the key tasks of customs and border security is verifying shipping container declarations to prevent illegal trafficking, and muon tomography offers a viable solution for this purpose. Common imaging methods using muons rely on data analysis of either muon scattering or absorption–transmission. We design a compact muon tomography system with dimensions of 3 × 3 × 3 m3, consisting of 2D position-sensitive detectors. These detectors include plastic scintillators, wavelength-shifting (WLS) fibers, and SiPMs. Through light transport modeling with GEANT4, we demonstrate that the proposed detector design—featuring 1 m × 1 m scintillator plates with 2 mm2 square-shaped WLS fibers—can achieve a spatial resolution of approximately 0.7–1.0 mm. Through Monte Carlo simulations, we demonstrate that combining muon scattering and absorption data enables the rapid and accurate identification of cargo materials. In a smuggling scenario where tobacco is falsely declared as paper towel rolls, this combined analysis distinguishes the two with 3 σ confidence at a spatial resolution of 1 mm (FWHM) for the muon detector, achieving results within a scanning time of 40 s for a 20-foot shipping container. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) GEANT4 visualization of a detector composed of a plastic scintillator readout provided by 2 × 2 mm<sup>2</sup> WLS fibers. (<b>b</b>) Simulated light photons transmitted through a plastic scintillator slab and WLS fibres with GEANT4. Generated scintillation photons, shown as light green lines. The SiPMs are depicted in red.</p>
Full article ">Figure 2
<p>Photon emission spectrum of BC-404 scintillator [<a href="#B29-instruments-09-00001" class="html-bibr">29</a>] (blue colour); absorption (red colour) and emission (green colour) spectra of BCF-91A [<a href="#B27-instruments-09-00001" class="html-bibr">27</a>] and SiPM PDE (dark blue colour).</p>
Full article ">Figure 3
<p>Illustration of uniformity of XY-plane positioning estimation (<b>a</b>). Histogram of distribution of detected photons (<b>b</b>).</p>
Full article ">Figure 4
<p>Comparison of position resolution in a central position (muon hit coordinates are <span class="html-italic">x</span><sub>1</sub> = 4 mm, <span class="html-italic">y</span><sub>1</sub> = 4 mm and <span class="html-italic">x</span><sub>2</sub> = 5, <span class="html-italic">y</span><sub>2</sub> = 5 mm) (<b>a</b>) and near the detector edge (muon hit coordinates are <span class="html-italic">x</span><sub>1</sub> = 493 mm, <span class="html-italic">y</span><sub>1</sub> = 3 mm and <span class="html-italic">x</span><sub>2</sub> = 494, <span class="html-italic">y</span><sub>2</sub> = 4 mm) (<b>b</b>).</p>
Full article ">Figure 5
<p>Reconstructed muon interaction position for the central detector area in 2D (<b>a</b>) and 3D (<b>d</b>) and profiles along the <span class="html-italic">x</span>-axis and <span class="html-italic">y</span>-axis for two beam positions: <span class="html-italic">x</span><sub>1</sub> = 4 mm (<b>b</b>), <span class="html-italic">y</span><sub>1</sub> = 4 mm (<b>c</b>) and <span class="html-italic">x</span><sub>2</sub> = 5 (<b>e</b>), <span class="html-italic">y</span><sub>2</sub> = 5 mm (<b>f</b>). The red line shows a Gaussian fit performed on each profile.</p>
Full article ">Figure 6
<p>Reconstructed muon interaction position for the central detector area in 2D (<b>a</b>) and 3D (<b>d</b>) and profiles along the <span class="html-italic">x</span>-axis and <span class="html-italic">y</span>-axis for two beam positions: <span class="html-italic">x</span><sub>1</sub> = 493 mm (<b>b</b>), <span class="html-italic">y</span><sub>1</sub> = 3 mm (<b>c</b>) and <span class="html-italic">x</span><sub>2</sub> = 494 (<b>e</b>), <span class="html-italic">y</span><sub>2</sub> = 4 mm (<b>f</b>). The red line shows a Gaussian fit performed on each profile.</p>
Full article ">Figure 7
<p>A GEANT4 model of the muon tomography station is illustrated, showing the front view (<b>a</b>) and the side view (<b>b</b>). The station consists of 9 muon trackers positioned above and 9 positioned below, each constructed from plastic scintillators measuring 1 m × 1 m × 1 cm with the WLS fiber readout. A space is allocated between the plastic scintillator detectors to accommodate the front-end electronics. The muon trajectories are shown as red (<math display="inline"><semantics> <msup> <mi>μ</mi> <mo>−</mo> </msup> </semantics></math>) and blue (<math display="inline"><semantics> <msup> <mi>μ</mi> <mo>+</mo> </msup> </semantics></math>) lines.</p>
Full article ">Figure 8
<p>(<b>a</b>) Simulated 200k muons (20 s) sampled on a 10 m × 10 m surface using a CRY generator. Three-dimensional image of the PoCA reconstruction of a cigarette smuggling scenario. Two- (<b>b</b>) and one-dimensional (<b>c</b>) profiles of the PoCA image of tobacco in the compact MTS (3 × 3 × 3 m<sup>3</sup>).</p>
Full article ">Figure 9
<p>(<b>a</b>) Scatter plot showing the distributions of scattering and scattering-to-stopping ratio for paper towel rolls and tobacco. Data points for paper towel rolls are marked in blue, and data points for tobacco are marked in red. (<b>b</b>) Scatter plot showing the distributions of scattering density versus stopped muons ratio for paper towel rolls and tobacco in small MTS for 20 s scanning time.</p>
Full article ">
20 pages, 4514 KiB  
Article
Extension of the Discrete Electron Transport Capabilities of the Geant4-DNA Toolkit to MeV Energies
by Ioanna Kyriakou, Hoang N. Tran, Laurent Desorgher, Vladimir Ivantchenko, Susanna Guatelli, Giovanni Santin, Petteri Nieminen, Sebastien Incerti and Dimitris Emfietzoglou
Appl. Sci. 2025, 15(3), 1183; https://doi.org/10.3390/app15031183 - 24 Jan 2025
Viewed by 512
Abstract
The discrete physics models available in the Geant4-DNA Monte Carlo toolkit are a subject of continuous evolution and improvement in order to meet the needs of state-of-the-art radiobiological research for medical and space applications. The current capabilities of Geant4-DNA for event-by-event electron transport [...] Read more.
The discrete physics models available in the Geant4-DNA Monte Carlo toolkit are a subject of continuous evolution and improvement in order to meet the needs of state-of-the-art radiobiological research for medical and space applications. The current capabilities of Geant4-DNA for event-by-event electron transport extend up to 1 MeV. In this work, Geant4-DNA’s most accurate electron inelastic model for sub-keV energies is improved and extended up to 10 MeV via the Relativistic Plane Wave Born Approximation and other theoretical considerations. Benchmark simulations of the electronic stopping power and range of electrons in liquid water using the new model show almost excellent agreement (at the few % level) with the recommendations of the International Commission on Radiation Units and Measurements (ICRU) up to 10 MeV, offering notable improvement (by a factor of ~2) over the default Geant4-DNA inelastic model and an order-of-magnitude higher electron limit. The present development will allow Geant4-DNA users to perform electron track-structure simulations up to 10 MeV, thus, covering a wider range of radiotherapeutic applications (including FLASH-RT) as well as space applications involving MeV electrons which are not currently reachable. Full article
(This article belongs to the Section Applied Physics General)
Show Figures

Figure 1

Figure 1
<p>Schematic of the different correction terms implemented into the present inelastic model (DNA-Opt4X). To reduce the computational burden, corrections are turned “ON/OFF” only within the energy range (regime) in which they have an effect larger/smaller than 1% to the electronic stopping power (SP).</p>
Full article ">Figure 2
<p>Electron inelastic cross sections for each (outer) ionization shell (<span class="html-italic">j</span> = 1–4) and excitation level (<span class="html-italic">k</span> = 1–5) of liquid water in the energy range from threshold (~10 eV) to 10 MeV calculated by the present model (DNA-Opt4X).</p>
Full article ">Figure 3
<p>The longitudinal and transverse contribution to the electron inelastic cross section for each outer ionization shell (<b>top</b> panel) and excitation level (<b>bottom</b> panel) of liquid water in the energy range from 100 keV to 10 MeV calculated by the present model (DNA-Opt4X).</p>
Full article ">Figure 4
<p>The contribution (in %) of the longitudinal and transverse terms to the total electron inelastic cross section of liquid water in the energy range from 100 keV to 10 MeV calculated by the present model (DNA-Opt4X).</p>
Full article ">Figure 5
<p>Comparison of the total ionization cross section (<b>top</b> panel) and the total excitation cross section (<b>bottom</b> panel) of liquid water as a function of electron energy calculated by different Geant4-DNA constructors, namely, DNA-Opt2 (default), DNA-Opt4X (present/new), and DNA-Opt6. Note that each constructor has a different upper energy limit of application (see text).</p>
Full article ">Figure 6
<p>(<b>Top</b> panel) Electronic stopping power (SP) of liquid water for electrons over the energy range from 20 eV to 10 MeV simulated by different Geant4-DNA constructors, namely, DNA-Opt2 (default), DNA-Opt4X (present/new), and DNA-Opt6, and compared against the SP values of ICRU Report 90 [<a href="#B59-applsci-15-01183" class="html-bibr">59</a>]. (<b>Bottom</b> panel) Percentage difference in the various Geant4-DNA constructors from ICRU [<a href="#B59-applsci-15-01183" class="html-bibr">59</a>].</p>
Full article ">Figure 7
<p>Percentage difference in the electronic stopping power (SP) calculated by the present model (DNA-Opt4X) with and without the various correction terms from the values of ICRU Report 90 [<a href="#B59-applsci-15-01183" class="html-bibr">59</a>]. The notation is as follows: “Mott-Co” denotes the Mott–Coulomb low-energy correction (see <a href="#sec2dot3-applsci-15-01183" class="html-sec">Section 2.3</a>), “Asy” refers to the high-energy asymptotic correction (see <a href="#sec2dot4-applsci-15-01183" class="html-sec">Section 2.4</a>), “Fermi” refers to the (relativistic) Fermi density correction (see <a href="#sec2dot4-applsci-15-01183" class="html-sec">Section 2.4</a>).</p>
Full article ">Figure 8
<p>The ratio of the electronic stopping power (SP) with the low-energy corrections (i.e., the Mott and/or the Coulomb corrections) to the SP without the corrections is calculated by the present model (DNA-Opt4X).</p>
Full article ">Figure 9
<p>(<b>Top</b> panel) range (or average pathlength) of electrons in liquid water over the energy range from 20 eV up to 10 MeV simulated by different Geant4-DNA constructors, namely, DNA-Opt2 (default), DNA-Opt4X (present/new), and DNA-Opt6, and compared against the range values of ICRU Report 90 [<a href="#B59-applsci-15-01183" class="html-bibr">59</a>]; (<b>Bottom</b> panel) Percentage difference in the various Geant4-DNA constructors from ICRU [<a href="#B59-applsci-15-01183" class="html-bibr">59</a>].</p>
Full article ">
26 pages, 6283 KiB  
Article
GEANT4 Simulation of the Gamma-Ray Total Absorption Facility
by Chong Zou, Guangyuan Luan, Haotian Luo, Qiwei Zhang, Jie Ren, Xichao Ruan, Hanxiong Huang, Zhaohui Wang, Guozhu He, Jie Bao, Qi Sun, Xiaoyu Wang, Mengxiao Kang, Jincheng Wang, Yingyi Liu, Haolan Yang and Xuanbo Chen
Symmetry 2025, 17(1), 92; https://doi.org/10.3390/sym17010092 - 9 Jan 2025
Viewed by 533
Abstract
To fulfill the needs of neutron capture reaction cross-section measurement in the keV energy region in the field of nuclear astrophysics and advanced nuclear energy system development, the 4π BaF_2 Gamma-Ray Total Absorption Facility (GTAF) developed by the Key Laboratory of Nuclear Data [...] Read more.
To fulfill the needs of neutron capture reaction cross-section measurement in the keV energy region in the field of nuclear astrophysics and advanced nuclear energy system development, the 4π BaF_2 Gamma-Ray Total Absorption Facility (GTAF) developed by the Key Laboratory of Nuclear Data of the China Institute of Atomic Energy (CIAE) was transplanted and installed at the Back-streaming White Neutron Source (Back-n) of the China Spallation Neutron Source (CSNS) in 2019. A series of results has been achieved and published based on the GTAF since then, and it has been identified that the need of reducing backgrounds is becoming increasingly urgent. In order to understand the origins of backgrounds and to optimize the facilities, a detailed simulation program using GEANT4 toolkits was established and is presented in this paper. The symmetry in the geometric arrangement of the 4π BaF2 detector array plays a critical role in ensuring uniform detection efficiency and accurate reconstruction of gamma-ray spectra, which is essential for neutron capture studies. To demonstrate the availability of the proven codes, several practical examples of assisting the process of experimental data and helping verify the optimization proposition are also shown in this paper. Full article
(This article belongs to the Section Physics)
Show Figures

Figure 1

Figure 1
<p>GTAF detector array and associated facilities installed in the Hall-2 of Back-n CSNS [<a href="#B6-symmetry-17-00092" class="html-bibr">6</a>].</p>
Full article ">Figure 2
<p>(<b>a</b>) Principle of isotope de-excitation: A captured neutron excites the nucleus, which de-excites through multiple gamma-ray emissions. (<b>b</b>) Schema of response multiplicity: gamma rays may be fully absorbed by one detector crystal (process 1) or undergo Compton scattering, triggering multiple crystals (processes 2,3).</p>
Full article ">Figure 3
<p>Principle of pile-up energy [<a href="#B6-symmetry-17-00092" class="html-bibr">6</a>].</p>
Full article ">Figure 4
<p>General data flow of simulation codes.</p>
Full article ">Figure 5
<p>Typical simulation reconstruction. (<b>a</b>) Mass plan of the Hall-2 geometry using CSG + CAD method; (<b>b</b>) central zone for the detector and its associated geometry using CSG + CAD method; (<b>c</b>) <math display="inline"><semantics> <mrow> <mi mathvariant="normal">B</mi> <mi mathvariant="normal">a</mi> <msub> <mrow> <mi mathvariant="normal">F</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> </mrow> </semantics></math> crystal geometry using CSG/CAD method; (<b>d</b>) sample tray/support geometry using CSG/CAD method; (<b>e</b>) bellow geometry using CSG + CAD method; (<b>f</b>) cage support geometry using CAD method; (<b>g</b>) <math display="inline"><semantics> <mrow> <msub> <mrow> <mi mathvariant="normal">B</mi> <mi mathvariant="normal">a</mi> <mi mathvariant="normal">F</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> </mrow> </semantics></math> pentagonal pyramid; (<b>h</b>) <math display="inline"><semantics> <mrow> <msub> <mrow> <mi mathvariant="normal">B</mi> <mi mathvariant="normal">a</mi> <mi mathvariant="normal">F</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> </mrow> </semantics></math> hexagonal pyramid.</p>
Full article ">Figure 5 Cont.
<p>Typical simulation reconstruction. (<b>a</b>) Mass plan of the Hall-2 geometry using CSG + CAD method; (<b>b</b>) central zone for the detector and its associated geometry using CSG + CAD method; (<b>c</b>) <math display="inline"><semantics> <mrow> <mi mathvariant="normal">B</mi> <mi mathvariant="normal">a</mi> <msub> <mrow> <mi mathvariant="normal">F</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> </mrow> </semantics></math> crystal geometry using CSG/CAD method; (<b>d</b>) sample tray/support geometry using CSG/CAD method; (<b>e</b>) bellow geometry using CSG + CAD method; (<b>f</b>) cage support geometry using CAD method; (<b>g</b>) <math display="inline"><semantics> <mrow> <msub> <mrow> <mi mathvariant="normal">B</mi> <mi mathvariant="normal">a</mi> <mi mathvariant="normal">F</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> </mrow> </semantics></math> pentagonal pyramid; (<b>h</b>) <math display="inline"><semantics> <mrow> <msub> <mrow> <mi mathvariant="normal">B</mi> <mi mathvariant="normal">a</mi> <mi mathvariant="normal">F</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> </mrow> </semantics></math> hexagonal pyramid.</p>
Full article ">Figure 6
<p>Deposited energy spectrum under different basic physical models.</p>
Full article ">Figure 7
<p>Simulation of neutron beam spot at Back-n of CSNS.</p>
Full article ">Figure 8
<p>GUI interface of GTAF simulation pre-processing program.</p>
Full article ">Figure 9
<p>Simulated calibration energy spectrum (pile-up energy) for different sources. The blue line represents <sup>137</sup>Cs, the black line represents <sup>60</sup>Co, and the red line represents <sup>22</sup>Na.</p>
Full article ">Figure 10
<p>Demonstration of multiplicity identification for <sup>60</sup>Co source simulated calibration experiment.</p>
Full article ">Figure 10 Cont.
<p>Demonstration of multiplicity identification for <sup>60</sup>Co source simulated calibration experiment.</p>
Full article ">Figure 11
<p>Simulation results of <sup>197</sup>Au sample response to the 4.9 eV monoenergetic neutron beam.</p>
Full article ">Figure 12
<p>Demonstrations of reaction channel discrimination.</p>
Full article ">Figure 13
<p>Demonstration of primary analysis of background.</p>
Full article ">Figure 14
<p>Background analysis and comparison.</p>
Full article ">Figure 15
<p>Mass plan with the optimized structure.</p>
Full article ">Figure 15 Cont.
<p>Mass plan with the optimized structure.</p>
Full article ">
17 pages, 11667 KiB  
Article
Silicon Drift Detectors for the Measurement and Reconstruction of Beta Spectra
by Andrea Nava, Leonardo Bernardini, Matteo Biassoni, Tommaso Bradanini, Marco Carminati, Giovanni De Gregorio, Carlo Fiorini, Giulio Gagliardi, Peter Lechner, Riccardo Mancino and Chiara Brofferio
Sensors 2024, 24(24), 8202; https://doi.org/10.3390/s24248202 - 22 Dec 2024
Viewed by 779
Abstract
The ASPECT-BET project, or An sdd-SPECTrometer for BETa decay studies, aims to develop a novel technique for the precise measurement of forbidden beta spectra in the 10 keV–1 MeV range. This technique employs a Silicon Drift Detector (SDD) as the main spectrometer with [...] Read more.
The ASPECT-BET project, or An sdd-SPECTrometer for BETa decay studies, aims to develop a novel technique for the precise measurement of forbidden beta spectra in the 10 keV–1 MeV range. This technique employs a Silicon Drift Detector (SDD) as the main spectrometer with the option of a veto system to reject events exhibiting only partial energy deposition in the SDD. A precise understanding of the spectrometer’s response to electrons is crucial for accurately reconstructing the theoretical shape of the beta spectrum. To compute this response, GEANT4 simulations optimized for low-energy electron interactions are used and validated with a custom-made electron gun. In this article we present the performance of these simulations in reconstructing the electron spectra measured with SDDs of a 109Cd monochromatic source, both in vacuum and in air. The allowed beta spectrum of a 14C source was also measured and analyzed, proving that this system is suitable for the application in ASPECT-BET. Full article
Show Figures

Figure 1

Figure 1
<p>47-pixel SDD matrix used for all the measurements here reported (<b>left</b>). Scheme of the 47 pixels: only the 7 red ones were acquired (<b>right</b>).</p>
Full article ">Figure 2
<p>Vacuum chamber in the Milano-Bicocca laboratory. The main detectors operated in this setup, a 47-pixel SDD matrix and a Pixet, are indicated. The e-gun, attached to a xy-movable stage, is also highlighted.</p>
Full article ">Figure 3
<p>(<b>Left</b>): picture of the e-gun. (<b>Center</b>): anode with 7 LEDs used to illuminate the cathode and the aluminum cylinder with UV light. (<b>Right</b>): gold-coated cathode used to collimate the electron beam. The aluminum cylinder, where electrons are produced, is visible in its center.</p>
Full article ">Figure 4
<p>(<b>Left</b>): CAD drawing of the e-gun. The most important components are highlighted. (<b>Right</b>): COMSOL simulation of the electron beam produced with the e-gun.</p>
Full article ">Figure 5
<p>(<b>Top</b>): measurement of the e-gun beam spot performed with the Pixet. A beam size of ∼0.5 mm and a rate <math display="inline"><semantics> <msup> <mn>10</mn> <mn>4</mn> </msup> </semantics></math> electrons/s were found. (<b>Bottom</b>): measurement of a 10 keV electron spectrum acquired with an SDD. Only the central part of the pixel was hit with the e-gun beam. The best fit of the spectrum done with the detector model is also shown.</p>
Full article ">Figure 6
<p>CAD schematics of the setup used in the measurements. Section of source and detector (<b>top</b>), and side view (<b>bottom</b>).</p>
Full article ">Figure 7
<p>M1 Data acquired with the SDD main pixel in a 3-h measurement. The energy of the X-ray peaks and the tag of the different IC electrons are shown.</p>
Full article ">Figure 8
<p>Data-MC comparison for different values of the effective Mylar thickness (<b>top</b>). Reduced <math display="inline"><semantics> <msup> <mi>χ</mi> <mn>2</mn> </msup> </semantics></math> as a function of the effective Mylar thickness (<b>bottom</b>). The best-fit value of 7.2 μm is extracted through a parabolic fit.</p>
Full article ">Figure 9
<p>Best fit of MC prediction to the data set acquired with the <sup>109</sup>Cd source in vacuum. The fit is done only in the non-shaded area.</p>
Full article ">Figure 10
<p>Comparison of MC prediction to the dataset acquired with the <sup>109</sup>Cd source in air. The comparison is done only in the non-shaded area.</p>
Full article ">Figure 11
<p>Data acquired with the main SDD in a 6 h measurement using a <sup>14</sup>C source (<b>top</b>). <math display="inline"><semantics> <msup> <mi>χ</mi> <mn>2</mn> </msup> </semantics></math> as a function of the effective Mylar thickness for the two models: the Fermi theory prediction and the one including the experimental shape factor (<b>bottom</b>). The best fit for the effective Mylar thickness is 4.5 μm.</p>
Full article ">Figure 12
<p>Fit to the data set acquired with the <sup>14</sup>C source in vacuum using the MC prediction. The fit is performed only in the non-shaded area. The theoretical input is also shown.</p>
Full article ">Figure 13
<p>Fits obtained by varying <math display="inline"><semantics> <mi>λ</mi> </semantics></math> (<b>top</b>), the baseline resolution <math display="inline"><semantics> <mi>σ</mi> </semantics></math> (<b>center</b>) and the charge cloud width in Si (<b>bottom</b>).</p>
Full article ">Figure 13 Cont.
<p>Fits obtained by varying <math display="inline"><semantics> <mi>λ</mi> </semantics></math> (<b>top</b>), the baseline resolution <math display="inline"><semantics> <mi>σ</mi> </semantics></math> (<b>center</b>) and the charge cloud width in Si (<b>bottom</b>).</p>
Full article ">Figure 14
<p>Fits obtained by varying the production cut for secondaries in GEANT4 (<b>top</b>) or the physics list used (<b>bottom</b>) in the simulations.</p>
Full article ">
13 pages, 2277 KiB  
Article
Modelling Potential Candidates for Targeted Auger Therapy
by Conor M. J. Buchanan, Eric O. Aboagye, Lee J. Evitts, Michael J. D. Rushton and Tim A. D. Smith
Biophysica 2024, 4(4), 711-723; https://doi.org/10.3390/biophysica4040046 - 18 Dec 2024
Viewed by 845
Abstract
Targeted Auger emitters are being considered as a cancer treatment owing to the high linear energy transfer of Auger electrons. When targeted to cancers, this allows for a highly efficient treatment with a low risk of damage to surrounding healthy tissue. The purpose [...] Read more.
Targeted Auger emitters are being considered as a cancer treatment owing to the high linear energy transfer of Auger electrons. When targeted to cancers, this allows for a highly efficient treatment with a low risk of damage to surrounding healthy tissue. The purpose of this study was to determine the most DNA-damaging Auger emitters from a range of radionuclides, some of which are clinically utilised. A Monte Carlo method-based software (Geant4-DNA version 10.7) was used to determine the energy deposition and number of DNA double-strand breaks from Auger (and internal conversion) electrons imposed on a tetranucleosome. The Auger emitters, 119Sb and 103Pd, have similar or slightly greater damaging properties compared to 123I, 111In, and 89Zr. 193mPt demonstrated the greatest therapeutic potency. Whilst 125I was highly damaging, its relatively long half-life (60 days) makes it less desirable for clinical use. Geant4-DNA modelling identified the radionuclide 193mPt as being highly favourable for use in radiotherapy. Full article
Show Figures

Figure 1

Figure 1
<p>A visual representation of targeted Auger therapy showing the radiopharmaceutical being taken up by a cancer cell into the cell nucleus where Auger electrons target the cell DNA.</p>
Full article ">Figure 2
<p>Detailed illustration of Auger electron emission following either electron capture or internal conversion processes.</p>
Full article ">Figure 3
<p>Track length, projected length, and penetration of incident electrons on a 1m radius sphere of liquid water.</p>
Full article ">Figure 4
<p>The simulated dose on the tetranucleosome from five prospective Auger-emitting radionuclides. Each is split to show the contribution from Auger electrons, conversion electrons, and β<sup>−</sup> particles.</p>
Full article ">Figure 5
<p>The equivalent dose applied on the tetranucleosome by four radionuclides currently used in nuclear medicine showing the contributions from Auger, conversion electrons, and β<sup>−</sup> particles.</p>
Full article ">Figure 6
<p>The number of double-strand breaks induced by low-energy electrons of increasing energy.</p>
Full article ">Figure 7
<p>The number of double-strand breaks induced in the tetranucleosome showing contributions from Auger electrons, conversion electrons, and β<sup>−</sup> particles emitted from novel radionuclides.</p>
Full article ">Figure 8
<p>The number of double-strand breaks induced on the tetranucleosome by radionuclides currently used in nuclear medicine showing contributions from Auger electrons, conversion electrons, and β<sup>−</sup> particles.</p>
Full article ">
16 pages, 6075 KiB  
Article
A Comparative Study of Neutron Shielding Performance in Al-Based Composites Reinforced with Various Boron-Containing Particles for Radiotherapy: A Monte Carlo Simulation
by Shiyan Yang, Yupeng Yao, Hanlong Wang and Hai Huang
Nanomaterials 2024, 14(21), 1696; https://doi.org/10.3390/nano14211696 - 23 Oct 2024
Viewed by 1194
Abstract
This study aimed to assess and compare the shielding performance of boron-containing materials for neutrons generated in proton therapy and used in boron neutron capture therapy (BNCT). Five composites, including AlB2, Al-B4C, Al-TiB2, Al-BN, and Al-TiB2 [...] Read more.
This study aimed to assess and compare the shielding performance of boron-containing materials for neutrons generated in proton therapy and used in boron neutron capture therapy (BNCT). Five composites, including AlB2, Al-B4C, Al-TiB2, Al-BN, and Al-TiB2-BN, were selected as shielding materials, with concrete used as a benchmark. The mass fraction of boron compounds in these materials ranged from 10% to 50%. The Monte Carlo toolkit Geant4 was employed to calculate shielding parameters, including neutron ambient dose equivalent, dose values, and macroscopic cross-section. Results indicated that, compared to concrete, these boron-containing materials more effectively absorb thermal neutrons. When the boron compound exceeds 30 wt.%, these materials exhibit better shielding performance than concrete of the same thickness for neutrons generated by protons. For a given material, its shielding capability increases with boron content. Among the five materials when the material thickness and boron compound content are the same, the shielding performance for neutrons generated by protons, from best to worst, is as follows: Al-TiB2, Al-B4C, AlB2, Al-TiB2-BN, and Al-BN. For BNCT, the shielding performance from best to worst is in the following order: Al-B4C, AlB2, Al-TiB2, Al-TiB2-BN, and Al-BN. The results of this study provide references and guidelines for the selection and optimization of neutron shielding materials in proton therapy and BNCT facilities. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Schematic diagram of the simulation setup for secondary neutrons produced by protons. (<b>b</b>) Three proton SOBPs with different modulation ranges were used as primary sources.</p>
Full article ">Figure 2
<p>Secondary neutron spectral fluence produced by protons with different energies recorded at various detector angles. All figures share the same legend.</p>
Full article ">Figure 3
<p>Neutron energy spectrum from the p-Li reaction at a proton energy of 2.8 MeV.</p>
Full article ">Figure 4
<p>The neutron spectral fluence recorded in the ICRU sphere after passing through different shielding materials of various thicknesses. The first row shows the energy spectra of the three studied neutrons at the source plan.</p>
Full article ">Figure 5
<p>Elastic and inelastic collision cross sections of neutrons in different elements that constitute the shielding materials.</p>
Full article ">Figure 6
<p>Neutron ambient dose equivalent values for different shielding materials. The left, middle, and right columns show the results for n-BNCT, n-SOBP1, and n-SOBP3, respectively. Each row uses the same legend.</p>
Full article ">Figure 7
<p>Dose distribution of different particles after passing through various thicknesses of concrete for the three studied neutron spectra.</p>
Full article ">Figure 8
<p>Comparison of total dose distribution for different materials at varying thicknesses and boron contents. The left, middle, and right columns present the results for n-BNCT, n-SOBP1, and n-SOBP3, respectively. Each row uses the same legend.</p>
Full article ">Figure 9
<p>Neutron macroscopic cross-section for various materials with different boron contents at 10 cm for n-BNCT (<b>top row</b>) and 50 cm thickness for n-SOBP1 (<b>middle row</b>) and n-SOBP3 (<b>bottom row</b>). The gray line represents the results of concrete at the corresponding thickness. Each column uses the same legend.</p>
Full article ">
11 pages, 1991 KiB  
Article
Research on B4C/PEEK Composite Material Radiation Shielding
by Hongxia Li, Hongping Guo, Hui Tu, Xiao Chen and Xianghua Zeng
Polymers 2024, 16(20), 2902; https://doi.org/10.3390/polym16202902 - 15 Oct 2024
Viewed by 965
Abstract
There are various types of charged particles in the space environment, which can cause different types of radiation damage to materials and devices, leading to on-orbit failures and even accidents for spacecraft. Developing lightweight and efficient radiation-shielding materials is an effective approach to [...] Read more.
There are various types of charged particles in the space environment, which can cause different types of radiation damage to materials and devices, leading to on-orbit failures and even accidents for spacecraft. Developing lightweight and efficient radiation-shielding materials is an effective approach to improving the inherent protection of spacecraft. The protective performance of different materials against proton and electron spectra in the Earth’s radiation belts is evaluated using a Geant4 simulation. Based on the simulation results, suitable hardening components were selected to design composite materials, and B4C/PEEK composites with different B4C contents were successfully prepared. The experimental results demonstrate that the simulated and experimental results for the electron, proton and neutron shielding performance of the B4C/PEEK composites are consistent. These composites exhibit excellent radiation shielding capabilities against electrons, protons and neutrons, and the radiation protection performance improves with increasing B4C content in the B4C/PEEK composite materials. Full article
(This article belongs to the Special Issue Advances in Functional Polymer Nanocomposites)
Show Figures

Figure 1

Figure 1
<p>Differential cumulative flux spectra of LEO orbital protons protected by different materials and the same mass thickness: (<b>a</b>) differential cumulative flux spectra after different material and thickness protection; (<b>b</b>) proton protection effect of different materials with the same mass thickness.</p>
Full article ">Figure 2
<p>Differential cumulative flux spectra of secondary neutrons and gamma rays after LEO orbital protons pass through protective layers of different thicknesses: (<b>a</b>) secondary neutrons; (<b>b</b>) <span class="html-italic">γ</span> rays.</p>
Full article ">Figure 3
<p>The differential cumulative flux spectra of the second gamma of GEO orbital electrons through 2 mm protective layer of different materials.</p>
Full article ">Figure 4
<p>Simulation calculation results of 1 MeV electron and 16 MeV proton protection rate of materials: (<b>a</b>) protective effect of different materials; (<b>b</b>) electronic protective effect of B<sub>4</sub>C/PEEK at different B<sub>4</sub>C addition levels; and (<b>c</b>) 16 MeV proton protective effect of B<sub>4</sub>C/PEEK at different B<sub>4</sub>C addition levels.</p>
Full article ">Figure 5
<p>Preparation flow diagram of B<sub>4</sub>C/PEEK composite materials.</p>
Full article ">Figure 6
<p>Distribution curves of film dose tablets along the thickness of B<sub>4</sub>C/PEEK under electron irradiation.</p>
Full article ">Figure 7
<p>Proton transmittance of B<sub>4</sub>C/PEEK composites changes with energy: (<b>a</b>) 10 wt% B<sub>4</sub>C/PEEK, (<b>b</b>) 20 wt% B<sub>4</sub>C/PEEK.</p>
Full article ">Figure 8
<p>Neutron protection curves of B<sub>4</sub>C/PEEK composite materials with different B<sub>4</sub>C addition amounts: (<b>a</b>) variation curves of neutron count and wavelength; (<b>b</b>) change curves of wavelength and protection effect.</p>
Full article ">
23 pages, 7773 KiB  
Article
Search for True Ternary Fission in Reaction 40Ar + 208Pb
by Md Ashaduzzaman, Antonio Di Nitto, Emanuele Vardaci, Giovanni La Rana, Pia Antonella Setaro, Tathagata Banerjee, Antonio Vanzanella and Giuseppe Alifano
Appl. Sci. 2024, 14(18), 8522; https://doi.org/10.3390/app14188522 - 21 Sep 2024
Viewed by 1000
Abstract
True ternary fission, the fission of a nucleus into three fragments of nearly equal mass, is an elusive and poorly known process influenced by shell effects. An increase in the probability of this process with respect to binary fission, which is very low [...] Read more.
True ternary fission, the fission of a nucleus into three fragments of nearly equal mass, is an elusive and poorly known process influenced by shell effects. An increase in the probability of this process with respect to binary fission, which is very low in spontaneous and neutron-induced fission, has been envisaged. Heavy-ion-induced reactions are adopted due to the possibility of an increase in the fissility parameter and the excitation energy of the compound nuclei. Nuclei with mass number around A = 250, accessible in heavy-ion-induced reactions, are favorable and should be investigated. It is still debated if the process takes place in a single step, direct ternary fission, or in a two step, sequential ternary fission. The purpose of this work is to define experimental conditions and observables that allow the disentangling of the products from the direct and sequential ternary fission, as well as from the usual most probable binary fission. This step is essential for gaining insights into the ternary fission dynamics and the binary to ternary fission competition. The method proposed here is for simulating the kinematics of the ternary and binary fission processes to compute the energy distributions and angular correlations of direct and sequential ternary fission products, as well as those of binary fission. The reaction taken as a benchmark is 40Ar + 208Pb at 230 MeV and is supposed to form the 248Fm* compound nucleus. The simulation results have been filtered by considering the response function of a multi-coincidence detection system virtually constructed using the Geant4 simulation toolkit. The simulations support the possibility of separating the products of different multimodal fission decays with the proposed setup that consequently represents an effective tool to obtain insights into ternary fission from the observables selected. Full article
(This article belongs to the Section Applied Physics General)
Show Figures

Figure 1

Figure 1
<p>In-plane kinematic plots of direct ternary fission (DTF-<b>top</b>) and sequential ternary fission(STF-<b>bottom</b>). In the DTF a simultaneous formation of the <math display="inline"><semantics> <mrow> <msub> <mi>A</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>A</mi> <mn>2</mn> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <msub> <mi>A</mi> <mn>3</mn> </msub> </semantics></math> fragments occurs. In STF the first binary scission produces the <math display="inline"><semantics> <msub> <mi>A</mi> <mn>1</mn> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mi>A</mi> <mn>23</mn> </msub> </semantics></math> fragments, the subsequent binary scission of the <math display="inline"><semantics> <msub> <mi>A</mi> <mn>23</mn> </msub> </semantics></math> fragment produces the final products <math display="inline"><semantics> <msub> <mi>A</mi> <mn>2</mn> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mi>A</mi> <mn>3</mn> </msub> </semantics></math>.</p>
Full article ">Figure 2
<p>Dalitz plot of the potential energies for collinear (<b>a</b>) and equatorial (<b>b</b>) ternary fragmentations of <sup>248</sup>Fm (Z = 100). The potential energies are calculated as a function of the charge numbers with the constraint for the fragment masses <math display="inline"><semantics> <mrow> <msub> <mi>A</mi> <mn>1</mn> </msub> <mo>≥</mo> <msub> <mi>A</mi> <mn>2</mn> </msub> <mo>≥</mo> <msub> <mi>A</mi> <mn>3</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 3
<p>Potential Energy of a ternary fragmentation involving two out of three double magic nuclei in collinear and equatorial configurations.</p>
Full article ">Figure 4
<p>Kinetic energies of light fragments <math display="inline"><semantics> <msub> <mi>E</mi> <mn>3</mn> </msub> </semantics></math> presented as a function of <math display="inline"><semantics> <msub> <mi>E</mi> <mn>1</mn> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mi>E</mi> <mn>2</mn> </msub> </semantics></math> by assuming DTF (<b>a</b>) and STF (<b>b</b>) mechanisms. The calculations have been performed with the conditions: <math display="inline"><semantics> <mrow> <msub> <mi>A</mi> <mn>1</mn> </msub> <mo>≥</mo> <msub> <mi>A</mi> <mn>2</mn> </msub> <mo>≥</mo> <msub> <mi>A</mi> <mn>3</mn> </msub> </mrow> </semantics></math>, emission of <math display="inline"><semantics> <msub> <mi>A</mi> <mn>1</mn> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mi>A</mi> <mn>2</mn> </msub> </semantics></math> on the opposite side with respect to the beam direction and wide range for the <math display="inline"><semantics> <msub> <mi>A</mi> <mn>1</mn> </msub> </semantics></math> energies. See the text for more details.</p>
Full article ">Figure 5
<p>Emission angles of light fragment <math display="inline"><semantics> <msub> <mi>θ</mi> <mn>3</mn> </msub> </semantics></math> as a function of the heavy (<math display="inline"><semantics> <msub> <mi>θ</mi> <mn>1</mn> </msub> </semantics></math>) and medium (<math display="inline"><semantics> <msub> <mi>θ</mi> <mn>2</mn> </msub> </semantics></math>) mass angles in the laboratory system calculated considering the DTF (<b>a</b>) and STF (<b>b</b>) mechanisms. In the calculations, <math display="inline"><semantics> <mrow> <msub> <mi>A</mi> <mn>1</mn> </msub> <mo>≥</mo> <msub> <mi>A</mi> <mn>2</mn> </msub> <mo>≥</mo> <msub> <mi>A</mi> <mn>3</mn> </msub> </mrow> </semantics></math> and two of the three fragments are doubly magic nuclei. See the text for more details.</p>
Full article ">Figure 6
<p>Kinetic energies of the product of <sup>132</sup>Sn + <sup>68</sup>Zn + <sup>48</sup>Ca tripartition. The products of the DTF and STF mechanisms are shown in the left and right panels, respectively.</p>
Full article ">Figure 7
<p>Angular correlations of the <sup>132</sup>Sn + <sup>68</sup>Zn + <sup>48</sup>Ca tripartition. The DTF (blue) and STF (pink) angular correlations <math display="inline"><semantics> <msub> <mi>θ</mi> <mn>2</mn> </msub> </semantics></math> vs. <math display="inline"><semantics> <msub> <mi>θ</mi> <mn>3</mn> </msub> </semantics></math> are presented for the fixed heavy fragment angle <math display="inline"><semantics> <msub> <mi>θ</mi> <mn>1</mn> </msub> </semantics></math> described at the top of each panel.</p>
Full article ">Figure 8
<p>DTF angular correlations of the <sup>132</sup>Sn + <sup>68</sup>Zn + <sup>48</sup>Ca tripartition for <math display="inline"><semantics> <msub> <mi>θ</mi> <mn>2</mn> </msub> </semantics></math> = <math display="inline"><semantics> <mrow> <mo>−</mo> <msub> <mi>θ</mi> <mn>1</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 9
<p>Energy correlations of the <sup>132</sup>Sn + <sup>68</sup>Zn + <sup>48</sup>Ca tripartition. The DTF (blue) and STF (pink) energy correlations <math display="inline"><semantics> <msub> <mi>E</mi> <mn>2</mn> </msub> </semantics></math> vs. <math display="inline"><semantics> <msub> <mi>E</mi> <mn>3</mn> </msub> </semantics></math> are presented for the fixed heavy fragment angle <math display="inline"><semantics> <msub> <mi>θ</mi> <mn>1</mn> </msub> </semantics></math> described at the top of each panel.</p>
Full article ">Figure 10
<p>Energies and angles of the light fragments from the <sup>132</sup>Sn + <sup>68</sup>Zn + <sup>48</sup>Ca tripartition. Under the condition of angle symmetric emission of the medium and heavy fragment at <math display="inline"><semantics> <mrow> <msub> <mi>θ</mi> <mn>1</mn> </msub> <mo>=</mo> <mo>−</mo> <msub> <mi>θ</mi> <mn>2</mn> </msub> <mo>=</mo> <mo>−</mo> <mn>70</mn> <mo>°</mo> </mrow> </semantics></math>, the angular and energy distributions of the <sup>48</sup>Ca fragments are shown in panels (<b>a</b>,<b>b</b>), while the angle vs. energy correlations are shown in panel (<b>c</b>), assuming both DTF (black) and STF (pink) processes. For details on energy, angular intervals, and steps considered in calculations, see the text.</p>
Full article ">Figure 11
<p>Angular correlations of BF events. The distributions of fragment 2 (<math display="inline"><semantics> <msub> <mi>θ</mi> <mn>2</mn> </msub> </semantics></math>) vs. fragment 1 (<math display="inline"><semantics> <msub> <mi>θ</mi> <mn>1</mn> </msub> </semantics></math>) for symmetric (blue line) and asymmetric fission (orange line). The green boxes represent the two angular ranges covered by the detection setup proposed in <a href="#sec2dot4-applsci-14-08522" class="html-sec">Section 2.4</a>.</p>
Full article ">Figure 12
<p>Schematic drawing of the detection apparatus consisting of three different detector types. The blue line superimposed on the z-axis (blue arrow) identifies the beam direction. See the text for more details.</p>
Full article ">Figure 13
<p>Velocity distributions of coincident events (<math display="inline"><semantics> <mrow> <msub> <mi>θ</mi> <mn>2</mn> </msub> <mo>=</mo> <mo>−</mo> <msub> <mi>θ</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>70</mn> <mo>°</mo> </mrow> </semantics></math>). The spectra at <math display="inline"><semantics> <msub> <mi>θ</mi> <mn>1</mn> </msub> </semantics></math> (<b>a</b>) and <math display="inline"><semantics> <msub> <mi>θ</mi> <mn>2</mn> </msub> </semantics></math> (<b>b</b>) correspond to the heavier and middle mass fragments in the case of TF, respectively. In asymmetric BF the heavier fragment is measured at <math display="inline"><semantics> <msub> <mi>θ</mi> <mn>1</mn> </msub> </semantics></math>. The location of the coincident events in the two-dimensional matrix <math display="inline"><semantics> <msub> <mi>v</mi> <mn>1</mn> </msub> </semantics></math> vs. <math display="inline"><semantics> <msub> <mi>v</mi> <mn>2</mn> </msub> </semantics></math> are shown in panel (<b>c</b>).</p>
Full article ">Figure 14
<p>Energy spectra measured with the two stages telescope detector at <math display="inline"><semantics> <mrow> <msub> <mi>θ</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>90</mn> <mo>°</mo> <mo>±</mo> <mn>15</mn> <mo>°</mo> </mrow> </semantics></math>. The spectra are simulated for the thin <math display="inline"><semantics> <mrow> <mo>Δ</mo> <mi>E</mi> </mrow> </semantics></math> (<b>a</b>) and thick <math display="inline"><semantics> <msub> <mi>E</mi> <mrow> <mi>r</mi> <mi>e</mi> <mi>s</mi> </mrow> </msub> </semantics></math> (<b>b</b>) detectors in the case of the 4 processes taken into account.</p>
Full article ">Figure 15
<p>Energy-velocity matrix of fragments emitted at <math display="inline"><semantics> <mrow> <mn>60</mn> <mo>°</mo> <mo>±</mo> <mn>5</mn> <mo>°</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 16
<p><math display="inline"><semantics> <mrow> <mo>Δ</mo> <mi>E</mi> <mo>−</mo> <mi>E</mi> </mrow> </semantics></math> matrix of nuclei with A = 48 and different atomic numbers. The width of the distribution in <math display="inline"><semantics> <mrow> <mo>Δ</mo> <mi>E</mi> </mrow> </semantics></math> depends mainly on the energy resolution of the detectors. Here we assume for both detectors an energy resolution (FWHM) of 900 keV.</p>
Full article ">
18 pages, 2273 KiB  
Article
Optimization of the Pixel Design for Large Gamma Cameras Based on Silicon Photomultipliers
by Carolin Wunderlich, Riccardo Paoletti and Daniel Guberman
Sensors 2024, 24(18), 6052; https://doi.org/10.3390/s24186052 - 19 Sep 2024
Viewed by 1029
Abstract
Most single-photon emission computed tomography (SPECT) scanners employ a gamma camera with a large scintillator crystal and 50–100 large photomultiplier tubes (PMTs). In the past, we proposed that the weight, size and cost of a scanner could be reduced by replacing the PMTs [...] Read more.
Most single-photon emission computed tomography (SPECT) scanners employ a gamma camera with a large scintillator crystal and 50–100 large photomultiplier tubes (PMTs). In the past, we proposed that the weight, size and cost of a scanner could be reduced by replacing the PMTs with large-area silicon photomultiplier (SiPM) pixels in which commercial SiPMs are summed to reduce the number of readout channels. We studied the feasibility of that solution with a small homemade camera, but the question on how it could be implemented in a large camera remained open. In this work, we try to answer this question by performing Geant4 simulations of a full-body SPECT camera. We studied how the pixel size, shape and noise could affect its energy and spatial resolution. Our results suggest that it would be possible to obtain an intrinsic spatial resolution of a few mm FWHM and an energy resolution at 140 keV close to 10%, even if using pixels more than 20 times larger than standard commercial SiPMs of 6 × 6 mm2. We have also found that if SiPMs are distributed following a honeycomb structure, the spatial resolution is significantly better than if using square pixels distributed in a square grid. Full article
Show Figures

Figure 1

Figure 1
<p>Scheme showing the different LASiP configurations (1–8).</p>
Full article ">Figure 2
<p>A simulated FFI using LASiP 2 without noise: (<b>a</b>) raw reconstruction, (<b>b</b>) after spatial linearity correction and (<b>c</b>) after uniformity correction. The blue dashed lines show the LASiP edges.</p>
Full article ">Figure 3
<p>Collected charge (relative to the total number of collected photons) as a function of the number of SiPMs used to collect the charge.</p>
Full article ">Figure 4
<p>Energy resolution as a function of the number of LASiPs used to collect the charge (LASiP 1–4) under different SiPM noise levels. The top <span class="html-italic">x</span> axis shows the total number of SiPMs employed.</p>
Full article ">Figure 5
<p>Reconstructed image of a capillary source at <math display="inline"><semantics> <mrow> <mi>x</mi> <mo>=</mo> <mo>−</mo> <mn>12</mn> <mo>,</mo> <mn>0</mn> <mo>,</mo> <mn>10</mn> </mrow> </semantics></math> mm for LASiP 1–4 at room-SiPM noise. The blue dashed lines mark the borders of the LASiPs.</p>
Full article ">Figure 6
<p>Normalized projections in the <span class="html-italic">x</span> axis of the reconstructed images of <a href="#sensors-24-06052-f005" class="html-fig">Figure 5</a>. The projections were performed in a narrow region around the center of the pixel.</p>
Full article ">Figure 7
<p>Mean FWHM of the reconstructed capillary as a function of the normalized distance to the pixel center <math display="inline"><semantics> <mover accent="true"> <mi>d</mi> <mo stretchy="false">^</mo> </mover> </semantics></math> (<math display="inline"><semantics> <mrow> <mover accent="true"> <mi>d</mi> <mo stretchy="false">^</mo> </mover> <mo>=</mo> <mn>0</mn> </mrow> </semantics></math> describes a capillary passing by the center of a pixel, <math display="inline"><semantics> <mrow> <mover accent="true"> <mi>d</mi> <mo stretchy="false">^</mo> </mover> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math> describes a capillary at the edge of a pixel).</p>
Full article ">Figure 8
<p>Reconstructed image of a simile-Derenzo source (the diameter of the circles are 1, 2, 3, 4, 6 and 8 mm) at room-SiPM noise. The dashed lines mark the edges of the LASiPs.</p>
Full article ">Figure 9
<p>Normalized projections in the <span class="html-italic">x</span> axis of the reconstructed images of capillaries passing through the center of the LASiP 3, LASiP 5 and LASiP 8 pixels. The projections were centered at different <span class="html-italic">y</span> values: (<b>a</b>) <math display="inline"><semantics> <mrow> <mi>y</mi> <mo>=</mo> <mn>5</mn> </mrow> </semantics></math> mm, (<b>b</b>) <math display="inline"><semantics> <mrow> <mi>y</mi> <mo>=</mo> <mn>10</mn> </mrow> </semantics></math> mm, (<b>c</b>) <math display="inline"><semantics> <mrow> <mi>y</mi> <mo>=</mo> <mn>20</mn> </mrow> </semantics></math> mm.</p>
Full article ">Figure 10
<p>Reconstructed images of a 1-D capillary source using LASiP 8 at room-SiPM noise. The dashed lines mark the edges of the LASiPs.</p>
Full article ">Figure 11
<p>Reconstructed image of the simile-Derenzo source at different noise levels using LASiP 3. The dashed lines mark the edges of the LASiPs.</p>
Full article ">
14 pages, 501 KiB  
Article
Microdosimetric Simulation of Gold-Nanoparticle-Enhanced Radiotherapy
by Maxim Azarkin, Martin Kirakosyan and Vladimir Ryabov
Int. J. Mol. Sci. 2024, 25(17), 9525; https://doi.org/10.3390/ijms25179525 - 2 Sep 2024
Cited by 1 | Viewed by 1140
Abstract
Conventional X-ray therapy (XRT) is commonly applied to suppress cancerous tumors; however, it often inflicts collateral damage to nearby healthy tissue. In order to provide a better conformity of the dose distribution in the irradiated tumor, proton therapy (PT) is increasingly being used [...] Read more.
Conventional X-ray therapy (XRT) is commonly applied to suppress cancerous tumors; however, it often inflicts collateral damage to nearby healthy tissue. In order to provide a better conformity of the dose distribution in the irradiated tumor, proton therapy (PT) is increasingly being used to treat solid tumors. Furthermore, radiosensitization with gold nanoparticles (GNPs) has been extensively studied to increase the therapeutic ratio. The mechanism of radiosensitization is assumed to be connected to an enhancement of the absorbed dose due to huge photoelectric cross-sections with gold. Nevertheless, numerous theoretical studies, mostly based on Monte Carlo (MC) simulations, did not provide a consistent and thorough picture of dose enhancement and, therefore, the radiosensitization effect. Radiosensitization by nanoparticles in PT is even less studied than in XRT. Therefore, we investigate the physics picture of GNP-enhanced RT using an MC simulation with Geant4 equipped with the most recent physics models, taking into account a wide range of physics processes relevant for realistic PT and XRT. Namely, we measured dose enhancement factors in the vicinity of GNP, with diameters ranging from 10 nm to 80 nm. The dose enhancement in the vicinity of GNP reaches high values for XRT, while it is very modest for PT. The macroscopic dose enhancement factors for realistic therapeutic GNP concentrations are rather low for all RT scenarios; therefore, other physico-chemical and biological mechanisms should be additionally invoked for an explanation of the radiosensitization effect observed in many experiments. Full article
(This article belongs to the Special Issue Nanoparticles in Nanobiotechnology and Nanomedicine: 2nd Edition)
Show Figures

Figure 1

Figure 1
<p>Spatial energy density enhancement factor (<math display="inline"><semantics> <mrow> <mi>D</mi> <mi>E</mi> <msub> <mi>F</mi> <mi>SE</mi> </msub> </mrow> </semantics></math>) as a function of distance from the center of the gold nanoparticle (GNP) immersed in a homogeneous water system. The <math display="inline"><semantics> <mrow> <mi>D</mi> <mi>E</mi> <msub> <mi>F</mi> <mi>SE</mi> </msub> </mrow> </semantics></math>s are measured for GNPs with diameters of 10 nm (<b>a</b>), 20 nm (<b>b</b>), 40 nm (<b>c</b>), and 80 nm (<b>d</b>). The vertical red line marks the GNP surface.</p>
Full article ">Figure 2
<p>A schematic (not to scale) layout of the simulated system: (<b>a</b>) The macroscopic setup represents a cube consisting of human soft tissue. The dark orange layer represents a tumor, the smaller blue cubes represent the microscopic volume with a gold nanoparticle inside, either on frontal or distal parts of the tumor layer. (<b>b</b>) A close-up of a microscopic volume represents a water cube with gold nanoparticles of various (10 to 80 nm) radii inside. The thick dark red arrows denote the primary beam particles, and the thin black arrows indicate the secondary particles that can interact with the GNP as well.</p>
Full article ">Figure 3
<p>Relative dose versus tissue depth for different types of radiation.</p>
Full article ">Figure 4
<p>Energy distribution of 140 kVp (<b>a</b>), 6 MV (<b>b</b>), and photons and protons (<b>c</b>) at different depths in tissue. These distributions are normalized by the number of initial beam photons (<math display="inline"><semantics> <msubsup> <mi>N</mi> <mrow> <mi>γ</mi> </mrow> <mi>init</mi> </msubsup> </semantics></math>) and protons (<math display="inline"><semantics> <msubsup> <mi>N</mi> <mrow> <mi mathvariant="normal">p</mi> </mrow> <mi>init</mi> </msubsup> </semantics></math>).</p>
Full article ">
21 pages, 1820 KiB  
Article
Enhanced Particle Classification in Water Cherenkov Detectors Using Machine Learning: Modeling and Validation with Monte Carlo Simulation Datasets
by Ticiano Jorge Torres Peralta, Maria Graciela Molina, Hernan Asorey, Ivan Sidelnik, Antonio Juan Rubio-Montero, Sergio Dasso, Rafael Mayo-Garcia, Alvaro Taboada, Luis Otiniano and for the LAGO Collaboration
Atmosphere 2024, 15(9), 1039; https://doi.org/10.3390/atmos15091039 - 28 Aug 2024
Cited by 1 | Viewed by 1103
Abstract
The Latin American Giant Observatory (LAGO) is a ground-based extended cosmic rays observatory designed to study transient astrophysical events, the role of the atmosphere on the formation of secondary particles, and space-weather-related phenomena. With the use of a network of Water Cherenkov Detectors [...] Read more.
The Latin American Giant Observatory (LAGO) is a ground-based extended cosmic rays observatory designed to study transient astrophysical events, the role of the atmosphere on the formation of secondary particles, and space-weather-related phenomena. With the use of a network of Water Cherenkov Detectors (WCDs), LAGO measures the secondary particle flux, a consequence of the interaction of astroparticles impinging on the atmosphere of Earth. This flux can be grouped into three distinct basic constituents: electromagnetic, muonic, and hadronic components. When a particle enters a WCD, it generates a measurable signal characterized by unique features correlating to the particle’s type and the detector’s specific response. The resulting charge histograms from these signals provide valuable insights into the flux of primary astroparticles and their key characteristics. However, these data are insufficient to effectively distinguish between the contributions of different secondary particles. In this work, we extend our previous research by using detailed simulations of the expected atmospheric response to the primary flux and the corresponding response of our WCDs to atmospheric radiation. This dataset, which was created through the combination of the outputs of the ARTI and Meiga simulation frameworks, simulated the expected WCD signals produced by the flux of secondary particles during one day at the LAGO site in Bariloche, Argentina, situated at 865 m above sea level. This was achieved by analyzing the real-time magnetospheric and local atmospheric conditions for February and March of 2012, where the resultant atmospheric secondary-particle flux was integrated into a specific Meiga application featuring a comprehensive Geant4 model of the WCD at this LAGO location. The final output was modified for effective integration into our machine-learning pipeline. With an implementation of Ordering Points to Identify the Clustering Structure (OPTICS), a density-based clustering algorithm used to identify patterns in data collected by a single WCD, we have further refined our approach to implement a method that categorizes particle groups using advanced unsupervised machine learning techniques. This allowed for the differentiation among particle types and utilized the detector’s nuanced response to each, thus pinpointing the principal contributors within each group. Our analysis has demonstrated that applying our enhanced methodology can accurately identify the originating particles with a high degree of confidence on a single-pulse basis, highlighting its precision and reliability. These promising results suggest the feasibility of future implementations of machine-leaning-based models throughout LAGO’s distributed detection network and other astroparticle observatories for semi-automated, onboard and real-time data analysis. Full article
Show Figures

Figure 1

Figure 1
<p>The diagonal shows the overall shape of the distribution of each PCA component while, the lower diagonal graphs show 2D projections between the different components obtained using PCA for the selected features in [<a href="#B18-atmosphere-15-01039" class="html-bibr">18</a>]. By visual exploration, it different ‘groups’ can be observed, aggregated within other groups of varying density, showing the complexity of the data.</p>
Full article ">Figure 2
<p>Construction of a reachability plot. Cluster 1 can be observed within the reachability plot as a valley (see the marked blue arrow). The algorithm can distinguish clusters with different densities of points in different hierarchies. Figure adapted from Wang et al., 2019 [<a href="#B43-atmosphere-15-01039" class="html-bibr">43</a>].</p>
Full article ">Figure 3
<p>General data pipeline for the ML modeling.</p>
Full article ">Figure 4
<p>Cross-correlation matrix of initial feature set. A darker shade of orange corresponds to a higher value for the cross-correlation between two features.</p>
Full article ">Figure 5
<p>Example of a reachability plot from one run of the ML pipeline. The x-axis shows every point in the dataset ordered such as from smallest to greatest reachability distance, <math display="inline"><semantics> <msub> <mi>ϵ</mi> <mi>r</mi> </msub> </semantics></math>, with respect to the point’s closest core group. The y-axis is the <math display="inline"><semantics> <msub> <mi>ϵ</mi> <mi>r</mi> </msub> </semantics></math> value. A consistent cut-off threshold of 0.08 was used in every run where each independent run produced a similar plot. Points belonging to a cluster are colored and labeled accordingly, while any point in black did not gain membership to any cluster. A total of eight clusters were found.</p>
Full article ">Figure 6
<p>Histogram of charge, in black, with resulting clusters labeled and colored. The Y-axis is in a logarithmic scale for clarity.</p>
Full article ">Figure 7
<p>Stacked bar chart showing the percentage particle compositions for each cluster. It can be seen that the algorithm is very accurate in grouping muonic contributions.</p>
Full article ">Figure 8
<p>Zoomed in version of the histogram of charge, in black, highlighting clusters 0, 1 and 2.</p>
Full article ">Figure 9
<p>Zoomed-in version of the histogram of charge, in black. Clusters have been combined into three groups: group 1 is clusters 0, 1, and 2; group 2 is clusters 3 and 4; while group 3 is clusters 5, 6, and 7.</p>
Full article ">
Back to TopTop