Exhaustive review of acceleration strategies for Monte Carlo simulations in photon transit

Monte Carlo simulation techniques have become the quintessence and a pivotal nexus of inquiry in the realm of simulating photon movement within biological fabrics.Through the stochastic sampling of tissue archetypes delineated by explicit optical characteristics,Monte Carlo simulations possess the theoretical capacity to render unparalleled accuracy in the depiction of exceedingly intricate phenomena.Nonetheless,the quintessential challenge associated with Monte Carlo simulation methodologies resides in their extended computational duration,which significantly impedes the refinement of their precision.Consequently,this discourse is specifically dedicated to exploring innovations in strategies and technologies aimed at expediting Monte Carlo simulations.It delves into the foundational concepts of various acceleration tactics,evaluates these strategies concerning their speed,accuracy,and practicality,and amalgamates a comprehensive overview and critique of acceleration methodologies for Monte Carlo simulations.Ultimately,the discourse envisages prospective trajectories for the employment of Monte Carlo techniques within the domain of tissue optics.

Health risk assessment of trace metal(loid)s in agricultural soils based on Monte Carlo simulation coupled with positive matrix factorization model in Chongqing, southwest China

This study aimed to investigate the pollution characteristics, source apportionment, and health risks associated with trace metal(loid)s(TMs) in the major agricultural producing areas in Chongqing, China. We analyzed the source apportionment and assessed the health risk of TMs in agricultural soils by using positive matrix factorization(PMF) model and health risk assessment(HRA) model based on Monte Carlo simulation. Meanwhile, we combined PMF and HRA models to explore the health risks of TMs in agricultural soils by different pollution sources to determine the priority control factors. Results showed that the average contents of cadmium(Cd), arsenic (As), lead(Pb), chromium(Cr), copper(Cu), nickel(Ni), and zinc(Zn) in the soil were found to be 0.26, 5.93, 27.14, 61.32, 23.81, 32.45, and 78.65 mg/kg, respectively. Spatial analysis and source apportionment analysis revealed that urban and industrial sources, agricultural sources, and natural sources accounted for 33.0%, 27.7%, and 39.3% of TM accumulation in the soil, respectively. In the HRA model based on Monte Carlo simulation, noncarcinogenic risks were deemed negligible(hazard index <1), the carcinogenic risks were at acceptable level(10^(-6)<total carcinogenic risk ≤ 10^(-4)), with higher risks observed for children compared to adults. The relationship between TMs, their sources, and health risks indicated that urban and industrial sources were primarily associated with As, contributing to 75.1% of carcinogenic risks and 55.7% of non-carcinogenic risks, making them the primary control factors. Meanwhile, agricultural sources were primarily linked to Cd and Pb, contributing to 13.1% of carcinogenic risks and 21.8% of non-carcinogenic risks, designating them as secondary control factors.

Machine learning-enhanced Monte Carlo and subset simulations for advanced risk assessment in transportation infrastructure

The maintenance of safety and dependability in rail and road embankments is of utmost importance in order to facilitate the smooth operation of transportation networks.This study introduces a comprehensive methodology for soil slope stability evaluation,employing Monte Carlo Simulation(MCS)and Subset Simulation(SS)with the"UPSS 3.0 Add-in"in MS-Excel.Focused on an 11.693-meter embankment with a soil slope(inclination ratio of 2H:1V),the investigation considers earthquake coefficients(kh)and pore water pressure ratios(ru)following Indian zoning requirements.The chance of slope failure showed a considerable increase as the Coefficient of Variation(COV),seismic coefficients(kh),and pore water pressure ratios(ru)experienced an escalation.The SS approach showed exceptional efficacy in calculating odds of failure that are notably low.Within computational modeling,the study optimized the worst-case scenario using ANFIS-GA,ANFIS-GWO,ANFIS-PSO,and ANFIS-BBO models.The ANFIS-PSO model exhibits exceptional accuracy(training R2=0.9011,RMSE=0.0549;testing R2=0.8968,RMSE=0.0615),emerging as the most promising.This study highlights the significance of conducting thorough risk assessments and offers practical insights into evaluating and improving the stability of soil slopes in transportation infrastructure.These findings contribute to the enhancement of safety and reliability in real-world situations.

Optimization of the Use of Spherical Targets for Point Cloud Registration Using Monte Carlo Simulation

Registrations based on the manual placement of spherical targets are still being employed by many professionals in the industry.However,the placement of those targets usually relies solely on personal experience without scientific evidence supported by numerical analysis.This paper presents a comprehensive investigation,based on Monte Carlo simulation,into determining the optimal number and positions for efficient target placement in typical scenes consisting of a pair of facades.It demonstrates new check-up statistical rules and geometrical constraints that can effectively extract and analyze massive simulations of unregistered point clouds and their corresponding registrations.More than 6×10^(7) sets of the registrations were simulated,whereas more than IOO registrations with real data were used to verify the results of simulation.The results indicated that using five spherical targets is the best choice for the registration of a large typical registration site consisting of two vertical facades and a ground,when there is only a box set of spherical targets available.As a result,the users can avoid placing extra targets to achieve insignificant improvements in registration accuracy.The results also suggest that the higher registration accuracy can be obtained when the ratio between the facade-to-target distance and target-to-scanner distance is approximately 3:2.Therefore,the targets should be placed closer to the scanner rather than in the middle between the facades and the scanner,contradicting to the traditional thought. Besides,the results reveal that the accuracy can be increased by setting the largest projected triangular area of the targets to be large.

Monte Carlo Simulations of Doping Properties of a Spin-3/2 Ising Nanotube

The effect of spin-1 impurities doping on the magnetic properties of a spin-3/2 Ising nanotube is investigated using Monte Carlo simulations within the Blume-Emery-Griffiths model in the presence of an external magnetic field. The thermal behaviors of the order parameters and different macroscopic instabilities as well as the hysteretic behavior of the material are examined in great detail as a function of the dopant density. It is found that the impurities concentration affects all the system magnetic properties generating for some specific values, compensation points and multi-cycle hysteresis. Doping conditions where the saturation/remanent magnetization and coercive field of the investigated material can be modified for permanent or soft magnets synthesis purpose are discussed.

Two Monte Carlo-based simulators for imaging-system modeling and projection simulation of flat-panel X-ray source

The advantages of a flat-panel X-ray source(FPXS)make it a promising candidate for imaging applications.Accurate imaging-system modeling and projection simulation are critical for analyzing imaging performance and resolving overlapping projection issues in FPXS.The conventional analytical ray-tracing approach is limited by the number of patterns and is not applicable to FPXS-projection calculations.However,the computation time of Monte Carlo(MC)simulation is independent of the size of the patterned arrays in FPXS.This study proposes two high-efficiency MC projection simulators for FPXS:a graphics processing unit(GPU)-based phase-space sampling MC(gPSMC)simulator and GPU-based fluence sampling MC(gFSMC)simulator.The two simulators comprise three components:imaging-system modeling,photon initialization,and physical-interaction simulations in the phantom.Imaging-system modeling was performed by modeling the FPXS,imaging geometry,and detector.The gPSMC simulator samples the initial photons from the phase space,whereas the gFSMC simulator performs photon initialization from the calculated energy spectrum and fluence map.The entire process of photon interaction with the geometry and arrival at the detector was simulated in parallel using multiple GPU kernels,and projections based on the two simulators were calculated.The accuracies of the two simulators were evaluated by comparing them with the conventional analytical ray-tracing approach and acquired projections,and the efficiencies were evaluated by comparing the computation time.The results of simulated and realistic experiments illustrate the accuracy and efficiency of the proposed gPSMC and gFSMC simulators in the projection calculation of various phantoms.

Chord length sampling correction analysis for dispersion fuel in Monte Carlo simulation

Dispersion fuels,knowned for their excellent safety performance,are widely used in advanced reactors,such as hightemperature gas-cooled reactors.Compared with deterministic methods,the Monte Carlo method has more advantages in the geometric modeling of stochastic media.The explicit modeling method has high computational accuracy and high computational cost.The chord length sampling(CLS)method can improve computational efficiency by sampling the chord length during neutron transport using the matrix chord length?s probability density function.This study shows that the excluded-volume effect in realistic stochastic media can introduce certain deviations into the CLS.A chord length correction approach is proposed to obtain the chord length correction factor by developing the Particle code based on equivalent transmission probability.Through numerical analysis against reference solutions from explicit modeling in the RMC code,it was demonstrated that CLS with the proposed correction method provides good accuracy for addressing the excludedvolume effect in realistic infinite stochastic media.

Enhancing Precision in Radiotherapy Delivery: Validating Monte Carlo Simulation Models for 6 MV Elekta Synergy Agility LINAC Photon Beam Using Two Models of the GAMOS Code

The most crucial requirement in radiation therapy treatment planning is a fast and accurate treatment planning system that minimizes damage to healthy tissues surrounding cancer cells. The use of Monte Carlo toolkits has become indispensable for research aimed at precisely determining the dose in radiotherapy. Among the numerous algorithms developed in recent years, the GAMOS code, which utilizes the Geant4 toolkit for Monte Carlo simula-tions, incorporates various electromagnetic physics models and multiple scattering models for simulating particle interactions with matter. This makes it a valuable tool for dose calculations in medical applications and throughout the patient’s volume. The aim of this present work aims to vali-date the GAMOS code for the simulation of a 6 MV photon-beam output from the Elekta Synergy Agility linear accelerator. The simulation involves mod-eling the major components of the accelerator head and the interactions of the radiation beam with a homogeneous water phantom and particle information was collected following the modeling of the phase space. This space was po-sitioned under the X and Y jaws, utilizing three electromagnetic physics mod-els of the GAMOS code: Standard, Penelope, and Low-Energy, along with three multiple scattering models: Goudsmit-Saunderson, Urban, and Wentzel-VI. The obtained phase space file was used as a particle source to simulate dose distributions (depth-dose and dose profile) for field sizes of 5 × 5 cm2 and 10 × 10 cm2 at depths of 10 cm and 20 cm in a water phantom, with a source-surface distance (SSD) of 90 cm from the target. We compared the three electromagnetic physics models and the three multiple scattering mod-els of the GAMOS code to experimental results. Validation of our results was performed using the gamma index, with an acceptability criterion of 3% for the dose difference (DD) and 3 mm for the distance-to-agreement (DTA). We achieved agreements of 94% and 96%, respectively, between simulation and experimentation for the three electromagnetic physics models and three mul-tiple scattering models, for field sizes of 5 × 5 cm2 and 10 × 10 cm2 for depth-dose curves. For dose profile curves, a good agreement of 100% was found between simulation and experimentation for the three electromagnetic physics models, as well as for the three multiple scattering models for a field size of 5 × 5 cm2 at 10 cm and 20 cm depths. For a field size of 10 × 10 cm2, the Penelope model dominated with 98% for 10 cm, along with the three multiple scattering models. The Penelope model and the Standard model, along with the three multiple scattering models, dominated with 100% for 20 cm. Our study, which compared these different GAMOS code models, can be crucial for enhancing the accuracy and quality of radiotherapy, contributing to more effective patient treatment. Our research compares various electro-magnetic physics models and multiple scattering models with experimental measurements, enabling us to choose the models that produce the most reli-able results, thereby directly impacting the quality of simulations. This en-hances confidence in using these models for treatment planning. Our re-search consistently contributes to the progress of Monte Carlo simulation techniques in radiation therapy, enriching the scientific literature.

Evaluation of the Reliability of a System: Approach by Monte Carlo Simulation and Application

The objective of this paper is to evaluate the reliability of a system in its different states (absence of failures, partial failure and total failure) and to propose actions to improve this reliability by an approach based on Monte Carlo simulation. It consists of a probabilistic evaluation based on Markov Chains. In order to achieve this goal, the functionalities of Markov Chains and Monte Carlo simulation steps are deployed. The application is made on a production system. .

Evaluation of Particle Properties of MgO/TiO2 Material by Monte Carlo Simulation Method

The simulation by the Monte Carlo method executed by the software PyPENELOPE proved effective to specify the particle propagation characteristics by calculating the absorption fractions, backscattering and transmission of electrons and secondary photons under the incidence of 0.5 to 20 KeV range of primary electrons. More than 99.9% of the primary electrons were transmitted in the 125 nm thick MgO/TiO2 material at 20 KeV. This occurred because several interactions took place in the transmitted primary irradiation such as characteristic, fluorescence, and bremsstrahlung produced when of the occupation of the KL3, KL2, KM3, and KM2 shell and sub-shell of titanium and magnesium which are the elements with a high atomic number in the material. The transmission particle characteristic of this material is therefore an indicator capable of improving the electrical performance and properties of the sensor.

Study on kinetics of propylene polymerizationat different temperatures via Monte Carlo simulation

The elementary reactions of propylene polymerization catalyzed by conventional Ziegler-Natta catalysts was proposed according to the comprehensive view and without considering the effect of any impurity in the material on propylene polymerization. The Monte Carlo simulation technique was employed to investigate the kinetics of propylene polymerization in order to determine the validity of the stationary state assumption and the effects of the polymerization temperature on the polymerization. The simulated total amount of active species, which only increases quickly at the beginning of the polymerization, indicates that the stationary state assumption in the studied system is valid. Moreover, significant effects of polymerization temperature on the polymerization conversion, and the molecular weight and its distribution were also analyzed. The simulated results show that the consumption rate of propylene increases with the increase of polymerization temperature; the maximum values of the number-average degree of polymerization are constant at different polymerization temperatures, however, the peak appears earlier with the higher temperature; as the polymerization temperature increases, the average molecular weight decreases and the molecular weight distribution changes greatly.

MONTE CARLO SIMULATION OF RADIATION FIELD OPTIMIZATION FOR MEDICAL LINAC

A method for designing an X-ray flatness filter for medical electron linac is developed. It is used in the optimization process in the electron beam radiation system. Monte Carlo simulation method is used and two examples of real radiation system optimization processes for China-made medical electron linac are provided： 15 MV X- ray system of BJ-20 linac, and 12 MeV electron system of BJ-14. Results are verified by using the traditional method.

Research and development of the automatic modeling system for Monte Carlo particle transport simulation

在广泛调研和分析现有几何建模方法特点的基础上研发了具有可视化用户界面的自动建模程序系统MCAM。它可以实现多种商用软件CAD模型与McNP模型之间的相互转换，且提供了模型建立、预处理、属性分析等基本功能和计算结果可视化及基于医学映像建模接口等扩展功能。全面系统地介绍了MCAM的设计思想与原理、总体结构、主要功能和国际合作协议框架下的应用测试等情况。实践表明，它是一个实用的MCNP计算辅助工具和核设计与核分析质量保证工具。

Comparison Between Three-Valley Model and Full Band Model in Monte Carlo Simulation of Bulk Wurtzite GaN

The Monte Carlo simulators with the three valley model and the full band Monte Carlo model are used to explore electron transport in bulk wurtzite gallium nitride (GaN).Comparison of the results based on the two models is made.The results based on both models are basically the same at the lower field region,but exhibit some differences at the higher field region.The electron average energy exhibits obvious difference at the high field region between the two models.This difference further causes several other differences of GaN properties,such as the drift velocity versus field characteristics,the repopulation.Because of the complicated energy band structures at the high energy region for wurtzite GaN,the analytical band structures in the three valley model can not cover all properties of the band structures of wurtzite GaN,so the results based on the full band Monte Carlo model should be more exact.

Monte Carlo numerical simulation and its application in probability analysis of long span bridges

To get the probability of long span bridges under the influence of external random factors, the Monte Carlo method using Latin hypercube sampling is applied. Combined with the condition assessment system on Runyang Suspension Bridge, which is the longest suspension bridge in China, the structural probabilities in normal and damaged situations are calculated with the external random factors considered including environmental temperature, wind load, load of vehicles, etc. The main assessment items contain the maximal vertical displacement of girder, the maximal stress of cables, the maximal horizontal displacement of towers etc. Finally, the probabilities and their cumulative distribution functions are provided. The analysis results can be plotted on line in a clear and vivid way, so the efficiency of assessment is increased and the decision-making of maintenance is more objective and accurate.

Monte Carlo方法在移动通信可靠性分析中的应用

文章主要探讨Monte Carlo方法在移动通信可靠性分析中的应用,并引入重要抽样技术对Monte Carlo方法进行优化。首先,文章介绍Monte Carlo方法的基本原理及其在移动通信可靠性分析中的应用。其次,引入重要抽样技术对该方法进行优化。最后,利用优化网络工程工具(Optimized Network Engineering Tools,OPNET)模拟平台构建相应的移动通信网络模型,并对比Monte Carlo方法和本文研究的优化方法在不同样本数量下的性能差异。结果表明,重要抽样技术的引入显著提高了Monte Carlo方法在移动通信可靠性分析中的准确性和效率,减少了方差,验证了该方法在移动通信可靠性分析中的有效性和实用性。

Full Band Monte Carlo Simulation of Electron Transport in Ge with Anisotropic Scattering Process

The electron transport properties in Ge are calculated by full band Monte Carlo technique with anisotropic scattering consideration.The calculation procedures are as follows:the full band structure is calculated by nonlocal empirical pseudopotential approach;the relative value of density of state (DOS) is computed by counting the number of states located in a certain region of the energy;the phonon dispersion curve is obtained from an adiabatic bond charge model;the electron phonon scattering rates are approximated by the nonparabolic model derived from Fermi’s golden rule at low energy region and scaled by DOS at higher energy region;the energy and momentum conservations are employed for choosing the final state after scattering.The validity of this Monte Carlo simulator and the physical models that are used is fully confirmed by comparing the program output to experimental results listed in references.As this Monte Carlo model can accurately reproduce the velocity and energy characteristics of electrons in Ge and the DOS scaled scattering rate can significantly reduce the computational cost for scattering rates,this approach is suitable for device simulation.

Monte Carlo Simulation of Impurity Scattering Effect in Resonant-Phonon-Assisted Terahertz Quantum-Cascade Lasers

We study the influence of ionized impurity scattering on the electron transport in resonant-phonon-assisted terahertz （THz） quantum-cascade lasers （QCLs）. We treat the ionized impurity scattering rates within the single subband static screening approximation. We find that the ionized impurity scattering supplies an additional current channel across the device,and affects the electrondistribution in different subbands. We conclude that the ionized impurity scattering should be taken into account in the study of the transport properties of resonant-phonon-assisted THz QCLs.

基于Monte-Carlo仿真的恒应力加速试验方案模拟评价

为了保证可靠性加速试验方案的正确性和可行性,提出基于仿真的模拟评价方法。以最具代表性的Weibull寿命分布产品的恒应力加速寿命试验为例,以代表试验结果准确性和稳定性的指标——中位寿命估计值渐进方差的均值和标准差作为评价指标,通过Monte-Carlo仿真模拟生成失效数据,对不同试验方案的统计精度和稳定性进行评价,为加速寿命试验优化设计的工程实现提供了思路和实例验证。

Using Monte Carlo Simulation Technology to Improve Intuitive Effect of Teaching Probability and Mathematical Statistics Course

With the illustration of a specific problem, this paper demonstrates that using Monte Carlo Simulation technology will improve intuitive effect of teaching Probability and Mathematical Statistics course, and save instructors' effort as well.And it is estimated that Monte Carlo Simulation technology will be one of the major teaching methods for Probability and Mathematical Statistics course in the future.