Ion Time-of-Flight Calculator

How fast an ion or electron moves after it is accelerated through a voltage, and how long it takes to fly down a tube to a detector. This is the arithmetic behind a time-of-flight (TOF) mass spectrometer, an electron energy analyzer, or any beamline timing. The calculator works relativistically and shows the classical answer beside it, models one- and two-field (Wiley–McLaren) ion sources with their space focus, and estimates the mass resolution.

Key terms
Accelerating voltage
The potential difference an ion falls through. An ion with charge number z gains z electronvolts of kinetic energy per volt: a singly charged ion accelerated through 2 kV carries 2 keV, a doubly charged one 4 keV.
m/z
The ion’s mass in atomic mass units divided by its charge number. Ions with the same m/z fly in exactly the same way through electric fields, so a TOF spectrum is a spectrum of m/z, not of mass.
Drift region
The field-free tube between the ion source and the detector, where ions coast at constant speed and separate by m/z. Flight time grows with the square root of m/z.
Wiley–McLaren source
An ion source with two accelerating fields in a row: a weak field where the ions are formed and a stronger one after it. Choosing the ratio of the two fields brings ions that started at different positions to the detector at the same time.
Space focus
The distance from the source at which ions of one m/z that started slightly closer to or farther from the exit arrive together, to first order. Ions that start farther back gain more energy and catch up with those that started ahead of them.
Turn-around time
An ion that starts moving away from the detector must first be stopped and turned around by the source field; it then trails an identical ion that started toward the detector by this time. It limits the resolution of any linear TOF.
Mass resolving power
m/Δm: how well neighboring masses are separated. In a TOF it is half the flight time divided by the width of a peak in time, t/(2Δt), because t grows as the square root of m.

Inputs

Ion source

Drift only: every ion enters the tube with the full energy z·e·U, and the time spent accelerating is ignored. One field or two fields: ions start at rest inside the source, and the time in the source and the energy each ion gains from its start position are included.

Particle

Fills in the mass and charge. To work from an m/z value, choose Custom and enter it as the mass with charge 1: flight times depend only on m/z.

u
e

Ion masses are atomic masses minus z electron masses (AME2020, CODATA 2022). Negative z (electrons, anions) only reverses the polarity of the voltages; enter the voltages as magnitudes.

Acceleration

The potential difference between where the ion starts and the drift tube, as a magnitude.

The field-free distance from the last grid of the source to the detector. Bench-top linear TOFs use 0.3–2 m.

Resolution (optional)

The width that the source pulse or laser pulse, the detector, and the digitizer add to every peak. A microchannel-plate detector with a fast digitizer adds about 0.5–2 ns.

Results

Flight time to the detectort—
Nonrelativistic flight time—
Speed in the drift regionv—
Kinetic energy in the drift regionK—
Lorentz factorγ—
Mass-to-charge ratiom/z—
Time in the source—
Time in the drift region—
Peak spacing for 1 u—The time between this peak and the one for an ion 1 u heavier with the same charge.
First-order space focus, drift lengthLf—
Extractor voltage for focus at L—
Arrival spread from Δs—
Turn-around timeΔtta—
Mass resolving powerm/Δm—

Layout

Time-of-flight layout

Flight time vs. m/z

How it is calculated

An ion of mass \(m\) and charge \(q = ze\) that falls through a potential difference \(U\) gains the kinetic energy \(K = |q|U\). Relativistically its Lorentz factor, speed, and momentum are

\[ \begin{gathered} \gamma = 1 + \frac{K}{mc^2}, \qquad v = c\sqrt{1 - \frac{1}{\gamma^2}}, \qquad p = \gamma m v = \frac{\sqrt{K^2 + 2Kmc^2}}{c}, \\[8pt] t_\text{drift} = \frac{L}{v}. \end{gathered} \]

The classical values are \(v = \sqrt{2K/m}\) and \(p = mv\). The two agree while \(K \ll mc^2\) (938 MeV for a proton, 511 keV for an electron): the relativistic flight time is longer by about \(\tfrac34 K/mc^2\), 0.04 ppm for Ar⁺ at 2 keV but 0.3 % for an electron at 2 keV.

In a uniform field \(E\) the force \(|q|E\) changes the momentum at a constant rate, both classically and relativistically, so the time to cross a field region is the momentum gained divided by the force. For the two-field source, with the back plate (repeller) at \(V_R\), the extraction grid at \(V_E\), and ground at the second grid,

\[ \begin{gathered} E_s = \frac{V_R - V_E}{S}, \qquad E_d = \frac{V_E}{d}, \qquad K_1 = |q| E_s s_0, \qquad K = K_1 + |q| V_E, \\[8pt] t = \frac{p_1}{|q| E_s} + \frac{p - p_1}{|q| E_d} + \frac{L}{v}, \end{gathered} \]

where \(p_1\) is the momentum at the extraction grid and \(p\), \(v\) the momentum and speed in the drift region. The one-field source is the same with no second gap: \(E_s = U/S\) and \(K = |q|E_s s_0\). The peak spacing for 1 u is the exact difference of two flight times, and every result uses \(|z|\).

Space focus and the Wiley–McLaren condition

An ion that starts farther from the exit grid (larger \(s_0\)) spends longer in the source but gains more energy, so it catches up with an ion that started ahead of it. The catch-up happens at the drift length where the arrival time does not change with \(s_0\) to first order, \(\partial t/\partial s_0 = 0\). The calculator finds it by differentiating the full relativistic \(t\) numerically; classically Wiley and McLaren’s result is

\[ L_f = 2 s_0\, k_0^{3/2} \left( 1 - \frac{d}{s_0\,(k_0 + \sqrt{k_0})} \right), \qquad k_0 = \frac{K}{K_1} = \frac{s_0 E_s + d E_d}{s_0 E_s}. \]

With one field (\(d = 0\), \(k_0 = 1\)) the focus is at \(L_f = 2s_0\), only a centimeter or two from the source: the reason a single-field source cannot space-focus a tube of useful length. With two fields, raising \(V_E\) toward \(V_R\) weakens the source field, raises \(k_0\), and pushes the focus out to the detector. The calculator solves for that extractor voltage by bisection. The arrival spread for a start-position spread \(\Delta s\) is the full range of arrival times over the interval \(s_0 \pm \Delta s/2\), evaluated at 41 points.

An ion born with energy \(\varepsilon_0\) moving away from the detector is stopped and turned around by the source field. It passes its starting point with the same speed as an identical ion that started toward the detector, but later by the turn-around time

\[ \Delta t_\text{ta} = \frac{2\sqrt{2 m \varepsilon_0}}{|q| E_s}. \]

No static field arrangement removes it; delayed extraction and reflectrons reduce its effect.

Mass resolution

Because \(t \propto \sqrt{m/z}\) in any electrostatic instrument, \(\Delta m/m = 2\Delta t/t\), so the resolving power is

\[ \frac{m}{\Delta m} = \frac{t}{2\,\Delta t}, \qquad \Delta t = \sqrt{\Delta t_\text{inst}^2 + \Delta t_s^2 + \Delta t_\text{ta}^2}. \]

The calculator adds the timing spread, the arrival spread from \(\Delta s\), and the turn-around time in quadrature, as if each were an independent peak width. That is an estimate: the start-position spread at the focus gives a skewed peak, and the turn-around time is the full separation of the two extreme ions, not a FWHM.

Typical values

Flight times through a 1 m drift region (drift only, singly charged), from the formulas above:

Particlem/zU = 2 kVU = 20 kV
Electron0.00054937.81 ns12.27 ns
Proton1.0071.616 µs0.511 µs
⁴⁰Ar⁺39.9610.18 µs3.218 µs
¹³²Xe⁺131.918.49 µs5.846 µs
C₆₀⁺720.043.19 µs13.66 µs
A 10 kDa protein, 1+10 000161.0 µs50.90 µs

Every factor of 10 in voltage shortens the flight time by \(\sqrt{10} \approx 3.16\); every factor of 4 in m/z doubles it. For the electron the relativistic correction is 0.29 % at 2 kV and 2.9 % at 20 kV.

Worked example

The default inputs: ⁴⁰Ar⁺ (m = 39.9618 u, z = 1) accelerated through 2000 V into a 1 m drift tube, with 2 ns of timing spread.

Switch to Two fields (Wiley–McLaren) to see the source defaults: VR = 2000 V, VE = 1830 V, S = 10 mm, s₀ = 5 mm, d = 10 mm. The source field is 170 V/10 mm = 17 V/mm, so an ion starting at s₀ gains 85 eV before the extraction grid and 1915 eV in all, and k₀ = 1915/85 = 22.5. The formula above gives Lf = 0.991 m, so the detector at 1 m sits 0.9 % beyond the focus; the calculator’s extractor voltage for focus at exactly 1 m is 1830.9 V. The ion arrives at 11.064 µs (0.494 µs in the first gap, 0.172 µs in the second). The 0.5 mm start-position spread adds only 0.24 ns, so with the 2 ns timing spread m/Δm ≈ 2750. An initial energy of 0.04 eV would add a turn-around time of 21 ns and cut m/Δm to about 260: in a linear TOF with a weak source field, thermal motion, not the start position, sets the resolution.

Assumptions and limits

References

Primary sources and standards

  1. W. C. Wiley and I. H. McLaren, “Time-of-flight mass spectrometer with improved resolution,” Rev. Sci. Instrum. 26, 1150–1157 (1955), doi:10.1063/1.1715212 (the two-field source, its space-focus condition, and the turn-around time).
  2. U. Boesl, “Time-of-flight mass spectrometry: Introduction to the basics,” Mass Spectrom. Rev. 36, 86–109 (2017), doi:10.1002/mas.21520 (one- and two-stage sources, flight-time broadening, and mass resolution).
  3. J. D. Jackson, Classical Electrodynamics, 3rd ed., Wiley (1999), ch. 11–12 (relativistic kinematics; motion in a uniform electric field, where the momentum grows linearly in time).
  4. P. J. Mohr, D. B. Newell, B. N. Taylor, and E. Tiesinga, “CODATA recommended values of the fundamental physical constants: 2022,” Rev. Mod. Phys. 97, 025002 (2025), doi:10.1103/RevModPhys.97.025002 (e, u, the electron, proton, deuteron, and alpha-particle masses, and c).
  5. M. Wang, W. J. Huang, F. G. Kondev, G. Audi, and S. Naimi, “The AME 2020 atomic mass evaluation (II). Tables, graphs and references,” Chin. Phys. C 45, 030003 (2021), doi:10.1088/1674-1137/abddaf (the atomic masses of the preset ions).

Free introductions