Grigor N. Grigorov

Monte Carlo simulations of dose near a nonradioactive gold seed

The relative doses and hot/cold spot positions around a non-radioactive gold seed, irradiated by a 6 or 18MV photon beam in water, were calculated using Monte Carlo simulation. Phase space files of 6 and 18MV photon beams with a field size of 1×1cm2 were generated by a Varian 21 EX linear accelerator using the EGSnrc and BEAMnrc code. The seed (1.2×1.2×3.2mm3) was positioned at the isocenter in a water phantom (20×20×20cm2) with source-to-axis distance=100cm. For the single beam geometry, the relative doses (normalized to the dose at 5mm distance above the isocenter) at the upstream seed surface were calculated to be 1.64 and 1.56 for the 6 and 18MV beams respectively when the central beam axis (CAX) is parallel to the width of the seed. These doses were slightly higher than those (1.58 and 1.52 for 6 and 18MV beams respectively) calculated when the CAX is perpendicular to the width of the seed. Compared to the relative dose profiles with the same beam geometry without the seed in the water phantom, the presence of the seed affects the dose distribution at about 3mm distance beyond both the upstream and downstream seed surface. For a pair of opposing beams with equal and unequal beam weight, the hot and cold spots of both opposing beams were mixed. For a 360degree photon arc around the longitudinal axis of the seed, the relative dose profile along the width of the seed was similar to that of the opposing beam pair, except the former geometry has a larger dose gradient near the seed surface. In this study, selected results from our simulation were compared to previous measurements using film dosimetry.

Dose measurements near a non-radioactive gold seed using radiographic film

The dose distribution near a non-radioactive gold seed under a 6 MV photon beam was measured using radiographic film, water equivalent bolus and solid water slabs. This type of small seed is typically used as a marker in target positional verification using a portal imager for conformal prostate treatment such as intensity modulated radiation therapy. A stack of three films was placed on top of the seed located on a soft bolus. Solid water slabs were then placed on top of the film. The films were exposed using a small 1x1 cm2 field. Then, using a similar experimental set-up and exposure, another stack of three films was placed under the seed, which was then covered by the soft bolus and solid water slabs. The cross-plane axial beam profiles at different depths, depending on the thickness of the film package, were measured. From the group of beam profiles above and below the seed, the dose distribution along a selected vertical line within the profiles was easily plotted. Compared to the dose with no seed at the isocentre and 5 cm of solid water, there was about a 21% increase in dose at 0.35 mm above the seed. On the other hand, there was about a 22% decrease in dose at the same distance below the seed. The dosimetry of the calibrated film was verified with a MOSFET detector. The change in dose due to the seed by varying the incident beam angles was also measured for this note.

Peripheral dose outside applicators in electron beams

The peripheral dose outside the applicators in electron beams was studied using a Varian 21 EX linear accelerator. To measure the peripheral dose profiles and point doses for the applicator, a solid water phantom was used with calibrated Kodak TL films. Peak dose spot was observed in the 4 MeV beam outside the applicator. The peripheral dose peak was very small in the 6 MeV beam and was ignorable at higher energies. Using the 10 x 10 cm(2) cutout and applicator, the dose peak for the 4 MeV beam was about 12 cm away from the field central beam axis (CAX) and the peripheral dose profiles did not change with depths measured at 0.2, 0.5 and 1 cm. The peripheral doses and profiles were further measured by varying the angle of obliquity, cutout and applicator size for the 4 MeV beam. The local peak dose was increased with about 3% per degree angle of obliquity, and was about 1% of the prescribed dose (angle of obliquity equals zero) at 1 cm depth in the phantom using the 10 x 10 cm(2) cutout and applicator. The peak dose position was also shifted 7 mm towards the CAX when the angle of obliquity was increased from 0 to 15 degrees.

Monte Carlo simulation of backscatter from lead for clinical electron beams using EGSnrc

In electron radiotherapy of superficial lesions in the eyelid, lip, buccal mucosa, ear, and nose, backscattered electrons are produced from the lead shield used to protect the critical tissue underneath the tumor. In this study, the backscattered electrons, produced by clinical electron beams using a Varian 21 EX linear accelerator, were studied using Monte Carlo simulations. The electron backscatter factor (EBF), defined as the ratio of dose at the tissue-lead interface to the dose at the same point without the presence of backscatter, was calculated using the Monte Carlo EGSnrc-based code. The calculated EBFs were verified with measurements using metal-oxide-semiconductor field effect transistor detectors. The effect of the (1) initial electron beam energy, (2) thickness of bolus over the lead shield, (3) beams angle of incidence, and (4) presence of an aluminum sheet used to absorb backscattered electrons, on the EBF, were studied. It is found that for lead shielding positioned at any fixed depth, the EBF decreases with an increase in initial electron beam energy (4-16 MeV). In addition, for depths within the electron practical range, Rp, and at a particular beam energy, the EBF increases with depth (or thickness of the treatment volume). When the electron beam angle increases from 0 degrees to 5 degrees, the EBF only decreases slightly (<4%) for all energies. The influence of the beam obliquity on the EBF is important when the treatment surface is not flat and perpendicular to the central beam axis. The use of an aluminum sheet to reduce backscattered electrons was also investigated. For a relatively low electron beam energy (4 MeV), a 2 mm aluminum sheet can reduce backscattering by 31%. While the electron beam energy increased, less backscattered electrons were produced and therefore removed by the same thickness of aluminum (only about 6% for 16 MeV). The Monte Carlo calculated EBFs from this study, characterized by the electron beam energy, depth of bolus above the lead shield, beam obliquity, and presence of an aluminum sheet, may provide important clinical information for radiation oncology staff when considering the effect of electron backscatter on radiotherapy using internal shielding.

Surface dosimetry for oblique tangential photon beams: a Monte Carlo simulation study.

The effect of beam obliquity on the surface relative dose profiles for the tangential photon beams was studied. The 6 and 15 MV photon beams with 4 x 4 and 10 x 10 cm2 field sizes produced by a Varian 21 EX linear accelerator were used. Phase-space models of the photon beams were created using Monte Carlo simulations based on the EGSnrc code, and were verified using film measurements. The relative dose profiles in the phantom skin, at 2 mm depth from the surface of the half-phantom geometry, or HPG, were calculated for increasing gantry angles from 270 to 280 deg clockwise. Relative dose profiles of a full phantom enclosing the whole tangential beam (full phantom geometry, or FPG) were also calculated using Monte Carlo simulation as a control for comparison. The results showed that, although the relative dose profiles in the phantom skin did not change significantly with an oblique beam using a FPG, the surface relative depth dose was increased for the HPG. In the HPG, with 6 MV photon beams and field size = 10 x 10 cm2, when the beam angle, starting from 270 deg, was increased from 1 to 3 deg, the relative depth doses in the phantom skin were increased from 68% to 79% at 10 cm depth. This increase in dose was slightly larger than the dose from 15 MV photon beams with the same field size and beam angles, where the relative depth doses in phantom skin were increased from 81% to 87% at 10 cm depth. A parameter called the percent depth dose (PDD) ratio, defined as the relative depth dose from the HPG to the relative depth dose from the FPG at a given depth along the phantom skin, was used to evaluate the effect of the phantom-air interface. It is found that the PDD ratio increased significantly when the beam angle was changed from zero to 1-3 degrees. Moreover, the PDD ratio, for a given field size, experienced a greater increase for 6 MV than for 15 MV. For the same photon beam energy, the PDD ratio increased more with a 4 x 4 cm2 field compared to 10 x 10 cm2. The results in this study will be useful for physicists and dosimetrists to predict the surface relative dose variations when using clinical tangential-like photon beams in radiation therapy.

SWIMRT: A graphical user interface using the sliding window algorithm to construct a fluence map machine file

A custom‐made computer program, SWIMRT, to construct “multileaf collimator (MLC) machine” file for intensity‐modulated radiotherapy (IMRT) fluence maps was developed using MATLAB® and the sliding window algorithm. The user can either import a fluence map with a graphical file format created by an external treatment‐planning system such as Pinnacle3 or create his or her own fluence map using the matrix editor in the program. Through comprehensive calibrations of the dose and the dimension of the imported fluence field, the user can use associated image‐processing tools such as field resizing and edge trimming to modify the imported map. When the processed fluence map is suitable, a “MLC machine” file is generated for our Varian 21 EX linear accelerator with a 120‐leaf Millennium MLC. This machine file is transferred to the MLC console of the LINAC to control the continuous motions of the leaves during beam irradiation. An IMRT field is then irradiated with the 2D intensity profiles, and the irradiated profiles are compared to the imported or modified fluence map. This program was verified and tested using film dosimetry to address the following uncertainties: (1) the mechanical limitation due to the leaf width and maximum traveling speed, and (2) the dosimetric limitation due to the leaf leakage/transmission and penumbra effect. Because the fluence map can be edited, resized, and processed according to the requirement of a study, SWIMRT is essential in studying and investigating the IMRT technique using the sliding window algorithm. Using this program, future work on the algorithm may include redistributing the time space between segmental fields to enhance the fluence resolution, and readjusting the timing of each leaf during delivery to avoid small fields. Possible clinical utilities and examples for SWIMRT are given in this paper. PACS numbers: 87.53.Kn, 87.53.St, 87.53.Uv

Electron radiotherapy: a study on dosimetric uncertainty using small cutouts

This note investigated the dosimetric uncertainties due to the positional error when centring a small cutout to the machine central beam axis (CAX) in electron radiotherapy. A group of six circular cutouts with 4 cm diameter were made with their centres shifting 0, 2, 4, 6, 8 and 10 mm from the machine CAX for the 6 x 6 cm(2) applicator. The per cent depth doses, beam profiles and output factors were measured using the 4, 9 and 16 MeV clinical electron beams produced by a Varian 21 EX linear accelerator. The 2D isodose distributions in the z-x (or cross-line) and z-y (or in-line) plane were calculated by Monte Carlo simulation using the EGSnrc system. When the cutout centre was shifted away from the machine CAX for the 4 MeV beam, the d(m), R(80) and R(90) at the machine CAX had no significant change (<0.1 mm). For higher energies of 9 and 16 MeV beams, the d(m) was reduced by 0.45 and 1.63 mm per mm, between the cutout centre and the machine CAX with off-axis shift <6 mm respectively. R(80) and R(90) were reduced by more than 0.3 mm per mm off-axis shift for both energies. The isodose coverage of the in-line axis beam profile was reduced when the cutout centre was shifted away from machine CAX. It is important for oncology staff to note such dosimetric changes in the clinical electron radiotherapy, particularly when a high energy electron beam is used for small cutout. Such positional uncertainty is unavoidable in fabricating an electron cutout in the mould room.

Effect of electron beam obliquity on lateral buildup ratio: a Monte Carlo dosimetry evaluation

The impact of the oblique electron beam on the lateral buildup ratio (LBR), used in the electron pencil beam model to predict the per cent depth dose (PDD) and dose per monitor unit (MU) for an irregular electron field, was examined using Monte Carlo simulation. The EGSnrc-based Monte Carlo code was used to model electron beams produced by a Varian 21 EX linear accelerator for different beam energies, angles of obliquity and field sizes. The Monte Carlo phase space model was verified by measurements using electron diode and radiographic film. For PDDs of oblique electron beams, it is found that the depth of maximum dose (d(m)) shifts towards the surface as the beam obliquity increases. Moreover, for increasing the beam angle of obliquity, the depth doses just beyond d(m) decrease with depth. The depth doses then increase eventually in a deeper depth close to the practical range. The LBRs and pencil beam radial spread function, calculated using PDDs with different field sizes, are found varying with electron beam energies, angles of obliquity and cutout diameters. It is found that LBR increases along the normalized depth when the beam angle of obliquity increases. This results in a decrease of the radial spread function with an increase of beam obliquity. When the size of the electron field increases, the variation of LBR with beam angle of obliquity decreases. It should be noted that when calculating dose per MU for an oblique electron beam with an irregular field misunderstanding and neglecting the effect of beam obliquity would lead to a significant deviation. A database of LBRs for oblique electron beams can be created using Monte Carlo simulation conveniently and is recommended when an oblique beam is used in electron radiotherapy.

Dosimetric dependence of the dimensional characteristics on a lead shield in electron radiotherapy: a Monte Carlo study

This study investigates the dosimetric dependence of the dimension of a lead (Pb) layer for shielding using clinical electron beams with different energies. Monte Carlo simulations were used to generate phase space files for the 4, 9 and 16 MeV electron beams produced by a Varian 21 EX linear accelerator using the EGSnrc‐based BEAMnrc code, and validated by measurements using films. Pb layers with different thicknesses (2, 4, 6 and 8 mm) and diameters (2.5, 3, 3.5 and 4 cm) were placed at the center of an electron field on a solid water phantom. Beam profiles were determined at the depth of maximum dose (dm) using Monte Carlo simulations. The dose profiles under the Pb layer at dm, including the penumbra at the edge of the layer and relative dose at the central beam axis (CAX), were studied with varying thicknesses and diameters of Pb. It is found that 2 mm of Pb is adequate to provide 5 half value layer (HVL) attenuation for the 4 MeV electron beams, and the beam profiles at dm are dependent on the diameter but not the thickness of the Pb. However, for the 9 and 16 MeV electron beams, the relative dose at the CAX and dm depends on both the thickness and diameter of the Pb layer. For 8 mm thickness of Pb, 4 and 5 HVL attenuation of electron beams with energies of 9 and 16 MeV can be achieved at dm, respectively. Moreover, the beam profile under the Pb layer at dm depends on: (1) the penumbra region at the edge of the Pb layer; (2) the beam attenuation varying with the thickness of the Pb layer; (3) the electron side scatter contributing to the CAX under the Pb layer; and (4) the photon contamination produced by the Pb layer. A parameter called “shielding area factor” (defined as the ratio of the length between two points of 50% relative doses in the beam profile at dm to the diameter of the Pb layer) is suggested to predict the required size and thickness of Pb for shielding a target with known dimension at dm. The dosimetric data calculated by Monte Carlo simulations in this study are useful to select the suitable thickness and size of Pb for the protection of critical tissue in electron radiotherapy. PACS number: 87.53.Bn; 87.55.kh and 87.55.km.

Measurement for the MLC leaf velocity profile by considering the leaf leakage using a radiographic film

A method to measure the velocity profile of a multi-leaf collimator (MLC) leaf along its travel range using a radiographic film is reported by considering the intra-leaf leakage. A specific dynamic MLC field with leaves travelling from the field edge to the isocentre line was designed. The field was used to expose a radiographic film, which was then scanned, and the dose profile along the horizontal leaf axis was measured. The velocity at a sampling point on the film can be calculated by considering the horizontal distance between the sampling point and the isocentre line, dose at the sampling point, dose rate of the linear accelerator, the total leaf travel time from the field edge to isocentre line and the pre-measured dose rate of leaf leakage. With the leaf velocities and velocity profiles for all MLC leaves measured routinely, a comprehensive and simple QA for the MLC can be set up to test the consistency of the leaf velocity performance which is essential to the IMRT delivery using a sliding window technique.