This post further builds on the rather remarkable results we got in our previous posts. Let us start with the basics once again. The elementary wavefunction is written as:

ψ = *a·e*^{−i[E·t − p∙x]/ħ} = *a·cos*(**p**∙**x**/ħ − E∙t/ħ) *+** i·a·sin*(**p**∙**x**/ħ − E∙t/ħ)

Of course, *Nature* (or God, as Einstein would put it) does not care about our conventions for measuring an angle (i.e. the phase of our wavefunction) clockwise or counterclockwise and, therefore, the ψ = *a·e*^{i}^{[E·t − p∙x]/ħ} function is also permitted. We know that *cos*(θ) = *cos*(−θ) and *sin*θ = −*sin*(*−*θ), so we can write: * *

ψ = *a·e*^{i}^{[E·t − p∙x]/ħ} = *a·cos*(E∙t/ħ − **p**∙**x**/ħ) *+** i·a·sin*(E∙t/ħ − **p**∙**x**/ħ)

*= **a·cos*(**p**∙**x**/ħ − E∙t/ħ) −* i·a·sin*(**p**∙**x**/ħ − E∙t/ħ)

The vectors **p** and **x** are the momentum and position vector respectively: **p** = (p_{x}, p_{y}, p_{z}) and **x** = (x, y, z). However, if we assume there is *no* uncertainty about **p** – not about the direction, and not about the magnitude – then the direction of **p** can be our x-axis. In this reference frame, **x** = (x, y, z) reduces to (x, 0, 0), and **p**∙**x**/ħ reduces to p∙x/ħ. This amounts to saying our particle is traveling along the x-axis or, if p = 0, that our particle is located somewhere on the x-axis. So we have an analysis in one dimension only then, which facilitates our calculations. The geometry of the wavefunction is then as illustrated below. The x-axis is the direction of propagation, and the y- and z-axes represent the real and imaginary part of the wavefunction respectively.

Note that, when applying the right-hand rule for the axes, the vertical axis is the y-axis, not the z-axis. Hence, we may associate the vertical axis with the cosine component, and the horizontal axis with the sine component. [You can check this as follows: if the origin is the (x, t) = (0, 0) point, then *cos*(θ) = *cos*(0) = 1 and *sin*(θ) = *sin*(0) = 0. This is reflected in both illustrations, which show a left- and a right-handed wave respectively.]

Now, you will remember that we speculated the two *polarizations* (left- versus right-handed) should correspond to the two possible values for the quantum-mechanical spin of the wave (+ħ/2 or −ħ/2). We will come back to this at the end of this post. Just focus on the essentials first: the cosine and sine components for the left-handed wave are shown below. Look at it carefully and try to understand. Needless to say, the cosine and sine function are the same, except for a phase difference of π/2: *sin*(θ) = *cos*(θ − π/2).

As for the wave velocity, and its direction of propagation, we know that the (phase) velocity of any waveform F(kx − ωt) is given by *v*_{p} = ω/k. In our case, we find that *v*_{p} = ω/k = (E/ħ)/(p/ħ) = E/p. Of course, the momentum might also be in the negative x-direction, in which case k would be equal to −p and, therefore, we would get a negative phase velocity: *v*_{p} = ω/k = (E/ħ)/(−p/ħ) = −E/p.

As you know, E/ħ = ω gives the *frequency in time* (expressed in radians per second), while p/ħ = k gives us the wavenumber, or the *frequency in space* (expressed in radians per meter). [If in doubt, check my post on essential wave math.] Now, you also know that *f* = ω/2π and λ = 2π/k, which gives us the two de Broglie relations:

- E = ħ∙ω = h∙
*f*
- p = ħ∙k = h/λ

The frequency in time (oscillations or radians per second) is easy to interpret. A particle will always have some mass and, therefore, some energy, and it is easy to appreciate the fact that the wavefunction of a particle with more energy (or more mass) will have a *higher density in time* than a particle with less energy.

However, the second *de Broglie *relation is somewhat harder to interpret. Note that the wavelength is *inversely *proportional to the momentum: λ = h/p. Hence, if p goes to zero, then the wavelength becomes infinitely long, so we write:

If p → 0 then λ* → ∞.*

For the limit situation, a particle with zero *rest *mass (m_{0} = 0), the velocity may be *c* and, therefore, we find that p = m_{v}∙*v* = m_{c}∙*c* = m∙*c *(all of the energy is kinetic) and, therefore, p∙*c *= m∙*c*^{2} = E, which we may also write as: E/p = *c*. Hence, for a particle with zero rest mass (m_{0} = 0), the wavelength can be written as:

λ = h/p = h*c*/E = h/m*c*

Of course, we are talking a *photon *here. We get the zero rest mass for a photon. In contrast, all *matter*-particles should have *some *mass[1] and, therefore, their velocity will *never* equal *c*.[2] The question remains: how should we interpret the inverse proportionality between *λ* and p?

Let us first see what this wavelength λ actually represents. If we look at the ψ = a·*cos*(p∙x/ħ − E∙t/ħ) − *i·a·sin*(p∙x/ħ – E∙t/ħ) once more, and if we write p∙x/ħ as Δ, then we can look at p∙x/ħ as a phase factor, and so we will be interested to know for what x this phase factor Δ = p∙x/ħ will be equal to 2π. So we write:

Δ =p∙x/ħ = 2π ⇔ x = 2π∙ħ/p = h/p = λ

So now we get a meaningful interpretation for that wavelength. It is the distance between the crests (or the troughs) of the wave, so to speak, as illustrated below. Of course, this two-dimensional wave has no real crests or troughs: they depend on your frame of reference.

So now we know what λ actually represent for our one-dimensional elementary wavefunction. Now, the time that is needed for one cycle is equal to T = 1/*f *= 2π·(ħ/E). Hence, we can now calculate the wave velocity:

*v* = λ/T = (h/p)/[2π·(ħ/E)] = E/p

Unsurprisingly, we just get the phase velocity that we had calculated already: *v *= *v*_{p} = E/p. It does not answer the question: what if p is zero? What if we are looking at some particle at rest? It is an intriguing question: we get an infinitely long wavelength, and an infinite phase velocity. Now, we know *phase *velocities can be superluminal, but they should not be infinite. So what does the mathematical inconsistency tell us? Do these infinitely long wavelengths and infinite wave velocities tell us that our particle *has to* move? Do they tell us our notion of a particle at rest is mathematically inconsistent?

Maybe. But maybe not. Perhaps the inconsistency just tells us our elementary wavefunction – or the concept of a *precise *energy, and a *precise *momentum – does not make sense. This is where the Uncertainty Principle comes in: stating that p = 0, implies zero uncertainty. Hence, the σ_{p} factor in the σ_{p}∙σ_{x} ≤ ħ/2 would be zero and, therefore, σ_{p}∙σ_{x} would be zero which, according to the Uncertainty Principle, it cannot be: it can be very small, but it cannot be zero.

It is interesting to note here that σ_{p} refers to the standard deviation from the mean, as illustrated below. Of course, the distribution may be or may not be normal – we don’t know – but a normal distribution makes intuitive sense, of course. Also, if we assume the mean is zero, then the uncertainty is basically about the *direction *in which our particle is moving, as the momentum might then be positive *or *negative.

The question of natural units may pop up. The Uncertainty Principle suggests a *numerical *value of the natural unit for momentum and distance that is equal to the *square root *of ħ/2, so that’s about 0.726×10^{−17} m for the distance unit and 0.726×10^{−17} N∙s for the momentum unit, as the product of both gives us ħ/2. To make this somewhat more real, we may note that 0.726×10^{−17} m is the attometer scale (1 am = 1×10^{−18} m), so that is very small but not unreasonably small.[3]

Hence, we need to superimpose a potentially infinite number of waves with energies and momenta centered on some mean value. It is only then that we get meaningful results. For example, the idea of a *group *velocity – which should correspond to the classical idea of the velocity of our particle – only makes sense in the context of wave *packet*. Indeed, the group velocity of a wave packet (*v*_{g}) is calculated as follows:

*v*_{g} = ∂ω_{i}/∂k_{i} = ∂(E_{i}/ħ)/∂(p_{i}/ħ) = ∂(E_{i})/∂(p_{i})

This assumes the existence of a *dispersion relation* which gives us ω_{i} as a function of k_{i} – what amounts to the same – E_{i} as a function of p_{i}. How do we get that? Well… There are a few ways to go about it but one interesting way of doing it is to re-write Schrödinger’s equation as the following *pair* of equations[4]:

*Re*(∂ψ/∂t) = −[ħ/(2m_{eff})]·*Im*(∇^{2}ψ) ⇔ ω·cos(kx − ωt) = k^{2}·[ħ/(2m_{eff})]·cos(kx − ωt)
*Im*(∂ψ/∂t) = [ħ/(2m_{eff})]·*Re*(∇^{2}ψ) ⇔ ω·sin(kx − ωt) = k^{2}·[ħ/(2m_{eff})]·sin(kx − ωt)

These equations imply the following dispersion relation:

ω = ħ·k^{2}/(2m)

Of course, we need to think about the subscripts now: we have ω_{i}, k_{i}, but… What about m_{eff} or, dropping the subscript, about m? Do we write it as m_{i}? If so, what is it? Well… It is the *equivalent* mass of E_{i} obviously, and so we get it from the mass-energy equivalence relation: m_{i} = E_{i}/*c*^{2}. It is a fine point, but one most people forget about: they usually just write m. However, if there is uncertainty in the energy, then Einstein’s mass-energy relation tells us we must have some uncertainty in the (equivalent) mass too, and the two will, obviously, be related as follows: σ_{m} = σ_{E}/*c*^{2}. We are tempted to do a few substitutions here. Let’s first check what we get when doing the m_{i} = E_{i}/*c*^{2} substitution:

ω_{i} = ħ·k_{i}^{2}/(2m_{i}) = (1/2)∙ħ·k_{i}^{2}∙*c*^{2}/E_{i} = (1/2)∙ħ·k_{i}^{2}∙*c*^{2}/(ω_{i}∙ħ) = (1/2)∙ħ·k_{i}^{2}∙*c*^{2}/ω_{i}

⇔ ω_{i}^{2}/k_{i}^{2} = *c*^{2}/2 ⇔ ω_{i}/k_{i} = *v*_{p} = *c*/2 !?

We get a very interesting but nonsensical *condition* for the dispersion relation here. I wonder what mistake I made. 😦

Let us try another substitution. The group velocity is what it is, right? It is the velocity of the group, so we can write: k_{i} = p/ħ = m_{i} ·*v*_{g}. This gives us the following result:

ω_{i} = ħ·(m_{i} ·*v*_{g})^{2}/(2m_{i}) = ħ·m_{i}·*v*_{g}^{2}/2

It is yet another interesting *condition *for the dispersion relation. Does it make any more sense? I am not so sure. That factor 1/2 troubles us. It only makes sense when we *drop *it. Now you will object that Schrödinger’s equation gives us the electron orbitals – and many other correct descriptions of quantum-mechanical stuff – so, surely, Schrödinger’s equation cannot be wrong. You’re right. It’s just that… Well… When we are splitting in up in two equations, as we are doing here, then we are looking at *one *of the two dimensions of the oscillation only and, therefore, it’s only *half *of the mass that counts. Complicated explanation but… Well… It should make sense, because the results that come out make sense. Think of it. So we write this:

*Re*(∂ψ/∂t) = −(ħ/m_{eff})·*Im*(∇^{2}ψ) ⇔ ω·cos(kx − ωt) = k^{2}·(ħ/m_{eff})·cos(kx − ωt)
*Im*(∂ψ/∂t) = (ħ/m_{eff})·*Re*(∇^{2}ψ) ⇔ ω·sin(kx − ωt) = k^{2}·(ħ/m_{eff})·sin(kx − ωt)

We then get the dispersion relation *without *that 1/2 factor:

ω_{i} = ħ·k_{i}^{2}/m_{i}

The m_{i} = E_{i}/*c*^{2} substitution then gives us the result we sort of expected to see:

ω_{i} = ħ·k_{i}^{2}/m_{i} = ħ·k_{i}^{2}∙*c*^{2}/E_{i} = ħ·k_{i}^{2}∙*c*^{2}/(ω_{i}∙ħ) ⇔ ω_{i}/k_{i} = *v*_{p} = *c*

Likewise, the other calculation also looks more meaningful now:

ω_{i} = ħ·(m_{i} ·*v*_{g})^{2}/m_{i} = ħ·m_{i}·*v*_{g}^{2}

Sweet ! 🙂

Let us put this aside for the moment and focus on something else. If you look at the illustrations above, you see we can sort of distinguish (1) a linear velocity – the speed with which those wave crests (or troughs) move – and (2) some kind of circular or tangential velocity – the velocity along the red contour line above. We’ll need the formula for a tangential velocity: *v*_{t} = *a*∙ω.

Now, if λ is zero, then *v*_{t} = *a*∙ω = *a*∙E/ħ is just all there is. We may double-check this as follows: the distance traveled in one period will be equal to 2π*a*, and the period of the oscillation is T = 2π·(ħ/E). Therefore, *v*_{t} will, effectively, be equal to *v*_{t} = 2π*a*/(2πħ/E) = *a*∙E/ħ. However, if λ is non-zero, then the distance traveled in one period will be equal to 2π*a *+ λ. The period remains the same: T = 2π·(ħ/E). Hence, we can write:

For an electron, we did this weird calculation. We had an angular momentum formula (for an electron) which we equated with the *real-life* +ħ/2 or −ħ/2 values of its spin, and we got a *numerical *value for *a*. It was the Compton radius: the scattering radius for an electron. Let us write it out:

Using the right numbers, you’ll find the *numerical* value for *a*: 3.8616×10^{−13} m. But let us just substitute the formula itself here:

This is fascinating ! And we just calculated that *v*_{p} is equal to *c*. For the elementary wavefunction, that is. Hence, we get this amazing result:

*v*_{t} = 2*c*

This *tangential *velocity is *twice *the *linear *velocity !

Of course, the question is: what is the *physical *significance of this? I need to further look at this. Wave velocities are, essentially, *mathematical *concepts only: the wave propagates through space, but *nothing else *is really moving. However, the geometric implications are obviously quite interesting and, hence, need further exploration.

One conclusion stands out: all these results reinforce our interpretation of the speed of light as a *property *of the vacuum – or of the fabric of spacetime itself. 🙂

[1] Even neutrinos should have some (rest) mass. In fact, the mass of the known neutrino flavors was estimated to be smaller than 0.12 eV/*c*^{2}. This mass combines the three known neutrino flavors.

[2] Using the Lorentz factor (γ), we can write the relativistically correct formula for the kinetic energy as KE = E − E_{0} = m_{v}*c*^{2} − m_{0}*c*^{2} = m_{0}γ*c*^{2} − m_{0}*c*^{2} = m_{0}*c*^{2}(γ − 1). As *v *approaches *c*, γ approaches infinity and, therefore, the kinetic energy would become infinite as well.

[3] It is, of course, extremely small, but 1 am is the current sensitivity of the LIGO detector for gravitational waves. It is also thought of as the upper limit for the length of an electron, for quarks, and for fundamental strings in string theory. It is, in any case, 1,000,000,000,000,000,000 times larger than the order of magnitude of the Planck length (1.616229(38)×10^{−35} m).

[4] The m_{eff} is the *effective* mass of the particle, which depends on the medium. For example, an electron traveling in a solid (a transistor, for example) will have a different effective mass than in an atom. In free space, we can drop the subscript and just write m_{eff} = m. As for the equations, they are easily derived from noting that two complex numbers a + *i*∙b and c + *i*∙d are equal if, and only if, their real and imaginary parts are the same. Now, the ∂ψ/∂t = *i*∙(ħ/m_{eff})∙∇^{2}ψ equation amounts to writing something like this: a + *i*∙b = *i*∙(c + *i*∙d). Now, remembering that *i*^{2} = −1, you can easily figure out that *i*∙(c + *i*∙d) = *i*∙c + *i*^{2}∙d = − d + *i*∙c.