Mecánica Estadística I - 2025
Departamento de Física - UNLP 06/06/2025
matera@fisica.unlp.edu.ar
rarias@fisica.unlp.edu.ar
Práctica 6: Magnetismo - Gases Reales - Cadenas de espines
Magnetismo
1. Un modelo sencillo de magnetismo en sólidos consiste en considerar a un sólido como un
sistema de N dipolos iguales arreglados en una determinada red. En estos dipolos, el
momento magnético de cada dipolo i viene dado por µ ⃗ i = −e eℏ
m Ŝi = 2 m σ̂i , donde Ŝi es el
1
operador momento angular de espín, representado para espín 2 por las matrices de Pauli σ̂.
El único tipo de energía existente en este sistema es la energía potencial magnética.
Considere este modelo bajo la acción de un campo magnético externo B ⃗ = (0, 0, |B|).
⃗
a) Escribir la función de partición del sistema. Discuta la base elegida para el cálculo.
b) Calcular el valor medio de la magnetización en el eje z en función del campo magnético
aplicado, la cantidad de dipolos y la temperatura. Interpretar físicamente el resultado
−24
en función de la magnitud x = µkB T , siendo µB = 9,21 × 10
B B
J.T −1 el magnetón de
−23 −1
Bohr y kB = 1,38 × 10 J.K la constante de Boltzmann.
c) Calcular la energía interna, la entropía, las probabilidades de que medir un espín up y
down y el calor específico. Graficar estas cantidades en función de la temperatura e
interpretar los límites T → 0 y T → ∞ en cada caso.
d) Mostrar que, para altas temperaturas, se recupera la ley de Curie para la
susceptibilidad magnética: χM ∝ C
T , con C una constante.
2. Considerar ahora el caso general, en donde el momento angular de espín total J puede
tomar cualquier valor entero o semi-entero. En este caso, la función de partición está dada
∑J mJ gµB B
por Z = mJ =−J e kB T donde g es el factor de Landé gJ = 3J(J+1)+S(S+1)−L(L+1)
2J(J+1) , que
para el caso anterior donde consideramos espín 2 resulta g ≈ 2.
1
Gases Reales
3. Utilizar el potencial de Sutherland
{
∞, r<D
u(r) = ( D )6
−ϵ r r>D
para calcular el segundo coeficiente del virial de un gas clásico. También determinar las
correcciones de primer orden a la ley del gas ideal para las distintas cantidades
termodinámicas del sistema.
1
4. Considerar un gas de partículas en d-dimensiones interactuando mediante el potencial
∞, 0 < r < a,
V (r) = −ϵ, a < r < b,
0, b < r < ∞.
a) Calcular el segundo coeficiente del virial y comentar sobre sus límites de temperaturas
altas y bajas.
b) Calcular la primer corrección a la compresibilidad isotérmica
( )
1 ∂V
κT = − .
V ∂P T,N
c) En el límite de alta temperatura llevar la ecuación de estado a la forma de van der
Waals e identificar los parámetros de van der Waals.
d) Para b = a y d = 1 calcular el tercer coeficiente del virial.
Cadenas de espines
5. Una red con N + 1 sitios tiene espines Si = ±1 en cada sitio, todos ellos afectados por un
campo magnético. Hay interacciones de igual intensidad entre uno de los espines, S0 , y el
resto. Luego, el Hamiltoniano es
∑
N ∑
N
H = −h Si − J Si S0 .
i=0 i=1
a) Encontrar la función de partición en el ensamble canónico Z(T, N ), la energía media
⟨E⟩ y, para i ̸= 0, ⟨Si ⟩ y ⟨S0 Si ⟩. Calcular los límites cuando h → 0 y J ̸= 0 y cuando
J → 0 con h ̸= 0.
b) Cuando h = 0 mostrar, para i ̸= j y j, i ̸= 0 que ⟨Si Sj ⟩ = ⟨S0 Si ⟩⟨S0 Sj ⟩.
6. Un modelo de Ising de spin 1/2 con interacción entre 3 espines en una red triangular tiene el
Hamiltoniano ∑
H=− JSi Sj Sk .
i,j,k
La suma es sobre los primeros vecinos y la interacción es ferromagnética, J > 0.
a) Usando a una expansión en las fluctuaciones a primer orden (es decir,
Si = mi + (Si − mi ) y expandir), encontrar el Hamiltoniano de campo medio.
(mi = ⟨Si ⟩ es la magnetización en el sitio i).
b) Usando el Hamiltoniano de campo medio y asumiendo un estado homogeneo en donde
mi = m, encontrar una expresión para la energía libre de Helmholtz y la ecuación de
campo medio para este problema.
