Independent symmetric three-level system
= Independent symmetric three-level system
For $N$ noninteracting units with energies $-\epsilon,0,\epsilon$, put $x=\beta\epsilon$. Then
$$
Z=(1+2\cosh x)^N,
\qquad
E=-\frac{2N\epsilon\sinh x}{1+2\cosh x},
$$
and
$$
S=Nk_B\left[
\log(1+2\cosh x)
-\frac{2x\sinh x}{1+2\cosh x}
\right].
$$