Is there any way to incorporate the intraband Drude-like term in optical spectrum computation? ATK does it via Kubo-Greenwood formula, where there should be one (-1/w^2) divergence at zero energy for real part of epsilon. I have tested it for all sorts of broadenings for simple metals like Au or Ag. However, the zero energy negative divergence is not produced. Can anybody help?