Skip to content

    Under active development. New calculators and features ship every week. Everything currently runs free while the site matures — pricing pages are a preview, checkout is not open yet.

    Engineering / Quantum Tunneling Probability Density Engine

    Quantum Tunneling Calculator

    Compute exact quantum tunneling transmission and reflection through a rectangular barrier for electrons, protons and alpha particles. Free, in-browser.

    ToolCalculatorExecutionIn browserCategoryEngineering
    ZorvLabs · Quantum tunneling

    Quantum Tunneling Probability Density Engine

    Exact finite square barrier scattering state solver with current-balance checks, probability-density profiles, and an energy-scan benchmark.
    Probability balance
    1.00000000
    R + T = 1.00000000
    Transmission
    0.04271914
    Exact current ratio
    Reflection
    0.95728086
    Exact current ratio
    T + R
    1.0000000000
    Probability conservation
    Barrier branch
    Below barrier
    Tunneling regime
    State controls
    View mode
    Particle preset
    Barrier height V0 (eV)
    12.00 eV
    Barrier width L (nm)
    0.180 nm
    Incident energy E (eV)
    6.00 eV
    Left x padding (nm)
    1.20 nm
    Right x padding (nm)
    1.20 nm
    Benchmark
    MIT OCW tunneling example: E = 6 eV, V0 = 12 eV, L = 0.18 nm → exact T = 0.042719.
    Closed-form T = 0.04271914
    Opaque-barrier WKB estimate = 0.04365672
    Source check uses the finite-barrier transmission relations and boundary matching.
    k (free region)
    12.549146 nm⁻¹
    Barrier wave number
    12.549146 nm⁻¹
    Current balance error
    0.0000000000
    Mass ratio m/mₑ
    1.000000
    2D scattering-state profile
    Potential barrier & probability density
    V(x) and normalized |ψ(x)|²0.0003.0006.0009.00012.000-1.20-0.68-0.170.350.861.38V(x)|ψ|² (scaled)
    Max density
    3.914086
    Max Re ψ
    1.399022
    Max Im ψ
    1.398920
    Transmission spectrum
    T(E) across the barrier
    Exact T(E) vs opaque-barrier WKB estimate0.0000.2500.5000.7501.0000.007.2014.4021.6028.8036.00T exactWKB approx
    Energy scan max
    36.00 eV
    Exact T at E
    0.04271914
    Exact R at E
    0.95728086
    3D surface view
    |ψ(x;E)|² probability-density surface
    3D probability density surface |ψ(x;E)|²x (nm)E (eV)|ψ|²surface grid42 × 68 cells · max |ψ|² = 3.99518
    Rows
    42
    Cols
    68
    Energy rows max
    36.00 eV
    Surface cap
    3.995177
    Equations used
    Schrödinger equation: iħ ∂Ψ/∂t = [−(ħ²/2m)∂²/∂x² + V(x)]Ψ
    Separation of variables: Ψ(x,t)=ψ(x)e^(−iEt/ħ)
    Probability density: ρ(x,t)=|Ψ(x,t)|²
    Boundary conditions: ψ continuous and ψ′ continuous at each interface
    Finite square barrier transmission: T = 1 / [1 + (V0²/(4E(V0−E)))sinh²(κL)] for E<V0, T = 1 / [1 + (V0²/(4E(E−V0)))sin²(qL)] for E>V0
    Top-of-barrier limit: T = 1 / (1 + k²L²/4) when E = V0
    Opaque-barrier estimate (optional overlay): T ≈ 4e^(−2κL)
    Exact matching solution
    Complex amplitudes from continuity at the interfaces
    QuantityRealImagMagnitude
    r0.00000000-0.978407310.97840731
    A0.010796350.010796350.01526834
    B0.98920365-0.989203651.39894522
    t-0.13125228-0.159662090.20668610
    Source notes
    The code comments and equations follow the standard Schrödinger formulation, the continuity conditions at the boundaries, the finite square barrier transmission coefficient, and the tunneling interpretation of finite barriers.
    Sources used for the equation set: Rutgers QM review for the Schrödinger equation, Born probability density, boundary continuity, and rectangular-barrier transmission; ETSU quantum-scattering notes for the piecewise finite-barrier wavefunctions and reflection/transmission definitions; MIT OCW for the opaque-barrier tunneling estimate.
    engineering methodology

    A clear path from inputs to results.

    Quantum particles can pass through a barrier they classically could not climb. This tool computes that tunneling probability for a particle hitting a rectangular potential barrier: enter the particle energy, the barrier height and width, and the particle type, and it returns the exact transmission and reflection probabilities from matching the Schrödinger-equation wavefunction across the barrier.

    Step-by-step methodWorked exampleAssumptions explained

    How to use

    1. Choose the particle: electron, proton, or alpha particle (this sets the mass used in the Schrödinger equation).
    2. Enter the incident energy E of the particle in electron-volts.
    3. Enter the barrier height V0 in electron-volts and the barrier width L in nanometres.
    4. Read the transmission probability T (the chance the particle tunnels through) and the reflection probability R = 1 − T.
    5. Compare against the built-in benchmark and the WKB approximation to see how the exact result behaves as you make the barrier taller, wider, or the energy lower.

    The formula

    T=[ 1+V02 sinh⁡2(κL)4E(V0−E) ]−1,κ=2m(V0−E)ℏ2,E<V0T = \left[\,1 + \dfrac{V_0^{2}\,\sinh^{2}(\kappa L)}{4E(V_0-E)}\,\right]^{-1}, \quad \kappa = \sqrt{\dfrac{2m(V_0-E)}{\hbar^{2}}}, \quad E < V_0

    Reference: Griffiths, Introduction to Quantum Mechanics, finite rectangular barrier transmission (exact boundary-matching solution); WKB opaque-barrier limit T ≈ 16·E(V0−E)/V0²·e^(−2κL).

    Worked example

    Electron through a 12 eV, 0.18 nm barrier at 6 eV

    Inputs
    • Particle: electron (mass ratio 1); incident energy E = 6 eV; barrier height V0 = 12 eV; barrier width L = 0.18 nm. The engine uses ħ²/2mₑ = 0.038099821161327 eV·nm².
    Calculation path
    1. 1

      Since E = 6 eV is below the 12 eV barrier, the wave is evanescent inside it with decay constant κ = √((V0−E)/(ħ²/2mₑ)) = √(6 / 0.0380998) = 12.5491 nm⁻¹. The dimensionless barrier thickness is κL = 12.5491 × 0.18 = 2.25885, giving sinh(κL) = 4.73378. The exact transmission is T = 1 / [1 + V0²·sinh²(κL) / (4E(V0−E))] = 1 / [1 + 144 × 22.4087 / (4 × 6 × 6)] = 1 / [1 + 22.4087] = 0.042719. Here V0²/(4E(V0−E)) = 144/144 = 1 exactly, so this also equals 1/cosh²(κL), an independent check that returns the same 0.042719.

    Result

    Transmission T = 0.04272 (about a 4.3% chance of tunneling through) and reflection R = 0.95728. The opaque-barrier WKB estimate 4·e^(−2κL) = 0.04366 is close, confirming the barrier is in the moderately opaque regime.

    Interpreting the result

    T is the probability that a particle in the incident beam emerges on the far side; R = 1 − T is the probability it is reflected, and the two always sum to one because the barrier neither creates nor absorbs particles. Tunneling falls off roughly exponentially with the barrier's opacity κL, so making the barrier a little wider or taller, or lowering the energy, drops T sharply — which is why scanning tunneling microscopes and alpha decay are so sensitive to distance and energy. When E exceeds V0 the particle is classically allowed over the barrier, and T oscillates toward 1 with resonances rather than decaying.

    Assumptions & limitations

    • The barrier is a one-dimensional rectangular step of constant height V0 and finite width L with flat potential on both sides; real barriers are smooth and this is an idealization.
    • The particle is a non-relativistic, spin-free plane wave described by the time-independent Schrödinger equation.
    • Energies and the barrier height are measured in electron-volts and the width in nanometres; the mass enters only through ħ²/2m, scaled by the particle's mass ratio to the electron.
    • The closed-form transmission applies for E < V0; the engine switches to the exact over-barrier solution for E ≥ V0, where resonant transmission appears.
    • The result is a steady-state scattering probability, not a time-resolved tunneling dynamics or a tunneling-time estimate.

    Frequently asked questions

    What does the transmission probability actually mean?

    It is the fraction of incident particles that make it through the barrier. A transmission of 0.04272 means that out of every thousand electrons that hit this 12 eV, 0.18 nm barrier at 6 eV, about 43 tunnel through and the other 957 are reflected. For a single electron it is the probability of being found on the far side.

    How does changing the barrier width affect tunneling?

    Very strongly. Transmission falls off roughly as e^(−2κL), so widening the barrier even slightly reduces T sharply. Doubling the 0.18 nm width to 0.36 nm in the worked example would cut the transmission by roughly the square of the current factor. This exponential sensitivity to distance is exactly what makes the scanning tunneling microscope able to resolve single atoms.

    Does the particle type change the answer?

    Yes, through the mass. The decay constant κ = √(2m(V0−E)/ħ²) grows with mass, so heavier particles tunnel far less. A proton (about 1836 times the electron mass) or an alpha particle (about 7294 times) sees a much larger κL for the same barrier and tunnels with a dramatically smaller probability than an electron — one reason electron tunneling dominates in solid-state devices.

    Why can a particle cross a barrier taller than its energy?

    In quantum mechanics the particle is a wave, and its wavefunction does not stop abruptly at the barrier — it decays exponentially inside the classically forbidden region. If the barrier is thin enough, a small but non-zero amplitude survives to the far side, giving a finite probability of transmission. This has no classical analogue; a classical ball with less energy than the barrier always bounces back.

    How does the WKB approximation compare to the exact result?

    The WKB opaque-barrier estimate 4·e^(−2κL) gives 0.04366 against the exact 0.04272 in this example — close but not identical. WKB is a semiclassical approximation that becomes accurate only when the barrier is thick and opaque (large κL). For thin or near-transparent barriers it can be off by a large factor, so the tool reports the exact boundary-matching result as authoritative.

    Methodology last reviewed 2026-10-06. Results are informational and not financial or engineering advice — see our disclaimer.