Thermodynamic functions computed from a partition function

x axis: k_BT/ε. Two-level system: ● excited population, ● entropy S/(Nk ln 2), ● heat capacity C/(Nk) with its Schottky maximum. Einstein solid: ● C/(3Nk), ● mean energy per oscillator in units of ε.