Choosing Molecular Coordinates

Introduction

A potential energy surface is independent of how we label a molecular geometry, but our calculations and explanations are not. Cartesian positions are universal and convenient for forces. Bond lengths, angles and torsions match chemical intuition. Collective coordinates compress many atomic movements into a progress variable. Choosing well can reveal a pathway; choosing badly can hide an intermediate or make a barrier appear at the wrong place.

Core explanation

For N atoms, Cartesian coordinates list 3N numbers: x, y and z for each nucleus. They are easy to use in gradients and molecular dynamics, but six combinations for a nonlinear free molecule represent overall translation and rotation rather than internal chemistry. Once these motions are removed, 3N−6 internal degrees of freedom remain. Cartesian displacements can be hard to interpret chemically because a bond stretch and rotation of the whole molecule may be mixed in a raw vector.

Internal coordinates use relative properties. Bond distances describe separation; bond angles describe bending; dihedral angles describe torsion around linked bonds. A redundant internal-coordinate set may include more entries than independent degrees of freedom yet make optimisation robust for a particular molecule. Internal coordinates can fail or become numerically awkward near a linear angle, bond breaking or a topology change. When a bond disappears, the originally chosen set may no longer describe the product naturally.

A collective coordinate is a function of many positions. For proton transfer from donor D to acceptor A, q = r(D–H) − r(A–H) can track which side holds the proton. Negative and positive q may identify donor- and acceptor-bound regions. Yet the donor–acceptor distance, solvent arrangement and other motions also affect the barrier. If those variables are fixed, a scan along q may overestimate or miss the real bottleneck. A good coordinate distinguishes states while allowing orthogonal degrees of freedom to relax or be sampled.

Mass weighting scales displacements by square roots of nuclear masses. It is useful for normal modes and for the conventional intrinsic reaction coordinate. In ordinary Cartesian space, equal geometric displacements of a hydrogen and a heavy atom have different kinetic implications; mass weighting reflects that distinction. A path's geometric shape can therefore depend on the metric used, even though the underlying stationary geometries of a fixed electronic PES do not change when isotope masses change in the Born–Oppenheimer approximation.

Coordinate choice affects visualisation and algorithms, but not the physical energy assigned to the same exact geometry. A reaction path drawn in one projection may appear to cross itself even when full-dimensional geometries differ. Likewise, a “flat” curve in one coordinate may conceal a steep orthogonal ridge. Always state what was held fixed, relaxed or averaged when showing a profile.

Step-by-step reasoning

Identify the chemically changing bonds and likely large-scale motions. Choose Cartesian coordinates for full dynamics, internal coordinates for chemically interpretable structural changes, or a collective coordinate for a targeted projection. Check coordinate independence and singularities. For a scan, decide which other variables relax. Compare candidate coordinates against configurations from both reactants and products, and test whether states overlap in the projected coordinate.

Visual explanation

Draw a three-atom D–H···A system with labelled Cartesian axes, two bond lengths and their difference q. On a second panel show that two configurations can share q=0 while having different D–A distances and energies. Draw a contour map in q and D–A distance, then project it onto q alone to illustrate information loss.

Real-world analogy

A journey can be described by latitude and longitude, road distance, or percentage of route completed. Each coordinate is useful for a different question; “50% completed” does not reveal whether the traveller is on a mountain or in a tunnel. A molecular collective coordinate likewise compresses information, but nuclear motion occurs in a much higher-dimensional space.

Real-world example

For a torsional conformational change, a dihedral angle can be a natural coordinate. A scan rotates the bond and relaxes other geometric variables, producing an energy profile with minima and a barrier. If the molecule also folds through a coupled ring motion, the same dihedral value may correspond to multiple distinct conformers. A two-dimensional scan or an unbiased sampling method can reveal the hidden coordinate.

Why?

Why not use Cartesian coordinates for every plot? They are high-dimensional and less chemically legible. Why not use only one bond length? Reactions often involve simultaneous bond making, breaking, solvent rearrangement or orientation change. Why mass-weight an IRC? It sets a metric connected to nuclear kinetic energy and vibrational modes rather than treating equal geometric displacements of all nuclei identically.

Common misconception

A reaction coordinate is not necessarily a literal distance, and a monotonic bond change is not proof of a complete reaction pathway. Another mistake is to assume stationary-point geometries differ by isotope solely because their mass-weighted paths differ; on an unchanged Born–Oppenheimer electronic surface, stationary geometries depend on potential energy, not nuclear mass, apart from corrections beyond that approximation.

Worked example

Question: For a proton between D and A, calculate q = r(D–H) − r(A–H) in two structures: (1.1 Å, 1.8 Å) and (1.5 Å, 1.5 Å). Interpret each value cautiously.

Reasoning: The first gives q = 1.1−1.8 = −0.7 Å, indicating H is closer to D. The second gives q=0, equal stated distances. The second could be near a transfer bottleneck, but q alone does not reveal D–A separation, solvent geometry or whether it is a saddle.

Answer: q is −0.7 Å and 0 Å; equal bond distances alone do not prove a transition state.

Quick check

1. What information does a one-dimensional proton-transfer coordinate omit if donor–acceptor distance is not included? Answer: It omits how far the donor and acceptor are apart and any coupled geometric or solvent rearrangement.

Exam focus

State coordinate definitions and units. Distinguish changing the labels of a geometry from changing its physical energy. For a constrained scan, report which degrees of freedom remain relaxed. Explain why two structures with the same projected coordinate can have different energies.

Advanced insight

Reaction coordinates can be tested dynamically through the committor, the probability that a configuration with sampled momenta reaches product before reactant. A good one-dimensional coordinate tends to separate reactant and product basins and locate a narrow bottleneck ensemble, but no simple coordinate is guaranteed. In solution, a solvent-polarisation coordinate can be more important than a solute bond length. Modern enhanced sampling often uses several collective variables to avoid concealing slow orthogonal motion.

Summary

Cartesian, internal and collective coordinates describe the same geometries from different viewpoints. Cartesian coordinates support general forces; internals express bonds and torsions; collective variables focus on a process but lose information. Mass weighting changes the path metric. A reliable profile states what was fixed, relaxed or averaged and tests whether its coordinates capture the actual bottleneck.

Practice questions

1. How many Cartesian nuclear coordinates are there for six atoms? Answer: 3×6 = 18 before removing overall translation and rotation.

2. Give one reason an internal-coordinate set may become awkward during bond breaking. Answer: A bond-length or angle definition tied to the reactant topology may become singular or chemically inappropriate as bonds disappear or rearrange.

3. Can two full-dimensional geometries have the same one-dimensional q value? Answer: Yes. They may differ in all other bond lengths, torsions or solvent positions and have different energies.

4. What is changed by mass weighting a path coordinate? Answer: The metric used to measure nuclear displacements, reflecting mass in vibrational or path calculations.

Sources: IUPAC, reaction coordinate; IUPAC, intrinsic reaction coordinate; Reaction dynamics and PES construction, Journal of Physical Chemistry A.