Monte Carlo simulation of exposure factor

In this paper, we simulate the exposure factor by a simple model of a free-air ionization chamber with the Monte Carlo programme Geant4. Special emphasis is placed on the discussion of the exposure factor related to parameters of the chamber model. The reason for the variation in exposure factor with incident ray energy is also analysed in terms of reaction cross section for different types of reactions. The obtained results indicate that our simulation is accurate in the calculation of the exposure factor and can serve as a reference in designing air ionization chambers.     

Variational and diffusion Monte Carlo simulations of a hydrogen molecular ion in a spherical box

The variational and diffusion Monte Carlo approaches are used to study the ground-state properties of a hydrogen molecular ion in a spheroidal box. In this work, we successfully treat the zero-point motion of protons in the same formalism with as of electrons and avoid the Born–Oppenheimer approximation in density function theory. The study shows that the total energy increases with the decrease in volume, and that the distance between protons decreases as the pressure increases.Considering the motion of protons, the kinetic energy of the electron is higher than that of the fixed model under the same conditions and increases by 5%. The kinetic energy of the proton is found to be small under high pressure, which is only a fraction of the kinetic energy of the electron.     

Study on gamma response function of E J301 organic liquid scintillator with GEANT4 and FLUKA

The gamma response function is required for energy calibration of EJ301 （5 cm in diameter and 20 cm in height） organic liquid scintillator detector by means of gamma sources. The GEANT4 and FLUKA Monte Carlo simulation packages were used to simulate the response function of the detector for standard 22Na, 60Co, 137Cs gamma sources. The simulated results showed a good agreement with experimental data by incorporating the energy resolution function to simulation codes. The energy resolution and the position of the maximum Compton electron energy were obtained by comparing measured light output distribution with simulated one. The energy resolution of the detector varied from 21.2% to 12.4% for electrons in the energy region from 0.341 MeV to 1.12 MeV. The accurate position of the maximum Compton electron energy was determined at the position 81% of maximum height of Compton edges distribution. In addition, the relation of the electron energy calibration and the effective neutron detection thresholds were described in detail. The present results indicated that both packages were suited for studying the gamma response function of EJ301 detector.     

The penetration, diffusion and energy deposition of high—energy photon＊

This paper presents a new theory for calculating the transport of high-energy photons and their secondary charged particles. We call this new algorithm characteristic line method, which is completely analytic. Using this new method we cannot only accurately calculate the transport behaviour of energetic photons, but also precisely describes the transport behaviour and energy deposition of secondary electrons, photoelectrons, Compton recoil electrons and positron－electron pairs. Its calculation efficiency is much higher than that of the Monte Carlo method. The theory can be directly applied to layered media situation and obtain a pencil-beam-modelled solution. Therefore, it may be applied to clinical applications for radiation therapy.     

A study of the calibration of wire positions for the BESⅢ drift chamber using physics data without magnetic field

A method to determine precisely three-dimensional wire positions in the BESⅢ drift chamber using physics events is introduced. In part the proposed technique takes advantage of the possibility that the magnetic field can be turned off, thereby providing a huge sample of straight tracks for which effects due to multiple scattering, energy loss and non-uniformity of the magnetic field are minimized during the calibration process. A toy Monte Carlo study is performed to demonstrate the viability of the method. As a result of the calibration process, the rms of the distribution of wire position deviation is reduced from 35μm to 10μm.     

Phase Transition of FeS in Terms of In-Situ Resistance Measurement

Electrical properties of stoichiometric iron sulfide （FeS） are investigated under high pressure with a designed diamond anvil cell. The process of phase transition is reflected by changing the electrical conductivity under high pressure, and the conductivity of FeS with the NiAs structure is found to be much smaller than other phases. Two new phase transitions without structural change are observed at 34.7 GPa and 61.3 GPa. The temperature dependence of the conductivity is found to be similar to that of a semiconductor when the pressure is higher than 35 GPa     

Design studies on the 4π γ-ray calorimeter for the ETF experiment at HIRFL-CSR

A high detection efficiency calorimeter which is used to detect γ-rays with energies from 1 MeV up to 10 MeV as well as light charged particles has been proposed. Design of the geometry, results of the crystal tests and Monte Carlo simulations are presented in this paper. The simulation results confirm that the calorimeter can obtain high detection efficiency and good energy resolution with the current designed geometry. And the calorimeter is competent for the future External Target Facility (ETF) experiments.     

Study of the background in the experiments of inner-shell ionization of atoms by positron impact

One of main difficulties in the experiments of inner-shell ionization of atoms by positron impact near threshold energy region is the relatively high low-energy background, which is caused by the deposited part of energy in semi-conductor X-ray detectors of 0.511 MeV γ rays that are produced by positron annihilations in targets and target chamber. In this paper, by using the Monte Carlo method, we simulated the backgrounds for the X-ray detectors with the sensitive layer thickness of 0.3 mm and 3 mm in the case of 0.511 MeV γ rays impacting vertically on a Ti plate of 0.2 mm in thickness, and compared the simulation results with the experimental observations of the other research group and our own. Moreover, we also simulated the backgrounds for a simplified experimental setup in the case of 20 keV positrons impacting vertically on a thick Ti target and observed that the backgrounds for the X-ray detectors with the sensitive layer thickness of 0.3 mm and 3 mm, are very similar.     

Monte Carlo simulation of prompt γ-ray spectra from depleted uranium under D-T neutron irradiation and electron recoil spectra in a liquid scintillator detector

To overcome the problem of inefficient computing time and unreliable results in MCNP5 calculation, a two-step method is adopted to calculate the energy deposition of prompt γ-rays in detectors for depleted uranium spherical shells under D-T neutron irradiation. In the first step, the γ-ray spectrum for energy below 7 MeV is calculated by MCNP5 code; secondly, the electron recoil spectrum in a BC501 A liquid scintillator detector is simulated based on EGSnrc Monte Carlo Code with the γ-ray spectrum from the first step as input. The comparison of calculated results with experimental ones shows that the simulations agree well with experiment in the energy region 0.4–3 MeV for the prompt γ-ray spectrum and below 4 MeVee for the electron recoil spectrum. The reliability of the two-step method in this work is validated.     

Effects due to a Pu-C source on a HPGe detector and the corresponding neutron shielding

A gamma spectrum of a Pu-C source is measured using a p-type HPGe detector, whose three peaks (full energy, single-escape and double-escape peak) can be used as a calibration source for the beam energy measurement system of BEPCII. The effect of fast neutron damage on the energy resolution of the HPGe detector is studied, which indicates that the energy resolution begins to deteriorate when the detector is subject to 2×10⁷ n/cm² fastneutrons. The neutron damage mechanism and detector repair methods are reviewed. The Monte Carlo simulation technique is utilized to study the shielding of the HPGe detector from the fast neutron radiation damage, which is of great significance for the future commissioning of the beam energy measurement system.     

Study of the energy response of high pressure ionization chamber for high energy gamma-ray

The energy response calibration of the commonly used high pressure ionization chamber is very difficult to obtain when the gamma-ray energy is more than 3 MeV.In order to get the calibration of the higher part of the high pressure ionization chamber,we use the Fluka Monte Carlo program to perfclrm the energy response in both the spherical and the cylindrical high pressure ionization chamber which are full of argon gas.The results compared with prior study when the gamma-ray energy is less than 1.25 MeV.Our result of Monte Carlo calculation shows agreement with those obtained by measurement within the uncertainty of the respective methods.The calculation of this study is significant for the high pressure ionization chamber to measure the high energy gamma-ray.     

能量响应平坦的康普顿探测器设计

以开展能量响应平坦的真空康普顿探测器为目的,采用厚薄材料叠加补偿设计思想,使用蒙特卡罗方法,对真空康普顿探测器的灵敏度能量响应进行了优化设计。优化的探测器系统采用0.01,1.00 mm厚Au叠加的发射极,3 mm厚Fe前窗以及3 mm挡铅。在f40 mm准直下,实现了在0.4~7.0 MeV 能区内,探测器本征灵敏度极值变化小于10.7%。该能量响应平坦性优于当前各种常见的探测器。     

利用252 Cf快裂变室测量BC501A液闪探测器的相对探测效率和响应函数

本文利用²⁵²Cf快裂变室和多参数数据采集系统,逐事例的同时记录了自发裂变中子和瞬发伽马的飞行时间(TOF),脉冲形状甄别(PSD)和反冲能量(RE,裂变中子是通过测量反冲质子;瞬发伽马是通过测量康普顿反冲电子)三维信息.详细介绍了通过离线数据分析完全扣除三维信息中的伽马事例贡献,以获得Φ50.8 mm×50.8 mm的BC501A液闪探测器的相对探测效率和响应函数的方法.在不通过探测器响应函数进行数据转换的条件下,利用中子的能量直接确定了中子的有效测量阈值.得到的中子 关键词： 252Cf快裂变室')" href="#">²⁵²Cf快裂变室 BC501A液闪探测器 相对探测效率 响应函数     

Investigation of NaI(Tl)-Crystal Having Anomalous Inhomogeneity of Light Yield Spatial Distribution

An experiment with a layer-by-layer irradiation of the large volume NaI(Tl)-crystal has been performed. The inhomogeneity of the light yield spatial distribution along the crystal’s axis was revealed. A new method of accounting for the revealed inhomogeneity in the processing of the gamma-ray spectra registered using a crystal with nonuniform response is proposed. This method is quite universal: the conventional stages of the detector investigation (the calibration with standard gamma-ray sources and the inhomogeneity estimation using a layer-by-layer irradiation) are complemented by another one stage with computing of the absorbed energy distribution using Monte Carlo simulation. Finally, based on the results of these stages, the combined broadening function is formed for calculation of the response function.     

Response of CaSO4:Dy phosphor based TLD badge system to high energy electron beams from medical linear accelerator and estimation of whole body dose and skin dose

Response of thermoluminescent dosimeter (TLD) badge based on CaSO₄:Dy phosphor to high energy electron beams from medical linear accelerators was simulated using FLUKA Monte Carlo code and experimentally verified. This study was carried out in order to determine the response of TL discs under different filter regions of the TLD badge and their ratios under different irradiation conditions in the energy range 6–20 MeV of high energy electron beams. It was found that the response of TL disc under metal filters of the TLD badge (D1) is always higher than the response of the TL discs under Polystyrene filter (D2) and open window (D3) and its response was found to decrease with increase in the energy of the electron beam. Estimation of whole body dose and skin dose including its error was carried out as an extension of the present study in case of accidental exposure of radiation worker to high energy electron beam.     

非全裸碳离子与氦原子碰撞过程中的多重电离研究

采用飞行时间符合技术测得１．５—７．５ＭｅＶ非全裸碳离子轰击氦原子，引起氦原子双重电离与单重电离的相对截面比。对碰撞过程进行多体经典轨道蒙特卡罗模拟，得到相应的Ｏｌｓｏｎ—Ｓｃｈｌａｃｈｔｅｒ曲线，结合本实验结果，给出非全裸碳离子引起氦原子双重电离和单重电离的绝对截面，并进行了讨论。     

Comparison of thermoluminescence energy response of optical fiber subject to photon irradiation between Monte Carlo simulation and experiments

This study is focused on energy response of Ge-doped silicon dioxide optical fiber subjected to photon irradiation. The TL responses for photon energies, ranging from 20 keV to 20MeV, were investigated as energy absorbed in the TL materials of Ge-doped silicon dioxide optical fiber. The simulation was performed using Monte Carlo N-particle transport code version 5 (MCNP5). The input parameters included in this study are source information, geometry specification, material information and tallies F6. Comparisons of energy response were done between simulation and previous results of experiments. A flat response can only be seen in the energy range of 200 keV to 10 MeV.     

Study of a sealed high gas pressure THGEM detector and response of alpha particle spectra

A sealed high gas pressure detector working in pure argon is assembled.It consists of a 5 cm x 5 cm PCB THGEM(THick Gaseous Electron Multiplier).The detector structure and experimental setup are described.The performance under high pressure(2 atm) is examined,selecting optimal voltages for the ionization region and induction region.The dependence of the shape of alpha particle spectra measured with relative gas gain on gas pressure(1.3-2.0 atm) has been studied.Eight data sets of relative gas gain versus working voltage of THGEM,expressed by weighting field E/P,are normalized,consistent with theory.The results show that the air tightness of the chamber is good,measured by a sensitive barometer and checked with gas gain.The experimental results are compared with Monte Carlo simulation of energy deposition without gas gain involved.     

井型高纯锗探测器体源样品无源效率刻度

采用蒙特卡罗程序MCNP建立物理模型，对井型高纯锗探测器的效率进行虚拟刻度。模拟计算在密度为0.4，1.2 gcm-3的6种不同成分环境样品中，探测器对射线的探测效率，当能量高于0.10 MeV时，体源样品的探测效率主要与样品密度、射线能量相关。以土壤、水和植物油样品为代表，结合所选取的函数模型，确定了密度范围0.1~1.6 gcm-3的固体源、密度1.0 gcm-3的水溶液和密度0.92 gcm-3的油溶液体源的效率函数及参数。实验采用标准源，对探测效率模拟结果进行了验证，探测效率的模拟值与实验值符合较好，二者误差均在3%以内。说明MCNP程序可以较为准确地模拟计算井型高纯锗探测器对射线的探测效率，验证了该无源效率刻度方法的准确性和可行性。