- Open Access
Charge density analysis for crystal engineering
© Krawczuk and Macchi; licensee Chemistry Central Ltd. 2014
- Received: 10 July 2014
- Accepted: 30 October 2014
- Published: 16 December 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.
- Charge density analysis
- Crystal engineering
- Supramolecular chemistry
- X-ray diffraction
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.
Overview of the most important interactions occurring between two c losed-shell electron density distributions (R is the distance between the two centers of masses)
Coulomb attraction/repulsion between unperturbed electron densities
Long range especially monopolar charges of ions (÷1/RL+1; where L is the sum of the multipole orders; L = 0 for charge-charge interactions)
Monopole-monopole interactions are not directional; increasing directionality for higher multipolar moments
Stabilization or destabilization, depending on the sign and orientation of the electrostatic moments of the interacting systems
Coulomb attraction between electron density of one molecule and field induced polarization of the other
Shorter range (÷1/R4)
Coulomb attraction between mutually polarized electron densities
Quite short range (÷1/R6)
Short range repulsion
The reduced probability of having two electrons with the same spin very close to each other (Fermi-hole)
Very short range (÷1/R12 or exponential)
Interaction between frontier orbitals of the interacting systems it implies partial covalence
Occurs only for contacts shorter than van der Waals distances
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.).
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.
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.
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.
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 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 .
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.
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.
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.
Acknowledgements and funding
We thank the Swiss national science foundation for financial support (project 200020–141271).
- Pepinsky R: Crystal Engineering: New Concepts in Crystallography. Phys Rev. 1955, 100: 971-Google Scholar
- Schmidt GMJ: Photodimerization in the solid state. Pure Appl Chem. 1971, 27: 647-678.Google Scholar
- Desiraju GR: Crystal engineering: A brief overview. J Chem Sci. 2010, 122: 667-675.Google Scholar
- Figgis BN, Hitchman MA: Ligand Field Theory and Its Applications. 2000, Wiley-VCH, New YorkGoogle Scholar
- Jeffrey GA: An Introduction to Hydrogen Bonding. 1997, Oxford University Press, OxfordGoogle Scholar
- Peresypkina EV, Blatov VA: Topology of molecular packings in organic crystals. Acta Cryst. 2000, B56: 1035-1045.Google Scholar
- Debye P: Zerstreuung von Röntgenstrahlen. Ann Phys. 1915, 46: 809-823.Google Scholar
- Coppens P: Comparative X-Ray and Neutron Diffraction Study of Bonding Effects in s-Triazine. Science. 1967, 158: 1577-1579.Google Scholar
- Coppens P: X-Ray charge densities and chemical bonding. 1997, Oxford University Press, OxfordGoogle Scholar
- Bader RFW: Atoms in Molecules: a quantum theory. 1990, Oxford University Press, OxfordGoogle Scholar
- 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.Google Scholar
- Stalke D: Meaningful Structural Descriptors from Charge Density. Chem. Eur J. 2011, 17: 9264-9278.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- Tsirelson VG, Ozerov RP: Electron Densities and Bonding in crystals. 1996, Institute of Physics, Bristol/PhiladelphiaGoogle Scholar
- Koritsanszky TS, Coppens P: Chemical applications of X-ray charge-density analysis. Chem Rev. 2001, 101: 1583-1627.Google Scholar
- Macchi P: Modern charge density studies: The entanglement of experiment and theory Cryst. Cryst. Rev. 2013, 19: 58-101.Google Scholar
- Stewart RF: Electron Population Analysis with Rigid Pseudoatoms. Acta Cryst. 1976, A32: 565-574.Google Scholar
- 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.Google Scholar
- Hirshfeld FL: Bonded-atom fragments for describing molecular charge-densities. Theor Chim Acta. 1977, 44: 129-138.Google Scholar
- Spackman MA, Byrom PG: A novel definition of a molecule in a crystal. Chem Phys Lett. 1997, 267: 215-220.Google Scholar
- McKinnon JJ, Spackmam MA, Mitchell AS: Novel tools for visualizing and exploring intermolecular interactions on molecular crystals. Acta CrystB. 2004, 60: 627-668.Google Scholar
- McKinnon JJ: Jayatilaka D. 2007, Towards quantitative analysis of intermolecular interactions with Hirshfeld surfaces Chem Commun, Spackman MAGoogle Scholar
- Spackman MA, Jayatilaka D: Hirshfeld Surface Analysis. CrystEngComm. 2009, 11: 19-32.Google Scholar
- 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 YorkGoogle Scholar
- Bader RFW: Why are there atoms in chemistry?. Can J Chem. 1998, 76: 973-988.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- Bader RFW, Stephens ME: Spatial localization of the electronic pair and number distributions in molecules. J Am Chem Soc. 1975, 97: 7391-7399.Google Scholar
- 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.Google Scholar
- Popelier PLA: Atoms in Molecules. An Introduction. 2000, Pearson Education, Harlow, Great BritainGoogle Scholar
- 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.Google Scholar
- Abramov YA: On the Possibility of Kinetic Energy Density Evaluation from the Experimental Electron-Density Distribution. Acta Cryst A. 1997, 53: 264-272.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- Becke AD, Edgecombe KE: A simple measure of electron localization in atomic and molecular systems. J Chem Phys. 1990, 92: 5397-5403.Google Scholar
- Tsirelson V, Stash A: Determination of the electron localization function from electron density Chem . Phys Lett. 2002, 351: 142-148.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- Hassel O: Structural aspects of interatomic charge-transfer bonding. Science. 1970, 170: 497-502.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- Stone AJ: Are Halogen Bonded Structures Electrostatically Driven?. J Am Chem Soc. 2013, 135: 7005-7009.Google Scholar
- Spackman MA: Rationalizing Molecular Crystal Structures using Hirshfeld Surfaces. 2014, XXIII Congress and General Assembly of the International Union of Crystallography, MontrealGoogle Scholar
- Spackman MA, McKinnon JJ: Fingerprinting Intermolecular Interactions in Molecular Crystals. CrystEngComm. 2002, 4: 378-392.Google Scholar
- Wolff SK, Grimwood DJ, McKinnon JJ, Turner MJ, Jayatilaka D, Spackman MA: CrystalExplorer (Version 3.1). 2012, University of Western Australia, PerthGoogle Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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, WeinheimGoogle Scholar
- Lipinski CA: Pharmaceutical Profiling in Drug Discovery for Lead Selection. 2004, AAPS Press Arlington, USAGoogle Scholar
- Trask AV, Motherwell WDS, Jones W: Pharmaceutical Cocrystalization: Engineering a Remedy for caffeine Hydration. Crystal Growth Design. 2005, 5: 1013-1021.Google Scholar
- 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.Google Scholar
- Steed JW: The role of co-crystals in pharmaceutical design. Trends Pharmacol Sci. 2013, 34: 185-193.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- Gryl M, Krawczuk A, Stadnicka K: Polymorphism of urea–barbituric acid co-crystals. Acta Cryst B. 2008, 64: 623-632.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- Mossotti OF: Memorie di Matematica e di Fisica della Società Italiana delle Scienze Residente in Modena. 1850, 24:49–74.Google Scholar
- Clausius R: Die mechanische Behandlung der Electricität. 1858, Vieweg, BraunschweigGoogle Scholar
- Dunmur DA: The local electric field in anisotropic molecular crystals. Mol Physics. 1972, 22: 109-115.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- Robinson FNH: Nonlinear optical coefficients. Bell Syst Tech J. 1967, 46: 913-956.Google Scholar
- 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.Google Scholar
- Grimwood DJ, Jayatilaka D: Wavefunctions derived from experiment. II. A wavefunction for oxalic acid dihydrate. Acta Cryst A. 2001, 57: 87-100.Google Scholar
- 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. 2005Google Scholar
- 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.Google Scholar
- Genoni A: Molecular orbitals Strictly Localized on Small Molecular Fragments from X-ray Diffraction Data. J Phys Chem Lett. 2013, 4: 1093-1099.Google Scholar
- 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.Google Scholar
- Yaghi OM, Li G, Li H: Selective binding and removal of gusts in a microporous metal-organic framework. Nature. 1995, 378: 703-706.Google Scholar
- Spackman MA, Jayatilaka D: Hirshfeld surface analysis. CrystEngComm. 2009, 11: 19-32.Google Scholar
- 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.Google Scholar
- Wheatley RJ: An overlap model for exchange-induction: application to alkali halides. Chem Phys Lett. 1998, 294: 487-492.Google Scholar
- Spackman MA: Atom-Atom potentials via electron gas theory. J Chem Phys. 1986, 85: 6579-6586.Google Scholar
- Spackman MA: A simple quantitative model of hydrogen bonding . J Chem Phys. 1986, 85: 6587-6601.Google Scholar
- Spackman MA: A simple quantitative model of hydrogen bonding. Applications to more complex systems . J Phys Chem. 1987, 91: 3179-3186.Google Scholar
- 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.Google Scholar
- Buckingham AD: Intermolecular Interactions: From Diatomics to Biopolymers. 1978, Wiley, WinchesterGoogle Scholar
- McWeeny R: Methods of molecular quantum mechanics. 1989, Academic Press, LondonGoogle Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- Komorowski L, Lipinski J, Szarek P: Polarization justified Fukui functions. J Chem Phys. 2009, 131: 124120-Google Scholar
- London F: On the Theory and Systematic of Molecular Forces. Z Physik. 1930, 63: 245-279.Google Scholar
- 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.Google Scholar
- Chimpri AS, Macchi P: Electron density building block approach for metal organic frameworks. Phys Scr. 2013, 87: 048105-048109.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- Jauch W, Reehuis M: Electron density distribution in paramagnetic and antiferromagnetic CoO: A γ-ray diffraction study. Phys Rev B. 2002, 65: 125111-125118.Google Scholar
- Jauch W, Reehuis M: Electron density distribution in paramagnetic and antiferromagnetic MnO: A γ-ray diffraction study. Phys Rev B. 2003, 67: 184420-184428.Google Scholar
- Schweizer J: Polarized neutrons and polarization analysis. Neutron scattering from magnetic crystals. Edited by: Chatterji T. 2006, Elsevier, Amsterdam, 153-213.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.