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

Resumen Técnico: Simulaciones Dinámicas de Ensamble de Espines Fuertemente Acoplados para Inferir la Naturaleza de las Correlaciones Electrónicas a partir de la Resonancia Magnética Nuclear

Planteamiento del Problema
La Resonancia Magnética Nuclear (RMN) es una herramienta crítica para sondear sistemas de electrones fuertemente correlacionados, tales como superconductores no convencionales y líquidos de espín, mediante la interacción hiperfina entre los espines nucleares y electrónicos. Sin embargo, los paquetes de simulación de RMN estándar (por ejemplo, PULSEE, SIMPSON) no son adecuados para estudiar estos materiales cuando las correlaciones electrónicas de largo alcance y fuerte intensidad inducen interacciones efectivas entre los espines nucleares. Las herramientas existentes suelen depender de la diagonalización exacta, la cual es computacionalmente intratable para sistemas que exceden los ~20 espines debido al crecimiento exponencial del espacio de Hilbert, o bien se centran en interacciones dipolares de corto alcance utilizando técnicas de sumación de Ewald. En consecuencia, existe una falta de herramientas numéricas eficientes capaces de simular grandes ensambles ($n > 100,000$) de espines nucleares interactuantes donde la interacción es mediada por un entorno electrónico de campo medio. Esta limitación dificulta la interpretación de firmas de RMN anómalas, como las formas de eco asimétricas y los desplazamientos espectrales dependientes de los pulsos, que pueden surgir de fases electrónicas exóticas como el estado de Fulde-Ferrell-Larkin-Ovchinnikov (FFLO).

Metodología
Los autores desarrollaron un paquete de software de código abierto, Spin Echo Sim, escrito en Julia, C++ y CUDA, para simular la dinámica de grandes ensambles de espines nucleares bajo la influencia de interacciones mediadas por electrones.

• Modelo Físico: La simulación emplea una aproximación de campo medio para integrar las libertades de grado electrónico, lo que resulta en una interacción núcleo-núcleo efectiva de largo alcance. El sistema se modela como una red cuadrada 2D de núcleos de espín-1/2. El Hamiltoniano efectivo para un espín nuclear $i$ incluye un término de Zeeman y un término de campo medio $M(i)$, el cual depende de la magnetización promediada espacialmente de los espines vecinos, ponderada por una función de localización $f(r, \xi)$. La función de localización (típicamente una Gaussiana) imita la extensión espacial de la susceptibilidad electrónica, mientras que los parámetros dependientes del eje ($\alpha_x, \alpha_y, \alpha_z$) capturan la anisotropía de la interacción electrónica.
• Dinámica y Disipación: La evolución temporal se resuelve en el marco rotatorio para eliminar la precesión de Larmor de alta frecuencia, lo que permite pasos de tiempo más grandes. La evolución de la matriz de densidad se trata en el espacio de Liouville para incorporar la disipación no unitaria (decoherencia) mediante operadores de salto de Lindblad, simulando los efectos de $T_1$ y $T_2$. La simulación maneja Hamiltonianos dependientes del tiempo aproximándolos como independientes del tiempo sobre pequeños pasos $dt$, validado mediante criterios de expansión de Magnus.
• Implementación Computacional:
  • Aceleración CUDA: La simulación central utiliza aceleración por GPU para manejar la masiva paralelización requerida para $n \sim 10^5$ espines.
  • Cálculo de Magnetización Local: Se utiliza un enfoque de "stencil" para calcular el campo medio local $M(i)$. En la implementación de CUDA, esto se optimiza mediante el aplanamiento de arreglos y el uso de bucles paralelos personalizados para realizar la multiplicación y suma elemento a elemento, evitando la sobrecarga de memoria de pre-calcular stencils desplazados requeridos en la versión de Julia.
  • Exponenciación de Matrices: Dado que CUDA carece de una función de exponencial de matriz integrada, los autores implementaron un algoritmo personalizado de escalado y cuadratura con un resolvedor de eliminación gaussiana escrito a mano para las matrices pequeñas (2x2 o 4x4), logrando una precisión comparable a las implementaciones de doble precisión en Julia.

Contribuciones Clave

1. Desarrollo de Software: El lanzamiento de Spin Echo Sim, un paquete altamente eficiente capaz de simular experimentos de eco de espín de RMN para ensambles que exceden los 100,000 espines. Se reporta que la implementación de CUDA es cientos de veces más rápida que las técnicas estándar de CPU o híbridas de CPU/GPU (demostrando una aceleración >100x sobre la línea base de Julia).
2. Validación de Aproximaciones: El artículo proporciona un análisis de convergencia riguroso respecto al tamaño del ensamble ($n$) y el paso de tiempo ($dt$). Demuestra que el modelo converge como $1/\sqrt{n}$ y exponencialmente con respecto a $1/(\Omega dt)$, estableciendo rangos de parámetros fiables para simulaciones precisas.
3. Identificación de Mecanismos: El trabajo elucida cómo las correlaciones electrónicas anisotrópicas y de largo alcance se manifiestan en los observables de RMN. Específicamente, identifica dos firmas primarias:
   • Asimetría Temporal: La inducción de estructuras de múltiples picos y "bloqueo de fase" en el dominio del tiempo, donde los espines se recoheren en tiempos distintos al tiempo de eco estándar debido a los torques inducidos por la interacción.
   • Desplazamientos Espectrales y Perforación de Huecos (Hole Burning): La emergencia de desplazamientos de frecuencia dependientes de los pulsos (variaciones del desplazamiento de Knight) y la "perforación de huecos" (desaparición del peso espectral en el centro de la línea) en el dominio de la frecuencia, impulsados por la interacción entre la fuerza de la interacción y los ángulos de giro de los pulsos.

Resultados

• Fuerza y Alcance de la Interacción: Las simulaciones revelan una transición suave desde un comportamiento de corto alcance (dominado por fluctuaciones locales) hacia un comportamiento de largo alcance (campo medio global). Estructuras de eco de múltiples picos fuertes emergen solo cuando el peso de la interacción $\omega$ excede aproximadamente $3\Gamma$ (donde $\Gamma$ es el ancho de línea).
• Dependencia del Pulso: El estudio demuestra que la posición espectral del eco es altamente sensible al ángulo de giro del pulso $\theta$. Las desviaciones de un pulso de $\pi/2$ perfecto rompen la simetría, induciendo un torque neto y un desplazamiento de frecuencia. Este desplazamiento es una sonda directa de la anisotropía entre las interacciones en el plano ($\alpha_{x,y}$) y fuera del plano ($\alpha_z$).
• Disipación: La inclusión del acoplamiento fuera del plano ($\omega_z$) combinado con interacciones de corto alcance conduce a una fuerte atenuación del eco, análoga al resultado de Pennington-Slichter para la relajación transversal en cupratos. Sin embargo, para rangos de interacción más largos, el sistema se vuelve más resistente a esta disipación.
• Aplicación a FFLO: Los autores aplican su modelo al contexto de la superconductividad FFLO (específicamente CeCoIn$_5$). Argumentan que las líneas de forma compuestas y los desplazantes anómalos observados experimentalmente no son meros artefactos del calentamiento por RF, sino que son indicativos de una transición a un estado con un incremento en la longitud de correlación electrónica y la anisotropía. La simulación sugiere que, a medida que el sistema entra en el estado FFLO, la longitud de correlación efectiva aumenta, lo que conduce a interacciones anisotrópicas de mayor alcance que inducen efectos de perforación de huecos y bloqueo de fase.

Significancia y Reivindicaciones
El artículo posiciona a Spin Echo Sim como una herramienta novedosa para extraer información microscópica de materiales fuertemente correlacionados a partir de datos de RMN. Los autores afirman que, mediante la medición cuidadosa del desplazamiento dependiente del pulso y la asimetría temporal de los ecos de espín, los investigadores pueden inferir:

• La anisotropía de la susceptibilidad de espín electrónica.
• El alcance espacial (longitud de correlación) de las correlaciones electrónicas.
• La fuerza total de la interacción mediada por electrones.

Este trabajo sirve como compañero de estudios experimentales y teóricos previos, proporcionando el marco computacional necesario para interpretar firmas complejas de RMN en fases exóticas de la materia. Los autores enfatizan que su enfoque de campo medio es apropiado para capturar la dinámica de orden principal de las interacciones mediadas por electrones, evitando la complejidad prohibitiva de los tratamientos totalmente correlacionados, permaneciendo al mismo tiempo computacionalmente factible para sistemas grandes. El artículo no pretende resolver el problema de muchos cuerpos en general, sino proporcionar un método numérico específico y eficiente para una clase distinta de problemas en la espectroscopia de RMN de electrones correlacionados.