Search

TLD-400 (CaF2:Mn) chips were applied for the gamma-ray dose measurement in a PMMA phantom exposed to a BNCT beam because of their very low neutron sensitivity. Since TLD-400 chips possess an adequate amount of Mn activator they have been employed in this work simultaneously for neuron activation measurement. The self-irradiation TL signals owing to the decay of the neutron induced 56Mn activity have been applied for a calibration of the TLD-400 chip in situ, where the activities were measured by an HPGe detector system and the energy deposition per disintegration of 56Mn was calculated by applying a Monte Carlo code. It was accidentally found that the irradiated TLD-400 chips were capable of emitting prominent scintillation lights owing to the induced 56Mn activity, which can easily be recorded by the TLD reader without heating and after a calibration can be used to determine the 56Mn activity.

In this work a modified absolute calibration method of neutron activation analysis (NAA) was proposed. Reaction rates per atom of elements of interest were calculated by using MCNPX Monte Carlo code with a well established disk surface neutron source. Time dependency of the absolute reaction rate per atom produced during the irradiation was achieved by utilizing a well calibrated neutron monitoring system. Three activation foils with given compositions, namely, Au(1%)/Al, Mn(88%)/Ni and Cu(100%) and a SS304 sample were measured to test the adequacy and performance of this proposed modified absolute calibration method.

The neutron beam monitoring system is indispensable to BNCT facility in order to achieve an accurate patient dose delivery. The neutron beam monitoring of a reactor-based BNCT (RB-BNCT) facility can be implemented through the instrumentation and control system of a reactor provided that the reactor power level remains constant during reactor operation. However, since the neutron flux in reactor core is highly correlative to complicated reactor kinetics resulting from such as fuel depletion, poison production, and control blade movement, some extent of variation may occur in the spatial distribution of neutron flux in reactor core. Therefore, a dedicated neutron beam monitoring system is needed to be installed in the vicinity of the beam path close to the beam exit of the RB-BNCT facility, where it can measure the BNCT beam intensity as closely as possible and be free from the influence of the objects present around the beam exit. In this study, in order to demonstrate the importance of a dedicated BNCT neutron beam monitoring system, the signals originating from the two in-core neutron detectors installed at THOR were extracted and compared with the three dedicated neutron beam monitors of the THOR BNCT facility. The correlation of the readings between the in-core neutron detectors and the BNCT neutron beam monitors was established to evaluate the improvable quality of the beam intensity measurement inferred by the in-core neutron detectors. In 29 sampled intervals within 16 days of measurement, the fluctuations in the mean value of the normalized ratios between readings of the three BNCT neutron beam monitors lay within 0.2%. However, the normalized ratios of readings of the two in-core neutron detectors to one of the BNCT neutron beam monitors show great fluctuations of 5.9% and 17.5%, respectively.

The paired ionization chambers (ICs) technique is commonly employed to determine neutron and photon doses in radiology or radiotherapy neutron beams, where neutron dose shows very strong dependence on the accuracy of accompanying high energy photon dose. During the dose derivation, it is an important issue to evaluate the photon and electron response functions of two commercially available ionization chambers, denoted as TE(TE) and Mg(Ar), used in our reactor based epithermal neutron beam. Nowadays, most perturbation corrections for accurate dose determination and many treatment planning systems are based on the Monte Carlo technique. We used general purposed Monte Carlo codes, MCNP5, EGSnrc, FLUKA or GEANT4 for benchmark verifications among them and carefully measured values for a precise estimation of chamber current from absorbed dose rate of cavity gas. Also, energy dependent response functions of two chambers were calculated in a parallel beam with mono-energies from 20 keV to 20 MeV photons and electrons by using the optimal simple spherical and detailed IC models. The measurements were performed in the well-defined (a) four primary M-80, M-100, M120 and M150 X-ray calibration fields, (b) primary 60Co calibration beam, (c) 6 MV and 10 MV photon, (d) 6 MeV and 18 MeV electron LINACs in hospital and (e) BNCT clinical trials neutron beam. For the TE(TE) chamber, all codes were almost identical over the whole photon energy range. In the Mg(Ar) chamber, MCNP5 showed lower response than other codes for photon energy region below 0.1 MeV and presented similar response above 0.2 MeV (agreed within 5% in the simple spherical model). With the increase of electron energy, the response difference between MCNP5 and other codes became larger in both chambers. Compared with the measured currents, MCNP5 had the difference from the measurement data within 5% for the 60Co, 6 MV, 10 MV, 6 MeV and 18 MeV LINACs beams. But for the Mg(Ar) chamber, the derivations reached 7.8–16.5% below 120 kVp X-ray beams. In this study, we were especially interested in BNCT doses where low energy photon contribution is less to ignore, MCNP model is recognized as the most suitable to simulate wide photon–electron and neutron energy distributed responses of the paired ICs. Also, MCNP provides the best prediction of BNCT source adjustment by the detector’s neutron and photon responses.

This study used simulations of galactic cosmic ray in the atmosphere to investigate the neutron background environment in Taiwan, emphasising its altitude dependence and spectrum variation near interfaces. The calculated results were analysed and compared with two measurements. The first measurement was a mobile neutron survey from sea level up to 3275 m in altitude conducted using a car-mounted high-sensitivity neutron detector. The second was a previous measured result focusing on the changes in neutron spectra near air/ground and air/water interfaces. The attenuation length of cosmic-ray neutrons in the lower atmosphere was estimated to be 163 g cm(-2) in Taiwan. Cosmic-ray neutron spectra vary with altitude and especially near interfaces. The determined spectra near the air/ground and air/water interfaces agree well with measurements for neutrons below 10 MeV. However, the high-energy portion of spectra was observed to be much higher than our previous estimation. Because high-energy neutrons contribute substantially to a dose evaluation, revising the annual sea-level effective dose from cosmic-ray neutrons at ground level in Taiwan to 35 μSv, which corresponds to a neutron flux of 5.30 × 10(-3) n cm(-2) s(-1), was suggested.

The (7)Li(p,xn)(7)Be nuclear reaction, based on the low-energy protons, could produce soft neutrons for accelerator-based boron neutron capture therapy (AB-BNCT). Based on the fact that the induced neutron field is relatively divergent, the relationship between the incident angle of proton beam and the neutron beam quality was evaluated in this study. To provide an intense epithermal neutron beam, a beam-shaping assembly (BSA) was designed. And a modified Snyder head phantom was used in the calculations for evaluating the dosimetric performance. From the calculated results, the intensity of epithermal neutrons increased with the increase in proton incident angle. Hence, either the irradiation time or the required proton current can be reduced. When the incident angle of 2.5-MeV proton beam is 120°, the required proton current is ∼13.3 mA for an irradiation time of half an hour.

Proton accelerators in the energy range of approximately 200 MeV have become increasingly popular for cancer treatment in recent years. These proton therapy facilities usually involve bulky concrete or iron in their shielding design or accelerator structure. Simple shielding data, such as source terms or attenuation lengths for various proton energies and materials are useful in designing accelerator shielding. Understanding the appropriateness or uncertainties associated with these data, which are largely generated from Monte Carlo simulations, is critical to the quality of a shielding design. This study demonstrated and investigated the problems of deep-penetration calculations on the estimation of shielding parameters through an extensive comparison between the FLUKA and MCNPX calculations for shielding against a 200-MeV proton beam hitting an iron target. Simulations of double-differential neutron production from proton bombardment were validated by comparison with experimental data. For the concrete shielding, the FLUKA calculated depth–dose distributions were consistent with the MCNPX results, except for some discrepancies in backward directions. However, for the iron shielding, if FLUKA is used inappropriately then overestimation of neutron attenuation can be expected as shown by this work because of the multigroup treatment for low-energy neutrons in FLUKA. Two neutron energy group structures, three degrees of self-shielding correction, and two iron compositions were considered in this study. Significant variation of the resulting attenuation lengths indicated the importance of problem-dependent multigroup cross sections and proper modeling of iron composition in deep-penetration calculations.

High energy photon (over 10 MeV) and neutron beams adopted in radiobiology and radiotherapy always produce mixed neutron/gamma-ray fields. The Mg(Ar) ionization chambers are commonly applied to determine the gamma-ray dose because of its neutron insensitive characteristic. Nowadays, many perturbation corrections for accurate dose estimation and lots of treatment planning systems are based on Monte Carlo technique. The Monte Carlo codes EGSnrc, FLUKA, GEANT4, MCNP5, and MCNPX were used to evaluate energy dependent response functions of the Exradin M2 Mg(Ar) ionization chamber to a parallel photon beam with mono-energies from 20 keV to 20 MeV. For the sake of validation, measurements were carefully performed in well-defined (a) primary M-100 X-ray calibration field, (b) primary 60Co calibration beam, (c) 6-MV, and (d) 10-MV therapeutic beams in hospital. At energy region below 100 keV, MCNP5 and MCNPX both had lower responses than other codes. For energies above 1 MeV, the MCNP ITS-mode greatly resembled other three codes and the differences were within 5%. Comparing to the measured currents, MCNP5 and MCNPX using ITS-mode had perfect agreement with the 60Co, and 10-MV beams. But at X-ray energy region, the derivations reached 17%. This work shows us a better insight into the performance of different Monte Carlo codes in photon–electron transport calculation. Regarding the application of the mixed field dosimetry like BNCT, MCNP with ITS-mode is recognized as the most suitable tool by this work.

This study aims to measure the neutronic characteristics of the 14-cm diameter boron neutron capture therapy beam at the Tsing Hua Open-pool Reactor and to provide essential information for the neutron source adjustment and validation. The measurements were performed both in air and inside a cubic PMMA phantom at the beam, with and without an 18-cm long extended PE collimator. The neutron intensity was determined by neutron activation analysis; the indirect neutron radiography and the cadmium difference method were coupled to provide the two-dimensional neutron flux distributions. According to the measurements, the angular distribution is highly forward for epithermal neutrons but much more divergent for thermal ones. A PE extended collimator will modify the beam characteristics and could benefit the treatment for head and neck tumors, which often needs higher therapeutical boron dose at the shallow region.

This paper presents a way to determine the angular and spatial distributions of the thermal neutron source strength of a boron neutron capture therapy (BNCT) beam. The experiments applied 1) the indirect neutron radiography, 2) the cadmium difference method, and 3) the instrumental neutron activation analysis. The measured data were processed by the spectrum deconvolution technique to resolve into a proper set of angular and spatial distributions. This paper took the epithermal neutron beam of the BNCT facility at the Tsing Hua Open-pool Reactor as an example.

This study is a part of the beam comparison campaign, inter-center dose comparison, between boron neutron capture therapy facilities at the Tsing Hua Open-pool Reactor and the High Flux Reactor. The clinical information exchange can improve the dosimetry uncertainty for medical physics in a mixed field. The method of paired Mg(Ar) and TE(TE) ionization chambers was used to determine the gamma-ray and neutron dose rates. Furthermore, activation foils, including gold, copper, and manganese, were employed to estimate the thermal and epithermal neutron fluxes. Measurements were performed free in air and also in a PMMA phantom. All the chambers were calibrated using a 60Co primary standard source at the Institute of Nuclear Energy Research, Taiwan. Spectrum dependent neutron sensitivity of TE(TE) chamber is one of the important parameters to evaluate dose components. The requested neutron spectra were calculated by the Monte Carlo code MCNP. The measured thermal neutron fluxes, gamma-ray and neutron dose rates of the THOR beam in the phantom were 2.6, 2.2, and 2.1 times of the HFR beam at 2.5-cm depth, respectively. The higher thermal neutron flux and neutron and gamma-ray dose rates are due to the higher epithermal neutron beam intensity of the THOR.

This study aims to determine the activation contamination in the metal-based ionization chambers, i.e. Mg(Ar)- and Al(Ar)-chambers. The contamination was estimated by directly measuring the signal caused by the radioactive nuclides after the irradiation. The INAA gamma ray spectra clearly showed that the ionization chambers were activated. The activation contamination could comprise more than 2% of the total measured current of the Mg(Ar)-chamber, and more than 50% of the Al(Ar)-chamber, both irradiated free-in-air at the THOR epithermal neutron beam. The activation contamination provides concrete evidence of the necessity of the thermal neutron sensitivity correction.

For the nuclear safeguard purpose, this work aims to nondestructively determine the plutonium content in the Taiwan Research Reactor (TRR) spent fuel rods in the storage pool before the stabilization process, which transforms the metal spent fuel rods into oxide powder. A SPent-fuel-Neutron-Counter (SPNC) system was designed and constructed to carry out underwater scan measurements of neutrons emitting from the spent fuel rod, from which the 240 Pu mass in the fuel rod will be determined. The SAS2 H control module of the SCALE 5.1 code package was applied to calculate the 240 Pu-to-Pu mass ratio in the TRR spent fuel rod according to the given power history. This paper presents the methodology and design of our detector system as well as the measurements of four TRR spent fuel rods in the storage pool and the comparison of the measured results with the facility declared values.

Activation foils, stacked in a package, together with a spectrum adjustment code are generally the means to characterize the neutron spectrum of a BNCT epithermal neutron beam. When using foil stacks, the phenomenon of neutron scattering is appreciably enhanced, especially in an epithermal neutron beam, as compared to a thermal neutron field. Contrary to the conventional idea of the self-shielding effect, self-shielding correction factors in most of the foils in the stack are larger than one. This paper demonstrates a way to correct the complex self-shielding effect using groupwise cross-sections modified for self-shielding, as calculated by the Monte Carlo neutron transport code, MCNP. This methodology can solve the perturbation problems caused by the foil stack applied in an epithermal neutron beam. In this study, the adjusted spectrum of the HB11 epithermal neutron beam, dedicated for BNCT research in Petten, The Netherlands, is utilized and two adjusted spectra are presented with and without the self-shielding correction. The difference between the spectra is highly significant and demonstrates that the self-shielding correction cannot be ignored.

In order to provide an improved and reliable neutron source description for treatment planning in boron neutron capture therapy (BNCT), a spectrum adjustment procedure named coarse-scaling adjustment has been developed and applied to the neutron spectrum measurements of both the Tsing Hua Open-pool Reactor (THOR) epithermal neutron beam in Taiwan and the High Flux Reactor (HFR) in The Netherlands, using multiple activation detectors. The coarse-scaling adjustment utilizes a similar idea as the well-known two-foil method, which adjusts the thermal and epithermal neutron fluxes according to the Maxwellian distribution for thermal neutrons and 1/ E distribution over the epithermal neutron energy region. The coarse-scaling adjustment can effectively suppress the number of oscillations appearing in the adjusted spectrum and provide better smoothness. This paper also presents a sophisticated 9-step process utilizing twice the coarse-scaling adjustment which can adjust a given coarse-group spectrum into a fine-group structure, i.e. 640 groups, with satisfactory continuity and excellently matched reaction rates between measurements and calculation. The spectrum adjustment algorithm applied in this study is the same as the well-known SAND-II.

The Cyclotron Center of Shin Kong Memorial Hospital is the most active [ 18 F]FDG producer in Taiwan. Any improvement of the production efficiency of the MINItrace™ cyclotron is valuable. A theoretical activity of [ 18 F]FDG has been calculated from the proton stopping power and the interaction cross-section of the 18 O(p,n) 18 F reaction to estimate the efficiency of the production process. A high-sensitivity moderated 3 He neutron detector was installed to investigate the real-time condition of the target irradiation by monitoring the associated neutron production. Based on the results of this study, a graphical user interface program has been developed for operators to arrange the operation conditions and to analyze the performance of the production process.

This study investigated the characteristics of gas bremsstrahlung and induced photoneutrons from the electron storage ring in National Synchrotron Radiation Research Center (NSRRC). Gas bremsstrahlung is generated in the vacuum chamber of an electron storage ring by interactions of the circulating electrons with the atoms of the residual gas. The emitted photons with sufficiently high energies can further produce neutrons by means of photonuclear interactions. Photoneutrons produced in a lead target struck by gas bremsstrahlung were measured as a function of target thickness by using high-efficiency neutron detectors in a well-defined configuration. A Monte Carlo simulation of the whole process including generation of gas bremsstrahlung, production of photoneutrons and subsequent neutron transport has been performed in an integrated manner and the results were directly compared with the measurements on an absolute basis. The agreement between them is very satisfactory and hence validates the evaluation of the characteristics of gas bremsstrahlung and induced photoneutrons in NSRRC.

This work presents an interactive code, WRNS, based on data base interpolation for preliminary neutron shielding design. A wide-range data base of transmission factors is also generated for nine shielding materials having thicknesses up to 200 cm and for neutron energies up to 400 MeV by using adjoint calculations of a one-dimensional discrete ordinates code. In addition, calculation details of the interpolation code are clearly interpreted along with WRNS code applications presented as well.

Standard Bonner spheres and proposed high-sensitivity Bonner cylinders were calibrated in a neutron calibration room, using a (252)Cf source. The Bonner sphere system consists of 11 polyethylene (PE) spheres of various diameters and 4 extended spheres that comprise embedded metal shells. Similar to the design of Bonner spheres, a set of Bonner cylinders was assembled using a large cylindrical (3)He tube as the central probe, which was wrapped using various thicknesses of PE. A layer of lead was employed inside one of the PE cylinders to increase the detection efficiency of high-energy neutrons. The central neutron probe used in the Bonner cylinders exhibited an efficiency of ∼17.9 times higher than that of the Bonner spheres. However, compared with the Bonner spheres, the Bonner cylinders are not fully symmetric in their geometry, exhibiting angular dependence in their responses to incoming neutrons. Using a series of calculations and measurements, this study presents a systematic comparison between Bonner spheres and cylinders in terms of their response functions, detection efficiencies, angular dependences and spectrum unfolding.

High energy proton beam (>8MeV) is favorable for producing neutrons with high yield. However, the produced neutrons are of high energies. These high energy neutrons can cause severe fast neutron contamination and degrade the BNCT treatment quality if they are not appropriately moderated. Hence, this study aims to briefly discuss the issue, from the viewpoint of fast neutron contamination control, whether high energy proton beam is ideal for AB-BNCT or not. In this study, D2O, PbF4, CaF2, and Fluental(™) were used standalone as moderator materials to slow down 1-, 6-, and 10-MeV parallelly incident neutrons. From the calculated results, we concluded that neutrons produced by high energy proton beam could not be easily moderated by a single moderator to an acceptable contamination level and still with reasonable epithermal neutron beam intensity. Hence, much more complicated and sophisticated designs of beam shaping assembly have to be developed when using high energy proton beams.

In order to monitor the neutron dose in the environment of the Synchrotron Radiation Research Center, the sensitivity of the monitor system to the natural background cosmic-ray neutrons must be established. Large-size BF/sub 3/ counters with and without an 8.5-cm polyethylene enclosure were used to measure the cosmic-ray neutrons in the Nuclear Engineering Building, NTHU. The response functions of the BF/sub 3/ counters were calculated by using ANISN code performed in adjoint mode. By applying the cosmic-ray neutron spectrum and the detector response functions the measured counting rates of the BF/sub 3/ counters were converted to deliver the annual dose of cosmic-ray neutrons. The shielding effect of the concrete slabs and windows of the building was determined by the measured data and compared to that predicted by neutron transport calculations. >

This paper aims to evaluate the performance of the source description of the THOR BNCT beam via different measurement techniques in different phantoms. The measurement included (1) the absolute reaction rate measurement of a set of triple activation foils, (2) the neutron and gamma-ray dose rates measured using the paired ionization chamber method, and (3) the relative reaction rate distributions obtained using the indirect neutron radiography. Three source descriptions, THOR-Y09, surface source file RSSA, and THOR-50C, were tested. The comparison results concluded that THOR-Y09 is a well-tested source description not only for neutron components, but also for gamma-ray component.

A refined dose assessment method has been used now in the THOR BNCT facility, which takes into account more delicate corrections, carefully handled calibration factors, and the spectrum- and kerma-weighted k t value. The refined method solved the previous problem of negative derived neutron dose in phantom at deeper positions. With the improved dose assessment, the calculated and measured gamma-ray dose rates match perfectly in a 15×15×15 cm 3 PMMA phantom.

This study aims to measure the two-dimensional (2D) neutron spatial distribution inside a cubic and a head PMMA phantom for the purpose of further comparison with the treatment planning. The measurements were made by using the indirect neutron radiography (INR), which utilized a thin copper foil and the imaging plate. The developed image provides satisfactory spatial resolution and very low statistical error (< 1%). As to the time cost, the whole procedure normally takes less than 3 h. The result shows that the indirect neutron radiography can be a quick and reliable method to provide a 2D neutron spatial distribution inside a phantom.

Full-scale Monte Carlo simulations of the cyclotron room of the Buddhist Tzu Chi General Hospital were carried out to improve the inadequate maze design. The double differential neutron source from the 1 8 O(p,n) 1 8 F reaction was adopted for the calculation. The weight window variance reduction technique, where the weight window was set by applying the adjoint flux, has been implemented in the MCNP run to facilitate the calculation of the dose rates outside the cyclotron room. Dose rates including neutron and gamma-ray components were calculated for some maze shielding modifications.