Interatomic Distances and Angles
This sheet: https://xraytools.com/geometry?contents=cell
Page: https://xraytools.com/geometry
The first thing anybody asks of a solved structure: how far apart are the atoms, and at what angles. Give a cell, a space group and the asymmetric unit — the page expands them by the symmetry, finds every neighbour inside a radius you choose, and names each one by the operation that produced it.
- Before this
- The atoms you supply are the asymmetric unit; the neighbours a distance is measured to are mostly symmetry copies the page generates. A distance to an atom you never typed is the normal case, not an error.
- You supply
- A unit cell, a space group and one atom per line of the asymmetric unit. Any constant or coordinate may carry its standard uncertainty in brackets, as a CIF writes it, and those are what the uncertainties below are propagated from.
- Reading it
- Uncertainties are propagated as if the parameters were uncorrelated, because a CIF does not carry the refinement’s variance–covariance matrix. That is not what a refinement program prints: for two atoms of one rigid group the correlation is usually positive, which makes the figure here larger. It is not an upper limit either.
Worked examples: Quartz, with uncertainties supplied · cubic – NaCl, Fm3m · cubic – Cu, Fm3m · cubic – α-Fe, Im3m · cubic – CsCl, Pm3m · cubic – ZnS, F43m · hexagonal – quartz, P3221 · tetragonal – cristobalite, P41212 · hexagonal – berlinite, P3121 · tetragonal – rutile, P42/mnm · orthorhombic – aragonite, Pmcn · monoclinic – ZrO2, P21/c · triclinic – albite, C1 · tetragonal – urea, P421m
Earlier on the path: CIF Parser Next on the path: Displacement Parameters and NPD Atoms On How to read a published structure, step 2 of 4
Notation here: U, B — what each one means here
Terms here: asymmetric unit · setting · zone
See also: CIF Parser · Displacement Parameters and NPD Atoms
What each input changes
- How far to look
- How far out neighbours are looked for. It decides which contacts are listed and nothing about their values; a bond does not become a bond by raising it.
- Atoms in the asymmetric unit
- The asymmetric unit. Most of the neighbours in the answer are symmetry copies the page generates, so a distance to an atom you never typed is the normal case.
Teaching with this page
- Objective
- After this page a learner can decide whether two reported distances are significantly different.
- Start from
- this worked example
- Ask first
- Two bonds are reported as 1.943(4) Å and 1.947(4) Å. Are they different?
- Watch for
- “No — the difference is inside the uncertainty, so they are equal”
- Then
- Displacement Parameters and NPD Atoms
Check yourself: Two bonds are reported as 1.943(4) Å and 1.947(4) Å. Are they different?
Not resolved, which is not the same as being equal No — the difference is inside the uncertainty, so they are equal
The difference is 0.004 Å against a combined uncertainty near 0.006 — well under the 3 s.u. that is conventionally taken as the threshold for calling a difference real — so the data cannot separate them. That is a statement about the data and not about the bonds: failing to resolve a difference is not evidence that there is none, and better data might resolve it. The combined figure is approximate as well — this page propagates as if the refined parameters were uncorrelated, and a CIF carries no variance–covariance matrix to do better.
These are distances and angles, not bonds. Every neighbour inside the radius is listed, whether or not anything holds the two atoms together, and a short contact between ions of the same charge is a repulsion rather than a bond. Contacts are found through every symmetry operation of the space group and across cell boundaries, so the list is the full environment and not only what lies inside one cell. Which of these contacts are bonds is asked separately below, by sorting them on how far each one exceeds the two covalent radii and cutting at the largest gap in that order — and it is that answer, not this table, that the picture draws.
The cell, filled
Si Z 14 · rcov 1.11 Å O Z 8 · rcov 0.66 Å
Hover an atom to name it. Click one to pin it and mark every copy of that site.
Every atom of the cell is drawn, including the symmetry-equivalent ones. An atom lying on a face, an edge or a corner belongs to each cell it touches and is drawn in all of them, so counting the spheres overcounts the contents of one cell. Colours are the standard CPK ones, and spheres are drawn at half the covalent radius (Cordero et al., 2008), reduced further only if that would make two of them run into each other — at full size two bonded atoms touch by definition. These are not ionic radii: in a salt the cation is drawn larger than the anion, which is the opposite of the ionic picture. A stick is drawn where both atoms' own bond ladders put the other below their cut, so every bond here is one the table above lists — and where the two ladders disagree the table is the fuller answer, not this picture. A sphere outside the cell is drawn only because a bond from inside reaches it, and its own bonds are not completed in turn.
The symmetry the atoms obey
The atoms obey exactly the 6 operations of the space group given, and no others that this cell would permit.
Distances
| atom | neighbour | distance | symmetry of the neighbour |
|---|---|---|---|
| Si | O | 1.6039 | x,y,z |
| O | 1.6039 | x-y,-y,-z+1/3 | |
| O | 1.6132 | -y+1,x-y,z-1/3 | |
| O | 1.6132 | -x+1,-x+y,-z+2/3 | |
| Si | 3.0575 | -y,x-y-1,z-1/3 | |
| Si | 3.0575 | -y+1,x-y,z-1/3 | |
| Si | 3.0575 | -x+y+1,-x,z+1/3 | |
| Si | 3.0575 | -x+y+1,-x+1,z+1/3 | |
| O | Si | 1.6039 | x,y,z |
| Si | 1.6132 | -x+y+1,-x+1,z+1/3 | |
| O | 2.6109 | x-y,-y,-z+1/3 | |
| O | 2.6157 | -x+1,-x+y,-z+2/3 | |
| O | 2.6157 | -x+1,-x+y+1,-z+2/3 | |
| O | 2.6306 | y,x,-z+1 | |
| O | 2.6438 | -y+1,x-y,z-1/3 | |
| O | 2.6438 | -x+y+1,-x+1,z+1/3 |
Which of these are bonds
| atom | neighbour | distance | radius sum | slack | gap |
|---|---|---|---|---|---|
| Si | 2 × O | 1.6039 | 1.7700 | −0.1661 | — |
| 2 × O | 1.6132 | 1.7700 | −0.1568 | 0.0093 | |
| 4 × Si | 3.0575 | 2.2200 | 0.8375 | 0.9943 | |
| Si: 4 neighbours below the cut — 4 × O. The gap is 0.9943 Å, 106.87 times the next largest. | |||||
| O | Si | 1.6039 | 1.7700 | −0.1661 | — |
| Si | 1.6132 | 1.7700 | −0.1568 | 0.0093 | |
| O | 2.6109 | 1.3200 | 1.2909 | 1.4477 | |
| 2 × O | 2.6157 | 1.3200 | 1.2957 | 0.0047 | |
| O | 2.6306 | 1.3200 | 1.3106 | 0.0150 | |
| O: 2 neighbours below the cut — 2 × Si. The gap is 1.4477 Å, 96.83 times the next largest; 2 further contacts inside this radius are not shown. | |||||
The cut is where the largest gap falls, not a criterion anyone chose. Slack is the contact's length less the sum of the two covalent radii, so a bond has little of it and a passing neighbour has a great deal; sorting by slack puts the coordination shell at the top whatever elements it is made of. Nothing here decides that a contact is a bond — it shows you where the evidence changes, and you decide.
A covalent radius is the wrong radius for a large ion, and that is where this fails. Where a published coordination number exists to compare against, the cut reproduces it for 26 of 35 atoms at a 5 Å radius, and 5 of the 9 misses have a soft cation — rock salt's sodium, caesium chloride's caesium, rutile's titanium, aragonite's calcium, baddeleyite's zirconium. The radius decides as much as the chemistry: the same 35 atoms give 23 right at 3.2 Å and 25 at 4 Å. A wrong cut usually has a gap barely bigger than the next one, so the ratio beside each cut is worth reading — but it is an association and not a test, and albite's sodium breaks it at 3.2 Å with a gap 3.4 times the next and the wrong answer.
Angles
| at | between | angle |
|---|---|---|
| Si | O and O | 108.962 |
| O and O | 110.529 | |
| O and O | 108.788 | |
| O and Si | 124.620 | |
| O and Si | 92.837 | |
| O and Si | 120.144 | |
| O and Si | 18.181 | |
| O and O | 108.788 | |
| O and O | 110.529 | |
| O and Si | 18.181 | |
| O and Si | 120.144 | |
| O and Si | 92.837 | |
| O and Si | 124.620 | |
| O and O | 109.241 | |
| O and Si | 93.021 | |
| O and Si | 18.072 | |
| O and Si | 113.632 | |
| O and Si | 109.144 | |
| O and Si | 109.144 | |
| O and Si | 113.632 | |
| O and Si | 18.072 | |
| O and Si | 93.021 | |
| Si and Si | 106.933 | |
| Si and Si | 91.197 | |
| Si and Si | 141.618 | |
| Si and Si | 123.308 | |
| Si and Si | 91.197 | |
| Si and Si | 106.933 | |
| O | Si and Si | 143.747 |
| Si and O | 35.519 | |
| Si and O | 35.724 | |
| Si and O | 152.341 | |
| Si and O | 139.454 | |
| Si and O | 34.850 | |
| Si and O | 109.593 | |
| Si and O | 163.412 | |
| Si and O | 111.015 | |
| Si and O | 35.487 | |
| Si and O | 35.380 | |
| Si and O | 130.975 | |
| Si and O | 34.621 | |
| O and O | 60.774 | |
| O and O | 157.969 | |
| O and O | 129.824 | |
| O and O | 59.702 | |
| O and O | 141.587 | |
| O and O | 139.844 | |
| O and O | 106.226 | |
| O and O | 60.019 | |
| O and O | 80.834 | |
| O and O | 60.521 | |
| O and O | 118.563 | |
| O and O | 59.524 | |
| O and O | 159.350 | |
| O and O | 59.460 | |
| O and O | 101.341 |
Every pair of neighbours inside the radius is here, which is more than the bond angles: two atoms that are not bonded to each other still subtend an angle at the middle one, and a small value usually means the two legs are of very different length. The distance table above is what says which neighbours are close enough to be bonds.
Torsions
| about | chain | torsion | sector | chains |
|---|---|---|---|---|
| Si-O | O–Si–O–Si | +−107.400 | −ac | 3 |
| O–Si–O–Si | +−88.973 | −sc | 3 | |
| O–Si–O–Si | +12.267 | sp | 3 | |
| O–Si–O–Si | +30.966 | +sc | 3 | |
| O–Si–O–Si | +131.865 | +ac | 3 | |
| O–Si–O–Si | +151.535 | ap | 3 |
A torsion angle is measured looking along the middle bond: it is the angle from the first atom to the fourth, projected onto the plane across that line. It is the first quantity on this page whose sign carries information no distance or angle does — mirror a crystal and every distance and every angle is unchanged, while every torsion changes sign. Quartz is the example: its two enantiomorphs give an identical set of bond lengths and the exactly opposite set of torsions. Two values are their own opposite and so carry no sign here: 0, where the chain is eclipsed, and 180, where it is anti. Each of those is superimposable on its own mirror image, so there is no hand to report.
The sign follows the convention of Klyne and Prelog, which IUPAC adopted as the standard for describing conformation: look from the first atom along the middle bond towards the fourth, and the torsion is positive when the near bond has to turn clockwise, through less than 180°, to eclipse the far one. Which end you look from does not change the answer — reading the chain backwards gives the same number, sign included, which is why a chain and its reverse are one row above. The form of it a reader can check against a picture rather than against algebra: a right-handed helix has positive torsions.
What this page computes, with b1 = B−A, b2 = C−B and b3 = D−C:
τ = atan2( |b2| b1·(b2×b3), (b1×b2)·(b2×b3) )
Two arguments rather than one, which is the whole point: an arc cosine of the angle between the two planes gives the size and throws the sign away, and a sign applied afterwards is a second convention to get wrong. Here it comes out of the arithmetic.
The sector column names the range, in IUPAC's terms:
- sp synperiplanar, 0–30°
- ±sc synclinal, 30–90°
- ±ac anticlinal, 90–150°
- ap antiperiplanar, 150–180°
IUPAC gives those ranges with shared endpoints — 0 to ±30 synperiplanar, 30 to 90 synclinal — so a torsion of exactly 30° is in two of them and the recommendation does not settle which. This page gives a boundary to the sector nearer zero; that is a choice, not a standard. The sign goes on sc and ac because +sc and −sc are two different sectors, while sp is one sector straddling zero and ap one straddling 180° — a sign on those would only repeat the number beside it.
Each of the three bonds in a chain comes from the ladder above, so every reservation there applies here three times over. Where a cut runs through a large soft cation the chains built on it are numerous and the least trustworthy on the page: read the ladder before the table.
Hydrogen bonds
This structure has no hydrogen in it, so it has no hydrogen bonds. Every other worked example below is an inorganic solid without hydrogen; load urea to see this table with something in it.
A hydrogen bond D–H···A is a hydrogen held between the atom it is covalently bonded to and a second one it is not. This table takes every hydrogen, finds its covalent bond from the ladder above, and keeps the contacts with the hydrogen genuinely between the two — which is the D–H···A angle exceeding 90°, and needs no cutoff to say. What survives is ranked by the same ladder, so the cut here and the cut above are one rule. One thing here is a convention rather than a measurement, and it is the elements: the donor must be N, O, F or S and the acceptor one of those or a halogen, because a hydrogen bond needs a polarised bond at one end and a lone pair at the other and neither is visible in a list of coordinates. That excludes C–H···O, which is a real if weaker interaction — the contact is still in the distance table above, it is only the name that is withheld. Urea is the example the page ships: each of its four N–H donates to a carbonyl oxygen, and each oxygen accepts four.
Every row rests on the ladder twice — once for the D–H bond it starts from and once for the cut among what is left — so each reservation there applies here twice. And the shortest contact of a hydrogen is often not its hydrogen bond: in urea the two closest neighbours of H1 are the other hydrogen on the same nitrogen and the carbon two bonds away, both of which sit beside it rather than in front of it.
Packing and voids
| probe radius / Å | void fraction | void volume / Å3 |
|---|---|---|
| 0.0 | 0.0958 | 10.8 |
| 0.2 | 0.0186 | 2.1 |
| 0.4 | 0.0000 | 0.0 |
| 0.6 | 0.0000 | 0.0 |
| 0.8 | 0.0000 | 0.0 |
| 1.2 | 0.0000 | 0.0 |
Packing fraction 0.9042 — that fraction of the cell is inside an atom.
The packing fraction is the row at a probe radius of zero: how much of the cell lies inside an atom. The rows below it ask a different question — how much space is left for a sphere of that radius to sit in without overlapping anything, which is what a crystallographer means by a void. A water molecule is usually given 1.2 Å.
Measured by sampling 110,592 points on a grid of 483, offset from the cell origin by an irrational fraction. The offset is not cosmetic: a grid that lines up with the lattice samples whole planes of points onto sphere boundaries, and its error then depends on the arithmetic relationship between the grid and the cell rather than on the resolution — refining it does not help. Against the four lattices whose packing fraction is an exact constant, this grid is right to about 5×10−4, so the fourth decimal above is the last one worth reading.
Radii are van der Waals radii from Alvarez (2013), the revision of Bondi's set derived from the Cambridge Structural Database. That choice is a convention and it is the only one on this table — the volumes themselves are geometry. A van der Waals radius describes how close a non-bonded neighbour comes, so in a metal or an ionic solid, where every contact is a bond, the spheres overlap and the packing fraction is 1 by construction. The number means what it says for molecular crystals.
About the uncertainties
No uncertainties are quoted, because none were given. Nothing in the cell or the coordinates above carries a bracket, so there is nothing to propagate. An absent uncertainty is not a zero one.
Write a constant as 4.9137(2) or a coordinate as
0.4697(1) and every distance and angle below will carry its own, split into the
part that comes from the cell and the part that comes from the coordinates.
Contacts were searched out to 3.2 Å. A neighbour is named by the operation that produces it from the atom in the list above, written out in full rather than as a numbered code.
Where this comes from
- Computer derivation of the symmetry elements implied in a structure description
Y. Le Page, J. Appl. Cryst. 1987, 20, 264–269 · doi:10.1107/S0021889887086710
The method behind the additional-symmetry panel: search the metric symmetry of the lattice, then test each candidate against the atoms. - MISSYM 1.1 — a flexible new release
Y. Le Page, J. Appl. Cryst. 1988, 21, 983–984 · doi:10.1107/S0021889888007022
The follow-up, and the source of the tolerance handling: how far an atom may sit from its image and still count. - Covalent radii revisited
B. Cordero et al., Dalton Trans. 2008, 2832 · doi:10.1039/b801115j
The radii the distance table draws its spheres from. It publishes standard deviations with them, which is why they can carry a stated tolerance rather than a hidden one. - Description of steric relationships across single bonds
W. Klyne, V. Prelog, Experientia 1960, 16, 521–523 · doi:10.1007/BF02158433
The sign convention the torsion column follows, and the source of the synperiplanar / synclinal / anticlinal / antiperiplanar names. A torsion is the one quantity on this page whose sign carries the handedness of the structure, so which of the two possible signs is printed is not a detail: mirror a crystal and every distance and every angle is unchanged while every torsion changes sign. - Basic terminology of stereochemistry (IUPAC Recommendations 1996)
G. P. Moss, Pure Appl. Chem. 1996, 68, 2193–2222 · doi:10.1351/pac199668122193
The recommendation that adopted the convention above and fixed the six sectors the table names. It gives their ranges with SHARED endpoints — 0 to ±30 synperiplanar, 30 to 90 synclinal — so a torsion of exactly 30 degrees falls in two of them and the recommendation does not settle which. This page gives a boundary to the sector nearer zero and says so, because that is a choice rather than a standard.