InternationalDefinition and exchange of crystallographic dataTables for Crystallography Volume G Edited by S. R. Hall and B. McMahon © International Union of Crystallography 2006 |
International Tables for Crystallography (2006). Vol. G, ch. 3.5, pp. 141-143
https://doi.org/10.1107/97809553602060000737 ## Chapter 3.5. Classification and use of electron density data
The electron density CIF dictionary (rhoCIF) provides terms suitable for a representation in software of the electron density distributions in crystals. It has been developed under the aegis of the IUCr Commission on Charge, Spin and Momentum Densities. The contents of the dictionary are adequate not only to record the results of experimental measurements, but also to represent a theoretical electron density that has been fitted with an appropriate model. |

Modern high-resolution X-ray diffraction methods enable the electron density distributions in many crystals to be accurately determined, allowing one to explore the interactions that hold crystals and molecules together (Coppens, 1997; Koritsanszky & Coppens, 2001). For example, Bader's (1990) theory of atoms in molecules gives a description of chemical bonding in terms of the topology of the electron density, but it requires an accurate knowledge of the electron density. [Note that the terms `electron density' (or distribution) and `charge density' tend to be used interchangeably in the literature. The term `electron density' is used here since `charge density' strictly includes both positive (nuclear) and negative (electronic) charge. Properties such as electrostatic potential, electric fields and their gradients, and molecular electric moments are derived from the charge density distribution, which includes the contribution of the nuclear charge (Spackman & Brown, 1994).] The advent of area detectors, more powerful laboratory computing facilities and developments in software have combined to produce a rapid increase in the number of accurate electron density determinations now available. While such studies are not yet routine, they are sufficiently common to justify defining a structure for a secure and systematic archive. This in turn requires a method of representing the electron density acceptable to the crystallographic community.

The CIF standard which has been adopted by the IUCr and is in common use in macromolecular and small-molecule crystallography is the obvious choice for such a file structure. At the Sagamore XIII conference held at Stare Jablonki, Poland, in 2000, the IUCr Commission on Charge, Spin and Momentum Densities decided to prepare a CIF dictionary for reporting electron densities in crystals. The result is the rhoCIF dictionary described in this chapter. It was prepared in conjunction with members of the Commission and other members of the community. [The name rhoCIF for this dictionary was chosen over the more obvious choice of edCIF to avoid any possible confusion with electron diffraction. The symbol conventionally represents the electron density in real space at the point defined by the position vector **r**.]

The distribution of electron density in a molecule or crystal can also be calculated theoretically, and the results are usually compared with experiment. The Hohenberg–Kohn theorem (Hohenberg & Kohn, 1964) shows that the ground-state energy is a unique functional of the electron density, thus providing a theoretical justification for relating chemical and physical properties directly to the observable electron density. Density functional theory has been applied to chemical problems for some time (see, for example, Labanowski & Andzelm, 1991; Ziegler, 1991; Kohn *et al.*, 1996; Nagy, 1998). Although the current version of the rhoCIF dictionary is intended to record the results of experimental measurements, it can also be used to represent a theoretical electron density that has been fitted with an appropriate model.

The rhoCIF dictionary contains the definitions of data items used to report accurate electron densities in crystals. The structure of these definitions conforms to the CIF rules and conventions, and selected data items reflect those currently needed in the field. In any field of science, the nature of the data used to describe the experiment and the results changes with time, and accurate electron density studies are no exception. The items defined in this dictionary describe the density representation approach used most often today – the atom-centred multipole expansion model introduced by Stewart (1973) and refined by Hansen & Coppens (1978). In this model, the electron density is expressed as the sum of three components. Two are the partial contributions of the core and valence shell of the spherical free atom, usually taken from the paper by Clementi & Roetti (1974). The third component is an aspherical multipole expansion of the remaining electron density. This is an atom-centred function, flexible enough to describe most electron density distributions likely to be encountered in a crystal, and with components that have the same symmetries as the atom-centred wave functions (orbitals) frequently used by chemists when describing chemical bonding.

The rhoCIF dictionary is intended to be used in combination with the core CIF dictionary, or, after adaptation to DDL2, the mmCIF dictionary. The two dictionaries can be combined to form a composite dictionary by concatenation. Alternatively, the merging procedures described in Section 3.1.9 can be used to create a virtual dictionary at run time. Because the dictionary is not self-contained, its structure can be kept very simple, comprising only two categories: one to define a local axis system around each atom and one to contain the various parameters of the Hansen & Coppens (1978) multipole function.

The two categories in the electron density CIF dictionary lie within the single ATOM category group and relate to the structural model. They are as follows:

Data items in this category are as follows:

ATOM_LOCAL_AXES
*The bullet (•) indicates a category key. The arrow is a reference to a parent data item.*

The ATOM_LOCAL_AXES category is used to define a local Cartesian axis system for each atom. Data items in this category are looped, with the category key ( ** _atom_local_axes_atom_label**) being a pointer to the

**defined in the ATOM_SITE category of the core CIF dictionary. The orientation of the Cartesian basis is defined in terms of vectors between atoms, and any arbitrary orientation can be described by including dummy atoms if necessary. The spherical coordinate system used in the description of the multipoles is defined in terms of the local Cartesian basis:**

`_atom_site_label`*z*is the polar axis from which the angle is measured, is measured from the

*x*axis in the

*xy*plane with the

*y*axis having a value of .

Example 3.5.3.1 shows how the category is used in part of an electron density CIF relating to a charge density study of the complex of 1,8-bis(dimethylamino)naphthalene with 1,2-dichloromaleic acid (a `proton sponge') (Mallinson *et al.*, 1997).

*Example 3.5.3.1. Definition of the local Cartesian axis system for the atoms in a proton sponge complex.*

Data items in this category are as follows:

ATOM_RHO_MULTIPOLE
*The bullet • indicates a category key. The arrow is a reference to a parent data item.*

Data items in this category are also looped. The category key ( ** _atom_rho_multipole_atom_label**) is also a child of

**in the ATOM_SITE category, thereby linking the electron density multipole coefficients to a specific atom site in the list of atom coordinates. This category contains all the parameters of the Hansen & Coppens (1978) functionwhere**

`_atom_site_label`**r**is a position vector with respect to the atomic nucleus, and and are the spherical core and valence electron densities, respectively. They are obtained from atomic orbital analytic wavefunctions such as those tabulated by Clementi & Roetti (1974). They are also the Fourier transforms of the X-ray scattering factors given in

**and**

`_atom_rho_multipole_scat_core`**.**

`_atom_rho_multipole_scat_valence` is the weight applied to the free-atom core (** _atom_rho_multipole_coeff_Pc**), is the weight applied to the free-atom valence shell (

**), is the number of remaining electrons (**

`_atom_rho_multipole_coeff_Pv`**) and (the atomic number) for a neutral atom. is the valence electron expansion factor (**

`_atom_rho_multipole_coeff_P00`**), is the radial function (Slater or equivalent) (**

`_atom_rho_multipole_kappa`**), is the multipole function expansion factor (**

`_atom_rho_multipole_radial_*`**), are the spherical harmonic coefficients (**

`_atom_rho_multipole_kappa_prime[l]`**) and is the spherical harmonic of order at the angle . The summations are performed over the index ranges , , where is the highest order of multipole applied.**

`_atom_rho_multiple_coeff_P[lm]`Example 3.5.3.2 demonstrates how the category is used in the proton sponge complex of Example 3.5.3.1. Only the first atom is shown in the example.

The first implementation of the dictionary appeared in a version of the *XD* program package (Koritsanszky *et al.*, 2003). Thus *XD* can read and write CIFs which include data items in the two categories described above (as well as, of course, items in the core CIF dictionary). It is envisaged that future developments to the rhoCIF dictionary will add features that are relevant to other widely used program packages. In the *XD* implementation, use was made of the *CIFtbx* library of Fortran functions for programming CIF applications (Hall & Bernstein, 1996).

### Acknowledgements

We thank the users of the program package *XD* (Koritsanszky *et al.*, 2003) who provided feedback during the development of the rhoCIF dictionary and Brian McMahon for helpful comments on its presentation.

### References

Bader, R. F. W. (1990).*Atoms in molecules: a quantum theory.*Oxford University Press.

Clementi, E. & Roetti, C. (1974).

*Roothan–Hartree–Fock atomic wavefunctions. Basis functions and their coefficients for ground and certain excited states of neutral and ionized atoms.*

*At. Data Nucl. Data Tables*,

**14**, 177–478.

Coppens, P. (1997).

*X-ray charge densities and chemical bonding.*Oxford University Press.

Hall, S. R. & Bernstein, H. J. (1996).

*CIF applications. V. CIFtbx2: extended tool box for manipulating CIFs.*

*J. Appl. Cryst.*

**29**, 598–603.

Hansen, N. K. & Coppens, P. (1978).

*Testing aspherical atom refinements on small-molecule data sets.*

*Acta Cryst.*A

**34**, 909–921.

Hohenberg, P. & Kohn, W. (1964).

*Inhomogeneous electron gas.*

*Phys. Rev. (Second Series)*,

**136**, B864–B871.

Kohn, W., Becke, A. D. & Parr, R. G. (1996).

*Density functional theory of electronic structure.*

*J. Phys. Chem.*

**100**, 12974–12980.

Koritsanszky, T., Richter, T., Macchi, P., Volkov, A., Gatti, C., Howard, S., Mallinson, P. R., Farrugia, L., Su, Z. & Hansen, N. K. (2003).

*XD*.

*Computer program package for multipole refinement and topological analysis of electron densities from diffraction data.*http://xd.chem.buffalo.edu .

Koritsanszky, T. S. & Coppens, P. (2001).

*Chemical applications of X-ray charge density analysis.*

*Chem. Rev.*

**101**, 1583–1627.

Labanowski, J. K. & Andzelm, J. W. (1991).

*Density-functional methods in chemistry.*New York: Springer-Verlag.

Mallinson, P. R., Wozniak, K., Smith, G. T. & McCormack, K. L. (1997).

*A charge density analysis of cationic and anionic hydrogen bonds in a `proton sponge' complex.*

*J. Am. Chem. Soc.*

**119**, 11502–11509.

Nagy, A. (1998).

*Density functional. Theory and application to atoms and molecules.*

*Phys. Rep.*

**298**, 2–79.

Spackman, M. A. & Brown, A. S. (1994).

*Charge densities from X-ray diffraction data.*

*R. Soc. Chem. Ann. Rep. Sect. C*, pp. 175–212.

Stewart, R. F. (1973).

*Electron population analysis with generalized X-ray scattering factors: higher multipoles.*

*J. Chem. Phys.*

**58**, 1668–1676.

Ziegler, T. (1991).

*Approximate density functional theory as a practical tool in molecular energetics and dynamics.*

*Chem. Rev.*

**91**, 651–667.