# Classical Wavefunctions

Right now I'm on the 4th floor of the physics building. If I walk down the hall to the balcony overlooking the atrium, I could drop a bouncy ball and watch its trajectory. It'll fall to the ground, bounce up to some fraction of its initial height, and repeat the process with a loss of energy each time until it finally comes to rest. In a perfect world with no dissipation of energy, the ball would always bounce exactly up to its initial height and repeat its trajectory forever.

In principle, I can drop the magic ball at a given time on the stopwatch and thenceforth predict exactly where it will be at any point in the future. On the other hand, I might just choose to tell you that I have in fact made the drop from the 4th floor without telling you when exactly that drop was. If you passed your physics classes this semester, you could still calculate things like the bounce rate, impact velocities, and what have you. But you couldn't say exactly where the ball is right now, because I haven't given you the time of drop.

On the other hand, you do know that the ball is moving slower when it's near the top of its trajectory and faster when it's near the bottom. Thus it spends more time near the top and less time near the ground. You could, if you wanted, calculate the probability that the ball is near some given height h at a randomly chosen instant. In some sense, this is what a quantum mechanical wavefunction is*. We calculated the quantum wavefunction for just this problem not too long ago, and we can compare it to the classical probability distribution for the classical situation I've described.

We'll do the comparison later. For now we'll just calculate the "classical wavefunction". We know that the probability P(x) that the ball will be found near a height x is proportional to how long it spends around x. How long it spends around x is inversely proportional to its velocity at x, because if it's going fast it's not going to be around long. So P(x) ~ 1/v(x). We should find v(x). From conservation of energy, we know that the total energy (kinetic and potential) of the ball is a constant. Call it E.

Remember that when I dropped the ball, its initial velocity is 0. This means that the actual numerical value of E is mgh, where h is the height from which I dropped it. Anyway, we can solve the above equation for v:

Since the probability is proportional to 1/v, we have:

But this isn't quite a probability density. For one thing the units are wrong, for another it doesn't integrate to 1. We need to find the normalization constant C such that:

Hey, it's exam time. Work that out yourselves for practice. Most of the constants will cancel, and we'll use some theorist-style simplification and define φ = x/h while picking units such that h = 1. In that case, the probability distribution is:

We can graph this:

From the graph, it's clear that you're a lot more likely to find the ball near the top - just as we expected. In fact, there's about a 71% chance that you'll find it above the halfway mark, and you're more than 4 times more likely to find it in the top 1/5th than you are in the bottom 1/5th.

But this doesn't look a whole lot like the Airy function, which is the quantum mechanical result. We'll take a closer look at the reason next time.

*In other very important senses, the two situations are very different. For instance, one's a stationary state and one isn't. The quantum analogue of the classical problem is really something more like a coherent state rather than a high-n eigenvalue.

Tags

### More like this

Last week, Rhett did a post on animating a bouncing ball in VPython. This was mostly making a point about the distinction between real simulation and animation, along the lines of yesterday's post on social construction of videogame reality. But, of course, my immediate reaction was, "That's not…
##### Clearing the air on the Airy fuction
There was some dissension in the comments of my post on solving the Schrodinger equation with a linear potential. What the post boiled down to was that the solution was Ai(u), where we found that u was: The point of the post was to work through and get that coefficient that's in front of the (Fx…
##### Energy Dissipation in a Physics Toy
A little while back, I used a photo of SteelyKid's toy Newton's cradle as the photo of the day, with a bonus video: I mentioned that I was going to do some analysis of this at some point, but didn't have time right then. I had a bit of time to poke at this yesterday, though, so here's some physics…
##### Quantum Bouncing Ball, Concluded
Well, last time we were looking about the classical probability density for a bouncing ball, and the quantum mechanical probability distribution for the same. They looked not even a little bit alike. This is kind of a problem, since we know from experience that classical physics works pretty well,…

I'm interested in seeing how far you go with this calculation.

I actually calculated the (approximate) coherent states for the quantum bouncer during my 4th year at uni. It's fairly special in that like the harmonic oscillator, the phase space trajectory of (, ) match the classical system quite well. This does not happen in other systems such as the infinite square well or the Coulomb problem.

I've planning on writing up a quick demonstrations.wolfram for it, when I get some time....

Damn html... it should be the expectation value of x and p inside those parentheses.

um... amusing spelling fail in title.

I assigned the probability distribution for a classical pendulum in my modern physics class this semester, and in looking for lecture material came across this interesting paper:

Am. J. Phys. 63 (9), 823, 1995: R.W. Robinett "Quantum and classical probability distributions for position and momentum"

Covers the harmonic oscillator, infinite well, and linear confining potential, as well as a couple of unbound systems.

If the probability is inversely related to the velocity, doesn't that mean that the probability of it being at the top of its arc is infinite?

Re: Paul

The pdf will approach infinity as the ball reaches the top of its arc, but the probability that it is exactly at the top at any given point in time is 0 since the distribution is continuous.

That's actually the density function not the distribution function.

A more direct derivation is as follows -

The random variable is f(t) = h -(1/2) g t^2

So for 0<=x<=h

P(x) = Prob({f(t) < x }) = Prob ({sqrt(2(h-x)/g < t })
= 1 - sqrt(2(h-x)/g)/sqrt(2h/g) = 1 - sqrt(1 - x/h)

So the distribution function is

0 x<=0
1 - sqrt(1-x/h) 0<=x<=h
1 h<=x

Then the density function is

0 x<=0
1/(2 h sqrt(1-x/h)) 0<=x<=h
0 h<=x

By Annonymous (not verified) on 15 Dec 2010 #permalink

My comment got garbled. The inequality signs show as
equal signs. The distribution function is 1 - sqrt(1-x/h)
for x between 0 and h, 0 for x less than 0 and 1 for x
greater than h.

The density function is 1/(2 h sqrt(1-x/h) ) for x between
0 and h and otherwise 0

By Annonymous (not verified) on 15 Dec 2010 #permalink

Some of the derivation got lost. The point is
that the distribution function for the random
variable h - (1/2) g t^2 comes immediately from
the definition of the distribution function using
just algebra. Then the density function is obtained by
differentiation.

By Annonymous (not verified) on 15 Dec 2010 #permalink

take your meds, you little fckers....

now we are going to bury you...

And the lesson from all of this? DOUBLE!
____________________________

What do you want, you little ****ers?

more of these idiots