
Advanced Quantum Mechanics course sequence · Interactive lesson
Vacuum fluctuations: ⟨E²⟩ and its spectrum
Topic 11.03 made E an operator; here we ask what it does in the state with no photons. Its mean is zero, its variance is not, and the variance at a point is infinite. We compute that variance and its spectrum, find where the infinity comes from, and sort what an instrument can register from what it never will.
UAQM-11.05 · Lesson 5 of 9 in this unit. Position is not a completion score.Build the model
No photons and no mean field, yet E has a variance, with a spectrum rising as ω³.
The key idea
Vacuum fluctuations of the electric field are its nonzero variance ⟨0|E⋅E|0⟩ in the photon vacuum, where the mean ⟨0|E|0⟩ is zero: each transverse mode contributes ħω/(2ε₀V), so the variance at a point diverges while any smooth average over a region is finite.
Example. One mode of wavelength 500 nm in a 1.00 cm³ box has a vacuum r.m.s. field of 0.150 V m⁻¹, yet the 6.7 × 10¹³ modes of that box with wavelengths of 500 nm or longer add in quadrature to 1.06 × 10⁶ V m⁻¹ at a point.
The complete physical model and assumptions
Topic 11.03 wrote the transverse field as a sum over modes of √(ħω/(2ε₀V)) times a combination of a and a†, and Unit 9 showed that every mode of a free field is an oscillator with ground-state energy ½ħω. In the photon vacuum the mean field vanishes, because ⟨0|a|0⟩ and ⟨0|a†|0⟩ both do, but the variance does not: in ⟨0|E⋅E|0⟩ the product a a† survives, and each mode and polarisation contributes ħω/(2ε₀V), which holds the electric half of that mode's zero-point energy. With the ω² density of modes the sum becomes a spectrum, ħω³/(2π²ε₀c³) per unit angular frequency: the input Quantum Mechanics II assumed for Welton's Lamb-shift estimate, now derived. Its integral grows as the fourth power of whatever cutoff is imposed, because the free field has no length of its own to stop it. That infinity belongs to the field at a mathematical point, which no instrument samples: weight the field with a Gaussian of size L and the variance is ħc/(π²ε₀L⁴), finite and cutoff-independent once the cutoff passes a few times 1/L. Two things are never seen directly. The constant Σ½ħω cancels from every energy difference, so normal ordering removes it at no cost unless the mode frequencies themselves move, as between the plates of topic 11.06; and a detector that works by absorption reads the normally ordered ⟨E⁽⁻⁾⋅E⁽⁺⁾⟩, which is zero in the vacuum however large ⟨E²⟩ is.
1Zero mean, nonzero variance
SI units throughout; where a four-vector appears the metric is (+, −, −, −). Quoted from topic 11.03: E(x, t) = i Σₖₛ √(ħωₖ/(2ε₀V)) [aₖₛêₖₛ e(ik⋅x − iωₖt) − aₖₛ†êₖₛ* e(−ik⋅x + iωₖt)], with unit polarisation vectors êₖₛ ⊥ k in a periodic box of volume V. Quoted from Unit 9: each mode is an oscillator whose ground state has energy ½ħω. Since aₖₛ|0⟩ = 0 and ⟨0|aₖₛ† = 0, ⟨0|E|0⟩ = 0: derived. In ⟨0|E⋅E|0⟩ every product with an a on the right or an a† on the left vanishes, leaving ⟨0|aₖₛaₖ'ₛ'†|0⟩ = δₖₖ'δₛₛ', and êₖₛ⋅êₖₛ* = 1 gives Σₖₛ ħωₖ/(2ε₀V), the same at every x and t: derived. The same steps give c²⟨B²⟩ = ⟨E²⟩. Check: one mode's electric energy (ε₀V/2)⟨E²⟩ₖₛ is ħω/4, its magnetic energy another ħω/4, total ½ħω, Unit 9's zero point. In the one-mode limit the variance is ħω/(2ε₀V), falling as 1/V: a bigger box gives each mode a smaller share and holds more modes to share it.
2From a sum over modes to a spectrum
A periodic box holds V d³k/(2π)³ wavevectors in d³k; with two polarisations and ω = ck that is Vω² dω/(π²c³) modes between ω and ω + dω. Each brings ħω/(2ε₀V), so ⟨E²⟩ = ∫ dω ħω³/(2π²ε₀c³): derived. Dimensions: ħω³ dω/(ε₀c³) is (J s⁻³)/(F m² s⁻³) = J/(F m²) = V² m⁻², and V has cancelled, as it must for a local quantity. The ω² mode density is the one in Planck's law, which the COBE FIRAS spectrum of the microwave background confirms to 50 parts per million of its peak; times ε₀, the spectrum is exactly the ½ in Planck's (n̄ + ½)ħω³/(π²c³). This is the zero-point spectrum Quantum Mechanics II assumed for Welton's estimate, now a consequence of quantising the field; Unit 13 puts it to use. Its shape matters more than its normalisation. The share of the variance below ω grows as ω⁴, so half of everything up to a cutoff ωc comes from above 0.841ωc: whatever cutoff is chosen controls the answer.
3Why the sum diverges, and what a cutoff stands for
Keep the modes up to ωc = cK and the integral gives ⟨E²⟩ = ħcK⁴/(8π²ε₀), derived. The number of modes grows as ω² and each one's share as ω, so nothing stops it. Dimensions say it had to. ħ, c and ε₀ combine into ħc/ε₀, in V² m², with no length left over, so the variance at a point can only be ħc/ε₀ times the fourth power of a wavenumber supplied from outside. A crystal has such a number, its lattice spacing, which ends the phonon spectrum at the Debye frequency; the free Maxwell field has none. Each halving of the shortest wavelength kept multiplies the r.m.s. field by four, from 2.65 × 10⁵ V m⁻¹ with every mode down to 1 µm to 4.5 × 10¹⁶ V m⁻¹ down to the electron's Compton wavelength, √(α/2π) times the Schwinger field m²c³/(eħ). A cutoff is not a property of the vacuum. It stands for whatever limits the question, the size of a probe, the bandwidth of a detector or the scale where the free-field model stops, and an answer that depends on it asked about the field at a point.
4Average over a region and the infinity goes away
No instrument reads E at a point. Let one register Ē = ∫f(x)E(x) d³x with a Gaussian weight f = (πL²)(−3/2) e(−|x|²/L²), normalised to 1. Each mode is then multiplied by f's Fourier transform, e(−k²L²/4), so ⟨Dz⟩ is ħc/(2π²ε₀) times the integral of k³ e(−k²L²/2) over k from 0 to ∞, which is 2/L⁴: ⟨Dz⟩ = ħc/(π²ε₀L⁴), derived, finite, and free of any cutoff. Modes with kL ≫ 1 oscillate many times across the weight and cancel; the weighted spectrum peaks at k = √3/L. With a sharp cutoff K as well, the result carries the factor 1 − (1 + K²L²/2)e(−K²L²/2), which passes 99% at KL = 3.6 and then stops changing. Limits: L → 0 returns the K⁴ divergence, and L → ∞ gives zero, as it must for a field with no mean averaged over all space. For L = 1.00 µm the r.m.s. field is 1.90 × 10⁴ V m⁻¹. Electro-optic sampling with femtosecond pulses (Riek et al., 2015) measured this variance directly and found it scaling inversely with the space-time volume sampled, as ħc/(ε₀L⁴) does when all four extents grow together.
5Normal ordering, and the constant no difference sees
The divergent piece is a number, not an operator. Normal ordering, every a† moved left of every a, writes E⋅E = :E⋅E: + ⟨0|E⋅E|0⟩ with ⟨0|:E⋅E:|0⟩ = 0; applied to the Hamiltonian it drops E₀ = Σ½ħω, as Unit 9 did. Nothing measured is lost. Every photon energy, every atomic transition frequency and every rate in this course is a difference between two states of the same field, and the constant cancels from it. It survives only where the set of mode frequencies changes with a parameter: move a conducting plate and every ωₖ shifts, so Σ½ħωₖ changes by an amount that, regulated, is finite and independent of the cutoff. That difference is the Casimir pressure, topic 11.06's calculation, measured to a few per cent since 1997. Gravity responds to absolute energy and so would see E₀ itself. Cut off at the Compton wavelength, ħcK⁴/(8π²) is already 1.8 × 10²² J m⁻³, against a measured dark-energy density near 5 × 10⁻¹⁰ J m⁻³; this course leaves that mismatch where Unit 9 left it.
6What a detector registers, and what it cannot
Split E = E⁽⁺⁾ + E⁽⁻⁾, with E⁽⁺⁾ holding the a's. In the vacuum ⟨0|E⋅E|0⟩ = ⟨0|E⁽⁺⁾⋅E⁽⁻⁾|0⟩: the whole variance is anti-normally ordered. A detector that works by absorption, an atom or a photodiode, gains energy, so by energy conservation in the golden rule, quoted, the field must lose it, which only E⁽⁺⁾ can arrange. Its rate is therefore proportional to ⟨E⁽⁻⁾⋅E⁽⁺⁾⟩ (Glauber, quoted), and E⁽⁺⁾|0⟩ = 0 makes it zero: a ground-state atom in the vacuum is never excited, and no photodiode counts zero-point photons. The variance shows itself wherever no energy is taken from the field: as the noise floor of a balanced homodyne detector with its signal port blocked, against which squeezed light has been measured 15 dB below (Vahlbruch et al., 2016); in level shifts such as the 1057.8 MHz Lamb shift, Unit 13; in forces between moving boundaries, topic 11.06; and in spontaneous emission, which Unit 12 derives.
Key equations
Read the units, assumptions and limits before substituting.
Vacuum variance as a sum over modes
⟨0|E⋅E|0⟩ = Σₖₛ ħωₖ/(2ε₀V)
Derived: only ⟨0|aₖₛaₖₛ†|0⟩ = 1 survives, and ⟨0|E|0⟩ = 0. Each term is one mode's vacuum variance, the value topic 11.09 uses for a single cavity mode.
E in V m⁻¹; ωₖ = c|k| in s⁻¹; V the box volume in m³; s = 1, 2 labels the two polarisations transverse to k; ε₀ in F m⁻¹.Zero-point spectral density
d⟨E²⟩/dω = ħω³/(2π²ε₀c³)
Derived: the input Welton's estimate assumed in Quantum Mechanics II. Times ε₀ it is the ½ in Planck's (n̄ + ½)ħω³/(π²c³), the zero-point energy per unit volume per unit ω.
Per unit angular frequency, in V² m⁻² s; it counts Vω²/(π²c³) modes per unit ω, both polarisations included.Variance at a point, with a sharp cutoff
⟨0|E⋅E|0⟩ = ħωc⁴/(8π²ε₀c³) = ħcK⁴/(8π²ε₀)
Derived, and divergent as K⁴: keeping every mode down to 100 nm already gives Eᵣₘₛ = 2.65 × 10⁷ V m⁻¹, and half the variance comes from the top 16% of the range.
ωc = cK is the highest angular frequency kept, K in m⁻¹; ħc/ε₀ = 3.571 × 10⁻¹⁵ V² m².Zero-point energy of one mode, split evenly
(ε₀V/2)⟨E²⟩ₖₛ + (V/2μ₀)⟨B²⟩ₖₛ = ħωₖ/4 + ħωₖ/4 = ħωₖ/2
Derived, and the check on the first card: the field energy returns Unit 9's ½ħω per oscillator. Summed to a cutoff K, the energy density is ħcK⁴/(8π²).
⟨E²⟩ₖₛ = ħωₖ/(2ε₀V) is one mode's share; B in T, μ₀ in H m⁻¹; c²⟨B²⟩ₖₛ = ⟨E²⟩ₖₛ.Variance of the field averaged over a region
⟨Dz⟩ = ħc/(π²ε₀L⁴)
Derived: f damps each mode by e(−k²L²/4), so the integral converges with no cutoff at all. The infinity belongs to the field at a point, not to anything an instrument reads.
Ē = ∫f E d³x, Gaussian weight f = (πL²)(−3/2) e(−r²/L²), L in m; L = 1.00 µm gives Ēᵣₘₛ = 1.90 × 10⁴ V m⁻¹.Absorption reads only the normally ordered part
⟨0|E⁽⁻⁾⋅E⁽⁺⁾|0⟩ = 0
Derived from aₖₛ|0⟩ = 0. An absorbing detector's rate is proportional to ⟨E⁽⁻⁾⋅E⁽⁺⁾⟩ (Glauber, quoted), so it is dark in the vacuum, where ⟨E²⟩ is all ⟨E⁽⁺⁾⋅E⁽⁻⁾⟩.
E⁽⁺⁾ is the part of E carrying the aₖₛ, E⁽⁻⁾ = (E⁽⁺⁾)† carries the aₖₛ†, and E = E⁽⁺⁾ + E⁽⁻⁾.Interactive lab
Predict first. Explore the model. Compare the result.
Explore the model
Try this. Hold L at 1.0 µm and push n from 5 to 8: the dashed point value climbs a factor of 100 per decade of K, while the solid curve bends near K ≈ 1.7/L and stops. Then shrink L to 0.1 µm: the plateau rises as 1/L², the bend moves right, and the two curves stay together longer.
Live results
- POINT FIELD r.m.s. TO K
- 67248 V m⁻¹
- AVERAGED r.m.s. TO K
- 18632 V m⁻¹
- AVERAGED r.m.s., ALL MODES
- 19021 V m⁻¹
- SHORTEST WAVELENGTH 2π/K
- 1.987 µm
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 K = 10ⁿ m⁻¹: exponent n: 6.5; Averaging radius L: 1.0 µm
This experiment is kept only while this page is open. It does not change saved course progress.
Catch the common trap
The tempting reading of ⟨E²⟩ = ħcK⁴/(8π²ε₀) is that quantum theory predicts an infinite field everywhere, so the cutoff must be a real feature of space, perhaps the Planck length, and every prediction must depend on it. It does not. The divergent quantity is the square of the field at a mathematical point, and E(x) is an operator-valued distribution: like a δ-function it has meaning only when integrated against a smooth weight, and no instrument does anything else. Average the field with a smooth weight over a region of size L, or over a time window of length L/c, and the variance is finite, of order ħc/(ε₀L⁴) and exactly ħc/(π²ε₀L⁴) for a Gaussian weight, and the same for every cutoff beyond a few times 1/L. The cutoff stands for the probe, not for the vacuum. What stays unexplained is the absolute zero-point energy density, which only gravity could sense, and that is a question about gravity, not about any measurement in this lesson.
Practice & worked examples
3 authored examples. Make an attempt, then reveal the reasoning one step at a time.
EasyExample 1
Keep only the vacuum modes with wavelengths between 400 nm and 700 nm. Find the variance ⟨0|E⋅E|0⟩ they give at a point, the r.m.s. field and their zero-point energy density. Compare that energy density with sunlight's at Earth, 1.0 kW m⁻² carried at c, and with what thermal photons add in a room at 300 K.
Sketch your answer or reasoning, then compare it with the authored solution. Your writing is not automatically graded or saved.
0 of 5 steps revealed.
- Integrate the spectral density across the band: ⟨E²⟩ = ∫ ħω³/(2π²ε₀c³) dω = (ħ/(8π²ε₀c³))(ω₂⁴ − ω₁⁴), with ω = 2πc/λ: ω₁ = 2.691 × 10¹⁵ s⁻¹ at 700 nm and ω₂ = 4.709 × 10¹⁵ s⁻¹ at 400 nm.
- ħ/(8π²ε₀c³) = 5.599 × 10⁻⁵¹ V² m⁻² s⁴ and ω₂⁴ − ω₁⁴ = 4.393 × 10⁶² s⁻⁴, so ⟨E²⟩ = 2.46 × 10¹² V² m⁻² and Eᵣₘₛ = 1.57 × 10⁶ V m⁻¹. As ω⁴ predicts, the 400 nm to 550 nm half of the band supplies 81% of it.
- In the vacuum c²⟨B²⟩ = ⟨E²⟩, so u = (ε₀/2)⟨E²⟩ + ⟨B²⟩/(2μ₀) = ε₀⟨E²⟩ = 21.8 J m⁻³. Sunlight carries I/c = 1.0 × 10³ W m⁻² ÷ 3.00 × 10⁸ m s⁻¹ = 3.3 × 10⁻⁶ J m⁻³, so the zero-point band holds 6.5 × 10⁶ times more.
- At 300 K each mode adds n̄ = 1/(eˣ − 1) quanta to its ½, with x = ħω/(kB T). At 700 nm, x = 1.771 eV ÷ 0.02585 eV = 68.51, so the thermal share is 2n̄ ≈ 2e(−68.51) = 3.5 × 10⁻³⁰, and less at shorter wavelengths: visible-band fluctuations in a warm room are zero-point to better than one part in 10²⁹.
- Checks: ħω⁴/(ε₀c³) has units (J s)(s⁻⁴)/(F m⁻¹ × m³ s⁻³) = J/(F m²) = V² m⁻². As the band closes, ω₁ → ω₂, the variance goes to zero; as ω₂ → ∞ it grows as ω₂⁴, the divergence of the field at a point.
Answer. ⟨E²⟩ = 2.46 × 10¹² V² m⁻², Eᵣₘₛ = 1.57 × 10⁶ V m⁻¹ and u = 21.8 J m⁻³: 6.5 × 10⁶ times sunlight's 3.3 × 10⁻⁶ J m⁻³, while thermal photons at 300 K add at most 3.5 parts in 10³⁰, at the 700 nm edge, and less across the rest of the band.
MediumExample 2
A probe registers the transverse vacuum field averaged with the Gaussian weight f = (πL²)(−3/2) e(−r²/L²), with L = 2.5 µm. Find the r.m.s. field it registers and the wavelength at which the averaged spectrum peaks. Then find the mode cutoff K that captures 99% of the averaged variance, and compare with a point calculation cut off at the same K.
Sketch your answer or reasoning, then compare it with the authored solution. Your writing is not automatically graded or saved.
0 of 5 steps revealed.
- The weight's Fourier transform is e(−k²L²/4), so each mode's ħωₖ/(2ε₀V) is multiplied by e(−k²L²/2). With two polarisations and V∫d³k/(2π)³ over the modes, ⟨Dz⟩ is ħc/(2π²ε₀) times the integral of k³ e(−k²L²/2) over k from 0 to ∞.
- The integral of k³ e(−ak²) from 0 to ∞ is 1/(2a²), which for a = L²/2 is 2/L⁴. So ⟨Dz⟩ = ħc/(π²ε₀L⁴) = 3.618 × 10⁻¹⁶ V² m² ÷ (2.5 × 10⁻⁶ m)⁴ = 9.26 × 10⁶ V² m⁻², and Ēᵣₘₛ = 3.04 × 10³ V m⁻¹. Dimensions: V² m² over m⁴ is V² m⁻².
- The integrand k³e(−k²L²/2) peaks where 3/k = kL², at k = √3/L = 6.93 × 10⁵ m⁻¹, a wavelength 2πL/√3 = 9.07 µm: this probe mostly senses mid-infrared vacuum modes.
- With a sharp cutoff K as well, the integral runs to K and gives the fraction 1 − (1 + K²L²/2)e(−K²L²/2): 59.4% at K = 2/L. Reaching 0.99 needs K²L²/2 = 6.64, so KL = 3.64 and K = 1.46 × 10⁶ m⁻¹, a shortest wavelength of 4.31 µm.
- A point calculation with the same cutoff gives √(ħc/(8π²ε₀)) K² = 6.725 × 10⁻⁹ V m × (1.46 × 10⁶ m⁻¹)² = 1.43 × 10⁴ V m⁻¹, 4.7 times larger, and still rising as K² while the averaged value has stopped. Limits: L → 0 restores the divergence, and L → ∞ sends Ēᵣₘₛ to zero.
Answer. Ēᵣₘₛ = 3.04 × 10³ V m⁻¹; the averaged spectrum peaks at 9.07 µm; 99% of the variance is in modes with K up to 1.46 × 10⁶ m⁻¹ (4.31 µm), where a point calculation already gives 1.43 × 10⁴ V m⁻¹ and is still rising.
HardExample 3
Multiply each mode's contribution by e(−ηk), with η a length, and find the equal-time vacuum correlation ⟨0|E(x)⋅E(y)|0⟩ as a function of r = |x − y|. For η = 10 nm evaluate it at r = 0 and at r = 1.00 µm, find where it changes sign, and say what the result shows about where the divergence of ⟨E²⟩ lives.
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.
- With Σₛ êₖₛ⋅êₖₛ* = 2 and the box sum replaced by V∫d³k/(2π)³, the steps that gave ⟨E²⟩ now give (ħc/ε₀)∫d³k/(2π)³ k e(ik⋅r) e(−ηk). The angular average of e(ik⋅r) is sin(kr)/(kr), leaving ħc/(2π²ε₀r) times the integral of k² sin(kr) e(−ηk) over k from 0 to ∞.
- That integral is the imaginary part of the one with e(−ηk + irk) in place of sin(kr) e(−ηk), which is 2/(η − ir)³ = 2(η + ir)³/(η² + r²)³; its imaginary part is 2(3η²r − r³)/(η² + r²)³. So ⟨E(x)⋅E(y)⟩ = (ħc/(π²ε₀))(3η² − r²)/(η² + r²)³: derived. Dimensions: (V² m²)(m²)/(m⁶) = V² m⁻².
- At r = 0 it is 3ħc/(π²ε₀η⁴) = 3 × 3.618 × 10⁻¹⁶ V² m² ÷ (1.0 × 10⁻⁸ m)⁴ = 1.09 × 10¹⁷ V² m⁻², an r.m.s. field of 3.3 × 10⁸ V m⁻¹, growing as η⁻⁴ as the regulator is removed: the point divergence again.
- At r = 1.00 µm ≫ η it is −3.616 × 10⁸ V² m⁻², within 0.06% of the regulator-free −ħc/(π²ε₀r⁴) = −3.618 × 10⁸ V² m⁻²; the difference, a fraction 6η²/r², vanishes as η → 0. At any nonzero separation the correlation is finite and, once η ≪ r, independent of the regulator.
- It changes sign at r = √3 η = 17.3 nm: a positive core of height ∝ η⁻⁴ and width ∝ η, then the negative tail −ħc/(π²ε₀r⁴). Integrated over all space the two cancel exactly, since that integral is (ħc/ε₀) times k e(−ηk) at k = 0, which is zero: a uniform average over all space does not fluctuate.
- So the divergence lives entirely at coincident points. An average over L ≫ η integrates the core and tail against a smooth weight, and the result, ħc/(π²ε₀L⁴), no longer depends on η.
Answer. ⟨E(x)⋅E(y)⟩ = (ħc/(π²ε₀))(3η² − r²)/(η² + r²)³: 1.09 × 10¹⁷ V² m⁻² at r = 0 (3.3 × 10⁸ V m⁻¹ r.m.s.), −3.62 × 10⁸ V² m⁻² at 1.00 µm, changing sign at 17.3 nm. The divergence sits only at r → 0.
Quick check
Select an answer and explain it before checking.
Review & continue
Vacuum fluctuations of the electric field are its nonzero variance ⟨0|E⋅E|0⟩ in the photon vacuum, where the mean ⟨0|E|0⟩ is zero: each transverse mode contributes ħω/(2ε₀V), so the variance at a point diverges while any smooth average over a region is finite.
- 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.