I have a 3D space (x, y, z) with an additional parameter at each point (energy), giving 4 dimensions of data in total.
I would like to find a set of x, y, z points w
Why not try quadlinear interpolation?
extend Trilinear interpolation by another dimension. As long as a linear interpolation model fits your data, it should work.