Principio de incertidumbre: incompatibilidad de observables y límites fundamentales

En el formalismo de la mecánica cuántica, toda magnitud física medible —posición, momento, energía, espín— se representa mediante un operador autoadjunto (también denominado hermítico) que actúa sobre el espacio de Hilbert \(\mathcal{H}\) del sistema. Un operador \(\hat{A}\) es autoadjunto si satisface \(\hat{A} = \hat{A}^\dagger\), donde \(\hat{A}^\dagger\) denota el adjunto definido por la relación \(\langle \phi | \hat{A}\psi \rangle = \langle \hat{A}^\dagger \phi | \psi \rangle\) para todos los vectores \(|\phi\rangle, |\psi\rangle \in \mathcal{H}\). Esta condición garantiza que todos los autovalores del operador son reales, requisito indispensable para que los resultados de medición correspondan a cantidades físicas observables.

El teorema espectral proporciona la piedra angular matemática: todo operador autoadjunto \(\hat{A}\) sobre un espacio de Hilbert separable admite una descomposición espectral de la forma

donde \(\sigma(\hat{A})\) es el espectro del operador y \(\{E(a)\}\) constituye una familia de proyectores espectrales (medida espectral con valores de proyección). En el caso de espectro discreto no degenerado, la descomposición se simplifica a \(\hat{A} = \sum_n a_n |a_n\rangle\langle a_n|\), donde \(\{|a_n\rangle\}\) es una base ortonormal de autoestados con autovalores \(a_n\). El postulado de medición de Born establece que, al medir \(\hat{A}\) sobre un estado \(|\psi\rangle\), la probabilidad de obtener el resultado \(a_n\) es \(p(a_n) = |\langle a_n | \psi \rangle|^2\), y el estado colapsa al autoestado correspondiente.

El valor esperado de un observable \(\hat{A}\) en el estado \(|\psi\rangle\) se define como

y la varianza queda determinada por \((\Delta A)^2 = \langle \hat{A}^2 \rangle - \langle \hat{A} \rangle^2\). La desviación estándar \(\Delta A = \sqrt{(\Delta A)^2}\) cuantifica la dispersión estadística intrínseca del observable en el estado dado, con independencia de cualquier limitación instrumental. Es fundamental subrayar que \(\Delta A = 0\) si y solo si \(|\psi\rangle\) es un autoestado de \(\hat{A}\); en cualquier otro caso, la medición exhibe una dispersión genuinamente cuántica. Nótese además que para un estado mixto descrito por un operador densidad \(\hat{\rho}\), el valor esperado se generaliza a \(\langle \hat{A} \rangle = \mathrm{Tr}(\hat{\rho}\,\hat{A})\) y la varianza adopta la forma \((\Delta A)^2 = \mathrm{Tr}(\hat{\rho}\,\hat{A}^2) - [\mathrm{Tr}(\hat{\rho}\,\hat{A})]^2\).

Para operadores sobre espacios de dimensión infinita —como el operador posición \(\hat{x}\) o el momento \(\hat{p}\) en \(L^2(\mathbb{R})\)— la autoadjunción debe distinguirse cuidadosamente de la mera simetría. Un operador simétrico densamente definido satisface \(\langle \phi | \hat{A}\psi \rangle = \langle \hat{A}\phi | \psi \rangle\) para todo \(|\phi\rangle, |\psi\rangle\) en su dominio, pero solo es autoadjunto cuando además \(\mathrm{Dom}(\hat{A}) = \mathrm{Dom}(\hat{A}^\dagger)\). La distinción es crucial: los operadores meramente simétricos pueden carecer de descomposición espectral completa. El teorema de Stone-von Neumann garantiza que, bajo condiciones de irreducibilidad, las relaciones canónicas de conmutación determinan esencialmente una única representación unitariamente equivalente.

El conmutador de dos operadores \(\hat{A}\) y \(\hat{B}\) se define como el operador

Este objeto algebraico codifica de manera precisa la noción de incompatibilidad entre observables. Cuando \([\hat{A}, \hat{B}] = 0\), los operadores conmutan y se dice que los observables correspondientes son compatibles: existe una base ortonormal de autoestados simultáneos de ambos operadores, de modo que pueden medirse conjuntamente sin perturbación mutua. En caso contrario, \([\hat{A}, \hat{B}] \neq 0\) señala una incompatibilidad fundamental que impide la existencia de autoestados comunes y, como veremos, impone cotas inferiores al producto de sus dispersiones.

El conmutador posee propiedades algebraicas notables. Es bilineal, \([\alpha\hat{A}+\beta\hat{B}, \hat{C}] = \alpha[\hat{A},\hat{C}] + \beta[\hat{B},\hat{C}]\), antisimétrico, \([\hat{A},\hat{B}] = -[\hat{B},\hat{A}]\), y satisface la identidad de Jacobi:

Estas tres propiedades confieren al espacio de operadores la estructura de un álgebra de Lie. En particular, las relaciones canónicas de conmutación del par posición-momento,

definen el álgebra de Heisenberg-Weyl, cuya representación irreducible es esencialmente única (módulo equivalencia unitaria) por el teorema de Stone-von Neumann. Para sistemas de espín, las componentes del momento angular satisfacen \([\hat{J}_i, \hat{J}_j] = i\hbar\,\epsilon_{ijk}\hat{J}_k\), constituyendo el álgebra de Lie \(\mathfrak{su}(2)\). El anticonmutador \(\{\hat{A}, \hat{B}\} = \hat{A}\hat{B} + \hat{B}\hat{A}\) complementa al conmutador y permite descomponer el producto de operadores como \(\hat{A}\hat{B} = \tfrac{1}{2}[\hat{A},\hat{B}] + \tfrac{1}{2}\{\hat{A},\hat{B}\}\), descomposición que resultará fundamental en la derivación de la relación de Robertson-Schrödinger.

La identidad de Leibniz para conmutadores, \([\hat{A}, \hat{B}\hat{C}] = [\hat{A},\hat{B}]\hat{C} + \hat{B}[\hat{A},\hat{C}]\), permite calcular conmutadores de productos de operadores de forma sistemática. Por ejemplo, para el hamiltoniano del oscilador armónico \(\hat{H} = \hat{p}^2/(2m) + m\omega^2\hat{x}^2/2\), se obtienen inmediatamente \([\hat{H}, \hat{x}] = -i\hbar\hat{p}/m\) y \([\hat{H}, \hat{p}] = i\hbar m\omega^2\hat{x}\), que reproducen las ecuaciones de movimiento de Heisenberg \(d\hat{x}/dt = \hat{p}/m\) y \(d\hat{p}/dt = -m\omega^2\hat{x}\), en perfecta correspondencia con la dinámica clásica.

En dimensión finita, especialmente relevante para qubits y registros cuánticos, los observables se representan mediante matrices hermitianas \(A, B \in \mathbb{C}^{d \times d}\). Las matrices de Pauli \(\sigma_x, \sigma_y, \sigma_z\) satisfacen \([\sigma_i, \sigma_j] = 2i\epsilon_{ijk}\sigma_k\), y junto con la identidad forman una base para el espacio de matrices hermitianas \(2 \times 2\). Para un sistema de \(n\) qubits, el espacio de operadores tiene dimensión \(4^n\), y la estructura de conmutación entre productos tensoriales de matrices de Pauli — el grupo de Pauli — gobierna las relaciones de incertidumbre de todos los observables del registro.

La desigualdad de Robertson (1929) establece la relación de incertidumbre más general para un par arbitrario de observables. Para operadores autoadjuntos \(\hat{A}\) y \(\hat{B}\) y cualquier estado normalizado \(|\psi\rangle\), se cumple

donde \(\Delta A = \sqrt{\langle \hat{A}^2 \rangle - \langle \hat{A} \rangle^2}\) y análogamente para \(\Delta B\). La demostración emplea la desigualdad de Cauchy-Schwarz aplicada a los vectores \(|\alpha\rangle = (\hat{A} - \langle\hat{A}\rangle)|\psi\rangle\) y \(|\beta\rangle = (\hat{B} - \langle\hat{B}\rangle)|\psi\rangle\). Se tiene \((\Delta A)^2(\Delta B)^2 = \langle\alpha|\alpha\rangle\langle\beta|\beta\rangle \geq |\langle\alpha|\beta\rangle|^2\). Puesto que \(\langle\alpha|\beta\rangle\) es en general complejo, se descompone en parte real e imaginaria:

donde hemos definido \(\hat{A}' = \hat{A} - \langle\hat{A}\rangle\) y \(\hat{B}' = \hat{B} - \langle\hat{B}\rangle\). Al tomar solo el módulo de la parte imaginaria se obtiene la desigualdad de Robertson. Schrödinger (1930) observó que, al retener también la parte real del producto interno, se obtiene una cota más fuerte:

Esta es la relación de Robertson-Schrödinger, que domina estrictamente la de Robertson siempre que el valor esperado del anticonmutador centrado no se anule. La cota de Schrödinger incluye la covarianza cuántica \(C_{AB} = \tfrac{1}{2}\langle\{\hat{A}',\hat{B}'\}\rangle = \tfrac{1}{2}\langle \hat{A}\hat{B} + \hat{B}\hat{A}\rangle - \langle\hat{A}\rangle\langle\hat{B}\rangle\), que captura las correlaciones clásicas entre ambos observables. En notación compacta:

Para comprender la saturación, consideremos cuándo se alcanza la igualdad. La desigualdad de Cauchy-Schwarz se satura si y solo si \(|\beta\rangle = \lambda|\alpha\rangle\) para algún \(\lambda \in \mathbb{C}\), es decir, cuando \((\hat{B} - \langle\hat{B}\rangle)|\psi\rangle = \lambda(\hat{A} - \langle\hat{A}\rangle)|\psi\rangle\). Los estados que saturan la relación de Robertson (sin el término del anticonmutador) satisfacen \(\lambda \in i\mathbb{R}\), mientras que la saturación completa de Robertson-Schrödinger admite \(\lambda\) complejo general. Los estados gaussianos comprimidos (squeezed states) proporcionan ejemplos paradigmáticos de saturación con correlaciones no nulas.

Resulta instructivo reformular la desigualdad en términos de la matriz de covarianza cuántica. Para un vector de operadores \(\hat{\mathbf{R}} = (\hat{A}_1, \hat{A}_2, \ldots, \hat{A}_n)^T\), se define la matriz de covarianza \(\Gamma_{ij} = \tfrac{1}{2}\langle\{\hat{A}_i', \hat{A}_j'\}\rangle\) y la matriz simpléctica \(\Omega_{ij} = -\tfrac{i}{2}\langle[\hat{A}_i, \hat{A}_j]\rangle\). La desigualdad de Robertson-Schrödinger se generaliza entonces a la condición matricial \(\Gamma + i\Omega \geq 0\), que constituye la base del formalismo de estados gaussianos en óptica cuántica y variables continuas.

La instancia más célebre del principio de incertidumbre surge al aplicar la relación de Robertson al par canónico posición-momento. Dado que \([\hat{x}, \hat{p}] = i\hbar\), el conmutador es proporcional a la identidad y su valor esperado es independiente del estado:

Esta desigualdad fue formulada originalmente por Heisenberg en 1927, aunque su argumento inicial invocaba un experimento mental —el microscopio de rayos gamma— que mezclaba la incertidumbre cuántica intrínseca con la perturbación causada por la medición. La derivación rigurosa a partir del formalismo de operadores, debida a Kennard (1927) y Robertson (1929), muestra que la cota es una propiedad puramente cinemática del espacio de Hilbert, independiente del aparato de medición empleado.

Los estados de mínima incertidumbre que saturan la desigualdad de Heisenberg son los paquetes de onda gaussianos. En la representación de posición, un estado de mínima incertidumbre centrado en \(x_0\) con momento medio \(p_0\) adopta la forma

con \(\Delta x = \sigma_x\) y \(\Delta p = \hbar/(2\sigma_x)\), de modo que \(\Delta x \cdot \Delta p = \hbar/2\) exactamente. La transformada de Fourier de este paquete es también una gaussiana con anchura \(\sigma_p = \hbar/(2\sigma_x)\), ilustrando la dualidad posición-momento como una manifestación del principio de incertidumbre.

Los estados coherentes \(|\alpha\rangle\) del oscilador armónico cuántico —autoestados del operador de aniquilación, \(\hat{a}|\alpha\rangle = \alpha|\alpha\rangle\)— son precisamente estados de mínima incertidumbre con \(\Delta x = \Delta p \cdot (m\omega)^{-1}\) en las variables adimensionales adecuadas, y con \(\Delta x \cdot \Delta p = \hbar/2\). Introducidos por Glauber en 1963, estos estados son los que más se asemejan a las ondas clásicas y desempeñan un papel central en la óptica cuántica. Evolucionan en el tiempo manteniendo su forma gaussiana, con el centro del paquete siguiendo la trayectoria clásica — una propiedad que los convierte en un puente natural entre la mecánica clásica y la cuántica.

Es posible construir estados comprimidos (squeezed states) que redistribuyen la incertidumbre entre las cuadraturas del campo, reduciendo \(\Delta x\) por debajo del nivel del estado coherente a costa de aumentar \(\Delta p\) (o viceversa), manteniendo siempre \(\Delta x \cdot \Delta p \geq \hbar/2\). Formalmente, un estado comprimido se obtiene aplicando el operador de compresión \(\hat{S}(\xi) = \exp[\tfrac{1}{2}(\xi^*\hat{a}^2 - \xi\hat{a}^{\dagger 2})]\) al estado de vacío. Estos estados encuentran aplicaciones cruciales en interferometría gravitacional (LIGO) y en la mejora de la sensibilidad metrológica más allá del límite cuántico estándar.

La relación de incertidumbre energía-tiempo, frecuentemente escrita como

requiere un tratamiento conceptual distinto al del par posición-momento, dado que en mecánica cuántica no relativista el tiempo no es un observable representado por un operador autoadjunto, sino un parámetro externo de la ecuación de Schrödinger. Por lo tanto, la relación no se deriva directamente de Robertson, y su interpretación ha sido objeto de extenso debate desde los primeros años de la teoría cuántica.

La formulación más rigurosa y ampliamente aceptada es la de Mandelstam y Tamm (1945). Para un observable arbitrario \(\hat{A}\) que no depende explícitamente del tiempo, la ecuación de movimiento de Heisenberg da \(d\langle\hat{A}\rangle/dt = (i/\hbar)\langle[\hat{H}, \hat{A}]\rangle\). Aplicando la relación de Robertson al par \((\hat{H}, \hat{A})\):

Definiendo el tiempo característico de evolución del observable como \(\tau_A = \Delta A / |d\langle\hat{A}\rangle/dt|\), se obtiene \(\Delta E \cdot \tau_A \geq \hbar/2\). Aquí \(\tau_A\) mide el tiempo que el sistema necesita para que el valor esperado de \(\hat{A}\) cambie en una cantidad comparable a su propia dispersión. Puesto que esta relación es válida para todo observable \(\hat{A}\), se define \(\Delta t = \min_{\hat{A}} \tau_A\), obteniendo la relación en su forma habitual.

La interpretación de Mandelstam-Tamm conduce directamente al concepto de límite cuántico de velocidad (quantum speed limit): existe una cota inferior al tiempo necesario para que un sistema cuántico evolucione desde un estado inicial \(|\psi_0\rangle\) hasta un estado ortogonal \(|\psi_\perp\rangle\). Si \(\Delta E\) es la dispersión de energía del estado inicial, el tiempo mínimo de tránsito satisface

Esta cota, demostrada originalmente por Mandelstam y Tamm y refinada por Margolus y Levitin, tiene consecuencias profundas: establece un límite fundamental a la velocidad de procesamiento de información cuántica, un bound termodinámico sobre la tasa de operaciones lógicas de cualquier computadora (cuántica o clásica) con energía finita, y conecta la mecánica cuántica con la teoría de la información a través de la noción de capacidad de canal.

Otra manifestación clásica de la relación energía-tiempo es la relación anchura de línea-vida media. Un estado excitado con vida media \(\tau\) presenta un espectro de emisión con anchura natural \(\Gamma \sim \hbar/\tau\), lo que implica que estados de vida larga emiten líneas espectrales estrechas y viceversa. En espectroscopía atómica, esta relación conecta directamente los tiempos de decaimiento radiativos con las anchuras de las líneas de emisión observadas, y en física de partículas relaciona los anchos de desintegración de resonancias con sus vidas medias a través de \(\Gamma \cdot \tau = \hbar\).

Dos bases ortonormales \(\mathcal{B}_1 = \{|e_i\rangle\}\) y \(\mathcal{B}_2 = \{|f_j\rangle\}\) de un espacio de Hilbert de dimensión \(d\) se denominan mutuamente insesgadas (mutually unbiased bases, MUB) si para todo par de vectores base se cumple

Esta condición significa que si un sistema se prepara en un autoestado de la primera base y se mide en la segunda, todos los resultados son equiprobables — se tiene ignorancia máxima. Las MUB constituyen la manifestación discreta más limpia del principio de incertidumbre: la certeza total en un observable implica la incertidumbre total en el observable conjugado.

El ejemplo paradigmático es el sistema de un qubit (\(d = 2\)). Las tres bases de Pauli conforman un conjunto maximal de MUB:

donde \(|{\pm}\rangle = (|0\rangle \pm |1\rangle)/\sqrt{2}\) y \(|{\pm i}\rangle = (|0\rangle \pm i|1\rangle)/\sqrt{2}\). En dimensión \(d\), el número máximo de MUB es a lo sumo \(d + 1\), y cuando \(d\) es potencia de un primo, se alcanza este máximo mediante construcciones explícitas basadas en cuerpos finitos \(\mathbb{F}_{d}\). Para dimensiones no primas, como \(d = 6\), la determinación del número máximo de MUB permanece como uno de los problemas abiertos más notables de la teoría de la información cuántica.

Las MUB desempeñan un papel central en la distribución cuántica de claves (QKD). En el protocolo BB84, Alice codifica bits aleatorios eligiendo entre dos bases mutuamente insesgadas (típicamente \(\mathcal{B}_Z\) y \(\mathcal{B}_X\)) y envía qubits a Bob, quien mide aleatoriamente en una de las dos bases. Cuando ambos eligen la misma base, obtienen resultados perfectamente correlacionados; cuando eligen bases distintas, los resultados de Bob son completamente aleatorios. Un espía (Eve) que intercepte y remida los qubits inevitablemente introduce errores detectables, precisamente porque la medición en una base insesgada destruye irremediablemente la información codificada en la otra. La seguridad del protocolo descansa directamente sobre el principio de incertidumbre.

Más allá de la criptografía, las MUB aparecen en la tomografía de estados cuánticos: un conjunto completo de \(d+1\) MUB permite reconstruir un operador densidad \(\hat{\rho}\) de dimensión \(d\) con el mínimo número de configuraciones de medición distintas. En la reconstrucción tomográfica, la probabilidad \(p_{b,k} = \langle f_k^{(b)} | \hat{\rho} | f_k^{(b)} \rangle\) medida para cada base \(b\) y cada resultado \(k\) proporciona información complementaria y maximamente informativa sobre el estado, haciendo que las MUB sean las bases óptimas para la tomografía.

Las relaciones de incertidumbre basadas en varianzas, como la de Robertson, presentan limitaciones: la cota depende del estado a través de \(\langle[\hat{A},\hat{B}]\rangle\), que puede anularse incluso para observables genuinamente incompatibles (por ejemplo, para el estado \(|0\rangle\) con \(\hat{A} = \sigma_x\) y \(\hat{B} = \sigma_y\), el conmutador es proporcional a \(\sigma_z\) cuyo valor esperado es \(\pm 1\), pero consideremos \(\hat{A} = \sigma_x\) y \(\hat{B} = \sigma_z\): para el estado \(|+\rangle\) se tiene \(\langle[\sigma_x, \sigma_z]\rangle = 0\)). Las relaciones de incertidumbre entrópicas superan esta deficiencia al cuantificar la incertidumbre mediante entropías en lugar de varianzas.

Para un observable discreto \(\hat{A}\) con autovalores \(\{a_k\}\) y autoestados \(\{|a_k\rangle\}\), la entropía de Shannon de la distribución de resultados de medición se define como

Deutsch (1983) fue el primero en conjeturar una relación de incertidumbre entrópica, que fue demostrada de forma óptima por Maassen y Uffink (1988):

donde \(c(\hat{A}, \hat{B}) = \max_{j,k} |\langle a_j | b_k \rangle|\) es la máxima superposición entre los autoestados de ambos observables y el logaritmo se toma en la base apropiada (base 2 para bits, base natural para nats). Para bases mutuamente insesgadas, \(c = 1/\sqrt{d}\), y la cota se convierte en \(H(\hat{A}) + H(\hat{B}) \geq \log d\), que es independiente del estado y estrictamente positiva siempre que \(d \geq 2\).

La superioridad de las relaciones entrópicas radica en varias propiedades. Primero, la cota de Maassen-Uffink depende solo de la estructura geométrica de las bases de medición, no del estado particular del sistema, a diferencia de la relación de Robertson. Segundo, las entropías capturan la forma completa de la distribución de probabilidad, no solo su segundo momento. Tercero, admiten generalizaciones naturales al caso de estados mixtos y mediciones con memoria cuántica: si el sistema \(A\) está entrelazado con una memoria cuántica \(B\), la relación de Maassen-Uffink se refuerza según

donde \(S(\hat{A}|B)_\rho\) es la entropía condicional de von Neumann tras la medición de \(\hat{A}\) en el subsistema, y \(S(A|B)_\rho\) es la entropía condicional del estado previo a la medición. Este resultado, demostrado por Berta, Christandl, Colbeck, Renes y Renner (2010), tiene consecuencias fundamentales: si \(A\) y \(B\) están maximamente entrelazados (\(S(A|B) = -\log d\)), la cota se anula, reflejando que el entrelazamiento puede compensar la incertidumbre.

La entropía de von Neumann \(S(\hat{\rho}) = -\mathrm{Tr}(\hat{\rho}\log\hat{\rho})\) generaliza la entropía de Shannon al caso de estados mixtos y proporciona la medida natural de incertidumbre en el formalismo de operadores densidad. Las relaciones de incertidumbre entrópicas formuladas en términos de entropías de Rényi, \(H_\alpha(\hat{A}) = (1-\alpha)^{-1}\log\sum_k p_k^\alpha\), ofrecen una familia parametrizada que interpola entre la min-entropía (\(\alpha \to \infty\), relevante para criptografía) y la entropía de Hartley (\(\alpha \to 0\), que cuenta el soporte efectivo de la distribución).

La metrología cuántica explota las propiedades del principio de incertidumbre para determinar los límites fundamentales de precisión en la estimación de parámetros físicos. Consideremos la tarea canónica de estimar una fase \(\theta\) acumulada durante la evolución unitaria \(\hat{U}(\theta) = e^{-i\theta\hat{G}}\), donde \(\hat{G}\) es el generador hermitiano de la transformación. La incertidumbre en la estimación de \(\theta\) está acotada inferiormente por la cota de Cramér-Rao cuántica:

donde \(\nu\) es el número de repeticiones del experimento y \(F_Q[\hat{\rho}, \hat{G}]\) es la información cuántica de Fisher (QFI). Para estados puros, \(F_Q = 4(\Delta G)^2\), lo que relaciona directamente la precisión metrológica con la dispersión del generador y, por extensión, con el principio de incertidumbre.

Si se emplean \(N\) partículas independientes (no entrelazadas) como sondas, cada una con dispersión \(\Delta G = \Delta G_0\), la QFI total escala como \(F_Q = N \cdot 4(\Delta G_0)^2\), conduciendo al límite cuántico estándar (SQL, standard quantum limit):

que escala como \(1/\sqrt{N}\) — el comportamiento de shot noise clásico. Sin embargo, utilizando estados entrelazados como el estado de Greenberger-Horne-Zeilinger (GHZ), \(|GHZ\rangle = (|0\rangle^{\otimes N} + |1\rangle^{\otimes N})/\sqrt{2}\), la QFI alcanza \(F_Q = 4N^2(\Delta G_0)^2\), lo que conduce al límite de Heisenberg:

que representa una mejora cuadrática sobre el SQL. El límite de Heisenberg es la cota fundamental permitida por la mecánica cuántica y no puede superarse sin violar el principio de incertidumbre. La transición del SQL al límite de Heisenberg mediante entrelazamiento constituye una de las demostraciones más directas de la ventaja cuántica en una tarea práctica.

Los estados comprimidos (squeezed states) ofrecen una vía práctica intermedia. En un interferómetro Mach-Zehnder alimentado con luz comprimida, la sensibilidad de fase mejora con el parámetro de compresión \(r\) según \(\delta\theta \sim e^{-r}/\sqrt{N}\), permitiendo superar el SQL sin requerir estados altamente entrelazados. LIGO utiliza esta técnica desde 2019, inyectando vacío comprimido (squeezed vacuum) en el puerto oscuro del interferómetro para reducir el ruido cuántico de disparo y acercarse al límite de Heisenberg.

Es importante señalar que, en presencia de decoherencia, el escalado de Heisenberg \(1/N\) generalmente se degrada. Para canales de pérdida con transmisividad \(\eta\), la cota de Cramér-Rao cuántica transita del escalado de Heisenberg al SQL según \(\delta\theta \geq 1/(\sqrt{N} \cdot 2\Delta G_0 \cdot \sqrt{\eta})\) en el régimen de alta pérdida. Este resultado tiene implicaciones directas para la metrología con procesadores NISQ, donde las tasas de error de compuerta limitan la profundidad de circuito utilizable para la preparación de estados sonda óptimos.

En la era de la computación cuántica de escala intermedia con ruido (NISQ, Noisy Intermediate-Scale Quantum), las relaciones de incertidumbre adquieren relevancia operativa directa. Los algoritmos variacionales — como VQE (Variational Quantum Eigensolver) y QAOA (Quantum Approximate Optimization Algorithm) — requieren la estimación de valores esperados \(\langle\hat{O}\rangle = \langle\psi(\boldsymbol{\theta})|\hat{O}|\psi(\boldsymbol{\theta})\rangle\) mediante repetición de mediciones sobre el circuito cuántico parametrizado. La incertidumbre estadística en esta estimación obedece al teorema del límite central:

donde \(\Delta O\) es la desviación estándar del observable en el estado preparado y \(N_{\mathrm{shots}}\) es el número de ejecuciones (disparos) del circuito. Esta relación cuantifica el ruido de disparo (shot noise) inherente a la medición cuántica y establece un compromiso fundamental: para reducir el error de estimación a la mitad, es necesario cuadruplicar el número de disparos. Para un hamiltoniano descompuesto como suma de \(M\) términos de Pauli, \(\hat{H} = \sum_{k=1}^M c_k \hat{P}_k\), el error total en la estimación de la energía escala como

lo que motiva técnicas de agrupación de términos conmutantes (grouping) y asignación óptima de disparos para minimizar la varianza total.

Los errores de compuerta en procesadores NISQ introducen una fuente adicional de incertidumbre que se superpone al ruido de disparo. Una compuerta de un qubit imperfecta puede modelarse como una rotación errónea \(\hat{U}_{\mathrm{real}} = e^{-i(\theta+\delta\theta)\hat{n}\cdot\hat{\boldsymbol{\sigma}}/2}\), donde \(\delta\theta\) fluctúa entre ejecuciones. La incertidumbre en el ángulo de rotación se propaga a los valores esperados de los observables posteriores a través de la cadena de compuertas del circuito. Para un circuito de profundidad \(L\) con tasa de error de compuerta \(\epsilon_g\), la fidelidad del estado final escala aproximadamente como \(\mathcal{F} \sim (1 - \epsilon_g)^L\), y la dispersión efectiva del observable medido se amplifica correspondientemente.

Los errores de lectura (readout errors) constituyen otra manifestación del principio de incertidumbre en la práctica. En procesadores superconductores, la discriminación entre los estados \(|0\rangle\) y \(|1\rangle\) se realiza mediante la medición de la señal reflejada por un resonador acoplado dispersivamente al qubit. La señal integrada durante un tiempo \(T_m\) presenta una relación señal-ruido que crece como \(\mathrm{SNR} \propto \sqrt{T_m}\), pero el tiempo de medición está acotado superiormente por el tiempo de relajación \(T_1\) del qubit. La tasa de error de clasificación típica en dispositivos actuales oscila entre el 0.5% y el 5%, y puede reducirse mediante protocolos de readout error mitigation basados en matrices de confusión calibradas.

El diseño de experimentos en procesadores NISQ requiere un equilibrio cuidadoso entre todas estas fuentes de incertidumbre. El presupuesto de error total para la estimación de un valor esperado combina el ruido de disparo (\(\propto 1/\sqrt{N_{\mathrm{shots}}}\)), los errores de compuerta (\(\propto L\epsilon_g\)), los errores de lectura y las contribuciones de decoherencia. La elección del número de disparos \(N_{\mathrm{shots}}\) se determina a partir de la tolerancia al error estadístico:

Para estimaciones de energía con precisión química (\(\epsilon \sim 1.6 \times 10^{-3}\) hartrees) en moléculas moderadas, el número de disparos requerido puede alcanzar \(10^6\)–\(10^9\), lo que constituye un cuello de botella significativo en la ejecución de algoritmos variacionales. Técnicas como la estimación clásica de sombras (classical shadows), el muestreo de Pauli adaptativo y los protocolos de mitigación de errores por extrapolación al ruido cero (zero-noise extrapolation) buscan aliviar estas demandas, redistribuyendo el presupuesto de incertidumbre entre las distintas contribuciones para maximizar la información extraíble por disparo.