Wigner crystallization in Bernal bilayer graphene

Resumen Técnico: Cristalización de Wigner en grafeno bicapa Bernal

Planteamiento del problema
El artículo investiga la naturaleza y la estabilidad del estado de cristal de Wigner (WC) en el grafeno bicapa Bernal (BBG) bajo la influencia de un campo de desplazamiento perpendicular ($U$). Si bien experimentos recientes en multicapas de grafeno rotadas y no rotadas han revelado fases correlacionadas ricas, la interacción específica entre la estructura de bandas única de BBG (específicamente la dispersión de tipo "sombrero mexicano" inducida por $U$), la curvatura de Berry y la deformación trigonal (trigonal warping) sobre la fase de WC sigue siendo una cuestión teórica abierta. Los autores pretenden determinar los límites de fase del WC en términos de densidad electrónica ($n$), campo de desplazamiento ($U$) y temperatura ($T$), así como caracterizar las propiedades del estado fundamental, incluyendo el magnetismo orbital y el ordenamiento de valles, distinguiéndolos de las fases de WC convencionales.

Metodología
Los autores emplean una descripción de oscilador armónico (HO) de partícula única para modelar el estado de WC. En este marco, cada electrón en la red de cristal de Wigner triangular es tratado como una partícula confinada en un pozo de potencial parabólico creado por la repulsión de Coulomb de sus vecinos. La fuerza de este confinamiento está determinada por la densidad electrónica.

El estudio procede en dos regímenes teóricos principales:

1. Modelo con simetría rotacional (sin deformación trigonal): Los autores primero descartan la deformación trigonal para tratar la dispersión de baja energía como rotacionalmente simétrica alrededor de los puntos $K$ y $K'$. Derivan una ecuación de Schrödinger de partícula única efectiva en el espacio de momentos, incorporando la conexión de Berry ($A(p)$) como un potencial vectorial. Esto permite soluciones analíticas para estimar los espectros de energía, las funciones de onda y el criterio de Lindemann para la fusión.
2. Modelo de deformación trigonal: Los autores incluyen posteriormente los efectos de la deformación trigonal, que rompe la simetría rotacional y divide la dispersión de "sombrero mexicano" en "mini-valles" discretos (un bolsillo central y tres bolsillos laterales). Para este caso, realizan la diagonalización numérica de la ecuación de Schrödinger bidimensional en una rejilla de espacio de momentos para calcular la razón de Lindemann y mapear el diagrama de fases.

La estabilidad del WC se evalúa utilizando el criterio de Lindemann, donde la fusión ocurre cuando el ancho de la raíz cuadrática media de la función de onda del electrón excede una fracción crítica ($\eta_c \approx 0.23$) de la constante de red. Los autores también calculan las integrales de intercambio para estimar las temperaturas de ordenamiento y evalúan el magnetismo mediante el valor de la esperanza del operador de momento angular orbital.

Contribuciones clave y resultados

• Dependencia del campo de desplazamiento y apantallamiento dieléctrico: El artículo establece que, con un campo de desplazamiento cero, el estado de WC queda excluido en BBG debido a la masa de banda evanescente a baja energía y al fuerte apantallamiento dieléctrico interbanda. Una fase de WC solo emerge cuando se aplica un campo de desplazamiento suficiente para abrir un hueco de banda y aplanar el fondo de la banda.
• Magnetismo orbital espontáneo (sin deformación trigonal): En ausencia de deformación trigonal, los autores demuestran que la curvatura de Berry impulsa una transición en el momento angular orbital ($\ell$) de los electrones del WC. A medida que el campo de desplazamiento $U$ aumenta, el flujo de Berry encerrado por la función de onda aumenta. Cuando este flujo excede un valor crítico ($\tilde{\Phi} > 1/2$), el estado fundamental transiciona de $\ell=0$ a $\ell=1$ (o $\ell=-1$ en el valle $K'$). Esta transición se manifiesta como un salto discontinuo en el magnetismo orbital del estado de WC.
• Fusión doblemente reentrante (con deformación trigonal): Al incluir la deformación trigonal, el diagrama de fases exhibe un comportamiento único de "doble reentrada" en función de la densidad electrónica a bajo $U$.
  1. A densidades muy bajas, los electrones se cristalizan en los "bolsillos laterales" de menor energía.
  2. A medida que la densidad aumenta, el WC se funde en un líquido de Fermi (FL) en una densidad crítica $n_{side}^c$.
  3. A una densidad más alta $n^*$, resulta energéticamente favorable que los electrones ocupen el "bolsillo central" de mayor masa, causando que el sistema se vuelva a congelar en un estado de WC.
  4. Finalmente, a una densidad aún más alta $n_{central}^c$, este WC del bolsillo central se funde nuevamente.
• Ordenamiento de mini-valles: La presencia de deformación trigonal rompe la simetría rotacional, haciendo que el momento angular $\ell$ sea mal definido. En su lugar, el estado fundamental se caracteriza por un patrón espacial no trivial de "ordenamiento de mini-valles". Basándose en un trabajo reciente (Ref. [102]), los autores predicen que los electrones se polarizan en mini-valles específicos en un patrón de franjas alternantes en la red de WC.
• Límites de fase: Los autores proporcionan estimaciones numéricas para los límites de fase. La densidad máxima para el estado de WC en BBG se estima en el orden de $1.2 \times 10^{10} \text{ cm}^{-2}$, ocurriendo a un campo de desplazamiento de aproximadamente $0.7 \text{ V/nm}$. La temperatura máxima de fusión se estima en el orden de $10 \text{ K}$ (disminuyendo con la densidad), con la región de reentrada del bolsillo central teniendo una temperatura máxima inferior de $\sim 1 \text{ K}$.

Significado y afirmaciones
El artículo afirma proporcionar un diagrama de fases teórico exhaustivo para el estado de WC en BBG, destacando cómo la geometría cuántica no trivial (curvatura de Berry) y la ruptura de simetría de la red (deformación trigonal) alteran fundamentalmente la naturaleza del cristal de Wigner.

• Efectos orbitales: El trabajo identifica un mecanismo para el magnetismo orbital espontáneo en un WC impulsado puramente por la curvatura de Berry y los campos de desplazamiento, distinto del ordenamiento basado en el espín.
• Comportamiento de fase complejo: La predicción de la fusión "doblemente reentrante" se presenta como una característica novedosa que surge de la competencia entre diferentes bolsillos de banda (central frente a laterales) con distintas masas efectivas y energías.
• Contexto experimental: Los autores señalan que el rango de densidad predicho para el WC en BBG ($n \lesssim 10^{10} \text{ cm}^{-2}$) es inferior a las densidades donde se han observado estados aislantes en experimentos recientes, lo que sugiere que las observaciones experimentales actuales en BBG pueden estar dominadas por el desorden o corresponder a regímenes diferentes. Sugieren que el comportamiento de reentrada y el ordenamiento de mini-valles podrían ser más fácilmente observables en grafeno de apilamiento rhomboédrico con un mayor número de capas ($N$), donde la estructura de bandas es más plana y el flujo de Berry es mayor ($N\pi$).
• Complementariedad: Se establece explícitamente que el artículo es complementario a un trabajo compañero (Ref. [1]) que se centra en el ordenamiento de espín; este trabajo se centra en los efectos orbitales y el diagrama de fases detallado incluyendo la deformación trigonal.

Los autores mantienen un tono modesto respecto a las limitaciones de la aproximación de HO, reconociendo que puede sobreestimar la estabilidad de la fase de WC en comparación con cálculos de Monte Carlo cuánticos más rigurosos, y que el límite de "fonón de Einstein" utilizado aquí no puede reproducir fases con números de Chern no triviales (como el cristal de efecto Hall anómalo cuántico).