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 ⁴⁰K, ²³⁸U, ²³²Th 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 ⁴⁰K, the ²³⁸U series, and the ²³²Th 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 for ⁴⁰K, the ²³⁸U series and the ²³²Th series, respectively.     

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.     

Monte Carlo simulation of magnetic nanostructured thin films

Using Monte Carlo simulation, we have compared the magnetic properties between nanostructured thin films and two-dimensional crystalline solids. The dependence of nanostructured properties on the interaction between particles that constitute the nanostructured thin films is also studied. The result shows that the parameters in the interaction potential have an important effect on the properties of nanostructured thin films at the transition temperatures.     

Monte Carlo simulation of grain growth in polycrystalline materials

Understanding the kinetics of grain growth, under the influence of second phase (such as impurities, voids and bubbles) is fundamental to advances in the control of microstructural evolution. As a precursor to this objective, we have investigated the grain growth kinetics in a polycrystalline material using a standard Q-state Potts’ model under Monte Carlo settings. Based on physical reasoning, new modifications are suggested to circumvent some of the disadvantages in the basic Potts model. The efficacy of these modifications vis-à-vis the basic model is verified. The influence of second phase particles on the impurity loaded grain boundaries is investigated for the study of grain growth kinetics.     

Monte Carlo simulation of glandular dose in a dedicated breast CT system

A dedicated breast CT system (DBCT) is a new method for breast cancer detection proposed in recent years. In this paper, the glandular dose in the DBCT is simulated using the Monte Carlo method. The phantom shape is half ellipsoid, and a series of phantoms with different sizes, shapes and compositions were constructed. In order to optimize the spectra, monoenergy X-ray beams of 5-80 keV were used in simulation. The dose distribution of a breast phantom was studied: a higher energy beam generated more uniform distribution, and the outer parts got more dose than the inner parts. For polyenergtic spectra, four spectra of Al filters with different thicknesses were simulated, and the polyenergtic glandular dose was calculated as a spectral weighted combination of the monoenergetic dose.     

Monte Carlo法模拟油滴的理论分布

在密立根油滴实验中,应用Monte Carlo法数值地模拟实验室条件下可能的油滴带电量Q、平衡电压U和匀速运动时间t的理论分布,分别从三维分布、时间分布及数量方面比较3种选择油滴的标准,给出较合理的结论.     

血液中Na活化放射性测量的蒙特卡罗模拟

针对超临界事故中人体受到中子辐照后感生的24Na活度测量,采用MCNP软件建立蒙特卡罗模拟模型,分别模拟不同类型NaI探测器对24Na衰变的两条射线全能峰的探测效率和塑料闪烁体探测器对24Na衰变的射线总计数的探测器效率。模拟结果表明：井型NaI探测器与圆柱型相比,24Na衰变的1.38 MeV和2.76 MeV 射线全能峰探测效率分别提高了4.30倍和4.11倍;塑料闪烁体探测器对24Na衰变的射线的探测效率是NaI探测器对24Na 射线的探测效率的1.72倍;同时粗略计算了探测器计数与人体所受中子辐照剂量之间的关系。     

Effective model for rare-earth Kitaev materials and its classical Monte Carlo simulation

Recently, the family of rare-earth chalcohalides were proposed as candidate compounds to realize the Kitaev spin liquid (KSL) [Chin. Phys. Lett. 38 047502 (2021)]. In the present work, we firstly propose an effective spin Hamiltonian consistent with the symmetry group of the crystal structure. Then we apply classical Monte Carlo simulations to preliminarily study the model and establish a phase diagram. When approaching to the low temperature limit, several magnetic long range orders are observed, including the stripe, the zigzag, the antiferromagnetic (AFM), the ferromagnetic (FM), the incommensurate spiral (IS), the multi-Q, and the 120° ones. We further calculate the thermodynamic properties of the system, such as the temperature dependence of the magnetic susceptibility and the heat capacity. The ordering transition temperatures reflected in the two quantities agree with each other. For most interaction regions, the system is magnetically more susceptible in the ab-plane than in the c-direction. The stripe phase is special, where the susceptibility is fairly isotropic in the whole temperature region. These features provide useful information to understand the magnetic properties of related materials.     

血液中Na活化放射性测量的蒙特卡罗模拟

针对超临界事故中人体受到中子辐照后感生的24Na活度测量,采用MCNP软件建立蒙特卡罗模拟模型,分别模拟不同类型NaI探测器对24Na衰变的两条射线全能峰的探测效率和塑料闪烁体探测器对24Na衰变的射线总计数的探测器效率。模拟结果表明:井型NaI探测器与圆柱型相比,24Na衰变的1.38 MeV和2.76 MeV 射线全能峰探测效率分别提高了4.30倍和4.11倍;塑料闪烁体探测器对24Na衰变的射线的探测效率是NaI探测器对24Na 射线的探测效率的1.72倍;同时粗略计算了探测器计数与人体所受中子辐照剂量之间的关系。     

Evaluation of medical exposure and exposure by the public in a typical scenario of examinations using mobile X-ray equipment through the Monte Carlo simulation

In this work irradiation scenarios that simulated chest and abdomen examinations involving mobile X-ray equipment in hospitals were modeled with the purpose of calculating conversion coefficient for effective dose (CC_E), normalized to entrance surface dose (ESD), applied to patients and public individuals. These coefficients can easily be used in this practice. Patients and public individuals were represented by a pair of anthropomorphic phantoms inserted in the MCNPX 2.7.0 radiation transport code. One of the phantoms (patient) was irradiated with the direct beam simulating examinations of the chest and abdomen, each with two fields of irradiation, ideal (IF) and extrapolated (EF). Using the software SPECGEN X-ray spectra from 60 to 100 kVp at 10 kVp intervals were generated and used in this work. The other phantom (public individual) was positioned 50–200 cm from the patient. In relation to the CC_E calculated in the patient, the average increase obtained between the irradiation fields was 62.4% for the chest examinations, and for the same conditions the CC_E was calculated for abdomen examinations and found to be 8.0%. Increasing the distance between public individual and patient, reductions of up to 81.7% in the CC_E in abdomen examinations and 83.4% in chest examinations were observed. Through the assessment of CC_E of these scenarios, it is possible to measure the damages relating to this practice for both patients and public individuals.     

Monte Carlo simulation using the Fourier transform of the interatomic potential

When doing Monte Carlo simulations using continuous potentials, the evaluation of the configurational potential energy ink-space by Fourier transformation is shown to be a computationally attractive scheme for systems where the long-range interatomic interaction spans a dimension comparable to the size of the simulated system.     

The Convergence of Markov Chain Monte Carlo Methods: From the Metropolis Method to Hamiltonian Monte Carlo

From its inception in the 1950s to the modern frontiers of applied statistics, Markov chain Monte Carlo has been one of the most ubiquitous and successful methods in statistical computing. The development of the method in that time has been fueled by not only increasingly difficult problems but also novel techniques adopted from physics. Here, the history of Markov chain Monte Carlo is reviewed from its inception with the Metropolis method to the contemporary state‐of‐the‐art in Hamiltonian Monte Carlo, focusing on the evolving interplay between the statistical and physical perspectives of the method.     

Monte Carlo simulation of the enantioseparation process

By means of Monte Carlo simulation, a study of enantioseparation by capillary electrophoresis has been carried out. A simplified system consisting of two enantiomers S (R) and a selector chiral C, which reacts with the enantiomers to form complexes RC (SC), has been considered. The dependence of Δμ$\Delta \mu$ (enantioseparation) with the concentration of chiral selector and with temperature have been analyzed by simulation. The effect of the binding constant and the charge of the complexes are also analyzed. The results are qualitatively satisfactory, despite the simplicity of the model.     

气泡幕后向散射光信号特性的蒙特卡罗方法研究

利用蒙特卡罗方法对水介质中的气泡幕的后向光散射回波信号进行了系统仿真。通过对计算结果与实验结果的比较,证明了蒙特卡罗方法的有效性。根据仿真结果,结合理论分析的方法,对气泡幕位置、厚度和接收器视场角等参数对回波信号的影响进行了分析,结果表明:回波信号出现时间与气泡幕位置存在一一对应的关系;在气泡幕的衰减系数ρσt不变的情况下,存在一个有效气泡幕厚度,当大于该厚度时,气泡幕的后向光散射回波信号基本不变;在系统各项参数不变的情况下,适当增加接收器的视场角,可有效地提高回波信号的信噪比。     

MacPherson-Srolovitz晶粒长大速率方程的仿真验证

2007年MacPherson和Srolovitz联合推导出一个三维个体晶粒长大的准确速率方程,但并未给出实验或计算机仿真的验证.采用Potts模型Monte Carlo方法对该速率方程进行了大尺度仿真验证.结果表明,仿真数据与MacPherson-Srolovitz速率方程符合很好,从而初步证实了该速率方程,即三维晶粒长大速率是晶粒棱长和晶粒平均宽度的函数. 关键词： 三维晶粒长大 速率方程 Monte Carlo仿真     

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

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

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

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

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.     

241Am-Be中子参考辐射的蒙特卡罗模拟研究

研究用于校准场所中子剂量监测仪表的241Am-Be中子参考辐射场计量特性。采用蒙特卡罗方法模拟了空气自由中子参考辐射(FRNR),GB/T 14055规定的最小尺寸中子参考辐射(SRNR)和实际中子参考辐射(ARNR)中不同检验点处中子周围剂量当量率、散射中子占比和能谱分布特征。研究结果表明,空气对FRNR中的剂量率和能谱分布影响小,近似为理想中子参考辐射;采用5%含硼聚乙烯作屏蔽的最小尺寸SRNR可减少热中子,降低散射中子占比,影锥法不适用于小尺寸中子参考辐射中对散射中子的修正;ARNR中的散射中子更少、占比更低,影锥法所得散射中子占比与理论值基本一致。     

Ni-15%Mo合金熔体热物理性质的Monte Carlo模拟

采用EAM势对6×6×6的Ni-15%Mo合金熔体进行Monte Carlo模拟,通过对不同温度下获得的NVT系统的平衡态统计分析得出Ni-15%Mo合金熔体在过冷态和过热态时的热物理性质.通过构造系统生成新表面,表面张力做功使系统能量发生改变,从而得到液态表面张力的模拟结果.Ni-15%Mo合金熔体的表面张力在1500—2000 K的温度范围内,随温度的变化规律为σ=1.918-1.130×10^-3（T-T_m） N/m 关键词： Monte Carlo模拟 表面张力 比热 Ni-15%Mo合金