1 SISTEMAS DE PARTÍCULAS INTERACTIVAS -El modelo de Ising. Magnetismo. -Aproximación del campo molecular de Weiss y aproximación de Bragg-Williams. -Fonones en sólidos. –Gases clásicos no ideales.
2 Ferromagnetismo. El modelo de Ising. Material ferromagnético. Ej: red de átomos con momento magnético. Este material puede magnetizarse aplicando un campo magnético H.
A T >T* todos los momento magnéticos están al azar. A T< T* pueden ordenarse, puede haber dominios. Hamiltoniano del sistema (sin campo): Simplificación de Ising: los spines sólo tienen componente z, y el campo externo será aplicado en dicha dirección. La suma sólo se hace a primeros vecinos.
3 Hamiltoniano del sistema ferromagnético uniaxial bajo campo magnético (Hz): Definimos la magnetización por spin:
4 Colectivo canónico y termodinámica Termodinámica: Energía libre de Helmholtz
5 Función de partición y variables termodinámicas Energía del sistema = Hamiltoniano El estado microscópico del sistema es la combinación de todos los spines: Y la función de partición es: Energía libre: Magnetización promedio: Energía: Calor específico:
6 Aplicaciones del Modelo de Ising Gas de red Sitios ocupados o vacíos. Interacción a primeros vecinos, -e Ising: se introduce un spin: Nº de partículas en una celda: Nº total de partículas: Interacción: Energía total:
7 Y la gran función de partición es: Esto es el Hamiltoniano de un sistema de Ising con: (Gp:) Ising Magnetic (Gp:) Lattice Gas (Gp:) Canonical ensemble (Gp:) Grand canonical ensemble (Gp:) Coupling constant J (Gp:) Interaction energy, e (Gp:) External field H (Gp:) Chemical potential, m (Gp:) Magnetization M (Gp:) Density (Gp:) Free energy A (Gp:) pressure P (Gp:) Susceptibility ? (Gp:) Compressibility a (Gp:) Isotropic phase (Gp:) Supercritical fluid (Gp:) Ordered phase (Gp:) Liquid or gas (Gp:) Curie point (Gp:) Critical point
Se usará el colectivo macrocanónico.
8 Aleaciones binarias Interacción a primeros vecinos, : vecinos AB 0 : vecinos AA o BB Spines: Interacción: Energía total: Nº de partículas: Aplicaciones del Modelo de Ising
9 Y la gran función de partición es: Esto es el Hamiltoniano de un sistema de Ising con: (no puedo usar GC. N cte, puedo cambiar nº A y nº B) (Gp:) Ising Magnetic (Gp:) Binary Mixture (Gp:) Canonical ensemble (Gp:) Semi-grand canonical ensemble (Gp:) Coupling constant J (Gp:) Interaction energy, 2e (Gp:) External field H (Gp:) Chemical potential diff. Dm (Gp:) Magnetization M (Gp:) Composition Df (Gp:) Isotropic phase (Gp:) Mixed phase (Gp:) Ordered phase (Gp:) Separated phase (Gp:) Curie point (Gp:) Critical mixing point
Se usará el colectivo “semi-macrocanónico”.
10 Modelo de Ising en 1D (cadena lineal de spines) Usamos el colectivo canónico: Si B=0 : Funciones termodinámicas:
11 Funciones termodinámicas: Modelo de Ising en 1D
12 kT/J kT/J kT/J
13 El modelo de Ising. (en general, 2D, 3D) 2D, red cuadrada, 4 primeros vecinos g : Nº de primeros vecinos. Vamos a reescribir el Hamiltoniano: Tipo de parejas de vecinos: N++, N–, N+- Esto permite escribir: Y el Hamiltoniano queda: Orden: a largo alcance: N+ , a corto alcance: N++ acoplamiento spin-spin
14 La aproximación de Bragg-Williams Usamos el colectivo canónico: La aproximación de Bragg-Williams consiste en desarrollar un método para evaluar el peso estadístico: Se cambian las variables por comodidad:
15 La aproximación de Bragg-Williams g(L) = g(N+) , por tanto es el nº de formas de tomar N+ números de entre N números La aproximación de Bragg-Williams propone la siguiente relación: Así el Hamiltoniano resulta: Si N+ es grande, N++ será grande. El orden a corto alcance surge del orden a largo alcance. Y pasamos de necesitar hallar a buscar cómo es g(L)
16 La aproximación de Bragg-Williams Y la función de partición es: L va de –1 a 1 en pasos de 2/N Si N es grande, Stirling. Si obtenemos lnZ ya podremos tener todas las funciones termodinámicas.
17 La aproximación de Bragg-Williams Recordamos el concepto de la distribución más probable. Habrá una configuración de N+,N++ con un valor de mucho mayor que en las demás. (H ? N, por lo tanto Z depende exponencialmente de N) Por tanto buscamos el L que maximice Z o ln Z. Lmax es la única variable independiente. H, T : parámetros externos del sistema J : interacción spin-spin. N : nº de spines Si N??, lnZ tiende a ser el valor del logaritmo del mayor sumando.
18 La aproximación de Bragg-Williams ¿Cuánto vale ese Lmax ? Solución gráfica de esa igualdad: Lmax f(Lmax)
19 La aproximación de Bragg-Williams. Aplicación al ferromagnetismo Sistema sin magnetizar: H=0. La solución gráfica da Lmax(H=0)=L0 Hay solución si la pendiente de en Lmax=0 es mayor que 1. Esto corresponde a T< TC. Temperatura de Curie: La solución es: Magnetización
20 La aproximación de Bragg-Williams. Aplicación al ferromagnetismo Por tanto, con H=0, La energía libre de Helmholtz: Si nos quedamos al orden más bajo en el desarrollo para H pequeño: No se consideran cambios de TC por efecto de H.
21 Magnetización: La aproximación de Bragg-Williams. Aplicación al ferromagnetismo Calor específico: Y podéis demostrar que,
22 Solución exacta del modelo de Ising en 2D: Lars Onsager Temperatura crítica: Comportamiento en Tc, exponentes críticos: H = M 15 M = t 1/8 X = t -7/4 (log) C = t 0 (log)
23 Variaciones del modelo de Ising Modelo de Potts: diferentes valores de s, para mezclas multicomponente Modelo de Heisenberg: considera como vector. Poner redes complejas: fcc, hexagonal, etc. Tener subredes Interacción a vecinos lejanos J(r) Hacer J aleatoria, para teoría de vidrios Introducir cinética y reorientación (dependencia con t)
24 Exponentes críticos
25 Teorías de campo medio y transiciones de fase.
26 Idea General Sabemos resolver problemas de 1 partícula o de muchas partículas sin interacción (gas ideal). Ahora hay que tratar con muchas partículas interaccionando Problema: Las interacciones entre partículas hacen que sea casi imposible resolver el cálculo de la función de partición Idea: Sustituir las fuerzas que actúan sobre una partícula dada por un campo externo efectivo
Esto es un teoría de campo medio.
27 Aproximación de campo medio Se asume que el papel de las partículas vecinas es crear un campo molecular promedio, que actúa sobre la partícula estudiada. La fuerza ejercida sobre si, debido a los vecinos y al campo externo es: El campo instantáneo que actúa sobre si es: Siendo su valor promediado:
28 Para resolver el problema usaremos la mecánica estadística de momentos magnéticos sin interacción (desacoplados) Número de vecinos Y se obtiene esta ecuación para la magnetización media por espín: Solución gráfica de esa igualdad: m0 -m0 Para h=0, bJz >1 o < 1 define las dos fases:
29 Comparación entre T. Campo Medio y la solución exacta: (para red cuadrada) ¿Por qué este desacuerdo? No se consideran las fluctuaciones
30 Exponentes críticos. Cerca de la transición m