Introducci´on a las transiciones de fase y a su simulaci´on

Introducci´on a las transiciones de fase y a su simulaci´on Introducci´on a las transiciones de fase y a su simulaci´on

sistemas.fciencias.unam.mx
from sistemas.fciencias.unam.mx More from this publisher
27.12.2013 Views

Figure9: Elprocesodenucleaciónycrecimientodegotasdeunafasemásestabledentrodeunafase metaestable en el modelo de Ising abajo de la temperatura crítica. El tiempo corre de izquierda a derecho. Varias gotas se pueden nuclear al mismo tiempo, como se percibe en las figuras. Al cabo delproceso, lanuevafasemás establedomina el sistema. 1 0.5 M 0 -0.5 -1 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 h Figure 10: Bucle de histéresis en el modelo de Ising. Se muestra la magnetización M en función del campoaplicadoh. Lasflechas indicanladirecciónenla cual h varía. hasta que se vuelva energéticamente favorable que se cree una gota de la fase nueva que ahora es más estable, es decir con la mayoría de espines para arriba. Se dice que una gota se nuclea. Despúes deesteprocesodenucleación–quepuedeinvolucrarmásdeunagota,comosemuestraenlafigura– las gotas que sehanformadocrecencon eltiempo, hasta que lanueva fasellene todo elsistema. Este proceso de nucleación siempre ocurre en las transiciones de primer orden, y del punto de vista simulacional se puede considerar como una muy fuerte indicación de la ocurrencia de una tal transicióndeprimerorden. Sinembargo,comolatransiciónnoocurredeinmediatocuandoelvalor delparámetrocruzaelvalorsupuestodelatransición, hayquetomarcuidadoconlaubicacióndela transición. Sisevuelveabajarelcampohaciavaloresnegativosdeh,entonceselsistemasequedaráenlafase magnetizada para arriba incluso para valores negativos pequeños de h. Por lo tanto, si dibujamos la magnetización en función del campo aplicado, vemos que el sistema no se regresa por la misma curva, como se demuestra en la figura 10. Este fenómeno se llama histéresis. Se observa en imanes reales,yengeneralencualquier sistema que exhibeuna transiciónde fasede primerorden. 20

(a) T = 1.3 < T c (b) T ≃ T c (c) T = 5 > T c Figure 11: Configuraciones representativas en el modelo de Ising antiferromagnético con J < 0 a distintastemperaturas T,enunaredcuadradadetamaño60 ×60: (a)abajodelatemperaturacrítica hay una fase ordenada antiferromagnéticamente, es decir alternando espines para arriba y espines para abajo; (b) cerca de la temperatura crítica se ven fluctuaciones de todos los tamaños; (c) arriba de latemperaturacrítica,enla fasedesordenada. Cuadradosrojos representanespines que apuntan paraarriba,yazulesparaabajo. 8.3 Antiferromagnetismo Lafigura11muestralatransicióndeunafasedesordenadaaunafaseantiferromagnética,variandoT con el campoexterno h = 1. Nohay ningún fenómeno denucleación, lo cualexcluye una transición deprimerorden,perosí sevencúmulos detodos los tamaños cercade T c ,ypor lotantoseconcluye tentativamentequehayunatransicióndesegundoorden. Estoseconfirma conunestudiodetallado de escalamiento de tamaño finito. Es interesante notar que todavía no se ha resuelto exactamentela termodinámica delmodelode Ising antiferromagnéticoenla presenciadeuncampo. 9 El algoritmode Wang–Landau Paraterminar elcurso, se presentaráel algoritmodeWang–Landau. Se tratade unmétodo tipoMonte Carlo desarrollado muy recientemente, que es mucho más eficiente que los métodos normales para ciertospropósitos, enparticularparacalcularlasfuncionestermodinámicasdeunsistema paratoda una gama de valores de la temperatura o de otro parámetro del sistema. Para hacerlo, el método calcula una cantidad llamada densidad de estados, que describe el sistema de una manera que no depende de la temperatura, pero de lo cual se pueden extraer las propiedades termodinámicas en funciónde latemperaturacon uncálculorelativamentesencillo. 9.1 La densidadde estados Como hemos visto, la cantidad fundamental que hay que calcular para obtener las cantidades termodinámicas de un sistema es la función de partición, Z(β). Pero la manera usual de calcularla no funciona, porque es una suma sobre demasiados estados microscópicos. Tomemos la definición de 21

(a) T = 1.3 < T c (b) T ≃ T c (c) T = 5 > T c<br />

Figure 11: Configuraciones representativas en el mo<strong>de</strong>lo <strong>de</strong> Ising antiferromagnético con J < 0 a<br />

distintastemperaturas T,enunaredcuadrada<strong>de</strong>tamaño60 ×60: (a)abajo<strong>de</strong>latemperaturacrítica<br />

hay una <strong>fase</strong> or<strong>de</strong>nada antiferromagnéticamente, es <strong>de</strong>cir alternando espines para arriba y espines<br />

para abajo; (b) cerca <strong>de</strong> la temperatura crítica se ven fluctuaciones <strong>de</strong> todos los tamaños; (c) arriba<br />

<strong>de</strong> latemperaturacrítica,enla <strong>fase</strong><strong>de</strong>sor<strong>de</strong>nada. Cuadradosrojos representanespines que apuntan<br />

paraarriba,yazulesparaabajo.<br />

8.3 Antiferromagnetismo<br />

Lafigura11muestralatransición<strong>de</strong>una<strong>fase</strong><strong>de</strong>sor<strong>de</strong>nadaauna<strong>fase</strong>antiferromagnética,variandoT<br />

con el campoexterno h = 1. Nohay ningún fenómeno <strong>de</strong>nucleación, lo cualexcluye una transición<br />

<strong>de</strong>primeror<strong>de</strong>n,perosí sevencúmulos <strong>de</strong>todos los tamaños cerca<strong>de</strong> T c ,ypor lotantoseconcluye<br />

tentativamentequehayunatransición<strong>de</strong>segundoor<strong>de</strong>n. Estoseconfirma conunestudio<strong>de</strong>tallado<br />

<strong>de</strong> escalamiento <strong>de</strong> tamaño finito. Es interesante notar que todavía no se ha re<strong>su</strong>elto exactamentela<br />

termodinámica <strong>de</strong>lmo<strong>de</strong>lo<strong>de</strong> Ising antiferromagnéticoenla presencia<strong>de</strong>uncampo.<br />

9 El algoritmo<strong>de</strong> Wang–Landau<br />

Paraterminar elcurso, se presentaráel algoritmo<strong>de</strong>Wang–Landau. Se trata<strong>de</strong> unmétodo tipoMonte<br />

Carlo <strong>de</strong>sarrollado muy recientemente, que es mucho más eficiente que los métodos normales para<br />

ciertospropósitos, enparticularparacalcular<strong>las</strong>funcionestermodinámicas<strong>de</strong>unsistema paratoda<br />

una gama <strong>de</strong> valores <strong>de</strong> la temperatura o <strong>de</strong> otro parámetro <strong>de</strong>l sistema. Para hacerlo, el método<br />

calcula una cantidad llamada <strong>de</strong>nsidad <strong>de</strong> estados, que <strong>de</strong>scribe el sistema <strong>de</strong> una manera que no<br />

<strong>de</strong>pen<strong>de</strong> <strong>de</strong> la temperatura, pero <strong>de</strong> lo cual se pue<strong>de</strong>n extraer <strong>las</strong> propieda<strong>de</strong>s termodinámicas en<br />

función<strong>de</strong> latemperaturacon uncálculorelativamentesencillo.<br />

9.1 La <strong>de</strong>nsidad<strong>de</strong> estados<br />

Como hemos visto, la cantidad fundamental que hay que calcular para obtener <strong>las</strong> cantida<strong>de</strong>s termodinámicas<br />

<strong>de</strong> un sistema es la función <strong>de</strong> partición, Z(β). Pero la manera u<strong>su</strong>al <strong>de</strong> calcularla no<br />

funciona, porque es una <strong>su</strong>ma sobre <strong>de</strong>masiados estados microscópicos. Tomemos la <strong>de</strong>finición <strong>de</strong><br />

21

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!