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.     

Study of Borromean halo nuclei by the neutron wall with simulation

The model of three-body Borromean halo nuclei breakup was described by using standard phase space distributions and the Monte Carlo simulation method was established to resolve the detection problem of two neutrons produced from breakup reaction on the neutron wall detector. For 6He case, overall resolution ~rEk for the Oaussiaal part of the detector response and the detection efficiency including solid angle acceptance with regard to the excitation energy Ek are obtained by the simulation of two neutrons from 6He breakup into the neutron wall. The effects of the algorithm on the angular and energy correlations of the fragments are briefly discussed.     

Simulation study using GEANT4 Monte Carlo code for a Gd-coated resistive plate chamber as a thermal neutron detector

The Resistive Plate Chamber (RPC) has been developed in many application areas ever since its introduction, from high energy physics experiments to positron emission tomography. Such detectors can be coated with a Gd layer that enables them to detect thermal neutrons. Consequently these RPCs can be utilized for industrial and medical purposes. Here, we present the configuration of a resistive plate chamber which is utilized to detect thermal neutrons by employing GEANT4 Monte Carlo code. The response of the RPC was evaluated as a function of neutron energy in the GEANT4 Monte Carlo code. The simulation results are taken for incident neutron energy in the energy range from 25 meV to 100 meV. The detection efficiency was found to be between 10% and 20%, depending on the detector configuration, for incident thermal neutrons of 25 meV energy.     

Background neutron in the endcap and barrel regions of resistive plate chamber for compact muon solenoid/large hadron collider using GEANT4

In this study the performance of double gap RPC has been tested by GEANT4 Monte Carlo simulation code. The detector response calculations taken as a function of the neutron energy in the range of 0.01 eV-1 GeV have been simulated through RPC set-up. In order to evaluate the response of detector in the LHC background environment, the neutron spectrum expected in the CMS muon endcap and barrel region were taken into account. A hit rate of about 165.5 Hz cm⁻², 34 Hz cm⁻², 33.6 Hz cm⁻², and 27.0 Hz cm⁻² due to an isotropic neutron source is calculated using GEANT4 standard electromagnetic package for a 20 × 20 cm² RPC in the ME1, ME2, ME3 and ME4, respectively. While for the same neutron source and using GEANT4 package a hit rate of about 0.42 Hz cm⁻², 0.7182 Hz cm⁻² was measured for the MB1 and MB4 stations respectively. Similar characteristics of hit rates have been observed for GEANT4 low electromagnetic package.      

Monte Carlo study for γ+N→π+N at a new compound target

An inbuilt compound target composed of carbon and tungsten is designed,and optimized by realistic GEANT4 Monte Carlo simulation.Also,we do a Monte Carlo study for single-pion photoproduction at the target.The results are presented from the simulation of pion yield,angular distribution and spectrum (at θ1ab,θcm=41°).These efforts are important to the coming measurement of the differential cross section for γ+N→π+N.     

低能中子探测的GEANT4模拟研究

研究了GEANT4蒙特卡罗模拟程序在低能中子核反应中应用的可行性, 新版的GEANT(4.9.4版本)升级了老版本在模拟6Li(n, α)3H核反应道的不足, 并且加入了最新截面数据库。 对多个核反应道进行了研究, 模拟结果显示核反应截面、 次级粒子能量和反应的分支比等是正确的; 比较了掺6Li材料热中子探测效率的实验结果与MCNP和GEANT4模拟结果, 它们基本一致。 研究表明, 利用GEANT4研究6Li(n, α)3H是可行的, 这为进一步研究闪烁体中子探测器的位置分辨率模拟提供了有力支持。 The feasible study of Monte Carlo simulation of low energy nuclear reaction was performed based on GEANT4. The reaction channel 6Li(n, α)3H, simulated insufficiently in the old version Geant4, can be given correct results in the new version (version 4.9.4). New cross section library was added into the program. The study of several nuclear reaction channels shows that the cross sections, the secondary particle energies and the branching ratios are consistent with the experimental values. For doped 6Li materials, the comparison of thermal neutrons detection efficiencies of experimental results, MCNP and GEANT4 simulations, shows that the simulations are reasonable. From the studies above, it is feasible to simulate the 6Li(n, α)3H reaction with GEANT4. It will provide helpful information for the further study of the position resolution of scintillation neutron detector.     

Influence of thyroid volume reduction on absorbed dose in 131I therapy studied by using Geant4 Monte Carlo simulation

A simulation study has been performed to quantify the effect of volume reduction on the thyroid absorbed dose per decay and to investigate the variation of energy deposition per decay due to β- and γ-activity of 131I with volume/mass of thyroid, for water, ICRP- and ICRU-soft tissue taken as thyroid material. A Monte Carlo model of the thyroid, in the Geant4 radiation transport simulation toolkit was constructed to compute the β- and γ-absorbed dose in the simulated thyroid phantom for various values of its volume. The effect of the size and shape of the thyroid on energy deposition per decay has also been studied by using spherical, ellipsoidal and cylindrical models for the thyroid and varying its volume in 1-25 cm3 range. The relative differences of Geant4 results for different models with each other and MCNP results lie well below 1.870%. The maximum relative difference among the Geant4 estimated results for water with ICRP and ICRU soft tissues is not more than 0.225%. S-values for ellipsoidal, spherical and cylindrical thyroid models were estimated and the relative difference with published results lies within 3.095%. The absorbed fraction values for beta particles show a good agreement with published values within 2.105% deviation. The Geant4 based simulation results of absorbed fractions for gammas again show a good agreement with the corresponding MCNP and EGS4 results （±6.667%） but have 29.032% higher values than that of MIRD calculated values. Consistent with previous studies, the reduction of the thyroid volume is found to have a substantial effect on the absorbed dose. Geant4 simulations confirm dose dependence on the volume/mass of thyroid in agreement with MCNP and EGS4 computed values but are substantially different from MIRD8 data. Therefore, inclusion of size/mass dependence is indicated for 131I radiotherapy of the thyroid.     

A geometric factor calculation method based on the isotropic flux assumption

One of the instruments onboard the China Seismic Electromagnetic Satellite （CSES） is the Low Energy Particle Detector （LEPD）. The primary objective of LEPD is to provide measurements of the fluxes, energy spectra and pitch angles of 100 keV to 10 MeV electrons and protons from 2 to 50 MeV in the Earth＇s magnetosphere. The geometric factor is one of the principle parameters of a detector, which converts the physical quantity-count rate to the particle quantity-flux. In this paper, we calculated the geometric factor of LEPD via computer modeling of an isotropic radiation environment. It was first demonstrated that the radiation intensity related should obey a cosine-law, then a general sampling method of generating this distribution via GPS of GEANT4 was explained. Furthermore, combined with flux normalization, a comparison of the geometric factor calculation of a set of 2-layer detectors with different shapes （cylinder, truncated cone and rectangle） was performed. Results show a generally good agreement between simulation and analytical calculations for the cylinder and truncated cone detectors, and the result of the rectangular one, for which there is no accurate analytical formula, is consistent with the previous simulated results by others. As a practical instance of the 2-layer rectangle detector, the geometric factor of LEPD is 10.336±0.036 m cm2·sr for 10 MeV proton and 8.211±0.032 m cm2·sr for 8 MeV electron.     

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.     

GEANT4 Monte Carlo simulation response of parallel plate avalanche counter for fast neutrons detection

In the present study, GEANT4 based Monte Carlo codes have been employed to evaluate parallel plate avalanche counter (PPAC) response for fast neutron detection. In order to detect fast neutrons, a thin polyethylene layer is coated on the surface of the electrode of the PPAC. Neutrons entering the converter produce protons which enter the counter and are detected. Fast neutrons in the energy range of 4.0 MeV–20.0 MeV have been transported onto the PPAC surface using GEANT4 MC code. The performance of the PPAC counter has been evaluated by means of simulation by employing QGSP_BERT_HP and QGSP_BIC_HP physics lists. The detection efficiencies of polyethylene-coated PPAC are 1.69 × 10^?2 and 1.86 × 10^?2 using converter thickness of 1 mm and 2 mm, respectively. The obtained results are discussed in detail.     

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.     

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.     

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.     

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.     

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.