Analytical Chemistry
The Karplus Equation: Reading Dihedral Angles from NMR 3J Coupling Constants
Twist two vicinal protons in a molecule to 180 degrees apart and their NMR signals split by a hefty 8 to 12 Hz; twist them to 90 degrees and that splitting nearly vanishes to almost 0 Hz. This dramatic, geometry-controlled swing is the heart of the Karplus equation, the single most-used bridge between an NMR number and a three-dimensional shape.
The Karplus equation relates the three-bond scalar coupling constant, written ³J (a "vicinal" coupling transmitted through three bonds, H–C–C–H), to the dihedral (torsion) angle θ between the two coupled nuclei. In its canonical form ³J = A cos²θ + B cosθ + C, it converts a measurable Hertz value into an angle, letting chemists deduce conformation, ring pucker, and protein backbone geometry directly from a spectrum.
- TypeEmpirical/theoretical NMR structure relation
- IntroducedMartin Karplus, 1959 (refined 1963)
- Key equation³J = A cos²θ + B cosθ + C
- Typical range0–12 Hz for vicinal ¹H–¹H
- Applies toVicinal 3-bond couplings (H–C–C–H, etc.)
- Measured by¹H NMR spectroscopy (multiplet splitting)
Interactive visualization
Press play, or step through manually. The visualization is yours to drive — try it before reading on.
Watch the 60-second explainer
A condensed visual walkthrough — narrated, captioned, under a minute.
What the Karplus Equation Is and Where It Applies
The Karplus equation is the quantitative link between a scalar coupling constant measured in an NMR spectrum and a molecular torsion angle. Specifically it governs ³J, the vicinal coupling transmitted over three bonds between two nuclei on adjacent atoms (classically H–C–C–H). Because ³J depends steeply on the H–C–C–H dihedral angle θ, a single splitting in Hertz reports on 3D geometry.
- Conformational analysis: distinguishing axial/equatorial protons, gauche vs. anti rotamers, and ring pucker in sugars and cyclohexanes.
- Structural biology: ³J(HN–Hα) values fix the protein backbone φ angle; ³J(H–H) couplings define side-chain χ angles.
- Nucleic acids: sugar-ring (deoxyribose/ribose) pseudorotation from ³J values.
It applies wherever a well-defined torsion connects two coupled spins. Analogous Karplus-type relations exist for ³J(C,H), ³J(H,P), ³J(C,C), and ³J(N,Cγ), each with its own coefficients, making the underlying idea one of NMR's most general structural tools.
The Derivation, Step by Step
In 1959, Martin Karplus used valence-bond theory to compute how the Fermi-contact interaction between two vicinal protons varies with rotation about the central C–C bond (J. Chem. Phys. 30, 11). The coupling is transmitted by electron spins in the bonding orbitals, and its size tracks the overlap of the two C–H bond orbitals across the intervening C–C bond.
- Step 1: Orbital overlap between the two C–H σ-bonds is maximal when the bonds are coplanar — that is, at θ = 0° (eclipsed) and θ = 180° (anti).
- Step 2: Overlap goes through a minimum near θ = 90°, where the bonds are orthogonal, so the Fermi-contact pathway is weakest.
- Step 3: A cos²θ term captures this twice-per-turn periodicity; a cosθ term makes anti (180°) larger than syn (0°); a constant C sets the baseline.
The result is ³J = A cos²θ + B cosθ + C. Karplus's original coefficients (A ≈ 8.5, B ≈ −0.28, C ≈ 0 for 0–90°) later got refined; in 1963 he stressed that A, B, C shift with substituent electronegativity, bond angles, and bond lengths.
Key Quantities and a Worked Example
A common textbook parameterization for hydrocarbons is ³J = 8.5 cos²θ − 0.28 (0°≤θ≤90°) and ³J = 9.5 cos²θ − 0.28 (90°≤θ≤180°). Ethane-derivative fits often use ³J = 12 cos²θ − cosθ + 2. The curve produces the diagnostic pattern chemists memorize:
- Anti (θ = 180°): ³J ≈ 9–12 Hz — large.
- Gauche (θ = 60°): ³J ≈ 2–4 Hz — small.
- Orthogonal (θ ≈ 90°): ³J ≈ 0 Hz — near-zero minimum.
Worked example: In a rigid trans-decalin-type ring, two adjacent protons are trans-diaxial (θ ≈ 180°). Plugging into 9.5cos²(180°) − 0.28 gives 9.5(1) − 0.28 ≈ 9.2 Hz. Measure that splitting, and you confirm both are axial. If instead you observe ~3 Hz, the geometry must be axial–equatorial (θ ≈ 60°): 8.5cos²(60°) − 0.28 = 8.5(0.25) − 0.28 ≈ 1.85 Hz. The ~5× difference is what makes ring stereochemistry so readable by NMR.
How It's Measured and Used in Practice
The input is simply the multiplet splitting in a ¹H NMR spectrum, read in Hertz (not ppm — coupling is field-independent). You measure the peak-to-peak spacing within a doublet, triplet, or more complex multiplet; that spacing is ³J.
- Extraction: for first-order multiplets, read splittings directly; for overlapping or strongly coupled systems, use spectral simulation or J-resolved 2D NMR.
- Angle read-out: invert the Karplus equation — but note it is many-to-one: a given ³J can map to two or more angles, so extra constraints (NOEs, rigid rings, additional couplings) are used to pick the right one.
- Averaging: for freely rotating bonds, the observed ³J is a population-weighted average over rotamers, letting you extract conformer populations rather than a single angle.
In protein NMR, Vuister and Bax (1993) calibrated ³J(HN–Hα) = 6.51 cos²(φ−60°) − 1.76 cos(φ−60°) + 1.60 to read backbone φ angles, a workhorse restraint in solution-structure determination.
Karplus vs. Related Relationships
The Karplus equation is one member of a family of angle–coupling relations, and it should not be confused with its cousins:
- Generalized (Haasnoot–Altona) equation: introduced in 1980, it augments the base curve with terms for substituent electronegativity and orientation, derived from 315 couplings and six parameters. It is essential for oxygen-rich systems like sugars where bare Karplus fails.
- Chemical shift (δ, in ppm): unrelated — δ reports electronic shielding, not geometry; Karplus concerns the splitting ³J in Hz.
- NOE (Nuclear Overhauser Effect): a through-space distance probe (1/r⁶), complementary to Karplus's through-bond angle information.
- Geminal (²J) and long-range (⁴J) couplings: different bond counts with their own, weaker angular dependences.
Karplus is uniquely a dihedral-angle tool: it answers "how are these two groups twisted?" where NOE answers "how far apart are they?" and δ answers "what electronic environment?"
Exceptions, Limits, and Significance
The Karplus equation is powerful but not exact. Its honest limits matter:
- Coefficient sensitivity: A, B, C are not universal — electronegative substituents (O, N, F, Cl) lower A and shift maxima, so a hydrocarbon curve mis-reads a C–O–C–H fragment by several Hz.
- Two-fold ambiguity: because cos²θ is symmetric, ³J alone often cannot distinguish, say, +60° from −60° or 30° from 150°.
- Motional averaging: flexible bonds give a Boltzmann-averaged ³J, not a snapshot; misreading an average as a fixed angle is a classic error.
- Bond-angle and electronegativity coupling: ring strain that changes H–C–C bond angles also perturbs J.
Its significance is hard to overstate: Karplus (later a 2013 Nobel laureate for multiscale modeling) gave chemists a direct ruler for molecular twist. From assigning steroid and carbohydrate stereochemistry to determining protein and RNA solution structures, the equation turns routine ¹H NMR into a conformational microscope — one of the most-cited relationships in all of chemistry.
| Dihedral angle θ | Conformation | Predicted ³J (Hz) | Where seen |
|---|---|---|---|
| 0° (syn/eclipsed) | Cis-eclipsed | ~8.5 | Rigid cyclic cis-vicinal H |
| 60° (gauche) | Synclinal | ~2–4 | Axial–equatorial, gauche rotamer |
| 90° | Orthogonal | ~0 | Locked bicyclic dihedrals |
| 120° | Anticlinal | ~2–4 | Equatorial–axial ring H |
| 180° (anti) | Antiperiplanar | ~9.5–12 | Trans-diaxial, extended chains |
Frequently asked questions
What does the Karplus equation actually tell you?
It converts a measured three-bond ³J coupling constant (in Hz) into the dihedral (torsion) angle between the two coupled nuclei. Using ³J = A cos²θ + B cosθ + C, you can infer whether two vicinal protons are anti (~180°, large J of 9–12 Hz) or gauche (~60°, small J of 2–4 Hz), revealing molecular conformation and stereochemistry directly from an NMR spectrum.
Why is the coupling large at 180° and near zero at 90°?
The coupling is transmitted by overlap of the two C–H bond orbitals across the central C–C bond via the Fermi-contact mechanism. Overlap is maximal when the bonds are coplanar (θ = 0° or 180°) and minimal when they are orthogonal (θ ≈ 90°). The cos²θ dependence in Karplus's valence-bond derivation captures exactly this geometry, driving J to a minimum near 90°.
Are the A, B, and C coefficients the same for every molecule?
No. Karplus himself warned in 1963 that A, B, and C depend on substituent electronegativity, bond angles, and bond lengths. Hydrocarbon fits use roughly A ≈ 8.5, B ≈ −0.28, C ≈ 0, but electronegative atoms like oxygen shrink the coupling and shift the curve. That is why generalized equations (Haasnoot–Altona, 1980) add explicit electronegativity terms for accurate work on sugars and other heteroatom-rich systems.
What is the difference between the original Karplus and the Haasnoot–Altona equation?
The original Karplus equation uses only the torsion angle. The Haasnoot–Altona generalized equation (1980) adds correction terms for the number, electronegativity, and orientation of substituents on the H–C–C–H fragment, parameterized from 315 experimental couplings with six fitted constants. It gives far better accuracy for conformational analysis of nucleosides, nucleotides, and carbohydrates.
Can one ³J value give more than one angle?
Yes — the relationship is many-to-one. Because cos²θ is symmetric, a single ³J can correspond to two or more dihedral angles (for example, a 4 Hz coupling might fit either ~40° or ~140°). Chemists resolve the ambiguity using additional couplings, NOE distance data, rigid ring constraints, or knowledge that a system is conformationally locked.
Is ³J measured in ppm or Hz, and does field strength matter?
Coupling constants are always measured in Hertz, taken as the spacing between lines within a multiplet. Unlike chemical shift (ppm), J is independent of the spectrometer's magnetic field strength — a 7 Hz coupling is 7 Hz at 300 MHz or 900 MHz. This field-independence is what makes ³J a reliable, transferable geometric probe.