Toward prompt gamma spectrometry for monitoring boron distributions during extra corporal treatment of liver metastases by boron neutron capture therapy: a Monte Carlo simulation study.

A Monte Carlo calculation was carried out for boron neutron capture therapy (BNCT) of extra corporal liver phantom. The present paper describes the basis for a subsequent clinical application of the prompt gamma spectroscopy set-up aimed at in vivo monitoring of boron distribution. MCNP code was used first to validate the homogeneity in thermal neutron field in the liver phantom and simulate the gamma ray detection system (collimator and detector) in the treatment room. The gamma ray of 478 keV emitted by boron in small specific region can be detected and a mathematical formalism was used for the tomography image reconstruction.

Validating a MCNPX model of Mg(Ar) and TE(TE) ionisation chambers exposed to 60CO gamma rays.

For an accurate determination of the absorbed doses in complex radiation fields (e.g. mixed neutron-gamma fields), a better interpretation of the response of ionisation chambers is required. This study investigates a model of the ionisation chambers using a different approach, analysing the collected charge per minute as a response of the detector instead of the dose. The MCNPX Monte Carlo code is used. In this paper, the model is validated using a well-known irradiation field only: a (60)Co source. The detailed MCNPX models of a Mg(Ar) and TE(TE) ionisation chamber is investigated comparing the measured charge per minute obtained free-in-air and in a water phantom with the simulated results. The difference between the calculations and the measurements for the TE(TE) chamber is within +/-2% whereas for the Mg(Ar) chamber is around +7%. The systematic discrepancy in the case of Mg(Ar) chamber is expected to be caused by an overestimation of the sensitive volume.

An international dosimetry exchange for BNCT part II: computational dosimetry normalizations.

The meaningful sharing and combining of clinical results from different centers in the world performing boron neutron capture therapy (BNCT) requires improved precision in dose specification between programs. To this end absorbed dose normalizations were performed for the European clinical centers at the Joint Research Centre of the European Commission, Petten (The Netherlands), Nuclear Research Institute, Rez (Czech Republic), VTT, Espoo (Finland), and Studsvik, Nyköping (Sweden). Each European group prepared a treatment plan calculation that was bench-marked against Massachusetts Institute of Technology (MIT) dosimetry performed in a large, water-filled phantom to uniformly evaluate dose specifications with an estimated precision of +/-2%-3%. These normalizations were compared with those derived from an earlier exchange between Brookhaven National Laboratory (BNL) and MIT in the USA. Neglecting the uncertainties related to biological weighting factors, large variations between calculated and measured dose are apparent that depend upon the 10B uptake in tissue. Assuming a boron concentration of 15 microg g(-1) in normal tissue, differences in the evaluated maximum dose to brain for the same nominal specification of 10 Gy(w) at the different facilities range between 7.6 and 13.2 Gy(w) in the trials using boronophenylalanine (BPA) as the boron delivery compound and between 8.9 and 11.1 Gy(w) in the two boron sulfhydryl (BSH) studies. Most notably, the value for the same specified dose of 10 Gy(w) determined at the different participating centers using BPA is significantly higher than at BNL by 32% (MIT), 43% (VTT), 49% (JRC), and 74% (Studsvik). Conversion of dose specification is now possible between all active participants and should be incorporated into future multi-center patient analyses.

Gel dosimetry in the BNCT facility for extra-corporeal treatment of liver cancer at the HFR Petten.

A thorough evaluation of the dose inside a specially designed and built facility for extra-corporeal treatment of liver cancer by boron neutron capture therapy (BNCT) at the High Flux Reactor (HFR) Petten (The Netherlands) is the necessary step before animal studies can start. The absorbed doses are measured by means of gel dosemeters, which help to validate the Monte Carlo simulations of the spheroidal liver holder that will contain the human liver for irradiation with an epithermal neutron beam. These dosemeters allow imaging of the dose due to gammas and to the charged particles produced by the (10)B reaction. The thermal neutron flux is extrapolated from the boron dose images and compared to that obtained by the calculations. As an additional reference, Au, Cu and Mn foil measurements are performed. All results appear consistent with the calculations and confirm that the BNCT liver facility is able to provide an almost homogeneous thermal neutron distribution in the liver, which is a requirement for a successful treatment of liver metastases.

Application of adjoint Monte Carlo to accelerate simulations of mono-directional beams in treatment planning for boron neutron capture therapy.

This paper deals with the application of the adjoint transport theory in order to optimize Monte Carlo based radiotherapy treatment planning. The technique is applied to Boron Neutron Capture Therapy where most often mixed beams of neutrons and gammas are involved. In normal forward Monte Carlo simulations the particles start at a source and lose energy as they travel towards the region of interest, i.e., the designated point of detection. Conversely, with adjoint Monte Carlo simulations, the so-called adjoint particles start at the region of interest and gain energy as they travel towards the source where they are detected. In this respect, the particles travel backwards and the real source and real detector become the adjoint detector and adjoint source, respectively. At the adjoint detector, an adjoint function is obtained with which numerically the same result, e.g., dose or flux in the tumor, can be derived as with forward Monte Carlo. In many cases, the adjoint method is more efficient and by that is much quicker when, for example, the response in the tumor or organ at risk for many locations and orientations of the treatment beam around the patient is required. However, a problem occurs when the treatment beam is mono-directional as the probability of detecting adjoint Monte Carlo particles traversing the beam exit (detector plane in adjoint mode) in the negative direction of the incident beam is zero. This problem is addressed here and solved first with the use of next event estimators and second with the application of a Legendre expansion technique of the angular adjoint function. In the first approach, adjoint particles are tracked deterministically through a tube to a (adjoint) point detector far away from the geometric model. The adjoint particles will traverse the disk shaped entrance of this tube (the beam exit in the actual geometry) perpendicularly. This method is slow whenever many events are involved that are not contributing to the point detector, e.g., neutrons in a scattering medium. In the second approach, adjoint particles that traverse an adjoint shaped detector plane are used to estimate the Legendre coefficients for expansion of the angular adjoint function. This provides an estimate of the adjoint function for the direction normal to the detector plane. In a realistic head model, as described in this paper, which is surrounded by 1020 mono-directional neutron/gamma beams and from which the best ones are to be selected, the example calculates the neutron and gamma fluxes in ten tumors and ten organs at risk. For small diameter beams (5 cm), and with comparable relative errors, forward Monte Carlo is seen to be 1.5 times faster than the adjoint Monte Carlo techniques. For larger diameter neutron beams (10 and 15 cm), the Legendre technique is found to be 6 and 20 times faster, respectively. In the case of gammas alone, for the 10 and 15 cm diam beams, both adjoint Monte Carlo Legendre and point detector techniques are respectively 2 and 3 times faster than forward Monte Carlo.

Design of a rotating facility for extracorporal treatment of an explanted liver with disseminated metastases by boron neutron capture therapy with an epithermal neutron beam.

In 2001, at the TRIGA reactor of the University of Pavia (Italy), a patient suffering from diffuse liver metastases from an adenocarcinoma of the sigmoid was successfully treated by boron neutron capture therapy (BNCT). The procedure involved boron infusion prior to hepatectomy, irradiation of the explanted liver at the thermal column of the reactor, and subsequent reimplantation. A complete response was observed. This encouraging outcome stimulated the Essen/Petten BNCT group to investigate whether such an extracorporal irradiation could be performed at the BNCT irradiation facility at the HFR Petten (The Netherlands), which has very different irradiation characteristics than the Pavia facility. A computational study has been carried out. A rotating PMMA container with a liver, surrounded by PMMA and graphite, is simulated using the Monte Carlo code MCNP. Due to the rotation and neutron moderation of the PMMA container, the initial epithermal neutron beam provides a nearly homogeneous thermal neutron field in the liver. The main conditions for treatment as reported from the Pavia experiment, i.e. a thermal neutron fluence of 4 x 10(12) +/- 20% cm(-2), can be closely met at the HFR in an acceptable time, which, depending on the defined conditions, is between 140 and 180 min.

A parameter study to determine the optimal source neutron energy in boron neutron capture therapy of brain tumours.

The values of the parameters used in boron neutron capture therapy (BNCT) to calculate a given dose to human tissue vary with patients due to different physical, biological and/or medical circumstances. Parameters include the tissue dimensions, the 10B concentration and the relative biological effectiveness (RBE) factors for the different dose components associated with BNCT. Because there is still no worldwide agreement on RBE values, more often than not, average values for these parameters are used. It turns out that the RBE-problem can be circumvented by taking into account all imaginable parameter values. Approaching this quest from another angle: the outcome will also provide the parameters (and values) which influence the optimal source neutron energy. For brain tumours it turns out that the 10B concentration, the RBE factors for 10B as well as fast neutrons, together with the dose limit set for healthy tissue, affect the optimal BNCT source neutron energy. By using source neutrons of a few keV together with neutrons of a few eV, it ensures that, under all imaginable circumstances, a maximum of alpha (and lithium) particles can be delivered in the tumour.