Superposición cuántica: estados, interferencia y medición

La mecánica cuántica se fundamenta axiomáticamente en cuatro postulados que relacionan los objetos matemáticos del formalismo con los resultados experimentales. El primero establece que todo sistema físico aislado se describe mediante un espacio de Hilbert \(\mathcal{H}\) —un espacio vectorial completo con producto interno— y que el estado del sistema en un instante dado corresponde a un vector de estado (o ket) \(|\psi\rangle\in\mathcal{H}\) normalizado: \(\langle\psi|\psi\rangle = 1\). Los vectores que difieren únicamente por una fase global \(e^{i\varphi}|\psi\rangle\) representan el mismo estado físico, de modo que los estados puros son en rigor rayos del espacio proyectivo \(\mathbb{P}(\mathcal{H})\).

Para un cúbit, \(\mathcal{H}=\mathbb{C}^2\). La base computacional canónica está formada por los vectores ortonormales \(|0\rangle\) y \(|1\rangle\), que satisfacen \(\langle 0|1\rangle=0\) y \(\langle j|j\rangle=1\). Cualquier estado puro se expresa entonces como una combinación lineal:

El segundo postulado dicta que la evolución temporal de un sistema cerrado es unitaria: existe un operador \(U\) tal que \(U^\dagger U = \mathbb{I}\) y \(|\psi(t)\rangle = U(t,t_0)\,|\psi(t_0)\rangle\). En el caso de evolución continua gobernada por un hamiltoniano independiente del tiempo \(H\), la unitaria es \(U(t)=e^{-iHt/\hbar}\), lo que garantiza la conservación de la norma —y, por tanto, de las probabilidades—. La linealidad de la ecuación de Schrödinger \(i\hbar\,\partial_t|\psi\rangle = H|\psi\rangle\) es precisamente lo que permite que una superposición de soluciones siga siendo solución: si \(|\phi_1\rangle\) y \(|\phi_2\rangle\) son soluciones, también lo es \(\alpha|\phi_1\rangle+\beta|\phi_2\rangle\).

El tercer postulado, relativo a la medición, introduce un quiebre con la evolución unitaria: la observación de un observable \(\hat{A}=\sum_a a\,|a\rangle\langle a|\) colapsa el estado a uno de los autoestados \(|a\rangle\) con probabilidad \(p(a)=|\langle a|\psi\rangle|^2\). Este proceso es irreversible, no unitario y estocástico. Finalmente, el cuarto postulado establece que el espacio de Hilbert de un sistema compuesto es el producto tensorial \(\mathcal{H}_A\otimes\mathcal{H}_B\) de los espacios de sus subsistemas, abriendo la puerta a estados entrelazados que no son factorizables.

La estructura de espacio de Hilbert impone restricciones fundamentales. La desigualdad de Cauchy–Schwarz \(|\langle\phi|\psi\rangle|^2 \le \langle\phi|\phi\rangle\langle\psi|\psi\rangle\) garantiza que las probabilidades de transición entre estados normalizados nunca exceden la unidad. El teorema espectral asegura que todo observable hermítico admite una descomposición en autovalores reales y autovectores ortogonales, lo cual proporciona la base sobre la que se proyecta al medir. Por último, la completitud del espacio —la convergencia de toda sucesión de Cauchy— garantiza que los límites de secuencias de estados permanezcan dentro de \(\mathcal{H}\), una propiedad esencial para definir operadores no acotados como el hamiltoniano de un oscilador armónico.

Puesto que la fase global es irrelevante, un estado puro de un cúbit posee únicamente dos grados de libertad reales. Eligiendo la fase de \(\alpha\) como referencia (haciéndola real y no negativa), la parametrización más natural emplea dos ángulos \(\theta\in[0,\pi]\) y \(\varphi\in[0,2\pi)\):

Esta expresión establece una correspondencia biyectiva entre estados puros del cúbit y puntos de la esfera de Bloch \(S^2\subset\mathbb{R}^3\). El polo norte \((\theta=0)\) representa \(|0\rangle\), el polo sur \((\theta=\pi)\) corresponde a \(|1\rangle\), y el ecuador \((\theta=\pi/2)\) alberga las superposiciones equilibradas con igual peso en ambas componentes, diferenciándose únicamente en la fase relativa \(\varphi\). Los estados canónicos \(|+\rangle=\tfrac{1}{\sqrt{2}}(|0\rangle+|1\rangle)\) y \(|-\rangle=\tfrac{1}{\sqrt{2}}(|0\rangle-|1\rangle)\) se encuentran en los extremos del eje \(x\), mientras que \(|{+i}\rangle\) y \(|{-i}\rangle\) residen sobre el eje \(y\).

Las coordenadas cartesianas del vector de Bloch \(\vec{r}=(r_x, r_y, r_z)\) se obtienen mediante los valores esperados de las matrices de Pauli:

Para estados puros \(|\vec{r}|=1\); para estados mixtos descritos por una matriz de densidad \(\rho=\tfrac{1}{2}(\mathbb{I}+\vec{r}\cdot\vec{\sigma})\), el vector de Bloch satisface \(|\vec{r}|<1\) y reside en el interior de la esfera. El centro \(\vec{r}=\vec{0}\) corresponde al estado máximamente mezclado \(\rho=\mathbb{I}/2\), que no contiene información alguna sobre el cúbit.

La representación de Bloch convierte las puertas unitarias de un cúbit en rotaciones de la esfera. Una rotación de ángulo \(\alpha\) alrededor del eje \(\hat{n}\) se expresa como \(R_{\hat{n}}(\alpha)=e^{-i\alpha\hat{n}\cdot\vec{\sigma}/2}=\cos(\alpha/2)\,\mathbb{I}-i\sin(\alpha/2)\,\hat{n}\cdot\vec{\sigma}\). La puerta Hadamard \(H\), por ejemplo, corresponde a una rotación de \(\pi\) alrededor del eje \((\hat{x}+\hat{z})/\sqrt{2}\), llevando \(|0\rangle\) al estado \(|+\rangle\) sobre el ecuador. Las puertas \(R_z(\varphi)\) son rotaciones alrededor del eje \(z\) que modifican la fase relativa sin alterar las probabilidades de medición en la base computacional, mientras que \(R_x(\theta)\) intercambia peso entre \(|0\rangle\) y \(|1\rangle\).

La esfera de Bloch proporciona además una interpretación geométrica de la fidelidad entre estados puros. Si dos estados puros \(|\psi\rangle\) y \(|\phi\rangle\) tienen vectores de Bloch separados por un ángulo \(\Theta\) en la esfera, la fidelidad es \(F=|\langle\phi|\psi\rangle|^2=\cos^2(\Theta/2)\). Para estados mixtos, la fidelidad generalizada de Uhlmann satisface \(F(\rho,\sigma)=\bigl(\mathrm{tr}\sqrt{\sqrt{\rho}\,\sigma\sqrt{\rho}}\bigr)^2\), una medida que se reduce al cuadrado del producto interno para estados puros y es ampliamente utilizada como figura de mérito en tomografía cuántica y benchmarks de procesadores cuánticos.

Los coeficientes \(\alpha\) y \(\beta\) en la expansión \(|\psi\rangle=\alpha|0\rangle+\beta|1\rangle\) son amplitudes de probabilidad. A diferencia de las probabilidades clásicas —reales y no negativas—, las amplitudes son números complejos que pueden sumarse, restarse e interferir entre sí. La conexión con el mundo experimental se establece mediante la regla de Born (1926):

donde \(|a_k\rangle\) es un autovector del observable medido. Para la base computacional del cúbit, esto se simplifica a \(p(0)=|\alpha|^2\) y \(p(1)=|\beta|^2\). La condición de normalización \(\sum_k p(a_k)=1\) es consecuencia directa de la normalización del estado y de la completitud de la base: \(\sum_k|a_k\rangle\langle a_k|=\mathbb{I}\).

La naturaleza compleja de las amplitudes tiene consecuencias profundas. Escribiendo \(\alpha=|\alpha|e^{i\phi_\alpha}\) y \(\beta=|\beta|e^{i\phi_\beta}\), la fase relativa \(\delta=\phi_\beta-\phi_\alpha\) es inaccesible mediante una medición en la base \(\{|0\rangle,|1\rangle\}\), pero se manifiesta inequívocamente al medir en una base complementaria. Por ejemplo, la probabilidad de obtener el resultado \(|+\rangle\) es:

El término \(2|\alpha||\beta|\cos\delta\) es el término de interferencia: depende exclusivamente de la fase relativa y no aparece en una mezcla clásica de \(|0\rangle\) y \(|1\rangle\) con las mismas probabilidades marginales. Este término puede variar continuamente entre \(-2|\alpha||\beta|\) y \(+2|\alpha||\beta|\), produciendo interferencia destructiva o constructiva según el valor de \(\delta\).

La regla de Born se extiende de manera natural a espacios de dimensión arbitraria. Para un sistema de \(n\) cúbits (\(\dim\mathcal{H}=2^n\)), un estado genérico se escribe \(|\psi\rangle=\sum_{k=0}^{2^n-1}c_k|k\rangle\) con \(\sum_k|c_k|^2=1\), y la probabilidad de observar la cadena de bits \(k\) es \(p(k)=|c_k|^2\). La riqueza computacional de la computación cuántica reside en la posibilidad de manipular \(2^n\) amplitudes simultáneamente mediante puertas unitarias que actúan sobre \(n\) cúbits —un paralelismo de amplitudes que no debe confundirse con el cómputo paralelo clásico—.

En el contexto de la tomografía de estado cuántico, la regla de Born fundamenta el problema inverso: a partir de estadísticas de medición en múltiples bases, se reconstruye la matriz de densidad \(\rho\) del sistema. Para un cúbit, bastan las estadísticas en tres bases mutuamente no sesgadas (por ejemplo \(\sigma_x\), \(\sigma_y\), \(\sigma_z\)) para determinar completamente \(\rho\). En general, para un sistema de \(d\) dimensiones se requieren mediciones en al menos \(d+1\) bases mutuamente no sesgadas (cuando existen) o, alternativamente, \(d^2-1\) parámetros reales independientes, lo cual escala exponencialmente con el número de cúbits y constituye uno de los principales desafíos de la verificación cuántica a gran escala.

La interferencia es la firma operacional de la superposición: cuando un sistema cuántico puede seguir dos o más caminos alternativos hacia un mismo estado final, las amplitudes asociadas a cada camino se suman antes de calcular la probabilidad. Si las amplitudes tienen la misma fase, se produce interferencia constructiva y la probabilidad se amplifica; si difieren en \(\pi\), la interferencia destructiva puede cancelar la probabilidad por completo. Este mecanismo no tiene análogo clásico, donde las probabilidades de caminos alternativos simplemente se suman.

El ejemplo paradigmático es el interferómetro de Mach–Zehnder. Un fotón incide sobre un divisor de haz equilibrado (50/50), que produce la superposición \(|0\rangle\mapsto\tfrac{1}{\sqrt{2}}(|0\rangle+i|1\rangle)\). Los dos caminos \(|0\rangle\) y \(|1\rangle\) acumulan fases distintas \(\phi_0\) y \(\phi_1\) al atravesar ramas ópticas de longitudes diferentes. Un segundo divisor de haz recombina las amplitudes:

La probabilidad de detección en el puerto \(|0\rangle\) es \(p(0)=\sin^2\!\bigl(\tfrac{\phi_0-\phi_1}{2}\bigr)\). Cuando \(\phi_0=\phi_1\) (desfase nulo), la interferencia es perfectamente destructiva en ese puerto y el fotón sale con certeza por \(|1\rangle\). Cuando \(\phi_0-\phi_1=\pi\), la situación se invierte. Variando la diferencia de fase continuamente, se observan las célebres franjas de interferencia sinusoidales.

La visibilidad de las franjas se define como \(\mathcal{V}=\frac{p_{\max}-p_{\min}}{p_{\max}+p_{\min}}\) y cuantifica el grado de coherencia del sistema. Para un estado puro en superposición equilibrada, \(\mathcal{V}=1\). La decoherencia, el ruido de fase y la distinción parcial de caminos reducen \(\mathcal{V}\) por debajo de la unidad. La visibilidad está directamente relacionada con las coherencias fuera de la diagonal de la matriz de densidad: en la base del interferómetro, \(\mathcal{V}=2|\rho_{01}|\).

En el contexto de la computación cuántica, la interferencia es el recurso que los algoritmos explotan para amplificar las amplitudes de las soluciones correctas y suprimir las incorrectas. El algoritmo de Deutsch–Jozsa constituye la demostración más simple: una sola evaluación del oráculo, seguida de puertas Hadamard, provoca que todas las amplitudes de resultados incorrectos se cancelen destructivamente si la función es constante, mientras que se refuercen constructivamente si es balanceada. El algoritmo de Grover generaliza este principio mediante el operador de inversión sobre la media: en cada iteración, la amplitud de la solución marcada crece en \(\sim 1/\sqrt{N}\), logrando una aceleración cuadrática respecto a la búsqueda clásica.

La interferencia cuántica se manifiesta también en el fenómeno de batido cuántico (quantum beating), observable en sistemas de dos niveles acoplados a campos coherentes. Si un cúbit se prepara en el estado \(|+\rangle\) y evoluciona bajo el hamiltoniano \(H=\tfrac{\omega}{2}\sigma_z\), la probabilidad de encontrarlo en \(|+\rangle\) tras un tiempo \(t\) oscila como \(p(+,t)=\cos^2(\omega t/2)\), describiendo las oscilaciones de Rabi cuando se incluye un campo transversal resonante. Estas oscilaciones son una manifestación directa de la evolución coherente de la superposición y su observación es rutinaria en laboratorios de cúbits superconductores, iones atrapados y puntos cuánticos.

La extensión de la superposición a registros de \(n\) cúbits revela una de las propiedades más extraordinarias de la mecánica cuántica: el crecimiento exponencial del espacio de estados. El espacio de Hilbert compuesto es el producto tensorial de los espacios individuales:

Un estado genérico de \(n\) cúbits se escribe como una superposición de \(2^n\) cadenas de bits: \(|\psi\rangle=\sum_{k=0}^{2^n-1}c_k|k\rangle\), donde cada \(|k\rangle\) representa una cadena binaria de \(n\) dígitos. Para \(n=300\), la dimensión del espacio de Hilbert \(2^{300}\approx 10^{90}\) supera el número estimado de partículas en el universo observable, lo cual ilustra por qué la simulación clásica exacta de un sistema cuántico genérico es intratable.

No todos los estados de un sistema compuesto son iguales en su estructura. Un estado producto tiene la forma \(|\psi\rangle=|\psi_A\rangle\otimes|\psi_B\rangle\), lo que implica que los subsistemas \(A\) y \(B\) no poseen correlaciones cuánticas y pueden describirse de manera independiente. Sin embargo, la mayoría de los estados en \(\mathcal{H}_A\otimes\mathcal{H}_B\) son estados entrelazados que no admiten tal factorización. El ejemplo prototípico es el estado de Bell:

En este estado, medir el primer cúbit y obtener \(|0\rangle\) colapsa instantáneamente el segundo a \(|0\rangle\), independientemente de la distancia espacial entre ambos. Las correlaciones resultantes violan las desigualdades de Bell, descartando toda teoría local de variables ocultas y confirmando la naturaleza no-local de la mecánica cuántica, como demostraron los experimentos de Aspect (1982) y los más recientes tests de Bell sin lagunas (2015).

La entropía de entrelazamiento cuantifica la cantidad de entrelazamiento en un estado bipartito puro. Se define como la entropía de von Neumann de la matriz de densidad reducida: \(S(\rho_A)=-\mathrm{tr}(\rho_A\log_2\rho_A)\), donde \(\rho_A=\mathrm{tr}_B(|\psi\rangle\langle\psi|)\). Para un estado producto, \(S=0\); para un estado de Bell, \(S=1\) ebit (bit de entrelazamiento). Esta cantidad es invariante bajo operaciones unitarias locales y constituye el recurso fundamental que los protocolos de teleportación cuántica, destilación de entrelazamiento y computación cuántica basada en mediciones consumen.

En la práctica, preparar y mantener superposiciones de muchos cúbits es un desafío formidable. Cada cúbit adicional duplica la dimensión del espacio de Hilbert y, con ella, el número de modos de decoherencia potenciales. Los registros cuánticos actuales —de decenas a centenares de cúbits— operan en el régimen NISQ (Noisy Intermediate-Scale Quantum), donde la profundidad máxima de circuito está limitada por las tasas de error de las puertas y la decoherencia. La corrección cuántica de errores, que codifica un cúbit lógico en múltiples cúbits físicos, es la estrategia fundamental para escalar hacia computación cuántica tolerante a fallos, pero exige una sobrecarga significativa en recursos físicos.

La medición proyectiva, o medición de von Neumann, constituye el formalismo estándar para describir la extracción de información de un sistema cuántico. Sea \(\hat{A}=\sum_k a_k\,P_k\) un observable con descomposición espectral en proyectores ortogonales \(P_k=|a_k\rangle\langle a_k|\) (asumiendo espectro no degenerado). Al medir \(\hat{A}\) sobre el estado \(|\psi\rangle\), se obtiene el resultado \(a_k\) con probabilidad:

y el estado post-medición colapsa a:

Este proceso es irreversible: la información contenida en las amplitudes relativas respecto a la base de medición se pierde irrecuperablemente. Si el observable tiene espectro degenerado, el proyector \(P_k\) proyecta sobre el subespacio propio asociado al autovalor \(a_k\), y el estado post-medición retiene la superposición dentro de ese subespacio.

La medición proyectiva satisface varias propiedades formales importantes: idempotencia (\(P_k^2=P_k\)), hermiticidad (\(P_k^\dagger=P_k\)), ortogonalidad (\(P_jP_k=\delta_{jk}P_k\)) y completitud (\(\sum_k P_k=\mathbb{I}\)). Estas condiciones garantizan que las probabilidades sumen uno y que una medición repetida inmediatamente después produzca el mismo resultado con certeza (reproducibilidad).

En sistemas compuestos, la medición parcial sobre un subsistema introduce correlaciones clásicas y modifica el estado del subsistema no medido. Si el estado bipartito es \(|\psi\rangle_{AB}=\sum_{j,k}c_{jk}|j\rangle_A|k\rangle_B\) y se mide el subsistema \(A\) obteniendo \(|j_0\rangle\), el estado del subsistema \(B\) colapsa condicionalmente a:

Para un estado de Bell \(|\Phi^+\rangle\), este colapso condicional es total: medir \(|0\rangle\) en \(A\) fuerza a \(B\) a colapsar a \(|0\rangle\) con certeza. Este fenómeno subyace a la teleportación cuántica, el intercambio de entrelazamiento (entanglement swapping) y los protocolos de distribución de claves cuánticas como BB84 y E91.

El problema de la medición sigue siendo una de las cuestiones fundacionales más debatidas en la interpretación de la mecánica cuántica. La ecuación de Schrödinger, que gobierna la evolución unitaria, no contiene ningún mecanismo para el colapso: el postulado de medición se introduce como una regla adicional. Las interpretaciones de Copenhague, de muchos mundos, de historias consistentes y la decoherencia ambiental ofrecen perspectivas distintas sobre esta aparente dualidad. Desde el punto de vista operacional —el relevante para la computación cuántica—, el postulado de medición es una regla empíricamente verificada que describe cómo las estadísticas de los detectores se relacionan con la descripción cuántica del sistema preparado.

En la práctica, ningún sistema cuántico está perfectamente aislado. La interacción con grados de libertad ambientales no controlados produce decoherencia: la pérdida gradual de las relaciones de fase entre las componentes de una superposición, convirtiendo un estado puro coherente en una mezcla estadística clásica. Formalmente, la decoherencia se describe mediante la evolución de la matriz de densidad reducida del sistema, obtenida al trazar sobre los grados de libertad del ambiente:

donde \(U(t)\) es la unitaria del sistema total (sistema + ambiente) y \(\rho_E\) es el estado inicial del ambiente. Esta evolución define un canal cuántico (mapa completamente positivo y que preserva la traza, CPTP), que se expresa mediante la representación de operadores de Kraus:

Para un cúbit, los canales de ruido más relevantes experimentalmente se caracterizan por dos escalas temporales: \(T_1\) (tiempo de relajación longitudinal, asociado a la pérdida de energía) y \(T_2\) (tiempo de coherencia transversal, asociado a la pérdida de fase). En general, \(T_2\le 2T_1\), y la igualdad se alcanza cuando la desfase pura (\(T_\varphi\)) es despreciable: \(1/T_2 = 1/(2T_1) + 1/T_\varphi\).

El canal de amortiguamiento de amplitud modela la relajación espontánea del estado excitado \(|1\rangle\) al fundamental \(|0\rangle\), con operadores de Kraus:

donde \(\gamma=1-e^{-t/T_1}\) es la probabilidad de decaimiento. Este canal describe la emisión espontánea en átomos, la relajación en cúbits superconductores y la pérdida de fotones en cavidades ópticas. A tiempo largo (\(\gamma\to 1\)), cualquier estado converge al fundamental \(|0\rangle\langle 0|\).

El canal de desfase (dephasing) preserva las poblaciones diagonales pero atenúa las coherencias fuera de la diagonal: \(\mathcal{E}_{\mathrm{deph}}(\rho)=(1-p)\rho+p\,\sigma_z\rho\,\sigma_z\), donde \(p=\tfrac{1}{2}(1-e^{-t/T_\varphi})\). Las coherencias \(\rho_{01}\) decaen como \(e^{-t/T_\varphi}\), convirtiendo la superposición en una mezcla clásica de \(|0\rangle\) y \(|1\rangle\) sin modificar las probabilidades de medición en la base computacional.

El canal despolarizante modela un ruido isotrópico que empuja cualquier estado hacia el estado máximamente mezclado:

En la esfera de Bloch, el canal despolarizante contrae uniformemente el vector de Bloch: \(\vec{r}\mapsto(1-4p/3)\vec{r}\). Cuando \(p=3/4\), el estado resultante es siempre \(\mathbb{I}/2\), independientemente del estado inicial. Este canal es particularmente útil como modelo teórico de ruido por su simetría, y aparece como modelo efectivo en escenarios de twirling aleatorio de puertas.

La plataforma Harmoniq materializa los conceptos de superposición desarrollados en las secciones anteriores a través de dos simuladores cuánticos de alto rendimiento: Hopper y Hermann. Ambos implementan la simulación de vector de estado completo, almacenando explícitamente las \(2^n\) amplitudes complejas del registro y aplicando puertas unitarias como multiplicaciones de matrices dispersas o transformaciones indexadas. Hopper está optimizado para circuitos de hasta 28 cúbits con ejecución en CPU multinúcleo, mientras que Hermann aprovecha aceleración en GPU para sistemas de hasta 34 cúbits, permitiendo explorar superposiciones en espacios de Hilbert de dimensión \(\sim 1.7\times 10^{10}\).

El compositor visual de Harmoniq permite construir circuitos arrastrando puertas sobre un diagrama de líneas de cúbits. Al aplicar una puerta Hadamard \(H\) sobre un cúbit inicializado en \(|0\rangle\), el simulador actualiza el vector de estado a \(|+\rangle=\tfrac{1}{\sqrt{2}}(|0\rangle+|1\rangle)\) en tiempo real, y el panel de probabilidades muestra inmediatamente la distribución \(p(0)=p(1)=0.5\). Cadenas de puertas \(R_z(\varphi)\) permiten explorar cómo la fase relativa modifica las distribuciones al medir en bases complementarias, haciendo visible la interferencia cuántica sin necesidad de ecuaciones explícitas.

Para circuitos más complejos, el SDK de Harmoniq (disponible en Python, Swift, Kotlin y React) expone una API programática que permite construir, parametrizar y ejecutar circuitos con shots configurables. La ejecución con múltiples shots reproduce el muestreo estocástico de la regla de Born, generando histogramas de frecuencia que convergen a las probabilidades teóricas conforme crece el número de repeticiones. Funcionalidades avanzadas incluyen la obtención del vector de estado completo (statevector), el cálculo de valores esperados de observables arbitrarios y la simulación de canales de ruido —amortiguamiento de amplitud, desfase y despolarización— con parámetros de \(T_1\) y \(T_2\) configurables, acercando la simulación a las condiciones de hardware real.

La integración de Harmoniq con los módulos de visualización permite representar el estado del cúbit en la esfera de Bloch en cada paso del circuito, ofreciendo una conexión directa entre la notación abstracta del espacio de Hilbert y la geometría del estado cuántico. Esta capacidad pedagógica es particularmente valiosa para instituciones educativas y equipos de investigación que buscan desarrollar intuición sobre superposición, interferencia y entrelazamiento antes de acceder a hardware cuántico real.

Los modelos de ruido disponibles en Harmoniq permiten además estudiar la degradación de la superposición bajo condiciones realistas. El usuario puede configurar tasas de error de puerta, tiempos \(T_1\) y \(T_2\) por cúbit y errores de lectura (readout errors), observando cómo las distribuciones de medición se alejan de las predicciones ideales conforme aumenta la profundidad del circuito. Esta funcionalidad es esencial para validar la viabilidad de algoritmos cuánticos antes de ejecutarlos en hardware real, donde cada ejecución consume recursos limitados y costosos.