Simulating Medical Imaging X-Ray Tubes with Various Parameters Using BEAMnrc Monte Carlo Software ()
1. Introduction
The spectrum distribution of photons produced by an x-ray tube is crucial for various medical operations, along with CT medical imaging [1], dosimetry [2], x-ray fluorescence [3], mammography [4], dental radiography [5], radiotherapy [6], and so on. The aim of our work is to analyze the x-ray tube photon production using BEAMnrc Monte Carlo software, because several other studies have been completed using other software as e.g., MCNP [7], GEANT4 [8].
The scope is to identify crucial parameters such as anode target materials [2], anode angles [7] and filtration materials [9] which, in combination with existing research with other Monte Carlo software packages, will allow us to determine a Medical Imaging tube-standard which can be utilized by other EGSnrc users providing, simultaneously, tube information not easily accessible. This is very important, because manufacturers and vendors do not provide construction and spectrum data for x-ray tubes at photon energies of 120 keV and 30 keV, making it, hence, difficult to uncover such characteristics.
Our intention is to demonstrate that BEAMnrc is capable in mimicking x-ray creation at low energies as good as other commercial software, with the benefit of its free distribution, capabilities and independent multilayers design.
Electrons emitted from a filament due to thermionic emission, in a vacuum tube are driven at high velocity and bombard an anode-target. Tungsten and gold are typically employed as target materials in medical imaging x-ray tubes because of properties such as their high atomic mass and melting point, which are significant for CT and basic radiography at energies between 80 - 150 keV [10].
The target-anode materials at energies between 10 - 30 keV, where most mammographic exams are conducted, employ mainly molybdenum and rhodium as target materials. A heated filament generates electrons, which are subsequently converted to x-ray photons through the emission of Bremsstrahlung and the creation of characteristic x-rays [11]. The anode’s inclination affects x-ray consumption inside the target and, thus, when the angle is increased, a higher number of photons can mix with primary beam as the anode’s self-consumption diminishes [12] [13].
To decrease radiation exposure, x-rays are created inside the anode material and various high Z-materials are employed as filters [14] [15]. Aluminium and copper are mostly used at energies of 80 - 150 keV because these preferentially filter out the low energy photons that contribute to the skin’s radiation exposure while, simultaneously, retain excellent picture contrast for clinical analysis. Because the energies in mammography are low, around 10 - 30 keV, and subsequently mammography has a potential of radiation damage and cancer [16], the use of filtration is extremely important. The most suitable metals that can be used as filters between 10 - 30 keV, are rhodium and molybdenum.
Even in complicated geometries, the Monte Carlo approach for simulating particle passage through matter, is the most precise method for x-ray spectrum production [17] [18] [19]. The transportation of electrons and photons within the target and filters, can be estimated by utilizing the Monte Carlo approach, through which, extensive information can be calculated regarding the components that contribute to the creation of the x-ray spectrum. BEAMnrc [20] [21] [22] [23], a part of the EGSnrcmp application code, is frequently used to simulate megavoltage clinical linear accelerators. However, via modifications and use of EGSnrcmp-layer modules, it can be used for simulating the generation of lower energy x-ray tubes, as effectively. Due to this, it is chosen as the appropriate code for this investigation. X-ray tube simulations utilizing the Monte Carlo approach have a significant cost in computational time, especially, when many bremsstrahlung photons are created by electrons in the keV energy range.
When running such simulations, the smallest decrease in calculation speed is crucial. The directional bremsstrahlung splitting approach is utilized in our procedures to reduce variance [24]. The necessity of applying variance reduction techniques to drastically reduce simulation time is emphasized. Very significant is also the investigation of the impact of the target’s angle in the spectra on all of the scenarios that we’ve chosen. Finally, it is illustrated how the various metals that act as filters, and, are placed right after the x-ray tube, can affect the spectral distributions.
After conducting research, we have managed to detect, the specifications of various materials that are used anodes, the operational angles of each employed tube and the utilized filtering materials. We wanted to examine x-ray tubes for CT, basic radiography and mammography, two distinct energies, e.g., 120 keV and 30 keV are chosen for that purpose. In each energy, two different materials for the anode are employed. The number of photons generated with and without variance reduction technique is calculated, showing how much computational time can be saved. Also as per previous studies it is shown that as the angle of the x-ray tube increases, the number of scored photons increases due to the anode’s self-absorption. Lastly, it is shown that in modern x-ray tube used for diagnostic purposes the use of filter reduces the low-energy bremsstrahlung photons that do not participate in the diagnostic process, lowering, hence, the dose in an examination. Various researchers [12] [13] chose other software for their Monte Carlo simulations and determine that the angle, anode material and filtration. These parameters play a very important role and needs to be studied further, in order to create more efficient x-rays tubes to use them in our everyday medical applications. It is hard to find via experimentation what is the best operational angle, the effect of operational angle in x-ray spectra, what are the suitable target materials and filter materials and what are the best combinations we can determine for the x-ray tube to operate and in what application, e.g. CT, dental radiology, mammography etc. That is the reason Monte Carlo can help us identify such parameters.
We contribute to the field through the open-source program BEAMnrc, which may be further used to conduct simulated experiments of modern x-ray tubes.
2. Materials and Methods
Modern medical imaging x-ray tube design has offered novel solutions to difficulties like rotating anode targets, improved circuit design, that provide power to the filament and anode-cathode circuits, improved tube housing, for more radiation shielding, etc. [25]. The method that we utilize focuses mostly on the creation of x-rays using various materials and characteristics within the production module, which for the BEAMnrc’s case, is the XTUBE module. The operational default values for each case (CT and mammography) are W-Au target materials, 22˚ operational angle and thin layers of 0.25 cm Al - 0.25 cm Al + 0.02 cm Cu to act as filters for the CT operation and 0.003 cm Mo - 0.003 cm Rh target materials, 17˚ operational angle and thin layers of Mo-Rh to act as filters. The above values were determined from past works and we chose these values as the research of x-ray tubes has been progressed over the years.
We are modelling a monoenergetic beam of electrons striking the anode at 120 keV and 30 keV respectively, simulating high quality DC voltage. Since x-ray tube filament-cathode relation is a computerized procedure that generates a single photon at the beginning of each Monte Carlo run, it is important to pre-determine the total number of runs, i.e., the total number of Monte Carlo generated x-ray photons. Towards this, we generate 40 million electron particles for the 120 keV tube and 400 million for the 30 keV tube. The anode-target materials of tungsten and gold with filters of aluminium and copper, at angles of 17˚, 22˚, 27˚ degrees [12] are modelled for the 120 keV x-ray tube [13]. For the 30 keV case, x-ray tube target materials are molybdenum and rhodium, with filters of rhodium and molybdenum [8], at angles of 12˚, 17˚, 22˚ degrees [13].
The energy transfer thresholds AE and AP are set at 512 keV and 1 keV, respectively, while the cut-off energies ECUT and PCUT are set to 660 keV and 150 keV. Because the utilized geometry is the x-ray tube component, we employ the BEAMnrc GUI environment’s components XTUBE and SLABS to describe our geometry. Then, several rectangular components are employed with varying thicknesses to configure the filters and the scoring plane. The thickness of the anode-target is set to 1 cm in the Z direction, underneath the XTUBE component we specify a 20 cm × 20 cm vacuum rectangle with 2 cm thickness on the Z axis, followed by a 20 cm × 20 cm air gap with 2 cm thickness on the Z axis, and finally, a 20 cm × 20 cm rectangle with the filters. Note that the thicknesses of the filters depend on the scenario we are exaggerating. We set the materials of the filters to air in the cases where filtering is air. Underneath the final slab component, the scoring plane is placed.
The world’s medium is set to AIR, the material in front of the XTUBE component and on the back is set to vacuum. The combination of materials for the target-anode, (W, Au, Mo, Rh) and the filters (Al, Al + Cu, Mo, Rh) [26] [27] can be found in Table 1. For some runs, we use the variance reduction technique and the exact same runs take place altering just the option of variance reduction to off. The corresponding results can be found in Table 2. The importance of the variance reduction technique (DBS) can be observed. It is very significant that the total simulation time is reduced, while, importantly, more particles hit and score in our phase space files.
Table 1. Target-anode, filter materials and target-anode angle combinations.
Table 2. Target-anode materials with and without DBS variance reduction technique.
The EGSnrc parameters that we use throughout this study are as follows. The algorithm EXACT is used for crossing the boundaries, the method PRESTA-II for the step of each electron particle, option on for the effects of the spin of the electron, the Koch-Motz distribution is utilized for the bremsstrahlung angular sampling (i.e. calculation of the emission angles of bremsstrahlung photons), the NRC cross sections database for bremsstrahlung is selected, bound Compton scattering is enabled, Koch-Motz pair angular sampling is selected, the pair cross sections of the NRC database are utilized, the sampling for angular photoelectrons is turned on, we take into account the scattering due to the Rayleigh phenomena, we account the atomic relaxations by switching the corresponding option to on and the photon cross-sections of the xcom database are chosen. The type of electron beam source is the parallel circular beam, i.e., 90 degrees’ incident angle between X-Y axis and zero angle between both X-Z and Y-Z axes. The beam source radius is set to 0.02 cm and the kinetic energy is 120 keV and 30 keV, as aforementioned. The value was chosen so that the focal spot of the beam produced to be the default option for the two energy cases. Over the year’s other research has been conducted for various sources for medical imaging examination [28], as for our research we are more focused on the typical model of the bombardment of high velocity electrons on a metal target and the production of characteristic x-rays.
The phase space files produced were from 2.2 Gigabytes to 7.5 Gigabytes based on the filtration substance and initial kinetic energies of the electrons. For the combination of x-ray tube at 120 keV electron beam energy, 40 × 106 electron particles were initially produced by the cathode towards the anode-target, and, for the 30 keV, 40 × 107 electron particles were used. Depending on the filters and the starting particles, each simulation took around 8 - 12 hours to complete so as to minimize the uncertainties. From all our phase space files we extract information using BEAMDP [29] application and present the generated spectral distributions via xmgrace software [30].
3. Results and Discussions
Our study is being divided into four processes, the first process is to measure the effect that the usage of directional bremsstrahlung splitting variance reduction technique has on the data and the calculation time. The second process is to validate our simulations of the four different anode-target materials used and evaluate the Spectral distributions produced with the theoretical x-ray characteristics of the materials. The third process is to evaluate the affect that the anode angle has on the Spectral distributions. The fourth and last process is to show the effect of filtering the x-rays produced on the x-ray spectra.
Firstly, we would like to emphasize to the negative effect that the lack of variance reduction techniques imposes. This is characteristically shown in Table 2, where the initial particles generate significantly lower scored particles both for the simulation of the materials Mo and Rh at 30 keV incident energy and the materials W and Au at 120 keV. It should be noted from Table 2, that the initial particles for materials W and Au are one order higher (40 × 107) compared to the ones of Mo and Rh materials. This was implemented because, for low energy x-ray beams, the bremsstrahlung phenomena in the electron beam hitting the anode, are significantly less probable and, in this manner, much more initial electrons are needed for the effective modelling of the bombardment of the anode-target material which generates the x-rays. The size of the phase space files was improved by using the DBS (Directional Bremsstrahlung Splitting) variance reduction technique, which also enables the much faster code execution and outcome scoring. For each case, we can clearly see that the scored particles and photons are increased by 2 - 3 orders of magnitude while, importantly, reducing the simulation by 20 - 30 times.
Next the approach of employing AIR as a filter is shown. In this way, the scored photons are simulated as being directly generated from the anode, at 22˚ for the 120 keV tube and at 17˚ for the 30 keV tube, as also, Adeli et al. [8] have implemented. However, in this publication the MCNP and Geant4 Monte Carlo software were used. As can be observed the BEAMnrc produce results comparable to MCNP and Geant4, for all anode-target materials and that is a very important finding. It implies, clearly, that the present model, the coding and the simulation are adequate and can be considered valid. The reader should emphasize here that both EGSnrcmp and BEAMnrc are codes well validated both for the high photon energy range and, most importantly, the low energy one [23] [24].
In Tables 3-6, the values of x-ray characteristic radiation for the Tungsten (W), Gold (Au), Molybdenum (Mo) and Rhodium (Rh) anode-target are presented from publication of Thomson et al. [31]. In Figures 1-4 the Spectral Distributions
Table 3. Tungsten (W) x-ray characteristic [31].
Table 4. Gold (Au) x-ray characteristic [31].
Table 5. Molybdenum (Mo) x-ray characteristic [31].
Table 6. Rhodium (Rh) x-ray characteristic [31].
Figure 1. Tungsten (W) spectral distribution.
Figure 2. Gold (Au) spectral distribution.
of each anode-target material, produced from BEAMnrc, are presented with the theoretical values of Tables 3-6. This is a very important finding because it
Figure 3. Molybdenum (Mo) spectral distribution.
Figure 4. Rhodium (Rh) spectral distribution.
implies, in an additional way, that our simulation model is valid and that the Monte Carlo software is verified. More analytically in Figure 1 we can see clearly the Ka1 = 0.057 MeV, Ka2 = 0.059 MeV, Kb1 = 0.067 MeV, Kb2 = 0.069 MeV and La,b,c = 0.012 MeV characteristic x-ray of the tungsten (W) target material. In Figure 2 we can see clearly the Ka1 = 0.068 MeV, Ka2 = 0.066 MeV, Kb1 = 0.077 MeV, La = 0.012 MeV, Lb = 0.011 MeV and Lc = 0.013 MeV characteristic x-ray of the gold (Au) anode-target material. The Ka1, Ka2, Kb1, Kb2 are the ones that contribute to the image analysis and photon detection process of a CT scanner. In Figure 3 we can see clearly the Ka = 0.017 MeV, Kb = 0.019 MeV and La,b,c = 0.002 MeV, characteristic x-ray. In Figure 4 we can clearly see the Ka = 0.020 MeV, Kb = 0.022 MeV and La,b,c = 0.002 MeV characteristic x-ray. The Ka and Kb are the characteristic x-rays that penetrate the soft tissue and contribute to the mammographic procedure.
Furthermore, similarly with the works of Mesbahi et al. [12] and Kim et al. [13], which report the effect of anode angle using MCNP, Figures 5-8 present comparable results, which are also derived by simulating a similar experimental design as the above publications. It may be observed that the anode angle has a noteworthy effect on the spectral distributions of each anode material. As shown in Figures 5-8, as the target’s angle increases, a higher number of photons can hit the scoring plane, and the anode-target self-absorption (filtration) decreases, this phenomenon is known as the heel effect and we can clearly see that in the spectral distribution produced. This effect is more intense in the spectral distribution of Figure 7 & Figure 8, for low energy photons, (30 keV photon energy). Based on the works of Verhaegen et al. [26] and Homolka et al. [9], Figures 9-12 present clearly that the use of filters can significantly reduce an important amount of bremsstrahlung photons from the low energy band, allowing, hence, the electronics to achieve a better image during an examination. By filtering, a very large
Figure 5. Tungsten (W) spectral distributions for angles 17˚ - 22˚ - 27˚.
Figure 6. Gold (Au) spectral distributions for angles 17˚ - 22˚ - 27˚.
Figure 7. Molybdenum (Mo) spectral distributions for angles 12˚ - 17˚ - 22˚.
amount of those photons will not reach the patient and will not produce extra radiation thus, the exposure and the dose of the patients, are significantly decreased. More analytically in Figure 9 we can see that the La,b,c characteristics has been cut from the spectral distribution with the use of Al and with an extra Cu layer even more bremsstrahlung photons have been cut off but the Ka1, Ka2, Kb1,
Figure 8. Rhodium (Rh) spectral distributions for angles 12˚ - 17˚ - 22˚.
Figure 9. Tungsten (W) with angle 22˚ spectral distributions using filters of Air, Al, Al + Cu.
Kb2, have sustained. Also in Figure 10 the same effect is taking place as the layers of Al and Cu have been applied the low energy bremsstrahlung photons don’t pass to the patient, but only the Ka1, Ka2, Kb1 are let through. The ideal situation would be for all the bremsstrahlung photons to cut off and only the K* characteristic
Figure 10. Gold (Au) with angle 22˚ spectral distributions using filters of Air, Al, Al + Cu.
Figure 11. Molybdenum (Mo) with angle 17˚ spectral distributions using filters of Air, Mo, Rh.
x-rays to pass. In Figure 11 we can observe that with the application of Mo filter a large amount of photons has been cut off and only the Ka and Kb characteristics are let through, in most mammographic examinations Mo target with Mo filter are being applied but in recent years the application of Rh filter is being used as
Figure 12. Rhodium (Rh) with angle 17˚ spectral distributions using filters of Air, Rh, Mo.
the intensity of the beam is lower but it lets other characteristics to be present. In Figure 12 we can observe the phenomena of lower energy photons to be cut-off letting only the Ka and Kb characteristic to pass, but in mammography only the Rh target with Rh filter is being used because if we apply Mo filtration the spectral distribution is being swift and the Mo filter is amplifying other characteristic x-ray production of the Rh target with increased intensity so this option is not ideal for the radiation exposure performance of the examination.
4. Conclusions
Because it’s difficult to find the exact features of commercial x-ray tubes for medical imaging purposes, and the existing technology is not easily accessible, designing x-ray tubes with Monte Carlo methods is challenging. In this paper, we manage to utilize a simple but flexible approach. We examine a wide range of topics regarding x-ray tubes. We investigate the use of variance reduction and demonstrate its importance in Monte Carlo simulation, as it yields even to 20 - 30 times faster running time, at least for a certain scenario. Secondly, we validate our model via the BEAMnrc application for all employed anode-target materials, by comparing outcomes to those of similar works using different Monte Carlo software. With all our combinations of anode-target material and angles, we manage to investigate the effects of the anode’s angle. We show that as the angle of the anode increases more photons are scored and this is a property connected to the heel effect and the self-absorption on the anode, we demonstrate these effects on the spectral distributions.
Lastly, we report that when photons are filtered in an examination, the dose of radiation received by a patient is significantly reduced. The findings obtained using our approach are comparable to the outcomes of recent studies.
After taking into consideration the theory and previous works on the subject of generation and creation of x-rays for medical imaging equipment we can support the view that BEAMnrc is a software that can be used to investigate efficiently and effectively x-ray tubes, as it is comparable to other commercial software, while it imposes great advantages.
Acknowledgements
The work of authors at the University of West Attica has been funded by The Special Account for Research Funds of ELKE/PADA.