三级近似束流光学计算程序

为了精确计算束流在离子光学系统中的传输,用Visual FORTRAN 6.5语言编写了一个计算程序,长约13000行. 此程序可以计算由三圆筒单透镜、三膜片单透镜、双元筒透镜、均匀场静电加速管、磁四极透镜、六极磁铁、静电四极透镜、偏转磁铁、螺线管透镜、ExB~正交电磁场分析器、静电偏转器、漂浮管、QWR(Quarter Wave Resonators)和SLR(Split Loop Resonators)射频加速元件等元件任意组成的离子光学系统. 粒子轨迹的计算可精确到三级近似. 粒子的分布类型也可以有多种选择. 程序具有最优化计算功能,即可以自动调整元件的参数,以实现所需要的光学条件. 各元件之后的横向和纵向相图以及系统的束流包络线以图形方式显示在屏幕上.     

ASYMPTOTIC APPROXIMATION METHOD AND ITS CONVERGENCE ON SEMI-INFINITE PROGRAMMING

The aim of this article is to discuss an asymptotic approximation model and its convergence for the minimax semi-infinite programming problem. An asymptotic surrogate constraints method for the minimax semi-infinite programming problem is presented by making use of two general discrete approximation methods. Simultaneously, the consistence and the epi-convergence of the asymptotic approximation problem are discussed.     

关于Tchebycheff—Fourier级数的（H,λ）平均的逼近

我们定义了（Ｈ，λ）求和法，它含有（Ｎ，ｐｎ），（Ｒ＾ｒｎ）和（Ｖｍｎ）求和法。讨论了函数ｆ（ｘ）∈Ｃ＾ｒ［－１，１］（ｒ∈Ｎ０）以及ｆ（ｘ）∈Ｗ＾ｒＨ＾ａ（ｒ∈Ｎ０，０＜ａ＜１）的切比晓夫－富里埃级数的逼近阶。     

The Micro-Theory of Ultra-Short Pulse Propagating in Nonlinear Fiber

The modified quantum nonlinear Schrodinger equation of ultra-short pulse propagating in fiber is derived using dispersion relation and the Hamiltonian of the transmission system. The derived equation is solved with linearization approximation, and modulation instability is analyzed. The equation is also solved with Hartree approximation. The results indicate that pulse power, loss and self-steeping effect change the critical frequency, the maximum gain and the gain spectrum range, but the third order dispersion has no effect on modulation instability. The expectation value of optical field is average of a set of modified classical solitons, and higher order terms change the amplitude, pulse position and phase of soliton.     

Duality in parameter space and approximation of measures for mixing repellers

For one-dimensional expanding mapsT with an invariant measure we consider, in a parameter space, the envelope _n of the real lines associated to any couple of points of the orbit, connected byn iterations ofT. If the map hass inverses and is piecewise linear, then the sets _n are just the union ofs ⁿ points and converge to the invariant Cantor set ofT. A correspondence between all the sets and their measures is established and allows one to associate the atomic measure on ₁ to the completly continuous measure on the Cantor set. If the map is nonlinear, hyperbolic, and hass inverses, the sets _n are homeomorphic to the Cantor set; they converge to the Cantor set ofT and their measures converge to the measure of the Cantor set whenn. The correspondence between the sets _n allows one to define converging approximation schemes for the map an its measure: one replaces each of thes ⁿ disjoint sets with a point in a convenient neighborhood and a probability equal to its measure and transforms it back to the original set ₁. All the approximations with linear Cantor systems previously proposed are recovered, the converging proprties being straightforward in the present scheme. Moreover, extensions to higher dimensionality and to nondisconnected repellers arte possible and are briefly examined.     

Asymptotic behavior of energy band associated with a negative energy level

The mathematical foundation of the tight binding approximation is given. If ₀ is a negative energy level of a real potentialq, then there exists an energy band for a one-dimensional chain with period 2T of the same atoms which lies near ₀. We study this band whenT tends to infinity.On leave of absence from the Department of Physics, Leningrad State University, Leningrad, USSR.     

Hereditary optimal control problems: Numerical method based upon a padé approximation

In this paper, we consider a particular approximation scheme which can be used to solve hereditary optimal control problems. These problems are characterized by variables with a time-delayed argumentx(t – ). In our approximation scheme, we first replace the variable with an augmented statey(t) x(t - ). The two-sided Laplace transform ofy(t) is a product of the Laplace transform ofx(t) and an exponential factor. This factor is approximated by a first-order Padé approximation, and a differential relation fory(t) can be found. The transformed problem, without any time-delayed argument, can then be solved using a gradient algorithm in the usual way. Four problems are solved to illustrate the validity and usefulness of this technique.This research was supported in part by the National Aeronautics and Space Administration under NASA Grant NCC-2-106.     

The projection approach to the Fokker-Planck equation. I. Colored Gaussian noise

It is shown that the Fokker-Planck operator can be derived via a projection-perturbation approach, using the repartition of a more detailed operator into a perturbation ₁ and an unperturbed part ₀. The standard Fokker-Planck structure is recovered at the second order in ₁, whereas the perturbation terms of higher order are shown to provoke the breakdown of this structure. To get rid of these higher order terms, a key approximation, local linearization (LL), is made. In general, to evaluate at the second order in ₁ the exact expression of the diffusion coefficient which simulates the influence of a Gaussian noise with a finite correlation time, a resummation up to infinite order in must be carried out, leading to what other authors call the best Fokker-Planck approximation (BFPA). It is shown that, due to the role of terms of higher order in ₁, the BFPA leads to predictions on the equilibrium distributions that are reliable only up to the first order in t. The LL, on the contrary, in addition to making the influence of terms of higher order in ₁ vanish, results in a simple analytical expression for the term of second order that is formally coincident with the complete resummation over all the orders in t provided by the Fox theory. The corresponding diffusion coefficient in turn is shown to lead in the limiting case to exact results for the steady-state distributions. Therefore, over the whole range 0 the LL turns out to be an approximation much more accurate than the global linearization proposed by other authors for the same purpose of making the terms of higher order in ₁ vanish. In the short- region the LL leads to results virtually coincident with those of the BFPA. In the large- region the LL is a more accurate approximation than the BFPA itself. These theoretical arguments are supported by the results of both analog and digital simulation.     

One-dimensional harmonic lattice caricature of hydrodynamics: Second approximation

The long-time behavior of an infinite chain of coupled harmonic oscillators is studied. In addition to a limiting hydrodynamic (Euler-type) equation, the next approximation is investigated. The corresponding equation is derived.     

Surface properties of finite classical Coulomb systems: Debye-Hückel approximation and computer simulations

We report analytical and numerical studies of surface correlations in finite, homogeneously polarizable, classical Coulomb systems placed in an insulating or conducting environment. Their purpose is to understand the phenomenological, shape-dependent laws of electrostatics, from the point of view of statistical mechanics; we focus on the knowledge of the dielectric susceptibility of the system, a quantity proportional to the equilibrium fluctuation of the system's instantaneous polarization per unit volume. This goal has been achieved for a system in a conducting state. The picture is that the shape-dependent part of the susceptibilities results from the action of unbounded observables (the second moments of the instantaneous polarization of the system) on long-range surface correlations and that the relations of electrostatics are verified by means of shape-dependent thermodynamic limits. This picture is supported (i) by exact solutions and asymptotic analysis of the Debye-Hückel approximation of multicomponent plasmas in disks and spheres with insulating and conducting environment and also in ellipses in a vacuum, and (ii) by computer simulations of a one-component plasma in a disk with different environments, notably a conducting environment with permeable and impermeable wall. These observations have revealed for the first time the reason why the susceptibility of a conducting disk in a conductor with impermeable walls diverges linearly with the radius of the disk: this is due to the occurrence of long-range radial correlations in the conductor. These findings are quantitatively interpreted in terms of a novel canonical Debye-Huckel approximation as contrasted to the ordinary grand canonical version. Lastly a fresh look at the problem of the surface correlations of a conductor in a vacuum, which places the observer close to the surface of the conductor but in the vacuum, is presented and applied to the disk, the ellipse, the cylinder, the sphere, and the wedge.