Is there a future for topological analysis in experimental charge-density research?

Understanding Hygroscopicity of Theophylline via a Novel Cocrystal Polymorph: A Charge Density Study

The charge density distribution in a novel cocrystal (1) complex of 1,3-dimethylxanthine (theophylline) and propanedioic acid (malonic acid) has been determined. The molecules crystallize in the triclinic, centrosymmetric space group P1̅, with four independent molecules (Z = 4) in the asymmetric unit (two molecules each of theophylline and malonic acid). Theophylline has a notably high hygroscopic nature, and numerous cocrystals have shown a significant improvement in stability to humidity. A charge density study of the novel polymorph has identified interesting theoretical results correlating the stability enhancement of theophylline via cocrystallization. Topological analysis of the electron density highlighted key differences (up to 17.8) in Laplacian (∇²ρ) between the experimental (EXP) and single-point (SP) models, mainly around intermolecular-bonded carbonyls. Further investigation via molecular electrostatic potential maps reaffirmed that the charge redistribution enhanced intramolecular hydrogen bonding, predominantly for N(2') and N(2) (61.2 and 61.8 kJ mol^-1, respectively). An overall weaker lattice energy of the triclinic form (-126.1 kJ mol^-1) compared to that of the monoclinic form (-133.8 kJ mol^-1) suggests a lower energy threshold to overcome to initiate dissociation. Future work via physical testing of the novel cocrystal in both dissolution and solubility will further solidify the correlation between theoretical and experimental results.

Topological analysis of chemical bonding in the layered FePSe3 upon pressure-induced phase transitions

Two pressure-induced phase transitions have been theoretically studied in the layered iron phosphorus triselenide (FePSe₃ ). Topological analysis of chemical bonding in FePSe₃ has been performed based on the results of first-principles calculations within the periodic linear combination of atomic orbitals (LCAO) method with hybrid Hartree-Fock-DFT B3LYP functional. The first transition at about 6 GPa is accompanied by the symmetry change from R 3 ¯ to C2/m, whereas the semiconductor-to-metal transition (SMT) occurs at about 13 GPa leading to the symmetry change from C2/m to P 3 ¯ 1 m . We found that the collapse of the band gap at about 13 GPa occurs due to changes in the electronic structure of FePSe₃ induced by relative displacements of phosphorus or selenium atoms along the c-axis direction under pressure. The results of the topological analysis of the electron density and its Laplacian demonstrate that the pressure changes not only the interatomic distances but also the bond nature between the intralayer and interlayer phosphorus atoms. The interlayer P-P interactions are absent in two non-metallic FePSe₃ phases while after SMT the intralayer P-P interactions weaken and the interlayer P-P interactions appear.

Identifying intermolecular atom⋯atom interactions that are not just bonding but also competitive

This highlight criticises the QTAIM method and discusses algorithms for identifying intermolecular interactions that are both bonding and competitive.

A Million Crystal Structures: The Whole Is Greater than the Sum of Its Parts

The founding in 1965 of what is now called the Cambridge Structural Database (CSD) has reaped dividends in numerous and diverse areas of chemical research. Each of the million or so crystal structures in the database was solved for its own particular reason, but collected together, the structures can be reused to address a multitude of new problems. In this Review, which is focused mainly on the last 10 years, we chronicle the contribution of the CSD to research into molecular geometries, molecular interactions, and molecular assemblies and demonstrate its value in the design of biologically active molecules and the solid forms in which they are delivered. Its potential in other commercially relevant areas is described, including gas storage and delivery, thin films, and (opto)electronics. The CSD also aids the solution of new crystal structures. Because no scientific instrument is without shortcomings, the limitations of CSD research are assessed. We emphasize the importance of maintaining database quality: notwithstanding the arrival of big data and machine learning, it remains perilous to ignore the principle of garbage in, garbage out. Finally, we explain why the CSD must evolve with the world around it to ensure it remains fit for purpose in the years ahead.

Measurement of Electric Fields Experienced by Urea Guest Molecules in the 18-Crown-6/Urea (1:5) Host-Guest Complex: An Experimental Reference Point for Electric-Field-Assisted Catalysis

High-resolution synchrotron and neutron single-crystal diffraction data of 18-crown-6/(pentakis)urea measured at 30 K are combined, with the aim of better appreciating the electrostatics associated with intermolecular interactions in condensed matter. With two 18-crown-6 molecules and five different urea molecules in the crystal, this represents the most ambitious combined X-ray/synchrotron and neutron experimental charge density analysis to date on a cocrystal or host-guest system incorporating such a large number of unique molecules. The dipole moments of the five urea guest molecules in the crystal are enhanced considerably compared to values determined for isolated molecules, and 2D maps of the electrostatic potential and electric field show clearly how the urea molecules are oriented with dipole moments aligned along the electric field exerted by their molecular neighbors. Experimental electric fields in the range of 10-19 GV m^-1, obtained for the five different urea environments, corroborate independent measurements of electric fields in the active sites of enzymes and provide an important experimental reference point for recent discussions focused on electric-field-assisted catalysis.

Electron-density critical points analysis and catastrophe theory to forecast structure instability in periodic solids

The critical points analysis of electron density, i.e. ρ(x), from ab initio calculations is used in combination with the catastrophe theory to show a correlation between ρ(x) topology and the appearance of instability that may lead to transformations of crystal structures, as a function of pressure/temperature. In particular, this study focuses on the evolution of coalescing non-degenerate critical points, i.e. such that ∇ρ(x_c) = 0 and λ₁, λ₂, λ₃ ≠ 0 [λ being the eigenvalues of the Hessian of ρ(x) at x_c], towards degenerate critical points, i.e. ∇ρ(x_c) = 0 and at least one λ equal to zero. The catastrophe theory formalism provides a mathematical tool to model ρ(x) in the neighbourhood of x_c and allows one to rationalize the occurrence of instability in terms of electron-density topology and Gibbs energy. The phase/state transitions that TiO₂ (rutile structure), MgO (periclase structure) and Al₂O₃ (corundum structure) undergo because of pressure and/or temperature are here discussed. An agreement of 3-5% is observed between the theoretical model and experimental pressure/temperature of transformation.

The future of topological analysis in experimental charge-density research

In a recent paper, Dittrich (2017) critically discussed the benefits of analysing experimental electron density within the framework of the quantum theory of atoms in molecules, often called simply the topological analysis of the charge density. The point he raised is important because it challenges the scientific production of a very active community. The question whether this kind of investigation is still sensible is intriguing and it fosters a multifaceted answer. Granted that none can predict the future of any field of science, but an alternative point of view emerges after answering three questions: Why should we investigate the electron charge (and spin) density? Is the interpretative scheme proposed by the quantum theory of atoms in molecules useful? Is an experimental charge density necessary?