Trimers in the Extended Hubbard Model

Resumen Técnico: Inestabilidad de Trímeros en el Modelo de Hubbard Extendido

Planteamiento del Problema
Este trabajo investiga la estabilidad de la fase aislante ferrimagnética de Lieb en el modelo de Hubbard extendido en una cadena de trímeros acoplados bipartita a medio llenado (un electrón por sitio). Si bien el modelo de Hubbard estándar en redes bipartitas con salto de vecino más cercano y repulsión en el sitio $U$ es conocido por exhibir un estado ferrimagnético robusto predicho por el teorema de Lieb, la introducción de una interacción Coulombiana de vecino más cercano $V$ se espera que induzca inestabilidades de carga. El problema central abordado es la evolución de este estado ferrimagnético a medida que $V$ aumenta, identificando específicamente las condiciones bajo las cuales la fase aislante homogénea se vuelve inestable ante la separación de fases (PS) y caracterizando la naturaleza de las fases coexistentes resultantes.

Metodología
Los autores emplean el método de la Densidad de Matriz de Renormalización (DMRG) para resolver el Hamiltoniano de Hubbard extendido para cadenas de trímeros acoplados con condiciones de contorno abiertas. El sistema consiste en $L$ trímeros, donde cada trímero contiene tres sitios etiquetados como $A$, $B_1$ y $B_2$. El Hamiltoniano incluye:

• Amplitud de salto intra-trímero $t$ (fijada en $t=1$).
• Amplitud de salto inter-trímero $t'$ (fijada en $t'=0.8$).
• Repulsión Coulombiana en el sitio $U$.
• Repulsión Coulombiana de vecino más cercano $V$.

Los cálculos se realizaron en cadenas de tamaños de hasta $L=60$ (total de sitios $N=180$), utilizando dimensiones de enlace máximas entre 512 y 2048, con un error de truncamiento de $\sim 1 \times 10^{-6}$. Para sondear las propiedades del estado fundamental, los autores analizaron:

1. Promedios espaciales y distribuciones locales de la ocupación de partículas ($\langle n_i \rangle$) y la ocupación de dobles ($\langle d_i \rangle = \langle n_{i\uparrow} n_{i\downarrow} \rangle$).
2. Espín total por trímero ($S/L$), calculado mediante funciones de correlación de espín en el sector $S_z=0$.
3. Compresibilidad ($\kappa$), derivada de las diferencias de energía de tamaño finito para identificar regímenes incomprensibles (con brecha/gap) versus compresibles (metálicos/separados de fase).
4. Métricas de homogeneidad, específicamente $\Delta_{max}$, definida como la desviación máxima de la ocupación de partículas en cualquier trímero respecto a la densidad homogénea de 3 electrones por trímero.

Contribuciones Clave y Resultados

• Análisis del Límite Atómico ($t=t'=0$): Los autores analizaron primero las configuraciones estáticas en el límite atómico. Identificaron que el aumento de $V$ desde cero favorece una transición de un estado estándar de medio llenado (un electrón por sitio) a un estado donde los sitios $A$ están vacíos y los sitios $B$ están doblemente ocupados. Esta transición ocurre en un valor crítico $V = U/4$, donde la energía del estado de medio llenado ($E/L = 4V$) es igual a la del estado de transferencia de carga ($E/L = U$).

• El Caso $U=0$: Para una repulsión en el sitio evanescente, cualquier $V$ finito induce separación de fases. El sistema se separa en dos fases distintas: una fase rica en dobles en la subred $B$ (aislante), y otra rica en dobles en la subred $A$ (metálica). El espín total por trímero es cero ($S/L = 0$) en este régimen.

• Estabilidad de la Fase Ferrimagnética de Lieb ($U \neq 0$):

  • Para $V=0$, el sistema exhibe el esperado estado aislante ferrimagnético de Lieb con $S/L = 1/2$.
  • A medida que $V$ aumenta, el sistema permanece en una fase aislante homogénea hasta un punto crítico. Sin embargo, hay una transferencia significativa de carga de los sitios $A$ a los sitios $B$, aumentando la densidad de dobles en la subred $B$.
  • Los autores identifican una línea de cruce (crossover) cerca de $V \approx U/4$. Por debajo de esta línea, el sistema es homogéneo. Entre la línea de cruce y la línea de separación de fases, el sistema permanece homogéneo pero exhibe un carácter "rico en dobles" con un aumento en la transferencia de carga, manteniendo el orden ferrimagnético ($S/L = 1/2$).

• Separación de Fases y Coexistencia:

  • Por encima de una línea crítica de separación de fases (PS) (observada en $V \gtrsim U/4$), el estado homogéneo se vuelve inestable. El sistema entra en una región de coexistencia donde existen espacialmente dos fases con distintas ocupaciones de trímero.
  • Para $U \gtrsim 1$: La región de coexistencia consiste en:
    1. Una fase aislante donde los sitios $B$ están predominantemente doblemente ocupados (aproximadamente 4 electrones por trímero), correspondiente a un estado rico en dobles.
    2. Una fase ferromagnética no saturada y metálica donde la distribución de carga es más uniforme pero carece de la condición estricta de medio llenado. Esta fase exhibe un momento magnético neto ($S/L > 0$), derivado de la interacción entre la itinerancia de la carga y el alineamiento ferrimagnético de los momentos en los sitios $A$ y $B$.
  • Para $U \lesssim 1$: La región de coexistencia se asemeja al caso $U=0$, consistiendo en dos fases paramagnéticas ($S/L \approx 0$), una rica en dobles en los sitios $B$ y la otra en los sitios $A$.

• Comportamiento de Alto $U$: Los autores estiman que la línea de cruce y la línea de PS se fusionan aproximadamente en $U \sim 16$. Para $U > 16$, el sistema transiciona directamente desde la fase ferrimagnética homogénea hacia la región de separación de fases, indicado por un salto en los promedios espaciales en lugar de un punto de inflexión.

Significado y Reivindicaciones
El artículo afirma demostrar que la fase aislante ferrimagnética de Lieb es notablemente robusta frente a la repulsión de vecino más cercano $V$, persistiendo más allá del punto donde la densidad de dobles aumenta significativamente. La contribución principal es el mapeo del diagrama de fases para cadenas de trímeros acoplados, distinguiendo entre un régimen ferrimagnético homogéneo "rico en dobles" y la subsiguiente región de separación de fases.

Los autores sostienen que sus resultados aclaran la naturaleza de la inestabilidad del estado fundamental en esta geometría de red específica. Destacan que, por encima de la línea de PS, el sistema soporta una coexistencia única de una fase paramagnética aislante y una fase ferromagnética no saturada metálica para $U$ intermedio o grande. El trabajo proporciona una estimación cuantitativa de las líneas de PS y de cruce utilizando DMRG, confirmando que la transición es impulsada por la competencia entre el costo energético de la doble ocupación ($U$) y la ganancia de energía al evitar la repulsión de vecino más cercano ($V$). Los autores concluyen que estos hallazgos proporcionan una base para futuras investigaciones sobre modelos dopados y los aspectos topológicos de tales transiciones de fase, sin proponer realizaciones experimentales específicas más allá del marco teórico.