Evaluation of the use of six diagnostic X-ray spectra computer codes

Characterization of backscatter factors for various tissue substitutes in diagnostic radiology: a Monte Carlo investigation

Accurate assessment of backscatter factors (BSFs) is critical in medical dosimetry to precisely quantify the increase in surface dose caused by photon scattering, particularly in the low-energy kilovoltage X-ray beams used in diagnostic radiology. This study aimed to conduct a comprehensive evaluation of BSF values for diagnostic X-ray beams through Monte Carlo simulations. The interactions of BSFs with widely used tissue substitutes, including water, ICRU tissue, polyester, polymethyl methacrylate (PMMA), and nylon, were examined across a range of conditions, including half-value layer (HVL), field size, and energy spectra. The results demonstrate that BSF values consistently increase with larger field sizes and higher beam energies/HVLs, highlighting the significant impact of these parameters on scatter contributions. Comparative analysis of the materials revealed that water most closely approximates the BSF behaviour of ICRU tissue, with deviations of -2.08-8% across the studied energy range and field sizes. Polyester and PMMA also showed promising agreement, converging to within ± 5% of ICRU tissue at higher energies and larger field sizes. In contrast, nylon exhibited more substantial deviations, particularly in smaller field sizes and lower energies. These findings provide essential insights to improve the accuracy of dosimetric models and enhance radiation safety in diagnostic radiology applications.

New opto-electro-mechanical sensor for two-dimensions dosimetry based on radiochromic films

This work presents the validation of a new Opto‒Electro-Mechanical (MOEM) system consisting of a matrix of photodetectors for two-dimensional dosimetry evaluation with radiochromic films. The proposed system is based on a 5 × 10 matrix of photodetectors controlled by both in-house electronic circuit and graphical user interface, which enables optical measurements directly. We present the first tests performed in an X-ray machine and 137Cs source with that array by using Gafchromic EBT3 films. We obtained similar results than with a standard method (e.g. flat-bed scanner). Results were compared with Monte Carlo simulations and very good agreement was found. Results show the feasibility of using this system for dose evaluations. To the best of our knowledge, this is the first MOEM sensor for radiotherapy. Further developments are ongoing to create an advanced 16 × 16 LDRs system covering 1.6 cm × 1.6 cm with a 1 mm of spatial resolution. We point to develop a portable dosimetry tool delivering dose maps in real time to improve the clinical application of radiochromic films.

Semiempirical, parameterized spectrum estimation for x-ray computed tomography

PURPOSE

To develop a tool to produce accurate, well-validated x-ray spectra for standalone use or for use in an open-access x-ray/CT simulation tool. Spectrum models will be developed for tube voltages in the range of 80 kVp through 140 kVp and for anode takeoff angles in the range of 5° to 9°.

METHODS

Spectra were initialized based on physics models, then refined using empirical measurements, as follows. A new spectrum-parameterization method was developed, including 13 spline knots to represent the bremsstrahlung component and 4 values to represent characteristic lines. Initial spectra at 80, 100, 120, and 140 kVp and at takeoff angles from 5° to 9° were produced using physics-based spectrum estimation tools XSPECT and SpekPy. Empirical experiments were systematically designed with careful selection of attenuator materials and thicknesses, and by reducing measurement contamination from scatter to <1%. Measurements were made on a 64-row CT scanner using the scanner's detector and using multiple layers of polymethylmethacrylate (PMMA), aluminum, titanium, tin, and neodymium. Measurements were made at 80, 100, 120, and 140 kVp and covering the entire 64-row detector (takeoff angles from 5° to 9°); a total of 6,144 unique measurements were made. After accounting for the detector's energy response, parameterized representations of the initial spectra were refined for best agreement with measurements using two proposed optimization schemes: based on modulation and based on gradient descent. X-ray transmission errors were computed for measurements vs calculations using the nonoptimized and optimized spectra. Half-value, tenth-value, and hundredth-value layers for PMMA, Al, and Ti were calculated.

RESULTS

Spectra before and after parameterization were in excellent agreement (e.g., R2 values of 0.995 and 0.997). Empirical measurements produced smoothly varying curves with x-ray transmission covering a range of up to 3.5 orders of magnitude. Spectra from the two optimization schemes, compared with the unoptimized physic-based spectra, each improved agreement with measurements by twofold through tenfold, for both postlog transmission data and for fractional value layers.

CONCLUSION

The resulting well-validated spectra are appropriate for use in the open-access x-ray/CT simulator under development, the x-ray-based Cancer Imaging Toolkit (XCIST), or for standalone use. These spectra can be readily interpolated to produce spectra at arbitrary kVps over the range of 80 to 140 kVp and arbitrary takeoff angles over the range of 5° to 9°. Furthermore, interpolated spectra over these ranges can be obtained by applying the standalone Matlab function available at https://github.com/xcist/documentation/blob/master/XCISTspectrum.m.

Dosimetric assessment of the exposure of radiotherapy patients due to cone-beam CT procedures

Cone-beam computed tomography (CBCT) is widely used for pre-treatment verification and patient setup in image-guided radiation therapy (IGRT). CBCT imaging is employed daily and several times per patient, resulting in potentially high cumulative imaging doses to healthy tissues that surround exposed target organs. Computed tomography dose index (CTDI) is the parameter used by CBCT equipment as indication of the radiation output to patients. This study aimed to increase the knowledge on the relation between CBCT organ doses and weighted CTDI (CTDI_W) for a thorax scanning protocol. A CBCT system was modelled using the Monte Carlo (MC) radiation transport program MCNPX2.7.0. Simulation results were validated against half-value layer (HVL), axial beam profile, patient skin dose (PSD) and CTDI measurements. For organ dose calculations, a male voxel phantom ("Golem") was implemented with the CBCT scanner computational model. After a successful MC model validation with measurements, a systematic comparison was performed between organ doses (and their distribution) and CTDI dosimetry concepts [CTDI_W and cumulative dose quantities f₁₀₀(150) and [Formula: see text]]. The results obtained show that CBCT organ doses vary between 1.2 ± 0.1 mGy and 3.3 ± 0.2 mGy for organs located within the primary beam. It was also verified that CTDI_W allows prediction of absorbed doses to tissues at distances of about 5 cm from the isocentre of the CBCT system, whereas f₁₀₀(150) allows prediction of organ doses at distances of about 10 cm from the isocentre, independently from its location. This study demonstrates that these dosimetric concepts are suitable methods that easily allow a good approximation of the additional CBCT imaging doses during a typical lung cancer IGRT treatment.

Estimating the absolute flux distribution for a synchrotron X-ray beam using ionization-chamber measurements with various filters

It is shown that an extensive set of accurate ionization-chamber measurements with a primary polychromatic synchrotron X-ray beam transmitted through various filter combinations/thicknesses can be used to quite effectively estimate the absolute flux distribution. The basic technique is simple but the `inversion' of the raw data to extract the flux distribution is a fundamentally ill-posed problem. It is demonstrated, using data collected at the Imaging and Medical Beamline (IMBL) of the Australian Synchrotron, that the absolute flux can be quickly and reliably estimated if a suitable choice of filters is made. Results are presented as a function of the magnetic field (from 1.40 to 4.00 T) of the superconducting multi-pole wiggler insertion device installed at IMBL. A non-linear least-squares refinement of the data is used to estimate the incident flux distribution and then comparison is made with calculations from the programs SPECTRA, XOP and spec.exe. The technique described is important not only in estimating flux itself but also for a variety of other, derived, X-ray properties such as beam quality, power density and absorbed-dose rate. The applicability of the technique with a monochromatic X-ray beam for which there is significant harmonic contamination is also demonstrated. Whilst absolute results can also be derived in this monochromatic beam case, relative (integrated) flux values are sufficient for our primary aim of establishing reliable determinations of the percentages of the various harmonic components.

A measurement-based generalized source model for Monte Carlo dose simulations of CT scans

The goal of this study is to develop a generalized source model for accurate Monte Carlo dose simulations of CT scans based solely on the measurement data without a priori knowledge of scanner specifications. The proposed generalized source model consists of an extended circular source located at x-ray target level with its energy spectrum, source distribution and fluence distribution derived from a set of measurement data conveniently available in the clinic. Specifically, the central axis percent depth dose (PDD) curves measured in water and the cone output factors measured in air were used to derive the energy spectrum and the source distribution respectively with a Levenberg-Marquardt algorithm. The in-air film measurement of fan-beam dose profiles at fixed gantry was back-projected to generate the fluence distribution of the source model. A benchmarked Monte Carlo user code was used to simulate the dose distributions in water with the developed source model as beam input. The feasibility and accuracy of the proposed source model was tested on a GE LightSpeed and a Philips Brilliance Big Bore multi-detector CT (MDCT) scanners available in our clinic. In general, the Monte Carlo simulations of the PDDs in water and dose profiles along lateral and longitudinal directions agreed with the measurements within 4%/1 mm for both CT scanners. The absolute dose comparison using two CTDI phantoms (16 cm and 32 cm in diameters) indicated a better than 5% agreement between the Monte Carlo-simulated and the ion chamber-measured doses at a variety of locations for the two scanners. Overall, this study demonstrated that a generalized source model can be constructed based only on a set of measurement data and used for accurate Monte Carlo dose simulations of patients' CT scans, which would facilitate patient-specific CT organ dose estimation and cancer risk management in the diagnostic and therapeutic radiology.

Evaluation of radiation dose to organs during kilovoltage cone-beam computed tomography using Monte Carlo simulation

Image-guided techniques for radiation therapy have improved the precision of radiation delivery by sparing normal tissues. Cone-beam computed tomography (CBCT) has emerged as a key technique for patient positioning and target localization in radiotherapy. Here, we investigated the imaging radiation dose delivered to radiosensitive organs of a patient during CBCT scan. The 4D extended cardiac-torso (XCAT) phantom and Geant4 Application for Tomographic Emission (GATE) Monte Carlo (MC) simulation tool were used for the study. A computed tomography dose index (CTDI) standard polymethyl methacrylate (PMMA) phantom was used to validate the MC-based dosimetric evaluation. We implemented an MC model of a clinical on-board imager integrated with the Trilogy accelerator. The MC model's accuracy was validated by comparing its weighted CTDI (CTDIw) values with those of previous studies, which revealed good agreement. We calculated the absorbed doses of various human organs at different treatment sites such as the head-and-neck, chest, abdomen, and pelvis regions, in both standard CBCT scan mode (125 kVp, 80 mA, and 25 ms) and low-dose scan mode (125 kVp, 40 mA, and 10 ms). In the former mode, the average absorbed doses of the organs in the head and neck and chest regions ranged 4.09-8.28 cGy, whereas those of the organs in the abdomen and pelvis regions were 4.30-7.48 cGy. In the latter mode, the absorbed doses of the organs in the head and neck and chest regions ranged 1.61-1.89 cGy, whereas those of the organs in the abdomen and pelvis region ranged between 0.79-1.85 cGy. The reduction in the radiation dose in the low-dose mode compared to the standard mode was about 20%, which is in good agreement with previous reports. We opine that the findings of this study would significantly facilitate decisions regarding the administration of extra imaging doses to radiosensitive organs.

Experimental and Monte Carlo-simulated spectra of standard mammography-quality beams

A spectrometric study of standard mammography-quality beams by using experimental and Monte Carlo simulation methods was carried out in this work. The qualities of these beams are described according to the International Electrotechical Commission 61267 standard and the Technical Report Series 457 International Atomic Energy Agency report. Specifically, the non-attenuated RQR-M beam series was studied.

METHODS

A Si-PIN diode-based spectrometer and the PENELOPE Monte Carlo code (v. 2008F1) were used for experiments and simulations, respectively. In addition, an ionization chamber was used to determine the half-value layers (HVLs) of each beam quality. The measurements were done in the mammography dosimeter calibration setup of our laboratory, and the Monte Carlo simulations reproduced such conditions.

RESULTS

The relative differences between the HVLs calculated from experimental and simulated spectra were lower than 2.4% for all the beam qualities studied. These differences are 1.2% and 3.1% when comparing the HVLs calculated from the experimental and simulated spectra to those determined by using the ionization chamber, respectively. A semi-empirical relation was found to obtain the nominal tube potential from the effective tube potential.

CONCLUSION

According to our results, the mammography beams used in this work have energy spectra similar to clinical beams.

CT scanner x-ray spectrum estimation from transmission measurements

PURPOSE

In diagnostic CT imaging, multiple important applications depend on the knowledge of the x-ray spectrum, including Monte Carlo dose calculations and dual-energy material decomposition analysis. Due to the high photon flux involved, it is difficult to directly measure spectra from the x-ray tube of a CT scanner. One potential method for indirect measurement involves estimating the spectrum from transmission measurements. The expectation maximization (EM) method is an accurate and robust method to solve this problem. In this article, this method was evaluated in a commercial CT scanner.

METHODS

Two step-wedges (polycarbonate and aluminum) were used to produce different attenuation levels. Transmission measurements were performed on the scanner and the measured data from the scanner were exported to an external computer to calculate the spectra. The EM method was applied to solve the equations that represent the attenuation processes of polychromatic x-ray photons. Estimated spectra were compared to the spectra simulated using a software provided by the manufacturer of the scanner. To test the accuracy of the spectra, a verification experiment was performed using a phantom containing different depths of water. The measured transmission data were compared to the transmission values calculated using the estimated spectra.

RESULTS

Spectra of 80, 100, 120, and 140 kVp from a dual-source CT scanner were estimated. The estimated and simulated spectra were well matched. The differences of mean energies were less than 1 keV. In the verification experiment, the measured and calculated transmission values were in excellent agreement.

CONCLUSIONS

Spectrum estimation using transmission data and the EM method is a quantitatively accurate and robust technique to estimate the spectrum of a CT system. This method could benefit studies relying on accurate knowledge of the x-ray spectra from CT scanner.

Monte Carlo simulator of realistic x-ray beam for diagnostic applications

PURPOSE

Monte Carlo simulation is a very useful tool for radiotherapy and diagnostic radiology. Yet even with the latest PCs, simulation of photon spectra emitted by an x-ray tube is a time-consuming task, potentially reducing the possibility to obtain relevant data such as dose evaluations, simulation of geometric settings, or monitor detector efficiency. This study developed and validated a method to generate random numbers for realistic beams in terms of photon spectrum and intensity to simulate x-ray tubes via Monte Carlo algorithms.

METHODS

Starting from literature data, the most common semiempirical models of bremsstrahlung are analyzed and implemented, adjusting their formulation to describe a large irradiation area (i.e., large field of view) and to take account of the heel effect as in common practice during patient examinations.

RESULTS

Simulation results show that Birch and Marshall's model is the fastest and most accurate for the aims of this work. Correction of the geometric size of the beam and validation of the intensity variation (heel effect) yielded excellent results with differences between experimental and simulated data of less than 6%.

CONCLUSIONS

The results of validation and execution time showed that the tube simulator calculates the x-ray photons quickly and efficiently and is perfectly capable of considering all the phenomena occurring in a real beam (total filtration, focal spot size, and heel effect), so it can be used in a wide range of applications such as industry, medical physics, or quality assurance.

Quantifying the effect of anode surface roughness on diagnostic x-ray spectra using Monte Carlo simulation

PURPOSE

The accurate prediction of x-ray spectra under typical conditions encountered in clinical x-ray examination procedures and the assessment of factors influencing them has been a longstanding goal of the diagnostic radiology and medical physics communities. In this work, the influence of anode surface roughness on diagnostic x-ray spectra is evaluated using MCNP4C-based Monte Carlo simulations.

METHODS

An image-based modeling method was used to create realistic models from surface-cracked anodes. An in-house computer program was written to model the geometric pattern of cracks and irregularities from digital images of focal track surface in order to define the modeled anodes into MCNP input file. To consider average roughness and mean crack depth into the models, the surface of anodes was characterized by scanning electron microscopy and surface profilometry. It was found that the average roughness (Ra) in the most aged tube studied is about 50 pm. The correctness of MCNP4C in simulating diagnostic x-ray spectra was thoroughly verified by calling its Gaussian energy broadening card and comparing the simulated spectra with experimentally measured ones. The assessment of anode roughness involved the comparison of simulated spectra in deteriorated anodes with those simulated in perfectly plain anodes considered as reference. From these comparisons, the variations in output intensity, half value layer (HVL), heel effect, and patient dose were studied.

RESULTS

An intensity loss of 4.5% and 16.8% was predicted for anodes aged by 5 and 50 microm deep cracks (50 kVp, 6 degrees target angle, and 2.5 mm A1 total filtration). The variations in HVL were not significant as the spectra were not hardened by more than 2.5%; however, the trend for this variation was to increase with roughness. By deploying several point detector tallies along the anode-cathode direction and averaging exposure over them, it was found that for a 6 degrees anode, roughened by 50 microm deep cracks, the reduction in exposure is 14.9% and 13.1% for 70 and 120 kVp tube voltages, respectively. For the evaluation of patient dose, entrance skin radiation dose was calculated for typical chest x-ray examinations. It was shown that as anode roughness increases, patient entrance skin dose decreases averagely by a factor of 15%.

CONCLUSIONS

It was concluded that the anode surface roughness can have a non-negligible effect on output spectra in aged x-ray imaging tubes and its impact should be carefully considered in diagnostic x-ray imaging modalities.

Quantitative contrast-enhanced mammography for contrast medium kinetics studies

Quantitative contrast-enhanced mammography, based on a dual-energy approach, aims to extract quantitative and temporal information of the tumour enhancement after administration of iodinated vascular contrast media. Simulations using analytical expressions and optimization of critical parameters essential for the development of quantitative contrast-enhanced mammography are presented. The procedure has been experimentally evaluated using a tissue-equivalent phantom and an amorphous silicon active matrix flat panel imager. The x-ray beams were produced by a tungsten target tube and spectrally shaped using readily available materials. Measurement of iodine projected thickness in mg cm(-2) has been performed. The effect of beam hardening does not introduce nonlinearities in the measurement of iodine projected thickness for values of thicknesses found in clinical investigations. However, scattered radiation introduces significant deviations from slope equal to unity when compared with the actual iodine projected thickness. Scatter correction before the analysis of the dual-energy images provides accurate iodine projected thickness measurements. At 10% of the exposure used in clinical mammography, signal-to-noise ratios in excess of 5 were achieved for iodine projected thicknesses less than 3 mg cm(-2) within a 4 cm thick phantom. For the extraction of temporal information, a limited number of low-dose images were used with the phantom incorporating a flow of iodinated contrast medium. The results suggest that spatial and temporal information of iodinated contrast media can be used to indirectly measure the tumour microvessel density and determine its uptake and washout from breast tumours. The proposed method can significantly improve tumour detection in dense breasts. Its application to perform in situ x-ray biopsy and assessment of the oncolytic effect of anticancer agents is foreseeable.

The development, validation and application of a multi-detector CT (MDCT) scanner model for assessing organ doses to the pregnant patient and the fetus using Monte Carlo simulations

The latest multiple-detector technologies have further increased the popularity of x-ray CT as a diagnostic imaging modality. There is a continuing need to assess the potential radiation risk associated with such rapidly evolving multi-detector CT (MDCT) modalities and scanning protocols. This need can be met by the use of CT source models that are integrated with patient computational phantoms for organ dose calculations. Based on this purpose, this work developed and validated an MDCT scanner using the Monte Carlo method, and meanwhile the pregnant patient phantoms were integrated into the MDCT scanner model for assessment of the dose to the fetus as well as doses to the organs or tissues of the pregnant patient phantom. A Monte Carlo code, MCNPX, was used to simulate the x-ray source including the energy spectrum, filter and scan trajectory. Detailed CT scanner components were specified using an iterative trial-and-error procedure for a GE LightSpeed CT scanner. The scanner model was validated by comparing simulated results against measured CTDI values and dose profiles reported in the literature. The source movement along the helical trajectory was simulated using the pitch of 0.9375 and 1.375, respectively. The validated scanner model was then integrated with phantoms of a pregnant patient in three different gestational periods to calculate organ doses. It was found that the dose to the fetus of the 3 month pregnant patient phantom was 0.13 mGy/100 mAs and 0.57 mGy/100 mAs from the chest and kidney scan, respectively. For the chest scan of the 6 month patient phantom and the 9 month patient phantom, the fetal doses were 0.21 mGy/100 mAs and 0.26 mGy/100 mAs, respectively. The paper also discusses how these fetal dose values can be used to evaluate imaging procedures and to assess risk using recommendations of the report from AAPM Task Group 36. This work demonstrates the ability of modeling and validating an MDCT scanner by the Monte Carlo method, as well as assessing fetal and organ doses by combining the MDCT scanner model and the pregnant patient phantom.

Estimation of patients' organ doses and conceptus doses from selected X-ray examinations in two Nigeria X-ray centres

In this study, organ and conceptus doses of patients undergoing chest, abdomen and skull radiograph examinations at two Nigeria X-ray centres, Niger State General Hospital (NGH) and Two-Tees (TTX), are reported. Air kerma was measured, and entrance surface dose (ESD) and half-value layer estimated for each set of tube potential (kV(p)), focus to skin distance and current-time product (mAs) used for each of the patients included in this study. Results show that the mean air kerma in the two centres are similar for the three projections considered in this study. Organ doses ranged from <0.01 to 2.18 mGy in NGH and from <0.01 to 1.29 mGy in TTX for examinations of the abdomen, from <0.01 to 0.20 mGy in NGH and from <0.01 to 0.13 mGy in TTX for examinations of the skull and from <0.01 to 3.90 mGy in NGH and from <0.01 to 1.96 mGy in TTX for examinations of the chest. Generally, no significant difference is seen between the organ doses of male and female patients. In NGH, organ doses are generally greater than those from TTX for the three examinations. The mean ESDs for examinations of the chest postero-anterior, abdomen antero-posterior (AP) and skull AP are, respectively, 5.37, 6.28 and 4.24 mGy in NGH, and 5.82, 5.33 and 4.76 mGy in TTX. The ESDs reported in this study, except for examinations of the chest, are generally lower than comparable values published in the literature. Conceptus doses were also estimated for female patients using normalised published conceptus dose data for abdomen examinations. The estimated conceptus doses were >1 mGy even when the conceptus was located 12 cm below the surface of the abdomen.

Study of influence of buildup factor form on simulated radiographic image

The objective of the study presented in this paper is the analysis of influence of different buildup factor forms on a simulated radiographic image. Simulated radiographic images are obtained by means of the ray-tracing technique. Scattered photons are modelled using the generally accepted geometric progression form, linear form and tabulated data of buildup factors. Simulated images were compared to the reference results obtained by Monte Carlo calculation. The best agreement to Monte Carlo simulated images is achieved for the geometric progression form of buildup factor.

Calculation of x-ray spectra emerging from an x-ray tube. Part II. X-ray production and filtration in x-ray targets

A new approach to the calculation of the x-ray spectrum emerging from an x-ray tube is proposed. Theoretical results for the bremsstrahlung cross section appearing in the literature are summarized. Four different treatments of electron penetration, based on the work presented in Part I, are then used to generate bremsstrahlung spectra. These spectra are compared to experimental data at 50, 80 and 100 kVp tube potentials. The most sophisticated treatment of electron penetration was required to obtain good agreement. With this treatment both the National Institute of Standards and Technology bremsstrahlung cross sections, based on accurate partial wave calculations, and the Bethe-Heitler cross section [H. A. Bethe and W. Heitler, Proc R. Soc. London, Ser. A. 146, 83-112 (1934)] corrected by a modified Elwert factor [G. Elwert, Ann. Phys. (Leipzig) 426, 178-208 (1939)], provided good agreement to measured data. An approximate treatment of the characteristic spectrum is suggested. The dependencies of the bremsstrahlung and characteristic outputs of an x-ray tube on tube potential are compared to experimentally derived data for 70-140 kVp potentials. Agreement is to within a few percent of the total output over the entire range. The spectral predictions of the semiempirical models of Birch and Marshall [R. Birch and M. Marshall, Phys. Med. Biol. 24, 505-513 (1979)] (IPEM Report 78) and of Tucker et al. [D. M. Tucker, G. T. Barnes, and D. P. Chakraborty, Med. Phys. 18, 211-218 (1991).] are also assessed. The predictions of Tucker et al. are very close to the model developed here. The predictions of IPEM Report 78 are similar, but consistently harder for the range of tube potentials examined (50-100 kV). Unlike the semiempirical models, the model proposed here requires the introduction of no empirical and unphysical parameters in the differential bremsstrahlung cross section, bar an overall normalization factor which is close to unity.

Calculation of x-ray spectra emerging from an x-ray tube. Part I. electron penetration characteristics in x-ray targets

The penetration characteristics of electron beams into x-ray targets are investigated for incident electron kinetic energies in the range 50-150 keV. The frequency densities of electrons penetrating to a depth x in a target, with a fraction of initial kinetic energy, u, are calculated using Monte Carlo methods for beam energies of 50, 80, 100, 120 and 150 keV in a tungsten target. The frequency densities for 100 keV electrons in Al, Mo and Re targets are also calculated. A mixture of simple modeling with equations and interpolation from data is used to generalize the calculations in tungsten. Where possible, parameters derived from the Monte Carlo data are compared to experimental measurements. Previous electron transport approximations in the semiempirical models of other authors are discussed and related to this work. In particular, the crudity of the use of the Thomson-Whiddington law to describe electron penetration and energy loss is highlighted. The results presented here may be used towards calculating the target self-attenuation correction for bremsstrahlung photons emitted within a tungsten target.

Study of increased radiation when an x-ray tube is placed in a strong magnetic field

When a fixed anode x-ray tube is placed in a magnetic field (B) that is parallel to the anode-cathode axis, the x-ray exposure increases with increasing B. It was hypothesized that the increase was caused by backscattered electrons which were constrained by B and reaccelerated by the electric field onto the x-ray tube target. We performed computer simulations and physical experiments to study the behavior of the backscattered electrons in a magnetic field, and their effects on the radiation output, x-ray spectrum, and off-focal radiation. A Monte Carlo program (EGS4) was used to generate the combined energy and angular distribution of the backscattered electrons. The electron trajectories were traced and their landing locations back on the anode were calculated. Radiation emission from each point was modeled with published data (IPEM Report 78), and thus the exposure rate and x-ray spectrum with the contribution of backscattered electrons could be predicted. The point spread function for a pencil beam of electrons was generated and then convolved with the density map of primary electrons incident on the anode as simulated with a finite element program (Opera-3d, Vector Fields, UK). The total spatial distribution of x-ray emission could then be calculated. Simulations showed that for an x-ray tube working at 65 kV, about 54% of the electrons incident on the target were backscattered. In a magnetic field of 0.5 T, although the exposure would be increased by 33%, only a small fraction of the backscattered electrons landed within the focal spot area. The x-ray spectrum was slightly shifted to lower energies and the half value layer (HVL) was reduced by about 6%. Measurements of the exposure rate, half value layer and focal spot distribution were acquired as functions of B. Good agreement was observed between experimental data and simulation results. The wide spatial distribution of secondary x-ray emission can degrade the MTF of the x-ray system at low spatial frequencies for B < 0.5 T.

Generation of calibrated tungsten target x-ray spectra: modified TBC model

In spite of the recent advances in the experimental detection of x-ray spectra, theoretical or semi-empirical approaches for determining realistic x-ray spectra in the range of diagnostic energies are important tools for planning experiments, estimating radiation doses in patients, and formulating radiation shielding models. The TBC model is one of the most useful approaches since it allows for straightforward computer implementation, and it is able to accurately reproduce the spectra generated by tungsten target x-ray tubes. However, as originally presented, the TBC model fails in situations where the determination of x-ray spectra produced by an arbitrary waveform or the calculation of realistic values of air kerma for a specific x-ray system is desired. In the present work, the authors revisited the assumptions used in the original paper published by . They proposed a complementary formulation for taking into account the waveform and the representation of the calculated spectra in a dosimetric quantity. The performance of the proposed model was evaluated by comparing values of air kerma and first and second half value layers from calculated and measured spectra by using different voltages and filtrations. For the output, the difference between experimental and calculated data was better then 5.2%. First and second half value layers presented differences of 23.8% and 25.5% in the worst case. The performance of the model in accurately calculating these data was better for lower voltage values. Comparisons were also performed with spectral data measured using a CZT detector. Another test was performed by the evaluation of the model when considering a waveform distinct of a constant potential. In all cases the model results can be considered as a good representation of the measured data. The results from the modifications to the TBC model introduced in the present work reinforce the value of the TBC model for application of quantitative evaluations in radiation physics.

A fit method for the determination of inherent filtration with diagnostic x-ray units

A method for the determination of total inherent filtration for clinical x-ray units using attenuation curves was devised. A model for the calculation of x-ray spectra is used to calculate kerma values which are then adjusted to the experimental data in minimizing the sum of the squared relative differences in kerma using a modified simplex fit process. The model considers tube voltage, voltage ripple, anode angle and additional filters. Fit parameters are the thickness of an additional inherent Al filter and a general normalization factor. Nineteen sets of measurements including attenuation data for three tube voltages and five Al-filter settings each were obtained. Relative differences of experimental and calculated kerma using the data for the additional filter thickness are within a range of -7.6% to 6.4%. Quality curves, i.e. the relationship of additional filtration to HVL, are often used to determine filtration but the results show that standard quality curves do not reflect the variety of conditions encountered in practice. To relate the thickness of the additional filter to the condition of the anode surface, the data fits were also made using tungsten as the filter material. These fits gave an identical fit quality compared to aluminium with a tungsten filter thickness of 2.12-8.21 microm which is within the range of the additional absorbing layers determined for rough anodes.