Cuantica 1 y 2

MECANICA CUANTICA I-II
Marzo 2004
Universidad de Concepcion
Departamento de Fsica
Universidad de Concepcion
Concepcion
Prefacio
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
Concepcion, Marzo 2004 LOS AUTORES
Indice general
Parte I Mecanica Cuantica I
Parte II Mecanica Cuantica II
1. Reglas de seleccion
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.1. Introduccion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.2. Rotaciones espaciales y transformaciones unitarias . . . . . . . . . . . . . . . 5
1.3. Operadores escalares . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.4. Operadores vectoriales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.5. Operadores Tensoriales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.6. Tarea 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2. Teora de Perturbaciones
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.1. Introduccion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.2. V independiente del tiempo. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.2.1.Caso de espectro no degenerado . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.2.2.Caso de espectro degenerado . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
2.3. Teora de perturbaciones dependiente del tiempo . . . . . . . . . . . . . . . . 26
2.3.1.Cuadro de Interaccion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
VIII

Indice general
2.3.2.Oscilaciones de Rabi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
2.3.3.Series de Dyson . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
2.3.4.Determinacion de probabilidades de transicion . . . . . . . . . . . . . . 31
2.4. Metodo variacional . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
3. Partculas identicas
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
3.1. Introduccion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
3.2. Simetra de permutacion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
3.3. Postulado de simetrizacion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
3.4. Sistema de dos electrones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
3.4.1.Estados de un sistema de tres electrones . . . . . . . . . . . . . . . . . . . 52
3.5. Atomo de Helio . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
3.5.1.Teora de perturbaciones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
3.5.2.Metodo variacional . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
3.5.3.Metodo de Hyllerras (1929) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
4. Cuantizacion del campo electromagnetico
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
4.1. Introduccion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
4.2. Ecuacion de ondas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
4.3. Campo electromagnetico libre . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
4.4. Estados de fotones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
4.5. Fotones como excitaciones mecanico cuanticas del CEM . . . . . . . . . . 70
4.6. Dependencia temporal de los operadores de creacion y destruccion . 72
4.7. Estados coherentes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
4.8. Operador densidad . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
Indice general IX
5. Interaccion de atomos con campos cuantizados
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
5.1. Introduccion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
5.2. Aproximacion Dipolar . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
5.3.

Atomo y campo cuantizado . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
5.4. Aproximacion de onda rotante . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
5.5.

Atomo de dos niveles con un unico modo del CEM . . . . . . . . . . . . . . 77
5.6. Colapos y resurguimientos de la poblacion atomica . . . . . . . . . . . . . . 79
5.7. Decaimiento atomico espontaneo en espacio libre . . . . . . . . . . . . . . . . 81
5.8. Wigner-Weisskopf teora del decaimiento espontaneo . . . . . . . . . . . . . 82
5.9. Coherencia atomica e interferencia . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
6. Teora de scattering
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
6.1. Introduccion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
6.2. Determinacion de la seccion transversal de dispersion . . . . . . . . . . . . 93
6.3. Ecuacion integral de dispersion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 94
6.4. Potenciales centrales V (r) = V (r). Ondas parciales. . . . . . . . . . . . . . . 96
Parte I
Mecanica Cuantica I
Parte II
Mecanica Cuantica II
1. Reglas de seleccion
Resumen En este captulo se estudia....
1.1. Introduccion
1.2. Rotaciones espaciales y transformaciones unitarias
Consideremos la rotacion un punto del espacio en un angulo . Esta
se describe por:
r
= exp()r
= r + r +
1
2!
r +
1
3!
r +. . .
Esto corresponde a una transformacion activa en que un punto P de coorden-
das r se desplaza hasta un punto P
de coordenadas r
. Antes de continuar
veremos que le sucede a la funcion de onda (r) ante una rotacion innite-
simal d (r
= r +d r)
(r
) = (r) +dr (r) = (r) +d r (r)

= (r) +d (r ) (r) = (r) +
i
d L(r)
(r
) = exp
_
i
L
_
(r).
Existe la posibilidad de efectuar una transformacion pasiva, en que el mismo
punto es denotado por los vectores r y r
debido a la rotacion del sistema

de referencia en un angulo , es decir, se describe el mismo punto en dos
sistemas de referencia diferentes r
= exp()r.
(r
) = (r)
=
(exp()r) = exp
_
L
_
(r) = (r)

(r) = exp
_
i
L
_
(r),
6
es decir,
(r) corresponde a la funcion de onda transformada unitaria en

el espacio de estados. As concluimos que: Una rotacion del sistema de
coordenadas en el espacio fsico real de un partcula produce un cambio de base
en el espacio de estados descrita por la transformacion U = exp
_
i
L
_
.
1.3. Operadores escalares
Los diversos observables asociados a un sistema fsico se pueden cla-
sicar en terminos de sus propiedades ante rotaciones. Esta clasicacion es
de gran utilidad al momento de determinar los elementos de matriz de los
operadores y a las reglas de seleccion que determinan cuales elementos son
no nulos. Un operador arbitrario bajo rotaciones transforma de acuerdo a
= UU
, U = exp(
i
J),
donde J es el momentum angular total.
Si los operadores son invariantes ante la rotacion se dice que son ope-
radores escalares. Un ejemplo es el Hamiltoniano de un sistema aislado. Estos
operadore satisfacen la relacion
USU
= S US SU = 0 = [S, U] [S, J] = 0.
Esto ultimo se demuestra al considerar innitesimal, por lo tanto
[S, J
z
] = 0, [S, J
2
] = 0.
El escribir estas expresiones en la base [j, m) se obtiene:
j
, m
[ SJ
z
J
z
S [j, m) = (mm
) j
, m
[ S [j, m) = 0
j
, m
[ SJ
2
J
2
S [j, m) =
2
(j(j + 1) j
(j
+ 1)) j
, m
[ S [j, m) = 0.
As
j
, m
[ S [j, m) =
_
,= 0 si m
= m, j
= j
0 en otros casos.
Esto nos ha permitido obtener las siguientes reglas de seleccion para opera-
dores escalares:
m = 0, j = 0.
Sin embargo, la invarianza ante rotaciones a un nos permite obtener mayor in-
formacion de los elementos de matriz del operador, ya que si s es un autovalor
de S en el estado [j, m) se cumple que:
S [j, m) = s [j, m) ,
SJ
+
[j, m) = J
+
S [j, m) = sJ
+
[j, m) ,
1. Reglas de seleccion 7
ya que [J
, S] = 0 y el valor propio s no depende se m, solo puede depender

de j, es decir, s = s(j) = j[ S [j). En resumen, podemos escribir
j
, m
[ S [j, m) =
j,j

m,m
j[ S [j) .
Si el sistema aislado depende de otros n umeros cuanticos que resumimos en
el ndice se tiene
j
, m
[ S [j, m; ) =
j,j

m,m
j;
[ S [j; ) .
1.4. Operadores vectoriales
Ciertamente, los operadores escalares es una clase restringida de ope-
radores y, en general, existe una transformacion del operador ante una rota-
cion. Por ejemplo, estan los operadores que poseen tres componentes, como
el momentum, que se pueden escribir como V [V
x
, V
y
, V
z
], que transforman
como las coordenadas de un punto en el espacio, es decir,
exp
_
i
J
_
V exp
_
J
_
= exp() V ,
donde representa el angulo de rotacion. En el caso innitesimal
exp
_
i
d J
_
V exp
_
d J
_
= V +d V
_
1 +
i
d J
_
V
_
1
i
d J
_
= V +d V
id J V iV d J = d V .
Escribiendo las relaciones para cada una de las componentes de J se obtiene
[J
x
, V
x
] = 0, [J
y
, V
x
] = iV
z
[J
z
, V
x
] = iV
y
[J
x
, V
y
] = iV
z
, [J
y
, V
y
] = 0 [J
z
, V
y
] = iV
x
[J
x
, V
z
] = iV
y
, [J
y
, V
z
] = iV
x
[J
z
, V
z
] = 0
[J
j
, V
k
] = i
jkl
V
l
,
donde
jkl
es el simbolo de Levi-Civita. En el caso de una rotacion nita se
puede obtener la transformacion del operador empleando los conmutadores
previos y la formula Baeker-Hausdorf-Campbell. Reescribiremos las expresio-
nes previas empleando los operadores de subida y bajada dados, respectiva-
mente, por V
= V
x
iV
y
, obteniendose:
[J
+
, V
+
] = 0, [J
, V
+
] = 2V
z
[J
z
, V
+
] = V
+
[J
+
, V
] = 2V
z
, [J
, V
] = 0 [J
z
, V
] = V
[J
+
, V
z
] = V
+
, [J
, V
z
] = V
[J
z
, V
z
] = 0
8
Empleando estos conmutadores en la base [j, m) podemos obtener infor-
macion de los elementos de matriz V . Por ejemplo, las reglas de seleccion
para V
z
se obtienen empleando las siguientes relaciones:
j
, m
[ J
z
V
+
V
+
J
z
[j, m) = (m
m) j
, m
[ V
+
[j, m) (1.1)
= j
, m
[ V
+
[j, m) ,
luego los elementos de matriz de V
+
son cero a menos que m
m = m = 1.
En forma analoga se obtienen las relaciones para operadores V
y V
z
, cuyo
resultado es:
j
, m
[ V
+
[j, m) =
_
0 si m ,= 1
,= 0 si m = 1
j
, m
[ V
[j, m) =
_
0 si m ,= 1
,= 0 si m = 1
j
, m
[ V
z
[j, m) =
_
0 si m ,= 0
,= 0 si m = 0
Debemos encontar la dependecia en j de coecientes del operador
vectorial V . Para esto emplearemos la siguiente relacion de conmutacion:
_
J
2
, V
l
= J
k
J
k
V
l
V
l
J
k
J
k
= (J
k
V
l
J
k
+i
klm
J
k
V
m
) (J
k
V
l
J
k
i
klm
V
m
J
k
)
= i
lmk
V
m
J
k
+i
lkm
J
k
V
m
= i
lmk
V
m
J
k
i
lkm
J
m
V
k
= i (V J J V )
l
,
_
J
2
, V
= i (V J J V ) = 2iK, K =
1
2
(V J J V )
(J V +V J)
m
=
mkl
J
k
V
l
+
mkl
V
k
J
l
=
mkl
J
k
V
l
+
mkl
(J
l
V
k
i
lkr
V
r
)
=
mkl
J
k
V
l
+
mlk
(J
k
V
l
i
klr
V
r
)
=
mkl
J
k
V
l
mkl
(J
k
V
l
i
klr
V
r
) = i
klm
klr
V
r
= i (
ll
mr
lr
ml
) V
r
= i (3V
m
lr
V
r
ml
) = 2iV
m
J V +V J = 2iV
K = V J iV = iV J V
Ahora, calcularemos J K:
(J K)
l
=
lmn
J
m
K
n
=
1
2
lmn
J
m
nst
(V
s
J
t
J
s
V
t
)
=
1
2
nlm
nst
J
m
(V
s
J
t
J
s
V
t
)
=
1
2
(
ls
mt
lt
ms
) J
m
(V
s
J
t
J
s
V
t
)
=
1
2
[J
m
V
l
J
m
J
m
J
l
V
m
J
m
V
m
J
l
+J
m
J
m
V
l
]
=
1
2
_
J
m
V
l
J
m
J
m
J
l
V
m
J V J
l
+J
2
V
l
=
1
2
_
J
m
J
m
V
l
i
mlu
J
m
V
u
J
m
V
m
J
l
i
lmu
J
m
V
u
J V J
l
+J
2
V
l
=
1
2
_
i
mlu
J
m
V
u
+i
mlu
J
m
V
u
2J V J
l
+ 2J
2
V
l
J K = J
2
V J V J.
En forma analoga:
K J = J V J V J
2
.
Con estos resultados podemos evaluar
_
J
2
,
_
J
2
, V
= i
_
J
2
, K
= i
2
(K J J K)
=
_
J V J V J
2
J
2
V +J V J
_
=
_
J
2
V 2J V J +V J
2
_
.
Al calcular los coecientes de estos operadores encontramos que:
j
, m
[
_
J
2
,
_
J
2
, V
[j, m) = j
, m
[ J
2
V 2J V J +V J
2
[j, m) ,
si j = j
nos permite obtener:

j, m
[ J
4
V
l
2J
2
V
l
J
2
+V J
4
[j, m) = 0
j, m
[ J
2
V
l
2J V J
l
+V
l
J
2
[j, m) = 0
j, m
[ V [j, m) =
j[ J V [j)
2
j(j + 1)
j, m
[ J [j, m) .
Observe que empleamos la condicion de ser un operador escalar para el pro-
ducto J V , por lo tanto es diagonal en la base [j, m) y su valor propio solo
depende de j. En el caso en que j
,= j en que j
, m
[ J V J [j, m) = 0 se
obtiene
_
j
2
(j
+ 1)
2
+j
2
(j + 1)
2
2j
(j
+ 1) j (j + 1)
2j
(j
+ 1) 2j (j + 1)
_
j
, m
[ V [j, m) = 0
((j
j)
2
1)((j
+j + 1)
2
1) j
, m
[ V [j, m) = 0.
10
Como (j
+j + 1)
2
1 ,= 0, ya que j
,= j se tiene que cuando

j = j
j = 1,
los coecientes j
, m
[ V [j, m) son no nulos. Al igual que en el caso de los

operadores escalares la dependencia explcita de los coecientes se obtiene
a partir de las relaciones de conmutacion de ls operadores V
. Asumiremos
que j = 1, es decir, j
= j +1. Recuerde que ademas deben satisfacerse las

reglas de seleccion para m.
j + 1, m+ 1[ [J
+
, V
+
] [j, m1) = 0
j + 1, m+ 1[ J
+
V
+
V
+
J
+
[j, m1) = 0
j + 1, m+ 1[ J
+
V
+
[j, m1) = j + 1, m+ 1[ V
+
J
+
[j, m1)
J
+
[j, m1) =
_
(j m+ 1) (j +m) [j, m)
J
[j + 1, m+ 1) =
_
(j +m+ 2) (j m+ 1) [j, m)
_
(j +m+ 2) (j m+ 1) j + 1, m[ V
+
[j, m1) =
_
(j m+ 1) (j +m) j + 1, m+ 1[ V
+
[j, m)
j + 1, m[ V
+
[j, m1)
_
(j +m) (j m+ 1)
=
j + 1, m+ 1[ V
+
[j, m)
_
(j +m+ 2) (j m+ 1)
(j)
j + 1, m[ V
+
[j, m1)
_
(j +m) (j m+ 1)
=
j + 1, m+ 1[ V
+
[j, m)
_
(j + 1 m) (j +m+ 2)
.
Esta razon es independiente del n umero cuantico m y la designaremos por
(j), as:
j + 1, m[ V
+
[j, m1) = (j)
_
(j +m) (j m+ 1).
En forma analoga se obtienen los coecientes para V
z
y V
empleando, res-
pectivamente, [J
, V
+
] = 2V
z
y [J
, V
z
] = V
. El resultado para j = 1
se resume en la siguiente tabla:
j = 1 :
j + 1, m+ 1[ V
+
[j, m) = (j)
_
(j +m+ 2) (j m+ 1)
j + 1, m[ V
z
[j, m) = (j)
_
(j +m+ 1) (j m+ 1)
j + 1, m1[ V
[j, m) = (j)
_
(j +m) (j m+ 1)
j = 0 :
j, m+ 1[ V
+
[j, m) = (j)
_
(j m) (j +m+ 1)
j, m[ V
z
[j, m) = (j)m
j, m1[ V
[j, m) = (j)
_
(j +m) (j m+ 1)
j = 1 :
j 1, m+ 1[ V
+
[j, m) = (j)
_
(j m1) (j m)
j 1, m[ V
z
[j, m) = (j)
_
(j m)(j +m)
j 1, m1[ V
[j, m) = (j)
_
(j +m) (j +m1).
1.5. Operadores Tensoriales
En el caso clasico sabemos que los tensores transforman como:
T
ijk...
R
ii
R
jj
. . . T
ijk...
ante una rotacion ortogonal dada por R, donde el n umero de ndices denota
el rango del tensor. En este caso las componentes del tensor se denominan
cartesianas. Mediante el empleo de estas componentes es posible construir
tensores de rango inferior mediante la contraccion de indices, como se realiza
al contruir un escalar. Ademas, mediante la mutiplicacion de tensores de
un determinado rango se pueden construir tensores de rango superior, sin
embargo, estos no poseen simetra denida ante las rotaciones.
En este sentido es conveniente denir operatores tensoriales esfericos
en analoga con los armonicos esfericos Y
m
l
( n) = n[l, m) de acuerdo a:
T
(k)
q
= Y
m=q
l=k
(V ),
j
r
V
j
, con j = x, y, z.
donde k es el rango y q es el n umero cuantico magnetico (m). Los opera-
dores tensoriales esfericos son irreducibles (poseen simetra denida ante R).
Veamos como transforman los armonicos esfericos, si [ n) = T(R) [ n),
n[j, m) = n[T
(R)[j, m)
=
m
j, m
[T
(R)[j, m) n[j, m
)
Y
m
l
( n) =
m
T
(l)
mm
(R
1
)Y
m
l
( n),
donde las funciones de Wigner estan dadas por:
T
(l)
mm
(R) = j, m[ exp
_
i
J n
_
[j, m) .
En forma analoga los operadores tensoriales esfericos (OTE) transforman
como:
T
(R)T
(k)
q
T(R) =
k
q=k
T
(k)
qq
(R)T
(k)
q
,
Y
m
l
( n) =
m
T
(l)
mm
(R
1
)Y
m
l
( n),
12
Si se asume que la rotacion es innitesimal los operadores tensoriales se de-
nen como:
_
J n, T
(k)
q
_
=
k
q=k
k, q[ J n[k, q) T
(k)
q
. (1.2)
_
J
z
, T
(k)
q
_
= qT
(k)
q
.
_
J
, T
(k)
q
_
=
_
(k q)(k q + 1)T
(k)
q1
.
La expresion (1.2) es la denicion alternativa para los OTE. Ejemplos de ellos
son:
T
(0)
0
=
U V
3
, T
(1)
q
=
(U V )
q
i
2
,
T
(2)
2
= U
1
V
1
, T
(2)
1
=
U
1
V
0
+U
0
V
1
2
,
T
(2)
0
=
U
+1
V
1
+ 2U
0
V
0
+U
1
V
1
6
Teorema: Si X
(k
1
)
q
1
y X
(k
2
)
q
2
son tensores irreducibles de rango k
1
y k
2
,
respectivamente. Entonces el tensor
T
(k)
q
=
q
1
,q
2
k
1
, k
2
; q
1
, q
2
[k
1
, k
2
; k, q) X
(k
1
)
q
1
X
(k
2
)
q
2
es irreducible de rango k. Los coecientes k
1
, k
2
; q
1
, q
2
[k
1
, k
2
; k, q) son los
coecientes de Clebsch-Gordan.
Teorema de Wigner-Eckart: Los elementos de matriz de un ope-
rador tensorial en la base de momentum angular satisfacen:
_
; jm[T
(k)
q
[; j, m
_
= j, k; m, q[j, k; jk)
_
; j[[T
(k)
q
[[; j
_
2j + 1
.
La demostracion es directa empleando la denicion de operador tensorial en
terminos de relaciones de conmutacion con las componentes de J. Ademas,
empleando la relacion de communtacion con J
z
obtenemos la regla de selec-
cion en m.
_
; jm[T
(k)
q
[; j, m
_
= 0 a menos que m= q +m.
1.6. Tarea 1
1. Demostrar que los elementos de matriz de matriz de un operador vecto-
rial, para los casos en que j = 0, 1, estan dados por:
j = 0 :
j, m+ 1[ V
+
[j, m) = (j)
_
(j m) (j +m+ 1)
j, m[ V
z
[j, m) = (j)m
j, m1[ V
[j, m) = (j)
_
(j +m) (j m+ 1)
j = 1 :
j 1, m+ 1[ V
+
[j, m) = (j)
_
(j m1) (j m)
j 1, m[ V
z
[j, m) = (j)
_
(j m)(j +m)
j 1, m1[ V
[j, m) = (j)
_
(j +m) (j +m1).
donde (j) y (j) son constantes que dependen solo del n umero cuantico
j, es decir, son independientes del valor de m.
Ademas, si V es hermtico, V
+
= V
, las funciones (j) y (j) estan

relacionadas por:
(j, , ) =
(j 1, , )
2. Problema 22, captulo 3, Sakurai, Modern Quantum Mechanics.
3. Problema 28, captulo 3, Sakurai, Modern Quantum Mechanics.
4. El operador s
r
= s
r
r
es la componente del spin del electron en la
direccion del vector r. Dado que s
r
es un operador escalar, conmuta
con cada componente del momentum angular total J = l + s, por lo
tanto [J, s
r
] = 0. Vericar la relacion de conmutacion anterior en forma
directa calculando los conmutadores [l, s
r
] y [s, s
r
] y las relaciones de
conmutacion de l y r.
2. Teora de Perturbaciones
Resumen En Mecanica Cuantica es posible resolver en forma analtica cerrada un
peque no n umero de problemas. Sin embargo, en una familia importante de proble-
mas se puede resolver en forma aproximada, estudiandolo como una perturbacion a
un problema conocido. En este captulo estudiaremos diversos metodos para abor-
dar en forma sistematica este tratamiento tanto para potenciales independientes
del tiempo como para aquellos en que depende explcitamente del tiempo.
2.1. Introduccion
Se desea determinar la evolucion de un sistema cuantico cuya solucion
exacta es conocida en ausencia de perturbaciones.
H = H
0
+V,
donde H
0
es el Hamiltoniano con solucion conocida, es decir, se han deter-
minado sus valores y vectores propios E
(0)
n
y
(0)
n
_
, respectivamente,
H
0
(0)
n
_
= E
(0)
n
(0)
n
_
y V es la perturbacion. En general, la pertubacion puede ser tiempo dependienteV =
V (t).
La perturbacion, en general, determina un desplazamiento de los nive-
les de energa del sistema original, obteniendose por lo tanto una nueva base,
asociada los nuevos niveles energeticos.
2.2. V independiente del tiempo
Al igual que en el caso no perturbado debemos resolver el problema
de valores y estados propios del sistema, es decir,
H = H
o
+V (2.1)
H[
n
i
) = E
n
[
n
i
) . (2.2)
2. Teora de Perturbaciones 15
Observe, que en general hemos asumido que puede existir mas de un estado
propio que tenga el mismo valor propio. Es conveniente reescribir la pertur-
bacion como
V = W, (2.3)
donde es el parametro perturbativo, es decir, 1 con lo cual podemos
permitir que W) H
o
). Por ejemplo, una partcula cargada en un potencial
armonico que es sometida a un campo electrico. El parametro no siempre
es evidente en un problema dado, en general, se debe construir la forma (2.3).
2.2.1. Caso de espectro no degenerado
Al expresar nuestro problema en la forma (2.2), podemos efectuar una
expansion en serie de potencias en ambos lados, en funcion del parametro
H[
n
) = E
n
[
n
) ,
n
= E
n
E
(0)
n
,
donde
n
representa el corrimiento de energa del nesimo nivel debido a
la perturbacion. Como los vectores propios de H
o
forman una base podemos
escribir:
[
n
) =
j
c
j
(0)
j
_
expansionenlabasenoperturbada (2.4)
[
n
) =
(0)
n
_
+
_
(1)
n
_
+
2
(2)
n
_
+. . .
_
(2.5)
n
=
(1)
n
+
2
(2)
n
+. . . (2.6)
Reemplazando estas deniciones en la ecuacion original se tiene que:
(H
o
+W) [
n
) =
_
E
(0)
n
+
n
_
[
n
)
(W
n
) [
n
) =
_
E
(0)
n
H
o
_
[
n
)
En esta ultima ecuacion el termino al lado izquierdo es ortogonal a
(0)
n
_
ya
que el menor termino es de orden , en forma alternativa podemos ver que
_
(1)
n
_
+
2
(2)
n
_
+. . .
_
=
j=n
c
j
(0)
j
_
.
Multiplicando este termino por la identidad menos el proyector
(0)
n
__
(0)
n
,
es decir
n
= I
(0)
n
__
(0)
n
=

j=n
(0)
j
__
(0)
j
, se obtiene que
n
(W
n
) [
n
) = (W
n
) [
n
) .
16
Esto nos permite eliminar las contribuciones del estado no perturbado
(0)
n
_
al estado perturbado luego [
n
), entonces podemos escribir:
[
n
) =
(0)
n
_
+
_
(1)
n
_
+
2
(2)
n
_
+. . .
_
=
(0)
n
_
+
1
E
(0)
n
H
o
n
(W
n
) [
n
) .
Observe que el estado anterior no esta normalizado, esto se debe efectuar al
terminar el calculo. Ademas, la correccion de primer orden en la energa se
obtiene directamente observando que:
_
(0)
n
(W
n
) [
n
) = 0
n
=
_
(0)
n
W [
n
) .
Esta ultima expresion determina los corrimientos de energa al reemplazar
las ecuaciones (2.5,2.6):
(1)
n
=
_
(0)
n
(0)
n
_
= W
nn
,
(2)
n
=
_
(0)
n
(1)
n
_
,
.
.
.
(m)
n
=
_
(0)
n
(m1)
n
_
.
Con esta informacion podemos determinar las correcciones a los estados, para
esto vemos que:
(0)
n
_
+
(1)
n
_
+
2
(2)
n
_
+. . .
=
(0)
n
_
+

E
(0)
n
H
o
n
_
W
(1)
n

(2)
n
. . .
__
(0)
n
_
+
(1)
n
_
+
2
(2)
n
_
+. . .
_
Agrupando terminos correspondientes a potencias iguales en se obtiene:
(1)
n
_
=
1
E
(0)
n
H
o
n
_
W
(1)
n
_
(0)
n
_
=
1
E
(0)
n
H
o
j=n
(0)
j
__
(0)
j
W
(1)
n
(0)
n
_
como
n
(0)
n
_
= 0,
(1)
n
_
=
1
E
(0)
n
H
o
n
W
(0)
n
_
=
j =n
_
(0)
j
(0)
n
_
E
(0)
n
E
(0)
j
(0)
j
_
=
j=n
W
jn
E
(0)
n
E
(0)
j
(0)
j
_
.
Ahora con esta correccion al estado determinamos en forma explcita la se-
gunda correccion a la energa como:
(2)
n
=
_
(0)
n
(1)
n
_
=
j=n
W
jn
E
(0)
n
E
(0)
j
_
(0)
n
(0)
j
_
=
j=n
[W
jn
[
2
E
(0)
n
E
(0)
j
,
n
= W
nn
+
2
j=n
W
jn
E
(0)
n
E
(0)
j
_
(0)
n
(0)
j
_
+...
= W
nn
+
2
j=n
[W
jn
[
2
E
(0)
n
E
(0)
j
+...
De la misma forma anterior obtenemos la segunda correccion al estado como:
(2)
n
_
=
1
E
(0)
n
H
o
n
_
W
(1)
n
_
(1)
n
(1)
n
_
(2)
n
(0)
n
__
.
El ultimo termino se cancela ya que
n
es un proyector ortogonal a
(0)
n
__
(0)
n
,
por lo tanto
(2)
n
_
=
1
E
(0)
n
H
o
n
_
W
1
E
(0)
n
H
o
n
W
W
nn
E
(0)
n
H
o
n
W
_
(0)
n
_
=
1
E
(0)
n
H
o
n
_
_
j=n
W
jn
E
(0)
n
E
(0)
j
W
W
nn
W
jn
E
(0)
n
H
o
_
_
(0)
j
_
=
k,j=n
W
k,j
W
jn
(0)
k
_
_
E
(0)
n
E
(0)
k
__
E
(0)
n
E
(0)
j
_
j=n
W
nn
W
jn
(0)
j
_
_
E
(0)
n
E
(0)
j
_
2
.
El procedimiento se puede continuar hasta obtener las soluciones al orden
deseado. Sin embargo, a un esta pendiente la normalizacion del estado [
n
)
18
que hemos obtenido:
[
n
)
N
= N
1/2
n
[
n
)
N
1
n
=
n
[
n
) = 1 +
2
j=n
[V
jn
[
2
_
E
(0)
n
E
(0)
j
_
2
+(
3
) +. . .
N
n
1
2
j=n
[V
jn
[
2
_
E
(0)
n
E
(0)
j
_
2
.
Problema 1. Oscilador armonico perturbado En el caso del oscilador armoni-
co perturbado consideraremos
H = H
o
+W,
H
o
=
1
2m
P
2
+
1
2
m
2
o
X
2
,
W =
1
2
m
2
0
X
2
.
Este problema resulta especialmente ilustrativo, pues obviamente se puede
resolver en forma exacta mediante una redenicion de la frecuencia del osci-
lador como
o
=
o
1 +,
lo que emplearemos posteriormente para comprobar nuestro resultado. En
la forma usual se denen operadores de subida y bajada para el oscilador
armonico
a =
_
m
o
2
_
X +
iP
m
o
_
, a
=
_
m
o
2
_
X
iP
m
o
_
H =
o
( n + 1/2) +
o
4
_
a +a
_
2
.
Analicemos el cambio del estado fundamental, as,
[0) =
0
(0)
_
+
j=n
W
jn
_
E
(0)
n
E
(0)
j
_
j
(0)
_
+. . .
=
0
(0)
_
+
j=n
W
jn
o
(n j)
j
(0)
_
+. . .
0
= W
00
+
2
j=n
[W
jn
[
2
_
E
(0)
n
E
(0)
j
_ +. . .
= W
00
+
2
j=n
[W
jn
[
2
o
(n j)
+. . .
W
00
=
_
0
(0)
0
(0)
_
=

o
4
_
0
(0)
_
a
2
+a
2
+ 2 n + 1
_
0
(0)
_
=

o
4
W
20
=
_
2
(0)
0
(0)
_
=

o
4
_
2
(0)
_
a
2
+a
2
+ 2 n + 1
_
0
(0)
_
=

o
4
2
[0) =
0
(0)
_
1
4
2
(0)
_
+. . . ,
0
=
o
_
4

2
16
+. . .
_
Como mencionamos anteriormente el sistema perturbado se puede resolver
en forma analitica cerrada con el cambio
o

H = ( n + 1/2) , E
n
= (n + 1/2)
=
o
_
1 +

2

2
8
+O
_
3
_
_
E
n
= E
(0)
n
+
o
_
1 +

2

2
8
+O
_
3
_
_
n
= E
n
E
(0)
n
=
o
(n + 1/2)
_
2

2
8
+O
_
3
_
_
0
=
o
_
4

2
16
+O
_
3
_
_
Problema 2. Determinar las correcciones hasta segundo orden de los estados
excitados del oscilador armonico cuando la perturbacion es del tipo W =
1
2
m
2
0
a
2
X
4
.
1. Como determinar el tiempo de validez de la aproximacion?
2. Es necesario determinar las correcciones al mismo orden en los estados
y los corrimientos en las energas?
Problema 3. Partcula en un campo electrico uniforme (atomo de Hidrogeno
o hidrogenoide) o efecto Stark cuadratico.
En este caso el Hamiltoniano del sistema no perturbado junto a la
peturbacion se pueden escribir como
H = H
o
+V
H
o
=
1
2m
P
2
+V
o
(r), V
o
(r) =
Ze
2
r
V = e [E[ z,
donde e es la carga electrica del electron, Z es el n umero de protones y
[E[ es la amplitud del campo electrico uniforme. En coordenadas esfericas el
20
problema no perturbado esta dado por:
2
2m
_
1
r
2
r
_
r
2
E
r
_
+
1
r
2
sin
_
sin
_
+
1
r
2
sin
2
2
_
Ze
2
r

E
= E
E
La funcion propia se descompone en su parte radial y angular
l,n,m
= R
n,l
(r)Y
m
l
(, ) = r[)
R
n,l
(r) =
_
_
2Z
na
o
_
3
(n l 1)!
2n[(n +l)!]
3
_
1/2
e
/2
L
2l+1
2+l
(), =
2Z
na
o
r
a
o
=

2
m
e
e
2
, E
n
=
Z
2
e
4
2n
2
a
o
, n l + 1.
No se considera degeneracion de los estados debido al spin del electron, el
cual se analizara posteriormente. De esta forma el corrimiento de energa de
los estados es:
n
= E
n
E
(0)
n
= e [E[ z
(n,l,m;n,l,m)
+e
2
[E[
2
j=n
z
(n,l,m;j,l
,m
2
_
E
(0)
n
E
(0)
j
_ +. . .
Los elementos matriciales del operador de posicion z son
z
(n,l,m;j,l
,m
)
= n, l, m[z[n
, l
, m
) =
_
_
_
,= 0 si
_
l
= l 1
m
= m
0 enotroscasos
Empleando la regla de seleccion anterior podemos ver que la primera correc-
cion a la energa se cancela y como menor orden se mantiene solo el segundo
termino. Consideremos el caso del estado fundamental [1, 0, 0), por lo tan-
to los estados que contribuyen seran aquellos del tipo [j, 1, 0) con lo que se
obtiene
1
e
2
[E[
2
j=n
z
(1,0,0;j,1,0)
2
_
E
(0)
1
E
(0)
j
_
Para obtener una expresion aproximada debemos reemplazar el termino de-
pendiente de la energa
E
(0)
1
+E
(0)
j
E
(0)
1
+E
(0)
2
=
e
2
2a
o
_
1
1
4
_
=
3e
2
8a
o
1

8
3
a
o
[E[
2
j=n
z
(1,0,0;j,1,0)
2
,
como ahora no hay problemas de divergencia en la suma
z
(1,0,0;j,1,0)
j
1, 0, 0[z[j
, l
, m
)j, l
, m
[z[1, 0, 0)
= 1, 0, 0[z
2
[1, 0, 0) =
1
3
1, 0, 0[r
2
[1, 0, 0) = a
2
o
.
Luego
1

8
3
a
3
o
[E[
2
.
Ademas, la polarizabilidad del atomo se dene en funcion de la respuesta de
un medio a la presencia del campo electrico, generando un corrimiento de
energa:
1

1
2
[E[
2
, =
16
3
a
3
o
.
Problema 4. Considere un sistema fsico de tres niveles descrito por
H =
0
_
_
5
3
1
_
_
+
0
_
_
0 1 2
1 0 1
2 1 0
_
_
,
de modo que 1. Determine las correciones de segundo en las energas pro-
pias y a primer orden en los estados propios. Alternativamente, el problema
se puede escribir como:
H =
0
(5 [3) 3[ + 3 [2) 2[ +[1) 1[) +
0
([3) 2[ +[2) 3[ + 2 ([3) 1[ +[1) 3[) +[2) 1[ +[1) 2[) ,
donde los estados son los elementos de la base no perturbada, por lo tanto
las correcciones para el estado fundamental estan dadas por:
(1)
1
_
=
W
21
E
(0)
1
E
(0)
2
2
(0)
_
+
W
31
E
(0)
1
E
(0)
3
3
(0)
_
=
1
2
_
2
(0)
_
+
3
(0)
__
[
1
)
N
=
1
_
1 +

2
4
)
_
1
(0)
_

2
_
2
(0)
_
+
3
(0)
__
_
(1)
1
= 1[ W [1) = 0
(2)
1
=
_
1
(0)
(1)
1
_
=
j=1
[W
j1
[
2
E
(0)
1
E
(0)
j
=
0
_
1
2
+
4
4
_
=
3
0
2
.
Problema 5. Comparar la evolucion del sistema descrito en el ejemplo ante-
rior empleando la solucion exacta y la aproximada, a diferentes ordenes en ,
22
para esto considere valores numericos de = 0,1 y 0,2. Por ejemplo, graque
la evolucion del operador Z
31
= [3) 3[ [1) 1[ . Que puede indicar sobre las
preguntas ?
H =
0
_
_
5 2
3
2 1
_
_
,
cuando el sistema se encuentra inicialmente en el estado [) =
_
1
(0)
_
+ 2
2
(0)
_
+
3
(0)
__
/2.
Que puede indicar sobre las preguntas a,b para este problema en particular?
2.2.2. Caso de espectro degenerado
El metodo desarrollado en el caso no degenerado considera la existen-
cia de un unico estado asociado a cada valor de la energa, por lo tanto es
posible calcular el inverso del operador
_
E
(0)
n
H
o
_
en un subespacio orto-
gonal, pues no apareceran divergencias. Aplicaremos esta misma idea en el
caso degenerado, para lo cual consideraremos los siguientes aspectos:
Asumiremos que el espectro de energa es degenerado en el subespacio
D, cuya dimension es g, es decir existen g estados con energa E
(0)
D
.
Sea
(0)
m
_
la base del subespacio degenerado D del problema no per-
turbado.
Sea
(0)
l,
_
la base del subespacio degenerado D del problema pertur-
bado. Naturalemente, se debe satisfacer la condicion
(0)
l,
_
lm0

(0)
l
_

(0)
m
_
.
La base
(0)
l
_
, en general, no coincide con la base
(0)
m
_
, mas bien
debe esperarse que sean combinaciones lineales de esta ultima base, es
decir,

(0)
l
_
=
mD
(0)
m
(0)
l
_
(0)
m
_
.
Observe, que en general los elementos
(0)
l,
_
tendran componentes fuera
del subespacio D.
Sea P
D
un operador de proyeccion en D
P
D
=
mD
(0)
m
__
(0)
m
,
por lo tanto la proyeccion ortogonal es
P
= 1 P
D
.
Al igual que en el caso no degenerado debemos resolver el problema
de valores y vectores propios que esta dado por:
(H
o
+W) [
l,
) = E
l
[
l,
) (E
l
H
o
W) [
l,
) = 0.
(E
l
H
o
W) (P
D
+P
) [
l,
) =
(E
l
H
o
W) P
D
[
l,
) + (E
l
H
o
W) P
[
l,
) = 0.
Es posible separar las contribuciones al estado perturbado en sus componen-
tes interior y exterior al subespacio D, mediante la aplicacion de los corres-
pondientes proyectores ortogonales, es decir,
P
D
(E
l
H
o
W) P
D
[
l,
) +P
D
(E
l
H
o
W) P
[
l,
) = 0
(P
D
(E
l
H
o
) P
D
W) P
D
[
l,
) +P
D
(E
l
H
o
) P
[
l,
) P
D
WP
[
l,
) = 0
__
E
l
E
(0)
D
_
P
D
P
D
W
_
P
D
[
l,
) +
_
E
l
E
(0)
D
_
P
D
P
[
l,
) P
D
WP
[
l,
) = 0
_
E
l
E
(0)
D
_
P
D
[
l,
) P
D
WP
D
[
l,
) P
D
WP
[
l,
) = 0 (2.7)
P
(E
l
H
o
W) P
D
[
l,
) +P
(E
l
H
o
W) P
[
l,
) = 0
(P
(E
l
H
o
) P
W) P
D
[
l,
) +P
(E
l
H
o
) P
[
l,
) P
WP
[
l,
) = 0
((E
l
H
o
) P
W) P
D
[
l,
) + (E
l
H
o
) P
2
[
l,
) P
WP
[
l,
) = 0
(E
l
H
o
P
W) P
[
l,
) P
WP
D
[
l,
) = 0 (2.8)
Lo que deseamos obtener es:
[
l,
) = (P
+P
D
) [
l,
) = P
[
l,
) +P
D
[
l,
)
=
(0)
l
_
+
(1)
l,
_
+
2
(2)
l,
_
+. . .
Empleando la ecuacion (2.8) se obtiene que
_
E
l
E
(0)
D
P
D
W
_
P
D
[
l,
) P
D
WP
[
l,
) = 0
(E
l
H
o
P
WP
) P
[
l,
) P
WP
D
[
l,
) = 0
P
[
l,
) = P
E
l
H
o
P
WP
WP
D
[
l,
)
En el caso de la correcion de primer orden se tiene que:
P
(1)
l,
_
=
j / D
W
jl
E
(0)
D
E
(0)
j
(0)
j
_
.
24
Con este resultado y empleando (2.7) podemos determinar la componente de
la primera correccion en el subespacio D.
_
E
l
E
(0)
D
P
D
WP
D
_
P
D
[
l,
) P
D
WP
[
l,
) = 0
__
E
l
E
(0)
D
_
P
D
WP
D
_
P
D
[
l,
) P
D
WP
E
l
H
o
P
WP
WP
D
[
l,
) = 0
_
_
E
l
E
(0)
D
_
P
D
WP
D
P
D
WP
E
l
H
o
P
WP
W
_
P
D
[
l,
) = 0 (2.9)
El ultimo termino en la expresion (2.9) es de orden
2
, por lo tanto a primer
orden se obtiene la siguiente ecuacion de estados y valores propios
__
E
l
E
(0)
D
_
P
D
WP
D
_
P
D
[
l,
) = 0,
l
= E
l
E
(0)
D
que determinan
(1)
l
y
(0)
l
_
a orden cero al resolver:
det
_
E
l
E
(0)
D
V
_
= 0.
En forma natural la perturbacion es diagonal en la nueva base:
(1)
l
=
_
(0)
l
(0)
l
_
.
Nuevamente volvemos a la ecuacion (2.9) para determinar la correccion de
primer orden:
_
_
E
l
E
(0)
D
_
P
D
WP
D
2
P
D
WP
1
E
l
H
o
P
WP
W
_
P
D
[
l,
) = 0
_
_
E
l
E
(0)
D
_
P
D
WP
D
2
P
D
WP
1
E
(0)
D
H
o
P
W
_
P
D
[
l,
) = 0
__
E
l
E
(0)
D
_
P
D
WP
D
_
P
D
[
l,
) =
2
P
D
WP
1
E
(0)
D
H
o
P
WP
D
[
l,
)
_
v
j
P
D
WP
D
_
P
D
(1)
l
j
_
= P
D
WP
1
E
(0)
D
H
o
P
l
j
_
P
D
(1)
j
_
=
1
_
v
j
P
D
WP
D
_P
D
WP
1
E
(0)
D
H
o
P
j
_
P
D
(1)
j
_
=
i=j
i
[ WP
1
E
(0)
D
H
o
P
j
_
_
v
j
v
i
_ [
i
) .
La correccion de segundo orden a la energa esta dada por:
(2)
l
=
_
(0)
l
(1)
l
_
=
j / D
[W
jl
[
2
E
(0)
D
E
(0)
j
.
Problema 6. Considere un sistema tridimensional sometido a una perturba-
cion
V = c [0) 1[ +c
[1) 0[ +a [0) 2[ +a
[2) 0[ .
Para este sistema se conoce que:
H
o
= E
0
[0) 0[ +E
0
[1) 1[ +E
2
[2) 2[ .
Se desea determinar la primera correccion a la energa y a los estados propios.
Solucion: Como sabemos este sistema se puede resolver en forma exac-
ta, pero emplearemos la teora de perturbaciones en el caso degenerado. Ob-
serve que
(H EI) =
_
_
E
0
E c a
c E
0
E 0
a 0 E
2
E
_
_
En este caso para el subespacio degenerado:
P
D
= [0) 0[ +[1) 1[ , P
= [2) 2[
(0)
1
_
=
1
2
([0) +[1)) , v
1
= E
0
+c
(0)
2
_
=
1
2
([0) +[1)) , v
2
= E
0
c
P
D
V P
D
= V, det (E
l
E
0
P
D
V P
D
) = 0
(0)
1
_
=
1
2
([0) +[1)) , v
1
= E
0
+c
(0)
2
_
=
1
2
([0) +[1)) , v
2
= E
0
c
P
(1)
1
_
=
V
21
E
0
E
2
[2) =
a/
2
E
0
E
2
[2) ,
P
(1)
2
_
=
V
22
E
0
E
2
[2) =
a/
2
E
0
E
2
[2) ,
P
D
(1)
1
_
=
_
(0)
2
V [2) 2[
1
E
(0)
D
H
o
[2) 2[ V
(0)
2
_
(v
1
v
2
)
(0)
2
_
=
a
2
2 2c (E
0
E
2
)
(0)
2
_
P
D
(1)
2
_
=
_
(0)
1
V [2) 2[
1
E
(0)
D
H
o
[2) 2[ V
(0)
1
_
(v
2
v
1
)
(0)
1
_
=
a
2
2 2c (E
0
E
2
)
(0)
1
_
26
Problema 7. Un sistema fsico puede ser modelado empleando un espacio de
Hilbert de 4 dimensiones cuyo Hamiltoniano no pertubado esta dado por:
H
o
= E
1
[1) 1[ +E
1
[2) 2[ +E
3
[3) 3[ +E
3
[4) 4[ . (2.10)
Este sistema es sometido a una perturbacion descrita por:
V = a [1) 2[ +a [2) 1[ +b [3) 4[ +b [4) 3[ , (2.11)
donde a y b son constantes reales.
Determine las correcciones de primer orden a las energas del sistema
y a los estados [1) y [4).
2.3. Teora de perturbaciones dependiente del tiempo
En este caso asumiremos que la perturbacion posee una dependencia
explcita del tiempo, es decir, V se puede escribir como:
V (t) = W(t).
Al igual que en el caso no dependiente del tiempo asumiremos conocida la
solucion del problema no perturbado. En el caso no perturbado la evolucion
temporal del vector de estado se escribe como:
H
o
[
n
) = E
n
[
n
) ,
[(t)) = U
o
(t) [(0)) = exp(i
Ht
) [(0)) ,
con [(0)) =

n
c
n
[
n
), donde c
n
son constantes tal que

n
[c
n
[
2
= 1. Por
lo tanto se obtiene que:
[(t)) =
n
c
n
e
iE
n
t/
[
n
) .
En el caso perturbado escribiremos la solucion de la siguiente forma:
[(t)) =
n
c
n
(t)e
iE
n
t/
[
n
) .
Observe que el termino exponencial e
iE
n
t/
permite eliminar las oscilaciones
rapidas en el coeciente asociado al estado [
n
). Se debe satisfacer:
i
[(t))
t
= (H
o
+V (t)) [(t))
i
_
_
c
n
(t) i
E
n
c
n
(t)
_
e
iE
n
t/
[
n
)
_
=
m
(E
m
+V (t)) c
m
(t)e
iE
m
t/
[
m
) /
n
[
i
_
n
_
c
n
(t) i
E
n
c
n
(t)
_
e
iE
n
t/
n
[ [
n
)
_
=
n
[ (E
m
+V (t)) c
m
(t)e
iE
m
t/
[
m
)
[i c
n
(t) +E
n
c
n
(t)] e
iE
n
t/
=
m
(E
m
n.m
+
n
[ V (t) [
m
)) c
m
(t)e
iE
m
t/
i c
n
(t) =
n
[ V (t) [
m
) c
m
(t)e
i(E
m
E
n
)t/
. (2.12)
Esto equivale a un sistema de ecuaciones lineales acopladas de primer orden.
Si denominamos frecuencias de transicion entre los niveles n m como
nm
= (E
n
E
m
) /. y se dene el vector x y la matriz A como:
x
j
= c
j
, A
nm=
n
[ V (t) [
m
) c
m
(t)e
i(E
m
E
n
)t/
,
nuestro problema se reduce a resolver la ecuacion
dx
dt
=
i
Ax.
2.3.1. Cuadro de Interaccion
El sistema de ecuaciones corresponde a la solucion en el cuadro de
Schrodinger (S). Sin embargo, en general para este tipo de problema es con-
veniente denir un cuadro de interaccion (I). Recuerde que en el cuadro de
Heisenberg (H) evolucionan los observables, mientras que en Schrodinger la
evolucion corresponde a los estados. El cuadro de interaccion, como veremos,
sera una cuadro intermedio en que ambas cantidades evolucionan.
La motivacion para introducir este cambio nuevamente considera el
conocimiento de la solucion de un problema cuantico en forma exacta, H
o
.
Esto nos permitira encontrar soluciones mas simples para el sistema de ecua-
ciones (2.12).
Asumiremos que a tiempo t = t
o
el sistema esta descrito por [)
S
=
[(t
o
)). As, se dene
[(t, t
o
))
I
= e
iH
o
t/
[(t, t
o
))
S
,
el estado del sistema en el cuadro de interaccion. Observe que el objetivo
que agregar U
o
al estado en el cuadro de Schrodinger es eliminar las osci-
laciones rapidas del sistema (considere en el caso que no existe perturba-
cion, [(t, t
o
))
I
= e
iH
o
t/
[(t, t
o
))
S
= e
iH
o
t/
e
iH
o
t/
[(t
o
))
S
= [(t
o
))
S
).
28
Emplearemos estas deniciones para derivar la ecuacion de movimiento de
[(t, t
o
))
I
. Observemos que:
i
[(t, t
o
))
I
t
= i

t
_
e
iH
o
t/
[(t, t
o
))
S
_
= e
iH
o
t/
H
o
[(t, t
o
))
S
+ie
iH
o
t/

t
[(t, t
o
))
S
= e
iH
o
t/
H
o
[(t, t
o
))
S
+e
iH
o
t/
(H
o
+V ) [(t, t
o
))
S
= e
iH
o
t/
V [(t, t
o
))
S
= e
iH
o
t/
V e
iH
o
t/
e
iH
o
t/
[(t, t
o
))
S
= V
I
(t, t
o
) [(t, t
o
))
I
,
donde
i
[(t, t
o
))
I
t
= V
I
(t, t
o
) [(t, t
o
))
I
, V
I
(t, t
o
) = e
iH
o
t/
V e
iH
o
t/
.
As, vemos que la interaccion evoluciona con H
o
y el estado evoluciona con
V
I
. Resulta directo demostrar que la ecuaciones diferenciales que satisfacen
los coecientes del estado en terminos de la base no perturbada es la misma
obtenida anteriormente, es decir, la ecuacion (2.12).
Podemos obtener una expresion para la evolucion de un operador del
sistema, generalizando la expresion obtenida para el operador de evolucion
comopara esto recordemos que:
A
I
(t, t
o
) = e
iH
o
t/
A
S
e
iH
o
t/
,
dA
I
dt
=
i
_
H
o
e
iH
o
t/
A
S
e
iH
o
t/
e
iH
o
t/
A
S
e
iH
o
t/
H
o
_
=
1
i
[A
I
, H
o
] .
2.3.2. Oscilaciones de Rabi
Este problema corresponde a un sistema de dos niveles acoplado a una
perturbacion que dependiente del tiempo:
H
o
= E
o
[0) 0[ +E
1
[1) 1[
V (t) = e
it
[0) 1[ +e
it
[1) 0[ .
Este potencial de interaccion induce transiciones entre los dos niveles del
sistema cuantizado. Observe que este problema corresponde a la interaccion
de un atomo con un campo electromagnetico clasico, cuya amplitud es y
frecuencia es . Se espera que la interaccion produzca transiciones cuando se
cumpla la condicion de resonancia:

E
1
E
0
.
El sistema de ecuaciones para los coecientes es:
i c
0
(t) = e
it
c
1
(t)e
i(E
1
E
0
)t/
,
i c
1
(t) = e
it
c
0
(t)e
i(E
1
E
0
)t/
.
Se efect ua el siguiente cambio de variables:
b
0
= c
0
e
i
2
(
E
1
E
0
)t
, b
1
= c
1
e
i
2
(
E
1
E
0
)t
b
0
=
i
2
_

E
1
E
0
_
b
0
i
b
1
b
1
=
i
2
_

E
1
E
0
_
b
1
i
b
0
b
1
=
_
_
(E
1
E
0
)
2
_
2
+

2
2
_
b
1
b
1
= Acos(t +), =
_
_
(E
1
E
0
)
2
_
2
+

2
2
_
1/2
.
Asumiremos que a tiempo t = 0 el sistema se encontraba en su estado fun-
damental, por lo tanto
c
1
(0) = 0 b
1
(0) = 0, = /2
b
1
=
i
sin(t).
[b
1
[
2
=
2
_
(E
1
E
0
)
2
_
2
+

2
2
sin
2
_
_
_
_
(E
1
E
0
)
2
_
2
+

2
2
_
1/2
t
_
_
Resulta instructivo escribir este problema empleando matrices de Pauli, ya
que el Hamiltoniano no perturbado es, esencialmente, proporcional a
z
:
H
o
=
(E
1
+E
0
)
2
I +
(E
1
E
0
)
2

z
,
pues el primer termino solo dene la energa de referecnia. Como la interac-
cion esta dada por:
x
= [1) 0[ +[0) 1[ ;
y
= i [1) 0[ +i [0) 1[
[1) 0[ =

x
+i
y
2
; [0) 1[ =

x
i
y
2
V (t) = e
it
[0) 1[ +e
it
[1) 0[ = e
it
x
i
y
2
+e
it
x
+i
y
2
=
_
e
it
x
i
y
2
+e
it
x
+i
y
2
_
= (cos (t)
x
+ sin(t)
y
) .
Problema 8. Empleando las propiedades de las matrices de Pauli determine
V
I
(t). Graque las probabilidades de que el sistema se encuentre en el nivel
inferior y superior. Interprete adecuadamente sus resultados.
30
Problema 9. Considere un sistema de tres niveles que es sometido a la una
interaccion del tipo:
H
o
= E
o
[0) 0[ +E
1
[1) 1[ +E
2
[2) 2[
V (t) =
_
e
i
1
t
[0) 1[ +e
i
1
t
[1) 0[
_
+
_
e
i
2
t
[0) 2[ +e
i
2
t
[2) 0[
_
.
Esto corresponde, por ejemplo, a un atomo de tres niveles interactuando con
dos campos elemtromagneticos clasicos. Determine el operador de evolucion
del sistema y discuta sobre posibles condiciones de resonancia.
2.3.3. Series de Dyson
Nuevamente consideramos una serie perturbativa, en este caso para
los coecientes que dependen del tiempo:
c
n
(t) = c
(0)
n
+c
(1)
n
+c
(2)
n
+. . .
La dinamica del operador de evolucion en el cuadro de interaccion se escribe
como:
i
dU
I
(t, t
o
)
dt
= V
I
(t, t
0
)U
I
(t, t
0
)
U
I
(t, t
o
) = 1 +
1
i
_
t
t
o
V
I
(t
1
)U
I
(t
1
, t
0
)dt
1
= 1 +
1
i
_
t
t
o
V
I
(t
1
)
_
1 +
1
i
_
t
1
t
o
V
I
(t
2
)U
I
(t
2
, t
0
)dt
2
_
dt
1
= 1 +
1
i
_
t
t
o
V
I
(t
1
)dt
1
+
1
(i)
2
_
t
t
o
_
t
1
t
o
V
I
(t
1
)V
I
(t
2
)dt
2
dt
1
. . . +
1
(i)
n
_
t
t
o
V
I
(t
1
)dt
1
_
t
1
t
o
V
I
(t
2
)dt
2
. . .
_
t
n1
t
o
V
I
(t
n1
)dt
n1
Nuevamente consideramos una serie perturbativa, en este caso para
los coecientes que dependen del tiempo:
c
n
(t) = c
(0)
n
+c
(1)
n
+c
(2)
n
+. . .
La dinamica del operador de evolucion en el cuadro de interaccion se escribe
como:
i
dU
I
(t, t
o
)
dt
= V
I
(t, t
0
)U
I
(t, t
0
)
U
I
(t, t
o
) = 1 +
1
i
_
t
t
o
V
I
(t
1
)U
I
(t
1
, t
0
)dt
1
= 1 +
1
i
_
t
t
o
V
I
(t
1
)
_
1 +
1
i
_
t
1
t
o
V
I
(t
2
)U
I
(t
2
, t
0
)dt
2
_
dt
1
= 1 +
1
i
_
t
t
o
V
I
(t
1
)dt
1
+
1
(i)
2
_
t
t
o
_
t
1
t
o
V
I
(t
1
)V
I
(t
2
)dt
2
dt
1
. . . +
1
(i)
n
_
t
t
o
V
I
(t
1
)dt
1
_
t
1
t
o
V
I
(t
2
)dt
2
. . .
_
t
n1
t
o
V
I
(t
n1
)dt
n1
2.3.4. Determinacion de probabilidades de transicion
Si estudiamos un sistema sujeto a una interaccion que depende en
forma explcita del tiempo entonces su evolucion se puede escribir como:
[(t))
I
= U
I
(t, t
o
) [(0)) =
m
_
(0)
m
U
I
(t, t
o
) [(0))
(0)
m
_
,
recuerde que
[(t))
I
= exp(iH
o
t/) [(t))
S
,
por lo tanto
_
(0)
m
U
I
(t, t
o
) [(0)) =
_
(0)
m
exp(iH
o
t/)U
S
(t, t
o
) exp(iH
o
t/) [(0))
=
n
_
(0)
m
U
S
(t, t
o
)
(0)
n
_
exp(i (E
m
E
n
) t/).
Consideremos el caso en que el sistema a t = 0 se encuentra en un estado
propio del sistema no perturbado, es decir, [(0)) =
(0)
n
_
, por lo tanto la
probabilidad de transicion del estado n al m esta dada por:
P
nm
=
(0)
m
U
I
(t, t
o
)
(0)
n
_
2
.
32
Para determinar esta cantidad emplearemos la serie de Dyson:
c
(0)
m
(t) =
mn
c
(1)
m
(t) =
1
i
_
t
t
o
_
(0)
m
V
I
(t)
(0)
n
_
dt
1
c
(2)
m
(t) =
1
(i)
2
k
_
t
t
o
_
(0)
m
V
I
(t)
(0)
k
_
dt
1
_
t
1
t
o
_
(0)
k
V
I
(t
2
)
(0)
n
_
dt
2
.
.
.
c
(0)
m
(t) =
mn
c
(1)
m
(t) =
1
i
_
t
t
o
exp(i (E
m
E
n
) t
1
/)
_
(0)
m
V (t
1
)
(0)
n
_
dt
1
c
(2)
m
(t) =
1
(i)
2
k
_
t
t
o
e
i(E
m
E
k
)t
1
/
_
(0)
m
V (t)
(0)
k
_
dt
1
_
t
1
t
o
e
i(E
k
E
n
)t
2
/
_
(0)
k
V (t
2
)
(0)
n
_
dt
2
.
.
.
Problema 10. Perturbacion constante: En este caso nos referimos como per-
turbacion a una interaccion que se conecta a un tiempo t = 0 y esta dada
por
V (t) =
_
0
V
o
t 0
t > 0
.
Observe que V
o
en esta expresion es un operador, por lo tanto las transicio-
nes inducidas por esta perturbacion se producen solo cuando los elementos
matriciales de V
o
sean no nulos. Para esto asumiremos que la base no pertur-
bada es [
n
) y el sistema se encuentra inicialmente en el estado [
N
), por
lo tanto:
c
N
(t) = 1 +
1
i
V
NN
t,
c
n
(t) =
V
nN
E
n
E
N
_
1 e
i
nN
t
_
.
La primera expresion indica como se modica el estado original, es decir,
[
N
)
_
1
i
V
NN
t
_
e
i
E
N
t
[
N
) e
i
(E
N
+V
NN
)t
[
N
) .
Esto corresponde, esencialmente, a un cambio de fase inducido por el corri-
miento de energa. La probabilidad de transcion a otro de los estados del
sistema es:
[c
n
(t)[
2
=
4 [V
nN
[
2
(E
n
E
N
)
2
sin
2
_
nN
2
t
_
,
nN
=
E
n
E
N
.
= [c
n
(t)[
2
max
sin
2
_
nN
2
t
_
, [c
n
(t)[
2
max
=
4 [V
nN
[
2
(E
n
E
N
)
2
.
Esta probabilidad adquiere valores signicativos solo para aquellos estados
con energas muy proximas a la energa del estado inicial. Esto requiere una
conservacion aproximada de la energa. Ademas, esta probabilidad es oscila-
toria y cambia signicativamente en tiempo que satisfacen
t =

E
n
E
N
tE .
Esto corresponde a decir que los unicos estados que estaran excitados por un
tiempo apreciable seran aquellos que aproximadamente conservan energa.
Problema 11. Espectro continuo: Una aplicacion interesante del ejemplo an-
terior es cuando el espectro del sistema es quasi continuo. Esto equivale
a que el espectro continua siendo discreto, pero sus niveles se encuentran
sucientemente cerca como para aproximar el sistema al continuo. Una apli-
cacion de esto es el problema de los modos de normales de una cavidad, al
imponer condiciones de borde perodicas y extender una de sus paredes hasta
el innito. Este sistema es bastante realista si consideramos una cavidad -
nita de largo l con una de sus paredes que transmite parcialmente emitiendo
luz, podemos imaginar el exterior de la cavidad como una segunda cavidad de
largo L , por lo tanto podemos describir el acoplamiento de la cavidad
original con el exterior.
Al tener un continuo lo relevante no es preguntarse por la probabilidad
de transicion a un estado dado, pues tiene medida nula. Por lo tanto debemos
preguntarnos por la probabilidad de transicion total a otros estados:
_
dn =
_
(E
n
)dE
n
.
Al mismo tiempo:
_
[c
n
(t)[
2
dE
n
=
4 [V
nN
[
2
2
_
sin
2
_
nN
2
t
_
2
nN
dE
n
=
2 [V
nN
[
2
t
[c
n
(t)[
2
=
2 [V
nN
[
2
t
(E
n
E
N
) .
34
Luego, la probabilidad de transicion total a un estado F con energa E
F
E
N
esta dada por:
_
[c
n
(t)[
2
(E
n
)dE
n
=
2t
_
[V
nN
[
2
(E
n
E
F
) (E
n
)dE
n
=
2t
[V
FN
[
2
(E
F
).
Como esta cantidad incrementa linealmente con el tiempo lo relevante es
determinar su tasa de cambio:
d
dt
__
[c
n
(t)[
2
(E
n
)dE
n
_
= w
FN
=
2
[V
FN
[
2
(E
F
).
Esta expresion es conocida como la regla de oro de Fermi.
Problema 12. Aproximacion de Born: Empleando el resultado anterior con-
sideraremos el problema de scaterring como un problema de teorade pertur-
baciones con el potencial siendo la perturbacion:
H
o
=
1
2m
P
2
.
Para construir un problema de espectro discreto se puede asumir que la
partcula se encuentra en una caja c ubica de grandes dimensiones, siendo L el
largo de cada uno de sus lados, por lo tanto los estados propios de H
o
son [k)
con r[k) =
1
L
3
exp(ik r). Por lo tanto el proceso de scattering es descrito
por el cambio [k) [k
). As los estados nales que debemos considerar son

aquellos que se encuentren en un angulo solido d en torno de la direccion
nal k
cuya magnitud [k
[ = k
adopta valores entre (k
, k
+dk
). El n ume-
ro total de estados que satisface esta condicion esta dado por el n umero de
puntos (n
x
, n
y
, n
z
) contenidos en el volumen n
2
dnd, con n
2
= n
2
x
+n
2
y
+n
2
z
.
Recuerde que la condicion de cuantizacion es
k
j
L = 2n
j
, j = x, y, z.
Empleando estas condiciones se tiene que:
(E
)dE
= n
2
dnd =
k
2
L
3
8
3
dk
d,
como E
=

2
2m
k
2
dE
=

2
k
m
dk
. Luego,
(E
) =
mL
3
2
k
d
8
3
.
As encontramos que:
w
kk
=
2mL
3
3
k
[k
[ V (r) [k)[
2
d
8
3
=
2m
3
L
3
k
_
exp(ik
r) V (r) exp(ik r) d
3
r
2
d
8
3
.
Como se tiene una partcula al interior de la caja moviendose con velocidad
v =
k
m
, el ujo incidente es S
i
= v/L
3
por unidad de area (cm
2
) y por unidad
de tiempo (s)
w
kk
= S
i
(k, k
)d =
k
mL
3
(k, k
)d
(k, k
) =
1
4
2m
2
_
V (r) exp(i (k k
) r) d
3
r
2
.
Como veremos mas adelante esta expresion es conocida como aproximacion
de Born a la amplitde de scattering.
Problema 13. Perturbacion perodica: Corresponde al caso de transciones ra-
diativas inducidas por un potencial del tipo
V (t) = V e
it
+V
e
it
,
donde V es el valor del operador al tiempo t = 0. La evolucion de los estados
diferentes al inicial se describe mediante:
c
n
(t) =
V
nN
+E
n
E
N
_
1 e
i(+
nN
)t
_
+
V
Nn
(E
n
E
N
)
_
1 e
i(
nN
)t
_
.
Empleando la discursion realizada en el primer ejemplo de esta seccion, la
condicion de cuasi conservacion de energa sera:
E
N
=
_
E
F
+
E
F

emision
absorcion
En el primero de estos, emision (+) el segundo termino en la ecuacion para
c
n
(t) se cancela, de la misma en la absorcion (-) se cancela el primer termino.
En forma compacta podemos escribir:
c
()
n
(t)
2
=
4
V
()
nN
2
( E
n
E
N
)
2
sin
2
_
(
nN
)
2
t
_
,
donde V
(+)
= V y V
(+)
= V
. Por lo tanto las probabilidades de transicion

son:
w
()
FN
=
2
V
()
FN
2
( E
N
E
F
).
Estas expresiones permiten justicar adecuadamente las suposiciones origi-
nales de Einstein al explicar la radiacion de cuerpo negro (conocidos como
coecientes de Einstein). Volveremos a este punto al estudiar interaccion ra-
diacion-materia.
36
Problema 14. Considere una partcula en una dimension moviendose bajo
la inuencia de alg un potencial independiente del tiempo. Los niveles de
energa y las correspondientes eigenfunciones para este problema se asumen
conocidas. Ahora, sometemos a la particula a un pulso viajero representado
por un potencial dependiente del tiempo:
V (t) = A(x ct)
1. Suponga que es sabido que en t = la partcula esta en el estado
fundamental de energa x[i) = u
i
(x). Obtenga la probabilidad de en-
contrar el sistema en alg un estado excitado con eigenfuncion de energa
x[f) = u
f
(x) en t = .
2. Interprete su resultado anterior fsicamente considerando el pulso de fun-
ciones delta como una superposicion de perturbaciones armonicas; re-
cuerde que:
(x ct) =
1
2c
_

de
i(
x
c
t)
Enfatice el rol de la conservacion de la energa, la cual es valida a un
cuando la perturbacion haya estado presente por un tiempo muy largo.
Solucion
Para este problema, se nos pide calcular la probabilidad de encontrar
al sistema en algun estado nal, conociendo el estado inicial de la partcula
sometida a un potencial de la forma:
V (t) = A(x ct)
Luego, utilizando la serie de Dyson, se tiene que:
c
(0)
n
(t) =
n,0
c
(1)
m
(t) =
1
i
_
t
t
0
(0)
m
[V
I
(t
1
)[
(0)
0
)dt
1
De este modo, el coeciente de orden 1 es:
c
(1)
m
=
1
i
_
t
t
0
m[V
I
(t
1
)[0)dt
1
=
1
i
_

m[e
iH
0
t
1
V (t
1
)e
iH
0
t
1
[0)dt
1
=
1
i
_

e
i(E
m
E
0
)t
1
m[V (t
1
)[0)dt
1
=
1
i
_

e
i
m0
t
1
m[
__

A(x ct)[x)x[dx
_
[0)dt
1
=
A
i
_

e
i
m0
t
1
_

(x ct)m[x)x[0)dxdt
1
=
A
i
_

e
i
m0
t
1
_

(x ct)u
m
(x)u
0
(x)dxdt
1
Antes de realizar la integracion espacial efectuareamos la integracion tempo-
ral con la ayuda de la transformada de Fourier de la funcion delta, con lo que
se obtiene:
_

e
i
m0
t
(x ct
) dt
=
1
2c
_

e
i
m0
t
dt
_
+
e
i
c
xit
d
=
1
c
_
+
e
i
c
x
d
1
2
_

e
i(
m0
)t
dt
=
1
c
_
+
e
i
c
x
(
m0
)d =
e
i
m0
c
x
c
.(2.13)
Por lo tanto, se obtiene:
c
(1)
m
=
A
ic
_

e
i
m0
c
x
u
m
(x)u
0
(x)dx (2.14)
Por ejemplo, para una partcula en una caja de largo L sometida a un po-
tencial, es conocida la solucion:
u
n
(x, t) =
_
2
L
e
ik
n
x
Donde k
n
=
2
n
=
2(n+2)
L
. Luego se tiene que:
c
(1)
m
=
2A
icL
_
L/2
L/2
exp [i(k
m0
+k
m
k
0
)x] dx
=
4A
icL
e
i(k
m0
+k
m
k
0
)L/2
e
i(k
m0
+k
m
k
0
)L/2
2i(k
m0
+k
m
k
0
)
=
2A
ic
sin((k
m0
+k
m
k
0
)L/2)
(k
m0
+k
m
k
0
)L/2
38
donde k
m0
=
m0
/c. La probabilidad de transicion desde el nivel fundamental
a un estado arbitrario es
P
0m
=
c
(1)
m
2
=
4A
2
2
c
2
_
sin((k
m0
+k
m
k
0
)L/2)
(k
m0
+k
m
k
0
)L/2
_
2
.
Los valores propios de la energa para una partcula en esta caja son:
E
n
=

2
k
2
n
2m
=
2
2
2
(n + 2)
2
mL
2
,
con lo que se puede simplicar la expresion anterior.
Para sumar todas las propbabilidades de transicion debemos normali-
zar adecuadamente el estado del sistema, despues de esto podemos analizar
la conservacion de energa sumando todos los terminos.
Problema 15. Un oscilador armonico unidimensional esta en su estado fun-
damental para t < 0. Para t 0 es sujeto a una fuerza (!no potencial)
dependiente del tiempo pero espacialmente uniforme:
F(t) = F
0
e
t/
1. Usando teora de perturbaciones dependientes del tiempo a primer orden,
obtener la probabilidad de encontrar el oscilador en su primer estado
excitado para t > 0. Muestre que el limite t ( nito) de su expresion
es independiente del tiempo. Es esto razonable o sorprendente?
2. Podemos encontrar estados excitados superiores?.
Ayuda: Puede usar:
n
[X[n) =
_

2m
_
n
n
,n1
+
n + 1
n
,n+1
_
(2.15)
Solucion: El termino forzante, puede ser escrito como un potencial recor-
dando que para fuerzas conservativas, se tiene que:
F =
dV
dx
De este modo se tiene que:
F = F
0
e
t/
=
dV
dx
V (x, t) =
_
F
0
e
t/
dx
V (x, t) = F
0
e
t/
x
De este modo, el Hamiltoniano del oscilador lo escribimos ahora como
(utilizando la notacion de operadores de creacion y destruccion):
H = H
0
+V
=
_
a
a +
1
2
_
+F
0
e
t/
x
Luego, utilizando las series de Dyson se tiene que:
c
(0)
n
=
n0
c
(1)
n
=
1
i
_
t
t
0
0
n
[V
I
(t
1
)[
0
0
)dt
1
=
1
i
_
t
0
n[e
iH
0
t
1
F
0
e
t
1
/
xe
iH
0
t
1
[0)dt
1
=
F
0
i
_
t
0
e
t
1
/
e
i
n0
t
1
n[x[0)dt
1
Donde
n0
=
E
n
E
0
. Utilizando la ecuacion 2.15 se tiene que:

c
(1)
n
=
F
0
i
_
t
0
e
(1/+i
n0
)t
1
_

2m
n,1
dt
1
=
F
0
i
_

2m
n,1
_
t
0
e
(1/+i
n0
)t
1
dt
1
=
F
0
i
_

2m
n,1
_
e
(1/+i
n0
)t
1
1/ +i
n0
c
n
=
nn
+
F
0
i
_

2m
n,1
_
e
(1/+i
n0
)t
1
1/ +i
n0
La probabilidad de que el sistema pase de el estado fundamental [0) a
uno [n) se obtiene como modulo cuadrado del coeciente c
n
, es decir:
P
0m
=
c
(1)
m
2
=
_
F
0
_

2m
_
2
n,1
_
1
1/
2
+
2
10
_
_
2e
t/
[cosh(t/) cos(
10
t)]
_
Esta probabilidad depende fuertemente de la amplitud de la fuerza, F
0
. El
termino e
t/
cosh(t/) converge rapidamente a 1/2 cuando t . Sin
embargo, el termino e
t/
cos(
10
t), en el lmite cuando t converge a
cero, luego la probabilidad para un tiempo muy grande es independiente del
tiempo y tiene un valor nito dado por:
P
0n
= c
n
c
n
=
_
F
0
2
_

2m
_
2
n,1
_
1
1/
2
+
2
10
_
.
40
Por lo tanto, dado que el acoplamiento con la fuerza es lineal, al encontrarse
el sistema en el estado fundamental solo puede alcanzar el primer estado
excitado (a primer orden en teora de perturbaciones). As se tiene que:
P
01
= c
n
c
n
=
_
F
0
2
_

2m
_
2 _
1
1/
2
+
2
_
.
Una graca de esta probabilidad P
10
, se muestra a continuacion en funcion
del tiempo: Vemos que este termino converge a un valor independiente del
0 10 20 30 40 50 60
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
Grafica de probabilidad de salto entre los estados0 y 1 versus tiempo
Tiempo
Probabilidad
Figura2.1. Graco de P
10
(probabilidad de salto del estado |0 al |1), en funcion
del tiempo
tiempo dado por:
P
01
=
F
2
0
2m
_
1
1/
2
+
2
10
_
Este resultado no debiera sorprendernos ya que es de esperar que el oscila-
dor tenga probabilidad de cambiar de estado al ser sometido a una fuerza
periodica.
2.4. Metodo variacional
Los metodos anteriores nos permiten obtener las soluciones perturba-
das en los casos en que se conoce la solucion a un problema similar al que se
desea resolver. Sin embargo, esto no siempre es as y es conveniente obtener
la energa del estado fundamental empleando estados (funciones) de prueba.
Para esto asumiremos que [) es un estado arbitrario normalizado, que
pertenece al espacio de estados del sistema. El valor medio del Hamiltoniano
en este estado es:
H) = [ H[)
0
[ H[
0
) = E
0
,
donde [
0
) es el estado fundamental. Para demostrar esto consideraremos un
estado no normalizado, que corresponde a una peque na desviacion del estado
fundamental
[) = [
0
) + [) ,
donde 1 y [) es un estado normalizado. Por lo tanto el valor medio de
H en este estado perturbado es:
H) =
[ H[)
[)
= E
[) = (
0
[ + [) ([
0
) + [)) = 1 +
_
0
[) +[
0
)) +
2
_
[ H[) = (
0
[ + [) H ([
0
) + [)) = E
0
+
_
0
[H[) +[H[
0
)) +
2
[H[)
_
.
Ademas, se asume que E = E
0
+E
(1)
0
+
2
E
(2)
0
, podemos igualar potencias
en el parametro peque no. El resultado es:
_
E
0
+E
(1)
0
+
2
E
(2)
0
_
_
1 +
_
(
0
[) +[
0
)) +
2
__
= E
0
+
_
(
0
[H[) +[H[
0
)) +
2
[H[)
_
E
(1)
0
+E
0
(
0
[) +[
0
)) = (
0
[H[) +[H[
0
))
= E
0
(
0
[) +[
0
)) E
(1)
0
= 0.
E
(2)
0
+E
0
= [H[)
E
(2)
0
= [H E
0
[) ,
como [) se puede expandir en la base de estados propios [
n
) de H como:
[) =
n
[
n
)
[H E
0
[) =
n=0
[
n
[ (E
n
E
0
)
2
0.
As, hemos obtenido que con un error de orden en el estado se obtiene un
error
2
en la energa del estado fundamental:
E = E
0
+
2
[H[) .
42
Problema 16. Energa deDeuteron: La fuerza entre el proton y el neutron
puede ser descrito por un potencial tipo Yukawa
V (r) = V
o
exp(r/r
o
)
r/r
o
,
donde r
o
representa el alcance de la fuerza y para los mesones de masa
esta dado por:
=

c
.
El Hamiltoniano para el estado del sistema es:
H =
1
2m
P
2
+V (r),
donde m es la masa reducida de un nucleon. La energa del estado funda-
mental corresponde a la energa de ligazon. La funcion de prueba a emplear
sera:
(r) = exp(r/r
o
),
donde es el parametro variacional. Al obtener una expresion para E
o
=
E
o
() minimizaremos sobre para obtener la energa aproximada del estado
fundamental (estado-s). La energa con la funcion de onda elegida es:
E =
2
2mr
2
o
_
4
2
(2 + 1)
2
1
_
, =
2mV
o
r
2
o
2
.
Problema 17. Demostrar que en el ejemplo anterior E esta dado por la expre-
sion obtenida. Determine el valor de que minimiza el valor de E. Finalmente,
muestre que si asume que E
0
= 2,226 Mev se obtiene que:
V
0
r
2
0
= 41,37
2
Mev(10
26
cm
2
)
Solucion
Recuerde que en coordenadas esfericas
2
=
1
r
2
r
r
2
r
+
1
r
2
sin
sin

+
1
r
2
sin
2
2
y
_
0
x
n
e
x
= n!. En este caso la funcion de prueba no es nor-
malizada, as es que determinamos esta constante: x = 2r/r
o
_

0
[(r)[
2
d
3
r = 4
_

0
e
2r/r
o
r
2
dr
=
r
3
o
2
3
_

0
e
x
x
2
dx =
r
3
o
3
.
Ademas,
_
(r)
1
2m
p
2
(r)d
3
r =

2
2m
_
(r)
2
(r)d
3
r
=
2
2
m
_

0
(r)
1
r
2
r
_
r
2
(r)
r
_
r
2
dr
=
2
2
m
_

0
(r)

r
_
r
2
r
o
(r)
_
dr
=
2
2
m
_

0
e
2
r
o
r
_
2r
r
o

2
r
2
r
2
o
_
dr
=

2
r
o
m
_

0
e
x
_
x
x
2
4
_
dx =

2
r
o
2m
_
(r)V (r)(r)d
3
r = 4V
o
_
(r)
exp(r/r
o
)
r/r
o
(r)r
2
dr
= 4V
o
r
o
_
e
(2+1)
r
o
r
rdr
=
4V
o
r
3
o
(2 + 1)
2
_
e
x
xdx = V
o
4r
3
o
(2 + 1)
2
Luego
E() = [ H[) / [ [)
E() =
2
r
o
2m
V
o
4r
3
o
(2+1)
2
r
3
o
3
,
2
=
2mV
o
r
2
o
2
E() =

2
2mr
2
o
3
_
1

2mr
2
o
V
o
2
4
(2 + 1)
2
_
E() =

2
2mr
2
o
_
4
2
3
(2 + 1)
2
_
dE()
d
=
_
2
12
2
2
(2 + 1)
2
+
16
2
3
(2 + 1)
3
_
= 0

2
6
(2 + 1)
2

8
(2 + 1)
3
= 1

2
=
(2 + 1)
3
2(2 + 3)
.
44
Reemplazando en la energa obtenida previamente:
E() =

2
2mr
2
o
2
(2 + 3)
(2 + 3 4 2)
=

2
2mr
2
o
2
(2 1)
(2 + 3)
.
Como se indica que experimentalmente se conoce que E
o
= 2,226 Mev, por
lo tanto:
E
o
o
=
2
(2 1)
(2 + 3)
V
o
o
=
(2 + 1)
3
2(2 + 3)
E
o
V
o
=
2
3
(2 1) (2 + 3)
(2 + 1)
4
Problema 18. Estime el menor valor propio () de la ecuacion diferencial
d
2
dx
2
+ ( [x[) = 0, (x ) = 0
usando el metodo variacional, con
(x) =
_
c ( [x[) , [x[ <
0, [x[ >
como funcion de prueba (Cuidado: d/dx es discontinua en x = 0). Puede
demostrarse que el valor exacto del menor eigenvalor es 1.019.
Solucion: Nuevamente la funcion de prueba no esta normalizada, por lo
tanto se requiere determinar:
[ [) =
_

(x)dx = c
2
_

( [x[)
2
dx = 2c
2
_

0
( x)
2
dx
= 2c
2
= 2c
2
3
3
c =
_
3
2
3
.
Reescribiremos el problema como:
2
2m
d
2
dx
2
+

2
2m
[x[ =
=

2
2m
H =
1
2m
p
2
+

2
2m
[x[
E() =

_
=
3
2
3
__

( [x[)
_
2
2m
d
2
dx
2
+

2
2m
[x[
_
( [x[) dx
_
=
3
2
3
_
_
_
0
( +x)
_
2
2m
d
2
dx
2

2
2m
x
_
( +x) dx
_
0
( x)
_
2
2m
d
2
dx
2
+

2
2m
x
_
( x) dx
_
_
Antes de continuar, debemos decir que hay que ser cuidados con el termino
relacionado con la segunda derivada. Como se menciono en el enunciado la
derivada de la funcion es discontinua, por lo tanto escribiremos:
_

(x)
d
2
(x)
dx
2
dx = lm
0
_

(x)
d
2
(x)
dx
2
dx = lm
0
_

(x)d
_
d(x)
dx
_
= lm
0
_
(x)d
_
d(x)
dx
_
_
d(x)
dx
_
2
dx
_
= lm
0
_
2c
2
( ) 2c
2
_
= 2
3
2
3
=
3
2
.
Con este resultado podemos reescribir el valor medio de la energa como:
E() =

2
m
3
2
+
3
2
3
2
m
_

0
_
2
x 2x
2
+x
3
_
dx
=

2
m
_
3
2
+

8
_
.
El mnimo se obtiene en forma numerica como
= 3,634. En forma alter-

nativa podemos minimizar la energa con respecto al parametro ,
dE()
d
=

2
m
_
3
+
1
8
_
= 0
3
+
1
8
= 0
= 3,6342
Al reemplazar
en la energa se obtiene que:

E( =
) =

2
m
0. 681 42 =

2
2m
1,3628.
Por lo tanto el valor obtenido para = 1,3628. Esto se puede mejorar en
forma notoria cambiando la funcion de prueba por una tipo exponencial.
3. Partculas identicas
3.1. Introduccion
3.2. Simetra de permutaci on
En mecanica clasica siempre es posible distinguir las trayectorias de
dos partculas:
Lo que usualmente se puede hacer en mecanica clasica es imaginar que
cada una de las partculas puede ser coloreada en forma diferente para
distinguirlas adecuadamente.
Una situacion distinta sucede con las partculas identicas en Mecanica
Cuantica, ya que estas no pueden ser rotuladas ni coloreadas, pues el
identicarlas completamente en su trayectoria corresponde a realizar per-
manentemente mediciones de posicion que nos hacen perder informacion
de su momentum.
Consideremos la gura que muestra dos fotones que llegan sim ultanea-
mente a un separador de haces
Las dos situaciones presentadas no pueden ser distinguidas, aun en el caso
que llegue un foton al separador de haces, pues podra ser uno transmitido
en el separador de haces o uno reejado en el.
En el caso mas simple tendremos solo dos partculas identicas
[k
) , [k
),
donde k
j
denota al conjunto de n umeros cuanticos que caracterizan a la
partcula.
El estado del sistema de dos partculas puede ser escrito como:
[k
) [k
) o [k
) [k
) .
3. Partculas identicas 47
0 10 20 30 40 50 60
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
Grafica de probabilidad de salto entre los estados0 y 1 versus tiempo
Tiempo
Probabilidad
Figura3.1. Realizar graco adecuado
Al efectuar mediciones sobre el sistema de partculas se obtendran k
y k
. Sin
embargo, este resultado se obtiene en ambos estados o en una superposicion
arbitraria de ellos
c [k
) [k
) +d [k
) [k
) .
Como vemos estos nos lleva a un nuevo tipo de degeneracion llamada degene-
racion de intercambio, es decir, que el CCOC no determina completamente
el estado de un sistema de partculas identicas.
Observemos que sucede con la operacion de intercambio de partculas,
la cual se dene mediante:
P
12
[k
) [k
) = [k
) [k
)
P
12
[k
) [k
) = [k
) [k
) P
2
12
= 1, P
12
= P
21
donde P
12
es el operador de intercambio entre las partculas 1 y 2. Veamos
que sucede con la accion de este operador sobre los operadores. Para lo cual
consideremos la accion de uno de los observables:
A
1
[a
) [a
) = a
[a
) [a
) ,
A
2
[a
) [a
) = a
[a
) [a
) ,
A
1
= A
1
1
2
, A
2
= 1
1
A
2
.
Aplicando el operador de intercambio se tiene que:
P
12
A
1
P
1
12
P
12
[a
) [a
) = a
P
12
[a
) [a
)
P
12
A
1
P
1
12
[a
) [a
) = a
[a
) [a
)
P
12
A
1
P
1
12
= A
2
,
48
por lo tanto el operador de intercambio modica los rotulos de las partculas
en los operadores.
Si deseamos denir estados propios de P
12
consideremos las siguientes
superposiciones:
[k
=
1
2
1/2
([k
) [k
) [k
) [k
)) ,
las cuales son conocidas como superposiciones simetricas (+) y antisimetricas
(-) que corresponden a valores propios 1 para P
12
. En forma natural se
pueden denir operadores de simetrizacion y antisimetrizacion como:
S
12

1
2
(1 +P
12
) , A
12

1
2
(1 P
12
) .
La accion de estos operadores sobre una superposicion arbitraria como la
escrita anteriomente nos da:
S
12
(c [k
) [k
) +d [k
) [k
)) =
c +d
2
([k
) [k
) +[k
) [k
)) ,
A
12
(c [k
) [k
) +d [k
) [k
)) =
c d
2
([k
) [k
) [k
) [k
)) .
Lo que hemos denido para el caso de dos partculas se puede genera-
lizar faclmente al caso de un sistema de muchas partculas identicas:
P
jk
[k
) [k
) . . .
k
j
_
. . .
k
k
_
. . .
k
l
_
. . . = [k
) [k
) . . .
k
k
_
. . .
k
j
_
. . .
k
l
_
. . .
P
2
jk
= 1,
y los valores propios de P
jk
son 1. Ademas, en general se cumple que:
[P
ij
, P
kl
] ,= 0
Problema 19. Considere un sistema de dos partculas descritas por el Hamil-
toniano
H =
p
2
1
2m
+
p
2
2
2m
+V ([x
1
x
2
[) +V
ext
(x
1
) +V
ext
(x
2
).
Que sucede con el sistema si inicialmente se encuentra en un estado simetrico
o antisimetrico?
Para responder determinaremos como se modica el Hamiltoniano
frente a la accion de P
12
(como operador):
P
12
HP
1
12
=
p
2
2
2m
+
p
2
1
2m
+V ([x
2
x
1
[) +V
ext
(x
2
) +V
ext
(x
1
)
=
p
2
1
2m
+
p
2
2
2m
+V ([x
1
x
2
[) +V
ext
(x
1
) +V
ext
(x
2
) = H
P
12
H HP
12
= 0, [H, P
12
] = 0.
Por lo tanto, si el sistema se encuentra incialmente en un estado simetrico o
antisimetrico lo mantendra durante la interaccion.
Problema 20. Considere un sistema de tres partculas identicas y determine
los estados simetricos y antisimetricos de ellas.
En este caso tendremos un estado del tipo
[k
) [k
) [k
) .
Observemos que existe degeneracion ante el intercambio de pares de partcu-
las, es decir, existen los operadores de intercambio P
12
, P
13
, P
23
y las su-
peropsiciones deben ser completamente simetricas o antisimetricas frente al
intercambio de partculas. De esta forma los estados resultantes son:
[k
=
1
(3!)
1/2
_
_
([k
) [k
) [k
) [k
) [k
) [k
)) +
([k
) [k
) [k
) [k
) [k
) [k
)) +
([k
) [k
) [k
) [k
) [k
) [k
))
_
_
.
Donde hemos asumido implcitamente que los n umeros cuanticos son todos
diferentes. Si se repite uno o mas de ellos no es posible construir superposi-
ciones antisimetricas. Consideremos el caso en que los dos primeros ndices
son iguales:
[k
)
+
=
_
2!
3!
_
1/2
[[k
) [k
) [k
) +[k
) [k
) [k
) +[k
) [k
) [k
)] .
En el caso general el coeciente de normalizacion es:
_
N
1
!N
2
! . . . N
n
!
N
_
1/2
, N
1
+N
2
+. . . +N
n
= N.
3.3. Postulado de simetrizacion
En la naturaleza los conjuntos de partculas identicas que se encuen-
tran en estados totalmente:
Simetricos ante el intercambio de cualquier par de partculas de
conocen como BOSONES, es decir, satisfacen la estadsitica de
Bose-Einstein.
Antisimetricos se conocen como fermiones, es decir, satisfacen
la estadstica de Fermi.
50
Hasta ahora no se ha observado una mezcla de simetras en un conjunto
de partculas identicas, aunque se ha intentado emplear para explicar algunos
fenemenos como superconductividad.
P
ij
[N bosones identicos) = +[N bosones identicos) ,
P
ij
[N fermiones identicos) = [N fermiones identicos) .
Existe una conexion entre el spin de la partcula y la estadstica que obedece:
Partculas de spin entero son bosones (
,
0
,
4
He)
Partculas de spin semientero son fermiones (e
, p
+
,
3
He)
Esta conexion entre spin y estadstica solo puede ser derivada en Mecani-
ca Cuantica Relativista. En terminos las clasicos la partculas satisfacen la
estadstica de Maxwell-Botzman, en que las partculas se pueden reconocer
en forma individual
Una consecuencia directa de la relacion entre spin y estadstica es el
Principio de exclusion de Pauli, que dice que dos electrones no pueden
tener todos sus n umeros cuanticos iguales. Esto resulta directo ya que los
electrones poseen spin 1/2 y por lo tanto son fermiones y deben ser anti-
simeticos ante el intercambio de partculas, pero no es posible construir un
estado antisimetrico con dos partculas con sus n umeros cuanticos identicos.
Problema 21. Consideremos dos electrones caracterizados por valores pro-
pios k
y k
. En este caso los estados posibles en cada una de las estadsticas

son:
Bose-Einstein [k
) [k
) [k
) [k
)
1
2
1/2
([k
) [k
) +[k
) [k
)) 3
Fermi-Dirac
1
2
1/2
([k
) [k
) [k
) [k
)) 1
Maxwell-Boltzman [k
) [k
) [k
) [k
) [k
) [k
) [k
) [k
) 4.
Las di-
ferencias entre estos sistemas se hacen evidentes a bajas temperaturas.
3.4. Sistema de dos electrones
Consideremos un sistema de dos electrones, que puede ser completa-
mente caracterizado por sus posiciones x
1
, x
2
y por sus componentes de spin
m
1
y m
2
. Por lo tanto la funcion de onda del sistema se puede escribir como:
(x
1
, x
2
) =
m
1
,m
2
c(m
1
, m
2
) x
1
, m
1
; x
2
, m
2
[) .
Asumieremos tambien que:
_
S
2
, H
= 0, por lo tanto la funcion de onda

esta dada por:
(x
1
, x
2
) = (x
1
, x
2
),
donde es la funcion de spin y que debe ser:
=
_
++
s
=
1
2
1/2
_
+
+
+
_
a
=
1
2
1/2
_
+
+
+
_
S = 1
S = 0
,
en que + () indica que m
j
= +1/2 (1/2). La accion del operador de
intercambio sobre el estado del sistema nos da:
x
1
, m
1
; x
2
, m
2
[P
12
[) = x
2
, m
2
; x
1
, m
1
[)
x
1
, m
1
; x
2
, m
2
[) = x
2
, m
2
; x
1
, m
1
[) .
Observe, ademas por la forma en que se descompone la funcion de onda del
sistema:
P
12
= P
(espacial)
12
P
(spin)
12
(x
1
, x
2
)
P
(espacial)
12
(x
2
, x
1
) (m
1
, m
2
)
P
(spin)
12
(m
2
, m
1
).
P
(spin)
12
=
1
2
_
1 +
4
2
S
1
S
2
_
,
S
1
S
2
=

2
4
_
1
3
triplete
singlete
,
por lo tanto, como el estado debe ser completamente antisimetrico, si la parte
espacial es simetrica la correspondiente a spin es antisimetrica y viceversa.
Un ejemplo simple de la aplicacion del operador de intercambio es la
accion sobre un sistema en que el interaccion mutua V ([x
1
x
2
[)y de spin
S
1
S
2
se pueden despreciar, es decir, el sistema puede ser descrito mediante
el Hamiltoniano:
H =
p
2
1
2m
+
p
2
2
2m
+V
ext
(x
1
) +V
ext
(x
2
),
por lo tanto la funcion de onda espacial se puede escribir como el producto
de las respectivas funciones de onda
A
(x
1
) y
B
(x
2
) :
(x
1
, x
2
) =
1
2
1/2
(
A
(x
1
)
B
(x
2
)
A
(x
2
)
B
(x
1
))
+
S
,
T
.
De esta forma la probabilidad de encontrar a la partcula 1 en el volumen
d
3
x
1
y a la partcula 2 en el volumen d
3
x
2
esta dada por:
P = [(x
1
, x
2
)[
2
d
3
x
1
d
3
x
2
=
1
2
_
[
A
(x
1
)[
2
[
B
(x
2
)[
2
+[
A
(x
2
)[
2
[
B
(x
1
)[
2
2 Re (
A
(x
2
)
B
(x
1
)
A
(x
2
)
B
(x
1
))
_
d
3
x
1
d
3
x
2
.
52
Observe que:
Si los electrones se encuentran en el estado singlete la probabilidad de
encontrarse en el mismo punto es distinta de cero.
Si los electrones se encuentran en un estado de triplete la probabilidad
de econtrarse en el mismo punto se cancela.
Los aspectos anteriores, de distinguibilidad, son relevantes si las funciones
de onda poseen un solapamiento apreciable. Si las partculas estan en una
region del espacio con solapamiento despreciable solo el primer termino
es relevante.
3.4.1. Estados de un sistema de tres electrones
3.5. Atomo de Helio
Un caso particular del tipo de Hamiltoniano que hemos discutido es
el atomo de Helio, que puede ser escrito como:
H (r
1
, r
2
) =
_
p
2
1
2m

Ze
2
r
1
_
+
_
p
2
2
2m

Ze
2
r
2
_
+
e
2
r
12
.
En este caso el problema no perturbado es el de cada uno de los electrones
junto a la repulsion coulombiana. Los operadors colectivos naturales de este
sistema son:
S
z
= s
z
1
+s
z
2
; L
z
= l
z
1
+l
z
2
; J
z
= L
z
+S
z
;
L = l
1
+l
2
; S = s
1
+s
2
; J = L+S.
Es facil demostrar que los operadores S
z
, L
z
, J
z
, L
2
, S
2
, J
2
conmutan entre
si y por lo tanto existe una ba com un en que H (r
1
, r
2
) es diagonal con valore
propios M
S
, M
L
, M, L(L+1), S(S+1), J(J +1), respectivamente. Algunos
comentarios:
H es degenerado en M
S
, M
L
, M.
H es degenerado en J, pues no depende del spin de los electrones y por lo
tanto es independiente de como se acople el momentum angular orbital
y el spin (acoplamiento spin-orbita).
Aunque H no dependa de S no podemos armar que los valores propios
de H son independientes de S, pues la funcion de onda total debe ser
antisimetrica.
Esta ultima armacion se puede demostrar de la siguiente forma (Hei-
senberg, 1926): Sea (r
1
, r
2
) una solucion aproximada de H (r
1
, r
2
) =
E (r
1
, r
2
), como el hamiltoniano es simetrico en r
1
y r
2
(r
2
, r
1
) tambien
es una solucion aproximada. Por lo tanto podemos construir las siguientes
superposiciones simetricas y antisimetricas:
s
=
1
(2(1 +))
1/2
((r
1
, r
2
) +(r
2
, r
1
))
a
=
1
(2(1 ))
1/2
((r
1
, r
2
) (r
2
, r
1
))
donde
=
_ _

(r
1
, r
2
)(r
2
, r
1
)d
3
r
1
d
3
r
2
.
Como se tiene dos electrones la funcion de onda total debe ser antisimetrica
y se puede construir como:
1
=
s

s
,
2
=
a

t
.
Con esta funciones podemos calcular la diferencia de energa:
E
1
E
2
=
_ _
(
s
(r
1
, r
2
)H
s
(r
2
, r
1
)
a
(r
1
, r
2
)H
a
(r
2
, r
1
)) d
3
r
1
d
3
r
2
=
1
2(1
2
)
_ _ _
(1 ) (
(r
1
, r
2
) +
(r
2
, r
1
)) H ((r
1
, r
2
) +(r
2
, r
1
))
(1 +) (
(r
1
, r
2
)
(r
2
, r
1
)) H ((r
1
, r
2
) (r
2
, r
1
))
_
d
3
r
1
d
3
r
2
=
1
(1
2
)
_ _ _
(
(r
2
, r
1
)H(r
1
, r
2
) +
(r
1
, r
2
)H(r
2
, r
1
))
(
(r
1
, r
2
)H(r
1
, r
2
) +
(r
2
, r
1
)H(r
2
, r
1
))
_
d
3
r
1
d
3
r
2
=
2 Re
(1
2
)
_ _
[
(r
1
, r
2
)H(r
2
, r
1
)
(r
1
, r
2
)H(r
1
, r
2
)] d
3
r
1
d
3
r
2
Es posible resolver este problema mediante diferentes esquemas: (i)
Teora de perturbaciones, (ii) metodo variacional, (iii) metodo de Hartree-
Fock y (iv) metodo de Hylleraas. Cada uno de los metodos nos permite un
grado diferente de aproximacion a la solucion para el caso de Z peque no. El
metodo de Hartree-Fock es tambien un metodo variacional, solo que en el no
se asuma una dependencia explcita de las funciones de onda facotrizadas y
se efectuan variaciones en la funcionas de onda mismas, con lo que resulta
un sistema de ecuaciones integro-diferenciales que son resueltos en forma
numerica. Sin embargo el hecho de emplear funciones de onda factorizadas
como funciones de prueba conlleva un inconveniente, pues pondera de igual
forma funcionas de onda que deben tener diferente probabilidad. Esto se
resuelve con el metodo de Hylleraas, quien introduce variables relativas y
as da cuenta de la correlacion de los electrones.
54
3.5.1. Teora de perturbaciones
Si despreciamos la interaccion entre electrones y como no hemos consi-
derado interaccion spin-spin la funcion de onda del sistema debe ser el de dos
atomos hidrogenoides con Z = 2. y la perturbacion la interaccion coulombia-
na entre los electrones, es decir,
H = H
o
1
+H
o
2
+V
12
H
o
j
=
_
p
2
j
2m

Ze
2
r
j
_
, V
12
=
e
2
r
12
.
Las soluciones al problema no perturbado son:
(0)
(r
1
, r
2
) =
n
1
(r
1
)
n
2
(r
2
),
n
1
(r
1
) = R
n
1
,l
1
(r
1
)Y
l
1
,m
1
(
1
,
1
)
H
o
j
(0)
(r
1
, r
2
) = (E
n
1
+E
n
2
)
(0)
(r
1
, r
2
)
E
n
j
=
mZ
2
e
4
2
2
n
2
1
2
.
donde n
j
n
j
l
j
m
j
. Por lo tanto las funciones de onda antisimetricas se
pueden escribir como:
1
=
1
2
_
n
1
(r
1
)
n
2
(r
2
) +
n
1
(r
2
)
n
2
(r
1
)
_

s
,
2
=
1
2
_
n
1
(r
1
)
n
2
(r
2
)
n
1
(r
2
)
n
2
(r
1
)
_

t
.
En este caso especco podemos determinar la diferencia de energa entre el
estado simetrico y el antisimetrico. Para se requiere determinar:
E
1
=
_ _

1
(r
1
, r
2
)H
1
(r
2
, r
1
)d
3
r
1
d
3
r
2
= (E
n
1
+E
n
2
) +n
1
; n
2
[
e
2
r
12
[n
1
; n
2
) +n
1
; n
2
[
e
2
r
12
[n
2
; n
1
)
E
2
=
_ _

2
(r
1
, r
2
)H
2
(r
2
, r
1
)d
3
r
1
d
3
r
2
= (E
n
1
+E
n
2
) +n
1
; n
2
[
e
2
r
12
[n
1
; n
2
) n
1
; n
2
[
e
2
r
12
[n
2
; n
1
)
Para determinar las energas en estas expresiones se requiere emplear la si-
guiente expansion:
1
r
12
=
1
(r
2
1
+r
2
2
2r
1r
2
cos )
1/2
=
l=0
r
l
<
r
l+1
>
4
2l + 1
l
m=l
Y
m
l
(
1
,
1
)Y
m
l
(
2
,
2
)
_
d
j
Y
m
l
(
j
,
j
) = 2
1/2
l,0
m,0
Por lo tanto,
n
1
; n
2
[
e
2
r
12
[n
1
; n
2
) =
l=0
_ _
r
2
1
dr
1
r
2
2
dr
2
R
2
n
1
,l
1
(r
1
)
r
l
<
r
l+1
>
R
2
n
2
,l
2
(r
2
)
4
2l + 1
m=l
__
d
1
Y
m
1
l
1
(
1
,
1
)Y
m
l
(
1
,
1
)Y
m
1
l
1
(
1
,
1
)
_
d
2
Y
m
2
l
2
(
2
,
2
)Y
m
l
(
2
,
2
)Y
m
2
l
2
(
2
,
2
)
_
=
l=0
_ _
r
2
1
dr
1
r
2
2
dr
2
R
2
n
1
,l
1
(r
1
)
r
l
<
r
l+1
>
R
2
n
2
,l
2
(r
1
)
4
2l + 1
__
d
1
Y
m
1
l
1
(
1
,
1
)Y
0
l
(
1
,
1
)Y
m
1
l
1
(
1
,
1
)
_
d
2
Y
m
2
l
2
(
2
,
2
)Y
0
l
(
2
,
2
)Y
m
2
l
2
(
2
,
2
)
_
.
Analogamente,
n
1
; n
2
[
e
2
r
12
[n
2
; n
1
) =
l=0
_ _
r
2
1
dr
1
r
2
2
dr
2
R
n
1
,l
1
(r
1
)R
n
2
,l
2
(r
1
)
r
l
<
r
l+1
>
R
2
n
2
,l
2
(r
1
)
4
2l + 1
__
d
1
Y
m
1
l
1
(
1
,
1
)Y
m
1
m
2
l
(
1
,
1
)Y
m
2
l
2
(
1
,
1
)
_
d
2
Y
m
2
l
2
(
2
,
2
)Y
m
1
m
2
l
(
2
,
2
)Y
m
1
l
1
(
2
,
2
)
_
.
l=0
_ _
r
2
1
dr
1
r
2
2
dr
2
R
n
1
,l
1
(r
1
)R
n
2
,l
2
(r
1
)
r
l
<
r
l+1
>
R
2
n
2
,l
2
(r
1
)
4
2l + 1
__
d
1
Y
m
1
l
1
(
1
,
1
)Y
m
1
m
2
l
(
1
,
1
)Y
m
2
l
2
(
1
,
1
)
_
2
Para esto recuerde que:
Y
m
l
(, ) e
im
Y
m
1
l
1
(
1
,
1
)Y
m
l
(
1
,
1
)Y
m
1
l
1
(
1
,
1
) e
i(mm
1
+m
1
)
m = 0
Y
m
1
l
1
(
1
,
1
)Y
m
l
(
1
,
1
)Y
m
2
l
2
(
1
,
1
) e
i(mm
1
+m
2
)
m = m
1
m
2
A modo de ejemplo considere que ambos electrones se encuentran en el estado
fundamental, es decir, la funcion de onda para ambos es
100
(x
j
) o el estado
(1s)
2
. Como este estado es necesariamente simetrico, pues l = 0, la funcion
de spin debe ser el singlete y la funcion total esta dada por:
(x
1
, x
2
) =
100
(x
1
)
100
(x
2
)
singlete
=
8
a
3
o
exp(2(r
1
+r
2
)/a
o
)
singlete
.
La energa de este estado es:
E = 8
e
2
2a
o
1,3E
exp
, E
exp
= 78,8eV
56
Para obtener una correccion a este valor de la energa podemos emplear teora
de perturbaciones, con lo cual
(1s)
2 =
_
e
2
r
12
_
=
64
2
a
6
o
e
2
_ _
exp(4(r
1
+r
2
)/a
o
)
1
r
12
d
3
r
1
d
3
r
2
.
_ _
e
(4(r
1
+r
2
)/a
o
)
1
r
12
dr
1
dr
2
=
_
_
_
r
1
0
e
(4(r
1
+r
2
)/a
o
)
r
2
2
r
1
dr
2
+
_
e
(4(r
1
+r
2
)/a
o
)
r
2
2
r
2
dr
2
_
r
2
1
dr
1
=
5
128
a
5
o
2
5

(1s)
2 =
5
2
e
2
2a
o
, E +
(1s)
2 = 74,8eV.
Si se asume que uno de los electrones se encuentra en el estado fundamental
y el otro en un estado excitado, la funcion de onda se puede escribir como:
(x
1
, x
2
) =
1
2
1/2
(
100
(x
1
)
nlm
(x
2
)
nlm
(x
2
)
100
(x
1
)) .
En este caso la energa se puede escribir como:
E = E
100
+E
nlm
+E,
donde E se calcula como en el caso anterior
E =
_
1, 0, 0; n, l, m
e
2
r
12
1, 0, 0; n, l, m
_
= I J,
I =
_
d
3
x
1
_
d
3
x
2
[
100
(x
1
)[
2
[
nlm
(x
2
)[
2
e
2
r
12
d
3
x
2
, .
J =
_
d
3
x
1
_
d
3
x
2
100
(x
1
)
100
(x
2
)
nlm
(x
2
)
nlm
(x
1
)
e
2
r
12
d
3
x
2
,
I, J > 0.
Si los electrones estan en singlete estos pueden estar mas cerca y por lo
tanto su energa es mayor, debido al termino de repulsion electroestatica.
PARAHELIO. Si los electrones se encuentran en el triplete deben estar mas
alejados y su energa es menor. ORTHOHELIO.
Problema 22. Determinar explcitamente una expresion reducida para los
terminos I y J de acuerdo al procedimiento anterior.
Solucion:
Las expresiones para I y J estan dados por:
I =
_
d
3
x
1
_
d
3
x
2
[
100
(
x
1
)[
2
[
nlm
(
x
2
)[
2
e
2
r
12
, (3.1)
J =
_
d
3
x
1
_
d
3
x
2
100
(
x
1
)
nlm
(
x
1
)
nlm
(

x
2
)
100
(
x
2
)
e
2
r
12
. (3.2)
El termino de interaccion entre los electrones se puede expandir como:
1
r
12
=
l=0
r
l
<
r
l+1
>
4
2l + 1
l
m=l
Y
m
l
(
1
,
1
)Y
m
l
(
2
,
2
). (3.3)
Las funciones de onda para un atomo Hidrogenoide con Z = 2 estan dadas
por:
nlm
(

x
j
) = R
nl
(r
j
)Y
m
l
(
j
,
j
), (3.4)
donde n representa el n umero cuantico principal y l y m son, respectivamente,
el momemtum angular orbital y su compoentne en el eje z. El ndice j denota
las coordenadas de la partcula j, con j = 1, 2. En forma expltica el estado
fundamental se representa como:
100
(

x
j
) = R
10
(r
j
)Y
0
0
(
j
,
j
). (3.5)
Debemos integrar en todo el espacio y los elementos de integracion son d
3
x
j
=
r
2
j
dr
j
d
j
.
Para el termino J al reemplazar las expresiones anteriores en (3.2)
J = 4e
2
_
r
2
1
dr
1
_
r
2
2
dr
2
R
10
(r
1
)R
10
(r
2
)R
n
1
l
1
(r
2
)R
n
1
l
1
(r
1
)
l=0
r
l
<
r
l+1
>
1
2l + 1
l
m=l
_
d
1
Y
m
1
l
1
(
1
,
1
)Y
m
l
(
1
,
1
)
_
d
2
Y
m
l
(
2
,
2
)Y
m
1
l
1
(
2
,
2
).
Recuerde que:
_
d
1
Y
m
1
l
1
(
1
,
1
)Y
m
l
(
1
,
1
) =
m
1
,m
l
1
,l
(3.6)
As, nuestra integral se nos reduce a
J =
4e
2
2l + 1
_
r
2
1
dr
1
_
r
2
2
dr
2
R
10
(r
2
)R
n
1
l
(r
2
)
r
l
<
r
l+1
>
=
4e
2
2l + 1
_

0
dr
1
R
10
(r
1
)R
n
1
l
(r
1
)
1
r
l1
1
_
r
1
0
dr
2
R
10
(r
2
)R
n
1
l
(r
2
)r
l+2
2
+
4e
2
2l + 1
_

0
dr
1
R
10
(r
1
)R
n
1
l
(r
1
)r
l+2
1
_
r
r
1
dr
2
R
10
(r
2
)R
n
1
l
(r
2
)
1
r
l1
2
58
En forma analoga podemos calcular el valor de I, con lo cual se obtiene
que:
I =
_
d
3
x
1
_
d
3
x
2
[
100
(
x
1
)[
2
[
n
1
l
1
m
1
(
x
2
)[
2
e
2
r
12
=
_
d
3
x
1
_
d
3
x
2
100
(
x
1
)
100
(
x
1
)
n
1
l
1
m
1
(
x
2
)
n
1
l
1
m
1
(
x
2
)
e
2
r
12
= e
2
_
r
2
1
dr
1
[R
10
(r
1
)]
2
_
r
2
2
dr
2
[R
n
1
l
1
(r
2
)]
2
l=0
r
l
<
r
l+1
>
4
2l + 1
l
m=l
_
d
1
Y
0
0
(
1
,
1
)Y
m
l
(
1
,
1
)
_
d
2
Y
m
1
l
1
(
2
,
2
)Y
m
1
l
1
(
2
,
2
)Y
m
l
(
2
,
2
)
= 4e
2
__

0
r
1
dr
1
[R
10
(r
1
)]
2
_
r
1
0
r
2
2
dr
2
[R
n
1
l
1
(r
2
)]
2
+
_

0
r
2
1
dr
1
[R
10
(r
1
)]
2
_

r
1
r
2
dr
2
[R
n
1
l
1
(r
2
)]
2
_
= 4e
2
_
1
a
3
_

0
r
1
dr
1
e
4r
1
a
_
r
1
0
r
2
2
dr
2
[R
n
1
l
1
(r
2
)]
2
+
_

r
1
r
2
dr
2
[R
n
1
l
1
(r
2
)]
2
_
Para evaluar en forma explcita las integrales se pueden emplear las
siguientes expresiones:
R
n,l
( =
2Z
na
r) =
(n l 1)!
2n[(n +l)!]
3
_
2Z
na
_
3
exp
_
2
_
l
L
2l+1
n+1
(), (3.7)
donde L
2l+1
n+1
() son los polinomios asociados de Laguerre.
3.5.2. Metodo variacional
En este caso solo calcularemos la correcion al estado 1s
2 1
S, para esto
se considerara como funcion de prueba:
(r) =
1
a
_
exp
_
a
r
_
.
Observe que el parametro corresponde a considerar una carga efectiva
diferente de Z = 2. Recuerde que para un atomo hidrogenoide el estado
fundamental esta dado por:
1s
(r) =
1
_
Z
a
_
exp
_
Z
a
r
_
, a =

2
me
2
(radiode Bohr) .
E
1s
=
Z
2
e
2
2a
Al emplear esta funcion de prueba se tiene que:
(1s)
2 =
_
e
2
r
12
_
=

6
Ze
2
2
a
6
o
_ _
exp(
a
o
(r
1
+r
2
))
1
r
12
d
3
r
1
d
3
r
2
=
16
6
Ze
2
a
6
o
_ _
exp(
a
o
(r
1
+r
2
))
1
r
12
r
2
1
dr
1
r
2
2
dr
2
.
La energa total obtenida en este caso es:
E
=
_
2
_
2
2Z
+
5
4
_
e
2
2a
E
mn
=
_
Z
2
5
4
Z +
25
128
_
e
2
2a
,
mn
= Z
5
16
.
Problema 23. Demuestre que la energa obtenida para el estado fundamental
del atomo de Helio empleando el metodo variacional es
E
mn
=
_
2Z
2
5
4
Z +
25
128
_
e
2
2a
, (3.8)
donde
mn
= Z
5
16
es el valor que minimiza la energa.
Solucion: En el calculo de la correccion al estado 1s
2 1
S emplearemos el
metodo variacional con una funcion de onda normalizada dada po:
(r) =
1
a
_
3/2
exp
_
a
r
_
(3.9)
donde es el parametro variacional. Recordemos que para un atomo de
Hidrogeno cuyo Hamiltoniano esta dado por:
H =
p
2
2m

Ze
2
r
2
1
, (3.10)
el estado fundamental es
1s
(r) =
1
_
Z
a
_
3/2
exp
_
Z
a
r
_
, (3.11)
donde a =

2
me
2
es el radio de Bohr del atomo y la energa de este estado
es E
1s
=
Z
2
e
2
2a
. Por lo tanto empleando el metodo variacional la energa se
determina como el mmino del valor medio del Hamiltoniano en el estado [)
E() = [H[) , (3.12)
ya que se encuentra normalizado y H esta dado por:
H = H
1
+H
2
+V
int
=
_
p
2
1
2m

Ze
2
r
2
1
_
+
_
p
2
2
2m

Ze
2
r
2
2
_
+
e
2
r
12
.
60
Aqu, los terminos
Ze
2
r
2
i
representean la accion coulumbiana entre cada partcu-
la y el n ucleo, y el termino
e
2
r
12
representa el potencial de interaccion mutua
entre los electrones. En el caso en consideracion el estado se escribe como:
(r
1
, r
2
) =
(r
2
)
(r
1
) =

3
a
3
exp
_
(r
1
+r
2
)
a
_
(3.13)
Calculando separadamente los promedios tenemos
Solucion: Recuerde que en coordenadas esfericas
2
=
1
r
2
r
r
2
r
y
_
0
x
n
e
x
= n!. En este caso la funcion de prueba no es normalizada, as es
que determinamos esta constante: x = 2r/r
o
[H[) =
_
r
2
1
dr
1
d
1
_
r
2
2
dr
2
d
2
r
1
r
2
[H[ r
1
r
2
)
T
=
_
r
2
1
dr
1
d
1
_
r
2
2
dr
2
d
2
_
r
1
r
2
_
p
2
1
2m
+
p
2
2
2m

Ze
2
r
2
1
Ze
2
r
2
2
+
e
2
r
12
_
r
1
r
2
_
Notando que en coordenadas esfericas

2
2m
2
1
=

2
2m
_
1
r
2
1
r
1
r
2
1
1
r
1
_
2
=

2
2m
a
_
1
r
2
1
_
r
2
1
1
_
r
1
_

2
2m
2
1
=

2
2m
a
_
a

2
r
1
_

2
2m
_
2
1
+
2
2
_
=
2
m
a
_
a

1
r
1
1
r
2
_

2
2m
_
r
2
1
dr
1
d
1
r
2
2
dr
2
d
2
2
1
+
2
2
_
=
2
m
_
a
_
2
+

2
m
a
_
r
2
1
dr
1
d
1
r
2
2
dr
2
d
2
_
1
r
1
+
1
r
2
_
2
=
2
m
_
a
_
2
+
2
2
m
_
a
_
2
,
donde empleamos que:
_
r
2
j
dr
j
d
j
2
j
r
j
= 4
1
a
_
3
_
a
2
_
2
_

0
x
j
e
x
j
dx
j
=

a
, j = 1, 2
_
r
2
1
dr
1
d
1
r
2
2
dr
2
d
2
_
1
r
1
+
1
r
2
_
2
=
2
a
.
Por lo tanto tambien se obtiene que:
Ze
2
_
r
2
1
dr
1
d
1
r
2
2
dr
2
d
2
_
1
r
1
+
1
r
2
_
2
=
e
2
2a
4Z (3.14)
Luego, la contribucion de los terminos libres a la energa es:
2
m
_
a
_
2
+ 2
2
m
_
a
_
2
e
2
2a
4Z =
e
2
2a
_
2
2
4Z
, (3.15)
donde reemplazamos

2
m
= e
2
a. Para el termino de interaccion entre los elec-
trones se tiene que:
[V
int
[) = e
2
_
r
2
1
dr
1
_
r
2
2
dr
2
l=0
r
l
<
r
l+1
>
4
2l + 1
l
m=l
_
Y
m
l
(
1
,
1
)d
1
_
Y
m
l
(
2
,
2
)d
2
= (4)
2
e
2
_

0
r
2
1
dr
1
_

0
r
2
2
dr
2
2
1
r
1
>
= (4)
2
e
2
__

0
r
1
2
1
dr
1
_
r
1
0
r
2
2
dr
2
2
2
+
_

0
r
2
1
2
1
dr
1
_

r
1
r
2
dr
2
2
2
_
=
e
2
2a
__

0
xe
x
_
x
2
e
x
2xe
x
2e
x
+ 2
_
dx +
_

0
x
2
e
x
e
x
(x + 1) dx
_
=
e
2
2a
5
4
En el calculo anterior se empleo:
_
Y
m
l
(
1
,
1
)d
1
=
4
l,0
m,0
1
a
_
3
_
r
1
0
r
2
e
2r
a
dr =
1
8
_
x
1
0
x
2
e
x
dx =
1
8
_
x
1
2
e
x
1
2x
1
e
x
1
2e
x
1
+ 2
_
1
a
_
3
_

r
1
re
2r
a
dr =
1
4
a
_

x
1
xe
x
dx =
1
4
a
_
e
x
1
(x
1
+ 1)
_
_

0
xe
x
_
x
2
e
x
2xe
x
2e
x
+ 2
_
dx +
_

0
x
2
e
2x
(x + 1) dx =
5
4
El resultado nal para la energa es:
E() =
e
2
2a
_
2
2
4Z +
5
4
_
. (3.16)
Al minimizar con respecto a se tiene que:
E()
=
e
2
2a
_
4 4Z +
5
4
_
= 0

=
_
Z
5
16
_
,
62
reemplazando en la energa
E(
) =
e
2
2a
_
2
_
Z
5
16
_
2
4Z
_
Z
5
16
_
+
5
4
_
Z
5
16
_
_
=
e
2
2a
_
2Z
2
5
4
Z +
25
128
_
3.5.3. Metodo de Hyllerras (1929)
Como mencionamos anteriormente, los otros metodos estan basados
en el hecho que la funciones de onda aproximadas se encuentran factorizadas
y por lo tanto no da cuenta adecuada de la correlacion entre los electrones.
El metodo de Hylleraas considera la denicion de las siguientes coordenadas
relativas:
s = r
1
+r
2
, t = r
1
r
2
, u = [r
12
[ ,
0 t s, 0 u s, 0 s .
Los angulos , , son los angulos de Euler del triangulo formado por el
n ucleo y los electrones. Para estados con momentum angular orbital nulo, la
funcion de onda no depende de los angulos y solo depende de los variables
s, t y u. El elemento de volumen es:
d = u(s
2
t
2
)dudtds.
La ecuacion de Schrodinger
_
2
2m
i
(E V )
_
= 0
se puede escribir en forma variacional como:
_
_
2
2m
(
i
)
2
(E V )
2
_
d = 0
,
donde el potencial es:
V = Ze
2
_
1
r
1
+
1
r
2
_
+
e
2
r
12
.
Este problema se puede reescribir como:
[M L NE] = 0,
M =
_

0
ds
_
s
0
du
_
u
0
dt
_
u(s
2
t
2
)
_
_
s
_
2
+
_
t
_
2
+
_
u
_
2
_
+2s(u
2
t
2
)
u
+ 2t(s
2
u
2
)
t
_
L =
_

0
ds
_
s
0
du
_
u
0
dt
_
4Zsu s
2
+t
2
_
2
N =
_

0
ds
_
s
0
du
_
u
0
dtu
_
s
2
t
2
_
2
.
Al efectuar un reescalamiento de todas las variables espaciales en el mismo
factor k la ecuacion
M L NE k
2
M kL NE.
Ahora es posible efectuar una variacion del parametro k, con lo que se obtiene:
_
k
2
M kL NE
_
= 0 E =
1
4NM
L
2
.
Como funcion de prueba al problema variacional total
(s, t, u) = e
s/2
C
n,2l,m
s
n
u
m
t
2l
.
El valor de la energa del estado 1s
2
es E = 5,80488
e
2
2a
o
y el valor expri-
mental es E = 5,810
e
2
2a
o
para la funcion:
(s, t, u) = e
s/2
_
1 + 0,08u + 0,01t
2
_
.
Problema 24. Sea H
0
el hamiltoniano de una particula de spin 1/2. Asumir
que H
0
corresponde solo a variables orbitales y posee un espectro de dimen-
sion N + 1, con energas distribuidas en forma equidistante entre 0 y N.
Estas energas no son degeneradas en el espacio orbital.
(a) Asuma que tiene un sistema de 4 electrones, cuyo hamiltoniano
total se puede escribir como
H = H
(1)
0
+H
(2)
0
+H
(3)
0
+H
(4)
0
(3.17)
donde el supra ndice denota al correspondiente electron. Encontrar el espec-
tro de energas de H y sus correspondiente degeneraciones, cuando N 4.
Soluci on:
Como se tiene un sistema de cuatro electrones, partculas de spin 1/2,
la funcion de onda total del sistema debe ser completamente antisimetrica.
64
Como la variables de spin del sistema no aparecen en forma expltica en el
Hamiltoniano para cada una de las partculas podemos escribir:
(i)
a
=
n
(r
i
)
(i)
m
s
, con i = 1, 2, 3, 4. (3.18)
donde a n, m
s
es el n umero cuantico total de la partcula. La funcion
de onda total ahora esta dada por:
total
=
P
(1)
P
(1)
a

(2)
b

(3)
c

(4)
d
. (3.19)
Aqu el ndice P denota el n umero de pertumaciones entre electrones, por
lo tanto (1)
P
posee el valor 1 (1) para un n umero par (impar) de per-
mutaciones. En forma alternativa la expresion anterior se puede representar
por:
total
=
1
4!
(1)
a

(1)
b

(1)
c

(1)
d
(2)
a

(2)
b

(2)
c

(2)
d
(3)
a

(3)
b

(3)
c

(3)
d
(3)
a

(3)
b

(3)
c

(3)
d
, (3.20)
donde [[ representa el determinante. Empleando propiedades de los deter-
minantes es directo probar que esta funcion de onda es antisimetrica ante el
intercambio de un par de electrones y que se anula si se tiene dos n umero
cuanticos totales iguales, por ejemplo a = b, pues el determinante tendra dos
columnas iguales.
Ahora, como ejercicio se puede ir conbiando los n umeros cuanticos
desde los mas bajos a los los mas altos para construir el espectro de este
sistema y determinar sus correspondientes degeneraciones.
Problema 25. Discuta acerca de que sucedera con los niveles de energa de
un atomo de Helio si los electrones fuesen bosones de spin 0. Sea lo mas
cuantitativo posible.
4. Cuantizacion del campo electromagnetico
4.1. Introduccion
4.2. Ecuacion de ondas
El campo electromagnetico (CEM) clasico se describe mediante las
ecuaciones de Maxwell
E =
B
t
, (4.1)
1
o
B =
o
E
t
+J (4.2)
E = , (4.3)
B = 0, (4.4)
donde y J son, respectivamente, las densidades de carga y corriente. Se ha
asumido que el medio es nomagnetico. Estas ecuaciones se pueden reescribir
en funcion de los potenciales escalar y vectorial, denidos como:
B = A,
E =
A
t
= 0.
Reemplazando estas deniciones en la ecuacion 2 y 3 se obtiene que:
( A)
2
A+
1
c
2
t
+
1
c
2
2
A
t
2
=
o
J, (4.5)
_
A
t
_
= . (4.6)
Recuerde que c = (
o
o
)
1/2
y A ( A)
2
A. Esto permite
obtener los potenciales en terminos de la distribucion de cargas y corriente.
66
Como sabemos existe un libertad adicional en la determinacion de los
potenciales, es decir, las ecuaciones anteriores permanecen invariantes ante
un cambio del tipo:
A = A
o
, =
o
+

t
,
donde = (r, t) es una funcion arbitraria de la posicion y el tiempo.
Esto se conoce invarianza de gauge. En particular, aqu elegiremos el gau-
ge de Coulomb (es el adecuado para el problema de propagacion de ondas
electromagneticas):
A = 0
2
= 0.
De esta forma las ecuaciones 5 y 6 se escriben como:
2
A+
1
c
2
t
+
1
c
2
2
A
t
2
=
o
J, (4.7)
2
= , (4.8)
y el potencial escalar es solucion de la ecuacion de Poisson:
(r) =
1
4
o
_
(r
)
[r r
[
d
3
r
.
Para simplicar la expresion para el potencial vectorial se descompone la den-
sidad de corriente en sus componentes transversal y longitudinal, de acuerdo
a:
J = J
T
+J
L
, J
T
= 0, J
L
= 0.
En este caso la ecuacion de continuidad se expresa como:
J
L
=
t
, J
L
= ,
con lo que se obtiene:
2
=
t
=
o
t
J
L
=
o
t
.
Ahora se puede determinar la dependencia en para el campo vectorial, el
resultado es:
2
A+
1
c
2
2
A
t
2
=
o
J
T
, (4.9)
Observe que la ecuacion (??) posee la forma de una ecuacion de onda. Su
solucion es:
A(r, t) =

o
4
_
J
T
(r
, t
)
[r r
[
d
3
r
, t
= t
[r r
[
c
,
donde t
es el tiempo retardado.
4. Cuantizaci on del CEM 67
Problema 26. Demostrar que la expresion anterior es la solucion de la ecua-
cion de ondas previa.
Al igual que la densidad de corriente el campo electrico puede ser
separado en sus componentes transversales y longitudinales como:
E = E
T
+E
L
, E
T
= 0, E
L
= 0,
E
T
=
A
t
, E
L
= .
Como el campo magnetico es naturalmente transversal, B = 0, obviare-
mos el subndice. Finalmente, las ecuaciones de Maxwell se pueden resumir
en:
E
T
=
B
t
,
1
o
B =
o
E
T
t
+J
T
E
T
= 0, B = 0.
Las ecuaciones transversales describen las ondas electromagneticas, que estan
solo determinadas por la componente transversal de la densidad de corriente.
4.3. Campo electromagnetico libre
Consideraremos solo una region del espacio en que no existen compo-
nentes transversales de corriente, es decir, J
T
= 0, por lo tanto
2
A =
1
c
2
2
A
t
2
.
Para resolver esta ecuacion consideramos una region c ubica del espacio de
longitud L, de longitud arbitraria. Nos referiremos a ella como cavidad, sin
embargo, no se requiere que las paredes de la cavidad sean reales. Como se
tiene una ecuacion de onda consideramos soluciones del tipo ondas viajeras
con condiciones de borde adecuado.
El potencial vectorial en la cavidad se puede expandir en serie de
Fourier como:
A =
k
_
A
k
(t)e
ikr
+A
k
(t)e
ikr
_
donde los componentes del vector de onda k satisfacen condiciones de borde
periodicas, es decir,
k
j
=
2
j
L
,
j
= 0, 1, 2, 3, . . .
68
para satisfacer el gauge de Coulomb se debe cumplir que A
k
(t) k = A
k
(t)
k = 0. Esto nos indica que existen dos direcciones perpendiculares de A
k
para cada k . Cada una de las componentes de Fourier satisface la ecuacion
de onda, por lo tanto se cumple que:
k
2
A
k
(t) +
1
c
2
2
A
k
(t)
t
2
= 0
2
k
A
k
(t) +

2
A
k
(t)
t
2
= 0,
con
k
= ck. Con el reemplazo A
k
(t) = A
k,
e
i
k
t
se obtiene
A =
_
A
k,
e
i
k
t
+A
k,
e
i
k
t
_
,
donde denota el grado de libertad de polarizacion u
k,
(r) =
()
e
ikr
, por
lo tanto A
k,
= A
k,
u
k,
(r). As se forma un sistema de ejes coordenados,
ortonormal, con los vectores (
(1)
,
(2)
, k/ [k[) Por lo tanto, cada componente
(modo) de A satisface la ecuacion de un oscilador armonico. Observe que se
satisfacen las siguientes propiedades:
1
V
_
d
3
r u
k,
(r) u
,
(r) =
k,k

,
.
1
V
_
d
3
r u
k,
(r) u
k
,
(r) =
k,k

,
Determinaremos el Hamiltoniano del campo usando las densidades de

energa de los campos:
H =

o
2
_
_
[E[
2
+c
2
[B[
2
_
d
3
r
E
k,
=
A
k
(t)
t
= i
k
_
A
k,
e
i
k
t+ikr
A
k,
e
i
k
tikr
_
E
k,
E
,
=
k
_
A
k,
u
k,
(r)e
i
k
t
A
k,
u
k,
e
i
k
t
_
_
A
,
(r)e
i
k
t
A
k
,
u
k
,
(r)e
i
k
t
_
=
k
_
_
_
_
_
A
k,
A
u
k,
(r) u
,
(r)e
i(
k
k
)t
+A
k
,
A
k,
u
k
,
(r) u
k,
e
i(
k
k
)t
A
k,
A
k
,
u
k
,
(r) u
k,
(r)e
i(
k
+
k
)t
A
k,
A
k,
(r) u
,
(r)e
i(
k
+
k
)t
_
_
_
_
_
o
2
_
[E[
2
d
3
r =

o
V
2
k,
2
k
2A
k,
A
k,
De la misma forma se tiene que:
B
k
= A
k
(t) =ik
_
A
k,
u
k,
(r)e
i
k
t
A
k,
u
k,
e
i
k
t
_
.
Por lo tanto, podemos escribir que:
H = 2
o
V
k,
2
k
A
k,
A
k,
.
Problema 27. Desarrollar los pasos intermedios para demostrar la expresion
anterior.
Al efectuar el cambio de variable desde
_
A
k,
, A
k,
_
(Q
k
, P
k
)
denida por
A
k,
=
1
(4
o
)
1/2
(
k
Q
k,
+iP
k,
)
()
k
,
A
k,
=
1
(4
o
)
1/2
(
k
Q
k,
iP
k,
)
()
k
,
H =
1
2
k,
_
P
2
k,
+
2
k
Q
2
k,
_
.
Por lo tanto la cuantizacion del campo corresponde a postular que P
k,
y
Q
k,
son operadores que satisfacen la relacion canonica de commutacion:
_
Q
k,
,
P
k
_
=
k,k

,
,
_
Q
k,
,
Q
k
_
= 0,
_
P
k,
,
P
k
_
= 0.
Ademas, podemos denir la siguiente combinacion lineal de estos operadores:
a
k,
=
1
(2)
1/2
_

Q
k,
+i
P
k,
_
,
a
k,
=
1
(2)
1/2
_

Q
k,
i
P
k,
_
,
los cuales satisfacen
_
a
k,
, a
_
=
k,k

,
,
[a
k,
, a
k
,
] =
_
a
k,
, a
_
= 0.
Al igual que en un oscilador armonico individual se puede denir:
N
k,
=a
k,
a
k,
,
que satisface las relaciones:
[a
k,
, N
k
,
] =
k,k

,
a
k,
,
_
a
k,
, N
k
_
=
k,k

,
a
k,
70
4.4. Estados de fotones
Los estados propios del operador N
k,
corresponden a un n umero da-
do de fotones con momentum y polarizacion dado por (k,). Sin embargo,
en general un estado que represente al campo electromagnetico contendra fo-
tones en diversos modos, por lo tanto un estado arbitrario se denota por el
producto tensorial de los estaods asociados a los diferentes modos:
n
k
1
,
1
, . . . , n
k
m
,
m
, . . .
_
=

n
k
1
,
1
_
. . .
n
k
m
,
m
_
. . .
La cantidad n
k
m
,
m
se denomina n umero de ocupacion del modo k
m
y pola-
rizacion
m
. Un ejemplo de estados de ocupacion es el estado de vaco
[0) =

0
k
1
,
1
_
0
k
2
,
2
_
0
k
3
,
3
_
. . .
0
k
m
,
m
_
. . .
Empleando las propiedades que conocemos para el oscilador armonico pode-
mos escribir el estado arbitrario como:
n
k
1
,
1
,n
k
2
,
2
,n
k
3
,
3
,...,
n
k
m
,
m
,
_
=
k
i
,
i
_
a
k
i
,
i
_
n
k
i
,
i
_
n
k
i
,
i
!
[0) ,
a
k
m
,
m
n
k
1
,
1
, . . . ,n
k
m
,
m
, . . .
_
=
_
n
k
m
,
m
+1
_
1/2
n
k
1
,
1
, . . . ,n
k
m
,
m
+1, . . .
_
a
k
m
,
m
n
k
1
,
1
, . . . ,n
k
m
,
m
, . . .
_
= n
1/2
k
m
,
m
n
k
1
,
1
, . . . ,n
k
m
,
m
1, . . .
_
.
4.5. Fotones como excitaciones mecanico cuanticas del
CEM
Al efectuar la cuantizacion de los coecientes de la expansion en serie
de Fourier del CEM se tiene que:
A =
1
V
1/2
k,
_
c
2
2
_
1/2
_
a
k,
(t)e
ikr
()
k
+a
k,
(t)e
ikr
()
k
_
.
Este operador es conocido como operador de campo o campo cuantizado.
De esta forma el operador Hamiltoniano del campo cuantizado se es-
cribe como
H =
1
2
_
_
E

E+
1
o

B

B
_
d
3
r,
reemplazando los correspondientes operadores se obtiene:
H =
k,
_
a
k,
a
k,
+ 1/2
_
.
Observe que el termino libre nos genera una contribucion innita en el lmite
en que L y dado que el nivel de referencia de la energa se puede elegir
en forma arbitraria, se dene:
H[0) = 0,
por lo tanto H =

k,
a
k,
a
k,
.
H
n
k
1
,
1
, . . . , n
k
m
,
m
, . . .
_
=
k
j
,
j
n
k
j
,
j
n
k
1
,
1
_
. . .
n
k
m
,
m
_
. . .
En el caso de radiacion clasica el momentum total esta dado por la integral
del vector de Poynting en el espacio, calculando la correspondiente cantidad
cuantizada se obtiene:
P =
1
c
_
_
E

B
_
d
3
r.
Observe que un termino tpico en esta expresion contiene factores
1
V

()
k
k
)
k
_
e
i(kk
)r
d
3
r = k
k,k

,
P =
k
j
,
j
k
_
a
k,
a
k,
+ 1/2
_
. =
k
j
,
j
ka
k,
a
k,
.
Problema 28. Demostrar las propiedades previas para el operador de momen-
tum.
El efecto de H y

P sobre el estado de un foton es:
Pa
k,
[0) = ka
k,
[0)
Ha
k,
[0) = a
k,
[0) = c [k[ a
k,
[0) ,
y empleando la relacion relativista para la masa y energa:
E
2
= m
2
c
4
+p
2
c
2
,
en este caso:
m
2
c
4
= ()
2
(k)
2
c
2
= 0.
Ademas, como vimos los fotones se caracterizan por su polarizacion
()
k
estos transforman como vectores y les corresponde un momentum angular
uno. Es posible denir las polarizaciones circulares como:
()
k
=

2
1/2
_
(1)
k
i
(2)
k
_
,
72
los cuales ante rotaciones innitesimales en torno a k se obtiene:
()
k
= i
()
k
m = 1.
Es necesario destacar que el valor m = 0 no existe debido a la condicion de
transversalidad. Observe que:
()
k

()
k
=
()
k

()
k
= 1
()
k

()
k
=
()
k

()
k
= 0
k
()
k
= 0.
4.6. Dependencia temporal de los operadores de
creacion y destruccion
A =
1
V
1/2
k,
_
c
2
2
_
1/2
_
a
k,
(t)e
ikr
()
k
+a
k,
(t)e
ikr
()
k
_
.
En el caso del vacio:
a
k,
(t) =
i
[H, a
k,
(t)] = ia
k,
(t)
a
k,
(t) = a
k,
(0)e
it
A =
1
V
1/2
k,
_
c
2
2
_
1/2
_
a
k,
(0)e
it+ikr
()
k
+a
k,
(0)e
itikr
()
k
_
.
4.7. Estados coherentes
En general, los estados relevantes fsicamente del CEM son superpo-
siciones de los estados de n umero de fotones. De las superposiciones conside-
raremos los estados coherentes denidos por:
[) = e
||
2
/2
n=0
n
(n!)
1/2
[n) .
Es directo demostrar que (ejercicio):
[) = 1
[) = exp
_
1
2
[[
2
1
2
[[
2
+
_
a [) = [)
[) = exp
_
a
1
2
[[
2
_
[0) = exp
_
a
a
_
[0) ,
donde = [[e
i
. Algunas de las propiedades de los estados coherentes son:
[ n[) = [[
2
[ n
2
[) = [[
4
+[[
2
n = [[
[
E [) = 2
_

2
o
V
_
1/2
[[ sin(k r t +)
[
E
2
[) =
_

2
o
V
_
1/2
_
4[[
2
sin
2
(k r t +) + 1
E =
_

2
o
V
_
1/2
Problema 29. Demostrar las propiedades previas de los estados coherentes.
4.8. Operador densidad
Hasta ahora solo hemos considerado casos en que el sistema fsico en
consideracion se encuentra en un estado puro. Sin embargo, generalmente lo
que se tiene es que solo existe una determinada probabilidad de encontrar al
sistema en un estado puro dado. Sean estos estados

j
_
y la probabilidad
de medir el sistema en ese estado es P
j
, en estos casos el sistema se describe
mediante el operador densidad denido por:
=
j
P
j
j
_
.
La matriz densidad posee toda la informacion del sistema fsico!
El valor medio de un observable esta dado por:
A) = Tr (A) =
k
k[ A [k)
=
j
P
j
j
_
,
74
donde [k) es una base ortonormal del sistema fsico. Algunas pro-
piedades de la matriz densidad:
1. Tr () = 1
2.
ij
=
ji
3. Tr
_
2
_
1
4.
ij
2

ii
jj
.
En el caso del campo electrico:
E) = 2
_

2
o
V
_
1/2
n
(n + 1)
1/2
Im(n + 1[ [n)) exp(ik r it +i) .
Problema 30. Distribucion de fotones a temperatura T:
=
n
P
n
[n) n[ , P
n
= [1 exp(/kT)] exp(n/kT)
5. Interaccion de atomos con campos
cuantizados
5.1. Introduccion
Consideraremos el Hamiltoniano de un campo y atomo cuantizados
H =
1
2m
j
_
p
j
+e

A(r
j
)
_
2
+
1
2
_
(r)(r)d
3
r
+
1
2
_
_
E

E+
1
o

B

B
_
d
3
r,
donde
(r) =
j
e(r r
j
) +Ze(r).
La interaccion entre el atomo y el campo corresponde al termino:
H
Inter
=
e
m
A(r
j
) p
j
+
e
2
2m
A(r
j
)
2
.
5.2. Aproximacion Dipolar
En el caso de un electron ligado por un potencial V (r) al n ucleo ubi-
cado en la posicion r
o
, la interaccion entre el atomo y el campo se puede
describir empleando la aproximacion dipolar, en donde el atomo esta inmer-
so en una region del espacio en que el campo esta descrito por el potencial
vectorial
V (r) = A(t) exp(ik (r
o
+r))
= A(t) exp(ik r
o
) (1 +ik r +. . .)
A(t) exp(ik r
o
) .
76
La ecuacion de Schrodinger para el electron en esta aproximacion esta dada
por:
_

2
2m
_
ie
A(t, r
o
)
_
2
+V (r)
_
(t, r) = i
(t, r)
t
.
Esto se ha escrito en el guage de Coulomb o gauge electromagnetico (U(t, r) =
0, A = 0). Efectuamos el siguiente cambio de variable:
(t, r) = exp
_
ie
A(t, r
o
) r
_
(t, r) .
Por lo tanto la nueva ecuacion para (t, r) es:
_

2
2m
2
+V (r) er E(t, r
o
)
_
(t, r) = i
(t, r)
t
,
H = H
o
+H
1
,
H
1
= er E(t, r
o
).
5.3.

Atomo y campo cuantizado
Nos restringiremos en forma exclusiva al caso de la interaccion de un
campo E con un atomo hidrogonoide en la aproximacion dipolar
H = H
a
+H
c
er E.
H
c
=
k
_
a
k
a
k
+
1
2
_
,
H
a
=
i
E
i
[i) i[ =
i
E
i
ii
,
ij
= [i) j[
er = e
i,j
[i) i[ r [j) j[ =
i,j
p
ij
ij
.
E =
k
E
k
_
a
k
+a
k
_
, E
k
=
_

k
3
o
V
_
1/2
.
Al reemplazar todas las deniciones es posible expresar el Hamiltoniano
como:
H =
k
a
k
a
k
+
i
E
i
ii
+
i,j
k
g
(ij)
k

ij
_
a
k
+a
k
_
,
donde la constante de acoplamiento se ha denido como g
(ij)
k
=
p
ij
k
E
k
.
5. Interacci on atomos-campos 77
Problema 31.

Atomo de dos niveles. Asumiremos que la constante de aco-
plamiento es real y que p
ab
= p
ba
(g
(ab)
k
= g
(ba)
k
= g
k
). El Hamiltoniano se
reduce a
H =
k
a
k
a
k
+E
a
aa
+E
b
bb
+
k
g
k
(
ab
+
ba
)
_
a
k
+a
k
_
,
=
k
a
k
a
k
+
1
2
z
+
k
g
k
(
+
+
)
_
a
k
+a
k
_
.
E
a
aa
+E
b
bb
=
1
2
k
(
aa
bb
) +
1
2
(E
a
+E
b
)
z
=
aa
bb
,
+
=
ab
,
=
ba
,
[
ij
,
kl
] =
il
jk
kj
il
.
5.4. Aproximacion de onda rotante
Observe el termino de interaccion en la expresion anterior se puede
escribir como
V
k
= g
k
_
+
a
k
+
+
a
k
+
k
+
a
k
_
.
En donde podemos observar la presencia de los siguientes procesos:
+
a
k
: absorcion del atomo y creacion de un foton,
a
k
: emision del atomo y destruccion de un foton,
+
a
k
: absorcion del atomo y destruccion de un foton,
k
: emision del atomo y creacion de un foton.
Los dos ultimos procesos convervan energa, en el primero se crea aproxima-
damente 2 y en el segundo se pierde esta misma cantidad. La aproximacion
de onda rotante se obtiene al despreciar los dos primero terminos:
H =
k
a
k
a
k
+
1
2
z
+
k
g
k
_
k
+
+
a
k
_
.
Esto se puede comprobar facilmente al calcular el termino de interaccion en
el cuadro de interaccion (ejercicio)
5.5.

Atomo de dos niveles con un unico modo del CEM
H = H
o
c
a
a +
1
2
z
. .
+V g
_
+
+
a
_
. .
.
78
Al describir la interaccion entre el atomo y el campo en el cuadro de inter-
accion obtenemos:
V
I
= e
iH
o
t/
V e
iH
o
t/
.
Los terminos que aparecen al efectuar el cambio de cuadro son del tipo:
e
i
c
a
at
ae
i
c
a
at
= ae
i
c
t
,
e
i
a
z
t
+
e
i
a
z
t
=
+
e
i
a
t
,
por lo tanto se obtiene nalmente:
V
I
= g
_
e
it
+
+
ae
it
_
, =
a
c
.
Determinacion del operador de evolucion: En este caso asumiremos adicio-
nalmente que = 0. Como en esta condicion el termino de interaccion es
independiente del tiempo el operador de evolucion se obtiene directamente
de:
U(t) = exp
_
i
V
I
t
_
= exp
_
igt
_
+
+
a
__
,
=
l=0
(igt)
l
l!
_
+
+
a
_
l
.
Los terminos en esta expansion se pueden separar en:
_
+
+
a
_
2m
=
_
aa
_
m
[a) a[ +
_
a
a
_
m
[b) b[ ,
_
+
+
a
_
2m+1
=
_
aa
_
m
a [a) b[ +a
_
aa
_
m
[b) a[ ,
por lo tanto:
U(t) = cos
_
gt
_
a
a + 1
_
1/2
_
[a) a[ + cos
_
gt
_
a
a
_
1/2
_
[b) b[
i
sin
_
gt
_
a
a + 1
_
1/2
_
(a
a + 1)
1/2
a [a) b[ ia
sin
_
gt
_
a
a + 1
_
1/2
_
(a
a + 1)
1/2
[b) a[ .
[(t)) = U(t) [(0)) .
En el caso que ,= 0 se puede emplear el metodo del operador de
evolucion, por razones didacticas consideraremos la derivacion enpleando una
expansion la base de estados de n umero para el campo:
i
[(t))
t
= V
I
[(t)) ,
[(t)) =
n
[c
a,n
(t) [a, n) +c
b,n
(t) [b, n)]
c
a,n
= ig (n + 1)
1/2
e
it
c
b,n+1
, c
b,n+1
= ig (n + 1)
1/2
e
it
c
a,n
.
Este es el mismo sistema de ecuaciones obtenidos en el caso de un sistema de
dos niveles interactuando con un campo clasico, cuya solucion es:
c
a,n
(t) =
_
_
c
a,n
(0)
_
cos
_
n
t
2
_
i

n
sin
_
n
t
2
_
_
2ig(n+1)
1/2
n
c
b,n+1
(0) sin
_
n
t
2
_
_
_
exp(it/2),
c
b,n+1
(t) =
_
_
c
b,n+1
(0)
_
cos
_
n
t
2
_
+i

n
sin
_
n
t
2
_
_
2ig(n+1)
1/2
n
c
a,n
(0) sin
_
n
t
2
_
_
_
exp(it/2).
n
=
2
+ 4g
2
(n + 1).
Problema 32. Determinar el operador de evolucion del modelo de Jaynes-
Cummings en el caso ,= 0.
Si inicialmente el atomo se encuentra en el estado excitado [a): c
a,n
(0) =
c
n
(0), c
b,n
(0) = 0.
c
a,n
(t) = c
n
(0)
_
cos
_
n
t
2
_
i

n
sin
_
n
t
2
__
exp(it/2),
c
b,n+1
(t) =
2ig(n + 1)
1/2
n
c
n
(0) sin
_
n
t
2
_
exp(it/2).
5.6. Colapos y resurguimientos de la poblacion atomica
Un efecto interesante de analizar en el caso de campo cuantizados es
la evolucion de la poblacion atomica. Recordemos que cuando la interaccion
del atomo de dos niveles se producen con un campos clasico, la evolucion de
la probabilidad de encontrar al atomo en el estado superior es:
P
(clasico)
a
(t) = cos
2
_
t
2
_
, (5.1)
donde =
_
2
+
2
o
, =
a

c
,
o
=
pE
o
. En esta descripcion se
considero una interaccion de dipo dipolar en la aproximacion de onda rotante
y se conoce como la frecuencia de Rabi del sistema.
En el caso de un modo del campo cuantizado y un atomo, en que el
campo se encuentra en un estado coherente de amplitud (el estado mas
clasico permitido por el principio de incertidumbre) y el atomo en el estado
superior esta probabilidad es:
P
(cuantico)
a
(t) =
n
p
n
[c
a,n
(t)[
2
= e
||
2
n=0
[[
2n
n!
cos
2
_
g
n + 1t
_
.
80
Si la amplitud del campo es sucientemente intensa, [[ gt, es posible
aproximar la expresion anterior por:
P
(cuantico)
a
(t)
1
2
_
1 + cos (2 [[ gt) e
g
2
t
2
_
. (5.2)
En terminos gracos podemos apreciar lo que sucede con la probabi-
lidad cuando [[ = 40 en funcion del parametro adimensional gt: Podemos
Figura5.1. Poblaci on del estado atomico superior a tiempos cortos en funcion de
gt. El atomo y el campo se asumieron inicialmente, en el estado superior y en
un estado coherente con = 40, respectivamente. Las oscilaciones de poblacion
at omica decaen con una envolvente del tipo exponencial cuadratica.
apreciar que la poblacion alcanza un estado estacionario cuando
gt
colapso
gt
n +n gt
n n = gt
_
[[
2
+[[ gt
_
[[
2
[[ (5.3)
= gt [[
_
1 +
1
[[

1
1
[[
_
= gt. (5.4)
Es decir, el colapso de la poblacion se produce debido a un amortiguamiento
exponencial que no depende de la intensidad del campo. Observe que este es
el tiempo caracterstico en que se produce el desfasamiento de las contribu-
ciones de las diferente frecuencias de Rabi del problema,
n
= g
n + 1. Este
fenomeno se conoce usualmente como Cummings collapse, podemos decir
que es un efecto de interferencia cuantica destructiva.
Un efecto notable de interferencia cuantica es el de resurgimiento,
revival, de las oscilaciones de la poblacion atomica. Esto se produce debido
a la naturaleza discreta de las oscilaciones de Rabi, inicialmente se produce
el colapso por el desfase. Sin embargo, en un tiempo posterior las oscilaciones
pueden contribuir en forma constructiva, que permiten que reaparezcan las
oscilaciones de la poblacion atomica. Este tiempo esta dado por:
gt
resurgimiento
4 (5.5)
A continuacion se muestra una situacion en que [[ = 4, en que el tiem-
po de resurgimiento es del orden de gt 30, la expresion anterior sera valida
cuando [[ gt. El resurgimiento se produce debido a la naturaleza discreta
Figura5.2. Poblaci on del estado atomico superior a tiempos largos en funcion de
gt. El atomo y el campo se asumieron inicialmente, en el estado superior y en un
estado coherente con = 4, respectivamente. Se puede observar que reaparecen las
oscilaciones de la poblacion atomica.
del campo y la evolucion atomica esta determinada por cuantos del campo.
Caractersticas esenciales de los resurgimientos se aprecia en la gura ante-
rior: (i) nunca alcanzan el valor maximo y (ii) son cada vez mas anchos. As,
el modelo de Jaynes-Cummings a un con su simpleza permite precedir (J. H.
Eberly, N. B. Narozny and J.J Sanchez-Mondragon, Phys. Rev. Lett 44, 1323
(1980) ) comportamientos no triviales en la interaccion radiacion-materia y
que han sido observados en forma experimental (en cavidades por G. Rempe
y H. Walther, Phys. Rev. Lett. 58, 353 (1987) ).
Si el campo se encuentra inicialmente en un estado de n umero, que
corresponde a un estado de intensidad denida, se recupera el comportaniento
clasico con
n
. Solo con la presencia de otras contribuciones de estados
de n umero inicial se pueden producir fenomenos de interferencia cuantica
destructiva y constructiva para diferentes tiempos.
5.7. Decaimiento atomico espontaneo en espacio libre
Para estudiar el decaimiento espontaneo de un atomo excitado hare-
mos uso de la descripcion cuantizada del campo electromagnetico, en donde
el atomo en el estado excitado decae a su estado fundamental entregando la
energa a un continuo de modos del campo electromagnetico a diferencia del
caso anterior en que la interaccion se produce con un unico modo. Para esto
haremos uso de la regla de oro de Fermi descrita en teora de terturbaciones
tiempo dependiente.
82
La evolucion del sistema atomo-campo esta dada por:
[(t)) =
n
_
C
a,n
(t)e
i(
a
+n
c
)t
[a, n) +C
b,n+1
(t)e
i(
b
+(n+1)
c
)t
[b, n + 1)
_
+C
b,0
(t)e
i
b
t
[b, 0) ,
[(0)) = [a, 0) .
Empleando teora de perturbaciones tiempo dependiente a primer orden se
obtiene que:
C
(1)
b,1
(t) = g
e
it
1
, =
a
c
.
C
(1)
b,1
(t)
2
= g
2
e
it
1
2
= 4g
2
sin
2 t
2
2
.
Por lo tanto, la probabilidad de decaer del atomo se obtiene sumando todas
las contribuciones de los diferentes modos (sin considerar polarizacion por el
momento)
P
ab
= 4
_

0
d t
2
(g())
2
T()
sin
2 (
a
)t
2
(
a
)
2
t
2
, (5.6)
donde D() es la densidad de estados, que vara como
2
. La densidad de
estados y la constante de acoplamiento son funciones que varan lentamente
en torno a la frecuencia atomica.
P
ab
=
_
2 (g())
2
T()
_
. .
t
. (5.7)
Esta derivacion simple nos permite visualizar algunas de las caractersticas
principales de la emision espontanea. Sin embargo, no da cuenta de que el
decamiento se produce acoplandose al campo electromagnetico tridimensional
y puede ser empleando para tiempo cortos, pues la regla de oro de Fermi no
describe el decaimiento exponencial de la probabilidad.
5.8. Wigner-Weisskopf teora del decaimiento
espontaneo
De acuerdo al resultado anterior, para efectuar una descripcion mas
preciso del decaimiento espontaneo necesitamos conocer la densidad de modos
del campo. Para esto efectuamos una extension al continuo de la distribucion
discreta de modos, en donde la condicion de cuantizacion es:
k
j
=
2n
j
L
, n
j
= 1, 2, . . . . dn
j
=
L
2
k
j
dn =
L
3
8
3
d
3
k.
Las sumas se modican de acuerdo a:
1
V
k,
f(k, )
_
dT() f(), T() =
1
2
c
2
. (5.8)
Si la funcion no depende de la polarizacion ni del angulo solido se obtiene
que:
1
V
k,
f(k, )
1
V
_
dnf(k, ) =
1
V
L
3
8
3
_
d
3
kf(k, )
=
1
8
3
_
d
2
c
3
_
2
0
d sin
_

0
d f(k, ).
En la descripcion de Wigner-Weisskopf (Z. Phys. 63, 54 (1930)) del
decaimiento espontaneo se considera la siguiente interaccion entre el atomo
y el campo:
H =
_
_
k,
a
k,
a
k,
+
a
aa
+
b
bb
_
_
+V
Int
,
V
Int
=
k,
g
k,
a
k,
+
+ adj.
La condicion inicial es [a, 0) por lo tanto el acoplamiento se produce con
los estados [b, 1
k,
), as el estado del sistema es:
[(t)) = C
a,{0}
(t)e
i
a
t
[a, 0)+
k,
C
b,{1
k,
}
(t)e
i(+
b
)t
[b, 1
k,
) . (5.9)
Al reemplazar esta expresion en la ecuacion de Schrodinger se obtiene:
C
a,{0}
(t) = i
k,
g
k,
e
i(
ab
)t
C
b,{1
k,
}
C
b,{1
k,
}
(t) = ig
k,
e
i(
ab
)t
C
a,{0}
.
C
a,{0}
(t) =
k,
[g
k,
[
2
=
V
8
3
c
3
_
d
2
_

0
d sin [g(, )[
2
_
2
0
d
_
t
0
dt
e
i(
ab
)(tt
)
C
a,{0}
(t
).
Observe, que la expresion para la evolucion del coeciente del estado [a, 0)
depende de tiempos anteriores. Para obtener una expresion aproximada se
84
emplea el concepto de aproximacion de grano grueso (coarse graining ap-
proximation), que corresponde a asumir que C
a,{0}
(t
) vara poco en un
intervalo de tiempo en que los otros coecientes contribuyen en forma apre-
ciable, es decir,
_
t
0
dt
e
i(
ab
)(tt
)
C
a,{0}
(t
) C
a,{0}
(t)
_
t
0
dt
e
i(
ab
)(tt
)
(5.10)
y
lm
t
_
t
0
dt
e
i(
ab
)(tt
)
= (
ab
) P
_
i

ab
_
. (5.11)
El lmite de t
1
qu se emplea la siguiente regularizacion:
lm
t
1 e
ixt
x
= lm
0
+
i
_

0
e
i(x+i)t
dt
= lm
0
+
1
x +i
= lm
0
+
_
x
x
2
+
2

i
x
2
+
2
_
=
1
x
i(x).
en esta expresion tiene el sentido de indicar que la evolucion del sistema
auxiliar (modos del campo en este caso) se produce en una escala de tiem-
po mucho menor que el atomo. El primer termino nos da la contribucion
al decaimiento y la parte principal produce un corrimiento de las energas
atomicas (Lamb shift), nos concentraremos solo en el efecto del termino de
decaimiento. Antes de continuar haremos la suma sobre las polarizaciones,
en donde el unico termino que tiene dependencia en ellas es:
[g(, )[
2
=
1
a[ er
()
[b)
_

2
o
V
2
=
p
2
ab
2
o
V
. (5.12)
Como k = k z,
(1)
= x,
(2)
= y y r = r sin cos x+r sin sin y+r cos z
se obtiene:
[g(, )[
2
=
p
2
ab
sin
2
2
o
V
_
cos
2
+ sin
2
_
=
p
2
ab
2
o
V
sin
2
_

0
d sin
3
=
4
3
C
a,{0}
(t) =
p
2
ab
6
o
c
3
_
d
3
(
ab
)C
a,{0}
(t) (5.13)
donde empleamos
_
0
d sin
3
_
2
0
d =
4
3
2. Finalmente, la ecuacion para
amplitud de probabilidad es:
C
a,{0}
(t) =
2
C
a,{0}
(t), =
1
4
o
4p
2
ab
3
ab
3c
3
. (5.14)
1
A
As, observamos que el acoplamiento del atomo al vaco del campo eletro-
magnetico nos permite obtener el decaimiento exponencial de la poblacion
atomica, P
a
(t) = e
t
. El resultado de esta interaccion es un proceso de tipo
irreversible, debido a la naturaleza continuo de las frecuencias. Si se modica
el vaco y solo se permite un n umero fnito de modos y las frecuencias son
conmensurables el proceso sera nuevamente reversible.
Problema 33. Considere un sistema fsico que puede ser descrito en forma
aproximada como un sistema de dos niveles, siendo estos [0) y [1). El sistema
se acopla a un ba no termico con un n umero nito de grados de libertad 2N+1.
El acoplamiento entre el ba no termico se asume valido en las aproximaciones
dipolar y de onda rotante.
Si inicialmente el sistema se encuentra en el estado excitado, [(0)) =
[1), y el ba no termico en el vaco, [0
n
) [0
N
, 0
N+1
, . . . , 0
n
, . . . , 0
N
), el
Hamiltoniano del sistema y ba no termico se puede expresar como:
H = H
(s)
0
+H
(bt)
0
+V
(sbt)
,
H
(s)
0
=
0
[0) 0[ +
1
[1) 1[ ,
H
(bt)
0
=
N
n=N
n
[1
n
) 1
n
[ ,
V
(sbt)
=
N
n=N
g
n
([0) 1[ [1
n
) 0
n
[ +[1) 0[ [0
n
) 1
n
[) ,
donde [1
n
) [0
N
, 0
N+1
, . . . , 1
n
, . . . , 0
N
) y [0
n
) [0
N
, 0
N+1
, . . . , 0
n
, . . . , 0
N
),
es decir, una excitacion en la frecuencia
n
y todos los modos en el vaco.
Ademas, se sabe que
n
= (
1
0
) n, donde es la constante de
decaimiento del sistema. Por ejemplo, este sistema permite modelar el aco-
plamiento de un atomo de dos-niveles y un conjunto nito de osciladores
armonicos.
De acuerdo a este modelo:
1. Escriba las ecuaciones que permiten determinar el estado del sistema
compuesto al tiempo t. Recuerde que el sistema fsico mas el ba no termico
forman un sistema cerrado y al tiempo inicial se encuentran en un estado
puro [(0)) = [(0)) [0
n
). Para esto asuma que el vector de estado del
sistema se puede escribir como:
[(t)) = A(t) [1) [0
n
) +
N
n=N
B
n
(t) [0) [1
n
) . (5.15)
Explique por que puede asumir esto.
86
2. Asuma que la constante de acoplamiento es del tipo g
n
= g (constante)
y obtenga un expresion aproximada para la probabilidad de encontrar al
sistema original al tiempo t.
3. Solucione en forma n umerica el sistema de ecuaciones y compare con el
resultado anterior. Asuma que N = 10, 100 y 1000. Ademas,

1
0
2
= 10
9
Hz y /2 = 20 10
3
Hz. Graque la probabilidad de encontrar al
sistema original en el estado exitado como una funcion del tiempo para
los distintos valores de N.
4. Resuelva lo anterior empleando g
n
= ge
(
n
(
1
0
))
2
/2
. Compare las
expresiones analiticas obtenidas e indique como ayuda esta suposicion a
la solucion numerica.
5. Graque los elementos de matriz de la matriz densidad del sistema al
tiempo t asumiendo que el sistema se encuentra inicialmente en el estado:
[(0)) = (c
o
[0) +c
1
[1)) [0
n
) , [c
0
[
2
+[c
1
[
2
= 1. (5.16)
Observe que el estado al tiempo t se puede escribir como:
[(t)) = (A
0
(t) [1) +A
0
(t) [1)) [0
n
) +
N
n=N
B
n
(t) [0) [1
n
) . (5.17)
5.9. Coherencia atomica e interferencia
En este seccion veremos como coherencia cuantica y correlaciones en
interaccion radiacion-materia es el origen de diversos fenomenos inesperados,
contraintuitivos.
a |
>
b |
>
c |
>
n
n
1
2
Figura5.3. Diagrama esquematico de un sistema atomico de tres niveles interac-
tuando con dos campos clasicos resonantes.
El primero de estos fenomenos es el atrapamiento coherente de po-
blaciones atomicas, en donde un atomo preparado en una superposicion de
sus niveles permite la cancelacion de la emision o absorcion en condiciones
especiales de los parametros del sistema. Para ello consideraremos un siste-
ma atomico de tres niveles interactuando con dos campos electromagneticos
clasicos. Escribiremos el Hamiltoniono en las apoximaciones dipolar y de onda
rotante, el sistema atomico se muestra en la gura,
H = H
o
+V,
H
o
= (
a
[a) a[ +
b
[b) b[ +
c
[c) c[)
V =
2
_
1
e
i
1
i
1
t
[a) b[ +
2
e
i
2
i
2
t
[a) c[
_
+ h.c.,
donde
j
e
i
j
t
y
j
corresponden a las frecuencia de Rabi compleja (asociada
al acoplamiento atomo-campo en la transicion a j) y la frecuencia del
campo, respectivamente, con j = 1, 2. Los campos se asumen resonantes a
las frecuencias atomicas, es decir,
aj
=
E
a
E
j
=
j
. La evolucion del vector
del estado del atomo esta dado por:
[(t)) = C
a
(t)e
i
a
t
[a) +C
b
(t)e
i
b
t
[b) +C
c
(t)e
i
c
t
[c) ,
C
a
=
i
2
_
1
e
i
1
C
b
+
2
e
i
2
C
c
_
C
b
=
i
2
1
e
i
1
C
a
,

C
c
=
i
2
2
e
i
2
C
a
.
Como condicion inicial emplearemos: [(t)) = cos

2
[b) + sin

2
e
i
[c), para
la cual la solucion de las ecuaciones anteriores es:
C
a
=
i sin(
2
)
1
e
i
1
cos

2
+
2
e
i(
2
+)
sin

2
_
C
b
=
1
2
__
2
2
+
2
1
cos
_
t
2
__
cos

2
2
1
2
e
i(
1
2
)
sin
2
_
t
4
_
sin

2
_
C
c
=
1
2
__
2
1
+
2
2
cos
_
t
2
__
e
i
sin

2
2
1
2
e
i(
1
2
)
sin
2
_
t
4
_
cos

2
_
,
donde
2
=
2
1
+
2
2
. El atrapamiento ocurre cuando:
1
=
2
, =

2
,
1
2
= . (5.18)
En estas condiciones las amplitudes de probabilidad son:
C
a
= 0, C
b
=
1
2
, C
c
=
1
2
e
i
. (5.19)
Observe que en este caso no se produce emision ni absorcion a un cuando exis-
ten campos electromagneticos presentes. Si los campos clasicos empleados son
reemplazados por campos cuantizados es posible describir adecuadamente
fenomenos asociados a luz lenta. Ademas, cuando los campos clasicos po-
seen una amplitud que depende del tiempo podemos efectuar transferencias
adiabaticas de las poblaciones atomicas de los niveles atomicos sin pobla-
miento del nivel superior.
88
Otro fenomeno de coherencia cuantica relacionado al anterior es el
transparencia eletromagneticamente inducida. En este caso nuevamente se
emplea un sistema atomico de tres niveles, en el cual los niveles superiores
son estimulados por un campo clasico intenso, esto permite que bajo ciertas
condiciones la otra transicion sea transparente a la presencia de un campo
electromagnetico resonante, como se muestra en la gura. Nuevamente se
considera un sistema de tres niveles cerrado:
a |
>
b |
>
c |
>
n
n
m
Probe
Drive
Figura5.4. Diagrama esquematico de los niveles atomicos y del campo de prueba
en la transicion |a |b y el campo intenso en la transicion |a |c.
H = H
o
+V +H
dis
,
H
o
= (
a
[a) a[ +
b
[b) b[ +
c
[c) c[)
V =
2
_
p
ab
c
e
it
[a) b[ +
e
i
t
[a) c[
_
+ h.c.,
= i [
a
[a) a[ +
b
[b) b[ +
c
[c) c[] .
El termino describe el decaimiento atomico esponstaneo en forma fenome-
nologica, ya que el termino es no hermtico. Las ecuaciones de movimiento pa-
ra los elementos de matriz de la matriz densidad son (i
t
= [H, ]+
1
2
, ):

ab
= (i
ab
+
ab
)
ab
i
2
p
ab
c
e
it
(
aa
bb
) +
i
2
e
i
cb

cb
= (i
cb
+
cb
)
cb
i
2
p
ab
c
e
it
ca
+
i
2
e
i
+i
ab

ac
= (i
ac
+
ac
)
ac
i
2
e
i
t
(
aa
cc
) +
i
2
p
ab
c
e
it
bc
.
En este problema el campo de prueba E puede ser considerado como una
perturbacion. Ademas, como el campo que estimula la transicion superior
es intenso y debe ser mantenido a todo orden. Se asume que toda la poblacion
inicialmente se encuentra en el estado inferior, es decir,
aa
=
cc
= 0,
bb
= 1.
Reemplazando esta condicion incial y manteniendo hasta primer orden en E
se obtiene que:

ab
= (i
ab
+
ab
)
ab
+
i
2
p
ab
c
e
it
+
i
2
e
i
cb

cb
= (i
cb
+
cb
)
cb
+
i
2
e
i
+i
ab
.
Ahora efectuamos el siguiente cambio de variables:
ab
= r
ab
e
it
,
cb
=
r
cb
e
i(+
ca
)t
. Asumiendo =
ab
y
=
ac
con lo cual se obtiene:
r
ab
= (i+
ab
)r
ab
+
i
2
p
ab
c
+
i
2
e
i
r
cb
r
cb
= (i+
cb
)r
cb
+
i
2
e
i
r
ab
.
As el problema consiste en resolver el sistema de ecuaciones acoplados, como
es de dismension dos su solucion exacta se conoce y es:
X = AX +B,
A =
_
i+
i
2
e
i
i
2
e
i
i+
_
, B =
_
i
2
p
ab
E
0
_
X(t) = A
1
B, A
1
=
_
_
_
_
_
i+
(
ab
+i)(
cb
+i)+
e
i
(
ab
+i)(
cb
+i)+
2i
e
i
(
ab
+i)(
cb
+i)+
i+
(
ab
+i)(
cb
+i)+
_
_
_
_
_
Luego,
ab
(t) =
ip
ab
ce
it
(
cb
+i)
2
_
(
ab
+i) (
cb
+i) +

2
4
_. (5.20)
Esto nos permite determinar la suceptibilidad compleja =
+i
, donde
el primer termino esta relacionado a la dispersion y el segundo a la absorcion.
Recuerde que en un medio lineal la polarizacion esta dada por P =
o
E =
2p
ab
ab
e
it
. Finalmente, se obtiene que:
=
N
a
[p
ab
[
2
__
ab
cb
4
_
+
cb
(
ab
+
cb
)
_
_
_
ab
cb
4
_
2
+
2
(
ab
+
cb
)
2
_
=
N
a
[p
ab
[
2
2
(
ab
+
cb
)
cb
_
ab
cb
4
__
_
_
ab
cb
4
_
2
+
2
(
ab
+
cb
)
2
_ ,
donde N
a
es la densidad de atomos. A continuacion se incluye un graco
en donde aparecen
(solida) y
(punteada) en funcion del parametro

90
Figura5.5. Gr aco de las componentes real e imaginaria de la susceptibilidad lineal
en funcion del parametro adimensional /
ab
x = /
ab
, para
= 2
ab
,
cb
= 10
4
ab
. Observe que en resonancia,
= 0, la parte real y imaginaria se cancelan. En esta situacion practicamente
no hay absorcion y por lo tanto el indice de refraccion es 1, es decir, el medio
atomico es transparente cuando la transicion superior es estimulada por un
campo intenso.
6. Teora de scattering
6.1. Introduccion
El fenomeno de stattering o dispersion de partculas es la forma en
que se puede estudiar un importante n umero de situaciones en fsica, en
que un haz de partculas (fuente) se dirige a otro conjunto de partculas
(blanco). Las partculas resultantes de este proceso de interaccion dependen
de la naturaleza de la interaccion misma. Considere por ejemplo:
1. Un haz de partculas producidas por un acelerador que inciden sobre
otro conjunto de partculas, las que se pueden encontrar tambien en mo-
vimiento en el sistema de laboratorio.
2. Un haz de luz laser incidiendo sobre un cristal o un conjunto de atomos.
La informacion que se puede obtener de esta interaccion es usualmente:
1. Direccion de deteccion
2. N umero de partculas detectadas
3. Energas de las partculas detectadas.
Un experimento tpico es mostrado en la gura adjunta:
Las caractersticas especcas del proceso de dispercion depende fuertemen-
te de la partculas y las energas involucradas en el proceso. En el caso de
partculas compuestas puede producirse desintegracion y/o en el caso de al-
tas energas se puede generar radiacion electromagnetica. Por lo tanto, en
general, el n umero de partculas iniciales es distinto al nal. Estos procesos
son conocidos como reacciones de dispersion.
El proceso de scattering puede ser elastico o inelastico. Desde el punto
de vista clasico en el primero de ellos se produce conservacion energa y mo-
mentum y en el segundo solo momemtum. Sin embargo, en el caso cuantico
92
Blanco
Detector
Detector
Hazincidente
q
q
1
2
Figura6.1. Funciones de Bessel.
estas no son las unicas variables que debemos tener en consideracion y na-
turalmente se debe producir conservacion de todos los n umeros cuanticos en
los procesos elasticos.
Consideraremos el caso de una partcula de masa que es dispersada
debido a la presencia de un potencial. Es este caso el Hamiltoniano del sistema
se puede escribir como:
H = H
o
+V, H
o
=
1
2
p
2
.
Recuerde que este caso corresponde a la interaccion efectiva de dos parculas
de masas m
1
y m
2
con 1/ = 1/m
1
+ 1/m
2
. Si el proceso de dispersion es
elastico la energa se conserva y el estado nal debe satisfacer que:
H
o
[) = E [) ,
(H
o
+V ) [) = E [) .
Podemos considerar el nuevo problema como una perturbacion del problema
original y por lo tanto [) V 0[) . Com este es un problema de energa
denida podemos efectuar el siguiente cambio: (r, t) = (r) exp(iEt/),
con lo que se obtiene:
_
2
2

2
+V (r)
_
(r) = E(r),
con el siguiente cambio de variables el problema se puede escribir en forma
mas compacta: E =

2
k
2
2
, V (r) =

2
2
U(r).
_
2
+k
2
U(r)
_
(r) = 0.
Esta ecuacion debe cumplir con los comportamiento asntoticos esperados, es
decir, para
t : partcula libre, representada por un onda plana ( exp(ikz)).
6. Teora de scattering 93
t : debe contener la parte transmitida, mas la parte desviada. Algu-
nas de las propiedades de esta onda desviada son.
En una direccion dada (coordenadas esfericas) debe tener la forma
exp(ikr)/r. La cual es solucion en la zona que no existe potencial.
Asumiento que existe una distancia r
o
en que U(r) 0.
_
2
+k
2
_
exp(ikr)
r
(
r ).
En general, la dispersion no es isotropica, depende de la direccion de
medicion.
Incorporando las condiciones previas es posible postular una solucion del
tipo:
v
k
(r)r e
ikz
+f
k
(, )
e
ikr
r
,
donde f
k
(, ) es la amplitude de desviacion, la cual naturalmente de-
pende del potencial. Por lo tanto, la funcion de onda tiempo dependiente
es:
(r, t) =
_
dkg(k)v
k
(r)e
iE
k
t/
r
_
dke
ikr
e
iE
k
t/
+
_
dkg(k)f
k
(, )
e
ikz
e
iE
k
t/
r
(r, t)r
_
dke
ikz
e
iE
k
t/
+
_
dkg(k)f
k
(, )
e
ikr
e
iE
k
t/
r
6.2. Determinacion de la seccion transversal de
dispersion
Para determinar la seccion transversal de dispersion emplearemos la
corriente de probabilidad. Recuerde que la corriente de probabilidad se dene
como:
J(r, t) =

2i
[
] =
1
Re
_
__
(r, t) = [(r, t)[
2
,
(r, t)
t
+ J(r, t) = 0.
Incialmente se tiene una onda plana exp(ikz) por lo tanto la corriente
esta dada por:
[J(r, t)[ =
k
.
94
La parte desviada se escribe en coordenadas esfericas y por lo tanto:
=

r
r +
1
r
+
1
r sin
.
La funcion desviada para grandes valores de r es f
k
(, )
e
ikz
r
, de donde
J
r
=
k
1
r
2
[f
k
(, )[
2
,
J
1
r
3
Re
_
1
i
f
k
(, )
f
k
(, )
_
,
J
1
r
3
sin
Re
_
1
i
f
k
(, )
f
k
(, )
_
,
J
, J
r J
r
,
con lo que la parte desviada es esencialmente radial a grandes distancias.
El n umero de partculas, por unidad de tiempo y de area que llegan a
un detector es proporcional a la corriente de probabilidad (la que es propor-
cional al ujo de partculas incidentes):
dn = CJ
d
dS = C(J
d
)
r
r
2
d = C
_
k
_
[f
k
(, )[
2
d
(, ) = [f
k
(, )[
2
(seccion transversal diferencial de dispersion)
6.3. Ecuacion integral de dispersion
En esta seccion emplearemos la funcion de Green asociada a la ecua-
cion para la funcion de onda:
_
2
+k
2
U(r)
_
(r) = 0,
_
2
+k
2
_
G(r) = (r)
(r) =
o
(r) +
_
d
3
r
G(r r
)U(r
)(r
).
_
2
+k
2
_
o
(r) = 0
Las soluciones para la funcion de Green son:
G
(r) =
1
4
e
ikr
r
2
_
e
ikr
r
_
=
1
r
2
_
e
ikr
_
+e
ikr
2
_
1
r
_
+ 2
_
1
r
_

_
e
ikr
_
_
e
ikr
_
= ike
ikr
r
r
,
2
_
e
ikr
_
= k
2
e
ikr
2ik
r
e
ikr
_
2
+k
2
_
G
(r) = 4e
ikr
(r).
Por lo tanto:
(r) = e
ikz
+
_
d
3
r
G
+
(r r
)U(r
)(r
).
Si la zona de accion del potencial esta restringida a una longitud caracterstica
L, entonces asumiremos que: r L y r
L. Observe la gura:
r
r r
r
-
u
[r r
[ r u r
G
+
(r r
) =
1
4
e
ik|rr
|
[r r
[
r
1
4
e
ikr
r
e
ikur
(r)r e
ikz
1
4
e
ikr
r
_
d
3
r
e
ikur
U(r
)(r
),
f
k
(, ) =
1
4
_
d
3
r
e
ikur
U(r
)(r
).
Como es las series de Dyson, empleando la solucion formal podemos escribir:
(r) = e
ik
i
r
+
_
d
3
r
G
+
(r r
)U(r
)(r
)
= e
ik
i
r
+
_
d
3
r
G
+
(r r
)U(r
)e
ik
i
r
+
_
d
3
r
_
d
3
r
G
+
(r r
)U(r
)G
+
(r
)U(r
)(r
).
Si el potencial se puede considerar como una perturbacion, entonces podremos
truncar a alg un orden la solucion formal y obtener una solucion aproximada.
96
En el caso que se desea mantener terminos de primer orden en U(r) se obtiene
el resultado conocido como la Aproximacion de Born:
(r) e
ik
i
r
+
_
d
3
r
G
+
(r r
)U(r
)e
ik
i
r
f
k
(, ) =
1
4
_
d
3
r
e
ikur
U(r
)e
ik
i
r
=
1
4
_
d
3
r
e
i(k
d
k
i
)r
U(r
)
=
1
4
_
d
3
r
e
iKr
U(r
), K = k
d
k
i
.
En esta aproximacion la seccion transversal de dispersion es:
(, ) =

2
4
2
_
d
3
r
e
iKr
V (r
2
.
6.4. Potenciales centrales V (r) = V (r). Ondas
parciales.
En los caso de potenciales centrales el momentum angular es una
constante de movimento y como tal conmuta con el Hamiltoniano, por lo
tanto es conveniente describir el estado de un sistema empleando la ba-
se de momentum angular, estados propios de H, L
2
, L
z
con valores pro-
pios
2
k
2
2
,
2
l(l + 1), m, respectivamente. En este caso el sistema, tam-
bien se puede describir empleando los operadores de momentum, es decir,
H, P
x
, P
y
, P
z
. Como el estado inicialmente se esscribe como una onda pla-
na y despues del proceso de dispersion lo natural es describirlo mediante los
estaods de momentum angular, harmonicos esfericos, es conveniente revisar
la forma de pasar de una descripcion a la otra.
momemtum denido momentum angular denido
P [p) = p[p) P [k) = k [k)
H
o
=
p
2
2
[p) H
o
[k) =

2
k
2
2
[k)
[p[ = (2E)
1/2
k[k
) = (k k
)
r [p) =
1
(2)
3/2
e
ipr/
r[k
) =
1
(2)
3/2
e
ipr/
k =
p
_
d
3
k [k) k[ = 1
[k) =
3/2
[p)
Al resolver el problema de ondas esfericas libres (V (r) = 0) se obtiene
la siguiente solucion:
(0)
k,l,m
(r) =
_
2k
2
_
1/2
j
l
(kr)Y
m
l
(, ),
k, l, m[k
, l
, m
) = (k k
)
l,l

m,m
_

0
dk
l=0
m=ll
[k, l, m)k, l, m[ = 1
donde j
l
(x) son las funciones de Bessel
(0)
k,l,m
(r)r
_
2k
2
_
1/2
Y
m
l
(, )
sin(kr l/2)
kr
=
_
2k
2
_
1/2
Y
m
l
(, )
e
ikr
e
i/2
e
ikr
e
i/2
2ikr
Problema 34. Scattering de un electron por un atomo de hidrogeno: El ato-
mo de hidrogeno posee un electron y por lo tanto el problema en general
debemos tratarlo como un problema de partculas identicas, para construir
los estados permitidos del sistema. Sin embargo, por ahora no consideraremos
esta aspecto de tal forma que el Hamiltoniano del sistema total es:
H = H
o
(r) +H
o
(R) +V (r, R)
H
o
(r) =
2
2
2
r
e
2
r
, H
o
(R) =
2
2
2
R
V (r, R) = e
2
_
1
R

1
[Rr[
_
.
98
Donde la variable r se reere al electron del atomo de Hidrogeno y R a al
electron que sufre el proceso scattering. Los problemas no perturbados se
escriben como:
(H
o
(r) E
n
) [n) = 0,
_
H
o
(R)

2
k
2
2
_
[k) = 0,
cuyas funciones de onda son:
n
(r) = funciones de onda atomo de Hidrogeno
1
e
ikR
= ondas planas normalizadas en el volumen .
Por lo tanto la seccion transversal de scattering en la aproximacion de Born
esta dada por:
d =
_

2
2
_
2
[n, k[ V (r, R) [n
, k
)[
2
sindd,
donde n
= n y k
= k corresponde al scattering elastico y n
,= n y k
2
< k
2
corresponden a scattering inelastico. En este ultimo caso el electron del atomo
se transere desde el nivel n al nivel n
.
En el caso de scattering elastico se satisfacen las siguientes relaciones:
E
n
+

2
k
2
2
= E
n
+

2
k
2
2
,
en terminos vectoriales:
Q
2
= k
2
+k
2
2kk
cos
Q = k
k,
es decir Q es el cambio de momentum que experimenta el electron durante
el proceso de scattering.
k[ V (r, R) [k
) =
e
2
__
1
R
e
iQR
R
2
dRd (cos ) d
_
1
[Rr[
e
iQR
R
2
dRd (cos ) d
_
=
2e
2
__
e
iQR
RdRd (cos )
_
1
[Rr[
e
iQR
R
2
dRd (cos )
_
=
2e
2
__
e
iQRcos
RdRd (cos ) e
iQr
_
e
iQR
dR
d (cos )
_
=
4e
2
Q
2
_
1 e
iQr
_
.
n, k[ V (r, R) [n
, k
) =
4e
2
Q
2
n[ 1 e
iQr
[n
)
=
4e
2
Q
2
_
n,n
n[ e
iQr
[n
)
_
.
Donde empleamos que:
_
e
iQRcos
RdRd(cos ) = 2
sinQR
QR
=
2i
Q
lm
0
_
_
e
i(Q+i)
e
i(Qi)R
_
dR
=
2
Q
lm
0
1
Q+i
+
1
Qi
=
2
Q
lm
0
Qi
Q
2
+
2
+
Q+i
Q
2
+
2
=
2
Q
lm
0
Qi + (Q+i)
Q
2
+
2
=
4
Q
lm
0
Q
Q
2
+
2
=
4
Q
2
Para evaluar el ultimo termino es conveniente emplear coordenadas
parabolicas en que:
x =
cos , y =
sin, z =
1
2
( ) ,
cuya transformacion inversa es:
= r z, = r +z, = tan
1
y
x
,
r =
1
2
( +)
d
3
r =
1
4
( +)ddd, Q r = Qz = Q
1
2
( ) .
Las funciones de onda del atomo de Hidrogeno en coordenadas parabolicas
es:
n,k,m
(, , ) =
1
2
e
im
2k
1
!k
2
!
a
3
o
n
4
[(k
1
+[m[)!(k
2
+[m[)!]
3
F(

na
o
)F(

na
o
),
F(u) = u
1/2|m|
e
u/2
L
|m|
k
1
+|m|
(u), k
1
+k
2
+[m[ + 1 = n
n
1
= k
1
+
1
2
[m[ +
1
2
, n
2
= k
2
+
1
2
[m[ +
1
2
, k
1
, k
2
0
100
con k
1
, k
2
0 enteros.
n[ e
iQr
[n
) = n[ e
iQ
1
2
()
[n
)
=
_

n,k,m
(, , )e
iQ
1
2
()
,k
,m
(, , )
1
4
( +)ddd
= 2
_

n,k,m
(, , )e
iQ
1
2
()
,k
,m
(, , )
1
4
( +)dd
I
j
=
_
e
1
2na
o
L
|m|
k
1
+|m|
(
1
na
o
)e
iQ
2
e
1
2n
a
o
L
|m|
k
1
+|m|
(
1
na
o
)
j
d
I
j
=
_
e
1
2na
o
+
1
2n
a
o
iQ
2
(j|m|)+|m|
L
|m|
k
1
+|m|
(
1
na
o
)L
|m|
k
1
+|m|
(
1
na
o
)d.

Cuantica 1 y 2

Cargado por

Copyright:

Formatos disponibles

Cuantica 1 y 2

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Cuantica 1 y 2

Cargado por

Copyright:

Formatos disponibles

MECANICA CUANTICA I-II

) = (r) +dr (r) = (r) +d r (r)

debido a la rotacion del sistema

(r) corresponde a la funcion de onda transformada unitaria en

, S] = 0 y el valor propio s no depende se m, solo puede depender

nos permite obtener:

,= j se tiene que cuando

[ V [j, m) son no nulos. Al igual que en el caso de los

= j +1. Recuerde que ademas deben satisfacerse las

, las funciones (j) y (j) estan

). As los estados nales que debemos considerar son

adopta valores entre (k

. Por lo tanto las probabilidades de transicion

. Utilizando la ecuacion 2.15 se tiene que:

= 3,634. En forma alter-

en la energa se obtiene que:

. En este caso los estados posibles en cada una de las estadsticas

= 0, por lo tanto la funcion de onda

Notando que en coordenadas esfericas

Determinaremos el Hamiltoniano del campo usando las densidades de

(punteada) en funcion del parametro

= k corresponde al scattering elastico y n

También podría gustarte