基于蒙特卡罗方法的几何因子计算程序

应用蒙特卡罗方法求解几何因子,基于蒙特卡罗方法的几何因子计算程序使用C++语言编写,可用于任意位姿的各种尺寸的圆面探测器对圆面源几何因子的计算。该程序使用了方差减小技巧。通过与国际通用蒙特卡罗计算程序（MCNP5）的计算结果对比,该方法具有结果准确（误差较小）、计算速度快、使用方便等优点。最终使用该程序计算几何因子,与实验数据进行对比,成功验证了中子深度分布分析（NDP）能谱测量系统探测器位姿的准确性（误差5%以内）,并对其移动位置进行修正,发现电机移动20 mm大约会产生1 mm的误差。     

基于蒙特卡罗方法的几何因子计算程序

应用蒙特卡罗方法求解几何因子,基于蒙特卡罗方法的几何因子计算程序使用C++语言编写,可用于任意位姿的各种尺寸的圆面探测器对圆面源几何因子的计算。该程序使用了方差减小技巧。通过与国际通用蒙特卡罗计算程序(MCNP5)的计算结果对比,该方法具有结果准确(误差较小)、计算速度快、使用方便等优点。最终使用该程序计算几何因子,与实验数据进行对比,成功验证了中子深度分布分析(NDP)能谱测量系统探测器位姿的准确性(误差5%以内),并对其移动位置进行修正,发现电机移动20mm大约会产生1mm的误差。     

Monte Carlo calculation using a personal computer for electron probe microanalysis

Monte Carlo simulation for calculating x‐ray spectra under electron impact is carried out by a personal computer. Reduction of the CPU time by reducing the number of incident electrons for calculation is performed. One of the methods used to reduce the number of incident electrons to obtain the results within an acceptable error is provided. We limit the present report to five typical compounds: an oxide (Al₂O₃), a nitride (TiN), a carbide (Al₄C₃), a light–light element compound (BN), and a heavy–light element compound (PbS). We calculate the relative x‐ray intensity generated from the five compounds listed above with the change of the number of the incident electrons. The calculated x‐ray intensity is then converted into the mass % by the ZAF method. Comparisons are made between the calculation and experiment for Al₄C₃ and PbS. It is found that the present Monte Carlo method is useful for investigating a relationship between the number of incident electrons and the acceptable error. The results of 150 incident electrons can be applied in practical analysis, and the calculation time is less than 2 min. Copyright © 2004 John Wiley & Sons, Ltd.     

Monte Carlo integration on GPU

We use a graphics processing unit (GPU) for fast computations of Monte Carlo integrations. Two widely used Monte Carlo integration programs, VEGAS and BASES, are parallelized for running on a GPU. By using W ⁺ plus multi-gluon production processes at LHC, we test the integrated cross sections and execution time for programs written in FORTRAN and running in the CPU and those running on a GPU. The integrated results agree with each other within statistical errors. The programs run about 50 times faster on the GPU than on the CPU.     

基于蒙特卡罗方法的6 MV Truebeam剂量计算

使用蒙特卡罗方法快速准确地模拟6 MV Varian Truebeam医用电子直线加速器射野剂量特性,探究厂家提供相位空间源的可用性及模拟方法的准确性。以Varian公司提供的出束窗口位置处相位空间文件作为Beamnrc/EGSnrc输入源,模拟射野形成结构并计算10 cm10 cm射野下的均匀水体模中的剂量分布,将计算结果与相同条件下的实验测量数据进行比较。同一计算机模拟时,传统完整机头模拟方法需4~5 h,文中所述方法模拟时间可缩减至48 min,剂量计算结果与相同条件下所得实际测量数据可较好地吻合,两者的百分深度剂量差异低于3%,且建成区吻合良好;不同深度处80%等剂量线所包含的射野区域内百分离轴剂量比差异均低于3%,且半影区效果好。使用厂家提供相位空间文件作为EGSnrc入射源,能快速模拟治疗头,得到剂量计算结果与实际测量值的差异满足剂量计算精度要求。     

Monte Carlo integration for lattice gauge theories with fermions

A Monte Carlo procedure is constructed for lattice gauge theories with fermions by replacing integration over fermion degrees of freedom in the path integral by conventional integration over effective boson degrees of freedom. The method is applied to gauge theories over two discrete subgroups of SU(2).     

Quantum Monte Carlo calculation for the neutron-rich Ca isotopes

We computed ground-state energies of calcium isotopes from ⁴²Ca to ⁴⁸Ca by means of the Auxiliary Field Diffusion Monte Carlo (AFDMC) method. Calculations were performed by replacing the ⁴⁰Ca core with a mean-field self-consistent potential computed using the Skyrme interaction. The energy of the external neutrons is calculated by projecting the ground state from a wave function built with the single-particle orbitals computed in the self-consistent external potential. The shells considered were the 1F _7/2 and the 1F _5/2 . The Hamiltonian employed is semi-realistic and includes tensor, spin-orbit and three-body forces. While absolute binding energies are too deep if compared with experimental data, the differences between the energies for nearly all isotopes are in very good agreement with the experimental data.     

Monte Carlo calculation of the configurational entropy and combinatory factor on binary alloys

A method is presented for the calculation of free energy and entropy of regular solutions. It is based on a calculation of the configurational density of states (combinatory factor), from Monte Carlo computations. Application to the planar square lattice shows an excellent agreement with exact results. Results on a bcc lattice are presented.     

蒙卡燃耗中的氙振荡现象分析及计算研究

根据氙振荡产生的机理,按照燃耗步长来划分：在短步长条件下会产生物理氙振荡,在长步长条件下会产生数值氙振荡。研究发现,当燃耗步长增大到19 h,物理氙振荡的振幅最小,随着燃耗步长继续增大,表现为数值氙振荡,当燃耗步长27 d,数值氙振荡振幅最小。采用平衡氙方法强制使中子通量与氙浓度平衡,在堆芯物理计算程序(RMC)的输运模块中运用平衡氙方法对碘氙的浓度进行不断的迭代更新,在保证收敛的情况下抑制了氙振荡现象的发生。     

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.     

基于蒙特卡罗正算输运的全局减方差方法

蒙特卡罗方法是当前形势下辐射屏蔽计算的首选分析工具。小概率深穿透问题则是屏蔽计算的关键与亟待解决的核心问题,需要使用有效的减方差技巧。针对全局问题,利用蒙特卡罗正算输运得到的粒子通量或探测响应来构建权重窗参数,将现有的粒子位置偏移拓展到位置和能量偏倚。利用国际屏蔽基准题进行测试验证,通过使用该方法,粒子被引导到模型的所有位置。平均相对误差降低到10%以下,几乎所有网格区域都有粒子统计。结果表明,基于蒙特卡罗正算输运的输运偏倚参数构建方法能够实现全局减方差。     

Monte Carlo calculation of hot electron mobility in PbS

Electron mobilities in PbS at a lattice temperature of 77°K have been calculated by the Monte Carlo method. The effects of the band nonparabolicity and of the deformation potential acoustic and polar optical phonon scattering have been included in the calculations. The results are compared with experiments and with the ones obtained on the basis of a displaced Maxwellian distribution function.     

Monte Carlo calculation of energy loss of channeled ions

Energy distribution of 2 MeV He ions channeled in 〈100〉 Si were determined using two different models of energy loss. A comparison of the calculated results with recent experimental data on the angular dependence of the energy loss is presented.