Scaling by Binormalization

We present an iterative algorithm (BIN) for scaling all the rows and columns of a real symmetric matrix to unit 2-norm. We study the theoretical convergence properties and its relation to optimal conditioning. Numerical experiments show that BIN requires 2–4 matrix–vector multiplications to obtain an adequate scaling, and in many cases significantly reduces the condition number, more than other scaling algorithms. We present generalizations to complex, nonsymmetric and rectangular matrices.     

Exact solutions and geometric phase factor of time-dependent three-generator quantum systems

There exist a number of typical and interesting systems and/or models, which possess three-generator Lie-algebraic structure, in atomic physics, quantum optics, nuclear physics and laser physics. The well-known fact that all simple 3-generator algebras are either isomorphic to the algebra sl (2, C) or to one of its real forms enables us to treat these time-dependent quantum systems in a unified way. By making use of both the Lewis-Riesenfeld invariant theory and the invariant-related unitary transformation formulation, the present paper obtains exact solutions of the time-dependent Schr?dinger equations governing various three-generator Lie-algebraic quantum systems. For some quantum systems whose time-dependent Hamiltonians have no quasialgebraic structures, it is shown that the exact solutions can also be obtained by working in a sub-Hilbert-space corresponding to a particular eigenvalue of the conserved generator (i.e., the time-independent invariant that commutes with the time-dependent Hamiltonian). The topological property of geometric phase factors and its adiabatic limit in time-dependent systems is briefly discussed. Received 6 July 2002 / Received in final form 21 October 2002 Published online 11 February 2003     

The numerical approximation of center manifolds in Hamiltonian systems

In this paper we develop a numerical method for computing higher order local approximations of center manifolds near steady states in Hamiltonian systems. The underlying system is assumed to be large in the sense that a large sparse Jacobian at the equilibrium occurs, for which only a linear solver and a low-dimensional invariant subspace is available. Our method combines this restriction from linear algebra with the requirement that the center manifold is parametrized by a symplectic mapping and that the reduced equation preserves the Hamiltonian form. Our approach can be considered as a special adaptation of a general method from Numer. Math. 80 (1998) 1-38 to the Hamiltonian case such that approximations of the reduced Hamiltonian are obtained simultaneously. As an application we treat a finite difference system for an elliptic problem on an infinite strip.     

Geometry and electronic structure of a heterometallic cluster Mo<Subscript>2</Subscript>Mg<Subscript>2</Subscript> in different oxidation states of Mo: a DFT study

Reduction of tetranuclear heterometallic complex Mo₂Mg₂ was simulated using the B3LYP and PBE density functional methods. The results of geometry calculations of the initial complex [Mo^VIO₂Mg(MeOH)₂(OMe)₄]₂ and a partially reduced Mo^V complex are in good agreement with experimental data. The reduced Mo^III complex is characterized by a decrease in the binding energy of aqua ligands. Structural rearrangement of the complex with release of a coordination position at the Mo atoms requires small energy expenditure. One can assume that the reduction of the polynuclear complex causes overcrowding of its coordination sphere, which favors formation of dinitrogen complexes. Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 3, pp. 441–457, March, 2008.     

Difference schemes on uniform grids performed by general discrete operators

Our aim is to set the foundations of a discrete vectorial calculus on uniform n-dimensional grids, that can be easily reformulated on general irregular grids. As a key tool we first introduce the notion of tangent space to any grid node. Then we define the concepts of vector field, field of matrices and inner products on the space of grid functions and on the space of vector fields, mimicking the continuous setting. This allows us to obtain the discrete analogous of the basic first order differential operators, gradient and divergence, whose composition define the fundamental second order difference operator. As an application, we show that all difference schemes, with constant coefficients, for first and second order differential operators with constant coefficients can be seen as difference operators of the form for suitable choices of q, and  . In addition, we characterize special properties of the difference scheme, such as consistency, symmetry and positivity in terms of q, and  .     

Indirect boundary element method for unsteady linearized flow over prolate and oblate spheroids and hemispheroidal protuberances

The indirect boundary element method was used to study the hydrodynamics of oscillatory viscous flow over prolate and oblate spheroids, and over hemispheroidal bodies hinged to a plate. Analytic techniques, such as spheroidal coordinates, method of images, and series representations, were used to make the numerical methods more efficient. A novel method for computing the hydrodynamic torque was used, since for oscillatory flow the torque cannot be computed directly from the weightings. Instead, a Green's function for torque was derived to compute the torque indirectly from the weightings. For full spheroids, the method was checked by comparing the results to exact solutions at low and high frequencies, and to results computed using the singularity method. For hemispheroids hinged to a plate, the method for low frequencies was checked by comparing the results to previous results, and to exact solutions at high frequencies. Copyright © 2004 John Wiley & Sons, Ltd.     

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).     

Kinetics and equilibrium of small metallic clusters: Ab initio confinement molecular dynamics study of 4

The ab initio molecular dynamics (AIMD) [1] is combined with the heuristic, successive confinement method of surveying a potential energy surface (PES) [2], thereby offering a framework for the simulation study of kinetics and equilibrium properties of metallic clusters. This approach is applied to the study of Au₄, a cluster possessing a simple but specific PES, which consists of very shallow and deep basins and due to this presents a challenge to the conventional AIMD methods. Among other things, the probabilities of the transitions between isomers have been found, and on this basis, both the time-dependent and equilibrium populations of the isomers have been calculated for the conditions typical of the NeNePo experiments [3] in the femtosecond pump-probe spectroscopy.     

延迟微分方程单支方法的非线性稳定性

本文讨论延迟微分方程单支方法的非线性稳定性 .对于 Kα,β,γ类非线性延迟微分方程 ,我们证明带有线性插值的 G( c,p,q) -代数稳定的单支方法当 c≤ 1时是 GR( p/2 ,q/2 ) -稳定及弱 GAR( p/2 ,q/2 ) -稳定的 ,当 c<1时是 GAR( p/2 ,q/2 ) -稳定的 .最后的数值试验表明了上述结论的正确性 .