Electromagnetism

The force whose sign a frequency chooses

A neutral particle in a non-uniform field is pulled towards strong field if it polarises more than the liquid around it and pushed away if it polarises less. With conduction in the picture, which of the two it does is decided by how fast the field alternates — and a living cell, a conductor wrapped in an insulator five nanometres thick, changes sign twice, where a dead one changes once. The same electrodes can then send the living cells one way and the dead ones the other.

Assumes: The force that lives where the model is not · The constant that depends on how fast it is asked

A drop of water with cells suspended in it, a pair of electrodes a few tens of micrometres apart, and an alternating voltage of a few volts. At one frequency of the voltage the cells crowd onto the edges of the electrodes. Turn the frequency down by a factor of ten and the same cells lift off and gather in the weak field above them. Nothing in the liquid has changed, the cells carry no net charge, and the force is not the one a battery exerts on an ion. It is the force on an induced dipole in a field that varies from place to place, and its sign is a property of the frequency.

A living cell and a dead one, told apart by a frequency. The real part of the dielectrophoretic factor for a living cell and a dead one, ten micrometres across, in a liquid of conductivity 0.01 S/m, against frequency. The cell is a conducting interior inside a membrane five nanometres thick. The living cell's intact membrane blocks slow fields, so at low frequency it behaves as an insulator and is pushed out of strong field; above 42 kHz the membrane is short-circuited by its own capacitance and the conducting interior is felt, so the cell is pulled in; above 140 MHz the permittivities take over and it is pushed again. The dead cell's membrane leaks, its interior has lost ions, and it changes sign only at 8.8 MHz. In 2 bands — below 42 kHz and from 8.9 MHz to 140 MHz — the same field pushes the two in opposite directions.
Fig. 1 The factor that sets the sign of the force on a living cell (solid) and a dead one (dashed), ten micrometres across, in a liquid of conductivity 0.01 S/m. Above zero a cell is pulled towards strong field; below zero it is pushed out of it. The living cell changes sign at 42 kHz and again at 140 MHz; the dead one only at 8.8 MHz. In the shaded bands the same field pushes the two in opposite directions.

A dipole in a gradient, and what it is compared with

The mechanism was set out in the argument about a glass slab drawn into a capacitor: a field polarises a neutral body, the induced dipole lines up with the field, and a dipole in a field that is stronger at one end than the other feels a net force towards the stronger end. The force goes as the gradient of the field’s square, so it does not care about the sign of the field. An alternating field works as well as a steady one, which is what makes the effect usable — a steady field in water would drive currents, electrolyse the electrodes and move every ion in the liquid.

For a sphere of radius aa in a liquid, the time-averaged force is

F=2πεma3ReK(ω)Erms2,K=εpεmεp+2εm,\mathbf{F} = 2\pi\varepsilon_m a^3\,\mathrm{Re}\,K(\omega)\,\nabla\left|\mathbf{E}_{\text{rms}}\right|^2, \qquad K = \frac{\varepsilon_p^* - \varepsilon_m^*}{\varepsilon_p^* + 2\varepsilon_m^*},

and nearly everything that is interesting is in KK. Without the stars it is the factor that decides how much field a sphere of dielectric keeps inside itself, and the numerator makes the essential point: the particle is not compared with empty space but with the liquid it displaces. A particle that polarises more than the liquid is pulled up the gradient; one that polarises less is pushed down it, for the same reason a bubble rises in water.

The stars carry the frequency. A material with conductivity σ\sigma in a field alternating at angular frequency ω\omega behaves exactly like a dielectric whose permittivity has an imaginary part, ε=εiσ/ω\varepsilon^* = \varepsilon - i\sigma/\omegaa conductor is a dielectric that is losing, and which description applies is a question of frequency. At low frequency the conductivity term dominates both starred permittivities and KK becomes (σpσm)/(σp+2σm)(\sigma_p - \sigma_m)/(\sigma_p + 2\sigma_m): the comparison is between how well the particle and the liquid conduct. At high frequency the conductivity terms have shrunk away and KK becomes (εpεm)/(εp+2εm)(\varepsilon_p - \varepsilon_m)/(\varepsilon_p + 2\varepsilon_m): the comparison is between how well they polarise. Nothing requires the two comparisons to agree.

Which way a neutral bead is pushed, against the frequency of the field. The real part of the factor K = (εp − εm)/(εp + 2εm) for three homogeneous beads in water of conductivity 0.01 S/m, against the frequency of the applied field. Positive K pulls a bead towards strong field and negative K pushes it away. At low frequency the factor is set by the conductivities alone and at high frequency by the permittivities alone, and one relaxation joins the two. A 200 nm latex bead runs from 0.25 to −0.48, changing sign at 3.3 MHz. A 2 µm latex bead runs from −0.36 to −0.48, never changing sign. A ceramic bead runs from −0.43 to 0.24, changing sign at 1.7 MHz. The bounds are +1 and −½ for any sphere whatever: a particle can be pulled twice as hard as it can be pushed.
Fig. 2 The sign-setting factor for three homogeneous beads in water of conductivity 0.01 S/m. A 200 nm latex bead, whose thin skin of counter-ions makes it conduct better than the water, is pulled at low frequency and pushed above 3.3 MHz. A 2 µm bead of the same latex conducts ten times less and is pushed at every frequency. A ceramic bead with a permittivity of 150 does the reverse of the small latex bead. The dotted lines at +1 and −½ bound every sphere.

Charge that cannot keep up

Why the two regimes are joined by a single step, and where the step falls, is a statement about how quickly charge can move.

In a slowly alternating field, free charge in the particle and in the liquid has time to flow and pile up at the particle’s surface. The accumulated surface charge makes the particle’s dipole whatever the conductivities dictate, and the permittivities are irrelevant. As the frequency rises, a point comes at which the field reverses before the charge has finished accumulating; above it the surface charge cannot form, and the dipole is whatever the bound charge in the two materials can manage. The crossover between the two is a relaxation — the Maxwell–Wagner relaxation, named for the two who worked out the polarisation of interfaces in layered dielectrics more than a century ago — and its characteristic frequency is the inverse of the charge-relaxation time of the combination, (σp+2σm)/2π(εp+2εm)(\sigma_p + 2\sigma_m)/2\pi(\varepsilon_p + 2\varepsilon_m).

The two latex beads in the second figure show where the conductivity of an insulator comes from, and it is a surprising place. Polystyrene conducts essentially nothing. What conducts is the thin cloud of counter-ions that every charged surface in water gathers around itself, a few nanometres thick, and its effect on the bead is that of a bulk conductivity equal to twice the surface conductance divided by the radius. The same surface on a smaller bead is spread over less volume, so the smaller bead behaves as the better conductor — ten times better for a bead ten times smaller. In water of 0.01 S/m the 200 nm bead out-conducts the liquid and is pulled at low frequency, while the 2 µm bead of the identical material does not and is pushed. A mixture of the two sizes separates in a field at a few hundred kilohertz, which is a way of sorting particles by size with no filter at all, using nothing but the fact that a surface and a volume scale differently.

A relaxation has two faces, and the second one turns out to be a separate force.

The two halves of the induced dipole: one pulls, one turns. The real and imaginary parts of the factor K for a bead of relative permittivity 2.55 and conductivity 0.02 S/m in water of 0.01 S/m. The real part is the share of the induced dipole in step with the field, and it sets the pull or push along a gradient. The imaginary part is the share a quarter of a cycle behind, and it sets the torque on the bead in a field that rotates — electrorotation. The imaginary part peaks where the real part changes most steeply, at 4.5 MHz, which is the Maxwell–Wagner frequency (σp + 2σm)/2π(εp + 2εm) at which charge can just keep up with the field at the bead's surface. The two curves are one analytic function seen from two sides.
Fig. 3 The real and imaginary parts of KK for the 200 nm latex bead. The real part is the share of the induced dipole in step with the field; it sets the pull or push along a gradient. The imaginary part is the share a quarter of a cycle behind; it sets the torque on the bead in a field whose direction rotates. The imaginary part peaks at 4.5 MHz, exactly where the real part falls most steeply, which is the Maxwell–Wagner frequency of the pair.

The imaginary part of KK is the part of the induced dipole that lags the field. In a field that only alternates in strength, the lagging part averages to no force at all. In a field that rotates — made with four electrodes driven a quarter of a cycle apart — the lagging dipole always points a little behind the field, and the field drags it round. The particle spins, at a rate proportional to ImK\mathrm{Im}\,K. This is electrorotation, and it is used to measure the same properties the translational force depends on, by counting revolutions under a microscope instead of watching drift.

The two parts are not independent measurements. They are the real and imaginary parts of one function of frequency that cannot respond before it is driven, and a function with that property has its two halves tied together: the step in one and the peak in the other are the same relaxation seen from two sides, and either can be computed from the other.

Why the bounds are +1 and −½

The factor KK has limits that no choice of materials can break, and they are worth seeing because they are lopsided. If the particle is enormously more polarisable or conductive than the liquid, K+1K \to +1. If it is enormously less, K1/2K \to -1/2. The asymmetry comes from the 2εm2\varepsilon_m in the denominator, which in turn comes from the geometry of the field round a sphere: the field that leaks around an insulating sphere and the field that is sucked into a conducting one are not mirror images.

The practical consequence is that the most a particle can be pulled is twice the most it can be pushed. A liquid can be made to repel almost anything by making it conduct far more than the particles — but not strongly, and never more strongly than a well-chosen particle can be attracted.

A cell is a conductor wrapped in an insulator

A living cell is not a homogeneous bead. Its interior is a salty solution that conducts nearly as well as seawater, and it is enclosed by a membrane of lipid about five nanometres thick that conducts almost not at all. A membrane that thin has an enormous capacitance per unit area, about a microfarad per square centimetre, and that capacitance is the whole story of the living cell’s curve in the first figure.

At low frequency the membrane blocks current, so the field cannot reach the conducting interior and the cell behaves as an insulating sphere. In a liquid that conducts at all, an insulating sphere has KK near 1/2-1/2, and the cell is pushed out of strong field. As the frequency rises, the membrane’s capacitance lets current through — a capacitor’s impedance falls as the frequency rises — and above a certain frequency the membrane is effectively short-circuited. The cell now presents its conducting interior to the field, and if the interior conducts better than the liquid, KK swings positive and the cell is pulled in. At still higher frequency the conductivities stop mattering altogether, the permittivities take over, and since the cytoplasm polarises a little less than water, KK ends slightly negative again. Two crossovers, and three regimes.

The first crossover is set by the membrane’s capacitance and the liquid’s conductivity together. For a cell of radius RR whose membrane has capacitance CmC_m per unit area, in a liquid of conductivity σm\sigma_m, it falls close to

f12σm2πRCm,f_1 \approx \frac{\sqrt2\,\sigma_m}{2\pi R\,C_m},

which for a membrane five nanometres thick with a relative permittivity of six — a microfarad per square centimetre — a radius of five micrometres and a liquid of 0.01 S/m comes to about 42 kHz, where the exact curve crosses. The formula says two things the curve alone does not. The crossover moves in proportion to the liquid’s conductivity, so raising it tenfold raises the crossover tenfold; and it moves inversely with the cell’s size, so large cells cross over sooner than small ones in the same liquid.

A living cell and a dead one, told apart by a frequency. The real part of the dielectrophoretic factor for a living cell and a dead one, ten micrometres across, in a liquid of conductivity 0.1 S/m, against frequency. The cell is a conducting interior inside a membrane five nanometres thick. The living cell's intact membrane blocks slow fields, so at low frequency it behaves as an insulator and is pushed out of strong field; above 400 kHz the membrane is short-circuited by its own capacitance and the conducting interior is felt, so the cell is pulled in; above 150 MHz the permittivities take over and it is pushed again. The dead cell's membrane leaks, its interior has lost ions, and it never changes sign. In one band — from 410 kHz to 150 MHz — the same field pushes the two in opposite directions.
Fig. 4 The same living and dead cells in a liquid ten times more conductive, 0.1 S/m. The living cell’s first crossover has moved from 42 kHz to 400 kHz, and its low-frequency push is as strong as before. The dead cell, whose interior now conducts less than the liquid, is pushed at every frequency; the band in which the two are driven apart runs from 410 kHz to 150 MHz.

A cell that has died has a membrane that leaks. Its conductivity rises by orders of magnitude, the ions inside escape towards the concentration outside, and the capacitive barrier that produced the living cell’s low-frequency push is gone. The dead cell in the first figure is pulled at low frequency, where the living one is pushed, and has only one crossover. Anywhere in the shaded bands, one field applied through one pair of electrodes sends the two kinds in opposite directions — which is how living and dead yeast, healthy and infected blood cells, and cancer cells from the blood cells around them have all been sorted, without labels and without touching them.

The crossover frequency is itself a measurement. Since it moves with the membrane’s capacitance per unit area, and that capacitance depends on how much the membrane is folded, a cell with a rougher, more folded surface crosses over at a lower frequency. Measuring the frequency at which a cell stops moving measures the area of membrane it carries per unit of apparent surface.

Pulled particles land on metal; pushed ones can float

The direction of the force has a consequence that is not obvious from the formula, and it connects the whole subject to a theorem about static fields.

The same field, and the two directions a neutral particle can be pushed in it. Two electrode strips on a floor held at opposite voltages, seen end on, with contours of the field's strength squared — closest together at the facing edges of the strips, where the field is strongest. The arrows are the direction of ∇|E|², which is the force on a particle that polarises more than the liquid around it (left), and the reverse for one that polarises less (right). The same electrodes collect the first kind at the gap between the strips and lift the second kind away into the weak field above them. Only the directions are drawn; the strips are modelled as two line charges just below the floor, which puts the field's maxima and minima in the right places and says nothing about their exact strength.
Fig. 5 Two electrode strips on a floor at opposite voltages, seen end on, with contours of the field’s strength. On the left, the direction of the force on a particle that polarises more than the liquid: towards the facing edges of the strips, where the field is strongest. On the right, the force on one that polarises less: up and away into the weak field above. Only directions are drawn; the strips are modelled as two line charges below the floor.

A particle pulled up the gradient goes where the field is strongest, and the field is strongest at the surface of the electrodes, at their edges. It cannot be held anywhere in between, and the reason is exact rather than practical. Each component of a static field in charge-free space satisfies Laplace’s equation, so the square of the field’s strength can have no maximum in the interior of the space: it is largest somewhere on the boundary. That is the same fact that makes it impossible for any arrangement of static charges to hold a charge still, transplanted from the potential to the field strength. A positively polarising particle always ends on an electrode.

A negatively polarising one is in the opposite situation. The square of the field’s strength can have a minimum in free space — it can be zero at a point surrounded by stronger field on every side — and a particle pushed down the gradient sits in that minimum stably. Eight electrodes at the corners of a small cube, driven in the right pattern, make a field with a zero at the centre, and a cell pushed out of strong field floats there in a cage with no walls. The alternating field and the time average are what make the dipole always align, but the reason the cage is possible at all is that the forbidden extremum is the maximum and not the minimum. It is the same asymmetry that lets a rapidly alternating force hold something up where no steady one could.

The grip goes as the volume

The force is proportional to a3a^3, which means dielectrophoresis is a tool for things of a particular size and becomes useless below it.

How deep a trap a field can dig, against the size of what it holds. The depth of the dielectrophoretic trap, 2π εm a³ |Re K| E², in units of the thermal energy kT at room temperature, against the particle's radius on logarithmic axes, for |Re K| = 0.5 and three field strengths. The depth goes as the cube of the radius, so a factor of ten in size is a factor of a thousand in grip. At 10⁵ V/m a trap is one kT deep for a radius of 57 nm; At 10⁶ V/m a trap is one kT deep for a radius of 12 nm; At 10⁷ V/m a trap is one kT deep for a radius of 3 nm. Below that radius the particle is shaken out of any trap the field can make — which is why a cell is held by a few volts across a gap of tens of micrometres, while a protein needs something like ten million volts per metre, and heats the liquid doing it.
Fig. 6 The depth of the trap a field can make — the energy 2πεma3ReKE22\pi\varepsilon_m a^3|\mathrm{Re}\,K|E^2 — in units of the thermal energy kTkT, against the particle’s radius, for three field strengths. The slope is three: ten times smaller is a thousand times weaker. At 10610^6 V/m a trap is one kTkT deep for a radius of 12 nm; below that, the particle’s own thermal motion carries it out of any trap the field can dig.

What any trap has to beat is the jiggling of thermal motion, which gives every particle an energy of order kTkT whatever its size. A cell ten micrometres across held by a field of 10510^5 V/m — a few volts across a few tens of micrometres — sits in a trap millions of kTkT deep and does not move at all except as the field directs. A virus fifty nanometres across needs fields ten times stronger to be held firmly. A protein a few nanometres across needs 10710^7 volts per metre, and at that field the current flowing through any real buffer heats the liquid, and the heating drives flows that push the protein about far harder than the dielectrophoretic force does. The size below which the method stops being worth using is set by that competition rather than by any limit on the field.

The same a3a^3 explains the optical version of the force. At the frequency of light the factor KK is set by the refractive indices, the gradient is the focus of a laser beam, and a microscope objective can make a gradient on a scale of a wavelength — which is the tweezer that holds a bead in a beam, dielectrophoresis at 101410^{14} Hz.

Where the single-shell cell stops being a cell

The field is taken to be uniform across the particle. The formula treats the particle as a point dipole in the local gradient, which requires the field to change little across it. Near the edge of an electrode, where the field changes on a scale of micrometres, a cell ten micrometres across violates that, and higher multipoles of its induced charge contribute forces the dipole term does not.

The cell is one shell and one interior. Real cells have a nucleus with its own envelope, organelles and a cell wall in plants and yeast, and each adds a relaxation to the factor. The single-shell model captures the first crossover well and the high-frequency behaviour only roughly.

Other forces are ignored. The field heats the liquid in proportion to its conductivity, and the heating makes flows; the field also acts on the charge in the thin layer at the electrode surfaces, making flows of its own at low frequency. Both can move particles faster than dielectrophoresis does, and much of the craft of using the effect lies in choosing the frequency and conductivity at which they are weak.

And the particles do not interact. A dense suspension of polarised particles attracts itself: neighbouring dipoles along the field attract, and particles form chains along the field lines — pearl chains — whose behaviour belongs to the physics of the suspension rather than of any one particle.

A spread of cells, not one crossover

Every curve here is a property of a single particle computed from bulk numbers — conductivities, permittivities and a membrane thickness — each of which varies from cell to cell by tens of per cent. A real population produces a spread of crossover frequencies, and the shaded bands of the first figure are narrower in practice than they look, because the edges of two distributions overlap. Whether a sorting succeeds is a question about those distributions, and a single curve cannot answer it.

Nor can the field picture show time. A particle in the right-hand panel does not arrive at the field’s minimum instantly; it drifts there against viscous drag at a speed proportional to a2a^2, so small particles are both weakly held and slow to arrive. The directions in the figure say where things go, not how long they take to get there.

Still open: whether the membrane’s electrical signature can diagnose a cell

That a cell’s crossover frequency shifts when it dies, when it becomes cancerous, when it is infected by a parasite, or when it differentiates into another type is well established; each changes the membrane’s capacitance per unit area or the interior’s conductivity. What is not established is whether those shifts are specific enough to be used as a diagnosis rather than as a preliminary sort. The electrical properties of a cell are a handful of numbers, and a great many different changes to a cell move the same handful. Studies that separate one cancer cell line from blood in a laboratory succeed; whether a crossover measurement on a patient’s cells can distinguish a disease from the ordinary variability between people is still being argued, and it depends on how much the few numbers the membrane can report really say about the cell behind it.

The next question the subject raises is what happens when the particles are no longer far apart. A suspension of polarisable particles in an alternating field stiffens from a liquid into a paste that will hold a load within milliseconds of the field being switched on, which is the basis of electrorheological fluids and the reason the physics of one polarised particle is only the beginning of the physics of many. The habit worth carrying from here is to ask of any polarisation what it is being compared with. The force on a neutral particle measures the difference between the particle and what it displaces, at the frequency being asked — change the liquid or change the frequency, and the same particle is pulled or pushed.

Part 5 of 5

This essay is one argument about Dielectrics. 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.

Brownian motionConductivityDielectricDipoleDispersionHarmonic functionPermittivityPolarisabilityRelaxation