?
GPU-Accelerated Matrix Exponent for Solving 1D Time-Dependent Schrödinger Equation
Non-adiabatic electron-ion quantum dynamics is still an area of many unresolved problems even for such simple systems as the H2+ molecular ion. Mathematical modelling based on time-dependent Schrödinger equation (TDSE) is an important method that can provide better understanding of these phenomena. In this work, we present TDSE solution for 1D TDSE that describes non-adiabatic electron-ion quantum dynamics for the simplified H2+ model. For solving TDSE, we use the real-space representation and the matrix exponent method that is quite computationally expensive but is free from usual symmetry-based simplifications. For this purpose, we make use of the very high performance of modern Nvidia V100/A100 GPUs and deploy our parallel multi-GPU matrix multiplication algorithm.