Física Estadística - (Joaquín J. Torres Agudo)
Física Estadística - (Joaquín J. Torres Agudo)
Física Estadística - (Joaquín J. Torres Agudo)
2
Bernoulli determin las correcciones a esta ley cuando uno
tena en cuenta el tamao de las molculas: aparece enton-
ces una dependencia con (V b), donde b es el volumen
que ocupan las molculas.
Joule demostr que P v 2, asumiendo que v = vi (para
todas las moculas).
Utilizando la misma hiptesis sobre v y suponiendo que
para todo tiempo, de acuerdo con las leyes de la probabi-
lidad, 1/6 de las molculas se mueve en cada una de las
6 direcciones independientes del espacio, Kroning (1856)
fue ms all y demostr que
1
P = nmv 2 (1)
3
(n densidad de molculas y m masa molecular).
En este punto aparece el trabajo de Clausius quien, sin
tener en cuenta las hiptesis tan restrictivas de Kroning,
deriva la ley de los gases ideales, y demuestra que (1) es
P
cierta con v 2 = N1 i vi2 , la velocidad cuadrtica media
de las molculas, teniendo en cuenta que las partculas se
mueven en lnea recta con velocidades vi al azar. Tambin
fu el primero en analizar fenmenos de transporte. Con l
se introducen de forma definitiva los conceptos de dinmi-
ca microscpica y estadstica en la teora (hay que hacer
promedios sobre la colectividad de molculas). Su trabajo
influy a Boltzmann; Maxwell lo considera el padre de la
teora cintica y Gibbs lo considera el padre de la ME.
Maxwell fue atrado por el trabajo de Clausius, y en 1860
publica su famosa ley de la distribucin de las velocidades
moleculares, derivada en base a principios elementales de
T. de probabilidades e inspirada en la ley Gaussiana de la
distribucin de errores aleatorios. Deriv una equacin del
transporte suponiendo que la distribucin de equilibrio de
velocidades debera permanecer invariante bajo colisiones
moleculares, y da los mismos resultados que las ecuaciones
ms fundamentales y posteriores de Boltzmann.
El trabajo de Boltzmann:
3
Entre 1868-71, Boltzmann generaliza el trabajo de Max-
well al caso de gases poliatmicos con posibilidad de fuer-
zas externas si stas dan lugar al famoso factor de Boltz-
mann e , con la energa total de la molcula.
el teorema de equiparticin que establece la contribucin
de cada grado de libertad a la energa total del sistema.
Muestra que esta distribucin generalizada (Maxwell-Boltzmann)
es estacionaria con respecto a las colisiones moleculares.
En 1872 deriva el Teorema H que da una base molecular
a la tendencia natural de los sistemas a alcanzar y perma-
necer en un estado de equilibrio, y establece de forma muy
fundamental la conexin entre la descripcin microscpica
propia de la ME y la fenomenolgica de la termodinmica.
Tambin establece un procedimiento para calcular la en-
tropa de un sistema fsico partiendo de un punto de vista
microscpico.
Usando el teorema H, muestra que cualquier distribucin
que venga afectada por colisiones relaja a la distribucin
de Maxwell-Boltzmann en el equilibrio.
En 1876 deriva su famosa ecuacin de transporte, de gran
utilidad para el estudio de los sistemas fuera del equilibrio.
El Teorema H y por tanto el carcter irreversible de los sis-
temas fsicos es criticado por Loschmidt (1876-7) y Zermelo
(1896) (que es un crtica a la teora cintica). Loschmidt se
pregunta, por ejemplo, cmo las consecuencias del teorema
puede conciliarse con el caracter reversible de las ecuaciones
de movimiento de las molculas. Zermelo, usando el teorema
de recurrencia de Poincar para sistemas cerrados con ener-
ga constante, que tienen un comportamiento dinmico cuasi-
peridico, muestra que las consecuencias del teorema H con-
tradicen la segunda ley de la termodinmica (pues a lo largo
de una de estas trayectorias quasiperidicas la entropa, que
es una funcin del microestado, no siempre crece si uno es-
pera un tiempo suficientemente largo). Boltzmann se defiende
aludiendo al papel de las fluctuaciones.
Al mismo tiempo que la crtica a la teora cintica aumenta,
4
se deserrolla la teora de colectividades (TC) que en s es la
ME propiamente dicha: En esta teora el estado del sistema,
que es funcin de sus coordenadas y momentos generalizados,
viene descrito por un punto en el espacio de las fases, (qi , pi)
o microestado, cuya evolucin viene descrita por las eqs. del
movimiento y est restringida a una cierta regin de acuerdo
con las ligaduras fsicas del sistema. Una colectividad es un
conjunto infinitamente grande de sistemas idnticos, tomados
en cualquier tiempo, con las mismas ligaduras fsicas, que di-
fieren slo en sus estados dinmicos ((qi, pi)). Uno tiene en la
regin accesible del espacio de las fases un enjambre de puntos
que se mueven con t.
La cantidad ms importante de la TC es la funcin densidad
(q, p, t) de los microestados. Cuando es estacionaria (/t=0)
caracteriza la colectividad estacionaria que representa a un sis-
tema en equilibrio. Tanto Maxwell como Boltzmann estudia-
ron este tipo de colectividades cuando depende slo de la
energa del sistema, que incluye a los sistemas ergdicos (son
aquellos tales que partiendo de un microestado compatible con
la energa E, si esperamos un tiempo suficientemente largo la
evolucin de dicho microestado pasa tan cerca como queramos
de cualquier otro microestado compatible con la energa E).
Se tiene entonces que para una cierta magnitud que depende
del microestado f (),
Z Z t0 +T
hf i = d ()f () = lm dt f (t) = f. (2)
T t0
5
Gibbs es tan general que se aplica a cualquier sistema fsico
tal que (i) sea mecnico en su estructura y (ii) obedezca las
ecs. del movimiento de Lagrange y Hamilton.
Estos desarrollos de Gibbs coinciden con la revolucin en la
fsica llevada por Planck (1900) con su hiptesis cuntica que
le sirvi para resolver el problema de la radiacin del cuerpo
negro, problema estudiado tanto por la mecnica, la termo-
dinmica y la electrodinmica, que eran las disciplinas mejor
establecidas del siglo XIX.
Tras el trabajo de Einstein (1905) sobre el efecto fotoelctrico y
el de Compton (1923) sobre dispersin de rayos X, se establece
la existencia de un cuanto de radiacin (el fotn). Se puede
entonces tratar el problema de la radiacin del cuerpo negro
mediante el estudio mecnico-estadstico de un gas de fotones
y derivar la frmula de Planck al igual que Maxwell deriv su
ley de distribucin de velocidades para un gas de molculas, y
as entender por qu son tan diferentes.
La respuesta a esta pregunta fue resuelta por Bose. En su fa-
moso trabajo de 1924, Bose deriva la frmula de Planck consi-
derando la radiacin del cuerpo negro como un gas de fotones,
pero en vez de considerar la posicin de cada fotn sobre los
estados de energa del sistema, consider el nmero de esta-
dos que contienen un nmero particular de fotones, que en
la prctica supona el considerar los fotones como partculas
indistinguibles.
Usando la hiptesis de De Broglie, Einstein argument que
lo que es vlido para fotones tambin lo sera para partcu-
las materiales (la indistinguibilidad proviene del carcter de
onda de las partculas) y en 1924-1925 aplica el mtodo de
Bose para el estudio de un gas ideal derivando lo que se co-
noce como estadstica Bose-Einstein, que difiere de la clsica
de Maxwell-Boltzmann por el carcter de indistinguibilidad de
las partculas. Se trata de ver de cuntas formas distintas gi es-
tados de energa pueden acomodar ni = 0, 1, 2... partculas (i)
distinguibles o (ii) indistinguibles. En estos trabajos Einstein
descubre el fenmeno de la condensacin de Bose-Einstein, que
6
fue utilizado por London para la comprensin microscpica de
las propiedades del He lquido a bajas T.
Utilizando el principio de exclusin de Pauli, Fermi (1926)
mostr que ciertos sistemas fsicos obedecen un tipo diferente
de estadstica, la estadstica de Fermi-Dirac, segn la cual no
ms de una partcula puede ocupar el mismo estado de energa
(ni = 0, 1). Esta estadstica fue aplicada por Fowler (1926) y
completada por Chandrasekhar (1939) para discutir los esta-
dos de equilibrio de las estrellas enanas blancas, y por Pauli
para explicar el dbil paramagnetismo, independiente de la
temperatura, observado en los metales alcalinos. Sommerfeld
desarroll la teora electrnica de los metales y Thomas (1927)
y Fermi (1928) estudiaron la distribucin de electrones en to-
mos pesados obteniendo estimaciones tericas para sus niveles
energticos, y desarrollando el llamado modelo de tomo de
Thomas-Fermi, aplicado a molculas, slidos y ncleos.
Toda la estructura de la mecnica estadstica tuvo que ser
reacondicionada para incluir el concepto de indistinguibilidad
de partculas. Al carcter estadstico debido al gran nmero
de partculas se le aade ahora los aspectos estadsticos debido
a la naturaleza ondulatoria de la descripcin mecnica. Esto
implica hacer un doble promedio de las variables dinmicas
sobre los estados del sistema.
Se hizo necesario reformular la teora de colectividades. Pri-
mero Landau y Von Neumann introducen la matriz densidad,
anlogo cuntico de la funcin densidad, tanto en la colec-
tividad cannica como en la microcannica. Segundo, Pauli
introduce la funcin de particin macrocannica (1927), que
permite tratar estadsticamente los nmero de ocupacin de
los estados del sistema.
Como ltimo hito importante Belifante (1939) y Pauli (1940)
descubren la conexin spn-estadstica por la que se explica por
qu algunos sistemas siguen la estadstica de Bose-Einstein
y otros la de Fermi-Dirac. Aquellos formados por partculas
cuyo espn es mltiplo entero de siguen la Bose-Einstein y
7
los formados por partculas con espn semientero impar la de
Fermi-Dirac.
Por ltimo decir que en los ltimos aos ha habido muchos
avances tericos que han servido para perfeccionar las tcnicas
matemticas usadas por la ME incluyendo su extensin a pro-
blemas fuera del equilibrio. Tambin el uso de los ordenadores
ha permitido la extensin del estudio mecnico-estadstico a
nuevos sistemas fsicos y otros sistemas complejos con un gran
nmero de constituyentes elementales.
8
9
Concepto y aplicabilidad de la Fsica Estadstica
Objetivo curso: Entender la Fsica Estadstica y dominar algu-
nas de sus tcnicas, imprescindible para comprender los sistemas
complejos y, en particular, las propiedades de la materia (incluso
plasmas) y la radiacin.
Como hemos dicho antes el curso de Fsica Estadstica es un curso
de introduccin a la Mecnica Estadstica, que es la mecnica de
sistemas fsicos formados por un conjunto grande de elementos sen-
cillos como: molculas en un gas, tomos en un cristal, fotones en
un haz de un lser, estrellas en una galaxia, coches en una carretera,
personas en un grupo social, rboles en un bosque, individuos en
un ecosistema, neuronas en el cerebro, etc.
Antes de seguir conviene recordar qu significa mecnica. La me-
cnica es la ciencia que se ocupa del estudio de los movimientos
conocidos (no slo los del espacio real, que los estudia la Mecnica
clsica):
Mecnica clsica: Ciencia que se ocupa del estudio del movi-
miento de los cuerpos macroscpicos, idealizados como puntos,
con muy pocos grados de libertad.
Mecnica cuntica: Mecnica ms precisa para describir mo-
vimientos del mismo tipo pero a escala muy pequea.
Estas mecnicas no son tiles cuando consideramos objetos con mu-
chos grados de libertad. El efecto de esta gran cantidad de grados
de libertad es las aparicin de otros movimientos o comporta-
mientos colectivos complejos. La ME es la ciencia que estudia estos
otros movimientos y las propiedades relacionadas con ellos como:
cambios de fase, autoorganizacin, formacin de patrones,
memoria asociativa, caos, criticalidad, etc.
Interesan sistemas con muchos grados de libertad, especialmente
grandes colecciones de objetos relativamente sencillos como, por
ejemplo:
muestra de gas, lquido o slido como conjunto de molculas
o de plasma como conjunto de iones y electrones, en cu-
yo caso interesan las consecuencias de los grados de libertad
10
de stas y del conjunto, incluyendo calores especficos, magne-
tizacin, viscosidad, susceptibilidades, y sus fluctuaciones, al
modificar los parmetros del sistema;
No necesariamente han de ser sistemas formados por muchas
partculas (tomos o molculas). Cualquier sistema con mu-
chos grados de libertad es suceptible de estudio por parte de
la FE:
Una molcula compleja (polmeros, DNA,...) como con-
junto de tomos;
Los tomos y ncleos muy pesados, como conjunto de
e, p, n y otras partculas elementales;
Un rayo lser y semejantes como conjunto de fotones;
Una galaxia como conjunto de estrellas;
Una comunidad como conjunto de indivduos, etc.
Veamos como ejemplo el caso de un fluido. En este caso la Fsica
Estadstica se pregunta:
Cada molcula del fluido tiene generalmente una di-
nmica (microscpica) casi trivial: interacciona con las
restantes segn las ecuaciones del movimiento de Newton.
Pero, a nivel macroscpico, se observa un comportamien-
to hidrodinmico que no comprendemos bien.
Esta paradoja es dramtica cuando, por alguna causa, esos movi-
mientos microscpicos se organizan para producir tornados o turbu-
lencia. Parece lgico pensar que la extraodinaria complejidad de un
flujo turbulento sea una manifestacin ms de la dinmica colectiva
emergente.
Cmo es esto posible? Cmo se produce esa complejidad?
La naturaleza tiene estructura jerrquica en capas (posibles
niveles de descripcin)!
Cmo ha de ser de grande un sistema es decir, qu nmero
de grados de libertad ha de contener para que, en un nivel
alto de descripcin pueda considerarse un medio continuo que
satisface las leyes fsicas (macroscpicas)? (en algunos casos
ha de ser del orden de NA )
11
Por qu se parecen tanto entre s fenmenos (macroscpicos)
aparentemente tan dispares como ferromagnetismo, condensa-
cin, superconductividad, turbulencia?
Cul es el origen ltimo de los principios de la termodinmica
o del comportamiento irreversible de muchos cambios natura-
les con el tiempo?
12
teora de procesos estocsticos: cuando la dinmica mi-
croscpica es desconocida o muy complicada (no conocemos el
H(q, p)), se hace la hiptesis de que los sistemas grandes rea-
lizan sus movimientos elementales al azar (es una descripcin
menos fina de los sistemas fsicos).
=
Simetra esfrica Simetra orthorrmbica
Sublimacin del yodo
irreversibilidad macroscpica: desaparece la inversin tem-
poral microscpica (las ecuaciones de Hamilton son reversibles
temporalmente).
13
Cuando tenemos un nmero pequeo de partculas se tiene:
~v ~v
=
cuando t t
14
Nota: Antes de profundizar en la naturaleza y mtodo de la FE,
repasaremos ciertas propiedades de las descricpciones microscpica
y macroscpica; lo hacemos con mucha brevedad, puesto que las
han de conocer bien, pero hay que insistir en su importancia para
lo que sigue. Los alumnos han de saber mecnica!
15
Tema 2: Construccin de la Fsica Estadstica:
Naturaleza y Mtodo
Descripcin microscpica clsica. Descripcin observable
(Balescu)
El estudio de la dinmica de los sistemas fsicos puede hacerse en
el marco de dos descripciones que proporcionan las Mecnica Cl-
sica y la Cuntica, respectivamente. Empecemos por estudiar la
descripcin microscpica clsica:
qs , p s ; s = 1, 2, ..., ,
3
Balescu, pginas 3 a 8, por ejemplo, para el rpido repaso que sigue de dinmica hamil-
toniana. El alumno ha de completarlo con sus libros de mecnica clsica y cuntica.
16
esto es:
17
no se corta a s misma, pues las ecs. de H. son diferenciales
de 1er orden.
Por la misma razn (como reflejo del determinismo cl-
sico):
0
t t
X
b H b H b b
= + |{b,{zH}} + ,
s=1
qs ps ps qs t t
Poisson
18
Casos particulares de (3):
Consecuencias de (3):
b = {b, H} Lb,
b (t) = U (t) b
19
Ms explcitamente, es
b (t ) = etL b (0 ) , b (0 ) = etL b (t )
nm q1n1 ...pm
(simplificando notacin).
(10)
n,m
con los mismos coeficientes nm . Pero, usando aqu (7) y (8) para
sustituir q1 (t) , ..., p (t) por sus valores iniciales y el t, se tiene:
X
b (t) = nm (t) q1n1 ...pm
,
n,m
20
esto es, una nueva funcin (pues tiene distintos coeficientes) de las
variables originales.
La misma funcin se puede definir en funcin de la base en t = 0
con coeficientes
nm (t), o en la base en t con coeficientes constantes
nm
21
+ nmeros molares, si hay varios componentes qumicos
+ momento de dipolo elctrico, si no es elctricamente neutro,
etc.
En general (ej, hidrodinmica) el macroestado es ms complejo
que el microestado: son necesarios campos B (~r, t) , es decir,
funciones contnuas (o contnuas a trozos) con dependencia
espacial y temporal gobernados por ecuaciones en derivadas
parciales o integrodiferenciales, con gradientes, que no tienen
invariancia temporal, es decir, son irreversibles; B (~r, t) =
B = const. es el caso particular de la Termodinmica del
equilibrio.
El macroestado no determina por completo al :
dos sistemas en el mismo macroestado pueden encontrarse en
distintos : existe degeneracin
llamaremos imagen dinmica de A:
(A) = conjunto de s compatibles con macroestado A
Sin embargo, supondremos que, en general, que lo contrario
s es verdad, el determina por completo el A, en el sentido
de que dos sistemas en el mismo han de encontrarse en el
mismo A.
A, como , puede cambiar con t, pero su evolucin NO es
determinista:
a) dos sistemas en mismo A en t0 (trivial: mujeres gestantes) pue-
den estar en distintos y conducir a distintos A en t > t0
(mujer+nio/a)
b) ejemplo sistema no conservativo:
en t0 , sistema 1 en macroestado energa E y microes-
tado 0
en t0 , sistema 2 en macroestado E pero microestado
0
t0 t : se sale de superficie de energa E,
0 t H (t ) = E
H (t ) = E 6= E
0 t
22
luego NO puede predecirse con confianza el A futuro,
Propugnamos:
reemplazar nuestro desconocimiento de condiciones ini-
ciales imposible y dudosamente conveniente
por estrategia consistente en mezclar teora de probabi-
lidades con el conocimiento de la forma genrica de las
ecs. de Hamilton para el problema
La FE ensea cmo hacerlo
23
Descripcin microscpica cuntica.
1111
0000
0000
1111
0000
1111 q
0000 s
1111
ps
24
esto es, la funcin de onda |i , ej, {q1 , ...q } u
otro conjunto que caracteriza microestado (no pode-
mos elegir ninguna pareja qs , ps dentro de este conjun-
to pues lo impide el PI), juega el papel del conjunto
(q, p) que tenamos en la descripcin clsica.
el espacio de Hilbert es ahora el espacio dinmico.
25
define su dimensin y permite expandir cualquier elemento
de H X
|i = an |ni (11)
n
{qs , ps} = 1
se concluye la analoga
1
{, } [, ]
i~
que sugiere que la ley de movimiento para operadores es
h i
bb = 1 bb, H
b (12)
i~
con Hb el operador hamiltoniano, que tambin juega aqu
papel privilegiado. De hecho esta ley es consecuencia de
la ec. de Schrdinger (Vamos a verlo):
b |i b
i~ |i = H |t i = eitH/~ |0 i ,
t
de donde, dada la hermiticidad de H b:
D E
b b
hbit ht | b | t i = e
b
itH/~ b
0 | be
b
itH/~
0 (13)
=h0 | eitH/~bbeitH/~ | 0 i =h0 | bbt | 0 i hbbt i0 (14)
b b
26
establece la forma de evolucionar los operadores. Derivan-
do, suponiendo bb sin dependencia explcita en t :
b b i +
bb t dbt = i He
b +be + e+ e e be H b (15)
dt ~ t
|{z} ~
=0
i b i b h
i b i 1 h bi
= Hbt bt H = H, bt = bt , H (16)
~ ~ ~ i~
Q.E.D.
27
28
Naturaleza y mtodo de la FE: Colectividad de Gibbs
Resumiendo: Como hemos dicho ya antes, la FE trata de relacionar
las descripciones microscpica (clsica o cuntica) y la descripcin
macroscpica (u observable), combinando la descripcin microsc-
pica con teora de probabilidades para llegar a la descripcin ma-
croscpica.
Problemas:
aplicacin restringida a gases muy diluidos o similares
para evitar correlaciones o interacciones entre los ele-
mentos; pero las interacciones no son despreciables, en
general;
los constituyentes no son idnticos a veces: ej, mezclas
Este es el mtodo de Maxwell y Boltzmann; es la base de
la pionera teora cintica de los gases que puede utilizar-
se para el estudio, por ejemplo, de un gas ideal mono-
componente muy diluido o un conjunto de osciladores
29
armnicos independientes; sin embargo el nnero de apli-
caciones fsicas es muy limitado y se utiliza rramente en
FE.
Para evitar estas limitaciones, Gibbs y Einstein fundaron la
mecnica o fsica estadstica propiamente dicha; proponen
2. Aplicar teora de probabilidades a un conjunto infinito de sis-
temas imaginarios, que se conoce como
colectividad de Gibbs:
es un conjunto infinitamente grande de copias macros-
cpicas imaginarias del sistema dado, todas en el
mismo macroestado, en particular, por ejemplo, con el
mismo hamiltoniano
todas ellas obedeciendo las mismas leyes dinmicas mi-
croscpicas.
cada una de ellas en distinto microestado compatible
con el macroestado dado (es decir lo que hemos lla-
mada imagen dinmica (A)); ej, correspondientes a
gran nmero de condiciones iniciales para las ecs. de
Hamilton (la colectividad de Gibbs va a ser una nu-
be de puntos en en t determinado, compatible con
el macroestado en cuestin; si el macroestado es una
determinada superficie de energa entonces la nube es-
tara sobre dicha superficie.)
A
1111
0000
A A A
0000
1111
0000
1111
0000
1111 1 2 3
0000
1111
4 5 6 7
A A A A i (A)
30
Ventajas:
Miembros cumplen sin dificultad condiciones para apli-
cacin ms eficaz de teora de probabilidades:
Los miembros de la colectividad son idnticos en
sentido macroscpico
son independientes, sin interacciones, por defini-
cin
su nmero es infinito por lo tanto vamos a tener
precisin infinita de los mtodos estadsticos. La
precisin estadstica de los resultados no depende
de nmero de constituyentes del sistema, sino del
nmero de miembros de la colectividad, que toma-
remos .
los miembros son macroscpicos, luego siempre distin-
guibles
= al contrario que en teora cintica, la FE es in-
sensible por completo al tipo de mecnica (clsica o
cuntica) que requieren los constituyentes pues el for-
malismo es idntico.
Consideraciones estadsticas son ms fundamentales
= menor necesidad de detalle microscpico (gene-
ralmente inaccesible); ej, comparar clculo P en teora
cintica siempre hay que hacer una gran descripcin
microscpica. En FE los detalles microscpicos van a
ser ms irrelevantes.
Colectividad es imaginaria no necesariamente realizable
luego hay gran libertad a la hora de elegirla segn las
necesidades matemticas que tengamos (puedo inven-
tarme una colectividad arbitraria para tratar un pro-
blema fsico en cuestin);
en principio, es un truco matemtico aunque (luego
profundizaremos en esta interpretacin fsica) trata de
reproducir gran nmero de experimentos realizados en
condiciones idnticas. Hacer en el papel un estudio es-
tadstico de la colectividad de Gibbs-Einstein va a ser
lo mismo que hacer medidas en el laboratorio de mag-
nitudes fsicas:
31
dada una situacin experimental (ej, gas en vasi-
ja, y nos interesa medir E), el experimentador no
conoce el microestado
si hacemos varias medidas (instantneas) de E,
encontraremos al sistema en distintos microestados
pues (A) es enorme, luego medimos Ei
E (i )
pero al observador le interesa el valor medio de Ei
y, quizs, p (Ei) propiedades estadsticas de las
fluctuaciones de Ei , correlaciones estadsticas en-
tre fluctuaciones de distintas magnitudes, etc.
FE supone que esta informacin puede obtenerse
sin realizar el experimento aplicando teora de
probabilidades a la colectividad (de Gibbs) cuyos
miembros distribuidos adecuadamente entre los mi-
croestados accesibles.
32
Postulados de la FE.
La FE tiene estructura de ciencia, pudiendo desarrollarse a partir
(por ejemplo) del siguiente conjunto de postulados:
Postulado Previo:
Postulado Primero:
I.A:
con ~r1 = {q1,x (t) , q1,y (t) , q1,z (t)} , etc., las posiciones de las
N molculas.
33
Notad que las energas y momentos son magnitudes fsicas mec-
nicas, mientras que lo afirmado en I.A falla, ej, para la T a
la que se le asocia una velocidad cuadrtica media, caracte-
rstica de todo el conjunto y para la S, pues S 6= S () sino
que S = S ( (A)).
I.B:
34
Falta por determinar cmo se reparten los miembros entre
los microestados accesibles; esto puede precisarse introdu-
ciendo, ej, la funcin densidad (de la mecnica clsica),
() definida en , que determina el peso de en la co-
lectividad, sujeta a la condicin:
Z
d () = 1.
y, en el cuntico:
hbi = Tr b
con el operador densidad, anlogo cuntico de la fun-
cin densidad, Tr (
) = 1.
35
Tolman (1938): argumentos convincentes, sin prueba rigurosa to-
dava de que poda ser un buen modelo de la realidad:
Prob [
/ M] = 0
Prob [ m M] medida (ej,vol.) de m
36
Dicho de otro modo, proponemos como segundo postulado:
Postulado II:
E energa E + E, E << E
37
As, al estudiar sistema cuntico en ME involucramos dos prome-
dios:
X D E
i | bi
lm
hbi =
i=1
| {z }
(1)
| {z }
(2)
38
Resumiendo:8
el mtodo mecnico-estadstico de la FE sugiere calcular (nos
referimos ahora, por sencillez, al caso clsico)
Z Z
hbi = d () b () , d () = 1,
M M2
M1
39
A (17) se le conoce como hiptesis ergdica; el estudio de si los
sistemas de inters tienen esta propiedad ha dado lugar a la teora
ergdica . Es decir, el sistema es ergdico cuando la trayectoria
sobre M llena densamente dicha superficie si esperamos un tiempo
suficientemente grande.
Si la regin accesible est formada por dos o ms subregiones disjun-
tas entonces el sistema no es ergdico. Partiendo de una condicin
inicial en una de ellas no se puede llegar a las otras: por ejemplo el
gas ideal, cada partcula de gas tiene una energa que es cte del mo-
vimiento entonces cara partcula me define una subregin accesible
al sistema. Un sistema es ergdico cuando = (H) slo depende
de las consts. del movimiento.
40
Funcin densidad, Teorema de Liouville clsico
Funcin densidad:
41
() valor finito bien definido en todo , lo que excluye
casos difciles, patolgicos, como () ( 0 ) (a
pesar de que son potencialmente interesantes)
42
As, la nube se comporta como un fluido, para el q la de-
rivada total (llamada derivada material o hidrodinmica
en dinmica de fluidos, sera la variacin que vera un ob-
servador que se mueve con el fluido) tiene contribucin
debida al movimiento de sus partculas y otra debida a los
cambios en el entorno de la partcula observadora (sera
la variacin que ve un observador fijo respecto al fluido)
Mt
t
0
M0
(t ) / (0 ) = 1.
43
Pero, si es cte. el no de puntos en M, d, y su medida, dV,
tambin densidad = d/dV = cte., luego
d
= 0,
dt
teorema de Liouville o, equivalentemente,
= {, H} iL,
t
ecuacin de Liouville.
t (Ut ) = 0 () , t,
t () = 0 () () , t,
t (Ut) = (Ut ) = 0 () = ()
11
Notad que esta propiedad es caracterstica de , para la que, en general:
d
=0 y 6= 0
dt t
mientras que se tiene, en general, cualquier otra funcin dinmica:
db b
6= 0 y =0
dt t
44
es decir,
(Ut ) = () , t.
Puesto que y Ut son dos puntos de la misma trayectoria, la
de una colectividad estacionaria es cte. a lo largo de cualquier
trayectoria que pueda recorrer el siste en su movimiento natural.
{, H} = 0,
45
tambin
H ()
() exp , (18)
con y parmetros independientes de para que sea =
f [H ()] .
= f [H () , L () , ...] ,
46
Matriz densidad, Teorema cuntico de Liouville:
1 X (i)
hbi = lm b . (20)
i=1
47
La ltima parte de (20) puede entonces escribirse:
" #
1 X
X
n
k 1 X
b(i) = b(ik ) ,
i=1 k i
k=1 k
Para mejorar analoga con caso clsico, podemos escribir este re-
sultado en trminos de la funcin de onda definida sobre el espacio
de configuracin . Sea k (q) la funcin de onda del estado |k i
y b k (q) la de b |k i , de modo que, si b se expresa en espacio de
configuracin mediante la matriz b (q; q ) hq| b |q i , se sigue:
Z
b k (q) = dq b (q; q ) k (q )
48
donde se ha definido la matriz densidad (en espacio de configu-
racin): X
(q ; q) k (q) fk k (q ) . (26)
k
Notemos la analoga de (25) con el caso clsico, es decir, hbi apa-
rece como el producto de dos funciones de 2 variables: la b (q; q )
determina la variable dinmica en cuestin, mientras que (q ; q)
determina la colect. Por supuesto, (26) no es la nica representacin
posible para la matriz densidad14
Puesto que fk son nmeros reales, la matriz densidad es hermtica,
(q ; q) = (q; q )
49
|2 i ,... Se tiene que16
X
|j i hj | = 1.
j
Se sigue sucesivamente:
X
hbi = fk hk | b |k i (27)
k
XX
= fk hk | j i hj | b|k i (28)
k j
X
= |j i traza b
hj | b (29)
j
con la suma extendida a todo el conjunto. Puesto que esta relacin es vlida | i , se sigue
lo indicado si el conjunto es completo.
50
La nonegatividad se corresponde con fk 0 k, y puede expre-
sarse:
Z Z
dq dq (q ) (q ; q) (q) 0 (q) .
d + i + i +
ut = + u+
t H = u
H,
dt ~ ~ t
17
|i no contiene otras dependencias funcionales, mientras que (q) hq | i depende
del espacio en cuestin.
51
+ = H.
pues H Entonces, derivando, se tiene sucesiva-
mente que
d + d
d
utu+
u ut = + t ut
u+ (31)
dt t dt t
dt
i + +
= ut H ut ut H ut = 0, (32)
~
luego u+ t = cte, de donde u+
t u t = 1 si tomamos u0 = 1.
t u
luego, en tiempo t,
X
t = ut |k i fk hk | u+ t 0 u+
t =u t
k
52
Derivando la ecuacin integrada:
t = u +
t0 u t 0 u
t +u
+
t
de donde, sustituyendo ut = ~i H
ut y u
+ i +
t H :
t = ~u
d
t i u u+ i u u+
= H t 0 t + t 0 t H (35)
dt ~ | {z } ~ | {z }
t t
i
i h i
= H t t H = H, t (36)
~ ~
que es la ec. de Von Neuman y es el anlogo cuntico a la ec. de
Liouville clsica.
Colectividades invariantes
53
Aplicamos la ec. de Von Neumann:
d i h i
0 = + H, ... hm | t | ni (37)
dt ~
d i h i En Em
= + H, ... exp i t hm | 0 | ni(38)
dt ~ ~
h i
Pero, H, 0 = 0 si la colectividad es estacionaria, y se sigue:
i En Em
= (En Em ) exp i t hm | 0 | ni = 0,
~ ~
donde 0 es para un tiempo inicial arbitrario. Es decir, para colec-
tividad estacionaria:
hm | | ni = 0 si Em 6= En :
ES LA COLECTIVIDAD MICROCANNICA!
54
queda determinada por normalizacin:
X am X X
traza
= = 1 = = am = 1
m
m E<Em <E+E
ES LA MACROCANNICA!
55
Tema 3: Colectividades
Colectividad microcannica. Entropa. Limitaciones.
Contacto con bao trmico. Potencial de van Hove. Colectivi-
dad cannica, . Funcin de particin. Degeneracin. Lmites
continuo y clsico. Operador densidad.
Colectividad macrocannica. Funcin de particin. Lmite cl-
sico.
Teoremas de existencia. Lmite termodinmico. Condiciones
sobre el hamiltoniano.
Densidad de energa libre. Estabilidad.
Generalizaciones: sistemas cunticos, plasmas, interacciones
gravitatorias.
Fluctuaciones. Frmula de Einstein. Divergencias. Opalescen-
cia crtica.
Equivalencia macroscpica de colectividades.
57
A lo largo de las lecciones de este bloque desarrollamos otras
colectividades a partir de la microcannica y establecemos
primera relacin con fsica macroscpica.
Tambin vamos a ver algunas propiedades fundamentales de
las colectividades.
H()= E + E
= (E, E; N, V ) = (E, E)
Viene a ser no de formas distintas de distribuir energa E entre N partculas, luego
19
P
depende de E y N (y de E). Tambin de V pues Em (o la suma m 1) es funcin de V.
59
Entropa
R
Determinada , se siguen promedios hbieq = d () b () . Pero
un mtodo ms directo de obtener informacin macroscpica a
partir de la microscpica contenida en o, equivalentemente, en :
60
Limitaciones:
X X
= am = 1 ,
m m
E Em E + E
61
fsicas:
62
Sistemas en equilibrio trmico: colectiviad cannica
Sea sistema aislado U (universo), en el que distinguimos parte, A,
que interacciona con el resto, A:
U = A A A A =
U
R0
A
RA A
63
Supondremos H AA suficiente para tener intercambio eficaz
de energa que permita equilibrio mutuo AA, pero relativa-
mente despreciable frente a EU , para tener EU EA + EA
Supondremos tambin
NA NA = EA EA (tienen la misma naturaleza)
XNX
H = (rij ) , que puede tomarse como definicin
i<j=1
de y representa las interacciones efectivas entre partculas
realista es el de Lenard-Jones, pero por sencillez supone-
mos un potencial de van Hove, que tiene las 2 propiedades
importantes de esta interaccin efectiva realista: ncleo duro
impenetrabilidad y atraccin de corto alcance:
LenardJones van Hove
r d0 R0 r
64
En estas condiciones, si
no
molculas mximo no de molculas dentro
HV en V alcance f atractiva de una dada |0 |
| {z } | 4 3
{z }
R 3
N=N V V = nV
3
4
0
d3
= Rd30
3 0 0
esto es,
R03
HV
0 nV 3 V
d0
Adems, H0 N nV, luego
E V , E V ,
A A A A
que justifica EA EA . Por otra parte, para estimar H AA no-
tamos que, para un potencial de alcance limitado, slo contribuyen
partculas en corredor de anchura R0 (ver figura) luego
H V 2
RA R0 R0 R0
AA
corredor = 0
E V R 3
R 1/3
A A A A V
A VA
65
Sabiendo que U puede representarse por la microcanni-
ca, y suponiendo que el sistema A tiene estados propios
caracterizados por r y valores propios Er , cul es la pro-
babilidad, pr , de encontrar a A en microestado con valor
propio de la energa Er ?20
Notamos:
U EU , EU + E (E EU ) :
probabilidad de que esto ocurra = U (EU ; E)
A Er : probabilidad = pr
A EA EU Er , EA + E
EA EU
EA EA = = E EA
(E EU )
probabilidad de que esto ocurra = A (EU Er ; E) que es
tambin el nmero de microestados de U en los que A est en
el estado con valor propio Er .
En consecuencia:
casos favorables A (EU Er ; E)
pr = =
casos posibles U (EU ; E)
66
ln A (E)
ln A (EU Er ) = ln A (EU ) + (Er ) + ...
E E=EU
ln A (EU ) Er
donde
ln A (E)
(39)
E E=EU
En consecuencia,
y se sigue
A (EU Er ; E) A (EU ; E) Er 1
pr = e eEr ,
U (EU ; E) U (EU ; E) Z
(40)
donde, segn las definiciones, los parmetros Z y son indepen-
dientes de Er , que se conoce como distribucin cannica (para la
energa) o factor de Boltzmann.
67
con la suma extendida a todos los microestados, sin restriccin
P
alguna, al contrario que en (E, E) = r 1, donde se suma
sobre microestados con energa entre E y E + E.
68
a un caso que tiene un rango de aplicaciones mucho ms amplio
que el de la colectividad microcannica introducida mediante un
postulado.
Los valores medios de operadores vienen dados en esta nueva colec-
tividad: D E
b = Tr b 1
= Tr b e H , (45)
Z
de acuerdo con el mtodo general.
Notemos que, para interpretar adecuadamente las frmulas anterio-
res, por ejemplo, las que se refieren a la representacin de energa,
hay que notar que los sistemas de inters tienen, generalmente, ni-
veles degenerados, es decir, r microestados con el mismo valor Er
de energa. Esta degeneracin ha de tenerse en cuenta en la su-
ma en (42) sobre todos los estados. Es decir, aquella frmula puede
escribirse ms explcitamente:
X
Z= r eEr ,
r
donde la suma es ahora sobre valores propios distintos Er , y la
probabilidad de que el sist est en cualquiera de estos estados es
Pr = Z 1r eEr , que juega el papel de peso del nivel degenerado
Er suponiendo que cada uno de los r tiene la misma probabilidad.
69
Esta expresin sugiere considerar formalmente a como una va-
riable compleja, con lo que la funcin de particin cannica se nos
muestra como la transformada de Laplace de la densidad de esta-
dos (E) . La importancia de este hecho es notable. Por una parte,
(46) establece una relacin formal entre las colectividades microca-
nnica y cannica, lo que resulta esencial en el contexto del llamado
problema del lmite termodinmico.24 Por otra parte, (46) permi-
te calcular a partir de Z que, generalmente, es un problema ms
sencillo.
24
De hecho, tendremos finalmente que
A = 1 ln Z,
S = k ln (c) ,
70
Lmite clsico
() = Z 1 exp [H ()]
71
notar que ste slo puede definirse con incertidumbres acordes con
Heisenberg, es decir, qi y pi tales que
Y
pi qi h .
j=1
Interpretacin de .
Tenemos
" #
e Er
1 X
pr = P = ln eEr (47)
r eEr Er r
72
puede estar justificado) y A1 +A2 +A pueda considerarse como un
sistema aislado describible en la microcannica:
A2
A1
n o n o
(1) (2)
Sean Er y Es los niveles accesibles a A1 y A2 , respectiva-
mente. La probabilidad conjunta de encontrar a A1 en r y a A2 en
s es
1 1 Er(1) 1 2 Es(2)
p(r,s) = e e
Z1 Z2
con X (i)
Zi = eiEr , i = 1, 2.
r
Pero tambin podemos considerar A = A1 + A2 con energas En =
(1) (2)
Er + Es (despreciando las energas de interaccin), en cuyo caso:
1 (Er(1) +Es(2) ) X (1) (2)
pn = e , Z= e (Er +Es ).
Z r,s
73
Potencial cannico.
Definimos
" #
X
A 1 ln eEr = 1 ln Z. (48)
r
(A) X (A)
d (A) = d + dEr usando (48):
r
E r
" # " #
X X X
= ln eEr d + ln eEr dEr
r r
Er r
X
d (A) = d (U ) dU + pr dEr
r
de donde:
X
dU
|{z}
pr dEr .
d [ (U A)] = (49)
| r {z }
aumento
medio de trabajo medio
sobre el exterior
energa
| {z }
sobre el sistema
| {z }
1er Principio: calor transferido desde el exterior dQ
74
Para interpretar este resultado, consideramos un proceso con varia-
Corolarios:
75
energa interna: U = N u = A + T S = kT 2 T ln Z =
ln Z
entalpa: H = A + T S + P V
y magnitudes intensivas:
A 2 a ln Z
presin: P = =n = kT
V N,T n N,T V
u 1 1 2 2
cV = = ln Z = k ln Z
T N,V N T N 2
1
1 V P
compresibilidad: T = = n
V P T n T
!
P 1 X
pr = eEr Er eEr
re = ln pr = Er ln .
| r
{z }
Z
Por otra parte,
X X
S
= (U A) = pr Er + ln Z pr
k
r | r{z }
X X =1
= pr (Er + ln Z) = pr ln pr .
r r
76
Es la entropa cannica o entropa de Gibbs:
X
S = k pr ln pr
r
77
78
Sistemas en equilibrio material: Colectividad Macrocan-
nica
La cannica (N, V, T fijos) es preferible a la microcannica (N, V, E)
por:
79
puede intercambiar E y N, estando el conjunto aislado:
A
V=cte
80
A est en un estado con (Es , Nr ) es
no config de U en las q
A est en (Es , Nr )
z }| {
(EU Es , NU Nr ; E)
pr,s = A o
n total microestados
accesibles a U
(E , N ; E)
| U U {z U }
ln A (EU Es , NU Nr ) ln A (EU , NU )
ln ln
Es A Nr A
E E = EU N E = EU
N = NU N = NU
ln
Sabemos que = (kT ) 1
= A , y definimos tal
E E = EU
N = NU
que
ln A
= ,
N E = EU
N = NU
P
con lo que, exigiendo normalizacin r,s pr,s = 1, se tiene
exp (Es + Nr )
pr,s = X
exp (Es + Nr )
r,s
81
Alternativamente, puede argumentarse (menos rigurosamente): Si
hacemos variaciones infinitesimales a V constante en = (N, V, E)
y usamos las definiciones arriba de T y :
ln ln 1
d (ln ) = dE + dN = dE + dN,
E N,V N V,E kT kT
pero S = k ln ,29 que implica k1dS = d (ln ) , luego
dS = T 1dE T 1dN dE = T dS + dN
que, comparada con la relacin termodinmica = es el potencial
qumico.
yN
ee, suma sobre los valores propios de H en la forma usual.
82
con Z (T, V, 0) 1, que establece relacin fundamental entre dos
funciones de particin: no es sino un desarrollo polinmico infinito
de potencias de la fugacidad y con coeficientes que son las funciones
de particin cannicas para sistemas con N = 0, 1, 2... partculas.
Lmite clsico
Potencial macrocannico.
J (V, z, T ) = V P (T, ) .
83
Este potencial satisface:
J J J
= S, = P, = hN i ,
T V, V T, T,V
de modo que
J 1
= z ln = ln
z
y se sigue la relacin fundamental30
kT
P (T, ) = ln (V, T, ) .
V
30
Puede verse que la constante aditiva de integracin no aparece.
84
Tema 4: El lmite termodinmico
Teoremas de existencia. Lmite termodinmico.
85
Entropa y energa libre para sistemas finitos con interaccio-
nes no presentan la aditividad que requiere la Termodinmica,
salvo en LT.
Teora ergdica: un modelo (mezcla) slo puede presentar com-
portamiento irreversible como el que caracteriza mundo real en
el LT.
Las distintas colectividades slo conducen a la misma descrip-
cin macroscpica en el LT.
86
densidad de energa libre:
)
kT
a
(n, T ) lm A(N,V,T = lm ln Z (N, V, T )
V V
|{z}
V
|{z}
V
Z (V, N ) = Q (V, N ) 3N ,
87
Para este (r) con alcance limitado, dada una partcula, un
no mximo de partculas que pueden interaccionar con ella,
d0
R0
88
de donde:
N ln N N
z }| {
ln Z 3N ln ln (N !) + N ln V + N
kT kT
ln Z (V, N ) N (+3 ln + ln N 1 ln V )
V V
1
A (V, nV ) nkT ln 3 n 1 . (54)
V
En definitiva, estabilidad = cota inferior para A/V, independiente
de V y N y, por tanto, finita en el LT.
89
Hove). Se tiene V1 + V2 < V, debido a la presencia del corredor. En-
tonces, Q (V, N ) Q (V1 + V2 , N ) ya que el integrando es siempre
definido positivo y se reduce el volumen de integracin; como 3N
no depende de V, sino de N que es el mismo, se tiene tambin
X Z Z Z
1 N!
Z (V, N ) d~q1 d~qN1 d~qN1 +1
3N N ! N1 ! (N N1)! V1 V1 V2
N1
Z
d~qN1 +N2 exp [H (~q1 , . . . , ~qN1 ) H (~qN1 +1, . . . , ~qN1+N2 ) N1N2 ]
V2
R R
V2 d~
qN1 +1 V2 d~qN1 +N2 exp [H (~qN1 +1, . . . , ~qN1+N2 )]
que conduce a
90
A (V, nV ) V1 A (V1 , n1V1 ) V2 A (V2 , n2V2 )
+ ,
V V V1 V V2
que establece una cota superior. Por otra parte,
N = N1 + N2 = nV = n1 V1 + n2V2 .
Nk , Vk , Nk /Vk n para k .
91
la secuencia estndar de cubos:
k+1
k
S k+1 + R
R
Sk + R
Sk
Sk + R, k = 1, 2, . . .
Esto es,
h i h i
3 3
A (Sk+1 + R) , Nk+1 A (Sk + R) , Nk
(Sk+1 + R)3 (Sk + R)3
92
(i)
puesto que Sk+1 + R = 2 (Sk + R) y todos los Nk son iguales por
construccin, de modo que, al pasar de un sistema al siguiente en
la secuencia, aumentando el volumen, la energa libre por unidad
de volumen no puede aumentar, sino, si acaso, disminuir.
QED33
93
Es interesante notar que el tratamiento de sistemas cunticos
en este contexto sigue lneas esencialmente paralelas a las del
caso clsico. Las diferencias son tcnicas:
hay que empezar demostrando que, partiendo de f.d.o. con
la simetra apropiada para bosones y fermiones, las condi-
ciones de contorno no afectan a las propiedades volum-
tricas en el LT. Aunque es obvio que ha de esperarse esta
propiedad, no es sencillo demostrarla en general (Ruelle
1969);
luego hay que probar que la operacin traza, que aparece
en la misma definicin de Z,
Z (V, N ) = Tr eH (57)
se mantiene bien definida (Ruelle 1969, p.60);
la condicin de estabilidad en este caso consiste en que
el valor propio ms bajo del H ha de tener cota inferior,
mn (n ) N . Se sigue entonces que la condicin cl-
sica de estabilidad (H N ) implica la cuntica.36
Es inmediato probar que V1 + V2 < V implica
Z (V, N ) Z (V1 , N1) Z (V2 , N2) .
En efecto, puesto que la energa de interaccin atractiva
es negativa, cuando desaparece la posibilidad de atraccin
entre dos partculas (una en 1 y otra en 2), es como si
apareciese un trmino positivo en el H, lo que aumenta la
magnitud de sus autovalores y, en consecuencia, disminuye
la traza (57).
Los restantes argumentos son geomtricos y, por tanto, no
se modifican respecto del caso clsico.
94
En qu medida los sistemas reales son estables y moderados?
95
Pero puede encontrarse una cota inferior para la energa
total como consecuencia directa de los principios cunti-
cos, y se sigue una cota inferior para la densidad de energa
libre igual que en el caso estndar.40
40
Por ejemplo, Dyson y Lenard (1967,1968) han mostrado que, para N fermiones en
presencia de cargas positivas, la energa del estado fundamental es 0 > N .
La prueba es interesante pues aparece, como condicin necesaria y suficiente, el que una
de las especies sea ferminica, lo que (a la vista del resultado de Onsager) sugiere que el
principio de exclusin de Pauli juega el papel de un ncleo rgido.
Para ste y otros detalles ver E.H. Lieb, Rev. Mod. Phys. 48, 553 (1976).
41
Si el sistema no es elctricamente neutro, aparecen trminos superficiales adicionales
que dependen de la forma del recipiente.
96
Estabilidad en el lmite termodinmico.
97
sigue un calor especfico definido positivo:42
dQ
C= 0
dT
Un razonamiento anlogo conduce a una compresibilidad:
1 dV
K= 0
V dP
98
2G V
= = V KT 0 = G (T, P ) es cncava en P.
P 2 T P T
(59)
(1) (2) 1
Nk+1 = 4Nk +4Nk = 4n1Vk +4n2Vk = 4 (n1 + n2) Vk = (n1 + n2) Vk+1.
2
99
Este resultado se generaliza inmediatamente a un no arbitrario , m,
de divisiones del sistema:
!
Xm Xm
a
i ni i a
(ni) , (60)
i=1 i=1
P
con i nonegativos tal que i i = 1, es decir, a (n) es completa-
mente convexa.
100
mos en (60):46
n1 = n; ni = 0, i = 1
= a (n) a (n) , 0 1.
1 = (; i = 0, i = 1)
(61)
De la definicin de a, se tiene:
a
(n) = na (v) , (62)
luego
1 v
a (n) = na (v ) con v = = v, (63)
n
esto es,
(63) (61) (62)
na (v ) = a (n) a (n) = na (v) ,
de donde:
a (v ) a (v) con v v,
es decir, a es funcin montonamente nocreciente del volumen
especfico.
El sistema es estable.
En realidad, el demostrar que la P es ciertamente funcin continua de v (su ya se sigue
de los Ts generales) requiere algo ms de esfuerzo, pues la derivada de una funcin convexa
puede tener discontinuidades, saltos. Ruelle prob inicialmente continuidad absoluta para
sistemas clsicos con la hiptesis adicional de que | (r)| fuese acotado. Este resultado se
extendi luego a potenciales estables moderados.
101
Es sencillo mostrar que a (v, T ) es cncava en T, para concluir que
el calor especfico es nonegativo. En realidadD esto ya se Esigue de
la relacin disipacin fluctuacin kT 2CV = (E hEi)2 que se
demuestra luego. Por otra parte, relaciones como las (58) y (59) se
siguen de las propiedades de a (v, T ) mediante transformacin de
Legendre.48
102
Esto es importante, como veremos en el contexto de la teora de
Yang y Lee.
103
104
Tema 5: Fluctuaciones y equivalencia entre colec-
tividades
Teora de fluctuaciones. Frmula de Einstein. Divergencias.
Opalescencia crtica.
Equivalencia macroscpica de colectividades.
105
Los trabajos mximos que se pueden obtener no van
a ser iguales en las tres
En el caso de la microcannica, la E es fija, mientras
que en las otras puede fluctuar. Tampoco las fluctua-
ciones van a ser las mismas.
Sin embargo, cuando uno puede olvidarse de estas dife-
rencias, conducen a resultados idnticos. Por ejemplo, se
obtiene la misma fsica al estudiar el gas ideal en cada
uno de los formalismos. Se plantea entonces:
106
P P P P
U ( r eEr )( r Er2 eEr )( r Er eEr )( r Er eEr )
= P Er 2
( r e )
P 2 Er P
Er 2
E e E r e
2
= P
r r
+ P
r
= E + hEi2
re re
Er Er
49
De hecho, Er = const., r sin trabajo exterior. Ms concretamente, se vi:
X hnr i
P dV = dEr
r
para el trabajo medio realizado por uno de los miembros de la colectividad en un proceso
en el que se alteren los valores Er .
50
Reservo calor especfico para la capacidad calorfica por partcula, CV = (u/T )v ,
v = V /N, u = U/N, que es magnitud intensiva. Calor molar es la capacidad calorfica
por mol.
107
D E
2
Por otra parte, CV = (U/T )V es extensiva, luego (E) N
y se tiene para la medida relativa de las fluctuaciones en la energa:
rD E
2
(E)
N 1/2 1
=
hEi N N
D E
2
Es decir, aunque la fluctuacin cuadrtica media (E) es muy
grande (macroscpica), la amplitud de la fluctuacin comparada
con el valor medio es despreciable, por ejemplo, de un 109 % si
N 1022, y tiende a cero en el LT.
Se sigue que, en el LT, el formalismo cannico conduce, de hecho,
a una situacin semejante a la representada por una colectividad
microcannica caracterizada por una energa U = hEi . Aunque E
es una variable en la cannica, sus miembros se agrupan estrecha-
mente alrededor de hEi .
Con objeto de profundizar en esta propiedad, adoptamos otro punto
de vista. Supongamos que E es una variable continua (lo que es
buena aproximacin en el LT) y escribamos, como obtuvimos antes:
p (E) dE = eE (E) dE.
Es decir, la densidad de probabilidad de que el sistema dado, como
miembro de una colectividad cannica, tenga energa en intervalo
(E, dE) es producto del factor de Boltzmann, que decrece mon-
tonamente con E, y de la densidad de estados o degeneracin, que
crece montonamente con E. Puesto que p (E) es integrable y aco-
tada, ha de presentar un extremo para cierto valor de E, E , esto
es
E
E
E
e (E) E=E
= 0 = E
ln e E=E
=0
= E [E + ln ]E=E = 0
ln (E)
= .
E E=E
108
Es decir, el extremo E es precisamente el mximo U del que ha-
blbamos antes, independientemente de la naturaleza del sistema
dado. Sabemos adems que este mximo es extraordinariamente
pronunciado, es decir, la energa E = U es muchos (depende de
N ) rdenes de magnitud ms probable que cualquier otro valor, por
poco que difiera de U.
Con esta informacin, podemos investigar la forma de p (E) . Desa-
rrollamos p (E) , o su logaritmo, alrededor de la posicin U de ese
mximo muy pronunciado:
ln p (E)
ln p (E) = ln eU (U ) + E E=U (E U )
| {z }
= 0 por ser U un extremo
1 2 ln p (E)
+ 2 (E U )2 +
2! E E=U
1 2
ln p (E)
ln eU (U ) + (E U )2
2! E 2
E=U
Adems, se tiene
S
ln eU (U ) = U + ln (U ) = U + = (U T S)
k
y51
2 E 1 1
ln e (E) = ,
E 2 kT 2 CV
E=U
51
En efecto, es
E 1 k ln (E)
ln e (E) = ln eE + ln (E) = + ,
E E E k E
| {z }
S(E)
E
E 2
h i1
2
de donde E 2ln e (E) = k1 E S(E)
2 . Notando entonces
x
y = y
x y que
z z
la definicin T1 S(E)
E
1
permite definir T (E) S(E)
E tal que T (U ) T, tenemos
E=U
2
E n o
E 2 ln e (E) = k1 E 1
E=U h i1
T (E)
E=U
1 1 E
= k T (E)2 T (E) = kT1 2 C1V , QED
E=U
109
luego:
1
ln p (E) (U T S) 2
(E U )2
2kT CV
de donde
" #
2
(E U )
p (E) = eE (E) e(U T S) exp
2kT 2CV
110
Equivalencia microcannica-macrocannica: A. Fluctua-
ciones en el nmero de partculas.
Por lo que respecta a la colectividad macrocannica, su diferencia
con las otras es tambin ms conceptual que prctica. Aunque N
y E son ahora variables, sus miembros se apilan de hecho, muy
estrechamente, con distribuciones tipo delta de Dirac en el LT, al-
rededor de los valores medios. En consecuencia, los resultados son
termodinmicamente equivalentes en el LT.
y
2 X X X
2
N = Nr pr,s = Nr pr,s + Nr (Nr 1) pr,s
r,s r,s r,s
2
1 2
= z +z
z z 2 T,V
y52
D E
1
z
(N )2 N 2 hN i2 = z .
z | {zz }
hN i T,V
2
2
1
52
En efeceto, se comprueba z z z z = 1
z 2 z2 + z
z
1
2 z
z
111
z
Finalmente, notando que dz = d, tenemos
kT
D E
2 hN i
(N ) = kT . (64)
V,T
Este resultado es paralelo al cannico
D Ede Einstein.
D Muestra
E que la
fluctuacin cuadrtica media (N )2 = (N hN i)2 N, de
modo que la fluctuacin relativa es
rD E
2
(N )
1
,
hN i N N
confirmando lo anunciado.
112
Pero ms importante es su implicacin en las proximidades del pun-
to crtico de la transicin gaslquido, por ejemplo.
P
T = Tc
v
Las isotermas son entonces prcticamente horizontales, de modo
que:
[P/V ]T TC 0,
KT TC y
las fluctuaciones relativas en la densidad dejan de ser despre-
ciables (se hacen de orden unidad)
113
B. Fluctuacin macrocannica de la energa
Tenemos sucesivamente:56
X
1 ln
hEi = Es pr,s = =
r,s
z,V z,V
2 X 2 1 2
E = Es pr,s =
r,s
2 z,V
D E
2
2 2 ln hEi
(E) E 2 hEi = =
2 z,V z,V
y, dado que d = dT /kT 2 y hEi U,
D E
U
(E)2 = kT 2 O (N ) .
T z,V
56
Puede evitarse esta seccin, o dar slo el resultado, que es intuitivo.
114
Usando aqu la relacin termodinmica57
" #2
U 1 N U
= CV + ,
T z,V T T,V N T,V
se sigue
D E " #2
N U
(E)2 = kT 2CV + kT .
T,V N T,V
luego, en particular:
U U U N U U
(*) = + = CV + .
T z,V T N N T T z V N T,V T z,V
Adems,
S= A
T N,V
= S
=
N T,V T
= A N,V
=
N T,V
U
S
dU = T dS P dV + dN = N T,V =+T N T,V
U
(**) =T
N T,V T N,V
115
luego
" #2
D E D E D E U
2 2 2
(E) = (E) + (N ) .
macro canon macro N T,V
116
Equivalencia de las colectividades (incluso en el punto cr-
tico)
Sea, como ejemplo, sistema clsico de N partculas confinadas en
regin 3D, , de volumen V , con energa potencial de interac-
cin H (~q1, ..., ~qN ) . Recordemos que la parte configuracional de la
funcin de particin cannica puede escribirse:
Z
1
Q (N ) = d~q1...d~qN eH = 3N eV a (N )
N!
n
Z = 3N Q
pues
= kT
a V ln Z
(66)
y que se tiene en la macrocannica:
X
= z N 3N Q (N ) = eV P , z 0, (67)
N =0
a (n, T ) = lm a , P (, T ) = lm P . (68)
117
ne Q (N ) 0 para N > Nm , luego, combinando (66) y (67):
X
Nm X
Nm
N V a
(N )
e V P
= z e = eN eV a (N )
N =0 N =0
X
Nm X
Nm
(N V a
) N
= e = exp V a
(N ) .
V
N =0 N =0
Definimos:
N
P () max a
(N ) (69)
N V
V P
eV P eV P n
| m{zV} e ,
Nm
luego
1
P P P + ln (nm V )
V
P () sup {n a
(n)}
0 n nm
118
P () < n1 a (n1 ) + ,
1
n1 a (n1) P () . (70)
2
N
Para , tomando una secuencia N tal que V n1, por el
teorema del LT, tenemos de (69) y (70):
N
lm nf por debajo P () = lm (nf) max a (N )
N V
= max {n a (n)}
1
n1 a (n1 ) P () .
2
P () = supn [n a
(n)]
a () = sup [n P (n)] ,
120
Transformadas de Legendre
Una funcin convexa es una funcin g(x) tal que g (x) > 0. Se
define entonces una transformada de Legendre g(p) de g como
a
(n) = n P ()
n
P()
121
~
a(n) P ()
n
P()
al calcular la inversa
a
(n) = n P ()
P
se tiene que n = , no est definida en ese punto.
T
debemos en estos caso utilizar una transformada de Legendre ge-
neralizada
P () = sup {n a
(n)}
n
Cuando no hay un cambio de fase coincide con la definicin usual
movindonos por diferentes n el mximo de n a(n) ocurre para
el n tal que la recta n es tangente a la curva a(n)
~
a(n) ~
a(n)
n* n
n
~ *)
a(n
n*
P()
En los cambios de fase con la nueva definicin tenemos que todas las
pendientes comprendidas entre esas pendientes lmites son vlidas
122
Seminario 1: Teora de la Informacin.
Informacin, incertidumbre e ignorancia.
Reglas para asignar probabilidades.
Informacin y Fsica estadstica.
Otras aplicaciones.
Bibliografa:
Conceptos relevantes.
123
El concepto de incertidumbre se sigue de consideraciones puramen-
te estadsticas. Por ejemplo, sea una situacin en la que pueden
presentarse n sucesos distintos de un experimento aleatorio, todos
igualmente probables a priori, es decir, de acuerdo con nuestro
conocimiento inicial. Podemos pensar en
124
I. Postulado
A A
valores observados hAi
II. Postulado
X
= Tr(
hAi
A) = |ip h|
cte |ri/E Er E + E
E < H < E + E pr =
0 resto
125
Dado un experimento aleatorio se dice que se puede dar un
esquema cuando se da los sucesos y sus probabilidades:
A1 A2 . . . An Ai sucesos del experimento A
Esquema A
p1 p2 . . . pn pi probabilidades o frequencias del suceso Ai
A1 . . . An B1 . . . Bm
A= B=
p1 . . . pn q1 . . . qm
A 1 B1 A 1 B2 . . . A n Bm
AB=
p 1 q1 p 1 q2 . . . p n qm
126
Intuimos que
Vamos a verlo:
X
r
1 2 r
U (A A . . . A ) = U (Ai)
i
U (nr ) = rU (n)
nr ms nr+1
r ln n s ln m (r + 1) ln n
y dividiendo por s ln n
r ln m r 1
+
s ln n s s
127
Las incertidumbres para ese tipo de esquemas verifican
(por ser creciente)
d
[U (xn ) = nU (x)] U (xn ) nxn1 = nU (x)
dx
d
[U (xn ) = nU (x)] U (xn ) xn ln x = U (x) .
dn
Combinando estas dos ecuaciones:
U (x) 1 1
U (x) = U (x) x ln x = ln U (x) = .
U (x) x ln x n x ln x
Integrando entre x y a :
128
distintas. Si el experimento se repite muchas, N, esperamos que
el suceso i ocurra N pi veces. La nica incertidumbre adicional
que tenemos es el orden que aparecen los N pi resultados del
suceso Ai. Los sucesos de este experimento son cadenas del
tipo por ejemplo
A1 A1 Ar . . . An
El nmero de cadenas distintas (teniendo en cuenta que no
vale permutar sucesos iguales pues me dan la misma cadena)
es
N!
W = Qn
r=1 (N pr )!
Podemos asociar a cada una de estas cadenas un suceso de
hacer N experimentos y todos ellos tienen asociada la misma
probabilidad pues en todas ellas el suceso A1 aparece N p1 ve-
ces, el A2 aparece N p2, etc. Luego aplicamos nuestro resultado
anterior
N! N!
UN Qn = k ln Qn ,
r=1 (N pr )! r=1 (N pr )!
En consecuencia,
" #
X
n
UN (W ) = k ln N ! ln (N pr )!
r=1
se sigue
X
n
UN (W ) = N k pr ln pr
r=1
129
UN (W ) Xn
U (n) = k pr ln pr . (74)
N r=1
Sea un problema con m sucesos posibles, cada uno con dos posi-
bilidades, 0 1. El nmero total de posibilidades es n = 2m y, en
consecuencia, la incertidumbre asociada es
U = k ln n = km ln 2.
130
U =m = K = 1/ ln 2 = log2 e,
luego se sigue:
U = log2 n
en lugar de la definicin de Shanon arriba; se dice entonces que
la incertidumbre est medida en bits (ejemplo, en Brillouin 1962,
p.2).65
131
En teora de la informacin, esta ltima condicin equi-
vale a
maximizar nuestra falta de informacin (o incertidumbre)
acerca del suceso o del conjunto de sucesos en cuestin.
132
Aplicacin de Teora de la informacin a la FE
Ejemplos de aplicacin del principio de mxima ignorancia o in-
certidumbre:
Antes de aplicarlo a las colectividades veamos un ejemplo sencillo
de un dado para ver el formalismo: Disponemos de un dado con
caras numeradas del 1 al 6 que, al ser tirado sobre superficie plana,
puede mostrar cualquiera de sus caras. Qu probabilidades, pi ,
hemos de asignar a cada uno de los 6 sucesos posibles?
ln pi + 1 = 0 i,
de donde
pi = e1 = const.
133
El parmetro puede determinarse usando la condicin de
P 1
normalizacin, e = 6e1 = 1 = e1 = 1/6, luego
1
pi = e1 = .
6
El resultado es razonable, aunque trivial: en ausencia de otra
info, lo mejor que podemos hacer es asignar probabilidades 1/6
a todas las caras.
Supongamos que arrojamos el dado muchas veces y obtenemos i0
como media de los nmeros que aparecen el la cara superior,
es decir,
X6
hii = ipi = i0 .
i=1
En el caso anterior, ha de ser
1
i0 = (1 + 2 + 3 + 4 + 5 + 6) = 3 5.
6
Supongamos, sin embargo, que el experimento arroja i0 = 45,
lo que delata cierta inhomogeneidad en el peso de las caras
que favorece los nmeros altos (dado trucado).
Cmo escoger pi si contamos con esta nueva informacin?
Usando multiplicadores de nuevo, la situacin actual requiere
que: P
U = k 6i=1 pi (ln pi + 1) = 0
P
i0 = 6i=1 i pi = 0
P
1 = 6i=1 pi = 0,
de donde, multiplicando por 1, k y k, respectivamente, y
sumando:
X
6
k pi (ln pi + 1 i ) = 0.
i=1
134
donde los parmetros (equivalentemente, Z e+1 ) y se
siguen de normalizacin:
X X
pi = 1 = Z = ei
i i
7 1
i0 = + 3 coth (3) coth . (78)
2 2 2
Notemos que
135
Fsica estadstica: colectividades.
En definitiva:
Z
S k hEi + k 1 = k d (ln + 1 + H ) = 0,
136
ln + 1 + H = 0
de donde
1 H
= e1 eH =
e (80)
Z
R
con Z = e1 = d eH , que es la distribucin cannica.
Notas:
De la condicin de normalizacin:
Z
Z = d eH . (81)
ln Z 1
S = k ln Z k = k 2 ln Z,
que coincide con la entropa cannica si interpretamos
1/kT y A = kT ln Z.68
Puede concluirse lo mismo en el caso de un sistema cuntico,
donde hay que obtener el operador densidad a partir de la
definicin
S = kTr ( ln )
68
De hecho, puede verse que esta interpretacin es correcta siguiendo el mtodo de la FE.
Por ejemplo, para que tenga sentido la definicin
Z
hEi = d H eH ,
he de ser = 1/kT.
137
y la informacin, por ejemplo, de que:
hEi = Tr (H ) , 1 = Tr (
) .
de donde
1 H
= e1 e H = e
Z
d = dE d (E0) ,
se tiene:
Z
S+k 1 = k dE d (ln + 1 ) = 0 = = e1 = cte.
138
y la cte. se sigue de (82):
R
1/, E0 < E < E0 + dE, = dE d (E0 )
=
0, en otro caso
139
Tema 6: Gases ideales clsicos
Sistemas ideales.
Gas de Boltzmann: rango de validez.
Gas ideal termodinmico.
Estructura molecular.
Definicin.
El sistema ms sencillo, matemtica y conceptualmente, es el sis-
tema ideal (SI), es decir, una clase de sistemas caracterizados por
hamiltonianos de la forma
X
N
H= Hj , (1)
j=1
Los N conjuntos {j} , cada uno de ellos una partcula o una mol-
cula, tienen interseccin nula (no interactan entre s). Esto incluye
muchos casos de inters:
2
SI = gas molecular, cuando {j} incluye adems GDL como
rotaciones y vibraciones de la molcula poliatmica
SI = slido, si las partculas son osciladores armnicos, fo-
nones en un cristal, etc.
Por supuesto, {j} puede incluir otros GDL relacionados con estruc-
turas intraatmicas (si interesa estudiar fenmenos de excitacin o
ionizacin, etc.) o con el espn, o puede representar cuasipartculas,
como veremos.
Lo fundamental en un SI es la propiedad aditiva (1) del H, que
refleja el que representa un conjunto de partculas independientes
entre s.
En consecuencia,
1
Un SI abandonado en un estado arbitrario es incapaz de tender al estado de equilibrio (no
es ergdico); la interaccin entre componentes es indispensable para que exista relajacin
cintica y se llegue a una posicin o estado de equilibrio.
3
modelo lmite
modelo de referencia
4
Los niveles energticos r que definen el microestado de cada par-
tcula estn asociados a su energa cintica y son
h2 2 2 2 h2 2
r (nx, ny , nz ) = (n +n +n ) n nx , ny , nz = 1, 2, 3
2mL2 x y z 2mL2 r
Luego (N, E, V ) es el nmero de soluciones positivas distintas
independientes de la ecuacin
X
3N
2mV 2/3E
n2r = 2
E
r=1
h
que no es sino la ec. de una esfera 3N-dimensional de radio E en
el espacio de las fases. Luego (N, E, V ) es el conjunto de puntos
del espacio de las fases que estn justo en la superficie de esa es-
fera. Vamos a considerar que tenemos una imprecisin en E, luego
nuestro sistema puede estar con energa E y E+E con E E,
entonces
(N, E, V )
(N, E, V ) = (N, E + E, V ) (N, E, V ) = E
E
donde (N, E, V ) es el volumen de espacio de las
fases de la esfera
3N-dimensional (en sus ejes positivos) de radio E (es decir todos
aquellos microestados con E E)
pero el volumen de una esfera d-dimensional de radio R es (ejercicio)
d/2
Vd = (d/2)! Rd de donde se tiene
3N E
(N, E, V ) (N, E, V )
2 E
y se tiene para la entropa
" 3/2#
V 4mE 3
S(N, E, V ) = k ln (N, E, V ) N k ln 3 + Nk
h 3N 2
No es extensiva.
2
Aqui hemos aplicado la frmula de Stirling para N 1, ln N ! N ln N N. Por otra
parte los trminos delante de (N, E, V ) dan contribuciones despreciables ( ln N respecto
a N ) en N
5
Da lugar a la paradoja de Gibbs segun la cual la entropa de
mezcla de dos partes de un gas ideal clsico (que es un proceso
reversible) es distinta de cero
La paradoja se resuelve metiendo un factor de indistinguibili-
dad (N !) dividiendo en el nmero de microestados accesibles
(N, E, V ) y que da una expresin diferente como veremos de
la entropa (Ec. de Sackur-Tetrode)
Gas de Boltzmann.
X
N
E= (j) , (2)
j=1
6
del momento de la partcula, los nmeros cunticos de rotacin y
vibracin, etc.; lo que llamamos nivel).3
Se sigue que todas esas posibilidades hay que contarlas como una
y que, en consecuencia, la funcin de particin cannica es
" #
X X n !nk !... X
Z= exp (j) (4)
(1) (N )
N! j
m m
7
luego: (j) 1 = kT = (j) 1 para la mayor
parte de los niveles, de modo que las exponenciales correspon-
dientes en (4) son muy pequeas, y slo contribuyen a Z los
niveles ms bajos (slo stos son, de hecho, accesibles). Las
partculas van a ir acumulndose en los niveles ms bajos de
forma que va a influir el principio de exclusin de Pauli: hay
que distinguir entre fermiones y bosones y son las leyes cun-
ticas las que determinan la Z.
A temperaturas altas: el nmero de niveles accesible de hecho
(que contribuyen a Z) crece, dndose una gran probabilidad
de que las N partculas se distribuyan entre niveles distintos:
8
El modelo de partculas indistinguibles definido como aqul que est
caracterizado por la igualdad en (5) resulta tener notable relevancia;
lo llamaremos gas de Boltzmann.
Ha quedado claro que se trata de un caso lmite (comn)
de las gases de bosones o de los de fermiones (o de mezclas
de ambos) en el lmite de T s altas.
Notemos que, al ser los niveles distinguibles, y al estar
ocupados en promedio por una nica partcula, resulta
que las partculas son, de hecho, distinguibles (por el nivel
que ocupan) en este lmite de T s altas, como en clsica,
aunque la indistinguibilidad ha dejado traza en el factor
1/N ! en el sentido que cada partcula puede ocupar uno
de los niveles y hay N ! formas de hacerlo.
Propiedad fundamental de Z : factoriza como consecuencia de
que los trminos en (2) son independientes:
( )
1 YN X h i
Z= exp (j) .
N ! j=1 (j)
m
Como las partculas son idnticas, los N factores aqu son iguales,
luego
1 N
Z= Z
N! 1
donde X
Z1 = exp ([])
m
con suma sobre todos los niveles de una partcula.
En consecuencia, la energa libre de este sistema es:
kT (N ln N N ) N kT ln Z1 (T, V )
N
= kT N ln N kT ln Z1 (T, V )
e
= N kT ln e N 1 Z1 (T, V )
9
Funcin de particin cannica, propiedades termodinmi-
cas; gases perfectos
Como hemos visto antes para un gas de Boltzmann la funcin de
particin queda
1 N
Z= Z
N! 1
donde X
Z1 = exp ([])
m
con suma sobre todos los niveles de una partcula. La energa libre
de un gas ideal de partculas idnticas es entonces
kT (N ln N N ) N kT ln Z1 (T, V )
N
= kT N ln N kT ln Z1 (T, V )
e
= N kT ln e N 1 Z1 (T, V )
Al final el estudio se reduce a calcular Z1 .
La energa de una partcula en un nivel dado es suma de las aso-
ciadas con el mov de traslacin del CM y con los grados internos
de libertad (rotacin, vibracin, excitacin electrnica, etc),
= tr + in .
10
X
Z1 = Z1,tr Z1,in , Z1,tr(in) = etr(in)
m
y
XXX h2
2 2 2
Z1,tr = exp n + ny + nz .
nx ny nz
2mL2 x
11
En consecuencia, sin menoscabo de la descripcin cuntica, el mo-
mento es aqu una variable continua a todos los efectos prcti-
cos, y es lcito reemplazar las sumas sobre n por integrales so-
bre p teniendo en cuenta que, como consecuencia de (6), es
d p = (h/L) d n :
!
3 Z p 2
+ p 2
+ p 2
L x y z
Z1,tr = 3 d px d py d pz exp .
h 2mkT
y
XN
p2j
H=
j=1
2m
se tiene
Z P Y
N
1
p2j
Z cl = 3N e 2m
d3 ~qj d3 ~pj .
h N! j=1
12
Z N
V 2
Z cl = ep /2mkT
4p2 d p
h3N N ! 0
N
1 V 3/2 1
3 N
= (2mkT ) = V .
N ! h3 N!
Termodinmica del gas de Boltzmann
La energa libre es
e Z1 V 3
A (T, V, N ) = N kT ln = N kT ln e Z1,in
N N
donde, en general:
X
in
Z1,in = exp
kT
{in }
8
Las variables internas (amplitud de vibracin de un oscilador, etc.) caracterizan mol-
culas aisladas con posicin dada para su CM, luego toda variacin se restringe a regin del
tamao de la molcula, pero no de tamao V.
13
donde
ain (T ) = kT ln Z1,in (T ) .
A (T, V, N )
a (T, n) lm
N N
1
= kT kT ln + 3kT ln (T ) + ain (T )
n
N
= kT ln n kT + 3kT ln (T ) + ain (T ) , n= .
V
14
a
Entropa por partcula: s = T n
y se sigue de arriba que
21 T
s = k ln n + k 3k ln (T ) 3kT ain (T )
5
=k ln n 3 ln (T ) ain (T ) .
2
3
u = a + T s = kT + ain (T ) T ain (T )
2
3
cV = k T ain (T )
2
15
gas perfecto = gas de Boltzmann sin GDL inter-
nos. En termo no se discute la posibilidad de GDL
internos y se define el gas perfecto como aquel para
el que cv = 32 k
slo los GDL internos pueden introducir dependen-
cia en la T en el cV .
h i1
1 V p
Compresibilidad isoterma: KT = V P T = n n
T
Fluctuaciones. La frmual de Einstein relaciona las fluc-
tuaciones de la energa con el cV por lo tanto slo cuando
haya GDL internos habr dependencia con la T. Tambin
otras derivadas de potenciales como la suceptibilidad estn
relacionada con propiedades microscpica de la estructura
interna de las partculas del sistema.
16
17
Estructura molecular, grados de libertad internos conge-
lados y no congelados, rotacin, vibracin y excitaciones
electrnicas: clculo de Zin
18
Pero antes de estudiar separadamente esto efectos, veamos una...
Propiedad general de los GDL:Cada GDL est asociado a un
espectro de niveles accesibles a la partcula cuando el GDL est
activado. Separacin entre niveles energa consecutivos correspon-
dientes a distintos procesos dinmicos asociados con los GDL inter-
nos depende de los parmetros mencionados (If , k , j ) y puede
ocurrir, de acuerdo con el valor de la T. Sea kT que es una medidad
de la excitacin trmica de la partcula:
19
sta es la forma de la contribucin tpica de un GDL a CV (T ) .
Puesto que suele ser muy distinta para GDL diferentes, la forma
CV (T ) resultante de la superposicin de las contribuciones de todos
los GDL puede ser muy complicada.
Excitacin electrnica:
Para tener en cuenta este efecto, resolver ec. Schrodinger para elec-
trones en presencia ncleo positivo (y/o hacer espectroscopa), lo
que d conjunto de estados electrnicos {j } , que conduce para
cada tomo o molcula:12
X
Zelectr = ej .
estados
electrnicos
12
En condiciones extremas de T y P, como en algunos sistemas astrofsicos (estrellas
interiores o de neutrones, etc) es necesario contar con reacciones nucleares, transformaciones
de partculas elementales, etc. No consideramos aqu estas complicaciones; ver ej., R.L.
Bowers, J.A. Campbell & R.L. Zimmerman, Phys. Rev. D7, 2289 (1973).
13
Esta separacin es, al menos, un orden de magnitud mayor que la que entre niveles
de rotacin y vibracin, aunque hay alguna excepcin (ej., si el estado fundamental es un
multiplete como en el NO).
14
Ms precisamente, slo la degeneracin 0 del estado fundamental contribuye:
X X
Zelectr e0 = 1 = 0 ,
estados de estados de
energa cero energa cero
20
Rotacin:
~2
rot = j (j + 1) , j = 0, 1, 2, ...,
2I
En este caso,
~2
rot = ,
2Ik
que generalmente tiene valor relativamente pequeo (ver luego),
de modo que este GDL suele contribuir clsicamente, con 21 k, a T
ambiente.
21
Esta suma no es fcil de calcular en general, pero puede apro-
ximarse:
y = j (j + 1)
haciendo
d y = (2j + 1) d j
Z
rot T
= d y exp y = , T rot
0 T rot
(8)
1 1 1 1
+ f (0) f (0) + f (0) f v (0) + ...
2 12 720 30240
donde, haciendo f (j) = (2j + 1) exp Trot j (j + 1) se
sigue la frmula de Mulholland (para molculas heteronu-
cleares):
2
T 1 1 rot 4 rot
Zrot = + + + + ..., T & rot
rot 3 15 T 315 T
22
De aqu se tiene para el cV :
" 2 3 #
1 rot 16 rot
cV,rot =k 1+ + + ... .
45 T 945 T
23
modo que Z ha de dividirse por 2.
En consecuencia, puede introducirse no de simetra
( = 1 para heteronucleares y = 2 para homonucleares)
y tenemos
cl T
Zrot = , T rot
rot
Se sigue que:
cl T
arot = kT ln Zrot = kT ln
rot
arot T 1/ rot
srot = = k ln + kT
T rot T / rot
T
= k ln +k
rot
urot
cV,rot = = k.
T
Es decir, en esta aproximacin clsica (T rot ), cV tiene
el valor que predice el T. de equiparticin, 2 12 k (el H
clsico para las alteras tiene dos tnos. cuadrticos, si gira
entorno al eje z, pues la componente x y la y del momento
angular). Por otra parte, no afecta a cV , de modo q
no se manifiesta diferencia alguna entre molculas lineales
(simtricas o antisimtricas) en esta aproximacin.
Pero a T s bajas, no todas las orientaciones son posibles,
se acoplan los estados rotacional y nuclear, luego hay que
tener en cuenta el carcter ferminico o bosnico de las
partculas, y el tratamiento de las homonucleares o sim-
tricas no es tan sencillo (Ver Pathria).
24
Vibracin:
de donde
Y
X
Zvib = Zi , avib = ai
i=1 i=1
16
Notemos: Si n = nmero de tomos en molcula, se tiene:
3n 5, molculas lineales
=
3n 6, en otro caso
25
La T caracterstica para vibracin se define i = hi /k. Zi se
calcula luego (tratamiento sistemtico de osciladores armnicos).
El resultado es
cV,i = kE (i /T ) ,
x
xe
donde E es funcin de Einstein (E(x) = (x1) 2 ), tal que cV,i
1
2 2 k para T , como se corresponde con el hecho de que el
hamiltoniano clsico de osc. tiene 2 trminos cuadrticos.
26
I ha de ser promedio sobre perodo vibracin pero, si hay
anarmonicidades, r0 crece con energa de vibracin.
27
Sistemas ideales en la macrocannica.
Z1 (V.T ) = V f (T )
Se sigue:
kT
P = ln = zkT f (T ) (9)
V
N = z ln = zV f (T ) (10)
z
17
Pathria p.104108.
28
y, combinando:
P V = N kT
es decir, ecuacin de estado es independiente de for-
ma de f (T ) , y concepto de gas ideal es muy general (ya lo
obtuvimos este resultado tambin en la microcannica)
De pasada, notamos otras propiedades de este modelo:
Notando U = ln y / = kT 2/T, y elimi-
nando luego z con (10), se tiene:
U = +kT 2 [zV f (T )] = zV kT 2 f (T )
T
2 f (T )
= N kT
f (T )
de donde
U f 2T f + f T 2f T 2 f 2
CV = = Nk
T f2
2T f f + T 2 f f f 2
= NK
f2
de modo que la forma de f (T ) es esencial para CV .
En casos sencillos, es f (T ) T , luego:
U = (N kT ) , CV = (N k) .
Tambin
U 1U
P V = N kT = = P =
V
es decir, la P es proporcional a la densidad de energa,
con cte. de proporcionalidad 1/, lo que es bien conocido
para el gas ideal norelativista con = 3/2 y para el gas
relativista extremo con = 3.
29
La funcin de particin es
X
X
(V, T, z) = N
z ZN (V, T ) = [z (T )]N [1 z (T )]1 ,
N =0 N =0
1 z
N =z ln = z ln = z ln (1 z) =
z z 1 z z 1 z
y, despejando z :
N N 1 1
z (T ) = =1 N 2 N 2 1
N +1 N N
19
Problema: Comprueben cmo se obtienen los resultados para un sistema de osciladores
cunticos en una dimensin si se hace aqu:
1
~
(T ) = 2 senh ,
2kT
y que se sigue el correspondiente resultado clsico haciendo:
kT
(T ) = .
~
Discutan este hecho.
30
luego
1
1 z (T )
N
y
1 1 1
z 1 . (11)
N
Equilibrio vaporslido:
kT kT ln N
P = ln = ln (1 z) kT
V V V
es decir, la Pslido 0 en el LT, como habra que esperar. Por otra parte,
zkT 2 (T )
U = kT 2 ln = kT 2 ln (1 z) = .
T T 1 z (T )
U (T )
kT 2 .
N (T )
En la misma aproximacin,
A ln N
kT ln (T ) + O
N N
S (T ) ln N
ln (T ) + T +O .
Nk (T ) N
31
Para el vapor, de (10):
Nvapor
zvapor = .
Vvapor f (T )
Luego, igualando:
Nvapor f (T )
= . (12)
Vvapor (T )
En consecuencia, si N = Nvapor +Nslido, la condicin necesaria
para la formacin de fase slida a partir del vapor es que:
f (T ) f (T )
N Nvapor = Vvapor =V ,
(T ) (T )
pues V Vvapor cuando empieza a formarse la nueva fase, es
decir:
N f (T )
.
V (T )
Si definimos la temperature de transicin, Ttrans. mediante:
f (Ttrans.) N
= ,
(Ttrans.) V
la condicin anterior puede expresarse:21
T Ttrans.(coexistencia de fases).
21
Si, como veremos con los teoremas de Yang-Lee, f (T ) y (T ) cumplen ciertas condi-
ciones matemticas que slo pueden demostrarse en presencia de interacciones.
32
33
Tema 7: Slidos
Pequeas oscilaciones. Modos normales. Fonones.
Osciladores armnicos, clsicos y cunticos, indepen-
dientes.
Calores especficos. Discusin. Ley DulongPetit.
Modelo continuo de Debye: discusin y variantes.
Modelo discreto de BornKarman.
34
Molculas de un slido no pueden moverse libremente,
como las del gas, sino que permanecen vibrando en
los nudos de una red.
X3N
1 2
+ (xi xi ) (xj xj ) + ...
i,j=1
2 xixj {xi }={
xi }
donde
24
Innecesario considerar GDL internos: como las FP factorizan, vale lo dicho para gas de
Boltzmann. Tambin olvidamos electrones libres, que estudiamos por separado.
35
0 ({xi}) es la mnima energa del slido (es cte.,
luego puede omitirse cambiando el cero de energa),
2o trmino se anula pues hay mnimo en {xi } = {
xi }
siguiente trmino representa pues las 1as contribuciones
no-nulas de vibraciones atmicas.
X
3N X
3N
1 1 2
H= mi2 + i,j i j , i,j .
i=1
2 i,j=1
2 xixj {xi }={
xi }
36
Fonones:
37
Sistema de osciladores independientes: caso clsico. Ley
de Doulong y Petit
En definitiva, teora de radiacin del cuerpo negro y teora de vi-
braciones en red pueden iniciarse siguiendo tratamiento paralelo a
gas ideal de Boltzmann, pero hay dos posibilidades de estudiarlos,
es decir, suponer al sistema como:
un conjunto de osciladores
una regin acotada que contiene un gas de los cuantos corres-
pondientes (fotones o fonones).
Estos cuantos son bosones en los dos casos, luego el 2o punto de vista
se desarrollar ms adelante. Desarrollamos ahora el 1er punto de
vista, separando tratamientos clsico y cuntico.
X
N
H= Hi
i=1
con
1 1 2
Hi = m 2 qi2 + p,
2 2m i
38
suponiendo de momento la misma frecuencia para todos los os-
ciladores. Se sigue, suponiendo que los elementos son en realidad
distinguibles entre s y = 1 (unidimensional):
QN () = [Q1 ()]N
con
Z
1
Q1 () = d qi d pi eHi .
h
Z Z
1 1 2 2 1 2
Q1 () = dq d p exp m q + p
h 2 2m
Z Z
1 1 2 2 1 2
= d q exp m q d p exp p
h 2 2m
1/2 1/2
1 2 2m 1 2 1
= = =
h m 2 h ~
luego:
QN () = (~)N
y
~
A kT ln QN () = N kT ln .
kT
26
Notamos que se permite que elongaciones y momentos tomen cualquier valor, aunque,
dado que las distribuciones son gaussianas, apenas contribuirn los extremos de los intervalos
de integracin. Por otra parte, en la integracin se usa
Z r
ax2
dx e =
a
39
Se sigue inmediatamente que la presin es
A
P = =0
V
A
S=
T
1
~ ~ ~k
= N k ln N kT 2 2
kT kT k T
| {z }
+N k
kT
= N k ln +1 ;
~
la energa interna es
U = A + T S = N kT ;
el capor especfico es
U U
CV = = = CP = N k.
T N,V T N,P
27
si las hubiese extendido a un cierto V habra una dependencia en V que desaparecera
en el LT.
28
Por ser P = 0, coinciden energa libre de Helmholtz y de Gibbs o entalpa libre, G
A A ~
N = A + P V, tenindose para el potencial qumico: = N = N = kT ln kT .
40
Problema: Conocida QN , calclese la densidad de esta-
dos (E) . Comprese las termodinmicas que se siguen
de QN y de (E) , discutiendo el resultado.29
con camino de integracin (que puede deformarse en tanto en cuanto la integral converja)
paralelo y a la derecha del eje imaginario, es decir, a lo largo de la lnea Re > 0. Pero
Z s +i
1 esx xn
I= d s =
2i s i sn+1 (n + 1)
S (N, E) = k ln (E)
E
k N ln ln E (N 1) ln (N 1) + (N 1)
~ | {z } | {z }
N ln N +ln N N
E
k N ln ln E + ln N + N
N ~
E
= N k ln +1 +k (ln N ln E)
N ~ | {z }
E N, luego este trmino es
despreciable frente anterior
E
N k ln +1
N ~
1 E
Tambin se tiene para la T del sistema: T = = Nk .
( E
S
)N
41
calor especfico por tomo = 3k
energa por mol (R = NAvogrado k): U = 3RT
calor especfico por mol (R 2 cal/mol grado): CV = 3R 6
c/m g.
42
43
Correcciones al modelo de Boltzmann. Sistema de oscila-
dores independientes: caso cuntico.
X X
r 1
Q1 () = e = exp r + ~
r r=0
2
X
ex/2 1
x/2 rx
=e e = =
r=0
1 ex ex/2 ex/2
1
1
= 2 senh ~
2
N N
1 exp ~
QN () = [Q1 ()]N = 2 senh ~ = 2
.
2 [1 exp (~)]N
44
1
A kT ln QN () = N kT ln 2 senh ~
2
1
= N ~ + kT ln (1 exp (~))
2
La entropa es
A ~
S= = Nk ln (1 exp (~))
T exp (~) 1
La energa interna es
U = A + TS =
1 ~
N ~ +
2 exp (~) 1
U U
CV = = = CP =
T N,V T N,P
exp (~)
N k (~)2
[exp (~) 1]2
45
clsicamente, U/N ~ es igual a kT /~, mientras que
cunticamente, U es siempre superior a kT :
46
Esto coincide globalmente con las observaciones experimenta-
les, aunque hay todava ciertos problemas tericos que discu-
timos ms adelante.
47
48
Problema: Clculo y discusin de la funcin densidad de estados
(E) para un conjunto de osciladores cunticos.
X X
N N +r1 r (N + r 1)! r
(1 x) = x = x,
r=0
r r=0
r! (N 1)!
tenemos:
N N
QN = 1 e~ e 2 ~
X
N + r 1 ( N2 ~+r~)
= e .
r=0
r
X
N +r1 1
(E) = E r + N ~ .
r=0
r 2
49
Problema: llegar al mismo resultado suponiendo que hay que dis-
tribuir r bolas indistinguibles en N cajas distinguibles de modo que
cada caja reciba 0,1,2,... bolas.
Fermi = , Planck = + /kT .
e/kT +1 2 e 1
50
La nica diferencia esencial es el signo del denominador, lo que se
corresponde con la diferencia formal entre las estadsticas cunticas
de FermiDirac y BoseEinstein, respectivamente, como veremos.
51
Modelos: Einstein versus Debye.
donde:
y el calor especfico
U X (~i)2 e~i
3N
CV (T ) = =k . (15)
T (e ~i 1)2
i=1
de donde se sigue que el potencial qumico, = (A/N )V,T de un sistema de fonones (como
en el caso de fotones, hacerlo como ejercicio) es idnticamente nulo.
52
como supuso Einstein por sencillez. Escribiendo la temperatura
caracterstica de Einstein, E ~E /k, se tiene
E
CV (T ) = 3N k E , (16)
T
donde se ha introducido la funcin de Einstein:
x2 ex
E (x) .
(ex 1)2
Se sigue de aqu lo discutido antes grficamente,32 que (tomando
E como parmetro fenomenolgico) describe las observaciones ex-
perimentales:
33
De hecho, este modelo dej claro en la poca cmo hay GDL (congelado) que no se
excitan hasta que sistema alcanza T caracterstica tal que kT sea comparable a separacin
entre los niveles correspondientes. Estas ideas fueron aplicadas enseguida a los GDL internos,
como se vi.
53
predicho (es decir va como una potencia T con 3).
La solucin no tard:
2
g () d = V 2 3 d ,
2 c
de donde
1
3 N 1 2
D = 18 2 + . (17)
V c3L c3T
54
En definitiva, Debye propone el espectro:
2
V 1 2 9N 2
g () = 2 c 3 + c3T
= 3 , si D
2 L D
0, si > D .
Z D
9N 2
CV = d 3 kE (x) , x ~
0 D
Z D
9N k x2 ex
2
3 d
D 0 (ex 1)2
= 3N kD (x0) (18)
donde
x0 ~D D /T,
D es la temperatura de Debye, y
se ha introducido la funcin de Debye:
Z x0
3 x4 ex d x
D (x0) 3
x0 0 (ex 1)2
Z x0
3x0 12 x3 d x
= x0 + 3 .
e 1 x0 0 ex 1
55
La situacin se resume en la figura:
56
Heat capacity of solids according to the Debye theory. The observed heat ca-
pacity of cooper is compared with the theory for a value of the characteristic
temperature = 308 K.
Pb Ag Zn Cu Al ClNa ClK
D de CV 88 215 308 345 398 308 233
D de cL , cT 73 214 305 332 402 320 240
donde las discrepancias son del orden del error experimental.
En consecuencia, el modelo ha encontrado innumerables aplicacio-
nes. Pero este xito fenomenolgico no debe ocultar una dificultad
bsica, cualitativa del modelo. Por ejemplo se ha supuesto que D
no tiene dependencia en la temperatura lo que supone que las cL
57
y cT no dependen de la temperatura, lo que no es cierto. Se pone
en evidencia, por ejemplo, al analizar la cuantizacin de las ondas
elsticas en experimentos de dispersin inelstica de rayos X y de
neutrones por slidos. Uno de estos experimentos mostr que g ()
tena la forma:
The normal-mode frequency distribution g() for aluminium. The solid curve
is derived from x-ray scattering measurements (Walker, 1956) while the dashed
curve represents the corresponding Debye approximation.
34
Notad que como max = 2c/min pero el valor minimo es el mismo para los modos
transversales y longitudinales L = T = ( 34 V /N )1/3
59
Modelo de Born von Karman
35
Evito muchos detalles. Vase Kittel p.168, Dekker p.47, Gopal p.147 y, un tratamiento
ms general, en L. Brillouin, Wave Propagation in Periodic Structures, Dover, N.Y. 1953.
36
No se determina el signo resultante, que depender del de los desplazamientos netos.
60
Estamos interesados en soluciones con forma de onda propa-
gndose, luego:
2
n (t) = ei(tqnd) , q= ,
donde
q = vector de onda, y
nd = posicin de equilibrio del tomo n respecto del ori-
gen.
61
correspondiente a las lineas rayadas en la figura (const. veloc.
propag. onda).
38
Dada la periodicidad de en q, el sistema tiene el mismo estado de vibracin para el
vector de onda q que para q 2s/d, con s = 1, 2, ... Con objeto de tener una relacin
nica entre y q, hay que limitar rango de variacin de q. Generalmente, se toma
q ,
d d
que se conoce como primera zona de Brillouin; la 2a est formada por los siguientes dos
medios perodos, uno a cada lado de q = 0, y as sucesivamente.
62
Otros hechos sobresalientes del modelo Born von Kar-
man
Condiciones peridicas en los lmites:
que implica
q = 1, 2, ..., N/2.
En consecuencia:
63
usando la representacin de ondas estacionarias 39 , en la que
0 q d y d g = (L/) d q, llevando a la expresin de CV :40
Z
L /d x2 ex ~
CV (T ) = k d q, x = .
0 (ex 1)2 kT
Falta figura
tomos distintos:
2 m = 2u cos qd 2u
2M = 2u cos qd 2u,
de donde
" 2 #1/2
m+M m+M 4 sin2 (qd)
2 = u u
mM mM mM
39
Para una onda estacionaria de g nodos (g = 1, 2, . . .) en una cuerda de longitud L su
g = 2L/g luego q= L g >0.
40
= mx sen qd 2 = q = 2
d arcsin
mx , luego d q = dd q d =
2r d 1 2 d
d 2 mx = d ( 2 2 )1/2 .
1 mx
mx
64
y tenemos dos posibles relaciones de dispersin
con r r r
2u 2u 2u (M + m)
1 = , 2 = , 3 = .
M m Mm
Es decir:
si q 0, hay dos soluciones, 4 q y 3 const
usando 4 en las soluciones 2n y 2n+1 se tiene , es decir,
los tomos se mueven conjuntamente, como en las ondas de
1/2
sonido ordinarias, con un velocidad c = 2ud2/ (M + m) .
A esta rama se la llama en consecuencia acstica, y se habla
de fonones acsticos.
usando 3, se llega a (M/m) , de modo que los tomos
adyacentes se mueven en sentido contrario, manteniendo fijo su
CDM. Si los dos tomos tienen cargas opuestas, como ClNa, se
origina dipolo fluctuante que puede interaccionar fuertemen-
te con campo electromagntico; se hable en consecuencia de
fonones pticos.
Con este modelo la distribucin g () muestra dos picos, de
acuerdo con observaciones experimentales.41
41
Gopal p.147152; tambin otras propiedades, p.152155.
65
Tema 8: Fenmenos magnticos I
En esta tema aplicamos de nuevo el concepto de sistema ideal para
empezar a comprender las propiedades magnticas de la materia.
Contenido principal:
42
De hecho, la masa del ncleo % 1000 veces la masa electrnica, y dado que el momento
e ~
magntico es ~ = g 2mc J, luego slo contribuye en la prctia a travs de su influencia en las
de los electrones y en la formacin de la del in.
66
Materiales diamagnticos:
negativa, < 0, y
pequea, tpicamente 106 / mol
depende poco de la T (salvo a temperaturas muy bajas,
ver ms abajo); por ejemplo:
T,P 0T,P 1 + (T T0) T P
dilatacin compresibilidad
Ejemplos:
67
A temperaturas muy bajas T s, los metales diamagnticos pre-
sentan fuertes variaciones de con T, y violentas oscilaciones
~ (cuyo periodo est relacionado con
al variar poquito valor de H
la energa de Fermi) que se llama efecto Hass van Alphen
que con otras propiedades diamagnticas estudiamos ms
adelante.
Materiales paramagnticos:
para nometales:
const 1
T,P , con v v0 [1 + (T T0) T P ] ; (, T 1)
v T
luego, en la prctica, se observa ley de CurieWeiss:
1
T,P
T
donde = Temperatura de Curie asociada al material. A
T s suficientemente altas, T , se tiene en la prctica
la ley de Curie 1/T (que la cumplen todos los
materiales paramagnticos metales o no a temperaturas
suficientemente altas).
los metales tienen un comportamiento ms complicado con
T.
Ejemplos:
68
Los tomos (Al, SO4Cu) tienen momento magntico perma-
nente y, en consecuencia, adems de efecto diamagntico, siem-
pre presente (aunque muy dbil generalmente), el momento
magntico puede interaccionar con el campo aplicado (sobre
todo si ste es muy intenso), lo que produce el efecto contrario
(esto es paramagnetismo y es ms intenso).43
Todos los materiales paramagnticos sufren cambio de fase a
T s suficientemente bajas (el paramagnetismo no es estable a
bajas temperaturas, es una fase desordenada, S 6= 0):
44
De hecho, como consecuencia de la ley de CurieWeiss:
1
T
2 9 0,
(T ) T 0
de modo que la de una fase paramagntica a T = 0 implicara la violacin del 3er Principio.
69
vidrios de espines y los sistemas con campos aleatorios. En
todos los casos estas fases son consecuencia de la interaccin
entre momentos magnticos permanentes (que son las relavan-
tes a temperaturas bajas)y no pueden estudiarse con un H
de sistema ideal, sin interacciones. A temperaturas altas, la
energa debido a la agitacin trmica es ms intensa y estos
materiales se comportan como paramagnticos.
J(r)
0
r
70
71
Dipolos en campo magntico.Teorema de Bohr van Leeu-
wen. Teora semiclsica de paramagnetismo de Langevin
Modelamos
P ahora cristal magntico como sistema ideal, es decir,
H = i Hi :45
45
El modelo slo tendr sentido cuando los elementos (ej, iones) magnticamente activos
estn suficientemente separados. Esta situacin es realizable experimentalmente; por ejem-
plo, diluyendo magnticamente sales de tierras raras, de modo que los iones magnticos
quedan en una subred separados por iones magnticamente inactivos. El resultado es prc-
ticamente un sistema ideal dado que las interacciones relevantes son de corto alcance y, a
distancias moderadas, despreciables frente a la agitacin trmica, que es notable a las T s
para las que este modelo es interesante. A bajas T s se hacen importantes las interacciones
frente a la agitacin trmica y se produce el cambio de fase a estados ferro y antiferromagn-
ticos, o a vidrio de espin, que no puede ser explicado con un H de elementos independientes;
las interacciones son esenciales para explicar este fenmeno.
72
En situaciones intermedia el campo consigue magnetizacin
parcial; la situacin est controlada entonces por alguna fun-
cin de la temperatura.
Respecto del H :
aparte de la agitacin trmica, los dipolos son prctica-
mente estticos, de modo que H no contiene trminos ci-
nticos46,
la energa potencial, debida a interaccin dipolos con cam-
po externo es, en buena aproximacin:47
X
N X
N
H= ~ = H
~i H cos i
i=1 i=1
XY
N
exp (H cos i) (22)
{i } i=1
73
P Q
Dado que los dipolos son independientes, conmutan y y
se tiene:
QN () = [Q1 ()]N
donde X
Q1 () = exp (H cos )
es la funcin de particin de cualquiera de los dipolos.
X
cos i eHi
X
N
{i }
= X .
i=1 eHi
{i }
P Q
Permutando y en (22), tomando logaritmo y derivando:
X N X
ln QN () = ln exp (H cos i)
H i=1
H
{i }
X
cos i eHi
X
N
{i }
= X = Mz ,
i=1 eHi
{i }
esto es,
74
1 N
Mz = ln QN () = ln Q1 () = N h cos i
H H
A
= pues A = kT ln QN ()
H
~ =
Un campo magntico H ~ A
~ aplicado al electrn transforma
su momento:
e~
p~j p~j A (~qj ) ,
c
75
~q ~q
e~
~p A (~qj ) ~k
c
1
Mz = ln QN () = 0.
H
QED
independiente de t.
49
q d p~ = |J| d ~q d ~k, se tiene
En efecto, para d ~
1 0 0 0 0 0
0 1 0 0 0 0
0 0 1 0 0 0
|J| = Ax Ax Ax 1 0 0 =1
x y z
Ay A A
x yy zy 0 1 0
Az
x A y
z
Az
z 0 0 1
(Ntese que todos los elementos por encima de la diagonal principal son nulos.) En todo
caso, la transformacin es cannica y el jacobiano de toda transformacin cannica es 1.
76
Dipolos: teora (semiclsica) de Langevin
Mz 1 1 1 Q1
z = ln Q1 () =
N H Q1 H
1
= coth (H) = L (H) ,
H
H
x
kT
Mz N 2 1
lm = (ley de Curie)
H0 H 3kT T
77
Dipolos: tratamiento cuntico
|~|z z = gB m
donde:
e~
B = es el magnetn de Bohr,
2me c
3 s (s + 1) ( + 1)
g= 2 + (j = + s)
2j (j + 1)
m es el no cuntico magntico: m = j, j + 1, ..., j 1, j
X
= exp (Hz ) .
posibles z
50
El momento magntico ~, tanto si el momento angular total J~ es debido por completo a
momento angular de espn, J~ = S,
~ a momento angular orbital, J~ = L,~ o a suma, J~ = L
~ + S,
~
satisface
e ~
~ =g
J,
2me c
donde el factor de Land (o giromagntico) tiene distinto valor en cada caso, g = 2 para
J~ = S
~ y g = 1 para J~ = L~ y (Ntese que J 2 = j(j + 1)~2 ; donde j = 1 , 3 , 5 , ... 0, 1, 2, ...
2 2 2
respectivamente).
Detalles en Pathria p.82, Kittel p.504 , mejor, EisbergResnick (Fsica cuntica, p.427-
429 de la edicin espaola, o bien sec. 10.6 sobre efecto Zeeman).
78
Si escribimos:
= max z = gB j
x H = HgB j,
de modo que
x
Hz = HgB m = m,
j
se tiene51
h i
X j 1
senh x 1 + 2j
mx
Q1 () = exp = ,
m=j
j x
senh 2j
luego
Mz 1 1 1 Q1
z = ln Q1 () =
N H Q1 H
1 1 1 x
= gB j 1+ coth x 1 + coth
2j 2j 2j 2j
= Bj (x) ,
donde
1 1 1 x
Bj (x) 1+ coth x 1 + coth
2j 2j 2j 2j
rn 1
51
Es progresin geomtrica, a1 + a1 r + a1 r2 + ... + a1 rn1 = Sn = a1 , con a1 = ex ,
r1
r = ex/j , n = 2j + 1, de modo que:
ex(2j+1)/j 1 ex(2j+1)/j 1 ex/j ex(2j+1)/2j ex(2j+1)/2j
S = ex x/j
= ex x/j
x/j = .
e 1 e 1 e ex/2j ex/2j
79
es la funcin de Brillouin de orden j:
Casos particulares:
3
H dbil y T alta (x 1) : escoth x = x1 + x3 x45 + ...
(0 < |x| < ), luego Bj (x) x3 1 + 1j (lneas discontinuas
en figura) y
x j + 1 2 1
z = 1+ H.
3 j 3kT j
80
N02 1
Mz = N0
z 1+ H
3kT j
Mz N02 1 Cj
lm 1+ =
H0 H 3kT j T
1/ versus T para una sal de gadolinio: Gd(C2 H5 SO4 )9H2 O; la lnea recta es la
ley de Curie.
N0 2 N0 g2 2B j(j+1)
53
Cj 3k 1 + 1j = 3k
81
1/ versus T para una muestra en polvo de copperpotasium sulphate hexahydrate
z /B versus H/T. Las curvas son la prediccin experimental. Los smbolos son datos para,
82
Lmite clsico.
e~
= gB j = g j ~j
2me c
1 1
1+ coth x 1 + coth (x)
2j 2j
1 x 1 2j x x3 1
coth + 3
+ ...
2j 2j 2j x 6j 8 45j x
1
lm Bj (x) = coth (x) = L (x)
j x
83
Campo inducido.
X
N
H = H cos i.
i=1
B = H = (1 + 4) H
se tiene
N0 2 1 1
= 1+ ,
3kT j 300T
luego
105 (no tiene dimensiones) a T ambiente, y
103 102, a muy bajas T s, ej, T 1o K.
En definitiva, el error por haber tomado H en lugar de B = (1 + 4) H
es ciertamente despreciable para sustancias paramagnticas, parti-
cularmente a T alta, que es cuando esta teora es de aplicacin.
84
85
Sistemas con temperaturas negativas.
TERMODINMICA:
S = S (U, X)
Vamos a ver que, sin embargo, unos pocos sistemas fsicos (sistemas
anormales) que violan este ltimo postulado, de modo que tienen
S capaz de disminuir con U y, en consecuencia, pueden presentar
T s negativas.
86
3. est esencialmente aislado de otros que no satisfagan estas con-
diciones.
54
Es cierto que las velocidades no pueden superar c, pero esta limitacin no afecta al
momento que aparece en el H, que no est acotado.
87
Esto se satisface tambin en algunos cristales, como F Li, para el
que:
espines
tenemos dos subsistemas
fonones (vibraciones en la red)
asociado interaccin espnfonn 5 minutos
asociado interaccin espnespn del orden de la precesin de
larmor de un ncleo en el campo de su vecino 105 segundos
88
y producen una contribucin neta (de la magnitud que sea):
U = N+ N .
con N es nmero de elementos con energa . Claramente este
sistema cumple la condicin 2, pues al tener slo dos estados ener-
gticos, la energa total est acotada. Nos preguntamos por el no
total de maneras de conseguir U, es decir, la degeneracin del
macroestado (N, U ) ; es claro que:
N!
(N, U ) = .
1 U 1 U
N+ ! N !
2 2
| {z } | {z }
N+ N
S = k ln (N, U )
1 U 1 U 1 U
k N ln N N N+ ln N+ + N+
2 2 2
1 U 1 U 1 U
N ln N + N
2 2 2
1 U 1 U 1 U 1 U
= k N ln N N+ ln N+ N ln N
2 2 2 2
y, derivando:
U
1 S k N
= = ln .
T U N 2 U
N+
Esto es,
89
U/
T > 0 : si N
N +U/ > 1, implicando U < 0, que ser el caso ms
frecuente (N > N+ : hay ms elementos en el estado que
en el +)
55
Acorde con hecho general: = N para sistemas fsicos con N constituyentes, cada uno
capaz de (valores) estados equiprobables, de modo que su entropa mxima es S = N k ln .
90
Tambin se tiene:
91
Termodinmica de sistema de espines nucleares.
" #N
X
QN () = [Q1 ()]N = eH
estados
= e e = [2 cosh ()]N
A n h i o
S= = N k ln 2 cosh tanh
T kT kT kT
92
Tambin podemos notar que el sistema tiene dos estados:
+ p2
p1
e e
p1 = , p2 = ,
Q1 Q1
y nmeros de ocupacin N p1 y N p2 .
el calor especfico:
U 2
2
CH = = Nk sech .
T kT kT
y la magnetizacin:56
M = N p1 N p2 = N tanh .
kT
S = N k (p1 ln p1 + p1 ln p1)
56
Ntese cmo estas expresiones implican U = HM, que corresponde a un sistema ideal.
Tambin se evidencia as otra vez cmo HM hace aqu el papel del trmino P V en fluidos.
93
Grficamente:57
94
simple monocomponente) cuando se mantienen constantes su
volumen y nmero molar:
Discusin.
58
Vase W.G. Proctor, en Scientific American Vol.239, No.2, pp.9099, August 1978.
95
Usando resonancia magntica nuclear pueden prepararse a T s
distintas, eventualmente , que se determinan con precisin
midiendo polarizacin.
+ p2 e/kT
p1 e/kT
96
aumento de energa en un sistema normal no acotado hace que
partculas se deslicen a niveles ms altos, aumentando T
sin lmite.
60
Problema: Sea una mquina termodinmica funcionando entre fuentes de calor a T1 y
T2 , ambas negativas, |T1 | < |T2 | , que toma calor Q2 de la fuente ms fra, convierte parte en
trabajo W, y cede el resto Q1 a la otra fuente. Convencerse de que nada hay en la Termo que
impida este funcionamiento y de que se consigue as el mvil perptuo de segunda especie.
Discutid esta aparente paradoja. (Libro de Biel.) De hecho no contradice el 2 principio pues
para mantener a los dos sistemas a temperaturas negativas hay que realizar un gran trabajo
sobre los espines.
97
imaginemos estado +H relativamente ms poblado que H
(que es el estable, en definitiva);
98
Tema 9: Gases ideales Cunticos
Indistinguibilidad cuntica: Conexin Espn-Estadstica.
Sistema ideal de bosones y fermiones. Funcin de parti-
cin cannica y macrocannica
Introduccin
99
(. . . , rk , . . . , rl , . . .) = (. . . , rl , . . . , rk , . . .) , = 1,
| (. . . rk . . . rl . . .)|2 = | (. . . rl . . . rk . . .)|2 ,
100
hay restricciones pero, pero el PP = que no puede haber dos fer-
miones en el sistema con los mismos nmeros cunticos (no se puede
construir dicho estado62. En consecuencia, es ms apropiado usar
el formalismo de segunda cuantizacin o de nmeros de ocupacin,
que incorpora naturalmente esta propiedad de los elementos.
con {mi (ri )} base completa para sistema de una partcula don-
de
c (. . . mk . . . ml . . .) = c (. . . ml . . . mk . . .) k, l,
de modo que
c (. . . mk . . . ml . . .) = 0 si mk = ml para = 1
101
El operador hamiltoniano para un sistema ideal se puede poner en
formalismo de segunda quantizacin como63
XX XX
H = p~ a+ a = p~ Np~,
p~, p~,
p~ p~
63
Aqu ap+ ~, ) es el operador creacin (destruccin) de una partcula de momento p
~, (ap ~y
+ + + + +
espn y son tal que [a , a] =1, a = (a ) y a = (a) , y actan de la forma
a|ni = n|n 1i; a+ |ni = n + 1|n + 1i; a+ a|ni = n|ni.
luego Np~, = a+ ap~, es el operador nmero de valores propios np = 0, 1, 2, . . .
p
~,
64
En el caso clsico no aparece la variable espn y las partculas son distinguibles
Q de forma
que cada configuracin o forma de ocupar los estados puede realizarse de N !/ np formas
distintas, que no es mas que el coeficiente multinomial de Newton, de donde:
N !N
X N ! Y p np 1 2 X
Zcl = Q e = e + e +... = exp [ p ] ,
P np ! p p
{np ; np =N }
resultado que ya obtuvimos antes salvo el factor 1/N ! (factor de indistinguibilidad cuntico)
65
Zcl factoriza pues no hay correlaciones entre elementos. Pero aqu, los p no son inde-
pendientes: se desarrollan correlaciones debidas a la indistinguibilidad.
102
Clculo en la macrocannica
(N ) |mi = Em (N )
donde Em , es el autovalor del hamiltoniano, H |mi ,
indica que cada estado del sistema con N elementos est caracte-
rizado por el conjunto de nmeros cunticos m. En nuestro caso
(N ) P P
m = {|np, i}p, y Em = p, p np con N = p, np luego
tenemos:
P P
X
X p np np
(T, V, ) = e p,
e p,
N =0 {np } !
P
np =N
p,
!
X
X XX
= exp ( p ) np
( )
N =0 P p
np ; np =N
p,
X
X
= exp [ (n1 + n2 )] exp [ (1 n1 + 2 n2)]
N =0 n1 ,n2
(n1 +n2 =N )
X XX X
X
N
= n1 n2
a b = aN n2 bn2
N n1 n2 N =0 n2 =0
(n1 +n2 =N )
103
fija N, suma para n2 desde 0 hasta N (n1 ya viene determinado
n1 = N n2 )
N =0 n2 = 0 1
N =1 n2 = 0, 1 a+b
N =2 n2 = 0, 1, 2 a2 + ab + b2
N =3 n2 = 0, 1, 2, 3 a3 + a2 b + ab2 + b3
.............. ..................... ................................
X
X
YX
n1 n2
= a b = n , = a, b, . . .
n1 =0 n2 =0 n=0
Luego
YX
= e(i)ni
i ni =0
104
que se puede extender a un conjunto de niveles etiquetados por
i = p, de donde se sigue:
YYX
= exp [ ( p ) np ] (25)
p np
Bosones: n = 0, 1, 2, . . . , luego
YY
(T, V, ) = 1 + e(p ) +2 e(p ) +
serie geomtrica
p
YY 1
(p )
= 1e
p
105
106
Termodinmica del Gas ideal cuntico: Ecuacin de esta-
do. Parmetro de degeneracin. Estadstica de los nmeros
de ocupacin
Una vez que ya tenemos la funcin de particin macrocannica
factorizada para el gas ideal cuntico vamos a obtener toda la ter-
modinmica.
El potencial macrocannico es
V P (T, ) = kT ln (T, V, )
XX n o
(p )
= kT ln 1 e
p
Pero:
P
p es independiente de , por lo que slo produce un factor
g, la degeneracin del nivel p , y
para sistemas macroscpicos, valores consecutivos
P de p~ estn R
suficientemente prximos para poder sustituir p por V h3 d~p
(salvo para bosones a bajas temperaturas, un caso que se trata
con detalle ms adelante).
luego:
Z
V p2 /2m
V P = gkT 3 d~p ln 1 exp
h kT
Z h
3/2 5/2 3
i
P (T, ) = gm (2kT ) h dx x ln 1 exp x +
0 kT
Integrando
por partes y usando la longitud de onda trmica, (T ) =
h/ 2mkT , se tiene
Z
2 2 x3/2
P (T, ) = kT g3 dx (26)
3 0 exp (x /kT )
107
Conviene reemplazar la dependencia en , por ejemplo por la den-
sidad de partculas,
Z
P 2 2 x3/2 exp (x /kT )
n= = g3 dx
T 3 0 [exp (x /kT ) ]2
que, integrada por partes es
Z
2 x1/2
n (T, ) = g3 dx (27)
0 exp (x /kT )
P = P (n, T )
U 3 2
u (T, ) = P P = u,
V 2 3
108
No es fcil obtener ecuaciones de estado explcitas a partir de las pa-
ramtricas (26) y (27); en lecciones posteriores veremos cmo hacer-
lo en casos particulares. Slo muy recientemente se han encontrado
algunas, por ejemplo para un gas de fermiones puede escribirse66
Z
3
P = nkT + 3/2 nP dy y 2 [n (y) P (y)]
2 g 1
donde
1 n (ln y)1/2
n (y) = tan Z
(ln )1/2
1 + n y d
1 ( y)
1P (ln y)3/2
P (y) = tan Z 3/2
(ln )
1 + P y d
1 ( y)
con
4 2g 8 2gkT
n = , P = .
n3 3P 3
que es bastante compleja y debe resolverse numricamente.
Para la mayora de los casos, como primer mtodo para obtener in-
formacin fsica relevante con poco esfuerzo, estudiamos la ecuacin
de estado para pequeas fugacidades,
z e/kT 1,
Z Z Z
x3/2
x3/2 1
dx = dx x 1 = dx x3/2 ex z 1 z ex
exp (x /kT ) e z
0 Z0 0
= dx x3/2 ex z 1 + ex z + e2x z 2 +
0
66
Hecho, por ejemplo, en M.M. Nieto, J.Math.Phys. 11, 1346 (1970), donde tambin se
trata la versin relativista de esta ecuacin de estado (ver luego para fermiones).
67
Para concluir esto es necesario conocer previamente la dependencia de en n y T ;
de todos modos, usando la expresin del potencial qumico para el gas de Boltzmann,
ya puede verse con facilidad que z 0 nos acerca al caso clsico pues para el gas de
Boltzmann z = n/f (T ) luego = kT ln z = kT ln n kT ln f (T ) < 0 pues n 1 para
preservar el caracter ideal y f(T) crece con T. Luego la fugacidad del gas de Boltzmann es
muy pequea.
109
que conduce a
3 5/2 5/2 2
P = gkT z 1 + 2 z+3 z +
3 3/2 3/2 2
n = g z 1 + 2 z+3 z +
Invirtiendo sta,68
1 3 3/2 1 3 1 1 2 6 2
z = g n 1 2 g n+ g n +
4 33/2
5/2 1 3 1 2 2 6 2
P = nkT 1 2 g n+ g n +
8 33/2
P = nkT 1 01768 g 1 3n 0 0033 g 2 6n2 + (28)
110
para fermiones ( = 1): la presin (para valores de n y T da-
dos) es mayor que la del gas perfecto: es como si existiera
una fuerza repulsiva adicional entre fermiones (consecuencia
del principio de exclusin en un entorno cooperativo) que des-
favoreciese configuraciones con partculas muy prximas entre
s;
para bosones ( = 1): la cooperatividad se manifiesta en un efecto
como el que producira una fuerza atractiva.
111
Nmeros de ocupacin
En formalismo de segunda cuantizacin, el valor medio, hnp i ,
del nmero de ocupacin del estado (~p, ) de momento lineal p~
y 3a componente de espn es
hnp i = Tr a+ a
p p
con el operador densidad macrocannico
1 XX
= exp ( p ) a+
p ap .
p
~
YX
= e(p )n
p~, n
y por definicin
P ( )
n exp [ ( p ) n] X
hnp i = Pn = ln exp [ ( p ) n]
n exp [ ( p ) n] () n
YX Yh i
(p )
exp [ ( p ) n] = 1e ,
~p, n ~p,
luego
h i e(p )
(p )
hnp i = ln 1 e = .
() 1 e(p )
112
Esto es,
1
hnp i =
e(p )
que me da el nmero de ocupacin promedio en el estado (~p, ),
con 69
+1 estadsitica o distribucin de Bose-Einstein
= 1 Fermi-Dirac
0 Maxwell-Boltzmann.
e(p )
1 + hnp i =
e(p )
que, para fermiones (es decir, n = 0, 1 de modo que 0 < hnp i < 1),
representa el nmero medio de huecos en el estado (~p, ) (densidad
de huecos). En el caso de bosones no tiene esta interpretacin pues
hni n no est acotado.
69
Puesto que hnp i es aqu independiente de , estamos suponiendo que hay el mismo
nmero medio de partculas en cada nivel con cada una de las posibles orientaciones de
espn, de modo que el estado de equilibrio del gas es no polarizado, con valor medio de
espn total nulo.
113
114
Tema 10: Estadstica de Fermi-Dirac
Gas ideal de fermiones degenerado; energa de Fermi,
potencial de Gibbs.
Gas relativista muy degenerado.
Enanas blancas: modelo de Chandrasekhar.
Magnetismo y gas de Fermi (paramagnetismo de Pau-
li, diamagnetismo de Landau, efecto Hass-van Alp-
hen).
Propiedades del gas electrnico en metales; coeficien-
tes de transporte.
3/2
1 h2
= ng 1.
2mkT
115
Energa de Fermi, F : la del nivel ms alto que se alcanza de
este modo suponiendo que el gas es ideal. Todo los niveles por
debajo de ste constituye el mar de Fermi.
1
hnp i = = (F p ) para (T 0o K) ,
e(p ) +1
(29)
hnp i
1
T = 0o K
con la funcin paso de Heaviside:
(gas ideal)
F
p
de modo que
< 0 si p < F
p es
> 0 si p > F ,
(T = 0o K) = F
De la ecuacin de estado para un gas ideal cuntico de fermiones:70
70
Se sigue haciendo x = en el resultado obtenido en una leccin anterior:
Z
2 x1/2
n (T, ) = g3 dx
0 exp (x /kT )
116
R 1/2 R
n= 0 d e() +1
= 0 d 1/2 (F p )
(30)
R F 3/2
= 0 d 1/2 = 32 F ,
2/3
6 2 ~2 2/3
F = n (31)
g 2m
r 1/3
F 6 2
pF = ~ n1/3
2m g
71
Es decir un gas de fermiones en el 0 absoluto es bastante energtico; por el PP, el
momento de cada fermin ser grande, en general, lo que hace que el gas de fermiones com-
pltamente degenerado sea un buen modelo, al menos cualitativo, para algunas situaciones
fsicas, como vemos a continuacin, es decir es un modelo con mucha fsica.
117
El nmero de ocupacin hnp i es entonces:72
exp [ ( )]
hnp i (1 hnp i) =
{exp [ ( )] + 1}2
kT
0
0 F p 0 0 F p 0 0 F p
118
Z
1
I= d f ()
0 exp [ ( )] + 1
Z
2 2 7 4
I= d f () + (kT ) f () + (kT )4 f () + , (32)
0 6 360
Demostracin:
Z Z
exp [ ( )]
I= dG () ; G () = d f () .
0 {exp [ ( )] + 1}2 0
Haciendo = x, y luego x = y :
Z Z
ex y ey
I= dx G (x + ) 2 = dy G + .
(ex +1) (ey +1)2
119
Z
X 1 yn
(n) ey
I= dy G ()
n0
n! n
(ey +1)2
X 1 1 Z
(n) n ey
= G () dy y
n0
n! n
(ey +1)2
X Z
1 1 (2n) 2n ey
I= G () 2 dy y .
n0
(2n)! 2n 0 (ey +1)2
R y
Para n = 0, es 0 dy (ey e+1)2 = 1
2
y, para n > 0, se tiene
integrando por partes:
Z Z
ey
dy y 2n
= 2n dy y 2n1 (e y + 1)1
0 (ey +1)2 0
= 2n 1 212n (2n 1)! (2n)
= 22n1 1 2n |B2n|
donde Bn son los nmeros de Bernouilli.75
Finalmente,
Z X
2
I d f () + 22n1 1 2n |B2n| (kT )2n f (2n1) ()
0 n1
(2n)!
2
73
En efecto, cambiando y y en ey / (ey +1) :
ey ey e2y ey ey
= = =
(ey +1)2 (ey +1)2 e2y [(ey +1) ey ]2 (1 + ey )2
74
La integral de un funcin impar en un intervalo intervalo simtrico simtrico es nula.
n1
75
Aqui hemos utilizado la conocida relacin B2n = (1) (2n)!
22n1 2n (2n) n 1.
120
cuyos primeros trminos. producen (32) QED.
" 2 #2/3
2 kT
F = 1 + +
8
" 2 #
2 kT
= F 1 + (33)
12 F
76
Haciendo f () = 1/2 :
Z
2 1
I1/2 = d 1/2 + (kT )2 1/2 +
0 6 2
" #
3/2
2 z }| { 2
2 2 kT 2 3/2 2 kT
= 3/2 + 2 1/2
+ = 1+ +
3 12 3 8
2/3
z }| { 2/3 " 2 #2/3
2 2/3
6 ~2 25/3 2/3 mg 2/3 2 2 kT
F = 1+ +
g 2m h2 3 8
121
que es un desarrollo vlido para bajas temperaturas.
1 u 2 kT
CV = = k,
n T n 2 F
77
u = U/V es la densidad de energa. En efecto
T S P V + N
= U/V =
u = T s P + n
V
pero P = B(T )I3/2 (ecuacin 26) luego derivando respecto a T sale s = S/V =
1 P V
= PT + 23 PT Tn donde para el tercer trmino hemos usado la expresin de
V
T ,V
n = P = 32 P. Se tiene P = 13 u
en funcin de y T. De ahi se sigue el resultado u si el
T
gas ideal es relativista.
78
En efecto, se tiene:
" 2 #
3 5 2 kT
u = nF 1 + +
5 12 F
79
De hecho vimos que para fonones (que son bosones) a T bajas el CV iba como T 3 .
122
Aplicacin del gas de fermiones muy degenerado
Gas relativista muy degenerado.
Enanas blancas: modelo de Chandrasekhar.
Introduccin
R.J. Tayler, Stellar Evolution, Rep. on Progr. in Physics 31, 167 (1968)
J. Nuckolls et al, Laser-induced thermonuclear fusion, Physics Today, August 1973, p. 46
Huang p. 230; Pathria p. 243
W. Herbst & G.E. Assousa, Supernovas and Star Formation, Sci. Am. 241, 138 (1979)
H.A. Bethe, Supernova Mechanisms, Rev. of Modern Physics 62, 801 (1990)
123
la escala temporal caractercitica para las variaciones en las
propiedades del Sol es de miles de millones de aos,
observaciones como las relacionadas con novas, supernovas,
etc. podran ser accidentes ms que etapas de una evolucin
estndar, etc.
pero hay observaciones, como el esquema de Hertzsprung-Rusell(HR),
a las que se les atribuye un gran valor indicativo del proceso seguido
por una estrella tpica:
124
10 % de excepciones:
prescindiendo de otros
detalles, estrellas gran
radio + brillantes de lo
que corresponde a su color
rojo, y otras pequeo
radio, gran densidad,
mates para su color blanco
125
Nube de hidrgeno homognea (las molculas se mueven al
azar). En un cierto momento se forma un ncleo ms denso
que vence el movimiento aleatorio de las molculas y empieza
un atraccin gravitatoria. En este momento la nube inicial
(transparente) se hace opaca y aparece una tendencia hacia el
colapso gravitatorio, al principio rpido, y despus ms lento
debido a la aparicin de gradientes de presin que se oponen
al colapso.
Aumenta la temperatura dendro del ncleo por bombardeo
de partculas que son aceleradas gravitatoriamente hacia el
ncleo. A la temperatura de 105 se produce la ionizacin del
H. El e exterior salta. El ncleo es entonces un gas mezcla de
protones y electrones.
Sigue la concentracin y aumentando la temperatura hasta
los 107 momento en el que la colisin entre protones es tan
violenta que se vence la repulsin culombiana y y aparecen las
interacciones fuertes originando la fusin nuclear H He,
que es un proceso muy exotrmico donde se emiten muchos
fotones y neutrinos (que no colapsan gravitatoriamente), la
estrella se hace muy brillante (visible). Es una fase estable con
una gran emisin de energa y la estrella pasa la mayor parte
de su vida.
Como consecuencia de la presin de radiacin, las capas ms
externas tienden a expandirse y aparecen como ms rojizas
(ms fras) pero muy luminosas dando lugar a una gijante ro-
ja, pues consecuencia de la expansin alcanzan un gran tama-
o. (Betelgeuse). La evolucin de este proceso depende mucho
del tamao de la estrella. Para tamaos grandes el proceso se
acelera, la estrella consume su combustible y muere antes
Cuando se ha consumido todo el H se frena la expansin y
vuelve la contraccin y no hay nada que la detenga (no hay
combustible nuclear) slo He. Se tiene entonces una enana
blanca (poco brillante)82
82
Suponiendo que la masa y la temperatura son tales que no se pueden iniciar procesos
de fusin de tomos de He
126
En este contexto se admite el siguiente modelo idealizado de
estrella enana blanca:
masa:
N
M Nm + 4mp = N (m + 2mp ) 2N mp ,
2
densidad de electrones:
N M/2mp 107gr/cm3
30 3
n= = = O 10 electrones/cm ,
V M/ 2mp 3 4 1024gr
momento de Fermi:
1/3
pF = 3 2 n ~ 1010 1027 O 1017gr cm/s ,
energa de Fermi:
F = p2F /2m O 107ergios = O 105eV ,
83
kT = (8,617385 105 eV K1 ) (107 K) = 862 eV
84
mn = 1,67493 1027 kg mp = 1,67262 1027 kg; me = 9,10939 1031 kg mp
127
y temperatura de Fermi:
TF F /k O 1010 K
128
Gas relativista de Fermi en el estado fundamental (es de-
cir, completamente degenerado)
y la velocidad:
2pc2 p/mc
u = p = cq = c tanh
p 2 m2 c4 + p2c2 2
1 + (p/mc)
X Z
2V pF
E0 = 2 3 4p2 d p
h 0
|~
p|<pF
Z
8V m4 c5 F
= 3
senh2 (cosh 1) cosh d
h 0
mu
85
el momento de una partcula relativista con velocidad u es p = mu y la
1(u/c)2
2
p
energa cintica = mc2 mc2 . De la definicin de se tiene (u)2 = c2 ( 2 1) = m 2 . De
2 2 4 2 2 2 2 4
esta ltima expresin multimplicando por c se tiene m c = p c + m c . tomando raiz
cuadrada y sustituyendo en la expresin de da el resultado requerido.
86
En espacio de momentos, para una partcula en una caja, la ec. de Dirac tiene la misma
solucin que la de Scrdinger, luego la densidad de estados para un gas relativista es la
misma que para el caso no-relativista.
129
y la presin (presin de degeneracin):87
Z pF Z pF
1N 2 8
P0 hp ui0 = 3 pu 4p2 d p = 3 u p3 d p
3V 3h 0 3h 0
4 5 Z F
8m c
= 3
senh4 d
3h 0
pF 1/3 ~
x senh F = = 3 2 n , (34)
mc mc
p
A (x) = x 2x 3 x2 + 1 + 3 senh1 x
2
88
p
2 3
B (x) = 8x x + 1 1 A (x)
se tiene
m4 c5 m4 c5
P0 = A (x) , E0 = V B (x) (35)
3h3 3h3
87
Puede obtenerse P0 = E0 /V notando que F (tq pF mc senh F ) depende de V ;
de hecho, pF puede definirse: 2 hV3 43 p3F = N. Alternativamente, recurdese que
1 1M 1 Nm 1N
P hu ui = hu ui = hu ui = hp ui
3 3V 3 V 3V
130
El comportamiento asinttico de las funciones A (x) y B (x) es:
8 4 1 5 no relativ.
x5 x7 + x9 x11 + , x1
A (x) =
5 7 3 22 mc pF
2x4 2x2 + 3 ln 2x 12
7
+ 4x5 2 + , x 1
relativista
mc pF
( 12 3 1 15 11
x5 x7 + x9 x + , x 1, no relativ.
B (x) = 5 7 6 176
6x4 8x3 + 6x2 3 ln 2x 14 + 4x3 2 + , x 1, relativista
es decir, una caja o un agente externo tienen que hacer el negativo de este trabajo para
comprimir la estrella desde un estado de densidad cero hasta un estado de densidad finita
(suponemos una esfera de radio R).
Es claro que este agente externo es la interaccin gravitatoria en la propia estrella, que
origina una energa de la forma GM 2 /R.
Notar que escribimos GM 2 /R para conseguir las dimensiones correctas: el sistema estrella
slo involucra una masa y una longitud. La teora no permite determinar que ha de
ser O (1) (no hay argumentos para esperar otra cosa); Galindo-Pascual escribe =
3/5. (Estrictamente hablando, (36) define ; su contenido fsico queda determinado por la
hiptesis de que es del orden de la unidad.)
131
Z R
2 G M2
P0 4r d r = (36)
R
1/3
pF N ~
x= = 3 2 (38)
mc V mc
132
m4 c5 8 5 (34) 5/3 (38) M 5/3
P0 x n .
3h3 5 R5
3 (9)2/3 ~2
R 5/3
M 1/3
40 G m mp M 1/3
1/2
9 3 (~c/G)3/2
M0 = .
64 3 m2p
Conclusin:
el tamao de una EEB, definida como se indic, disminuye al
aumentar su masa (las pequeas son ms densas),93 tanto si
el gas es relativista como si no;
no es posible tener una EEB con M > M0 (la ec. de equilibrio
no tiene soluciones reales).
133
no relativista
relativista
M0
M
Conclusiones
Dependencia de R de la EEB con M Las ms pequeas son
las ms densas.
El hecho de que el gas sea relativista o no, no aparta nada
importante cualitativamente.
134
Hay efectos cunticos importantes alta degeneracin.
No es posible tener una EEB con M > M0 pues no hay solucin
para R. Si M > M0 la presin de degeneracin no es suficiente
para evitar el colapso gravitacional por esto no tiene soluciones
reales
Hiptesis simplificatorias
Nuestro modelo es un sistema de fermiones relativista muy de-
generado. Hemos supuesto que los ncleos de He se encuentran
en reposo y que por tanto su movimiento no va a influir en la
dinmica del sistema. Sim embargo las temperaturas son tan
altas que la energa cintica de estos ncleos ser muy alta
luego modelos ms realistan deben incluir el efecto del movi-
miento de los ncleos.
135
Figura 1: Masa versus radio. La curva discont. se sigue de las consideraciones
de degeracin no-relat. de Fowler: el radio para la configuracin de equilibrio
disminuye, pero slo se hace cero para masa infinita. En este caso, una estrella
de cualquier masa puede confortablemente equilibrarse en el estado EEB. Esto
no ocurre para configuraciones completamente degeneradas. (Dibujo del propio
Chandrasekhar.)
136
mxime cuando estamos viendo que la densidad es muy gran-
de.
137
dad muy grande. Estos procesos son bastantes corrientes. En
el siglo pasado se han fotografiado unas 50 supernovas.
138
139
140
Otras aplicaciones del gas de electrones
Magnetismo y gas de Fermi: paramagnetismo de Pau-
li, diamagnetismo de Landau, efecto Hass-van Alphen.
Propiedades del gas electrnico en metales: calores es-
pecficos; efectos magnticos; coeficientes de transpor-
te, nmero de Lorentz.
Introduccin
Mz
lm ,
H0 H
~ = Mx , My , Mx M ~ es la magnetizacin con
negativa; M H
que responde el sistema al campo aplicado;
~ originando paramag-
sus espines tienden a alinearse con H,
netismo, con > 0.
94
El ncleo influye poco, salvo en formacin de la electrnica, pues resulta pesado para
tener mm orbital apreciable y el intrnseco es unas 103 veces el del electrn.
141
(Adems, si H muy intenso y baja T, manifiesta comportamiento
oscilatorio.)
H
Mz = n
BJ (x) , x
kT
142
Mz = n
= ngB J = = 0
1 e ~ ~
p~ + A ~ H
2m c
p2 ~
= ~ H
2m
~ luego ~H
En estas condiciones, ~ paralelo o antiparalelo a H, ~ =
~
95
El momento medio de los electrones es siempre muy alto, propiedad del gas de Fermi
incluso en el 0 K, consecuencia del PP (tenemos como mucho que la mitad de los espines
apuntan en la direccin del campo), por lo que no puede alcanzarse saturacin magntica: si
dos electrones estn en el mismo estado, es imposible que tengan sus espines en la direccin
~ Por lo tanto, ha de esperarse una susceptibilidad lmite, 0 , que ser independiente de T
H.
aunque probablemente con fuerte dependencia en la densidad del gas (F ).
143
~ (donde los valores propios de ~ H
H ~ son H con = 1 ) y pueden
2
distinguirse dos tipos de partculas en el gas con energa respectivas
p2 paralelos
=
H
2m +antiparalelos
e. total F
p2
e. cintica 0 (F
H)
2m
momento H)]1/2
0 p [2m (F
+
V (2m)3/2 h
4 3/2 3/2
i
M =
N N = (F +
H) (F
H)
3h3
96
En el tomo de T-F , escribamos que la mx energa de un electrn en el nivel de Fermi
es max = F e, mientras que aqu escribimos max = F = p2 2m H. Esto es, un
factor externo (H) cambia el nivel de Fermi, mientras que el interno debido al ncleo no
lo cambia. (De hecho, F la tenemos definida para un sistema sin interacciones mutuas, y
las contiene.)
97
Recordad que, puesto que los valores propios del momento son px = (h/L) nx , con
nx = 0, 1, 2, . . . , a cada celda de una red cbica le corresponden L3 /h3 = V /h3 estados.
P R R p 2
Luego N = p 1 N = hV3 d~ p = 4V
h3 0
p dp, que da el resultado requerido.
144
que implica la susceptibilidad por unidad de volumen:98
Mz
0 = lm
H0 H V
i
V (2m)3/2
4 3 h 1/2 1/2
= lm
(F +
H) + (F
H)
3h3 H0 2
V (2m)3/2 1/2 3 2 n
4
= F =
h3 2 F
que no depende de la temperatura, slo de la densidad y de F
y por lo tanto el PP da un lmite inferior a la suceptibilidad
como en experimentos reales.
a temperaturas altas, sin embargo, sistema se comporta
como gas de Boltzmann (estadstica de Maxwell-Boltzmann),
y se sigue lo obtanido ya en contexto
de paramagnetismo de
2
Langevin, es decir, Mz = N /3kT (1 + 1/J) H luego, para
1 2
J = 2 , es Mz = N HkT, de donde:
2
n C
= = , 2k.
C = n
kT T
Termodinmica
Sea np = # partculas con momento p ~ cuyo momento magntico
es paralelo (antiparalelo) al campo. La energa total del gas es
98
En la ltima igualdad aqu hemos usado que
2 2/3 2 )
F = 6g ~
n 2/3 3/2 6 2 ~3
2m = F = n =
g=2 g (2m)3/2
1/2 n 2 ~3 3 n h3
F = 3 3/2
=
F (2m) 2 F 4 (2m)3/2
145
X p2 2
p
E= H n+
p + H n
+ p
2m 2m
p~
X p2 X
= n+ + n
H n +
n
2m p p p p
p~ p~
X p2
= n+
H N + N
p + np
2m
p~
X N
X
X
H(2N + N ) P p2 + P p2
Z= e e ~ 2m np e ~ 2m np
p p
N + =0 n+p
np
P P
( p =N )
n+ +
( p =N )
n
N h
HN
X
HN +
2 +
+
i
=e e Z0 N Z0 N N
N + =0
donde
X P p2
~ 2m np
Z0 (N ) = e p exp [A0 (N )]
np =0,1
P
( np =N )
146
partculas sin espn de masa m, y A0 (N ) es la energa libre de
este sistema ficticio.
Tomando logaritmos:
1 1 X
N
ln Z (N ) =
H+ ln HN + A0 N + A0 N N +
exp 2
N N +
| {z }
N =0
f (N + )
0 N + 0 N N + = 2
H
99
El logaritmo de esta suma es igual al logaritmo del trmino mayor salvo contribuciones
de O (ln N ) , luego despreciables para sistema macroscpico (N 1022 ) y en el lmite-
Termodinmico.
PN De hecho, puede interpretarse el sistema en la macrocannica debido a la
N + =0 . En estas condiciones, hemos mostrado que la aproximacin es rigurosamente vlida
para sistemas macroscpicos.
100
Puede comprobarse que esta ecuacin slo tiene una solucin real entre 0 y N, que
corresponde a un mximo.
101
En efecto:
A0 (N N + ) A0 (N N + ) (N N + ) A0 (N N + )
= =
N + (N N + ) N + (N N + )
147
con 0 (N ) el potencial qumico del sistema ficticio.102
N+ N
M = N r, (con el parmetro 0 r 1)
es decir,
1+r 1r
N+ = N, N = N,
2 2
y se tiene
1+r 1r
0 N 0 N = 2
H
2 2
102
Se est imaginando el sistema real como una mezcla (termodinmica) de componentes
+ y con intercambio continuo de partculas. Se puede llegar a la misma conclusin por
tanto igualando los potenciales qumico en el equilibrio: + = . Pero
p2 +H (antiparalelo) (N partculas)
=
2m H (paralelo) (N + partculas)
sinespn
partculas
+ = 0 N +
H (paralelas)
= 0 N + H (antiparalelas)
y se sigue lo indicado
148
de donde103
1 M 22 n
= lm =h i
V H0 H 0 (xN )
x x= 21
Para hacer explcito este resultado, hay que usar aqu la expresin
adecuada para el potencial qumico del sistema de partculas sin
espn; en particular, se tiene
luego
2 N H
2
N r h
M = i para H 0
0 (xN )
x x= 12
y se sigue lo indicado.
104
En efecto, es
2/3
3N h2
F =
4gV 2m
luego se tiene para el sistema ficticio (g = 1, y cambia N xN ):
2/3
2/3 3N h2
0 (xN ) = x
4V 2m
As: 2/3
0 (xN ) 2 h 1/3 i 3N h2 22
= x = real
x 3 1/2 4V 2m 3 F
x= 12 | {z }
24/3 /3
y se sigue:
2 n3
2 3 2 n
0 = real
=
2
2 F 2 real
F
149
es un resultado que ya habamos obtenido antes mediante otras
consideraciones.
a T s bajas, pero finitas, es
" 2 #
2 kT
F 1 +...
12 F
y se sigue que
" 2 #
2
kT
0 1 +...
12 F
g z2 z3
= 3 z 3/2 + 3/2 + ,
2 3
si T :
3 N N 3
z=e = kT ln .
g V gV
150
0 (xN ) 1 N 3
= kT = 2kT
x x= 12 xN 3 gV V x= 12
2
n
=
kT
n3
= 1 5/2 + ,
2
Diamagnetismo de Landau
El diamagnetismo es consecuencia de la libertad de movimiento que
tienen los electrones en presencia de un campo magntico. Pero (T.
de Bohr-van Leeuwen) la teora clsica es incapaz de explicar efec-
tos magnticos, luego ha de estar asociado a fenmenos cunticos.
Landau (1930): asociado con cuantizacin rbitas de partculas car-
gadas en presencia de H,~ es decir, con el hecho de que la Ecintica de
las partculas asociada con el movimiento a H ~ no puede tomar
cualquier valor. Esta circunstancia = adicional a la paramag-
ntica, cualitativamente semejante a sta, pero negativa.
~ sin
Sea un gas de electrones libres, salvo que interaccionan con H,
espn (ya conocemos las consecuencias de este grado de libertad),
confinado en volumen V. Con notacin usual, el hamiltoniano es
1 e ~ 2
H= p~ A + (~r) .
2m c
151
Resolviendo ec. de Schrdinger(ver [1] en la bibliografia), sus valores
propios son105
e~H 1 p2z
(j, pz ) = j+ + , J = 0, 1, 2, . . . (39)
mc 2 2m
con degeneracin
eH
g = Lx Ly
hc
Puede estimarse:
105
La interpretacin de (39) es que el movimiento del electrn puede descomponerse en:
plano (x, y) : describe crculos bajo accin de la fuerza de Lorentz que origina velocidad
angular constante = eHmc. La energa asocida est cuantizada en unidades ~ =
e~Hmc.
152
Z Z Z Z
Lx Ly
g= d nx d ny = 2
d px d py
h
p= L n h
en un crculo e~H p2 2
x +py e~H
mc j 2m mc (j+1)
Lx Ly e~H eH
= 2m {(j + 1) j} = L x L y ,
h2 mc hc
X X
X h i
(j,pz )
ln = ln 1 + z e = g ln 1 + z e
j=0 pz
Z
eH X h i
Lz
= d pz Lx Ly ln 1 + z e(j,pz )
h hc j=0
Z X
h i
V eH (j,pz )
= 2 d pz ln 1 + z e (40)
hc j=0
Implicaciones:
106
donde la densidad de estados en el espacio de momentos en el eje z es Lz /h.
153
T s altas: N es finito luego, a T s altas, cuando tienden a
ocuparse los estados con no cuntico j alto, es necesario z 1
para asegurar convergencia.107 Entonces108
Z X
V eH
ln 2 dp z exp [ (j, pz )]
hc j=0
Z X
zV eH p2 e~H 1
= dp exp j+
h2 c j=0
2m mc 2
Z
zV eH
p2 X e~H 1
= d p exp exp j+
h2 c 2m mc 2
| {z }|j=0 {z }
integral gausiana
serie geomtrica
integral serie
}| {z
z
1/2 }| 1{
zV eH 2m e~H zV x
= 2 senh = ,
h2 c 2mc 3 senh x
con x B H (B = magnetn de Bohr)
y tambin que (z 1)
zV x
N ln ,
3 senh x
y109
1
M ln = B N L (x)
H
107
De hecho, ya sabemos que z 1 corresponde al comportamiento gas de Boltzmann
caracterstico de T s altas.
108
Si z e con > 0 (e < 1) y z 1, se tiene 1, luego podemos usar aqu
que ln (1 + ) .
109
En efecto, se tiene sucesivamente:
154
1
donde L (x) coth x es la funcin de Langevin que ya
x
apareci antes. De hecho, este resultado coincide con el para-
magntico, salvo que ahora es M < 0 (independientemente
del signo de la carga eltrica de la partcula!), como se observa
en situaciones diamagnticas. Es la tendencia a que se presen-
ten momentos magnticos opuestos a la direccin del campo
x 2 H
M B N = N B
3 3kT
y
n2B
=
3kT
= Pauli
+ Landau
N
z }| {
zV 1 x cosh x zV x 1
= B = B 3 coth x
3 senh x senh2 x senh x x
110 1 x x3 2x5 x x3
coth x = x + 3 45 + 945 = L (x) = 3 45 +
155
M ~2 + O ~4 ,
F
z exp = exp + 1.
kT kT
1
0 = n2B F ,
2
kT B H F
2
111 2 kT
F 1 12 F
112
Se usa la frmula de Euler,
X Z
1 1
f j+ f (x) d x + f (0)
j=0
2 0 24
156
la funcin de particin macrocannica tiene una parte oscila-
toria de la forma:114
X
(1) cos B H F 4
(ln )osc (B H) 3/2
3/2 senh
=1 B H kT
114
Algunos detalles en Pathria p. 231.
115
En consecuencia, F puede determinarse experimentalmente, de manera relativamente
sencilla, por este mtodo.
157
158
Propiedades del gas electrnico en metales
Calores especficos
159
el valor experimental es aproximadamente 3R por mol para
slidos monoatmicos, en lugar de la prediccin 92 R117
si la T no es tan alta, CV depende fuertemente en T.
Sommerfeld (1928) not que era muy importante tratar el gas elec-
trnico como un gas de Fermi degenerado. En efecto, tomando
g = 2, tenemos
2/3
1 3N h2
electr
F = ,
22/3 4V 2m
luego la prediccin es CV = 92 N k = 29 R
118
Aqu, y en otras muchas aplicaciones del modelo gas de Fermi, se est despreciando
la posible energa de interaccin, en particular, la debida a la presencia de un potencial
producido por los iones (ej. los localizados en los nudos de la red), que produce efectos
importantes como comentamos luego, y las interacciones entre electrones. Respecto de
stas, hay que notar que la interaccin entre dos electrones aislados es coulombiana, nada
160
F 5146 1012 erg = 3 21 eV
TF F /k 33 104 K
para el sodio.120
161
2 kT
CV = k . (42)
2 F
Pero, a T s orinarias, es (kT F ) T TF O 102 , es decir,
esta contribucin es del orden de 102 R por mol, despreciable frente
a la de la red, 3R por mol. En consecuencia, los metales a T s
ordinarias han de presentar la contribucin bosnica (41), lo que
confirman los experimentos.
CV = T + T 3
122
Se confirma comportamiento lineal de CV T en funcin de T 2 en el rango 1-4 K. Por
ejemplo, Corak et al, Phys. Rev. 98, 1699 (1955); ver tambin Pathria p. 234.
162
Comparando con (41) y (42): la pendiente de estas rectas (a muy
baja T ) permite determinar D (temp de Debye), mientras que el
corte con el eje vertical est relacionado con F .123
Efectos magnticos
Pauli (1930): tiene que ver con el PP, pues dos electrones en el
mismo estado no pueden alinear sus espines en direccin H, ~ lo
que impide saturacin a bajas T s (es decir, T TF ), cuando
hay mucha degeneracin. En consecuencia, la de un metal ha de
ser suma de la contribucin de Pauli (debida al espn de electrones
123
Para el cobre se obtiene as D = (3438 0 5) K, que se compara muy bien con D =
345 K obtenido a partir de las constantes elsticas segn procedimiento descrito en otra
parte y = 7 104 J / mol grado2 . A partir de F , suponiendo que el gas es ideal, se sigue
= 5 4 104 J / mol grado2. La diferencia desaparece suponiendo una masa efectiva para
los electrones del cobre 13 veces la electrnica. (Insisto, sin embargo, que este concepto de
masa efectiva parece tener una validez bastante restringida.)
124
El mismo resultado se obtiene como lmite T de la susceptibilidad de un gas de
Fermi, como se vi.
163
degenerados) y de la de Landau (originada por el momento orbital
~
de los electrones en presencia de H):
3n2 n2B
Pauli = , Landau =
2F 2F
Coeficientes de transporte
164
slidos.125 Una propiedad de y es la ley (fenomenolgica) de
Wiedemann-Franz (1853) que establece que L = /T es indepen-
diente de T y del metal en cuestin.126
2
k
L=3 = 223 108 W K2 ,
e
0C 100 C
Pb 247 256
Pt 251 260
Ag 231 237
Au 235 240
165
Si el gas est sometido a campo elctrico E, cada partcula sufre
aceleracin a = eE/m, viajando una distancia 12 at2 en cualquier
tiempo t, o
1
2
una distancia media a t en el tiempo caracterstico hti
2
1
2
a t es la distancia media que recorren
2
los electrones antes de sufrir otra colisin
y
1
hvi = a t2 hti es la velocidad media de arrastre,
2
166
luego
2
1
2a t eE eEF
hvi = = aF = F =
hti m mvF
Puesto que
eF
movilidad velocidad media de arrastre por u. de campo = mvF
se tiene finalmente:
ne2 F
= (43)
mvF
167
Suponiendo que la direccin de la velocidad es aleatoria, puesto
que hay 6 direcciones, la superior P ser atravesada en promedio
por 16 nvF electrones por u. de tiempo viniendo de A y otros tantos
viniendo de B.
luego
1 1 2 k2T
= nvF F = nvF F ,
3 T 3 2 F
1
6 nvF F 2 k 2 T n 2 k 2 F T
= 1 2
= .
2 mvF
3mvF
2
2 k
L= = = 2 44 108 W K2
T 3 e
168
Bibliografa
1 e ~ 2
H= p~ A + (~r)
2m c
Ax = 0, Ay = Hx, Az = 0.
Se sigue la ec de Schrodinger:
2
2 ieHx 2 2mE
+ + + 2 = 0. (44)
x2 y ~c z 2 ~
y se sigue para u :
" 2 #
2mE1 eH
u + 2
ky x u = 0, (45)
~ ~c
128
que origina un campo uniforme de intensidad H en la direccin z
169
2
donde E1 E 2m ~
kz2 es la enrrga asociada con el mov en
el plano transversal. Pero esta ecuacin es formalmente similar
a la ecuacin de onda de un oscilador armnico de frecuencia
eHmc y centro en x0 = ~cky eH, de lo que se sigue (la
frecuencia es 2L = eHmc y se tiene E1 = (2j + 1) ~L)
que los vp son
1 eH
E1 = j + ~
2 mc
y
1 e~H ~2kz2
E= j+ +
2 mc 2m
QED.
170
No existen soluciones si x0 queda fuera del volumen. Para un
pequeo rango de x0 cerca de las paredes, a uno y otro lado, la
presencia de la pared modifica la funcin de onda y E1. Para
un campo de 103 Gauss y un electrn de 2 eV, el radio es de
unos 5103 cm, luego esta regin puede generalmente despre-
ciarse.]
El intervalo permitido de ky es eH
~c L1 . Se siguen entonces
L1L2eH2~c estados para cada energa, QED.
171
172
Tema 11: BOSONES
Gas ideal de bosones degenerado
Condensacin de Bose-Einstein
Aplicacin del gas de bosones: superfluidez del helio lquido,
etc..
Introduccin
X
p
V P (T, ) = kT ln (T, V, ) = kT ln 1 exp .
p,
kT
(46)
donde se tiene = 1 para bosones.
ln X
N = kT = hnp, i
p,
1
hnp, i =
e(p ) 1
173
que expresa una cualidad esencial de la distribucin de Bose-Einstein.129
P (T, z) = kT 3Y (z)
n (T, z) = 3 y (z)
donde z = e/kT , = h ,
2xmkT
Z
4 3/2
Y (z) = d 1 ,
3 0 z e 1
y Z
2 1/2
y (z) = d 1
0 z e 1
con
p2
=
kT 2mkT
129
Puesto que este hnp, i es independiente de , estamos suponiendo que el estado de
equilibrio es no polarizado, es decir, con el mismo # partculas en promedio con igual
orientacin de espn en cada nivel, de modo que el espn total en el gas es cero.
130
Y (Z) y y(z) se definen en trminos de las funciones de Bose-Einstein definidas como
Z
1 x1
g (z) = dx
() 0 z ex 1
1
174
Gas ideal de bosones degenerado
y (z) = zY (z) .
X
X
5/2
Y (z) = m m
z , y (z) = m3/2 z m
m=1 m=1
y, sustituyendo aqu n = 3 y y P = kT 3 Y,
Y z Y
3 y = kT 3 = y= z
z kT z
132
La fugacidad tiene un distinto significado para fermiones: los trminos de esas series
tienen entonces signo alternado de modo que, aunque el radio de convergencia es z = 1
tambin, este punto no es una singularidad esencial, pues el denominador e +z no se anula,
y las funciones para fermiones pueden continuarse analticamente ms all de z = 1 y existen
para 0 z .
175
0z1
P 5/2 5
Y (1) = m=1 m 2
1342
P 3/2 3
y (1) = m=1 m 2
2612
y (z = 1) = .
puesto que
1 1 1 1 1
= 1/2 + 1/2 + 1/2 + 1/2 +
2 1 2 3 4
1 1 1 1
> + + + +
1 2 3 4
es mayor que la armnica, puesto que 0 < x < x. y la armnica diverge (a pesar de que
an 0)
176
tiene la forma de la figura:134
177
tal que la ecuacin de estado no tiene soluciones con significado
fsico para T < Tc (n)
Para reiterar este comportamiento, partimos de n3 1 , man-
tenemos T =cte, y aumentamos la densidad n
Llegamos ahora nc 3 y (1) que permite definir una densidad
o volumen especfico crticos
3 h3
n1
c = vc = = 3/2
T 3/2,
y (1) 2 612 (2mk)
Es decir,
Notamos entonces en
Z
2 1/2 p2
n3 = y (z) = d , =
0 z 1 e 1 2mkT
178
puede ser), es decir, aquella extraa desaparicin de partculas
refleja que el formalismo no da cuenta de ellas.
Z
V p2/2m
V P = kT 3 d ~p g ln 1 exp .
h 0 kT
Pero este paso slo est justificado cuando son finitos todos los
trminos en la suma sobre ~p de (46), lo que no ocurre para bosones;
es decir, despreciamos el trmino p = 0 que contribuye a P V de
acuerdo con el no medio de ocupacin del estado fundamental, que
resulta ser hn0 i N, de modo que, para un sist macroscpico (N
), produce una singularidad tipo delta que ha sido despreciada.
X
p2/2m
V P = kT ln 1 exp
p
kT
X
= kT ln (1 z) + ln 1 z e/kT
p6=0
179
P R
Aqu ya puede reemplazarse p6=0 por la integral 0 d p~; el
resultado es el de antes135 pero con las correcciones ln (1 z)
y z (1 z) , es decir,
P (T, z) = kT V1 ln (1 z) + 3Y (z)
1 z
n (T, z) = + 3 y (z)
V 1z
180
A cualquier valor de la ordenada le corresponde un z <1.
= z0 (T, n) para n3 . y (1) =
2 612
z T < Tc (n)
1 para n3 y (1)
n > nc (T )
ahora si metemos partculas la densidad sigue aumentando.
181
182
Condensacin de Bose-Einstein136
1e z
hn0i = ( ) = = .
e 0 1 1 e 1z
3
Usando esto en (47): n3 = y (z) + V hn0 i , es decir,
V y (z) y (z)
hn0i = nV =N 1
3 n3
" 3/2#
3c T
hn0i = N 1 3 = N 1 , T < Tc .
Tc
y (z)
hn0i = N 1 0, T > Tc .
n3
En definitiva:
3/2
N 1 TTc para T < Tc
hn0 i =
0 para T > Tc :
136
Tratamiento alternativo en Huang p.270.
183
Es decir, para T < Tc se produce acumulacin de partculas en el
estado fundamental p~ = 0 que, eventualmente, habr capturado a
todas para T 0. Se conoce a este fenmeno como condensacin
de Bose-Einstein (que ocurre en el espacio de momentos, es decir,
todas las partculas pasan a tener p = 0).137
Es una condensacin en el espacio de momentos (no en el real)
consecuencia directa de la simetra de la funcin de onda (no de las
interacciones entre partculas como en la condensacin real), pero
tiene una notable similitud formal con la condensacin ordinaria en
la transicin lquido-vapor (o la superfluidez, la superconductividad
o el cambio de fase paramagnetico-ferromagntico)138.
De hecho, el comportamiendo de hn0i N es como el del parmetro
de orden L G y otras magnitudes semejan a la situacin
137
Algunas referencias interesantes:
no puede ocurrir en d = 1, 2 : Hohenberg, Phys. Rev 158, 383 (1967)
puede presentarse en d = 1, 2 si rotaciones o campos gravitatorios: Widom, Phys.
Rev. 176, 254 (1968)
d = 1 con atraccin tipo : Ioriatti et al., Am. J. Phys. 44, 744 (1976)
138
En todas ellas la fase condensada tiene menos simetras que la anterior: hay una ruptura
de la simetra. Esta variacin de la simetra viene gobernada por un parmetro de orden.
184
139
en un cambio de fase de primer orden.
Por ejemplo, recordemos
1
la 1 ecuacin de estado, P = kT V ln (1 z) + 3 Y donde
a
1
V ln (1 z) es despreciable para z 6= 1 pero tambin si 1 z =
O (1/N ) , pues se hace O V1 ln N 0, de manera que se tiene
V
3 Y [z0 (T, n)] para n < nc
P (T, n) = kT
Y (1) para n > nc
Y [z0 (T, v)] para v > vc
P (T, v) = kT 3
Y (1) para v < vc ,
luego las isotermas presentan una seccin plana para 0 < v < vc ,
como en condensacin ordinaria:
139
Me refiero a hechos como el que describo a continuacin, y no a que un gas de bosones
bajo la accin de un campo gravitatorio experimenta una separacn de fases similar a la
condensacin ordinaria (Lamb & Nordsieck, Phys.Rev. 59, 677).
185
La seccin plana refleja que las partculas en el estado fundamental
(fase condensada) no contribuyen a la P 140 Notar que de llevar
la analoga hasta sus ltimas consecuencias, se seguira que estas
isotermas describen equilibrio entre tal fase condensada y una fase
gaseosa de volumen vc (T ) , y que la fase condensada tiene v = 0,
es decir, densidad infinita!
140
A medida que van cayendo en el estado fundamental, esto es para n > nc , no implica
un aumento de la presin.
141
Ej, a partir de la expresin para la entropa que luego deducimos.
186
El comportamiento a bajas T s (T 3/2) contrasta con el tpico
de fermiones (Cv T, T 0) y con el de fotones y fonones
(Cv T 3, T 0).142
de donde
Y (z0 )
para
n3 ln z0 T > Tc
s2
=
k5
Y (1) para T < Tc ,
n3
142
Resumiendo:
a a
P = n2 , s=
n T T n
notando que
a n T
= 1
n T T a a n
187
s 0 para T 0, de acuerdo con el postulado de Planck
(propiedad que no tiene el modelo clsico gas de Boltzmann
(cuya validez no incluye bajas T s);
Discusin
Sin embargo, no est claro que este punto de vista que encontraris
en algunos libros pueda seguir defendindose:
188
mediante una ingeniosa trampa magntica:145
145
Esta fenomenologa se describe, por ejemplo, en Physics Today, Agosto 1995, p.17.
189
Termodinmica del He lquido
Termodinmica del helio
Modelo de Tisza.
Introduccin
El resto de esta leccin tiene motivacin mltiple:
discutir la aplicacin del modelo gas ideal de bosones pues el
He a bajas temperaturas delata a nivel macroscpico los efectos
cunticos que se estn produciendo a bajas temperaturas en
su seno.
describir un tipo comportamiento general que tienden a pre-
sentar muchos sistemas materiales cerca 0K, y repasar as la
fsica de un sistema singular, 4 He
introducir un interesante mtodo para analizar la consecuen-
cia de las interacciones entre partculas, basado en el concepto
de cuasi-partcula
hacer evidente que no todava teora microscpica para el
comportamiento superfluido del 4He, a pesar de las maravillo-
sas ideas de Landau,
190
y de que el fenmeno es conocido hace ms de 50 aos
146
Para otros detalles, puede verse, por ejemplo: Special Issue: Helium-3 and Helium-
4, Phys. Today, February (1987); Superfluid 3 He: Theory and Recent Experiments, E.V.
Thuneberg & J.P. Pekola, Europhys. News 22, 3 (1991)
191
Termodinmica del helio
Puede verse una introduccin a las propiedades del He en Thermodynamics and Statistical
Mechanis de Finkelstein, p.64, y tratamientos ms detallados en J. Wilks, The properties
of Liquid and Solid helium, Clarendon Press 1967 (502-WIL).
192
El gas ms difcil de licuar y de solidificar: No se hace lquido hasta
enfriar por debajo de unos 5K a P normal, y permanece lquido a
T s menores (extrapolando, incluso en 0K); es necesario aplicarle
fuertes presiones para transformarlo en slido.
Se entiende, ya que:
6
6 r
(r) = 5 67 10 exp 21 5 1 08 ;
r
193
el de masa ms ligera entre los nobles, que implica mayor mo-
vilidad y dificultad de localizacin en una red;149
principio variacional = configuracin ms probable en 0K, la
que minimice E total (posiciones de equilibrio), pero principio
incertidumbre
x p ~ (x)2 (p)2 ~2
E = (p)2 2m
2
1 ~
E 10 K,
2m x x anchura pozo
194
Diagrama de fases del istopo 4 He:150
150
En realidad es algo ms complejo, pues tambin se presentan tres fases slidas (HCP,
hexagonal con mximo empaquetamiento, la nica estable a bajas T s, BCC, que ocupa una
pequea regin junto a la curva de fusin, y la FCC que aparece a presiones superiores a las
mostradas en la figura):
HCP solid
BCC solid
Pressure (bars)
25
Liquid helium I
Liquid helium II
critical
point
Gas
Temperature (K)
195
la curva de presin de vapor presenta un punto crtico para
Tc = 4 21 K y Pc = 2 6 atm (aproximadamente)
196
La forma de CV recuerda la discutida en la condensacin B-E:
197
Podemos definir
hn0i +
Tr eH a a
0 0
= lm = lm H
,
N N N N Tr (e )
que es una medida del nmero de ocupacin del estado funda-
mental. Para 6= 0 dara cuenta de una ocupacin macrosc-
pica del estado fundamental como en la condensacin B-E.
Hacer Mecnica Estadstica de bosones con interacciones!(*)
198
que se extendera por todo el fluido; los tomos se mo-
veran entonces coherentemente, sin poder perder E por
choques con otros, pues estn en el estado de mnima E.
Hoy sabemos que la de superfluidez requiere ciertas in-
teracciones entre tomos.
199
Modelo de Tisza
Modelo de los dos fluidos
Segn el modelo:
entropa nula
viscosidad nula (fluir sin resistencia por capilares, r
102 cm)
200
superfluido: formado por partculas en estado fundamental =
perfectamente ordenadas, luego no transportan entropa,
y no pueden interaccionar (es decir, intercambiar momen-
to), a menos que sean excitadas, luego no hay viscosidad. Lo
peculiar es tener partculas en el estado fundamental sin estar
en una red.
componente normal: formado por partculas excitadas, pueden
chocar con paredes, intercambiar momento y frenarse, como
fluidos ordinarios
Es forma (imprecisa) de hablar: no tiene sentido divisin
real en partculas normales y superfluidas; slo se propone
considerar dos componentes, como matiza luego la teora
microscpica.
Sin otras hiptesis, pueden entenderse cualitativamente
(incluso cuantitativamente a veces) algunas de las propie-
dades ms extraas del He-II:
Observacin:
201
capilares (como, por ej, en 1), y el movimiento de un objeto en el
seno del fluido (2)
(1) (2)
202
la ley de Poiseuille,
flujo de masa 1 P,
203
1111111111111
0000000000000 0000000000000
1111111111111
0000000000000
1111111111111
0000000000000
1111111111111 0000000000000
1111111111111
0000000000000
1111111111111 0000000000000
1111111111111
0000000000000
1111111111111
0000000000000
1111111111111
0000000000000
1111111111111 0000000000000
1111111111111
0000000000000
1111111111111 0000000000000
1111111111111
0000000000000
1111111111111 0000000000000
1111111111111
0000000000000
1111111111111
0000000000000
1111111111111 0000000000000
1111111111111
0000000000000
1111111111111
0000000000000
1111111111111 0000000000000
1111111111111
0000000000000
1111111111111 0000000000000
1111111111111
0000000000000
1111111111111
000000000
111111111
0000000000000
1111111111111
TB
0000000000000
1111111111111
0000000000000
1111111111111 TA
0000000000000
1111111111111
000000000
111111111
0000000000000
1111111111111
0000000000000
1111111111111
11111111
00000000
0000000000000
1111111111111 0000000000000
1111111111111
0000000000000
1111111111111 0000000000000
1111111111111
0000000000000
1111111111111
0000000000000
1111111111111 0000000000000
1111111111111
0000000000000
1111111111111 0000000000000
1111111111111
0000000000000
1111111111111 0000000000000
1111111111111
0000000000000
1111111111111
0000000000000
1111111111111 0000000000000
1111111111111
0000000000000
1111111111111
0000000000000
1111111111111 0000000000000
1111111111111
0000000000000
1111111111111 0000000000000
1111111111111
0000000000000
1111111111111 0000000000000
1111111111111
0000000000000
1111111111111
TB < TA
PB < PA
Lo contrario:
204
(1) (2) (3)
156
J.F. Allen & H. Jones, Nature 141, 243 (1938)
205
Las medidas relacionadas con los efectos anteriores pueden
explicarse cuantitativamente con facilidad:
206
3. tercer sonido: excitaciones en las que slo contribuye el super-
fluido, ej, perturbacin entre lminas paralelas muy prximas
de forma que slo tenemos superfluido entre ellas. stas se
propagan a velocidad c3 y sirven para confirmar la teora de
Tisza.
4. cuarto sonido: Cuando el He II se mueve sobre una sustancia
slida muy porosa que slo permita movimiento del superfluido
2 1 2
2 2 + 1 2 T = 0
c1 t
2 1 2 2
T + 2 2 = 0
2 2
c2 t t
c1
T 0 K: 1 , 2 0 c2 , luego se desacoplan
3
y, en efecto, 1o y 2o sonidos son ondas de densidad y trmica,
respectivamente.
a T s finitas, las ecuaciones estn acopladas; para resolverlas po-
demos probar:
h i h i
~
= 0 +1 exp i t k ~r , T = T0 +T1 exp i t k ~r ~
207
Experimento de Andronikashvili: determinacin emprica de
ecuacin de estado n /s a lo largo de la curva de presin de vapor:
208
209
Excitaciones elementales
Excitaciones elementales. Caso del 4 He.
Fotones. Gas en movimiento uniforme.
Rotones. Gas de rotones en movimiento.
Velocidad crtica de superflujo.
Vrtices. Teorema de Feynman.
Superconductividad.
Modelo de vrtices.
Introduccin
210
El estudio de este concepto y aplicaciones es el objetivo del
resto de esta leccin160
Excitaciones elementales
211
= puede hablarse de excitacin elemental161:
E = p F
excitacin asociada con
part. exc. de energa p
E = F p
hueco correspondiente
212
Para hacer ME, es necesario el espectro de las excitaciones ele-
mentales, que depender del estado fundamental y de las in-
teracciones
Ejemplo de antes vale para 3 He (ferminico), pero concepto apli-
cado en:
semiconductores y aislantes: espectro presenta saltos ener-
ga
superconductores (ver luego analoga superfluidez-superconducividad)
sustancias ferromagnticas: ondas de espn
solitones: excitaciones no-lineales de gran amplitud que
mantienen identidad incluso en colisiones
defectos de la materia, sean ferminicos o bosnicos, ...
Primera aplicacin: al 4He por Landau, que argumenta:
excitaciones en lquido dbilmente excitado han de poder
aparecer y desaparecer de una en una,
el momento angular en sistema de bosones slo puede va-
riar en no entero, luego
las excitaciones elementales en sistema de bosones
han de ser bosones
Se observa en el 4 He que: CV T 3 , como corresponde a
fonones, para T 0 K, luego
Las excitaciones elementales a muy bajas T s, es decir, asociadas
a pequeos p o grandes = h/p : son
fonones,162 esto es, ley de dispersin o espectro = cp
(c = velocidad propagacin del sonido en helio a 0K)
En definitiva, Landau propone:
fluido normal (en primera aproximacin) = gas ideal de
fonones en movimiento relativo al superfluido, siendo ste
ltimo el responsable del resto de masa y momento del
lquido He-II.
162
El razonamiento lleva a esperar que excitaciones elementales ferminicas en 3 He, donde
CV T
213
1a consecuencia: calores especficos
214
Gas en movimiento uniforme
Qu propiedades hay asociadas al helio movindose respecto de las
partculas en estado fundamental.164 Cual es la relacin entre el
momento total del gas y la velocidad v del movimiento de su masa.
215
una partcula con energa (~p) es:
1 +1 bosones
hn (~p)i = , = 1 fermiones
exp ( + + ~ p~)
0 M-B
(48)
V
RR 1
h3 0 0 m (c cos ) p2 d p 2 sen d
v= e 1
V
R R 1 2
3
h 0 0 em 1 p d p 2 sen d
donde m + ~ ~p = pc 1 + c cos .
Haciendo aqu cos = y p 1 + c
= p, se sigue:
R1
3
1 1 + c d
v=c R 3 =
1
1 1 + c d
166
es el ngulo que forman las velocidades de las partculas con el eje z, esto es con la
direccin que marca ~v
216
que atribuye significado a .
1
multiplicadores de Lagrange: = = ~ = ~v
kT kT
asociados con conservacin de no part energa m. total
P P
nii = E, nip~i = P~ ,
queda
1
hn (~p)i =
exp [ ( ~v p~)] 1
A partir de la distribucin en K:
notamos que
hni > 0 = ( ~v p~) > 0
= v < c,
= pc
217
X X
P = Pz = hn (~p)i pz = hn (~p)i p cos
p
~ ~p
Z Z
V p cos
= 3 v
p2 d p 2 sen d ,
h 0 0 exp pc 1 c cos 1
cos = , p 1 + = p
c
de modo que
Z Z h i4
2V p3 d p v
P = 3 1 cos cos sen d
h 0 ep c 1 0 c
16 5 v/c2
=V
45h3c3 4 (1 v 2/c2 )3
y la energa total:
X
E= hn (~p)i pc
p~
Z Z h i4
2V c p3 d p v
= 1 cos sen d
h3 0 ep c 1 0 c
2
4 5 1 + 31 (v/c)
=V
15h3c3 4 (1 v 2/c2 )3
218
Luego en el gas tenemos bien definido un momento y una velocidad,
luego parece lgico definir una masa inercial del gas de fonones.
Interpretemos P/v como masa inercial del gas de fonones;
se sigue la densidad de masa inercial asociada:
P 16 5 1
= = .
vV 45h3c5 4 (1 v 2/c2 )3
16 5k 4 4
0 = T
45h3c5
El de Andronikashvili establece
n = He (T T )5 6 , T & 1K.
219
No es posible hacer medidas a T inferior con este mtodo (basado
en viscosidad), pero medidas indirectas (basadas en la entropa y
en la velocidad del segundo sonido: Klerk, Hudson & Pellman 1953)
conducen a buen acuerdo:
Consecuencias:
otra vez, la idea de Landau (slo fonones) conduce a buen
acuerdo con datos a muy bajas T s
hay que explicar ese cambio que (como para el CV ) se produce
a unos 05K
probablemente: la ley de dispersin = pc deja de ser lineal en
p, pasando a tener importancia otros tipos de excitaciones
elementales (desaparece el comportamiento fonnico que
tiene ese espectro)
220
Otras propiedades de un gas en movimiento
y se sigue la presin
kT 4 5 k 4T 4 1
P= ln =
V 45h3c3 (1 v 2/c2 )2
y la entropa170
16 5 k 4T 3 1
S=V
45h c (1 v 2 /c2)2
3 3
170
Usando la expresin de la energa obtenida anteriormente, tambin se tiene
16 5 k 4 T 4 1
E + PV = V
45h3 c3 (1 v 2 /c2 )3
221
Otra consecuencia es que
1 1
hn (~p)i = h i , hn0 (~p0 )i = h i
~v~
p 0
exp + kT exp 0 + kT0
~v p~ 0
+ 0 + .
kT kT0
Puesto que:
~v p~
p = 0
1 v 2/c2
se sigue que
p
T = T0 1 v 2/c2 (Planck 1917, Einstein 1917),
222
223
Otras excitaciones elementales a temperaturas y momen-
tos mayores: Rotones
Hechos:
Rotones
n es
no
Fo
p p
0
224
A estas excitaciones nuevas que aparecen entorno a p0 Landau las
llam rotones, pues imagin que eran debidas a perturbaciones lo-
cales de carcter rotacional en el fluido.
1
hn (~p)i = .
e 1
171
Landau no us este argumento; de hecho, su naturaleza bosnica resulta irrelevante.
225
A las T s de inters, e 1,
experimentalmente, /k 86 K,de modo que e > e8 6 =
5432 1
X
X
X
[Q1 (V, T )]N
N
= z QN (V, T ) = QN (V, T ) =
0 N!
N =0 N =0 N =0
" #
X
= exp [Q1 (V, T )] = exp exp ()
p
pero z = e = 1, luego172 N = ln .
kT
Es P (V, T ) V ln , luego
P X
= N = ln = exp ()
kT p
X
172
N es controlado por T, no por . Adems N = hn(~
p)i = ln .
p
~
226
P V
R
Usando aqu la parte rotn, y ~p : h3 d~p :
Z ( " #)
2
P V 1 (p p0)
=N = 3 exp + 4p2 d p
kT h 0 kT 2
Z
4V
notar: esto es lo mismo que 3 n () p2 d p
h 0
Z " #2
1/2
4p20 V (2kT )
e/kT (2kT )1/2
2
= d x ex 1+ x
h3 p0 (2kT )1/2 p0
R 2
En definitiva, la integral es d x ex = 1/2 , de modo que
la contribucin de los rotones en el rango de T s de inters es
227
Se sigue que:
As:
228
la integral sobre n () conduce a P~ = 0 (pues se refiere a K0
donde el gas est en reposo global)
para la otra integral, poniendo pz = p cos , ~v ~p = vp cos :
Z Z
1 1 n 3
0 3 vp cos p cos d p
vh
Z Z
1 n 2
= 3 p2 cos2 p d p 2 sen d
h
Z
4 n 4
= 3 p dp
3h 0
173
Comprobar que conduce a la expresin que obtuvimos anteriormente para fonones por
un mtodo distinto.
174
Antes vimos al calcular N que
Z Z
4V 4V 2 4V p20
N= 3 n()p2 dp = 3 e/kT e(pp0 ) /2kT p2 dp = 3
(2kT )1/2 e/kT
h 0 h 0 h
donde la ltima igualdad se obtiene viendo que el rango de parmetros de y p0 y
2
T en el que trabajamos hace que n() sea muy picuda, n() = e/kT e(pp0 ) /2kT
/kT 1/2
e (2kT ) (p p0 ). Esto mismo se puede usar ahora dando el resultado requerido
229
Usando datos, confirmamos la situacin imaginada arriba:
roton. Como la masa efectiva del fluido normal va como 1/T entonces hp2 i = hp20 i =
6 f (T )
= hp2 i/2m 6= f (T ) pues m = cte bien definida.
230
desaparezca la superfluidez:
231
En otro caso, es decir, si se cumple el criterio de Landau,
v < (/p)mn vc,
no puede aparecer componente normal, luego ha de mantenerse
el carcter superfluido.
Implicaciones:
232
Landau (1941-47) no toma en serio esta posibilidad, conduciendo
sus estudios en otra direccin sin caer en la cuenta de que sus
excitaciones elementales han de justificarse como consecuencia
de esas interacciones
1
p/~ - 026 : parte lineal con pendiente 237 m/s, pr-
xima a velocidad sonido en He c 238 m/s;
179
Ver Huang p.383 y siguientes!
233
luego mximo (maxn) para /k = 1392 K y p/~ =
1
1 12 , transicin suave hacia:
p
= 8 65 K, = 192 1, = 0 16 mHe (T = 11 K )
k ~
que confirma la hiptesis de Landau.
luego cambio en 2a derivada que sugiere otro mximo...
234
otras observaciones sugieren falta de viscosidad (sin la
que es difcil intuir la turbulencia en un fluido ordina-
rio)
El rgimen turbulento en este caso debe ser compatible
con la condicin de Landau.
235
Bibliografa
~ es
accesible en el macroestado N, V, E, P
X
~ =
N, V, E, P W ({ni})
{ni }
donde
W {ni } = no microest distintos asociado con conjunto dado
de nos de ocupacin {ni } ,
la suma esta restringida por las leyes de conservacin:
X X X
ni = N, ni i = E, nip~i = P~
Entonces,
X
S = k ln W ({ni })
{ni }
236
[2] En lugar de volver a calcular directamente esta integral, recor-
demos que Z
4 2 N
n () p d p =
h3 0 V
237
a esas T s, pueden ser importantes otras exc elem, como ve-
remos
238
Tema 12: Fenmenos magnticos II
Transiciones orden-desorden: aleaciones binarias
Modelo Lenz-Ising y variantes: Modelo de Heisen-
berg, cristal magtico, gas reticular, red de neuronas
Soluciones -matricial y combinatorial- en dimensin
uno
Teorema de Peierls: existencia de cambio de fase en
dimensin dos.
Solucin en dimensin dos.
Soluciones en aproximacin campo medio (d = ).
Teora Curie-Weiss del ferromagnetismo.
Introduccin
En temas anteriores, al estudiar el paramagnetismo, vimos que esta
fase no es estable a bajas temperaturas y aparece un cambio de fa-
se hacia estados ferromagnticos, antiferromagnticos o vidrios de
espn. Esto ocurre porque a bajas temperaturas las interacciones
van a ser muy importantes luego el tratamiento de un sistema de
espines como sistema ideal no es posible. En este tema vamos a
estudiar esta transicin y nos va a servir para introducir ciertas
tcnicas matemticas utilizadas en el modelado de sistemas com-
plejos (soluciones exactas, aproximaciones campo medio, etc.) que
han jugado papel esencial en teora de cambios de fase.
La fase paramagnetica es una fase desordenada con los espines
orientados en todas las direcciones mientras que la fase ferro o an-
tiferro son fases ordenadas con los espines orientados en ciertas
direcciones. Comenzamos primero estudiando de forma genrica la
transicin de fase orden-desorden.
2
Sea red regular cuyos nudos (intersecciones) estn ocupados
por objetos de dos especies distintas, ej, in tipo A o in tipo
B. En dimensin dos, una configuracin posible es:
| | | |
A A B A
| | | |
B A B B
| | | |
A B A B
| | | |
B B A A
| | | |
3
1. Al-Zn, Cu-Ti, Ni-Si, Au-Pt,... tendencia agrupamiento, como
modelo 1 (pero complicada pues la red cbica-simple en la
que pensamos es casi una excepcin, hay nudos vacos, y
deformaciones y defectos donde quizs la tendencia es a otro
tipo de orden, o los segundos vecinos pueden interaccionar,
etc.)
2. Au-Cu, Au-Cd, Mg-Cd, Cu-Zn,... presentan tendencia a for-
mar super-redes (tambin complicada por otros factores)
4
No trozos en que se divide fase rica en As y la irregularidad
tpica de stos puede caracterizarse mediante longitud de corre-
lacin, , que ser f (T ) : es pequea, igual o proporcional
al radio de la esfera para T = 0K, y crece para T TC (de
hecho, diverge en TC ).1
Para T > TC , slo hay una fase homognea, mezcla de As
y Bs, pero orden de corto alcance, es decir, agrupamientos
microscpicos (ej, del orden de unos pocos espaciados de la
red)
5
donde A (~ri) = potencial qumico local asociado con la especie
A. La primera suma se extiende sobre qN/2 parejas de vecinos
prximos donde q es el nmero de coordinacin de la red.
Introducimos variables ocupacin tales que:
1
+1 (A) cA
i = 2 (1 + si )
si = = B 1
1 (B) ci = 2 (1 si )
donde
1
2Jij AB (rij ) [AA (rij ) + BB (rij )]
2
y
X1
Hi [BB (rij ) AA (rij )] + A (~ri ) B (~ri )
4
j(6=i)
X
qN/2
X
H = J si sj H si ,
i,j i
i6=j: vp
6
H0 = 0 (reescalamiento de la energa)
isotropa, de modo que Jij = J i, j, con
hMi = (A B ) N,
donde
X
N
M= si
i=1
4
Esto es as porque la P = 0 y en el hamiltoniano hay un trmino restando proporcional a
lo que hemos llamado H. Si esto lo interpretamos como el campo magntico y lo consideramos
variable independiente se tiene algo como A HM que no es ms que la energa de Gibbs de
un sistema magntico. En el caso de que hablemos de mezcla binaria H no es una variable
independiente y por lo tanto tenemos A.
5
Al escribir ZN (T, H) estamos implicando, adems del H, una energa trmica kB T
que puede interpretarse consecuencia de las vibraciones de la red (o sistema de fonones)
coexistiendo con los grados de libertad de ocupacin, representable como bao trmico, de
modo que tendencia al orden, medida por J y H, y tendencia a la aleatoriedad, medida
por T.
7
es no relativo de nudos ocupados por As y Bs en una confi-
guracin. Es decir, hMi es el promedio de M sobre todas las
configuraciones, cada una con su peso:
X X
hMi = M exp [ H ({si })] = kB T ln ZN (T, H)
s1 =1 sN =1
H
hMi =
6 0 (A B 6= 0) para H 0, T < TC .
8
NAB = NBA = no total de parejas AB vecinos proximos
Tenemos inmediatamente que6
qNA = 2NAA + NAB
qNB = 2NBB + NAB
Se sigue que:
NB = N NA (de N = NA + NB )
NAB = qNA 2NAA (de la 1a igualdad)
NBB = 21 qNB 12 NAB (de la 2a igualdad)
= 21 q (N NA ) 21 (qNA 2NAA ) (de las anteriores)
= 21 qN 21 qNA 21 qNA + NAA
= 12 qN qNA + NAA
de donde:
P 1
hiji si sj = NAA + NBB NAB = 4NAA 2qNA + 2 qN
P
i si = NA NB = 2NA N
luego
1
H (NA , NAA ) = J 4NAA 2qNA + qN H (2NA N )
2
1
= 4JNAA + 2 (qJ H) NA qJ H N
2
6
En efecto, tomemos un nudo A cualquiera y lanzamos uniones a sus vecinos prximos:
habremos dibujado q lneas. Repetimos con todos los As: habremos dibujado qNA lneas.
Pero stas pueden tambin contarse notando que:
entre cada AA habr dos lneas 2NAA
entre cada AB habr una lnea NAB
no habr lnea entre las parejas BB
Esto justifica la 1a igualdad; un argumento similar lleva a la 2a
9
y, finalmente, la funcin de particin es7:
X
G
ZN (T, H) = e = g (NA , NAA) exp [H (NA , NAA)]
NA ,NAA
X
N
( 12 qJH)N
=e e2(qJH)NA
NA =0
X
g (NA , NAA) e4JNAA
NAA
7
La suma es sobre todos los valores posibles de NA , primero, y luego sobre todos los NAA
compatibles con el NA en cuestin.
10
11
Modelo Lenz-Ising y variantes
con la interpretacin
+1 espn up
si = ,
1 espn down
12
M 6= 0) o antiferromagnticas a bajas T s.
Este modelo puede resolverse exactamente para d = 1 y 2.
Ising (1925) lo resolvi para d = 1, con el resultado de que no hay
transicin a T finita, es decir, TC = 0.
13
Nota: La ausencia de transicin en d = 1 puede entenderse: Sea
estado ordenado, con todos los espines down. Es posible que una
fluctuacin trmica del equilibrio invierta una serie de espines con-
secutivos en la cadena. Esta inversin no puede ser contrarrestada
por efecto de los espines no invertidos en los extremos, pues hay
tendencias misma intensidad y distinto signo:
14
cuadrada, con vecinos prximos, H = 0. Hoy se considera una de
las contribuciones ms importantes en fsica terica de las ltimas
dcadas: por primera vez y casi nica se llegaba a reproducir
un cambio de fase y fenmenos crticos en el lmite N para
una interaccin razonable; adems se sabe ahora que es extraordi-
nariamente realista para multitud de fenmenos.
15
(D)
Las variables ~si pueden interpretarse como espines cuanti-
zados,10 capaces de presentarse en 2s + 1 estados, con s =
no cuntico de espn; el modelo requiere entonces tratamiento
cuntico, con lmite clsico para s .
Modelo de Heisenberg: d = D = 3.
La versin clsica (s ) , fue estudiada por Heller & Kra-
mers (1934) a bajas T s.
La versin cuntica (s finito; generalmente,
p se normaliza la
magnitud de cada espn dividiendo por s (s + 1)), y repre-
senta aceptablemente algunos materiales,11 como compuestos
ferromagnticos del europio (EuO, EuS) y los materiales anti-
ferromagnticos RbMnF3 y KMnF3 (incluso MnO, CoO y otros
xidos) tomando Jij < 0 i, j, que favorece la alineacin anti-
paralela de espines (formacin de super-redes entremezcladas
tal que dado un espin en una red, todos sus vecinos prximos
estn en la otra red) como en estados antiferromagnticos.
16
P
restringida por la condicin s2i = N
d = 1 para cualquier D
Cristal Magntico
Nos podemos preguntar el por qu el hamiltoniano de Ising o Hei-
senberg es una buena hiptesis para describir un cristal magntico,
en qu pudo basar Lenz su intuicin para proponerselo a Ising.
17
Cmo es la interaccin entre ~si y ~sj vps? La mecnica cuntica
nos dice que es:
Kij = energa directa o culombiana entre los espines
Jij = energa de intercambio, de naturaleza cuntica
para espines antiparalelos, es decir, para estado
Kij Jij , donde +:
singlete , antisimtrico, de espn total S = 0
para espines paralelos, es decir, para el estado
:
triplete , simtrico, de espn tot. S = 1
18
triplete s = 21 , S = 1 : = 41 , = 2Jij 14
singlete s = 21 , S = 0 : = 43 , = 2Jij 34
luego
1 3
= 2Jij + 2Jij = 2Jij
4 4
QED15
En este contexto, Ising equivale a reemplazar el producto
Por ltimo, resulta que Jij , aunque es del mismo orden de magni-
tud que Kij , decrece mucho ms rpidamente con la distancia y, de
hecho, es relativamente pequea para vecinos segundos.
15
De hecho, esto es una prueba de consistencia, pero puede demostrarse, salvo
una cte. aditiva (K. Stevens, Phys.Rep. 24C, 1 (1976) y apndice 1 de tesina de
Parra), partiendo del hamiltoniano para dos tomos en el q figuren los tnos. de enrga
correspondientesa los electrones de las capas atmicas incompletas. En definitiva,
se entiende el hamiltoniano de Heisenberg.
16
Se sabe (Matsubara & Matsuda 1956, Betts et al. 1968-70) que tambin es interesante el
modelo con dos trminos en el hamiltoniano, six sjx + siy sjy , pues simula ciertas propiedades
del helio y de sustancias ferromagnticas aislantes.
19
Otras variantes: Gas reticular.
Propuesto por Yang & Lee (1952) para estudiar fluidos y sus tran-
siciones:
20
y el clculo se reduce, como en Ising, al de la degeneracin g (NA , NAA).17
(T, N, ) = eN P (P = presin)
X
N
= z NA ZNA (T, N )
NA =0
X
N X
NA
= e g (NA , NAA ) e0 NAA .
NA =0 NAA
Para acabar de explotar esta relacin, notemos que, para una mezcla
binaria (como vimos):
17
N hace el papel de volumen total accesible al gas (en unidades de volumen de una celda
de la red);
NAA = no total de parejas de nudos vecinos prximos ocupados;
g (NA , NAA ) = no de formas distintas de distribuir NA partculas indistinguibles en N nudos,
de modo que se obtengan NAA parejas de partculas vp;
la suma extendida a todos los posibles valores de NAA compatibles con los valores dados de
NA y N
21
En definitiva:
mezcla binaria gas reticular Ising ferromagn.
NA NA N+
NB N NA N
AA 2AB + BB 0 4J
P N + 12 qJ H
G
22
23
Solucin Ising d=1: mtodo matricial.
El caso d = 1 es muy interesante, a pesar de no presentar cambio
de fase (como veremos), puesto que tiene solucin exacta y per-
mite estudiar la naturaleza de mtodos para tratar problemas ms
complejos.
de modo que:
24
X
N
1 X
N
H (~ ) = J ii+1 H (i + i+1)
i=1
2 i=1
X N
1
= Jii+1 + H (i + i+1)
i=1
2
y se tiene:
X
Z (H, T ) = eH(~ )
{~ }
N
( )
X X X 1
= exp Jii+1 + H (i + i+1) .
1 =1 N =1 i=1
2
e(J+H) eJ
P= ,
eJ e(JH)
19
se sigue:
X
Z (H, T ) = 1 PN 1 = traza PN = N N
+ + ,
1 =1
18
Es decir, los posibles elementos de matriz son
h+1 |P| + 1i = e(J+H) ; h1 |P| 1i = e(JH)
19
En efecto, teniendo en cuenta la definicin de producto de matrices:
X X
Z (H, T ) = h1 | P|2 i h2 |P|3 i h3 |
| {z } | {z }
1 =1 N =1
=1 =1
|N 1 i hN 1 |P|N i hN |P |1 i
| {z } | {z }
=1 =1
25
donde son los valores propios de la matriz P, esto es las solu-
ciones de la ec secular:
(J+H)
e eJ
= 0,
eJ e(JH)
de donde:
q
= e J
cosh (H) e2J cosh2 (H) 2 senh (2J)
o, equivalentemente:
q
= eJ cosh (H) e4J + senh2 (H) .
1
= ln + + ln 1 + ln +
N + N
26
En definitiva:
q
1 J
ln Z (H, T ) +ln cosh (H) + e4J + senh2 (H) ,
N N kT
1
a (H, T ) = kT lm ln Z (H, T )
N N
q
= J kT ln cosh (H) + e4J + senh2 (H)
a (0, T ) = J kT ln 1 + e2J
= kT ln [2 cosh (J)]
calor especfico:
2a
c0V = T = k (J)2 sech2 (J) ,
T 2 N,V
que se comporta:
27
c0V /k 0.4
0.2
0
0 1 2
kT/J
es decir, nunca singularmente (como se esperara en cambio de
fase), sino que tiene pico ancho para kT /J 1, lo que indica
transicin gradual desde orden completo cuando T = 0 hacia
desorden completo a altas T s.:
a 1
m (H, T ) = = senh (H)
H T +
senh (H)
=q
e4J + senh2 (H)
que se comporta:
28
m(H ,T)
+1
T1
T2
T1 < T2
-1
m (H = 0, T ) 6= 0,
es decir, hay magnetizacin espontnea.
29
C. Itzykson & J.M. Drouffe, Statistical Field Theory, vol.
1, Cambridge Monographs on Mathematical Physics, 1990.
C.J. Thompson, Algebraic Derivation of the Partition
Function of a Two-dimensional Ising Model, Journal of
Mathematical Physics 6, 1392-1395 (1965)
T.D. Schultz, D.C. Mattis & E.H. Lieb, Two-dimensional
Ising Model as a Soluble Problem of Many Fermions, Re-
view of Modern Physics 36, 856-871 (1964)
T.T. Wu, Theory of Toeplitz Determinats and the Spin
Correlations of the Two-dimensional Ising Model, I, Phy-
sical Review 149, 380-401 (1966)
McKoy & Wu, The Two-dimensional Ising Model, Har-
ward Univ. Press 1973.
L.P. Kadanoff, Spin-spin Correlations in the Two-dimensional
Ising Model, Il Nuovo Cimento XLIVB, 276-304 (1966)
R.P. Feymann, Statistical Mechanics, Benjamin (1981),
p.127.
Hamiltoniano:
X
N X
N
H = 2J ~si ~si1 g ~ ~si,
H |~s|
i=1 i=1
30
Cadena unidimensional con espines de dimensin ar-
bitraria
Interaccin coulombiana en d = 1
31
J. Joffrin, Disordered Systems Experimental Viewpoint,
Ill-condensed Matter, Les Houches 1978, p. 63, North-
Holland.
Redes de neuronas
1 X
P
Jij =
N =1 i j
32
Adems todas neuronas j actan sobre una dada i con
intensidad Jij Puede resolverse por completo con tcnica
de rplicas.
En d = 2 :
X X
H = J1 sj,k sj.k+1 Jj sj,k sj+1.k
jk jk
Hamiltoniano: X
H= Jij si sj
hiji
33
prximos (bond dilution). Modela impurezas en sistemas
reales y es tambin muy interesante por la dependencia de
su comportamiento crtico con d, por ser una especie de
percolacin a T finita, etc.
Gas de Takahashi
Modelo X-Y
X
H= Jij sxi sxj + syi syj
ij
34
Inters: sx sy no conmutan. Si parmetro de orden es m =
P x
i si , no conmuta con H = estado fundamental no es
con todos los espines alineados.
35
Teorema de Peierls.
Veamos los argumentos de Peierls21 que, aparte de un inters hist-
rico, como se discuti, se han usado con xito en otros problemas.
Demuestra que debe haber cambio de fase en el modelo de Ising en
d=2
a
m= ,
H
36
luego a es simtrica dada T,
a (H) = a (H) ,
y se tiene la antisimetra
m (H) = m (H) .
En definitiva, ha de ser:
Figura 1
37
Figura 2
0, con independiente de
m
1
25
hi i = Tr i eH Tr eH
38
Parte 1: Sea un cuadrado, y la condicin lmite especial i =
+1 en todo el contorno del cuadrado,26 como en este ejemplo:
39
En efecto, para formar estos polgonos de b lados, podemos partir de un nudo
cualquiera de la red, y podemos seguir en cualquiera de las direcciones del
espacio, luego hay 4N formas distintas de empezar en una red cuadrada. Cada
uno de los restantes lados puede trazarse, como mucho, en 3 direcciones, pues
no puede volverse por el camino recorrido, luego pueden formarse 4N 3b1
fronteras con longitud b siguiendo este mtodo. Puesto que cualquiera de los b
vrtices del polgono podra haberse tomado como punto de partida, se sigue
b 4N 3b1 /b, QED. (Como b 4, tambin puede escribirse b N 3b1 ,
pero no es necesario para el argumento.)
Para demostrar
N
hNi (1 ) , 0:
2
Sea {C} el conjunto de todas las ~ en las que aparece una
frontera dada, (j, b)
Sea C la frontera que se obtiene de C invirtiendo (es decir, i i)
todos los espines dentro de esa frontera.29
Entonces, para H = 0 (el caso en cuestin):
H (C ) = H (C) 2bJ,
P
eH(C)
{C}
= P = e2bJ
e+2bJ eH(C) hay correspondencia uno a uno entre {C} y {C }
{C }
29
Es decir, en la frontera grande del ejemplo se invierte tambin el espn central + .
41
En consecuencia:
" #" #
X b 2 X b hD
(j)
Ei X b2 4N3b1
hN i Xb e2bJ
4 j=1
16 b
b=4,6,... b=4,6,...
N X b
= b 3 e2bJ (a comparar con la expresin de arriba)
12
b=4,6,...
lm a
(H) lm a
(0) H.
30 1
Esto es, la serie converge para x 3 exp (2J) 1, es decir, T 2J (k ln 3) . En
este caso, puede hacerse la suma y se tiene
hN i x4 2 x2
2.
N 6 (1 x2 )
42
Pero sabemos que lm a
(H) = a (H) , luego
a (H) a (0) H, H 0,
a (H) a (0+ )
H
luego
a (H) a (0+ ) a
lm = lm+ lm+ m (H) me .
H0 H H0 H H0
43
Ms compleja es la extensin al caso de magnetizacin cons-
tante (como conviene a las mezclas binarias) con fases separa-
das. El formalismo adecuado entonces es el cannico (en lugar
del macrocannico que hemos usado hasta ahora). En la prc-
tica, slo hay que aplicar las relaciones generales entre los dos
formalismos.Por ej, la probabilidad de que magnetizacin =
m es
X X
m Pm eH = eN a (m )
1 N
donde P
1 si i i = m N
Pm =
0 en otro caso
44
espn entero, s = n : se llega a este mismo resultado, pero
la demostracin es ms complicada (Lebowitz & Gallavot-
ti 1971).
Se ha extendido al modelo de Heisenberg, con interaccin
isotrpica Eij = J ~si ~sj /s2, si los componentes de ~s no to-
man valores cuantizados, se sigue una solucin equivalente a
tomar s , y (para d = 2 y J > 0 interaccin ferro-
magntica) puede demostrarse que no hay transisin (Mermin
& Wagner 1966). Podemos plantear aqu un problema intere-
sante: cunta anisotropa es necesaria para que este modelo
presente magnetizacin espontnea:
45
Bibliografa
46
cin = m es
X X
m Pm eH = eN a (m )
1 N
donde P
1 si i i = m N
Pm =
0 en otro caso
Por ej, se demuestra, con espines + en los bordes, que los espi-
nes tienden a agruparse dentro de una gran frontera que, para
T s suficientemente bajas, tiene aproximadamente la forma de
un cuadrado.
47
(hacer) que cota superior para TC tiende a cero cuando n .
(Luego esta extensin no es muy interesante, salvo por mostrar
una de las limitaciones del TP: el modelo tiene transicin, pero
el TP establece una cota demasiado baja para TC .)
48
antiferromagntico (Griffiths, en Domb & Gree II, p. 5-13),
modelos no-reticulares y modelos cunticos, y ciertas propieda-
des especiales de modelos familiares:
49
Solucin del modelo de Ising en dimensin 2.
La ventaja del mtodo de las matrices de tranferencia es que se pue-
de generalizar fcilmente a dos o ms dimensiones. Consideremos
el caso de dimensin d=2. Tomamos una red bidimensional (con
m filas y n columnas) sobre un cilindro (condiciones de contorno
periodicas en columnas pero no en filas). La energa configuracional
del sistema de espines es:
XX
m1 n X
m X
n X
m X
n
H(S) = J si,j si+1,j J si,j si,j+1 H si,j
i=1 j=1 i=1 j=1 i=1 j=1
(3)
con si,n+1 = si,1.
50
La funcin de particin queda entonces como
X
Zn,m = exp(H(s))
{s}
" !#
X X
n
= exp {V1 (j ) + V2 (j , j+1)}
(8)
X 1 ,...,n j=1
= T (1, 2)T (2, 3) . . . T (n1, n)T (n, 1)
1 ,...,n X
= T n (1, 1)
1
T (, ) = exp[V1()] exp[V2(, )]
! !
X
m1 X
m X
m
= exp J si si+1 + H si exp J si si
i=1 i=1 i=1
(9)
que se puede poner en forma simtrica como (Ejercicio)
T (, ) = exp[ V1 ()] exp[V2(, )] exp[ V1 ( )]
2 2
! !
J X
m1
H X
m X m
= exp si si+1 + si exp J si si
2 i=1 2 i=1
! i=1
J X
m1
H X m
exp si si+1 + si
2 i=1 2 i=1
(10)
que es una matriz 2 2 . Entonces la funcin de particin queda:
m m
X
m
2
Zn,m = Tr (T)n = nj (11)
j=1
51
m (aunque no es una condicin necesaria) entonces
1 1
a(H) = lm lm ln Zn,m
m n mn " !#
1 1 X2m
= lm ln 1 + lm lm ln 1 + (j /1)n
m m m n mn
j=2
1
= lm ln 1 .
m m
(12)
El problema se reduce tambin a encontrar el autovalor ms grande
y = J.
52
La energa libre por espn es entonces
1 X
m1
1
a(H = 0) = ln(2 sinh 2) + lm 2k+1. (15)
2 m 2m
k=0
= J coth 2
ZZ
1 d1d2
1 + (sinh2
2 1)
2 cosh 2 (sinh 2)(cos 1 + cos 2)
2
0
(20)
53
Aqu la integral diverge logartmicamente en el origen (1 = 2 = 0)
cuando cosh2 2 = 2 sinh 2 (Ejercicio: demostrarlo). Hay una
singularidad, o cambio de fase, por lo tanto cuando
sinh 2c = 1 (22)
54
0
u(H=0)/J
-1
-2
0 1 2 3 4
kT/J
En la vecindad de k1 = 1 se tiene
55
2
c0V /k
0
1 2 3 4
kT/J
56
1
m0=m(H=0)
0.5
0
1 2 3 4
kT/J
57
Teoras de Campo medio
Los modelos de magnetismo descritos anteriormente son, en gene-
ral, difciles de resolver de forma analtica exacta. El modelo de
Ising slo est resuelto en 1 y 2 dimensiones. Cuando no podemos
trabajar de forma analtica podemos hacer simulaciones numricas
del sistema, o bien utilizar aproximaciones. Entre stas, estn las
teoras de campo medio y teora de perturbaciones.
c F
c I
[ \
C = I F , I F =
I = {i C /si |i j| = 1 j C }. (32)
58
La aproximacin de campo medio considera de forma exacta, en el
hamiltoniano, la interaccin entre espines en el cluster. El resto de
las interacciones se incluye en un campo medio en el hamiltoniano,
desconocido a priori, que se calcula de forma autoconsistente.
59
Interaccin a vecinos prximos:
Suponemos que, dado un espn en la posicin i, se tiene para la
interaccin de intercambio Ji,j = J, |i j| = 1. La hiptesis de
campo medio consiste en suponer
X
H(s) = Hi si (34)
i
con X
Hi = H + Ji,j sj H + Jqm, (35)
|ij|=1
60
Interaccin de largo alcance:
Consideramos ahora Ji,j de largo alcance: dos espines situados en
posiciones lejanas sienten la interaccin de intercambio. Sea Ji,j =
J
N i, j , siendo J intensivo ( que H(s)/N no diverje X
en el
lmite termodinmico N ) . Si llamamos x(s) = N 1 si, el
i
hamiltoniano queda
J JN
H(s) = x(s)2 HN x(s). (40)
2 2
La funcin de particin es
X n o
J X JN
Z= eH(s) = e 2 exp 2 x(s) 2
+ HN x(s) .
s s
(41)
Utilizando la transformacin gaussiana
Z + y2
a2 1 2 + 2ay
e = dy e , (42)
2
con 1
JN 2
a= x(s), (43)
2
la funcin de particin queda
J 1 Z +
y2 X 1
Z=e 2 dy e 2 exp{(JN ) 2 x(s)y + HN x(s)}.
2 s
(44)
1
Haciendo el cambio y = (JN ) 2 m, y realizando la suma sobre
configuraciones de espines se tiene
J JN 1 Z +
N
Z=2 e 2 (
)2 dm eN f (m), (45)
2
con
Jm2
f (m) = ln{cosh(Jm + H)}. (46)
2
En el lmite termodinmico, el nico trmino que va a contribuir a
Z es el trmino dominante en la integral, que va a ser el punto de
silla, determinado por la ecuacin
f (m)
= 0, (47)
m
61
que da como resultado:
m = tanh(Jm + H). (48)
Para campo magntico nulo (H = 0), la magnetizacin espontnea
es
m0 = tanh(Jm0). (49)
Esta ecuacin tiene solucin no trivial para T < Tc = J , de forma
k
que tenemos, al igual que para el caso de vecinos prximos, un
cambio de fase ferromagntico-paramagntico.
B. Aproximacin de orden 1
Esta aproximacin se conoce tambin con el nombre de aproxima-
cin de Bethe-Peierls. Trata con ms precisin la interaccin de un
espn dado con sus vecinos prximos. Ms concretamente considera
de forma exacta la interaccin de un espn si0 con sus q vecinos
prximos, y el resto de las interacciones en la red a travs de un
campo medio, que de nuevo obtenemos de forma autoconsistente.
Vamos a considerar de nuevo interacciones de intercambio slo a
vecinos prximos. Llamamos F al conjunto de espines que son ve-
cinos prximos a si0 y que constituyen su frontera. El hamiltoniano
de Ising queda de la forma
1 X X
H(s) = Ji,j sisj H si HBP (sc)
2
|ij|=1 i
X X
HBP (sc) = Jsi0 sj Hsi0 (H + h) sj , (50)
jF jF
62
Para determinar el campo h, imponemos que el promedio del espn
ha de ser igual en el centro y la frontera del cluster, es decir,
1X
hsi0 i = h sj i. (53)
q
jF
Por definicin,
1 X Z+ Z
hsi0 i si0 eHBP (sc) = (54)
ZBP s ZBP
c
63
Esta ecuacin tiene solucin no trivial h0 y por tanto m0 (dada por
Eq. 57) para T < Tc , con
2J
Tc = . (59)
q
kB ln q 2
64
0
Onsager
Bethe-Peierls
Bragg-Willians
u(H=0)/J
-1
-2
0 1 2 3 4
kT/J
Onsager
Bethe-Peierls
1 Bragg-Willians
m0
0.5
0
0 1 2 3 4
kT/J
Onsager
Bethe-Peierls
Bragg-Willians
2
/k V
0
c
0
0 1 2 3 4
kT/J
65
Tema 13: Introduccin a la teora de los cambios
de fase: teoremas de Yang y Lee
Fase termodinmica
Ceros de la funcin de particin.
Teoremas de Yang y Lee.
Consecuencias. Cambios de fase y clasificacin.
66
El estudio experimental de estos cambios de fase nos ha ensea-
do que vienen asociadas, casi por definicin, con discontinuidades,
divergencias u otros comportamientos noanalticos de ciertas mag-
nitudes termodinmicas tales como calores especficos y susceptibi-
lidades.
37
Vanse los libros de Thompson, Ruelle, Ma, Pathria y Huang.
67
Una duda razonable
H = H0 + (1, ..., N )
con
P p2i
H0 = i 2m es la energa cintica de las molculas libres
es la parte potencial que suponemos slo depende de po-
P
siciones y puede escribirse = i<j ij con ij potenciales
(efectivos) entre dos cuerpos.
Z p
1
QN (V, T ) = d~r1...d~rN e , = h2 /2mkB T ,
N !3N
X
(z, V, T ) = z N QN (V, T ) , z e/kT .
N 0
68
La presin y el volumen especfico vienen entonces dados por:38
1
P = ln (z, V, T )
V
1 1
= ln (z, V, T ) .
v ln z V
En consecuencia,
QN (V, T ) = 0, N > Nm (V )
(es decir, esos estados tienen energa infinita) y, a todos los efectos
prcticos, es un polinomio en z de grado mximo Nm (V ) :
X
Nm
(z, V, T ) = z N QN (V, T )
N =0
38
Notemos que v = V / hN i es aqu un parmetro independiente de V. La segunda ecuacin
se sigue entonces de que hN i = 1 z z
y v1 = n = hN i /V.
39
Es suficiente con que Nm sea creciente con V.
69
cuyos coeficientes son definidos reales y positivos.40
Sean
zk = zk (V, T ) , k = 1, 2, ..., Nm
X
Nm
QN (V, T ) z N = 0. (62)
N =0
Se tiene que:
40
Notemos que estos coeficientes slo pueden determinarse conociendo (si sabemos cal-
cular la integral QN ). En lo que sigue, determinamos propiedades generales importantes sin
llegar a determinarlos.
41
Por ejemplo, si xk = , , se tiene:
(x ) (x ) = x2 ( + ) x + ,
70
Pero los ceros de son los puntos singulares del potencial ter-
modinmico, V1 ln (z, V, T ) , y, en consecuencia, de P y de
1/v,
luego, tanto P como 1/v, son funciones analticas de z en una
regin del plano que incluye el eje real positivo y, en conse-
cuencia,
P es funcin analtica de v en una regin del plano v que incluye
el eje real positivo;
notad que, dado su significado fsico, los parmetros v y z = e/kT
son definidos reales y positivos en Fsica, luego
la presin
o cualquier otra magnitud, que siempre puede obtenerse a partir de sta en la
colectividad macrocannica,
71
se tiene (por ejemplo; hay ahora otros formalismos):
Z X
QN (V, T ) = d~r1...d~rN (1...N ) eH (1...N )
Conclusin:
72
Adems, el sistema ha de encontrarse, para valores dados de
V y T, en una fase bien definida a lo largo del eje real positivo,
no pudiendo aparecer en ste un cambio de fase o coexistencia
de fases.
Pero (63) sugiere que podra tenerse P/v O (1/V ) en al-
guna regin del diagrama termodinmico, de modo que P/v
sera ciertamente pequea
para un sistema macroscpico, para
el que V O 1024 tpicamente, y se seguira P prcticamen-
te constante (en escala macroscpica), tal como se observa en
una coexistencia de fases en el laboratorio.
Es decir, se plantea la posibilidad de que el formalismo refleje
cambios de fase slo en el caso de sistemas macroscpicos que,
como vimos, requieren ser representados en este contexto con
la condicin V (para evitar influencia de las condiciones
en los lmites o efectos superficiles). Ms explcito, los modelos
han de considerarse en el LT:
V
V , N , finito e igual al volumen especfico dado.
N
En estas condiciones, la ecucacin de estado paramtrica es:
1
P = lm ln (z, V, T ) (64)
V V
1 1
= lm ln (z, V, T ) (65)
v V ln z V
73
Sabemos que, para potenciales muy generales, Nm V / 3 o, al
menos, crece con V, y que los ceros son zk = zk (V, T ) , de modo
que se movern continuamente por el plano z al variar V.
Es decir, al crecer V, la distribucin inicial de ceros se va haciendo
ms densa. En estas condiciones, puede ocurrir que, para V ,
acabe eventualmente hacindose continua, de modo que aparezca
un punto z0 de acumulacin de ceros en el eje real, es decir, una
singularidad esencial aislada:43
74
Sea una regin R en el plano z que contiene un segmento del
eje real positivo libre de ceros de (z, V ) para cualquier valor
de V. En consecuencia, P 0 y P/v < 0 en R, de modo
que el sistema se encuentra en una fase homognea nica por
todo R.
Si admitimos la posibilidad indicada ms arriba, podramos
definir varias regiones R, R , ... de este tipo, sin solapamiento
entre ellas, cada una correspondiente al sistema en cada una
de sus posibles fases:
44
Un punto de acumulacin de ceros slo puede conducir a:
una funcin idnticamente nula, que no es nuestro caso, o
una singularidad esencial aislada como, por ejemplo, la funcin sen (1/z) , para la que
z = 0, es el lmite de la sucesin de ceros z = 1/n (n = 1, 2, ...)
(Por el contrario, tan (1/z) tiene polos en los puntos z = 2/n, n = 1, 3, ..., y, por
tanto, z = 0, punto de acumulacin de la sucesin de polos, es un punto singular esencial
no aislado) pues en un entorno de z = 0 hay otro tipo de singularidades (los puntos de la
sucesin anterior con n impar)
75
Teoremas de Yang-Lee
Primer Teorema de Yang-Lee, que tiene dos partes:
Primera parte:
Segunda parte:
76
Para que F (z) sea funcin continua de z slo nos falta probar que
su derivada (logartmica) permanece acotada V :
P
N m
z =z 1 + zQ1 + z 2 Q2 + ... = zQ1 + 2z 2 Q2 + ... = N z N QN
z z N =1
QN >0 P
N m P
N m
Nm z N QN < Nm 1 + z N QN = Nm
N =1 N =1
1 1 1 Nm por ej 1 V 1
ln = z < = const 3 = const 3 < .
ln z V V z V V
Sea una regin R del plano complejo z, tal que incluya un segmen-
to del eje real positivo y no contenga ceros para cualquier V
(incluyendo V ).
Se sigue que, z R, la funcin
1
F (z, V ) ln (z, V ) (67)
V
converge uniformemente para V hacia un lmite que
es funcin analtica en R.
77
En efecto, por ser QN 0, se tiene:
1 h Nm
i
|F (z, V )| ln 1 + |z| Q1 (V ) + + |z| QNm (V ) .
V
Pero, en el crculo D, es |z| < = no real, luego:
1
|F (z, V )| ln 1 + Q1 + + Nm QNm .
V
Ahora bien, la derecha de esta expresin est acotada
pues lo est para cualquier V finito y, por el primer
teorema (ya que QN y son reales), tiene lmite para
V por lo que |F (z, V )| est acotada.
Por ltimo, el primer teorema implica que F (z, V )
tiene un lmite dentro de D que es precisamente F (z)
(que es real y positivo).
Pero stas son las condiciones del Teorema de Vitali que
implica que, en estas condiciones, se verifica el segundo
Teorema de Y-L (convergencia uniforme hacia una funcin
analtica) en D.
A continuacin, podemos contruir en R otro crculo D con
centro en D, las condiciones de Vitali tambin se verifican
en D ;
as, sucesivamente, puede extenderse la validez del teore-
ma a toda la regin R.
78
Algunas consecuencias:
79
Como vimos, esto quiere decir que el sistema se encuentra
en una fase nica z y T que se correspondan con esta
situacin.
80
Se sigue de los teoremas Yang-Lee que el sistema puede
presentarse en dos fases, para z < z0 y z > z0 , respecti-
vamente.
P (z) ha de ser continua en z = z0 , como consecuencia
del 1er teorema, pero su primera derivada, 1/v, puede ser
discontinua:46
46
Problema: Si hay cambio en el volumen especfico, tal como se implica en este caso,
ha de ser tq 1/v (z) crezca cuando z aumenta pasando por z0 .
En efecto, de relaciones discutidas al estudiar el formalismo macrocannico, se tiene:
D E
(N hN i)2
1 1
= ln = 0
ln z v (z) ln z ln z V V
81
Eliminando z entre la ec. para P (z) y la ec. para 1/v(z) se
sigue la relacin P (v) con la forma indicada en el grfico
anterior.
A este tipo de CF se le llama de primer orden, con una
regin (vL , vG) de coexistencia de las fases.
Esta situacin la presentan tpicamente los fluidos para
T < TC .
82
Puede tambin ocurrir que la singularidad entre R1 y R2 sea
de tal naturaleza que 1/v, 1a derivada de P (v) , permanezca
continua pero la segunda, 2P/z 2 , presente discontinuidad:
83
Habr que definir varias regiones, R1 , R2, R3,... corres-
pondientes a distintas fases del sistema;
z0 corresponder a la coexistencia de dos fases, z0 a la de
otras dos, etc.
84
En definitiva, hemos mostrado que:
ej, podra ocurrir que al acercarse los ceros al eje real para
V no lo hicieran hacia un punto nico, sino hacia
un segmento;
todava se tendra un cambio de fase desde un punto de
vista matemtico en este caso, aunque la interpretacin
fsica de tal anomala sera comprometida;
Utilidad:
85
Se sigue de lo anterior que esas curvas de coexistencia
son el lugar geomtrico de las singularidades esenciales
aisladas z0 (T ) o, ms concretamente, P (z0 ), al variar T.
86
Lamentablemente, no es fcil en la prctica determinar z0 (T )
puede ser tan complicado como calcular la funcin de parti-
cin.
87
Fsica Estadstica. Relaci
on 1
October 9, 2009
4. Estime el tiempo necesario para que una molecula de aire en una habitacion
llegue a una distancia de L = 5 metros de donde esta. Considere que no hay
movimiento macrosc opico del aire, que la temperatura y la presion son uni-
formes, y que el camino libre medio y la velocidad media de las moleculas son
l = 5 106 m y v = 500 m/s, respectivamente.
b) Aproximadamente, cuanto tiempo habra que esperar para que este re-
sultado fuera cierto?
1
6. Una caja de volumen 2V es dividida en mitades por una na particion. El
lado de la izquierda contiene un gas perfecto a presion p0 y el lado de la derecha
esta inicialmente vaco. Se perfora la particion con un peque
no agujero de area
A. Cu al es la presi
on p en el lado de la izquierda como funcion del tiempo?
Considere que la temperatura es constante a ambos lados de la caja, y exprese
el resultado en terminos de la velocidad media v.
2
Fsica Estadstica.
Soluciones relaci on 1
1
con m algun n
umero natural. Por lo tanto, en este caso n/m = 1 /2 sera un
n
umero racional. La condicion necesaria y suficiente para que se cubra densa-
mente el espacio de las fases (y el sistema sea ergodico) es, pues, que 1 /2 sea
un n
umero irracional; es decir, que las velocidades sean inconmensurables.
4. Estime el tiempo necesario para que una molecula de aire en una habitaci
on
llegue a una distancia de L = 5 metros de donde est a. Considere que no hay
movimiento macrosc opico del aire, que la temperatura y la presi on son uni-
formes, y que el camino libre medio y la velocidad media de las moleculas son
l = 5 106 m y v = 500 m/s, respectivamente.
l L2
t=n = = 104 s.
v lv
2
coma aleatoriamente es equivalente a seleccionar un numero n = 10000 como
supervivientes. Como N n, podemos considerar que todas las selecciones
posibles son identicas. Hay 2n formas de elegir n bacterias. Para cada se-
n
leccion, hay Cm formas de que m sean rojas. Por lo tanto, la distribucion de
probabilidad es
1 n 1 n!
p(m) = n
Cm = n , m = 0, 1, ...n.
2 2 m!(n m)!
6. Una caja de volumen 2V es dividida en mitades por una fina partici on. El
lado de la izquierda contiene un gas perfecto a presi
on p0 y el lado de la derecha
est
a inicialmente vaco. Se perfora la partici
on con un pequeno agujero de area
A. Cu al es la presi
on p en el lado de la izquierda como funci on del tiempo?
Considere que la temperatura es constante a ambos lados de la caja, y exprese
el resultado en terminos de la velocidad media v.
Como el agujero es peque no, podemos suponer que los gases a ambos lados
se encuentran en equlibrio termico en cualquier momento. Si el n umero de
partculas por unidad de volumen en el lado izquierdo en el tiempo t = 0 es n0 ,
los numeros de partculas por unidad de volumen en el tiempo t en los lados
izquierdo y derecho son n(t) y no n(t), respectivamente. Por lo tanto tenemos
que
dn(t) A A
V = n(t)v + [n0 n(t)]v,
dt 6 6
q
donde v = 8kT m es la velocidad media de las part culas. Esto se puede simpli-
ficar a
dn(t) A A
+ n(t)v = n0 v.
dt 3V 6V
Con la condicion n(0) = n0 , tenemos
n0 Av
n(t) = 1 + exp t
2 3V
y, para la presion,
p0 Av
p(t) = 1 + exp t .
2 3V
3
Fsica Estadstica. Relaci
on 2
a) el sistema y
b) la reserva.
c) A partir de a) y b), derive la relacion de Boltzmann para n1 /n2 .
1
2 .
b) Derive la aproximacion de Stirling ln n! ' n ln n n, escribiendo ln n! como
una integral.
c) Reescriba la solucion de a) utilizando el resultado de b). Encuentre el valor
de n que maximiza S(n).
d) Tratando la energa total del sistema, E, como continua, demuestre que este
sistema puede exhibir temperaturas negativas.
e) Por que es posible una temperatura negativa aqu pero no para un gas en
una caja?
2
Fsica Estadstica.
Soluciones relaci on 2
Consideremos una caja con forma de prisma de longitud L y lados con area
A, en cuyo interior se mueve una sola particula a lo largo de la direccion de la
longitud con velocidad vx . Cada vez que llega a una pared, rebota elasticamente,
luego el tiempo necesario para que choque dos veces contra la misma superficie
es
2L
t= ,
vx
es decir, que choca contra esta pared con frecuencia
vx
f= .
2L
La partcula ejercera una fuerza sobre la pared igual a la variacion de su mo-
mento por unidad de tiempo. En cada choque el momento cambia en
p = 2mvx ,
luego la fuerza es
mvx2
F = f p = .
L
La presion es la fuerza media por unidad de area, que en el caso de solo una
partcula es por tanto
mvx2
P1 = ,
V
donde V = LA es el volumen de la caja. La presion ejercida por N partculas
es por tanto
N
P = mvx2 ,
V
2 1
P 2
con vx N i vx,i la velocidad cuadr
atica media en la direccion x.
Si las partculas pueden moverse en tres dimensiones, la velocidad cuadratica
media total sera v 2 = vx2 + vy2 + vz2 . Si hay hay muchas partculas y no hay
1
ninguna predisposicion por moverse en una direccion sobre otra, tendremos
vx2 = vy2 = vz2 . Por lo tanto,
N mv 2
P = ,
3V
o, en terminos de la energa cinetica media Ec = 12 mv 2 ,
2N Ec
P = .
3V
Asociando la temperatura, T , a un promedio de la energa cinetica de las
partculas, y definiendo la constante k tal que
2Ec
= kT,
3V
tenemos
P V = N kT.
En terminos del n
umero de moles n = N/NA (con NA el n
umero de Avogadro),
y definiendo la constante R NA k, esto es
P V = nRT.
PV R
R= nT = 0.082 l atm/K/mol y k = NA = 1.38 1025 l atm/K.
d
= 0,
dt
es decir, que la funcion densidad es constante a lo largo de cualquier trayectoria
en el espacio de las fases. Desarrollando la derivada, podemos escribir:
d X X H H
= qi + qi + = + ,
dt i
qi pi t i
qi pi pi qi t
2
donde para la segunda igualdad hemos utilizado las ecuaciones de Hamilton.
Esto se puede escribir mas condensadamente utilizando corchetes de Poisson,
X H H
{, H} ,
i
qi pi pi qi
como
d
= {, H} + . (1)
dt t
Notese que esta expresion es valida para describir la evolucion temporal de
cualquier ecuacion dinamica f = f (qi , pi , t), no solo la del caso particular f = .
d
= 0,
dt
y queda demostrado el teorema.
Para un sistema en equilibrio, la densidad no puede depender explcitamente
del tiempo, luego
{, H} = 0.
Este resultado justifica el postulado fundamental de la mecanica estadstica de
equilibrio: en un sistema cerrado, a energa constante, todos los microestados
3
accesibles son igual de probables.
a) el sistema y
b) la reserva.
c) A partir de a) y b), derive la relaci
on de Boltzmann para n1 /n2 .
N! N! n2 n2
S1 = k ln k ln = k ln ' k ln
(n2 1)!(n1 + 1)! n1 !n2 ! n1 + 1 n1
E2 E1
S2 = .
T
c) De S1 + S2 = 0 tenemos
n2 E2 E1
= exp .
n1 kT
4
paralelos a la cadena. El n
umero de estados compatibles con esta condicion
es
n!
= .
m!(n m)!
Por tanto, la entropa es
n!
S = k ln = k ln x
a n ! n xa n !
a) El n
umero de microestados compatibles con n es
N!
(n) = ,
n!(N n)!
5
y la entropa
N!
S(n) = k ln(n) = k ln .
n!(N n)!
b)
n
X Z n
ln n! = ln m ' ln xdx = n ln n n + 1 ' n ln n n.
m=1 1
c)
S(n) N n
' N ln n ln .
k N n N n
La entropa maxima ocurre cuando
dS N n
=0 1 ln n + ln(N n) = 0.
dn N n N n
Por lo tanto, S = Smax cuando n = N/2.
e) La razon es que aqu la energa de una sola partcula tiene un lmite su-
perior. Para un gas, no existe este lmite y la entropa es una funcion creciente
de E, por lo que no puede darse una temperatura negativa. Desde el punto de
vista de la energa, un sistema con temperatura negativa estara mas caliente
que cualquier sistema con temperatura positiva.
6
Por lo tanto,
kT n0 RT n0
z= ln = ln .
mg n g n
Utilizando = 28 g/mol, g = 9.8m/s2 , R = 8.31J/Kmol, T = 300 K, tenemos
que n0 /n = 2 ocurre para z = 6297 m.
7
Fsica Estadstica. Relaci
on 3
November 4, 2009
XN
1 2 X
H= p + V (|ri rj |).
i=1
2m i i<j
1
c) Encuentre la capacidad calorfica a H constante.
j(j + 1)h2
j = ,
8 2 ma2
donde j = 0, 1, 2, .., y m y a son constantes positivas. La degeneracion de cada
nivel es gj = 2j + 1:
2
Fsica Estadstica.
Soluciones relaci on 3
November 4, 2009
XN
1 2 X
H= p + V (|ri rj |).
i=1
2m i i<j
QN = V N .
1
As, la funcion de particion canonica queda
3N/2
2m
Zc = V N,
de manera que
1 3 2m
ln Zc = ln + ln V.
N 2
Al igual que en la colectividad microcanonica, para que la entropa sea extensiva
en el lmite termodinamico hay que introducir la correccion de Boltzmann N !;
ademas del factor h3N . (Estos terminos pueden obtenerse como el lmite clasico
de un gas cuantico monoatomico.) Teniendo en cuenta estos factores, la funcion
de particion queda
3N/2
1 1 1 2m
Z= Zc = QN ,
N ! h3N N! h2
a) La funcion de particion es
X
Z= exp(En ),
n
donde la suma es sobre los posibles estados del sistema. Por tanto, la energa
interna es P
n En exp(En )
U = hEi = P = ln Z,
n exp(En )
donde consideramos que los {En } son constantes.
b) Para hE 2 i, tenemos
P 2
2
n En exp(En ) 1 2
hE 2 i = P = Z = ln Z + ln Z .
n exp(En ) Z 2
2
por otro,
U 2
= 2 ln Z = U 2 hE 2 i.
Luego, como mantener el conjunto {En } constante es equivalente a mantener
constante el volumen, la capacidad calorfica a volumen constante viene dada
por la formula de Einstein:
1
CV = (hE 2 i U 2 ).
kT 2
a) Cu
al es la mnima energa posible del sistema?
b) Determine la funci
on de partici
on.
c) Cu
al es la energa m
as probable del sistema?
b) La funcion de particion es
Z = (e/kT + 1 + e/kT )N .
a a1 2N sinh kT
U = N 1 = ,
a+1+ a 1 + 2 cosh kT
donde a exp(/kT ).
a) Determine la magnetizaci
on total inducida en el sistema, y obtenga aproxi-
maciones para temperaturas altas y bajas.
3
a) Definiendo a H/kT , el momento magnetico medio es
R
cos exp(a cos )d
hi = R ,
exp(a cos )d
1 a a3
coth(a) = + ...
a 3 45
y quedarnos en el termino lineal. As tenemos
N 2 H
hM i|T ' .
3kT
Por el contrario, a temperaturas bajas tenemos que coth(a) 1/a ' 1, luego
hM i|T 0 = N .
U hM i
C= = H = N k(1 x2 csch2 x2 ).
T T
c) Utilizando las expresiones del apartado a), tenemos que la susceptibilidad
magnetica a temperaturas altas es ' N 2 /(kT ), mientras que a temperat-
uras bajas, . Vemos que en el lmite de bajas temperaturas, habra
magnetizacion espontanea.
j(j + 1)h2
j = ,
8 2 ma2
donde j = 0, 1, 2, .., y m y a son constantes positivas. La degeneraci
on de cada
nivel es gj = 2j + 1:
a) Encuentre la expresi
on general para la funci
on de partici
on.
4
c) Determine la energa interna U y la capacidad calorfica CV para temper-
aturas altas.
d) Obtanga un aproximaci
on para Z a temperaturas bajas, y obtenga U y CV
para este caso.
a) Definiendo
h2
,
8 2 ma2 k
la funcion de particion sera
X
j(j + 1)
Z= (2j + 1) exp .
j=0
T
h i
b) A temperaturas altas, /T 1, y exp j(j+1)T cambia lentamente con
h i
j(j+1)
j, de modo que podemos tratar (2j + 1) exp T como una funcion con-
tinua de j. Escribiendo para este caso x j(j + 1), y dx = 2j + 1, Z adquiere
la forma integral
Z
T 8 2 ma2 kT
Z= ex/T dx = = .
0 h2
ln Z kT 2 Z
U = = = kT,
Z T
y la capacidad calorfica es CV = k.
Z ' 1 + 3e2/T ,
6k 2/T
U' e ,
Z
12k2 2/T
CV ' e .
Z 2T 2
5
Fsica Estadstica. Relaci
on 4
December 2, 2009
1. Suponga que la energa de una partcula puede ser escrita como E(z) =
az 2 , donde a es constante y z es una variable que puede representar tanto una
posicion como un momento, y que va de a . Utilizando estadstica de
Boltzamann, cual es el valor esperado de la energa?
4. Una partcula de polvo larga y fina (es decir, con forma de aguja) flota
en una caja llena de gas a temperatura constante T . En media, es el vector de
momento angular casi paralelo o casi perpendicular al eje largo de la partcula?
1
6. El circuito de la figura esta en equilibrio termico con su entorno a tem-
peratura T . Encuentre la expresion clasica para la corriente media que pasa por
el inductor.
2
Fsica Estadstica.
Soluciones relaci on 4
December 2, 2009
1. Suponga que la energa de una partcula puede ser escrita como E(z) =
az 2 , donde a es constante y z es una variable que puede representar tanto una
posici on como un momento, y que va de a . Utilizando estadstica de
Boltzamann, cu al es el valor esperado de la energa?
Seg
un la distribucion de Boltzmann, f (z) exp[E(z)/(kT )]. Por lo tanto,
el valor esperado de la energa es
R
Z az 2
exp az 2
dz
kT
E= f (z)E(z)dz = R az2 .
exp kT dz
Utilizando la identidad
Z
m cx2 m+12
x e dx = (m+1)/2 ,
0 2c
tenemos E = 12 kT .
R H
H xi x eH d
xi = R j ,
xj eH d
1
La parte integrada desaparece ya que el hamiltoniano diverge para valores
extremos. Teniendo en cuenta que xi /xj = ij , el termino restante es
Z
1
ij eH d.
Por lo tanto, el valor esperado que buscamos es
H
xi = ij kT.
xj
Para xi = xj = pi , tenemos
H
pi = hpi qi i = kT,
pi
e, igualmente, para xi = xj = qi ,
H
qi = hqi pi i = kT.
qi
Sumando sobre los 3N componentes,
* + * +
X H X
pi = pi qi = 3N kT,
i
pi i
y * + * +
X H X
qi = qi pi = 3N kT.
i
qi i
La u
ltima igualdad constituye el teorema del virial.
En muchas situaciones fsicas, el hamiltoniano del sistema es una funcion
cuadratica de sus coordenadas, luego, a traves de una transformacion canonica,
puede ser escrita en la forma
X X
H= Aj Pj2 + Bj Q2j , (1)
j j
2
cierto, se requiere que estos grados de libertad puedan ser excitados libremente.
En algunas situaciones, sin no hay suficiente energa disponible, algunos grados
de libertad pueden estar congelados y no contribuir significativamente a la
energa ni al calor especfico. En general, mientras mas alta la temperatura del
sistema, mas valido el teorema.
P
Seg
un el teorema del virial, V h i qi Fi i = 3N kT , con las fuerzas Fi =
pi . En este caso, las partculas solo interactuan con la pared, luego solo en esta
superficie existiran fuerzas. La fuerza asociada al elemento de superficie dS sera
P dS, donde P es la presion. En el lmite continuo, tenemos
I Z
V = P r S = P ( r)dV = 3P V,
S V
Digamos que el eje largo es el eje z. Los principales momentos de inercia satis-
faceran por tanto Iz < Ix , Iy . En equilibrio termico tendremos
1 1 1
Iz z2 = Ix x2 = Iy y2 ,
2 2 2
de manera que r r
Ix Iy
|z | = |x | = |y |.
Iz Iz
Por lo tanto, r
Ix
|Iz z | = |Ix x | < |Ix x |;
Iz
o, igualmente, r
Iy
|Iz z | = |Iy y | < |Iy y |.
Iz
Luego el vector momento angular es casi perpendicular al eje largo de la partcula.
3
a) La velocidad media de traslaci on de las moleculas.
b) La velocidad media de rotaci on de las moleculas alrededor de un eje perpen-
dicular al bisector de la lnea que une los dos atomos.
c) La capacidad calorfica CV .
4
Fsica Estadstica.
Prueba de clase 1
Nombre y apellido(s):
1
Fsica Estadstica.
Prueba de clase 1
U = U0 + nE.
b) Hay CnN formas de elegir n atomos de los N vertices, y CnN formas de situarlas
en N sitios entre vertices, luego = (CnN )2 . La entropa es por tanto:
N!
S = k ln = 2k ln .
n!(N n)!
N
n= .
eE/(2kT ) + 1
1
3. Utilizando la funci on de partici
on canonica, demuestre que la capacidad
calorfica a volumen constante esta relacionada con las fluctuaciones de la en-
erga, E, a traves de la f
ormula de Einstein (3 puntos):
1
CV = (hE 2 i hEi2 ).
kT 2
La funcion de particion es
X
Z= exp(En ),
n
donde la suma es sobre los posibles estados del sistema. Por tanto, la energa
interna es P
n En exp(En )
U = hEi = P = ln Z,
n exp(En )
donde consideramos que los {En } son constantes. Para hE 2 i, tenemos
P 2
2
2 n En exp(En ) 1 2
hE i = P = Z= ln Z + ln Z .
n exp(En ) Z 2
2
Fsica Estadstica. Relaci
on 5
December 9, 2009
1
Fsica Estadstica.
Soluciones relaci on 5
December 9, 2009
n
Como Vn sera proporcional a R , se puede escribir
Vn (R) = Cn Rn ,
1
Haciedo uso de la identidad
Z
1 +1
exp(y 2 )y dy = , > 1,
0 2(+1)/2 2
con = 1, tenemos que la expresion anterior vale
1 n n
n/2 = nCn = !Cn ,
2 2 2
luego la constante es
n/2 n/2
Cn = n
n
= n .
2 2 2 !
As, tenemos que el volumen de una esfera n-dimensional de radio R es
n/2
Vn (R) = n Rn ,
2 !
y su superficie,
2 n/2 n1
Sn (R) = R .
n2
2
Energa (eV) poblacion
30.1 103 3.1%
21.5 103 8.5%
12.9 103 23%
4.3 103 63%
T = 99.4K.
3
Si kT i , podemos aproximar
2
i
ci ' 4k e2i /(kT ) .
kT
4
Fsica Estadstica.
Examen parcial de febrero
Nombre y apellido(s):
U (n) es funci
on creciente de su argumento
Para r esquemas mutuamente independientes de un experimento aleatorio
se tiene U (nr ) = rU (n)
(1 punto).
3. Una cremallera tiene N enlaces. Cada enlace tiene un estado en que est a
cerrado, con energa 0, y otro en que est
a abierto, con energa . La cremallera
solo puede abrirse desde la izquierda, y un enlace solo puede abrirse si todos
los enlaces a su izquierda estan abiertos. (A veces se utiliza este modelo para
describir moleculas de ADN.)
a) Encuentre la funcion de particion (1 punto).
b) Obtenga el n umero medio de enlaces abiertos, hni, y demuestre que a tem-
peraturas bajas (T ), hni es independiente de N (1 punto).
p2x p2y p2
E(px , py , pz ) = + + z + k ,
2m 2m 2m
donde k representa los niveles de energa internos de las moleculas del gas.
a) Eval
ue la energa libre F . Muestre explcitamente la dependencia con el vol-
umen (1 punto).
Considere ahora otro reciepiente, tambien a temperatura T , con el mismo n umero
de moleculas de un gas identico y a presi
on P2 .
b) Obtenga la entropa total de los dos gases en funci on de P1 , P2 , T y N (1
punto).
c) Conectamos ahora los dos recipientes de manera que sus gases se pueden
mezclar sin hacer trabajo. Eval ue el cambio de entropa del sistema (compruebe
que su resultado tiene sentido para el caso P1 = P2 ) (1 punto).
1
Fsica Estadstica.
Soluciones examen parcial de febrero
1. Una cremallera tiene N enlaces. Cada enlace tiene un estado en que esta
cerrado, con energa 0, y otro en que esta abierto, con energa . La cremallera
solo puede abrirse desde la izquierda, y un enlace solo puede abrirse si todos
los enlaces a su izquierda estan abiertos. (A veces se utiliza este modelo para
describir moleculas de ADN.)
El n
umero medio de enlaces abiertos es
PN n
PN
n=0 ne n n
hni = = Pn=0N
Z n=0
n
1 N +1
= ln Z = ln
1
(N + 1) N +1
= .
1 1 N +1
Si kT = 1 tenemos que 1, luego
hni ' e ,
que no depende de N .
1
donde k representa los niveles de energa internos de las moleculas del gas.
a) Eval
ue la energa libre F . Muestre explcitamente la dependencia con el
volumen.
La energa libre es
!N
1 X
F = kT ln Z = kT ln eEk
N!
k
!
X
Ek
= N kT ln e + kT ln N !
k
Podemos calcular la entropa total S de ambos gases (utilizando que los gases
son identicos de manera que f (kT ) es la misma para ambos contenedores):
F
S= .
(kT ) V
La entropa de cada contenedor es
eV1
S1 = N ln N f 0 (kT ),
N
2
eV2
S2 = N ln N f 0 (kT ).,
N
y la entropa total,
eV1 eV1
S = S1 + S2 = N ln N f 0 (kT ) + N ln N f 0 (kT )
N N
kT kT
= N ln + N ln 2N f 0 (kT ) = N ln(P1 P2 ) + 2N ln(kT ) 2N f 0 (kT ).
P1 P2
c) Conectamos ahora los dos recipientes de manera que sus gases se pueden
mezclar sin hacer trabajo. Eval ue el cambio de entropa del sistema (com-
pruebe que su resultado tiene sentido para el caso P1 = P2 ).
e(V1 + V2 )
S = 2N ln 2N f 0 (kT )
2N
= 2N ln P + 2N ln(kT ) 2N f 0 (kT ).
La presion es ahora
2P1 P2
P = ,
P1 + P2
de manera que
2P1 P2
S = 2N ln + 2N ln(kT ) 2N f 0 (kT )
P1 + P2
" 2 #
2P1 P2
= N ln + 2N ln(kT ) 2N f 0 (kT ).
P1 + P2
4P1 P2 (P1 + P2 )2
S = S S = N ln = N ln .
(P1 + P2 )2 4P1 P2
(P1 + P2 )2
1 (P1 P2 )2 0,
4P1 P2
que siempre se cumple. Cuando P1 = P2 (V1 = V2 ), S = 0.
3
Fsica Estadstica.
Relaci
on 6
1
de adoptar el estado si = +1, donde es una temperatura inversa, y
1 1
P (si 1) = 1 P (si +1) = tanh (hi )
2 2
de adoptar el estado si = 1. Se puede guardar una serie de M patrones de
actividad i = 1 en los pesos sinapticos de acuerdo con la regla de Hebb,
M
1 X
ij = .
N =1 i j
Para T = 0 todos los estados con < 0 estan ocupados y los demas vacos;
mientras que para T > 0, algunos estados con > 0 estaran ocupados y asimismo
algunos con < 0 estaran vacos. Si a = b, cuanto vale ? Discuta como cambia
cuando a > b y cuando b > a..
c) Si, a T = 0, existe un exceso de nd electrones que no pueden ser acomodados
a < 0, cuanto valdra ahora ?
2
Fsica Estadstica.
Soluciones relaci on 6
S(N/2) = Smax = kN ln 2.
dS
c) De C(T ) = T dT se tiene que
Z
C(T )
Smax Smin = dT = c1 (1 ln 2).
0 T
Utilizando Smax = kN ln 2 y Smin = S(0) = 0, tenemos
ln 2
c1 = kN .
1 ln 2
1
d) La energa del sistema es E = n2 . Por tanto, la temperatura cumple
1 S S
= = ,
T E 2 n
Si K J, para n = 1 la suma es
X
eKs1 s2 = eK + eK
s1 =1
luego tenemos
Z = (eK + eK )N .
2
b) La energa libre es
F = kT ln Z = kT N ln(eK + eK ),
3
si = hsi = P (s +1) P (s 1) = 2P (s +1) 1, y obtenga:
a) Una ecuacion para el valor estacionario del solapamiento, m .
b) La temperatura crtica por encima de la cual no hay memoria. (A temper-
auras altas, puede usar la aproximacion tanh(x) = x 13 x3 + O(x5 ).)
Para T = 0 todos los estados con < 0 estan ocupados y los dem
as vacos; mien-
tras que para T > 0, algunos estados con > 0 estar an ocupados y asimismo
algunos con < 0 estar an vacos. Si a = b, cu anto vale ? Discuta c omo
cambia cuando a > b y cuando b > a..
c) Si, a T = 0, existe un exceso de nd electrones que no pueden ser acomodados
a < 0, cu
anto valdra ahora ?
a) Los electrones son fermiones (tienen espn 1/2) de manera que siguen la dis-
tribucion de Fermi, seg
un la cual la probabilidad de ocupar un nivel de energa
es
1
F () = .
1 + e()
Por lo tanto, la probabilidad de encontrar un electron con energa + es
1
F ( + ) = .
1 + e
La probabilidad de no encontrar ning
un electron a energa = es
e 1
1 F ( ) = = = F ( + ).
1 + e 1 + e
El n
umero de ausencias de electrones (huecos) a < 0 es
Z 0 Z 0 Z
1 1
np = D()[1F ()]d = b()1/2 ()
d = b(0 )1/2 (0 +)
d0 .
1 + e 0 1 + e
Ya que ne = np , tenemos,
a 1 + e(+g )
= .
b 1 + e(+)
4
Si a = b, entonces = 21 g . Si a > b, se tiene que < 1
2 2 , mientras que si
b > a, entonces > 21 2 .
c) Se tiene que Z
nd = a( g )1/2 d,
g
de donde 2/3
3nd
= g + .
2a
5
Fsica Estadstica.
Relaci
on 7
1
de un gas ideal de Fermi a T = 0 con densidad n y energa de Fermi F .
2
Fsica Estadstica.
Soluciones relaci on 7
a) La distribucion es equivalente a la de Boltzmann si e kT 1 (lmite de
gas diluido). En este caso podemos obtener el n umero de partculas como
Z
2V (2m)3/2 ()/kT
N= e d = V 3 e/kT ,
h3 0
donde = h/ 2mkT es la longitud de onda de de Broglie (longitud de onda
termica). Luego la condicion e kT 1 es equivalente a
3
V 3 d
e kT = = 1.
N
Es decir, que las partculas se encuentren separadas por distancias mucho may-
ores que las longitudes de onda asociadas a su energa cinetica.
b) La energa interna es
Z Z
2V (2m)3/2 /kT 3/2 /kT 2/kT 3/2 2/kT
U= e e d + e e d
h3 0 0
3 1 3
= N kT 1 + 3 .
2 4 2d
1
2. Sea un gas de fotones a una temperatura T en un volumen V . Sabiendo
2
que la densidad de estados es 2d
c3 :
a) Cuantos fotones hay?
b) Calcule la energa interna.
c) Que forma tiene la densidad espectral de energa, ()?
b) La energa interna es
Z 4 Z 3
V ~ 3 d ~ kT x dx
U= 2 3 = 2 3 x1
.
c 0 e~/kT 1 c ~ 0 e
R
c) De u = U/V = 0
()d, tenemos que
~ 3
() = .
2 c3 e~/kT 1
2
4. Encuentre la energa de Fermi para un gas de electrones mono-dimensional.
El n
umero de estados en un intervalo dk es
Ldk
dN = g ,
2~
donde g = 2s + 1 = 2 y L es la longitud del metal.El n
umero total de electrones
N (igual al de atomos) es
Z kF Z kF
L 2LkF
N= dN = g dk = .
kF 2~ kF ~
Por lo tanto,
N ~
pF = ~ = ,
2L 2d
donde d es la distancia entre atomos.La energa de Fermi es
kF2 (~)2
F = = .
2m 8md2
Por lo tanto,
P 5
V = P,
V T 3
o bien
1 V 3
= .
V P T 5P
A T = 0, tenemos
Z
2E 2 V ~2 k 2 ~2 kF5
P = = 2 d3 k = .
3V 3V (2)3 k<kF 2m 15m 2
3
Por otra parte, de Z
V
nV = N = 2 d3 k
(2)3 k<kF
obtenemos
kF3
n= .
3 2
En un gas ideal, la energa de una partcula es
~2 k 2
(k) = ,
2m
luego
~2 kF2
F = .
2m
As tenemos que, a T = 0,
2
p= nF ,
5
y
3
= .
2nF
4
Fsica Estadstica.
Relaci
on 8
1
Fsica Estadstica.
Soluciones relaci on 8
c) Hay magnetizaci
on expont
anea (ferromagnetismo)?
a) Cada espn (momento) puede adoptar uno de dos estados posibles, alin-
eado o antialineado, contribuyendo una energa H o H, respectivamente. La
magnetizacion, en funcion de los n
umeros de espies alineados y antialineados,
es por tanto,
1
m = (n n ).
N
n
La probabilidad de estar alineado con el campo es N eH , mientras que
n
la de estar antialineado es N eH , con (kT )1 . Por lo tanto, la
magnetizacion es
eH eH
m = H = tanh(H).
e + eH
b)
m
(T, H) = 2 .
T T cosh (H)
Si H = 0, 0 = , luego 0 = C/T , con C = k.
1
ideal.)
a) Cu
al es ahora la ecuaci
on de estado, m(T, H)?
m = tanh(H + m),
2
luego
1 = kT (1 + m2 ) .
Si H = 0 y T > Tc , tenemos que m = 0, as que
1
0 = (T Tc )1 .
k
En el caso T . Tc , podemos usar que m2 ' 3( 1), con lo que
1
' (T Tc )1 .
2k
El exponente vale, por tanto, 1.