Differential algebraic method for computing the high order aberrations of practical electron lenses

Differential algebraic method is of an effective technique in computer numerical analysis. It implements exactly differentiation up to arbitrary high order based on the nonstandard analysis. Some complicated nonlinear dynamics problems including high order aberrations of electron optics systems can be solved by mapping properties of differential algebraic quantities. However, the existing electron optical differential algebraic method can only be applied to those problems where the electric and/or magnetic fields are expressed in analytic forms. In this paper, the principle of differential algebraic method is applied to practical electron lenses whose electric/magnetic fields are in the forms of discrete arrays, for example, the files computed by FEM or FDM method. Thus the developed new differential algebraic method is applicable to engineering design. The programs were written for computing the high order aberrations of practical electron lenses. An example is given to show that the numerical results have good accuracy compared with the results computed by using the electric fields expressed in analytical form; the accuracy is limited only by the accuracy of the numerical computation of the fields and the arithmetic errors, and it is enough for engineering design.     

Differential algebraic method for arbitrary-order chromatic aberration analysis of electrostatic lenses

In this paper, differential algebra is applied to calculate arbitrary high-order chromatic aberrations of electrostatic electron lenses. Expressions of differential algebraic form of high-order combined chromatic aberration coefficients are obtained and arbitrary-order chromatic aberrations can be calculated numerically. As an example, a typical Schiske's electrostatic lens has been studied. All the first- to third-order chromatic aberration coefficients of the lens have been calculated, and the pattern of the first-order chromatic aberration has been given as well.     

Aberration analysis of electron mirrors by a differential algebraic method

A differential algebraic (DA) method has been developed for the aberration analysis of electron mirrors. Since large ray slopes occur near the turning points in mirrors, the axial position is no longer suitable as the independent variable and the electron trajectory equation used in conventional lens theory is no longer feasible. A DA solution of the electron motion equation, wherein a single DA ray trace is performed on a non-standard extension of real number space called _nD_v, enables the aberrations of a mirror system to be obtained, in principle up to arbitrary order n, and with very high accuracy, due to the remarkable algebraic properties of _nD_v. With the DA method, the enormous effort to derive explicit formulae for the aberration coefficients of electron mirrors is avoided. A software package MIRROR_DA has been developed for the aberration analysis of electron mirrors, based on the DA method. Two examples of electron mirrors are presented. For the first example, for which the electrostatic and magnetic fields are represented by analytical models, the results computed with MIRROR_DA were shown to be in good agreement with those extracted by direct ray tracing, with relative deviations of less than 0.065% for all the primary aberration coefficients. The second example consists of a real magnetic lens and electrostatic mirror, with numerically computed fields, and from the results of MIRROR_DA, the spherical aberration coefficient C_s3 is almost cancelled out because of the correction effect of the mirror. The MIRROR_DA software is a novel, effective and precise tool for the aberration analysis of electron mirrors, capable of handling realistic and complicated systems of electron lenses and electron mirrors.     

Computer optimization of retarding lens systems for ESCA spectrometers

The performance of four-element electrostatic lenses as retarding systems between source and analyzer in ESCA spectrometers is calculated. The potential distribution in the lens is defined by an axial potential of the type (?(z) = V_o + Σ(V_i - V_{i- 1})/2 · tanh (ω/a_i (z - z_i)). For a given general shape of the lens and a given retardation ratio, the potentials of the two middle electrodes are fitted to give a paraxial image with a prescribed magnification at the exit slit of the lens system. The equipotential surfaces forming the electrodes are found by calculating the potential in an off-axis point, using the series expansion. All third-order geometrical and first-order chromatic aberrations of the lenses are calculated and used together with the second-order aberrations of the analyzer to calculate optimum dimensions of the lens elements and of the emittance-defining slits. A computer program, of which one part calculates the lens properties and one the properties of the entire system lens-analyzer, is described.Two lens systems are presented in some detail. The first one is intended for use with a hemispherical electrostatic analyzer. The angular acceptance is here defined by an aperture stop inside the lens. In this system, the image position and magnification can be kept constant for retardation ratios at least between 1:2 and 60:1, with moderate potentials on the middle electrodes. The second lens system is designed for a magnetic spectrometer of the π √2-type. Here, the central trajectory in the lens is slightly curved by the magnetic field, and the angular acceptance is defined by a baffle after the lens. This system is optimized for a constant retardation ratio of 5:1.     

Different algebraic method for computing the high-order aberrations of practical combined focusing-deflection systems

Differential algebraic (DA) method is an effective technique in computer numerical analysis. It implements conveniently differentiation up to arbitrary high order, based on the nonstandard analysis. Some complicated nonlinear dynamics problem including high-order aberrations of electron optical systems can be solved effectively by mapping properties of DA quantities. However, the existing study applying DA method to the practical system is limited to simple electron lenses. In this paper, the application of DA method is extended to practical combined focusing-deflection systems, and the aberrations up to fifth order are calculated. The electric and magnetic fields of the electron lenses and deflectors, which are computed by finite element method (FEM) or finite difference method (FDM), are in the form of discrete arrays. Local analytical expressions for electric and magnetic field quantities are constructed from these arrays for arbitrary place where electron is traced in DA computation. Thus the developed DA method is applicable for engineering design and computing problem. Based on the study of the expressions and the structure of the aberration coefficients for the fifth-order aberrations of the combined focusing-deflection system and the local analytic expression of the practical electromagnetic field, correlative computer software was developed, whose interface is convenient to calculate the high-order aberrations of the practical systems. The new DA method and software are tested with a uniform magnetic field deflection having an analytic solution. And the results show that the accuracy is sufficiently high, only limited by machine precision. Finally, as an example, a practical combined magnetic focusing and magnetic deflection system is analyzed and discussed.     

Mobility matrix for arbitrary spherical particles in solution

We review the theory of the extended mobility matrix for N arbitrary, spherically symmetric particles immersed in an incompressible fluid. The two-particle mobility functions can be evaluated to any desired order in the inverse interparticle distance by means of an algebraic computer program implementing exact recursion relations. We correct some earlier published expressions and summarize known results for the scattering coefficients which characterize the hydrodynamic properties of the particles. Explicit results are presented for stick and slip hard spheres, for permeable spheres and for fluid droplets.     

Study on differential algebraic method for electron optical properties of wide electron beams in immersion electrostatic-magnetic lenses

Differential algebraic method for calculation of electron optical properties is extended to wide electron beams in combined immersion electrostatic-magnetic lenses, based on Taylor series expansion of the differential algebraic data structure of wide electron beam focusing systems. Higher order aberrations referred to off-axis central trajectories are calculated and this makes the evaluation of wide electron beam focusing properties much simplified with high-enough precision. A program was written, and a specific combined immersion electrostatic-magnetic lens is calculated as an example. Computed results show that the second aperture aberrations referred to off-axis central trajectories of the system of the primary one, whereas the third- and over third-order ones only play a secondary role.     

Veselago’s lens consisting of left-handed materials with arbitrary index of refraction

We study Veselago’s lens with arbitrary index of refraction and characteristic impedance. Using a full wave optics calculation, we show that this lens can be considered as an imaging system and we derive the appropriate lens formula. The lens with arbitrary index and impedance retains some of the properties of the matched lens, such as the invariance of its optical axis, three-dimensional imaging and easy manufacturing, but it loses the property of sub-wavelength resolution. We also show that identical results can be obtained for the impedance matched lens in the framework of paraxial geometrical optics, from which it can be inferred that optical systems containing such a lens can be studied and designed using traditional ray-tracing tools.     

Lie algebraic analysis for the nonlinear transportof intense pulsed beams in electrostatics lenses

The Lie algebraic method is applied to the analysis of the nonlinear transport of an intense pulsed beam in cylindrically symmetrical electrostatic lenses, and particle orbits in a six-dimensional phase space (x, p_x, y, p_y, τ, p_τ) are obtained in the second order approximation. They can also be acquired in the third or higher order approximation if needed. In the analysis, we divide the electrostatic lenses into several segments. Each segment is considered as a uniform accelerating field, and each dividing point is treated as a thin lens. The particle distribution in a three-dimensional ellipsoid is of Gaussian type.     

Computing the Abel map

We present the next step in an ongoing research program to allow for the black-box computation of the so-called finite-genus solutions of integrable differential equations. This next step consists of the black-box computation of the Abel map from a Riemann surface to its Jacobian. Using a plane algebraic curve representation of the Riemann surface, we provide an algorithm for the numerical computation of this Abel map. Since our plane algebraic curves are of arbitrary degree and may have arbitrary singularities, the Abel map of any connected compact Riemann surface may be obtained in this way. This generality is necessary in order for these algorithms to be relevant for the computation of the finite-genus solutions of any integrable equation.     

Algebraic treatment of quantum-mechanical models with modified Coulomb potentials

Nonrelativistic models with modified Coulomb potentials are solved by an algebraic method based on SO(4,2) dynamical group. The nonrelativistic model with Yukawa long-range potential and the Stark effect in the nonrelativistic hydrogen atom in the homogeneous and inhomogeneous external electrostatic fields are studied in details. The algebraic method for the Dirac and Klein-Gordon relativistic hydrogen atom as well as relativistic models with the rotationally symmetric modified Coulomb potentials are also discussed.     

对流扩散方程的涡区分离解法

对流扩散方程是流体计算中一个基本方程,常用的数值方法导至解一个高阶的代数方程组,要求较大的存贮量和较长的计算时间。本文提出一种涡区分离解法,它利用对流扩散方程的迎风性质,把涡区从对流支配区分离出来,仅在各个涡区建立代数方程组并求解。而在对流支配区,则充分利用其抛物性,只需采用显式格式进行计算。由于在各涡区建立的这些方程组阶数和带宽都较小,因此要求存贮量较小,计算速度较快。对于雷诺数较大,涡区范围较小的问题,该方法特别有效。     

The Newest Release of the Ortocartan Set of Programs for Algebraic Calculations in Relativity

The program Ortocartan for algebraic calculations in relativity has just been implemented in the Codemist Standard Lisp and can now be used under the Windows 98 and Linux operating systems. The paper describes the new facilities and subprograms that have been implemented since the previous release in 1992. These are: the possibility to write the output as Latex input code and as Ortocartan's input code, the calculation of the Ellis evolution equations for the kinematic tensors of flow, the calculation of the curvature tensors from given (torsion-free) connection coefficients in a manifold of arbitrary dimension, the calculation of the lagrangian from a given metric by the Landau-Lifshitz method, the calculation of the Euler–Lagrange equations from a given lagrangian (only for sets of ordinary differential equations) and the calculation of first integrals of sets of ordinary differential equations of second order (the first integrals are assumed to be polynomials of second degree in the first derivatives of the functions).     

Bethe ansatz for the XXX-S chain with non-diagonal open boundaries

We consider the algebraic Bethe ansatz solution of the integrable and isotropic XXX-S Heisenberg chain with non-diagonal open boundaries. We show that the corresponding K-matrices are similar to diagonal matrices with the help of suitable transformations independent of the spectral parameter. When the boundary parameters satisfy certain constraints we are able to formulate the diagonalization of the associated double-row transfer matrix by means of the quantum inverse scattering method. This allows us to derive explicit expressions for the eigenvalues and the corresponding Bethe ansatz equations. We also present evidences that the eigenvectors can be build up in terms of multiparticle states for arbitrary S.     

On the integration of some classes of (<Emphasis Type="Italic">n</Emphasis>+1)-dimensional nonlinear equations of mathematical physics

This paper presents a method for construction of the families of particular solutions to some new classes of (n+1)-dimensional nonlinear partial differential equations (PDE). The method is based on the specific link between algebraic matrix equations and PDEs. Admittable solutions depend on arbitrary functions of n variables. Examples of deformed Burgers-type equations are given.     

On the polarizability of a pair of cylinders in a transverse electric field

A sequential scheme for calculating the polarizability of a pair of parallel cylinders with an arbitrary (but fairly symmetric) form is suggested. An infinite set of algebraic equations is derived for coefficients involved in an expression for the potential. The numerical solution of this set makes it possible to find the polarizability with any degree of accuracy. This method is unrelated to the cylinder shape. The electrostatic properties of the cylinders are described in terms of the multipolar polarizability matrix.     

Focusing ground-state atoms with an electrostatic field

A general formula for the trajectory of atoms in an arbitrary potential with axial symmetry is derived. We apply this formula to show that an axially symmetric electrostatic imaging lens for ground-state neutral atoms is possible. Because of its simple construction and capability of focusing atoms of any species, such a lens system can be used in a variety of applications. PACS 03.75.Be; 39.25.+k; 32.80.Pj     

Elliptic PDE formulation and boundary conditions of the spherical harmonics method of arbitrary order for general three-dimensional geometries

The inherent complexity of the radiative transfer equation makes the exact treatment of radiative heat transfer impossible even for idealized situations and simple boundary conditions. Therefore, a wide variety of efficient solution methods have been developed for the RTE. Among these solution methods the spherical harmonics method, the moment method, and the discrete ordinates method provide means to obtain higher-order approximate solutions to the equation of radiative transfer. Although the assembly of the governing equations for the spherical harmonics method requires tedious algebra, their final form promises great accuracy for any given order, since it is a spectral method (rather than finite difference/finite volume in the case of discrete ordinates). In this study, a new methodology outlined in a previous paper on the spherical harmonics method (P_N) is further developed. The new methodology employs successive elimination of spherical harmonic tensors, thus reducing the number of first-order partial differential equations needed to be solved simultaneously by previous P_N approximations (=(N+1)²). The result is a relatively small set (=N(N+1)/2) of second-order, elliptic partial differential equations, which can be solved with standard PDE solution packages. General boundary conditions and supplementary conditions using rotation of spherical harmonics in terms of local coordinates are formulated for the general P_N approximation for arbitrary three-dimensional geometries. Accuracy of the P_N approximation can be further improved by applying the “modified differential approximation” approach first developed for the P₁-approximation. Numerical computations are carried out with the P₃ approximation for several new two-dimensional problems with emitting, absorbing, and scattering media. Results are compared to Monte Carlo solutions and discrete ordinates simulations and a discussion of ray effects and false scattering is provided.     

New Exact Travelling Wave Solutions for Generalized Zakharov-Kuzentsov Equations Using General Projective Riccati Equation Method

Applying the generalized method, which is a direct and unified algebraic method for constructing multiple travelling wave solutions of nonlinear partial differential equations (PDEs), and implementing in a computer algebraic system, we consider the generalized Zakharov-Kuzentsov equation with nonlinear terms of any order. As a result, we can not only successfully recover the previously known travelling wave solutions found by existing various tanh methods and other sophisticated methods, but also obtain some new formal solutions. The solutions obtained include kink-shaped solitons, bell-shaped solitons, singular solitons, and periodic solutions.     

How behavior of systems with sparse spectrum can be predicted on a quantum computer

Call a spectrum of Hamiltonian H sparse if each eigenvalue can be quickly restored within ε from its rough approximation within ε₁ by means of some classical algorithm. It is shown how the behavior of a system with a sparse spectrum up to time T=(1?ρ)/14ε can be predicted on a quantum computer with the time complexity t=4/(1?ρ)ε₁ plus the time of classical algorithm, where ρ is the fidelity. The quantum knowledge of Hamiltonian eigenvalues is considered as the new Hamiltonian W _H whose action on each eigenvector of H gives the corresponding eigenvalue. Speedup of evolution for systems with a sparse spectrum is possible, because, for such systems, the Hamiltonian W _H can be quickly simulated on the quantum computer. For an arbitrary system (even in the classical case), its behavior cannot be predicted on a quantum computer even for one step ahead. By this method, we can also restore the history with the same efficiency.