• Progress in Medical Physics
• /
• v.19 no.1
• /
• pp.21-34
• /
• 2008

The parallel Monte Carlo electron and photon transport (PMCEPT) code [Kum and Lee, J. Korean Phys. Soc. 47, 716 (2006)] for calculating electron and photon beam doses has been developed based on the three dimensional geometry defined by computed tomography (CT) images and implemented on the Beowulf PC cluster. Understanding the limitations of Monte Carlo codes is useful in order to avoid systematic errors in simulations and to suggest further improvement of the codes. We evaluated the PMCEPT code by comparing its normalized depth doses for electron and photon beams with those of MCNP5, EGS4, DPM, and GEANT4 codes, and with measurements. The PMCEPT results agreed well with others in homogeneous and heterogeneous media within an error of $1{\sim}3%$ of the dose maximum. The computing time benchmark has also been performed for two cases, showing that the PMCEPT code was approximately twenty times faster than the MCNP5 for 20-MeV electron beams irradiated on the water phantom. For the 18-MV photon beams irradiated on the water phantom, the PMCEPT was three times faster than the GEANT4. Thus, the results suggest that the PMCEPT code is indeed appropriate for both fast and accurate simulations.

• Radiation Oncology Journal
• /
• v.16 no.1
• /
• pp.71-79
• /
• 1998

Purpose : The objective of this study is to introduce our installation of a non-commercial 3D Planning system, Plunc and confirm it's clinical applicability in various treatment situations. Materials and Methods : We obtained source codes of Plunc, offered by University of North Carolina and installed them on a Pentium Pro 200MHz (128MB RAM, Millenium VGA) with Linux operating system. To examine accuracy of dose distributions calculated by Plunc, we input beam data of 6MV Photon of our linear accelerator(Siemens MXE 6740) including tissue-maximum ratio, scatter-maximum ratio, attenuation coefficients and shapes of wedge filters. After then, we compared values of dose distributions(Percent depth dose; PDD, dose profiles with and without wedge filters, oblique incident beam, and dose distributions under air-gap) calculated by Plunc with measured values. Results : Plunc operated in almost real time except spending about 10 seconds in full volume dose distribution and dose-volume histogram(DVH) on the PC described above. As compared with measurements for irradiations of 90-cm 550 and 10-cm depth isocenter, the PDD curves calculated by Plunc did not exceed $1\%$ of inaccuracies except buildup region. For dose profiles with and without wedge filter, the calculated ones are accurate within $2\%$ except low-dose region outside irradiations where Plunc showed $5\%$ of dose reduction. For the oblique incident beam, it showed a good agreement except low dose region below $30\%$ of isocenter dose. In the case of dose distribution under air-gap, there was $5\%$ errors of the central-axis dose. Conclusion : By comparing photon dose calculations using the Plunc with measurements, we confirmed that Plunc showed acceptable accuracies about $2-5\%$ in typical treatment situations which was comparable to commercial planning systems using correction-based a1gorithms. Plunc does not have a function for electron beam planning up to the present. However, it is possible to implement electron dose calculation modules or more accurate photon dose calculation into the Plunc system. Plunc is shown to be useful to clear many limitations of 2D planning systems in clinics where a commercial 3D planning system is not available.

• Progress in Medical Physics
• /
• v.10 no.1
• /
• pp.55-62
• /
• 1999

In an effort to study the characteristics of x-rays utilized in radiation therapy, we calculated the energy distribution and the mean energy of x-rays generated from a tungsten target bombarded by 6, 10, and 15 MeV electron beams, using a Monte Carlo technique. The average photon energies calculated as a function of the beam radius lied in 1.4 ∼ 1.6, 2.1 ∼ 2.5 and 2.8 ∼ 3.3 MeV ranges for 4, 10, and 15 MV electron beams, respectively, which turned out to have no strong dependence on the radius. Using the energy distributions of 6,10, and 15 MV x-rays obtained for the target distance of 100 cm, percentage depth doses were determined using Monte Carlo calculations. For the case 10 MV, a comparison was made between our calculation and measurement performed by others. The calculated percentage depth dose appeared somewhat smaller than the measured one except in the surface region. We conclude that this is due to the fact that the beam hardening effect resulting from the flattening filter was not properly allowed for in our Monte Carlo calculations.

• Progress in Medical Physics
• /
• v.13 no.3
• /
• pp.149-155
• /
• 2002

Recent advances in radiation transport algorithms, computer hardware performance, and parallel computing make the clinical use of Monte Carlo based dose calculations possible. To compare the speed and accuracies of dose calculations between different developed codes, a benchmark tests were proposed at the XIIth ICCR (International Conference on the use of Computers in Radiation Therapy, Heidelberg, Germany 2000). A Monte Carlo treatment planning comprised of 28 various Intel Pentium CPUs was implemented for routine clinical use. The purpose of this study was to evaluate the performance of our system using the above benchmark tests. The benchmark procedures are comprised of three parts. a) speed of photon beams dose calculation inside a given phantom of 30.5 cm$\times$39.5 cm $\times$ 30 cm deep and filled with 5 ㎣ voxels within 2% statistical uncertainty. b) speed of electron beams dose calculation inside the same phantom as that of the photon beams. c) accuracy of photon and electron beam calculation inside heterogeneous slab phantom compared with the reference results of EGS4/PRESTA calculation. As results of the speed benchmark tests, it took 5.5 minutes to achieve less than 2% statistical uncertainty for 18 MV photon beams. Though the net calculation for electron beams was an order of faster than the photon beam, the overall calculation time was similar to that of photon beam case due to the overhead time to maintain parallel processing. Since our Monte Carlo code is EGSnrc, which is an improved version of EGS4, the accuracy tests of our system showed, as expected, very good agreement with the reference data. In conclusion, our Monte Carlo treatment planning system shows clinically meaningful results. Though other more efficient codes are developed such like MCDOSE and VMC＋＋, BEAMnrc based on EGSnrc code system may be used for routine clinical Monte Carlo treatment planning in conjunction with clustering technique.

• Proceedings of the Korean Society of Medical Physics Conference
• /
• 2002.09a
• /
• pp.113-115
• /
• 2002

The Monte Carlo simulation method is a numerical solution to a problem that models objects interacting with other objects or their environment based upon simple object-object or object-environment relationships. In spite of its great accuracy, It was turned away because of long calculation time to simulate a model. But, it is used to simulate a linear accelerator frequently with the advance of computer technology. To simulate linear accelerator in Monte Carlo simulations, there are many parameters needed to input to Monte Carlo code. These data can be supported by a linear accelerator manufacturer. Although the model of a linear accelerator is the same, a different characteristic property can be found. Thus, we performed a commissioning process of 6MV photon beam in Varian 2300C/D model with BEAM/EGS4 Monte Carlo code. The head geometry data were put into BEAM/EGS4 data. The mean energy and energy spread of the electron beam incident on the target were varied to match Monte Carlo simulations to measurements. TLDs (thermoluminescent dosimeter) and radiochromic films were employed to measure the absorbed dose in a water phantom. Beam profile was obtained in 40cm${\times}$40cm field size and Depth dose was in 10cm${\times}$10cm. At first, we compared the depth dose between measurements and Monte Carlo simulations varying the mean energy of an incident electron beam. Then, we compared the beam profile with adjusting the beam radius of the incident electron beam in Monte Carlo simulation. The results were found that the optimal mean energy was 6MV and beam radius of 0.1mm was well matched to measurements.

• Journal of Biosystems Engineering
• /
• v.38 no.2
• /
• pp.87-94
• /
• 2013

Purpose: Pineapple is now the third most important tropical fruit in world production after banana and citrus. Phytosanitary irradiation is recognized as a promising alternative treatment to chemical fumigation. However, most of the phytosanitary irradiation studies have dealt with physiochemical properties and its efficacy. Accurate dose calculation is crucial for ensuring proper process control in phytosanitary irradiation. The objective of this study was to optimize phytosanitary irradiation treatment of pineapple in various radiation sources using Monte Carlo simulation. Methods: 3-D geometry and component densities of the pineapple, extracted from CT scan data, were entered into a radiation transport Monte Carlo code (MCNP5) to obtain simulated dose distribution. Radiation energy used for simulation were 2 MeV (low-energy) and 10 MeV (high-energy) for electron beams, 1.25 MeV for gamma-rays, and 5 MeV for X-rays. Results: For low-energy electron beam simulation, electrons penetrated up to 0.75 cm from the pineapple skin, which is good for controlling insect eggs laid just below the fruit surface. For high-energy electron beam simulation, electrons penetrated up to 4.5 cm and the irradiation area occupied 60.2% of the whole area at single-side irradiation and 90.6% at double-side irradiation. For a single-side only gamma- and X-ray source simulation, the entire pineapple was irradiated and dose uniformity ratios (Dmax/Dmin) were 2.23 and 2.19, respectively. Even though both sources had all greater penetrating capability, the X-ray treatment is safer and the gamma-ray treatment is more widely used due to their availability. Conclusions: These results are invaluable for optimizing phytosanitary irradiation treatment planning of pineapple.

• Progress in Medical Physics
• /
• v.21 no.3
• /
• pp.304-310
• /
• 2010

Less execution of the electron arc treatment could in large part be attributed to the lack of an adequate planning system. Unlike most linear accelerators providing the electron arc mode, no commercial planning systems for the electron arc plan are available at this time. In this work, with the expectation that an easily accessible planning system could promote electron arc therapy, a commercial planning system was commissioned and evaluated for the electron arc plan. For the electron arc plan with use of a Varian 21-EX, Pinnacle3 (ver. 7.4f), with an electron pencil beam algorithm, was commissioned in which the arc consisted of multiple static fields with a fixed beam opening. Film dosimetry and point measurements were executed for the evaluation of the computation. Beam modeling was not satisfactory with the calculation of lateral profiles. Contrary to good agreement within 1% of the calculated and measured depth profiles, the calculated lateral profiles showed underestimation compared with measurements, such that the distance-to-agreement (DTA) was 5.1 mm at a 50% dose level for 6 MeV and 6.7 mm for 12 MeV with similar results for the measured depths. Point and film measurements for the humanoid phantom revealed that the delivered dose was more than the calculation by approximately 10%. The electron arc plan, based on the pencil beam algorithm, provides qualitative information for the dose distribution. Dose verification before the treatment should be mandatory.

• The Journal of Korean Society for Radiation Therapy
• /
• v.26 no.2
• /
• pp.163-170
• /
• 2014

Purpose : If electron, chosen for superficial oncotherapy, was applied with bolus, it could work as an important factor to a therapy result by showing a drastic change in surface dose. Hence the calculation value and the actual measurement value of surface dose of Treatment Planning System (TPS) according to four variables influencing surface dose when using bolus on an electron therapy were compared and analyzed in this paper. Materials and Methods : Four variables which frequently occur during the actual therapies (A: bolus thickness - 3, 5, 10 mm, B: field size - $6{\time}6$, $10{\time}10$, $15{\time}15cm2$, C: energy - 6, 9, 12 MeV, D: gantry angle - $0^{\circ}$, $15^{\circ}$) were set to compare the actual measurement value with TPS(Pinnacle 9.2, philips, USA). A computed tomography (lightspeed ultra 16, General Electric, USA) was performed using 16 cm-thick solid water phantom without bolus and total 54 beams where A, B, C, and D were combined after creating 3, 5 and 10 mm bolus on TPS were planned for a therapy. At this moment SSD 100 cm, 300 MU was investigated and measured twice repeatedly by placing it on iso-center by using EBT3 film(International Specialty Products, NJ, USA) to compare and analyze the actual measurement value and TPS. Measured film was analyzed with each average value and standard deviation value using digital flat bed scanner (Expression 10000XL, EPSON, USA) and dose density analyzing system (Complete Version 6.1, RIT, USA). Results : For the values according to the thickness of bolus, the actual measured values for 3, 5 and 10 mm were 101.41%, 99.58% and 101.28% higher respectively than the calculation values of TPS and the standard deviations were 0.0219, 0.0115 and 0.0190 respectively. The actual values according to the field size were $6{\time}6$, $10{\time}10$ and $15{\time}15cm2$ which were 99.63%, 101.40% and 101.24% higher respectively than the calculation values and the standard deviations were 0.0138, 0.0176 and 0.0220. The values according to energy were 6, 9, and 12 MeV which were 99.72%, 100.60% and 101.96% higher respectively and the standard deviations were 0.0200, 0.0160 and 0.0164. The actual measurement value according to beam angle were measured 100.45% and 101.07% higher at $0^{\circ}$ and $15^{\circ}$ respectively and standard deviations were 0.0199 and 0.0190 so they were measured 0.62% higher at $15^{\circ}$ than $0^{\circ}$. Conclusion : As a result of analyzing the calculation value of TPS and measurement value according to the used variables in this paper, the values calculated with TPS on 5 mm bolus, $6{\time}6cm2$ field size and low-energy electron at $0^{\circ}$ gantry angle were closer to the measured values, however, it showed a modest difference within the error bound of maximum 2%. If it was beyond the bounds of variables selected in this paper using electron and bolus simultaneously, the actual measurement value could differ from TPS according to each variable, therefore QA for the accurate surface dose would have to be performed.

• Progress in Medical Physics
• /
• v.33 no.4
• /
• pp.142-149
• /
• 2022

Purpose: This study seeks to compare the dosimetric parameters of the bulk electron density (ED) approach and synthetic computed tomography (CT) image in terms of position variation of the air cavity in magnetic resonance-guided radiotherapy (MRgRT) for patients with pancreatic cancer. Methods: This study included nine patients that previously received MRgRT and their simulation CT and magnetic resonance (MR) images were collected. Air cavities were manually delineated on simulation CT and MR images in the treatment planning system for each patient. The synthetic CT images were generated using the deep learning model trained in a prior study. Two more plans with identical beam parameters were recalculated with ED maps that were either manually overridden by the cavities or derived from the synthetic CT. Dose calculation accuracy was explored in terms of dose-volume histogram parameters and gamma analysis. Results: The D_95% averages were 48.80 Gy, 48.50 Gy, and 48.23 Gy for the original, manually assigned, and synthetic CT-based dose distributions, respectively. The greatest deviation was observed for one patient, whose D_95% to synthetic CT was 1.84 Gy higher than the original plan. Conclusions: The variation of the air cavity position in the gastrointestinal area affects the treatment dose calculation. Synthetic CT-based ED modification would be a significant option for shortening the time-consuming process and improving MRgRT treatment accuracy.

• Progress in Medical Physics
• /
• v.13 no.1
• /
• pp.21-26
• /
• 2002

A practical calculation algorithm which calculates the relative output factor(ROF) for irregular shaped electron field has been developed and evaluated the accuracy of the algorithm. The algorithm adapted two-source model, which assumes that the electron dose can be express as sum of the primary source component and the scattered component from the shielding block. Original two-source model has been modified in order to make the algorithm simpler and to reduce the number of parameters needed in the calculation, while the calculation error remains within clinical tolerance range. The primary source is assumed to have Gaussian distribution, while the scattered component follows the inverse square law. Depth and angular dependency of the primary and the scattered are ignored ROF can be calculated with three parameters such as, the effective source distance, the variance of primary source, and the scattering power of the block. The coefficients are obtained from the square shaped-block measurements and the algorithm is confirmed from the rectangular or irregular shaped-fields used in the clinic. The results showed less than 1.0 ％ difference between the calculation and measurements for most cases. None of cases which have bigger than 2.1 ％ have been found. By improving the algorithm for the aperture region which shows the largest error, the algorithm could be practically used in the clinic, since one can acquire the 1011 parameter's with minimum measurements(5∼6 measurements per cones) and generates accurate results within the clinically acceptable range.