Characterization of prompt gamma-ray emission with respect to the Bragg peak for proton beam range verification: A Monte Carlo study

In this paper we report a Geant4 simulation study to investigate the characteristic prompt gamma (PG) emission in a water phantom for real-time monitoring of the Bragg peak (BP) during proton beam irradiation. The PG production, emission spatial correlation with the BP, and position preference for detection with respect to the BP have been quantified in different PG energy windows as a function of proton pencil-beam energy from 100 to 200 MeV. The PG response to small BP shifts was evaluated using a 2 cm-thick slab with different human body materials embedded in a water phantom. Our results show that the prominent characteristic PG emissions of 4.44, 5.21 and 6.13 MeV exhibit distinctive correlation with the dose deposition curve. The accuracy in BP position identification using these characteristic PG rays is highly consistent as the beam energy increases from 100 to 200 MeV. There exists a position preference for PG detection with respect to the BP position, which has a strong dependence on the proton beam energy and PG energies. It was also observed that a submillimeter shift of the BP position can be realized by using PG signals. These results indicate that the characteristic PG signal is sensitive and reliable for BP tracking. Although the maximization of the PG measurement associated with the BP is difficult, it can be optimized with energy and detection position preferences.     

Influence of beam incidence and irradiation parameters on stray neutron doses to healthy organs of pediatric patients treated for an intracranial tumor with passive scattering proton therapy

PurposeIn scattering proton therapy, the beam incidence, i.e. the patient’s orientation with respect to the beam axis, can significantly influence stray neutron doses although it is almost not documented in the literature.MethodsMCNPX calculations were carried out to estimate stray neutron doses to 25 healthy organs of a 10-year-old female phantom treated for an intracranial tumor. Two beam incidences were considered in this article, namely a superior (SUP) field and a right lateral (RLAT) field. For both fields, a parametric study was performed varying proton beam energy, modulation width, collimator aperture and thickness, compensator thickness and air gap size.ResultsUsing a standard beam line configuration for a craniopharyngioma treatment, neutron absorbed doses per therapeutic dose of 63 μGy Gy⁻¹ and 149 μGy Gy⁻¹ were found at the heart for the SUP and the RLAT fields, respectively. This dose discrepancy was explained by the different patient’s orientations leading to changes in the distance between organs and the final collimator where external neutrons are mainly produced. Moreover, investigations on neutron spectral fluence at the heart showed that the number of neutrons was 2.5 times higher for the RLAT field compared against the SUP field. Finally, the influence of some irradiation parameters on neutron doses was found to be different according to the beam incidence.ConclusionBeam incidence was thus found to induce large variations in stray neutron doses, proving that this parameter could be optimized to enhance the radiation protection of the patient.     

A deep learning approach for converting prompt gamma images to proton dose distributions: A Monte Carlo simulation study

PurposeIn proton therapy, imaging prompt gamma (PG) rays has the potential to verify proton dose (PD) distribution. Despite the fact that there is a strong correlation between the gamma-ray emission and PD, they are still different in terms of the distribution and the Bragg peak (BP) position. In this work, we investigated the feasibility of using a deep learning approach to convert PG images to PD distributions.MethodsWe designed the Monte Carlo simulations using 20 digital brain phantoms irradiated with a 100-MeV proton pencil beam. Each phantom was used to simulate 200 pairs of PG images and PD distributions. A convolutional neural network based on the U-net architecture was trained to predict PD distributions from PG images.ResultsOur simulation results show that the pseudo PD distributions derived from the corresponding PG images agree well with the simulated ground truths. The mean of the BP position errors from each phantom was less than 0.4 mm. We also found that 2000 pairs of PG images and dose distributions would be sufficient to train the U-net. Moreover, the trained network could be deployed on the unseen data (i.e. different beam sizes, proton energies and real patient CT data).ConclusionsOur simulation study has shown the feasibility of predicting PD distributions from PG images using a deep learning approach, but the reliable prediction of PD distributions requires high-quality PG images. Image-degrading factors such as low counts and limited spatial resolution need to be considered in order to obtain high-quality PG images.     

Fiber background rejection and crystal over-response compensation for GaN based in vivo dosimetry

For dosimetric measurements using an implantable optical fiber probe with GaN (Gallium Nitride) scintillator as radioluminescence (RL) transducer, a bi-channel method is proposed to reject the background contribution of the irradiated fiber segment. It is based on spectral differences between the narrow-band light emission from GaN and the large-band background from the irradiated optical fiber. Experimental validation of this method using 6 MV photon beam has shown that the remaining background contribution after subtraction is below 1.2% for square field sizes ranging from 3 cm to 20 cm. Furthermore, a compensation method for the over-response of GaN is also proposed, since GaN is not tissue equivalent. The over-response factor of GaN exhibits a linear increase with square field aperture and depends on depth from phantom surface. This behaviour is modelled to allow compensation in specific conditions. The proposed method has been evaluated and has shown a maximum deviation of 3% for a 6 MV photon beam and 1% for an 18 MV photon beam at a depth beyond the build-up region.     

An inhomogeneous most likely path formalism for proton computed tomography

PurposeMultiple Coulomb scattering (MCS) poses a challenge in proton CT (pCT) image reconstruction. The assumption of straight paths is replaced with Bayesian models of the most likely path (MLP). Current MLP-based pCT reconstruction approaches assume a water scattering environment. We propose an MLP formalism based on accurate determination of scattering moments in inhomogeneous media.MethodsScattering power relative to water (RScP) was calculated for a range of human tissues and investigated against relative stopping power (RStP). Monte Carlo simulation was used to compare the new inhomogeneous MLP formalism to the water approach in a slab geometry and a human head phantom. An MLP-Spline-Hybrid method was investigated for improved computational efficiency.ResultsA piecewise-linear correlation between RStP and RScP was shown, which may assist in iterative pCT reconstruction. The inhomogeneous formalism predicted Monte Carlo proton paths through a water cube with thick bone inserts to within 1.0 mm for beams ranging from 210 to 230 MeV incident energy. Improvement in accuracy over the conventional MLP ranged from 5% for a 230 MeV beam to 17% for 210 MeV. There was no noticeable gain in accuracy when predicting 200 MeV proton paths through a clinically relevant human head phantom. The MLP-Spline-Hybrid method reduced computation time by half while suffering negligible loss of accuracy.ConclusionsWe have presented an MLP formalism that accounts for material composition. In most clinical cases a water scattering environment can be assumed, however in certain cases of significant heterogeneity the proposed algorithm may improve proton path estimation.     

The optimal balance between quality and efficiency in proton radiography imaging technique at various proton beam energies: A Monte Carlo study

Proton radiography is a novel imaging modality that allows direct measurement of the proton energy loss in various tissues. Currently, due to the conversion of so-called Hounsfield units from X-ray Computed Tomography (CT) into relative proton stopping powers (RPSP), the uncertainties of RPSP are 3–5% or higher, which need to be minimized down to 1% to make the proton treatment plans more accurate.In this work, we simulated a proton radiography system, with position-sensitive detectors (PSDs) and a residual energy detector (RED). The simulations were built using Geant4, a Monte Carlo simulation toolkit. A phantom, consisting of several materials was placed between the PSDs of various Water Equivalent Thicknesses (WET), corresponding to an ideal detector, a gaseous detector, silicon and plastic scintillator detectors. The energy loss radiograph and the scattering angle distributions of the protons were studied for proton beam energies of 150 MeV, 190 MeV and 230 MeV. To improve the image quality deteriorated by the multiple Coulomb scattering (MCS), protons with small angles were selected. Two ways of calculating a scattering angle were considered using the proton’s direction and position.A scattering angle cut of 8.7 mrad was applied giving an optimal balance between quality and efficiency of the radiographic image. For the three proton beam energies, the number of protons used in image reconstruction with the direction method was half the number of protons kept using the position method.     

Hot and cold contrasts in high-resolution Tc-99m planar scintigraphy: A survey of fifty-two camera heads using the PICKER thyroid phantom

This study was aimed at comparing the sensitivity and hot and cold contrasts obtained when imaging the PICKER thyroid phantom using gamma cameras fitted with either their ultra-high or high-resolution low-energy parallel hole collimator.Seventeen camera models from Elscint, General Electric, Siemens and Sopha Medical Vision were involved in the study for a total of 30 cameras and 52 camera heads. A single operator conducted the study in order to minimize the impact of human factors. The phantom contained about 74 MBq ^99mTc and was imaged at 10 cm from the collimator face with the energy window that are recommended by the camera manufacturer. A total of 1 million counts were accumulated.Hot and cold contrasts were in mean of about 0.05 higher when using an ultra-high-resolution than when using a high-resolution low-energy collimator. This higher contrast was obtained at the expense of a mean reduction in sensitivity of 30%. In particular, Elscint cameras demonstrated a 30% lower sensitivity whatever the collimator type. The Sopha Medical Vision DST and DSX cameras and the General Electric Magicam camera offered the lowest contrasts among the cameras with a high-resolution collimator. Although this was accompanied by a higher than the mean sensitivity for the DST and DSX, the Magicam demonstrated sensitivity roughly identical to the mean of all the cameras with a high-resolution collimator.     

Development and validation of a new set-up simulator dedicated to ocular proton therapy at CNAO

An Eye Tracking System (ETS) is used at CNAO for providing a stable and reproducible ocular proton therapy (OPT) set-up, featuring a fixation light (FL) and monitoring stereo-cameras embedded in a rigid case. The aim of this work is to propose an ETS set-up simulation algorithm, that automatically provides the FL positioning in space, according to patient-specific gaze direction and avoiding interferences with patient, beam and collimator.Two configurations are provided: one in the CT room for acquiring images required for treatment planning with the patient lying on a couch, and one related to the treatment room with the patient sitting in front of the beam. Algorithm validation was performed reproducing ETS simulation (CT) and treatment (room) set-up for 30 patients previously treated at CNAO. The positioning accuracy of the device was quantified through a set of 14 control points applied to the ETS case and localizable both in the CT volume and in room X-ray images.Differences between the position of ETS reference points estimated by the algorithm and those measured by imaging systems are reported. The corresponding gaze direction deviation is on average 0.2° polar and 0.3° azimuth for positioning in CT room and 0.1° polar and 0.4° azimuth in the treatment room.The simulation algorithm was embedded in a clinically usable software application, which we assessed as capable of ensuring ETS positioning with an average accuracy of 2 mm in CT room and 1.5 mm in treatment room, corresponding to gaze direction deviations consistently lower than 1°.     

GEANT4 models for the secondary radiation flux in the collimation system of a 300 MeV proton microbeam

In Harbin, we are developing a 300 MeV proton microbeam for many applications in space science including upset studies in microelectronic devices, radiation hardness of materials for satellites and radiation effects in human tissues. There are also applications of this facility proposed for proton therapy. The microbeam system will employ a purpose-built proton synchrotron to provide the beam. However there are many challenges to be addressed in the design, construction and operation of this facility. Here we address two important design aspects for which we apply GEANT4 modeling. First, the high energy proton beam interacts strongly with beam line materials, especially the collimation slits, to produce showers of secondary particles which could introduce significant background signals and degrade the resolution of the proton microbeam. Second, the beam transport within the residual vacuum of the beam line may also introduce undesirable background radiation. In both cases mitigation strategies need to be incorporated during the design phase of the new system. We study the use of a dipole magnet following the aperture collimator to reduce the flux of secondary particles incident on the analysis chamber. Monte Carlo simulations are performed using GEANT4 and SRIM. By inserting the dipole magnet, we find as expected a significant reduction in the scattering of protons and other particles, such as neutrons and gamma rays, at the collimation system exit position. Secondary radiation from the residual gas pressure within the beam line vacuum system are also modelled and found to be negligible under the standard operating conditions.     

Evaluations of a flat-panel based compact daily quality assurance device for proton pencil beam scanning (PBS) system

PurposeTo evaluate the flat-panel detector quenching effect and clinical usability of a flat-panel based compact QA device for PBS daily constancy measurements.Materials & MethodThe QA device, named Sphinx Compact, is composed of a 20x20 cm² flat-panel imager mounted on a portable frame with removable plastic modules for constancy checks of proton energy (100 MeV, 150 MeV, 200 MeV), Spread-Out-Bragg-Peak (SOBP) profile, and machine output. The potential quenching effect of the flat-panel detector was evaluated. Daily PBS QA tests of X-ray/proton isocenter coincidence, the constancy of proton spot position and sigma as well as the energy of pristine proton beam, and the flatness of SOBP proton beam through the 'transformed' profile were performed and analyzed. Furthermore, the sensitivity of detecting energy changes of pristine proton beam was also evaluated.ResultsThe quenching effect was observed at depths near the pristine peak regions. The flat-panel measured range of the distal 80% is within 0.9 mm to the defined ranges of the delivered proton beams. X-ray/proton isocenter coincidence tests demonstrated maximum mismatch of 0.3 mm between the two isocenters. The device can detect 0.1 mm change of spot position and 0.1 MeV energy changes of pristine proton beams. The measured transformed SOBP beam profile through the wedge module rendered as flat.ConclusionsEven though the flat-panel detector exhibited quenching effect at the Bragg peak region, the proton range can still be accurately measured. The device can fulfill the requirements of the daily QA tests recommended by the AAPM TG224 Report.     

Proton microbeam radiotherapy with scanned pencil-beams – Monte Carlo simulations

Irradiation, delivered by a synchrotron facility, using a set of highly collimated, narrow and parallel photon beams spaced by 1 mm or less, has been termed Microbeam Radiation Therapy (MRT). The tolerance of healthy tissue after MRT was found to be better than after standard broad X-ray beams, together with a more pronounced response of malignant tissue. The microbeam spacing and transverse peak-to-valley dose ratio (PVDR) are considered to be relevant biological MRT parameters. We investigated the MRT concept for proton microbeams, where we expected different depth-dose profiles and PVDR dependences, resulting in skin sparing and homogeneous dose distributions at larger beam depths, due to differences between interactions of proton and photon beams in tissue. Using the FLUKA Monte Carlo code we simulated PVDR distributions for differently spaced 0.1 mm (sigma) pencil-beams of entrance energies 60, 80, 100 and 120 MeV irradiating a cylindrical water phantom with and without a bone layer, representing human head. We calculated PVDR distributions and evaluated uniformity of target irradiation at distal beam ranges of 60–120 MeV microbeams. We also calculated PVDR distributions for a 60 MeV spread-out Bragg peak microbeam configuration. Application of optimised proton MRT in terms of spot size, pencil-beam distribution, entrance beam energy, multiport irradiation, combined with relevant radiobiological investigations, could pave the way for hypofractionation scenarios where tissue sparing at the entrance, better malignant tissue response and better dose conformity of target volume irradiation could be achieved, compared with present proton beam radiotherapy configurations.     

Incorrect dosimetric leaf separation in IMRT and VMAT treatment planning: Clinical impact and correlation with pretreatment quality assurance

PurposeDynamic treatment planning algorithms use a dosimetric leaf separation (DLS) parameter to model the multi-leaf collimator (MLC) characteristics. Here, we quantify the dosimetric impact of an incorrect DLS parameter and investigate whether common pretreatment quality assurance (QA) methods can detect this effect.Methods16 treatment plans with intensity modulated radiation therapy (IMRT) or volumetric modulated arc therapy (VMAT) technique for multiple treatment sites were calculated with a correct and incorrect setting of the DLS, corresponding to a MLC gap difference of 0.5 mm. Pretreatment verification QA was performed with a bi-planar diode array phantom and the electronic portal imaging device (EPID). Measurements were compared to the correct and incorrect planned doses using gamma evaluation with both global (G) and local (L) normalization. Correlation, specificity and sensitivity between the dose volume histogram (DVH) points for the planning target volume (PTV) and the gamma passing rates were calculated.ResultsThe change in PTV and organs at risk DVH parameters were 0.4–4.1%. Good correlation (>0.83) between the PTV_mean dose deviation and measured gamma passing rates was observed. Optimal gamma settings with 3%L/3 mm (per beam and composite plan) and 3%G/2 mm (composite plan) for the diode array phantom and 2%G/2 mm (composite plan) for the EPID system were found. Global normalization and per beam ROC analysis of the diode array phantom showed an area under the curve <0.6.ConclusionsA DLS error can worsen pretreatment QA using gamma analysis with reasonable credibility for the composite plan. A low detectability was demonstrated for a 3%G/3 mm per beam gamma setting.     

Correction of measured Gamma-Knife output factors for angular dependence of diode detectors and PinPoint ionization chamber

Dosimetry for Gamma-Knife requires detectors with high spatial resolution and minimal angular dependence of response. Angular dependence and end effect time for p-type silicon detectors (PTW Diode P and Diode E) and PTW PinPoint ionization chamber were measured with Gamma-Knife beams. Weighted angular dependence correction factors were calculated for each detector. The Gamma-Knife output factors were corrected for angular dependence and end effect time. For Gamma-Knife beams angle range of 84°–54°. Diode P shows considerable angular dependence of 9% and 8% for the 18 mm and 14, 8, 4 mm collimator, respectively. For Diode E this dependence is about 4% for all collimators. PinPoint ionization chamber shows angular dependence of less than 3% for 18, 14 and 8 mm helmet and 10% for 4 mm collimator due to volumetric averaging effect in a small photon beam. Corrected output factors for 14 mm helmet are in very good agreement (within ±0.3%) with published data and values recommended by vendor (Elekta AB, Stockholm, Sweden). For the 8 mm collimator diodes are still in good agreement with recommended values (within ±0.6%), while PinPoint gives 3% less value. For the 4 mm helmet Diodes P and E show over-response of 2.8% and 1.8%, respectively. For PinPoint chamber output factor of 4 mm collimator is 25% lower than Elekta value which is generally not consequence of angular dependence, but of volumetric averaging effect and lack of lateral electronic equilibrium. Diodes P and E represent good choice for Gamma-Knife dosimetry.     

Development of an x-ray-opaque-marker system for quantitative phantom positioning in patient-specific quality assurance

PurposeWe developed an x-ray-opaque-marker (XOM) system with inserted fiducial markers for patient-specific quality assurance (QA) in CyberKnife (Accuray) and a general-purpose linear accelerator (linac). The XOM system can be easily inserted or removed from the existing patient-specific QA phantom. Our study aimed to assess the utility of the XOM system by evaluating the recognition accuracy of the phantom position error and estimating the dose perturbation around a marker.MethodsThe recognition accuracy of the phantom position error was evaluated by comparing the known error values of the phantom position with the values measured by matching the images with target locating system (TLS; Accuray) and on-board imager (OBI; Varian). The dose perturbation was evaluated for 6 and 10 MV single-photon beams through experimental measurements and Monte Carlo simulations.ResultsThe root mean squares (RMSs) of the residual position errors for the recognition accuracy evaluation in translations were 0.07 mm with TLS and 0.30 mm with OBI, and those in rotations were 0.13° with TLS and 0.15° with OBI. The dose perturbation was observed within 1.5 mm for 6 MV and 2.0 mm for 10 MV from the marker.ConclusionsSufficient recognition accuracy of the phantom position error was achieved using our system. It is unnecessary to consider the dose perturbation in actual patient-specific QA. We concluded that the XOM system can be utilized to ensure quantitative and accurate phantom positioning in patient-specific QA with CyberKnife and a general-purpose linac.     

Development of Monte Carlo based real-time treatment planning system with fast calculation algorithm for boron neutron capture therapy

PurposeWe simulated the effect of patient displacement on organ doses in boron neutron capture therapy (BNCT). In addition, we developed a faster calculation algorithm (NCT high-speed) to simulate irradiation more efficiently.MethodsWe simulated dose evaluation for the standard irradiation position (reference position) using a head phantom. Cases were assumed where the patient body is shifted in lateral directions compared to the reference position, as well as in the direction away from the irradiation aperture.For three groups of neutron (thermal, epithermal, and fast), flux distribution using NCT high-speed with a voxelized homogeneous phantom was calculated. The three groups of neutron fluxes were calculated for the same conditions with Monte Carlo code. These calculated results were compared.ResultsIn the evaluations of body movements, there were no significant differences even with shifting up to 9 mm in the lateral directions. However, the dose decreased by about 10% with shifts of 9 mm in a direction away from the irradiation aperture.When comparing both calculations in the phantom surface up to 3 cm, the maximum differences between the fluxes calculated by NCT high-speed with those calculated by Monte Carlo code for thermal neutrons and epithermal neutrons were 10% and 18%, respectively. The time required for NCT high-speed code was about 1/10th compared to Monte Carlo calculation.ConclusionsIn the evaluation, the longitudinal displacement has a considerable effect on the organ doses.We also achieved faster calculation of depth distribution of thermal neutron flux using NCT high-speed calculation code.     

A novel method for dose distribution registration using fiducial marks made by a megavoltage beam in film dosimetry for intensity-modulated radiation therapy quality assurance

PurposePhotographic film is widely used for the dose distribution verification of intensity-modulated radiation therapy (IMRT). However, analysis for verification of the results is subjective. We present a novel method for marking the isocenter using irradiation from a megavoltage (MV) beam transmitted through slits in a multi-leaf collimator (MLC).MethodsWe evaluated the effect of the marking irradiation at 500 monitor units (MU) on the total transmission through the MLC using an ionization chamber and Radiochromic Film. Film dosimetry was performed for quality assurance (QA) of IMRT plans. Three methods of registration were used for each film: marking by irradiating with an MV beam through slits in the MLC (MLC-IC); marking with a fabricated phantom (Phantom-IC); and a subjective method based on isodose lines (Manual). Each method was subjected to local γ-analysis.ResultsThe effect of the marking irradiation on the total transmission was 0.16%, as measured by a ionization chamber at a 10-cm depth in a solid phantom, while the inter-leaf transmission was 0.3%, determined from the film. The mean pass rates for each registration method agreed within ±1% when the criteria used were a distance-to-agreement (DTA) of 3 mm and a dose difference (DD) of 3%. For DTA/DD criteria of 2 mm/3%, the pass rates in the sagittal plane were 96.09 ± 0.631% (MLC-IC), 96.27 ± 0.399% (Phantom-IC), and 95.62 ± 0.988% (Manual).ConclusionThe present method is a versatile and useful method of improving the objectivity of film dosimetry for IMRT QA.     

Chirp phantom for MTF calculations. A study of its precision in noisy environments

Bar pattern phantoms are used to determine the maximum number of line-pairs per mm that an imaging system can resolve. In some cases, a numerical determination of the modulation transfer function (MTF) can also be carried out. However, calculations can only be performed in a relatively small number of frequencies because of the small number of bar groups in the phantom. In this work, a new bar pattern phantom has been simulated. This phantom consists of 66 pairs of lines of different periods and these periods vary exponentially with spatial position, like in a chirp wave. An oversampling procedure has been implemented to obtain the pre-sampled MTF of the system and the results obtained have been compared with those obtained with the edge method, recommended by the IEC. Monte Carlo simulations were carried out for three different levels of noise aimed at investigating the effect of noise on the uncertainties of the MTF determination. In addition, using the analytic expressions for the MTF calculation, statistical fluctuations of noise in phantom images were propagated to MTF values. Despite the smaller size of the chirp phantom, uncertainties in the chirp method are smaller than those of the edge method. For the edge image, the standard deviation of the MTF is proportional to the frequency f, whereas for the chirp method it is proportional to its square root. It is shown that applying an oversampling method allows the use of a single line pair per period without compromising the precision in noisy environments.     

Experimental assessment of the influence of beam hardening filters on image quality and patient dose in volumetric 64-slice X-ray CT scanners

Beam hardening filters have long been employed in X-ray Computed Tomography (CT) to preferentially absorb soft and low-energy X-rays having no or little contribution to image formation, thus allowing the reduction of patient dose and beam hardening artefacts. In this work, we studied the influence of additional copper (Cu) and aluminium (Al) flat filters on patient dose and image quality and seek an optimum filter thickness for the GE LightSpeed VCT 64-slice CT scanner using experimental phantom measurements. Different thicknesses of Cu and Al filters (0.5–1.6 mm Cu, 0.5–4 mm Al) were installed on the scanner’s collimator. A planar phantom consisting of 13 slabs of Cu having different thicknesses was designed and scanned to assess the impact of beam filtration on contrast in the intensity domain (CT detector’s output). To assess image contrast and image noise, a cylindrical phantom consisting of a polyethylene cylinder having 16 holes filled with different concentrations of K₂HPO₄ solution mimicking different tissue types was used. The GE performance and the standard head CT dose index (CTDI) phantoms were also used to assess image resolution characterized by the modulation transfer function (MTF) and patient dose defined by the weighted CTDI. A 100 mm pencil ionization chamber was used for CTDI measurement. Finally, an optimum filter thickness was determined from an objective figure of merit (FOM) metric. The results show that the contrast is somewhat compromised with filter thickness in both the planar and cylindrical phantoms. The contrast of the K₂HPO₄ solutions in the cylindrical phantom was degraded by up to 10% for a 0.68 mm Cu filter and 6% for a 4.14 mm Al filter. It was shown that additional filters increase image noise which impaired the detectability of low density K₂HPO₄ solutions. It was found that with a 0.48 mm Cu filter the 50% MTF value is shifted by about 0.77 lp/cm compared to the case where the filter is not used. An added Cu filter with approximately 0.5 mm thickness accounts for 50% reduction in radiation-absorbed dose as measured by the weighted CTDI. The FOM results indicate that with an additional filter of 0.5 mm Cu or minimum 4 mm Al, a good compromise between image quality and patient dose is achieved for CT images acquired at tube voltages of 120 and 140 kVp. The results seem to indicate that an optimum filter for high kVp acquisitions, routinely used in cardiovascular imaging, should be 0.5 mm copper or 4 mm aluminium minimum.     

Monte Carlo simulation and measurement of radiation leakage from applicators used in external electron radiotherapy

External electron radiotherapy is performed using a cone or applicator to collimate the beam. However, because of a trade-off between collimation and scattering/bremsstrahlung X-ray production, applicators generate a small amount of secondary radiation (leakage). We investigate the peripheral dose outside the radiation field of a Varian-type applicator. The dose and fluence outside the radiation field were analyzed in a detailed Monte Carlo simulation. The differences between the calculation results and data measured in a water phantom in an ionization chamber were less than ±1% in regions more than 3 mm below the surface of the phantom and at the depth of dose maximum. The calculated fluence was analyzed inside and outside the radiation field on a plane just above the water phantom surface. Changing the electron energy affected the off-axis fluence distribution outside the radiation field; however, the size of the applicator had little effect on this distribution. For each energy, the distributions outside the radiation field were similar to the dose distribution at shallow depths in the water phantom. The effect of secondary electrons generation by photon transmission through the alloy making up the lowest scraper was largest in the region from the field edge to directly below the cutout and at higher beam energies. The results of the Monte Carlo simulation confirm that the peripheral dose outside the field is significantly affected by radiation scattered or transmitted from the applicator, and the effect increases with the electron energy.     

Visualization of air and metal inhomogeneities in phantoms irradiated by carbon ion beams using prompt secondary ions

PurposeNon-invasive methods for monitoring of the therapeutic ion beam extension in the patient are desired in order to handle deteriorations of the dose distribution related to changes of the patient geometry. In carbon ion radiotherapy, secondary light ions represent one of potential sources of information about the dose distribution in the irradiated target. The capability to detect range-changing inhomogeneities inside of an otherwise homogeneous phantom, based on single track measurements, is addressed in this paper.MethodsAir and stainless steel inhomogeneities, with PMMA equivalent thickness of 10 mm and 4.8 mm respectively, were inserted into a PMMA-phantom at different positions in depth. Irradiations of the phantom with therapeutic carbon ion pencil beams were performed at the Heidelberg Ion Beam Therapy Center. Tracks of single secondary ions escaping the phantom under irradiation were detected with a pixelized semiconductor detector Timepix. The statistical relevance of the found differences between the track distributions with and without inhomogeneities was evaluated.ResultsMeasured shifts of the distal edge and changes in the fragmentation probability make the presence of inhomogeneities inserted into the traversed medium detectable for both, 10 mm air cavities and 1 mm thick stainless steel. Moreover, the method was shown to be sensitive also on their position in the observed body, even when localized behind the Bragg-peak.ConclusionsThe presented results demonstrate experimentally, that the method using distributions of single secondary ion tracks is sensitive to the changes of homogeneity of the traversed material for the studied geometries of the target.