Tuesday, 6 October 2020

quantum mechanics - How do I determine the location of a free particle with Schrödinger's equation?


I'm trying to get to grips with the Schrödinger equation by looking at a free particle. I'm certain at some point I massively misunderstood something.


According to a textbook and a lecture the free particle moving in positive x direction can be described by


$$ \Psi(x,t) = A e^{i(kx - \omega t)} = \psi(x) \cdot e^{-i\omega t} $$


Classically, I would expect the particle to move along its path with some constant velocity $v = (x-x_0)/t$, so I would like to determine the probability $P(x,\Delta x,t)$ to find the particle between $x$ and $\Delta x$ at time $t$ in order to compare it with the classical location $x(t) = x_0 + v\cdot t$.


Since I'm looking for a location, I have to use the (trivial) location operator $\hat{r} = r$ and I get:


$$ P(x,\Delta x,t) = \int_x^{x+\Delta x}\Psi(x,t)^* \hat{\Psi(x,t)} dx\\ = \int_x^{x+\Delta x}\Psi(x,t)^* \Psi(x,t) dx\\ = A^2\int_x^{x+\Delta x} {e^{i(kx - \omega t)}}^* {e^{i(kx - \omega t)}} dx \\ = A^2\int_x^{x+\Delta x} dx = A^2 \Delta x $$


Which doesn't make sense to me at all, since it doesn't depend on either $x$ or $t$. The particle is not at all locations along the x-axis with the same probability at all times; I would rather expect it to move along the axis with the velocity $v$ but regardless which constants I shove into $A$ by using normalisation constraints, $P$ will never depend on $t$. But according to my understanding it should.


Obviously, my understanding is wrong and/or I made some mistakes in my calculation. Where am I going wrong?



Answer




When you solve the Schrodinger equation for a free particle you get a family of solutions of the form $\Psi(x,t) = A e^{i(kx - \omega t)}$ and all superpositions of these functions. So just solving the Schrodinger equation doesn't give you a solution for a specific particle. For that you need to specify the initial conditions.


If you take the solution to be $\Psi(x,t) = A e^{i(kx - \omega t)}$ then you are (without realising it) specifying the initial condition to be a completely delocalised particle i.e. one for which we have precise knowledge of the momentum but no knowledge of the position. That's why when you attempt to calculate the position you get a silly answer.


If you specify the initial conditions as $\Psi(x, 0)$ then you have effectively created a wavepacket describing your particle, so it does have a finite uncertainty in position, and of course now a finite uncertainty in momentum. You can now calculate the expectation value of position as a function of time.


Your $\Psi(x, 0)$ will probably be expressed as a linear superposition of the plane wave solutions. To calculate the superposition just Fourier transform your $\Psi(x, 0)$.


Response to comment:


In your comment you ask:



I would have to plug in some constraints first and get a value for A or how would I go about doing that?



but it isn't just a matter of choosing the value of $A$ in $A e^{i(kx - \omega t)}$ because this doesn't describe a localised particle whatever value you choose for $A$.



Suppose at time $t = 0$ we know the position of the particle precisely, $x = 0$. This means the initial wavefunction is a delta function:


$$ \Psi(x, 0) = \delta(x) e^{-i\omega t} $$


i.e. $\Psi(x, 0)$ is zero everywhere except at $x = 0$. The position of this particle is obviously $x = 0$.


The trouble is that it isn't obvious how this wavefunction evolves in time. We know how the plane waves evolve in time, so we can easily calculate the time evolution if we could express the $\delta(x)$ function as a sum of plane waves:


$$ \delta(x) = \sum\limits_i A_i e^{i(kx - \omega t)} $$


The problem is working out how to do this sum, i.e. what are the values of the coefficients $a_i$ and how many terms we need in the sum. We can work this out by Fourier transforming our $\delta$ function, because this is exactly what a Fourier transform does. It expresses any function as an integral of plane waves. The Wikipedia article I've linked goes into more detail on this. The answer is that:


$$ \delta(x) = \frac{1}{2\pi} \int_{-\infty}^{\infty} e^{i(kx - \omega t)} dk $$


In fact choosing a $\delta$ function as the initial conditions isn't helpful because if we have an exact position we have infinite uncertainty in momentum, and if the momentum is infinitely uncertain we can't calculate the future position. If you're trying to describe a real system you would choose something like a Gaussian:


$$ \Psi(x, 0) = k \space e^{-(x^2/\Delta x^2)} $$


This describes a particle with the expectation value of $x = 0$ and the uncertainty in $x = \Delta x$. You can now Fourier transform your $\Psi(x, 0)$ to express it as an integral of plane waves, and you can now calculate the expectation value of $x$ as a function of time.



No comments:

Post a Comment

Understanding Stagnation point in pitot fluid

What is stagnation point in fluid mechanics. At the open end of the pitot tube the velocity of the fluid becomes zero.But that should result...