Finite-temperature formation of magnetic plateaus and simplex liquid states on the frustrated ruby lattice

Resumo Técnico: Formação de Patamares Magnéticos e Estados de Líquido de Símplex em Temperatura Finita na Rede Ruby Frustrada

Enunciado do Problema
A frustração geométrica em sistemas quânticos é conhecida por estabilizar fases não convencionais, como líquidos de spin quânticos, que evitam o ordenamento magnético tradicional em baixas temperaturas. No entanto, compreender a formação desses estados em temperaturas finitas e suas assinaturas termodinâmicas permanece um desafio significativo. Especificamente, o estado fundamental a temperatura zero ($T=0$) do antiferromagneto de Heisenberg de spin-1/2 na rede ruby, com interações de vizinho próximo, foi mostrado anteriormente por exibir patamares magnéticos fracionários ($m_z = 0, 1/3, 1/2, 2/3$). Esses patamares foram interpretados como fases cristalinas discretas que quebram simetrias de rede. Questões em aberto permaneciam em relação à degenerescência exata do estado fundamental, a estabilidade desses patamares em temperatura finita e se o sistema sofre transições de fase convencionais para alcançá-los. Resolver esses problemas requer simulações que unam resultados computacionais com realizações experimentais onde flutuações térmicas e quânticas competem.

Metodologia
Os autores empregam um ansatz de estado de rede de tensores infinito purificado (iTNS) para representar a raiz quadrada não normalizada da matriz de densidade de equilíbrio, $\rho^{1/2}(\beta) = \exp(-\beta H/2)$, onde $\beta = 1/T$. O estudo foca no caso isotrópico ($J_t = J_h = J_d = J$) do Hamiltoniano de Heisenberg na rede ruby.

• Representação de Rede: Os vértices das unidades triangulares são agrupados em triângulos (símplex), mapeando a rede ruby para uma geometria de colmeia (honeycomb) de escala macroscópica com número de coordenação $z=3$. O iTNS utiliza uma célula unitária de dois sítios, onde cada tensor possui seis índices físicos (representando os três spins em cada triângulo) e três índices virtuais.
• Algoritmo: O sistema é inicializado na temperatura infinita ($\rho^{1/2}(0) = I$) e resfriado via evolução no tempo imaginário usando uma Trotterização de segunda ordem do propagador.
• Observáveis: Quantidades termodinâmicas são calculadas usando duas abordagens:
  1. Propagação de Crença (BP): A contração padrão de BP é aplicada à rede de norma para calcular valores esperados.
  2. BP Corrigida por Loop: Para melhorar a precisão além da aproximação de campo médio da BP, os autores aplicam correções de cluster de primeira ordem (correções de loop) à densidade de energia livre. Isso contabiliza as correlações ao redor dos hexágonos da rede.
• Escala: As simulações utilizam hardware de GPU para atingir dimensões de ligação até $\chi \sim 200$, permitindo a medição precisa de observáveis termodinâmicos no limite termodinâmico através de uma ampla faixa de temperatura.

Principais Resultados

1. Estado de Líquido de Símplex em Campo Zero:

   • Ao resfriar a partir da temperatura infinita, o sistema não sofre uma transição de fase para um ordenamento cristalino. Em vez disso, forma um "estado de líquido de símplex" desordenado, caracterizado por uma entropia residual não nula ($s_0 \approx \frac{1}{3}\lambda$, onde $\lambda$ é a entropia residual de coberturas de dímeros em uma rede de colmeia).
   • A capacidade térmica ($C_v$) permanece finita e contínua em todas as temperaturas, sem picos divergentes indicativos de uma transição de fase.
   • A análise da matriz de densidade reduzida revela um autovalor principal $\lambda_1 \approx 1/3$, correspondendo a um estado onde dois símplex estão fortemente pareados. Os demais autovalores correspondem a estados fracamente correlacionados. O sistema é uma mistura desordenada perfeita de um número exponencialmente grande de configurações cristalinas $|\psi_d\rangle$ envolvendo pares de spins de símplex.
   • A densidade de energia deste estado líquido ($e/J \approx -0.498$) coincide com a energia do estado fundamental de configurações cristalinas específicas encontradas em estudos anteriores de $T=0$, mas o sistema permanece em um estado líquido com invariância de translação em vez de selecionar um único cristal.

2. Patamares Magnéticos e Comportamento Induzido por Campo:

   • Na presença de um campo magnético externo, o sistema exibe patamares magnéticos estáveis nas densidades de magnetização $m_z = 0, 1/3, 1/2$ e $2/3$.
   • Assinaturas Termodinâmicas: A formação desses patamares ocorre suavemente, sem transições de fase convencionais. A capacidade térmica permanece contínua, embora picos rasos sejam observados.
   • Entropia Residual: Os patamares abrigam "estados de líquido de símplex" com degenerescência extensa. O patamar $m_z = 1/2$, em particular, exibe uma grande entropia residual (aproximadamente o dobro do patamar $m_z=0$).
   • Estabilidade: O patamar $m_z = 0$ é o mais estável contra flutuações térmicas. O patamar $m_z = 1/3$ é menos estável, aparecendo apenas em uma estreita faixa de campo em temperaturas muito baixas ($T/J \le 10^{-2}$). O patamar $m_z = 2/3$ é observado em cálculos de estado fundamental, mas requer temperaturas menores do que as atualmente acessíveis na simulação de temperatura finita para ser totalmente resolvido.
   • Troca de Fase: Na interface entre os patamares $m_z = 1/2$ e $m_z = 1/3$, o sistema exibe uma troca termicamente impulsionada entre estados líquidos. Isso é marcado por um pico "lambda" assimétrico na capacidade térmica e um aumento abrupto na suscetibilidade magnética, provavelmente refletindo flutuações críticas próximas a uma transição de fase quântica de temperatura zero.

Significância e Alegações
O artigo demonstra que técnicas de redes de tensores, especificamente estados de redes de tensores infinitos otimizados com propagação de crença e correções de loop, fornecem uma rota poderosa e escalável para compreender magnetos quânticos frustrados em temperaturas finitas.

• Validação do Método: O trabalho valida a abordagem de redes de tensores baseada em BP como capaz de capturar propriedades termodinâmicas com precisão ao longo de muitos ordens de magnitude de temperatura, incluindo a detecção de gaps de energia e entropias residuais.
• Natureza do Estado Fundamental: O estudo esclarece que os patamares magnéticos na rede ruby não são simples ordenamentos cristalinos, mas sim "estados de líquido de símplex" — misturas desordenadas de exponencialmente muitas configurações cristalinas que retêm entropia residual não nula e não quebram simetrias de rede.
• Ausência de Transições de Fase: Uma descoberta fundamental é que o sistema atinge esses estados magnéticos estáveis sem passar por transições de fase de temperatura finita, conforme evidenciado pela capacidade térmica contínua e finita.
• Direções Futuras: Os autores sugerem que sua metodologia pode ser estendida para estudar sistemas frustrados tridimensionais, como a rede pirocloro, e para investigar perturbações que possam quebrar a degenerescência do estado fundamental em favor de superposições de ordem topológica.