We prove that any linear multi-step method G1^T of the form ∑k=0^mαkZk = T∑k=0^mβkJ^-1↓ΔH(Zk) with odd order u (u≥ 3) cannot be conjugate to a symplectic method G2^T of order w (w 〉 u) via any generalize...We prove that any linear multi-step method G1^T of the form ∑k=0^mαkZk = T∑k=0^mβkJ^-1↓ΔH(Zk) with odd order u (u≥ 3) cannot be conjugate to a symplectic method G2^T of order w (w 〉 u) via any generalized linear multi-step method G3^T of the form ∑k=0^mαkZk = T∑k=0^mβkJ^-1↓ΔH(∑l=0^mγklZl). We also give a necessary condition for this kind of generalized linear multi-step methods to be conjugate-symplectic. We also demonstrate that these results can be easily extended to the case when G3^T is a more general operator.展开更多
By performing a particular spatial discretization to the nonlinear Schrodinger equation(NLSE),we obtain a non-integrable Hamiltonian system which can be decomposed into three integrable parts(L-L-N splitting).We integ...By performing a particular spatial discretization to the nonlinear Schrodinger equation(NLSE),we obtain a non-integrable Hamiltonian system which can be decomposed into three integrable parts(L-L-N splitting).We integrate each part by calculating its phase flow,and develop explicit symplectic integrators of different orders for the original Hamiltonian by composing the phase flows.A 2nd-order reversible constructed symplectic scheme is employed to simulate solitons motion and invariants behavior of the NLSE.The simulation results are compared with a 3rd-order non-symplectic implicit Runge-Kutta method,and the convergence of the formal energy of this symplectic integrator is also verified.The numerical results indicate that the explicit symplectic scheme obtained via L-L-N splitting is an effective numerical tool for solving the NLSE.展开更多
基金Acknowledgements. We would like to thank the editors for their valuable suggestions and corrections. This research is supported by the National Natural Science Foundation of China (Grant Nos. 10471145 and 10672143), and by Morningside Center of Mathematics, Chinese Academy of Sciences.
文摘We prove that any linear multi-step method G1^T of the form ∑k=0^mαkZk = T∑k=0^mβkJ^-1↓ΔH(Zk) with odd order u (u≥ 3) cannot be conjugate to a symplectic method G2^T of order w (w 〉 u) via any generalized linear multi-step method G3^T of the form ∑k=0^mαkZk = T∑k=0^mβkJ^-1↓ΔH(∑l=0^mγklZl). We also give a necessary condition for this kind of generalized linear multi-step methods to be conjugate-symplectic. We also demonstrate that these results can be easily extended to the case when G3^T is a more general operator.
基金This research is partially supported by the Informatization Construction of Knowledge Innovation Projects of the Chinese Academy of Sciences“Supercomputing En-vironment Construction and Application”(INF105-SCE)National Natural Science Foundation of China(Grant Nos.10471145 and 10672143).
文摘By performing a particular spatial discretization to the nonlinear Schrodinger equation(NLSE),we obtain a non-integrable Hamiltonian system which can be decomposed into three integrable parts(L-L-N splitting).We integrate each part by calculating its phase flow,and develop explicit symplectic integrators of different orders for the original Hamiltonian by composing the phase flows.A 2nd-order reversible constructed symplectic scheme is employed to simulate solitons motion and invariants behavior of the NLSE.The simulation results are compared with a 3rd-order non-symplectic implicit Runge-Kutta method,and the convergence of the formal energy of this symplectic integrator is also verified.The numerical results indicate that the explicit symplectic scheme obtained via L-L-N splitting is an effective numerical tool for solving the NLSE.