Clase 1
Clase 1
Hamilton-Jacobi , Schrödinger
Ecuaciónes de Schrödinger dependiente e independiente del tiempo
Ecuación de continuidad
Partícula libre
Por otra parte, el patrón de interferencia formado por una partícula con…rma
la hipótesis de Louis De Broglie sobre el comportamiento simétrico de la natu-
raleza. Erwin Schrödinger convencido del comportamiento ondulatorio de una
partícula microscópica considera la formulación clásica de Hamilton - Jacobi
que describe el movimiento de una partícula. Esto llevó a Schrödinger a buscar
una ecuación de movimiento para la función de onda (r; t) de una partícula
cuántica que sea compatible con la ecuación de Hamilton –Jacobi, en algún
límite.
1
1 Ecuación de Hamilton - Jacobi
@S
pi = :
@qi
La acción S se considera como un funcional de las coordenadas generalizadas
qi y del tiempo t. En consecuencia, la derivada temporal total de la acción se
puede obtener de dos formas. La primera es la siguiente,
dS @S X @S
= + q_i ;
dt @t i
@qi
@S X
= + pi q_i : (2)
@t i
dS
= L: (3)
dt
Igualando ambas expresiones, (2) y (3), y considerando que el hamiltoniano
H (qi ; pi ) está dado por,
X
H (qi ; pi ) = pi q_i L;
i
obtenemos la relación,
@S X
+ pi q_i = L;
@t i
@S X
+ pi q_i L = 0;
@t i
@S
+ H (qi ; pi ) = 0;
@t
1- Luis Roa Oppliger, Universidad de Concepción, Chile -
2
teniendo en cuenta ahora que los momentos generalizados pi son las derivadas
parciales con respecto a la respectiva coordenada qi de la acción, @S=@qi , obten-
emos la ecuación de Hamilton –Jacobi para la acción S,
@S @S
+ H qi ; = 0: (4)
@t @qi
Por otra parte, el diferencial total dS la podemos obtenter de (2), esto es,
X @S
dS = pi dqi + dt;
i
@t
esta acción coincide con la fase de una onda plana. Debido a este simple
resultado Schödinger propuso que la partícula libre es descrita por una función
de onda (r; t) donde S(r; t) es su fase.
3
donde ~ = h=2 , h es la constante de Planck, r = @ @ @
@x ; @y ; @z y r2 = r r =
@2 @2 @2
@x2 + @y 2 + @z 2 . Reemplazando (6) en cada término de (7) obtenermos,
i
r (r; t)= AeiS(r;t)=~ rS;
~
1 2 i
r2 (r; t) = AeiS(r;t)=~ 2 (rS) + AeiS(r;t)=~ r2 S;
~ ~
@ (r; t) i @S
= A eiS(r;t)=~ ;
@t ~ @t
esto es,
~2 1 2 i i @S
AeiS(r;t)=~ (rS) + AeiS(r;t)=~ r2 S i~A eiS(r;t)=~ = 0;
2m ~2 ~ ~ @t
1 2 ~ 2 @S
AeiS(r;t)=~ (rS) i r S+ = 0;
2m 2m @t
como hemos considerado una partícula libre, entonces podemos reemplazar el
primer término del parentesis cuadrado por el Hamiltoniano H y el término
imaginario se anula, debido a que el momento es constante, esto es,
rS = p;
2
(rS) p2
= = H;
2m 2m
r2 S = r p = 0;
por lo tanto, en el parentesis cuadrado tenemos,
@S
H+ = 0;
@t
que corresponde a la ecuación de Hamilton –Jacobi, sin la necesidad de consid-
erar algún límite.
4
es decir, la amplitud depende de la posición r y del tiempo t. En este caso la
acción S(r; t) satisface la ecuación de Hamilton - Jacobi si la función (r; t)
satisface la siguiente ecuación de onda,
~2 2 @ (r; t)
r (r; t) + V (r) (r; t) i~ = 0:
2m @t
5
la expresión del lado derecho es idénticamente cero 8r y 8t, si todo el término
que está dentro del parentesis cuadrado grande es cero, esto es,
~2 1 2 1 2 i~ 1
0 = r A+ (rS) r A2 (rS) + V (r)
2m A 2m 2m A2
i~ @A @S
+ ;
A @t @t
@S 1 2 ~2 1 2
= + (rS) + V (r) r A
@t 2m 2m A
i~ @A 1
2
2A + r A2 (rS) ;
2A @t m
@S 1 2 ~2 1 2
= + (rS) + V (r) r A
@t 2m 2m A
i~ @A2 1
2
+ r A2 (rS) :
2A @t m
El lado derecho de esta ecuación se anula cuando su parte real y su parte imag-
inaria se anulan idénticamente. Al anular la parte real tenemos,
@S 1 2 ~2 1 2
+ (rS) + V (r) = r A;
@t 2m 2m A
que corresponde a la ecuación de Hamilton – Jacobi en el límite ~2 =m ! 0.
Como vemos este límite da cuenta del principio de correspondencia establecido
por Niels Bohr en 1923; La física clásica debe emerger como una aproximación
de la física cuántica cuando los sistemas aumentan de tamaño. Este enunciado
de Niels Bohr no es un principio, es más bien un área de investigación que busca
determinar los estados cuánticos o condiciones físicas en los cuales los sistemas
cuánticos emulan el comportamiento clásico.
6
1 ~
= Re r ; = complejo conjugado de ,
m i
~
= Re [ i( r )] ;
m
~
= Re f i [Re ( r ) + i Im ( r )]g ;
m
~
= Re [ i Re ( r ) + Im ( r )] ;
m
~
= Im ( r ); (10)
m
~
= ( r r ); (11)
2mi
las expresiones (10) y (11) son las más usadas al momento de calcular explici-
tamente j (r; t) y solucionar problemas.
En consecuencia, de este formalismo, además de la ecuación de Hamilton –
Jacobi con su término no clásico proporcional a ~2 =m, obtenemos la ecuación
de continuidad (9) y la ecuación de Schrödinger dependiente del tiempo para la
función de onda (r; t) de la partícula sometida a un potencial V (r),
~2 2 @ (r; t)
r (r; t) + V (r) (r; t) = i~ . (12)
2m @t
Esta es una ecuación diferencial lineal y determinista para (r; t), esto es, dada
la condición inicial (r; 0) se obtiene la solución (r; t).
Notemos que (12) tiene tres términos, un término con segundas derivadas espa-
ciales, otro término con derivada temporal y un tercer término que da cuenta
del efecto que produce el protencial V (r) al cual está sometida la partícula.
Ahora podemos anunciar una primera versión de los dos primeros principios de
la Mecánica Cuántica:
Primer principio: Toda la información de una partícula o sistema aislado
está contenida en la función de onda (r; t).
Segundo principio: La evolución de (r; t) está dada por la ecuación de
Schrödinger (12).
Max Born en 1926 interpretó a (r; t) como una amplitud de densidad de prob-
abilidad en la posición r y tiempo t. Para ser más preciso, la probabilidad
prob( ; t) de encontrar a la partícula dentro de una región …nita r 2 en un
2
instante t está dada por la suma continua de probabilidades j (r; t)j d3 r 2 ,
R 2
prob( ; t) = r2 j (r; t)j d3 r:
En consecuencia, la suma de las probabilidades sobre todo el espacio cartesiano
trimimensional r = (x; y; z) tiene que ser 1 en 8t, esto es,
R1 R1 R1 2
1 1 1
j (r; t)j dxdydz = 1;
1- Luis Roa Oppliger, Universidad de Concepción, Chile -
7
esto es conocido como la normalización de (r; t) y por esta razón la función
de onda (r; t) debe ser de módulo cuadrado integrable.
Como vimos en el experimento de la doble rendija y como aclararemos aun más
en las próximas clases, el indeterminismo se presenta al observar una variable
de la partícula, por ejemplo, la posición X en el caso del experimento de la doble
rendija, la componente Sz sobre el eje z del momento angular en el experimento
de Stern-Gerlach.
~2 2 @ (t)
r (r) + V (r) (r) (t) = i~ (r) ;
2m @t
1 ~2 2 i~ @ (t)
r (r) + V (r) (r) = ;
(r) 2m (t) @t
ambos lados dependen de variables independientes, por lo tanto se debe satis-
facer que,
i~ @ (t)
= E; (13a)
(t) @t
1 ~2 2
r (r) + V (r) (r) = E; (13b)
(r) 2m
donde E es una constante con unidades de energía y que corresponde a la energía.
De la ecuación (13a) se obtiene la siguiente ecuación y solución,
@ (t) iE iEt=~
= (t) ; ! E (t) = e ;
@t ~
la constante de normalización de E (t) la consideramos en (r). De (13b) se
obtiene,
~2 2
r + V (r) E (r) = E E (r) ; (14)
2m
que corresponde a una ecuación de funciones propias y respectivos valores
propios f E (r) ; E 2 Rg. Esta ecuación (14) es conocida como la ecuación de
Schrödinger independiente del tiempo. Claramente es V (r) quien determina
las soluciones propias f E (r) ; E 2 Rg.
1- Luis Roa Oppliger, Universidad de Concepción, Chile -
8
Para potenciales V (r) que con…nan a la partícula en una región …nita del es-
pacio, el conjunto de soluciones f E (r) ; E 2 Rg, resulta ser discreto, esto es
f n (r) ; En 2 R; con n enterog. Por la linealidad de la ecuación (14) una
solución general (r; t) siempre puede ser escrita como la superposición de
iEn t=~
las soluciones particulares n (r) e , esto es,
X
(r; t) = Cn n (r) e iEn t=~ ; (15)
n
donde las constantes Cn pueden ser complejas. Como veremos más adelante
f n (r) = hrj n ig son las representaciones de los estados propios fj n ig del
operador hamiltoniano H y los valores propios fEn g es el conjunto de valores de
energía permitidos, conocido también como el espectro de energías del sistema.
Algo que demostraremos en la 2da clase de álgebra lineal es que las funciones
propias f n (r)g tiene la propiedad de ser ortogonales entre sí, esto es,
ZZZ
n (r) m (r) dr = n;m delta de Kronecker,
R3
así en t = 0 tenemos,
ZZZ
Cn = n (r) (r; 0) dr; (16)
R3
por lo tanto, dada la función de onda inicial (r; 0) esta relación (16) permite
obtener los coe…cientes Cn y en consecuencia obtenemos (r; t) de la expresión
(15).
Ahora podemos enuncial una 1ra versión del tercer principio de la mecánica
cuántica: Al medir u observar una variable, en este caso la energía de la
1- Luis Roa Oppliger, Universidad de Concepción, Chile -
9
partícula, su función de onda (r; t) se proyecta en forma completamente aleato-
ria a una de las funciones de onda propia n (r) y el resultado que se observa
en el medidor o detector es el valore propio En asociado a la función propia
n (r) favorecida en el proceso aleatorio. Esquemáticamente esto es,
prob(En ;t)
Menerg{a [ (r; t)] ) (r; t) ! n (r) :
En
10
Representación geométrica en el caso que (r; t) involucre solo tres funciones
propias f n 1 (r) ; n (r) ; n+1 (r) g.
11
3
d2psi(x)/dx2
2
0
0 1 2 3 4 5 6
x
dpsi(x)/dx 8
6
0
0 1 2 3 4 5 6
x
psi(x) 20
15
10
0
0 1 2 3 4 5 6
x
1- Luis Roa Oppliger, Universidad de Concepción, Chile -
12
Si para x 2 [x0 ; x1 ] el potencial V (x) = 1, entonces la función de onda
(x) = 0 en [x0 ; x1 ]. De esta forma d2 (x) =dx2 y d (x) =dx se anulan
en [x0 ; x1 ]. Esto garantiza que (x) sea normalizable.
Si V (x) es un potencial con una barrera del tipo delta de Dirac en x0 ,
entonce de la ecuación de Schrödinger independiente del tiempo tenemos
que,
d2 (x) 2m
= 2 [ (x x0 ) E] (x) ; 2 R;
dx2 ~
d (x) d (x) 2m
= (x0 ) ; (18)
dx x=x0+ dx x=x0 ~2
13
Así, d (x) =dx tendrá la discontinuidad (18) en cada punto donde haya
un potencial delta de Dirac. Note que pude ser positivo o negativo,
> 0 corresponde a una barrera con altura in…nita y < 0 corresponde
a un pozo con profundidad in…nita, ambos con ancho in…nitesimal.
5 Partícula libre
Para una partícula libre el potencial V (x) = 0, 8 x 2 [ 1; 1] y la ecuación
de Schrödinger independiente del tiempo se reduce a la forma,
d2 (x) 2mE
= (x) ;
dx2 ~2
2
p
= (x) ; con p2 = 2mE;
~2
cuyas soluciones o funciones propias son,
eixp=~ p
p (x) = p ; con p = 2mE; (19)
2 ~
así de…nimos,
p2
Ep =; con 1 < p < 1; (20)
2m
note que en este caso el espectro de valores de penergías permitidas Ep es con-
tinuo. La variable continua p es el momento y 2 ~ es una constante de nor-
malización a la delta de Dirac. Esto signi…ca lo siguiente,
Z 1 Z 1
1 1 0
p (x) p 0 (x) dx = p e ixp=~ p eixp =~ dx;
1 1 2 ~ 2 ~
Z 1
1 0
= e ix(p p )=~ dx;
2 ~ 1
| {z }
2 ~ (p p0 )
1
= 2 ~ (p p0 ) ;
2 ~
= (p p0 ) ; delta de Dirac
0; si p 6= p0 ;
= (21)
1; si p = p0 :
Considerando la parte temporal, en este caso las soluciones particulares son las
ondas planas con momento p bien determinado,
1
p (x; t) = p eixp=~ e iEp t=~ ;
2 ~
1 2
= p ei(xp p t=2m)=~ ; con 1 < p < 1:
2 ~
1- Luis Roa Oppliger, Universidad de Concepción, Chile -
14
Una onda plana positiva (p > 0) se propagan en el sentido positivo del eje x y
una onda plana negativa (p < 0) se propaga en el sentido negativo del eje x. La
velocidad de fase de una onda plana depende de p y es vp = @S=@x
m = p=m.
2
La siguiente …gura muestra Re[ei(px p t=2m)=~ ] = cos[(px p2 t=2m)=~] para
dos ondas planas positivas, una con p = 1 (negra) y otra con p = 2 (gris) en
diferentes tiempos t=2m~ = 0, 0:5 y 1.
-1
-5 -4 -3 -2 -1 0 1 2 3 4 5
t---> x/hbar
La …gura claramente muestra que la onda plana con p = 2 hse propaga al doble
i
2
de la velocidad de la onda plana con p = 1. La parte Im ei(px p t=2m)=~ =
sin px p2 t=2m =~ tiene el mismo comportamiento temporal. En conse-
cuencia las ondas planas se dispersan o separan entre sí con el paso del tiempo.
Por linealidad de la ecuación de Schrödinger, cualquier solución (x; t) puede
2 p
ser representada como una suma de ondas planas ei(xp p t=2m)=~ = 2 ~, esto
es,
Z 1
(x; t) = (p) p (x) e iEt=~ dp; (22a)
1
Z 1
2
= p21 ~ (p) ei(xp p t=2m)=~ dp: (22b)
1
15
Z 1 Z 1
= (p0 ) p (x) p0 (x) dx dp0 ;
1 1
Z 1
0
= (p ) (p p0 ) dp0 ;
1
= (p) ;
esto es, Z 1
p1 ixp=~
(p) = 2 ~
(x; 0) e dx (23)
1
6 Ejemplo
La función de onda inicial (x; 0) de una partícula está dada por la siguiente
expresión,
1 x2 =4 2
(x; 0) = 1=4 p e x : (24)
(2 ) x
Note que la densidad de probabilidad es una distribución gaussiana normalizada,
2 2
e x =22 x
j (x; 0)j = p ; (25)
2 x
con,
Z 1
2 1 R1 x2 =2 2
j (x; 0)j dx = p 1
e x dx;
1 2 x| {z }
p
2 x
1 p
= p 2 x;
2 x
= 1:
En este caso, inicialmente la posición promedio hXi = 0 ya que la distribución
2
o densidad de probabilidad j (x; 0)j es una función par, por lo tanto,
R1 2
hXi = 1
x j (x; 0)j dx;
R 1 2 2
= p2 1 xe x{z=2 x} dx;
1 |
x
función impar
= 0:
Por otra parte, en t = 0 el promedio del cuadrado de la posición hX 2 i = 2
x,
esto es,
1- Luis Roa Oppliger, Universidad de Concepción, Chile -
16
R1 2
hX 2 i = x2 j (x; 0)j dx;
1
R 1 2 x2 =2 2
= p 1 x e {z x dx;
2 x 1| }
p
2 3
x
p 3
= p 1 2
2 x
x;
2
= x;
2
La siguiente …gura muestra j (x; 0)j como función de x para diferentes valores
de x , 0:3 (sólida), 1 (segmento) y 3 (punto).
1.0
0.5
0.0
-12 -10 -8 -6 -4 -2 0 2 4 6 8 10 12
x
x2 =2 2 p
En el límite x ! 0 la distribución gaussiana e x = 2 x ! (x).
17
Tarea 1
Una partícula libre inicialmente es descrita por la función de onda (x; 0) dada
por la siguiente expresión,
1 x2 =4 2
(x; 0) = 1=4 p e x :
(2 ) x
Recuerde que, p
X (t) = hX 2 (t)i hX (t)i2 ;
con, Z 1
2
hX 2 (t)i = x2 j (x; t)j dx:
1
———————————————————————————————————
18