Variational optimization of projected entangled-pair states on the triangular lattice

Resumen Técnico: Optimización Variacional de Estados de Pares Entrelazados Proyectados en la Red Triangular

Planteamiento del Problema
La simulación de sistemas de muchos cuerpos cuánticos frustrados, particularmente en redes triangulares y de kagome, sigue siendo un desafío significativo en la física de la materia condensada. Si bien los Estados de Pares Entrelazados Proyectados (PEPS) se han convertido en un marco estándar para simular sistemas cuánticos bidimensionales, su optimización variacional en geometrías triangulares ha estado históricamente constreñida. Los enfoques convencionales suelen mapear la red triangular en una red cuadrada, introduciendo interacciones de siguiente vecino diagonal. Este mapeo obliga a que las correlaciones a lo largo de la dirección diagonal se redirijan a través de las direcciones restantes de la red, lo que constituye una representación ineficiente de la estructura de entrelazamiento intrínseca del sistema. Además, los métodos existentes a menudo dependen de mapeos indirectos que no logran capturar plenamente la simetría rotacional y la conectividad específica de la red de Bravais triangular, lo que potencialmente limita la precisión de las minimizaciones de energía variacional para una dimensión de enlace fija.

Metodología
Los autores presentan un algoritmo de Grupo de Renormalización de Matriz de Transferencia de Esquina (CTMRG) de propósito general diseñado específicamente para PEPS infinitos en la red triangular. Los componentes metodológicos clave incluyen:

• Formulación Triangular Nativa: En lugar de mapear a una red cuadrada, el algoritmo opera directamente sobre la geometría triangular. El ansatz PEPS se define utilizando tensores de siete índices (un índice físico y seis índices virtuales que conectan con todos los vecinos más cercanos), reflejando el número de coordinación de la red triangular.
• Diferenciación Automática (AD): El enfoque integra la Diferenciación Automática (AD) para calcular los gradientes de energía de manera eficiente, permitiendo la minimización directa de la energía variacional sin la necesidad de derivar manualmente reglas de actualización complejas.
• Procedimiento CTMRG Radial: A diferencia del CTMRG de la red cuadrada, que utiliza direcciones ortogonales, la variante triangular emplea un esquema de absorción radial. El entorno consiste en seis tensores de esquina ($C$) y doce tensores de borde ($T$). El algoritmo realiza métodos de potencia iterativos para calcular los tensores del entorno de punto fijo.
• Absorción y Renormalización: La absorción de los tensores PEPS en el entorno se realiza de forma radial. Debido a que los vectores de la red no son ortogonales, las actualizaciones de los tensores de esquina y de borde son interdependientes. El método utiliza un proceso de renormalización de dos pasos que involucra dos conjuntos de proyectores ($P$ y $Q$) para truncar las dimensiones de enlace ampliadas de vuelta a una dimensión de enlace del entorno fija $\chi_E$, asegurando la estabilidad numérica.
• Manejo de Simetrías: El algoritmo no asume la simetría rotacional de los tensores de entrada, lo que le permite acomodar estados arbitrarios, incluyendo estados de rotura de simetría y espirales. También puede extenderse para incorporar simetrías físicas para mejorar la eficiencia.

Contribuciones Clave

1. Desarrollo Algorítmico: El artículo presenta un método computacionalmente eficiente y preciso para contraer y optimizar variacionalmente PEPS directamente en la red triangular, evitando la necesidad de mapeos a redes cuadradas.
2. Representabilidad Mejorada: Al utilizar la geometría nativa, el ansatz de PEPS triangular (TL-PEPS) captura las correlaciones cuánticas a lo largo de tres direcciones de red independientes de manera más fiel que sus contrapartes de red cuadrada (SL-PEPS). El análisis teórico muestra que, para una dimensión de enlace dada $\chi_B$, el límite superior de la entropía de entrelazamiento para TL-PEPS es significativamente mayor (escalando como $(8L-2)\log\chi_B$) en comparación con SL-PEPS (que escala como $4L\log\chi_B$).
3. Aumento de Parámetros Variacionales: Con una dimensión de enlace fija, el ansatz TL-PEPS posee un mayor número de parámetros variacionales ($p\chi_B^6$) en comparación con el SL-PEPS ($p\chi_B^4$), ofreciendo un mayor poder de expresión para los mismos recursos computacionales en términos de dimensión de enlace.
4. Marco de Referencia (Benchmarking): Los autores proporcionan un protocolo de evaluación riguroso que compara el TL-PEPS nativo frente al SL-PEPS mapeado, asegurando comparaciones justas al igualar los parámetros variacionales o analizar la convergencia con respecto a las dimensiones de enlace del entorno.

Resultados
El método se aplicó a dos modelos de espín-1/2 de Heisenberg antiferromagnético paradigmáticos:

• Modelo de Heisenberg de Red Triangular:

  • El enfoque TL-PEPS produjo energías variacionales más bajas que el SL-PEPS para las mismas dimensiones de enlace.
  • Con $\chi_B = 5$, el TL-PEPS alcanzó una energía del estado fundamental de $E = -0.55136$, que es el valor de energía más bajo de red de tensores reportado hasta la fecha para este modelo, superando los resultados previos de extrapolación de DMRG e iPEPS.
  • La extrapolación al límite de dimensión de enlace infinita produjo $E = -0.55164$, consistente con los resultados de SL pero obtenido con propiedades de convergencia mejoradas.
  • El parámetro de orden magnético $m$ mostró una reducción sistemática en TL-PEPS ($m \approx 0.157$) en comparación con SL-PEPS ($m \approx 0.182$), atribuido a la mayor sensibilidad del ansatz a la estructura de entrelazamiento.

• Modelo de Heisenberg de Red de Kagome:

  • El TL-PEPS proporcionó consistentemente energías variacionales más bajas en todas las dimensiones de enlace accesibles en comparación con la formulación SL.
  • La energía del estado fundamental extrapolada fue $E = -0.43929$, mejorando ligeramente los resultados previos de DMRG y de la diagonalización exacta.
  • Los datos de TL-PEPS exhibieron un comportamiento de escalamiento más limpio con respecto al inverso de la longitud de correlación ($1/\xi$), lo que facilita extrapolaciones más fiables incluso en dimensiones de enlace bajas donde los datos de SL-PEPS mostraban inexactitudes.
  • Los resultados apoyan un estado fundamental no magnético consistente con un líquido de espines cuánticos (QSL), aunque la naturaleza con brecha (gapped) o sin brecha (gapless) sigue siendo una cuestión abierta en la literatura.

Significancia y Reivindicaciones
El artículo afirma que la formulación nativa de la red triangular ofrece una ventaja significativa en la precisión variacional para imanes cuánticos frustrados. Las mejoras se atribuyen a:

1. Representación Fiel del Entrelazamiento: El ansatz respeta naturalmente la geometría de la red y las simetrías de rotación, lo que conduce a una representación más precisa de la estructura de entrelazamiento.
2. Eficiencia de Parámetros: El mayor número de parámetros variacionales a una dimensión de enlace fija permite al ansatz explorar el espacio de Hilbert de manera más efectiva.
3. Viabilidad: A pesar del mayor costo computacional por paso (escalando como $O(\chi_B^{13})$ en comparación con la red cuadrada), el método sigue siendo factible y competitivo, particularmente cuando se compara con simulaciones de red cuadrada con un número equivalente de parámetros variacionales.

Los autores concluyen que este método proporciona una perspectiva imparcial sobre problemas abiertos en el magnetismo cuántico frustrado y allana el camino para el estudio de otras redes basadas en la estructura triangular (por ejemplo, panal, dice, maple-leaf) y modelos con interacciones de siguiente vecino sin rutas de conexión complejas. La implementación está disponible en la biblioteca de código abierto variPEPS.