Electromagnetism

The static field a metal may or may not reflect

The Casimir force between two gold plates at room temperature has been computed two ways for twenty-five years, and the two answers differ by forty per cent at three micrometres. Both use the same measured optical properties of gold. They part company over a single number: whether gold reflects a magnetic field that does not change in time. Common sense and every electrical engineer say it does not, since a static field passes straight through a copper sheet. The most precise force measurements say it does. The disagreement is not about the Casimir effect at all; it is about what a metal's electrons do with a fluctuation slower than their collisions.

Assumes: The force that forgets Planck's constant · The attraction that weakens when light is too slow

The force that forgets Planck’s constant followed the Casimir pressure between two plates out from the quantum regime, where it falls as the inverse fourth power of the gap, to the thermal one, where it falls as the inverse cube and contains Planck’s constant not at all. Far out, only the static fluctuations of the field survive, and the pressure is a purely classical ζ(3)kT/4πa3\zeta(3)kT/4\pi a^3 — if both polarisations of those static fluctuations are reflected. For plates of real metal, that essay found, one of them may not be, and the force is then half as large. It drew the two possibilities with perfect mirrors and left the question of which is right.

This essay puts real gold in the gap. Gold’s response to light is measured, from the infrared to the ultraviolet, to a few per cent, and it is the input to the calculation Evgeny Lifshitz devised in 1955 for the force between any two materials. With real gold the question of the static term does not go away. It becomes sharper, and it turns out to be a question about one point on a curve that nobody can measure.

Two descriptions of gold that agree everywhere but at zero. How strongly gold answers a slowly varying field, (ε(iξ) − 1)ξ²/ωₚ² in the two forms of its permittivity at imaginary frequency ξ, against ξ on a logarithmic axis: the plasma model (dashed), which ignores the electrons' collisions and gives exactly 1, and the Drude model (solid), which includes them at the measured rate ħγ = 35 meV and falls to zero below about γ = 5.3·10¹³ per second. The Lifshitz force samples this curve only at the Matsubara frequencies 2πlkT/ħ (dots, 300 K): the first, 2.47·10¹⁴ per second, lies above γ, and there the two descriptions differ by 18 per cent, at the next by 10, and less beyond. At l = 0, ξ = 0, they part completely: 1 against 0. Everything the two models disagree about is in that one point.
Fig. 1 Gold’s response to a slowly varying field, (ε(iξ)−1)ξ2/ωp2(\varepsilon(i\xi) - 1)\xi^2/\omega_p^2, against imaginary frequency ξ\xi: the plasma description (dashed), which ignores the electrons’ collisions and gives exactly 1, and the Drude description (solid), with the measured collision rate ℏγ\hbar\gamma = 35 meV, falling to zero below γ=5.3×1013\gamma = 5.3 \times 10^{13} per second. The Lifshitz force samples it only at the Matsubara frequencies (dots, 300 K): at the first the two differ by 18 per cent, at the second by 10, and less beyond. At ξ=0\xi = 0 they part completely.

The force adds up a list of frequencies

The Lifshitz formula at a temperature TT is a sum. Each term belongs to one of the Matsubara frequencies, ξl=2πlkT/ℏ\xi_l = 2\pi l kT/\hbar, which are the frequencies at which a thermal field at temperature TT can be decomposed, and each term depends on how strongly the plates reflect fluctuations of that frequency. At room temperature the first nonzero frequency, ξ1\xi_1, is 2.5×10142.5 \times 10^{14} per second, an infrared frequency; the others are its multiples. The term with l=0l = 0 is the static field, and it is weighted by one half.

The frequencies are imaginary: the formula needs the permittivity ε(iξ)\varepsilon(i\xi) on the imaginary axis, which is obtained from the measured real-frequency absorption by the dispersion relation that the answer that cannot come first derived from causality. That step is standard and uncontroversial. The difficulty is that the measurements stop at low frequency — optical data for gold run down to a few tenths of an electronvolt — and the dispersion relation needs an extrapolation below that. For a metal the extrapolation is the free-electron response, and there are two forms of it.

The plasma description treats gold’s conduction electrons as a gas that accelerates freely in a field: ε(iξ)=1+ωp2/ξ2\varepsilon(i\xi) = 1 + \omega_p^2/\xi^2, with the plasma frequency ℏωp=9.0\hbar\omega_p = 9.0 eV fixed by the density of electrons. The Drude description adds the collisions that give a real metal its resistance: ε(iξ)=1+ωp2/ξ(ξ+γ)\varepsilon(i\xi) = 1 + \omega_p^2/\xi(\xi+\gamma), with a collision rate γ\gamma of about 5×10135 \times 10^{13} per second for gold at room temperature, the number that sets its electrical resistivity. The frequency below which nothing gets in used the plasma form for the ionosphere, where collisions are rare; a metal at room temperature is a plasma whose electrons collide every twenty femtoseconds.

The first figure plots the combination that enters the reflection of a transverse-electric fluctuation, (ε−1)ξ2(\varepsilon - 1)\xi^2, scaled by ωp2\omega_p^2. For the plasma description it is exactly one at every frequency. For the Drude description it follows the plasma line at high frequency and falls away below the collision rate. The Matsubara frequencies are marked, and they all lie above the fall: at the first the two descriptions differ by 18 per cent, at the second by 10, and at higher ones by less. Then there is the point at zero, which a logarithmic axis cannot show and the Lifshitz sum cannot omit. There the plasma description gives one and the Drude description gives zero.

A static field and a copper sheet

Whether a metal reflects a static magnetic fluctuation. The reflection coefficient of gold for a transverse-electric fluctuation at zero frequency — a static magnetic field varying across the surface with wavenumber q — against q in units of ωₚ/c, the inverse of gold's 22 nm plasma length. The plasma model (dashed) reflects it almost perfectly wherever the field varies slowly compared with that length: −0.993 at the wavenumber a 3 μm gap selects (dot). This is a superconductor's behaviour — field expelled by currents that never decay. The Drude model (solid) reflects nothing at all: in a metal whose currents are damped by collisions, a static magnetic field soaks through, as it does through a copper sheet. Every difference between the two predictions of the force comes from this zero-frequency number.
Fig. 2 The reflection coefficient of gold for a transverse-electric fluctuation at zero frequency — a static magnetic field varying across the surface with wavenumber qq — against qq in units of ωp/c\omega_p/c, the inverse of gold’s 22 nm plasma length. The plasma description (dashed) reflects it almost perfectly where the field varies slowly on that scale: −0.993 at the wavenumber a 3 μm gap selects (dot), a superconductor’s behaviour. The Drude description (solid) reflects nothing: the field soaks through.

The physical content of the zero-frequency point is easier to see than the formula. A transverse-electric fluctuation at zero frequency is a static magnetic field, varying in strength across the surface of the plate and unchanging in time. The second figure shows how strongly each description reflects it. The plasma description reflects it almost perfectly: for a field varying over the micrometre scale of the gap, the reflection coefficient is −0.993. The Drude description reflects it not at all.

Here the Drude description has the whole of classical electromagnetism behind it. A static magnetic field passes through a sheet of copper or gold as if the sheet were not there. When the field was switched on, it did induce currents in the metal that opposed it — the circuit that fights its own change is the rule — but those currents met resistance and died away, and with them the opposition. The magnet that falls slowly down a copper tube is slowed by exactly those currents, and a magnet held still inside the tube feels nothing. Only a superconductor, whose currents never decay, keeps a static field out, which is what the field that is pushed out made the defining property of superconductivity. The plasma description, with no collisions, describes a metal whose electrons are dissipationless and so gives an ordinary metal a superconductor’s response to static magnetic fields.

So the Drude description looks correct and the plasma description looks like an error. At the frequencies that matter for the force at short distances that error is harmless — the plasma line and the Drude line agree there — and the plasma form is often used as a convenient approximation. The trouble starts at the zero-frequency term, where the error, if it is one, is total.

Forty per cent at three micrometres

The force the two descriptions predict, against the gap. The Casimir pressure between two gold plates at 300 K computed with the Drude description, divided by the same pressure computed with the plasma description, against their separation, from the Lifshitz sum with gold's measured plasma frequency (9 eV) and collision rate (35 meV). At 0.1 μm the two agree to 2.1 per cent; at 1 μm the Drude pressure is 16 per cent lower; at 3 μm, 41 per cent; at 10 μm it is 0.50 of the plasma value, approaching one half, where only the static terms survive and the Drude metal has lost its TE half. For comparison, perfect mirrors with that one term removed (dashed) give nearly the same curve: the finite plasma frequency matters only at the shortest gaps.
Fig. 3 The Casimir pressure between gold plates at 300 K in the Drude description divided by that in the plasma description, against separation, from the Lifshitz sum with gold’s plasma frequency (9 eV) and collision rate (35 meV). At 0.1 μm they agree to 2.1 per cent; at 1 μm the Drude pressure is 16 per cent lower; at 3 μm, 41; at 10 μm it is 0.50 of the plasma value. Perfect mirrors with only the static TE term removed (dashed) give nearly the same curve.

The third figure computes the full Lifshitz pressure both ways, with gold’s measured plasma frequency and collision rate, and divides one by the other. At a tenth of a micrometre the two agree to two per cent. The gap is small compared with the thermal wavelength ℏc/kT\hbar c/kT, so many Matsubara frequencies contribute, and the static term is a small share of the total. As the gap widens, fewer frequencies matter and the static term grows in relative importance. At a micrometre the Drude pressure is sixteen per cent lower; at three, forty-one per cent; at ten, half. Far out, the pressure is the static term alone, and in the Drude description the static term has lost its TE half.

The dashed curve is the perfect-mirror model of the earlier essay, with only that one term removed, and it lies almost on top of the gold curve. The comparison makes the point precisely. Everything about real gold — its finite plasma frequency, its collisions at the nonzero Matsubara frequencies, the penetration of fields into its surface — changes the ratio by a few per cent at the shortest gaps and by nothing further out. The whole of the forty per cent at three micrometres is the static TE term.

Where the two predictions part, term by term. The Casimir pressure between gold plates 3 μm apart at 300 K, in micropascals, split into the parts the Lifshitz sum is made of: the static (zero-frequency) term for each polarisation and the sum over every other Matsubara frequency, for the plasma description (upper bar) and the Drude description (lower bar). The static TM terms are equal, 7.337 μPa, because both models make gold a perfect mirror for a static electric field. The sums over the other frequencies differ by 0.8 per cent. The static TE term is 7.025 μPa in the plasma description and zero in the Drude, and it accounts for 99.7 per cent of the 7.048 μPa between the two totals.
Fig. 4 The pressure between gold plates 3 μm apart at 300 K, in micropascals, split into its parts for the plasma description (upper bars) and the Drude description (lower): the static TM terms are equal, 7.337 μPa, since both make gold a perfect mirror for a static electric field; the sums over all other Matsubara frequencies differ by 0.8 per cent; the static TE term is 7.025 μPa in one and zero in the other, and makes up 99.7 per cent of the 7.048 μPa between the totals.

The fourth figure takes the three-micrometre pressure apart. The static TM term is the response to a static electric field, and both descriptions make gold a perfect mirror for that: the electrons rearrange until the field inside is zero, whatever their collisions, as the inside of a conductor showed. The sum over all the nonzero frequencies differs between the descriptions by less than one per cent. The static TE term is seven micropascals in the plasma description and nothing in the Drude. The disagreement between two calculations that share every measured input is one term, of one sign, at one frequency.

Why that term grows in importance with the gap has a simple reason. A fluctuation at Matsubara frequency ξl\xi_l contributes to the force only if it can reach across the gap, and its contribution falls as e−2ξla/ce^{-2\xi_l a/c}. At a tenth of a micrometre, 2ξ1a/c2\xi_1 a/c is 0.16 and dozens of frequencies contribute nearly equally, so the static term is one of many. At three micrometres the same exponent is 4.9: the first nonzero frequency contributes less than a hundredth of what it would at zero gap, and the higher ones nothing. The force has become a sum of one term and a correction, and whatever is uncertain about that one term is uncertain about the force. The crossover is set by the thermal wavelength ℏc/kT\hbar c/kT, 7.6 micrometres at room temperature: at gaps well below it the force is quantum and the dispute is a small correction; at gaps above it the force is thermal and the dispute is the whole answer.

The same fluctuations do other work, and there the dissipative picture is not in doubt. The heat that crosses a gap too narrow for light found that two surfaces closer than a thermal wavelength exchange heat far faster than blackbody radiation allows, through evanescent fields that tunnel across the gap. For two metal surfaces, calculations find that this heat is carried largely by the transverse-electric fluctuations — slowly varying magnetic fields that induce eddy currents in the far surface and are absorbed by its resistance. Those fluctuations exist only because the electrons collide; a collisionless plasma neither emits nor absorbs them. The metal that must dissipate to carry heat across a gap is, in the force measurements at short range, behaving as though its static magnetic response did not dissipate at all.

The measurements that disagree with each other

Where each kind of experiment can tell the descriptions apart. How much smaller the Drude description makes the Casimir pressure between gold surfaces than the plasma description, in per cent, against separation at 300 K, with the ranges explored by two kinds of experiment shaded. Between 0.16 and 0.75 μm, where micromechanical oscillators measure the force gradient to a fraction of a per cent, the difference is 3.1 to 11.9 per cent; those experiments agree with the plasma description. Between 0.7 and 7 μm, where a torsion pendulum measured the force in 2011, the difference grows from 11 to 49 per cent — large, but in a range where the electrostatic force of uneven surface potentials is comparable and had to be fitted and subtracted; that experiment agreed with the Drude description. The first sees a small difference measured cleanly; the second a large one measured beside a large background.
Fig. 5 How much lower the Drude description puts the pressure than the plasma description, in per cent, against separation at 300 K, with the ranges of two kinds of experiment shaded. Between 0.16 and 0.75 μm, where micromechanical oscillators measure the force gradient to a fraction of a per cent, the difference is 3.1 to 11.9 per cent, and those experiments agree with the plasma description. Between 0.7 and 7 μm, where a torsion pendulum measured the force in 2011, it is 11 to 49 per cent, beside an electrostatic background that had to be fitted; that experiment agreed with the Drude description.

The question would be academic if it could not be tested, and it has been tested repeatedly, with results that do not agree. The last figure places the experiments on the curve of the difference.

At the short end, Ricardo Decca and colleagues measured the force between a gold-coated sphere and a gold-coated plate oscillating on a micromachined torsional oscillator, from 0.16 to 0.75 micrometres, repeatedly between 2003 and 2016 with improving precision. In that range the two descriptions differ by three to twelve per cent, and the measurements, precise to a fraction of a per cent, sit on the plasma prediction and exclude the Drude one by many standard deviations. Other measurements since — with an atomic-force microscope and surfaces of nickel, a magnetic metal for which the two descriptions differ in a different way, and with a differential arrangement designed so that the two predictions for the measured quantity differ enormously — have found the plasma description again.

At the long end, Alex Sushkov, Woo-Joong Kim, Diego Dalvit and Steve Lamoreaux measured the force between a gold-coated lens and plate on a torsion pendulum, from 0.7 to 7 micrometres, and published in 2011 a result that agreed with the Drude description. Here the difference between the descriptions is large, tens of per cent, but so is a competing force. Real metal surfaces are not equipotentials: grains of different crystal orientation have slightly different work functions, and the patches of different potential produce an electrostatic attraction that, at a few micrometres, is comparable with the Casimir force itself. The experiment measured a total force many times larger than the predicted thermal Casimir force and extracted it by fitting and subtracting the patch contribution. Critics have argued that the subtraction carries more uncertainty than was assigned; the authors that the short-range experiments neglect the patches, which matter less there but not nothing.

The standoff is the reverse of what one would expect. The experiments that see a small difference cleanly favour the description that appears physically wrong. The experiment that sees a large difference, but beside a large background, favours the description that appears physically right.

The entropy at absolute zero

There is a theoretical argument against the Drude description as well, and it is about thermodynamics rather than electromagnetism. The Casimir energy between two plates is a free energy, and its derivative with respect to temperature is an entropy. For the Drude description with a collision rate that vanishes as the temperature falls — as it does in an ideal crystal with no impurities, where the only scattering is by lattice vibrations that freeze out — the Casimir entropy at absolute zero comes out not zero but a finite negative number that depends on the separation. That would violate the third law of thermodynamics, which the two laws absolute zero will not allow took apart, and it was pointed out by Vladimir Mostepanenko and colleagues in 2002 as evidence that the Drude description is being misapplied.

The answer offered by others is that real metals always have impurities, so the collision rate does not fall to zero but levels off at a residual value set by the defects, and with a residual rate the Drude entropy does go to zero, though only at temperatures so low that no experiment reaches them. Whether an argument about a perfect crystal at absolute zero can rule out a description of an imperfect crystal at room temperature is itself disputed. What the argument shows is that the static term is delicate at both ends: what a metal does with a fluctuation at zero frequency depends on processes slower than anything an optical measurement sees.

What the dispute is really about

The Lifshitz formula treats the plates as having a local permittivity, a number at each frequency that relates the current at a point to the field at that point. At zero frequency and over micrometre distances that description may simply be wrong. A static magnetic field varying over a micrometre across a gold surface is not a uniform field, and the electrons that respond to it travel tens of nanometres between collisions at room temperature. Whether their response is local — the assumption behind both the Drude and plasma forms — or whether the spatial variation matters, as it does in the anomalous skin effect at low temperature, has been examined with nonlocal permittivities that make the two pictures blend, and the conclusions depend on details of the surface that are hard to know.

There is also the question of what the electrons’ collisions are. In the Drude description a collision randomises an electron’s momentum and dissipates energy, and a field that varies more slowly than the collision rate drives a current limited by resistance. A static fluctuation of the thermal field is a fluctuation, not a field applied by an outside agent, and some theorists have argued that dissipative and dissipationless responses enter the fluctuation–dissipation relation differently at zero frequency, so that the static term should be computed in a way neither standard description does. Half a kT in a piece of wire found the fluctuation–dissipation relation fixing the noise in a resistor from its resistance alone; whether the same reasoning applies without modification to the zero-frequency fluctuations of the field between two resistors is the point at issue.

What a reflection coefficient leaves out

The figures show the pressure between two infinite, flat, smooth plates of gold in thermal equilibrium at 300 K, and every experiment differs from that in every respect. The surfaces are a sphere or a lens against a plate, handled by the proximity-force approximation, which is accurate to a fraction of a per cent at the gaps involved but is still an approximation. The surfaces are rough at the scale of tens of nanometres, the gold films are a few hundred nanometres thick and have their own grain structure, and their optical properties differ from tabulated bulk values by a few per cent. Each of these corrections is computed and applied, and the experimental groups have argued over each. None of them is drawn here, and none of them affects the main point the figures make: that the two descriptions share everything except the static TE term, so any measurement that separates them is a measurement of that one term.

Still open: how a superconductor would vote

The most direct test proposed so far would change the static TE term while holding everything else fixed. A plate of a superconducting metal, cooled through its transition temperature, changes from a metal that lets a static magnetic field through to one that expels it. In the Drude description the static TE reflection would jump from zero to nearly one at the transition, adding a measurable step to the force between a plate and a sphere at a few micrometres; in the plasma description, which already treats the metal as reflecting static fields, nothing should change. Platforms that measure the force between superconducting surfaces at a few kelvin have been built, and the step has been looked for, but none has yet reached the precision needed to see it or to rule it out.

The habit worth carrying away is to locate a disagreement exactly. Two calculations of the Casimir force between gold plates that share every measured property of gold differ by 16 per cent at a micrometre and 41 per cent at three — and the whole difference is the reflection of a static magnetic fluctuation, about −1 if the electrons never collide and exactly 0 if they do, at a frequency no optical measurement reaches. When a prediction rests on an extrapolation, the part of it that depends on the extrapolation is the part a measurement is testing.

Part 7 of 7

This essay is one argument about Retardation. The others:

The objects named here

The third axis, after the field and the reading path: the things themselves, and every essay that touches each one.

Casimir effectDrude modelEddy currentsLifshitz theoryMatsubara frequencyPlasma frequencyReflection coefficientThermal fluctuations