Modelo de Ising

Mecánica estadística
Termodinámica Teoría cinética
Estadísticas de partículas Teorema de estadística de espín Partículas indistinguibles Maxwell-Boltzmann Bose-Einstein Fermi-Dirac Paraestadística Estadísticas de Anyonic Estadísticas de trenzas
Conjuntos termodinámicos NVE Microcanónico NVT canónico µVT Gran canónico NPH isoentálpico-isobárico NPT isotérmico-isobárico
Modelos Debye Einstein Isando Potts
Potenciales Energía interna Entalpía Energía libre de Helmholtz Energía libre de Gibbs Gran potencial / Energía libre de Landau
Científicos Maxwell Boltzman Helmholtz Bosé Gibbs Einstein Dirac Fiesta de los testigos de Neumann Tolmán Debye Fermi Sincronizar Isando Landó
en a mi

El modelo de Ising (o modelo de Lenz-Ising ), llamado así por los físicos Ernst Ising y Wilhelm Lenz , es un modelo matemático del ferromagnetismo en mecánica estadística . El modelo consiste en variables discretas que representan momentos dipolares magnéticos de "espines" atómicos que pueden estar en uno de dos estados (+1 o −1). Los espines están dispuestos en un gráfico, generalmente una red (donde la estructura local se repite periódicamente en todas las direcciones), lo que permite que cada espín interactúe con sus vecinos. Los espines vecinos que concuerdan tienen una energía menor que los que no concuerdan; el sistema tiende a la energía más baja pero el calor perturba esta tendencia, creando así la posibilidad de diferentes fases estructurales. El modelo permite la identificación de transiciones de fase como un modelo simplificado de la realidad. El modelo de Ising de red cuadrada bidimensional es uno de los modelos estadísticos más simples para mostrar una transición de fase . ^[1]

El modelo de Ising fue inventado por el físico Wilhelm Lenz  (1920), quien lo propuso como problema a su alumno Ernst Ising. El modelo unidimensional de Ising fue resuelto por Ising (1925) en solitario en su tesis de 1924; ^[2] no tiene transición de fase. El modelo de Ising bidimensional de red cuadrada es mucho más difícil y solo recibió una descripción analítica mucho más tarde, por Lars Onsager  (1944). Generalmente se resuelve mediante un método de matriz de transferencia , aunque existen diferentes enfoques, más relacionados con la teoría cuántica de campos .

En dimensiones mayores de cuatro, la transición de fase del modelo de Ising se describe mediante la teoría del campo medio . El modelo de Ising para dimensiones mayores también se exploró con respecto a varias topologías de árboles a fines de la década de 1970, lo que culminó en una solución exacta del modelo de Barth (1981) independiente del tiempo y de campo cero para árboles de Cayley cerrados con una relación de ramificación arbitraria y, por lo tanto, una dimensionalidad arbitrariamente grande dentro de las ramas del árbol. La solución de este modelo exhibió un comportamiento de transición de fase nuevo e inusual, junto con correlaciones de espín-espín de largo alcance y vecino más cercano que no se desvanecen, consideradas relevantes para redes neuronales grandes como una de sus posibles aplicaciones.

El problema de Ising sin un campo externo se puede formular de manera equivalente como un problema de corte máximo de grafo (Max-Cut) que se puede resolver mediante optimización combinatoria .

Considere un conjunto de sitios de red, cada uno con un conjunto de sitios adyacentes (por ejemplo, un grafo ) que forman una red de dimensiones . Para cada sitio de red hay una variable discreta tal que , que representa el espín del sitio. Una configuración de espín , es una asignación de valor de espín a cada sitio de red.${\displaystyle \Lambda }$${\displaystyle d}$${\displaystyle k\in \Lambda }$${\displaystyle \sigma _{k}}$${\displaystyle \sigma _{k}\in \{-1,+1\}}$${\displaystyle {\sigma }=\{\sigma _{k}\}_{k\in \Lambda }}$

Para dos sitios adyacentes cualesquiera existe una interacción . Además, un sitio tiene un campo magnético externo que interactúa con él. La energía de una configuración está dada por la función hamiltoniana${\displaystyle i,j\in \Lambda }$ ${\displaystyle J_{ij}}$${\displaystyle j\in \Lambda }$ ${\displaystyle h_{j}}$${\displaystyle {\sigma }}$

$${\displaystyle H(\sigma )=-\sum _{\langle ij\rangle }J_{ij}\sigma _{i}\sigma _{j}-\mu \sum _{j}h_{j}\sigma _{j},}$$

donde la primera suma es sobre pares de espines adyacentes (cada par se cuenta una vez). La notación indica que los sitios y son vecinos más cercanos. El momento magnético está dado por . Nótese que el signo en el segundo término del hamiltoniano anterior debería ser en realidad positivo porque el momento magnético del electrón es antiparalelo a su espín, pero el término negativo se usa convencionalmente. ^[3] La probabilidad de configuración está dada por la distribución de Boltzmann con temperatura inversa :${\displaystyle \langle ij\rangle }$${\displaystyle i}$${\displaystyle j}$${\displaystyle \mu }$ ${\displaystyle \beta \geq 0}$

$${\displaystyle P_{\beta }(\sigma )={\frac {e^{-\beta H(\sigma )}}{Z_{\beta }}},}$$

donde , y la constante de normalización${\displaystyle \beta =1/(k_{\text{B}}T)}$

$${\displaystyle Z_{\beta }=\sum _{\sigma }e^{-\beta H(\sigma )}}$$

es la función de partición . Para una función de los espines ("observable"), se denota por${\displaystyle f}$

$${\displaystyle \langle f\rangle _{\beta }=\sum _{\sigma }f(\sigma )P_{\beta }(\sigma )}$$

el valor esperado (medio) de .${\displaystyle f}$

Las probabilidades de configuración representan la probabilidad de que (en equilibrio) el sistema esté en un estado con configuración .${\displaystyle P_{\beta }(\sigma )}$${\displaystyle \sigma }$

El signo menos en cada término de la función hamiltoniana es convencional. Utilizando esta convención de signos, los modelos de Ising pueden clasificarse según el signo de la interacción: si, para un par i ,  j${\displaystyle H(\sigma )}$

El sistema se denomina ferromagnético o antiferromagnético si todas las interacciones son ferromagnéticas o todas son antiferromagnéticas. Los modelos originales de Ising eran ferromagnéticos y todavía se suele suponer que "modelo de Ising" significa un modelo de Ising ferromagnético.

En un modelo ferromagnético de Ising, los espines tienden a estar alineados: las configuraciones en las que los espines adyacentes son del mismo signo tienen mayor probabilidad. En un modelo antiferromagnético, los espines adyacentes tienden a tener signos opuestos.

La convención de signos de H (σ) también explica cómo un sitio de espín j interactúa con el campo externo. Es decir, el sitio de espín quiere alinearse con el campo externo. Si:

Los modelos de Ising suelen examinarse sin un campo externo que interactúe con la red, es decir, h  = 0 para todos los j en la red Λ. Utilizando esta simplificación, el hamiltoniano se convierte en

$${\displaystyle H(\sigma )=-\sum _{\langle i~j\rangle }J_{ij}\sigma _{i}\sigma _{j}.}$$

Cuando el campo externo es cero en todas partes, h  = 0, el modelo de Ising es simétrico al cambiar el valor del espín en todos los sitios de la red; un campo distinto de cero rompe esta simetría.

Otra simplificación común es suponer que todos los vecinos más cercanos ⟨ ij ⟩ tienen la misma fuerza de interacción. Entonces podemos establecer J _ij = J para todos los pares i ,  j en Λ. En este caso, el hamiltoniano se simplifica aún más a

$${\displaystyle H(\sigma )=-J\sum _{\langle i~j\rangle }\sigma _{i}\sigma _{j}.}$$

Un subconjunto S del conjunto de vértices V(G) de un grafo no dirigido ponderado G determina un corte del grafo G en S y su subconjunto complementario G\S. El tamaño del corte es la suma de los pesos de las aristas entre S y G\S. Un tamaño máximo de corte es al menos el tamaño de cualquier otro corte, que varíe S.

Para el modelo de Ising sin un campo externo en un grafo G, el hamiltoniano se convierte en la siguiente suma sobre los bordes del grafo E(G)

${\displaystyle H(\sigma )=-\sum _{ij\in E(G)}J_{ij}\sigma _{i}\sigma _{j}}$.

Aquí cada vértice i del grafo es un sitio de espín que toma un valor de espín . Una configuración de espín dada divide el conjunto de vértices en dos subconjuntos dependientes de , aquellos con espín hacia arriba y aquellos con espín hacia abajo . Denotamos por el conjunto de aristas dependiente de que conecta los dos subconjuntos de vértices complementarios y . El tamaño del corte para biparticionar el grafo no dirigido ponderado G se puede definir como${\displaystyle \sigma _{i}=\pm 1}$${\displaystyle \sigma }$${\displaystyle V(G)}$${\displaystyle \sigma }$${\displaystyle V^{+}}$${\displaystyle V^{-}}$${\displaystyle \delta (V^{+})}$${\displaystyle \sigma }$${\displaystyle V^{+}}$${\displaystyle V^{-}}$ ${\displaystyle \left|\delta (V^{+})\right|}$${\displaystyle \delta (V^{+})}$

$${\displaystyle \left|\delta (V^{+})\right|={\frac {1}{2}}\sum _{ij\in \delta (V^{+})}W_{ij},}$$

donde denota un peso del borde y se introduce la escala 1/2 para compensar el conteo doble de los mismos pesos .${\displaystyle W_{ij}}$${\displaystyle ij}$${\displaystyle W_{ij}=W_{ji}}$

Las identidades

$${\displaystyle {\begin{aligned}H(\sigma )&=-\sum _{ij\in E(V^{+})}J_{ij}-\sum _{ij\in E(V^{-})}J_{ij}+\sum _{ij\in \delta (V^{+})}J_{ij}\\&=-\sum _{ij\in E(G)}J_{ij}+2\sum _{ij\in \delta (V^{+})}J_{ij},\end{aligned}}}$$

donde la suma total en el primer término no depende de , implica que minimizar en es equivalente a minimizar . Definir el peso del borde convierte así el problema de Ising sin un campo externo en un problema de corte máximo de grafo ^[4] que maximiza el tamaño del corte , que está relacionado con el hamiltoniano de Ising de la siguiente manera,${\displaystyle \sigma }$${\displaystyle H(\sigma )}$${\displaystyle \sigma }$${\displaystyle \sum _{ij\in \delta (V^{+})}J_{ij}}$${\displaystyle W_{ij}=-J_{ij}}$${\displaystyle \left|\delta (V^{+})\right|}$

$${\displaystyle H(\sigma )=\sum _{ij\in E(G)}W_{ij}-4\left|\delta (V^{+})\right|.}$$

Una cantidad significativa de preguntas estadísticas que se pueden plantear sobre este modelo se encuentran en el límite de un gran número de giros:

• En una configuración típica, ¿la mayoría de los espines son +1 o −1, o están divididos equitativamente?
• Si un giro en cualquier posición i es 1, ¿cuál es la probabilidad de que el giro en la posición j también sea 1?
• Si se cambia β , ¿hay una transición de fase?
• En una red Λ, ¿cuál es la dimensión fractal de la forma de un gran cúmulo de espines +1?

El caso más estudiado del modelo de Ising es el modelo de campo cero ferromagnético invariante a la traslación en una red de dimensión  d , es decir, Λ = Z ^d , J _ij  = 1, h  = 0.

En su tesis doctoral de 1924, Ising resolvió el modelo para el caso d  = 1, que puede considerarse como una red horizontal lineal donde cada sitio solo interactúa con su vecino izquierdo y derecho. En una dimensión, la solución no admite transición de fase . ^[5] Es decir, para cualquier β positivo, las correlaciones ⟨σ _i σ _j ⟩ decaen exponencialmente en | i  −  j |:$${\displaystyle \langle \sigma _{i}\sigma _{j}\rangle _{\beta }\leq C\exp \left(-c(\beta )|i-j|\right),}$$

y el sistema está desordenado. Sobre la base de este resultado, concluyó incorrectamente ^{[ cita requerida ]} que este modelo no exhibe comportamiento de fase en ninguna dimensión.

El modelo de Ising experimenta una transición de fase entre una fase ordenada y una desordenada en dos dimensiones o más. Es decir, el sistema está desordenado para valores β pequeños, mientras que para valores β grandes el sistema exhibe un orden ferromagnético:

$${\displaystyle \langle \sigma _{i}\sigma _{j}\rangle _{\beta }\geq c(\beta )>0.}$$

Esto fue demostrado por primera vez por Rudolf Peierls en 1936, ^[6] utilizando lo que ahora se llama un argumento de Peierls .

El modelo de Ising en una red cuadrada bidimensional sin campo magnético fue resuelto analíticamente por Lars Onsager  (1944). Onsager demostró que las funciones de correlación y la energía libre del modelo de Ising están determinadas por un fermión de red no interactuante. Onsager anunció la fórmula para la magnetización espontánea para el modelo bidimensional en 1949, pero no proporcionó una derivación. Yang (1952) dio la primera prueba publicada de esta fórmula, utilizando una fórmula límite para los determinantes de Fredholm , demostrada en 1951 por Szegő en respuesta directa al trabajo de Onsager. ^[7]

Se han derivado rigurosamente varias desigualdades de correlación para las correlaciones de espín de Ising (para estructuras reticulares generales), lo que ha permitido a los matemáticos estudiar el modelo de Ising tanto dentro como fuera de criticidad.

Dado cualquier subconjunto de espines y en la red, se cumple la siguiente desigualdad,${\displaystyle \sigma _{A}}$${\displaystyle \sigma _{B}}$

$${\displaystyle \langle \sigma _{A}\sigma _{B}\rangle \geq \langle \sigma _{A}\rangle \langle \sigma _{B}\rangle ,}$$

dónde . ${\displaystyle \langle \sigma _{A}\rangle =\langle \prod _{j\in A}\sigma _{j}\rangle }$

Con , resulta el caso especial .${\displaystyle B=\emptyset }$${\displaystyle \langle \sigma _{A}\rangle \geq 0}$

Esto significa que los espines están correlacionados positivamente en el ferroimán de Ising. Una aplicación inmediata de esto es que la magnetización de cualquier conjunto de espines aumenta con respecto a cualquier conjunto de constantes de acoplamiento .${\displaystyle \langle \sigma _{A}\rangle }$${\displaystyle J_{B}}$

La desigualdad de Simon-Lieb ^[8] establece que para cualquier conjunto que se desconecte de (por ejemplo, el límite de una caja con estar dentro de la caja y estar fuera),${\displaystyle S}$${\displaystyle x}$${\displaystyle y}$${\displaystyle x}$${\displaystyle y}$

$${\displaystyle \langle \sigma _{x}\sigma _{y}\rangle \leq \sum _{z\in S}\langle \sigma _{x}\sigma _{z}\rangle \langle \sigma _{z}\sigma _{y}\rangle .}$$

Esta desigualdad se puede utilizar para establecer la nitidez de la transición de fase para el modelo de Ising. ^[9]

Esta desigualdad se demuestra primero para un tipo de modelo de percolación correlacionado positivamente , que incluye una representación del modelo de Ising. Se utiliza para determinar las temperaturas críticas del modelo de Potts planar utilizando argumentos de percolación (que incluye el modelo de Ising como un caso especial). ^[10]

Uno de los argumentos de Demócrito en apoyo del atomismo era que los átomos explican naturalmente los límites de fase nítidos observados en los materiales ^{[ cita requerida ]} , como cuando el hielo se derrite en agua o el agua se convierte en vapor. Su idea era que pequeños cambios en las propiedades a escala atómica conducirían a grandes cambios en el comportamiento agregado. Otros creían que la materia es inherentemente continua, no atómica, y que las propiedades a gran escala de la materia no son reducibles a propiedades atómicas básicas.

Aunque las leyes de los enlaces químicos dejaron claro a los químicos del siglo XIX que los átomos eran reales, entre los físicos el debate continuó hasta bien entrado el siglo XX. Los atomistas, en particular James Clerk Maxwell y Ludwig Boltzmann , aplicaron la formulación de Hamilton de las leyes de Newton a sistemas grandes y descubrieron que el comportamiento estadístico de los átomos describe correctamente los gases a temperatura ambiente. Pero la mecánica estadística clásica no explicaba todas las propiedades de los líquidos y los sólidos, ni de los gases a baja temperatura.

Una vez formulada la mecánica cuántica moderna , el atomismo ya no estaba en conflicto con la experimentación, pero esto no condujo a una aceptación universal de la mecánica estadística, que iba más allá del atomismo. Josiah Willard Gibbs había dado un formalismo completo para reproducir las leyes de la termodinámica a partir de las leyes de la mecánica. Pero muchos argumentos erróneos sobrevivieron desde el siglo XIX, cuando la mecánica estadística se consideraba dudosa. Los lapsos de intuición se derivaban principalmente del hecho de que el límite de un sistema estadístico infinito tiene muchas leyes cero-uno que están ausentes en los sistemas finitos: un cambio infinitesimal en un parámetro puede conducir a grandes diferencias en el comportamiento general, agregado, como esperaba Demócrito.

A principios del siglo XX, algunos creían que la función de partición nunca podría describir una transición de fase, basándose en el siguiente argumento:

1. La función de partición es una suma de e ^{−β E} sobre todas las configuraciones.
2. La función exponencial es en todas partes analítica como función de β.
3. La suma de funciones analíticas es una función analítica.

Este argumento funciona para una suma finita de exponentes y establece correctamente que no hay singularidades en la energía libre de un sistema de tamaño finito. Para sistemas que están en el límite termodinámico (es decir, para sistemas infinitos), la suma infinita puede conducir a singularidades. La convergencia al límite termodinámico es rápida, de modo que el comportamiento de fase es evidente ya en una red relativamente pequeña, aunque las singularidades se suavizan por el tamaño finito del sistema.

Esto fue establecido por primera vez por Rudolf Peierls en el modelo de Ising.

Poco después de que Lenz e Ising construyeran el modelo de Ising, Peierls pudo demostrar explícitamente que una transición de fase ocurre en dos dimensiones.

Para ello, comparó los límites de alta y baja temperatura. A temperatura infinita (β = 0) todas las configuraciones tienen la misma probabilidad. Cada giro es completamente independiente de cualquier otro, y si se representan gráficamente las configuraciones típicas a temperatura infinita de modo que más/menos se representen en blanco y negro, parecen nieve de televisión . Para temperatura alta, pero no infinita, hay pequeñas correlaciones entre posiciones vecinas, la nieve tiende a agruparse un poco, pero la pantalla sigue viéndose aleatoriamente y no hay un exceso neto de negro o blanco.

Una medida cuantitativa del exceso es la magnetización , que es el valor medio del espín:

$${\displaystyle M={\frac {1}{N}}\sum _{i=1}^{N}\sigma _{i}.}$$

Un argumento falso análogo al argumento de la sección anterior establece ahora que la magnetización en el modelo de Ising es siempre cero.

1. Cada configuración de espines tiene la misma energía que la configuración con todos los espines invertidos.
2. Así que para cada configuración con magnetización M hay una configuración con magnetización − M con igual probabilidad.
3. Por lo tanto , el sistema debería pasar la misma cantidad de tiempo en la configuración con magnetización M que con magnetización − M.
4. Por lo tanto, la magnetización promedio (a lo largo del tiempo) es cero.

Como antes, esto solo demuestra que la magnetización promedio es cero en cualquier volumen finito. Para un sistema infinito, las fluctuaciones podrían no ser capaces de empujar el sistema desde un estado mayormente positivo a uno mayormente negativo con una probabilidad distinta de cero.

Para temperaturas muy altas, la magnetización es cero, como lo es a temperatura infinita. Para ver esto, note que si el espín A tiene solo una pequeña correlación ε con el espín B, y B solo está débilmente correlacionado con C, pero C es independiente de A, la cantidad de correlación de A y C es como ε ² . Para dos espines separados por una distancia L , la cantidad de correlación es como ε ^L , pero si hay más de un camino por el cual las correlaciones pueden viajar, esta cantidad aumenta por el número de caminos.

El número de caminos de longitud L en una red cuadrada en d dimensiones es ya que hay 2 d opciones para ir en cada paso.$${\displaystyle N(L)=(2d)^{L},}$$

Un límite en la correlación total está dado por la contribución a la correlación al sumar todos los caminos que unen dos puntos, que está limitada arriba por la suma de todos los caminos de longitud L dividida por que tiende a cero cuando ε es pequeño.$${\displaystyle \sum _{L}(2d)^{L}\varepsilon ^{L},}$$

A bajas temperaturas (β ≫ 1) las configuraciones están cerca de la configuración de energía más baja, aquella en la que todos los espines son positivos o negativos. Peierls se preguntó si es estadísticamente posible a baja temperatura, comenzando con todos los espines negativos, fluctuar hasta un estado en el que la mayoría de los espines sean positivos. Para que esto suceda, las gotitas de espín positivo deben poder solidificarse para formar el estado positivo.

La energía de una gota de espines positivos en un fondo negativo es proporcional al perímetro de la gota L, donde los espines positivos y negativos son vecinos entre sí. Para una gota con perímetro L , el área está en algún lugar entre ( L  − 2)/2 (la línea recta) y ( L /4) ² (la caja cuadrada). El costo de probabilidad para introducir una gota tiene el factor e ^{−β L} , pero esto contribuye a la función de partición multiplicada por el número total de gotas con perímetro L , que es menor que el número total de caminos de longitud L : De modo que la contribución total de espines de las gotas, incluso contando en exceso al permitir que cada sitio tenga una gota separada, está limitada por encima de$${\displaystyle N(L)<4^{2L}.}$$$${\displaystyle \sum _{L}L^{2}4^{2L}e^{-4\beta L},}$$

que tiende a cero en valores grandes de β. Para valores de β suficientemente grandes, esto suprime exponencialmente los bucles largos, de modo que no pueden ocurrir y la magnetización nunca fluctúa demasiado lejos de −1.

Así, Peierls estableció que la magnetización en el modelo de Ising finalmente define sectores de superselección , dominios separados no vinculados por fluctuaciones finitas.

Kramers y Wannier demostraron que la expansión a alta temperatura y la expansión a baja temperatura del modelo son iguales hasta que se produce un reajuste general de la energía libre. Esto permitió determinar con exactitud el punto de transición de fase en el modelo bidimensional (suponiendo que existe un único punto crítico).

Después de la solución de Onsager, Yang y Lee investigaron la forma en que la función de partición se vuelve singular a medida que la temperatura se acerca a la temperatura crítica.

El modelo de Ising puede resultar a menudo difícil de evaluar numéricamente si hay muchos estados en el sistema. Consideremos un modelo de Ising con

L = |Λ|: el número total de sitios en la red,

σ _j ∈ {−1, +1}: un sitio de espín individual en la red, j  = 1, ..., L ,

S ∈ {−1, +1} ^L : estado del sistema.

Dado que cada sitio de espín tiene ±1 espín, hay 2 ^L estados diferentes que son posibles. ^[11] Esto motiva la razón por la que el modelo de Ising se simula utilizando métodos de Monte Carlo . ^[11]

El hamiltoniano que se utiliza habitualmente para representar la energía del modelo cuando se utilizan métodos de Monte Carlo es Además, el hamiltoniano se simplifica aún más al suponer un campo externo h nulo , ya que muchas de las cuestiones que se plantean para ser resueltas utilizando el modelo se pueden responder en ausencia de un campo externo. Esto nos lleva a la siguiente ecuación de energía para el estado σ: Dado este hamiltoniano, se pueden calcular cantidades de interés como el calor específico o la magnetización del imán a una temperatura dada. ^[11]$${\displaystyle H(\sigma )=-J\sum _{\langle i~j\rangle }\sigma _{i}\sigma _{j}-h\sum _{j}\sigma _{j}.}$$$${\displaystyle H(\sigma )=-J\sum _{\langle i~j\rangle }\sigma _{i}\sigma _{j}.}$$

El algoritmo Metropolis-Hastings es el algoritmo de Monte Carlo más utilizado para calcular las estimaciones del modelo de Ising. ^[11] El algoritmo elige primero las probabilidades de selección g (μ, ν), que representan la probabilidad de que el algoritmo seleccione el estado ν de entre todos los estados, dado que uno está en el estado μ. Luego utiliza las probabilidades de aceptación A (μ, ν) de modo que se satisfaga el equilibrio detallado . Si se acepta el nuevo estado ν, entonces nos movemos a ese estado y repetimos con la selección de un nuevo estado y la decisión de aceptarlo. Si no se acepta ν, entonces nos quedamos en μ. Este proceso se repite hasta que se cumple algún criterio de detención, que para el modelo de Ising suele ser cuando la red se vuelve ferromagnética , lo que significa que todos los sitios apuntan en la misma dirección. ^[11]

Al implementar el algoritmo, uno debe asegurarse de que g (μ, ν) se seleccione de manera que se cumpla la ergodicidad . En equilibrio térmico, la energía de un sistema solo fluctúa dentro de un rango pequeño. ^[11] Esta es la motivación detrás del concepto de dinámica de inversión de espín único , ^[12] que establece que en cada transición, solo cambiaremos uno de los sitios de espín en la red. ^[11] Además, al usar la dinámica de inversión de espín único, uno puede pasar de cualquier estado a cualquier otro estado invirtiendo cada sitio que difiere entre los dos estados uno a la vez.

La cantidad máxima de cambio entre la energía del estado actual, H _μ y la energía de cualquier nuevo estado posible H _ν (usando la dinámica de inversión de espín único) es 2 J entre el espín que elegimos "invertir" para pasar al nuevo estado y el vecino de ese espín. ^[11] Por lo tanto, en un modelo de Ising 1D, donde cada sitio tiene dos vecinos (izquierdo y derecho), la diferencia máxima en energía sería 4 J.

Sea c el número de coordinación de red ; el número de vecinos más cercanos que tiene cualquier sitio de red. Suponemos que todos los sitios tienen el mismo número de vecinos debido a las condiciones de contorno periódicas . ^[11] Es importante señalar que el algoritmo Metropolis-Hastings no funciona bien alrededor del punto crítico debido a la desaceleración crítica. Se requieren otras técnicas como los métodos multigrid, el algoritmo de Niedermayer, el algoritmo de Swendsen-Wang o el algoritmo de Wolff para resolver el modelo cerca del punto crítico; un requisito para determinar los exponentes críticos del sistema.

Hay paquetes de código abierto disponibles que implementan estos algoritmos. ^[13]

Específicamente para el modelo de Ising y utilizando la dinámica de giro único, se puede establecer lo siguiente.

Dado que hay L sitios en total en la red, utilizando un solo giro como la única forma de realizar la transición a otro estado, podemos ver que hay un total de L estados nuevos ν a partir de nuestro estado actual μ. El algoritmo supone que las probabilidades de selección son iguales a los L estados: g (μ, ν) = 1/ L . El balance detallado nos dice que la siguiente ecuación debe cumplirse:

$${\displaystyle {\frac {P(\mu ,\nu )}{P(\nu ,\mu )}}={\frac {g(\mu ,\nu )A(\mu ,\nu )}{g(\nu ,\mu )A(\nu ,\mu )}}={\frac {A(\mu ,\nu )}{A(\nu ,\mu )}}={\frac {P_{\beta }(\nu )}{P_{\beta }(\mu )}}={\frac {{\frac {1}{Z}}e^{-\beta (H_{\nu })}}{{\frac {1}{Z}}e^{-\beta (H_{\mu })}}}=e^{-\beta (H_{\nu }-H_{\mu })}.}$$

Por lo tanto, queremos seleccionar la probabilidad de aceptación para que nuestro algoritmo satisfaga$${\displaystyle {\frac {A(\mu ,\nu )}{A(\nu ,\mu )}}=e^{-\beta (H_{\nu }-H_{\mu })}.}$$

Si H _ν > H _μ , entonces A (ν, μ) > A (μ, ν). Metropolis establece el mayor de A (μ, ν) o A (ν, μ) en 1. Según este razonamiento, el algoritmo de aceptación es: ^[11]

$${\displaystyle A(\mu ,\nu )={\begin{cases}e^{-\beta (H_{\nu }-H_{\mu })},&{\text{if }}H_{\nu }-H_{\mu }>0,\\1&{\text{otherwise}}.\end{cases}}}$$

La forma básica del algoritmo es la siguiente:

1. Elija un sitio de espín utilizando la probabilidad de selección g (μ, ν) y calcule la contribución a la energía que involucra este espín.
2. Invierte el valor del giro y calcula la nueva contribución.
3. Si la nueva energía es menor, mantenga el valor invertido.
4. Si la nueva energía es mayor, solo quédate con probabilidad${\displaystyle e^{-\beta (H_{\nu }-H_{\mu })}.}$
5. Repetir.

El cambio de energía H _ν  −  H _μ depende únicamente del valor del espín y de sus vecinos gráficos más cercanos. Por lo tanto, si el gráfico no está demasiado conectado, el algoritmo es rápido. Este proceso producirá eventualmente una selección de la distribución.

Es posible ver el modelo de Ising como una cadena de Markov , ya que la probabilidad inmediata P _β (ν) de transición a un estado futuro ν solo depende del estado presente μ. El algoritmo de Metropolis es en realidad una versión de una simulación de Monte Carlo de cadena de Markov , y dado que utilizamos dinámicas de cambio de espín único en el algoritmo de Metropolis, cada estado puede verse como si tuviera enlaces a exactamente L otros estados, donde cada transición corresponde a cambiar un solo sitio de espín al valor opuesto. ^[14] Además, dado que el cambio de la ecuación de energía H _σ solo depende de la fuerza de interacción del vecino más cercano J , el modelo de Ising y sus variantes, como el modelo de Sznajd, pueden verse como una forma de modelo de votante para la dinámica de la opinión.

El límite termodinámico existe mientras la desintegración de la interacción sea con α > 1. ^[15]${\displaystyle J_{ij}\sim |i-j|^{-\alpha }}$

• En el caso de la interacción ferromagnética con 1 < α < 2, Dyson demostró, por comparación con el caso jerárquico, que hay transición de fase a una temperatura suficientemente pequeña. ^[16]${\displaystyle J_{ij}\sim |i-j|^{-\alpha }}$
• En el caso de la interacción ferromagnética , Fröhlich y Spencer demostraron que hay una transición de fase a una temperatura suficientemente pequeña (en contraste con el caso jerárquico). ^[17]${\displaystyle J_{ij}\sim |i-j|^{-2}}$
• En el caso de interacción con α > 2 (que incluye el caso de interacciones de rango finito), no hay transición de fase a ninguna temperatura positiva (es decir, β finito), ya que la energía libre es analítica en los parámetros termodinámicos. ^[15]${\displaystyle J_{ij}\sim |i-j|^{-\alpha }}$
• En el caso de interacciones con vecinos más próximos , E. Ising proporcionó una solución exacta del modelo. A cualquier temperatura positiva (es decir, β finito), la energía libre es analítica en los parámetros termodinámicos y la correlación de espín truncada de dos puntos decae exponencialmente rápido. A temperatura cero (es decir, β infinito), hay una transición de fase de segundo orden: la energía libre es infinita y la correlación de espín truncada de dos puntos no decae (permanece constante). Por lo tanto, T = 0 es la temperatura crítica de este caso. Se satisfacen las fórmulas de escala. ^[18]

En el caso del vecino más próximo (con condiciones de contorno periódicas o libres) se dispone de una solución exacta. El hamiltoniano del modelo de Ising unidimensional en una red de sitios L con condiciones de contorno libres es donde J y h pueden ser cualquier número, ya que en este caso simplificado J es una constante que representa la fuerza de interacción entre los vecinos más próximos y h es el campo magnético externo constante aplicado a los sitios de la red. Entonces la energía libre es y la correlación espín-espín (es decir, la covarianza) es donde C (β) y c (β) son funciones positivas para T > 0. Sin embargo, para T → 0, la longitud de correlación inversa c (β) se desvanece.$${\displaystyle H(\sigma )=-J\sum _{i=1,\ldots ,L-1}\sigma _{i}\sigma _{i+1}-h\sum _{i}\sigma _{i},}$$$${\displaystyle f(\beta ,h)=-\lim _{L\to \infty }{\frac {1}{\beta L}}\ln Z(\beta )=-{\frac {1}{\beta }}\ln \left(e^{\beta J}\cosh \beta h+{\sqrt {e^{2\beta J}(\sinh \beta h)^{2}+e^{-2\beta J}}}\right),}$$$${\displaystyle \langle \sigma _{i}\sigma _{j}\rangle -\langle \sigma _{i}\rangle \langle \sigma _{j}\rangle =C(\beta )e^{-c(\beta )|i-j|},}$$

La prueba de este resultado es un cálculo simple.

Si h = 0, es muy fácil obtener la energía libre en el caso de condición de contorno libre, es decir cuando Entonces el modelo se factoriza bajo el cambio de variables.$${\displaystyle H(\sigma )=-J\left(\sigma _{1}\sigma _{2}+\cdots +\sigma _{L-1}\sigma _{L}\right).}$$$${\displaystyle \sigma '_{j}=\sigma _{j}\sigma _{j-1},\quad j\geq 2.}$$

Esto da$${\displaystyle Z(\beta )=\sum _{\sigma _{1},\ldots ,\sigma _{L}}e^{\beta J\sigma _{1}\sigma _{2}}e^{\beta J\sigma _{2}\sigma _{3}}\cdots e^{\beta J\sigma _{L-1}\sigma _{L}}=2\prod _{j=2}^{L}\sum _{\sigma '_{j}}e^{\beta J\sigma '_{j}}=2\left[e^{\beta J}+e^{-\beta J}\right]^{L-1}.}$$

Por lo tanto, la energía libre es

$${\displaystyle f(\beta ,0)=-{\frac {1}{\beta }}\ln \left[e^{\beta J}+e^{-\beta J}\right].}$$

Con el mismo cambio de variables

$${\displaystyle \langle \sigma _{j}\sigma _{j+N}\rangle =\left[{\frac {e^{\beta J}-e^{-\beta J}}{e^{\beta J}+e^{-\beta J}}}\right]^{N},}$$

Por lo tanto, decae exponencialmente tan pronto como T ≠ 0; pero para T = 0, es decir, en el límite β → ∞ no hay decaimiento.

Si h ≠ 0 necesitamos el método de la matriz de transferencia. Para las condiciones de contorno periódicas el caso es el siguiente. La función de partición es Los coeficientes pueden verse como las entradas de una matriz. Hay diferentes opciones posibles: una conveniente (porque la matriz es simétrica) es o En el formalismo matricial donde λ ₁ es el valor propio más alto de V , mientras que λ ₂ es el otro valor propio: y |λ ₂ | < λ ₁ . Esto da la fórmula de la energía libre.$${\displaystyle Z(\beta )=\sum _{\sigma _{1},\ldots ,\sigma _{L}}e^{\beta h\sigma _{1}}e^{\beta J\sigma _{1}\sigma _{2}}e^{\beta h\sigma _{2}}e^{\beta J\sigma _{2}\sigma _{3}}\cdots e^{\beta h\sigma _{L}}e^{\beta J\sigma _{L}\sigma _{1}}=\sum _{\sigma _{1},\ldots ,\sigma _{L}}V_{\sigma _{1},\sigma _{2}}V_{\sigma _{2},\sigma _{3}}\cdots V_{\sigma _{L},\sigma _{1}}.}$$${\displaystyle V_{\sigma ,\sigma '}}$$${\displaystyle V_{\sigma ,\sigma '}=e^{{\frac {\beta h}{2}}\sigma }e^{\beta J\sigma \sigma '}e^{{\frac {\beta h}{2}}\sigma '}}$$$${\displaystyle V={\begin{bmatrix}e^{\beta (h+J)}&e^{-\beta J}\\e^{-\beta J}&e^{-\beta (h-J)}\end{bmatrix}}.}$$$${\displaystyle Z(\beta )=\operatorname {Tr} \left(V^{L}\right)=\lambda _{1}^{L}+\lambda _{2}^{L}=\lambda _{1}^{L}\left[1+\left({\frac {\lambda _{2}}{\lambda _{1}}}\right)^{L}\right],}$$$${\displaystyle \lambda _{1}=e^{\beta J}\cosh \beta h+{\sqrt {e^{2\beta J}(\sinh \beta h)^{2}+e^{-2\beta J}}},}$$

La energía del estado más bajo es − JL , cuando todos los espines son iguales. Para cualquier otra configuración, la energía adicional es igual a 2 J multiplicada por el número de cambios de signo que se producen al recorrer la configuración de izquierda a derecha.

Si designamos el número de cambios de signo en una configuración como k , la diferencia de energía con respecto al estado de energía más bajo es 2 k . Dado que la energía es aditiva en el número de cambios de espín, la probabilidad p de tener un cambio de espín en cada posición es independiente. La relación entre la probabilidad de encontrar un cambio de espín y la probabilidad de no encontrarlo es el factor de Boltzmann:

$${\displaystyle {\frac {p}{1-p}}=e^{-2\beta J}.}$$

El problema se reduce a lanzamientos de moneda independientes y sesgados . Esto básicamente completa la descripción matemática.

A partir de la descripción en términos de lanzamientos independientes, se pueden entender las estadísticas del modelo para líneas largas. La línea se divide en dominios. Cada dominio tiene una longitud promedio exp(2β). La longitud de un dominio se distribuye exponencialmente, ya que existe una probabilidad constante en cualquier paso de encontrar un lanzamiento. Los dominios nunca se vuelven infinitos, por lo que un sistema largo nunca está magnetizado. Cada paso reduce la correlación entre un giro y su vecino en una cantidad proporcional a p , por lo que las correlaciones caen exponencialmente.

$${\displaystyle \langle S_{i}S_{j}\rangle \propto e^{-p|i-j|}.}$$

La función de partición es el volumen de configuraciones, cada una ponderada por su peso de Boltzmann. Como cada configuración se describe mediante los cambios de signo, la función de partición factoriza:

$${\displaystyle Z=\sum _{\text{configs}}e^{\sum _{k}S_{k}}=\prod _{k}(1+p)=(1+p)^{L}.}$$

El logaritmo dividido por L es la densidad de energía libre:

$${\displaystyle \beta f=\log(1+p)=\log \left(1+{\frac {e^{-2\beta J}}{1+e^{-2\beta J}}}\right),}$$

que es analíticamente alejada de β = ∞. Un signo de una transición de fase es una energía libre no analítica, por lo que el modelo unidimensional no tiene una transición de fase.

Para expresar el hamiltoniano de Ising utilizando una descripción mecánico cuántica de los espines, reemplazamos las variables de espín con sus respectivas matrices de Pauli . Sin embargo, dependiendo de la dirección del campo magnético, podemos crear un hamiltoniano de campo transversal o de campo longitudinal. El hamiltoniano de campo transversal viene dado por

$${\displaystyle H(\sigma )=-J\sum _{i=1,\ldots ,L}\sigma _{i}^{z}\sigma _{i+1}^{z}-h\sum _{i}\sigma _{i}^{x}.}$$

El modelo de campo transversal experimenta una transición de fase entre un régimen ordenado y desordenado en J  ~  h . Esto se puede demostrar mediante una aplicación de matrices de Pauli

$${\displaystyle \sigma _{n}^{z}=\prod _{i=1}^{n}T_{i}^{x},}$$

$${\displaystyle \sigma _{n}^{x}=T_{n}^{z}T_{n+1}^{z}.}$$

Al reescribir el hamiltoniano en términos de estas matrices de cambio de base, obtenemos

$${\displaystyle H(\sigma )=-h\sum _{i=1,\ldots ,L}T_{i}^{z}T_{i+1}^{z}-J\sum _{i}T_{i}^{x}.}$$

Dado que los roles de h y J se intercambian, el hamiltoniano experimenta una transición en J = h . ^[19]

Cuando no hay un campo externo, podemos derivar una ecuación funcional que satisfaga usando renormalización. ^[20] Específicamente, sea la función de partición con sitios. Ahora tenemos: donde . Sumamos sobre cada uno de , para obtener Ahora, como la función cosh es par, podemos resolver como . Ahora tenemos una relación de autosimilitud: Tomando el límite, obtenemos donde .${\displaystyle f(\beta ,0)=f(\beta )}$${\displaystyle Z_{N}(\beta ,J)}$${\displaystyle N}$$${\displaystyle Z_{N}(\beta ,J)=\sum _{\sigma }e^{K\sigma _{2}(\sigma _{1}+\sigma _{3})}e^{K\sigma _{4}(\sigma _{3}+\sigma _{5})}\cdots }$$${\displaystyle K:=\beta J}$${\displaystyle \sigma _{2},\sigma _{4},\cdots }$$${\displaystyle Z_{N}(\beta ,J)=\sum _{\sigma }(2\cosh(K(\sigma _{1}+\sigma _{3})))\cdot (2\cosh(K(\sigma _{3}+\sigma _{5})))\cdots }$$${\displaystyle Ae^{K'\sigma _{1}\sigma _{3}}=2\cosh(K(\sigma _{1}+\sigma _{3}))}$${\textstyle A=2{\sqrt {\cosh(2K)}},K'={\frac {1}{2}}\ln \cosh(2K)}$$${\displaystyle {\frac {1}{N}}\ln Z_{N}(K)={\frac {1}{2}}\ln \left(2{\sqrt {\cosh(2K)}}\right)+{\frac {1}{2}}{\frac {1}{N/2}}\ln Z_{N/2}(K')}$$$${\displaystyle f(\beta )={\frac {1}{2}}\ln \left(2{\sqrt {\cosh(2K)}}\right)+{\frac {1}{2}}f(\beta ')}$$${\displaystyle \beta 'J={\frac {1}{2}}\ln \cosh(2\beta J)}$

Cuando es pequeño, tenemos , por lo que podemos evaluar numéricamente iterando la ecuación funcional hasta que sea pequeño.${\displaystyle \beta }$${\displaystyle f(\beta )\approx \ln 2}$${\displaystyle f(\beta )}$${\displaystyle K}$

• En el caso ferromagnético hay una transición de fase. A baja temperatura, el argumento de Peierls prueba la magnetización positiva para el caso del vecino más próximo y luego, por la desigualdad de Griffiths , también cuando se agregan interacciones de mayor alcance. Mientras tanto, a alta temperatura, la expansión del cúmulo proporciona analiticidad de las funciones termodinámicas.
• En el caso del vecino más próximo, Onsager calculó con exactitud la energía libre mediante la equivalencia del modelo con fermiones libres en la red. McCoy y Wu calcularon las funciones de correlación espín-espín.

Onsager (1944) obtuvo la siguiente expresión analítica para la energía libre del modelo de Ising en la red cuadrada anisotrópica cuando el campo magnético en el límite termodinámico en función de la temperatura y las energías de interacción horizontal y vertical y , respectivamente${\displaystyle h=0}$${\displaystyle J_{1}}$${\displaystyle J_{2}}$

$${\displaystyle -\beta f=\ln 2+{\frac {1}{8\pi ^{2}}}\int _{0}^{2\pi }d\theta _{1}\int _{0}^{2\pi }d\theta _{2}\ln[\cosh(2\beta J_{1})\cosh(2\beta J_{2})-\sinh(2\beta J_{1})\cos(\theta _{1})-\sinh(2\beta J_{2})\cos(\theta _{2})].}$$

A partir de esta expresión para la energía libre, se pueden calcular todas las funciones termodinámicas del modelo utilizando una derivada apropiada. El modelo de Ising 2D fue el primer modelo que exhibió una transición de fase continua a una temperatura positiva. Ocurre a la temperatura que resuelve la ecuación${\displaystyle T_{c}}$

$${\displaystyle \sinh \left({\frac {2J_{1}}{kT_{c}}}\right)\sinh \left({\frac {2J_{2}}{kT_{c}}}\right)=1.}$$

En el caso isótropo cuando las energías de interacción horizontal y vertical son iguales , la temperatura crítica ocurre en el siguiente punto${\displaystyle J_{1}=J_{2}=J}$${\displaystyle T_{c}}$

$${\displaystyle T_{c}={\frac {2J}{k\ln(1+{\sqrt {2}})}}=(2.269185\cdots ){\frac {J}{k}}}$$

Cuando las energías de interacción , son ambas negativas, el modelo de Ising se convierte en un antiferromagnético. Dado que la red cuadrada es bipartita, es invariante bajo este cambio cuando el campo magnético , por lo que la energía libre y la temperatura crítica son las mismas para el caso antiferromagnético. Para la red triangular, que no es bipartita, el modelo de Ising ferromagnético y antiferromagnético se comportan de manera notablemente diferente. Específicamente, alrededor de un triángulo, es imposible hacer que los 3 pares de espines sean antiparalelos, por lo que el modelo de Ising antiferromagnético no puede alcanzar el estado de energía mínima. Este es un ejemplo de frustración geométrica .${\displaystyle J_{1}}$${\displaystyle J_{2}}$${\displaystyle h=0}$

Empecemos con una analogía con la mecánica cuántica. El modelo de Ising en una red periódica larga tiene una función de partición

$${\displaystyle \sum _{\{S\}}\exp {\biggl (}\sum _{ij}S_{i,j}\left(S_{i,j+1}+S_{i+1,j}\right){\biggr )}.}$$

Piense en la dirección i como el espacio y en la dirección j como el tiempo . Esta es una suma independiente de todos los valores que pueden tomar los espines en cada intervalo de tiempo. Este es un tipo de integral de trayectoria , es la suma de todos los historiales de espines.

Una integral de trayectoria se puede reescribir como una evolución hamiltoniana. El hamiltoniano recorre el tiempo realizando una rotación unitaria entre el tiempo t y el tiempo t + Δ t :$${\displaystyle U=e^{iH\Delta t}}$$

El producto de las matrices U, una tras otra, es el operador de evolución temporal total, que es la integral de trayectoria con la que comenzamos.

$${\displaystyle U^{N}=(e^{iH\Delta t})^{N}=\int DXe^{iL}}$$

donde N es el número de porciones de tiempo. La suma de todas las trayectorias se obtiene mediante un producto de matrices, cada elemento de la matriz es la probabilidad de transición de una porción a la siguiente.

De manera similar, se puede dividir la suma de todas las configuraciones de la función de partición en porciones, donde cada porción es la configuración unidimensional en el momento 1. Esto define la matriz de transferencia :$${\displaystyle T_{C_{1}C_{2}}.}$$

La configuración en cada porción es una colección unidimensional de espines. En cada porción temporal, T tiene elementos de matriz entre dos configuraciones de espines, una en el futuro inmediato y otra en el pasado inmediato. Estas dos configuraciones son C ₁ y C ₂ , y todas son configuraciones de espines unidimensionales. Podemos pensar en el espacio vectorial sobre el que actúa T como todas las combinaciones lineales complejas de estas. Utilizando la notación de la mecánica cuántica:$${\displaystyle |A\rangle =\sum _{S}A(S)|S\rangle }$$

donde cada vector base es una configuración de espín de un modelo de Ising unidimensional.${\displaystyle |S\rangle }$

Al igual que el hamiltoniano, la matriz de transferencia actúa sobre todas las combinaciones lineales de estados. La función de partición es una función matricial de T, que se define por la suma de todas las historias que vuelven a la configuración original después de N pasos:$${\displaystyle Z=\mathrm {tr} (T^{N}).}$$

Como se trata de una ecuación matricial, se puede evaluar en cualquier base. Por lo tanto, si podemos diagonalizar la matriz T , podemos hallar Z.

La contribución a la función de partición para cada par de configuraciones pasadas/futuras en una porción es la suma de dos términos. Existe la cantidad de cambios de espín en la porción pasada y existe la cantidad de cambios de espín entre la porción pasada y la futura. Defina un operador en las configuraciones que cambie el espín en el sitio i:

$${\displaystyle \sigma _{i}^{x}.}$$

En la base de Ising habitual, al actuar sobre cualquier combinación lineal de configuraciones pasadas, produce la misma combinación lineal pero con el espín en la posición i de cada vector base invertido.

Defina un segundo operador que multiplique el vector base por +1 y −1 según el espín en la posición i :

$${\displaystyle \sigma _{i}^{z}.}$$

T se puede escribir en términos de estos:

$${\displaystyle \sum _{i}A\sigma _{i}^{x}+B\sigma _{i}^{z}\sigma _{i+1}^{z}}$$

donde A y B son constantes que se deben determinar para reproducir la función de partición. La interpretación es que la configuración estadística en esta porción contribuye de acuerdo con el número de cambios de espín en la porción y si el espín en la posición i se ha invertido o no.

Al igual que en el caso unidimensional, trasladaremos nuestra atención de los espines a los cambios de espín. El término σ ^z en T cuenta el número de cambios de espín, que podemos escribir en términos de operadores de creación y aniquilación de cambios de espín:

$${\displaystyle \sum C\psi _{i}^{\dagger }\psi _{i}.\,}$$

El primer término cambia de dirección, por lo que, dependiendo del estado base, puede decirse:

1. Mueve un giro una unidad hacia la derecha
2. Mueve un giro una unidad hacia la izquierda
3. produce dos giros en sitios vecinos
4. Destruye dos spin-flips en sitios vecinos.

Escribiendo esto en términos de operadores de creación y aniquilación:$${\displaystyle \sigma _{i}^{x}=D{\psi ^{\dagger }}_{i}\psi _{i+1}+D^{*}{\psi ^{\dagger }}_{i}\psi _{i-1}+C\psi _{i}\psi _{i+1}+C^{*}{\psi ^{\dagger }}_{i}{\psi ^{\dagger }}_{i+1}.}$$

Ignore los coeficientes constantes y concentre su atención en la forma. Todos son cuadráticos. Como los coeficientes son constantes, esto significa que la matriz T se puede diagonalizar mediante transformadas de Fourier.

Al realizar la diagonalización se produce la energía libre de Onsager.

Onsager anunció la famosa siguiente expresión para la magnetización espontánea M de un ferroimán de Ising bidimensional en la red cuadrada en dos conferencias diferentes en 1948, aunque sin pruebas ^[7] donde y son las energías de interacción horizontal y vertical.$${\displaystyle M=\left(1-\left[\sinh 2\beta J_{1}\sinh 2\beta J_{2}\right]^{-2}\right)^{\frac {1}{8}}}$$${\displaystyle J_{1}}$${\displaystyle J_{2}}$

En 1951, Yang (1952) realizó una derivación completa utilizando un proceso de limitación de valores propios de la matriz de transferencia. La prueba fue simplificada en gran medida en 1963 por Montroll, Potts y Ward ^[7] utilizando la fórmula límite de Szegő para los determinantes de Toeplitz al tratar la magnetización como el límite de las funciones de correlación.

En el punto crítico, el modelo de Ising bidimensional es una teoría de campos conforme bidimensional . Las funciones de correlación de espín y energía se describen mediante un modelo mínimo , que se ha resuelto con exactitud.

Tanto en tres como en dos dimensiones, el caso más estudiado del modelo de Ising es el modelo invariante en la traslación sobre una red cúbica con acoplamiento de vecinos más próximos en el campo magnético cero. Muchos teóricos buscaron durante muchas décadas una solución analítica tridimensional, que fuera análoga a la solución de Onsager en el caso bidimensional. ^{[21] [22]} Tal solución no se ha encontrado hasta ahora, aunque no hay pruebas de que no exista.

En tres dimensiones, Alexander Polyakov y Vladimir Dotsenko demostraron que el modelo de Ising tiene una representación en términos de cuerdas fermiónicas que no interactúan . Esta construcción se ha llevado a cabo en la red y se desconoce el límite del continuo , que describe conjeturalmente el punto crítico.

En tres como en dos dimensiones, el argumento de Peierls muestra que hay una transición de fase. Esta transición de fase se sabe rigurosamente que es continua (en el sentido de que la longitud de correlación diverge y la magnetización tiende a cero), y se llama punto crítico . Se cree que el punto crítico puede describirse mediante un punto fijo del grupo de renormalización de la transformación del grupo de renormalización de Wilson-Kadanoff. También se cree que la transición de fase puede describirse mediante una teoría de campo conforme unitaria tridimensional, como lo evidencian las simulaciones de Monte Carlo , ^{[23] [24]} los resultados de diagonalización exactos en modelos cuánticos, ^[25] y los argumentos teóricos de campo cuántico. ^[26] Aunque es un problema abierto establecer rigurosamente la imagen del grupo de renormalización o la imagen de la teoría de campo conforme, los físicos teóricos han utilizado estos dos métodos para calcular los exponentes críticos de la transición de fase, que concuerdan con los experimentos y con las simulaciones de Monte Carlo.

Esta teoría de campo conforme que describe el punto crítico de Ising tridimensional está bajo investigación activa utilizando el método del bootstrap conforme . ^{[27] [28] [29] [30]} Este método actualmente produce la información más precisa sobre la estructura de la teoría crítica (ver exponentes críticos de Ising ).

En 2000, Sorin Istrail de Sandia National Laboratories demostró que el modelo de Ising del vidrio de espín en una red no plana es NP-completo . Es decir, suponiendo que P ≠ NP, el modelo general de Ising del vidrio de espín es exactamente solucionable solo en casos planos , por lo que las soluciones para dimensiones superiores a dos también son intratables. ^[31] El resultado de Istrail solo se refiere al modelo de vidrio de espín con acoplamientos que varían espacialmente y no dice nada sobre el modelo ferromagnético original de Ising con acoplamientos iguales.

En cualquier dimensión, el modelo de Ising se puede describir productivamente mediante un campo medio que varía localmente. El campo se define como el valor de espín promedio en una región grande, pero no tan grande como para incluir todo el sistema. El campo aún tiene variaciones lentas de un punto a otro, a medida que se mueve el volumen promedio. Estas fluctuaciones en el campo se describen mediante una teoría de campo continuo en el límite del sistema infinito.

El campo H se define como los componentes de Fourier de longitud de onda larga de la variable de espín, en el límite en que las longitudes de onda son largas. Hay muchas formas de tomar el promedio de longitud de onda larga, dependiendo de los detalles de cómo se cortan las longitudes de onda altas. Los detalles no son demasiado importantes, ya que el objetivo es encontrar las estadísticas de H y no los espines. Una vez que se conocen las correlaciones en H , las correlaciones de larga distancia entre los espines serán proporcionales a las correlaciones de larga distancia en H.

Para cualquier valor del campo de variación lenta H , la energía libre (probabilidad logarítmica) es una función analítica local de H y sus gradientes. La energía libre F ( H ) se define como la suma de todas las configuraciones de Ising que son consistentes con el campo de longitud de onda larga. Dado que H es una descripción burda, hay muchas configuraciones de Ising consistentes con cada valor de H , siempre que no se requiera demasiada exactitud para la coincidencia.

Dado que el rango permitido de valores del espín en cualquier región depende únicamente de los valores de H dentro de un volumen promedio de esa región, la contribución de energía libre de cada región depende únicamente del valor de H allí y en las regiones vecinas. Por lo tanto, F es una suma de todas las regiones de una contribución local, que depende únicamente de H y sus derivadas.

Por simetría en H , solo contribuyen las potencias pares. Por simetría de reflexión en una red cuadrada, solo contribuyen las potencias pares de gradientes. Escribiendo los primeros términos en la energía libre:

$${\displaystyle \beta F=\int d^{d}x\left[AH^{2}+\sum _{i=1}^{d}Z_{i}(\partial _{i}H)^{2}+\lambda H^{4}+\cdots \right].}$$

En una red cuadrada, las simetrías garantizan que los coeficientes Z _i de los términos derivados sean todos iguales. Pero incluso para un modelo de Ising anisotrópico, donde los Z i _en diferentes direcciones son diferentes, las fluctuaciones en H son isotrópicas en un sistema de coordenadas donde las diferentes direcciones del espacio se reescalan.

En cualquier red, el término derivado es una forma cuadrática definida positiva y se puede utilizar para definir la métrica del espacio. Por lo tanto, cualquier modelo de Ising invariante en la traslación es invariante en la rotación a grandes distancias, en coordenadas que hacen que Z _ij = δ _ij . La simetría rotacional surge espontáneamente a grandes distancias simplemente porque no hay muchos términos de orden bajo. En puntos multicríticos de orden superior, esta simetría accidental se pierde.$${\displaystyle Z_{ij}\,\partial _{i}H\,\partial _{j}H}$$

Dado que β F es una función de un campo que varía lentamente en el espacio, la probabilidad de cualquier configuración de campo es (omitiendo los términos de orden superior):

$${\displaystyle P(H)\propto e^{-\int d^{d}x\left[AH^{2}+Z|\nabla H|^{2}+\lambda H^{4}\right]}=e^{-\beta F[H]}.}$$

El promedio estadístico de cualquier producto de H términos es igual a:

$${\displaystyle \langle H(x_{1})H(x_{2})\cdots H(x_{n})\rangle ={\int DH\,e^{-\int d^{d}x\left[AH^{2}+Z|\nabla H|^{2}+\lambda H^{4}\right]}H(x_{1})H(x_{2})\cdots H(x_{n}) \over \int DH\,e^{-\int d^{d}x\left[AH^{2}+Z|\nabla H|^{2}+\lambda H^{4}\right]}}.}$$

El denominador en esta expresión se llama función de partición : y la integral sobre todos los valores posibles de H es una integral de trayectoria estadística. Integra exp(β F ) sobre todos los valores de H , sobre todos los componentes de Fourier de longitud de onda larga de los espines. F es un lagrangiano "euclidiano" para el campo H . Es similar al lagrangiano en de un campo escalar en la teoría cuántica de campos , la diferencia es que todos los términos derivados entran con un signo positivo, y no hay un factor global de i (por lo tanto "euclidiano").$${\displaystyle Z=\int DH\,e^{-\int d^{d}x\left[AH^{2}+Z|\nabla H|^{2}+\lambda H^{4}\right]}}$$

La forma de F se puede utilizar para predecir qué términos son los más importantes mediante el análisis dimensional. El análisis dimensional no es completamente sencillo, porque es necesario determinar la escala de H.

En el caso genérico, elegir la ley de escala para H es fácil, ya que el único término que contribuye es el primero,

$${\displaystyle F=\int d^{d}x\,AH^{2}.}$$

Este término es el más significativo, pero da lugar a un comportamiento trivial. Esta forma de energía libre es ultralocal, lo que significa que es una suma de una contribución independiente de cada punto. Esto es como los cambios de espín en el modelo unidimensional de Ising. Cada valor de H en cualquier punto fluctúa de forma completamente independiente del valor en cualquier otro punto.

La escala del campo puede redefinirse para absorber el coeficiente A , y entonces queda claro que A solo determina la escala general de fluctuaciones. El modelo ultralocal describe el comportamiento de alta temperatura de longitud de onda larga del modelo de Ising, ya que en este límite los promedios de fluctuación son independientes de un punto a otro.

Para encontrar el punto crítico, baje la temperatura. A medida que la temperatura baja, las fluctuaciones en H aumentan porque las fluctuaciones están más correlacionadas. Esto significa que el promedio de un gran número de espines no se vuelve pequeño tan rápidamente como si no estuvieran correlacionados, porque tienden a ser iguales. Esto corresponde a la disminución de A en el sistema de unidades donde H no absorbe A . La transición de fase solo puede ocurrir cuando los términos secundarios en F pueden contribuir, pero como el primer término domina a largas distancias, el coeficiente A debe ajustarse a cero. Esta es la ubicación del punto crítico:

$${\displaystyle F=\int d^{d}x\left[tH^{2}+\lambda H^{4}+Z(\nabla H)^{2}\right],}$$

donde t es un parámetro que pasa por cero en la transición.

Como t se desvanece, fijar la escala del campo utilizando este término hace que los otros términos se disparen. Una vez que t es pequeño, la escala del campo se puede configurar para fijar el coeficiente del término H ⁴ o el término (∇ H ) ² en 1.

Para hallar la magnetización, fijamos la escala de H de modo que λ sea uno. Ahora el campo H tiene dimensión − d /4, de modo que H ⁴ d ^d x es adimensional y Z tiene dimensión 2 −  d /2. En esta escala, el término de gradiente solo es importante a grandes distancias para d ≤ 4. Por encima de cuatro dimensiones, en longitudes de onda largas, la magnetización general solo se ve afectada por los términos ultralocales.

Hay un detalle sutil: el campo H fluctúa estadísticamente y las fluctuaciones pueden desplazar el punto cero de t . Para ver cómo, considere la división de H ⁴ de la siguiente manera:

$${\displaystyle H(x)^{4}=-\langle H(x)^{2}\rangle ^{2}+2\langle H(x)^{2}\rangle H(x)^{2}+\left(H(x)^{2}-\langle H(x)^{2}\rangle \right)^{2}}$$

El primer término es una contribución constante a la energía libre y se puede ignorar. El segundo término es un desplazamiento finito en t . El tercer término es una cantidad que escala hasta cero en distancias largas. Esto significa que al analizar el escalamiento de t mediante análisis dimensional, lo importante es el desplazamiento de t . Esto era históricamente muy confuso, porque el desplazamiento de t en cualquier λ finito es finito, pero cerca de la transición t es muy pequeño. El cambio fraccionario en t es muy grande y en unidades donde t es fijo el desplazamiento parece infinito.

La magnetización se encuentra en el mínimo de la energía libre, y ésta es una ecuación analítica. En términos del desplazamiento t ,

$${\displaystyle {\partial \over \partial H}\left(tH^{2}+\lambda H^{4}\right)=2tH+4\lambda H^{3}=0}$$

Para t < 0, los mínimos están en H proporcionales a la raíz cuadrada de t . Por lo tanto, el argumento de la catástrofe de Landau es correcto en dimensiones mayores que 5. El exponente de magnetización en dimensiones mayores que 5 es igual al valor del campo medio.

Cuando t es negativo, las fluctuaciones en torno al nuevo mínimo se describen mediante un nuevo coeficiente cuadrático positivo. Como este término siempre predomina, a temperaturas inferiores a la transición las fluctuaciones vuelven a ser ultralocales a grandes distancias.

Para encontrar el comportamiento de las fluctuaciones, reescalar el campo para fijar el término de gradiente. Entonces la dimensión de escala de longitud del campo es 1 −  d /2. Ahora el campo tiene fluctuaciones espaciales cuadráticas constantes a todas las temperaturas. La dimensión de escala del término H ² es 2, mientras que la dimensión de escala del término H ⁴ es 4 −  d . Para d < 4, el término H ⁴ tiene una dimensión de escala positiva. En dimensiones mayores que 4 tiene dimensiones de escala negativas.

Esta es una diferencia esencial. En dimensiones superiores a 4, fijar la escala del término de gradiente significa que el coeficiente del ^término H4 es cada vez menos importante a medida que las longitudes de onda son cada vez más largas. La dimensión en la que las contribuciones no cuadráticas comienzan a contribuir se conoce como dimensión crítica. En el modelo de Ising, la dimensión crítica es 4.

En dimensiones superiores a 4, las fluctuaciones críticas se describen mediante una energía libre puramente cuadrática en longitudes de onda largas. Esto significa que todas las funciones de correlación se pueden calcular a partir de promedios gaussianos :

$${\displaystyle \langle S(x)S(y)\rangle \propto \langle H(x)H(y)\rangle =G(x-y)=\int {dk \over (2\pi )^{d}}{e^{ik(x-y)} \over k^{2}+t}}$$

válido cuando x  −  y es grande. La función G ( x  −  y ) es la continuación analítica al tiempo imaginario del propagador de Feynman , ya que la energía libre es la continuación analítica de la acción del campo cuántico para un campo escalar libre. Para dimensiones 5 y superiores, todas las demás funciones de correlación a largas distancias se determinan entonces por el teorema de Wick . Todos los momentos impares son cero, por simetría ±. Los momentos pares son la suma sobre toda la partición en pares del producto de G ( x  −  y ) para cada par.

$${\displaystyle \langle S(x_{1})S(x_{2})\cdots S(x_{2n})\rangle =C^{n}\sum G(x_{i1},x_{j1})G(x_{i2},x_{j2})\ldots G(x_{in},x_{jn})}$$

donde C es la constante de proporcionalidad, por lo que basta con conocer G , que determina todas las correlaciones multipunto del campo.

Para determinar la forma de G , considere que los campos en una integral de trayectoria obedecen las ecuaciones clásicas de movimiento derivadas al variar la energía libre:

$${\displaystyle {\begin{aligned}&&\left(-\nabla _{x}^{2}+t\right)\langle H(x)H(y)\rangle &=0\\\rightarrow {}&&\nabla ^{2}G(x)+tG(x)&=0\end{aligned}}}$$

Esto es válido únicamente en puntos no coincidentes, ya que las correlaciones de H son singulares cuando los puntos chocan. H obedece a las ecuaciones clásicas de movimiento por la misma razón que los operadores mecánicos cuánticos las obedecen: sus fluctuaciones están definidas por una integral de trayectoria.

En el punto crítico t = 0, esta es la ecuación de Laplace , que se puede resolver mediante el método de Gauss de la electrostática. Defina un campo eléctrico análogo mediante

$${\displaystyle E=\nabla G}$$

Lejos del origen:

$${\displaystyle \nabla \cdot E=0}$$

ya que G es esféricamente simétrico en d dimensiones, y E es el gradiente radial de G . Integrando sobre una esfera grande de dimensión d  − 1,

$${\displaystyle \int d^{d-1}SE_{r}=\mathrm {constant} }$$

Esto da como resultado:

$${\displaystyle E={C \over r^{d-1}}}$$

y G se puede encontrar integrando con respecto a r .

$${\displaystyle G(r)={C \over r^{d-2}}}$$

La constante C fija la normalización general del campo.

Cuando t no es igual a cero, de modo que H fluctúa a una temperatura ligeramente alejada de la crítica, la función de dos puntos decae a grandes distancias. La ecuación a la que obedece se modifica:

$${\displaystyle \nabla ^{2}G+tG=0\to {1 \over r^{d-1}}{d \over dr}\left(r^{d-1}{dG \over dr}\right)+tG(r)=0}$$

Para r pequeño en comparación con , la solución diverge exactamente de la misma manera que en el caso crítico, pero se modifica el comportamiento a larga distancia.${\displaystyle {\sqrt {t}}}$

Para ver cómo, es conveniente representar la función de dos puntos como una integral, introducida por Schwinger en el contexto de la teoría cuántica de campos:

$${\displaystyle G(x)=\int d\tau {1 \over \left({\sqrt {2\pi \tau }}\right)^{d}}e^{-{x^{2} \over 4\tau }-t\tau }}$$

Esto es G , ya que la transformada de Fourier de esta integral es fácil. Cada contribución fija de τ es una gaussiana en x , cuya transformada de Fourier es otra gaussiana de ancho recíproco en k .

$${\displaystyle G(k)=\int d\tau e^{-(k^{2}-t)\tau }={1 \over k^{2}-t}}$$

Esta es la inversa del operador ∇ ²  −  t en el espacio k , que actúa sobre la función unidad en el espacio k , que es la transformada de Fourier de una fuente de función delta localizada en el origen. Por lo tanto, satisface la misma ecuación que G con las mismas condiciones de contorno que determinan la intensidad de la divergencia en 0.

La interpretación de la representación integral sobre el tiempo propio τ es que la función de dos puntos es la suma de todos los caminos de caminatas aleatorias que unen la posición 0 con la posición x a lo largo del tiempo τ. La densidad de estos caminos en el tiempo τ en la posición x es gaussiana, pero los caminantes aleatorios desaparecen a una tasa constante proporcional a t, de modo que la gaussiana en el tiempo τ disminuye en altura por un factor que disminuye constantemente de manera exponencial. En el contexto de la teoría cuántica de campos, estos son los caminos de los cuantos localizados relativísticamente en un formalismo que sigue los caminos de partículas individuales. En el contexto estadístico puro, estos caminos todavía aparecen por la correspondencia matemática con los campos cuánticos, pero su interpretación es menos directamente física.

La representación integral muestra inmediatamente que G ( r ) es positiva, ya que se representa como una suma ponderada de gaussianas positivas. También proporciona la tasa de decaimiento en r grande, ya que el tiempo adecuado para que un paseo aleatorio alcance la posición τ es r ² y en este tiempo, la altura gaussiana ha decaído en . Por lo tanto, el factor de decaimiento apropiado para la posición r es .${\displaystyle e^{-t\tau }=e^{-tr^{2}}}$${\displaystyle e^{-{\sqrt {t}}r}}$

Una aproximación heurística para G ( r ) es:

$${\displaystyle G(r)\approx {e^{-{\sqrt {t}}r} \over r^{d-2}}}$$

Esta no es una forma exacta, excepto en tres dimensiones, donde las interacciones entre trayectorias cobran importancia. Las formas exactas en dimensiones superiores son variantes de las funciones de Bessel .

La interpretación de las correlaciones como cuantos de tamaño fijo que recorren caminos aleatorios permite entender por qué la dimensión crítica de la interacción H ^{4 es 4. El término} H ⁴ puede considerarse como el cuadrado de la densidad de los caminantes aleatorios en cualquier punto. Para que un término de este tipo altere las funciones de correlación de orden finito, que solo introducen unos pocos nuevos caminos aleatorios en el entorno fluctuante, los nuevos caminos deben intersecarse. De lo contrario, el cuadrado de la densidad es simplemente proporcional a la densidad y solo desplaza el coeficiente H ² en una constante. Pero la probabilidad de intersección de los caminos aleatorios depende de la dimensión, y los caminos aleatorios de dimensión superior a 4 no se intersecan.

La dimensión fractal de un paseo aleatorio ordinario es 2. El número de bolas de tamaño ε necesarias para cubrir el camino aumenta a medida que ε ⁻² . Dos objetos de dimensión fractal 2 se intersectarán con una probabilidad razonable solo en un espacio de dimensión 4 o menor, la misma condición que para un par genérico de planos. Kurt Symanzik argumentó que esto implica que las fluctuaciones críticas de Ising en dimensiones mayores que 4 deberían describirse mediante un campo libre. Este argumento finalmente se convirtió en una prueba matemática.

El modelo de Ising en cuatro dimensiones se describe mediante un campo fluctuante, pero ahora las fluctuaciones interactúan. En la representación del polímero, las intersecciones de recorridos aleatorios son marginalmente posibles. En la continuación del campo cuántico, los cuantos interactúan.

El logaritmo negativo de la probabilidad de cualquier configuración de campo H es la función de energía libre

$${\displaystyle F=\int d^{4}x\left[{Z \over 2}|\nabla H|^{2}+{t \over 2}H^{2}+{\lambda \over 4!}H^{4}\right]\,}$$

Los factores numéricos están ahí para simplificar las ecuaciones de movimiento. El objetivo es comprender las fluctuaciones estadísticas. Como cualquier otra integral de trayectoria no cuadrática, las funciones de correlación tienen una expansión de Feynman como partículas que viajan a lo largo de caminos aleatorios, dividiéndose y volviéndose a unir en los vértices. La fuerza de interacción está parametrizada por la cantidad clásicamente adimensional λ.

Aunque el análisis dimensional muestra que tanto λ como Z son adimensionales, esto es engañoso. Las fluctuaciones estadísticas de longitud de onda larga no son exactamente invariantes de escala y solo se vuelven invariantes de escala cuando la fuerza de interacción desaparece.

La razón es que se utiliza un valor de corte para definir H y este valor de corte define la longitud de onda más corta. Las fluctuaciones de H en longitudes de onda cercanas al valor de corte pueden afectar las fluctuaciones de longitud de onda más larga. Si el sistema se escala junto con el valor de corte, los parámetros se escalarán mediante análisis dimensional, pero luego la comparación de parámetros no compara el comportamiento porque el sistema reescalado tiene más modos. Si el sistema se reescala de tal manera que el valor de corte de longitud de onda corta permanece fijo, las fluctuaciones de longitud de onda larga se modifican.

Una forma heurística rápida de estudiar el escalamiento es cortar los números de onda H en un punto λ. Los modos de Fourier de H con números de onda mayores que λ no pueden fluctuar. Un reescalamiento de la longitud que haga que todo el sistema sea más pequeño aumenta todos los números de onda y desplaza algunas fluctuaciones por encima del límite.

Para restablecer el antiguo límite, realice una integración parcial sobre todos los números de onda que solían estar prohibidos, pero que ahora fluctúan. En los diagramas de Feynman, la integración sobre un modo fluctuante en el número de onda k vincula las líneas que llevan el momento k en una función de correlación en pares, con un factor del propagador inverso.

En el reescalado, cuando el sistema se reduce por un factor de (1+ b ), el coeficiente t aumenta por un factor de (1+ b ) ² según el análisis dimensional. El cambio en t para un valor infinitesimal b es 2 bt . Los otros dos coeficientes son adimensionales y no cambian en absoluto.

El efecto de orden más bajo de la integración se puede calcular a partir de las ecuaciones de movimiento:

$${\displaystyle \nabla ^{2}H+tH=-{\lambda \over 6}H^{3}.}$$

Esta ecuación es una identidad dentro de cualquier función de correlación, alejada de otras inserciones. Después de integrar los modos con Λ < k < (1+ b )Λ, será una identidad ligeramente diferente.

Como la forma de la ecuación se conservará, para encontrar el cambio en los coeficientes es suficiente analizar el cambio en el término H ^{3 . En una expansión del diagrama de Feynman, el término} H ³ en una función de correlación dentro de una correlación tiene tres líneas colgantes. Uniendo dos de ellas en un número de onda grande k se obtiene un cambio H ³ con una línea colgante, por lo que es proporcional a H :

$${\displaystyle \delta H^{3}=3H\int _{\Lambda <|k|<(1+b)\Lambda }{d^{4}k \over (2\pi )^{4}}{1 \over (k^{2}+t)}}$$

El factor 3 proviene del hecho de que el bucle se puede cerrar de tres maneras diferentes.

La integral debe dividirse en dos partes:

$${\displaystyle \int dk{1 \over k^{2}}-t\int dk{1 \over k^{2}(k^{2}+t)}=A\Lambda ^{2}b+Bbt}$$

La primera parte no es proporcional a t y, en la ecuación de movimiento, puede ser absorbida por un desplazamiento constante en t . Esto se debe a que el término H ³ tiene una parte lineal. Solo el segundo término, que varía de t a t , contribuye al escalamiento crítico.

Este nuevo término lineal se suma al primer término del lado izquierdo, lo que cambia t en una cantidad proporcional a t . El cambio total en t es la suma del término del análisis dimensional y este segundo término de los productos de operadores :

$${\displaystyle \delta t=\left(2-{B\lambda \over 2}\right)bt}$$

Entonces t se reescala, pero su dimensión es anómala , se modifica en una cantidad proporcional al valor de λ.

Pero λ también cambia. El cambio en λ requiere considerar que las líneas se dividen y luego se vuelven a unir rápidamente. El proceso de orden más bajo es aquel en el que una de las tres líneas de H ³ se divide en tres, que rápidamente se unen con una de las otras líneas del mismo vértice. La corrección del vértice es

$${\displaystyle \delta \lambda =-{3\lambda ^{2} \over 2}\int _{k}dk{1 \over (k^{2}+t)^{2}}=-{3\lambda ^{2} \over 2}b}$$

El factor numérico es tres veces mayor porque hay un factor adicional de tres al elegir cuál de las tres nuevas líneas contratar. Por lo tanto

$${\displaystyle \delta \lambda =-3B\lambda ^{2}b}$$

Estas dos ecuaciones juntas definen las ecuaciones del grupo de renormalización en cuatro dimensiones:

$${\displaystyle {\begin{aligned}{dt \over t}&=\left(2-{B\lambda \over 2}\right)b\\{d\lambda \over \lambda }&={-3B\lambda \over 2}b\end{aligned}}}$$

El coeficiente B se determina mediante la fórmula$${\displaystyle Bb=\int _{\Lambda <|k|<(1+b)\Lambda }{d^{4}k \over (2\pi )^{4}}{1 \over k^{4}}}$$

y es proporcional al área de una esfera tridimensional de radio λ, multiplicada por el ancho de la región de integración b Λ dividido por Λ ⁴ :$${\displaystyle B=(2\pi ^{2}\Lambda ^{3}){1 \over (2\pi )^{4}}{b\Lambda }{1 \over b\Lambda ^{4}}={1 \over 8\pi ^{2}}}$$

En otras dimensiones, la constante B cambia, pero la misma constante aparece tanto en el flujo t como en el flujo de acoplamiento. La razón es que la derivada con respecto a t del bucle cerrado con un solo vértice es un bucle cerrado con dos vértices. Esto significa que la única diferencia entre el escalado del acoplamiento y el t son los factores combinatorios de unión y división.

Debería ser posible investigar tres dimensiones partiendo de la teoría de cuatro dimensiones, porque las probabilidades de intersección de los recorridos aleatorios dependen continuamente de la dimensionalidad del espacio. En el lenguaje de los grafos de Feynman, el acoplamiento no cambia mucho cuando se cambia la dimensión.

El proceso de continuar alejándose de la dimensión 4 no está completamente bien definido sin una prescripción sobre cómo hacerlo. La prescripción solo está bien definida en los diagramas. Reemplaza la representación de Schwinger en la dimensión 4 con la representación de Schwinger en la dimensión 4 − ε definida por:$${\displaystyle G(x-y)=\int d\tau {1 \over t^{d \over 2}}e^{{x^{2} \over 2\tau }+t\tau }}$$

En la dimensión 4 − ε, el acoplamiento λ tiene una dimensión de escala positiva ε, y esta debe agregarse al flujo.

$${\displaystyle {\begin{aligned}{d\lambda \over \lambda }&=\varepsilon -3B\lambda \\{dt \over t}&=2-\lambda B\end{aligned}}}$$

El coeficiente B depende de la dimensión, pero se cancelará. El punto fijo para λ ya no es cero, sino: donde las dimensiones de escala de t se modifican en una cantidad λ B = ε/3.$${\displaystyle \lambda ={\varepsilon \over 3B}}$$

El exponente de magnetización se altera proporcionalmente a:$${\displaystyle {\tfrac {1}{2}}\left(1-{\varepsilon \over 3}\right)}$$

que es .333 en 3 dimensiones (ε = 1) y .166 en 2 dimensiones (ε = 2). Esto no está tan lejos del exponente medido .308 y el exponente bidimensional de Onsager .125.

El comportamiento de un modelo de Ising en un grafo completamente conectado puede entenderse completamente mediante la teoría del campo medio . Este tipo de descripción es apropiada para redes cuadradas de dimensiones muy altas, porque entonces cada sitio tiene una cantidad muy grande de vecinos.

La idea es que si cada espín está conectado a un gran número de espines, solo es importante la relación media de espines + a espines −, ya que las fluctuaciones sobre esta media serán pequeñas. El campo medio H es la fracción media de espines que son + menos la fracción media de espines que son −. El coste energético de invertir un único espín en el campo medio H es ±2 JNH . Es conveniente redefinir J para absorber el factor N , de modo que el límite N → ∞ sea suave. En términos del nuevo J , el coste energético de invertir un espín es ±2 JH .

Este coste de energía da la relación entre la probabilidad p de que el espín sea + y la probabilidad 1− p de que el espín sea −. Esta relación es el factor de Boltzmann:$${\displaystyle {p \over 1-p}=e^{2\beta JH}}$$

de modo que$${\displaystyle p={1 \over 1+e^{-2\beta JH}}}$$

El valor medio del espín se da promediando 1 y −1 con los pesos p y 1 −  p , por lo que el valor medio es 2 p  − 1. Pero este promedio es el mismo para todos los espines y, por lo tanto, es igual a H .$${\displaystyle H=2p-1={1-e^{-2\beta JH} \over 1+e^{-2\beta JH}}=\tanh(\beta JH)}$$

Las soluciones de esta ecuación son los posibles campos de medias consistentes. Para β J < 1 solo hay una solución en H = 0. Para valores mayores de β hay tres soluciones y la solución en H = 0 es inestable.

La inestabilidad significa que aumentar un poco el campo medio por encima de cero produce una fracción estadística de espines que son + que es mayor que el valor del campo medio. Por lo tanto, un campo medio que fluctúa por encima de cero producirá un campo medio aún mayor y, finalmente, se asentará en la solución estable. Esto significa que, para temperaturas por debajo del valor crítico β J = 1, el modelo de Ising de campo medio experimenta una transición de fase en el límite de N grande .

Por encima de la temperatura crítica, las fluctuaciones de H se amortiguan porque el campo medio restablece la fluctuación al campo cero. Por debajo de la temperatura crítica, el campo medio se lleva a un nuevo valor de equilibrio, que es la solución H positiva o H negativa de la ecuación.

Para β J = 1 + ε, justo por debajo de la temperatura crítica, el valor de H se puede calcular a partir de la expansión de Taylor de la tangente hiperbólica:$${\displaystyle H=\tanh(\beta JH)\approx (1+\varepsilon )H-{(1+\varepsilon )^{3}H^{3} \over 3}}$$

Dividiendo por H para descartar la solución inestable en H = 0, las soluciones estables son:$${\displaystyle H={\sqrt {3\varepsilon }}}$$

La magnetización espontánea H crece cerca del punto crítico como la raíz cuadrada del cambio de temperatura. Esto es cierto siempre que H pueda calcularse a partir de la solución de una ecuación analítica que sea simétrica entre valores positivos y negativos, lo que llevó a Landau a sospechar que todas las transiciones de fase de tipo Ising en todas las dimensiones deberían seguir esta ley.

El exponente del campo medio es universal porque los cambios en el carácter de las soluciones de ecuaciones analíticas siempre se describen mediante catástrofes en la serie de Taylor , que es una ecuación polinómica. Por simetría, la ecuación para H solo debe tener potencias impares de H en el lado derecho. Cambiar β solo debería cambiar suavemente los coeficientes. La transición ocurre cuando el coeficiente de H en el lado derecho es 1. Cerca de la transición:$${\displaystyle H={\partial (\beta F) \over \partial h}=(1+A\varepsilon )H+BH^{3}+\cdots }$$

Sean cuales sean A y B , siempre que ninguno de ellos esté ajustado a cero, la magnetización espontánea crecerá como la raíz cuadrada de ε. Este argumento sólo puede fallar si la energía libre β F no es analítica o no genérica en el β exacto donde ocurre la transición.

Sin embargo, la magnetización espontánea en sistemas magnéticos y la densidad en gases cerca del punto crítico se miden con mucha precisión. La densidad y la magnetización en tres dimensiones tienen la misma dependencia de ley de potencia con respecto a la temperatura cerca del punto crítico, pero el comportamiento de los experimentos es:$${\displaystyle H\propto \varepsilon ^{0.308}}$$

El exponente también es universal, ya que es el mismo en el modelo de Ising que en el imán y el gas experimentales, pero no es igual al valor del campo medio. Esto fue una gran sorpresa.

Esto también es cierto en dos dimensiones, donde$${\displaystyle H\propto \varepsilon ^{0.125}}$$

Pero allí no hubo ninguna sorpresa, porque fue predicho por Onsager .

En tres dimensiones, la serie perturbativa de la teoría de campos es una expansión en una constante de acoplamiento λ que no es particularmente pequeña. El tamaño efectivo del acoplamiento en el punto fijo es uno sobre el factor de ramificación de las trayectorias de las partículas, por lo que el parámetro de expansión es aproximadamente 1/3. En dos dimensiones, el parámetro de expansión perturbativa es 2/3.

Pero la renormalización también se puede aplicar de forma productiva a los espines directamente, sin pasar a un campo promedio. Históricamente, este enfoque se debe a Leo Kadanoff y es anterior a la expansión perturbativa ε.

La idea es integrar iterativamente los espines reticulares, generando un flujo en los acoplamientos. Pero ahora los acoplamientos son coeficientes de energía reticular. El hecho de que exista una descripción continua garantiza que esta iteración convergerá a un punto fijo cuando la temperatura se ajuste a la criticidad.

Escriba el modelo de Ising bidimensional con un número infinito de posibles interacciones de orden superior. Para mantener la simetría de reflexión de espín, solo contribuyen las potencias pares:$${\displaystyle E=\sum _{ij}J_{ij}S_{i}S_{j}+\sum J_{ijkl}S_{i}S_{j}S_{k}S_{l}\ldots .}$$

Por invariancia de traslación, J _ij es solo una función de ij. Por la simetría rotacional accidental, en i y j grandes su tamaño solo depende de la magnitud del vector bidimensional i  −  j . Los coeficientes de orden superior también están restringidos de manera similar.

La iteración de renormalización divide la red en dos partes: espines pares y espines impares. Los espines impares se encuentran en las posiciones de la red del tablero de ajedrez impar, y los pares en las posiciones del tablero de ajedrez par. Cuando los espines se indexan por la posición ( i , j ), los sitios impares son aquellos con i  +  j impar y los sitios pares aquellos con i  +  j par, y los sitios pares solo están conectados a los sitios impares.

Los dos valores posibles de los espines impares se integrarán mediante la suma de ambos valores posibles. Esto producirá una nueva función de energía libre para los espines pares restantes, con nuevos acoplamientos ajustados. Los espines pares vuelven a estar en una red, con ejes inclinados 45 grados respecto a los antiguos. Al desrotar el sistema se restablece la configuración anterior, pero con nuevos parámetros. Estos parámetros describen la interacción entre espines a distancias mayores.${\textstyle {\sqrt {2}}}$

Partiendo del modelo de Ising y repitiendo esta iteración, eventualmente se modifican todos los acoplamientos. Cuando la temperatura es más alta que la temperatura crítica, los acoplamientos convergerán a cero, ya que los espines a grandes distancias no están correlacionados. Pero cuando la temperatura es crítica, habrá coeficientes distintos de cero que vincularán espines en todos los órdenes. El flujo se puede aproximar considerando solo los primeros términos. Este flujo truncado producirá aproximaciones cada vez mejores a los exponentes críticos cuando se incluyan más términos.

La aproximación más simple es mantener únicamente el término J habitual y descartar todo lo demás. Esto generará un flujo en J , análogo al flujo en t en el punto fijo de λ en la expansión ε.

Para hallar el cambio en J , considere los cuatro vecinos de un sitio impar. Estos son los únicos espines que interactúan con él. La contribución multiplicativa a la función de partición de la suma de los dos valores del espín en el sitio impar es:$${\displaystyle e^{J(N_{+}-N_{-})}+e^{J(N_{-}-N_{+})}=2\cosh(J[N_{+}-N_{-}])}$$

donde N _± es el número de vecinos que son ±. Ignorando el factor 2, la contribución de energía libre de este sitio extraño es:$${\displaystyle F=\log(\cosh[J(N_{+}-N_{-})]).}$$

Esto incluye interacciones entre el vecino más cercano y el vecino más cercano siguiente, como se esperaba, pero también una interacción de cuatro espines que debe descartarse. Para truncar las interacciones entre el vecino más cercano, considere que la diferencia de energía entre todos los espines iguales y con números iguales de + y – es:$${\displaystyle \Delta F=\ln(\cosh[4J]).}$$

A partir de los acoplamientos de vecinos más próximos, la diferencia de energía entre todos los espines iguales y espines escalonados es 8 J . La diferencia de energía entre todos los espines iguales y no escalonados pero con espín neto cero es 4 J . Ignorando las interacciones de cuatro espines, un truncamiento razonable es el promedio de estas dos energías o 6 J . Dado que cada enlace contribuirá a dos espines impares, el valor correcto para comparar con el anterior es la mitad:$${\displaystyle 3J'=\ln(\cosh[4J]).}$$

Para J pequeñas , esto fluye rápidamente a un acoplamiento cero. Para J grandes , fluye a acoplamientos grandes. El exponente de magnetización se determina a partir de la pendiente de la ecuación en el punto fijo.

Las variantes de este método producen buenas aproximaciones numéricas para los exponentes críticos cuando se incluyen muchos términos, tanto en dos como en tres dimensiones.

La motivación original del modelo fue el fenómeno del ferromagnetismo . El hierro es magnético; una vez magnetizado, permanece magnetizado durante mucho tiempo en comparación con cualquier tiempo atómico.

En el siglo XIX se pensaba que los campos magnéticos se debían a corrientes en la materia, y Ampère postuló que los imanes permanentes eran causados ​​por corrientes atómicas permanentes. Sin embargo, el movimiento de partículas cargadas clásicas no podía explicar las corrientes permanentes, como demostró Larmor . Para que haya ferromagnetismo, los átomos deben tener momentos magnéticos permanentes que no se deban al movimiento de cargas clásicas.

Una vez descubierto el espín del electrón, quedó claro que el magnetismo debía deberse a un gran número de espines de electrones, todos ellos orientados en la misma dirección. Era natural preguntarse cómo los espines de los electrones sabían en qué dirección apuntar, porque los electrones de un lado de un imán no interactúan directamente con los electrones del otro lado. Sólo pueden influir en sus vecinos. El modelo de Ising fue diseñado para investigar si una gran fracción de los espines de los electrones podían orientarse en la misma dirección utilizando sólo fuerzas locales.

El modelo de Ising puede reinterpretarse como un modelo estadístico del movimiento de los átomos. Dado que la energía cinética depende únicamente del momento y no de la posición, mientras que la estadística de las posiciones depende únicamente de la energía potencial, la termodinámica del gas depende únicamente de la energía potencial para cada configuración de átomos.

Un modelo burdo consiste en hacer del espacio-tiempo una red e imaginar que cada posición contiene un átomo o no. El espacio de configuración es el de los bits independientes B _i , donde cada bit es 0 o 1 dependiendo de si la posición está ocupada o no. Una interacción atractiva reduce la energía de dos átomos cercanos. Si la atracción es solo entre vecinos más cercanos, la energía se reduce en −4 JB _i B _j por cada par vecino ocupado.

La densidad de los átomos se puede controlar añadiendo un potencial químico , que es un coste de probabilidad multiplicativo por añadir un átomo más. Un factor multiplicativo en la probabilidad se puede reinterpretar como un término aditivo en el logaritmo: la energía. La energía extra de una configuración con N átomos se cambia por μN . El coste de probabilidad de un átomo más es un factor de exp(− βμ ).

Entonces la energía del gas reticular es:$${\displaystyle E=-{\frac {1}{2}}\sum _{\langle i,j\rangle }4JB_{i}B_{j}+\sum _{i}\mu B_{i}}$$

Reescribiendo los bits en términos de espines,${\displaystyle B_{i}=(S_{i}+1)/2.}$$${\displaystyle E=-{\frac {1}{2}}\sum _{\langle i,j\rangle }JS_{i}S_{j}-{\frac {1}{2}}\sum _{i}(4J-\mu )S_{i}}$$

Para redes donde cada sitio tiene un número igual de vecinos, este es el modelo de Ising con un campo magnético h = ( zJ  −  μ )/2, donde z es el número de vecinos.

En los sistemas biológicos, se han utilizado versiones modificadas del modelo de gas reticular para comprender una variedad de comportamientos de unión, entre los que se incluyen la unión de ligandos a receptores en la superficie celular ^[32] , la unión de proteínas de quimiotaxis al motor flagelar ^{[33] y la condensación del ADN [34]} .

La actividad de las neuronas en el cerebro se puede modelar estadísticamente. Cada neurona en cualquier momento está activa + o inactiva −. Las neuronas activas son aquellas que envían un potencial de acción a través del axón en cualquier ventana de tiempo dada, y las inactivas son aquellas que no lo hacen.

Siguiendo el enfoque general de Jaynes, ^{[35] [36]} una interpretación posterior de Schneidman, Berry, Segev y Bialek, ^[37] es que el modelo de Ising es útil para cualquier modelo de función neuronal, porque un modelo estadístico para la actividad neuronal debe elegirse utilizando el principio de máxima entropía . Dada una colección de neuronas, un modelo estadístico que puede reproducir la tasa de disparo promedio para cada neurona introduce un multiplicador de Lagrange para cada neurona: Pero la actividad de cada neurona en este modelo es estadísticamente independiente. Para permitir correlaciones de pares, cuando una neurona tiende a disparar (o no disparar) junto con otra, introduzca multiplicadores de Lagrange por pares: donde no están restringidos a los vecinos. Tenga en cuenta que esta generalización del modelo de Ising a veces se denomina distribución binaria exponencial cuadrática en estadística. Esta función de energía solo introduce sesgos de probabilidad para un espín que tiene un valor y para un par de espines que tienen el mismo valor. Las correlaciones de orden superior no están restringidas por los multiplicadores. Un patrón de actividad obtenido a partir de esta distribución requiere la mayor cantidad de bits para almacenarse en una computadora, en el esquema de codificación más eficiente imaginable, en comparación con cualquier otra distribución con la misma actividad promedio y correlaciones por pares. Esto significa que los modelos de Ising son relevantes para cualquier sistema que se describa mediante bits que sean lo más aleatorios posible, con restricciones en las correlaciones por pares y el número promedio de 1, lo que ocurre con frecuencia tanto en las ciencias físicas como en las sociales.$${\displaystyle E=-\sum _{i}h_{i}S_{i}}$$$${\displaystyle E=-{\tfrac {1}{2}}\sum _{ij}J_{ij}S_{i}S_{j}-\sum _{i}h_{i}S_{i}}$$${\displaystyle J_{ij}}$

Con el modelo de Ising, los llamados vidrios de espín también pueden describirse mediante el hamiltoniano habitual , donde las variables S describen los espines de Ising, mientras que las J _i,k se toman de una distribución aleatoria. Para los vidrios de espín, una distribución típica elige enlaces antiferromagnéticos con probabilidad p y enlaces ferromagnéticos con probabilidad 1 −  p (también conocido como el modelo de Ising de enlace aleatorio). Estos enlaces permanecen fijos o "extinguidos" incluso en presencia de fluctuaciones térmicas. Cuando p  = 0 tenemos el modelo de Ising original. Este sistema merece interés por sí mismo; en particular, uno tiene propiedades "no ergódicas" que conducen a un extraño comportamiento de relajación. También se ha prestado mucha atención al modelo de Ising relacionado con el enlace y el sitio diluido, especialmente en dos dimensiones, lo que conduce a un comportamiento crítico intrigante. ^[38]${\textstyle H=-{\frac {1}{2}}\,\sum J_{i,k}\,S_{i}\,S_{k},}$

El modelo de Ising fue instrumental en el desarrollo de la red de Hopfield . El modelo original de Ising es un modelo para el equilibrio. Roy J. Glauber en 1963 estudió el modelo de Ising evolucionando en el tiempo, como un proceso hacia el equilibrio térmico ( dinámica de Glauber ), agregando el componente del tiempo. ^[39] (Kaoru Nakano, 1971) ^{[40] [41]} y ( Shun'ichi Amari , 1972), ^[42] propusieron modificar los pesos de un modelo de Ising por la regla de aprendizaje hebbiana como un modelo de memoria asociativa. La misma idea fue publicada por (William A. Little  [de] , 1974), ^[43] quien fue citado por Hopfield en su artículo de 1982.

El modelo de Sherrington-Kirkpatrick del vidrio de espín, publicado en 1975, ^[44] es la red de Hopfield con inicialización aleatoria. Sherrington y Kirkpatrick descubrieron que es muy probable que la función de energía del modelo SK tenga muchos mínimos locales. En el artículo de 1982, Hopfield aplicó esta teoría recientemente desarrollada para estudiar la red de Hopfield con funciones de activación binarias. ^[45] En un artículo de 1984, la extendió a funciones de activación continua. ^[46] Se convirtió en un modelo estándar para el estudio de redes neuronales a través de la mecánica estadística. ^{[47] [48]}

El estanque de fusión se puede modelar mediante el modelo de Ising; los datos de topografía del hielo marino influyen bastante en los resultados. La variable de estado es binaria para una aproximación 2D simple, ya sea agua o hielo. ^[49]

Con el fin de investigar un modelo de Ising con relevancia potencial para redes neuronales grandes (por ejemplo, con o interacciones por nodo), por sugerencia de Krizan en 1979, Barth (1981) obtuvo la expresión analítica exacta para la energía libre del modelo de Ising en el árbol de Cayley cerrado (con una relación de ramificación arbitrariamente grande) para un campo magnético externo cero (en el límite termodinámico) aplicando las metodologías de Glasser (1970) y Jellito (1979)${\displaystyle 10^{4}}$${\displaystyle 10^{5}}$