On the symplectic integration of Hamiltonian systems

Os sistemas Hamiltonianos formam uma das classes mais importantes de equações diferenciais. Além de constituírem o formalismo central da física clássica, sua aplicação se estende a uma grande variedade de outros campos de estudo. Esses sistemas possuem uma característica notória do ponto de vista da matemática, a saber, que a sua ação sobre seus estados iniciais preserva uma estrutura geométrica conhecida como simpleticidade. Este fato tem importantes consequências sobre as características qualitativas do comportamento do sistema, em especial no longo prazo. Neste trabalho, são estudados métodos numéricos para obter soluções aproximadas para sistemas Hamiltonianos (já que, via de regra, soluções exatas não podem ser encontradas) que preservem a estrutura simplética das equações originais. Para tal, é feita uma revisão da teoria clássica da integração numérica de equações diferenciais, bem como de temas mais recentes como os integradores exponenciais. Além de expor a literatura mais recente sobre integradores simpléticos do tipo Runge-Kutta Exponencial, o trabalho propõe um algoritmo para o cálculo computacionalmente eficientes de integrais envolvendo exponenciais de matrizes, que são centrais para a integração simplética estável de ordem alta.Hamiltonian systems form one of the most important classes of differential equations describing the evolution of physical phenomena. They are the backbone of classical mechanics and their application covers many different areas such as molecular dynamics, hydrodynamics, celestial and statistical mechanics, just to mention a few of them. A noteworthy feature of Hamiltonian systems is that their flow possesses a geometric property -known as symplecticity- which has a major impact on the long-time behavior of the solution. Since in general closed-form solutions can be found only in few particular cases, the construction and analysis of numerical integrators -able to produce discrete approximations that are also symplecticity preserving- is crucial for studying these systems. In this work we present the key ideas about Hamiltonian systems and their theoretical properties. We also review the main numerical methods and techniques to design and analyze symplectic integrators. Special attention is given to the stability and dynamical properties of the methods, as well as their effectiveness for long-term simulations. Finally, we propose an algorithm to improve the computational implementation of the family of exponential-based symplectic integrators recently found in the literature