This post is intended to help you think about, and work with, those mysterious amplitudes. More in particular, I’ll explore how *potential *differences change amplitudes. But let’s first recapitulate the basics.

In my previous post, I explained why the young French *Comte* Louis de Broglie, when writing his PhD thesis back in 1924, i.e. before Schrödinger, Born, Heisenberg and others had published their work, boldly proposed to the ω·t − **k**·**x** argument in the wavefunction of a particle with the relativistic invariant product of the momentum and position four-vectors p_{μ }= (E, **p**) = (E, p_{x}, p_{y}, p_{z},) and x_{μ }= (t, **x**) = (t, x, y, z), provided the energy and momentum are re-scaled in terms of ħ. Hence, he wrote:

θ = ω·t − **k**·**x** = (p_{μ}x_{μ})/ħ = (E∙t − **p**∙**x**)/ħ = (E/ħ)∙t − (**p**/ħ)∙**x**

As it’s usually instructive to do a quick dimensional analysis, let’s do one here too. Energy is expressed in joule, and dividing it by the quantum of action, which is expressed in joule·seconds (J·s) gives us the dimension of an (angular) frequency indeed, which, in turns, yields a pure number. Likewise, linear momentum can be expressed in newton·seconds which, when divided by joule·seconds (J·s), yields a quantity expressed *per meter. *Hence, the dimension of p/ħ is m^{–1}, which again yields a pure number when multiplied with the dimension of the coordinates x, y or z.

In the mentioned post, I also gave an unambiguous answer to the question as to what energy concept should be used in the equation: it is the *total *energy of the particle we are trying to describe, so that includes its kinetic energy, its rest mass energy and, finally, its potential energy in whatever force field it may find itself, such as a gravitational and/or electromagnetic force field. Now, while we know that, when talking potential energy, we have some liberty in choosing the zero point of our energy scale, this issue is easily overcome by noting that we are always talking about the amplitude to go from one state to another, or to go from one point in spacetime to another. Hence, what matters is the potential *difference*, really.

Feynman, in his description of the conservation of energy in a quantum-mechanical context, distinguishes:

- The rest energy m
_{0}∙*c*^{2}, which he describes as the rest energy ‘of the parts of the particle’. [One should remember he wrote this*before*the existence of quarks and the associated theory of matter was confirmed.] - The energy ‘over and above’ the rest energy, which includes both the kinetic energy, i.e. m∙
*v*^{2}/2 = p^{2}/(2m), as well as the ‘binding and/or excitation energy’, which he refers to as ‘internal energy’. - Finally, there is the potential energy, which we’ll denote by U.

In my previous post, I also gave you the relativistically correct formula for the energy of a particle with momentum p:

However, we will follow Feynman in his description, who uses the *non*-relativistic formula E_{p }= E_{int }+ p^{2}/(2m) + U. This is quite OK if we assume that the classical velocity of our particle does not approach the speed of light, so that covers a rather large number of real-life situations. Also, to make the example more real, we will assume the potential energy is electrostatic, and given by the formula U = q·Φ, with Φ the electrostatic potential (so just think of a number expressed in *volt*). Of course, q·Φ will be negative if the signs of q (i.e. the electric charge of our particle) and Φ are opposite, and positive if both have the same sign, as opposites *attract *and like *repel *when it comes to electric charge.

The illustration below visualizes the situation for Φ_{2 }< Φ_{1}. For example, we may assume Φ_{1} is zero, that Φ_{2 }is negative, and that our particle is positively charged, so U_{2 }= qΦ_{2 }< 0. So it’s all rather simple really: we have two areas with a potential equal to U_{1 }= qΦ_{1 }and U_{2 }= qΦ_{2 }< 0 respectively. Hence, we need to use E_{1 }= E_{int }+ p_{1}^{2}/(2m) + U_{1} to substitute ω_{1 }for E_{1}/ħ in the first area, and then E_{2 }= E_{int }+ p_{2}^{2}/(2m) + U_{2 }to substitute ω_{2 }for E_{2}/ħ in the second area, which U_{2 }– U_{1 }< 0.

The corresponding amplitudes, or wavefunctions, are:

- Ψ
_{1}(θ_{1}) = Ψ_{1}(x, t) =*a*·*e*^{−i}^{θ}^{1}*= a*·*e*^{−i[(Eint + p12/(2m) + U1)·t − p1∙x]/ħ } - Ψ
_{2}(θ_{2}) = Ψ_{2}(x, t) =*a·e*^{−i}^{θ}^{2}*= a*·*e*^{−i[(Eint + p22/(2m) + U2)·t − p2∙x]/ħ }

Now how should we *think *about these two equations? We are definitely talking *different *wavefunctions. However, having said that, there is no reason to assume the different potentials would have an impact on the *temporal *frequency. Therefore, we can boldly equate ω_{1 }and ω_{2 }and, therefore, write that:

E_{int }+ p_{1}^{2}/(2m) + U_{1 }= E_{int }+ p_{2}^{2}/(2m) + U_{2 }⇔ p_{1}^{2}/(2m) − p_{2}^{2}/(2m) = U_{2 }– U_{1 }< 0

⇒ p_{1}^{2 }− p_{2}^{2 }< 0 ⇔ p_{2}^{ }> p_{1}

What this says is that the kinetic energy, and/or the momentum, of our particle is greater in the second area, which is what we would classically expect, as a positive charged particle will pick up speed – and, therefore, momentum and kinetic energy – as it moves from an area with zero potential to an area with negative potential. However, the λ = h/p relation then implies that λ_{2}^{ }= h/p_{2}^{ }is *smaller *than λ_{1}^{ }= h/p_{2}, which is what is illustrated by the dashed lines in the illustration above – which represent surfaces of equal phase, or wavefronts – and also by the second diagram in the illustration, which shows the real part of the complex-valued amplitude and compares the wavelengths λ_{1}^{ }and λ_{2}. [As you know, the imaginary part is just like the real part but with a phase shift equal to π/2. Ideally, we should show both, but you get the idea.]

To sum it all up, **the classical statement energy conservation principle is equivalent to the quantum-mechanical statement that the temporal frequency f or ω**, i.e. the

*time-rate of change*of the phase of the wavefunction,

**does**– as long as the conditions do not change with time, of course –

*not*change**but that the**.

*spatial*frequency, i.e. the wave number k or the wavelength λ – changes as the potential energy and/or kinetic energy change**Tunneling**

The p_{1}^{2}/(2m) − p_{2}^{2}/(2m) = U_{2 }– U_{1 }equation may be re-written to illustrate the quantum-mechanical effect of tunneling, i.e. the *penetration of a potential barrier*. Indeed, we can re-write p_{1}^{2}/(2m) − p_{2}^{2}/(2m) = U_{2 }– U_{1 }as

p_{2}^{2} = 2m·[p_{1}^{2}/(2m) − (U_{2 }– U_{1})]

and, importantly, try to analyze what happens if U_{2 }– U_{1} is *larger* than p_{1}^{2}/(2m), so we get a *negative *value for p_{2}^{2}. Just imagine that Φ_{1} is zero again, and that our particle is positively charged, but that Φ_{2 }is also positive (instead of negative, as in the example above), so our particle is being repelled. In practical terms, it means that our particle just doesn’t have enough energy to “climb the potential hill”. Quantum-mechanically, however, the amplitude is still given by that equation above, and we have a purely imaginary number for p_{2}, as the square root of a *negative *number is a purely imaginary number, just like √−4 = 2*i*. So let’s denote p_{2} as *i*·p’ and let’s analyze what happens by breaking our *a·e ^{i}*

^{θ}

^{2}

*function up in two separate parts by writing:*

*a·e*

^{−i}^{θ}

^{2}

*=*

*a·e*

^{−i}^{[}

^{(E2/ħ)∙t − (i·p’/ħ)x]}=

*a·e*

^{−i}^{(E2/ħ)∙t}

*·e*

^{i2·p’·x/ħ}=

*a·e*

^{−i}^{(E2/ħ)∙t}

*·e*

^{−p’·x/ħ}.

Now, the *e*^{−p’·x/ħ} factor in our formula for *a·e ^{−i}*

^{θ}

^{2}

*is a*

*real-valued*exponential function, and it’s a decreasing function, with the same shape as the general

*e*

^{−x}function, which I depict below.

This *e*^{−p’·x/ħ} basically ‘kills’ our wavefunction as we move in the positive x-direction, past the potential barrier, which is what is illustrated below.

However, the story doesn’t finish here. We may imagine that the region with the prohibitive potential is rather small—like a few wavelengths only—and that, past that region, we’ve got another region where p_{2}^{2} = 2m·[p_{1}^{2}/(2m) − (U_{2 }– U_{1})] is *not *negative. That’s the situation that’s depicted below, which also shows what *might *happen: the amplitude decays exponentially, but does *not *reach zero and, hence, there is a possibility that a particle might make it through the barrier, and that it will be found on the other side, with a real-valued and positive momentum and, hence, with a regular wavefunction.

Feynman gives a very interesting example of this: *alpha*-decay. Alpha decay is a type of radioactive decay in which an atomic nucleus emits an α-particle (so that’s a helium nucleus, really), thereby transforming or ‘decaying’ into an atom with a reduced mass and atomic number. The Wikipedia article on it hais not bad, but Feynman’s explanation is more to the point, especially when you’ve understood all of the above. The graph below illustrates the basic idea as it shows the potential energy U of an α-particle as a function of the distance from the center. As Feynman puts it: “If one tried to shoot an α-particle with the energy E *into* the nucleus, it would feel an electrostatic repulsion from the nucleus and would, classically, get no closer than the distance r_{1}, where its total energy is equal to U. Closer in, however, the potential energy is much lower because of the strong attraction of the short-range nuclear forces. How is it then that in radioactive decay we find α-particles which started out inside the nucleus coming out with the energy E? Because they start out with the energy E inside the nucleus and “leak” through the potential barrier.”

As for the numbers involved, the mean life of an α-particle in the uranium nucleus is as long as 4.5 billion years, according to Feynman, whereas the oscillations inside the nucleus are in the range of 10^{22} cycles per second! So how can one get a number like 10^{9} years from 10^{−22} seconds? The answer, as Feynman notes, is that that exponential gives a factor of about *e*^{−45}. So that gives the very small but definite probability of leakage. Once the α-particle is in the nucleus, there is almost no amplitude at all for finding it outside. However, if you take many nuclei and wait long enough, you’ll find one. 🙂

Now, that should be it for today, but let me end this post with something I should have told you a while ago, but then I didn’t, because I thought it would distract you from the essentials. If you’ve read my previous post carefully, you’ll note that I wrote the wavefunction as Ψ(θ) = *a·e*^{i}^{θ}, rather as* a·e ^{−i}*

^{θ}, with the minus sign in front of the complex exponent. So why is that?

There is a long and a short answer to that. I’ll give the short answer. You’ll remember that the phase of our wavefunction is like the hand of a stopwatch. Now we could imagine a stopwatch going counter-clockwise, and we could actually make one. Now, there is no arbitrariness here: it’s one way or the other, depending on our other conventions, and the phase of our complex-valued wavefunction does actually turn *clockwise* if we write things the way we’re writing them, rather than anti-clockwise. That’s a direction that’s actually *not* as per the usual mathematical convention: an angle in the unit circle is usually measured *counter*-clockwise. If you’d want it that way, we can fix easily by reversing the signs inside of the bracket, so we could write θ = **k**·**x** − ω·t, which is actually what you’ll often see. But so there’s only way to get it right: there’s a direction to it, and if we use the θ = ω·t − **k**·**x**, then we need the minus sign in the Ψ(θ) = *a·e ^{−i}*

^{θ}equation.

It’s just one of those things that is easy to state, but actually gives us a lot of food for thought. Hence, I’ll probably come back to this one day. As for now, however, I think you’ve had enough. Or I’ve had enough, at least. 🙂 I hope this was not too difficult, and that you enjoyed it.

## 6 thoughts on “Potential energy and amplitudes: energy conservation and tunneling effects”