均匀磁场下碳离子笔形束的剂量变化分析及位置修正方法

分析外加均匀磁场对于碳离子笔形束剂量分布的影响,并考虑修正这种影响,为磁共振成像引导碳离子放射治疗的临床应用提供指导。本文利用蒙特卡罗方法模拟计算了不同能量碳离子笔形束在不同强度磁场下的剂量分布情况,发现垂直于碳离子束入射方向的均匀磁场对于碳离子笔形束射程缩短的影响很小,磁场对碳离子束的主要影响是引起束流横向偏转,特别是碳离子束布拉格峰位置的横向侧移。横向侧移程度与碳离子束的能量和磁场强度相关,根据模拟结果,得到了一个计算碳离子束布拉格峰在磁场中相对横向偏转的方程,并提出一种校正外加磁场引起的碳离子束布拉格峰横移的角度修正方法。这些结果可用于指导磁共振图像引导碳离子放射治疗计划系统的研发。     

80.5 MeV/u碳离子诱发铜靶的放射性剩余产物测量

高能重带电粒子能直接穿透靶原子核外电子层,与原子核发生直接碰撞,发生散裂反应,产生一系列具有放射性的剩余产物核.重带电粒子诱发靶材放射性剩余核与辐射防护和人员安全有着密切联系,当前,大部分剩余核产额主要依靠蒙特卡罗粒子输运程序进行模拟计算,其准确程度亟需通过实验测量进行准确评估.本文利用能量为80.5 MeV/u的(12)^C⁶⁺粒子对薄铜靶开展了辐照实验与伽玛射线测量,结合伽玛谱学分析方法,得出了辐照产生的18种放射性剩余产物的初始活度和产生截面值,并与PHITS模拟结果进行对比.结果表明,PHITS模拟程序对放射性剩余核种类的估计具有较高可靠性,在其绝对产额方面,与实验测量仍具有较大偏差.     

Monte Carlo simulation of electric conductivity for pure and doping fullerene (C60)

In this paper, we have used Monte Carlo (MC) method to simulate and study the temperature and doping effects on the electric conductivity of fullerene (C60). The results show that the band gap has reduced by the doping and the charge carrier transport is facilitated from valence band to conduction band by the temperature where is touched a 300 K. In this case, the conductivity reached a value of $4\times {10}^{-7}\text{S}{\text{cm}}^{-1}$. The electric conductivity of C60 can increase by the triphenylmethane dye crystal violet (CV) alkali metal to reach $4\times {10}^{-3}\text{S}{\text{cm}}^{-1}$ at 303 K. Our results of MC simulation have a good agreement with those extracted from literature [10], [33].     

Carathéodory's method for a class of dynamic games

The method of equivalent variational methods, originally due to Carathéodory for free problems in the calculus of variations is extended to investigate feedback Nash equilibria for a class of n-person differential games. Both the finite-horizon and infinite-horizon cases are considered. Examples are given to illustrate the presented results.     

Monte Carlo study of interfacial silicon suboxide layers and oxidation kinetics

A simple simulation scheme that simultaneously describes the growth kinetics of SiO₂ films at the nanometer scale and the SiO_x/Si interface dynamics (its extent, and spatial/temporal evolution) is presented. The simulation successfully applies to experimental data in the region above and below 10 nm, reproduces the Deal and Grove linear-parabolic law and the oxide growth rate enhancement in the very thin film regime (the so-called anomalous region). According to the simulation, the oxidation is governed mainly by two processes: (a) the formation of a transition suboxide layer and (b) its subsequent drift towards the silicon bulk. We found that it is the superposition of these two processes that produces the crossover from the anomalous oxidation region behavior to the linear-parabolic law.     

An optical measurement of vortex shape at a free surface

We have proposed an optical method of vortex shape measurement based on Fourier transform profilometry (FTP) and verified it by experiment. The results of our experiment proposed in this paper show that FTP can efficiently reconstruct the vortex shape at a free surface and this method is suitable for wide use in studying such problems as liquid shear flow, wake of an object, flow behind a bluff body, and wetting angle.     

锥形光纤中光波传输特性的蒙特卡罗模拟

为了分析锥形光纤中光子流的传输特性和空间分布特性，采用了蒙特卡罗算法对其进行了模拟仿真，得到了大量光子经过锥形光纤后的光子分布图.模拟实验结果表明，通过锥形光纤的大量光子集中在探测器中心，到达探测器某一位置的光子数随着该位置到探测器中心距离的增加而减少，并且在锥形光纤锥度比一定的条件下，改变锥形光纤长度，光强度随着长度的增加而下降，但通过锥形光纤的光子分布是一样的.这些研究结果对锥形光纤的设计和实际应用提供了重要依据.     

Non-universal critical behaviour of heteronuclear dimers on a square lattice––A Monte Carlo study

Monte Carlo simulation within the grand canonical ensemble, the histogram reweighting technique, and finite size scaling analysis are used to explore the phase behaviour of heteronuclear dimers, composed of A and B type atoms, on a square lattice. We have found that for the models with attractive BB and AB nearest-neighbour energy, u_BB=u_AB=−1, and for non-repulsive energy between AA nearest-neighbour sites, u_AA<0, the system belongs to the universality class of the two-dimensional Ising model. However, when u_AA>0, the system exhibits a non-universal critical behaviour. We have evaluated the dependences of the critical point characteristics on the value of u_AA.     

电子回旋共振放电的电离特性PIC/MCC模拟(Ⅰ)——物理模型与理论方法

采用粒子模拟与蒙特卡罗相结合(PIC/MCC)的方法对电子回旋共振(ECR)放电中的电离过程进行了模拟，其中带电粒子与微波的相互作用由PIC方法的电磁模型描述，粒子间的碰撞过程由MCC方法描述.考虑的碰撞类型有电子与中性粒子的弹性、激发、电离碰撞，离子与中性粒子的弹性、电荷交换碰撞，碰撞截面均依赖于能量而变化.阐述了理论分析的过程，为数值模拟ECR放电奠定了基础. 关键词： 电子回旋共振放电 粒子模拟 蒙特卡罗 电离