Práctica Computacional: La Ecuación de Schrödinger
Partícula confinada en un pozo de potencial unidimensional

Instructor: Darío Mitnik


Parte teórica


La figura 1 representa un ``pozo de potencial" unidimensional (1-$d$ potential well)
$\displaystyle V(x)$ $\textstyle =$ $\displaystyle 0 \hspace{2.69cm} (\vert x\vert>a)$  
$\displaystyle V(x)$ $\textstyle =$ $\displaystyle -\vert E_0\vert \hspace{2cm} (\vert x\vert<a ).$  

Resolver la Ecuación de Schrödinger y obtener el espectro de energías posibles del sistema.

Figure 1: Pozo de potencial unidimensional, de ancho $a$.
\begin{figure}\vspace{1cm}
\hspace{5 cm}
\psfig{figure=potentialwell.eps,width=5.5cm,angle=-0}\vspace{2pt}
\end{figure}


Procedimiento de solucion: separar el problema en distintos sectores,
solucionarlos independientemente y hacer coincidir las funciones (y sus 
derivadas) en los bordes. 
Se debe solucionar separadamente el caso de soluciones pares e impares.


En el sector I (adentro del pozo) reemplazar:
$k_i^2 = 2m(E_0 + E)/\hbar^2 = 2m(E_0 - E_b)/\hbar^2 $ y $p_i = \hbar k_i $
En el sector II (afuera del pozo) reemplazar:
$ \alpha^2 = 2mE_b/\hbar^2 $ y $ \rho^2 = 2mE_0/\hbar^2 $
la continuidad de las funciones y sus derivadas llevaran a 
obtener ecuaciones trascendentales de la forma:
$k_i a ~\tan k_i a = \alpha a$ para el caso par
$k_i a ~\cot k_i a = -\alpha a$ para el caso impar
NO INTENTE RESOLVERLAS !

Una forma de resolver este problema es graficando:
$k_i a ~ \tan k_i a $ ó $-k_i a ~ \cot k_i a $
Las soluciones posibles se encuentran en las intersecciones de esta funcion con
$(k_i a)^2 + (\alpha a)^2 = (\rho a)^2$

Preguntas


  1. Utilizar el gráfico dado en la Figura 2 y calcular el primer estado ligado impar, para el caso en que $V_0=1$ a.u. y $a_0=1.75$ a.u..
  2. Trazar un boceto de las primeras soluciones y compararlas con las del pozo infinito.
  3. Analice el problema de una partícula en un pozo de potencial $V = 0$ ($\vert x-x_0\vert>a$) y $V = -\vert E_0\vert$ ($\vert x-x_0\vert<a$).
  4. Cómo cambian las energías si en lugar de $V$ se utiliza $V+\epsilon$?
  5. Cómo cambian las energías si en lugar de tener un potencial de rango $a$, este es $a+L$ ?
  6. Las expresiones de las energías que se obtienen en estos ejercicios son exactas (no hay $\Delta_E$). Qué se puede decir acerca del tiempo que tarda una partícula en decaer desde un estado excitado?

  7. En las siguientes figuras están representados un potencial unidimensional, la función de onda radial correspondiente al 5$^{to}$ nivel de una partícula en ese potencial, y su energía correspondiente. Sólo uno de los gráficos es correcto. Identificar cuál es.

    \begin{figure}\centerline{\psfig{figure=potquestion1.ps,width=9.0 cm,angle=-90}}\end{figure}

  8. En las figuras siguientes están representados potenciales unidimensionales, y sus estados fundamentales, desplazados verticalmente con su energía correspondiente. Cada función de onda es solución de alguno de los potenciales dibujados en la figura. Identificar cuál de estos gráficos es el correcto.

    \begin{figure}\vspace{.75cm}
\centerline{\psfig{figure=potquestion2.eps,width=9.0 cm,angle=-0}}\end{figure}



Figure: Curva $-x \cot(x)$.
\begin{figure}\centerline{\psfig{figure=deuteron.eps,width=14.0 cm,angle=-90}}\vspace{.25cm}
\end{figure}


Parte computacional

Solución de la caja de potencial utilizando Mathematica

1) Autofunciones del potencial

Definicion del ancho y la profundidad y la ubicacion del pozo

RowBox[{RowBox[{Clear[V0, R0, Del, x0] ;, , RowBox[{RowBox[{V0, =,  , 5.}], ;}], ɯ ... #62371;, Poten[r_] = If[r> (R0 - Del/2), -V0, 0] * If[r< (R0 + Del/2), 1, 0] ;}], }]

Plot[Poten[r], {r, 0, 10}] ;

[Graphics:HTMLFiles/well_3.gif]

Soluciones pares

Busqueda de las energias

(* Cell -> Convert To = Input Form *)

Plot[{ξ*Tan[ξ], Sqrt[γ^2 - ξ^2]}, {ξ, 0, γ}, PlotRange -> {0, γ}];

[Graphics:HTMLFiles/well_4.gif]

(* Cell -> Convert To = Traditional Form *)

Plot[{ξ tan(ξ), (γ^2 - ξ^2)^(1/2)}, {ξ, 0, γ}, PlotRange ... 371;BackgroundRGBColor[0, 1, 1], AxesLabel {ξ, "γ"}] ;

[Graphics:HTMLFiles/well_6.gif]

RowBox[{RowBox[{FindRoot, [, RowBox[{{ξ Tan[ξ] η, ξ^2 + η^2γ^2}, ,, RowBox[{{, RowBox[{ξ, ,, 1.4}], }}], ,, {η, .5}}], ]}],  }]

RowBox[{{, RowBox[{RowBox[{ξ, , 1.3549}], ,, RowBox[{η, , 6.17772}]}], }}]

Definicion de los coeficientes

Clear[kr, κ] RowBox[{k, =,  , RowBox[{2,  , *,  , RowBox[{1.3549, /, Del}]}]}] RowBox[{κ,  , =,  , RowBox[{2, *,  , RowBox[{6.17772, /, Del}]}]}]

Clear[A] ; A = ^(-κ (Del/2))/Cos[k (Del/2)] ;

Definicion de las soluciones pares

Clear[solp, C2, Normapar] solp[r_] = If[r<R0 - Del/2, ^(κ (r - R0)), If[r<R0 + Del/2, A Cos[k (r - R0)], ^(-κ (r - R0))]] 

C2 = NIntegrate[ (solp[r])^2  , {r, 0, Infinity}] ; Normapar = 1/Sqrt[C2] ;

Plot[{Poten[r], Normapar * solp[r]}, {r, 0, 10}, PlotRange {-1, 1}] ;

[Graphics:HTMLFiles/well_15.gif]

Soluciones impares

Busqueda de las energias

Plot[{-ξ Cot[ξ], (γ^2 - ξ^2)^(1/2)}, {ξ, 0, γ}, PlotRange {0, γ}] ;

[Graphics:HTMLFiles/well_17.gif]


Created by Mathematica



Darío Mitnik
U.B.A.