⚛️ quantum physics

qiskit-qudits: A Qiskit Extension for Simulating Qudit Circuits

Este artículo presenta **qiskit-qudits**, una extensión de Qiskit que simula qudits arbitrarios de nivel d codificándolos en qubits con relleno de identidad para manejar dimensiones que no son potencias de dos, ofreciendo tanto descomposiciones de compuertas estándar como síntesis unitaria densa exacta para dimensiones de hasta 16, al tiempo que valida su implementación mediante la verificación numérica de la transformada cuántica de Fourier y las descomposiciones de compuertas.

Autores originales: Francesco Pudda

Publicado 2026-08-25

📖 6 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Francesco Pudda

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 el mundo de la computación cuántica, el bloque de construcción estándar es el cúbit, un sistema diminuto que puede existir en dos estados a la vez, a menudo descritos como cero y uno. Pero la naturaleza no se limita a solo dos opciones. Muchos sistemas físicos, como ciertos átomos o iones, poseen naturalmente tres, cuatro o incluso más niveles de energía distintos. Los científicos llaman a estos sistemas multinivel "qudits". Aunque los cúbits son la base de las computadoras cuánticas actuales, los qudits ofrecen un panorama más rico para la computación, permitiendo potencialmente algoritmos más eficientes y un mejor uso del hardware físico. Sin embargo, un obstáculo significativo ha persistido: las herramientas de software más populares para diseñar programas cuánticos fueron construidas exclusivamente para cúbits de dos niveles. Esto significaba que los investigadores que deseaban experimentar con qudits debían escribir su propio software complejo desde cero o abandonar las poderosas herramientas que la comunidad en general había desarrollado.

Francesco Pudda ha cerrado esta brecha con una nueva extensión de software llamada qiskit-qudits. Esta herramienta permite a los científicos diseñar circuitos cuánticos utilizando qudits multinivel mientras siguen utilizando la suite de software estándar Qiskit, la cual es ampliamente utilizada para construir y probar programas cuánticos. El logro central de este trabajo es un método de traducción ingenioso. Dado que el hardware subyacente y el software estándar piensan en términos de cúbits de dos niveles, el nuevo sistema codifica cada qudit multinivel en un pequeño grupo de cúbits. Por ejemplo, un sistema de tres niveles se representa mediante dos cúbits, y un sistema de cinco niveles mediante tres. Este enfoque crea una ligera ineficiencia porque el grupo de cúbits puede representar más estados de los que el qudit realmente utiliza, pero el software está diseñado para ignorar estos estados adicionales no utilizados. Al tratar estos estados no utilizados como si no hicieran nada, el sistema asegura que la lógica compleja del qudit funcione perfectamente dentro del lenguaje familiar de los cúbits.

Los investigadores no solo construyeron una idea teórica; crearon una biblioteca plenamente funcional que se integra perfectamente con los flujos de trabajo cuánticos existentes. El software define un nuevo conjunto de instrucciones, o puertas, que actúan sobre estos sistemas multinivel. Estas puertas pueden realizar operaciones estándar como desplazar niveles o crear superposiciones complejas. Una característica clave del diseño es su capacidad para manejar diferentes tamaños de qudits simultáneamente. Un solo circuito puede contener un sistema de tres niveles, un sistema de cuatro niveles y un sistema de cinco niveles trabajando juntos, y el software gestiona la traducción para cada uno de forma automática. Cuando el número de niveles es una potencia de dos, el software descompone la operación en una serie de pasos de cúbits estándar que las herramientas de optimización de la computadora pueden procesar fácilmente. Cuando el número de niveles no es una potencia de dos, el software utiliza un método diferente, calculando la transformación matemática exacta para todo el grupo de cúbits a la vez. Esto asegura que el sistema funcione correctamente para cualquier tamaño hasta dieciséis niveles, un rango que cubre muchos sistemas físicos prácticos.

Para demostrar que su método de traducción era preciso, el equipo realizó rigurosas pruebas numéricas. Compararon el comportamiento de su software contra las reglas matemáticas conocidas para un tipo específico de operación cuántica llamada transformada de Fourier cuántica, que es una herramienta fundamental en muchos algoritmos cuánticos. Probaron esto a través de varios tamaños, incluyendo aquellos que no encajan ordenadamente en potencias de dos. Los resultados mostraron que la salida del software coincidía con las expectativas teóricas con una precisión extrema, con errores tan pequeños que eran efectivamente cero para propósitos prácticos. También verificaron que cada tipo de puerta que crearon se comportaba exactamente como pretendían, ya fuera un simple desplazamiento o una interacción multinivel compleja. Las pruebas confirmaron que el software maneja correctamente las partes "no utilizadas" del espacio codificado, asegurando que ninguna información se filtre o se corrompa durante el proceso.

El artículo también destaca las compensaciones involucradas en este enfoque. Debido a que el sistema utiliza cúbits adicionales para representar cada qudit, requiere más memoria para simularse en una computadora clásica de lo que requeriría un simulador de qudits nativo. Para sistemas muy grandes, esta sobrecarga se vuelve significativa, y los investigadores señalan que para la simulación pura de sistemas masivos, una herramienta construida específicamente para qudits podría seguir siendo más eficiente. Sin embargo, para la gran mayoría de los casos de uso actuales, el beneficio de utilizar el ecosistema estándar de Qiskit supera este costo. El software permite a los investigadores dibujar circuitos al nivel del qudit, ejecutarlos a través de pasadas de optimización estándar y prepararlos para su ejecución en hardware cuántico real o simuladores sin necesidad de reescribir el código manualmente.

Uno de los aspectos más prácticos de este trabajo es cómo maneja la medición. Cuando se mide un sistema cuántico, el resultado debe convertirse de vuelta de los datos brutos de cúbits a los niveles significativos del qudit. El nuevo software incluye una capa de decodificación que realiza esta traducción automáticamente. También incluye un mecanismo de seguridad para detectar la "fuga" (leakage), una situación en la que un sistema termina accidentalmente en un estado que no corresponde a ningún nivel válido del qudit. En una simulación perfecta, esto nunca debería suceder, pero en el hardware real o cuando se introduce ruido, puede ocurrir. El software permite a los usuarios elegir cómo manejar estos eventos raros, ya sea descartando los datos erróneos o marcándolos para su revisión, haciendo que los resultados de la simulación sean más fiables y transparentes.

El trabajo representa un paso significativo hacia la realización de que la computación cuántica sea más accesible y versátil. Al convertir los conceptos teóricos de los sistemas cuánticos multinivel en una capa de software funcional y transducible, elimina una barrera importante para los investigadores. El equipo verificó que su método funciona para dimensiones que van desde dos hasta dieciséis, cubriendo las necesidades de muchas plataformas experimentales actuales. Si bien el software depende de una estrategia de codificación específica que añade una capa de complejidad, el autor demuestra que esta complejidad se gestiona de manera efectiva, permitiendo que las herramientas estándar de la comunidad cuántica operen en una gama mucho más amplia de sistemas físicos. El resultado es una herramienta que permite a los científicos explorar todo el potencial de la mecánica cuántica multinivel sin tener que abandonar la infraestructura poderosa y establecida en la que ya confían.

Resumen Técnico: QISKIT-QUDITS

Planteamiento del Problema

Los qudits, que generalizan los qubits de dos niveles a $d$ niveles computacionales, ofrecen ventajas potenciales en el procesamiento de información cuántica. Sin embargo, el marco de trabajo Qiskit carece de soporte nativo para qudits, tratándolos como una brecha de extensión de terceros conocida. Las herramientas existentes enfrentan una dicotomía: los simuladores de qudits nativos (por ejemplo, Cirq, QuTiP) no pueden aprovechar la extensa pila de compilación, las pasadas del transpilador y los backends de hardware de Qiskit, mientras que las herramientas orientadas a qubits carecen de una superficie de programación para operaciones a nivel de qudit. Actualmente, los usuarios deben implementar manualmente esquemas de codificación, matrices de puertas, descomposiciones y decodificación clásica, un proceso propenso a errores con respecto al tratamiento del espacio de estado no físico y la decodificación de la medición.

Metodología

El artículo presenta qiskit-qudits, una extensión de Qiskit que simula qudits de $d$ niveles codificando cada uno en $m = \lceil \log_2 d \rceil$ qubits. El software se construye sobre las clases públicas Instruction, Gate y ControlledGate de Qiskit, asegurando la compatibilidad con el ecosistema existente.

Decisiones de Diseño Centrales

Codificación y Relleno de Identidad (Identity Padding):
Dado que $d$ no es necesariamente una potencia de dos, el espacio de Hilbert codificado ( $D = 2^N$ ) es generalmente mayor que el subespacio lógico ( $L = \prod d_i$ ). La librería adopta la convención de relleno de identidad: cada puerta actúa como la identidad en el complemento ortogonal del subespacio lógico. Formalmente, para una puerta lógica $G$ , la unitaria codificada es $U = E G E^\dagger + (I - \Pi)$ , donde $E$ es la isometría de codificación y $\Pi$ es el proyector sobre el subespacio lógico. Esta elección evita la fuga (leakage) durante la simulación sin ruido y preserva las propiedades de permutación.
Jerarquía de Puertas y Descomposición:
La librería define tres clases base: QuditGate (qudit único), QuditControlledGate (multicontrol) y QuditMultiGate.
- Ruta de Potencia de Dos: Cuando todas las dimensiones de los operandos son potencias de dos, las puertas se descomponen en cascadas fijas de puertas de qubits estándar de Qiskit (por ejemplo, mcx, cp, synth_qft_full). Esto permite que el transpilador de Qiskit las optimice.
- Respaldo de Síntesis Densa (Dense-Synthesis Fallback): Para dimensiones que no son potencias de dos, se recurre a la síntesis unitaria densa exacta. La librería genera la matriz unitaria completa $D \times D$ (vía el protocolo __array__ de NumPy) y la incrusta como una instrucción unitaria opaca única, que la pasada UnitarySynthesis luego descompone.
Modelo de Circuito y Medición:
El QuditQuantumCircuit mantiene dos vistas sincronizadas: un registro de instrucciones a nivel de qudit (.data) y un circuito codificado a nivel de qubit (.circuit). La medición se gestiona mediante QuditMeasure, que mapea los qubits de codificación a bits clásicos (ClByte). Una capa de decodificación dedicada traduce los conteos del backend de vuelta a niveles de qudit, detectando explícitamente eventos de "fuga" (resultados $\ge d$ ), los cuales pueden manejarse mediante políticas configurables (keep, drop, raise).

Contribuciones Clave

Implementación de una Capa de Qudit Transpilable: El trabajo convierte las definiciones teóricas de puertas del trabajo previo [1] en una extensión funcional y transpilable de Qiskit. Soporta dimensiones $2 \le d \le 16$ y registros de dimensiones mixtas.
Simulación Exacta para Dimensiones que no son Potencias de Dos: A diferencia de los enfoques previos de simulación de qubits limitados a $d=2^m$ , esta librería soporta síntesis exacta para $d$ arbitrarios dentro del rango soportado utilizando respaldos unitarios densos.
Verificación de Teoría e Implementación:
- Verificación A: Se verificó que las matrices unitarias densas reproducen la Transformada Discreta de Fourier (DFT) en el subespacio lógico tanto para dimensiones de potencia de dos como para las que no lo son, extendiendo la verificación de [1].
- Verificación B: Se validó que los circuitos descompuestos (tanto las rutas de cascada como las densas) coinciden con las definiciones unitarias densas en 138 instancias de puertas, con desviaciones dentro de la precisión de la máquina ( $\approx 10^{-15}$ ).
- Invariante de Relleno (Padding Invariant): Se confirmó que el invariante de relleno de identidad se mantiene exactamente, asegurando que no haya mezcla entre los subespacios válidos e inválidos.
Detección de Fuga (Leakage): Proporciona un mecanismo explícito para detectar y manejar eventos de fuga en las simulaciones, distinguiendo entre el comportamiento ideal y los errores inducidos por ruido.

Resultados

La verificación numérica confirma la corrección de la implementación:

QFT vs. DFT: La implementación de QFT coincide con la matriz DFT teórica con un margen de $10^{-15}$ para dimensiones $d \in \{2, 3, 4, 5, 8\}$ .
Descomposición de Puertas: Todas las puertas, incluyendo SUMX, SUMP, SWAP y QFT, producen descomposiciones consistentes con sus definiciones unitarias densas.
Análisis de Recursos:
- Potencia de Dos: Las descomposiciones son eficientes y amigables para el transpilador. Por ejemplo, una QFT de 3 qutrits ( $d=3$ ) requiere 1,783 puertas CX, mientras que una QFT de 3 ququart ( $d=4$ ) requiere solo 56 puertas CX.
- No Potencia de Dos: El respaldo denso incurre en un overhead significativo. Una QFT de 3 qutrits ( $d=3$ ) se sintetiza en 1,783 puertas CX, mientras que una QFT de 3 ququint ( $d=5$ ) requiere 119,383 puertas CX. Esto resalta el costo de la síntesis densa para dimensiones más grandes que no son potencias de dos.
- Restricciones de Memoria: El camino denso está limitado por la memoria; una puerta en ocho qudits de 16 niveles ( $N=32$ ) requeriría $\approx 256$ EiB de memoria, lo que hace que el límite práctico sea $N \approx 13\text{--}14$ en hardware convencional.

Significado y Limitaciones

El artículo posiciona a qiskit-qudits como un puente que permite a los usuarios escribir algoritmos de qudits a nivel lógico mientras retienen el acceso total a la pila de compilación, backends y primitivas de Qiskit. Aborda explícitamente la "parte no física" del espacio codificado y la decodificación de la medición, que suelen ser propensas a errores en implementaciones manuales.

Limitaciones y Compensaciones:

Límite de Dimensión: La librería impone un límite de $2 \le d \le 16$ y un máximo de 8 qudits por puerta. Estos son límites nominales; el límite práctico está determinado por la memoria requerida para la síntesis unitaria densa ( $N \lesssim 13\text{--}14$ ).
Semántica de Fuga: La convención de relleno de identidad congela la amplitud filtrada en lugar de modelar la fuga física. No conmuta con el producto tensorial, lo que significa que la equivalencia del circuito debe verificarse en el subespacio proyectado ( $\Pi U \Pi$ ).
Optimización: La implementación actual no explota las filas/columnas de "no importa" (don't-care) en el transpilador, y el respaldo denso para dimensiones que no son potencias de dos es computacionalmente costoso en comparación con los simuladores de qudits nativos o la síntesis de lógica reversible optimizada.

El trabajo concluye que, si bien el respaldo denso permite la simulación exacta para dimensiones que no son potencias de dos, las mejoras futuras deberían centrarse en la síntesis de lógica reversible agnóstica a la dimensión y en exponer el subespacio no físico al transpilador para reducir el overhead de recursos.