Isotope Patterns and Elemental Clues

A mass spectrum is a bar-code of isotopes before it is anything else. Look at the peak cluster near the molecular ion of bromobenzene (C₆H₅Br). You see two nearly equal-height bars two mass units apart. This is not measurement noise. It is the fingerprint of bromine itself. Nature ships bromine as almost exactly a 1:1 mix of ⁷⁹Br and ⁸¹Br. Long before a chemist runs a high-resolution instrument or a fragmentation study, that simple pattern already answers a question. Does this molecule contain bromine? The answer comes with near certainty. This matters because molecular formula is the first fact any structure elucidation needs. Isotope ratios hand it over almost for free. Every element on the periodic table carries its own signed isotopic ratio. Carbon has a small but reliable 1.1% ¹³C tax on every atom. Chlorine shows a 3:1 split between mass 35 and 37. Sulfur has a characteristic 4.4% bump two mass units up. These ratios are constants of nature. They stay essentially unchanged by chemical environment. So they survive intact into any spectrum, regardless of ionization method. This holds for electron ionization, electrospray, or MALDI. It works as long as the instrument resolves the isotopologues. The big picture is this: a mass spectrum's M, M+1, and M+2 peaks encode a small system of linear equations. Solve them and you recover useful information, with quantifiable uncertainty. You learn how many carbons are present. You learn whether halogens or sulfur are present, and how many. This is the same logic used in The Mass Spectrum and the m/z Axis to define the axis itself. Now it becomes a chemical detective tool instead. Intensities become data. Natural abundance tables become the model. Formula deduction becomes a fitting problem with real statistical weight behind it. This article builds the isotope-pattern toolkit from first principles. It covers the binomial and multinomial statistics behind isotope clusters. It gives the working formulas for %M+1 and %M+2 relative to M. It presents the diagnostic doublets and triplets for Cl, Br, and S. It also works through three fully worked examples. These range from a simple hydrocarbon, to a brominated drug fragment, to a sulfur/chlorine-containing exam problem. It sets up naturally for interpreting fragmentation patterns. It also leads into the accurate-mass methods in high-resolution and accurate mass measurements. There, isotope reasoning and mass-defect reasoning merge into one unified craft.

Key takeaways
  • Every element is a fixed natural mixture of isotopes; these ratios do not change with chemical bonding and are tabulated constants.
  • The M+1 peak height (relative to M = 100%) is built mainly from ¹³C content: %M+1 ≈ 1.1 × (number of carbons), plus small ¹⁵N/²H/¹⁷O terms.
  • The M+2 peak is diagnostic: large M+2 (~33% of M) signals one chlorine; ~98% of M signals one bromine; ~4.4% signals one sulfur.
  • Two chlorines give a 9:6:1 triplet; two bromines give a 1:2:1 triplet; mixed Cl/Br give characteristic multiplets — all derivable from the binomial expansion of isotope fractions.
  • Isotope-pattern fitting is a statistical inference problem: measured intensities are compared to a multinomial model, and formula candidates are ranked by chi-square-like residuals.
  • These patterns are read directly off unit-resolution spectra, no high-resolution instrument required, making them the fastest formula clue available.
  • Key idea: the height and spacing of the M+1/M+2 satellite peaks are a direct, calculable readout of an element's own natural isotope ratios, turning a spectrum into an elemental census.

Where this lives today

Isotope-pattern reading remains the fastest, cheapest formula clue in a working mass spectrometry laboratory, and it has only grown more central as instruments and databases have matured. In pharmaceutical impurity and metabolite identification, LC-MS screening routinely uses the bromine and chlorine doublet/triplet signatures to flag halogenated degradants or metabolites within seconds, well before any high-resolution accurate-mass formula search is attempted — the isotope pattern narrows the candidate list before the accurate mass narrows it further, exactly the workflow anticipated in high-resolution and accurate mass methods. In environmental and forensic chemistry, brominated flame retardants and chlorinated pesticides are screened in complex matrices largely by recognizing their unmistakable M/M+2/M+4 fingerprints against a noisy chemical background, because the pattern survives even when exact mass accuracy is compromised by low ion counts. In proteomics and peptide mass fingerprinting, software that assigns peptide identities from tandem mass spectra explicitly models the expected isotope envelope (dominated by ¹³C at high carbon counts) to distinguish real peptide monoisotopic peaks from noise and to verify charge state, since a +2 charge state halves the apparent isotope spacing to 0.5 Da rather than 1 Da — a direct extension of the binomial carbon model built here into the multiply-charged regime, a natural bridge toward interpreting fragmentation patterns in tandem MS. In isotope-ratio mass spectrometry (a distinct but related discipline), the same natural-abundance bookkeeping is inverted: instead of using known ratios to deduce formula, tiny deviations from standard terrestrial ratios (in carbon, nitrogen, oxygen, sulfur) are measured with exquisite precision to trace metabolic pathways, authenticate food origin, and reconstruct paleoclimate records — proof that the humble M+1/M+2 satellite peak, first explained by Aston a century ago, still carries information at both the routine-identification and frontier-research ends of modern chemistry.

Ask the tutor: If a spectrum shows M+2 at 45% and M+4 at 3.5%, could this be one bromine plus something else, or does it demand a different explanation entirely?

The physical picture

A mass spectrometer separates ions by mass-to-charge ratio, but "mass" here means the mass of a particular isotopic composition of a molecule, not some averaged textbook molecular weight. A molecule of chlorobenzene is not one species but effectively two: C₆H₅³⁵Cl (mass 112) and C₆H₅³⁷Cl (mass 114), present in nature's fixed ratio because every chlorine atom, wherever it sits in a molecule, is drawn from the same terrestrial isotopic reservoir. The reservoir ratios are known to high precision: ³⁵Cl is 75.77% and ³⁷Cl is 24.23% of all chlorine atoms — call it 3.05 : 1, conventionally quoted as 100 : 32.5 relative to the lighter isotope taken as 100. Because chlorination doesn't discriminate isotopes (chemistry cares about protons and electrons, isotopes differ only in neutrons), every C–Cl bond in every molecule on Earth reproduces this same ratio. The consequence for a mass spectrum is immediate: a molecular ion containing one chlorine atom does not appear as a single peak. It appears as two peaks, at M and M+2, in a fixed area ratio of about 3:1. This is not fragmentation, not a contaminant, not noise — it is the isotope pattern of chlorine itself, printed onto the spectrum. Carbon does something similar but far more subtly: 98.93% of carbon atoms are ¹²C, but 1.07% are ¹³C, one mass unit heavier. A molecule with many carbons has many independent chances to contain one ¹³C, so the M+1 peak grows in proportion to the carbon count — small for a single carbon, substantial for a fifteen-carbon steroid. The clean way to see all of this is the binomial/multinomial model: each atom of a polyisotopic element is an independent trial with a fixed probability of "success" (being the heavy isotope). Summing over all atoms of that element in the molecule, and convolving over multiple elements, generates the entire isotope cluster from a small number of natural constants.

Figure
Chlorine's isotope doublet versus bromine's near-1:1 doublet, drawn to scale
Cl-containing ion: M and M+2 (≈100:32.5) M (35Cl) M+2 (37Cl) 100% 32.5% Br-containing ion: M and M+2 (≈100:97.3) M (79Br) M+2 (81Br)

The bromine case is even more visually striking: since ⁷⁹Br (50.69%) and ⁸¹Br (49.31%) are almost equally abundant, any bromine-containing ion shows an almost exactly 1:1 doublet at M and M+2 — one of the most recognizable signatures in all of mass spectrometry, visible by eye without any calculation. Sulfur adds a third flavor: ³²S is 94.99% abundant, but ³⁴S sits at 4.25% (with a negligible 0.75% ³³S folded mostly into rounding). This produces a modest but diagnostic M+2 peak around 4–4.5% of M — far smaller than chlorine's 32.5% or bromine's 97%, but distinguishable from the "no heavy-halogen, no sulfur" background once carbon count is accounted for.

Figure
Building the M, M+1, M+2 cluster for a small organic ion, showing the carbon (¹³C) contribution stacking with a heteroatom contribution
C₈ fragment (no heteroatom) vs. same fragment +1 S M M+1 (8.6%) M+2 (~0%) M M+1 (8.6%) M+2 (~4.5% from S)

The key visual lesson: M+1 tracks carbon count almost regardless of what else is in the molecule, while M+2 is where halogens and sulfur announce themselves, because ¹³C alone essentially never produces meaningful M+2 intensity in typical organic ions (it would take two ¹³C atoms simultaneously, a much rarer event captured by the binomial cross-term).

The core ideas

For a single element with a light isotope of abundance aa (set to 100) and a heavy isotope of abundance bb (relative to a=100a=100), nn atoms of that element in a molecule generate a binomial isotope pattern. The relative intensity of the peak with kk heavy atoms out of nn is

Ik=(nk) a n−k b kI_k = \binom{n}{k}\,a^{\,n-k}\,b^{\,k}

where:

  • nn = number of atoms of that element in the ion,
  • aa = relative abundance of the light isotope (normalized, e.g. 100 for ³⁵Cl),
  • bb = relative abundance of the heavy isotope (e.g. 32.5 for ³⁷Cl relative to a=100),
  • IkI_k = relative intensity of the peak kk mass units above the all-light isotopologue (before normalizing the whole set to 100 at the base peak).

Reading: this is exactly the term-by-term binomial expansion of (a+b)n(a+b)^n; each term corresponds to a specific count of heavy isotopes, and the spacing between terms is fixed by the mass difference between isotopes (2 amu for Cl, Br, S). For carbon specifically, with 13C^{13}\mathrm{C} natural abundance p=0.0107p = 0.0107 (1.07%) and nCn_C carbon atoms, the practical approximate formula (valid because pnC≪1p n_C \ll 1 for typical molecules) is

%M+1≈1.1 nC  +  0.36 nN  +  0.02 nO  +  0.016 nS\%M{+}1 \approx 1.1\, n_C \;+\; 0.36\, n_N \;+\; 0.02\, n_O \;+\; 0.016 \, n_S

where:

  • %M+1\%M{+}1 = intensity of the M+1 peak as a percentage of the M peak (M ≡ 100%),
  • nC,nN,nO,nSn_C, n_N, n_O, n_S = number of carbon, nitrogen, oxygen, and sulfur atoms respectively,
  • the coefficients are 100×100\times the natural abundance ratio of each element's next-heaviest isotope to its main isotope (¹³C/¹²C = 1.08%, ¹⁵N/¹⁴N = 0.37%, ¹⁷O/¹⁶O ≈ 0.04% giving the small oxygen term, ³³S/³²S ≈ 0.75%).

Reading: M+1 is almost entirely a carbon-counting tool; nitrogen, oxygen, and sulfur nudge it slightly but rarely change the carbon count you'd infer by rounding. This is the standard rule chemists use to estimate nCn_C directly from %M+1/1.1. For the M+2 peak, the dominant contributions are the heavy isotopes with a 2-amu mass defect — ³⁷Cl, ⁸¹Br, ³⁴S, ¹⁸O — plus the small binomial cross-term from two ¹³C atoms simultaneously:

%M+2  ≈  (nC2)(1.1100)2×100  +  0.20 nO  +  4.4 nS  +  32.5 nCl  +  97.3 nBr\%M{+}2 \;\approx\; \binom{n_C}{2}\left(\frac{1.1}{100}\right)^{2}\times 100 \;+\; 0.20\,n_O \;+\; 4.4\,n_S \;+\; 32.5\,n_{Cl} \;+\; 97.3\,n_{Br}

where:

  • the first term is the two-¹³C binomial cross-term (usually negligible unless nCn_C is large, e.g. > 20),
  • 0.20 nO0.20\,n_O comes from ¹⁸O at 0.205% natural abundance,
  • 4.4 nS4.4\,n_S and 32.5 nCl32.5\,n_{Cl} and 97.3 nBr97.3\,n_{Br} come directly from the ³⁴S/³²S, ³⁷Cl/³⁵Cl, ⁸¹Br/⁷⁹Br natural abundance ratios (each multiplied by the number of atoms, to first order, when only one halogen atom is present — for two or more, the full binomial/multinomial expansion below is needed).

Reading: for molecules without heavy halogens or excess sulfur, M+2 is tiny (a fraction of a percent) and easily lost in noise; the moment M+2 is a clearly resolvable few percent or more of M, a halogen or sulfur atom is essentially certain to be present.

Derivation AThe chlorine doublet from the binomial isotope model

Goal. Derive the exact M : M+2 ratio for a molecular ion containing exactly one chlorine atom, and show how it generalizes to two chlorines.

Step 1 — set up the two-state distribution. Chlorine exists in nature as ³⁵Cl (abundance p35=0.7577p_{35} = 0.7577) and ³⁷Cl (abundance p37=0.2423p_{37} = 0.2423). For a molecule containing exactly one Cl atom, the ion is either the "all-³⁵Cl" isotopologue or the "one-³⁷Cl" isotopologue — there is no intermediate, because n=1n=1 atom can only be light or heavy, nothing else. So the intensity ratio is simply

I(M+2)I(M)=p37p35=0.24230.7577=0.3199\frac{I(M+2)}{I(M)} = \frac{p_{37}}{p_{35}} = \frac{0.2423}{0.7577} = 0.3199

Step 2 — express as a percentage. Multiplying by 100 gives %M+2 ≈ 32.0% of the M peak height, consistent with the commonly quoted 100:32.5 (small variation in tabulated abundance values, e.g. IUPAC vs. CRC tables, shifts the last digit).

Step 3 — extend to two chlorine atoms via the binomial expansion. With n=2n=2 independent Cl atoms, each drawn from the same (p35,p37)(p_{35}, p_{37}) distribution, the probability of kk heavy atoms out of 2 is binomial:

Ik=(2k) p35 2−k p37 k,k=0,1,2I_k = \binom{2}{k}\,p_{35}^{\,2-k}\,p_{37}^{\,k}, \qquad k = 0,1,2

Expanding term by term (before normalizing to I0=100I_0 = 100):

I0:I1:I2=p352  :  2 p35p37  :  p372I_0 : I_1 : I_2 = p_{35}^2 \;:\; 2\,p_{35}p_{37} \;:\; p_{37}^2

Step 4 — plug in numbers. With p35=0.7577p_{35}=0.7577, p37=0.2423p_{37}=0.2423:

p352=0.5741,2p35p37=0.3672,p372=0.0587p_{35}^2 = 0.5741,\quad 2p_{35}p_{37} = 0.3672,\quad p_{37}^2 = 0.0587

Normalizing to the M peak (p352→100p_{35}^2 \to 100):

I0:I1:I2=100:64.0:10.2I_0 : I_1 : I_2 = 100 : 64.0 : 10.2

This is the well-known 9:6:1 pattern (rounded ratios), the classic signature read directly off a spectrum for a dichloro-compound: three peaks at M, M+2, M+4 with heights roughly in the ratio 9:6:1 (equivalently 100:64:10.2, matching the more careful binomial values above).

Step 5 — limiting case check. As the number of chlorine atoms n→0n \to 0, all binomial weight collapses onto I0I_0 and the M+2/M+4 peaks vanish — consistent with a chlorine-free molecule showing no such satellites. For n=1n=1, the formula correctly reduces to the simple two-term ratio of Step 1, since (10)=(11)=1\binom{1}{0}=\binom{1}{1}=1 and there is no k=2k=2 term available.

Go deeper · Derivation B · 1 minGeneral multinomial formula for mixed halogen/sulfur/carbon clusters
Derivation BGeneral multinomial formula for mixed halogen/sulfur/carbon clusters

Goal. Build the complete isotope cluster for a molecule containing several polyisotopic elements simultaneously (e.g., C, Cl, Br, S together), and show it reduces correctly to the single-element cases already derived.

Step 1 — treat each element as an independent generating function. Represent each element's isotope distribution as a polynomial in a bookkeeping variable xx, where the exponent of xx tracks mass units above the lightest isotope. For carbon (2 isotopes, spacing 1 amu):

GC(x)=(0.9893+0.0107 x)nCG_C(x) = (0.9893 + 0.0107\,x)^{n_C}

For chlorine (2 isotopes, spacing 2 amu, tracked as x2x^2):

GCl(x)=(0.7577+0.2423 x2)nClG_{Cl}(x) = (0.7577 + 0.2423\,x^{2})^{n_{Cl}}

Analogous generating functions apply for Br (0.5069+0.4931 x20.5069 + 0.4931\,x^2) raised to nBrn_{Br}, and S using its three isotopes (0.9499+0.0075 x+0.0425 x2)nS(0.9499 + 0.0075\,x + 0.0425\,x^{2})^{n_S} (the ³³S term is small but included for completeness).

Step 2 — the whole-molecule pattern is the product of the generating functions. Because isotope substitution at each element is statistically independent, the combined distribution is

Gtotal(x)=GC(x)⋅GCl(x)⋅GBr(x)⋅GS(x)⋯G_{\text{total}}(x) = G_C(x)\cdot G_{Cl}(x)\cdot G_{Br}(x)\cdot G_S(x)\cdots

and the relative intensity of the M+mm peak is the coefficient of xmx^m in the expanded polynomial, all normalized so the largest term (usually x0x^0, the "all-light-isotope" ion) is scaled to 100%.

Step 3 — worked reduction for CH₃Br (methyl bromide), a clean 2-element case. Here nC=1n_C = 1, nBr=1n_{Br}=1, no Cl or S. The relevant product to second order is

G(x)=(0.9893+0.0107x)(0.5069+0.4931x2)G(x) = (0.9893 + 0.0107x)(0.5069 + 0.4931x^2)

Expanding:

G(x)=(0.9893)(0.5069)+(0.0107)(0.5069) x+(0.9893)(0.4931) x2+(0.0107)(0.4931) x3G(x) = (0.9893)(0.5069) + (0.0107)(0.5069)\,x + (0.9893)(0.4931)\,x^2 + (0.0107)(0.4931)\,x^3
G(x)=0.5015+0.00542 x+0.4879 x2+0.00528 x3G(x) = 0.5015 + 0.00542\,x + 0.4879\,x^2 + 0.00528\,x^3

Step 4 — normalize to M = 100. Divide every coefficient by 0.5015 and multiply by 100:

M:M+1:M+2:M+3=100:1.08:97.3:1.05M : M{+}1 : M{+}2 : M{+}3 = 100 : 1.08 : 97.3 : 1.05

This matches the textbook expectation exactly: CH₃Br shows a near 1:1 M/M+2 doublet (bromine's signature) with a barely-visible 1% M+1 shoulder from the single carbon — small enough to often be folded into peak-width noise at unit resolution, but present.

Step 5 — limiting case check. Setting nBr=0n_{Br}=0 collapses GBr(x)→1G_{Br}(x) \to 1, and the whole pattern reduces to the pure-carbon M+1 pattern (1.1% per carbon), recovering the standard carbon-only rule. Setting nC=0n_C=0 instead recovers the pure single-bromine doublet from the earlier derivation, confirming the generating-function method contains both single-element results as special cases.

Worked examples

Example 1 (foundational) — Reading the carbon count from M+1 in a simple ketone. A compound believed to be a small aliphatic ketone shows, in its EI mass spectrum, a molecular ion at m/z 86 with relative intensity 100%, and an M+1 peak at m/z 87 with relative intensity 6.6%. Estimate the number of carbon atoms, and check consistency with a proposed formula C₅H₁₀O (MW 86.13, nominal mass 86). Using the carbon-dominated approximation (ignoring the small O contribution, since 0.02×1=0.02%0.02 \times 1 = 0.02\% is far below the measurement precision):

nC≈%M+11.1=6.61.1=6.0n_C \approx \frac{\%M{+}1}{1.1} = \frac{6.6}{1.1} = 6.0

This is a slight surprise: the simple estimate suggests 6 carbons, not 5. Checking further, C₆H₁₄ has nominal mass 86 too (14×6 + 14 = 100, not matching — recompute: C₆H₁₄ = 6(12)+14(1) = 72+14 = 86, yes it matches at nominal mass 86 also). Both C₅H₁₀O and C₆H₁₄ are nominally mass-86 formulas ("isobaric" at unit resolution), and the M+1 intensity is exactly the tool that discriminates them: 5 carbons would predict %M+1 ≈ 5.5% (plus ~0.02% from one oxygen, negligible), while 6 carbons predicts %M+1 ≈ 6.6%. The observed 6.6% points cleanly to the six-carbon, oxygen-free hydrocarbon, C₆H₁₄ (2-methylpentane or a related isomer), not the ketone originally guessed.

Limiting-case sanity check: if the compound had zero carbons (an impossible but instructive limit for an organic ion), the model predicts %M+1 → 0, correctly reflecting that only heteroatoms like N or O could produce any M+1 signal at all — a check that confirms the formula's structure (proportionality to nCn_C, vanishing at nC=0n_C=0) rather than an accidental fit to this one data point.

Example 2 (applied) — Identifying a brominated fragment ion in a pharmaceutical impurity study. A liquid chromatography–mass spectrometry (LC-MS) trace of a drug degradation product shows a fragment ion cluster with peaks at m/z 261 (100%), m/z 262 (12.9%), and m/z 263 (98.1%). No peak is observed at m/z 264 above baseline. Determine the likely elemental composition clue (halogen content and approximate carbon count).

Step 1 — read M+2. %M+2 = 98.1%, essentially matching the bromine signature (theoretical single-Br doublet ≈ 97.3%, allowing for ±2% measurement scatter typical of ion-counting statistics at these intensities). This is far too large to be sulfur (expected ~4.4%) and not the 32.5% expected for a single chlorine — it is unambiguously one bromine atom.

Step 2 — read M+1 to extract carbon count, correcting for the small carbon contribution to bromine's own isotope pattern (Br itself contributes zero M+1, since ⁷⁹Br → ⁸¹Br is a 2-amu jump, not 1-amu). So the entire 12.9% M+1 signal is attributable to carbon (plus any N/O, assumed absent or small here):

nC≈12.91.1=11.7≈12n_C \approx \frac{12.9}{1.1} = 11.7 \approx 12

Step 3 — propose a formula. A fragment with 12 carbons, one bromine, at nominal mass 261: subtracting Br (mass 79 for the lighter isotope) leaves 261−79=182261 - 79 = 182 mass units for the C₁₂ + remaining H/N/O skeleton. With 12 carbons accounting for 12×12=14412\times12 = 144 mass units, the remaining 182−144=38182-144=38 mass units must be distributed among H, N, O — consistent, for instance, with C₁₂H₁₄NOBr (144 + 14 + 14 + 16 + 79 = 267, needs adjustment) or more simply C₁₂H₂₂Br (144+22+79 = 245, too low) — the precise remaining-atom assignment requires the accurate mass or further fragmentation data, but the halogen and carbon-count clues are already secured from the isotope pattern alone, exactly the handoff point where high-resolution and accurate mass measurement would pin down the exact formula.

Limiting-case check: if the observed M+2 had instead been near 32%, the same reasoning machinery (Step 1) would have flagged chlorine instead of bromine, showing the method's diagnostic power depends entirely on which characteristic ratio the data matches, not on any assumption baked in beforehand.

Example 3 (exam/competition level) — Disentangling a mixed Cl₂S-containing ion. A mass spectrum shows a cluster of four peaks with the following relative intensities: M (m/z 200) = 100%, M+2 = 69.5%, M+4 = 16.4%, M+6 = 1.3%. A student proposes the ion contains either (a) two chlorine atoms and one sulfur atom, or (b) three chlorine atoms. Use the multinomial isotope model to decide which is more consistent with the data.

Step 1 — model (b): three chlorines only. Using the binomial expansion for n=3n=3 with p35=0.7577p_{35}=0.7577, p37=0.2423p_{37}=0.2423:

Ik=(3k) p35 3−kp37 kI_k = \binom{3}{k}\,p_{35}^{\,3-k}p_{37}^{\,k}
I0=0.75773=0.4349,I1=3(0.7577)2(0.2423)=0.4175I_0 = 0.7577^3 = 0.4349,\quad I_1 = 3(0.7577)^2(0.2423)=0.4175
I2=3(0.7577)(0.2423)2=0.1336,I3=0.24233=0.01423I_2 = 3(0.7577)(0.2423)^2 = 0.1336,\quad I_3 = 0.2423^3 = 0.01423

Normalizing to I0=100I_0=100: M:M+2:M+4:M+6=100:96.0:30.7:3.3M{:}M{+}2{:}M{+}4{:}M{+}6 = 100:96.0:30.7:3.3. This predicts M+2 ≈ 96%, far above the observed 69.5%. Three chlorines is inconsistent with the data.

Step 2 — model (a): two chlorines plus one sulfur. Build the combined generating function using the Cl₂ binomial result already derived (100:64.0:10.2 for M:M+2:M+4, from the first derivation) convolved with sulfur's own small three-isotope distribution (0.9499+0.0075x+0.0425x2)(0.9499 + 0.0075x + 0.0425x^2). To leading order, sulfur mainly contributes its own M+2 term (4.25%/94.99% ≈ 4.47% relative), which convolves into the chlorine pattern by shifting a small fraction of each chlorine-cluster peak two mass units higher. Approximating the combination as a product of the two normalized distributions:

%M+2combined≈%M+2Cl2+%M+2S×(I0Cl2100)≈64.0+4.47=68.5\%M{+}2_{\text{combined}} \approx \%M{+}2_{Cl_2} + \%M{+}2_S \times \left(\frac{I_0^{Cl_2}}{100}\right) \approx 64.0 + 4.47 = 68.5
%M+4combined≈%M+4Cl2+%M+2Cl2×4.47100≈10.2+64.0×0.0447≈10.2+2.86=13.1\%M{+}4_{\text{combined}} \approx \%M{+}4_{Cl_2} + \%M{+}2_{Cl_2}\times \frac{4.47}{100} \approx 10.2 + 64.0\times0.0447 \approx 10.2+2.86=13.1

(A fully rigorous treatment multiplies the two generating polynomials term by term rather than adding leading contributions, which would raise these estimates by a percent or two — but even this simplified convolution already lands close to the data.)

Step 3 — compare. Predicted (a): M:M+2:M+4 ≈ 100:68.5:13.1, versus observed 100:69.5:16.4. This is a far closer match than model (b)'s 100:96:30.7. The data strongly favor two chlorines plus one sulfur.

Step 4 — sanity check on M+6. Two chlorines alone predict essentially zero M+6 (the I3I_3 term requires three heavy Cl atoms, impossible with only two chlorines present), while the observed small 1.3% M+6 is naturally explained as a higher-order cross-term (one ³⁷Cl and the rarer sulfur isotope combination), confirming that a mixed-element multinomial model, not a pure-chlorine model, is required.

Limiting-case check: if the sulfur atom were removed from model (a) entirely (nS→0n_S \to 0), the prediction collapses exactly onto the pure-Cl₂ result (100:64.0:10.2) from the first derivation — confirming internal consistency between the single-element and multi-element treatments.

Limiting cases, scope & common traps

SituationWhat the pattern showsCommon pitfall
Pure hydrocarbon, few carbons (nC<10n_C < 10)Small M+1 (~1–5%), negligible M+2Mistaking natural M+1 noise for a heteroatom signal
One Cl presentM+2 ≈ 32–33% of MConfusing with two ¹³C atoms (which give <0.1% at this size)
One Br presentM+2 ≈ 97–98% of M, nearly equal doubletMisreading the doublet as two co-eluting compounds in LC-MS
One S presentM+2 ≈ 4.4% of MAttributing to noise, or confusing with two ¹³C in a 20+ carbon molecule
Two Cl present9:6:1 ratio at M, M+2, M+4Assuming a single Cl and missing the M+4 peak entirely
Very large carbon count (nC>30n_C > 30, e.g. steroids, peptides)M+1 can exceed 30%; M+2 from ¹³C₂ term becomes non-negligibleWrongly attributing large M+2 to a halogen when it is simply many carbons
Low ion statistics / weak molecular ionPoisson counting noise inflates or deflates M+1/M+2Over-interpreting a single low-intensity spectrum without averaging scans

The chemist behind it

Francis William Aston built the first high-resolution mass spectrograph in 1919 at the Cavendish Laboratory, and his most consequential discovery was not a new element but a new kind of fact about old ones: elements are mixtures of isotopes in fixed proportions, and those proportions could be read directly off photographic plates as line intensities. Aston's "whole number rule" — that isotopic masses are close to whole numbers while chemical atomic weights are not — followed directly from recognizing that chlorine's famously non-integer atomic weight (35.45) was simply the abundance-weighted average of ³⁵Cl and ³⁷Cl, present in the very ratio still used today to read the M/M+2 doublet. He was awarded the Nobel Prize in Chemistry in 1922 for the discovery, by mass spectrograph, of isotopes in a large number of nonradioactive elements, and for his enunciation of the whole-number rule. Aston's meticulous plate-intensity measurements — essentially the same kind of relative-peak-height bookkeeping used in every isotope-pattern calculation in this article — established both the technique and the philosophy: an element's isotopic signature is a measurable, tabulable constant, not an incidental detail.

Origins · 1 minFrancis William Aston and how this was discovered
FW
Francis William Aston
1877–1945 · Cavendish Laboratory, Cambridge, developer of the mass spectrograph and discoverer of isotopic natural abundance patterns

History
- 1913 — J. J. Thomson's positive-ray parabola apparatus first hints at neon having two distinct masses.
- 1919 — Aston builds his first true mass spectrograph, resolving ²⁰Ne and ²²Ne cleanly and establishing isotopes as a general phenomenon.
- 1920 — Aston measures isotopic masses for dozens of elements, formulating the whole-number rule linking isotope mass to atomic weight non-integrality.
- 1922 — Aston awarded the Nobel Prize in Chemistry for isotope discovery via mass spectrography.
- 1930s–1950s — systematic tabulation of natural isotopic abundances (Cl, Br, S, and others) becomes standard reference data for organic mass spectrometry.
Why does the M+1 peak mainly reflect carbon and not other elements?

Carbon's ¹³C isotope is both reasonably abundant (1.07%) and ubiquitous — organic molecules typically contain far more carbon atoms than nitrogen or oxygen atoms, so even though nitrogen's ¹⁵N (0.37%) and oxygen's ¹⁷O (0.04%) exist, their small abundance and typically low atom count mean their contribution to M+1 is a minor correction rather than the dominant term. Since M+1 scales linearly with atom count for each element, and carbon usually dominates atom count in organic molecules, the M+1 peak becomes, to good approximation, a direct carbon counter: %M+1 ≈ 1.1 × n_C. This approximation breaks down only for nitrogen-rich or highly oxygenated small molecules where the correction terms become comparable to the carbon term itself.

Can isotope patterns alone give an exact molecular formula?

Not by themselves — isotope-pattern reasoning at unit mass resolution narrows down elemental content (carbon count estimate, presence/absence/count of Cl, Br, S) but cannot distinguish, for instance, C₆H₁₄ from C₅H₁₀O at nominal mass 86 without the M+1 clue used together with a nominal-mass formula search, and even then several isobaric candidates often survive. The definitive step is high-resolution accurate-mass measurement, which resolves the tiny mass defects between formulas of the same nominal mass (a task belonging to accurate-mass instrumentation), using the isotope-count clues developed here as a first filter that dramatically shrinks the candidate list before that finer measurement is applied.

Why do two chlorine atoms give a 9:6:1 ratio rather than simply doubling the single-chlorine doublet?

Two independent chlorine atoms create three possible combinations — both light (³⁵Cl,³⁵Cl), one of each, or both heavy (³⁷Cl,³⁷Cl) — and the "one of each" outcome can occur in two equivalent ways (first atom heavy/second light, or vice versa), giving it double statistical weight relative to the two homogeneous outcomes. This is exactly the binomial coefficient (21)=2\binom{2}{1}=2 appearing in the expansion of (p35+p37)2(p_{35}+p_{37})^2, and plugging in natural abundances yields the characteristic 100:64:10.2 ratio (colloquially rounded to 9:6:1), not a simple doubling of the single-atom 100:32.5 doublet.

How does the isotope pattern method extend to multiply-charged ions in electrospray mass spectrometry?

For an ion carrying charge zz, the mass-to-charge spacing between isotope peaks becomes 1/z1/z instead of 1 Da, because the isotope mass difference (still 1 amu in absolute mass) is divided by the charge on the ion. A +2 charge state shows isotope peaks spaced 0.5 m/z units apart, and a +3 state shows 0.33 m/z spacing; recognizing this spacing is itself the standard method for determining the charge state of large biomolecular ions in electrospray spectra, extending the same carbon-driven binomial isotope envelope logic developed here into the proteomics regime.

Is the natural isotopic abundance of an element ever significantly different from the standard tabulated value?

For the vast majority of routine organic and pharmaceutical mass spectrometry, standard terrestrial abundances (as tabulated by IUPAC) are accurate to the precision needed for formula deduction, and any small natural variation (a fraction of a percent, arising from mass-dependent geological or biological fractionation) is far below the resolution needed to affect elemental-composition conclusions. Such small deviations are significant in the separate discipline of isotope-ratio mass spectrometry, where they are the measured quantity of interest rather than a nuisance, but for the formula-deduction purposes described here, tabulated average abundances are treated as fixed constants.

Why doesn't nitrogen's M+1 contribution get its own famous "nitrogen rule" pattern the way Cl, Br, and S do?

Nitrogen's isotope contribution to M+1 (0.37% per atom) is real but too small and non-distinctive by itself to serve as a standalone diagnostic the way Cl's 32.5% or Br's 97% M+2 does; nitrogen's true diagnostic signature in mass spectrometry is instead the "nitrogen rule" governing nominal mass parity (an odd number of nitrogen atoms gives an odd nominal molecular mass, for typical C,H,N,O compounds), which is a completely separate piece of reasoning from the isotope-abundance patterns discussed here, though both are read from the same spectrum.

What happens to the isotope pattern if the molecule contains two different halogens, say one Cl and one Br?

The two halogens' generating functions multiply together rather than simply add, producing a more complex multi-peak cluster: the M peak (both light isotopes), an M+2 peak (either the Cl or the Br atom is heavy, but not both — two contributing pathways), and an M+4 peak (both heavy simultaneously). Working the arithmetic through the product (0.7577+0.2423x2)(0.5069+0.4931x2)(0.7577+0.2423x^2)(0.5069+0.4931x^2) gives a three-peak pattern with intensities roughly 100:73:24, a shape distinct from either single-halogen doublet and itself a recognizable diagnostic for the presence of both a chlorine and a bromine in the same fragment.

Before this