Geodesic
On the possibility of constructing a conservative numerical method of~solution of the initial value problem for Hamiltonian systems on the basis of~the two-stage symmetric symplectic Runge--Kutta methods
Matematičeskoe modelirovanie, Tome 26 (2014) no. 10, pp. 47-63.

Voir la notice de l'article provenant de la source Math-Net.Ru

The question about the possibility of working out a Hamiltonian system initial value problem solution numerical method which gives an approximate solution that conforms the total energy conservation law is studied in the work. The method is constructed on basis of the family of two-stage symmetric symplectic Runge–Kutta methods. The properties of the suggested method are studied by the example of the model problem about the particle motion in the cubic potential field. The possibility of working out the method giving a numerical solution which preserves the total energy on the period of the problem finite solution, except for small neighbourhood of return points, is shown. The time dependencies of the symplecticity and reversibility defects on the numerical solution obtained by the worked out method are investigated.
Keywords: molecular dynamics, Hamiltonian systems, numerical methods for solving initial value problem, energy conservation, Runge–Kutta methods. 
@article{MM_2014_26_10_a3,
     author = {P. A. Aleksandrov and G. G. Elenin},
     title = {On the possibility of constructing a conservative numerical method of~solution of the initial value problem for {Hamiltonian} systems on the basis of~the two-stage symmetric symplectic {Runge--Kutta} methods},
     journal = {Matemati\v{c}eskoe modelirovanie},
     pages = {47--63},
     publisher = {mathdoc},
     volume = {26},
     number = {10},
     year = {2014},
     language = {ru},
     url = {http://geodesic.mathdoc.fr/item/MM_2014_26_10_a3/}
}
TY  - JOUR
AU  - P. A. Aleksandrov
AU  - G. G. Elenin
TI  - On the possibility of constructing a conservative numerical method of~solution of the initial value problem for Hamiltonian systems on the basis of~the two-stage symmetric symplectic Runge--Kutta methods
JO  - Matematičeskoe modelirovanie
PY  - 2014
SP  - 47
EP  - 63
VL  - 26
IS  - 10
PB  - mathdoc
UR  - http://geodesic.mathdoc.fr/item/MM_2014_26_10_a3/
LA  - ru
ID  - MM_2014_26_10_a3
ER  - 
%0 Journal Article
%A P. A. Aleksandrov
%A G. G. Elenin
%T On the possibility of constructing a conservative numerical method of~solution of the initial value problem for Hamiltonian systems on the basis of~the two-stage symmetric symplectic Runge--Kutta methods
%J Matematičeskoe modelirovanie
%D 2014
%P 47-63
%V 26
%N 10
%I mathdoc
%U http://geodesic.mathdoc.fr/item/MM_2014_26_10_a3/
%G ru
%F MM_2014_26_10_a3
P. A. Aleksandrov; G. G. Elenin. On the possibility of constructing a conservative numerical method of~solution of the initial value problem for Hamiltonian systems on the basis of~the two-stage symmetric symplectic Runge--Kutta methods. Matematičeskoe modelirovanie, Tome 26 (2014) no. 10, pp. 47-63. http://geodesic.mathdoc.fr/item/MM_2014_26_10_a3/

[1] Allen M. P., Tildesley D. J., Computer Simulation of Liquids, Clarendon Press, Oxford, 1991

[2] Deuflhard P., Hermans J., Leimkuhler B. et. al. (eds.), Computational Molecular Dynamics: Challenges, Methods, Ideas, Lecture Notes in Computational Science and Engineering, 4, Berlin, 1998, 494 pp. | MR | Zbl

[3] Hairer E., Lubich C., Wanner G., Geometric Numerical Integration, second revised edition, Springer, Berlin, 2006, 644 pp. | MR | Zbl

[4] Hairer E., Norsett S. P., Wanner G., Solving Ordinary Differential Equations, v. I, Nonstiff Problems, second revised edition, Springer, 1993, 528 pp. | MR | Zbl

[5] Landau L. D., Lifshits E. M., Teoreticheskaya fizika, v. 1, Mekhanika, Nauka, M., 1965

[6] Lasagni F. M., “Canonical Runge–Kutta methods”, Journal of Applied Mathematics and Physics (ZAMP), 39 (1988), 952–953 | DOI | MR | Zbl

[7] Sanz-Serna J. M., “Runge–Kutta schemes for Hamiltonian systems”, BIT, 28:4 (1988), 877–883 | DOI | MR | Zbl

[8] Suris Y. B., “On the conservation of the symplectic structure in the numerical solution of Hamiltonian systems”, Numerical solution of ordinary differential equations, Keldysh Institute of Applied Mathematics, USSR Academy of Sciences, M., 1988, 148–160 | MR | Zbl

[9] Elenin G. G., Aleksandrov P. A., Tochnye resheniya sistemy razreshayuschikh uravnenii se-meistva simmetrichno-simplekticheskikh metodov Runge-Kutty dlya zadachi o dvizhenii materialnoi tochki v pole kubicheskogo potentsiala, preprint, MAKS Press, M., 2011, 58 pp.

[10] Elenin G. G., Shlyakhov P. I., “O konservativnosti dvukhstadiinykh simmetrichno-simplekticheskikh metodov Runge–Kutty i metoda Shtermera–Verle”, Differentsialnye uravneniya, 46:7 (2010), 983–989 | MR | Zbl

[11] Elenin G. G., Aleksandrov P. A., “O konservativnosti dvukhparametricheskogo semeistva trekh-stadiinykh simmetrichno-simplekticheskikh metodov Runge–Kutty”, Differentsialnye uravneniya, 48:7 (2012), 981–989 | Zbl

[12] Oevel W., Sofroniou M., Symplectic Runge–Kutta schemes. II: Classification of symmetric methods, Preprint, University of Paderborn, 1996, 47 pp. | MR

[13] Aleksandrov P. A., Elenin G. G., Chislo reshenii sistemy razreshayuschikh uravnenii semei-stva simmetrichno-simplekticheskikh metodov Runge-Kutty dlya zadachi o dvizhenii materialnoi tochki v pole kubicheskogo potentsiala, preprint, MAKS Press, M., 2010, 38 pp.

[14] Aleksandrov P. A., Elenin G. G., Tochnye resheniya sistemy razreshayuschikh uravnenii metoda srednei tochki dlya zadachi o dvizhenii materialnoi tochki v pole kubicheskogo potentsiala, preprint, MAKS Press, M., 2011, 28 pp.

[15] Attetkov A. V., Galkin S. V., Zarubin V. S., Metody optimizatsii, Ucheb. dlya vuzov, eds. V. S. Zarubin, A. P. Krischenko, Izd-vo MGTU im. N. E. Baumana, M., 2003, 440 pp.

[16] Vasilev F. P., Chislennye metody resheniya ekstremalnykh zadach, Nauka, M., 1988, 552 pp.

[17] Korneenko V. P., Metody optimizatsii, Uchebnik, Vyssh. shk., M., 2007, 664 pp.

[18] Muller D. E., “A method for solving algebraic equations using an automatic computer”, Mathematical Tables and Other Aids to Computation, 10:56 (1956), 208–215 | DOI | MR | Zbl

[19] Samarskii A. A., Gulin A. V., Chislennye metody, Ucheb. posobie dlya vuzov, Nauka, M., 1989, 432 pp.

[20] Elenin G. G., Shlyakhov P. I., “Geometricheskaya struktura prostranstva parametrov trekhstadii-nykh simplekticheskikh metodov Runge–Kutty”, Matematicheskoe modelirovanie, 23:5 (2011), 16–34 | MR | Zbl

[21] Hunter J. D., “Matplotlib: A 2D graphics environment”, Computing in Science Engineering, 9:3 (2007), 90–95 | DOI