⚛️ física cuántica

Realistic simulation of universal multiport interferometers with losses within a Gaussian formalism

Este artículo propone un formalismo gaussiano para simular de manera realista interferómetros multiport de uso universal con pérdidas ópticas no uniformes, permitiendo la caracterización y predicción precisa de los resultados de algoritmos cuánticos en plataformas fotónicas NISQ ruidosas.

Autores originales: Sebastiano Corli, Dario Chemoli, Samuele Altilia, Edoardo Suerra, Dario Tamascelli, Simone Cialdi, Enrico Prati

Publicado 2026-10-14

📖 6 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Sebastiano Corli, Dario Chemoli, Samuele Altilia, Edoardo Suerra, Dario Tamascelli, Simone Cialdi, Enrico Prati

Artículo original bajo licencia CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Esta es una explicación generada por IA del artículo a continuación. No ha sido escrita ni avalada por los autores. Para mayor precisión técnica, consulte el artículo original. Leer descargo de responsabilidad completo

En la carrera por construir una máquina que pueda resolver problemas fuera del alcance de las computadoras actuales, la luz ha emergido como un contendiente principal. A diferencia de otros enfoques que dependen de diminutos circuitos electrónicos o átomos atrapados, estas computadoras fotónicas utilizan partículas individuales de luz, llamadas fotones, para transportar información. La ventaja reside en cómo se comporta la luz de forma natural; viaja rápido, genera poco calor y puede manipularse con espejos y lentes para realizar cálculos complejos. Sin embargo, este camino está plagado de un obstáculo único y persistente: la pérdida. A medida que los fotones viajan a través de los delicados canales de vidrio de un chip informático, algunos inevitablemente desaparecen, absorbidos por el material o dispersos en el aire. Esta desaparición no es solo un fallo menor; es la barrera principal que impide que estas máquinas alcancen todo su potencial. Para entender qué tan cerca estamos de un dispositivo funcional, los científicos deben ser capaces de simular estas máquinas en computadoras clásicas, pero hacerlo es increíblemente difícil cuando el sistema es imperfecto. Si una simulación asume que la luz se comporta perfectamente, puede predecir un resultado que la máquina real y defectuosa nunca podrá producir.

Un equipo de investigadores ha desarrollado una nueva forma de modelar estas máquinas imperfectas, creando un plano realista de cómo se comporta la luz cuando se escapa del sistema. Su trabajo se centra en un tipo específico de computadora óptica que utiliza una red de interferómetros diminutos —dispositivos que dividen y recombinan haces de luz— para procesar información. En un mundo ideal, estos dispositivos guiarían cada fotón exactamente hacia donde necesita ir. En la realidad, el vidrio y las conexiones entre ellos son imperfectos, lo que provoca que una fracción de la luz se pierda en cada paso. Los investigadores se dieron cuenta de que intentar rastrear cada fotón perdido y el entorno al que escapa haría que la simulación por computadora fuera demasiado pesada para ejecutarse. En su lugar, idearon un método para condensar todas estas pérdidas en una forma matemática manejable. Al tratar toda la red de caminos de luz como un único sistema interactuando con un entorno simplificado, pudieron calcular el estado final de la luz sin necesidad de simular las infinitas posibilidades de la energía perdida.

El núcleo de su descubrimiento es una técnica que les permite tomar una red compleja y desordenada de caminos de luz y describirla mediante un mapa único y unificado de cómo cambia la luz. Demostraron que incluso cuando la luz se pierde de manera desigual —es decir, cuando algunos caminos son más sucios o están más dañados que otros—, el efecto general puede capturarse mediante un conjunto específico de números que describen la "imprecisión" o incertidumbre del sistema. Este enfoque funciona para dos formas distintas de construir estas computadoras: una donde la luz fluye a través de muchos cables paralelos a la vez, y otra donde la luz viaja en bucles, rebotando de un lado a otro. Los investigadores probaron su método ejecutando miles de simulaciones en una computadora clásica, imitando el comportamiento de estos dispositivos ópticos bajo diversas condiciones. Encontraron que, si bien el comportamiento de la luz se ve afectado por la disposición específica de la máquina, el nuevo modelo podía predecir con precisión el resultado del experimento, incluyendo la probabilidad de detectar un número específico de fotones al final.

Uno de los hallazgos más significativos de este estudio desafía un atajo común utilizado por los científicos. Durante años, muchos investigadores han asumido que las pérdidas ópticas son uniformes, es decir, que ocurren exactamente al mismo ritmo en todas las partes de la máquina, para facilitar sus cálculos. Las nuevas simulaciones muestran que esta suposición suele ser incorrecta. En uno de los diseños comunes, donde la luz viaja a través de una malla triangular de componentes, las pérdidas no son uniformes; varían dependiendo de dónde se encuentre la luz en la red. Los investigadores demostraron que asumir pérdidas uniformes puede conducir a predicciones erróneas sobre la salida de la máquina, especialmente a medida que la máquina se hace más grande. En contraste, su nuevo modelo, que tiene en cuenta la naturaleza específica y desigual de las pérdidas, proporciona una imagen mucho más precisa de lo que la máquina realmente hará. Esto es crucial porque si los científicos quieren construir sistemas de corrección de errores para solucionar estos fallos, primero deben saber exactamente cómo son esos fallos.

El equipo también exploró cómo estas pérdidas afectan el resultado final de un cálculo conocido como muestreo de bosones gaussianos, una tarea diseñada para demostrar que las computadoras cuánticas pueden superar a las clásicas. Simularon el proceso de enviar luz comprimida —un estado especial donde la incertidumbre de la luz se reduce en una dirección— a través de una red con pérdidas. Observaron que, a medida que la luz viaja a través de la máquina imperfecta, la probabilidad de detectar fotones cambia drásticamente. El estado de vacío, donde no se detectan fotones, se vuelve mucho más común en un entorno ruidoso que en uno perfecto. Además, descubrieron que el patrón específico de pérdidas cambia la probabilidad de detectar ciertos números de fotones, rompiendo la simetría que existiría en una máquina perfecta. Al comparar su modelo realista contra la suposición simplificada de pérdida uniforme, pudieron ver exactamente dónde divergen ambos. La diferencia aumenta a medida que aumenta el número de caminos de luz y a medida que la pérdida se vuelve mayor, confirmando que la visión simplificada es insuficiente para dispositivos a gran escala.

Este trabajo proporciona una herramienta vital para el futuro de la computación cuántica óptica. Ofrece una forma de tomar datos del mundo real sobre cómo un chip específico pierde luz e integrarlos directamente en una simulación para ver cómo se desempeñaría un algoritmo cuántico en ese hardware exacto. Esto significa que, en lugar de adivinar cómo se comportará una máquina, los ingenieros ahora pueden probar sus diseños con un alto grado de realismo incluso antes de que sean construidos. Los investigadores sugieren que este marco puede adaptarse a cualquier tipo de arquitectura óptica de variables continuas, ya sea que utilice trayectorias espaciales o bucles basados en el tiempo. Al caracterizar con precisión estas pérdidas, la comunidad científica puede acercarse al desarrollo de técnicas efectivas de corrección de errores, las cuales son esenciales para construir una computadora cuántica tolerante a fallos. El estudio no pretende haber resuelto el problema de la pérdida de luz, pero ha proporcionado un mapa claro y realista del terreno, mostrando exactamente dónde se encuentran los obstáculos y cómo navegar por ellos. Esta claridad es un paso necesario para convertir la promesa teórica de la computación cuántica fotónica en una realidad práctica.

Título: Simulación realista de interferómetros multiport universales con pérdidas dentro de un formalismo Gaussiano

Planteamiento del problema
Las tecnologías fotónicas son contendientes principales para demostrar la ventaja cuántica, pero las pérdidas de fotones representan el obstáculo primario para la realización de computadoras cuánticas ópticas tolerantes a fallos. La caracterización precisa de estas pérdidas es crítica para las simulaciones clásicas de algoritmos cuánticos en backends ruidosos. Los enfoques de simulación actuales suelen depender de la hipótesis a priori de que las pérdidas ópticas son uniformes en todos los canales. Los autores argumentan que esta suposición puede conducir a predicciones incorrectas sobre los resultados de los algoritmos en hardware específico. Además, aunque existen códigos de corrección de errores para dispositivos basados en qubits, no existe una técnica establecida para corregir las pérdidas de fotones en la computación de variables continuas (CV). Las bibliotecas de software existentes para simulaciones de CV en regímenes ruidosos son insuficientes en comparación con las de qubits basados en materia (por ejemplo, Qiskit, Cirq). Modelar las pérdidas mediante la introducción de modos ambientales externos para cada interacción conduce a un aumento polinomial en el número de modos, lo que hace que las simulaciones del espacio de Fock sean computacionalmente inviables.

Metodología
Los autores proponen un formalismo gaussiano para caracterizar y simular las pérdidas ópticas en interferómetros multiport, aplicable tanto a arquitecturas de codificación de compartimentos espaciales (spatial-bin) como de compartimentos temporales (time-bin). La metodología procede de la siguiente manera:

Marco Teórico: El sistema se describe utilizando la base de Heisenberg (cuadraturas $\hat{q}$ y $\hat{p}$ ) en lugar de la base de Fock. La dinámica de los estados gaussianos está gobernada por el vector de media y la matriz de covarianza. Las pérdidas ópticas se modelan como interacciones entre modos físicos y modos ambientales externos, representadas como operaciones de divisor de haz (beam-splitter).
Caracterización de Interferómetros Imperfectos: Los autores se basan en la caracterización de interferómetros Mach-Zehnder (MZI) imperfectos descrita por Qi et al. [52]. Modelan un interferómetro con pérdidas $T_\eta$ incorporando factores de pérdida en la entrada, la salida y dentro del propio divisor de haz.
Descomposición en Valores Singulares (SVD): Para evitar la explosión computacional de rastrear numerosos modos ambientales externos, los autores emplean SVD. El operador no unitario $A$ resultante de una malla de interferómetros imperfectos se descompone como $A = U_2 \sqrt{H} U_1$ , donde $U_1$ y $U_2$ son matrices unitarias y $\sqrt{H}$ es una matriz diagonal que contiene los factores disipativos (valores singulares). Esto condensa todos los canales de pérdida en $M$ modos externos (donde $M$ es el número de modos físicos), reduciendo el tamaño del sistema de un aumento polinomial a una matriz de covarianza manejable de $2M \times 2M$ .
Derivación de la Matriz de Covarianza: Los autores derivan una expresión analítica para la matriz de covarianza $\sigma_S$ de un proceso de Muestreo de Bosones Gaussianos (GBS) sujeto a estas pérdidas. Esta matriz se obtiene trazando parcialmente fuera (partially tracing out) los modos ambientales de la transformación simpléctica global.
Simulación Numérica: Los autores implementan un flujo de trabajo utilizando un backend gaussiano (biblioteca de Python PNN) para simular GBS. Comparan:
- Resultados ideales de GBS.
- Resultados realistas de GBS utilizando la matriz de pérdida no uniforme $\sqrt{H}$ derivada.
- Resultados de GBS asumiendo pérdidas uniformes (donde $\sqrt{H} \propto I$ ).
  Prueban los esquemas de descomposición de Clements y Reck para la malla del interferómetro.

Contribuciones Clave

Formalismo Analítico: El artículo proporciona una expresión de forma cerrada (Ecuación 24) para la matriz de covarianza de un proceso GBS realizado con interferómetros imperfectos, derivada mediante SVD. Esto permite la simulación de GBS ruidoso y algoritmos generales de tipo dyne en backends clásicos sin recurrir al espacio de Fock.
Reducción Eficiente de Modos: Al aprovechar la SVD, el método reduce el número de modos auxiliares requeridos para modelar las pérdidas de una dependencia polinomial de la profundidad del circuito a una dependencia lineal ( $2M$ ), disminuyendo significativamente los costes computacionales.
Validación de la Suposición de Uniformidad: El trabajo evalúa numéricamente la validez de asumir pérdidas ópticas uniformes. Demuestra que esta suposición no es universalmente cierta y depende fuertemente del esquema de descomposición del interferómetro (Clements vs. Reck) y del número de modos.
Caracterización Específica del Dispositivo: El marco permite la integración de datos experimentales relativos a imperfecciones específicas del hardware (por ejemplo, tasas de pérdida específicas para MZIs individuales) directamente en el flujo de simulación.

Resultos

Desplazamientos de Distribución: Las simulaciones numéricas muestran que las pérdidas ópticas aumentan la población del estado de vacío y permiten la detección de conteos de fotones impares (que están prohibidos en entradas de estado comprimido o squeezed ideales), alterando significativamente la distribución de probabilidad en comparación con el caso ideal.
Dependencia de la Descomposición: El estudio encuentra que la desviación de la matriz de pérdida $\sqrt{H}$ $H$ respecto a la matriz identidad (es decir, la no uniformidad de las pérdidas) se comporta de manera diferente según la arquitectura:
- En el esquema Reck (equivalente a la codificación de tiempo en compartimentos basada en bucles), la distancia de la identidad aumenta con el número de modos $M$ . Los autores ajustan esta tendencia a una ley de potencia y a una función logarítmica, encontrando que los ajustes de ley de potencia son ligeramente mejores.
- En el esquema Clements (codificación espacial), la distancia oscila y converge a un rango acotado después de $M=5$ , lo que sugiere una desviación más estable de la uniformidad.
Divergencia del Hafnian: Al simular el muestreo del Hafnian de una matriz, los autores observan que los resultados de un modelo de pérdida no uniforme realista divergen más rápidamente del caso ideal que la aproximación de pérdida uniforme. Esta discrepancia se vuelve más pronunciada a medida que aumentan las pérdidas ópticas (menor transmisividad $\eta$ ) y a medida que crece el número de modos.

Significancia y Reivindicaciones
Los autores afirman que este trabajo proporciona un "plano" (blueprint) para futuras caracterizaciones de dispositivos y arquitecturas específicas. El marco propuesto permite simulaciones realistas de algoritmos cuánticos en plataformas ópticas de la era NISQ (Escala Intermedia con Ruido). Al modelar con precisión las pérdidas no uniformes, el método tiene como objetivo prevenir predicciones erróneas sobre los resultados de los algoritmos que podrían surgir de asumir distribuciones de pérdida uniformes. Los autores posicionan este trabajo como un paso fundamental hacia la integración de modelos de ruido realistas con técnicas de mitigación de errores y el posible desarrollo de futuros códigos de corrección de errores para la computación de variables continuas. Declaran explícitamente que el trabajo no propone nuevos códigos de corrección de errores, sino una herramienta de simulación para comprender mejor y caracterizar el panorama del ruido.

Más como este