Physics for Information Science. Module 2, matter waves to quantum tunnelling.

GAPHT121 · 15 pages · Your place is saved in this browser

Contents · Where each topic sits
  1. What you will be able to do
  2. How to use these notes
  3. Where each topic sits
  4. The idea
  5. Deriving the de Broglie wavelength (route 1: from the photon's energy)
  6. Deriving it (route 2: from the photon's momentum)
  7. The de Broglie wavelength of a particle
  8. An electron accelerated through a voltage
  9. Some properties of matter waves
  10. Experimental evidence
  11. de Broglie wavelength: worked examples
  12. The idea
  13. Physical significance (Born's interpretation)
  14. A one-dimensional picture
  15. Normalisation
  16. Why a wave function must be normalisable
  17. Conditions for a well-behaved wave function
  18. Worked example
  19. Setting it up
  20. Expressing $\omega$ and $k$ through energy and momentum
  21. A second way to reach the same form
  22. Its probability density is constant
  23. A single wave cannot represent a localised particle
  24. A wave packet
  25. Phase velocity and group velocity
  26. From wave packets to uncertainty
  27. Statement
  28. Other conjugate pairs
  29. What the principle does and does not say
  30. Why we never notice it in everyday life
  31. 1. Electrons cannot exist inside the nucleus
  32. 2. Natural width of spectral lines (natural line broadening)
  33. The idea
  34. Derivation, route 1 (differentiating twice)
  35. Derivation, route 2 (energy and momentum operators)
  36. Three dimensions
  37. What $V$ does
  38. The idea
  39. Derivation
  40. Three dimensions
  41. Stationary states
  42. The set-up
  43. The equation inside the box
  44. Method A: sine and cosine form
  45. Method B: exponential form
  46. Why $n=0$ is not allowed
  47. Allowed energy levels
  48. Normalised wave functions
  49. What the solution tells us
  50. Worked examples
  51. The idea
  52. The wave function in each region
  53. Transmission coefficient
  54. Applications
  55. Numerical problems
  56. Questions asked in KTU exams
  57. Formula sheet
Text size18

Page 1 of 15

Module 2 at a glance

Quantum mechanics is the branch of mechanics that gives the mathematical description of the motion and interaction of very small particles such as electrons, protons and atoms. Classical (Newtonian) mechanics works perfectly for cricket balls and planets, but it fails at the atomic scale. It cannot explain why atoms are stable, why their spectra have sharp lines, or why electrons are diffracted by crystals. Quantum mechanics takes its place there. It is based on three ideas that have no place in everyday physics:

  1. Quantisation of energy. A confined particle cannot have any energy it likes. Only certain discrete energies are allowed.
  2. Wave–particle duality. Every moving particle has a wave associated with it, called the matter wave. The wave is described by a wave function Ψ\Psi, and ∣Ψ∣2|\Psi|^2 tells where the particle is likely to be, not where it definitely is.
  3. The uncertainty principle. Position and momentum (and also energy and time) cannot both be known exactly at the same time.

From these ideas we build the Schrödinger equation, the equation that the wave function follows. We solve it for two standard problems: a particle trapped in a box (which explains why energy comes in discrete levels and why a confined particle can never be at rest) and a particle meeting a barrier that it does not have enough energy to climb (which explains tunnelling).

What you will be able to do

  • Derive the de Broglie relation λ=h/p\lambda=h/p and calculate the wavelength of any particle, and of an electron accelerated through a voltage.
  • Explain the physical meaning of the wave function (Born's interpretation), write the conditions for a well-behaved wave function, and normalise a wave function.
  • Write the wave function of a free particle and show that ω=E/ℏ\omega=E/\hbar and k=p/ℏk=p/\hbar.
  • Explain wave packets, phase and group velocity, and how a wave packet leads to the uncertainty principle.
  • State the uncertainty principle for the three conjugate pairs and use it to show that electrons cannot exist inside the nucleus and that spectral lines have a natural width.
  • Derive the time-dependent Schrödinger equation by two methods, and derive the time-independent equation from it.
  • Solve the Schrödinger equation for a particle in a one-dimensional box: find the energy eigenvalues and the normalised eigenfunctions, and interpret them.
  • Explain quantum mechanical tunnelling, write the wave function in each region, define the transmission coefficient and list its applications.

How to use these notes

  • Each page follows idea → explanation or derivation → worked example → common mistakes.
  • Green boxes show key ideas and definitions. Amber boxes marked Watch out show mistakes that make you lose marks. Boxes marked Note on the source show where the lecture notes were corrected.
  • The derivations of the Schrödinger equation and the particle in a box are the most frequently asked essay questions. Read them, then close the notes and write them yourself, step by step, until you can do it without looking.
  • The last page has practice problems with full answers, the past KTU exam questions and a formula sheet.

Where each topic sits

PartPages
Matter wavesMatter waves, de Broglie wavelength examples
The wave functionWave function and probability, normalisation, the free-particle wave function, wave packets
UncertaintyUncertainty principle, its applications
Schrödinger equationTime-dependent form (two derivations), time-independent form
ApplicationsParticle in a box (solution, then energies and wave functions), tunnelling
RevisionPractice questions, exam questions and formula sheet

Before you start. You should know wavelength λ\lambda, frequency ν\nu, period T=1/νT=1/\nu, angular frequency ω=2πν=2π/T\omega=2\pi\nu=2\pi/T and wave number k=2π/λk=2\pi/\lambda; the photon relations E=hνE=h\nu and p=h/λp=h/\lambda; complex numbers, with eiθ=cos⁡θ+isin⁡θe^{i\theta}=\cos\theta+i\sin\theta, the conjugate z∗z^* and ∣z∣2=z∗z|z|^2=z^*z; partial derivatives and simple integration, including sin⁡2θ=12(1−cos⁡2θ)\sin^2\theta=\tfrac12(1-\cos2\theta).

Page 2 of 15

Matter waves

The idea

By the early twentieth century it was clear that light has a dual nature. Interference, diffraction and polarisation show that light is a wave. The photoelectric effect and the Compton effect show that light also behaves as a stream of particles, called photons, each carrying energy E=hνE=h\nu and momentum p=h/λp=h/\lambda.

In 1924 the French physicist Louis de Broglie argued by analogy and symmetry: if waves (light) can behave like particles, then particles should also be able to behave like waves. He proposed that every moving object has a wave associated with it under suitable conditions. These waves are called matter waves or de Broglie waves.

Matter waves. The waves associated with a moving particle are called matter waves or de Broglie waves. Their wavelength is called the de Broglie wavelength.

Deriving the de Broglie wavelength (route 1: from the photon's energy)

For a photon the energy can be written in two ways, one from its wave nature and one from its particle nature:

E=hν=hcλ(wave nature)E=h\nu=\frac{hc}{\lambda}\qquad\text{(wave nature)} E=mc2(particle nature)E=mc^2\qquad\text{(particle nature)}

Since both describe the same photon, they must be equal. This shows its wave–particle dual nature:

mc2=hcλmc^2=\frac{hc}{\lambda}

Cancel one factor of cc:

mc=hλmc=\frac{h}{\lambda}

Here mcmc (mass times velocity) is the momentum pp of the photon. So

λ=hp\lambda=\frac{h}{p}

Deriving it (route 2: from the photon's momentum)

The momentum of a photon is known directly from experiment and theory to be p=h/λp=h/\lambda, which can be written as λ=h/p\lambda=h/p. De Broglie's step was to assume that this equation is completely general. It applies not only to photons but to electrons, protons, atoms and every other particle.

The de Broglie wavelength of a particle

A particle of mass mm moving with speed vv has momentum p=mvp=mv. Its kinetic energy is KE=12mv2KE=\tfrac12mv^2, which can be written as

KE=12mv2=(mv)22m=p22m⟹p=2m KEKE=\frac12mv^2=\frac{(mv)^2}{2m}=\frac{p^2}{2m}\qquad\Longrightarrow\qquad p=\sqrt{2m\,KE}

so the de Broglie wavelength can be written in three equivalent ways:

de Broglie wavelength.

λ=hp=hmv=h2m KE\lambda=\frac{h}{p}=\frac{h}{mv}=\frac{h}{\sqrt{2m\,KE}}

An electron accelerated through a voltage

An electron of charge ee and mass mm, starting from rest and accelerated through a potential difference of VV volts, gains kinetic energy KE=eVKE=eV. So

12mv2=eV⟹v=2eVm\frac12mv^2=eV\qquad\Longrightarrow\qquad v=\sqrt{\frac{2eV}{m}}

Then the wavelength of the electron is

λ=hmv=hmm2eV=h2meV\lambda=\frac{h}{mv}=\frac{h}{m}\sqrt{\frac{m}{2eV}}=\frac{h}{\sqrt{2meV}}

This is the same as putting KE=eVKE=eV in λ=h/2m KE\lambda=h/\sqrt{2m\,KE}. Now put h=6.626×10−34 J sh=6.626\times10^{-34}\ \text{J s}, m=9.11×10−31 kgm=9.11\times10^{-31}\ \text{kg} and e=1.6×10−19 Ce=1.6\times10^{-19}\ \text{C}:

λ=6.626×10−342×9.11×10−31×1.6×10−19×V=1.227×10−9V m\lambda=\frac{6.626\times10^{-34}}{\sqrt{2\times9.11\times10^{-31}\times1.6\times10^{-19}\times V}}=\frac{1.227\times10^{-9}}{\sqrt V}\ \text{m}

Since 1 A˚=10−10 m1\ \text{Å}=10^{-10}\ \text{m}, the same result in ångström is

λ=12.27V A˚≈150V A˚\lambda=\frac{12.27}{\sqrt V}\ \text{Å}\approx\sqrt{\frac{150}{V}}\ \text{Å}

(because 12.272≈15012.27^2\approx150). Both forms appear in textbooks; they are the same formula.

Watch out. These shortcut formulas are only for electrons, because the mass and charge of the electron are already inside the number 1.227×10−91.227\times10^{-9}. For a proton, neutron, alpha particle or anything else, go back to λ=h/2m KE\lambda=h/\sqrt{2m\,KE} with the correct mass.

Some properties of matter waves

  • The lighter the particle, the longer its wavelength for a given speed (λ∝1/m\lambda\propto1/m). This is why wave effects are noticeable for electrons but not for everyday objects.
  • The faster the particle, the shorter its wavelength (λ∝1/v\lambda\propto1/v).
  • A particle at rest (v=0v=0) would have an infinite wavelength. Matter waves are associated only with moving particles.
  • Matter waves are not electromagnetic waves. They are produced by any moving particle, charged or neutral (neutrons are also diffracted).
  • The wave is not a physical disturbance of a medium. As the wave function pages explain, it is a probability wave.

Experimental evidence

De Broglie had no experimental evidence when he proposed the idea. But he could show that it explained the stationary orbits and the quantised energy levels that Bohr had assumed for the hydrogen atom in 1913. An electron orbit is allowed only if a whole number of de Broglie wavelengths fits around it, 2πr=nλ2\pi r=n\lambda, which is exactly Bohr's condition on angular momentum.

The direct proof came in 1927. Davisson and Germer in the USA, and independently G. P. Thomson in Britain, showed that electrons are diffracted by crystals exactly as X-rays are, with the wavelength given by λ=h/p\lambda=h/p. Today, electron diffraction and the electron microscope depend on this.

Key idea. Wave and particle are two aspects of the same object. Which aspect we see depends on the experiment. Diffraction shows the wave, and a collision with a detector shows the particle.

Page 3 of 15

de Broglie wavelength: worked examples

Constants: h=6.626×10−34 J sh=6.626\times10^{-34}\ \text{J s}, me=9.11×10−31 kgm_e=9.11\times10^{-31}\ \text{kg}, mp=1.673×10−27 kgm_p=1.673\times10^{-27}\ \text{kg}, mn=1.675×10−27 kgm_n=1.675\times10^{-27}\ \text{kg}, e=1.6×10−19 Ce=1.6\times10^{-19}\ \text{C}.

Example 1. Why we never notice matter waves. Find the de Broglie wavelength of (a) a ball of mass 500 g500\ \text{g} moving at 50 km/h50\ \text{km/h}, and (b) an electron accelerated through 100 V100\ \text{V}.

(a) First convert the speed: v=50 km/h=50×10003600=13.9 m/sv=50\ \text{km/h}=\dfrac{50\times1000}{3600}=13.9\ \text{m/s}. Then

λ=hmv=6.626×10−340.5×13.9=9.5×10−35 m\lambda=\frac{h}{mv}=\frac{6.626\times10^{-34}}{0.5\times13.9}=\mathbf{9.5\times10^{-35}\ \text{m}}

This is roughly 10−34 m10^{-34}\ \text{m}, which is extremely small compared with the size of the ball (and far smaller than even a nucleus). No experiment can ever detect a wave aspect in the behaviour of the ball, so for everyday objects classical mechanics is perfectly good.

(b) Using the electron formula,

λ=1.227×10−9100=1.23×10−10 m=1.23 A˚\lambda=\frac{1.227\times10^{-9}}{\sqrt{100}}=\mathbf{1.23\times10^{-10}\ \text{m}}=1.23\ \text{Å}

This is about the size of an atom, and it is much larger than the classical radius of the electron (≈10−15 m\approx10^{-15}\ \text{m}). So the wave nature of the electron cannot be neglected. It is also comparable to the spacing between atoms in a crystal, which is why crystals diffract electrons.

Note on the source. One lecture set gave the wavelength of the electron in part (b) as 0.123×10−10 m0.123\times10^{-10}\ \text{m}. The correct value is 1.23×10−10 m1.23\times10^{-10}\ \text{m}, ten times larger, as the formula λ=1.227×10−9/V\lambda=1.227\times10^{-9}/\sqrt V shows for V=100V=100.

Example 2. The Davisson–Germer experiment. In the original experiment the electrons were accelerated through 54 V54\ \text{V}. Find their wavelength.

λ=1.227×10−954=1.227×10−97.35=1.67×10−10 m=1.67 A˚\lambda=\frac{1.227\times10^{-9}}{\sqrt{54}}=\frac{1.227\times10^{-9}}{7.35}=\mathbf{1.67\times10^{-10}\ \text{m}}=1.67\ \text{Å}

This matched the wavelength measured from the diffraction pattern of a nickel crystal, and this is how the experiment confirmed de Broglie's hypothesis.

Example 3. A slow (thermal) neutron. Find the de Broglie wavelength of a neutron with kinetic energy 0.025 eV0.025\ \text{eV}.

The electron formula cannot be used. Change the energy to joules, KE=0.025×1.6×10−19=4.0×10−21 JKE=0.025\times1.6\times10^{-19}=4.0\times10^{-21}\ \text{J}, and use the general form:

λ=h2mnKE=6.626×10−342×1.675×10−27×4.0×10−21=1.81×10−10 m\lambda=\frac{h}{\sqrt{2m_nKE}}=\frac{6.626\times10^{-34}}{\sqrt{2\times1.675\times10^{-27}\times4.0\times10^{-21}}}=\mathbf{1.81\times10^{-10}\ \text{m}}

Again about 1.8 A˚1.8\ \text{Å}. Slow neutrons are diffracted by crystals and are used to study crystal structures.

Example 4. A proton versus an electron. A proton and an electron are each accelerated through 100 V100\ \text{V}. Compare their wavelengths.

Both gain the same kinetic energy eVeV, so λ=h/2meV∝1/m\lambda=h/\sqrt{2meV}\propto1/\sqrt m:

λp=6.626×10−342×1.673×10−27×1.6×10−19×100=2.86×10−12 m\lambda_p=\frac{6.626\times10^{-34}}{\sqrt{2\times1.673\times10^{-27}\times1.6\times10^{-19}\times100}}=\mathbf{2.86\times10^{-12}\ \text{m}} λeλp=mpme=1836≈43\frac{\lambda_e}{\lambda_p}=\sqrt{\frac{m_p}{m_e}}=\sqrt{1836}\approx43

The proton's wavelength is about 43 times shorter, because it is 1836 times heavier.

Example 5. Reverse problem. Through what voltage must an electron be accelerated to have a wavelength of 1 A˚1\ \text{Å}?

From λ=150/V A˚\lambda=\sqrt{150/V}\ \text{Å}: 1=150/V1=\sqrt{150/V}, so V=150 VV=\mathbf{150\ \text{V}}.

Watch out. Always work in SI units: mass in kg (500 g is 0.50.5 kg), speed in m/s (divide km/h by 3.6) and energy in joules (multiply eV by 1.6×10−191.6\times10^{-19}). Most wrong answers here come from a unit that was not converted.

Page 4 of 15

The wave function and probability

The idea

In every wave, something varies periodically with space and time:

WaveWhat varies
Water wavethe height of the water surface
Sound wavethe pressure of the air
Light (electromagnetic) wavethe electric and magnetic fields
Matter wavethe wave function Ψ\Psi

The quantity whose variations make up matter waves is called the wave function, written Ψ\Psi (the Greek capital letter psi, pronounced "sigh"). In general it is a function of all three space coordinates and of time:

Ψ=Ψ(x,y,z,t)\Psi=\Psi(x,y,z,t)

The wave function is a mathematical function that represents the state of the particle. The mathematical theory of quantum mechanics tells how to use Ψ(x,y,z,t)\Psi(x,y,z,t) to find the position, velocity, momentum, energy, angular momentum and so on of the particle. Everything that can be known about the particle is contained in Ψ\Psi.

Physical significance (Born's interpretation)

  1. Ψ\Psi itself cannot be measured. The wave function is in general a complex quantity (it contains i=−1i=\sqrt{-1}). A complex number cannot be the result of a measurement, so Ψ\Psi has no direct physical meaning by itself.
  2. ∣Ψ∣2|\Psi|^2 is real and measurable. The square of the absolute value of Ψ\Psi,
∣Ψ∣2=Ψ∗Ψ|\Psi|^2=\Psi^*\Psi

where Ψ∗\Psi^* is the complex conjugate of Ψ\Psi, is always real and positive.

  1. ∣Ψ∣2|\Psi|^2 is a probability density. ∣Ψ(x,y,z,t)∣2|\Psi(x,y,z,t)|^2 represents the probability of finding the particle at the point (x,y,z)(x,y,z) at time tt, per unit volume. It is called the probability density or the probability distribution function. Where ∣Ψ∣2|\Psi|^2 is large, the particle is likely to be found. Where it is small, the particle is unlikely to be found.
  2. Probability in a region. The probability of finding the particle in a small volume dVdV around (x,y,z)(x,y,z) at time tt is ∣Ψ∣2 dV|\Psi|^2\,dV. For a finite region it is the integral
P=∫∣Ψ∣2 dVP=\int|\Psi|^2\,dV

PP always lies between 00 (never there) and 11 (certainly there). For example, if P=0.4P=0.4 for a region, there is a 40 % chance of finding the particle in that region.

  1. Probability amplitude. Ψ\Psi itself is not a probability, but its square is. So Ψ\Psi is also called the probability amplitude.
  2. Statistical link between wave and particle. The wave function connects the particle nature and the wave nature statistically. It does not tell exactly where the particle is. It tells where the particle is likely to be. If we repeat the same experiment many times, the particle is found most often where ∣Ψ∣2|\Psi|^2 is largest.

Key idea. Quantum mechanics does not predict the exact location of a particle. It predicts the probability of finding it in each region of space. The matter wave is a probability wave.

A one-dimensional picture

For a particle moving only along xx, the wave function is Ψ(x,t)\Psi(x,t) and ∣Ψ(x,t)∣2 dx|\Psi(x,t)|^2\,dx is the probability of finding the particle between xx and x+dxx+dx. The probability of finding it between x=ax=a and x=bx=b is

P(a≤x≤b)=∫ab∣Ψ∣2 dxP(a\le x\le b)=\int_a^b|\Psi|^2\,dx

Watch out. The probability density is ∣Ψ∣2=Ψ∗Ψ|\Psi|^2=\Psi^*\Psi, not Ψ2\Psi^2. For a complex Ψ\Psi, Ψ2\Psi^2 is itself complex and cannot be a probability. Only when Ψ\Psi is real (as for the particle in a box) are the two the same.

Page 5 of 15

Normalisation and well-behaved wave functions

Normalisation

A particle that exists in a particular region must be found somewhere in that region. So the probability of finding it anywhere in the region must be exactly 1 (100 %):

∫V∣Ψ∣2 dV=1\int_V|\Psi|^2\,dV=1

A wave function that satisfies this is said to be normalised (normalised to unity).

Now suppose we have a wave function Ψ\Psi that solves the Schrödinger equation but is not normalised. If we multiply a solution by a constant AA (not depending on xx, yy, zz), we get another solution, ΨN=AΨ\Psi_N=A\Psi. We choose AA so that the new wave function is normalised. Adding the probability at all points in the volume VV,

∑all points in V∣ΨN∣2=1or∫V∣ΨN∣2 dV=1\sum_{\text{all points in }V}|\Psi_N|^2=1\qquad\text{or}\qquad\int_V|\Psi_N|^2\,dV=1 that is∣A∣2∫V∣Ψ∣2 dV=1\text{that is}\qquad |A|^2\int_V|\Psi|^2\,dV=1

Normalisation. The process of multiplying a wave function by a suitable constant so that the probability of finding the particle within the region becomes exactly 1 is called normalisation. The constant is called the normalisation constant.

Why a wave function must be normalisable

∣Ψ∣2|\Psi|^2 is supposed to be a probability density. A probability can only be interpreted if the total adds up to 1, so the integral ∫∣Ψ∣2 dV\int|\Psi|^2\,dV must be finite. A wave function whose integral is infinite cannot be normalised and cannot describe a real particle.

Conditions for a well-behaved wave function

The physically acceptable ("well-behaved") wave functions satisfy these conditions:

  1. Ψ\Psi is finite everywhere. An infinite value would mean an infinite probability.
  2. Ψ\Psi is single-valued. At any point and time there must be only one value of the probability.
  3. Ψ\Psi is continuous. The probability cannot jump suddenly from place to place.
  4. The first derivatives ∂Ψ∂x\dfrac{\partial\Psi}{\partial x}, ∂Ψ∂y\dfrac{\partial\Psi}{\partial y}, ∂Ψ∂z\dfrac{\partial\Psi}{\partial z} are continuous and single-valued, wherever the potential is finite. These derivatives decide the momentum of the particle, which must change smoothly.
  5. Ψ\Psi is normalisable, so Ψ→0\Psi\to0 as x,y,z→±∞x,y,z\to\pm\infty.

Worked example

Example 1. A particle confined between x=0x=0 and x=Lx=L has the wave function ψ=A x(L−x)\psi=A\,x(L-x) inside, and ψ=0\psi=0 outside. Find the normalisation constant AA.

Normalisation requires

∫0L∣A∣2x2(L−x)2 dx=1\int_0^L|A|^2x^2(L-x)^2\,dx=1

Expand x2(L−x)2=L2x2−2Lx3+x4x^2(L-x)^2=L^2x^2-2Lx^3+x^4 and integrate term by term:

∫0L(L2x2−2Lx3+x4)dx=L53−2L54+L55=L5(10−15+630)=L530\int_0^L\left(L^2x^2-2Lx^3+x^4\right)dx=\frac{L^5}{3}-\frac{2L^5}{4}+\frac{L^5}{5}=L^5\left(\frac{10-15+6}{30}\right)=\frac{L^5}{30}

So ∣A∣2L530=1|A|^2\dfrac{L^5}{30}=1, giving

A=30L5A=\mathbf{\sqrt{\frac{30}{L^5}}}

This function is finite, single-valued, continuous and zero at both walls, so it is well behaved.

Watch out. When normalising, square the absolute value of the constant (∣A∣2|A|^2), integrate over the whole region where the particle can be, and remember to take the square root at the end.

Page 6 of 15

The wave function of a free particle

Setting it up

Consider a particle moving freely along the xx direction with energy EE and momentum pp. Its matter wave is a travelling wave, which we write in complex form as

Ψ=a e−i(ωt−kx)(1)\Psi=a\,e^{-i(\omega t-kx)}\qquad(1)

where aa is the amplitude, ω\omega is the angular frequency and kk is the wave number:

ω=2πT=2πνk=2πλ\omega=\frac{2\pi}{T}=2\pi\nu\qquad\qquad k=\frac{2\pi}{\lambda}

Expressing ω\omega and kk through energy and momentum

Angular frequency. Use Planck's relation E=hνE=h\nu, so ν=E/h\nu=E/h:

ω=2πT=2πν=2πhνh=2πh E=Eℏ(2)\omega=\frac{2\pi}{T}=2\pi\nu=\frac{2\pi h\nu}{h}=\frac{2\pi}{h}\,E=\frac{E}{\hbar}\qquad(2)

where ℏ=h2π=1.055×10−34 J s\hbar=\dfrac{h}{2\pi}=1.055\times10^{-34}\ \text{J s} (read "h-bar").

Wave number. Use the de Broglie relation λ=h/p\lambda=h/p:

k=2πλ=2πph=pℏ(3)k=\frac{2\pi}{\lambda}=\frac{2\pi p}{h}=\frac{p}{\hbar}\qquad(3)

Put (2) and (3) in (1):

Ψ=a e−i(Eℏt−pℏx)\Psi=a\,e^{-i\left(\frac{E}{\hbar}t-\frac{p}{\hbar}x\right)}  Ψ=a e−iℏ(Et−px) \boxed{\,\Psi=a\,e^{-\frac{i}{\hbar}(Et-px)}\,}

This is the wave function of a free particle of energy EE and momentum pp. It is the starting point for the Schrödinger equation.

A second way to reach the same form

The same result can be reached through frequency and wavelength (this is the method used in the second lecture set). Write the wave as

Ψ=A e−iω(t−xv)\Psi=A\,e^{-i\omega\left(t-\frac xv\right)}

where vv is the wave speed. With ω=2πν\omega=2\pi\nu and v=νλv=\nu\lambda:

Ψ=A e−i2πν(t−xνλ)=A e−2πi(νt−xλ)\Psi=A\,e^{-i2\pi\nu\left(t-\frac{x}{\nu\lambda}\right)}=A\,e^{-2\pi i\left(\nu t-\frac{x}{\lambda}\right)}

Now use E=hνE=h\nu (so ν=E/h\nu=E/h) and λ=h/p\lambda=h/p (so 1/λ=p/h1/\lambda=p/h):

Ψ=A e−2πi(Eth−pxh)=A e−2πih(Et−px)\Psi=A\,e^{-2\pi i\left(\frac{Et}{h}-\frac{px}{h}\right)}=A\,e^{-\frac{2\pi i}{h}(Et-px)}

Since 2πh=1ℏ\dfrac{2\pi}{h}=\dfrac1\hbar,

Ψ=A e−iℏ(Et−px)\Psi=A\,e^{-\frac{i}{\hbar}(Et-px)}

exactly as before.

Its probability density is constant

∣Ψ∣2=Ψ∗Ψ=A e+iℏ(Et−px)⋅A e−iℏ(Et−px)=A2e0=A2|\Psi|^2=\Psi^*\Psi=A\,e^{+\frac{i}{\hbar}(Et-px)}\cdot A\,e^{-\frac{i}{\hbar}(Et-px)}=A^2e^{0}=A^2

This is a constant. It does not depend on position or time. The particle is equally likely to be found at any point along the xx axis. This makes sense, because this wave has one exact wavelength and so describes a particle with an exactly known momentum but a completely unknown position. The next page shows how a particle that is located somewhere is described.

Watch out. The sign convention matters only for consistency. These notes always use e−iℏ(Et−px)e^{-\frac{i}{\hbar}(Et-px)}, which represents a wave moving in the +x+x direction. Use the same form in the whole derivation.

Page 7 of 15

Wave packets

A single wave cannot represent a localised particle

The wave function of a moving body shows the probability of finding it at a particular place at a particular time. Consider the free-particle wave function Ψ=Ae−i(ωt−kx)\Psi=Ae^{-i(\omega t-kx)}. Its real part, Acos⁡(ωt−kx)A\cos(\omega t-kx), is a sine-shaped wave of constant amplitude that goes on for ever in both directions.

As the previous page showed, its probability density is

∣Ψ∣2=A2|\Psi|^2=A^2

which is the same everywhere. So this wave says the particle is equally likely to be found at any point along the xx axis. That describes a particle with an exact wavelength, and hence an exact momentum, but with no position at all.

A wave packet

A real particle is more likely to be found in some places than in others at a given time. Its wave function must therefore be large in a small region of space and nearly zero everywhere else. Graphically this is a short burst of waves, a wave packet (also called a wave group).

Top: a single wavelength extends everywhere. Bottom: a wave packet confined to a region of width Δx

How a wave packet is formed. A wave packet is formed by the superposition (adding together) of a large number of waves of slightly different wavelengths and frequencies. In one small region the waves are in step and reinforce each other, which gives a large amplitude. Away from that region they go out of step and cancel, so the amplitude dies away.

  • The larger the number of component waves (the wider the range of wavelengths), the smaller the length of the wave packet.
  • The smaller the number of component waves, the longer the packet.

When the particle moves, the wave packet moves as a whole.

Phase velocity and group velocity

Two different velocities describe a wave packet:

  • The phase velocity vpv_p is the velocity of the individual component waves (the crests inside the packet).
  • The group velocity vgv_g is the velocity with which the wave packet as a whole (its envelope) moves.

The group velocity is equal to the velocity of the particle, so it is also called the particle velocity. It is the group, not the individual waves, that carries the particle along.

From wave packets to uncertainty

This picture leads directly to the uncertainty principle:

  • If the wave function is a single sine wave of wavelength λ\lambda, the momentum p=h/λp=h/\lambda is known exactly. But ∣Ψ∣2|\Psi|^2 has the same value all over space, so the position at a given time cannot be determined. The particle is "spread out" in space: exact momentum means a very large uncertainty in position.
  • If the wave function is a narrow packet of width Δx\Delta x, then ∣Ψ∣2|\Psi|^2 is zero everywhere except in that region, so the position is known to within Δx\Delta x. But the packet is formed by the interference of a large number of waves of slightly different wavelengths. Since λ\lambda can take any value in that range, the momentum p=h/λp=h/\lambda can also take a range of values. So there is an uncertainty in momentum Δp\Delta p.
  • For a narrower packet the number of component waves must be larger, so Δp\Delta p is larger.

Key idea. You cannot make Δx\Delta x and Δp\Delta p both small at the same time. This trade-off is not a fault of our instruments. It is part of the wave nature of matter. It is Heisenberg's uncertainty principle (next page).

Watch out. The individual waves inside a packet (phase velocity) and the packet itself (group velocity) move at different speeds. When an exam asks for the speed of the particle, it is the group velocity.

Page 8 of 15

Heisenberg's uncertainty principle

Statement

When we treat a moving particle as a wave group, there are basic limits to the accuracy with which we can measure particle properties such as position and momentum.

Heisenberg's uncertainty principle. It is impossible to know both the position and the momentum of a particle exactly at the same time. If Δx\Delta x is the uncertainty in the position and Δp\Delta p the uncertainty in the corresponding component of momentum (pxp_x), then

Δx Δp ≥ ℏ2\Delta x\,\Delta p\ \ge\ \frac{\hbar}{2}

where ℏ=h2π=1.055×10−34 J s\hbar=\dfrac{h}{2\pi}=1.055\times10^{-34}\ \text{J s}, so that ℏ2=h4π\dfrac\hbar2=\dfrac{h}{4\pi}.

In simple words: the product of the uncertainty in position and the uncertainty in momentum can never be smaller than ℏ/2\hbar/2. The more precisely one is known, the less precisely the other can be known.

Other conjugate pairs

A pair of variables connected by an uncertainty relation is called a conjugate pair. Besides position and momentum, the important pairs are energy and time, and angular momentum and angular position:

ΔE Δt ≥ ℏ2\Delta E\,\Delta t\ \ge\ \frac{\hbar}{2} ΔL Δθ ≥ ℏ2\Delta L\,\Delta\theta\ \ge\ \frac{\hbar}{2}

Here ΔE\Delta E and Δt\Delta t are the uncertainties in finding the energy and time, and ΔL\Delta L and Δθ\Delta\theta are the uncertainties in the angular momentum and the angle.

There is no uncertainty relation between Δx\Delta x and a component of momentum other than the xx component. For example, xx and pyp_y can both be known exactly at the same time.

Note on the source. One lecture set wrote these relations with ℏ\hbar on the right (Δx Δp≥ℏ\Delta x\,\Delta p\ge\hbar, ΔE Δt≥ℏ\Delta E\,\Delta t\ge\hbar, ΔJ Δθ≥ℏ\Delta J\,\Delta\theta\ge\hbar) and the other with ℏ/2\hbar/2. The exact lower limit is ℏ/2\hbar/2, used in all these notes. Some textbooks write ℏ\hbar or even hh because they make an estimate of the order of magnitude. Then the answers change by a factor of 2 or 4π4\pi. Write which form you use.

What the principle does and does not say

  • It is not about poor instruments or experimental error. Even with perfect instruments, Δx Δp\Delta x\,\Delta p cannot be smaller than ℏ/2\hbar/2. It is a property of nature, coming from the wave nature of matter.
  • It applies to simultaneous knowledge. You can measure position very precisely, but then the momentum becomes very uncertain, and the other way round.
  • Because ℏ\hbar is tiny, the limits matter only on the atomic scale. For everyday objects the uncertainties are far too small to notice.

Why we never notice it in everyday life

Example 1. A bullet of mass 25 g25\ \text{g} moves with a speed of 400 m/s400\ \text{m/s}. The speed is measured to an accuracy of 0.02 %0.02\,\%. How precisely can the position of the bullet be found?

The uncertainty in speed is 0.02 %0.02\,\% of 400 m/s400\ \text{m/s}:

Δv=400×0.02100=0.08 m/s\Delta v=400\times\frac{0.02}{100}=0.08\ \text{m/s}

The uncertainty in momentum is

Δp=m Δv=0.025×0.08=2×10−3 kg m/s\Delta p=m\,\Delta v=0.025\times0.08=2\times10^{-3}\ \text{kg m/s}

So

Δx ≥ ℏ2 Δp=1.055×10−342×2×10−3=2.64×10−32 m\Delta x\ \ge\ \frac{\hbar}{2\,\Delta p}=\frac{1.055\times10^{-34}}{2\times2\times10^{-3}}=\mathbf{2.64\times10^{-32}\ \text{m}}

This uncertainty is far smaller than even the size of an atom (about 10−10 m10^{-10}\ \text{m}), so for a bullet it is completely negligible.

Note on the source. The lecture notes gave 2.626×10−32 m2.626\times10^{-32}\ \text{m}. When we calculate h/(4π Δp)h/(4\pi\,\Delta p) with h=6.625×10−34h=6.625\times10^{-34}, we get 2.64×10−32 m2.64\times10^{-32}\ \text{m}.

Example 2. An electron in an atom. An electron is confined within an atom of size about 10−10 m10^{-10}\ \text{m}. Find the minimum uncertainty in its speed.

Δp ≥ ℏ2 Δx=1.055×10−342×10−10=5.3×10−25 kg m/s\Delta p\ \ge\ \frac{\hbar}{2\,\Delta x}=\frac{1.055\times10^{-34}}{2\times10^{-10}}=5.3\times10^{-25}\ \text{kg m/s} Δv=Δpm=5.3×10−259.11×10−31=5.8×105 m/s\Delta v=\frac{\Delta p}{m}=\frac{5.3\times10^{-25}}{9.11\times10^{-31}}=\mathbf{5.8\times10^{5}\ \text{m/s}}

An uncertainty of more than half a million metres per second. For an electron in an atom, the uncertainty principle is very large, and the idea of a definite orbit with a definite speed breaks down.

Watch out. Change a percentage accuracy into an actual Δv\Delta v before multiplying by the mass, and use the mass in kilograms. Also remember that Δx\Delta x comes out as a minimum (≥\ge), not an exact value.

Page 9 of 15

Applications of the uncertainty principle

Planck's constant is so small that the limitations given by the uncertainty principle are important only at the subatomic level. On that scale, the principle is a powerful tool for understanding many phenomena.

1. Electrons cannot exist inside the nucleus

In the early days, people thought that the electrons coming out of a nucleus as β\beta rays in radioactivity actually existed inside the nucleus before they were emitted. The uncertainty principle shows that this is impossible.

Step 1: the uncertainty in position. The diameter of a nucleus is typically about 10−14 m10^{-14}\ \text{m}. If an electron existed inside it, the maximum uncertainty in its position would be the size of the nucleus:

Δx=10−14 m\Delta x=10^{-14}\ \text{m}

Step 2: the minimum uncertainty in momentum.

Δp=ℏ2 Δx=h4π Δx=6.625×10−344×3.14×10−14=5.3×10−21 kg m/s\Delta p=\frac{\hbar}{2\,\Delta x}=\frac{h}{4\pi\,\Delta x}=\frac{6.625\times10^{-34}}{4\times3.14\times10^{-14}}=5.3\times10^{-21}\ \text{kg m/s}

Step 3: the minimum momentum. If this is the uncertainty in the momentum, the momentum itself must be at least as large:

pmin⁡≈5.3×10−21 kg m/sp_{\min}\approx5.3\times10^{-21}\ \text{kg m/s}

Step 4: the minimum energy. At such a large momentum the electron moves at nearly the speed of light, so its rest-mass energy can be neglected and its kinetic energy is about pcpc:

KE≈pc=5.3×10−21×3×108=1.6×10−12 J=1.6×10−121.6×10−19 eV≈10 MeVKE\approx pc=5.3\times10^{-21}\times3\times10^{8}=1.6\times10^{-12}\ \text{J}=\frac{1.6\times10^{-12}}{1.6\times10^{-19}}\ \text{eV}\approx\mathbf{10\ \text{MeV}}

Step 5: compare with experiment. An electron inside the nucleus would need a kinetic energy of at least about 10 MeV10\ \text{MeV}. But the electrons emitted in β\beta decay have energies of only a fraction of this, a few MeV at most. Hence electrons do not exist inside the nucleus. The electrons emitted in β\beta decay are created at the moment of decay, when a neutron changes into a proton.

Note on the source. One lecture set used Δx Δp≥ℏ\Delta x\,\Delta p\ge\hbar with Δx=10−14 m\Delta x=10^{-14}\ \text{m} taken as the nuclear radius, and got Δp=1.1×10−20 kg m/s\Delta p=1.1\times10^{-20}\ \text{kg m/s} and about 20 MeV20\ \text{MeV}. The other used ℏ/2\hbar/2 and got 10 MeV10\ \text{MeV}. In both cases the conclusion is the same, because both values are far above the energies seen in β\beta decay.

2. Natural width of spectral lines (natural line broadening)

When excited atoms return to the ground state, they emit radiation. The emitted lines do not have a perfectly sharp frequency. Each line has a small spread of frequencies, called its natural width. The uncertainty principle explains this.

Step 1: the uncertainty in time. The average time between the excitation of an atom and the moment it radiates is called the lifetime of the excited state. It is typically about 10−8 s10^{-8}\ \text{s}. This is the uncertainty in the time of the de-excitation:

Δt=10−8 s\Delta t=10^{-8}\ \text{s}

Step 2: the uncertainty in energy. From ΔE Δt≈ℏ/2\Delta E\,\Delta t\approx\hbar/2,

ΔE≈ℏ2 Δt=6.625×10−342×3.14×2×10−8≈5.3×10−27 J\Delta E\approx\frac{\hbar}{2\,\Delta t}=\frac{6.625\times10^{-34}}{2\times3.14\times2\times10^{-8}}\approx5.3\times10^{-27}\ \text{J}

Step 3: the uncertainty in frequency. Since E=hνE=h\nu, ΔE=h Δν\Delta E=h\,\Delta\nu:

Δν=ΔEh=ℏ2h Δt=14π Δt=14π×10−8≈8×106 Hz\Delta\nu=\frac{\Delta E}{h}=\frac{\hbar}{2h\,\Delta t}=\frac{1}{4\pi\,\Delta t}=\frac{1}{4\pi\times10^{-8}}\approx\mathbf{8\times10^{6}\ \text{Hz}}

So the emitted frequency is not sharp: the spectral line has a finite width of a few megahertz. The shorter the lifetime of the excited state, the broader the line.

Note on the source. The lecture set that used ΔE Δt≥ℏ\Delta E\,\Delta t\ge\hbar got Δν=1/(2πΔt)≈16 MHz\Delta\nu=1/(2\pi\Delta t)\approx16\ \text{MHz} for the same lifetime. It differs from 8 MHz8\ \text{MHz} only by the factor of 2 in the convention.

Key idea. Both applications use the same idea: a quantity confined to a small range (Δx\Delta x or Δt\Delta t) forces a large spread in its conjugate (Δp\Delta p or ΔE\Delta E).

Watch out. Use Δν=ΔE/h\Delta\nu=\Delta E/h (with hh), not ΔE/ℏ\Delta E/\hbar. If you combine it with ΔE=ℏ/(2Δt)\Delta E=\hbar/(2\Delta t), you get the neat result Δν=1/(4πΔt)\Delta\nu=1/(4\pi\Delta t), which is worth remembering.

Page 10 of 15

The time-dependent Schrödinger equation

The idea

How can we find Ψ\Psi for a particle? The answer is the Schrödinger equation. It is to quantum mechanics what Newton's second law is to classical mechanics: given the forces acting (described by the potential energy VV), it tells how the wave function behaves in space and time.

The Schrödinger equation cannot be derived from older physics. We get it by starting from the wave function of a free particle and using Planck's relation E=hνE=h\nu and de Broglie's relation p=h/λp=h/\lambda. The equation we reach is then taken as the general law, valid for any potential. Two equivalent methods are given below. Learn one thoroughly. The second one is useful because it also gives the energy and momentum operators.

Derivation, route 1 (differentiating twice)

Consider a particle moving in the xx direction with momentum pp and energy EE. Its wave function is (see the free-particle page)

Ψ=A e−iℏ(Et−px)(1)\Psi=A\,e^{-\frac{i}{\hbar}(Et-px)}\qquad(1)

Step 1: differentiate (1) with respect to tt. Only the exponent depends on tt, and its derivative is −iEℏ-\dfrac{iE}{\hbar}:

∂Ψ∂t=−iEℏ A e−iℏ(Et−px)=−iEℏ Ψ(using (1))\frac{\partial\Psi}{\partial t}=-\frac{iE}{\hbar}\,A\,e^{-\frac{i}{\hbar}(Et-px)}=-\frac{iE}{\hbar}\,\Psi\qquad\text{(using (1))}

Multiply both sides by iℏi\hbar. Since iℏ⋅(−iEℏ)=−i2E=Ei\hbar\cdot\left(-\dfrac{iE}{\hbar}\right)=-i^2E=E:

EΨ=iℏ ∂Ψ∂t(2)E\Psi=i\hbar\,\frac{\partial\Psi}{\partial t}\qquad(2)

Step 2: differentiate (1) twice with respect to xx. Each differentiation brings down a factor ipℏ\dfrac{ip}{\hbar}:

∂Ψ∂x=ipℏ Ψ\frac{\partial\Psi}{\partial x}=\frac{ip}{\hbar}\,\Psi ∂2Ψ∂x2=ipℏ⋅ipℏ Ψ=i2p2ℏ2 Ψ=−p2ℏ2 Ψ\frac{\partial^2\Psi}{\partial x^2}=\frac{ip}{\hbar}\cdot\frac{ip}{\hbar}\,\Psi=\frac{i^2p^2}{\hbar^2}\,\Psi=-\frac{p^2}{\hbar^2}\,\Psi

Multiply both sides by −ℏ2-\hbar^2:

p2Ψ=−ℏ2 ∂2Ψ∂x2(3)p^2\Psi=-\hbar^2\,\frac{\partial^2\Psi}{\partial x^2}\qquad(3)

Step 3: total energy. The total energy of the particle is the sum of its kinetic energy and potential energy VV. With KE=12mv2=(mv)22m=p22mKE=\tfrac12mv^2=\dfrac{(mv)^2}{2m}=\dfrac{p^2}{2m}:

E=p22m+VE=\frac{p^2}{2m}+V

Multiply by Ψ\Psi:

EΨ=p2Ψ2m+VΨ(4)E\Psi=\frac{p^2\Psi}{2m}+V\Psi\qquad(4)

Step 4: put (2) and (3) in (4).

 iℏ ∂Ψ∂t=−ℏ22m ∂2Ψ∂x2+VΨ (5)\boxed{\,i\hbar\,\frac{\partial\Psi}{\partial t}=-\frac{\hbar^2}{2m}\,\frac{\partial^2\Psi}{\partial x^2}+V\Psi\,}\qquad(5)

This is the time-dependent Schrödinger equation in one dimension. It is the most general form of the Schrödinger equation. VV and Ψ\Psi may depend on both position and time, which is why it is called time-dependent.

Derivation, route 2 (energy and momentum operators)

Start from the same wave function, Ψ=Ae−iℏ(Et−px)\Psi=Ae^{-\frac{i}{\hbar}(Et-px)}.

Differentiating partially with respect to tt:

∂Ψ∂t=Ae−iℏ(Et−px)⋅(−iEℏ)=−iEℏ Ψ⟹EΨ=−ℏi∂Ψ∂t=iℏ ∂Ψ∂t(6)\frac{\partial\Psi}{\partial t}=Ae^{-\frac{i}{\hbar}(Et-px)}\cdot\left(\frac{-iE}{\hbar}\right)=\frac{-iE}{\hbar}\,\Psi \qquad\Longrightarrow\qquad E\Psi=-\frac{\hbar}{i}\frac{\partial\Psi}{\partial t}=i\hbar\,\frac{\partial\Psi}{\partial t}\qquad(6)

(using −1i=i-\dfrac1i=i). Differentiating partially with respect to xx:

∂Ψ∂x=Ae−iℏ(Et−px)⋅(ipℏ)=ipℏ Ψ⟹pΨ=ℏi∂Ψ∂x=−iℏ ∂Ψ∂x(7)\frac{\partial\Psi}{\partial x}=Ae^{-\frac{i}{\hbar}(Et-px)}\cdot\left(\frac{ip}{\hbar}\right)=\frac{ip}{\hbar}\,\Psi \qquad\Longrightarrow\qquad p\Psi=\frac{\hbar}{i}\frac{\partial\Psi}{\partial x}=-i\hbar\,\frac{\partial\Psi}{\partial x}\qquad(7)

Equations (6) and (7) say that, when they act on Ψ\Psi, energy and momentum can be replaced by operators:

Energy and momentum operators.

E → iℏ ∂∂tp → −iℏ ∂∂xE\ \to\ i\hbar\,\frac{\partial}{\partial t}\qquad\qquad p\ \to\ -i\hbar\,\frac{\partial}{\partial x}

These are used to find the energy and momentum of a particle when its wave function is known.

Now write the total energy, E=p22m+VE=\dfrac{p^2}{2m}+V, multiply by Ψ\Psi, and replace EE and pp by their operators (p2p^2 means applying the momentum operator twice):

iℏ ∂Ψ∂t=12m(−iℏ ∂∂x)2Ψ+VΨi\hbar\,\frac{\partial\Psi}{\partial t}=\frac{1}{2m}\left(-i\hbar\,\frac{\partial}{\partial x}\right)^2\Psi+V\Psi

Since (−iℏ)2=i2ℏ2=−ℏ2(-i\hbar)^2=i^2\hbar^2=-\hbar^2,

iℏ ∂Ψ∂t=−ℏ22m ∂2Ψ∂x2+VΨi\hbar\,\frac{\partial\Psi}{\partial t}=-\frac{\hbar^2}{2m}\,\frac{\partial^2\Psi}{\partial x^2}+V\Psi

which is equation (5) again.

Three dimensions

If Ψ\Psi depends on xx, yy and zz, the second derivative is replaced by the sum of all three:

iℏ ∂Ψ∂t=−ℏ22m(∂2Ψ∂x2+∂2Ψ∂y2+∂2Ψ∂z2)+VΨi\hbar\,\frac{\partial\Psi}{\partial t}=-\frac{\hbar^2}{2m}\left(\frac{\partial^2\Psi}{\partial x^2}+\frac{\partial^2\Psi}{\partial y^2}+\frac{\partial^2\Psi}{\partial z^2}\right)+V\Psi oriℏ ∂Ψ∂t=−ℏ22m∇2Ψ+VΨ\text{or}\qquad i\hbar\,\frac{\partial\Psi}{\partial t}=-\frac{\hbar^2}{2m}\nabla^2\Psi+V\Psi

where ∇2=∂2∂x2+∂2∂y2+∂2∂z2\nabla^2=\dfrac{\partial^2}{\partial x^2}+\dfrac{\partial^2}{\partial y^2}+\dfrac{\partial^2}{\partial z^2} is the Laplacian operator.

What VV does

The potential energy function VV represents the effect of the surroundings on the particle. Any restriction on the motion of the particle (walls, barriers, the attraction of a nucleus) appears in VV. Once VV is known, the Schrödinger equation can be solved for the wave function Ψ\Psi of the particle, and from Ψ\Psi we get ∣Ψ∣2|\Psi|^2 and all the other physical quantities.

Watch out. Be careful with the signs: there is a minus sign in front of ℏ22m∂2Ψ∂x2\dfrac{\hbar^2}{2m}\dfrac{\partial^2\Psi}{\partial x^2}, and iℏi\hbar (not −iℏ-i\hbar) on the left. They come from (ipℏ)2=−p2ℏ2\left(\tfrac{ip}{\hbar}\right)^2=-\tfrac{p^2}{\hbar^2} and from multiplying −iEℏΨ-\tfrac{iE}{\hbar}\Psi by iℏi\hbar. In an exam, write each step. Examiners give marks for the steps, not only for the final equation.

Page 11 of 15

The time-independent Schrödinger equation

The idea

In many situations, the potential energy of the particle does not depend on time. The forces acting on the particle, and hence VV, change only with the position of the particle, V=V(x)V=V(x). The particle is then in a steady state. When this is true, the Schrödinger equation can be simplified by separating the wave function into a part that depends only on position and a part that depends only on time. The result is the steady-state or time-independent form.

Derivation

Step 1: separate the wave function. Start from the free-particle wave function and split the exponential into two factors:

Ψ=A e−iℏ(Et−px)=A eiℏpx  e−iℏEt\Psi=A\,e^{-\frac{i}{\hbar}(Et-px)}=A\,e^{\frac{i}{\hbar}px}\;e^{-\frac{i}{\hbar}Et}

Let ψ=A eiℏpx\psi=A\,e^{\frac{i}{\hbar}px}, a function of position only (it is sometimes written Ψ0\Psi_0 or ψ0\psi_0). Then

Ψ=ψ e−iℏEt(1)\Psi=\psi\,e^{-\frac{i}{\hbar}Et}\qquad(1)

So the wave function Ψ\Psi has a time-independent part ψ\psi and a time-dependent part e−iEt/ℏe^{-iEt/\hbar}.

Step 2: differentiate (1) with respect to tt. ψ\psi does not depend on tt, so only the exponential is differentiated:

∂Ψ∂t=−iEℏ ψ e−iℏEt(2)\frac{\partial\Psi}{\partial t}=-\frac{iE}{\hbar}\,\psi\,e^{-\frac{i}{\hbar}Et}\qquad(2)

Step 3: differentiate (1) twice with respect to xx. The exponential does not depend on xx, so only ψ\psi is differentiated:

∂2Ψ∂x2=d2ψdx2 e−iℏEt(3)\frac{\partial^2\Psi}{\partial x^2}=\frac{d^2\psi}{dx^2}\,e^{-\frac{i}{\hbar}Et}\qquad(3)

Step 4: put in the time-dependent equation. The time-dependent Schrödinger equation is

iℏ ∂Ψ∂t=−ℏ22m ∂2Ψ∂x2+VΨ(4)i\hbar\,\frac{\partial\Psi}{\partial t}=-\frac{\hbar^2}{2m}\,\frac{\partial^2\Psi}{\partial x^2}+V\Psi\qquad(4)

Putting (1), (2) and (3) in (4):

iℏ(−iEℏ)ψ e−iℏEt=−ℏ22m d2ψdx2 e−iℏEt+Vψ e−iℏEti\hbar\left(-\frac{iE}{\hbar}\right)\psi\,e^{-\frac{i}{\hbar}Et}=-\frac{\hbar^2}{2m}\,\frac{d^2\psi}{dx^2}\,e^{-\frac{i}{\hbar}Et}+V\psi\,e^{-\frac{i}{\hbar}Et}

On the left, iℏ⋅(−iEℏ)=Ei\hbar\cdot\left(-\dfrac{iE}{\hbar}\right)=E, so

Eψ e−iℏEt=−ℏ22m d2ψdx2 e−iℏEt+Vψ e−iℏEtE\psi\,e^{-\frac{i}{\hbar}Et}=-\frac{\hbar^2}{2m}\,\frac{d^2\psi}{dx^2}\,e^{-\frac{i}{\hbar}Et}+V\psi\,e^{-\frac{i}{\hbar}Et}

Step 5: cancel the time factor. Divide throughout by e−iEt/ℏe^{-iEt/\hbar}, which is never zero:

Eψ=−ℏ22m d2ψdx2+VψE\psi=-\frac{\hbar^2}{2m}\,\frac{d^2\psi}{dx^2}+V\psi

Step 6: rearrange. Take everything to one side:

ℏ22m d2ψdx2+(E−V) ψ=0\frac{\hbar^2}{2m}\,\frac{d^2\psi}{dx^2}+(E-V)\,\psi=0

and multiply through by 2mℏ2\dfrac{2m}{\hbar^2}:

 d2ψdx2+2mℏ2 (E−V) ψ=0 \boxed{\,\frac{d^2\psi}{dx^2}+\frac{2m}{\hbar^2}\,(E-V)\,\psi=0\,}

This is the time-independent Schrödinger equation in one dimension. Since ψ\psi depends only on xx, there is no need for partial derivatives. Ordinary derivatives d2ψdx2\dfrac{d^2\psi}{dx^2} are used.

Three dimensions

∂2ψ∂x2+∂2ψ∂y2+∂2ψ∂z2+2mℏ2(E−V)ψ=0\frac{\partial^2\psi}{\partial x^2}+\frac{\partial^2\psi}{\partial y^2}+\frac{\partial^2\psi}{\partial z^2}+\frac{2m}{\hbar^2}(E-V)\psi=0 or∇2ψ+2mℏ2(E−V)ψ=0\text{or}\qquad\nabla^2\psi+\frac{2m}{\hbar^2}(E-V)\psi=0

Stationary states

For a steady state, the probability density is

∣Ψ∣2=∣ψ∣2 ∣e−iEt/ℏ∣2=∣ψ∣2|\Psi|^2=|\psi|^2\,\big|e^{-iEt/\hbar}\big|^2=|\psi|^2

because ∣e−iEt/ℏ∣=1|e^{-iEt/\hbar}|=1. So the probability density does not change with time. This is why these solutions are called stationary states. All the problems in this module (the box and the barrier) use the time-independent equation, because their potentials do not change with time.

Key idea. The time-independent equation can also be written as −ℏ22md2ψdx2+Vψ=Eψ-\dfrac{\hbar^2}{2m}\dfrac{d^2\psi}{dx^2}+V\psi=E\psi: some operator acting on ψ\psi gives a number EE times ψ\psi. Only certain values of EE (the energy eigenvalues) allow acceptable solutions ψ\psi (the eigenfunctions). The particle in a box shows this clearly.

Watch out. An equivalent and equally correct way to write the result is ℏ22md2ψdx2+(E−V)ψ=0\dfrac{\hbar^2}{2m}\dfrac{d^2\psi}{dx^2}+(E-V)\psi=0 (multiply through by ℏ2/2m\hbar^2/2m). Do not mix the two forms in the middle of a derivation.

Page 12 of 15

Particle in a one-dimensional box: solving the equation

The set-up

Consider a particle of mass mm moving inside a box along the xx direction, bouncing back and forth between two walls a distance LL apart. The walls at x=0x=0 and x=Lx=L are infinitely high (potential barriers that cannot be crossed), so the particle can never get past them, however much energy it has. This is the one-dimensional infinite square well potential, also called a particle in a one-dimensional box. The potential energy varies as:

V=0for 0<x<L(region II, inside the box)V=0\quad\text{for }0<x<L\qquad\text{(region II, inside the box)} V=∞for x≤0 and x≥L(regions I and III, outside)V=\infty\quad\text{for }x\le0\text{ and }x\ge L\qquad\text{(regions I and III, outside)}

Infinite square well: V = 0 between x = 0 and x = L, and V infinite outside (regions I, II and III)

The position of the particle at any moment lies in 0<x<L0<x<L. Since the particle cannot exist outside the box, the probability of finding it there is zero, so Ψ∗Ψ=0\Psi^*\Psi=0 and hence ψ=0\psi=0 for x≤0x\le0 and x≥Lx\ge L. Our task is to find ψ\psi inside the box. Since the potential does not depend on time, we use the time-independent Schrödinger equation.

The equation inside the box

The time-independent Schrödinger equation is

d2ψdx2+2mℏ2(E−V)ψ=0\frac{d^2\psi}{dx^2}+\frac{2m}{\hbar^2}(E-V)\psi=0

Inside the box V=0V=0, so

d2ψdx2+2mEℏ2 ψ=0\frac{d^2\psi}{dx^2}+\frac{2mE}{\hbar^2}\,\psi=0

Put

k2=2mEℏ2,that isk=2mEℏ2(1)k^2=\frac{2mE}{\hbar^2},\qquad\text{that is}\qquad k=\sqrt{\frac{2mE}{\hbar^2}}\qquad(1)

so the equation becomes

d2ψdx2+k2ψ=0(2)\frac{d^2\psi}{dx^2}+k^2\psi=0\qquad(2)

This is the equation of simple harmonic motion. Its general solution can be written in two equivalent ways.

Method A: sine and cosine form

The general solution of (2) is

ψ=Asin⁡kx+Bcos⁡kx(3)\psi=A\sin kx+B\cos kx\qquad(3)

with two constants AA and BB fixed by the boundary conditions. The wave function must be continuous, and it is zero outside the box, so it must vanish at both walls:

ψ=0  at x=0andψ=0  at x=L\psi=0\ \text{ at }x=0\qquad\text{and}\qquad\psi=0\ \text{ at }x=L

Boundary condition at x=0x=0:

ψ(0)=Asin⁡(k⋅0)+Bcos⁡(k⋅0)=A⋅0+B⋅1=B\psi(0)=A\sin(k\cdot0)+B\cos(k\cdot0)=A\cdot0+B\cdot1=B

So B=0B=0 and

ψ=Asin⁡kx\psi=A\sin kx

Boundary condition at x=Lx=L:

ψ(L)=Asin⁡kL=0\psi(L)=A\sin kL=0

We cannot have A=0A=0, because then ψ=0\psi=0 everywhere and there would be no particle. So

sin⁡kL=0⟹kL=nπ⟹k=nπL,n=1,2,3,…(4)\sin kL=0\qquad\Longrightarrow\qquad kL=n\pi\qquad\Longrightarrow\qquad k=\frac{n\pi}{L},\qquad n=1,2,3,\dots\qquad(4)

Method B: exponential form

The general solution of (2) can also be written as

ψ=Aeikx+Be−ikx\psi=Ae^{ikx}+Be^{-ikx}

At x=0x=0: ψ(0)=A+B=0\psi(0)=A+B=0, so B=−AB=-A and

ψ=A(eikx−e−ikx)=A⋅2isin⁡kx=2iAsin⁡kx\psi=A\left(e^{ikx}-e^{-ikx}\right)=A\cdot2i\sin kx=2iA\sin kx

using eikx−e−ikx=2isin⁡kxe^{ikx}-e^{-ikx}=2i\sin kx.

At x=Lx=L: 2iAsin⁡kL=02iA\sin kL=0. Since 2iA≠02iA\ne0, sin⁡kL=0\sin kL=0, so again kL=nπkL=n\pi and k=nπ/Lk=n\pi/L. The wave function is

ψ=2iAsin⁡nπxL\psi=2iA\sin\frac{n\pi x}{L}

which is the same as Method A with the constant C=2iAC=2iA in place of AA. Both methods lead to

ψn=Csin⁡nπxL,n=1,2,3,…\psi_n=C\sin\frac{n\pi x}{L},\qquad n=1,2,3,\dots

The wave function has several expressions, one for each value of nn. The nnth wave function is ψn\psi_n.

Why n=0n=0 is not allowed

For n=0n=0, k=0k=0 and ψ=Csin⁡0=0\psi=C\sin0=0 everywhere inside the box. A wave function that is zero everywhere means the particle is not present at all, so n=0n=0 is not allowed. Negative values of nn give the same functions with a minus sign (the same physical state), so we take n=1,2,3,…n=1,2,3,\dots

Watch out. The general solution has two constants. The boundary condition at x=0x=0 removes the cosine (or makes B=−AB=-A). The condition at x=Lx=L quantises kk. Do not drop a constant at the start.

Note on the source. One lecture set wrote the general solution as "Acos⁡kx+sin⁡kxA\cos kx+\sin kx", where the second constant is missing. It is corrected above.

Page 13 of 15

Particle in a box: energy levels and wave functions

Allowed energy levels

The value of nn depends on the energy of the particle. From (1) and (4) on the previous page,

k=2mEℏ2andkL=nπk=\sqrt{\frac{2mE}{\hbar^2}}\qquad\text{and}\qquad kL=n\pi ⟹2mEℏ2  L=nπ\Longrightarrow\qquad\sqrt{\frac{2mE}{\hbar^2}}\;L=n\pi

Since nn must be a whole number, the energy EE can take only certain values. Squaring both sides:

2mEℏ2 L2=n2π2\frac{2mE}{\hbar^2}\,L^2=n^2\pi^2  En=n2π2ℏ22mL2=n2h28mL2 n=1,2,3,…\boxed{\,E_n=\frac{n^2\pi^2\hbar^2}{2mL^2}=\frac{n^2h^2}{8mL^2}\,}\qquad n=1,2,3,\dots

The second form uses ℏ=h/2π\hbar=h/2\pi, so π2ℏ2=π2h2/4π2=h2/4\pi^2\hbar^2=\pi^2h^2/4\pi^2=h^2/4. Each allowed energy EnE_n is called an energy eigenvalue.

n=1: E1=h28mL2n=2: E2=4h28mL2=4E1n=3: E3=9h28mL2=9E1n=1:\ E_1=\frac{h^2}{8mL^2}\qquad n=2:\ E_2=\frac{4h^2}{8mL^2}=4E_1\qquad n=3:\ E_3=\frac{9h^2}{8mL^2}=9E_1

and in general En=n2E1E_n=n^2E_1: E4=16E1E_4=16E_1, E5=25E1E_5=25E_1, and so on. The number nn is called the quantum number.

Normalised wave functions

The wave function for level nn is ψn=Csin⁡nπxL\psi_n=C\sin\dfrac{n\pi x}{L}. The square of the absolute value of the wave function gives the probability of finding the particle. Since the particle is certainly somewhere inside the box,

∫0Lψn∗ψn dx=1(for all values of n)\int_0^L\psi_n^*\psi_n\,dx=1\qquad\text{(for all values of }n)

This is the normalisation condition. Putting ψn\psi_n:

CC∗∫0Lsin⁡2(nπxL)dx=1CC^*\int_0^L\sin^2\left(\frac{n\pi x}{L}\right)dx=1

Use sin⁡2θ=1−cos⁡2θ2\sin^2\theta=\dfrac{1-\cos2\theta}{2}:

∣C∣2∫0L1−cos⁡(2nπxL)2 dx=1|C|^2\int_0^L\frac{1-\cos\left(\frac{2n\pi x}{L}\right)}{2}\,dx=1

Integrate each term:

∫0L12 dx=L2∫0L12cos⁡2nπxL dx=L4nπ[sin⁡2nπxL]0L=L4nπ(sin⁡2nπ−sin⁡0)=0\int_0^L\frac12\,dx=\frac L2\qquad\qquad\int_0^L\frac12\cos\frac{2n\pi x}{L}\,dx=\frac{L}{4n\pi}\Big[\sin\frac{2n\pi x}{L}\Big]_0^L=\frac{L}{4n\pi}\left(\sin2n\pi-\sin0\right)=0

So

∣C∣2⋅L2=1⟹∣C∣2=2L⟹∣C∣=2L|C|^2\cdot\frac{L}{2}=1\qquad\Longrightarrow\qquad|C|^2=\frac2L\qquad\Longrightarrow\qquad|C|=\sqrt{\frac2L}

The normalised wave function is

 ψn=2L sin⁡nπxL n=1,2,3,…\boxed{\,\psi_n=\sqrt{\frac2L}\,\sin\frac{n\pi x}{L}\,}\qquad n=1,2,3,\dots

Each ψn\psi_n is an energy eigenfunction. For each value of nn there is an energy level EnE_n and a corresponding wave function ψn\psi_n.

Wave functions and probability densities for n = 1, 2, 3, and the energy levels E1, 4E1, 9E1, 16E1

The figure shows the first three wave functions ψn\psi_n (left), the corresponding probability densities ∣ψn∣2|\psi_n|^2 (middle) and the energy levels (right). As xx goes from 00 to LL, the angle nπxL\dfrac{n\pi x}{L} goes from 00 to nπn\pi, so ψn\psi_n contains exactly nn half-wavelengths: one half-wave for n=1n=1, two for n=2n=2, three for n=3n=3.

What the solution tells us

  1. Energy is quantised. Inside the potential well, the particle can have only the discrete energies EnE_n, not any arbitrary value. This is the origin of energy levels in atoms.
  2. Zero-point energy. The lowest energy is E1=h28mL2=π2ℏ22mL2E_1=\dfrac{h^2}{8mL^2}=\dfrac{\pi^2\hbar^2}{2mL^2}, not zero. This minimum is called the zero-point energy. A confined particle can never be completely at rest. (This agrees with the uncertainty principle: if the particle is confined to Δx≈L\Delta x\approx L, there is a spread in momentum and hence some kinetic energy.)
  3. The levels are not equally spaced. En∝n2E_n\propto n^2, so the gaps increase: E2−E1=3E1E_2-E_1=3E_1, E3−E2=5E1E_3-E_2=5E_1, E4−E3=7E1E_4-E_3=7E_1, …
  4. Preferred positions. In the ground state (n=1n=1) the probability of finding the particle is greatest at the centre (x=L/2x=L/2). For n=2n=2 the maxima are at L/4L/4 and 3L/43L/4, a quarter and three quarters of the way across, and the probability at the centre is zero (a node). For n=3n=3 the maxima are at L/6L/6, L/2L/2 and 5L/65L/6. Classically every position would be equally likely.
  5. Smaller box, larger energies. En∝1/L2E_n\propto1/L^2 and En∝1/mE_n\propto1/m. For a large box or a heavy particle, the levels are so close together that the quantisation cannot be noticed, and the classical picture is recovered.

Key idea. When a particle is confined in a potential well, we get (i) certain preferred positions of maximum probability, (ii) only certain discrete energy levels, and (iii) a minimum energy that the particle must have.

Worked examples

Example 1. Find the ground-state energy of an electron in an infinitely deep box of width 1 A˚1\ \text{Å}.

E1=h28mL2=(6.626×10−34)28×9.11×10−31×(10−10)2=6.02×10−18 JE_1=\frac{h^2}{8mL^2}=\frac{(6.626\times10^{-34})^2}{8\times9.11\times10^{-31}\times(10^{-10})^2}=6.02\times10^{-18}\ \text{J} E1=6.02×10−181.6×10−19 eV=37.6 eVE_1=\frac{6.02\times10^{-18}}{1.6\times10^{-19}}\ \text{eV}=\mathbf{37.6\ \text{eV}}

The next levels are E2=4E1=150 eVE_2=4E_1=150\ \text{eV} and E3=9E1=338 eVE_3=9E_1=338\ \text{eV}.

Example 2. In Example 1 the electron drops from n=2n=2 to n=1n=1 and emits a photon. Find the wavelength of the photon.

ΔE=E2−E1=3E1=3×37.6=112.8 eV\Delta E=E_2-E_1=3E_1=3\times37.6=112.8\ \text{eV} λ=hcΔE=1240 eV nm112.8 eV=11.0 nm\lambda=\frac{hc}{\Delta E}=\frac{1240\ \text{eV nm}}{112.8\ \text{eV}}=\mathbf{11.0\ \text{nm}}

(using hc=1240 eV nmhc=1240\ \text{eV nm}). This is in the far ultraviolet.

Example 3. Find the ground-state energy of an electron in a box of width 1 nm1\ \text{nm}.

The box is 10 times wider than in Example 1, and E1∝1/L2E_1\propto1/L^2, so the energy is 100 times smaller:

E1=37.6100=0.376 eVE_1=\frac{37.6}{100}=\mathbf{0.376\ \text{eV}}

Example 4. An electron in a box of width LL is in its first excited state. Find the probability of finding it in the central half of the box.

The first excited state is n=2n=2. The central half runs from L/4L/4 to 3L/43L/4:

P=∫L/43L/42Lsin⁡22πxL dx=1L∫L/43L/4(1−cos⁡4πxL)dx=1L[x−L4πsin⁡4πxL]L/43L/4P=\int_{L/4}^{3L/4}\frac2L\sin^2\frac{2\pi x}{L}\,dx=\frac1L\int_{L/4}^{3L/4}\left(1-\cos\frac{4\pi x}{L}\right)dx=\frac1L\left[x-\frac{L}{4\pi}\sin\frac{4\pi x}{L}\right]_{L/4}^{3L/4}

The sine term is sin⁡3π−sin⁡π=0\sin3\pi-\sin\pi=0, so

P=1L(3L4−L4)=12P=\frac1L\left(\frac{3L}{4}-\frac L4\right)=\mathbf{\frac12}

For comparison, in the ground state (n=1n=1) the same central half has probability 12+1π≈0.82\tfrac12+\tfrac1\pi\approx0.82, because the particle is concentrated near the centre.

Example 5. A particle is in the ground state of a box of width LL. Find the probability of finding it between x=0x=0 and x=L/4x=L/4.

P=∫0L/42Lsin⁡2πxL dx=1L[x−L2πsin⁡2πxL]0L/4=1L(L4−L2π)=14−12π=0.091P=\int_0^{L/4}\frac2L\sin^2\frac{\pi x}{L}\,dx=\frac1L\left[x-\frac{L}{2\pi}\sin\frac{2\pi x}{L}\right]_0^{L/4}=\frac1L\left(\frac L4-\frac{L}{2\pi}\right)=\frac14-\frac{1}{2\pi}=\mathbf{0.091}

Only about 9 %, much less than the 25 % a classical particle would have, because the ground-state wave function is small near the walls.

Watch out. The quantum number starts at n=1n=1, so the "first excited state" is n=2n=2, not n=1n=1. Also, the energy depends on n2n^2: if nn is doubled, the energy becomes four times.

Page 14 of 15

Quantum mechanical tunnelling

The idea

Consider a stream of particles, each of energy EE, falling on a potential barrier of height V0V_0 and width LL, with the energy less than the barrier height (E<V0E<V_0):

V=0for x<0(region I)V=0\quad\text{for }x<0\qquad\text{(region I)} V=V0for 0≤x≤L(region II, the barrier)V=V_0\quad\text{for }0\le x\le L\qquad\text{(region II, the barrier)} V=0for x>L(region III)V=0\quad\text{for }x>L\qquad\text{(region III)}
  • Classically, when a particle of lower kinetic energy approaches a higher potential barrier, it bounces back. It does not have enough energy to climb over, so it can never appear on the other side.
  • Quantum mechanically, the wave function representing the particle does not vanish at the barrier. If we solve the time-independent Schrödinger equation in the three regions, the wave function is non-zero in all three regions I, II and III. So there is a finite probability that the particle goes through the barrier and is found beyond it.

Quantum mechanical tunnelling. The phenomenon of particles passing through potential barriers higher than their own energy is known as quantum mechanical tunnelling. It is purely a quantum effect, showing the wave property of quantum objects, and has no classical counterpart.

Tunnelling: the incident wave is partly reflected, decays inside the barrier and emerges as a smaller transmitted wave

The wave function in each region

Region I (x<0x<0): incident and reflected waves. Here V=0V=0, so the solution is a travelling wave in each direction:

ψI=Aeikx+Be−ikx,k=2mEℏ\psi_{\text{I}}=Ae^{ikx}+Be^{-ikx},\qquad k=\frac{\sqrt{2mE}}{\hbar}

AA is the amplitude of the incident wave (moving right) and BB the amplitude of the reflected wave (moving left). The non-zero value of the wave function in region I is due to the incident particles and the particles reflected from the barrier.

Region II (0≤x≤L0\le x\le L): inside the barrier. Here E<V0E<V_0, so E−V0E-V_0 is negative and the solution no longer oscillates. It decays exponentially with distance into the barrier, roughly as

ψII≈Ce−κx,κ=2m(V0−E)ℏ\psi_{\text{II}}\approx Ce^{-\kappa x},\qquad\kappa=\frac{\sqrt{2m(V_0-E)}}{\hbar}

The wave function is small at the far side of the barrier, but not zero. Since the probability of finding the particle is non-zero in region II, the particle can be inside the barrier, which is classically forbidden.

Region III (x>Lx>L): the transmitted wave. Here V=0V=0 again, so in general

ψIII=Feikx+Ge−ikx\psi_{\text{III}}=Fe^{ikx}+Ge^{-ikx}

where FF is the amplitude of the transmitted wave. But nothing comes back from the right, so there is no reflected wave in region III: G=0G=0 and

ψIII=Feikx\psi_{\text{III}}=Fe^{ikx}

Interpretation. A non-zero wave function in region III means there is a finite probability of finding particles there. The incident particles have less energy than the barrier height, so they could not have climbed over it. Also, the wave function is non-zero inside the barrier. So we conclude that the particles have reached region III by penetrating, or "tunnelling" through, the barrier.

Transmission coefficient

The probability of tunnelling is measured by the transmission coefficient TT. Since the probability of finding the particle is proportional to the square of the amplitude,

T=probability density of the transmitted waveprobability density of the incident wave=∣F∣2∣A∣2T=\frac{\text{probability density of the transmitted wave}}{\text{probability density of the incident wave}}=\frac{|F|^2}{|A|^2}

In classical physics T=0T=0 when E<V0E<V_0. In quantum mechanics it is a finite, non-zero quantity. For a barrier that is not too thin, it is approximately

T≈e−2κL,κ=2m(V0−E)ℏT\approx e^{-2\kappa L},\qquad\kappa=\frac{\sqrt{2m(V_0-E)}}{\hbar}

so it falls off very rapidly:

  • The wider the barrier (LL larger), the smaller TT.
  • The higher the barrier compared with the energy (V0−EV_0-E larger), the smaller TT.
  • The heavier the particle (mm larger), the smaller TT. So tunnelling is important for light particles such as electrons.

Example 1. An electron has energy 1 eV1\ \text{eV} less than the height of a barrier. Find the approximate transmission coefficient for a barrier of width (a) 0.5 nm0.5\ \text{nm}, (b) 1 nm1\ \text{nm}.

κ=2×9.11×10−31×1×1.6×10−191.055×10−34=5.1×109 m−1\kappa=\frac{\sqrt{2\times9.11\times10^{-31}\times1\times1.6\times10^{-19}}}{1.055\times10^{-34}}=5.1\times10^{9}\ \text{m}^{-1}

(a) 2κL=2×5.1×109×0.5×10−9=5.12\kappa L=2\times5.1\times10^{9}\times0.5\times10^{-9}=5.1, so T≈e−5.1≈6×10−3T\approx e^{-5.1}\approx\mathbf{6\times10^{-3}}.

(b) 2κL=10.22\kappa L=10.2, so T≈e−10.2≈3.5×10−5T\approx e^{-10.2}\approx\mathbf{3.5\times10^{-5}}.

If the width is doubled, the tunnelling probability becomes about 170 times smaller. The scanning tunnelling microscope uses this extreme sensitivity to width.

Applications

  • Tunnel diode (Module 4): electrons tunnel through a very thin depletion layer, giving a negative-resistance region used in high-frequency oscillators.
  • Zener breakdown: in heavily doped p-n junctions, electrons tunnel through the thin depletion region at a sharp reverse voltage.
  • Alpha decay: alpha particles inside a nucleus have less energy than the nuclear barrier, yet they escape by tunnelling.
  • Scanning tunnelling microscope (STM): the tunnelling current between a sharp tip and a surface depends so strongly on the gap that it can map individual atoms.
  • Josephson effect and SQUIDs: Cooper pairs tunnel through a thin insulating layer between two superconductors.

Watch out. Tunnelling does not mean the particle gains energy to go over the barrier. Its energy stays EE throughout. It appears on the far side with the same energy, only with a much smaller probability (amplitude). Also, T≈e−2κLT\approx e^{-2\kappa L} is an approximation for a fairly wide barrier, not an exact result.

Page 15 of 15

Practice questions and formula sheet

Constants: h=6.626×10−34 J sh=6.626\times10^{-34}\ \text{J s}, ℏ=1.055×10−34 J s\hbar=1.055\times10^{-34}\ \text{J s}, me=9.11×10−31 kgm_e=9.11\times10^{-31}\ \text{kg}, e=1.6×10−19 Ce=1.6\times10^{-19}\ \text{C}, c=3×108 m/sc=3\times10^{8}\ \text{m/s}.

Numerical problems

1. A certain excited state of a helium atom has an average lifetime of 2.11×10−8 s2.11\times10^{-8}\ \text{s}. What is the minimum uncertainty with which the frequency of the emitted radiation can be measured?

Δν=14π Δt=14π×2.11×10−8=3.8×106 Hz=3.8 MHz\Delta\nu=\frac{1}{4\pi\,\Delta t}=\frac{1}{4\pi\times2.11\times10^{-8}}=\mathbf{3.8\times10^{6}\ \text{Hz}}=3.8\ \text{MHz}

2. Calculate the uncertainty in the frequency of the emitted radiation if the uncertainty in time of an excited atom is 5×10−8 s5\times10^{-8}\ \text{s}.

Δν=14π×5×10−8=1.6×106 Hz=1.6 MHz\Delta\nu=\frac{1}{4\pi\times5\times10^{-8}}=\mathbf{1.6\times10^{6}\ \text{Hz}}=1.6\ \text{MHz}

3. The time gap between the excitation of an atom and the emission of radiation is Δt=10−8 s\Delta t=10^{-8}\ \text{s}. Find the uncertainty in the frequency of the radiation.

Δν=14π×10−8=8.0×106 Hz=8.0 MHz\Delta\nu=\frac{1}{4\pi\times10^{-8}}=\mathbf{8.0\times10^{6}\ \text{Hz}}=8.0\ \text{MHz}

4. Find the energy of the ground state of an electron moving in an infinitely high potential box of width 1 A˚1\ \text{Å}.

E1=h28mL2=(6.626×10−34)28×9.11×10−31×10−20=6.02×10−18 J=37.6 eVE_1=\frac{h^2}{8mL^2}=\frac{(6.626\times10^{-34})^2}{8\times9.11\times10^{-31}\times10^{-20}}=6.02\times10^{-18}\ \text{J}=\mathbf{37.6\ \text{eV}}

5. An electron confined in a one-dimensional box of width LL is known to be in its first excited state. Determine the probability of finding the electron in the central half.

n=2n=2, and integrating ∣ψ2∣2=2Lsin⁡22πxL|\psi_2|^2=\frac2L\sin^2\frac{2\pi x}{L} from L/4L/4 to 3L/43L/4 gives P=12P=\mathbf{\tfrac12} (worked on the energy-levels page).

6. Find the de Broglie wavelength of an electron accelerated through 100 V100\ \text{V}.

λ=1.227×10−9100=1.23×10−10 m\lambda=\frac{1.227\times10^{-9}}{\sqrt{100}}=\mathbf{1.23\times10^{-10}\ \text{m}}

7. A bullet of mass 25 g25\ \text{g} moving at 400 m/s400\ \text{m/s} has its speed measured to 0.02 %0.02\,\%. Find the minimum uncertainty in its position.

Δv=0.08 m/s\Delta v=0.08\ \text{m/s}, Δp=0.025×0.08=2×10−3 kg m/s\Delta p=0.025\times0.08=2\times10^{-3}\ \text{kg m/s}, so Δx≥ℏ2Δp=2.64×10−32 m\Delta x\ge\dfrac{\hbar}{2\Delta p}=\mathbf{2.64\times10^{-32}\ \text{m}}.

8. Find the de Broglie wavelength of a neutron with kinetic energy 0.025 eV0.025\ \text{eV} (mn=1.675×10−27 kgm_n=1.675\times10^{-27}\ \text{kg}).

λ=h2mnKE=6.626×10−342×1.675×10−27×4.0×10−21=1.81×10−10 m\lambda=\frac{h}{\sqrt{2m_nKE}}=\frac{6.626\times10^{-34}}{\sqrt{2\times1.675\times10^{-27}\times4.0\times10^{-21}}}=\mathbf{1.81\times10^{-10}\ \text{m}}

9. An electron is confined to a region of size 10−10 m10^{-10}\ \text{m}. Find the minimum uncertainty in its speed.

Δp≥ℏ2Δx=5.3×10−25 kg m/s\Delta p\ge\dfrac{\hbar}{2\Delta x}=5.3\times10^{-25}\ \text{kg m/s}, so Δv=Δp/m=5.8×105 m/s\Delta v=\Delta p/m=\mathbf{5.8\times10^{5}\ \text{m/s}}.

10. For an electron in a box of width 1 A˚1\ \text{Å}, find the wavelength of the photon emitted in the transition n=2→n=1n=2\to n=1.

ΔE=3E1=112.8 eV\Delta E=3E_1=112.8\ \text{eV}, so λ=1240112.8 nm=11.0 nm\lambda=\dfrac{1240}{112.8}\ \text{nm}=\mathbf{11.0\ \text{nm}}.

Note on the source. The lecture sheets used ΔE Δt≈ℏ/2\Delta E\,\Delta t\approx\hbar/2 in one place and ℏ\hbar in another. With ℏ/2\hbar/2, Δν=1/(4πΔt)\Delta\nu=1/(4\pi\Delta t) as above. With ℏ\hbar every frequency width would be twice as large.

Questions asked in KTU exams

Short answer

  1. What is tunnel effect? What is quantum mechanical tunnelling?
  2. State and explain Heisenberg's uncertainty principle.
  3. What are the characteristics of a well-defined wave function?
  4. A well-behaved wave function is normalised. Why?
  5. Give the probability interpretation of the wave function.
  6. Explain the absence of electrons in the nucleus on the basis of the uncertainty principle.
  7. What is meant by a wave function? Write its normalisation condition.
  8. How do you account for natural line broadening on the basis of Heisenberg's uncertainty principle?

Essay

  1. Write the Schrödinger equation for a particle trapped in a one-dimensional box of width LL and solve it to obtain the energy eigenvalues and wave functions.
  2. State the uncertainty principle. Explain the absence of electrons inside the nucleus using this principle.
  3. Starting from the time-dependent equation, derive Schrödinger's time-independent wave equation.
  4. Formulate Schrödinger's time-dependent equation starting from a plane wave equation, using de Broglie's formula and Einstein's relation for photon energy.

Formula sheet

Matter waves and the wave function

de Broglie wavelengthλ=hp=hmv=h2m KE\lambda=\dfrac hp=\dfrac{h}{mv}=\dfrac{h}{\sqrt{2m\,KE}}
Electron through VV voltsλ=h2meV=1.227×10−9V m≈150V A˚\lambda=\dfrac{h}{\sqrt{2meV}}=\dfrac{1.227\times10^{-9}}{\sqrt V}\ \text{m}\approx\sqrt{\dfrac{150}{V}}\ \text{Å}
Free-particle wave functionΨ=Ae−iℏ(Et−px)\Psi=Ae^{-\frac{i}{\hbar}(Et-px)}, with ω=E/ℏ\omega=E/\hbar, k=p/ℏk=p/\hbar
Probability density∣Ψ∣2=Ψ∗Ψ\lvert\Psi\rvert^2=\Psi^*\Psi
Probability in a regionP=∫∣Ψ∣2 dVP=\displaystyle\int\lvert\Psi\rvert^2\,dV
Normalisation∫∣Ψ∣2 dV=1\displaystyle\int\lvert\Psi\rvert^2\,dV=1

Uncertainty principle

Position and momentumΔx Δp≥ℏ/2\Delta x\,\Delta p\ge\hbar/2
Energy and timeΔE Δt≥ℏ/2\Delta E\,\Delta t\ge\hbar/2
Angular momentum and angleΔL Δθ≥ℏ/2\Delta L\,\Delta\theta\ge\hbar/2
Natural line widthΔν=14π Δt\Delta\nu=\dfrac{1}{4\pi\,\Delta t}
Electron in nucleusΔx=10−14 m⇒KE≳10 MeV\Delta x=10^{-14}\ \text{m}\Rightarrow KE\gtrsim10\ \text{MeV}

Schrödinger equation

Time-dependentiℏ∂Ψ∂t=−ℏ22m∂2Ψ∂x2+VΨi\hbar\dfrac{\partial\Psi}{\partial t}=-\dfrac{\hbar^2}{2m}\dfrac{\partial^2\Psi}{\partial x^2}+V\Psi
Time-independentd2ψdx2+2mℏ2(E−V)ψ=0\dfrac{d^2\psi}{dx^2}+\dfrac{2m}{\hbar^2}(E-V)\psi=0
Three dimensions∇2ψ+2mℏ2(E−V)ψ=0\nabla^2\psi+\dfrac{2m}{\hbar^2}(E-V)\psi=0
OperatorsE→iℏ∂∂tE\to i\hbar\dfrac{\partial}{\partial t},  p→−iℏ∂∂x\ p\to-i\hbar\dfrac{\partial}{\partial x}

Particle in a box and tunnelling

Energy levelsEn=n2h28mL2=n2π2ℏ22mL2E_n=\dfrac{n^2h^2}{8mL^2}=\dfrac{n^2\pi^2\hbar^2}{2mL^2},  n=1,2,3,…\ n=1,2,3,\dots
Wave functionsψn=2Lsin⁡nπxL\psi_n=\sqrt{\dfrac2L}\sin\dfrac{n\pi x}{L}
Zero-point energyE1=h28mL2E_1=\dfrac{h^2}{8mL^2}
Level ratiosEn=n2E1E_n=n^2E_1: E1,4E1,9E1,16E1,…E_1, 4E_1, 9E_1, 16E_1,\dots
Transmission coefficientT=∣F∣2∣A∣2≈e−2κLT=\dfrac{\lvert F\rvert^2}{\lvert A\rvert^2}\approx e^{-2\kappa L},  κ=2m(V0−E)ℏ\ \kappa=\dfrac{\sqrt{2m(V_0-E)}}{\hbar}

You’ve reached the end of S1 CSE Physics · Module 2.

Source: KTU GAPHT121 study pack, Module 2: Quantum mechanics. Reorganised and rewritten from the lecture notes, with every worked example re-checked. Check the original material if a step looks unclear.