Saved in:
| Main Authors: | , |
|---|---|
| Format: | Preprint |
| Published: |
2023
|
| Subjects: | |
| Online Access: | https://arxiv.org/abs/2310.05633 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866910618849640448 |
|---|---|
| author | Fereidani, Roya Moghaddasi Vaníček, Jiří JL |
| author_facet | Fereidani, Roya Moghaddasi Vaníček, Jiří JL |
| contents | Gaussian wavepacket dynamics has proven to be a useful semiclassical approximation for quantum simulations of high-dimensional systems with low anharmonicity. Compared to Heller's original local harmonic method, the variational Gaussian wavepacket dynamics is more accurate, but much more difficult to apply in practice because it requires evaluating the expectation values of the potential energy, gradient, and Hessian. If the variational approach is applied to the local cubic approximation of the potential, these expectation values can be evaluated analytically, but still require the costly third derivative of the potential. To reduce the cost of the resulting local cubic variational Gaussian wavepacket dynamics, we describe efficient high-order geometric integrators, which are symplectic, time-reversible, and norm-conserving. For small time steps, they also conserve the effective energy. We demonstrate the efficiency and geometric properties of these integrators numerically on a multi-dimensional, nonseparable coupled Morse potential. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2310_05633 |
| institution | arXiv |
| publishDate | 2023 |
| record_format | arxiv |
| spellingShingle | High-order geometric integrators for the local cubic variational Gaussian wavepacket dynamics Fereidani, Roya Moghaddasi Vaníček, Jiří JL Numerical Analysis Mathematical Physics Quantum Physics Gaussian wavepacket dynamics has proven to be a useful semiclassical approximation for quantum simulations of high-dimensional systems with low anharmonicity. Compared to Heller's original local harmonic method, the variational Gaussian wavepacket dynamics is more accurate, but much more difficult to apply in practice because it requires evaluating the expectation values of the potential energy, gradient, and Hessian. If the variational approach is applied to the local cubic approximation of the potential, these expectation values can be evaluated analytically, but still require the costly third derivative of the potential. To reduce the cost of the resulting local cubic variational Gaussian wavepacket dynamics, we describe efficient high-order geometric integrators, which are symplectic, time-reversible, and norm-conserving. For small time steps, they also conserve the effective energy. We demonstrate the efficiency and geometric properties of these integrators numerically on a multi-dimensional, nonseparable coupled Morse potential. |
| title | High-order geometric integrators for the local cubic variational Gaussian wavepacket dynamics |
| topic | Numerical Analysis Mathematical Physics Quantum Physics |
| url | https://arxiv.org/abs/2310.05633 |