ISSN: 2320-2459

**Habibeh Ghasemi ^{1}, Zohreh Azma^{2}, Ali Jabbary Arfaee^{3} and Mahdi Sadeghi^{4}^{*}**

^{1}Medical Radiation Engineering Department, Science and Research Branch, Tehran, Iran

^{2}Medical Radiation Engineering Department, Shahid Beheshti University, Tehran, Iran

^{3}Shohada Tajrish Hospital Radiation Oncology, Shahid Beheshti University, Tehran, Iran

^{4}Medical Physics Department, School of Medicine, Iran University of Medical Science, Tehran, Iran

- *Corresponding Author:
- Mahdi Sadeghi

Medical physics department,

School of Medicine, Iran University of Medical Science,

Tehran, Iran

**Tel:**+989123595322

**E-mail:**mahdisadeghi2003@yahoo.com

**Received Date**: 10/07/2018; **Accepted Date:** 26/09/2018; **Published Date**: 04/10/2018

**Visit for more related articles at** Research & Reviews: Journal of Pure and Applied Physics

Background: In the treatment with the electron beams to get the appropriate dose distribution in the tumor region, the distance between virtual electron source position and patient’s body should be correctly identified. This study is an attempt to find the effective SSD (Source-surface distance)by using the Monte Carlo N-particle (MCNP v1. 51) code to substitute the measurements in the clinic. Materials and Methods: The MC (Monte Carlo) simulation of Oncor Linac (Siemens Co. Germany) was done based on manufacturer data in 9 MeV electron beam energy. In order to obtain the values of effective SSD using the inverse slope (IS) technique, point doses are calculated at dmax (the depth of maximum dose on the central axis), inside of the water phantom model, with the phantom simulated first at the standard SSD (5 cm air gap) and then with an increasing air gap of 3, 6, 9 and 12 cm between the distal end of applicator and the phantom surface. Results: The measured and MC-calculated effective SSDs using IS method are presented and compared with those reported from other works. The effective SSDs were found to be 95.68 cm and 96.66 cm from the MC tallies of *F8 and *F4, respectively. In addition, the effective SSDs were found to be 95.59 cm from ion chamber measurement. Conclusion: Considering the results from three methods of calculation, and comparing the results with experimental methods, it was shown that Monte Carlo simulation is a useful method in determination of electron virtual source position.

Electron beams have an important role in the treatment of superficial and shallow cancers in modern radiotherapy [1].
Because of the electron interactions with the accelerator head materials including scattering foils, monitoring chambers, photon
jaws, and the applicators, they may be thought of as emanating from the virtual source that is not in the real position of the
accelerator exit window. The International Commission on Radiation Units and Measurements (ICRU) Report No. 35 [2], defines
an effective extended electron source as the source which when placed in a vacuum at some distance SSD from the phantom
surface (Z=0) would produce exactly the same electron fluence at Z=0 as the actual beam. The distance from this source position
to the patient’s skin surface is called the effective SSD. To obtain an accurate dose calculation of electron beams, correction
for measured effective SSD must be applied. In addition, effective SSD may be used to calculate output at extended SSD for
corrections of central axis percentage depth dose (PDD) [3], off-axis dose (OAD) values, output factors, and X-ray contamination [4].
Hence, it is essential that these corrections should be examined for each treatment unit prior to referral to treat at extended SSD
using electron beams. The SSD_{eff} is dependent on the electron energy and the field size used.

There are three independent experimental techniques to ascertain the position of the effective electron source. These are projection of the wires grid image to a photographic film, variation of field size with the nominal source-to-film distance, and variation of dose rate with the distance [5-7]. The last method, namely the ‘inverse slope’ proposed by Khan [8,9] was chosen in the present work since it simulates the situation in which we need to know the dose rate from the electron beam output of a Siemens Oncor Linac machine (Siemens Co. Germany) at different extended SSD. Khan et al. [10] have recommended an inverse slope (IS) method that simulates as closely as possible the clinical circumstance. In this regard, by taking:

(1)

Which g is the air gap between the distal end of the applicator and the surface of the water phantom, *f* is the effective SSD
value, *I _{0}* is dose at d

f=(1/*slope*) – *d _{m}* (2)

Herein, we attempted to estimate the effective electron source position using measurements and Monte Carlo simulation from a Siemens Oncor LINAC at the Alborz Hospital (Karaj, Iran).

**Measurement**

In the current work, to verification of the simulated beam's dosimetric properties, the PDD along clinical axis data and the
OAD values both were measured within a 50 × 50 × 50 cm^{3} PTW scanning water phantom tank (MP3-m, PTW) according to
recommendations of the Technical Reports series No. 398 protocol, reported by International Atomic Energy Agency (IAEA) [11].
Measurements were performed along the central axis by using a 0.02 cc nominal sensitive volume Advanced Markus^{®} electron
chamber (TM 34045; PTW, Freiburg, Germany) and acquired data were processed using the MEPHYSTO (Medical Physics Tools)
mc2 software v7.20 (PTW-beam analyzer) at the Alborz Hospital (Karaj, Iran). The Advanced Markus chamber is a vented parallel
plate with an approximate moving speed of 50 mm/sec and had a 0.87 mm thick PMMA (Poly Methyl Meth Acrylate) protective cap
for in-water measurements. Chamber positions inside the phantom were automatically controlled using a TBA control unit (PTW).
The PDD step sizes were 2 mm, 5 mm and 10 mm for depths shallower than 3 cm, 4 cm and 8 cm, respectively. The profile step
sizes were set to 2 mm in the penumbra and 5mm everywhere else.

In order to obtain the values of effective SSD, point doses are measured by radiating 100 monitor unit (MU) at near the
depth of maximum dose, d_{max} (i.e. 2.4 cm) inside of the water phantom, with the phantom first at the standard SSD (no air gap)
and then with an increasing air gap of 3, 6, 9 and 12 cm between the distal end of applicator and the phantom surface. Although
the effective SSD is acquired through doing measurements at the d_{max}, its value does not change considerably with the depth of
measurement [10]. This procedure was carried out for electron beam energy of 9 MeV using applicator sizes of 10 × 10 cm^{2}. The
uncertainty in the measurement of dose values is obtained less than 2%.

**Monte Carlo Simulation**

In this work, the MCNP5 v1.51 Monte Carlo code [12-18] was employed to simulate the radiation transport of electron beams through the Siemens ONCOR Linac (Siemens Co. Germany) treatment head and water phantom based on the dimensions and materials provided by the manufacturer. The accelerator exit window, primary collimators, flattening filter, monitor chambers, Y-jaws, and X-jaws were simulated in the geometric section of the code. The electron source was simulated with a 0.1 mm radius in disk form, which produced 9 MeV electron beams and had a Gaussian energy distribution with a full width at half maximum (FWHM) of 1.007 MeV and centered at 10.06 MeV for the 9 MeV electron beam. A view of simulated geometry for the treatment head of Siemens Oncor LINAC is shown in Figure 1.

**Dose Calculations**

The percentage depth dose along central axis and off-axis dose values acquired in a water phantom with dimensions of 50
× 50 × 50 cm^{3} with an SSD of 100 cm. Since no mesh tally is available for pulse height tally *F8 with full physics transport (p,
e-mode) in the MCNP5 code, a three dimension (3D) tally distribution of cells were applied in the simulation. To achieve depth
dose calculations within various depths of the water phantom, a narrow rectangle region was defined with 5 mm lateral width and
160 mm height along the beam central axis. The narrow rectangle was then split into tiny scoring voxels with size varied between
5.0 × 5.0 × 1.0 and 5.0 × 5.0 × 7.0, depending on the spatial resolution required. Usually, a smaller voxel size was chosen for the
depth dose buildup region and the beam penumbra region [19]. Applying the same procedure was followed to calculate transverse
dose values, except that the main axes of scoring narrow rectangle was considered perpendicular to the central axis of the beam.
To plot PDD and OAD curves, all absorbed dose values were normalized to the maximum dose in *d _{max}*. In the simulations for PDD
and OAD calculations, the cut-off energy value was defined as 0.01 and 0.05 MeV for the photons and electrons, respectively. A
total number of 1 × 10

In order to obtain the values of effective SSD, point doses are calculated at *d _{max}*, inside of the 50 × 50 × 50 cm

**Depth Dose Distribution and Dose Profile**

To validate our simulation results, the MC-calculated PDD and OAD values in a water phantom were firstly obtained for a
9-MeV electron beam using open 10 × 10 cm^{2} applicator from Siemens ONCOR Linac and then compared with the experimental
data measured at the Alborz Hospital, Karaj, Iran. Figures 2 and 3 illustrate the compare of experimental versus MC-based
theoretical relative depth dose curves from 9-MeV electron beams for the field size of 10 × 10 cm^{2} in water. As it can be seen
from Figure 2 the MC-calculated PDD overestimates the dose in the buildup region with local differences up to 4.5%. The larger
deviation in the build-up region is due to a high gradient of dose distribution, which makes ionization chamber measurements
unreliable as well as the waterproof enclosure of the chamber, which affects the PDD at the surface of the phantom [21]. Also, the
ﬁnite size of the chamber, which perturbs the absorbed dose, may be another reason [22]. Beyond *d _{max}*, differences are less than
1% from R

Off-axis dose values were determined at the depths of *d _{max}* within the water phantom. The measured and calculated profiles
for 9-MeV electron beams are presented in Figure 4. The relative statistical errors on the MC results were less than 2% for dose
fall-off region and reach a maximum of 6% outside the radiation field. As Figure 5 depicts, the agreement between measurement
and MC calculations was within 1.7% for plateau region, but it increased to 4.6% for the region located outside the radiation ﬁeld,
which can be ensure deviation limits [23,24]. These differences occurred because the effects of electronic components and the lead
blocks shielding the target, the monitor chamber, and the jaw inside the treatment head were negligible and not included in the
design drawing. The design drawing of the electron applicators does not describe an accurate width or thickness of each scraper
in detail, although the thickness of each scraper is given. The error in the thickness measured using a digital caliper affected the
results calculated outside the applicator [25]. The precision of the caliper measurements was about 1 mm. The Figure 5 exhibits
a difference of less than 0.6 mm for 80-20% penumbra width between measured and MC calculated results.

**Table 1.** Differences for measured and calculated percentage depth dose.

Energy (MeV) | Differences up to d (%)_{max} |
Differences at d (mm)_{max} |
Differences at R_{50} (mm) |
Differences at R_{p} (mm) |
---|---|---|---|---|

9 | 1.95 | 0.5 | 1 | 1.5 |

**Table 2.** Differences for measured and calculated off-axis dose values.

Energy (MeV) | Differences at plateau area (%) | Differences at D_{50} (mm) |
Differences a bremsstrahlung tail (%) |
---|---|---|---|

9 | 1.7 | 0.7 | 3.07 |

**The Effective SSD**

To ﬁnd the effective source position, measurement and calculations were carried out at various SSD values ranging from
100 to 112 cm for a radiation field of 10 × 10 cm^{2} at 9-MeV electron beam energy. The inverse slope (IS) technique was used to
ﬁnd the effective SSD as recommended by Khan et al. [26]. For clinical requirements, the IS technique is recommended, because
it is more clinically related as well as give accurate inverse-square corrections [26]. For instance, the usage of IS technique is
shown in Figure 6. It can be seen that the graphs were drawn by indicating the extended SSDs on the z-axis, and square-root
value of *I _{0}/I_{g}* at the

**Table 3. **Comparison of calculated and measured extended source-to-surface distance (SSD) factors.

SSD (cm) | Measured† | Calculated (*F8) | Calculated (*F4) | Uncertainty, *F8 (%) | Uncertainty, *F4 (%) |
---|---|---|---|---|---|

100 | 651 | 1.418e-14 | 8.605e-5 | 1.8 | 1.5 |

103 | 615.33 | 1.325e-14 | 8.31e-5 | 2 | 1.8 |

106 | 580.33 | 1.277e-14 | 7.766e-5 | 2.1 | 1.8 |

109 | 545.67 | 1.197e-14 | 7.355e-5 | 2.1 | 1.9 |

112 | 517.81 | 1.12e-14 | 6.86e-5 | 2.2 | 2 |

**Table 3** compares measured and calculated extended SSD factors for the radiation field of 10 × 10 cm^{2} at 9-MeV electron
beam energy included in the study. The maximum differences of 0.05% and 1.06% were observed between measurement and MC-
*F8 and MC-*F4 calculation, respectively. Absorbed doses resulting from the use of an *F4 volume flux tally method, compared
to an *F8 tally using both the photon-electron mode (p, e-mode) were all in agreement.

The results of the measurements, obtained with the inverse slope method, have been compared with those reported from
the literature for a 10 × 10 cm^{2} applicator and the selected 9 MeV energy in **Table 4**. in early studies of effective SSD, the amount
of effective SSD in the mentioned position calculated 83.2 Cm by khan et al. [31] in Toshiba (LMR 13) machine. In the latest studies
Rajasekar et al. [32] in Mitsubishi machine and Kim et al. [31] in Varian, calculated SSD_{eff}, respectively 87.2 and 80.6 Cm. As **Table 4** depicts, the comparison of virtual SSD obtained in this study with those obtained by other authors as well as suggest that SSD_{eff} is nm dependent and should be measured for each accelerator [33]. It should be pointed out that these differences may be caused
by differences within the collimation systems of LINACs as outlined previously. Because collimation system has various beam
scattering components and as a result, the SSD_{eff} varies with collimation system in a complex manner [34].

**Table 4. **Comparison of the measured and calculated effective SSDs of the present study with those of other workers for the 10 × 10 cm^{2} applicator in 9-MeV energy.

Authors | Method, devices | Accelerator type | SSD_{eff} |
---|---|---|---|

Jamshidi et al. | ISL†, 0.6-cm3 Farmer-type ionization chamber with a ^{60}Co cap |
Varian Clinac-2500 | 72.7 |

Faermann et al. | IS††, 0.6 cc cylindrical ion chamber, model 2571 NE | Varian Clinac 18 | 76.18 |

Kim et al. | IS, Farmer type chamber (FC 65-G) | iX (Varian) | 80.6 |

Roback et al. | IS, 0.6-cm^{3} Farmer ionization chamber |
Varian Clinac 2500 | 82.9 |

Khan et al. | IS, 0.6-ml Farmer ion chamber | Toshiba (LMR 13) | 83.2 |

Rajasekar et al. | IS, radiation ﬁeld analyser (RFA-300, Scanditronics) with a p-semiconductor probe | Mitsubishi | 87.2 |

Roback et al. | IS, 0.6-cm^{3} Farmer ionization chamber |
Varian Clinac 2100C | 88.3 |

Ravindran et al. | IS, 0.5 cm^{3} Capintec parallel plate chamber |
Siemens Mevatron | 92.2 |

Cygler et al. | ISL, RK cylindrical chamber (0.12 cc) | Siemens MD-2 | 94 |

This work | IS, 0.02cc Advanced Markus® electron chamber (TM 34045; PTW) | Siemens ONCOR | 95.59 |

This work | IS, MCNP5-*F8 | Siemens ONCOR | 95.68 |

This work | IS, MCNP5-*F4 | Siemens ONCOR | 96.66 |

In the present study, we have simulated 9 MeV electron beams generated from the Siemens Oncore Linac treatment head by MCNP5 transport code. Considering the good agreements between the results of MC-calculated PDD and OAD values with the corresponding measured values and the model simulation was veriﬁed. The simulation results show that the differences between measured and calculated values of effective SSD within 1%. Considering the results from three methods of calculation, and comparing the results with experimental methods, it was shown that Monte Carlo simulation is a useful method in determination of electron virtual source position and can be used safely in upcoming dose calculations in uncommon SSD distances encountered in the clinic.

Author’s would like to thank the “Alborz Hospital, Karaj, Iran” for their support.

- OShea T, Foley MJ, Rajasekar D, Downes PA, Van der Putten W, et al. Electron beam therapy at extended source-to-surface distance: a Monte Carlo investigation. J Appl Clin Med Phys 2008;9:2811.
- Ibbott GS. Radiation Dosimetry: Electron Beams with Energies between 1 and 50 MeV (ICRU Report No. 35). Med Phys 1985;12:813.
- Gibbons JP, Antolak JA, Followill DS, Huq MS, Klein EE, et al. Monitor unit calculations for external photon and electron beams: Report of the AAPM Therapy Physics Committee Task Group No. 71. Med Phys 2014;41.
- Laughlin JS. High energy electron treatment planning for inhomogeneities. Br J Radiol 1965;38:143-147.
- Wolfgang P, Julius K. Dosimetrie zur Betatrontherapie. Thieme, West Germany 1965.
- Fehrentz D, Liebig B, Schröder-Babo P. Consideration of large in homogeneous regions of the computation of electron dose distributions. Strahlentherapie 1976;151:423-442.
- Briot E, Dutreix A. Dosimetry of high energy electron beams from a linear accelerator (author's transl). J Radiol Electrol Med Nucl 1976;57:447-454.
- Khan FM. The Physics of Radiation Therapy. Baltimore: Williams & Wilkins 1984.
- Khan FM. Clinical Electron Beam Dosimetry. In: Keriakes JG, elson hr, Born CG, eds. Radiation Oncology Physics. AAPM Monograph No. 15. New York, NY: American Institute of physics 1986:211-264.
- Khan FM, Sewchand W, Levitt SH. Effect of air space on depth dose in electron beam therapy. Radiology 1978;126:249
- Andreo P, Burns DT, Hohlfeld K, Huq MS, Kanai T, et al. IAEA TRS-398 Absorbed dose determination in external beam radiotherapy: An International code of practice for dosimetry based on standards of absorbed dose to water 2000.
- Monte Carlo Team. Monte Carlo N-particle extended Los Alamos National Laboratory Report 2008.
- Saidi P, Sadeghi M, Shirazi A, Tenreiro C. ROPES eye plaque brachytherapy dosimetry for two models of 103Pd seeds. Australas Phys Eng Sci Med 2011;34:223-231.
- Hadadi A, Sadeghi M, Sardari D, Khanchi A, Shirazi A, et al. Monte Carlo characterization of biocompatible beta-emitting 90Y glass seed incorporated with the radionuclide 153Sm as a SPECT marker for brachytherapy applications. Journal of Applied Clinical Medical Physics 2013;14:90-103.
- Hosseini H, Sadeghi M, Ataeinia V. Dosimetric comparison of four new design 103Pd brachytherapy sources: optimal design using silver and copper rod cores. Med Phys 2009;36:3080-3085.
- Saidi P, Sadeghi M, Shirazi A, Tenreiro C. Dosimetric parameters of the new design 103Pd brachytherapy source based on Monte Carlo study. Physica Medica 2011;28:13-18.
- Saidi P, Sadeghi M, Shirazi A, Tenreiro C. Monte Carlo calculation of dosimetry parameters for the IR08-103Pd brachytherapy source. Medical Physics 2010;37:2509.
- Sadeghi M, Hosseini SH, Raisali G. Experimental measurements and Monte Carlo calculations of dosimetric parameters of the IRA1-103Pd brachytherapy source. Applied Radiation and Isotopes 2008;66:1431-1437
- Cho SH, Vassiliev ON, Lee S, Liu H H, Ibbott GS, et al. Reference photon dosimetry data and reference phase space data for the 6 MV photon beam from Varian Clinac 2100 series linear accelerators. Med Phys 2005;32:137-148
- http://physics.nist.gov/PhysRefData/XrayMassCoef/tab4.html
- Mohammadi N, Miri-Hakimaba H, Rafat-Motavali L, Akbari F, Abdollahi S. Neutron spectrometry and determination of neutron contamination around the 15 MV Siemens Primus LINAC. J Radioanal Nucl Chem 2015;304:1001-1008.
- Siebers JV, Keall PJ, Libby B, Mohan R. Comparison of EGS4 and MCNP4b Monte Carlo codes for generation of photon phase space distributions for a Varian 2100C. Phys Med Biol 1999;44:3009-3026.
- Brahme A, Chavaudra J, Landberg T, McCullough EC, Nusslin F, et al. Accuracy requirements and quality assurance of external beam therapy with photons and electrons 1988.
- Wieslander E, Knoos T. A virtual linear accelerator for veriﬁcation of treatment planning systems. Phys Med Biol 2000;45:2887-2896.
- Shimozato T, Okudaira K, Fuse H, Tabushi K. Monte Carlo simulation and measurement of radiation leakage from applicators used in external electron radiotherapy. Med Phys 2013:388-396.
- Khan FM, Sewchand W, Levitt SH. Effect of airspace on depth dose in electron beam therapy. Radiology 1978;126:249-251.
- Faermann S, Lavie K, Datz H. Determination of the Virtual Source Position for the Electron Beams of a Varian Clinac 18 Linear Accelerator: a comparison of experimental methods. Annual Meeting of the Israel Nuclear Societies 2002.
- Jamshidi A, Kuchnir FT, Reft CS. Determination of the source position for the electron beams from a high energy linear accelerator. Med Phys 1986;13:942-945.
- Roback DM, Khan FM, Gibbons JP, Sethi A. Effective SSD for electron beams as a function of energy and beam collimation. Med Phys 1995;22:2093-2095.
- Faermann S, Krutman M, Weinstein M. A comparison of the methods of determination of virtual sources for electron beams Proc 1983.
- Kim MT, Lee HK, Heo YC, Cho JH. A Study on Effective Source-Skin Distance using Phantom in Electron Beam Therapy. Journal of Magnetics 2014;19:15-19.
- Rajasekar D, Datta NR, Das KJ, Ayyagari S. Electron beam therapy at extended SSDs: an analysis of output correction factors for a Mitsubishi linear accelerator. Phys Med Biol 2002;47:3301-3311.
- Ravindran BP. A study on virtual source position for electron beams from a Mevatron MD linear accelerator. Phys Med Biol 1999;44:1309-1315.
- Cygler J, Li XA, Ding GX, Lawrence E. Practical approach to electron beam dosimetry at extended SSD. Phys Med Biol 1997;42:1505-1514.