A simulation study of the e+/e－ sensitivity of a low resistive phosphate glasselectrode in a RPC using GEANT MC

The resistivity of conventional glass is quite high and is unacceptable in a high rate environment. Low resistive glass-electrodes could be a solution for this problem. The present study reports the e^＋/e^- simulation results of an RPC detector made from low resistive phosphate glass electrodes. The detailed geometrical configuration of the content materials which are the essential components of the glass of the RPC detector have been created with the GEANT4 simulation code. Two different types of particle sources, i.e. for e^＋/e^- , have been located on the detectors surface to evaluate the performance of the phosphate glass RPC. Both of the particles have been simulated as a function of their respective energies in the range of 0.1 MeV 1.0 GeV. The present simulation work has shown that the resistive electrode plays an important role for the particle production in the RPC configuration.     

Monte Carlo simulation of the property of a scintillation bar in the multi-neutron correlation spectrometer

To perform a kinematically complete measurement of the dissociation reaction for neutron-rich nuclei, a multi-neutron correlation spectrometer is proposed at Peking University. A Monte Carlo simulation code based on GEANT4 is developed for a single scintillation bar which processes not only the energy deposition but also the light propagation in the scintillator and the light collection and conversion to signal at the end of the bar in a realistic way. The simulating method is described in detail in this paper, and the timing and position resolutions and detector efficiency are studied based on the simulation and compared with the experimental results. A new method of crosstalk rejection has been demonstrated to be important for the design of the whole spectrometer.     

An efficient numerical tool for dose deposition prediction applied to synchrotron medical imaging and radiation therapy

Medical imaging and radiation therapy are widely used synchrotron‐based techniques which have one thing in common: a significant dose delivery to typically biological samples. Among the ways to provide the experimenters with image guidance techniques indicating optimization strategies, Monte Carlo simulation has become the gold standard for accurately predicting radiation dose levels under specific irradiation conditions. A highly important hampering factor of this method is, however, its slow statistical convergence. A track length estimator (TLE) module has been coded and implemented for the first time in the open‐source Monte Carlo code GATE/Geant4. Results obtained with the module and the procedures used to validate them are presented. A database of energy‐absorption coefficients was also generated, which is used by the TLE calculations and is now also included in GATE/Geant4. The validation was carried out by comparing the TLE‐simulated doses with experimental data in a synchrotron radiation computed tomography experiment. The TLE technique shows good agreement versus both experimental measurements and the results of a classical Monte Carlo simulation. Compared with the latter, it is possible to reach a pre‐defined statistical uncertainty in about two to three orders of magnitude less time for complex geometries without loss of accuracy.     

Finite Size Effect in Path Integral Monte Carlo Simulations of ^4He Systems

Path integral Monte Carlo （PIMC） simulations are a powerful computational method to study interacting quantum systems at finite temperatures. In this work, PIMC has been applied to study the finite size effect of the simulated systems of ^4He. We determine the energy as a function of temperature at saturated-vapor-pressure （SVP） conditions in the temperature range of T ∈ [1.0 K,4.0 K], and the equation of state （EOS） in the grmmd state For systems consisted of 32, 64 and 128 ^4He atoms, respectively, We find that the energy at SVP is influenced significantly by the size of the simulated system in the temperature range of T ∈ [2.1 K, 3.0 K] and the larger the system is, the better results are obtained in comparison with the experimental values; while the EOS appeared to be unrelated to it.     

Study of neutron response for two hybrid RPC setups using the GEANT4 MC simulation approach

The present article describes a detailed neutron simulation study in the energy range 10^-10 MeV to 1.0 GeV for two different RPC configurations. The simulation studies were taken by using the GEANT4 MC code. Aluminum was utilized on the GND and readout strips for the （a） Bakelite-based and （b） glass-based RPCs. For the former type of RPC setup the neutron sensitivity for the isotropic source was Sn = 2.702 × 10^-2 at En = 1.0 GeV, while for the latter type of RPC, the neutron sensitivity for the same source was evaluated as Sn = 4.049 × 10^-2 at En = 1.0 GeV. These results were further compared with the previous RPC configuration in which copper was used for ground and pickup pads. Additionally A1 was employed at （GND＋strips） of the phosphate glass RPC setup and compared with the copper-based phosphate glass RPC. Good agreement with sensitivity values was obtained with the current and previous simulation results.     

Monte Carlo simulation of indoor external exposure due to gamma-emitting radionuclides in building materials

The use of building materials containing naturally occurring radionuclides such as^40K,^238U,^232Th and their progeny, could lead to external exposures to the residents of such buildings. In this paper, a set of models are constructed to calculate the specific effective dose rates（the effective dose rate per Bq/kg of ^40K, the ^238U series,and the ^232Th series） imposed on residents by building materials with the MCNPX code. The effect of chemical composition, position concerned in the room and thickness as well as density of material is analyzed. In order to facilitate more precise assessment of indoor external dose due to gamma-emitting radionuclides in building materials,three regressive expressions are proposed and validated by measured data to calculate specific effective rates for40 K,the238U series and the232 Th series, respectively.     

Ep≤63 GeV 质子入射铅、铋引起的中子产生双微分截面的模拟(英文)

加速器驱动次临界系统(ADS) 液态Pb-Bi 散裂靶的设计中,需要可靠的理论计算工具精确地预言几个GeV 能量范围的质子引起的散裂反应产生的各种粒子和核素。利用蒙特卡罗模拟软件包Geant4 计算研究了800 MeV至3 GeV 质子入射铅、铋材料引起的中子产生双微分截面。比较了Geant4 不同物理模型得到的模拟结果与现有的实验数据。其中,Geant4 的QGSP BERT和QGSP INCL ABLA 物理模型模拟结果很好地再现了实验数据。本工作证实了Geant4 蒙特卡罗模拟软件包适合用于能量高达3 GeV 的质子入射铅、铋引起的中子产生双微分截面的模拟计算。A detailed design of the liquid Pb-Bi spallation target of the Accelerator Driven Systems (ADS) requires powerful and reliable computational tools that can accurately predict particles and nuclides production by the proton induced spallation reactions in the energy range of a few GeV. In this paper, the neutron production double-differential cross sections for Pb and Bi target materials at incident proton kinetic energies between 800 MeV and 3 GeV are studied by calculations with Monte Carlo simulation package Geant4. The simulated results of Geant4 with several physics models are compared with available experimental data. The simulated results generated by QGSP BERT and QGSP INCL ABLA physics models of Geant4 well reproduce the available experimental data. The present results validated that Geant4 Monte Carlo simulation package is suitable for simulations of neutron production double-differential cross sections of proton induced reaction on Pb and Bi targets in the incident energy range up to 3 GeV.     

News on PHOTOS Monte Carlo： γ＊→π^＋π^-（γ） and K^±→π^＋π^-e^±υ（γ）

PHOTOS Monte Carlo is widely used for simulating QED effects in decay of intermediate particles and resonances. It can be easily connected to other main process generators. In this paper we consider decaying processes γ＊→π^＋π^-（γ） and K^±→π^＋π^-e^±υ（γ） in the framework of Scalar QED. These two processes are interesting not only for the technical aspect of PHOTOS Monte Carlo, but also for precision measurement of αED（Mz）, g--2, as well as ππ scattering lengths.     

A Study of Surfactant-Induced Growth of Ag on Ag （111）

A new growth model is introduced to describe surfactant-induced growth of Ag on Ag （111） with realistic physical parameters. In this model, the A-S exchange mechanism is considered for the first time. Using the Monte Carlo simulations, the influence of exchange mechanism, surface temperature T, the exchange barrier Eεx, and the coverage of surfactant θM on the growth mode and morphology during multilayer film growth of Ag/Ag （111） are studied in detail Both the referenced value of surfactant coverage and the method to obtain perfect layer-by-layer film in surfactantinduced Ag/Ag （111） system are provided. Our simulation results are consistent with many experimental observations for surfactant-induced growth of Ag on Ag （111）.     

Coverage Dependence of Growth Mode in Heteroepitaxy of Ni on Cu（100） Surface

A realistic kinetic Monte Carlo （KMC） simulation model with physical parameters is developed, which well reproduces the heteroepitaxial growth of multilayered Ni thin film on Cu（100） surfaces at room temperature. The effects of mass transport between interlayers and edge diffusion of atoms along the islands are included in the simulation model, and the surface roughness and the layer distribution versus total coverage are calculated. Speeially, the simulation model reveals the transition of growth mode with coverage and the difference between the Ni heteroepitaxy on Cu（100） and the Ni homoepitaxy on Ni（100）. Through comparison of KMC simulation with the real scanning tunneling microscopy （STM） experiments, the Ehrlich-Schwoebel （ES） barrier Ees is estimated to be 0.18±0.02 eV for Ni/Cu（100） system while 0.28 eV for Ni/Ni（100）. The simulation also shows that the growth mode depends strongly on the thickness of thin film and the surface temperature, and the critical thickness of growth mode transition is dependent on the growth condition such as surface temperature and deposition flux as well.     

Speckle intensity images of target based on Monte Carlo method

Speckle intensity in the detector plane is deduced in the free-space optical system and imaging system based on Van Cittert-Zemike theorem. The speckle intensity images of plane target and conical target are obtained by using the Monte Carlo method and measured experimentally. The results show that when the range extent of target is smaller, the speckle size along the same direction become longer, and the speckle size increase with increasing incident light wavelengths. The speckle size increases and the speckle intensity images of target is closer to the actual object when the aperture scale augments. These findings are useful to access the target information by speckle in laser radar systems.     

Theoretical study of molecular hydrogen and spiltover hydrogen storage on two-dimensional covalent-organic frameworks

Molecular hydrogen and spiltover hydrogen storages on five two-dimensional （2D） covalent-organic frameworks （COFs） （PPy-COF, TP-COF, BTP-COF, COF-18 A, and HHTP-DPB COF） are investigated using the grand canonical Monte Carlo （GCMC） simulations and the density functional theory （DFT）, respectively. The GCMC simulated results show that HHTP-DPB COF has the best performance for hydrogen storage, followed by BTP-COF, TP-COF, COF-18 A, and PPy-COE However, their adsorption amounts at room temperature are all too low to meet the uptake target set by US Department of Energy （US-DOE） and enable practical applications. The effects of pore size, surface area, and isosteric heat of hydrogen on adsorption amount are considered, which indicate that these three factors are all the important factors for determining the H2 adsorption amount. The chemisorptions of spiltover hydrogen atoms on these five COFs represented by the cluster models are investigated using the DFT method. The saturation cluster models are constructed by considering all possible adsorption sites for these cluster models. The average binding energy of a hydrogen atom and the saturation hydrogen storage density are calculated. The large average binding energy indicates that the spillover process may pro- ceed smoothly and reversibly. The saturation hydrogen storage density is much larger than the physisorption uptake of H2 molecules at 298 K and 100 bar （1 bar = 105 Pa）, and is close to or exceeds the 2010 US-DOE target of 6 wt% for hydrogen storage. This suggests that the hydrogen storage capacities of these COFs by spillover may be significantly enhanced. Thus 2D COFs studied in this paper are suitable hydrogen storage media by spillover.     

蒙特卡罗技术在放射诊断剂量快速计算中的应用

采用蒙特卡罗（MC）技术研究放射性核素131I的放疗特性。给出了131I衰变过程的MC抽样表达式，实现了含电子和光子的混合抽样模拟。对水客体中131I衰变的MC模拟结果表明:131I能够有效地治疗甲状腺肿瘤而几乎不影响正常组织。变密度客体的MC模拟结果表明:沉积能量与面密度之间的关系与密度无关，满足卷积核自适应的要求；并且MC模拟能够提供剂量卷积核函数，为快速剂量计算提供数据。     

Large energy-loss straggling of swift heavy ions in ultra-thin active silicon layers

Monte Carlo simulations reveal considerable straggling of energy loss by the same ions with the same energy in fully-depleted silicon-on-insulator (FDSOI) devices with ultra-thin sensitive silicon layers down to 2.5 nm. The absolute straggling of deposited energy decreases with decreasing thickness of the active silicon layer. While the relative straggling increases gradually with decreasing thickness of silicon films and exhibits a sharp rise as the thickness of the silicon film descends below a threshold value of 50 nm, with the dispersion of deposited energy ascending above ±10%. Ion species and energy dependence of the energy-loss straggling are also investigated. For a given beam, the dispersion of deposited energy results in large uncertainty on the actual linear energy transfer (LET) of incident ions, and thus single event effect (SEE) responses, which pose great challenges for traditional error rate prediction methods.     

Strong and electromagnetic J/ψ and ψ（2S） decays into pion and kaon pairs

We discuss interference effects important for the form factors extraction in the vicinity of J/ψ and ψ（2S） resonances in combination with resonance parameters determination. The implementation to the Monte Carlo event generator PHOKHARA of the J/ψ and ψ（2S） contributions to the muon, pion and kaon pairs production associated with a photon at next-to-leading order is also described.     

Energy of Long-Lifetime Configurations in Zero-Temperature Dynamics

Using Monte Carlo method with zero-temperature dynamics, we investigate energy evolution of Ising spin configuration on a square lattice. The energies of some configurations exhibit long duration before those configurations reach the final state -- ground state or frozen stripe state. For ground-state dynamical realization, the duration occurs when the energy per spin is 4/L, where L is the lattice size. For stripe-state dynamical realization, the energy is slightly higher than 2/L when the duration appears in the last evolution stage. In addition, it is found that the average energy per spin in final state is approximately 2/3L.     

Design and simulations for the detector based on DSSSD

The present paper describes the design and simulation results of a position-sensitive charged particle detector based on the Double Sided Silicon Strip Detector （DSSSD）. Also, the characteristics of the DSSSD and its testing result were are discussed. With the application of the DSSSD, the position-sensitive charged particle detector can not only give particle flux and energy spectra information and identify different types of charged particles, but also measure the location and angle of incident particles. As the detector can make multi-parameter measurements of charged particles, it is widely used in space detection and exploration missions, such as charged particle detection related to earthquakes, space environment monitoring and solar activity inspection.     

基于Geant4的反应堆临界及动态计算方法

Geant4是一款基于C++面向对象技术的蒙特卡罗开发程序包,可以模拟各种已知粒子与物质之间的相互作用。然而该程序包没有提供临界源功能,无法直接用于反应堆物理计算。因此,利用Geant4提供的基础物理模型和粒子跟踪控制等功能,用两种不同方法实现了临界源的设置,实现了基于Geant4的反应堆静态计算程序G4-RSM和反应堆动态计算程序G4-RDM。两个程序均可用于反应堆临界计算,与MCNP计算结果相对误差在5%以内。G4-RDM程序除可用于临界计算外,还可用于模拟堆内事故工况下的中子学瞬态变化。     

Absorbed dose estimations of 131I for critical organs using the GEANT4 Monte Carlo simulation code

The aim of this study is to compare the absorbed doses of critical organs of ¹³¹I using the MIRD (Medical Internal Radiation Dose) with the corresponding predictions made by GEANT4 simulations. S-values (mean absorbed dose rate per unit activity) and energy deposition per decay for critical organs of ¹³¹I for various ages, using standard cylindrical phantom comprising water and ICRP soft-tissue material, have also been estimated. In this study the effect of volume reduction of thyroid, during radiation therapy, on the calculation of absorbed dose is also being estimated using GEANT4. Photon specific energy deposition in the other organs of the neck, due to ¹³¹I decay in the thyroid organ, has also been estimated. The maximum relative difference of MIRD with the GEANT4 simulated results is 5.64% for an adult's critical organs of ¹³¹I. Excellent agreement was found between the results of water and ICRP soft tissue using the cylindrical model. S-values are tabulated for critical organs of ¹³¹I, using 1, 5, 10, 15 and 18 years (adults) individuals. S-values for a cylindrical thyroid of different sizes, having 3.07% relative differences of GEANT4 with Siegel & Stabin results. Comparison of the experimentally measured values at 0.5 and 1 m away from neck of the ionization chamber with GEANT4 based Monte Carlo simulations results show good agreement. This study shows that GEANT4 code is an important tool for the internal dosimetry calculations.     

Reverse electric field Monte Carlo simulation for vector radiative transfer in the atmosphere

In this paper, a reverse electric field Monte Carlo （REMC） method is proposed to study the vector radiation transfer in the atmosphere. The REMC is based on tracing the multiply scattered electric field to simulate the vector transmitted radiance. The reflected intensities with different total optical depth values are obtained, which accord well with the results in the previous research. Stokes vector and the degree of polarization are numerically investigated. The simulation result shows that when the solar zenith angle is determined, the zenith angle of detector has two points, of which the degree of polarization does not change with the ground albedo and the optical depth. The two points change regularly with the solar zenith angle. Moreover, our REMC method can be applied to the vector radiative transfer in the atmosphere-ocean system.