Módulo II.6 · Artículo 04

Suma de momentos

Dos espines ½ dan un triplete y un singlete separados por exactamente 1 ħ² de ⟨S₁·S₂⟩: +¼ arriba y −¾ abajo, por dos caminos que no comparten cuenta. Ese único ħ², metido en el fundamental del hidrógeno con el término de contacto de Fermi, sale convertido en 1422,808 MHz; corregir el radio de Bohr por la masa reducida —y a₀ entra al cubo— lo baja a 1420,485, a 56 ppm de la raya de 21 cm medida. De paso, ℓ ⊗ ½ se resuelve bajando con J₋ y sin abrir ninguna tabla. Cierra el módulo con su resumen en frío.

Entre los dos niveles hiperfinos del hidrógeno hay exactamente 1 ħ²: S1S2\langle \mathbf{S}_1\cdot\mathbf{S}_2\rangle vale +0,250 ħ² en el triplete y −0,750 ħ² en el singlete, y la diferencia es 1 ħ² sin redondear nada. Este artículo hace dos cosas con ese número. La primera es deducirlo dos veces por caminos que no comparten ninguna cuenta —despejándolo de S2\mathbf{S}^2 y sumando las tres componentes S1kS2k\langle S_{1k}S_{2k}\rangle— y ver de paso que S1zS2z\langle S_{1z}S_{2z}\rangle vale −0,250 ħ² en los dos estados: el eje z no los distingue. La segunda es cobrárselo. Metido en el fundamental del hidrógeno con el término de contacto de Fermi, ese 1 ħ² se convierte en 1422,808 MHz; corregir el radio de Bohr por la masa reducida —y entra al cubo, no al cuadrado— lo deja en 1420,485 MHz, a 56 ppm de los 1420,405 751 7667 MHz que se miden en el cielo. Los 56 ppm que quedan no son redondeo: son relatividad, QED y estructura del protón, y se dicen con su nombre.

Necesitas: este módulo entero, y de él tres cosas concretas. Del artículo 01, el operador escalera J±=Jx±iJyJ_\pm = J_x \pm iJ_y con su coeficiente j(j+1)m(m±1)\hbar\sqrt{j(j+1)-m(m\pm1)}, que aquí se usa en cada deducción. Del artículo 02, que el espín ½ se representa con las matrices de Pauli y que S=σ/2\mathbf{S} = \hbar\boldsymbol{\sigma}/2. Del artículo 03, los dos factores gg y las dos escalas magnéticas — geg_e, gpg_p, μB y μN—, que son los cuatro números que entran en la raya de 21 cm. De fuera del módulo hacen falta dos: del II.5 artículo 04, que el nivel nn tiene n2n^2 orbitales; y del I.4 artículo 01, el valor ψ100(0)2=2,15×1030|\psi_{100}(0)|^2 = 2{,}15\times10^{30} m⁻³ y que la raya de 21 cm es proporcional a él. Lo que NO hace falta: la desigualdad de Bell. El singlete aparece aquí como estado acoplado, no como argumento sobre el realismo local; eso es el I.5 artículo 03.
Dos espines ½: el triplete, el singlete y ⟨S₁·S₂⟩

Cuatro estados producto —|↑↑⟩, |↑↓⟩, |↓↑⟩, |↓↓⟩— y dos etiquetas nuevas, S y M. Los de M = ±1 no tienen elección: son |1,+1⟩ y |1,−1⟩. En el plano M = 0 sí la hay, y el mando la recorre: |ψ⟩ = cos χ |↑↓⟩ + sen χ |↓↑⟩. En χ = +45° las dos amplitudes son iguales y sale |1,0⟩; en χ = −45° son opuestas y sale |0,0⟩; en χ = 0 el estado es el producto |↑↓⟩, que no es autoestado de y por eso reparte. El mando llega exactamente a los dos extremos con Inicio y Fin.

M =
pesosS = 1 · tripletesinglete · S = 0⟨S₁·S₂⟩−¾−½−¼0autovalor del singleteautovalor del tripleteCzzCxxCyy−¼ ħ²0+¼ ħ²

El dibujo sale de los mismos números que las casillas de abajo, no de una cuenta aparte. Y no afirma nada por su cuenta: un píxel de la barra de pesos vale unas dos décimas de punto porcentual, así que en χ = +44° el trozo de singlete es más fino que un píxel y aun así no es cero — el peso exacto es siempre la cifra impresa. Las dos marcas de la regla son los dos caminos a ⟨S₁·S₂⟩: el punto es la suma de las tres correlaciones y el anillo es (S²−S₁²−S₂²)/2. Si alguna vez se separan, una de las dos cuentas está mal.

P(triplete) 50.0 % S = 1
P(singlete) 50.0 % S = 0
⟨S₁·S₂⟩ -0.2500 en ħ² · sólo +¼ y −¾ son autovalores
Czz -0.2500 ⟨S₁zS₂z
Cxx 0.0000 ⟨S₁xS₂x
Cyy 0.0000 ⟨S₁yS₂y
(S²−S₁²−S₂²)/2 -0.2500 el otro camino a ⟨S₁·S₂⟩

χ = 0°, o sea el producto |↑↓⟩: es autoestado de SzM = 0 — y no de . Reparte 50.0 % y 50.0 %: aquí la mitad y la mitad es literal, porque las dos amplitudes son +0.7071 y +0 en la base producto y eso da ±0.7071 en la base |S, M⟩. ⟨S₁·S₂⟩ = -0.2500 ħ², que es la media de +¼ y −¾ pesada a partes iguales — y no es autovalor de nada.

Censo que hace este panel al cargarse, no cableado: de los 93 estados que los mandos alcanzan —91 posiciones del deslizador más los dos sectores M = ±1—, 3 son triplete puro (los dos de M = ±1 y el χ = +45°), 1 es singlete puro y 1 reparte exactamente mitad y mitad (el χ = 0, o sea |↑↓⟩). Ninguna otra posición recibe esas etiquetas, y ningún rótulo de arriba imprime «100 %» ni «0 %» sobre un peso que no lo sea: en χ = +44° el panel escribe 99,97 %, que es lo que hay. Los dos caminos a ⟨S₁·S₂⟩ —las tres correlaciones por un lado, el peso del triplete por otro— se separan como mucho 2.2e-16 en los 93 estados, y la probabilidad cierra: P(triplete) + P(singlete) vale 1 exacto en 58 de los 93 y se aparta como mucho 2.2e-16 en el resto. Los pesos son cocientes de Rayleigh, ⟨ψ|A|ψ⟩/⟨ψ|ψ⟩, y por eso los «100 %» y los «0 %» de los cuatro estados con nombre son el peso y no el redondeo.

Sumar no es apilar: los números cuánticos buenos cambian

Dos electrones, cada uno con su espín. El espacio de estados es el producto tensorial de los dos, así que tiene 2×2=42\times2 = 4 dimensiones, y la base más obvia es la de los cuatro estados producto:  ⁣|\!\uparrow\uparrow\rangle,  ⁣|\!\uparrow\downarrow\rangle,  ⁣|\!\downarrow\uparrow\rangle y  ⁣|\!\downarrow\downarrow\rangle. Cada uno de ellos tiene m1m_1 y m2m_2 perfectamente definidos. Parece que no hay nada más que decir.

Lo hay, y es lo que decide el resto del artículo. El operador que gobierna la energía de dos espines que interactúan casi nunca es S1zS_{1z} ni S2zS_{2z}: es S1S2\mathbf{S}_1\cdot\mathbf{S}_2, porque la interacción no distingue direcciones del espacio. Y ese operador no conmuta con S1zS_{1z}. Es fácil de ver sin calcular nada: S12\mathbf{S}_1^2 y S22\mathbf{S}_2^2 son múltiplos de la identidad en este problema, así que [S1S2,S1z]=12[S2,S1z][\mathbf{S}_1\cdot\mathbf{S}_2, S_{1z}] = \tfrac12[\mathbf{S}^2, S_{1z}], y eso no es cero. El verificador lo comprueba en las nueve series de acoplamiento que recorre: en todas, J2\mathbf{J}^2 conmuta con JzJ_z y no conmuta con J1zJ_{1z}. Traducido: si un hamiltoniano contiene S1S2\mathbf{S}_1\cdot\mathbf{S}_2, entonces m1m_1 y m2m_2 dejan de ser números cuánticos buenos —dejan de conservarse— y hay que cambiar de base.

La base nueva la fabrica el momento angular total J=J1+J2\mathbf{J} = \mathbf{J}_1 + \mathbf{J}_2. El primer resultado es puramente algebraico y no cuesta nada: la suma de dos momentos angulares vuelve a cumplir el álgebra del artículo 01, porque J1\mathbf{J}_1 y J2\mathbf{J}_2 actúan sobre factores distintos del producto y por tanto conmutan entre sí. De ahí que J2\mathbf{J}^2 y JzJ_z vuelvan a tener autovalores 2J(J+1)\hbar^2 J(J+1) y M\hbar M, con M=JJM = -J\ldots J. Lo que no es gratis es saber qué valores de JJ aparecen. La respuesta es la serie de Clebsch-Gordan:

j1j2=j1j2    j1j2+1        j1+j2,j_1 \otimes j_2 = |j_1-j_2| \;\oplus\; |j_1-j_2|+1 \;\oplus\; \cdots \;\oplus\; j_1+j_2 ,

y su primera comprobación no es una cita, es una cuenta: las dimensiones tienen que cuadrar. Comprobado en nueve series, desde 1212\tfrac12\otimes\tfrac12 hasta 522\tfrac52\otimes2, diagonalizando J2\mathbf{J}^2 en la base producto. En las nueve, (2j1+1)(2j2+1)=J(2J+1)(2j_1+1)(2j_2+1) = \sum_J (2J+1): 4 = 3+1 para dos espines ½, 6 = 4+2 para 1121\otimes\tfrac12, 9 = 5+3+1 para 111\otimes1, 30 = 2+4+6+8+10 para 522\tfrac52\otimes2. Ninguna de esas igualdades es una identidad aritmética que se pueda escribir a ojo: es el mismo número de estados contado en dos bases.

½ ⊗ ½: la escalera se cierra y sobra un estado

El caso de dos espines ½ se resuelve entero con el operador escalera del artículo 01, sin abrir ninguna tabla y sin diagonalizar ninguna matriz. Se empieza por arriba. El estado con M=+1M = +1 es único —sólo  ⁣|\!\uparrow\uparrow\rangle tiene m1+m2=+1m_1+m_2 = +1— así que tiene que ser el tope de un multiplete, y como MM llega a 1 ese multiplete es S=1S=1:

1,+1= ⁣.|1,+1\rangle = |\!\uparrow\uparrow\rangle .

Ahora se baja. S=S1+S2S_- = S_{1-} + S_{2-} actúa sobre los dos factores, y sobre un espín ½ vale S ⁣= ⁣S_-|\!\uparrow\rangle = \hbar|\!\downarrow\rangle. Por tanto S ⁣=( ⁣+ ⁣)S_-|\!\uparrow\uparrow\rangle = \hbar(|\!\downarrow\uparrow\rangle + |\!\uparrow\downarrow\rangle). Por el otro lado, el coeficiente de escalera del artículo 01 dice que S1,+1=12101,0=21,0S_-|1,+1\rangle = \hbar\sqrt{1\cdot2 - 1\cdot0}\,|1,0\rangle = \hbar\sqrt2\,|1,0\rangle, y 2=1,4142136\sqrt2 = 1{,}4142136 es exactamente la norma que se mide del vector de la izquierda antes de normalizarlo. Igualando:

1,0= ⁣+ ⁣2,1,1= ⁣.|1,0\rangle = \frac{|\!\uparrow\downarrow\rangle + |\!\downarrow\uparrow\rangle}{\sqrt2}, \qquad |1,-1\rangle = |\!\downarrow\downarrow\rangle .

Los dos coeficientes valen +0,70710678, con el mismo signo. Eso deja tres estados colocados y uno sin colocar: en el plano M=0M=0 hay dos estados producto y sólo uno de ellos se ha usado. El que falta es lo ortogonal a 1,0|1,0\rangle dentro de ese plano, y con M=0M=0 y sin sitio en el multiplete S=1S=1 tiene que ser S=0S=0:

0,0= ⁣ ⁣2.|0,0\rangle = \frac{|\!\uparrow\downarrow\rangle - |\!\downarrow\uparrow\rangle}{\sqrt2} .

El menos va en el singlete. Es el error más frecuente del tema y no chirría al escribirlo, así que conviene tener tres razones independientes para el signo, y las tres están medidas. La primera es la construcción de arriba: el triplete sale de bajar, y bajar suma los dos términos con el mismo signo. La segunda es la diagonalización directa de S2\mathbf{S}^2: sus cuatro autovalores son 222\hbar^2 tres veces y 00 una, y el autovector del 0 tiene componentes +0,707107 y −0,707107 sobre  ⁣|\!\uparrow\downarrow\rangle y  ⁣|\!\downarrow\uparrow\rangle, y exactamente 0 sobre los otros dos. La tercera es la tabla de Clebsch-Gordan de sympy, que con el convenio de Condon-Shortley da 12,+12;12,120,0=+0,70710678\langle\tfrac12,+\tfrac12;\tfrac12,-\tfrac12|0,0\rangle = +0{,}70710678 y 12,12;12,+120,0=0,70710678\langle\tfrac12,-\tfrac12;\tfrac12,+\tfrac12|0,0\rangle = -0{,}70710678. Los tres caminos coinciden en el signo.

Y hay una cuarta razón, que no es de cálculo sino de física, y que es la que hay que recordar: bajo el intercambio de las dos partículas el triplete es simétrico —sus tres estados— y el singlete es antisimétrico. Es la única combinación de signos que puede producir eso, porque  ⁣|\!\uparrow\uparrow\rangle y  ⁣|\!\downarrow\downarrow\rangle son simétricos sin remedio y el multiplete entero tiene que serlo con ellos. De esa distinción vive el II.7, el módulo de partículas idénticas: con espín ½, la antisimetría total obliga a que un singlete de espín lleve una parte espacial simétrica y un triplete una antisimétrica, y ahí es donde nace la interacción de intercambio.

⟨S₁·S₂⟩ por dos caminos, y el único ħ² que hay entre los dos niveles

El observable que interesa no es S2\mathbf{S}^2 sino S1S2\mathbf{S}_1\cdot\mathbf{S}_2, porque es él quien aparece en el hamiltoniano. Se despeja del cuadrado de la suma:

S2=(S1+S2)2=S12+S22+2S1S2S1S2=S2S12S222.\mathbf{S}^2 = (\mathbf{S}_1+\mathbf{S}_2)^2 = \mathbf{S}_1^2 + \mathbf{S}_2^2 + 2\,\mathbf{S}_1\cdot\mathbf{S}_2 \quad\Longrightarrow\quad \mathbf{S}_1\cdot\mathbf{S}_2 = \frac{\mathbf{S}^2 - \mathbf{S}_1^2 - \mathbf{S}_2^2}{2} .

Con S12=S22=342\mathbf{S}_1^2 = \mathbf{S}_2^2 = \tfrac34\hbar^2 —que aquí no se teclea: sale de las matrices j=12j=\tfrac12 construidas por la escalera— los dos autovalores son inmediatos. En el triplete S2=22\mathbf{S}^2 = 2\hbar^2, luego

S1S2S=1=2343422=+24=+0,250 2,\langle \mathbf{S}_1\cdot\mathbf{S}_2\rangle_{S=1} = \frac{2 - \tfrac34 - \tfrac34}{2}\,\hbar^2 = +\frac{\hbar^2}{4} = +0{,}250\ \hbar^2 ,

y en el singlete S2=0\mathbf{S}^2 = 0, luego

S1S2S=0=0343422=324=0,750 2.\langle \mathbf{S}_1\cdot\mathbf{S}_2\rangle_{S=0} = \frac{0 - \tfrac34 - \tfrac34}{2}\,\hbar^2 = -\frac{3\hbar^2}{4} = -0{,}750\ \hbar^2 .

Ese 34-\tfrac34 es el número que más veces se escribe mal. El error habitual es poner 2/4-\hbar^2/4, que es un número perfectamente creíble, tiene el signo bueno y existe en este problema —es S1zS2z\langle S_{1z}S_{2z}\rangle, que aparece en la sección siguiente—, pero no es esto. La diferencia entre los dos autovalores es 14(34)=1\tfrac14 - (-\tfrac34) = 1, exactamente 1 ħ², y de ese 1 sale el tamaño entero de la raya de 21 cm. Con el −¼ saldría cuatro veces menor.

El segundo camino no comparte nada con el primero: se escriben las tres componentes y se suman. S1S2=S1xS2x+S1yS2y+S1zS2z\mathbf{S}_1\cdot\mathbf{S}_2 = S_{1x}S_{2x} + S_{1y}S_{2y} + S_{1z}S_{2z}, con cada S1kS2k=24σkσkS_{1k}S_{2k} = \tfrac{\hbar^2}{4}\,\sigma_k\otimes\sigma_k. Contraídos los tres productos de Kronecker sobre los autovectores de S2\mathbf{S}^2 y sumados, los dos caminos dan lo mismo, y el panel de arriba imprime los dos resultados en dos casillas distintas para que se vean moverse juntos.

Las tres correlaciones, y el eje que no distingue

La suma de las tres componentes esconde algo que la suma sola no deja ver, y es la parte del artículo que más se usa después. Llamemos Ckk=S1kS2kC_{kk} = \langle S_{1k}S_{2k}\rangle a cada correlación por separado, en unidades de 2\hbar^2. Para los dos estados de M=0M=0 salen así:

EstadoCzzCxxCyySuma = ⟨S₁·S₂⟩
|1,0⟩ · triplete−0,250+0,250+0,250+0,250
|0,0⟩ · singlete−0,250−0,250−0,250−0,750

Mírese la primera columna. CzzC_{zz} vale −0,250 ħ² en los dos estados. No se parece: es el mismo número. Un experimento que mida S1zS_{1z} y S2zS_{2z} a la vez y calcule la correlación de los dos resultados no distingue el triplete del singlete, por muchas veces que se repita. La razón es geométrica y se ve en la construcción: los dos estados son combinaciones de los mismos dos productos  ⁣|\!\uparrow\downarrow\rangle y  ⁣|\!\downarrow\uparrow\rangle, y los dos tienen los espines antiparalelos en z con probabilidad 1. Lo único que los separa es el signo relativo, y el signo relativo no lo ve un observable diagonal en esa base.

Lo que sí lo ve son CxxC_{xx} y CyyC_{yy}, que cambian de signo: +0,250 ħ² en el triplete y −0,250 ħ² en el singlete. En el eje z los dos estados están anticorrelacionados; en los ejes x e y, el triplete está correlacionado y el singlete anticorrelacionado. Ésa es toda la diferencia entre los dos, y es lo que hay que medir para verla. La lección general, que vale mucho más allá de este problema: un observable que devuelve el mismo valor en dos estados no es una mala medida, es una medida de otra cosa.

Y por qué la suma sigue saliendo distinta. Si CzzC_{zz} es igual en los dos y la suma no lo es, el desacuerdo tiene que estar entero en las otras dos, y está: −0,250 + 0,250 + 0,250 = +0,250 frente a −0,250 − 0,250 − 0,250 = −0,750. Conviene también no confundir CzzC_{zz} con S1S2\langle \mathbf{S}_1\cdot\mathbf{S}_2\rangle: en el singlete valen −¼ ħ² y −¾ ħ², un factor 3 de diferencia, y los dos son números de este mismo problema. Escribir el primero donde va el segundo es el error del párrafo anterior visto desde el otro lado.

Ejemplo resuelto 1 · El estado |↑↓⟩, que no es autoestado de nada interesante

Problema. Un electrón apunta arriba y el otro abajo: ψ= ⁣|\psi\rangle = |\!\uparrow\downarrow\rangle. Es el estado que cualquiera dibujaría al decir «espines opuestos». (a) ¿Es autoestado de SzS_z? ¿Y de S2\mathbf{S}^2? (b) Si se mide S2\mathbf{S}^2, ¿qué valores salen y con qué probabilidad? (c) ¿Cuánto vale S1S2\langle \mathbf{S}_1\cdot\mathbf{S}_2\rangle, y es un resultado posible de una medida?

Solución. (a) Sz=S1z+S2zS_z = S_{1z}+S_{2z} da (+22)ψ=0(+\tfrac\hbar2 - \tfrac\hbar2)|\psi\rangle = 0: sí es autoestado, con M=0M=0. Con S2\mathbf{S}^2 no hay tanta suerte. La forma rápida de verlo es invertir las dos definiciones de arriba, que es un sistema lineal de dos por dos:

 ⁣=1,0+0,02, ⁣=1,00,02.|\!\uparrow\downarrow\rangle = \frac{|1,0\rangle + |0,0\rangle}{\sqrt2}, \qquad |\!\downarrow\uparrow\rangle = \frac{|1,0\rangle - |0,0\rangle}{\sqrt2} .

Como 1,0|1,0\rangle y 0,0|0,0\rangle tienen autovalores distintos de S2\mathbf{S}^2222\hbar^2 y 00—, la mezcla no puede ser autoestado de nada.

(b) La regla de Born sobre esa descomposición: 1,0 ⁣2=0,5|\langle 1,0|\!\uparrow\downarrow\rangle|^2 = 0{,}5 y 0,0 ⁣2=0,5|\langle 0,0|\!\uparrow\downarrow\rangle|^2 = 0{,}5, las dos medidas y no supuestas. Medir S2\mathbf{S}^2 devuelve 222\hbar^2 la mitad de las veces y 00 la otra mitad. Nótese que las dos mitades son exactas y no aproximadas: el coeficiente es 1/21/\sqrt2 por construcción.

(c) El valor esperado es la media pesada: 0,5(+14)+0,5(34)=140{,}5\cdot(+\tfrac14) + 0{,}5\cdot(-\tfrac34) = -\tfrac14, o sea −0,250 ħ², que es también lo que da la suma de las tres correlaciones —Czz=0,250C_{zz} = -0{,}250 y Cxx=Cyy=0C_{xx} = C_{yy} = 0 en este estado, como imprime el panel—. Y no es un resultado posible de ninguna medida: los autovalores de S1S2\mathbf{S}_1\cdot\mathbf{S}_2 son +¼ ħ² y −¾ ħ² y no hay más.

Resultado. «Los dos espines apuntan en sentidos opuestos» no describe un estado de espín total definido. Describe una superposición a partes iguales de los dos que sí lo tienen, y esa superposición es la que un hamiltoniano con S1S2\mathbf{S}_1\cdot\mathbf{S}_2 partirá en dos energías distintas en cuanto se encienda. El estado que la intuición clásica considera el más definido de los cuatro es justo el que no tiene definida la magnitud que gobierna su energía — y el 0,250 2-0{,}250\ \hbar^2 del apartado (c) es, además, la trampa: coincide con CzzC_{zz} y con el 2/4-\hbar^2/4 que se escribe por error en el singlete, y los tres son cosas distintas.

ℓ ⊗ ½: los Clebsch del 2p sin abrir una tabla

El caso que hace falta para un átomo real no es 1212\tfrac12\otimes\tfrac12 sino 12\ell\otimes\tfrac12: el momento angular orbital de un electrón sumado a su propio espín. Es lo que el acoplamiento espín-órbita pone en el hamiltoniano, y el II.8 lo encenderá. Aquí se resuelve el 2p, o sea =1\ell=1, y se resuelve del mismo modo que antes: bajando.

La serie da 112=32121\otimes\tfrac12 = \tfrac32 \oplus \tfrac12, con 4 + 2 = 6 estados, que es la dimensión del producto. Para j=32j = \tfrac32 el autovalor de J2\mathbf{J}^2 es 2j(j+1)=1542=3,752\hbar^2 j(j+1) = \tfrac{15}{4}\hbar^2 = 3{,}75\,\hbar^2. Arriba del todo sólo hay un estado con mj=m+ms=32m_j = m_\ell + m_s = \tfrac32, así que

32,32=1,+1 ⁣.\left|\tfrac32,\tfrac32\right\rangle = |1,+1\rangle\,|\!\uparrow\rangle .

Se aplica J=L+SJ_- = L_- + S_-. Por la izquierda, el coeficiente de escalera del artículo 01 con j=32j=\tfrac32 y m=32m=\tfrac32 da 1543212=3\hbar\sqrt{\tfrac{15}{4}-\tfrac32\cdot\tfrac12} = \hbar\sqrt3. Por la derecha, LL_- sobre 1,+1|1,+1\rangle da 21,0\hbar\sqrt2\,|1,0\rangle y SS_- sobre  ⁣|\!\uparrow\rangle da  ⁣\hbar|\!\downarrow\rangle. Igualando y dividiendo por 3\hbar\sqrt3:

32,+12=23  1,0 ⁣  +  13  1,+1 ⁣.\left|\tfrac32,+\tfrac12\right\rangle = \sqrt{\tfrac23}\;|1,0\rangle|\!\uparrow\rangle \;+\; \sqrt{\tfrac13}\;|1,+1\rangle|\!\downarrow\rangle .

Los dos coeficientes de Clebsch-Gordan son 2/3=\sqrt{2/3} = 0,81649658 y 1/3=\sqrt{1/3} = 0,57735027, y sus cuadrados suman 1 sin ayuda: el estado sale normalizado del propio procedimiento, no hay que normalizarlo a mano. El verificador los cruza contra la tabla de Clebsch-Gordan de sympy, calculada por la fórmula de Racah, y contra la diagonalización directa de J2\mathbf{J}^2 en la base producto: las tres vías coinciden con desviación medida 0. La misma bajada, un peldaño más, da

32,12=23  1,0 ⁣  +  13  1,1 ⁣,32,32=1,1 ⁣.\left|\tfrac32,-\tfrac12\right\rangle = \sqrt{\tfrac23}\;|1,0\rangle|\!\downarrow\rangle \;+\; \sqrt{\tfrac13}\;|1,-1\rangle|\!\uparrow\rangle , \qquad \left|\tfrac32,-\tfrac32\right\rangle = |1,-1\rangle|\!\downarrow\rangle .

Con eso está entero el multiplete de cuatro. El de dos sale de lo mismo que antes: en el plano mj=+12m_j = +\tfrac12 hay dos estados producto, uno se ha gastado en j=32j=\tfrac32, y el que queda —lo ortogonal— es 12,+12\left|\tfrac12,+\tfrac12\right\rangle. No hace falta ninguna tabla en ningún momento. Lo único que se usa es el coeficiente de escalera y la ortogonalidad, y las dos cosas están deducidas en el artículo 01.

Que el procedimiento produzca una tabla y no un montón de números sueltos también es comprobable, y está comprobado en cuatro series: la matriz de coeficientes es cuadrada, cumple UU=UU=1UU^\dagger = U^\dagger U = 1 y diagonaliza J2\mathbf{J}^2. Es decir: el cambio de base producto a base acoplada es un giro, no una aproximación, y por eso no se pierde ni se inventa ningún estado por el camino.

La comprobación que siempre tiene que salir

Cualquier cuenta de acoplamiento tiene una comprobación gratuita, y conviene hacerla siempre porque caza la mayoría de los errores de coeficiente:

Lz+Sz=mj,\langle L_z\rangle + \langle S_z\rangle = m_j\hbar ,

porque Jz=Lz+SzJ_z = L_z + S_z y el estado es autoestado de JzJ_z con autovalor mjm_j\hbar. Lo bonito es que los dos sumandos no son autovalores de nada: son valores esperados, y salen fraccionarios. En 32,12\left|\tfrac32,-\tfrac12\right\rangle, con pesos 23\tfrac23 sobre 1,0 ⁣|1,0\rangle|\!\downarrow\rangle y 13\tfrac13 sobre 1,1 ⁣|1,-1\rangle|\!\uparrow\rangle:

Lz=23(0)+13()=3=0,333,Sz=23 ⁣(2)+13 ⁣(+2)=6=0,167,\langle L_z\rangle = \tfrac23(0) + \tfrac13(-\hbar) = -\frac{\hbar}{3} = -0{,}333\,\hbar , \qquad \langle S_z\rangle = \tfrac23\!\left(-\tfrac\hbar2\right) + \tfrac13\!\left(+\tfrac\hbar2\right) = -\frac{\hbar}{6} = -0{,}167\,\hbar ,

y 1316=12-\tfrac13 - \tfrac16 = -\tfrac12, que es mjm_j. Sale. En 32,+12\left|\tfrac32,+\tfrac12\right\rangle los dos cambian de signo, +0,333+0{,}333\,\hbar y +0,167+0{,}167\,\hbar, y suman +12+\tfrac12\hbar. En el tope, 32,+32\left|\tfrac32,+\tfrac32\right\rangle, no hay mezcla y los dos valen sus autovalores enteros: +1,000+1{,}000\,\hbar y +0,500+0{,}500\,\hbar, que suman 32\tfrac32\hbar. El verificador hace esta comprobación en los seis estados de esta serie y nunca falla. Y por eso es útil: cuando falla, el error está en el coeficiente que se acaba de escribir.

Hay que resistir aquí una tentación concreta: leer Sz=/6\langle S_z\rangle = -\hbar/6 como «el espín apunta un poco hacia abajo». No apunta. Cada vez que se mida SzS_z de ese electrón saldrá +/2+\hbar/2 o /2-\hbar/2, nunca /6-\hbar/6; lo que pasa es que sale hacia abajo dos de cada tres veces. El /6-\hbar/6 es una media, y la media de un observable con dos autovalores no tiene por qué parecerse a ninguno de los dos.

El mecanismo que hace cambiar de signo a ⟨Sz⟩. Fíjese en de dónde salió el /6-\hbar/6: de restar dos pesos, 23\tfrac23 con el espín abajo y 13\tfrac13 con el espín arriba. Ahora piense en el otro multiplete de la misma serie, el j=12j=\tfrac12. Sus estados están construidos sobre los mismos dos productos —no hay otros con ese mjm_j— y son ortogonales a los de j=32j=\tfrac32. Ortogonal, con dos componentes reales, significa una sola cosa: los dos pesos se intercambian y uno de los coeficientes cambia de signo. Si en j=32j=\tfrac32 el peso grande está sobre el espín abajo, en j=12j=\tfrac12 el peso grande está sobre el espín arriba — y por tanto Sz\langle S_z\rangle sale con el signo contrario, con el mismo mjm_j, el mismo \ell, el mismo electrón y el mismo espín. La suma Lz+Sz\langle L_z\rangle + \langle S_z\rangle sigue valiendo mjm_j\hbar en los dos, porque Lz\langle L_z\rangle se mueve al revés para compensar. Ponerle números a esto es el problema 5 de la hoja, y la cuenta cabe en tres líneas con lo que hay aquí. Lo que no cabe es la intuición: nada en «un electrón con mj=+12m_j = +\tfrac12» dice hacia dónde apunta su espín, y la respuesta depende de un número — jj— que no es una propiedad del espín ni de la órbita, sino de cómo están acoplados.

El 2 de 2n², contado por j

El II.5 dejó el nivel nn del hidrógeno con n2n^2 orbitales, contados dos veces y sin espín. El I.4 llevaba años publicando la capacidad de las capas —2, 8, 18, 32— y el 2 de 2n22n^2 se daba como un hecho de la tabla periódica. Aquí se cierra el círculo: ese 2 es la degeneración del espín, y en el lenguaje de este módulo es la dimensión 2s+1=22s+1 = 2 de la representación j=12j = \tfrac12 del álgebra de momento angular. De ninguna otra cosa.

Y se puede contar de dos maneras que no se hablan entre sí. La primera es tensorial: cada nm|n\ell m\rangle se multiplica por los dos estados de espín, así que hay 2n22n^2 estados. La segunda usa la base acoplada: para cada \ell se forman los multipletes j=±12j = \ell\pm\tfrac12 y se suman sus 2j+12j+1 estados. Para n=2n=2: el =0\ell=0 da un solo j=12j=\tfrac12 con 2 estados, y el =1\ell=1 da j=32j=\tfrac32 con 4 y j=12j=\tfrac12 con 2. Total 2 + 4 + 2 = 8, que es 2n22n^2. El verificador hace las dos cuentas para n=1,2,3n = 1, 2, 3 y 4 y las compara con las cuatro filas de capacidad que el I.4 publica: 2 · 8 · 18 · 32, las cuatro coinciden.

La segunda cuenta parece un rodeo y no lo es: es la que sobrevive cuando la degeneración se rompe. En el momento en que el II.8 encienda LS\mathbf{L}\cdot\mathbf{S}, los ocho estados de n=2n=2 dejarán de tener la misma energía y se repartirán precisamente por multipletes: un 2S1/2^2S_{1/2}, un 2P1/2^2P_{1/2} y un 2P3/2^2P_{3/2}. La primera cuenta habrá dejado de significar algo; la segunda dirá exactamente cuántos estados hay en cada raya. Es el mismo total escrito en la base buena.

La raya de 21 cm: el contacto de Fermi, montado del todo

Ya está todo. El protón tiene espín ½ y momento magnético; el electrón también; y en el fundamental del hidrógeno los dos ocupan el mismo sitio. El acoplamiento de sus dos espines parte el nivel 1s en dos —el triplete arriba, el singlete abajo— y la transición entre ellos es la raya de 21 cm, la señal con la que se ha cartografiado el hidrógeno neutro de la Galaxia. El I.4 la nombra y dice que su tamaño es proporcional a ψ(0)2|\psi(0)|^2. Aquí se calcula.

El mecanismo es el término de contacto de Fermi. Un dipolo magnético puntual dentro de una densidad de carga produce, además del campo dipolar habitual, un término proporcional a δ3(r)\delta^3(\mathbf{r}) que sólo actúa donde el electrón está encima del núcleo. Su valor esperado es, por tanto, proporcional a la densidad de probabilidad del electrón en el origen, y el operador de espín que lleva es exactamente el del artículo:

ΔE=2μ03gegpμBμNψ(0)2  SeSpS=1SeSpS=02.\Delta E = \frac{2\mu_0}{3}\,|g_e|\,g_p\,\mu_B\,\mu_N\,|\psi(0)|^2\; \frac{\langle \mathbf{S}_e\cdot\mathbf{S}_p\rangle_{S=1} - \langle \mathbf{S}_e\cdot\mathbf{S}_p\rangle_{S=0}}{\hbar^2} .

El último factor es el 1 de la sección anterior —la diferencia +¼ − (−¾)— y no hay más espín en la fórmula. Los otros cuatro números vienen del artículo 03: ge=2,002319|g_e| = 2{,}002319, gp=5,586g_p = 5{,}586, μB=57,8838\mu_B = 57{,}8838 µeV/T y μN=31,5245\mu_N = 31{,}5245 neV/T. Y ψ(0)2|\psi(0)|^2 viene del I.4, que publica 2,15×10302{,}15\times10^{30} m⁻³ para el 1s; recalculado aquí como 1/πa031/\pi a_0^3 da 2,14806×10302{,}14806\times10^{30} m⁻³, que coincide con la cifra publicada dentro de su tercera cifra significativa (9,0×1049{,}0\times10^{-4} relativo, que es lo que cuesta redondear a tres cifras).

Sustituyendo sale ΔE=5,88426\Delta E = 5{,}88426 µeV, o sea

f=ΔEh=1422,808 MHz.f = \frac{\Delta E}{h} = 1422{,}808\ \text{MHz} .

La medida es 1420,405 751 7667 MHz (Essen, Donaldson, Bangham y Hope, Nature 229, 110, 1971; el año anterior Hellwig y colaboradores habían dado 1420,405 751 768 MHz, y las dos siguen citándose. Es una de las frecuencias mejor conocidas de la física, y por eso se escribe con todas sus cifras y no con tres). El cálculo se pasa un 0,17 %. Para una cuenta que no lleva ni un parámetro ajustable ni una constante de acoplamiento inventada, es un resultado espectacular. Y para lo que se sabe del hidrógeno, es un desacuerdo vergonzoso: es un 0,17 % sobre una raya que se conoce con doce cifras.

Y por dos caminos, porque un número calculado una vez no está verificado. La misma ΔE\Delta E se puede escribir sin μ0\mu_0, sin a0a_0 y sin μB\mu_B, como 43(ge/2)gp(me/mp)α2Eh\tfrac43\,(|g_e|/2)\,g_p\,(m_e/m_p)\,\alpha^2 E_h con EhE_h la energía de Hartree. Las dos formas no comparten ninguna entrada de la tabla de constantes salvo los dos factores gg, y coinciden en 1,2×10111{,}2\times10^{-11}. El 1422,808 no es un número que se haya sacado de un sitio: es un número que sale de dos sitios.

El factor perdido: a₀ entra al cubo

El 0,17 % no es misterioso, y encontrarlo es la lección más transferible de todo el módulo. El cálculo usó ψ(0)2=1/πa03|\psi(0)|^2 = 1/\pi a_0^3 con el radio de Bohr de siempre, a0=52,9177a_0 = 52{,}9177 pm. Pero a0a_0 se define con la masa desnuda del electrón, y el hidrógeno de verdad tiene un protón de masa finita que también se mueve. El II.5 ya lo pagó en las energías con la masa reducida; en el tamaño se paga igual, y hacia arriba:

aH=a0meμ=a0(1+memp)=52,9465 pm.a_H = a_0\,\frac{m_e}{\mu} = a_0\left(1 + \frac{m_e}{m_p}\right) = 52{,}9465\ \text{pm} .

Es un átomo mayor, y en una fracción exactamente igual a me/mpm_e/m_p, o sea el inverso de los mp/me=1836,15m_p/m_e = 1836{,}15 del artículo 03. Parece nada. Pero en ψ(0)2|\psi(0)|^2 el radio entra al cubo, y ahí está el factor perdido: la corrección no es (1+me/mp)1(1+m_e/m_p)^{-1} ni (1+me/mp)2(1+m_e/m_p)^{-2}, sino

(1+memp)3=11,632×103,\left(1 + \frac{m_e}{m_p}\right)^{-3} = 1 - 1{,}632\times10^{-3} ,

o sea −1,632 ‰. Aplicado al resultado anterior:

f=1420,485 MHz.f = 1420{,}485\ \text{MHz} .

Se pasa ahora 56,1 ppm: el desacuerdo cae de 1,69 × 10⁻³ a 5,61 × 10⁻⁵. La masa reducida se ha llevado el 97 % del error, y ha bastado con no equivocarse en un exponente. Ésa es la moraleja, y no es retórica: con el radio al cuadrado la corrección habría sido sólo dos tercios de ésta —el exponente es la corrección—, el resultado habría mejorado bastante menos y habría seguido pareciendo un éxito. Los errores de exponente no producen disparates: producen resultados ligeramente peores que se aceptan sin más.

Y los 56 ppm que quedan tampoco son ruido, ni redondeo, ni «el error del método». Son enormes al lado de la incertidumbre con la que se conoce la raya —doce cifras significativas— y tienen nombre: correcciones relativistas, correcciones de electrodinámica cuántica del mismo tipo que la anomalía aea_e del artículo 03, y —la más difícil— la estructura interna del protón, que no es un dipolo puntual sino una distribución de corriente de tamaño finito. Decir «coincide con la medida» sobre 56 ppm sería falso y además perezoso: lo que coincide es lo que el modelo contiene, y lo que sobra es la lista de lo que no contiene.

De la frecuencia medida —no de la calculada, que daría 21,1049 cm— salen las otras tres cifras con las que se maneja la raya en radioastronomía. La longitud de onda es λ=c/f=\lambda = c/f = 21,1061 cm —el nombre «21 cm» redondea y pierde la tercera cifra, lo cual está bien para un nombre y mal para una cuenta—. La energía del fotón es 5,874 µeV, siete órdenes de magnitud por debajo de las energías ópticas del mismo átomo. Y la tercera es la que explica todo lo demás.

Ejemplo resuelto 2 · Por qué una nube a 100 K emite una raya de 68 mK

Problema. El hidrógeno neutro del medio interestelar está típicamente entre 50 y 200 K. La raya de 21 cm se observa en emisión desde ese gas, y su brillo se usa para pesar el hidrógeno de una galaxia entera. (a) ¿Qué temperatura corresponde a la energía de la raya? (b) ¿Qué fracción de átomos está en el nivel de arriba a 100 K? (c) ¿Qué mide entonces el brillo de la raya: la temperatura del gas o la cantidad de gas?

Solución. (a) La escala térmica de una transición es hf/kBhf/k_B. Con f=1420,4057517667f = 1420{,}405\,751\,766\,7 MHz sale

hfkB=0,06817 K=68,2 mK.\frac{hf}{k_B} = 0{,}068\,17\ \text{K} = 68{,}2\ \text{mK} .

(b) El cociente de poblaciones de dos niveles a temperatura TT es (garriba/gabajo)ehf/kBT(g_{\text{arriba}}/g_{\text{abajo}})\,e^{-hf/k_BT}. Las degeneraciones las da este artículo: el nivel de arriba es el triplete, con tres estados, y el de abajo es el singlete, con uno. A 100 K el exponente vale 68,2 mK/100 K=6,8×10468{,}2\ \text{mK}/100\ \text{K} = 6{,}8\times10^{-4}, así que la exponencial es 1 con cuatro cifras y el cociente es simplemente 3:1. Tres cuartas partes de los átomos están arriba, y lo estarían igual a 50 K o a 200 K.

(c) Como el reparto no depende de TT en todo el rango del medio interestelar, el brillo de la raya no mide la temperatura: mide cuántos átomos hay en la línea de visión. Es un contador de hidrógeno, no un termómetro.

Resultado. Los 68,2 mK son el número que hace útil a la raya de 21 cm, y lo son por ser pequeños. Cualquier transición con hf/kBhf/k_B comparable a la temperatura del gas tendría el nivel de arriba casi vacío en el gas frío y sólo se vería en absorción, o se vería en emisión con un brillo que dependería exponencialmente de TT y que habría que corregir antes de contar nada. Aquí no: el gas más frío de la Galaxia está tres órdenes de magnitud por encima de la escala de la transición, el triplete está poblado a tope por pura contabilidad de estados, y la señal es proporcional a la densidad de columna y a nada más. Que la raya sea débil —es una transición prohibida a orden dipolar eléctrico, y por átomo es rarísima— deja de importar en cuanto hay bastante hidrógeno delante, y en una galaxia lo hay.

Hasta dónde vale el «campo débil»: 7,17 mT

Falta una honestidad. Todo lo anterior supone que los cuatro estados hiperfinos son el triplete y el singlete, o sea que el acoplamiento de los dos espines domina sobre cualquier otra cosa. En cuanto hay un campo magnético externo eso deja de ser cierto, porque el campo acopla cada espín por separado —y con fuerzas muy distintas, ya que el electrón es 658 veces más magnético que el protón— y por tanto empuja al sistema de vuelta hacia la base producto. El problema entero se resuelve diagonalizando una matriz 4 × 4, y para la transición central hay forma cerrada, la fórmula de Breit-Rabi:

f(B)=Ah1+x2,x(geμB+gpμN)BA,f(B) = \frac{A}{h}\sqrt{1+x^2} , \qquad x \equiv \frac{(|g_e|\mu_B + g_p\mu_N)\,B}{A} ,

con AA la separación hiperfina a campo cero. El verificador comprueba la forma cerrada contra la diagonalización numérica en siete campos, de 1 mT a 200 mT, con desviación 00 en todos.

La pregunta interesante es cuándo deja de valer el cálculo de arriba, y la respuesta tiene una trampa incorporada. Lo que se mira, y lo que dibujan casi todos los paneles del tema, es el cruce: el campo al que la energía Zeeman del electrón iguala a la hiperfina, o sea x=1x = 1. Ese campo vale 50,6 mT —equivalentemente, geμBB/h=1420|g_e|\mu_B B/h = 1420 MHz da 50,68 mT—. Parece razonable decir «por debajo de 50 mT el campo es débil». Es falso.

Porque lo que hay que comparar no es xx con 1, sino el corrimiento con la precisión que se pida. Con 1+x21+x2/2\sqrt{1+x^2} \approx 1 + x^2/2, un corrimiento del 1 % llega ya en x0,14x \approx 0{,}14, y resolviendo la ecuación exacta con brentq\texttt{brentq} sobre la diagonalización 4 × 4 —y cruzándolo con la forma cerrada, que coinciden a 7×10147\times10^{-14}— sale

B(1%)=7,17 mT.B(1\,\%) = 7{,}17\ \text{mT} .

Siete veces antes del cruce: el cociente medido es 7,05. El campo débil se rompe mucho antes de donde uno miraría, y ésa es la forma general del error, no una anécdota del hidrógeno: la escala en la que un término empieza a notarse no es la escala en la que empieza a dominar, y entre las dos hay tanto espacio como raíz de la precisión que se exija. La tabla lo hace visible:

Bf de la transición centralQué es ese campo
01420,405 751 7667 MHzla raya medida
1 mT1420,68 MHzel campo de un imán de nevera
5 mT1427,32 MHz
7,16 mT1434,55 MHzaquí ya se ha corrido un 1 % (7,17 mT)
10 mT1447,87 MHz
50,6 mT2008,62 MHzel cruce, x = 1
200 mT5790,42 MHzun electroimán pequeño: ×4,08

Y lo que esto abre. En el cruce de 50,6 mT la raya ya está en 2008,62 MHz: no es que el campo débil esté «cerca del límite», es que el límite ya está muy atrás. A 200 mT —un electroimán pequeño, dos órdenes de magnitud por debajo del imán de resonancia magnética del problema 2 de la hoja— la transición central está en 5790,42 MHz, cuatro veces más arriba, y hablar de «la raya de 21 cm» dentro de ese imán no significa nada. El II.8 recoge exactamente esto: el Zeeman de campo débil y el de campo fuerte no son dos aproximaciones a elegir por gusto, son dos regímenes con una frontera calculable, y calcular la frontera —en vez de suponer que está donde se cruzan las dos energías— es la mitad del trabajo. La otra mitad es aceptar que entre 7 y 50 mT no vale ninguna de las dos formas simples.

Lo que este módulo deja hecho

El módulo entró prometiendo lo que el II.3 había dejado escrito: que (+1)\ell(\ell+1) saldría de un conmutador igual que nn había salido de una cadena. Sale, y de paso el álgebra produce algo que la esfera del II.5 no admitía —los semienteros—, que es todo el espín. Con eso, el artículo 02 construyó la representación de dimensión 2 y su rotación a mitad de ángulo, el 03 le puso el factor gg y sus escalas magnéticas, y este artículo ha cerrado la operación que faltaba: sumar dos de ellos.

Lo que queda abierto, dicho con su nombre. La estructura fina no está calculada: aquí se ha montado LS\mathbf{L}\cdot\mathbf{S} como problema de acoplamiento, pero no se ha escrito el hamiltoniano que lo produce ni se ha calculado el desdoblamiento; eso es el II.8 y es teoría de perturbaciones. El signo del corrimiento Zeeman se ha usado sin discutir qué estado sube y cuál baja en cada régimen. Y los 56 ppm de la raya de 21 cm no se cierran en este sitio: hacen falta la ecuación de Dirac —el módulo III.7— y datos de estructura del protón que no son de mecánica cuántica no relativista. Es la misma clase de deuda que el II.5 dejó con los 0,147 meV de la ionización del hidrógeno: el modelo llega hasta donde llega, y se dice hasta dónde.

Ejercicios

Ejercicio 1 · Bajar dos veces, y comprobar que la escalera se cierra

Partiendo de 1,+1= ⁣|1,+1\rangle = |\!\uparrow\uparrow\rangle, aplica S=S1+S2S_- = S_{1-}+S_{2-} dos veces. (a) Comprueba en cada paso que la norma del vector obtenido coincide con el coeficiente S(S+1)M(M1)\hbar\sqrt{S(S+1)-M(M-1)} del artículo 01. (b) Comprueba que S1,1=0S_-|1,-1\rangle = 0, o sea que la escalera se cierra por abajo. (c) Comprueba sobre 1,0|1,0\rangle que S1z+S2z=M\langle S_{1z}\rangle + \langle S_{2z}\rangle = M\hbar, y di por qué esa comprobación no distingue 1,0|1,0\rangle de 0,0|0,0\rangle.

Solución

(a) Primer paso: S ⁣=( ⁣+ ⁣)S_-|\!\uparrow\uparrow\rangle = \hbar(|\!\downarrow\uparrow\rangle+|\!\uparrow\downarrow\rangle), cuya norma es 2=1,4142136\hbar\sqrt2 = 1{,}4142136\,\hbar. El coeficiente pedido, con S=1S=1 y M=+1M=+1, es 210=2\hbar\sqrt{2-1\cdot0} = \hbar\sqrt2. Coinciden, y dividiendo por esa norma sale 1,0|1,0\rangle con sus dos coeficientes +0,70710678. Segundo paso: sobre 1,0|1,0\rangle, SS_- da 2 ⁣\hbar\sqrt2\,|\!\downarrow\downarrow\rangle, y el coeficiente con M=0M=0 es 20(1)=2\hbar\sqrt{2-0\cdot(-1)} = \hbar\sqrt2. Otra vez.

(b) S ⁣=0S_-|\!\downarrow\downarrow\rangle = 0 porque S ⁣=0S_-|\!\downarrow\rangle = 0 en cada factor, y el coeficiente lo predice: con M=1M=-1 vale 2(1)(2)=0\hbar\sqrt{2-(-1)(-2)} = 0. Esto no es un detalle decorativo. Es el mismo mecanismo que en el artículo 01 fuerza a que jj sea entero o semientero: la escalera tiene que cerrarse por los dos extremos, y sólo se cierra si el número de peldaños es entero. En el oscilador del II.3 no pasaba —allí aa^\dagger no anula nunca a ningún estado— y por eso la torre no tiene tope.

(c) Sobre 1,0|1,0\rangle, cada término contribuye +22=0+\tfrac\hbar2 - \tfrac\hbar2 = 0 con peso ½, así que S1z=S2z=0\langle S_{1z}\rangle = \langle S_{2z}\rangle = 0 y la suma es 0=M0 = M\hbar. Sale. Y no distingue nada: el singlete tiene exactamente las mismas dos componentes con los mismos módulos, y también da 0 + 0. La segunda lección del apartado, que es la misma que la de CzzC_{zz}: una comprobación que pasa siempre es útil para cazar errores y completamente inútil para distinguir estados. Confundir las dos cosas es cómo se acaba creyendo que un estado está identificado porque ha superado un test.

Ejercicio 2 · Diseñar la medida que sí distingue

Un laboratorio prepara pares de espines y quiere saber si su fuente produce 1,0|1,0\rangle o 0,0|0,0\rangle. Su aparato mide S1zS_{1z} y S2zS_{2z} a la vez, muchas veces, y calcula la correlación. (a) ¿Qué obtendrá con cada una de las dos fuentes? (b) Demuéstralo escribiendo CzzC_{zz} sobre los dos estados, sin citar la tabla del artículo. (c) Propón la modificación mínima del aparato que sí resuelve la duda, y di qué valores tiene que ver en cada caso. (d) ¿Cambia algo si en vez de la correlación mide sólo S1z\langle S_{1z}\rangle?

Solución

(a) Lo mismo con las dos: −0,250 ħ². El experimento no puede distinguirlas, y ninguna cantidad de estadística lo arregla, porque no es un problema de ruido.

(b) S1zS2zS_{1z}S_{2z} es diagonal en la base producto, con autovalor +2/4+\hbar^2/4 sobre  ⁣|\!\uparrow\uparrow\rangle y  ⁣|\!\downarrow\downarrow\rangle y 2/4-\hbar^2/4 sobre los otros dos. Los dos estados en cuestión sólo tienen componentes sobre esos otros dos, cada una con peso ½, así que Czz=12(24)+12(24)=24C_{zz} = \tfrac12(-\tfrac{\hbar^2}{4}) + \tfrac12(-\tfrac{\hbar^2}{4}) = -\tfrac{\hbar^2}{4} en ambos. Obsérvese qué es lo que no ha intervenido: el signo relativo entre las dos componentes. Un operador diagonal no lo puede ver, y el signo relativo es lo único que separa los dos estados.

(c) Basta girar los dos analizadores 90° y medir CxxC_{xx} con el mismo montaje. Sale +0,250 ħ² con el triplete y −0,250 ħ² con el singlete: dos números que ni siquiera tienen el mismo signo. Con CyyC_{yy} pasa lo mismo. No hace falta un aparato nuevo, hace falta otro eje.

(d) No, y esto es lo más instructivo del ejercicio. S1z=0\langle S_{1z}\rangle = 0 en los dos estados —cada uno tiene la mitad de su peso con el espín 1 arriba y la mitad con él abajo—, así que tampoco los distingue. Y no es que haya que buscar otro observable de una partícula: ninguno sirve. Si se traza sobre el segundo espín, tanto 1,0|1,0\rangle como 0,0|0,0\rangle dejan el mismo estado reducido —la mezcla a partes iguales 121\tfrac12\,\mathbf{1}, que no tiene ninguna dirección preferida—, y dos estados con el mismo estado reducido son indistinguibles con cualquier medida local. Toda la información que los separa es conjunta, y por eso hay que medir los dos a la vez. Ése es exactamente el ingrediente que el I.5 artículo 02 llama entrelazamiento, aquí visto desde el lado del álgebra: el singlete no es «un espín arriba y otro abajo», es un estado que no tiene partes.

Ejercicio 3 · Qué orbitales pueden tener hiperfino de contacto

El término de contacto es proporcional a ψ(0)2|\psi(0)|^2. (a) Cerca del origen la parte radial de un orbital se comporta como Rn(r)rR_{n\ell}(r)\sim r^{\ell}; deduce de ahí qué orbitales tienen ψ(0)2=0|\psi(0)|^2 = 0. (b) De los ocho estados de n=2n=2 del hidrógeno, ¿cuántos sienten el término de contacto? (c) El desdoblamiento del 2s por contacto, ¿será mayor o menor que el del 1s, y en qué proporción aproximada? (d) ¿Quiere eso decir que el 2p no tiene estructura hiperfina?

Solución

(a) Con RrR\sim r^\ell, sólo =0\ell=0 sobrevive en r=0r=0: para 1\ell\ge1 la función radial se anula allí y con ella ψ(0)2|\psi(0)|^2. Sólo los orbitales s tienen densidad en el núcleo. Es el mismo factor rr^\ell que produce la barrera centrífuga del II.5, visto en su consecuencia menos obvia.

(b) Los ocho estados de n=2n=2 son 2n2=82\cdot n^2 = 8, repartidos en 2 del 2s y 6 del 2p. Sólo los dos del 2s sienten el contacto. Contado por multipletes, como en el artículo: el 2S1/2^2S_{1/2} (2 estados) sí, y el 2P1/2^2P_{1/2} (2) y el 2P3/2^2P_{3/2} (4) no.

(c) Mucho menor. ψn00(0)2|\psi_{n00}(0)|^2 escala como 1/n31/n^3, así que el 2s tiene ocho veces menos densidad en el origen que el 1s y su desdoblamiento por contacto es unas ocho veces menor. Como en el 1s vale 1420 MHz, en el 2s cabe esperar una octava parte de los 1420 MHz. La comprobación de orden de magnitud que conviene hacerse es que la raya sigue estando en microondas y no se va a otra región del espectro.

(d) No, y confundir «sin contacto» con «sin hiperfino» es el error que este apartado existe para romper. El contacto es uno de los términos del acoplamiento entre el espín nuclear y el electrón; los otros son la interacción dipolo-dipolo entre los dos momentos magnéticos —que no necesita que el electrón esté en el origen, sólo que esté cerca— y el acoplamiento del espín nuclear con la corriente orbital del electrón, que existe precisamente cuando 0\ell\neq0. En un orbital s esos dos se anulan por simetría esférica y queda sólo el contacto; en un p pasa lo contrario. La segunda lección: que un término se anule no significa que el efecto se anule, significa que hay que ir a buscar el siguiente, y el siguiente casi nunca es cero.

Ejercicio 4 · Dieciocho estados, contados por j

El nivel n=3n=3 del hidrógeno. (a) Cuenta sus estados como 2n22n^2. (b) Cuéntalos otra vez formando, para cada \ell, los multipletes j=±12j = \ell\pm\tfrac12 y sumando sus 2j+12j+1 estados; escribe la lista de términos espectroscópicos. (c) ¿Por qué el =0\ell=0 aporta un solo multiplete y no dos? (d) Cuando el II.8 rompa la degeneración, ¿cuál de las dos cuentas seguirá diciendo algo?

Solución

(a) 2n2=29=2n^2 = 2\cdot9 = 18, que es la tercera fila de la tabla de capacidades que publica el I.4: 2, 8, 18, 32.

(b) Con =0,1,2\ell = 0, 1, 2:

Multipletes2j+1Subtotal
0 (s)²S1/222
1 (p)²P1/2 · ²P3/22 · 46
2 (d)²D3/2 · ²D5/24 · 610

2 + 6 + 10 = 18. Coincide, y no era obvio que coincidiera: la primera cuenta multiplica por 2 y la segunda reparte 2(2+1)2(2\ell+1) en dos trozos de tamaños 22\ell y 2+22\ell+2.

(c) Porque j=12j = \ell - \tfrac12 daría j=12j = -\tfrac12, que no existe. La serie de Clebsch-Gordan empieza en j1j2|j_1-j_2|, y con j1=0j_1 = 0 el único valor posible es j2j_2. Dicho de otro modo: sin momento angular orbital no hay nada con lo que acoplar el espín, y el multiplete es el del espín solo.

(d) La segunda. En cuanto el hamiltoniano contenga LS\mathbf{L}\cdot\mathbf{S}, la energía dependerá de jj y los 18 estados se repartirán en cinco grupos con energías distintas, uno por fila de la tabla y dos por las filas con dos multipletes. La cuenta 2n22n^2 seguirá dando 18 y no dirá dónde están. La lección para llevarse, que es la del artículo entero: elegir la base no cambia cuántos estados hay, pero decide si los puedes nombrar de uno en uno cuando la degeneración se rompa. Y hay una comprobación gratis que conviene hacer siempre: los subtotales de la segunda cuenta tienen que sumar 2n22n^2 exactamente; si no suman, hay un multiplete de más o de menos.

Resumen en frío

Módulo II.6 · Momento angular y espín, en una tabla

Todo lo que este módulo deja utilizable, con dónde se dedujo cada cosa. Las filas se han rehecho desde el verificador, no copiadas del párrafo que las enuncia: scripts/verificar-espin.py remonta cada cifra por un segundo camino, falla si alguna no cuadra y sale con código 0. Convenios: la coma es decimal y el separador de millares es el espacio fino, sólo desde cinco cifras (1420 MHz va sin él y 28 025 MHz/T con él); \hbar se omite donde la unidad de la columna ya lo dice; las cifras con más de tres dígitos son las que la fuente justifica —redondear la raya de 21 cm a tres cifras borra justo la comparación que este módulo hace—. Pensada para leerse dentro de seis meses sin releer nada.

QuéFórmula o valorDónde
El álgebra del momento angular[Li, Lj] = iħ εijk Lk, deducido de x y p sin matricesart. 01
Lo que conmuta con todo[L², Lz] = 0, y L² no es el operador nuloart. 01
Operadores escaleraL± = Lx ± iLy, con [Lz, L±] = ±ħL±art. 01
… y su coeficienteL±|ℓ,m⟩ = ħ√(ℓ(ℓ+1) − m(m±1))|ℓ,m±1⟩, comprobado sobre 12 armónicos con sympyart. 01
La identidad que lo cierraL² = LL± + Lz² ± ħLz; y [L₊,L₋] = 2ħLzart. 01
Por qué ℓ(ℓ+1) y no ℓ²Lx²+Ly² ≥ 0 da |m| ≤ √(ℓ(ℓ+1)) = 1,4142 para ℓ = 1, mayor que ℓart. 01
… y qué añade el enteroLa cota sola no da |m| ≤ ℓ; lo da que la escalera se cierre en un número entero de pasosart. 01
El semientero, en el álgebraj = ½, 3/2, 5/2, 7/2 cumplen [Ji,Jj] y J² = j(j+1) con desviación ≤ 1,6 × 10⁻¹⁵art. 01
… y lo que la esfera diceℓ = ½ hace estallar Legendre: pendiente 0,3183 = |sen ℓπ|/π, el mismo 0,3183 del II.5art. 01
… y qué lo impone de verdadei2πm = −1 con m = ½: es la univaluación en el ángulo, no el álgebraart. 01
La promesa del II.3, cumplida[N,a] = −a es [Lz,L₋]/ħ = −L₋; la diferencia es que L₊|ℓ,ℓ⟩ = 0 y a† no anula nuncaart. 01
Matrices de PauliS = ħσ/2; σiσj = δij + iεijkσk; det σ = −1, tr σ = 0art. 02
… y de dónde salenSon 2Jk/ħ de la representación j = ½ de la escalera, no una definiciónart. 02
S² del espín ½¾ħ²·1, o sea s(s+1) con s = ½art. 02
Rotación de un espinore−iθ σ·n/2: serie y forma cerrada coinciden a 9 × 10⁻¹⁶ en 73 ángulosart. 02
La regla de Stern-GerlachP(+n|+z) = cos²(θ/2), no cos²θ; a 60° dan 0,75 y 0,25art. 02
El ángulo a la mitadU(2π) = −1 y U(4π) = +1 en j = ½; con j = 1 ya es +1 a 2πart. 02
El espinor gira la mitad que el vectorEl espinor gira θ/2 y el vector de Bloch θ, con la misma R: la cubierta es dobleart. 02
Tres imanesSGz → SGx → SGz da 1/8; sin el de en medio, 1/2art. 02
Cadena de N imanescos2N(90°/N): 25 % con 2, 78,1 % con 10, 96,0 % con 60art. 02
… y no es una analogíaSon las MISMAS cifras que los N polarizadores del I.5, porque el ángulo del espinor entra a la mitadart. 02
El haz de Stern-GerlachEfusivo: v³e−mv²/2kT, v más probable = √(3kBT/m) = 543 m/s para Ag a 1000 °Cart. 02
… y no es la del gas√(2kBT/m) daría 443 m/s, un 18 % menosart. 02
La deflexiónΔz = μz(∂B/∂z)L²/(6kBT): no depende de la masa; Ag y H se desvían igualart. 02
… y «la separación» son tres númerosmoda 0,0539 · vp 0,108 · media 0,162 mm, en razón ½ : 1 : 3/2art. 02
Lo que la placa de 1922 exigeCon 10 T/cm, μz = 0,928 μB (7,2 %, dentro del ±10 % declarado); μ ∝ 1/Gart. 02
Dos manchas y no tres2ℓ+1 es impar para todo ℓ entero: dos manchas no salen de ningún momento orbitalart. 02
g del electrónCODATA publica ge = −2,002 319 304 360 92: el signo va en g o en la fórmula, no en las dosart. 03
… y su consecuencia medibleμz(ms = +½) = −1,0012 μB: con ∂Bz/∂z > 0 el espín arriba se desvía hacia abajoart. 03
Magnetón de Bohreħ/2me = 57,8838 µeV/T ≙ 0,672 K/T ≙ 13,9962 GHz/T; derivado, no pegadoart. 03
Zeeman a 1 TΔE = |geBB = 115,90 µeV; a 300 K, kBT/μBB = 447, y a 1 K, 1,49art. 03
Frecuencia de Larmor|geB/h = 28 025 MHz/T (cinco cifras: a tres se pierde la comparación)art. 03
Frecuencia de ciclotróne/2πme = 27 992 MHz/T; fL/fc = g/2 exactamenteart. 03
… y su diferencia ES el experimentofL − fc = 32,5 MHz/T, o sea ae entero: es la trampa de Penningart. 03
Anomalía del electrónae = (|ge|−2)/2 = 1,159 65 × 10⁻³art. 03
… y el término de Schwingerα/2π = 1,161 41 × 10⁻³, que acierta ae con un 0,15 % (1516 ppm)art. 03
… en energíaaeμB = 67,1 neV/T: el «extra» de QED, medibleart. 03
g del protón5,586, que no es 2: el protón no es un Dirac puntual. Y gn = −3,826 sin cargaart. 03
Magnetón nucleareħ/2mp = 31,5245 neV/T = μB·me/mpart. 03
RMN del protóngpμN/h = 42,58 MHz/Tart. 03
Por qué el hiperfino es pequeñoep| = 658,2, y mp/me = 1836,15art. 03
Serie de Clebsch-Gordanj₁ ⊗ j₂ = |j₁−j₂| ⊕ … ⊕ j₁+j₂; (2j₁+1)(2j₂+1) = Σ(2J+1), medido en 9 seriesart. 04
Por qué hay que acoplarJ² conmuta con Jz y NO con J1z: el número cuántico bueno cambiaart. 04
½ ⊗ ½4 = 3 + 1: triplete S = 1 y singlete S = 0; S² tiene tres autovalores 2ħ² y uno 0art. 04
El triplete de M = 0|1,0⟩ = (|↑↓⟩+|↓↑⟩)/√2, los dos coeficientes +0,707107art. 04
El singlete|0,0⟩ = (|↑↓⟩−|↓↑⟩)/√2: el menos va aquí, por tres caminos independientesart. 04
… y qué significa el menosTriplete simétrico, singlete antisimétrico bajo intercambioart. 04
S₁·S₂(S² − S₁² − S₂²)/2, con S₁² = S₂² = ¾ħ²art. 04
⟨S₁·S₂⟩ del triplete+0,250 ħ², por los dos caminos, con desviación medida 0art. 04
⟨S₁·S₂⟩ del singlete−0,750 ħ², y NO −0,250: ése es Czz, un factor 3 menorart. 04
… y su diferencia1 ħ² exacto: es el factor entero que fija el tamaño de la raya de 21 cmart. 04
Correlaciones de |1,0⟩Czz = −0,250 · Cxx = Cyy = +0,250 ħ²; suman +0,250art. 04
Correlaciones de |0,0⟩Las tres −0,250 ħ²; suman −0,750art. 04
El eje que no distingueCzz vale −0,250 ħ² en los dos: hay que medir en x o en yart. 04
|↑↓⟩ en la base acoplada(|1,0⟩+|0,0⟩)/√2: autoestado de Sz con M = 0 y de S² de nada; ⟨S₁·S₂⟩ = −0,250 ħ²art. 04
1 ⊗ ½6 = 4 + 2: j = 3/2 con j(j+1) = 3,75 y j = ½art. 04
Los Clebsch del 2p, sin tablaBajar con J₋ desde |3/2,3/2⟩ da √⅔ = 0,81650 y √⅓ = 0,57735art. 04
La comprobación gratis⟨Lz⟩ + ⟨Sz⟩ = mjħ siempre, aunque ninguno de los dos sea autovalorart. 04
… en |3/2,−½⟩⟨Lz⟩ = −0,333 ħ y ⟨Sz⟩ = −0,167 ħ, que suman −½ħart. 04
… y en el tope |3/2,+3/2⟩+1,000 ħ y +0,500 ħ: sin mezcla, los dos son autovaloresart. 04
La tabla de Clebsch es un giroUU† = U†U = 1 a 2 × 10⁻¹⁶ y diagonaliza J²: no se pierde ni se inventa ningún estadoart. 04
El 2 de 2n²Es 2s+1 = 2, la dimensión de j = ½. Da 2 · 8 · 18 · 32, las cuatro filas del I.4art. 04
… contado por jΣℓ,j(2j+1) = 2n²: 2 · 8 · 18 · 32. Es la cuenta que sobrevive cuando el II.8 rompa la degeneraciónart. 04
|ψ₁₀₀(0)|²1/πa₀³ = 2,148 × 10³⁰ m⁻³, o 2,15 × 10³⁰ como lo publica el I.4art. 04
… y ℓ > 0|ψ(0)|² = 0 por el factor r: sólo los orbitales s tienen hiperfino de contactoart. 04
Contacto de FermiΔE = (2μ₀/3)|ge|gpμBμN|ψ(0)|² × 1, donde el 1 es la diferencia de ⟨S₁·S₂⟩art. 04
… con a₀ desnudoΔE = 5,884 µeV ⇒ 1422,808 MHz, un 0,17 % altoart. 04
… por un segundo camino(4/3)(|ge|/2)gp(me/mp)α²Eh, sin μ₀ ni a₀: coinciden a 1,2 × 10⁻¹¹art. 04
El factor perdidoaH = a₀(1+me/mp) = 52,9465 pm, y entra al CUBO: −1,632 ‰art. 04
… con la masa reducida1420,485 MHz, a 56,1 ppm de la medida: la μ se lleva el 97 % del errorart. 04
La raya, medida1420,405 751 7667 MHz (Essen et al., 1971)art. 04
… y lo que sobra56 ppm: relatividad, QED y estructura del protón. No es redondeoart. 04
λ de la rayac/f = 21,1061 cm: «21 cm» pierde la tercera cifraart. 04
Energía de la raya5,874 µeVart. 04
Escala térmicahf/kB = 68,2 mK ⇒ poblaciones 3:1 a cualquier T del medio interestelarart. 04
Breit-Rabif(B) = (A/h)√(1+x²), x = (|geB+gpμN)B/A; exacta contra la 4×4 en 7 camposart. 04
Dónde se rompe el campo débil7,17 mT para un 1 % de corrimientoart. 04
… frente al crucex = 1 en 50,6 mT: 7,05 veces más tarde de donde empieza a notarseart. 04
… y en un electroimán pequeño200 mT ⇒ 5790,42 MHz: cuatro veces la raya. Ahí no hay «21 cm»art. 04
ConstantesCODATA 2022 vía src/data/constants.ts: μB = 9,274 010 066 × 10⁻²⁴ J/T · μN = 5,050 784 × 10⁻²⁷ J/T · a₀ = 52,9177 pmtodo