A New 3D Potential-Density Basis Set

A set of bi-orthogonal potential-density basis functions is introduced to model the density and its associated gravitational field of three dimensional stellar systems. Radial components of our basis functions are weighted integral forms of spherical Bessel functions. We discuss the properties of our basis functions and demonstrate their shapes for the latitudinal Fourier number $l=2$.