Analytical parameterization of Bragg curves for proton beams in muscle, bone, and polymethylmethacrylate.

Proton dose calculation in media other than water may be of interest for either research purposes or clinical practice. Current study aims to quantify the required parameters for analytical proton dosimetry in muscle, bone, and PMMA. Required analytical dosimetry parameters were extracted from ICRU-49 report and Janni study. Geant4 Toolkit was also used for Bragg curve simulation inside the investigated media at different proton energies. Calculated and simulated dosimetry data were compared using gamma analysis. Simulated and calculated Bragg curves are consistent, a fact that confirms the validity of reported parameters for analytical proton dosimetry inside considered media. Furthermore, derived analytical parameters for these media are different from those of water. Listed parameters can be reliably utilized for analytical proton dosimetry inside muscle, bone, and PMMA. Furthermore, accurate proton dosimetry inside each medium demands dedicated analytical parameters and one is not allowed to use the water coefficients for non-water media.

Dynamic anthropomorphic thorax phantom for quality assurance of motion management in radiotherapy.

Background and purpose: Motion management techniques are important to spare the healthy tissue adequately. However, they are complex and need dedicated quality assurance. The aim of this study was to create a dynamic phantom designed for quality assurance and to replicate a patient's size, anatomy, and tissue density. Materials and methods: A computed tomography (CT) scan of a cancer patient was used to create molds for the lungs, heart, ribs, and vertebral column via additive manufacturing. A pump system and software were developed to simulate respiratory dynamics. The extent of respiratory motion was quantified using a 4DCT scan. End-to-end tests were conducted to evaluate two motion management techniques for lung stereotactic body radiotherapy (SBRT). Results: The chest wall moved between 4 mm and 13 mm anteriorly and 2 mm to 7 mm laterally during the breathing. The diaphragm exhibited superior-inferior movement ranging from 5 mm to 16 mm in the left lung and 10 mm to 36 mm in the right lung. The left lung tumor displaced ± 7 mm superior-inferiorly and anterior-posteriorly. The CT numbers were for lung: -716 ± 108 HU (phantom) and -713 ± 70 HU (patient); bone: 460 ± 20 HU (phantom) and 458 ± 206 HU (patient); soft tissue: 92 ± 9 HU (phantom) and 60 ± 25 HU (patient). The end-to-end testing showed an excellent agreement between the measured and the calculated dose for ion chamber and film dosimetry. Conclusions: The phantom is recommended for quality assurance, evaluating the institution's specific planning and motion management strategies either through end-to-end testing or as an external audit phantom.

Erratum to: Surface guided 3DCRT in deep-inspiration breath-hold for left sided breast cancer radiotherapy: implementation and first clinical experience in Iran.

[This corrects the article DOI: 10.5603/RPOR.a2022.0103.].

Application of Localized Surface Plasmon Resonance of Conjugated Gold Nanoparticles in Spectral Diagnosis of SARS-CoV-2: A Numerical Study.

Severe respiratory syndrome COVID-19 (SARS-CoV-2) outbreak has became the most important global health issue, and simultaneous efforts to fast and low-cost diagnosis of this virus were performed by researchers. One of the most usual tests was colorimetric methods based on the change of color of gold nanoparticles in the presence of viral antibodies, antigens, and other biological agents. This spectral change can be due to the aggregation of the particles or the shift of localized surface plasmon resonance due to the electrical interactions of surface agents. It is known that surface agents could easily shift the absorption peak of metallic nanocolloids which is attributed to the localized surface plasmon resonance. Experimental diagnosis assays for colorimetric detection of SARS-CoV-2 using Au NPs were reviewed, and the shift of absorption peak was studied from the viewpoint of numerical analysis. Using the numerical method, the refractive index and real and imaginary parts of the effective relative permittivity of the viral biological shell around Au NPs were obtained. This model gives a quantitative description of colorimetric assays of the detection of SARS-CoV-2 using Au NPs.

A dose planning study for cardiac and lung dose sparing techniques in left breast cancer radiotherapy: Can free breathing helical tomotherapy be considered as an alternative for deep inspiration breath hold?

Purpose: To investigate the possibility to be able to offer left sided breast cancer patients, not suitable for DIBH, an organ at risk saving treatment. Materials and Methods: Twenty patients receiving radiotherapy for left breast cancer in DIBH were enrolled in the study. Planning CT scans were acquired in the same supine treatment position in FB and DIBH. 3DCRT_DIBH plans were designed and optimized using two parallel opposed tangent beams (with some additional segments) for the breast and chest wall and anterior-posterior fields for regional lymph nodes irradiation. Additionally, FB helical tomotherapy plans were optimized to minimize heart and lung dose. All forty plans were optimized with at least 95% of the total CTV covered by the 95% of prescribed dose of 50 Gy in 25 fractions. Results: HT_FB plans showed significantly better dose homogeneity and conformity compared to the 3DCRT_DIBH specially for regional nodal irradiation. The heart mean dose was almost comparable in 3DCRT_DIBH and HT_FB while the volume (%) of the heart receiving 25 Gy had a statistically significant reduction from 7.90 ± 3.33 in 3DCRT_DIBH to 0.88 ± 0.66 in HT_FB. HT_FB was also more effective in left descending artery (LAD) mean dose reduction about 100% from 30.83 ± 9.2 Gy to 9.7 ± 3.1. The ipsilateral lung volume receiving 20 Gy has a further reduction of 43 % in HT_FB compared with 3DCRT_DIBH. For low dose comparison, 3DCRT_DIBH was superior for contralateral organ sparing compared to the HT_FB due to the limited angle for dose delivery. Conclusion: For patients who cannot be a candidate for DIBH for any reason, HT in free breathing may be a good alternative and provides heart and ipsilateral lung dose sparing, however with the cost of increased dose to contralateral breast and lung.

Performance evaluation of buildup bolus during external radiotherapy of mastectomy patients: treatment planning and film dosimetry.

A buildup bolus is used during the post-mastectomy radiotherapy (PMRT) to overcome under-dosage issues in the chest wall. The current study is aimed at evaluating the performance of a bolus in dose enhancement through both film dosimetry and treatment planning approaches. Twenty patients were enrolled in current research. The received dose by the skin at the lateral and medial regions of the chest wall in the presence and absence bolus was evaluated. Film dosimetry results showed that the presence of the bolus can averagely increase the skin dose by about 80% (P value < 0.001) and 92% (P value < 0.001) in lateral and medial regions, respectively. No significant difference was observed between the measured and treatment planning system (TPS)-calculated dose values in the presence of bolus. The presence of the bolus can considerably increase the absorbed dose by superficial chest wall regions. The TPS shows a favorable performance in superficial dose calculations in the presence of the buildup bolus. Hosseini et al.: demonstration of implemented research in the current study.

Surface guided 3DCRT in deep-inspiration breath-hold for left sided breast cancer radiotherapy: implementation and first clinical experience in Iran.

Background: The aim of the study is to evaluate the overall accuracy of the surface-guided radiotherapy (SGRT) workflow through a comprehensive commissioning and quality assurance procedures and assess the potential benefits of deep-inspiration breath-hold (DIBH) radiotherapy as a cardiac and lung dose reduction approach for left-sided breast cancer irradiation. Materials and methods: Accuracy and reproducibility of the optical surface scanner used for DIBH treatment were evaluated using different phantoms. Patient positioning accuracy and reproducibility of DIBH treatment were evaluated. Twenty patients were studied for treatment plan quality in target dose coverage and healthy organ sparing for the two different treatment techniques. Results: Reproducibility tests for the surface scanner showed good stability within 1 mm in all directions. The maximum position variation between applied shifts on the couch and the scanner measured offsets is 1 mm in all directions. The clinical study of 200 fractions showed good agreement between the surface scanner and portal imaging with the isocenter position deviation of less than 3 mm in each lateral, longitudinal, and vertical direction. The standard deviation of the DIBH level showed a value of < 2 mm during all evaluated DIBHs. Compared to the free breathing (FB) technique, DIBH showed significant reduction of 48% for heart mean dose, 43% for heart V25, and 20% for ipsilateral lung V20. Conclusion: Surface-guided radiotherapy can be regarded as an accurate tool for patient positioning and monitoring in breast radiotherapy. DIBH treatment are considered to be effective techniques in heart and ipsilateral lung dose reductions for left breast radiotherapy.

Comparing the inter-seed effect for some iodine-125 brachytherapy sources through a Monte Carlo simulation approach.

BACKGROUND AND OBJECTIVE: One of the principal deficits of TG-43(U1) dosimetry protocol is ignoring the attenuation effect of each seed during the multiseed brachytherapy implants. To take into account this shadowing effect, a parameter known as the inter-seed effect (ISE) has been introduced. Current study aims to evaluate and compare the ISE for some I-125 brachytherapy seeds through a Monte Carlo (MC) simulation approach. METHODS: Five different models of I-125 seeds including 6711, ProstaSeed, STM 1251, 3500, and IAI-125A were simulated by MCNPX MC Code. Validity of simulated seed models was confirmed through comparing the corresponding dosimetric parameters such as radial dose function and anisotropy function with those reported by relevant literature. Then, the ISE for each seed at different distances was determined in three seeds implant configuration. Additionally, the relevant attenuation factors (AF) were also determined for each brachytherapy source. RESULTS: The obtained results demonstrated that minimum inter-seed attenuation belongs to the ProstaSeed. In return, the maximum inter-seed attenuation was found for 3500 brachytherapy seed. The mean attenuation factor for 6711, ProstaSeed, STM 1251, 3500, and IAI-125A seed models was equal to 0.764, 0.829, 0.785, 0.594, and 0.772, respectively. CONCLUSIONS: It can be concluded that the considered I-125 brachytherapy seeds show different inter-seed effects and mean AF values for the same multiseed implant arrangement. This finding can be attributed to discrepancies in internal seed design and configuration.

Effective energy assessment during breast cancer intraoperative radiotherapy by low-energy X-rays: A Monte Carlo study.

The study reported in the present paper aimed to evaluate the effective energy (Eeff) of X-rays emitted from the surface of a bare X-ray probe and from different spherical applicators with various diameters, which are widely employed for low kV intraoperative radiotherapy (IORT) of breast cancer. A previously validated Monte Carlo model of the INTRABEAM system along with applicator diameters of 1.5-5 cm (with 0.5 cm increments) was employed for this purpose. The results show that the presence of the applicator can considerably harden the X-rays produced by the bare probe so that Eeff increases by a factor of about 2.6. Variations of applicator diameter also affects the X-ray effective energy. Specifically, increasing the applicator diameter from 1.5 to 3 cm and 3.5-5 cm resulted in an increase in the Eeff by 8.8% and 14.6%, respectively. The validity of the calculated Eeff values was confirmed by a reasonable agreement between the obtained probability density distributions (PDDs) for the full X-ray energy spectrum and those for the corresponding single effective energies, for different applicator diameters. The Eeff values obtained for different applicator diameters and the bare probe alone can be used as an alternative for the corresponding full energy spectra, in Monte Carlo-based dosimetry simulations of low-energy therapeutic X-rays, as well as for determining quality conversion factors of any ion chambers employed for low kV-IORT absolute dosimetry.

Impact of spherical applicator diameter on relative biologic effectiveness of low energy IORT X-rays: A hybrid Monte Carlo study.

INTRODUCTION: Low-kV IORT (Low kilovoltage intraoperative radiotherapy) using INTRABEAM machine and dedicated spherical applicators is a candidate modality for breast cancer treatment. The current study aims to quantify the RBE (relative biologic effectiveness) variations of emitted X-rays from the surface of different spherical applicators and bare probe through a hybrid Monte Carlo (MC) simulation approach. MATERIALS AND METHODS: A validated MC model of INTRABEAM machine and different applicator diameters, based on GEANT4 Toolkit, was employed for RBE evaluation. To doing so, scored X-ray energy spectra at the surface of each applicator diameter/bare probe were used to calculate the corresponding secondary electron energy spectra at various distances inside the water and breast tissue. Then, MCDS (Monte Carlo damage simulation) code was used to calculate the RBE values according to the calculated electron spectra. RESULTS: Presence of spherical applicators can increase the RBE of emitted X-rays from the bare probe by about 22.3%. In return, changing the applicator diameter has a minimal impact (about 3.2%) on RBE variation of emitted X-rays from each applicator surface. By increasing the distance from applicator surface, the RBE increments too, so that its value enhances by about 10% with moving from 2 to 10 mm distance. Calculated RBE values within the breast tissue were higher than those of water by about 4% maximum value. CONCLUSION: Ball section of spherical IORT applicators can affect the RBE value of the emitted X-rays from INTRABEAM machine. Increased RBE of breast tissue can reduce the prescribed dose for breast irradiation if INTRABEAM machine has been calibrated inside the water.

Optimization of an ultra-fast silicon detector for proton and carbon beams using GEANT4 Monte Carlo toolkit.

AIM: The purpose of this study was to investigate the crosstalk effects between adjacent pixels in a thin silicon detector with 50 um thickness. BACKGROUND: There are some limitations in the applications of detectors in hadron therapy. So it is necessary to have a detector with concurrent excellent time and resolution. In this work, the GEANT4 toolkit was applied to estimate the best value for energy cutoff in the thin silicon detector in order to optimize the detector. MATERIALS AND METHODS: GEANT4 toolkit was applied to simulate the transport and interactions of particles. Calculations were performed for a thin silicon detector (2 cm × 2 cm×0.005 cm) irradiated by proton and carbon ion beams. A two-dimensional array of silicon pixels in the x-y plane with 100 um × 100 um × 50 um dimensions build the whole detector. In the end, the ROOT package is used to interpret and analyze the results. RESULTS: It is seen that by the presence of energy cutoff, pixels with small deposited energy are ignored. The best values for energy cutoff are 0.01 MeV and 0.7 MeV for proton and carbon ion beams, respectively. By applying these energy cutoff values, efficiency and purity values are maximized and also minimum output errors are achieved. CONCLUSIONS: The results are reasonable, good and useful to optimize the geometry of future silicon detectors in order to be used as beam monitoring in hadron therapy applications.

Monte Carlo based analysis and evaluation of energy spectrum for low-kV IORT spherical applicators.

Low-kV IORT is an increasing modality for breast cancer treatment. Soft X-rays from INTRABEAM (Carl Zeiss Meditec AG, Oberkochen, Germany), a dedicated IORT device along with special spherical applicators are employed for this purpose. A Monte Carlo model of INTRAMBEAM and spherical applicators are introduced in the current study to evaluate the dosimetric and physical characteristics of emitted X-rays from the bare probe and different applicator diameters. X-ray probe and different applicator diameters of 1.5cm to 5cm were simulated by GEANT4 Monte Carlo Toolkit. Then, the validity of the simulated model was evaluated by comparing the Monte Carlo based PDD (percentage depth dose) and anisotropy data with those reported by the manufacturer. Finally, the physical characteristics of X-rays such as the mean and most probable energy as well as the LET of secondary electrons were obtained and analyzed. There was a good agreement between the simulated and reported PDDs for bare probe and different applicator diameters. The anisotropy values were also within an exemplary range reported by the manufacturer. The X-ray mean energy shifts from 25.6keV to 28.6keV with variations of applicator diameter. The maximum variation of the secondary electron LET was 9% with changing the applicator diameter. The usefulness of GEANT4 Toolkit for Monte Carlo based commissioning of INTRABEAM machine was confirmed in the current study. The variations of the applicator diameter slightly changed the physical characteristics of low-kV X-rays and LET of secondary electrons which cannot considerably affect the relative biologic effectiveness (RBE) for different applicator diameters.

A Monte Carlo study on dose perturbation due to dental restorations in a 15 MV photon beam.

AIM: The main purpose of this study is to evaluate the effect of dose perturbation due to common dental restoration materials in the head and neck radiotherapy with a 15 MV external photon beam. SETTING AND DESIGN: Teeth with three dental restorations such as tooth filled with Amalgam, Ni-Cr alloy, and Ceramco were simulated by MCNPX Monte Carlo code. In this simulation, the dental materials were exposed by a 15 MV photon beam from a Siemens Primus linac, inside a water phantom. MATERIALS AND METHODS: A Siemens Primus linear accelerator and a phantom including: tooth only, tooth with Amalgam, tooth with Ni-Cr alloy, and tooth with Ceramco were simulated by MCNPX Monte Carlo code, separately. The percentage dose change was evaluated relative to dose in water versus depth for these samples on the beam's central axis. The absolute dose by prescription of 100 cGy dose in water phantom at 3.0 cm depth was calculated for water, tooth, tooth with Amalgam, tooth with Ni-Cr alloy, and tooth with Ceramco. RESULTS: The maximum percentage dose change is related to tooth with Ni-Cr alloy, tooth, tooth with Ceramco, and tooth with Amalgam with amounts of 7.73%, 6.95%, 4.7%, and 3.06% relative to water at 0.75 cm depth, respectively. When 100.0 cGy dose was prescribed at 3.1 cm, the maximum absolute dose was 201.0% in the presence of tooth with Ni-Cr alloy at 0.75 cm. CONCLUSION: Introduction of the compositions of dental restorations can improve the accuracy of dosimetric calculations in treatment planning and protect the healthy tissues surrounding teeth from a considerable overdose.

Experimental study of the influence of dental restorations on thermal and fast photo-neutron production in radiotherapy with a high-energy photon beam.

In head and neck radiation therapy, the presence of dental restorations can increase unwanted neutron dose to the patient. This study aimed at the measurement of secondary neutron production induced by irradiation of a healthy tooth, Amalgam, Ni-Cr alloy and Ceramco with a photon beam generated in the treatment head of a Siemens Primus linac at a voltage of 15 MV. The irradiation field amounted to 10 × 10 cm2. The measurements of thermal and fast-neutron equivalent doses were performed by means of CR-39 detectors positioned in various depths of a Perspex (polymethyl methacrylate) phantom as at open field as at presence of corresponding dental restorations. The general trend of thermal neutron as well as fast-neutron equivalent dose behind the denture samples reveals their reduction with increasing depth. The maximum values of thermal-neutron dose related to Amalgam, Ceramco and Ni-Cr alloy amount to 1.45 mSv/100 MU, 1.38 mSv/100 MU and 1.32 mSv/100 MU, whereas the corresponding maximum values of fast-neutron dose at the depth of 1.8 cm amount to 0.19 mSv/100 MU, 1.04 mSv/100 MU and 0.97 mSv/100 MU, respectively. The present study investigates the neutron dose accompanied with radiotherapy. It is recommended that attempts have to be made to ensure that dental restorations are not in the path of the primary high-energy photon beam. Considering treatment planning, the guidelines of radiation protection should be improved.

Monte Carlo based determination of dose distribution for some patch sources employed for radionuclide skin therapy.

Several methods are introduced for skin cancer treatment. An encouraging method is radionuclide skin therapy, where high-energy beta emitting radionuclides such as 32P, 90Y, 188Re and 166Ho are employed for skin irradiation. This study aimed to calculate the dose distribution for mentioned radionuclides at different layers of skin phantom through Monte Carlo simulation. Depth dose distribution, transverse dose profile and isodose curves related to the patch sources under investigation were calculated by MCNPX code. All of calculations were performed inside a developed skin phantom. Obtained results were compared with those reported by other studies to evaluate the validity of simulations. The results showed that the 166Ho and 32P have steeper dose gradient within the depth which can lead to the better normal tissue sparing. Simulated depth dose distributions were in a good agreement with other published studies and confirmed the validity of performed simulations. The obtained transverse dose profiles at 0.2 mm depth had acceptable symmetry and flatness that can improve the dose uniformity within the target area. Calculated isodose curves showed that the 90% isodose level covers a circular area with the diameter of around 8 mm for all studied beta sources. From the results, it can be concluded that 166Ho and 32P are more effective in treatment of superficial skin lesions, while, 90Y and 188Re are more rational choices in treatment of deeply distributed skin tumors. Size of employed patch source should be based on the target area to minimize the delivered dose to the adjacent tissues.

Effect of various dental restorations on dose distribution of 6 MV photon beam.

AIM: The purpose of this study is to evaluate the effect of various dental restoration materials on dose distribution in radiotherapy of head and neck cancer with 6 MV photon beam of a medical linac. SETTING AND DESIGN: The dental restorations include tooth, tooth with amalgam, tooth with Ni-Cr alloy, and tooth with Ceramco. Dose perturbation due to the dental restorations on 6 MV beam of Siemens Primus linac was calculated by MCNPX Monte Carlo code. These dental materials were separately simulated in a cubic water phantom. MATERIALS AND METHODS: Photon percentage dose change in the presence of tooth, tooth with amalgam, tooth with Ni-Cr alloy, and tooth with Ceramco was calculated at various depths on the central axis of the beam relative to the dose in water. In another evaluation, the absolute dose (cGy) for water, tooth, tooth with amalgam, tooth with Ni-Cr alloy, and tooth with Ceramco was obtained by calculating 100 cGy dose at 0.75 cm depth in the water phantom. RESULTS: Based on the calculations performed, maximum percentage dose change due to backscattering was 11%, 8%, 6%, and 4% for amalgam, Ni-Cr alloy, Ceramco and tooth at depth of 0.75 cm, respectively. The maximum dose perturbation by amalgam is due to its higher mass density and atomic number. CONCLUSION: Ignoring the effect of dental restoration compositions on dose distribution causes discrepancy in clinical treatment planning system calculations.

Simulation of positron emitters for monitoring of dose distribution in proton therapy.

AIM: The purpose of this work was to estimate the dependency between the produced positron emitters and the proton dose distribution as well as the dependency between points of annihilation and the proton dose distribution. BACKGROUND: One important feature of proton therapy is that, through the non-elastic nuclear interaction of protons with the target nuclei such as 12C, 14N and 16O, it produces a small number of positron-emitting radioisotopes along the beam-path. These radioisotopes allow imaging the Bragg peak position which is related to the proton dose distribution by using positron emission tomography. METHODS: In this study, the GEANT4 toolkit was applied to simulate a soft and bone tissue phantom in proton therapy to evaluate the positron emitter productions and the actual annihilation points of ß+. Simulation was done by delivering pencil and spread-out Bragg peak (SOBP) proton beams. RESULTS: The findings showed that (15O, 11C, 13N) and (11C, 15O, 38K, 30P, 39Ca, 13N) are the most suitable positron emitters in the soft and bone tissue respectively. By increasing the proton energy, the distance between the peak of annihilation profile and Bragg peak is almost constant, but the distance between the Bragg peak position and positron annihilation point peak in bone tissue is smaller than that in the soft tissue. The peak of ß+ activity distribution becomes sharper at higher proton energies. CONCLUSIONS: There is a good relationship between the positions of positron annihilation profile and positron emitters radioactive decay. Also, GEANT4 is a powerful and suitable tool for simulation of nuclear interactions and positron emitters in tissues.

Assessment of Mean Glandular Dose in Mammography System with Different Anode-Filter Combinations Using MCNP Code.

BACKGROUND: X-ray mammography is one of the general methods for early detection of breast cancer. Since glandular tissue in the breast is sensitive to radiation and it increases the risk of cancer, the given dose to the patient is very important in mammography. OBJECTIVES: The aim of this study was to determine the average absorbed dose of X-ray radiation in the glandular tissue of the breast during mammography examinations as well as investigating factors that influence the mean glandular dose (MGD). One of the precise methods for determination of MGD absorbed by the breast is Monte Carlo simulation method which is widely used to assess the dose. MATERIALS AND METHODS: We studied some different X-ray sources and exposure factors that affect the MGD. "Midi-future" digital mammography system with amorphous-selenium detector was simulated using the Monte Carlo N-particle extended (MCNPX) code. Different anode/filter combinations such as tungsten/silver (W/Ag), tungsten/rhodium (W/Rh), and rhodium/aluminium (Rh/Al) were simulated in this study. The voltage of X-ray tube ranged from 24 kV to 32 kV with 2 kV intervals and the breast phantom thickness ranged from 3 to 8 cm, and glandular fraction g varied from 10% to 100%. RESULTS: MGD was measured for different anode/filter combinations and the effects of changing tube voltage, phantom thickness, combination and glandular breast tissue on MGD were studied. As glandular g and X-ray tube voltage increased, the breast dose increased too, and the increase of breast phantom thickness led to the decrease of MGD. The obtained results for MGD were consistent with the result of Boone et al. that was previously reported. CONCLUSION: By comparing the results, we saw that W/Rh anode/filter combination is the best choice in breast mammography imaging because of the lowest delivered dose in comparison with W/Ag and Rh/Al. Moreover, breast thickness and g value have significant effects on MGD.

Ion therapy for uveal melanoma in new human eye phantom based on GEANT4 toolkit.

Radiotherapy with ion beams like proton and carbon has been used for treatment of eye uveal melanoma for many years. In this research, we have developed a new phantom of human eye for Monte Carlo simulation of tumors treatment to use in GEANT4 toolkit. Total depth-dose profiles for the proton, alpha, and carbon incident beams with the same ranges have been calculated in the phantom. Moreover, the deposited energy of the secondary particles for each of the primary beams is calculated. The dose curves are compared for 47.8MeV proton, 190.1MeV alpha, and 1060MeV carbon ions that have the same range in the target region reaching to the center of tumor. The passively scattered spread-out Bragg peak (SOBP) for each incident beam as well as the flux curves of the secondary particles including neutron, gamma, and positron has been calculated and compared for the primary beams. The high sharpness of carbon beam׳s Bragg peak with low lateral broadening is the benefit of this beam in hadrontherapy but it has disadvantages of dose leakage in the tail after its Bragg peak and high intensity of neutron production. However, proton beam, which has a good conformation with tumor shape owing to the beam broadening caused by scattering, can be a good choice for the large-size tumors.

Evaluation of energy deposition and secondary particle production in proton therapy of brain using a slab head phantom.

AIM: Evaluation of energy deposition of protons in human brain and calculation of the secondary neutrons and photons produced by protons in proton therapy. BACKGROUND: Radiation therapy is one of the main methods of treating localized cancer tumors. The use of high energy proton beam in radiotherapy was proposed almost 60 years ago. In recent years, there has been a revival of interest in this subject in the context of radiation therapy. High energy protons suffer little angular deflection and have a well-defined penetration range, with a sharp increase in the energy loss at the end of their trajectories, namely the Bragg peak. MATERIALS AND METHODS: A slab head phantom was used for the purpose of simulating proton therapy in brain tissue. In this study simulation was carried out using the Monte Carlo MCNPX code. RESULTS: By using mono energetic proton pencil beams, energy depositions in tissues, especially inside the brain, as well as estimating the neutron and photon production as a result of proton interactions in the body, together with their energy spectra, were calculated or obtained. The amount of energy escaped from the head by secondary neutrons and photons was determined. CONCLUSIONS: It was found that for high energy proton beams the amount of escaped energy by neutrons is almost 10 times larger than that by photons. We estimated that at 110 MeV beam energy, the overall proton energy "leaked" from the head by secondary photons and neutrons to be around 1%.