Spin-imbalanced fermion on a dynamic lattice

Resumen Técnico: Fermiones con Desequilibrio de Espín en una Red Dinámica

Planteamiento del Problema
Este estudio investiga las propiedades magnéticas del estado fundamental de un modelo de red dinámica de fermiones de espín-1/2 en una dimensión (1D). Mientras que los sistemas 1D convencionales (como el modelo de Fermi-Hubbard extendido) han establecido que los vectores de orden de la onda de densidad de espín (SDW) siguen la relación $k = 2k_F$ en sistemas equilibrados, el comportamiento de estos sistemas bajo un desequilibrio de espín de fermiones ($N_\uparrow \neq N_\downarrow$) dentro de una red dinámica permanece en gran medida inexplorado. Específicamente, el artículo aborda cómo la interacción entre fermiones itinerantes y espines localizados centrados en los enlaces, mediada por el salto dependiente del espín, influye en la modulación magnética cuando la población de fermiones está parcialmente polarizada.

Metodología
Los autores emplean el método de la Renormalización de Matriz de Densidad (DMRG) para simular el estado fundamental del sistema. El Hamiltoniano del modelo describe a los fermiones itinerantes saltando entre los sitios de la red con una amplitud $t$, la cual es modulada dinámicamente por los espines localizados que residen en los enlaces mediante un término de acoplamiento $\alpha$. Los espines localizados están sujetos a un campo longitudinal $\Delta$ y un campo transversal $\beta$. Los fermiones itinerantes experimentan un campo Zeeman externo $h$ (que controla el desequilibrio de espín) y una interacción de Hubbard repulsiva en el sitio $U$.

Las simulaciones se centran en el régimen quasiadiabático ($\beta \ll t$), fijando $\beta = 0.02$ y $t=1$. Se analizan tamaños de sistema que van desde $L=24$ hasta $L=60$ con condiciones de contorno abiertas. Las observables clave incluyen el factor de estructura de espín localizado $S_z(k)$, la función de correlación en el espacio real $S_z(r)$ y las funciones de distribución de momento resueltas por espín $n_\sigma(k)$. Se utiliza el escalamiento de tamaño finito para extrapolar los resultados al límite termodinámico ($L \to \infty$).

Contribuciones Clave y Resultados

1. Identificación de Fases Magnéticas:
   El estudio mapea el diagrama de fases, identificando dos fases primarias: una fase paramagnética (PM) y una fase de onda de densidad de espín (SDW). La transición entre estas fases es impulsada por la competencia entre el acoplamiento de salto dependiente del espín ($\alpha$) y el campo longitudinal ($\Delta$). Cuando $\alpha$ y $\Delta$ son comparables, emergen órdenes de SDW espacialmente modulados. Cuando un parámetro domina, el sistema entra en una fase PM caracterizada por un alineamiento de espín uniforme y la ausencia de correlaciones magnéticas de largo alcance.

2. Fenómenos de Ordenamiento Dual de SDW (SDW-I y SDW-II):
   En el régimen de espín parcialmente polarizado, los autores descubren dos patrones de ordenamiento SDW distintos, los cuales se desvían fundamentalmente de la expectativa estándar de $k=2k_F$ para sistemas equilibrados:

   • SDW-I: El vector de onda de ordenamiento $k_{max}$ exhibe una dependencia lineal de la densidad de espín mayoritario ($\rho_\uparrow$), siguiendo $k_{max} \approx 2\pi\rho_\uparrow$. Esto corresponde al vector de anidamiento (nesting) de la superficie de Fermi de espín hacia arriba.
   • SDW-II: Por el contrario, $k_{max}$ exhibe una dependencia lineal de la densidad de espín minoritario ($\rho_\downarrow$), siguiendo $k_{max} \approx 2\pi\rho_\downarrow$. Esto corresponde al vector de anidamiento de la superficie de Fermi de espín hacia abajo.
     La transición entre estos dos fenómenos es controlada mediante el ajuste de los parámetros de acoplamiento ($\alpha, \Delta$) y el factor de llenado.

3. Origen Microscópico:
   El análisis de las funciones de distribución de momento resueltas por espín revela que el pico de la SDW dominante está dictado no por la densidad de llenado total, sino por las ocupaciones relativas de los componentes de espín. La emergencia de estos patrones duales se interpreta como un análogo magnético de una transición tipo Peierls, impulsada por el anidamiento de la superficie de Fermi específico de los sectores resueltos por espín de los fermiones itinerantes.

4. Robustez frente a la Interacción de Hubbard:
   Se encuentra que la inclusión de interacciones de Hubbard repulsivas en el sitio ($U$) potencia la tendencia hacia la polarización de espín de los fermiones (reduciendo el campo crítico requerido para la polarización completa) al suprimir la doble ocupación. Sin embargo, la interacción repulsiva no altera cualitativamente la naturaleza de las fases SDW-I y SDW-II; estos patrones de ordenamiento distintos permanecen robustos dentro de sus respectivas regiones de parámetros.

Significancia
El artículo afirma revelar una nueva ruta para sintonizar las modulaciones magnéticas en sistemas correlacionados unidimensionales. Al demostrar que el desequilibrio de espín de los fermiones puede impulsar el sistema hacia distintas fases de SDW caracterizadas por diferentes vectores de anidamiento ($2k_{F\uparrow}$ vs. $2k_{F\downarrow}$), este trabajo enriquece la comprensión microscópica del magnetismo de red dinámica. Los resultados sugieren que el orden magnético en tales sistemas no es meramente una propiedad estática, sino que está acoplado dinámicamente a la configuración de los espines localizados y a la población de espín específica de los fermiones itinerantes. Esto proporciona un marco teórico para comprender cómo el desequilibrio de espín moldea las fases cuánticas en sistemas de baja dimensión más allá del límite convencional equilibrado.