Introduction

This program calculates the entropy of two Einstein solids that can exchange energy. The plot shows the entropy
of the two systems, S_{A} (blue curve) and S_{B} (green curve), and the entropy S (red curve) of the composite system. N_{A} and N_{B} are the number of particles in each system. The total energy is given by E = E_{A} + E_{B}.

The entropy is computed by using the binomial coefficients to determine the number of states as a function of energy. In each
case the entropy is plotted versus E_{A}, the energy of system A. The entropy of system A is given by

S(E_{A}) = k ln Ω(E_{A}),

where Ω_{A}(E_{A}) is the number of microstates of system A with energy E_{A} and k is Boltzmann's constant. We will choose units such that k = 1. The entropy of system B is given by a similar expression. The entropy of the composite system is given by

S(E_{A}) = k ln [Ω_{A}(E_{A}) Ω_{B}(E_{B})] = S_{A} + S_{B}.

The total entropy of the system is given by

.

We will see that as the total number of particles increases, the total entropy can be approximated by

,

where is the most probable value of E_{A}; that is, the value of E_{A} for which S(E_{A}) is a maximum.

Problems

- Choose various values of the parameters and describe the qualitative behavior of S
_{A}, S_{B}, and S as a function of E_{A} - The message box gives the relative error in the total entropy when it is approximated by its maximum value as a function of E
_{A}. Describe the behavior of this error as you increase the size of the system and keep the energy per particle a constant. - Discuss why the total entropy has a maximum. What can you say about the slopes of S
_{A}and S_{B}at the value of E_{A}where S(E_{A}) is a maximum? - The values of the inverse slopes dS
_{A}/dE_{A}and dS_{B}/dE_{B}can be obtained by clicking your mouse on the corresponding curves. The values of the (inverse) slopes are given in the lower right corner. We know that the two systems are in equilibrium when their temperatures are equal. Discuss the relation of the slopes to the temperature.

Java Classes

- EntropyApp

Updated 28 December 2009.