Let’s play a bit with the stuff we found in our previous post. This is going to be unconventional, or *experimental*, if you want. The idea is to give you… Well… Some ideas. So you can play yourself. 🙂 Let’s go.

Let’s first look at Feynman’s (simplified) formula for the amplitude of a photon to go from point *a* to point *b*. If we identify point *a *by the position vector ** r_{1}** and point

*b*by the position vector

**, and using Dirac’s fancy**

*r*_{2}*bra-ket*notation, then it’s written as:

So we have a vector dot product here: **p**∙*r*_{12} = |**p**|∙|*r*_{12}|· *cos*θ = p∙*r*_{12}·*cos*α. The angle here (α) is the angle between the **p **and *r*_{12} vector. All good. Well… No. We’ve got a problem. When it comes to calculating *probabilities*, the α angle doesn’t matter: |*e*^{i·θ}/*r*|^{2} = 1/*r*^{2}. Hence, for the probability, we get: P = | 〈**r**_{2}|**r**_{1}〉 |^{2} = 1/*r*_{12}^{2}. ** Always !** Now that’s strange. The θ =

**p**∙

*r*_{12}/

*ħ*argument gives us a different phase depending on the angle (α) between

**p**and

*r*_{12}. But… Well… Think of it:

*cos*α goes from 1 to 0 when α goes from 0 to ±90° and, of course, is

*negative*when

**p**and

*r*_{12}have opposite directions but… Well… According to this formula, the

*probabilities*do

*not*depend on the direction of the momentum. That’s just weird, I think. Did Feynman, in his iconic

*Lectures*, give us a meaningless formula?

Maybe. We may also note this function looks like the elementary wavefunction for any particle, which we wrote as:

ψ(** x**,

*t*) =

*a·e*

^{−i}^{∙θ}=

*a·e*

^{−i}^{∙}

^{(E∙t − p∙x)/ħ}=

*a*·

*e*

^{−i}

^{∙}

^{(E∙t)/ħ}·

*e*

^{i}^{∙}

^{(p∙x)/ħ}

The only difference is that the 〈**r**_{2}|**r**_{1}〉 sort of abstracts away from time, so… Well… Let’s get a feel for the quantities. Let’s think of a photon *carrying *some typical amount of energy. Hence, let’s talk visible light and, therefore, photons of a few eV only – say 5.625 eV = 5.625×1.6×10^{−19} J = 9×10^{−19} J. Hence, their momentum is equal to p = E/*c* = (9×10^{−19} N·m)/(3×10^{5} m/s) = 3×10^{−24} N·s. That’s tiny but that’s only because *newtons *and *seconds *are enormous units at the (sub-)atomic scale. As for the distance, we may want to use the thickness of a playing card as a starter, as that’s what Young used when establishing the *experimental *fact of light interfering with itself. Now, playing cards in Young’s time were obviously rougher than those today, but let’s take the smaller distance: modern cards are as thin as 0.3 mm. Still, that distance is associated with a value of θ that is equal to 13.6 *million*. Hence, the *density *of our wavefunction is enormous at this scale, and it’s a bit of a miracle that Young could see any interference at all ! As shown in the table below, we only get meaningful values (remember: θ is a phase *angle*) when we go down to the *nanometer* scale (10^{−9} m) or, even better, the *angstroms *scale ((10^{−9} m).* *

So… Well… Again: what can we do with Feynman’s formula? Perhaps he didn’t give us a *propagator* function but something that is more general (read: more *meaningful*) at our (limited) level of knowledge. As I’ve been reading Feynman for quite a while now – like three or four years 🙂 – I think… Well… Yes. That’s it. Feynman wants us to *think* about it. 🙂 *Are you joking again, Mr. Feynman? *🙂 So let’s assume the reasonable thing: let’s assume it gives us the amplitude to go from point *a *to point *b *by the position vector ** along some path r**. So, then, in line with what we wrote in our previous post, let’s say p·

*r*(momentum over a distance) is the action (

*S*) we’d associate with this particular path (

*r*) and then see where we get. So let’s write the formula like this:

ψ = *a*·*e*^{i·θ} = (1/*r*)·*e*^{i·S/ħ} = *e*^{i·p∙r/ħ}/*r*

We’ll use an index to denote the various paths: *r*_{0} is the straight-line path and *r*_{i} is any (other) path. Now, quantum mechanics tells us we should calculate this amplitude *for every possible path*. The illustration below shows the straight-line path and two nearby paths. So each of these paths is associated with some amount of action, which we measure in *Planck units*: θ = *S*/*ħ*.

The time interval is given by *t *= *t*_{0 }= *r*_{0}/*c*, *for all paths*. Why is the time interval the same *for all paths*? Because we think of a photon going from some *specific *point in space *and in time* to some other *specific *point in space *and in time*. Indeed, when everything is said and done, we do think of light as traveling from point *a *to point *b *at the speed of light (*c*). In fact, all of the weird stuff here is all about trying to explain *how *it does that. 🙂

Now, if we would think of the photon *actually *traveling along this or that path, then this implies its velocity along any of the nonlinear paths will be *larger *than *c*, which is OK. That’s just the weirdness of quantum mechanics, and you should actually *not *think of the photon actually traveling along one of these paths anyway although we’ll often put it that way. Think of something fuzzier, whatever that may be. 🙂

So the action is energy times time, or momentum times distance. Hence, the difference in action between two paths *i *and *j* is given by:

δ*S *= p·*r*_{j} − p·*r*_{i} = p·(*r*_{j} − *r*_{i}) = p·Δ*r*

I’ll explain the δ*S < *2π*ħ*/3 thing in a moment. Let’s first pause and think about the *uncertainty *and how we’re modeling it. We can effectively think of the variation in *S *as some uncertainty in the action: δ*S *= Δ*S *= p·Δ*r*. However, if *S* is also equal to energy times time (*S *= E·*t*), and we insist *t *is the same for all paths, then we must have some uncertainty in the energy, right? Hence, we can write δ*S *as Δ*S *= ΔE·*t*. But, of course, E = E = *m*·*c*^{2} = p·*c*, so we will have an uncertainty in the momentum as well. Hence, the variation in *S *should be written as:

δ*S *= Δ*S* = Δp·Δ*r*

That’s just logical thinking: if we, somehow, entertain the idea of a photon going from some *specific *point in spacetime to some other *specific *point in spacetime along various paths, then the variation, or *uncertainty*, in the action will effectively combine some uncertainty in the momentum and the distance. We can calculate Δp as ΔE/*c*, so we get the following:

δ*S *= Δ*S* = Δp·Δ*r = *ΔE·Δ*r*/*c = *ΔE·Δ*t* with Δ*t** = *Δ*r*/*c*

So we have the two expressions for the Uncertainty Principle here: Δ*S* = Δp·Δ*r* = ΔE·Δ*t*. Just be careful with the interpretation of Δ*t*: it’s just the equivalent of Δ*r*. We just express the uncertainty in distance in *seconds *using the (absolute) speed of light. We are *not *changing our spacetime interval: we’re still looking at a photon going from *a *to *b *in *t *seconds, *exactly*. Let’s now look at the δ*S < *2π*ħ*/3 thing. If we’re adding *two *amplitudes (two *arrows* or *vectors*, so to speak) and we want the magnitude of the result to be larger than the magnitude of the two contributions, then the angle between them should be smaller than 120 degrees, so that’s 2π/3 *rad*. The illustration below shows how you can figure that out geometrically.Hence, if *S*_{0} is the action for *r*_{0}, then *S*_{1} = *S*_{0} + *ħ *and *S*_{2} = *S*_{0} + 2·*ħ *are still good, but *S*_{3} = *S*_{0} + 3·*ħ* is *not *good. Why? Because the difference in the phase angles is Δθ = *S*_{1}/*ħ* − *S*_{0}/*ħ* = (*S*_{0} + *ħ*)/*ħ* − *S*_{0}/*ħ* = 1 and Δθ = *S*_{2}/*ħ* − *S*_{0}/*ħ* = (*S*_{0} + 2·*ħ*)/*ħ* − *S*_{0}/*ħ* = 2 respectively, so that’s 57.3° and 114.6° respectively and that’s, effectively, *less *than 120°. In contrast, for the next path, we find that Δθ = *S*_{3}/*ħ* − *S*_{0}/*ħ* = (*S*_{0} + 3·*ħ*)/*ħ* − *S*_{0}/*ħ* = 3, so that’s 171.9°. So that amplitude gives us a *negative *contribution.

Let’s do some calculations using a spreadsheet. To simplify things, we will assume we measure everything (time, distance, force, mass, energy, action,…) in Planck units. Hence, we can simply write: *S _{n}* =

*S*

_{0}+

*n*. Of course,

*n*= 1, 2,… etcetera, right? Well… Maybe not. We are

*measuring*action in units of

*ħ*, but do we actually think action

*comes*in units of

*ħ*? I am not sure. It would make sense, intuitively, but… Well… There’s uncertainty on the energy (E) and the momentum (

*) of our photon, right? And how accurately can we measure the distance? So there’s some randomness everywhere. 😦 So let’s leave that question open as for now.*

**p**We will also assume that the phase angle for *S*_{0} is equal to 0 (or some multiple of 2π, if you want). That’s just a matter of choosing the origin of time. This makes it really easy: Δ*S*_{n} = *S _{n}* −

*S*

_{0}=

*n*, and the associated phase angle θ

_{n}= Δθ

_{n}is the same. In short, the amplitude for each path reduces to ψ

_{n}=

*e*

^{i·n}/

*r*

_{0}. So we need to add these

*first*and

*then*calculate the magnitude, which we can then square to get a probability. Of course, there is also the issue of normalization (probabilities have to add up to one) but let’s tackle that later. For the calculations, we use Euler’s

*r*·

*e*

^{i·θ}=

*r*·(

*cos*θ +

*i*·

*sin*θ) =

*r*·

*cos*θ +

*i*·

*r*·

*sin*θ formula. Needless to say, |

*r*·

*e*

^{i·θ}|

^{2}= |

*r|*

^{2}·

*|*

*e*

^{i·θ}|

^{2}

*=*|

*r*|

^{2}·(

*cos*

^{2}θ +

*sin*

^{2}θ) =

*r*. Finally, when adding complex numbers, we add the real and imaginary parts respectively, and we’ll denote the ψ

_{0}+ ψ

_{1}+ψ

_{2}+ … sum as Ψ.

Now, we also need to see how our Δ*S* = Δp·Δ*r* works out. We may want to assume that the uncertainty in p and in *r *will both be proportional to the overall uncertainty in the action. For example, we could try writing the following: Δ*S*_{n} = Δp* _{n}*·Δ

*r*=

_{n}*n*·Δp

_{1}·Δ

*r*

_{1}. It also makes sense that you may want Δp

*and Δ*

_{n}*r*to be proportional to Δp

_{n}_{1}and Δ

*r*

_{1}respectively. Combining both, the assumption would be this:

Δp* _{n}* = √

*n*·Δp

_{1 }and Δ

*r*= √

_{n}*n*·Δ

*r*

_{1}

So now we just need to decide how we will distribute Δ*S*_{1} = *ħ *= 1 over Δp_{1} and Δ*r*_{1} respectively. For example, if we’d assume Δp_{1} = 1, then Δ*r*_{1} = *ħ*/Δp_{1} = 1/1 = 1. These are the calculations. I will let you analyze them. 🙂Well… We get a weird result. It reminds me of Feynman’s explanation of the partial reflection of light, shown below, but… Well… That doesn’t make much sense, does it?

Hmm… Maybe it does. 🙂 Look at the graph more carefully. The *peaks *sort of oscillate out so… Well… That might make sense… 🙂

Does it? Are we doing something wrong* *here? These amplitudes should reflect the ones that are reflected in those nice animations (like this one, for example, which is part of that’s part of the Wikipedia article on Feynman’s path integral formulation of quantum mechanics). So what’s wrong, if anything? Well… Our paths differ by some fixed amount of action, which doesn’t quite reflect the geometric approach that’s used in those animations. The graph below shows how the distance *r *varies as a function of *n*.

If we’d use a model in which the distance would *increase* linearly or, preferably, exponentially, then we’d get the result we want to get, right?

Well… Maybe. Let’s try it. Hmm… We need to think about the geometry here. Look at the triangle below. If *b *is the straight-line path (*r*_{0}), then *ac *could be one of the crooked paths (*r _{n}*). To simplify, we’ll assume

*isosceles*triangles, so

*a*equals

*c*and, hence,

*r*= 2·

_{n}*a*= 2·

*c*. We will also assume the successive paths are separated by the same vertical distance (

*h = h*

_{1}) right in the middle, so

*h*=

_{b}*h*

_{n}=

*n*·

*h*

_{1}. It is then easy to show the following:This gives the following graph for

*r*= 10 and

_{n}*h*

_{1 }= 0.01.

Is this the right step increase? Not sure. We can vary the values in our spreadsheet. Let’s first build it. The photon will have to travel *faster* in order to cover the extra distance in the same time, so its momentum will be higher. Let’s think about the velocity. Let’s start with the first path (*n *= 1). In order to cover the *extra *distance Δ*r*_{1}, the velocity *c*_{1} must be equal to (*r*_{0} + Δ*r*_{1})/*t *= *r*_{0}/*t *+ Δ*r*_{1}/*t =* *c *+ Δ*r*_{1}/*t *= *c*_{0}* *+ Δ*r*_{1}/*t*. We can write *c*_{1} as *c*_{1} = *c*_{0}* *+ Δ*c*_{1}, so Δ*c*_{1} = Δ*r*_{1}/*t**.* Now, the *ratio *of p_{1} and p_{0} will be equal to the *ratio *of *c*_{1} and *c*_{0} because p_{1}/p_{0 }= (m*c*_{1})/m*c*_{0}) = *c*_{1}/*c*_{0}. Hence, we have the following formula for p_{1}:

p_{1} = p_{0}·*c*_{1}/*c*_{0} = p_{0}·(*c*_{0}* *+ Δ*c*_{1})/*c*_{0} = p_{0}·[1 + Δ*r*_{1}/(*c*_{0}*·t*) = p_{0}·(1 + Δ*r*_{1}/*r*_{0})

For p* _{n}*, the logic is the same, so we write:

p* _{n}* = p

_{0}·

*c*/

_{n}*c*

_{0}= p

_{0}·(

*c*

_{0}

*+ Δ*

*c*)/

_{n}*c*

_{0}= p

_{0}·[1 + Δ

*r*/(

_{n}*c*

_{0}

*·t*) = p

_{0}·(1 + Δ

*r*/

_{n}*r*

_{0})

Let’s do the calculations, and let’s use meaningful values, so the nanometer scale and actual values for Planck’s constant and the photon momentum. The results are shown below.

Pretty interesting. In fact, this looks *really *good. The *probability *first swings around wildly, because of these zones of constructive and destructive interference, but then stabilizes. [Of course, I would need to normalize the probabilities, but you get the idea, right?] So… Well… I think we get a *very* meaningful result with this model. Sweet ! 🙂 ** I’m lovin’ it !** 🙂 And, here you go, this is (part of) the calculation table, so you can see what I am doing. 🙂

The graphs below look even better: I just changed the *h*_{1}/*r*_{0} ratio from 1/100 to 1/10. The probability stabilizes almost immediately. 🙂 So… Well… It’s not as fancy as the referenced animation, but I think the educational value of this thing here is at least as good ! 🙂

🙂 This is good stuff… 🙂

**Post scriptum **(19 September 2017): There is an obvious inconsistency in the model above, and in the calculations. We assume there is a path r_{1} = , *r*_{2}, *r*_{2},etcetera, and then we calculate the action for it, and the amplitude, and then we add the amplitude to the sum. But, surely, we should count these paths *twice*, in two-dimensional space, that is. Think of the graph: we have positive and negative interference zones that are sort of layered around the straight-line path, as shown below.

In three-dimensional space, these lines become surfaces. Hence, rather than adding *one* arrow for every δ* * having *one *contribution only, we may want to add… Well… In three-dimensional space, the formula for the surface around the straight-line path would probably look like π·*h*_{n}·*r*_{1}, right? Hmm… Interesting idea. I changed my spreadsheet to incorporate that idea, and I got the graph below. It’s a nonsensical result, because the probability does swing around, but it gradually spins out of control: it never stabilizes.That’s because we increase the *weight *of the paths that are further removed from the center. So… Well… We shouldn’t be doing that, I guess. 🙂 I’ll you look for the right formula, OK? Let me know when you found it. 🙂