Krylov integrators for Hamiltonian systems
Permalink
Citation
Eirola , T & Koskela , A 2019 , ' Krylov integrators for Hamiltonian systems ' , BIT Numerical Mathematics , vol. 59 , no. 1 , pp. 57–76 . https://doi.org/10.1007/s10543-018-0732-y
| Title: | Krylov integrators for Hamiltonian systems |
| Author: | Eirola, Timo; Koskela, Antti |
| Contributor organization: | Department of Mathematics and Statistics |
| Date: | 2019-03 |
| Language: | eng |
| Number of pages: | 20 |
| Belongs to series: | BIT Numerical Mathematics |
| ISSN: | 0006-3835 |
| DOI: | https://doi.org/10.1007/s10543-018-0732-y |
| URI: | http://hdl.handle.net/10138/265089 |
| Abstract: | We consider Arnoldi-like processes to obtain symplectic subspaces for Hamiltonian systems. Large dimensional systems are locally approximated by ones living in low dimensional subspaces, and we especially consider Krylov subspaces and some of their extensions. These subspaces can be utilized in two ways: by solving numerically local small dimensional systems and then mapping back to the large dimension, or by using them for the approximation of necessary functions in exponential integrators applied to large dimensional systems. In the former case one can expect an excellent energy preservation and in the latter this is so for linear systems. We consider second order exponential integrators which solve linear systems exactly and for which these two approaches are in a certain sense equivalent. We also consider the time symmetry preservation properties of the integrators. In numerical experiments these methods combined with symplectic subspaces show promising behavior also when applied to nonlinear Hamiltonian problems. |
| Subject: |
111 Mathematics
Hamiltonian systems Exponential integrators Krylov subspace method Symplectic integrators Hamiltonian Lanczos algorithm EXPONENTIAL INTEGRATORS EQUATIONS Hamiltonian systems Exponential integrators Krylov subspace methods Symplectic integrators Symmetric integrators Hamiltonian Lanczos algorithm EXPONENTIAL INTEGRATORS EQUATIONS |
| Peer reviewed: | Yes |
| Rights: | unspecified |
| Usage restriction: | openAccess |
| Self-archived version: | publishedVersion |
Files in this item
Total number of downloads: Loading...
| Files | Size | Format | View |
|---|---|---|---|
| Eirola_Koskela2 ... tegratorsForHamiltonia.pdf | 1.042Mb |
View/ |
This item appears in the following Collection(s)
-
Articles from TUHAT CRIS [46485]
Related items
Showing items related by title, author, creator and subject.
-
(Helsingin yliopisto, 2019)The spin-orbit-coupled insulator Sr 3 NiIrO 6 is a strongly correlated transition metal compound, where an interplay of geometric frustration and spin anisotropy gives rise to novel magnetic phases. Resonant inelastic x-ray scattering (RIXS) is a powerful probe of the low-lying quasi-particle excitations that underpin these emergent properties. In this work, we partition the active space into approximately non-interacting parts in order to introduce a tight-binding single-particle model Hamiltonian describing the distorted IrO6 octahedra in Sr3NiIrO6. We then use this model to calculate its RIXS spectrum at the Ir L3-edge in the sub-electronvolt range. The results of this calculation are compared with experiments performed at the European Synchrotron Radiation Facility, and with a multiplet crystal field model calculation. We find that this one electron model largely agrees with the full-multiplet model and describes the d-d excitations observed in experiment. The addition of an exchange field term explains the low-lying temperature-dependent magnetic feature, disambiguating the sign of the crystal-field term, and suggesting that the feature is well localized at low temperatures, and is best described as an orbitally- entangled local spin-flip excitation. However, the correspondence at room temperature diminishes, suggesting that dispersive description is necessary to model this regime. The drastic reduction in active space entailed by this model facilitates the creation of extended non-collinear Heisenberg-like models, which can be calculated at a lower computational cost than full multiplet extended models.
-
(2020)We investigate the properties of the excited spectra of the even-even isotopes of krypton using a Generalised Bohr Hamiltonian with three different Skyrme functionals. In particular, we investigate the evolution of the low-lying states and their associated electromagnetic transitions. The model reproduces quite nicely the energy trends apart from 88Kr, where none of the interactions used here are able to grasp a sudden change in the energy spectrum. Additionally, we explore the neutron deficient region 72-76Kr which is a proposed region for shape coexistence. We observe that the model can reproduce the structure of the experimental spectrum of 72Kr exceedingly well. © Published under licence by IOP Publishing Ltd.
-
(Helsingin yliopisto, 2021)HMC is a computational method build to efficiently sample from a high dimensional distribution. Sampling from a distribution is typically a statistical problem and hence a lot of works concerning Hamiltonian Monte Carlo are written in the mathematical language of probability theory, which perhaps is not ideally suited for HMC, since HMC is at its core differential geometry. The purpose of this text is to present the differential geometric tool's needed in HMC and then methodically build the algorithm itself. Since there is a great introductory book to smooth manifolds by Lee and not wanting to completely copy Lee's work from his book, some basic knowledge of differential geometry is left for the reader. Similarly, the author being more comfortable with notions of differential geometry, and to cut down the length of this text, most theorems connected to measure and probability theory are omitted from this work. The first chapter is an introductory chapter that goes through the bare minimum of measure theory needed to motivate Hamiltonian Monte Carlo. Bulk of this text is in the second and third chapter. The second chapter presents the concepts of differential geometry needed to understand the abstract build of Hamiltonian Monte Carlo. Those familiar with differential geometry can possibly skip the second chapter, even though it might be worth while to at least flip through it to fill in on the notations used in this text. The third chapter is the core of this text. There the algorithm is methodically built using the groundwork laid in previous chapters. The most important part and the theoretical heart of the algorithm is presented here in the sections discussing the lift of the target measure. The fourth chapter provides brief practical insight to implementing HMC and also discusses quickly how HMC is currently being improved.