Square shape wave packet spreading

[youngurlee]

for a 1D free particle with initial wave function \phi(x') square shaped(e.g. \phi(x')=1,x'\in [a,b],otherwise it vanishes),
my question is: how does it evolve with time t?


if we deal with it in P basis, it is easily solved, using the propagator U(t)=∫|p'><p'|e^{-\frac{ip'^2 t}{2m\hbar}}dp';

but if we directly solve SE in X basis, where P must be written as -i\frac{∂}{∂x'}, the initial wavefunction is not continous, so the equation becomes improper at the ends of the interval[a,b],


so why dose the SE equation seems so distinct in these 2 representations? what goes wrong in X representation?

 

[youngurlee]

my question comes thus:

suppose we set up a device to detect a particle, it can detect the particle when the particle occur with in the region [a,b],
so when the device really detected the particle, the wave function of the particle must collaspe to one vanishes without the region, in a special case, to a square shaped one.

I want to know its latter evolution, but I got confused when I tried to solve the SE directly in the Schrödinger representation .

Thanks for any tip.

 

[vanhees71]

You can simply use the solution in momentum space and transform it to position space. That's also nicely done numerically if you can't find an analytic expression. In fact you have

\psi(t,x)=\langle x|\psi(t) \rangle=\int \mathrm{d} p \langle x | p \rangle \langle p|\psi> = \int \mathrm{d} p \frac{\exp(\mathrm{i} p x)}{\sqrt{2 \pi}} \tilde{\psi}(t,p).

Now you have

\tilde{\psi}(t,p)=\langle p |\exp[-\mathrm{i} \hat{p}^2/(2m)] \psi(t=0) \rangle = \exp[-\mathrm{i} p^2/(2m)] \tilde{\psi}_0(p).

The wave function at t=0 in the momentum-space representation is of course

\tilde{\psi}_0(p)=\int \mathrm{d} x \frac{\exp(-\mathrm{i} p x)}{(2 \pi)^{1/2}} \psi_0(x).

Now you have just to plug everything together.

 

[youngurlee]

vanhees71, thank you for your reply.

Your method is exactly what I mean by solving it using the propagator U(t) in terms of P's eigenbras and eigenkets.

What I wonder is why can't I solve it in coordinate space.

For example, after infinitesimal interval \Delta t, by -i\hbar\frac{\partial}{\partial t}|\psi\rangle=H|\psi\rangle, we have
\Delta|\psi\rangle=\frac{\Delta t}{i\hbar}H|\psi\rangle.

if we work in P basis, that's no problem, since then
\Delta\psi(p)=\frac{p^2\Delta t}{i 2m\hbar}\psi(p), and \psi(p) behaves very well.

however, if we work in X basis,
\Delta\psi(x)=i\hbar\Delta t\frac{\partial^2}{\partial t^2}\psi(x), but look at how \frac{\partial^2}{\partial t^2}\psi(x) behaves here, it just becomes improper.

with this perspective, how could these two methods lead to the same result?