I am interested in evaluating the following sum
$I\left(s\right)=\sum_{\left\{ m,n\right\} \neq\left\{ 0,0\right\} }\frac{1}{\left(l^{2}+lm+m^{2}\right)^{s}}.$
This is a Madelung Constants type sum arises in triangular/hexagonal lattice in two dimension.
(For a reference: The corresponding expression for a square lattice, similar sum is given by
$\sum_{\left\{ m,n\right\} \neq\left\{ 0,0\right\} }\frac{1}{\left(l^{2}+m^{2}\right)^{2}}$
This has a closed form answer/expression given in