Yo no estoy seguro de cómo obtuvo la última expresión. El estándar de Sommerfeld de expansión (para más detalles, véase, por ejemplo, Ashcroft Y Mermin) da un poco diferente resultado, que es
EF≈μ[1+π28(kBTμ)2)2/3≈μ[1+π212(kBTμ)2]
al líder de la no-trivial orden en kBT/μ, es decir, O((kBT/μ)2).
(He puesto E0=0 en la expresión de la Wikipedia.)
Podemos invertir esta relación mediante la sustitución de μ=EF[1+c2(kBT/EF)2+⋯] en el de arriba. Es decir,
EF=EF[1+c2(kBT/EF)2+⋯]{1+π212(kBTEF)2[1+c2(kBT/EF)2+⋯)−2}.
Comparando el orden cero de términos en kBT/EF a ambos lados de la ecuación anterior, simplemente obtener EF=EF. La comparación de los términos de segundo orden, tenemos 0=π212+c2. Por lo tanto
μ=EF[1−π212(kBTEF)2]
hasta el O((kBT/EF)2).
Para determinar el siguiente orden de corrección, a μ, se debe incluir un término de orden superior en la Sommerfeld de expansión, lo que da
EF≈μ[1+π28(kBTμ)2+7π4640(kBTμ)4]2/3.
La expansión de este a O((kBT/μ)4), la sustitución de μ=EF[1+c2(kBT/EF)2+c4(kBT/EF)4+⋯] (a los que ya hemos determinado c2 antes) en ella, y luego la coincidencia de los coeficientes de ambos lados hasta el O((kBT/EF)4) conducirá al resultado deseado.