Módulo II.4 · Hoja de problemas

Hoja de problemas del módulo II.4

Seis problemas sobre el formalismo, con pistas y sin soluciones. La brújula: en este módulo casi ninguna propiedad se lee, se calcula — y la cuenta que la decide casi nunca es la que parece.

La matriz (1102)\bigl(\begin{smallmatrix}1&1\\0&2\end{smallmatrix}\bigr) tiene por autovalores 1 y 2 —reales, redondos, respetables— y con sus autovectores la suma de la regla de Born no vale 1: recorre de 0,293 a 1,707 según cuál sea el estado. Ahí está, en una línea, el error que más veces se comete en este tema: leer una propiedad en la superficie del objeto en vez de calcular la cuenta que la decide. El espectro no certifica que algo sea un observable; lo certifica un conmutador. Una cota que sale cero no es una cota pequeña. Una suma truncada que no llega a 1 no mide la base, mide el corte. Los seis problemas de abajo pasan por los tres ejes.

El primero: qué compra cada hipótesis, exactamente. La hermiticidad compra autovalores reales; la ortogonalidad de los autovectores la compra ser normal, que es otra condición y se comprueba con [A, A†]. Hay matrices no hermíticas con los autovectores perfectamente perpendiculares, y el problema 1 tiene una en la mano. Antes de escribir «como no es hermítica…», di qué se sigue de eso y qué no.

El segundo: la cota es de la pareja y del estado. Lo dijo el artículo 03 con el panel: de los 65 341 estados de su rejilla saturan 1439, y no son dos puntos sino dos circunferencias. En una desviación típica eso se traduce en un hábito: si una cota te sale cero, mira si el producto también lo es —en cuyo caso no has encontrado un contraejemplo, has encontrado una identidad— y si te sale positiva, mira cuánto le sobra al producto, porque el sobrante tiene dos orígenes distintos y se pueden separar.

El tercero: convergencia no es una sola cosa. Una función puede desarrollarse entera en una base —Parseval exacto, suma 1— y tener a la vez energía media infinita. Las dos series son distintas y hay funciones perfectamente respetables en las que una converge y la otra no. El problema 4 tiene una, y es la más sencilla que se puede dibujar.

Y una brújula de método, que ahorra páginas: en este módulo casi todo se decide con una cuenta de una línea, y casi nunca es la cuenta larga. ¿Son ortogonales estos autovectores? Un conmutador, no un problema de autovalores. ¿Vale cero esta cota? Mira si el estado es real antes de calcular ⟨p⟩. ¿Es completa esta base? Ninguna suma finita lo demuestra, y una sola suma que falle lo refuta. ¿Depende del orden esta pareja de medidas? El conmutador otra vez. Lo caro es lanzarse a integrar lo que se contesta mirando una simetría.

Necesitas: el módulo II.4 completo. En concreto, el 1 continúa el ejercicio 5 del artículo 02, que clasifica su matriz y no mira sus autovectores, y su último apartado vuelve sobre la derivada discretizada del ejercicio 3 del mismo artículo. El 2 continúa el ejercicio 6 del artículo 01. El 3 pide la relación de Robertson del artículo 03 sobre el operador número del II.3, artículo 03, del que hacen falta el álgebra de a y a† y los conmutadores [N, a] = −a y [N, a†] = a†; su apartado (c) arranca donde acaba el ejercicio 4 del artículo 04, que mide N sobre ese mismo estado y de cuyo reparto 0,5 / 0,5 sale σN en un renglón. El 4 continúa el ejercicio 4 del artículo 01 con una función que el panel no puede generar. El 5 continúa el ejercicio 6 del artículo 03 y repite, sobre otro sistema y otros pesos, la cadena ⟨H⟩ → σE → τ → T/τ del ejercicio 3 de ese mismo artículo; compara además con la mezcla del pozo de 1 nm del II.2. El 6 es el artículo 04 entero: los postulados III y IV aplicados en cadena.
Los datos van impresos en cada enunciado, con las cifras que hacen falta para comprobarse: no hay que ir a buscarlos a los artículos. Las constantes son las de CODATA 2022, y en todo el módulo ħ = 0,658212 eV·fs.
Nada de esta hoja necesita el II.5, y nada se contesta arrastrando un mando de los dos paneles: el del artículo 01 aparece en el problema 4 como objeto de estudio, no como calculadora.
Problema 1 ●○○ · Una matriz que finge, y unos autovectores que no

El ejercicio 5 del artículo 02 clasificó A3=(0110)A_3 = \bigl(\begin{smallmatrix}0&-1\\1&0\end{smallmatrix}\bigr) en tres casillas —normal, espectro real, hermítica— y ahí la dejó. Aquí se toma su traspuesta y se le pregunta lo que aquella clasificación no pregunta: dónde están sus autovectores.

Datos. A=(0110)A = \bigl(\begin{smallmatrix}0&1\\-1&0\end{smallmatrix}\bigr), las matrices de Pauli σx=(0110)\sigma_x = \bigl(\begin{smallmatrix}0&1\\1&0\end{smallmatrix}\bigr), σy=(0ii0)\sigma_y = \bigl(\begin{smallmatrix}0&-i\\i&0\end{smallmatrix}\bigr), σz=(1001)\sigma_z = \bigl(\begin{smallmatrix}1&0\\0&-1\end{smallmatrix}\bigr). Para el apartado (e): un electrón en un anillo de perímetro L = 1,000 nm, con m = me; el II.2 publica que su primer modo está en 1,504 eV, cuatro veces la E1 = 0,376 eV del pozo del mismo tamaño. La derivada se discretiza en 64 nodos equiespaciados con diferencias centradas, (Dψ)j=(ψj+1ψj1)/2Δx(D\psi)_j = (\psi_{j+1}-\psi_{j-1})/2\,\Delta x, y los índices se cierran en círculo.

(a) Escribe A† y di si A es hermítica. Mide cuánto no lo es con ‖A − A†‖, tomando como norma el mayor elemento en valor absoluto. ¿Puede A ser un observable? Da el nombre de lo que sí es. (b) Sus autovalores, por el polinomio característico, y compruébalos con la traza y el determinante, que tienen que devolverte su suma y su producto. Ahora los autovectores normalizados y su solape |⟨v+|v⟩|. Sale ortogonal — y la parte que hay que contestar es ésta: ¿qué comprobación de una línea lo predecía antes de calcular ningún autovector? Y entonces, ¿qué le compra exactamente la hermiticidad a un observable, si no es la ortogonalidad? (c) Ahora iA: ¿es hermítica? ¿Cuál es su espectro? ¿Y sus autovectores? Compáralos uno a uno con los de A y di qué ha movido la i y qué ha dejado quieto. Identifica iA con una matriz de Pauli, y localiza sus dos autoestados en el panel del artículo 03: ¿tienen algo de particular allí? (d) Calcula A†A y A². Además de normal, A es otra cosa: dile el nombre. Con A² = −1̂, suma la serie de eθAe^{\theta A} agrupando términos pares e impares, di qué transformación del plano es y cuánto vale e2πAe^{2\pi A}. Escribe en una frase la regla general de la que esto es el caso 2 × 2, y señala en qué renglón del postulado V se usa. (e) La misma estructura, en grande. La matriz D del anillo es real y antisimétrica: sin diagonalizar nada, ¿qué forma tiene su espectro y qué forma tiene el de −iħD? Prueba la onda plana ψj=eikxj\psi_j = e^{ikx_j} con k=2πm/Lk = 2\pi m/L y obtén el autovalor exacto de la malla. Evalúalo para m = 1 en kg·m/s y compáralo con el ħk del continuo: ¿por cuánto se queda corta la malla, y en qué signo? ¿Qué le pasa a ese autovalor cuando k crece hasta el borde —qué tiene el operador discreto que el continuo no tiene—?

Pista

El camino de (b) no pasa por los autovectores: pasa por [A, A†]. El teorema espectral general del artículo 02 dice «normal ⟺ diagonalizable en base ortonormal», y A conmuta con su adjunta de la manera más barata posible. En (e), una matriz circulante tiene por autovectores las ondas planas que caben en el anillo, y no hay que resolver ningún sistema: basta sustituir.

Ojo con la implicación que este problema existe para romper: «no hermítica ⇒ autovectores no ortogonales» es falsa, y el artículo 02 tiene el contraejemplo dicho al revés —la triangular de la entradilla es no normal, y ésa sí pierde la ortogonalidad—. En (c), multiplicar por i no es un cambio de base: mira dónde acaban los autovectores antes de escribir nada. Y en (e), el factor que separa la malla del continuo no es un error de redondeo: es una función de k, y tiene máximo.

Comprobación: A² tiene que salirte proporcional a la identidad, y con eso eθAe^{\theta A} se cierra en senos y cosenos sin sumar ninguna serie a mano. Comprueba entonces tu eθAe^{\theta A} en θ = 0 y en θ = π/2 aplicándolo a un vector cualquiera y mirando dónde acaba: si la transformación no es la que has dicho que es, has perdido un signo en los términos pares. En (e), tu autovalor tiene que tender a ħk cuando Δx → 0 con k fija: si no, has confundido Δx con L.

Problema 2 ●○○ · Tres direcciones, dos repartos y un trío que no es base

El ejercicio 6 del artículo 01 recortó la resolución de la identidad quitándole un término y midió lo que se perdía. Aquí no se le quita nada a la base: se cambia de base, y luego se usa una que no lo es.

Datos. En ℂ³, el estado ψ=(1,1,1)/3|\psi\rangle = (1, 1, 1)/\sqrt3. Tres tercias:

{ek}: (1,1,0)2, (1,1,0)2, (0,0,1);{uk}: (1,0,0), (0,1,0), (0,0,1);{gk}: (1,1,0)2, (1,0,0), (0,0,1).\{e_k\}:\ \tfrac{(1,1,0)}{\sqrt2},\ \tfrac{(1,-1,0)}{\sqrt2},\ (0,0,1); \qquad \{u_k\}:\ (1,0,0),\ (0,1,0),\ (0,0,1); \qquad \{g_k\}:\ \tfrac{(1,1,0)}{\sqrt2},\ (1,0,0),\ (0,0,1).

(a) Comprueba que los eke_k son una base ortonormal y que cumplen kekek=1^\sum_k |e_k\rangle\langle e_k| = \hat{\mathbb 1} —las dos cosas, porque son dos comprobaciones distintas y sólo la segunda es la resolución de la identidad—. (b) Calcula los tres ck=ekψc_k = \langle e_k|\psi\rangle y los tres |ck|², y después los del mismo |ψ⟩ en la base canónica {uk}\{u_k\}. Los dos repartos son distintos y las dos sumas valen lo mismo. Escribe en una frase qué mide entonces |ck|² y de qué depende, y qué parte de un valor esperado se lleva la base y qué parte se lleva el estado. (c) ¿Qué ángulo forma |ψ⟩ con e1e_1, y cuál con e3e_3? Comprueba que suman 90° y explica por qué tenían que sumar eso y no otra cosa. Uno de los dos tiene nombre propio en resonancia magnética y se define por una condición de una línea sobre su coseno: encuéntrala. (d) Ahora {gk}\{g_k\}. Comprueba que los tres son unitarios y que generan ℂ³ —el determinante no se anula—, y calcula kgkψ2\sum_k |\langle g_k|\psi\rangle|^{2}. ¿Qué acaba de detectar Parseval, y cuál de las dos hipótesis de (a) es la que falla? Intenta reconstruir |ψ⟩ con kgkψgk\sum_k \langle g_k|\psi\rangle\,|g_k\rangle y mide el error. Existe una tercia asociada con la que la reconstrucción sí funciona: di cómo se calcula y por qué en una base ortonormal coincide con la de partida.

Pista

Dos de los tres coeficientes de (b) salen sin calcular nada: uno por simetría entre las dos primeras componentes de |ψ⟩, y otro porque e3e_3 es un eje canónico. En (d), la tercia {gk}\{g_k\} se distingue de {ek}\{e_k\} en un solo vector, y ese vector no es perpendicular a otro de la lista.

Ojo: «normalizados» y «linealmente independientes» son las dos propiedades que {gk}\{g_k\} sí tiene, y no bastan. La trampa consiste en leer gkψ\langle g_k|\psi\rangle como «el coeficiente de |ψ⟩ en la dirección k»: proyección y coeficiente son la misma cosa sólo si la base es ortonormal, y la suma de (d) es lo que mide la diferencia.

Comprobación: la suma de (d) tiene que salirte mayor que 1 y con denominador 3. Si te sale menor que 1, has proyectado sobre vectores sin normalizar. Y en (c), el ángulo con e3e_3 tiene que ser el mismo que |ψ⟩ forma con cada uno de los tres ejes canónicos —es la simetría de (1,1,1)—: si no lo es, revisa el coseno.

Problema 3 ●●○ · Robertson donde el conmutador no es un número

Kennard es el caso en que [x̂, p̂] = iħ es un múltiplo de la identidad y la cota no depende del estado. Aquí se toma la otra pareja natural del oscilador —el número de cuantos y la posición—, cuyo conmutador es un operador, y se mira qué le pasa a la cota estado por estado.

Datos. El electrón del II.3 con x0 = 1,000 nm: ħω = 76,20 meV, k = 12,21 mN/m, m = me. Del mismo módulo, x^=x02(a^+a^)\hat x = \tfrac{x_0}{\sqrt2}(\hat a + \hat a^{\dagger}), p^=ix02(a^a^)\hat p = \tfrac{i\hbar}{x_0\sqrt2}(\hat a^{\dagger} - \hat a), [a^,a^]=1[\hat a, \hat a^{\dagger}] = 1, [N,a^]=a^[N, \hat a] = -\hat a y [N,a^]=+a^[N, \hat a^{\dagger}] = +\hat a^{\dagger} con N = a†a. Del artículo 03 de este módulo, σAσB12[A^,B^]\sigma_A\sigma_B \ge \tfrac12|\langle[\hat A,\hat B]\rangle| y su versión fina, que añade el término del anticonmutador.

(a) Calcula [N,x^][N, \hat x] usando sólo los dos conmutadores de arriba, sin escribir ninguna matriz, y déjalo en función de p̂. ¿Es hermítico ese conmutador? Deduce de ahí en qué eje del plano complejo vive ⟨[N, x̂]⟩, y por qué la cota de Robertson lleva un módulo y un medio. Escribe la cota en la forma más corta que puedas: sale una expresión con ⟨p̂⟩ dentro. (b) En el autoestado |n⟩: da σN\sigma_N, da σx\sigma_x en nanómetros para n = 0, 1, 2, 3, 4 y 5, y da la cota. Robertson queda 0 ≥ 0. El ejercicio 2 del artículo 03 separó la saturación trivial —cota cero y producto cero— de la cota que no prohíbe nada: ¿cuál de las dos es ésta, y qué habría que enseñar para convencer a alguien de que la desigualdad aquí no dice nada? (c) El estado (0+i1)/2(|0\rangle + i|1\rangle)/\sqrt2, que es la mezcla al 50 % un cuarto de periodo después de (0+1)/2(|0\rangle+|1\rangle)/\sqrt2. Calcula σN\sigma_N, σx\sigma_x en nm, la cota en nm y el producto en nm. ¿Cuánto le sobra al producto sobre la cota? (d) Demuestra que ese sobrante es el más pequeño de todo el ciclo: escribe la razón producto/cota en función de ωt en forma cerrada y minimízala. ¿Se satura Robertson en algún instante? La desigualdad tiene dos maneras independientes de no ser tensa —que Cauchy-Schwarz no sature y que el anticonmutador no se anule—: en el instante de (c), calcula cuál de las dos está actuando, con el número que lo demuestre. El ejercicio 6 del artículo 03 explicó sin calcular por qué x̂ no satura en un oscilador; aquí toca ponerle la cifra.

Pista

Todo (a) sale de escribir x̂ como suma de a y a† y aplicar los dos conmutadores dados por separado: no hace falta ninguna matriz. El resultado se parece muchísimo a p̂, y ese parecido es el problema entero. En (d), la razón sólo depende de ωt a través de sen ωt, y conviene escribirla antes de evaluarla en ningún sitio.

Ojo con (b): σN = 0 en |n⟩ hace que el producto sea cero y que la cota lo sea, y las dos cosas hay que decirlas. Y ojo con el estado de (c): la fase i no cambia |cn|², así que σN es la misma que sin ella — lo que cambia es ⟨p̂⟩, que en el estado real valía cero. Ése es el motivo de que el enunciado ponga la fase y no la olvide.

Comprobación: en (c), el producto y la cota tienen que salirte los dos proporcionales a x0, y su razón tiene que ser un número adimensional sin rastro de la ω ni de ħ. En (d), evalúa tu forma cerrada en ωt = π/2 y compárala con lo que te dio (c): si no coinciden, has puesto un coseno donde va un seno.

Problema 4 ●●○ · Una caja, no una gaussiana

El ejercicio 4 del artículo 01 recorrió el panel con la gaussiana desplazada, que cierra en trece términos. Aquí se desarrolla en la misma base la función más sencilla que se puede dibujar, y no cierra: converge a 1, converge despacio, y su energía media no converge.

Datos. En la variable adimensional ξ = x/x0 del oscilador —el mismo electrón del problema 3, x0 = 1,000 nm, ħω = 76,20 meV—, la caja

f(ξ)=12 si ξ<1,f(ξ)=0 fuera,f(\xi) = \tfrac{1}{\sqrt2}\ \text{si } |\xi| < 1, \qquad f(\xi) = 0 \ \text{fuera},

con ψ0=π1/4eξ2/2\psi_0 = \pi^{-1/4}e^{-\xi^{2}/2}, ψ1=2ξψ0\psi_1 = \sqrt2\,\xi\,\psi_0 y ψ2=(2ξ21)eξ2/2/2π\psi_2 = (2\xi^{2}-1)e^{-\xi^{2}/2}/\sqrt{2\sqrt\pi}. Del artículo 01, para comparar: la gaussiana desplazada un x0 tiene |c0|² = 0,6065 y |c1|² = 0,3033, y su suma llega a 1 en n = 12, con un resto de 1,2329 × 10⁻¹⁴ — que es el término de comparación que (d) necesita. Del panel, que hace falta en (d): corta en N = 12, integra sobre una rejilla de 401 nodos en ξ ∈ [−8, 8] —o sea dξ = 0,04— sumando ψnfdξ\psi_n f\,d\xi nodo a nodo, y suma la cola hasta n = 40 en vez de restarla de 1. Y un detalle suyo que hay que reproducir para contestar (d): normaliza f sobre esa misma rejilla, dividiendo por ifi2dξ\sqrt{\sum_i f_i^{2}\,d\xi}, no por la norma analítica.

(a) Comprueba que f está normalizada. Di cuánto vale c1 sin integrar nada, y por qué. Calcula c0 por cuadratura y en forma cerrada —sale una función error—, y da |c0|² y cuánto le falta a 1 con n ≤ 1. (b) Sigue sumando. ¿Converge Σ|cn|² a 1? ¿A qué ritmo? Ajusta el exponente de cn2nα|c_n|^{2}\propto n^{-\alpha} sobre los pares, di de qué rasgo de f sale ese exponente —no es la anchura— y cuenta cuántos términos hacen falta para llegar al 99 % y al 99,5 %. Compara ξ4/ξ22\langle\xi^{4}\rangle/\langle\xi^{2}\rangle^{2} de la caja con el 3 de cualquier gaussiana y di qué te dice ese número sobre la forma. (c) Ahora la energía, H=ωncn2(n+12)\langle H\rangle = \hbar\omega\sum_n |c_n|^{2}(n+\tfrac12). Da la suma parcial hasta n = 12 y hasta n = 4000, en meV, y di qué está pasando y con qué ley crece. Demuéstralo sin una sola función de Hermite: suaviza el borde de la caja con una rampa de anchura w y calcula ⟨p̂²⟩ en función de w. Enuncia en una frase la conclusión —hay dos conjuntos de funciones y f está en uno y no en el otro— y da además σx de la caja en nm: ¿es más ancha o más estrecha que el estado fundamental? (d) El panel del artículo 01 con esta f. Si se le pudiera cambiar la función y nada más: ¿qué marcaría su casilla de Σ con el corte en N = 12, y qué concluiría un lector que sólo mirase esa casilla? ¿Y qué c0 daría su cuadratura de 401 nodos, y con qué error relativo frente al exacto de (a)? Los dos números son defectos de cosas distintas: di de qué es defecto cada uno, y cuál de los dos se arregla moviendo un mando.

Pista

En (a), 11eξ2/2dξ\int_{-1}^{1} e^{-\xi^{2}/2}d\xi es la única integral del problema y ya tiene nombre. En (b), la ley de potencias se ve mejor en escala doble logarítmica y sobre los pares solos, porque los impares son todos cero. En (c), la rampa convierte una derivada que no existe en una que vale 1/w en un tramo de anchura w, y el cuadrado de eso integrado sobre el tramo es lo que se dispara.

Ojo con el orden de las dos conclusiones de (b) y (c): que Σ|cn|² converja a 1 no autoriza a sumar la serie de la energía término a término y esperar un número. Son dos series distintas, con el mismo |cn|² y pesos distintos. Y ojo en (d): uno de los dos números que se piden depende del corte y el otro no —depende de dónde caen los nodos respecto del salto de f—, y confundirlos lleva a «subiendo N se arregla», que es falso para uno de los dos.

Comprobación cruzada, y es bonita: la constante que gobierna el resto de Parseval y la constante que gobierna el crecimiento de ⟨H⟩ tienen que ser la misma, porque las dos salen del mismo exponente. Ajusta una y predice la otra; si no coinciden en las dos primeras cifras, tu exponente está mal. Y para el c0 del panel son dos pasos, no uno: cuenta primero qué nodos de su rejilla recogen algo —el salto de f cae justo encima de uno, y la suma deja fuera una franja de anchura dξ/2 a cada lado—, y acuérdate después de que el panel divide por la norma que él mismo mide sobre esa rejilla, que con la caja ya no vale 1. Los dos efectos tienen el mismo origen y se cancelan en parte: si te saltas el segundo paso, tu discrepancia sale unas tres veces mayor que la que el panel enseña de verdad. Y si te sale en el sexto decimal, has integrado la caja como si fuera continua.

Problema 5 ●●● · La cota es del estado; saturarla, del observable

El ejercicio 6 del artículo 03 preparó un oscilador en la mezcla al 50 % de sus dos primeros niveles y encontró un observable que satura Mandelstam-Tamm en todo instante: el pseudoespín 01+10|0\rangle\langle1| + |1\rangle\langle0|, que dentro de ese plano es σ^x\hat\sigma_x. Aquí se rompe la simetría de los pesos, y con ella se mueve el observable que satura.

Datos. El mismo electrón de los problemas 3 y 4: x0 = 1,000 nm, ħω = 76,20 meV, k = 12,21 mN/m, En=ω(n+12)E_n = \hbar\omega(n+\tfrac12). El estado es

ψ=130+231.|\psi\rangle = \sqrt{\tfrac13}\,|0\rangle + \sqrt{\tfrac23}\,|1\rangle .

Del artículo 03, τBσB/dB/dt/2σE\tau_B \equiv \sigma_B/|d\langle B\rangle/dt| \ge \hbar/2\sigma_E, la condición de saturación en un qubitσ^xσ^z=0\langle\hat\sigma_x\rangle\langle\hat\sigma_z\rangle = 0 para la pareja (σ^z,σ^x)(\hat\sigma_z, \hat\sigma_x)— y, para comparar en (c) y en (d), la mezcla al 50 % de los dos primeros niveles del pozo de 1 nm: ΔE = 1,128 eV, σE = 0,564 eV, cota 0,583 fs, batido T = h/ΔE = 3,666 fs = 2πτ y un mínimo de τx de 1,262 cotas. En todo el módulo ħ = 0,658212 eV·fs y h = 4,13567 × 10⁻¹⁵ eV·s.

(a) Calcula ⟨H⟩ y la dispersión de energía σE en meV, por dos caminos que no compartan cuenta —la media pesada de los dos niveles y H2H2\sqrt{\langle H^{2}\rangle - \langle H\rangle^{2}}—, y da la cota de Mandelstam-Tamm en femtosegundos. Escribe también ⟨H⟩/ħω y σE/ħω como fracciones exactas: las dos son bonitas y la segunda lleva una raíz. (b) Dentro del plano de |0⟩ y |1⟩ esto es un qubit con H^=Eˉ+ω2σ^z\hat H = \bar E + \tfrac{\hbar\omega}{2}\hat\sigma_z. Sitúa el estado en la esfera de Bloch en t = 0 y comprueba que el pseudoespín del artículo 03 —el σ^x\hat\sigma_x del plano— ya no satura. Encuentra el observable del plano que sí lo hace en t = 0, escríbelo en la base de |0⟩ y |1⟩ e identifícalo con un operador físico del oscilador: no es la posición, y tampoco es exactamente el momento. Di qué le falta al momento de verdad para ser ése. (c) Ponle cifra a lo que el artículo 03 dejó dicho de palabra. Calcula el mínimo sobre el ciclo de τx=σx/dx/dt\tau_x = \sigma_x/|d\langle x\rangle/dt| en femtosegundos y en unidades de la cota, y di qué porcentaje le sobra a la posición. Explica el mecanismo con un número, y sin derivar nada: descompón x̂|ψ⟩ en la parte que vive en el plano de |0⟩ y |1⟩ y la que se ha ido al nivel n = 2, y saca de ahí el mismo cociente que acabas de medir. ¿Le pasa lo mismo al momento entero, y por qué? (d) Compara la cota de (a) con los 0,583 fs de la mezcla del pozo. ¿Cuál es más larga, por qué factor, y qué magnitud fija ese factor? Un lector escribe «el pozo va más deprisa porque es más pequeño»: corrígelo con la magnitud que de verdad manda, y di qué tendría que cambiarse del oscilador —y en qué proporción— para igualar las dos cotas.

Pista

Para una mezcla de dos niveles con pesos p0 y p1, σE=p0p1ΔE\sigma_E = \sqrt{p_0p_1}\,\Delta E, y conviene comprobarlo con la definición antes de usarlo. En (b), la condición de saturación del artículo 03 se anula por uno de sus dos factores; mira cuál puede anularse aquí, porque el otro ya no.

Ojo con el reflejo de copiar el 2π del artículo 03: ese cociente entre el periodo y la cota vale 2π sólo en la mezcla al 50 %, y ésta no lo es. Y ojo en (c): σx de este estado no es la del fundamental ni la de |1⟩, porque el término cruzado entre los dos niveles no se anula — es justamente lo que hace que ⟨x⟩ se mueva.

Comprobación: el mínimo de τx tiene que salirte por encima de la cota, nunca por debajo; si te sale menor, has derivado ⟨x⟩ respecto de ωt y no de t. Y en (b) hay una comprobación que se hace sin calcular nada: en t = 0 el valor esperado del pseudoespín del artículo 03 está en un extremo de su oscilación, así que su ritmo es cero — mira qué le hace eso a su τB y tendrás, de un vistazo, lo lejos que está de saturar. El observable que encuentres tiene que darte, en t = 0, la cota de (a) exacta y no aproximada.

Problema 6 ●●● · Cuánto destruye una medida, y cómo se le da la vuelta

El ejercicio 5 del artículo 04 midió σz y luego σx sobre una superposición, y encontró cuatro parejas al 25 %. Ese caso es un extremo: aquí se pide la familia entera de la que es el extremo, y después lo que ocurre cuando el extremo se recorre a pasos pequeños.

Datos. Un qubit preparado en |0⟩, autoestado de σz con autovalor +1. Se mide primero

σ^θ=cosθσ^z+sinθσ^x,\hat\sigma_\theta = \cos\theta\,\hat\sigma_z + \sin\theta\,\hat\sigma_x ,

sin mirar el resultado, y después σz. Del artículo 04, el postulado III —p(a)=Paψ2p(a) = \|P_a|\psi\rangle\|^{2}— y el postulado IV, que deja el estado en Paψ/PaψP_a|\psi\rangle/\|P_a|\psi\rangle\|. Para el apartado (b), un hamiltoniano H^=(ε/2)σ^x\hat H = (\varepsilon/2)\hat\sigma_x y otro H^=(ε/2)σ^y\hat H = (\varepsilon/2)\hat\sigma_y, aplicados entre las dos medidas durante un tiempo t.

(a) Calcula P(z = +1) al final, en función de θ, aplicando IV a cada rama y sumando las dos contribuciones. Evalúala en θ = 30°, 45°, 60° y 75°. ¿Dónde es máximo el daño, y cuánto vale? Comprueba que θ = 0° y θ = 90° reproducen los dos casos que el artículo 04 ya tiene, y di cuál es cuál. (b) Con θ = 90°, deja pasar un tiempo t entre las dos medidas bajo H ∝ σx. ¿Cambia P(z = +1)? Ahora bajo H ∝ σy: calcula las dos probabilidades condicionadas a cada resultado de la primera medida, en función de φ = εt/ħ, y vuelve a sumar. La respuesta débil a este apartado es «no cambia porque el estado es autoestado de H»; la fuerte es otra y es mucho más fuerte. Enúnciala, y di qué la protege — qué propiedad de las dos ramas hace que sus inclinaciones se cancelen exactamente. (c) La vuelta de tuerca. En vez de una medida a 90°, haz N medidas escalonadas, la k-ésima a θk = kπ/2N, y llévate sólo la cadena en la que todas dan +1. Calcula la probabilidad de esa cadena en forma cerrada y evalúala para N = 1, 2, 5, 10 y 100. ¿Hacia dónde tiende con N → ∞, y en qué estado acaba el sistema? Esto tiene nombre, y el glosario lo trae con su versión óptica —una pila de polarizadores—. El factor que separa las dos versiones es el medio ángulo de la esfera de Bloch: un giro de Δθ en la esfera mueve el estado Δθ/2. Con eso como dato, comprueba tu fórmula contra los tres porcentajes ópticos que el glosario publica y explica por qué N = 2 no da el 25 % que dan dos polarizadores cruzados por la mitad. (d) Con las tres partes delante, escribe qué postulado es el responsable de cada cosa: del daño de (a), de la invariancia de (b) y del arrastre de (c). Y contesta a esto: ¿es «medir perturba» un enunciado del formalismo? Si lo es, di dónde está escrito; si no lo es, di qué lo sustituye y qué predice de más.

Pista

En (a) hay que sumar sobre las dos ramas de la primera medida —es una medida no selectiva— y cada rama contribuye con el producto de dos probabilidades de Born. Sale un polinomio en cos θ de grado dos. En (c), el paso entre dos medidas consecutivas es siempre el mismo ángulo, así que la probabilidad de la cadena es una potencia N-ésima.

Ojo con el factor 2 de (c): el enunciado te da qué pasa —el estado lleva medio ángulo— y lo que se pide es de dónde sale. Míralo en la parametrización del qubit: |0⟩ y |1⟩ son ortogonales y en la esfera están a 180°, así que la esfera no puede estar midiendo el ángulo del espacio de Hilbert. Escribe el solape de dos estados cuyos puntos de Bloch difieren en Δθ antes de elevar nada al cuadrado, y verás en qué renglón entra el 2 y por qué la luz, que no lleva medio ángulo, da otro número. En (b), la trampa es contestar con el caso particular que el enunciado te pone delante: prueba el segundo hamiltoniano antes de generalizar, y fíjate en que las condicionales sí se mueven.

Comprobación: tu fórmula de (a) tiene que valer 1 en θ = 0 y coincidir en θ = 90° con la suma de dos de las cuatro parejas del ejercicio 5 del artículo 04. En (c), tu forma cerrada con N = 1 tiene que devolverte exactamente el caso θ = 90° de (a): si no, has contado mal los pasos o has puesto el ángulo entero donde va el medio.

Fin del módulo II.4. Lo que se ha ganado tiene un nombre preciso: los tres módulos anteriores usaban un formalismo, y ahora está escrito. El producto escalar es una definición con tres axiomas y no una integral, y de ellos salen Cauchy-Schwarz y Parseval, que convierten «esta base es completa» en una suma que se mide. Un observable es un operador hermítico, y se sabe qué compra esa hipótesis —autovalores reales— y qué no compra, porque hay contraejemplos con número. La relación de Robertson se deduce en cinco líneas de Cauchy-Schwarz, y con ella el ħ/2 de Kennard y la cota de Mandelstam-Tamm que el II.1 usó a crédito dejan de ser fórmulas prestadas. Y los cinco postulados están escritos seguidos, con la parte incómoda señalada y medida en vez de disimulada.

Y esta hoja ha entrenado lo que de verdad se transfiere: qué cuenta decide qué. Ninguno de los seis problemas se resuelve por parecido. La ortogonalidad se decide con un conmutador y no resolviendo un problema de autovalores; que una tercia sea base se decide con Parseval y no con un determinante; que una función se pueda desarrollar y que tenga energía media son dos preguntas distintas sobre la misma serie; y que una medida «perturbe» resulta ser una función continua de un ángulo, con un extremo en el que arrastra el estado entero sin romperlo.

Lo que falta. Todo el teorema espectral de aquí es de dimensión finita o de funciones bien portadas: el caso infinito con dominios, la diferencia entre simétrico y autoadjunto —que el artículo 02 deja abierta a propósito, y que el problema 1 no puede cerrar por mucho que discretice—, el espectro continuo y el teorema de Gleason son el III.1. La medida generalizada, cuando el aparato no es un proyector, tampoco está. Y la degeneración se ha tratado como posibilidad formal, sin un solo sistema de este módulo donde ocurra de verdad.

Ésa es justamente la deuda que manda, y por eso el siguiente es el II.5, «Cuántica en tres dimensiones: el átomo de hidrógeno». Al pasar a tres dimensiones el potencial de Coulomb entrega un espectro con degeneración de verdad —varios estados con la misma energía— y entonces el ejercicio 6 del artículo 02 deja de ser una curiosidad de laboratorio: hará falta un conjunto completo de observables que conmuten, tres de ellos, para nombrar un estado sin ambigüedad, y ahí es donde los números cuánticos dejan de ser una convención de notación. El resumen en frío del artículo 04 está escrito para que ese módulo se pueda empezar sin volver a leer éste.