2025-11-16T02:31:12.228614

Tensor network method for solving the Ising model with a magnetic field

Chung
We study the two-dimensional square lattice Ising ferromagnet and antiferromagnet with a magnetic field by using tensor network method. Focusing on the role of guage fixing, we present the partition function in terms of a tensor network. The tensor has a different symmetry property for ferromagnets and antiferromagnets. The tensor network of the partition function is interpreted as a multiple product of the one-dimensional quantum Hamiltonian. We perform infinite density matrix renormalization group to contract the two-dimensional tensor network. We present the numerical result of magnetization and entanglement entropy for the Ising ferromagnet and antiferromagnet side by side. In order to determine the critical line in the parameter space of temperature and magnetic field, we use the half-chain entanglement entropy of the one-dimensional quantum state. The entanglement entropy precisely indicates the critical line forming the parabolic shape for the antiferromagnetic case, but shows the critical point for the ferromagnetic case.
academic

Método de red tensorial para resolver el modelo de Ising con campo magnético

Información Básica

  • ID del Artículo: 2501.01098
  • Título: Método de red tensorial para resolver el modelo de Ising con campo magnético
  • Autor: Myung-Hoon Chung (College of Science and Technology, Hongik University)
  • Clasificación: cond-mat.stat-mech, quant-ph
  • Fecha de Publicación: 3 de enero de 2025
  • Enlace del Artículo: https://arxiv.org/abs/2501.01098

Resumen

Este artículo utiliza el método de red tensorial para estudiar ferromagnetos y antiferromagnetos de Ising con campo magnético en redes cuadradas bidimensionales. Se enfatiza el papel de la fijación de calibre, representando la función de partición en forma de red tensorial. Los tensores del ferromagneto y antiferromagneto poseen propiedades de simetría distintas. La red tensorial de la función de partición se interpreta como un producto múltiple de hamiltonianos cuánticos unidimensionales. Se utiliza el grupo de renormalización de matriz de densidad infinita (iDMRG) para contraer la red tensorial bidimensional. El artículo presenta resultados numéricos comparativos de magnetización y entropía de entrelazamiento para ferromagnetos y antiferromagnetos de Ising. Para determinar las líneas críticas en el espacio de parámetros de temperatura y campo magnético, se utiliza la entropía de entrelazamiento de media cadena de estados cuánticos unidimensionales. La entropía de entrelazamiento indica precisamente una línea crítica de forma parabólica en el caso antiferromagnético, mientras que muestra un punto crítico en el caso ferromagnético.

Antecedentes y Motivación de la Investigación

Contexto del Problema

  1. Desafíos de Complejidad Computacional: La complejidad computacional de sistemas cuánticos de muchos cuerpos crece exponencialmente con el tamaño del sistema y la dimensionalidad, requiriendo métodos numéricos efectivos
  2. Desarrollo de Redes Tensoriales: Desde los años 1990, las técnicas de red tensorial se han convertido en herramientas fundamentales para tratar sistemas de muchos cuerpos, descubiertas independientemente en física estadística y sistemas cuánticos de muchos cuerpos
  3. Importancia del Modelo de Ising: El modelo de Ising bidimensional es un sistema prototípico para investigaciones teóricas de correlaciones cuánticas y mecánica estadística clásica

Motivación de la Investigación

  1. Innovación Metodológica: Explorar la aplicación del método de red tensorial en el modelo de Ising bidimensional con campo magnético
  2. Análisis de Simetría: Investigar el papel de la fijación de calibre en la construcción de redes tensoriales, particularmente las diferentes simetrías en casos ferromagnético y antiferromagnético
  3. Detección de Transiciones de Fase: Utilizar la entropía de entrelazamiento como marcador para detectar transiciones de fase cuántica, una aplicación importante de conceptos de información cuántica en física estadística

Contribuciones Principales

  1. Marco Unificado: Estudiar simultáneamente la representación de red tensorial de modelos de Ising ferromagnético y antiferromagnético, revelando diferencias esenciales en simetría entre ambos
  2. Análisis de Fijación de Calibre: Análisis profundo del papel crucial de la fijación de calibre en la construcción de redes tensoriales, demostrando la equivalencia de diferentes opciones de calibre
  3. Aplicación de Entropía de Entrelazamiento: Uso exitoso de la entropía de entrelazamiento de media cadena para determinar líneas críticas, particularmente descubriendo regiones críticas de forma parabólica en el caso antiferromagnético
  4. Optimización de Métodos Numéricos: Adopción del método iDMRG para contraer efectivamente redes tensoriales bidimensionales, evitando dificultades del método de grano grueso en el manejo de tensores defectuosos

Explicación Detallada de Métodos

Definición de Tareas

Investigar el modelo de Ising en redes cuadradas bidimensionales, incluyendo:

  • Entrada: Parámetro de temperatura β, intensidad de campo magnético B, constante de acoplamiento J
  • Salida: Magnetización, entropía de entrelazamiento, ubicación de líneas críticas
  • Restricciones: Manejo de casos ferromagnético (J>0) y antiferromagnético (J>0 pero interacción repulsiva)

Construcción de Red Tensorial

Representación de Función de Partición

Función de partición del modelo de Ising:

Z = Σ{σ} e^(-βH({σ}))

Donde el hamiltoniano es:

  • Caso ferromagnético: Hf=JijσiσjBiσiH_f = -J\sum_{\langle ij\rangle} \sigma_i\sigma_j - B\sum_i \sigma_i
  • Caso antiferromagnético: Ha=JijσiσjBiσiH_a = J\sum_{\langle ij\rangle} \sigma_i\sigma_j - B\sum_i \sigma_i

Descomposición de Valores Singulares

Descomposición de términos de interacción mediante descomposición de valores singulares:

e^(βJσ_iσ_j) = Σ_m W_σim W_σjm
e^(-βJσ_iσ_j) = Σ_m V_σim W_σjm

Donde:

W = (√cosh(βJ)    √sinh(βJ))
    (√cosh(βJ)   -√sinh(βJ))

V = (√cosh(βJ)   -√sinh(βJ))
    (√cosh(βJ)    √sinh(βJ))

Construcción de Tensor de Cuatro Patas

  • Tensor Ferromagnético: Truldf=σeβBσWσrWσuWσlWσdT^f_{ruld} = \sum_σ e^{βBσ} W_{σr}W_{σu}W_{σl}W_{σd}
  • Tensor Antiferromagnético: Trulda=σeβBσWσrVσuVσlWσdT^a_{ruld} = \sum_σ e^{βBσ} W_{σr}V_{σu}V_{σl}W_{σd}

Observación clave: TruldfT^f_{ruld} posee simetría completa de permutación de índices, mientras que TruldaT^a_{ruld} no posee esta simetría.

Implementación del Algoritmo iDMRG

Método de Matriz de Transferencia

Mapeo de modelo estadístico bidimensional a modelo cuántico unidimensional:

... T|A⟩ T|B⟩ T|A⟩ T|B⟩ ... = λ_max × ... |A⟩|B⟩|A⟩|B⟩ ...

Donde T es el operador de matriz de transferencia de fila a fila.

Representación de Estado de Producto Matricial

El estado de frontera se asume como estado de producto matricial (MPS), determinado por iDMRG con tensores de tres patas A, B y coeficientes de Schmidt λ_i.

Transformación de Calibre

Para el caso antiferromagnético, mediante transformación de calibre:

W_σr V_σu V_σl W_σd → V_σr W_σu W_σl V_σd

Se asegura la consistencia del proceso de contracción superior e inferior.

Configuración Experimental

Configuración de Parámetros

  • Constante de Acoplamiento: J = 1 (fijo)
  • Dimensión de Enlace Ambiental: χ = 60
  • Rango de Temperatura: T ∈ 2.10, 2.40 (ferromagnético), T ∈ 2.18, 2.34 (antiferromagnético)
  • Rango de Campo Magnético: B ∈ -0.10, 0.10 (ferromagnético), B ∈ -1.5, 1.5 (antiferromagnético)

Métricas de Evaluación

  1. Magnetización: Calculada mediante la razón de red tensorial defectuosa a función de partición
  2. Entropía de Entrelazamiento de Media Cadena: Sh=iλi2logλi2S_h = -\sum_i \lambda_i^2 \log \lambda_i^2
  3. Comparación de Temperatura Crítica: Comparación con la solución exacta de Onsager Tc=2.269185314...T_c = 2.269185314...

Comparativas de Referencia

  • Solución exacta de Onsager para modelo de Ising bidimensional (caso de campo cero)
  • Magnetización exacta de Yang para caso ferromagnético
  • Resultados de línea crítica antiferromagnética de Müller-Hartmann y Zittartz

Resultados Experimentales

Resultados Principales

Magnetización

  1. Caso Ferromagnético: Coincide bien con resultados exactos de Yang en campo cero
  2. Caso Antiferromagnético: La magnetización total desaparece por encima de T=2.27 en campo cero, mostrando patrón de tablero por debajo pero con magnetización total cero

Análisis de Entropía de Entrelazamiento

  1. Caso Ferromagnético: Picos agudos cerca de T≈2.27 y campo cero, correspondiendo a punto crítico
  2. Caso Antiferromagnético: Muestra región parabólica, todos los valores de entropía grande casi idénticos, formando línea crítica

Comparativa de Comportamiento Crítico

  • Ferromagneto: Exhibe comportamiento crítico puntual (punto crítico)
  • Antiferromagneto: Exhibe comportamiento crítico lineal (línea crítica de forma parabólica)

Verificación de Precisión Numérica

La temperatura crítica en campo cero T_c = 2.269185314 es altamente consistente con el valor teórico, verificando la precisión del método numérico.

Trabajo Relacionado

Desarrollo de Redes Tensoriales

  1. Desarrollo Histórico: Desde modelo de Baxter, estado AKLT hasta algoritmo DMRG de White
  2. Unificación Teórica: Östlund y Rommer demostraron que el estado base obtenido por DMRG es precisamente un estado de producto matricial
  3. Perspectiva de Información Cuántica: Reinterpretación de DMRG desde perspectiva de entrelazamiento, desarrollo de marcos MPO y MPS

Métodos de Contracción de Red Tensorial Bidimensional

  1. Métodos de Grano Grueso: Grupo de renormalización tensorial, grupo de renormalización tensorial de orden superior, renormalización de red tensorial
  2. Métodos de Tratamiento de Frontera: Estado de producto matricial uniforme variacional, grupo de renormalización de matriz de transferencia angular, iDMRG

Investigación del Modelo de Ising

  • Solución exacta de Onsager (1944)
  • Aplicación de entropía de entrelazamiento en detección de transiciones de fase
  • Investigación de comportamiento crítico en modelo de Ising antiferromagnético

Conclusiones y Discusión

Conclusiones Principales

  1. Diferencias de Simetría: Las representaciones de red tensorial de modelos de Ising ferromagnético y antiferromagnético poseen simetrías esencialmente diferentes
  2. Entropía de Entrelazamiento como Marcador: La entropía de entrelazamiento actúa exitosamente como marcador de transiciones de fase, indicando comportamiento crítico más claramente que la magnetización
  3. Diferencias en Comportamiento Crítico: El ferromagneto muestra punto crítico puntual, el antiferromagneto muestra línea crítica de forma parabólica

Limitaciones

  1. Limitación de Dimensión de Enlace: La dimensión de enlace χ = 60 puede limitar la estabilidad numérica en regiones de campo magnético alto
  2. Complejidad Computacional: El método iDMRG aún enfrenta limitaciones de recursos computacionales al procesar sistemas grandes
  3. Rango de Aplicabilidad: El método es principalmente aplicable a funciones de partición representables como redes tensoriales

Direcciones Futuras

  1. Dimensiones de Enlace Más Altas: Aumentar la dimensión de enlace para mejorar precisión numérica y estabilidad
  2. Extensión a Otros Modelos: Generalizar el método a otros modelos de física estadística
  3. Investigación de Transiciones de Fase Cuántica: Exploración profunda de la aplicación de entropía de entrelazamiento en detección de diversos tipos de transiciones de fase cuántica

Evaluación Profunda

Fortalezas

  1. Innovación Metodológica: Uso ingenioso de fijación de calibre para revelar diferencias esenciales entre modelos ferromagnético y antiferromagnético
  2. Alta Precisión Numérica: Verificación mediante comparación con soluciones exactas conocidas demuestra confiabilidad del método
  3. Perspectiva Física Profunda: La aplicación de entropía de entrelazamiento como marcador de transiciones de fase demuestra integración profunda de información cuántica y física estadística
  4. Presentación Clara: Lógica clara del artículo con detalles técnicos suficientes

Deficiencias

  1. Rango de Parámetros Limitado: Datos inestables en región de campo magnético grande para caso antiferromagnético, sugiriendo necesidad de tiempos de simulación más largos
  2. Análisis Teórico Insuficiente: Falta de explicación teórica sobre por qué la entropía de entrelazamiento puede indicar precisamente líneas críticas
  3. Análisis de Eficiencia Computacional: Comparación insuficiente de eficiencia con otros métodos de contracción de red tensorial bidimensional

Impacto

  1. Contribución Académica: Proporciona caso importante para aplicación de método de red tensorial en física estadística
  2. Valor Metodológico: El método propuesto de detección de transiciones de fase mediante entropía de entrelazamiento posee universalidad
  3. Reproducibilidad: Detalles técnicos suficientes facilitan reproducción y extensión por otros investigadores

Escenarios de Aplicabilidad

  1. Modelos Estadísticos Clásicos Bidimensionales: Aplicable a diversos modelos de física estadística bidimensional representables como redes tensoriales
  2. Investigación de Transiciones de Fase Cuántica: Proporciona herramienta numérica efectiva para investigar transiciones de fase cuántica
  3. Aplicaciones en Ciencia de Materiales: Puede utilizarse para investigar comportamiento de transiciones de fase en materiales magnéticos

Referencias

El artículo cita 26 referencias importantes, abarcando trabajos clásicos en teoría de redes tensoriales, método DMRG, soluciones exactas del modelo de Ising y otros campos clave, proporcionando base teórica sólida para la investigación.


Evaluación General: Este es un artículo de física teórica de alta calidad que aplica exitosamente el método de red tensorial al estudio del modelo de Ising bidimensional con campo magnético, con contribuciones importantes tanto en metodología como en perspectiva física. Particularmente, el enfoque de detectar transiciones de fase mediante entropía de entrelazamiento posee valor académico importante y perspectivas de aplicación prometedoras.