Linearization-preserving self-adjoint and symplectic integrators

This article is concerned with geometric integrators which are linearization-preserving, i.e. numerical integrators which preserve the exact linearization at every fixed point of an arbitrary system of ODEs. For a canonical Hamiltonian system, we propose a new symplectic and self-adjoint B-series method which is also linearization-preserving. In a similar fashion, we show that it is possible to construct a self-adjoint and linearization-preserving B-series method for an arbitrary system of ODEs. Some numerical experiments on Hamiltonian ODEs are presented to test the behaviour of both proposed methods. This work was supported by the Australian Research Council and by the Marsden Fund of the Royal Society of New Zealand.     

High-Order Symplectic and Symmetric Composition Methods for Multi-Frequency and Multi-Dimensional Oscillatory Hamiltonian Systems

The multi-frequency and multi-dimensional adapted Runge-Kutta-Nyström (ARKN) integrators, and multi-frequency and multi-dimensional extended Runge-Kutta-Nyström(ERKN) integrators have been developed to efficiently solve multi-frequency oscillatory Hamiltonian systems. The aim of this paper is to analyze and derive high-order symplectic and symmetric composition methods based on the ARKN integrators and ERKN integrators. We first consider the symplecticity conditions for the multi-frequency and multi-dimensional ARKN integrators. We then analyze the symplecticity of the adjoint integrators of the multi-frequency and multi-dimensional symplectic ARKN integrators and ERKN integrators, respectively. On the basis of the theoretical analysis and by using the idea of composition methods, we derive and propose four new high-order symplectic and symmetric methods for the multi-frequency oscillatory Hamiltonian systems. The numerical results accompanied in this paper quantitatively show the advantage and efficiency of the proposed high-order symplectic and symmetric methods.     

Invariant norm quantifying nonlinear content of Hamiltonian systems

Given a Hamiltonian system, one can represent it using a symplectic map. This symplectic map is specified by a set of homogeneous polynomials which are uniquely determined by the Hamiltonian. In this paper, we construct an invariant norm in the space of homogeneous polynomials of a given degree. This norm is a function of parameters characterizing the original Hamiltonian system. Such a norm has several potential applications.     

Preservation of equilibria for symplectic methods applied to Hamiltonian systems

In this paper, the linear stability of symplectic methods for Hamiltonian systems is studied. In par- ticular, three classes of symplectic methods are considered： symplectic Runge-Kutta （SRK） methods, symplectic partitioned Runge-Kutta （SPRK） methods and the composition methods based on SRK or SPRK methods. It is shown that the SRK methods and their compositions preserve the ellipticity of equilibrium points uncondi- tionally, whereas the SPRK methods and their compositions have some restrictions on the time-step.     

Exponential Rosenbrock integrators for option pricing

In this paper, we are concerned with the time integration of differential equations modeling option pricing. In particular, we consider the Black-Scholes equation for American options. As an alternative to existing methods, we present exponential Rosenbrock integrators. These integrators require the evaluation of the exponential and related functions of the Jacobian matrix. The resulting methods have good stability properties. They are fully explicit and do not require the numerical solution of linear systems, in contrast to standard integrators. We have implemented some numerical experiments in Matlab showing the reliability of the new method.     

New open modified Newton Cotes type formulae as multilayer symplectic integrators

New modified open Newton Cotes integrators are introduced in this paper. For the new proposed integrators the connection between these new algorithms, differential methods and symplectic integrators is studied. Much research has been done on one step symplectic integrators and several of them have obtained based on symplectic geometry. However, the research on multistep symplectic integrators is very poor. Zhu et al. [1] studied the well known open Newton Cotes differential methods and they presented them as multilayer symplectic integrators. Chiou and Wu [2] studied the development of multistep symplectic integrators based on the open Newton Cotes integration methods. In this paper we introduce a new open modified numerical method of Newton Cotes type and we present it as symplectic multilayer structure. The new obtained symplectic schemes are applied for the solution of Hamilton’s equations of motion which are linear in position and momentum. An important remark is that the Hamiltonian energy of the system remains almost constant as integration proceeds. We have applied also efficiently the new proposed method to a nonlinear orbital problem and an almost periodic orbital problem.     

Conservation properties of numerical integrators for highly oscillatory Hamiltonian systems

^** Email: David.Cohen{at}math.unige.ch. Present address: Mathematisches Institut, Universität Tübingen, D-72076 Tübingen, Germany (cohen{at}na.uni-tuebingen.de) Modulated Fourier expansion is used to show long-time near-conservationof the total and oscillatory energies of numerical methods forHamiltonian systems with highly oscillatory solutions. The numericalmethods considered are an extension of the trigonometric methods.A brief discussion of conservation properties in the continuousproblem and in the multi-frequency case is also given.     

Preservation of stability properties near fixed points of linear Hamiltonian systems by symplectic integrators

Based on reasonable testing model problems, we study the preservation by symplectic Runge-Kutta method (SRK) and symplectic partitioned Runge-Kutta method (SPRK) of structures for fixed points of linear Hamiltonian systems. The structure-preservation region provides a practical criterion for choosing step-size in symplectic computation. Examples are given to justify the investigation.     

A molecular dynamics model for symplectic intergrators

Mechanical systems in the very large scale like in celestial mechanics or in the very small scale like in the molecular dynamics can be modelled without dissipation. The resulting Hamiltonian systems possess conservation properties, which are characterized with the term symplecticness, Numerical integration schemes should preserve the symplecticness. Different methods are introduced and their performance is studied for constant and variable step size. As test examples two systems from molecular dynamics are introduced.     

Two Novel Classes of Arbitrary High-Order Structure-Preserving Algorithms for Canonical Hamiltonian Systems

In this paper, we systematically construct two classes of structure-preserving schemes with arbitrary order of accuracy for canonical Hamiltonian systems. The one class is the symplectic scheme, which contains two new families of parameterized symplectic schemes that are derived by basing on the generating function method and the symmetric composition method, respectively. Each member in these schemes is symplectic for any fixed parameter. A more general form of generating functions is introduced, which generalizes the three classical generating functions that are widely used to construct symplectic algorithms. The other class is a novel family of energy and quadratic invariants preserving schemes, which is devised by adjusting the parameter in parameterized symplectic schemes to guarantee energy conservation at each time step. The existence of the solutions of these schemes is verified. Numerical experiments demonstrate the theoretical analysis and conservation of the proposed schemes.     

Rational Krylov methods in exponential integrators for European option pricing

The aim of this paper is to analyze efficient numerical methods for time integration of European option pricing models. When spatial discretization is adopted, the resulting problem consists of an ordinary differential equation that can be approximated by means of exponential Runge–Kutta integrators, where the matrix‐valued functions are computed by the so‐called shift‐and‐invert Krylov method. To our knowledge, the use of this numerical approach is innovative in the framework of option pricing, and it reveals to be very attractive and efficient to solve the problem at hand. In this respect, we propose some a posteriori estimates for the error in the shift‐and‐invert approximation of the core‐functions arising in exponential integrators. The effectiveness of these error bounds is tested on several examples of interest. They can be adopted as a convenient stopping criterion for implementing the exponential Runge–Kutta algorithm in order to perform time integration. Copyright © 2013 John Wiley & Sons, Ltd.     

Weyl-Titchmarsh theory for symplectic difference systems

In this work, we establish Weyl-Titchmarsh theory for symplectic difference systems. This paper extends classical Weyl-Titchmarsh theory and provides a foundation for studying spectral theory of symplectic difference systems.     

Hamiltonian torus actions on symplectic orbifolds and toric varieties

In the first part of the paper, we build a foundation for further work on Hamiltonian actions on symplectic orbifolds. Most importantly, we prove the orbifold versions of the abelian connectedness and convexity theorems. In the second half, we prove that compact symplectic orbifolds with completely integrable torus actions are classified by convex simple rational polytopes with a positive integer attached to each open facet and that all such orbifolds are algebraic toric varieties.

     

Regular reduction of controlled Hamiltonian system with symplectic structure and symmetry

In this paper, our goal is to study the regular reduction theory of regular controlled Hamiltonian (RCH) systems with symplectic structure and symmetry, and this reduction is an extension of regular symplectic reduction theory of Hamiltonian systems under regular controlled Hamiltonian equivalence conditions. Thus, in order to describe uniformly RCH systems defined on a cotangent bundle and on the regular reduced spaces, we first define a kind of RCH systems on a symplectic fiber bundle. Then we introduce regular point and regular orbit reducible RCH systems with symmetry by using momentum map and the associated reduced symplectic forms. Moreover, we give regular point and regular orbit reduction theorems for RCH systems to explain the relationships between RpCH-equivalence, RoCH-equivalence for reducible RCH systems with symmetry and RCH-equivalence for associated reduced RCH systems. Finally, as an application we regard rigid body and heavy top as well as them with internal rotors as the regular point reducible RCH systems on the rotation group SO(3) and on the Euclidean group SE(3), as well as on their generalizations, respectively, and discuss their RCH-equivalence. We also describe the RCH system and RCH-equivalence from the viewpoint of port Hamiltonian system with a symplectic structure.     

辛对称Hamilton算子值域的闭性

利用扰动理论和算子矩阵的因式分解,研究了辛对称Hamilton算子值域的闭性.针对对角占优、上行占优等情形,在一定条件下给出了值域闭性的若干描述,并得到了一般情形的结果.     

Numerical invariant tori of symplectic integrators for integrable Hamiltonian systems

In this paper, we study the persistence of invariant tori of integrable Hamiltonian systems satisfying Rssmann's non-degeneracy condition when symplectic integrators are applied to them. Meanwhile, we give an estimate of the measure of the set occupied by the invariant tori in the phase space. On an invariant torus,numerical solutions are quasi-periodic with a diophantine frequency vector of time step size dependence. These results generalize Shang's previous ones(1999, 2000), where the non-degeneracy condition is assumed in the sense of Kolmogorov.     

Computation of periodic solutions of Hamiltonian systems

This paper attempts to give a practical method to compute global periodic solutions of autonomous Hamiltonian systems of arbitrary finite order. The proposed numerical method is based on continuation of solutions branching from equlibrium points and requires no iterations. Moreover, during computation of one-parameter families of periodic orbits, their possible bifurcations are determined as well.     

BACKWARD ERROR ANALYSIS OF SYMPLECTIC INTEGRATORS FOR LINEAR SEPARABLE HAMILTONIAN SYSTEMS

Symplecticness, stability, and asymptotic properties of Runge-Kutta, partitioned Runge-Kutta, and Runge-Kutta-Nystrom methods applied to the simple Hamiltonian system p = -vg, q = kp are studied. Some new results in connection with P-stability are presented. The main part is focused on backward error analysis. The numerical solution produced by a symplectic method with an appropriate stepsize is the exact solution of a perturbed Hamiltonian system at discrete points. This system is studied in detail and new results are derived. Numerical examples are presented.     

Homoclinic solutions for a class of subquadratic second-order Hamiltonian systems

In this paper, we study the existence of infinitely many homoclinic solutions for a class of subquadratic second-order Hamiltonian systems. By using the variant fountain theorem, we obtain a new criterion for guaranteeing that second-order Hamiltonian systems has infinitely many homoclinic solutions. Recent results from the literature are generalized and significantly improved. An example is also given in this paper to illustrate our main results.