Pure adaptive search in global optimization

Pure adaptive seach iteratively constructs a sequence of interior points uniformly distributed within the corresponding sequence of nested improving regions of the feasible space. That is, at any iteration, the next point in the sequence is uniformly distributed over the region of feasible space containing all points that are strictly superior in value to the previous points in the sequence. The complexity of this algorithm is measured by the expected number of iterations required to achieve a given accuracy of solution. We show that for global mathematical programs satisfying the Lipschitz condition, its complexity increases at mostlinearly in the dimension of the problem.This work was supported in part by NATO grant 0119/89.     

硼中子俘获治疗的蒙特卡罗方法模拟

用通用蒙特卡罗程序MCNP模拟了粒子在人脑中的输运过程. 吸收剂量率主要来自以下四个反应:¹⁰B(n,α)⁷Li,¹⁴N(n,p)¹⁴C,¹H(n,γ)²D,快中子弹性散射反应.对肿瘤区的贡献主要来自硼中子吸收反应.结果表明,超热中子比热中子适合于深肿瘤的治疗,而热中子对浅肿瘤的治疗有优越性,比如皮肤癌.同确定论方法的结果相比,蒙特卡罗方法不失为一种模拟中子俘获治疗的好工具.     

Azine bridged silver coordination polymers: Powder X-ray diffraction route to crystal structure determination of silver benzotriazole

In continuation of our interest in solid-state structures of silver complexes of photographic importance, the structure for silver benzotriazole (AgBZT), has now been obtained. The preferred method for solving crystal structures is via single-crystal X-ray diffraction (XRD). However, for some materials, growing single crystals of appropriate size and quality is often difficult or even impossible. AgBZT is an example of such a silver complex with poor solubility. The usual routes to preparing single crystals using recrystallization from a cooperating solvent resulted in polycrystalline powder samples. We propose a crystal structure for AgBZT, solved from synchrotron X-ray powder diffraction data, using a direct-space Monte Carlo simulated annealing approach. AgBZT crystals are monoclinic, (P2₁/c), with unit cell dimensions, a=14.8052(3) Å, b=3.7498(4) Å, c=12.3495(12) Å, and β=114.200(6)°. The AgBZT complex is constructed from all three of the Benzotriazole (BZT) nitrogens bonding to a separate silver atom. As a consequence of this bonding mode, the structure is a highly cross-linked, coordination polymer.     

蒙特卡罗哈密顿新方案及其检验

蒙特卡罗哈密顿方法(MCH)是研究量子理论的数值模拟新方法, 其优点是可求出超出基态的能谱和波函数. 旧MCH方案需要自由粒子的信息, 较难推广应用于格点规范理论. 本文提出克服这个困难的新方案. 首先介绍这一方案的思想, 并以1维量子力学模型V(x)=μ²x²+λx⁴(其中μ²< 0,λ>0)为例说明实现这一新方案的具体计算步骤和方法.     

Global optimization by random perturbation of the gradient method with a fixed parameter

The paper deals with the global minimization of a differentiable cost function mapping a ball of a finite dimensional Euclidean space into an interval of real numbers. It is established that a suitable random perturbation of the gradient method with a fixed parameter generates a bounded minimizing sequence and leads to a global minimum: the perturbation avoids convergence to local minima. The stated results suggest an algorithm for the numerical approximation of global minima: experiments are performed for the problem of fitting a sum of exponentials to discrete data and to a nonlinear system involving about 5000 variables. The effect of the random perturbation is examined by comparison with the purely deterministic gradient method.     

Pervaporation of benzene/cyclohexane mixtures using ion-exchange membrane containing copper ions

The sorption, diffusion, and pervaporation (PV) properties of benzene/cyclohexane (Bz/Cx) mixtures on cation-exchange membranes containing copper ions (Cu(II)) were investigated. The equilibrium sorption isotherms of pure vapors in the membranes and the partial solubility of binary solutions in the membranes were described using the UNIQUAC model. The τ_iM and τ_Mi values were 0.978 and 0.591 for Bz, and 0.922 and 0.475 for Cx. The transient regimes of vapor sorption were employed to calculate the concentration-dependent diffusion coefficients. Long’s model sufficiently explained the diffusivity of Bz and Cx in the membranes. The pre-exponential factors were 3×10⁻¹³ m²/s and the plasticization factors were 3.0 and 3.6 for Bz and Cx, respectively. Excellent agreement was found with the experimental results applying the solubility and diffusivity data to simulate the pervaporation performance (flux and selectivity) using the modified Maxwell–Stefan equation. The membrane containing Cu(II) demonstrates better facilitating capability for Bz transport than that with Na(I), mainly due to its preferential sorption property toward Bz. Replacing Na(I) with Cu(II) into a Neosepta membrane resulted in better separation efficiency and higher Bz flux throughout the entire Bz concentration range.     

A kinetic Monte Carlo method on super-lattices for the study of the defect formation in the growth of close packed structures

A modified Kinetic Lattice Monte Carlo model has been developed to predict growth rate regimes and defect formation in the case of the homo-epitaxial growth of close packed crystalline structures. The model is an improvement over standard Monte Carlo algorithms, which usually retain fixed atom positions and bond partners indicative of perfect crystal lattices. Indeed, we extend the concepts of Monte Carlo growth simulations on super-lattices containing additional sites (defect sites) with respect to those of the reference material. This extension implies a reconsideration of the energetic mapping, which is extensively presented, and allows to describe a complex phenomenology that is out of accessibility of standard stochastic approaches. Results obtained using the Kawasaki and the Bond-Counting rules for the transition probability of the Monte Carlo event are discussed in details. These results demonstrate how the defect types (local or extended), the formation mechanisms and the defect generation regimes can be characterized using our approach.     

Electrode “edge effects” analyzed by the Green function method

An exact method based on Green's equation is used to find the diffusion-controlled faradaic current for certain electrode geometries that incorporate edges and vertices. Thereby the magnitudes of the time-independent current density associated with angled electrode/electrode and electrode/insulator junctions are calculated. As well, the square-root-of-time-dependent currents associated with vertices, receive attention. These terms extend to longer times, the Cottrell formulation appropriate for short times. Though most of the problems solved here have been tackled previously, the novel Green function approach is shown to be straightforward and intuitive.     

Rate of Convergence of a Stochastic Particle System for the Smoluchowski Coagulation Equation

By continuing the probabilistic approach of Deaconu et al. (2001), we derive a stochastic particle approximation for the Smoluchowski coagulation equations. A convergence result for this model is obtained. Under quite stringent hypothesis we obtain a central limit theorem associated with our convergence. In spite of these restrictive technical assumptions, the rate of convergence result is interesting because it is the first obtained in this direction and seems to hold numerically under weaker hypothesis. This result answers a question closely connected to the Open Problem 16 formulated by Aldous (1999).     

Monte Carlo Simulation of Electron Transport and X-Ray Generation. II. Radiative Processes and Examples in Electron Probe Microanalysis

Theoretical models for Monte Carlo simulation of radiative processes, i.e. bremsstrahlung and characteristic x-ray emission, are presented. Possible strategies for simulating electron transport are briefly described. For mechanisms involving energy loss and angular deflections, difficulties for strict implementation of accurate numerical differential cross sections still remain due to the strong correlations between these variables. Practical solutions for the case of inelastic collisions and bremsstrahlung emission are described. Comparisons of simulation results with experimental data for several problems of interest in electron probe microanalysis are presented.