Dynamic Simulations of Strongly Coupled Spin Ensembles for Inferring Nature of Electronic Correlations from Nuclear Magnetic Resonance

Resumo Técnico: Simulações Dinâmicas de Ensembles de Spins Fortemente Acoplados para Inferir a Natureza das Correlações Eletrônicas a partir de Ressonância Magnética Nuclear

Definição do Problema
A Ressonância Magnética Nuclear (RMN) é uma ferramenta crítica para investigar sistemas eletrônicos fortemente correlacionados, como supercondutores não convencionais e líquidos de spin, por meio da interação hiperfina entre spins nucleares e eletrônicos. No entanto, os pacotes de simulação de RMN padrão (ex: PULSEE, SIMPSON) são inadequados para o estudo desses materiais quando correlações eletrônicas fortes e de longo alcance induzem interações efetivas entre spins nucleares. As ferramentas existentes dependem tipicamente de diagonalização exata, que é computacionalmente intratável para sistemas que excedem ~20 spins devido ao crescimento exponencial do espaço de Hilbert, ou focam em interações dipolares de curto alcance usando técnicas de soma de Ewald. Consequentemente, há uma carência de ferramentas numéricas eficientes capazes de simular grandes ensembles ($n > 100.000$) de spins nucleares interagentes onde a interação é mediada por um ambiente eletrônico de campo médio. Essa limitação dificulta a interpretação de assinaturas anômalas de RMN, como formas de eco assimétricas e deslocamentos espectrais dependentes de pulso, que podem surgir de fases eletrônicas exóticas como o estado de Fulde-Ferrell-Larkin-Ovchinnikov (FFLO).

Metodologia
Os autores desenvolveram um pacote de software de código aberto, o Spin Echo Sim, escrito em Julia, C++ e CUDA, para simular a dinâmica de grandes ensembles de spins nucleares sob a influência de interações mediadas por elétrons.

• Modelo Físico: A simulação utiliza uma aproximação de campo médio para integrar fora os graus de liberdade eletrônicos, resultando em uma interação núcleo-núcleo efetiva de longo alcance. O sistema é modelado como uma rede quadrada 2D de núcleos spin-1/2. O Hamiltoniano efetivo para um spin nuclear $i$ inclui um termo de Zeeman e um termo de campo médio $M(i)$, que depende da magnetização espacialmente média dos spins vizinhos, ponderada por uma função de localização $f(r, \xi)$. A função de localização (tipicamente uma Gaussiana) mimetiza a extensão espacial da suscetibilidade eletrônica, enquanto parâmetros dependentes do eixo ($\alpha_x, \alpha_y, \alpha_z$) capturam a anisotropia da interação eletrônica.
• Dinâmica e Dissipação: A evolução temporal é resolvida no referencial rotativo para eliminar a precessão de Larmor de alta frequência, permitindo passos de tempo maiores. A evolução da matriz densidade é tratada no espaço de Liouville para incorporar dissipação não unitária (decoerência) via operadores de salto de Lindblad, simulando efeitos de $T_1$ e $T_2$. A simulação lida com Hamiltonianos dependentes do tempo aproximando-os como independentes do tempo sobre pequenos passos $dt$, validado via critérios de expansão de Magnus.
• Implementação Computacional:
  • Aceleração CUDA: A simulação central utiliza aceleração por GPU para lidar com a massiva paralelização necessária para $n \sim 10^5$ spins.
  • Cálculo de Magnetização Local: Uma abordagem de "stencil" é usada para calcular o campo médio local $M(i)$. Na implementação CUDA, isso é otimizado através do achatamento de arrays e uso de loops paralelos customizados para realizar multiplicação elemento a elemento e soma, evitando o overhead de memória de pré-calcular stencils deslocados exigido na versão em Julia.
  • Exponenciação de Matriz: Como o CUDA carece de uma função de exponencial de matriz integrada, os autores implementaram um algoritmo customizado de escalonamento e quadratura (scaling-and-squaring) com um solver de eliminação gaussiana escrito à mão para as pequenas matrizes densidade (2x2 ou 4x4), alcançando precisão comparável às implementações de dupla precisão em Julia.

Principais Contribuições

1. Desenvolvimento de Software: O lançamento do Spin Echo Sim, um pacote altamente eficiente capaz de simular experimentos de eco de spin de RMN para ensembles que excedem 100.000 spins. A implementação CUDA é relatada como centenas de vezes mais rápida que as técnicas padrão de CPU ou híbridas CPU/GPU (demonstrando um aumento de velocidade >100x sobre a linha de base em Julia).
2. Validação de Aproximações: O artigo fornece uma análise de convergência rigorosa em relação ao tamanho do ensemble ($n$) e ao passo de tempo ($dt$). Demonstra que o modelo converge como $1/\sqrt{n}$ e exponencialmente em relação a $1/(\Omega dt)$, estabelecendo intervalos de parâmetros confiáveis para simulações precisas.
3. Identificação de Mecanismos: O trabalho elucida como correlações eletrônicas de longo alcance e anisotrópicas se manifestam em observáveis de RMN. Especificamente, identifica duas assinaturas primárias:
   • Assimetria Temporal: A indução de estruturas multi-pico e "travamento de fase" (phase locking) no domínio do tempo, onde os spins recorreem em tempos distintos do tempo de eco padrão devido a torques induzidos pela interação.
   • Deslocamentos Espectrais e Buraco de Queima (Hole Burning): A emergência de deslocamentos de frequência dependentes de pulso (variações do deslocamento de Knight) e "buraco de queima" (desaparecimento do peso espectral no centro da linha) no domínio da frequência, impulsionados pelo interplay entre a força de interação e os ângulos de rotação do pulso.

Resultos

• Força e Alcance da Interação: As simulações revelam um crossover suave de comportamento de curto alcance (dominado por flutuações locais) para comportamento de longo alcance (campo médio global). Estruturas de eco multi-pico fortes emergem apenas quando o peso da interação $\omega$ excede aproximadamente $3\Gamma$ (onde $\Gamma$ é a largura de linha).
• Dependência de Pulso: O estudo demonstra que a posição espectral do eco é altamente sensível ao ângulo de rotação do pulso $\theta$. Desvios de um pulso de $\pi/2$ perfeito quebram a simetria, induzindo um torque líquido e um deslocamento de frequência. Este deslocamento é uma sonda direta da anisotropia entre as interações no plano ($\alpha_{x,y}$) e fora do plano ($\alpha_z$).
• Dissipação: A inclusão do acoplamento fora do plano ($\omega_z$) combinado com interações de curto alcance leva a uma forte atenuação do eco, análogo ao resultado de Pennington-Slichter para relaxação transversal em cupratos. No entanto, para alcances de interação mais longos, o sistema torna-se mais resiliente a essa dissipação.
• Aplicação ao FFLO: Os autores aplicam seu modelo ao contexto da supercondutividade FFLO (especificamente CeCoIn$_5$). Eles argumentam que as linhas de forma compostas e os deslocamentos anômalos observados experimentalmente não são meros artefatos de aquecimento de RF, mas indicativos de uma transição para um estado com aumento do comprimento de correlação eletrônica e anisotropia. A simulação sugere que, conforme o sistema entra no estado FFLO, o comprimento de correlação efetivo aumenta, levando a interações anisotrópicas de maior alcance que induzem efeitos de buraco de queima e travamento de fase.

Significância e Alegações
O artigo posiciona o Spin Echo Sim como uma ferramenta inovadora para extrair informações microscópicas de materiais fortemente correlacionados a partir de dados de RMN. Os autores alegam que, ao medir cuidadosamente o deslocamento dependente de pulso e a assimetria temporal dos ecos de spin, pesquisadores podem inferir:

• A anisotropia da suscetibilidade de spin eletrônica.
• O alcance espacial (comprimento de correlação) das correlações eletrônicas.
• A força total de interação mediada por elétrons.

O trabalho serve como um complemento a estudos experimentais e teóricos anteriores, fornecendo o arcabouço computacional necessário para interpretar assinaturas complexas de RMN em fases exóticas da matéria. Os autores enfatizam que sua abordagem de campo médio é apropriada para capturar a dinâmica de ordem principal das interações mediadas por elétrons, evitando a complexidade proibitiva de tratamentos totalmente correlacionados, mantendo-se, contudo, computacionalmente viável para grandes sistemas. O artigo não pretende resolver o problema geral de muitos corpos, mas sim fornecer um método numérico específico e eficiente para uma classe distinta de problemas em espectroscopia de RMN de elétrons correlacionados.