We report a Monte Carlo algorithm
for calculation of occupied (“network”)
and unoccupied (“void”) space in crystal structures.
The variation of the volumes of the voids and the network of intermolecular
contacts with pressure sensitively reveals discontinuities associated
with first- and second-order phase transitions, providing insights
into the effect of compression (and, in principle, other external
stimuli) at a level between those observed in individual contact distances
and the overall unit cell dimensions. The method is shown to be especially
useful for the correlation of high-pressure crystallographic and spectroscopic
data, illustrated for naphthalene, where a phase transition previously
detected by vibrational spectroscopy, and debated in the literature
for over 80 years, has been revealed unambiguously in crystallographic
data for the first time. Premonitory behavior before a phase transition
and crystal collapse at the end of a compression series has also been
detected. The network and void volumes for 129 high-pressure studies
taken from the Cambridge Structural Database (CSD) were fitted to
equation of state to show that networks typically have bulk moduli
between 40 and 150 GPa, while those of voids fall into a much smaller
range, 2–5 GPa. These figures are shown to reproduce the narrow
range of overall bulk moduli of molecular solids (
ca.
5–20 GPa). The program, called CellVol, has been written
in Python using the CSD Python API and can be run through the command
line or through the Cambridge Crystallographic Data Centre’s
Mercury interface.