Quantum Mechanics – How to Treat 1D Infinite Square Well When Box Size Increases Suddenly

homework-and-exercisesquantum mechanicsquantum-chemistryschroedinger equationwavefunction

I am currently working my way through John S. Townsend book "A Fundamental Approach to Modern Physics" (ISBN: 978-1-891389-62-7). Exercise 3.12 (p.111) is about the 1D infinite square well. The box has the potential barriers at $x=0$ and $x=L$.
$$
V\left(x\right) = \begin{cases} \infty & x < 0 \\ 0 & 0 \leq x \leq L \\ \infty & x > L \end{cases}
$$

The text states the following:

A particle of mass $m$ is in the lowest energy (ground) state of the infinite potential energy well. At time $t=0$ the wall located at $x=L$ is suddenly pulled back to a position at $x=2L$. This change occurs so rapidly that instantaneously the wave function does not change. ($a$) Calculate the probability that a measurement of the energy will yield the ground-state energy of the new well.

The answer to this question involves re-solving the T.I.S.E., applying new boundary conditions, and re-normalizing the wave function. Then we project the new wave-function onto the old one, and calculate the overlap of the two wave-functions (initial and final wave-function).

I cannot find any discussion about this particular problem that would make me able to reason my way to calculate an overlap between the initial and final wave-function. No way. Unless I have a background in linear algebra, and would be able to foresee the need to project the one onto the other, and sort of find my way logically to the answer, I see no way any student without any experience with quantum physics would be able to solve this on their own. I had to use Google until I found the solution to a similar problem.

Or, maybe QM just does not come naturally to me. A lot of problem-solving is required in these types of courses to for us to build up some "intuition" about what to expect as results, what to do, etc. I still do not fully understand the logic behind the solution to the problem.

Phenomenologically, there is a particle inside the box. Then, suddenly, the box widens to twice its size. However, the wave-function does not change. Even though the particle has a bigger box in which to move around, the initial wave-function does not take this extra room into account, and would equal zero when $x>L$. Hence, I would not expect for the particle to venture out to the new region much, if we used only the initial wave-function. Or, is the problem that the particle will indeed venture out there (since $V(x)=0$, and it is physically allowed to), but we have not considered this in our normalizing of the function? Hence, it makes no sense to calculate any property of the particle with the "initial" wave-function, since this is simply the incorrect wave-function for the new well?

New ideas and thoughts

Okay, so I have to evaluate the integral

$$
c_1 = \int\limits_{- \infty}^{\infty} \psi_1 (x) \Psi (x) dx
$$

where $\psi_1$ is the theoretical wave-function for a particle in the entire box, that is from $0$ to $2L$, and $\Psi$ is the wave-function for the actual particle, that is from $0$ to $L$.

Now, I can split this integral into two parts:

$$
c_1 = \int\limits_{- \infty}^{\infty} \psi_1 (x) \Psi (x) dx = \int\limits_{0}^{L} \psi_1 (x) \Psi (x) dx + \int\limits_{L}^{2L} \psi_1 (x) \Psi (x) dx
$$

We see that since the actual particle's wave-function is not defined when $x>L$, the second term will be zero – that is, the wave-functino is normalized for $0<x<L$, and the boundary conditions ensures that $\Psi \rightarrow 0$ when $ x=0$ and $x=L$. Thus, the the integral over is reduced to

\begin{align*}
c_1 &= \int\limits_{0}^{L} \psi_1 (x) \Psi (x) dx \\
&= \int\limits_{0}^{L} \sqrt{\frac{1}{L}} \sin \left( \frac{\pi x}{2L} \right) \sqrt{\frac{2}{L}} \sin \left( \frac{\pi x}{L} \right) dx \\
&= \frac{\sqrt{2}}{L} \int\limits_{0}^{L} \sin \left( \frac{\pi x}{2L} \right) \sin \left( \frac{\pi x}{L} \right) dx \\
&= \frac{4 \sqrt{2}}{3 \pi} \\
\Rightarrow c_1^2 &= \frac{32}{9 \pi ^2} \\
c_1^2 & \approx 0.36
\end{align*}

I think maybe my difficulty with "visualzing" the problem is that I do not fully understand the expression for $c_n^2$, and how this gives the probability for measuring the energy level $n$.

I know this seems like I am rambling, and maybe I am, but I hope you understand my confusion. I appreciate any help!

Best Answer

...it makes no sense to calculate any property of the particle with the "initial" wave-function, since this is simply the incorrect wave-function for the new well?

The wavefunction can't be "incorrect for the well". Your wavefunction is just an initial condition for time-dependent Schrödinger equation. Here's how it would evolve if you solve the time-dependent equation (I ignore normalization here):

enter image description here

We see that since the actual particle's wave-function is not defined when $x>L$...

The wavefunction is defined for all $x\in\mathbb R$. It's just zero outside of $(0,L)$ because the potential is infinite there.

the second term will be zero

This remains true, however, because of what I said above.

the wave-functino is normalized for $0<x<L$

Actually, again, the wavefunction is normalized period. It's defined for the whole real line, and zero outside the well, so when you normalized using the integral over the well, it's the same as if you integrated over $\mathbb R$. If it were not this way, your "normalization for some domain" would not make any sense, i.e. it'd not be a normalization at all.

Your further calculation looks OK to me.

I think maybe my difficulty with "visualzing" the problem is that I do not fully understand the expression for $c^2_n$, and how this gives the probability for measuring the energy level $n$.

That this gives the probability for measuring energy level $n$ is known as Born rule. You find projection of your actual wave function onto the eigenstate of energy, namely on state $n$. By Born rule, its magnitude squared is the probability of measuring the system to appear in that eigenstate.

The fact that you used only the unchanged original wavefunction for calculations, even though it almost immediately drastically changes as time goes on, is that despite its change in shape, its coefficients $c_n$ actually only change their phase like $$c_n\propto\exp\left(-\frac i\hbar E_nt\right),$$

but remain the same in magnitude — because the potential is independent of time. Thus you can measure the energy after some time passes, and will still get the same result.

Related Question