Problemas para exámen sobre ensambles microcanónicos y

canónicos

Fı́sica Estadı́stica 2019B

Profesor: Dr. Jorge Emmanuel Sánchez Rodrı́guez

Jorge Octavio Vizcaı́no Guajardo

1 Problemas
1.1 Problema 1
Considere una cadena unidimensional conformada por n(>> 1) elementos. Sea a la dis-
tancia entre cada elemento y x la distancia entre los extremos. Determine la entropı́a de
esta cadena como función de x y obtenga la relación entre la temperatura de la cadena y
la fuerza de tensión que es necesaria para mantenerla una distancia x, asumiendo que las
uniones pueden girar libremente.

Solución.

Primero encontremos los estados posibles para este sistema fı́sico. Como describe el
planteamiento, cada elemento de la cadena puede girar y moverse libremente, y se desea
medir la longitud en el eje horizontal entre sus extremos bajo distintas configuraciones, dado
un sistema coordenado centrado en su extremo izquierdo. Evidentemente, dependiendo de
la posición de cada elemento, la longitud x en el eje horizontal variará entre 0 ≤ x ≤ na.
Si clasificamos las posiciones horizontales de cada elemento en dos grupos, tendremos
los que se encuentran a la derecha con respecto a su eje de giro d y los que se dirigen a la
izquierda de su eje de giro i. Ası́ tendremos nd elementos que apuntan a la derecha y ni
elementos hacia la izquierda, de tal forma que
n = nd + ni (1)
x = (nd − ni ) a (2)
Si se efectúa (a ∗ 1 + 2) y (a ∗ 1 − 2), se obtienen las siguientes relaciones:
na + x na − x
nd = , nd = (3)
2a 2a
1
 Ahora, podemos apreciar que la cantidad de estados posibles que puede tomar el sistema
estará dado por una combinatoria entre el número total de elementos de la cadena (n) y la
cantidad de elementos que están, ya sea a la izquierda o a la derecha (nd , ni ), de la siguiente
forma
n!
W = . (4)
nd !ni !
Por lo tanto, la entropı́a del sistema será

S = k ln W = k [n ln n + nd ln nd + ni ln ni ]
 
na + x na + x na − x na − x
= k ln n ln n − ln − ln
2a 2a 2a 2a
h n x  n x  n x  n x i
= k ln n ln n − 1+ ln 1+ − 1− ln 1−
2 na 2 na 2 na 2 na
h n x  x  n x  x i
= k n ln n + n ln 2 − n ln n − 1+ ln 1 + − 1− ln 1 −
 2 na na 2 na na
1 x  x  1 x  x
= kn ln 2 − 1+ ln 1 + − 1− ln 1 − (5)
2 na na 2 na na

Vemos pues, que la entropı́a de este sistema no depende de la temperatura de la cadena,

sino únicamente de su configuración.
Ahora bien,
Debido a que no se menciona un potencial o fuerza externa que influya sobre el sistema,
podemos proponer el hamiltoniano de la cadena como una combinación lineal de la energı́a
cinética para cada elemento
n
X p2i
H= =E (6)
i=1
2m

Por lo tanto, la energı́a no dependerá de la longitud x de la cadena.

A partir de la ecuación para la energı́a libre de Helmholtz
X X
dF = E − SdT − Xi dxi + µj dNj , (7)

es posible obtener la tensión del sistema, como su variación con respecto a la longitud de la
cadena a temperatura constante (si consideramos que para este caso no existe contribución
del potencial quı́mico)
     
∂F ∂F ∂S ∂S
X = = = −T (8)
∂x T ∂S x ∂x T ∂x
−1   x  x  kT  ln(1 + x/na)  kT x
= −nkT ln 1 + − 1 − ln 1 − +1 = = arctanh (9)
2na na na 2a ln(1 − x/na) a na

Por lo tanto, hemos encontrado una relación entre la tensión de la cadena y su temper-
atura.
Sin embargo, podemos realizar un análisis más detallado si aproximamos la función
tanh−1 (y), y = x/na por serie de Taylor alrededor de y = 0, que corresponde fı́sicamente

2
a cuando x << na. Puede verse que
 
kT −1 d −1 2 d −1
X = tanh (y = 0) + y tanh (y = 0) + y tanh (y = 0) + ...
a dy dy
kT kT
= (y + O(3)) = 2 x + ... (10)
a na
Por lo tanto, cuando la longitud horizontal de la cadena es considerablemente pequeña,
su tensión varı́a de forma lineal con respecto a la temperatura, y cuando se comienza a
x
extender, su tensión aumenta rápidamente en función de la tangente hiperbólica de na .

1.2 Problema 2
Sea ps la probabilidad de que un sistema esté en un estado s con energı́a Es , muestre que si
la entropı́a se define como
N
X
S = −k ps ln ps , (11)
s=1

los valores de ps que hacen a S un máximo bajo la condición de que la energı́a promedio del
sistema es E, sigue la distribución canónica.

Solución.

Tenemos la condición de que la entropı́a promedio debe ser E y que la probabilidad total
de todos los N estados, es decir
N
X
Es ps = Ē = E (12)
s=1
N
X
ps = 1 (13)
s=1

Para resolver este problema con restricciones, empleamos los multiplicadores de Lagrange,
con el fin de encontrar los valores para los que se tiene la entropı́a máxima. Por lo tanto las
ecuaciones para la energı́a y probabilidad serán las restricciones a las cuales asignaremos las
incógnitas λ y κ y las sumamos a la entropı́a (S + κE + λ).
Para maximizar la entropı́a, hacemos una expansión en serie de potencias alrededor de
la variación de cada ps e igualamos a cero. Entonces procedemos con la ecuación ?? como
N
X
δ (S + κE + λ) = δ (−kps ln ps + κEs ps + λps ) = 0
s=1
N N
X 1X k
= (−k ln ps − k + κEs + λ) δps − (δps )2 + ... = 0
s=1
2 s=1 ps
(14)

3
Por lo tanto, para los términos lineales de δps , obtenemos que

−k ln ps − k + κEs + λ
κ λ
ln ps = −1 + Es +
k k
−1+ κ Es + λ
ps = e k k (15)

Por lo tanto, ahora debemos determinar los coeficientes κ y λ para los cuales se satisfagan
las condiciones para la probabilidad total y el valor esperado de la energı́a. Sustituyendo lo
que hemos encontrado para ps , en dichas ecuaciones, tenemos que
N N
−1+ κ E +λ λ κ
X X
−1
E = Es e k s k =e k Es e k Es (16)
s=1 s=1
N N
−1+ κ E +λ λ κ
X X
−1
1 = e k s k =e k e k Es (17)
s=1 s=1

Ahora bien, como Es puede tomar valores tanto pequeños como grandes, para que la
primer ecuación converja, debe pasar que el término κ < 0, por lo que podemos definir
κ
k
= −β, β > 0. Ası́
N
λ
X
−1
E = e k Es e−βEs (18)
s=1
N
λ
X
1 = e k −1 e−βEs (19)
s=1
(20)
PN
La sumatoria en el último término es equivalente a la distribución Z(β) = s=1 e−βEs
para ensambles canónicos. Ası́, podemos ver que
λ
e k −1 e−βEs e−βEs
ps = λ = (21)
e k −1 Z(β) Z(β)

1.3 Problema 3
Cuando una partı́cula de spin 1/2 se sitúa en un campo magnético H, su nivel de energı́a varı́a
entre −µH y +µH, con un momento magnético −µ y +µ en dirección al campo magnético,
respectivamente. Suponga que un sistema de N número de estas partı́culas está en medio de
un campo magnético H y se mantiene a una temperatura T . Encuentre la energı́a interna,
entropı́a, calor especı́fico y el momento magnético total M de este sistema con la ayuda de
la distribución canónica.

Solución.

El espı́n asociado a cada partı́cula es independiente, por lo que la función de partición

para las N partı́culas puede considerarse como el producto de las funciones asociadas a cada

4
espı́n. Es decir

ZN = Z1 · Z2 · · · Zn = ZnN , donde Zi = e−βE = eβµH + e−βµH

= 2cosh (βµH) (22)

Por lo tanto
 
N N N N µH
ZN = 2 cosh (βµH) = 2 cosh (23)
kT

Ası́, para la energı́a tendremos

  N µH


∂ ln ZN ∂ µH sinh kT
Ē = − =− ln 2cosh = −N µH µH

= −N µH tanh (βµH)(24)
∂β ∂β kT cosh kT

Para la energı́a libre de Helmholtz, tenemos

F = −kT ln ZN = −N kT ln (2cosh (βµH)) (25)

Y la entropı́a
 
∂F µH
S=− = −N k ln (2cosh (βµH)) − tanh (βµH) (26)
∂T kT

Por otro lado, para el momento magnético tenemos que

∂F sinh (βµH)
M =− = N kT βµ ln ()
∂H cosh (βµH)
= N µtanh (βµH) (27)

Y finalmente para el calor especı́fico a volumen constante

∂E N µH 2 N k (βµH)2
CV = = sech (βµH) = (28)
∂T kT 2 cosh2 (βµH)

1.4 Problema 4
¯
2
Encuentre la fructuación M − ¯M del momento total M del sistema de espines que se
discutió en el problema 3.

Solución.

Podemos abordar este problema definiendo primero el valor promedio del momento
magnético, el cual será
n n
X X µe−βEn eβµH − e−βµH
M̄ = µPn = = Nµ
i=1 i=1
Z(β) eβµH + e−βµH
1 ∂ ln z
= N µtanh(βµH) = (29)
H ∂β

5
Asimismo
n n
X µ2 e−βEn 1 X 1 ∂ 2 −βEn 1 ∂ 2z
M̄ 2 = = e = (30)
i=1
Z(β) Z i=1 H 2 ∂β 2 H 2 Z ∂β 2
Entonces
2
1 ∂ 2z 1 ∂ 2 ln z

¯
2 ¯ 2 1 ∂ ln z
M − ¯M = M 2 − 2M M̄ + M̄M = M̄ 2 − M̄ 2 = − =
H 2 Z ∂β 2 H ∂β H 2 ∂β 2
= −N µ2 sech2 (βµH) (31)

1.5 Problema 5
Considere la distribución espacial de partı́culas con carga e dadas por el número de densidad
n(r). Si el potencial de un campo externo es Φext , la energı́a potencial total es
n(r)n(r0 )
Z Z Z
1 2 0
U= e drdr + e n(r)Φext dr, (32)
2 |r − r0 |
Asuma que la energı́a del sistema es
Z
S = −k n(r) (ln n(r) − 1) dr (33)

y encuentre las ecuaciones donde n (r) y el potencial estático φ(r) se satisfacen en el

estado de equilibrio.

Solución.

Para este problema debemos encontrar un n(r) tal que se encuentre un mı́nimo para la
energı́a libre de Helmholtz. Por lo tanto, emplearemos el método de los multiplicadores de
Lagrange.
Veamos las variaciones respecto a r. Tenemos que
n(r0 )dr0
Z Z Z
2
δU = e δn(r)dr + e Φext δn(r)dr (34)
|r − r0 |
Además
Z
T δS = − kT ln nδn(r)dr (35)

Y para el número total de partı́culas

Z
δN = δn(r)dr (36)

Ahora igualamos la variación de la energı́a libre de Helmholtz a cero, para obtener su

máximo y agregamos un multiplicador de Lagrange para N
δF = δU − T δS − µδN = 0
n(r0 )dr0
Z Z Z Z Z
2
= e δn(r)dr + e Φext δn(r)dr kT ln nδn(r)dr − µ δn(r)dr
|r − r0 |
n(r0 )dr0
 Z Z Z Z 
2
= e dr + e Φext dr kT ln ndr − µ δn(r) = 0
|r − r0 |
(37)

6
Ası́, para el término lineal δn resulta

n(r0 )dr0
Z Z
2
e + eΦext + kT ln n − µ = 0 (38)
|r − r0 |

El potencial eléctrico en este caso se puede considerar como una suma entre el potencial
externo y el potencial debido a las partı́culas, es decir

n(r0 )
Z
φ(r) = φe xt(r) + dr0 + eext (39)
|r − r0 |

Si sustituimos esto en la ecuación previa, vemos que

eφ(r) + kT ln n − µ = 0 (40)

Y despejando para n
1 eφ µ
n = e− kT (eφ−µ) = n0 e− kT , con n0 = e kT (41)