Surrogate density function