Monte Carlo model of the temperature rise at a GaAs surface under an electron beam

Scanning electron microscopy (SEM) has frequently been used to study semiconductor materials. It offers the possibility of obtaining reliable qualitative and quantitative information on relevant local material parameters. The temperature rise due to electron‐beam bombardment can influence some semiconductor parameters, which then will influence the SEM information. In this work we propose a model calculation based on the Monte Carlo (MC) method to calculate the temperature rise due to electron‐beam heating. The results show that the temperature rise increases with increasing numbers of electrons (electron‐beam current), and the inverse behavior is observed with respect to the electron energy (electron‐beam voltage). The decrease in temperature rise with depth is also obtained. Copyright © 2006 John Wiley & Sons, Ltd.     

The features of electron dose distributions in circular objects: Comparison of Monte Carlo calculation predictions with dosimetry

The features of electron dose field formation in the multi-layer circular objects are related with its surface irregularity such as convexity, concavity and roundness of inner and outer layers. The simulation of dose distributions in multi-layer tubes irradiated with a scanned electron beam (EB) was carried out by Monte Carlo (MC) method with utilization of the software ModeCEB. The effects of mutual influence on dose field formation in contacting multi-layers tubes irradiated with EB were MC simulated and measured with a film dosimetry. An experimental validation of the obtained simulation predictions for dose distributions in multi-layer tubes irradiated with 10 MeV electrons was performed on radiation facility with linear electron accelerator LAE 13/9, INCT, Warsaw. Comparison of MC simulation results with a film dosimetry is discussed in the report.     

Quantum Monte Carlo Calculation of Correlation Effects on Bond Orders

We present a computational approach, using quantum Monte Carlo, that provides some insight into the effect of electron correlation on chemical bonding between individual pairs of atoms. Our approach rests upon a recently suggested relation between the bond order and charge fluctuations with respect to atomic domains. Within the present implementation we have taken a compromise between conceptual rigour and computational simplicity. In a first step atomic domains were obtained from Hartree-Fock (HF) densities, using Bader’s definition of atoms in molecules. These domains were used in a second step in quantum Monte Carlo calculations to determine bond orders for pairs of atoms. Correlation effects have been studied by comparison of HF bond orders with those obtained from pure diffusion quantum Monte Carlo calculations. We illustrate this concept for C–O and C–S bonds in different molecular environments. Our results suggest an approximate linear relation between bond order and bond length for these kinds of bonds.     

Monte Carlo analysis of the electron thermalization process in the afterglow of a microsecond dc pulsed glow discharge

A Monte Carlo model is utilized for studying the behavior of electrons in the afterglow of an analytical microsecond dc pulsed glow discharge. This model uses several quantities as input data, such as electric field and potential, ion flux at the cathode, the fast argon ion and atom impact ionization rates, slow electron density, the electrical characterization of the pulse (voltage and current profiles) and temperature profile. These quantities were obtained by earlier Monte Carlo — fluid calculations for a pulsed discharge. Our goal is to study the behavior of the so-called Monte Carlo electrons (i.e., those electrons created at the cathode or by ionization collisions in the plasma which are followed by using the Monte Carlo model) from their origin to the moment when they are absorbed at the cell walls or when they have lost their energy by collisions (being transferred to the group of slow electrons) in the afterglow of the pulsed discharge. The thermalization of the electrons is a phenomenon where the electron-electron Coulomb collisions acquire a special importance. Indeed, in the afterglow the cross sections of the other electron reactions taken into account in the model are very low, because of the very low electron energy. We study the electron energy distributions at several times during and after the pulse and at several positions in the plasma cell, focusing on the thermalization and on the behavior of the electrons in the afterglow. Also, the time evolution of the rates of the various collision processes, the average electron energy, the densities of Monte Carlo and slow electrons and the ionization degree are investigated.     

A simple approach to treat anisotropic elastic collisions in Monte Carlo calculations of the electron energy distribution function in cold plasmas

A simplified method allowing one to treat anisotropic electron heavy species elastic scattering in MonteCarlo models of gas discharges with the proper value for collision frequency is proposed The method is applied to an electric discharge in a Ne · Xe/HCl mixture, and the results are compared with the solution of the two-term expansion of the Boltzmann equation under the same conditions. Methods for reduction of computational time in Monte Carlo codes and the use of the Monte Carlo flux method are also discussed.     

Monte Carlo study of electron correlation functions for small molecules

A study is made of electron-electron correlation functions for use in trial wave functions for small molecules. New forms are proposed that have only a few variational parameters, and these parameters have physical meanings that are easily discerned. Total energies for H₂, LiH and Li₂ computed using these correlation functions are presented, and comparison is made with previous forms, including the Jastrow-Pade form often used in Monte Carlo studies. We further treat the possibility that correlation depends not only on the separation of a pair of electrons but also on the location of the electron pair relative to the nuclei — indicative of a density-dependent or many body correlation effect. Our results indicate that such a many-body correlation effect is weakly present.This work was supported by the Director, Office of Energy Research, Office of Basic Energy Sciences, Chemical Sciences Division of the U.S. Department of Energy under Contract No. DE-AC03-76SF00098     

Monte Carlo方法在高分子科学中的应用

介绍了Monte Carlo方法的历史及其特点,并描述了它在现代高分子科学研究中的广泛应用情况,并对其前景作了一些展望。     

表面活性剂与高分子链混合体系的模拟

计算机模拟了高分子链对表面活性剂胶束形成过程的影响，以及高分子链构象性质随胶束化过程的变化.结果表明,当高分子链与表面活性剂之间的相互作用强度超过临界值后，高分子链的存在有利于表面活性剂胶束的形成.临界聚集浓度（CAC）与临界胶束浓度（CMC）的比值CAC/CMC随高分子链长的增大和相互吸引作用的增强而减小.在CAC之前，高分子链与表面活性剂分子只有动态的聚集；但在CAC之后，表面活性剂胶束随表面活性剂浓度X的增加而增大，并静态地吸附在高分子链上，形成表面活性剂/高分子聚集体.随着表面活性剂分子的加入，高分子链的均方末端距和平均非球形因子先保持恒定；从X略小于CAC开始， 和快速减小，至极小值后又逐渐增大.模拟结果支持高分子链包裹在胶束表面的实验模型.     

Monte Carlo Simulations of Polymers in Nanoslits

Summary: Monte Carlo computer simulations have been performed for model polymers confined in slits of thickness comparable to the transverse diameter of the chains. The density of polymer within the slits is allowed to vary with the slit thickness in such a way that the content of the slits is always in equilibrium with a large reservoir of bulk polymer. The calculations reveal the presence of polymer‐mediated attractive or repulsive interactions between the slit plates, oscillating with the slit thickness in good agreement with experimental results.

The base cell used in the simulations.     

Monte Carlo calculations of hole size distributions: Simulation of positron annihilation spectroscopy

Consequences are explored of a hole size distribution in an amorphous polymer for the ortho-positronium (o-Ps) lifetime (τ₃) and intensity (I₃), determined by positron annihilation lifetime spectroscopy. The disordered lattice model, with a vacancy fraction h as a central quantity, is used to represent the equation-of-state behavior of the polymer. By means of Monte Carlo simulations, we obtain the cluster size distribution as a function of h and hence temperature. The predicted average cluster size and the cluster concentration are compared to τ₃ and I₃ data, respectively, for bisphenol-a polycarbonate. Furthermore, the influence of an o-Ps lifetime distribution on the experimental mean τ₃ is investigated. By mimicking the computational methods used in experimental analysis, agreement between experiment and theory in respect to τ₃ and to I₃ in the melt ensues. In the glass, however, the experimental I₃ becomes increasingly smaller with decreasing temperature than is computed. These deviations may result from a distortion of the equilibrium free volume. © 1992 John Wiley & Sons, Inc.     

Optimized Monte Carlo sampling in forward-backward semiclassical dynamics

Forward-backward semiclassical dynamics (FBSD) provides a rigorous and powerful methodology for calculating time correlation functions in condensed phase systems characterized by substantial quantum mechanical effects associated with zero-point motion, quantum dispersion, or identical particle exchange symmetries. The efficiency of these simulations arises from the use of classical trajectories to capture all dynamical information. However, full quantization of the density operator makes these calculations rather expensive compared to fully classical molecular dynamics simulations. This article discusses the convergence properties of various correlation functions and introduces an optimal Monte Carlo sampling scheme that leads to a significant reduction of statistical error. A simple and efficient procedure for normalizing the FBSD results is also discussed. Illustrative examples on model systems are presented.     

Measurements and calculations of absorbed dose in electron beam radiation processing

An electron beam current densimeter has been described and used for dose measurement in EB radiation processing. An improved bipartition model of electron transport is applied to calculate the reflection coefficients and energy deposition distributions produced by 0.2–3 MeV electrons in semi-infinite media of Al, C, MAR and nylon, and the energy deposition produced by 2 MeV electrons in the two-layer medium consisting of copper and polystyrene. In addition, the depth dose distribution of 300 keV electrons in Ti-air-nylon composite-layer has been evaluated. The calculation results are in good agreement with the measurement data.     

Monte Carlo global sensitivity analysis of Boltzmann equation for electron kinetics in H2 discharges

Global sensitivity analysis with the Monte Carlo method is applied to the Boltzmann equation for the electron energy distribution function (eedf). The results show the sensitivity of eedf and related quantities on the global variation of cross sections set. A new indicator of global sensitivity is used, which appears to be of easy applicability.     

端点附壁的高分子链形状的Monte Carlo模拟

分别基于简立方格点和四面体格点模型对一端吸附在无限大平面的高分子链（平面接枝高分子链）的形状进行了Ｍｏｎｔｅ Ｃａｒｌｏ模拟,结果表明,接枝高分子链的形状更偏离球形,〈Ｌ＾２１〉：〈Ｌ＾２２〉：〈Ｌ＾２３〉的极限值约为１：２．７５：１２．５,其中〈Ｌ＾２１〉,〈Ｌ＾２２〉和〈Ｌ＾２３〉分别为回转半径张量的本征值Ｌ＾２１,Ｌ＾２２和Ｌ＾２３（Ｌ＾２１〈Ｌ＾２２〈Ｌ＾２３）的统计平均;链长相同时,接枝     

分形表面上CO氧化反应的蒙特卡罗模拟

众所周知,人们在理想完好的单晶表面上的动力学研究已取得了很大的成功[1],然而对许多实际催化剂来说,传统催化理论和实验之间不符的情况仍然俯拾皆是．现在人们对造成这种情况的原因已经有了更加深入的认识：大量的实验事实表明,在表面催化反应体系中,不但催化剂表面具有复杂的分形结构,而且催化剂表面上的活性中心分布也具有复杂的分形特性[2]．“分形”是指那些具有分数维数的几何对象．这些对象往往是不规则的,不能用通常的欧氏几何来描述[3]．将分形引入催化科学中最早的是Pfeifer和Avnir等[4]人,至今已有近二十年的历史,现…     

Monte Carlo simulations of biomolecules: The MC module in CHARMM

We describe the implementation of a general and flexible Monte Carlo (MC) module for the program CHARMM, which is used widely for modeling biomolecular systems with empirical energy functions. Construction and use of an almost arbitrary move set with only a few commands is made possible by providing several predefined types of moves that can be combined. Sampling can be enhanced by noncanonical acceptance criteria, automatic optimization of step sizes, and energy minimization. A systematic procedure for improving MC move sets is introduced and applied to simulations of two peptides. The resulting move sets allow MC to sample the configuration spaces of these systems much more rapidly than Langevin dynamics. The rate of convergence of the difference in free energy between ethane and methanol in explicit solvent is also examined, and comparable performances are observed for MC and the Nosé-Hoover algorithm. Its ease of use combined with its sampling efficiency make the MC module in CHARMM an attractive alternative for exploring the behavior of biomolecular systems.     

Diffusion Monte Carlo study of correlation in the hydrogen molecule

The diffusion Monte Carlo (DMC) method shows that correlation in H₂ produces a set of three spatial changes: (i) an enhancement in the electron density distribution n( r ) in the left and right anti‐binding regions that include separately the immediate vicinity of each of the two nuclei, (ii) a reduction in n( r ) in the binding region intervening between the two nuclei as a counterbalance, and (iii) a concomitant increase in the equilibrium internuclear separation. It is stressed that the correlation energy E^c (= T^c + V^c) for diatomic molecules be defined by the difference in the total energy between the exact and the Hartree–Fock (HF) variational calculations that are performed at individually optimized internuclear separations. It is this definition that makes it possible to involve a significant contribution from a correlation‐induced change in the equilibrium internuclear separation as part of the correlation energy and to relate (i) and (ii) to (iii) in consistency with the electrostatic theorem. The present calculations fulfill the virial theorem to an accuracy of ?V/T = 2.00 for DMC and ?V^HF/T^HF = 2.000 for HF. The present correlation energy E^c = ?0.0408 hartree is not only in good agreement with the most accurate value previously reported, but also can be analyzed into all its components in accordance with the correlational virial theorem 2T^c + V^c = 0. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Monte Carlo Simulation of Electron Transport and X-Ray Generation. I. Electron Elastic and Inelastic Scattering

We present different theoretical approaches to determine differential cross sections for elastic and inelastic interactions of electrons. These cross sections are the basic ingredients for accurate Monte Carlo simulation of electron transport in matter. The considered models range from simple analytical approximations employed in early calculations to purely numerical differential cross sections described by large databases calculated with state-of-the-art theory.     

A sparse algorithm for the evaluation of the local energy in quantum Monte Carlo

A new algorithm is presented for the sparse representation and evaluation of Slater determinants in the quantum Monte Carlo (QMC) method. The approach, combined with the use of localized orbitals in a Slater-type orbital basis set, significantly extends the size molecule that can be treated with the QMC method. Application of the algorithm to systems containing up to 390 electrons confirms that the cost of evaluating the Slater determinant scales linearly with system size.     

Monte Carlo simulation on layered polymeric films

Thin films of polymer blends composed of alternating copolymer, diblock copolymer and/or homopolymer are studied using Monte Carlo simulation. A multilayer morphology is observed in the film, that is, the blended polymers assemble into individual domains arranged from interior to the surfaces of the film. The coexisting components residing throughout the neighboring domains in the film make no distinguishable interface between any neighboring domains. By this means, it forms a vertical composition gradient in the polymeric film. Being different from layer-by-layer deposition of polyelectrolyte or hydrogen bonding approach etc., the layered structure in this study is formed by polymer blending in one step. Alternating copolymers are found to be essential components to form vertical composition gradient （layered structure） in thin films.