Skip to main content

Surface Properties

de and di

ded_e is the distance from the Hirshfeld surface to the nearest nucleus outside the surface. did_i is the corresponding distance to the nearest nucleus inside the surface.

ded_e and did_i were the original surface properties mapped onto the Hirshfeld surface. Normally, we only map ded_e, since this property provides information about the crystal packing environment around the surface.

The simplest and most immediately useful property to map onto the surface is the distance from the surface to the nearest nucleus external to the surface, which we call ded_e. This property provides an immediate picture of the nature of intermolecular contacts in the crystal.

In a similar fashion, we can define did_i, the distance from the surface to the nearest nucleus internal to the surface. We have not yet utilised this property as an independent property mapped on the surface, however it is valuable when used in conjunction with ded_e.

The range of ded_e and did_i across the Hirshfeld surface varies considerably depending on the atoms in the molecule (size dependence) and the particular type of intermolecular interaction experienced (interaction dependence). Mapping de over the same range of all molecules would reduce the colour contrast for molecules which feature only a small range of contact distances (such as hydrocarbons). To gain the maximum benefit from the colour mapping on the surface, we choose a ded_e range most suited to the each group of molecules under direct comparison.

dnorm

dnormd_\text{norm} is a normalised contact distance. did_i is normalised by the van der Waals radius of the atom involved; ded_e is similarly normalised, and the sum of these two normalised quantities is the dnormd_\text{norm} property:

dnorm=de+did_{\text{norm}} = d_{|e|} + d_{|i|}

where

de=dervdwrvdwd_{|e|} = \frac{d_e - r_{\text{vdw}}}{r_{\text{vdw}}}
di=dirvdwrvdwd_{|i|} = \frac{d_i - r_{\text{vdw}}}{r_{\text{vdw}}}\\

Where atoms make intermolecular contacts closer than the sum of their van der Waals radii, these contacts will be highlighted in red on the dnormd_\text{norm} surface. Longer contacts are blue, and contacts around the sum of van der Waals radii are white, as shown in the diagram on the right.

Experimental features

We understand that the van der Waals radius is an arbitrary value for choosing to highlight intermolecular contacts (it is a sensible arbitrary value, though). So, in the case that you might want to highlight contacts that are somewhat greater than the this distance, we have been experimenting with two alternative methods:

  1. Secondary Highlight: Change the colour scheme so that contacts at the van der Waals separation are yellow, and contacts 5% greater than the van der Waals separation are white.
  2. Sensitivity:Drag the slider bar to change the white point of the dnormd_\text{norm} property between 0 and 10% greater than van der Waals separation. This should not be used for publication purposes, as there is no way to reproduce the exact setting you use.

Curvedness and Shape Index

A molecular surface defines the shape of a molecule. Because the Hirshfeld surface defines the shape of the molecule in terms of its surrounding crystalline environment, we think that the local shape of the surface may provide some chemical insight.

At any point pp on the surface, we can determine the outward normal, and there exist two principal directions u\mathbf{u} and v\mathbf{v} (see right) along which the principal curvatures κ1\kappa_1 and κ2\kappa_2 are calculated

The two conventional measures of curvature, the mean curvature, HH, and the gaussian curvature, KK, do not provide much physical information. Koenderink12 has introduced two more useful measures of surface curvature, the curvedness, CC, and the shape index, SS as follows:

C=κ12+κ222C = \sqrt{\frac{\kappa_1^2 + \kappa_2^2}{2}}
S=2πarctan(κ2+κ1κ2κ1)S = \frac{2}{\pi} \arctan(\frac{\kappa_2 + \kappa_1}{\kappa_2 - \kappa_1})

with κ2κ1\kappa_2 \ge \kappa_1.

The curvedness is a function of the root-mean-square curvature of the surface, with flat areas of the surface having a low curvedness and areas of sharp curvature having a high curvedness. Areas on the Hirshfeld surface with high curvedness tend to divide the surface into contact patches with each neighbouring molecule, so that the curvedness of the Hirshfeld surface could be used to define a coordination number in the crystal.

The shape index is a qualitative measure of shape and can be sensitive to very subtle changes in surface shape, particularly in regions where the total curvature (or the curvedness) is very low. One important attribute of the shape index is that two shapes where the shape index differs only by a sign represent complementary "stamp" and "mould" pairs. This means that maps of shape index on the Hirshfeld surface can be used to identify complementary hollows (with S<1S \lt 1) and bumps (with S>1S \gt 1).


  1. J.J. Koenderink, Solid Shape,Cambridge MA, MIT Press, 1990
  2. J.J. Koenderink, A.J. Van Doorn, Image and Vision Computing, 1992, 10, 557.