Física Estadística II – 2024

Unidad 2 - Parte 2

2. Estadística de Fermi-Dirac y Bose-Einstein. Límite de Maxwell-Boltzmann. Energía de

Fermi. Gas ideal de Fermi - Dirac; electrones en metales. Calor específico electrónico. El gas
ideal de bosones. La condensación de Bose-Einstein.

En esta parte veremos: El gas ideal de bosones. La condensación de Bose-Einstein.

Ahora resolvamos el gas de bosones, que tiene una particularidad !!

Comencemos con algo que está en la clase de Zanette que me resultó muy interesante e
ilustrativo, y es un ejemplo de un sistema de bosones con solamente 2 estados de energía
disponibles (primera parte de la Clase 16,
https://www.youtube.com/watch?v=eElzO2Swhd0&list=PLv0-vSkgrxj-
5sd5SmNX3h4ZgqGw5xmTl&index=16).

Supongamos un sistema de bosones con dos niveles de energía:

o = 0
1 = 
Como son bosones:
1
1 1 1
N=∑ = +
eβ (εl − μ) − 1 e−β μ − 1 eβ (ε− μ) − 1
𝑙=0

Escribiendo esto en función de z, en vez de , me queda

𝑧 z
N= +
1 − z eβ ε − z
En el límite de N  , como

eβ ε > 1 y z<1
para que z sea una solución en la ecuación de N, entonces z debe tender a 1, de manera que el
primer término de la derecha tiende a , como lo hace N.

Pueden probar con MatLab o Python de resolver la ecuación y les da dos soluciones:
1
z=(N + eβ ε + N*eβ ε  (N2*e2β ε - 2*N2*eβ ε + N2 + 2*N*e2β ε - 4*N*eβ ε + 2*N +
e2β ε + 2*eβ ε + 1)^(1/2) + 1)/(2*(N + 2))
Si hacen tender N a infinito, queda z  1.
En realidad, trabajando un poco más (ver video del Dr. Zanette):

1
z1-N
Entonces:

No  N
y

1
N  eβ ε −1

Esto quiere decir que a medida que agrego partículas, casi todas van al estado fundamental
independientemente de la temperatura !!!
Entonces, a diferencia de los fermiones, cuando agrego bosones a un sistema de muy muchos
bosones (N muy muy grande), el estado de energía  se satura, y todo lo que agrego va a
parar al estado fundamental !!!

Qué diferente al comportamiento de los fermiones !!!

Retomemos:
Volvamos al gas donde, en vez de 2 estados tengo infinitos estados.
Empecemos con N.

1 z
N=∑ =∑
eβ (εl − μ) − 1 eβ εl − z
l l

Al pasar la sumatoria a integral, en este caso el término que corresponde a l=0 diverge porque
z1 y exp(0)=1.

Además, la densidad de estados que voy a utilizar que es la que corresponde a partículas con
energía cinética en una caja de 3 dimensiones:

2
 3
dϕ π 8m 1 2
ω(ε) dε = dε = ( 2 ) V ε2 dε
dε 4 h
en función de p queda:

3 1
2 2
π 8m p p2 4πV
( 2 ) V ( ) ( ) dp = 3 p2 dp
4 h 2m m h

4πV ∞ 2
∑→ ∫ p dp
h3 0
l

Es decir que para p=0, que sería equivalente a l=0 (estado fundamental), el peso que le da la
densidad de estados es 0 !!!.
Entonces lo tengo que “aislar” antes de pasar a una integral 

z z z
N=∑ β pl 2
= + ∑ β pl 2
1−z
l e 2m −z l (sin l=0) e 2m −z
y ahora si paso a integral.

Si no lo hiciera, dan mal los resultados. De esto se dio cuenta Einstein, quien hizo esta
corrección. Por eso creo que se llama estadísitica de Bose-Einstein.
Hacer esto, equivale a modificar la densidad de estados:

3
Como el término de l0 en la integral es muy pequeño, puedo igual integrar desde cero, a
pesar de haberlo puesto ya explícitamente fuera de la integral:

z 4πV ∞ 2 z
N= + 3 ∫ p β p2
dp
1−z h 0
e 2m −z
Algunos libros tienen en cuenta este “cachitín” de más porque a la integral la puse desde
cero, y para ser exactos deberíamos escribir, por ejemplo, para N:

z 4πV ∞ 2 z
N= + 3 ∫ p β p2
dp − IN (z, … )
1−z h 0
e 2m −z
que es muy pequeño y tiende a cero para z1

Hacemos lo mismo ahora para el potencial Grandcanónico donde teníamos que:

BE = - kT ln ZGBE = kT ∑l ln[1 − e−β (εl − μ) ]

O sea:

BE = kT ∑ ln[1 − z e−β εl ]

l

Aquí también para el nivel fundamental, es decir l=0, tenemos 0=0, y como para esta
situación z1 tendríamos ln(1-1e0) o sea ln(0), que diverge y la integral se comería a este
término como 0 porque la densidad de estados que estoy usando vale 0 en l=0. Entonces, lo
“aíslo” de nuevo antes de pasar a integral.
Me queda así:

BE = kT ln(1 − z) + kT ∑ ln[1 − z e−β εl ]

l (sin l=0)

y haciendo el cambio de variable de  a p, y de suma a integral:

4
 4πkT V ∞ 2 −
β p2
BE = kT ln(1 − z) + 3
∫ p ln [1 − z e 2m ] dp
h 0

Hacemos como en FD el cambio de variable a x2=p2/2m y nos queda:

∞
4kT V 2
BE = kT ln(1 − z) + ∫ x 2 ln[1 − z e−x ] dx
Λ3 √ 𝜋 0

o sea para P que es -/V:

∞
kT 4kT 2
P = − ln(1 − z) − 3 ∫ x 2 ln[1 − z e−x ] dx
𝑉 Λ √𝜋 0
(Aquí también estoy cometiendo un pequeño error al integrar desde cero que algunos libros
lo consideran restando un término IP)

y para N queda:
∞
z 4V 2
z
N= + ∫ x x2 dx
1 − z Λ3 √ π 0 e −z

Usando las funciones g5/2 y g3/2

∞ ∞
4 2 −x2
zα
g 5 (z) = − ∫ x ln(1 − z e ) dx = ∑ 5/2
2 √π 0 α
α=1

∞ ∞
d 4 2
z zα
g 3 (z) = z g 5 (z) = ∫ x x2 dx = ∑ 3/2
2 dz 2 √π 0 e − z α
α=1

y además para z1 tienden a las funciones zeta de Riemann ():

5
 5 3
g 5 (1) = ζ ( ) ≈ 1.342 … y g 3 (1) = ζ ( ) ≈ 2.612 …
2 2 2 2

Como para amigarnos un poco, aquí están las gráficas de estas funciones en el rango de su
dominio para nuestro análisis (recordemos que z está acotado entre 0 y 1 para los bosones):

Entonces nos queda:

kT kT
P = − ln(1 − z) + 3 g 5 (z)
V Λ 2

z V
N= + 3 g 3 (z)
1−z Λ 2

O sea que para la densidad n (=N/V) queda:

1 z 1
n= + 3 g 3 (z)
V1 − z Λ 2

Veamos que ocurre con todo esto en el límite V  , que sería para nosotros el límite
termodinámico.

6
Examinemos las funciones que aparecen como primer término en P y n, es decir:

kT
− ln(1 − z) en el caso de P
V
1 z
en el caso de n
V1− z
Si fijamos la densidad de nuestro gas de bosones, es decir n=constante, y fijamos también T,
y agrandamos el volumen (V y z1), tenemos:

kT
lim (− ln(1 − z)) = 0
𝑉→∞ V

1 z
lim ( ) = 𝑛0
𝑉→∞ V 1 − z

z
(Recordemos que es la cantidad de partículas en el estado fundamental, entonces dividido
1−z
en V, da la densidad en el estado fundamental: n0. Ojo aquí porque veníamos usando n0
precisamente para la cantidad (absoluta) de partículas en l=0 ¡!!)

(z en realidad es función de V: z(V) tal que z(V)1. Se podría pensar que si n es fijo y
agrando mucho el volumen, el recipiente está como vacío así que no debería costarme nada
agregar una partícula al gas ¡!!! Les parece razonable ¿?)

Es decir que en el límite de V podemos despreciar la contribución a la presión del primer
término pero en el caso de la densidad no ¡!
1 z
Veamos cómo se comporta en las cercanías de z=1 a medida que agrandamos V:
V 1−z

7
 1 z
Vemos que para z1 el factor siempre se va a .
V 1−z
1 z
A medida que agrandamos V, z necesita estar más cerca de 1 para que sea significativo.
V 1−z
Si no es casi cero ¡!! Hasta que de repente se va a  ¡!!

Ahora, veamos el límite V y N de manera que la densidad n=N/V permanezca finita.
Me queda para la presión:

kT
g 5 (z) si z < 1
Λ3 2
P= y
kT
3
g 5 (1) si z = 1
{Λ 2

1
3
g 3 (z) si z < 1
𝑁 Λ 2
n= = y
𝑉 1
𝑛 + g 3 (1) si z = 1
{ 0 Λ3 2

Aclaremos esto con un gráfico (T  ), donde en vez de la densidad n, se grafica n3, es
decir:
8
 𝑛 Λ3 = 𝑛0 Λ3 + g 3 (𝑧)
2
o
𝑛 λ 𝑇 3 = 𝑛0 λ 𝑇 3 + g 3 (𝑧)
2

En general, y lejos de z=1, el mayor aporte a la densidad de partículas proviene de g 3 (𝑧).

2
A medida que z tiende a 1 g 3 (𝑧) tiende a una constante: g 3 (1) = 2.612
2 2

Entonces ya este término no puede aportar al crecimiento de n y todo debe venir de 𝑛0 λ 𝑇 3 .

Esto se denomina “condensación de Bose-Einstein”.
En otras palabras: a medida que z1 se da una condensación macroscópica en el estado
fundamental del gas.
O: la densidad numérica de partículas en el estado fundamental 0, es decir n0, se vuelve una
cantidad macroscópica ¡!!

En el límite termodinámico la condensación de Bose-Einstein comienza cuando z=1 (o, lo

que es lo mismo cuando =0). Es decir, cuando

𝑛 λ 𝑇 3 = g 3 (1) ≈ 2.612
2
Este punto se denomina punto crítico que es justo cuando comienza la condensación.
Entonces los valores de las variables termodinámicas justo en este punto reciben el nombre
de variables “críticas”.

9
Para la densidad crítica, nc, en función de T tenemos entonces:

g 3 (1) 3
1 √2π m kT
nc = = 2 3 ≈ 2.612 ( )
vc λT h

donde vc es el volumen crítico por partícula.

La temperatura crítica Tc (es decir la temperatura a la cual se da la condensación de BE en

función de n) se obtiene de

g 3 (1)
λT𝑐 3 = 2
𝑛
es decir:
2
3 2
2 2
h 𝑛 h 𝑛 3
Tc = ( ) ≈ ( )
2π m k g 3 (1) 2π m k 2.612
2
El parámetro de orden (que veremos más adelante cuando veamos fenómenos críticos, pero
como para que tengan una idea es un parámetro que se hace cero para T mayores a Tc, o para
n mayor a nc) se define como

g3 (1) 3
𝑛𝑜 2 𝑇 2
𝜂= =1− =1− ( )
𝑛
𝑛 λT
3 𝑇𝑐

La gráfica de este parámetro queda en función de T/Tc queda:

10
La región por debajo de la curva muestra la región de co-existencia entre las dos fases: la
“normal” y la “condensada”

Veamos ahora P. No entraremos en detalle, simplemente veamos las isotermas en la figura de

P vs. v (el volumen por partícula)

Para la entropía tenemos:

Y para el calor específico a densidad constante (que si suponemos que N casi no varía es
como CV), se tiene:
11
La transición de fase del gas de Bose-Einstein es el resultado de la estadística.
El gas de Fermi-Dirac no presenta esta transición ¡!!

En el límite de T altas tenemos que

z0
 𝑔5 (𝑧) ≈ 𝑔3 (𝑧) ≈ 𝑔1 (𝑧) ≈ 𝑧
2 2 2

 Para n promedio que teníamos dos expresiones, tendríamos la de z<1

3
1 z √2𝜋𝑚𝑘𝑇 𝜇
n = 3 g 3 (z) → 3 = ( ) 𝑒 𝑘𝑇
Λ 2 Λ ℎ

Y lo mismo para P. Parto de la expresión para z<1:

kT kT
𝑃 = 3 g 5 (z) → 3 z
Λ 2 Λ
Y ohhhhh ¡!!! z/3 es n promedio 

12
 kT 𝑁𝑘𝑇
𝑃= z = nkT =
Λ3 𝑉

Para Cn teníamos en el caso de z<1:

g3 (z)
15 1 9 2
𝐶𝑛 = 𝑘 g (z) − 𝑛 𝑘
3 5
4 Λ 2 4 g1 (z)
2
15 1 9 z 15 9 3 3 𝑁𝑘
𝐶𝑛 → 𝑘 z−𝑛𝑘 = 𝑛𝑘 − 𝑛𝑘 = 𝑛𝑘 =
4 Λ 3 4 z 4 4 2 2 𝑉

Conclusión: el gas de Bose Einstein se comporta como el gas clásico en el límite de T alta.

-----

Comparemos ahora las presiones (en realidad su variación con v por ejemplo) para los tres
tipos de gases: FD, BE y clásico (Cl o MB).

v: volumen por partícula (en realidad es volumen promedio por partícula N/V, y N es
promedio ¡!!! V no es promedio. Es una constante. Recuerden que vemos recipientes de
volumen fijo ¡!!!

P para BE es muchísimo menor que P de un gas de fermiones o uno clásico ¡!!

13
Por qué ¿????
Porque cuando estamos a v<vc hay una avalancha de partículas hacia el estado fundamental y
por lo tanto esas partículas NO contribuyen a P ¡!!!!

P para FD es siempre mayor que en BE y en el gas clásico.

Esto ocurre porque solo una partícula puede estar en el estado fundamental mientras que en el
clásico no hay límite ¡!! (en realidad ya estamos “violando” aquí los supuestos que hicimos
para poder usar tranquilos MB ¡!! ojo ¡!!!)

-------
Volvamos un poco la distribución de cada estadística:

1
𝐵𝐸
eβ (εl − μ) − 1

𝑛𝑙 = e−β (εl − μ) 𝑀𝐵

1
𝐹𝐷
{ eβ (εl − μ) + 1

14
Volvamos al número promedio de ocupación por microestado de energía de una partícula
para fermiones (estadística de FD):

1
𝑛(ε) 𝑜 𝑛ε =
eβ (ε− μ) + 1

Para cualquier , y una T dada, veamos que pasa para distintos :

* si -  n  0 porque eβ (ε− μ)  
 N/V  0

* si    n  1 porque eβ (ε− μ)  0
 N/V   que en realidad no es así porque tengo un número finito de estados 

1
𝑛ε =
eβ (ε− μ) + 1
Ahora para cualquier , y una  dada, veamos que pasa para distintas T:

* si T 0 y >  n  0 porque eβ (ε− μ)  

* si T 0 y <  n  1 porque eβ (ε− μ)  0
15
Es decir que a 0 K (0 absoluto) las partículas se acomodan cada una en un microestado
comenzando desde el estado fundamental (=0) y así siguiendo hasta acomodarse las N que
haya. El último  en llenarse es =.
Todos los microestados con > quedarán vacíos
El hecho de que solo una partícula pueda estar en un microestado hace que la energía del
punto 0 sea grande.
Para T mayor que 0K la función n deja de ser una función escalón y comienzan a ocuparse
los niveles con > por las partículas que estaban por debajo y comienzan a excitarse con el
kT que recibe el sistema cuando aumento T.

Tanto en el caso de BE como en FD cuando el número de estados disponibles excede al N,

ambas tienden a MB ¡!

Esto implica n  0 para todo 

En otras palabras, casi todos los microestados están vacíos excepto unos cuantos que tiene
n=1.

Esta condición se cumple cuando N/V  0 (a esto ya los sabíamos de gas clásico: debe tener
densidad baja cuando tengo T fija, o para densidad baja si agrando T puedo llegar al límite
cásico también).

Densidad baja implica z 0, por lo siguiente:

1 z
𝑛ε = =
eβ (ε− μ) + 1 eβ ε + z
Ahora, para z chico:
𝑛ε → z e−β ε

porque en el denominador de n z se hace despreciable frente al exponencial ¡!

Y ya sabemos del gas clásico que N= eβ μ Z (donde Z es la función de partición canónica de
una única partícula en el gas) o sea

𝑛ε e−β ε
=
𝑁 𝑍
16
𝑛ε
es la probabilidad de que una partícula del gas esté en el microestado  para un gas con
𝑁
distribución de MB

Cómo se modifica la curva de n para fermiones si en vez de microestados la hacemos en

función de niveles ?
En el caso de partículas libres en una caja quedaría:

Les parece bien ???

Ahora BE:
1
𝑛ε = ε− μ
e 𝑘𝑇 −1
En el caso de BE para T0 , si  es menor que el estado fundamental, entonces - es
siempre finito y tendrá n0, es decir que N promedio es también cero.

Si  es casi del valor del estado fundamental, entonces para ese estado tendré n y el
resto vacío (tendría daría un 0/0 en la exponencial, que es como si tuviera un 0, exp(0)=1 y m
queda un 1/0 en n)

17
En otras palabras, en el límite de T0 el único valor termodinámicamente significativo de 
es el que es igual a la energía del estado fundamental, que en general se considera como 0.
En este límite todas las partículas están en estado fundamental.
Aquí también para z0

𝑛ε → z e−β ε

Como en FD ¡!!!

-----------
Vamos ahora a un ejemplo de BE:

En el caso del gas de electrones los efectos cuánticos se ven a T ambiente gracias a que la
masa del electrón es muy pequeña (además de la gran densidad del gas que vimos como
ejemplo).

Dado que en el caso de los bosones la T a la cual ocurre la condensación de Bose está dada
por:
2
h2 𝑛 3
Tc ≈ ( )
2π m k 2.612
convendría utilizar el elemento más liviano posible para no tener que bajar tanto T. Pero con
el Hidrógeno (que es el elemento de menor masa atómica de la Tabla Periódica) no se
puede !!!

El experimento que demostró por primera vez la existencia de la condensación de BE se

realizó en 1995 (o por lo menos es la fecha del paper) con átomos de Rubidio(87): ~2000
átomos y una Tc de ~200 nK !!!

Veamos el caso del Helio 4: 4He, que es otro elemento liviano.

La masa del 4He, si bien es pequeña, es mucho más grande que la de los electrones, entonces
los efectos cuánticos se empiezan a notar a T muy bajas ¡!!!
Para el Helio
V
= 0.67 recien para el Helio líquido a T = 4.2 K
N Λ3

18
En el gas de Helio las fuerzas intermoleculares son débiles aún a T baja, por lo tanto se hacen
visibles los efectos cuánticos que acabamos de ver sin que la “interacción” entorpezca
nuestra teoría.
Recordemos que todo lo anterior lo deducimos para gas ideal: interacción despreciable ¡!!!
Sin embargo, las propiedades tan interesantes que tiene el Helio cuando se produce la
condensación de Bose resultan de una mezcla de la condensación misma + los efectos de la
interacción. En realidad la superfluidez no es consecuencia de la condensación de Bose.
Algunos aspectos se parecen a los de un condensado de Bose, y por eso se dice que se trata
de una condensación de Bose. Pero muchos otros no.
Cuales son estos aspectos que se parecen a una condensación de Bose ?
De acuerdo a lo que vimos, calculemos a que T se produciría la condensación de Bose-
Einstein en un gas de Helio.

2
h2 𝑛 3
Tc = ( )
2π m k 2.612
El Helio líquido tiene densidad 27.6 cm3/mol
m = masa de un átomo de Helio 4 = 6.7× 10-27 kg (lo puedo sacar a partir de la masa atómica
o como ~2 protones + 2 neutrones)
h = constante de Planck = 6.6× 10-34 m2 kg s-1
k = constante de Boltzmann = 1.38 × 10-23 m2 kg s-2 K-1
N/V = n = 0.036 mol/cm3

Esto da TC = 3.13K que es próxima a la TC de superfluidez: TC= 2.19K. O al menos se

parecen !!

El Cv se comporta también medio parecido al de la condensación de Bose:

19
Se lo denomina transición  (lambda) por la forma de la curva de CV.

Algo importante de notar:

Esta condensación es una condensación en el espacio de las fases (no en el espacio de
coordenadas): las partículas pasan a tener energía y momento cero.

20