On a General Method for Resolving Integrals of Multiple Spherical Bessel Functions Against Power Laws into Distributions

We here present a method of performing integrals of products of spherical Bessel functions (SBFs) weighted by a power-law. Our method, which begins with double-SBF integrals, exploits a differential operator $\hat{D}$ defined via Bessel's differential equation. Application of this operator raises the power-law in steps of two. We also here display a suitable base integral expression to which this operator can be applied for both even and odd cases. We test our method by showing that it reproduces previously-known solutions. Importantly, it also goes beyond them, offering solutions in terms of singular distributions, Heaviside functions, and Gauss's hypergeometric,$\;_2{\rm F}_1$ for $all$ double-SBF integrals with positive semi-definite integer power-law weight. We then show how our method for double-SBF integrals enables evaluating $arbitrary$ triple-SBF overlap integrals, going beyond the cases currently in the literature. This in turn enables reduction of arbitrary quadruple, quintuple, and sextuple-SBF integrals and beyond into tractable forms.