Higher Berry curvature, second Chern numbers and magnetoelectric coupling in crystalline insulators

Resumen Técnico: Mayor Curvatura de Berry, Segundos Números de Chern y Acoplamiento Magnetoeléctrico en Aislantes Cristalinos

Planteamiento del Problema
Los invariantes topológicos son esenciales para clasificar los estados fundamentales de muchos cuerpos cuánticos. Mientras que la curvatura de Berry (una 2-forma) y el número de Chern asociado caracterizan con éxito familias de estados fundamentales de dimensión cero (por ejemplo, estados de Bloch en la zona de Brillouin), desarrollos teóricos recientes sugieren una generalización a sistemas espacialmente extendidos. En este marco, las familias de estados de entrelazamiento de corto alcance en $d$ dimensiones espaciales se caracterizan por formas diferenciales de grado superior (por ejemplo, una 3-forma para sistemas 1D). Específicamente, para sistemas unidimensionales parametrizados por un espacio tridimensional, el invariante relevante es el número de Dixmier–Douady–Kapustin–Spodyneiko (DDKS), $\nu$, derivado de una curvatura de Berry de orden superior.

Existe una brecha significativa en la computación explícita y manifiestamente cuantizada de estos invariantes superiores para sistemas genéricos de fermiones libres o interactuantes. Aunque existen expresiones analíticas para modelos específicos de fermiones libres (como el modelo de Dirac 4D), los métodos numéricos para computar el segundo número de Chern ($ch_2$) en cuatro dimensiones suelen depender de la teoría de gauge en red o de respuestas de campos magnéticos en el espacio real. Estos métodos estándar carecen a menudo de "cuantización manifiesta", lo que significa que convergen a valores enteros solo en el límite continuo o requieren un refinamiento extensivo de la malla, en lugar de devolver enteros por construcción. Además, la interpretación física del número DDKS en tres dimensiones espaciales, particularmente su relación con el acoplamiento magnetoeléctrico orbital y el ángulo axión de Chern-Simons ($\theta_{CS}$), sigue siendo una cuestión abierta.

Metodología
Los autores abordan estos problemas reformulando el modelo de red de un aislante de Chern 4D como una familia de cadenas unidimensionales (1D) invariantes por traslación infinitas parametrizadas por la zona de Brillouin tridimensional ($T^3$).

1. Reformulación del Modelo: Partiendo del Hamiltoniano de Dirac 4D, los autores realizan una transformada de Fourier parcial sobre un componente del momento ($k_4$), lo que resulta en una familia de Hamiltonianos 1D $\{H_\lambda\}$ que dependen de los parámetros $\lambda = (\lambda_1, \lambda_2, \lambda_3)$. Estos Hamiltonianos describen cadenas 1D con dos orbitales por celda unidad.
2. Marco Numérico (iMPS): Los estados fundamentales de estas cadenas 1D se aproximan mediante estados de producto de matrices infinitos (iMPS). Los autores utilizan el algoritmo propuesto en la Ref. [13], que generaliza el método de Fukui-Hatsugai-Suzuki para computar la curvatura de Berry de orden superior.
3. Curvatura de Berry Superior Discreta:
   • El espacio de parámetros $T^3$ se discretiza en una rejilla y se triangula en tetraedros.
   • Para cada vértice, se obtiene un estado fundamental iMPS mediante iDMRG.
   • Se computan las funciones de transición entre tensores iMPS vecinos utilizando los autovectores dominantes de matrices de transferencia mixtas.
   • Se define una fase invariante de gauge $\phi(\Delta_2)$ para las caras triangulares, ponderada por el espectro de Schmidt para manejar la variación de la dimensión de enlace.
   • La curvatura de Berry de orden superior $F(\Delta_3)$ para un tetraedro se define como la suma orientada de estas fases sobre su frontera.
   • El número DDKS $\nu$ se obtiene sumando $F(\Delta_3)$ sobre toda la triangulación de $T^3$, produciendo un entero manifiestamente cuantizado.
4. Análisis Comparativo: Los autores computan el segundo número de Chern $ch_2$ para el mismo modelo utilizando un enfoque estándar de teoría de gauge en red (Ref. [34]) y comparan los diagramas de fase resultantes.
5. Acoplamiento Magnetoeléctrico: El estudio investiga la relación entre el número DDKS y el ángulo axión de Chern-Simons $\theta_{CS}$, que gobierna el acoplamiento magnetoeléctrico orbital en 3D. Comparan la variación de $\theta_{CS}$ (derivada analíticamente de la segunda forma de Chern) con la variación de una "fase de Berry superior" $\theta_{HBP}$ (derivada de la construcción iMPS) a través de un ciclo de bombeo.

Contribuciones Clave y Resultados

• Computación Manifiestamente Cuantizada de $ch_2$: Los autores demuestran que el número DDKS, computado mediante el algoritmo de curvatura de Berry superior basado en iMPS, genera un diagrama de fase topológico para el aislante de Chern 4D que es exactamente congruente con el diagrama de fase del segundo número de Chern $ch_2$. Esta concordancia se mantiene a través de diferentes regímenes de masa ($m$), incluyendo regiones donde los cálculos estándar de teoría de gauge en red de $ch_2$ sufren errores significativos de rejilla finita. El método iMPS produce enteros exactos incluso en rejillas gruesas (p. ej., $4^3$), mientras que el método de gauge en red requiere rejillas mucho más finas ($10^4$ puntos) para converger.
• Conexión Teórica vía Fórmula de Künneth: El artículo establece un vínculo riguroso entre el segundo número de Chern 4D y el número DDKS 3D utilizando la fórmula de Künneth. Al descomponer el toro 4D $T^4$ como $S^1 \times T^3$, los autores muestran que la segunda clase de Chern en $T^4$ corresponde al producto del generador de $H^1(S^1, \mathbb{Z})$ y la clase en $H^3(T^3, \mathbb{Z})$ medida por $\nu$. Así, para esta descomposición específica, $ch_2 = \nu$.
• Acoplamiento Magnetoeléctrico en el Modelo de Dirac 4D: Para el modelo de Dirac 4D, los autores encuentran que el cambio en el ángulo axión de Chern-Simons $\Delta\theta_{CS}$ a lo largo de un ciclo de bombeo es reproducido precisamente por el cambio en la fase de Berry superior $\Delta\theta_{HBP}$. Esto sugiere una correspondencia directa entre la fase de Berry superior y la respuesta magnetoeléctrica en sistemas donde los números de Chern de primer orden de las secciones de la zona de Brillouin (BZ) se anulan.
• Contraejemplo en el Aislante de Hopf: Los autores prueban esta correspondencia en el aislante de Hopf, un modelo de dos bandas con un invariante de Hopf no trivial. En este caso, aunque el ángulo de Chern-Simons está cuantizado por el invariante de Hopf, la fase de Berry superior $\theta_{HBP}$ computada no es ni cuantizada ni igual a $\theta_{CS}$. Los autores atribuyen esta discrepancia a la naturaleza "delicada" del invariante de Hopf, el cual es estrictamente bien definido solo para modelos de dos bandas en el espacio de momentos, una condición que puede no preservarse en la descripción iMPS en el límite termodinámico del espacio real.

Significado y Reivindicaciones
El artículo afirma proporcionar la primera computación numérica, manifiestamente cuantizada, del segundo número de Chern para un modelo de red 4D utilizando curvatura de Berry superior e iMPS. Esto valida el marco teórico de que las curvaturas de Berry de orden superior pueden servir como sondas para invariantes topológicos de grado superior en familias de estados espacialmente extendidos.

Los autores posicionan sus resultados dentro del contexto de la geometría cuántica y las propiedades topológicas, señalando que, si bien la correspondencia entre las fases de Berry superiores y el acoplamiento magnetoeléctrico se mantiene para el modelo de Dirac 4D, no es universal, como lo evidencia el caso del aislante de Hopf. Enfatizan que su método es particularmente adecuado para sistemas interactuantes donde las descripciones de un solo cuerpo de Bloch fallan, aunque actualmente se centran en referencias de fermiones libres. El trabajo destaca el potencial de las dimensiones sintéticas y el accionamiento cuasi-periódico en plataformas experimentales para realizar y sondear estos fenómenos topológicos de mayor dimensión.

El artículo concluye modestamente, afirmando que la relación entre la curvatura de Berry superior y la respuesta física en tres dimensiones sigue siendo una cuestión abierta, particularmente respecto a la universalidad de la conexión con $\theta_{CS}$. Sugiere que las fases de Berry superiores pueden capturar información geométrica distinta de las respuestas de Chern-Simons estándar, lo que justifica una mayor investigación sobre su papel en el transporte no lineal y la física de los estados fundamentales.