Disorder in crystal structure¶
Introduction¶
The class ccdc.crystal.Crystal.Disorder
describes the disorder in a
crystal structure. We can investigate the disorder and make selection of the
available disorder groups.
Loading a structure¶
>>> from functools import reduce
>>> import operator
>>> from ccdc import crystal
>>> from ccdc import io
Here we load a cif file for a crystal structure with multiple disorder assemblies.
>>> cry_reader = io.CrystalReader(disorder_crystal_cif_file)
>>> cry = cry_reader[0]
>>> cry.has_disorder
True
>>> cry.disorder.is_suppressed
False
Investigate crystal disorder¶
The disorder in this structure is not suppressed, which means the disorder is fully represented, and we are able to work with individual disorder groups.
>>> len(cry.disorder.assemblies)
3
>>> [len(assembly.groups) for assembly in cry.disorder.assemblies]
[2, 2, 2]
There are three disorder assemblies, each with two disorder groups. By default, the first disorder group of each disorder assembly is active. Let’s select the second disorder group of the second disorder assembly as active. The void volume of the crystal structure changes as a result.
>>> [assembly.active.id for assembly in cry.disorder.assemblies]
[0, 0, 0]
>>> round(cry.void_volume(probe_radius=0.6, grid_spacing=0.3), 2)
23.88
>>> cry.disorder.assemblies[1].groups[1].activate()
>>> [assembly.active.id for assembly in cry.disorder.assemblies]
[0, 1, 0]
>>> round(cry.void_volume(probe_radius=0.6, grid_spacing=0.3), 2)
24.23
Going through all combinations of the disorder¶
The ccdc.crystal.Crystal.Disorder.combinations
attribute returns an
iterator that yields all combinations of the disorder groups. The crystal is
updated accordingly with the yielded disorder groups. This provides a
convenient way to calculate crystal property taking into account the disorder
occupancies.
>>> void_vol = 0
>>> for active_groups in cry.disorder.combinations:
... vol = cry.void_volume(probe_radius=0.6, grid_spacing=0.3)
... occ = reduce(operator.mul, [group.occupancy for group in active_groups], 1)
... void_vol += vol * occ
>>> print(f'The disorder averaged void volume is {void_vol:.2f} cubic angstroms')
The disorder averaged void volume is 24.04 cubic angstroms