← SciSim / Chemistry

§1 Interactive Simulation

Three.js r128
E_n (eV)
0.038
n (level)
1
L (nm)
1.0
Nodes
0
λ (nm)
2.0
ΔE→n+1 (eV)
0.114
⚙ Controls
Quantum Number n
1110
Box Length L | nm
0.11.020
Particle mass (×mₑ)
111836
nx (for 2D/3D)
115
ny (for 2D/3D)
115
Animation Speed
0.25×
Show nodes
Show probability
Animate phase

§2 The Idea, Step by Step

Pluck a guitar string and you don't get just any pitch — you get a fundamental note and its overtones, never anything in between. A jump rope is the same: you can shake it into one hump, or two, or three, but never two-and-a-half. An electron trapped between two hard walls behaves exactly like that jump rope. Only certain standing-wave shapes "fit," so the electron can only have certain energies. Energy comes in steps, not a smooth ramp — and that is the whole meaning of the word quantized.

Three things set the size of those steps: the level number $n$ (1 for the fundamental shape, 2 for the next, and so on), the box length $L$, and the particle's mass $m$. The allowed energies are

$$E_n = \frac{n^2 h^2}{8mL^2}$$

where $h$ is Planck's constant. Two big lessons hide in that one line: energy climbs like $n^2$, so the gaps get wider as you go up; and energy shrinks like $1/L^2$, so a bigger box gives lower, more tightly packed levels. For an electron in a $1\text{ nm}$ box the lowest rung is already $E_1 = 0.38\text{ eV}$ — small, but not zero. Even the ground state is never perfectly still.

Going deeper, each level carries a wave shape $\psi_n(x)=\sqrt{2/L}\,\sin(n\pi x/L)$, and the thing you can actually measure is $|\psi_n|^2$ — the probability of finding the particle at each point. Level $n$ has exactly $n-1$ nodes, spots where the particle is never found, just like the still points on a vibrating string. The "$n^2$ ladder" is forced by the walls: the wave must drop to zero at both ends, so only whole numbers of half-wavelengths can fit. In the sim, the $n$ slider picks the shape, the $L$ slider stretches the box, and the mass slider lets you swap the electron for something heavier (a proton is about $1836\times$ heavier).

Try this in the sim above: (1) Slide $n$ from 1 to 5 and count the nodes — you always get $n-1$. (2) Open the "E vs L" graph and drag $L$ wider; watch every level slide down toward zero. (3) Push the mass slider up: the heavier the particle, the smaller the energy gaps, until the quantum steps blur into the smooth classical world.

§3 Equation Derivation

Particle in a Box — 1D Schrödinger Equation Solution
\[E_n = \frac{n^2 h^2}{8mL^2} \qquad \psi_n(x)=\sqrt{\frac{2}{L}}\sin\!\left(\frac{n\pi x}{L}\right)\]
SymbolMeaningUnit
\(E_n\)Energy of level nJ or eV
\(n\)Quantum number = 1,2,3,…dimensionless
\(h\)Planck's constant = 6.626×10⁻³⁴J·s
\(m\)Particle masskg
\(L\)Box lengthm
\(\psi_n\)Wavefunction (normalized)m⁻¹/²
\(|\psi_n|^2\)Probability densitym⁻¹
Step-by-Step Derivation
Step 1 — Set up the Schrödinger equation Inside box (0 ≤ x ≤ L), V = 0: \(-\dfrac{\hbar^2}{2m}\dfrac{d^2\psi}{dx^2} = E\psi\)
Outside box, V = ∞, so ψ = 0.
Step 2 — General solution \(\psi(x) = A\sin(kx) + B\cos(kx)\), where \(k = \sqrt{2mE}/\hbar\)
Step 3 — Apply boundary conditions ψ(0) = 0 → B = 0. ψ(L) = 0 → A sin(kL) = 0 → kL = nπ → \(k_n = n\pi/L\)
Step 4 — Quantized energy \(E_n = \dfrac{\hbar^2 k_n^2}{2m} = \dfrac{n^2\pi^2\hbar^2}{2mL^2} = \dfrac{n^2h^2}{8mL^2}\)
Key: energy ∝ n², ∝ 1/L², ∝ 1/m. Ground state E₁ > 0 (zero-point energy!)
Step 5 — Normalize the wavefunction \(\int_0^L |\psi|^2 dx = 1 \Rightarrow A = \sqrt{2/L}\)
\(\psi_n(x) = \sqrt{2/L}\,\sin(n\pi x/L)\) — n-1 nodes (zeros) inside box
Step 6 — 2D and 3D extensions \(E_{n_x,n_y,n_z} = \dfrac{h^2}{8m}\!\left(\dfrac{n_x^2}{L_x^2}+\dfrac{n_y^2}{L_y^2}+\dfrac{n_z^2}{L_z^2}\right)\) Degenerate levels occur when n_x²/L_x² + n_y²/L_y² + n_z²/L_z² are equal by different (nx,ny,nz)
Worked Example — Electron in 1nm Box

m = mₑ = 9.109×10⁻³¹ kg, L = 1.0×10⁻⁹ m, n = 1

\(E_1 = \dfrac{(6.626\times10^{-34})^2}{8\times9.109\times10^{-31}\times(10^{-9})^2} = 6.02\times10^{-20}\text{ J} = 0.376\text{ eV}\)

Energy spacing E₂ − E₁ = E₁(4−1) = 3×0.376 = 1.13 eV → photon λ = hc/ΔE = 1100 nm (near-IR)

Reference: Atkins & de Paula — Physical Chemistry, 11th Ed., §7C "Translational Motion: The Particle in a Box"

§4 Frequently Asked Questions

🔬 SimulationWhat does each simulation mode show?
The 1D Wavefunction mode shows ψn(x) = √(2/L)sin(nπx/L) as an animated oscillating wave. The 1D Probability mode shows |ψn|² — where the particle is most likely to be found, with n−1 nodes (zero probability regions). 2D Box renders a color-coded |ψ|² heatmap for (nx,ny) states, revealing nodal planes. 3D Box shows an isosurface visualization. Superposition shows what happens when multiple states are mixed — the probability oscillates in time (wave packet). Key: The wavefunction ψ is not directly observable — only |ψ|² (probability density) has physical meaning.
🌍 Real LifeWhere is the particle-in-a-box model used in real chemistry?
The particle-in-a-box model explains the color of conjugated dye molecules like β-carotene (the orange pigment in carrots). The π electrons are delocalized over the conjugated chain of length L; using the PIB formula gives absorption wavelengths that match experiment surprisingly well. Quantum dots (semiconductor nanocrystals used in QLED TVs) work on this principle — smaller dots = smaller L = higher energy = bluer color. The model also explains why metals conduct electricity (free electron model) and provides the foundation for understanding band gaps in semiconductors. Key: QLED TVs, solar cells, and biological pigments all rely on the quantum confinement described by particle-in-a-box physics.
💡 Non-ObviousWhy can't the particle have zero energy (why is E₁ > 0)?
If E = 0, then ψ = 0 everywhere — the particle doesn't exist. More fundamentally, the Heisenberg uncertainty principle requires Δx·Δp ≥ ℏ/2. Confining the particle to a box of length L means Δx ≤ L, so Δp ≥ ℏ/2L > 0. Nonzero Δp means nonzero kinetic energy. This minimum energy E₁ = h²/8mL² is called the zero-point energy — even at absolute zero, the particle has kinetic energy. This is not a thermal effect but a quantum confinement effect. Key: Zero-point energy is a direct consequence of the uncertainty principle — confinement in space forces nonzero momentum, hence nonzero energy.
🧮 MathematicalHow do you calculate the wavelength of light absorbed by a PIB system?
Calculate ΔE = E_{n+1} − E_n = (2n+1)h²/(8mL²), then λ = hc/ΔE. Example: electron in 2nm box, n=1→2 transition: ΔE = 3h²/(8mL²) = 3×(6.626e-34)²/(8×9.11e-31×(2e-9)²) = 4.52×10⁻²⁰ J = 0.282 eV. λ = hc/ΔE = (6.626e-34 × 3e8)/(4.52e-20) = 4400 nm (mid-IR). Longer box → smaller ΔE → longer wavelength absorbed — this is why longer conjugated molecules absorb at longer (redder) wavelengths. Key: Larger box = smaller energy gaps = longer wavelength absorption — this directly explains the color tuning of quantum dots.
🧪 ConceptualWhat do the nodes in the wavefunction physically mean?
Nodes are positions where ψ = 0, which means |ψ|² = 0 — the particle has zero probability of being found there. For the nth level, there are n−1 nodes inside the box. The ground state (n=1) has no nodes — the particle is most likely to be found in the middle of the box. This seems paradoxical: for n=2, there are two regions where the particle can be found (left and right halves), separated by a node at x=L/2, yet the particle must somehow "cross" this node without being there. This is resolved by recognizing the particle doesn't have a classical trajectory. Key: More nodes = higher energy = more oscillations of the wavefunction — this is the quantum analogue of higher harmonics in a standing wave on a string.
🌌 DeepHow does the PIB model connect to molecular orbital theory?
In a linear conjugated molecule (like ethylene, butadiene, hexatriene), the π electrons move along the carbon chain — a nearly one-dimensional potential well. The PIB wavefunctions approximate the molecular orbitals: lowest MO has no nodes (like ψ₁), next has one node (like ψ₂), etc. Hückel MO theory refines this. The HOMO-LUMO gap corresponds to ΔE in the PIB, which determines the absorption wavelength. This connects directly to why π→π* transitions shift to longer wavelengths in longer conjugated systems. Key: PIB is the simplest model for molecular orbitals in conjugated systems — it provides the physical intuition behind MO theory.

Reference: LibreTexts Chemistry — "Particle in a Box" https://chem.libretexts.org/Bookshelves/Physical_and_Theoretical_Chemistry_Textbook_Maps/Physical_Chemistry_(LibreTexts)/03%3A_The_Schrödinger_Equation_and_a_Particle_in_a_Box | MIT OCW 5.61 Physical Chemistry https://ocw.mit.edu/courses/5-61-physical-chemistry-fall-2017/

§5 Common Misconceptions

❌ Misconception: "The particle bounces back and forth like a classical ball inside the box."
✅ Correction: The particle does not have a trajectory — it exists as a standing wave described by ψn(x). The particle doesn't "move" from one end to the other; its wavefunction is simultaneously spread throughout the box. The probability of finding the particle at any position x is |ψn(x)|² dx. Classical particle-like behavior (bouncing) is only recovered in the limit of very large quantum numbers (n → ∞), called the correspondence principle.
📖 Reference: Atkins & de Paula — Physical Chemistry, 11th Ed., §7C.1 — "The particle does not have a definite position"
❌ Misconception: "n=0 is the ground state with zero energy."
✅ Correction: n=0 is forbidden: ψ₀ = √(2/L)sin(0) = 0 everywhere, meaning the particle has zero probability of existing — it's not a physical state. The minimum allowed quantum number is n=1, giving the zero-point energy E₁ = h²/8mL² > 0. This is fundamentally different from classical mechanics, where a particle at rest has E=0.
📖 Reference: McQuarrie & Simon — Physical Chemistry: A Molecular Approach, Chapter 3.2: "Quantization Conditions" — n=0 excluded explicitly
❌ Misconception: "Larger boxes always have lower energies at every level."
✅ Correction: For a given n, En ∝ 1/L² — yes, larger L gives lower En. However, all levels decrease, not just the lowest. Importantly, for a macroscopic box (L=1m), E₁ ≈ 6×10⁻³⁸ J, which is so tiny that quantum effects vanish — this is why we don't see quantum confinement in everyday objects. Only when L is on the nanometer scale do quantum energy gaps become chemically meaningful (≥ 0.01 eV).
📖 Reference: Levine — Physical Chemistry, 6th Ed., Chapter 2.4: "The Particle in a Box" — box size comparison table
❌ Misconception: "The node means the particle can never cross that point — so the two halves are isolated regions."
✅ Correction: The node at x=L/2 (n=2 state) means the probability of detecting the particle there is zero, but this doesn't mean the particle is localized in only one half. The particle's wavefunction exists on both sides simultaneously — quantum mechanics forbids asking "which side is it on" without making a measurement. This is quantum superposition: the n=2 state is a coherent combination spanning the entire box.
📖 Reference: Atkins & de Paula — Physical Chemistry, 11th Ed., §7C.2: "The interpretation of wavefunctions"
❌ Misconception: "The particle in a box is just a toy model with no real applications."
✅ Correction: The PIB model quantitatively predicts the absorption spectra of conjugated dye molecules to within 10–20% with just one parameter (L). Quantum dots used in QLED displays (Samsung, LG) are explicitly engineered using quantum confinement. The free electron model for metals (PIB in 3D) explains metallic conductivity and optical properties. The PIB also forms the foundation for the Schrödinger equation solutions used in all of molecular orbital theory.
📖 Reference: Skoog et al. — Fundamentals of Analytical Chemistry, 9th Ed., Chapter 6 on spectroscopy | Quantum dot QLED technology

Section 4 reference: Taber, K.S. — Chemical Misconceptions (RSC, 2002) | Tsaparlis, G. & Papaphotis, G. — Quantum chemical concepts: are they helpful to students? Chem. Educ. Res. Pract. 2002, 3, 129