Validation of mammographic x-ray spectra generated using Particle and Heavy Ion Transport code System

A Monte Carlo (MC) code is a robust method to generate a mammographic x-ray spectrum because the geometry of a mammography system can be flexible and directly modeled in MC simulation. However, simulations from MC code need to be validated before it can be reliably used for specific applications. This study aimed to generate and validate the x-ray spectra of relevant anodes used in mammography and breast tomosynthesis using Particle and Heavy Ion Transport code System (PHITS). PHITS version 3.08 was used to generate the x-ray spectra of molybdenum (Mo), rhodium (Rh), and tungsten (W) anodes. The Mo anode spectrum derived using PHITS was compared with those obtained using other MC codes. The generated spectra of all anodes were compared with the literature. Parameters including spectral shape, K characteristic x-ray yield, heel effect, and half-value layer (HVL) were used for a comparative assessment. The differences in these assessment parameters conducted by PHITS and PHITSEGS5 simulations were studied. Regarding the comparative parameters, PHITSEGS5 simulation improved the accuracy of mammographic x-ray generation compared to PHITS simulation; K x-ray and bremsstrahlung yields of the Mo anode spectrum generated by PHITSEGS5 simulation were a better agreement with those generated by other MC code simulations. The PHITSEGS5 spectra overestimated K x-ray and low-energy bremsstrahlung photons in comparison with measured spectra. Subsequently, HVLs calculated from PHITSEGS5 spectra were 1.0% (Mo/Mo) and 7.0% (W/Al) lower than those derived from measured spectra. For Mo and Rh anodes, relative difference of HVLs calculated from PHITSEGS5 spectra and those obtained from literature and measurement were within the TRS 457 acceptance criteria (+/- 0.02 mm Al). The observed difference exceeded the acceptance criteria for W anode. Regarding existed consistency in HVL between simulation and measurement, PHITSEGS5 simulation can be reliably used to generate x-ray spectra of Mo and Rh anodes. However, its accuracy should be improved for generating W anode spectrum.