Numerical methods for finding multiple eigenvalues of matrices depending on parameters

Summary. Let be a square matrix dependent on parameters and , of which we choose as the eigenvalue parameter. Many computational problems are equivalent to finding a point such that has a multiple eigenvalue at . An incomplete decomposition of a matrix dependent on several parameters is proposed. Based on the developed theory two new algorithms are presented for computing multiple eigenvalues of with geometric multiplicity . A third algorithm is designed for the computation of multiple eigenvalues with geometric multiplicity but which also appears to have local quadratic convergence to semi-simple eigenvalues. Convergence analyses of these methods are given. Several numerical examples are presented which illustrate the behaviour and applications of our methods. Received December 19, 1994 / Revised version received January 18, 1996     

A posteriori error estimation with finite element methods of lines for one-dimensional parabolic systems

Summary Consider the solution of one-dimensional linear initial-boundary value problems by a finite element method of lines using a piecewiseP ^th -degree polynomial basis. A posteriori estimates of the discretization error are obtained as the solutions of either local parabolic or local elliptic finite element problems using piecewise polynomial corrections of degreep+1 that vanish at element ends. Error estimates computed in this manner are shown to converge in energy under mesh refinement to the exact finite element discretization error. Computational results indicate that the error estimates are robust over a wide range of mesh spacings and polynomial degrees and are, furthermore, applicable in situations that are not supported by the analysis.This research was partially supported by the U.S. Air Force Office of Scientific Research, Air Force Systems Command, USAF, under Grant Number AFOSR 90-0194; by the U.S. Army Research Office under Contract Number DAAL03-91-G-0215; and by the National Science Foundation under Institutional Infrastructure Grant Number CDA-8805910     

On the Moments of Overflow and Freed Carried Traffic for the GI/M/C/0 System

In circuit-switched networks call streams are characterized by their mean and peakedness (two-moment method). The GI/M/C/0 system is used to model a single link, where the GI-stream is determined by fitting moments appropriately. For the moments of the overflow traffic of a GI/M/C/0 system there are efficient numerical algorithms available. However, for the moments of the freed carried traffic, defined as the moments of a virtual link of infinite capacity to which the process of calls accepted by the link (carried arrival process) is virtually directed and where the virtual calls get fresh exponential i.i.d. holding times, only complex numerical algorithms are available. This is the reason why the concept of the freed carried traffic is not used. The main result of this paper is a numerically stable and efficient algorithm for computing the moments of freed carried traffic, in particular an explicit formula for its peakedness. This result offers a unified handling of both overflow and carried traffics in networks. Furthermore, some refined characteristics for the overflow and freed carried streams are derived.     

Surface grating formation and erasing on a-Se films

Surface relief formation at holographic recording on amorphous selenium films was demonstrated and investigated. The presence of this optical phase modulation component is essential for ensuring significant, stable and erasable optical recording in a-Se films at 290–320 K temperatures, where conventional photodarkening was known as insignificant and unstable. Photocrystallization can only be observed in super-exposed a-Se films at the given experimental conditions of hologram recording. Erasing behavior of surface relief gratings under heat treatment was also investigated in order to reveal further details of the mechanism. Photoinduced structural transformations within the amorphous phase, connected to local ordering under the condition of light-induced fluidity, are proposed as an explanation for the relief formation and erasing. The observed reversible optical recording process may be useful for the various optoelectronic applications of photoconductive a-Se layers. Received: 12 June 2000 / Accepted: 6 June 2001 / Published online: 30 August 2001     

An algebraic characterization ofB-convergent Runge-Kutta methods

Summary In the analysis of discretization methods for stiff intial value problems, stability questions have received most part of the attention in the past.B-stability and the equivalent criterion algebraic stability are well known concepts for Runge-Kutta methods applied to dissipative problems. However, for the derivation ofB-convergence results — error bounds which are not affected by stiffness — it is not sufficient in many cases to requireB-stability alone. In this paper, necessary and sufficient conditions forB-convergence are determined.This paper was written while J. Schneid was visiting the Centre for Mathematics and Computer Science with an Erwin-Schrödinger stipend from the Fonds zur Förderung der wissenschaftlichen Forschung     

Point approximation of a space-homogeneous transport equation

Summary We present an approximation method of a space-homogeneous transport equation which we prove is convergent. The method is very promising for numerical computation. Comparison of a numerical computation with an exact solution is given for the Master equation.     

C72,La2@C72几何结构和电子性质的计算研究

采用密度泛涵理论(density functional theory,DFT)中的广义梯度近似(generally gradient approximation,GGA)对富勒烯C72和内掺金属La富勒烯La2@C72三种同分异构体的几何结构和电子结构进行研究.发现在C72的三种同分异构体中,满足独立五边形规则(isolated-pentagon-rule,IPR)的C72(D6d)结构最为稳定;在La2@C72三种同分异构体中,有着两对两两相邻五元环(twofused-pentagon)的La2@C72(#10611)结构最为稳定,而满足IPR的La2@C72(D6d)的结构变成了最不稳定结构.从能级图和态密度图分析得知,笼子稳定性的变化与La原子的原子轨道与C72原子轨道之间的杂化有关.Mulliken电荷分析得知,La2@C72(#10611)的两个La原子共转移了约3个电子给C72,并且,它们几乎分布在整个C笼上,形成的电子结构为La23+@C3-72.净自旋分析得知,La2@C72(#10611)中La原子磁性完全淬灭.     

Ag+掺杂的立方相Y2O3:Eu纳米晶体粉末发光强度研究

采用化学自燃烧法制备了不同Ag⁺掺杂浓度的Y₂O₃:Eu纳米晶体粉末样品(［Y³⁺］∶［Eu³⁺］∶［Ag⁺］=99∶1∶X，X=0—3.5×10^-2)，以及通过退火处理得到了相应的体材料.根据X射线衍射谱确定所得纳米和体材料样品均为纯立方相.实验表明在纳米尺寸样品中随着Ag离子浓度的增加，荧光发射强度随之增加，当X=2×10^-2时达到最大值，其发光强度比X=0时提高了近50%.当Ag离子浓度继续增加，样品发光强度保持不变.在相应的体材料样品中则没有观察到此现象.通过对各样品的发射光谱，激发光谱，X射线衍射图谱，透射电镜(TEM)照片和荧光衰减曲线的研究，分析了引起纳米样品荧光强度变化的原因是由于Ag离子与表面悬键氧结合，从而使这一无辐射通道阻断，使发光中心Eu³⁺的量子效率提高；Ag⁺的引入所带来的另一个效应是使激发更为有效.这两方面原因使发光效率得到了提高.