According to the second law of thermodynamics, a system assumes a configuration of maximum entropy at thermodynamic equilibrium. We seek a probability distribution of states ρ ρ i that maximizes the discrete Gibbs entropy S = − − k B ∑ ∑ i ρ ρ i ln ⁡ ⁡ ρ ρ i subject to two physical constraints: