Charge density analysis for crystal engineering
Chemistry Central Journalvolume 8, Article number: 68 (2014)
This review reports on the application of charge density analysis in the field of crystal engineering, which is one of the most growing and productive areas of the entire field of crystallography.
While methods to calculate or measure electron density are not discussed in detail, the derived quantities and tools, useful for crystal engineering analyses, are presented and their applications in the recent literature are illustrated. Potential developments and future perspectives are also highlighted and critically discussed.
In modern crystallography, a crucial issue is the understanding of interactions that enable the assembly of molecules and the fabrication of flexible or rigid organic or metal-organic polymers.
The supra-molecular paradigm is often associated with crystal engineering. This name was originally introduced by Pepinsky , but later used by Schmidt  with a different meaning namely the usage of crystals for controlled stereospecific chemical reactions. In Schmidt’s view, the crystal is a matrix in which the reaction occurs and, at the same time, a precursor of the desired material. On the other hand, crystal engineering has later evolved towards the rationalization of binding motifs and their usage to create crystalline materials with specific structural or functional features : the crystal and its structure have become the subject themselves of the speculation and the target of the research. Crystal engineering is the initial and fundamental step leading to the fabrication of a material and it implies the design, the preparation and the characterization of crystalline species.
In this context, the accurate analysis of those linkages that build up the desired structural motifs, are extremely important. Most of these bonds are, however, more elusive than typical chemical bonds of organic molecules, whose nature is known and well rationalized since decades.
Coordinative bonds in metal organic frameworks are most of the time well known because identical to those typical of simple complexes and often understood within the ligand field theory . On the other hand, it is the regio-selectivity in multi-dentate organic linkers to be more intriguing and sometimes difficult to predict.
Even more complicated is understanding the nature and the role of various intermolecular non-covalent interactions in crystals based only on weaker forces, see Table 1 for a summary. This field has attracted enormous attention, starting from the most well-known of these interactions, namely the hydrogen bond  (HB).
Recognition and classification of intermolecular bonding features is important not only to understand the key factors that promote aggregation, but also to enable the classification of solids through topological analysis , which is a method to rationalize both the structural motifs and, at least in principle, the resulting material properties, thus the fundamental steps of a proper material design.
Since the early days of X-ray diffraction, it became clear that it was in principle possible not only to ascertain the positions of atoms in crystals, but also to observe the distribution of electrons  and therefore to “visualize” the chemical bonding. This became really feasible much later  and it is nowadays quite common to analyze molecules and crystals in terms of electron density partitioning . Among the most relevant achievements, important is the analysis of chemical bonding, through the quantum theory of atoms in molecules (QTAIM) , which has been successfully applied to coordinative bonds , as well as to most of the known intermolecular interactions –. Moreover, electron density partitioning enables the evaluation of electrostatic interactions between molecules, therefore provides quantitative measures of involved energies.
Methods to obtain the electron density experimentally or to calculate it by first principles are well known and explained in textbooks , and review articles , and we will not focus on that in this paper.
Here it is important to recall the following concepts and notions:
The electron density (ED, ρ(r)) is a quantum mechanical observable, that can be measured, for example, through scattering experiments, in particular X-ray diffraction from crystals. Although it can be directly calculated by Fourier summations of structure factors, the electron density is better obtained as a three-dimensional function fitted against the measured structure factors, which enables a deconvolution of the atomic thermal motion from the (static) electron density distribution.
The most adopted method to reconstruct the electron density is the multipolar model ,, where ρ(r) is expanded into atomic - or better pseudo-atomic - multipolar functions, based on a radial function centered at the nuclear site and an angular function (spherical harmonics, usually truncated at hexadecapolar level).
While the multipolar electron density is not a true quantum mechanical function, it can be compared to those computed ab initio by quantum chemical methods that use various degrees of approximation to solve the Schrödinger equation.
From the electron density some important properties are straightforwardly calculated, like the electrostatic potential, field and field gradients or the electrostatic moments of an atom, a functional group, a molecule or a monomeric unit of a polymer. These partial quantities require that an assumption is made on how to recognize an atom in a crystal (and therefore a functional group or a molecule). The most adopted scheme is offered by QTAIM, but other recent applications make use of Hirshfeld “stockholder” partitioning, as for example the Hirshfeld atom  or the Hirshfeld molecular surfaces –.
Other important properties cannot be obtained from the electron density, because they would require not only the trace of the first order density matrix (from which the Bragg scattering of X-rays depends), but also the out of diagonal component or the second order density matrix. These quantities, albeit connected to observables and experimentally available quantities, are very difficult to measure and more often they are obtained only via theoretical calculations.
In the past two decades, methods have been proposed that directly refine elements of density matrices or coefficients of a quantum mechanical wavefunction, including information from scattering experiments (X-ray diffraction, Compton scattering, polarized neutron scattering), see  for a comprehensive review on the subject. These approaches, albeit less straightforward than the traditional multipolar expansion, are extremely appealing because they combine theory and experiment and offer a wider spectrum of properties, because the full density matrix becomes available.
In this paper, we present some of the many tools offered by electron density analysis for crystal engineering studies and we will show some applications reported so far in the literature, giving some perspectives for future developments in this field.
Electron densities of studies of organic crystals
Characterization of Intra- and Inter-molecular interactions
As introduced above, at the basis of crystal engineering is the understanding how molecules interact with each other to form a three-dimensional structure in the solid state. The more insight we get into the nature of these weaker, intermolecular bonding, the more effective materials we can obtain.
With no doubts, QTAIM is one of the most powerful tool for evaluating the interactions within crystal structures, because it analyzes the gradient field of the electron density, hence it enables visualizing its concentrations and depletions, knowing that electrons are in fact the “glue” that stick atoms and then molecules together. QTAIM is grounded on the idea that atoms can still be identified in molecules and provides the quantum mechanical bases for that ,. This justifies the hard space partitioning of ρ(r) into atomic basins (Ω) used to quantify atomic volumes and electron populations. The inter-atomic surface (IAS) shared by two bonded atoms enables to evaluate the nature of the bonding between them, especially analyzing the electron density properties at the bond critical point (BCP), a point on IAS where the gradient of ED is equal to zero (∇ρ(r ) = 0). In order to extract chemical information on the bond, such as its strength, order, polarity etc., properties evaluated at BCPs are crucial. One of the most important electronic property at BCP is the Laplacian of electron density, ∇2 ρ(r ). Bader et al.  noted that covalent bonds are typically associated with the approach of the valence shell charge concentrations of the bonded atoms, producing a local accumulation of charge at the BCP, thus characterized by a negative ∇2 ρ(r ). On the contrary, a positive Laplacian indicates the local depletion of electron density, typical of closed-shell interactions, i.e. interaction between two electronic systems with the outermost electronic shell filled, as it occurs in ionic bonds, or any other interaction between molecules (van der Waals, medium-weak hydrogen bonding etc.).
This paradigm works well for most organic compounds but it fails when heavier atoms (e.g. transition metals) are concerned . In fact, the rather elusive outermost shell of these elements, makes the sign of ∇2 ρ(r ) no longer discriminating. For almost all bonds to a transition metal, the corresponding BCPs are found in regions of charge depletion ,, thus producing a kind of “Hegelian night”. For this reason, other indicators were found to be more useful, for example the energy densities and the electron delocalization indices that however require the entire first order density matrix to be calculated, therefore they cannot be retrieved just from the electron density (trace of the first order density matrix), which is a more straightforward observable. A local kinetic G(r ) and potential V(r ) energy density functions can be defined from the first order density matrix. Cremer and Kraka  were the first to introduce the idea that the total energy density H(r ) (=G(r ) + V(r)) reflects a dominant covalence when H(r ) <0 (i.e. when the potential energy density is in excess). As a matter of fact, the total energy density, better than the Laplacian, defines sensible boarders of a molecule, see for example Figure 1, in which the H(r ) distribution of two approaching glycine molecules is drawn. When a strong hydrogen bond is eventually formed, the valence regions of the two molecules belong to the same synaptic domain of negative energy density, thus H(r ) at the H…O BCP is negative, in keeping with most of the current consensus that for such short distances the bond must contain significant amount of covalence. Electron delocalization indices, instead, measure the number of electron pairs shared by two atomic basins  and can also be used to reveal the degree of covalence in intermolecular bonding. In addition, they may anticipate exchange paths in magnetic frameworks, especially in metal based materials. The nature of metal-ligand bonds is extremely important in one sector of crystal engineering, namely that of coordination polymers, and will be discussed in section 3. Other interactions between molecules are more genuinely classified into the closed-shell class, although some amount of covalence might be present, sometime playing a fundamental role.
In order to characterize genuine hydrogen bonds and to differentiate them from pure van der Waals interactions, Koch and Popelier , proposed eight conditions to be fulfilled. The first four criteria can be easily checked also from experimental data, because relying only on the electron density function and its derivatives in the crystal, whereas others comparisons of quantum mechanical calculations for the HB aggregate and the isolated molecules. The presence of a BCP between donor and acceptor atoms linked through a bond path and the presence of charge density at the BCP is the basis of first two criteria. Positive value of the Laplacian at BCP and its correlation with interaction energies constitutes the third condition. Noteworthy, this condition can be controversial because very strong, symmetric HB’s are associated with negative Laplacian, indicating even large stabilization energy. The fourth criterion, considered as “necessary and sufficient”, concerns the mutual penetration of the hydrogen (H) and the acceptor (A) atomic basins. The following relation must be fulfilled:
where are non-bonded radii of hydrogen and the acceptor atom taken as the gas phase van der Waals radii and r H , r A are corresponding bonding radii taken as distances from BCP to the nuclei. Any violation of the above condition indicates van der Waals nature of the considered contact. Other criteria express a loss of electrons and energetic destabilization of H-atom resulting from increased net charge of the atom, as well as a decrease of dipolar polarization and volume depletion of H-atom.
Another convenient classification of weak electrostatic interactions is based on the electronic energy densities, introduced by Espinosa et al. , who extended the idea by Cremer and Kraka. In fact, weak electrostatic interactions, can be classified in terms of kinetic energy density G(r BCP) and potential energy density V(r BCP) at the bond critical point. The relationship between those two functions reflects how electrons around BCPs are affected by the formation of a hydrogen bond (HB). As mentioned above, energy densities in principle require the full density matrix to be computed, however Abramov  proposed a functional to estimate the kinetic energy density, based only on the electron density and its derivatives, therefore making it available to experimental determinations as well. In particular, at points where ∇ρ(r ) vanishes (like all the critical points of ρ(r ), CP’s):
In turn, the potential energy density V(r CP) is then obtained applying the local virial theorem:
Following Cremer and Kraka, in closed-shell interactions, the local kinetic energy density G(r CP) (everywhere positive) is in excess of local potential energy density V(r CP) (everywhere negative), thus H(r CP) >0. Furthermore, the larger is |V(r CP)∣, the larger is the shared character of the interaction and the electronic stabilization of the structure. It is also observed that in closed-shell interactions the amount of kinetic energy per electron is large, typically G(r CP)/ρ(r CP) >1 (in atomic units). Because at BCP the kinetic and potential energy densities depend exponentially on the distance between hydrogen atom and the HB acceptor, a correlation was found between the energy of the hydrogen bond and potential energy density :
which can be interpreted as the energetic response of the hydrogen bond to the force exerted on the electrons around BCPs. Note that the ½ coefficient in equation (4) is not dimensionless. Spackman  has shown that this correlation can be actually predicted even from the pro-crystal ED, i.e. from the simple summation of atomic spherical electron density functions that are easily calculated once the crystal structure is known.
The above mentioned quantitative indicators can be used to analyze any weak interactions and were extensively applied by Munshi and Guru Row . They first reported on a comparison between experimental and theoretical electron density studies for three bioactive molecules: 2-thiouracil, cytosine monohydrate and salicylic acid. They gave a quantitative description of all the interactions and could clearly differentiate strong and weak contacts. Moreover, they showed that the nature of weak interactions is not lost in the presence of strong hydrogen bonds. Those studies contributed to evaluation of preferred orientations at the protein binding sites. The same group studied for the first time the differences of the electron density in polymorphs of 3-acetylcoumarin . This research clearly indicated that for the purpose of “quantitative crystal engineering” conventional crystal structure analysis based only on geometrical features is insufficient and inadequate. Only detailed ED analysis can justify the occurrence of any interaction in the crystalline state and therefore provide a useful input to design new materials. Extensive studies on aliphatic dicarboxylic acids  revealed interesting systematics in the topology of ED. The electron density associated with the side-chain interactions, as a fraction of the total intermolecular density, plotted against the number of methylene groups revealed an alternating behavior. Acids with even numbers of carbons exhibit higher ρ(r) values at bond critical points, compared to their odd neighbors. This explains the relatively higher melting points in the even-member acids since side-chain interactions play a major role for the cohesion of acid molecules in the solid state. Howard and co-workers , based on experimental ED studies of trans-cinnamic acid and coumarin-3-carboxylic acid, postulated that the presence of strong interactions not maintained in the irradation products “may influence the ability of a compound to undergo a solid-state [2,+2] cycloaddition reaction”.
A general hypothesis concerning azide building blocks was proposed by Bushmarinov et al. . Based on the QTAIM and the electron localization function (ELF ,), geometrical preferences in favor of hydrogen bond formation were explained. They proved that the number of interactions to the terminal nitrogen atoms of the azide only depends on steric effects, thus supramolecular systems based on hydrogen bonds to an azide will be independent from the torsions involving terminal atoms of the azide.
Important insight into the understanding the crystal structure –property correlation was given by Gopalan et al. . They studied organic crystals exhibiting NLO properties and showed how the non-centric nature of the crystal field affects molecular dipole moment and therefore optical properties of the solid.
Increasing attention is attracted by halogen bonded – crystals. In terms of charge density analysis, Bianchi et al.  reported on the investigation of the co-crystal of 1,2-bis(4-pyridil)ethylene with 1,4-diiodotetrafluorobenzene. Based on QTAIM topology, they classified the interaction between the pyridyl donor N and the di-iodobenzene I acceptor as a closed-shell interaction. Indeed this is characterized by a positive Laplacian at the intermolecular BCP, although accompanied by a negative energy density as the authors also pointed out. A clear manifestation of the “Hegelian night” is that C-I interactions would appear, at first sight, “similar to those of metal-metal and metal-ligand bonds in organometallic compounds” . Interestingly, the authors proved that equation (4) remains substantially valid for this interaction, by comparing the ab initio interaction energies with the empirical derivation from the kinetic energy density.
Bui et al.  developed a model to rationalize halogen bonding based on accurate studies on hexa-halobenzene molecular crystals. The deformation density, i.e. the difference between the total electron density and the superposition of spherical atomic densities (hereinafter called the promolecule), enabled the visualization of the so-called σ-hole, first anticipated by Politzer et al. , using quantum chemical calculations. Other studies have followed –, again stressing on the visualization of the σ-hole by means of the electrostatic potential or the Laplacian distribution, and therefore addressing the overwhelming contribution of the electrostatic term. However, recent work by Stone  has demonstrated that some stereochemical features of the halogen-bonded packing originate from the necessity to minimize the inter-atomic repulsion term, rather than from a stabilizing, though weak, electrostatic interaction. Accordingly, Spackman has very recently shown that in many cases the interaction between halogen bonded molecules is associated only with a small or even negligible stabilization . Therefore, further investigations are expected in the next future on this topic.
To facilitate the discussion of all intermolecular contacts in molecular crystals it is very useful to introduce Hirshfeld surface (HS) analysis –. The Hirshfeld molecule is an extension of the concept of Hirshfeld atom , which is not based on a quantum mechanical definition, as QTAIM, but on a rather simplified interpretation of a multivariate function, like the electron density when a breakdown into atomic terms is adopted. Hirshfeld defined the atom as a “stockholder”, who receives from the “asset” an “equity” proportional to the “investment”. In this naive example, the asset is the molecular electron density (computed or measured), whereas the investment is the electron density of the isolated atom, calculated in its ground state and spherically averaged. The equity, evaluated at each point and integrated over whole space, can be positive or negative, leading to a negatively or positively charged atom respectively a.
In case of a molecule in a crystal, mutatis mutandis, the same concept can be applied. However, Spackman realized that the fuzzy partitioning of the Hirshfeld approach (each point in space belongs to many atoms, with its own share) was not very useful for crystal engineering. A hard space definition of the building blocks is much preferable. Therefore, he defined a molecule in a crystal as the, unique, region of space whose procrystal density has at least 50% share from the given pro-molecule. Noteworthy, a tessellation of space is not complete with this partitioning, because regions without a dominant pro-molecule are in principle possible, albeit in general extremely small. The Hirshfeld surface gives a unique signature of a molecule in a crystal, because it strongly depends on the surrounding, so the same molecule in different crystal packing looks different. On the Hirshfeld surface, some functions can be mapped, as for example d norm , which combines the internal di and external de distances from the surface to the nearest nucleus. On Figure 2a, the Hirshfeld surface of L-aspartic acid (L-Asp) is shown: contact zones shorter than van der Waals radii are marked as red areas and highlight hydrogen bond sites of the molecule. Hirshfeld surfaces are very often accompanied by 2D fingerprints ,, Figure 2b, scatter-plots of de and di that uniquely identify each type of interaction in the crystal structure. In case of L-Asp, the strongest interactions are those of O…H type constituting the highest fraction of 72.7%. Other close contacts are also present, including very weak C…H interactions (2.8%) and non-directional H…H contacts contributing in 18.9%.
Beside the numerous applications of this methodology and its growing appeal in crystal engineering, it should be stressed that this analysis does not rely on quantum mechanics and therefore its predictive power is based only on empirical evidences.
A recent and alternative way of quantifying non-covalent interactions (NCIs) between molecules was introduced by Johnson et al.  and Contreras-García et al. . The NCI descriptor enables visualizing regions of space involved in either attractive or repulsive interactions. The NCI index depends on the reduced electron density gradient (RDG):
Scatterplots of s(r) against ρ(r) address non-covalent interactions. In fact, in the low-gradient and low-density regions characteristic spikes occur which are not observed for covalent bonds. If we only consider ED/RDG regions, the information about the nature of the interaction would be lost, since different types of interactions appear in the same very narrow range. However, the sign of the second eigenvalue λ 2 of the Hessian matrix of ρ(r) (with ∇2 ρ(r ) = λ 1 + λ 2 + λ 3); λ 1 ≤ λ 2 ≤ λ 3) indicates whether the interaction is stabilizing (λ 2 < 0) or destabilizing (λ 2 > 0). Therefore, diagrams of s(ρ(r) ⋅ sign(λ 2)) allow recognizing the type of NCI, whereas the amount of density itself issues the strength of that interaction. A spike in the low-gradient, low-density area at negative λ 2 indicates stabilizing interactions like hydrogen bonds, a smaller spike and slightly negative λ 2 is the fingerprint of a weaker stabilizing interaction, and a spike associated with positive λ 2 indicates that no NCI is active. The shape of RDG surfaces also allows for qualitative description of interactions strength. Small disc-shaped RDG domains denote stronger interactions whereas broad multiform domains refer to much weaker interactions. NCI approach can be applied to experimental or theoretical electron density distributions, as for example shown by Saleh et al.  or by Hey et al. .
The importance of intermolecular interactions can be evaluated also through the analysis of atomic polarizabilities, in particular their deformation with respect to non-interacting molecules. Recently, we have developed a program, PolaBer , which enables to calculate distributed atomic polarizabilities based on a partitioning of the electron density. The advantage of this approach is the definition of atomic contribution to a molecular property (the molecular polarizability) or a crystalline property (the linear susceptibility), which enables to identify the key-features for large polarizabilities. Therefore, this approach might be useful for crystal engineering purposes.
The ED partitioning follows QTAIM, although other schemes could be adopted. The main advantage of QTAIM is that it is based on quantum mechanical ground, therefore together with atoms in molecules one consistently define bonds as well. Moreover, it ensures a maximal transferability between different systems as already demonstrated by Matta and Bader . Within this approach atomic properties such as charges Q(Ω), energies E(Ω) and, in particular, dipole moments μ(Ω) can be calculated by integrating their corresponding operators over the volume of atomic basin Ω. Atomic polarizability tensors are obtained from numerical derivatives of atomic dipole moments with respect to external electric field:
where is the atomic dipolar component along the i direction computed with a given electric field (0 or ε) in j direction. Full description of the procedure is given in Krawczuk et al. , based on the theory developed by Keith . For crystal engineering purposes, it is essential to differentiate weak non-covalent intra-molecular or intermolecular interactions from covalent bonds inside the molecule. In fact, the partitioning scheme distinguishes two contributions to the atomic dipole: one is due to the polarization inside the atomic basin, the other originates from distributing the atomic charge over all the bonds to the atom creating a bond dipole. These quantities are easily computed from a system of equations involving all bonds and all atomic charges, however an ambiguity occurs when a ring is present. Keith  suggested including an additional condition to enable solution of the system of equations: the sum of ring bond charges should be zero. However, if all bonds are taken as equivalent in the ring, an anomalous importance is attributed to weaker interactions, producing mathematically correct but physically unrealistic atomic polarizabilities. Therefore, a weighting scheme is applied in PolaBer: in the ring conditions, bond dipoles are multiplied to the inverse of their strength, measured by the electron density at the BCP. This avoid drastic changes of the atomic polarizabilities, if a weak BCP generates a ring in the molecular graph. On the other hand, rings made of strong covalent bonds (like those of aryls) truly affect the atomic polarizabilities; accordingly all bonds have similar or even identical weight if symmetric. On Figure 3 atomic polarizabilities in L-valine are visualized. In the zwitterionic form, an intramolecular weak hydrogen bond of C − H…O type is present. If no weighting scheme is applied, the polarizabilities of oxygen and hydrogen atoms are substantially different than those of the same molecule in a conformation where no intramolecular HB occurs.
Since atomic polarizabilities are second order positive tensors, they are easily visualized as ellipsoids with main axes having dimensions of volumes. The visualization is done in the same real space as the molecule assuming that 1Å3 (unit of polarizability tensor) =1Å (unit of atomic coordinates), though for visualization purposes a scaling factor is necessary to reduce the size of ellipsoids (typically 1Å3 = 0.4Å for atoms of the second period). The size of the ellipsoid is proportional to the total atomic polarizability, whereas the ellipsoid axes indicate the anisotropy of the polarizability, thus the directions along which the atomic electron density is more or less polarizable. Although weaker than covalent bonds, hydrogen bonds may affect the polarizabilities. The perturbation is mainly due to the electrostatic interaction occurring between the donor and the acceptor atoms in the hydrogen bond system. On Figure 4, a comparison between isolated molecule of oxalic acid and a dimer is shown. In general, polarizabilities are larger along covalent bond directions, especially towards atoms with high polarizabilities (see carboxylic groups). When a hydrogen bond is formed, the oxygen atoms are slightly modified in orientation and are stretched along HB direction (compare O1 and O2 ellipsoids in both pictures), due to the perturbation produced by the incoming donor atom. The increased polarizability along the direction of the HB can be measured by the bond polarizability defined as the projection of the atomic ellipsoids on the bond vector:
where r ΩΩ' is a unit vector in the direction of Ω-Ω’ bond. α Ω − Ω Ω' is a scalar which reflects how feasible is the polarization of the electron density along the bond, upon application of an electric field in the same direction. Values of bond polarizabilities for carboxylic groups of oxalic acid are also given on Figure 4. Larger values of bond polarizabilities of O-H bond in dimer confirm the elongation of hydrogen polarizability along the donor-acceptor path.
The design of co-crystals for multifunctional materials has brought lots of attention in last few years, especially in the field of pharmaceutical compounds – where at least one of the components is an active pharmaceutical ingredient (API). So obtained co-crystals gain new chemical and physical properties (i.e. solubility, density, hygroscopic abilities, melting point etc.), usually drastically different from individual components and ideally tunable in order to obtain the desired functionality. The crucial point in crystal engineering of drugs is to understand and evaluate potential intermolecular interactions that a given molecule may exhibit and rationalize the consequences for the supramolecular architecture.
One of the first charge density analysis on API co-crystals was presented by Hathwar and co-workers . The main goal of the study was to quantitatively describe differences between a co-crystal of nicotinamide (API component) with salicilic acid and the salt formed by nicotinamide and oxalic acid. The region of main interest was the proton transfer path to the nitrogen atom on the pyridine ring of nicotinamide. Topological analysis revealed bonding features associated with N…H − O and N − H…O hydrogen bonds for the salt and the co-crystal, respectively. A similar picture was obtained from the electrostatic potential maps where the electropositive region on oxygen atom of salicilic acid indicated close-shelled interaction whereas electronegative region of oxygen atom on oxalic acid suggested covalent bond with H atom. All above observations confirmed earlier hypothesis that a formation of a co-crystal is strongly dependent on pKa of the individual constituents. Those studies offered very convenient way of verifying the continuum from co-crystal to salt by assessing interaction energies in terms of charge transfer character at the critical point.
Hathwar et al.  proposed also a library of transferable multipolar parameters for structural fragments representing supramolecular synthons. Based on the high resolution X-ray diffraction datasets, the library would provide the criteria to design and fabricate new synthons and therefore mimic the 3D formation based on a given hydrogen bond system. Since it was already proven that multipolar parameters are transferable for molecules or molecular fragments – authors wanted to test if this is also true for supramolecular synthons. The transferability was tested on methoxy-benzoic acid, acetanilide, 4-(acetylamino)benzoic acid, 4-methylbenzoic acid, and 4-methylacetanilide. Electron density features derived with the supramolecular synthon based fragments approach (SBFA), were compared to experimentally obtained values and showed a very good agreement, except for some discrepancies in monopole parameters. The SBFA can be successfully applied for essential topological features of ED for intra- and intermolecular interactions (synthons) in molecular crystals, especially when no good quality crystals can be obtained and therefore no high-resolution data can be gathered. SBFA model was also applied for synthons with C − H…O, C − H…F and F…F interactions  proving that this model can also be applied for weak interactions.
Charge density studies provide valuable information on subtle features in case of polymorphism in co-crystals. In our best knowledge, so far only couple of papers reported on charge density studies for polymorphs of co-crystals. Gryl et al.  confirmed earlier hypothesis  that the polymorphic forms of barbituric acid and urea originate from the existence of resonance structures of the barbituric acid molecule. Both, experimental and theoretical charge density studies indicated characteristic features of two, among six, possible mesomeric forms, see Scheme one in . It was possible to recognize electron density displacement in barbituric acid molecule towards those two resonant forms, which influence the type of hydrogen bonds formed in each polymorphic form and therefore results in different packing topology.
Schmidtmann et al.  studied short, strong hydrogen bonds (SSHB) of O − H…N type in two polymorphic forms of isonicotinamide-oxalic acid crystallizing in C 2/c (I) and P‾1 (II) space groups. It was the first case where topological analysis of ED confirmed formation of rather unusual centered heteronuclear intermolecular SSHB of O…H…N type between oxalic acid and isonicotinamide in polymorph II. The presence of such an interaction, where carboxylic H atom is equally shared between O and N atom, raises the question whether this compound should be considered as a co-crystal or a salt, however this is beyond the scope of this paper.
Dubey et al.  applied supramolecular synthon based fragments approach  to study polymorphism of orcinol:4,4′-bipyridine co-crystals and showed the transferability of multipole parameters of O − H…N synthon in the polymorphic forms of the studied compound.
Although a rather small number of studies have reported on co-crystals from charge density point of view, all above examples clearly indicate the importance of the topology of electron density towards understanding the correlation between crystal structure and the physical property of the materials.
Optical properties from electron density studies
One of the most challenging part of crystal engineering is to design an efficient material which will exhibit desired physical properties. To accomplish this goal, a deep understanding of the connection between molecular structure and the solid-state property is needed. In this chapter we want to focus on optical properties and attempt to estimate those properties from electron density studies, both experimental and theoretical.
The optical response of the material is determined by its electric susceptibilities, that depends on the dipole (hyper)polarizabilities of individual atoms and therefore molecules. The electric (hyper)polarizabilities of a molecule have long been used to predict and understand their chemical reactivity, intermolecular interactions and physical properties. For example, the first order polarizabilities enable to calculate refractive indices, applying the Clausius-Mossotti equation , or, taking into account long range interactions, the anisotropic Lorentz field factor approach .
There were several attempts to obtain molecular (hyper)polarizabilities from electric moments using experimental charge density studies of materials with potential non-linear optical properties and therefore to predict the size of the NLO effect. Within first reports –, the estimation of NLO properties was done from a Robinson model  allowing to connect molecular polarizabilities to multipolar moments of the electron density distribution. First results were very promising, especially in case of molecular polarizability, but failed when estimating the hyperpolarizabilities. Explanation for this was given later by Whitten et al.  who stated that the electron density obtained from regular multipole model does not carry sufficient information to determine reliably molecular hyperpolarizability and therefore to estimate correctly the non-linear optical effect. Instead, the X-ray constrained wavefunction (XCW, ) approach, implemented in the software package TONTO , provides a pseudo-quantum mechanical wavefunction constrained to reproduce the experimentally observed structure factors. The wavefunction calculations minimize the sum of Hartee-Fock energy and the χ2 function which defines the precision level of the experimental structure factor against the calculated ones:
where N r and N p are numbers of reflections and parameters, respectively, F and F * are calculated and experimental structure factors and σ 2 is the uncertainty of each structure factor. This leads to an “experimental wavefunction” from which physical properties can be easily calculated, thanks to the fact that a pseudo density matrix can be defined. In a recent paper , the XCW model was successfully applied to calculate molecular dipole moments, polarizabilities and refractive indices of small organic compounds (which find application in laser dyes) using coupled perturbed Hartree-Fock (CPHF) approach, therefore calculating the field induced perturbation to the molecule embedded in the crystal. A new XCW approach was recently proposed by Genoni , who implemented the method for extremely localized molecular orbital wave functions. An alternative way to estimate the crystal properties is however that of calculating with high accuracy the molecular quantities, then applying corrections for the perturbation of the crystal packing. For example, quite effective is the distributed atomic polarizability approach to estimate linear optical properties implemented in PolaBer . The atomic polarizabilities can be used to calculate the electric susceptibility, through the anisotropic Lorentz approximation. The advantage of using atomic polarizabilities rather than the whole molecular ones, is that we can extract separate information about the role of each functional group in the molecule, which is very important to design new molecules. One nice example is the calculations of refractive indices recently reported for the L-histidinium hydrogen oxalate crystal structure . Obtained values of refractive indices were comparable with the ones obtained from couple-perturbed Kohn-Sham theory, although slightly different from the experimental ones. In this kind of comparison, one should consider that calculations are generally carried out at zero frequency instead of finite one, and therefore refractive indices are underestimated. Nonetheless, these results are promising and could open a new field in applications of electron density partitioning for material properties.
Metal organic materials
As discussed in the previous sections, the electron density analysis offers many tools to investigate materials, in particular the stereo-electronic features that enable understanding the robustness of a given type of aggregation or the breakdown of a crystal property in terms of atoms, functional groups or molecular building blocks.
For metal organic materials, and in particular porous metal organic frameworks (MOFs) , quite useful is the possibility to “observe” interactions occurring in cavities, channels or layers, where guest molecules or counter ions can be hosted and could diffuse, for example during ion exchange processes.
For example, Hirshfeld surfaces have been adopted not only to define molecules in crystals, but also as a qualitative tool to investigate mutual relationships between building blocks of materials  and to find possible exchange channels for ions . In fact, HS could be used, on one hand, to visualize the complementarity between functional groups in building blocks and therefore to visually address potentially robust synthons. This relatively simple approach has been enthusiastically received in the crystal engineering community, so that HS plots usually accompany many papers in this field, although, as already discussed, the quantum mechanical information therein is sometime overestimated. On the other hand, the procrystal electron densities enable the visualization of sites available for guest molecules and therefore potentially usable channels in porous frameworks. An available site is expected when the procrystal electron density is below a given threshold. Albeit heuristic, this concept implies that a guest molecule can be hosted in a framework only in regions where the short range repulsion associated with the Pauli exclusion-principle (Table 1) is small enough. In fact, it is demonstrated that this repulsion is proportional to the overlapping density , therefore a region of small electron density of the framework should be more accessible. Moreover, as discussed in section 2, the electron density is proportional to the amount of kinetic energy density, which would produce destabilization, therefore the criterion is actually grounded on energy considerations.
This qualitative picture calls for more accurate evaluation of the stabilization or destabilization produced when two molecules interact. Since the 1980s’, Spackman – has proposed simple models to evaluate the electrostatic energy of two sets of multipolar distributions, based on the classical point-multipole approach by Buckingham , though including corrections for the diffuse nature of the electron density distributions that could also penetrate one into the other. Moreover, he proposed a set of atom based parameters to estimate the repulsion as well as the dispersion. In the classical McWeeny approach , the electrostatic energy is simply the zero order energy of the Coulombic interaction between the two electron density distributions, see also Table 1. At this point, it should be reminded that the actually observed electron density in crystals implicitly includes the polarization induced by the electric field experienced by the molecule in the crystal, therefore the electrostatic energy computed from experimentally derived multipolar parameters includes the first order perturbed polarization energy . The dispersion energy, instead, comes from the mutual perturbation of excited configurations and it is often regarded as the electrostatic interaction occurring between instantaneous multipoles of the two interacting systems.
Repulsion and dispersion are normally estimated by atom parameterized functionals, that depends on the nth of the distance between atoms (6th and 12th power for dispersion and repulsion, respectively, in classical Lennard-Jones potential ) or otherwise have exponential decay.
After the pioneering work by Spackman, other approaches have been proposed, which are especially trying to overcome the inaccuracy of a point-multipole approximation. For example, Gavezzotti  has proposed the well-known PIXEL approach, where the electrostatic interaction is evaluated in terms of Coulomb interactions between the charge of the interacting molecules at each point (pixel) in space. This is nothing else than a numerical solution to the Coulomb integral of two charge distributions, that Volkov and Coppens ,, instead solved with a more refined quadrature, obtaining a considerable speed up in the calculation. The more interesting feature of Gavezzotti’s approach was the estimation of polarization energy, genuinely deconvoluted from zero order electrostatic interaction if starting from molecular electron densities evaluated in the gas phase, as well as the estimation of repulsion and dispersion. Repulsion follows the idea of Wheatly , through the evaluating of the overlap integral between two electron density distributions, whereas the functionals for polarization and dispersion  make use of the concept of polarizability functions . Thus, polarization is the energy obtained by summing at each point the interaction between the electric field of one molecule and the local field induced dipole of the other which depends on the molecular (or atomic) polarizability. Dispersion instead is the interaction between both polarizabilities, following the original proposal by London .
The use of the PIXEL semi-classical sums, although very easy to parallelize and therefore suited for modern GPU’s clusters, is not so convenient for larger systems, where faster algorithms like Volkov’s ,, are more efficient. Moreover, some approximation could be adopted to simplify the calculation: in fact it is not necessary to use the exact electron density functions for molecules or atoms which are quite distant from each other. In this case, the point multipole approximation  is sufficiently accurate and much faster to evaluate. For even longer range interaction, the molecular moments, instead of the atomic ones, could be used.
Despite the theoretical background is in fact available, interaction energies are rarely used to investigate MOFs or in general of metal organic materials. Adams and Rao  proposed an energetic criterion based on Morse-type functionals and the bond-valence mismatch, which seems to be promising and quite simple to use, but limited to inorganic systems. Instead, an interaction energy mapping was proposed by Chimpri and Macchi , that while preserving the same simple visualization of threshold density approaches, it contains much more quantitative information because the energy is in fact evaluated in more accurate way using Volkov’s method and therefore it is the total interaction energy to be computed, not only the repulsion. Moreover, the adopted threshold is more physically grounded, because an available s ite is simply defined as the one for which the interaction energy between the guest molecule and the framework is stabilizing. It is very important to stress that this energy is calculated from the electron density of the host and the guest, therefore at variance from the criterion of minimal electron density of the host, the interaction energy map is able to screen the capability of the framework to host a particular guest. Some examples of the interaction energy mapping are shown in Figure 5, for some classical metal-oxalate honeycomb structures.
Quite important MOFs are those which maximize the porosity, without losing chemical and mechanical stability. One famous example is MOF5  (Zn4O(1,3,5-benzenetribenzoate)2), which is often taken as a benchmark for its extraordinary porous properties. Civalleri et al.  have analyzed the electronic structure and the electrostatic potential φ(r) in MOF-5, demonstrating the feasibility of a fully ab initio simulation and the useful approach based on the electrostatic potential. It is worth pointing out that while the electron density threshold criterion only addresses the short-range repulsion, φ(r) maps only address the electrostatic attraction/repulsion of a point charge, therefore their interpretation should be considered very carefully.
Among metal framework structures, of particular interest are nowadays those carrying magnetic properties, in particular able to build up 1D, 2D or 3D magnetism. For this reason, ED studies on magnetic frameworks are becoming more common, sometime coupling the electron density with the spin density analysis.
Iversen et al. – have deeply investigated the connection between electron density topology and framework structures in magnetic materials, analyzing in particular the stereoelectronics at the metal site. In fact, the multipolar expansion enables to roughly estimate the electronic population of transition metal d-orbitals , therefore having at least a confirmation of the electronic structures predicted from the magnetic measurements. In some cases, it is even possible to estimate the magnetic moments ,, by careful consideration of the electronic population. However, the full information on the spin density is available only if using a radiation able to interact with the electronic spin, as for example polarized neutrons –. In particular, recent work by Deutsch et al.  opens the possibility to simultaneously determine the charge and spin density, from combination of X-ray diffraction and polarized neutron diffraction. Very important is in fact the possibility to determine the paths of magnetic exchange, in view of the growing interest for metal framework materials with tunable magnetism.
In this mini-review, we have focused on the application of charge density analysis in the vast but still growing field of crystal engineering, showing that crystallographic information, which goes beyond the determination of the molecular geometry or the supramolecular packing, can be of enormous impact for this kind of research.
The charge density is inherently connected with X-ray diffraction, given that it is in fact the electron distribution which is visualized after tentatively assigning the unknown phases to the measured structure factors. However, crystallographers rarely push the structure modeling much further and do not investigate the finer details of electron density in between atoms and molecules, although this is nowadays much more easily available thanks to more powerful radiation sources (including synchrotron radiation) and more sensitive and faster detectors. Pushing this limit would enable extracting more information on the robustness of a given synthon, of a given metal-ligand coordination or a host-guest interaction, which are quite important in crystal engineering.
In parallel to the experimental side, also the ab initio quantum chemical analysis of these materials made enormous progresses in the past decade, which enable calculation of a number of properties that are not straightforwardly available from measured scattering intensities.
In the future, one may expect that the whole process of material design and fabrication could be integrated, having as initial step predictions based on calculated or measured electron density models of the molecular building blocks, that could offer an initial guess of the expected property of a crystalline solid.
aIt is always important to remember that the electron density is a probability function to find one electron in space, therefore it is always positive, despite the charge of the electron is actually negative. Thus an excess of electron density in an atom means a negative charge.
Pepinsky R: Crystal Engineering: New Concepts in Crystallography. Phys Rev. 1955, 100: 971-
Schmidt GMJ: Photodimerization in the solid state. Pure Appl Chem. 1971, 27: 647-678.
Desiraju GR: Crystal engineering: A brief overview. J Chem Sci. 2010, 122: 667-675.
Figgis BN, Hitchman MA: Ligand Field Theory and Its Applications. 2000, Wiley-VCH, New York
Jeffrey GA: An Introduction to Hydrogen Bonding. 1997, Oxford University Press, Oxford
Peresypkina EV, Blatov VA: Topology of molecular packings in organic crystals. Acta Cryst. 2000, B56: 1035-1045.
Debye P: Zerstreuung von Röntgenstrahlen. Ann Phys. 1915, 46: 809-823.
Coppens P: Comparative X-Ray and Neutron Diffraction Study of Bonding Effects in s-Triazine. Science. 1967, 158: 1577-1579.
Coppens P: X-Ray charge densities and chemical bonding. 1997, Oxford University Press, Oxford
Bader RFW: Atoms in Molecules: a quantum theory. 1990, Oxford University Press, Oxford
Macchi P, Sironi A: Chemical bonding in transition metal carbonyl clusters: complementary analysis of theoretical and experimental electron densities. Coord Chem Rev. 2003, 238: 383-412.
Stalke D: Meaningful Structural Descriptors from Charge Density. Chem. Eur J. 2011, 17: 9264-9278.
Flaig R, Koritsanszky T, Dittrich B, Wagner A, Luger P: Intra- and intermolecular topological properties of amino acids:A comparative study of experimental and theoretical results. J Am Chem Soc. 2002, 124: 3407-3417.
Kubicki M, Borowiak T, Dutkiewicz G, Souhassou M, Jelsch C, Lecomte C: Experimental electron density of 1-phenyl-4-nitroimidazole and topological analysis of C––H–-O and C––H–-N hydrogen bonds. J Phys Chem B. 2002, 106: 3706-3714.
May E, Destro R, Gatti C: The unexpected and large enhancement of the dipole moment in the 3,4-bis(dimethylamino)-3-cyclobutene-1,2-dione (DMACB) molecule upon crystallization: A new role of the intermolecular CH–-O interactions. J Am Chem Soc. 2001, 123: 12248-12254.
Gatti C, May E, Destro R, Cargnoni F: Fundamental properties and nature of C-H–-O interactions in crystals on the basis of experimental and theoretical charge densities. The case of 3,4-bis(dimethylamino)-3-cyclobutene-1,2-dione (DMACB) crystal. J Phys Chem A. 2002, 106: 2707-2720.
Tsirelson VG, Ozerov RP: Electron Densities and Bonding in crystals. 1996, Institute of Physics, Bristol/Philadelphia
Koritsanszky TS, Coppens P: Chemical applications of X-ray charge-density analysis. Chem Rev. 2001, 101: 1583-1627.
Macchi P: Modern charge density studies: The entanglement of experiment and theory Cryst. Cryst. Rev. 2013, 19: 58-101.
Stewart RF: Electron Population Analysis with Rigid Pseudoatoms. Acta Cryst. 1976, A32: 565-574.
Hansen NK, Coppens P: Electron population analysis of accurate diffraction data. 6 Testing aspherical atom refinements on small-molecule data sets. Acta Cryst. 1978, A34: 909-921.
Hirshfeld FL: Bonded-atom fragments for describing molecular charge-densities. Theor Chim Acta. 1977, 44: 129-138.
Spackman MA, Byrom PG: A novel definition of a molecule in a crystal. Chem Phys Lett. 1997, 267: 215-220.
McKinnon JJ, Spackmam MA, Mitchell AS: Novel tools for visualizing and exploring intermolecular interactions on molecular crystals. Acta CrystB. 2004, 60: 627-668.
McKinnon JJ: Jayatilaka D. 2007, Towards quantitative analysis of intermolecular interactions with Hirshfeld surfaces Chem Commun, Spackman MA
Spackman MA, Jayatilaka D: Hirshfeld Surface Analysis. CrystEngComm. 2009, 11: 19-32.
Gillet JM, Koritsanszky T: Past, Present and Future of Charge Density and Density Matrix Refinements. Modern Charge Density Analysis. Edited by: Gatti C, Macchi P. 2012, Springer, New York
Bader RFW: Why are there atoms in chemistry?. Can J Chem. 1998, 76: 973-988.
Bader RFW, MacDougall PJ, Lau CDH: Bonded and nonbonded charge concentrations and their relation to molecular geometry and reactivity. J Am Chem Soc. 1984, 106: 1594-1605.
Macchi P, Proserpio DM, Sironi A: Experimental Electron Density in a Transition Metal Dimer: Metal-Metal and Metal-Ligand Bonds. J Am Chem Soc. 1998, 120: 13429-13435.
Cremer D, Kraka E: Chemical Bonds without Bonding Electron Density – Does the Difference Electron-Density Analysis Suffice for a Description of the Chemical Bond?. Angew Chem. 1984, 23: 627-628.
Bader RFW, Stephens ME: Spatial localization of the electronic pair and number distributions in molecules. J Am Chem Soc. 1975, 97: 7391-7399.
Koch U, Popelier P: Characterization of C-H-O Hydrogen-Bonds on the Basis of the Charge-Density. J Physical Chem. 1995, 99: 9747-9754.
Popelier PLA: Atoms in Molecules. An Introduction. 2000, Pearson Education, Harlow, Great Britain
Espinosa E, Molins E, Lecomte C: Hydrogen bond strengths revealed by topological analyses of experimentally observed electron densities. Chem Phys Lett. 1998, 285: 170-173.
Abramov YA: On the Possibility of Kinetic Energy Density Evaluation from the Experimental Electron-Density Distribution. Acta Cryst A. 1997, 53: 264-272.
Spackman MA: Hydrogen bond energetics from topological analysis of experimental electron densities: Recognising the importance of the promolecule. Chem Phys Lett. 1999, 301: 425-429.
Munshi P, Guru Row TN: Intra- and intermolecular interactions in small bioactive molecules: cooperative features from experimental and theoretical charge-density analysis. Acta Cryst. B . 2006, 62: 612-626.
Munshi P, Guru Row TN: Topological Analysis of Charge Density Distribution in Concominant Polymorphs of 3-Acetylcoumarin, A Case of Packing Polymorphism. Cryst Growth Des. 2006, 6: 708-718.
Srinivasa Gopalan R, Kumaradhas P, Kulkarni GU, Rao CNR: An experimental charge density study of aliphatic dicarboxylic acids. J Mol Structure. 2000, 521: 97-106.
Howard JAK, Mahon MF, Raithby PR, Sparkes HA: Trans-cinnamic acid and coumarin-3-carboxylic acid: experimental charge-density studies to shed light on [2 + 2] cycloaddition reactions. Acta Cryst B. 2009, 65: 230-237.
Bushmarinov IS, Nabiev OG, Kostyanovsky RG, Antipin MY, Lyssenko KA: The azide anion as a building block in crystal engineering from a charge density point of view. Cryst Eng Comm. 2011, 13: 2930-2934.
Becke AD, Edgecombe KE: A simple measure of electron localization in atomic and molecular systems. J Chem Phys. 1990, 92: 5397-5403.
Tsirelson V, Stash A: Determination of the electron localization function from electron density Chem . Phys Lett. 2002, 351: 142-148.
Gopalan RS, Kulkarni GU, Rao CNR: An Experimental Charge Density Study of the Effect of the Noncentric Crystal Field on the Molecular Properties of Organic NLO Materials. ChemPhysChem. 2000, 1: 127-135.
Hassel O, Romming C: Direct structural evidence for weak charge-transfer bonds in solids containing chemically saturated molecules. Q Rev Chem Soc. 1962, 16: 1-18.
Hassel O: Structural aspects of interatomic charge-transfer bonding. Science. 1970, 170: 497-502.
Legon AC: Prereactive Complexes of Dihalogens XY with Lewis Bases B in the Gas Phase: A Systematic Case for the Halogen Analogue B–-XY of the Hydrogen Bond B–-HX. Angew Chem Int Ed Eng. 1999, 38: 2686-2714.
Bianchi R, Forni A, Pilati T: The Experimental Electron Density Distribution in the Complex of (E)-1,2-Bis(4-pyridyl)ethylene with 1,4-Diiodotetrafluorobenzene at 90K. Chem Eur J. 2003, 9: 1631-1638.
Bui TTT, Dahoui S, Lecomte C, Desiraju GR, Espinosa E: The Nature of Halogen · · · Halogen Interactions: A Model Derived from Experimental Charge-Density Analysis. Angew Chem Int Ed Eng. 2009, 48: 3838-3841.
Brinck T, Murray JS, Politzer P: Surface electrostatic potentials of halogenated methanes as indicators of directional intermolecular interactions. Int J Quantum Chem. 1992, 44: 57-64.
Brinck T, Murray JS, Politzer P: Molecular-surface electrostatic potentials and local ionization energies of group-v-vii hydrides and their anions - relationships for aqueous and gas-phase acidities. Int J Quantum Chem. 1993, 48: 73-88.
Nelyubina YV, Antipin MY, Dunin DS, Kotov VY, Lyssenko KA: Unexpected “amphoteric” character of the halogen bond: The charge density study of the co-crystal of N-methylpyrazine iodide with I2. Chem Comm. 2010, 46: 5325-5327.
Brezgunova ME, Aubert E, Dahaoui S, Fertey P, Lebègue S, Jelsch C, Ángyán JG, Espinosa E: Charge density analysis and topological properties of Hal 3-synthons and their comparison with competing hydrogen bonds. Cryst Growth Des. 2012, 12: 5373-5386.
Pavan MS, Durga Prasad K, Guru Row TN: Halogen bonding in fluorine: Experimental charge density study on intermolecular F–-F and F–-S donor-acceptor contacts. Chem Comm. 2013, 49: 7558-7560.
Stone AJ: Are Halogen Bonded Structures Electrostatically Driven?. J Am Chem Soc. 2013, 135: 7005-7009.
Spackman MA: Rationalizing Molecular Crystal Structures using Hirshfeld Surfaces. 2014, XXIII Congress and General Assembly of the International Union of Crystallography, Montreal
Spackman MA, McKinnon JJ: Fingerprinting Intermolecular Interactions in Molecular Crystals. CrystEngComm. 2002, 4: 378-392.
Wolff SK, Grimwood DJ, McKinnon JJ, Turner MJ, Jayatilaka D, Spackman MA: CrystalExplorer (Version 3.1). 2012, University of Western Australia, Perth
Johnson ER, Keinan S, Mori-Sánchez P, Contreras-García J, Cohen AJ, Yan W: Revealing Noncovalent Interactions. J Am Chem Soc. 2010, 132: 6498-6506.
Contreras-García J, Johnson ER, Keinan S, Chaudret R, Piquemal J-P, Beratan DN, Yang W: NCIPLOT: A Program for Plotting Noncovalent Interaction Regions. J Chem Theory Comput. 2011, 7: 625-632.
Saleh G, Gatti C, Lo Presti L, Contreras-García J: Revealing non-covalent Interactions in Molecular Crystals through Their Experimental Electron densities. Chem Eur J. 2012, 18: 15523-15536.
Hey J, Leusser D, Kratzert D, Fliegl H, Dieterich JM, Mata RA, Stalke D: Heteroaromaticity approached by charge density investigations and electronic structure calculations. Phys Chem Chem Phys. 2013, 15: 20600-20610.
Krawczuk A, Pérez D, Macchi P: PolaBer: a program to calculate and visualize distributed atomic polarizabilities based on electron density partitioning. J Appl Cryst. 2014, 47: 1452-1458.
Matta CF, Bader RFW: An Atoms-In-Molecules study of the genetically-encoded amino acids: I. Effects of conformation and of tautomerization on geometric, atomic, and bond properties. Proteins. 2000, 40: 310-329.
Krawczuk-Pantula A, Pérez D, Macchi P: Distributed atomic polarizabilities from electron density. 1. Motivations and Theory. Trans Amer Cryst Ass. 2012, 42: 1-25.
Keith TA: Atomic Response Properties. The Quantum Theory of Atoms in Molecules: from Solid State to DNA and Drug Design. Edited by: Matta CF, Boyd RJ. 2007, Viley-VCH, Weinheim
Lipinski CA: Pharmaceutical Profiling in Drug Discovery for Lead Selection. 2004, AAPS Press Arlington, USA
Trask AV, Motherwell WDS, Jones W: Pharmaceutical Cocrystalization: Engineering a Remedy for caffeine Hydration. Crystal Growth Design. 2005, 5: 1013-1021.
Aakeröy CB, Forbes S, Desper J: Using Cocrystals To Systematically Modulate Aqueous Solubility and Melting Behavior of an Anticancer Drug. J Am Chem Soc. 2009, 131: 17048-17049.
Steed JW: The role of co-crystals in pharmaceutical design. Trends Pharmacol Sci. 2013, 34: 185-193.
Hathwar VR, Pal R, Guru Row TN: Charge Density Analysis of Crystals of Nicotinamide with Salicylic Acid and Oxalic Acid: An Insight into the Salt to Cocrystal Continuum. Cryst Growth Des. 2010, 10: 3306-3310.
Hathwar VR, Thakur TS, Guru Row TN, Desiraju GR: Transferability of multipole charge density parameters for supramolecular synthons: A new tool for quantitative crystal engineering. Cryst Growth Des. 2011, 11: 616-623.
Volkov A, Messerschmidt M, Coppens P: Improving the scattering-factor formalism in protein refinement: application of the University at Buffalo Aspherical-Atom Databank to polypeptide structures. Acta Cryst D. 2007, 63: 160-170.
Volkov A, Li X, Koritsanszky T, Coppens P: Ab Initio Quality Electrostatic Atomic and Molecular Properties Including Intermolecular Energies from a Transferable Theoretical Pseudoatom Databank J . Phys Chem A. 2004, 108: 4283-4300.
Dominiak PM, Volkov A, Li X, Messerschmidt M, Coppens P: A Theoretical Databank of Transferable Aspherical Atoms and Its Application to Electrostatic Interaction Energy Calculations of Macromolecules J . Chem Theory Comput. 2007, 3: 232-247.
Dittrich B, Hubschle CB, Luger P, Spackman MA: Introduction and validation of an invariom database for amino-acid, peptide and protein molecules. Acta Cryst D. 2006, 62: 1325-1335.
Hubschle CB, Luger P, Dittrich BJ: Automation of invariom and of experimental charge density modelling of organic molecules with the preprocessor program InvariomTool J . Appl Cryst. 2007, 40: 623-627.
Dittrich B, Strumpel M, Spackman MA, Koritsanszky T: Invarioms for improved absolute structure determination of light-atom crystal structures. Acta Cryst A. 2006, 62: 217-223.
Zarychta B, Pichon-Pesme V, Guillot B, Lecomte C, Jelsch C: On the application of an experimental multipolar pseudo-atom library for accurate refinement of small-molecule and protein crystal structures. Acta CrystA. 2007, 63: 108-125.
Jelsch C, Pichon-Pesme V, Lecomte C, Aubry A: Transferability of Multipole Charge-Density Parameters: Application to Very High Resolution Oligopeptide and Protein Structures Acta Cryst . Acta Cryst D. 1998, 54: 1306-1318.
Guillot B, Jelsch C, Podjarny A, Lecomte C: Charge-density analysis of a protein structure at subatomic resolution: the human aldose reductase case. Acta Cryst D. 2008, 64: 567-588.
Hathwar VR, Thakur TS, Dubey R, Pavan MS, Guru Row TN, Desiraju GR: Extending the supramolecular synthon based fragment approach (SBFA) for transferability of multipole charge density parameters to monofluorobenzoic acids and their co-crystals with isonicotinamide: importance of C–H · · · O, C–H · · · F and F · · · F intermolecular regions. J Phys Chem A. 2011, 115: 12852-12863.
Gryl M, Krawczuk-Pantula A, Stadnicka K: Charge-density analysis in polymorphs of urea-barbituric acid co-crystals Acta Cryst . Acta Cryst. B. 2011, 67: 144-154.
Gryl M, Krawczuk A, Stadnicka K: Polymorphism of urea–barbituric acid co-crystals. Acta Cryst B. 2008, 64: 623-632.
Schmidtmann M, Farrugia LJ, Middlemiss DS, Gutmann MJ, McIntyre GJ, Wilson CC: Experimental and Theoretical Charge Density Study of Polymorphic isonicotinamide-Oxalic Acid Molecular Complexes with Strong O…H…N Hydrogen Bonds. J Phys Chem A. 2009, 113: 13985-13997.
Dubey R, Pavan MS, Guru Row TN, Desiraju GR: Crystal landscape in the orcinol: 4,4′-bipyridine system: synthon modularity, polymorphism and transferability of multipole charge density parameters. IUCrJ. 2014, 1: 8-18.
Mossotti OF: Memorie di Matematica e di Fisica della Società Italiana delle Scienze Residente in Modena. 1850, 24:49–74.
Clausius R: Die mechanische Behandlung der Electricität. 1858, Vieweg, Braunschweig
Dunmur DA: The local electric field in anisotropic molecular crystals. Mol Physics. 1972, 22: 109-115.
Fkyerat A, Guelzim A, Baert F, Zyss J, Perigaud A: Assessment of the polarizabilities of a nonlinear optical compound [N-(4-nitrophenyl)-(L)-prolinol] from an experimental electronic density study. Phys Rev. 1996, 53: 16236-16246.
Hamzaoui F, Zanoun A, Vergoten G: The molecular linear polarizability from X-ray diffraction study. The case of 3-methyl 4-nitropyridine N-oxide (POM). J Mol Struct. 2004, 697: 17-22.
Chouaih A, Hamzaoui F, Vergoten G: Capability of X-ray diffraction to the determination of the macroscopic linear susceptibility in a crystalline environment: the case of 3-Methyl 4-Nitropyridine N-oxide (POM) J . Mol Struct. 2005, 738: 33-38.
Robinson FNH: Nonlinear optical coefficients. Bell Syst Tech J. 1967, 46: 913-956.
Whitten AE, Jayatilaka D, Spackman MA: Effective molecular polarizabilities and crystal refractive indices estimated from x-ray diffraction data. J Chem Phys. 2006, 125: 174505-174514.
Grimwood DJ, Jayatilaka D: Wavefunctions derived from experiment. II. A wavefunction for oxalic acid dihydrate. Acta Cryst A. 2001, 57: 87-100.
Jayatilaka D, Grimewood DJ, Lee A, Lemay A, Russell AJ, Taylor C, Wolff SK, Cassam-Chenai P, Whitten A: TONTO, a system for computational chemistry. 2005
Hickstein DD, Cole JM, Turner MJ, Jayatilaka D: Modeling electron density distributions from X-ray diffraction to derive optical properies: Constrained wavefunction versus multipole refinement. J Chem Phys. 2013, 139: 064108-064114.
Genoni A: Molecular orbitals Strictly Localized on Small Molecular Fragments from X-ray Diffraction Data. J Phys Chem Lett. 2013, 4: 1093-1099.
Chimpri AS, Gryl M, Dos Santos LHR, Krawczuk A: Correlation between accurate Electron density and Linear Optical Properties in Amino acid derivatives: L-Histidinium Hydrogen Oxalate. Crystal Growth Design. 2013, 13: 2995-3010.
Yaghi OM, Li G, Li H: Selective binding and removal of gusts in a microporous metal-organic framework. Nature. 1995, 378: 703-706.
Spackman MA, Jayatilaka D: Hirshfeld surface analysis. CrystEngComm. 2009, 11: 19-32.
Filsø MO, Turner MJ, Gibbs GV, Adams S, Spackman MA, Iversen BB: Visualizing Lithium-Ion Migration Pathways in Battery Materials. Chemistry . Eur J. 2013, 19: 15535-15544.
Wheatley RJ: An overlap model for exchange-induction: application to alkali halides. Chem Phys Lett. 1998, 294: 487-492.
Spackman MA: Atom-Atom potentials via electron gas theory. J Chem Phys. 1986, 85: 6579-6586.
Spackman MA: A simple quantitative model of hydrogen bonding . J Chem Phys. 1986, 85: 6587-6601.
Spackman MA: A simple quantitative model of hydrogen bonding. Applications to more complex systems . J Phys Chem. 1987, 91: 3179-3186.
Spackman MA: The use of the promolecular charge density to approximate the penetration contribution to intermolecular electrostatic energies. Chem Phys Lett. 2006, 418: 158-162.
Buckingham AD: Intermolecular Interactions: From Diatomics to Biopolymers. 1978, Wiley, Winchester
McWeeny R: Methods of molecular quantum mechanics. 1989, Academic Press, London
Coppens P, Abramov Y, Carducci M, Korjov B, Novozhilova I, Alhambra C, Pressprich MR: Experimental Charge Densities and Intermolecular Interactions: Electrostatic and Topological Analysis of DL-Histidine. J Am Chem Soc. 1999, 121: 2585-2593.
Lennard-Jones JE: On the determination of molecular fields. II From the equation of state of a gas. Proc R Soc Lond A. 1924, 106: 463-477.
Gavezzotti A: Calculation of Intermolecular Interaction Energies by Direct Numerical Integration over Electron Densities. I. Electrostatic and Polarization Energies in Molecular Crystals. J Phys Chem B. 2002, 106: 4145-4154.
Volkov A, Koritsanszky T, Coppens P: Combination of the exact potential and multipole methods (EP/MM) for evaluation of intermolecular electrostatic interaction energies with pseudoatom representation of molecular electron densities. Chem Phys Lett. 2004, 391: 170-175.
Volkov A, King H, Coppens P: Dependence of the Intermolecular Electrostatic Interaction Energy on the Level of Theory and the Basis Set. J Chem Theor Comput. 2006, 2: 81-89.
Gavezzotti A: Calculation of Intermolecular Interaction Energies by Direct Numerical Integration over Electron Densities. 2. An Improved Polarization Model and the Evaluation of Dispersion and Repulsion Energies. J Phys Chem B. 2003, 107: 2344-2353.
Komorowski L, Lipinski J, Szarek P: Polarization justified Fukui functions. J Chem Phys. 2009, 131: 124120-
London F: On the Theory and Systematic of Molecular Forces. Z Physik. 1930, 63: 245-279.
Adams S, Rao RP: Transport pathways for mobile ions in disordered solids from the analysis of energy-scaled bond-valence mismatch landscapes. Phys Chem Chem Phys. 2009, 11: 3210-3216.
Chimpri AS, Macchi P: Electron density building block approach for metal organic frameworks. Phys Scr. 2013, 87: 048105-048109.
Li H, Eddaoudi M, O’Keeffe M, Yaghi O: Design and synthesis of an exceptionally stable and highly porous metal-organic framework. Nature. 1999, 402: 276-279.
Civalleri B, Napoli F, Nöel Y, Roetti C, Dovesi R: Ab-initio prediction of materials properties with CRYSTAL: MOF-5 as a case study. CrystEngComm. 2006, 8: 364-371.
Poulsen RD, Bentien A, Graber T, Iversen BB: Synchrotron charge-density studies in materials chemistry: 16 K X-ray charge density of a new magnetic metal-organic framework material [Mn2(C8H4O4)2(C3H7NO)2]. Acta Cryst. 2004, A60: 382-389.
Poulsen RD, Bentien A, Chevalier M, Iversen BB: Synthesis, physical properties, multitemperature crystal structure, and 20 K synchrotron X-ray charge density of a magnetic metal organic framework structure [Mn2(C8H4O4)2(C3H7NO)2]. J Am Chem Soc. 2005, 127: 9156-9166.
Poulsen RD, Jørgensen MRV, Overgaard J, Larsen FK, Morgenroth W, Graber T, Chen YS, Iversen BB: Synchrotron X-Ray Charge-Density Study of Coordination Polymer [Mn(HCOOH) 2(H 2O) 2] ∞ Chemistry . Eur J. 2007, 13: 9775-9979.
Holladay A, Leung P, Coppens P: Generalized relations between d-orbital occupancies of transition metal atoms and electron density multipole population parameters from X-ray diffraction data. Acta Cryst. 1983, A39: 377-387.
Jauch W, Reehuis M: Electron density distribution in paramagnetic and antiferromagnetic CoO: A γ-ray diffraction study. Phys Rev B. 2002, 65: 125111-125118.
Jauch W, Reehuis M: Electron density distribution in paramagnetic and antiferromagnetic MnO: A γ-ray diffraction study. Phys Rev B. 2003, 67: 184420-184428.
Schweizer J: Polarized neutrons and polarization analysis. Neutron scattering from magnetic crystals. Edited by: Chatterji T. 2006, Elsevier, Amsterdam, 153-213.
Pillet S, Souhassou M, Pontillon Y, Caneschi A, Gatteschi D, Lecomte C: Investigation of magnetic interaction pathways by experimental electron density measurements: application to an organic free radical, p(methylthio)phenyl nitronyl nitroxide. New J Chem. 2001, 25: 131-143.
Pillet S, Souhassou M, Mathonière C, Lecomte C: Electron density distribution of an oxamato bridged Mn(II)-Cu(II) bimetallic chain and correlation to magnetic properties. J Am Chem Soc. 2004, 126: 1219-1228.
Claiser N, Souhassou M, Lecomte C, Gillon B, Carbonera C, Caneschi A, Dei A, Gatteschi D, Bencini A, Pontillon Y, Lelièvre-Berna E: Combined charge and spin density experimental study of the yttrium(III) semiquinonato complex Y(HBPz3)2(DTBSQ) and DFT calculations. J Phys Chem B. 2005, 109: 2723-2732.
Deutsch M, Gillon B, Claiser N, Gillet JM, Lecomte C, Souhassou M: First spin-resolved electron distributions in crystals from combined polarized neutron and X-ray diffraction experiments. IUCr J. 2014, 1: 194-199.
Acknowledgements and funding
We thank the Swiss national science foundation for financial support (project 200020–141271).
The authors declare that they have no competing interests.
Both authors have equally contributed in writing this review and have given final approval of the version to be published.