MEMBER'S PAPER




Measurement and Calculation of Neutron and Gamma Absorbed Dose in a Head Phantom in Boron Neutron Caputure Theory

A Pazirandeh and M E Fazilat1

Physics Department, University of Tehran
PO Box 1943,  Tehran  19395,  Iran

email:  paziran@chamran.ut.ac.ir

1Tarbiat Moddaress University, Tehran, Iran

ABSTRACT

In Boron Neutron Capture Therapy, neutron and gamma dose determination is of prime importance. The neutron and gamma absorbed doses were measured in a cylindrical head phantom, 18cm in diameter and 20cm in height. For this purpose gold foils (bare and cadmium covered) and TLD-700 were used for neutron flux and gamma dose measurements respectively. The contributions to absorbed dose were due to 10B(n,a)7Li, 1H(n,g)2D, 14N(n,p)14C and recoil proton reactions. The KERMA factors for different energy groups were used to convert the neutron flux to the relevant absorbed dose. There was good agreement between experimental values and the calculated dose.

  Keywords:       Boron Neutron Capture Therapy, Phantom, TLD, Dose components, Spectrum unfolding, Neutron Filter.

INTRODUCTION

Pioneering work of boron neutron capture therapy (BNCT) as a  potential, clinical technique to destroy malignant tumors - especially glioblastoma- is usually attributed to two individuals:  William H. Sweet and Hiroshi Hatanaka. Their work started in the early 1950s. In this technique, a neutron beam from a research reactor core passing through a proper moderator and a filter reduce the core neutron energy to epithermal range and a collimator converts them to an epithermal neutron beam. Prior to neutron bombardment,  the patient is administered a boron-10 compound such as BSH (Na12B12 H11 SH). Then the epithermal neutron beam impinges on the patient head. The epithermal neutron beam is normally used for deep-seated brain tumors since the epithermal neutrons are more or less slowed down and turn to thermal neutrons before reaching the tumor.  It is worth mentioning that, because of its anti-body properties, BSH overwhelms the tumor.  BSH is a stable and non-toxic compound.1 Therefore it can be applied without hesitation. The neutrons during passage through the brain undergo several reactions including 10B(n,a)7Li, 1H(n,g)2D, 14N(n,p)14C and elastic scattering from protons. Because the cross section for 10B(n,a)7Li reaction is rather high,3838 barns, Boron is very suitable for tumor therapy. In this reaction 2.79 MeV is released in the  form of kinetic energy of two ions, 7Li and 4He,  in a proportion reciprocal to their atomic masses. It should be pointed out that in 94% of reactions, a gamma ray with 478KeV is emitted. The energy of the alpha particles released in this reaction is enough to travel about 10mm in tissue. Therefore, each alpha particle because of its high LET can destroy a cancer cell. Because of the high Boron-10 concentration in the tumor, 40-50ppm, as compared to 10ppm in other parts of the brain, the absorbed dose due to this reaction in the tumor is several times higher than that in normal tissue. This is enhanced by the  peak thermal neutron flux in the tumor. 

  Success in BNCT lends itself in precise dose determination and it plays a key role in this technique. Since direct dose measurement in a human head is not practical, a head phantom is normally used instead. Ordinary distilled water is a suitable material that has very similar properties to normal brain tissue with regard to neutron interaction and energy absorption.  Different shaped head phantoms namely cylindrical, elliptical or cubic are used. We used a cylindrically-shaped phantom, 18cm in diameter and 20cm high filled with distilled water. The measurements were carried out in the head phantom containing a small tumor-shaped cylinder containing boric acid, 50ppm B-10. The phantom was placed in front of the neutron beam tube of  the Tehran Research Reactor (TRR) in such a way as to receive the emerging epithermal neutrons from the collimator. The beam was directed at the small cylinder.  The gamma ray absorbed dose component due to 10B(n,a)7Li, 1H(n,g)2D, and 14N(n,p)14C reactions and core gamma rays was measured using TLD-700 and the neutron flux was measured using gold foils. To convert neutron flux to absorbed dose the relevant KERMA factors were applied2.

  MATERIALS AND METHODS

  In 1993 there was an attempt to investigate the possibility of using the BNCT technique at the Tehran Research Reactor3 .  A collimator and an Al-Fe filter were designed, constructed and installed in a beam tube of TRR, see Fig. 1. The filter was composed of 24 pieces of alternating aluminum and iron blocks, 5cm each. To stop core gamma rays from reaching the target, a 15cm bismuth block was placed behind the filter.  

Fig. 1  Collimator and Al-Fe filter detailed design

To determine the central beam hole axis, a radiographic cassette containing an x-ray film was placed in front of the beam hole, once with a gadolinium converter and second time without a converter. The irradiation period extended 24 hours. The x-ray film was scanned by a laser-scanner to determine the center of the neutron beam.

The leakage neutron spectrum was measured at the beam outlet using an He-3 counter and the unfolding technique3,4. The spectrum from a bare Cf-252 neutron source, a distribution of well thermalized neutrons and the BNCT beam port neutron spectrum are shown in Fig. 2. 

The dose measurements were carried out in a cylindrical head phantom made of Plexiglas, filled with distilled water.  Bare and cadmium covered gold foils, 1cm diameter and 5mm, were used in three radial rows to measure thermal and epicadmium fluxes.  The irradiated gold foils were counted on a HPGe detector. The absolute net area under the prominent Au-198 photo-peak, 412Kev was determined utilizing a detector efficiency factor.    Using the following relations, thermal and epithermal neutron fluxes were determined.5

                                                       (1)

 

                                                                  (2)

        where

M      molecular weight of Au, u,

Na   =   Avogadro Number,

As     saturated activity, dps,

m       gold foil weight, g,

FCd   = cadmium self-absorption coefficient6,

Ac   =   cadmium covered gold foil activity, dps,

sth   =   thermal neutron capture cross section, 99b,

RI    =   integral resonance of Au-197, 1550 b,

Fres     resonance self-absorption coefficient6,

e             Detector efficiency.

Fig.3  shows thermal and epicadmium fluxes along three lines inside the head phantom.

For  gamma ray  dose  measurements the TLD-700s, 4mm in diameter and 1.5mm thick were used.    The TLDs were calibrated prior to irradiation.

DOSE DETERMINATION

Major reactions in BNCT treatment are due to 10B(n,a)7Li, 14N(n,p)14C , 1H(n,g)2D, recoil protons and (n,n’) reactions. The absorbed dose in brain is as a result of thermal and epithermal neutron reactions as follows:

1.      10B(n,a)7Li reaction

         10B + n à  a + 7Li  +2.79Mev        6%

         10B + n à  a + 7Li  +2.31Mev      94%

                     7Li à g + 0.478 Mev

 

2.    14N(n,p)14C  reaction

         14N + n à p + 14C + 0.62Mev

In this reaction 0.580 Mev energy is given off as kinetic energy of the proton and the remaining 40Kev, is taken away by 14C as kinetic energy.

3.  1H(n,g)2D reaction.

     Following neutron capture by a hydrogen nucleus, a 2.25 Mev gamma ray is emitted. Due to the high proportion of hydrogen in the  brain, absorbed dose of this reaction is rather high. TLD-700 chips measured the absorbed dose of gamma rays of these reactions and core gamma rays. Ten TLD-700 chips were used in two radial rows. They were irradiated for 24 hours in the head phantom. Table 1 shows measured dose along the path to the small cylinder in the phantom.

4.  Recoil proton and (n,n’) reactions.

     The determination of absorbed dose of epithermal neutrons is rather complicated. Scattering of epithermal neutrons from hydrogen nuclei resulted in proton recoil or elastic scattering. To estimate the absorbed dose due to these two reactions, the slowing down neutron spectrum was determined pointwise using the transport equation with anisotropic scattering. 

      We employed the ANISN/PC code along with a  17-group cross section library [IRAN.LIB7].  We calculated from the computational output the space dependant spectrum along the neutron trajectory towards the tumor.  At each point inside the head phantom, the absorbed dose of ith energy group was calculated from equations (3) and (4) using calculated fluxes by ANISN/PC.

                                                                       (3)
  where

Fn(Ei) =   KERMA (Kinetic Energy Released in
                  Medium) factor for ith group

f(Ei)  =   Group neutron flux, n/cm2.s

The KERMA factor for ith group, Fn (Ei) was calculated using the following relation8:

                                  Fn(Ei)  =  1.602*10-8 sN m-1Etr                    (4)

      where

 

Fn (Ei)  =   KERMA factor as a function of neutron energy, in J/kg or Grey

s            neutron-isotope cross section reaction for specified neutron energy,b

          Number of atoms in the sample

          mass of sample, g

Etr          total kinetic energy transferred to medium.

The incident epithermal neutrons on the human head phantom, undergo elastic scattering from hydrogen nuclei which may lead to recoil protons. As a result, the neutron spectrum extends to the thermal region as they travel towards the small cylinder tumor inside the phantom.


RESULTS AND DISCUSSIONS

  Table 2 and Fig. 4  give the total and the contribution of each component in building up RBE absorbed dose rate in the tumor along its central line. The result is in good agreement with those of Raaijimakers et.al.9. In spite of the fact that the neutron flux level is too low to be feasible in the clinical application. However, this facility can be utilized for in-vitro and in-vivo experiments. The results of the project, especially the neutron beam energy and dose distribution in tumor and surrounding tissues encourages us for future planning. The data obtained from this project has provided enough information to extend the project to practical application.

ACKNOWLEDGEMENT

Thanks are due to Tehran Research Reactor (TRR) operating staff in cooperating during implementation of the project and thanks to Neutron Physics Group, specially H.Ghods and H. Zandi, for assisting radiography and foil activation analysis measurements.

Table 1

Absorbed Dose of Gamma-rays, mrad/hr

Distance from surface of phantom, cm

 

Central axis

 

Left side

 

Right side

            0

   174  ± 26

  160  ± 24

  150  ± 22.5

            2

   260  ± 39

  235  ± 35

  220  ± 33

            4

   220  ± 33

  210  ± 31

  205  ± 31

            6

   165  ± 25

  160  ± 24

  150  ± 22.5

            8

   130  ± 20

  125  ± 19

  120  ± 18

  

Table 2

Absorbed dose due to different reactions along central axis (RBE-mrad/hr)

Distance in   Phantom,
cm

Absorbed Dose of Healthy Tissues

Absorbed Dose of Tumor

Absorbed Dose due to 40ppm 10B

Gamma-ray Absorbed Dose

Epithermal Neutron     Absorbed Dose

Thermal
 Neutron
Absorbed Dose

0.0

235 ± 19.9

 304 ± 22.6

   92 ± 11.2

87 ± 13

123 ± 14.8

2.3 ±  0.27

2.0

285 ± 23.3

 561 ± 48.5

368 ± 44

 130 ± 19.5

54 ± 6.4

9 ± 1.1

3.5

221 ± 19.4

 440 ± 39.2

  292 ± 35.2

   114 ± 17

27 ± 3.2

  7 ± 0.84

5.0

171 ± 15.6

 334 ± 29.6

216 ± 26

97 ± 14

   15.4 ± 1.8

5 ± 0.6

6.5

 126.4 ± 13.2

 225 ± 20.4

132 ± 16

  80 ± 12.5

9.9 ± 1.6

3.5 ± 0.42

8.0

      96 ± 10.47

163.5 ± 14

    90 ± 10.4

65 ± 10

6.2 ± 1.7

2.3 ± 0.27

 

REFERENCES

1.    Soloway A.H., Alarn F., Barth R.F., Anisuzzaman A.K.M. and Bapat B.V. Tumor Targeting Agents for Neutron Capture Therapy, Basic Life Science, Vol. 54, 1990.

2.    Caswell R.S and Coyne J.J., KERMA factor for neutron energies below 30Mev, Radiation Research Vol. 83, 1980.

3.    Marashi M.K., Investigation and production on an epithermal neutron\beam for boron neutron capture therapy, Ph.D thesis, Physics Faculty,    Amir Kabir University, 1993, (supervisor
Ali Pazirandeh)

4.    Matzke Manfred,  Unfolding of  Pulse Height  Spectra, The HEPRO  Program System,
PTB- N-19 Braunschweig, 1994.

5.    Khanbabaii Sh., Unfolding of Neutron Field Spectrum using HEPRO CODES, Physics Department, University of Tehran 1999, (supervisor Ali Pazirandeh).

6.    Beckurts K.H., Wirtz K., Neutron Physics, Springer-Verlag, 1964.

7.    Marashi M.K., Maiorino, J.R., Mendonca A.G., Santos A., IRAN-LIB (Improved Range of ANISN/PC Library), A P3 Coupled Neutron Gamma Cross Section Library in ISOTXS foramat to be used by ANISN/PC, CCC-0514/02), Ann. Nucl. Energy, Vol. 18, No. 10, 597-602, 1991.

8.    Rogus R.D.  , Harling O.K. Yanch J.C.,(1994) Mixed field dosimetry of epithermal neutron for boron capture therapy at the MTR-II research reactor, Med. Phys. 21(10), October 1994.

9.    Raaijimakers C.P.J. and Konijnenberg M.W., (1995) Determination of dose components in phantoms irradiated with an epithermal neutron beam for boron neutron capture therapy, Med. Phys. 22(3), March 1995.

 

Further study:

Wheeler F.J. and Nigg W., Three dimensional radiation distribution analysis for boron capture therapy, Nuc.Sci and Eng., 110(1992) 16-31, 1992

Brugger R.M. , Shin J-L. A. and H.B. Liu,  An epithermal  neutron beam for neutron capture therapy at the Missouri University Research Reactor, Nuclear Technology, 98, June 1992.

Barthy R.F., Soloway A.H., and Fairchild R.G., Boron neutron capture therapy for cancerCancer 7 0,12(1992) December 1992.

Khan Faizm, The Physics of Radiotherapy;
Willaim and Wilkins, Baltimore, Maryland, USA, 1994.


Previous Home Page Next