Skip to main content

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.
Starting model: Vacuum fluctuations: ⟨E²⟩ and its spectrumLog-log plot of r.m.s. vacuum field against mode cutoff K = 10ⁿ m⁻¹, n = 6.5. Dashed: the field at a point, 67248 V m⁻¹, rising as K² without limit. Solid: averaged over a Gaussian of radius L = 1.0 µm, 18632 V m⁻¹, levelling onto the faint line at 19021 V m⁻¹, its all-mode value.r.m.s. vacuum field / V m⁻¹, log scaledashed: E at a point, grows as K²solid: averaged over L = 1.0 µm10²10⁴10⁶10⁵10⁶10⁷10⁸mode cutoff K / m⁻¹, log scale
Start with a picture. Change the model in the lab below.
1

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.

2

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⁽⁻⁾.
3

Interactive lab

Predict first. Explore the model. Compare the result.

Explore the model

6.5
5 to 8; steps of 0.1.
1.0 µm
0.1 to 5 µm; steps of 0.1.

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.

Interactive Advanced Quantum Mechanics modelLog-log plot of r.m.s. vacuum field against mode cutoff K = 10ⁿ m⁻¹, n = 6.5. Dashed: the field at a point, 67248 V m⁻¹, rising as K² without limit. Solid: averaged over a Gaussian of radius L = 1.0 µm, 18632 V m⁻¹, levelling onto the faint line at 19021 V m⁻¹, its all-mode value.r.m.s. vacuum field / V m⁻¹, log scaledashed: E at a point, grows as K²solid: averaged over L = 1.0 µm10²10⁴10⁶10⁵10⁶10⁷10⁸mode cutoff K / m⁻¹, log scale
Log-log plot of r.m.s. vacuum field against mode cutoff K = 10ⁿ m⁻¹, n = 6.5. Dashed: the field at a point, 67248 V m⁻¹, rising as K² without limit. Solid: averaged over a Gaussian of radius L = 1.0 µm, 18632 V m⁻¹, levelling onto the faint line at 19021 V m⁻¹, its all-mode value.

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

Baseline 6.5; choose 5 to 8, in steps of 0.1.
4. Compared with the saved baseline, this result will…

This experiment is kept only while this page is open. It does not change saved course progress.

4

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.

5

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.

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.

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.

6

Quick check

Select an answer and explain it before checking.

In the vacuum, the modes with wavelengths between 400 nm and 700 nm give ⟨0|E⋅E|0⟩ ≈ 2.5 × 10¹² V² m⁻² at a point, an energy density millions of times that of sunlight. Why does a silicon photodiode, band gap 1.12 eV, sealed in a dark box at absolute zero, register no photocurrent from it?
7

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.

Extra challenging questions

Connect several ideas, explore supplied data, and compare your reasoning with a worked solution.