The geometry of the wavefunction (2)

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·ei[E·t − px]/ħa·cos(px/ħ − E∙t/ħ) + i·a·sin(px/ħ − 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·ei[E·t − px]/ħ function is also permitted. We know that cos(θ) = cos(−θ) and sinθ = −sin(θ), so we can write:    

ψ = a·ei[E·t − p∙x]/ħa·cos(E∙t/ħ − px/ħ) + i·a·sin(E∙t/ħ − px/ħ)

= a·cos(px/ħ − E∙t/ħ) − i·a·sin(px/ħ − E∙t/ħ)

The vectors p and x are the momentum and position vector respectively: p = (px, py, pz) 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 px/ħ 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).

circular polarizaton with components

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 vp = ω/k. In our case, we find that vp = ω/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: vp = ω/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:

  1. E = ħ∙ω = h∙f
  2. 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 (m0 = 0), the velocity may be c and, therefore, we find that p = mvv = mcc = m∙c (all of the energy is kinetic) and, therefore, p∙c = m∙c2 = E, which we may also write as: E/p = c. Hence, for a particle with zero rest mass (m0 = 0), the wavelength can be written as:

λ = h/p = hc/E = h/mc

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.

wavelength

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 = vp = 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.

Standard_deviation_diagram

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 (vg) is calculated as follows:

vg = ∂ωi/∂ki = ∂(Ei/ħ)/∂(pi/ħ) = ∂(Ei)/∂(pi)

This assumes the existence of a dispersion relation which gives us ωi as a function of ki – what amounts to the same – Ei as a function of pi. 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]:

  1. Re(∂ψ/∂t) = −[ħ/(2meff)]·Im(∇2ψ) ⇔ ω·cos(kx − ωt) = k2·[ħ/(2meff)]·cos(kx − ωt)
  2. Im(∂ψ/∂t) = [ħ/(2meff)]·Re(∇2ψ) ⇔ ω·sin(kx − ωt) = k2·[ħ/(2meff)]·sin(kx − ωt)

These equations imply the following dispersion relation:

ω = ħ·k2/(2m)

Of course, we need to think about the subscripts now: we have ωi, ki, but… What about meff or, dropping the subscript, about m? Do we write it as mi? If so, what is it? Well… It is the equivalent mass of Ei obviously, and so we get it from the mass-energy equivalence relation: mi = Ei/c2. 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/c2. We are tempted to do a few substitutions here. Let’s first check what we get when doing the mi = Ei/c2 substitution:

ωi = ħ·ki2/(2mi) = (1/2)∙ħ·ki2c2/Ei = (1/2)∙ħ·ki2c2/(ωi∙ħ) = (1/2)∙ħ·ki2c2i

⇔ ωi2/ki2 = c2/2 ⇔ ωi/ki = vp = 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: ki = p/ħ = mi ·vg. This gives us the following result:

ωi = ħ·(mi ·vg)2/(2mi) = ħ·mi·vg2/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) = −(ħ/meffIm(∇2ψ) ⇔ ω·cos(kx − ωt) = k2·(ħ/meff)·cos(kx − ωt)
  • Im(∂ψ/∂t) = (ħ/meffRe(∇2ψ) ⇔ ω·sin(kx − ωt) = k2·(ħ/meff)·sin(kx − ωt)

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

ωi = ħ·ki2/mi

The mi = Ei/c2 substitution then gives us the result we sort of expected to see:

ωi = ħ·ki2/mi = ħ·ki2c2/Ei = ħ·ki2c2/(ωi∙ħ) ⇔ ωi/ki = vp = c

Likewise, the other calculation also looks more meaningful now:

ωi = ħ·(mi ·vg)2/mi = ħ·mi·vg2

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: vt = a∙ω.

Now, if λ is zero, then vt = 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, vt will, effectively, be equal to vt = 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:

F1

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:

F2

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: F3

This is fascinating ! And we just calculated that vp is equal to c. For the elementary wavefunction, that is. Hence, we get this amazing result:

vt = 2c

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/c2. 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 − E0 = mvc2 − m0c2 = m0γc2 − m0c2 = m0c2(γ − 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 meff 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 meff = 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∙(ħ/meff)∙∇2ψ equation amounts to writing something like this: a + i∙b = i∙(c + i∙d). Now, remembering that i2 = −1, you can easily figure out that i∙(c + i∙d) = i∙c + i2∙d = − d + i∙c.

Advertisements

The geometry of the wavefunction

My posts and article on the wavefunction as a gravitational wave are rather short on the exact geometry of the wavefunction, so let us explore that a bit here. By now, you know the formula for the elementary wavefunction by heart:

ψ = a·ei[E·t − px]/ħa·cos(px/ħ − E∙t/ħ) + i·a·sin(px/ħ − E∙t/ħ)

If we assume the momentum p is all in the x-direction, then the p and x vectors will have the same direction, and px/ħ reduces to p∙x/ħ. This amounts to saying our particle is traveling along the x-axis. The geometry of the wavefunction is 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. 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. I am convinced these correspond to the two possible values for the quantum-mechanical spin of the wave: +ħ/2 or −ħ/2. But… Well… Who am I? The cosine and sine components are shown below. Needless to say, the cosine and sine function are the same, except for a phase difference of π/2: sin(θ) = cos(θ − π/2)  circular polarizaton with components

Surely, Nature doesn’t care a hoot about our conventions for measuring the phase angle clockwise or counterclockwise and therefore, the ψ = a·ei[E·t − px]/ħ function should, effectively, also be permitted. We know that cos(θ) = cos(θ) and sinθ = sin(θ), so we can write:    

ψ = a·ei[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/ħ)

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). Of course, we may write: f = ω/2π  and λ = 2π/k, which gives us the two de Broglie relations:

  1. E = ħ∙ω = h∙f
  2. p = ħ∙k = h/λ

The frequency in time is easy to interpret (a particle will always have some mass and, therefore, some energy), but the wavelength is inversely proportional to the momentum: λ = h/p. Hence, if p goes to zero, then the wavelength becomes infinitely long: if p → 0, then λ → ∞. For the limit situation, a particle with zero rest mass (m0 = 0), the velocity may be c and, therefore, we find that p = mvv = m∙c  and, therefore, p∙c = m∙c2 = E, which we may also write as: E/p = c. Hence, for a particle with zero rest mass, the wavelength can be written as:

λ = h/p = hc/E = h/mc

However, we argued that the physical dimension of the components of the wavefunction may be usefully expressed in N/kg units (force per unit mass), while the physical dimension of the electromagnetic wave are expressed in N/C (force per unit charge). This, in fact, explains the dichotomy between bosons (photons) and fermions (spin-1/2 particles). Hence, all matter-particles should have some mass.[1] But how we interpret the inverse proportionality between λ and p?

We should probably first ask ourselves what wavelength we are talking about. The wave only has a phase velocity here, which is equal to vp = ω/k = (E/ħ)/(p/ħ) = E/p. Of course, we know that, classically, the momentum will be equal to the group velocity times the mass: p = m·vg. However, when p is zero, we have a division by zero once more: if p → 0, then vp = E/p → ∞. Infinite wavelengths and infinite phase velocities probably tell us that our particle has to move: our notion of a particle at rest is mathematically inconsistent. If we associate this elementary wavefunction with some particle, and if we then imagine it to move, somehow, then we get an interesting relation between the group and the phase velocity:

vp = ω/k = E/p = E/(m·vg) = (m·c2)/(m·vg) = c2/vg

We can re-write this as vp·vg = c2, which reminds us of the relationship between the electric and magnetic constant (1/ε0)·(1/μ0) = c2. But what is the group velocity of the elementary wavefunction? Is it a meaningful concept?

The phase velocity is just the ratio of ω/k. In contrast, the group velocity is the derivative of ω with respect to k. So we need to write ω as a function of k. Can we do that even if we have only one wave? We do not have a wave packet here, right? Just some hypothetical building block of a real-life wavefunction, right? Right. So we should introduce uncertainty about E and p and build up the wave packet, right? Well… Yes. But let’s wait with that, and see how far we can go in our interpretation of this elementary wavefunction. Let’s first get that ω = ω(k) relation. You’ll remember we can write Schrödinger’s equation – the equation that describes the propagation mechanism for matter-waves – as the following pair of equations:

  1. Re(∂ψ/∂t) = −[ħ/(2m)]·Im(∇2ψ) ⇔ ω·cos(kx − ωt) = k2·[ħ/(2m)]·cos(kx − ωt)
  2. Im(∂ψ/∂t) = [ħ/(2m)]·Re(∇2ψ) ⇔ ω·sin(kx − ωt) = k2·[ħ/(2m)]·sin(kx − ωt)

This tells us that ω = ħ·k2/(2m). Therefore, we can calculate ∂ω/∂k as:

∂ω/∂k = ħ·k/m = p/m = vg

We learn nothing new. We are going round and round in circles here, and we always end up with a tautology: as soon as we have a non-zero momentum, we have a mathematical formula for the group velocity – but we don’t know what it represents – and a finite wavelength. In fact, using the p = ħ∙k = h/λ relation, we can write one as a function of the other:

λ = h/p = h/mvg ⇔ vg = h/mλ

What does this mean? It resembles the c = h/mλ relation we had for a particle with zero rest mass. Of course, it does: the λ = h/mc relation is, once again, a limit for vg going to c. By the way, it is interesting to note that the vp·vg = c2 relation implies that the phase velocity is always superluminal. That’ easy to see when you re-write the equation in terms of relative velocities: (vp/c)·(vg/c) = βphase·βgroup = 1. Hence, if βgroup < 1, then βphase > 1.

So what is the geometry, really? Let’s look at the ψ = a·cos(p∙x/ħ – E∙t/ħ) i·a·sin(p∙x/ħ – E∙t/ħ) formula once more. If we write p∙x/ħ as Δ, then we will be interested to know for what x this phase factor 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’s that distance between the crests of the wave, so to speak, as illustrated below.

wavelength

Can we now find a meaningful (i.e. geometric) interpretation for the group and phase velocity? If you look at the illustration above, you see we can sort of distinguish (1) a linear velocity (the speed with which those wave crests move) and (2) some kind of circular or tangential velocity (the velocity along the red contour line above). We’ll probably need the formula for the tangential velocity: v = a∙ω. If p = 0 (so we have that weird infinitesimally long wavelength), then we have two velocities:

  1. The tangential velocity around the a·ei·E·t  circle, so to speak, and that will just be equal to v = a∙ω = a∙E/ħ.
  2. The red contour line sort of gets stretched out, like infinitely long, and the velocity becomes… What does it do? Does it go to ∞ , or to c?

Let’s think about this. For a particle at rest, we had 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 so we got a numerical value for a. It was the Compton radius: the scattering radius for an electron. Let me copy it once again:

Compton radius formula

Just to bring this story a bit back to Earth, you should note the calculated value: = 3.8616×10−13 m. We did then another weird calculation. We said all of the energy of the electron had to be packed in this cylinder that might of might not be there. The point was: the energy is finite, so that elementary wavefunction cannot have an infinite length in space. Indeed, assuming that the energy was distributed uniformly, we jotted down this formula, which reflects the formula for the volume of a cylinder:

E = π·a2·l ⇔ = E/(π·a2)

Using the value we got for the Compton scattering radius (= 3.8616×10−13 m), we got an astronomical value for l. Let me write it out:

= (8.19×10−14)/(π·14.9×10−26) ≈ 0.175×1012 m

It is, literally, an astronomical value: 0.175×1012 m is 175 million kilometer, so that’s like the distance between the Sun and the Earth. We wrote, jokingly, that such space is too large to look for an electron and, hence, that we should really build a proper packet by making use of the Uncertainty Principle: allowing for uncertainty in the energy should, effectively, reduce the uncertainty in position.

But… Well… What if we use that value as the value for λ? We’d get that linear velocity, right? Let’s try it. The period is equal to T = T = 2π·(ħ/E) = h/E and λ = E/(π·a2), so we write:formula for vWe can write this as a function of m and the and ħ constants only:velocitiy 2

A weird formula but not necessarily nonsensical: we get a finite speed. Now, if the wavelength becomes somewhat less astronomical, we’ll get different values of course. I have a strange feeling that, with these formula, we should, somehow, be able to explain relativistic length contraction. But I will let you think about that as for now. Here I just wanted to show the geometry of the wavefunction a bit more in detail.

[1] The discussions on the mass of neutrinos are interesting in this regard. Scientists all felt the neutrino had to have some (rest) mass, so my instinct on this is theirs. In fact, only recently experimental confirmation came in, and the mass of the known neutrino flavors was estimated to be something like 0.12 eV/c2. This mass combines the three known neutrino flavors. To understand this number, you should note it is the same order of magnitude of the equivalent mass of low-energy photons, like infrared or microwave radiation.