Física Cuántica

UNIVERSIDAD DE MURCIA Departamento de Física Miguel Ortuño Ortín [email protected] Miguel Albaladejo Serrano [email protected]

Índice general

Índice

iii

Figuras

iv

Bibliografía

iv

Notas del autor

vi

1. Orígenes de la Física Cuántica

1

1.1. Radiación del cuerpo negro. Postulado de Planck . . . . . . . . . . . . .

1

1.2. Efecto fotoeléctrico. Efecto Compton . . . . . . . . . . . . . . . . . . . .

2

1.3. Postulado de De Broglie . . . . . . . . . . . . . . . . . . . . . . . . . . .

3

1.4. Modelo atómico de Bohr . . . . . . . . . . . . . . . . . . . . . . . . . . .

4

1.5. Experimento de la doble rendija . . . . . . . . . . . . . . . . . . . . . . .

5

1.5.1. Experimento con proyectiles (clásicos) . . . . . . . . . . . . . . .

6

1.5.2. Experimento con ondas . . . . . . . . . . . . . . . . . . . . . . . .

7

1.5.3. Experimentos con proyectiles (cuánticos): electrones . . . . . . . .

8

Átomo de Rutherford . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

10

A.1.

Potencia radiada . . . . . . . . . . . . . . . . . . . . . . . . . . .

10

A.2.

Conclusión . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

13

A.

2. Ecuación de Schrödinger

14

i

Índice general

Índice general

1.

Paquetes de ondas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

14

2.

Ecuación de Schrödinger . . . . . . . . . . . . . . . . . . . . . . . . . . .

15

3.

Interpretación probabilística . . . . . . . . . . . . . . . . . . . . . . . . .

17

4.

Operadores. Valores esperados . . . . . . . . . . . . . . . . . . . . . . . .

18

5.

Representación de momentos . . . . . . . . . . . . . . . . . . . . . . . . .

22

5.1.

Función de onda en p como transformada de Fourier . . . . . . . .

22

5.2.

Teorema de Parseval . . . . . . . . . . . . . . . . . . . . . . . . .

22

6.

Autovalores y autovectores . . . . . . . . . . . . . . . . . . . . . . . . . .

22

7.

Evolución temporal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

23

3. Problemas unidimensionales 1.

25

Introducción. Generalidades . . . . . . . . . . . . . . . . . . . . . . . . . 1.1.

25

Ecuación de Schrödinger para potenciales constantes. Continuidades y discontinuidades . . . . . . . . . . . . . . . . . . . . . . . . . .

25

Simetrías. Operador paridad . . . . . . . . . . . . . . . . . . . . .

27

Potencial escalón . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

28

2.1.

Reflexión parcial . . . . . . . . . . . . . . . . . . . . . . . . . . .

28

2.2.

Reflexión total . . . . . . . . . . . . . . . . . . . . . . . . . . . .

30

Pozo de potencial finito . . . . . . . . . . . . . . . . . . . . . . . . . . . .

31

3.1.

Estados ligados . . . . . . . . . . . . . . . . . . . . . . . . . . . .

31

3.2.

Coeficientes de transmisión . . . . . . . . . . . . . . . . . . . . . .

33

4.

Barreras de potencial . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

36

5.

Pozo de potencial infinito . . . . . . . . . . . . . . . . . . . . . . . . . . .

37

5.1.

Centrado en el origen . . . . . . . . . . . . . . . . . . . . . . . . .

37

5.2.

Pozo asimétrico . . . . . . . . . . . . . . . . . . . . . . . . . . . .

39

Potenciales delta . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

42

1.2. 2.

3.

6.

4. Estructura y postulados 1.

43

Preliminares. Funciones de onda . . . . . . . . . . . . . . . . . . . . . . . ii

44

Índice general

Índice general

2.

Espacio de estados. Notación de Dirac . . . . . . . . . . . . . . . . . . .

47

A.

Resumen de los postulados . . . . . . . . . . . . . . . . . . . . . . . . . .

49

A.1.

Descripción del estado de un sistema . . . . . . . . . . . . . . . .

49

A.2.

Descripción de magnitudes físicas . . . . . . . . . . . . . . . . . .

49

A.3.

Teoría de la medida . . . . . . . . . . . . . . . . . . . . . . . . . .

49

A.4.

Valores posibles de una medida . . . . . . . . . . . . . . . . . . .

50

5. Espín. Momentos Angulares

51

1.

Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

51

2.

Evidencia experimental . . . . . . . . . . . . . . . . . . . . . . . . . . . .

52

2.1.

Estructura fina de las líneas espectrales . . . . . . . . . . . . . . .

52

2.2.

Efecto Zeeman anómalo . . . . . . . . . . . . . . . . . . . . . . .

52

2.3.

Existencia de momentos angulares semienteros . . . . . . . . . . .

53

3.

Postulados de Pauli . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

53

4.

Espín 1/2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

54

5.

Descripción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

55

6.

Momento angular total . . . . . . . . . . . . . . . . . . . . . . . . . . . .

58

7.

Suma de espines . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

59

8.

Método general . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

62

8.1.

Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

62

8.2.

Cambio de base . . . . . . . . . . . . . . . . . . . . . . . . . . . .

63

8.3.

Autovalores de Jz . . . . . . . . . . . . . . . . . . . . . . . . . . .

64

8.4.

Autovalores de J2 . . . . . . . . . . . . . . . . . . . . . . . . . . .

64

8.5.

Autovectores de J2 y Jz . . . . . . . . . . . . . . . . . . . . . . .

66

iii

Índice de figuras

1.1. Ondas estacionarias del electrón . . . . . . . . . . . . . . . . . . . . . . .

5

1.2. Dispositivo para el experimento de la doble rendija . . . . . . . . . . . .

6

1.3. Experimento de la doble rendija con proyectiles macroscópicos . . . . . .

7

1.4. Experimento de la doble rendija para un fenómeno ondulatorio . . . . . .

8

1.5. Experimento de la doble rendija con los electrones. . . . . . . . . . . . .

9

1.6. Curvas con interferencias suavizadas . . . . . . . . . . . . . . . . . . . .

11

3.1. Potencial escalón . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

28

3.2. Representación del coeficiente de transmisión T en función de α = V0 /E.

30

3.3. Pozo cuadrado de potencial . . . . . . . . . . . . . . . . . . . . . . . . .

32

3.4. Representación de estados ligados en un pozo cuadrado . . . . . . . . . .

34

3.5. Barrera cuadrada de potencial . . . . . . . . . . . . . . . . . . . . . . . .

36

3.6. Pozo de potencial infinito centrado en el origen

. . . . . . . . . . . . . .

37

3.7. Primeros estados separados según su forma: cos kn x o sin kn x . . . . . . .

40

3.8. Los cuatro primeros estados . . . . . . . . . . . . . . . . . . . . . . . . .

41

5.1. Valores de M y autovectores asociados. . . . . . . . . . . . . . . . . . . .

65

5.2. Valores de M y degeneración asociada. . . . . . . . . . . . . . . . . . . .

65

iv

Bibliografía [1] C. Cohen-Tannoudji. Quantum Mechanics, volume I. [2] C. Cohen-Tannoudji. Quantum Mechanics, volume II. [3] R. Ferrer. Mecánica Cuántica I. Apuntes Universidad de Chile. [4] R. P. Feynman, R. B. Leighton, and M. Sands. The Feynman lectures on Physics, volume I. Pearson, 1963. [5] R. P. Feynman, R. B. Leighton, and M. Sands. The Feynman lectures on Physics, volume III. Pearson, 1963. [6] F. J. Ynduráin Muñoz. Mecánica Cuántica. Alianza Editorial, 2003.

v

Notas del autor Estas notas de clase están sacadas de la asignatura Física Cuántica, del 3ž curso de la Licenciatura en Física de la Universidad de Murcia, impartida por Miguel Ortuño Ortín, al que el autor está agradecidísimo. Concretamente, se basan en el desarrollo de la asignatura 2003-2004. Por supuesto, están sin terminar, están en continuo desarrollo, y agradeceré todas las sugerencias y comentarios que me mandéis a la dirección de correo de la portada. Sin más, os dejo con Feynman.

Dijo una vez un poeta: “El Universo entero está contenido en un vaso de vino”. Probablemente nunca sabremos lo que quería decir, ya que los poetas no escriben para ser comprendidos. Pero es cierto que si miramos un vaso de vino lo suficientemente cerca, vemos el Universo entero. Ahí estan las cosa de la física: el líquido que se arremolina y se evapora dependiendo del viento y del tiempo, las reflexiones en el vidrio, y nuestra imaginación agrega los átomos. El vidrio es un destilado de las rocas terrestres y en su composición vemos los secretos de la edad del universo y la evolución de las estrellas. ¿Qué extraño arreglo de elementos químicos hay en el vino? ¿Cómo llegaron a ser? Están los fermentos, las enzimas, los sustratos y los productos. Allí en el vino se encuentra la gran generalización: toda vida es fermentación. Nadie puede descubrir la química del vino sin descubrir, como lo hizo Louis Pasteur, la causa de muchas enfermedades. ¡Cuán vívido es el vino tinto que imprime su existencia dentro del conocimiento de quien lo observa! ¡Si nuestras pequeñas mentes, por alguna conveniencia, dividen este vaso de vino, este universo, en partes – física, biología, geología, astronomía, psicología, etc. . . , recuerden que la naturaleza no lo sabe! Así, reunamos todo de nuevo sin olvidar en última instancia para qué sirve. Dejemos que nos dé un último placer: ¡bébanselo y olvídense de todo! Richard P. Feynman, 1918-1988

vi

Capítulo 1

Orígenes de la Física Cuántica Vamos a comenzar describiendo los experimentos que hicieron patente la necesidad de una nueva teoría física, que hoy en día conocemos como Física Cuántica, para explicar los fenómenos de la naturaleza.

1.1.

Radiación del cuerpo negro. Postulado de Planck

Un cuerpo negro es, por definición, una cavidad que absorbe completamente toda la luz que llega a ella. Se trata de un concepto ideal, pero una buena realización del mismo es una cavidad con las paredes de negro y un pequeño orificio al exterior. Imaginemos dicho cuerpo negro con las paredes a una temperatura absoluta T . Podemos medir la intensidad de la radiación que sale de él por el pequeño orificio en función de la temperatura T y la frecuencia ν. En muy buena aproximación, dicha intensidad será igual a la que incide en las paredes del cuerpo negro. A través de dos conceptos clásicos muy generales podemos predecir cómo debería de ser la intensidad en el marco de la teoría clásica. En primer lugar, el teorema de equipartición nos dice que la energía media por grado de libertad ha de ser igual a kB T , siendo kB la constante de Boltzmann. En segundo lugar, el número de grados de libertad del campo electromagnético en el interior de un cuerpo negro viene dado por el número de modos estacionarios. Teniendo en cuenta estos dos ingredientes, Rayleigh-Jeans dedujeron la ley que lleva su nombre y que nos dice que la intensidad de la radiación en función de la frecuencia es 8πν 2 (1.1) u(ν, T ) = 3 kB T c Esta teoría coincidía con los experimentos para frecuencias bajas (T /νgj1 ,j2 10−10 K s), pero fallaba para frecuencias altas. Además, la teoría tenía que ser incorrecta ya que predecía una energía total infinita, pues el número de modos diverge en el límite de altas frecuencias, mientras que la energía por modo se mantiene constante. Este resultado se conoce como divergencia o catástrofe ultravioleta. 1

1.2. Efecto fotoeléctrico. Efecto Compton Capítulo 1. Orígenes de la Física Cuántica Planck (1900,1901) propuso una fórmula de la intensidad de la radiación del cuerpo negro que se ajustaba perfectamente a los datos experimentales. Estaba basada en la proposición de que el intercambio energético entre las paredes del cuerpo negro y el campo electromagnético sólo podía realizarse en cuantos de energía, no de forma continua. Planck no llegó a proponer que la radiación electromagnética en sí estuviera cuantizada. El cuanto de energía o unidad energética elemental sería proporcional a la frecuencia ν y vendría dado por ∆E = hν (1.2) en donde h es una constante universal conocida como constante de Planck, cuyo valor es h = 6,626 · 10−34 Js. Cada modo propio de oscilación del campo electromagnético puede vibrar con una energía múltiplo entero del cuanto de energía. A una temperatura T , la probabilidad de que un modo tenga una energía E es proporcional al factor de Boltzmann exp{−E/kB T } de modo que su energía media será: P∞ −En /kb T n=0 En e , En = nhν E(ν) = N (ν) P ∞ −En /kb T n=0 e Sumando la serie llegamos a: E(ν) = N (ν)

hν ehν/KB T − 1

Esta expresión está de acuerdo con la teoría clásica, ecuación (1.1), en el límite de ν pequeño y, como hemos dicho, reproduce perfectamente los resultados experimentales.

1.2.

Efecto fotoeléctrico. Efecto Compton

Fue Einstein (1905) quien realmente propuso la cuantización de las ondas electromagnéticas para explicar el efecto fotoeléctrico. Si se ilumina un metal (hacemos incidir radiación electromagnética sobre él), podemos observar que se “arrancan” electrones, es decir, son emitidos. Podemos variar la frecuencia con la que la luz–radiación incide y medir la energía y la dependencia de la energía con la frecuencia. Experimentalmente se observan los siguientes hechos: Existe una frecuencia umbral, más hacia abajo de la cual no se emiten electrones El número de electrones que se emiten, una vez superada dicha frecuencia umbral, es proporcional a la intensidad de la radiación incidente, pero la energía de los electrones emitidos no. 2

1.3. Postulado de De Broglie

Capítulo 1. Orígenes de la Física Cuántica

Estos dos hechos serían lógicos si la luz estuviese compuesta de entidades individuales, partículas (llamadas fotones) de energía E = hν cada una. Esta fue, precisamente, la propuesta de Einstein. Esta propuesta era revolucionaria1 . Las interacciones se producirían así de manera individual. Cuánticamente, se tiene I = N hν. Una intensidad alta puede deberse a que N o ν sean altos, pero clásicamente no se sabía distinguir esto. Tenemos entonces los siguientes resultados: Si lanzamos fotones de muy baja energía, da igual cuantos lancemos (es decir, el tiempo que transcurra desde que empecemos a lanzarlos), ninguno de ellos conseguirá que sean arrancados electrones. Para una intensidad dada, la energía del electrón emitido depende proporcionalmente de la frecuencia. 1 2 mv = hν − W, W = hνc donde νc es la frecuencia umbral 2 max Otra de las principales evidencias del carácter corpuscular de la luz fue el efecto Compton. Al “golpear” con fotones la materia (generalmente electrones) se obtienen unos resultados que no son explicables clásicamente. Si lo son, sin embargo, al considerar que la interacción fotón–electrón se produce como un choque elástico (este punto es el que pone de manifiesto que la luz está compuesta por partículas). Se tiene, para la energía del fotón, según la teoría de la relatividad restringida: q hν h = E = m20 c4 + p2 c2 = pc = hν ⇒ p = c λ Compton usó lo que Einstein ya había dicho para explicar su experimento. A las siguientes relaciones se las llama relaciones de Planck–Einstein: E = hν (1.3a) h p = (1.3b) λ En la física clásica se distinguía perfectamente entre las partículas y las ondas. Con estos experimentos mencionados se puso de manifiesto que, de algún modo, las ondas electromagnéticas estaban compuestas por partículas. Es decir, las ondas electromagnéticas (la luz) tienen ambos caracteres: ondulatorio y corpuscular.

1.3.

Postulado de De Broglie

Puesto que el campo electromagnético, que desde el siglo XIX (a partir de las leyes de Maxwell) se había considerado como un fenómeno ondulatorio, mostraba ahora 1

Desde los experimentos de Young, que estudiamos más adelante, en torno a 1800, nadie dudaba acerca de la naturaleza ondulatoria de la luz

3

1.4. Modelo atómico de Bohr

Capítulo 1. Orígenes de la Física Cuántica

también un carácter corpuscular, de Broglie postuló (1923) que las partículas materiales pueden tener también un comportamiento ondulatorio. Propuso que toda partícula tiene una onda asociada cuya longitud de onda es λ = h/p, en analogía con la fórmula de Compton–relaciones de Planck–Einstein. Por ejemplo, un electrón, bajo determinadas circunstancias podría tener comportamientos ondulatorios tales como la difracción, que fueron observados poco después de haberlos predicho de Broglie.

1.4.

Modelo atómico de Bohr

Los experimentos de Ruhterford, en 1911, demostraron que la materia debe estar compuesta por átomos, y también permitieron discriminar la naturaleza de estos2 . Bombardeando con partículas α unas láminas de oro muy delgadas, observó que algunas se desviaban, debido a la existencia de fuertes campos eléctricos, mientras que la gran mayoría de las partículas α atravesaban las láminas de oro. Esto le permitió deducir la siguiente estructura en el átomo: Núcleo El núcleo central debería estar compuesto por una carga positiva, y en el debería concentrarse la mayor parte de la masa. Corteza electrónica Los electrones, de carga negativa (en total, el átomo ha de ser neutro) que se encontrarían orbitando alrededor del núcleo. Una amplia zona de vacío entre ambos. Esta imagen estaría en principio de acuerdo con la imagen clásica (newtoniana) de la mecánica: los electrones orbitan alrededor del núcleo como los planetas alrededor del Sol. El desacuerdo surgía con las leyes de Maxwell: según éstas, una carga en movimiento acelerado (y el movimiento circular lo es) debería radiar energía de forma continua, y acabar precipitando en el núcleo, es decir, el átomo debería ser inestable, cosa que es a todas luces falso3 . Además, en los espectros de absorción–emisión se observaban líneas discretas, es decir, que las energías eran “tomadas” o emitidas en forma de cuantos, y no de manera continua. Esto condujo a Bohr a postular la existencia de órbitas estables para los electrones en los átomos, sin saber muy bien a que podría deberse la existencia de dichas órbitas. Al pasar un electrón de una órbita a otra lo hace absorbiendo o emitiendo un cuanto de energía, un fotón, lo que permitiría explicar las líneas espectrales discretas observadas. Las 2

Él mismo afirma: “Un cálculo simple demuestra que el átomo debe de ser asiento de un intenso campo eléctrico para que se produzca una gran desviación en una colisión simple...”

3

Para un cálculo más detallado ver apéndice A, pág. 10.

4

1.5. Experimento de la doble rendija

Capítulo 1. Orígenes de la Física Cuántica

órbitas son estables cuando la onda asociada al electrón es estable, esto es, estacionaria: 2πr = nλ

(1.4)

Los múltiplos de la longitud de onda han de “caber” exactamente en la longitud de la circunferencia de la órbita, y esto implica la cuantización del momento angular: h h deBroglie h L = pr = r= = n = n~ ⇒ L = n~ (1.5) λ 2πr/n 2π Así pues, añadiendo este ingrediente a la fórmula clásica de un movimiento orbital, se tienen ciertas órbitas posibles, con sus correspondientes energías, y no cualquier órbita.

Figura 1.1: Ondas estacionarias del electrón Pero extender esto a otros sistemas resultaba muy complicado. Pronto quedó patente que eran necesarias ciertas renovaciones profundas a nivel fundamental en la física, y no sólo la adición de una serie de recetas como podía ser la ecuación (1.5). Se hizo necesario la formulación axiomática de una nueva Mecánica, la Mecánica Cuántica, que revolucionaría (y aún hoy se siente ese revuelo) los conceptos y las aplicaciones de la Física.

1.5.

Experimento de la doble rendija

En esta sección seguiremos básicamente el libro de Feynman [4]. Para tratar de entender el comportamiento cuántico de los electrones (usaremos los electrones para nuestros ejemplos), compararemos su comportamiento en un dispositivo experimental especial, con el comportamiento habitual que tendrían otros fenómenos físicos de carácter ondulatorio (las ondas en el mar, o cosas similares). Dicho experimento se encuentra esquematizado en la figura 1.2. 5

1.5. Experimento de la doble rendija

Capítulo 1. Orígenes de la Física Cuántica

Figura 1.2: Dispositivo para el experimento de la doble rendija

1.5.1.

Experimento con proyectiles (clásicos)

Para este experimento, nuestro dispositivo consiste en un cañón que dispara proyectiles. Entre el cañón que dispara y la pared que detecta los proyectiles hay otra pared con dos rendijas intermedias R1 y R2 , que podemos tener abiertas y cerradas, independientes entre sí. En la pared final, de contención, disponemos de un detector, algo que permite controlar cuántos proyectiles y donde exactamente golpean la pared, o, lo que es lo mismo, determinar la probabilidad de que un proyectil golpee la pared en un punto concreto, en función de posición en la pared, digamos, la coordenada x. Trataremos ahora con cuerpos macroscópicos, como pudieran ser pelotas de tenis o balas, que no puedan partirse por la mitad, ni hacer cosas extrañas. Observaremos entonces, en nuestro experimento, que los proyectiles llegan al detector individualmente, y pasan, como esperábamos, por una y solo una de las rendijas intermedias. Si hacemos el experimento con R1 abierta y R2 cerrada, encontraremos un diagrama de probabilidad como el que se recoge en la figura 1.3. Si lo hacemos al revés, con R1 abierta y R2 cerrada, el diagrama es equivalente pero desplazado. No nos extraña, tampoco, que al hacer el experimento con ambas rendijas abiertas, el diagrama que obtenemos se parece mucho a los anteriores, pero ahora está centrado, y en realidad, es la suma de ambos resultados. Es decir, se tiene P12 (x) = P1 (x) + P2 (x). Decimos entonces que en este experimento no hay interferencia.

6

1.5. Experimento de la doble rendija

Capítulo 1. Orígenes de la Física Cuántica

Figura 1.3: Experimento de la doble rendija con proyectiles macroscópicos

1.5.2.

Experimento con ondas

Ahora deseamos hacer esencialmente el mismo experimento, pero con ondas de agua, o algo parecido. Evidentemente, ya no tendremos un cañón, sino una fuente puntual de ondas, que se propaguen por igual en todas las direcciones, en particular hacia nuestro sistema. Ahora ya no tendremos un detector que nos sirva para contar el número de proyectiles que chocan con la pared, sino un detector de intensidad, o algo que nos permita medir la intensidad de la onda al llegar a la pared. De nuevo hacemos el experimento con una y después otra rendijas abiertas, pero cuando hacemos el experimento con ambas rendijas abiertas, el resultado que obtenemos no es la suma de los dos resultados anteriores, por lo que decimos que existe interferencia entre ambas ondas. £Qué es lo que ocurre? En algunos lugares, donde I12 tiene sus máximos, las ondas provinientes de una y otra rendija están en fase, y las ondas se suman para dar una gran amplitud. Por contra, hay algunos lugares, justo allí donde la amplitud de I12 es mínima, en los que las ondas van desfasadas con una diferencia justamente de π. La intensidad es proporcional al cuadrado de la amplitud de la onda, es decir: I(x) = |A(x)|2 = |A1 (x) + A2 (x)|2 = |A1 (x)2 | + |A2 (x)|2 + 2|A1 (x)||A2 (x)| = p =I1 (x) + I2 (x) + 2 I1 (x)I2 (x) 6= I1 (x) + I2 (x) En definitiva, en los experimentos con ondas, la intensidad puede tomar cualquier valor, no llega en cuantos, y se obtienen interferencias.

7

1.5. Experimento de la doble rendija

Capítulo 1. Orígenes de la Física Cuántica

Figura 1.4: Experimento de la doble rendija para un fenómeno ondulatorio

1.5.3.

Experimentos con proyectiles (cuánticos): electrones

Ahora imaginamos un experimento similar al de la sección 1.5.1. Construimos un cañón de electrones, de modo que conseguimos (imaginariamente)4 que podamos disparar los electrones de uno en uno. Al tener cada una de las rendijas abiertas por separado, al llegar los electrones a la pantalla, estos golpean por separado, golpean en un solo sitio: no son como una onda, o eso parece. Pero si abrimos ambas rendijas a la vez, entonces no obtenemos una gráfica para la probabilidad que sea la suma de ambos resultados por separado, sino que obtenemos una gráfica en la que se observa interferencia: P12 6= P1 + P2 . La situación se encuentra representada en la figura 1.5. Esto no es lo que esperaríamos si cada electrón pasase por una y sólo una de las rendijas R1 o R2 . Podríamos pensar entonces que los electrones no van por una y sólo una de las dos rendijas, que a lo mejor van en formas más complicadas: a lo mejor el electrón se divide en dos mitades... aunque esto es imposible, porque hemos visto que llegan de uno en uno a la pantalla... Podemos intentar explicar, y de hecho mucha gente lo ha intentado ya, la curva de probabilidad para el golpeo del electrón, y no tendremos éxito. Sin embargo, la matemática necesaria para explicarla es muy simple, pues es la misma que la de la figura 1.4, y aquella figura era simple, pues describía una onda. Concluimos entonces que los electrones llegan a la pantalla como 4

Este experimento, descrito como lo hemos hecho, es en realidad impensable, pues el dispositivo debería ser muy pequeño para apreciar los efectos que queremos ver. Sabemos lo que pasará porque se han hecho experimentos distintos, pero esencialmente idénticos, en los que suceden estos fenómenos que describimos

8

1.5. Experimento de la doble rendija

Capítulo 1. Orígenes de la Física Cuántica

Figura 1.5: Experimento de la doble rendija con los electrones. partículas, pero la probabilidad de llegar a un punto u otro de la pantalla se distribuye como si fuera una onda. Por esto decimos que los electrones a veces se comportan como una onda y otras veces como partículas, y esto es llamado por todos y con razón, la dualidad onda-partícula. Ahora podremos interesarnos un poco más acerca de por donde pasa el electrón en realidad: por la rendija R1 , por R2 , por ambas, o por ninguna. Podríamos poner a nuestro aparato un detector más, una fuente de luz, cuya intensidad podemos variar, que nos ayude a decir por donde ha pasado el electrón. Así, si ponemos dicha fuente de luz entre ambas rendijas, al pasar un electrón por la rendija R1 obtendremos un destello luminoso en sus alrededores, si pasa por R2 , un destello cerca de dicha rendija, si el electrón se divide, dos destellos, o algo así, lo importante es que esto nos ayudará a saber que hace realmente el electrón. Si llevamos a cabo este experimento, resulta que justo antes de cada golpeo de un electrón en la pantalla obtenemos un destello luminoso que llega desde la rendija R1 o desde R2 , pero nunca de ambas a la vez, ni nada por el estilo. Debemos concluir, en contra de lo que dijimos anteriormente, que cada electrón pasa por una y sólo una de las rendijas, pero esto ya es desconcertante, pues hemos obtenido dos conclusiones contradictorias. Si hacemos el experimento como hemos dicho, con este otro detector, cerrando una y luego otra rendija, obtenemos los diagramas de probabilidad que obtuvimos en la figura 1.5. Ahora bien, £que obtenemos para la probabilidad total, con ambas rendijas?. Ya no tenemos lo que obtuvimos en esa figura, ahora tenemos algo como la figura 1.3, con proyectiles macroscópicos. Ahora que observamos los electrones detenidamente, parecen 9

A. Átomo de Rutherford

Capítulo 1. Orígenes de la Física Cuántica

comportarse bien, ya no generan interferencias. Y si apagamos la luz, volvemos a obtener interferencia. Como vimos, la luz efectivamente lleva un cantidad de movimiento, de valor p = h/λ. Podríamos intentar usar luz con una longitud de onda mayor, intentando evitar perturbar tanto los electrones, de modo que obtengamos interferencia, pero a la vez sepamos perfectamente por donde pasó el electrón. Pero no obtendremos éxito, pues encontraremos que es imposible conseguir que la luz no perturbe los electrones y no perturbar el diagrama de interferencia, y a la vez determinar con exactitud a través de cual de las dos rendijas ha pasado todos y cada uno de los electrones. Y básicamente en esto consiste el principio de incertidumbre de Heisenberg, que imponen de algún modo una limitación a nuestra capacidad experimental, no una limitación de nuestros aparatos o de nuestras capacidades, sino una limitación intrínseca a la naturaleza: hay cosas que no se pueden determinar con exactitud absoluta al mismo tiempo. Por el momento, nadie ha encontrado una manera de esquivar este principio de incertidumbre o de indeterminación propio de la naturaleza. Quedan algunas preguntas a las que no hemos respondido totalmente, y ahora responderemos a ellas, al menos, tan “totalmente” como podamos. Los electrones, £pasan por una y solo una de las dos rendijas, o siguen otros caminos alternativos?. Lo único que podemos decir para no caer en contradicción con los experimentos es que, cuando los observamos de cerca, el electrón pasa efectivamente por una y sólo una de las rendijas, y entonces se observa interferencia. Pero si no los observamos, si no colocamos un aparato que nos pueda decir esto, entonces no sabemos si el electrón atraviesa una rendija, las dos o ninguna, es más no sabemos ni siquiera que camino sigue de forma general, y se obtiene una interferencia. Por otro lado, si todo en la naturaleza tiene un carácter ondulatorio, £por qué al hacer el experimento con proyectiles grandes, sin ningún detector, no obteníamos interferencias?. Bueno, es que en realidad si las obteníamos, pero ni nosotros ni nuestros detectores éramos capaces de detectarlas, por eso creíamos que no teníamos interferencias. Las longitudes de las ondas asociadas a tales objetos son tan pequeñas que ningún detector sería capaz de distinguir los máximos y los mínimos, sino que presentaría en realidad un promedio, una curva de interferencia muy suavizada, algo parecido a lo que se representa en la figura 1.6.

A.

A.1.

Decaimiento del electrón al núcleo en el átomo de Rutherford Potencia radiada

De las leyes de Maxwell se deduce que una carga en movimiento acelerado irradia una energía por unidad de tiempo (potencia) que viene dada por la fórmula de Larmor. 10

A. Átomo de Rutherford

Capítulo 1. Orígenes de la Física Cuántica

Figura 1.6: Curvas con interferencias suavizadas Para el caso de la carga del electrón, esta es: dW e2 a2 = dt 6πε0 c3

(1.6a)

Aquí a representa la aceleración de la carga, c la velocidad de la luz, y e la carga del electrón. Usaremos esta fórmula más tarde para considerarla en la conservación de la energía. Por otro lado las fuerzas que intervienen son las de la propia atracción eléctrica (según la Ley de Coulomb) y la “fuerza centrífuga” para un movimiento circular, como se supone en este caso. Es decir: Fcent acent Fcoul

me v 2 = r v2 = r 1 e2 = − 4 πε0 r2

(1.6b)

El equilibrio que se supone implica que el balance de fuerzas sea nulo: Fcent + Fcoul =

me v 2 1 e2 − =0 r 4πε0 r2 11

(1.6c)

A. Átomo de Rutherford

Capítulo 1. Orígenes de la Física Cuántica

De aquí se deduce que 1 e2 (1.6d) 4πε0 r Además, la energía para una órbita de este tipo, análogamente a las órbitas planetarias, viene dada por: mv 2 =

1 e2 1 1 e2 1 e2 1 1 e2 1 = − =− E = me v 2 − 2 4πε0 r 4πε0 2 r 4πε0 r 4πε0 2 r

(1.6e)

Esta ecuación, derivada, nos da la variación de la energía orbital de la partícula: dE 1 1 e2 dr = dt 2 4πε0 r2 dt

(1.6f)

Conviene definir una constante, k0 , de la forma k0 =

1 4πε0

para mayor simplicidad en los cálculos posteriores. Si tenemos en cuenta que debe cumplirse: dW dE + =0 dt dt

(1.6g)

llegamos, dados los valores para a y v dados respectivamente por las condiciones (1.6b) y (1.6d), a la siguiente ecuación: −

3 m2e c3 2 r dr = dt 4 k02 e4

Ahora integramos esta ecuación: Z Z Z τ 0 3 m2e c3 0 2 3 m2e c3 r=r 2 − r dr = + 2 4 r dr = dt 4 k02 e4 r=r0 4 k0 e 0 0

(1.6h)

(1.6i)

Esto significa que: τ=

3 m2e c3 3 r 4 k02 e4 0

(1.6j)

Ahora tomaremos unos valores para los datos que aparecen en la ecuación (1.6j). Estos valores serán aproximados, pues tampoco queremos un valor exacto de este tiempo, sino simplemente demostrar que es de un orden de magnitud tan pequeño, que si este cálculo se correspondiera con la realidad, no existirían átomos estables. Así pues, aplicando dichos valores a la expresión (1.6j), tendremos que el tiempo en que se esperaría que el electrón precipite sobre el núcleo es: 12

A. Átomo de Rutherford

Capítulo 1. Orígenes de la Física Cuántica

Cuadro 1.1: Valores de las magnitudes Magnitud

Valor

Unidades

Masa del electrón me

9,1 · 10−31

kg

Velocidad de la luz c

3 · 108

m/s

Constante culombiana k0

9 · 109

Nm2 C−2

1,6 · 10−19

C

' 10−10

m

Carga del electrón e Radio de la órbita del electrón r0

τ ' 3 · 10−10

s

Efectivamente, el tiempo en que supuestamente el electrón caería sobre el núcleo no nos da un valor infinito, o similar, sino un valor concreto, finito, y muy pequeño, en contra de la evidencia experimental.

A.2.

Conclusión

De este sencillo cálculo se deduce que aunque las predicciones hechas por Rutherford acerca de la estructura de la materia son en general acertadas, el modelo que se plantea según la física clásica es incompleto, esto es, necesita ser revisado con una nueva teoría, que es la mecánica cuántica. La hipótesis de Bohr acerca de la existencia de órbitas estables cuantizadas (no puede existir cualquier órbita) en las cuales se puede localizar el electrón, y la emisión discreta de energía en forma de cuantos al pasar aquellos de una órbita a otra puso de manifiesto, junto a otros hechos experimentales, la validez de la mecánica cuántica a nivel atómico.

13

Capítulo 2

Mecánica ondulatoria. Ecuación de Schrödinger 1.

Paquetes de ondas Ahora veremos algunas propiedades generales de las ondas, no importará de que tipo.

La exponencial compleja eikx es una onda plana. Según la fórmula de Euler, eikx = cos kx + i sin kx. kx ha de ser adimensional, luego [k] = L−1 . k es el número de ondas, es inversamente a la longitud de onda. El módulo de la exponencial compleja ikx 2 proporcional −ikx ikx es 1, e = e e = 1. El uso de números complejos es una conveniencia matemática útil en algunas ramas de la física clásica, pero es que resulta ser fundamental, una propiedad intrínseca, en la física cuántica1 . Una onda plana está totalmente extendida. Si sumamos dos ondas (o un número finito de ellas) de distinta frecuencia, parece que se vayan localizando. Para localizarlas totalmente, sumamos un número infinito de ellas, es decir, las integramos. Formamos entonces un paquete de ondas. Tenemos: Z +∞ g(k)eikx dk (2.1) −∞

En la anterior expresión, g(k) es una función de peso, que tomará valores en una zona determinada, normalmente (Esto hará que se pueda encontrar la onda localizada). Lo más normal es que g(k) sea una gaussiana, teniendo entonces un paquete de ondas gaussiano: Z +∞

2

e−α(k−k0 ) eikx dk

Ψ(x) =

(2.2)

−∞ 1

Por ejemplo, el estado de un sistema se describe mediante un vector de un espacio vectorial cuyo cuerpo de escalares es el plano complejo, C

14

2. Ecuación de Schrödinger

Capítulo 2. Ecuación de Schrödinger

La función Ψ(x) es una función de onda (no necesariamente cuántica). La constante α estará relacionada con la anchura del paquete gaussiano. Para hacer esta integral, hacemos la siguiente transformación: k − k0 = k 0 , por lo que dk = dk 0 , y queda: Z +∞ 02 0 ik0 x e−αk eik x dk 0 Ψ(x) = e −∞

Como aparece un término en k 0 y otro en k 02 , lo que intentamos es completar el cuadrado perfecto. Z +∞ ix 2 x2 0 ik0 x Ψ(x) = e e−α(k − 2α ) e− 4α dk 0 −∞

El término en x vuelve a salir, y hacemos otro cambio de variables, k 0 − dk 0 = dk 00 , y obtenemos: r π ik0 x − x2 Ψ(x) = e e 4α α

ix 2α

= k 00 , luego (2.3)

Por tanto: π −x2 |Ψ(x)|2 = e 2σ α p √ var(x) = σx = α

2.

Ecuación de Schrödinger

Vimos en el tema anterior que la función de onda para el electrón libre viene dada por:   Z ∞ i(px − Et) dp ψ(x, t) = g(p) exp ~ −∞ En general, se tendrá exp(ikx − wt), solo que hemos sustituido las propiedades particulares (cuánticas) que hemos visto. Por motivos de interés matemático, y sin ninguna implicación física, hacemos:   Z ∞ 1 i(px − Et) ψ(x, t) = √ φ(p) exp dp ~ 2π~ −∞ Así conseguiremos que ψ(x, t) y φ(p, t) serán las funciones de onda en x y p, respectivamente.

15

2. Ecuación de Schrödinger

Capítulo 2. Ecuación de Schrödinger

Vamos a ver la ecuación diferencial que satisface la anterior ecuación:   Z ∞ iE ∂ψ 1 i(px − Et) − φ(p) exp = √ dp ∂t ~ 2π~ −∞ ~   Z ∞ ip ∂ψ 1 i(px − Et) = √ φ(p) exp dp ∂x ~ 2π~ −∞ ~   Z ∞ 2 p i(px − Et) ~2 ∂ 2 ψ 1 φ(p) exp dp − = √ 2m ∂x2 ~ 2π~ −∞ 2m

Si consideramos que el electrón tiene solo energía cinética, E = p2 /2m, se cumple entonces que: ~2 ∂ 2 ψ ∂ψ =− (2.4) i~ ∂t 2m ∂x2 Schrödinger generalizó2 2.4, pensando que debería existir también un término que contribuyera con energía potencial, V , de la forma: i~

∂ψ(x, t) ~2 ∂ 2 ψ(x, t) =− + V (x)ψ(x, t) ∂t 2m ∂x2

(2.5)

De esta ecuación son de destacar varios hechos: La energía total aparece como una derivada temporal de un función, la función de onda. El momento, sin embargo, es introducido como una variación espacial de dicha función. El tiempo, t, es un simple parámetro, se introduce exactamente igual que en la física clásica, y no como en teorías como la Relatividad, general o especial, en las que aparece como una coordenada más, intrínsecamente relacionado con la estructura del espacio-tiempo. Además, estos pasos que hemos presentado no son una demostración, sino una motivación, si acaso. La ecuación de Schrödinger no se deduce, se postula. 2

En realidad publicó primero una versión relativista de la ecuación, la de Klein-Gordon ??

16

3. Interpretación probabilística

3.

Capítulo 2. Ecuación de Schrödinger

Interpretación probabilística de la mecánica cuántica

La interpretación estadística de la mecánica cuántica3 data de los inicios de ésta. Hemos visto que toda estado de una partícula viene definido por una función de onda, que representa su carácter ondulatorio. £Cómo se da cuenta del comportamiento material?. Para esto se interpreta la cantidad: ρ(r, t) = |ψ(r, t)|2 = ψ(r, t)∗ ψ(r, t)

(2.6)

como la densidad de probabilidad de encontrar la partícula en el instante de tiempo t en la posición r. La función de onda permite pues determinar completamente el estado de un sistema, tanto su carácter material como su carácter ondulatorio. Este concepto es, sin duda, una de las grandes revoluciones de la mecánica cuántica: no somos capaces de afirmar que la partícula esté en un punto concreto en un cierto instante, sino que lo único que sabemos – y no siempre – la probabilidad de que una partícula esté aquí y la probabilidad de que esté allá4 . Vamos a demostrar ahora que la probabilidad total se conserva. Para ello vamos a establecer una ecuación de continuidad de forma análoga a las que encontramos en electromagnetismo o en la mecánica de los fluidos. Definimos la corriente de probabilidad como: j(r, t) =

i~ (ψ(r, t)∇ψ ∗ (r, t) − ψ ∗ (r, t)∇ψ(r, t)) 2m

(2.7)

Para obtener la evolución temporal de la densidad de probabilidad debemos hacer su derivada (omitimos las dependencias temporales por simplicidad), y tener en cuenta que la evolución temporal de la función de onda viene dada en función de su hamiltoniano por la ecuación (2.5), donde las derivadas de segundo orden en x han de ser generalizadas 3

Veáse el discurso de Max Born al recibir el premio Nobel (11 de Diciembre de 1954): “The work, for which I have had the honour to be awarded the Nobel Prize for 1954 contains no discovery of a fresh natural phenomenon, but rather the basis for a new mode of thought in regard to natural phenomena.” “El trabajo por el que he tenido el honor de recibir el Premio Nobel de 1954 no presenta el descubrimiento de un nuevo fenómeno natural, sino más bien continene las bases de un nuevo modo de pensamiento en relación a los fenómenos naturales.”

4 Este concepto tan anti-intuitivo es lo que produjo más estupor entre los científicos de la época, hasta el punto de que el mismo Einstein, que propuso la necesaria cuantización de la radiación electromagnética en fotones, escribiera su famosa frase, asegurando que “Dios no puede estar jugando a los dados”, a lo que Bohr respondió: “Deja de decirle a Dios lo que tiene que hacer”.

17

4. Operadores. Valores esperados

Capítulo 2. Ecuación de Schrödinger

a operadores laplacianos:  
∂ρ ∂ψ ∗ ~2 ∗ ∂ψ i~ = i~ ψ +ψ =− ψ ∗ ∇2 ψ − ψ∇2 ψ ∗ ∂t ∂t ∂t 2m

(2.8)

Ahora bien, de la definición de j, y teniendo en cuenta que f ∇2 f ∗ − f ∗ ∇2 f = ∇ · (f ∇f ∗ − f ∗ ∇f ), podemos llegar a la ecuación: ∂ρ(r, t) + ∇ · j(r, t) = 0 (2.9) ∂t Esta es la ecuación de continuidad que andábamos buscando. Podemos integrarla sobre un volumen V encerrado en una superficie (real o imaginaria) ∂V para obtener: Z Z Z Z ∂ρ(r, t) ∂ρ(r, t) dv + ∇ · j(r, t)dv = 0 = dv + j(r, t) · ds = 0 ∂t ∂t V V V ∂V Esto nos permite interpretar j(r, t) como la probabilidad por unidad de tiempo de que la partícula atraviese la superficie ∂V . Así pues, podemos visualizar la partícula como si fuera un fluido con densidad ρ(r, t), que, como tal, cumple una/la ecuación de continuidad. Ahora, para demostrar, como anunciábamos, que la probabilidad total se conserva, pasaremos al caso unidimensional, por ser más cómodo e intuitivo, aunque la solución será general. Para ello, escribimos la ecuación de continuidad en una dimensión (x) y la integramos: Z Z b ∂j ∂ b ρ(x, t)dx = − dx = j(a, t) − j(b, t) t a a ∂x Pero, en el límite en que a → −∞ y b → ∞, como la función de onda no puede diverger (se debe hacer cero en los infinitos), tenemos que j(−∞, t) = j(∞, t) = 0, y por tanto: Z ∞ Z ∞ 2 |ψ(x, 0)| dx = 1 ⇒ |ψ(x, t)|2 dx = 1 ∀t −∞

−∞

Este hecho, como hemos dicho varias veces, se puede enunciar diciendo que la probabilidad total se conserva.

4.

Operadores. Valores esperados

Un operador es un objeto matemático que actúa sobre un vector devolviendo otro. En el caso de la mecánica cuántica, dichos operadores serán funciones, las funciones de onda5 Para calcular magnitudes, es decir, para obtener algún resultado tangible (un 5

Todo esto se verá con mayor detenimiento en el capítulo 4. Aquí solo presentaremos algunos resultados útiles que nos permitan trabajar sobre ciertos problemas sencillos, como haremos en el capítulo 3.

18

4. Operadores. Valores esperados

Capítulo 2. Ecuación de Schrödinger

número, algún valor de una cantidad física), £qué debemos hacer en mecánica cuántica?. Bien, las magnitudes llevan asociados operadores, que actúan sobre las funciones de onda. Un ejemplo de magnitudes calculables son los valores medios. Por ejemplo el valor medio de la magnitud x, la posición, que lleva asociada el operador X, que consiste simplemente en multiplicar la función por x, vendía dado por: Z ∞ Z ∞ 2 ψ ∗ (x)xψ(x)dx (2.10) x |ψ(x)| dx = hx i = −∞

−∞

Obsérvese la abismal diferencia entre la física clásica y la nuevea física, la física cuántica. Con la operación anteriormente indicada, una integral de una función, tan complicada como se requiera en cada caso, multiplicada por x, nos da simplemente un valor medio de x. En la física clásica, una simple sustitución de algún parámetro, como el tiempo, nos daría inmediatamente el valor de la posición exacto, no en promedio. En general, los operadores no consistirán en una simple multiplicación, sino que implicarán, dentro de este ámbito, derivadas, integraciones, etc. Entonces, el valor medio de una magnitud no se calcularía como hemos anunciado anteriormente, sino de una forma similar, pero con una diferencia sutil. El valor medio de la magnitud A, cuyo operador asociado es A, viene dado por: Z ∞ ψ ∗ (x, t)Aψ(x, t)dx (2.11) hA i = inf ty

Hasta ahora hemos visto, sin darnos cuenta de ello, dos operadores: el operador posición, x, y el operador energía potencial, V (x). Hemos trabajado con ellos sin saber nada sobre como operar con ellos precisamente porque solo implican multiplicar la función de onda por las magnitudes necesarias. El siguiente operador que vamos a ver, sin embargo, ya no es tan trivial; se trata del operador momento, p, que está relacionado con el momento cinético clásico. Su expresión la podemos intuir de la forma general del hamiltoniano (viendo cual es la aportación de la energía cinéticaa éste), según la ecuación (2.5): −

p2 ~2 ∂ 2 ψ(x, t) ≡ ψ(x, t) 2m ∂x2 2m

o viendo que (como caso particular) para una onda plana:   i(px−Et) ∂exp ~ i ~ ∂ = p exp(. . .) ⇒ p = ∂x ~ i ∂x Como hicimos con la ecuación de Schrödinger, debemos señalar que esto no es ni mucho menos una deducción, solo una orientación, o una motivación. Definimos, en general, el operador momento (para la variable x, solamente) como: p=

~ ∂ i ∂x

operador momento lineal 19

(2.12)

4. Operadores. Valores esperados

Capítulo 2. Ecuación de Schrödinger

Así, por ejemplo, el valor medio del momento lineal es: Z ∞ ~ ∂ ψ ∗ (x, t) (ψ(x, t))dx hp i = i ∂x −∞ De los operadores de la energía, ya conocemos el operador energía potencial, y tenemos cierto conocimiento del operador hamiltoniano, pero no conocemos la energía cinética. Por la analogía con la física clásica, donde Ec = π 2 /2m, con π el momento cinético, se tendrá que el operador energía cinética está dado por el término 1/2m multiplicado por p2 . Pero ahora p es un operador. Surge la pregunta £qué es el “cuadrado” de un operador?. El “cuadrado” de un operador es este mismo actuando sobre el resultado de hacer actuar dicho operador sobre una función, por ejemplo, T 2 (ψ) = T (T (ψ)). Así, para el caso de p, se tiene: ~2 ∂ 2 p2 Ec = =− 2m 2m ∂x2 Podemos escribir la ecuación de Schrödinger de forma más compacta definiendo H, el operador hamiltoniano, como: ~2 ∂ 2 + V (x) 2m ∂x2 de modo que la ecuación fundamental de la mecánica cuántica queda como: H=−

i~

∂ψ(x, t) = Hψ(x, t) ∂t

Ecuación de Schrödinger

(2.13)

Es importante observar que la dependencia temporal está contenida de forma explícita en la función de onda, y no en los operadores. Así, los valores medios evolucionarán con el tiempo no por la magnitud que representan, sino por como dependa del tiempo la función de onda. Es por eso que, macroscópicamente, encontramos una depencia de x, por ejemplo, con el tiempo: es porque lo que en realidad medimos en la física clásica son promedios temporales de las magnitudes. Por otro lado, como seguramente habrá percibido el lector, parece posible que el valor esperado de una magnitud sea un número complejo, atendiendo a la definición dada en (2.11). Esto no sucede así en la realidad, y se debe a que todo operador asociado a una magnitud física observable es hermítico, es decir, cumple la relación: Z b ∗ Z b Z b ∗ ψ2 Aψ1 dx = ψ1 Aψ2 dx = ψ1 (Aψ2 )∗ dx (2.14) a

a

a

Es directo entonces comprobar que el valor esperado de todo operador hermítico es real. Para ello, tomamos sobre la definición de hermiticidad las funciones ψ1 = ψ2 = ψ, por lo que tenemos: Z ∞ ∗ Z ∞ ∗ ∗ hA i = ψ (x, t)Aψ(x, t)dx = ψ (x, t)Aψ(x, t) = hA i∗ −∞

−∞

20

4. Operadores. Valores esperados

Capítulo 2. Ecuación de Schrödinger

Es decir, se cumple hA i = hA i∗ , por lo que hA i es real si A es hermítico. Como ejercicio, demuestrese que los operadores x, V y p son hermíticos. Téngase en cuenta que los operadores x y V solo involucran multiplicaciones por números reales. En general, se observa que son hermíticos operadores en los que aparezcan: (multiplicaciones por) Funciones reales. ∂ Expresiones del tipo i ∂x , por la aparición de i.

Expresiones del tipo

∂2 , ∂x2

a través de la integración por partes.

Por último, aparece entre los operadores una relación muy importante, la realación de conmutación. Es evidente que en esta mecánica que estamos desarrollando los operadores no tienen por qué conmutar entre sí, es decir, no es indiferente el orden en que se apliquen dos operadores sobre una misma función. Esta mecánica de funciones de onda es la mecánica de Schrödinger. Existe otra formulación, desarrollada en paralelo y al mismo tiempo que aquella, la formulación matricial de Heisenberg, que se demuestra totalmente equivalente a la primera, en la que también aparece de forma natural la no conmutación de los operadores, ya que en esta vienen dados por matrices, que, como se sabe de las nociones fundamental del Álgebra Lineal, no conmutan entre sí. Se define el conmutador de dos operadores A y B, como: [A, B] = AB − BA

(2.15)

El conmutador es un nuevo operador. Aplicado sobre una función, se tiene: [A, B] ψ = A(Bψ) − B(Aψ) Como ejemplo, calculemos el conmutador de los operadores posición y     ∂ψ ∂xψ ~ ∂ψ ∂x ∂ψ ~ x − = x −ψ −x = [x, p] (ψ) = i ∂x ∂x i ∂x ∂x ∂x

momento: ~ (−ψ) i

Por tanto, se tiene: [x, p] = i~

(2.16)

De esta relación tan simple depende gran parte de la teoría de la mecánica cuántica, y del asombro que esta produce, pues sobre ella está fundamentado el famoso principio de incertidumbre de Heisenberg.

21

5. Representación de momentos

5.

Capítulo 2. Ecuación de Schrödinger

Representación de momentos

5.1.

Función de onda en p como transformada de Fourier

5.2.

Teorema de Parseval

6.

Autovalores y autovectores

Supongamos que tenemos un operador A, asociado a una magnitud física A. Supongamos que, a aplicarlo sobre una función, se tiene lo siguientes: Aψn,i (x) = λn · ψn,i (x)

(2.17)

Las funciones que cumplen la anterior relación, para algún valor λ, se llaman autofunciones o funciones propias del operador A, así como los λ se llaman autovalores o valores propios. Pueden existir muchos autovalores de un operador A (por ello colocamos un subíndice n en los autovalores). Además, el conjunto de los autovalores puede ser continuo o degenerado. Para un mismo valor de λ pueden existir varias funciones que cumplan (2.17): decimos entonces que el espectro de A es degenerado (Por esto hay un subíndice i en las funciones). Hay varias propiedades generales, que serán demostradas en el capítulo 4, que resultan interesantes y que por ello adelantamos aquí: Los valores propios de un operador hermítico son reales. Si un operador tiene dos o más valores propios distintos, se cumple que las funciones asociadas a cada uno de dichos autovalores son ortogonales entre sí (incluso las que corresponden a un mismo autovalor). Es decir, λn 6= λm ⇒ hψn |ψm i (La expresión h|i denota el producto escalar). El conjunto de las funciones propias de casi6 cualquier operador hermítico forma una base en el espacio vectorial correspondiente. Los postulados de la mecánica cuántica afirman que, al realizar una medida de un observable, los posibles valores que se pueden obtener son los valores propios del operador asociado. Si tenemos un operador cuyas funciones propias forman una base, podemos escribir cualquier función de estado del sistema en esa base de la forma: X ψ(x) = cn ψn (x) n 6

En los espacios vectoriales de dimensión infinita, que aparecen a menudo en mecánica cuántica, existe una sutil diferencia entre operadores hermíticos y operadores auto–conjugados, y de este hecho proviene este casi.

22

7. Evolución temporal

Capítulo 2. Ecuación de Schrödinger

de manera que los coeficientes cn adquieren una especial relevancia: están relacionados con las probabilidades de obtener los distintos valores propios (recordemos, los valores posibles de una medida) del operador asociado a la base en la que está desarrollada la función de onda del sistema. Es decir: P(λn ) = |cn |2 De nuevo nos encontramos con el hecho de que en mecánica cuántica solo podemos dar probabilidades de que un cierto acontecimiento suceda. Vemos aquí una nueva manera de obtener valores esperados. En la sección 4 estudiábamos como obtener valores esperados para una distribución de probabilidad continua. Con la ecuación anterior tenemos una distribución de probabilidad discreta. Sustituyendo las integrales por sumas, podemos obtener el valor esperado de la magnitud A de la siguiente manera: X hA i = λn |cn |2 (2.18) n

Esta ecuación sería la análoga a la ecuación (2.10).

7.

Evolución temporal

Vamos a estudiar aquí la solución más general de la ecuación de Schrödinger. Partimos de la ecuación: ∂ψ(x, t) ~2 ∂ 2 ψ(x, t) i~ =− +V (x)ψ(x, t) Ecuación de Schrödinger independiente del tiempo ∂t 2m ∂x2 (2.19) Lo primero que podemos hacer es intentar probar con un cambio de variables, de la forma: ψ(x, t) = ξ(t)ϕ(x) Si aplicamos lo anterior a la ecuación y dividimos por la función dada en la forma anterior, obtenemos: ~2 1 d2 ϕ(x) i~ ∂ξ(t) =− + V (x) ≡ cte. = En ξ(t) ∂t 2m ϕ(x) dx2 Como el primer miembro depende exclusivamente de t y el segundo de x, es lógico que para que ambos términos sean iguales deberán ser una constante, que no dependa ni de x ni de t. Esa constante viene dada por la energía E. Esto se debe a que los autovalores del operador hamiltoniano, como se muestra en la siguiente ecuación, vienen dados por las posibles energías del sistema: Hψ(x, t) = En ψ(x, t)

(2.20)

Hemos pasado de tener una ecuación con dos incógnitas a tener dos ecuaciones de una incógnita. La ecuación que queda para ξ(t) tiene como solución: ξ(t) = Ce− 23

iEn t ~

7. Evolución temporal

Capítulo 2. Ecuación de Schrödinger

Para ϕ(x) no sabemos nada, solo que en general será una solución de la ecuación de autovalores: Hψn (x) = En ψn (x) Ecuación de Schrödinger independiente del tiempo

(2.21)

Esta es la ecuación de Schrödinger independiente del tiempo. La solución mas general vendrá dada por el una combinación lineal del producto de estas dos funciones, combinación lineal que dependerá del índice n: ψ(x, t) =

X

αn e−

n

24

iEn t ~

ϕn (x)

(2.22)

Capítulo 3

Problemas unidimensionales 1. 1.1.

Introducción. Generalidades Ecuación de Schrödinger para potenciales constantes. Continuidades y discontinuidades

Vamos a estudiar ahora problemas unidimensionales, para resolver la ecuación de Schrödinger en casos sencillos, ya que fuera de estos es de difícil solución analítica. Además muchos de estos casos sencillos son modelos que sirven para aproximar casos físicos reales, como los núcleos atómicos. El hecho de que se estudien en una sola dimensión es porque son de fácil generalización a más dimensiones, y también porque a veces la dependencia es efectivamente unidimensional (radial, por ejemplo). Si consideramos la ecuación de Schrödinger en varias dimensiones, dependiente del tiempo (por ejemplo, (2.5)), al aplicarla al caso unidimensional tendremos que el operd2 ador ∇2 pasa a ser simplemente dx 2 . Además buscaremos soluciones estacionarias de la ecuación, es decir, independientes del tiempo. Tenemos pues la siguiente ecuación:

i~

~2 d2 ψ(x, t) ∂ψ(x, t) + ψ(x, t)V (x) =− ∂t 2m dx2

(3.1)

Si introducimos una separación de variables en la anterior ecuación de la forma ψ(x, t) = ξ(t)ϕ(x), y dividimos toda la ecuación por ψ(x, t) en la forma dada, se obtiene (omitimos las dependencias respecto a las variables por simplicidad): i~ ∂ξ ~2 1 d2 ϕ =− + V (x) ≡ cte. ξ ∂t 2m ϕ dx2 La dos miembros de la ecuación anterior tienen que ser iguales a una constante porque el miembro de la izquierda depende solo de t, mientras que el segundo depende solo de 25

1. Introducción. Generalidades

Capítulo 3. Problemas unidimensionales

x. Así pues, llamando a esa constante E, que será la energía (ahora veremos por qué), tenemos:

−

~2 d2 ϕ(x) + V (x)ϕ(x) = Eϕ(x) 2m dx2

(3.2)

Esta es la ecuación de Schrödinger unidimensional independiente del tiempo. Decimos que esa constante E representa la energía porque es el autovalor (en realidad podrá haber b = Eϕ, que es una ecuación de autovalores, y donde H b es el varios) de la ecuación Hϕ operador hamiltoniano. Todo esto se verá con más detalle y quedará mejor entendido en el capítulo dedicado a la estructura y postulados de la mecánica cuántica. Si consideramos, como haremos en este capítulo, potenciales en los que, de una manera u otra, el potencial es constante, V (x) = V , podemos resolver directamente la ecuación, considerando los casos en los que E > V ó E < V . Podemos escribir la ecuación de la forma (omitiendo las dependencias1 ): d2 ψ 2m + 2 (E − V )ψ = 0 dx2 ~

(3.3)

Según los casos, la solución será: E − V > 0 En este caso la solución será una combinación lineal de exponenciales imaginarias: ψ(x) = A exp (ikx) + B exp (−ikx) (3.4a) definiendo k como:

r k=

2m (E − V ) ~2

(3.4b)

E − V < 0 En este otro caso la solución será una combinación lineal de exponenciales normales: ψ(x) = A exp (qx) + B exp (−qx) (3.5a) donde hemos definido q como: r q=

−

2m (E − V ) ~2

(3.5b)

También es importante como ha de ser la función de onda en las discontinuidades. Deberá de ser siempre continua, y su derivada también2 , salvo en el caso de un potencial 1

Escribiremos a partir de ahora la parte espacial de la función de onda como ψ = ψ(x), y no como ϕ, como hicimos antes en la separación de variables, y no tomaremos ψ como la función de onda total, ψ = ψ(x, t) 2 Una forma de entender esto es que, si no lo fuera, tendríamos que su derivada segunda tendría un valor infinito, mientras que en la ecuación de Schrödinger aparece igualada a algo finito.

26

1. Introducción. Generalidades

Capítulo 3. Problemas unidimensionales

delta (definido mediante una delta de Dirac, como veremos más adelante). Esto se puede ver integrando la ecuación de Schrödinger en las cercanías de una discontinuidad, situada arbitrariamente, por comodidad en la notación, en x = 0. Se tiene, para el caso de un potencial constante, V (x) = V :     Z ε Z ε 2 dψ(x) dψ(x) d ψ(x) ε→0 (E − V )ψ(x)dx −→ 0 ⇒ = dx = − dx2 dx dx −ε −ε ε −ε Sin embargo, si el potencial viene definido no por una constante sino por una función delta, de intensidad λ, V (x) = λδ(x), se tendrá, usando los mismos cálculos que en el apartado anterior:     Z ε dψ(x) 2mλ dψ(x) 2mλ − = δ(x)ψ(0)dx = 2 ψ(0) 2 dx dx ~ −ε ~ +ε −ε Esto deriva en la siguiente condición, para una discontinuidad situada en un punto a:     dψ(x) dψ(x) 2mλ (3.6) − = 2 ψ(a) dx dx ~ a+ε a−ε Otra cosa importante que tendremos que tener en cuenta serán las condiciones de contorno, es decir, las condiciones específicas que tendrá que cumplir la ecuación de onda, que serán precisamente lo que dará origen a fenómenos interesantes, tales como la cuantización de la energía y la aparición de estados discretos, el efecto túnel, resonancias, estados ligados, etc.

1.2.

Simetrías. Operador paridad

Definimos el operador paridad, Π, mediante la siguiente ecuación: Π | ψ(x) i = | ψ(−x) i Para ver los posibles autovalores de Π en la ecuación Πϕ(x) = λπ | ϕ i sólo tenemos que hacer actuar dos veces sobre un vector propio ϕ(x): √ Π2 | ϕ(x) i = λ2π | ϕ i = Π | Πϕ(x) i = Π | ϕ(−x) i = | ϕ(x) i ⇒ λπ = 1 = ±1 Por otro lado, veremos más adelante (capítulo 4) que la variación del promedio temporal de una magnitud asociada a un operador A varía como ∂t hA i ∝ [H, A]. Es decir, una cantidad se conserva si su operador asociado conmuta con el hamiltoniano. Para el caso de nuestro operador Π se tiene que:  2  P ,Π = 0 27

2. Potencial escalón

Capítulo 3. Problemas unidimensionales 

[V (x), Π] =

= 0, si V (x) es par, V (x) = V (−x); 6= 0, si V (x) no es par.

Es decir, si el potencial es par, V (x) = V (−x), entonces [H, Π] = 0, y la paridad se conserva. Esto es importante, porque, si el potencial es par, permite elegir como funciones propias del hamiltoniano funciones que sean pares o impares, que permitirán resolver con más facilidad los distintos problemas, según como sea el espectro de H: Estados no degenerados La función ha de ser necesariamente impar o impar. Estados degenerados Se pueden elegir funciones propias pares o impares (y con ellas establecer una combinación lineal para obtener una solución general)

2.

Potencial escalón

Vamos a estudiar el comportamiento de una partícula en un potencial escalón (véase figura 3.1), definido matemáticamente de la siguiente forma:  0 si x < 0; V (x) = V0 si x > 0. I

V (x)

6

II

V0

O -

x

Figura 3.1: Potencial escalón Habrá que considerar soluciones en uno y otro lado del salto de potencial, y después unirlas. Según que la energía sea E < V0 ó E > V0 tendremos distintos tipos de solución.

2.1.

E > V0 . Reflexión parcial

En este caso, la solución en uno y otro lado del salto de potencial es exponencial imaginaria, y escribimos: ψI (x) = AI eikI x + BI e−ikI x ψII (x) = AII eikII x + BII e−ikII x 28

2. Potencial escalón r kI =

2mE ~2

kII

Capítulo 3. Problemas unidimensionales r 2m(E − V0 ) = ~2

Elegimos arbitrariamente BII = 0. Esto lo hacemos queriendo indicar que no hay ondas que provengan desde +∞. Estableciendo las condiciones de que la función de onda y su derivada sean continuas en x = 0 obtenemos estas condiciones: AI + BI = AII ikI (AI − BI ) = ikII AII De aquí, y atendiendo a las definiciones de kI y kII , se tiene:

AII AI BI AI

2kI kI + kII kI − kII = kI + kII

=

Ahora vamos a definir dos magnitudes y conceptos que usaremos con frecuencia. El coeficiente de transmisión, T , representa la probabilidad de que una onda incidente sea transmitida a través del salto de potencial, mientras que el coeficiente de reflexión, R, es la probabilidad de que una onda incidente se refleje. Si tenemos en cuenta como hemos escrito las funciones de onda, podemos definirlos de la siguiente manera: 2 kII AII T = k I AI 2 BI R = AI

kII 4kI2 kI (kI + kII )2 (kI − kII )2 R = (kI + kII )2 T =

Es directo comprobar que se cumple T +R = 1, es decir, la onda incidente ciertamente se transmite o se refleja. Esto es un resultado que resulta chocante en comparación con la física clásica: un partícula que llegase a un salto de potencial de esta manera, siendo su energía mayor que el potencial, jamás rebotaría, siempre sería transmitida. Cuánticamente cabe la posibilidad de que la onda/partícula se refleje. En la figura 3.2 hay una representación del coeficiente T frente a α, donde hemos definido E − V0 = E(1 − α), es decir, α = V0 /E. 29

2. Potencial escalón

Capítulo 3. Problemas unidimensionales

Figura 3.2: Representación del coeficiente de transmisión T en función de α = V0 /E.

2.2.

E < V0 . Reflexión total

En este caso se tiene que la solución en un lado (izquierdo) se corresponde con combinación de exponenciales imaginarias, mientras que en el otro la solución viene dada por una exponencial real. Escribimos: ψI (x) = AI eikx + BI e−ikx ψII (x) = AII eqx + BII e−qx r r 2mE 2m(E − V0 ) q= − k= 2 ~ ~2 Para que la función de onda sea normalizable, y la integral no diverja, es necesario hacer AII = 0. En estas circunstancias, y estableciendo las misma condiciones de empalme que antes, se llega a un resultado muy parecido al anterior (solo hay que cambiar kI → k y kII → iq):

BII AI BI AI

2k k + iq k − iq = k + iq =

Si intentamos hallar el coeficiente de reflexión, R, y operamos con él, obtenemos R = 1. Esto quiere decir que la partícula siempre se refleja. Sin embargo, dado que la función de onda es no nula en la zona II, existe una cierta probabilidad de encontrar a la onda “dentro” del escalón. Es inmediato comprobar que la corriente de probabilidad se conserva. Si definimos je 30

3. Pozo de potencial finito

Capítulo 3. Problemas unidimensionales

como la corriente entrante, y js como corriente saliente, se tendrá: ~k (|AI |2 + |BI |2 ) m ~q 2 = |C| m

je = js

Si operamos y atendemos a los resultados obtenidos anteriormente:   (k − q)2 ~q 4k 2 ~q 2 ~k je = |AI | 1− = = |AII |2 = js . 2 2 m (k + q) m (k + q) m Además, se puede observar, teniendo en cuenta los resultados de los coeficientes que acompañan a las exponenciales complejas, que la reflexión produce un cierto desfase. Si hallamos las constantes en el límite V0 → ∞ encontramos: BI V0 →∞ −→ −1 ⇒ BI −→ −AI AI BII V0 →∞ −→ −1 ⇒ BII −→ 0 AI Esto es, cuando el potencial se hace muy grande, toda la onda se refleja, si desfase alguno, y además no hay probabilidad de que la onda/partícula penetre en el escalón.

3.

Pozo de potencial finito

Vamos a estudiar un pozo cuadrado de potencial, cuya definición matemática está dada en (3.7). Una representación esquemática viene dada en la figura 3.3.  0, si |x| ≥ a; V (x) = (3.7) −V0 , si |x| ≤ a .

3.1.

−V0 < E < 0. Estados ligados

Las soluciones de este caso vienen dadas por: ψ1 (x) = A exp (qx) + B exp (−qx) ψ2 (x) = C exp (ikx) + D exp (−ikx) ψ3 (x) = F exp (qx) + G exp (−qx) 31

3. Pozo de potencial finito

Capítulo 3. Problemas unidimensionales

V =0 I

V =0 II

−a

III a

−V0

Figura 3.3: Pozo cuadrado de potencial

r

r 2m|E| 2m(V0 − |E|) q= k = ~2 ~2 Para evitar la divergencia en la normalización tendremos que tomar G = A = 0.

(3.8)

A los estados que aparecerán en la zona II se les llama estados ligados: el pozo cuadrado atrae, liga a los electrones. Si éstos van con mucha energía pasan por él sin ser muy perturbados, pero si van con poca energía, pueden quedar ligados al pozo. La complicación que puede aparecer al solucionar el sistema de ecuaciones anterior se resuelve mediante las consideraciones de paridad del potencial. Funciones propias pares En este caso se tiene B = F y C = D. Esto implica que las consideraciones de empalme serán las mismas en −a y en a. Se tiene: C (exp (ika) + exp (−ika)) = B exp (−qx) ikC (exp (ika) − exp (−ika)) = −qB exp (qx) Las anteriores ecuaciones nos dicen implícitamente para que energías hay solución. Atendiendo a la fórmula de Euler sobre la exponencial compleja, se tiene la siguiente condición: q condición de estado ligado (3.9) tan (ka) = k Pero la anterior ecuación es una ecuación implícita, no podemos despejar analíticamente E para saber para que energías existe solución. Vamos a hallarlas mediante un método gráfico. Atendiendo a la definición de k y q, en la ecuación (3.8), y si hacemos estas otras definiciones: 2mV0 2 y =k·a λ = a2 (k 2 + q 2 ) = a ~2 32

3. Pozo de potencial finito

Capítulo 3. Problemas unidimensionales

llegamos fácilmente a la expresión: p λ − y2 tan(y) = y

(3.10)

Los niveles de un pozo dependen del producto de a2 y V0 , es decir, de la profundidad energética y la anchura espacial del pozo. Funciones propias impares Si realizamos las mismas operaciones que en el caso anterior, obtenemos: p λ − y2 cot(−y) = y

(3.11)

Representación Ahora representamos las posibles soluciones, para soluciones pares e impares, en una única gráfica. En√ realidad, en la gráfica están representadas las funciones cot(−πx), √ n2 −x2 λ tan(πx), f (x) = , hemos hecho y = πx, y n = π . Las soluciones –los posibles esx tados ligados– vienen dados por los cortes de cualquiera de las funciones trigonométricas con la función f (x). Podemos ver que el número de soluciones es: √ ! 2 λ = 1 + int(2n) ν = 1 + int π

3.2.

E > 0 > −V0 . Coeficientes de transmisión

Lo que se pretende en este caso es calcular los coeficientes de transmisión y reflexión. Tenemos las funciones de onda:

ψ1 = A exp(ikx) + B exp(−ikx) ψ2 = C exp(iqx) + D exp(−iqx) ψ3 = E exp(ikx) + F exp(−ikx) r k=

2mE ~2

r q= 33

2m(E − |V0 | ~2

3. Pozo de potencial finito

Capítulo 3. Problemas unidimensionales

Figura 3.4: Representación de estados ligados en un pozo cuadrado En este caso hay que usar también la paridad, pero es algo más complicado. Podemos elegir funciones propias del hamiltoniano pares o impares, pero no quiere decir que todas las soluciones tengan que serlo. Lo mejor será encontrar soluciones pares e impares, y establecer una combinación lineal que nos dé la solución general. Funciones propias pares En este caso, consideramos A = G = A+ = 1, B = F = B+ , y C = D = C+ . Las condiciones de empalme en −a, que serán las mismas que en a, son: B+ exp(ika) + exp(−ika) = C+ (exp(iqa) + exp(−iqa)) k (B+ exp(ika) − exp(−ika)) = qC+ (exp(iqa) − exp(−iqa)) Estas condiciones nos dan: B+ = exp(−2ika)

34

i kq cot(qa) − 1 i kq cot(qa) + 1

(3.12)

3. Pozo de potencial finito

Capítulo 3. Problemas unidimensionales

Funciones propias impares Ahora consideramos A = −G = A− = 1, B = −F = B− y C = −D = −C− . Haciendo los mismos pasos que en el apartado anterior, obtenemos: B− = exp(−2ika)

i kq tan(qa) − 1 i kq tan(qa) + 1

(3.13)

Solución general Ahora, como dijimos con anterioridad, tenemos que encontrar una combinación lineal que nos dé una solución general con las condiciones que nosotros le impongamos. Como hay una constante global que siempre permanecerá por determinar, podemos considerar arbitrariamente A = 13 . Además imponemos que G = 1, queriendo significar con ello que no consideramos ondas que provengan desde +∞. En definitiva, queremos una solución ψ = Cpar ψpar + Cimpar ψimpar , donde las ψpar , ψimpar vienen dadas, en cada parte del potencial, mediante las ecuaciones generales con los coeficientes vistos en cada uno de los apartados anteriores, y los coeficientes Cpar y Cimpar se obtendrán imponiendo las condiciones anteriormente mencionadas: Cpar A+ + Cimpar A− = A = 1 → Cpar + Cimpar = 1 Cpar G+ + Cimpar G− = G = 0 → Cpar − Cimpar = 0 En definitiva, se obtiene: Cpar = Cimpar =

1 2

Como los coeficientes se obtendrán a partir de B y F , nos interesa expresar estos en función de B− y B+ : B+ + B− B = 2 B+ − B− F = 2 Coeficientes de transmisión y reflexión Como de costumbre, hallamos: 1 (B+ + B− ) (B+ + B− )∗ 4 1 = (B+ − B− ) (B+ − B− )∗ 4

R = |B|2 = T = |F |2 3

Otra alternativa sería calcular los cocientes A/A, B/A, C/A, etc. . . como hicimos con anterioridad. Pero estos dos métodos son totalmente equivalentes

35

4. Barreras de potencial

Capítulo 3. Problemas unidimensionales

Desarrollando los cálculos anteriores, podemos llegar a las siguientes expresiones: (q 2 − k 2 )2 sin2 (2qa) 4k 2 q 2 + (q 2 − k 2 )2 sin2 (2qa) 4q 2 k 2 T = 4k 2 q 2 + (q 2 − k 2 )2 sin2 (2qa)

R =

Es directo comprobar que R + T = 1. Decimos que hay resonancias cuando el coeficiente de transmisión es máximo, esto es, cuando sin2 (2qa) = 0 ⇔ 2qa = nπ, n ∈ N.

4.

Barreras de potencial

Vamos a estudiar ahora las soluciones de la ecuación de Schrödinger para los potenciales conocidos como barreras de potencial, definidos como (3.14), cuya representación esquemática se encuentra en la figura 3.5.   0, si x < −a; V0 , si −a ≤ x ≤ a; V (x) =  0, si x > a.

(3.14)

6

V0 I

II

III

V =0 -

−a

a

x

Figura 3.5: Barrera cuadrada de potencial DESARROLLAR EL EFECTO TÚNEL.

36

5. Pozo de potencial infinito

5.

Capítulo 3. Problemas unidimensionales

Pozo de potencial infinito

Vamos a estudiar este caso particularmente sencillo, un caso en el que, en realidad, lo único que se considera es que la función de onda no se extiende a todo el espacio, sino que se localiza en una región concreta del espacio. Fuera de ella se considera un potencial infinito, es decir, un potencial que la partícula no puede superar nunca, cualquiera que sea su energía.

5.1.

Centrado en el origen

En la figura 3.6 podemos observar como es este potencial. Matemáticamente, viene descrito por las partes:  V (x) =

0 si |x| < a +∞ si |x| > a

(3.15)

Figura 3.6: Pozo de potencial infinito centrado en el origen

V (x) V (x) = +∞ V (x) = +∞ V (x) = 0 x

La forma de potencial nos asegura tajantemente que la función de onda se anulará en los extremos del potencial, es decir, tendremos las condiciones de contorno: ψ(−a) = ψ(a) = 0 ψ(x; |x| > a) = 0

(3.16a) (3.16b)

En la zona central, como el potencial es cero, podremos escribir la ecuación de Schrödinger como ecuación diferencial ordinaria de segundo grado con coeficientes constantes de la manera siguiente: d2 ψ + k2ψ = 0 dx2 37

(3.17a)

5. Pozo de potencial infinito

Capítulo 3. Problemas unidimensionales

Siendo k2 =

2mE ~2

(3.17b)

Las autofunciones de esta ecuación las podemos tomar como: ψ(x) = Ae−ikx + Beikx = C1 eik(x−a) + C2 e−ik(x−a)

(3.18)

El cambio en la segunda igualdad lo podemos hacer sin pérdida de generalidad ninguna sin más que reabsorber la constante e±ika en las distintas constantes A y B. Aplicando las condiciones de contorno dadas por 3.16, obtenemos:
ψ(a) = C1 + C2 = 0 ⇒ C1 = −C2 ⇒ ψ(x) = C1 eik(x−a) − e−ik(x−a) (3.19a) Podemos reescribir la anterior expresión en forma sinusoidal: ψ(x) = 2iC1 sin (k(x − a)) = C 1 sin (k(x − a))

(3.19b)

Aplicando la otra condición de contorno, hallamos: ψ(−a) = C 1 sin(−2ka) = 0 ⇒ 2ka = nπ,

n = 1, 2, . . .

(3.19c)

Esto quiere decir que la energía está cuantizada, dada la relación que existe entre k y E, según 3.17b:

E=

~2 π 2 2 n, 8ma2

n = 1, 2, . . .

(3.20)

Ahora normalizaremos la función de onda para obtener el valor de la constante C 1 :

Z

a 2

2

2

Z

0

|C 1 | sin (k(x − a))dx = |C 1 | −a

= |C 1 |2

sin2 (k(z − a))dz

−2a



z sin(2kz) − 2 4k

0 −2a

1 = |C 1 |2 a = 1 ⇒ |C 1 | = √ a

(3.21)

Por lo tanto, podemos escribir la solución como sigue:

 nπ  1 ψ(x) = C 1 sin(k(x − a)) = √ sin (x − a) 2a a 38

n = 1, 2, 3, . . .

(3.22)

5. Pozo de potencial infinito

Capítulo 3. Problemas unidimensionales

Sin embargo, es más cómodo desarrollar el seno de la diferencia como una diferencia de senos y cosenos, pues la situación se simplifica. En efecto, tenemos:   nπ  nπ 1 nπ  1 (x − a) = √ sin x− ψ(x) = √ sin 2a 2a 2 a a     nπ    nπ 1 1 nπ nπ  x cos − √ cos x sin = √ sin (3.23) 2a 2 2a 2 a a Ahora tenemos en cuenta que: nπ =1 2 nπ n par ⇒ cos =1 2

n impar ⇒ sin

nπ =0 2 nπ sin =0 2

cos

Podemos entonces separar la solución en los dos casos posibles señalados arriba, n par o impar, y obtener, como solución final (ver representación en las figuras 3.7 y 3.8):  1  nπ   x n = 1, 3, 5, . . .  √ cos 2a a  ψn (x) = 1 nπ    √ sin x n = 2, 4, 6, . . . 2a a

(3.24)

Podemos calcular el valor del incremento relativo de energía entre dos niveles consecutivos: En+1 − En = En

~2 π 2 (n 8ma2

+ 1)2 − ~2 π 2 2 n 8ma2

~2 π 2 2 n 8ma2

=

2n − 1 2 =n→∞ 2 n n

(3.25)

Es decir, para números cuánticos suficientemente grandes, el espectro de energías parece casi continuo (ley de los números grandes).

5.2.

Pozo asimétrico

El potencial ahora es de la forma:  V (x) =

0 si 0 < x < 2a +∞ si x < 0 ó x > 2a

(3.26)

Las autofunciones son las mismas que antes, salvo que ahora no conviene tomar como argumentos de las exponenciales complejas los valores x − a, sino simplemente x: ψ(x) = C1 e−ikx + C2 eikx 39

(3.27)

5. Pozo de potencial infinito

Capítulo 3. Problemas unidimensionales

Figura 3.7: Primeros estados separados según su forma: cos kn x o sin kn x

Las condiciones de contorno son ahora: ψ(0) = ψ(L) = 0 ψ(x; 0 < x < 2a) = 0

(3.28a) (3.28b)

Aplicándolas, obtenemos de nuevo C1 = −C2 , y por tanto: ψ(x) = 2iC1 sin kx = C 1 sin kx

(3.29a)

Y también: ψ(2a) = C 1 sin(2ka) = 0 ⇒ 2ka = nπ,

n = 1, 2, . . .

(3.29b)

Evidentemente, la energía también está cuantizada y los posibles valores que obtenemos para la energía son los mismos que en el caso anterior: E=

~2 π 2 2 n, 8ma2

n = 1, 2, . . .

La constante de normalización es también la misma, |C 1 | = 40

(3.29c)

√1 . a

Y obtenemos, como

5. Pozo de potencial infinito

Capítulo 3. Problemas unidimensionales

Figura 3.8: Los cuatro primeros estados

resultado final:  nπ  1 ψn (x) = √ sin x 2a a

n = 1, 2, 3, . . .

(3.30)

Todas las consideraciones hechas respecto al caso anterior son válidas en este caso, con la salvedad de la aparición de los estados “pares” e “impares”. Aquí obtenemos solo un tipo de resultados. En realidad, el resultado tampoco es tan distinto del que obtuvimos en el apartado 5.1. La equivalencia la tenemos que buscar entre esta expresión 3.30 y la ecuación 3.22 y no la 3.24. De hecho, si hacemos un cambio de variable x → x−a volvemos a obtener el caso anterior. Todo esto se debe a que en el primer caso el potencial tiene simetría par, por lo que el operador potencial conmuta con el operador hamiltoniano, y entonces podemos elegir (en ciertos casos y bajo ciertas condiciones) las soluciones como pares o impares. En este segundo caso, el potencial asimétrico, su operador no conmuta con el hamiltoniano, y no podemos hacer esta elección.

41

6. Potenciales delta

6.

Capítulo 3. Problemas unidimensionales

Potenciales delta

En muchas situaciones la forma matemática de un potencial viene dado por una función delta de Dirac, de la forma V (x) = λδ(x − a), diciéndose entonces que existe un potencial delta en a de intensidad λ. Vamos a estudiar en detalle el caso de una barrera de potencial cuadrada (ver sección 4), en el que el valor de la “altura” del potencial tiende a infinito (consideramos E < V0 ), y la anchura del mismo tiende a cero, manteniéndose constante el producto de ambos, es decir:  V0 , |x| ≤ a; V (x) = 0, |x| ≥ a. cumpliéndose: l´ım a→0 V0 → ∞

2aV0 = V

Trataremos de comprobar que en este límite la función de onda en las dos partes es continua, mientras que su derivada es discontinua y satisface: 2m dψ1 dψ2 − = − 2 V ψ1 dx dx ~

(3.31)

Por analogía con la sección 3, las soluciones en las regiones son4 : ψ1 (x) = exp(ikx) +

B+ + B− exp(−ikx) 2

B+ − B− exp(ikx) 2 ψ3 (x) = C exp(qx) + D exp(−qx) ψ2 (x) =

(3.32a) (3.32b) (3.32c)

donde los coeficientes B+ y B− son: B+ = exp(−2ika) B− = exp(−2ika) r q=

4

2m(E + V0 ) ~2

r k=

Para las expresiones de C y D ver 3

42

i kq coth(qa) − 1 i kq coth(qa) + 1 i kq tanh(qa) − 1 i kq tanh(qa) + 1

2mE ~2

(3.32d) (3.32e) (3.32f)

Capítulo 4

Estructura y postulados de la mecánica cuántica Hemos avanzado en la descripción de la perspectiva de la física cuántica. Introducido el concepto de función de onda, y unas reglas para manejarse con ella, hemos resuelto varios problemas relativos a movimiento unidimensional de una partícula. Ahora vamos hacer una pausa para introducir las herramientas matemáticas y la notación necesarias para profundizar en la estructura de la mecánica cuántica. Principalmente, todo consiste en hacer abstracción. Hemos considerado que una partícula venía descrita por una función de onda ψ(x), que representaba como se distribuye la probabilidad de encontrar a ésta en el espacio. También hablamos de la función en el espacio de momentos φ(p), que nos daba lo equivalente para la densidad de probabilidad de su momento. Por último, veíamos que cuando tenemos un pozo infinito, o en los estados ligados de un pozo finito, obtenemos un conjunto de funciones propias. Puede escribirse entonces cualquier estado (ligado en el segundo caso) como combinación lineal de éstas. En este caso, basta una sucesión de números ci para caracterizar el estado por completo. Todas estas maneras de describir un sistema son igualmente válidas, y aquí subyace la motivación de la abstracción. Dada una de éstas caracterizaciones del estado del sistema, siempre es posible en principio encontrar las otras (haciendo la transformada de Fourier, obteniendo funciones propias, etc.). Es por eso que debemos intentar extraer lo que subyace en común a todas ellas al hablar del mismo estado físico.

43

1. Preliminares. Funciones de onda

1.

Capítulo 4. Estructura y postulados

Preliminares. Funciones de onda

Una función de onda básicamente ha de ser continua, derivable, y para guardar consistencia con su interpretación probabilística, ha de verificarse: Z ∞ |ψ(x)|2 d3 x = 1 (caso unidimensional) (4.1) −∞

El conjunto de todas las funciones de Rn en C tal que existe la integral de su cuadrado extendida a todo Rn (esto es, converge a un valor no infinito) se denomina el espacio de funciones de cuadrado integrable L2C (Rn ). En realidad una función de onda estará contenida en el subespacio de L2C (Rn ) de funciones continuas y derivables. Para una partícula en el espacio tridimensional en principio tomaríamos n = 3, ya que la función de onda se evalúa en cada punto de éste, con cordenadas (x, y, z). Sin embargo, continuaremos sin concretar este valor, ya que abre la posibilidad a grados de libertad inferiores (en problemas uni o bidimensionales p.ej.) o superiores (varias partículas). Denotaremos este espacio de funciones de onda L0 . Al ser L0 un espacio vectorial, la combinación lineal de funciones de onda es una función de onda (sin tener en cuenta que la integral anterior valga 1 o no)1 . En este conjunto existe además un producto escalar definido por: Z ∞ Z ∞ Z ∗ ψ ∗ (x)φ(x)dn x (4.2) ... ψ φ= (ψ, φ) = Rn

−∞

−∞

Este producto escalar cumple las propiedades que son de esperar en un espacio euclídeo complejo: (ψ, φ) = (φ, ψ)∗ (Hermiticidad) (4.3) (ψ, αφ + βχ) = α(ψ, φ) + β(ψ, χ) (αψ + βζ, φ) = α∗ (ψ, φ) + β ∗ (ζ, φ) (Linealidad)

(4.4)

(ψ, ψ) = |ψ|2 ≥ 0 (Norma)

(4.5)

(ψ, ψ) = 0 ⇔ ψ ≡ 0

(4.6)

Así, diremos que ψ y φ son ortogonales si (ψ, φ) = 0. A partir de aquí podemos definir el concepto de base ortonormal. Un conjunto de funciones de onda {ui }I (donde el conjunto de índices I en principio es infinito numerable o, de otra forma, N) es una base ortonormal si: 1 Si no vale 1, no sucede nada, salvo que para calcular probabilidades, aparecerán multiplicadas por un factor igual al valor de esta integral. Si van a hacerse otras manipulaciones, lo que procede es normalizarla a 1 dividiendo por una constante real (la raíz cuadrada de su norma, véase más adelante).

44

1. Preliminares. Funciones de onda

Capítulo 4. Estructura y postulados

i) Toda función de onda puede desarrollarse como combinación lineal de las funciones de la base: ∞ X ψ= ci ui (4.7) i=1

ii) El conjunto es ortonormal, esto es : (ui , uj ) = δij

(4.8)

El espacio de funciones de onda se dice que es un espacio de Hilbert separable, y esto conlleva varias consecuencias interesantes. Una de ellas es el hecho de que las bases ortonormales de L0 sean numerables (discretas), como hemos supuesto antes. Veamos la forma en que pueden calcularse los coeficientes ci : (uj , ψ) = (uj ,

∞ X

ci ui ) =

i=1

∞ X

ci (uj , ui ) =

i=1

∞ X

ci δji = cj

i=1

Así: ci = (ui , ψ)

(4.9)

Esto nos lleva inmediatamente a que el producto escalar puede expresarse como: ψ=

∞ X

ai u i

,

φ=

i=1

(ψ, φ) =

∞ X ∞ X

∞ X

bj uj

j=1

a∗i bj (ui , uj ) =

i=1 j=1

∞ X ∞ X

a∗i bj δij =

i=1 j=1

∞ X

a∗i bi

(4.10)

i=1

Y también tenemos una expresión, pues, para la norma al cuadrado de una función: |ψ|2 =

∞ X

|ci |2 =

i=1

∞ X

|(ui , ψ)|2

(4.11)

i=1

A modo formal, se hace necesario introducir el concepto de función de onda generalizada y base generalizada. Como puede verse en el apéndice, mediante la distribución δ de Dirac, puede escribirse: Z ψ(x) = ψ(y) δ(y − x) dn y (4.12) Rn

Tambien, dada la transformada de Fourier φ de ψ, ésta verifica por definición: Z 1 ψ(x) = φ(p)eip·x~ dn p (4.13) n/2 (2π~) Rn 45

1. Preliminares. Funciones de onda

Capítulo 4. Estructura y postulados

Estas expresiones nos inducen a considerar los conjuntos: {δy }y∈Rn , {vp }p∈Rn ,

δy = δ(x − y)

vp =

1 eip·x~ (2π~)n/2

(4.14) (4.15)

como bases del espacio de funciones de onda. El problema es que esto en sí no tiene sentido: las bases de L0 son numerables (éstas son no numerables), y lo que es más importante, están formadas por vectores de L0 (en el caso de las δy ni siquiera tenemos funciones). Sin embargo, aceptaremos estas bases generalizadas porque tienen su utilidad para la comprensión formal de la teoría. Esto es idéntico al proceso de introducción de la distribución δ : tiene un cierto significado y puede implementarse de un modo más riguroso a un nivel mucho más profundo. De hecho esto ya lo hemos utilizado, en la resolución de problemas unidimensionales. Allí obteníamos como funciones propias del hamiltoniano combinaciones de ondas planas que no pertenecían a L0 . Simplemente teníamos los ladrillos con los que se construyen los estados, una base generalizada asociada a la energía. Para estas bases generalizadas, el conjunto de coeficientes es continuo, y la combinación lineal debe ser entendida como integral. Así, reproduciendo los resultados anteriores para bases ortonormales en sentido estricto: Z Base generalizada {ui (x)}I (ut , us ) = u∗t (x)us (x) ds = δ(t − s) (4.16) Rn

Z

Z φ(t)ut (x) dt ,

ψ(x) =

χ(s)us (x) ds

ζ(x) =

(4.17)

I

I

Z φ(t) = (ut , ψ) =

ut (x)ψ(x) dn x

(4.18)

Rn

Z (ψ, ζ) =

φ∗ (t)χ(t) dt

(4.19)

|φ(t)|2 dt

(4.20)

I 2

Z

|ψ| = I

Las dos primeras líneas son la definición definitiva de base ortonormal generalizada: la ortonormalidad y la posibilidad de desarrollar cualquier función de onda como una combinación (integral) de las funciones de la base. I es un conjunto de índices continuo cualquiera, en general, Rm (no tiene por qué ser m = n)2 . El último concepto fundamental para el desarrollo de esta teoría es el de operador. Ya hemos visto que podemos obtener valores medios de magitudes como el momento o En el caso de índices “continuos”, nos referimos ciertamente a un subconjunto I de Rn . Es importante destacar que todos los espacios Rn tienen idéntico cardinal (número de elementos). En principio si la base no es numerable, la consideraremos continua. 2

46

2. Espacio de estados. Notación de Dirac

Capítulo 4. Estructura y postulados

la posición de la partícula, transformando la función de alguna forma (y luego haciendo una integral). En general, trabajaremos exclusivamente con operadores lineales, esto es, aplicaciones lineales en el espacio de funciones de onda3 : A : L0 → L0 A(αψ + βφ) = αAψ + βAφ , ∀α, β ∈ C ψ, φ ∈ L0

(4.21)

Los operadores permiten manipular las funciones de onda para, entre otras cosas, obtener información sobre el estado del sistema. Esto ya se ha sugerido antes, pero avanzaremos más sobre ello. Notemos ahora que el valor medio de un operador A para un estado ψ puede escribirse ahora como: Z ψ ∗ Aψ = (ψ, Aψ) (4.22) hA i = Rn

Llegados a este punto hemos de mirar atrás y darnos cuenta de a dónde hemos llegado. Dado un estado mediante la función de onda ψ, ésta contiene toda la información sobre el sistema. Si especificamos una base {ϕi }N , conocer los coeficientes ci del desarrollo de ψ nos proporciona el mismo conocimiento. Tambien vale esto para una base continua {uα }A . De hecho, dada {δy }y∈Rn , los coeficientes de ψ con respecto a esta base son precisamente ψ(x), con lo que la propia función de onda se representa a sí misma. Los valores de la función de onda del espacio de momentos φ(p) no son sino los coeficientes respecto a la base de ondas planas {vp }p∈Rn . Un mismo estado se ve representado por distintas funciones de onda (según trabajemos en el espacio de momentos, en el espacio natural de la posición,etc.), sucesiones de escalares complejos ci ,. . . Debemos encontrar alguna manera entonces de reunir todas estas interesantes propiedades para acabar con la ambigüedad que encontramos.

2.

Espacio de estados. Notación de Dirac

Las distintas formas en las que un estado queda caracterizado tienen algo en común: todas se refieren a un espacio de Hilbert separable de dimensión infinita. L0 lo es ya, y a él pertenecen las funciónes de onda, funciónes de onda en el espacio de momentos, los coeficientes referentes a una base continua. . . Por otro lado, las sucesiones de escalares ∞ X complejos ci tales que: |ci |2 < ∞ son el marco natural donde se trabaja cuando i=1

tenemos una base discreta, y también forman un espacio de este tipo. 3 Así como diremos operador en vez de aplicación lineal, utilizaremos una notación distinta a la usual f (x) para las funciones: Aψ . Nótese que Aψ(x) es el valor de la función Aψ en x

47

2. Espacio de estados. Notación de Dirac

Capítulo 4. Estructura y postulados

El paso que vamos a dar se apoya sobre un importante hecho: todos estos espacios son perfectamente isomorfos, esto es, existe (al menos) una correspondencia biunívoca entre cada pareja de ellos, que ademas conserva el producto escalar. . . . ↔ ψ(x) ↔ {ci }N ↔ φ(p) ↔ . . . . . . ↔ ψ 0 (y) ↔ {c0j }N ↔ φ0 (q) ↔ . . . Z ... = Rn

............... Z ∞ X ∗ 0 n ∗ 0 ψ (x)ψ (x) d x = ci ci = i=1

φ∗ (p)φ0 (p) dn p = . . .

Rn

Así, simplemente diremos que el estado de un sistema viene representado por un vector de un espacio E, llamado espacio de estados. Este espacio de estados es un espacio de Hilbert separable genérico, sin ninguna de las connotaciones propias de los vistos anteriormente. A los vectores de este espacio se les denomina kets, y se les denota mediante el símbolo | i. Así, usualmente convendremos que el vector de E correspondiente a una función de onda ψ (del espacio natural de la posición) se escribe | ψ i. En este espacio existe también un producto escalar que denotaremos de momento por ( | ψ i, | φ i). El espacio dual E ∗ Un funcional lineal ϕ en un espacio vectorial como E es una función que asocia a cada vector | ψ i un escalar ϕ( | ψ i) ∈ C, de forma que se verifique: ϕ(λ | ψ i + λ0 | ψ i0 ) = λϕ( | ψ i) + λ0 ϕ( | ψ i0 ),

(4.23)

dados λ, λ0 ∈ C, ψ, ψ 0 ∈ E cualesquiera. Al conjunto de todos los funcionales lineales se le denota por E ∗ , y es un espacio vectorial de Hilbert (no exactamente isomorfo a E, ahora lo veremos). El interés de considerar esta definición es descubrir el potencial de la notación de Dirac. Un famoso teorema del análisis funcional4 afirma que a cada funcional lineal ϕ le corresponde un vector | χ i, de forma que ϕ( | ψ i) = ( | χ i, | ψ i) , ∀ | ψ i ∈ E. Denotaremos entonces el funcional lineal asociado al vector | ψ i, que actuando sobre otro vector | φ i, da ( | ψ i, | φ i), por h ψ | . Los vectores del espacio E ∗ los llamaremos entonces vectores bra, y el producto escalar de | ψ i por | φ i se denotará de ahora en adelante como el funcional h ψ | actuando sobre | φ i, escrito h ψ | φ i. Esto justifica los nombres de bra y ket de la notación de Dirac5 . 4 5

El teorema de brackets: del inglés, paréntesis h | i

48

A. Resumen de los postulados

A.

Capítulo 4. Estructura y postulados

Resumen de los postulados

A.1.

Descripción del estado de un sistema

Postulado 1 (Estado de un sistema). En un instante de tiempo concreto t0 , el estado de un sistema físico viene descrito por un vector | ψ(t0 ) i, llamado ket, de un espacio de Hilbert H El hecho de que especifiquemos los estados como vectores de espacios de Hilbert implica la linealidad, lo que en la práctica implica superposición: si dos vectores | ϕ1 i y | ϕ2 i son posibles estados de un sistema, también lo será la combinación lineal | ψ i = α | ϕ1 i + β | ϕ2 i ∀α, β ∈ C.

A.2.

Descripción de magnitudes físicas

Postulado 2 (Magnitudes físicas). Toda magnitud físicamente medible A viene descrita b que actúa en H. Este operador es un observable. por un operador A Con estos dos postulados se observa la radical diferencia entre la mecánica clásica y la cuántica: en la primera el estado de un sistema en un instante de tiempo se describe mediante seis cantidades correspondientes a dos vectores tridimensionales (r y p) y las magnitudes mediante cantidades escalares o vectoriales; en la mecánica cuántica el estado de un sistema está dado por una única función o vector de un espacio vectorial abstracto, y las magnitudes por operadores.

A.3.

Teoría de la medida

Postulado 3 (Posibles valores de una medida). Los únicos posibles resultados de una medida realizada sobre un sistema respecto a una magnitud A es uno de los autovalores b asociado a dicha magnitud. del operador A Los posibles valores de una medida, como es de esperar, son siempre reales, ya que b A es hermítico por definición. b es discreto, entonces la magnitud Si el espectro de los autovalores de un operador A A está cuantizada.

49

A. Resumen de los postulados

A.4.

Capítulo 4. Estructura y postulados

Valores posibles de una medida

Postulado 3.1 (Caso discreto no degenerado). La probabilidad de obtener como resultado de una medida de A sobre un sistema en el estado normalizado | ψ i el autovalor no degenerado λn correspondiente al autovector normalizado a la unidad | ϕn i viene dada por: Pψ (λn ) = |h ϕn | ψ i|2 (4.24) Postulado 3.2 (Caso discreto degenerado). La probabilidad de obtener como resultado de una medida de A sobre un sistema en el estado normalizado | ψ i el autovalor b viene dada por: degenerado λn del operador asociado A gn X i h ϕn | ψ i 2 Pψ (λn ) =

(4.25)

i=1

gn es la degeneración del autovalor λn y los | ϕin i forman una base ortonormal del subespacio de H asociado al autovalor λn . Postulado 3.3 (Caso continuo). La probabilidad de obtener como resultado de una medida de A sobre un sistema en el estado normalizado | ψ i un resultado entre λ y λ + dλ viene dada por: dPψ (λn ) = |h ϕα | ψ i|2 dα (4.26) b asociado al donde | ϕα i es el autovector correspondiente al autovalor α del operador A observable A.

50

Capítulo 5

Espín electrónico. Composición de momentos angulares 1.

Introducción

Hasta ahora hemos considerado las partículas, en particular los electrones, como puntuales y con tres grados de libertad asociados a las posiciones x, y y z, es decir, hemos caracterizado los electrones mediante una función de onda ψ(x, y, z), normalmente. Con estas herramientas de trabajo hemos conseguido muchas descripciones, que ampliaban sin duda los contenidos de la mecánica clásica, explicando fenómenos atómicos con el átomo de hidrógeno, sus niveles de energía, etc. . . ; además, estos modelos sirven de base para modelar otras estructuras de mayor complejidad (eso sí, a través de aproximaciones, pues una resolución analítica de la ecuación de Schr´’odinger suele ser imposible). Pero la realidad, los experimentos, van más allá de lo que hemos estudiado hasta ahora. Cuando se estudian ciertos experimentos en profundidad, aparecen ciertos fenómenos que nuestra teoría, en el estado en que se encuentra hasta nuestro desarrollo actual, no puede explicar. Este resultado no debería ser sorprendente, ya que desde el primer momento, como tuvimos a bien avisar, hemos despreciado toda aportación que la teoría de la Relatividad pudiera hacer: son necesarias sin duda alguna gran cantidad de correcciones relativistas que hemos despreciado: así, la variación de la masa (no la masa en reposo) con la velocidad, efectos magnéticos, etc. . . Estos efectos son pequeños; son pequeños, sí, pero pueden ser medidos. La ecuación de Dirac plantea una descripción cuántica y relativista del electrón, pero sus implicaciones suponen profundas modificaciones respecto a lo que significa la ecuación de Schr´’odinger. Además, de esta teoría emerge con total naturalidad una propiedad fundamental del electrón, y de todas las partículas, el espín, con la misma relevancia que puedan tener la masa en reposo o la carga de la partícula. Sin embargo, al contrario de lo que sucede hoy en día, el espín del electrón fue descubierto experimentalmente mucho antes de que fuera predicho teóricamente. Es más, Pauli desarrolló una teoría en la que el espín se introduce de una forma simple de manera no-relativista, a través de unos 51

2. Evidencia experimental

Capítulo 5. Espín. Momentos Angulares

nuevos postulados que se añaden -en ningún momento sustituyen- a los anteriores. La teoría demuestra estar en perfecto acuerdo con los experimentos. Nosotros estudiaremos aquí esta teoría de Pauli, que no incluye elementos relativistas y es bastante más simple.

2.

Evidencia experimental del espín electrónico

Son los experimentos los que mandan sobre nuestros modelos físicos y no al revés. La evidencia experimental de la existencia de una propiedad sin ninguna analogía clásica, y con comportamiento de momento angular, al que llamamos espín, es numerosa. Nosotros nos fijaremos solo en aquellos experimentos que son más clásicos e intuitivos, y que nos permitirán extraer las consecuencias necesarias para justificar la existencia del espín.

2.1.

Estructura fina de las líneas espectrales

Cuando hablábamos del átomo de hidrógeno afirmábamos sin miedo que la teoría desarrollada entonces explicaba perfectamente las líneas espectrales observadas (las frecuenciasenergías observadas en las transiciones de un estado a otro de los electrones en los distintos niveles energéticos). Bien, esto es cierto, pero solo hasta un grado de precisión determinado. Más allá de esta precisión, estas líneas se desdoblan en otras que tienen frecuencias casi idénticas, pero que se pueden distinguir con los aparatos adecuados. Es decir, existen niveles atómicos que son muy parecidos, pero distintos. Este desdoblamiento, conocido como estructura fina del átomo de hidrógeno, no puede ser explicado con las herramientas y postulados desarrollados hasta ahora. Se hace necesaria la introducción de ciertas hipótesis adicionales.

2.2.

Efecto Zeeman anómalo

Al colocar un átomo en el seno de un campo magnético las lineas espectrales anteriormente mencionadas se subdividen a su vez en un cierto número de líneas equiespaciadas, siendo el intervalo proporcional al campo magnético. Esto se entiende si se tiene en cuenta que aparece un momento magnético asociado al momento angular orbital de la forma: µB q~ M= L µB = ~ 2me El factor µB es conocido como magnetón de Bohr. Esta teoría es confirmada por los experimentos en ciertos casos –efecto Zeeman normal–, pero en otras ocasiones se demuestra incapaz de dar cuenta cuantitativamente de todas estas subdivisiones –efecto Zeeman anómalo–. Se hace necesaria la introducción de ciertas hipótesis adicionales. 52

3. Postulados de Pauli

2.3.

Capítulo 5. Espín. Momentos Angulares

Existencia de momentos angulares semienteros

El experimento de Stern-Gerlach, que no describiremos aquí, revela la existencia de valores semienteros del momento angular. Sin embargo, como veremos más adelante, en este mismo capítulo, la composición de momentos angulares enteros da como resultado otro momento angular de posibles valores distintos, pero siempre enteros. Y los momentos angulares orbitales, como vimos en capítulos anteriores, son siempre enteros. Es decir, existen ciertos momentos angulares, no orbitales, cuyos valores son semienteros. Se hace necesaria la introducción de ciertas hipótesis adicionales.

3.

Postulados de Pauli

Para resolver estas y otras dificultades, Uhlenbeck y Goudsmit propusieron en 1925 la hipótesis de que el electrón gira1 sobre sí mismo y esto le confiere un momento angular intrínseco que llamamos espín. (MOMENTO MAGNÉTICO ??). Posteriormente Pauli propuso estas hipótesis más precisamente, dando una descripición del espín que es el límite del caso relativista, añadiendo una serie de postulados a los ya conocidos. Hasta ahora hemos considerado la cuantización de variables orbitales como p y r, asociando a ellas los observables P y R, que actúan en el espacio de estados Er , que llamaremos espacio de estados orbitales. A estas variables orbitales debemos añadir las variables de espín, que satisfagan los siguientes postulados: 1. El operador de espín S es un momento angular, con lo que –como vimos en el capítulo dedicado al momento angular– esto implica; es decir, satisfacen las siguientes identidades: [Si , Sj ] = i~εijk Sk (5.1) 2. Los operadores de espín actúan sobre un nuevo espacio, el espacio de estados de espín, Es , donde S2 y Sz consituyen un sistema completo de observables compatibles. Es está generado por un conjunto de vectores propios | s, m i comunes a S2 y Sz que satisfacen: S2 | s, m i = s(s + 1)~2 | s, m i Sz | s, m i = m~ | s, m i

(5.2) (5.3)

3. El espacio de estados total, E es el producto tensorial de Er y Es : E = Er ⊗ Es 1

En inglés el verbo girar es spin, de donde viene el nombre spin, castellanizado como espín.

53

(5.4)

4. Espín 1/2

Capítulo 5. Espín. Momentos Angulares

Sabemos que todo valor (que llamábamos j, en general) de un momento angular ha de ser entero o semi–entero. Una partícula está caracterizada por un valor único de s; se dice que la partícula tiene un espín s. El espacio de estados de espín es siempre de dimensión finita de valor 2s + 1, y todos los estados de espín de una partícula son vectores propios de S2 con autovalor s(s + 1)~2 . Concretamente, el electrón es una partícula de espín 1/2, como lo son también el neutrón y el protón. Se conocen partículas con espín desde 0 hasta 11/2, pasando por todos los enteros y semienteros. Con estos nuevos postulados (a los que habremos de añadir los que se refieren a la simetrización o antisimetrización de funciones de onda de partículas idénticas) queda patente que, excepto para el caso particular en que s = 0, no basta con especificar un ket del espacio Er , es decir, los sistemas {X, Y, Z}, o {PX , PY , PZ }, etc. . . no constituyen un sistema completo de observables compatibles. Todo vector de estado de una partícula es una combinación lineal de kets que son productos tensoriales de un ket de Er y otro de Es .

4.

Propiedades del momento angular 1/2

Ahora resumiremos algunas de las propiedades básicas para un momento angular s = 1/2 (no necesariamente un espín), que se deducen directamente de las propiedades generales estudiadas para el momento angular. Tomaremos como base para nuestra descripción una base de Es formada por los vectores { | + i, | − i}, desentendiéndonos por completo de la parte angular. Los vectores | + i y | − i satisfacen las propiedades:

S2 | ± i = | S iz | ± i = h+|−i = h+|+i = b = 1

3 2 ~ |±i 4 1 ± ~| ± i 2 0 h−|−i = 1 | + ih + | + | − ih − |

54

(5.5) (5.6) (5.7) (5.8) (5.9) (5.10)

5. Descripción

Capítulo 5. Espín. Momentos Angulares

Los operadores S± se definen a partir de Sx y Sy , cumpliendo las siguientes propiedades: S± S+ | + i S+ | − i S− | + i S− | − i S+2

= = = = = =

Sx ± iSy 0 ~| + i ~| − i 0 S−2 = 0

(5.11) (5.12) (5.13) (5.14) (5.15) (5.16)

La expresión más general para un vector de estado de espín se describe mediante un vector arbitrario de Es : | χ i = c+ | + i + c− | − i (5.17) Los operadores Sx , Sy u y Sz tienen las siguientes expresiones matriciales en la base { | + i, | − i}: S = ~σ  0 σx = 1  1 σy = 0  0 σz = i

(5.18)  1 0  0 −1  −i 0

(5.19) (5.20) (5.21)

Estas matrices poseen las siguientes propiedades, que pueden deducirse directamente de su expresión explícita anterior: σi2 σi σj + σj σi [σi , σj ] = 2iεijk σi σj Trσi = 0 det σi

5.

= 1 = 0 = iσz = −1

(5.22) (5.23) (5.24) (5.25) (5.26) (5.27)

Descripción no-relativista de una partícula de espín 1/2

Si se han de tener en cuenta todas las propiedades de una partícula (todos sus grados de libertad), el estado de una partícula viene descrito por un vector de E = Er ⊗ Es . Para 55

5. Descripción

Capítulo 5. Espín. Momentos Angulares

ello, habremos de extender (en el sentido matemático de la palabra) los operadores que antes habíamos definido para los dos espacios de estado por separado al espacio total E. Usaremos como sistema de observables el formado por {X, Y, Z, S2 , Sz }. Como base, tomaremos los kets obtenidos como producto tensorial de vectores de las bases de Er , | r i ≡ | x, y, z i, y vectores de la base de Es , | ε i, donde ε puede referirse a | + i o | − i. Es decir, por definición, los vectores de la base anteriormente indicada son de la forma: | r, ε i ≡ | x, y, z, ε i = | r i ⊗ | ε i

(5.28)

Estos vectores, por definición, son autovectores de los operadores X, Y , Z, S 2 y Sz , con autovalores: X | r, ε i = x | r, ε i Y | r, ε i = y | r, ε i Z | r, ε i = z | r, ε i 3 2 S 2 | r, ε i = ~ | r, ε i 4 ~ Sz | r, ε i = ε | r, ε i 2

(5.29) (5.30) (5.31) (5.32) (5.33)

Además, los autovectores cumplen las relaciones de ortonormalidad y cierre: h r0 , ε0 | r, ε i = δε0 ,ε δ(r0 − r) XZ b dr | r, ε ih r, ε | = 1

(5.34a) (5.34b)

ε

Todo vector puede ser expresado en la base { | r, ε i}, de la forma: XZ |ψi = dr | r, ε ih r, ε | ψ i

(5.35)

ε

Los coeficientes h r, ε | ψ i son las coordenadas del vector | ψ i en la base { | r, ε i}, es decir, existen dos funciones que especifican el estado del sistema, que son: ψ+ (r) = h r, + | ψ i ψ− (r) = h r, − | ψ i

(5.36) (5.37)

Estas dos funciones son las dos componentes de una estructura llamada espinor, que escribimos de la forma:   ψ+ (r) [ψ] (r) = (5.38) ψ− (r) La representación análoga a esta del bra h ψ | asociado al estado | ψ i es la matriz adjunta de esta, es decir, un espinor de la forma:
? ? (r) ψ− (r) [ψ]† (r) = ψ+ (5.39) 56

5. Descripción

Capítulo 5. Espín. Momentos Angulares

Esta notación compacta permite simplificar la expresión del producto escalar de dos vectores de estado, que podemos escribir como: Z XZ hψ |ϕi = drh ψ | r, ε ih r, ε | ϕ i = dr [ψ]† (r) [ϕ] (r) (5.40) ε

En la notación anterior, primero se realiza el producto de matrices y después la integración. De los vectores de el espacio de estados, algunos (los de la base, entre otros) serán producto vectorial de la base de vectores correspondiente a las partes orbital y de espín, respectivamente: |ψi = |ϕi ⊗ |χi , (5.41) con | ϕ i y | χ i dados por: h r | ϕ i = ϕ(r) | χ i = c+ | + i + c− | − i El espinor tiene entonces una forma más simple que la general, ya que:   c [ψ] (r) = ψ(r) + c−

(5.42) (5.43)

(5.44)

Esto equivale a decir que las funciones del espinor ϕ± (r) son ϕ± (r) = R simplemente 3 2 2 ∗ c± ϕ(r). Por tanto, h ψ | ψ i = h ϕ | ϕ ih χ | χ i = (|c+ | + |c− | ) V d rϕ (r)ϕ(r). Discusiones similares se pueden hacer para conocer la forma de los operadores que actúan sobre una partícula, ya no solo en el espacio orbital, sino sobre el espacio vectorial de estados, que es producto vectorial de los espacios de espín y orbital, E = Er ⊗ Es , recordando que cada operador de un espacio no afecta sobre las componentes del vector que corresponden al otro espacio. Por ejemplo, el operador Sz , que antes considerábamos que actuaba solo sobre Es , ahora actúa sobre E, pero sobre un vector | r, ± i, actuará solo sobre la parte del espín, no sobre la otra parte del vector. Aunque también puede haber operadores mixtos, como productos de operadores en los que un factor actúe sobre el espacio Er y otro factor sobre el espacio Es . Por ejemplo, el operador Sz en una representación | r, ± i, viene dado por:   1 0 Sz ≡ ~ 0 −1

(5.45)

Que es la misma forma que tenía en Es , evidentemente. Pero operadores como X o Px ya no son simples operaciones, sino matrices cuyos elementos consisten en operaciones:   ∂   0 x 0   (5.46) bXc ≡ bPy c ≡ −i~  ∂y ∂  0 x 0 ∂y 57

6. Momento angular total

Capítulo 5. Espín. Momentos Angulares

Operadores mixtos, más complicados, pueden ser:   ∂ ~2  0 − ∂z  bPz Sy c ≡  ∂ 2 0 ∂z

(5.47)

La representación espinorial también permite expresar el valor esperado de un operador en un estado | ψ i que tenga en cuenta el espín, de la forma: Z hA i = h ψ | A | ψ i = d3 r [ψ]† (r) bAc [ψ]† (r) (5.48) V

6.

Momento angular total en mecánica clásica y en mecánica cuántica

El momento angular total de un sistema de N partículas clásicas respecto a un punto PN O es la suma de los momentos individuales, L = i=1 Li . Cuando las fuerzas externas son cero o paralelas a los vectores de posición, entonces el momento total de las fuerzas es cero y el momento angular total es una (en realidad, tres) constante de movimiento. No lo serán en general, sin embargo, los momentos angulares individuales de cada partícula. Supongamos ahora un sistema cuántico de dos partículas, con un hamiltoniano H 0 = H1 + H2 , siendo Hi : ~2 2 ∇ + V (ri ) (5.49) Hi = − 2mi Es claro entonces que los momentos angulares cumplen la relación de conmutación [Li , H 0 ] = 0, que es el equivalente cuántico de decir que Li “es una constante de movimiento”. Pero supongamos ahora que en el sistema hay también una interacción dada por un potencial central de la forma U (r), con r la distancia entre ambas partículas. Es decir, el hamiltoniano es de la forma H = H 0 + U . Ahora las relaciones de conmutación de los momentos angulares son de la forma [Li , H] = [Li , U ]. Para el caso de la componente z del momento angular de la partícula 1 se tiene:   ~ ∂ ∂ [Liz , U ] = x1 − y1 , (5.50) i ∂y1 ∂x1 que, en general, será distinto de cero. Es decir, debido a U (r), las componentes de los Li dejan de ser constantes de movimiento. Pero podemos intentar ver que sucede con el momento angular total, L = L1 + L2 . Si tomamos la componente z de este operador vectorial, y calculamos su conmutador con H, tenemos:   ~ ∂U ∂U ∂U ∂U [Lz , H] = [Lz , U ] = [L1z , U ] + [L2z , U ] = x1 − y1 + x2 − y2 = i ∂y1 ∂x1 ∂y2 ∂x2 ~ v0 = (x1 (y1 − y2 ) − y1 (x1 − x2 ) − x2 (y1 − y2 ) + y2 (x1 − x2 )) = 0 (5.51) i r 58

7. Suma de espines

Capítulo 5. Espín. Momentos Angulares

Vemos entonces que el momento angular total Li es una constante de movimiento, como sucede en mecánica cuántica. Este resultado lo hemos obtenido para el caso particular en el que sumamos dos momentos angulares orbitales. A menudos se presentarán en mecánica cuántica el problema de sumar dos momentos angulares cualesquiera (sean orbitales o de espín, o combinaciones de ellos). Por ello debemos desarrollar un método que nos permite hacer esto fácilmente. Para ello, comenzaremos con el ejemplo sencillo de la suma de dos espines s = 1/2, en el apartado 7, del que sacaremos algunas directrices, para luego aplicarlas al caso general, en la sección 8.

7.

Suma de espines s = 1/2

Para describir un sistema compuesto de dos espines, tomamos como base el conjunto de vectores obtenidos como producto tensorial de cada uno de los vectores que describen a una partícula de espín s = 1/2, los vectores { | ± i}. La base que decimos será entonces { | ε1 i ⊗ | ε2 i} ≡ { | ε1 i | ε2 i}, donde los εi = ±, y que escribiremos de forma más abreviada como { | ε1 ε2 i}. Estos vectores son autovectores de los operadores Si2 y de los Siz 2 : 3 2 ~ | ε1 ε2 i 4 Siz | ε1 ε2 i = εi ~ | ε1 ε2 i Si2 | ε1 ε2 i =

(5.52) (5.53)

Como queremos expresar estos vectores en otra base, en la que entre el operador S2 , debemos definir el espín total, S = S1 + S2 , comprobando que S es un momento angular, ya que obedece las relaciones de conmutación que definen un momento angular. Tenemos que realizar un cambio de base. Partimos de una base dada por un conjunto completo de observables compatibles, dado mediante {S12 , S22 , S1z , S2z }, y debemos cambiar a otro, dado por {S12 , S22 , S 2 , Sz }, para lo que construiremos una base cuyos vectores denotaremos por { | SM i}, en vista de que queremos obtener: S 2 | SM i = S(S + 1)~2 | SM i Sz | SM i = M ~ | SM i 3 2 Si2 | SM i = ~ | SM i 4

(5.54) (5.55) (5.56)

Estas son las relaciones habituales de vectores y valores propios de un momento angular S. Está claro que los | ε1 ε2 i son autovectores de la componente z del espín total, ya que Sz | ε1 ε2 i = ~2 (ε1 + ε2 ) | ε1 ε2 i. Como εi = ±1, entonces los valores posibles de M son 2

Los índices i hacen referencia a las partículas 1 y 2

59

7. Suma de espines

Capítulo 5. Espín. Momentos Angulares

−1, 0, 1, que es lo que esperaríamos para un momento angular con S = 1 (M = −1, 0, 1) ó S = 0 (M = 0). Para el caso M = 1, la única posibilidad es que ε1 = ε2 = 1, que se corresponde a | + + i (ambos espines “hacia arriba”). Análogamente, para el caso M = −1, el único autovector posible es el | − − i (ambos espines “hacia abajo”). Estos dos autovalores de M son, por tanto, no degenerados. Pero el caso M = 0 es distinto ya que podemos tener ε1 = 1 y ε2 = −1, o bien ε1 = −1 y ε2 = 1. Por lo tanto, M = 0 es doblemente degenerado. La matriz que representa el operador  1 ~ 0 Sz ≡   2 0 0

Sz es: 0 0 0 0

 0 0 0 0  0 0 0 −1 | ε1 ε2 i

(5.57)

El subíndice hace referencia a la base en la que se escribe la matriz. El operador S2 se puede obtener de varias formas, ya que: S2 = S21 + S22 + 2S1 S2 S1 S2 = S1x S2x + S1y S2y + S1z S2z = S1z S2z +

1 (S1+ S2− + S1− S2+ ) 2

La matriz que representa entonces al operador S2 es:   2 0 0 0 ~ 0 1 1 0  S2 ≡  2 0 1 1 0 0 0 0 2 | ε1 ε2 i

(5.58) (5.59)

(5.60)

Sz ya es diagonal en esta base. Por los teoremas que vimos en los primeros capítulos, ahora solo resta diagonalizar la restricción de S2 al subespacio generado por { | + − i, | − + i} (el subespacio correspondiente a M = 0), es decir, la parte central de S2 , ya que es la que debería romper la degeneración de Sz y no lo hace. En definitiva, hay que diagonalizar la submatriz:   2 1 1 2 S0 ≡ ~ (5.61) 1 1 Si designamos con λ~2 , los autovalores, se tiene que λ = 2 o bien λ = 0. Los correspondientes autovectores son: 1 λ = 2 −→ √ [ | + − i + | − + i] 2 1 λ = 0 −→ √ [ | + − i − | − + i] 2 60

(5.62) (5.63)

7. Suma de espines

Capítulo 5. Espín. Momentos Angulares

Encontramos entonces, en definitiva, que, en la base del espín total, el operador S2 tiene dos autovalores posibles, 2~2 y 0. Como S 2 | SM i = S(S + 1)~2 | SM i, tenemos que los valores posibles de S son, respectivamente, S = 1 y S = 0. Los autovectores asociados a S = 1 son: | + +i

1 √ ( | + − i + | − + i) 2

| − −i

A su vez, éstos son autovectores de Sz , con autovalores: | + + i −→ M = 1 1 √ ( | + − i + | − + i) −→ M = 0 2 | − − i −→ M = −1 Por otro lado, asociado a S = 0, solo existe un autovector, autovalor de Sz dado por M = 0.

√1 2

( | + − i + | − + i), con

Tenemos, tal y como cabría esperar, que, para cada S, los autovalores posibles de Sz son de la forma −S, −S + 1, . . . 0 . . . , S − 1, S. Por otro lado, hemos visto que los vectores que obtenemos como base son de la forma: singlete

triplete

1 | 00 i = √ ( | + − i − | − + i) 2   | 11 i = | + + i | 10 i = √12 ( | + − i + | − + i)  | 11 i = | − − i

Por último, destacar algunas propiedades de simetría de la base anterior. El vector | 00 i es lo que se conoce como singlete (dado que aparece solo), y es antisimétrico respecto del intercambio de dos partículas. El conjunto de vectores | 1M i, al estar compuesto por tres elementos, es lo que se conoce como triplete, y es simétrico respecto del intercambio de dos espines. Estas propiedades de simetría, que podrían considerarse fruto de la casualidad, serán de gran importancia al tratar partículas idénticas, y aparecerán con frecuencia en la mecánica cuántica. Además, ofrece una regla para recordar la forma de los vectores de la base. De los vectores | 1M i han de formar parte, claramente, los vectores | + + i (con M = 1, pues los espines están “hacia arriba”, y | − − i (con M = −1, ya que los espines se encuentran “hacia abajo”). El otro vector vector que falta para completar el triplete ha de tener M = 0 y por tanto ser una combinación lineal de los vectores | + − i y | − + i (que, además, son los que faltan aparte de los | + + i y | − − i). Pero como ha de ser simétrico, como el resto del triplete, y estar normalizado, entonces está claro que ha de ser | 10 i = √12 ( | + − i + | − + i). Finalmente, el estado singlete, con S = M = 0, ha de ser también una combinación lineal de los | + − i y | − + i. Pero como ha de ser antisimétrico, entonces | 00 i = √12 ( | + − i − | − + i). 61

8. Método general

8. 8.1.

Capítulo 5. Espín. Momentos Angulares

Método general de suma de momentos angulares Introducción

Sea un momento angular J y una base { | k, j, m i} asociada a J. Se cumplen entonces las siguientes ecuaciones de valores propios: J2 | k, j, m i = j(j + 1)~2 | k, j, m i Jz | k, j, m i = m~ | k, j, m i Además, se tienen los operadores J± , cuya acción viene dada por: p J± | k, j, m i = ~ j(j + 1) − m(m ± 1) | k, j, m ± 1 i

(5.64) (5.65)

(5.66)

Llamemos ε(k, j) a los subespacios generados por los vectores de la base correspondiente a k y j fijos. Estos subespacios: Son (2j + 1)-dimensionales. Son invariantes bajo la acción de J2 , Jz y F (J en general, es decir, los elementos de matriz no nulos de estos operadores se encuentran solo en cada subespacio ε(k, j). Sus elementos de matriz, en cada subespacio ε(k, j), son independientes de k. Este pequeño resumen anterior contiene la información básica que necesitaremos de la teoría del momento angular de un sistema. Ahora consideremos dos subsistemas distintos (que indicaremos a partir de ahora como 1 y 2). El espacio de estados del sistema conjunto es el producto tensorial de los espacios de estados de los sistemas por separado, ε = ε1 ⊗ ε2 . Tomamos como base de estados del sistema los vectores de la forma | k1 , k2 ; j1 , j2 ; m1 , m2 i, con: | k1 , k2 ; j1 , j2 ; m1 , m2 i = | k1 , j1 , m1 i ⊗ | k2 , j2 , m2 i El espacio de estados se puede descomponer en la suma directa: X X ε= ε(k1 , k2 ; j1 , j2 ) = ε1 (k1 ; j1 ) ⊗ ε(k2 ; j2 ) ⊕

(5.67)

(5.68)

⊕

Los subespacios ε(k1 , k2 ; j1 , j2 ) son (2j1 + 1)(2j2 + 1)-dimensionales. El momento angular total del sistema, J, es la suma de los momentos angulares3 J1 y J2 , es decir: J = J1 + J2 (5.69) 3

En realidad, es la suma de las extensiones de los operadores J1 y J2 , desde los espacios originales en que actuaban, ε1 y ε2 , respectivamente, al espacio de estados ε.

62

8. Método general

Capítulo 5. Espín. Momentos Angulares

Se puede comprobar, de la definición anterior, que se cumplen las siguientes relaciones de conmutación:         Jz , J21 = Jz , J22 = J2 , J21 = J2 , J22 = [J1z , Jz ] = [J2z , Jz ] = [J1 , J2 ] = 0 (5.70) El operador J2 se puede expresar como: J2 = J21 + J22 + 2J1z J2z + J1+ J2− + J1− J2+

8.2.

(5.71)

Cambio de base

Los vectores { | k1 , k2 ; j1 , j2 ; m1 , m2 i} son autovectores simultáneamente de los operadores: J21 , J22 , J1z y J2z , (5.72) con sus respectivos autovalores. También los operadores J21 ,

J22 ,

J2

y Jz

(5.73)

conmutan entre sí; nosotros buscamos una base basada en los operadores 5.73 para el estudio del momento angular total, por lo que habrá que realizar un cambio de base de 5.72 a 5.73. Es decir, buscamos vectores, que denotaremos (para dos valores j1 y j2 fijos) como | J, M i que satisfagan simultáneamente: J2 | J, M i = J(J + 1)~2 | J, M i Jz | J, M i = M ~ | J, M i

(5.74) (5.75)

Además, por supuesto, deberán ser autovectores de J21 y J22 , con autovalores asociados a j1 y a j2 , respectivamente. Como los subespacios ε(k1 , k2 ; j1 , j2 ) tienen elementos de matriz no nulos solo dentro de cada subespacio, el problema consistirá “solo” en realizar un cambio de base dentro de cada subespacio. Además, como los elementos de matriz de cualquier operador función de J21 , J22 , J 2 yJz son independientes de los valores de k1 y k2 , el problema de la diagonalización será el mismo en cada subespacio. Por tanto, adoptaremos la siguiente notación reducida: ε(k1 , k2 ; j1 , j2 ) −→ ε(j1 , j2 ) | k1 , k2 ; j1 , j2 ; m1 , m2 i −→ | j1 , j2 ; m1 , m2 i Por otro lado, como J es un momento angular, y los ε(j1 , j2 ) son invariantes bajo la acción de cualquier función de J, se tiene que ε(j1 , j2 ) se puede descomponer en suma directa de los distintos subespacios asociados a J (que tendrán que ver con los autovalores de este operador), llamémoslos ε(k, J), que son invariantes bajo la acción de J2 , Jz y J± , esto es: X ε(j1 , j2 ) = ε(k, J) (5.76) ⊕

Tenemos, por ahora, un doble objetivo: 63

8. Método general

Capítulo 5. Espín. Momentos Angulares

Dados j1 y j2 , encontrar los valores posibles de J, y cuantos subespacios ε(k, J) existen para cada par de valores j1 y j2 . Expresar los autovectores de J2 y Jz en función de los vectores de la “antigua” base, es decir, 5.72.

8.3.

Autovalores de Jz

Asumimos, sin pérdida de generalidad alguna, que j1 ≥ j2 . Los vectores | j1 , j2 ; m1 , m2 i son autovectores de Jz , luego Jz | j1 , j2 ; m1 , m2 i = (m1 + m2 )~ | j1 , j2 ; m1 , m2 i,

(5.77)

y sus autovalores son M = m1 + m2 . Como los valores máximos que pueden alcanzar m1 y m2 son j1 y j2 , respectivamente, y los mínimos −j1 y −j2 , entonces los valores posibles de M son: M = j1 + j2 ,

j1 + j2 − 1,

...

,

−j1 − j2 + 1,

−(j1 + j2 )

(5.78)

Ahora estudiaremos la degeneración4 gj1 ,j2 (M ) de cada autovalor M , es decir, el número de vectores que hay con el mismo autovalor. Para ello consideremos el gráfico 5.1, en el que se muestran los valores posibles de m1 y m2 y el valor M al que dan lugar, para dos valores concretos de j1 y j2 (j1 = 3 y j2 = 2). Con cada punto va asociado un vector | j1 , j2 ; m1 , m2 i. Todos los puntos situados sobre la misma línea dan origen a un mismo valor de M (como se puede comprobar matemáticamente de manera fácil), luego el número de tales puntos es el número de vectores asociados al mismo autovalor M , es decir, la degeneración de dicho autovalor gj1 ,j2 (M ). Así, gj1 ,j2 (j1 +j2 ) = 1, gj1 ,j2 (j1 + j2 − 1) = 2, . . . , gj1 ,j2 (j1 − j2 ) = 2j2 − 1. Conforme M disminuye en una unidad, su degeneración aumenta (también en una unidad), hasta llegar a j1 − j2 . A partir de ahí toma un valor constante hasta llegar a M = −(j1 − j2 ), valor tras el cual la degeneración empieza a disminuir de forma simétrica con respecto a como subía, es decir, se tiene gj1 ,j2 (−M ) = gj1 ,j2 (M ), de donde, concretando, gj1 ,j2 (−(j1 + j2 )) = 1. Esta información se resume en el gráfico 5.2. Por último, resaltar que los valores de M serán enteros si j1 y j2 son a la vez enteros o semienteros, mientras que si uno de los dos es entero y el otro semientero, entonces los valores de M serán semienteros.

8.4.

Autovalores de J2

Comentamos en el apartado anterior que el máximo valor que puede tomar M es M = j1 + j2 . Por otro lado, como los valores de M de un momento angular van desde 4

Explicarlo más matemáticamente

64

8. Método general

Capítulo 5. Espín. Momentos Angulares m2 6

x

x

x j2 = 2

x

x

(3, 2)

x

x

M = −1

M =5

x

x

x

x

x

x

x

M =4 M = −3

x

x

x

x

x

x

x - m1

j1 = 3 x

x

x

x

x

x

x

x

x

x

x

x

x

x

(−3, −2)

M =0 Figura 5.1: Valores de M y autovectores asociados.

gj1 ,j2 (M )

5 4 3 2 1 M −5

−4

−3

−2

−1

0

1

2

3

Figura 5.2: Valores de M y degeneración asociada. 65

4

5

8. Método general

Capítulo 5. Espín. Momentos Angulares

− hasta J, deducimos que no hay valores de J mayores que j1 + j2 . Con J = j1 + j2 se asocia un subespacio único, ε(J = j1 + j2 ), en el que está el único vector con M = j1 + j2 . Este subespacio también contiene un único vector (porque J es un momento angular) con autovalor de Jz asociado a M = j1 + j2 − 1. Ahora bien, este autovalor es doblemente degenerado, luego existe otro vector en ε(j1 , j2 ). Dicho vector está en otro subespacio, ε(J = j1 + j2 − 1), y además, está en él de forma no degenerada. Formalicemos esto un poco más. Sea pj1 ,j2 (J) el número de subespacios de un autovalor asociado a J. Entonces se tiene: gj1 ,j2 (M ) = pj1 ,j2 (|M |) + pj1 ,j2 (|M | + 1) + . . . + pj1 ,j2 (j1 + j2 ) (5.79) Podemos invertir la relación para obtener pj1 ,j2 (J): pj1 ,j2 (J) = gj1 ,j2 (M = J) − gj1 ,j2 (M = J + 1)

(5.80)

De aquí podemos observar que: pj1 ,j2 (J > j1 + j2 ) = 0, ya que gj1 ,j2 (M ) = 0 para |M | > j1 + j2 . pj1 ,j2 (j1 − j2 ≤ J ≤ j1 + j2 ) = 1, porque gj1 ,j2 (M ) disminuye en una unidad al aumentar el valor de M en una unidad. pj1 ,j2 (J < j1 − j2 ) = 0, porque no podemos aceptar valores de pj1 ,j2 (J) negativos, no tendría sentido. Es decir, hemos llegado a la conclusión de que: J = j1 + j2 ,

j1 + j2 − 1,

...

,

|j1 − j2 |

(5.81)

Hemos puesto |j1 − j2 | para generalizar al caso en que fuera j1 < j2 , aunque no tiene mayor trascendencia, claro está. Otra forma de ver que el valor mínimo de J ha de ser j1 − j2 es ver la dimensión de ε(j1 , j2 ). Llamemos x a ese valor mínimo de J. Para cada valor de J hay un subespacio de dimensión 2J + 1 debida a los distintos valores posibles de M . También sabemos que el valor máximo de J es J = j1 + j2 . Como la dimensión de ε(j1 , j2 ) es (2j1 + 1)(2j2 + 1), se debe cumplir: j1 +j2 X 2J + 1 = (2j1 + 1)(2j2 + 1) (5.82) J=x

Desarrollando la suma, se encuentra que x = |j1 − j2 |.

8.5.

Autovectores de J2 y Jz

66

Formulario general Mecánica ondulatoria

Z

∞

Z ∞ 1 ψ(p)eipx/~ dp ψ(x) = √ 2π~ −∞ Z ∞ 1 ψ(x)e−ipx/~ dx ψ(p) = √ 2π~ −∞ √ x2 e− 2σ2 dx = 2π σ

−∞

1 2π

Z

∞

0

e−ik(x−x ) dk = δ(x − x0 )

−∞

i~ j(x) = 2m



∂ψ ∗ ∂ψ ψ − ψ∗ ∂x ∂x



ρ = |ψ|2 ∂ρ + divj(x) = 0 ∂t dhA i i = h[H, A] i dt ~  Potencial delta:

dψ dx



 − 0+

67

dψ dx

 = 0−

2m αψ(0) ~2

i~∂t | ψ i = H | ψ i H | ϕin i = En | ϕin i X | ψ(t) i = cn e−iEn t/~ | ϕn (0) i n

~ ∂ i ∂x P2 H = + V (x) (Expresión general) 2m Z ∞ p2 1 dp0 V (p)ψ(p − p0 ) Eψ(p) = ψ(p) + √ 2m 2π~ −∞ P =

Problemas unidimensionales Pozo infinito Centrado en origen, anchura 2a

~2 π 2 2 n 8ma2   nπ  1  x , n = 2, 4, 6, . . .  √ sin 2a  a  ψn (x) = 1 nπ   √ cos x , n = 1, 3, 5, . . . 2a a En =

Centrado en L/2, anchura L

~2 π 2 2 En = n 2mL2 r  nπ  2 ψn (x) = sin x L L

68

Oscilador armónico H=

1 P2 + mωX 2 2m 2

ˇ = X

r

mω X ~

r

1 P mω~

Pˇ =

ˇ = 1 (Pˇ 2 + X ˇ 2) H 2 ˇ H = ~ω H

1 ˇ a = √ (X + iPˇ ) 2 1 ˇ a† = √ (X − iPˇ ) 2 ˇ = a† a + 1 = N + 1 H 2 2 √ n + 1 | ϕn+1 i a† | ϕn i = √ a | ϕn i = n | ϕn−1 i 1 | ϕn i = √ (a† )n | ϕ0 i n! √ h ϕn0 | a | ϕn i = nδ 0 √ n ,n−1 † h ϕn0 | a | ϕn i = n + 1 δn0 ,n+1 r i √ ~ h√ h ϕn0 | X | ϕn i = n + 1 δn0 ,n+1 + n δn0 ,n−1 r2mω h i √ mω~ √ h ϕn0 | P | ϕn i = i n + 1 δn0 ,n+1 − n δn0 ,n−1 2

69

 n 1/2   n mω 2 ~ mω 1/4 mω d 1 x− e− 2~ x n 2 n! mω π~ ~ dx  mω 1/4 mω 2 · e− 2~ x π~   1/4 mω 2 4 mω 3 xe− 2~ x π ~ i mω 2  mω 1/4 h mω 2 x2 − 1 e− 2~ x 4π~ ~ 

ϕn (x) = ϕ0 (x) = ϕ1 (x) = ϕ2 (x) =

Momento angular General

J [J , J ]  i2 j J ,J [Jz , J] [J+ , J− ] JJ∓ J2

J 2 | k, j, m i Jz | k, j, m i 2 J J | k, j, m i Jz J | k, j, m i

Jx ± iJy i~εi,j,k Jk 0 ±~J 2~Jz Jx2 + Jy2 ± ~Jz = J 2 − Jz2 ± ~Jz 1 = (J+ J− + J− J+ ) + Jz2 2

= = = = = =

j(j + 1)~2 | k, j, m i m~ | k, j, m i j(j + 1)~2 | k, j, m i (m ± 1)~ | k, j, m i p J | k, j, m i = ~ j(j + 1) − m(m ± 1) | k, j, m ± 1 i = = = =

h k, j, m | Jz | k 0 , j 0 , m0 i = m~ δk,k0 δj,j 0 δm,m0 h k, j, m | J 2 | k 0 , j 0 , m0 i = j(j + 1)~2 δk,k0 δj,j 0 δm,m0 p h k, j, m | J | k 0 , j 0 , m0 i = ~ j(j + 1) − m0 (m0 ± 1) δk,k0 δj,j 0 δm,m0 ±1 70

~ h k, j, m | Jx | k 0 , j 0 , m0 i = δk,k0 δj,j 0 2 hp i p 0 0 0 0 j(j + 1) − m (m + 1) δm,m0 +1 + j(j + 1) − m (m − 1) δm,m0 −1 ~ h k, j, m | Jy | k 0 , j 0 , m0 i = δk,k0 δj,j 0 2i hp i p j(j + 1) − m0 (m0 + 1) δm,m0 +1 − j(j + 1) − m0 (m0 − 1) δm,m0 −1

(1/2)

(Jz )

~ = 2

(1/2)

(Jx )



~ = 2

 1 0 0 −1



2 (1/2)

(J )

 0 1 1 0

(1/2)

(Jy )



(Jz )(1)

(Jx )(1)

 1 0 0 = ~ 0 0 0  0 0 −1

  0 1 0 ~ = √ 1 0 1 2 0 1 0

3~2 = 4 ~ = 2





 1 0 0 1

0 −i i 0





(J 2 )(1)

(Jy )(1)

 1 0 0 = 2~2 0 1 0 0 0 1

  0 −i 0 ~ = √  i 0 −i 2 0 i 0

Momento angular orbital

 cos ϕ ∂ ∂ i~ sin ϕ + ∂θ tan θ ∂ϕ   ∂ sin ϕ ∂ i~ − cos ϕ + ∂θ tan θ ∂ϕ ∂ −i~ ∂ϕ  2  ∂ 1 ∂ 1 ∂2 2 −~ + + ∂θ2 tan θ ∂θ sin2 θ ∂ϕ2 

Lx = Ly = Lz = L2 =

L2 Ylm (θ, ϕ) = l(l + 1)~2 Ylm (θ, ϕ) Lz Ylm (θ, ϕ) = m~Ylm (θ, ϕ) 71

  ∂ ∂ L = ~e ± + i cot θ ∂θ ∂ϕ l l L+ Yl (θ, ϕ) = 0 ⇒ Yl (θ, ϕ) = cl (sin θ)l eilϕ r (−1)l (2l + 1)! cl = 2l l! 4π ±iϕ

1 Y00 (θ, ϕ) = √ 4π r

3 Y1±1 (θ, ϕ) = ∓ sin θe±iϕ 8π r 3 cos θ Y10 (θ, ϕ) = ∓ 4π r 15 Y2±2 (θ, ϕ) = sin2 θe±2iϕ 32π r 15 ±1 Y2 (θ, ϕ) = ∓ sin θ cos θe±iϕ 8π r
5 0 3 cos2 θ − 1 Y2 (θ, ϕ) = ∓ 16π

Composición de espines



0  ~ 0 S1x =  2 1 0  0 0 ~ 0 0 S1y =   2 i 0 0 i  1 0 ~ 0 1 S1z =  2 0 0 0 0

0 0 0 1 −i 0 0 0 0 0 −1 0

1 0 0 0

 0 1  0 0  0 −i  0 0  0 0  0 −1

 0  ~ 1 S2x =  2 0 0  0 ~ i S2y =   2 0 0  1 ~ 0 S2z =  2 0 0

72

1 0 0 0

0 0 0 1

 0 0  1 0

 0 0 0 0  0 −i i 0  0 0 0 −1 0 0   0 1 0 0 0 −1 −i 0 0 0

s P21 | SM i = (−1)S+1 | SM i
1 2 S − S12 − S22 S1 S2 = 2

Teoría de perturbaciones General

ˆ H = H 0 + W = H 0 + λW | ϕ n i = | 0 i + λ | 1 i + λ2 | 2 i + · · · En (λ) = En0 + λε1 + λ2 ε2 + · · · ...

Nivel no degenerado

ˆ | ϕn i ε1 = h ϕ n | W ˆ | ϕn i X X h ϕip | W |1i = | ϕip i 0 − E0 E n p p6=n i 2 ˆ | ϕi i X X h ϕip | W p ε2 = 0 0 En − Ep p6=n i

Nivel degenerado ˆ n | 0 i = ε1 | 0 i W

73