Mass of an exponential spherical atmosphere ($M=4\pi {\rho }_0{h}_0(R^2+2R{h}_0+2h_0^2)$)

For a spherical surface of radius $R$ with exterior mass density $\rho (r)={\rho }_0{e}^{-(r-R)/h_0}$, integrate the spherical volume element to obtain $M=4\pi {\rho }_0{\int }_0^{\infty }(R+h{)}^2{e}^{-h/h_0}dh$. The moments ${\int }_0^{\infty }{h}^n{e}^{-h/h_0}dh=n!h_0^{n+1}$ give the displayed expression. Its leading thin-atmosphere approximation is surface area times column mass, $4\pi R^2{\rho }_0{h}_0$; the two positive correction terms account for the growth of shell area with altitude.