Entanglement dynamics of monitored noninteracting fermions on graphics processing units

Mediante el uso de técnicas de GPU para simular sistemas de fermiones no interactuantes monitoreados en tamaños de red sin precedentes, este trabajo confirma la ausencia de transiciones de fase inducidas por mediciones en una dimensión y descubre una transición crítica en dos dimensiones con exponentes que desafían las predicciones del modelo sigma no lineal.

Lo esencial

Imagina que tienes un enorme tablero de ajedrez donde cada casilla puede tener una ficha (un electrón) o estar vacía. Ahora, imagina que este tablero está "vivo": las fichas se mueven, saltan y bailan entre sí siguiendo las reglas de la mecánica cuántica. Este es nuestro sistema de fermiones no interactuantes (electrones que no se pelean entre sí, solo se mueven).

El problema que resuelve este artículo es como intentar entender cómo se "entrelazan" (se conectan mágicamente) estas fichas cuando alguien las está espiando constantemente.

Aquí te explico la historia, los descubrimientos y por qué fue tan difícil, usando analogías sencillas:

1. El Gran Dilema: ¿Espejo o Laberinto?

En el mundo cuántico, hay dos formas en las que la información puede comportarse:

• La Ley del Volumen (El Laberinto): Si el sistema es libre, la información se mezcla por todo el tablero. Es como si todas las fichas estuvieran conectadas entre sí; para saber qué pasa en una esquina, necesitas mirar todo el tablero.
• La Ley del Área (El Espejo): Si alguien mide (espiando) el sistema muy a menudo, rompe esas conexiones mágicas. La información se queda "atrapada" en pequeñas zonas locales. Es como si cada ficha solo pudiera hablar con sus vecinos inmediatos.

Los científicos llevaban años discutiendo: ¿Cuánto hay que espiar para que el sistema pase de ser un laberinto gigante a ser un conjunto de espejos pequeños? A esto le llaman una "transición de fase inducida por medición".

2. El Problema de las "Gafas de Realidad Virtual" (El Hardware)

Antes de este estudio, los científicos intentaban simular esto en computadoras normales (CPU). El problema es que para ver la respuesta correcta, necesitas un tablero enorme.

• Imagina que intentas ver si un río se seca o no. Si solo miras un charco pequeño (un tablero de 100x100), el agua parece que se va a secar porque se evapora rápido. Pero si miras el río entero (un tablero gigante), ves que el agua fluye eternamente.
• Los estudios anteriores usaban "charcos" (tableros pequeños, menos de 1000 casillas). Por eso, algunos decían que el río se secaba (había una transición) y otros decían que no. Los resultados eran confusos.

La solución de este equipo: Usaron GPUs (tarjetas gráficas, como las de las consolas de videojuegos o para renderizar películas).

• La analogía: Si una computadora normal es como un solo trabajador que mueve fichas una por una, una GPU es como un ejército de 100 trabajadores moviendo fichas al mismo tiempo.
• Gracias a esto, pudieron simular tableros gigantescos:
  • En 1D (una línea): ¡Más de 16,000 casillas! (Antes solo llegaban a 1,000).
  • En 2D (un cuadrado): Un tablero de 160 x 160.

3. Los Descubrimientos: Lo que pasó en 1D y 2D

En una Dimensión (Una línea recta)

• La pregunta: Si espiamos a los electrones en una línea, ¿se rompe el entrelazamiento?
• El resultado: No.
• La explicación: Con los tableros pequeños, parecía que sí se rompía. Pero al usar el tablero gigante (16,000 casillas), vieron que la "correlación" (la conexión entre fichas) nunca se corta completamente, solo se hace muy larga. Es como si el hilo que une a las fichas fuera tan largo que parecía que no existía, pero en realidad, si esperas lo suficiente, siempre hay una conexión.
• Conclusión: En 1D, no importa cuánto espiemos, nunca hay una transición de fase. El sistema siempre mantiene su "volumen" de información.

En Dos Dimensiones (Un cuadrado)

• La pregunta: ¿Qué pasa si el tablero es un cuadrado?
• El resultado: Sí, hay una transición.
• La explicación: Aquí el sistema es más complejo. Si espiamos con una intensidad moderada, el sistema cambia drásticamente.
  • Si espiamos poco: Es un laberinto gigante (Ley del Volumen).
  • Si espiamos mucho: Se convierte en pequeños espejos aislados (Ley del Área).
  • El punto crítico: Existe un momento exacto (un "punto de inflexión") donde el sistema se vuelve caótico y las conexiones se vuelven "invariantes de escala" (se ven iguales sin importar cuánto te alejes).
• Lo sorprendente: Los matemáticos (usando un modelo llamado NLSM) habían predicho que esto pasaría, pero fallaron al predecir los números exactos. Dijeron que el punto de cambio sería uno, pero los datos reales mostraron que es diferente. Necesitamos la simulación gigante para ver la verdad.

4. ¿Por qué es importante esto?

Imagina que estás construyendo una computadora cuántica. Quieres que la información se mantenga entrelazada para hacer cálculos mágicos, pero el entorno (el calor, las vibraciones) actúa como un "espía" que mide tu sistema y destruye la magia.

Este estudio nos dice:

1. En sistemas pequeños (1D): No te preocupes tanto por el espionaje, el sistema es muy robusto.
2. En sistemas grandes (2D): ¡Cuidado! Hay un umbral exacto. Si cruzas esa línea de "espionaje", pierdes la magia cuántica de golpe.

Resumen en una frase

Los científicos usaron la potencia bruta de las tarjetas gráficas para simular sistemas cuánticos gigantescos, demostrando que en una línea el entrelazamiento nunca muere, pero en un plano sí existe un punto de no retorno, corrigiendo así teorías matemáticas que antes parecían correctas pero que fallaban por no tener "lente de aumento" suficiente.

La moraleja: A veces, para ver la verdad de la naturaleza, no necesitas nuevas teorías, sino simplemente más poder de cómputo para mirar el problema en su verdadera escala.

Resumen técnico

1. El Problema

La dinámica del entrelazamiento en sistemas cuánticos de muchos cuerpos sometidos a mediciones continuas es un área de investigación activa, centrada en la existencia de Transiciones de Fase Inducidas por Medición (MIPT, por sus siglas en inglés).

• Contexto: En sistemas de fermiones libres (no interactuantes) con simetría U(1), los resultados teóricos y numéricos previos han sido contradictorios.
• La Controversia:
  • La teoría del Modelo Sigma No Lineal (NLSM) predice que para fermiones de Dirac en 1D, no existe MIPT (el entrelazamiento sigue una ley de área para cualquier tasa de monitoreo). Sin embargo, para fermiones de Majorana o en 2D, el NLSM predice una transición.
  • Estudios numéricos anteriores reportaron resultados mixtos: algunos confirmaron la ausencia de transición en 1D, mientras que otros (usando protocolos de "salto cuántico") reportaron una transición.
• Limitación Crítica: La mayoría de los estudios numéricos previos se limitaron a tamaños de red pequeños ($L \lesssim 1000$ en 1D y $\sim 60 \times 60$ en 2D). Dado que la longitud de correlación en la fase de ley de área crece exponencialmente con la inversa de la tasa de monitoreo ($\gamma$), los sistemas pequeños no pueden distinguir entre una caída exponencial real (ley de área) y una caída de ley de potencias (ley de volumen) causada por efectos de tamaño finito.

2. Metodología

Los autores abordaron el problema mediante simulaciones numéricas de alta precisión utilizando Unidades de Procesamiento Gráfico (GPU) para superar las limitaciones de tamaño de sistema.

• Modelo: Fermiones de Dirac no interactuantes en 1D y 2D con simetría U(1).
  • Hamiltoniano: Cadena de hopping con condiciones de frontera periódicas y llenado de medio electrón.
  • Protocolos de Medición:
    1. Medición Projectiva (PM): Mediciones fuertes y aleatorias de la ocupación de sitios individuales.
    2. Difusión de Estado Cuántico (QSD): Mediciones débiles y continuas en todos los sitios, gobernadas por una ecuación de Schrödinger estocástica.
• Técnica Computacional:
  • Implementación masiva en GPUs (tarjetas NVIDIA A100).
  • Se aprovecharon las propiedades cuadráticas del sistema (teorema de Wick) para trabajar exclusivamente con la matriz de correlación ($L \times L$) en lugar de la función de onda completa, lo que reduce la complejidad computacional.
  • Escalabilidad: Se alcanzaron tamaños de red sin precedentes:
    • 1D: Hasta $L = 16,384$ (comparado con $L \lesssim 2000$ en estudios anteriores).
    • 2D: Hasta $160 \times 160$ (comparado con $\sim 60 \times 60$).
• Observables:
  • Se calculó la Entropía de Entrelazamiento (EE) a través de la función de correlación de densidad $C(r, t)$.
  • En 2D, se utilizó la Información Mutua ($I_2$) y la covarianza de número de partículas ($G_{AB}$) como indicadores de la transición, ya que son invariantes de escala en el punto crítico.

3. Contribuciones Clave

• Superación de la Barrera de Tamaño: Demostraron que el uso de GPUs permite acceder a regímenes de tamaño de sistema donde los efectos de borde desaparecen, resolviendo la ambigüedad de estudios anteriores.
• Validación Cuantitativa: Proporcionaron una descripción cuantitativa completa de la dinámica de entrelazamiento, validando o refutando predicciones analíticas del NLSM con datos numéricos de alta precisión.
• Unificación de Protocolos: Compararon sistemáticamente los protocolos de medición projectiva y difusión, mostrando que la física subyacente es robusta independientemente del método de medición, aunque los valores críticos difieren.

4. Resultados Principales

A. En 1 Dimensión (1D)

• Conclusión: No existe transición de fase inducida por medición (MIPT) para fermiones de Dirac con simetría U(1), ni para el protocolo PM ni para QSD.
• Evidencia:
  • Para confirmar la ley de área (decaimiento exponencial de $C(r)$), es necesario alcanzar tamaños de red $L \sim 10,000$.
  • En sistemas pequeños ($L < 4000$), la longitud de correlación es tan grande que el sistema parece estar en una fase de ley de volumen (decaimiento de ley de potencias), lo que llevó a falsos positivos en estudios anteriores.
  • Al ajustar la longitud de correlación ($l_{cor}$) con la tasa de monitoreo ($\gamma$), se encuentra que $l_{cor}$ crece exponencialmente y la tasa crítica $\gamma_c$ es consistentemente cero ($\gamma_c = 0.0 \pm 0.1$).

B. En 2 Dimensiones (2D)

• Conclusión: Sí existe una MIPT a una tasa de monitoreo crítica finita ($\gamma_c$).
• Características de la Transición:
  • La transición se caracteriza por una información mutua invariante de escala (independiente del tamaño del sistema $L$) en el punto crítico.
  • Valores Críticos:
    • Protocolo QSD: $\gamma_c \approx 4.77 \pm 0.01$.
    • Protocolo PM: $\gamma_c \approx 5.72 \pm 0.02$.
  • Exponente Crítico: El exponente crítico $\nu$ que gobierna la divergencia de la longitud de correlación es $\nu \approx 1.3$ para ambos protocolos.
• Implicación Teórica:
  • El valor $\nu \approx 1.3$ no coincide con la predicción analítica del NLSM ($\nu = 1$) ni con la clase de universalidad BDI de la transición de Anderson ($\nu \approx 1.06$).
  • El resultado coincide con la transición de Anderson no hermitiana en 3D de la clase de universalidad AI†, sugiriendo que la simetría efectiva del Hamiltoniano no hermitiano es la que controla la física crítica.

5. Significado e Impacto

• Resolución de Controversias: El trabajo aclara definitivamente que la ausencia de MIPT en 1D para fermiones libres es un hecho robusto, pero solo es visible en el límite termodinámico realizable gracias a la potencia de cómputo de las GPUs.
• Precisión Cuantitativa: Establece que las predicciones analíticas basadas en NLSM son cualitativamente correctas (existencia de transición en 2D, ausencia en 1D) pero cuantitativamente incorrectas en cuanto a los exponentes críticos y la ubicación exacta de la transición.
• Nueva Clase de Universalidad: Sugiere que la física de las transiciones inducidas por medición en sistemas fermiónicos monitoreados pertenece a una clase de universalidad diferente (posiblemente AI†) a la de los sistemas de Anderson tradicionales, lo que requiere un replanteamiento de los modelos teóricos efectivos.
• Metodología: Demuestra que la simulación de sistemas cuánticos de muchos cuerpos en GPUs es esencial para estudiar fenómenos de fase donde las longitudes de correlación divergen exponencialmente, abriendo la puerta a estudios cuantitativos precisos en sistemas de monitoreo más complejos.