Successive column correction algorithms for solving sparse nonlinear systems of equations

This paper presents two algorithms for solving sparse nonlinear systems of equations: the CM-successive column correction algorithm and a modified CM-successive column correction algorithm. Aq-superlinear convergence theorem and anr-convergence order estimate are given for both algorithms. Some numerical results and the detailed comparisons with some previously established algorithms show that the new algorithms have some promise of being very effective in practice.This research was partially supported by contracts and grants: DOE DE-AS05-82ER1-13016, AFOSR 85-0243 at Rice University, Houston, U.S.A. and Natural Sciences and Engineering Research Council of Canada grant A-8639.     

Nonlinear Time Series Analysis Since 1990:Some Personal Reflections

Abstract I reflect upon the development of nonlinear time series analysis since 1990 by focusing on five majorareas of development. These areas include the interface between nonlinear time series analysis and chaos,thenonparametric/semiparametric approach,nonlinear state space modelling,financial time series and nonlinearmodelling of panels of time series.     

New feature in the Auger peak of adsorbed oxygen

Recent joint theoretical and experimental investigations of Auger core-core-valence spectra of alkali adatoms on simple metals have revealed that such technique is capable to ascertain contributions from different adsorption environments in the signal [M.I. Trioni, S. Caravati, G.P. Brivio, L. Floreano, F. Bruno, A. Morgante, Phys. Rev. Lett. 93 (2004) 206802]. Consequently, to verify if such an effect is present also for other chemical species, we study theoretically the KLV transition of oxygen either as a bulk impurity or as an adsorbate in/on Al and Ag (jellium-like). We make use of the Fermi golden rule in which the matrix elements of the interaction are calculated within DFT. We verify that the relevant physical quantity of this phenomenon is the excited local density of states (LDOS), calculated within a region centered on the core ionized atom. The Auger rate for oxygen in Ag bulk displays a single asymmetric peak, while for adsorbed oxygen a second smaller feature at lower energies, and very close to the first one, appears. This unexpected result follows from the removal of the degeneracy of the m quantum number of the 2p states of oxygen at the surface. It is only displayed on the electronically less dense metal (Ag), but not on Al.     

差分流线扩散法的局域稳定性分析

本文讨论差分流线扩散法的局部稳定性分析，在理论上证明了该方法依旧保持着流线扩散法扰动局部影响的特性。     

Finite valued feedback laws and piecewise classical solutions

Often, in engineering literature, we find control systems in which the open loop inputs are piecewise constant and take values in a finite set. Such open loop inputs cause the system to have fairly regular solutions. On the other hand, when acting in closed loop, feedback laws taking values in a finite set may not be reinterpreted as open loop inputs of the considered type. In fact, pathological behaviours such as the accumulation of discontinuities may appear (Zeno phenomenon). We give some conditions which can be used as tools for building finite valued feedback laws not causing such pathological behaviours.     

Preliminary results of a non-spherical aerosol method for the retrieval of the atmospheric aerosol optical properties

There is experimental evidence that the non-sphericity of certain atmospheric particles can cause scattering properties different from those predicted by standard Mie theory. Numerous studies indicate the need to consider the presence of non-spherical particles in modeling the optical properties of atmospheric aerosols. On the other hand, natural aerosols show a great variety of shapes, making difficult a realistic choice of a particle shape (or shape mixture) model. In this paper, we test a parameterization of the particle shape in the retrieval of size distribution, phase function, single scattering albedo and asymmetry parameter from direct and sky-radiance measurements. For this purpose we have substituted the Kernel based on the Mie theory included in the model SKYRAD.PACK by one derived for non-spherical particles. The method is applied under different atmospheric conditions, including Saharan dust outbreak, polluted and local mineral episodes. We compare the results with those obtained by the well known spheroids algorithm used in the AERONET network.     

The six-vertex model eigenvectors as critical limit of the eight-vertex model bethe ansatz

The critical limit of the eight-vertex model eigenvectors obtained by means of the generalized Bethe Ansatz is shown to give the six-vertex eigenvectors as constructed in a previous paper by two of the authors. Furthermore, an explicit mapping is established between these eigenvectors and the usual Bethe Ansatz eigenvectors of the six-vertex model. This allows one to show that the indexv labeling the eight-vertex eigenstates becomes exactly the third component of the total spin in the critical limit.     

Optimized analysis of the critical behavior in polymer mixtures from Monte Carlo simulations

A complete outline is given for how to determine the critical properties of polymer mixtures with extrapolation methods similar to the Ferrenberg-Swendsen techniques recently devised for spin systems. By measuring not only averages but the whole distribution of the quantities of interest, it is possible to extrapolate the data obtained in only a few simulations nearT _c over the entire critical region, thereby saving at least 90% of the computer time normally needed to locate susceptibility peaks or cumulant intersections and still getting more precise results. A complete picture of the critical properties of polymer mixtures in the thermodynamic limit is then obtained with finite-size scaling functions. Since the amount of information extracted from a simulation in this way is drastically increased as compared to conventional methods, the investigation of mixtures with long chains or built-in asymmetries is now possible. As an example, the critical points, exponents, and amplitudes of dense, symmetric polymer mixtures with chain lengths ranging fromN=16 up toN=256 are determined within the framework of the 3D bond fluctuation model using grand canonical simulation techniques. As an example for an asymmetry, the generalization of the method to asymmetric monomer potentials is briefly discussed.     

Application of the finite Hankel transform to a diffusion problem without azimuthal symmetry

The problem treated in this paper concerns calculating the evolution of the pressure in a single-phase, slightly compressible fluid in a porous medium consisting of communicating layers. The fluid is produced through a point sink located on the side of an otherwise sealed cylindrical wellbore. This location of the sink causes the flow around the wellbore to be azimuthally asymmetric.The problem is solved through successive application of Laplace, finite Fourier and finite Hankel transforms. Although apparently straightforward, this approach leads to serious numerical difficulties. The published form of the inversion formula for the finite Hankel transform leads to inaccurate computation for the higher azimuthal modes even with 128 bit arithmetic. An alternative form is developed which enables accurate evaluation of the solution with the more practical 64 bit arithmetic. The technique for two-layer solution presented here can be directly extended to a problem with a larger number of communicating layers. This is the first instance of successful application of the finite Hankel transform to an azimuthally asymmetric diffusion problem.