Abstract
We present software for calculating the polarisability and atom-photon scattering rates of hydrogen atoms in S-states (orbital angular momentum l = 0). It can also calculate magic wavelengths between pairs of hydrogen S states. Such calculations are essential to the potential optical trapping of atomic hydrogen. Calculations proceed numerically via implicit summation over a basis of radial Sturmian functions. This naturally includes all contributions from both bound and unbound states. The software was developed with low-lying states in mind but applies to all S states and any hydrogen isotope. Extensions to other hydrogen-like systems or different orbital angular momentum states are possible. It is written in Python and uses NumPy and SciPy extensively. It is available via GitHub.
