9.6: Particle in a Semi-infinite Potential Well
- Page ID
- 135868
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
Numerical Solutions to Schrödinger's Equation for the Particle in the Semi-infinite Box
Parameters go here: \( x_{min} = 0\) \( x_{max} = 5\) \( m = 1\) \( lb = 2\)
Potential energy \( V(x) = if[(x \geq lb), V_{0}, 0]\)
Given \( \frac{d^{2}}{dx^{2}} \psi (x) = 2 m (V(x) -E) \psi (x)\) \( \psi (x_{min}) = 0\) \( \psi '(0) = 0.1\)
\( \psi := Odesolve (x, x_{max})\) \( \psi = \frac{ \psi (x)}{ \sqrt{ \int_{x_{min}}^{x_{max}}} \psi(x)^{2} dx}\)
Enter energy guess: E = 0.766
Calculate the probability that the particle is in the barrier: \(\int_{2}^{5} \psi (x)^{2}dx = 0.092\)
Calculate the probability that the particle is not in the barrier: \(\int_{0}^{2} \psi (x)^{2}dx = 0.908\)
Calculate and display the momentum distribution:
Fourier transform: p = -10,-9.9 .. 10 \( \Phi (p) = \int_{ -x_{min}}^{ x_{max}} exp(-1 p x) \psi (x)~dx\)