That looks like the sum you need to do if you were to calculate the Partition function for a quantum mechanical dumbell. That sum cannot be computed in closed form. The typical approximations are to assume either low temperature ( [itex]C \rightarrow \infty[/itex]), in which case the sum is approximately
[tex]1 + 3e^{-2C} + \ldots .[/tex]
The other limit is high temperature, [itex]C \rightarrow 0[/itex], in which case the sum is approximately an integral,
[tex]\int_0^\infty dn~(2n+1)e^{-Cn(n+1)},[/itex]<br />
<br />
which is easily solved by substitution.<br />
<br />
Systematic corrections to the integral form can be computed using the Euler-Maclaurin formula:<br />
<br />
<a href="http://en.wikipedia.org/wiki/Euler-Maclaurin_formula" target="_blank" class="link link--external" rel="nofollow ugc noopener">http://en.wikipedia.org/wiki/Euler-Maclaurin_formula</a>[/tex]