El potencial químico se define como: $$ \mu = -T\frac{\partial{S(N,V,E)}}{\partial{N}} $$ A mí me parece que esto es completamente independiente de donde puedo poner el punto de referencia de la energía, debido a la diferencia de entropía es relevante (y también de la temperatura se define como una diferencia en la entropía).
Sin embargo, la de Fermi-Dirac distribución es: $$<n_r> = \frac{1}{exp(\beta(\epsilon_r-\mu))+1}$$
Pero si puedo cambiar el valor del punto de referencia de la energía, el valor de $<n_r>$ cambios, lo que provoca una contradicción. En mi libro acerca de la mecánica estadística, que el estado a pesar de que $\epsilon_r-\mu$ es independiente del punto de referencia de la energía, pero no veo por qué, porque se $\mu$ parece ser independiente, mientras que $\epsilon_r$ no parece independiente del punto de referencia.