| Both sides previous revisionPrevious revisionNext revision | Previous revision |
| general_considerations [2020/07/16 19:09] – [4.iii.4 Boundary Conditions] admin | general_considerations [2020/07/21 05:39] (current) – [4.iii.4 Boundary Conditions] admin |
|---|
| The notation $\lim_{\epsilon\rightarrow 0_+}$ means that we take the limit from the positive side, i.e. we start with a small positive value of $\epsilon > 0$ and then consider what happens as it gets smaller and smaller. Similarly, the notation $\lim_{\epsilon\rightarrow 0_-}$ means taking the limit from the negative side with $\epsilon < 0$. | The notation $\lim_{\epsilon\rightarrow 0_+}$ means that we take the limit from the positive side, i.e. we start with a small positive value of $\epsilon > 0$ and then consider what happens as it gets smaller and smaller. Similarly, the notation $\lim_{\epsilon\rightarrow 0_-}$ means taking the limit from the negative side with $\epsilon < 0$. |
| |
| The free particle solutions satisfy both of these conditions everywhere, so the boundary conditions do not constrain these solutions further. For the infinite square well potential, we started with the general solution to the TISE $\psi(x)=A\sin kx + B\cos k x$ and then applied the first boundary condition. We noted that $\psi(x)$ must be zero in the regions of infinite potential to avoid our solutions having infinite energy. Then, the first boundary condition implies that the solutions must go to zero at the ends of the well, so we applied $\psi(0)=\psi(a)=0$. This did two things. First, it gave us $B=0$, which removes the degeneracy. We now only have one independent solution $A\sin k$ for each energy. The second thing it did was to constrain the possible values of $k$& to be $\frac{n\pi}{a}$ for $n=1,2,\cdots$, which gave us a discrete energy spectrum. Both the nondegeneracy and the discreteness come from the boundary conditions. The fact that the solutions have to be zero outside the potential also made the solutions normalizable. | The free particle solutions satisfy both of these conditions everywhere, so the boundary conditions do not constrain these solutions further. For the infinite square well potential, we started with the general solution to the TISE $\psi(x)=A\sin kx + B\cos k x$ and then applied the first boundary condition. We noted that $\psi(x)$ must be zero in the regions of infinite potential to avoid our solutions having infinite energy. Then, the first boundary condition implies that the solutions must go to zero at the ends of the well, so we applied $\psi(0)=\psi(a)=0$. This did two things. First, it gave us $B=0$, which removes the degeneracy. We now only have one independent solution $A\sin k$ for each energy. The second thing it did was to constrain the possible values of $k$ to be $\frac{n\pi}{a}$ for $n=1,2,\cdots$, which gave us a discrete energy spectrum. Both the nondegeneracy and the discreteness come from the boundary conditions. The fact that the solutions have to be zero outside the potential also made the solutions normalizable. |
| |
| For the infinite square well potential, the second boundary condition holds inside the well, but we did not apply it at $x=0$ or $x=a$ because the potential jumps to infinity at these points. The solutions | For the infinite square well potential, the second boundary condition holds inside the well, but we did not apply it at $x=0$ or $x=a$ because the potential jumps to infinity at these points. The solutions |
| do not have a continuous derivative at $x=0$ or $x=a$. For example, | do not have a continuous derivative at $x=0$ or $x=a$. For example, |
| |
| \[\lim_{\epsilon \rightarrow 0_+} \left . \frac{\mathrm{d}\psi}{\mathrm{d}x}\right |_{\epsilon} = \lim_{\epsilon \rightarrow 0} \frac{n\pi}{a}\cos \frac{n\pi\epsilon}{a} = \frac{n\pi}{a},\] | \[\lim_{\epsilon \rightarrow 0_+} \left . \frac{\mathrm{d}\psi}{\mathrm{d}x}\right |_{x=\epsilon} = \lim_{\epsilon \rightarrow 0} \frac{n\pi}{a}\cos \frac{n\pi\epsilon}{a} = \frac{n\pi}{a},\] |
| |
| \[\lim_{\epsilon \rightarrow 0_-} \left . \frac{\mathrm{d}\psi}{\mathrm{d}x}\right |_{\epsilon} = \lim_{\epsilon \rightarrow 0} 0 = 0.\] | \[\lim_{\epsilon \rightarrow 0_-} \left . \frac{\mathrm{d}\psi}{\mathrm{d}x}\right |_{x=\epsilon} = \lim_{\epsilon \rightarrow 0} 0 = 0.\] |
| |
| In fact, there are no solutions with finite energy that have continuous derivatives at the two ends of the well. | In fact, there are no solutions with finite energy that have continuous derivatives at the two ends of the well. |
| \[\int_{x'-\epsilon}^{x'+\epsilon}\frac{\mathrm{d}^2\psi}{\mathrm{d}x^2}\,\mathrm{d}x=\frac{2m}{\hbar^2}\int_{x'-\epsilon}^{x'+\epsilon}V(x)\psi(x) \,\mathrm{d}x-\frac{2mE}{\hbar^2}\int_{x'-\epsilon}^{x'+\epsilon} \psi(x)\,\mathrm{d}x\] | \[\int_{x'-\epsilon}^{x'+\epsilon}\frac{\mathrm{d}^2\psi}{\mathrm{d}x^2}\,\mathrm{d}x=\frac{2m}{\hbar^2}\int_{x'-\epsilon}^{x'+\epsilon}V(x)\psi(x) \,\mathrm{d}x-\frac{2mE}{\hbar^2}\int_{x'-\epsilon}^{x'+\epsilon} \psi(x)\,\mathrm{d}x\] |
| |
| Note that $\int_{x'-\epsilon}^{x'+\epsilon}\frac{\mathrm{d}^2\psi}{\mathrm{d}x^2}\,\mathrm{d}x= \left . \frac{\mathrm{d}\psi}{\mathrm{d}x} \right |_{x'+\epsilon} - \left . \frac{\mathrm{d}\psi}{\mathrm{d}x} \right |_{x'-\epsilon}$, so taking the limit of the right hand side as $\epsilon \rightarrow 0_+$ will give us the boundary condition on the derivative. Now, if $\psi(x)$ is continuous then the second integral on the right hand side vanishes, so we are left with | Note that $\int_{x'-\epsilon}^{x'+\epsilon}\frac{\mathrm{d}^2\psi}{\mathrm{d}x^2}\,\mathrm{d}x= \left . \frac{\mathrm{d}\psi}{\mathrm{d}x} \right |_{x=x'+\epsilon} - \left . \frac{\mathrm{d}\psi}{\mathrm{d}x} \right |_{x=x'-\epsilon}$, so taking the limit of the right hand side as $\epsilon \rightarrow 0_+$ will give us the boundary condition on the derivative. Now, if $\psi(x)$ is continuous then the second integral on the right hand side vanishes, so we are left with |
| |
| \[\lim_{\epsilon \rightarrow 0_+} \left [ \left . \frac{\mathrm{d}\psi}{\mathrm{d}x} \right |_{x'+\epsilon} - \left . \frac{\mathrm{d}\psi}{\mathrm{d}x} \right |_{x'-\epsilon}\right ]=\frac{2m}{\hbar^2}\int_{x'-\epsilon}^{x'+\epsilon}V(x)\psi(x)\,\mathrm{d}x\] | \[\lim_{\epsilon \rightarrow 0_+} \left [ \left . \frac{\mathrm{d}\psi}{\mathrm{d}x} \right |_{x=x'+\epsilon} - \left . \frac{\mathrm{d}\psi}{\mathrm{d}x} \right |_{x=x'-\epsilon}\right ]=\frac{2m}{\hbar^2} \lim_{\epsilon \rightarrow 0_+} \left [ \int_{x'-\epsilon}^{x'+\epsilon}V(x)\psi(x)\,\mathrm{d}x \right ]\] |
| |
| We will see how this can be used to determine the boundary condition for a delta function potential in section 4.vi. | We will see how this can be used to determine the boundary condition for a delta function potential in section 4.vi. |
| {{ :gensteppotential.png?600 |}} | {{ :gensteppotential.png?600 |}} |
| |
| Note that the potential has a minimum $V_0$ and values $V_+$ and $V_-$, which extend to $x = +\infty$ and $x = -\infty$ respectively. The features that the energy eigenstates have depend on how the energy $E$& is related to these three values. | Note that the potential has a minimum $V_0$ and values $V_+$ and $V_-$, which extend to $x = +\infty$ and $x = -\infty$ respectively. The features that the energy eigenstates have depend on how the energy $E$ is related to these three values. |
| |
| It is helpful to start by rewriting the TISE in the form | It is helpful to start by rewriting the TISE in the form |
| This is consistent with the intuition from classical mechanics, since the particle would have to have a negative kinetic energy everywhere to exist in this potential, which is not possible. | This is consistent with the intuition from classical mechanics, since the particle would have to have a negative kinetic energy everywhere to exist in this potential, which is not possible. |
| |
| ===== V_0 < E < V_+: Bound States ===== | ===== $V_0 < E < V_+$: Bound States ===== |
| |
| In this case, $k_-$ and $k_+$ are real and positive, but $k_0$ is imaginary. Define $k_0' = -ik_0 so that $k_0'= \sqrt{\frac{2m}{\hbar^2} \left ( E -V_0\right )}$ is real and positive. As before, we can discard the terms that diverge as $x\rightarrow\pm\infty$ and we are left with the solution | In this case, $k_-$ and $k_+$ are real and positive, but $k_0$ is imaginary. Define $k_0' = -ik_0$ so that $k_0'= \sqrt{\frac{2m}{\hbar^2} \left ( E -V_0\right )}$ is real and positive. As before, we can discard the terms that diverge as $x\rightarrow\pm\infty$ and we are left with the solution |
| |
| \[ \psi(x)=\begin{cases} \psi_-(x) = Ae^{k_-x} , & x<a \\ \psi_0(x)=Ce^{ik_0'x} + De^{-ik_0' x}, & a\leq x \leq b \\ \psi_+(x) = Ge^{-k_+x}, & x>b.\end{cases}\] | \[ \psi(x)=\begin{cases} \psi_-(x) = Ae^{k_-x} , & x<a \\ \psi_0(x)=Ce^{ik_0'x} + De^{-ik_0' x}, & a\leq x \leq b \\ \psi_+(x) = Ge^{-k_+x}, & x>b.\end{cases}\] |
| In summary, for bound states the eigenstates are normalizable and nondegenerate, and there is a discrete energy spectrum. | In summary, for bound states the eigenstates are normalizable and nondegenerate, and there is a discrete energy spectrum. |
| |
| ===== E>V_+: Scattering States ===== | ===== $E>V_+$: Scattering States ===== |
| |
| The cases $V_+ < E < V_-$ and $E>V_-$ are both called //scattering// states. There does not seem to be a standard terminology to distinguish these two cases, but I will call $V_+ < E < V_-$ //semi-scattering// states and $E>V_+$ //full-scattering// states. | The cases $V_+ < E < V_-$ and $E>V_-$ are both called //scattering// states. There does not seem to be a standard terminology to distinguish these two cases, but I will call $V_+ < E < V_-$ //semi-scattering// states and $E>V_+$ //full-scattering// states. |
| Since this solution oscillates as $x\rightarrow - \infty$, it will not be normalizable. We can still choose to set one of the five unknown constants arbitrarily, since multiplying a wavefunction by an overall constant does not change the physical state. As before, the boundary conditions give us four equations but, in this case, we have four unknown constants, so there is a unique solution with no further constraints on the energy $E$. Therefore, there will be a continuous energy spectrum: all possible values of $E$ satisfying $V_+ < E < V_-$ are allowed. Since the boundary conditions uniquely determine the four unknown constants, the solutions are nondegenerate. | Since this solution oscillates as $x\rightarrow - \infty$, it will not be normalizable. We can still choose to set one of the five unknown constants arbitrarily, since multiplying a wavefunction by an overall constant does not change the physical state. As before, the boundary conditions give us four equations but, in this case, we have four unknown constants, so there is a unique solution with no further constraints on the energy $E$. Therefore, there will be a continuous energy spectrum: all possible values of $E$ satisfying $V_+ < E < V_-$ are allowed. Since the boundary conditions uniquely determine the four unknown constants, the solutions are nondegenerate. |
| |
| This type of solution is called semi-scattering because the particle can escape to $x = -\infty$, but not $x=+\infty$. As we have seen, semi-scattering states are unnormalizable, nondegenerate, and have a continuous energy spectrum. | This type of solution is called semi-scattering because the particle can escape to $x = +\infty$, but not $x=-\infty$. As we have seen, semi-scattering states are unnormalizable, nondegenerate, and have a continuous energy spectrum. |
| |
| In the full scattering case, all of $k_-$, $k_0$ and $k_+$ are imaginary, so we replace them by the real positive constants $k_-'= \sqrt{\frac{2m}{\hbar^2} \left ( E -V_-\right )}$, $k_0'= \sqrt{\frac{2m}{\hbar^2} \left ( E -V_0\right )}$ and $k_+'= \sqrt{\frac{2m}{\hbar^2} \left ( E -V_+\right )}$. There are no divergent terms in the general solution, so we have | In the full scattering case, all of $k_-$, $k_0$ and $k_+$ are imaginary, so we replace them by the real positive constants $k_-'= \sqrt{\frac{2m}{\hbar^2} \left ( E -V_-\right )}$, $k_0'= \sqrt{\frac{2m}{\hbar^2} \left ( E -V_0\right )}$ and $k_+'= \sqrt{\frac{2m}{\hbar^2} \left ( E -V_+\right )}$. There are no divergent terms in the general solution, so we have |
| * Normalizable | * Normalizable |
| * Nondegenerate | * Nondegenerate |
| * Have a discrete energy spectrum</span> | * Have a discrete energy spectrum |
| * The solution will be a //semi-scattering state// if the region in which $E > V(x)$ extends to either $x\rightarrow +\infty$ or $x\rightarrow-\infty$ but not both. Semi-scattering states are: | * The solution will be a //semi-scattering state// if the region in which $E > V(x)$ extends to either $x\rightarrow +\infty$ or $x\rightarrow-\infty$ but not both. Semi-scattering states are: |
| * Unnormalizable | * Unnormalizable |
| A node of a wavefunction is a point at which it crosses the $x$-axis, from negative to positive or vice versa. In the solutions for an infinite square well potential | A node of a wavefunction is a point at which it crosses the $x$-axis, from negative to positive or vice versa. In the solutions for an infinite square well potential |
| |
| \[\psi(x)=\begin{cases} 0, x\leq 0 \\ \sqrt{\frac{2}{a}}\sin\left ( \frac{n\pi x}{a}\right ), & 0<x<a \\ 0, & x\geq a, \end{cases}\] | \[\psi(x)=\begin{cases} 0, & x\leq 0 \\ \sqrt{\frac{2}{a}}\sin\left ( \frac{n\pi x}{a}\right ), & 0<x<a \\ 0, & x\geq a, \end{cases}\] |
| |
| if we exclude the end points $x=0$ and $x=a$, the ground state has no nodes, the first excited state has one node, the second excited state has two nodes, etc. This behavior is generally true for the bound states of all potentials. A proof is beyond the scope of this course, but the intuition is the same as for the infinite square well. In the regions where $E > V(x)$, the solution is oscillatory. The more times the solution oscillates in this region the higher the energy. Roughly, more oscillations indicate a smaller "wavelength" and hence a larger momentum. "Wavelength" is in scare quotes here because the solution need not be periodic, so the concept does not strictly apply, but the general idea that a larger number of nodes gives a larger mean momentum, and hence energy, is correct. Therefore, the ground state will have the smallest possible number of nodes, i.e. zero, the first excited state the next smallest, i.e. one, and so on. | if we exclude the end points $x=0$ and $x=a$, the ground state has no nodes, the first excited state has one node, the second excited state has two nodes, etc. This behavior is generally true for the bound states of all potentials. A proof is beyond the scope of this course, but the intuition is the same as for the infinite square well. In the regions where $E > V(x)$, the solution is oscillatory. The more times the solution oscillates in this region the higher the energy. Roughly, more oscillations indicate a smaller "wavelength" and hence a larger momentum. "Wavelength" is in scare quotes here because the solution need not be periodic, so the concept does not strictly apply, but the general idea that a larger number of nodes gives a larger mean momentum, and hence energy, is correct. Therefore, the ground state will have the smallest possible number of nodes, i.e. zero, the first excited state the next smallest, i.e. one, and so on. |
| At the end of section 4.ii we considered the symmetric infinite potential well | At the end of section 4.ii we considered the symmetric infinite potential well |
| |
| \[V(x) = \begin{cases} \infty, & x<-\frac{a}{2}, \\ 0, & -\frac{a}{2} \leq x \leq +\frac{a}{2}, \\ \infty, x>+\frac{a}{2}, \end{cases}\] | \[V(x) = \begin{cases} \infty, & x<-\frac{a}{2}, \\ 0, & -\frac{a}{2} \leq x \leq +\frac{a}{2}, \\ \infty, & x>+\frac{a}{2}, \end{cases}\] |
| |
| as illustrated below. | as illustrated below. |
| We noted that the eigenstates | We noted that the eigenstates |
| |
| \[\psi_n(x) = \begin{cases} \sqrt{\frac{2}{a}} \cos \left ( \frac{n\pi x}{a}\right ), & n=1,3,5,7,\cdots \sqrt{\frac{2}{a}} \sin \left ( \frac{n\pi x}{a} \right ), & n=2,4,6,8,\cdots,\\ \end{cases}\] | \[\psi_n(x) = \begin{cases} \sqrt{\frac{2}{a}} \cos \left ( \frac{n\pi x}{a}\right ), & n=1,3,5,7,\cdots \\ \sqrt{\frac{2}{a}} \sin \left ( \frac{n\pi x}{a} \right ), & n=2,4,6,8,\cdots,\\ \end{cases}\] |
| |
| were a sequence of alternating odd and even functions. The ground state is odd, the first excited state is even, the second excited state is odd, etc. | were a sequence of alternating odd and even functions. The ground state is odd, the first excited state is even, the second excited state is odd, etc. |
| \[\hat{\mathcal{P}}|\psi\rangle = -\int_{+\infty}^{-\infty} \mathrm{d}x'\, \psi(-x') |-x'\rangle = \int_{-\infty}^{+\infty} \mathrm{d}x\, \psi(-x) |x\rangle, \] | \[\hat{\mathcal{P}}|\psi\rangle = -\int_{+\infty}^{-\infty} \mathrm{d}x'\, \psi(-x') |-x'\rangle = \int_{-\infty}^{+\infty} \mathrm{d}x\, \psi(-x) |x\rangle, \] |
| |
| and hence $\hat{\mathcal{P}}$ maps a wavefunction $\psi(x)$ to $\psi(-x)$. | where, in the second step we got rid of the minus sign by switching the limits and we replaced $x'$ by $x$ because $x'$ is a dummy variable. Hence, $\hat{\mathcal{P}}$ maps a wavefunction $\psi(x)$ to $\psi(-x)$. |
| |
| In the in class activities, you will show that $\hat{\mathcal{P}}$ is a Hermitian operator, which implies that it has real eigenvalues. Note also that | In the in class activities, you will show that $\hat{\mathcal{P}}$ is a Hermitian operator, which implies that it has real eigenvalues. Note also that |
| \[\hat{\mathcal{P}} \hat{p}^2 = -\hat{p} \hat{\mathcal{P}} \hat{p} = +\hat{p}^2 \hat{\mathcal{P}}.\] | \[\hat{\mathcal{P}} \hat{p}^2 = -\hat{p} \hat{\mathcal{P}} \hat{p} = +\hat{p}^2 \hat{\mathcal{P}}.\] |
| |
| Similarly, if the potential is even, $V(-\hat{x}) = V(\hat{x})$ the its Taylor expansion only includes even powers of $\hat{x}$, | Similarly, if the potential is even, $V(-\hat{x}) = V(\hat{x})$, its Taylor expansion only includes even powers of $\hat{x}$, |
| |
| \[V(\hat{x}) = a_0 + a_2 \hat{x}^2 + a_3 \hat{x}^4 + \cdots.\] | \[V(\hat{x}) = a_0 + a_2 \hat{x}^2 + a_3 \hat{x}^4 + \cdots.\] |
| The same type of argument works for the rest of the excited states. If they have an even number of nodes then there cannot be a node at $x=0$ and hence the wavefunction must be even. If there are an odd number of nodes then one of them must be at $x=0$ meaning that the wavefunction must be odd. Since the number of nodes increases by one each time we increase the energy to the next allowed value in the spectrum, we get an alternating sequence of even and odd wavefunctions. | The same type of argument works for the rest of the excited states. If they have an even number of nodes then there cannot be a node at $x=0$ and hence the wavefunction must be even. If there are an odd number of nodes then one of them must be at $x=0$ meaning that the wavefunction must be odd. Since the number of nodes increases by one each time we increase the energy to the next allowed value in the spectrum, we get an alternating sequence of even and odd wavefunctions. |
| |
| {{:question-mark.png?nolink&50 |}}====== In Class Activities ====== | {{:question-mark.png?nolink&50 |}} |
| | ====== In Class Activities ====== |
| |
| - Prove that the parity operator $\hat{\mathcal{P}}$ is Hermitian, i.e. | - Prove that the parity operator $\hat{\mathcal{P}}$ is Hermitian, i.e. |
| HINT: Use the fact that $\hat{\mathcal{P}}$ maps $\psi(x)$ to $\psi(-x)$ and the definition of the inner product | HINT: Use the fact that $\hat{\mathcal{P}}$ maps $\psi(x)$ to $\psi(-x)$ and the definition of the inner product |
| \[\langle \phi | \psi \rangle = \int_{-\infty}^{+\infty} \mathrm{d}x\, \phi^*(x) \psi(x).\] | \[\langle \phi | \psi \rangle = \int_{-\infty}^{+\infty} \mathrm{d}x\, \phi^*(x) \psi(x).\] |
| - Show that $\hat{x}$ is an odd operator. It is sufficient to show that $\hat{x}\hat{\mathcal{P}}|x\rangle = -\hat{\mathcal{P}}\hat{x}|x\rangle$. Why? | - |
| - Show that $\hat{\mathcal{P}}|p\rangle = |-p\rangle$ using $\langle x | p \rangle = \frac{1}{\sqrt{2\pi \hbar}} e^{ipx/\hbar}$. | - Show that $\hat{x}$ is an odd operator. It is sufficient to show that $\hat{x}\hat{\mathcal{P}}|x\rangle = - |
| - Show that $\hat{p}$ is an odd operator. | \hat{\mathcal{P}}\hat{x}|x\rangle$. Why? |
| | - Show that $\hat{\mathcal{P}}|p\rangle = |-p\rangle$ using $\langle x | p \rangle = \frac{1}{\sqrt{2\pi \hbar}} e^{ipx/\hbar}$. |
| | - Show that $\hat{p}$ is an odd operator. |
| |
| |