Symmetric and symplectic exponential integrators for nonlinear Hamiltonian systems
This letter studies symmetric and symplectic exponential integrators when applied to numerically computing nonlinear Hamiltonian systems. We first establish the symmetry and symplecticity conditions of exponential integrators and then show that these conditions are extensions of the symmetry and symplecticity conditions of Runge-Kutta methods. Based on these conditions, some symmetric and symplectic exponential integrators up to order four are derived. Two numerical experiments are carried out and the results demonstrate the remarkable numerical behavior of the new exponential integrators in comparison with some symmetric and symplectic Runge-Kutta methods in the literature.
Keywords: exponential integrators; symmetric methods; symplectic methods; Hamiltonian systems
MSC (2000): 65L05, 65P10
In this letter, we explore efficient symmetric and symplectic methods for solving the initial value problems expressed in the following from
where is assumed to be a linear operator on a Banach space with a norm , is the infinitesimal generator of a strongly continuous semigroup on and the function is analytic (see, e.g. ). It follows from the assumption of that there exist two constants and such that
We note that the linear operator can be a matrix if is chosen as or . Under this situation, is accordingly the matrix exponential function. It is known that the exact solution of (1) can be represented by the variation-of-constants formula
Problems of the form (1) often arise in a wide range of practical applications such as quantum physics, engineering, flexible body dynamics, mechanics, circuit simulations and other applied sciences (see, e.g. [5, 7, 18, 21]). Some highly oscillatory problems, Schrödinger equations, parabolic partial differential equations with their spatial discretisations all fit the form. In order to solve (1) effectively, many researches have been done and the readers are referred to [9, 10, 11, 12, 13] for example. Among them, a standard form of exponential integrators is formulated and these integrators have been studied by many researchers. We refer the reader to [1, 2, 3, 4, 5, 8, 15, 17, 19, 20] for some examples on this topic and a systematic survey of exponential integrators is referred to .
On the other hand, it can be observed that the problem (1) can become a nonlinear Hamiltonian system if
where is a smooth potential function, is a symmetric matrix, and with the identity . The energy of this Hamiltonian system is . For this system, symplectic exponential integrators (see ) are strongly recommended since they can preserve the symplecticity of the original problems and provide good long time energy preservation and stability. Besides, it is shown in  that symmetric methods also have excellent long time behaviour when applied to reversible differential equations and symmetric exponential integrators have been considered for solving Schrödinger equations in . However, it seems that symmetric exponential integrators have not been used for ODEs and moreover symmetric and symplectic exponential integrators have never been studied so far, which motives this letter.
The main contribution of this letter is to analyse and derive symmetric and symplectic exponential integrators. The integrators have symmerty and symplecticity simultaneously. The letter is organized as follows. In Section 2, we present the scheme of exponential integrators and derive its properties including the symmetry and symplecticity conditions. Then we are devoted to the construction of some practical symmetric and symplectic exponential integrators in Section 3. In Section 4, we carry out two numerical experiments and the numerical results demonstrate the remarkable efficiency of the new integrators in comparison with some existing methods in the scientific literature. The last section is concerned with conclusions.
2 Exponential integrators and their properties
In this section, we first present the scheme of exponential integrators and then analyze their symmetry and symplecticity conditions.
It is worth mentioning that if an exponential integrator reduces to a classical Runge-Kutta (RK) method with the coefficients for .
The next theorem gives the symmetry conditions of the EI method.
An s-stage exponential integrator is symmetric if and only if its coefficients satisfy the following conditions:
Proof For the exponential integrator , exchanging and replacing by yields
Then we have
For the second formula of (7) and the first formula of , it is required that the following conditions are true
Based on these conditions, we obtain that the following two formulae are equal
This implies that
It is noted that when , these symmetry conditions become
which are the exact symmetry conditions of -stage RK methods.
About the symplecticity conditions of exponential integrators, we have the following result.
(See .) If the coefficients of an -stage exponential integrator satisfy
where and for then the integrator is symplectic.
We also remark that when , these conditions reduce to
which are the exact symplecticity conditions of -stage RK methods.
3 Symmetric and symplectic EI
In this section, we derive a class of symmetric and symplectic exponential integrators. To this end, we consider the following special exponential integrators.
(See ) Define a special kind of -stage exponential integrators by
are the coefficients of an -stage RK method. We denote this class of exponential integrators by SEI.
For these special exponential integrators, the following properties can be derived.
which can be simplified as the symmetry conditions (10) of RK methods. Thus the first statement is true. The second statement can be obtained immediately by considering Theorem 3.2 of . Based on the above two results, the third one holds. The last result comes from Theorem 3.1 of .
In what follows, based on Theorem 3.2 we construct some practical symmetric and symplectic SEI integrators.
3.1 One-stage symmetric and symplectic SEI
First consider an one-stage RK method with the coefficients:
From these formulae, it follows that
This gives a class of symmetric and symplectic exponential integrators by considering (13). As an example, we choose and denote the method as SSSEI1s2. It can be checked that this RK method is implicit midpoint rule. Thus the symmetric and symplectic SEI is of order two.
3.2 Two-stage symmetric and symplectic SEI
Consider a two-stage RK method whose coefficients are given by a Butcher tableau:
The symmetry conditions of this method are
The RK method is symplectic if
In the light of the third-order conditions of RK methods (see )
we choose the parameters by
This choice as well as (19) and (13) gives an symmetric and symplectic SEI. Moreover, it can be seen that the corresponding RK method is Gauss method of order four (see ). Therefore, the SEI is also of order four, which is denoted by SSSEI2s4.
3.3 Three-stage symmetric and symplectic SEI
We turn to considering three-stage symmetric and symplectic integrators. The following Butcher tableau describes a three-stage RK method:
This method is symmetric if the following conditions are true
The RK methods are symplectic if
and denote the method as SSSEI3s4. According to the fourth-order conditions of RK methods (see )
it can be checked that the coefficients of the RK method satisfy all the conditions. Thus this RK method is of order four and the symmetric and symplectic SEI has same order. This symmetric and symplectic SEI is denoted by SSSEI3s4.
4 Numerical experiments
This section presents two numerical experiments to show the remarkable efficiency of the new integrators as compared with some existing RK methods. The integrators for comparisons are chosen as:
SSSEI1s2: the one-stage symmetric and symplectic EI of order two presented in this letter;
SSSEI2s4: the two-stage symmetric and symplectic EI of order four presented in this letter;
SSSEI3s4: the three-stage symmetric and symplectic EI of order four presented in this letter;
SSRK1s2: the one-stage symmetric and symplectic RK method of order two obtained by letting for SSSEI1s2 (implicit midpoint rule);
SSRK2s4: the two-stage symmetric and symplectic RK method of order four obtained by letting for SSSEI2s4 (Gauss method of order four);
SSRK3s4: the three-stage symmetric and symplectic RK method of order four obtained by letting for SSSEI3s4 (the method was given in ).
Problem 1. As the first numerical example, we consider the Duffing equation defined by
It is a Hamiltonian system with the Hamiltonian The exact solution of this system is with the Jacobi elliptic function . For this problem, we choose , , and for The efficiency curves are shown in Figure 1 (i). We integrate this problem with a fixed stepsize in the interval for . The results of energy conservation are presented in Figure 1 (ii).
Problem 2. The second numerical example is the following averaged system in wind-induced oscillation
where is a detuning parameter and is a damping factor with The first integral (when ) or Lyapunov function (when ) of this system is
We choose the initial values Firstly we consider and solve the problem on the interval with for . The efficiency curves are shown in Figure 2 (i). Then this problem is integrated with on the interval See Figure 2 (ii) for the energy conservation for different methods. Secondly we choose and the efficiency curves are shown in Figure 2 (iii) on with
From the numerical results, it follows clearly that the symmetric and symplectic exponential integrators behave much better than symmetric and symplectic RK methods.
5 Conclusions and discussions
In this letter, in order to solve the differential equations (1) by using symmetric and symplectic methods, we present the symmetry and symplecticity conditions for exponential integrators. Then based on these conditions, we consider a special kind of exponential integrators and construct some practical symmetric and symplectic exponential integrators. The remarkable efficiency of the new integrators is shown by the numerical results from two numerical experiments in comparison with some existing RK methods in the literature.
-  H. Berland, B. Owren, B. Skaflestad, B-series and order conditions for exponential integrators, SIAM J. Numer. Anal. 43 (2005) 1715-1727.
-  M. Caliari, A. Ostermann, Implementation of exponential Rosenbrock-type integrators, Appl. Numer. Math. 59 (2009) 568-581.
-  M.P. Calvo, C. Palencia, A class of explicit multistep exponential integrators for semilinear problems, Numer. Math. 102 (2006) 367-381.
-  E. Celledoni, D. Cohen, B. Owren, Symmetric exponential integrators with an application to the cubic Schrödinger equation, Found. Comput. Math. 8 (2008) 303-317.
-  V. Grimm, M. Hochbruck, Error analysis of exponential integrators for oscillatory second-order differential equations, J. Phys. A: Math. Gen. 39 (2006) 5495-5507.
-  E. Hairer, C. Lubich, G. Wanner, Geometric Numerical Integration: Structure-Preserving Algorithms for Ordinary Differential Equations, 2nd edn. Springer-Verlag, Berlin, Heidelberg, 2006.
-  M. Hochbruck, A. Ostermann, Exponential integrators, Acta Numer. 19 (2010) 209-286.
-  M. Hochbruck, A. Ostermann, J. Schweitzer, Exponential rosenbrock-type methods, SIAM J. Numer. Anal. 47 (2009) 786-803.
-  A. Iserles, On the global error of discretization methods for highly-oscillatory ordinary differential equations, BIT Num. Math. 42 (2002) 561-599.
-  A. Iserles, Think globally, act locally: solving highly-oscillatory ordinary differential equations, Appl. Num. Anal. 43 (2002) 145-160.
-  A.K. Kassam, L.N.Trefethen, Fourth-order time-stepping for stiff PDEs, SIAM J. Sci. Comput. 26 (2005) 1214-1233.
-  M. Khanamiryan, Quadrature methods for highly oscillatory linear and nonlinear systems of ordinary differential equations: part I, BIT Num. Math. 48 (2008) 743-762.
-  S. Krogstad, Generalized integrating factor methods for stiff PDEs, J. Comput. Phys. 203 (2005) 72-88.
-  L. Mei, X. Wu, Symplectic exponential Runge-Kutta methods for solving nonlinear Hamiltonian systems, J. Comput. Phys. 338 (2017) 567–584.
-  A. Ostermann, M. Thalhammer, W.M. Wright, A class of explicit exponential general linear methods, BIT Numer. Math. 46 (2006) 409-431.
-  J.M. Sanz-Serna, L. Abia, Order conditions for canonical Runge–Kutta schemes, SIAM J. Numer. Analy. 28 (1991) 1081-1096.
-  B. Wang, A. Iserles, X. Wu, Arbitrary–order trigonometric Fourier collocation methods for multi-frequency oscillatory systems, Found. Comput. Math. 16 (2016) 151-181.
-  B. Wang, T. Li, Y. Wu, Arbitrary-order functionally fitted energy-diminishing methods for gradient systems, Appl. Math. Lett. 83 (2018) 130-139
-  B. Wang, X. Wu, The formulation and analysis of energy-preserving schemes for solving high-dimensional nonlinear Klein-Gordon equations. IMA. J. Numer. Anal. DOI: 10.1093/imanum/dry047 (2018)
-  B. Wang, H. Yang, F. Meng, Sixth order symplectic and symmetric explicit ERKN schemes for solving multi-frequency oscillatory nonlinear Hamiltonian equations, Calcolo, 54 (2017) 117-140
-  X. Wu, B. Wang, Recent Developments in Structure-Preserving Algorithms for Oscillatory Differential Equations. Springer Nature Singapore Pte Ltd, 2018