IUCr journals news
Fine-tuning the rigid pseudoatom model
A. Volkov & P. Coppens, Acta Cryst. (2001). A57, 395–405
The recent adoption of high-speed data-collection methods, especially CCD detectors, and the wider availability of user-friendly software such as XD (Koritsanszky, Howard, Mallinson, Su, Richter & Hansen, 1995) have attracted new workers to charge-density studies from X-ray diffraction data, extending the variety of the systems studied, as well as leading to greater penetration of the results in mainstream chemistry and physics journals. However, it is increasingly clear that limitations of the rigid pseudoatom model (multipole refinement) upon which most studies are based may be compromising the more ambitious studies, especially topological analyses, where the objective is to determine the fine details of the curvature of the electron density.
There seem to be two ways forward: (i) fine-tune the rigid pseudoatom model, especially the radial functions, or (ii) abandon the model completely. Examples of the latter approach include maximum entropy methods and the derivation of 'experimental wavefunctions' from X-ray diffraction data. However, the vast majority of charge-density determinations still use the multipole refinement procedure, and numerous model studies have appeared recently that explore its strengths and limitations. Especially active in recent years has been the Coppens group in Buffalo, and the latest publication by Volkov & Coppens [Acta Cryst. (2001). A57, 395–405] explores the use of additional radial functions on non-H atoms.
A more flexible pseudoatom may be obtained by adding a second radial function per multipole, with a different set of n and α. But what are the optimum values to choose for each atom? To answer this question, Volkov and Coppens performed several thousand model refinements, optimizing α for choices of n between 4 and 22 for a second set of radial functions on first-row atoms. Inclusion of additional radial functions led to significant improvements in the agreement with the theoretical electron density and topological properties along the bonds, but discrepancies still remained, especially for polar bonds such as N–O and C–O. The figure illustrates the improvement in the electron density for an O–C bond in methyl carbamate. Surprisingly, optimum values of n for the second deformation set were found to be as high as 20 and 22 in several instances, corresponding to quite narrow radial functions. The authors did not explore why this might be the case for some bond types, and do not provide a standard recipe for incorporating additional radial functions in multipole refinements. They conclude that application of this more flexible model will require experimental structure factors of exceptional accuracy, which is undoubtedly true if one aims to refine n and α simultaneously. However, it is worth remembering that Hirshfeld’s multipole model (Hirshfeld, 1977), seldom used today, has always incorporated two-term radial functions for some higher multipoles, albeit in a constrained fashion, and its demands on the X-ray data are much more modest. It would seem desirable to compare its performance with the one- and two-term Hansen–Coppens models. This paper provides some valuable guidance and much-needed impetus for future studies aimed at fine-tuning the multipole model.
Hirshfeld, F. L. (1977). Isr. J. Chem. 16, 226–229.
Koritsanszky, T., Howard, S. T., Mallinson, P. R., Su, Z., Richter, T. & Hansen, N. K. (1995). XD. A Computer Program Package for Multipole Refinement and Analysis of Charge Densities from X-ray Diffraction Data. Free University of Berlin, Germany.
School of Physical Sciences & Engineering,
U. of New England, Armidale, NSW, Australia