At a glance
A steady beam of particles with energy E meets a rectangular barrier of height U and full width 2a, and the sim plots the fraction T that gets through, for every energy at once. You change the barrier and the particle and watch how T responds.
What you'll learn
- How the classical and the wave picture of the barrier compare.
- What the wave does inside the barrier, and how to read T and R = 1 − T.
- How T depends on the barrier's width and height and on the particle's mass.
- What sets the special energies where T = 1, and where they can occur.
What's on screen
- Transmission T (blue curve): the fraction that gets through.
- Classical limit (dashed grey): the classical prediction.
- Reflection R = 1 − T (pink curve, off until ticked).
- E = U boundary (pink vertical line): the barrier top on the energy axis.
- Resonance markers (magenta dots on T = 1, labelled n = 1, 2, …).
- Scanning marker (amber dot): sets the energy that Current Values reports.
- Current Values: T, R, E, k and η or κ at the marker.
- Barrier diagram (under the plot): a schematic of the waves at the marker's energy.
The controls
Press ⌖ on a row to see where that control sits in the sim.
| Control | What it changes | Range | |
|---|---|---|---|
| View mode | Plots T against the energy E, or against the wave number k = √(2mE). | T vs E / T vs k opens on T vs E | |
| Half-width a | The barrier's half-width. The full width is 2a. | 0.05–2.5 Å opens at 1.00 | |
| Height U | The height of the barrier top. | 0.5–15 eV opens at 5.00 | |
| Mass m | The mass of the particle, in electron masses. | 0.1–1.0 mₑ opens at 1.00 | |
| Scanner speed | How fast the amber marker sweeps along the curve. | 0.25–1.00× opens at 1.00 | |
| Current Values | T, R, E, k and η or κ at the scanning marker. With the scanner paused, they sit at the middle of the energy or k axis. | panel | |
| The plot | Hover over the curve to read T, E, k and η or κ at that point. | E from 0 to 3U | |
| Reflection R (Display) | Adds the pink curve R = 1 − T. | off at the start | |
| Classical prediction (Display) | Shows or hides the dashed grey line. | on at the start | |
| Region shading (Display) | Tints the plot left of E = U ("Tunneling") and right of it ("Propagating"). | on at the start | |
| Resonance markers (Display) | Shows or hides the magenta dots at T = 1. | on at the start | |
| Animated scanner (Display) | Turns the sweeping marker on or off. It is the same switch as Start / Pause. | on at the start | |
| Show barrier diagram (Display) | Shows or hides the schematic under the plot. | on at the start | |
| Start / Pause | Runs or freezes the scanning marker, so you can read T at one energy. | top bar | |
| Reset | Returns a, U, m and the scanner speed to their opening values. | top bar |
Fixed: ħ = 1. The sim reads the slider numbers directly, so Å and eV are labels only; real units would need ħ²/2mₑ = 3.81 eV·Å².
A wave at a wall
A classical particle with E < U is turned back at the near edge, and with E > U it passes. Classically T is 0 below the top and 1 above: the dashed grey step. The sim instead treats the beam as a steady wave ψ(x) that solves the time-independent Schrödinger equation, with k = √(2mE)/ħ outside the barrier.
- ψ and its slope must be continuous at x = −a and x = +a. A wave that is exactly zero beyond x = +a could not match this slope, so some wave always gets through.
- T is the transmitted fraction of the beam and R the reflected fraction. Nothing is absorbed, so R = 1 − T.
- T and R are fractions of a steady beam. The plot does not follow a single particle across the wall.
Decay inside the barrier
For E < U the equation inside the barrier loses its oscillating solutions:
- The amplitude of the decaying part falls as e−ηx and its probability density as e−2ηx.
- The barrier is only 2a thick, so the decay is cut short. What is left at x = +a leaves as the transmitted wave.
- Above the top the inside oscillates, with κ in place of η. Since κ < k, its wavelength is longer than outside.
- The barrier diagram is a schematic. The reflected wave is not drawn, the interior curve simply falls to √T of the incident height, and nothing in it is a particle sitting inside the barrier.
What the width does
Below the top the sim uses the exact result, with k and η as above:
When 2ηa ≫ 1, sinh grows like e2ηa/2 and the result simplifies:
- The exponent 4ηa is the density decay e−2ηx carried across the full width 2a. The width adds to it directly.
- At E = U/2 with the opening settings the exponent is 8.94 and T = 5.2 × 10⁻⁴. At a = 2 Å it is 17.9 and T = 6.8 × 10⁻⁸.
- The estimate needs 2ηa ≫ 1. At a = 0.35 Å the exponent is only 3.13, and the estimate gives 0.175 against the exact 0.160.
These sliders move the sim too. The numbers come from the sim's own formula (ħ = 1).
What the mass and height do
Mass and height reach T through η = √(2m(U − E))/ħ, which sits under a square root inside the exponent.
- Mass. Quartering m from 1.00 to 0.25 mₑ halves η, so the exponent 4ηa goes from 8.94 to 4.47. T at E = U/2 rises from 5.2 × 10⁻⁴ to 4.5 × 10⁻², about 86 times.
- Height. Raising U from 5 to 10 eV with E held at U/2 multiplies η by √2. The exponent becomes 12.6 and T falls to 1.3 × 10⁻⁵.
- Mass and width together. For a given E and U, T depends on m and a only through m a². So a = 1.00 Å at m = 0.25 mₑ gives the same T as a = 0.50 Å at m = 1.00 mₑ.
Above the top: perfect transmission and dips
For E > U the sinh turns into a sine, and a sine has zeros:
- Each edge of the barrier reflects part of the wave, even where the potential steps down. When a whole number of half-wavelengths of the interior wave fits in 2a, that is 2κa = nπ, the two reflections cancel: R = 0 and T = 1.
- At the opening settings the markers sit at E₁ = 6.23 eV and E₂ = 9.93 eV. Between them T falls to 0.76 near 7.4 eV, and the next dip, near 12.3 eV, only reaches 0.94.
- Below the top sinh² is positive, so T < 1 at every E < U. Only a very thin wall comes close: T = 0.95 at U/2 for a = 0.05 Å.
Where it matters
- Scanning tunnelling microscope. Electrons tunnel across a vacuum gap against a work function of about 4 eV, so η ≈ 1.0 Å⁻¹. The current falls roughly eightfold for every extra ångström of gap, which is why single atoms can be resolved.
- Alpha decay. An α particle leaks through a Coulomb barrier. A factor of about 2 in decay energy (8.95 MeV for ²¹²Po, 4.08 MeV for ²³²Th) spans about 24 orders of magnitude in half-life, from 0.3 μs to 10¹⁰ years.
- Fusion in the Sun. Core thermal energies are around 1 keV, far below the proton–proton Coulomb barrier of around 1 MeV. The Sun burns because a tail of protons tunnels.
These use SI constants and textbook values, not the sim's ħ = 1 numbers. Real barriers are not rectangles.
Go deeper: a barrier of any shape
- This is the WKB estimate, good when the potential changes slowly compared with the wavelength.
Summary
Key points
- Below the top T ≈ 16(E/U)(1 − E/U) e−4ηa for a thick wall, with η = √(2m(U − E))/ħ.
- Width enters the exponent directly; mass and height enter through a square root.
- Above the top T = 1 at 2κa = nπ, with dips between.
Common mix-ups
- "A thin enough wall would let a classical particle through too." The classical line ignores a.
- "The particle borrows energy to climb over." E is the same on both sides; what crosses is wave amplitude.
- "Halve the mass and T doubles." The mass sits inside the exponent: a quarter of the mass gives 86 times the T.
- "Above the top, everything gets through." T dips below 1 between the markers, because the edges reflect.
Things to try on your own
- Raise a to 2 Å and count the resonance markers, then do the same at a = 0.5 Å.
- Hover over the first magenta marker and check that E₁ − U = π²/(8ma²) for your settings.