应用蒙特卡罗模拟进行正电子发射断层成像仪散射特性分析

通过使用基于Geant4的蒙特卡罗模拟代码GATE对全身用全3D正电子发射断层成像仪和含有挡板的2D/3D小动物用正电子发射断层成像仪进行建模,系统地分析了3D采集条件下正电子发射断层成像仪的散射分数、散射分布、多次散射、视野外散射四个主要方面和2D采集条件下挡板对散射分数和散射分布的影响.针对全3D散射校正的难点: 多次散射和视野外散射,设计了附加实验,拟合得到了多次散射光子的百分比随体模横截面积变化的关系和不同环的位置受到视野外散射光子的影响;针对2D散射校正,对挡板引入的散射光子进行分离,单独分析, 关键词： 正电子发射断层成像仪（PET） 蒙特卡罗模拟 散射特性 散射校正     

用于HIRFL重离子治癌的PET成像实验研究与蒙特卡罗模拟

为了HIRFL重离子治癌的需要,使用两个位置灵敏闪烁体探测器组成了一个简单的PET成像系统.对PET成像进行了实验研究,实验测量得到了物体成像以及γ射线能谱,对于511keV全能峰处的能量分辨率为186%,峰总比为52.4%.采用GEANT3程序对该系统进行了模拟计算,并与实验进行比较,最后对该系统进行了优化设计     

光学相干层析信号的模拟分析与计算

对传统的模拟介质散射响应函数的蒙特卡罗算法的计算过程进行了改进，提出光子不可分割的假设，这样可以避免人为选定光子被吸收阈值的不确定性；同时采用双曲线法模拟入射高斯光束，解决了对表面为非平面介质的传输信号进行模拟的问题.通过对盖玻片（平面）与眼角膜（曲面）的光学相干层析探测信号及层析图的数值模拟，并与相应实验结果进行比较，证实了该方法与假设的可行性. 关键词： 光学相干层析 蒙特卡罗法模拟 散射     

Thermodynamic properties of ternary alloys from Monte Carlo simulations

Methods are presented for obtaining both the integral and partial molar energies, entropies and free energies of the components in rigid ternary substitutional alloy systems from the one Monte Carlo simulation. Tests of the methods for some model systems are presented and discussed.     

材料中He深度分布演化的Monte Carlo模拟研究

基于He泡生长的迁移-合并机理,用Monte Carlo方法模拟了对材料进行等温退火过程中He深度分布的演化,探讨了不同参数对这一演化的影响.研究表明:材料中He泡的初始浓度和尺寸将影响He深度分布的变化,而退火温度则对演化速率起重要作用但对最终的He深度分布影响较小;随着反应的进行,整个系统的演化是逐渐趋缓的. 关键词： He 深度分布 Monte Carlo模拟     

First-order transition in a 2D classical XY-model using microcanonical Monte Carlo simulations

We have simulated two-dimensional classical XY-model in a microcanonical ensemble using the Monte Carlo technique. Simulations were carried out on a square lattice having 25, 100 or 900-spins with periodic boundary conditions. The nearest neighbour interaction potential was taken to beV(θ)=2J[1−cos¹⁰⁰(θ/2)]. The system energy, mean square magnetization and vortex-density were calculated as functions of temperature. A sudden change in the system energy, vortex density and mean square magnetization was observed at the first-order transition which is associated with this choice of the nearest neighbour interaction potential. The transition temperature increases with decrease in the system size. It is found that the creation of a vortex-antivortex pair costs more energy during the first-order transition than the energy associated with a tightly bound vortex-antivortex pair.     

蒙特卡罗法计算无序激光器的阈值

提出了计算无序激光器阈值增益的模型,提出了已知平均自由程和散射强度分布得到分布函数及相应的随机变量的方法.对各向同性的光散射,导出了各分布函数,并在每次散射中得到相应的随机变量.运用蒙特卡罗法模拟光子在无序介质中的随机行走,得到了无序激光器的阈值.讨论了阈值增益与光子平均自由程、增益区半径、增益区深度、非增益区光吸收系数的关系,提出了降低阈值的方法.理论结果与实验数据是符合的. 关键词： 激光物理 阈值 蒙特卡罗 散射     

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.     

Monte Carlo simulations of jet quenching in heavy ion collisions

Jet quenching is one of the major discoveries of the heavy-ion program at Rhic. While there is a wealth of data from Rhic that will soon be supplemented with measurements at the Lhc, on the theoretical side the situation is less clear. A thorough understanding of jet quenching is, however, beneficial, as it is expected that medium-induced modifications of jets allow one to characterise properties of the QCD matter produced in heavy ion collisions. This talk aims at summarising the main ideas and concepts of the currently available Monte Carlo models for jet quenching.     

聚合物中多重光散射传导的Monte Carlo数值模拟

基于米氏散射(Mie scattering)理论,建立填充分散粒子群的聚合物对光散射传导的Monte Carlo数学模型.在此基础上,编写了一套仿真模拟程序.通过模拟单个光子在聚合物中的多重散射运动过程,把问题扩展到以激光束或线状光为入射源,得到在聚合物板块内的光传导情况,并且在计算机上图像化地重现整个物理过程,对输出光强的分布情况进行模拟统计分析.模拟结果表明,利用体散射机制,可以将点光源和线光源转换为平面光输出,输出光的状态可以通过对比计算结果实施有效控制.     

TLD measurements and Monte Carlo simulations for glandular dose and scatter fraction assessment in mammography: A comparative study

The main purpose of this study was to validate and compare Mean Glandular Dose (MGD) values obtained using Monte Carlo simulations with experimental values obtained from Entrance Surface Dose (ESD) and depth dose measurements performed in a Hospital mammography unit. ESD and depth dose were measured using ThermoLuminescent Dosimeters (TLDs), and a tissue equivalent mammography phantom recommended by the American College of Radiology (ACR). Measurements and Monte Carlo simulations were also compared with the MGD calculated using the Automatic Exposure Control (AEC) system of the mammographic unit. In the simulations the Doppler energy broadening effect was also taken into account. The simulated ESD are about 5%–10% higher than the measured ESD values. The deviation between the measured and simulated MGD values in the phantom is of about 15%. The MGD evaluated using the AEC system is smaller both with respect to the Monte Carlo simulation and experimental result by a factor of about 15% and 25% respectively. Moreover the BackScatter Factor (BSF) estimated by Monte Carlo simulations was used for the MGD calculation according to the Wu’s method. Finally the inclusion of the energy broadening effect on MGD calculation produces negligible variations on the simulated results.     

基于蒙特卡洛和器件仿真的单粒子翻转计算方法

文章提出了一种基于蒙特卡洛和器件仿真的存储器单粒子翻转截面获取方法,可以准确计算存储器单粒子效应,并定位单粒子翻转的灵敏区域.基于该方法,计算了国产静态存储器和现场可编程门阵列(FPGA)存储区的单粒子效应的截面数据,仿真结果和重离子单粒子效应试验结果符合较好.仿真计算揭示了器件单粒子翻转敏感程度与器件n,p截止管区域面积相关的物理机理,并获得了不同线性能量转移(LET)值下单粒子翻转灵敏区域分布.采用蒙特卡洛方法计算了具有相同LET、不同能量的离子径迹分布,结果显示高能离子的电离径迹半径远大于低能离子,而低能离子径迹中心的能量密度却要高约两到三个数量级.随着器件特征尺寸的减小,这种差别的影响将会越来越明显,阈值LET和饱和截面将不能完全描述器件单粒子效应结果.     

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.     

逆Compton散射的MonteCarlo模拟

研究模拟均匀、各向同性介质内的逆康普顿散射过程的Monte Carlo方法,分析电子速度分布乘抽样法和光子散射方向分布乘抽样法的抽样效率,完善电子速度分布的直接抽样法、光子散射能量分布的乘加抽样法.通过对这四种抽样方法抽样费用的分析和比较,得出各种方法的最适条件,该方法可以模拟更高能的辐射.     

Guaranteeing total balance in Metropolis algorithm Monte Carlo simulations

The condition of detailed balance has long been used as a proxy for the more difficult-to-prove condition of total balance, which along with ergodicity is required to guarantee convergence of a Markov Chain Monte Carlo (MCMC) simulation to the correct probability distribution. However, some simple-to-program update schemes such as the sequential and checkerboard Metropolis algorithms are known not to satisfy detailed balance for such common systems as the Ising model.     

Two-sided bounds on the free energy from local states in Monte Carlo simulations

We show that a precise assessment of free energy estimates in Monte Carlo simulations of lattice models is possible by using cluster variation approximations in conjunction with the local states approximations proposed by Meirovitch. The local states method (LSM) utilizes entropy expressions which recently have been shown to correspond to a converging sequence of upper bounds on the thermodynamic limit entropy density (i.e., entropy per lattice site), whereas the cluster variation method (CVM) supplies formulas that in some cases have been proven to be, and in other cases are believed to be, lower bounds. We have investigated CVM-LSM combinations numerically in Monte Carlo simulations of the two-dimensional Ising model and the two-dimensional five-states ferromagnetic Potts model. Even in the critical region the combination of upper and lower bounds enables an accurate and reliable estimation of the free energy from data of a single run. CVM entropy approximations are therefore useful in Monte Carlo simulation studies and in establishing the reliability of results from local states methods.     

Monte Carlo simulation of light scattering from size distributed sub-micron spherical CdS particles in a volume element

Simulation of polarized light scattering by spherical particles having modal radius of 180 nm is presented in this paper. A Monte Carlo method which is based on the Stokes-Mueller formalism developed in ANSI Standard C-language is used for simulation. Single scattering is considered in our program with monodispersed sub-micron sized spherical CdS particles. We have considered only θ dependent scattering as described by Mie theory for spherical CdS particles. The experiments for studying light scattering properties of these particles were conducted in a designed and developed laser based light scattering studies setup. The simulation results were compared with experimental results and theoretical results obtained purely from Mie theory. The closeness of agreement or disagreement between these results is discussed in this paper.     

分子串模型中空间弛豫模式的弛豫动力学的蒙特卡罗模拟

基于玻璃化转变的分子串模型的分子串弛豫方程,提出了更为精确的模拟分子串中所有空间弛豫模式(SRM)的蒙特卡罗模拟方案. 模拟得出各个SRM的弛豫时间随温度和分子串长度的变化结果与分子串模型中分子串弛豫方程的预言完全一致,即理论预期和模拟结果相互印证. 应当指出,分子串能否作为液态中集体单元的必要条件是在考虑到分子串之间的不均匀随机相互作用后,分子串的所有SRM的定性特征是不能改变的,这就需要对不同分子串的SRM之间的耦合进行研究. 但是迄今为止,仍未发现相关的严格解,仅有近似的自洽弛豫平均场方法. 由此可知,所提出的模拟方案为研究不同分子串的SRM之间的耦合(包括上述自洽场的可行性)提供了必要的基础. 关键词： 玻璃化转变 弛豫动力学 蒙特卡罗模拟 分子串     

低能电子在固体表面背散射系数的直接Monte Carlo方法模拟

应用单次碰撞的直接Monte Carlo方法计算能量范围从100 eV~10 keV的电子在固体Al,Si,Au表面的背散射系数,其中低能电子在固体中的弹性散射和非弹性散射截面分别应用Mott散射截面和Born近似下的广义振子强度计算模型得到.通过与压缩历史Monte Carlo方法的模拟计算结果及实验值的比较,结果表明,对于100 eV~10 keV范围的低能区电子,采用直接方法计算得到的电子背散射系数与实验值符合较好,直接方法比压缩历史方法更适合于能量在10 keV以下的电子输运计算.     

Monte Carlo simulations including energy from an entropic force

Several experimental techniques have shown that the primary response of many materials comes from a heterogeneous distribution of independently relaxing nanoscale regions; but most Monte Carlo simulations have homogeneous correlations. Resolving this discrepancy may require including the energy needed to change the configurational entropy, which is often used in theoretical treatments of thermal fluctuations, but not in computer simulations. Here the local configurational entropy is shown to give a nonlinear correction to the Metropolis algorithm that restores conservation of energy, maintains maximum entropy, and yields heterogeneous correlations. The nonlinear correction also improves agreement between Monte Carlo simulations of the Ising model and measurements of specific heat and structural correlations from the Jahn–Teller distortion in LaMnO₃.