Probability a particle is in a certain region

AI Thread Summary
The discussion revolves around calculating the probability of a particle being in a certain region given a specific potential. The user initially sets up wave functions for two regions but suspects an error when applying boundary conditions at the interface. They explore the implications of the energy being less than zero and the necessity for continuity in the wave function's derivative. The conversation highlights the importance of correctly interpreting the conditions for sinusoidal and exponential solutions based on the potential's value. Ultimately, the user acknowledges a misunderstanding regarding the relationship between the wave function's parameters and the continuity requirements, suggesting a need for numerical solutions to resolve the issue.
etotheipi
Homework Statement
Given the particle is in its ground state with ##E <0##, calculate ##\mathbb{P}(x > \xi)## for the following potential:\begin{align*} V(x) =
\begin{cases}



\infty & x < 0 \\

- \lambda & 0\leq x\leq \xi \\

0 & x > \xi

\end{cases}

\end{align*}where ##\lambda > 0##.
Relevant Equations
N/A
I think I made an error somewhere. In ##[0,a]## I let ##\varphi(x) = \varphi_1(x) := p\sin{kx}## whilst in ##(\xi, \infty)## I let ##\varphi(x) = \varphi_2 (x) := re^{-\gamma x}##, and the constraints at ##x=\xi## are \begin{align*}

\varphi_1'(\xi) = \varphi_2'(\xi) &\implies pk\cos{k\xi} = -\gamma re^{-\gamma \xi} \\
\varphi_1(\xi) = \varphi_2(\xi) &\implies p\sin{k\xi} = re^{-\gamma \xi}

\end{align*}which implies ##\gamma = -k/\tan{k\xi}##. Then we can normalise,\begin{align*}

\int_0^\xi dx \, \varphi_1^2(x) + \int_\xi^{\infty} dx \, \varphi_2^2(x) &= \frac{(2k \xi - \sin{2 k \xi})p^2}{4k} + \frac{r^2 e^{-2\gamma \xi}}{2\gamma} \overset{!}{=} 1

\end{align*}and, since ##p^2 \sin^2{k \xi} = r^2 e^{-2\gamma \xi}## it follows that\begin{align*}

\mathbb{P}(x > \xi) = \frac{r^2 e^{-2\gamma \xi}}{2\gamma} = \frac{2\sin^2{k \xi} \tan{k \xi}}{2\sin^2{k \xi} \tan{k \xi} + \sin{2k \xi} - 2k \xi}\end{align*}That the particle is in its ground state implies ##k \xi = \pi / 2##, such that the sinusoidal part joins nicely onto the decaying exponential solution in the ##x > \xi## region. But that gives a non-sense answer; so I wonder where I made my mistake(s). Thanks!
 
Last edited by a moderator:
Physics news on Phys.org
etotheipi said:
Homework Statement:: Given the particle is in its ground state with ##E <0##, calculate ##\mathbb{P}(x > \xi)## for the following potential:\begin{align*} V(x) =
\begin{cases}
\infty & x < 0 \\

- \lambda & 0\leq x\leq \xi \\

0 & x > \xi

\end{cases}

\end{align*}where ##\lambda > 0##.
Relevant Equations:: N/A

I think I made an error somewhere. In ##[0,a]## I let ##\varphi(x) = \varphi_1(x) := p\sin{kx}## whilst in ##(\xi, \infty)## I let ##\varphi(x) = \varphi_2 (x) := re^{-\gamma x}##, and the constraints at ##x=\xi## are \begin{align*}

\varphi_1'(\xi) = \varphi_2'(\xi) &\implies pk\cos{k\xi} = -\gamma re^{-\gamma \xi} \\
\varphi_1(\xi) = \varphi_2(\xi) &\implies p\sin{k\xi} = re^{-\gamma \xi}

\end{align*}which implies ##\gamma = -k/\tan{k\xi}##. Then we can normalise,\begin{align*}

\int_0^\xi dx \, \varphi_1^2(x) + \int_\xi^{\infty} dx \, \varphi_2^2(x) &= \frac{(2k \xi - \sin{2 k \xi})p^2}{4k} + \frac{r^2 e^{-2\gamma \xi}}{2\gamma} \overset{!}{=} 1

\end{align*}and, since ##p^2 \sin^2{k \xi} = r^2 e^{-2\gamma \xi}## it follows that\begin{align*}

\mathbb{P}(x > \xi) = \frac{r^2 e^{-2\gamma \xi}}{2\gamma} = \frac{2\sin^2{k \xi} \tan{k \xi}}{2\sin^2{k \xi} \tan{k \xi} + \sin{2k \xi} - 2k \xi}\end{align*}That the particle is in its ground state implies ##k \xi = \pi / 2##, such that the sinusoidal part joins nicely onto the decaying exponential solution in the ##x > \xi## region. But that gives a non-sense answer; so I wonder where I made my mistake(s). Thanks!
##k \xi = \pi / 2## implies ##\varphi_1'(\xi) =0##, no? Which is a bit of a problem for ##\varphi_2'(\xi)## with the chosen exponential form.
Not knowing anything much about QM, I do not know how you are able to assume those forms for the two regions.
 
haruspex said:
Not knowing anything much about QM, I do not know how you are able to assume those forms for the two regions.
They're the solutions to the Schrödinger equation that satisfy the boundary conditions at ##x=0## and ##x \to \infty##.
 
vela said:
They're the solutions to the Schrödinger equation that satisfy the boundary conditions at ##x=0## and ##x \to \infty##.
Ok, I've done a little reading..
The sign of E-V(x) is critical, yes? Where positive, we should get a trig solution and where negative an exponential one.
But V is given as ##-\lambda## in the central range, and we are told E<0 so it should be exponential solutions for all x>0?
 
haruspex said:
##k \xi = \pi / 2## implies ##\varphi_1'(\xi) =0##, no? Which is a bit of a problem for ##\varphi_2'(\xi)## with the chosen exponential form.
Yeah, that's the non-sense bit that I can't seem to figure out 😄. That seems to be the false assumption (i.e. I believe we are instead only allowed to write ##k\xi > \pi/2##), but the question remains to then find the actual value of ##k \xi## in the ground state...

haruspex said:
Ok, I've done a little reading..
The sign of E-V(x) is critical, yes? Where positive, we should get a trig solution and where negative an exponential one.
But V is given as ##-\lambda## in the central range, and we are told E<0 so it should be exponential solutions for all x>0?
How I read it was that ##- \lambda < E < 0##, so ##E > - \lambda## means you have sinusoidal solutions in the central region and ##E < 0## means you have a decaying exponential in the ##x > \xi## region. The question doesn't explicitly state that, but it's the only way the situation seems to make sense in the first place (as far as I can tell... ☺️)
 
Last edited by a moderator:
etotheipi said:
How I read it was that ##- \lambda < E < 0##, so ##E > - \lambda## means you have sinusoidal solutions in the central region and ##E < 0## means you have a decaying exponential in the ##x > \xi## region. The question doesn't explicitly state that, but it's the only way the situation seems to make sense in the first place (as far as I can tell... ☺️)
Just realized I tripped over the double negative.
##E-V=E+\lambda##, so it is unclear which sign it will have.
You are saying we should take it as positive... ok.
That leaves two things to question:
- does the derivative need to be continuous? I note that in the infinite well we have zero outside the well and ##\sin(nx\pi/L) ## inside. That does not have a continuous derivative at the boundaries.
- why should ##k\xi=\pi/2##?
 
  • Like
Likes etotheipi
haruspex said:
That leaves two things to question:
- does the derivative need to be continuous? I note that in the infinite well we have zero outside the well and ##\sin(nx\pi/L) ## inside. That does not have a continuous derivative at the boundaries.
Infinite potential barriers are a bit weird; from the Schroedinger equation it's possible to show that solutions have continuous derivatives so long as ##V(x)## does not have any infinite discontinuities. So whilst the discontinuity in ##\varphi'## at ##x=0## can be tolerated, we still need to impose continuity in ##\varphi'## everywhere else.

haruspex said:
- why should ##k\xi=\pi/2##?
I think that is my mistake. I edited my previous post #5 just a short while ago to reflect that the actual statement should be ##k \xi > \pi/2##; that's because it's only possible to satisfy the continuity conditions at ##x = \xi## if the gradient of the wave-function is negative. And without loss of generality we can take ##k## to be positive (because the wave-function itself is not physically meaningful, only it's [modulus] square).
 
etotheipi said:
Yeah, that's the non-sense bit that I can't seem to figure out 😄. That seems to be the false assumption (i.e. I believe we are instead only allowed to write ##k\xi > \pi/2##), but the question remains to then find the actual value of ##k \xi## in the ground state...
You have to solve it numerically.

You can also look at the system graphically. Let ##k_0^2 = 2m\lambda/\hbar^2##. Then you have
$$k^2 = \frac{2m(E+\lambda)}{\hbar^2} = -\gamma^2 + k_0^2.$$ So the solutions to the system correspond to the intersection of ##\gamma = -k \cot k\xi## with the circle ##\gamma^2 + k^2 = k_0^2## in the ##k\gamma## plane.
 
  • Like
Likes etotheipi
Back
Top