Adding energy output field for density fitting