A numerical algorithm for solving the coupled Schrödinger equations using inverse power method
Abstract
The inverse power method is a numerical algorithm to obtain the eigenvectors of a matrix. In this work, we develop an iteration algorithm, based on the inverse power method, to numerically solve the Schrödinger equation that couples an arbitrary number of components. Such an algorithm can also be applied to the multi-body systems. To show the power and accuracy of this method, we also present an example of solving the Dirac equation under the presence of an external scalar potential and a constant magnetic field, with source code publicly available.
- Publication:
-
Computer Physics Communications
- Pub Date:
- October 2024
- DOI:
- 10.1016/j.cpc.2024.109284
- arXiv:
- arXiv:2403.02747
- Bibcode:
- 2024CoPhC.30309284Z
- Keywords:
-
- Coupled Schrödinger-like equations;
- Eigenvalue and eigenstate;
- Inverse power method;
- Physics - Computational Physics;
- Nuclear Theory;
- Quantum Physics
- E-Print:
- 8 pages, 2 figures