Early stages of silicon oxidation

Experiments have shown that the early stages of silicon oxidation proceed layer by layer, so that one layer is essentially complete before another develops. Other experiments show that the mechanism does not involve step growth, the most obvious mechanism. We use a new approach to modelling the growth to show that these two observations can be understood when there is a rate-determining step which depends strongly on the local oxide thickness. The rate in question might be the sticking probability, or the rate of incorporation of adsorbed oxygen species into the oxide network. Such mechanisms are possible when transport by an ionic species dominates, contrary to the situation for thicker films. Our modelling suggests the mechanisms are driven by the image interaction, as in earlier suggestions by Stoneham and Tasker, rather than an effect of the electric field central to the Mott-Cabrera mechanism.     

SU(2)纯格点规范场热力学非对称点阵Monte-Carlo模拟

在非对称点阵(ξ=2,3)上模拟了SU(2)纯格点规范场系统的热力学行为,结果表明,非对称点阵Monte-Carlo模拟是实现快速趋向连续物理极限的有效途径.     

The radiative transfer in cylindrical medium and partition allocation method by overlap regions

In this paper the radiative transfer in one-dimensional and two-dimensional cylindrical medium is simulated by the Monte-Carlo (M-C) method. Our results agree with the previous ones very well. It indicates that our cylindrical M-C model is creditable. In this paper, the partition allocation method of radiative heat transfer in participating cylindrical media is presented, in which every sub-domain is isolated by imaginary black wall at certain equivalent temperature and overlaps each other. The stitched results of all sub-domains can predict the results of a whole zone accurately. The partitioned treatment by overlap regions can achieve reasonable result, save memory efficiently, and compute parallel.     

Weibull分布损伤失效率模型常应力下的参数估计

Bhattacharyya和Soejoeti(1989)对步进应力加速寿命试验提出损伤失效率模型(TFR模型)．本文在全加速步加试验场合下指出文献[1]，[2]，[5]中用通常的回归分析方法求取常应力下参数的估计是不合理的，同时给出了如何求取常应力下参数估计的一种方法．     

X-Ray Emission from Thin Films on a Substrate – Calculation and Experiments

 Monte Carlo simulation of electron transport in solids is widely used in electron microscopy, spectroscopy and microanalysis. The reliability of physical models incorporated in a Monte Carlo code is usually checked by comparing with experimental results. Elastic or inelastic collisions are usually considered as the basic interactions of electrons with atoms. In our Monte Carlo code the single scattering model is employed for simulation of X-ray emission from thin films of Au on the Si substrate. The electron beam energy was in the range 10–30 keV, the take-off-angle was 40°. The simulated values of X-ray production were calculated in our Monte Carlo code using several models of ionisation cross-sections. For the emitted intensities the depths of inelastic collision and X-ray absorption were taken into account, then the k-ratios were calculated. These data were compared with experimental values of k-ratios calculated from X-ray intensities of Au M and Au L characteristic lines. We followed mainly the dependence of the k-ratios of the film element on film thickness. The film thickness was in the range 0.05–1 μm. Reasonably good agreement was found for dependences of X-ray intensity on film thickness in the whole energy range and for both lines, especially for Powell’s model of the ionisation cross-section.     

Computer models for electrical conductivity of multicomponent solid-state systems

Models for the computer simulation of electrical conductivity of polycrystalline materials and mixtures are proposed. The effect of model parameters on the concentration behavior of conductivity is considered. The correspondence between the model parameters and the characteristics of physical objects is established, thus making it possible to analyze the spatial relations in solid-state systems.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 11, pp. 1868–1873, November, 1993.     

二维指数信号模型中参数Bootstrap估计的强相合性

本文主要研究了在统计信号处理当中具有广泛应用的二维带白噪声指数信号模型中参数估计的Bootstrap逼近, 借助于回归模型中Bootstrap逼近的构造方法, 给出了二维指数信号模型参数的自助估计, 并证明了自助估计具有强相合性. 最后采用了Monte-Carlo法对所提的方法进行随机模拟, 模拟的结果表明当噪声不服从正态分布时, Bootstrap方法的估计效果优于最小二乘估计.     

Particles movement and surface quality in PLD/PR systems

A three-dimensional model based on Monte-Carlo and Finite Elements techniques has been used for simulating plume behavior, ‘micron-sized particles’ movement and interaction with obstacles in a Pulsed Laser Deposition with Plasma Reflection (PLD/PR) system. Have been investigated the influences of mass, surface size and emission time on trajectory and film surface quality as well. Droplet and ‘big-size particles’ deposition statistics are presented and a comparison between theoretical and experimental results upon thin film surface quality as well. One can observe that particles mass and surface size have a strong influence on the particles trajectory by affecting the collisions parameters during the entire propagation process. The emission time should influence the particles trajectory by affecting the probability of interaction with other particles. By making a 10,000 particles statistic for a normal distribution of these investigated parameters, we obtain reasonable good results in modeling ‘big-size particles’ tendency to be deposited at lower reflection angles. These results sustain assumption of ‘big particles’ deflection by plume fine particles during the propagation process.     

Monte-Carlo法模拟二维Ising模型 ——Metropolis、Swendsen-Wang与Wolff算法的对比

Ising模型是一种应用广泛的磁自旋相互作用模型,其二维情况严格求解极为复杂,实际应用中通常利用Wolff算法进行模拟.Wolff算法目前被认为是最好的聚类翻转Monte-Carlo算法.Metropolis和Swendsen-Wang算法同Wolff算法类似,理论上也适用于Ising模型的模拟,却未有文章将三者系统对比来说明Wolff算法的优越性,本科课程对于Monte-Carlo算法的介绍也较少.本文分别利用三种算法模拟了二维Ising模型,介绍了其算法原理、参数选择及实现方式,分析对比了三种算法的模拟效果和适用范围,从而总结说明在二维Ising模型的模拟中Wolff算法效果更好的原因.