Stability in Debye series calculation for light scattering by absorbing particles and bubbles

The Debye-series decomposition is of importance for understanding of light scattering features and for the validity of the geometrical optics approximation to light scattering. The numerical stability and accuracy for calculating light scattering with Debye series is studied and an improved algorithm is proposed in this work. The ratios of the Riccati-Bessel functions and the logarithmic derivatives of the Riccati-Bessel functions are employed and calculated with proper recurrences. Exemplifying results are provided to show the improvement of the algorithm.