High-accuracy numerical calculations of the bound states of a hydrogen atom in a constant magnetic field with arbitrary strength

We develop a simple and effective method for solving the Schrödinger equation of a hydrogen atom in a constant magnetic field with arbitrary strength. Energies are obtained not only for the ground and low-lying states but also for highly excited states with precision from 12 up to 20 decimal digits. The calculations are performed for an entire range of magnetic field intensity up to 9.4 x 10^8 Tesla, the strongest field ever observed. The strong point of the development of the method is the construction of an anharmonic oscillator model for a hydrogen atom in a constant magnetic field via the Kustaanheimo–Stiefel transformation. This model allows the use of purely algebraic calculations and the Feranchuk–Komarov (FK) operator method for effectively solving the Schrödinger equation. The advantages of the basis set in this work are also discussed to extend its application to other problems, such as multi-electron atoms in a constant magnetic field. We also provide a program written by FORTRAN for the solutions mentioned above.