Density expansion based on Hermite functions