Say I have two identical (fermionic) non-interacting particles in a 1D harmonic oscillator. I would like to compute the entropy of the system as the temperature $T$ varies, for which I need the partition function of this system, then I would calculate the entropy like
$$ F = -k_B T \log(Z) \Longrightarrow S = -\frac{\partial F}{\partial T} $$
However I am stuck when trying to compute $Z$, can anyone help?
My attempt
The particles are non-interacting, therefore the energy of the system is given by $E_{n, m} = E_n + E_m= \hbar \omega(1+n+m)$ so to compute the (canonical) partition function we need to compute
$$ Z = \sum_{n, m}e^{-\beta E_{n, m}} = e^{-\beta \hbar \omega}\sum_{n, m}\Big(e^{-\beta\hbar \omega}\Big)^n\Big(e^{-\beta\hbar \omega}\Big)^m = e^{-\beta \hbar \omega}\sum_n\Big(e^{-\beta\hbar \omega}\Big)^n\Big(e^{-\beta\hbar \omega}\Big)^n + e^{-\beta \hbar \omega}\sum_{n<m}\Big(e^{-\beta\hbar \omega}\Big)^n\Big(e^{-\beta\hbar \omega}\Big)^m = \frac{e^{-\beta \hbar \omega}}{1-e^{-2\beta \hbar \omega}} + \sum_{n<m}\Big(e^{-\beta\hbar \omega}\Big)^n\Big(e^{-\beta\hbar \omega}\Big)^m $$
at this point I am stuck. How can I compute the second part of $Z$?