Skip to content

Ions Fading by Pairs

Irodov 3.215DifficultyDifficulty 2 of 5UniquenessUniqueness 3 of 5

Problem Irodov 3.215

An ionizer that makes n˙i=3.5×109\dot n_i = 3.5\times10^{9} ion pairs per cm3^3 of air per second has run long enough for the ion concentration to settle, and is then switched off. The ions are lost only by recombination, with coefficient r=1.67×10−6r = 1.67\times10^{-6} cm3^3/s. How long after the switch-off does the ion concentration fall by a factor η=2.0\eta = 2.0?

Given
  • n˙i=3.5×109\dot n_i = 3.5\times10^{9} cm−3^{-3}s−1^{-1}
  • r=1.67×10−6r = 1.67\times10^{-6} cm3^3/s
  • η=2.0\eta = 2.0
Find
the time tt for the concentration to fall η\eta times

Solution

Step 1The starting concentration

Having run long enough, the ionizer has brought the air to equilibrium, where the pairs it makes are lost by recombination at the same rate, rn02r n_0^2 per unit volume per second:

n˙i=rn02\dot n_i = r n_0^2

So the concentration at the moment of switching off is

n0=n˙ir(1) n_0 = \sqrt{\frac{\dot n_i}{r}} \tag{1}
Step 2The decay by recombination

After the switch-off nothing is produced. Positive and negative ions have the same concentration nn and recombine in pairs, so

dndt=−rn2\frac{dn}{dt} = -r n^2

Separating the variables and integrating from the switch-off, where the concentration is n0n_0:

∫n0ndn′n′2=−r∫0tdt′⟹1n−1n0=rt\int_{n_0}^{n}\frac{dn'}{n'^2} = -r\int_0^t dt' \quad\Longrightarrow\quad \frac{1}{n} - \frac{1}{n_0} = r t

Solving for the concentration:

n(t)=n01+rn0t(2) n(t) = \frac{n_0}{1 + r n_0 t} \tag{2}
The ion concentration after the switch-off (solid): it falls to n_0/2 at t = τ, with τ = 1/(r n_0). Dashed: an exponential with the same initial slope, which the recombination decay lags behind.

The ion concentration after the switch-off (solid): it falls to n0/2n_0/2 at t=τt = \tau, with τ=1/(rn0)\tau = 1/(r n_0). Dashed: an exponential with the same initial slope, which the recombination decay lags behind.

Step 3The time to fall η\eta times

Setting n=n0/ηn = n_0/\eta in (2):

1+rn0t=η⟹t=η−1rn01 + r n_0 t = \eta \quad\Longrightarrow\quad t = \frac{\eta - 1}{r n_0}

Substituting n0n_0 from (1) into the product rn0r n_0:

rn0=rn˙ir=rn˙ir n_0 = r\sqrt{\frac{\dot n_i}{r}} = \sqrt{r\dot n_i}

Numerically, with η=2.0\eta = 2.0:

t=2.0−11.67×10−6×3.5×109 st = \frac{2.0 - 1}{\sqrt{1.67\times10^{-6}\times 3.5\times10^{9}}}\ \text{s}
t=η−1rn˙i=13 mst = \frac{\eta - 1}{\sqrt{r\dot n_i}} = 13\ \text{ms}
(3)
Check

Units: cm3s−1⋅cm−3s−1\text{cm}^3\text{s}^{-1}\cdot\text{cm}^{-3}\text{s}^{-1} is s−2\text{s}^{-2}, whose root is a rate, so tt is a time. Limits: η=1\eta = 1 gives t=0t = 0, as it must; and the time grows only linearly with η\eta, so a hundredfold fall takes 99/rn˙i99/\sqrt{r\dot n_i}, about 1.31.3 s, about a hundred times the halving time, where an exponential decay would need fewer than seven half lives.

Remark

Recombination is a two-body loss, so the rate per ion, rnr n, drops as the ions thin out: each halving takes twice as long as the one before. There is no fixed half life, unlike a radioactive decay or an RCRC discharge, and the slow tail in the graph is why a gas stays weakly conducting long after its ionizer stops.