A discrete formulation for three-dimensional winding number

Este artículo presenta una formulación discreta robusta para calcular el número de enrollamiento tridimensional en sistemas con degeneraciones, definiendo dos versiones de flujo discreto: una práctica que casi siempre está cuantizada y otra modificada que garantiza estrictamente la cuantización entera.

Lo esencial

Imagina que estás intentando contar cuántas veces una pelota de goma se ha enrollado alrededor de un poste. En el mundo de la física, esto se llama un "número de enrollamiento" (winding number). Si la pelota da una vuelta completa, el número es 1; si da dos, es 2. Este concepto es crucial para entender materiales extraños, como superconductores o aislantes topológicos, que tienen propiedades mágicas protegidas por su forma geométrica interna.

El problema es que en la vida real (y en las computadoras), no podemos ver el espacio como una superficie lisa y perfecta. Lo vemos como una cuadrícula de puntos, como los píxeles en una pantalla o las casillas de un tablero de ajedrez.

El desafío: El "puzzle" de los puntos
El autor de este artículo, Ken Shiozaki, se enfrenta a un problema clásico: ¿Cómo contamos esas vueltas si solo tenemos una cuadrícula de puntos y no una superficie continua?

Anteriormente, los científicos intentaban resolver esto siguiendo a cada "banda" de energía (como si fueran hilos de colores) de un punto a otro, asegurándose de que el hilo rojo del punto A se conectara con el hilo rojo del punto B. Pero, ¡ay! A veces, los hilos se cruzan, se mezclan o se vuelven indistinguibles (esto se llama "degeneración"). Cuando eso pasa, el método antiguo se rompe, como intentar seguir una sola hebra en un ovillo de lana enredado.

La solución: El "hueco" mágico (θ-gap)
Shiozaki propone una idea brillante y más sencilla. En lugar de seguir a los hilos individuales, mira el "cielo" donde están los hilos.

Imagina que los valores de energía de tu sistema son como estrellas en un círculo. A veces, hay un espacio vacío entre las estrellas, un "hueco" o gap.

1. La regla del hueco: El autor dice: "Elegimos un hueco en el círculo de estrellas y lo llamamos nuestro punto de referencia".
2. No necesitas seguir los hilos: En lugar de preguntar "¿A dónde fue el hilo rojo?", preguntamos "¿Cuántas estrellas hay entre el hueco A y el hueco B?".
3. Robustez: Esto es genial porque si los hilos se cruzan o se mezclan, el número total de estrellas en ese espacio sigue siendo el mismo. No te importa si se enredaron; solo te importa cuántas hay. Esto hace que el cálculo sea muy resistente a errores y a situaciones caóticas.

Dos versiones de la receta
El autor ofrece dos formas de hacer este cálculo en una computadora:

• La versión "rápida y práctica" (Flux no modificado): Es como usar una regla simple. Mides el espacio localmente en cada cuadradito de tu cuadrícula. En la mayoría de los casos, si tu cuadrícula es lo suficientemente fina (muchos píxeles), esta regla simple te dará el número entero correcto casi siempre. Es como estimar el área de un lago mirando solo las orillas; suele salir bien si miras de cerca.
• La versión "matemáticamente perfecta" (Flux modificado): A veces, la regla simple puede fallar y darte un número como 2.99 en lugar de 3. Para arreglar esto, el autor crea una versión "corregida" que mira un poco más allá, revisando cómo se conectan los cuadraditos vecinos. Esta versión garantiza matemáticamente que el resultado sea siempre un número entero exacto (3, no 2.99). Es como tener un GPS que no solo mira el camino inmediato, sino que verifica la ruta completa para asegurarse de que no te has perdido.

¿Por qué importa esto?
Este método es como pasar de intentar contar las vueltas de una cuerda enredada a mano (difícil y propenso a errores) a usar un escáner que cuenta automáticamente los nudos sin importar cómo se enrede la cuerda.

Permite a los físicos y a los ingenieros:

1. Simular materiales cuánticos complejos en computadoras de manera mucho más rápida y fiable.
2. Estudiar sistemas donde las partículas se comportan de manera caótica o donde hay "cruces" inesperados, sin que el cálculo se rompa.
3. Descubrir nuevos materiales con propiedades topológicas que podrían usarse en computadoras cuánticas futuras.

En resumen, Shiozaki ha creado una nueva "regla de conteo" para el universo cuántico que es tan inteligente que no se confunde cuando las cosas se ponen feas o enredadas, asegurando que siempre obtengamos la respuesta correcta: un número entero que describe la forma oculta de la materia.

Resumen técnico

A continuación presento un resumen técnico detallado del artículo "A discrete formulation for three-dimensional winding number" (Una formulación discreta para el número de enrollamiento tridimensional) de Ken Shiozaki, publicado en el Yukawa Institute for Theoretical Physics.

1. Planteamiento del Problema

El número de enrollamiento tridimensional ($W_3$) es un invariante topológico fundamental definido para aplicaciones de mapas suaves $g: X \to U(N)$, donde $X$ es una variedad tridimensional cerrada y orientada. En física, este invariante es crucial en la teoría de bandas topológicas (por ejemplo, para superconductores tridimensionales con simetría de inversión temporal) y en la teoría de gauge no abeliana (cálculo del número de instantones).

La definición continua clásica es:
$$W_3[g] = \frac{1}{24\pi^2} \int_X \text{Tr}\left[(g^{-1}dg)^3\right]$$

El desafío: En análisis numéricos, la función $g$ a menudo solo se define en un conjunto finito de puntos de una red (lattice). Existen formulaciones discretas para el número de Chern en 2D, pero para $W_3$ en 3D, los métodos existentes (como el de Höckendorf et al.) requieren rastrear y emparejar bandas individuales (autovalores y autovectores) entre puntos adyacentes. Este proceso se vuelve intrincado y propenso a errores cuando el mapa unitario $g$ presenta degeneraciones (accidentales o impuestas por simetría), ya que el ordenamiento de los autovalores puede volverse ambiguo.

2. Metodología Propuesta

El autor propone una formulación discreta alternativa basada en el concepto de "gaps $\theta$" ($\theta$-gaps), evitando la necesidad de rastrear bandas individuales a través de la red.

Conceptos Clave:

• Gap $\theta$: Una matriz $g \in U(N)$ tiene un gap $\theta$ si ninguno de sus autovalores es $e^{i\theta}$. Esto permite definir una rama específica del logaritmo complejo, $\log_\theta z$.
• Construcción Local: En parches donde $g(x)$ mantiene un gap $\theta$, se puede construir una 2-forma $B_\theta$ tal que $H = dB_\theta$ (donde $H$ es la 3-forma de la densidad de winding).
• Invariancia: La diferencia entre las formas $B_{\theta_1}$ y $B_{\theta_2}$ es una derivada total, lo que permite expresar el número de enrollamiento global como una suma de integrales de línea sobre los bordes de las celdas de la red.

Algoritmo Discreto:

1. Discretización: Se aproxima la variedad $X$ mediante una red cúbica $L$.
2. Selección de Parámetros $\theta$: Para cada cubo $c$ de la red, se selecciona un parámetro $\theta_c$ tal que todos los autovalores de $g$ en los vértices de ese cubo mantengan un gap en $\theta_c$. Esto se logra "difuminando" (smearing) los autovalores para encontrar un intervalo libre de espectro.
3. Definición de Marcos de Autoestados: Se definen subespacios de autoestados $\gamma_{\theta_1, \theta_2}(v)$ en cada vértice $v$, que contienen los autovectores cuyos autovalores caen entre los ángulos $\theta_1$ y $\theta_2$.
4. Cálculo del Flujo: El número de enrollamiento se calcula sumando contribuciones de flujo sobre las caras (plaquettes) de la red.

3. Contribuciones Clave y Resultados

El artículo introduce dos versiones de la formulación discreta:

A. Flujo No Modificado ($\Phi_p$)

• Definición: Se calcula utilizando únicamente la información local de la cara (plaquette) y los marcos de autoestados en sus vértices.
  $$\Phi_p \approx \text{Arg} \left[ \prod_{\text{bordes } \partial p} \det(\gamma^\dagger \gamma) \right]$$
• Ventaja: Es altamente práctico y computacionalmente eficiente.
• Limitación: No garantiza estrictamente la cuantización entera en todas las situaciones, aunque el autor demuestra que para mallas suficientemente finas, el valor es casi siempre un entero (la desviación es pequeña).

B. Flujo Modificado ($\tilde{\Phi}_p$)

• Definición: Para garantizar la cuantización estricta, se redefine la conexión de Berry en los bordes. Esta versión incorpora información de los cuatro cubos adyacentes a cada borde de la red para reorganizar la conexión discreta y eliminar elementos fuera de la diagonal que surgen cuando los marcos de autoestados se superponen de manera compleja.
• Resultado: La suma de estos flujos modificados, $\tilde{W}_3^{\text{dis}}[g] = \frac{1}{2\pi} \sum_p \tilde{\Phi}_p$, es estrictamente un número entero.

Validación Numérica:

El autor verificó la fórmula mediante dos modelos:

1. Modelo de Hamiltoniano de Bloch: Un modelo de superconductor tridimensional en un toro. La fórmula discreta reprodujo correctamente los valores de $W_3$ analíticos (0, $\pm 1$, $\pm 2$) dependiendo de los parámetros de masa y acoplamiento.
2. Modelos Aleatorios: Se generaron matrices aleatorias con hopping de corto alcance.
   • Se observó que el flujo no modificado ($\Phi_p$) a menudo no daba un entero exacto en mallas gruesas o con perturbaciones.
   • El flujo modificado ($\tilde{\Phi}_p$) siempre produjo un entero, convergiendo al valor topológico correcto ($W_3 = -1$ en el ejemplo) a medida que aumentaba la resolución de la red.

4. Significado e Impacto

• Robustez ante Degeneraciones: A diferencia de métodos previos que fallan o requieren un emparejamiento manual complejo de bandas en presencia de degeneraciones (puntos excepcionales o simetrías), este método es robusto porque no depende del ordenamiento global de las bandas, sino de la existencia de un gap angular local.
• Simplicidad de Implementación: Elimina la necesidad de algoritmos de seguimiento de bandas (band tracking), simplificando significativamente el código numérico para sistemas de muchos cuerpos o materiales complejos.
• Versatilidad: La formulación se extiende naturalmente a matrices invertibles $GL_N(\mathbb{C})$ mediante descomposición de valores singulares (SVD), permitiendo el cálculo de invariantes topológicos en sistemas donde la unitariedad no está garantizada a priori (como en puntos excepcionales).
• Aplicabilidad Práctica: Aunque la versión modificada es teóricamente necesaria para la cuantización estricta, la versión no modificada es suficiente para la mayoría de las evaluaciones numéricas prácticas con mallas finas, ofreciendo un equilibrio óptimo entre precisión y costo computacional.

En resumen, este trabajo proporciona una herramienta matemática y numérica robusta para calcular invariantes topológicos tridimensionales en sistemas discretos, superando las limitaciones de los métodos de emparejamiento de bandas tradicionales y abriendo la puerta a simulaciones más eficientes de fases topológicas complejas.