KAHIBARO
8.4. Detector Materials
Table of Contents
LYSO
LYSO, or lutetium yttrium oxyorthosilicate, is one of the most common scintillator materials in modern PET scanners and in many gamma imaging detectors. In GATE, LYSO is typically defined as a compound material with elements lutetium, yttrium, silicon, and oxygen, together with a specific density and fractional composition that match published data or manufacturer specifications.
A typical LYSO composition in medical imaging simulations contains a large fraction of lutetium. This is important, because natural lutetium contains the radioisotope Lu 176, which is slightly radioactive and produces an intrinsic background in LYSO crystals. When you simulate PET detectors with LYSO, you may or may not want to include this intrinsic radioactivity, depending on your study. For many basic imaging or geometry studies it can be ignored, but for low activity or background studies it can be relevant.
In practice, you normally do not write a completely new LYSO definition by hand. The Geant4 material database or GATE example material files often provide a standard LYSO material that you can reuse. If you do define it yourself, you must ensure:
- Correct element fractions, usually given as mass fractions or number of atoms per formula unit.
- A realistic density, commonly around 7.1 g/cm$^3$.
- Consistency with the rest of your simulation, for example the same LYSO properties for all detector crystals.
When you are interested only in interaction and energy deposition, the elemental composition and density are enough. If you want to simulate optical photons later, you will also need optical properties, such as the refractive index and the light yield. These are usually attached as additional properties to the same LYSO material.
LYSO has a relatively high effective atomic number and a high density, so it has a good stopping power for 511 keV photons. This is one reason why it is widely used in PET. In GATE, using LYSO instead of a lower density material will clearly change the detection efficiency and the ratio between photoelectric and Compton interactions inside the crystal. You can verify this by comparing energy spectra from simulations with different crystal materials.
GATE does not force a single canonical definition for LYSO. For reproducible work you should always document the exact density and composition that you used. If you plan to compare with another study, try to match their definition as closely as possible.
For realistic PET simulations with LYSO, always define a density close to 7.1 g/cm$^3$ and use a consistent composition throughout all detector crystals.
BGO
BGO, bismuth germanate (Bi$_4$Ge$_3$O$_{12}$), is another classical scintillator for gamma and PET detectors. Compared to LYSO it has higher density but slower scintillation. BGO is still useful in simulations that focus on high stopping power or when you want to reproduce older or specific detector systems that used BGO crystals.
In GATE, BGO is typically created as a compound containing bismuth, germanium, and oxygen in the correct stoichiometric ratio. The high atomic number of bismuth gives BGO excellent gamma stopping efficiency, which is especially relevant in thick PET detector rings or in high energy gamma experiments. A typical density is around 7.1 g/cm$^3$ for BGO as well, but you should check reference data.
As with LYSO, you can often reuse a pre defined BGO material from a material database or from example scripts. If you define BGO yourself, you specify each element, then build the compound with its formula and density. Because BGO has a high effective atomic number, it produces a relatively strong photopeak at 511 keV and at other diagnostic energies in simulated spectra. This can be useful when you want to study energy resolution or detector performance.
BGO scintillation is slower than LYSO, which matters for time of flight PET and very high count rate simulations. Basic energy deposition simulations in GATE do not care about this timing detail, because they work at the level of energy deposited per step. Only when you add optical photon simulation and realistic electronics models do the timing properties of BGO become important.
In medical imaging educational simulations, BGO is often used to illustrate how different detector materials change sensitivity and count rates. By keeping the geometry and sources identical and changing only the material between LYSO and BGO, you can compare detection efficiency and the relative height of photopeaks and Compton continua.
When comparing detector materials in GATE, change only the material definition, not the crystal size or geometry, if you want to isolate the effect of material properties such as density and effective atomic number.
NaI
NaI, sodium iodide, is a very common scintillator, especially in gamma cameras and single photon detectors. In medical imaging, NaI(Tl) crystals are standard in many SPECT systems. In GATE, NaI is usually modelled as NaI or NaI doped with thallium, written NaI_Tl or similar, depending on the material library. For basic interaction and energy deposition, the small concentration of thallium does not significantly change the bulk interaction properties, so NaI alone is often sufficient.
NaI has a lower density than LYSO or BGO, around 3.7 g/cm$^3$, and a moderately high atomic number due to iodine. It has good light yield and has historically been very important in nuclear medicine imaging. In simulations, NaI crystals are usually larger in thickness than LYSO crystals to achieve similar stopping power, because of the lower density.
To use NaI in GATE, you either select a built in material or define a compound out of sodium and iodine with the correct stoichiometry and density. If you want to mimic a commercial gamma camera, you also need to use realistic crystal thickness and crystal area, but these belong to the geometry part, not the material part.
NaI is hygroscopic, which influences real detector design but does not affect Monte Carlo interaction physics. In simulation you do not need to worry about encapsulation or hermetic sealing unless you want to model the detector housing explicitly. However, housing materials such as aluminum or stainless steel can introduce additional attenuation and scattering, so for detailed studies you may include them as separate volumes around the NaI crystal.
Because NaI is often associated with gamma cameras, its typical photon energy range is between about 100 keV and 300 keV. In this range the balance between photoelectric and Compton interactions in NaI differs from that in LYSO or BGO. Energy spectra from SPECT type simulations will show a different photopeak shape and Compton tail depending on your choice of NaI thickness and geometry.
To ensure reproducibility, you should record whether you used NaI or NaI(Tl) in your simulation report. Even if the macroscopic interaction cross sections are almost the same, the notation helps others understand that you intend the material to represent a real thallium doped scintillator crystal.
For SPECT and gamma camera simulations, NaI with density around 3.7 g/cm$^3$ and realistic crystal thickness is important to obtain correct detection efficiency and energy spectra.
CsI
CsI, cesium iodide, is another scintillator material used in some medical imaging detectors, small gamma probes, and flat panel detectors. CsI can be doped with thallium or sodium, but in many Monte Carlo use cases you treat it simply as CsI with a given density and stoichiometric composition.
In GATE, you define CsI as a compound of cesium and iodine, or you use a predefined CsI material from a library. Its density is typically around 4.5 g/cm$^3$, higher than NaI but lower than LYSO and BGO. Because of its relatively high effective atomic number and light yield, CsI can be interesting in simulations where you compare alternative scintillator materials or where you mimic specific commercial detectors.
CsI is sometimes used in structured or pixelated form, for example needle like crystals grown on a substrate, particularly in digital X ray imaging. For such detectors, both material and geometry are important. The material defines how photons interact, while the geometry defines the effective fill factor and spatial resolution. In a beginner level GATE simulation, you often approximate these complex structures by uniform CsI slabs or simple pixel arrays.
For energy deposition simulations, you define CsI with its composition and density and assign it to the detector volume. The result will be a different balance between absorbed and transmitted photons compared to NaI or LYSO. In the diagnostic X ray energy range, CsI has good absorption efficiency and can be used as the sensitive layer of flat panel detectors. If you run a CT or X ray simulation with CsI detectors, you should set realistic thicknesses so that the detector does not unrealistically absorb all photons or let almost all pass.
If you plan to extend your simulations to optical photon transport and photodetector response, CsI requires additional optical parameters similar to LYSO or NaI(Tl). These include refractive index, light yield, emission spectrum, and surface properties. The detector geometry chapter will cover how to attach these properties to the correct surfaces.
For many educational examples, CsI can act as a second scintillator material to compare with NaI. By changing only the material data and keeping detector thickness fixed, you can study how different densities and atomic numbers modify spatial resolution and detection efficiency.
When you model CsI detectors in GATE, always couple the chosen material definition with a realistic thickness, because the combination of density and thickness determines the overall detection efficiency for your X ray or gamma energy range.
Silicon
Silicon plays a different role from the scintillator crystals described above. It is not usually a bulk scintillator material in medical imaging, but it is a fundamental material for semiconductor detectors and photodetectors. In GATE you can use silicon to model solid state detectors directly or as part of photodetectors such as silicon photomultipliers.
Silicon as a material is usually available directly in the Geant4 material database as a pure element, so you do not need to create a custom compound. You simply assign silicon to a detector volume that represents, for example, a semiconductor pixel detector or a silicon diode array. Its density is about 2.33 g/cm$^3$, and its interaction properties with photons and charged particles are well defined.
In energy deposition simulations, silicon behaves differently from high Z scintillators. For gamma detection, its lower atomic number and density mean that interaction probabilities at typical SPECT or PET energies are lower. Silicon based gamma detectors are therefore often thinner and designed for specific energy ranges. On the other hand, silicon is excellent for charged particle detection, including electrons and protons in some beam monitoring devices, because of its good signal to noise properties and the ability to fabricate segmented structures.
In optical photon simulations, silicon is used to represent photodetectors that detect scintillation light. Here the primary interest is not gamma interaction in silicon, but the conversion of optical photons to electronic signals. This is modeled through quantum efficiency or detection efficiency properties attached to silicon based surfaces or volumes. The underlying material is still silicon, but most of the relevant behavior is captured in additional optical and detector response parameters.
If you are building a mixed detector in GATE, for example a scintillator crystal coupled to a silicon photodetector, you would typically define two main materials: one for the scintillator, such as LYSO or CsI, and one for the photodetector, silicon. The geometry defines how they are coupled, and the physics and digitizer configuration define how interactions in each part are recorded and converted into hits and signals.
Because silicon is an elemental material with broad use throughout physics, you should avoid redefining it from scratch. Rely on the standard material definition provided by Geant4 or by your GATE distribution, so that its interaction cross sections and density remain consistent with other simulations and with published data.
Use the standard Geant4 silicon material for detector volumes, and only add extra properties, such as optical or electronic response, instead of redefining silicon itself.
Views: 11