Efficient geometric integrators for nonadiabatic quantum dynamics. I. The adiabatic representation
Fuente:
arXiv
Saved in:
| Main Authors: | , |
|---|---|
| Format: | Preprint |
| Published: |
2019
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866910618731151360 |
|---|---|
| author | Choi, Seonghoon Vaníček, Jiří |
| author_facet | Choi, Seonghoon Vaníček, Jiří |
| contents | Geometric integrators of the Schrödinger equation conserve exactly many invariants of the exact solution. Among these integrators, the split-operator algorithm is explicit and easy to implement, but, unfortunately, is restricted to systems whose Hamiltonian is separable into a kinetic and potential terms. Here, we describe several implicit geometric integrators applicable to both separable and non-separable Hamiltonians, and, in particular, to the nonadiabatic molecular Hamiltonian in the adiabatic representation. These integrators combine the dynamic Fourier method with recursive symmetric composition of the trapezoidal rule or implicit midpoint method, which results in an arbitrary order of accuracy in the time step. Moreover, these integrators are exactly unitary, symplectic, symmetric, time-reversible, and stable, and, in contrast to the split-operator algorithm, conserve energy exactly, regardless of the accuracy of the solution. The order of convergence and conservation of geometric properties are proven analytically and demonstrated numerically on a two-surface NaI model in the adiabatic representation. Although each step of the higher order integrators is more costly, these algorithms become the most efficient ones if higher accuracy is desired; a thousand-fold speedup compared to the second-order trapezoidal rule (the Crank-Nicolson method) was observed for wavefunction convergence error of $10^{-10}$. In a companion paper [J. Roulet, S. Choi, and J. Van\'ıček (2019)], we discuss analogous, arbitrary-order compositions of the split-operator algorithm and apply both types of geometric integrators to a higher-dimensional system in the diabatic representation. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_1902_04661 |
| institution | arXiv |
| publishDate | 2019 |
| record_format | arxiv |
| spellingShingle | Efficient geometric integrators for nonadiabatic quantum dynamics. I. The adiabatic representation Choi, Seonghoon Vaníček, Jiří Chemical Physics Computational Physics Quantum Physics Geometric integrators of the Schrödinger equation conserve exactly many invariants of the exact solution. Among these integrators, the split-operator algorithm is explicit and easy to implement, but, unfortunately, is restricted to systems whose Hamiltonian is separable into a kinetic and potential terms. Here, we describe several implicit geometric integrators applicable to both separable and non-separable Hamiltonians, and, in particular, to the nonadiabatic molecular Hamiltonian in the adiabatic representation. These integrators combine the dynamic Fourier method with recursive symmetric composition of the trapezoidal rule or implicit midpoint method, which results in an arbitrary order of accuracy in the time step. Moreover, these integrators are exactly unitary, symplectic, symmetric, time-reversible, and stable, and, in contrast to the split-operator algorithm, conserve energy exactly, regardless of the accuracy of the solution. The order of convergence and conservation of geometric properties are proven analytically and demonstrated numerically on a two-surface NaI model in the adiabatic representation. Although each step of the higher order integrators is more costly, these algorithms become the most efficient ones if higher accuracy is desired; a thousand-fold speedup compared to the second-order trapezoidal rule (the Crank-Nicolson method) was observed for wavefunction convergence error of $10^{-10}$. In a companion paper [J. Roulet, S. Choi, and J. Van\'ıček (2019)], we discuss analogous, arbitrary-order compositions of the split-operator algorithm and apply both types of geometric integrators to a higher-dimensional system in the diabatic representation. |
| title | Efficient geometric integrators for nonadiabatic quantum dynamics. I. The adiabatic representation |
| topic | Chemical Physics Computational Physics Quantum Physics |
| url | https://arxiv.org/abs/1902.04661 |