This commit is contained in:
pkirwin 2021-03-05 15:41:55 -07:00
parent 05a3d276e9
commit a09d982b30
4 changed files with 29 additions and 4 deletions

Binary file not shown.

Binary file not shown.

View File

@ -680,9 +680,28 @@ axis([0 1e-9 0 0.04]);
\item %4a
\begin{figure}[H]
\centering
\begin{subfigure}{0.5\textwidth}
\centering
\includegraphics[width=\textwidth]{q4a_e1.png}
\caption{}
\end{subfigure}%
\begin{subfigure}{0.5\textwidth}
\centering
\includegraphics[width=\textwidth]{q.png}
\caption{}
\end{subfigure}
\caption{(a),(b) Previous plots but with the classical momentum marked in red.}
\end{figure}
% TODO I really don't like this explanation...
An \( a \) value of \( 0.62 \angstrom \) was chosen, in order to provide an adequately
shaped graph without sacrificing too much computation time.
An \( a \) value of \( 0.53 \angstrom \) was chosen in order to provide an adequately
shaped graph without sacrificing too much computation time and ensure that the first two
numerical energies correspond to the given experimental results. The experimental results are
\(\SI{0.14395}{\electronvolt}\) and \(\SI{0.43185}{\electronvolt} \) for the first and second energy levels
respectively, and the numerical results with our chosen \(a\) are \(\SI{0.0.14386}{\electronvolt}\) and \(\SI{0.43140}{\electronvolt} \)
% maybe talk about normalization?
@ -692,7 +711,13 @@ axis([0 1e-9 0 0.04]);
\lipsum[1]
\item %4bii
Between \( 0V \) and \( 0.25V \), only the first energy level is carrying any current.
This current drops to 0 above \( 0.25V\) because the coupling between the contacts and that
energy level drops to 0, meaning no electrons can transfer.
Between \(0.4V\) and \(0.65V\), only the second energy level is carrying current.
This energy level stops dropping current because it's shifted energy drops below
the threshold where the contacts have any coupling with it.
\item %4biii

View File

@ -18,7 +18,7 @@ w0 = (0.2879*q/hbar);
% YOU MUST ENTER AN APPROPRIATE VALUE OF "a"
a = 12 / sqrt(m*w0/hbar) / 100;
a = 12 / sqrt(m*w0/hbar) / 120;
beta = sqrt(m*w0/hbar);
x = [-6/beta:a:6/beta];
N = length(x);