$${\rho_f=\frac{A_{21}}{B_{12} \exp\big(\frac{h f}{k_B T}\big)-B_{21}}}$$