Scattering of N2 from Ni(111)

A cold (T_rot<10 K) beam of N₂ with an initial translational energy of 0.40 eV strikes an Ni(111) surface at surface temperatures from 300 to 873 K at several incident angles from 15 to 60°. The rotational energy and angular distributions of the scattered molecules are probed using (2+1) resonance-enhanced multiphoton ionization. Molecules scattered in the specular direction have mean rotational energies that are independent of surface temperature, whereas those scattered at angles far from the specular show a dependence on surface temperature, caused likely by multiple collisions with the surface before escape. A rotational rainbow, seen in systems such as CO–Ni(111) and N₂–Ag(111), is not seen in this system. For molecules that scatter close to the specular direction, approximately 10% of the initial translational energy is converted into rotational energy of the scattered N₂. For surface temperatures above room temperature, the angular distributions indicate that molecules that scatter into low-J states also tend to exit in a broad peak (10–20° FWHM) near the specular, and this peak is broadened with increasing incident angle. The molecules that scatter into high-J states have a much broader distribution, indicating that they are trapped rotationally during the initial collision and suffer multiple collisions before leaving the surface.     

Mega-Electron-Volt Electron Scattering and Radiography of Plasma

A Monte Carlo code is developed to study mega-electron-volt(MeV) electron scattering and transport in plasma based on multiple scattering.A scaling law relating the angular width of a scattered beam to the incident electron energy and the areal density of plasma is found,which may provide a method of MeV electron radiography for diagnosing the areal density of high-temperature,dense plasma under fusion conditions.The study on the MeV electron beam radiography also shows that plasma density interfaces could be discriminated by electron scattering.     

基于蒙特卡罗方法的XHA600D医用电子直线加速器的入射电子束参数研究

使用蒙特卡罗方法研究入射电子束参数对XHA600D医用电子直线加速器产生的剂量分布的影响,并确定优化的入射电子束参数。根据厂商提供的XHA600D加速器治疗头的几何、材料参数,使用蒙特卡罗程序EGSnrc对不同的入射电子束参数进行模拟并记录其在水模体中产生的剂量分布,将模拟结果与测量结果进行比较。模拟的入射电子束参数包括平均能量、径向强度分布、角度展宽和能量展宽;实验测量数据包括4 cm×4 cm、10 cm×10 cm、30 cm×30 cm射野条件下的百分深度剂量与离轴剂量。结果表明当入射电子束的平均能量为6 MeV、径向强度的半高宽(Full Width at Half Maximum, FWHM)为0.25 cm、角度展宽为0.15°时,模拟结果和测量结果吻合非常好。这些参数可以作为建立适用于XHA600D加速器的TPS(Treatment Planning System)剂量计算模型的基础参数。     

Frictional cooling: Latest experimental results

Frictional cooling — that is cooling a beam of very low-energetic charged particles by moderation in matter and simultaneous acceleration in an electrostatic field — has been shown to be feasible during our experiments in 1994/1995 at PSI. In agreement with our previous closed-form and Monte Carlo calculations we found a significant increase in spectral density and a decrease in the angular spread in the case of a beam of negative muons.     

Validation of Monte Carlo simulation of 6 MV photon beam produced by Varian Clinac 2100 linear accelerator using BEAMnrc code and DOSXYZnrc code

The Monte Carlo model for the photon-beam output from the Varian Clinac 2100 linear accelerator was validated to compare the calculated to measured PDD and beam dose profiles The Monte Carlo calculation method is considered to be the most accurate method for dose calculation in radiotherapy. The objective of this study is to build a Monte Carlo geometry of Varian Clinac 2100 linear accelerator as realistically as possible. The Monte Carlo codes used in this work were the BEAMnrc code to simulate the photons beam and the DOSXYZnrc code to examinate the absorbed dose in the water phantom. We have calculated percentage depth dose (PDD) and beam profiles of the 6 MV photon beam for the 6 × 6 cm², 10 × 10 cm² and 15 × 15 cm² field sizes. We have used the gamma index technique for the quantitative evaluation to compare the measured and calculated distributions. Good agreement was found between calculated PDD and beam profile compared to measured data. The comparison was evaluated using the gamma index method and the criterions were 3% for dose difference and 3 mm for distance to agreement. The gamma index acceptance rate was more than 97% of both distribution comparisons PDDs and dose profiles and our results were more developed and accurate. The Varian Clinac 2100 linear accelerator was accurately modeled using Monte Carlo codes: BEAMnrc and DOSXYZnrc codes package.     

Study on the radiation problem caused by electron beam loss in accelerator tubes

The beam dynamic code PARMELA was used to simulate the transportation process of accelerating electrons in S-band SW linacs with different energies of 2.5, 6 and 20 MeV. The results indicated that in the ideal condition, the percentage of electron beam loss was 50% in accelerator tubes. Also we calculated the spectrum, the location and angular distribution of the lost electrons. Calculation performed by Monte Carlo code MCNP demonstrated that the radiation distribution of lost electrons was nearly uniform along the tube axis, the angular distributions of the radiation dose rates of the three tubes were similar, and the highest leaking dose was at the angle of 160° with respect to the axis. The lower the energy of the accelerator, the higher the radiation relative leakage. For the 2.5 MeV accelerator, the maximum dose rate reached 5% of the main dose and the one on the head of the electron gun was 1%, both of which did not meet the eligible protection requirement for accelerators. We adopted different shielding designs for different accelerators. The simulated result showed that the shielded radiation leaking dose rates fulfilled the requirement.     

MONTE CARLO STUDIES OF RADIATION DAMAGE INDUCED BY FUSION DYNAMIC SIMULATION OF SPUTTERING

In this paper, the sputtering caused by fusion of ⁴He particles in the first wall materials is investigated by using the dynamic Monte Carlo simulation, which is based on the binary collision approximation. The dependences of sputtering yields on the incident energy and angle, as well as the comparisons of results calculated using the Monte Carlo methed with results from experiments, are discussed. Energy spectrum and angular distributions of sputtered species, the depth of origin of the sputtered particles as well as range distributions of incident ions are given. From a comparison between the related experimental data and the calculated results with static and dynamic simulation it has been found that a better agreement is obtained for dynamic simulation.     

考虑中子泄漏修正的快堆蒙特卡罗少群截面参数制作与验证方法

一种基于B1均匀化方程的泄漏修正模型在连续能量蒙特卡罗程序TRIPOLI4中得以实现并且用于制作少群截面参数。此蒙卡泄漏修正模型通过在连续能量的蒙卡模拟以及求解B1均匀化方程之间迭代,最终得到蒙卡模拟下的临界状态。通过此方法得到的少群截面参数较其他蒙卡以及确定论方法有两点显著优势:用于求解B1均匀化方程的少群常数是用通过临界状态的通量谱得到的;考虑了泄漏效应的蒙卡模拟可以更真实地反映组件计算时的能谱状态。为验证此泄漏修正模型,一个由连续能量的TRIPOLI4模拟而得到的数值临界实验被用于分析与比较。通过与其他蒙卡程序SERPENT以及确定论程序ECCO进行结果对比,可证明此B1泄漏修正方法能够给出更精确的用于堆芯计算的少群截面参数。     

Energy exchange in structure scattering: a molecular dynamics study for Cs from Pt(1 1 1)

We have employed a classical molecular dynamics simulation to investigate the energy transfer of a heavy projectile ion to a surface, i.e. Cs⁺ impacting onto Pt(1 1 1), for incidence energies between 25 and 100 eV and an incidence angle of 45°. The in-plane scattering results show a continuous increase of the final energy with increasing scattering angle. All scattering intensities have a main supraspecular peak and scattering into subspecular angles increases with increasing incidence energy. The large projectile/target mass ratio causes a high energy loss and a strong angular dependence of the final energy distribution. The trends of the energy transfer and its angular dependence can be understood in terms of a binary collision model, augmented with double collisions and an the image charge correction. Backscattering at high incidence energies leads to a distribution of very low final energies, indicating the onset of surface sputtering. Peaks in the energy spectra arise from impact site dependent scattering and can be assigned to single, double, triple or sputtering type collisions.     

HIAF高能辐照终端感生放射性

本工作是基于蒙特卡罗模拟软件FLUKA对高能强流重离子加速器（HIAF）高能辐照终端感生放射性进行初步研究。该终端可运行质子最高能量为9.3 GeV,最大流强是1.45×1012 pps（particle per second）。研究内容包括:（1）预测高能辐照终端内活化物质的放射性活度特性;（2）预测不同冷却时间高能辐照终端内残余剂量率分布。研究结果表明,HIAF正常运行时高能辐照终端内的感生放射性主要受束流垃圾桶活化产生的放射性核素影响。当加速器连续运行100天冷却4小时,垃圾桶表面残余剂量率为2.375 mSv·h-1。终端内空气中13N和15O动态饱和比浓度大于其对应的导出空气浓度。冷却水中13N和15O的活度大于对应的ALI_min。该研究是HIAF辐射防护基础研究以及加速器环境影响评价的一项重要内容。The Monte Carlo code FLUKA was used to predict the induced radioactivity of high-energy irradiation terminal of HIAF. The maximum energy of proton is 9.3 GeV, and the maximum current is 1.45×1012 pps (particle per second). In this study we were to predict:(1) the activity properties of activated substances in the experimental terminal; (2) the residual dose rate distribution in the experimental terminal at different cooling time. The results indicate that the induced radioactivity in the high energy irradiation terminal of the HIAF is mainly affected by the radionuclide induced in the beam dump. The residual dose rate on the surface of the beam dump is 2.375 mSv·h-1, after 100 d irradiation and 4 h cooling. The dynamic saturation ratio of 13N and 15O induced in the air inside the terminal is higher than its corresponding derived air concentration. The activity of 13 N and15O induced in cooling water is higher than its ALI_min. This study is a part of radiation protection basic research and environmental impact assessment for HIAF.     

Particle in Cell/Monte Carlo Collision Method for Simulation of RF Glow Discharges: Effect of Super Particle Weighting

A parallel Particle in Cell/Monte Carlo Collision (PIC/MCC) numerical code for glow discharge plasma simulations is developed and verified. This method is based on simultaneous solution of the Lorentz equations of motion of super particles, coupled with the Poisson's equation for electric field. Collisions between the particles are modelled by the Monte Carlo method. Proper choice of particle weighting is critically important in order to perform adequate and efficient PIC simulations of plasma. Herein, effects of particle weighting on the simulations of capacitive radio‐frequency argon plasma discharges are studied in details. (© 2014 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

粗糙海面对高斯分布激光光束的反射模型推导

研究激光光束海面反射光强的方向分布特性对海上光电对抗等领域的工程实践具有重要意义.本文采用分形方法模拟粗糙海面,并在海面基准坐标系中建立起描述粗糙海面几何特征的数学方程,然后基于蒙特卡罗方法模拟高斯光束,依据几何光学原理在基准坐标系下推导了高斯光束的海面反射模型,采用该模型可以编程计算激光光束海面反射光强的方向分布.将模拟计算结果与实验结果进行了对比分析,结果表明该模型可以较好地反映激光光束海面反射光强的分布趋势,验证了模型的有效性.     

In vitro cell irradiation systems based on 210Po alpha source: construction and characterisation

One way of studying the risk to human health of low-level radiation exposure is to make biological experiments on living cell cultures.

Two ²¹⁰Po -particle emitting devices, with 0.5 and 100 MBq activity, were designed and constructed to perform such experiments irradiating monolayers of cells. Estimates of dose rate at the cell surface were obtained from measurements by a PIPS -particle spectrometer and from calculations by the SRIM 2000, Monte Carlo charged particle transport code. Particle fluence area distributions were measured by solid state nuclear track detectors. The design and dosimetric characterisation of the devices are discussed.     

Scattering and recoiling mapping of the Kr–Pt(1 1 1) system by SARIS

The technique of angle resolved mapping of scattering and recoiling imaging spectra (SARIS) combined with computer simulations is demonstrated to be a valuable tool for characterization of atomic collision events on surfaces. The energy distributions of scattered Kr and fast recoiled Pt atoms from a Pt(1 1 1) surface were measured as a function of exit angle. The use of a large area microchannel plate detector and time-of-flight techniques decreases the collection time and increases the number of detected trajectories above that of other designs. Classical ion trajectory simulations using the three-dimensional scattering and recoiling imaging code are used to simulate the kinematics of the scattering and recoiling particles. It is shown that SARIS mapping allows one to probe the kinematics of both scattered and recoiled particles, the probability for their occurrence in specific trajectories, their detection probabilities, and their threshold detection velocity. The measured and simulated energy distributions agree quantitatively if the detection efficiency is taken into account. The observed value of the threshold detection velocity for Pt atoms, ν^th=3.78(5)×10⁴ m/s, is in good agreement with previous studies.     

Monte Carlo simulations of oscillations, chaos and pattern formation in heterogeneous catalytic reactions

Experimental studies employing surface science methods indicate that kinetic oscillations, chaos, and pattern formation in heterogeneous catalytic reactions often result from the interplay of rapid chemical reaction steps and relatively slow complementary processes such as oxide formation or adsorbate-induced surface restructuring. In general, the latter processes should be analysed in terms of theory of phase transitions. Therefore, the conventional mean-field reaction–diffusion equations widely used to describe oscillations in homogeneous reactions are strictly speaking not applicable. Under such circumstances, application of the Monte Carlo method becomes almost inevitable. In this review, we discuss the advantages and limitations of employing this technique and show what can be achieved in this way. Attention is focused on Monte Carlo simulations of CO oxidation on (1 0 0) and (1 1 0) single-crystal Pt and polycrystal Pt, Pd and Ir surfaces and of NO reduction by CO and H₂ on Pt(1 0 0). CO oxidation on supported nanometre-sized catalyst particles and NO reduction on composite catalysts are also discussed. The results show that with current computer facilities the MC technique has become an effective tool for analysing temporal oscillations and pattern formation on the nanometre scale in catalytic reactions occurring on both single crystals and supported particles.     

蒙特卡罗与离散纵标耦合方法在压水堆堆腔漏束计算中的应用

为克服蒙特卡罗(MC)方法计算时间长和离散纵标(SN)方法复杂几何描述不精确的困难,采用SN-MC耦合计算流程,研究了基于蒙特卡罗方法和离散纵标方法的耦合方法。耦合方法的主要思想是根据离散纵标程序提供的中子角通量,利用接口程序计算出面源的分布概率,然后由修改后的源抽样子程序生成包括上下圆面源和圆柱面源的组合源,提供给蒙特卡罗程序进行计算分析。初步计算结果表明,该耦合方法是正确的,可用于压水堆堆腔漏束的计算分析。     

Observation of transient adsorption or skipping motion in small angle ion-surface scattering

Transient adsorption or skipping motion of silicon ions scattered under specular conditions from a Cu(111) surface is reported. The phenomena is identified on the basis of the observation of discrete energy loss peaks in the scattered positive and negative ion kinetic energy distributions, together with their behaviour upon variation of the crystal azimuth, incidence angle of the beam to the surface, and scattering angular distributions. Up to five reflections in a binding potential were seen. We propose that the beam is trapped following an energy loss, by electron-hole pair or plasmon excitation, which exceeds the incident beam energy component normal to the surface. The results lend strong support to the thesis that ion beam techniques can probe the family of potential surfaces upon which adsorption and reaction, at thermal energies, proceed.     

Energy transfer in hyperthermal Xe-graphite surface scattering

The scattering of a hyperthermal Xe from a graphite (0001) surface has been studied using a molecular beam-surface scattering technique and molecular dynamics (MD) simulations. The angular and velocity distributions of scattered Xe atoms were measured at incidence energies from 0.45 to 3.5 eV, three incidence angles of 15°, 35° and 60° and the surface temperatures of 300 K and 550 K. The observed time-of-flight spectra exhibit a sharp velocity distribution with only one velocity component, which is ascribed to the direct inelastic scattering process. The angle-resolved energy ratios of the mean final translational energy over the mean incidence energy E_f/E_i agree well with those predicted based on the assumption of the conservation of the momentum parallel to the surface. The Hard-Cube model, where the mass of the cube is approximately 310 u, has reproduced the angle-resolved flux distributions of scattered Xe atoms. In the Hard-Cube model almost 80% of the normal component of the incidence translational energy is found to be lost in collision. The MD simulations reproduce well the experimental results by using the Brenner potential for intralayer C atoms and a Lennard-Jones potential for interlayer C–C pair interactions.     

圆柱腔体光电输运的蒙特卡罗模拟

应用三维并行蒙特卡罗程序JMCT, 计算了特征温度分别为1, 3, 5, 8keV的黑体谱X射线入射到铝、二氧化硅、金的表面的背散射光电产额和电子能谱, 并与文献结果进行对比, 验证了程序的正确性, 进而针对系统电磁脉冲(SGEMP)研究中的典型几何结构——金属圆柱腔体, 模拟计算了其在黑体谱X射线照射下的光电输运过程, 用温度为1keV、注量为1J/m2的黑体谱X射线平行照射圆柱腔侧面, 半个侧面发射光电子, 计算得出了和方位角相关的不同面上发射光电子的光电产额、能谱分布和角分布, 结果表明掠入射的X射线会产生更高的光电产额; 光电子的发射角分布都基本符合余弦角分布的规律。     

X射线与CsI:Na转换屏作用规律的研究

本文用MonteCarlo方法模拟了轫致辐射X光子与CsI:Na转换屏的作用过程.给出了正向散射光子的角分布以及正向出射光子的能谱分布,指出正向散射光子的主要成分是CsI:Na的特征辐射光子,其能量为33.17和35.98keV.此外,还对吸收光子数与管电压和CsI:Na层厚度的关系进行了说明.