
 cancer is a major social problem , and it is the main cause of death between the ages 4565 years . in the treatment of cancer ,
hadron therapy , thus , has great prospects for being used in early stages of tumor disease not amenable to surgery .
light particles are often used for radiation therapy because they have a well - defined penetration in tissue , the depth being dependent on the incident energy of the particles and the nature of the irradiated tissue .
the other main reason for using these particles in radiation therapy lies in their physical , the ability to deliver the dose to the interest target .
beams of charged particles have specific dose distribution , exhibiting a flat entrance dose relatively ( plateau ) , followed by a sharp dose peak , the bragg peak , in which the particles lose rest of their energy . due to their physical and radiobiological properties ,
they allow one to obtain a more conformal treatment , sparing better the healthy tissues located in proximity of the tumor and allowing a higher control of the disease .
indeed , its relative biological effectiveness ( rbe ) does depend on the linear energy transfer ( let ) and can become substantially large at the falling edge of the bragg peak [ 3 , 4 ] .
treatment planning in radiation therapy uses mathematical and physical formalisms to optimize the delivering a high and conformal dose to the target and limiting the doses to critical structures .
the dose tolerance levels for critical structures , as well as the required doses for various tumor types , are typically defined on the basis of decades of clinical experience .
then by proper choosing of different types of particles like h and the other light particles can be reduced the relevant dose to critical organs .
it is well recognized that h are extremely valuable to treat tumors close to critical structures such head and neck , brain stem , prostate , spinal cord , eye   or optic nerves . 
 

the main types of radiation therapy induced fatalities that have been widely reported are secondary cancers .
so , the most efficient way to prevent these secondary cancers is reducing the amount of dose scattered to the internal organs ; for example choosing a radiation technique or proper particle that minimizes the scattered dose to the other organs . in the regards to the secondary cancers ,
a recent review showed that secondary tumors occur more frequently in organs that are close to radiation fields , in the high / intermediate dose areas , and that is important to evaluate the scattered dose to those the internal organs along with their secondary cancer .
all these light particles deliver different levels of scattered dose to the internal organs and hence may induce different risks of secondary cancers .
the aim of this study is to evaluate the amount of the scattered dose to the internal organs situated in the intermediate / high dose region including the spine and nerves , etc . 
 

recently , orecchia has been reported that the potential indications for proton therapy to treatment of thyroid cancer is about 50 patients per year in the national center for oncological hadron therapy in pavia , italy .
taheri has investigated the potential advantages of intensity modulated h therapy ( impt ) compared with intensity modulated radiation therapy ( imrt ) in nasopharyngeal carcinoma .
they found that impt plans reduced the averaged mean dose of the thyroid gland by a factor of two related to imrt . 
 
 in this study , we have calculated the light particles energy deposition in the thyroid gland , for different particles such h , h , he , and he by mcnpx monte carlo code , and estimation of the total photon and neutron production due to interaction of incident particles by tissue , that is very important in order to evaluate the risk of secondary cancers .
one of the major concerns in hadron therapy is due to the neutrons and photons . because of their long rang
, they can damage the critical organs and increase the probability the secondary cancers . 


 mcnpx is a general purpose monte carlo radiation transport code designed to track many particles over broad ranges of energies .
the code deals with transport of more than 40 particles , and coupled transport , such transport of secondary gamma rays and neutrons resulting from different interactions in hadron therapy .
it is important that a dose calculation algorithm to be able to model complex and heterogeneous density mediums , as well as , calculating of the primary and secondary particles dose .
a cylindrical geometry was used to describe the neck phantom , having 6 cm of length and 6 cm of diameter .
the layers of the phantom include 0.12 cm skin , 0.4 cm adipose and 0.6 cm skeletal muscle , after these layers we consider the thyroid gland as cylindrical with 0.6 cm of diameter , vertebrae and spine with 0.4 cm diameters , respectively .
figure 1 shows the geometry of the neck phantom and we derived their material compositions and densities from icrp publications .
four light particles including h , h , he and he are used in our simulation .
geometry of the neck phantom including layers from left to right : 0.12 cm skin , 0.4 adipose , 0.6 cm skeletal muscle , 0.6 thyroid , 1 cm vertebrae , 0.4 spine .

 we have calculated the energy deposition variation with depth for the same range particles to compare their bragg peaks in the target organ . as well as ,
the suitable energy interval for scanning of the thyroid organ with step of 1 mev have been obtained : 37 - 53 mev for h , 50 - 71 mev for h , 128 - 190 mev for he and finally 146 - 214 mev for he .
light particles slow down in tissue , mainly through myriad coulomb interaction with electrons of atoms . as it is well known , the dominant interactions in tissue are ionization and excitation processes in more 99.9% of energy loss and the rest is belong to the nuclear interactions . in the treatment plan of hadron therapy ,
the energy of light particles with the same range of 50 mev for h ( 2.14 cm ) in thyroid organ are , 68 mev , 176mev , and 200 mev for h , he , he , respectively .
figure 2 shows the variation energy deposition with depth in the neck phantom for these particles . 
 
variation of the energy deposition with depth in the neck phantom for any particles with the same range is shown in figure 2 .
the sharper peak with less width , have more concentration for adjusting beam on the tumor and also beam energy spread .
it is well known that the fwhm and peak to plateau ratio ( ptpr ) are two significant factors that can be studied in treatment planning .
these results have been listed in table 1 and show that he with the highest value for ptpr and the minimum value for fwhm factor is a good choice . according to the figure 2 , with comparing the particles with same range can be understand the energy deposition in peak rises by increasing of the particle mass and atomic number . 
 
peak to plateau ratio and fwhm . 
 in hadron therapy of thyroid gland , vertebrae , spine and nerve can be considered as the critical organs which maybe damaged by the exposure of neutrons and photons secondary particles .
we have evaluated the total flux of the secondary particles as recorded in table 2 . 
  total flux for photons and neutrons secondary particles per one particle of beam . 

as it is well known , presence of neutron and photon secondary particles in hadron therapy , a fraction of dose can be deposited inside the thyroid gland and other far organs .
tables 3 and 4 indicate that how much energy fractions of neutrons and photons deposited in the tumor and other important healthy organs in the phantom .
the obtained results in table 3 and 4 indicated that secondary particles of h beam deposits maximum percentage dose in thyroid gland but the maximum damage for other sensitive organs is related to the h beam .
the results illustrated that we have the maximum absorption in thyroid region and then in spine and vertebrae , respectively . 
 
( a ) dose curves for neutron in sensitive organs and ( b ) dose curves for photons in sensitive organs .
consequently , when assessing the impact of neutron and photon doses , not only the absorbed dose is important but also the neutron and photon energy distribution spectra are essential in radiobiological effects evaluation .
the secondary neutrons may have energies up to the primary particle beams energy . based on the hadroinc model , neutrons with energy in excess of 10 mev produced by an intranuclear cascade processes are mainly forward - peaked and below 10 mev neutrons produced by an evaporation process , and are emitted more isotropically . 

figure 4 shows the simulated energy distribution of neutrons entering the neck by light particle beams .
the results show that neutrons produced by h interactions in the patient body have an average energy lower than the neutrons produced by other hardron beams with also wider energy distribution .
because of neutron elastic scattering in the soft tissue , there is a prevailing field of low energy .
furthermore of low energy neutron dose , a large fraction of the neutron dose is induced by fast neutrons . 
 
neutron spectra produced by h , h , he , and he beams . 
 for low energies neutrons ,
the most likely interaction in the body , because of the presence of hydrogen , is elastic scattering . in the low / thermal energy region
, there is a decreasing probability of neutrons slowing down , as low - energy neutrons sparsely interact with the body s material . as a result ,
an extensive part of the patient s body may be exposed to the secondary photon and neutron radiation fields .
the small amount of dose from uncharged particles - neutrons and photons- measurement or simulation is difficult and may be time consuming . in this research ,
two types of have been done : first , calculation of the secondary particle flux according to the energy of secondary particle , and the second is evaluation of the secondary particle flux according to the primary particle energy .
the neutron and photon energy spectra have been illustrated in figures 4 and 5 , respectively . 
 
spectra of photon flux for ( a ) h , ( b ) h , ( c ) he , and ( d ) he . 

the neutron spectra for these particles have significant high intensity smooth peak in low energy region .
thermal neutrons have a different and often much larger effective neutron absorption cross section for a nuclides in tissue than fast neutrons , and can therefore often be absorbed more easily by an atomic nucleus , creating a heavier , often unstable isotope as a neutron activation process .
after the peak , there is a continuum flat that the intensity is decreased smoothly by increase in neutron energy .
the photons flux spectra are important remarkable quantity need to study , deeply . according to the results of simulation
, there are some peaks which related to elements that have been specified in figure 5 .
some peaks are known very well , 4.43 mev from c excitation , 5.22 from cl , 6.13 and 6.92 from o. 
 
 total neutron and photon production for different energies of h , h , he , and he primary beams have been illustrated in figures 6 - 9 .
photon production result is fitted by polynomial third order function for h energy , whereas for other beams it is described by a linear function .
the figures show that the neutron total flux increased for he , he , h and h , respectively .
the total photon flux production due to h , he , he and h are increasing . 
 
variation of ( a ) neutron and ( b ) photon flux in terms of h energy .
variation of ( a ) neutron and ( b ) photon flux in terms of h energy .
variation of ( a ) neutron and ( b ) photon flux in terms of he energy .
variation of ( a ) neutron and ( b ) photon flux in terms of he energy .

 in this research , we simulated simplified neck phantom involving a tumor placed in 2.14 cm depth with mcnpx monte carlo code for h , h , he and he .
the results show that by increasing the atomic number , the bragg peak becomes sharper and the incident beam can deposit more energy in tumor . as mentioned before ,
one of the important factor in external therapy is secondary particles that produced by nuclear interaction with tissue .
we evaluated the flux and energy deposition of the secondary particles including neutron and photon in critical organs such spine and vertebrae placed after thyroid gland .
based on the dose of secondary particles , it can be concluded that h and he are the best therapy choices for thyroid glands whereas h is the worse particle . 

 
 background : hadron therapy is a novel technique of cancer radiation therapy which employs charged particles beams , 1h and light ions in particular . due to their physical and radiobiological properties , they allow one to obtain a more conformal treatment , sparing better the healthy tissues located in proximity of the tumor and allowing a higher control of the disease . 

 objective : as it is well known , these light particles can interact with nuclei in the tissue , and produce the different secondary particles such as neutron and photon . these particles can damage specially the critical organs behind of thyroid gland . 

 
 methods : in this research , we simulated neck geometry by mcnpx code and calculated the light particles dose at distance of 2.14 cm in thyroid gland , for different particles beam : 1h , 2h , 3he , and 4he . 
 thyroid treatment is important because the spine and vertebrae is situated right behind to the thyroid gland on the posterior side . 

 
 results : the results show that 2h has the most total flux for photon and neutron , 1.944e-3 and 1.7666e-2 , respectively . whereas 1h and 3he have best conditions , 8.88609e-4 and 1.35431e-3 for photon , 4.90506e-4 and 4.34057e-3 for neutron , respectively . 
 the same calculation has obtained for energy depositions for these particles . 

 
 conclusion : in this research , we investigated that which of these light particles can deliver the maximum dose to the normal tissues and the minimum dose to the tumor . by comparing these results for the mentioned light particles , 
 we find out 1h and 3he is the best therapy choices for thyroid glands whereas 2h is the worst . 
 
Introduction
Material and Methods
Results and Discussion
Conclusion
None