
Advanced Quantum Mechanics course sequence · Interactive lesson
The Casimir force between conducting plates
This is the one calculation in the course where a change in the field's zero-point energy accounts for a force you can put on a balance. We count the modes two mirrors allow, subtract two infinite sums the only honest way, arrive at 1.3 mPa across a micrometre, then ask what the cutoff really is and what the measured force does and does not prove.
UAQM-11.06 · Lesson 6 of 9 in this unit. Position is not a completion score.Build the model
Two mirrors change which modes the vacuum may hold, and a plate feels the difference.
The key idea
The Casimir force is the attraction between two neutral, parallel conducting plates in vacuum, π²ħc/(240d⁴) per unit area for ideal plates a distance d apart, which arises because the plates change the zero-point energies of the field modes between them.
Example. Ideal mirrors 1.00 µm apart are pressed together by 1.30 mPa; on a pair 1.0 cm × 1.0 cm that is 0.13 µN, the weight of 13 µg, and halving the gap multiplies it by 16.
The complete physical model and assumptions
Two parallel mirrors add nothing to the vacuum; they restrict it. Between perfectly conducting plates a distance d apart, the transverse field of topic 11.03 may only have kz = nπ/d, and every allowed mode is still an oscillator holding ½ħω. Moving a plate shifts the frequency of every mode with n ≥ 1, so Σ½ħω, which normal ordering could discard while the modes stayed put, now depends on d. The sum is infinite, and so is the free-space energy it replaces, but the two are infinite in the same way. Regulate both with one smooth cutoff and every term that grows without limit as the cutoff is removed is either proportional to d, the free-space energy density times the gap's volume, which the outside vacuum takes back as the plate moves, or independent of d, which pushes on nothing; the remaining cutoff-dependent terms vanish with the cutoff. What survives is finite and universal: −π²ħc/(720d³) per unit area, an attraction of π²ħc/(240d⁴), with no charge, no mass and no cutoff in it. The model stops where its idealisation does. A real metal reflects only below its plasma frequency, so the result needs d well above c/ωₚ, about 22 nm for gold, and at room temperature thermal photons take over beyond a few micrometres. It also stops short of a philosophical verdict: Lifshitz's theory reaches the same force from the fluctuating charges inside the plates, so the measured force confirms the quantised field, not the reality of its zero-point energy.
1The modes two mirrors allow
SI units throughout; the metric signature would be (+, −, −, −), though no four-vector is needed here. Put perfectly conducting plates of area A = L² at z = 0 and z = d, with L ≫ d and periodic conditions in x and y. At each plate the transverse field of topic 11.03, quoted, must have zero tangential E and zero normal B, so its z-dependence is sin(nπz/d) or cos(nπz/d) with kz = nπ/d and ω = c√(k∥² + n²π²/d²): derived. For n ≥ 1 two polarisations survive, TE with E parallel to the plates and TM with E in the plane of k and the normal. For n = 0 only TM does, a field Ez uniform across the gap; a TE field uniform in z would be tangential at the plates and must vanish. Each family holds A d²k∥/(2π)² modes. For d = 1.00 µm the n = 1 family starts at πc/d = 9.42 × 10¹⁴ s⁻¹, a free-space wavelength of 2d = 2.00 µm, yet the n = 0 family carries every longer wavelength too.
2Two infinite energies, and the difference a plate feels
Each mode is an oscillator with ground-state energy ½ħω, quoted from Unit 9. Summing over the modes above, with s = √(k∥² + n²π²/d²) as the variable, the gap holds (ħc/4π)Σ′ₙ F(nπ/d) per unit area, where F(κ) = 2∫ s² f(s/K) ds from κ to ∞ and Σ′ gives n = 0 half weight for its single polarisation: derived. The comparison is free space in the same volume A d, (ħc/4π)(d/π)∫ F(κ) dκ, and the reason is mechanical. Move a plate out by δd and the gap gains A δd of volume while the vacuum outside loses exactly that much, so the only energy that can change is the gap's minus free space's. Casimir's 1948 set-up, a plate inside a large conducting box, makes the same subtraction explicit. The cutoff f must be the same function on both sides, with f(0) = 1 and f → 0 well above K. With f = e(−λs) the free-space term is 3ħcd/(π²λ⁴) per unit area, 1.54 mJ m⁻² for λ = 50 nm and d = 1.00 µm, millions of times the answer we are after.
3Euler-Maclaurin turns the difference into one number
Write G(n) = F(nπ/d). The difference is a sum minus the integral of the same smooth function, and the Euler-Maclaurin formula, quoted, gives it as −G′(0)/12 + G‴(0)/720 − …, with the ½ on n = 0 exactly the trapezoid weight it needs. Differentiate F: F′(κ) = −2κ² f(κ/K), so F′(0) = F″(0) = 0 and F‴(0) = −4f(0) = −4, while F⁽⁵⁾(0) carries a factor 1/K². Then G‴(0) = (π/d)³F‴(0) = −4π³/d³ and ΔE/A = (ħc/4π)(−4π³/(720d³)) = −π²ħc/(720d³): derived. Only f(0) = 1 entered, which is why every smooth cutoff gives the same limit; that the remainder vanishes for each such f is stated, not proved here. A sharp cutoff fails that condition. With a step at |k| = K, the same difference swings from +480 to −302 times the Casimir value as Kd/π runs from 4 to 4.5, and its swings grow as K²: no limit exists.
4The pressure, checked three ways
Differentiate: P = −∂(ΔE/A)/∂d = −π²ħc/(240d⁴), derived, the minus sign meaning attraction. With π²ħc/240 = 1.300 × 10⁻²⁷ J m, P is 1.30 mPa at 1.00 µm, 13.0 Pa at 100 nm and 1.3 × 10⁵ Pa, about an atmosphere, at 10 nm, where the ideal model has already failed. Dimensions: J m ÷ m⁴ = J m⁻³ = Pa. Limits: as ħ → 0 it vanishes, as it must for a classical field at zero temperature, which has no zero-point motion; as d → ∞ it dies as d⁻⁴, faster than the d⁻² pressure of any residual voltage. Measured, quoted: Lamoreaux (1997), torsion pendulum, sphere and plate, 0.6 to 6 µm, agreement to about 5%; Mohideen and Roy (1998), atomic-force microscope, 0.1 to 0.9 µm, about 1%; Bressi and co-workers (2002), parallel plates themselves, 0.5 to 3 µm, the coefficient to 15%. Sphere-plate data are compared through F = 2πR ΔE/A, the medium example.
5What the cutoff stands for
The cutoff is not an embarrassment; it is the metal. A real conductor reflects well only below its plasma frequency ωₚ and turns transparent above it, so modes much shorter than c/ωₚ never feel the plates. For gold ħωₚ ≈ 9.0 eV and c/ωₚ ≈ 22 nm. The ideal result is the limit d ≫ c/ωₚ, where the modes that decide the force, with wavelengths near d, are reflected almost perfectly. The figure shows what happens away from it: two smooth cutoffs that agree as λ/d → 0 split in opposite directions as λ grows, so the departure belongs to the metal, not the vacuum. Lifshitz's theory (1956), quoted, computes it from the metal's dielectric function; for a plasma-model metal P ≈ P₀[1 − 16δ₀/(3d) + …] with δ₀ = c/ωₚ, about 10% below ideal for gold at 1 µm. Also quoted: at gaps of a few tens of nanometres the force crosses over to the unretarded van der Waals form ∝ d⁻³, and at room temperature thermal photons become comparable near 3 µm, where whether to treat the metal with dissipation or without is still argued between experiments.
6What the measured force does and does not show
The derivation used only how Σ½ħω changes with d, never Σ½ħω itself, so it says nothing about the absolute zero-point energy that gravity would feel. Nor does the force alone prove the zero-point energy real. Lifshitz's theory, and later Schwinger's source theory, obtain the same pressure as the retarded van der Waals attraction between fluctuating charges and currents in the two plates, never mentioning the field's zero-point energy, a point Jaffe (2005) made explicit; it is named here, not derived. The coupling gives it away. A real metal has ωₚ² = ne²/(ε₀m), which vanishes with e, and then so does the force; the ideal result, free of e, is the limit in which the plates respond infinitely strongly at every relevant frequency. The two accounts are one physics in two bookkeepings: modes of the field with boundary conditions, or correlated fluctuations of the charges that impose them. What the measurements do establish is that the quantised field gets the force between neutral bodies right where retardation matters.
Key equations
Read the units, assumptions and limits before substituting.
Modes between ideal plates
ω = c√(k∥² + n²π²/d²), n = 0, 1, 2, …
Derived from zero tangential E at both plates. The n = 0 family fits every wavelength, so the gap is never empty of long waves.
k∥ in-plane wavevector (m⁻¹), d plate gap (m), ω in s⁻¹; two polarisations for n ≥ 1, one (TM, E along z) for n = 0.Regulated energy difference per unit area
ΔE/A = (ħc/4π)[Σ′ₙ F(nπ/d) − (d/π)∫₀^∞ F(κ) dκ]
Derived. Both pieces grow as K⁴d; only their difference changes the total energy when a plate moves, because the outside vacuum takes back the volume the gap gives up.
F(κ) = 2∫ s² f(s/K) ds from s = κ to ∞, in m⁻³; cutoff f with f(0) = 1, K in m⁻¹; Σ′ halves the n = 0 term; ΔE/A in J m⁻².Euler-Maclaurin for a half-weighted sum
Σ′ₙ G(n) − ∫₀^∞ G(n) dn = −G′(0)/12 + G‴(0)/720 − G⁽⁵⁾(0)/30240 + …
Quoted from analysis. The ½ on n = 0 is exactly the trapezoid weight it needs, so no G(0) term survives and a single derivative carries the answer.
G(n) = F(nπ/d), smooth and decaying; here G′(0) = 0 and G‴(0) = −4π³/d³, while later terms fall as 1/K².Casimir energy per unit area
ΔE/A = −π²ħc/(720 d³)
Derived for ideal plates as the cutoff is removed, and stated to be the same for every smooth cutoff. Neither e nor K appears: only ħ, c and d.
ħc = 3.1615 × 10⁻²⁶ J m, d in m, ΔE/A in J m⁻²; −4.33 × 10⁻¹⁰ J m⁻² at d = 1.00 µm.Casimir pressure between ideal plates
P = −∂(ΔE/A)/∂d = −π²ħc/(240 d⁴)
Derived. Halving the gap multiplies it by 16; Bressi and co-workers (2002) measured the coefficient of d⁻⁴ between parallel plates to 15%.
P in Pa, negative meaning attraction; π²ħc/240 = 1.300 × 10⁻²⁷ J m gives 1.30 mPa at 1.00 µm and 13.0 Pa at 100 nm.Finite conductivity, plasma-model metal
P ≈ P₀[1 − 16δ₀/(3d) + 24δ₀²/d²]
Quoted from Lifshitz's theory. For gold ħωₚ ≈ 9.0 eV gives δ₀ ≈ 22 nm, and at 1 µm the pressure is about 10% below the ideal value.
δ₀ = c/ωₚ (m), the depth the field enters the metal; ωₚ the plasma frequency (s⁻¹); P₀ the ideal pressure; valid for d ≫ δ₀.Interactive lab
Predict first. Explore the model. Compare the result.
Explore the model
Try this. Start at λ = 50 nm and d = 1000 nm, where both cutoffs match the ideal value within 0.4%. Slide d to 200 nm: the markers ride out to where the curves fan apart, 0.915 and 1.078, as the pressure climbs 625-fold. Then cut λ to 5 nm: both curves flatten onto 1, as the derivation promised.
Live results
- IDEAL PRESSURE π²ħc/240d⁴
- 1.300 mPa
- CASIMIR ENERGY PER AREA
- -0.4334 nJ m⁻²
- SUBTRACTED 3ħcd/π²λ⁴
- 1.538 mJ m⁻²
- LOG₁₀ (SUBTRACTED ÷ CASIMIR)
- 6.55
Results follow this lesson’s model and assumptions. Displayed values are rounded.
Predict → change one variable → check
Plan a change before running it. Every other input stays at the saved baseline.
Saved baseline. Cutoff length λ: 50 nm; Plate separation d: 1000 nm
This experiment is kept only while this page is open. It does not change saved course progress.
Catch the common trap
The usual picture is that the gap shuts out every wave longer than 2d, so the outside vacuum, holding all its modes, pushes harder than the inside one, and the force could be estimated by adding up ½ħω for the missing long waves. The calculation does not work that way. Long waves are not shut out: the n = 0 family fits any wavelength, and counting every mode below a frequency ω, the gap holds more modes than the same volume of free space for every ω below 3πc/(4d). Above that the counts trade places, the gap ahead at every half-integer value of ωd/(πc) and behind at every integer, with swings that grow with ω, so any count-and-stop estimate depends on where it stops: the sharp-cutoff failure again. The attraction comes from comparing the whole sum with the whole integral through a smooth cutoff, and its sign is fixed by the G‴(0)/720 term, not by a tally of absent modes.
Practice & worked examples
3 authored examples. Make an attempt, then reveal the reasoning one step at a time.
EasyExample 1
Two ideal, uncharged conducting plates, each 1.0 cm × 1.0 cm, face each other across a vacuum gap of 0.50 µm. Find the Casimir pressure and the force on one plate. Then find the potential difference V between the plates whose electrostatic pressure ε₀V²/(2d²) would equal it, and say what that means for an experiment.
Sketch your answer or reasoning, then compare it with the authored solution. Your writing is not automatically graded or saved.
0 of 4 steps revealed.
- Derived in this lesson: P = π²ħc/(240d⁴). With π²ħc/240 = 1.300 × 10⁻²⁷ J m and d⁴ = (5.0 × 10⁻⁷ m)⁴ = 6.25 × 10⁻²⁶ m⁴, P = 2.08 × 10⁻² Pa = 20.8 mPa, attractive. Dimensions: J m ÷ m⁴ = J m⁻³ = Pa.
- Force on A = 1.0 × 10⁻⁴ m²: F = PA = 2.1 × 10⁻⁶ N, the weight of about 0.21 mg. Because P ∝ d⁻⁴, the same plates 1.0 µm apart feel 1/16 of this, 1.3 × 10⁻⁷ N.
- Equal electrostatic pressure: ε₀V²/(2d²) = P gives V = d√(2P/ε₀) = 5.0 × 10⁻⁷ m × √(2 × 2.08 × 10⁻² Pa ÷ 8.854 × 10⁻¹² F m⁻¹) = 5.0 × 10⁻⁷ m × 6.85 × 10⁴ V m⁻¹ = 34 mV. Check: (F m⁻¹)(V²)/m² = J m⁻³ = Pa.
- The two laws scale differently, ε₀V²/(2d²) as d⁻² and the Casimir pressure as d⁻⁴, so at 1.0 µm the same 34 mV gives 5.2 mPa, four times the Casimir pressure. Residual potentials of tens of millivolts are routine between real metal surfaces, so every measurement since 1997 has had to find and cancel this voltage before it could see the Casimir force.
Answer. P = 20.8 mPa (attractive) and F = 2.1 × 10⁻⁶ N on 1.0 cm²; a residual 34 mV would match it at 0.50 µm and outweigh it at any larger gap.
MediumExample 2
A metal-coated sphere of radius R = 150 µm hangs d = 0.40 µm above a large flat plate, both treated as ideal conductors. In the proximity-force approximation each ring of the sphere interacts as if it faced a parallel plate at its own local gap. Derive the force from ΔE/A = −π²ħc/(720d³), evaluate it and its gradient, and find the radius of the patch that supplies 7/8 of the force.
Sketch your answer or reasoning, then compare it with the authored solution. Your writing is not automatically graded or saved.
0 of 6 steps revealed.
- A ring of radius ρ on the sphere sits at local gap s = d + ρ²/(2R) for ρ ≪ R, so ds = ρ dρ/R and the interaction energy is U(d) = ∫ 2πρ ΔE/A(s) dρ = 2πR ∫ ΔE/A(s) ds over s from d to ∞.
- F = −dU/dd = 2πR ΔE/A(d) = −π³ħcR/(360d³): derived within the approximation, whose relative error is of order d/R (quoted), here 3 × 10⁻³. The minus sign is attraction.
- Numbers: π³ħc = 31.006 × 3.1615 × 10⁻²⁶ J m = 9.803 × 10⁻²⁵ J m. Times R = 1.50 × 10⁻⁴ m and divided by 360d³ = 360 × 6.40 × 10⁻²⁰ m³ = 2.304 × 10⁻¹⁷ m³, |F| = 6.38 × 10⁻¹² N. Dimensions: J m × m ÷ m³ = J m⁻¹ = N.
- Gradient: F ∝ d⁻³, so |dF/dd| = 3|F|/d = 3 × 6.38 × 10⁻¹² N ÷ 4.0 × 10⁻⁷ m = 4.8 × 10⁻⁵ N m⁻¹, the quantity a dynamic measurement reads from the shift of an oscillator's frequency.
- Rings out to local gap s₀ supply the fraction [ΔE/A(d) − ΔE/A(s₀)]/ΔE/A(d) = 1 − (d/s₀)³, which is 7/8 at s₀ = 2d, that is at ρ₀ = √(2Rd) = √(2 × 1.50 × 10⁻⁴ m × 4.0 × 10⁻⁷ m) = 1.1 × 10⁻⁵ m. The force comes from a patch of radius 11 µm, and no two surfaces need be held parallel, which is why most measurements use this geometry.
- Limits: as R → ∞ the force grows without bound, the infinite-plate case, and as d → ∞ it falls as d⁻³, one power slower than the plate pressure because F = P(d) × 2πRd/3, an effective area 2πRd/3 that grows with d.
Answer. |F| = π³ħcR/(360d³) = 6.4 × 10⁻¹² N (6.4 pN), attractive, with gradient 4.8 × 10⁻⁵ N m⁻¹; 7/8 of it comes from within ρ₀ = √(2Rd) = 11 µm of the closest point.
HardExample 3
With the cutoff f = e(−λk) the regulated difference can be summed exactly, without Euler-Maclaurin. Using F(κ) = 2∫ s² e(−λs) ds from κ to ∞, evaluate Σ′ₙ F(nπ/d) − (d/π)∫₀^∞ F(κ) dκ in closed form and expand it in a = πλ/d. For d = 1.00 µm and λ = 50 nm, find the free-space term the subtraction removes, the Casimir energy with its first cutoff correction, and the change in the pressure.
Sketch your answer or reasoning, then compare it with the authored solution. Your writing is not automatically graded or saved.
0 of 6 steps revealed.
- Integrating by parts, F(κ) = (2/λ³)e(−λκ)(λ²κ² + 2λκ + 2), so G(n) = F(nπ/d) = (2/λ³)qⁿ(a²n² + 2an + 2) with q = e(−a). With S₀ = Σqⁿ = 1/(1 − q), S₁ = Σnqⁿ = q/(1 − q)² and S₂ = Σn²qⁿ = q(1 + q)/(1 − q)³, the sum is (2/λ³)(a²S₂ + 2aS₁ + 2S₀ − 1), the −1 removing half of G(0) = 4/λ³, and the integral is (2/λ³)(6/a).
- Expand for small a: 2S₀ = 2/a + 1 + a/6 − a³/360 + a⁵/15120, 2aS₁ = 2/a − a/6 + a³/120 − a⁵/3024 and a²S₂ = 2/a − a³/120 + a⁵/1512. The total is 6/a + 1 − a³/360 + a⁵/2520 + …: the 6/a cancels the integral, the +1 cancels the half-weight −1, and the terms in a cancel among themselves.
- Multiplying the remainder (2/λ³)(−a³/360 + a⁵/2520) by ħc/4π, with a³/λ³ = π³/d³, gives ΔE/A = −(π²ħc/720d³)(1 − π²λ²/(7d²) + …): the Euler-Maclaurin result again, with its first cutoff correction, which vanishes as λ → 0. The subtracted free-space term is (ħc/4π)(2/λ³)(6/a) = 3ħcd/(π²λ⁴).
- Numbers: 3ħcd/(π²λ⁴) = 9.485 × 10⁻³² J m² ÷ (9.870 × 6.25 × 10⁻³⁰ m⁴) = 1.54 × 10⁻³ J m⁻². The correction is π²λ²/(7d²) = 3.52 × 10⁻³, so ΔE/A = −4.334 × 10⁻¹⁰ × (1 − 0.00352) J m⁻² = −4.318 × 10⁻¹⁰ J m⁻²: the answer is 1 part in 3.6 × 10⁶ of the term subtracted.
- Pressure: from ΔE/A = −C/d³ + Cπ²λ²/(7d⁵) with C = π²ħc/720, P = −3C/d⁴ + 5Cπ²λ²/(7d⁶) = −(π²ħc/240d⁴)(1 − 5π²λ²/(21d²)), 0.59% weaker: 1.292 mPa instead of 1.300 mPa.
- The cutoff (1 + λk)e(−λk) gives the same leading term with the correction's sign reversed, as the figure shows. So the size of the departure belongs to the cutoff chosen, not to nature; a real metal's departure, first order in c/(ωₚd) in Lifshitz's theory, has to be quoted from a model of the metal.
Answer. ΔE/A = −(π²ħc/720d³)(1 − π²λ²/7d²) = −4.318 × 10⁻¹⁰ J m⁻², against a subtracted 1.54 × 10⁻³ J m⁻² (3.6 × 10⁶ times larger); the pressure is 0.59% below ideal, 1.292 mPa.
Quick check
Select an answer and explain it before checking.
Review & continue
The Casimir force is the attraction between two neutral, parallel conducting plates in vacuum, π²ħc/(240d⁴) per unit area for ideal plates a distance d apart, which arises because the plates change the zero-point energies of the field modes between them.
- Explain the relationship without looking at the equation.
- Check a sign, a unit and a model assumption in your working.
- Return to the lab or worked example if a step is still unclear.