Optimal search tree algorithm for points inside spherical regions

Question

I have a large number of 3d spheres, each with some individual radius. I need to find, for each sphere S, all other spheres whose center is within the radius of S. Right now