Crossover in the Ordered Phase in the Non-Mermin-Wagner-Hohenberg Regime of Spin Models with Long-Range Coupling

Resumen Técnico: Crossover en la Fase Ordenada en el Régimen No-Mermin-Wagner-Hohenberg de Modelos de Espín con Acoplamiento de Largo Alcance

Planteamiento del Problema
Los modelos de espín clásicos con interacciones de largo alcance que decaen como $r^{-\sigma}$ exhiben comportos de fase complejos que difieren significendamente de sus contrapartes de corto alcance. Mientras que el teorema de Mermin-Wagner-Hohenberg (MWH) prohíbe el verdadero orden de largo alcance (LRO) en sistemas de baja dimensión con interacciones de corto alcance, la presencia de interacciones de largo alcance permite un LRO verdadero en el "régimen no-MWH" ($D < \sigma < 2D$, donde $D$ es la dimensión espacial). La investigación previa se ha centrado principalmente en la transición entre las fases ordenada y desordenada o en valores específicos de $\sigma$. Sin embargo, ha faltado una investigación exhaustiva sobre los comportamientos de escala dentro de la verdadera fase ordenada de largo alcance. Específicamente, no está claro si la magnetización en este régimen converge monótonamente a un valor finito o si exhibe regímenes de escala distintos dependiendo del tamaño del sistema y el rango de interacción.

Metodología
Los autores emplearon simulaciones de Monte Carlo utilizando el algoritmo de Metropolis de actualización única para estudiar modelos XY y Heisenberg clásicos en una dimensión (1D) y dos dimensiones (2D) con acoplamiento de largo alcance $J/|i-j|^\sigma$.

• Modelos: Las simulaciones se realizaron en cadenas 1D y redes cuadradas 2D con condiciones de contorno periódicas. El Hamiltoniano incluye interacciones que decaen como $r^{-\sigma}$.
• Parámetros de Simulación: Los tamaños del sistema ($L$) variaron hasta $L=4000$ (1D) y $L=120$ (2D). Las simulaciones consistieron en $10^5$ barridos de calentamiento seguidos de $5 \times 10^6$ barridos de muestreo. Se tuvo especial cuidado a bajas temperaturas ($T < 0.1$) para mitigar la metaestabilidad restringiendo los giros de espín propuestos.
• Observables: El estudio analizó la energía ($E$), la magnetización ($M$), la susceptibilidad magnética ($\chi$) y los cúmulos magnéticos (razones de Binder, $U_n$). Crucialmente, los autores utilizaron las primeras derivadas de los cúmulos y de la magnetización con respecto a la temperatura inversa ($\beta$) para localizar las transiciones de fase con mayor precisión que los métodos tradicionales de intersección.
• Análisis de Escala: Se aplicó el escalamiento de tamaño finito (FSS) para extraer exponentes críticos ($\nu, \gamma$) y temperaturas críticas ($T_c$). Los autores examinaron específicamente la dependencia de la magnetización con $\ln L$ y el comportamiento de las funciones de correlación espín-espín $G(r)$ para distinguir entre diferentes regímenes ordenados.

Contribuciones Clave y Resultados

1. Identificación de Distintos Regímenes de Escala (EnLRO y ReLRO):
   El hallazgo principal es que la verdadera fase de orden de largo alcance no es uniforme, sino que comprende dos regímenes de escala distintos:

   • Orden de Largo Alcance Mejorado (EnLRO): En este régimen, la magnetización $M$ aumenta con el tamaño del sistema $L$. Esto ocurre en valores de $\sigma$ más pequeños (interacciones de largo alcance más fuertes).
   • Orden de Largo Alcance Reducido (ReLRO): En este régimen, la magnetización $M$ disminuye con el tamaño del sistema $L$ (aunque permanece finita), ocurriendo en valores de $\sigma$ más grandes (interacciones de largo alcance más débiles).
   • Crossover: Un crossover entre estos regímenes ocurre en un $\sigma_c$ crítico. Cerca de este punto, $M$ exhibe un comportamiento no monótono con $L$ (disminuyendo y luego aumentando).

2. Cuantificación de los Puntos de Crossover:

   • Modelo XY 1D: El crossover ocurre en $\sigma_c \approx 1.575$ a bajas temperaturas ($T=0.01$). La temperatura de crossover aumenta con $\sigma$.
   • Modelo XY 2D: El crossover ocurre en $\sigma_c \approx 3.2$ a bajas temperaturas.
   • Modelos de Heisenberg: Se observaron crossovers similares en modelos de Heisenberg 1D y 2D en $\sigma_c \approx 1.575$ y $\sigma_c \approx 3.22$, respectivamente. Los autores sugieren que esta similitud indica que el fenómeno es una característica genérica de los modelos de espín continuos con interacciones de largo alcance, probablemente debido a la simetría de rotación continua compartida.

3. Comportamiento de la Función de Correlación:

   • En el régimen ReLRO, las funciones de correlación espín-espín decaen algebraicamente ($G(r) \sim r^{-\kappa}$) hacia un valor no nulo.
   • En el régimen EnLRO, el decaimiento se describe mejor mediante una función exponencial ($G(r) \sim k_1 + k_2 e^{-r/\xi}$) que se aproxima a una constante no nula. La longitud de correlación $\xi$ muestra un cambio discontinuo cerca de $\sigma_c$.

4. Diagramas de Fase y Exponentes Críticos:
   Se construyeron diagramas de fase exhaustivos en el plano $T-\sigma$ para ambos modelos 1D y 2D. Los autores extrajeron los exponentes críticos $\nu$ y $\gamma$ para las transiciones de orden-a-desorden. Por ejemplo, en el modelo 1D a $\sigma=1.1$, $\nu \approx 1.525$, mientras que a $\sigma=1.8$, $\nu \approx 2.65$. La temperatura crítica $T_c$ diverge cuando $\sigma \to D$.

5. Mecanismo de Crossover:
   El artículo propone un mecanismo basado en la competencia entre el alineamiento de espines, las excitaciones de ondas de espín y los defectos locales:

   • Baja Temperatura: La disminución de $\sigma$ suprime la excitación y propagación de ondas de espín, impulsando el sistema de ReLRO a EnLRO.
   • Alta Temperatura: Los defectos locales (que rompen la cadena o la red en subsistemas) bloquean la propagación de las ondas de espín. A medida que el tamaño del sistema aumenta, los defectos proliferan más rápido que las perturbaciones de las ondas de espín, aislando efectivamente las regiones ordenadas y mejorando el orden local, lo que conduce a EnLRO.

Significancia y Reivindicaciones
El artículo afirma refinar la comprensión de los modelos de espín clásicos en el régimen no-MWH al identificar que la fase ordenada contiene distintos comportamientos de escala previamente pasados por alto. Los autores postulan que la distinción entre EnLRO y ReLRO es una característica genérica de los modelos de espín continuos.

El estudio destaca que, si bien el tránsito EnLRO-ReLRO es un cambio en el comportamiento de escala y no una transición de fase termodinámica tradicional (con calor latente o ruptura de simetría), representa un cambio significativo en la respuesta del sistema a los efectos de tamaño finito. Los autores señalan que, aunque actualmente no está claro si este crossover constituye una verdadera transición de fase con un parámetro de orden aún no identificado, la caracterización de estas diferencias de escala proporciona puntos de referencia valiosos para sistemas experimentales, como las matrices de átomos de Rydberg, que operan en regímenes de tamaño finito donde estos efectos son prominentes. El trabajo cierra la brecha entre los modelos teóricos y las restricciones de tamaño finito observadas en las realizaciones experimentales actuales de sistemas con interacciones de largo alcance.